Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

抗原变异可使得多种致病微生物易于逃避宿主免疫应答。表达在感染红细胞表面的恶性疟原虫红细胞表面蛋白1（PfPMP1）与感染红细胞、内皮细胞、树突状细胞以及胎盘的单个或多个受体作用，在黏附及免疫逃避中起关键的作用。每个单倍体基因组var基因家族编码约60种成员，通过启动转录不同的var基因变异体为抗原变异提供了分子基础。

疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A

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.

One of the most pervasive problems afflicting people throughout the world is inadequate access to clean water and sanitation. Problems with water are expected to grow worse in the coming decades, with water scarcity occurring globally, even in regions currently considered water-rich. Addressing these problems calls out for a tremendous amount of research to be conducted to identify robust new methods of purifying water at lower cost and with less energy, while at the same time minimizing the use of chemicals and impact on the environment. Here we highlight some of the science and technology being developed to improve the disinfection and decontamination of water, as well as efforts to increase water supplies through the safe re-use of wastewater and efficient desalination of sea and brackish water.

Science and technology for water purification in the coming decades

Department of Materials Science, University of Patras, 26504 Rio Patras, Greece, Theoretical and Physical Chemistry Institute, National Hellenic Research Foundation, 48 Vass. Constantinou Avenue, 116 35 Athens, Greece, Institut de Biologie Moleculaire et Cellulaire, UPR9021 CNRS, Immunologie et Chimie Therapeutiques, 67084 Strasbourg, France, and Dipartimento di Scienze Farmaceutiche, Universita di Trieste, Piazzale Europa 1, 34127 Trieste, Italy

Chemistry of Carbon Nanotubes

We use molecular dynamics simulations to study osmotically driven transport of water molecules through hexagonally packed carbon nanotube membranes. Our simulation setup comprises two such semipermeable membranes separating compartments of pure water and salt solution. The osmotic force drives water flow from the pure-water to the salt-solution compartment. Monitoring the flow at molecular resolution reveals several distinct features of nanoscale flows. In particular, thermal fluctuations become significant at the nanoscopic length scales, and as a result, the flow is stochastic in nature. Further, the flow appears frictionless and is limited primarily by the barriers at the entry and exit of the nanotube pore. The observed flow rates are high (5.8 water molecules per nanosecond and nanotube), comparable to those through the transmembrane protein aquaporin-1, and are practically independent of the length of the nanotube, in contrast to predictions of macroscopic hydrodynamics. All of these distinct characteristics of nanoscopic water flow can be modeled quantitatively by a 1D continuous-time random walk. At long times, the pure-water compartment is drained, and the net flow of water is interrupted by the formation of structured solvation layers of water sandwiched between two nanotube membranes. Structural and thermodynamic aspects of confined water monolayers are studied.

https://www.pnas.org/content/pnas/100/18/10175.full.pdf

Osmotic water transport through carbon nanotube membranes

Water molecules confined to nonpolar pores and cavities of nanoscopic dimensions exhibit highly unusual properties. Water filling is strongly cooperative, with the possible coexistence of filled and empty states and sensitivity to small perturbations of the pore polarity and solvent conditions. Confined water molecules form tightly hydrogen-bonded wires or clusters. The weak attractions to the confining wall, combined with strong interactions between water molecules, permit exceptionally rapid water flow, exceeding expectations from macroscopic hydrodynamics by several orders of magnitude. The proton mobility along 1D water wires also substantially exceeds that in the bulk. Proteins appear to exploit these unusual properties of confined water in their biological function (e.g., to ensure rapid water flow in aquaporins or to gate proton flow in proton pumps and enzymes). The unusual properties of water in nonpolar confinement are also relevant to the design of novel nanofluidic and molecular separation devices or fuel cells.

/pdf/water-in-nonpolar-confinement-from-nanotubes-to-proteins-and-412uecrowv.pdf

Water in Nonpolar Confinement: From Nanotubes to Proteins and Beyond *

Proteins can be denatured by pressures of a few hundred MPa. This finding apparently contradicts the most widely used model of protein stability, where the formation of a hydrophobic core drives protein folding. The pressure denaturation puzzle is resolved by focusing on the pressure-dependent transfer of water into the protein interior, in contrast to the transfer of nonpolar residues into water, the approach commonly taken in models of protein unfolding. Pressure denaturation of proteins can then be explained by the pressure destabilization of hydrophobic aggregates by using an information theory model of hydrophobic interactions. Pressure-denatured proteins, unlike heat-denatured proteins, retain a compact structure with water molecules penetrating their core. Activation volumes for hydrophobic contributions to protein folding and unfolding kinetics are positive. Clathrate hydrates are predicted to form by virtually the same mechanism that drives pressure denaturation of proteins.

https://www.pnas.org/content/pnas/95/4/1552.full.pdf

The pressure dependence of hydrophobic interactions is consistent with the observed pressure denaturation of proteins.

A molecular model of poorly understood hydrophobic effects is heuristically developed using the methods of information theory. Because primitive hydrophobic effects can be tied to the probability of observing a molecular-sized cavity in the solvent, the probability distribution of the number of solvent centers in a cavity volume is modeled on the basis of the two moments available from the density and radial distribution of oxygen atoms in liquid water. The modeled distribution then yields the probability that no solvent centers are found in the cavity volume. This model is shown to account quantitatively for the central hydrophobic phenomena of cavity formation and association of inert gas solutes. The connection of information theory to statistical thermodynamics provides a basis for clarification of hydrophobic effects. The simplicity and flexibility of the approach suggest that it should permit applications to conformational equilibria of nonpolar solutes and hydrophobic residues in biopolymers.

https://www.pnas.org/content/pnas/93/17/8951.full.pdf

An information theory model of hydrophobic interactions.

A theoretical approach is developed to quantify hydrophobic hydration and interactions on a molecular scale, with the goal of insight into the molecular origins of hydrophobic effects. The model is based on the fundamental relation between the probability for cavity formation in bulk water resulting from molecular-scale density fluctuations and the hydration free energy of the simplest hydrophobic solutes, hard particles. This probability is estimated using an information theory (IT) approach, incorporating experimentally available properties of bulk water:  the density and radial distribution function. The IT approach reproduces the simplest hydrophobic effects:  hydration of spherical nonpolar solutes, the potential of mean force (PMF) between methane molecules, and solvent contributions to the torsional equilibrium of butane. Applications of this approach to study temperature and pressure effects provide new insights into the thermodynamics and kinetics of protein folding. The IT model relates the hydr...

Shekhar Garde

Papers

Osmotic water transport through carbon nanotube membranes

Water in Nonpolar Confinement: From Nanotubes to Proteins and Beyond *

The pressure dependence of hydrophobic interactions is consistent with the observed pressure denaturation of proteins.

An information theory model of hydrophobic interactions.

Hydrophobic Effects on a Molecular Scale