Massive sequence perturbation of the Raf ras binding domain reveals relationships between sequence conservation, secondary structure propensity, hydrophobic core organization and stability

Insertion state of modular protein nanopores into a membrane

In the past decade, significant progress has been made in the development of new protein nanopores. Despite these advancements, there is a pressing need for the creation of nanopores equipped with relatively large functional groups for the sampling of biomolecular events on their extramembranous side. Here, we designed, produced, and analyzed protein nanopores encompassing a robust truncation of a monomeric β-barrel membrane protein. An exogenous stably folded protein was anchored within the aqueous phase via a flexible peptide tether of varying length. We have extensively examined the pore-forming properties of these modular protein nanopores using protein engineering and single-molecule electrophysiology. This study revealed distinctions in the nanopore conductance and current fluctuations that arose from tethering the exogenous protein to either the N terminus or the C terminus. Remarkably, these nanopores insert into a planar lipid membrane with one specific conductance among a set of three substate conductance values. Moreover, we demonstrate that the occurrence probabilities of these insertion substates depend on the length of the peptide tether, the orientation of the exogenous protein with respect to the nanopore opening, and the molecular mass of tethered protein. In addition, the three conductance values remain unaltered by major changes in the composition of modular nanopores. The outcomes of this work serve as a platform for further developments in areas of protein engineering of transmembrane pores and biosensor technology.

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

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 (KRas^GTP). 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 KRas^GTP, 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.

A Phage-Assisted Continuous Selection Approach for Deep Mutational Scanning of Protein-Protein Interactions

Protein-protein interactions (PPIs) are critical for organizing molecules in a cell and mediating signaling pathways. Dysregulation of PPIs is often a key driver of disease. To better understand the biophysical basis of such disease processes-and to potentially target them-it is critical to understand the molecular determinants of PPIs. Deep mutational scanning (DMS) facilitates the acquisition of large amounts of biochemical data by coupling selection with high throughput sequencing (HTS). The challenging and labor-intensive design and optimization of a relevant selection platform for DMS, however, limits the use of powerful directed evolution and selection approaches. To address this limitation, we designed a versatile new phage-assisted continuous selection (PACS) system using our previously reported proximity-dependent split RNA polymerase (RNAP) biosensors, with the aim of greatly simplifying and streamlining the design of a new selection platform for PPIs. After characterization and validation using the model KRAS/RAF PPI, we generated a library of RAF variants and subjected them to PACS and DMS. Our HTS data revealed positions along the binding interface that are both tolerant and intolerant to mutations, as well as which substitutions are tolerated at each position. Critically, the "functional scores" obtained from enrichment data through continuous selection for individual variants correlated with K_D values measured in vitro, indicating that biochemical data can be extrapolated from sequencing using our new system. Due to the plug and play nature of RNAP biosensors, this method can likely be extended to a variety of other PPIs. More broadly, this, and other methods under development support the continued development of evolutionary and high-throughput approaches to address biochemical problems, moving toward a more comprehensive understanding of sequence-function relationships in proteins.

Crystal structure of an ancient protein: evolution by conformational epistasis

The structural mechanisms by which proteins have evolved new functions are known only indirectly. We report x-ray crystal structures of a resurrected ancestral protein-the approximately 450 million-year-old precursor of vertebrate glucocorticoid (GR) and mineralocorticoid (MR) receptors. Using structural, phylogenetic, and functional analysis, we identify the specific set of historical mutations that recapitulate the evolution of GR's hormone specificity from an MR-like ancestor. These substitutions repositioned crucial residues to create new receptor-ligand and intraprotein contacts. Strong epistatic interactions occur because one substitution changes the conformational position of another site. "Permissive" mutations-substitutions of no immediate consequence, which stabilize specific elements of the protein and allow it to tolerate subsequent function-switching changes-played a major role in determining GR's evolutionary trajectory.

Universal strategies in research and drug discovery based on protein-fragment complementation assays

Changes in the interactions among proteins that participate in a biochemical pathway can reflect the immediate regulatory responses to intrinsic or extrinsic perturbations of the pathway. Thus, methods that allow for the direct detection of the dynamics of protein-protein interactions can be used to probe the effects of any perturbation on any pathway of interest. Here we describe experimental strategies - based on protein-fragment complementation assays (PCAs) - that can achieve this. PCA-based strategies can be used with or instead of traditional target-based drug discovery strategies to identify novel pathway-component proteins of therapeutic interest, to increase the quantity and quality of information about the actions of potential drugs, and to gain insight into the intricate networks that make up the molecular machinery of living cells.

The transition state of the ras binding domain of Raf is structurally polarized based on Phi-values but is energetically diffuse

The ras binding domain (RBD) of the Ser/Thr kinase c-Raf/Raf-1 spans 78 residues and adopts a structure characteristic of the beta-grasp ubiquitin-like topology. Recently, the primary sequence of Raf RBD has been nearly exhaustively mutated experimentally by insertion of stretches of degenerate codons, which revealed sequence conservation and hydrophobic core organization similar to that found in an alignment of beta-grasp ubiquitin-like proteins. These results now allow us to examine the relationship between sequence conservation and the folding process, particularly viewed through the analysis of transition state (TS) structure. Specifically, we present herein a protein engineering study combining classic truncation (Ala/Gly) and atypical mutants to predict folding TS ensemble properties. Based on classical Phi-value analysis, Raf RBD TS structure is particularly polarized around the N-terminal beta-hairpin. However, all residues constituting the inner layer of the hydrophobic core are involved in TS stabilization, although they are clearly found in a less native-like environment. The TS structure can also be probed by a direct measure of its destabilization upon mutation, DeltaDeltaG(U-++). Viewed through this analysis, Raf RBD TS is a more diffuse structure, in which all residues of the hydrophobic core including beta-strands 1, 2, 3 and 5 and the major alpha-helix play similar roles in TS stabilization. In addition, Phi-values and DeltaDeltaG(U-++) reveal striking similarities in the TS of Raf RBD and ubiquitin, a structural analogue displaying insignificant sequence identity (<12%). However, ubiquitin TS appears more denatured-like and polarized around the N-terminal beta-hairpin. We suggest that analysis of Phi-values should also consider the direct impact of mutations on differences in free energy between the unfolded and TS (DeltaDeltaG(U-++)) to ensure that the description of TS properties is accurate. Finally, the impact of these findings on the modeling of protein folding is discussed.