A Structural Model for the DEAD Box Helicase YxiN in Solution: Localization of the RNA Binding Domain

DEAD box proteins consist of a common helicase core formed by two globular RecA domains that are separated by a cleft. The helicase core acts as a nucleotide-dependent switch that alternates between open and closed conformations during the catalytic cycle of duplex separation, thereby providing basic helicase activity. Flanking domains can direct the helicase core to a specific RNA substrate by mediating high-affinity or high-specificity RNA binding. In addition, they may position RNA for the helicase core or may directly contribute to unwinding. While structures of different helicase cores have been determined previously, little is known about the orientation of flanking domains relative to the helicase core.YxiN is a DEAD box protein that consists of a helicase core and a C-terminal RNA binding domain (RBD) that mediates specific binding to hairpin 92 in 23S rRNA. To provide a framework for understanding the functional cooperation of the YxiN helicase core and the RBD, we mapped the orientation of the RBD in single-molecule fluorescence resonance energy transfer experiments. We present a model for the global conformation of YxiN in which the RBD lies above a slightly concave patch that is formed by flexible loops on the surface of the C-terminal RecA domain. The orientation of the RBD is different from the orientations of flanking domains in the Thermus thermophilus DEAD box protein Hera and in Saccharomyces cerevisiae Mss116p, in line with the different functions of these DEAD box proteins and of their RBDs. Interestingly, the corresponding patch on the C-terminal RecA domain that is covered by the YxiN RBD is also part of the interface between the translation factors eIF4A and eIF4G. Possibly, this region constitutes an adaptable interface that generally allows for the interaction of the helicase core with additional domains or interacting factors.     

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 Ras‐binding domain region of RGS14 regulates its functional interactions with heterotrimeric G proteins

RGS14 is a 60 kDa protein that contains a regulator of G protein signaling (RGS) domain near its N‐terminus, a central region containing a pair of tandem Ras‐binding domains (RBD), and a GPSM (G protein signaling modulator) domain (a.k.a. Gi/o‐Loco binding [GoLoco] motif) near its C‐terminus. The RGS domain of RGS14 exhibits GTPase accelerating protein (GAP) activity toward Gαi/o proteins, while its GPSM domain acts as a guanine nucleotide dissociation inhibitor (GDI) on Gαi1 and Gαi3. In the current study, we investigate the contribution of different domains of RGS14 to its biochemical functions. Here we show that the full‐length protein has a greater GTPase activating activity but a weaker inhibition of nucleotide dissociation relative to its isolated RGS and GPSM regions, respectively. Our data suggest that these differences may be attributable to an inter‐domain interaction within RGS14 that promotes the activity of the RGS domain, but simultaneously inhibits the activity of the GPSM domain. The RBD region seems to play an essential role in this regulatory activity. Moreover, this region of RGS14 is also able to bind to members of the B/R4 subfamily of RGS proteins and enhance their effects on GPCR‐activated Gi/o proteins. Overall, our results suggest a mechanism wherein the RBD region associates with the RGS domain region, producing an intramolecular interaction within RGS14 that enhances the GTPase activating function of its RGS domain while disfavoring the negative effect of its GPSM domain on nucleotide dissociation. J. Cell. Biochem. 114: 1414–1423, 2013. © 2012 Wiley Periodicals, Inc.     

Sam68 associates with the SH3 domains of Grb2 recruiting GAP to the Grb2-SOS complex in insulin receptor signaling

The 68 kDa Src substrate associated during mitosis (Sam68) is an RNA binding protein with Src homology (SH) 2 and 3 domain binding sites. We have recently found that Sam68 is a substrate of the insulin receptor (IR) that translocates from the nucleus to the cytoplasm and that Tyr-phosphorylated Sam68 associates with the SH2 domains of p85 PI3K and GAP, in vivo and in vitro. In the present work, we have further demonstrated the cytoplasmic localization of Sam68, which is increased in cells overexpressing IR. Besides, we sought to further study the association of Sam68 with the Ras-GAP pathway by assessing the interactions with SH3 domains of Grb2. We employed GST-fusion proteins containing the SH3 domains of Grb2 (N or C), and recombinant Sam68 for in vitro studies. In vivo studies of protein-protein interaction were assessed by co-immunoprecipitation experiments with specific antibodies against Sam68, GAP, Grb2, SOS, and phosphotyrosine; and by affinity precipitation with the fusion proteins (SH3-Grb2). Insulin stimulation of HTC-IR cells promotes phosphorylation of Sam68 and its association with the SH2 domains of GAP. Sam68 is constitutively associated with the SH3 domains of Grb2 and it does not change upon insulin stimulation, but Sam68 is Tyr-phosphorylated and promotes the association of GAP with the Grb2-SOS complex. In vitro studies with fusion proteins showed that Sam68 association with Grb2 is preferentially mediated by the C-terminal SH3 domains of Grb2. In conclusion, Sam68 is a substrate of the IR and may have a role as a docking protein in IR signaling, recruiting GAP to the Grb2-SOS complex, and in this way it may modulate Ras activity.     

