Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “directed evolution”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87Linked to original sources

Dikaryons of the basidiomycete fungus Schizophyllum commune: evolution in long-term culture.

The impact of ploidy on adaptation is a central issue in evolutionary biology. While many eukaryotic organisms exist as diploids, with two sets of gametic genomes residing in the same nucleus, most basidiomycete fungi exist as dikaryons in which the two genomes exist in separate nuclei that are physically paired and that divide in a coordinated manner during hyphal extension. To determine if haploid monokaryotic and dikaryotic mycelia adapt to novel environments under natural selection, we serially transferred replicate populations of each ploidy state on minimal medium for 18 months (approximately 13,000 generations). Dikaryotic mycelia responded to selection with increases in growth rate, while haploid monokaryotic mycelia did not. To determine if the haploid components of the dikaryon adapt reciprocally to one another's presence over time, we recovered the intact haploid components of dikaryotic mycelia at different time points (without meiosis) and mated them with nuclei of different evolutionary histories. We found evidence for coadaptation between nuclei in one dikaryotic line, in which a dominant deleterious mutation in one nucleus was followed by a compensatory mutation in the other nucleus; the mutant nuclei that evolved together had the best overall fitness. In other lines, nuclei had equal or higher fitness when paired with nuclei of other histories, indicating a heterozygote advantage. To determine if genetic exchange occurs between the two nuclei of a dikaryon, we developed a 24-locus genotyping system based on single nucleotide polymorphisms to monitor somatic exchange. We observed genetic exchange and recombination between the nuclei of several different dikaryons, resulting in genotypic variation in these mitotic cell lineages.

Base Sequence↗

[In vitro molecular evolution of humanized mouse ScFv specific against human cross-linked fibrin].

In early work at our laboratory, a mouse ScFv All specific against human cross-linked fibrin was picked out by using phage displaying technique, and it was useful in prepareing targeting thrombolytic agent. To reduce its immunogenicity, All was humanized according to the surface reshaping approach. In order to increase the affinity and specifity of humanized ScFv, the Mixed HCDR3 and LCDR3 library was constructed and selected, and five different humanized ScFv with better affinity or specifity than humanized ScFv were identified, but the affinity of these five humanized scFvs was still lower than parent mouse scFv All. In this study, in order to further increase the affinity of humanized scFv, the identified five different CDR3 mutant of humanized ScFvs were mutated and reassembled with the methods of error-prone PCR and DNA shuffling, then the shuffled ScFvs were cloned into plasmid pHB-1 HSCFV to construct ScFv library of 10(5). The ScFv library was displayed on the surface of phage and panned with antigen DD, then five ScFvs with better affinity or specifity were identified by using ScFv-alkaline phosphatase detection system which was established in our laboratory. Based on identified five better ScFv, the second round error-prone PCR and DNA shuffling program was taken out, a ScFv library of 10(6) was constructed and selected as before. And four humanized ScFv that had evidently better affinity than parent mouse ScFv was identified. This work laid a foundation for further research of targeting thrombolitic agent with low immunogenicity.

Animals↗

Enhanced synthesis of poly(3-hydroxybutyrate) in recombinant Escherichia coli by means of error-prone PCR mutagenesis, saturation mutagenesis, and in vitro recombination of the type II polyhydroxyalkanoate synthase gene.

Type II synthase (PhaC1(Ps)) for polyhydroxyalkanoate (PHA) from Pseudomonas sp. 61-3 was subjected to an in vitro evolution system including PCR-mediated mutagenesis in order to improve the function of PhaC1(Ps) in terms of its ability to produce poly(3-hydroxybutyrate) [P(3HB)] in recombinant Escherichia coli. Based on our established in vivo assay system, two positions (Ser325 and Gln481) where mutations provided remarkable increases in P(3HB) synthesis were identified. Saturation mutagenesis at these positions was carried out to explore whether there might be more beneficial sequences for P(3HB) synthesis than those identified in the point mutation library. As a result, five single mutants [S325C (T) and Q481M (K, R)] gave rise to highly enhanced P(3HB) synthesis. Drastically enhanced P(3HB) synthesis (up to 340- to 400-fold the amount of that of the wild type) was further achieved by generation of all five variants of the double mutants combining the codons for residues 325/481. It is feasible that the replacement of Ser (specific for type II synthase) by Thr (specific for type I synthase) at position 325 resulted in acquiring greater P(3HB) synthesis ability as exhibited by type I synthases. The other hot spot, 481, that positively contributes to enhanced P(3HB) synthesis is located adjacent to a His479, a residue that forms a putative catalytic diad that can be inferred by sequence alignment.

