Thermodynamics of Ras/effector and Cdc42/effector interactions probed by isothermal titration calorimetry

Multiplexed imaging for probing RAS-RAF interactions in living cells

GTP-bound RAS interacts with its protein effectors in response to extracellular stimuli, leading to chemical inputs for downstream pathways. Significant progress has been made in measuring these reversible protein-protein interactions (PPIs) in various cell-free environments. Yet, acquiring high sensitivity in heterogeneous solutions remains challenging. Here, using an intermolecular fluorescence resonance energy transfer (FRET) biosensing approach, we develop a method to visualize and localize HRAS-CRAF interactions in living cells. We demonstrate that the EGFR activation and the HRAS-CRAF complex formation can be concurrently probed in a single cell. This biosensing strategy discriminates EGF-stimulated HRAS-CRAF interactions at the cell and organelle membranes. In addition, we provide quantitative FRET measurements for assessing these transient PPIs in a cell-free environment. Finally, we prove the utility of this approach by showing that an EGFR-binding compound is a potent inhibitor of HRAS-CRAF interactions. The outcomes of this work form a fundamental basis for further explorations of the spatiotemporal dynamics of various signaling networks.

Structural and biophysical properties of farnesylated KRas interacting with the chaperone SmgGDS-558

KRas is a small GTPase and membrane-bound signaling protein. Newly synthesized KRas is post-translationally modified with a membrane-anchoring prenyl group. KRas chaperones are therapeutic targets in cancer due to their participation in trafficking oncogenic KRas to membranes. SmgGDS splice variants are chaperones for small GTPases with basic residues in their hypervariable domain (HVR), including KRas. SmgGDS-607 escorts pre-prenylated small GTPases, while SmgGDS-558 escorts prenylated small GTPases. We provide a structural description of farnesylated and fully processed KRas (KRas-FMe) in complex with SmgGDS-558 and define biophysical properties of this interaction. Surface plasmon resonance measurements on biomimetic model membranes quantified the thermodynamics of the interaction of SmgGDS with KRas, and small-angle x-ray scattering was used to characterize complexes of SmgGDS-558 and KRas-FMe structurally. Structural models were refined using Monte Carlo and molecular dynamics simulations. Our results indicate that SmgGDS-558 interacts with the HVR and the farnesylated C-terminus of KRas-FMe, but not its G-domain. Therefore, SmgGDS-558 interacts differently with prenylated KRas than prenylated RhoA, whose G-domain was found in close contact with SmgGDS-558 in a recent crystal structure. Using immunoprecipitation assays, we show that SmgGDS-558 binds the GTP-bound, GDP-bound, and nucleotide-free forms of farnesylated and fully processed KRas in cells, consistent with SmgGDS-558 not engaging the G-domain of KRas. We found that the dissociation constant, Kd, for KRas-FMe binding to SmgGDS-558 is comparable with that for the KRas complex with PDEδ, a well-characterized KRas chaperone that also does not interact with the KRas G-domain. These results suggest that KRas interacts in similar ways with the two chaperones SmgGDS-558 and PDEδ. Therapeutic targeting of the SmgGDS-558/KRas complex might prove as useful as targeting the PDEδ/KRas complex in KRas-driven cancers.

Structure-based design of CDC42 effector interaction inhibitors for the treatment of cancer

CDC42 family GTPases (RHOJ, RHOQ, CDC42) are upregulated but rarely mutated in cancer and control both the ability of tumor cells to invade surrounding tissues and the ability of endothelial cells to vascularize tumors. Here, we use computer-aided drug design to discover a chemical entity (ARN22089) that has broad activity against a panel of cancer cell lines, inhibits S6 phosphorylation and MAPK activation, activates pro-inflammatory and apoptotic signaling, and blocks tumor growth and angiogenesis in 3D vascularized microtumor models (VMT) in vitro. Additionally, ARN22089 has a favorable pharmacokinetic profile and can inhibit the growth of BRAF mutant mouse melanomas and patient-derived xenografts in vivo. ARN22089 selectively blocks CDC42 effector interactions without affecting the binding between closely related GTPases and their downstream effectors. Taken together, we identify a class of therapeutic agents that influence tumor growth by modulating CDC42 signaling in both the tumor cell and its microenvironment.

Multiple Mechanisms Drive the Evolutionary Adaptation of Phytophthora infestans Effector Avr1 to Host Resistance

Effectors, a group of small proteins secreted by pathogens, play a central role in antagonistic interactions between plant hosts and pathogens. The evolution of effector genes threatens plant disease management and sustainable food production, but population genetic analyses to understand evolutionary mechanisms of effector genes are limited compared to molecular and functional studies. Here we investigated the evolution of the Avr1 effector gene from 111 Phytophthora infestans isolates collected from six areas covering three potato cropping regions in China using a population genetic approach. High genetic variation of the effector gene resulted from diverse mechanisms including base substitution, pre-termination, intragenic recombination and diversifying selection. Nearly 80% of the 111 sequences had a point mutation in the 512th nucleotide (T512G), which generated a pre-termination stop codon truncating 38 amino acids in the C-terminal, suggesting that the C-terminal may not be essential to ecological and biological functions of P. infestans. A significant correlation between the frequency of Avr1 sequences with the pre-termination and annual mean temperature in the collection sites suggests that thermal heterogeneity might be one of contributors to the diversifying selection, although biological and biochemical mechanisms of the likely thermal adaptation are not known currently. Our results highlight the risk of rapid adaptation of P. infestans and possibly other pathogens as well to host resistance, and the application of eco-evolutionary principles is necessary for sustainable disease management in agricultural ecosystems.

The intramolecular allostery of GRB2 governing its interaction with SOS1 is modulated by phosphotyrosine ligands

Growth factor receptor-bound protein 2 (GRB2) is a trivalent adaptor protein and a key element in signal transduction. It interacts via its flanking nSH3 and cSH3 domains with the proline-rich domain (PRD) of the RAS activator SOS1 and via its central SH2 domain with phosphorylated tyrosine residues of receptor tyrosine kinases (RTKs; e.g. HER2). The elucidation of structural organization and mechanistic insights into GRB2 interactions, however, remain challenging due to their inherent flexibility. This study represents an important advance in our mechanistic understanding of how GRB2 links RTKs to SOS1. Accordingly, it can be proposed that (1) HER2 pYP-bound SH2 potentiates GRB2 SH3 domain interactions with SOS1 (an allosteric mechanism); (2) the SH2 domain blocks cSH3, enabling nSH3 to bind SOS1 first before cSH3 follows (an avidity-based mechanism); and (3) the allosteric behavior of cSH3 to other domains appears to be unidirectional, although there is an allosteric effect between the SH2 and SH3 domains.

Evaluating the Efficiency of the Martini Force Field to Study Protein Dimerization in Aqueous and Membrane Environments

