Studies on evolutionary and selective properties of hypercycles using a Monte Carlo method

Summary The most relevant properties of hypercycles were previously studied mainly from a theoretical point of view. We have developed a Monte Carlo method simulating hypercyclic organization to obtain information about the dynamics of this prebiotic organization. Nucleation, growth, and selective properties have been tested and the results obtained are in good agreement with those of the theoretical predictions. The influence of hypercyclic organization of the error threshold has also been studied. As a consequence of the emergence of a hypercycle, the value of this threshold decreases. The amount of this decrease depends on the population size. Moreover, for some interval of quality factor values, either the hypercycle organization or an error catastrophe can be produced, depending on the initial conditions. The influence of these phenomena on both the dynamic behavior and evolutionary advantages of the hypercycle, as well as their decisive roles on genome size, are discussed.Presented at the FEBS Symposium on Genome Organization and Evolution, held in Crete, Greece, September 1–5, 1986     

"Living" under the challenge of information decay: the stochastic corrector model vs. hypercycles

The combined problem of having a large genome size when the accuracy of replication was a limiting factor is probably the most difficult transition to explain at the late stages of RNA world. One solution has been to suggest the existence of a cyclically coupled system of autocatalytic and cross-catalytic molecular mutualists, where each member helps the following member and receives help from the preceding one (i.e., a "hypercycle"). However, such a system is evolutionarily unstable when mutations are taken into account because it lacks individuality. In time, the cooperating networks of genes should have been encapsulated in a cell-like structure. But once the cell was invented, it closely aligned genes' common interests and helped to reduce gene selfishness, so there was no need for hypercycles. A simple package of competing genes, described by the "stochastic corrector model" (SCM), could have provided the solution. Until now, there is no clear demonstration that the proposed mechanisms (compartmentalized hypercycles and the stochastic corrector model) do in fact solve the error threshold problem. Here, we present a Monte Carlo model to test the viability of protocell populations that enclose a hypercyclic (HPC) or a non-hypercyclic (SCM) system when faced with realistic mutation rates before the evolution of efficient enzymic machinery for replication. The numerical results indicate that both systems are efficient information integrators and are able to overcome the danger of information decay in the absence of accurate replication. However, a population of SCM protocells can tolerate higher deleterious mutation rates and reaches an equilibrium mutational load lower than that in a population of HPC protocells.     

Growth of a hypercycle and comparison with conventional autocatalysis

The hypercycle is a chemical model for reproduction which has been advocated for early stages of biological evolution. Its advantage is thought to lie in the high rate of growth conferred by hyperbolic kinetics. Earlier work has also indicated a saturation phase at large constituent concentrations. The present paper shows that both the saturation phase and the phase of hyperbolic growth have been introduced into the kinetics by making some of the reactions reversible. Reversibility is not essential to the operation of the hypercycle and the system with irreversible reactions grows faster. However, even the fastest hypercycle does not grow as fast as the simple auto-catalyst obtained by removing that reaction which is characteristic of the hypercycle. Also, both because the hypercycle is not a coherent system and because its growth requires reactions among separated constituents in the medium, it is more subject to decay than a simple autocatalytic particle. With greater complexity, slower growth, and more severe decay, the hypercycle is not a satisfactory alternative to conventional models of reproduction.     

Group selection of early replicators and the origin of life

A major problem of the origin of life has been that of information integration. As Eigen (1971) has shown, a mutant distribution of RNAs replicating without the aid of a replicase cannot integrate sufficient information for the functioning of a higher-level unit utilizing several types of encoded enzymes. He proposed the hypercycle model to bridge this gap in prebiology. It can be shown by a nonlinear game model, incorporating mutation of a hypercycle, that the selection properties of hypercycles make them inefficient information integrators as they cannot compete favourably with all kinds of less efficient information carriers or mutationally coupled hypercycles. The stochastic corrector model is presented as an alternative resolution of Eigen's paradox. It assumes that replicative templates are competing within replicative compartments, whose selective values depend on the internal template composition via a catalytic acid in replication and "metabolism". The dynamics of template replication are analyzed by numerical simulation of master equations. Due to the stochasticity in replication and compartment fission the best compartment types recur. An Eigen equation at the compartment level is set up and calculated. Even selfish template mutants cannot destroy the system though they make it less efficient. The genetic information of templates is evaluated at both levels, and the higher (compartment) level successfully constrains the lower (template) one. Compartmentation together with stochastic effects is sufficient to integrate information dispersed in competitive replicators. Compartment selection is considered to be group selection of replicators. Implications for the origin of life are discussed.     

