Mark, thanks for the prompt response!<br><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im"><blockquote style="border-left:1px solid #00f;padding-left:13px;margin-left:0" type="cite">
<div>I have done Normal Mode Analysis and have calculated partial charges and the optimized geometry of a few compounds using high-level QM calculations. Now I want to see (if possible) how well GROMACS can reproduce the normal modes if I start from the same optimized geometry and use the same partial charges.</div>
</blockquote><div><br></div></div><div>In general for NMA to make sense you need to be at a stationary point w.r.t. the atomic degrees of freedom of the model being used. That won't be quite true at a QM geometry, so there's a sense of apples-vs-oranges comparison.</div>
<div class="im"><div><br></div></div></blockquote><div> </div><div>If I get it right you mean that NMA in GROMACS must start from an energy minimum (stationary point) w.r.t the ff used (GAFF in my case), which means that an energy minimization is neccessary ever if I use an QM optimum geometry and the respective partial charges. Namely there is no way to reproduce the normal modes I obtained from QM calculations, correct?<br>
<br></div><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="im"><div><br>An obvious problem is that the starting compound geometry is not in full precision</div>
<div><br></div></div><div>The starting geometry is in full precision if it's the same as that used for the QM calculation. That is quite possible to achieve with .pdb or .gro input.</div><div class="im"><div><br></div>
</div></blockquote><div> </div><div>The same as the starting geometry or as the optimized geometry?<br> </div><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div class="im"><div></div><blockquote style="border-left:1px solid #00f;padding-left:13px;margin-left:0" type="cite"><div> as highlighted in the documentation:<br>
<br><a href="http://www.gromacs.org/Documentation/How-tos/Normal_Mode_Analysis" target="_blank">http://www.gromacs.org/Documentation/How-tos/Normal_Mode_Analysis</a><br><br>Is it possible to create a full precision .trr coordinate file from a .gro or any other structure file with modified 8-decimal point coordinates?</div>
</blockquote><div><br></div></div><div>I think you are misunderstanding the use of the word "precision" here. In general, the same configuration will be represented differently in .trr and .gro formats, with the former being a closer approximation. Accordingly, one will get a different result for NMA on the endpoint of GROMACS EM as observed in the .trr file and as observed in the .gro file. The former will be closer to the stationary point, and so lead to more acceptable estimates of the normal modes. However, here you want to do NMA on the same coordinates with two programs, so it is up to you to represent the coordinates in a way that the two programs can compute on the same approximation to the coordinates of the stationary point. There's no need to convert to .trr (or the QM binary format), because all that does is treat 2.613 as 2.6130000.</div>
<div><br></div></blockquote><div><br>So a command line like this will do the job, right?<br><br>grompp_d4.5.5 -f nm.mdp -c ${ligand}_8_decimal_points.gro -p ${ligand}.top -o nm.tpr<br><br><br><br>Thomas<br><br><br><br> </div>
</div><p style="margin-bottom:0cm" align="LEFT">======================================================================</p>
<p style="margin-bottom:0cm" align="LEFT">Thomas Evangelidis</p>
<p style="margin-bottom:0cm" align="LEFT">PhD student</p>
<p style="margin-bottom:0cm" align="LEFT">Biomedical Research
Foundation, Academy of Athens</p>
<p style="margin-bottom:0cm" align="LEFT">4 Soranou Ephessiou , 115 27
Athens, Greece<br><br>email: <a href="mailto:tevang@bioacademy.gr" target="_blank">tevang@bioacademy.gr</a></p>
<p style="margin-bottom:0cm" align="LEFT"> <a href="mailto:tevang3@gmail.com" target="_blank">tevang3@gmail.com</a></p>
<p style="margin-bottom:0cm" align="LEFT"><br>website:
<a href="https://sites.google.com/site/thomasevangelidishomepage/" target="_blank">https://sites.google.com/site/thomasevangelidishomepage/</a></p>
<p style="margin-bottom:0cm" align="LEFT"><br>
</p>
<br>