Molecular dynamics simulation on a parallel computer.

For the purpose of molecular dynamics simulations of large biopolymers we have built a parallel computer with a systolic loop architecture, based on Transputers as computational units, and have programmed it in Occam 11. The computational nodes of the computer are linked together in a systolic ring. The program based on this .topology for large biopolymers increases its computational throughput nearly linearly with the number of computational nodes. The program developed is closely related to the simulation programs CHARMM and XPLOR, the input files required (force field, protein structure file, coordinates) and output files generated (sets of atomic coordinates representing dynamic trajectories and energies) are compatible with the corresponding files of these programs. Benchmark results of simulations of biopolymers comprising 66, 568, 3 634, 5 797 and 12 637 atoms are compared with XPLOR simulations on conventional computers (Cray, Convex, Vax). These results demonstrate that the software and hardware developed provide extremely cost effective biopolymer simulations. We present also a simulation (equilibrium of X-ray structure) of the complete photosynthetic reaction center of Rhodopseudomonus viridis (12 637 atoms). The simulation accounts for the Coulomb forces exactly, i.e. no cut-off had been assumed

Choice of fluorophore affects dynamic DNA nanostructures

The ability to dynamically remodel DNA origami structures or functional nanodevices is highly desired in the field of DNA nanotechnology. Concomitantly, the use of fluorophores to track and validate the dynamics of such DNA-based architectures is commonplace and often unavoidable. It is therefore crucial to be aware of the side effects of popular fluorophores, which are often exchanged without considering the potential impact on the system. Here, we show that the choice of fluorophore can strongly affect the reconfiguration of DNA nanostructures. To this end, we encapsulate a triple-stranded DNA (tsDNA) into water-in-oil compartments and functionalize their periphery with a single-stranded DNA handle (ssDNA). Thus, the tsDNA can bind and unbind from the periphery by reversible opening of the triplex and subsequent strand displacement. Using a combination of experiments, molecular dynamics (MD) simulations, and reaction-diffusion modelling, we demonstrate for 12 different fluorophore combinations that it is possible to alter or even inhibit the DNA nanostructure formation—without changing the DNA sequence. Besides its immediate importance for the design of pH-responsive switches and fluorophore labelling, our work presents a strategy to precisely tune the energy landscape of dynamic DNA nanodevices

Orientation effects in reactions of allenyl cations with styrene

Allenyl cations, generated from allenyl or alkynyl halides and Ag+, attack styrene at the side chain or at the aromatic nucleus. The allenyl/alkynyl product ratio is dependent on the structure of the precursor halide except for highly substituted systems

Escape from a metastable well under a time-ramped force

Thermally activated escape of an over-damped particle from a metastable well under the action of a time-ramped force is studied. We express the mean first passage time (MFPT) as the solution to a partial differential equation, which we solve numerically for a model case. We discuss two approximations of the MFPT, one of which works remarkably well over a wide range of loading rates, while the second is easy to calculate and can provide a valuable first estimate.Comment: 9 pages, including 2 figure

Membrane protein sequestering by ionic protein–lipid interactions.

Neuronal exocytosis is catalyzed by the SNARE protein syntaxin-1A(1). Syntaxin-1A is clustered in the plasma membrane at sites where synaptic vesicles undergo exocytosis(2,3). However, how syntaxin-1A is sequestered is unknown. Here, we show that syntaxin clustering is mediated by electrostatic interactions with the strongly anionic lipid phosphatidylinositol-4,5-bisphosphate (PIP2). We found with super-resolution STED microscopy on the plasma membrane of PC12 cells that PIP2 is the dominant inner-leaflet lipid in ~73 nm-sized microdomains. This high accumulation of PIP2 was required for syntaxin-1A sequestering, as destruction of PIP2 by the phosphatase synaptojanin-1 reduced syntaxin-1A clustering. Furthermore, co-reconstitution of PIP2 and the C-terminal part of syntaxin-1A in artificial giant unilamellar vesicles resulted in segregation of PIP2 and syntaxin-1A into distinct domains even when cholesterol was absent. Our results demonstrate that electrostatic protein-lipid interactions can result in the formation of microdomains independent of cholesterol or lipid phases

Action-derived molecular dynamics in the study of rare events

We present a practical method to generate classical trajectories with fixed initial and final boundary conditions. Our method is based on the minimization of a suitably defined discretized action. The method finds its most natural application in the study of rare events. Its capabilities are illustrated by non-trivial examples. The algorithm lends itself to straightforward parallelization, and when combined with molecular dynamics (MD) it promises to offer a powerful tool for the study of chemical reactions.Comment: 7 Pages, 4 Figures (3 in color), submitted to Phys. Rev. Let

Inter-subunit coupling enables fast CO2-fixation by reductive carboxylases

Enoyl-CoA carboxylases/reductases (ECRs) are some of the most efficient CO2-fixing enzymes described to date. However, the molecular mechanisms underlying the extraordinary catalytic activity of ECRs on the level of the protein assembly remain elusive. Here we used a combination of ambient-temperature X-ray free electron laser (XFEL) and cryogenic synchrotron experiments to study the structural organization of the ECR from Kitasatospora setae. The K. setae ECR is a homotetramer that differentiates into a pair of dimers of open- and closed-form subunits in the catalytically active state. Using molecular dynamics simulations and structure-based mutagenesis, we show that catalysis is synchronized in the K. setae ECR across the pair of dimers. This conformational coupling of catalytic domains is conferred by individual amino acids to achieve high CO2-fixation rates. Our results provide unprecedented insights into the dynamic organization and synchronized inter- and intrasubunit communications of this remarkably efficient CO2-fixing enzyme during catalysis.