Light on the structural communication in Ras GTPases

Structural communication between the GTPase Sec4p and its activator Sec2p: Determinants of GEF activity and early deformations to nucleotide release

Ras GTPases are molecular switches that cycle between OFF and ON states depending on the bound nucleotide (i.e. GDP-bound and GTP-bound, respectively).

The Rab GTPase, Sec4p, plays regulatory roles in multiple steps of intracellular vesicle trafficking. Nucleotide release is catalyzed by the Guanine Nucleotide Exchange Factor (GEF) Sec2p.

Here, the integration of structural information with molecular dynamics (MD) simulations addressed a number of questions concerning the intrinsic and stimulated dynamics of Sec2p and Sec4p as well as the chain of structural deformations leading to GEF-assisted activation of the Rab GTPase.

Sec2p holds an intrinsic ability to adopt the conformation found in the crystallographic complexes with Sec4p, thus suggesting that the latter selects and shifts the conformational equilibrium towards a pre-existing bound-like conformation of Sec2p.

The anchoring of Sec4p to a suitable conformation of Sec2p favors the Sec2p-assisted pulling on itself of the α1/switch 1 (SWI) loop and of SWI, which loose any contact with GDP. Those deformations of Sec4p would occur earlier. Formation of the final Sec2p-Sec4p hydrophobic interface, accomplishes later. Disruption of the nucleotide cage would cause firstly loss of interactions with the guanine ring and secondly loss of interactions with the phosphates.

The ease in sampling the energy landscape and adopting a bound-like conformation likely favors the catalyzing ability of GEFs for Ras GTPases.

PSNtools for standalone and web-based structure network analyses of conformational ensembles

Structure graphs, in which interacting amino acids/nucleotides correspond to linked nodes, represent cutting-edge tools to investigate macromolecular function.

The graph-based approach defined as Protein Structure Network (PSN) was initially implemented in the Wordom software and subsequently in the webPSN server. PSNs are computed either on a molecular dynamics (MD) trajectory (PSN-MD) or on a single structure. In the latter case, information on atomic fluctuations is inferred from the Elastic Network Model-Normal Mode Analysis (ENM-NMA) (PSN-ENM). While Wordom performs both PSN-ENM and PSN-MD analyses but without output post-processing, the webPSN server performs only single-structure PSN-EMN but assisting the user in input setup and output analysis.

Here we release for the first time the standalone software PSNtools, which allows calculation and post-processing of PSN analyses carried out either on single structures or on conformational ensembles. Relevant unique and novel features of PSNtools are either comparisons of two networks or computations of consensus networks on sets of homologous/analogous macromolecular structures or conformational ensembles. Network comparisons and consensus serve to infer differences in functionally different states of the same system or network-based signatures in groups of bio-macromolecules sharing either the same functionality or the same fold.

In addition to the new software, here we release also an updated version of the webPSN server, which allows performing an interactive graphical analysis of PSN-MD, following the upload of the PSNtools output.

PSNtools, the auxiliary binary version of Wordom software, and the WebPSN server are freely available at http://webpsn.hpc.unimo.it/wpsn3.php.

Dynamics and structural communication in the ternary complex of fully phosphorylated V2 vasopressin receptor, vasopressin, and β-arrestin 1

