Regulation of GTPase function by autophosphorylation.

A unifying feature of the RAS superfamily is a conserved GTPase cycle by which these proteins transition between active and inactive states. We demonstrate that autophosphorylation of some GTPases is an intrinsic regulatory mechanism that reduces nucleotide hydrolysis and enhances nucleotide exchange, altering the on/off switch that forms the basis for their signaling functions. Using X-ray crystallography, nuclear magnetic resonance spectroscopy, binding assays, and molecular dynamics on autophosphorylated mutants of H-RAS and K-RAS, we show that phosphoryl transfer from GTP requires dynamic movement of the switch II region and that autophosphorylation promotes nucleotide exchange by opening the active site and extracting the stabilizing Mg2+. Finally, we demonstrate that autophosphorylated K-RAS exhibits altered effector interactions, including a reduced affinity for RAF proteins in mammalian cells. Thus, autophosphorylation leads to altered active site dynamics and effector interaction properties, creating a pool of GTPases that are functionally distinct from their non-phosphorylated counterparts.

Raf promotes dimerization of the Ras G-domain with increased allosteric connections.

Ras dimerization is critical for Raf activation. Here we show that the Ras binding domain of Raf (Raf-RBD) induces robust Ras dimerization at low surface densities on supported lipid bilayers and, to a lesser extent, in solution as observed by size exclusion chromatography and confirmed by SAXS. Community network analysis based on molecular dynamics simulations shows robust allosteric connections linking the two Raf-RBD D113 residues located in the Galectin scaffold protein binding site of each Raf-RBD molecule and 85 Å apart on opposite ends of the dimer complex. Our results suggest that Raf-RBD binding and Ras dimerization are concerted events that lead to a high-affinity signaling complex at the membrane that we propose is an essential unit in the macromolecular assembly of higher order Ras/Raf/Galectin complexes important for signaling through the Ras/Raf/MEK/ERK pathway.

Mechanisms of isoform-specific residue influence on GTP-bound HRas, KRas, and NRas.

HRas, KRas, and NRas are GTPases with a common set of effectors that control many cell-signaling pathways, including proliferation through Raf kinase. Their G-domains are nearly identical in sequence, with a few isoform-specific residues that have an effect on dynamics and biochemical properties. Here, we use accelerated molecular dynamics (aMD) simulations consistent with solution x-ray scattering experiments to elucidate mechanisms through which isoform-specific residues associated with each Ras isoform affects functionally important regions connected to the active site. HRas-specific residues cluster in loop 8 to stabilize the nucleotide-binding pocket, while NRas-specific residues on helix 3 directly affect the conformations of switch I and switch II. KRas, the most globally flexible of the isoforms, shows greatest fluctuations in the switch regions enhanced by a KRas-specific residue in loop 7 and a highly dynamic loop 8 region. The analysis of isoform-specific residue effects on Ras proteins is supported by NMR experiments and is consistent with previously published biochemical data.

DRoP: Automated detection of conserved solvent-binding sites on proteins.

Water and ligand binding play critical roles in the structure and function of proteins, yet their binding sites and significance are difficult to predict a priori. Multiple solvent crystal structures (MSCS) is a method where several X-ray crystal structures are solved, each in a unique solvent environment, with organic molecules that serve as probes of the protein surface for sites evolved to bind ligands, while the first hydration shell is essentially maintained. When superimposed, these structures contain a vast amount of information regarding hot spots of protein-protein or protein-ligand interactions, as well as conserved water-binding sites retained with the change in solvent properties. Optimized mining of this information requires reliable structural data and a consistent, objective analysis tool. Detection of related solvent positions (DRoP) was developed to automatically organize and rank the water or small organic molecule binding sites within a given set of structures. It is a flexible tool that can also be used in conserved water analysis given multiple structures of any protein independent of the MSCS method. The DRoP output is an HTML format list of the solvent sites ordered by conservation rank in its population within the set of structures, along with renumbered and recolored PDB files for visualization and facile analysis. Here, we present a previously unpublished set of MSCS structures of bovine pancreatic ribonuclease A (RNase A) and use it together with published structures to illustrate the capabilities of DRoP.

