In Vitro Directed Evolution of a Restriction Endonuclease with More Stringent Specificity.

Restriction endonuclease (REase) specificity engineering is extremely difficult. Here we describe a multistep protocol that helps to produce REase variants that have more stringent specificity than the parental enzyme. The protocol requires the creation of a library of expression selection cassettes (ESCs) for variants of the REase, ideally with variability in positions likely to affect DNA binding. The ESC is flanked on one side by a sequence for the restriction site activity desired and a biotin tag and on the other side by a restriction site for the undesired activity and a primer annealing site. The ESCs are transcribed and translated in a water-in-oil emulsion, in conditions that make the presence of more than one DNA molecule per droplet unlikely. Therefore, the DNA in each cassette molecule is subjected only to the activity of the translated, encoded enzyme. REase variants of the desired specificity remove the biotin tag but not the primer annealing site. After breaking the emulsion, the DNA molecules are subjected to a biotin pulldown, and only those in the supernatant are retained. This step assures that only ESCs for variants that have not lost the desired activity are retained. These DNA molecules are then subjected to a first PCR reaction. Cleavage in the undesired sequence cuts off the primer binding site for one of the primers. Therefore, PCR amplifies only ESCs from droplets without the undesired activity. A second PCR reaction is then carried out to reintroduce the restriction site for the desired specificity and the biotin tag, so that the selection step can be reiterated. Selected open reading frames can be overexpressed in bacterial cells that also express the cognate methyltransferase of the parental REase, because the newly evolved REase targets only a subset of the methyltransferase target sites.

Crystal Structure and Directed Evolution of Specificity of NlaIV Restriction Endonuclease.

Specificity engineering is challenging and particularly difficult for enzymes that have the catalytic machinery and specificity determinants in close proximity. Restriction endonucleases have been used as a paradigm for protein engineering, but successful cases are rare. Here, we present the results of a directed evolution approach to the engineering of a dimeric, blunt end cutting restriction enzyme NlaIV (GGN/NCC). Based on the remote similarity to EcoRV endonuclease, regions for random mutagenesis and in vitro evolution were chosen. The obtained variants cleaved target sites with an up to 100-fold kcat/KM preference for AT or TA (GGW/WCC) over GC or CG (GGS/SCC) in the central dinucleotide step, compared to the only ~17-fold preference of the wild-type enzyme. To understand the basis of the increased specificity, we determined the crystal structure of NlaIV. Despite the presence of DNA in the crystallization mix, the enzyme crystallized in the free form. We therefore constructed a computational model of the NlaIV-DNA complex. According to the model, the mutagenesis of the regions that were in the proximity of DNA did not lead to the desired specificity change, which was instead conveyed in an indirect manner by substitutions in the more distant regions.

Identification of protein structural elements responsible for the diversity of sequence preferences among Mini-III RNases.

Many known endoribonucleases select their substrates based on the presence of one or a few specific nucleotides at or near the cleavage site. In some cases, selectivity is also determined by the structural features of the substrate. We recently described the sequence-specific cleavage of double-stranded RNA by Mini-III RNase from Bacillus subtilis in vitro. Here, we characterized the sequence specificity of eight other members of the Mini-III RNase family from different bacterial species. High-throughput analysis of the cleavage products of Φ6 bacteriophage dsRNA indicated subtle differences in sequence preference between these RNases, which were confirmed and characterized by systematic analysis of the cleavage kinetics of a set of short dsRNA substrates. We also showed that the sequence specificities of Mini-III RNases are not reflected by different binding affinities for cognate and non-cognate sequences, suggesting that target selection occurs predominantly at the cleavage step. We were able to identify two structural elements, the α4 helix and α5b-α6 loop that were involved in target selection. Characterization of the sequence specificity of the eight Mini-III RNases may provide a basis for better understanding RNA substrate recognition by Mini-III RNases and adopting these enzymes and their engineered derivatives as tools for RNA research.

Sequence-specific endoribonucleases.

