Many-particle Brownian and Langevin Dynamics Simulations with the Brownmove package

<p>Abstract</p> <p>Background</p> <p>Brownian Dynamics (BD) is a coarse-grained implicit-solvent simulation method that is routinely used to investigate binary protein association dynamics, but due to its efficiency in handling large simulation volumes and particle numbers it is well suited to also describe many-protein scenarios as they often occur in biological cells.</p> <p>Results</p> <p>Here we introduce our "brownmove" simulation package which was designed to handle many-particle problems with varying particle numbers and allows for a very flexible definition of rigid and flexible protein and polymer models. Both a Brownian and a Langevin dynamics (LD) propagation scheme can be used and hydrodynamic interactions are treated efficiently with our recently introduced TEA-HI ansatz [Geyer, Winter, JCP 130 (2009) 114905]. With simulations of constrained polymers and flexible models of spherical proteins we demonstrate that it is crucial to include hydrodynamics when multi-bead models are used in BD or LD simulations. Only then both the translational and the rotational diffusion coefficients and the timescales of the internal dynamics can be reproduced correctly. In the third example project we show how constant density boundary conditions [Geyer et al, JCP 120 (2004) 4573] can be used to set up a non-equilibrium simulation of diffusional transport across an array of fixed obstacles. Finally, we demonstrate how the agglomeration dynamics of multiple particles with attractive patches can be analysed conveniently with the help of a dynamic interaction network.</p> <p>Conclusions</p> <p>Combining BD and LD propagation, fast hydrodynamics, a flexible protein model, and interfaces for "open" simulation settings, our freely available "brownmove" simulation package constitutes a new platform for coarse-grained many-particle simulations of biologically relevant diffusion and transport processes.</p

Diffusion in crowded biological environments: applications of Brownian dynamics

Biochemical reactions in living systems occur in complex, heterogeneous media with total concentrations of macromolecules in the range of 50 - 400 mgml. Molecular species occupy a significant fraction of the immersing medium, up to 40% of volume. Such complex and volume-occupied environments are generally termed 'crowded' and/or 'confined'. In crowded conditions non-specific interactions between macromolecules may hinder diffusion - a major process determining metabolism, transport, and signaling. Also, the crowded media can alter, both qualitatively and quantitatively, the reactions in vivo in comparison with their in vitro counterparts. This review focuses on recent developments in particle-based Brownian dynamics algorithms, their applications to model diffusive transport in crowded systems, and their abilities to reproduce and predict the behavior of macromolecules under in vivo conditions

Diffusion of hydrophobin proteins in solution and interactions with a graphite surface

<p>Abstract</p> <p>Background</p> <p>Hydrophobins are small proteins produced by filamentous fungi that have a variety of biological functions including coating of spores and surface adhesion. To accomplish these functions, they rely on unique interface-binding properties. Using atomic-detail implicit solvent rigid-body Brownian dynamics simulations, we studied the diffusion of HFBI, a class II hydrophobin from <it>Trichoderma reesei</it>, in aqueous solution in the presence and absence of a graphite surface.</p> <p>Results</p> <p>In the simulations, HFBI exists in solution as a mixture of monomers in equilibrium with different types of oligomers. The oligomerization state depends on the conformation of HFBI. When a Highly Ordered Pyrolytic Graphite (HOPG) layer is present in the simulated system, HFBI tends to interact with the HOPG layer through a hydrophobic patch on the protein.</p> <p>Conclusions</p> <p>From the simulations of HFBI solutions, we identify a tetrameric encounter complex stabilized by non-polar interactions between the aliphatic residues in the hydrophobic patch on HFBI. After the formation of the encounter complex, a local structural rearrangement at the protein interfaces is required to obtain the tetrameric arrangement seen in HFBI crystals. Simulations performed with the graphite surface show that, due to a combination of a geometric hindrance and the interaction of the aliphatic sidechains with the graphite layer, HFBI proteins tend to accumulate close to the hydrophobic surface.</p

Cotranslational protein assembly imposes evolutionary constraints on homomeric proteins

Cotranslational protein folding can facilitate rapid formation of functional structures. However, it might also cause premature assembly of protein complexes, if two interacting nascent chains are in close proximity. By analyzing known protein structures, we show that homomeric protein contacts are enriched towards the C-termini of polypeptide chains across diverse proteomes. We hypothesize that this is the result of evolutionary constraints for folding to occur prior to assembly. Using high-throughput imaging of protein homomers in vivo in E. coli and engineered protein constructs with N- and C-terminal oligomerization domains, we show that, indeed, proteins with C-terminal homomeric interface residues consistently assemble more efficiently than those with N-terminal interface residues. Using in vivo, in vitro and in silico experiments, we identify features that govern successful assembly of homomers, which have implications for protein design and expression optimization

Diffusion, Crowding & Protein Stability in a Dynamic Molecular Model of the Bacterial Cytoplasm

A longstanding question in molecular biology is the extent to which the behavior of macromolecules observed in vitro accurately reflects their behavior in vivo. A number of sophisticated experimental techniques now allow the behavior of individual types of macromolecule to be studied directly in vivo; none, however, allow a wide range of molecule types to be observed simultaneously. In order to tackle this issue we have adopted a computational perspective, and, having selected the model prokaryote Escherichia coli as a test system, have assembled an atomically detailed model of its cytoplasmic environment that includes 50 of the most abundant types of macromolecules at experimentally measured concentrations. Brownian dynamics (BD) simulations of the cytoplasm model have been calibrated to reproduce the translational diffusion coefficients of Green Fluorescent Protein (GFP) observed in vivo, and “snapshots” of the simulation trajectories have been used to compute the cytoplasm's effects on the thermodynamics of protein folding, association and aggregation events. The simulation model successfully describes the relative thermodynamic stabilities of proteins measured in E. coli, and shows that effects additional to the commonly cited “crowding” effect must be included in attempts to understand macromolecular behavior in vivo

Parametrization of Backbone Flexibility in a Coarse-Grained Force Field for Proteins (COFFDROP) Derived from All-Atom Explicit-Solvent Molecular Dynamics Simulations of All Possible Two-Residue Peptides

Recently, we reported the parametrization of a set of coarse-grained (CG) nonbonded potential functions, derived from all-atom explicit-solvent molecular dynamics (MD) simulations of amino acid pairs and designed for use in (implicit-solvent) Brownian dynamics (BD) simulations of proteins; this force field was named COFFDROP (COarse-grained Force Field for Dynamic Representations Of Proteins). Here, we describe the extension of COFFDROP to include bonded backbone terms derived from fitting to results of explicit-solvent MD simulations of all possible two-residue peptides containing the 20 standard amino acids, with histidine modeled in both its protonated and neutral forms. The iterative Boltzmann inversion (IBI) method was used to optimize new CG potential functions for backbone-related terms by attempting to reproduce angle, dihedral, and distance probability distributions generated by the MD simulations. In a simple test of the transferability of the extended force field, the angle, dihedral, and distance probability distributions obtained from BD simulations of 56 three-residue peptides were compared to results from corresponding explicit-solvent MD simulations. In a more challenging test of the COFFDROP force field, it was used to simulate eight intrinsically disordered proteins and was shown to quite accurately reproduce the experimental hydrodynamic radii (<i>R</i>_hydro), provided that the favorable nonbonded interactions of the force field were uniformly scaled downward in magnitude. Overall, the results indicate that the COFFDROP force field is likely to find use in modeling the conformational behavior of intrinsically disordered proteins and multidomain proteins connected by flexible linkers