Classification of GTP-dependent K-Ras4B active and inactive conformational states.

Classifying reliably active and inactive molecular conformations of wildtype (WT) and mutated oncogenic proteins is a key, ongoing challenge in molecular cancer studies. Here, we probe the GTP-bound K-Ras4B conformational dynamics using long-time atomistic molecular dynamics (MD) simulations. We extract and analyze the detailed underlying free energy landscape of WT K-Ras4B. We use two key reaction coordinates, labeled d1 and d2 (i.e., distances coordinating the Pß atom of the GTP ligand with two key residues, T35 and G60), shown to correlate closely with activities of WT and mutated K-Ras4B. However, our new K-Ras4B conformational kinetics study reveals a more complex network of equilibrium Markovian states. We show that a new reaction coordinate is required to account for the orientation of acidic K-Ras4B sidechains such as D38 with respect to the interface with binding effector RAF1 and rationalize the activation/inactivation propensities and the corresponding molecular binding mechanisms. We use this understanding to unveil how a relatively conservative mutation (i.e., D33E, in the switch I region) can lead to significantly different activation propensities compared with WT K-Ras4B. Our study sheds new light on the ability of residues near the K-Ras4B-RAF1 interface to modulate the network of salt bridges at the binding interface with the RAF1 downstream effector and, thus, to influence the underlying GTP-dependent activation/inactivation mechanism. Altogether, our hybrid MD-docking modeling approach enables the development of new in silico methods for quantitative assessment of activation propensity changes (e.g., due to mutations or local binding environment). It also unveils the underlying molecular mechanisms and facilitates the rational design of new cancer drugs.

MST2-RASSF protein-protein interactions through SARAH domains.

The detailed, atomistic-level understanding of molecular signaling along the tumor-suppressive Hippo signaling pathway that controls tissue homeostasis by balancing cell proliferation and death through apoptosis is a promising avenue for the discovery of novel anticancer drug targets. The activation of kinases such as Mammalian STE20-Like Protein Kinases 1 and 2 (MST1 and MST2)-modulated through both homo- and heterodimerization (e.g. interactions with Ras association domain family, RASSF, enzymes)-is a key upstream event in this pathway and remains poorly understood. On the other hand, RASSFs (such as RASSF1A or RASSF5) act as important apoptosis activators and tumor suppressors, although their exact regulatory roles are also unclear. We present recent molecular studies of signaling along the Ras-RASSF-MST pathway, which controls growth and apoptosis in eukaryotic cells, including a variety of modern molecular modeling and simulation techniques. Using recently available structural information, we discuss the complex regulatory scenario according to which RASSFs perform dual signaling functions, either preventing or promoting MST2 activation, and thus control cell apoptosis. Here, we focus on recent studies highlighting the special role being played by the specific interactions between the helical Salvador/RASSF/Hippo (SARAH) domains of MST2 and RASSF1a or RASSF5 enzymes. These studies are crucial for integrating atomistic-level mechanistic information about the structures and conformational dynamics of interacting proteins, with information available on their system-level functions in cellular signaling.

Conformational analysis of replica exchange MD: Temperature-dependent Markov networks for FF amyloid peptides.

Recent molecular modeling methods using Markovian descriptions of conformational states of biomolecular systems have led to powerful analysis frameworks that can accurately describe their complex dynamical behavior. In conjunction with enhanced sampling methods, such as replica exchange molecular dynamics (REMD), these frameworks allow the systematic and accurate extraction of transition probabilities between the corresponding states, in the case of Markov state models, and of statistically-optimized transition rates, in the case of the corresponding coarse master equations. However, applying automatically such methods to large molecular dynamics (MD) simulations, with explicit water molecules, remains limited both by the initial ability to identify good candidates for the underlying Markovian states and by the necessity to do so using good collective variables as reaction coordinates that allow the correct counting of inter-state transitions at various lag times. Here, we show that, in cases when representative molecular conformations can be identified for the corresponding Markovian states, and thus their corresponding collective evolution of atomic positions can be calculated along MD trajectories, one can use them to build a new type of simple collective variable, which can be particularly useful in both the correct state assignment and in the subsequent accurate counting of inter-state transition probabilities. In the case of the ubiquitously used root-mean-square deviation (RMSD) of atomic positions, we introduce the relative RMSD (RelRMSD) measure as a good reaction coordinate candidate. We apply this method to the analysis of REMD trajectories of amyloid-forming diphenylalanine (FF) peptides-a system with important nanotechnology and biomedical applications due to its self-assembling and piezoelectric properties-illustrating the use of RelRMSD in extracting its temperature-dependent intrinsic kinetics, without a priori assumptions on the functional form (e.g., Arrhenius or not) of the underlying conformational transition rates. The RelRMSD analysis enables as well a more objective assessment of the convergence of the REMD simulations. This type of collective variable may be generalized to other observables that could accurately capture conformational differences between the underlying Markov states (e.g., distance RMSD, the fraction of native contacts, etc.).

