Analysis of the effect of electrostatic energy truncation in molecular dynamics simulations of immunoglobulin G light chain dimer

Molecular dynamics (MD) simulations of immunoglobulin G (IgG) light chain dimer using particle mesh Ewald (PME) and cutoff methods of treating electrostatic interactions were performed. The results indicate that structural parameters (RMSD, radius of gyration, solvent accessible surface) are very similar for both schemes; however, PME simulation shows increased mobility of side chains. This leads to larger fluctuations in the distance between the monomers in the dimer molecule, and, as a consequence, results in decreased number of interactions across the dimer interface. The wall clock time of the simulations was also compared. It was shown that the PME method is approximately 30% faster than the cutoff method for the system studied on a single processor.Figure Backbone order parameters for PME (red) and cutoff (green) calculations. Thick, horizontal lines show stable secondary structures     

Monte Carlo Simulations of the Chemical Potential and Free Energy for Trimer and Hexamer Rings

Abstract

The chemical potential of a trimer and hexamer model ring system was determined by computer simulation over a range of temperatures and densities. Such ring molecules are important as model aromatic and naphthenic hydrocarbons. Thermodynamic integration of the pressure along a reversible path, Widom's ghost particle insertion method and Kirkwood's charging parameter method were used over a molecular density range of 0.05 to 0.30. Data were obtained by Monte Carlo simulation of a 96 molecule system that was modelled with a Lennard-Jones 6-12 truncated potential. The original insertion method, which does not take into account the orientation of the molecule when it is inserted, gives results for the chemical potential which deviate from that obtained using the thermodynamic pressure integration. At high density or temperature the deviation is significant. We have modified the Widom insertion technique to account for this short range orientation and find good agreement between this technique and the thermodynamic integration method for the chemical potential. We also calculated the free energy difference between our model ring molecules and ring molecules made up of hard spheres.     

Monte Carlo simulations in the isothermal-isobaric ensemble: use of a ‘shell’ particle for simulating polyatomic fluids

The recent reformulation of the isothermal-isobaric ensemble requires the use of a ‘shell’ particle to define uniquely the volume of the system, thereby avoiding the redundant counting of configurations. A previous modification of the Monte Carlo method, in which trial moves are generated and accepted consistent with the correct constant pressure partition function, is extended here to the case of polyatomic fluids. With a ‘shell’ molecule, either the centre of mass of the molecule or the location of any one of the atoms within the molecule can be chosen to define the system volume. Ensemble averages obtained with the use of the shell molecule differ from ensemble averages determined with the old (i.e. no shell particle) Monte Carlo algorithm, specifically for small system sizes, although both sets of averages become equal, as they must, in the thermodynamic limit. Monte Carlo simulations in the constant pressure ensemble for various Lennard-Jones polyatomic fluids, both for pure component and binary mixtures, demonstrate these differences for small systems. For mixtures, Monte Carlo simulations may include attempted identity swaps for the shell molecule, as the choice of which component serves as the shell molecule is arbitrary when periodic boundary conditions are applied.     

Multiple time step diffusive Langevin dynamics for proteins

We present an algorithm for simulating the long time scale dynamics of proteins and other macromolecules. Our method applies the concept of multiple time step integration to the diffusive Langevin equation, in which short time scale dynamics are replaced by friction and noise. The macromolecular force field is represented at atomic resolution. Slow motions are modeled by constrained Langevin dynamics with very large time steps, while faster degrees of freedom are kept in local thermal equilibrium. In the limit of a sufficiently large molecule, our algorithm is shown to reduce the CPU time required by two orders of magnitude. We test the algorithm on two systems, alanine dipeptide and bovine pancreatic trypsin inhibitor (BPTI), and find that it accurately calculates a variety of equilibrium and dynamical properties. In the case of BPTI, the CPU time required is reduced by nearly a factor of 60 compared to a conventional, unconstrained Langevin simulation using the same force field. Proteins 30:215–227, 1998. © 1998 Wiley-Liss, Inc.     

基于EM算法的遗传重组率估计方法