Ribonucleases are nucleolytic enzymes that commonly occur in living organisms and act by cleaving RNA molecules. These enzymes are involved in basic cellular processes, including the RNA maturation that accompanies the formation of functional RNAs, as well as RNA degradation that enables removal of defective or dangerous molecules or ones that have already fulfilled their cellular functions. RNA degradation is also one of the main processes that determine the amount of transcripts in the cell and thus it makes an important element of the gene expression regulation system. Ribonucleases can catalyse reactions involving RNA molecules containing specific sequences, structures or sequences within a specific structure, they can also cut RNAs non-specifically. In this article, we discuss ribonucleases cleaving the phosphodiester bond inside RNA molecules within or close to particular sequences. We also present examples of protein engineering of ribonucleases towards the development of molecular tools for sequence-specific cleavage of RNA.

Sequence-specific cleavage of dsRNA by Mini-III RNase.

Ribonucleases (RNases) play a critical role in RNA processing and degradation by hydrolyzing phosphodiester bonds (exo- or endonucleolytically). Many RNases that cut RNA internally exhibit substrate specificity, but their target sites are usually limited to one or a few specific nucleotides in single-stranded RNA and often in a context of a particular three-dimensional structure of the substrate. Thus far, no RNase counterparts of restriction enzymes have been identified which could cleave double-stranded RNA (dsRNA) in a sequence-specific manner. Here, we present evidence for a sequence-dependent cleavage of long dsRNA by RNase Mini-III from Bacillus subtilis (BsMiniIII). Analysis of the sites cleaved by this enzyme in limited digest of bacteriophage Φ6 dsRNA led to the identification of a consensus target sequence. We defined nucleotide residues within the preferred cleavage site that affected the efficiency of the cleavage and were essential for the discrimination of cleavable versus non-cleavable dsRNA sequences. We have also determined that the loop α5b-α6, a distinctive structural element in Mini-III RNases, is crucial for the specific cleavage, but not for dsRNA binding. Our results suggest that BsMiniIII may serve as a prototype of a sequence-specific dsRNase that could possibly be used for targeted cleavage of dsRNA.

Sequence-specific cleavage of the RNA strand in DNA-RNA hybrids by the fusion of ribonuclease H with a zinc finger.

Ribonucleases (RNases) are valuable tools applied in the analysis of RNA sequence, structure and function. Their substrate specificity is limited to recognition of single bases or distinct secondary structures in the substrate. Currently, there are no RNases available for purely sequence-dependent fragmentation of RNA. Here, we report the development of a new enzyme that cleaves the RNA strand in DNA-RNA hybrids 5 nt from a nonanucleotide recognition sequence. The enzyme was constructed by fusing two functionally independent domains, a RNase HI, that hydrolyzes RNA in DNA-RNA hybrids in processive and sequence-independent manner, and a zinc finger that recognizes a sequence in DNA-RNA hybrids. The optimization of the fusion enzyme's specificity was guided by a structural model of the protein-substrate complex and involved a number of steps, including site-directed mutagenesis of the RNase moiety and optimization of the interdomain linker length. Methods for engineering zinc finger domains with new sequence specificities are readily available, making it feasible to acquire a library of RNases that recognize and cleave a variety of sequences, much like the commercially available assortment of restriction enzymes. Potentially, zinc finger-RNase HI fusions may, in addition to in vitro applications, be used in vivo for targeted RNA degradation.

Delineation of structural domains and identification of functionally important residues in DNA repair enzyme exonuclease VII.

Exonuclease VII (ExoVII) is a bacterial nuclease involved in DNA repair and recombination that hydrolyses single-stranded DNA. ExoVII is composed of two subunits: large XseA and small XseB. Thus far, little was known about the molecular structure of ExoVII, the interactions between XseA and XseB, the architecture of the nuclease active site or its mechanism of action. We used bioinformatics methods to predict the structure of XseA, which revealed four domains: an N-terminal OB-fold domain, a middle putatively catalytic domain, a coiled-coil domain and a short C-terminal segment. By series of deletion and site-directed mutagenesis experiments on XseA from Escherichia coli, we determined that the OB-fold domain is responsible for DNA binding, the coiled-coil domain is involved in binding multiple copies of the XseB subunit and residues D155, R205, H238 and D241 of the middle domain are important for the catalytic activity but not for DNA binding. Altogether, we propose a model of sequence-structure-function relationships in ExoVII.

Tau inhibits tubulin oligomerization induced by prion protein.

