The GAP arginine finger movement into the catalytic site of Ras increases the activation entropy

Allosteric modulation of NF1 GAP: Differential distributions of catalytically competent populations in loss-of-function and gain-of-function mutants

Neurofibromin (NF1), a Ras GTPase-activating protein (GAP), catalyzes Ras-mediated GTP hydrolysis and thereby negatively regulates the Ras/MAPK pathway. NF1 mutations can cause neurofibromatosis type 1 manifesting tumors, and neurodevelopmental disorders. Exactly how the missense mutations in the GAP-related domain of NF1 (NF1GRD) allosterically impact NF1 GAP to promote these distinct pathologies is unclear. Especially tantalizing is the question of how same-domain, same-residue NF1GRD variants exhibit distinct clinical phenotypes. Guided by clinical data, we take up this dilemma. We sampled the conformational ensembles of NF1GRD in complex with GTP-bound K-Ras4B by performing molecular dynamics simulations. Our results show that mutations in NF1GRD retain the active conformation of K-Ras4B but with biased propensities of the catalytically competent populations of K-Ras4B-NF1GRD complex. In agreement with clinical depiction and experimental tagging, compared to the wild type, NF1GRD E1356A and E1356V mutants effectively act through loss-of-function and gain-of-function mechanisms, leading to neurofibromatosis and developmental disorders, respectively. Allosteric modulation of NF1GRD GAP activity through biasing the conformational ensembles in the different states is further demonstrated by the diminished GAP activity by NF1GRD isoform 2, further manifesting propensities of conformational ensembles as powerful predictors of protein function. Taken together, our work identifies a NF1GRD hotspot that could allosterically tune GAP function, suggests targeting Ras oncogenic mutations by restoring NF1 catalytic activity, and offers a molecular mechanism for NF1 phenotypes determined by their distinct conformational propensities.

Mediating kinase activity in Ras-mutant cancer: potential for an individualised approach?

It is widely acknowledged that there is a considerable number of oncogenic mutations within the Ras superfamily of small GTPases which are the driving force behind a multitude of cancers. Ras proteins mediate a plethora of kinase pathways, including the MAPK, PI3K, and Ral pathways. Since Ras was considered undruggable until recently, pharmacological targeting of pathways downstream of Ras has been attempted to varying success, though drug resistance has often proven an issue. Nuances between kinase pathway activation in the presence of various Ras mutants are thought to contribute to the resistance, however, the reasoning behind activation of different pathways in different Ras mutational contexts is yet to be fully elucidated. Indeed, such disparities often depend on cancer type and disease progression. However, we are in a revolutionary age of Ras mutant targeted therapy, with direct-targeting KRAS-G12C inhibitors revolutionising the field and achieving FDA-approval in recent years. However, these are only beneficial in a subset of patients. Approximately 90% of Ras-mutant cancers are not KRAS-G12C mutant, and therefore raises the question as to whether other distinct amino acid substitutions within Ras may one day be targetable in a similar manner, and indeed whether better understanding of the downstream pathways these various mutants activate could further improve therapy. Here, we discuss the favouring of kinase pathways across an array of Ras-mutant oncogenic contexts and assess recent advances in pharmacological targeting of various Ras mutants. Ultimately, we will examine the utility of individualised pharmacological approaches to Ras-mediated cancer.

Differential roles of putative arginine fingers of AAA+ ATPases Rvb1 and Rvb2

The evolutionarily conserved AAA+ ATPases Rvb1 and Rvb2 proteins form a heteromeric complex (Rvb1/2) required for assembly or remodeling of macromolecular complexes in essential cellular processes ranging from chromatin remodeling to ribosome biogenesis. Rvb1 and Rvb2 have a high degree of sequence and structural similarity, and both contain the classical features of ATPases of their clade, including an N-terminal AAA+ subdomain with the Walker A motif, an insertion domain that typically interacts with various binding partners, and a C-terminal AAA+ subdomain containing a Walker B motif, the Sensor I and II motifs, and an arginine finger. In this study, we find that despite the high degree of structural similarity, Rvb1 and Rvb2 have distinct active sites that impact their activities and regulation within the Rvb1/2 complex. Using a combination of biochemical and genetic approaches, we show that replacing the homologous arginine fingers of Rvb1 and Rvb2 with different amino acids not only has distinct effects on the catalytic activity of the complex, but also impacts cell growth, and the Rvb1/2 interactions with binding partners. Using molecular dynamics simulations, we find that changes near the active site of Rvb1 and Rvb2 cause long-range effects on the protein dynamics in the insertion domain, suggesting a molecular basis for how enzymatic activity within the catalytic site of ATP hydrolysis can be relayed to other domains of the Rvb1/2 complex to modulate its function. Further, we show the impact that the arginine finger variants have on snoRNP biogenesis and validate the findings from molecular dynamics simulations using a targeted genetic screen. Together, our results reveal new aspects of the regulation of the Rvb1/2 complex by identifying a relay of long-range molecular communication from the ATPase active site of the complex to the binding site of cofactors. Most importantly, our findings suggest that despite high similarity and cooperation within the same protein complex, the two proteins have evolved with unique properties critical for the regulation and function of the Rvb1/2 complex.

Structural basis for the distinct roles of non-conserved Pro116 and conserved Tyr124 of BCH domain of yeast p50RhoGAP

p50RhoGAP is a key protein that interacts with and downregulates the small GTPase RhoA. p50RhoGAP is a multifunctional protein containing the BNIP-2 and Cdc42GAP Homology (BCH) domain that facilitates protein-protein interactions and lipid binding and the GAP domain that regulates active RhoA population. We recently solved the structure of the BCH domain from yeast p50RhoGAP (YBCH) and showed that it maintains the adjacent GAP domain in an auto-inhibited state through the β5 strand. Our previous WT YBCH structure shows that a unique kink at position 116 thought to be made by a proline residue between alpha helices α6 and α7 is essential for the formation of intertwined dimer from asymmetric monomers. Here we sought to establish the role and impact of this Pro116. However, the kink persists in the structure of P116A mutant YBCH domain, suggesting that the scaffold is not dictated by the proline residue at this position. We further identified Tyr124 (or Tyr188 in HBCH) as a conserved residue in the crucial β5 strand. Extending to the human ortholog, when substituted to acidic residues, Tyr188D or Tyr188E, we observed an increase in RhoA binding and self-dimerization, indicative of a loss of inhibition of the GAP domain by the BCH domain. These results point to distinct roles and impact of the non-conserved and conserved amino acid positions in regulating the structural and functional complexity of the BCH domain.