Acyltransferases↗

The adaptive dynamics of altruism in spatially heterogeneous populations.

We study the spatial adaptive dynamics of a continuous trait that measures individual investment in altruism. Our study is based on an ecological model of a spatially heterogeneous population from which we derive an appropriate measure of fitness. The analysis of this fitness measure uncovers three different selective processes controlling the evolution of altruism: the direct physiological cost, the indirect genetic benefits of cooperative interactions, and the indirect genetic costs of competition for space. In our model, habitat structure and a continuous life cycle makes the cost of competing for space with relatives negligible. Our study yields a classification of adaptive patterns of altruism according to the shape of the costs of altruism (with decelerating, linear, or accelerating dependence on the investment in altruism). The invasion of altruism occurs readily in species with accelerating costs, but large mutations are critical for altruism to evolve in selfish species with decelerating costs. Strict selfishness is maintained by natural selection only under very restricted conditions. In species with rapidly accelerating costs, adaptation leads to an evolutionarily stable rate of investment in altruism that decreases smoothly with the level of mobility. A rather different adaptive pattern emerges in species with slowly accelerating costs: high altruism evolves at low mobility, whereas a quasi-selfish state is promoted in more mobile species. The high adaptive level of altruism can be predicted solely from habitat connectedness and physiological parameters that characterize the pattern of cost. We also show that environmental changes that cause increased mobility in those highly altruistic species can beget selection-driven self-extinction, which may contribute to the rarity of social species.

Adaptation, Physiological↗

Fending off decay: a combinatorial approach in intact cells for identifying mRNA stability elements.

The strategy of systematic evolution, whereby nucleic acid sequences or conformers can be selected and amplified from a randomized population, has been exploited by many research groups for numerous purposes. It is, however, a technique largely performed in vitro, under nonphysiological conditions. We have now modified this in vitro approach to accomplish selection in growing cells. Here, we report that this new methodology has been used in vivo to select RNA elements that confer increased transcript stability. A randomized cassette was embedded in a 3'-untranslated region (UTR), downstream from the luciferase reporter open reading frame. A heterogeneous population of capped luciferase mRNA was then generated by in vitro transcription. Human liver Hep G2 cells were electroporated with this population of luciferase mRNA and total cytoplasmic RNA was isolated after varying lengths of incubation. Following RT-PCR, the 3' UTR was used to reconstruct a new population of luciferase templates, permitting subsequent cycles of in vitro transcription, electroporation, RNA isolation, and RT-PCR. Increasing the incubation time at each cycle before RNA isolation imposed selection for stable transcripts. The functional half-life of the luciferase mRNA population increased from 55 to 140 min after four cycles. Subsequent sequencing of the selected 3' UTRs revealed G-U rich elements in clones with extended chemical and functional half-lives.

3' Untranslated Regions↗

Evolutional design of a hyperactive cysteine- and methionine-free mutant of Escherichia coli dihydrofolate reductase.

We developed a strategy for finding out the adapted variants of enzymes, and we applied it to an enzyme, dihydrofolate reductase (DHFR), in terms of its catalytic activity so that we successfully obtained several hyperactive cysteine- and methionine-free variants of DHFR in which all five methionyl and two cysteinyl residues were replaced by other amino acid residues. Among them, a variant (M1A/M16N/M20L/M42Y/C85A/M92F/C152S), named as ANLYF, has an approximately seven times higher k(cat) value than wild type DHFR. Enzyme kinetics and crystal structures of the variant were investigated for elucidating the mechanism of the hyperactivity. Steady-state and transient binding kinetics of the variant indicated that the kinetic scheme of the catalytic cycle of ANLYF was essentially the same as that of wild type, showing that the hyperactivity was brought about by an increase of the dissociation rate constants of tetrahydrofolate from the enzyme-NADPH-tetrahydrofolate ternary complex. The crystal structure of the variant, solved and refined to an R factor of 0.205 at 1.9-angstroms resolution, indicated that an increased structural flexibility of the variant and an increased size of the N-(p-aminobenzoyl)-L-glutamate binding cleft induced the increase of the dissociation constant. This was consistent with a large compressibility (volume fluctuation) of the variant. A comparison of folding kinetics between wild type and the variant showed that the folding of these two enzymes was similar to each other, suggesting that the activity enhancement of the enzyme can be attained without drastic changes of the folding mechanism.

