Quantum simulation of ferrocytochrome c

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.

Simulation as experiment: a philosophical reassessment for biological modeling

Some scientific modelers suggest that complex simulation models that mimic biological processes should have a limited place in ecological and evolutionary studies. However, complex simulation models can have a role that is different from that of simpler models that are designed to be fit to data. Simulation can be viewed as another kind of experimental system and should be analyzed as such. Here, I argue that current discussions in the philosophy of science and in the physical sciences fields about the use of simulation as an experimental system have important implications for biology, especially complex sciences such as evolution and ecology. Simulation models can be used to mimic complex systems, but unlike nature, can be manipulated in ways that would be impossible, too costly or unethical to do in natural systems. Simulation can add to theory development and testing, can offer hypotheses about the way the world works and can give guidance as to which data are most important to gather experimentally.

The statistical-thermodynamic basis for computation of binding affinities: a critical review

Although the statistical thermodynamics of noncovalent binding has been considered in a number of theoretical papers, few methods of computing binding affinities are derived explicitly from this underlying theory. This has contributed to uncertainty and controversy in certain areas. This article therefore reviews and extends the connections of some important computational methods with the underlying statistical thermodynamics. A derivation of the standard free energy of binding forms the basis of this review. This derivation should be useful in formulating novel computational methods for predicting binding affinities. It also permits several important points to be established. For example, it is found that the double-annihilation method of computing binding energy does not yield the standard free energy of binding, but can be modified to yield this quantity. The derivation also makes it possible to define clearly the changes in translational, rotational, configurational, and solvent entropy upon binding. It is argued that molecular mass has a negligible effect upon the standard free energy of binding for biomolecular systems, and that the cratic entropy defined by Gurney is not a useful concept. In addition, the use of continuum models of the solvent in binding calculations is reviewed, and a formalism is presented for incorporating a limited number of solvent molecules explicitly.

A priori evaluation of aqueous polarization effects through Monte Carlo QM-MM simulations

A Monte Carlo quantum mechanical-molecular mechanical (QM-MM) simulation method was used to determine the contributions of the solvent polarization effect to the total interaction energies between solute and solvent for amino acid side chains and nucleotide bases in aqueous solution. In the present AM1-TIP3P approach, the solute molecule is characterized by valence electrons and nucleus cores with Hartree-Fock theory incorporating explicit solvent effects into the total Hamiltonian, while the solvent is approximated by the three-point charge TIP3P model. The polarization energy contributes 10 to 20 percent of the total electrostatic energy in these systems. The performance of the hybrid AM1-TIP3P model was further validated by consideration of bimolecular complexes with water and by computation of the free energies of solvation of organic molecules using statistical perturbation theory. Excellent agreement with ab initio 6-31G(d) results and experimental solvation free energies was obtained.