Mapping of the auto-inhibitory interactions of protein kinase R by nuclear magnetic resonance

The dsRNA-dependent protein kinase (PKR) is a key mediator of the anti-viral and anti-proliferative effects of interferon. Unphosphorylated PKR is characterized by inhibitory interactions between the kinase and RNA binding domains (RBDs), but the structural details of the latent state and its unraveling during activation are not well understood. To study PKR regulation by NMR we assigned a large portion of the backbone resonances of the catalytically inactive K296R kinase domain, and performed (15)N-heteronuclear single quantum coherence (HSQC) titrations of this kinase domain with the RBDs. Chemical shift perturbations in the kinase indicate that RBD2 binds to the substrate eIF2alpha docking site in the kinase C-lobe. Consistent with these results, a mutation in the eIF2alpha docking site, F495A, displays weaker interactions with the RBD. The full-length RBD1+2 binds more strongly to the kinase domain than RBD2 alone. The observed chemical shift changes extend from the eIF2alpha binding site into the kinase N-lobe and inside the active site, consistent with weak interactions between the N-terminal part of the RBD and the kinase.     

Structural insights into the role of the WW2 domain on tandem WW–PPxY motif interactions of oxidoreductase WWOX

Class I WW domains are present in many proteins of various functions and mediate protein interactions by binding to short linear PPxY motifs. Tandem WW domains often bind peptides with multiple PPxY motifs, but the interplay of WW–peptide interactions is not always intuitive. The WW domain–containing oxidoreductase (WWOX) harbors two WW domains: an unstable WW1 capable of PPxY binding and stable WW2 that cannot bind PPxY. The WW2 domain has been suggested to act as a WW1 domain chaperone, but the underlying mechanism of its chaperone activity remains to be revealed. Here, we combined NMR, isothermal calorimetry, and structural modeling to elucidate the roles of both WW domains in WWOX binding to its PPxY-containing substrate ErbB4. Using NMR, we identified an interaction surface between these two domains that supports a WWOX conformation compatible with peptide substrate binding. Isothermal calorimetry and NMR measurements also indicated that while binding affinity to a single PPxY motif is marginally increased in the presence of WW2, affinity to a dual-motif peptide increases 10-fold. Furthermore, we found WW2 can directly bind double-motif peptides using its canonical binding site. Finally, differential binding of peptides in mutagenesis experiments was consistent with a parallel N- to C-terminal PPxY tandem motif orientation in binding to the WW1–WW2 tandem domain, validating structural models of the interaction. Taken together, our results reveal the complex nature of tandem WW-domain organization and substrate binding, highlighting the contribution of WWOX WW2 to both protein stability and target binding.     

Contribution of the tyrosines to the structure and function of the human U1A N-terminal RNA binding domain.

RNA binding domains (RBDs) are members of a large family of proteins that share minimal sequence conservation but adopt an alpha beta sandwich global fold. Defining the contributions of specific amino acids to RBD structure and RNA binding is critical to understanding the functions of these proteins. In these experiments with the human U1A N-terminal RNA binding domain (RBD1), the contributions from each of its four tyrosines to protein structure, stability, and RNA binding were measured. Each tyrosine was substituted with phenylalanine and one other selected residue, and the resulting proteins were characterized by chemical denaturation to measure their unfolding free energy, by binding free energies to the wild-type RNA hairpin, and by 19F NMR to probe for structural changes. Features of the protein identified in these experiments include a possible tyrosine/lysine contact in an alpha-helix, which may be an example of an energetically favorable aromatic/amino side chain interaction. One long loop of the protein, which shows unusual 15N backbone and tyrosine side-chain dynamics, is implicated in protein:protein association. The diverse interactions of the four tyrosine residues in the organization of RBD1 illustrate how each member of this family of proteins will have unique molecular details that contribute to function.     

Finding biologically relevant protein domain interactions: conserved binding mode analysis

