Re: [AMBER] Hydration Free Enery by TI

From: Kepa K. Burusco <kekoburgo.yahoo.es>
Date: Tue, 19 Mar 2013 19:27:18 +0000 (GMT)

[19-III-2013] Hi Thomas and Ignacio, I have tried to create the input files for the boxes V0 = [solvent+solute], and V1 = [solvent] according to the method proposed by Tom: creating system crd + top files, then using "delete system system.1" and saving again: 1) I have been completely unable to do this process with AmberTools12 "tleap" using the "delete" command (actually I cannot find it in the manual) BUT using "remove" I can delete the solute from the solvent box using the sequence: ********************************************************** source leaprc.ff99SB loadamberparams frcmod.tip4pew loadoff solvents.lib loadoff anq.lib solvatebox ANQ TIP4PEWBOX 10.0 setbox ANQ vdw check ANQ set default PdbWriteCharges on savePDB ANQ anqWAT.pdb saveMol2 ANQ anqWAT.mol2 1 saveamberparm ANQ anqWAT.top anqWAT.crd remove ANQ ANQ.1 savePDB ANQ anqBOX.pdb saveMol2 ANQ anqBOX.mol2 1 saveamberparm ANQ anqBOX.top anqBOX.crd quit ********************************************************** Do they (delete remove) do exactly the same? 2) Then I face a new problem. The cartesian coordinates of the water molecules in the V0 and V1 states are exactly the same, BUT, the atom numbers and the the residue numbers in the solvent box have been shifted after the solute removal: New atom number = old atom number - number of atoms in the solute molecule Number of residue = old residue number - 1. And then, according to the amber12 manual for the Thermodynamic Integration I read: "The input (mdin) and starting coordinate files must be the same for the two groups. Furthermore, the two prmtop files must have the same number number of atoms, in the same order (since one common set of coordinates will be used for both.)..." [pg. 114] ... so... what can I do? 3) Thomas, could you please recommend me any reference about the 1- and 2- step transformations "in explicit solvent" and "in gas phase"? I am very curious about the practical and theoretical aspects of the transformation when in some cases you need the "in vacuo" step to close the cycle and in others no (regarding the charge removal) Actually I was thinking precisely in: > You dont need to do this anymore. The tranformation keeps internal > energies intact, so the dG for solute->nothing always equals 0.0. > > You do need to do the vacuum charge removal for the solute though. That > can be avoided by switching to one-step transformations if you want to. Sorry for being so persistent with this doubt. I send attached the PDB files after solvent removal and the library of the organic molecule in case is of any help.   Thank you very much!! Kepa K. ******************************************************* ________________________________ De: "steinbrt.rci.rutgers.edu" <steinbrt.rci.rutgers.edu> Para: AMBER Mailing List <amber.ambermd.org> Enviado: Martes 19 de marzo de 2013 11:07 Asunto: Re: [AMBER] Hydration Free Enery by TI Hi, > But this may be a problem when you are disappearing, say, one residue in a > stretch of RNA, since then you will end up with two separate pieces of > RNA. but, if you are removing an RNA base in the middle of a strand, you do want to end up with two separate pieces of DNA, dont you? Removing one base and subsequently closing the gap would be two big changes to simulate at once and I am not sure if a typical TI calculation would converge well on this. You could do it by deleting residue n in v1 and adding a (very long) bond between residues n-1 and n+1 at the appropriate place in leap. This transformation sounds like horribly difficult sampling to me, though. > So, I am wondering, is it possible to do in the VDW step what one usually > does for the charge, i.e, > > crgmask_v0=:10 > scmask_v0='' > crgmask_v1=:10 > scmask_v1=:10 > > while using the same prm and rst file for both V0 and V1? (in V0 you have > a > residue 10 without charge but with VDW, while in V1 you have residue 10 > without charge and without VDW). No, you cant set things up this way. The syntax you describe would make sense for disappearing residues, but that is not how the current code works. Kind Regards, Thomas Dr. Thomas Steinbrecher formerly at the BioMaps Institute Rutgers University 610 Taylor Rd. Piscataway, NJ 08854 _______________________________________________ AMBER mailing list AMBER.ambermd.org http://lists.ambermd.org/mailman/listinfo/amber

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber

Received on Tue Mar 19 2013 - 13:00:02 PDT
Custom Search