ATP binding and ATP hydrolysis in full-length MsbA monitored via time-resolved Fourier transform infrared spectroscopy

The essential Escherichia coli ATPase MsbA is a lipid flippase that serves as a prototype for multi drug resistant ABC transporters. Its physiological function is the transport of lipopolisaccharides to build up the outer membranes of gram negative bacteria. Although several structural and biochemical studies of MsbA have been conducted previously, a detailed picture of the dynamic processes that link ATP hydrolysis to allocrit transport remains elusive. We report here for the first time time-resolved Fourier transform infrared (FTIR) spectroscopic measurements of the ATP binding and ATP hydrolysis reaction of full-length MsbA and determined reaction rates at 288 K of k ₁ = 0.49 ± 0.28 s^-1 and k ₂ = 0.014 ± 0.003 s^-1, respectively. We further verified these rates with photocaged NPEcgAppNHp where only nucleotide binding was observable and the negative mutant MsbA-H537A that showed slow hydrolysis (k ₂ < 2 × 10^-4 s^-1). Besides single turnover kinetics, FTIR measurements also deliver IR signatures of all educts, products and the protein. ADP remains protein-bound after ATP hydrolysis. In addition, the spectral changes observed for the two variants MsbA-S378A and MsbA-S482A correlated with the loss of hydrogen bonding to the γ-phosphate of ATP. This study paves the way for FTIR-spectroscopic investigations of allocrite transport in full-length MsbA.

Lifetime of actin-dependent protein nanoclusters

Protein nanoclusters (PNCs) are dynamic collections of a few proteins that spatially organize in nanometer-length clusters. PNCs are one of the principal forms of spatial organization of membrane proteins, and they have been shown or hypothesized to be important in various cellular processes, including cell signaling. PNCs show remarkable diversity in size, shape, and lifetime. In particular, the lifetime of PNCs can vary over a wide range of timescales. The diversity in size and shape can be explained by the interaction of the clustering proteins with the actin cytoskeleton or the lipid membrane, but very little is known about the processes that determine the lifetime of the nanoclusters. In this paper, using mathematical modeling of the cluster dynamics, we model the biophysical processes that determine the lifetime of actin-dependent PNCs. In particular, we investigated the role of actin aster fragmentation, which had been suggested to be a key determinant of the PNC lifetime, and we found that it is important only for a small class of PNCs. A simple extension of our model allowed us to investigate the kinetics of protein-ligand interaction near PNCs. We found an anomalous increase in the lifetime of ligands near PNCs, which agrees remarkably well with experimental data on RAS-RAF kinetics. In particular, analysis of the RAS-RAF data through our model provides falsifiable predictions and novel hypotheses that will not only shed light on the role of RAS-RAF kinetics in various cancers, but also will be useful in studying membrane protein clustering in general.

BY-kinases: Protein tyrosine kinases like no other

BY-kinases (for bacterial tyrosine kinases) constitute a family of protein tyrosine kinases that are highly conserved in the bacterial kingdom and occur most commonly as essential components of multicomponent assemblies responsible for the biosynthesis, polymerization, and export of complex polysaccharides involved in biofilm or capsule formation. BY-kinase function has been attributed to a cyclic process involving formation of an oligomeric species, its disassembly into constituent monomers, and subsequent reassembly, depending on the overall phosphorylation level of a C-terminal cluster of tyrosine residues. However, the relationship of this process to the active/inactive states of the enzyme and the mechanism of its integration into the polysaccharide production machinery remain unclear. Here, we synthesize the substantial body of biochemical, cell-biological, structural, and computational data, acquired over the nearly 3 decades since the discovery of BY-kinases, to suggest means by which they fulfill their physiological function. We propose a mechanism involving temporal coordination of the assembly/disassembly process with the autokinase activity of the enzyme and its ability to be dephosphorylated by its counteracting phosphatase. We speculate that this temporal control enables BY-kinases to function as molecular timers that coordinate the diverse processes involved in the synthesis, polymerization, and export of complex sugar derivatives. We suggest that BY-kinases, which deploy distinctive catalytic domains resembling P-loop nucleoside triphosphatases, have uniquely adapted this ancient fold to drive functional processes through exquisite spatiotemporal control over protein-protein interactions and conformational changes. It is our hope that the hypotheses proposed here will facilitate future experiments targeting these unique protein kinases.

Millisecond molecular dynamics simulations of KRas-dimer formation and interfaces

Ras dimers have been proposed as building blocks for initiating the extracellular signal-regulated kinase (ERK)/mitogen-activated protein kinase (MAPK) cellular signaling pathway. To better examine the structure of possible dimer interfaces, the dynamics of Ras dimerization, and its potential signaling consequences, we performed molecular dynamics simulations totaling 1 ms of sampling, using an all-atom model of two full-length, farnesylated, guanosine triphosphate (GTP)-bound, wild-type KRas4b proteins diffusing on 29%POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine)-mixed POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membranes. Our simulations unveil an ensemble of thermodynamically weak KRas dimers spanning multiple conformations. The most stable conformations, having the largest interface areas, involve helix α2 and a hypervariable region (HVR). Among the dimer conformations, we found that the HVR of each KRas has frequent interactions with various parts of the dimer, thus potentially mediating the dimerization. Some dimer configurations have one KRas G-domain elevated above the lipid bilayer surface by residing on top of the other G-domain, thus likely contributing to the recruitment of cytosolic Raf kinases in the context of a stably formed multi-protein complex. We identified a variant of the α4-α5 KRas-dimer interface that is similar to the interfaces obtained with fluorescence resonance energy transfer (FRET) data of HRas on lipid bilayers. Interestingly, we found two arginine fingers, R68 and R149, that directly interact with the beta-phosphate of the GTP bound in KRas, in a manner similar to what is observed in a crystal structure of GAP-HRas complex, which can facilitate the GTP hydrolysis via the arginine finger of GTPase-activating protein (GAP).

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.

Common Patterns of Hydrolysis Initiation in P-loop Fold Nucleoside Triphosphatases