SARAH Domain-Mediated MST2-RASSF Dimeric Interactions.

RASSF enzymes act as key apoptosis activators and tumor suppressors, being downregulated in many human cancers, although their exact regulatory roles remain unknown. A key downstream event in the RASSF pathway is the regulation of MST kinases, which are main effectors of RASSF-induced apoptosis. The regulation of MST1/2 includes both homo- and heterodimerization, mediated by helical SARAH domains, though the underlying molecular interaction mechanism is unclear. Here, we study the interactions between RASSF1A, RASSF5, and MST2 SARAH domains by using both atomistic molecular simulation techniques and experiments. We construct and study models of MST2 homodimers and MST2-RASSF SARAH heterodimers, and we identify the factors that control their high molecular stability. In addition, we also analyze both computationally and experimentally the interactions of MST2 SARAH domains with a series of synthetic peptides particularly designed to bind to it, and hope that our approach can be used to address some of the challenging problems in designing new anti-cancer drugs.

Peptide dimerization-dissociation rates from replica exchange molecular dynamics.

We show how accurate rates of formation and dissociation of peptide dimers can be calculated using direct transition counting (DTC) from replica-exchange molecular dynamics (REMD) simulations. First, continuous trajectories corresponding to system replicas evolving at different temperatures are used to assign conformational states. Second, we analyze the entire REMD data to calculate the corresponding rates at each temperature directly from the number of transition counts. Finally, we compare the kinetics extracted directly, using the DTC method, with indirect estimations based on trajectory likelihood maximization using short-time propagators and on decay rates of state autocorrelation functions. For systems with relatively low-dimensional intrinsic conformational dynamics, the DTC method is simple to implement and leads to accurate temperature-dependent rates. We apply the DTC rate-extraction method to all-atom REMD simulations of dimerization of amyloid-forming NNQQ tetrapetides in explicit water. In an assessment of the REMD sampling efficiency with respect to standard MD, we find a gain of more than a factor of two at the lowest temperature.

Thermal and aqueous stability improvement of graphene oxide enhanced diphenylalanine nanocomposites.

Nanocomposites of diphenylalanine (FF) and carbon based materials provide an opportunity to overcome drawbacks associated with using FF micro- and nanostructures in nanobiotechnology applications, in particular their poor structural stability in liquid solutions. In this study, FF/graphene oxide (GO) composites were found to self-assemble into layered micro- and nanostructures, which exhibited improved thermal and aqueous stability. Dependent on the FF/GO ratio, the solubility of these structures was reduced to 35.65% after 30 min as compared to 92.4% for pure FF samples. Such functional nanocomposites may extend the use of FF structures to e.g. biosensing, electrochemical, electromechanical or electronic applications.

Phosphorylation of RAF Kinase Dimers Drives Conformational Changes that Facilitate Transactivation.

