Molecular Mechanism for Conformational Dynamics of Ras·GTP Elucidated from In-Situ Structural Transition in Crystal

High Pressure Promotes Binding of the Allosteric Inhibitor Zn2+-Cyclen in Crystals of Activated H-Ras

In this work, we experimentally investigate the potency of high pressure to drive a protein toward an excited state where an inhibitor targeted for this state can bind. Ras proteins are small GTPases cycling between active GTP-bound and inactive GDP-bound states. Various states of GTP-bound Ras in active conformation coexist in solution, amongst them, state 2 which binds to effectors, and state 1, weakly populated at ambient conditions, which has a low affinity for effectors. Zn2+-cyclen is an allosteric inhibitor of Ras protein, designed to bind specifically to the state 1. In H-Ras(wt).Mg2+.GppNHp crystals soaked with Zn2+-cyclen, no binding could be observed, as expected in the state 2 conformation which is the dominant state at ambient pressure. Interestingly, Zn2+-cyclen binding is observed at 500 MPa pressure, close to the nucleotide, in Ras protein that is driven by pressure to a state 1 conformer. The unknown binding mode of Zn2+-cyclen to H-Ras can thus be fully characterized in atomic details. As a more general conjunction from our study, high pressure x-ray crystallography turns out to be a powerful method to induce transitions allowing drug binding in proteins that are in low-populated conformations at ambient conditions, enabling the design of specific inhibitors.

Assessing the mechanism of fast-cycling cancer-associated mutations of Rac1 small Rho GTPase

Rho-GTPases proteins function as molecular switches alternating from an active to an inactive state upon Guanosine triphosphate (GTP) binding and hydrolysis to Guanosine diphosphate (GDP). Among them, Rac subfamily regulates cell dynamics, being overexpressed in distinct cancer types. Notably, these proteins are object of frequent cancer-associated mutations at Pro29 (P29S, P29L, and P29Q). To assess the impact of these mutations on Rac1 structure and function, we performed extensive all-atom molecular dynamics simulations on wild-type (wt) and oncogenic isoforms of this protein in GDP- and GTP-bound states. Our results unprecedentedly elucidate that P29Q/S-induced structural and dynamical perturbations of Rac1 core domain weaken the binding of the catalytic site Mg2+ ion, and reduce the GDP residence time within protein, enhancing the GDP/GTP exchange rate and Rac1 activity. This broadens our knowledge of the role of cancer-associated mutations on small GTPases mechanism supplying valuable information for future drug discovery efforts targeting specific Rac1 isoforms.

GTP-Bound N-Ras Conformational States and Substates Are Modulated by Membrane and Point Mutation

Oncogenic Ras proteins are known to present multiple conformational states, as reported by the great variety of crystallographic structures. The GTP-bound states are grouped into two main states: the "inactive" state 1 and the "active" state 2. Recent reports on H-Ras have shown that state 2 exhibits two substates, directly related to the orientation of Tyr32: toward the GTP-bound pocket and outwards. In this paper, we show that N-Ras exhibits another substate of state 2, related to a third orientation of Tyr32, toward Ala18 and parallel to the GTP-bound pocket. We also show that this substate is highly sampled in the G12V mutation of N-Ras and barely present in its wild-type form, and that the G12V mutation prohibits the sampling of the GTPase-activating protein (GAP) binding substate, rendering this mutation oncogenic. Furthermore, using molecular dynamics simulations, we explore the importance of the membrane on N-Ras' conformational state dynamics and its strong influence on Ras protein stability. Moreover, the membrane has a significant influence on the conformational (sub)states sampling of Ras. This, in turn, is of crucial importance in the activation/deactivation cycle of Ras, due to the binding of guanine nucleotide exchange factor proteins (GEFs)/GTPase-activating proteins (GAPs).

Beyond structural bioinformatics for genomics with dynamics characterization of an expanded KRAS mutational landscape

Current capabilities in genomic sequencing outpace functional interpretations. Our previous work showed that 3D protein structure calculations enhance mechanistic understanding of genetic variation in sequenced tumors and patients with rare diseases. The KRAS GTPase is among the critical genetic factors driving cancer and germline conditions. Because KRAS-altered tumors frequently harbor one of three classic hotspot mutations, nearly all studies have focused on these mutations, leaving significant functional ambiguity across the broader KRAS genomic landscape observed in cancer and non-cancer diseases. Herein, we extend structural bioinformatics with molecular simulations to study an expanded landscape of 86 KRAS mutations. We identify multiple coordinated changes strongly associated with experimentally established KRAS biophysical and biochemical properties. The patterns we observe span hotspot and non-hotspot alterations, which can all dysregulate Switch regions, producing mutation-restricted conformations with different effector binding propensities. We experimentally measured mutation thermostability and identified shared and distinct patterns with simulations. Our results indicate mutation-specific conformations, which show potential for future research into how these alterations reverberate into different molecular and cellular functions. The data we present is not predictable using current genomic tools, demonstrating the added functional information derived from molecular simulations for interpreting human genetic variation.

Beyond Structural Bioinformatics for Genomics with Dynamics Characterization of an Expanded KRAS Mutational Landscape

