Directed evolution of a (βα)8-barrel enzyme to catalyze related reactions in two different metabolic pathways

Enzymes participating in different metabolic pathways often have similar catalytic mechanisms and structures, suggesting their evolution from a common ancestral precursor enzyme. We sought to create a precursor-like enzyme for N′-[(5′-phosphoribosyl)formimino]-5-aminoimidazole-4-carboxamide ribonucleotide (ProFAR) isomerase (HisA; EC 5.3.1.16) and phosphoribosylanthranilate (PRA) isomerase (TrpF; EC 5.3.1.24), which catalyze similar reactions in the biosynthesis of the amino acids histidine and tryptophan and have a similar (βα)8-barrel structure. Using random mutagenesis and selection, we generated several HisA variants that catalyze the TrpF reaction both in vivo and in vitro, and one of these variants retained significant HisA activity. A more detailed analysis revealed that a single amino acid exchange could establish TrpF activity on the HisA scaffold. These findings suggest that HisA and TrpF may have evolved from an ancestral enzyme of broader substrate specificity and underscore that (βα)8-barrel enzymes are very suitable for the design of new catalytic activities.

Directed evolution of polymerase function by compartmentalized self-replication

We describe compartmentalized self-replication (CSR), a strategy for the directed evolution of enzymes, especially polymerases. CSR is based on a simple feedback loop consisting of a polymerase that replicates only its own encoding gene. Compartmentalization serves to isolate individual self-replication reactions from each other. In such a system, adaptive gains directly (and proportionally) translate into genetic amplification of the encoding gene. CSR has applications in the evolution of polymerases with novel and useful properties. By using three cycles of CSR, we obtained variants of Taq DNA polymerase with 11-fold higher thermostability than the wild-type enzyme or with a >130-fold increased resistance to the potent inhibitor heparin. Insertion of an extra stage into the CSR cycle before the polymerase reaction allows its application to enzymes other than polymerases. We show that nucleoside diphosphate kinase and Taq polymerase can form such a cooperative CSR cycle based on reciprocal catalysis, whereby nucleoside diphosphate kinase produces the substrates required for the replication of its own gene. We also find that in CSR the polymerase genes themselves evolve toward more efficient replication. Thus, polymerase genes and their encoded polypeptides cooperate to maximize postselection copy number. CSR should prove useful for the directed evolution of enzymes, particularly DNA or RNA polymerases, as well as for the design and study of in vitro self-replicating systems mimicking prebiotic evolution and viral replication.

Computer simulations of enzyme catalysis: Finding out what has been optimized by evolution

The origin of the catalytic power of enzymes is discussed, paying attention to evolutionary constraints. It is pointed out that enzyme catalysis reflects energy contributions that cannot be determined uniquely by current experimental approaches without augmenting the analysis by computer simulation studies. The use of energy considerations and computer simulations allows one to exclude many of the popular proposals for the way enzymes work. It appears that the standard approaches used by organic chemists to catalyze reactions in solutions are not used by enzymes. This point is illustrated by considering the desolvation hypothesis and showing that it cannot account for a large increase in kcat relative to the corresponding kcage for the reference reaction in a solvent cage. The problems associated with other frequently invoked mechanisms also are outlined. Furthermore, it is pointed out that mutation studies are inconsistent with ground state destabilization mechanisms. After considering factors that were not optimized by evolution, we review computer simulation studies that reproduced the overall catalytic effect of different enzymes. These studies pointed toward electrostatic effects as the most important catalytic contributions. The nature of this electrostatic stabilization mechanism is far from being obvious because the electrostatic interaction between the reacting system and the surrounding area is similar in enzymes and in solution. However, the difference is that enzymes have a preorganized dipolar environment that does not have to pay the reorganization energy for stabilizing the relevant transition states. Apparently, the catalytic power of enzymes is stored in their folding energy in the form of the preorganized polar environment.

Stepwise Photoinhibition of Photosystem II1 : Studies with Synechocystis Species PCC 6803 Mutants with a Modified D-E Loop of the Reaction Center Polypeptide D1