RAF kinases are key players in the MAPK signaling pathway and are important targets for personalized cancer therapy. RAF dimerization is part of the physiological activation mechanism, together with phosphorylation, and is known to convey resistance to RAF inhibitors. Herein, molecular dynamics simulations are used to show that phosphorylation of a key N-terminal acidic (NtA) motif facilitates RAF dimerization by introducing several interprotomer salt bridges between the αC-helix and charged residues upstream of the NtA motif. Additionally, we show that the R-spine of RAF interacts with a conserved Trp residue in the vicinity of the NtA motif, connecting the active sites of two protomers and thereby modulating the cooperative interactions in the RAF dimer. Our findings provide a first structure-based mechanism for the auto-transactivation of RAF and could be generally applicable to other kinases, opening new pathways for overcoming dimerization-related drug resistance.

Membrane Protein Structure, Function, and Dynamics: a Perspective from Experiments and Theory.

Membrane proteins mediate processes that are fundamental for the flourishing of biological cells. Membrane-embedded transporters move ions and larger solutes across membranes; receptors mediate communication between the cell and its environment and membrane-embedded enzymes catalyze chemical reactions. Understanding these mechanisms of action requires knowledge of how the proteins couple to their fluid, hydrated lipid membrane environment. We present here current studies in computational and experimental membrane protein biophysics, and show how they address outstanding challenges in understanding the complex environmental effects on the structure, function, and dynamics of membrane proteins.

Exposure of chick embryos to cadmium changes the extra-embryonic vascular branching pattern and alters expression of VEGF-A and VEGF-R2.

Cadmium (Cd) has several industrial applications, and is found in tobacco products, a notable source of human exposure. Vascular endothelial cells are key targets of Cd toxicity. Here, we aim to quantify the alteration to vascular branching pattern following Cd exposure in the chick extra-embryonic membrane (EEM) using fractal analysis, and explore molecular cues to angiogenesis such as VEGF-A and VEGF-R2 expression following Cd treatment. Chicken embryos were incubated for 60 h to Hamburger-Hamilton developmental stage 16-17, then explanted and treated with 50 µL of 50 µmol cadmium acetate (CdAc) or an equivalent volume of equimolar sodium acetate (NaAc). Images of embryos and their area vasculosa (AV) were captured and analyzed at 4 different time points (4, 8, 24 and 48 h) following treatment. Vascular branching in the AV was quantified using its fractal dimension (Df), estimated using a box counting method. Gallinaceous VEGF ELISA was used to measure the VEGF-A concentration in the EEM following treatment, with determination of the relative expression of VEGF-A and VEGF-R2 using quantitative real-time RT-PCR. Vascular branching increased monotonically in the control group at all time points. The anti-angiogenic effect of Cd exposure on the AV was reflected by a significant reduction in Df when compared with controls. Df was more markedly reduced in cultures with abnormal embryos. The expression of VEGF-A protein, and VEGF-A and VEGF-R2 mRNA were reduced in Cd-exposed EEMs. Both molecules contribute to growth, vessel sprouting and branching processes, which supports our findings using fractal analysis.

Structure and dynamics of the fibronectin-III domains of Aplysia californica cell adhesion molecules.

Due to their homophilic and heterophilic binding properties, cell adhesion molecules (CAMs) such as integrin, cadherin and the immunoglobulin superfamily CAMs are of primary importance in cell-cell and cell-substrate interactions, signalling pathways and other crucial biological processes. We study the molecular structures and conformational dynamics of the two fibronectin type III (Fn-III) extracellular domains of the Aplysia californica CAM (apCAM) protein, by constructing and probing an atomically-detailed structural model based on apCAM's homology with other CAMs. The stability and dynamic properties of the Fn-III domains, individually and in tandem, are probed and analysed using all-atom explicit-solvent molecular dynamics (MD) simulations and normal mode analysis of their corresponding elastic network models. The refined structural model of the Fn-III tandem of apCAM reveals a specific pattern of amino acid interactions that controls the stability of the ß-sheet rich structure and could affect apCAM's response to physical or chemical changes of its environment. It also exposes the important role of several specific charged residues in modulating the structural properties of the linker segment connecting the two Fn-III domains, as well as of the inter-domain interface.

Molecular mechanisms of asymmetric RAF dimer activation.