Water in Ras Superfamily Evolution.

The Ras GTPase superfamily of proteins coordinates a diverse set of cellular outcomes, including cell morphology, vesicle transport, and cell proliferation. Primary amino acid sequence analysis has identified Specificity determinant positions (SDPs) that drive diversified functions specific to the Ras, Rho, Rab, and Arf subfamilies (Rojas et al. 2012, J Cell Biol 196:189-201). The inclusion of water molecules in structural and functional adaptation is likely to be a major response to the selection pressures that drive evolution, yet hydration patterns are not included in phylogenetic analysis. This article shows that conserved crystallographic water molecules coevolved with SDP residues in the differentiation of proteins within the Ras superfamily of small GTPases. The patterns of water conservation between protein subfamilies parallel those of sequence-based evolutionary trees. Thus, hydration patterns have the potential to help elucidate functional significance in the evolution of amino acid residues observed in phylogenetic analysis of homologous proteins. © 2019 Wiley Periodicals, Inc.

Diets containing cunhã (Clitoria ternatea L.) hay and forage cactus (Opuntia sp.) meal on production and meat quality of Boer crossbred goat.

The aim of this study was to evaluate the effect of association of forage cactus meal (CM) and cunhã hay (CH) on carcass traits, meat quality, and fatty acid (FA) profile of goats. Twenty-four goats were finished on diets composed with the association of CM and CH (83% CH + 17% CM, 67% CH + 33% CM, or 50% CH + 50% CM), corresponding to 60% of the total diet; the remainder was composed of 10% elephant grass and 30% concentrate. The control treatment consisted of 70% elephant grass and 30% concentrate. Carcass yield and commercial cut weight were higher in the animals fed CH and CM at any concentration (P < 0.05). Animals fed with intermediate and high CM content resulted in meat with better proportions of unsaturated FA (P = 0.0281), desirable FA (P = 0.0024), omega 6/omega 3 ratio (P = 0.0340), and hypo-/hypercholesterolemic index (P = 0.0140). The inclusion of CH and CM provided carcass yield and FA profile more favorable for consumer market and human health.

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.

Neutron Crystal Structure of RAS GTPase Puts in Question the Protonation State of the GTP γ-Phosphate.

RAS GTPase is a prototype for nucleotide-binding proteins that function by cycling between GTP and GDP, with hydrogen atoms playing an important role in the GTP hydrolysis mechanism. It is one of the most well studied proteins in the superfamily of small GTPases, which has representatives in a wide range of cellular functions. These proteins share a GTP-binding pocket with highly conserved motifs that promote hydrolysis to GDP. The neutron crystal structure of RAS presented here strongly supports a protonated γ-phosphate at physiological pH. This counters the notion that the phosphate groups of GTP are fully deprotonated at the start of the hydrolysis reaction, which has colored the interpretation of experimental and computational data in studies of the hydrolysis mechanism. The neutron crystal structure presented here puts in question our understanding of the pre-catalytic state associated with the hydrolysis reaction central to the function of RAS and other GTPases.

Tyrosine phosphorylation of RAS by ABL allosterically enhances effector binding.

