Search

We report the measurement and analysis of sulfonium-π, thioether-π, and ammonium-π interactions in a β-hairpin peptide model system, coupled with computational investigation and PDB analysis. These studies indicated that the sulfonium-π interaction is the strongest and that polarizability contributes to the stronger interaction with sulfonium relative to ammonium. Computational studies demonstrate that differences in solvation of the trimethylsulfonium versus the trimethylammonium group also contribute to the stronger sulfonium-π interaction. In comparing sulfonium-π versus sulfur-π interactions in proteins, analysis of SAM- and SAH-bound enzymes in the PDB suggests that aromatic residues are enriched in close proximity to the sulfur of both SAM and SAH, but the populations of aromatic interactions of the two cofactors are not significantly different, with the exception of the Me-π interactions in SAM, which are the most prevalent interaction in SAM but are not possible for SAH. This suggests that the weaker interaction energies due to loss of the cation-π interaction in going from SAM to SAH may contribute to turnover of the cofactor.

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.

Recombinant T cell receptors (TCRs) can be used to redirect naïve T cells to eliminate virally infected or cancerous cells; however, they are plagued by low stability and uneven expression. Here, we use molecular modeling to identify mutations in the TCR constant domains (Cα/Cβ) that increase the unfolding temperature of Cα/Cβ by 20 °C, improve the expression of four separate α/β TCRs by 3- to 10-fold, and improve the assembly and stability of TCRs with poor intrinsic stability. The stabilizing mutations rescue the expression of TCRs destabilized through variable domain mutation. The improved stability and folding of the TCRs reduces glycosylation, perhaps through conformational stabilization that restricts access to N-linked glycosylation enzymes. The Cα/Cβ mutations enables antibody-like expression and assembly of well-behaved bispecific molecules that combine an anti-CD3 antibody with the stabilized TCR. These TCR/CD3 bispecifics can redirect T cells to kill tumor cells with target HLA/peptide on their surfaces in vitro., Recombinant T-cells receptors can redirect naïve T cells but often have low or uneven expression. Here, the authors model mutations in TCR constant domains to increase T cell receptor expression, assembly, and stability and generate bispecific molecules.

Biomolecular structure drives function, and computational capabilities have progressed such that the prediction and computational design of biomolecular structures is increasingly feasible. Because computational biophysics attracts students from many different backgrounds and with different levels of resources, teaching the subject can be challenging. One strategy to teach diverse learners is with interactive multimedia material that promotes self-paced, active learning. We have created a hands-on education strategy with a set of fifteen modules that teach topics in biomolecular structure and design, from fundamentals of conformational sampling and energy evaluation to applications like protein docking, antibody design, and RNA structure prediction. Our modules are based on PyRosetta, a Python library that encapsulates all computational modules and methods in the Rosetta software package. The workshop-style modules are implemented as Jupyter Notebooks that can be executed in the Google Colaboratory, allowing learners access with just a web browser. The digital format of Jupyter Notebooks allows us to embed images, molecular visualization movies, and interactive coding exercises. This multimodal approach may better reach students from different disciplines and experience levels as well as attract more researchers from smaller labs and cognate backgrounds to leverage PyRosetta in their science and engineering research. All materials are freely available at https://github.com/RosettaCommons/PyRosetta.notebooks.

SxIP-iLID is a novel optogenetic tool designed to assess the temporal role of proteins on microtubule dynamics. The authors establish that optogenetic cross-linking of microtubule and actin networks decreases MT growth velocities and increases the cell area void of microtubules., We developed a novel optogenetic tool, SxIP–improved light-inducible dimer (iLID), to facilitate the reversible recruitment of factors to microtubule (MT) plus ends in an end-binding protein–dependent manner using blue light. We show that SxIP-iLID can track MT plus ends and recruit tgRFP-SspB upon blue light activation. We used this system to investigate the effects of cross-linking MT plus ends and F-actin in Drosophila melanogaster S2 cells to gain insight into spectraplakin function and mechanism. We show that SxIP-iLID can be used to temporally recruit an F-actin binding domain to MT plus ends and cross-link the MT and F-actin networks. Cross-linking decreases MT growth velocities and generates a peripheral MT exclusion zone. SxIP-iLID facilitates the general recruitment of specific factors to MT plus ends with temporal control enabling researchers to systematically regulate MT plus end dynamics and probe MT plus end function in many biological processes., Graphical Abstract

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.