Current capabilities in genomic sequencing outpace functional interpretations. Our previous work showed that 3D protein structure calculations enhance mechanistic understanding of genetic variation in sequenced tumors and patients with rare diseases. The KRAS GTPase is among the critical genetic factors driving cancer and germline conditions. Because KRAS-altered tumors frequently harbor one of three classic hotspot mutations, nearly all studies have focused on these mutations, leaving significant functional ambiguity across the broader KRAS genomic landscape observed in cancer and non-cancer diseases. Herein, we extend structural bioinformatics with molecular simulations to study an expanded landscape of 86 KRAS mutations. We identify multiple coordinated changes strongly associated with experimentally established KRAS biophysical and biochemical properties. The patterns we observe span hotspot and non-hotspot alterations, which can all dysregulate Switch regions, producing mutation-restricted conformations with different effector binding propensities. We experimentally measured mutation thermostability and identified shared and distinct patterns with simulations. Our results indicate mutation-specific conformations which show potential for future research into how these alterations reverberate into different molecular and cellular functions. The data we present is not predictable using current genomic tools, demonstrating the added functional information derived from molecular simulations for interpreting human genetic variation.

GAP positions catalytic H-Ras residue Q61 for GTP hydrolysis in molecular dynamics simulations, complicating chemical rescue of Ras deactivation

Functional interaction of Ras signaling proteins with upstream, negative regulatory GTPase activating proteins (GAPs) represents a crucial step in cellular decision making related to growth and survival. Key components of the catalytic transition state for Ras deactivation by GAP-accelerated hydrolysis of Ras-bound guanosine triphosphate (GTP) are thought to include an arginine residue from the GAP (the arginine finger), a glutamine residue from Ras (Q61), and a water molecule that is likely coordinated by Q61 to engage in nucleophilic attack on GTP. Here, we use in-vitro fluorescence experiments to show that 0.1-100 mM concentrations of free arginine, imidazole, and other small nitrogenous molecule fail to accelerate GTP hydrolysis, even in the presence of the catalytic domain of a mutant GAP lacking its arginine finger (R1276A NF1). This result is surprising given that imidazole can chemically rescue enzyme activity in arginine-to-alanine mutant protein tyrosine kinases (PTKs) that share many active site components with Ras/GAP complexes. Complementary all-atom molecular dynamics (MD) simulations reveal that an arginine finger GAP mutant still functions to enhance Ras Q61-GTP interaction, though less extensively than wild-type GAP. This increased Q61-GTP proximity may promote more frequent fluctuations into configurations that enable GTP hydrolysis as a component of the mechanism by which GAPs accelerate Ras deactivation in the face of arginine finger mutations. The failure of small molecule analogs of arginine to chemically rescue catalytic deactivation of Ras is consistent with the idea that the influence of the GAP goes beyond the simple provision of its arginine finger. However, the failure of chemical rescue in the presence of R1276A NF1 suggests that the GAPs arginine finger is either unsusceptible to rescue due to exquisite positioning or that it is involved in complex multivalent interactions. Therefore, in the context of oncogenic Ras proteins with mutations at codons 12 or 13 that inhibit arginine finger penetration toward GTP, drug-based chemical rescue of GTP hydrolysis may have bifunctional chemical/geometric requirements that are more difficult to satisfy than those that result from arginine-to-alanine mutations in other enzymes for which chemical rescue has been demonstrated.

Common Mechanism of Activated Catalysis in P-loop Fold Nucleoside Triphosphatases-United in Diversity

To clarify the obscure hydrolysis mechanism of ubiquitous P-loop-fold nucleoside triphosphatases (Walker NTPases), we analysed the structures of 3136 catalytic sites with bound Mg-NTP complexes or their analogues. Our results are presented in two articles; here, in the second of them, we elucidated whether the Walker A and Walker B sequence motifs-common to all P-loop NTPases-could be directly involved in catalysis. We found that the hydrogen bonds (H-bonds) between the strictly conserved, Mg-coordinating Ser/Thr of the Walker A motif ([Ser/Thr]WA) and aspartate of the Walker B motif (AspWB) are particularly short (even as short as 2.4 ångströms) in the structures with bound transition state (TS) analogues. Given that a short H-bond implies parity in the pKa values of the H-bond partners, we suggest that, in response to the interactions of a P-loop NTPase with its cognate activating partner, a proton relocates from [Ser/Thr]WA to AspWB. The resulting anionic [Ser/Thr]WA alkoxide withdraws a proton from the catalytic water molecule, and the nascent hydroxyl attacks the gamma phosphate of NTP. When the gamma-phosphate breaks away, the trapped proton at AspWB passes by the Grotthuss relay via [Ser/Thr]WA to beta-phosphate and compensates for its developing negative charge that is thought to be responsible for the activation barrier of hydrolysis.

Delineating The RAS Conformational Landscape

