Zinc Protein Simulations Using the Cationic Dummy Atom (CaDA) Approach

  • Amber Force Field Parameters for Tetrahedron-Shaped Zinc Divalent Cation Used by the Cationic Dummy Atom Approach to Zinc Protein Simulations

    Table 1. The Most Recent Amber Force Field Parameters for Tetrahedron-Shaped Zinc Divalent Cation Used by the Cationic Dummy Atom (CaDA) Approach to Zinc Protein Simulations (Feb 25, 2005)

Table 1

Amber all-atom force field parameters for the CaDA approach using an explicit solvent model

bond K [ kcal/(mol Å2)] Req (Å)
DZ-ZN 640 0.90
DZ-DZ 640 1.470
angle K [kcal/(mol radian2)] Teq (deg.)
DZ-ZN-DZ 55 109.50
DZ-DZ-DZ 55 60.0
DZ-DZ-ZN 55 35.25
dihedral IDIVF Vn/2 (kcal/mol) γ(deg.) N
ZN-DZ-DZ-DZ 1 0 35.3 2
DZ-ZN-DZ-DZ 1 0 120.0 2
DZ-DZ-DZ-DZ 1 0 70.5 2
VDW mass r*(Å) Eps (kcal/mol) VDW mass r*(Å) Eps (kcal/mol)
ZN 53.38 3.1 1E-6 DZ 3.0 0 0

Table 2. Amber Atom Types and Charges of Histidinate, Hydroxide and Tetrahedron-Shaped Zinc Divalent Cation Used by the CaDA Approach

Atom name Atom type Charge
N N -0.5641
H H 0.2469
CA CT 0.3171
HA H1 0.0096
CB CT -0.1347
HB2 HC 0.0083
HB3 HC 0.0381
CG CC 0.1504
ND1 NA -0.7626
CE1 CR 0.4994
HE1 H5 -0.0295
NE2 NB -0.7656
CD2 CV 0.0405
HD2 H4 0.0525
C C 0.4588
O O -0.5653
H1 HO 0.2049
O OH -1.2049
Tetrahedron-Shaped Zinc Divalent Cation    
DZ DZ 0.5

Generating Topology and Coordinates Files Using xLeap (AmberTools V1.5)
Parameter and example files:

Download and copy the parameter files to respective directories (see commands below).

cp leaprc.zinc $AMBERHOME/dat/leap/cmd
cp znb.lib hin.lib hydroxide.lib $AMBERHOME/dat/leap/lib cp frcmod.zinc $AMBERHOME/dat/leap/parm

Launch xLeap by typing the following command in c shell of X11:
$AMBERHOME/exe/xleap –s –f $AMBERHOME/dat/leap/cmd/leaprc.ff99SB

Type the following commands in xLeap:
source $AMBERHOME/dat/leap/cmd/leaprc.zinc
m = loadpdb revised_zinc.pdb
addions m counterion 0
solvatebox m TIP3PBOX 10.0
saveamberparm m prmtop prmcrd

xLeap does not recognize $AMBERHOME and you need to replace it with an actual path to the amber directory (e.g., /Applications/amber11).

"revised_zinc.pdb" is a pdb file that has been revised as follows:

Change the residue and atom names of the zinc divalent cation to ZNB and Zn, respectively. ZNB is the residue name of the tetrahedron-shaped zinc divalent cation and has 4 peripheral dummy atoms attached to the central zinc atom. xLeap will automatically add these dummies if they are not defined in the pdb file.

Change the residue name of the zinc-coordinating water molecule to HO- (hydroxide).

Change HIS and CYS to HIN and CYM, respectively, if these residues coordinate the zinc ion. HIN and CYM are anionic forms of histidine and cysteine, respectively.

Change GLU and ASP to GLH and ASH, respectively, if these residues form a hydrogen bond with HIN, HO-, or CYM. GLH and ASH are the neutral forms of glutamic and aspartic acids, respectively. Swap the coordinates of OE1 and OE2 of GLH or ASH if OE1 is closer to the hydrogen bond acceptor of HIN, HO-, or CYM than OE2 (see below).

ATOM 1622 OE1 GLU C 94 30.263 40.009 33.982 1.00 31.93 O
ATOM 1623 OE2 GLU C 94 30.708 38.077 34.933 1.00 33.98 O1-
ATOM 1622 OE1 GLH C 94 30.708 38.077 34.933 1.00 33.98 O
ATOM 1623 OE2 GLH C 94 30.263 40.009 33.982 1.00 31.93 O1-

"counterion" is either Na+ or Cl-.

"prmtop" and "prmcrd" are topology and coordinate files, respectively, for sander and/or pmemd