Protein phosphorylation is one of the most common post-translational modifications in cell regulatory mechanisms. Dimerization plays also a crucial role in the kinase activity of many kinases, including RAF, CDK2 (cyclin-dependent kinase 2) and EGFR (epidermal growth factor receptor), with heterodimers often being the most active forms. However, the structural and mechanistic details of how phosphorylation affects the activity of homo- and hetero-dimers are largely unknown. Experimentally, synthesizing protein samples with fully specified and homogeneous phosphorylation states remains a challenge for structural biology and biochemical studies. Typically, multiple changes in phosphorylation lead to activation of the same protein, which makes structural determination methods particularly difficult. It is also not well understood how the occurrence of phosphorylation and dimerization processes synergize to affect kinase activities. In the present article, we review available structural data and discuss how MD simulations can be used to model conformational transitions of RAF kinase dimers, in both their phosphorylated and unphosphorylated forms.

Dynamics of protein folding: probing the kinetic network of folding-unfolding transitions with experiment and theory.

The problem of spontaneous folding of amino acid chains into highly organized, biologically functional three-dimensional protein structures continues to challenge the modern science. Understanding how proteins fold requires characterization of the underlying energy landscapes as well as the dynamics of the polypeptide chains in all stages of the folding process. In recent years, important advances toward these goals have been achieved owing to the rapidly growing interdisciplinary interest and significant progress in both experimental techniques and theoretical methods. Improvements in the experimental time resolution led to determination of the timescales of the important elementary events in folding, such as formation of secondary structure and tertiary contacts. Sensitive single molecule methods made possible probing the distributions of the unfolded and folded states and following the folding reaction of individual protein molecules. Discovery of proteins that fold in microseconds opened the possibility of atomic-level theoretical simulations of folding and their direct comparisons with experimental data, as well as of direct experimental observation of the barrier-less folding transition. The ultra-fast folding also brought new questions, concerning the intrinsic limits of the folding rates and experimental signatures of barrier-less "downhill" folding. These problems will require novel approaches for even more detailed experimental investigations of the folding dynamics as well as for the analysis of the folding kinetic data. For theoretical simulations of folding, a main challenge is how to extract the relevant information from overwhelmingly detailed atomistic trajectories. New theoretical methods have been devised to allow a systematic approach towards a quantitative analysis of the kinetic network of folding-unfolding transitions between various configuration states of a protein, revealing the transition states and the associated folding pathways at multiple levels, from atomistic to coarse-grained representations. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Computer Simulations of the Dissociation Mechanism of Gleevec from Abl Kinase with Milestoning.

Gleevec (a.k.a., imatinib) is an important anticancer (e.g., chronic myeloid leukemia) chemotherapeutic drug due to its inhibitory interaction with the Abl kinase. Here, we use atomically detailed simulations within the Milestoning framework to study the molecular dissociation mechanism of Gleevec from Abl kinase. We compute the dissociation free energy profile, the mean first passage time for unbinding, and explore the transition state ensemble of conformations. The milestones form a multidimensional network with average connectivity of about 2.93, which is significantly higher than the connectivity for a one-dimensional reaction coordinate. The free energy barrier for Gleevec dissociation is estimated to be ∼10 kcal/mol, and the exit time is ∼55 ms. We examined the transition state conformations using both, the committor and transition function. We show that near the transition state the highly conserved salt bridge K217 and E286 is transiently broken. Together with the calculated free energy profile, these calculations can advance the understanding of the molecular interaction mechanisms between Gleevec and Abl kinase and play a role in future drug design and optimization studies.

Replica Exchange Molecular Dynamics of Diphenylalanine Amyloid Peptides in Electric Fields.

