Structure of a transient intermediate for GTP hydrolysis by ras

Unraveling the Significance of Mg2+ Dependency and Nucleotide Binding Specificity of H-RAS

RAS is a small GTPase and acts as a binary molecular switch; the transition from its active to inactive state plays a crucial role in various cell signaling processes. Molecular dynamics simulations at the atomistic level suggest that the absence of cofactor Mg2+ ion generally leads to pronounced structural changes in the Switch-I than Switch-II regions and assists GTP binding. The presence of the Mg2+ ion also restricts the rotation of ϒ phosphate and enhances the hydrolysis rate of GTP. Further, the simulations reveal that the stability of the protein is almost uncompromised when Mg2+ is replaced with Zn2+ and not the Ca2+ ion. The specificity of H-RAS to GTP was evaluated by substituting with ATP and CTP, which indicates that the binding pocket tolerates purine bases over pyrimidine bases. However, the D119 residue specifically interacts with the guanine base and serves as one of the primary interactions that leads to the selectivity of GTP over ATP. The ring displacement of 32Y serves as gate dynamics in H-RAS which are important for its interaction with GAP for the nucleotide exchange and is restricted in the presence of ATP. Finally, the point mutations 61, 16, and 32 influence the structural changes, specifically in the Switch-II region, which are expected to impact the GTP hydrolysis and thus are termed oncogenic mutations.

A RhoA structure with switch II flipped outward revealed the conformational dynamics of switch II region

Small GTPase RhoA switches from GTP-bound state to GDP-bound state by hydrolyzing GTP, which is accelerated by GTPases activating proteins (GAPs). However, less study of RhoA structural dynamic changes was conducted during this process, which is essential for understanding the molecular mechanism of GAP dissociation. Here, we solved a RhoA structure in GDP-bound state with switch II flipped outward. Because lacking the intermolecular interactions with guanine nucleotide, we proposed this conformation of RhoA could be an intermediate after GAP dissociation. Further molecular dynamics simulations found the conformational changes of switch regions are indeed existing in RhoA and involved in the regulation of GAP dissociation and GEF recognition. Besides, the guanine nucleotide binding pocket extended to switch II region, indicating a potential "druggable" cavity for RhoA. Taken together, our study provides a deeper understanding of the dynamic properties of RhoA switch regions and highlights the direction for future drug development.

Equilibria between conformational states of the Ras oncogene protein revealed by high pressure crystallography

In this work, we experimentally investigate the allosteric transitions between conformational states on the Ras oncogene protein using high pressure crystallography. Ras protein is a small GTPase involved in central regulatory processes occurring in multiple conformational states. Ras acts as a molecular switch between active GTP-bound, and inactive GDP-bound states, controlling essential signal transduction pathways. An allosteric network of interactions between the effector binding regions and the membrane interacting regions is involved in Ras cycling. The conformational states which coexist simultaneously in solution possess higher Gibbs free energy than the ground state. Equilibria between these states can be shifted by applying pressure favouring conformations with lower partial molar volume, and has been previously analyzed by high-pressure NMR spectroscopy. High-pressure macromolecular crystallography (HPMX) is a powerful tool perfectly complementary to high-pressure NMR, allowing characterization at the molecular level with a high resolution the different allosteric states involved in the Ras cycling. We observe a transition above 300 MPa in the crystal leading to more stable conformers. Thus, we compare the crystallographic structures of Ras(wt)·Mg²⁺·GppNHp and Ras(D33K)·Mg²⁺·GppNHp at various high hydrostatic pressures. This gives insight into per-residue descriptions of the structural plasticity involved in allosteric equilibria between conformers. We have mapped out at atomic resolution the different segments of Ras protein which remain in the ground-state conformation or undergo structural changes, adopting excited-energy conformations corresponding to transient intermediate states. Such in crystallo phase transitions induced by pressure open the possibility to finely explore the structural determinants related to switching between Ras allosteric sub-states without any mutations nor exogenous partners.

The equilibria between structural states induced by pressure within the crystal structure of Ras are illustrated with different colors corresponding to different Ras substates.

Conformational Features of Ras: Key Hydrogen-Bonding Interactions of Gln61 in the Intermediate State during GTP Hydrolysis

The Ras protein is one of the most important drug targets for battling cancers. To effectively design novel drugs of Ras, we characterize here its conformational ensembles for the hydrolysis intermediate state RasGDP·Pi and the product state RasGDP by extensive replica-exchange molecular dynamics simulations. Several substates for RasGDP·Pi have been identified, while structural analyses have revealed an unrecognized hydrogen-bonding network that stabilizes the hydrolysis intermediate state. More interestingly, Gln61, which is involved in numerous oncogenic mutations, was found to be engaged in this hydrogen-bonding network, adopting a specific conformation that always points to Pi in contrast to that in the RasGTP state. The simulations also reveal that RasGDP has more than one substate, suggesting a conformational selection mechanism for the interaction between Ras and the guanine nucleotide exchange factors (GEFs). These findings offer new opportunities for the drug design of Ras by stabilizing the hydrolysis intermediate or disrupting its interaction with the GEFs.

Mechanistic insights into the effect of phosphorylation on Ras conformational dynamics and its interactions with cell signaling proteins

Ras undergoes interconversion between the active GTP-bound state and the inactive GDP-bound state. This GTPase cycle, which controls the activities of Ras, is accelerated by Ras GTPase-activating proteins (GAPs) and guanine nucleotide exchange factors (SOS). Oncogenic Ras mutations could affect the GTPase cycle and impair Ras functions. Additionally, Src-induced K-Ras Y32/64 dual phosphorylation has been reported to disrupt GTPase cycle and hinder Ras downstream signaling. However, the underlying mechanisms remain unclear. To address this, we performed molecular dynamics simulations (~30 μs in total) on unphosphorylated and phosphorylated K-Ras4B in GTP- and GDP-bound states, and on their complexes with GTPase cycle regulators (GAP and SOS) and the effector protein Raf. We found that K-Ras4B dual phosphorylation mainly alters the conformation at the nucleotide binding site and creates disorder at the catalytic site, resulting in the enlargement of GDP binding pocket and the retard of Ras-GTP intrinsic hydrolysis. We observed phosphorylation-induced shift in the distribution of Ras-GTP inactive-active sub-states and recognized potential druggable pockets in the phosphorylated Ras-GTP. Moreover, decreased catalytic competence or signal delivery abilities due to reduced binding affinities and/or distorted catalytic conformations of GAP, SOS and Raf were observed. In addition, the allosteric pathway from Ras/Raf interface to the distal Raf L4 loop was compromised by Ras phosphorylation. These results reveal the mechanisms by which phosphorylation influences the intrinsic or GAP/SOS catalyzed transformations between GTP- and GDP-bound states of Ras and its signal transduction to Raf. Our findings project Ras phosphorylation as a target for cancer drug discovery.