Protein-protein complex assembly is one of the major drivers of biological response. Understanding the mechanisms of protein oligomerization/dimerization would allow one to elucidate how these complexes participate in biological activities and could ultimately lead to new approaches in designing novel therapeutic agents. However, determining the exact association pathways and structures of such complexes remains a challenge. Here, we use parallel tempering metadynamics simulations in the well-tempered ensemble to evaluate the performance of Martini 2.2P and Martini open-beta 3 (Martini 3) force fields in reproducing the structure and energetics of the dimerization process of membrane proteins and proteins in an aqueous solution in reasonable accuracy and throughput. We find that Martini 2.2P systematically overestimates the free energy of association by estimating large barriers in distinct areas, which likely leads to overaggregation when multiple monomers are present. In comparison, the less viscous Martini 3 results in a systematic underestimation of the free energy of association for proteins in solution, while it performs well in describing the association of membrane proteins. In all cases, the near-native dimer complexes are identified as minima in the free energy surface albeit not always as the lowest minima. In the case of Martini 3, we find that the spurious supramolecular protein aggregation present in Martini 2.2P multimer simulations is alleviated and thus this force field may be more suitable for the study of protein oligomerization. We propose that the use of enhanced sampling simulations with a refined coarse-grained force field and appropriately defined collective variables is a robust approach for studying the protein dimerization process, although one should be cautious of the ranking of energy minima.

An Effective MM/GBSA Protocol for Absolute Binding Free Energy Calculations: A Case Study on SARS-CoV-2 Spike Protein and the Human ACE2 Receptor

The binding free energy calculation of protein-ligand complexes is necessary for research into virus-host interactions and the relevant applications in drug discovery. However, many current computational methods of such calculations are either inefficient or inaccurate in practice. Utilizing implicit solvent models in the molecular mechanics generalized Born surface area (MM/GBSA) framework allows for efficient calculations without significant loss of accuracy. Here, GBNSR6, a new flavor of the generalized Born model, is employed in the MM/GBSA framework for measuring the binding affinity between SARS-CoV-2 spike protein and the human ACE2 receptor. A computational protocol is developed based on the widely studied Ras-Raf complex, which has similar binding free energy to SARS-CoV-2/ACE2. Two options for representing the dielectric boundary of the complexes are evaluated: one based on the standard Bondi radii and the other based on a newly developed set of atomic radii (OPT1), optimized specifically for protein-ligand binding. Predictions based on the two radii sets provide upper and lower bounds on the experimental references: -14.7(ΔGbindBondi)<-10.6(ΔGbindExp.)<-4.1(ΔGbindOPT1) kcal/mol. The consensus estimates of the two bounds show quantitative agreement with the experiment values. This work also presents a novel truncation method and computational strategies for efficient entropy calculations with normal mode analysis. Interestingly, it is observed that a significant decrease in the number of snapshots does not affect the accuracy of entropy calculation, while it does lower computation time appreciably. The proposed MM/GBSA protocol can be used to study the binding mechanism of new variants of SARS-CoV-2, as well as other relevant structures.

Biophysical and Structural Characterization of Novel RAS-Binding Domains (RBDs) of PI3Kα and PI3Kγ

Phosphatidylinositol-3-kinases (PI3Ks) are lipid kinases that phosphorylate phosphatidylinositol 4,5-bisphosphate to generate a key lipid second messenger, phosphatidylinositol 3,4,5-bisphosphate. PI3Kα and PI3Kγ require activation by RAS proteins to stimulate signaling pathways that control cellular growth, differentiation, motility and survival. Intriguingly, RAS binding to PI3K isoforms likely differ, as RAS mutations have been identified that discriminate between PI3Kα and PI3Kγ, consistent with low sequence homology (23%) between their RAS binding domains (RBDs). As disruption of the RAS/PI3Kα interaction reduces tumor growth in mice with RAS- and epidermal growth factor receptor driven skin and lung cancers, compounds that interfere with this key interaction may prove useful as anti-cancer agents. However, a structure of PI3Kα bound to RAS is lacking, limiting drug discovery efforts. Expression of full-length PI3K isoforms in insect cells has resulted in low yield and variable activity, limiting biophysical and structural studies of RAS/PI3K interactions. This led us to generate the first RBDs from PI3Kα and PI3Kγ that can be expressed at high yield in bacteria and bind to RAS with similar affinity to full-length PI3K. We also solved a 2.31 Å X-ray crystal structure of the PI3Kα-RBD, which aligns well to full-length PI3Kα. Structural differences between the PI3Kα and PI3Kγ RBDs are consistent with differences in thermal stability and may underly differential RAS recognition and RAS-mediated PI3K activation. These high expression, functional PI3K RBDs will aid in interrogating RAS interactions and could aid in identifying inhibitors of this key interaction.

MMGB/SA Consensus Estimate of the Binding Free Energy Between the Novel Coronavirus Spike Protein to the Human ACE2 Receptor

The ability to estimate protein-protein binding free energy in a computationally efficient via a physics-based approach is beneficial to research focused on the mechanism of viruses binding to their target proteins. Implicit solvation methodology may be particularly useful in the early stages of such research, as it can offer valuable insights into the binding process, quickly. Here we evaluate the potential of the related molecular mechanics generalized Born surface area (MMGB/SA) approach to estimate the binding free energy ΔGbind between the SARS-CoV-2 spike receptor-binding domain and the human ACE2 receptor. The calculations are based on a recent flavor of the generalized Born model, GBNSR6. Two estimates of ΔGbind are performed: one based on standard bondi radii, and the other based on a newly developed set of atomic radii (OPT1), optimized specifically for protein-ligand binding. We take the average of the resulting two ΔGbind values as the consensus estimate. For the well-studied Ras-Raf protein-protein complex, which has similar binding free energy to that of the SARS-CoV-2/ACE2 complex, the consensus ΔGbind = -11.8 ± 1 kcal/mol, vs. experimental -9.7 ± 0.2 kcal/mol. The consensus estimates for the SARS-CoV-2/ACE2 complex is ΔGbind = -9.4 ± 1.5 kcal/mol, which is in near quantitative agreement with experiment (-10.6 kcal/mol). The availability of a conceptually simple MMGB/SA-based protocol for analysis of the SARS-CoV-2 /ACE2 binding may be beneficial in light of the need to move forward fast.

Conformation-specific inhibitors of activated Ras GTPases reveal limited Ras dependency of patient-derived cancer organoids

The small GTPases H, K, and NRAS are molecular switches indispensable for proper regulation of cellular proliferation and growth. Several mutations in the genes encoding members of this protein family are associated with cancer and result in aberrant activation of signaling processes caused by a deregulated recruitment of downstream effector proteins. In this study, we engineered variants of the Ras-binding domain (RBD) of the C-Raf proto-oncogene, Ser/Thr kinase (CRAF). These variants bound with high affinity with the effector-binding site of Ras in an active conformation. Structural characterization disclosed how the newly identified RBD mutations cooperate and thereby enhance affinity with the effector-binding site in Ras compared with WT RBD. The engineered RBD variants closely mimicked the interaction mode of naturally occurring Ras effectors and acted as dominant-negative affinity reagents that block Ras signal transduction. Experiments with cancer cells showed that expression of these RBD variants inhibits Ras signaling, reducing cell growth and inducing apoptosis. Using these optimized RBD variants, we stratified patient-derived colorectal cancer organoids with known Ras mutational status according to their response to Ras inhibition. These results revealed that the presence of Ras mutations was insufficient to predict sensitivity to Ras inhibition, suggesting that not all of these tumors required Ras signaling for proliferation. In summary, by engineering the Ras/Raf interface of the CRAF-RBD, we identified potent and selective inhibitors of Ras in its active conformation that outcompete binding of Ras-signaling effectors.