Binding Sites↗

Phage display selection of hairpin loop soyacystatin variants that mediate high affinity inhibition of a cysteine proteinase.

Two hairpin-loop domains in cystatin family proteinase inhibitors form an interface surface region that slots into the active site cleft of papain-like cysteine proteinases, and determine binding affinity. The slot region surface architecture of the soybean cysteine proteinase inhibitor (soyacystatin N, scN) was engineered using techniques of in vitro molecular evolution to define residues that facilitate interaction with the proteinase cleft and modulate inhibitor affinity and function. Combinatorial phage display libraries of scN variants that contain mutations in the essential motifs of the first (QVVAG) and second (EW) hairpin-loop regions were constructed. Approximately 1010-1011 phages expressing recombinant scN proteins were subjected to biopanning selection based on binding affinity to immobilized papain. The QVVAG motif in the first hairpin loop was invariant in all functional scN proteins. All selected variants (30) had W79 in the second hairpin-loop motif, but there was diversity for hydrophobic and basic amino acids in residue 78. Kinetic analysis of isolated scN variants identified a novel scN isoform scN(LW) with higher papain affinity than the wild-type molecule. The variant contained an E78L substitution and had a twofold lower Ki (2.1 pM) than parental scN, due to its increased association rate constant (2.6 +/- 0.09 x 107 M-1sec-1). These results define residues in the first and second hairpin-loop regions which are essential for optimal interaction between phytocystatins and papain, a prototypical cysteine proteinase. Furthermore, the isolated variants are a biochemical platform for further integration of mutations to optimize cystatin affinity for specific biological targets.

Base Sequence↗

Laboratory evolution of P450 BM-3 for mediated electron transfer.

Preparative synthesis with P450 monooxygenases is hampered in cell-free systems by the requirement for cofactors such as NAD(P)H as reduction equivalents. A validated medium-throughput screening system was designed for improving P450 monooxygenases by mediated electron transfer with zinc/cobalt(III)sepulchrate (Zn/Co(III)sep) as an alternative and cost-effective cofactor system. The monooxygenase P450 BM-3 F87A was used as a model system for developing the screening system in a 96-well format. A coefficient of variation of less than 10% was achieved under optimized screening conditions. The mediator evolution screen was validated by comparing the activity of P450 BM-3 to P450 BM-3 F87A and by screening a saturation mutagenesis library at amino acid position R47. For mediated electron transfer, two double mutants P450 BM-3(F87A R47F) and P450 BM-3 (F87A R47Y) were identified with a two-threefold increased catalytic efficiency (up to 32 microM(-1) min(-1) for P450 BM-3(F87A R47F) and 34 microM(-1) min(-1) for P450 BM-3 (F87A R47Y)) compared to P450 BM-3 F87A. The kinetic constants of the double mutants are, in contrast to those of P450 BM-3 F87A, dependent on Co(III)sep concentration in the presence of NADPH. kcat increases from 145 min(-1) (0.25 mM Co(III)sep) to 197 min(-1) (0.5 mM Co(III)sep), and Km decreases simultaneously from 7.0 microM to 3.7 microM, for P450 BM-3 (F87A R47F). For P450 BM-3 (F87A R47Y), kcat increases from 138 min(-1) (0.25 mM Co(III)sep) up to 187 min(-1) (0.5 mM Co(III)sep), and Km decreases from 8.2 microM to 4.2 microM. Due to lower Km values, the catalytic efficiencies were improved six times for P450 BM-3 (F87A R47F) and three times for P450 BM-3 (F87A R47Y), when comparing catalytic efficiencies of the mediated electron-transfer system to the natural reduction equivalent NADPH.

