Systolic Loop Methods for Molecular Dynamics Simulation,Generalised for Macromolecules

Abstract

Systolic loop programs have been shown to be very efficient for molecular dynamics simulations of liquid systems on distributed memory parallel computers. The original methods address the case where the number of molecules simulated exceeds the number of processors used. Simulations of large flexible molecules often do not meet this condition, requiring the three- and four-body terms used to model chemical bonds within a molecule to be distributed over several processors. This paper discusses how the systolic loop methods may be generalised to accommodate such systems, and describes the implementation of a computer program for simulation of protein dynamics. Performance figures are given for this program running typical simulations on a Meiko Computing Surface using different number of processors.     

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 II. 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 Rhodopseudomonas viridis (12 637 atoms). The simulation accounts for the Coulomb forces exactly, i.e. no cut-off had been assumed.     

Mobile Macromolecules in Plant Development

Plant cells transmit developmental signals and distribute nutrients via dynamic intercellular channels termed plasmodesmata (PD). Multidisciplinary inquiries have provided evidence that plasmodesmatal regulation is critical to fundamental plant functions, such as development, host–pathogen interactions, and systemic RNA silencing. This review focuses on macromolecules that transport cell-to-cell through PD and describes their implications on plant development.     

Study of Liquid Dimethyl Sulfoxide by Computer Simulation

The paper reports Monte Carlo and molecular dynamicsresults for pure liquid dimethyl sulfoxide (DMSO) at298 K and 1 atm. The classical 6–12 Lennard–Jones plusCoulomb pairwise potential was used to calculateintermolecular interaction energy. Potentialparameters for the liquid were optimized in this work.Some thermodynamic and dynamical properties obtained,such as heat of vaporization, density and diffusioncoefficient, are in good agreement with theexperimental values. The present model is comparedwith other models for DMSO reported previously. It isshown to be an improvement over earlier potentials.The structure factors and the radial distributionfunctions (rdf), are compared with experimentalresults for the liquid. The analysis shows that thestructure of DMSO is not completely understood yet anddeserves deeper investigation. The geometry of thedimer that corresponds to the rdf plots obtained, isreported. The results suggest that the dipole momentof this dimer plays an important role in the structureof the liquid.     

A Molecular Dynamics Computer Simulation Performance Comparison of Java Versus C

Abstract

The relative performance of molecular dynamics (MD) computer simulations of fluids written in ANSI C is compared to that achieved by a comparable program written in Java. The performance of the Java program is shown to be dependent upon its runtime environment. The Java Runtime Environment (JRE) from the Java Development Kit (JDK) 1.2 provides a Just-In-Time (JIT) compiler option on Solaris and Windows 95 platforms which decreases the execution time by approximately 4–10× compared to the standard Java interpreter. The compiled Java implementation of the MD computer simulation runs between 30–100% slower. depending on the platform, compared to the equivalent C implementation. The stability of the two simulations, as measured by conservation of energy is shown to be identical to within ～ 1% over 10⁵ time steps.     

Molecular Dynamic Simulation of Chaperonin-Mediated Protein Folding

We use molecular dynamics methods to simulate chaperonin-mediated refolding of barnase. A chaperonin term is added to the force field in order to simulate the hydrophobic environment in the central cavity of the chaperonins. Two aspects of our simulation results are consistent with experiments: (1) The hydrophobic environment of the central cavity of the chaperonin is an advantageous condition for the refolding of the misfolded intermediates. (2) One cycle of binding and release is not enough for the successful folding. Chaperonin-assisted protein folding maybe a procedure of multiple cycles of binding and release from the chaperonin.     

Diffusive Dynamics in a Detailed Potential: Application to Biological Macromolecules

Abstract

The local dynamics of macromolecules is obtained to second-order in the mode-coupling expansion of the Smoluchowski diffusion theory. The NMR spin-lattice relaxation times of different ¹³C or ¹⁵N nuclei along the chains are calculated and compared to experimental data from the literature. The macromolecules are considered as fluctuating 3D structures undergoing rotational diffusion. The fluctuations can be evaluated with any technique for sampling the configurational space. In the presented test cases Molecular Dynamics simulations have been applied to a DNA fragment and to the NK-2 homeodomain. In the case of the double-stranded DNA fragment d(TpCpGpCpG)₂, second and even first order theories are found to be in close agreement with experimental results. The major advantage of the diffusion technique is that only a good statistics is important as input while the solvent dynamic effects enter through hydrodynamic theory. Application based on Hybrid Monte Carlo schemes coupled with J-walking, are now in progress.     

