10. W. Kabsch, Acta Cryst., A32, 922 (1976), and A34, 827 (1978). 11. S. J. Weiner, P. A. Kollman, D. Case, P. Weiner, J. Am. Chem. Soc., 106, 765 (1984); (b) P. Weiner, P. A. Kollman, J. Comput. Chem., 2, 287 (1981). 12. A. Perry, I nt. J. Computer Math, Sec. B6, 327 (1978). 13. U. Burkert, N. L. Allinger, 11 Molecular Mechanics", ACS Monograph, 177, 21 (1982). 14. A. J. Irwin, J. B. Jones, J. Am. Chem. Soc., 98, 8476 (1976). 15. J. R. Miller, J. V. Beitz, R. K. Huddleston, / Chem. Soc.f 106, 5057 (1984). 16 F. M. Menger, F. F. Chow, H. Kaiserman, P. C. VesQuez, J. Am. Chem. Soc., 105, 4996 (1983). 17. (a) G. Wipff, A. Dearing, P. K. Weiner, J. M. Blaney, P.A. Kollman,/ Am. Chem. Soc.r 105, 997 (1983); (b) P. Kollman, Acc. Chem. Res., 18, 105 (1985). 1& M. W. Makinen, W. Maret, M. B. Yim, Proc. Natl. Acad. Sci. USA. 80, 2584 (1983).

The concept of a quantum chemical system (QCS) as a viable and useful tool in quantum chemical calculations requires that the system be continuously modified and expanded to handle new developments. The MOLE QCS, as originally conceived, was designed to handle the study of molecular electronic structure by means of the following types of calculations: (1) semiempirical (CNDO/2, Hückel, extended Hückel); (2) ab initio SCF; (3) configuration interaction (CI); (4) natural orbital analysis; and (5) expectation values (electronic moments and point properties). Since 1968, MOLE has been continuously modified to improve the efficiency of the algorithms which accomplish the above. Recently major additions have altered MOLE sufficiently to justify reporting its new capabilities. The additions are as follows: (1) extended CI capability; (2) perturbation theory CI for the choosing of optimal configurations; (3) configuration generation capability; (4) improved SCF proficiency; and (5) valence-only electronic structure capability based on the use of model potentials for core electrons. The present paper describes the current status of the MOLE system.

In molecular mechanics (MM) studies, atom types and/or bond types of molecules are needed to determine prior to energy calculations. We present here an automatic algorithm of perceiving atom types that are defined in a description table, and an automatic algorithm of assigning bond types just based on atomic connectivity. The algorithms have been implemented in a new module of the AMBER packages. This auxiliary module, antechamber (roughly meaning "before AMBER"), can be applied to generate necessary inputs of leap-the AMBER program to generate topologies for minimization, molecular dynamics, etc., for most organic molecules. The algorithms behind the manipulations may be useful for other molecular mechanical packages as well as applications that need to designate atom types and bond types.

The complexation of an adamantyl-phosphate derivative with one beta-cyclodextrin, with two beta-cyclodextrins, and with two beta-cyclodextrins dimerized with a disulfide bridge was studied by computational methods (MD, MM/PBSA, and MM/GBSA) to analyze and rationalize the chelate effect. Although this effect is usually explained by invoking favorable entropy contribution due to the preorganization of the ligand, it has been determined experimentally that in this case it is enthalpy-driven. The computational results are in accord with this finding, although the entropy contribution due to the solvent structural organization around the complex is crucial for the final estimates of the free energy of complexation.

