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

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.

Solution structure and pressure response of thioredoxin-1 of Plasmodium falciparum

We present here the solution structures of the protein thioredoxin-1 from Plasmodium falciparum (PfTrx-1), in its reduced and oxidized forms. They were determined by high-resolution NMR spectroscopy at 293 K on uniformly 13C-, 15N-enriched, matched samples allowing to identification of even small structural differences. PfTrx-1 shows an α/β-fold with a mixed five-stranded β-sheet that is sandwiched between 4 helices in a β1 α1 β2 α2 β3 α3 β4 β5 α4 topology. The redox process of the CGPC motif leads to significant structural changes accompanied by larger chemical shift changes from residue Phe25 to Ile36, Thr70 to Thr74, and Leu88 to Asn91. By high-field high-pressure NMR spectroscopy, rare conformational states can be identified that potentially are functionally important and can be used for targeted drug development. We performed these experiments in the pressure range from 0.1 MPa to 200 MPa. The mean combined, random-coil corrected B1* values of reduced and oxidized thioredoxin are quite similar with -0.145 and -0.114 ppm GPa-1, respectively. The mean combined, random-coil corrected B2* values in the reduced and oxidized form are 0.179 and 0.119 ppm GPa-2, respectively. The mean ratios of the pressure coefficients B2/B1 are -0.484 and -0.831 GPa-1 in the reduced and oxidized form respectively. They differ at some points in the structure after the formation of the disulfide bond between C30 and C33. The thermodynamical description of the pressure dependence of chemical shifts requires the assumption of at least three coexisting conformational states of PfTrx-1. These three conformational states were identified in the reduced as well as in the oxidized form of the protein, therefore, they represent sub-states of the two main oxidation states of PfTrx-1.

The switch states of the GDP-bound HRAS affected by point mutations: a study from Gaussian accelerated molecular dynamics simulations and free energy landscapes

Point mutations play a vital role in the conformational transformation of HRAS. In this work, Gaussian accelerated molecular dynamics (GaMD) simulations followed by constructions of free energy landscapes (FELs) were adopted to explore the effect of mutations D33K, A59T and L120A on conformation states of the GDP-bound HRAS. The results from the post-processing analyses on GaMD trajectories suggest that mutations alter the flexibility and motion modes of the switch domains from HRAS. The analyses from FELs show that mutations induce more disordered states of the switch domains and affect interactions of GDP with HRAS, implying that mutations yield a vital effect on the binding of HRAS to effectors. The GDP-residue interaction network revealed by our current work indicates that salt bridges and hydrogen bonding interactions (HBIs) play key roles in the binding of GDP to HRAS. Furthermore, instability in the interactions of magnesium ions and GDP with the switch SI leads to the extreme disorder of the switch domains. This study is expected to provide the energetic basis and molecular mechanism for further understanding the function of HRAS.Communicated by Ramaswamy H. Sarma.

Mechanism of GTPase activation of a prokaryotic small Ras-like GTPase MglA by an asymmetrically interacting MglB dimer

Cell polarity oscillations in Myxococcus xanthus motility are driven by a prokaryotic small Ras-like GTPase, mutual gliding protein A (MglA), which switches from one cell pole to the other in response to extracellular signals. MglA dynamics is regulated by MglB, which functions both as a GTPase activating protein (GAP) and a guanine nucleotide exchange factor (GEF) for MglA. With an aim to dissect the asymmetric role of the two MglB protomers in the dual GAP and GEF activities, we generated a functional MglAB complex by coexpressing MglB with a linked construct of MglA and MglB. This strategy enabled us to generate mutations of individual MglB protomers (MglB1 or MglB2 linked to MglA) and delineate their role in GEF and GAP activities. We establish that the C-terminal helix of MglB1, but not MglB2, stimulates nucleotide exchange through a site away from the nucleotide-binding pocket, confirming an allosteric mechanism. Interaction between the N-terminal β-strand of MglB1 and β0 of MglA is essential for the optimal GEF activity of MglB. Specific residues of MglB2, which interact with Switch-I of MglA, partially contribute to its GAP activity. Thus, the role of the MglB2 protomer in the GAP activity of MglB is limited to restricting the conformation of MglA active site loops. The direct demonstration of the allosteric mechanism of GEF action provides us new insights into the regulation of small Ras-like GTPases, a feature potentially present in many uncharacterized GEFs.

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.