Analysis of Ras-effector interaction competition in large intestine and colorectal cancer context

Cancer is the second leading cause of death globally, and colorectal cancer (CRC) is among the five most common cancers. The small GTPase KRAS is an oncogene that is mutated in ~30% of all CRCs. Pharmacological treatments of CRC are currently unsatisfactory, but much hope rests on network-centric approaches to drug development and cancer treatment. These approaches, however, require a better understanding of how networks downstream of Ras oncoproteins are connected in a particular tissue context – here colon and CRC. Previously we have shown that competition for binding to a ‘hub’ protein, such as Ras, can induce a rewiring of signal transduction networks. In this study, we analysed 56 established and predicted effectors that contain a structural domain with the potential ability to bind to Ras oncoproteins and their link to pathways coordinating intestinal homoeostasis and barrier function. Using protein concentrations in colon tissue and Ras-effector binding affinities, a computational network model was generated that predicted how effectors differentially and competitively bind to Ras in colon context. The model also predicted both qualitative and quantitative changes in Ras-effector complex formations with increased levels of active Ras – to simulate its upregulation in cancer – simply as an emergent property of competition for the same binding interface on the surface of Ras. We also considered how the number of Ras-effector complexes at the membrane can be increased by additional domains present in some effectors that are recruited to the membrane in response to specific conditions (inputs/stimuli/growth factors) in colon context and CRC.

A systems mechanism for KRAS mutant allele-specific responses to targeted therapy

Cancer treatment decisions are increasingly guided by which specific genes are mutated within each patient's tumor. For example, agents inhibiting the epidermal growth factor receptor (EGFR) benefit many colorectal cancer (CRC) patients, with the general exception of those whose tumor includes a KRAS mutation. However, among the various KRAS mutations, that which encodes the G13D mutant protein (KRASG13D) behaves differently; for unknown reasons, KRASG13D CRC patients benefit from the EGFR-blocking antibody cetuximab. Controversy surrounds this observation, because it contradicts the well-established mechanisms of EGFR signaling with regard to RAS mutations. Here, we identified a systems-level, mechanistic explanation for why KRASG13D cancers respond to EGFR inhibition. A computational model of RAS signaling revealed that the biophysical differences between the three most common KRAS mutants were sufficient to generate different sensitivities to EGFR inhibition. Integrated computation with experimentation then revealed a nonintuitive, mutant-specific dependency of wild-type RAS activation by EGFR that is determined by the interaction strength between KRAS and the tumor suppressor neurofibromin (NF1). KRAS mutants that strongly interacted with and competitively inhibited NF1 drove wild-type RAS activation in an EGFR-independent manner, whereas KRASG13D weakly interacted with and could not competitively inhibit NF1 and, thus, KRASG13D cells remained dependent on EGFR for wild-type RAS activity. Overall, our work demonstrates how systems approaches enable mechanism-based inference in genomic medicine and can help identify patients for selective therapeutic strategies.

The protonation states of GTP and GppNHp in Ras proteins

The small GTPase Ras transmits signals in a variety of cellular signaling pathways, most prominently in cell proliferation. GTP hydrolysis in the active center of Ras acts as a prototype for many GTPases and is the key to the understanding of several diseases, including cancer. Therefore, Ras has been the focus of intense research over the last decades. A recent neutron diffraction crystal structure of Ras indicated a protonated γ-guanylyl imidodiphosphate (γ-GppNHp) group, which has put the protonation state of GTP in question. A possible protonation of GTP was not considered in previously published mechanistic studies. To determine the detailed prehydrolysis state of Ras, we calculated infrared and NMR spectra from quantum mechanics/molecular mechanics (QM/MM) simulations and compared them with those from previous studies. Furthermore, we measured infrared spectra of GTP and several GTP analogs bound to lipidated Ras on a membrane system under near-native conditions. Our findings unify results from previous studies and indicate a structural model confirming the hypothesis that γ-GTP is fully deprotonated in the prehydrolysis state of Ras.

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.

Deconstruction of the Ras switching cycle through saturation mutagenesis

Ras proteins are highly conserved signaling molecules that exhibit regulated, nucleotide-dependent switching between active and inactive states. The high conservation of Ras requires mechanistic explanation, especially given the general mutational tolerance of proteins. Here, we use deep mutational scanning, biochemical analysis and molecular simulations to understand constraints on Ras sequence. Ras exhibits global sensitivity to mutation when regulated by a GTPase activating protein and a nucleotide exchange factor. Removing the regulators shifts the distribution of mutational effects to be largely neutral, and reveals hotspots of activating mutations in residues that restrain Ras dynamics and promote the inactive state. Evolutionary analysis, combined with structural and mutational data, argue that Ras has co-evolved with its regulators in the vertebrate lineage. Overall, our results show that sequence conservation in Ras depends strongly on the biochemical network in which it operates, providing a framework for understanding the origin of global selection pressures on proteins.

The RAS-Effector Interface: Isoform-Specific Differences in the Effector Binding Regions

RAS effectors specifically interact with the GTP-bound form of RAS in response to extracellular signals and link them to downstream signaling pathways. The molecular nature of effector interaction by RAS is well-studied but yet still incompletely understood in a comprehensive and systematic way. Here, structure-function relationships in the interaction between different RAS proteins and various effectors were investigated in detail by combining our in vitro data with in silico data. Equilibrium dissociation constants were determined for the binding of HRAS, KRAS, NRAS, RRAS1 and RRAS2 to both the RAS binding (RB) domain of CRAF and PI3Kα, and the RAS association (RA) domain of RASSF5, RALGDS and PLCε, respectively, using fluorescence polarization. An interaction matrix, constructed on the basis of available crystal structures, allowed identification of hotspots as critical determinants for RAS-effector interaction. New insights provided by this study are the dissection of the identified hotspots in five distinct regions (R1 to R5) in spite of high sequence variability not only between, but also within, RB/RA domain-containing effectors proteins. Finally, we propose that intermolecular β-sheet interaction in R1 is a central recognition region while R3 may determine specific contacts of RAS versus RRAS isoforms with effectors.

Computing Protein-Protein Association Affinity with Hybrid Steered Molecular Dynamics

