Structure, substrate binding and activity of a unique AAA+ protein: the BrxL phage restriction factor.

Bacteriophage exclusion ('BREX') systems are multi-protein complexes encoded by a variety of bacteria and archaea that restrict phage by an unknown mechanism. One BREX factor, termed BrxL, has been noted to display sequence similarity to various AAA+ protein factors including Lon protease. In this study we describe multiple CryoEM structures of BrxL that demonstrate it to be a chambered, ATP-dependent DNA binding protein. The largest BrxL assemblage corresponds to a dimer of heptamers in the absence of bound DNA, versus a dimer of hexamers when DNA is bound in its central pore. The protein displays DNA-dependent ATPase activity, and ATP binding promotes assembly of the complex on DNA. Point mutations within several regions of the protein-DNA complex alter one or more in vitro behaviors and activities, including ATPase activity and ATP-dependent association with DNA. However, only the disruption of the ATPase active site fully eliminates phage restriction, indicating that other mutations can still complement BrxL function within the context of an otherwise intact BREX system. BrxL displays significant structural homology to MCM subunits (the replicative helicase in archaea and eukaryotes), implying that it and other BREX factors may collaborate to disrupt initiation of phage DNA replication.

De novo design of a fluorescence-activating ß-barrel.

The regular arrangements of ß-strands around a central axis in ß-barrels and of α-helices in coiled coils contrast with the irregular tertiary structures of most globular proteins, and have fascinated structural biologists since they were first discovered. Simple parametric models have been used to design a wide range of α-helical coiled-coil structures, but to date there has been no success with ß-barrels. Here we show that accurate de novo design of ß-barrels requires considerable symmetry-breaking to achieve continuous hydrogen-bond connectivity and eliminate backbone strain. We then build ensembles of ß-barrel backbone models with cavity shapes that match the fluorogenic compound DFHBI, and use a hierarchical grid-based search method to simultaneously optimize the rigid-body placement of DFHBI in these cavities and the identities of the surrounding amino acids to achieve high shape and chemical complementarity. The designs have high structural accuracy and bind and fluorescently activate DFHBI in vitro and in Escherichia coli, yeast and mammalian cells. This de novo design of small-molecule binding activity, using backbones custom-built to bind the ligand, should enable the design of increasingly sophisticated ligand-binding proteins, sensors and catalysts that are not limited by the backbone geometries available in known protein structures.

Identification and characterization of the WYL BrxR protein and its gene as separable regulatory elements of a BREX phage restriction system.

Bacteriophage exclusion ('BREX') phage restriction systems are found in a wide range of bacteria. Various BREX systems encode unique combinations of proteins that usually include a site-specific methyltransferase; none appear to contain a nuclease. Here we describe the identification and characterization of a Type I BREX system from Acinetobacter and the effect of deleting each BREX ORF on growth, methylation, and restriction. We identified a previously uncharacterized gene in the BREX operon that is dispensable for methylation but involved in restriction. Biochemical and crystallographic analyses of this factor, which we term BrxR ('BREX Regulator'), demonstrate that it forms a homodimer and specifically binds a DNA target site upstream of its transcription start site. Deletion of the BrxR gene causes cell toxicity, reduces restriction, and significantly increases the expression of BrxC. In contrast, the introduction of a premature stop codon into the BrxR gene, or a point mutation blocking its DNA binding ability, has little effect on restriction, implying that the BrxR coding sequence and BrxR protein play independent functional roles. We speculate that elements within the BrxR coding sequence are involved in cis regulation of anti-phage activity, while the BrxR protein itself plays an additional regulatory role, perhaps during horizontal transfer.

Rational design of α-helical tandem repeat proteins with closed architectures.

Tandem repeat proteins, which are formed by repetition of modular units of protein sequence and structure, play important biological roles as macromolecular binding and scaffolding domains, enzymes, and building blocks for the assembly of fibrous materials. The modular nature of repeat proteins enables the rapid construction and diversification of extended binding surfaces by duplication and recombination of simple building blocks. The overall architecture of tandem repeat protein structures--which is dictated by the internal geometry and local packing of the repeat building blocks--is highly diverse, ranging from extended, super-helical folds that bind peptide, DNA, and RNA partners, to closed and compact conformations with internal cavities suitable for small molecule binding and catalysis. Here we report the development and validation of computational methods for de novo design of tandem repeat protein architectures driven purely by geometric criteria defining the inter-repeat geometry, without reference to the sequences and structures of existing repeat protein families. We have applied these methods to design a series of closed α-solenoid repeat structures (α-toroids) in which the inter-repeat packing geometry is constrained so as to juxtapose the amino (N) and carboxy (C) termini; several of these designed structures have been validated by X-ray crystallography. Unlike previous approaches to tandem repeat protein engineering, our design procedure does not rely on template sequence or structural information taken from natural repeat proteins and hence can produce structures unlike those seen in nature. As an example, we have successfully designed and validated closed α-solenoid repeats with a left-handed helical architecture that--to our knowledge--is not yet present in the protein structure database.

