Darrin M. York

Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems

An N⋅log(N) method for evaluating electrostatic energies and forces of large periodic systems is presented. The method is based on interpolation of the reciprocal space Ewald sums and evaluation of the resulting convolutions using fast Fourier transforms. Timings and accuracies are presented for three large crystalline ionic systems.

CHARMM: The biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model‐building capabilities. The CHARMM program is applicable to problems involving a much broader class of many‐particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical‐molecular mechanical force fields, to all‐atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

A chemical potential equalization method for molecular simulations

A formulation of the chemical potential (electronegativity) equalization principle is presented from the perspective of density‐functional theory. The resulting equations provide a linear‐response framework for describing the redistribution of electrons upon perturbation by an applied field. The method has two main advantages over existing electronegativity equalization and charge equilibration methods that allow extension to accurate molecular dynamics simulations. Firstly, the expansion of the energy is taken about the molecular ground state instead of the neutral atom ground states; hence, in the absence of an external field, the molecular charge distribution can be represented by static point charges and dipoles obtained from fitting to high‐level ab initio calculations without modification. Secondly, in the presence of applied fields or interactions with other molecules, the density response can be modeled accurately using basis functions. Inclusion of basis functions with dipolar or higher order mul...

An Efficient Linear-Scaling Ewald Method for Long-Range Electrostatic Interactions in Combined QM/MM Calculations

A method is presented for the efficient evaluation of long-range electrostatic forces in combined quantum mechanical and molecular mechanical (QM/MM) calculations of periodic systems. The QM/MM-Ewald method is a linear-scaling electrostatic method that utilizes the particle mesh Ewald algorithm for calculation of point charge interactions of molecular mechanical atoms and a real-space multipolar expansion for the quantum mechanical electrostatic terms plus a pairwise periodic correction factor for the QM and QM/MM interactions that does not need to be re-evaluated during the self-consistent field procedure. The method is tested in a series of molecular dynamics simulations of the ion-ion association of ammonium chloride and ammonium metaphosphate and the dissociative phosphoryl transfer of methyl phosphate and acetyl phosphate. Results from periodic boundary molecular dynamics (PBMD) simulations employing the QM/MM-Ewald method are compared with corresponding PBMD simulations using electrostatic cutoffs and with results from nonperiodic stochastic boundary molecular dynamics (SBMD) simulations, with cutoffs and with full electrostatics (no cutoff). The present method allows extension of linear-scaling Ewald methods to molecular simulations of enzyme and ribozyme reactions that use combined QM/MM potentials.

Linear‐scaling semiempirical quantum calculations for macromolecules

A linear‐scaling method to carry out semiempirical quantum mechanical calculations for large systems has been developed based on the density matrix version of the divide‐and‐conquer approach. The method has been tested and demonstrated to be accurate and efficient. With this implementation, semiempirical quantum mechanical calculations are made possible for large molecules over 9000 atoms on a typical workstation. For biological macromolecules, solvent effects are included with a dielectric continuum model.

Specific Reaction Parametrization of the AM1/d Hamiltonian for Phosphoryl Transfer Reactions: H, O, and P Atoms

A semiempirical AM1/d Hamiltonian is developed to model phosphoryl transfer reactions catalyzed by enzymes and ribozymes for use in linear-scaling calculations and combined quantum mechanical/molecular mechanical simulations. The model, designated AM1/d-PhoT, is parametrized for H, O, and P atoms to reproduce high-level density-functional results from a recently constructed database of quantum calculations for RNA catalysis ( http://theory.chem.umn.edu/Database/QCRNA ), including geometries and relative energies of minima, transition states and reactive intermediates, dipole moments, proton affinities, and other relevant properties. The model is tested in the gas phase and in solution using a QM/MM potential. The results indicate that the method provides significantly higher accuracy than MNDO/d, AM1, and PM3 methods and, for the transphosphorylation reactions, is in close agreement with the density-functional calculations at the B3LYP/6-311++G(3df,2p) level with a reduction in computational cost of 3-4 orders of magnitude. The model is expected to have considerable impact on the application of semiempirical QM/MM methods to transphosphorylation reactions in solution, enzymes, and ribozymes and to ultimately facilitate the design of improved next-generation multiscale quantum models.

The fast Fourier Poisson method for calculating Ewald sums

The conventional Ewald expression for the electrostatic energy and forces is recast in a form that can be evaluated to high accuracy in order N log(N) steps using fast Fourier transforms. The fast Fourier Poisson method does not rely on interpolation approaches or Taylor/multipole expansions, and can be easily integrated with conventional molecular dynamics algorithms.

Accurate Proton Affinity and Gas-Phase Basicity Values for Molecules Important in Biocatalysis

Benchmark quantum calculations of proton affinities and gas-phase basicities of molecules relevant to biochemical processes, particularly acid/base catalysis, are presented and compared for a variety of multilevel and density functional quantum models. Included are nucleic acid bases in both keto and enol tautomeric forms, ribose in B-form and A-form sugar pucker conformations, amino acid side chains and backbone molecules, and various phosphates and phosphoranes, including thio substitutions. This work presents a high-level thermodynamic characterization of biologically relevant protonation states and provides a benchmark database for development of next-generation semiempirical and approximate density functional quantum models and parametrization of methods to predict pK(a) values and relative solvation energies.

Constant pH Replica Exchange Molecular Dynamics in Explicit Solvent Using Discrete Protonation States: Implementation, Testing, and Validation

By utilizing Graphics Processing Units, we show that constant pH molecular dynamics simulations (CpHMD) run in Generalized Born (GB) implicit solvent for long time scales can yield poor pKa predictions as a result of sampling unrealistic conformations. To address this shortcoming, we present a method for performing constant pH molecular dynamics simulations (CpHMD) in explicit solvent using a discrete protonation state model. The method involves standard molecular dynamics (MD) being propagated in explicit solvent followed by protonation state changes being attempted in GB implicit solvent at fixed intervals. Replica exchange along the pH-dimension (pH-REMD) helps to obtain acceptable titration behavior with the proposed method. We analyzed the effects of various parameters and settings on the titration behavior of CpHMD and pH-REMD in explicit solvent, including the size of the simulation unit cell and the length of the relaxation dynamics following protonation state changes. We tested the method with the amino acid model compounds, a small pentapeptide with two titratable sites, and hen egg white lysozyme (HEWL). The proposed method yields superior predicted pKa values for HEWL over hundreds of nanoseconds of simulation relative to corresponding predicted values from simulations run in implicit solvent.

Quantum mechanical/molecular mechanical simulation study of the mechanism of hairpin ribozyme catalysis.

The molecular mechanism of hairpin ribozyme catalysis is studied with molecular dynamics simulations using a combined quantum mechanical and molecular mechanical (QM/MM) potential with a recently developed semiempirical AM1/d-PhoT model for phosphoryl transfer reactions. Simulations are used to derive one- and two-dimensional potentials of mean force to examine specific reaction paths and assess the feasibility of proposed general acid and base mechanisms. Density-functional calculations of truncated active site models provide complementary insight to the simulation results. Key factors utilized by the hairpin ribozyme to enhance the rate of transphosphorylation are presented, and the roles of A38 and G8 as general acid and base catalysts are discussed. The computational results are consistent with available experimental data, provide support for a general acid/base mechanism played by functional groups on the nucleobases, and offer important insight into the ability of RNA to act as a catalyst without explicit participation by divalent metal ions.