Coarse-grained modeling of the actin filament derived from atomistic-scale simulations

Response of an actin filament network model under cyclic stretching through a coarse grained Monte Carlo approach

Cells are complex, dynamic systems that actively adapt to various stimuli including mechanical alterations. Central to understanding cellular response to mechanical stimulation is the organization of the cytoskeleton and its actin filament network. In this manuscript, we present a minimalistic network Monte Carlo based approach to model actin filament organization under cyclic stretching. Utilizing a coarse-grained model, a filament network is prescribed within a two-dimensional circular space through nodal connections. When cyclically stretched, the model demonstrates that a perpendicular alignment of the filaments to the direction of stretch emerges in response to nodal repositioning to minimize net nodal forces from filament stress states. In addition, the filaments in the network rearrange and redistribute themselves to reduce the overall stress by decreasing their individual stresses. In parallel, we cyclically stretch NIH 3T3 fibroblasts and find a similar cytoskeletal response. With this work, we test the hypothesis that a first-principles mechanical model of filament assembly in a confined space is by itself capable of yielding the remodeling behavior observed experimentally. Identifying minimal mechanisms sufficient to reproduce mechanical influences on cellular structure has important implications in a diversity of fields, including biology, physics, medicine, computer science, and engineering.

"Fluctuograms" reveal the intermittent intra-protein communication in subtilisin Carlsberg and correlate mechanical coupling with co-evolution

The mechanism of intra-protein communication and allosteric coupling is key to understanding the structure-property relationship of protein function. For subtilisin Carlsberg, the Ca²⁺-binding loop is distal to substrate-binding and active sites, yet the serine protease function depends on Ca²⁺ binding. The atomic molecular dynamics (MD) simulations of apo and Ca²⁺-bound subtilisin show similar structures and there is no direct evidence that subtilisin has alternative conformations. To model the intra-protein communication due to Ca²⁺ binding, we transform the sequential segments of an atomic MD trajectory into separate elastic network models to represent anharmonicity and nonlinearity effectively as the temporal and spatial variation of the mechanical coupling network. In analogy to the spectrogram of sound waves, this transformation is termed the “fluctuogram” of protein dynamics. We illustrate that the Ca²⁺-bound and apo states of subtilisin have different fluctuograms and that intra-protein communication proceeds intermittently both in space and in time. We found that residues with large mechanical coupling variation due to Ca²⁺ binding correlate with the reported mutation sites selected by directed evolution for improving the stability of subtilisin and its activity in a non-aqueous environment. Furthermore, we utilize the fluctuograms calculated from MD to capture the highly correlated residues in a multiple sequence alignment. We show that in addition to the magnitude, the variance of coupling strength is also an indicative property for the sequence correlation observed in a statistical coupling analysis. The results of this work illustrate that the mechanical coupling networks calculated from atomic details can be used to correlate with functionally important mutation sites and co-evolution.

A hallmark of protein molecules is their machine-like behaviors while carrying out biological functions. At the molecular level, molecular signals such as binding a metal ion at an action site can cause long-range effects and alter protein function. Such phenomena are often referred to as intra-protein communication or allosteric coupling. Elucidating the underlying mechanisms could lead to novel discovery of molecular modulators to regulate protein function in a more specific and effective manner. A long-standing puzzle is the roles of the anharmonicity and nonlinearity in protein dynamics. To incorporate these characters in modeling intra-protein communication, we devise a “fluctuogram” analysis to record the choreography of allosteric coupling in an atomic molecular dynamics simulation. We show that fluctuogram analysis can bridge the results of physics-based simulation and sequence alignment in bioinformatics by capturing the residues that exhibit high correlation in a multiple sequence alignment. We also show that the fluctuograms calculated from atomic details have the potential to be applied as a tool to select mutation sites for modulating protein function.

Origin of twist-bend coupling in actin filaments