Structure, subunit organization and behavior of the asymmetric Type IIT restriction endonuclease BbvCI.

BbvCI, a Type IIT restriction endonuclease, recognizes and cleaves the seven base pair sequence 5'-CCTCAGC-3', generating 3-base, 5'-overhangs. BbvCI is composed of two protein subunits, each containing one catalytic site. Either site can be inactivated by mutation resulting in enzyme variants that nick DNA in a strand-specific manner. Here we demonstrate that the holoenzyme is labile, with the R1 subunit dissociating at low pH. Crystallization of the R2 subunit under such conditions revealed an elongated dimer with the two catalytic sites located on opposite sides. Subsequent crystallization at physiological pH revealed a tetramer comprising two copies of each subunit, with a pair of deep clefts each containing two catalytic sites appropriately positioned and oriented for DNA cleavage. This domain organization was further validated with single-chain protein constructs in which the two enzyme subunits were tethered via peptide linkers of variable length. We were unable to crystallize a DNA-bound complex; however, structural similarity to previously crystallized restriction endonucleases facilitated creation of an energy-minimized model bound to DNA, and identification of candidate residues responsible for target recognition. Mutation of residues predicted to recognize the central C:G base pair resulted in an altered enzyme that recognizes and cleaves CCTNAGC (N = any base).

Computational design of ligand-binding proteins with high affinity and selectivity.

The ability to design proteins with high affinity and selectivity for any given small molecule is a rigorous test of our understanding of the physiochemical principles that govern molecular recognition. Attempts to rationally design ligand-binding proteins have met with little success, however, and the computational design of protein-small-molecule interfaces remains an unsolved problem. Current approaches for designing ligand-binding proteins for medical and biotechnological uses rely on raising antibodies against a target antigen in immunized animals and/or performing laboratory-directed evolution of proteins with an existing low affinity for the desired ligand, neither of which allows complete control over the interactions involved in binding. Here we describe a general computational method for designing pre-organized and shape complementary small-molecule-binding sites, and use it to generate protein binders to the steroid digoxigenin (DIG). Of seventeen experimentally characterized designs, two bind DIG; the model of the higher affinity binder has the most energetically favourable and pre-organized interface in the design set. A comprehensive binding-fitness landscape of this design, generated by library selections and deep sequencing, was used to optimize its binding affinity to a picomolar level, and X-ray co-crystal structures of two variants show atomic-level agreement with the corresponding computational models. The optimized binder is selective for DIG over the related steroids digitoxigenin, progesterone and ß-oestradiol, and this steroid binding preference can be reprogrammed by manipulation of explicitly designed hydrogen-bonding interactions. The computational design method presented here should enable the development of a new generation of biosensors, therapeutics and diagnostics.

The effect of increasing numbers of repeats on TAL effector DNA binding specificity.

Transcription activator-like effectors (TALEs) recognize their DNA targets via tandem repeats, each specifying a single nucleotide base in a one-to-one sequential arrangement. Due to this modularity and their ability to bind long DNA sequences with high specificity, TALEs have been used in many applications. Contributions of individual repeat-nucleotide associations to affinity and specificity have been characterized. Here, using in vitro binding assays, we examined the relationship between the number of repeats in a TALE and its affinity, for both target and non-target DNA. Each additional repeat provides extra binding energy for the target DNA, with the gain decaying exponentially such that binding energy saturates. Affinity for non-target DNA also increases non-linearly with the number of repeats, but with a slower decay of gain. The difference between the effect of length on affinity for target versus non-target DNA manifests in specificity increasing then diminishing with increasing TALE length, peaking between 15 and 19 repeats. Modeling across different hypothetical saturation levels and rates of gain decay, reflecting different repeat compositions, yielded a similar range of specificity optima. This range encompasses the mean and median length of native TALEs, suggesting that these proteins as a group have evolved for maximum specificity.

CSAR Benchmark Exercise 2013: Evaluation of Results from a Combined Computational Protein Design, Docking, and Scoring/Ranking Challenge.