Several mutant strains of Synechocystis sp. PCC 6803 with large deletions in the D-E loop of the photosystem II (PSII) reaction center polypeptide D1 were subjected to high light to investigate the role of this hydrophilic loop in the photoinhibition cascade of PSII. The tolerance of PSII to photoinhibition in the autotrophic mutant ΔR225-F239 (PD), when oxygen evolution was monitored with 2,6-dichloro-p-benzoquinone and the equal susceptibility compared with control when monitored with bicarbonate, suggested an inactivation of the QB-binding niche as the first event in the photoinhibition cascade in vivo. This step in PD was largely reversible at low light without the need for protein synthesis. Only the next event, inactivation of QA reduction, was irreversible and gave a signal for D1 polypeptide degradation. The heterotrophic deletion mutants ΔG240-V249 and ΔR225-V249 had severely modified QB pockets, yet exhibited high rates of 2,6-dichloro-p-benzoquinone-mediated oxygen evolution and less tolerance to photoinhibition than PD. Moreover, the protein-synthesis-dependent recovery of PSII from photoinhibition was impaired in the ΔG240-V249 and ΔR225-V249 mutants because of the effects of the mutations on the expression of the psbA-2 gene. No specific sequences in the D-E loop were found to be essential for high rates of D1 polypeptide degradation.

Duplication and functional divergence in the chalcone synthase gene family of Asteraceae: evolution with substrate change and catalytic simplification.

Plant-specific polyketide synthase genes constitute a gene superfamily, including universal chalcone synthase [CHS; malonyl-CoA:4-coumaroyl-CoA malonyltransferase (cyclizing) (EC 2.3.1.74)] genes, sporadically distributed stilbene synthase (SS) genes, and atypical, as-yet-uncharacterized CHS-like genes. We have recently isolated from Gerbera hybrida (Asteraceae) an unusual CHS-like gene, GCHS2, which codes for an enzyme with structural and enzymatic properties as well as ontogenetic distribution distinct from both CHS and SS. Here, we show that the GCHS2-like function is encoded in the Gerbera genome by a family of at least three transcriptionally active genes. Conservation within the GCHS2 family was exploited with selective PCR to study the occurrence of GCHS2-like genes in other Asteraceae. Parsimony analysis of the amplified sequences together with CHS-like genes isolated from other taxa of angiosperm subclass Asteridae suggests that GCHS2 has evolved from CHS via a gene duplication event that occurred before the diversification of the Asteraceae. Enzyme activity analysis of proteins produced in vitro indicates that the GCHS2 reaction is a non-SS variant of the CHS reaction, with both different substrate specificity (to benzoyl-CoA) and a truncated catalytic profile. Together with the recent results of Durbin et al. [Durbin, M. L., Learn, G. H., Jr., Huttley, G. A. & Clegg, M. T. (1995) Proc. Natl. Acad. Sci. USA 92, 3338-3342], our study confirms a gene duplication-based model that explains how various related functions have arisen from CHS during plant evolution.

Crystal structure of D-amino acid oxidase: a case of active site mirror-image convergent evolution with flavocytochrome b2.