In previous studies we have demonstrated that prion protein (PrP) interacts with tubulin and disrupts microtubular cytoskeleton by inducing tubulin oligomerization. These observations may explain the molecular mechanism of toxicity of cytoplasmic PrP in transmissible spongiform encephalopathies (TSEs). Here, we check whether microtubule associated proteins (MAPs) that regulate microtubule stability, influence the PrP-induced oligomerization of tubulin. We show that tubulin preparations depleted of MAPs are more prone to oligomerization by PrP than those containing traces of MAPs. Tau protein, a major neuronal member of the MAPs family, reduces the effect of PrP. Importantly, phosphorylation of Tau abolishes its ability to affect the PrP-induced oligomerization of tubulin. We propose that the binding of Tau stabilizes tubulin in a conformation less susceptible to oligomerization by PrP. Since elevated phosphorylation of Tau leading to a loss of its function is observed in Alzheimer disease and related tauopathies, our results point at a possible molecular link between these neurodegenerative disorders and TSEs.

Interactions between subunits in heterodimeric Ncd molecules.

The nonprocessive minus-end-directed kinesin-14 Ncd is involved in the organization of the microtubule (MT) network during mitosis. Only one of the two motor domains is involved in the interaction with the MT. The other head is tethered to the bound one. Here we prepared, purified, and characterized mutated Ncd molecules carrying point mutations in one of the heads, thus producing heterodimeric motors. The mutations tested included substitutions in Switch I and II: R552A, E585A, and E585D; the decoupling mutant N600K; and a deletion in the motor domain in one of the subunits resulting in a single-headed molecule (NcN). These proteins were isolated by two sequential affinity chromatography steps, followed by measurements of their affinities to MT, enzymatic properties, and the velocity of the microtubule gliding test in vitro. A striking observation is a low affinity of the single-headed NcN for MT both without nucleotides and in the presence of 5'-adenylyl-beta,gamma-imidodiphosphate, implying that the tethered head has a profound effect on the structure of the Ncd-MT complex. Mutated homodimers had no MT-activated ATPase and no motility, whereas NcN had motility comparable with that of the wild type Ncd. Although the heterodimers had one fully active and one inactive head, the ATPase and motility of Ncd heterodimers varied dramatically, clearly demonstrating that interactions between motor domains exist in Ncd. We also show that the bulk property of dimeric proteins that interact with the filament with only one of its heads depends also on the distribution of the filament-interacting subunits.

Prion protein region 23-32 interacts with tubulin and inhibits microtubule assembly.

In previous studies we have demonstrated that prion protein (PrP) binds directly to tubulin and this interaction leads to the inhibition of microtubule formation by inducement of tubulin oligomerization. This report is aimed at mapping the regions of PrP and tubulin involved in the interaction and identification of PrP domains responsible for tubulin oligomerization. Preliminary studies focused our attention to the N-terminal flexible part of PrP encompassing residues 23-110. Using a panel of deletion mutants of PrP, we identified two microtubule-binding motifs at both ends of this part of the molecule. We found that residues 23-32 constitute a major site of interaction, whereas residues 101-110 represent a weak binding site. The crucial role of the 23-32 sequence in the interaction with tubulin was confirmed employing chymotryptic fragments of PrP. Surprisingly, the octarepeat region linking the above motifs plays only a supporting role in the interaction. The binding of Cu(2+) to PrP did not affect the interaction. We also demonstrate that PrP deletion mutants lacking residues 23-32 exhibit very low efficiency in the inducement of tubulin oligomerization. Moreover, a synthetic peptide corresponding to this sequence, but not that identical with fragment 101-110, mimics the effects of the full-length protein on tubulin oligomerization and microtubule assembly. At the cellular level, peptide composed of the PrP motive 23-30 and signal sequence (1-22) disrupted the microtubular cytoskeleton. Using tryptic and chymotryptic fragments of alpha- and beta-tubulin, we mapped the docking sites for PrP within the C-terminal domains constituting the outer surface of microtubule.

The mitotic kinesin-14 Ncd drives directional microtubule-microtubule sliding.

