Effects of frustration, confinement, and surface interactions on the dimerization of an off-lattice β-barrel protein

How Efficient Is Replica Exchange Molecular Dynamics? An Analytic Approach

Replica exchange molecular dynamics (REMD) has become a standard technique for accelerating relaxation in biosimulations. Despite its widespread use, questions remain about its efficiency compared with conventional, constant temperature molecular dynamics (MD). An analytic approach is taken to describe the relative efficiency of REMD with respect to MD. This is applied to several simple two-state models and to several real proteins-protein L and the B domain of protein A-to predict the relative efficiency of REMD with respect to MD in actual applications. In agreement with others, we find the following: as long as there is a positive activation energy for folding, REMD is more efficient than MD; the effectiveness of REMD is strongly dependent on the activation enthalpy; and the efficiency of REMD for actual proteins is a strong function of the maximum temperature. Choosing the maximum temperature too high can result in REMD becoming significantly less efficient than conventional MD. A good rule of thumb appears to be to choose the maximum temperature of the REMD simulation slightly above the temperature at which the enthalpy for folding vanishes. Additionally, we find that the number of replicas in REMD, while important for simulations shorter than one or two relaxation times, has a minimal effect on the asymptotic efficiency of the method.

Reaching new levels of realism in modeling biological macromolecules in cellular environments

An increasing number of studies are aimed at modeling cellular environments in a comprehensive and realistic fashion. A major challenge in these efforts is how to bridge spatial and temporal scales over many orders of magnitude. Furthermore, there are additional challenges in integrating different aspects ranging from questions about biomolecular stability in crowded environments to the description of reactive processes on cellular scales. In this review, recent studies with models of biomolecules in cellular environments at different levels of detail are discussed in terms of their strengths and weaknesses. In particular, atomistic models, implicit representations of cellular environments, coarse-grained and spheroidal models of biomolecules, as well as the inclusion of reactive processes via reaction-diffusion models are described. Furthermore, strategies for integrating the different models into a comprehensive description of cellular environments are discussed.

Insights from coarse-grained Gō models for protein folding and dynamics

Exploring the landscape of large scale conformational changes such as protein folding at atomistic detail poses a considerable computational challenge. Coarse-grained representations of the peptide chain have therefore been developed and over the last decade have proved extremely valuable. These include topology-based Gō models, which constitute a smooth and funnel-like approximation to the folding landscape. We review the many variations of the Gō model that have been employed to yield insight into folding mechanisms. Their success has been interpreted as a consequence of the dominant role of the native topology in folding. The role of local contact density in determining protein dynamics is also discussed and is used to explain the ability of Gō-like models to capture sequence effects in folding and elucidate conformational transitions.

Self-assembly of β-sheet forming peptides into chiral fibrillar aggregates

The authors introduce a novel mid-resolution off-lattice coarse-grained model to investigate the self-assembly of beta-sheet forming peptides. The model retains most of the peptide backbone degrees of freedom as well as one interaction center describing the side chains. The peptide consists of a core of alternating hydrophobic and hydrophilic residues, capped by two oppositely charged residues. Nonbonded interactions are described by Lennard-Jones and Coulombic terms. The influence of different levels of "hydrophobic" and "steric" forces between the side chains of the peptides on the thermodynamics and kinetics of aggregation was investigated using Langevin dynamics. The model is simple enough to allow the simulation of systems consisting of hundreds of peptides, while remaining realistic enough to successfully lead to the formation of chiral, ordered beta tapes, ribbons, as well as higher order fibrillar aggregates.

The Physics of the Interactions Governing Folding and Association of Proteins

The review discusses the molecular origins of the forces and free energies that determine several things about proteins, and how experiment and theory reveal this information. The first subject is the stability of the folded, native structures. The second is the range of molecular mechanisms by which proteins find their way to those folded structures in laboratory environments. The third is the much more complex problem of how folding occurs in the cellular environment. This topic includes a discussion of crowding and of the roles of chaperone molecules. The review concludes with a discussion of protein aggregation and fibril formation and of misfolding and therapies associated with it.