Amino Acid Substitution↗

Overexpression of calsequestrin in L6 myoblasts: formation of endoplasmic reticulum subdomains and their evolution into discrete vacuoles where aggregates of the protein are specifically accumulated.

Calsequestrin (CSQ), the major low-affinity Ca(2+)-binding glycoprotein of striated muscle fibers, is concentrated to yield aggregates that occupy the lumen of the terminal cisternae of the sarcoplasmic reticulum (SR). When infected or transfected into L6 myoblast, the protein is also concentrated, however, in dense vacuoles apparently separate from the endoplasmic reticulum (ER). CSQ-rich cells appear otherwise normal; in particular, neither other proteins involved in Ca2+ homeostasis nor ER chaperones are increased. The CSQ dense vacuoles are shown herein to be specialized ER subdomains as demonstrated by 1) the endoglycosidase H sensitivity of their CSQ and 2) two markers, calreticulin and calnexin (but not others, protein disulfide isomerase and BiP), intermixed with the vacuole content. Their formation is shown to start with the aggregation of CSQ at discrete sites of the ER lumen. When cells were transfected with both CSQ and calreticulin, only the first gave rise to vacuoles; the second remained diffusely distributed within the ER lumen. The possibility that CSQ aggregation is an artifact of overexpression appears unlikely because 1) within dense vacuoles CSQ molecules are not disulfide cross-linked, 2) their turnover is relatively slow (t = 12 h), and 3) segregated CSQ is bound to large amounts of Ca2+. Transfection of a tagged CSQ into cells already overexpressing the protein revealed the continuous import of the newly synthesized protein into preassembled vacuoles. The tendency to aggregation appears, therefore, as a property contributing to the segregation of CSQ within the ER lumen and to its accumulation within specialized subdomains. The study of L6 cells expressing CSQ-rich vacuoles might thus ultimately help to unravel mechanisms by which the complexity of the sarcoplasmic reticulum is established in muscle fibers.

Animals↗

New NTP analogs: the synthesis of 4'-thioUTP and 4'-thioCTP and their utility for SELEX.

The synthesis of the triphosphates of 4'-thiouridine and 4'-thiocytidine, 4'-thioUTP (7; thioUTP) and 4'-thioCTP (8; thioCTP), and their utility for SELEX (systematic evolution of ligands by exponential enrichment) is described. The new nucleoside triphosphate (NTP) analogs 7 and 8 were prepared from appropriately protected 4'-thiouridine and -cytidine derivatives using the one-pot method reported by J. Ludwig and F. Eckstein [(1989) J. Org. Chem., 54, 631-635]. Because SELEX requires both in vitro transcription and reverse transcription, we examined the ability of 7 and 8 for SELEX by focusing on the two steps. Incorporation of 7 and 8 by T7 RNA polymerase to give 4'-thioRNA (thioRNA) proceeded well and was superior to those of the two sets of frequently used modified NTP analogs for SELEX (2'-NH2dUTP and 2'-NH2dCTP; 2'-FdUTP and 2'-FdCTP), when an adequate leader sequence of DNA template was selected. We revealed that a leader sequence of about +15 of DNA template is important for the effective incorporation of modified NTP analogs by T7 RNA polymerase. In addition, reverse transcription of the resulting thioRNA into the complementary DNA in the presence of 2'-deoxynucleoside triphosphates (dNTPs) also proceeded smoothly and precisely. The stability of the thioRNA toward RNase A was 50 times greater than that of the corresponding natural RNA. With these successful results in hand, we attempted the selection of thioRNA aptamers to human alpha-thrombin using thioUTP and thioCTP, and found a thioRNA aptamer with high binding affinity (K(d) = 4.7 nM).

Base Sequence↗

Encoding of target direction and speed during visual instruction and arm tracking in dorsal premotor and primary motor cortical neurons.

