Characterization of norovirus RNA replicase for in vitro amplification of RNA.

BACKGROUND: The isothermal amplification of RNA in vitro has been used for the study of in vitro evolution of RNA. Although Qß replicase has been traditionally used as an enzyme for this purpose, we planned to use norovirus replicase (NV3D(pol)) due to its structural simplicity in the scope of in vitro autonomous evolution of the protein. Characteristics of the enzyme NV3D(pol) in vitro were re-evaluated in this context. RESULTS: NV3D(pol), synthesized by using a cell-free translation system, represented the activities which were reported in the previous several studies and the reports were not fully consistent each other. The efficiency of the initiation of replication was dependent on the 3'-terminal structure of single-stranded RNA template, and especially, NV3D(pol) preferred a self-priming small stem-loop. In the non-self-priming and primer-independent replication reaction, the presence of -CCC residues at the 3'-terminus increased the initiation efficiency and we demonstrated the one-pot isothermal RNA (even dsRNA) amplification by 16-fold. NV3D(pol) also showed a weak activity of elongation-reaction from a long primer. Based on these results, we present a scheme of the primer-independent isothermal amplification of RNA with NV3D(pol) in vitro. CONCLUSIONS: NV3D(pol) can be used as an RNA replicase in in vitro RNA + protein evolution with the RNA of special terminal sequences.

Biophysical connection between evolutionary dynamics and thermodynamics in in vitro evolution.

We analyzed a mathematical model of in vitro evolution conducted by repetition of mutagenesis and selection processes. The selection process consists of the selective enrichment and subsequent sampling as follows: each mutant with fitness W is amplified by the Boltzmann factor exp(rW/k(B)T(the)), where the fitness W is defined as the negative Gibbs free energy (-ΔG) in a reaction of the phenotypic molecules and r is the round number of the selective enrichment; then, an arbitrary mutant is randomly chosen from the resulting mutant population and it becomes a new parent in the next generation. As a result, we found that the evolutionary dynamics is described in a mathematical framework similar to thermodynamics: the "evolution constant" k(E) and "evolutionary temperature" T(evo) play key roles similar to the Boltzmann constant k(B) and thermodynamic temperature T(the), respectively. In the stationary state of the evolutionary dynamics, the attractor of the fitness is in inverse proportion to k(E)T(evo). Furthermore, beyond the mathematical analogy, we obtained a biophysical connection between evolutionary dynamics and thermodynamics. Particularly, we found that T(evo) and T(the) are connected by k(E)T(evo)≈k(B)T(the)/2r. These results suggest that we can predict the fitness value in the stationary state by the thermodynamic temperature T(the) in the experimental setup.

Toward an evolutionary containment of evolving pathogen-receptors by using an ensemble of multiple mutant ligands: from the viewpoint of fitness landscape in sequence space.

It is known that even if a ligand peptide is designed to bind to a target receptor on the surface of a pathogen such as viruses, bacteria or cancer cells, it is likely that some receptors are subject to random mutation and thus the ligand has a reduced ability to bind to these receptors. This issue is known as drug-resistant or escape mutants. In this paper, we present an idea to inhibit the evolving receptors by using an ensemble of all possible single- or double-point mutant sequences of the ligand peptide. Several mutant ligands in the ensemble are expected to bind to the mutant receptors, and then the ensemble may create a defensive wall surrounding the target receptors in receptor-sequence space. We examined the effectiveness of this "evolutionary containment" of the evolving receptors through eight peptide-protein complex systems, which were retrieved from the Protein Data Bank (PDB). As a result, we obtained a suggestion that the original (or parent) ligand sequence should be designed to have as high fitness as possible but to be not local optima, in order to maximize the rate of the evolutionary containment. This may be a strategy of the drug-design against evolving pathogens.

Gel shift selection of translation enhancer sequences using messenger RNA display.