Alteration of the quaternary structure of human UDP-glucose dehydrogenase by a double mutation

There are conflicting views for the polymerization process of human UDP-glucose dehydrogenase (UGDH) and no clear evidence has been reported yet. Based on crystal coordinates for Streptococcus pyogenes UGDH, we made double mutant A222Q/S233G. The double mutagenesis had no effects on expression, stability, and secondary structure. Interestingly, A222Q/S233G was a dimeric form and showed an UGDH activity, although it showed increased Km values for substrates. These results suggest that Ala222 and Ser233 play an important role in maintaining the hexameric structure and the reduced binding affinities for substrates are attributable to its altered subunit communication although quaternary structure may not be critical for catalysis.     

Phase diagrams for DNA crystallization systems

Phase diagrams for several oligonucleotide duplex-spermine systems have been constructed. These diagrams characterize the duplex and spermine concentrations ranges in which crystalline precipitates are formed. All of them are wedge-like form. The slope of the upper branch of the diagram is determined by the oligonucleotide length. The position of the lower branch depends on both the nucleotide sequence and its length. The position of the lower branch depends on both the nucleotide sequence and its length. It has been shown that the addition to the system of MgCl2 and NaCl salts and MPD results in specific changes in the diagrams. A model for oligonucleotide duplex-spermine system has been suggested which explains the main characteristic features of the obtained phase diagrams. The experimental phase diagrams for the (pGpT)n (pApC)n-spermine system (n = 2,3,4) have been analyzed ion terms of this model and the values of the binding constants of spermine and Mg2+ ions binding to duplexes have been determined. It permitted to identify the complexes that precipitated in different regions of the phase diagrams under various conditions. The diagram obtained in the presence of a cobalt hexammine counterion is also considered. It has been shown that this phase diagram, in general, is similar to those obtained for the oligonucleotide duplex-spermine system.     

Stability in simple grazing models: effects of explicit functions

A class of simple models of a grazed pasture is based on the assumption that the dynamics of green biomass can be expressed as a balance between rates of growth and consumption, both of which are functions only of the total amount of green biomass present. Graphical analysis, using only the general shape of these functions, has shown that the pasture may be either continuously or discontinuously stable under grazing. Further results require explicit specification of these functions. Four growth functions and four consumption functions are considered here and their mathematical and biological properties discussed. The stability properties of the 16 models resulting from combinations of these functions are compared by plotting the isoclines of zero pasture growth, and by obtaining expressions in terms of the parameters for some critical points and quantities. Both methods show considerable similarity between the models in their qualitative behaviour and in the role of certain parameters in determining critical values. There are substantial differences in quantitative predictions of critical values and particularly of the boundary between continuous and discontinuous stability. In the absence of an ungrazable residual, some models are discontinuously stable in the whole parameter space, others only in a limited region of it-albeit the region in which most real systems may be.     

The correlated actions of vitamins C,K, and Q

Three vitamins (C, K, and Q), two of which are unequivocally established and the existence of the third supported by both experimental and clinical evidence, are needed to prevent hemorrhagic states and therefore can be designated the hemostatic vitamins. The coordinated actions of these vitamins can be epitomized by a diagram. Vitamin K is responsible for the synthesis of four basic clotting factors that function in two distinct pathways for the activation of prothrombin to thrombin. Vitamin Q also functions in these pathways. In the intrinsic, it supplies by means of platelets the Q factor that with factors VIII and IX generates intrinsic (plasma) thromboplastin. In the extrinsic pathway, it is related to tissue thromboplastin which has the Q factor as a part of its structure. It appears to be a phospholipid obtained from exogenous sources. Both vitamins C and K have a potential redox mechanism in their structure which can be hypothecated to function in the synthesis and maintenance of mesenchymal tissue.     

