AMBER Archive (2009)

Subject: Re: [AMBER] GAFF and NME, ACE terminal res

From: Jio M (jiomm_at_yahoo.com)
Date: Wed Jul 15 2009 - 23:29:08 CDT

Dear AMBER users,
 
Thanks Francois for reply .
 
one thing more I want to know that, I didnt mention clearly earlier (sorry for that), can I use the constraint charges in first RESP fitting mentioned in tutorial for NMe and ACE i.e.
 
NMe-------->
 
N -.4157
H .2719
C -.1490
H .0976
H .0976
H .0976
 net charge 0
similarly for ACE
 
after this the charges so derived on residues I want to use the residues to be later saved as prmtop file using GAFF.
 
thanks
 
JIomm                      

 
 
 
 

> I am new to AMBER though I was going through Advanced Amber tutorial  Tutorial 1 using dye and linker.

Be aware there are several inconsistencies in this tutorial.

> I want to make same type of residues to be protected by NMe (for C  terminal) or ACE(N terminal). pls clarify the following:
>
> 1) can I use GAFF for these residues to be used later and making prmtop file.

This depends on the structure you are going to study. If you want to develop a fragment // a component of a protein or a nucleic acid (i. e. a modified nucleotide or unusual amino acid), you should not use GAFF & you should use ff99SB (or ff03; for proteins only).

If your structure is organic you could indeed use GAFF.

> 2) I cannot use other than HF/631-G* method to optimise. the charges
> that I will use, during RESP fit, for NMe and ACE (as given in
> tutorial) are actually derived by optimising amino acids by HF/631-G*
> method also.

Just forget it: the HF/6-31G* or HF/6-31G** theory level (TL) was used in the Amber force field topology databse construction. You do not need MP2 correction or DFT in the particular cases presented.

> but for GAFF parameters have been derived by using MP2/631-G (will this
> effect the charges so derived by HF method during simulation). this is
> confusion pls clarify it whether I can use GAFF or not or  I have
> misunderstood the things

Do not mix up the TL used in force field paremeter developement & that used in charge derivation. Here, we talk about charge derivation and using HF/6-31G*//HF/6-31G* (or HF/6-31G*//HF/6-31G**) is perfect for GAFF or ff99SB; [for ff03 you need to use B3LYP/cc-pVTZ SCRF(IEFPCM,solvent=Ether)//HF/6-31G**]

For ff99SB or more generally ff94, please read:
Cieplak et al. Application of the multimolecule and multiconformational RESP methodology to biopolymers: Charge derivation for DNA, RNA, and proteins. J. Comput. Chem. 1995, 16, 1357.

For ff03, see:
Duan et al. A Point-Charge Force Field for Molecular Mechanics Simulations of Proteins Based on Condensed-Phase Quantum Mechanical Calculations. J. Comput. Chem. 2003, 24, 1999.

Finally, I suggest you to read as well:
http://q4md-forcefieldtools.org/Tutorial/ where charge derivation is carried out accordingly to what was done in the Amber force field topology database (Cieplak et al. J. Comput. Chem. 1995, 16, 1357).

regards, Francois

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