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

/pdf/self-assembled-monolayers-of-thiolates-on-metals-as-a-form-1neb0hj3k4.pdf

Self-assembled monolayers of thiolates on metals as a form of nanotechnology.

Molecular mechanics is powerful for its speed in atomistic simulations, but an accurate force field is required. The Amber ff99SB force field improved protein secondary structure balance and dynamics from earlier force fields like ff99, but weaknesses in side chain rotamer and backbone secondary structure preferences have been identified. Here, we performed a complete refit of all amino acid side chain dihedral parameters, which had been carried over from ff94. The training set of conformations included multidimensional dihedral scans designed to improve transferability of the parameters. Improvement in all amino acids was obtained as compared to ff99SB. Parameters were also generated for alternate protonation states of ionizable side chains. Average errors in relative energies of pairs of conformations were under 1.0 kcal/mol as compared to QM, reduced 35% from ff99SB. We also took the opportunity to make empirical adjustments to the protein backbone dihedral parameters as compared to ff99SB. Multiple sm...

/pdf/ff14sb-improving-the-accuracy-of-protein-side-chain-and-2f0ci0lvo5.pdf

ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB

In just three years, the green fluorescent protein (GFP) from the jellyfish Aequorea victoria has vaulted from obscurity to become one of the most widely studied and exploited proteins in biochemistry and cell biology. Its amazing ability to generate a highly visible, efficiently emitting internal fluorophore is both intrinsically fascinating and tremendously valuable. High-resolution crystal structures of GFP offer unprecedented opportunities to understand and manipulate the relation between protein structure and spectroscopic function. GFP has become well established as a marker of gene expression and protein targeting in intact cells and organisms. Mutagenesis and engineering of GFP into chimeric proteins are opening new vistas in physiological indicators, biosensors, and photochemical memories.

/pdf/the-green-fluorescent-protein-3cl09aqrj2.pdf

The green fluorescent protein

The green fluorescent protein (GFP) of the jellyfish Aequorea Victoria has attracted widespread interest since the discovery that its chromophore is generated by the autocatalytic, posttranslational cyclization and oxidation of a hexapeptide unit. This permits fusion of the DNA sequence of GFP with that of any protein whose expression or transport can then be readily monitored by sensitive fluorescence methods without the need to add exogenous fluorescent dyes. The excited state dynamics of GFP were studied following photo-excitation of each of its two strong absorption bands in the visible using fluorescence upconversion spectroscopy (about 100 fs time resolution). It is shown that excitation of the higher energy feature leads very rapidly to a form of the lower energy species, and that the excited state interconversion rate can be markedly slowed by replacing exchangeable protons with deuterons. This observation and others lead to a model in which the two visible absorption bands correspond to GFP in two ground-state conformations. These conformations can be slowly interconverted in the ground state, but the process is much faster in the excited state. The observed isotope effect suggests that the initial excited state process involves a proton transfer reaction that is followed by additional structural changes. These observations may help to rationalize and motivate mutations that alter the absorption properties and improve the photo stability of GFP.

https://www.pnas.org/content/pnas/93/16/8362.full.pdf

Ultra-fast excited state dynamics in green fluorescent protein: multiple states and proton transfer.

The fusion and spreading of phospholipid bilayers on glass surfaces was investigated as a function of pH and ionic strength. Membrane fusion to the support was favorable at high ionic strength and low pH for vesicles containing a net negative charge; however, neutral and positively charged vesicles fused under all conditions attempted. This result suggests that van der Waals and electrostatic interactions govern the fusion process. Membrane spreading over a planar surface was favorable at low pH regardless of the net charge on the bilayer, and the process is driven by van der Waals forces. On the other hand membrane propagation is impeded at high pH or on highly curved surfaces. In this case a combination of hydration and bending interactions is primarily responsible for arresting the spreading process. These results provide a framework for understanding many of the factors that influence the effectiveness of scratches on planar supported bilayers as barriers to lateral diffusion and lead to a simple meth...

Formation and spreading of lipid bilayers on planar glass supports

Lithographically patterned grids of photoresist, aluminum oxide, or gold on oxidized silicon substrates were used to partition supported lipid bilayers into micrometer-scale arrays of isolated fluid membrane corrals. Fluorescently labeled lipids were observed to diffuse freely within each membrane corral but were confined by the micropatterned barriers. The concentrations of fluorescent probe molecules in individual corrals were altered by selective photobleaching to create arrays of fluid membrane patches with differing compositions. Application of an electric field parallel to the surface induced steady-state concentration gradients of charged membrane components in the corrals. In addition to producing patches of membrane with continuously varying composition, these gradients provide an intrinsically parallel means of acquiring information about molecular properties such as the diffusion coefficient in individual corrals.

https://science.sciencemag.org/content/sci/275/5300/651.full.pdf

Micropatterning fluid lipid bilayers on solid supports.

Stark spectroscopy has been applied to a wide range of molecular systems and materials. A generally useful method for obtaining electronic and vibrational Stark spectra that does not require sophisticated equipment is described. By working with frozen glasses it is possible to study nearly any molecular system, including ions and proteins. Quantitative analysis of the spectra provides information on the change in dipole moment and polarizability associated with a transition. The change in dipole moment reflects the degree of charge separation for a transition, a quantity of interest to a variety of fields. The polarizability change describes the sensitivity of a transition to an electrostatic field such as that found in a protein or an ordered synthetic material. Applications to donor-acceptor polyenes, transition metal complexes (metal-to-ligand and metal-to-metal mixed valence transitions), and nonphotosynthetic biological systems are reviewed.

Steven G. Boxer

Papers

Ultra-fast excited state dynamics in green fluorescent protein: multiple states and proton transfer.

Formation and spreading of lipid bilayers on planar glass supports

Micropatterning fluid lipid bilayers on solid supports.

STARK SPECTROSCOPY: Applications in Chemistry, Biology, and Materials Science

Model membrane systems and their applications