Actin filaments are semiflexible polymers that display large-scale conformational twisting and bending motions. Modulation of filament bending and twisting dynamics has been linked to regulatory actin-binding protein function, filament assembly and fragmentation, and overall cell motility. The relationship between actin filament bending and twisting dynamics has not been evaluated. The numerical and analytical experiments presented here reveal that actin filaments have a strong intrinsic twist-bend coupling that obligates the reciprocal interconversion of bending energy and twisting stress. We developed a mesoscopic model of actin filaments that captures key documented features, including the subunit dimensions, interaction energies, helicity, and geometrical constraints coming from the double-stranded structure. The filament bending and torsional rigidities predicted by the model are comparable to experimental values, demonstrating the capacity of the model to assess the mechanical properties of actin filaments, including the coupling between twisting and bending motions. The predicted actin filament twist-bend coupling is strong, with a persistence length of 0.15-0.4 μm depending on the actin-bound nucleotide. Twist-bend coupling is an emergent property that introduces local asymmetry to actin filaments and contributes to their overall elasticity. Up to 60% of the filament subunit elastic free energy originates from twist-bend coupling, with the largest contributions resulting under relatively small deformations. A comparison of filaments with different architectures indicates that twist-bend coupling in actin filaments originates from their double protofilament and helical structure.

Characterizing protein energy landscape by self-learning multiscale simulations: application to a designed β-hairpin

Characterizing the energy landscape of proteins at atomic resolution is still a very challenging problem, since it simultaneously requires high accuracy in estimating specific interactions and high efficiency in conformational sampling. Here, for these two requirements to meet, we extended the self-learning multiscale simulation (SLMS) method developed recently and applied it to the designed β-hairpin CLN025. The SLMS integrates all-atom and coarse-grained (CG) models in an iterative way such that the conformational sampling is performed by the CG model, the AA energy is used to calibrate the energy landscape, and the CG model is improved by the calibrated energy landscape. We extended the SLMS in two aspects, use of the energy decomposition for self-learning of the CG potential and a two-bead/residue CG model. The results show that the self-learning greatly improved the CG potential, and with the derived CG potential, the β-hairpin CLN025 robustly folded to the native structure. The self-learning iteration progressively enhanced the context dependence in the CG potential and increased the energy gap between the native and the denatured states of the CG model, leading to a funnel-like energy landscape. By using the SLMS method, without prior knowledge of the native structure but with the help of the AA energy, we can obtain a tailor-made CG potential specific to the target protein. The method can be useful for de novo structure prediction as well.

Coarse-Grained Representations of Large Biomolecular Complexes from Low-Resolution Structural Data

High-resolution atomistic structures of many large biomolecular complexes have not yet been solved by experiments, such as X-ray crystallography or NMR. Often however low-resolution information is obtained by alternative techniques, such as cryo-electron microscopy or small-angle X-ray scattering. Coarse-grained (CG) models are an appropriate choice to computationally study these complexes given the limited resolution experimental data. One of the important questions therefore is how to define CG representations from these low-resolution density maps. This work provides a space-based essential dynamics coarse-graining (ED-CG) method to define a CG representation from a density map without detailed knowledge of its underlying atomistic structure and primary sequence information. This method is demonstrated on G-actin (both the atomic structure and its density map). It is then applied to the density maps of the Escherichia coli 70S ribosome and the microtubule. The results indicate that the method can define highly CG models that still preserve functionally important dynamics of large biomolecular complexes.

Multiscale methods for protein folding simulations

Inherently hierarchic nature of proteins makes multiscale computational methods especially useful in the studies of folding and other functional dynamics. With the multiscale strategies, one can achieve improved accuracy and efficiency by coupling the atomistic and the coarse grained simulations. Depending on the problems studied, very different implementation protocols can be used to realize the multiscale idea. Here, we give detailed introductions to the currently used multiscale protocols, together with some recent applications to the protein folding simulations in our group. The advantages and weakness, as well as the application scopes of these multiscale protocols are discussed. The directions for the future developments are also proposed.

Actin filament remodeling by actin depolymerization factor/cofilin

We investigate, using molecular dynamics, how the severing protein, actin depolymerization factor (ADF)/cofilin, modulates the structure, conformational dynamics, and mechanical properties of actin filaments. The actin and cofilactin filament bending stiffness and corresponding persistence lengths obtained from all-atom simulations are comparable to values obtained from analysis of thermal fluctuations in filament shape. Filament flexibility is strongly affected by the nucleotide-linked conformation of the actin subdomain 2 DNase-I binding loop and the filament radial mass density distribution. ADF/cofilin binding between subdomains 1 and 3 of a filament subunit triggers reorganization of subdomain 2 of the neighboring subunit such that the DNase-I binding loop (DB-loop) moves radially away from the filament. Repositioning of the neighboring subunit DB-loop significantly weakens subunit interactions along the long-pitch helix and lowers the filament bending rigidity. Lateral filament contacts between the hydrophobic loop and neighboring short-pitch helix monomers in native filaments are also compromised with cofilin binding. These works provide a molecular interpretation of biochemical solution studies documenting the disruption of filament subunit interactions and also reveal the molecular basis of actin filament allostery and its linkage to ADF/cofilin binding.

