Search

Many proteins have low thermodynamic stability, which can lead to low expression yields and limit functionality in research, industrial and clinical settings. This article introduces two, web-based tools that use the high-resolution structure of a protein along with the Rosetta molecular modeling program to predict stabilizing mutations. The protocols were recently applied to three genetically and structurally distinct proteins and successfully predicted mutations that improved thermal stability and/or protein yield. In all three cases, combining the stabilizing mutations raised the protein unfolding temperatures by more than 20°C. The first protocol evaluates point mutations and can generate a site saturation mutagenesis heatmap. The second identifies mutation clusters around user-defined positions. Both applications only require a protein structure and are particularly valuable when a deep multiple sequence alignment is not available. These tools were created to simplify protein engineering and enable research that would otherwise be infeasible due to poor expression and stability of the native molecule.

Using the molecular modeling program Rosetta, we designed a de novo protein, called SEWN0.1, which binds the heterotrimeric G protein Gα(q.) The design is helical, well‐folded, and primarily monomeric in solution at a concentration of 10 μM. However, when we solved the crystal structure of SEWN0.1 at 1.9 Å, we observed a dimer in a conformation incompatible with binding Gα(q). Unintentionally, we had designed a protein that adopts alternate conformations depending on its oligomeric state. Recently, there has been tremendous progress in the field of protein structure prediction as new methods in artificial intelligence have been used to predict structures with high accuracy. We were curious if the structure prediction method AlphaFold could predict the structure of SEWN0.1 and if the prediction depended on oligomeric state. When AlphaFold was used to predict the structure of monomeric SEWN0.1, it produced a model that resembles the Rosetta design model and is compatible with binding Gα(q), but when used to predict the structure of a dimer, it predicted a conformation that closely resembles the SEWN0.1 crystal structure. AlphaFold's ability to predict multiple conformations for a single protein sequence should be useful for engineering protein switches.

The FastDesign protocol in the molecular modeling program Rosetta iterates between sequence optimization and structure refinement to stabilize de novo designed protein structures and complexes. FastDesign has been used previously to design novel protein folds and assemblies with important applications in research and medicine. To promote sampling of alternative conformations and sequences, FastDesign includes stages where the energy landscape is smoothened by reducing repulsive forces. Here, we discover that this process disfavors larger amino acids in the protein core because the protein compresses in the early stages of refinement. By testing alternative ramping strategies for the repulsive weight, we arrive at a scheme that produces lower energy designs with more native-like sequence composition in the protein core. We further validate the protocol by designing and experimentally characterizing over 4000 proteins and show that the new protocol produces higher stability proteins. This article is protected by copyright. All rights reserved.

Promiscuous enzymes that modify peptides and proteins are powerful tools for labeling biomolecules; however, directing these modifications to desired substrates can be challenging. Here, we use computational interface design to install a substrate recognition domain adjacent to the active site of a promiscuous enzyme, catechol O-methyltransferase. This design approach effectively decouples substrate recognition from the site of catalysis and promotes modification of peptides recognized by the recruitment domain. We determined the crystal structure of this novel multidomain enzyme, SH3-588, which shows that it closely matches our design. SH3-588 methylates directed peptides with catalytic efficiencies exceeding the wild-type enzyme by over 1,000-fold, whereas peptides lacking the directing recognition sequence do not display enhanced efficiencies. In competition experiments, the designer enzyme preferentially modifies directed substrates over undirected substrates, suggesting that we can use designed recruitment domains to direct post-translational modifications to specific sequence motifs on target proteins in complex multisubstrate environments.

Here we introduce Z-lock, an optogenetic approach for reversible, light-controlled steric inhibition of protein active sites. The light oxygen voltage (LOV) domain and Zdk, a small protein that binds LOV selectively in the dark, are appended to the protein of interest where they sterically block the active site. Irradiation causes LOV to change conformation and release Zdk, exposing the active site. Computer-assisted protein design was used to optimize linkers and Zdk-LOV affinity, for both effective binding in the dark, and effective light-induced release of the intramolecular interaction. Z-lock cofilin was shown to have actin severing ability in vitro, and in living cancer cells it produced protrusions and invadopodia. An active fragment of the tubulin acetylase αTAT was similarly modified and shown to acetylate tubulin on irradiation.

Dengue viruses (DENV serotypes 1-4) and Zika virus (ZIKV) are related flaviviruses that continue to be a public health concern, infecting hundreds of millions of people annually. The traditional live-attenuated virus vaccine approach has been challenging for the four DENV serotypes because of the need to achieve balanced replication of four independent vaccine components. Subunit vaccines represent an alternative approach that may circumvent problems inherent with live-attenuated DENV vaccines. In mature virus particles, the envelope (E) protein forms a homodimer that covers the surface of the virus and is the major target of neutralizing antibodies. Many neutralizing antibodies bind to quaternary epitopes that span across both E proteins in the homodimer. For soluble E (sE) protein to be a viable subunit vaccine, the antigens should be easy to produce and retain quaternary epitopes recognized by neutralizing antibodies. However, WT sE proteins are primarily monomeric at conditions relevant for vaccination and exhibit low expression yields. Previously, we identified amino acid mutations that stabilize the sE homodimer from DENV2 and dramatically raise expression yields. Here, we tested whether these same mutations raise the stability of sE from other DENV serotypes and ZIKV. We show that the mutations raise thermostability for sE from all the viruses, increase production yields from 4-fold to 250-fold, stabilize the homodimer, and promote binding to dimer-specific neutralizing antibodies. Our findings suggest that these sE variants could be valuable resources in the efforts to develop effective subunit vaccines for DENV serotypes 1 to 4 and ZIKV.