Mutations in RAS isoforms (KRAS, NRAS, and HRAS) are among the most frequent oncogenic alterations in many cancers, making these proteins high priority therapeutic targets. Effectively targeting RAS isoforms requires an exact understanding of their active, inactive, and druggable conformations. However, there is no structural catalog of RAS conformations to guide therapeutic targeting or examining the structural impact of RAS mutations. Here we present an expanded classification of RAS conformations based on analyses of the catalytic switch 1 (SW1) and switch 2 (SW2) loops. From 721 human KRAS, NRAS, and HRAS structures available in the Protein Data Bank (206 RAS-protein co-complexes, 190 inhibitor-bound, and 325 unbound, including 204 WT and 517 mutated structures), we created a broad conformational classification based on the spatial positions of Y32 in SW1 and Y71 in SW2. Clustering all well-modeled SW1 and SW2 loops using a density-based machine learning algorithm defined additional conformational subsets, some previously undescribed. Three SW1 conformations and nine SW2 conformations were identified, each associated with different nucleotide states (GTP-bound, nucleotide-free, and GDP-bound) and specific bound proteins or inhibitor sites. The GTP-bound SW1 conformation could be further subdivided based on the hydrogen bond type made between Y32 and the GTP γ-phosphate. Further analysis clarified the catalytic impact of G12D and G12V mutations and the inhibitor chemistries that bind to each druggable RAS conformation. Overall, this study has expanded our understanding of RAS structural biology, which could facilitate future RAS drug discovery.

TRPM8-Rap1A Interaction Sites as Critical Determinants for Adhesion and Migration of Prostate and Other Epithelial Cancer Cells

Emerging evidence indicates that the TRPM8 channel plays an important role in prostate cancer (PCa) progression, by impairing the motility of these cancer cells. Here, we reveal a novel facet of PCa motility control via direct protein-protein interaction (PPI) of the channel with the small GTPase Rap1A. The functional interaction of the two proteins was assessed by active Rap1 pull-down assays and live-cell imaging experiments. Molecular modeling analysis allowed the identification of four putative residues involved in TRPM8-Rap1A interaction. Point mutations of these sites impaired PPI as shown by GST-pull-down, co-immunoprecipitation, and PLA experiments and revealed their key functional role in the adhesion and migration of PC3 prostate cancer cells. More precisely, TRPM8 inhibits cell migration and adhesion by trapping Rap1A in its GDP-bound inactive form, thus preventing its activation at the plasma membrane. In particular, residues E207 and Y240 in the sequence of TRPM8 and Y32 in that of Rap1A are critical for the interaction between the two proteins not only in PC3 cells but also in cervical (HeLa) and breast (MCF-7) cancer cells. This study deepens our knowledge of the mechanism through which TRPM8 would exert a protective role in cancer progression and provides new insights into the possible use of TRPM8 as a new therapeutic target in cancer treatment.

Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.

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.

Elucidating microscopic events driven by GTP hydrolysis reaction in the Ras-GAP system with semi-reactive molecular dynamics simulations: the alternative role of a phosphate binding loop for mechanical energy storage

ATPase and GTPase have been widely found as chemical energy-mechanical work transducers, whereas the physicochemical mechanisms are not satisfactorily understood. We addressed the problem by examining John Ross' conjecture that repulsive Coulomb interaction between ADP/GDP and inorganic phosphate (P_i) does the mechanical work upon the system. We effectively simulated the consequence of a GTP hydrolysis reaction in a complex system of Rat sarcoma (Ras) and GTPase activation protein (GAP) in the framework of classical molecular dynamics by switching force field parameters between the reactant and product systems. We then observed a ca. 5 kcal mol^-1 increase of potential energy about the phosphate-binding loop (P-loop) in the Ras protein, indicating that the mechanical work generated via the GTP hydrolysis is converted into the local interaction energy and stored in the P-loop. Interestingly, this local energy storage in the P-loop depends on neither impulsive nor consecutive collisions of GDP and P_i with the P-loop. Instead, GTP-GDP conversion itself does work on the Ras system, elevating the potential energy. These observations encourage us to challenge a conjecture previously given by Ross. We assert that triphosphate nucleotide hydrolyses do mechanical work by producing emergent steric interaction accompanied by relaxation, namely, a shift of the biomolecular system to the non-equilibrium state on the reshaped potential energy landscape. Recalling the universality of the P-loop motif among GTPases and ATPases, the observations that we obtained through this study would progress the physicochemical understanding of the operating principles of GTP/ATP hydrolysis-driven biological nano-machines.

Structural basis for late maturation steps of the human mitoribosomal large subunit

Mitochondrial ribosomes (mitoribosomes) synthesize a critical set of proteins essential for oxidative phosphorylation. Therefore, mitoribosomal function is vital to the cellular energy supply. Mitoribosome biogenesis follows distinct molecular pathways that remain poorly understood. Here, we determine the cryo-EM structures of mitoribosomes isolated from human cell lines with either depleted or overexpressed mitoribosome assembly factor GTPBP5, allowing us to capture consecutive steps during mitoribosomal large subunit (mt-LSU) biogenesis. Our structures provide essential insights into the last steps of 16S rRNA folding, methylation and peptidyl transferase centre (PTC) completion, which require the coordinated action of nine assembly factors. We show that mammalian-specific MTERF4 contributes to the folding of 16S rRNA, allowing 16 S rRNA methylation by MRM2, while GTPBP5 and NSUN4 promote fine-tuning rRNA rearrangements leading to PTC formation. Moreover, our data reveal an unexpected involvement of the elongation factor mtEF-Tu in mt-LSU assembly, where mtEF-Tu interacts with GTPBP5, similar to its interaction with tRNA during translational elongation.