G protein-coupled receptors (GPCRs) are critically regulated by arrestins, which not only desensitize G-protein signaling but also initiate a G protein-independent wave of signaling. The information from structure determination was herein exploited to build a structural model of the ternary complex, comprising fully phosphorylated V2 vasopressin receptor (V2R), the agonist arginine vasopressin (AVP), and β-arrestin 1 (β-arr1). Molecular simulations served to explore dynamics and structural communication in the ternary complex. Flexibility and mechanical profiles reflect fold of V2R and β-arr1. Highly conserved amino acids tend to behave as hubs in the structure network and contribute the most to the mechanical rigidity of V2R seven-helix bundle and of β-arr1. Two structurally and dynamically distinct receptor-arrestin interfaces assist the twist of the N- and C-terminal domains (ND and CD, respectively) of β-arr1 with respect to each other, which is linked to arrestin activation. While motion of the ND is essentially assisted by the fully phosphorylated C-tail of V2R (V2RCt), that of CD is assisted by the second and third intracellular loops and the cytosolic extensions of helices 5 and 6. In the presence of the receptor, the β-arr1 inter-domain twist angle correlates with the modes describing the essential subspace of the ternary complex. β-arr1 motions are also influenced by the anchoring to the membrane of the C-edge-loops in the β-arr1-CD. Overall fluctuations reveal a coupling between motions of the agonist binding site and of β-arr1-ND, which are in allosteric communication between each other. Mechanical rigidity points, often acting as hubs in the structure network and distributed along the main axis of the receptor helix bundle, contribute to establish a preferential communication pathway between agonist ligand and the ND of arrestin. Such communication, mediated by highly conserved amino acids, involves also the first amino acid in the arrestin C-tail, which is highly dynamic and is involved in clathrin-mediated GPCR internalization.

PRECOG: PREdicting COupling probabilities of G-protein coupled receptors

G-protein coupled receptors (GPCRs) control multiple physiological states by transducing a multitude of extracellular stimuli into the cell via coupling to intra-cellular heterotrimeric G-proteins. Deciphering which G-proteins couple to each of the hundreds of GPCRs present in a typical eukaryotic organism is therefore critical to understand signalling. Here, we present PRECOG (precog.russelllab.org): a web-server for predicting GPCR coupling, which allows users to: (i) predict coupling probabilities for GPCRs to individual G-proteins instead of subfamilies; (ii) visually inspect the protein sequence and structural features that are responsible for a particular coupling; (iii) suggest mutations to rationally design artificial GPCRs with new coupling properties based on predetermined coupling features.

Interconnecting Flexibility, Structural Communication, and Function in RhoGEF Oncoproteins

Dbl family Rho guanine nucleotide exchange factors (RhoGEFs) play a central role in cell biology by catalyzing the exchange of guanosine 5'-triphosphate for guanosine 5'-diphosphate (GDP) on RhoA. Insights into the oncogenic constitutive activity of the Lbc RhoGEF were gained by analyzing the structure and dynamics of the protein in different functional states and in comparison with a close homologue, leukemia-associated RhoGEF. Higher intrinsic flexibility, less dense and extended structure network, and less stable allosteric communication pathways in Lbc, compared to the nonconstitutively active homologue, emerged as major determinants of the constitutive activity. Independent of the state, the essential dynamics of the two RhoGEFs is contributed by the last 10 amino acids of Dbl homology (DH) and the whole pleckstrin homology (PH) domains and tends to be equalized by the presence of RhoA. The catalytic activity of the RhoGEF relies on the scaffolding action of the DH domain that primarily turns the switch I (SWI) of RhoA on itself through highly conserved amino acids participating in the stability core and essential for function. Changes in the conformation of SWI and disorganization of the RhoA regions deputed to nucleotide binding are among the major RhoGEF effects leading to GDP release. Binding of RhoA reorganizes the allosteric communication on RhoGEF, strengthening the communication among the canonical RhoA binding site on DH, a secondary RhoA binding site on PH, and the binding site for heterotrimeric G proteins, suggesting dual roles for RhoA as a catalysis substrate and as a regulatory protein. The structure network-based analysis tool employed in this study proved to be useful for predicting potentially druggable regulatory sites in protein structures.

Illuminating G-Protein-Coupling Selectivity of GPCRs