Unveiling conformational dynamics changes of H-Ras induced by mutations based on accelerated molecular dynamics

Uncovering molecular basis with regard to the conformational change of two switches I and II in the GppNHp (GNP)-bound H-Ras is highly significant for the understanding of Ras signaling. For this purpose, accelerated molecular dynamics (aMD) simulations and principal component (PC) analysis are integrated to probe the effect of mutations G12V, T35S and Q61K on conformational transformation between two switches of the GNP-bound H-Ras. The RMSF and cross-correlation analyses suggest that three mutations exert a vital effect on the flexibility and internal dynamics of two switches in the GNP-bound H-Ras. The results stemming from PC analysis indicate that two switches in the GNP-bound WT H-Ras tend to form a closed state in most conformations, while those in the GNP-bound mutated H-Ras display transformation between different states. This conclusion is further supported by free energy landscapes constructed by using the distances of residues 12 away from 35 and 35 away from 61 as reaction coordinates and different experimental studies. Interaction scanning is performed on aMD trajectories and the information shows that conformational transformations of two switches I and II induced by mutations extremely affect the GNP-residue interactions. Meanwhile, the scanning results also signify that residues G15, A18, F28, K117, A146 and K147 form stable contacts with GNP, while residues D30, E31, Y32, D33, P34 and E62 in two switches I and II produce unstable contacts with GNP. This study not only reveals dynamic behavior changes of two switches in H-Ras induced by mutations, but also unveils general principles and mechanisms with regard to functional conformational changes of H-Ras.

Comparative MD simulations and advanced analytics based studies on wild-type and hot-spot mutant A59G HRas

The Ras family of proteins is known to play an important role in cellular signal transduction. The oncoprotein Ras is also found to be mutated in ~90% of the pancreatic cancers, of which G12V, G13V, A59G and Q61L are the known hot-spot mutants. These ubiquitous proteins fall in the family of G-proteins, and hence switches between active GTP bound and inactive GDP bound states, which is hindered in most of its oncogenic mutant counterparts. Moreover, Ras being a GTPase has an intrinsic property to hydrolyze GTP to GDP, which is obstructed due to mutations and lends the mutants stuck in constitutively active state leading to oncogenic behavior. In this regard, the present study aims to understand the dynamics involved in the hot-spot mutant A59G-Ras using long 10μs classical MD simulations (5μs for each of the wild-type and mutant systems) and comparing the same with its wild-type counterpart. Advanced analytics using Markov State Model (MSM) based approach has been deployed to comparatively understand the transition path for the wild-type and mutant systems. Roles of crucial residues like Tyr32, Gln61 and Tyr64 have also been established using multivariate PCA analyses. Furthermore, this multivariate PCA analysis also provides crucial features which may be used as reaction coordinates for biased simulations for further studies. The absence of formation of pre-hydrolysis network is also reported for the mutant conformation, using the distance-based analyses (between crucial residues) of the conserved regions. The implications of this study strengthen the hypothesis that the disruption of the pre-hydrolysis network in the mutant A59G ensemble might lead to permanently active oncogenic conformation in the mutant conformers.

Review: Precision medicine and driver mutations: Computational methods, functional assays and conformational principles for interpreting cancer drivers

At the root of the so-called precision medicine or precision oncology, which is our focus here, is the hypothesis that cancer treatment would be considerably better if therapies were guided by a tumor’s genomic alterations. This hypothesis has sparked major initiatives focusing on whole-genome and/or exome sequencing, creation of large databases, and developing tools for their statistical analyses—all aspiring to identify actionable alterations, and thus molecular targets, in a patient. At the center of the massive amount of collected sequence data is their interpretations that largely rest on statistical analysis and phenotypic observations. Statistics is vital, because it guides identification of cancer-driving alterations. However, statistics of mutations do not identify a change in protein conformation; therefore, it may not define sufficiently accurate actionable mutations, neglecting those that are rare. Among the many thematic overviews of precision oncology, this review innovates by further comprehensively including precision pharmacology, and within this framework, articulating its protein structural landscape and consequences to cellular signaling pathways. It provides the underlying physicochemical basis, thereby also opening the door to a broader community.

Tyrosyl phosphorylation of KRAS stalls GTPase cycle via alteration of switch I and II conformation

Deregulation of the RAS GTPase cycle due to mutations in the three RAS genes is commonly associated with cancer development. Protein tyrosine phosphatase SHP2 promotes RAF-to-MAPK signaling pathway and is an essential factor in RAS-driven oncogenesis. Despite the emergence of SHP2 inhibitors for the treatment of cancers harbouring mutant KRAS, the mechanism underlying SHP2 activation of KRAS signaling remains unclear. Here we report tyrosyl-phosphorylation of endogenous RAS and demonstrate that KRAS phosphorylation via Src on Tyr32 and Tyr64 alters the conformation of switch I and II regions, which stalls multiple steps of the GTPase cycle and impairs binding to effectors. In contrast, SHP2 dephosphorylates KRAS, a process that is required to maintain dynamic canonical KRAS GTPase cycle. Notably, Src- and SHP2-mediated regulation of KRAS activity extends to oncogenic KRAS and the inhibition of SHP2 disrupts the phosphorylation cycle, shifting the equilibrium of the GTPase cycle towards the stalled ‘dark state’.

Deregulation of the RAS GTPase cycle due to mutations in RAS genes is commonly associated with cancer development. Here authors use NMR and mass spectrometry to shows that KRAS phosphorylation via Src alters the conformation of switch I and II regions and thereby impacts the GTPase cycle.

Atypical Switch-I Arginine plays a catalytic role in GTP hydrolysis by Rab21 from Entamoeba histolytica

Entamoeba histolytica, the causative agent of amoebic dysentery, liver abscess and colitis, exploits its vesicular trafficking machinery for survival and virulence. Rab family of small GTPases play a key role in the vesicular transport by undergoing the GTP/GDP cycle which is central to the biological processes. Amoebic genome encodes several atypical Rab GTPases which are unique due to absence of conserved sequence motif(s) or atypical residues in their catalytic site [Saito-Nakano et al., 2005 ]. Previously, EhRab21 has been reported to involve in amoebic invasion and migration [Emmanuel et al., 2015 ]. The conserved Glutamine of switch-II region is universally accepted to be crucial for GTP hydrolysis. Mutations that reduce the sidechain polarity of Glutamine render the protein GTPase activity deficient [Krengel et al., 1990]. Here, we report a catalytic role of atypical switch-I Arginine (R36) in intrinsic GTP hydrolysis catalysed by EhRab21. Unlike the GTPase activity deficient QL mutants, the GTPase activity of EhRab21Q64L was found to be marginally enhanced compared to the wild-type protein. Although EhRab21R36L mutant showed normal GTPase activity, the double mutant (R36L/Q64L) was found to be GTPase deficient. Thus, EhRab21 is a unique member of small GTPase family in which an atypical switch-I Arginine is capable of driving GTP hydrolysis independent of the conserved switch-II Glutamine.