Oncogenic mutations Q61L and Q61H confer active form-like structural features to the inactive state (state 1) conformation of H-Ras protein

GTP-bound forms of Ras proteins (Ras•GTP) assume two interconverting conformations, "inactive" state 1 and "active" state 2. Our previous study on the crystal structure of the state 1 conformation of H-Ras in complex with guanosine 5'-(β, γ-imido)triphosphate (GppNHp) indicated that state 1 is stabilized by intramolecular hydrogen-bonding interactions formed by Gln61. Since Ras are constitutively activated by substitution mutations of Gln61, here we determine crystal structures of the state 1 conformation of H-Ras•GppNHp carrying representative mutations Q61L and Q61H to observe the effect of the mutations. The results show that these mutations alter the mode of hydrogen-bonding interactions of the residue 61 with Switch II residues and induce conformational destabilization of the neighboring regions. In particular, Q61L mutation results in acquirement of state 2-like structural features. Moreover, the mutations are likely to impair an intramolecular structural communication between Switch I and Switch II. Molecular dynamics simulations starting from these structures support the above observations. These findings may give a new insight into the molecular mechanism underlying the aberrant activation of the Gln61 mutants.

Normal Mode Analysis of KRas4B Reveals Partner Specific Dynamics

Ras GTPase interacts with its regulators and downstream effectors for its critical function in cellular signaling. Targeting the disrupted mechanisms in Ras-related human cancers requires understanding the distinct dynamics of these protein-protein interactions. We performed normal mode analysis (NMA) of KRas4B in wild-type or mutant monomeric and neurofibromin-1 (NF1), Son of Sevenless 1 (SOS1) or Raf-1 bound dimeric conformational states to reveal partner-specific dynamics of the protein. Gaussian network model (GNM) analysis showed that the known KRas4B lobes further partition into subdomains upon binding to its partners. Furthermore, KRas4B interactions with different partners suppress the flexibility in not only their binding sites but also distant residues in the allosteric lobe in a partner-specific way. The conformational changes can be driven by intrinsic residue fluctuations of the open state KRas4B-GDP, as we illustrated with anisotropic network model (ANM) analysis. The allosteric paths connecting the nucleotide binding residues to the allosteric site at α3-L7 portray differences in the inactive and active states. These findings help in understanding the partner-specific KRas4B dynamics, which could be utilized for therapeutic targeting.

Selective and noncovalent targeting of RAS mutants for inhibition and degradation

Activating mutants of RAS are commonly found in human cancers, but to date selective targeting of RAS in the clinic has been limited to KRAS(G12C) through covalent inhibitors. Here, we report a monobody, termed 12VC1, that recognizes the active state of both KRAS(G12V) and KRAS(G12C) up to 400-times more tightly than wild-type KRAS. The crystal structures reveal that 12VC1 recognizes the mutations through a shallow pocket, and 12VC1 competes against RAS-effector interaction. When expressed intracellularly, 12VC1 potently inhibits ERK activation and the proliferation of RAS-driven cancer cell lines in vitro and in mouse xenograft models. 12VC1 fused to VHL selectively degrades the KRAS mutants and provides more extended suppression of mutant RAS activity than inhibition by 12VC1 alone. These results demonstrate the feasibility of selective targeting and degradation of KRAS mutants in the active state with noncovalent reagents and provide a starting point for designing noncovalent therapeutics against oncogenic RAS mutants.

Mechanistic Differences of Activation of Rac1P29S and Rac1A159V

Rac1 is a small GTPase that plays key roles in actin reorganization, cell motility, and cell survival/growth as well as in various cancer types and neurodegenerative diseases. Similar to other Ras superfamily GTPases, Rac1 switches between active GTP-bound and inactive GDP-bound states. Switch I and II regions open and close during GDP/GTP exchange. P29S and A159V (paralogous to K-RasA146) mutations are the two most common somatic mutations of Rac1. Rac1P29S is a known hotspot for melanoma, whereas Rac1A159V most commonly occurs in head and neck cancer. To investigate how these substitutions induce the Rac1 dynamics, we used atomistic molecular dynamics simulations on the wild-type Rac1 and two mutant systems (P29S and A159V) in the GTP bound state, and on the wild-type Rac1 and P29S mutated system in the GDP bound state. Here, we show that P29S and A159V mutations activate Rac1 with different mechanisms. In Rac1P29S-GTP, the substitution increases the flexibility of Switch I based on RMSF and dihedral angle calculations and leads to an open conformation. We propose that the open Switch I conformation is one of the underlying reasons for rapid GDP/GTP exchange of Rac1P29S. On the other hand, in Rac1A159V-GTP, some of the contacts of the guanosine ring of GTP with Rac1 are temporarily lost, enabling the guanosine ring to move toward Switch I and subsequently close the switch. Rac1A159V-GTP adopts a Ras state 2 like conformation, where both switch regions are in closed conformation and Thr35 forms a hydrogen bond with the nucleotide.

Engineering subtilisin proteases that specifically degrade active RAS