The encoding of direction and speed in the discharge of dorsal premotor (PMd) and primary motor (MI) neurons was studied during two-dimensional visually-instructed pursuit arm movements in which eight directions and four constant speeds were independently manipulated. Each trial consisted of equal durations of visual observation of target movement without hand movement (cue) and visual pursuit-tracking of the target with the hand (track). A total of 240 neurons was recorded from PMd and MI in two Macaca mulatta monkeys. Two classes of regression analyses were used to relate neuronal firing during the cue and track periods to direction and speed. First, the average firing from each period was fitted to target direction or speed. Period-averaged firing significantly correlated with direction more frequently in the track than in the cue period. Conversely, correlations with speed (with or without direction) were more common in the cue than in the track period. Secondly, a binwise regression evaluated the temporal evolution of firing correlations with direction and speed. Supporting the period-based results, significant binwise correlations of the discharge with speed occurred preferentially during the cue period when there was no hand movement. Prior to movement, correlations of the firing with direction became significant and continued through the movement. Both analyses demonstrated a distinct tendency for neurons to be modulated by speed information early and by direction information later. This temporal parcellation reflects both the sequential demands of the task and constraints placed on the neural computations. The early representation of target speed is hypothesized to reflect the need to calculate a 'go signal' for the initiation of movement.

Action Potentials↗

In vitro evolution of the hammerhead ribozyme to a purine-specific ribozyme using mutagenic PCR with two nucleotide analogues.

The conventional hammerhead ribozyme cleaves RNA 3' to nucleotide triplets with the general formula NUH, where N is any nucleotide, U is uridine and H is any nucleotide except guanosine. In order to isolate hammerhead ribozyme sequences capable of cleaving 3' to the GUG triplet, we performed a mutagenic selection protocol starting with the conventional sequence of an NUH-cleaving ribozyme. The 22 nucleotides in the core and the stem-loop II region were subjected to mutagenic PCR using the two nucleotide analogues 6-(2-deoxy-beta-d-ribofuranosyl)-3,4-dihydro-8H-pyrimido-[4,5-C)][1, 2] oxazin-7-one and of 8-oxo-2'-deoxyguanosine. After five repetitions of the selection cycle, several clones showed cleavage activity. One sequence, having one deletion, showed at least a 90 times higher in trans cleavage rate than the starting ribozyme. It cleaved 3' to GUG and GUA. The sequence of this ribozyme is essentially identical with that obtained previously by selection for AUG cleavage starting with a randomised core and stem-loop II region. This identical result of two independent selection procedures supports the notion that sequences for NUR cleavage, where R is a purine nucleotide, are not compatible with the classical hammerhead structure, and that the sequence space for this cleavage specificity is very limited. The cleavage of NUR triplets is not restricted to the sequence of the substrate that was used for selection but is sequence-independent for in trans cleavage, although the sequence context influences the value for the cleavage rate somewhat. Analysis of cleavage activities indicates the importance of A at position L2.5 in loop II.

8-Hydroxy-2'-Deoxyguanosine↗

A cold-adapted protease engineered by experimental evolution system.

A new cold-adapted protease subtilisin BPN' mutant, termed m-51, was successfully isolated by use of an evolutionary program consisting of two-step in vitro random mutagenesis, which we developed for the screening of mutant subtilisins with increased activity at low temperature. The m-51 mutant showed 70% higher catalytic efficiency, expressed by the k(cat)/K(m) value, than the wild-type at 10 degrees C against N-succinyl-L-Ala-L-Ala-L-Pro-L-Phe-p-nitroanilide as a synthetic substrate. This cold-adaptation was achieved mainly by the increase in the k(cat) value in a temperature-dependent manner. Genetic analysis revealed that m-51 had three mutations, Ala-->Thr at position -31 (A-31T) in the prodomain, Ala-->Val at position 88 (A88V), and Ala-->Thr at position 98 (A98T). From kinetic parameters of the purified mutant enzymes, it was found that the A98T mutation led to 30% activity increase, which was enhanced up to 70% by the accompanying neutral mutation A88V. The A-31T mutation severely constrained the autoprocessing-mediated maturation of the pro-subtilisin in the Escherichia coli expression system, thus probably causing an activity-non-detectable mutation in the first step of mutagenesis. No distinct change was observed in the thermal stability of any mutant or in the substrate specificity for m-51. In the molecular models of the two single mutants (A88V and A98T), relatively large displacements of alpha carbon atoms were found around the mutation points. In the model of the double mutant (A88V/A98T), on the other hand, the structural changes around the mutation point counterbalanced each other, and thus no crucial displacements occurred. This mutual effect may be related to the enhanced activity of the double mutant.