Engineered, light-sensitive protein switches are used to interrogate a broad variety of biological processes. These switches are typically constructed by genetically fusing naturally occurring light-responsive protein domains with functional domains from other proteins. Protein activity can be controlled using a variety of mechanisms including light-induced colocalization, caging, and allosteric regulation. Protein design efforts have focused on reducing background signaling, maximizing the change in activity upon light stimulation, and perturbing the kinetics of switching. It is common to combine structure-based modeling with experimental screening to identify ideal fusion points between domains and discover point mutations that optimize switching. Here, we introduce commonly used light-sensitive domains and summarize recent progress in using them to regulate protein activity.

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 A 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.

Many scientific disciplines rely on computational methods for data analysis, model generation, and prediction. Implementing these methods is often accomplished by researchers with domain expertise but without formal training in software engineering or computer science. This arrangement has led to underappreciation of sustainability and maintainability of scientific software tools developed in academic environments. Some software tools have avoided this fate, including the scientific library Rosetta. We use this software and its community as a case study to show how modern software development can be accomplished successfully, irrespective of subject area. Rosetta is one of the largest software suites for macromolecular modeling, with 3.1 million lines of code and many state-of-the-art applications. Since the mid 1990s, the software has been developed collaboratively by the RosettaCommons, a community of academics from over 60 institutions worldwide with diverse backgrounds including chemistry, biology, physiology, physics, engineering, mathematics, and computer science. Developing this software suite has provided us with more than two decades of experience in how to effectively develop advanced scientific software in a global community with hundreds of contributors. Here we illustrate the functioning of this development community by addressing technical aspects (like version control, testing, and maintenance), community-building strategies, diversity efforts, software dissemination, and user support. We demonstrate how modern computational research can thrive in a distributed collaborative community. The practices described here are independent of subject area and can be readily adopted by other software development communities.

During signal transduction, the G protein, Gαq, binds and activates phospholipase C-β isozymes. Several diseases have been shown to manifest upon constitutively activating mutation of Gαq, such as uveal melanoma. Therefore, methods are needed to directly inhibit Gαq. Previously, we demonstrated that a peptide derived from a helix-turn-helix (HTH) region of PLC-β3 (residues 852-878) binds Gαq with low micromolar affinity and inhibits Gαq by competing with full-length PLC-β isozymes for binding. Since the HTH peptide is unstructured in the absence of Gαq, we hypothesized that embedding the HTH in a folded protein might stabilize the binding-competent conformation and further improve the potency of inhibition. Using the molecular modeling software Rosetta, we searched the Protein Data Bank for proteins with similar HTH structures near their surface. The candidate proteins were computationally docked against Gαq and their surfaces were redesigned to stabilize this interaction. We then used yeast surface display to affinity mature the designs. The most potent design bound Gαq/i with high affinity in vitro (KD = 18 nM) and inhibited activation of PLC-β isozymes in HEK293 cells. We anticipate that our genetically encoded inhibitor will help interrogate the role of Gαq in healthy and disease model systems. Our work demonstrates that grafting interaction motifs into folded proteins is a powerful approach for generating inhibitors of protein-protein interactions.

Light-sensitive proteins can be used to perturb signaling networks in living cells and animals with high spatiotemporal resolution. We recently engineered a protein heterodimer that dissociates when irradiated with blue light and demonstrated that by fusing each half of the dimer to termini of a protein that it is possible to selectively block binding surfaces on the protein when in the dark. Upon activation with light, the dimer dissociates and exposes the binding surface, allowing the protein to bind its partner. Critical to the success of this system, called Z-lock, is that the linkers connecting the dimer components to the termini are engineered so that the dimer forms over the appropriate binding surface. Here, we develop and test a protocol in the Rosetta molecular modeling program for designing linkers for Z-lock. We show that the protocol can predict the most effective linker sets for three different light-sensitive switches, including a newly designed switch that binds the Rho-family GTPase Cdc42 upon stimulation with blue light. This protocol represents a generalized computational approach to placing a wide variety of proteins under optogenetic control with Z-lock.

Proteins perform an amazingly diverse set of functions in all aspects of life. Critical to the function of many proteins are the highly specific three-dimensional structures they adopt. For this reason, there is strong interest in learning how to rationally design proteins that adopt user-defined structures. Over the last 25 years, there has been significant progress in the field of computational protein design as rotamer-based sequence optimization protocols have enabled accurate design of protein tertiary and quaternary structure. In this award article, I will summarize how the molecular modeling program Rosetta is used to design new protein structures and describe how we have taken advantage of this capability to create proteins that have important applications in research and medicine. I will highlight three protein design stories: the use of protein interface design to create therapeutic bispecific antibodies, the engineering of light-inducible proteins that can be used to recruit proteins to specific locations in the cell, and the de novo design of new protein structures from pieces of naturally occurring proteins.

The prediction of protein three-dimensional structure from amino acid sequence has been a grand challenge problem in computational biophysics for decades, owing to its intrinsic scientific interest and also to the many potential applications for robust protein structure prediction algorithms, from genome interpretation to protein function prediction. More recently, the inverse problem - designing an amino acid sequence that will fold into a specified three-dimensional structure - has attracted growing attention as a potential route to the rational engineering of proteins with functions useful in biotechnology and medicine. Methods for the prediction and design of protein structures have advanced dramatically in the past decade. Increases in computing power and the rapid growth in protein sequence and structure databases have fuelled the development of new data-intensive and computationally demanding approaches for structure prediction. New algorithms for designing protein folds and protein-protein interfaces have been used to engineer novel high-order assemblies and to design from scratch fluorescent proteins with novel or enhanced properties, as well as signalling proteins with therapeutic potential. In this Review, we describe current approaches for protein structure prediction and design and highlight a selection of the successful applications they have enabled.

