Relation between the conformational heterogeneity and reaction cycle of Ras: molecular simulation of Ras

Ras functions as a molecular switch by cycling between the active GTP-bound state and the inactive GDP-bound state. It is known experimentally that there is another GTP-bound state called state 1. We investigate the conformational changes and fluctuations arising from the difference in the coordinations between the switch regions and ligands in the GTP- and GDP-bound states using a total of 830 ns of molecular-dynamics simulations. Our results suggest that the large fluctuations among multiple conformations of switch I in state 1 owing to the absence of coordination between Thr-35 and Mg²⁺ inhibit the binding of Ras to effectors. Furthermore, we elucidate the conformational heterogeneity in Ras by using principal component analysis, and propose a two-step reaction path from the GDP-bound state to the active GTP-bound state via state 1. This study suggests that state 1 plays an important role in signal transduction as an intermediate state of the nucleotide exchange process, although state 1 itself is an inactive state for signal transduction.     

Dynamic polymorphism of Ras observed by single molecule FRET is the basis for molecular recognition

Ras regulates signal transduction pathway function by dynamically interacting with various effectors. To understand the basis for Ras function, its conformational dynamics were measured in the absence and presence of effectors using single molecule fluorescence resonance energy transfer (FRET) between probes located on the Switch II region and GTP. The time trajectories of FRET efficiency from GTP-bound Ras showed that this conformation spontaneously varies among multiple states. Among them, a low FRET state was identified as an inactive state. The transition involving the inactive conformational state occurred in the time range of seconds. In contrast, fluctuation occurring most probably between multiple active high FRET conformational states lasted approximately 30 ms but converged to a specific conformational state upon binding to an effector. Thus, Ras conformation spontaneously fluctuates to readily interact with various effectors.     

Active and Inactive Cdc42 Differ in Their Insert Region Conformational Dynamics

Cell division control protein 42 homolog (Cdc42) protein, a Ras superfamily GTPase, regulates cellular activities, including cancer progression. Using all-atom molecular dynamics (MD) simulations and essential dynamic analysis, we investigated the structure and dynamics of the catalytic domains of GDP-bound (inactive) and GTP-bound (active) Cdc42 in solution. We discovered substantial differences in the dynamics of the inactive and active forms, particularly in the “insert region” (residues 122–135), which plays a role in Cdc42 activation and binding to effectors. The insert region has larger conformational flexibility in the GDP-bound Cdc42 than in the GTP-bound Cdc42. The G2 loop and switch I at the effector lobe of the catalytic domain exhibit large conformational changes in both the GDP- and the GTP-bound systems, but in the GTP-bound Cdc42, the switch I interactions with GTP are retained. Oncogenic mutations were identified in the Ras superfamily. In Cdc42, the G12V and Q61L mutations decrease the GTPase activity. We simulated these mutations in both GDP- and GTP-bound Cdc42. Although the overall structural organization is quite similar between the wild type and the mutants, there are small differences in the conformational dynamics, especially in the two switch regions. Taken together, the G12V and Q61L mutations may play a role similar to their K-Ras counterparts in nucleotide binding and activation. The conformational differences, which are mainly in the insert region and, to a lesser extent, in the switch regions flanking the nucleotide binding site, can shed light on binding and activation. We propose that the differences are due to a network of hydrogen bonds that gets disrupted when Cdc42 is bound to GDP, a disruption that does not exist in other Rho GTPases. The differences in the dynamics between the two Cdc42 states suggest that the inactive conformation has reduced ability to bind to effectors.     

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

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

Interconversion of active and inactive 30 S ribosomal subunits is accompanied by a conformational change in the decoding region of 16 S rRNA