Structural insight into the rearrangement of the switch I region in GTP-bound G12A K-Ras

K-Ras, a molecular switch that regulates cell growth, apoptosis and metabolism, is activated when it undergoes a conformation change upon binding GTP and is deactivated following the hydrolysis of GTP to GDP. Hydrolysis of GTP in water is accelerated by coordination to K-Ras, where GTP adopts a high-energy conformation approaching the transition state. The G12A mutation reduces intrinsic K-Ras GTP hydrolysis by an unexplained mechanism. Here, crystal structures of G12A K-Ras in complex with GDP, GTP, GTPγS and GppNHp, and of Q61A K-Ras in complex with GDP, are reported. In the G12A K-Ras-GTP complex, the switch I region undergoes a significant reorganization such that the Tyr32 side chain points towards the GTP-binding pocket and forms a hydrogen bond to the GTP γ-phosphate, effectively stabilizing GTP in its precatalytic state, increasing the activation energy required to reach the transition state and contributing to the reduced intrinsic GTPase activity of G12A K-Ras mutants.

Mechanistic Insights into the Differential Catalysis by RheB and Its Mutants: Y35A and Y35A-D65A

RheB GTPase is a Ras-related molecular switch, which regulates the mTOR signaling pathway by cycling between the active [guanosine triphosphate (GTP)] state and inactive [guanine diphosphate (GDP)] state. Impairment of GTPase activity because of mutations in several small GTPases is known to be associated with several cancers. The conventional GTPase mechanism such as in H-Ras requires a conserved glutamine (Q64) in the switch-II region of RheB to align the catalytic water molecule for efficient GTP hydrolysis. The conformation of this conserved glutamine is different in RheB, resulting in an altered conformation of the entire switch-II region. Studies on the atypical switch-II conformation in RheB revealed a distinct, noncanonical mode of GTP hydrolysis. An RheB mutant Y35A was previously shown to exclusively enhance the intrinsic GTPase activity of RheB, whereas the Y35A-D65A double mutant was shown to reduce the elevated GTPase activity. Here, we have used all-atom molecular dynamics (MD) simulations for comprehensive understanding of the conformational dynamics associated with the fast (Y35A) and slow (Y35A-D65A) hydrolyzing mutants of RheB. Using a combination of starting models from PDB structures and in-silico generated mutant structures, we discuss the observed conformational deviations in wild type (WT) versus mutants. Our results show that a number of interactions of RheB with phosphates of GTP as well as Mg²⁺ are destabilized in Y35A mutant in the switch-I region. We report distinct water dynamics at the active site of WT and mutants. Furthermore, principal component analysis showed significant differences in the conformational space sampled by the WT and mutants. Our observations provide improved understanding of the noncanonical GTP hydrolysis mechanism adopted by RheB and its modulation by Y35A and Y35A-D65A mutants.

Electric-Field-Induced Protein Translocation via a Conformational Transition in SecDF: An MD Study

SecDF is an important component of the Sec protein translocation machinery embedded in the bacterial membrane, which is associated with many functions, such as stabilizing other Sec translocon components within the membrane, maintaining the transmembrane (TM) potential, and facilitating the ATP-independent stage of the translocation mechanism. Related studies suggest that SecDF undergoes functionally important conformational changes that involve mainly its P1-head domain and that these changes are coupled with the proton motive force (Δp). However, there still is not a clear understanding of how SecDF functions, its exact role in the translocation machinery, and how its function is related to Δp. Here, using all-atom molecular dynamics simulations combined with umbrella sampling, we study the P1-head conformational change and how it is coupled to the proton motive force. We report potentials of mean force along a root-mean-square-distance-based reaction coordinate obtained in the presence and absence of the TM electrical potential. Our results show that the interaction of the P1 domain dipole moment with the TM electrical field considerably lowers the free-energy barrier in the direction of F-form to I-form transition.

The small GTPases K-Ras, N-Ras, and H-Ras have distinct biochemical properties determined by allosteric effects

H-Ras, K-Ras, and N-Ras are small GTPases that are important in the control of cell proliferation, differentiation, and survival, and their mutants occur frequently in human cancers. The G-domain, which catalyzes GTP hydrolysis and mediates downstream signaling, is 95% conserved between the Ras isoforms. Because of their very high sequence identity, biochemical studies done on H-Ras have been considered representative of all three Ras proteins. We show here that this is not a valid assumption. Using enzyme kinetic assays under identical conditions, we observed clear differences between the three isoforms in intrinsic catalysis of GTP by Ras in the absence and presence of the Ras-binding domain (RBD) of the c-Raf kinase protein (Raf-RBD). Given their identical active sites, isoform G-domain differences must be allosteric in origin, due to remote isoform-specific residues that affect conformational states. We present the crystal structure of N-Ras bound to a GTP analogue and interpret the kinetic data in terms of structural features specific for H-, K-, and N-Ras.

Differentiating the pre-hydrolysis states of wild-type and A59G mutant HRas: An insight through MD simulations

The most representative member of the Ras subfamily is its HRas isoform. Ras proteins being GTPases, possess an intrinsic activity to hydrolyze the GTP molecule to GDP. During the transition phases, between active and inactive states, P-loop and switch regions show maximum variations. Various hot-spot Ras mutants (G12V, A59G, Q61L etc) have been reported, that limit the protein's conformation in the permanent active state. In the present study, we aim to explore the structural dynamics of one such crucial mutant of Ras namely A59G which belongs to the conserved Switch II region of the protein. Approximately ∼15μs of Classical Molecular Dynamics (CMD) simulations have been carried out on the mutant and wild-type complexes. Further, a metadynamics simulation of 500ns was also carried out, which suggests an energy barrier of ∼9.56kcal/mol between wild-type and mutant conformation. We demonstrate the role of water molecule in maintaining the required interaction networks in the pre-hydrolysis state, its impact on A59G mutation, distinct orientation of the Gln61 residue in two conformations, disruption of crucial Gly60 and γ phosphate and the change in the Switch II region. The outcome of our study captures the pre-hydrolysis state of the HRas protein. It also establishes the fact that this mutation makes the movement of Switch II region and the conserved DXXGQ motif highly constrained, which is known to be an important requirement for hydrolysis. This suggests that the A59G mutation may decrease the rate of intrinsic hydrolysis as well as GAP-mediated hydrolysis.

Distinct dynamics and interaction patterns in H- and K-Ras oncogenic P-loop mutants