Projection Methods for the Analysis of Complex Motions in Macromolecules

Abstract

In studies of macromolecular dynamics it is often desirable to analyze complex motions in terms of a small number of coordinates. Only for simple types of motion, e.g., rigid-body motions, these coordinates can be easily constructed from the Cartesian atomic coordinates. This article presents an approach that is applicable to infinitesimal or approximately infinitesimal motions, e.g., Cartesian velocities, normal modes, or atomic fluctuations. The basic idea is to characterize the subspace of interesting motions by a set of (possibly linearly dependent) vectors describing elementary displacements, and then project the dynamics onto this subspace. Often the elementary displacements can be found by physical intuition. The restriction to small displacements facilitates the study of complicated coupled motions and permits the construction of collective-motion subspaces that do not correspond to any set of generalized coordinates.

As an example for this technique, we analyze the low-frequency normal modes of proteins up to ≈ 20 THz (600 cm^?1) in order to see what kinds of motions occupy which frequency range. This kind of analysis is useful for the interpretation of spectroscopic measurements on proteins, e.g., inelastic neutron scattering experiments.     

Computer Simulation by Molecular Dynamics as a Tool for Modelling of Molecular Systems

Abstract

A survey is given of methods for simulation of molecular systems on a computer. The various assumptions, approximations and limitations are discussed and the possibility of making comparisons with experimental quantities is assessed. Finally, a number of practical applications of molecular dynamics simulation techniques in chemistry are reviewed.     

The Lennard-Jones Fluid Revisited: Computer Simulation Results

Abstract

Pseudoexperimental data of high accuracy on the pressure and the internal energy of the Lennard-Jones fluid have been generated both by the Monte Carlo and molecular dynamics methods for five subcritical and three supercritical isotherms. Values of the chemical potential of the Lennard-Jones fluid computed by a new version of the gradual insertion particle method for two isotherms up to very high densities are also reported and discussed, and compared with existing data.     

Computer Simulation of Liquid Methanol II. System Size Effects

Abstract

A series of molecular dynamics simulations of liquid methanol has been carried out on a supernode transputer array. Four system sizes from 125 to 512 molecules have been considered, in order to study the effect of system size on the calculated structural, orientational and dynamic properties. The dielectric constant and the dielectric relaxation time are compared with experimental data.     

Modification of a Computer Simulation Model for a Plant-Nematode System

New data on egg development and death rates, and refinements of logic concerning interaction of the nematode and host, were incorporated into a simulation model of a Meloidogyne arenaria and grapevine system. Simulations of field data improved but other areas of weakness in the model were discovered. Two peaks in the egg population curve suggested that the nematode was able to complete two life cycles before host dormancy and declining temperatures limited physiological activity.     

Computer Simulation of Isothermal Mass Transport in Graphite Slit Pores

Abstract

Results are presented from a simulation study of the mass transport of oxygen and nitrogen through graphite slit pores. The work is motivated by an attempt to understand the molecular origins of the kinetic selectivity displayed when air is separated into its major components using pressure swing adsorption. A combination of non-equilibrium molecular dynamics (NEMD), equilibrium molecular dynamics (EMD) and grand canonical Monte Carlo methods has been employed in our study to extract the maximum information. Transport diffusivities, self-diffusivities, permeabilities and Darken thermodynamic factors have been calculated as a function of pore width and temperature for pure component oxygen and nitrogen. In addition, new EMD simulation data for an 80:20 mixture of nitrogen and oxygen is reported, including a direct calculation of the Stefan-Maxwell coefficients. The results are discussed in terms of the oxygen selectivity and the possible mechanisms, which increase or decrease this quantity.

We find that the pore width behaviour of the diffusion coefficients consists of three distinct regimes: a regime at larger pore widths in which single component diffusion coefficients are largely independent of pore width, an optimum pore width at which both diffusivities increase substantially but the slit pore is selective towards nitrogen, and a regime at very low pore widths at which the diffusivities decrease sharply, but the slits are selective towards oxygen. The mechanism behind each of these regimes is discussed in terms of “entropic” effects and potential barrier heights.