Zamir, Elson and their co-workers have shown that 30 S ribosomal subunits are reversibly inactivated by depletion of monovalent or divalent cations. We have re-investigated the conformation of 16 S rRNA in the active and inactive forms of the 30 S subunit, using a strategy that is designed to eliminate reversible ion-dependent conformational effects that are unrelated to the heat-dependent Zamir-Elson transition. A combination of structure-specific chemical probes enables us to monitor the accessibility of pyrimidines at N-3 and purines at N-1 and N-7. Chemically modified bases are identified by end-labeling followed by analine-induced strand scission (in some cases preceded by hybrid selection), or by primer extension using synthetic DNA oligomers. These studies show the following: The transition from the active to the inactive state cannot be described as a simple loosening or unfolding of native structure, such as that which is observed under conditions of more severe ion depletion. Instead, it has the appearance of a reciprocal interconversion between two differently structured states; some bases become more reactive toward the probes, whilst others become less reactive as a result of inactivation. Changes in reactivity are almost exclusively confined to the "decoding site" centered at positions 1400 and 1500, but significant differences are also detected at U723 and G791 in the central domain. This may reflect possible structural and functional interactions between the central and 3' regions of 16 S rRNA. The inactive form also shows significantly decreased reactivity at positions 1533 to 1538 (the Shine-Dalgarno region), in agreement with earlier findings. The principal changes in reactivity involve the universally conserved nucleotides G926, C1395, A1398 and G1401. The three purines show reciprocal behavior at their N-1 versus N-7 positions. G926 loses its reactivity at N-1, but becomes highly reactive at N-7 as a result of the transition of the inactive state. In contrast, A1398 and G1401 become reactive at N-1, but lose their hyper-reactivity at N-7. The possible structural and functional implications of these findings are discussed.     

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

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

Characterization of a guanine nucleotide dissociation stimulator for a ras-related GTPase.

ras-related GTPases participate in signaling for a variety of cellular processes. The GTPases cycle between a GTP-bound active state and a GDP-bound inactive state. This cycling is partially controlled by guanine nucleotide dissociation stimulators (GDS, also known as exchange factors). We report on the molecular cloning of cDNAs encoding a new mammalian GDS protein, using sequences derived from the yeast ras GDS proteins as probes. The encoded protein stimulates the dissociation of guanine nucleotides from the ras-related ralA and ralB GTPases at a rate at least 30-fold faster than the intrinsic nucleotide dissociation rate. This new GDS, ralGDS, is at least 20-fold more active on the ralA and ralB GTPases than on any other GTPase tested, including other members of the ras family (H-ras, N-ras, K-ras, R-ras, rap1a and rap2), members of the rho family (rhoA, rhoB and CDC42-Hs) and members of the rab family (rab3a and ypt1). While the ralGDS protein is phosphorylated on serine residues, we find no evidence that phosphorylation affects the activity of insect cell-expressed ralGDS towards the ralA or ralB GTPase. The 3600 nucleotide ralGDS mRNA and the 115 kDa protein were found in all tissues and cell lines examined.     

Direct visualization of Ras proteins in spatially distinct cell surface microdomains

Localization of signaling complexes to specific microdomains coordinates signal transduction at the plasma membrane. Using immunogold electron microscopy of plasma membrane sheets coupled with spatial point pattern analysis, we have visualized morphologically featureless microdomains, including lipid rafts, in situ and at high resolution. We find that an inner-plasma membrane lipid raft marker displays cholesterol-dependent clustering in microdomains with a mean diameter of 44 nm that occupy 35% of the cell surface. Cross-linking an outer-leaflet raft protein results in the redistribution of inner leaflet rafts, but they retain their modular structure. Analysis of Ras microlocalization shows that inactive H-ras is distributed between lipid rafts and a cholesterol-independent microdomain. Conversely, activated H-ras and K-ras reside predominantly in nonoverlapping, cholesterol-independent microdomains. Galectin-1 stabilizes the association of activated H-ras with these nonraft microdomains, whereas K-ras clustering is supported by farnesylation, but not geranylgeranylation. These results illustrate that the inner plasma membrane comprises a complex mosaic of discrete microdomains. Differential spatial localization within this framework can likely account for the distinct signal outputs from the highly homologous Ras proteins.     

Conformation and stability of alpha-helical membrane proteins. 1. Influence of salts on conformational equilibria between active and Inactive states of rhodopsin

We studied the influence of salts on the pH-dependent conformational equilibria between the active and the inactive photoproduct states of rhodopsin, Meta II and Meta I, respectively, and between the active and inactive conformations of the apoprotein opsin. In both equilibria, the active species is favored in the presence of medium to high concentration of salt. The ion selectivity for the Meta I/Meta II equilibrium is particularly pronounced for the anions and follows the series trichloroacetate > thiocyanate > iodide > bromide > sulfate > chloride > acetate. The Hill coefficient of this salt-induced transition is close to 2.0. Both ion selectivity and Hill coefficient suggest that the transition is mainly regulated by ion binding to two specific charged binding sites in the protein with smaller contributions being due to the Hofmeister effect. We propose that these putative ion binding sites are identical to those sites that are titrated in the corresponding pH-dependent conformational transition. They presumably function as ionic locks, which keep the receptor in an inactive conformation, and which may be disrupted either by pH-dependent protonation or by salt-dependent ion binding.     