Evolution trains a from-scratch catalyst Metal-bound peptides can catalyze simple reactions such as ester hydrolysis and may have been the starting point for the evolution of modern enzymes. Studer et al. selected progressively more-proficient variants of a small protein derived from a computationally designed zinc-binding peptide. The resulting enzyme could perform the trained reaction at rates typical for naturally evolved enzymes and serendipitously developed a strong preference for a single enantiomer of the substrate. A structure of the final catalyst highlights how small, progressive changes can remodel both catalytic residues and protein architecture in unpredictable ways. Science , this issue p. 1285

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

Summary Dendritic spines are the major loci of synaptic plasticity and are considered as possible structural correlates of memory. Nonetheless, systematic manipulation of specific subsets of spines in the cortex has been unattainable, and thus, the link between spines and memory has been correlational. We developed a novel synaptic optoprobe, AS-PaRac1 (activated synapse targeting photoactivatable Rac1), which can label recently potentiated spines specifically, and induce the selective shrinkage of AS-PaRac1-containing spines. In vivo imaging of AS-PaRac1 revealed that a motor learning induced substantial synaptic remodelling in a small subset of neurons. The acquired motor learning was disrupted by the optical shrinkage of the potentiated spines, whereas it was not affected by the identical manipulation of spines evoked by a distinct motor task in the same cortical region. Taken together, our results demonstrate that a newly acquired motor skill depends on the formation of a task-specific dense synaptic ensemble.

Over the past decade, the Rosetta biomolecular modeling suite has informed diverse biological questions and engineering challenges ranging from interpretation of low-resolution structural data to design of nanomaterials, protein therapeutics, and vaccines. Central to Rosetta’s success is the energy function: amodel parameterized from small molecule and X-ray crystal structure data used to approximate the energy associated with each biomolecule conformation. This paper describes the mathematical models and physical concepts that underlie the latest Rosetta energy function, beta_nov15. Applying these concepts,we explain how to use Rosetta energies to identify and analyze the features of biomolecular models.Finally, we discuss the latest advances in the energy function that extend capabilities from soluble proteins to also include membrane proteins, peptides containing non-canonical amino acids, carbohydrates, nucleic acids, and other macromolecules.