Computing protein-protein association affinities is one of the fundamental challenges in computational biophysics/biochemistry. The overwhelming amount of statistics in the phase space of very high dimensions cannot be sufficiently sampled even with today's high-performance computing power. In this article, we extend a potential of mean force (PMF)-based approach, the hybrid steered molecular dynamics (hSMD) approach we developed for ligand-protein binding, to protein-protein association problems. For a protein complex consisting of two protomers, P1 and P2, we choose m (≥3) segments of P1 whose m centers of mass are to be steered in a chosen direction and n (≥3) segments of P2 whose n centers of mass are to be steered in the opposite direction. The coordinates of these m + n centers constitute a phase space of 3(m + n) dimensions (3(m + n)D). All other degrees of freedom of the proteins, ligands, solvents, and solutes are freely subject to the stochastic dynamics of the all-atom model system. Conducting SMD along a line in this phase space, we obtain the 3(m + n)D PMF difference between two chosen states: one single state in the associated state ensemble and one single state in the dissociated state ensemble. This PMF difference is the first of four contributors to the protein-protein association energy. The second contributor is the 3(m + n - 1)D partial partition in the associated state accounting for the rotations and fluctuations of the (m + n - 1) centers while fixing one of the m + n centers of the P1-P2 complex. The two other contributors are the 3(m - 1)D partial partition of P1 and the 3(n - 1)D partial partition of P2 accounting for the rotations and fluctuations of their m - 1 or n - 1 centers while fixing one of the m/n centers of P1/P2 in the dissociated state. Each of these three partial partitions can be factored exactly into a 6D partial partition in multiplication with a remaining factor accounting for the small fluctuations while fixing three of the centers of P1, P2, or the P1-P2 complex, respectively. These small fluctuations can be well-approximated as Gaussian, and every 6D partition can be reduced in an exact manner to three problems of 1D sampling, counting the rotations and fluctuations around one of the centers as being fixed. We implement this hSMD approach to the Ras-RalGDS complex, choosing three centers on RalGDS and three on Ras (m = n = 3). At a computing cost of about 71.6 wall-clock hours using 400 computing cores in parallel, we obtained the association energy, -9.2 ± 1.9 kcal/mol on the basis of CHARMM 36 parameters, which well agrees with the experimental data, -8.4 ± 0.2 kcal/mol.

Molecular determinants of drug-specific sensitivity for epidermal growth factor receptor (EGFR) exon 19 and 20 mutants in non-small cell lung cancer

We hypothesized that aberrations activating epidermal growth factor receptor (EGFR) via dimerization would be more sensitive to anti-dimerization agents (e.g., cetuximab). EGFR exon 19 abnormalities (L747_A750del; deletes amino acids LREA) respond to reversible EGFR kinase inhibitors (TKIs). Exon 20 in-frame insertions and/or duplications (codons 767 to 774) and T790M mutations are clinically resistant to reversible/some irreversible TKIs. Their impact on protein function/therapeutic actionability are not fully elucidated.In our study, the index patient with non-small cell lung cancer (NSCLC) harbored EGFR D770_P772del_insKG (exon 20). A twenty patient trial (NSCLC cohort) (cetuximab-based regimen) included two participants with EGFR TKI-resistant mutations ((i) exon 20 D770>GY; and (ii) exon 19 LREA plus exon 20 T790M mutations). Structural modeling predicted that EGFR exon 20 anomalies (D770_P772del_insKG and D770>GY), but not T790M mutations, stabilize the active dimer configuration by increasing the interaction between the kinase domains, hence sensitizing to an agent preventing dimerization. Consistent with predictions, the two patients harboring D770_P772del_insKG and D770>GY, respectively, responded to an EGFR antibody (cetuximab)-based regimen; the T790M-bearing patient showed no response to cetuximab combined with erlotinib. In silico modeling merits investigation of its ability to optimize therapeutic selection based on structural/functional implications of different aberrations within the same gene.

Understanding and Manipulating Electrostatic Fields at the Protein-Protein Interface Using Vibrational Spectroscopy and Continuum Electrostatics Calculations

Biological function emerges in large part from the interactions of biomacromolecules in the complex and dynamic environment of the living cell. For this reason, macromolecular interactions in biological systems are now a major focus of interest throughout the biochemical and biophysical communities. The affinity and specificity of macromolecular interactions are the result of both structural and electrostatic factors. Significant advances have been made in characterizing structural features of stable protein-protein interfaces through the techniques of modern structural biology, but much less is understood about how electrostatic factors promote and stabilize specific functional macromolecular interactions over all possible choices presented to a given molecule in a crowded environment. In this Feature Article, we describe how vibrational Stark effect (VSE) spectroscopy is being applied to measure electrostatic fields at protein-protein interfaces, focusing on measurements of guanosine triphosphate (GTP)-binding proteins of the Ras superfamily binding with structurally related but functionally distinct downstream effector proteins. In VSE spectroscopy, spectral shifts of a probe oscillator's energy are related directly to that probe's local electrostatic environment. By performing this experiment repeatedly throughout a protein-protein interface, an experimental map of measured electrostatic fields generated at that interface is determined. These data can be used to rationalize selective binding of similarly structured proteins in both in vitro and in vivo environments. Furthermore, these data can be used to compare to computational predictions of electrostatic fields to explore the level of simulation detail that is necessary to accurately predict our experimental findings.

Structural basis for the interaction of the adaptor protein grb14 with activated ras

Grb14, a member of the Grb7-10-14 family of cytoplasmic adaptor proteins, is a tissue-specific negative regulator of insulin signaling. Grb7-10-14 contain several signaling modules, including a Ras-associating (RA) domain, a pleckstrin-homology (PH) domain, a family-specific BPS (between PH and SH2) region, and a C-terminal Src-homology-2 (SH2) domain. We showed previously that the RA and PH domains, along with the BPS region and SH2 domain, are necessary for downregulation of insulin signaling. Here, we report the crystal structure at 2.4-Å resolution of the Grb14 RA and PH domains in complex with GTP-loaded H-Ras (G12V). The structure reveals that the Grb14 RA and PH domains form an integrated structural unit capable of binding simultaneously to small GTPases and phosphoinositide lipids. The overall mode of binding of the Grb14 RA domain to activated H-Ras is similar to that of the RA domains of RalGDS and Raf1 but with important distinctions. The integrated RA-PH structural unit in Grb7-10-14 is also found in a second adaptor family that includes Rap1-interacting adaptor molecule (RIAM) and lamellipodin, proteins involved in actin-cytoskeleton rearrangement. The structure of Grb14 RA-PH in complex with H-Ras represents the first detailed molecular characterization of tandem RA-PH domains bound to a small GTPase and provides insights into the molecular basis for specificity.

Subcellular control of Rac-GTPase signalling by magnetogenetic manipulation inside living cells

Many cell functions rely on the coordinated activity of signalling pathways at a subcellular scale. However, there are few tools capable of probing and perturbing signalling networks with a spatial resolution matching the intracellular dimensions of their activity patterns. Here we present a generic magnetogenetic approach based on the self-assembly of signalling complexes on the surface of functionalized magnetic nanoparticles inside living cells. The nanoparticles act as nanoscopic hot spots that can be displaced by magnetic forces and trigger signal transduction pathways that bring about a cell response. We applied this strategy to Rho-GTPases, a set of molecular switches known to regulate cell morphology via complex spatiotemporal patterns of activity. We demonstrate that the nanoparticle-mediated activation of signalling pathways leads to local remodelling of the actin cytoskeleton and to morphological changes.

The PaaX-type repressor MeqR2 of Arthrobacter sp. strain Rue61a, involved in the regulation of quinaldine catabolism, binds to its own promoter and to catabolic promoters and specifically responds to anthraniloyl coenzyme A