The P-loop fold nucleoside triphosphate (NTP) hydrolases (also known as Walker NTPases) function as ATPases, GTPases, and ATP synthases, are often of medical importance, and represent one of the largest and evolutionarily oldest families of enzymes. There is still no consensus on their catalytic mechanism. To clarify this, we performed the first comparative structural analysis of more than 3100 structures of P-loop NTPases that contain bound substrate Mg-NTPs or their analogues. We proceeded on the assumption that structural features common to these P-loop NTPases may be essential for catalysis. Our results are presented in two articles. Here, in the first, we consider the structural elements that stimulate hydrolysis. Upon interaction of P-loop NTPases with their cognate activating partners (RNA/DNA/protein domains), specific stimulatory moieties, usually Arg or Lys residues, are inserted into the catalytic site and initiate the cleavage of gamma phosphate. By analyzing a plethora of structures, we found that the only shared feature was the mechanistic interaction of stimulators with the oxygen atoms of gamma-phosphate group, capable of causing its rotation. One of the oxygen atoms of gamma phosphate coordinates the cofactor Mg ion. The rotation must pull this oxygen atom away from the Mg ion. This rearrangement should affect the properties of the other Mg ligands and may initiate hydrolysis according to the mechanism elaborated in the second article.

Genetically manipulating endogenous Kras levels and oncogenic mutations in vivo influences tissue patterning of murine tumorigenesis

Despite multiple possible oncogenic mutations in the proto-oncogene KRAS, unique subsets of these mutations are detected in different cancer types. As KRAS mutations occur early, if not being the initiating event, these mutational biases are ostensibly a product of how normal cells respond to the encoded oncoprotein. Oncogenic mutations can impact not only the level of active oncoprotein, but also engagement with proteins. To attempt to separate these two effects, we generated four novel Cre-inducible (LSL) Kras alleles in mice with the biochemically distinct G12D or Q61R mutations and encoded by native (nat) rare or common (com) codons to produce low or high protein levels. While there were similarities, each allele also induced a distinct transcriptional response shortly after activation in vivo. At one end of the spectrum, activating the KrasLSL-natG12D allele induced transcriptional hallmarks suggestive of an expansion of multipotent cells, while at the other end, activating the KrasLSL-comQ61R allele led to hallmarks of hyperproliferation and oncogenic stress. Evidence suggests that these changes may be a product of signaling differences due to increased protein expression as well as the specific mutation. To determine the impact of these distinct responses on RAS mutational patterning in vivo, all four alleles were globally activated, revealing that hematolymphopoietic lesions were permissive to the level of active oncoprotein, squamous tumors were permissive to the G12D mutant, while carcinomas were permissive to both these features. We suggest that different KRAS mutations impart unique signaling properties that are preferentially capable of inducing tumor initiation in a distinct cell-specific manner.

Mechanistic Insights into the Long-range Allosteric Regulation of KRAS Via Neurofibromatosis Type 1 (NF1) Scaffold Upon SPRED1 Loading

Allosteric regulation is the most direct and efficient way of regulating protein function, wherein proteins transmit the perturbations at one site to another distinct functional site. Deciphering the mechanism of allosteric regulation is of vital importance for the comprehension of both physiological and pathological events in vivo as well as the rational allosteric drug design. However, it remains challenging to elucidate dominant allosteric signal transduction pathways, especially for large and multi-component protein machineries where long-range allosteric regulation exits. One of the quintessential examples having long-range allosteric regulation is the ternary complex, SPRED1-RAS-neurofibromin type 1 (NF1, a RAS GTPase-activating protein), in which SPRED1 facilitates RAS-GTP hydrolysis by interacting with NF1 at a distal, allosteric site from the RAS binding site. To address the underlying mechanism, we performed extensive Gaussian accelerated molecular dynamics simulations and Markov state model analysis of KRAS-NF1 complex in the presence and absence of SPRED1. Our findings suggested that SPRED1 loading allosterically enhanced KRAS-NF1 binding, but hindered conformational transformation of the NF1 catalytic center for RAS hydrolysis. Moreover, we unveiled the possible allosteric pathways upon SPRED1 binding through difference contact network analysis. This study not only provided an in-depth mechanistic insight into the allosteric regulation of KRAS by SPRED1, but also shed light on the investigation of long-range allosteric regulation among complex macromolecular systems.

Identification of functional substates of KRas during GTP hydrolysis with enhanced sampling simulations

As the hub of major signaling pathways, Ras proteins are implicated in 19% of tumor-caused cancers due to perturbations in their conformational and/or catalytic properties. Despite numerous studies, the functions of the conformational substates for the most important isoform, KRas, remain elusive. In this work, we perform an extensive simulation analysis on the conformational landscape of KRas in its various chemical states during the GTP hydrolysis cycle: the reactant state KRasGTP·Mg²⁺, the intermediate state KRasGDP·Pi·Mg²⁺ and the product state KRasGDP·Mg²⁺. The results from enhanced sampling simulations reveal that State 1 of KRasGTP·Mg²⁺ has multiple stable substates in solution, one of which might account for interacting with GEFs. State 2 of KRasGTP·Mg²⁺ features two substates "Tyr32_in" and "Tyr32_out", which are poised to interact with effectors and GAPs, respectively. For the intermediate state KRasGDP·Pi·Mg²⁺, Gln61 and Pi are found to assume a broad set of conformations, which might account for the weak oncogenic effect of Gln61 mutations in KRas in contrast to the situation in HRas and NRas. Finally, the product state KRasGDP·Mg²⁺ has more than two stable substates in solution, pointing to a conformation-selection mechanism for complexation with GEFs. Based on these results, some specific inhibition strategies for targeting the binding sites of the high-energy substates of KRas during GTP hydrolysis are discussed.

An ATPase with a twist: A unique mechanism underlies the activity of the bacterial tyrosine kinase, Wzc

BY-kinases constitute a protein tyrosine kinase family that encodes unique catalytic domains that deviate from those of eukaryotic kinases resembling P-loop nucleotide triphosphatases (NTPases) instead. We have used computational and supporting biochemical approaches using the catalytic domain of the Escherichia coli BY-kinase, Wzc, to illustrate mechanistic divergences between BY-kinases and NTPases despite their deployment of similar catalytic motifs. In NTPases, the “arginine finger” drives the reactive conformation of ATP while also displacing its solvation shell, thereby making favorable enthalpic and entropic contributions toward βγ-bond cleavage. In BY-kinases, the reactive state of ATP is enabled by ATP·Mg²⁺-induced global conformational transitions coupled to the conformation of the Walker-A lysine. While the BY-kinase arginine finger does promote the desolvation of ATP, it does so indirectly by generating an ordered active site in combination with other structural elements. Bacteria, using these mechanistic variations, have thus repurposed an ancient fold to phosphorylate on tyrosine.

Divergent Mechanisms Activating RAS and Small GTPases Through Post-translational Modification