Community Structure-Activity Resource (CSAR) conducted a benchmark exercise to evaluate the current computational methods for protein design, ligand docking, and scoring/ranking. The exercise consisted of three phases. The first phase required the participants to identify and rank order which designed sequences were able to bind the small molecule digoxigenin. The second phase challenged the community to select a near-native pose of digoxigenin from a set of decoy poses for two of the designed proteins. The third phase investigated the ability of current methods to rank/score the binding affinity of 10 related steroids to one of the designed proteins (pKd = 4.1 to 6.7). We found that 11 of 13 groups were able to correctly select the sequence that bound digoxigenin, with most groups providing the correct three-dimensional structure for the backbone of the protein as well as all atoms of the active-site residues. Eleven of the 14 groups were able to select the appropriate pose from a set of plausible decoy poses. The ability to predict absolute binding affinities is still a difficult task, as 8 of 14 groups were able to correlate scores to affinity (Pearson-r > 0.7) of the designed protein for congeneric steroids and only 5 of 14 groups were able to correlate the ranks of the 10 related ligands (Spearman-ρ > 0.7).

Comprehensive homing endonuclease target site specificity profiling reveals evolutionary constraints and enables genome engineering applications.

Homing endonucleases (HEs) promote the evolutionary persistence of selfish DNA elements by catalyzing element lateral transfer into new host organisms. The high site specificity of this lateral transfer reaction, termed homing, reflects both the length (14-40 bp) and the limited tolerance of target or homing sites for base pair changes. In order to better understand molecular determinants of homing, we systematically determined the binding and cleavage properties of all single base pair variant target sites of the canonical LAGLIDADG homing endonucleases I-CreI and I-MsoI. These Chlorophyta algal HEs have very similar three-dimensional folds and recognize nearly identical 22 bp target sites, but use substantially different sets of DNA-protein contacts to mediate site-specific recognition and cleavage. The site specificity differences between I-CreI and I-MsoI suggest different evolutionary strategies for HE persistence. These differences also provide practical guidance in target site finding, and in the generation of HE variants with high site specificity and cleavage activity, to enable genome engineering applications.

The structural basis of mRNA recognition and binding by yeast pseudouridine synthase PUS1.

The chemical modification of RNA bases represents a ubiquitous activity that spans all domains of life. Pseudouridylation is the most common RNA modification and is observed within tRNA, rRNA, ncRNA and mRNAs. Pseudouridine synthase or 'PUS' enzymes include those that rely on guide RNA molecules and others that function as 'stand-alone' enzymes. Among the latter, several have been shown to modify mRNA transcripts. Although recent studies have defined the structural requirements for RNA to act as a PUS target, the mechanisms by which PUS1 recognizes these target sequences in mRNA are not well understood. Here we describe the crystal structure of yeast PUS1 bound to an RNA target that we identified as being a hot spot for PUS1-interaction within a model mRNA at 2.4 Å resolution. The enzyme recognizes and binds both strands in a helical RNA duplex, and thus guides the RNA containing the target uridine to the active site for subsequent modification of the transcript. The study also allows us to show the divergence of related PUS1 enzymes and their corresponding RNA target specificities, and to speculate on the basis by which PUS1 binds and modifies mRNA or tRNA substrates.

De novo design of knotted tandem repeat proteins.

De novo protein design methods can create proteins with folds not yet seen in nature. These methods largely focus on optimizing the compatibility between the designed sequence and the intended conformation, without explicit consideration of protein folding pathways. Deeply knotted proteins, whose topologies may introduce substantial barriers to folding, thus represent an interesting test case for protein design. Here we report our attempts to design proteins with trefoil (31) and pentafoil (51) knotted topologies. We extended previously described algorithms for tandem repeat protein design in order to construct deeply knotted backbones and matching designed repeat sequences (N = 3 repeats for the trefoil and N = 5 for the pentafoil). We confirmed the intended conformation for the trefoil design by X ray crystallography, and we report here on this protein's structure, stability, and folding behaviour. The pentafoil design misfolded into an asymmetric structure (despite a 5-fold symmetric sequence); two of the four repeat-repeat units matched the designed backbone while the other two diverged to form local contacts, leading to a trefoil rather than pentafoil knotted topology. Our results also provide insights into the folding of knotted proteins.

Computational design of sequence-specific DNA-binding proteins.