The genes coding for quinaldine catabolism in Arthrobacter sp. strain Rue61a are clustered on the linear plasmid pAL1 in two upper pathway operons (meqABC and meqDEF) coding for quinaldine conversion to anthranilate and a lower pathway operon encoding anthranilate degradation via coenzyme A (CoA) thioester intermediates. The meqR2 gene, located immediately downstream of the catabolic genes, codes for a PaaX-type transcriptional repressor. MeqR2, purified as recombinant fusion protein, forms a dimer in solution and shows specific and cooperative binding to promoter DNA in vitro. DNA fragments recognized by MeqR2 contained a highly conserved palindromic motif, 5'-TGACGNNCGTcA-3', which is located at positions -35 to -24 of the two promoters that control the upper pathway operons, at positions +4 to +15 of the promoter of the lower pathway genes and at positions +53 to +64 of the meqR2 promoter. Disruption of the palindrome abolished MeqR2 binding. The dissociation constants (K(D)) of MeqR2-DNA complexes as deduced from electrophoretic mobility shift assays were very similar for the four promoters tested (23 nM to 28 nM). Anthraniloyl-CoA was identified as the specific effector of MeqR2, which impairs MeqR2-DNA complex formation in vitro. A binding stoichiometry of one effector molecule per MeqR2 monomer and a K(D) of 22 nM were determined for the effector-protein complex by isothermal titration calorimetry (ITC). Quantitative reverse transcriptase PCR analyses suggested that MeqR2 is a potent regulator of the meqDEF operon; however, additional regulatory systems have a major impact on transcriptional control of the catabolic operons and of meqR2.

Specific Conformational States of Ras GTPase upon Effector Binding

To uncover the structural and dynamical determinants involved in the highly specific binding of Ras GTPase to its effectors, the conformational states of Ras in uncomplexed form and complexed to the downstream effectors Byr2, PI3Kγ, PLCε, and RalGDS were investigated using molecular dynamics and cross-comparison of the trajectories. The subtle changes in the dynamics and conformations of Ras upon effector binding require an analysis that targets local changes independent of global motions. Using a structural alphabet, a computational procedure is proposed to quantify local conformational changes. Positions detected by this approach were characterized as either specific for a particular effector, specific for an effector domain type, or as effector unspecific. A set of nine structurally connected residues (Ras residues 5–8, 32–35, 39–42, 55–59, 73–78, and 161–165), which link the effector binding site to the distant C-terminus, changed dynamics upon effector binding, indicating a potential effector-unspecific signaling route within the Ras structure. Additional conformational changes were detected along the N-terminus of the central β-sheet. Besides the Ras residues at the effector interface (e.g., D33, E37, D38, and Y40), which adopt effector-specific local conformations, the binding signal propagates from the interface to distant hot-spot residues, in particular to Y5 and D57. The results of this study reveal possible conformational mechanisms for the stabilization of the active state of Ras upon downstream effector binding and for the structural determinants responsible for effector specificity.

Electrostatic effects of mutations of Ras glutamine 61 measured using vibrational spectroscopy of a thiocyanate probe

Mutations of human oncoprotein p21(Ras) (hereafter Ras) at glutamine 61 are known to slow the rate of guanosine triphosphate (GTP) hydrolysis and transform healthy cells into malignant cells. It has been hypothesized that this glutamine plays a role in the intrinsic mechanism of GTP hydrolysis by interacting with an active site water molecule that electrostatically stabilizes the formation of the charged transition state at the γ-phosphate during hydrolysis. We have tested the interactions between amino acids at this position and water by measuring changes in the electrostatic field experienced by a nitrile probe positioned near Ras Q61 using vibrational Stark effect (VSE) spectroscopy. We mutated this glutamine to every amino acid except cysteine and proline and then incubated these mutants with a Ral guanine nucleotide dissociation stimulator (Ral) containing the I18C mutation that was chemically labeled with a thiocyanate vibrational spectroscopic probe. The formation of the docked Ras Q61X-labeled Ral complex was confirmed by measurement of the dissociation constant of the interaction. We measured the absorption energy of this nitrile to determine any differences in electrostatic environment in the immediate vicinity of the thiocyanate probe between wild type and mutants of Ras. For each Ras Q61X mutant, we correlate the change in electrostatic field at position 61 with the solvent accessible surface area of polar components of the mutant side chain determined from a Boltzmann-weighted ensemble of structures, as well as the residue's hydration potential. These results support the hypothesis that the role of Ras Q61 is to stabilize water in or near the active site during GTP hydrolysis. The substantial effect that nonpolar side chains of Ras Q61X have on the absorption energy of the thiocyanate must be investigated with further experiments.

Isothermal titration calorimetry - a new method for the quantification of microbial degradation of trace pollutants

The environmental fate and, in particular, biodegradation rates of hydrophobic organic compounds (HOC) are of high interest due to the ubiquity, persistence, and potential health effects of these compounds. HOC tend to interact with bioreactor materials and sampling devices and are frequently volatile, so that conventionally derived degradation parameters are often biased. We report on the development and validation of a novel calorimetric approach that serves to gain real time information on the kinetics and the physiology of HOC bioconversion in aqueous systems while overcoming weaknesses of conventional biodegradation experiments. Soil bacteria Mycobacterium frederiksbergense LB501T, Rhodococcus erythropolis K2-3 and Pseudomonas putida G7 were exposed to pulsed titrations of dissolved anthracene, 4-(2,4-dichlorophenoxy)butyric acid or naphthalene, respectively, and the thermal responses were monitored. The combinations of strains and pollutants were selected as examples for complete and partial biodegradation and complete degradation with storage product formation, respectively. Heat production signals were interpreted thermodynamically and in terms of Michaelis-Menten kinetics. The half-saturation constant k(D) and the degradation rate r(D)(Max) were derived. Comparison with conventional methods shows the suitability to extract kinetic degradation parameters of organic trace pollutants from simple ITC experiments, while thermodynamic interpretation provided further information about the metabolic fate of HOC compounds.

What makes Ras an efficient molecular switch: a computational, biophysical, and structural study of Ras-GDP interactions with mutants of Raf

Ras is a small GTP-binding protein that is an essential molecular switch for a wide variety of signaling pathways including the control of cell proliferation, cell cycle progression and apoptosis. In the GTP-bound state, Ras can interact with its effectors, triggering various signaling cascades in the cell. In the GDP-bound state, Ras looses its ability to bind to known effectors. The interaction of the GTP-bound Ras (Ras(GTP)) with its effectors has been studied intensively. However, very little is known about the much weaker interaction between the GDP-bound Ras (Ras(GDP)) and Ras effectors. We investigated the factors underlying the nucleotide-dependent differences in Ras interactions with one of its effectors, Raf kinase. Using computational protein design, we generated mutants of the Ras-binding domain of Raf kinase (Raf) that stabilize the complex with Ras(GDP). Most of our designed mutations narrow the gap between the affinity of Raf for Ras(GTP) and Ras(GDP), producing the desired shift in binding specificity towards Ras(GDP). A combination of our best designed mutation, N71R, with another mutation, A85K, yielded a Raf mutant with a 100-fold improvement in affinity towards Ras(GDP). The Raf A85K and Raf N71R/A85K mutants were used to obtain the first high-resolution structures of Ras(GDP) bound to its effector. Surprisingly, these structures reveal that the loop on Ras previously termed the switch I region in the Ras(GDP).Raf mutant complex is found in a conformation similar to that of Ras(GTP) and not Ras(GDP). Moreover, the structures indicate an increased mobility of the switch I region. This greater flexibility compared to the same loop in Ras(GTP) is likely to explain the natural low affinity of Raf and other Ras effectors to Ras(GDP). Our findings demonstrate that an accurate balance between a rigid, high-affinity conformation and conformational flexibility is required to create an efficient and stringent molecular switch.

