CHARMM (Chemistry at HARvard Macromolecular Mechanics) is a highly flexible computer program which uses empirical energy functions to model macromolecular systems. The program can read or model build structures, energy minimize them by first- or second-derivative techniques, perform a normal mode or molecular dynamics simulation, and analyze the structural, equilibrium, and dynamic properties determined in these calculations. The operations that CHARMM can perform are described, and some implementation details are given. A set of parameters for the empirical energy function and a sample run are included.

CHARMM: A program for macromolecular energy, minimization, and dynamics calculations

New protein parameters are reported for the all-atom empirical energy function in the CHARMM program. The parameter evaluation was based on a self-consistent approach designed to achieve a balance between the internal (bonding) and interaction (nonbonding) terms of the force field and among the solvent−solvent, solvent−solute, and solute−solute interactions. Optimization of the internal parameters used experimental gas-phase geometries, vibrational spectra, and torsional energy surfaces supplemented with ab initio results. The peptide backbone bonding parameters were optimized with respect to data for N-methylacetamide and the alanine dipeptide. The interaction parameters, particularly the atomic charges, were determined by fitting ab initio interaction energies and geometries of complexes between water and model compounds that represented the backbone and the various side chains. In addition, dipole moments, experimental heats and free energies of vaporization, solvation and sublimation, molecular volume...

/pdf/all-atom-empirical-potential-for-molecular-modeling-and-2nds0dxw12.pdf

All-atom empirical potential for molecular modeling and dynamics studies of proteins.

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecu- lar simulation program. It has been developed over the last three decades with a primary focus on molecules of bio- logical 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 estima- tors, 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 numer- ous 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.

/pdf/charmm-the-biomolecular-simulation-program-4kexem0xkx.pdf

CHARMM: the biomolecular simulation program.

Motivation: Molecular simulation has historically been a low-throughput technique, but faster computers and increasing amounts of genomic and structural data are changing this by enabling large-scale automated simulation of, for instance, many conformers or mutants of biomolecules with or without a range of ligands. At the same time, advances in performance and scaling now make it possible to model complex biomolecular interaction and function in a manner directly testable by experiment. These applications share a need for fast and efficient software that can be deployed on massive scale in clusters, web servers, distributed computing or cloud resources.

Results: Here, we present a range of new simulation algorithms and features developed during the past 4 years, leading up to the GROMACS 4.5 software package. The software now automatically handles wide classes of biomolecules, such as proteins, nucleic acids and lipids, and comes with all commonly used force fields for these molecules built-in. GROMACS supports several implicit solvent models, as well as new free-energy algorithms, and the software now uses multithreading for efficient parallelization even on low-end systems, including windows-based workstations. Together with hand-tuned assembly kernels and state-of-the-art parallelization, this provides extremely high performance and cost efficiency for high-throughput as well as massively parallel simulations.

Availability: GROMACS is an open source and free software available from http://www.gromacs.org.

Contact: erik.lindahl@scilifelab.se

Supplementary information:Supplementary data are available at Bioinformatics online.

Gromacs 4.5

For molecular dynamics simulations of hydrated proteins a simple yet reliable model for the intermolecular potential for water is required. Such a model must be an effective pair potential valid for liquid densities that takes average many-body interactions into account. We have developed a three-point charge model (on hydrogen and oxygen positions) with a Lennard-Jones 6–12 potential on the oxygen positions only. Parameters for the model were determined from 12 molecular dynamics runs covering the two-dimensional parameter space of charge and oxygen repulsion. Both potential energy and pressure were required to coincide with experimental values. The model has very satisfactory properties, is easily incorporated into protein-water potentials, and requires only 0.25 sec computertime per dynamics step (for 216 molecules) on a CRAY-1 computer.

Interaction Models for Water in Relation to Protein Hydration

The dynamics of a folded globular protein (bovine pancreatic trypsin inhibitor) have been studied by solving the equations of motion for the atoms with an empirical potential energy function. The results provide the magnitude, correlations and decay of fluctuations about the average structure. These suggest that the protein interior is fluid-like in that the local atom motions have a diffusional character.

Dynamics of folded proteins

