Synthetic auxotrophy remains stable after continuous evolution and in coculture with mammalian cells.

Understanding the evolutionary stability and possible context dependence of biological containment techniques is critical as engineered microbes are increasingly under consideration for applications beyond biomanufacturing. While synthetic auxotrophy previously prevented Escherichia coli from exhibiting detectable escape from batch cultures, its long-term effectiveness is unknown. Here, we report automated continuous evolution of a synthetic auxotroph while supplying a decreasing concentration of essential biphenylalanine (BipA). After 100 days of evolution, triplicate populations exhibit no observable escape and exhibit normal growth rates at 10-fold lower BipA concentration than the ancestral synthetic auxotroph. Allelic reconstruction reveals the contribution of three genes to increased fitness at low BipA concentrations. Based on its evolutionary stability, we introduce the progenitor strain directly to mammalian cell culture and observe containment of bacteria without detrimental effects on HEK293T cells. Overall, our findings reveal that synthetic auxotrophy is effective on time scales and in contexts that enable diverse applications.

EGFR-Binding Peptides: From Computational Design towards Tumor-Targeting of Adeno-Associated Virus Capsids.

The epidermal growth factor receptor (EGFR) plays a central role in the progression of many solid tumors. We used this validated target to analyze the de novo design of EGFR-binding peptides and their application for the delivery of complex payloads via rational design of a viral vector. Peptides were computationally designed to interact with the EGFR dimerization interface. Two new peptides and a reference (EDA peptide) were chemically synthesized, and their binding ability characterized. Presentation of these peptides in each of the 60 capsid proteins of recombinant adeno-associated viruses (rAAV) via a genetic based loop insertion enabled targeting of EGFR overexpressing tumor cell lines. Furthermore, tissue distribution and tumor xenograft specificity were analyzed with systemic injection in chicken egg chorioallantoic membrane (CAM) assays. Complex correlations between the targeting of the synthetic peptides and the viral vectors to cells and in ovo were observed. Overall, these data demonstrate the potential of computational design in combination with rational capsid modification for viral vector targeting opening new avenues for viral vector delivery and specifically suicide gene therapy.

Computational design of a modular protein sense-response system.

Sensing and responding to signals is a fundamental ability of living systems, but despite substantial progress in the computational design of new protein structures, there is no general approach for engineering arbitrary new protein sensors. Here, we describe a generalizable computational strategy for designing sensor-actuator proteins by building binding sites de novo into heterodimeric protein-protein interfaces and coupling ligand sensing to modular actuation through split reporters. Using this approach, we designed protein sensors that respond to farnesyl pyrophosphate, a metabolic intermediate in the production of valuable compounds. The sensors are functional in vitro and in cells, and the crystal structure of the engineered binding site closely matches the design model. Our computational design strategy opens broad avenues to link biological outputs to new signals.

A general strategy to construct small molecule biosensors in eukaryotes.

Biosensors for small molecules can be used in applications that range from metabolic engineering to orthogonal control of transcription. Here, we produce biosensors based on a ligand-binding domain (LBD) by using a method that, in principle, can be applied to any target molecule. The LBD is fused to either a fluorescent protein or a transcriptional activator and is destabilized by mutation such that the fusion accumulates only in cells containing the target ligand. We illustrate the power of this method by developing biosensors for digoxin and progesterone. Addition of ligand to yeast, mammalian, or plant cells expressing a biosensor activates transcription with a dynamic range of up to ~100-fold. We use the biosensors to improve the biotransformation of pregnenolone to progesterone in yeast and to regulate CRISPR activity in mammalian cells. This work provides a general methodology to develop biosensors for a broad range of molecules in eukaryotes.

Biocontainment of genetically modified organisms by synthetic protein design.

Genetically modified organisms (GMOs) are increasingly deployed at large scales and in open environments. Genetic biocontainment strategies are needed to prevent unintended proliferation of GMOs in natural ecosystems. Existing biocontainment methods are insufficient because they impose evolutionary pressure on the organism to eject the safeguard by spontaneous mutagenesis or horizontal gene transfer, or because they can be circumvented by environmentally available compounds. Here we computationally redesign essential enzymes in the first organism possessing an altered genetic code (Escherichia coli strain C321.ΔA) to confer metabolic dependence on non-standard amino acids for survival. The resulting GMOs cannot metabolically bypass their biocontainment mechanisms using known environmental compounds, and they exhibit unprecedented resistance to evolutionary escape through mutagenesis and horizontal gene transfer. This work provides a foundation for safer GMOs that are isolated from natural ecosystems by a reliance on synthetic metabolites.

A mechanism for tunable autoinhibition in the structure of a human Ca2+/calmodulin- dependent kinase II holoenzyme.

Calcium/calmodulin-dependent kinase II (CaMKII) forms a highly conserved dodecameric assembly that is sensitive to the frequency of calcium pulse trains. Neither the structure of the dodecameric assembly nor how it regulates CaMKII are known. We present the crystal structure of an autoinhibited full-length human CaMKII holoenzyme, revealing an unexpected compact arrangement of kinase domains docked against a central hub, with the calmodulin-binding sites completely inaccessible. We show that this compact docking is important for the autoinhibition of the kinase domains and for setting the calcium response of the holoenzyme. Comparison of CaMKII isoforms, which differ in the length of the linker between the kinase domain and the hub, demonstrates that these interactions can be strengthened or weakened by changes in linker length. This equilibrium between autoinhibited states provides a simple mechanism for tuning the calcium response without changes in either the hub or the kinase domains.

