6.2.2. Definition of Effective Properties 3064 6.3. Response Properties to Magnetic Fields 3066 6.3.1. Nuclear Shielding 3066 6.3.2. Indirect Spin−Spin Coupling 3067 6.3.3. EPR Parameters 3068 6.4. Properties of Chiral Systems 3069 6.4.1. Electronic Circular Dichroism (ECD) 3069 6.4.2. Optical Rotation (OR) 3069 6.4.3. VCD and VROA 3070 7. Continuum and Discrete Models 3071 7.1. Continuum Methods within MD and MC Simulations 3072

/pdf/quantum-mechanical-continuum-solvation-models-3edv963g14.pdf

Quantum mechanical continuum solvation models.

We describe the development, current features, and some directions for future development of the Amber package of computer programs. This package evolved from a program that was constructed in the late 1970s to do Assisted Model Building with Energy Refinement, and now contains a group of programs embodying a number of powerful tools of modern computational chemistry, focused on molecular dynamics and free energy calculations of proteins, nucleic acids, and carbohydrates.

/pdf/the-amber-biomolecular-simulation-programs-e2nyq1ef44.pdf

The Amber biomolecular simulation programs

Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects

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.

AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules

Reformulation et extension d'une methode d'estimation de la difference d'entropie de configuration entre les conformations moleculaires de proteines a partir de simulations de dynamique moleculaire

Evaluation of the configurational entropy for proteins: application to molecular dynamics simulations of an α-helix

The dielectric properties of proteins are central to their stability and activity. We use the Frohlich-Kirkwood theory of dielectrics to analyze two 1-ns molecular dynamics simulations of ferro- and ferricytochrome c in spherical droplets of 1400 water molecules. Protein and solvent are idealized as a series of concentric, spherical, dielectric media. Analysis results depend strongly on the treatment of the charged protein side chains at the protein/solvent interface. If charged side chains are viewed as part of the protein medium, then the protein dipole fluctuations are dominated by large, mutually uncorrelated, anisotropic, motions of the charged side chains. It is then incorrect to view the protein region as a single, homogeneous dielectric material. If one does take this view, estimates of the protein "dielectric constant" vary from 16 to 37, depending on the exact choice of model parameters. In contrast, if the charged portions of the charged side chains are viewed as part of the solvent medium, then theory and simulation are consistent: the protein dipole fluctuations excluding charged side chains are roughly those of a homogeneous, isotropic dielectric medium, with a dielectric constant of 4.7 +/- 1.0 (ferro) or 3.4 +/- 1.0 (ferri), in agreement with powder experiments. Statistical uncertainty and sensitivity to model parameters are small. Analysis of the radial dependence of the dipole fluctuations suggests that the inner half of the protein has a somewhat lower dielectric constant of 1.5-2, consistent with its biological function in electron transfer. These results suggest that Poisson-Boltzmann models could treat the protein bulk as a low-dielectric medium and the charged surface groups as part of the solvent region.

Internal and Interfacial Dielectric Properties of Cytochrome c from Molecular Dynamics in Aqueous Solution

Interdomain motion in liver alcohol dehydrogenase. Structural and energetic analysis of the hinge bending mode.

Microscopic Dielectric Properties of Cytochrome c from Molecular Dynamics Simulations in Aqueous Solution

The mean square amplitudes of atomic fluctuations for a polypeptide (decaglycine) α-helix evaluated from molecular dynamics simulations at seven temperatures between 5 and 300 K are compared with analytic harmonic results and with experimental values. Above 100 K the harmonic approximation significantly underestimates the amplitudes of the displacements. Analysis of the time dependence of the fluctuations shows that low-frequency modes ( 10 ps). Quantum corrections to the amplitude of the fluctuations are found to be small above 50 K. The mean square amplitudes obtained from the molecular dynamics simulations are compared with the values derived from x-ray temperature (Debye-Waller) factors for metmyoglobin (80, 250, and 300 K) and ferrocytochrome c (300 K).

David Perahia

Papers

Evaluation of the configurational entropy for proteins: application to molecular dynamics simulations of an α-helix

Internal and Interfacial Dielectric Properties of Cytochrome c from Molecular Dynamics in Aqueous Solution

Interdomain motion in liver alcohol dehydrogenase. Structural and energetic analysis of the hinge bending mode.

Microscopic Dielectric Properties of Cytochrome c from Molecular Dynamics Simulations in Aqueous Solution

Molecular dynamics of an α-helical polypeptide: Temperature dependence and deviation from harmonic behavior