Sory.<br>My ligand is Citalopram(a drug)<br><br><div class="gmail_quote">On Mon, Jan 17, 2011 at 1:35 PM, mohsen ramezanpour <span dir="ltr"><<a href="mailto:ramezanpour.mohsen@gmail.com">ramezanpour.mohsen@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear justin<br>I generated it by PRODRG server,and modifying Protein.top as it is in ENZYME/DRUG tutorial.<br>
Actually I don't know where is problem.<br>Ok,I will modify my EM parameter to do a good energy minimization.<br>
Thanks in advance<div><div></div><div class="h5"><br><br><div class="gmail_quote">On Sat, Jan 15, 2011 at 6:10 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div><br>
<br>
mohsen ramezanpour wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
I checked this one but it did not solve the problem.<br>
Actually I did in your way and I found the force is very high on one atom of my ligand,C12.<br>
I checked it's structure with pymol,it 's situation was normal.<br>
can I change it's coordinate a few?<br>
I think it can make force less.<br>
please let me know how can i do<br>
</blockquote>
<br></div>
Making ad hoc changes to coordinates is a bad idea. What you might gain by relaxing nonbonded forces you might strain within the molecule (bonded interactions). Proper energy minimization should resolve any high forces.<br>
<br>
What is your ligand? How did you generate its topology? Poor parameters can also give bad contacts and forces.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
thanks in advance<div><br>
<br>
On Wed, Jan 12, 2011 at 3:22 PM, mohsen ramezanpour <<a href="mailto:ramezanpour.mohsen@gmail.com" target="_blank">ramezanpour.mohsen@gmail.com</a> <mailto:<a href="mailto:ramezanpour.mohsen@gmail.com" target="_blank">ramezanpour.mohsen@gmail.com</a>>> wrote:<br>
<br>
<br>
<br>
<br>
On Wed, Jan 12, 2011 at 2:59 PM, Mark Abraham<br></div><div><div></div><div>
<<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a> <mailto:<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a>>> wrote:<br>
<br>
On 12/01/2011 9:37 PM, mohsen ramezanpour wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear Dr,Tsjerk<br>
I want to estimate protein-drug binding free energy.<br>
I am using umbrella sampling for this mean.<br>
my drug is inside of a hole in protein. then I have to rotate<br>
my system to can fit the pulling line along one box axis.<br>
besides I have to pull drug not at direction which connect COM<br>
of protein and drug,but it is better to pull it along line<br>
which connects drug to<br>
a residue inside of hole.<br>
</blockquote>
<br>
I would<br>
<br>
1) take my starting bound configuration,<br>
2) strip away anything except the protein complex,<br>
3) delete the box information,<br>
4) rotate the complex with editconf until I was happy with its<br>
orientation,<br>
5) then generate a suitable box around that orientation,<br>
6) do EM<br>
7) solvate and neutralize<br>
8) do EM<br>
9) etc.<br>
<br>
Mark<br>
Dear Mark<br>
<br>
thanks for your reply<br>
I will check this ways too.<br>
<br>
<br>
</div></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div><div></div><div>
<br>
I rotated box with editconf ,solvated system with<br>
genbox,neutralized with genion,<br>
now I want to generate NPT and then generating configuration<br>
as umbrella sampling tutorial.<br>
<br>
<br>
<br>
On Wed, Jan 12, 2011 at 1:53 PM, Tsjerk Wassenaar<br></div></div><div>
<<a href="mailto:tsjerkw@gmail.com" target="_blank">tsjerkw@gmail.com</a> <mailto:<a href="mailto:tsjerkw@gmail.com" target="_blank">tsjerkw@gmail.com</a>>> wrote:<br>
<br>
Hi Mohsen,<br>
<br>
You're doing something terribly wrong. But why you want to<br>
do what you<br>
attempt eludes me. Maybe it helps if you give an<br>
explanation of what<br>
you want, in stead of what doesn't work. In addition, give<br>
the set of<br>
commands that bring you up to this point, and not only the<br>
output of<br>
mdrun. That way we can probably see where you go astray.<br>
<br>
Cheers,<br>
<br>
Tsjerk<br>
<br>
<br>
On Wed, Jan 12, 2011 at 10:55 AM, mohsen ramezanpour<br>
<<a href="mailto:ramezanpour.mohsen@gmail.com" target="_blank">ramezanpour.mohsen@gmail.com</a><br></div><div><div></div><div>
<mailto:<a href="mailto:ramezanpour.mohsen@gmail.com" target="_blank">ramezanpour.mohsen@gmail.com</a>>> wrote:<br>
><br>
> Dear Dr.Tsjerk<br>
><br>
> Before doing md for generating NPT, I did an EM,the<br>
result was:<br>
> poteintial energy:-2.2611160*10^(6)<br>
> Max F=4.8960352*10^(4) on atom 5289<br>
> Besides I had done EM before on the same system,I just<br>
add solvent by genbox<br>
> and Ions by genion.<br>
> the above result is for Energy minimization after adding<br>
ions by genion.<br>
><br>
> I did what you said.but when I was generating NPT<br>
equilibration<br>
> I recieved this massage:<br>
><br>
><br>
> Step 0, time 0 (ps) LINCS WARNING<br>
> relative constraint deviation after LINCS:<br>
> rms 8.387059, max 321.381958 (between atoms 5289 and 5290)<br>
> bonds that rotated more than 30 degrees:<br>
> atom 1 atom 2 angle previous, current, constraint length<br>
> 5293 5294 74.4 0.1137 12.9330 0.1090<br>
> 5291 5293 59.1 0.1422 32.1605 0.1390<br>
> 5291 5292 103.5 0.1096 9.8747 0.1090<br>
> 5289 5291 89.5 0.1383 39.2022 0.1390<br>
> 5289 5290 85.7 0.1410 43.8439 0.1360<br>
> 5289 5287 85.4 0.1426 44.0196 0.1390<br>
> 5287 5288 90.2 0.1091 1.5178 0.1090<br>
> 5285 5287 88.7 0.1391 1.9186 0.1390<br>
> 5285 5286 41.8 0.1092 0.1322 0.1090<br>
> 5284 5293 87.2 0.1429 21.7425 0.1390<br>
> 5284 5285 92.4 0.1393 2.0880 0.1390<br>
> 5277 5278 79.1 0.1532 0.5213 0.1530<br>
> 5276 5284 63.0 0.1394 1.8946 0.1390<br>
> 5276 5277 79.2 0.1534 0.7278 0.1530<br>
> 5276 5275 79.9 0.1432 0.6036 0.1430<br>
> 5274 5275 82.0 0.1431 0.3189 0.1430<br>
> 5276 5272 82.1 0.1393 0.6113 0.1390<br>
> 5272 5273 74.2 0.1333 0.2251 0.1330<br>
> 5272 5270 77.7 0.1332 0.2428 0.1330<br>
> 5270 5271 46.5 0.1091 0.1682 0.1090<br>
> 5268 5270 36.9 0.1391 0.1854 0.1390<br>
> 5273 5266 41.4 0.1334 0.1933 0.1330<br>
> Wrote pdb files with previous and current coordinates<br>
> step 0Warning: 1-4 interaction between 5272 and 5293 at<br>
distance 22.443<br>
> which is larger than the 1-4 table size 2.400 nm<br>
> These are ignored for the rest of the simulation<br>
> This usually means your system is exploding,<br>
> if not, you should increase table-extension in your mdp file<br>
> or with user tables increase the table size<br>
><br>
> Please let me know what can I do.<br>
> thanks in advance<br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
> On Thu, Jan 6, 2011 at 11:20 PM, Tsjerk Wassenaar<br></div></div><div>
<<a href="mailto:tsjerkw@gmail.com" target="_blank">tsjerkw@gmail.com</a> <mailto:<a href="mailto:tsjerkw@gmail.com" target="_blank">tsjerkw@gmail.com</a>>> wrote:<br>
>><br>
>> Hi Mohsen,<br>
>><br>
>> I think rotating a molecule with editconf will not<br>
rotate the box. Then<br>
>> again, if it did, it would result in a box violating<br>
Gromacs requirements.<br>
>> Either way, it's not going to work like that. Build a<br>
new box after<br>
>> rotation... And have a good look at what you're<br>
actually trying now by<br>
>> taking the rotated system and stack it a few times<br>
using genconf -nbox 2 2 2<br>
>><br>
>> Cheers,<br>
>><br>
>> Tsjerk<br>
>><br>
>> On Jan 6, 2011 7:22 PM, "mohsen ramezanpour"<br>
>> <<a href="mailto:ramezanpour.mohsen@gmail.com" target="_blank">ramezanpour.mohsen@gmail.com</a><br></div><div>
<mailto:<a href="mailto:ramezanpour.mohsen@gmail.com" target="_blank">ramezanpour.mohsen@gmail.com</a>>> wrote:<br>
>><br>
>> Dear Amit<br>
>> I entered these commands for rotating box:<br>
>> editconf -f conf.gro -o output.pdb -rotate 0 0 25.4<br>
>> and then:<br>
>> editconf -f output.pdb -o newbox.pdb -rotate 0 127.67548 0<br>
>> as a result my molecul is located out of box totally,of<br>
course drug and<br>
>> protein are bind to eachother yet.<br>
>> thanks in advance for your attention and reply<br>
>><br>
>> On Tue, Jan 4, 2011 at 1:25 PM, Amit Choubey<br></div><div>
<<a href="mailto:kgp.amit@gmail.com" target="_blank">kgp.amit@gmail.com</a> <mailto:<a href="mailto:kgp.amit@gmail.com" target="_blank">kgp.amit@gmail.com</a>>> wrote: ><br>
>> > Could you post the e...<br>
>><br>
>> --<br>
>> gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><div><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<br>
>> <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a><br>
before posting!<br>
>> Please don't post (un)subscribe requests to the list.<br>
Use the<br>
>> www interface or send it to<br>
<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<div><br>
>> Can't post? Read<br>
<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 <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><div><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<br>
>> <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a><br>
before posting!<br>
>> Please don't post (un)subscribe requests to the list.<br>
Use the<br>
>> www interface or send it to<br>
<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<div><br>
>> Can't post? Read<br>
<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
><br>
><br>
> --<br>
> gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><div><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<br>
> <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a><br>
before posting!<br>
> Please don't post (un)subscribe requests to the list.<br>
Use the<br>
> www interface or send it to<br>
<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<div><br>
> Can't post? Read<br>
<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
><br>
<br>
<br>
<br>
--<br>
Tsjerk A. Wassenaar, Ph.D.<br>
<br>
post-doctoral researcher<br>
Molecular Dynamics Group<br>
* Groningen Institute for Biomolecular Research and<br>
Biotechnology<br>
* Zernike Institute for Advanced Materials<br>
University of Groningen<br>
The Netherlands<br>
--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><div><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<br>
<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<div><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>
<br>
</div></blockquote><div>
<br>
<br>
--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><div><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<br>
<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 <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<div><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>
<br>
<br>
</div></blockquote>
<br>
-- <br>
========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br><font color="#888888">
-- <br></font><div><div></div><div>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><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 www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<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>
</div></div></blockquote></div><br>
</div></div></blockquote></div><br>