Ignore the "4n unknown element" message, wherein n is the number of the zinc ion in the protein, after "m = loadpdb revised_zinc.pdb." This benign message is caused by the lack of dummy atoms defined in the pdb file.

The CaDA Approach References

  • Pang, Y.-P., Novel zinc protein molecular dynamics simulations: Steps toward antiangiogenesis for cancer treatment. J. Mol. Model. 1999, 5, 196-202.
  • Pang, Y.-P.; Xu, K.; El Yazal, J.; Prendergast, F. G., Successful molecular dynamics simulation of the zinc-bound farnesyltransferase using the cationic dummy atom approach. Protein Sci. 2000, 9, 1857-1865.
  • Pang, Y.-P., Successful molecular dynamics simulation of two zinc complexes bridged by a hydroxide in phosphotriesterase using the cationic dummy atom method. Proteins. 2001, 45, 183-189.
  • Oelschlaeger, P.; Schmid, R. D.; Pleiss, J., Insight into the mechanism of the IMP-1 metallo-beta-lactamase by molecular dynamics simulations. Protein Eng. 2003, 16, (5), 341-350.
  • Oelschlaeger, P.; Schmid, R. D.; Pleiss, J., Modeling domino effects in enzymes: Molecular basis of the substrate specificity of the bacterial metallo-beta-lactamases IMP-1 and IMP-6. Biochemistry 2003, 42, (30), 8945-8956.
  • Serotype-selective, small-molecule inhibitors of the zinc endopeptidase of botulinum neurotoxin serotype A. JG Park, PC Sill, EF Makiyi, AT Garcia-Sosa, CB Millard, JJ Schmidt, YP Pang, Bioorg. Med. Chem., 2006, 14, 395-408.
  • Computer-aided Lead Optimization: Improved Small-Molecule Inhibitor of the Zinc Endopeptidase of Botulinum Neurotoxin Serotype A. Jing Tang, Jewn Giew Park, Charles B. Millard, James J. Schmidt, and Yuan-Ping Pang, PLoS ONE, 2007, 2(8): e761.
  • Potent New Small-Molecule Inhibitor of Botulinum Neurotoxin Serotype A Endopeptidase Developed by Synthesis-Based Computer-Aided Molecular Design. Yuan-Ping Pang, Anuradha Vummenthala, Rajesh K. Mishra, Jewn Giew Park, Shaohua Wang, Jon Davis, Charles B. Millard, and James J. Schmidt, PLoS ONE, 2009, 4(11): e7730.
  • Small Molecules Showing Significant Protection of Mice against Botulinum Neurotoxin Serotype A. Yuan-Ping Pang, Jon Davis, Shaohua Wang, Jewn Giew Park, Madhusoodana P. Nambiar, J. Schmidt, and Charles B. Millard, PLoS ONE, 2010, 5(4): e10129.

Q & A

Q: Energy minimization failed independent of what zinc protein was used.

A: This is likely due to special characters in leaprc.zinc, hin.lib, znb.lib, hydroxide.lib, or frcmod.zinc. These characters are introduced after copying the files from one operating system to the other. If you initially download the files to a Windows system, you need to follow the procedure below to make these files work properly in a Unix or Linux system.

  1. Download and save the .zip files to a windows system
  2. Transfer the .zip files to a Unix or Linux system
  3. Change the ".zip" extension to ".gz"
  4. Unzip the .gz files with the "gunzip" command
  5. Add ".lib" extension to the hin, znb, and hydroxide files
  6. Add ".zinc" extension to the leaprc and frcmod files

Another possibility is that the residue or atom name is not labeled correctly.

Q: The four dummy atoms do not point towards their respective coordinates after energy minimization.

A: This happens when DIELC is not 1.0 as the CaCA parameters are developed for DIELC = 1.0.

Q: The hydroxide always fused with one of the dummy atoms leading to an infinite electrostatic energy, when minimizing a protein containing a hydroxide-bound, tetrahedron-shaped zinc divalent cation.

A: This happens when the atom type of the dummy atom is mistakenly defined as "H" as it should be "DZ."

Q: The xLeap procedure failed when using 2CAB.pdb (carbonic anhydrase form B). However, the same procedure worked when using 1G54.pdb (carbonic anhydrase II).

A: This problem is caused by the 4 lines after the zinc atom line in the 2CAB.pdb (see Figure 1). Simply delete these lines and the procedure should work with 2CAB.pdb as well.

Q: Can the CaDA parameters be used with the ff03 force field?

A: The CaDA parameters work with all AMBER all-atom force fields excluding implicit solvent models.