Functional genetic and biophysical analyses of membrane disruption by human adenovirus

The identification of the adenovirus (AdV) protein that mediates endosome penetration during infection has remained elusive. Several lines of evidence from previous studies suggest that the membrane lytic factor of AdV is the internal capsid protein VI. While these earlier results imply a role for protein VI in endosome disruption, direct evidence during cell entry has not been demonstrated. To acquire more definitive proof, we engineered random mutations in a critical N-terminal amphipathic α-helix of VI in an attempt to generate AdV mutants that lack efficient membrane penetration and infection. Random mutagenesis within the context of the AdV genome was achieved via the development of a novel technique that incorporates both error-prone PCR and recombineering. Using this system, we identified a single mutation, L40Q, that significantly reduced infectivity and selectively impaired endosome penetration. Furthermore, we obtained biophysical data showing that the lack of efficient endosomalysis is associated with reduced insertion of the L40Q mutation in protein VI (VI-L40Q) into membranes. Our studies indicate that protein VI is the critical membrane lytic factor of AdV during cellular entry and reveal the biochemical basis for its membrane interactions.     

Does a General Temperature-Dependent Q10 Model of Soil Respiration Exist at Biome and Global Scale?

Soil respiration (SR) is commonly modeled by a Q10 (an indicator of temperature sensitivity)function in ecosystem models. Q10is usually treated as a constant of 2 in these models, although Q10 value of SR often decreases with increasing temperatures. It remains unclear whether a general temperaturedependent Q10 model of SR exists at biome and global scale. In this paper, we have compiled the long-term Q10 data of 38 SR studies ranging from the Boreal, Temperate, to Tropical/Subtropical biome on four continents.Our analysis indicated that the general temperature-dependent biome Q10 models of SR existed, especially in the Boreal and Temperate biomes. A single-exponential model was better than a simple linear model in fitting the average Q10 values at the biome scale. Average soil temperature is a better predictor of Q10 value than average air temperature in these models, especially in the Boreal biome. Soil temperature alone could explain about 50% of the Q10 variations in both the Boreal and Temperate biome single-exponential Q10 model. Q10 value of SR decreased with increasing soil temperature but at quite different rates among the three biome Q10 models. The k values (Q10 decay rate constants) were 0.09, 0.07, and 0.02/℃ in the Boreal, Temperate, and Tropical/Subtropical biome, respectively, suggesting that Q10 value is the most sensitive to soil temperature change in the Boreal biome, the second in the Temperate biome, and the least sensitive in the Tropical/Subtropical biome. This also indirectly confirms that acclimation of SR in many soil warming experiments probably occurs. The k value in the "global" single-exponential Q10 model which combined both the Boreal and Temperate biome data set was 0.08/℃. However, the global general temperature-dependent Q10model developed using the data sets of the three biomes is not adequate for predicting Q10 values of SR globally.The existence of the general temperature-dependent Q10 models of SR in the Boreal and Temperate biome has important implications for modeling SR, especially in the Boreal biome. More detail model runs are needed to exactly evaluate the impact of using a fixed Q10 vs a temperature-dependent Q10 on SR estimate in ecosystem models (e.g., TEM, Biome-BGC, and PnET).     

Defects in antigen presentation of mutant influenza haemagglutinins are reversed by mutations in the MHC class II molecule