Characterization of a Ras mutant with identical GDP- and GTP-bound structures

We previously characterized the G60A mutant of Ras and showed that the switch regions of the GTP-bound but not the GDP-bound form of this mutant adopt an "open conformation" similar to that seen in nucleotide-free Ras. Here, we mutate Lys147 of the conserved (145)SAK(147) motif in the G60A background and characterize the resulting double mutant (DM). We show that RasDM is the first structure of a Ras protein with identical GDP- and GTP-bound structures. Both structures adopt the open conformation of the active form of RasG60A. The increase in the accessible surface area of the nucleotide is consistent with a 4-fold increase in its dissociation rate. Stopped-flow experiments show no major difference in the two-step kinetics of association of GDP or GTP with the wild type, G60A, or RasDM. Addition of Sos fails to accelerate nucleotide exchange. Overexpression of the G60A or double mutant of Ras in COS-1 cells fails to activate Erk and shows a strong dominant negative effect. Our data suggest that flexibility at position 60 is required for proper Sos-catalyzed nucleotide exchange and that structural information is somehow shared among the switch regions and the different nucleotide binding motifs.

Improved binding of raf to Ras.GDP is correlated with biological activity

The GTP-binding protein Ras plays a central role in the regulation of various cellular processes, acting as a molecular switch that triggers signaling cascades. Only Ras bound to GTP is able to interact strongly with effector proteins like Raf kinase, phosphatidylinositol 3-kinase, and RalGDS, whereas in the GDP-bound state, the stability of the complex is strongly decreased, and signaling is interrupted. To determine whether this process is only controlled by the stability of the complex, we used computer-aided protein design to improve the interaction between Ras and effector. We challenged the Ras.Raf complex in this study because Raf among all effectors shows the highest Ras affinity and the fastest association kinetics. The proposed mutations were characterized as to their changes in dynamics and binding strength. We demonstrate that Ras-Raf interaction can only be improved at the cost of a loss in specificity of Ras.GTP versus Ras.GDP. As shown by NMR spectroscopy, the Raf mutation A85K leads to a shift of Ras switch I in the GTP-bound as well as in the GDP-bound state, thereby increasing the complex stability. In a luciferase-based reporter gene assay, Raf A85K is associated with higher signaling activity, which appears to be a mere matter of Ras-Raf affinity.

Thermodynamic characterization of two homologous protein complexes: associations of the semaphorin receptor plexin-B1 RhoGTPase binding domain with Rnd1 and active Rac1

Plexin receptors function in response to semaphorin guidance cues in a variety of developmental processes involving cell motility. Interactions with Rho, as well as Ras family small GTPases are critical events in the cell signaling mechanism. We have recently determined the structure of a cytoplasmic domain (RBD) of plexin-B1 and mapped its binding interface with several Rho-GTPases, Rac1, Rnd1, and RhoD. All three GTPases associate with a similar region of this plexin domain, but show different functional behavior in cells. To understand whether thermodynamic properties of the GTPase-RBD interaction contribute to such different behavior, we have examined the interaction at different temperatures, buffer, and pH conditions. Although the binding affinity of both Rnd1 and Rac1 with the plexin-B1 RBD is similar, the detailed thermodynamic properties of the interactions are considerably different. These data suggest that on Rac1 binding to the plexin-B1 RBD, the proteins become more rigid in the complex. By contrast, Rnd1 binding is consistent with unchanged or slightly increased flexibility in one or both proteins. Both GTPases show an appreciable reduction in affinity for the dimeric plexin-B1 RBD indicating that GTPase binding is not cooperative with dimer formation, but that a partial steric hindrance destabilizes the dimer. However, a reduced affinity binding mode to a disulphide stabilized model for the dimeric RBD is also possible. Consistent with cellular studies, the interaction thermodynamics imply that further levels of regulation involving additional binding partners and/or regions outside of the RhoGTPase binding domain are required for receptor activation.

A systems perspective of ras signaling in cancer

The development of cancer reflects the complex interactions and properties of many proteins functioning as part of large biochemical networks within the cancer cell. Although traditional experimental models have provided us with wonderful insights on the behavior of individual proteins within a cancer cell, they have been deficient in simultaneously keeping track of many proteins and their interactions in large networks. Computational models have emerged as a powerful tool for investigating biochemical networks due to their ability to meaningfully assimilate numerous network properties. Using the well-studied Ras oncogene as an example, we discuss the use of models to investigate pathologic Ras signaling and describe how these models could play a role in the development of new cancer drugs and the design of individualized treatment regimens.

Dissecting the thermodynamics of GAP-RhoA interactions

We describe a detailed study of the RhoA-binding epitope of the GAP domain of Graf, including the determination of the thermodynamic and kinetic parameters of the interaction of wild-type domain, and of its 15 single-site mutants, with cognate GTPases. We show that residues important for the structural integrity of the Arg-finger loop are critical for binding Rho and for the catalytic activity of GAP, but GTPase selectivity appears to be modulated by a much more subtle interplay of electrostatic and hydrophobic interactions involving residues on the periphery of the main interface. The eight residues targeted in this study are involved in three distinct patches on the surface, two of which appear to interact with highly conserved regions of the GTPase, while the third plays a role in GTPase selectivity.

Mapping the nucleotide and isoform-dependent structural and dynamical features of Ras proteins

Ras GTPases are conformational switches controlling cell proliferation, differentiation, and development. Despite their prominent role in many forms of cancer, the mechanism of conformational transition between inactive GDP-bound and active GTP-bound states remains unclear. Here we describe a detailed analysis of available experimental structures and molecular dynamics simulations to quantitatively assess the structural and dynamical features of active and inactive states and their interconversion. We demonstrate that GTP-bound and nucleotide-free G12V H-ras sample a wide region of conformational space, and show that the inactive-to-active transition is a multiphase process defined by the relative rearrangement of the two switches and the orientation of Tyr32. We also modeled and simulated N- and K-ras proteins and found that K-ras is more flexible than N- and H-ras. We identified a number of isoform-specific, long-range side chain interactions that define unique pathways of communication between the nucleotide binding site and the C terminus.

Molecular dynamics-solvated interaction energy studies of protein-protein interactions: the MP1-p14 scaffolding complex