During mitosis and meiosis, the bipolar spindle facilitates chromosome segregation through microtubule sliding as well as microtubule growth and shrinkage. Kinesin-14, one of the motors involved, causes spindle collapse in the absence of kinesin-5 (Refs 2, 3), participates in spindle assembly and modulates spindle length. However, the molecular mechanisms underlying these activities are not known. Here, we report that Drosophila melanogaster kinesin-14 (Ncd) alone causes sliding of anti-parallel microtubules but locks together (that is, statically crosslinks) those that are parallel. Using single molecule imaging we show that Ncd diffuses along microtubules in a tail-dependent manner and switches its orientation between sliding microtubules. Our results show that kinesin-14 causes sliding and expansion of an anti-parallel microtubule array by dynamic interactions through the motor domain on the one side and the tail domain on the other. This mechanism accounts for the roles of kinesin-14 in spindle organization.

Subunits interactions in kinesin motors.

Kinesins form a large and diverse superfamily of proteins involved in numerous important cellular processes. The majority of them are molecular motors moving along microtubules. Conversion of chemical energy into mechanical work is accomplished in a sequence of events involving both biochemical and conformational alternation of the motor structure called the mechanochemical cycle. Different members of the kinesin superfamily can either perform their function in large groups or act as single molecules. Conventional kinesin, a member of the kinesin-1 subfamily, exemplifies the second type of motor which requires tight coordination of the mechanochemical cycle in two identical subunits to accomplish processive movement toward the microtubule plus end. Recent results strongly support an asymmetric hand-over-hand model of "walking" for this protein. Conformational strain between two subunits at the stage of the cycle where both heads are attached to the microtubule seems to be a major factor in intersubunit coordination, although molecular and kinetic details of this phenomenon are not yet deciphered. We discuss also current knowledge concerning intersubunit coordination in other kinesin subfamilies. Members of the kinesin-3 class use at least three different mechanisms of movement and can translocate in monomeric or dimeric forms. It is not known to what extent intersubunit coordination takes place in Ncd, a dimeric member of the kinesin-14 subfamily which, unlike conventional kinesin, exercises a power-stroke toward the microtubule minus end. Eg5, a member of the kinesin-5 subfamily is a homotetrameric protein with two kinesin-1-like dimeric halves controlled by their relative orientation on two microtubules. It seems that diversity of subunit organization, quaternary structures and cellular functions in the kinesin superfamily are reflected also by the divergent extent and mechanism of intersubunit coordination during kinesin movement along microtubules.

Identification of a new subfamily of HNH nucleases and experimental characterization of a representative member, HphI restriction endonuclease.

The restriction endonuclease (REase) R. HphI is a Type IIS enzyme that recognizes the asymmetric target DNA sequence 5'-GGTGA-3' and in the presence of Mg(2+) hydrolyzes phosphodiester bonds in both strands of the DNA at a distance of 8 nucleotides towards the 3' side of the target, producing a 1 nucleotide 3'-staggered cut in an unspecified sequence at this position. REases are typically ORFans that exhibit little similarity to each other and to any proteins in the database. However, bioinformatics analyses revealed that R.HphI is a member of a relatively big sequence family with a conserved C-terminal domain and a variable N-terminal domain. We predict that the C-terminal domains of proteins from this family correspond to the nuclease domain of the HNH superfamily rather than to the most common PD-(D/E)XK superfamily of nucleases. We constructed a three-dimensional model of the R.HphI catalytic domain and validated our predictions by site-directed mutagenesis and studies of DNA-binding and catalytic activities of the mutant proteins. We also analyzed the genomic neighborhood of R.HphI homologs and found that putative nucleases accompanied by a DNA methyltransferase (i.e. predicted REases) do not form a single group on a phylogenetic tree, but are dispersed among free-standing putative nucleases. This suggests that nucleases from the HNH superfamily were independently recruited to become REases in the context of RM systems multiple times in the evolution and that members of the HNH superfamily may be much more frequent among the so far unassigned REase sequences than previously thought.

Theoretical model of restriction endonuclease HpaI in complex with DNA, predicted by fold recognition and validated by site-directed mutagenesis.

