----- Original Message ----- From: Binh Khanh Mai <khanhbinh...@yahoo.com> Date: Tuesday, May 18, 2010 14:10 Subject: [gmx-users] Computational Alanine Scanning To: gmx-users@gromacs.org
----------------------------------------------------------- | > > > Dear all, > I'm trying use GROMACS to calculate binding free energy with MM-PBSA > method and do Computational Alanine Scanning. After run MD and obtain pdb > files, I use PYMOL to mutate one residue to Alanine to estimate the energy > difference between 2 structures. My issue is that I don't know how to > generate topology file (.top file) after mutation but not ignore hydrogen in > initial file. If using pdb2gmx command, it needs me ignore hydrogen and I > think hydrogen after using this command is not in the same place with > hydrogen at first. This will my my results unbelievable. Do your mutation with PYMOL, and later plan to use those coordinates as inputs to grompp. Use pdb2gmx however you need to in order to create the new .top. Assuming PYMOL is not messing around with atom/residue naming and ordering except for the mutated residue, all you have to do is edit a copy of the PYMOL coordinate file so that the new residue has atom names, residue names and atom ordering that match the .top. Then supply the .top with the fixed coordinate file to grompp. Alternatively, work out why pdb2gmx requires you to ignore and regenerate hydrogens, and fix the issue there. Mark -- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/mailing_lists/users.php