Base Sequence↗

In vitro selection of DNA aptamers to anthrax spores with electrochemiluminescence detection.

Systematic evolution of ligands by exponential enrichment (SELEX) was used to select and PCR amplify DNA sequences (aptamers) capable of binding to and detecting nonpathogenic Sterne strain Bacillus anthracis spores. A simplified affinity separation approach was employed, in which autoclaved anthrax spores were used as the separation matrix. An aptamer-magnetic bead-electrochemiluminescence (AM-ECL) sandwich assay scheme was devised for detecting anthrax spores. Using a low SELEX DNA to spore ratio (154 ng DNA/10(6) spores), at least three distinct populations of single-stranded DNA aptamers, having varied affinities for anthrax spores, were noted by the AM-ECL assay. Results reflect detection of spore components with a dynamic range equivalent to < 10- > 6 x 10(6) anthrax spores. In the low DNA to spore ratio experiments, aptamers could be liberated from spore pellets by heating at 96 degrees C for 5 min after each round of SELEX. When a much higher DNA to spore ratio (10,256 ng DNA/10(6) spores) was used for SELEX development, a higher affinity set of aptamers was selected that could not be heat-eluted even at 99 degrees C for 5 min following round four of SELEX. However, high affinity spore surface bound aptamers were detectable via their 5'-biotinylated tails using labeled avidin and could be eluted in deionized water. Aptamers have potential for use as inexpensive, in vitro-generated receptors for biosensors in biological warfare detection and other areas.

Bacillus anthracis↗

In vitro evolution of the binding specificity of neocarzinostatin, an enediyne-binding chromoprotein.

Neocarzinostatin is the most studied member of the enediyne-chromoprotein family, and is clinically used as an antitumoral agent. Neocarzinostatin could be a promising drug delivery vehicle if new binding specificities could be conferred to its protein scaffold. We used in vitro evolution methods to demonstrate that this approach is feasible. We created large libraries containing between 1.7 x 10(8) and 1.4 x 10(9) independent clones, where up to 13 side chains pointing toward the binding crevice were randomly substituted. We then used phage display to select variants that bind to a model ligand (testosterone) which is unrelated to the natural ligand of neocarzinostatin. Several different binders were selected from each library. The corresponding proteins were expressed in Escherichia coli and their affinities and specificities were characterized in detail. K(D) values of about 20 nM were obtained for streptavidin-bound testosterone. The K(D) of selected proteins for free soluble testosterone are between 7 and 55 microM and therefore higher than the K(D) for streptavidin-bound testosterone. The spacer and streptavidin used during selection contributed to the high affinity of the selected binders for the target. Binding studies of 15 different steroids related to testosterone allowed us to determine that C3, 4, 5, 6, and 7 on cycles A and B and the conjugated 3 oxo group of the steroid molecule were essential for molecular recognition. Other testosterone analogues substituted on C1, 2, 9, 11, 15, and 17 were not discriminated from testosterone. These results demonstrate that the binding specificity of this protein family can be extended to compounds that are completely unrelated to the natural enediyne chromophore family. This type of highly expressed, stable proteins with tailored binding properties have a wide potential range of applications.

Amino Acid Sequence↗

A trans acting ribozyme that phosphorylates exogenous RNA.