RAS proteins are signal transduction gatekeepers that mediate cell growth, survival, and differentiation through interactions with multiple effector proteins. The RAS effector RAS- and RAB-interacting protein 1 (RIN1) activates its own downstream effectors, the small GTPase RAB5 and the tyrosine kinase Abelson tyrosine-protein kinase (ABL), to modulate endocytosis and cytoskeleton remodeling. To identify ABL substrates downstream of RAS-to-RIN1 signaling, we examined human HEK293T cells overexpressing components of this pathway. Proteomic analysis revealed several novel phosphotyrosine peptides, including Harvey rat sarcoma oncogene (HRAS)-pTyr(137). Here we report that ABL phosphorylates tyrosine 137 of H-, K-, and NRAS. Increased RIN1 levels enhanced HRAS-Tyr(137) phosphorylation by nearly 5-fold, suggesting that RAS-stimulated RIN1 can drive ABL-mediated RAS modification in a feedback circuit. Tyr(137) is well conserved among RAS orthologs and is part of a transprotein H-bond network. Crystal structures of HRAS(Y137F) and HRAS(Y137E) revealed conformation changes radiating from the mutated residue. Although consistent with Tyr(137) participation in allosteric control of HRAS function, the mutations did not alter intrinsic GTP hydrolysis rates in vitro. HRAS-Tyr(137) phosphorylation enhanced HRAS signaling capacity in cells, however, as reflected by a 4-fold increase in the association of phosphorylated HRAS(G12V) with its effector protein RAF proto-oncogene serine/threonine protein kinase 1 (RAF1). These data suggest that RAS phosphorylation at Tyr(137) allosterically alters protein conformation and effector binding, providing a mechanism for effector-initiated modulation of RAS signaling.

Structurally conserved Nop56/58 N-terminal domain facilitates archaeal box C/D ribonucleoprotein-guided methyltransferase activity.

Box C/D RNA-protein complexes (RNPs) guide the 2'-O-methylation of nucleotides in both archaeal and eukaryotic ribosomal RNAs. The archaeal box C/D and C'/D' RNP subcomplexes are each assembled with three sRNP core proteins. The archaeal Nop56/58 core protein mediates crucial protein-protein interactions required for both sRNP assembly and the methyltransferase reaction by bridging the L7Ae and fibrillarin core proteins. The interaction of Methanocaldococcus jannaschii (Mj) Nop56/58 with the methyltransferase fibrillarin has been investigated using site-directed mutagenesis of specific amino acids in the N-terminal domain of Nop56/58 that interacts with fibrillarin. Extensive mutagenesis revealed an unusually strong Nop56/58-fibrillarin interaction. Only deletion of the NTD itself prevented dimerization with fibrillarin. The extreme stability of the Nop56/58-fibrillarin heterodimer was confirmed in both chemical and thermal denaturation analyses. However, mutations that did not affect Nop56/58 binding to fibrillarin or sRNP assembly nevertheless disrupted sRNP-guided nucleotide modification, revealing a role for Nop56/58 in methyltransferase activity. This conclusion was supported with the cross-linking of Nop56/58 to the target RNA substrate. The Mj Nop56/58 NTD was further characterized by solving its three-dimensional crystal structure to a resolution of 1.7 Å. Despite low primary sequence conservation among the archaeal Nop56/58 homologs, the overall structure of the archaeal NTD domain is very well conserved. In conclusion, the archaeal Nop56/58 NTD exhibits a conserved domain structure whose exceptionally stable interaction with fibrillarin plays a role in both RNP assembly and methyltransferase activity.

Allosteric modulation of Ras positions Q61 for a direct role in catalysis.

Ras and its effector Raf are key mediators of the Ras/Raf/MEK/ERK signal transduction pathway. Mutants of residue Q61 impair the GTPase activity of Ras and are found prominently in human cancers. Yet the mechanism through which Q61 contributes to catalysis has been elusive. It is thought to position the catalytic water molecule for nucleophilic attack on the gamma-phosphate of GTP. However, we previously solved the structure of Ras from crystals with symmetry of the space group R32 in which switch II is disordered and found that the catalytic water molecule is present. Here we present a structure of wild-type Ras with calcium acetate from the crystallization mother liquor bound at a site remote from the active site and likely near the membrane. This results in a shift in helix 3/loop 7 and a network of H-bonding interactions that propagates across the molecule, culminating in the ordering of switch II and placement of Q61 in the active site in a previously unobserved conformation. This structure suggests a direct catalytic role for Q61 where it interacts with a water molecule that bridges one of the gamma-phosphate oxygen atoms to the hydroxyl group of Y32 to stabilize the transition state of the hydrolysis reaction. We propose that Raf together with the binding of Ca(2+) and a negatively charged group mimicked in our structure by the acetate molecule induces the ordering of switch I and switch II to complete the active site of Ras.