Despite years of study, the structural or dynamical basis for the differential reactivity and oncogenicity of Ras isoforms and mutants remains unclear. In this study, we investigated the effects of amino acid variations on the structure and dynamics of wild type and oncogenic mutants G12D, G12V, and G13D of H- and K-Ras proteins. Based on data from µs-scale molecular dynamics simulations, we show that the overall structure of the proteins remains similar but there are important differences in dynamics and interaction networks. We identified differences in residue interaction patterns around the canonical switch and distal loop regions, and persistent sodium ion binding near the GTP particularly in the G13D mutants. Our results also suggest that different Ras variants have distinct local structural features and interactions with the GTP, variations that have the potential to affect GTP release and hydrolysis. Furthermore, we found that H-Ras proteins and particularly the G12V and G13D variants are significantly more flexible than their K-Ras counterparts. Finally, while most of the simulated proteins sampled the effector-interacting state 2 conformational state, G12V and G13D H-Ras adopted an open switch state 1 conformation that is defective in effector interaction. These differences have implications for Ras GTPase activity, effector or exchange factor binding, dimerization and membrane interaction. Proteins 2017; 85:1618-1632. © 2017 Wiley Periodicals, Inc.

The Structural Basis of Oncogenic Mutations G12, G13 and Q61 in Small GTPase K-Ras4B

Ras mediates cell proliferation, survival and differentiation. Mutations in K-Ras4B are predominant at residues G12, G13 and Q61. Even though all impair GAP-assisted GTP → GDP hydrolysis, the mutation frequencies of K-Ras4B in human cancers vary. Here we aim to figure out their mechanisms and differential oncogenicity. In total, we performed 6.4 μs molecular dynamics simulations on the wild-type K-Ras4B (K-Ras4B(WT)-GTP/GDP) catalytic domain, the K-Ras4B(WT)-GTP-GAP complex, and the mutants (K-Ras4B(G12C/G12D/G12V)-GTP/GDP, K-Ras4B(G13D)-GTP/GDP, K-Ras4B(Q61H)-GTP/GDP) and their complexes with GAP. In addition, we simulated 'exchanged' nucleotide states. These comprehensive simulations reveal that in solution K-Ras4B(WT)-GTP exists in two, active and inactive, conformations. Oncogenic mutations differentially elicit an inactive-to-active conformational transition in K-Ras4B-GTP; in K-Ras4B(G12C/G12D)-GDP they expose the bound nucleotide which facilitates the GDP-to-GTP exchange. These mechanisms may help elucidate the differential mutational statistics in K-Ras4B-driven cancers. Exchanged nucleotide simulations reveal that the conformational transition is more accessible in the GTP-to-GDP than in the GDP-to-GTP exchange. Importantly, GAP not only donates its R789 arginine finger, but stabilizes the catalytically-competent conformation and pre-organizes catalytic residue Q61; mutations disturb the R789/Q61 organization, impairing GAP-mediated GTP hydrolysis. Together, our simulations help provide a mechanistic explanation of key mutational events in one of the most oncogenic proteins in cancer.

Ras Conformational Ensembles, Allostery, and Signaling

Ras proteins are classical members of small GTPases that function as molecular switches by alternating between inactive GDP-bound and active GTP-bound states. Ras activation is regulated by guanine nucleotide exchange factors that catalyze the exchange of GDP by GTP, and inactivation is terminated by GTPase-activating proteins that accelerate the intrinsic GTP hydrolysis rate by orders of magnitude. In this review, we focus on data that have accumulated over the past few years pertaining to the conformational ensembles and the allosteric regulation of Ras proteins and their interpretation from our conformational landscape standpoint. The Ras ensemble embodies all states, including the ligand-bound conformations, the activated (or inactivated) allosteric modulated states, post-translationally modified states, mutational states, transition states, and nonfunctional states serving as a reservoir for emerging functions. The ensemble is shifted by distinct mutational events, cofactors, post-translational modifications, and different membrane compositions. A better understanding of Ras biology can contribute to therapeutic strategies.

Mapping the Conformation Space of Wildtype and Mutant H-Ras with a Memetic, Cellular, and Multiscale Evolutionary Algorithm

An important goal in molecular biology is to understand functional changes upon single-point mutations in proteins. Doing so through a detailed characterization of structure spaces and underlying energy landscapes is desirable but continues to challenge methods based on Molecular Dynamics. In this paper we propose a novel algorithm, SIfTER, which is based instead on stochastic optimization to circumvent the computational challenge of exploring the breadth of a protein's structure space. SIfTER is a data-driven evolutionary algorithm, leveraging experimentally-available structures of wildtype and variant sequences of a protein to define a reduced search space from where to efficiently draw samples corresponding to novel structures not directly observed in the wet laboratory. The main advantage of SIfTER is its ability to rapidly generate conformational ensembles, thus allowing mapping and juxtaposing landscapes of variant sequences and relating observed differences to functional changes. We apply SIfTER to variant sequences of the H-Ras catalytic domain, due to the prominent role of the Ras protein in signaling pathways that control cell proliferation, its well-studied conformational switching, and abundance of documented mutations in several human tumors. Many Ras mutations are oncogenic, but detailed energy landscapes have not been reported until now. Analysis of SIfTER-computed energy landscapes for the wildtype and two oncogenic variants, G12V and Q61L, suggests that these mutations cause constitutive activation through two different mechanisms. G12V directly affects binding specificity while leaving the energy landscape largely unchanged, whereas Q61L has pronounced, starker effects on the landscape. An implementation of SIfTER is made available at http://www.cs.gmu.edu/~ashehu/?q=OurTools. We believe SIfTER is useful to the community to answer the question of how sequence mutations affect the function of a protein, when there is an abundance of experimental structures that can be exploited to reconstruct an energy landscape that would be computationally impractical to do via Molecular Dynamics.

Specific cancer-associated mutations in the switch III region of Ras increase tumorigenicity by nanocluster augmentation

Hotspot mutations of Ras drive cell transformation and tumorigenesis. Less frequent mutations in Ras are poorly characterized for their oncogenic potential. Yet insight into their mechanism of action may point to novel opportunities to target Ras. Here, we show that several cancer-associated mutations in the switch III region moderately increase Ras activity in all isoforms. Mutants are biochemically inconspicuous, while their clustering into nanoscale signaling complexes on the plasma membrane, termed nanocluster, is augmented. Nanoclustering dictates downstream effector recruitment, MAPK-activity, and tumorigenic cell proliferation. Our results describe an unprecedented mechanism of signaling protein activation in cancer.

DOI:http://dx.doi.org/10.7554/eLife.08905.001

Cancer is a disease that develops when cells within the body acquire genetic mutations that allow them to grow and divide rapidly. Many human cancers have mutations in a gene that encodes a protein called Ras, which promotes cell growth and division by controlling the activities of other proteins.

