Development of the CHARMM Force Field for Lipids

The development of the CHARMM additive all-atom lipid force field (FF) is traced from the early 1990's to the most recent version (C36) published in 2010. Though simulations with early versions yielded useful results, they failed to reproduce two important quantities: a zero surface tension at the experimental bilayer surface area, and the signature splitting of the deuterium order parameters in the glycerol and upper chain carbons. Systematic optimization of parameters based on high level quantum mechanical data and free energy simulations have resolved these issues, and bilayers with a wide range of lipids can be simulated in tensionless ensembles using C36. Issues associated with other all-atom lipid FFs, success and limitations in the C36 FF and ongoing developments are also discussed.

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.

Human aldehyde dehydrogenase catalytic activity and structural interactions with coenzyme analogs

K(m) and V(max) values for 10 coenzyme analogs never previously studied with any aldehyde dehydrogenase and NADP(+) were compared with those for NAD(+) for three human aldehyde dehydrogenases (EC 1.2.1.3); the cytoplasmic E1 (the product of the aldh1 gene), the mitochondrial E2 (the product of the aldh2 gene) and the cytoplasmic E3 (the product of the aldh9 gene) isozymes. Structural information on changes in coenzyme-protein interactions were obtained via molecular dynamics (MD) studies with the E2 isozyme and quantum mechanical (QM) calculations were used to study changes in charge distribution of the pyridine ring and relative free energies of solvation of the purine ring in the analogs. E1 showed the broadest substrate specificity and was the only isozyme subject to substrate inhibition, both of which are suggested to be due to the two coenzyme conformations observed previously in the sheep crystal structure. NADP(+) selectivity is indicated to be influenced by Glu195 in E1 and E2. Substitutions in the purine ring affected K(m) but not V(max), with the changes in K(m) being dominated by the hydrophobicity of the purine ring as indicted by the QM calculations. Substitutions in the pyridine ring sometimes rendered the coenzymes inactive, with no consistent pattern observed for the three coenzymes. Structural analysis of the coenzyme analog-E2 MD simulations revealed different structural perturbations of the surrounding active site, though interactions with Asn169 and Glu399 were preserved in all cases.