The Rosetta software for macromolecular modeling, docking and design is extensively used in laboratories worldwide. During two decades of development by a community of laboratories at more than 60 institutions, Rosetta has been continuously refactored and extended. Its advantages are its performance and interoperability between broad modeling capabilities. Here we review tools developed in the last 5 years, including over 80 methods. We discuss improvements to the score function, user interfaces and usability. Rosetta is available at ., This Perspective reviews tools developed over the past five years in the macromolecular modeling, docking and design software Rosetta.

The spread of dengue (DENV) and Zika virus (ZIKV) is a major public health concern. The primary target of antibodies that neutralize DENV and ZIKV is the envelope (E) glycoprotein, and there is interest in using soluble recombinant E (sRecE) proteins as subunit vaccines. However, the most potent neutralizing antibodies against DENV and ZIKV recognize epitopes on the virion surface that span two or more E proteins. Therefore, to create effective DENV and ZIKV vaccines, presentation of these quaternary epitopes may be necessary. The sRecE proteins from DENV and ZIKV crystallize as native-like dimers, but studies in solution suggest that these dimers are marginally stable. To better understand the challenges associated with creating stable sRecE dimers, we characterized the thermostability of sRecE proteins from ZIKV and three DENV serotypes, DENV2-4. All four proteins irreversibly unfolded at moderate temperatures (46-53 °C). At 23 °C and low micromolar concentrations, DENV2 and ZIKV were primarily dimeric, and DENV3-4 were primarily monomeric, whereas at 37 °C, all four proteins were predominantly monomeric. We further show that the dissociation constant for DENV2 dimerization is very temperature-sensitive, ranging from

β‐sheets often have one face packed against the core of the protein and the other facing solvent. Mutational studies have indicated that the solvent‐facing residues can contribute significantly to protein stability, and that the preferred amino acid at each sequence position is dependent on the precise structure of the protein backbone and the identity of the neighboring amino acids. This suggests that the most advantageous methods for designing β‐sheet surfaces will be approaches that take into account the multiple energetic factors at play including side chain rotamer preferences, van der Waals forces, electrostatics, and desolvation effects. Here, we show that the protein design software Rosetta, which models these energetic factors, can be used to dramatically increase protein stability by optimizing interactions on the surfaces of small β‐sheet proteins. Two design variants of the β‐sandwich protein from tenascin were made with 7 and 14 mutations respectively on its β‐sheet surfaces. These changes raised the thermal midpoint for unfolding from 45°C to 64°C and 74°C. Additionally, we tested an empirical approach based on increasing the number of potential salt bridges on the surfaces of the β‐sheets. This was not a robust strategy for increasing stability, as three of the four variants tested were unfolded.

We engineered a photoactivatable system for rapidly and reversibly exporting proteins from the nucleus by embedding a nuclear export signal in the LOV2 domain from phototropin 1. Fusing the chromatin modifier Bre1 to the photoswitch, we achieved light-dependent control of histone H2B monoubiquitylation in yeast, revealing fast turnover of the ubiquitin mark. Moreover, this inducible system allowed us to dynamically monitor the status of epigenetic modifications dependent on H2B ubiquitylation.

Cellular signal transduction is often regulated at multiple steps in order to achieve more complex logic or precise control of a pathway. For instance, some signaling mechanisms couple allosteric activation with localization to achieve high signal to noise. Here, we create a system for light activated nuclear import that incorporates two levels of control. It consists of a nuclear import photoswitch, Light Activated Nuclear Shuttle (LANS), and a protein engineered to preferentially interact with LANS in the dark, Zdk2. First, Zdk2 is tethered to a location in the cytoplasm, which sequesters LANS in the dark. Second, LANS incorporates a nuclear localization signal (NLS) that is sterically blocked from binding to the nuclear import machinery in the dark. When activated with light, LANS both dissociates from its tethered location and exposes its NLS, which leads to nuclear accumulation. We demonstrate that this coupled system improves the dynamic range of LANS in mammalian cells, yeast, and C. elegans and provides tighter control of transcription factors that have been fused to LANS.

The use of bispecific antibodies (BsAbs) to treat human diseases is on the rise. Increasingly complex and powerful therapeutic mechanisms made possible by BsAbs are spurring innovation of novel BsAb formats and methods for their production. The long-lived in vivo pharmacokinetics, optimal biophysical properties and potential effector functions of natural IgG monoclonal (and monospecific) antibodies has resulted in a push to generate fully IgG BsAb formats with the same quaternary structure as monoclonal IgGs. The production of fully IgG BsAbs is challenging because of the highly heterogeneous pairing of heavy chains (HCs) and light chains (LCs) when produced in mammalian cells with two IgG HCs and two LCs. A solution to the HC heterodimerization aspect of IgG BsAb production was first discovered two decades ago; however, addressing the LC mispairing issue has remained intractable until recently. Here, we use computational and rational engineering to develop novel designs to the HC/LC pairing issue, and particularly for κ LCs. Crystal structures of these designs highlight the interactions that provide HC/LC specificity. We produce and characterize multiple fully IgG BsAbs using these novel designs. We demonstrate the importance of specificity engineering in both the variable and constant domains to achieve robust HC/LC specificity within all the BsAbs. These solutions facilitate the production of fully IgG BsAbs for clinical use.

