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.

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.

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.

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.

Structural Modulation of Human Amylin Protofilaments by Naturally Occurring Mutations.

Human islet amyloid polypeptide (hIAPP), also known as amylin, is a 37-amino-acid peptide, co-secreted with insulin, and widely found in fibril form in type-2 diabetes patients. By using all-atom molecular dynamics simulations, we study hIAPP fibril segments (i.e., fibrillar oligomers) formed with sequences of naturally occurring variants from cat, rat, and pig, presenting different aggregation propensities. We characterize the effect of mutations on the structural dynamics of solution-formed hIAPP fibril models built from solid-state NMR data. Results from this study are in agreement with experimental observations regarding their respective relative aggregation propensities. We analyze in detail the specific structural characteristics and infer mechanisms that modulate the conformational stability of amylin fibrils. Results provide a platform for further studies and the design of new drugs that could interfere with amylin aggregation and its cytotoxicity. One particular mutation, N31K, has fibril-destabilizing properties, and could potentially improve the solubility of therapeutic amylin analogs.

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.

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.

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.

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.

Nanoscale Piezoelectric Properties of Self-Assembled Fmoc-FF Peptide Fibrous Networks.

Fibrous peptide networks, such as the structural framework of self-assembled fluorenylmethyloxycarbonyl diphenylalanine (Fmoc-FF) nanofibrils, have mechanical properties that could successfully mimic natural tissues, making them promising materials for tissue engineering scaffolds. These nanomaterials have been determined to exhibit shear piezoelectricity using piezoresponse force microscopy, as previously reported for FF nanotubes. Structural analyses of Fmoc-FF nanofibrils suggest that the observed piezoelectric response may result from the noncentrosymmetric nature of an underlying ß-sheet topology. The observed piezoelectricity of Fmoc-FF fibrous networks is advantageous for a range of biomedical applications where electrical or mechanical stimuli are required.

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.

Conformational dynamics of human IAPP monomers.

We study the conformational dynamics of the human Islet Amyloid Polypeptide (hIAPP) molecule - a 37 residue-long peptide associated to type 2 diabetes - using molecular dynamics (MD) simulations. We identify partially structured conformational states of the hIAPP monomer, categorized by both end-to-end distance and secondary structure, as suggested by previous experimental and computational studies. The MD trajectories of hIAPP are analyzed using data-driven methods, in particular principal component analysis, in order to identify preferred conformational states of the amylin monomer and to discuss their relative stability as compared to corresponding states in the amylin dimer. These potential hIAPP conformational states could be further tested and described experimentally, or in conjunction with modern computational analysis tools such as Markov state-based methods for extracting kinetics and thermodynamics from atomistic MD trajectories.

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.

Peptide folding kinetics from replica exchange molecular dynamics.

We show how accurate kinetic information, such as the rates of protein folding and unfolding, can be extracted from replica-exchange molecular dynamics (REMD) simulations. From the brief and continuous trajectory segments between replica exchanges, we estimate short-time propagators in conformation space and use them to construct a master equation. For a helical peptide in explicit water, we determine the rates of transitions both locally between microscopic conformational states and globally for folding and unfolding. We show that accurate rates in the approximately 1(100 ns) to approximately 1(1 ns) range can be obtained from REMD with exchange times of 5 ps , in excellent agreement with results from long equilibrium molecular dynamics.

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.