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

Comparative MD simulations and advanced analytics based studies on wild-type and hot-spot mutant A59G HRas

The Ras family of proteins is known to play an important role in cellular signal transduction. The oncoprotein Ras is also found to be mutated in ~90% of the pancreatic cancers, of which G12V, G13V, A59G and Q61L are the known hot-spot mutants. These ubiquitous proteins fall in the family of G-proteins, and hence switches between active GTP bound and inactive GDP bound states, which is hindered in most of its oncogenic mutant counterparts. Moreover, Ras being a GTPase has an intrinsic property to hydrolyze GTP to GDP, which is obstructed due to mutations and lends the mutants stuck in constitutively active state leading to oncogenic behavior. In this regard, the present study aims to understand the dynamics involved in the hot-spot mutant A59G-Ras using long 10μs classical MD simulations (5μs for each of the wild-type and mutant systems) and comparing the same with its wild-type counterpart. Advanced analytics using Markov State Model (MSM) based approach has been deployed to comparatively understand the transition path for the wild-type and mutant systems. Roles of crucial residues like Tyr32, Gln61 and Tyr64 have also been established using multivariate PCA analyses. Furthermore, this multivariate PCA analysis also provides crucial features which may be used as reaction coordinates for biased simulations for further studies. The absence of formation of pre-hydrolysis network is also reported for the mutant conformation, using the distance-based analyses (between crucial residues) of the conserved regions. The implications of this study strengthen the hypothesis that the disruption of the pre-hydrolysis network in the mutant A59G ensemble might lead to permanently active oncogenic conformation in the mutant conformers.

Differentiating the pre-hydrolysis states of wild-type and A59G mutant HRas: An insight through MD simulations

The most representative member of the Ras subfamily is its HRas isoform. Ras proteins being GTPases, possess an intrinsic activity to hydrolyze the GTP molecule to GDP. During the transition phases, between active and inactive states, P-loop and switch regions show maximum variations. Various hot-spot Ras mutants (G12V, A59G, Q61L etc) have been reported, that limit the protein's conformation in the permanent active state. In the present study, we aim to explore the structural dynamics of one such crucial mutant of Ras namely A59G which belongs to the conserved Switch II region of the protein. Approximately ∼15μs of Classical Molecular Dynamics (CMD) simulations have been carried out on the mutant and wild-type complexes. Further, a metadynamics simulation of 500ns was also carried out, which suggests an energy barrier of ∼9.56kcal/mol between wild-type and mutant conformation. We demonstrate the role of water molecule in maintaining the required interaction networks in the pre-hydrolysis state, its impact on A59G mutation, distinct orientation of the Gln61 residue in two conformations, disruption of crucial Gly60 and γ phosphate and the change in the Switch II region. The outcome of our study captures the pre-hydrolysis state of the HRas protein. It also establishes the fact that this mutation makes the movement of Switch II region and the conserved DXXGQ motif highly constrained, which is known to be an important requirement for hydrolysis. This suggests that the A59G mutation may decrease the rate of intrinsic hydrolysis as well as GAP-mediated hydrolysis.

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(2+) and water-forming stable interactions.

Tryptophan rotamer distribution revealed for the α-helix in tear lipocalin by site-directed tryptophan fluorescence

Rotamer libraries are a valuable tool for protein structure determination, modeling, and design. Site-directed tryptophan fluorescence (SDTF) was used in combination with the rotamer model for the fluorescence intensity decays to solve α-helical conformations of proteins in solution. Single Trp mutations located in an α-helical segment of human tear lipocalin were explored for structure assignment. Along with fluorescence λ(max) values, the rotamer model assignment of fluorescence lifetimes fits the backbone conformation. Typically, Trp fluorescence in proteins shows three lifetimes. However, for the α-helix, two lifetimes assigned to t and g(-) rotamers were satisfactory to describe Trp fluorescence intensity decays. The g(+) rotamer is not feasible in the α-helix due to steric restriction. Trp rotamer distributions obtained by fluorescence were compared with the rotamer library derived from X-ray crystallography data of proteins. The Trp rotamer distributions vary for solvent exposed and buried (tertiary interaction) sites. A new strategy using the rotamer distribution with SDTF (RD-SDTF) removes the limitation of regular SDTF and other labeling techniques, in which site-specific differences, e.g., accessibility, are presumed. The RD-SDTF technique does not rely on environmental differences of side chains and is able to detect α-helical structure where all side chains are exposed to solvent. Potentially, this technique is applicable to various proteins including membrane proteins, which are rich in α-helix motif.