Sequence-specific DNA-binding proteins (DBPs) play critical roles in biology and biotechnology, and there has been considerable interest in the engineering of DBPs with new or altered specificities for genome editing and other applications. While there has been some success in reprogramming naturally occurring DBPs using selection methods, the computational design of new DBPs that recognize arbitrary target sites remains an outstanding challenge. We describe a computational method for the design of small DBPs that recognize specific target sequences through interactions with bases in the major groove, and employ this method in conjunction with experimental screening to generate binders for 5 distinct DNA targets. These binders exhibit specificity closely matching the computational models for the target DNA sequences at as many as 6 base positions and affinities as low as 30-100 nM. The crystal structure of a designed DBP-target site complex is in close agreement with the design model, highlighting the accuracy of the design method. The designed DBPs function in both Escherichia coli and mammalian cells to repress and activate transcription of neighboring genes. Our method is a substantial step towards a general route to small and hence readily deliverable sequence-specific DBPs for gene regulation and editing.

High-resolution profiling of homing endonuclease binding and catalytic specificity using yeast surface display.

Experimental analysis and manipulation of protein-DNA interactions pose unique biophysical challenges arising from the structural and chemical homogeneity of DNA polymers. We report the use of yeast surface display for analytical and selection-based applications for the interaction between a LAGLIDADG homing endonuclease and its DNA target. Quantitative flow cytometry using oligonucleotide substrates facilitated a complete profiling of specificity, both for DNA-binding and catalysis, with single base pair resolution. These analyses revealed a comprehensive segregation of binding specificity and affinity to one half of the pseudo-dimeric interaction, while the entire interface contributed specificity at the level of catalysis. A single round of targeted mutagenesis with tandem affinity and catalytic selection steps provided mechanistic insights to the origins of binding and catalytic specificity. These methods represent a dynamic new approach for interrogating specificity in protein-DNA interactions.

Design of functionalised circular tandem repeat proteins with longer repeat topologies and enhanced subunit contact surfaces.

Circular tandem repeat proteins ('cTRPs') are de novo designed protein scaffolds (in this and prior studies, based on antiparallel two-helix bundles) that contain repeated protein sequences and structural motifs and form closed circular structures. They can display significant stability and solubility, a wide range of sizes, and are useful as protein display particles for biotechnology applications. However, cTRPs also demonstrate inefficient self-assembly from smaller subunits. In this study, we describe a new generation of cTRPs, with longer repeats and increased interaction surfaces, which enhanced the self-assembly of two significantly different sizes of homotrimeric constructs. Finally, we demonstrated functionalization of these constructs with (1) a hexameric array of peptide-binding SH2 domains, and (2) a trimeric array of anti-SARS CoV-2 VHH domains. The latter proved capable of sub-nanomolar binding affinities towards the viral receptor binding domain and potent viral neutralization function.

Incorporation of sensing modalities into de novo designed fluorescence-activating proteins.

Through the efforts of many groups, a wide range of fluorescent protein reporters and sensors based on green fluorescent protein and its relatives have been engineered in recent years. Here we explore the incorporation of sensing modalities into de novo designed fluorescence-activating proteins, called mini-fluorescence-activating proteins (mFAPs), that bind and stabilize the fluorescent cis-planar state of the fluorogenic compound DFHBI. We show through further design that the fluorescence intensity and specificity of mFAPs for different chromophores can be tuned, and the fluorescence made sensitive to pH and Ca2+ for real-time fluorescence reporting. Bipartite split mFAPs enable real-time monitoring of protein-protein association and (unlike widely used split GFP reporter systems) are fully reversible, allowing direct readout of association and dissociation events. The relative ease with which sensing modalities can be incorporated and advantages in smaller size and photostability make de novo designed fluorescence-activating proteins attractive candidates for optical sensor engineering.

Engineering and functionalization of large circular tandem repeat protein nanoparticles.

Protein engineering has enabled the design of molecular scaffolds that display a wide variety of sizes, shapes, symmetries and subunit compositions. Symmetric protein-based nanoparticles that display multiple protein domains can exhibit enhanced functional properties due to increased avidity and improved solution behavior and stability. Here we describe the creation and characterization of a computationally designed circular tandem repeat protein (cTRP) composed of 24 identical repeated motifs, which can display a variety of functional protein domains (cargo) at defined positions around its periphery. We demonstrate that cTRP nanoparticles can self-assemble from smaller individual subunits, can be produced from prokaryotic and human expression platforms, can employ a variety of cargo attachment strategies and can be used for applications (such as T-cell culture and expansion) requiring high-avidity molecular interactions on the cell surface.

