Computational redesign of endonuclease DNA binding and cleavage specificity

Enzymatic excision of uracil residues in nucleosomes depends on the local DNA structure and dynamics

The excision of uracil bases from DNA is accomplished by the enzyme uracil DNA glycosylase (UNG). Recognition of uracil bases in free DNA is facilitated by uracil base pair dynamics, but it is not known whether this same mechanistic feature is relevant for detection and excision of uracil residues embedded in nucleosomes. Here we investigate this question using nucleosome core particles (NCPs) generated from Xenopus laevis histones and the high-affinity "Widom 601" positioning sequence. The reactivity of uracil residues in NCPs under steady-state multiple-turnover conditions was generally decreased compared to that of free 601 DNA, mostly because of anticipated steric effects of histones. However, some sites in NCPs had equal or even greater reactivity than free DNA, and the observed reactivities were not readily explained by simple steric considerations or by global DNA unwrapping models for nucleosome invasion. In particular, some reactive uracils were found in occluded positions, while some unreactive uracils were found in exposed positions. One feature of many exposed reactive sites is a wide DNA minor groove, which allows penetration of a key active site loop of the enzyme. In single-turnover kinetic measurements, multiphasic reaction kinetics were observed for several uracil sites, where each kinetic transient was independent of the UNG concentration. These kinetic measurements, and supporting structural analyses, support a mechanism in which some uracils are transiently exposed to UNG by local, rate-limiting nucleosome conformational dynamics, followed by rapid trapping of the exposed state by the enzyme. We present structural models and plausible reaction mechanisms for the reaction of UNG at three distinct uracil sites in the NCP.

Recent advances in targeted genome engineering in mammalian systems

Targeted genome engineering enables researchers to disrupt, insert, or replace a genomic sequence precisely at a predetermined locus. One well-established technology to edit a mammalian genome is known as gene targeting, which is based on the homologous recombination (HR) mechanism. However, the low HR frequency in mammalian cells (except for mice) prevents its wide application. To address this limitation, a custom-designed nuclease is used to introduce a site-specific DNA double-strand break (DSB) on the chromosome and the subsequent repair of the DSB by the HR mechanism or the non-homologous end joining mechanism results in efficient targeted genome modifications. Engineered homing endonucleases (also called meganucleases), zinc finger nucleases, and transcription activator-like effector nucleases represent the three major classes of custom-designed nucleases that have been successfully applied in many different organisms for targeted genome engineering. This article reviews the recent developments of these genome engineering tools and highlights a few representative applications in mammalian systems. Recent advances in gene delivery strategies of these custom-designed nucleases are also briefly discussed.

An energy-based conformer library for side chain optimization: improved prediction and adjustable sampling

Side chain optimization is a fundamental component of protein modeling applications such as docking, structural prediction, and design. In these applications side chain flexibility is often provided by rotamer or conformer libraries, which are collections of representative side chain conformations. Here we demonstrate that the sampling provided by the library can be substantially improved by adding an energetic criterion to its creation. The result of the new procedure is the Energy-Based library, a conformer library selected according to the propensity of its elements to fit energetically into natural protein environments. The new library performs outstandingly well in side chain optimization, producing structures with significantly lower energies and resulting in improved side chain conformation prediction. In addition, because the library was created as an ordered list, its size can be adjusted to any desired level. This feature provides unprecedented versatility in tuning sampling. It allows to precisely balance the number of conformers required by each amino acid type, equalizing their chances to fit into structural environments. It also allows to scale the amount of sampling to the specific requirement of any given side optimization problem. A rotameric version of the library was also produced with the same method to support applications that require a dihedral-only description of side chain conformation. The libraries are available at http://seneslab.org/EBL.

Role of the biomolecular energy gap in protein design, structure, and evolution

The folding of natural biopolymers into unique three-dimensional structures that determine their function is remarkable considering the vast number of alternative states and requires a large gap in the energy of the functional state compared to the many alternatives. This Perspective explores the implications of this energy gap for computing the structures of naturally occurring biopolymers, designing proteins with new structures and functions, and optimally integrating experiment and computation in these endeavors. Possible parallels between the generation of functional molecules in computational design and natural evolution are highlighted.

Engineering domain fusion chimeras from I-OnuI family LAGLIDADG homing endonucleases

Although engineered LAGLIDADG homing endonucleases (LHEs) are finding increasing applications in biotechnology, their generation remains a challenging, industrial-scale process. As new single-chain LAGLIDADG nuclease scaffolds are identified, however, an alternative paradigm is emerging: identification of an LHE scaffold whose native cleavage site is a close match to a desired target sequence, followed by small-scale engineering to modestly refine recognition specificity. The application of this paradigm could be accelerated if methods were available for fusing N- and C-terminal domains from newly identified LHEs into chimeric enzymes with hybrid cleavage sites. Here we have analyzed the structural requirements for fusion of domains extracted from six single-chain I-OnuI family LHEs, spanning 40–70% amino acid identity. Our analyses demonstrate that both the LAGLIDADG helical interface residues and the linker peptide composition have important effects on the stability and activity of chimeric enzymes. Using a simple domain fusion method in which linker peptide residues predicted to contact their respective domains are retained, and in which limited variation is introduced into the LAGLIDADG helix and nearby interface residues, catalytically active enzymes were recoverable for ∼70% of domain chimeras. This method will be useful for creating large numbers of chimeric LHEs for genome engineering applications.

Rewritable digital data storage in live cells via engineered control of recombination directionality