The de novo design of proteins is a rigorous test of our understanding of the key determinants of protein structure. The helix bundle is an interesting de novo design model system due to the diverse topologies that can be generated from a few simple α-helices. Previously, noncomputational studies demonstrated that connecting amphipathic helices together with short loops can sometimes generate helix bundle proteins, regardless of the bundle's exact sequence. However, using such methods, the precise positions of helices and side chains cannot be predetermined. Since protein function depends on exact positioning of residues, we examined if sequence design tools in the program Rosetta could be used to design a four-helix bundle with a predetermined structure. Helix position was specified using a folding procedure that constrained the design model to a defined topology, and iterative rounds of rotamer-based sequence design and backbone refinement were used to identify a low energy sequence for characterization. The designed protein, DND_4HB, unfolds cooperatively (Tm >90°C) and a NMR solution structure shows that it adopts the target helical bundle topology. Helices 2, 3, and 4 agree very closely with the design model (backbone RMSD = 1.11 A) and >90% of the core side chain χ1 and χ2 angles are correctly predicted. Helix 1 lies in the target groove against the other helices, but is displaced 3 A along the bundle axis. This result highlights the potential of computational design to create bundles with atomic-level precision, but also points at remaining challenges for achieving specific positioning between amphipathic helices.

We computationally designed a de novo protein-protein interaction between wild-type ubiquitin and a redesigned scaffold. Our strategy was to incorporate zinc at the designed interface to promote affinity and orientation specificity. A large set of monomeric scaffold surfaces were computationally engineered with three-residue zinc coordination sites, and the ubiquitin residue H68 was docked to the open coordination site to complete a tetrahedral zinc site. This single coordination bond was intended as a hotspot and polar interaction for ubiquitin binding, and surrounding residues on the scaffold were optimized primarily as hydrophobic residues using a rotamer-based sequence design protocol in Rosetta. From thousands of independent design simulations, four sequences were selected for experimental characterization. The best performing design, called Spelter, binds tightly to zinc (Kd

Hypoxia-inducible factors (HIFs) are heterodimeric transcription factors central to hypoxia response and cancer development. Within the HIF-2 complex, one domain (HIF-2α PAS-B) contains a large (290 A3) buried cavity filled with water molecules within its hydrophobic core. Such cavities are uncommon except in the case of ligand-binding proteins, leading to the hypothesis that HIF-2α can be regulated by small molecules. The development of artificial HIF-2α inhibitors validates this hypothesis but raises questions about the impact of this cavity on HIF-2α PAS-B structure and function. To answer these points, we used computational methods to construct a repacked protein containing a smaller cavity within the native fold. Experimental validation of a five-mutation variant confirms achieving these objectives and stabilizing the folded structure. Complementary functional data establish that ligands cannot bind this variant although heterodimerization remains unchanged. Altogether, our strategy innovatively addresses the roles of solvated cavities in maintaining protein stability and function.

Structure-based protein design tests our understanding of the minimal determinants of protein structure and function. Previous studies have demonstrated that placing zinc binding amino acids (His, Glu, Asp or Cys) near each other in a folded protein in an arrangement predicted to be tetrahedral is often sufficient to achieve binding to zinc. However, few designs have been characterized with high-resolution structures. Here, we use X-ray crystallography, binding studies and mutation analysis to evaluate three alternative strategies for designing zinc binding sites with the molecular modeling program Rosetta. While several of the designs were observed to bind zinc, crystal structures of two designs reveal binding configurations that differ from the design model. In both cases, the modeling did not accurately capture the presence or absence of second-shell hydrogen bonds critical in determining binding-site structure. Efforts to more explicitly design second-shell hydrogen bonds were largely unsuccessful as evidenced by mutation analysis and low expression of proteins engineered with extensive primary and secondary networks. Our results suggest that improved methods for designing interaction networks will be needed for creating metal binding sites with high accuracy.

Lysine selectivity is of critical importance during polyubiquitin chain formation because the identity of the lysine controls the biological outcome. Ubiquitins are covalently linked in polyubiquitin chains through one of seven lysine residues on its surface and the C terminus of adjacent protomers. Lys 48-linked polyubiquitin chains signal for protein degradation; however, the structural basis for Lys 48 selectivity remains largely unknown. The ubiquitin-conjugating enzyme Ube2R1/2 has exquisite specificity for Lys 48, and computational docking of Ube2R1/2 and ubiquitin predicts that Lys 48 is guided to the active site through a key electrostatic interaction between Arg 54 on ubiquitin and Asp 143 on Ube2R1/2. The validity of this interaction was confirmed through biochemical experiments. Since structural examples involving Arg 54 in protein-ubiquitin complexes are exceedingly rare, these results provide additional insight into how ubiquitin-protein complexes can be stabilized. We discuss how these findings relate to how other ubiquitin-conjugating enzymes direct the lysine specificity of polyubiquitin chains.

Cell growth and differentiation are controlled by growth factor receptors coupled to the GTPase Ras. Oncogenic mutations disrupt GTPase activity, leading to persistent Ras signaling and cancer progression. Recent evidence indicates that monoubiquitination of Ras leads to Ras activation. Mutation of the primary site of monoubiquitination impairs the ability of activated K-Ras (one of the three mammalian isoforms of Ras) to promote tumor growth. To determine the mechanism of human Ras activation, we chemically ubiquitinated the protein and analyzed its function by NMR, computational modeling and biochemical activity measurements. We established that monoubiquitination has little effect on the binding of Ras to guanine nucleotide, GTP hydrolysis or exchange-factor activation but severely abrogates the response to GTPase-activating proteins in a site-specific manner. These findings reveal a new mechanism by which Ras can trigger persistent signaling in the absence of receptor activation or an oncogenic mutation.

SummaryProtein design tests our understanding of protein stability and structure. Successful design methods should allow the exploration of sequence space not found in nature. However, when redesigning naturally occurring protein structures, most fixed backbone design algorithms return amino acid sequences that share strong sequence identity with wild-type sequences, especially in the protein core. This behavior places a restriction on functional space that can be explored and is not consistent with observations from nature, where sequences of low identity have similar structures. Here, we allow backbone flexibility during design to mutate every position in the core (38 residues) of a four-helix bundle protein. Only small perturbations to the backbone, 1–2 Å, were needed to entirely mutate the core. The redesigned protein, DRNN, is exceptionally stable (melting point >140°C). An NMR and X-ray crystal structure show that the side chains and backbone were accurately modeled (all-atom RMSD = 1.3 Å).

Pathways of ubiquitin-like (UBL) molecule transfer regulate a myriad of cellular cascades. Here, we report a high-throughput assay that correlates catalytic human NEDD8 transfer to bacterial survival. The assay was utilized to screen mutant NEDD8 and NAE (NEDD8-activating enzyme) libraries to engineer a more stable NEDD8 and redesign the NEDD8-NAE interface. This approach will be useful in understanding the specificities underlying UBL pathways.

Summary Mutation of surface residues to charged amino acids increases resistance to aggregation and can enable reversible unfolding. We have developed a protocol using the Rosetta computational design package that "supercharges" proteins while considering the energetic implications of each mutation. Using a homology model, a single-chain variable fragment antibody was designed that has a markedly enhanced resistance to thermal inactivation and displays an unanticipated ≈30-fold improvement in affinity. Such supercharged antibodies should prove useful for assays in resource-limited settings and for developing reagents with improved shelf lives.

Fluorescent biosensors for living cells currently require laborious optimization and a unique design for each target. They are limited by the availability of naturally occurring ligands with appropriate target specificity. Here we describe a biosensor based on an engineered fibronectin monobody scaffold that can be tailored to bind different targets via high throughput screening. This Src family kinase (SFK) biosensor was made by derivatizing a monobody specific for activated SFK with a bright dye whose fluorescence increases upon target binding. We identified sites for dye attachment and alterations to eliminate vesiculation in living cells, providing a generalizable scaffold for biosensor production. This approach minimizes cell perturbation because it senses endogenous, unmodified target, and because sensitivity is enhanced by direct dye excitation. Automated correlation of cell velocities and SFK activity revealed that SFK are activated specifically during protrusion. Activity correlates with velocity, and peaks 1–2 microns from the leading edge.

We substantiate our preliminary description of the class I tryptophanyl-tRNA synthetase minimal catalytic domain with details of its construction, structure, and steady-state kinetic parameters. Generating that active fragment involved deleting 65% of the contemporary enzyme, including the anticodon-binding domain and connecting peptide 1, CP1, a 74-residue internal segment from within the Rossmann fold. We used protein design (Rosetta), rather than phylogenetic sequence alignments, to identify mutations to compensate for the severe loss of modularity, thus restoring stability, as evidenced by renaturation described previously and by 70-ns molecular dynamics simulations. Sufficient solubility to enable biochemical studies was achieved by expressing the redesigned Urzyme as a maltose-binding protein fusion. Michaelis-Menten kinetic parameters from amino acid activation assays showed that, compared with the native full-length enzyme, TrpRS Urzyme binds ATP with similar affinity. This suggests that neither of the two deleted structural modules has a strong influence on ground-state ATP binding. However, tryptophan has 103 lower affinity, and the Urzyme has comparably reduced specificity relative to the related amino acid, tyrosine. Molecular dynamics simulations revealed how CP1 may contribute significantly to cognate amino acid specificity. As class Ia editing domains are nested within the CP1, this finding suggests that this module enhanced amino acid specificity continuously, throughout their evolution. We call this type of reconstructed protein catalyst an Urzyme (Ur prefix indicates original, primitive, or earliest). It establishes a model for recapitulating very early steps in molecular evolution in which fitness may have been enhanced by accumulating entire modules, rather than by discrete amino acid sequence changes.

DNA methylation patterns regulate gene expression programs and are maintained through a highly coordinated process orchestrated by the RING E3 ubiquitin ligase UHRF1. UHRF1 controls DNA methylation inheritance by reading epigenetic modifications to histones and DNA to activate histone H3 ubiquitylation. Here, we find that all five domains of UHRF1, including the previously uncharacterized ubiquitin-like domain (UBL), cooperate for hemi-methylated DNA-dependent H3 ubiquitin ligation. Our structural and biochemical studies, including mutations found in cancer genomes, reveal a bifunctional requirement for the UBL in histone modification: (1) the UBL makes an essential interaction with the backside of the E2 and (2) the UBL coordinates with other UHRF1 domains that recognize epigenetic marks on DNA and histone H3 to direct ubiquitin to H3. Finally, we show UBLs from other E3s also have a conserved interaction with the E2, Ube2D, highlighting a potential prevalence of interactions between UBLs and E2s.

We recently described two protein G variants (NuG1 and NuG2) with redesigned first hairpins that were almost twice as stable, folded 100-fold faster, and had a switched folding mechanism relative to the wild-type protein. To test the structural accuracy of our design algorithm and to provide insights to the dramatic changes in the kinetics and thermodynamics of folding, we have now determined the crystal structures of NuG1 and NuG2 to 1.8 A and 1.85 A, respectively. We find that they adopt hairpin structures that are closer to the computational models than to wild-type protein G; the RMSD of the NuG1 hairpin to the design model and the wild-type structure are 1.7 A and 5.1 A, respectively. The crystallographic B factor in the redesigned first hairpin of NuG1 is systematically higher than the second hairpin, suggesting that the redesigned region is somewhat less rigid. A second round of structure-based design yielded new variants of NuG1 and NuG2, which are further stabilized by 0.5 kcal/mole and 0.9 kcal/mole.

A challenge in the structure-based design of specificity is modeling the negative states, i.e., the complexes that you do not want to form. This is a difficult problem because mutations predicted to destabilize the negative state might be accommodated by small conformational rearrangements. To overcome this challenge, we employ an iterative strategy that cycles between sequence design and protein docking in order to build up an ensemble of alternative negative state conformations for use in specificity prediction. We have applied our technique to the design of heterodimeric CH3 interfaces in the Fc region of antibodies. Combining computationally and rationally designed mutations produced unique designs with heterodimer purities greater than 90%. Asymmetric Fc crystallization was able to resolve the interface mutations; the heterodimer structures confirmed that the interfaces formed as designed. With these CH3 mutations, and those made at the heavy-/light-chain interface, we demonstrate one-step synthesis of four fully IgG-bispecific antibodies.

SummaryB. stearothermophilus tryptophanyl-tRNA synthetase catalysis proceeds via high-energy protein conformations. Unliganded MD trajectories of the pretransition-state complex with Mg2+•ATP and the (post) transition-state analog complex with adenosine tetraphosphate relax rapidly in opposite directions, the former regressing, the latter progressing along the structural reaction coordinate. The two crystal structures (rmsd 0.7 Å) therefore lie on opposite sides of a conformational free-energy maximum as the chemical transition state forms. SNAPP analysis illustrates the complexity of the associated long-range conformational coupling. Switching interactions in four nonpolar core regions are locally isoenergetic throughout the transition. Different configurations, however, propagate their effects to unfavorable, longer-range interactions at the molecular surface. Designed mutation shows that switching interactions enhance the rate, perhaps by destabilizing the ground state immediately before the transition state and limiting nonproductive diffusion before and after the chemical transition state, thereby reducing the activation entropy. This paradigm may apply broadly to energy-transducing enzymes.

The ability to manipulate protein binding affinities is important for the development of proteins as biosensors, industrial reagents, and therapeutics. We have developed a structure-based method to rationally predict single mutations at protein-protein interfaces that enhance binding affinities. The protocol is based on the premise that increasing buried hydrophobic surface area and/or reducing buried hydrophilic surface area will generally lead to enhanced affinity if large steric clashes are not introduced and buried polar groups are not left without a hydrogen bond partner. The procedure selects affinity enhancing point mutations at the protein-protein interface using three criteria: 1) the mutation must be from a polar amino acid to a non-polar amino acid or from a non-polar amino acid to a larger non-polar amino acid, 2) the free energy of binding as calculated with the Rosetta protein modeling program should be more favorable than the free energy of binding calculated for the wild type complex and 3) the mutation should not be predicted to significantly destabilize the monomers. The Rosetta energy function emphasizes short-range interactions: steric repulsion, Van der Waals forces, hydrogen bonding, and an implicit solvation model that penalizes placing atoms adjacent to polar groups. The performance of the computational protocol was experimentally tested on two separate protein complexes; Gαi1 from the heterotrimeric G-protein system bound to the RGS14 GoLoco motif, and the E2, UbcH7, bound to the E3, E6AP from the ubiquitin pathway. 12 single-site mutations that were predicted to be stabilizing were synthesized and characterized in the laboratory. 9 of the 12 mutations successfully increased binding affinity with 5 of these increasing binding by over 1.0 kcal/mol. To further assess our approach we searched the literature for point mutations that pass our criteria and have experimentally determined binding affinities. Of the 8 mutations identified, 5 were accurately predicted to increase binding affinity, further validating the method as a useful tool to increase protein-protein binding affinities.