Proteins evolved through the shuffling of functional domains, and therefore, the same domain can be found in different proteins and species. Interactions between such conserved domains often involve specific, well-determined binding surfaces reflecting their important biological role in a cell. To find biologically relevant interactions we developed a method of systematically comparing and classifying protein domain interactions from the structural data. As a result, a set of conserved binding modes (CBMs) was created using the atomic detail of structure alignment data and the protein domain classification of the Conserved Domain Database. A conserved binding mode is inferred when different members of interacting domain families dock in the same way, such that their structural complexes superimpose well. Such domain interactions with recurring structural themes have greater significance to be biologically relevant, unlike spurious crystal packing interactions. Consequently, this study gives lower and upper bounds on the number of different types of interacting domain pairs in the structure database on the order of 1000-2000. We use CBMs to create domain interaction networks, which highlight functionally significant connections by avoiding many infrequent links between highly connected nodes. The CBMs also constitute a library of docking templates that may be used in molecular modeling to infer the characteristics of an unknown binding surface, just as conserved domains may be used to infer the structure of an unknown protein. The method's ability to sort through and classify large numbers of putative interacting domain pairs is demonstrated on the oligomeric interactions of globins.     

Solution structure of the Ras binding domain of the protein kinase Byr2 from Schizosaccharomyces pombe.

BACKGROUND: After activation, small GTPases such as Ras transfer the incoming signal to effectors by specifically interacting with the binding domain of these proteins. Structural details of the binding domain of different effectors determine which pathway is predominantly activated. Byr2 from fission yeast is a functional homolog of Raf, which is the direct downstream target of Ras in mammalians that initiates a protein kinase cascade. The amino acid sequence of Byr2's Ras binding domain is only weakly related to that of Raf, and Byr2's three-dimensional structure is unknown. RESULTS: We have solved the 3D structure of the Ras binding domain of Byr2 (Byr2RBD) from Schizosaccharomyces pombe in solution. The structure consists of three alpha helices and a mixed five-stranded beta pleated sheet arranged in the topology betabetaalphabetabetaalphabetaalpha with the first seven canonic secondary structure elements forming a ubiquitin superfold. 15N-(1)H-TROSY-HSQC spectroscopy of the complex of Byr2RBD with Ras*Mg(2+)*GppNHp reveals that the first and second beta strands and the first alpha helix of Byr2 are mainly involved in the protein-protein interaction as observed in other Ras binding domains. Although the putative interaction site of H-Ras from human and Ras1 from S. pombe are identical in sequence, binding to Byr2 leads to small but significant differences in the NMR spectra, indicating a slightly different binding mode. CONCLUSIONS: The ubiquitin superfold appears to be the general structural motif for Ras binding domains even in cases with vanishing sequence identity. However, details of the 3D structure and the interacting interface are different, thereby determining the specifity of the recognition of Ras and Ras-related proteins.     

gamma-Syntrophin scaffolding is spatially and functionally distinct from that of the alpha/beta syntrophins

The syntrophins are a family of scaffolding proteins with multiple protein interaction domains that link signaling proteins to dystrophin family members. Each of the three most characterized syntrophins (alpha, beta1, beta2) contains a PDZ domain that binds a unique set of signaling proteins including kinases, ion and water channels, and neuronal nitric oxide synthase (nNOS). The PDZ domains of the gamma-syntrophins do not bind nNOS. In vitro pull-down assays show that the gamma-syntrophins can bind dystrophin but have unique preferences for the syntrophin binding sites of dystrophin family members. Despite their ability to bind dystrophin in vitro, neither gamma-syntrophin isoform co-localizes with dystrophin in skeletal muscle. Furthermore, gamma-syntrophins do not co-purify with dystrophin isolated from mouse tissue. These data suggest that the interaction of gamma-syntrophin with dystrophin is transient and potentially subject to regulatory mechanisms. gamma1-Syntrophin is highly expressed in brain and is specifically localized in hippocampal pyramidal neurons, Purkinje neurons in cerebellum, and cortical neurons. gamma2-Syntrophin is expressed in many tissues including skeletal muscle where it is found only in the subsynaptic space beneath the neuromuscular junction. In both neurons and muscle, gamma-syntrophin isoforms localize to the endoplasmic reticulum where they may form a scaffold for signaling and trafficking.     

The Ras/p120 GTPase-activating protein (GAP) interaction is regulated by the p120 GAP pleckstrin homology domain

Pleckstrin homology domains are structurally conserved functional domains that can undergo both protein/protein and protein/lipid interactions. Pleckstrin homology domains can mediate inter- and intra-molecular binding events to regulate enzyme activity. They occur in numerous proteins including many that interact with Ras superfamily members, such as p120 GAP. The pleckstrin homology domain of p120 GAP is located in the NH(2)-terminal, noncatalytic region of p120 GAP. Overexpression of the noncatalytic domains of p120 GAP may modulate Ras signal transduction pathways. Here, we demonstrate that expression of the isolated pleckstrin homology domain of p120 GAP specifically inhibits Ras-mediated signaling and transformation but not normal cellular growth. Furthermore, we show that the pleckstrin homology domain binds the catalytic domain of p120 GAP and interferes with the Ras/GAP interaction. Thus, we suggest that the pleckstrin homology domain of p120 GAP may specifically regulate the interaction of Ras with p120 GAP via competitive intra-molecular binding.     