Ras congregates at the membrane that surrounds the cell and can assemble into clusters (called nanoclusters) that each contain six to eight Ras proteins. The tight packing of the proteins in these nanoclusters increases the amount of Ras in the membrane locally, which allows Ras to interact with other proteins more efficiently to promote growth and cell division.

In normal cells, other proteins control when Ras is active. However, in many cancer cells, Ras is active all the time due to mutations that occur in three ‘hotspots’ within its gene. Other mutations in the gene that encodes Ras are also found in cancer cells, but these are less common and it is not clear how they alter the activity of the protein.

Here, Solman et al. used microscopy and biochemical techniques to study the effects of some of the less common mutations on Ras activity in human cells. The experiments show that several mutations that alter a region of Ras called the ‘switch III region’ moderately increase the activity of Ras. The mutations probably alter the way that Ras sits in the membrane, which in turn changes the way it interacts with other proteins and the membrane so that more Ras nanoclusters form.

Solman et al.'s findings reveal a new way that Ras can be activated in cancer cells. The next challenge is to develop drugs that block the formation of Ras nanoclusters and to find out if they have the potential to be used to treat cancer.

DOI:http://dx.doi.org/10.7554/eLife.08905.002

Differential dynamics of RAS isoforms in GDP- and GTP-bound states

RAS subfamily proteins regulates cell growth promoting signaling processes by cycling between active (GTP-bound) and inactive (GDP-bound) states. Different RAS isoforms, though structurally similar, exhibit functional specificity and are associated with different types of cancers and developmental disorders. Understanding the dynamical differences between the isoforms is crucial for the design of inhibitors that can selectively target a particular malfunctioning isoform. In this study, we provide a comprehensive comparison of the dynamics of all the three RAS isoforms (HRAS, KRAS, and NRAS) using extensive molecular dynamics simulations in both the GDP- (total of 3.06 μs) and GTP-bound (total of 2.4 μs) states. We observed significant differences in the dynamics of the isoforms, which rather interestingly, varied depending on the type of the nucleotide bound and the simulation temperature. Both SwitchI (Residues 25-40) and SwitchII (Residues 59-75) differ significantly in their flexibility in the three isoforms. Furthermore, Principal Component Analysis showed that there are differences in the conformational space sampled by the GTP-bound RAS isoforms. We also identified a previously unreported pocket, which opens transiently during MD simulations, and can be targeted to regulate nucleotide exchange reaction or possibly interfere with membrane localization. Further, we present the first simulation study showing GDP destabilization in the wild-type RAS protein. The destabilization of GDP/GTP occurred only in 1/50 simulations, emphasizing the need of guanine nucleotide exchange factors (GEFs) to accelerate such an extremely unfavorable process. This observation along with the other results presented in this article further support our previously hypothesized mechanism of GEF-assisted nucleotide exchange.

GTP-specific fab fragment-based GTPase activity assay

GTPases are central cellular signaling proteins, which cycle between a GDP-bound inactive and a GTP-bound active conformation in a controlled manner. Ras GTPases are frequently mutated in cancer and so far only few experimental inhibitors exist. The most common methods for monitoring GTP hydrolysis rely on luminescent GDP- or GTP-analogs. In this study, the first GTP-specific Fab fragment and its application are described. We selected Fab fragments using the phage display technology. Six Fab fragments were found against 2'/3'-GTP-biotin and 8-GTP-biotin. Selected antibody fragments allowed specific detection of endogenous, free GTP. The most potent Fab fragment (2A4(GTP)) showed over 100-fold GTP-specificity over GDP, ATP, or CTP and was used to develop a heterogeneous time-resolved luminescence based assay for the monitoring of GTP concentration. The method allows studying the GEF dependent H-Ras activation (GTP binding) and GAP-catalyzed H-Ras deactivation (GTP hydrolysis) at nanomolar protein concentrations.

Mechanism of the exchange reaction in HRAS from multiscale modeling

HRAS regulates cell growth promoting signaling processes by cycling between active (GTP-bound) and inactive (GDP-bound) states. Understanding the transition mechanism is central for the design of small molecules to inhibit the formation of RAS-driven tumors. Using a multiscale approach involving coarse-grained (CG) simulations, all-atom classical molecular dynamics (CMD; total of 3.02 µs), and steered molecular dynamics (SMD) in combination with Principal Component Analysis (PCA), we identified the structural features that determine the nucleotide (GDP) exchange reaction. We show that weakening the coupling between the SwitchI (residues 25–40) and SwitchII (residues 59–75) accelerates the opening of SwitchI; however, an open conformation of SwitchI is unstable in the absence of guanine nucleotide exchange factors (GEFs) and rises up towards the bound nucleotide to close the nucleotide pocket. Both I21 and Y32, play a crucial role in SwitchI transition. We show that an open SwitchI conformation is not necessary for GDP destabilization but is required for GDP/Mg escape from the HRAS. Further, we present the first simulation study showing displacement of GDP/Mg away from the nucleotide pocket. Both SwitchI and SwitchII, delays the escape of displaced GDP/Mg in the absence of GEF. Based on these results, a model for the mechanism of GEF in accelerating the exchange process is hypothesized.

A homogeneous quenching resonance energy transfer assay for the kinetic analysis of the GTPase nucleotide exchange reaction

A quenching resonance energy transfer (QRET) assay for small GTPase nucleotide exchange kinetic monitoring is demonstrated using nanomolar protein concentrations. Small GTPases are central signaling proteins in all eukaryotic cells acting as a "molecular switches" that are active in the GTP-state and inactive in the GDP-state. GTP-loading is highly regulated by guanine nucleotide exchange factors (GEFs). In several diseases, most prominently cancer, this process in misregulated. The kinetics of the nucleotide exchange reaction reports on the enzymatic activity of the GEF reaction system and is, therefore, of special interest. We determined the nucleotide exchange kinetics using europium-labeled GTP (Eu-GTP) in the QRET assay for small GTPases. After GEF catalyzed GTP-loading of a GTPase, a high time-resolved luminescence signal was found to be associated with GTPase bound Eu-GTP, whereas the non-bound Eu-GTP fraction was quenched by soluble quencher. The association kinetics of the Eu-GTP was measured after GEF addition, whereas the dissociation kinetics could be determined after addition of unlabeled GTP. The resulting association and dissociation rates were in agreement with previously published values for H-Ras(Wt), H-Ras(Q61G), and K-Ras(Wt), respectively. The broader applicability of the QRET assay for small GTPases was demonstrated by determining the kinetics of the Ect2 catalyzed RhoA(Wt) GTP-loading. The QRET assay allows the use of nanomolar protein concentrations, as more than 3-fold signal-to-background ratio was achieved with 50 nM GTPase and GEF proteins. Thus, small GTPase exchange kinetics can be efficiently determined in a HTS compatible 384-well plate format.

Probing the wild-type HRas activation mechanism using steered molecular dynamics, understanding the energy barrier and role of water in the activation