We describe the design, kinetic properties, and structures of engineered subtilisin proteases that degrade the active form of RAS by cleaving a conserved sequence in switch 2. RAS is a signaling protein that, when mutated, drives a third of human cancers. To generate high specificity for the RAS target sequence, the active site was modified to be dependent on a cofactor (imidazole or nitrite) and protease sub-sites were engineered to create a linkage between substrate and cofactor binding. Selective proteolysis of active RAS arises from a 2-step process wherein sub-site interactions promote productive binding of the cofactor, enabling cleavage. Proteases engineered in this way specifically cleave active RAS in vitro, deplete the level of RAS in a bacterial reporter system, and also degrade RAS in human cell culture. Although these proteases target active RAS, the underlying design principles are fundamental and will be adaptable to many target proteins.

Structure of an inactive conformation of GTP-bound RhoA GTPase

By using ³¹P NMR, we present evidence that the Rho family GTPase RhoA, similar to Ras GTPases, exists in an equilibrium of conformations when bound to GTP. High-resolution crystal structures of RhoA bound to the GTP analog GMPPNP and to GDP show that they display a similar overall inactive conformation. In contrast to the previously reported crystal structures of GTP analog-bound forms of two RhoA dominantly active mutants (G14V and Q63L), GMPPNP-bound RhoA assumes an open conformation in the Switch I loop with a previously unseen interaction between the γ-phosphate and Pro36, instead of the canonical Thr37. Molecular dynamics simulations found that the oncogenic RhoA^G14V mutant displays a reduced flexibility in the Switch regions, consistent with a crystal structure of GDP-bound RhoA^G14V. Thus, GDP- and GTP-bound RhoA can present similar inactive conformations, and the molecular dynamics in the Switch regions are likely to have a role in RhoA activation.

Structural basis of ubiquitination mediated by protein splicing in early Eukarya

BACKGROUND

Inteins are intervening proteins, which are known to perform protein splicing. The reaction results in the production of an intein domain and an inteinless protein, which shows no trace of the insertion. BIL2 is part of the polyubiquitin locus of Tetrahymena thermophila (BUBL), where two bacterial-intein-like (BIL) domains lacking the C + 1 nucleophile, are flanked by two independent ubiquitin-like domains (ubl4/ubl5).

METHODS

We solved the X-ray structures of BIL2 in both the inactive and unprecedented, zinc-induced active, forms. Then, we characterized by mass spectrometry the BUBL splicing products in the absence and in the presence of T.thRas-GTPase. Finally, we investigated the effect of ubiquitination on T.thRas-GTPase by molecular dynamics simulations.

RESULTS

The structural analysis demonstrated that zinc-induced conformational change activates protein splicing. Moreover, mass spectrometry characterization of the splicing products shed light on the possible function of BIL2, which operates as a "single-ubiquitin-dispensing-platform", allowing the conjugation, via isopeptide bond formation (K(εNH₂)-C-ter), of ubl4 to either ubl5 or T.thRas-GTPase. Lastly, we demonstrated that T.thRas-GTPase ubiquitination occurs in proximity of the nucleotide binding pocket and stabilizes the protein active state.

CONCLUSIONS

We demonstrated that BIL2 is activated by zinc and that protein splicing induced by this intein does not take place through classical or aminolysis mechanisms but via formation of a covalent isopeptide bond, causing the ubiquitination of endogenous substrates such as T.thRas-GTPase.

GENERAL SIGNIFICANCE

In this "enzyme-free" ubiquitination mechanism the isopeptide formation, which canonically requires E1-E2-E3 enzymatic cascade and constitutes the alphabet of ubiquitin biology, is achieved in a single, concerted step without energy consumption.

Dynamically encoded reactivity of Ras enzymes: opening new frontiers for drug discovery

Decoding molecular flexibility in order to understand and predict biological processes-applying the principles of dynamic-structure-activity relationships (DSAR)-becomes a necessity when attempting to design selective and specific inhibitors of a protein that has overlapping interaction surfaces with its upstream and downstream partners along its signaling cascade. Ras proteins are molecular switches that meet this definition perfectly. The close-lying P-loop and the highly flexible switch I and switch II regions are the site of nucleotide-, assisting-, and effector-protein binding. Oncogenic mutations that also appear in this region do not cause easily characterized overall structural changes, due partly to the inherent conformational heterogeneity and pliability of these segments. In this review, we present an overview of the results obtained using approaches targeting Ras dynamics, such as nuclear magnetic resonance (NMR) measurements and experiment-based modeling calculations (mostly molecular dynamics (MD) simulations). These methodologies were successfully used to decipher the mutant- and isoform-specific nature of certain transient states, far-lying allosteric sites, and the internal interaction networks, as well as the interconnectivity of the catalytic and membrane-binding regions. This opens new therapeutic potential: the discovered interaction hotspots present hitherto not targeted, selective sites for drug design efforts in diverse locations of the protein matrix.

Structural impact of GTP binding on downstream KRAS signaling

