<html><head></head><body><div style="font-family: Verdana;font-size: 12.0px;"><div>
<div>Has anyone done an all atom 53a6 or 54a7 .itp/rtp that can be downloaded (reguarding protein/amono acids or and nucleic? I can make one in a week (for proteins) just using a topology builder and modifying things, there are already lipids around 100 with all atom in either field, but after reviewing liturature, even the hybrid 53a6 is better than opls or charm all atom as far as fitting computational against real measured biochem. I assume if one did this, it would have to be tested first, or in parallel to any work, however I have found multiple relics and hints of people whom obviously already did this, but can not find an .itp or .rtp for this? this is also in reguards to membrane simulations, as most of the comparisons I have seen show even berger hybrid UA/charged hydrogens in 53a6 as much (by about 1-5% difference) vs around 10-20% difference in biochem Vs comp. for all other force fields.</div>
<div> </div>
<div>Sincerely,</div>
<div> </div>
<div>Stephan Watkins</div>
<div>
<div name="quote" style="margin:10px 5px 5px 10px; padding: 10px 0 10px 10px; border-left:2px solid #C3D9E5; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div style="margin:0 0 10px 0;"><b>Gesendet:</b> Donnerstag, 05. Juni 2014 um 23:31 Uhr<br/>
<b>Von:</b> "Piggot T." <T.Piggot@soton.ac.uk><br/>
<b>An:</b> "gmx-users@gromacs.org" <gmx-users@gromacs.org><br/>
<b>Betreff:</b> Re: [gmx-users] Need some help with membrane simulation</div>
<div name="quoted-content">I haven't tried it at all, but in the release notes for 5.0.x there is a LJ force switching option mentioned. This would probably be preferable for membrane simulations using CHARMM36, rather than switching the potential.<br/>
<br/>
Cheers<br/>
<br/>
Tom<br/>
________________________________________<br/>
From: gromacs.org_gmx-users-bounces@maillist.sys.kth.se [gromacs.org_gmx-users-bounces@maillist.sys.kth.se] on behalf of Justin Lemkul [jalemkul@vt.edu]<br/>
Sent: 05 June 2014 22:12<br/>
To: gmx-users@gromacs.org<br/>
Subject: Re: [gmx-users] Need some help with membrane simulation<br/>
<br/>
On 6/5/14, 4:57 PM, Ali Khan wrote:<br/>
> Justin,<br/>
><br/>
> 4.6.5 does not have Potential-switch-Verlet. When I added the other<br/>
> changes, I still had the same error.<br/>
><br/>
<br/>
Likely something more fundamental is wrong with the system. What was the<br/>
outcome of EM? Exact numbers produced by mdrun, please.<br/>
<br/>
> Are these recommendations ones that I should generally use when<br/>
> implementing charmm in gromacs?<br/>
<br/>
Yes. I'm still experimenting with all the new Gromacs features, but as best I<br/>
can tell, in terms of reproducing a "normal" CHARMM setup:<br/>
<br/>
cutoff-scheme = Verlet<br/>
vdwtype = cutoff<br/>
vdw-modifier = Potential-switch<br/>
coulombtype = PME<br/>
rlist = 1.2<br/>
rvdw = 1.2<br/>
rvdw-switch = 1.0<br/>
rcoulomb = 1.2<br/>
<br/>
-Justin<br/>
<br/>
><br/>
><br/>
> On Thu, Jun 5, 2014 at 4:29 PM, Justin Lemkul <jalemkul@vt.edu> wrote:<br/>
><br/>
>><br/>
>><br/>
>> On 6/5/14, 2:47 PM, Ali Khan wrote:<br/>
>><br/>
>>> I am using the charmm36 force field.<br/>
>>><br/>
>>> *My energy minimization .mdp file:*<br/>
>>><br/>
>>> ; minim.mdp - used as input into grompp to generate em.tpr<br/>
>>> ; Parameters describing what to do, when to stop and what to save<br/>
>>> integrator = steep ; Algorithm (steep = steepest descent minimization)<br/>
>>> emtol = 1000.0 ; Stop minimization when the maximum force < 1000.0<br/>
>>> kJ/mol/nm<br/>
>>> emstep = 0.01 ; Energy step size<br/>
>>> nsteps = 50000000 ; Maximum number of (minimization) steps to perform<br/>
>>><br/>
>>> ; Parameters describing how to find the neighbors of each atom and how to<br/>
>>> calculate the interactions<br/>
>>> nstlist = 1 ; Frequency to update the neighbor list and long range forces<br/>
>>> ns_type = grid ; Method to determine neighbor list (simple, grid)<br/>
>>> rlist = 1.2 ; Cut-off for making neighbor list (short range forces)<br/>
>>> coulombtype = PME ; Treatment of long range electrostatic interactions<br/>
>>> rcoulomb = 1.2 ; Short-range electrostatic cut-off<br/>
>>> rvdw = 1.2 ; Short-range Van der Waals cut-off<br/>
>>> pbc = xyz ; Periodic Boundary Conditions (yes/no)<br/>
>>> vdw-type = cutoff<br/>
>>> vdw-modifier = Potential-shift-Verlet<br/>
>>> cutoff-scheme = Verlet<br/>
>>><br/>
>>><br/>
>>><br/>
>>><br/>
>>><br/>
>>> *Temperature equilibration .mdp file:*<br/>
>>><br/>
>>> title = NVT equilibration<br/>
>>> define = -DPOSRES ; position restrain the protein<br/>
>>> ; Run parameters<br/>
>>> integrator = md ; leap-frog integrator<br/>
>>> nsteps = 500000 ; 2 * 500000 = 1000 ps<br/>
>>> dt = 0.002 ; 2 fs<br/>
>>> ; Output control<br/>
>>> nstxout = 0<br/>
>>> nstvout = 0<br/>
>>> nstenergy = 500<br/>
>>> nstlog = 0<br/>
>>> ; Bond parameters<br/>
>>> continuation = no ; first dynamics run<br/>
>>> constraint_algorithm = lincs ; holonomic constraints<br/>
>>> constraints = all-bonds ; all bonds (even heavy atom-H bonds) constrained<br/>
>>> lincs_iter = 1 ; accuracy of LINCS<br/>
>>> lincs_order = 4 ; also related to accuracy<br/>
>>> ; Neighborsearching<br/>
>>> ns_type = grid ; search neighboring grid cells<br/>
>>> nstlist = 5 ; 10 fs<br/>
>>> rlist = 1.2 ; short-range neighborlist cutoff (in nm)<br/>
>>> rcoulomb = 1.2 ; short-range electrostatic cutoff (in nm)<br/>
>>> rvdw = 1.2 ; short-range van der Waals cutoff (in nm)<br/>
>>> rcoulomb-switch = 0<br/>
>>> rvdw-switch = 0.8<br/>
>>><br/>
>><br/>
>> Use rvdw-switch = 1.0.<br/>
>><br/>
>><br/>
>> ; long-range cut-off for switched potentials<br/>
>>> rlistlong = 1.4<br/>
>>><br/>
>><br/>
>> Probably never had the desired effect in terms of using CHARMM force<br/>
>> fields and isn't necessary. The Verlet scheme with rlist = 1.2 should<br/>
>> provide adequate buffer space.<br/>
>><br/>
>><br/>
>> ; Electrostatics<br/>
>>> coulombtype = PME ; Particle Mesh Ewald for long-range electrostatics<br/>
>>> pme_order = 4 ; cubic interpolation<br/>
>>> fourierspacing = 0.16 ; grid spacing for FFT<br/>
>>> ; Relative dielectric constant for the medium and the reaction field<br/>
>>> epsilon_r = 1<br/>
>>> epsilon_rf = 1<br/>
>>> ; Temperature coupling is on<br/>
>>> tcoupl = V-rescale ; modified Berendsen thermostat<br/>
>>> tc-grps = POPC Water_and_ions ; two coupling groups - more accurate<br/>
>>> tau_t = 0.1 0.1 ; time constant, in ps<br/>
>>> ref_t = 310 310 ; reference temperature, one for each group, in K<br/>
>>> ; Pressure coupling is off<br/>
>>> pcoupl = no ; no pressure coupling in NVT<br/>
>>> ; Periodic boundary conditions<br/>
>>> pbc = xyz ; 3-D PBC<br/>
>>> ; Dispersion correction<br/>
>>> DispCorr = No ; account for cut-off vdW scheme<br/>
>>> ; Method for doing Van der Waals<br/>
>>> vdw-type = cutoff<br/>
>>> vdw-modifier = Potential-shift-Verlet<br/>
>>><br/>
>><br/>
>> Use Potential-switch-Verlet. I don't remember if this was added to the<br/>
>> code in time for 4.6.5, but it works well in any of the 5.0 releases.<br/>
>><br/>
>> -Justin<br/>
>><br/>
>><br/>
>> cutoff-scheme = Verlet<br/>
>>> ; Velocity generation<br/>
>>> gen_vel = yes ; assign velocities from Maxwell distribution<br/>
>>> gen_temp = 310 ; temperature for Maxwell distribution<br/>
>>> gen_seed = -1 ; generate a random seed<br/>
>>><br/>
>>><br/>
>>><br/>
>>> On Thu, Jun 5, 2014 at 2:36 PM, Justin Lemkul <jalemkul@vt.edu> wrote:<br/>
>>><br/>
>>><br/>
>>>><br/>
>>>> On 6/5/14, 2:11 PM, Ali Khan wrote:<br/>
>>>><br/>
>>>> Hi all,<br/>
>>>>><br/>
>>>>> I don't know if anyone has experienced this problem before, but I<br/>
>>>>> figure I<br/>
>>>>> ask. I have several times run a simulation on a 135A x 135A POPC<br/>
>>>>> membrane<br/>
>>>>> acquired fro the charmm-gui website. I wanted now to simulate one that<br/>
>>>>> was<br/>
>>>>> slightly larger (150A x 150A). I made the larger membrane on charmm-gui,<br/>
>>>>> and I did everything the same way as I did before (even used the same<br/>
>>>>> .mdp<br/>
>>>>> files). After energy minimization, I get an error during my temperature<br/>
>>>>> equilibration as follows:<br/>
>>>>><br/>
>>>>> *Back Off! I just backed up step0c_n1.pdb to ./#step0c_n1.pdb.4#*<br/>
>>>>> *Wrote pdb files with previous and current coordinates*<br/>
>>>>><br/>
>>>>> *WARNING: Listed nonbonded interaction between particles 31617 and<br/>
>>>>> 31619*<br/>
>>>>> *at distance 3f which is larger than the table limit 3f nm.*<br/>
>>>>><br/>
>>>>> *This is likely either a 1,4 interaction, or a listed interaction<br/>
>>>>> inside*<br/>
>>>>> *a smaller molecule you are decoupling during a free energy<br/>
>>>>> calculation.*<br/>
>>>>> *Since interactions at distances beyond the table cannot be computed,*<br/>
>>>>> *they are skipped until they are inside the table limit again. You will*<br/>
>>>>> *only see this message once, even if it occurs for several<br/>
>>>>> interactions.*<br/>
>>>>><br/>
>>>>> *IMPORTANT: This should not happen in a stable simulation, so there is*<br/>
>>>>> *probably something wrong with your system. Only change the<br/>
>>>>> table-extension*<br/>
>>>>><br/>
>>>>> *distance in the mdp file if you are really sure that is the reason.*<br/>
>>>>><br/>
>>>>><br/>
>>>>><br/>
>>>>> *-------------------------------------------------------*<br/>
>>>>> *Program mdrun, VERSION 4.6.5*<br/>
>>>>> *Source code file: /home/ali/Downloads/gromacs-4.6.5/src/mdlib/pme.c,<br/>
>>>>> line:<br/>
>>>>> 851*<br/>
>>>>><br/>
>>>>> *Fatal error:*<br/>
>>>>> *5 particles communicated to PME node 1 are more than 2/3 times the<br/>
>>>>> cut-off<br/>
>>>>> out of the domain decomposition cell of their charge group in dimension<br/>
>>>>> x.*<br/>
>>>>><br/>
>>>>> *This usually means that your system is not well equilibrated.*<br/>
>>>>> *For more information and tips for troubleshooting, please check the<br/>
>>>>> GROMACS*<br/>
>>>>> *website at <a href="http://www.gromacs.org/Documentation/Errors" target="_blank">http://www.gromacs.org/Documentation/Errors</a><br/>
>>>>> <<a href="http://www.gromacs.org/Documentation/Errors" target="_blank">http://www.gromacs.org/Documentation/Errors</a>>*<br/>
>>>>> *-------------------------------------------------------*<br/>
>>>>><br/>
>>>>><br/>
>>>>> I referred to the gromacs website, and it stated that this error means<br/>
>>>>> that<br/>
>>>>> my system is blowing up. I have tried several things to fix it<br/>
>>>>> (including:<br/>
>>>>> remaking the system on charmm-gui, running energy minimization with<br/>
>>>>> small<br/>
>>>>> energy steps [emstep = 0.01, 0.001, 0.0001], energy minimizing to Fmax<br/>
>>>>> of<br/>
>>>>> 100 kJ kJ/mol/nm instead of 1000, and I have done all these combinations<br/>
>>>>> in<br/>
>>>>> both single and double precision) and none of these methods worked. I<br/>
>>>>> also<br/>
>>>>> should add that sometimes I am able to perform the temperature<br/>
>>>>> equilibration, but then I would get the same error during the pressure<br/>
>>>>> equilibration. I doubt it is charmm-gui, because they are usually very<br/>
>>>>> reliable in generating membranes.<br/>
>>>>><br/>
>>>>><br/>
>>>>> On the contrary, we've had numerous reports of crashing membranes, due<br/>
>>>> to<br/>
>>>> some clashes in the starting structures. Manual modifications to the<br/>
>>>> coordinates of the offending atoms (fractions of an Angstrom, really)<br/>
>>>> usually solves it. The output of EM tells you where the maximum force<br/>
>>>> is.<br/>
>>>> How big is this force? What atom is it acting upon? That's where you<br/>
>>>> start your investigation.<br/>
>>>><br/>
>>>><br/>
>>>> Does anyone have any suggestions? I can post any of my input files if<br/>
>>>><br/>
>>>>> requested.<br/>
>>>>><br/>
>>>>><br/>
>>>>> Always post an .mdp file when a run is crashing. Membranes can be<br/>
>>>> particularly challenging to get right; small errors in settings can have<br/>
>>>> catastrophic consequences.<br/>
>>>><br/>
>>>> -Justin<br/>
>>>><br/>
>>>> --<br/>
>>>> ==================================================<br/>
>>>><br/>
>>>> Justin A. Lemkul, Ph.D.<br/>
>>>> Ruth L. Kirschstein NRSA Postdoctoral Fellow<br/>
>>>><br/>
>>>> Department of Pharmaceutical Sciences<br/>
>>>> School of Pharmacy<br/>
>>>> Health Sciences Facility II, Room 601<br/>
>>>> University of Maryland, Baltimore<br/>
>>>> 20 Penn St.<br/>
>>>> Baltimore, MD 21201<br/>
>>>><br/>
>>>> jalemkul@outerbanks.umaryland.edu | (410) 706-7441<br/>
>>>> <a href="http://mackerell.umaryland.edu/~jalemkul" target="_blank">http://mackerell.umaryland.edu/~jalemkul</a><br/>
>>>><br/>
>>>> ==================================================<br/>
>>>><br/>
>>>> --<br/>
>>>> Gromacs Users mailing list<br/>
>>>><br/>
>>>> * Please search the archive at <a href="http://www.gromacs.org/" target="_blank">http://www.gromacs.org/</a><br/>
>>>> Support/Mailing_Lists/GMX-Users_List before posting!<br/>
>>>><br/>
>>>> * Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br/>
>>>><br/>
>>>> * For (un)subscribe requests visit<br/>
>>>> <a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users</a> or<br/>
>>>> send a mail to gmx-users-request@gromacs.org.<br/>
>>>><br/>
>>>><br/>
>> --<br/>
>> ==================================================<br/>
>><br/>
>> Justin A. Lemkul, Ph.D.<br/>
>> Ruth L. Kirschstein NRSA Postdoctoral Fellow<br/>
>><br/>
>> Department of Pharmaceutical Sciences<br/>
>> School of Pharmacy<br/>
>> Health Sciences Facility II, Room 601<br/>
>> University of Maryland, Baltimore<br/>
>> 20 Penn St.<br/>
>> Baltimore, MD 21201<br/>
>><br/>
>> jalemkul@outerbanks.umaryland.edu | (410) 706-7441<br/>
>> <a href="http://mackerell.umaryland.edu/~jalemkul" target="_blank">http://mackerell.umaryland.edu/~jalemkul</a><br/>
>><br/>
>> ==================================================<br/>
>> --<br/>
>> Gromacs Users mailing list<br/>
>><br/>
>> * Please search the archive at <a href="http://www.gromacs.org/" target="_blank">http://www.gromacs.org/</a><br/>
>> Support/Mailing_Lists/GMX-Users_List before posting!<br/>
>><br/>
>> * Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br/>
>><br/>
>> * For (un)subscribe requests visit<br/>
>> <a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users</a> or<br/>
>> send a mail to gmx-users-request@gromacs.org.<br/>
>><br/>
<br/>
--<br/>
==================================================<br/>
<br/>
Justin A. Lemkul, Ph.D.<br/>
Ruth L. Kirschstein NRSA Postdoctoral Fellow<br/>
<br/>
Department of Pharmaceutical Sciences<br/>
School of Pharmacy<br/>
Health Sciences Facility II, Room 601<br/>
University of Maryland, Baltimore<br/>
20 Penn St.<br/>
Baltimore, MD 21201<br/>
<br/>
jalemkul@outerbanks.umaryland.edu | (410) 706-7441<br/>
<a href="http://mackerell.umaryland.edu/~jalemkul" target="_blank">http://mackerell.umaryland.edu/~jalemkul</a><br/>
<br/>
==================================================<br/>
--<br/>
Gromacs Users mailing list<br/>
<br/>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List</a> before posting!<br/>
<br/>
* Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br/>
<br/>
* For (un)subscribe requests visit<br/>
<a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users</a> or send a mail to gmx-users-request@gromacs.org.<br/>
--<br/>
Gromacs Users mailing list<br/>
<br/>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List</a> before posting!<br/>
<br/>
* Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br/>
<br/>
* For (un)subscribe requests visit<br/>
<a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users</a> or send a mail to gmx-users-request@gromacs.org.</div>
</div>
</div>
</div></div></body></html>