We designed a new approach for selection of translation enhancer sequences that enables efficient protein synthesis in cell-free systems. The selection is based on a gel shift assay of a messenger RNA (mRNA)-protein fusion product that is synthesized in a cell-free translation system using an mRNA display method. A library of randomized 20-nt-long sequences, with all possible combinations of the four nucleotides, upstream of a coding region was screened by successive rounds of screening in which the translation time of the succeeding round was reduced compared with the previous round. An efficient translation enhancer sequence capable of more rapid initiation of cell-free protein synthesis, with a minimal translation time of 5 min, than a natural longer enhancer sequence (Xenopus ß-globin 5'UTR) was selected using rabbit reticulocyte extract as a model cell-free translation system. Furthermore, a successful screening of cap-independent translation enhancer sequence and a significant sequence similarity of the selected candidates validated the efficiency of the combined mRNA display and gel shift assay method for the rapid development of advanced cell-free translation systems.

cDNA display: a novel screening method for functional disulfide-rich peptides by solid-phase synthesis and stabilization of mRNA-protein fusions.

We report a robust display technology for the screening of disulfide-rich peptides, based on cDNA-protein fusions, by developing a novel and versatile puromycin-linker DNA. This linker comprises four major portions: a 'ligation site' for T4 RNA ligase, a 'biotin site' for solid-phase handling, a 'reverse transcription primer site' for the efficient and rapid conversion from an unstable mRNA-protein fusion (mRNA display) to a stable mRNA/cDNA-protein fusion (cDNA display) whose cDNA is covalently linked to its encoded protein and a 'restriction enzyme site' for the release of a complex from the solid support. This enables not only stabilizing mRNA-protein fusions but also promoting both protein folding and disulfide shuffling reactions. We evaluated the performance of cDNA display in different model systems and demonstrated an enrichment efficiency of 20-fold per selection round. Selection of a 32-residue random library against interleukin-6 receptor generated novel peptides containing multiple disulfide bonds with a unique linkage for its function. The peptides were found to bind with the target in the low nanomolar range. These results show the suitability of our method for in vitro selections of disulfide-rich proteins and other potential applications.

Novel concept microarray enabling PCR and multistep reactions through pipette-free aperture-to-aperture parallel transfer.

BACKGROUND: The microarray has contributed to developing the omic analysis. However, as it depends basically on the surface reaction, it is hard to perform bulk reactions and sequential multistep reactions. On the other hand, the popular microplate technology, which has a great merit of being able to perform parallel multistep reactions, has come to its limit in increasing the number of wells (currently, up to 9600) and reducing the volume to deal with due to the difficulty in operations. RESULTS: Here, we report a novel microarray technology which enables us to explore advanced applications, termed microarray-with-manageable volumes (MMV). The technical essence is in the pipette-free direct parallel transfer from well to well performed by centrifugation, evading the evaporation and adsorption-losses during handling. By developing the MMV plate, accompanying devices and techniques, generation of multiple conditions (256 kinds) and performance of parallel multistep reactions, including PCR and in vitro translation reactions, have been made possible. These were demonstrated by applying the MMV technology to searching lysozyme-crystallizing conditions and selecting peptides aimed for Aß-binding or cathepsin E-inhibition. CONCLUSIONS: With the introduction of a novel concept microarray (MMV) technology, parallel and multistep reactions in sub-µL scale have become possible.

Quantum dot FRET biosensors that respond to pH, to proteolytic or nucleolytic cleavage, to DNA synthesis, or to a multiplexing combination.

Fluorescent acceptors have been immobilized on nanoparticulate quantum dots (QDs), which serve in turn as their FRET donors. The broad excitation and narrow emission bands of QDs mark them as having excellent potential as donors for FRET and, in principle, differently colored QDs could be excited simultaneously. The present work describes the preparation and operation of FRET-based QD bioprobes individually able to detect the actions of protease, deoxyribonuclease, DNA polymerase, or changes in pH. In addition, two such QD-mounted biosensors were excited at a single wavelength, and shown to operate simultaneously and independently of each other in the same sample solution, allowing multiplex detection of the action of a protease, trypsin, in the presence of deoxyribonuclease.

Fitting protein-folding free energy landscape for a certain conformation to an NK fitness landscape.

The NK fitness landscape is a mathematical landscape model with a parameter k that governs the degree of ruggedness of the landscape. We presented a procedure to fit a given landscape to the NK fitness landscape by introducing the "apparent k-value"k(app). In this paper, we defined the protein free energy (DeltaG) landscape in amino acid sequence space, where DeltaG is the folding free energy from a random coil to a "certain conformation". Applying this landscape to our fitting procedure, we examined the statistical properties of the landscape. For calculation of a conformation energy, amino acid residues are represented by points, and interaction energies among amino acid residues are given as (1+K)-body interactions, that is, an unit of interacting (1+K) residues cooperatively contribute a single energy value to the conformational energy. Our results suggest that the apparent k-value of the free energy landscape is k(app) approximately K, and that the number of possible interactions among residues is unrelated to the k(app) value. This leads to the inference that k(app) takes values about 1-3 in real landscapes, if nature adopts two-body approximately four-body interaction energies.

Solid-phase translation and RNA-protein fusion: a novel approach for folding quality control and direct immobilization of proteins using anchored mRNA.

A novel cell-free translation system is described in which template-mRNA molecules were captured onto solid surfaces to simultaneously synthesize and immobilize proteins in a more native-state form. This technology comprises a novel solid-phase approach to cell-free translation and RNA-protein fusion techniques. A newly constructed biotinylated linker-DNA which enables puromycin-assisted RNA-protein fusion is ligated to the 3' ends of the mRNA molecules to attach the mRNA-template on a streptavidin-coated surface and further to enable the subsequent reactions of translation and RNA-protein fusion on surface. The protein products are therefore directly immobilized onto solid surfaces and furthermore were discovered to adopt a more native state with proper protein folding and superior biological activity compared with conventional liquid-phase approaches. We further validate this approach via the production of immobilized green fluorescent protein (GFP) on microbeads and by the production and assay of aldehyde reductase (ALR) enzyme with 4-fold or more activity. The approach developed in this study may enable to embrace the concept of the transformation of 'RNA chip-to-protein chip' using a solid-phase cell-free translation system and thus to the development of high-throughput microarray platform in the field of functional genomics and in vitro evolution.

Directed evolution by accumulating tailored mutations: thermostabilization of lactate oxidase with less trade-off with catalytic activity.

We assumed that adverse effects posed by introducing multiple mutations could be decomposed into those of each of the component mutations and that the risk could be reduced by the accumulation of mutations that were finely tuned for directed improvement of a specific property. We propose here a directed evolution strategy for improving a specific property with less effect on other ones. This strategy is composed of fine-tuning of mutations and their accumulation by our original mutation-assembling method. In this study, we selected lactate oxidase (LOX) as a model enzyme, because its directed evolution had showed a trade-off between thermostability and catalytic activity. Mutation profiling at each of the sites found by error-prone PCR revealed a strong inverse relationship between the two properties. Thermostable mutations with less effect on catalytic activity were selected at each site and accumulated with ideal combinations by our method. The resultant multiple mutants exhibited 5- to 10-fold superior catalytic activity and comparable thermostability with those created by accumulating thermostable mutations, which were not tuned for catalytic activity. This result demonstrates that the accumulation of fine-tuned mutations is an advantageous approach to reduce the risk of adverse effects posed by accumulating multiple mutations.

An RNA Binding Peptide Consisting of Four Types of Amino Acid by in Vitro Selection Using cDNA Display.

RNA-protein interactions have a central role in the living world. In this article, we examined whether primitive peptides (30 residues) consisting of four types of amino acid (Gly, Ala, Asp, and Val) could interact with tRNA as a model of primitive RNAs in the RNA world. By in vitro selection of binding peptides using the cDNA display method, a characteristic peptide was selected from a random peptide library and assayed by electrophoretic mobility shift and pull-down assays. Interestingly, the selected peptide bound to a single-stranded region including a loop structure of an RNA molecule with some sequence specificity.

Protease-sensitive signalling by chemically engineered intramolecular fluorescent resonance energy transfer mutants of green fluorescent protein.

The native cysteine residues of green fluorescent protein (GFP) at positions 48 and 70 were replaced by non-thiolic amino acids, and new cysteine sites were introduced at specific, surface positions. Based on molecular modeling of the GFP structure, the sites chosen for mutagenesis to Cys were glutamic acid at position 6 and isoleucine at position 229. These new, unique cysteine sites provided reactive thiol groups suitable for site-specific chemical modification by eosin-based fluorescence labels. The new constructs were designed to serve as the basis of proof of principle for fluorescence resonance energy transfer (FRET) using an enzyme-activated (trypsin) intervening sequence between native and chemically conjugated fluorophores. These eosin moieties provided chemical FRET partners for the native GFP chromophore. On excitation, these GFP-eosin constructs exhibited strong intramolecular FRET, with quenching of the native GFP (511 nm) fluorophore emission and emission around 540 nm, corresponding to eosin. GFP mutants engineered with trypsin-sensitive sequences close to the eosin site, so that on trypsinolysis FRET was destroyed, the emission wavelength switching from that of the chemical FRET partner back to that of the native GFP fluorophore, providing efficient, ratio-based detection. This protein engineering provides the basis for novel bioprobes for enzymatic triggering using intramolecular FRET between GFP and carefully sited chemical labels.

Biased mutation-assembling: an efficient method for rapid directed evolution through simultaneous mutation accumulation.

We have developed an efficient optimization technique, 'biased mutation-assembling', for improving protein properties such as thermostability. In this strategy, a mutant library is constructed using the overlap extension polymerase chain reaction technique with DNA fragments from wild-type and phenotypically advantageous mutant genes, in which the number of mutations assembled in the wild-type gene is stochastically controlled by the mixing ratio of the mutant DNA fragments to wild-type fragments. A high mixing ratio results in a mutant composition biased to favor multiple-point mutants. We applied this strategy to improve the thermostability of prolyl endopeptidase from Flavobacterium meningosepticum as a case study and found that the proportion of thermostable mutants in a library increased as the mixing ratio was increased. If the proportion of thermostable mutants increases, the screening effort needed to find them should be reduced. Indeed, we isolated a mutant with a 1200-fold longer activity half-life at 60 degrees C than that of wild-type prolyl endopeptidase after screening only 2000 mutants from a library prepared with a high mixing ratio. Our results indicate that an aggressive accumulation of advantageous mutations leads to an increase in the quality of the mutant library and a reduction in the screening effort required to find superior mutants.

Multi-line split DNA synthesis: a novel combinatorial method to make high quality peptide libraries.

BACKGROUND: We developed a method to make a various high quality random peptide libraries for evolutionary protein engineering based on a combinatorial DNA synthesis. RESULTS: A split synthesis in codon units was performed with mixtures of bases optimally designed by using a Genetic Algorithm program. It required only standard DNA synthetic reagents and standard DNA synthesizers in three lines. This multi-line split DNA synthesis (MLSDS) is simply realized by adding a mix-and-split process to normal DNA synthesis protocol. Superiority of MLSDS method over other methods was shown. We demonstrated the synthesis of oligonucleotide libraries with 1016 diversity, and the construction of a library with random sequence coding 120 amino acids containing few stop codons. CONCLUSIONS: Owing to the flexibility of the MLSDS method, it will be able to design various "rational" libraries by using bioinformatics databases.

An Efficient Ligation Method in the Making of an in vitro Virus for in vitro Protein Evolution.

The "in vitro virus" is a molecular construct to perform evolutionary protein engineering. The "virion (=viral particle)" (mRNA-peptide fusion), is made by bonding a nascent protein with its coding mRNA via puromycin in a test tube for in vitro translation. In this work, the puromycin-linker was attached to mRNA using the Y-ligation, which was a method of two single-strands ligation at the end of a double-stranded stem to make a stem-loop structure. This reaction gave a yield of about 95%. We compared the Y-ligation with two other ligation reactions and showed that the Y-ligation gave the best productivity. An efficient amplification of the in vitro virus with this "viral genome" was demonstrated.

Estimation of statistical binding properties of ligand population during in vitro selection based on population dynamics theory.

During in vitro selection process, it is very valuable to monitor the binding properties of the ligand population in real time, particularly the population average of the association constant in the population. If this monitoring can be realized, the selection process can be controlled in a rational way. In this paper, we present a simple method to estimate the binding properties of the ligand population during in vitro selection. The framework of the method is as follows. First, the number of all the collected ligand molecules, which are eluted after incubation and washing, is measured. Ideally, this number corresponds to the number of all the ligand molecules bound with the target-receptor or other materials in a test tube. This measurement is performed through several successive rounds of selection. Second, the measured numbers of molecules are subjected to a theoretical analysis, based on the mathematical theory of population dynamics in the selection process. Then, we can estimate the probability density of the binding free energy in the ligand population. The validity of our method was confirmed by several computer simulations based on a physicochemical model.

Optimal terminal sequences for continuous or serial isothermal amplification of dsRNA with norovirus RNA replicase.

The norovirus RNA replicase (NV3D(pol), 56 kDa, single chain monomeric protein) can amplify double-stranded (ds) RNA isothermally. It will play an alternative role in the in vitro evolution against traditional Qß RNA replicase, which cannot amplify dsRNA and consists of four subunits, three of which are borrowed from host E.coli. In order to identify the optimal 3'-terminal sequence of the RNA template for NV3D(pol), an in vitro selection using the serial transfer was performed for a random library having the 3'-terminal sequence of ---UUUUUUNNNN-3'. The population landscape on the 4-dimensional sequence space of the 17(th) round of transfer gave a main peak around ---CAAC-3'. In the preceding studies on the batch amplification reaction starting from a single-stranded RNA, a template with 3'-terminal C-stretch was amplified effectively. It was confirmed that in the batch amplification the ---CCC-3' was much more effective than the ---CAAC-3', but in the serial transfer condition in which the ----CAAC-3' was sustained stably, the ---CCC-3' was washed out. Based on these results we proposed the existence of the "shuttle mode" replication of dsRNA. We also proposed the optimal terminal sequences of RNA for in vitro evolution with NV3D(pol).

Theoretical consideration of selective enrichment in in vitro selection: optimal concentration of target molecules.

We considered an in vitro selection system composed of a peptide-ligand library and a single target protein receptor, and examined effective strategies to realize maximum efficiency in selection. In the system, a ligand molecule with sequence s binds to a target receptor with probability of [R]/(K(ds)+[R]) (specific binding) or binds to non-target materials with probability of q (non-specific binding), where [R] and K(ds) represent the free target-receptor concentration at equilibrium and dissociation constant K(d) of the ligand sequence s with the receptor, respectively. Focusing on the fittest sequence with the highest affinity (represented by K(d1) ≡ min{K(ds)|s=1,2,…,M}) in the ligand library with a library size N and diversity M, we examined how the target concentration [R] should be set in each round to realize the maximum enrichment of the fittest sequence. In conclusion, when N >> M (that realizes a deterministic process), it is desirable to adopt [R]=K(d1), and when N=M (that realizes a stochastic process), [R]=[Formula in text] only in the first round (where * represents the population average) and [R]=K(d1) in the subsequent rounds. Based on this strategy, the mole fraction of the fittest increases by (2q)(-r) times after the rth round. With realistic parameters, we calculated several quantities such as the optimal [R] values and number of rounds needed. These values were quite reasonable and consistent with observations, suggesting the validity of our theory.

A mathematical consideration of the word-composition vector method in comparison of biological sequences.

To measure the similarity or dissimilarity between two given biological sequences, several papers proposed metrics based on the "word-composition vector". The essence of these metrics is as follows. First, we count the appearance frequencies of all the K-tuple words throughout each of two given sequences. Then, the two given sequences are transformed into their respective word-composition vectors. Next, the distance metrics, for example the angle between the two vectors, are calculated. A significant issue is to determine the optimal word size K. With a mathematical model of mutational events (including substitutions, insertions, deletions and duplications) that occur in sequences, we analyzed how the angle between the composition vectors depends on the mutational events. We also considered the optimal word size (=resolution) from our original approach. Our results were verified by computational experiments using artificially generated sequences, amino acid sequences of hemoglobin and nucleotide sequences of 16S ribosomal RNA.

One-pot preparation of mRNA/cDNA display by a novel and versatile puromycin-linker DNA.

A rapid, easy, and robust preparation method for mRNA/cDNA display using a newly designed puromycin-linker DNA is presented. The new linker is structurally simple, easy to synthesize, and cost-effective for use in "in vitro peptide and protein selection". An introduction of RNase T1 nuclease site to the new linker facilitates the easy recovery of mRNA/cDNA displayed protein by an improvement of the efficiency of ligating the linker to mRNAs and efficient release of mRNA/cDNA displayed protein from the solid-phase (magnetic bead). For application demonstration, affinity selections were successfully performed. Furthermore, we introduced a "one-pot" preparation protocol to perform mRNA display easy. Unlike conventional approaches that require tedious and downstream multistep process including purification, this protocol will make the mRNA/cDNA display methods more practical and convenient and also facilitate the development of next-generation, high-throughput mRNA/cDNA display systems amenable to automation.