Coarse-grained dynamics of the receiver domain of NtrC: fluctuations, correlations and implications for allosteric cooperativity

Receiver domains are key molecular switches in bacterial signaling. Structural studies have shown that the receiver domain of the nitrogen regulatory protein C (NtrC) exists in a conformational equilibrium encompassing both inactive and active states, with phosphorylation of Asp54 allosterically shifting the equilibrium towards the active state. To analyze dynamical fluctuations and correlations in NtrC as it undergoes activation, we have applied a coarse-grained dynamics algorithm using elastic network models. Normal mode analysis reveals possible dynamical pathways for the transition of NtrC from the inactive state to the active state. The diagonalized correlation between the inactive and the active (phosphorylated) state shows that most correlated motions occur around the active site of Asp54 and in the region Thr82 to Tyr101. This indicates a coupled correlation of dynamics in the "Thr82-Tyr101" motion. With phosphorylation inducing significant flexibility changes around the active site and alpha3 and alpha4 helices, we find that this activation makes the active-site region and the loops of alpha3/beta4 and alpha4/beta5 more stable. This means that phosphorylation entropically favors the receiver domain in its active state, and the induced conformational changes occur in an allosteric manner. Analyses of the local flexibility and long-range correlated motion also suggest a dynamics criterion for determining the allosteric cooperativity of NtrC, and may be applicable to other proteins.     

Rho小G蛋白介导的细胞周期调控

Rho小G蛋白作为一个信号分子家族具有多样化的功能, 可以调节细胞骨架重排 、细胞迁移、细胞极性、基因表达、细胞周期调控等. Rho小G蛋白家族对细胞周期 调控的研究主要集中在其对于有丝分裂期细胞的调节作用,包括调节有丝分裂期前 期细胞趋圆化、后期染色体排列及收缩环的收缩作用.近期的研究显示,Rho小G蛋白及其效应分子对于细胞周期G1、S、G2期的调控主要是通过影响细胞周期的正调控因子细胞周期蛋白D1 (cyclin D1) 和负调控因子细胞周期蛋白依赖型激酶相互作用蛋白1及细胞周期蛋白依赖型激酶抑制蛋白27 (p21cip1/p27kip1) 进行的.本文总结了Rho小G蛋白及其效应分子在细胞周期调控,尤其是对G1/S期调控的研究进展,并简要阐述了Rho小G蛋白介导的细胞周期调控异常与癌症发生的关系.     

Conformational states of ADP ribosylation factor 1 complexed with different guanosine triphosphates as studied by 31P NMR spectroscopy

Guanine nucleotide binding proteins (GNB-proteins) play an essential role in cellular signaling, acting as molecular switches, cycling between the inactive, GDP-bound form and the active, GTP-bound form. It has been shown that conformational equilibria also exist within the active form of GNB-proteins between conformational states with different functional properties. Here we present (31)P NMR data on ADP ribosylation factor 1 (Arf1), a GNB-protein involved in Golgi traffic, promoting the coating of secretory vesicles. To investigate conformational equilibria in active Arf1, the wild type and switch I mutants complexed with GTP and a variety of commonly used GTP analogues, namely, GppCH(2)p, GppNHp, and GTPγS, were analyzed. To gain deeper insight into the conformational state of active Arf1, we titrated with Cu(2+)-cyclen and GdmCl and formed the complex with the Sec7 domain of nucleotide exchange factor ARNO and an effector GAT domain. In contrast to the related proteins Ras, Ral, Cdc42, and Ran, from (31)P NMR spectroscopic view, Arf1 exists predominantly in a single conformation independent of the GTP analogue used. This state seems to correspond to the so-called state 2(T) conformation, according to Ras nomenclature, which is interacting with the effector domain. The exchange of the highly conserved threonine in position 48 with alanine led to a shift of the equilibrium toward a conformational state with typical properties obtained for state 1(T) in Ras, such as interaction with guanine nucleotide exchange factors, a lower affinity for nucleoside triphosphates, and greater sensitivity to chaotropic agents. In active Arf1(wt), the effector interacting conformation is strongly favored. These intrinsic conformational equilibria of active GNB-proteins could be a fine-tuning mechanism of regulation and thereby an interesting target for the modulation of protein activity.     

Analysis of Ras:RasGEF interactions by phage display and static multi-angle light scattering

