Search

We used the molecular modeling program Rosetta to identify clusters of amino acid substitutions in antibody fragments (scFvs and scAbs) that improve global protein stability and resistance to thermal deactivation. Using this methodology, we increased the melting temperature (T(m)) and resistance to heat treatment of an antibody fragment that binds to the Clostridium botulinum hemagglutinin protein (anti-HA33). Two designed antibody fragment variants with two amino acid replacement clusters, designed to stabilize local regions, were shown to have both higher T(m) compared to the parental scFv and importantly, to retain full antigen binding activity after 2 hours of incubation at 70 °C. The crystal structure of one thermostabilized scFv variants was solved at 1.6 Å and shown to be in close agreement with the RosettaAntibody model prediction.

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

Genetic circuits implement computational operations within a cell. Debugging them is difficult because their function is defined by multiple states (e.g., combinations of inputs) that vary in time. Here, we develop RNA‐seq methods that enable the simultaneous measurement of: (i) the states of internal gates, (ii) part performance (promoters, insulators, terminators), and (iii) impact on host gene expression. This is applied to a three‐input one‐output circuit consisting of three sensors, five NOR/NOT gates, and 46 genetic parts. Transcription profiles are obtained for all eight combinations of inputs, from which biophysical models can extract part activities and the response functions of sensors and gates. Various unexpected failure modes are identified, including cryptic antisense promoters, terminator failure, and a sensor malfunction due to media‐induced changes in host gene expression. This can guide the selection of new parts to fix these problems, which we demonstrate by using a bidirectional terminator to disrupt observed antisense transcription. This work introduces RNA‐seq as a powerful method for circuit characterization and debugging that overcomes the limitations of fluorescent reporters and scales to large systems composed of many parts.

DNAplotlib (www.dnaplotlib.org) is a computational toolkit for the programmable visualization of highly customizable, standards-compliant genetic designs. Functions are provided to aid with both visualization tasks and to extract and overlay associated experimental data. High-quality output is produced in the form of vector-based PDFs, rasterized images, and animated movies. All aspects of the rendering process can be easily customized or extended by the user to cover new forms of genetic part or regulation. DNAplotlib supports improved communication of genetic design information and offers new avenues for static, interactive and dynamic visualizations that map and explore the links between the structure and function of genetic parts, devices and systems; including metabolic pathways and genetic circuits. DNAplotlib is crossplatform software developed using Python and released under the MIT license.

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

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.

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.

Computational design of proteins with altered ligand specificity is an emerging method for the creation of new biosensing systems. In this work, we investigated the outcome of site-directed mutagenesis on the Escherichia coli ribose binding protein (RBP), which is frequently used as a design scaffold for computational searches. A ribose biosensor was first constructed whereby an environmentally sensitive fluorescent probe was covalently attached to RBP at position S265C. This protein conjugate displayed a 54% decrease in emission intensity upon the addition of saturating ribose concentrations and exhibited an apparent dissociation constant (K(d) ) of 3.4 microM. Site-directed mutants within the RBP binding pocket were created and examined for ribose binding ability and overall structural stability. Because as many as 12 mutations are needed to alter ligand specificity in RBP, we measured the effect of single and multiple alanine mutations on stability and signal transduction potential of the ribose biosensor. Single alanine mutations had significant impact on both stability and signaling. Mutations of N190A and F214A each produced melting temperatures8 degrees C below those observed for the wild-type protein. Residue Q235, located in the hinge region of RBP, appeared to be a hot spot for global protein stability as well. Additional single alanine mutations demonstrated as much as 200-fold increase in apparent K(d) but retained overall protein stability. The data collected from this study may be incorporated into design algorithms to help create more stable biosensors and optimize signal transduction properties for a variety of important analytes.

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 10 nM) and binds ubiquitin with a Kd of 20 µM in the presence of zinc and 68 µM in the absence of zinc. Mutagenesis studies and nuclear magnetic resonance chemical shift perturbation experiments indicate that Spelter interacts with H68 and the target surface on ubiquitin; however, H68 does not form a hotspot as intended. Instead, mutation of H68 to alanine results in tighter binding. Although a 3/1 zinc coordination arrangement at an interface cannot be ruled out as a means to improve affinity, our study led us to conclude that 2/2 coordination arrangements or multiple-zinc designs are more likely to promote high-affinity protein interactions.

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.

Protein structures evolved through a complex interplay of cooperative interactions and it is still very challenging to design new protein folds de novo. Here, we present a strategy to design self-assembling polypeptide nanostructured polyhedra, based on modularization using orthogonal dimerizing segments. We designed end experimentally demonstrated formation of the tetrahedron that self-assembles from a single polypeptide chain comprising 12 concatenated coiled-coil-forming segments separated by flexible peptide hinges. Path of the polypeptide chain is guided by the defined order of segments that traverse each of the 6 edges of the tetrahedron exactly twice, forming coiled-coil dimers with their corresponding partners. Coincidence of the polypeptide termini in the same vertex is demonstrated by reconstitution of the split fluorescent protein by the polypeptide with the correct tetrahedral topology, while polypeptides with a deleted or scrambled segment order fail to self-assemble correctly. This design platform provides the basis for construction of new topological polypeptide folds based on the set of orthogonal interacting polypeptide segments.