Conformational States of the GDP- and GTP-Bound HRAS Affected by A59E and K117R: An Exploration from Gaussian Accelerated Molecular Dynamics

The HRAS protein is considered a critical target for drug development in cancers. It is vital for effective drug development to understand the effects of mutations on the binding of GTP and GDP to HRAS. We conducted Gaussian accelerated molecular dynamics (GaMD) simulations and free energy landscape (FEL) calculations to investigate the impacts of two mutations (A59E and K117R) on GTP and GDP binding and the conformational states of the switch domain. Our findings demonstrate that these mutations not only modify the flexibility of the switch domains, but also affect the correlated motions of these domains. Furthermore, the mutations significantly disrupt the dynamic behavior of the switch domains, leading to a conformational change in HRAS. Additionally, these mutations significantly impact the switch domain's interactions, including their hydrogen bonding with ligands and electrostatic interactions with magnesium ions. Since the switch domains are crucial for the binding of HRAS to effectors, any alterations in their interactions or conformational states will undoubtedly disrupt the activity of HRAS. This research provides valuable information for the design of drugs targeting HRAS.

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).

Pushed to extremes: distinct effects of high temperature versus pressure on the structure of STEP

Protein function hinges on small shifts of three-dimensional structure. Elevating temperature or pressure may provide experimentally accessible insights into such shifts, but the effects of these distinct perturbations on protein structures have not been compared in atomic detail. To quantitatively explore these two axes, we report the first pair of structures at physiological temperature versus. high pressure for the same protein, STEP (PTPN5). We show that these perturbations have distinct and surprising effects on protein volume, patterns of ordered solvent, and local backbone and side-chain conformations. This includes interactions between key catalytic loops only at physiological temperature, and a distinct conformational ensemble for another active-site loop only at high pressure. Strikingly, in torsional space, physiological temperature shifts STEP toward previously reported active-like states, while high pressure shifts it toward a previously uncharted region. Altogether, our work indicates that temperature and pressure are complementary, powerful, fundamental macromolecular perturbations.

Exploring the structural dynamics of proteins by pressure perturbation using macromolecular crystallography

High pressure is a convenient thermodynamic parameter to probe the dynamics of proteins as it is intimately related to volume which is essential for protein function. To be biologically active, a protein fluctuates between different substates. Pressure perturbation can promote some hidden substates by modifying the population between them. High pressure macromolecular crystallography (HPMX) is a perfect tool to capture and to characterize such substates at a molecular level providing new insights on protein dynamics. The present chapter describes the use of the diamond anvil cell to perform HPMX experiments. It also provides tips on sample preparation and optimal data collection as well as on efficient analysis of the resulting high-pressure structures.

Reduced dynamic complexity allows structure elucidation of an excited state of KRASG13D

Localized dynamics of RAS, including regions distal to the nucleotide-binding site, is of high interest for elucidating the mechanisms by which RAS proteins interact with effectors and regulators and for designing inhibitors. Among several oncogenic mutants, methyl relaxation dispersion experiments reveal highly synchronized conformational dynamics in the active (GMPPNP-bound) KRAS^G13D, which suggests an exchange between two conformational states in solution. Methyl and ³¹P NMR spectra of active KRAS^G13D in solution confirm a two-state ensemble interconverting on the millisecond timescale, with a major P^γ atom peak corresponding to the dominant State 1 conformation and a secondary peak indicating an intermediate state different from the known State 2 conformation recognized by RAS effectors. High-resolution crystal structures of active KRAS^G13D and KRAS^G13D-RAF1 RBD complex provide snapshots of the State 1 and 2 conformations, respectively. We use residual dipolar couplings to solve and cross-validate the structure of the intermediate state of active KRAS^G13D, showing a conformation distinct from those of States 1 and 2 outside the known flexible switch regions. The dynamic coupling between the conformational exchange in the effector lobe and the breathing motion in the allosteric lobe is further validated by a secondary mutation in the allosteric lobe, which affects the conformational population equilibrium.