Heterotrimetic G proteins consist of four subfamilies (G_s, G_i/o, G_q/11, and G_12/13) that mediate signaling via G-protein-coupled receptors (GPCRs), principally by receptors binding Gα C termini. G-protein-coupling profiles govern GPCR-induced cellular responses, yet receptor sequence selectivity determinants remain elusive. Here, we systematically quantified ligand-induced interactions between 148 GPCRs and all 11 unique Gα subunit C termini. For each receptor, we probed chimeric Gα subunit activation via a transforming growth factor-α (TGF-α) shedding response in HEK293 cells lacking endogenous G_q/11 and G_12/13 proteins, and complemented G-protein-coupling profiles through a NanoBiT-G-protein dissociation assay. Interrogation of the dataset identified sequence-based coupling specificity features, inside and outside the transmembrane domain, which we used to develop a coupling predictor that outperforms previous methods. We used the predictor to engineer designer GPCRs selectively coupled to G₁₂. This dataset of fine-tuned signaling mechanisms for diverse GPCRs is a valuable resource for research in GPCR signaling.

Uncovering GPCR and G Protein Function by Protein Structure Network Analysis

Protein structure network (PSN) analysis is one of the graph theory-based approaches currently used for investigating structural communication in biomolecular systems. Information on the system's dynamics can be provided by atomistic molecular dynamics (MD) simulations or coarse grained elastic network models paired with normal mode analysis (ENM-NMA). This chapter reports on selected applications of PSN analysis to uncover the structural communication in G protein coupled receptors (GPCRs) and G proteins. Strategies to highlight changes in structural communication caused by mutations, ligand and protein binding are described. Conserved amino acids, sites of misfolding mutations, or ligands acting as functional switches tend to behave as hubs in the native structure networks. Densely linked regions in the protein structure graphs could be identified as playing central roles in protein stability and function. Changes in the communication pathway fingerprints depending on the bound ligand or following amino acid mutation could be highlighted as well. A bridge between misfolding and misrouting could be established in rhodopsin mutants linked to inherited blindness. The analysis of native network perturbations by misfolding mutations served to infer key structural elements of protein responsiveness to small chaperones with implications for drug discovery.

Structure network analysis to gain insights into GPCR function

G protein coupled receptors (GPCRs) are allosteric proteins whose functioning fundamentals are the communication between the two poles of the helix bundle. Protein structure network (PSN) analysis is one of the graph theory-based approaches currently used to investigate the structural communication in biomolecular systems. Information on system's dynamics can be provided by atomistic molecular dynamics (MD) simulations or coarse grained elastic network models paired with normal mode analysis (ENM–NMA). The present review article describes the application of PSN analysis to uncover the structural communication in G protein coupled receptors (GPCRs). Strategies to highlight changes in structural communication upon misfolding, dimerization and activation are described. Focus is put on the ENM–NMA-based strategy applied to the crystallographic structures of rhodopsin in its inactive (dark) and signalling active (meta II (MII)) states, highlighting changes in structure network and centrality of the retinal chromophore in differentiating the inactive and active states of the receptor.

Protein-Protein Interface and Disease: Perspective from Biomolecular Networks

Protein-protein interactions are involved in many important biological processes and molecular mechanisms of disease association. Structural studies of interfacial residues in protein complexes provide information on protein-protein interactions. Characterizing protein-protein interfaces, including binding sites and allosteric changes, thus pose an imminent challenge. With special focus on protein complexes, approaches based on network theory are proposed to meet this challenge. In this review we pay attention to protein-protein interfaces from the perspective of biomolecular networks and their roles in disease. We first describe the different roles of protein complexes in disease through several structural aspects of interfaces. We then discuss some recent advances in predicting hot spots and communication pathway analysis in terms of amino acid networks. Finally, we highlight possible future aspects of this area with respect to both methodology development and applications for disease treatment.

Catching Functional Modes and Structural Communication in Dbl Family Rho Guanine Nucleotide Exchange Factors