Protein-protein interactions are critical for many biological processes, and over the past several decades, this importance has prompted researchers to investigate the physical and chemical bases of protein binding. How much do we now understand about how proteins interact? Perhaps the most rigorous way to answer this question is to endeavor to design, at an atomic level of detail, such an interaction from scratch. On page 816 of this issue, Fleishman et al. ( 1 ) take on this challenge. They used a computational method that enabled them to design two proteins that bind to a preselected surface on an influenza virus. This work demonstrates how far we have come in understanding and being able to predict protein-protein interactions, and the technique could prove useful for developing biosensors, reagents, and more effective drugs.

Reengineering protein surfaces to exhibit high net charge, referred to as "supercharging", can improve reversibility of unfolding by preventing aggregation of partially unfolded states. Incorporation of charged side chains should be optimized while considering structural and energetic consequences, as numerous mutations and accumulation of like-charges can also destabilize the native state. A previously demonstrated approach deterministically mutates flexible polar residues (amino acids DERKNQ) with the fewest average neighboring atoms per side chain atom (AvNAPSA). Our approach uses Rosetta-based energy calculations to choose the surface mutations. Both protocols are available for use through the ROSIE web server. The automated Rosetta and AvNAPSA approaches for supercharging choose dissimilar mutations, raising an interesting division in surface charging strategy. Rosetta-supercharged variants of GFP (RscG) ranging from -11 to -61 and +7 to +58 were experimentally tested, and for comparison, we re-tested the previously developed AvNAPSA-supercharged variants of GFP (AscG) with +36 and -30 net charge. Mid-charge variants demonstrated ∼3-fold improvement in refolding with retention of stability. However, as we pushed to higher net charges, expression and soluble yield decreased, indicating that net charge or mutational load may be limiting factors. Interestingly, the two different approaches resulted in GFP variants with similar refolding properties. Our results show that there are multiple sets of residues that can be mutated to successfully supercharge a protein, and combining alternative supercharge protocols with experimental testing can be an effective approach for charge-based improvement to refolding.

The Rosetta molecular modeling software package provides experimentally tested and rapidly evolving tools for the 3D structure prediction and high-resolution design of proteins, nucleic acids, and a growing number of non-natural polymers. Despite its free availability to academic users and improving documentation, use of Rosetta has largely remained confined to developers and their immediate collaborators due to the code’s difficulty of use, the requirement for large computational resources, and the unavailability of servers for most of the Rosetta applications. Here, we present a unified web framework for Rosetta applications called ROSIE (Rosetta Online Server that Includes Everyone). ROSIE provides (a) a common user interface for Rosetta protocols, (b) a stable application programming interface for developers to add additional protocols, (c) a flexible back-end to allow leveraging of computer cluster resources shared by RosettaCommons member institutions, and (d) centralized administration by the RosettaCommons to ensure continuous maintenance. This paper describes the ROSIE server infrastructure, a step-by-step ‘serverification’ protocol for use by Rosetta developers, and the deployment of the first nine ROSIE applications by six separate developer teams: Docking, RNA de novo, ERRASER, Antibody, Sequence Tolerance, Supercharge, Beta peptide design, NCBB design, and VIP redesign. As illustrated by the number and diversity of these applications, ROSIE offers a general and speedy paradigm for serverification of Rosetta applications that incurs negligible cost to developers and lowers barriers to Rosetta use for the broader biological community. ROSIE is available at http://rosie.rosettacommons.org.

Recent modeling of filopodia—the actin-based cell organelles employed for sensing and motility—reveals that one of the key limiting factors of filopodial length is diffusional transport of G-actin monomers to the polymerizing barbed ends. We have explored the possibility of active transport of G-actin by myosin motors, which would be an expected biological response to overcome the limitation of a diffusion-based process. We found that in a straightforward implementation of active transport the increase in length was unimpressive, ≤30%, due to sequestering of G-actin by freely diffusing motors. However, artificially removing motor sequestration reactions led to approximately threefold increases in filopodial length, with the transport being mainly limited by the motors failing to detach from the filaments near the tip, clogging the cooperative conveyer belt dynamics. Making motors sterically transparent led to a qualitative change of the dynamics to a different regime of steady growth without a stationary length. Having identified sequestration and clogging as ubiquitous constraints to motor-driven transport, we devised and tested a speculative means to sidestep these limitations in filopodia by employing cross-linking and putative scaffolding roles of Ena/VASP proteins. We conclude that a naïve design of molecular-motor-based active transport would almost always be inefficient—an intricately organized kinetic scheme, with finely tuned rate constants, is required to achieve high-flux transport.