Conservation and divergence of Grb7 family of Ras-binding domains

Ras proteins are signal-transducing GTPases that cycle between inactive GDP-bound and active GTP-bound forms. Ras is a prolific signaling molecule interacting with a spectrum of effector molecules and acting through more than one signaling pathway. The Ras-effector proteins contain a Ras-associating (RA) domain through which these associate with Ras in a GTP-dependent manner. The RA domain is highly conserved among the members of the growth factor receptor-bound (Grb) 7 family of proteins which includes Grb7, Grb10 and Grb14. Our laboratory has reported an unusual observation that RA domain of Grb14 binds to the C-terminal nucleotide binding site of cyclic nucleotide gated channel (CTR-CNGA1) and inhibits the channel activity. Molecular modeling of the CTR-CNGA1 displays 50%–70% tertiary structural similarity towards Ras proteins. We named this region as Ras-like domain (RLD). The interaction between RA-Grb14 and RLD-CNGA1 is mediated through a simple protein-protein interaction temporally and spatially regulated by light and cGMP. It is interesting to note that Grb14 binds to GTPase-mutant Rab5, a Ras-related small GTPase whereas Grb10 binds only to GTP-bound form of active Rab5 but not to GTPase-defective mutant Rab5. These results suggest that Grb14 might have been evolved later in the evolution that binds to both Ras and nucleotide binding proteins such as CNGA1. Our studies also suggest that eukaryotic CNG channels could be evolved through a gene fusion between prokaryotic ion channels and cyclic nucleotide binding proteins, both of which might have undergone several sequence variations for functional adaptation during evolution.     

Distinct dynamics and interaction patterns in H‐ and K‐Ras oncogenic P‐loop mutants

Despite years of study, the structural or dynamical basis for the differential reactivity and oncogenicity of Ras isoforms and mutants remains unclear. In this study, we investigated the effects of amino acid variations on the structure and dynamics of wild type and oncogenic mutants G12D, G12V, and G13D of H‐ and K‐Ras proteins. Based on data from µs‐scale molecular dynamics simulations, we show that the overall structure of the proteins remains similar but there are important differences in dynamics and interaction networks. We identified differences in residue interaction patterns around the canonical switch and distal loop regions, and persistent sodium ion binding near the GTP particularly in the G13D mutants. Our results also suggest that different Ras variants have distinct local structural features and interactions with the GTP, variations that have the potential to affect GTP release and hydrolysis. Furthermore, we found that H‐Ras proteins and particularly the G12V and G13D variants are significantly more flexible than their K‐Ras counterparts. Finally, while most of the simulated proteins sampled the effector‐interacting state 2 conformational state, G12V and G13D H‐Ras adopted an open switch state 1 conformation that is defective in effector interaction. These differences have implications for Ras GTPase activity, effector or exchange factor binding, dimerization and membrane interaction. Proteins 2017; 85:1618–1632. © 2017 Wiley Periodicals, Inc.     

Unraveling protein interaction networks with near-optimal efficiency

The functional characterization of genes and their gene products is the main challenge of the genomic era. Examining interaction information for every gene product is a direct way to assemble the jigsaw puzzle of proteins into a functional map. Here we demonstrate a method in which the information gained from pull-down experiments, in which single proteins act as baits to detect interactions with other proteins, is maximized by using a network-based strategy to select the baits. Because of the scale-free distribution of protein interaction networks, we were able to obtain fast coverage by focusing on highly connected nodes (hubs) first. Unfortunately, locating hubs requires prior global information about the network one is trying to unravel. Here, we present an optimized 'pay-as-you-go' strategy that identifies highly connected nodes using only local information that is collected as successive pull-down experiments are performed. Using this strategy, we estimate that 90% of the human interactome can be covered by 10,000 pull-down experiments, with 50% of the interactions confirmed by reciprocal pull-down experiments.     

WW and SH3 domains, two different scaffolds to recognize proline-rich ligands