Type II restriction enzymes are commercially important deoxyribonucleases and very attractive targets for protein engineering of new specificities. At the same time they are a very challenging test bed for protein structure prediction methods. Typically, enzymes that recognize different sequences show little or no amino acid sequence similarity to each other and to other proteins. Based on crystallographic analyses that revealed the same PD-(D/E)XK fold for more than a dozen case studies, they were nevertheless considered to be related until the combination of bioinformatics and mutational analyses has demonstrated that some of these proteins belong to other, unrelated folds PLD, HNH, and GIY-YIG. As a part of a large-scale project aiming at identification of a three-dimensional fold for all type II REases with known sequences (currently approximately 1000 proteins), we carried out preliminary structure prediction and selected candidates for experimental validation. Here, we present the analysis of HpaI REase, an ORFan with no detectable homologs, for which we detected a structural template by protein fold recognition, constructed a model using the FRankenstein monster approach and identified a number of residues important for the DNA binding and catalysis. These predictions were confirmed by site-directed mutagenesis and in vitro analysis of the mutant proteins. The experimentally validated model of HpaI will serve as a low-resolution structural platform for evolutionary considerations in the subgroup of blunt-cutting REases with different specificities. The research protocol developed in the course of this work represents a streamlined version of the previously used techniques and can be used in a high-throughput fashion to build and validate models for other enzymes, especially ORFans that exhibit no sequence similarity to any other protein in the database.

Direct interaction between prion protein and tubulin.

Recently published data show that the prion protein in its cellular form (PrP(C)) is a component of multimolecular complexes. In this report, zero-length cross-linking with 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC) allowed us to identify tubulin as one of the molecules interacting with PrP(C) in complexes observed in porcine brain extracts. We found that porcine brain tubulin added to these extracts can be cross-linked with PrP(C). Moreover, we observed that the 34 kDa species identified previously as full-length diglycosylated prion protein co-purifies with tubulin. Cross-linking of PrP(C) species separated by Cu(2+)-loaded immobilized metal affinity chromatography confirmed that only the full-length protein but not the N-terminally truncated form (C1) binds to tubulin. By means of EDC cross-linking and cosedimentation experiments, we also demonstrated a direct interaction of recombinant human PrP (rPrP) with tubulin. The stoichiometry of cosedimentation implies that rPrP molecules are able to bind both the alpha- and beta-isoforms of tubulin composing microtubule. Furthermore, prion protein exhibits higher affinity for microtubules than for unpolymerized tubulin.

A theoretical model of restriction endonuclease NlaIV in complex with DNA, predicted by fold recognition and validated by site-directed mutagenesis and circular dichroism spectroscopy.

Restriction enzymes (REases) are commercial reagents commonly used in DNA manipulations and mapping. They are regarded as very attractive models for studying protein-DNA interactions and valuable targets for protein engineering. Their amino acid sequences usually show no similarities to other proteins, with rare exceptions of other REases that recognize identical or very similar sequences. Hence, they are extremely hard targets for structure prediction and modeling. NlaIV is a Type II REase, which recognizes the interrupted palindromic sequence GGNNCC (where N indicates any base) and cleaves it in the middle, leaving blunt ends. NlaIV shows no sequence similarity to other proteins and virtually nothing is known about its sequence-structure-function relationships. Using protein fold recognition, we identified a remote relationship between NlaIV and EcoRV, an extensively studied REase, which recognizes the GATATC sequence and whose crystal structure has been determined. Using the 'FRankenstein's monster' approach we constructed a comparative model of NlaIV based on the EcoRV template and used it to predict the catalytic and DNA-binding residues. The model was validated by site-directed mutagenesis and analysis of the activity of the mutants in vivo and in vitro as well as structural characterization of the wild-type enzyme and two mutants by circular dichroism spectroscopy. The structural model of the NlaIV-DNA complex suggests regions of the protein sequence that may interact with the 'non-specific' bases of the target and thus it provides insight into the evolution of sequence specificity in restriction enzymes and may help engineer REases with novel specificities. Before this analysis was carried out, neither the three-dimensional fold of NlaIV, its evolutionary relationships or its catalytic or DNA-binding residues were known. Hence our analysis may be regarded as a paradigm for studies aiming at reducing 'white spaces' on the evolutionary landscape of sequence-function relationships by combining bioinformatics with simple experimental assays.

Inference of relationships in the 'twilight zone' of homology using a combination of bioinformatics and site-directed mutagenesis: a case study of restriction endonucleases Bsp6I and PvuII.

