Role of tyrosine hot‐spot residues at the interface of colicin E9 and immunity protein 9: A comparative free energy simulation study

The endonuclease activity of the bacterial colicin 9 enzyme is controlled by the specific and high‐affinity binding of immunity protein 9 (Im9). Molecular dynamics simulation studies in explicit solvent were used to investigate the free energy change associated with the mutation of two hot‐spot interface residues [tyrosine (Tyr): Tyr54 and Tyr55] of Im9 to Ala. In addition, the effect of several other mutations (Leu33Ala, Leu52Ala, Val34Ala, Val37Ala, Ser48Ala, and Ile53Ala) with smaller influence on binding affinity was also studied. Good qualitative agreement of calculated free energy changes and experimental data on binding affinity of the mutations was observed. The simulation studies can help to elucidate the molecular details on how the mutations influence protein–protein binding affinity. The role of solvent and conformational flexibility of the partner proteins was studied by comparing the results in the presence or absence of solvent and with or without positional restraints. Restriction of the conformational mobility of protein partners resulted in significant changes of the calculated free energies but of similar magnitude for isolated Im9 and for the complex and therefore in only modest changes of binding free energy differences. Although the overall binding free energy change was similar for the two Tyr–Ala mutations, the physical origin appeared to be different with solvation changes contributing significantly to the Tyr55Ala mutation and to a loss of direct protein–protein interactions dominating the free energy change due to the Tyr54Ala mutation. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Molecular dynamics simulation of PNPLA3 I148M polymorphism reveals reduced substrate access to the catalytic cavity

A missense mutation I148M in PNPLA3 (patatin‐like phospholipase domain‐containing 3 protein) is significantly correlated with nonalcoholic fatty liver disease (NAFLD). To glean insights into mutation's effect on enzymatic activity, we performed molecular dynamics simulation and flexible docking studies. Our data show that the size of the substrate‐access entry site is significantly reduced in mutants, which limits the access of palmitic acid to the catalytic dyad. Besides, the binding free energy calculations suggest low affinity for substrate to mutant enzyme. The substrate‐bound system simulations reveal that the spatial arrangement of palmitic acid is distinct in wild‐type from that in mutant. The substrate recognition specificity is lost due to the loop where the I148M mutation was located. Our results provide strong evidence for the mechanism by which I148M affects the enzyme activity and suggest that mediating the dynamics may offer a potential avenue for NAFLD. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Vapour-Liquid Phase Equilibria of n-alkanes by Direct Monte Carlo Simulations

Abstract

We report results of direct Monte Carlo simulations of n-pentane and n-decane at the liquidvapour interface for a number of temperatures. The intermolecular interactions are modeled using the last version of the anisotropic united atom model (AUA4). We have used the local long range correction energy and an algorithm allowing to select randomly with equal probability two different displacements. The liquid and vapour densities are in excellent agreement with experimental data and with those previously calculated using the GEMC method.     

Electroosmotic Flow in Nanoscale Parallel-plate Channels: Molecular Simulation Study and Comparison with Classical Poisson–Boltzmann Theory

Molecular dynamics simulations have been carried out for simple electrolyte systems to study the electrokinetically driven osmotic flow in parallel-plate channels of widths ～10–120?nm. The results are compared with the classical theory predictions based on the solution to the Poisson–Boltzmann equation. We find that despite some of the limitations in the Poisson–Boltzmann equation, such as assumption of the Boltzmann distribution for the ions, the classical theory captures the general trend of the variations of the osmotic flow with channel width, as characterized by the mobility of the fluid in channels between ～10 and 120?nm at moderate to low ion concentration. At moderate concentration (corresponding to relatively low surface potential), the classical theory is almost quantitative. The theory and simulation show more disagreement at low concentration, primarily caused by the high surface potential where the assumption of Boltzmann distribution becomes inaccurate. We discuss the limitations of the Poisson–Boltzmann equation as applied to the nanoscale channels.     