RAS is a founding member of the RAS superfamily of GTPases. These small 21 kDa proteins function as molecular switches to initialize signaling cascades involved in various cellular processes, including gene expression, cell growth, and differentiation. RAS is activated by GTP loading and deactivated upon GTP hydrolysis to GDP. Guanine nucleotide exchange factors (GEFs) and GTPase-activating proteins (GAPs) accelerate GTP loading and hydrolysis, respectively. These accessory proteins play a fundamental role in regulating activities of RAS superfamily small GTPase via a conserved guanine binding (G)-domain, which consists of five G motifs. The Switch regions lie within or proximal to the G2 and G3 motifs, and undergo dynamic conformational changes between the GDP-bound "OFF" state and GTP-bound "ON" state. They play an important role in the recognition of regulatory factors (GEFs and GAPs) and effectors. The G4 and G5 motifs are the focus of the present work and lie outside Switch regions. These motifs are responsible for the recognition of the guanine moiety in GTP and GDP, and contain residues that undergo post-translational modifications that underlie new mechanisms of RAS regulation. Post-translational modification within the G4 and G5 motifs activates RAS by populating the GTP-bound "ON" state, either through enhancement of intrinsic guanine nucleotide exchange or impairing GAP-mediated down-regulation. Here, we provide a comprehensive review of post-translational modifications in the RAS G4 and G5 motifs, and describe the role of these modifications in RAS activation as well as potential applications for cancer therapy.

Mechanism of Guanosine Triphosphate Hydrolysis by the Visual Proteins Arl3-RP2: Free Energy Reaction Profiles Computed with Ab Initio Type QM/MM Potentials

We report the results of calculations of the Gibbs energy profiles of the guanosine triphosphate (GTP) hydrolysis by the Arl3-RP2 protein complex using molecular dynamics (MD) simulations with ab initio type QM/MM potentials. The chemical reaction of GTP hydrolysis to guanosine diphosphate (GDP) and inorganic phosphate (Pi) is catalyzed by GTPases, the enzymes, which are responsible for signal transduction in live cells. A small GTPase Arl3, catalyzing the GTP → GDP reaction in complex with the activating protein RP2, constitute an essential part of the human vision cycle. To simulate the reaction mechanism, a model system is constructed by motifs of the crystal structure of the Arl3-RP2 complexed with a substrate analog. After selection of reaction coordinates, energy profiles for elementary steps along the reaction pathway GTP + H2O → GDP + Pi are computed using the umbrella sampling and umbrella integration procedures. QM/MM MD calculations are carried out, interfacing the molecular dynamics program NAMD and the quantum chemistry program TeraChem. Ab initio type QM(DFT)/MM potentials are computed with atom-centered basis sets 6-31G** and two hybrid functionals (PBE0-D3 and ωB97x-D3) of the density functional theory, describing a large QM subsystem. Results of these simulations of the reaction mechanism are compared to those obtained with QM/MM calculations on the potential energy surface using a similar description of the QM part. We find that both approaches, QM/MM and QM/MM MD, support the mechanism of GTP hydrolysis by GTPases, according to which the catalytic glutamine side chain (Gln71, in this system) actively participates in the reaction. Both approaches distinguish two parts of the reaction: the cleavage of the phosphorus-oxygen bond in GTP coupled with the formation of Pi, and the enzyme regeneration. Newly performed QM/MM MD simulations confirmed the profile predicted in the QM/MM minimum energy calculations, called here the pathway-I, and corrected its relief at the first elementary step from the enzyme-substrate complex. The QM/MM MD simulations also revealed another mechanism at the part of enzyme regeneration leading to pathway-II. Pathway-II is more consistent with the experimental kinetic data of the wild-type complex Arl3-RP2, whereas pathway-I explains the role of the mutation Glu138Gly in RP2 slowing down the hydrolysis rate.

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.

Targeting KRAS mutant cancers by preventing signaling transduction in the MAPK pathway

KRAS genes are the most commonly mutated oncogenes in cancer. Unfortunately, effective therapeutic strategies for targeting KRAS mutant cancers have proven to be difficult to obtain. A key reason for this setback is due to the lack of success direct KRAS mutant inhibitors have received. Researchers have turned their efforts away from targeting the KRAS nucleotide-binding site directly and towards targeting other areas of the MAPK signaling pathway to block KRAS function. Researchers found that inhibiting enzymes and protein-protein interactions involved in the MAPK signaling pathway inhibit the activation of KRAS mutant therefore can lead to a potential therapeutic for KRAS mutated cancers. Throughout the past two decades, various indirect inhibitors have been designed and tested. EGFR and MEK inhibitors have presented with less success; however, significant advances have been made when targeting the plasma membrane localization process and the allosteric site of KRAS mutant. Farnesyltransferase and allosteric inhibitors have both advanced to human clinical trials. This comprehensive review presents the most recent developments of direct and indirect KRAS mutant inhibitors. This review summarizes published data on the inhibitory and anti-cancer activity of compounds that target KRAS activation as well as highlights the most promising strategies for targeting KRAS mutant cancers.

Understanding the Biological Activities of Vitamin D in Type 1 Neurofibromatosis: New Insights into Disease Pathogenesis and Therapeutic Design

Vitamin D is a fat-soluble steroid hormone playing a pivotal role in calcium and phosphate homeostasis as well as in bone health. Vitamin D levels are not exclusively dependent on food intake. Indeed, the endogenous production-occurring in the skin and dependent on sun exposure-contributes to the majority amount of vitamin D present in the body. Since vitamin D receptors (VDRs) are ubiquitous and drive the expression of hundreds of genes, the interest in vitamin D has tremendously grown and its role in different diseases has been extensively studied. Several investigations indicated that vitamin D action extends far beyond bone health and calcium metabolism, showing broad effects on a variety of critical illnesses, including cancer, infections, cardiovascular and autoimmune diseases. Epidemiological studies indicated that low circulating vitamin D levels inversely correlate with cutaneous manifestations and bone abnormalities, clinical hallmarks of neurofibromatosis type 1 (NF1). NF1 is an autosomal dominant tumour predisposition syndrome causing significant pain and morbidity, for which limited treatment options are available. In this context, vitamin D or its analogues have been used to treat both skin and bone lesions in NF1 patients, alone or combined with other therapeutic agents. Here we provide an overview of vitamin D, its characteristic nutritional properties relevant for health benefits and its role in NF1 disorder. We focus on preclinical and clinical studies that demonstrated the clinical correlation between vitamin D status and NF1 disease, thus providing important insights into disease pathogenesis and new opportunities for targeted therapy.

Infrared Difference Spectroscopy of Proteins: From Bands to Bonds