Ras is one of the most common oncogenes in human cancers. It belongs to a family of GTPases that functions as binary conformational switches by timely switching of their conformations from GDP to GTP and vice versa. It attains the final active state structure via an intermediate GTP-bound state. The transition between these states is a millisecond-time-scale event. This makes studying this mechanism beyond the scope of classical molecular dynamics. In the present study, we describe the activation pathway of the HRas protein complex along the distance-based reaction coordinate using steered molecular dynamics. Approximately ~720 ns of MD simulations using CMD and SMD was performed. We demonstrated the change in orientation and arrangement of the two switch regions and the role of various hydrogen bonds during the activation process. The weighted histogram analysis method was also performed, and the potential of mean force was calculated between the inactive and active via the intermediate state (state 1) of HRas. The study indicates that water seems to play a crucial role in the activation process and to transfer the HRas protein from its intermediate state to the fully active state. The implications of our study hereby suggest that the HRas activation mechanism is a multistep process. It starts from the inactive state to an intermediate state 1 followed by trapping of water molecules and flipping of the Thr35 residue to form a fully active state (state 2). This state 2 also comprises Gly60, Thr35, GTP, Mg(2+) and water-forming stable interactions.

Lessons from computer simulations of Ras proteins in solution and in membrane

BACKGROUND

A great deal has been learned over the last several decades about the function of Ras proteins in solution and membrane environments. While much of this knowledge has been derived from a plethora of experimental techniques, computer simulations have also played a substantial role.

SCOPE OF REVIEW

Our goal here is to summarize the contribution of molecular simulations to our current understanding of normal and aberrant Ras function. We focus on lessons from molecular dynamics simulations in aqueous and membrane environments.

MAJOR CONCLUSIONS

The central message is that a close interaction between theory and simulation on the one hand and cell-biological, spectroscopic and other experimental approaches on the other has played, and will likely continue to play, a vital role in Ras research.

GENERAL SIGNIFICANCE

Atomistic insights emerging from detailed simulations of Ras in solution and in bilayers may be the key to unlock the secret that to date prevented development of selective anti-Ras inhibitors for cancer therapy.

A small-molecule induces apoptosis and suppresses metastasis in pancreatic cancer cells

Pancreatic cancer is one of the most malignant tumor diseases with the characters of aggressive growth and metastasis. With the inefficiency of the current therapeutics, new potential targets and new therapeutic agents for healing of pancreatic cancer are critically needed. We have previously found a small molecule, named 4-tert-butyl-2-[(cyclohexylamino) methyl]-6-methylphenol (TBMMP, NSC number: 48160), which can freeze the intermediate of Ras-GTP hydrolysis in the open non-signaling conformation with high affinity and high specificity in silico. In this work, we studied the effect and mechanism of TBMMP on two pancreatic cancer cell lines, CFPAC-1 and BxPC-3. The results showed that TBMMP could restrain the growth of the pancreatic cancer cells with IC(50) value 84.3 μM for CPFAC-1 and 94.5 μM for BxPC-3, respectively. Additionally, TBMMP increased cytochrome c release, reduced mitochondrial membrane potential, activated caspase-3, -9, elevated ROS and increased expression of the Bax in the pancreatic cancer cell lines. The results indicated that TBMMP induced the apoptosis of pancreatic cancer cells through the mitochondrial pathway. Further, we also found that TBMMP could suppress the metastasis of both pancreatic cancer cells in vitro. Taken together, we proposed that TBMMP might be a therapeutic potential lead for treating patients with pancreatic cancer.

Ras inhibition via direct Ras binding--is there a path forward?

Three decades after identification of the Ras oncogene, no effective treatments for Ras mutant tumors are available despite intensive drug discovery efforts. Here we critically review the attempts to inhibit Ras function via direct binding of small molecules at the Ras surface with the aim to disrupt its interaction with other proteins. Multiple binders at different binding sites have been discovered, and recent efforts afforded crystal structures of Ras-binder complexes. Albeit with low affinities, many of the binders were shown to impart inhibitory activities, and inhibition of nucleotide exchange as a consequence of disrupting the Ras-SOS interaction has been the most commonly identified mode of action. We see two key challenges in the development of these early starting points: Enhancing binding affinities and achieving selectivity, both against other GTPases and for mutant Ras over the wildtype form. In light of the large unmet medical need, we encourage the continued search for functionally active Ras binders, and we believe that integrated use of biophysical and biochemical tools will provide the highest chances for success. Given the failures experienced in the past and the significant hurdles ahead, we propose that this challenge be tackled through alliances between industry and academia.

Electrostatic effects of mutations of Ras glutamine 61 measured using vibrational spectroscopy of a thiocyanate probe

Mutations of human oncoprotein p21(Ras) (hereafter Ras) at glutamine 61 are known to slow the rate of guanosine triphosphate (GTP) hydrolysis and transform healthy cells into malignant cells. It has been hypothesized that this glutamine plays a role in the intrinsic mechanism of GTP hydrolysis by interacting with an active site water molecule that electrostatically stabilizes the formation of the charged transition state at the γ-phosphate during hydrolysis. We have tested the interactions between amino acids at this position and water by measuring changes in the electrostatic field experienced by a nitrile probe positioned near Ras Q61 using vibrational Stark effect (VSE) spectroscopy. We mutated this glutamine to every amino acid except cysteine and proline and then incubated these mutants with a Ral guanine nucleotide dissociation stimulator (Ral) containing the I18C mutation that was chemically labeled with a thiocyanate vibrational spectroscopic probe. The formation of the docked Ras Q61X-labeled Ral complex was confirmed by measurement of the dissociation constant of the interaction. We measured the absorption energy of this nitrile to determine any differences in electrostatic environment in the immediate vicinity of the thiocyanate probe between wild type and mutants of Ras. For each Ras Q61X mutant, we correlate the change in electrostatic field at position 61 with the solvent accessible surface area of polar components of the mutant side chain determined from a Boltzmann-weighted ensemble of structures, as well as the residue's hydration potential. These results support the hypothesis that the role of Ras Q61 is to stabilize water in or near the active site during GTP hydrolysis. The substantial effect that nonpolar side chains of Ras Q61X have on the absorption energy of the thiocyanate must be investigated with further experiments.

The role of conserved waters in conformational transitions of Q61H K-ras