A functional analysis was undertaken of the effects of mutating single amino acid residues in the alpha chain of the I-Ak molecule (to alanine; residues 50-79) on the ability of I-Ak transfectants to process and present influenza haemagglutinin to CD4+ T cell clones specific for two major antigenic sites of the HA1 subunit. In each instance, T cells were insensitive to a majority of substitutions in Ak with the exception of a few critical residues that differed for individual T cell clones. But more significantly, the failure of T cell clones to respond to mutant influenza viruses, containing drift substitutions within a T cell recognition site, in association with wild type I-Ak, could be reversed by single substitutions in Ak alpha. A T cell clone specific for HA1 120-139 failed to respond to a laboratory mutant virus (HA1 135 Gly----Arg) whereas optimal responses were observed with a mutant Ak transfectant (Ak alpha 56 Arg----Ala). Similarly, mutant transfectant 62 (Ak alpha 62 Gly----Ala) was able to present a natural variant virus A/TEX/77 to a T cell clone specific for HA1 48-67. We propose that Ak alpha 56 and Ak alpha 62 increase the affinity of association of mutant HA1 peptides for class II and therefore confer T cell recognition of variant viruses.     

Phase Diagrams for DNA Crystallization Systems

Abstract

Phase diagrams for several oligonucleotide duplex -spermine systems have been constructed. These diagrams characterize the duplex and spermine concentrations ranges in which crystalline precipitates are formed. All of them are wedge-like form. The slope of the upper branch of the diagram is determined by the oligonucleotide length. The position of the lower branch depends on both the nucleotide sequence and its length. The position of the lower branch depends on both the nucleotide sequence and its length. It has been shown that the addition to the system ofMgCl₂ and NaCl salts and MPD results in specific changes in the diagrams. A model for oligonucleotide duplex-spermine system has been suggested which explains the main characteristic features of the obtained phase diagrams. The experimental phase diagrams for the (pGpT)_n · (pApC)_n-spermine system (n = 2,3,4) have been analyzed ion terms of this model and the values of the binding constants of spermine and Mg²⁺ions binding to duplexes have been determined. It permitted to identify the complexes that precipitated in different regions of the phase diagrams under various conditions. The diagram obtained in the presence of a cobalt hexammine counterion is also considered. It has been shown that this phase diagram, in general, is similar to those obtained for the oligonucleotide duplex-spermine system.     

Noxa1 as a moderate activator of Nox2-based NADPH oxidase

Noxa1 was discovered as an activating factor for Nox1, an O(2)(-)-generating enzyme. Subsequent studies have shown that Noxa1 is colocalized with Nox2 in several cell types, including vascular cells. Nox2 activation by Noxa1 has been examined in reconstituted model cells. However, little is known about the kinetic properties of Noxa1 in Nox2 activation. In the present study, we used purified cyt.b(558) (Nox2 plus p22(phox)), Rac(Q61L), and Noxo1 to examine the ability of Noxa1 to activate Nox2. In the pure reconstitution system, Noxa1 activated Nox2 with lower efficiency than p67(phox), a canonical activator of Nox2. The EC(50) value of Noxa1 was considerably higher than that of p67(phox). The V(max) value with Noxa1 and Noxo1 was one-third of that with p67(phox) and p47(phox). The EC(50) value of Noxo1 or Rac(Q61L) was also higher when Noxa1 was used. The affinity of FAD for the oxidase and the stability of the active complex were remarkably low when Noxa1 and Noxo1 were used compared with p67(phox) and p47(phox). The stability was not improved by fusion of Noxa1 with Rac(Q61L). These findings show that Noxa1 has quite different kinetic properties from p67(phox) and suggest that Noxa1 may function as a moderate activator of Nox2.     

The temperature-composition phase diagram and mesophase structure characterization of monopentadecenoin in water.

The temperature-composition phase diagram of monopentadecenoin, a monoacylglycerol with a cis monounsaturated fatty acid 15 carbon atoms long (C15:1c10) in water was constructed using x-ray diffraction. Low- and wide-angle diffraction patterns were collected from samples of fixed hydration as a function of temperature in the heating direction on x-ray-sensitive film. The temperature and hydration ranges investigated were 0-104 degrees C and 0-60% (w/w) water, respectively. The phases identified in the system include the lamellar crystalline phase, the lamellar liquid crystalline phase, the fluid isotropic phase, and two inverted cubic phases belonging to space groups la3d (Q230) and Pn3m (Q244). Particular attention has been devoted to the issues of phase equilibrium, phase boundary verification, and structure characterization. The phase diagrams of monopentadecenoin, monomyristolein (C14:1c9), and monoolein (C18:1c9) are compared, and the impact of molecular structure on mesophase stability and structure is discussed.     