Infrared difference spectroscopy probes vibrational changes of proteins upon their perturbation. Compared with other spectroscopic methods, it stands out by its sensitivity to the protonation state, H-bonding, and the conformation of different groups in proteins, including the peptide backbone, amino acid side chains, internal water molecules, or cofactors. In particular, the detection of protonation and H-bonding changes in a time-resolved manner, not easily obtained by other techniques, is one of the most successful applications of IR difference spectroscopy. The present review deals with the use of perturbations designed to specifically change the protein between two (or more) functionally relevant states, a strategy often referred to as reaction-induced IR difference spectroscopy. In the first half of this contribution, I review the technique of reaction-induced IR difference spectroscopy of proteins, with special emphasis given to the preparation of suitable samples and their characterization, strategies for the perturbation of proteins, and methodologies for time-resolved measurements (from nanoseconds to minutes). The second half of this contribution focuses on the spectral interpretation. It starts by reviewing how changes in H-bonding, medium polarity, and vibrational coupling affect vibrational frequencies, intensities, and bandwidths. It is followed by band assignments, a crucial aspect mostly performed with the help of isotopic labeling and site-directed mutagenesis, and complemented by integration and interpretation of the results in the context of the studied protein, an aspect increasingly supported by spectral calculations. Selected examples from the literature, predominately but not exclusively from retinal proteins, are used to illustrate the topics covered in this review.

Diversity of mechanisms in Ras-GAP catalysis of guanosine triphosphate hydrolysis revealed by molecular modeling

The mechanism of the deceptively simple reaction of guanosine triphosphate (GTP) hydrolysis catalyzed by the cellular protein Ras in complex with the activating protein GAP is an important issue because of the significance of this reaction in cancer research. We show that molecular modeling of GTP hydrolysis in the Ras-GAP active site reveals a diversity of mechanisms of the intrinsic chemical reaction depending on molecular groups at position 61 in Ras occupied by glutamine in the wild-type enzyme. First, a comparison of reaction energy profiles computed at the quantum mechanics/molecular mechanics (QM/MM) level shows that an assignment of the Gln61 side chain in the wild-type Ras either to QM or to MM parts leads to different scenarios corresponding to the glutamine-assisted or the substrate-assisted mechanisms. Second, replacement of Gln61 by the nitro-analog of glutamine (NGln) or by Glu, applied in experimental studies, results in two more scenarios featuring the so-called two-water and the concerted-type mechanisms. The glutamine-assisted mechanism in the wild-type Ras-GAP, in which the conserved Gln61 plays a decisive role, switching between the amide and imide tautomer forms, is consistent with the known experimental results of structural, kinetic and spectroscopy studies. The results emphasize the role of the Ras residue Gln61 in Ras-GAP catalysis and explain the retained catalytic activity of the Ras-GAP complex towards GTP hydrolysis in the Gln61NGln and Gln61Glu mutants of Ras.

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.

All Atom Motion Tree detects side chain-related motions and their coupling with domain motion in proteins

Structural changes of proteins are closely related with their molecular function. We previously developed a computational tool, Motion Tree (MT), to compare protein structures and describe structural changes using solely the Cα atoms. Here, we have extended MT to incorporate all heavy atoms to analyze side chain-related (SCR) motions. All Atom Motion Tree (AAMT) was applied to 76 proteins that exhibited a simple domain motion identified by MT. AAMT also detected 921 SCR motions. We examined the coupling of domain and SCR motions and classified the structural changes in terms of coupling. The statistical results indicated that it is common for coupled SCR motions to also couple with the domain motion. The classification correlates properties of domain motions and SCR motions. The AAMT results suggest that a large domain motion with a sizable domain boundary is accompanied by SCR motions composed of more than a single residue, which induces further couplings of SCR motions.

KRAS G13D sensitivity to neurofibromin-mediated GTP hydrolysis

KRAS mutations occur in ∼35% of colorectal cancers and promote tumor growth by constitutively activating the mitogen-activated protein kinase (MAPK) pathway. KRAS mutations at codons 12, 13, or 61 are thought to prevent GAP protein-stimulated GTP hydrolysis and render KRAS-mutated colorectal cancers unresponsive to epidermal growth factor receptor (EGFR) inhibitors. We report here that KRAS G13-mutated cancer cells are frequently comutated with NF1 GAP but NF1 is rarely mutated in cancers with KRAS codon 12 or 61 mutations. Neurofibromin protein (encoded by the NF1 gene) hydrolyzes GTP directly in complex with KRAS G13D, and KRAS G13D-mutated cells can respond to EGFR inhibitors in a neurofibromin-dependent manner. Structures of the wild type and G13D mutant of KRAS in complex with neurofibromin (RasGAP domain) provide the structural basis for neurofibromin-mediated GTP hydrolysis. These results reveal that KRAS G13D is responsive to neurofibromin-stimulated hydrolysis and suggest that a subset of KRAS G13-mutated colorectal cancers that are neurofibromin-competent may respond to EGFR therapies.

A GAP-GTPase-GDP-Pi Intermediate Crystal Structure Analyzed by DFT Shows GTP Hydrolysis Involves Serial Proton Transfers

Cell signaling by small G proteins uses an ON to OFF signal based on conformational changes following the hydrolysis of GTP to GDP and release of dihydrogen phosphate (Pi ). The catalytic mechanism of GTP hydrolysis by RhoA is strongly accelerated by a GAP protein and is now well defined, but timing of inorganic phosphate release and signal change remains unresolved. We have generated a quaternary complex for RhoA-GAP-GDP-Pi . Its 1.75 Å crystal structure shows geometry for ionic and hydrogen bond coordination of GDP and Pi in an intermediate state. It enables the selection of a QM core for DFT exploration of a 20 H-bonded network. This identifies serial locations of the two mobile protons from the original nucleophilic water molecule, showing how they move in three rational steps to form a stable quaternary complex. It also suggests how two additional proton transfer steps can facilitate Pi release.

GTP Hydrolysis Without an Active Site Base: A Unifying Mechanism for Ras and Related GTPases

GTP hydrolysis is a biologically crucial reaction, being involved in regulating almost all cellular processes. As a result, the enzymes that catalyze this reaction are among the most important drug targets. Despite their vital importance and decades of substantial research effort, the fundamental mechanism of enzyme-catalyzed GTP hydrolysis by GTPases remains highly controversial. Specifically, how do these regulatory proteins hydrolyze GTP without an obvious general base in the active site to activate the water molecule for nucleophilic attack? To answer this question, we perform empirical valence bond simulations of GTPase-catalyzed GTP hydrolysis, comparing solvent- and substrate-assisted pathways in three distinct GTPases, Ras, Rab, and the G_αi subunit of a heterotrimeric G-protein, both in the presence and in the absence of the corresponding GTPase activating proteins. Our results demonstrate that a general base is not needed in the active site, as the preferred mechanism for GTP hydrolysis is a conserved solvent-assisted pathway. This pathway involves the rate-limiting nucleophilic attack of a water molecule, leading to a short-lived intermediate that tautomerizes to form H₂PO₄^- and GDP as the final products. Our fundamental biochemical insight into the enzymatic regulation of GTP hydrolysis not only resolves a decades-old mechanistic controversy but also has high relevance for drug discovery efforts. That is, revisiting the role of oncogenic mutants with respect to our mechanistic findings would pave the way for a new starting point to discover drugs for (so far) "undruggable" GTPases like Ras.

