Hello Dr. Gase,
I’ve tried a few different approaches to calculate inter-domain non-bonded interactions.
1) In cpptraj:
lie LIE :1-425 :426-850 out LIEanalysis_NoProt_Del4.dat cutvdw 12.0 cutelec 12.0 diel 4
Where the two masks correspond to all the residues in the two protein domains that are in contact.
Or,
2) Again in cpptraj,
esander XXXX out XXXX.dat parmname XXXX.parm7 cut 9999.0 igb 2 extdiel 78.5 intdiel 4.0 ntb 0 ntf 2 ntc 2
Where XXXX is either DM1, DM2, or DM1DM2 for domain #1, #2, or both #1 and #2. I saved separate trajectories in a prior step with one or the other or both of the domains from the original trajectory.
(I haven’t used the “pairwise” function in cpptraj, just because my understanding is that the electrostatics are not screened by the environment in any way.)
The output in mdout reads (for one time step):
NSTEP = 18000000 TIME(PS) = 54000.000 TEMP(K) = 299.00 PRESS = 0.0
Etot = -444472.8967 EKtot = 109676.4062 EPtot = -554149.3030
BOND = 2571.3718 ANGLE = 6977.2214 DIHED = 8519.2870
1-4 NB = 2865.7846 1-4 EEL = 38657.1336 VDWAALS = 74142.2504
EELEC = -687882.3518 EHBOND = 0.0000 RESTRAINT = 0.0000
------------------------------------------------------------------------------
So, even in the simulation, the VDW term is positive. This puzzles me.
Did this help to clarify my question? Please let me know if I should provide more information.
Best,
Matthew
> On Jun 30, 2021, at 9:01 AM, David A Case <david.case.rutgers.edu> wrote:
>
> *Message sent from a system outside of UConn.*
>
>
> On Mon, Jun 28, 2021, Matthew Guberman-Pfeffer wrote:
>
>>
>> I am trying to compute how the stability of a protein-protein interface
>> changes with protonation of the acidic groups within a 10 Å distance of
>> the interface. I’ve computed the non-bonded interactions for MD snapshots
>> using either the linear interaction method with some dielectric constant
>> specified (4 or 80) or using esander. One aspect that confuses me is that
>> the VDW energies are positive. I was expecting strongly negative values. Is
>> there just a sign convention I’m missing?
>
> It's not clear what you are computing, or exactly how. How do the numbers
> you get compare to the values that are printed in the mdout file? Are you
> only computing a subset of the LJ interactions?
>
> ....dac
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&data=04%7C01%7C%7Cd2b3944ace404d8a2da008d93bc74574%7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C637606549307248214%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=NhAukwOWBi0uUY6L8IWkN9ZHCNpZnSFSuRdYIu26TA4%3D&reserved=0
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jun 30 2021 - 07:00:02 PDT