The use of synthetic biological systems in research, healthcare, and manufacturing often requires autonomous history-dependent behavior and therefore some form of engineered biological memory. For example, the study or reprogramming of aging, cancer, or development would benefit from genetically encoded counters capable of recording up to several hundred cell division or differentiation events. Although genetic material itself provides a natural data storage medium, tools that allow researchers to reliably and reversibly write information to DNA in vivo are lacking. Here, we demonstrate a rewriteable recombinase addressable data (RAD) module that reliably stores digital information within a chromosome. RAD modules use serine integrase and excisionase functions adapted from bacteriophage to invert and restore specific DNA sequences. Our core RAD memory element is capable of passive information storage in the absence of heterologous gene expression for over 100 cell divisions and can be switched repeatedly without performance degradation, as is required to support combinatorial data storage. We also demonstrate how programmed stochasticity in RAD system performance arising from bidirectional recombination can be achieved and tuned by varying the synthesis and degradation rates of recombinase proteins. The serine recombinase functions used here do not require cell-specific cofactors and should be useful in extending computing and control methods to the study and engineering of many biological systems.

LAHEDES: the LAGLIDADG homing endonuclease database and engineering server

LAGLIDADG homing endonucleases (LHEs) are DNA cleaving enzymes, also termed ‘meganucleases’ that are employed as gene-targeting reagents. This use of LHEs requires that their DNA specificity be altered to match sequences in genomic targets. The choice of the most appropriate LHE to target a particular gene is facilitated by the growing number of such enzymes with well-characterized activities and structures. ‘LAHEDES’ (The LAGLIDADG Homing Endonuclease Database and Engineering Server) provides both an online archive of LHEs with validated DNA cleavage specificities and DNA-binding interactions, as well as a tool for the identification of DNA sequences that might be targeted by various LHEs. Searches can be performed using four separate scoring algorithms and user-defined choices of LHE scaffolds. The webserver subsequently provides information regarding clusters of amino acids that should be interrogated during engineering and selection experiments. The webserver is fully open access and can be found at http://homingendonuclease.net.

Targeted gene therapies: tools, applications, optimization

Many devastating human diseases are caused by mutations in a single gene that prevent a somatic cell from carrying out its essential functions, or by genetic changes acquired as a result of infectious disease or in the course of cell transformation. Targeted gene therapies have emerged as potential strategies for treatment of such diseases. These therapies depend upon rare-cutting endonucleases to cleave at specific sites in or near disease genes. Targeted gene correction provides a template for homology-directed repair, enabling the cell's own repair pathways to erase the mutation and replace it with the correct sequence. Targeted gene disruption ablates the disease gene, disabling its function. Gene targeting can also promote other kinds of genome engineering, including mutation, insertion, or gene deletion. Targeted gene therapies present significant advantages compared to approaches to gene therapy that depend upon delivery of stably expressing transgenes. Recent progress has been fueled by advances in nuclease discovery and design, and by new strategies that maximize efficiency of targeting and minimize off-target damage. Future progress will build on deeper mechanistic understanding of critical factors and pathways.

Improved modeling of side-chain--base interactions and plasticity in protein--DNA interface design

Combinatorial sequence optimization for protein design requires libraries of discrete side-chain conformations. The discreteness of these libraries is problematic, particularly for long, polar side chains, since favorable interactions can be missed. Previously, an approach to loop remodeling where protein backbone movement is directed by side-chain rotamers predicted to form interactions previously observed in native complexes (termed "motifs") was described. Here, we show how such motif libraries can be incorporated into combinatorial sequence optimization protocols and improve native complex recapitulation. Guided by the motif rotamer searches, we made improvements to the underlying energy function, increasing recapitulation of native interactions. To further test the methods, we carried out a comprehensive experimental scan of amino acid preferences in the I-AniI protein-DNA interface and found that many positions tolerated multiple amino acids. This sequence plasticity is not observed in the computational results because of the fixed-backbone approximation of the model. We improved modeling of this diversity by introducing DNA flexibility and reducing the convergence of the simulated annealing algorithm that drives the design process. In addition to serving as a benchmark, this extensive experimental data set provides insight into the types of interactions essential to maintain the function of this potential gene therapy reagent.

Restriction endonucleases: natural and directed evolution