Computational approaches such as Principal Component Analysis (PCA) and Elastic Network Model-Normal Mode Analysis (ENM-NMA) are proving to be of great value in investigating relevant biological problems linked to slow motions with no demand in computer power. In this study, these approaches have been coupled to the graph theory-based Protein Structure Network (PSN) analysis to dissect functional dynamics and structural communication in the Dbl family of Rho Guanine Nucleotide Exchange Factors (RhoGEFs). They are multidomain proteins whose common structural feature is a DH-PH tandem domain deputed to the GEF activity that makes them play a central role in cell and cancer biology. While their common GEF action is accomplished by the DH domain, their regulatory mechanisms are highly variegate and depend on the PH and the additional domains as well as on interacting proteins. Major evolutionary-driven deformations as inferred from PCA concern the α6 helix of DH that dictates the orientation of the PH domain. Such deformations seem to depend on the mechanisms adopted by the GEF to prevent Rho binding, i.e. functional specialization linked to autoinhibition. In line with PCA, ENM-NMA indicates α6 and the linked PH domain as the portions of the tandem domain holding almost the totality of intrinsic and functional dynamics, with the α6/β1 junction acting as a hinge point for the collective motions of PH. In contrast, the DH domain holds a static scaffolding and hub behavior, with structural communication playing a central role in the regulatory actions by other domains/proteins. Possible allosteric communication pathways involving essentially DH were indeed found in those RhoGEFs acting as effectors of small or heterotrimeric RasGTPases. The employed methodology is suitable for deciphering structure/dynamics relationships in large sets of homologous or analogous proteins.

The Ras-Membrane Interface: Isoform-specific Differences in The Catalytic Domain

The small GTPase Ras is mutated in about 20% of human cancers, primarily at active site amino acid residues G12, G13, and Q61. Thus, structural biology research has focused on the active site, impairment of GTP hydrolysis by oncogenic mutants, and characterization of protein-protein interactions in the effector lobe half of the protein. The C-terminal hypervariable region has increasingly gained attention due to its importance in H-Ras, N-Ras, and K-Ras differences in membrane association. A high-resolution molecular view of the Ras-membrane interaction involving the allosteric lobe of the catalytic domain has lagged behind, although evidence suggests that it contributes to isoform specificity. The allosteric lobe has recently gained interest for harboring potential sites for more selective targeting of this elusive "undruggable" protein. The present review reveals critical insight that isoform-specific differences appear prominently at these potentially targetable sites and integrates these differences with knowledge of Ras plasma membrane localization, with the intent to better understand the structure-function relationships needed to design isoform-specific Ras inhibitors.

Network and atomistic simulations unveil the structural determinants of mutations linked to retinal diseases

A number of incurable retinal diseases causing vision impairments derive from alterations in visual phototransduction. Unraveling the structural determinants of even monogenic retinal diseases would require network-centered approaches combined with atomistic simulations. The transducin G38D mutant associated with the Nougaret Congenital Night Blindness (NCNB) was thoroughly investigated by both mathematical modeling of visual phototransduction and atomistic simulations on the major targets of the mutational effect. Mathematical modeling, in line with electrophysiological recordings, indicates reduction of phosphodiesterase 6 (PDE) recognition and activation as the main determinants of the pathological phenotype. Sub-microsecond molecular dynamics (MD) simulations coupled with Functional Mode Analysis improve the resolution of information, showing that such impairment is likely due to disruption of the PDEγ binding cavity in transducin. Protein Structure Network analyses additionally suggest that the observed slight reduction of theRGS9-catalyzed GTPase activity of transducin depends on perturbed communication between RGS9 and GTP binding site. These findings provide insights into the structural fundamentals of abnormal functioning of visual phototransduction caused by a missense mutation in one component of the signaling network. This combination of network-centered modeling with atomistic simulations represents a paradigm for future studies aimed at thoroughly deciphering the structural determinants of genetic retinal diseases. Analogous approaches are suitable to unveil the mechanism of information transfer in any signaling network either in physiological or pathological conditions.

Lessons from computer simulations of Ras proteins in solution and in membrane

BACKGROUND

A great deal has been learned over the last several decades about the function of Ras proteins in solution and membrane environments. While much of this knowledge has been derived from a plethora of experimental techniques, computer simulations have also played a substantial role.

SCOPE OF REVIEW