Multiscale network model for large protein dynamics

Protein dynamics is essential for gaining insight into biological functions of proteins. Although protein dynamics is well delineated by molecular model, the molecular model is computationally prohibited for simulating large protein structures. In this work, we provide a multiscale network model (MNM) that allows the efficient computation on low-frequency normal modes related to structural deformation of proteins as well as dynamic behavior of functional sites. Specifically, MNM consists of two regions, one of which is described as a low-resolution structure, while the other is dictated by a high-resolution structure. The high-resolution regions using all alpha carbons of the protein are mainly binding site parts, which play a critical function in molecules, while the low-resolution parts are constructed from a further coarse-grained model (not using all alpha carbons). The feasibility of MNM to observe the cooperative motion of a protein structure was validated. It was shown that the MNM enables us to understand functional motion of proteins with computational efficiency.

Hierarchical coarse-graining strategy for protein-membrane systems to access mesoscopic scales

An overall multiscale simulation strategy for large scale coarse-grain simulations of membrane protein systems is presented. The protein is modeled as a heterogeneous elastic network, while the lipids are modeled using the hybrid analytic-systematic (HAS) methodology, where in both cases atomistic level information obtained from molecular dynamics simulation is used to parameterize the model. A feature of this approach is that from the outset liposome length scales are employed in the simulation (i.e., on the order of 1/2 a million lipids plus protein). A route to develop highly coarse-grained models from molecular-scale information is proposed and results for N-BAR domain protein remodeling of a liposome are presented.

Defining coarse-grained representations of large biomolecules and biomolecular complexes from elastic network models

Coarse-grained (CG) models of large biomolecular complexes enable simulations of these systems over long timescales that are not accessible for atomistic molecular dynamics (MD) simulations. A systematic methodology, called essential dynamics coarse-graining (ED-CG), has been developed for defining coarse-grained sites in a large biomolecule. The method variationally determines the CG sites so that key dynamic domains in the protein are preserved in the CG representation. The original ED-CG method relies on a principal component analysis (PCA) of a MD trajectory. However, for many large proteins and multi-protein complexes such an analysis may not converge or even be possible. This work develops a new ED-CG scheme using an elastic network model (ENM) of the protein structure. In this procedure, the low-frequency normal modes obtained by ENM are used to define dynamic domains and to define the CG representation accordingly. The method is then applied to several proteins, such as the HIV-1 CA protein dimer, ATP-bound G-actin, and the Arp2/3 complex. Numerical results show that ED-CG with ENM (ENM-ED-CG) is much faster than ED-CG with PCA because no MD is necessary. The ENM-ED-CG models also capture functional essential dynamics of the proteins almost as well as those using full MD with PCA. Therefore, the ENM-ED-CG method may be better suited to coarse-grain a very large biomolecule or biomolecular complex that is too computationally expensive to be simulated by conventional MD, or when a high resolution atomic structure is not even available.

Structure and dynamics of the actin filament

We used all-atom molecular dynamics simulations to investigate the structure and properties of the actin filament, starting with either the recent Oda model or the older Holmes model. Simulations of monomeric and polymerized actin show that polymerization changes the nucleotide-binding cleft, bringing together the Q137 side chain and bound ATP in a way that may enhance the ATP hydrolysis rate in the filament. Simulations with different bound nucleotides and conformations of the DNase I binding loop show that the persistence length of the filament depends only on loop conformation. Computational modeling reveals how bound phalloidin stiffens actin filaments and inhibits the release of gamma-phosphate from ADP-P(i) actin.

Multiscale simulation of protein mediated membrane remodeling

Proteins interacting with membranes can result in substantial membrane deformations and curvatures. This effect is known in its broadest terms as membrane remodeling. This review article will survey current multiscale simulation methodologies that have been employed to examine protein mediated membrane remodeling.

REACH coarse-grained normal mode analysis of protein dimer interaction dynamics