To investigate the stability and functional role of long-residence water molecules in the Q61H variant of the signaling protein K-ras, we analyzed all available Ras crystal structures and conformers derived from a series of independent explicit solvent molecular dynamics (MD) simulations totaling 1.76 µs. We show that the protein samples a different region of phase space in the presence and absence of several crystallographically conserved and buried water molecules. The dynamics of these waters is coupled with the local as well as the global motions of the protein, in contrast to less buried waters whose exchange with bulk is only loosely coupled with the motion of loops in their vicinity. Aided by two novel reaction coordinates involving the distance (d) between the C_α atoms of G60 at switch 2 and G10 at the P-loop and the N-C_α-C-O dihedral (ξ) of G60, we further show that three water molecules located in lobe1, at the interface between the lobes and at lobe2, are involved in the relative motion of residues at the two lobes of Q61H K-ras. Moreover, a d/ξ plot classifies the available Ras x-ray structures and MD-derived K-ras conformers into active GTP-, intermediate GTP-, inactive GDP-bound, and nucleotide-free conformational states. The population of these states and the transition between them is modulated by water-mediated correlated motions involving the functionally critical switch 2, P-loop and helix 3. These results suggest that water molecules act as allosteric ligands to induce a population shift among distinct switch 2 conformations that differ in effector recognition.

K-ras belongs to the Ras family of G-proteins that regulate cell proliferation and development. To execute its function, K-ras adopts different conformational states when it is active and inactive. In addition to these two states, it samples many transient intermediate conformations as it makes the transition from one state to the other. Mutations that affect the population of these states can cause cancer or developmental disorder. Using simulation approaches, here we show that a number of water molecules buried within the structure of an oncogenic K-ras protein modulate the distribution of its conformational states. Moreover, a detailed analysis based on two novel structural parameters revealed the existence of long-range water-mediated interactions that facilitate a dynamic coupling between the two lobes of the protein. These findings pave the way for a dynamics-guided strategy to inhibit abnormal Ras signaling.

Fluorescence detection of GDP in real time with the reagentless biosensor rhodamine-ParM

The development of novel fluorescence methods for the detection of key biomolecules is of great interest, both in basic research and in drug discovery. Particularly relevant and widespread molecules in cells are ADP and GDP, which are the products of a large number of cellular reactions, including reactions catalysed by nucleoside triphosphatases and kinases. Previously, biosensors for ADP were developed in this laboratory, based on fluorophore adducts with the bacterial actin homologue ParM. It is shown in the present study that one of these biosensors, tetramethylrhodamine–ParM, can also monitor GDP. The biosensor can be used to measure micromolar concentrations of GDP on the background of millimolar concentrations of GTP. The fluorescence response of the biosensor is fast, the response time being <0.2 s. Thus the biosensor allows real-time measurements of GTPase and GTP-dependent kinase reactions. Applications of the GDP biosensor are exemplified with two different GTPases, measuring the rates of GTP hydrolysis and nucleotide exchange.

Interconversion of two GDP-bound conformations and their selection in an Arf-family small G protein

ADP-ribosylation factor (Arf) and other Arf-family small G proteins participate in many cellular functions via their characteristic GTP/GDP conformational cycles, during which a nucleotide(∗)Mg(2+)-binding site communicates with a remote N-terminal helix. However, the conformational interplay between the nucleotides, the helix, the protein core, and Mg(2+) has not been fully delineated. Herein, we report a study of the dynamics of an Arf-family protein, Arl8, under various conditions by means of NMR relaxation spectroscopy. The data indicated that, when GDP is bound, the protein core, which does not include the N-terminal helix, reversibly transition between an Arf-family GDP form and another conformation that resembles the Arf-family GTP form. Additionally, we found that the N-terminal helix and Mg(2+), respectively, stabilize the aforementioned former and latter conformations in a population-shift manner. Given the dynamics of the conformational changes, we can describe the Arl8 GTP/GDP cycle in terms of an energy diagram.

Germline KRAS mutations cause aberrant biochemical and physical properties leading to developmental disorders

The KRAS gene is the most common locus for somatic gain-of-function mutations in human cancer. Germline KRAS mutations were shown recently to be associated with developmental disorders, including Noonan syndrome (NS), cardio-facio-cutaneous syndrome (CFCS), and Costello syndrome (CS). The molecular basis of this broad phenotypic variability has in part remained elusive so far. Here, we comprehensively analyzed the biochemical and structural features of ten germline KRAS mutations using physical and cellular biochemistry. According to their distinct biochemical and structural alterations, the mutants can be grouped into five distinct classes, four of which markedly differ from RAS oncoproteins. Investigated functional alterations comprise the enhancement of intrinsic and guanine nucleotide exchange factor (GEF) catalyzed nucleotide exchange, which is alternatively accompanied by an impaired GTPase-activating protein (GAP) stimulated GTP hydrolysis, an overall loss of functional properties, and a deficiency in effector interaction. In conclusion, our data underscore the important role of RAS in the pathogenesis of the group of related disorders including NS, CFCS, and CS, and provide clues to the high phenotypic variability of patients with germline KRAS mutations.

Structure of the dominant negative S17N mutant of Ras

The use of the dominant negative mutant of Ras has been crucial in elucidating the cellular signaling of Ras in response to the activation of various membrane-bound receptors. Although several point mutants of Ras exhibit a dominant negative effect, the asparagine to serine mutation at position 17 (S17N) remains the most popular and the most effective at inhibiting the activation of endogenous Ras. It is now widely accepted that the dominant negative effect is due to the ability of the mutant to sequester upstream activators and its inability to activate downstream effectors. Here, we present the crystal structure of RasS17N in the GDP-bound form. In the three molecules that populate the asymmetric unit, the Mg(2+) ion that normally coordinates the beta-phosphate is absent because of steric hindrance from the Asn17 side chain. Instead, a Ca(2+) ion is coordinating the alpha-phosphate. Also absent from one molecule is electron density for Phe28, a conserved residue that normally stabilizes the nucleotide's guanine base. Except for Phe28, the nucleotide makes conserved interactions with Ras. Combined, the inability of Phe28 to stabilize the guanine base and the absence of a Mg(2+) ion to neutralize the negative charges on the phosphates explain the weaker affinity of GDP for Ras. Our data suggest that the absence of the Mg(2+) should also dramatically affect GTP binding to Ras and the proper positioning of Thr35 necessary for the activation of switch 1 and the binding to downstream effectors, a prerequisite for the triggering of signaling pathways.

Functional study on GTP hydrolysis by the GTP-binding protein from Sulfolobus solfataricus, a member of the HflX family

GTPase domains from members of the HflX protein family have their catalytic glutamine residue of the DxxGQ motif substituted by phenylalanine, while they are still able to hydrolyse GTP. This appears to challenge the traditional view of GTP hydrolysis mechanism of Ras-like GTPases. SsGBP from the hyperthermophilic archaeon Sulfolobus solfataricus provided the first crystal structure of the HflX family. Here, we report structure-based mutagenesis analyses on SsGBP. Six-point mutations were individually introduced in the Ras-like GTPase domain including regions of P-loop, switches I and II. Intrinsic GTPase activities and thermal stabilities of these variants together with the wild-type full-length SsGBP and its isolated GTPase domain were analysed. Both functional and structural analyses of G235P and G235S mutants, which showed total and partial loss of the GTP hydrolyzing activity, respectively, support our hypothesis that the role of aligning a nucleophilic water molecule by the Ras Gln60 residue is replaced by the backbone amide group of Gly235 in SsGBP. Together with functional studies of other mutants, we conclude that the classical view of GTP hydrolysis mechanism likely remains the same in the HflX family with a twist in the entity of the nucleophilic alignment.

