<html><head></head><body><div style="font-family: Verdana;font-size: 12.0px;"><div>
<div>You can also look them up (angles bond distances) in the CRC handbook or online and just put the angle(distances direct into the .itp file for your ligand, as well as impropers...it works for all the force fields, unless you use hybrids (CH3=1 representation)...it comes out wierd as the different force fields put the same numbers different. Such as distance 1.18, 118.0 to 11800000.0 depending on representation is all...</div>
<div> </div>
<div>Stephan</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> Mittwoch, 17. Juli 2013 um 19:33 Uhr<br/>
<b>Von:</b> "gromacs query" <gromacsquery@gmail.com><br/>
<b>An:</b> "Discussion list for GROMACS users" <gmx-users@gromacs.org><br/>
<b>Betreff:</b> Re: [gmx-users] defining impropers necessary?</div>
<div name="quoted-content">Dear Justin,<br/>
<br/>
Thanks for reply and explanation, and..:<br/>
<br/>
>> you've got an amide flanked by two methylene groups as the repeat unit<br/>
>> All the amino acids in the aminoacids.rtp file specify impropers<br/>
centered on the C and N atoms of the peptide bond.<br/>
<br/>
I meant to say its easy to define (in general) impropers for chiral<br/>
centers, but not sure about bonds (here peptide). In oplsaa.ff I just<br/>
looked ALA in aminoacids.rtp for which there are two impropers.<br/>
<br/>
ALA impropers:<br/>
<br/>
-C CA N H improper_Z_N_X_Y<br/>
CA +N C O improper_O_C_X_Y<br/>
<br/>
which should be as follows:<br/>
<br/>
H CH3 H<br/>
I I I<br/>
(-)C---------N--CA--C-------N(+)<br/>
I I I<br/>
O H O<br/>
^ ^ ^<br/>
^ ^ ^<br/>
[-C]---------[ALA]-------[N+]<br/>
<br/>
<br/>
Does this mean both improper terms are making peptide bond planar (two<br/>
peptide bonds: left side and right side of ALA), but in my case as you<br/>
noticed there are two methylene groups. So is it correct if I define<br/>
improper for my polymer as:<br/>
<br/>
C C2 N H<br/>
C1 N C O<br/>
<br/>
<br/>
thanks<br/>
<br/>
regards,<br/>
Jiom<br/>
<br/>
<br/>
On Wed, Jul 17, 2013 at 6:51 PM, Justin Lemkul <jalemkul@vt.edu> wrote:<br/>
<br/>
><br/>
><br/>
> On 7/17/13 11:39 AM, gromacs query wrote:<br/>
><br/>
>> Dear Justin,<br/>
>><br/>
>> 1) I can understand the improper for stereocenters (chiral) easily but<br/>
>> with<br/>
>><br/>
><br/>
> OPLS doesn't use impropers for chiral centers.<br/>
><br/>
><br/>
> bonds I am confused (I have tried to explain below please let me know if I<br/>
>> have defined impropers correctly for peptide to keep it planar).<br/>
>><br/>
>> H1 O H3<br/>
>> \ I /<br/>
>> --C1---C--N---C2--<br/>
>> / I \<br/>
>> H2 H H4<br/>
>><br/>
>><br/>
>><br/>
> OK, funky text drawing :) But it looks like you've got an amide flanked<br/>
> by two methylene groups as the repeat unit.<br/>
><br/>
><br/>
> So I think improper (being zero) to keep peptide planar should be defined<br/>
>> (and sufficiently) as: O C N C1<br/>
>> and dihedral should be which is simple as: O C N H<br/>
>><br/>
>> But it confuses me, why one needs to define improper for peptide. If I<br/>
>> just<br/>
>> keep dihedral angle strict (high force constant) then it will make peptide<br/>
>> freezed! so it will remain planar.<br/>
>><br/>
>><br/>
> The OPLS approach for amides is to define normal dihedrals as well as<br/>
> impropers. Keep in mind that rotations about a bond and out-of-plane<br/>
> bending are different types of motion. I suppose you could achieve what<br/>
> you're thinking of by very strong force constants for the dihedrals, but I<br/>
> would wonder if you're then throwing the balance of 1-4 interactions<br/>
> (especially charges) out of whack.<br/>
><br/>
><br/>
> 2) Also, I am not sure whether OPLS needs this peptide bond to be defined<br/>
>> with improper, someone have experienced this please suggest?<br/>
>><br/>
>><br/>
> All the amino acids in the aminoacids.rtp file specify impropers centered<br/>
> on the C and N atoms of the peptide bond. I would stick with that approach.<br/>
><br/>
><br/>
> -Justin<br/>
><br/>
> --<br/>
> ==============================**====================<br/>
><br/>
> Justin A. Lemkul, Ph.D.<br/>
> 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 <jalemkul@outerbanks.umaryland.edu> | (410)<br/>
> 706-7441<br/>
><br/>
> ==============================**====================<br/>
> --<br/>
> gmx-users mailing list gmx-users@gromacs.org<br/>
> <a href="http://lists.gromacs.org/**mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/**mailman/listinfo/gmx-users</a><<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a>><br/>
> * Please search the archive at <a href="http://www.gromacs.org/**" target="_blank">http://www.gromacs.org/**</a><br/>
> Support/Mailing_Lists/Search<<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a>>before posting!<br/>
> * Please don't post (un)subscribe requests to the list. Use the www<br/>
> interface or send it to gmx-users-request@gromacs.org.<br/>
> * Can't post? Read <a href="http://www.gromacs.org/**Support/Mailing_Lists" target="_blank">http://www.gromacs.org/**Support/Mailing_Lists</a><<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a>><br/>
><br/>
--<br/>
gmx-users mailing list gmx-users@gromacs.org<br/>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br/>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br/>
* Please don't post (un)subscribe requests to the list. Use the<br/>
www interface or send it to gmx-users-request@gromacs.org.<br/>
* Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a></div>
</div>
</div>
</div></div></body></html>