Oncogenic RAS proteins, involved in ∼30% of human tumors, are molecular switches of various signal transduction pathways. Here we apply a new protocol for the NMR study of KRAS in its (inactive) GDP- and (activated) GTP-bound form, allowing a comprehensive analysis of the backbone dynamics of its WT-, G12C- and G12D variants. We found that Tyr32 shows opposite mobility with respect to the backbone of its surroundings: it is more flexible in the GDP-bound form while more rigid in GTP-complexes (especially in WT- and G12D-GTP). Using the G12C/Y32F double mutant, we showed that the presence of the hydroxyl group of Tyr32 has a marked effect on the G12C-KRAS-GTP system as well. Molecular dynamics simulations indicate that Tyr32 is linked to the γ-phosphate of GTP in the activated states – an arrangement shown, using QM/MM calculations, to support catalysis. Anchoring Tyr32 to the γ-phosphate contributes to the capture of the catalytic waters participating in the intrinsic hydrolysis of GTP and supports a simultaneous triple proton transfer step (catalytic water → assisting water → Tyr32 → O1G of the γ-phosphate) leading to straightforward product formation. The coupled flip of negatively charged residues of switch I toward the inside of the effector binding pocket potentiates ligand recognition, while positioning of Thr35 to enter the coordination sphere of the Mg²⁺ widens the pocket. Position 12 mutations do not disturb the capture of Tyr32 by the γ-phosphate, but (partially) displace Gln61, which opens up the catalytic pocket and destabilizes catalytic water molecules thus impairing intrinsic hydrolysis.

Nucleotide exchange to the physiological, activated, GTP-bound form of KRAS results in the anchoring of Tyr32 within the active site.

The Plasma Membrane as a Competitive Inhibitor and Positive Allosteric Modulator of KRas4B Signaling

Mutant Ras proteins are important drivers of human cancers, yet no approved drugs act directly on this difficult target. Over the last decade, the idea has emerged that oncogenic signaling can be diminished by molecules that drive Ras into orientations in which effector-binding interfaces are occluded by the cell membrane. To support this approach to drug discovery, we characterize the orientational preferences of membrane-bound K-Ras4B in 1.45-ms aggregate time of atomistic molecular dynamics simulations. Individual simulations probe active or inactive states of Ras on membranes with or without anionic lipids. We find that the membrane orientation of Ras is relatively insensitive to its bound guanine nucleotide and activation state but depends strongly on interactions with anionic phosphatidylserine lipids. These lipids slow Ras' translational and orientational diffusion and promote a discrete population in which small changes in orientation control Ras' competence to bind multiple regulator and effector proteins. Our results suggest that compound-directed conversion of constitutively active mutant Ras into functionally inactive forms may be accessible via subtle perturbations of Ras' orientational preferences at the membrane surface.

Exploiting RAS Nucleotide Cycling as a Strategy for Drugging RAS-Driven Cancers

Oncogenic mutations in RAS genes result in the elevation of cellular active RAS protein levels and increased signal propagation through downstream pathways that drive tumor cell proliferation and survival. These gain-of-function mutations drive over 30% of all human cancers, presenting promising therapeutic potential for RAS inhibitors. However, many have deemed RAS “undruggable” after nearly 40 years of failed drug discovery campaigns aimed at identifying a RAS inhibitor with clinical activity. Here we review RAS nucleotide cycling and the opportunities that RAS biochemistry presents for developing novel RAS inhibitory compounds. Additionally, compounds that have been identified to inhibit RAS by exploiting various aspects of RAS biology and biochemistry will be covered. Our current understanding of the biochemical properties of RAS, along with reports of direct-binding inhibitors, both provide insight on viable strategies for the discovery of novel clinical candidates with RAS inhibitory activity.

A temperature-controlled cold-gas humidifier and its application to protein crystals with the humid-air and glue-coating method

The room-temperature experiment has been revisited for macromolecular crystallography. Despite being limited by radiation damage, such experiments reveal structural differences depending on temperature, and it is expected that they will be able to probe structures that are physiologically alive. For such experiments, the humid-air and glue-coating (HAG) method for humidity-controlled experiments is proposed. The HAG method improves the stability of most crystals in capillary-free experiments and is applicable at both cryogenic and ambient temperatures. To expand the thermal versatility of the HAG method, a new humidifier and a protein-crystal-handling workbench have been developed. The devices provide temperatures down to 4°C and successfully maintain growth at that temperature of bovine cytochrome c oxidase crystals, which are highly sensitive to temperature variation. Hence, the humidifier and protein-crystal-handling workbench have proved useful for temperature-sensitive samples and will help reveal temperature-dependent variations in protein structures.

E. coli elongation factor Tu bound to a GTP analogue displays an open conformation equivalent to the GDP-bound form

According to the traditional view, GTPases act as molecular switches, which cycle between distinct ‘on’ and ‘off’ conformations bound to GTP and GDP, respectively. Translation elongation factor EF-Tu is a GTPase essential for prokaryotic protein synthesis. In its GTP-bound form, EF-Tu delivers aminoacylated tRNAs to the ribosome as a ternary complex. GTP hydrolysis is thought to cause the release of EF-Tu from aminoacyl-tRNA and the ribosome due to a dramatic conformational change following P_i release. Here, the crystal structure of Escherichia coli EF-Tu in complex with a non-hydrolysable GTP analogue (GDPNP) has been determined. Remarkably, the overall conformation of EF-Tu·GDPNP displays the classical, open GDP-bound conformation. This is in accordance with an emerging view that the identity of the bound guanine nucleotide is not ‘locking’ the GTPase in a fixed conformation. Using a single-molecule approach, the conformational dynamics of various ligand-bound forms of EF-Tu were probed in solution by fluorescence resonance energy transfer. The results suggest that EF-Tu, free in solution, may sample a wider set of conformations than the structurally well-defined GTP- and GDP-forms known from previous X-ray crystallographic studies. Only upon binding, as a ternary complex, to the mRNA-programmed ribosome, is the well-known, closed GTP-bound conformation, observed.