We have developed a new-generation Amber united-atom force field for simulations involving highly demanding conformational sampling such as protein folding and protein-protein binding. In the new united-atom force field, all hydrogens on aliphatic carbons in all amino acids are united with carbons except those on C alpha. Our choice of explicit representation of all protein backbone atoms aims at minimizing perturbation to protein backbone conformational distributions and to simplify development of backbone torsion terms. Tests with dipeptides and solvated proteins show that our goal is achieved quite successfully. The new united-atom force field uses the same new RESP charging scheme based on B3LYP/cc-pVTZ//HF/6-31g** quantum mechanical calculations in the PCM continuum solvent as that in the Duan et al. force field. van der Waals parameters are empirically refitted starting from published values with respect to experimental solvation free energies of amino acid side-chain analogues. The suitability of mixing new point charges and van der Waals parameters with existing Amber covalent terms is tested on alanine dipeptide and is found to be reasonable. Parameters for all new torsion terms are refitted based on the new point charges and the van der Waals parameters. Molecular dynamics simulations of three small globular proteins in the explicit TIP3P solvent are performed to test the overall stability and accuracy of the new united-atom force field. Good agreements between the united-atom force field and the Duan et al. all-atom force field for both backbone and side-chain conformations are observed. In addition, the per-step efficiency of the new united-atom force field is demonstrated for simulations in the implicit generalized Born solvent. A speedup around two is observed over the Duan et al. all-atom force field for the three tested small proteins. Finally, the efficiency gain of the new united-atom force field in conformational sampling is further demonstrated with a well-known toy protein folding system, an 18 residue polyalanine in distance-dependent dielectric. The new united-atom force field is at least a factor of 200 more efficient than the Duan et al. all-atom force field for ab initio folding of the tested peptide.

The open nucleotide pocket conformation of actin in the profilin:actinCaATP x-ray structure has been hypothesized to be a crucial intermediate for nucleotide exchange in the actin depolymerization/polymerization cycle. The requirement for ancillary modification of actin for crystallization leads to ambiguities in this interpretation, however. We have used molecular dynamics simulations to model the thermodynamic properties of the actin x-ray structure, outside the crystal lattice, in an aqueous environment with profilin removed. Our simulations show that the open-nucleotide-pocket, profilin-free structure is actually unstable, and closes. The coordination of actin to the nucleotide in the molecular-dynamics-derived closed structure is virtually identical to that in the closed profilin:actinSrATP x-ray structure. Thus, there is currently no thermodynamically stable structure representing the open-nucleotide-pocket state of actin.

AbstractThe original article to which this Erratum refers was published in Journal of Computational Chemistry (2004) 25(9) 1157–1174

In this work, an efficient strategy was presented to search drug leads for human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) using hierarchical database screenings, which included a pharmacophore model, multiple-conformation rigid docking, solvation docking, and molecular mechanics-Poisson-Boltzmann/surface area (MM-PB/SA) sequentially. Encouraging results were achieved in searching a refined available chemical directory (ACD) database: the enrichment factor after the first three filters was estimated to be 25-fold; the hit rate for all the four filters was predicted to be 41% in a control test using 37 known HIV-1 non-nucleoside reverse transcriptase inhibitors; 10 out of 30 promising solvation-docking hits had MM-PB/SA binding free energies better than -6.8 kcal/mol and the best one, HIT15, had -17.0 kcal/mol. In conclusion, the hierarchical multiple-filter database searching strategy is an attractive strategy in drug lead exploration.

We have investigated the sampling efficiency in molecular dynamics with the PB implicit solvent when self-guiding forces are added. Compared with a high-temperature dynamics simulation, the use of self-guiding forces in room-temperature dynamics is found to be rather efficient as measured by potential energy fluctuation, gyration radius fluctuation, backbone RMSD fluctuation, number of unique clusters, and distribution of low RMSD structures over simulation time. Based on the enhanced sampling method, we have performed ab initio folding simulations of two small proteins, betabetaalpha1 and villin headpiece. The preliminary data for the folding simulations is presented. It is found that betabetaalpha1 folding proceeds by initiation of the turn and the helix. The hydrophobic collapse seems to be lagging behind or at most concurrent with the formation of the helix. The hairpin stability is weaker than the helix in our simulations. Its role in the early folding events seems to be less important than the more stable helix. In contrast, villin headpiece folding proceeds first by hydrophobic collapse. The formation of helices is later than the collapse phase, different from the betabetaalpha1 folding.