Type II restriction endonucleases (REs) are highly sequence-specific compared with other classes of nucleases. PD-(D/E)XK nucleases, initially represented by only type II REs, now comprise a large and extremely diverse superfamily of proteins and, although sharing a structurally conserved core, typically display little or no detectable sequence similarity except for the active site motifs. Sequence similarity can only be observed in methylases and few isoschizomers. As a consequence, REs are classified according to combinations of functional properties rather than on the basis of genetic relatedness. New alignment matrices and classification systems based on structural core connectivity and cleavage mechanisms have been developed to characterize new REs and related proteins. REs recognizing more than 300 distinct specificities have been identified in RE database (REBASE: http://rebase.neb.com/cgi-bin/statlist ) but still the need for newer specificities is increasing due to the advancement in molecular biology and applications. The enzymes have undergone constant evolution through structural changes in protein scaffolds which include random mutations, homologous recombinations, insertions, and deletions of coding DNA sequences but rational mutagenesis or directed evolution delivers protein variants with new functions in accordance with defined biochemical or environmental pressures. Redesigning through random mutation, addition or deletion of amino acids, methylation-based selection, synthetic molecules, combining recognition and cleavage domains from different enzymes, or combination with domains of additional functions change the cleavage specificity or substrate preference and stability. There is a growing number of patents awarded for the creation of engineered REs with new and enhanced properties.

Endonucleases: tools to correct the dystrophin gene

BACKGROUND

Various endonucleases can be engineered to induce double-strand breaks (DSBs) in chosen DNA sequences. These DSBs are spontaneously repaired by nonhomologous-end-joining, resulting in micro-insertions or micro-deletions (INDELs). We detected, characterized and quantified the frequency of INDELs produced by one meganuclease (MGN) targeting the RAG1 gene, six MGNs targeting three introns of the human dystrophin gene and one pair of zinc finger nucleases (ZFNs) targeting exon 50 of the human dystrophin gene. The experiments were performed in human cells (i.e. 293 T cells, myoblasts and myotubes).

METHODS

To analyse the INDELs produced by the endonucleases the targeted region was polymerase chain reaction amplified and the amplicons were digested with the Surveyor enzyme, cloned in bacteria or deep sequenced.

RESULTS

Endonucleases targeting the dystrophin gene produced INDELs of different sizes but there were clear peaks in the size distributions. The positions of these peaks were similar for MGNs but not for ZFNs in 293 T cells and in myoblasts. The size of the INDELs produced by these endonucleases in the dystrophin gene would have permitted a change in the reading frame. In a subsequent experiment, we observed that the frequency of INDELs was increased by re-exposition of the cells to the same endonuclease.

CONCLUSIONS

Endonucleases are able to: (i) restore the normal reading of a gene with a frame shift mutation; (ii) delete a nonsense codon; and (iii) knockout a gene. Endonucleases could thus be used to treat Duchenne muscular dystrophy and other hereditary diseases that are the result of a nonsense codon or a frame shift mutation.

Chemical and biological approaches to improve the efficiency of homologous recombination in human cells mediated by artificial restriction DNA cutter

A chemistry-based artificial restriction DNA cutter (ARCUT) was recently prepared from Ce(IV)/EDTA complex and a pair of pseudo-complementary peptide nucleic acids. This cutter has freely tunable scission-site and site specificity. In this article, homologous recombination (HR) in human cells was promoted by cutting a substrate DNA with ARCUT, and the efficiency of this bioprocess was optimized by various chemical and biological approaches. Of two kinds of terminal structure formed by ARCUT, 3′-overhang termini provided by 1.7-fold higher efficiency than 5′-overhang termini. A longer homology length (e.g. 698 bp) was about 2-fold more favorable than shorter one (e.g. 100 bp). When the cell cycle was synchronized to G2/M phase with nocodazole, the HR was promoted by about 2-fold. Repression of the NHEJ-relevant proteins Ku70 and Ku80 by siRNA increased the efficiency by 2- to 3-fold. It was indicated that appropriate combination of all these chemical and biological approaches should be very effective to promote ARCUT-mediated HR in human cells.

Genome engineering with zinc-finger nucleases

Zinc-finger nucleases (ZFNs) are targetable DNA cleavage reagents that have been adopted as gene-targeting tools. ZFN-induced double-strand breaks are subject to cellular DNA repair processes that lead to both targeted mutagenesis and targeted gene replacement at remarkably high frequencies. This article briefly reviews the history of ZFN development and summarizes applications that have been made to genome editing in many different organisms and situations. Considerable progress has been made in methods for deriving zinc-finger sets for new genomic targets, but approaches to design and selection are still being perfected. An issue that needs more attention is the extent to which available mechanisms of double-strand break repair limit the scope and utility of ZFN-initiated events. The bright prospects for future applications of ZFNs, including human gene therapy, are discussed.

An enumerative stepwise ansatz enables atomic-accuracy RNA loop modeling

Atomic-accuracy structure prediction of macromolecules should be achievable by optimizing a physically realistic energy function but is presently precluded by incomplete sampling of a biopolymer's many degrees of freedom. We present herein a working hypothesis, called the "stepwise ansatz," for recursively constructing well-packed atomic-detail models in small steps, enumerating several million conformations for each monomer, and covering all build-up paths. By making use of high-performance computing and the Rosetta framework, we provide first tests of this hypothesis on a benchmark of 15 RNA loop-modeling problems drawn from riboswitches, ribozymes, and the ribosome, including 10 cases that are not solvable by current knowledge-based modeling approaches. For each loop problem, this deterministic stepwise assembly method either reaches atomic accuracy or exposes flaws in Rosetta's all-atom energy function, indicating the resolution of the conformational sampling bottleneck. As a further rigorous test, we have carried out a blind all-atom prediction for a noncanonical RNA motif, the C7.2 tetraloop/receptor, and validated this model through nucleotide-resolution chemical mapping experiments. Stepwise assembly is an enumerative, ab initio build-up method that systematically outperforms existing Monte Carlo and knowledge-based methods for 3D structure prediction.

Towards artificial metallonucleases for gene therapy: recent advances and new perspectives

The process of DNA targeting or repair of mutated genes within the cell, induced by specifically positioned double-strand cleavage of DNA near the mutated sequence, can be applied for gene therapy of monogenic diseases. For this purpose, highly specific artificial metallonucleases are developed. They are expected to be important future tools of modern genetics. The present state of art and strategies of research are summarized, including protein engineering and artificial ‘chemical’ nucleases. From the results, we learn about the basic role of the metal ions and the various ligands, and about the DNA binding and cleavage mechanism. The results collected provide useful guidance for engineering highly controlled enzymes for use in gene therapy.

RosettaRemodel: a generalized framework for flexible backbone protein design

We describe RosettaRemodel, a generalized framework for flexible protein design that provides a versatile and convenient interface to the Rosetta modeling suite. RosettaRemodel employs a unified interface, called a blueprint, which allows detailed control over many aspects of flexible backbone protein design calculations. RosettaRemodel allows the construction and elaboration of customized protocols for a wide range of design problems ranging from loop insertion and deletion, disulfide engineering, domain assembly, loop remodeling, motif grafting, symmetrical units, to de novo structure modeling.

Divalent metal ion differentially regulates the sequential nicking reactions of the GIY-YIG homing endonuclease I-BmoI

Homing endonucleases are site-specific DNA endonucleases that function as mobile genetic elements by introducing double-strand breaks or nicks at defined locations. Of the major families of homing endonucleases, the modular GIY-YIG endonucleases are least understood in terms of mechanism. The GIY-YIG homing endonuclease I-BmoI generates a double-strand break by sequential nicking reactions during which the single active site of the GIY-YIG nuclease domain must undergo a substantial reorganization. Here, we show that divalent metal ion plays a significant role in regulating the two independent nicking reactions by I-BmoI. Rate constant determination for each nicking reaction revealed that limiting divalent metal ion has a greater impact on the second strand than the first strand nicking reaction. We also show that substrate mutations within the I-BmoI cleavage site can modulate the first strand nicking reaction over a 314-fold range. Additionally, in-gel DNA footprinting with mutant substrates and modeling of an I-BmoI-substrate complex suggest that amino acid contacts to a critical GC-2 base pair are required to induce a bottom-strand distortion that likely directs conformational changes for reaction progress. Collectively, our data implies mechanistic roles for divalent metal ion and substrate bases, suggesting that divalent metal ion facilitates the re-positioning of the GIY-YIG nuclease domain between sequential nicking reactions.

Mining endonuclease cleavage determinants in genomic sequence data

Homing endonucleases have great potential as tools for targeted gene therapy and gene correction, but identifying variants of these enzymes capable of cleaving specific DNA targets of interest is necessary before the widespread use of such technologies is possible. We identified homologues of the LAGLIDADG homing endonuclease I-AniI and their putative target insertion sites by BLAST searches followed by examination of the sequences of the flanking genomic regions. Amino acid substitutions in these homologues that were located close to the target site DNA, and thus potentially conferring differences in target specificity, were grafted onto the I-AniI scaffold. Many of these grafts exhibited novel and unexpected specificities. These findings show that the information present in genomic data can be exploited for endonuclease specificity redesign.

RosettaScripts: a scripting language interface to the Rosetta macromolecular modeling suite

Macromolecular modeling and design are increasingly useful in basic research, biotechnology, and teaching. However, the absence of a user-friendly modeling framework that provides access to a wide range of modeling capabilities is hampering the wider adoption of computational methods by non-experts. RosettaScripts is an XML-like language for specifying modeling tasks in the Rosetta framework. RosettaScripts provides access to protocol-level functionalities, such as rigid-body docking and sequence redesign, and allows fast testing and deployment of complex protocols without need for modifying or recompiling the underlying C++ code. We illustrate these capabilities with RosettaScripts protocols for the stabilization of proteins, the generation of computationally constrained libraries for experimental selection of higher-affinity binding proteins, loop remodeling, small-molecule ligand docking, design of ligand-binding proteins, and specificity redesign in DNA-binding proteins.

Computational design of virus-like protein assemblies on carbon nanotube surfaces

There is a general need for the engineering of protein-like molecules that organize into geometrically specific superstructures on molecular surfaces, directing further functionalization to create richly textured, multilayered assemblies. Here we describe a computational approach whereby the surface properties and symmetry of a targeted surface define the sequence and superstructure of surface-organizing peptides. Computational design proceeds in a series of steps that encode both surface recognition and favorable intersubunit packing interactions. This procedure is exemplified in the design of peptides that assemble into a tubular structure surrounding single-walled carbon nanotubes (SWNTs). The geometrically defined, virus-like coating created by these peptides converts the smooth surfaces of SWNTs into highly textured assemblies with long-scale order, capable of directing the assembly of gold nanoparticles into helical arrays along the SWNT axis.

Computational protein design with a generalized Born solvent model: application to Asparaginyl-tRNA synthetase

Computational Protein Design (CPD) is a promising method for high throughput protein and ligand mutagenesis. Recently, we developed a CPD method that used a polar-hydrogen energy function for protein interactions and a Coulomb/Accessible Surface Area (CASA) model for solvent effects. We applied this method to engineer aspartyl-adenylate (AspAMP) specificity into Asparaginyl-tRNA synthetase (AsnRS), whose substrate is asparaginyl-adenylate (AsnAMP). Here, we implement a more accurate function, with an all-atom energy for protein interactions and a residue-pairwise generalized Born model for solvent effects. As a first test, we compute aminoacid affinities for several point mutants of Aspartyl-tRNA synthetase (AspRS) and Tyrosyl-tRNA synthetase and stability changes for three helical peptides and compare with experiment. As a second test, we readdress the problem of AsnRS aminoacid engineering. We compare three design criteria, which optimize the folding free-energy, the absolute AspAMP affinity, and the relative (AspAMP-AsnAMP) affinity. The sequences and conformations are improved with respect to our previous, polar-hydrogen/CASA study: For several designed complexes, the AspAMP carboxylate forms three interactions with a conserved arginine and a designed lysine, as in the active site of the AspRS:AspAMP complex. The conformations and interactions are well maintained in molecular dynamics simulations and the sequences have an inverted specificity, favoring AspAMP over AsnAMP. The method is not fully successful, since experimental measurements with the seven most promising sequences show that they do not catalyze at a detectable level the adenylation of Asp (or Asn) with ATP. This may be due to weak AspAMP binding and/or disruption of transition-state stabilization.

A synthetic homing endonuclease-based gene drive system in the human malaria mosquito

Genetic methods of manipulating or eradicating disease vector populations have long been discussed as an attractive alternative to existing control measures because of their potential advantages in terms of effectiveness and species specificity^1–3. The development of genetically engineered malaria-resistant mosquitoes has shown, as a proof-of-principle, the possibility of targeting the mosquito’s ability to serve as a disease vector^4–7. The translation of these achievements into control measures requires an effective technology to spread a genetic modification from laboratory mosquitoes to field populations⁸. We have previously suggested that homing endonuclease genes (HEGs), a class of simple selfish genetic elements, could be exploited for this purpose⁹. Here we demonstrate that a synthetic genetic element, consisting of mosquito regulatory regions¹⁰ and the homing endonuclease gene I-SceI^11–13, can substantially increase its transmission to the progeny in transgenic mosquitoes of the human malaria vector Anopheles gambiae. We show that the I-SceI element is able to rapidly invade receptive mosquito cage populations, validating mathematical models for the transmission dynamics of HEGs. Molecular analyses confirm that expression of I-SceI in the male germline induces high rates of site-specific chromosomal cleavage and gene conversion, which results in the gain of the I-SceI gene, and underlies the observed genetic drive. These findings demonstrate a new mechanism by which genetic control measures can be implemented. Our results also show in principle how sequence-specific genetic drive elements like HEGs could be used to take the step from the genetic engineering of individuals to the genetic engineering of populations.

Homing endonucleases: from microbial genetic invaders to reagents for targeted DNA modification

Homing endonucleases are microbial DNA-cleaving enzymes that mobilize their own reading frames by generating double strand breaks at specific genomic invasion sites. These proteins display an economy of size, and yet recognize long DNA sequences (typically 20 to 30 base pairs). They exhibit a wide range of fidelity at individual nucleotide positions in a manner that is strongly influenced by host constraints on the coding sequence of the targeted gene. The activity of these proteins leads to site-specific recombination events that can result in the insertion, deletion, mutation, or correction of DNA sequences. Over the past fifteen years, the crystal structures of representatives from several homing endonuclease families have been solved, and methods have been described to create variants of these enzymes that cleave novel DNA targets. Engineered homing endonucleases proteins are now being used to generate targeted genomic modifications for a variety of biotech and medical applications.

Context dependence between subdomains in the DNA binding interface of the I-CreI homing endonuclease

Homing endonucleases (HE) have emerged as precise tools for achieving gene targeting events. Redesigned HEs with tailored specificities can be used to cleave new sequences, thereby considerably expanding the number of targetable genes and loci. With HEs, as well as with other protein scaffolds, context dependence of DNA/protein interaction patterns remains one of the major limitations for rational engineering of new DNA binders. Previous studies have shown strong crosstalk between different residues and regions of the DNA binding interface. To investigate this phenomenon, we systematically combined mutations from three groups of amino acids in the DNA binding regions of the I-CreI HE. Our results confirm that important crosstalk occurs throughout this interface in I-CreI. Detailed analysis of success rates identified a nearest-neighbour effect, with a more pronounced level of dependence between adjacent regions. Taken together, these data suggest that combinatorial engineering does not necessarily require the identification of separable functional or structural regions, and that groups of amino acids provide acceptable building blocks that can be assembled, overcoming the context dependency of the DNA binding interface. Furthermore, the present work describes a sequential method to engineer tailored HEs, wherein three contiguous regions are individually mutated and assembled to create HEs with engineered specificity.

Analysis of the bacterial luciferase mobile loop by replica-exchange molecular dynamics

Bacterial luciferase contains an extended 29-residue mobile loop. Movements of this loop are governed by binding of either flavin mononucleotide (FMNH2) or polyvalent anions. To understand this process, loop dynamics were investigated using replica-exchange molecular dynamics that yielded conformational ensembles in either the presence or absence of FMNH2. The resulting data were analyzed using clustering and network analysis. We observed the closed conformations that are visited only in the simulations with the ligand. Yet the mobile loop is intrinsically flexible, and FMNH2 binding modifies the relative populations of conformations. This model provides unique information regarding the function of a crystallographically disordered segment of the loop near the binding site. Structures at or near the fringe of this network were compatible with flavin binding or release. Finally, we demonstrate that the crystallographically observed conformation of the mobile loop bound to oxidized flavin was influenced by crystal packing. Thus, our study has revealed what we believe are novel conformations of the mobile loop and additional context for experimentally determined structures.

Insect population control by homing endonuclease-based gene drive: an evaluation in Drosophila melanogaster

Insects play a major role as vectors of human disease as well as causing significant agricultural losses. Harnessing the activity of customized homing endonuclease genes (HEGs) has been proposed as a method for spreading deleterious mutations through populations with a view to controlling disease vectors. Here, we demonstrate the feasibility of this method in Drosophila melanogaster, utilizing the well-characterized HEG, I-SceI. In particular, we show that high rates of homing can be achieved within spermatogonia and in the female germline. We show that homed constructs continue to exhibit HEG activity in the subsequent generation and that the ectopic homing events required for initiating the strategy occur at an acceptable rate. We conclude that the requirements for successful deployment of a HEG-based gene drive strategy can be satisfied in a model dipteran and that there is a reasonable prospect of the method working in other dipterans. In characterizing the system we measured repair outcomes at the spermatogonial, spermatocyte, and spermatid stages of spermatogenesis. We show that homologous recombination is restricted to spermatogonia and that it immediately ceases when they become primary spermatocytes, indicating that the choice of DNA repair pathway in the Drosophila testis can switch abruptly during differentiation.

Determinants of precatalytic conformational transitions in the DNA cytosine methyltransferase M.HhaI

The DNA methyltransferase M.HhaI is an excellent model for understanding how recognition of a nucleic acid substrate is translated into site-specific modification. In this study, we utilize direct, real-time monitoring of the catalytic loop position via engineered tryptophan fluorescence reporters to dissect the conformational transitions that occur in both enzyme and DNA substrate prior to methylation of the target cytosine. Using nucleobase analogues in place of the target and orphan bases, the kinetics of the base flipping and catalytic loop closure rates were determined, revealing that base flipping precedes loop closure as the rate-determining step prior to methyl transfer. To determine the mechanism by which individual specific hydrogen bond contacts at the enzyme-DNA interface mediate these conformational transitions, nucleobase analogues lacking hydrogen bonding groups were incorporated into the recognition sequence to disrupt the major groove recognition elements. The consequences of binding, loop closure, and catalysis were determined for four contacts, revealing large differences in the contribution of individual hydrogen bonds to DNA recognition and conformational transitions on the path to catalysis. Our results describe how M.HhaI utilizes direct readout contacts to accelerate extrication of the target base that offer new insights into the evolutionary history of this important class of enzymes.

Successful targeting and disruption of an integrated reporter lentivirus using the engineered homing endonuclease Y2 I-AniI

Current antiviral therapy does not cure HIV-infected individuals because the virus establishes lifelong latent infection within long-lived memory T cells as integrated HIV proviral DNA. Here, we report a new therapeutic approach that aims to cure cells of latent HIV infection by rendering latent virus incapable of replication and pathogenesis via targeted cellular mutagenesis of essential viral genes. This is achieved by using a homing endonuclease to introduce DNA double-stranded breaks (dsb) within the integrated proviral DNA, which is followed by triggering of the cellular DNA damage response and error-prone repair. To evaluate this concept, we developed an in vitro culture model of viral latency, consisting of an integrated lentiviral vector with an easily evaluated reporter system to detect targeted mutagenesis events. Using this system, we demonstrate that homing endonucleases can efficiently and selectively target an integrated reporter lentivirus within the cellular genome, leading to mutation in the proviral DNA and loss of reporter gene expression. This new technology offers the possibility of selectively disabling integrated HIV provirus within latently infected cells.

Computational design of protein-ligand interfaces: potential in therapeutic development

Computational design of protein-ligand interfaces finds optimal amino acid sequences within a small-molecule binding site of a protein for tight binding of a specific small molecule. It requires a search algorithm that can rapidly sample the vast sequence and conformational space, and a scoring function that can identify low energy designs. This review focuses on recent advances in computational design methods and their application to protein-small molecule binding sites. Strategies for increasing affinity, altering specificity, creating broad-spectrum binding, and building novel enzymes from scratch are described. Future prospects for applications in drug development are discussed, including limitations that will need to be overcome to achieve computational design of protein therapeutics with novel modes of action.

Comprehensive computational design of mCreI homing endonuclease cleavage specificity for genome engineering

Homing endonucleases (HEs) cleave long (∼ 20 bp) DNA target sites with high site specificity to catalyze the lateral transfer of parasitic DNA elements. In order to determine whether comprehensive computational design could be used as a general strategy to engineer new HE target site specificities, we used RosettaDesign (RD) to generate 3200 different variants of the mCreI LAGLIDADG HE towards 16 different base pair positions in the 22 bp mCreI target site. Experimental verification of a range of these designs demonstrated that over 2/3 (24 of 35 designs, 69%) had the intended new site specificity, and that 14 of the 15 attempted specificity shifts (93%) were achieved. These results demonstrate the feasibility of using structure-based computational design to engineer HE variants with novel target site specificities to facilitate genome engineering.

Meganucleases and other tools for targeted genome engineering: perspectives and challenges for gene therapy

The importance of safer approaches for gene therapy has been underscored by a series of severe adverse events (SAEs) observed in patients involved in clinical trials for Severe Combined Immune Deficiency Disease (SCID) and Chromic Granulomatous Disease (CGD). While a new generation of viral vectors is in the process of replacing the classical gamma-retrovirus-based approach, a number of strategies have emerged based on non-viral vectorization and/or targeted insertion aimed at achieving safer gene transfer. Currently, these methods display lower efficacies than viral transduction although many of them can yield more than 1% of engineered cells in vitro. Nuclease-based approaches, wherein an endonuclease is used to trigger site-specific genome editing, can significantly increase the percentage of targeted cells. These methods therefore provide a real alternative to classical gene transfer as well as gene editing. However, the first endonuclease to be in clinic today is not used for gene transfer, but to inactivate a gene (CCR5) required for HIV infection. Here, we review these alternative approaches, with a special emphasis on meganucleases, a family of naturally occurring rare-cutting endonucleases, and speculate on their current and future potential.

Meganuclease-mediated Inhibition of HSV1 Infection in Cultured Cells

Herpes simplex virus type 1 (HSV1) is a major health problem. As for most viral diseases, current antiviral treatments are based on the inhibition of viral replication once it has already started. As a consequence, they impair neither the viral cycle at its early stages nor the latent form of the virus, and thus cannot be considered as real preventive treatments. Latent HSV1 virus could be addressed by rare cutting endonucleases, such as meganucleases. With the aim of a proof of concept study, we generated several meganucleases recognizing HSV1 sequences, and assessed their antiviral activity in cultured cells. We demonstrate that expression of these proteins in African green monkey kidney fibroblast (COS-7) and BSR cells inhibits infection by HSV1, at low and moderate multiplicities of infection (MOIs), inducing a significant reduction of the viral load. Furthermore, the remaining viral genomes display a high rate of mutation (up to 16%) at the meganuclease cleavage site, consistent with a mechanism of action based on the cleavage of the viral genome. This specific mechanism of action qualifies meganucleases as an alternative class of antiviral agent, with the potential to address replicative as well as latent DNA viral forms.

De novo design of peptide-calcite biomineralization systems

Many organisms produce complex, hierarchically structured, inorganic materials via protein-influenced crystal growth--a process known as biomineralization. Understanding this process would shed light on hard-tissue formation and guide efforts to develop biomaterials. We created and tested a computational method to design protein-biomineralization systems. The algorithm folds a protein from a fully extended structure and simultaneously optimizes the fold, orientation, and sequence of the protein adsorbed to a crystal surface. We used the algorithm to design peptides (16 residues) to modify calcite (CaCO(3)) crystallization. We chemically synthesized six peptides that were predicted to bind different states of a calcite growth plane. All six peptides dramatically affected calcite crystal growth (as observed by scanning electron microscopy), and the effects were dependent on the targeted state of the {001} growth plane. Additionally, we synthesized and assayed scrambled variants of all six designed peptides to distinguish cases where sequence composition determines the interactions versus cases where sequence order (and presumably structure) plays a role. Scrambled variants of negatively charged peptides also had dramatic effects on calcite crystallization; in contrast, scrambled variants of positively charged peptides had a variable effect on crystallization, ranging from dramatic to mild. Special emphasis is often placed on acidic protein residues in calcified tissue mineralization; the work presented here suggests an important role for basic residues as well. In particular, this work implicates a potential role for basic residues in sequence-order specificity for peptide-mineral interactions.

The I-CreI meganuclease and its engineered derivatives: applications from cell modification to gene therapy

Meganucleases (MNs) are highly specific enzymes that can induce homologous recombination in different types of cells, including mammalian cells. Consequently, these enzymes are used as scaffolds for the development of custom gene-targeting tools for gene therapy or cell-line development. Over the past 15 years, the high resolution X-ray structures of several MNs from the LAGLIDADG family have improved our understanding of their protein-DNA interaction and mechanism of DNA cleavage. By developing and utilizing high-throughput screening methods to test a large number of variant-target combinations, we have been able to re-engineer scores of I-CreI derivatives into custom enzymes that target a specific DNA sequence of interest. Such customized MNs, along with wild-type ones, have allowed for exploring a large range of biotechnological applications, including protein-expression cell-line development, genetically modified plants and animals and therapeutic applications such as targeted gene therapy as well as a novel class of antivirals.

Evolution of I-SceI homing endonucleases with increased DNA recognition site specificity

Elucidating how homing endonucleases undergo changes in recognition site specificity will facilitate efforts to engineer proteins for gene therapy applications. I-SceI is a monomeric homing endonuclease that recognizes and cleaves within an 18-bp target. It tolerates limited degeneracy in its target sequence, including substitution of a C:G(+4) base pair for the wild-type A:T(+4) base pair. Libraries encoding randomized amino acids at I-SceI residue positions that contact or are proximal to A:T(+4) were used in conjunction with a bacterial one-hybrid system to select I-SceI derivatives that bind to recognition sites containing either the A:T(+4) or the C:G(+4) base pairs. As expected, isolates encoding wild-type residues at the randomized positions were selected using either target sequence. All I-SceI proteins isolated using the C:G(+4) recognition site included small side-chain substitutions at G100 and either contained (K86R/G100T, K86R/G100S and K86R/G100C) or lacked (G100A, G100T) a K86R substitution. Interestingly, the binding affinities of the selected variants for the wild-type A:T(+4) target are 4- to 11-fold lower than that of wild-type I-SceI, whereas those for the C:G(+4) target are similar. The increased specificity of the mutant proteins is also evident in binding experiments in vivo. These differences in binding affinities account for the observed ∼36-fold difference in target preference between the K86R/G100T and wild-type proteins in DNA cleavage assays. An X-ray crystal structure of the K86R/G100T mutant protein bound to a DNA duplex containing the C:G(+4) substitution suggests how sequence specificity of a homing enzyme can increase. This biochemical and structural analysis defines one pathway by which site specificity is augmented for a homing endonuclease.

Develop reusable and combinable designs for transcriptional logic gates

One limit on developing complex synthetic gene circuits is the lack of basic components such as transcriptional logic gates that can process combinatorial inputs. Here, we propose a strategy to construct such components based on reusable designs and convergent reengineering of well-studied natural systems. We demonstrated the strategy using variants of the transcription factor (TF) LacI and operator Olac that form specifically interacting pairs. Guided by a mathematical model derived from existing quantitative knowledge, rational designs of transcriptional NAND, NOR and NOT gates have been realized. The NAND gates have been designed based on direct protein-protein interactions in coupling with DNA looping. We demonstrated that the designs are reusable: a multiplex of logic devices can be readily created using the same designs but different combinations of sequence variants. The designed logic gates are combinable to form compound circuits: a demonstration logic circuit containing all three types of designed logic gates has been synthesized, and the circuit truthfully reproduces the pre-designed input-output logic relations.

Computational protein design using flexible backbone remodeling and resurfacing: case studies in structure-based antigen design

Computational protein design has promise for vaccine design and other applications. We previously transplanted the HIV 4E10 epitope onto non-HIV protein scaffolds for structural stabilization and immune presentation. Here, we developed two methods to optimize the structure of an antigen, flexible backbone remodeling and resurfacing, and we applied these methods to a 4E10 scaffold. In flexible-backbone remodeling, an existing backbone segment is replaced by a de novo designed segment of prespecified length and secondary structure. With remodeling, we replaced a potentially immunodominant domain on the scaffold with a helix-loop segment that made intimate contact to the protein core. All three domain trim designs tested experimentally had improved thermal stability and similar binding affinity for the 4E10 antibody compared to the parent scaffold. A crystal structure of one design had a 0.8 Å backbone RMSD to the computational model in the rebuilt region. Comparison of parent and trimmed scaffold reactivity to anti-parent sera confirmed the deletion of an immunodominant domain. In resurfacing, the surface of an antigen outside a target epitope is redesigned to obtain variants that maintain only the target epitope. Resurfaced variants of two scaffolds were designed in which 50 positions amounting to 40% of the protein sequences were mutated. Surface-patch analyses indicated that most potential antibody footprints outside the 4E10 epitope were altered. The resurfaced variants maintained thermal stability and binding affinity. These results indicate that flexible-backbone remodeling and resurfacing are useful tools for antigen optimization and protein engineering generally.

A protein sequence meta-functional signature for calcium binding residue prediction

The diversity of characterized protein functions found amongst experimentally interrogated proteins suggests that a vast array of unknown functions remains undiscovered. These protein functions are imparted by specific geometric distributions of amino acid residue chemical moieties, each contributing a functional interaction. We hypothesize that individual residue function contributions are predictable through sequence analytic knowledge based algorithms, and that they can be recombined to understand composite protein function by predicting spatial relation in tertiary structure. We assess the former by training a meta-functional signature algorithm to specifically predict calcium ion binding residues from protein sequence. We estimate the latter by testing for match between predictive contribution of positions in predicted secondary structures and patterns of side chain proximity forced by secondary structure moieties. Specific training for calcium binding results in 83% area under the receiver operator characteristic curve added value over random (AUCoR) and p<10(-300) significance as measured by Kendall's τ in ten fold cross validation for parallel sets of 811 residues in 336 proteins and 696 residues in 299 proteins. Training for generalized function results in 63% AUCoR and p≅10(-221) for the same tests. Including inference of side chain proximity improves predictive ability by 2% AUCoR consistently. The results demonstrate that protein meta-functional signatures can be trained to predict specific protein functions by considering amino acid identity and structural features accessible from sequence, laying the groundwork for composite sequence based function site prediction.

Improving computational protein design by using structure-derived sequence profile

Designing a protein sequence that will fold into a predefined structure is of both practical and fundamental interest. Many successful, computational designs in the last decade resulted from improved understanding of hydrophobic and polar interactions between side chains of amino acid residues in stabilizing protein tertiary structures. However, the coupling between main-chain backbone structure and local sequence has yet to be fully addressed. Here, we attempt to account for such coupling by using a sequence profile derived from the sequences of five residue fragments in a fragment library that are structurally matched to the five-residue segments contained in a target structure. We further introduced a term to reduce low complexity regions of designed sequences. These two terms together with optimized reference states for amino-acid residues were implemented in the RosettaDesign program. The new method, called RosettaDesign-SR, makes a 12% increase (from 34 to 46%) in fraction of proteins whose designed sequences are more than 35% identical to wild-type sequences. Meanwhile, it reduces 8% (from 22% to 14%) to the number of designed sequences that are not homologous to any known protein sequences according to psi-blast. More importantly, the sequences designed by RosettaDesign-SR have 2-3% more polar residues at the surface and core regions of proteins and these surface and core polar residues have about 4% higher sequence identity to wild-type sequences than by RosettaDesign. Thus, the proteins designed by RosettaDesign-SR should be less likely to aggregate and more likely to have unique structures due to more specific polar interactions.

Molecular basis of engineered meganuclease targeting of the endogenous human RAG1 locus

Homing endonucleases recognize long target DNA sequences generating an accurate double-strand break that promotes gene targeting through homologous recombination. We have modified the homodimeric I-CreI endonuclease through protein engineering to target a specific DNA sequence within the human RAG1 gene. Mutations in RAG1 produce severe combined immunodeficiency (SCID), a monogenic disease leading to defective immune response in the individuals, leaving them vulnerable to infectious diseases. The structures of two engineered heterodimeric variants and one single-chain variant of I-CreI, in complex with a 24-bp oligonucleotide of the human RAG1 gene sequence, show how the DNA binding is achieved through interactions in the major groove. In addition, the introduction of the G19S mutation in the neighborhood of the catalytic site lowers the reaction energy barrier for DNA cleavage without compromising DNA recognition. Gene-targeting experiments in human cell lines show that the designed single-chain molecule preserves its in vivo activity with higher specificity, further enhanced by the G19S mutation. This is the first time that an engineered meganuclease variant targets the human RAG1 locus by stimulating homologous recombination in human cell lines up to 265 bp away from the cleavage site. Our analysis illustrates the key features for à la carte procedure in protein–DNA recognition design, opening new possibilities for SCID patients whose illness can be treated ex vivo.

Specificity in computational protein design

A long-standing goal of computational protein design is to create proteins similar to those found in Nature. One motivation is to harness the exquisite functional capabilities of proteins for our own purposes. The extent of similarity between designed and natural proteins also reports on how faithfully our models represent the selective pressures that determine protein sequences. As the field of protein design shifts emphasis from reproducing native-like protein structure to function, it has become important that these models treat the notion of specificity in molecular interactions. Although specificity may, in some cases, be achieved by optimization of a desired protein in isolation, methods have been developed to address directly the desire for proteins that exhibit specific functions and interactions.

DNA technologies: what's next applied to microbiology research?

This perspective discusses current DNA technologies used in basic and applied microbiology research and speculates on possible new future technologies. DNA remains one of the most fascinating molecules known to humans and will continue to revolutionize many areas ranging from medicine, food and forensics to robotics and new industrial bioproducts/biofuel from waste materials. What's next with DNA is not always obvious, but history shows the international microbiology research community will readily adopt it.

Sub-AQUA: real-value quality assessment of protein structure models

Computational protein tertiary structure prediction has made significant progress over the past years. However, most of the existing structure prediction methods are not equipped with functionality to predict accuracy of constructed models. Knowing the accuracy of a structure model is crucial for its practical use since the accuracy determines potential applications of the model. Here we have developed quality assessment methods, which predict real value of the global and local quality of protein structure models. The global quality of a model is defined as the root mean square deviation (RMSD) and the LGA score to its native structure. The local quality is defined as the distance between the corresponding Calpha positions of a model and its native structure when they are superimposed. Three regression methods are employed to combine different types of quality assessment measures of models, including alignment-level scores, residue-position level scores, atomic-detailed structure level scores and composite scores. The regression models were tested on a large benchmark data set of template-based protein structure models of various qualities. In predicting RMSD and the LGA score, a combination of two terms, length-normalized SPAD, a score that assesses alignment stability by considering suboptimal alignments, and Verify3D normalized by the square of the model length shows a significant performance, achieving 97.1 and 83.6% accuracy in identifying models with an RMSD of <2 and 6 A, respectively. For predicting the local quality of models, we find that a two-step approach, in which the global RMSD predicted in the first step is further combined with the other terms, can dramatically increase the accuracy. Finally, the developed regression equations are applied to assess the quality of structure models of whole E. coli proteome.