The REACH (realistic extension algorithm via covariance Hessian) coarse-grained biomolecular simulation method is a self-consistent multiscale approach directly mapping atomistic molecular dynamics simulation results onto a residue-scale model. Here, REACH is applied to calculate the dynamics of protein-protein interactions. The intra- and intermolecular fluctuations and the intermolecular vibrational densities of states derived from atomistic molecular dynamics are well reproduced by the REACH normal modes. The phonon dispersion relations derived from the REACH lattice dynamics model of crystalline ribonuclease A are also in satisfactory agreement with the corresponding all-atom results. The REACH model demonstrates that increasing dimer interaction strength decreases the translational and rotational intermolecular vibrational amplitudes, while their vibrational frequencies are relatively unaffected. A comparative study of functionally interacting biological dimers with crystal dimers, which are formed artificially via crystallization, reveals a relation between their static structures and the interprotein dynamics: i.e., the consequence of the extensive interfaces of biological dimers is reduction of the intermonomer translational and rotational amplitudes, but not the frequencies.

Nucleotide-dependent conformational states of actin

The influence of the state of the bound nucleotide (ATP, ADP-Pi, or ADP) on the conformational free-energy landscape of actin is investigated. Nucleotide-dependent folding of the DNase-I binding (DB) loop in monomeric actin and the actin trimer is carried out using all-atom molecular dynamics (MD) calculations accelerated with a multiscale implementation of the metadynamics algorithm. Additionally, an investigation of the opening and closing of the actin nucleotide binding cleft is performed. Nucleotide-dependent free-energy profiles for all of these conformational changes are calculated within the framework of metadynamics. We find that in ADP-bound monomer, the folded and unfolded states of the DB loop have similar relative free-energy. This result helps explain the experimental difficulty in obtaining an ordered crystal structure for this region of monomeric actin. However, we find that in the ADP-bound actin trimer, the folded DB loop is stable and in a free-energy minimum. It is also demonstrated that the nucleotide binding cleft favors a closed conformation for the bound nucleotide in the ATP and ADP-Pi states, whereas the ADP state favors an open confirmation, both in the monomer and trimer. These results suggest a mechanism of allosteric interactions between the nucleotide binding cleft and the DB loop. This behavior is confirmed by an additional simulation that shows the folding free-energy as a function of the nucleotide cleft width, which demonstrates that the barrier for folding changes significantly depending on the value of the cleft width.

Models of protein linear molecular motors for dynamic nanodevices

Protein molecular motors are natural nano-machines that convert the chemical energy from the hydrolysis of adenosine triphosphate into mechanical work. These efficient machines are central to many biological processes, including cellular motion, muscle contraction and cell division. The remarkable energetic efficiency of the protein molecular motors coupled with their nano-scale has prompted an increasing number of studies focusing on their integration in hybrid micro- and nanodevices, in particular using linear molecular motors. The translation of these tentative devices into technologically and economically feasible ones requires an engineering, design-orientated approach based on a structured formalism, preferably mathematical. This contribution reviews the present state of the art in the modelling of protein linear molecular motors, as relevant to the future design-orientated development of hybrid dynamic nanodevices.

Bending flexibility of actin filaments during motor-induced sliding

Muscle contraction and other forms of cell motility occur as a result of cyclic interactions between myosin molecules and actin filaments. Force generation is generally attributed to ATP-driven structural changes in myosin, whereas a passive role is ascribed to actin. However, some results challenge this view, predicting structural changes in actin during motor activity, e.g., when the actin filaments slide on a myosin-coated surface in vitro. Here, we analyzed statistical properties of the sliding filament paths, allowing us to detect changes of this type. It is interesting to note that evidence for substantial structural changes that led to increased bending flexibility of the filaments was found in phalloidin-stabilized, but not in phalloidin-free, actin filaments. The results are in accordance with the idea that a high-flexibility structural state of actin is a prerequisite for force production, but not the idea that a low-to-high flexibility transition of the actin filament should be an important component of the force-generating step per se. Finally, our data challenge the general view that phalloidin-stabilized filaments behave as native actin filaments in their interaction with myosin. This has important implications, since phalloidin stabilization is a routine procedure in most studies of actomyosin function.

Molecular dynamics simulation and coarse-grained analysis of the Arp2/3 complex