Molecular switches such as small GTPases of the Ras family cycle between inactive GDP-bound and active GTP-bound states. Their essential role in controlling development and cell homeostasis requires mechanisms which determine amplitude and timing of activation. This is achieved in part by the action of guanine nucleotide exchange factors, which function as highly controlled enzymes whose activity relies on spatial segregation and intra- and intermolecular regulation. Here, we describe two experimental methodologies that permit the identification and characterization of GTPase binding sites on activators by assaying complex formation within a broad range of affinities. In the first assay system, proteins presented on the surface of filamentous phage are used to probe affinity determinants of protein-protein interactions. In this application, a protein-displayed phage library is generated by random mutagenesis and a plate-based selection is performed to identify mutations that confer higher binding affinity to an immobilized target. The second method uses light scattering as a tool for measuring the molecular weight, stoichiometry, and polydispersity of protein complexes in solution. In this application, conventional gel filtration chromatography provides initial fractionation, and in-line light scattering measurements allow accurate determination of molar masses of the eluent. This technique also provides information about conformational homogeneity which can be used as a quality     

Putative Active States of a Prototypic G-Protein-Coupled Receptor from Biased Molecular Dynamics

A major current focus of structural work on G-protein-coupled receptors (GPCRs) pertains to the investigation of their active states. However, for virtually all GPCRs, active agonist-bound intermediate states have been difficult to characterize experimentally owing to their higher conformational flexibility, and thus intrinsic instability, as compared to inactive inverse agonist-bound states. In this work, we explored possible activation pathways of the prototypic GPCR bovine rhodopsin by means of biased molecular dynamics simulations. Specifically, we used an explicit atomistic representation of the receptor and its environment, and sampled the conformational transition from the crystal structure of a photoactivated deprotonated state of rhodopsin to the low pH crystal structure of opsin in the presence of 11-trans-retinal, using adiabatic biased molecular dynamics simulations. We then reconstructed the system free-energy landscape along the predetermined transition trajectories using a path collective variable approach based on metadynamics. Our results suggest that the two experimental endpoints of rhodopsin/opsin are connected by at least two different pathways, and that the conformational transition is populated by at least four metastable states of the receptor, characterized by a different amplitude of the outward movement of transmembrane helix 6.     

Mapping l1 ligase ribozyme conformational switch

L1 ligase (L1L) molecular switch is an in vitro optimized synthetic allosteric ribozyme that catalyzes the regioselective formation of a 5′-to-3′ phosphodiester bond, a reaction for which there is no known naturally occurring RNA catalyst. L1L serves as a proof of principle that RNA can catalyze a critical reaction for prebiotic RNA self-replication according to the RNA world hypothesis. L1L crystal structure captures two distinct conformations that differ by a reorientation of one of the stems by around 80 Å and are presumed to correspond to the active and inactive state, respectively. It is of great interest to understand the nature of these two states in solution and the pathway for their interconversion. In this study, we use explicit solvent molecular simulation together with a novel enhanced sampling method that utilizes concepts from network theory to map out the conformational transition between active and inactive states of L1L. We find that the overall switching mechanism can be described as a three‐state/two‐step process. The first step involves a large-amplitude swing that reorients stem C. The second step involves the allosteric activation of the catalytic site through distant contacts with stem C. Using a conformational space network representation of the L1L switch transition, it is shown that the connection between the three states follows different topographical patterns: the stem C swing step passes through a narrow region of the conformational space network, whereas the allosteric activation step covers a much wider region and a more diverse set of pathways through the network.     

Molecular dynamic simulations of the N-terminal receiver domain of NtrC reveal intrinsic conformational flexibility in the inactive state

The N-terminal receiver domain of NtrC is the molecular switch in the two-component signal transduction. It is the first protein where structures of both the active (phosphyroylated) and inactive (unphosphyroylated) states are determined experimentally. Phosphorylation of the NtrC at the active site induces large structural change. NMR experiments suggested that the wild type unphosphorylated NtrC adopts both the active and the inactive conformations and the phosphorylation stabilizes the active conformations. We applied free (unconstrained) molecular dynamic (MD) simulation to examine the intrinsic flexibilities and stabilities of the NtrC receiver domain in both the active and inactive conformations. Molecular dynamic simulations showed that the inactive state of NtrC receiver domain is more flexible than the active state. There were large movements in helix 4 and loop beta3-alpha3 which coincide with major structural differences between the inactive and active states. We observed large root-mean-square deviations from the initial starting structure and the large root-mean-square fluctuations during MD simulation for the inactive state. We then investigated the activation pathway with Targeted MD simulation. We show that the intrinsic flexibility in the loop beta3-alpha3 plays an important role in triggering the conformational change. Phosphorylation at the active site may serve to stabilize the conformational change. These results together suggest that the unphosphorylated NtrC receiver domain could be involved in a conformational equilibrium between two different states.     