Allosteric site variants affect GTP hydrolysis on Ras.

RAS GTPases are proto-oncoproteins that regulate cell growth, proliferation, and differentiation in response to extracellular signals. The signaling functions of RAS, and other small GTPases, are dependent on their ability to cycle between GDP-bound and GTP-bound states. Structural analyses suggest that GTP hydrolysis catalyzed by HRAS can be regulated by an allosteric site located between helices 3, 4, and loop 7. Here we explore the relationship between intrinsic GTP hydrolysis on HRAS and the position of helix 3 and loop 7 through manipulation of the allosteric site, showing that the two sites are functionally connected. We generated several hydrophobic mutations in the allosteric site of HRAS to promote shifts in helix 3 relative to helix 4. By combining crystallography and enzymology to study these mutants, we show that closure of the allosteric site correlates with increased hydrolysis of GTP on HRAS in solution. Interestingly, binding to the RAS binding domain of RAF kinase (RAF-RBD) inhibits GTP hydrolysis in the mutants. This behavior may be representative of a cluster of mutations found in human tumors, which potentially cooperate with RAF complex formation to stabilize the GTP-bound state of RAS.

Shift in the equilibrium between on and off states of the allosteric switch in Ras-GppNHp affected by small molecules and bulk solvent composition.

Ras GTPase cycles between its active GTP-bound form promoted by GEFs and its inactive GDP-bound form promoted by GAPs to affect the control of various cellular functions. It is becoming increasingly apparent that subtle regulation of the GTP-bound active state may occur through promotion of substates mediated by an allosteric switch mechanism that induces a disorder to order transition in switch II upon ligand binding at an allosteric site. We show with high-resolution structures that calcium acetate and either dithioerythritol (DTE) or dithiothreitol (DTT) soaked into H-Ras-GppNHp crystals in the presence of a moderate amount of poly(ethylene glycol) (PEG) can selectively shift the equilibrium to the "on" state, where the active site appears to be poised for catalysis (calcium acetate), or to what we call the "ordered off" state, which is associated with an anticatalytic conformation (DTE or DTT). We also show that the equilibrium is reversible in our crystals and dependent on the nature of the small molecule present. Calcium acetate binding in the allosteric site stabilizes the conformation observed in the H-Ras-GppNHp/NOR1A complex, and PEG, DTE, and DTT stabilize the anticatalytic conformation observed in the complex between the Ras homologue Ran and Importin-ß. The small molecules are therefore selecting biologically relevant conformations in the crystal that are sampled by the disordered switch II in the uncomplexed GTP-bound form of H-Ras. In the presence of a large amount of PEG, the ordered off conformation predominates, whereas in solution, in the absence of PEG, switch regions appear to remain disordered in what we call the off state, unable to bind DTE.

Allosteric modulation of Ras-GTP is linked to signal transduction through RAF kinase.

Ras is a key signal transduction protein in the cell. Mutants of Gly(12) and Gln(61) impair GTPase activity and are found prominently in cancers. In wild type Ras-GTP, an allosteric switch promotes disorder to order transition in switch II, placing Gln(61) in the active site. We show that the "on" and "off" conformations of the allosteric switch can also be attained in RasG12V and RasQ61L. Although both mutants have similarly impaired active sites in the on state, RasQ61L stabilizes an anti-catalytic conformation of switch II in the off state of the allosteric switch when bound to Raf. This translates into more potent activation of the MAPK pathway involving Ras, Raf kinase, MEK, and ERK (Ras/Raf/MEK/ERK) in cells transfected with RasQ61L relative to RasG12V. This differential is not observed in the Raf-independent pathway involving Ras, phosphoinositide 3-kinase (PI3K), and Akt (Ras/PI3K/Akt). Using a combination of structural analysis, hydrolysis rates, and experiments in NIH-3T3 cells, we link the allosteric switch to the control of signaling in the Ras/Raf/MEK/ERK pathway, supporting a GTPase-activating protein-independent model for duration of the Ras-Raf complex.