We describe here a general Amber force field (GAFF) for organic molecules. GAFF is designed to be compatible with existing Amber force fields for proteins and nucleic acids, and has parameters for most organic and pharmaceutical molecules that are composed of H, C, N, O, S, P, and halogens. It uses a simple functional form and a limited number of atom types, but incorporates both empirical and heuristic models to estimate force constants and partial atomic charges. The performance of GAFF in test cases is encouraging. In test I, 74 crystallographic structures were compared to GAFF minimized structures, with a root-mean-square displacement of 0.26 A, which is comparable to that of the Tripos 5.2 force field (0.25 A) and better than those of MMFF 94 and CHARMm (0.47 and 0.44 A, respectively). In test II, gas phase minimizations were performed on 22 nucleic acid base pairs, and the minimized structures and intermolecular energies were compared to MP2/6-31G* results. The RMS of displacements and relative energies were 0.25 A and 1.2 kcal/mol, respectively. These data are comparable to results from Parm99/RESP (0.16 A and 1.18 kcal/mol, respectively), which were parameterized to these base pairs. Test III looked at the relative energies of 71 conformational pairs that were used in development of the Parm99 force field. The RMS error in relative energies (compared to experiment) is about 0.5 kcal/mol. GAFF can be applied to wide range of molecules in an automatic fashion, making it suitable for rational drug design and database searching.

The prediction of pK a shifts of ionizable groups in proteins is of great relevance for a number of important biological phenomena. We present an implementation of the MM‐GBSA approach, which combines molecular mechanical (MM) and generalized Born (GB) continuum solvent energy terms, to the calculation of pK a values of a panel of nine proteins, including 69 individual comparisons with experiment. While applied so far mainly to the calculation of biomolecular binding free energies, we show that this method can also be used for the estimation of protein pK a shifts, with an accuracy around 1 pK a unit, even for strongly shifted residues. Our analysis reveals that the nonelectrostatic terms that are part of the MM‐GBSA free energy expression are important contributors to improved prediction accuracy. This suggests that most of the previous approaches that focus only on electrostatic interactions could be improved by adding other nonpolar energy terms to their free energy expression. Interestingly, our method yields best accuracy at protein dielectric constants of ε int = 2–4, which is in contrast to previous approaches that peak at higher ε int ≥ 8. An important component of our procedure is an intermediate minimization step of each protonation state involving different rotamers and tautomers as a way to explicitly model protein relaxation upon (de)protonation. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1865–1872, 2004

Development of the new generation of therapeutics against the influenza viral coat protein neuraminidase is a response to the continuing threat of influenza epidemics. A variety of structurally similar compounds have been reported that vary greatly in their ability to inhibit neuraminidase, a critical enzyme that cleaves sialic acid and promotes virion release. To determine how neuraminidase exhibits this wide range of affinities with structurally similar compounds, molecular dynamic simulations, coupled with free energy calculations, were used to determine the binding components of a series of neuraminidase inhibitors. Using four cocrystal structures of neuraminidase-inhibitor complexes, we examined the structural and energetic components of ligand potency and selectivity. An in-depth energetic analysis, including internal energy, entropy, and nonbonded interactions, reveals that potency of ligand binding is governed by nonpolar contacts. Electrostatic components generally oppose binding, although two of the best inhibitors use electrostatic interactions to orient the ligand. This investigation suggests that the enhanced selectivity and potency of the better ligands may arise from an improved positioning of their ligand atoms in the active site due to polar and hydrophobic functionalities. Simulations that included crystal water molecules in the active site indicate that the more potent ligands make less use of water-mediated interactions.

Direct hydroxide attack on the scissile carbonyl of the substrate has been suggested as a likely mechanism for esterase antibodies elicited by phosphonate haptens, which mimic the transition states for the alkaline hydrolysis of esters. 1 The unique amidase activity of esterase antibody 43C9 has been attributed to nucleophilic attack by an active‐site histidine residue. 2 Yet, the active site of 43C9 is strikingly similar to those of other esterase antibodies, particularly 17E8. We have carried out quantum mechanical calculations, molecular dynamics simulations, and free energy calculations to assess the mechanism involving direct hydroxide attack for 43C9. Results support this mechanism and suggest that the mechanism is plausible for other antiphosphonate antibodies that catalyze the hydrolysis of ( p ‐nitro)phenyl esters. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 1371–1377, 2003

