Spectroscopy and Quantum Chemical Studies

Universal Solvation Model Based on Solute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions

We present a new continuum solvation model based on the quantum mechanical charge density of a solute molecule interacting with a continuum description of the solvent. The model is called SMD, where the "D" stands for "density" to denote that the full solute electron density is used without defining partial atomic charges. "Continuum" denotes that the solvent is not represented explicitly but rather as a dielectric medium with surface tension at the solute-solvent boundary. SMD is a universal solvation model, where "universal" denotes its applicability to any charged or uncharged solute in any solvent or liquid medium for which a few key descriptors are known (in particular, dielectric constant, refractive index, bulk surface tension, and acidity and basicity parameters). The model separates the observable solvation free energy into two main components. The first component is the bulk electrostatic contribution arising from a self-consistent reaction field treatment that involves the solution of the nonhomogeneous Poisson equation for electrostatics in terms of the integral-equation-formalism polarizable continuum model (IEF-PCM). The cavities for the bulk electrostatic calculation are defined by superpositions of nuclear-centered spheres. The second component is called the cavity-dispersion-solvent-structure term and is the contribution arising from short-range interactions between the solute and solvent molecules in the first solvation shell. This contribution is a sum of terms that are proportional (with geometry-dependent proportionality constants called atomic surface tensions) to the solvent-accessible surface areas of the individual atoms of the solute. The SMD model has been parametrized with a training set of 2821 solvation data including 112 aqueous ionic solvation free energies, 220 solvation free energies for 166 ions in acetonitrile, methanol, and dimethyl sulfoxide, 2346 solvation free energies for 318 neutral solutes in 91 solvents (90 nonaqueous organic solvents and water), and 143 transfer free energies for 93 neutral solutes between water and 15 organic solvents. The elements present in the solutes are H, C, N, O, F, Si, P, S, Cl, and Br. The SMD model employs a single set of parameters (intrinsic atomic Coulomb radii and atomic surface tension coefficients) optimized over six electronic structure methods: M05-2X/MIDI!6D, M05-2X/6-31G, M05-2X/6-31+G, M05-2X/cc-pVTZ, B3LYP/6-31G, and HF/6-31G. Although the SMD model has been parametrized using the IEF-PCM protocol for bulk electrostatics, it may also be employed with other algorithms for solving the nonhomogeneous Poisson equation for continuum solvation calculations in which the solute is represented by its electron density in real space. This includes, for example, the conductor-like screening algorithm. With the 6-31G basis set, the SMD model achieves mean unsigned errors of 0.6-1.0 kcal/mol in the solvation free energies of tested neutrals and mean unsigned errors of 4 kcal/mol on average for ions with either Gaussian03 or GAMESS.

The Fluid Mosaic Model of the Structure of Cell Membranes

A fluid mosaic model is presented for the gross organization and structure of the proteins and lipids of biological membranes. The model is consistent with the restrictions imposed by thermodynamics. In this model, the proteins that are integral to the membrane are a heterogeneous set of globular molecules, each arranged in an amphipathic structure, that is, with the ionic and highly polar groups protruding from the membrane into the aqueous phase, and the nonpolar groups largely buried in the hydrophobic interior of the membrane. These globular molecules are partially embedded in a matrix of phospholipid. The bulk of the phospholipid is organized as a discontinuous, fluid bilayer, although a small fraction of the lipid may interact specifically with the membrane proteins. The fluid mosaic structure is therefore formally analogous to a two-dimensional oriented solution of integral proteins (or lipoproteins) in the viscous phospholipid bilayer solvent. Recent experiments with a wide variety of techniqes and several different membrane systems are described, all of which abet consistent with, and add much detail to, the fluid mosaic model. It therefore seems appropriate to suggest possible mechanisms for various membrane functions and membrane-mediated phenomena in the light of the model. As examples, experimentally testable mechanisms are suggested for cell surface changes in malignant transformation, and for cooperative effects exhibited in the interactions of membranes with some specific ligands. Note added in proof: Since this article was written, we have obtained electron microscopic evidence (69) that the concanavalin A binding sites on the membranes of SV40 virus-transformed mouse fibroblasts (3T3 cells) are more clustered than the sites on the membranes of normal cells, as predicted by the hypothesis represented in Fig. 7B. T-here has also appeared a study by Taylor et al. (70) showing the remarkable effects produced on lymphocytes by the addition of antibodies directed to their surface immunoglobulin molecules. The antibodies induce a redistribution and pinocytosis of these surface immunoglobulins, so that within about 30 minutes at 37 degrees C the surface immunoglobulins are completely swept out of the membrane. These effects do not occur, however, if the bivalent antibodies are replaced by their univalent Fab fragments or if the antibody experiments are carried out at 0 degrees C instead of 37 degrees C. These and related results strongly indicate that the bivalent antibodies produce an aggregation of the surface immunoglobulin molecules in the plane of the membrane, which can occur only if the immunoglobulin molecules are free to diffuse in the membrane. This aggregation then appears to trigger off the pinocytosis of the membrane components by some unknown mechanism. Such membrane transformations may be of crucial importance in the induction of an antibody response to an antigen, as well as iv other processes of cell differentiation.

Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models