Thus far, identification of functionally important residues in Type II restriction endonucleases (REases) has been difficult using conventional methods. Even though known REase structures share a fold and marginally recognizable active site, the overall sequence similarities are statistically insignificant, unless compared among proteins that recognize identical or very similar sequences. Bsp6I is a Type II REase, which recognizes the palindromic DNA sequence 5'GCNGC and cleaves between the cytosine and the unspecified nucleotide in both strands, generating a double-strand break with 5'-protruding single nucleotides. There are no solved structures of REases that recognize similar DNA targets or generate cleavage products with similar characteristics. In straightforward comparisons, the Bsp6I sequence shows no significant similarity to REases with known structures. However, using a fold-recognition approach, we have identified a remote relationship between Bsp6I and the structure of PvuII. Starting from the sequence-structure alignment between Bsp6I and PvuII, we constructed a homology model of Bsp6I and used it to predict functionally significant regions in Bsp6I. The homology model was supported by site-directed mutagenesis of residues predicted to be important for dimerization, DNA binding and catalysis. Completing the picture of sequence-structure-function relationships in protein superfamilies becomes an essential task in the age of structural genomics and our study may serve as a paradigm for future analyses of superfamilies comprising strongly diverged members with little or no sequence similarity.

A homology model of restriction endonuclease SfiI in complex with DNA.

BACKGROUND: Restriction enzymes (REases) are commercial reagents commonly used in recombinant DNA technologies. They are attractive models for studying protein-DNA interactions and valuable targets for protein engineering. They are, however, extremely divergent: the amino acid sequence of a typical REase usually shows no detectable similarities to any other proteins, with rare exceptions of other REases that recognize identical or very similar sequences. From structural analyses and bioinformatics studies it has been learned that some REases belong to at least four unrelated and structurally distinct superfamilies of nucleases, PD-DxK, PLD, HNH, and GIY-YIG. Hence, they are extremely hard targets for structure prediction and homology-based inference of sequence-function relationships and the great majority of REases remain structurally and evolutionarily unclassified. RESULTS: SfiI is a REase which recognizes the interrupted palindromic sequence 5'GGCCNNNN--NGGCC3' and generates 3 nt long 3' overhangs upon cleavage. SfiI is an archetypal Type IIF enzyme, which functions as a tetramer and cleaves two copies of the recognition site in a concerted manner. Its sequence shows no similarity to other proteins and nothing is known about the localization of its active site or residues important for oligomerization. Using the threading approach for protein fold-recognition, we identified a remote relationship between SfiI and BglI, a dimeric Type IIP restriction enzyme from the PD-DxK superfamily of nucleases, which recognizes the 5'GCCNNNN--NGGC3' sequence and whose structure in complex with the substrate DNA is available. We constructed a homology model of SfiI in complex with its target sequence and used it to predict residues important for dimerization, tetramerization, DNA binding and catalysis. CONCLUSIONS: The bioinformatics analysis suggest that SfiI, a Type IIF enzyme, is more closely related to BglI, an "orthodox" Type IIP restriction enzyme, than to any other REase, including other Type IIF REases with known structures, such as NgoMIV. NgoMIV and BglI belong to two different, very remotely related branches of the PD-DxK superfamily: the alpha-class (EcoRI-like), and the beta-class (EcoRV-like), respectively. Thus, our analysis provides evidence that the ability to tetramerize and cut the two DNA sequences in a concerted manner was developed independently at least two times in the evolution of the PD-DxK superfamily of REases. The model of SfiI will also serve as a convenient platform for further experimental analyses.

Laboratory scale structural genomics.

At Lawrence Livermore National Laboratory, the development of the TB structural genomics consortium crystallization facility has paralleled several local proteomics research efforts that have grown out of gene expression microarray and comparative genomics studies. Collective experience gathered from TB consortium labs and other centers involved in the NIH-NIGMS protein structure initiative allows us to explore the possibilities and challenges of pursuing structural genomics on an academic laboratory scale. We discuss our procedures and protocols for genomic targeting approaches, primer design, cloning, small scale expression screening, scale-up and purification, through to automated crystallization screening and data collection. The procedures are carried out by a small group using a combination of traditional approaches, innovative molecular biochemistry approaches, software automation, and a modest investment in robotic equipment.