Recent efforts to design de novo or redesign the sequence and structure of proteins using computational techniques have met with significant success. Most, if not all, of these computational methodologies attempt to model atomic-level interactions, and hence high-resolution structural characterization of the designed proteins is critical for evaluating the atomic-level accuracy of the underlying design force-fields. We previously used our computational protein design protocol RosettaDesign to completely redesign the sequence of the activation domain of human procarboxypeptidase A2. With 68% of the wild-type sequence changed, the designed protein, AYEdesign, is over 10 kcal/mol more stable than the wild-type protein. Here, we describe the high-resolution crystal structure and solution NMR structure of AYEdesign, which show that the experimentally determined backbone and side-chains conformations are effectively superimposable with the computational model at atomic resolution. To isolate the origins of the remarkable stabilization, we have designed and characterized a new series of procarboxypeptidase mutants that gain significant thermodynamic stability with a minimal number of mutations; one mutant gains more than 5 kcal/mol of stability over the wild-type protein with only four amino acid changes. We explore the relationship between force-field smoothing and conformational sampling by comparing the experimentally determined free energies of the overall design and these focused subsets of mutations to those predicted using modified force-fields, and both fixed and flexible backbone sampling protocols.

Cereblon (CRBN), a substrate receptor for the cullin-RING ubiquitin ligase 4 (CRL4) complex, is a direct protein target for thalidomide teratogenicity and antitumor activity of immunomodulatory drugs (IMiDs). Here we report that glutamine synthetase (GS) is an endogenous substrate of CRL4^(CRBN). Upon exposing cells to high glutamine concentration, GS is acetylated at lysines 11 and 14, yielding a degron that is necessary and sufficient for binding and ubiquitylation by CRL4^(CRBN) and degradation by the proteasome. Binding of acetylated degron peptides to CRBN depends on an intact thalidomide-binding pocket but is not competitive with IMiDs. These findings reveal a feedback loop involving CRL4^(CRBN) that adjusts GS protein levels in response to glutamine and uncover a new function for lysine acetylation.