Controlling the Revolving and Rotating Motion Direction of Asymmetric Hexameric Nanomotor by Arginine Finger and Channel Chirality

Nanomotors in nanotechnology are as important as engines in daily life. Many ATPases are nanoscale biomotors classified into three categories based on the motion mechanisms in transporting substrates: linear, rotating, and the recently discovered revolving motion. Most biomotors adopt a multisubunit ring-shaped structure that hydrolyzes ATP to generate force. How these biomotors control the motion direction and regulate the sequential action of their multiple subunits is intriguing. Many ATPases are hexameric with each monomer containing a conserved arginine finger. This review focuses on recent findings on how the arginine finger controls motion direction and coordinates adjacent subunit interactions in both revolving and rotating biomotors. Mechanisms of intersubunit interactions and sequential movements of individual subunits are evidenced by the asymmetrical appearance of one dimer and four monomers in high-resolution structural complexes. The arginine finger is situated at the interface of two subunits and extends into the ATP binding pocket of the downstream subunit. An arginine finger mutation results in deficiency in ATP binding/hydrolysis, substrate binding, and transport, highlighting the importance of the arginine finger in regulating energy transduction and motor function. Additionally, the roles of channel chirality and channel size are discussed as related to controlling one-way trafficking and differentiating the revolving and rotating mechanisms. Finally, the review concludes by discussing the conformational changes and entropy conversion triggered by ATP binding/hydrolysis, offering a view different from the traditional concept of ATP-mediated mechanochemical energy coupling. The elucidation of the motion mechanism and direction control in ATPases could facilitate nanomotor fabrication in nanotechnology.

Biochemical and kinetic properties of the complex Roco G-protein cycle

Roco proteins have come into focus after mutations in the gene coding for the human Roco protein Leucine-rich repeat kinase 2 (LRRK2) were discovered to be one of the most common genetic causes of late onset Parkinson's disease. Roco proteins are characterized by a Roc domain responsible for GTP binding and hydrolysis, followed by a COR dimerization device. The regulation and function of this RocCOR domain tandem is still not completely understood. To fully biochemically characterize Roco proteins, we performed a systematic survey of the kinetic properties of several Roco protein family members, including LRRK2. Together, our results show that Roco proteins have a unique G-protein cycle. Our results confirm that Roco proteins have a low nucleotide affinity in the micromolar range and thus do not strictly depend on G-nucleotide exchange factors. Measurement of multiple and single turnover reactions shows that neither Pi nor GDP release are rate-limiting, while this is the case for the GAP-mediated GTPase reaction of some small G-proteins like Ras and for most other high affinity Ras-like proteins, respectively. The KM values of the reactions are in the range of the physiological GTP concentration, suggesting that LRRK2 functioning might be regulated by the cellular GTP level.

Multi-target, ensemble-based virtual screening yields novel allosteric KRAS inhibitors at high success rate

RAS mutations account for >15% of all human tumors, and of these ~85% are due to mutations in a particular RAS gene: KRAS. Recent studies revealed that KRAS harbors four druggable allosteric sites. Here, we have (a) used molecular simulations to generate ensembles of wild type and four major oncogenic KRAS mutants (G12V, G12D, G13D, and Q61H); (b) characterized the druggability of each allosteric pocket in each protein; (c) conducted extensive ensemble-based virtual screening using pocket-tailored ligand libraries; (d) prioritized hits through hierarchical postdocking analysis; and (e) validated predicted hits with NMR. Of the 785 diverse potential hits identified by our in silico analysis, we tested 90 for their ability to bind KRAS using NMR and found that nine cause backbone amide chemical shift perturbations of residues near the functionally responsive switch loops, suggesting potential binding. We conducted detailed biophysical analyses on a novel indole-based compound to demonstrate the potential of our workflow to yield lead compounds. We believe the detailed information documented in this work regarding the druggability profile of each allosteric site and the chemical fingerprints of compounds that target them will serve as vital resources for future structure-based drug design efforts against KRAS, a high-value target for cancer therapy.

Ras-Specific GTPase-Activating Proteins-Structures, Mechanisms, and Interactions

Ras-specific GTPase-activating proteins (RasGAPs) down-regulate the biological activity of Ras proteins by accelerating their intrinsic rate of GTP hydrolysis, basically by a transition state stabilizing mechanism. Oncogenic Ras is commonly not sensitive to RasGAPs caused by interference of mutants with the electronic or steric requirements of the transition state, resulting in up-regulation of activated Ras in respective cells. RasGAPs are modular proteins containing a helical catalytic RasGAP module surrounded by smaller domains that are frequently involved in the subcellular localization or contributing to regulatory features of their host proteins. In this review, we summarize current knowledge about RasGAP structure, mechanism, regulation, and dual-substrate specificity and discuss in some detail neurofibromin, one of the most important negative Ras regulators in cellular growth control and neuronal function.

Human Miro Proteins Act as NTP Hydrolases through a Novel, Non-Canonical Catalytic Mechanism

Mitochondria are highly dynamic organelles that play a central role in multiple cellular processes, including energy metabolism, calcium homeostasis and apoptosis. Miro proteins (Miros) are “atypical” Ras superfamily GTPases that display unique domain architecture and subcellular localisation regulating mitochondrial transport, autophagy and calcium sensing. Here, we present systematic catalytic domain characterisation and structural analyses of human Miros. Despite lacking key conserved catalytic residues (equivalent to Ras Y32, T35, G60 and Q61), the Miro N-terminal GTPase domains display GTPase activity. Surprisingly, the C-terminal GTPase domains previously assumed to be “relic” domains were also active. Moreover, Miros show substrate promiscuity and function as NTPases. Molecular docking and structural analyses of Miros revealed unusual features in the Switch I and II regions, facilitating promiscuous substrate binding and suggesting the usage of a novel hydrolytic mechanism. The key substitution in position 13 in the Miros leads us to suggest the existence of an “internal arginine finger”, allowing an unusual catalytic mechanism that does not require GAP protein. Together, the data presented here indicate novel catalytic functions of human Miro atypical GTPases through altered catalytic mechanisms.