A molecular dynamics investigation and coarse-grained analysis of inactivated actin-related protein (Arp) 2/3 complex is presented. It was found that the nucleotide binding site within Arp3 remained in a closed position with bound ATP or ADP, but opened when simulation with no nucleotide was performed. In contrast, simulation of the isolated Arp3 subunit with bound ATP, showed a fast opening of the nucleotide binding cleft. A homology model for the missing subdomains 1 and 2 of Arp2 was constructed, and it was also found that the Arp2 binding cleft remained closed with bound nucleotide. Within the nucleotide binding cleft a distinct opening and closing period of 10 ns was observed in many of the simulations of Arp2/3 as well as isolated Arp3. Substitution studies were employed, and several alanine substitutions were found to induce a partial opening of the ATP binding cleft in Arp3 and Arp2, whereas only a single substitution was found to induce opening of the ADP binding cleft. It was also found that the nucleotide type did not cause a substantial change on interfacial contacts between Arp3 and the ArpC2, ArpC3 and ArpC4 subunits. Nucleotide-free Arp3 had generally less stable contacts, but the overall contact architecture was constant. Finally, nucleotide-dependent coarse-grained models for Arp3 are developed that serve to further highlight the structural differences induced in Arp3 by nucleotide hydrolysis.

A systematic methodology for defining coarse-grained sites in large biomolecules

Coarse-grained (CG) models of biomolecules have recently attracted considerable interest because they enable the simulation of complex biological systems on length-scales and timescales that are inaccessible for atomistic molecular dynamics simulation. A CG model is defined by a map that transforms an atomically detailed configuration into a CG configuration. For CG models of relatively small biomolecules or in cases that the CG and all-atom models have similar resolution, the construction of this map is relatively straightforward and can be guided by chemical intuition. However, it is more challenging to construct a CG map when large and complex domains of biomolecules have to be represented by relatively few CG sites. This work introduces a new and systematic methodology called essential dynamics coarse-graining (ED-CG). This approach constructs a CG map of the primary sequence at a chosen resolution for an arbitrarily complex biomolecule. In particular, the resulting ED-CG method variationally determines the CG sites that reflect the essential dynamics characterized by principal component analysis of an atomistic molecular dynamics trajectory. Numerical calculations illustrate this approach for the HIV-1 CA protein dimer and ATP-bound G-actin. Importantly, since the CG sites are constructed from the primary sequence of the biomolecule, the resulting ED-CG model may be better suited to appropriately explore protein conformational space than those from other CG methods at the same degree of resolution.

Systematic multiscale parameterization of heterogeneous elastic network models of proteins

We present a method to parameterize heterogeneous elastic network models (heteroENMs) of proteins to reproduce the fluctuations observed in atomistic simulations. Because it is based on atomistic simulation, our method allows the development of elastic coarse-grained models of proteins under different conditions or in different environments. The method is simple and applicable to models at any level of coarse-graining. We validated the method in three systems. First, we computed the persistence length of ADP-bound F-actin, using a heteroENM model. The value of 6.1 +/- 1.6 microm is consistent with the experimentally measured value of 9.0 +/- 0.5 microm. We then compared our method to a uniform elastic network model and a realistic extension algorithm via covariance Hessian (REACH) model of carboxy myoglobin, and found that the heteroENM method more accurately predicted mean-square fluctuations of alpha-carbon atoms. Finally, we showed that the method captures critical differences in effective harmonic interactions for coarse-grained models of the N-terminal Bin/amphiphysin/Rvs (N-BAR) domain of amphiphysin, by building models of N-BAR both bound to a membrane and free in solution.

Coarse-grained force field for the nucleosome from self-consistent multiscaling

A coarse-grained simulation model for the nucleosome is developed, using a methodology modified from previous work on the ribosome. Protein residues and DNA nucleotides are represented as beads, interacting through harmonic (for neighboring) or Morse (for nonbonded) potentials. Force-field parameters were estimated by Boltzmann inversion of the corresponding radial distribution functions obtained from a 5-ns all-atom molecular dynamics (MD) simulation, and were refined to produce agreement with the all-atom MD simulation. This self-consistent multiscale approach yields a coarse-grained model that is capable of reproducing equilibrium structural properties calculated from a 50-ns all-atom MD simulation. This coarse-grained model speeds up nucleosome simulations by a factor of 10(3) and is expected to be useful in examining biologically relevant dynamical nucleosome phenomena on the microsecond timescale and beyond.

Cofilin increases the bending flexibility of actin filaments: implications for severing and cell mechanics

