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

Early Steps of Supported Bilayer Formation Probed by Single Vesicle Fluorescence Assays

Synthetic lipid-oligonucleotide conjugates inserted into lipid vesicles mediate fusion when one population of vesicles displays the 5'-coupled conjugate and the other the 3'-coupled conjugate, so that anti-parallel hybridization allows the membrane surfaces to come into close proximity. Improved assays show that lipid mixing proceeds more quickly and to a much greater extent than content mixing, suggesting the latter is rate limiting. To test the effect of membrane-membrane spacing on fusion, a series of conjugates was constructed by adding 2-24 noncomplementary bases at the membrane-proximal ends of two complementary sequences. Increasing linker lengths generally resulted in progressively reduced rates and extents of lipid and content mixing, in contrast to higher vesicle docking rates. The relatively flexible, single-stranded DNA linker facilitates docking but allows greater spacing between the vesicles after docking, thus making the transition into fusion less probable, but not preventing it altogether. These experiments demonstrate the utility of DNA as a model system for fusion proteins, where sequence can easily be modified to systematically probe the effect of distance between bilayers in the fusion reaction.

https://www.pnas.org/content/pnas/106/4/979.full.pdf

Effects of linker sequences on vesicle fusion mediated by lipid-anchored DNA oligonucleotides.

Lateral variations in membrane composition are postulated to play a central role in many cellular events, but it has been difficult to probe membrane composition and organization on length scales of tens to hundreds of nanometers. We present a high-resolution imaging secondary ion mass spectrometry technique to reveal the lipid distribution within a phase-separated membrane with a lateral resolution of approximately 100 nanometers. Quantitative information about the chemical composition within small lipid domains was obtained with the use of isotopic labels to identify each molecular species. Composition variations were detected within some domains.

https://science.sciencemag.org/content/sci/313/5795/1948.full.pdf

Phase Separation of Lipid Membranes Analyzed with High-Resolution Secondary Ion Mass Spectrometry

Brownian ratchets use a time-varying asymmetric potential that can be applied to separate diffusing particles or molecules. A new type of Brownian ratchet, a geometrical Brownian ratchet, has been realized. Charged, fluorescently labeled phospholipids in a two-dimensional fluid bilayer were driven in one direction by an electric field through a two-dimensional periodic array of asymmetric barriers to lateral diffusion fabricated from titanium oxide on silica. Diffusion spreads the phospholipid molecules in the orthogonal direction, and the asymmetric barriers rectify the Brownian motion, causing a directional transport of molecules. The geometrical ratchet can be used as a continuous molecular sieve to separate mixtures of membrane-associated molecules that differ in electrophoretic mobility and diffusion coefficient.

https://science.sciencemag.org/content/sci/285/5430/1046.full.pdf

Brownian Ratchets: Molecular Separations in Lipid Bilayers Supported on Patterned Arrays

Water-soluble organic fluorophores are widely used as labels in biological systems. However, in many cases these fluorophores can interact strongly with lipid bilayers, influencing the interaction of the target with the bilayer and/or leading to misleading fluorescent signals. Here, we quantify the interaction of 32 common water-soluble dyes with model lipid bilayers to serve as an additional criterion when selecting a dye label.

/pdf/choose-your-label-wisely-water-soluble-fluorophores-often-4ylmpfxn8y.pdf

Steven G. Boxer

Papers

Early Steps of Supported Bilayer Formation Probed by Single Vesicle Fluorescence Assays

Effects of linker sequences on vesicle fusion mediated by lipid-anchored DNA oligonucleotides.

Phase Separation of Lipid Membranes Analyzed with High-Resolution Secondary Ion Mass Spectrometry

Brownian Ratchets: Molecular Separations in Lipid Bilayers Supported on Patterned Arrays

Choose your label wisely: water-soluble fluorophores often interact with lipid bilayers.