Amide-imide tautomerization in the glutamine side chain in enzymatic and photochemical reactions in proteins

Amide-imide tautomerization presents a pervasive class of chemical transformations in organic chemistry of natural compounds. In this Perspective, we describe two distinctively different protein systems, in which the amide-imide tautomerization in the glutamine side chain takes place in enzymatic or photochemical reactions. First, hydrolysis of guanosine triphosphate (GTP) catalyzed by the Ras-GAP protein complex suggests the occurrence of the imide tautomer of glutamine in reaction intermediates. Second, photoexcitation of flavin-binding protein domains (BLUFs) initiates a chain of reactions in the chromophore-binding pocket, including amide-imide tautomerization of glutamine. Mechanisms of these reactions at the atomic level have been revealed in quantum mechanics/molecular mechanics (QM/MM) simulations. To reinforce conclusions on the critical role of amide-imide tautomerization of glutamine in these reactions we describe results of new quantum chemistry and QM/MM calculations for relevant molecular model systems. We reexamine results of the recent IR spectroscopy studies of BLUF domains, which provide experimental evidences of Gln tautomerization in proteins. We also propose to validate the glutamine-assisted mechanism of enzymatic GTP hydrolysis by using IR spectroscopy in a proper range of wavenumbers.

Water Distribution within Wild-Type NRas Protein and Q61 Mutants during Unrestrained QM/MM Dynamics

Point mutations in p21^ras are associated with ∼30% of human tumors by disrupting its GTP hydrolysis cycle, which is critical to its molecular switch function in cellular signaling pathways. In this work, we investigate the impact of Gln 61 substitutions in the structure of the p21^N-ras active site and particularly focus on water reorganization around GTP, which appears to be crucial to evaluate favorable and unfavorable hydration sites for hydrolysis. The NRas-GTP complex is analyzed using a hybrid quantum mechanics/molecular mechanics approach, treating for the first time to our knowledge transient water molecules at the ab initio level and leading to results that account for the electrostatic coupling between the protein complex and the solvent. We show that for the wild-type protein, water molecules are found around the GTP γ-phosphate group, forming an arch extended from residues 12 to 35. Two density peaks are observed, supporting previous results that suggest the presence of two water molecules in the active site, one in the vicinity of residue 35 and a second one stabilized by hydrogen bonds formed with nitrogen backbone atoms of residues 12 and 60. The structural changes observed in NRas Gln 61 mutants result in the drastic delocalization of water molecules that we discuss. In mutants Q61H and Q61K, for which water distribution is overlocalized next to residue 60, the second density peak supports the hypothesis of a second water molecule. We also conclude that Gly 60 indirectly participates in GTP hydrolysis by correctly positioning transient water molecules in the protein complex and that Gln 61 has an indirect steric effect in stabilizing the preorganized catalytic site.

Specific Substates of Ras To Interact with GAPs and Effectors: Revealed by Theoretical Simulations and FTIR Experiments

The oncogenic Ras protein adopts various specific conformational states to execute its function in signal transduction. The large number of Ras structures obtained from X-ray and NMR experiments illustrates the diverse conformations that Ras adopts. It is difficult, however, to connect specific structural features with Ras functions. We report the free-energy landscape of Ras·GTP based on extensive explicit solvent simulations. The free-energy map clearly shows that the functional state 2 of Ras·GTP in fact has two distinct substates, denoted here as "Tyr32in" and "Tyr32out". Unbiased MD simulations show that the two substrates interconvert on the submicrosecond scale in solution, pointing to a novel mechanism for Ras·GTP to selectively interact with GAPs and effectors. This proposal is further supported by time-resolved FTIR experiments, which demonstrate that Tyr32 destabilizes the Ras·GAP complex and facilitates an efficient termination of Ras signaling.

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.

Importance of an Orchestrate Participation of each Individual Residue Present at a Catalytic Site

GTP hydrolysis is indispensable to keep a living cell healthy. Nature has evolved so many enzymes to enhance the slow GTP hydrolysis. Rab GTPases are evolved to regulate vesicle trafficking. GTPase activating proteins (GAPs) accelerates their intrinsic slow GTP hydrolysis in order to maintain the sustainability between cellular events. Any malfunction/interference in this hydrolysis disrupts normal cellular events and causes severe diseases. In this study, GTP hydrolysis mechanism of Rab33B catalyzed by TBC-domain GAP protein Gyp1p has been decoded using extensive ab initio QM/MM metadynamics simulations. An organized coupled movement of individual residues present at the catalytic site is found to be the key factor for this reaction. An unorganized coupled movement leads the hydrolysis through very high energy pathways. This also reveals that the chemical transformations occurring at a catalytic site are residue specific.

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.

Common mechanisms of catalysis in small and heterotrimeric GTPases and their respective GAPs

GTPases are central switches in cells. Their dysfunctions are involved in severe diseases. The small GTPase Ras regulates cell growth, differentiation and apoptosis by transmitting external signals to the nucleus. In one group of oncogenic mutations, the 'switch-off' reaction is inhibited, leading to persistent activation of the signaling pathway. The switch reaction is regulated by GTPase-activating proteins (GAPs), which catalyze GTP hydrolysis in Ras, and by guanine nucleotide exchange factors, which catalyze the exchange of GDP for GTP. Heterotrimeric G-proteins are activated by G-protein coupled receptors and are inactivated by GTP hydrolysis in the Gα subunit. Their GAPs are called regulators of G-protein signaling. In the same way that Ras serves as a prototype for small GTPases, Gαi1 is the most well-studied Gα subunit. By utilizing X-ray structural models, time-resolved infrared-difference spectroscopy, and biomolecular simulations, we elucidated the detailed molecular reaction mechanism of the GTP hydrolysis in Ras and Gαi1. In both proteins, the charge distribution of GTP is driven towards the transition state, and an arginine is precisely positioned to facilitate nucleophilic attack of water. In addition to these mechanistic details of GTP hydrolysis, Ras dimerization as an emerging factor in signal transduction is discussed in this review.

Mechanism of the intrinsic arginine finger in heterotrimeric G proteins