Current knowledge of the large RhoGAP family of proteins

The Rho GTPases are implicated in almost every fundamental cellular process. They act as molecular switches that cycle between an active GTP-bound and an inactive GDP-bound state. Their slow intrinsic GTPase activity is greatly enhanced by RhoGAPs (Rho GTPase-activating proteins), thus causing their inactivation. To date, more than 70 RhoGAPs have been identified in eukaryotes, ranging from yeast to human, and based on sequence homology of their RhoGAP domain, we have grouped them into subfamilies. In the present Review, we discuss their regulation, biological functions and implication in human diseases.     

The topology design principles that determine the spatiotemporal dynamics of G-protein cascades

Small monomeric G-proteins control cellular behavior, cycling between inactive GDP-bound and active GTP-bound states. Activating and deactivating transitions are regulated by guanine nucleotide exchange factors (GEFs) and GTPase activating proteins (GAPs), respectively. G-proteins can control different GEF and GAP activities, thereby creating GTPase signaling cascades. Here, we characterize all 128 different wiring topologies of two-layer cascades, which include feedforward/feedback interactions and an auto-regulatory loop. Exclusion of "mirror" designs leaves 64 topologies, which are classified into eight groups. We demonstrate that eight different cascades in each group generate the same number of steady states and similar spatiotemporal dynamics. Two groups (featuring 16 topologies) can generate three distinct dynamics: (i) bistable switches, (ii) excitable behavior, and (iii) sustained oscillations, giving rise to propagating waves of G-protein activation switches and pulses. Four other groups can produce switch-like, bistable behaviors and trigger waves. The remaining two groups have a single steady state. This first, complete classification of all possible interaction circuitries systematically links topological design to the spatiotemporal dynamics of G-protein cascades, predicting and explaining experimentally observed behavior.     

Structural insights into Rab21 GTPase activation mechanism by molecular dynamics simulations

Rab proteins belong to the family of monomeric GTPases which are involved in the cellular membrane trafficking. Rab21 protein exists in interchangeable GTP- and GDP-bound states. Rabs switch between two active and inactive conformations like other GTPases. The inactive form of Rab is bound to GDP while its active form is bounded with the GTP. Interexchange between active and inactive form is mediated by the GDP/GTP exchange factor (GEF) which catalyses the conversion from GDP-bound to GTP-bound form, thereby activating the Rab. While the GTP hydrolysis of Rabs is regulated by a GTPase-activating protein (GAP) which causes Rab inactivation. Here, we report the structural flexibility of the Rab21-GTP and Rab21-GDP complexes by docking and molecular dynamics (MD) simulations. Structural analysis of exchange mechanisms of the co-factors complexed with Rab21 reveals that Cys29, Thr33, His48, Gln78 and Lys133 are essentially important in the activation of proteins. Furthermore, a significant change in the orientation of the interacting co-factors, with slight variation in the free energy of binding was observed. Complexation of GEF with Rab21-GTP and Rab21-GDP reveal a flipping of the switchable residues. Finally, 50 ns MD simulations confirm that the GTP-bound Rab21 complex is thermodynamically more favoured than the corresponding GDP-bound complex. This study provides a detailed understanding of the structural elements involved in the conformational changes of Rab21.     

Conformational states of the switch I region of Ha-ras-p21 in hinge residue mutants studied by fluorescence lifetime and fluorescence anisotropy measurements

The hinge residues (Val29 and Ile36) of the switch I region (also known as the effector loop) of the Ha-ras-p21 protein have been mutated to glycines to accelerate the conformational changes typical for the effector loop. In this work, we have studied the influence of the combined mutations on the steady-state structure of the switch I region of the protein in both the inactive GDP-bound conformation as in the active GTP-bound conformation. Here, we use the fluorescence properties of the single tryptophan residue in the Y32W mutant of Ha-ras-p21. This mutant has already been used extensively as a reference form of the protein. Reducing the size of the side chains of the hinge residues not only accelerates the conformational changes but also affects the steady-state structures of the effector loop as indicated by the changes in the fluorescence properties. A thorough analysis of the fluorescence changes (quantum yield, lifetimes, etc.) proves that these changes are from a reshuffling between the rotamer populations of Trp. The population reshuffling is caused by the overall structural rearrangement along the switch I region. The effects are clearly more pronounced in the inactive GDP-bound conformation than in the active GTP-bound conformation. The effect of both mutations seems to be additive in the GDP-bound state, but cooperative in the GTP-bound state.