Our goal here is to summarize the contribution of molecular simulations to our current understanding of normal and aberrant Ras function. We focus on lessons from molecular dynamics simulations in aqueous and membrane environments.

MAJOR CONCLUSIONS

The central message is that a close interaction between theory and simulation on the one hand and cell-biological, spectroscopic and other experimental approaches on the other has played, and will likely continue to play, a vital role in Ras research.

GENERAL SIGNIFICANCE

Atomistic insights emerging from detailed simulations of Ras in solution and in bilayers may be the key to unlock the secret that to date prevented development of selective anti-Ras inhibitors for cancer therapy.

A Mixed Protein Structure Network and Elastic Network Model Approach to Predict the Structural Communication in Biomolecular Systems: The PDZ2 Domain from Tyrosine Phosphatase 1E As a Case Study

Graph theory is being increasingly used to study the structural communication in biomolecular systems. This requires incorporating information on the system's dynamics, which is time-consuming and not suitable for high-throughput investigations. We propose a mixed Protein Structure Network (PSN) and Elastic Network Model (ENM)-based strategy, i.e., PSN-ENM, for fast investigation of allosterism in biological systems. PSN analysis and ENM-Normal Mode Analysis (ENM-NMA) are implemented in the structural analysis software Wordom, freely available at http://wordom.sourceforge.net/ . The method performs a systematic search of the shortest communication pathways that traverse a protein structure. A number of strategies to compare the structure networks of a protein in different functional states and to get a global picture of communication pathways are presented as well. The approach was validated on the PDZ2 domain from tyrosine phosphatase 1E (PTP1E) in its free (APO) and peptide-bound states. PDZ domains are, indeed, the systems whose structural communication and allosteric features are best characterized both in vitro and in silico. The agreement between predictions by the PSN-ENM method and in vitro evidence is remarkable and comparable to or higher than that reached by more time-consuming computational approaches tested on the same biological system. Finally, the PSN-ENM method was able to reproduce the salient communication features of unbound and bound PTP1E inferred from molecular dynamics simulations. High speed makes this method suitable for high throughput investigation of the communication pathways in large sets of biomolecular systems in different functional states.

Modeling the structural communication in supramolecular complexes involving GPCRs

This article describes a computational strategy aimed at studying the structural communication in G-Protein Coupled Receptors (GPCRs) and G proteins. The strategy relies on comparative Molecular Dynamics (MD) simulations and analyses of wild-type (i.e., reference state) vs. mutated (i.e., perturbed state), or free (i.e., reference state) vs. bound (i.e., perturbed state) forms of a GPCR or a G protein. Bound forms of a GPCR include complexes with small ligands and/or receptor dimers/oligomers, whereas bound forms of heterotrimeric GDP-bound G proteins concern the complex with a GPCR. The computational strategy includes structure prediction of a receptor monomer (in the absence of high-resolution structure), a receptor dimer/oligomer, and a receptor-G protein complex, which constitute the inputs of MD simulations. Finally, the analyses of the MD trajectories are instrumental in inferring the structural/dynamics differences between reference and perturbed states of a GPCR or a G protein. In this respect, focus will be put on the analysis of protein structure networks and communication paths.

Network analysis to uncover the structural communication in GPCRs

Protein structure network (PSN) analysis is one of the graph theory-based approaches currently used to investigate the structural communication in biomolecular systems. Information on system dynamics can be provided by atomistic molecular dynamics simulations or coarse-grained Elastic Network Models paired with Normal Mode Analysis (ENM-NMA). This chapter describes the application of PSN analysis to uncover the structural communication in G protein-coupled receptors (GPCRs). Strategies to highlight changes in structural communication upon misfolding mutations, dimerization, and activation are described. Focus is put on the ENM-NMA-based strategy applied to the crystallographic structures of rhodopsin in its inactive (dark) and signaling active (meta II (MII)) states, highlighting clear changes in the PSN and the centrality of the retinal chromophore in differentiating the inactive and active states of the receptor.