Unintended specificity of an engineered ligand-binding protein facilitated by unpredicted plasticity of the protein fold.

Attempts to create novel ligand-binding proteins often focus on formation of a binding pocket with shape complementarity against the desired ligand (particularly for compounds that lack distinct polar moieties). Although designed proteins often exhibit binding of the desired ligand, in some cases they display unintended recognition behavior. One such designed protein, that was originally intended to bind tetrahydrocannabinol (THC), was found instead to display binding of 25-hydroxy-cholecalciferol (25-D3) and was subjected to biochemical characterization, further selections for enhanced 25-D3 binding affinity and crystallographic analyses. The deviation in specificity is due in part to unexpected altertion of its conformation, corresponding to a significant change of the orientation of an α-helix and an equally large movement of a loop, both of which flank the designed ligand-binding pocket. Those changes led to engineered protein constructs that exhibit significantly more contacts and complementarity towards the 25-D3 ligand than the initial designed protein had been predicted to form towards its intended THC ligand. Molecular dynamics simulations imply that the initial computationally designed mutations may contribute to the movement of the helix. These analyses collectively indicate that accurate prediction and control of backbone dynamics conformation, through a combination of improved conformational sampling and/or de novo structure design, represents a key area of further development for the design and optimization of engineered ligand-binding proteins.

Sampling and energy evaluation challenges in ligand binding protein design.

The steroid hormone 17α-hydroxylprogesterone (17-OHP) is a biomarker for congenital adrenal hyperplasia and hence there is considerable interest in development of sensors for this compound. We used computational protein design to generate protein models with binding sites for 17-OHP containing an extended, nonpolar, shape-complementary binding pocket for the four-ring core of the compound, and hydrogen bonding residues at the base of the pocket to interact with carbonyl and hydroxyl groups at the more polar end of the ligand. Eight of 16 designed proteins experimentally tested bind 17-OHP with micromolar affinity. A co-crystal structure of one of the designs revealed that 17-OHP is rotated 180° around a pseudo-two-fold axis in the compound and displays multiple binding modes within the pocket, while still interacting with all of the designed residues in the engineered site. Subsequent rounds of mutagenesis and binding selection improved the ligand affinity to nanomolar range, while appearing to constrain the ligand to a single bound conformation that maintains the same "flipped" orientation relative to the original design. We trace the discrepancy in the design calculations to two sources: first, a failure to model subtle backbone changes which alter the distribution of sidechain rotameric states and second, an underestimation of the energetic cost of desolvating the carbonyl and hydroxyl groups of the ligand. The difference between design model and crystal structure thus arises from both sampling limitations and energy function inaccuracies that are exacerbated by the near two-fold symmetry of the molecule.

Precise assembly of complex beta sheet topologies from de novo designed building blocks.

Design of complex alpha-beta protein topologies poses a challenge because of the large number of alternative packing arrangements. A similar challenge presumably limited the emergence of large and complex protein topologies in evolution. Here, we demonstrate that protein topologies with six and seven-stranded beta sheets can be designed by insertion of one de novo designed beta sheet containing protein into another such that the two beta sheets are merged to form a single extended sheet, followed by amino acid sequence optimization at the newly formed strand-strand, strand-helix, and helix-helix interfaces. Crystal structures of two such designs closely match the computational design models. Searches for similar structures in the SCOP protein domain database yield only weak matches with different beta sheet connectivities. A similar beta sheet fusion mechanism may have contributed to the emergence of complex beta sheets during natural protein evolution.

A synthetic sex ratio distortion system for the control of the human malaria mosquito.

It has been theorized that inducing extreme reproductive sex ratios could be a method to suppress or eliminate pest populations. Limited knowledge about the genetic makeup and mode of action of naturally occurring sex distorters and the prevalence of co-evolving suppressors has hampered their use for control. Here we generate a synthetic sex distortion system by exploiting the specificity of the homing endonuclease I-PpoI, which is able to selectively cleave ribosomal gene sequences of the malaria vector Anopheles gambiae that are located exclusively on the mosquito's X chromosome. We combine structure-based protein engineering and molecular genetics to restrict the activity of the potentially toxic endonuclease to spermatogenesis. Shredding of the paternal X chromosome prevents it from being transmitted to the next generation, resulting in fully fertile mosquito strains that produce >95% male offspring. We demonstrate that distorter male mosquitoes can efficiently suppress caged wild-type mosquito populations, providing the foundation for a new class of genetic vector control strategies.