Using the MP1-p14 scaffolding complex from the mitogen-activated protein kinase signaling pathway as model system, we explored a structure-based computational protocol to probe and characterize binding affinity hot spots at protein-protein interfaces. Hot spots are located by virtual alanine-scanning consensus predictions over three different energy functions and two different single-structure representations of the complex. Refined binding affinity predictions for select hot-spot mutations are carried out by applying first-principle methods such as the molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy (SIE) to the molecular dynamics (MD) trajectories for mutated and wild-type complexes. Here, predicted hot-spot residues were actually mutated to alanine, and crystal structures of the mutated complexes were determined. Two mutated MP1-p14 complexes were investigated, the p14(Y56A)-mutated complex and the MP1(L63A,L65A)-mutated complex. Alternative ways to generate MD ensembles for mutant complexes, not relying on crystal structures for mutated complexes, were also investigated. The SIE function, fitted on protein-ligand binding affinities, gave absolute binding affinity predictions in excellent agreement with experiment and outperformed standard MM-GBSA predictions when tested on the MD ensembles of Ras-Raf and Ras-RalGDS protein-protein complexes. For wild-type and mutant MP1-p14 complexes, SIE predictions of relative binding affinities were supported by a yeast two-hybrid assay that provided semiquantitative relative interaction strengths. Results on the MP1-mutated complex suggested that SIE predictions deteriorate if mutant MD ensembles are approximated by just mutating the wild-type MD trajectory. The SIE data on the p14-mutated complex indicated feasibility for generating mutant MD ensembles from mutated wild-type crystal structure, despite local structural differences observed upon mutation. For energetic considerations, this would circumvent costly needs to produce and crystallize mutated complexes. The sensitized protein-protein interface afforded by the p14(Y56A) mutation identified here has practical applications in screening-based discovery of first-generation small-molecule hits for further development into specific modulators of the mitogen-activated protein kinase signaling pathway.

Choose your own path: specificity in Ras GTPase signaling

The Ras superfamily of small G proteins contributes importantly to numerous cellular and physiological processes (M. F. Olsen and R. Marais, Semin. Immunol., 2000, 12, 63). This family comprises a large class of proteins (more than 150) which all share a common enzymatic function: hydrolysis of the gamma-phosphate of guanosine triphosphate (GTP) to create the products guanosine diphosphate (GDP) and inorganic phosphate (Y. Takai, T. Sasaki and T. Matozaki, Physiol. Rev., 2001, 81, 153). For this reason Ras family proteins, which include the Ras, Rho, Arf/Sara, Ran and Rab subfamilies, are classified as GTPases (G. W. Reuther and C. J. Der, Curr. Opin. Cell Biol., 2000, 12, 157). Guanine nucleotide coupling is a key regulator of enzymatic function; thus, Ras family GTPases participate in signal transduction. Ras signaling depends on binding to effectors. Many of the known effectors can bind to multiple Ras isotypes, often leading to common cellular outcomes, but each Ras isotype also engages specific effector pathways to mediate unique functions. Further, each Ras isotype can propagate multiple signaling pathways, indicating the presence of cellular determinants which allow for promiscuity in Ras-effector interactions while also maintaining specificity. Small distinctions in sequence, structure, and/or cellular regulation contribute to these differences in Ras-effector binding and subsequent cellular effects. A major focus of investigation in the Ras signaling field is identifying the determinants of these individualized functions. This review will attempt to summarize the current state of understanding of this question (with a particular focus on the Ras subfamily) and the approaches being taken to address it, and will discuss prospective areas for future investigation.

Real-time observation of temperature-dependent protein–protein interactions using real-time dual-color detection system

This study examined the ability of a real-time dual-color detection system to allow direct observations of the kinetics of temperature-dependent protein-protein interaction at a single-molecule level. The primary target protein was an Alexa Fluor 488-labeled actin conjugate, which had been pre-incubated with an unlabeled rabbit anti-actin antibody (IgG). The complementary fluorescent protein was Alexa Fluor 633-labeled goat anti-rabbit IgG antibody, which interacts with the rabbit anti-actin antibody (IgG) bound to the Alexa Fluor 488-labeled actin conjugate. The individual protein molecules labeled with different fluorescent dyes in solution were effectively focused, interacted with the other protein molecules at 500 aM, and detected directly in real-time using the dual-wavelength (lambda(ex)=488 and 635 nm) laser-induced fluorescence detection system. The kinetics of the protein-protein interactions were examined at different temperatures (12-32 degrees C). At concentrations in the aM range, the number of bound complex molecules through the protein-protein interaction decreased gradually with time at a given temperature, and increased with decreasing temperature at a set time. A high concentration (above 500 pM) of the protein sample caused aggregation and nonspecific binding of the protein molecules, even though the protein molecules were not an example of complementary binding. The results demonstrated that the real-time kinetics of a protein-protein interaction could be analyzed effectively at the single-molecule level without any time delay using the real-time dual-color detection system.

Structure of a transient intermediate for GTP hydrolysis by ras

The flexibility of the conserved 57DTAGQ61 motif is essential for Ras proper cycling in response to growth factors. Here, we increase the flexibility of the 57DTAGQ61 motif by mutating Gln61 to Gly. The crystal structure of the RasQ61G mutant reveals a new conformation of switch 2 that bears remarkable structural homology to an intermediate for GTP hydrolysis revealed by targeted molecular dynamics simulations. The mutation increased retention of GTP and inhibited Ras binding to the catalytic site, but not to the distal site of Sos. Most importantly, the thermodynamics of RafRBD binding to Ras are altered even though the structure of switch 1 is not affected by the mutation. Our results suggest that interplay and transmission of structural information between the switch regions are important factors for Ras function. They propose that initiation of GTP hydrolysis sets off the separation of the Ras/effector complex even before the GDP conformation is reached.

Crystal Structure of the Human Retinitis Pigmentosa 2 Protein and Its Interaction with Arl3

The crystal structure of human retinitis pigmentosa 2 protein (RP2) was solved to 2.1 angstroms resolution. It consists of an N-terminal beta helix and a C-terminal ferredoxin-like alpha/beta domain. RP2 is functionally and structurally related to the tubulin-specific chaperone cofactor C. Seven of nine known RP2 missense mutations identified in patients are located in the beta helix domain, and most of them cluster to the hydrophobic core and are likely to destabilize the protein. Two residues, Glu138 and the catalytically important Arg118, are solvent-exposed and form a salt bridge, indicating that Glu138 might be critical for positioning Arg118 for catalysis. RP2 is a specific effector protein of Arl3. The N-terminal 34 residues and beta helix domain of RP2 are required for this interaction. The abilitities of RP2 to bind Arl3 and cause retinitis pigmentosa seem to be correlated, since both the R118H and E138G mutants show a drastically reduced affinity to Arl3.

NMR studies of Escherichia coli acyl carrier protein: dynamic and structural differences of the apo- and holo-forms

Two indicators of conformational variability of Escherichia coli acyl carrier protein (ACP) have been investigated, namely backbone dynamics and chemical shift variations of ACP. Hydrophobic interactions between the 4'-PP prosthetic group and the hydrophobic pocket enclosed by the amphipathic helices resulted in chemical shift perturbations in the residues near the prosthetic group binding sites and contact sites in the hydrophobic pockets upon conversion from apo- to holo-forms. At pH 7.9, destabilization of ACP due to negative charge repulsions and the deprotonated state of His 75 resulted in observed chemical shift changes in the C-terminal region. Model-free analysis showed that the alpha(1)alpha(2) loop region near the prosthetic group binding site in ACP shows the greatest flexibility (lowest S(2) values) and this result may suggest these flexibilities are required for structural rearrangements when the acyl chain binds to the prosthetic group of ACP. Flexibility of ACP shown in this study is essential for its ability to interact with functionally different enzyme partners specifically and weakly in the rapid delivery of acyl chain from one partner to another.

