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.

A new program package, XEASY, was written for interactive computer support of the analysis of NMR spectra for three-dimensional structure determination of biological macromolecules. XEASY was developed for work with 2D, 3D and 4D NMR data sets. It includes all the functions performed by the precursor program EASY, which was designed for the analysis of 2D NMR spectra, i.e., peak picking and support of sequence-specific resonance assignments, cross-peak assignments, cross-peak integration and rate constant determination for dynamic processes. Since the program utilizes the X-window system and the Motif widget set, it is portable on a wide range of UNIX workstations. The design objective was to provide maximal computer support for the analysis of spectra, while providing the user with complete control over the final resonance assignments. Technically important features of XEASY are the use and flexible visual display of 'strips', i.e., two-dimensional spectral regions that contain the relevant parts of 3D or 4D NMR spectra, automated sorting routines to narrow down the selection of strips that need to be interactively considered in a particular assignment step, a protocol of resonance assignments that can be used for reliable bookkeeping, independent of the assignment strategy used, and capabilities for proper treatment of spectral folding and efficient transfer of resonance assignments between spectra of different types and different dimensionality, including projected, reduced-dimensionality triple-resonance experiments.

A new adaptive umbrella sampling technique for molecular dynamics simulations is described. The high efficiency of the technique renders multidimensional adaptive umbrella sampling possible and thereby enables uniform sampling of the conformational space spanned by several degrees of freedom. The efficiency is achieved by using the weighted histogram analysis method to combine the results from different simulations, by a suitable extrapolation scheme to define the umbrella potential for regions that have not been sampled, and by a criterion to identify simulations during which the system was not in equilibrium. The technique is applied to two test systems, the alanine dipeptide and the threonine dipeptide, to sample the configurational space spanned by one or two dihedral angles. The umbrella potentials applied at the end of each adaptive umbrella sampling run are equal to the negative of the corresponding potentials of mean force. The trajectories obtained in the simulations can be used to calculate dynamical variables that are of interest. An example is the distribution of the distance between the HN and the Hβ proton that can be important for the interpretation of NMR experiments. Factors influencing the accuracy of the calculated quantities are discussed. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1450–1462, 1997

No abstract