LOVTRAP is an optogenetic approach for reversible light-induced protein dissociation using protein A fragments that bind to the LOV domain only in the dark, with tunable kinetics and a >150-fold change in the dissociation constant (Kd). By reversibly sequestering proteins at mitochondria, we precisely modulated the proteins' access to the cell edge, demonstrating a naturally occurring 3-mHz cell-edge oscillation driven by interactions of Vav2, Rac1, and PI3K proteins.

The structurally defined ubiquitin-like homology fold (UBL) can engage in several unique protein–protein interactions and many of these complexes have been characterized with high-resolution techniques. Using Rosetta's structural classification tools, we have created the Ubiquitin Structural Relational Database (UbSRD), an SQL database of features for all 509 UBL-containing structures in the PDB, allowing users to browse these structures by protein–protein interaction and providing a platform for quantitative analysis of structural features. We used UbSRD to define the recognition features of ubiquitin (UBQ) and SUMO observed in the PDB and the orientation of the UBQ tail while interacting with certain types of proteins. While some of the interaction surfaces on UBQ and SUMO overlap, each molecule has distinct features that aid in molecular discrimination. Additionally, we find that the UBQ tail is malleable and can adopt a variety of conformations upon binding. UbSRD is accessible as an online resource at rosettadesign.med.unc.edu/ubsrd.