D-amino acid oxidase is the prototype of the FAD-dependent oxidases. It catalyses the oxidation of D-amino acids to the corresponding alpha-ketoacids. The reducing equivalents are transferred to molecular oxygen with production of hydrogen peroxide. We have solved the crystal structure of the complex of D-amino acid oxidase with benzoate, a competitive inhibitor of the substrate, by single isomorphous replacement and eightfold averaging. Each monomer is formed by two domains with an overall topology similar to that of p-hydroxybenzoate hydroxylase. The benzoate molecule lays parallel to the flavin ring and is held in position by a salt bridge with Arg-283. Analysis of the active site shows that no side chains are properly positioned to act as the postulated base required for the catalytic carboanion mechanism. On the contrary, the benzoate binding mode suggests a direct transfer of the substrate alpha-hydrogen to the flavin during the enzyme reductive half-reaction.The active site Of D-amino acid oxidase exhibits a striking similarity with that of flavocytochrome b2, a structurally unrelated FMN-dependent flavoenzyme. The active site groups (if these two enzymes are in fact superimposable once the mirror-image of the flavocytochrome b2 active site is generated with respect to the flavin plane. Therefore, the catalytic sites of D-amino acid oxidase and flavocytochrome b2 appear to have converged to a highly similar but enantiomeric architecture in order to catalvze similar reactions (oxidation of alpha-amino acids or alpha-hydroxy acids), although with opposite stereochemistry.

Rates of decomposition of ribose and other sugars: implications for chemical evolution.

The existence of the RNA world, in which RNA acted as a catalyst as well as an informational macromolecule, assumes a large prebiotic source of ribose or the existence of pre-RNA molecules with backbones different from ribose-phosphate. The generally accepted prebiotic synthesis of ribose, the formose reaction, yields numerous sugars without any selectivity. Even if there were a selective synthesis of ribose, there is still the problem of stability. Sugars are known to be unstable in strong acid or base, but there are few data for neutral solutions. Therefore, we have measured the rate of decomposition of ribose between pH 4 and pH 8 from 40 degrees C to 120 degrees C. The ribose half-lives are very short (73 min at pH 7.0 and 100 degrees C and 44 years at pH 7.0 and 0 degrees C). The other aldopentoses and aldohexoses have half-lives within an order of magnitude of these values, as do 2-deoxyribose, ribose 5-phosphate, and ribose 2,4-bisphosphate. These results suggest that the backbone of the first genetic material could not have contained ribose or other sugars because of their instability.

Kinetic Analysis and Pollutant Evolution During Thermal Degradation of Printed Circuit Boards Considering the Effect of Metals

Paper submitted to the 7th International Symposium on Feedstock Recycling of Polymeric Materials (7th ISFR 2013), New Delhi, India, 23-26 October 2013.

The chances of death, and other studies in evolution.

"Errata" slip tipped in at v. 1, p. [1]; Addenda slip tipped in at v. 2, p. [413]

A Lanczos-powered implementation of the Faber polynomial quantum time propagator for reaction probabilities

Complementing our recent work on subspace wavepacket propagation [Chem. Phys. Lett. 336 (2001) 149], we introduce a Lanczos-based implementation of the Faber polynomial quantum long-time propagator. The original version [J. Chem. Phys. 101 (1994) 10493] implicitly handles non-Hermitian Hamiltonians, that is, those perturbed by imaginary absorbing potentials to handle unwanted reflection effects. However, like many wavepacket propagation schemes, it encounters a bottleneck associated with dense matrix-vector multiplications. Our implementation seeks to reduce the quantity of such costly operations without sacrificing numerical accuracy. For some benchmark scattering problems, our approach compares favourably with the original. (C) 2004 Elsevier B.V. All rights reserved.

Bioactive glass sol-gel foam scaffolds:evolution of nanoporosity during processing and in situ monitoring of apatite layer formation using small- and wide-angle X-ray scattering

Recent work has highlighted the potential of sol-gel-derived calcium silicate glasses for the regeneration or replacement of damaged bone tissue. The work presented herein provides new insight into the processing of bioactive calcia-silica sol-gel foams, and the reaction mechanisms associated with them when immersed in vitro in a simulated body fluid (SBF). Small-angle X-ray scattering and wide-angle X-ray scattering (diffraction) have been used to study the stabilization of these foams via heat treatment, with analogous in situ time-resolved data being gathered for a foam immersed in SBF. During thermal processing, pore sizes have been identified in the range of 16.5-62.0 nm and are only present once foams have been heated to 400 degrees C and above. Calcium nitrate crystallites were present until foams were heated to 600 degrees C; the crystallite size varied from 75 to 145 nm and increased in size with heat treatment up to 300 degrees C, then decreased in size down to 95 rim at 400 degrees C. The in situ time-resolved data show that the average pore diameter decreases as a function of immersion time in SBF, as calcium phosphates grow on the glass surfaces. Over the same time, Bragg peaks indicative of tricalcium phosphate were evident after only 1-h immersion time, and later, hydroxycarbonate apatite was also seen. The hydroxycarbonate apatite appears to have preferred orientation in the (h,k,0) direction.

Persistent localized states for a chaotically mixed bistable reaction

We describe the evolution of a bistable chemical reaction in a closed two-dimensional chaotic laminar flow, from a localized initial disturbance. When the fluid mixing is sufficiently slow, the disturbance may spread and eventually occupy the entire fluid domain. By contrast, rapid mixing tends to dilute the initial state and so extinguish the disturbance. Such a dichotomy is well known. However, we report here a hitherto apparently unremarked intermediate case, a persistent highly localized disturbance. Such a localized state arises when the Damkoehler number is great enough to sustain a "hot spot," but not so great as to lead to global spread. We show that such a disturbance is located in the neighborhood of an unstable periodic orbit of the flow, and we describe some limited aspects of its behavior using a reduced, lamellar model. Copyright American Physical Society (APS) 2006.

A bistable reaction-diffusion system in a stretching flow

We examine the evolution of a bistable reaction in a one-dimensional stretching flow, as a model for chaotic advection. We derive two reduced systems of ordinary differential equations (ODEs) for the dynamics of the governing advection-reaction-diffusion partial differential equations (PDE), for pulse-like and for plateau-like solutions, based on a non-perturbative approach. This reduction allows us to study the dynamics in two cases: first, close to a saddle-node bifurcation at which a pair of nontrivial steady states are born as the dimensionless reaction rate (Damkoehler number) is increased, and, second, for large Damkoehler number, far away from the bifurcation. The main aim is to investigate the initial-value problem and to determine when an initial condition subject to chaotic stirring will decay to zero and when it will give rise to a nonzero final state. Comparisons with full PDE simulations show that the reduced pulse model accurately predicts the threshold amplitude for a pulse initial condition to give rise to a nontrivial final steady state, and that the reduced plateau model gives an accurate picture of the dynamics of the system at large Damkoehler number. Published in Physica D (2006)

Chaotic mixing of a competitive-consecutive reaction

The evolution of a competitive-consecutive chemical reaction is computed numerically in a two-dimensional chaotic fluid flow with initially segregated reactants. Results from numerical simulations are used to evaluate a variety of reduced models commonly adopted to model the full advection-reaction-diffusion problem. Particular emphasis is placed upon fast reactions, where the yield varies most significantly with Peclet number (the ratio of diffusive to advective time scales). When effects of the fluid mechanical mixing are strongest, we find that the yield of the reaction is underestimated by a one-dimensional lamellar model that ignores the effects of fluid mixing, but overestimated by two other lamellar models that include fluid mixing.

Evolution of genome sequencing techniques

The quality and the speed for genome sequencing has advanced at the same time that technology boundaries are stretched. This advancement has been divided so far in three generations. The first-generation methods enabled sequencing of clonal DNA populations. The second-generation massively increased throughput by parallelizing many reactions while the third-generation methods allow direct sequencing of single DNA molecules. The first techniques to sequence DNA were not developed until the mid-1970s, when two distinct sequencing methods were developed almost simultaneously, one by Alan Maxam and Walter Gilbert, and the other one by Frederick Sanger. The first one is a chemical method to cleave DNA at specific points and the second one uses ddNTPs, which synthesizes a copy from the DNA chain template. Nevertheless, both methods generate fragments of varying lengths that are further electrophoresed. Moreover, it is important to say that until the 1990s, the sequencing of DNA was relatively expensive and it was seen as a long process. Besides, using radiolabeled nucleotides also compounded the problem through safety concerns and prevented the automation. Some advancements within the first generation include the replacement of radioactive labels by fluorescent labeled ddNTPs and cycle sequencing with thermostable DNA polymerase, which allows automation and signal amplification, making the process cheaper, safer and faster. Another method is Pyrosequencing, which is based on the “sequencing by synthesis” principle. It differs from Sanger sequencing, in that it relies on the detection of pyrophosphate release on nucleotide incorporation. By the end of the last millennia, parallelization of this method started the Next Generation Sequencing (NGS) with 454 as the first of many methods that can process multiple samples, calling it the 2º generation sequencing. Here electrophoresis was completely eliminated. One of the methods that is sometimes used is SOLiD, based on sequencing by ligation of fluorescently dye-labeled di-base probes which competes to ligate to the sequencing primer. Specificity of the di-base probe is achieved by interrogating every 1st and 2nd base in each ligation reaction. The widely used Solexa/Illumina method uses modified dNTPs containing so called “reversible terminators” which blocks further polymerization. The terminator also contains a fluorescent label, which can be detected by a camera. Now, the previous step towards the third generation was in charge of Ion Torrent, who developed a technique that is based in a method of “sequencing-by-synthesis”. Its main feature is the detection of hydrogen ions that are released during base incorporation. Likewise, the third generation takes into account nanotechnology advancements for the processing of unique DNA molecules to a real time synthesis sequencing system like PacBio; and finally, the NANOPORE, projected since 1995, also uses Nano-sensors forming channels obtained from bacteria that conducts the sample to a sensor that allows the detection of each nucleotide residue in the DNA strand. The advancements in terms of technology that we have nowadays have been so quick, that it makes wonder: ¿How do we imagine the next generation?