The epigenetic inheritance of DNA methylation requires UHRF1, a histone- and DNA-binding RING E3 ubiquitin ligase that recruits DNMT1 to sites of newly replicated DNA through ubiquitylation of histone H3. UHRF1 binds DNA with selectivity towards hemi-methylated CpGs (HeDNA); however, the contribution of HeDNA sensing to UHRF1 function remains elusive. Here, we reveal that the interaction of UHRF1 with HeDNA is required for DNA methylation but is dispensable for chromatin interaction, which is governed by reciprocal positive cooperativity between the UHRF1 histone- and DNA-binding domains. HeDNA recognition activates UHRF1 ubiquitylation towards multiple lysines on the H3 tail adjacent to the UHRF1 histone-binding site. Collectively, our studies are the first demonstrations of a DNA-protein interaction and an epigenetic modification directly regulating E3 ubiquitin ligase activity. They also define an orchestrated epigenetic control mechanism involving modifications both to histones and DNA that facilitate UHRF1 chromatin targeting, H3 ubiquitylation, and DNA methylation inheritance. DOI: http://dx.doi.org/10.7554/eLife.17101.001, eLife digest Cells are able to regulate the activity of their genes in response to different cues. Genetic information is encoded in DNA and one way to regulate gene activity is to modify the DNA by attaching chemical “epigenetic” markers to it. When a cell divides, these epigenetic markers can be inherited by the daughter cells so that they share the same patterns of gene activity as the parent cell. When the DNA of the parent cell is copied prior to cell division, the epigenetic markers are also copied onto the new DNA. Mistakes in this process are linked to a wide range of diseases in humans, such as cancer and neurological disorders. One type of epigenetic marker is known as a methyl tag and it is added to DNA by certain enzymes in a process called DNA methylation. A protein called UHRF1 is required for human cells to inherit patterns of DNA methylation through cell division. This protein binds to newly copied DNA that lacks some methyl tags as well as to another protein associated with DNA called histone H3. UHRF1 modifies histone H3 by attaching a small protein molecule called ubiquitin to it. This helps to recruit a DNA methylation enzyme to place methyl tags on the newly copied DNA. However, it was not clear how the various properties of UHRF1 allow it to control how DNA methylation is inherited. Harrison et al. addressed this question by studying purified proteins and DNA fragments outside of living cells. The results show that UHRF1 binding to DNA and histone H3 work together to bring UHRF1 to the sites on DNA that require methylation. Further experiments revealed that the methylation pattern on newly copied DNA is able to activate the ability of UHRF1 to place ubiquitin on histone H3. The findings of Harrison et al. reveal a new mechanism by which dividing cells control how DNA methylation is inherited by their daughter cells. A future challenge will be to find out how attaching ubiquitin to histone H3 activates DNA methylation. DOI: http://dx.doi.org/10.7554/eLife.17101.002

Haptotaxis is the process by which cells respond to gradients of substrate-bound cues, such as extracellular matrix proteins (ECM); however, the cellular mechanism of this response remains poorly understood and has mainly been studied by comparing cell behavior on uniform ECMs with different concentrations of components. To study haptotaxis in response to gradients, we utilized microfluidic chambers to generate gradients of the ECM protein fibronectin, and imaged the cell migration response. Lamellipodia are fan-shaped protrusions that are common in migrating cells. Here, we define a new function for lamellipodia and the cellular mechanism required for haptotaxis – differential actin and lamellipodial protrusion dynamics lead to biased cell migration. Modest differences in lamellipodial dynamics occurring over time periods of seconds to minutes are summed over hours to produce differential whole cell movement towards higher concentrations of fibronectin. We identify a specific subset of lamellipodia regulators as being crucial for haptotaxis. Numerous studies have linked components of this pathway to cancer metastasis and, consistent with this, we find that expression of the oncogenic Rac1 P29S mutation abrogates haptotaxis. Finally, we show that haptotaxis also operates through this pathway in 3D environments.

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.

The precise spatio-temporal dynamics of protein activity are often critical in determining cell behaviour, yet for most proteins they remain poorly understood; it remains difficult to manipulate protein activity at precise times and places within living cells. Protein activity has been controlled by light, through protein derivatization with photocleavable moieties or using photoreactive small-molecule ligands. However, this requires use of toxic ultraviolet wavelengths, activation is irreversible, and/or cell loading is accomplished via disruption of the cell membrane (for example, through microinjection). Here we have developed a new approach to produce genetically encoded photoactivatable derivatives of Rac1, a key GTPase regulating actin cytoskeletal dynamics in metazoan cells. Rac1 mutants were fused to the photoreactive LOV (light oxygen voltage) domain from phototropin, sterically blocking Rac1 interactions until irradiation unwound a helix linking LOV to Rac1. Photoactivatable Rac1 (PA-Rac1) could be reversibly and repeatedly activated using 458- or 473-nm light to generate precisely localized cell protrusions and ruffling. Localized Rac activation or inactivation was sufficient to produce cell motility and control the direction of cell movement. Myosin was involved in Rac control of directionality but not in Rac-induced protrusion, whereas PAK was required for Rac-induced protrusion. PA-Rac1 was used to elucidate Rac regulation of RhoA in cell motility. Rac and Rho coordinate cytoskeletal behaviours with seconds and submicrometre precision. Their mutual regulation remains controversial, with data indicating that Rac inhibits and/or activates Rho. Rac was shown to inhibit RhoA in mouse embryonic fibroblasts, with inhibition modulated at protrusions and ruffles. A PA-Rac crystal structure and modelling revealed LOV-Rac interactions that will facilitate extension of this photoactivation approach to other proteins.

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.