Abstract An analytical algorithm, called SETTLE, for resetting the positions and velocities to satisfy the holonomic constraints on the rigid water model is presented. This method is still based on the Cartesian coordinate system and can be used in place of SHAKE and RATTLE. We implemented this algorithm in the SPASMS package of molecular mechanics and dynamics. Several series of molecular dynamics simulations were carried out to examine the performance of the new algorithm in comparison with the original RATTLE method. It was found that SETTLE is of higher accuracy and is faster than RATTLE with reasonable tolerances by three to nine times on a scalar machine. Furthermore, the performance improvement ranged from factors of 26 to 98 on a vector machine since the method presented is not iterative. © 1992 by John Wiley & Sons, Inc.

A full coupled-cluster singles and doubles model: The inclusion of disconnected triples

The coupled-cluster singles and doubles model (CCSD) is derived algebraically, presenting the full set of equations for a general reference function explicitly in spin–orbital form. The computational implementation of the CCSD model, which involves cubic and quartic terms, is discussed and results are reported and compared with full CI calculations for H2O and BeH2. We demonstrate that the CCSD exponential ansatz sums higher-order correlation effects efficiently even for BeH2, near its transition state geometry where quasidegeneracy efforts are quite large, recovering 98% of the full CI correlation energy. For H2O, CCSD plus the fourth-order triple excitation correction agrees with the full CI energy to 0.5 kcal/mol. Comparisons with low-order models provide estimates of the effect of the higher-order terms T1T2, T21T2, T31, and T41 on the correlation energy.

Efficacious Form for Model Pseudopotentials

A simple way has been discovered to put model pseudopotentials, $V(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})={\ensuremath{\Sigma}}_{\mathrm{lm}}|{Y}_{\mathrm{lm}}〉{V}_{l}(r)\ifmmode\times\else\texttimes\fi{}〈{Y}_{\mathrm{lm}}|$, into a form which reduces the number of integrals of $V(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$ required for an energyband calculation from $\frac{\mathrm{mn}(n+1)}{2}$ to $\mathrm{mn}$ for each $l$ in the sum (where $n$ is the number of plane waves used in the expansion and $m$ the number of points in the Brillouin zone at which the calculation is performed). The new form may be chosen to improve the accuracy of the pseudopotential when used in other chemical environments.

Interaction Models for Water in Relation to Protein Hydration

Classical Electrostatics in Biology and Chemistry

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.

Topics in Applied Physics

"Topics in Applied Physics." Nuclear Technology, 39(3), pp. 327–328

The future of quantum biology

Biological systems are dynamical, constantly exchanging energy and matter with the environment in order to maintain the non-equilibrium state synonymous with living. Developments in observational techniques have allowed us to study biological dynamics on increasingly small scales. Such studies have revealed evidence of quantum mechanical effects, which cannot be accounted for by classical physics, in a range of biological processes. Quantum biology is the study of such processes, and here we provide an outline of the current state of the field, as well as insights into future directions.