1
|
Lowegard AU, Frenkel MS, Holt GT, Jou JD, Ojewole AA, Donald BR. Novel, provable algorithms for efficient ensemble-based computational protein design and their application to the redesign of the c-Raf-RBD:KRas protein-protein interface. PLoS Comput Biol 2020; 16:e1007447. [PMID: 32511232 PMCID: PMC7329130 DOI: 10.1371/journal.pcbi.1007447] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2019] [Revised: 07/01/2020] [Accepted: 05/13/2020] [Indexed: 11/25/2022] Open
Abstract
The K* algorithm provably approximates partition functions for a set of states (e.g., protein, ligand, and protein-ligand complex) to a user-specified accuracy ε. Often, reaching an ε-approximation for a particular set of partition functions takes a prohibitive amount of time and space. To alleviate some of this cost, we introduce two new algorithms into the osprey suite for protein design: fries, a Fast Removal of Inadequately Energied Sequences, and EWAK*, an Energy Window Approximation to K*. fries pre-processes the sequence space to limit a design to only the most stable, energetically favorable sequence possibilities. EWAK* then takes this pruned sequence space as input and, using a user-specified energy window, calculates K* scores using the lowest energy conformations. We expect fries/EWAK* to be most useful in cases where there are many unstable sequences in the design sequence space and when users are satisfied with enumerating the low-energy ensemble of conformations. In combination, these algorithms provably retain calculational accuracy while limiting the input sequence space and the conformations included in each partition function calculation to only the most energetically favorable, effectively reducing runtime while still enriching for desirable sequences. This combined approach led to significant speed-ups compared to the previous state-of-the-art multi-sequence algorithm, BBK*, while maintaining its efficiency and accuracy, which we show across 40 different protein systems and a total of 2,826 protein design problems. Additionally, as a proof of concept, we used these new algorithms to redesign the protein-protein interface (PPI) of the c-Raf-RBD:KRas complex. The Ras-binding domain of the protein kinase c-Raf (c-Raf-RBD) is the tightest known binder of KRas, a protein implicated in difficult-to-treat cancers. fries/EWAK* accurately retrospectively predicted the effect of 41 different sets of mutations in the PPI of the c-Raf-RBD:KRas complex. Notably, these mutations include mutations whose effect had previously been incorrectly predicted using other computational methods. Next, we used fries/EWAK* for prospective design and discovered a novel point mutation that improves binding of c-Raf-RBD to KRas in its active, GTP-bound state (KRasGTP). We combined this new mutation with two previously reported mutations (which were highly-ranked by osprey) to create a new variant of c-Raf-RBD, c-Raf-RBD(RKY). fries/EWAK* in osprey computationally predicted that this new variant binds even more tightly than the previous best-binding variant, c-Raf-RBD(RK). We measured the binding affinity of c-Raf-RBD(RKY) using a bio-layer interferometry (BLI) assay, and found that this new variant exhibits single-digit nanomolar affinity for KRasGTP, confirming the computational predictions made with fries/EWAK*. This new variant binds roughly five times more tightly than the previous best known binder and roughly 36 times more tightly than the design starting point (wild-type c-Raf-RBD). This study steps through the advancement and development of computational protein design by presenting theory, new algorithms, accurate retrospective designs, new prospective designs, and biochemical validation. Computational structure-based protein design is an innovative tool for redesigning proteins to introduce a particular or novel function. One such function is improving the binding of one protein to another, which can increase our understanding of important protein systems. Herein we introduce two novel, provable algorithms, fries and EWAK*, for more efficient computational structure-based protein design as well as their application to the redesign of the c-Raf-RBD:KRas protein-protein interface. These new algorithms speed-up computational structure-based protein design while maintaining accurate calculations, allowing for larger, previously infeasible protein designs. Additionally, using fries and EWAK* within the osprey suite, we designed the tightest known binder of KRas, a heavily studied cancer target that interacts with a number of different proteins. This previously undiscovered variant of a KRas-binding domain, c-Raf-RBD, has potential to serve as a tool to further probe the protein-protein interface of KRas with its effectors and its discovery alone emphasizes the potential for more successful applications of computational structure-based protein design.
Collapse
Affiliation(s)
- Anna U. Lowegard
- Program in Computational Biology and Bioinformatics, Duke University Medical Center, Durham, North Carolina, United States of America
- Department of Computer Science, Duke University, Durham, North Carolina, United States of America
| | - Marcel S. Frenkel
- Department of Biochemistry, Duke University Medical Center, Durham, North Carolina, United States of America
| | - Graham T. Holt
- Program in Computational Biology and Bioinformatics, Duke University Medical Center, Durham, North Carolina, United States of America
- Department of Computer Science, Duke University, Durham, North Carolina, United States of America
| | - Jonathan D. Jou
- Department of Computer Science, Duke University, Durham, North Carolina, United States of America
| | - Adegoke A. Ojewole
- Program in Computational Biology and Bioinformatics, Duke University Medical Center, Durham, North Carolina, United States of America
- Department of Computer Science, Duke University, Durham, North Carolina, United States of America
| | - Bruce R. Donald
- Department of Computer Science, Duke University, Durham, North Carolina, United States of America
- Department of Biochemistry, Duke University Medical Center, Durham, North Carolina, United States of America
- * E-mail:
| |
Collapse
|
2
|
Filchtinski D, Sharabi O, Rüppel A, Vetter IR, Herrmann C, Shifman JM. What makes Ras an efficient molecular switch: a computational, biophysical, and structural study of Ras-GDP interactions with mutants of Raf. J Mol Biol 2010; 399:422-35. [PMID: 20361980 DOI: 10.1016/j.jmb.2010.03.046] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2009] [Revised: 03/19/2010] [Accepted: 03/23/2010] [Indexed: 11/16/2022]
Abstract
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.
Collapse
Affiliation(s)
- Daniel Filchtinski
- Physikalische Chemie I, Fakultät für Chemie und Biochemie, Ruhr-Universität-Bochum, Universitätstr. 150, 44780 Bochum, Germany
| | | | | | | | | | | |
Collapse
|
3
|
Fischer A, Hekman M, Kuhlmann J, Rubio I, Wiese S, Rapp UR. B- and C-RAF display essential differences in their binding to Ras: the isotype-specific N terminus of B-RAF facilitates Ras binding. J Biol Chem 2007; 282:26503-16. [PMID: 17635919 DOI: 10.1074/jbc.m607458200] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Recruitment of RAF kinases to the plasma membrane was initially proposed to be mediated by Ras proteins via interaction with the RAF Ras binding domain (RBD). Data reporting that RAF kinases possess high affinities for particular membrane lipids support a new model in which Ras-RAF interactions may be spatially restricted to the plane of the membrane. Although the coupling features of Ras binding to the isolated RAF RBD were investigated in great detail, little is known about the interactions of the processed Ras with the functional and full-length RAF kinases. Here we present a quantitative analysis of the binding properties of farnesylated and nonfarnesylated H-Ras to both full-length B- and C-RAF in the presence and absence of lipid environment. Although isolated RBD fragments associate with high affinity to both farnesylated and nonfarnesylated H-Ras, the full-length RAF kinases revealed fundamental differences with respect to Ras binding. In contrast to C-RAF that requires farnesylated H-Ras, cytosolic B-RAF associates effectively and with significantly higher affinity with both farnesylated and nonfarnesylated H-Ras. To investigate the potential farnesyl binding site(s) we prepared several N-terminal fragments of C-RAF and found that in the presence of cysteine-rich domain only the farnesylated form of H-Ras binds with high association rates. The extreme N terminus of B-RAF turned out to be responsible for the facilitation of lipid independent Ras binding to B-RAF, since truncation of this region resulted in a protein that changed its kinase properties and resembles C-RAF. In vivo studies using PC12 and COS7 cells support in vitro results. Co-localization measurements using labeled Ras and RAF documented essential differences between B- and C-RAF with respect to association with Ras. Taken together, these data suggest that the activation of B-RAF, in contrast to C-RAF, may take place both at the plasma membrane and in the cytosolic environment.
Collapse
Affiliation(s)
- Andreas Fischer
- Institut für Medizinische Strahlenkunde und Zellforschung, University of Wuerzburg, 97078 Wuerzburg, Germany
| | | | | | | | | | | |
Collapse
|
4
|
Campbell-Valois FX, Tarassov K, Michnick SW. Massive sequence perturbation of a small protein. Proc Natl Acad Sci U S A 2005; 102:14988-93. [PMID: 16214894 PMCID: PMC1257691 DOI: 10.1073/pnas.0500465102] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Most protein topologies rarely occur in nature, thus limiting our ability to extract sequence information that could be used to predict structure, function, and evolutionary constraints on protein folds. In principle, the sequence diversity explored by a given protein topology could be expanded by introducing sequence perturbations and selecting variant proteins that fold correctly. However, our capacity to explore sequence space is intrinsically limited by the enormous number of sequences generated from the 20 amino acids and the limited number of variants likely to fold. Here we sought to test whether the sequence space for naturally existing proteins can be explored by simple, sequential degeneration of a complete set of short sequence segments of a model protein, without long-range covariation. Using the Raf ras binding domain as a model of a small protein capable of autonomous folding, we degenerated 72 of 76 positions of the primary structure for the 20 amino acids in segments of four to seven residues defined by secondary structure and selected the folded species for interaction with h-ras by using an in vivo survival-selection assay. The methodology presented allowed for rigorous statistical analysis and comparison of sequence diversity. The ensemble of sequence variants of Raf ras binding domain obtained have recaptured the diversity observed for the ubiquitin-roll topology. A signature sequence for this fold and the implication of this strategy to protein design and structure prediction are discussed.
Collapse
Affiliation(s)
- F-X Campbell-Valois
- Département de Biochimie, Université de Montréal, C.P. 6128, Succ. Centre-ville, Montréal, QC, Canada H3C 3J7
| | | | | |
Collapse
|
5
|
Abstract
Ras effectors have convergently developed a common subdomain in their otherwise unrelated protein body for their interaction with Ras. Structural analysis revealed that the mode of interaction is highly similar for all Ras effectors, but is completely different from that of effectors of other subfamilies of small GTPases. Whereas the molecular mechanism of effector activation is still elusive, detailed knowledge about the thermodynamics and dynamics of the interaction with Ras has accumulated.
Collapse
Affiliation(s)
- Christian Herrmann
- Max-Planck-Institute for Molecular Physiology, Otto-Hahn-Strasse 11, 44227 Dortmund, Germany.
| |
Collapse
|
6
|
Quilliam LA, Rebhun JF, Castro AF. A growing family of guanine nucleotide exchange factors is responsible for activation of Ras-family GTPases. PROGRESS IN NUCLEIC ACID RESEARCH AND MOLECULAR BIOLOGY 2003; 71:391-444. [PMID: 12102558 DOI: 10.1016/s0079-6603(02)71047-7] [Citation(s) in RCA: 222] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
GTPases of the Ras subfamily regulate a diverse array of cellular-signaling pathways, coupling extracellular signals to the intracellular response machinery. Guanine nucleotide exchange factors (GEFs) are primarily responsible for linking cell-surface receptors to Ras protein activation. They do this by catalyzing the dissociation of GDP from the inactive Ras proteins. GTP can then bind and induce a conformational change that permits interaction with downstream effectors. Over the past 5 years, approximately 20 novel Ras-family GEFs have been identified and characterized. These data indicate that a variety of different signaling mechanisms can be induced to activate Ras, enabling tyrosine kinases, G-protein-coupled receptors, adhesion molecules, second messengers, and various protein-interaction modules to relocate and/or activate GEFs and elevate intracellular Ras-GTP levels. This review discusses the structure and function of the catalytic or CDC25 homology domain common to almost all Ras-family GEFs. It also details our current knowledge about the regulation and function of this rapidly growing family of enzymes that include Sos1 and 2, GRF1 and 2, CalDAG-GEF/GRP1-4, C3G, cAMP-GEF/Epac 1 and 2, PDZ-GEFs, MR-GEF, RalGDS family members, RalGPS, BCAR3, Smg GDS, and phospholipase C(epsilon).
Collapse
Affiliation(s)
- Lawrence A Quilliam
- Department of Biochemistry and Molecular, Biology and Walther Oncology Center, Indiana University School of Medicine, Indianapolis 46202, USA
| | | | | |
Collapse
|
7
|
Catimel B, Rothacker J, Nice E. The use of biosensors for microaffinity purification: an integrated approach to proteomics. JOURNAL OF BIOCHEMICAL AND BIOPHYSICAL METHODS 2001; 49:289-312. [PMID: 11694286 DOI: 10.1016/s0165-022x(01)00205-6] [Citation(s) in RCA: 19] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- B Catimel
- Ludwig Institute for Cancer Research, Melbourne Tumour Biology Branch, The Royal Melbourne Hospital, Post Office Box 2008, Victoria 3050, Australia
| | | | | |
Collapse
|
8
|
Catimel B, Teh T, Fontes MR, Jennings IG, Jans DA, Howlett GJ, Nice EC, Kobe B. Biophysical characterization of interactions involving importin-alpha during nuclear import. J Biol Chem 2001; 276:34189-98. [PMID: 11448961 DOI: 10.1074/jbc.m103531200] [Citation(s) in RCA: 134] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Proteins containing the classical nuclear localization sequences (NLSs) are imported into the nucleus by the importin-alpha/beta heterodimer. Importin-alpha contains the NLS binding site, whereas importin-beta mediates the translocation through the nuclear pore. We characterized the interactions involving importin-alpha during nuclear import using a combination of biophysical techniques (biosensor, crystallography, sedimentation equilibrium, electrophoresis, and circular dichroism). Importin-alpha is shown to exist in a monomeric autoinhibited state (association with NLSs undetectable by biosensor). Association with importin-beta (stoichiometry, 1:1; K(D) = 1.1 x 10(-8) m) increases the affinity for NLSs; the importin-alpha/beta complex binds representative monopartite NLS (simian virus 40 large T-antigen) and bipartite NLS (nucleoplasmin) with affinities (K(D) = 3.5 x 10(-8) m and 4.8 x 10(-8) m, respectively) comparable with those of a truncated importin-alpha lacking the autoinhibitory domain (T-antigen NLS, K(D) = 1.7 x 10(-8) m; nucleoplasmin NLS, K(D) = 1.4 x 10(-8) m). The autoinhibitory domain (as a separate peptide) binds the truncated importin-alpha, and the crystal structure of the complex resembles the structure of full-length importin-alpha. Our results support the model of regulation of nuclear import mediated by the intrasteric autoregulatory sequence of importin-alpha and provide a quantitative description of the binding and regulatory steps during nuclear import.
Collapse
Affiliation(s)
- B Catimel
- Ludwig Institute for Cancer Research, Royal Melbourne Hospital, Victoria 3050, Australia
| | | | | | | | | | | | | | | |
Collapse
|
9
|
Abstract
We have compiled a comprehensive list of the articles published in the year 2000 that describe work employing commercial optical biosensors. Selected reviews of interest for the general biosensor user are highlighted. Emerging applications in areas of drug discovery, clinical support, food and environment monitoring, and cell membrane biology are emphasized. In addition, the experimental design and data processing steps necessary to achieve high-quality biosensor data are described and examples of well-performed kinetic analysis are provided.
Collapse
Affiliation(s)
- R L Rich
- Center for Biomolecular Interaction Analysis, University of Utah, Salt Lake City, UT 84132, USA
| | | |
Collapse
|
10
|
Use of thiazolidine-mediated ligation for site specific biotinylation of mouse EGF for biosensor immobilisation. ACTA ACUST UNITED AC 2001. [DOI: 10.1007/bf02446519] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|