Inducible dimers are powerful tools for controlling biological processes through colocalizing signaling molecules. To be effective, an inducible system should have a dissociation constant in the "off" state that is greater (i.e., weaker affinity) than the concentrations of the molecules that are being controlled, and in the "on" state a dissociation constant that is less (i.e., stronger affinity) than the relevant protein concentrations. Here, we reengineer the interaction between the light inducible dimer, iLID, and its binding partner SspB, to better control proteins present at high effective concentrations (5-100 μM). iLID contains a light-oxygen-voltage (LOV) domain that undergoes a conformational change upon activation with blue light and exposes a peptide motif, ssrA, that binds to SspB. The new variant of the dimer system contains a single SspB point mutation (A58V), and displays a 42-fold change in binding affinity when activated with blue light (from 3 ± 2 μM to 125 ± 40 μM) and allows for light-activated colocalization of transmembrane proteins in neurons, where a higher affinity switch (0.8-47 μM) was less effective because more colocalization was seen in the dark. Additionally, with a point mutation in the LOV domain (N414L), we lengthened the reversion half-life of iLID. This expanded suite of light induced dimers increases the variety of cellular pathways that can be targeted with light.

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.

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.

Degenerate codon (DC) libraries efficiently address the experimental library-size limitations of directed evolution by focusing diversity toward the positions and toward the amino acids (AAs) that are most likely to generate hits; however, manually constructing DC libraries is challenging, error prone and time consuming. This paper provides a dynamic programming solution to the task of finding the best DCs while keeping the size of the library beneath some given limit, improving on the existing integer-linear programming formulation. It then extends the algorithm to consider multiple DCs at each position, a heretofore unsolved problem, while adhering to a constraint on the number of primers needed to synthesize the library. In the two library-design problems examined here, the use of multiple DCs produces libraries that very nearly cover the set of desired AAs while still staying within the experimental size limits. Surprisingly, the algorithm is able to find near-perfect libraries where the ratio of amino-acid sequences to nucleic-acid sequences approaches 1; it effectively side-steps the degeneracy of the genetic code. Our algorithm is freely available through our web server and solves most design problems in about a second.

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.

A myriad of innovative bispecific antibody (BsAb) platforms have been reported. Most require significant protein engineering to be viable from a development and manufacturing perspective. Single-chain variable fragments (scFvs) and diabodies that consist only of antibody variable domains have been used as building blocks for making BsAbs for decades. The drawback with Fv-only moieties is that they lack the native-like interactions with CH1/CL domains that make antibody Fab regions stable and soluble. Here, we utilize a redesigned Fab interface to explore 2 novel Fab-based BsAbs platforms. The redesigned Fab interface designs limit heavy and light chain mixing when 2 Fabs are co-expressed simultaneously, thus allowing the use of 2 different Fabs within a BsAb construct without the requirement of one or more scFvs. We describe the stability and activity of a HER2×HER2 IgG-Fab BsAb, and compare its biophysical and activity properties with those of an IgG-scFv that utilizes the variable domains of the same parental antibodies. We also generated an EGFR × CD3 tandem Fab protein with a similar format to a tandem scFv (otherwise known as a bispecific T cell engager or BiTE). We show that the Fab-based BsAbs have superior biophysical properties compared to the scFv-based BsAbs. Additionally, the Fab-based BsAbs do not simply recapitulate the activity of their scFv counterparts, but are shown to possess unique biological activity.