Computational studies of the principle of dynamic-change-driven protein interactions.

Dynamic allostery emphasizes a role of entropy change manifested as a sole change in protein fluctuations without structural changes. This kind of entropy-driven effect remains largely understudied. The most significant examples involve protein-ligand interactions, leaving protein-protein interactions, which are critical in signaling and other cellular events, largely unexplored. Here we study an example of how protein-protein interaction (binding of Ras to the Ras binding domain [RBD] of the effector protein Raf) affects a subsequent protein association process (Ras dimerization) by quenching Ras internal motions through dynamic allostery. We also investigate the influence of point mutations or ambient temperature, respectively, on the protein dynamics and interaction of two other systems: in adenylate kinase (ADK) and in the EphA2 SAM:Ship2 SAM complex. Based on these examples, we postulate that there are different ways in which dynamic-change-driven protein interactions are manifested and that it is likely a general biological phenomenon.

Comparative analysis of the 15.5kD box C/D snoRNP core protein in the primitive eukaryote Giardia lamblia reveals unique structural and functional features.

Box C/D ribonucleoproteins (RNP) guide the 2'-O-methylation of targeted nucleotides in archaeal and eukaryotic rRNAs. The archaeal L7Ae and eukaryotic 15.5kD box C/D RNP core protein homologues initiate RNP assembly by recognizing kink-turn (K-turn) motifs. The crystal structure of the 15.5kD core protein from the primitive eukaryote Giardia lamblia is described here to a resolution of 1.8 Å. The Giardia 15.5kD protein exhibits the typical α-ß-α sandwich fold exhibited by both archaeal L7Ae and eukaryotic 15.5kD proteins. Characteristic of eukaryotic homologues, the Giardia 15.5kD protein binds the K-turn motif but not the variant K-loop motif. The highly conserved residues of loop 9, critical for RNA binding, also exhibit conformations similar to those of the human 15.5kD protein when bound to the K-turn motif. However, comparative sequence analysis indicated a distinct evolutionary position between Archaea and Eukarya. Indeed, assessment of the Giardia 15.5kD protein in denaturing experiments demonstrated an intermediate stability in protein structure when compared with that of the eukaryotic mouse 15.5kD and archaeal Methanocaldococcus jannaschii L7Ae proteins. Most notable was the ability of the Giardia 15.5kD protein to assemble in vitro a catalytically active chimeric box C/D RNP utilizing the archaeal M. jannaschii Nop56/58 and fibrillarin core proteins. In contrast, a catalytically competent chimeric RNP could not be assembled using the mouse 15.5kD protein. Collectively, these analyses suggest that the G. lamblia 15.5kD protein occupies a unique position in the evolution of this box C/D RNP core protein retaining structural and functional features characteristic of both archaeal L7Ae and higher eukaryotic 15.5kD homologues.

Thermodynamic, enzymatic and structural effects of removing a salt bridge at the base of loop 4 in (pro)caspase-3.

Interactions between loops 2, 2' and 4, known as the loop bundle, stabilize the active site of caspase-3. Loop 4 (L4) is of particular interest due to its location between the active site and the dimer interface. We have disrupted a salt bridge between K242 and E246 at the base of L4 to determine its role in overall conformational stability and in maintaining the active site environment. Stability measurements show that only the K242A single mutant decreases stability of the dimer, whereas both single mutants and the double mutant demonstrate much lower activity compared to wild-type caspase-3. Structural studies of the caspase-3 variants show the involvement of K242 in hydrophobic interactions that stabilize helix 5, near the dimer interface, and the role of E246 appears to be to neutralize the positive charge of K242 within the hydrophobic cluster. Overall, the results suggest E246 and K242 are important in procaspase-3 for their interaction with neighboring residues, not with one another. Conversely, formation of the K242-E246 salt bridge in caspase-3 is needed for an accurate, stable conformation of loop L4 and proper active site formation in the mature enzyme.