Aminoacyl-tRNA synthetases (aaRS) catalyze both chemical steps that translate the universal genetic code. Rodin and Ohno offered an explanation for the existence of two aaRS classes, observing that codons for the most highly conserved Class I active-site residues are anticodons for corresponding Class II active-site residues. They proposed that the two classes arose simultaneously, by translation of opposite strands from the same gene. We have characterized wild-type 46-residue peptides containing ATP-binding sites of Class I and II synthetases and those coded by a gene designed by Rosetta to encode the corresponding peptides on opposite strands. Catalysis by WT and designed peptides is saturable, and the designed peptides are sensitive to active-site residue mutation. All have comparable apparent second-order rate constants 2.9-7.0E-3 M(-1) s(-1) or ∼750,000-1,300,000 times the uncatalyzed rate. The activities of the two complementary peptides demonstrate that the unique information in a gene can have two functional interpretations, one from each complementary strand. The peptides contain phylogenetic signatures of longer, more sophisticated catalysts we call Urzymes and are short enough to bridge the gap between them and simpler uncoded peptides. Thus, they directly substantiate the sense/antisense coding ancestry of Class I and II aaRS. Furthermore, designed 46-mers achieve similar catalytic proficiency to wild-type 46-mers by significant increases in both kcat and Km values, supporting suggestions that the earliest peptide catalysts activated ATP for biosynthetic purposes.

We recently used computational protein design to create an extremely stable, globular protein, Top7, with a sequence and fold not observed previously in nature. Since Top7 was created in the absence of genetic selection, it provides a rare opportunity to investigate aspects of the cellular protein production and surveillance machinery that are subject to natural selection. Here we show that a portion of the Top7 protein corresponding to the final 49 C-terminal residues is efficiently mis-translated and accumulates at high levels in Escherichia coli. We used circular dichroism, size-exclusion chromatography, small-angle X-ray scattering, analytical ultra-centrifugation, and NMR spectroscopy to show that the resulting C-terminal fragment (CFr) protein adopts a compact, extremely stable, homo-dimeric structure. Based on the solution structure, we engineered an even more stable variant of CFr by disulfide-induced covalent circularisation that should be an excellent platform for design of novel functions. The accumulation of high levels of CFr exposes the high error rate of the protein translation machinery. The rarity of correspondingly stable fragments in natural proteins coupled with the observation that high quality ribosome binding sites are found to occur within E. coli protein-coding regions significantly less often than expected by random chance implies a stringent evolutionary pressure against protein sub-fragments that can independently fold into stable structures. The symmetric self-association between two identical mis-translated CFr sub-domains to generate an extremely stable structure parallels a mechanism for natural protein-fold evolution by modular recombination of protein sub-structures.

Protein conformational switches are ubiquitous in nature and often regulate key biological processes. To design new proteins that can switch conformation, protein designers have focused on the two key components of protein switches: the amino acid sequence must be compatible with the multiple target states and there must be a mechanism for perturbing the relative stability of these states. Proteins have been designed that can switch between folded and disordered states, between distinct folded states and between different aggregation states. A variety of trigger mechanisms have been used, including pH shifts, post-translational modification and ligand binding. Recently, computational protein design methods have been applied to switch design. These include algorithms for designing novel ligand-binding sites and simultaneously optimizing a sequence for multiple target structures.

Loss of side-chain conformational entropy is an important force opposing protein folding and the relative preferences of the amino acids for being buried or solvent exposed may be partially determined by which amino acids lose more side-chain entropy when placed in the core of a protein. To investigate these preferences, we have incorporated explicit modeling of side-chain entropy into the protein design algorithm, RosettaDesign. In the standard version of the program, the energy of a particular sequence for a fixed backbone depends only on the lowest energy side-chain conformations that can be identified for that sequence. In the new model, the free energy of a single amino acid sequence is calculated by evaluating the average energy and entropy of an ensemble of structures generated by Monte Carlo sampling of amino acid side-chain conformations. To evaluate the impact of including explicit side-chain entropy, sequences were designed for 110 native protein backbones with and without the entropy model. In general, the differences between the two sets of sequences are modest, with the largest changes being observed for the longer amino acids: methionine and arginine. Overall, the identity between the designed sequences and the native sequences does not increase with the addition of entropy, unlike what is observed when other key terms are added to the model (hydrogen bonding, Lennard-Jones energies, and solvation energies). These results suggest that side-chain conformational entropy has a relatively small role in determining the preferred amino acid at each residue position in a protein.