We determined the flexural (bending) rigidities of actin and cofilactin filaments from a cosine correlation function analysis of their thermally driven, two-dimensional fluctuations in shape. The persistence length of actin filaments is 9.8 microm, corresponding to a flexural rigidity of 0.040 pN microm(2). Cofilin binding lowers the persistence length approximately 5-fold to a value of 2.2 microm and the filament flexural rigidity to 0.0091 pN microm(2). That cofilin-decorated filaments are more flexible than native filaments despite an increased mass indicates that cofilin binding weakens and redistributes stabilizing subunit interactions of filaments. We favor a mechanism in which the increased flexibility of cofilin-decorated filaments results from the linked dissociation of filament-stabilizing ions and reorganization of actin subdomain 2 and as a consequence promotes severing due to a mechanical asymmetry. Knowledge of the effects of cofilin on actin filament bending mechanics, together with our previous analysis of torsional stiffness, provide a quantitative measure of the mechanical changes in actin filaments associated with cofilin binding, and suggest that the overall mechanical and force-producing properties of cells can be modulated by cofilin activity.

REACH coarse-grained biomolecular simulation: transferability between different protein structural classes

Coarse graining of protein interactions provides a means of simulating large biological systems. The REACH (Realistic Extension Algorithm via Covariance Hessian) coarse-graining method, in which the force constants of a residue-scale elastic network model are calculated from the variance-covariance matrix obtained from atomistic molecular dynamics (MD) simulation, involves direct mapping between scales without the need for iterative optimization. Here, the transferability of the REACH force field is examined between protein molecules of different structural classes. As test cases, myoglobin (all alpha), plastocyanin (all beta), and dihydrofolate reductase (alpha/beta) are taken. The force constants derived are found to be closely similar in all three proteins. An MD version of REACH is presented, and low-temperature coarse-grained (CG) REACH MD simulations of the three proteins are compared with atomistic MD results. The mean-square fluctuations of the atomistic MD are well reproduced by the CGMD. Model functions for the CG interactions, derived by averaging over the three proteins, are also shown to produce fluctuations in good agreement with the atomistic MD. The results indicate that, similarly to the use of atomistic force fields, it is now possible to use a single, generic REACH force field for all protein studies, without having first to derive parameters from atomistic MD simulation for each individual system studied. The REACH method is thus likely to be a reliable way of determining spatiotemporal motion of a variety of proteins without the need for expensive computation of long atomistic MD simulations.

A finite element framework for computation of protein normal modes and mechanical response

A computational framework based on the Finite Element Method is presented to calculate the normal modes and mechanical response of proteins and their supramolecular assemblies. Motivated by elastic network models, proteins are treated as continuum elastic solids with molecular volume defined by their solvent-excluded surface. The discretized Finite Element representation is obtained using a surface simplification algorithm that facilitates the generation of models of arbitrary prescribed spatial resolution. The procedure is applied to a mutant of T4 phage lysozyme, G-actin, syntenin, cytochrome-c', beta-tubulin, and the supramolecular assembly filamentous actin (F-actin). Equilibrium thermal fluctuations of alpha-carbon atoms and their inter-residue correlations compare favorably with all-atom-based results, the Rotational-Translational Block procedure, and experiment. Additionally, the free vibration and compressive buckling responses of F-actin are in quantitative agreement with experiment. The proposed methodology is applicable to any protein or protein assembly and facilitates the incorporation of specific atomic-level interactions, including aqueous-electrolyte-mediated electrostatic effects and solvent damping. The procedure is equally applicable to proteins with known atomic coordinates as it is to electron density maps of proteins, protein complexes, and supramolecular assemblies of unknown atomic structure.

Coarse-grained models reveal functional dynamics--II. Molecular dynamics simulation at the coarse-grained level--theories and biological applications

Molecular dynamics (MD) simulation has remained the most indispensable tool in studying equilibrium/non-equilibrium conformational dynamics since its advent 30 years ago. With advances in spectroscopy accompanying solved biocomplexes in growing sizes, sampling their dynamics that occur at biologically interesting spatial/temporal scales becomes computationally intractable; this motivated the use of coarse-grained (CG) approaches. CG-MD models are used to study folding and conformational transitions in reduced resolution and can employ enlarged time steps due to the absence of some of the fastest motions in the system. The Boltzmann-Inversion technique, heavily used in parameterizing these models, provides a smoothed-out effective potential on which molecular conformation evolves at a faster pace thus stretching simulations into tens of microseconds. As a result, a complete catalytic cycle of HIV-1 protease or the assembly of lipid-protein mixtures could be investigated by CG-MD to gain biological insights. In this review, we survey the theories developed in recent years, which are categorized into Folding-based and Molecular-Mechanics-based. In addition, physical bases in the selection of CG beads/time-step, the choice of effective potentials, representation of solvent, and restoration of molecular representations back to their atomic details are systematically discussed.

Unique elastic properties of the spectrin tetramer as revealed by multiscale coarse-grained modeling