Pushed to extremes: distinct effects of high temperature vs. pressure on the structure of an atypical phosphatase

Protein function hinges on small shifts of three-dimensional structure. Elevating temperature or pressure may provide experimentally accessible insights into such shifts, but the effects of these distinct perturbations on protein structures have not been compared in atomic detail. To quantitatively explore these two axes, we report the first pair of structures at physiological temperature vs. high pressure for the same protein, STEP (PTPN5). We show that these perturbations have distinct and surprising effects on protein volume, patterns of ordered solvent, and local backbone and side-chain conformations. This includes novel interactions between key catalytic loops only at physiological temperature, and a distinct conformational ensemble for another active-site loop only at high pressure. Strikingly, in torsional space, physiological temperature shifts STEP toward previously reported active-like states, while high pressure shifts it toward a previously uncharted region. Together, our work argues that temperature and pressure are complementary, powerful, fundamental macromolecular perturbations.

Harnessing Pressure-Axis Experiments to Explore Volume Fluctuations, Conformational Substates, and Solvation of Biomolecular Systems

Intrinsic thermodynamic fluctuations within biomolecules are crucial for their function, and flexibility is one of the strategies that evolution has developed to adapt to extreme environments. In this regard, pressure perturbation is an important tool for mechanistically exploring the causes and effects of volume fluctuations in biomolecules and biomolecular assemblies, their role in biomolecular interactions and reactions, and how they are affected by the solvent properties. High hydrostatic pressure is also a key parameter in the context of deep-sea and subsurface biology and the study of the origin and physical limits of life. We discuss the role of pressure-axis experiments in revealing intrinsic structural fluctuations as well as high-energy conformational substates of proteins and other biomolecular systems that are important for their function and provide some illustrative examples. We show that the structural and dynamic information obtained from such pressure-axis studies improves our understanding of biomolecular function, disease, biological evolution, and adaptation.

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.

Structural communication between the GTPase Sec4p and its activator Sec2p: Determinants of GEF activity and early deformations to nucleotide release

Ras GTPases are molecular switches that cycle between OFF and ON states depending on the bound nucleotide (i.e. GDP-bound and GTP-bound, respectively).

The Rab GTPase, Sec4p, plays regulatory roles in multiple steps of intracellular vesicle trafficking. Nucleotide release is catalyzed by the Guanine Nucleotide Exchange Factor (GEF) Sec2p.

Here, the integration of structural information with molecular dynamics (MD) simulations addressed a number of questions concerning the intrinsic and stimulated dynamics of Sec2p and Sec4p as well as the chain of structural deformations leading to GEF-assisted activation of the Rab GTPase.

Sec2p holds an intrinsic ability to adopt the conformation found in the crystallographic complexes with Sec4p, thus suggesting that the latter selects and shifts the conformational equilibrium towards a pre-existing bound-like conformation of Sec2p.

The anchoring of Sec4p to a suitable conformation of Sec2p favors the Sec2p-assisted pulling on itself of the α1/switch 1 (SWI) loop and of SWI, which loose any contact with GDP. Those deformations of Sec4p would occur earlier. Formation of the final Sec2p-Sec4p hydrophobic interface, accomplishes later. Disruption of the nucleotide cage would cause firstly loss of interactions with the guanine ring and secondly loss of interactions with the phosphates.

The ease in sampling the energy landscape and adopting a bound-like conformation likely favors the catalyzing ability of GEFs for Ras GTPases.