Structure of doubly prenylated Ypt1:GDI complex and the mechanism of GDI-mediated Rab recycling

In eukaryotic cells Rab/Ypt GTPases represent a family of key membrane traffic controllers that associate with their targeted membranes via C-terminally conjugated geranylgeranyl groups. GDP dissociation inhibitor (GDI) is a general and essential regulator of Rab recycling that extracts prenylated Rab proteins from membranes at the end of their cycle of activity and facilitates their delivery to the donor membranes. Here, we present the structure of a complex between GDI and a doubly prenylated Rab protein. We show that one geranylgeranyl residue is deeply buried in a hydrophobic pocket formed by domain II of GDI, whereas the other lipid is more exposed to solvent and is skewed across several atoms of the first moiety. Based on structural information and biophysical measurements, we propose mechanistic and thermodynamic models for GDI and Rab escort protein-mediated interaction of RabGTPase with intracellular membranes.

Nucleotide binding and filament assembly of recombinant yeast septin complexes

Septins are filament-forming GTPases involved in cytokinesis and cortical organization. In the yeast Saccharomyces cerevisiae, the septins encoded by CDC3, CDC10, CDC11, and CDC12 form a high-molecular-weight complex, localized at the cytoplasmic face of the plasma membrane in the mother-bud neck. While septin function at the cellular level is fairly well understood, progress on structure-function analysis of these proteins has been slow and limited by the lack of large amounts of pure complex. While monomeric septins form apparently non-native aggregates, stable recombinant complexes of two, three, or four yeast septins can be produced by co-expression from bi-cistronic vectors in E. coli. The septin polypeptides show various degrees of saturation with guanine nucleotides in different complexes. The binary core Cdc3p-Cdc12p complex contains no bound nucleotide. While ternary complexes are partially saturated and can bind extraneously added nucleotide with micromolar affinity, only the complete four-component septin complex is fully coordinated with tightly bound GDP/GTP after chromatographic purification. We show here that the nucleotide-binding sites of the septins show drastic changes on formation of higher oligomers. Although the binary core Cdc3p-Cdc12p complex does not form filaments, the ternary and quaternary complexes form bundles of paired filaments. In the case of ternary complexes, filament formation is stimulated by guanine nucleotide, but is not dependent on the presence or absence of the gamma-phosphate.

Recognizing and defining true Ras binding domains I: biochemical analysis

Common domain databases contain sequence motifs which belong to the ubiquitin fold family and are called Ras binding (RB) and Ras association (RalGDS/AF6 Ras associating) (RA) domains. The name implies that they bind to Ras (or Ras-like) GTP-binding proteins, and a few of them have been documented to qualify as true Ras effectors, defined as binding only to the activated GTP-bound form of Ras. Here we have expressed a large number of these domains and investigated their interaction with Ras, Rap and M-Ras. While their (albeit weak) sequence homology suggest that the domains adopt a common fold, not all of them bind to Ras proteins, irrespective of whether they are called RB or RA domains. We used fluorescence spectroscopy and isothermal titration calorimetry to show that the binding affinities vary over a large range, and are usually specific for either Ras or Rap. Moreover, the specificity is dictated by a set of key residues in the interface. Stopped-flow kinetic analysis showed that the association rate constants determine the different affinities of effector binding, while the dissociation rate constants are in a similar range. Manual sequence analysis allowed us to define positively charged sequence epitopes in certain secondary structure elements of the ubiquitin fold (beta1, beta2 and alpha1) which are located at similar positions and comprise the hot spots of the binding interface. These residues are important to qualify an RA/RB domain as a true candidate Ras or Rap effector.

Recognizing and Defining True Ras Binding Domains II: In Silico Prediction Based on Homology Modelling and Energy Calculations

Considering the large number of putative Ras effector proteins, it is highly desirable to develop computational methods to be able to identify true Ras binding molecules. Based on a limited sequence homology among members of the Ras association (RA) and Ras binding (RB) sub-domain families of the ubiquitin super-family, we have built structural homology models of Ras proteins in complex with different RA and RB domains, using the FOLD-X software. A critical step in our approach is to use different templates of Ras complexes, in order to account for the structural variation among the RA and RB domains. The homology models are validated by predicting the effect of mutating hot spot residues in the interface, and residues important for the specificity of interaction with different Ras proteins. The FOLD-X calculated energies of the best-modelled complexes are in good agreement with previously published experimental data and with new data reported here. Based on these results, we can establish energy thresholds above, or below which, we can predict with 96% confidence that a RA/RB domain will or will not interact with Ras. This study shows the importance of in depth structural analysis, high quality force-fields and modelling for correct prediction. Our work opens the possibility of genome-wide prediction for this protein family and for others, where there is enough structural information.

The nucleotide switch in Cdc42 modulates coupling between the GTPase-binding and allosteric equilibria of Wiskott-Aldrich syndrome protein

The GTP/GDP nucleotide switch in Ras superfamily GTPases generally involves differential affinity toward downstream effectors, with the GTP-bound state having a higher affinity for effector than the GDP-bound state. We have developed a quantitative model of allosteric regulation of the Wiskott-Aldrich syndrome protein (WASP) by the Rho GTPase Cdc42 to better understand how GTPase binding is coupled to effector activation. The model accurately predicts WASP affinity for Cdc42, activity toward Arp2/3 complex, and activation by Cdc42 as functions of a two-state allosteric equilibrium in WASP. The ratio of GTPase affinities for the inactive and active states of WASP is appreciably larger for Cdc42-GTP than for Cdc42-GDP. The greater ability to distinguish between the two states of WASP makes Cdc42-GTP a full WASP agonist, whereas Cdc42-GDP is only a partial agonist. Thus, the nucleotide switch controls not only the affinity of Cdc42 for its effector but also the efficiency of coupling between the Cdc42-binding and allosteric equilibria in WASP. This effect can ensure high fidelity and specificity in Cdc42 signaling in crowded membrane environments.

Structural Basis of Rho GTPase-Mediated Activation of the Formin mDia1

Diaphanous-related formins (DRFs) regulate dynamics of unbranched actin filaments during cell contraction and cytokinesis. DRFs are autoinhibited through intramolecular binding of a Diaphanous autoinhibitory domain (DAD) to a conserved N-terminal regulatory element. Autoinhibition is relieved through binding of the GTPase RhoA to the N-terminal element. We report the crystal structure of the dimeric regulatory domain of the DRF, mDia1. Dimerization is mediated by an intertwined six-helix bundle, from which extend two Diaphanous inhibitory domains (DIDs) composed of five armadillo repeats. NMR and biochemical mapping indicate the RhoA and DAD binding sites on the DID partially overlap, explaining activation of mDia1 by the GTPase. RhoA binding also requires an additional structurally independent segment adjacent to the DID. This regulatory construction, involving a GTPase binding site spanning a flexibly tethered arm and the inhibitory module, is observed in many autoinhibited effectors of Ras superfamily GTPases, suggesting evolutionary pressure for this design.