Assessment of flexible backbone protein design methods for sequence library prediction in the therapeutic antibody Herceptin-HER2 interface.

Computational protein design methods can complement experimental screening and selection techniques by predicting libraries of low-energy sequences compatible with a desired structure and function. Incorporating backbone flexibility in computational design allows conformational adjustments that should broaden the range of predicted low-energy sequences. Here, we evaluate computational predictions of sequence libraries from different protocols for modeling backbone flexibility using the complex between the therapeutic antibody Herceptin and its target human epidermal growth factor receptor 2 (HER2) as a model system. Within the program RosettaDesign, three methods are compared: The first two use ensembles of structures generated by Monte Carlo protocols for near-native conformational sampling: kinematic closure (KIC) and backrub, and the third method uses snapshots from molecular dynamics (MD) simulations. KIC or backrub methods were better able to identify the amino acid residues experimentally observed by phage display in the Herceptin-HER2 interface than MD snapshots, which generated much larger conformational and sequence diversity. KIC and backrub, as well as fixed backbone simulations, captured the key mutation Asp98Trp in Herceptin, which leads to a further threefold affinity improvement of the already subnanomolar parental Herceptin-HER2 interface. Modeling subtle backbone conformational changes may assist in the design of sequence libraries for improving the affinity of antibody-antigen interfaces and could be suitable for other protein complexes for which structural information is available.

ROSETTA3: an object-oriented software suite for the simulation and design of macromolecules.

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

Computer-aided design of functional protein interactions.

Predictive methods for the computational design of proteins search for amino acid sequences adopting desired structures that perform specific functions. Typically, design of 'function' is formulated as engineering new and altered binding activities into proteins. Progress in the design of functional protein-protein interactions is directed toward engineering proteins to precisely control biological processes by specifically recognizing desired interaction partners while avoiding competitors. The field is aiming for strategies to harness recent advances in high-resolution computational modeling-particularly those exploiting protein conformational variability-to engineer new functions and incorporate many functional requirements simultaneously.

Backbone flexibility in computational protein design.

The field of computational protein design has produced striking successes, including the engineering of novel enzymes. Many of these achievements employed methodologies that sample amino acid side-chains on a fixed backbone, while methods that explicitly model backbone flexibility have so far largely focused on the design of new structures rather than functions. Recent methodological improvements in conformational sampling techniques, some borrowed from the field of robotics to model mechanically accessible conformations, now provide exciting opportunities to explore amino acid sequences and backbone structures simultaneously. Incorporating functional constraints into flexible backbone design should help to achieve challenging engineering goals that exploit the role of conformational variability in controlling biological processes, while more generally advancing biophysical understanding of the relationship between variations in protein sequence, structure, dynamics, and function.

Strengths of hydrogen bonds involving phosphorylated amino acid side chains.

Post-translational phosphorylation plays a key role in regulating protein function. Here, we provide a quantitative assessment of the relative strengths of hydrogen bonds involving phosphorylated amino acid side chains (pSer, pAsp) with several common donors (Arg, Lys, and backbone amide groups). We utilize multiple levels of theory, consisting of explicit solvent molecular dynamics, implicit solvent molecular mechanics, and quantum mechanics with a self-consistent reaction field treatment of solvent. Because the approximately 6 pKa of phosphate suggests that -1 and -2 charged species may coexist at physiological pH, hydrogen bonds involving both protonated and deprotonated phosphates for all donor-acceptor pairs are considered. Multiple bonding geometries for the charged-charged interactions are also considered. Arg is shown to be capable of substantially stronger salt bridges with phosphorylated side chains than Lys. A pSer hydrogen-bond acceptor tends to form more stable interactions than a pAsp acceptor. The effect of phosphate protonation state on the strengths of the hydrogen bonds is remarkably subtle, with a more pronounced effect on pAsp than on pSer.

Genome-wide analysis of p53 under hypoxic conditions.

Hypoxia is an important nongenotoxic stress that modulates the tumor suppressor activity of p53 during malignant progression. In this study, we investigated how genotoxic and nongenotoxic stresses regulate p53 association with chromatin, p53 transcriptional activity, and p53-dependent apoptosis. We found that genotoxic and nongenotoxic stresses result in the accumulation and binding of the p53 tumor suppressor protein to the same cognate binding sites in chromatin. However, it is the stress that determines whether downstream signaling is mediated by association with transcriptional coactivators. In contrast to p53 induced by DNA-damaging agents, hypoxia-induced p53 has primarily transrepression activity. Using extensive microarray analysis, we identified families of repressed targets of p53 that are involved in cell signaling, DNA repair, cell cycle control, and differentiation. Following our previous study on the contribution of residues 25 and 26 to p53-dependent hypoxia-induced apoptosis, we found that residues 25-26 and 53-54 and the polyproline- and DNA-binding regions are also required for both gene repression and the induction of apoptosis by p53 during hypoxia. This study defines a new role for residues 53 and 54 of p53 in regulating transrepression and demonstrates that 25-26 and 53-54 work in the same pathway to induce apoptosis through gene repression.