Quantitative Proteomics Combined with Affinity MS Revealed the Molecular Mechanism of Ginsenoside Antitumor Effects

Ginsenosides have previously been demonstrated to effectively inhibit cancer cell growth and survival in both animal models and cell lines. However, the specific ginsenoside component that is the active ingredient for cancer treatment through interaction with a target protein remains unknown. By an integrated quantitative proteomics approach via affinity mass spectrum (MS) technology, we deciphered the core structure of the ginsenoside active ingredient derived from crude extracts of ginsenosides and progressed toward identifying the target protein that mediates its anticancer activity. The Tandem Mass Tag (TMT) labeling quantitative proteomics technique acquired 55620 MS/MS spectra that identified 5499 proteins and 3045 modified proteins. Of these identified proteins, 224 differentially expressed proteins and modified proteins were significantly altered in nonsmall cell lung cancer cell lines. Bioinformatics tools for comprehensive analysis revealed that the Ras protein played a general regulatory role in many functional pathways and was probably the direct target protein of a compound in ginsenosides. Then, affinity MS screening based on the Ras protein identified 20(s)-protopanaxadiol, 20(s)-Ginsenoside Rh2, and 20(s)-Ginsenoside Rg3 had affinity with Ras protein under different conditions. In particular, 20(s)-protopanaxadiol, whose derivatives are the reported antitumor compounds 20(s)-Ginsenoside Rh2 and 20(s)-Ginsenoside Rg3 that have a higher affinity for Ras via a low KD of 1.22 μM and the mutation sites of G12 and G60, was demonstrated to play a core role in those interactions. Moreover, the molecular mechanism and bioactivity assessment results confirmed the identity of the chemical ligand that was directly acting on the GTP binding pocket of Ras and shown to be effective in cancer cell bioactivity profiles.

Methionine 170 is an Environmentally Sensitive Membrane Anchor in the Disordered HVR of K-Ras4B

Ras protein colocalization at the plasma membrane is implicated in the activation of signaling cascades that promote cell growth, survival, and motility. However, the mechanisms that underpin Ras self-association remain unclear. We use molecular dynamics simulations to show how basic and hydrophobic components of the disordered C-terminal membrane tether of K-Ras4B combine to regulate its membrane interactions. Specifically, anionic lipids attract lysine residues to the membrane surface, thereby splitting the peptide population into two states that exchange on the microsecond time scale. These states differ in the membrane insertion of a methionine residue, which is influenced by local membrane composition. As a result, these states may impose context-dependent biases on the disposition of Ras' signaling domain, with possible implications for the accessibility of its effector binding surfaces. We investigate Ras' ability to nanocluster by fly-casting for patches of anionic lipids and find that while anionic lipids promote the intermolecular association of K-Ras4B membrane tethers, at short range this appears to be a passive process in which anionic lipids electrostatically screen these cationic peptides to mitigate their natural repulsion. Together with the sub-microsecond stability of interpeptide contacts, this result suggests that experimentally observed K-Ras4B nanoclustering is not driven by direct intermolecular contact of its membrane tethers.

Allostery and dynamics in small G proteins

The Ras family of small guanine nucleotide-binding proteins behave as molecular switches: they are switched off and inactive when bound to GDP but can be activated by GTP binding in response to signal transduction pathways. Early structural analysis showed that two regions of the protein, which change conformation depending on the nucleotide present, mediate this switch. A large number of X-ray, NMR and simulation studies have shown that this is an over-simplification. The switch regions themselves are highly dynamic and can exist in distinct sub-states in the GTP-bound form that have different affinities for other proteins. Furthermore, regions outside the switches have been found to be sensitive to the nucleotide state of the protein, indicating that allosteric change is more widespread than previously thought. Taken together, the accrued knowledge about small G protein structures, allostery and dynamics will be essential for the design and testing of the next generation of inhibitors, both orthosteric and allosteric, as well as for understanding their mode of action.

Molecular Basis for Allosteric Inhibition of GTP-Bound H-Ras Protein by a Small-Molecule Compound Carrying a Naphthalene Ring