The structural complexity required for substrate recognition within an active site constrains the evolution of novel catalytic functions. To evaluate those constraints within populations of incipient ribozymes, we performed a selection for kinase ribozymes under conditions that allowed competition for phosphorylation at nine candidate sites. Two candidate sites are the hydroxyl groups on a "quasi-diffusible" chloramphenicol (Cam) moiety tethered to the evolving library through an inert, flexible linker. A subtractive step was included to allow only seven ribose 2' hydroxyls to compete with the two Cam hydroxyls for phosphorylation. After the library was incubated with gamma-thio-ATP (ATPgammaS), active species were recovered from a polyacrylamide gel containing [(N-acryloylamino)phenyl] mercury (APM) and amplified for further cycles of selection. Activity assays on selected isolates and truncated derivatives identified the essential secondary structure of the dominant RNA motif. Phosphorylation was independent of the Cam moiety, indicating ribose 2' phosphorylation. The dominant motif was separated into catalytic "ribozyme" and "substrate" strands. Partial alkaline digestion of the substrate strand before and after phosphorylation identified the precise modification site as the first purine (R) within the required sequence 5'-RAAAANCG-3'. The reaction shows approximately 10-fold preference for ATPgammaS over ATP and is independent of pH over a wide range (5.5-8.9), consistent with a dissociative reaction mechanism that is rate-limited by formation of a metaphosphate transition state. Divalent metal ions are required, with a slight preference of Mn(2+) > Mg(2+) > Ca(2+). Lack of reactivity in [Co(NH(3))(6)](3+) indicates a requirement for inner sphere contact with the metal ion, either for structural stabilization, catalysis, or both.

Adenosine Triphosphate↗

A stable disulfide-free gene-3-protein of phage fd generated by in vitro evolution.

Disulfide bonds provide major contributions to the conformational stability of proteins, and their cleavage often leads to unfolding. The gene-3-protein of the filamentous phage fd contains two disulfides in its N1 domain and one in its N2 domain, and these three disulfide bonds are essential for the stability of this protein. Here, we employed in vitro evolution to generate a disulfide-free variant of the N1-N2 protein with a high conformational stability. The gene-3-protein is essential for the phage infectivity, and we exploited this requirement for a proteolytic selection of stabilized protein variants from phage libraries. First, optimal replacements for individual disulfide bonds were identified in libraries, in which the corresponding cysteine codons were randomized. Then stabilizing amino acid replacements at non-cysteine positions were selected from libraries that were created by error-prone PCR. This stepwise procedure led to variants of N1-N2 that are devoid of all three disulfide bonds but stable and functional. The best variant without disulfide bonds showed a much higher conformational stability than the disulfide-containing wild-type form of the gene-3-protein. Despite the loss of all three disulfide bonds, the midpoints of the thermal transitions were increased from 48.5 degrees C to 67.0 degrees C for the N2 domain and from 60.0 degrees C to 78.7 degrees C for the N1 domain. The major loss in conformational stability caused by the removal of the disulfides was thus over-compensated by strongly improved non-covalent interactions. The stabilized variants were less infectious than the wild-type protein, probably because the domain mobility was reduced. Only a small fraction of the sequence space could be accessed by using libraries created by error-prone PCR, but still many strongly stabilized variants could be identified. This is encouraging and indicates that proteins can be stabilized by mutations in many different ways.

Bacteriophage M13↗

Experimental rugged fitness landscape in protein sequence space.

The fitness landscape in sequence space determines the process of biomolecular evolution. To plot the fitness landscape of protein function, we carried out in vitro molecular evolution beginning with a defective fd phage carrying a random polypeptide of 139 amino acids in place of the g3p minor coat protein D2 domain, which is essential for phage infection. After 20 cycles of random substitution at sites 12-130 of the initial random polypeptide and selection for infectivity, the selected phage showed a 1.7x10(4)-fold increase in infectivity, defined as the number of infected cells per ml of phage suspension. Fitness was defined as the logarithm of infectivity, and we analyzed (1) the dependence of stationary fitness on library size, which increased gradually, and (2) the time course of changes in fitness in transitional phases, based on an original theory regarding the evolutionary dynamics in Kauffman's n-k fitness landscape model. In the landscape model, single mutations at single sites among n sites affect the contribution of k other sites to fitness. Based on the results of these analyses, k was estimated to be 18-24. According to the estimated parameters, the landscape was plotted as a smooth surface up to a relative fitness of 0.4 of the global peak, whereas the landscape had a highly rugged surface with many local peaks above this relative fitness value. Based on the landscapes of these two different surfaces, it appears possible for adaptive walks with only random substitutions to climb with relative ease up to the middle region of the fitness landscape from any primordial or random sequence, whereas an enormous range of sequence diversity is required to climb further up the rugged surface above the middle region.

Amino Acid Sequence↗