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.

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

A major revival in the use of classical electrostatics as an approach to the study of charged and polar molecules in aqueous solution has been made possible through the development of fast numerical and computational methods to solve the Poisson-Boltzmann equation for solute molecules that have complex shapes and charge distributions. Graphical visualization of the calculated electrostatic potentials generated by proteins and nucleic acids has revealed insights into the role of electrostatic interactions in a wide range of biological phenomena. Classical electrostatics has also proved to be successful quantitative tool yielding accurate descriptions of electrical potentials, diffusion limited processes, pH-dependent properties of proteins, ionic strength-dependent phenomena, and the solvation free energies of organic molecules.

https://science.sciencemag.org/content/sci/268/5214/1144.full.pdf

Classical Electrostatics in Biology and Chemistry

Publisher Summary This chapter discusses the integrated methods for the construction of three-dimensional models and computational probing of structure–function relations in G protein-coupled receptors (GPCR). The rapid pace of cloning and expression of G protein-coupled receptors offers attractive opportunities to probe the structural basis of signal transduction mechanisms at the level of these cell-surface receptors. Major insights have emerged from comparisons and classifications of the amino acid sequences of GPCRs into families defined by evolutionary developments and adapted to perform selective functions. Structural data on GPCRs, based on biochemical, immunological, and biophysical approaches have validated consensus architecture of GPCRs with an extracellular N-terminus, a cytoplasmic C-terminus, and a transmembrane portion comprised of seven-transmembrane helical domains connected by loops. Developments in the molecular modeling and computational exploration of GPCR proteins indicate a tantalizing potential to alleviate some of these difficulties. These expectations are based on the increased rate of success achieved by molecular modeling and computational simulation methods in providing structural insights relevant to the functions of biological molecules.

[19] Integrated methods for the construction of three-dimensional models and computational probing of structure-function relations in G protein-coupled receptors

Prediction of pH-dependent properties of proteins.

A new method of simultaneously scaling and shifting the Lennard‐Jones (LJ) potential in molecular dynamics (MD) and thermodynamic integration (TI) simulations is presented. The approach allows the smooth creation or annihilation of atoms or molecules in an ensemble of solvent molecules during a molecular simulation. By scaling and shifting the LJ potential in the direction of the interatomic distance between particles, the method eliminates the problem of the creation or annihilation of a large repulsive LJ potential at the initial or final state of a TI. The optimal degree of shifting and scaling the LJ potential as a function of a control variable λ was studied for the annihilation and creation of neon in aqueous solution. The method was further tested on the calculation of the free energy of aqueous solvation of a small molecule, ethanol. In contrast to linear scaling of the LJ potential during TI, the calculated free energies using the new separation‐shifted scaling approach are reasonably well converged after 200–500 ps of simulation and show smaller hysteresis comparing forward and reverse TI.

Separation‐shifted scaling, a new scaling method for Lennard‐Jones interactions in thermodynamic integration

Theoretical chemistry, as implemented on fast computers, is beginning to yield accurate predictions of the thermodynamic and kinetic properties of large molecular assemblies. In addition to providing detailed insights into the origins of molecular activity, theoretical calculations can be used to design new molecules with specific properties. This article describes two types of calculations that show special promise as design tools, the thermodynamic cycle-perturbation method and the Brownian reactive dynamics method. These methods can be applied to calculate equilibrium and rate constants that describe many aspects of molecular recognition, stability, and reactivity.

https://science.sciencemag.org/content/sci/238/4826/486.full.pdf

Computer-aided molecular design.

The enzyme acetylcholinesterase generates a strong electrostatic field that can attract the cationic substrate acetylcholine to the active site. However, the long and narrow active site gorge seems inconsistent with the enzyme's high catalytic rate. A molecular dynamics simulation of acetylcholinesterase in water reveals the transient opening of a short channel, large enough to pass a water molecule, through a thin wall of the active site near tryptophan-84. This simulation suggests that substrate, products, or solvent could move through this "back door," in addition to the entrance revealed by the crystallographic structure. Electrostatic calculations show a strong field at the back door, oriented to attract the substrate and the reaction product choline and to repel the other reaction product, acetate. Analysis of the open back door conformation suggests a mutation that could seal the back door and thus test the hypothesis that thermal motion of this enzyme may open multiple routes of access to its active site.

http://science.sciencemag.org/content/sci/263/5151/1276.full.pdf

Open "Back Door" in a Molecular Dynamics Simulation of Acetylcholinesterase

The electrostatic contributions to free energies of solvation of several small molecules have been calculated, treating the solvent as a statistical continuum. The computational method is based on solving the linearized Poisson-Boltzmann equation for the electrostatic potentials using the finite-difference scheme. A careful study of convergence indicates the importance of a fine grid spacing, as well as the short comings of rotational averaging. The computed free energies of solvation are in excellent agreement with the experimental results as well as the free energy perturbation calculations

J. A. McCammon

Papers

Prediction of pH-dependent properties of proteins.

Separation‐shifted scaling, a new scaling method for Lennard‐Jones interactions in thermodynamic integration

Computer-aided molecular design.

Open "Back Door" in a Molecular Dynamics Simulation of Acetylcholinesterase

Continuum model calculations of solvation free energies: accurate evaluation of electrostatic contributions