Development and Evaluation of Emulsions from <Emphasis Type="Italic">Carapa guianensis</Emphasis> (Andiroba) Oil

Carapa guianensis, a popular medicinal plant known as “Andiroba” in Brazil, has been used in traditional medicine as an insect repellent and anti-inflammatory product. Additionally, this seed oil has been reported in the literature as a repellent against Aedes aegypti. The aim of this work is to report on the emulsification of vegetable oils such as “Andiroba” oil by using a blend of nonionic surfactants (Span 80® and Tween 20®), using the critical hydrophilic–lipophilic balance (HLB) and pseudo-ternary diagram as tools to evaluate the system’s stability. The emulsions were prepared by the inverse phase method. Several formulations were made according to a HLB spreadsheet design (from 4.3 to 16.7), and the products were stored at 25°C and 4°C. The emulsion stabilities were tested both long- and short-term, and the more stable one was used for the pseudo-ternary diagram study. The emulsions were successfully obtained by a couple of surfactants, and the HLB analysis showed that the required HLB of the oil was 16.7. To conclude, the pseudo-ternary diagram identified several characteristic regions such as emulsion, micro-emulsion, and separation of phases.     

Identification and Structural Analysis of Amino Acid Substitutions that Increase the Stability and Activity of Aspergillus niger Glucose Oxidase

Glucose oxidase is one of the most conspicuous commercial enzymes due to its many different applications in diverse industries such as food, chemical, energy and textile. Among these applications, the most remarkable is the manufacture of glucose biosensors and in particular sensor strips used to measure glucose levels in serum. The generation of ameliorated versions of glucose oxidase is therefore a significant biotechnological objective. We have used a strategy that combined random and rational approaches to isolate uncharacterized mutations of Aspergillus niger glucose oxidase with improved properties. As a result, we have identified two changes that increase significantly the enzyme''s thermal stability. One (T554M) generates a sulfur-pi interaction and the other (Q90R/Y509E) introduces a new salt bridge near the interphase of the dimeric protein structure. An additional double substitution (Q124R/L569E) has no significant effect on stability but causes a twofold increase of the enzyme''s specific activity. Our results disclose structural motifs of the protein which are critical for its stability. The combination of mutations in the Q90R/Y509E/T554M triple mutant yielded a version of A. niger glucose oxidase with higher stability than those previously described.     

Aspergillus fumigatus来源植酸酶的Q23L和Q23LG272E突变研究

Aspergillus fumigatus来源的植酸酶具有热稳定性好、pH作用范围广的优点,但其比活性很低。设计的植酸酶Q23L突变能在pH4.5～7.0范围内大幅提高比活性,但pH稳定性却显著下降,为了进一步改良Q23L的pH稳定性,在Q23L分子上加入了G272E突变。将原酶、突变酶Q23L和突变酶Q23LG272E分别在毕赤酵母GS115中表达,表达酶经纯化后进行酶学性质比较分析,结果表明:突变酶Q23L的比活性比原酶显著提高,在pH5.5比活性由51u/mg提高到109u/mg,但其pH稳定性,尤其是在pH3.0～4.0酸性条件下的稳定性却显著降低,低于80%。突变酶Q23LG272E在pH3.0～4.5和pH6.5～7.0时的稳定性比Q23L有所提高,恢复到原酶的水平,而比活性基本维持在Q23L的水平。通过一级序列和三维结构比较,分析了可能影响Q23LG272E酶学性质的因素,为进一步研究植酸酶的结构与功能提供了材料。     

Functional significance of a highly conserved glutamate residue of the human noradrenaline transporter