The self-assembling propensity of amyloid peptides such as diphenylalanine (FF) allows them to form ordered, nanoscale structures, with biocompatible properties important for biomedical applications. Moreover, piezoelectric properties allow FF molecules and their aggregates (e.g., FF nanotubes) to be aligned in a controlled way by the application of external electric fields. However, while the behavior of FF nanostructures emerges from the biophysical properties of the monomers, the detailed responses of individual peptides to both temperature and electric fields are not fully understood. Here, we study the temperature-dependent conformational dynamics of FF peptides solvated in explicit water molecules, an environment relevant to biomedical applications, by using an enhanced sampling method, replica exchange molecular dynamics (REMD), in conjunction with applied electric fields. Our simulations highlight and overcome possible artifacts that may occur during the setup of REMD simulations of explicitly solvated peptides in the presence of external electric fields, a problem particularly important in the case of short peptides such as FF. The presence of the external fields could overstabilize certain conformational states in one or more REMD replicas, leading to distortions of the underlying potential energy distributions observed at each temperature. This can be overcome by correcting the REMD initial conditions to include the lower-energy conformations induced by the external field. We show that the converged REMD data can be analyzed using a Markovian description of conformational states and show that a rather complex, 3-state, temperature-dependent conformational dynamics in the absence of electric fields collapses to only one of these states in the presence of the electric fields. These details on the temperature- and electric-field-dependent thermodynamic and kinetic properties of small FF amyloid peptides can be useful in understanding and devising new methods to control their aggregation-prone biophysical properties and, possibly, the structural and biophysical properties of FF molecular nanostructures.

Long-time methods for molecular dynamics simulations: Markov State Models and Milestoning.

Molecular dynamics (MD) studies of biomolecules require the ability to simulate complex biochemical systems with an increasingly larger number of particles and for longer time scales, a problem that cannot be overcome by computational hardware advances alone. A main problem springs from the intrinsically high-dimensional and complex nature of the underlying free energy landscape of most systems, and from the necessity to sample accurately such landscapes for identifying kinetic and thermodynamic states in the configurations space, and for accurate calculations of both free energy differences and of the corresponding transition rates between states. Here, we review and present applications of two increasingly popular methods that allow long-time MD simulations of biomolecular systems that can open a broad spectrum of new studies. A first approach, Markov State Models (MSMs), relies on identifying a set of configuration states in which the system resides sufficiently long to relax and loose the memory of previous transitions, and on using simulations for mapping the underlying complex energy landscape and for extracting accurate thermodynamic and kinetic information. The Markovian independence of the underlying transition probabilities creates the opportunity to increase the sampling efficiency by using sets of appropriately initialized short simulations rather than typically long MD trajectories, which also enhances sampling. This allows MSM-based studies to unveil bio-molecular mechanisms and to estimate free energy barriers with high accuracy, in a manner that is both systematic and relatively automatic, which accounts for their increasing popularity. The second approach presented, Milestoning, targets accurate studies of the ensemble of pathways connecting specific end-states (e.g., reactants and products) in a similarly systematic, accurate and highly automatic manner. Applications presented range from studies of conformational dynamics and binding of amyloid-forming peptides, cell-penetrating peptides and the DFG-flip dynamics in Abl kinase. As highlighted by the increasing number of studies using both methods, we anticipate that they will open new avenues for the investigation of systematic sampling of reactions pathways and mechanisms occurring on longer time scales than currently accessible by purely computational hardware developments.

The transition between active and inactive conformations of Abl kinase studied by rock climbing and Milestoning.

BACKGROUND: Kinases are a family of enzymes that catalyze the transfer of the ɤ-phosphate group from ATP to a protein's residue. Malfunctioning kinases are involved in many health problems such as cardiovascular diseases, diabetes, and cancer. Kinases transitions between multiple conformations of inactive to active forms attracted considerable interest. METHOD: A reaction coordinate is computed for the transition between the active to inactive conformation in Abl kinase with a focus on the DFG-in to DFG-out flip. The method of Rock Climbing is used to construct a path locally, which is subsequently optimized using a functional of the entire path. The discrete coordinate sets along the reaction path are used in a Milestoning calculation of the free energy landscape and the rate of the transition. RESULTS: The estimated transition times are between a few milliseconds and seconds, consistent with simulations of the kinetics and with indirect experimental data. The activation requires the transient dissociation of the salt bridge between Lys271 and Glu286. The salt bridge reforms once the DFG motif is stabilized by a locked conformation of Phe382. About ten residues are identified that contribute significantly to the process and are included as part of the reaction space. CONCLUSIONS: The transition from DFG-in to DFG-out in Abl kinase was simulated using atomic resolution of a fully solvated protein yielding detailed description of the kinetics and the mechanism of the DFG flip. The results are consistent with other computational methods that simulate the kinetics and with some indirect experimental measurements. GENERAL SIGNIFICANCE: The activation of kinases includes a conformational transition of the DFG motif that is important for enzyme activity but is not accessible to conventional Molecular Dynamics. We propose a detailed mechanism for the transition, at a timescale longer than conventional MD, using a combination of reaction path and Milestoning algorithms. The mechanism includes local structural adjustments near the binding site as well as collective interactions with more remote residues.