The force-extension profile of tetrameric spectrin is determined by using multiscale computer simulation. Fluctuation results of atomistic simulations of double spectrin repeat units (DSRU) are used to systematically build a coarse-grained (CG) model for the tetrameric form of spectrin. It is found that the spectrin tetramer can be modeled as a soft polymer with a unique flat force-extension profile over the range of biologically important lengths. It is also concluded that in the cytoskeletal network of the red blood cell the tetramer is in an "overcompressed" state. These findings are in contrast to the commonly used models of spectrin tetramer elasticity, namely the "entropic spring" polymer models. From these results, it is concluded that stable intact helical linker regions are needed to maintain the soft elasticity of the spectrin tetramer.

Microseconds dynamics simulations of the outer-membrane protease T

Conformational fluctuations of enzymes may play an important role for substrate recognition and/or catalysis, as it has been suggested in the case of the protease enzymatic superfamily. Unfortunately, theoretically addressing this issue is a problem of formidable complexity, as the number of the involved degrees of freedom is enormous: indeed, the biological function of a protein depends, in principle, on all its atoms and on the surrounding water molecules. Here we investigated a membrane protease enzyme, the OmpT from Escherichia coli, by a hybrid molecular mechanics/coarse-grained approach, in which the active site is treated with the GROMOS force field, whereas the protein scaffold is described with a Go-model. The method has been previously tested against results obtained with all-atom simulations. Our results show that the large-scale motions and fluctuations of the electric field in the microsecond timescale may impact on the biological function and suggest that OmpT employs the same catalytic strategy as aspartic proteases. Such a conclusion cannot be drawn within the 10- to 100-ns timescale typical of current molecular dynamics simulations. In addition, our studies provide a structural explanation for the drop in the catalytic activity of two known mutants (S99A and H212A), suggesting that the coarse-grained approach is a fast and reliable tool for providing structure/function relationships for both wild-type OmpT and mutants.

Brushes, cables, and anchors: recent insights into multiscale assembly and mechanics of cellular structural networks

The remarkable ability of living cells to sense, process, and respond to mechanical stimuli in their environment depends on the rapid and efficient interconversion of mechanical and chemical energy at specific times and places within the cell. For example, application of force to cells leads to conformational changes in specific mechanosensitive molecules which then trigger cellular signaling cascades that may alter cellular structure, mechanics, and migration and profoundly influence gene expression. Similarly, the sensitivity of cells to mechanical stresses is governed by the composition, architecture, and mechanics of the cellular cytoskeleton and extracellular matrix (ECM), which are in turn driven by molecular-scale forces between the constituent biopolymers. Understanding how these mechanochemical systems coordinate over multiple length and time scales to produce orchestrated cell behaviors represents a fundamental challenge in cell biology. Here, we review recent advances in our understanding of these complex processes in three experimental systems: the assembly of axonal neurofilaments, generation of tensile forces by actomyosin stress fiber bundles, and mechanical control of adhesion assembly.

Coarse-grained free energy functions for studying protein conformational changes: a double-well network model

In this work, a double-well network model (DWNM) is presented for generating a coarse-grained free energy function that can be used to study the transition between reference conformational states of a protein molecule. Compared to earlier work that uses a single, multidimensional double-well potential to connect two conformational states, the DWNM uses a set of interconnected double-well potentials for this purpose. The DWNM free energy function has multiple intermediate states and saddle points, and is hence a "rough" free energy landscape. In this implementation of the DWNM, the free energy function is reduced to an elastic-network model representation near the two reference states. The effects of free energy function roughness on the reaction pathways of protein conformational change is demonstrated by applying the DWNM to the conformational changes of two protein systems: the coil-to-helix transition of the DB-loop in G-actin and the open-to-closed transition of adenylate kinase. In both systems, the rough free energy function of the DWNM leads to the identification of distinct minimum free energy paths connecting two conformational states. These results indicate that while the elastic-network model captures the low-frequency vibrational motions of a protein, the roughness in the free energy function introduced by the DWNM can be used to characterize the transition mechanism between protein conformations.

Coarse-grained biomolecular simulation with REACH: realistic extension algorithm via covariance Hessian

