[AMBER] Inputs QM/MM MD and MMGBSA

From: Soumendranath Bhakat <bhakatsoumendranath.gmail.com>
Date: Mon, 12 May 2014 01:27:16 +0530

Dear Amberists;

I was just testing some inputs for QM/MM MD for a ligand bound protein
complex. Here are some inputs for different stages of MD as well as MMGBSA.
Just for clarification of MD please let me know if these inputs are correct
or not.

Full_Mini.in

Initial minimization of QMMM: solvent molecules and added ions
 &cntrl
  imin = 1,
  maxcyc = 2500,
  ncyc = 750,
  ntb = 1,
  ntr = 1,
  cut = 12.0,
  ifqnt = 1,
 /
&qmmm
qmmask= '250,238.75,86,89,75',
qmcharge=0,
qmtheory=1,
qmshake=1,
qm_ewald=1,
qm_pme=1,
/

Hold the Protein fixed
10.0
RES 1 557

Heating.in


































*Heating Step of QMMM: stage-2 &cntrl imin= 0, irest=0, NTX=1, ntb=
1, NTPR=500, NTWX=500, NTWR=500, ntr=1, Tempi=0.0, Temp0=300.0,
ifqnt=1, NTT=3, gamma_ln=1.0, NTC=2, NTF=2, cut= 12.0,
nstlim=2500, dt=0.002,/&qmmmqmmask=
'250,238.75,86,89,75',qmcharge=0,qmtheory=1,qmshake=1,qm_ewald=1,qm_pme=1,/Keep
Protein and inhibitor fixed with weak restraints10.0RES 1 557ENDEND*

equil.in

Equilibration Step of MMP3 (MMMM): stage-1
 &cntrl
  imin= 0,
  irest=1,
  NTX=7,
  ntb=2,
  ntp=1,
  PRES0=1.0,
  TAUP=2.0,
  NTPR=500,
  NTWX=500,
  ntr=0,
  Tempi=300.0,
  Temp0=300.0,
  NTT=3,
  ifqnt=1,
  gamma_ln=1.0,
  NTC=2,
  NTF=2,
  cut=12.0,
  nstlim=250000,
  dt=0.002
 /
&qmmm
qmmask= '250,238.75,86,89,75',
qmcharge=0,
qmtheory=1,
qmshake=1,
qm_ewald=1,
qm_pme=1,
/
*md.in <http://md.in>*
Equilibration Step of MMP3 (MMMM): stage-1
 &cntrl
  imin= 0,
  irest=1,
  NTX=7,
  ntb=2,
  ntp=1,
  PRES0=1.0,
  TAUP=2.0,
  NTPR=500,
  NTWX=500,
  ntr=0,
  Tempi=300.0,
  Temp0=300.0,
  NTT=3,
  ifqnt=1,
  gamma_ln=1.0,
  NTC=2,
  NTF=2,
  cut=12.0,
  nstlim=2500000,
  dt=0.002
 /
&qmmm
qmmask= '250,238.75,86,89,75',
qmcharge=0,
qmtheory=1,
qmshake=1,
qm_ewald=1,
qm_pme=1,
/

*MMGBSA script*

Input file for running PB and GB in serial
&general
 startframe=1, endframe=5000, interval=5verbose=2,entropy=1,keep_file=0,
/
&gb
igb=5,saltcon=0.15,ifqnt=1,qmmask= '250,238.75,86,89,75',
qmcharge=0,qmtheory=1,
/

Please confirm once if these scripts are well to perform QM/MM MD at PM3
level and calculate post MD QM/MM MMGBSA.

Thanks a lot in advance

-- 
Thanks & Regards;
Soumendranath Bhakat
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun May 11 2014 - 13:00:02 PDT
Custom Search