WW domains are small protein modules composed of approximately 40 amino acids. These domains fold as a stable, triple stranded beta-sheet and recognize proline-containing ligands. WW domains are found in many different signaling and structural proteins, often localized in the cytoplasm as well as in the cell nucleus. Based on analyses of seven structures of WW domains, we discuss their diverse binding preferences and sequence conservation patterns. While modeling WW domains for which structures have not been determined we uncovered a case of potential molecular and functional convergence between WW and SH3 domains. The binding surface of the modeled WW domain of Npw38 protein shows a remarkable similarity to the SH3 domain of Sem5 protein, confirming biochemical data on similar binding predilections of both domains.     

The ubiquitin domain superfold: structure-based sequence alignments and characterization of binding epitopes

Ubiquitin-like domains are present, apart from ubiquitin-like proteins themselves, in many multidomain proteins involved in different signal transduction processes. The sequence conservation for all ubiquitin superfold family members is rather poor, even between subfamily members, leading to mistakes in sequence alignments using conventional sequence alignment methods. However, a correct alignment is essential, especially for in silico methods that predict binding partners on the basis of sequence and structure. In this study, using 3D-structural information we have generated and manually corrected sequence alignments for proteins of the five ubiquitin superfold subfamilies. On the basis of this alignment, we suggest domains for which structural information will be useful to allow homology modelling. In addition, we have analysed the energetic and electrostatic properties of ubiquitin-like domains in complex with various functional binding proteins using the protein design algorithm FoldX. On the basis of an in silico alanine-scanning mutagenesis, we provide a detailed binding epitope mapping of the hotspots of the ubiquitin domain fold, involved in the interaction with different domains and proteins. Finally, we provide a consensus fingerprint sequence that identifies all sequences described to belong to the ubiquitin superfold family. It is possible that the method that we describe may be applied to other domain families sharing a similar fold but having low levels of sequence homology.     

The Ras-Byr2RBD complex: structural basis for Ras effector recognition in yeast.

BACKGROUND: The small GTP binding protein Ras has important roles in cellular growth and differentiation. Mutant Ras is permanently active and contributes to cancer development. In its activated form, Ras interacts with effector proteins, frequently initiating a kinase cascade. In the lower eukaryotic Schizosaccharomyces pombe, Byr2 kinase represents a Ras target that in terms of signal-transduction hierarchy can be considered a homolog of mammalian Raf-kinase. The activation mechanism of protein kinases by Ras is not understood, and there is no detailed structural information about Ras binding domains (RBDs) in nonmammalian organisms. RESULTS: The crystal structure of the Ras-Byr2RBD complex at 3 A resolution shows a complex architecture similar to that observed in mammalian homologous systems, with an interprotein beta sheet stabilized by predominantly polar interactions between the interacting components. The C-terminal half of the Ras switch I region contains most of the contact anchors, while on the Byr2 side, a number of residues from topologically distinct regions are involved in complex stabilization. A C-terminal helical segment, which is not present in the known mammalian homologous systems and which is part of the auto-inhibitory region, has an additional binding site outside the switch I region. CONCLUSIONS: The structure of the Ras-Byr2 complex confirms the Ras binding module as a communication element mediating Ras-effector interactions; the Ras-Byr2 complex is also conserved in a lower eukaryotic system like yeast, which is in contrast to other small GTPase families. The extra helical segment might be involved in kinase activation.     

Prediction of Ras-effector interactions using position energy matrices

MOTIVATION: One of the more challenging problems in biology is to determine the cellular protein interaction network. Progress has been made to predict protein-protein interactions based on structural information, assuming that structural similar proteins interact in a similar way. In a previous publication, we have determined a genome-wide Ras-effector interaction network based on homology models, with a high accuracy of predicting binding and non-binding domains. However, for a prediction on a genome-wide scale, homology modelling is a time-consuming process. Therefore, we here successfully developed a faster method using position energy matrices, where based on different Ras-effector X-ray template structures, all amino acids in the effector binding domain are sequentially mutated to all other amino acid residues and the effect on binding energy is calculated. Those pre-calculated matrices can then be used to score for binding any Ras or effector sequences. RESULTS: Based on position energy matrices, the sequences of putative Ras-binding domains can be scanned quickly to calculate an energy sum value. By calibrating energy sum values using quantitative experimental binding data, thresholds can be defined and thus non-binding domains can be excluded quickly. Sequences which have energy sum values above this threshold are considered to be potential binding domains, and could be further analysed using homology modelling. This prediction method could be applied to other protein families sharing conserved interaction types, in order to determine in a fast way large scale cellular protein interaction networks. Thus, it could have an important impact on future in silico structural genomics approaches, in particular with regard to increasing structural proteomics efforts, aiming to determine all possible domain folds and interaction types. AVAILABILITY: All matrices are deposited in the ADAN database (http://adan-embl.ibmc.umh.es/). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.