A COMPLETE description of an enzyme requires a knowledge of its structure and the dynamics of its function. From the crystal structures of enzymes and enzyme–inhibitor complexes and the known chemistry of model systems, it has been possible in some cases to draw reasonable inferences concerning the mechanisms of enzyme-catalysed reactions. Little has been done so far, however, to supplement such essentially static results by an investigation of the reaction dynamics. This requires an understanding of the internal motions of the enzyme, as well as those of the substrate, since both are likely to be essential to the function. Here we present a theoretical study of a low frequency vibration involving the two globular lobes of lysozyme between which the cleft containing the active site is located1. Any motion involving this cleft could play a part in the catalytic activity; in fact, atom displacements of up to 0.75 A found in a comparison of the free enzyme and the enzyme-inhibitor complex indicate that the cleft has closed down in the latter2. The force constant for the low frequency bending vibration corresponding to the opening and closing of the cleft is obtained from empirical energy functions3. Because the protein surface moves appreciably during the vibration, damping effects resulting from the viscous dissipation in the solvent are included in the calculation4,5.

The hinge-bending mode in lysozyme

Side-chain torsional potentials in the bovine pancreatic trypsin inhibitor are calculated from empirical energy functions by use of the known X-ray structure of the protein and the rigid-geometry mapping technique. The potentials are analyzed to determine the roles and relative importance of contributions from the dipeptide backbone, the protein, and the crystalline environment of solvent and other protein molecules. The structural characteristics of the side chains determine two major patterns of energy surfaces, E(X1,X2): a gamma-branched pattern and a pattern for longer, straight side chains (Arg, Lys, Glu, and Met). Most of the dipeptide potential curves and surfaces have a local minimum corresponding to the side-chain torsional angles in the X-ray structure. Addition of the protein forces sharpens and/or selects from these minima, providing very good agreement with the experimental conformation for most side chains at the surface or in the core of the protein. Inclusion of the crystalline environment produces still better results, especially for the side chains extending away from the protein. The results are discussed in terms of the details of the interactions due to the surrounding, calculated solvent-accessibility figures and the temperature factors derived from the crystallographic refinement of the pancreatic trypsin inhibitor.

Side-chain torsional potentials: effect of dipeptide, protein, and solvent environment.

A reaction path is presented by which the effects of oxygen binding in hemoglobin are transmitted from a heme group to the surface of its subunit. Starting from the known deoxy geometry, it is shown by calculations with empirical energy functions and comparisons with available data how the change in heme geometry on ligation introduces a perturbation that leads to the tertiary structural alterations essential for cooperatively. It is found that there is little strain on the unliganded heme; instead, the reduced oxygen affinity of hemoglobin results from the strain on the liganded subunit in a tetramer with the deoxy quarternary structure.

https://www.pnas.org/content/pnas/74/3/801.full.pdf

Mechanism of tertiary structural change in hemoglobin.

Conformational potentials of sidechains in the bovine pancreatic trypsin inhibitor have been studied with an empirical energy function. Calculated minimumenergy positions are in excellent agreement with the x-ray structure for sidechains in the core or at the surface of the protein; as expected, angles for sidechains that are directed out into the solvent do not agree with the calculated values. The contributions to the potentials are analyzed and compared with the potentials for the free amino acid. Although there is a large restriction in the available conformational space due to nonbonded interactions, the minimum energy positions in the protein are close to those of the free amino acid; the significance of this result is discussed. To estimate the effective barriers for rotation of the aromatic rings (tyrosine and phenylalanine), calculations are done in which the protein is permitted to relax as a function of the ring orientation. Thr resulting barriers, which are much lowere than the rigid rotation barriers, are used to evaluate the rotation rates; comparison is made with the available nuclear magnetic resonance data.

https://www.pnas.org/content/pnas/72/6/2002.full.pdf

Bruce R. Gelin

Papers

Dynamics of folded proteins

The hinge-bending mode in lysozyme

Side-chain torsional potentials: effect of dipeptide, protein, and solvent environment.

Mechanism of tertiary structural change in hemoglobin.

Sidechain torsional potentials and motion of amino acids in porteins: bovine pancreatic trypsin inhibitor