A previously developed computer program for protein design, RosettaDesign, was used to predict low free energy sequences for nine naturally occurring protein backbones. RosettaDesign had no knowledge of the naturally occurring sequences and on average 65% of the residues in the designed sequences differ from wild-type. Synthetic genes for ten completely redesigned proteins were generated, and the proteins were expressed, purified, and then characterized using circular dichroism, chemical and temperature denaturation and NMR experiments. Although high-resolution structures have not yet been determined, eight of these proteins appear to be folded and their circular dichroism spectra are similar to those of their wild-type counterparts. Six of the proteins have stabilities equal to or up to 7 kcal/mol greater than their wild-type counterparts, and four of the proteins have NMR spectra consistent with a well-packed, rigid structure. These encouraging results indicate that the computational protein design methods can, with significant reliability, identify amino acid sequences compatible with a target protein backbone.

We have improved the original Rosetta centroid/backbone decoy set by increasing the number of proteins and frequency of near native models and by building on sidechains and minimizing clashes. The new set consists of 1,400 model structures for 78 different and diverse protein targets and provides a challenging set for the testing and evaluation of scoring functions. We evaluated the extent to which a variety of all-atom energy functions could identify the native and close-to-native structures in the new decoy sets. Of various implicit solvent models, we found that a solvent-accessible surface area-based solvation provided the best enrichment and discrimination of close-to-native decoys. The combination of this solvation treatment with Lennard Jones terms and the original Rosetta energy provided better enrichment and discrimination than any of the individual terms. The results also highlight the differences in accuracy of NMR and X-ray crystal structures: a large energy gap was observed between native and non-native conformations for X-ray structures but not for NMR structures.

Protein –protein docking algorithms provide a means to elucidate structural details for presently unknown complexes. Here, we present and evaluate a new method to predict protein –protein complexes from the coordinates of the unbound monomer components. The method employs a low-resolution, rigid-body, Monte Carlo search followed by simultaneous optimization of backbone displacement and side-chain conformations using Monte Carlo minimization. Up to 10 5 independent simulations are carried out, and the resulting “decoys” are ranked using an energy function dominated by van der Waals interactions, an implicit solvation model, and an orientation-dependent hydrogen bonding potential. Top-ranking decoys are clustered to select the final predictions. Small-perturbation studies reveal the formation of binding funnels in 42 of 54 cases using coordinates derived from the bound complexes and in 32 of 54 cases using independently determined coordinates of one or both monomers. Experimental binding affinities correlate with the calculated score function and explain the predictive success or failure of many targets. Global searches using one or both unbound components predict at least 25% of the native residue –residue contacts in 28 of the 32 cases where binding funnels exist. The results suggest that the method may soon be useful for generating models of biologically important complexes from the structures of the isolated components, but they also highlight the challenges that must be met to achieve consistent and accurate prediction of protein –protein interactions. q 2003 Elsevier Ltd. All rights reserved

Light-activatable proteins allow precise spatial and temporal control of biological processes in living cells and animals. Several approaches have been developed for controlling protein localization with light, including the conditional inhibition of a nuclear localization signal (NLS) with the Light Oxygen Voltage (AsLOV2) domain of phototropin 1 from Avena sativa. In the dark, the switch adopts a closed conformation that sterically blocks the NLS motif. Upon activation with blue light the C-terminus of the protein unfolds, freeing the NLS to direct the protein to the nucleus. A previous study showed that this approach can be used to control the localization and activity of proteins in mammalian tissue culture cells. Here, we extend this result by characterizing the binding properties of a LOV/NLS switch and demonstrating that it can be used to control gene transcription in yeast. Additionally, we show that the switch, referred to as LANS (light-activated nuclear shuttle), functions in the C. elegans embryo and allows for control of nuclear localization in individual cells. By inserting LANS into the C. elegans lin-1 locus using Cas9-triggered homologous recombination, we demonstrated control of cell fate via light-dependent manipulation of a native transcription factor. We conclude that LANS can be a valuable experimental method for spatial and temporal control of nuclear localization in vivo.

In the ubiquitin-proteasome system, protein substrates are degraded via covalent modification by a polyubiquitin chain. The polyubiquitin chain must be assembled rapidly in cells, because a chain of at least four ubiquitins is required to signal for degradation, and chain-editing enzymes in the cell may cleave premature polyubiquitin chains before achieving this critical length. The ubiquitin-conjugating enzyme Cdc34 and ubiquitin ligase SCF are capable of building polyubiquitin chains onto protein substrates both rapidly and processively; this may be explained at least in part by the atypically fast rate of Cdc34 and SCF association. This rapid association has been attributed to electrostatic interactions between the acidic C-terminal tail of Cdc34 and a feature on SCF called the basic canyon. However, the structural aspects of the Cdc34-SCF interaction and how they permit rapid complex formation remain elusive. Here, we use protein cross-linking to demonstrate that the Cdc34-SCF interaction occurs in multiple conformations, where several residues from the Cdc34 acidic tail are capable of contacting a broad region of the SCF basic canyon. Similar patterns of cross-linking are also observed between Cdc34 and the Cul1 paralog Cul2, implicating the same mechanism for the Cdc34-SCF interaction in other members of the cullin-RING ubiquitin ligases. We discuss how these results can explain the rapid association of Cdc34 and SCF.