AbstractFor Abstract see ChemInform Abstract in Full Text.

We have applied the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method (J. Srinivasan, T. E. Cheatham, P. Cieplak, P. A. Kollman, and D. A. Case, Journal of the American Chemical Society, 1998, Vol. 120, pp. 9401-9409) to study the interaction of an RNA aptamer with theophylline and its analogs. The MM-PBSA free energy analysis provides a reasonable absolute binding free energy for the RNA aptamer-theophylline complex formation. Energetic analysis reveals that the van der Waals interaction and the nonpolar contribution to solvation provide the basis for the favorable absolute free energy of complex. This trend is similar to other protein-ligand interactions studied previously. The MM-PBSA method also ranks the relative binding energies of five theophylline analogs approximately correctly, but not as well as the more conventional thermodynamic integration calculations, which were carried out to convert theophylline into its analogs. The comparison of MM-PBSA with TI suggests that the MM-PBSA method has some difficulties with the first-solvation-shell energetics.

Abstract This article discusses the basic physical chemical principles that determine the strength of non‐covalent interactions between ligands and receptors. We consider the basic thermodynamic principles involved in association and consider how molecular mechanics force fields can be helpful in understanding such phenomena. Three examples are discussed in some detail: biotin‐avidin, trimethoprim‐dihydrofolate reductase, and intercalator‐DNA interactions.

Molecular mechanics models have been applied extensively to study the dynamics of proteins and nucleic acids. Here we report the development of a third-generation point-charge all-atom force field for proteins. Following the earlier approach of Cornell et al., the charge set was obtained by fitting to the electrostatic potentials of dipeptides calculated using B3LYP/cc-pVTZ//HF/6-31G** quantum mechanical methods. The main-chain torsion parameters were obtained by fitting to the energy profiles of Ace-Ala-Nme and Ace-Gly-Nme di-peptides calculated using MP2/cc-pVTZ//HF/6-31G** quantum mechanical methods. All other parameters were taken from the existing AMBER data base. The major departure from previous force fields is that all quantum mechanical calculations were done in the condensed phase with continuum solvent models and an effective dielectric constant of epsilon = 4. We anticipate that this force field parameter set will address certain critical short comings of previous force fields in condensed-phase simulations of proteins. Initial tests on peptides demonstrated a high-degree of similarity between the calculated and the statistically measured Ramanchandran maps for both Ace-Gly-Nme and Ace-Ala-Nme di-peptides. Some highlights of our results include (1) well-preserved balance between the extended and helical region distributions, and (2) favorable type-II poly-proline helical region in agreement with recent experiments. Backward compatibility between the new and Cornell et al. charge sets, as judged by overall agreement between dipole moments, allows a smooth transition to the new force field in the area of ligand-binding calculations. Test simulations on a large set of proteins are also discussed.

We used classical molecular mechanics (MM) simulations and quantum mechanical (QM) structural relaxations to examine the active site of myosin when bound to ATP. Two conformations of myosin have been determined by x-ray crystallography. In one, there is no direct interaction between switch 2 and the nucleotide (open state). In the other (closed state), the universally conserved switch 2 glycine forms a hydrogen bond with a γ-phosphate oxygen. MM simulations indicate that the two states are thermodynamically stable and allow us to investigate the extent to which the P-loop, switch 1, and switch 2 are involved in hydrolysis. We find that the open structure has a higher affinity for ATP than the closed structure, and that ATP is distorted toward a transition state by interactions with the protein. We also examine how the structure of the binding site changes with either MgATP or CaATP as the nucleotide in myosin in the open conformer. Our analyses suggest that higher CaATPase rates occur because the leaving phosphate (Pi) group is more weakly bound and dissociation occurs faster. Finally, we validate the use of a particular formulation of a QM methodology (Car-Parrinello) to further refine the structures of the active site.