Coarse-graining of protein interactions provides a means of simulating large biological systems. Here, a coarse-graining method, REACH, is introduced, in which the force constants of a residue-scale elastic network model are calculated from the variance-covariance matrix obtained from atomistic molecular dynamics (MD) simulation. In test calculations, the C(alpha)-atoms variance-covariance matrices are calculated from the ensembles of 1-ns atomistic MD trajectories in monomeric and dimeric myoglobin, and used to derive coarse-grained force constants for the local and nonbonded interactions. Construction of analytical model functions of the distance-dependence of the interresidue force constants allows rapid calculation of the REACH normal modes. The model force constants from monomeric and dimeric myoglobin are found to be similar in magnitude to each other. The MD intra- and intermolecular mean-square fluctuations and the vibrational density of states are well reproduced by the residue-scale REACH normal modes without requiring rescaling of the force constant parameters. The temperature-dependence of the myoglobin REACH force constants reveals that the dynamical transition in protein internal fluctuations arises principally from softening of the elasticity in the nonlocal interactions. The REACH method is found to be a reliable way of determining spatiotemporal protein motion without the need for expensive computations of long atomistic MD simulations.

Bridging from molecular simulation to biochemical networks

How can we make the connection between the three-dimensional structures of individual proteins and understanding how complex biological systems involving many proteins work? The modelling and simulation of protein structures can help to answer this question for systems ranging from multimacromolecular complexes to organelles and cells. On one hand, multiscale modelling and simulation techniques are advancing to permit the spatial and temporal properties of large systems to be simulated using atomic-detail structures. On the other hand, the estimation of kinetic parameters for the mathematical modelling of biochemical pathways using protein structure information provides a basis for iterative manipulation of biochemical pathways guided by protein structure. Recent advances include the structural modelling of protein complexes on the genomic level, novel coarse-graining strategies to increase the size of the system and the time span that can be simulated, and comparative molecular field analyses to estimate enzyme kinetic parameters.

Multiscale modeling of biomolecular systems: in serial and in parallel

Considerable progress has been recently achieved in the multiscale modeling of complex biological processes. Multiscale models have now investigated the structure and dynamics of lipid membranes, proteins, peptides and DNA over length and time scales ranging from the atomic to the macroscopic. Serial multiscale methods that parameterize low-resolution coarse-grained models with data from high-resolution models have studied long time or length scale phenomena that cannot be investigated with atomically detailed models. Parallel multiscale methods that directly couple high- and low-resolution models have efficiently explored slow structural transitions and the importance of long-wavelength fluctuations for biological molecules. The success of such models relies upon new theories and methods for constructing accurate multiscale bridges that transfer information between models with different resolutions.

Atomistic and Coarse-grained Analysis of Double Spectrin Repeat Units: The Molecular Origins of Flexibility

Spectrin is an ubiquitous protein in metazoan cells, and its flexibility is one of the keys to maintaining cellular structure and organization. Both alpha-spectrin and beta-spectrin polypeptides consist primarily of triple coiled-coil modular repeat units, and two important factors that determine spectrin flexibility are the bending flexibility between two consecutive repeat units and the conformational flexibility of individual repeat units. Atomistic molecular dynamics (MD) simulations are used here to study double spectrin repeat units (DSRUs) from the human erythrocyte beta-spectrin (HEbeta89) and the chicken brain alpha-spectrin (CBalpha1617). From the results of MD simulations, a highly conserved Trp residue in the A-helix of most repeat units that has been suggested to be important in conferring stability to the coiled-coil structures is found not to have a significant effect on the conformational flexibility of individual repeat units. Characterization of the bending flexibility for two consecutive repeats of spectrin via atomistic simulations and coarse-grained (CG) modeling indicate that the bending flexibility is governed by the interactions between the AB-loop of the first repeat unit, the BC-loop of the second repeat unit and the linker region. Specifically, interactions between residues in these regions can lead to a strong directionality in the bending behavior of two repeat units. The biological implications of these finding are discussed.

Bending elasticity of anti-parallel beta-sheets

Using a coarse-grained elastic model, we examine the bending properties of anti-parallel beta-sheets comprised of uniform amino-acid residues in vacuum as well as in explicit solvent. By comparing the conformational probability of the beta-sheet from molecular dynamics simulations with the same quantities obtained from the coarse-grained model, we compute the elastic bending constant, kappa. Equilibrium fluctuations of the beta-sheet and its response to external forces are well reproduced by a model with a uniform isotropic bending constant. An anisotropic bending model is also investigated, although the computed anisotropy is relatively weak and most of the observed properties are well described by an isotropic model. The presence of explicit solvent also lowers the bending constant. The sequence dependence of our result and its implications in protein conformational dynamics are discussed.