Crystal Structure Reveals the Full Ras-Raf Interface and Advances Mechanistic Understanding of Raf Activation.

Ras and Raf-kinase interact through the Ras-binding (RBD) and cysteine-rich domains (CRD) of Raf to signal through the mitogen-activated protein kinase pathway, yet the molecular mechanism leading to Raf activation has remained elusive. We present the 2.8 Å crystal structure of the HRas-CRaf-RBD_CRD complex showing the Ras-Raf interface as a continuous surface on Ras, as seen in the KRas-CRaf-RBD_CRD structure. In molecular dynamics simulations of a Ras dimer model formed through the α4-α5 interface, the CRD is dynamic and located between the two Ras protomers, poised for direct or allosteric modulation of functionally relevant regions of Ras and Raf. We propose a molecular model in which Ras binding is involved in the release of Raf autoinhibition while the Ras-Raf complex dimerizes to promote a platform for signal amplification, with Raf-CRD centrally located to impact regulation and function.

Fragment-based identification of druggable 'hot spots' of proteins using Fourier domain correlation techniques.

MOTIVATION: The binding sites of proteins generally contain smaller regions that provide major contributions to the binding free energy and hence are the prime targets in drug design. Screening libraries of fragment-sized compounds by NMR or X-ray crystallography demonstrates that such 'hot spot' regions bind a large variety of small organic molecules, and that a relatively high 'hit rate' is predictive of target sites that are likely to bind drug-like ligands with high affinity. Our goal is to determine the 'hot spots' computationally rather than experimentally. RESULTS: We have developed the FTMAP algorithm that performs global search of the entire protein surface for regions that bind a number of small organic probe molecules. The search is based on the extremely efficient fast Fourier transform (FFT) correlation approach which can sample billions of probe positions on dense translational and rotational grids, but can use only sums of correlation functions for scoring and hence is generally restricted to very simple energy expressions. The novelty of FTMAP is that we were able to incorporate and represent on grids a detailed energy expression, resulting in a very accurate identification of low-energy probe clusters. Overlapping clusters of different probes are defined as consensus sites (CSs). We show that the largest CS is generally located at the most important subsite of the protein binding site, and the nearby smaller CSs identify other important subsites. Mapping results are presented for elastase whose structure has been solved in aqueous solutions of eight organic solvents, and we show that FTMAP provides very similar information. The second application is to renin, a long-standing pharmaceutical target for the treatment of hypertension, and we show that the major CSs trace out the shape of the first approved renin inhibitor, aliskiren. AVAILABILITY: FTMAP is available as a server at http://ftmap.bu.edu/.

A constitutively active and uninhibitable caspase-3 zymogen efficiently induces apoptosis.

The caspase-3 zymogen has essentially zero activity until it is cleaved by initiator caspases during apoptosis. However, a mutation of V266E in the dimer interface activates the protease in the absence of chain cleavage. We show that low concentrations of the pseudo-activated procaspase-3 kill mammalian cells rapidly and, importantly, this protein is not cleaved nor is it inhibited efficiently by the endogenous regulator XIAP (X-linked inhibitor of apoptosis). The 1.63 A (1 A = 0.1 nm) structure of the variant demonstrates that the mutation is accommodated at the dimer interface to generate an enzyme with substantially the same activity and specificity as wild-type caspase-3. Structural modelling predicts that the interface mutation prevents the intersubunit linker from binding in the dimer interface, allowing the active sites to form in the procaspase in the absence of cleavage. The direct activation of procaspase-3 through a conformational switch rather than by chain cleavage may lead to novel therapeutic strategies for inducing cell death.