EM算法是在不完全信息资料下实现参数极大似然估计的一种通用方法．本文导出了双位点不同标记类型,包括共显性-共显性,共显性-显性和显性-显性3种模式下,估计遗传重组率的EM算法,以及获得重组率抽样方差的Bootstrap方法;并将之推广到部分个体缺失标记基因型(未检测到电泳谱带)下的重组率估计．通过大量Monte Carlo模拟研究发现: (1)连锁紧密时,样本容量对重组率的估计影响不大;连锁松散时,需要较大样本容量才可检测到连锁以及实现重组率的较精确估计．(2)用包含缺失标记的所有个体估计重组率比仅用其中的非缺失标记个体估计更准确,且可显著提高连锁检测的统计功效．     

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.     

Sex and recombination in the Hötzel aging model

Why sex evolved and it prevails in nature remain one of thegreat puzzles of evolution. Most biologists would explain that it promotes genetic variability, however this explanation suffers from several difficulties. What advantages might sex confer? The present communication aims at certain investigations related to this question, in this way we introduce sexual recombination on the Hötzel model (with males and females) and wecompare these results with those from asexual reproduction without recombination.     

Predator group size distributions in predator–prey systems

The grouping behavior is common in nature, e.g., fish school, bird flocks and insects swarms. Indeed, numerous theoretical and empirical predator-prey models have demonstrated the impact of group-living animals on ecosystems. To examine the interactions between individuals in the same group or competition between groups, we introduced different models based on Monte Carlo simulation and mean-field theory and found that the predator group sizes follow the geometric distribution and logarithmic distribution, as in previous empirical and theoretical cases. Our models also provide an intuitive explanation for these distributions. A new distribution based on the Holling-III functional response is presented; this distribution is heavy tailed in some specific cases.     

MolMod – an open access database of force fields for molecular simulations of fluids

The MolMod database is presented, which is openly accessible at http://molmod.boltzmann-zuse.de and contains intermolecular force fields for over 150 pure fluids at present. It was developed and is maintained by the Boltzmann-Zuse Society for Computational Molecular Engineering (BZS). The set of molecular models in the MolMod database provides a coherent framework for molecular simulations of fluids. The molecular models in the MolMod database consist of Lennard-Jones interaction sites, point charges, and point dipoles and quadrupoles, which can be equivalently represented by multiple point charges. The force fields can be exported as input files for the simulation programmes ms2 and ls1 mardyn, GROMACS, and LAMMPS. To characterise the semantics associated with the numerical database content, a force field nomenclature is introduced that can also be used in other contexts in materials modelling at the atomistic and mesoscopic levels. The models of the pure substances that are included in the database were generally optimised such as to yield good representations of experimental data of the vapour–liquid equilibrium with a focus on the vapour pressure and the saturated liquid density. In many cases, the models also yield good predictions of caloric, transport, and interfacial properties of the pure fluids. For all models, references to the original works in which they were developed are provided. The models can be used straightforwardly for predictions of properties of fluid mixtures using established combination rules. Input errors are a major source of errors in simulations. The MolMod database contributes to reducing such errors.     

Assessment of effects of experimental imprecision on optimized biochemical systems.

Metabolic pathways may be optimized with S-system models that prescribe profiles of control variables leading to optimal output while keeping metabolites and enzyme activities within predefined ranges. Monte Carlo simulations show how much the yield and the corresponding metabolite concentrations would be affected by inaccuracies in the experimental implementation of the prescribed profiles. For a recent model of citric acid production in Aspergillus niger, the yield is roughly normally distributed, whereas the distributions of metabolite concentrations differ greatly in shape and statistical characteristics. Even moderate inaccuracies may lead to constraint violations, which appear to be correlated with high logarithmic gains.     

Test of the Inverse Monte Carlo Method for the Calculation of Interatomic Potential Energies in Atomic Liquids

Abstract

The principle purpose of this paper is to demonstrate the use of the Inverse Monte Carlo technique for calculating pair interaction energies in monoatomic liquids from a given equilibrium property. This method is based on the mathematical relation between transition probability and pair potential given by the fundamental equation of the “importance sampling” Monte Carlo method. In order to have well defined conditions for the test of the Inverse Monte Carlo method a Metropolis Monte Carlo simulation of a Lennard Jones liquid is carried out to give the equilibrium pair correlation function determined by the assumed potential. Because an equilibrium configuration is prerequisite for an Inverse Monte Carlo simulation a model system is generated reproducing the pair correlation function, which has been calculated by the Metropolis Monte Carlo simulation and therefore representing the system in thermal equilibrium. This configuration is used to simulate virtual atom displacements. The resulting changes in atom distribution for each single simulation step are inserted in a set of non-linear equations defining the transition probability for the virtual change of configuration. The solution of the set of equations for pair interaction energies yields the Lennard Jones potential by which the equilibrium configuration has been determined.     