Computational Opportunities and Challenges in Finding Cyclic Peptide Modulators of Protein-Protein Interactions.

Peptide cyclization can improve stability, conformational constraint, and compactness. However, apart from beta-turn structures, which are well incorporated into cyclic peptides (CPs), many primary peptide structures and functions are markedly altered by cyclization. Accordingly, to mimic linear peptide interfaces with cyclic peptides, it can be beneficial to screen combinatorial cyclic peptide libraries. Computational methods have been developed to screen CPs, but face a number of challenges. Here, we review methods to develop in silico computational libraries, and the potential for screening naturally occurring libraries of CPs. The simplest and most rapid computational pharmacophore methods that estimate peptide three-dimensional structures to be screened versus targets are relatively easy to implement, and while the constraint on structure imposed by cyclization makes them more effective than the same approaches with linear peptides, there are a large number of limiting assumptions. In contrast, full molecular dynamics simulations of cyclic peptide structures not only are costly to implement, but also require careful attention to interpretation, so that not only is the computation time rate limiting, but the interpretation time is also rate limiting due to the analysis of the typically complex underlying conformational space of CPs. A challenge for the field of computational cyclic peptide screening is to bridge this gap effectively. Natural compound libraries of short cyclic peptides, and short cyclized regions of proteins, encoded in the genomes of many organisms present a potential treasure trove of novel functionality which may be screened via combined computational and experimental screening approaches.

Are current molecular dynamics force fields too helical?

Accurate force fields are essential for the success of molecular dynamics simulations. In apparent contrast to the conformational preferences of most force fields, recent NMR experiments suggest that short polyalanine peptides in water populate the polyproline II structure almost exclusively. To investigate this apparent contradiction, with its ramifications for the assessment of molecular force fields and the structure of unfolded proteins, we performed extensive simulations of Ala(5) in water ( approximately 5 micros total time), using twelve different force fields and three different peptide terminal groups. Using either empirical or density-functional-based Karplus relations for the J-couplings, we find that most current force fields do overpopulate the alpha-region, with quantitative results depending on the choice of Karplus relation and on the peptide termini. Even after reweighting to match experiment, we find that Ala(5) retains significant alpha- and beta-populations. In fact, several force fields match the experimental data well before reweighting and have a significant helical population. We conclude that radical changes to the best current force fields are not necessary, based on the NMR data. Nevertheless, experiments on short peptides open the way toward the systematic improvement of current simulation models.

Coarse master equations for peptide folding dynamics.

We construct coarse master equations for peptide folding dynamics from atomistic molecular dynamics simulations. A maximum-likelihood propagator-based method allows us to extract accurate rates for the transitions between the different conformational states of the small helix-forming peptide Ala5. Assigning the conformational states by using transition paths instead of instantaneous molecular coordinates suppresses the effects of fast non-Markovian dynamics. The resulting master equations are validated by comparing their analytical correlation functions with those obtained directly from the molecular dynamics simulations. We find that the master equations properly capture the character and relaxation times of the entire spectrum of conformational relaxation processes. By using the eigenvectors of the transition rate matrix, we are able to systematically coarse-grain the system. We find that a two-state description, with a folded and an unfolded state, roughly captures the slow conformational dynamics. A four-state model, with two folded and two unfolded states, accurately recovers the three slowest relaxation process with time scales between 1.5 and 7 ns. The master equation models not only give access to the slow conformational dynamics but also shed light on the molecular mechanisms of the helix-coil transition.