Excited protein states of human tear lipocalin for low- and high-affinity ligand binding revealed by functional AB loop motion

Human tear lipocalin (TL), a prominent member of lipocalin family, exhibits functional and structural promiscuity. The plasticity of loop regions modulates entry to the ligand pocket at the "open" end of the eight-stranded beta-barrel. Site-directed multi-distance measurements using fluorescence resonance energy transfer between functional loops register two excited protein states for low- and high-affinity ligand binding. At low pH, the longest loop AB adopts the conformation of the low-affinity excited protein state that matches the crystal structure of holo-TL at pH 8. A "crankshaft" like movement is detected for the loop AB in a low pH transition. At pH 7.3 the holo-protein assumes a high-affinity excited protein state, in which the loop AB is more compact (RMS=3.1A). In the apo-holo transition, the reporter Trp 28 moves about 4.5A that reflects a decrease in distance between Glu27 and Lys108. This interaction fixes the loop AB conformation for the high-affinity mode. No such movement is detected at low pH, where Glu27 is protonated. Data strongly indicate that the protonation state of Glu27 modulates the conformation of the loop AB for high- and low-affinity binding.

Reversible movement of switch 1 loop of myosin determines actin interaction

The conserved switch 1 loop of P-loop NTPases is implicated as a central element that transmits information between the nucleotide-binding pocket and the binding site of the partner proteins. Recent structural studies have identified two states of switch 1 in G-proteins and myosin, but their role in the transduction mechanism has yet to be clarified. Single tryptophan residues were introduced into the switch 1 region of myosin II motor domain and studied by rapid reaction methods. We found that in the presence of MgADP, two states of switch 1 exist in dynamic equilibrium. Actin binding shifts the equilibrium towards one of the MgADP states, whereas ATP strongly favors the other. In the light of electron cryo-microscopic and X-ray crystallographic results, these findings lead to a specific structural model in which the equilibrium constant between the two states of switch 1 is coupled to the strength of the actin-myosin interaction. This has implications for the enzymatic mechanism of G-proteins and possibly P-loop NTPases in general.

Cdc42p GDP/GTP cycling is necessary for efficient cell fusion during yeast mating

The highly conserved small Rho G-protein, Cdc42p plays a critical role in cell polarity and cytoskeleton organization in all eukaryotes. In the yeast Saccharomyces cerevisiae, Cdc42p is important for cell polarity establishment, septin ring assembly, and pheromone-dependent MAP-kinase signaling during the yeast mating process. In this study, we further investigated the role of Cdc42p in the mating process by screening for specific mating defective cdc42 alleles. We have identified and characterized novel mating defective cdc42 alleles that are unaffected in vegetative cell polarity. Replacement of the Cdc42p Val36 residue with Met resulted in a specific cell fusion defect. This cdc42[V36M] mutant responded to mating pheromone but was defective in cell fusion and in localization of the cell fusion protein Fus1p, similar to a previously isolated cdc24 (cdc24-m6) mutant. Overexpression of a fast cycling Cdc42p mutant suppressed the cdc24-m6 fusion defect and conversely, overexpression of Cdc24p suppressed the cdc42[V36M] fusion defect. Taken together, our results indicate that Cdc42p GDP-GTP cycling is critical for efficient cell fusion.

Synthesis of GTP-derived Ras ligands

A practical and convenient method for the synthesis of acid- and base-sensitive GTP analogues carrying a further substituent at the terminal phosphate has been developed. Key to the successful synthesis of these potential ligands of the Ras protein is the use of Pd0-sensitive allyl protecting groups in a one-pot synthesis that avoids evaporation steps. Initial biochemical analysis of a representative compound revealed that such GTP analogues can bind to Ras and might open up the possibility of developing small molecules that can act as deactivators of oncogenic Ras.