A Comparative Molecular Dynamics Study of Thermophilic and Mesophilic Ribonuclease HI Enzymes

Abstract

We studied a pair of homologous thermophilic and mesophilic ribonuclease HI enzymes by molecular dynamics simulations. Each protein was subjected to three 5 ns simulations in explicit water at both 310 K and 340 K. The thermophilic enzyme showed larger overall positional fluctuations at both temperatures, while only the mesophilic enzyme at the higher temperature showed significant instability. When the temperature is changed, the relative flexibility of different local segments on the two proteins changed differently. Principal component analysis showed that the simulations of the two proteins explored largely overlapping regions in the conformational space. However, at 340 K, the collective structure variations of the thermophilic protein are different from those of the mesophilic protein. Our results, although not in accordance with the view that hyperthermostability of proteins may originate from their conformational rigidity, are consistent with several recent experimental and simulation studies which showed that thermophilic proteins may be conformationally more flexible than their mesophilic counterparts. The decorrelation between conformational rigidity and hyperthermostability may be attributed to the temperature dependence and long range nature of electrostatic interactions that play more important roles in the structural stability of thermophilic proteins.     

“Nanostructures in Thin-film Epitaxy: Exploring and Exploiting Substrate-mediated Interactions”

We review recent advances made in understanding the ramifications of substrate-mediated interactions for thin-film growth. Experimental studies and first-principles calculations with density-functional theory (DFT) indicate that substrate-mediated interactions can significantly influence thin-film growth. We review the findings from our kinetic Monte Carlo simulations used to model the growth of thin films, both with and without substrate-mediated interactions. For Ag heteroepitaxy on Pt(1?1?1), the pair interaction energies and adsorbate diffusion barriers were obtained from DFT calculations. Island densities for this system show significant deviations from what is predicted by classical nucleation theory. The electronic interactions created by the adsorbed atoms lead to the formation of repulsive barriers surrounding small islands and, as a result, sharp island-size distributions are produced. The island-size distributions can be manipulated by changing the growth conditions to yield desirable island sizes and shapes.     

A Study of Conformational Equilibria in Chain Molecules. I. Liquid n-Butane

Abstract

We have performed molecular dynamics simulations for liquid n-butane in order to understand liquid structures in terms of both inter- and intra-molecular interactions. Each n-butane molecule consists of four sites interacting with LJ potential and only a dihedral angle is taken into account as the internal degree of freedom. The population of gauche conformations with respect to the ideal gas state is found to increase in the liquid state. To investigate how the intermolecular interaction affects the dihedral angle distribution, we also adopt the repulsive LJ potential (RLJ) model. It is found that the nearest neighbor packing of the methyl and/or methylene groups can be approximately represented by using only the repulsive interaction. From the dihedral angle distribution, however, the rate of the shift of RLJ model to gauche is larger than that of LJ model and the attractive force also plays a significant role in the conformational equilibrium.     

Integrity of H1 helix in prion protein revealed by molecular dynamic simulations to be especially vulnerable to changes in the relative orientation of H1 and its S1 flank

In the template-assistance model, normal prion protein (PrPC), the pathogenic cause of prion diseases such as Creutzfeldt-Jakob in human, bovine spongiform encephalopathy in cow, and scrapie in sheep, converts to infectious prion (PrPSc) through an autocatalytic process triggered by a transient interaction between PrPC and PrPSc. Conventional studies suggest the S1-H1-S2 region in PrPC to be the template of S1-S2 β-sheet in PrPSc, and the conformational conversion of PrPC into PrPSc may involve an unfolding of H1 in PrPC and its refolding into the β-sheet in PrPSc. Here we conduct a series of simulation experiments to test the idea of transient interaction of the template-assistance model. We find that the integrity of H1 in PrPC is vulnerable to a transient interaction that alters the native dihedral angles at residue Asn¹⁴³, which connects the S1 flank to H1, but not to interactions that alter the internal structure of the S1 flank, nor to those that alter the relative orientation between H1 and the S2 flank.