[AMBER] PMF calculations

From: Mahmoud Soliman <mahmoudelkot.gmail.com>
Date: Fri, 15 Apr 2011 20:06:11 +0200

   Dear users,
   I am doing umbrella sampling/wham calculation to get a free energy profile
   of some nucleophilic substitution reaction between the enzyme and substrate.
   I am using bond distance restraints from 3.2 A to 1.0 A. I use the same
   run_script provided by Wang here
   [1]http://projects.eml.org/mcm/people/wang/PMF/tutorial.html (below). My
   concern is that he uses a fixed starting coordinate file (btn.crd) for all
   simulations even at different restraint values. My understanding is,
   particularly with my case, if I need to put a distance restraint at value
   2.0 A and my starting coordinate (equilibrated starting structure) is about
   3.2A, I assume that is a big jump and it will be very hard to get a smooth
   SCF convergence. This will not be the case if I start from the previous
   window which is say 1.8 A rather than from 3.2 A. Am I right??
   The only difference between my case and his case is that he is using a set
   of dihedral angle restraints but I think the same concept applies
   #!/bin/tcsh
   set number=(-180 -160 -140 -120 -100 -80 -60 -40 -20 0 20 40 60 80 100 120
   140 160 180)
   foreach N ($number)
    sed -e 's/minimum/'$N'/g' umbrella.in > input.in
    $AMBERHOME/exe/sander -O -i input.in -o btn_$N.out -p btn.top -c btn.crd -r
   btn_$N.xyz -x btn_$N.trj
   end
   My script
   #!/bin/bash
   Number=(3.2 3.0 2.8 2.6 2.4 2.2 2.0 1.8 1.6 1.4 1.2 1.0)
   for i in ${Number}
   do
   sed -e 's/minimum/'$i'/g' umbrella.in > input.in
   /export/home/nmonama/amber/bin/sander -O -i input.in -o dist_$i.out -p
   com_solvated.top -c com_solvated_equil.qmmm.rst -r dist_$i.xyz -x
   dist_$i.trj
   done
   Thanks
   Best wishes
   Mahmoud

   --

   *************************************************

   Mahmoud E. Soliman

   Computational Chemistry & Modeling (PhD)

   Department of Chemistry

   University of Bath

   Bath

   BA2 7AY

   United Kingdom

   [2]http://people.bath.ac.uk/mess20/

   [3]http://www.bath.ac.uk/person/812559


   *********************************************

   Mahmoud E. Soliman

   Lecturer of pharmaceutical organic chemistry

   Pharmaceutical Organic Chemistry Dept.

   Faculty of pharmacy

   Zagazig University

   Zagazig

   Egypt

   **********************************************

   Email:

   [4]mess20.bath.ac.uk

   [5]meelkot.zu.edu.eg

   [6]mahmoudelkot.gmail.com

References

   1. http://projects.eml.org/mcm/people/wang/PMF/tutorial.html
   2. http://people.bath.ac.uk/mess20/
   3. http://www.bath.ac.uk/person/812559
   4. mailto:mess20.bath.ac.uk
   5. mailto:meelkot.zu.edu.eg
   6. mailto:mahmoudelkot.gmail.com
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Apr 15 2011 - 11:00:03 PDT
Custom Search