Characterization of a Ras mutant with identical GDP- and GTP-bound structures

We previously characterized the G60A mutant of Ras and showed that the switch regions of the GTP-bound but not the GDP-bound form of this mutant adopt an "open conformation" similar to that seen in nucleotide-free Ras. Here, we mutate Lys147 of the conserved (145)SAK(147) motif in the G60A background and characterize the resulting double mutant (DM). We show that RasDM is the first structure of a Ras protein with identical GDP- and GTP-bound structures. Both structures adopt the open conformation of the active form of RasG60A. The increase in the accessible surface area of the nucleotide is consistent with a 4-fold increase in its dissociation rate. Stopped-flow experiments show no major difference in the two-step kinetics of association of GDP or GTP with the wild type, G60A, or RasDM. Addition of Sos fails to accelerate nucleotide exchange. Overexpression of the G60A or double mutant of Ras in COS-1 cells fails to activate Erk and shows a strong dominant negative effect. Our data suggest that flexibility at position 60 is required for proper Sos-catalyzed nucleotide exchange and that structural information is somehow shared among the switch regions and the different nucleotide binding motifs.

Real-time NMR study of three small GTPases reveals that fluorescent 2'(3')-O-(N-methylanthraniloyl)-tagged nucleotides alter hydrolysis and exchange kinetics

The Ras family of small GTPases control diverse signaling pathways through a conserved "switch" mechanism, which is turned on by binding of GTP and turned off by GTP hydrolysis to GDP. Full understanding of GTPase switch functions requires reliable, quantitative assays for nucleotide binding and hydrolysis. Fluorescently labeled guanine nucleotides, such as 2'(3')-O-(N-methylanthraniloyl) (mant)-substituted GTP and GDP analogs, have been widely used to investigate the molecular properties of small GTPases, including Ras and Rho. Using a recently developed NMR method, we show that the kinetics of nucleotide hydrolysis and exchange by three small GTPases, alone and in the presence of their cognate GTPase-activating proteins (GAPs) and guanine nucleotide exchange factors, are affected by the presence of the fluorescent mant moiety. Intrinsic hydrolysis of mantGTP by Ras homolog enriched in brain (Rheb) is approximately 10 times faster than that of GTP, whereas it is 3.4 times slower with RhoA. On the other hand, the mant tag inhibits TSC2GAP-catalyzed GTP hydrolysis by Rheb but promotes p120 RasGAP-catalyzed GTP hydrolysis by H-Ras. Guanine nucleotide exchange factor-catalyzed nucleotide exchange for both H-Ras and RhoA was inhibited by mant-substituted nucleotides, and the degree of inhibition depends highly on the GTPase and whether the assay measures association of mantGTP with, or dissociation of mantGDP from the GTPase. These results indicate that the mant moiety has significant and unpredictable effects on GTPase reaction kinetics and underscore the importance of validating its use in each assay.

CHARMM: the biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

Ras conformational switching: simulating nucleotide-dependent conformational transitions with accelerated molecular dynamics

Ras mediates signaling pathways controlling cell proliferation and development by cycling between GTP- and GDP-bound active and inactive conformational states. Understanding the complete reaction path of this conformational change and its intermediary structures is critical to understanding Ras signaling. We characterize nucleotide-dependent conformational transition using multiple-barrier-crossing accelerated molecular dynamics (aMD) simulations. These transitions, achieved for the first time for wild-type Ras, are impossible to observe with classical molecular dynamics (cMD) simulations due to the large energetic barrier between end states. Mapping the reaction path onto a conformer plot describing the distribution of the crystallographic structures enabled identification of highly populated intermediate structures. These structures have unique switch orientations (residues 25-40 and 57-75) intermediate between GTP and GDP states, or distinct loop3 (46-49), loop7 (105-110), and alpha5 C-terminus (159-166) conformations distal from the nucleotide-binding site. In addition, these barrier-crossing trajectories predict novel nucleotide-dependent correlated motions, including correlations of alpha2 (residues 66-74) with alpha3-loop7 (93-110), loop2 (26-37) with loop10 (145-151), and loop3 (46-49) with alpha5 (152-167). The interconversion between newly identified Ras conformations revealed by this study advances our mechanistic understanding of Ras function. In addition, the pattern of correlated motions provides new evidence for a dynamic linkage between the nucleotide-binding site and the membrane interacting C-terminus critical for the signaling function of Ras. Furthermore, normal mode analysis indicates that the dominant collective motion that occurs during nucleotide-dependent conformational exchange, and captured in aMD (but absent in cMD) simulations, is a low-frequency motion intrinsic to the structure.

Mapping the nucleotide and isoform-dependent structural and dynamical features of Ras proteins

Ras GTPases are conformational switches controlling cell proliferation, differentiation, and development. Despite their prominent role in many forms of cancer, the mechanism of conformational transition between inactive GDP-bound and active GTP-bound states remains unclear. Here we describe a detailed analysis of available experimental structures and molecular dynamics simulations to quantitatively assess the structural and dynamical features of active and inactive states and their interconversion. We demonstrate that GTP-bound and nucleotide-free G12V H-ras sample a wide region of conformational space, and show that the inactive-to-active transition is a multiphase process defined by the relative rearrangement of the two switches and the orientation of Tyr32. We also modeled and simulated N- and K-ras proteins and found that K-ras is more flexible than N- and H-ras. We identified a number of isoform-specific, long-range side chain interactions that define unique pathways of communication between the nucleotide binding site and the C terminus.

Transformation efficiency of RasQ61 mutants linked to structural features of the switch regions in the presence of Raf

Transformation efficiencies of Ras mutants at residue 61 range over three orders of magnitude, but the in vitro GTPase activity decreases 10-fold for all mutants. We show that Raf impairs the GTPase activity of RasQ61L, suggesting that the Ras/Raf complex differentially modulates transformation. Our crystal structures show that, in transforming mutants, switch II takes part in a network of hydrophobic interactions burying the nucleotide and precatalytic water molecule. Our results suggest that Y32 and a water molecule bridging it to the gamma-phosphate in the wild-type structure play a role in GTP hydrolysis in lieu of the Arg finger in the absence of GAP. The bridging water molecule is absent in the transforming mutants, contributing to the burying of the nucleotide. We propose a mechanism for intrinsic hydrolysis in Raf-bound Ras and elucidate structural features in the Q61 mutants that correlate with their potency to transform cells.