The aim of the study was to investigate the role of glutamate residue 113 in transmembrane domain 2 of the human noradrenaline transporter in determining cell surface expression and functional activity. This residue is absolutely conserved in all members of the Na+- and Cl--dependent transporter family. Mutations to alanine (hE113A), aspartate (hE113D) and glutamine (hE113Q) were achieved by site-directed mutagenesis and the mutants were expressed in transfected COS-7 or HEK-293 cells. Cell surface expression of hE113A and hE113D, but not hE113Q, was markedly reduced compared with wild type, and functional noradrenaline uptake was detected only for the hE113Q mutant. The pharmacological properties of the hE113Q mutant showed very little change compared with wild type, except for a decrease in Vmax values for noradrenaline and dopamine uptake of 2-3-fold. However, the hE113D mutant showed very marked changes in its properties, compared with wild type, with 82-260-fold decreases in the affinities of the substrates, noradrenaline, dopamine and MPP+, and increased Na+ affinity for stimulation of nisoxetine binding. The results of the study show that the size and not the charge of the 113 glutamate residue of the noradrenaline transporter seems to be the most critical factor for maintenance of transporter function and surface expression.     

Stabilization of yeast hexokinase A by polyol osmolytes: correlation with the physicochemical properties of aqueous solutions

Osmolytes of the polyol series are known to accumulate in biological systems under stress and stabilize the structures of a wide variety of proteins. While increased surface tension of aqueous solutions has been considered an important factor in protein stabilization effect, glycerol is an exception, lowering the surface tension of water. To clarify this anomalous effect, the effect of a series of polyols on the thermal stability of a highly thermolabile two domain protein yeast hexokinase A has been investigated by differential scanning calorimetry and by monitoring loss in the biological activity of the enzyme as a function of time. A larger increase in the T(m) of domain 1 compared with that of domain 2, varying linearly with the number of hydroxyl groups in polyols, has been observed, sorbitol being the best stabilizer against both thermal as well as urea denaturation. Polyols help retain the activity of the enzyme considerably and a good correlation of the increase in T(m) (DeltaT(m)) and the retention of activity with the increase in the surface tension of polyol solutions, except glycerol, which breaks this trend, has been observed. However, the DeltaT(m) values show a linear correlation with apparent molal heat capacity and volume of aqueous polyol solutions including glycerol. These results suggest that while bulk solution properties contribute significantly to protein stabilization, interfacial properties are not always a good indicator of the stabilizing effect. A subtle balance of various weak binding and exclusion effects of the osmolytes mediated by water further regulates the stabilizing effect. Understanding these aspects is critical in the rational design of stable protein formulations.     

Three amino acids that are critical to formation and stability of the P22 tailspike trimer

The P22 tailspike protein folds by forming a folding competent monomer species that forms a dimeric, then a non-native trimeric (protrimer) species by addition of folding competent monomers. We have found three residues, R549, R563, and D572, which play a critical role in both the stability of the native tailspike protein and assembly and maturation of the protrimer. King and colleagues reported previously that substitution of R563 to glutamine inhibited protrimer formation. We now show that the R549Q and R563K variants significantly delay the protrimer-to-trimer transition both in vivo and in vitro. Previously, variants that destabilize intermediates have shown wild-type chemical stability. Interestingly, both the R549Q and R563K variants destabilize the tailspike trimer in guanidine denaturation studies, indicating that they represent a new class of tailspike folding variants. R549Q has a midpoint of unfolding at 3.2M guanidine, compared to 5.6M for the wild-type tailspike protein, while R563K has a midpoint of unfolding of 1.8 M. R549Q and R563K also denature over a broader pH range than the wild-type tailspike protein and both proteins have increased sensitivity to pH during refolding, suggesting that both residues are involved in ionic interactions. Our model is that R563 and D572 interact to stabilize the adjacent turn, aiding the assembly of the dimer and protrimer species. We believe that the interaction between R563 and D572 is also critical following assembly of the protrimer to properly orient D572 in order to form a salt bridge with R549 during protrimer maturation.