We have also found that permeability selectivity is substantially reduced in a mixture of the two gases with a composition similar to that of air. Cross diffusion coefficients in the mixture have been calculated and shown to be non-negligible.     

Parameterization and Application of an Implicit Solvent Model for Macromolecules

Abstract

As the field of theoretical biophysics begins to recognize systems of longer timescales and larger magnitude, rapid approaches for investigating these systems are required. One promising simplification of the typical system of a solute surrounded by water is the use of implicit solvation models. The generalized Born implicit solvent offers a rapid approach for computing the electrostatic effects of bulk solvent without the explicit representation of water molecules. This report describes the parameterization of a generalized Born (GB) model for protein and nucleic acid structures. As a demonstration of the usefulness of this approach, the GB model is applied toward the discrimination of misfolded and properly folded protein structures. This study attempts to illustrate the potential of the GB model for molecular dynamics simulations over longer timescales as well as for screening large structural databases.     

Molecular Dynamics Simulation on the Connection Machine

An algorithm for the efficient calculation of macromolecular force fields on the Connection Machine is described. The full force field is separated into bond interactions and non-bonding interactions. Only the latter are implemented on the Connection Machine, the former, less computationally intensive tasks are performed by an existing, conventional molecular dynamics code on the front end. Parallelization of the evaluation of non-bonding interactions is achieved by the Replicated Systolic Loop algorithm introduced in this paper. The algorithm is a variant of the Systolic Loop scheme often used for the computation of 2-particle forces for the classical N-particle problem.     

灰树花深层培养的生长动力学与计算机模拟

研究了灰树花的生长特性,建立了灰树花的生长速率模型和基质消耗动力学模型;并运用SIMULINK仿真环境,建立了灰树花在培养过程中的主要因素通气量和pH对菌丝生物影响的模型,结果表明,模拟结果较理想。     

Analysis of the Disk-to-Disk Energy Transfer Processes in C-Phycocyanin Complexes by Computer Simulation Technique

Based on the crystal structure and spectral properties of C-phycocyanin (C-PC) from cyanobacteria, models for complexes with 2 and 3 C-PC hexamer disks were built and the energy transfer dynamic properties were studied by the use of stochastic computer simulation approach. In addition, an experimental parameter of 0.056 ps^–1, corresponding to a time constant of 18 ps, derived from the previous time-resolved measurement, was used for simulation of the energy transfer process from the three terminal symmetrically equivalent ₈₄ chromophores of the core-linked disk to an ₈₄ chromophore of the allophycocyanin (APC) core. The simulation showed: (1) The disk-to-disk energy transfer can be as fast as several picoseconds. (2) The energy transfer efficiencies from the first disk to the core would depend on the length of the rod (i.e. the number of disks). Efficiencies of 0.95, 0.87, and 0.75 were found for the rods with 1, 2 and 3 hexamer disks, respectively. (3) The energy transfer along a rod in a native phycobilisome (PBS) is probably very close to the one-way manner. It is the core of PBS that makes the excitation energy be transferred fast in a nearly one-way manner.     

一种水蛭素类融合蛋白与凝血酶作用的动力学模拟

通过将凝血酶抑制剂水蛭素的C端20肽片段嫁接到血小板结合蛋白AnnexinⅤ上,可以期望获得既具有抗凝血活性,同时又具有导向性的新型工程蛋白质分子.利用计算机辅助分子设计手段模拟了该融合蛋白的分子结构,并对该融合蛋白对凝血酶的抑制作用进行了分子动力学模拟,得到了支持上述想法的结果.     

太赫兹(THz)光谱在生物大分子研究中的应用

太赫兹(THz)辐射是一种新型的远红外相干辐射源,近年来,在生物大分子研究中得到了广泛的应用,特别是在生物分子的结构和动力学特性等方面有着巨大的应用潜力.结合THz光谱的特点,介绍了利用THz光谱对蛋白质、糖类及DNA等生物大分子的探索研究,以及THz技术在测定水环境与生物分子相互作用等方面的应用.探讨了该技术在生物学领域应用中有待解决的问题及发展前景.     

Molecular Dynamics Simulation of Protein Folding with Supersecondary Structure Constraints

We integrate molecular dynamics simulation methods with a newly developed supersecondary structure prediction method and compute the structure of a protein molecule, crambin. The computed structure is similar to the crystal structure with an rms error of 3.94 Å.