Heterotrimeric G proteins are crucial molecular switches that maintain a large number of physiological processes in cells. The signal is encoded into surface alterations of the Gα subunit that carries GTP in its active state and GDP in its inactive state. The ability of the Gα subunit to hydrolyze GTP is essential for signal termination. Regulator of G protein signaling (RGS) proteins accelerates this process. A key player in this catalyzed reaction is an arginine residue, Arg178 in Gα_i1, which is already an intrinsic part of the catalytic center in Gα in contrast to small GTPases, at which the corresponding GTPase-activating protein (GAP) provides the arginine "finger." We applied time-resolved FTIR spectroscopy in combination with isotopic labeling and site-directed mutagenesis to reveal the molecular mechanism, especially of the role of Arg178 in the intrinsic Gα_i1 mechanism and the RGS4-catalyzed mechanism. Complementary biomolecular simulations (molecular mechanics with molecular dynamics and coupled quantum mechanics/molecular mechanics) were performed. Our findings show that Arg178 is bound to γ-GTP for the intrinsic Gα_i1 mechanism and pushed toward a bidentate α-γ-GTP coordination for the Gα_i1·RGS4 mechanism. This movement induces a charge shift toward β-GTP, increases the planarity of γ-GTP, and thereby catalyzes the hydrolysis.

Pressure and Temperature Effects on the Activity and Structure of the Catalytic Domain of Human MT1-MMP

Membrane type 1-matrix metalloproteinase (MT1-MMP or MMP-14) is a zinc-transmembrane metalloprotease involved in the degradation of extracellular matrix and tumor invasion. While changes in solvation of MT1-MMP have been recently studied, little is known about the structural and energetic changes associated with MT1-MMP while interacting with substrates. Steady-state kinetic and thermodynamic data (including activation energies and activation volumes) were measured over a wide range of temperatures and pressures by means of a stopped-flow fluorescence technique. Complementary temperature- and pressure-dependent Fourier-transform infrared measurements provided corresponding structural information of the protein. MT1-MMP is stable and active over a wide range of temperatures (10-55 °C). A small conformational change was detected at 37 °C, which is responsible for the change in activity observed at the same temperature. Pressure decreases the enzymatic activity until complete inactivation occurs at 2 kbar. The inactivation is associated with changes in the rate-limiting step of the reaction caused by additional hydration of the active site upon compression and/or minor conformational changes in the active site region. Based on these data, an energy and volume diagram could be established for the various steps of the enzymatic reaction.

Theoretical vibrational spectroscopy of intermediates and the reaction mechanism of the guanosine triphosphate hydrolysis by the protein complex Ras-GAP

The structures and vibrational spectra of the reacting species upon guanosine triphosphate (GTP) hydrolysis to guanosine diphosphate and inorganic phosphate (Pi) trapped inside the protein complex Ras-GAP were analyzed following the results of QM/MM simulations. The frequencies of the phosphate vibrations referring to the reactants and to Pi were compared to those observed in the experimental FTIR studies. A good correlation between the theoretical and experimental vibrational data provides a strong support to the reaction mechanism of GTP hydrolysis by the Ras-GAP enzyme system revealed by the recent QM/MM modeling. Evolution of the vibrational bands associated with the inorganic phosphate Pi during the elementary stages of GTP hydrolysis is predicted.

Protein Regulation in Signal Transduction

SUMMARYCells must respond to a diverse, complex, and ever-changing mix of signals, using a fairly limited set of parts. Changes in protein level, protein localization, protein activity, and protein-protein interactions are critical aspects of signal transduction, allowing cells to respond highly specifically to a nearly limitless set of cues and also to vary the sensitivity, duration, and dynamics of the response. Signal-dependent changes in levels of gene expression and protein synthesis play an important role in regulation of protein levels, whereas posttranslational modifications of proteins regulate their degradation, localization, and functional interactions. Protein ubiquitylation, for example, can direct proteins to the proteasome for degradation or provide a signal that regulates their interactions and/or location within the cell. Similarly, protein phosphorylation by specific kinases is a key mechanism for augmenting protein activity and relaying signals to other proteins that possess domains that recognize the phosphorylated residues.

Rab32/38 and the xenophagic restriction of intracellular bacteria replication

Rab GTPases' subversion by intracellular pathogens during infection has been extensively documented. Recent findings have implicated a key intracellular bacterial restriction/containment function for Rab32/38 in Salmonella species in macrophages and Listeria monocytogenes in dendritic cells. Rab32/38 aids the phagolysosome maturation, and mediates a parallel xenophagy mechanism by engaging prohibitins.

Two-Photon Activation of p-Hydroxyphenacyl Phototriggers: Toward Spatially Controlled Release of Diethyl Phosphate and ATP

Two-photon activation of the p-hydroxyphenacyl (pHP) photoactivated protecting group is demonstrated for the first time using visible light at 550 nm from a pulsed laser. Broadband two-photon absorption measurements reveal a strong two-photon transition (>10 GM) near 4.5 eV that closely resembles the lowest-energy band at the same total excitation energy in the one-photon absorption spectrum of the pHP chromophore. The polarization dependence of the two-photon absorption band is consistent with excitation to the same S3 ((1)ππ*) excited state for both one- and two-photon activation. Monitoring the progress of the uncaging reaction under nonresonant excitation at 550 nm confirms a quadratic intensity dependence and that two-photon activation of the uncaging reaction is possible using visible light in the range 500-620 nm. Deprotonation of the pHP chromophore under mildly basic conditions shifts the absorption band to lower energy (3.8 eV) in both the one- and two-photon absorption spectra, suggesting that two-photon activation of the pHP chromophore may be possible using light in the range 550-720 nm. The results of these measurements open the possibility of spatially and temporally selective release of biologically active compounds from the pHP protecting group using visible light from a pulsed laser.

(19)F NMR and DFT Analysis Reveal Structural and Electronic Transition State Features for RhoA-Catalyzed GTP Hydrolysis

Molecular details for RhoA/GAP catalysis of the hydrolysis of GTP to GDP are poorly understood. We use (19)F NMR chemical shifts in the MgF3(-) transition state analogue (TSA) complex as a spectroscopic reporter to indicate electron distribution for the γ-PO3(-) oxygens in the corresponding TS, implying that oxygen coordinated to Mg has the greatest electron density. This was validated by QM calculations giving a picture of the electronic properties of the transition state (TS) for nucleophilic attack of water on the γ-PO3(-) group based on the structure of a RhoA/GAP-GDP-MgF3(-) TSA complex. The TS model displays a network of 20 hydrogen bonds, including the GAP Arg85' side chain, but neither phosphate torsional strain nor general base catalysis is evident. The nucleophilic water occupies a reactive location different from that in multiple ground state complexes, arising from reorientation of the Gln-63 carboxamide by Arg85' to preclude direct hydrogen bonding from water to the target γ-PO3(-) group.

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.