The non-homologous end joining (NHEJ) pathway is used in diverse species to repair chromosome breaks, and is defined in part by a requirement for Ku. We previously demonstrated mammalian Ku has intrinsic 5' deoxyribosephosphate (5'dRP) and apurinic/apyrimidinic (AP) lyase activity, and showed this activity is important for excising abasic site damage from ends. Here we employ systematic mutagenesis to clarify the protein requirements for this activity. We identify lysine 31 in the 70 kD subunit (Ku70 K31) as the primary candidate nucleophile required for catalysis, but additional mutation of Ku70 K160 and six other lysines within Ku80 were required to eliminate all activity. Ku from Saccharomyces cerevisiae also possesses 5'dRP/AP lyase activity, and robust activity was also reliant on lysines in Ku70 analogous to K31 and K160. By comparison, these lysines are not conserved in Xenopus laevis Ku, and Ku from this species has negligible activity. A role for residues flanking Ku70 K31 in expanding the range of abasic site contexts that can be used as substrate was also identified. Our results suggest an active site well located to provide the substrate specificity required for its biological role.

We have developed a protocol for identifying proteins that are predisposed to bind linear epitopes on target proteins of interest. The protocol searches through the protein database for proteins (scaffolds) that are bound to peptides with sequences similar to accessible, linear epitopes on the target protein. The sequence match is considered more significant if residues calculated to be important in the scaffold-peptide interaction are present in the target epitope. The crystal structure of the scaffold-peptide complex is then used as a template for creating a model of the scaffold bound to the target epitope. This model can then be used in conjunction with sequence optimization algorithms or directed evolution methods to search for scaffold mutations that further increase affinity for the target protein. To test the applicability of this approach we targeted three disease-causing proteins: a tuberculosis virulence factor (TVF), the apical membrane antigen (AMA) from malaria, and hemagglutinin from influenza. In each case the best scoring scaffold was tested, and binders with Kds equal to 37 μM and 50 nM for TVF and AMA, respectively, were identified. A web server (http://rosettadesign.med.unc.edu/scaffold/) has been created for performing the scaffold search process with user-defined target sequences.

Accurate energy functions are critical to macromolecular modeling and design. We describe new tools for identifying inaccuracies in energy functions and guiding their improvement, and illustrate the application of these tools to the improvement of the Rosetta energy function. The feature analysis tool identifies discrepancies between structures deposited in the PDB and low-energy structures generated by Rosetta; these likely arise from inaccuracies in the energy function. The optE tool optimizes the weights on the different components of the energy function by maximizing the recapitulation of a wide range of experimental observations. We use the tools to examine three proposed modifications to the Rosetta energy function: improving the unfolded state energy model (reference energies), using bicubic spline interpolation to generate knowledge-based torisonal potentials, and incorporating the recently developed Dunbrack 2010 rotamer library (Shapovalov & Dunbrack, 2011). © 2013 Elsevier Inc.

Here we show that a recent computationally designed zinc-mediated protein interface is serendipitously capable of catalyzing carboxyester and phosphoester hydrolysis. Although the original motivation was to design a de novo zinc-mediated protein-protein interaction (called MID1-zinc), we observed in the homodimer crystal structure a small cleft and open zinc coordination site. We investigated if the cleft and zinc site at the designed interface were sufficient for formation of a primitive active site that can perform hydrolysis. MID1-zinc hydrolyzes 4-nitrophenyl acetate with a rate acceleration of 10(5) and a k(cat)/K(M) of 630 M(-1) s(-1) and 4-nitrophenyl phosphate with a rate acceleration of 10(4) and a k(cat)/K(M) of 14 M(-1) s(-1). These rate accelerations by an unoptimized active site highlight the catalytic power of zinc and suggest that the clefts formed by protein-protein interactions are well-suited for creating enzyme active sites. This discovery has implications for protein evolution and engineering: from an evolutionary perspective, three-coordinated zinc at a homodimer interface cleft represents a simple evolutionary path to nascent enzymatic activity; from a protein engineering perspective, future efforts in de novo design of enzyme active sites may benefit from exploring clefts at protein interfaces for active site placement.

SummaryPhotocontrol of functional peptides is a powerful tool for spatial and temporal control of cell signaling events. We show that the genetically encoded light-sensitive LOV2 domain of Avena Sativa phototropin 1 (AsLOV2) can be used to reversibly photomodulate the affinity of peptides for their binding partners. Sequence analysis and molecular modeling were used to embed two peptides into the Jα helix of the AsLOV2 domain while maintaining AsLOV2 structure in the dark but allowing for binding to effector proteins when the Jα helix unfolds in the light. Caged versions of the ipaA and SsrA peptides, LOV-ipaA and LOV-SsrA, bind their targets with 49- and 8-fold enhanced affinity in the light, respectively. These switches can be used as general tools for light-dependent colocalization, which we demonstrate with photo-activable gene transcription in yeast.

International audience; The CAPRI (Critical Assessment of Predicted Interactions) and CASP (Critical Assessment of protein Structure Prediction) experiments have demonstrated the power of community-wide tests of methodology in assessing the current state of the art and spurring progress in the very challenging areas of protein docking and structure prediction. We sought to bring the power of community-wide experiments to bear on a very challenging protein design problem that provides a complementary but equally fundamental test of current understanding of protein-binding thermodynamics. We have generated a number of designed protein-protein interfaces with very favorable computed binding energies but which do not appear to be formed in experiments, suggesting that there may be important physical chemistry missing in the energy calculations. A total of 28 research groups took up the challenge of determining what is missing: we provided structures of 87 designed complexes and 120 naturally occurring complexes and asked participants to identify energetic contributions and/or structural features that distinguish between the two sets. The community found that electrostatics and solvation terms partially distinguish the designs from the natural complexes, largely due to the nonpolar character of the designed interactions. Beyond this polarity difference, the community found that the designed binding surfaces were, on average, structurally less embedded in the designed monomers, suggesting that backbone conformational rigidity at the designed surface is important for realization of the designed function. These results can be used to improve computational design strategies, but there is still much to be learned; for example, one designed complex, which does form in experiments, was classified by all metrics as a nonbinder.

The inhibitory switch (IS) domain of p21-activated kinase 1 (PAK1) stabilizes full-length PAK1 in an inactive conformation by binding to the PAK1 kinase domain. Competitive binding of small guanosine triphosphatases to the IS domain disrupts the autoinhibitory interactions and exposes the IS domain binding site on the surface of the kinase domain. To build an affinity reagent that selectively binds the activated state of PAK1, we used molecular modeling to reengineer the isolated IS domain so that it was soluble and stable, did not bind to guanosine triphosphatases and bound more tightly to the PAK1 kinase domain. Three design strategies were tested: in the first and second cases, extension and redesign of the N-terminus were used to expand the hydrophobic core of the domain, and in the third case, the termini were redesigned to be adjacent in space so that the domain could be stabilized by insertion into a loop in a host cyan fluorescent protein (CFP). The best-performing design, called CFP-PAcKer, was based on the third strategy and bound the kinase domain of PAK1 with an affinity of 400 nM. CFP-PAcKer binds more tightly to a full-length variant of PAK1 that is stabilized in the “open” state (Kd = 3.3 μM) than to full-length PAK1 in the “closed” state (undetectable affinity), and binding can be monitored with fluorescence by placing an environmentally sensitive fluorescence dye on CFP-PAcKer adjacent to the binding site.

Background Few existing protein-protein interface design methods allow for extensive backbone rearrangements during the design process. There is also a dichotomy between redesign methods, which take advantage of the native interface, and de novo methods, which produce novel binders. Methodology Here, we propose a new method for designing novel protein reagents that combines advantages of redesign and de novo methods and allows for extensive backbone motion. This method requires a bound structure of a target and one of its natural binding partners. A key interaction in this interface, the anchor, is computationally grafted out of the partner and into a surface loop on the design scaffold. The design scaffold's surface is then redesigned with backbone flexibility to create a new binding partner for the target. Careful choice of a scaffold will bring experimentally desirable characteristics into the new complex. The use of an anchor both expedites the design process and ensures that binding proceeds against a known location on the target. The use of surface loops on the scaffold allows for flexible-backbone redesign to properly search conformational space. Conclusions and Significance This protocol was implemented within the Rosetta3 software suite. To demonstrate and evaluate this protocol, we have developed a benchmarking set of structures from the PDB with loop-mediated interfaces. This protocol can recover the correct loop-mediated interface in 15 out of 16 tested structures, using only a single residue as an anchor.

The de novo design of protein-binding peptides is challenging because it requires the identification of both a sequence and a backbone conformation favorable for binding. We used a computational strategy that iterates between structure and sequence optimization to redesign the C-terminal portion of the RGS14 GoLoco motif peptide so that it adopts a new conformation when bound to Gα(i1). An X-ray crystal structure of the redesigned complex closely matches the computational model, with a backbone root-mean-square deviation of 1.1 Å.

We have recently completed a full re-architecturing of the Rosetta molecular modeling program, generalizing and expanding its existing functionality. The new architecture enables the rapid prototyping of novel protocols by providing easy to use interfaces to powerful tools for molecular modeling. The source code of this rearchitecturing has been released as Rosetta3 and is freely available for academic use. At the time of its release, it contained 470,000 lines of code. Counting currently unpublished protocols at the time of this writing, the source includes 1,285,000 lines. Its rapid growth is a testament to its ease of use. This document describes the requirements for our new architecture, justifies the design decisions, sketches out central classes, and highlights a few of the common tasks that the new software can perform.

During ubiquitin conjugation, the thioester bond that links ‘donor’ ubiquitin to ubiquitin-conjugating enzyme (E2) undergoes nucleophilic attack by the ε-amino group of an acceptor lysine, resulting in formation of an isopeptide bond. Models of ubiquitination have envisioned the donor ubiquitin to be a passive participant in this process. However, we show here that the I44A mutation in ubiquitin profoundly inhibits its ability to serve as a donor for ubiquitin chain initiation or elongation, but can be rescued by computationally-predicted compensatory mutations in the E2 Cdc34. The donor defect of ubiquitin-I44A can be partially suppressed either by using a low pKa amine (hydroxylamine) as the acceptor or by performing reactions at higher pH, suggesting that the discharge defect arises in part due to inefficient deprotonation of the acceptor lysine. We propose that interaction between Cdc34 and the donor ubiquitin organizes the active site to promote efficient ubiquitination of substrate.

GoLoco motif proteins bind to the inhibitory G(i) subclass of G-protein α subunits and slow the release of bound GDP; this interaction is considered critical to asymmetric cell division and neuro-epithelium and epithelial progenitor differentiation. To provide protein tools for interrogating the precise cellular role(s) of GoLoco motif/Gα(i) complexes, we have employed structure-based protein design strategies to predict gain-of-function mutations that increase GoLoco motif binding affinity. Here, we describe fluorescence polarization and isothermal titration calorimetry measurements showing three predicted Gα(i1) point mutations, E116L, Q147L, and E245L; each increases affinity for multiple GoLoco motifs. A component of this affinity enhancement results from a decreased rate of dissociation between the Gα mutants and GoLoco motifs. For Gα(i1)(Q147L), affinity enhancement was seen to be driven by favorable changes in binding enthalpy, despite reduced contributions from binding entropy. The crystal structure of Gα(i1)(Q147L) bound to the RGS14 GoLoco motif revealed disorder among three peptide residues surrounding a well defined Leu-147 side chain. Monte Carlo simulations of the peptide in this region showed a sampling of multiple backbone conformations in contrast to the wild-type complex. We conclude that mutation of Glu-147 to leucine creates a hydrophobic surface favorably buried upon GoLoco peptide binding, yet the hydrophobic Leu-147 also promotes flexibility among residues 511-513 of the RGS14 GoLoco peptide.

The importance of a protein-protein interaction to a signaling pathway can be established by showing that amino acid mutations that weaken the interaction disrupt signaling, and that additional mutations that rescue the interaction recover signaling. Identifying rescue mutations, often referred to as second-site suppressor mutations, controls against scenarios in which the initial deleterious mutation inactivates the protein or disrupts alternative protein-protein interactions. Here, we test a structure-based protocol for identifying second-site suppressor mutations that is based on a strategy previously described by Kortemme and Baker. The molecular modeling software Rosetta is used to scan an interface for point mutations that are predicted to weaken binding but can be rescued by mutations on the partner protein. The protocol typically identifies three types of specificity switches: knob-in-to-hole redesigns, switching hydrophobic interactions to hydrogen bond interactions, and replacing polar interactions with nonpolar interactions. Computational predictions were tested with two separate protein complexes; the G-protein Galpha(i1) bound to the RGS14 GoLoco motif, and UbcH7 bound to the ubiquitin ligase E6AP. Eight designs were experimentally tested. Swapping a buried hydrophobic residue with a polar residue dramatically weakened binding affinities. In none of these cases were we able to identify compensating mutations that returned binding to wild-type affinity, highlighting the challenges inherent in designing buried hydrogen bond networks. The strongest specificity switches were a knob-in-to-hole design (20-fold) and the replacement of a charge-charge interaction with nonpolar interactions (55-fold). In two cases, specificity was further tuned by including mutations distant from the initial design. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Emerging evidence suggests that ubiquitination serves as a protein trafficking signal in addition to its well characterized role in promoting protein degradation. The yeast G protein α subunit Gpa1 represents a rare example of a protein that undergoes both mono- and poly-ubiquitination. Whereas mono-ubiquitinated Gpa1 is targeted to the vacuole, poly-ubiquitinated Gpa1 is directed instead to the proteasome. Here we investigate the structural requirements for mono- and poly-ubiquitination of Gpa1. We find that variants of Gpa1 engineered to be unstable are more likely to be poly-ubiquitinated and less likely to be mono-ubiquitinated. In addition, mutants that cannot be myristoylated are no longer mono-ubiquitinated but are still polyubiquitinated. Finally, we show that the ubiquitin ligase Rsp5 is necessary for Gpa1 mono-ubiquitination in vivo and that the purified enzyme is sufficient to catalyze Gpa1 mono-ubiquitination in vitro. Taken together, these data indicate that mono- and poly-ubiquitination have distinct enzyme and substrate recognition requirements; whereas poly-ubiquitination targets misfolded protein for degradation, a distinct ubiquitination apparatus targets the fully mature, fully myristoylated G protein for mono-ubiquitination and delivery to the vacuole.

The computer-based design of protein-protein interactions is a rigorous test of our understanding of molecular recognition and an attractive approach for creating novel tools for cell and molecular research. Considerable attention has been placed on redesigning the affinity and specificity of naturally occurring interactions. Several studies have shown that reducing the desolvation costs for binding while preserving shape complimentarity and hydrogen bonding is an effective strategy for improving binding affinities. In favorable cases specificity has been designed by focusing only on interactions with the target protein, while in cases with closely related off-target proteins, it has been necessary to explicitly disfavor unwanted binding partners. The rational design of protein-protein interactions from scratch is still an unsolved problem, but recent developments in flexible backbone design and energy functions hold promise for the future.