Lipid Membrane Structure and Dynamics Studied by All-Atom Molecular Dynamics Simulations of Hydrated Phospholipid Bilayers

Abstract

The structure and dynamics of phosphatidylcholine bilayers are examined by reviewing the results of several nanoseconds of molecular dynamics simulations on a number of bilayer and monolayer models. The lengths of these simulations, the longest single one of which was 2 nanoseconds, were sufficiently long to effectively sample many of the longer-scale motions governing the behaviour of biomembranes. These simulations reproduce many experimental observables well and provide a degree of resolution currently unavailable experimentally.     

Controllers for imposing continuum-to-molecular boundary conditions in arbitrary fluid flow geometries

We present a new parallelised controller for steering an arbitrary geometric region of a molecular dynamics (MD) simulation towards a desired thermodynamic and hydrodynamic state. We show that the controllers may be applied anywhere in the domain to set accurately an initial MD state, or solely at boundary regions to prescribe non-periodic boundary conditions (PBCs) in MD simulations. The mean molecular structure and velocity autocorrelation function remain unchanged (when sampled a few molecular diameters away from the constrained region) when compared with those distributions measured using PBCs. To demonstrate the capability of our new controllers, we apply them as non-PBCs in parallel to a complex MD mixing nano-channel and in a hybrid MD continuum simulation with a complex coupling region. The controller methodology is easily extendable to polyatomic MD fluids.     

Molecular dynamics simulations of submonolayer hexane and pentane films on graphite

We present results of molecular dynamics computer simulations of hexane (C₆H₁₄ or C6) and pentane (C₅H₁₂ or C5) adlayers physisorbed onto a graphite substrate, for various submonolayer coverages. The hexane and pentane molecules incorporate explicit hydrogens and the graphite is modelled as a six-layer all-atom structure. Even though C6 and C5 have different structures at monolayer completion, both systems generally behave similarly in the submonolayer regime and results are in reasonable agreement with experiment for both systems. Specifically, there are four distinct topological regimes involving empty space: at densities closest to full coverage, there are large domains with individual vacancies, then with decreasing density, large vacancy patches appear first, followed by the formation of connected networks of smaller domains with multiple orientations that ultimately separate into individual patches. The energetics and melting behaviour of all systems are readily understood within the framework of the topology presented at various densities.     

Study of the stability of aluminium trimeric clusters in aqueous solutions

Ab initio molecular dynamics (AIMD) based on density functional theory has been used to study small aluminium–oxygen complexes in water. Such Al–O clusters have been seen in several recent mass spectrometry studies. In this study, we have focused on trimeric Al–O clusters. The initial very compact trimeric Al–O structures opened up and formed linear Al–O chains. The typical Al–O coordination number in these chain structures was 5. We have performed long (up to 200 ps) AIMD simulations and these chain structures are stable on the nanosecond time scale. We have also studied the reactivity of the Al–O dimer and solvated Al. We found a formation path for the trimeric cluster, which has a action barrier (0.04 eV) and a reaction free energy of ? 0.55 eV. This suggests that the association of a dimer and a monomer Al–O species is fast and thermodynamically a very favourable process.     

The Glassy Wall Boundary for Simulations of Inhomogeneous Systems

Abstract

For simulations on systems that invlove a physical inhomogeneity at infinite dilution, the usual periodic boundary conditions are inappropriate. In previous studies on such systems, surfaces have been re-introduced in order to contain the particles of the sample. When smooth surfaces are used to contain liquids, undesirable structural artifacts in the liquid are promoted by these surfaces. In this study a rigid but structurally liquid-like containing surface, that we call a glassy wall boundary, is introduced, in an attempt to solve this problem. It serves as a containing surface, but “looks” nearly liquid-like to molecules in the mobile liquid near it. We illustrate the properties of this boundary for a system that consists of an isolated polyion surrounded by SPC water at 300K and ～ 1 gcm^?3. We show that this boundary reduces or eliminates some of the problems caused by a smooth surface.