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. In addition, 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 paper provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM paper in 1983.

Abstract:A hybrid quantum mechanical/molecular mechanical (QM/MM) potential energy function with HartreeFock, density functional theory (DFT), and post-HF (RIMP2, MP2, CCSD) capability has been implemented in the CHARMM and Q-Chem software packages. In addition, we have modified CHARMM and Q-Chem to take advantage of the newly introduced replica path and the nudged elastic band methods, which are powerful techniques for studying reaction pathways in a highly parallel (i.e., parallel/parallel) fashion, with each pathway point being distributed to a different node of a large cluster. To test our implementation, a series of systems were studied and comparisons were made to both full QM calculations and previous QM/MM studies and experiments. For instance, the differences between HF, DFT, MP2, and CCSD QM/MM calculations of H2O· · · H2O, H2O· · · Na + , and H2O· · · Cl − complexes have been explored. Furthermore, the recently implemented polarizable Drude water model was used to make comparisons to the popular TIP3P and TIP4P water models for doing QM/MM calculations. We have also computed the energetic profile of the chorismate mutase catalyzed Claisen rearrangement at various QM/MM levels of theory and have compared the results with previous studies. Our best estimate for the activation energy is 8.20 kcal/mol and for the reaction energy is −23.1 kcal/mol, both calculated at the MP2/6-31+G(d)//MP2/6-31+G(d)/C22 level of theory.

Two new techniques for modeling chemical processes in condensed phases with combined quantum mechanical and molecular mechanical (QM/MM) potentials are introduced and tested on small, model compounds. The first technique, the double link atom (DLA) method, is an extension of the traditional, single link atom (SLA) method to avoid some of the problems with the latter method. These problems are primarily electrostatic, as the SLA method can produce an unphysical overall charge or dipole. The second technique, the delocalized Gaussian MM charge (DGMM) method, is an empirical way to include the delocalized character of the electron density of atoms in the MM region. This can be important for the electrostatic interaction of the QM region with nearby atoms in the MM region, and it can simplify the rules governing which classical interactions are included in the energies and forces. Even for very short distances, the DGMM method does not require the neglect of the MM host in the QM calculation. The DGMM method can be used for modeling reactions in solution, and it can be combined with methods such as the link atom, frozen orbital, or pseudopotential methods for terminating the QM region at a covalent bond. The DLA and the DGMM methods have been combined effectively. Presented here are tests on small, model systems that mimic properties important for reactions in proteins, in particular rotational barriers, proton affinities, and deprotonation energies. The new methods yield improved energetics for model compounds, vis-à-vis a point-MM-charge and SLA treatment.

rnThe application of hybrid quantum mechanical and molecular mechanical (QM/MM) potentials to the study of chemical reactions in enzymes is outlined. The discussion is general and addresses the difficulties encountered in an enzyme QM/MM study. First, general criteria for determining whether a particular enzyme is an appropriate candidate for a QM/MM approach are outlined. Methods for obtaining starting structures are detailed. The importance of choosing appropriate levels of ab initio or semiempirical theory is emphasized. Approaches for interfacing the QM and MM regions are briefly discussed, with greater detail given to describing our CHARMM-GAMESS interface. Techniques for partitioning the system into QM and MM regions are explored. Link atom placement, as distant from reacting atoms as possible within the confines of computational efficiency, is examined in some detail. Methods for determining reaction paths are also discussed.

Distributed Replica (REPDSTR) is a powerful parallelization technique enabling simulations of a group of replicas in a parallel/parallel fashion, where each replica is distributed to different nodes of a large cluster [Theor. Chem. Acc. 109: 140 (2003)]. Here, we use the framework provided by REPDSTR to combine a staged free energy perturbation protocol with different values of the thermodynamic coupling parameters with replica-exchange molecular dynamics (FEP/REMD. The structure of REPDSTR, which allows multiple parallel input/output (I/O), facilitates the treatment of replica-exchange to couple the N window simulations. As a result, each of the N synchronous window simulations benefit from the sampling carried out by the N-1 others. As illustrative examples of the FEP/REMD strategy, calculations of the absolute hydration and binding free energy of small molecules were performed using the biomolecular simulation program CHARMM adapted for the IBM Blue Gene/P platform. The computations show that a FEP/REMD strategy significantly improves the sampling and accelerate the convergence of absolute free energy computations.

A hybrid quantum mechanical−molecular mechanical (QM−MM) potential energy function with ab initio and density functional capabilities has been implemented in the CHARMM program. It makes use of the quantum mechanical program CADPAC and the CHARMM molecular mechanics energy function; a GAMESS(US) interface to the CHARMM program was already available. To test the methodology, a series of relatively small systems are studied and comparisons are made of full QM calculations with those from various QM−MM partitions. Both density functional and Hartree−Fock calculations for the quantum region are presented and, where possible, compared with results from previous AM1−MM calculations. For the density functional based QM−MM calculations, the LDA and BLYP functionals were used. The performances of both the density functional and Hartree−Fock based QM−MM calculations compare well with pure quantum calculations. The link atom method was tested by performing a number of QM−MM simulations on the complexes of metal cations with model ligands of biological interest. It was found that it gave good results for the structures, binding energies, and charge distributions.

A new web portal for the CHARMM macromolecular modeling package, CHARMMing (CHARMM interface and graphics, http://www.charmming.org), is presented. This tool provides a user friendly interface for the preparation, submission, monitoring, and visualization of molecular simulations (i.e., energy minimization, solvation, and dynamics). The infrastructure used to implement the web application is described. Two additional programs have been developed and integrated with CHARMMing: GENRTF, which is employed to define structural features not supported by the standard CHARMM force field, and a job broker, which is used to provide a portable method for using grid and cluster computing with CHARMMing. The use of the program is described with three proteins: 1YJP, 1O1O and 1UFY. Source code is provided allowing CHARMMing to be downloaded, installed, and used by supercomputing centers and research groups that have a CHARMM license. Although no software can replace a scientist’s own judgment and experience, CHARMMing eases the introduction of newcomers to the molecular modeling discipline by providing a graphical method for running simulations.

No abstract