The ras oncogene products (H-Ras, K-Ras, and N-Ras) have been regarded as some of the most promising targets for anticancer drug discovery because their activating mutations are frequently found in human cancers. Nonetheless, molecular targeted therapy for them is currently unavailable. Here, we report the discovery of a small-molecule compound carrying a naphthalene ring, named KBFM123, which binds to the GTP-bound form of H-Ras. The solution structure of its complex with the guanosine 5'-(β,γ-imide) triphosphate-bound form of H-RasT35S (H-RasT35S·GppNHp) indicates that the naphthalene ring of KBFM123 interacts directly with a hydrophobic pocket located between switch I and switch II and allosterically inhibits the effector interaction by inducing conformational changes in switch I and its flanking region in strand β2, which are directly involved in recognition of the effector molecules, including c-Raf-1. In particular, Asp38 of H-Ras, a crucial residue for the interaction with c-Raf-1 via the formation of a salt bridge with Arg89 of the Ras-binding domain (RBD) of c-Raf-1, shows a drastic conformational change: its side chain orients toward the opposite direction. Consistent with these results, KBFM123 exhibits an activity to inhibit, albeit weakly, the association of H-RasG12V·GppNHp with the c-Raf-1 RBD. The binding of the naphthalene ring to the hydrophobic pocket of H-RasT35S·GppNHp is further supported by nuclear magnetic resonance analyses showing that two other naphthalene-containing compounds with distinct structures also exhibit similar binding properties with KBFM123. These results indicate that the naphthalene ring could become a promising scaffold for the development of Ras inhibitors.

Ras Binder Induces a Modified Switch-II Pocket in GTP and GDP States

Covalent inhibitors of K-Ras(G12C) have been reported that exclusively recognize the GDP state. Here, we utilize disulfide tethering of a non-natural cysteine (K-Ras(M72C)) to identify a new switch-II pocket (S-IIP) binding ligand (2C07) that engages the active GTP state. Co-crystal structures of 2C07 bound to H-Ras(M72C) reveal binding in a cryptic groove we term S-IIG. In the GppNHp state, 2C07 binding to a modified S-IIP pushes switch I away from the nucleotide, breaking the network of polar contacts essential for adopting the canonical GTP state. Biochemical studies show that 2C07 alters nucleotide preference and inhibits SOS binding and catalyzed nucleotide exchange. 2C07 was converted to irreversible covalent analogs, which target both nucleotide states, inhibit PI3K activation in vitro, and function as occupancy probes to detect reversible engagement in competition assays. Targeting both nucleotide states opens the possibility of inhibiting oncogenic mutants of Ras, which exist predominantly in the GTP state in cells.

Structural basis for intramolecular interaction of post-translationally modified H-Ras•GTP prepared by protein ligation

Ras undergoes post-translational modifications including farnesylation, proteolysis, and carboxymethylation at the C terminus, which are necessary for membrane recruitment and effector recognition. Full activation of c-Raf-1 requires cooperative interaction of the farnesylated C terminus and the activator region of Ras with its cysteine-rich domain (CRD). However, the molecular basis for this interaction remains unclear because of difficulties in preparing modified Ras in amounts sufficient for structural studies. Here, we use Sortase A-catalyzed protein ligation to prepare modified Ras in sufficient amounts for NMR and X-ray crystallographic analyses. The results show that the farnesylated C terminus establishes an intramolecular interaction with the catalytic domain and brings the farnesyl moiety to the proximity of the activator region, which may be responsible for their cooperative recognition of c-Raf-1-CRD.

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.

Nucleotide Dependent Switching in Rho GTPase: Conformational Heterogeneity and Competing Molecular Interactions

Ras superfamily of GTPases regulate myriad cellular processes through a conserved nucleotide (GTP/GDP) dependent switching mechanism. Unlike Ras family of GTPases, for the Rho GTPases, there is no clear evidence for the existence of “sub-states” such as state 1 & state 2 in the GTP bound form. To explore the nucleotide dependent conformational space of the Switch I loop and also to look for existence of state 1 like conformations in Rho GTPases, atomistic molecular dynamics and metadynamics simulations on RhoA were performed. These studies demonstrate that both the nucleotide-free state and the GDP bound “OFF” state have very similar conformations, whereas the GTP bound “ON” state has unique conformations with signatures of two intermediate states. The conformational free energy landscape for these systems suggests the presence of multiple intermediate states. Interestingly, the energetic penalty of exposing the non-polar residues in the GTP bound form is counter balanced by the favourable hydrogen bonded interactions between the γ-phosphate group of GTP with the highly conserved Tyr34 and Thr37 residues. These competing molecular interactions lead to a tuneable energy landscape of the Switch I conformation, which can undergo significant changes based on the local environment including changes upon binding to effectors.

Dynamic regulation of GDP binding to G proteins revealed by magnetic field-dependent NMR relaxation analyses

Heterotrimeric guanine-nucleotide-binding proteins (G proteins) serve as molecular switches in signalling pathways, by coupling the activation of cell surface receptors to intracellular responses. Mutations in the G protein α-subunit (Gα) that accelerate guanosine diphosphate (GDP) dissociation cause hyperactivation of the downstream effector proteins, leading to oncogenesis. However, the structural mechanism of the accelerated GDP dissociation has remained unclear. Here, we use magnetic field-dependent nuclear magnetic resonance relaxation analyses to investigate the structural and dynamic properties of GDP bound Gα on a microsecond timescale. We show that Gα rapidly exchanges between a ground-state conformation, which tightly binds to GDP and an excited conformation with reduced GDP affinity. The oncogenic D150N mutation accelerates GDP dissociation by shifting the equilibrium towards the excited conformation.

Heterotrimeric guanine-nucleotide-binding proteins (G proteins) act as molecular switches. Here the authors use NMR relaxation analyses, which reveal the dynamics of G protein alpha subunit binding to GDP on a microsecond timescale.