Horizontal gene transfer among genomes: The complexity hypothesis

Increasingly, studies of genes and genomes are indicating that considerable horizontal transfer has occurred between prokaryotes. Extensive horizontal transfer has occurred for operational genes (those involved in housekeeping), whereas informational genes (those involved in transcription, translation, and related processes) are seldomly horizontally transferred. Through phylogenetic analysis of six complete prokaryotic genomes and the identification of 312 sets of orthologous genes present in all six genomes, we tested two theories describing the temporal flow of horizontal transfer. We show that operational genes have been horizontally transferred continuously since the divergence of the prokaryotes, rather than having been exchanged in one, or a few, massive events that occurred early in the evolution of prokaryotes. In agreement with earlier studies, we found that differences in rates of evolution between operational and informational genes are minimal, suggesting that factors other than rate of evolution are responsible for the observed differences in horizontal transfer. We propose that a major factor in the more frequent horizontal transfer of operational genes is that informational genes are typically members of large, complex systems, whereas operational genes are not, thereby making horizontal transfer of informational gene products less probable (the complexity hypothesis).

Quantitative modeling of stochastic systems in molecular biology by using stochastic Petri nets

An integrated understanding of molecular and developmental biology must consider the large number of molecular species involved and the low concentrations of many species in vivo. Quantitative stochastic models of molecular interaction networks can be expressed as stochastic Petri nets (SPNs), a mathematical formalism developed in computer science. Existing software can be used to define molecular interaction networks as SPNs and solve such models for the probability distributions of molecular species. This approach allows biologists to focus on the content of models and their interpretation, rather than their implementation. The standardized format of SPNs also facilitates the replication, extension, and transfer of models between researchers. A simple chemical system is presented to demonstrate the link between stochastic models of molecular interactions and SPNs. The approach is illustrated with examples of models of genetic and biochemical phenomena where the UltraSAN package is used to present results from numerical analysis and the outcome of simulations.

Phospholipase D activity is required for suppression of yeast phosphatidylinositol transfer protein defects

Yeast phosphatidylinositol transfer protein (Sec14p) function is essential for production of Golgi-derived secretory vesicles, and this requirement is bypassed by mutations in at least seven genes. Analyses of such ‘bypass Sec14p’ mutants suggest that Sec14p acts to maintain an essential Golgi membrane diacylglycerol (DAG) pool that somehow acts to promote Golgi secretory function. SPO14 encodes the sole yeast phosphatidylinositol-4,5-bisphosphate-activated phospholipase D (PLD). PLD function, while essential for meiosis, is dispensable for vegetative growth. Herein, we report specific physiological circumstances under which an unanticipated requirement for PLD activity in yeast vegetative Golgi secretory function is revealed. This PLD involvement is essential in ‘bypass Sec14p’ mutants where normally Sec14p-dependent Golgi secretory reactions are occurring in a Sec14p-independent manner. PLD catalytic activity is necessary but not sufficient for ‘bypass Sec14p’, and yeast operating under ‘bypass Sec14p’ conditions are ethanol-sensitive. These data suggest that PLD supports ‘bypass Sec14p’ by generating a phosphatidic acid pool that is somehow utilized in supporting yeast Golgi secretory function.

The roles of the two proton input channels in cytochrome c oxidase from Rhodobacter sphaeroides probed by the effects of site-directed mutations on time-resolved electrogenic intraprotein proton transfer

The crystal structures of cytochrome c oxidase from both bovine and Paracoccus denitrificans reveal two putative proton input channels that connect the heme-copper center, where dioxygen is reduced, to the internal aqueous phase. In this work we have examined the role of these two channels, looking at the effects of site-directed mutations of residues observed in each of the channels of the cytochrome c oxidase from Rhodobacter sphaeroides. A photoelectric technique was used to monitor the time-resolved electrogenic proton transfer steps associated with the photo-induced reduction of the ferryl-oxo form of heme a3 (Fe4+ = O2−) to the oxidized form (Fe3+OH−). This redox step requires the delivery of a “chemical” H+ to protonate the reduced oxygen atom and is also coupled to proton pumping. It is found that mutations in the K channel (K362M and T359A) have virtually no effect on the ferryl-oxo-to-oxidized (F-to-Ox) transition, although steady-state turnover is severely limited. In contrast, electrogenic proton transfer at this step is strongly suppressed by mutations in the D channel. The results strongly suggest that the functional roles of the two channels are not the separate delivery of chemical or pumped protons, as proposed recently [Iwata, S., Ostermeier, C., Ludwig, B. & Michel, H. (1995) Nature (London) 376, 660–669]. The D channel is likely to be involved in the uptake of both “chemical” and “pumped” protons in the F-to-Ox transition, whereas the K channel is probably idle at this partial reaction and is likely to be used for loading the enzyme with protons at some earlier steps of the catalytic cycle. This conclusion agrees with different redox states of heme a3 in the K362M and E286Q mutants under aerobic steady-state turnover conditions.

Charge transfer and transport in DNA

We explore charge migration in DNA, advancing two distinct mechanisms of charge separation in a donor (d)–bridge ({Bj})–acceptor (a) system, where {Bj} = B1,B2, … , BN are the N-specific adjacent bases of B-DNA: (i) two-center unistep superexchange induced charge transfer, d*{Bj}a → d∓{Bj}a±, and (ii) multistep charge transport involves charge injection from d* (or d+) to {Bj}, charge hopping within {Bj}, and charge trapping by a. For off-resonance coupling, mechanism i prevails with the charge separation rate and yield exhibiting an exponential dependence ∝ exp(−βR) on the d-a distance (R). Resonance coupling results in mechanism ii with the charge separation lifetime τ ∝ Nη and yield Y ≃ (1 + δ̄ Nη)−1 exhibiting a weak (algebraic) N and distance dependence. The power parameter η is determined by charge hopping random walk. Energetic control of the charge migration mechanism is exerted by the energetics of the ion pair state d∓B1±B2 … BNa relative to the electronically excited donor doorway state d*B1B2 … BNa. The realization of charge separation via superexchange or hopping is determined by the base sequence within the bridge. Our energetic–dynamic relations, in conjunction with the energetic data for d*/d− and for B/B+, determine the realization of the two distinct mechanisms in different hole donor systems, establishing the conditions for “chemistry at a distance” after charge transport in DNA. The energetic control of the charge migration mechanisms attained by the sequence specificity of the bridge is universal for large molecular-scale systems, for proteins, and for DNA.

Conformation of coenzyme pyrroloquinoline quinone and role of Ca2+ in the catalytic mechanism of quinoprotein methanol dehydrogenase

The ab initio structures of 2,7,9-tricarboxypyrroloquinoline quinone (PQQ), semiquinone (PQQH), and dihydroquinone (PQQH2) have been determined and compared with ab initio structures of the (PQQ)Ca2+, (PQQH)Ca2+, and (PQQH2)Ca2+ complexes as well as the x-ray structure of (PQQ)Ca2+ bound at the active site of the methanol dehydrogenase (MDH) of methyltropic bacteria. Plausible mechanisms for the MDH oxidation of methanol involving the (PQQ)Ca2+ complex are explored via ab initio computations and discussed. Considering the reaction of methanol with PQQ in the absence of Ca2+, nucleophilic addition of methanol to the PQQ C-5 carbonyl followed by a retro-ene elimination is deemed unlikely due to large energy barrier. A much more favorable disposition of the methanol C-5 adduct to provide formaldehyde involves proton ionization of the intermediate followed by elimination of methoxide concerted with hydride transfer to the oxygen of the C-4 carbonyl. Much the same transition state is reached if one searches for the transition state beginning with Asp-303–CO2−general-base removal of the methanol proton of the (PQQ)Ca2+O(H)CH3 complex concerted with hydride transfer to the oxygen at C-4. For such a mechanism the role of the Ca2+ moiety would be to (i) contribute to the formation of the ES complex (ii) provide a modest decrease in the pKa of methanol substrate,; and (iii) polarize the oxygen at C-5.

Disruption of the telomerase catalytic subunit gene from Arabidopsis inactivates telomerase and leads to a slow loss of telomeric DNA

Telomerase is an essential enzyme that maintains telomeres on eukaryotic chromosomes. In mammals, telomerase is required for the lifelong proliferative capacity of normal regenerative and reproductive tissues and for sustained growth in a dedifferentiated state. Although the importance of telomeres was first elucidated in plants 60 years ago, little is known about the role of telomeres and telomerase in plant growth and development. Here we report the cloning and characterization of the Arabidopsis telomerase reverse transcriptase (TERT) gene, AtTERT. AtTERT is predicted to encode a highly basic protein of 131 kDa that harbors the reverse transcriptase and telomerase-specific motifs common to all known TERT proteins. AtTERT mRNA is 10–20 times more abundant in callus, which has high levels of telomerase activity, versus leaves, which contain no detectable telomerase. Plants homozygous for a transfer DNA insertion into the AtTERT gene lack telomerase activity, confirming the identity and function of this gene. Because telomeres in wild-type Arabidopsis are short, the discovery that telomerase-null plants are viable for at least two generations was unexpected. In the absence of telomerase, telomeres decline by approximately 500 bp per generation, a rate 10 times slower than seen in telomerase-deficient mice. This gradual loss of telomeric DNA may reflect a reduced rate of nucleotide depletion per round of DNA replication, or the requirement for fewer cell divisions per organismal generation. Nevertheless, progressive telomere shortening in the mutants, however slow, ultimately should be lethal.

Growth suppression of glioma cells by PTEN requires a functional phosphatase catalytic domain

Deletions of all or part of chromosome 10 are the most common genetic alterations in high-grade gliomas. The PTEN gene (also called MMAC1 and TEP1) maps to chromosome region 10q23 and has been implicated as a target of alteration in gliomas and also in other cancers such as those of the breast, prostate, and kidney. Here we sought to provide a functional test of its candidacy as a growth suppressor in glioma cells. We used a combination of Northern blot analysis, protein truncation assays, and sequence analysis to determine the types and frequency of PTEN mutations in glioma cell lines so that we could define appropriate recipients to assess the growth suppressive function of PTEN by gene transfer. Introduction of wild-type PTEN into glioma cells containing endogenous mutant alleles caused growth suppression, but was without effect in cells containing endogenous wild-type PTEN. The ectopic expression of PTEN alleles, which carried mutations found in primary tumors and have been shown or are expected to inactivate its phosphatase activity, caused little growth suppression. These data strongly suggest that PTEN is a protein phosphatase that exhibits functional and specific growth-suppressing activity.

Conformational changes between the active-site and regulatory light chain of myosin as determined by luminescence resonance energy transfer: The effect of nucleotides and actin

Myosin is thought to generate movement of actin filaments via a conformational change between its light-chain domain and its catalytic domain that is driven by the binding of nucleotides and actin. To monitor this change, we have measured distances between a gizzard regulatory light chain (Cys 108) and the active site (near or at Trp 130) of skeletal myosin subfragment 1 (S1) by using luminescence resonance energy transfer and a photoaffinity ATP-lanthanide analog. The technique allows relatively long distances to be measured, and the label enables site-specific attachment at the active-site with only modest affect on myosin’s enzymology. The distance between these sites is 66.8 ± 2.3 Å when the nucleotide is ADP and is unchanged on binding to actin. The distance decreases slightly with ADP-BeF3, (−1.6 ± 0.3 Å) and more significantly with ADP-AlF4 (−4.6 ± 0.2 Å). During steady-state hydrolysis of ATP, the distance is temperature-dependent, becoming shorter as temperature increases and the complex with ADP⋅Pi is favored over that with ATP. We conclude that the distance between the active site and the light chain varies as Acto-S1-ADP ≈ S1-ADP > S1-ADP-BeF3 > S1-ADP-AlF4 ≈ S1-ADP-Pi and that S1-ATP > S1-ADP-Pi. The changes in distance are consistent with a substantial rotation of the light-chain binding domain of skeletal S1 between the prepowerstroke state, simulated by S1-ADP-AlF4, and the post-powerstroke state, simulated by acto-S1-ADP.

Correct protein folding in glycerol

Water is the natural medium for protein folding, which is also used in all in vitro studies. In the present work, we posed, and answered affirmatively, a question of whether it is possible to fold correctly a typical protein in a nonaqueous solvent. To this end, unfolded and reduced hen egg-white lysozyme was refolded and reoxidized in glycerol containing varying amounts of water. The unfolded/reduced enzyme was found to regain spontaneously substantial catalytic activity even in the nearly anhydrous solvent; for example, the refolding yield in 99% glycerol was still some one-third of that in pure water, and one-half of that was regained even in 99.8% glycerol. The less than full recovery of the enzymatic activity in glycerol is, as in water, because of competing protein aggregation during the refolding. Lysozyme reoxidation in glycerol was successfully mediated by two dissimilar oxidizing systems, and the refolding yield was markedly affected by the pH of the last aqueous solution before the transfer into glycerol. No recovery of the lysozyme activity was observed when the refolding/reoxidation reaction was carried out in the denaturing solvent dimethyl sulfoxide. This study paves the way for a systematic investigation of the solvent effect on protein folding and demonstrates that water is not a unique milieu for this process.

Comparative genomics of the restriction-modification systems in Helicobacter pylori

Helicobacter pylori is a Gram-negative bacterial pathogen with a small genome of 1.64–1.67 Mb. More than 20 putative DNA restriction-modification (R-M) systems, comprising more than 4% of the total genome, have been identified in the two completely sequenced H. pylori strains, 26695 and J99, based on sequence similarities. In this study, we have investigated the biochemical activities of 14 Type II R-M systems in H. pylori 26695. Less than 30% of the Type II R-M systems in 26695 are fully functional, similar to the results obtained from strain J99. Although nearly 90% of the R-M genes are shared by the two H. pylori strains, different sets of these R-M genes are functionally active in each strain. Interestingly, all strain-specific R-M genes are active, whereas most shared genes are inactive. This agrees with the notion that strain-specific genes have been acquired more recently through horizontal transfer from other bacteria and selected for function. Thus, they are less likely to be impaired by random mutations. Our results also show that H. pylori has extremely diversified R-M systems in different strains, and that the diversity may be maintained by constantly acquiring new R-M systems and by inactivating and deleting the old ones.

Toward a quantum-mechanical description of metal-assisted phosphoryl transfer in pyrophosphatase

The wealth of kinetic and structural information makes inorganic pyrophosphatases (PPases) a good model system to study the details of enzymatic phosphoryl transfer. The enzyme accelerates metal-complexed phosphoryl transfer 1010-fold: but how? Our structures of the yeast PPase product complex at 1.15 Å and fluoride-inhibited complex at 1.9 Å visualize the active site in three different states: substrate-bound, immediate product bound, and relaxed product bound. These span the steps around chemical catalysis and provide strong evidence that a water molecule (Onu) directly attacks PPi with a pKa vastly lowered by coordination to two metal ions and D117. They also suggest that a low-barrier hydrogen bond (LBHB) forms between D117 and Onu, in part because of steric crowding by W100 and N116. Direct visualization of the double bonds on the phosphates appears possible. The flexible side chains at the top of the active site absorb the motion involved in the reaction, which may help accelerate catalysis. Relaxation of the product allows a new nucleophile to be generated and creates symmetry in the elementary catalytic steps on the enzyme. We are thus moving closer to understanding phosphoryl transfer in PPases at the quantum mechanical level. Ultra-high resolution structures can thus tease out overlapping complexes and so are as relevant to discussion of enzyme mechanism as structures produced by time-resolved crystallography.

On the role of the K-proton transfer pathway in cytochrome c oxidase

Cytochrome c oxidase is a membrane-bound enzyme that catalyzes the four-electron reduction of oxygen to water. This highly exergonic reaction drives proton pumping across the membrane. One of the key questions associated with the function of cytochrome c oxidase is how the transfer of electrons and protons is coupled and how proton transfer is controlled by the enzyme. In this study we focus on the function of one of the proton transfer pathways of the R. sphaeroides enzyme, the so-called K-proton transfer pathway (containing a highly conserved Lys(I-362) residue), leading from the protein surface to the catalytic site. We have investigated the kinetics of the reaction of the reduced enzyme with oxygen in mutants of the enzyme in which a residue [Ser(I-299)] near the entry point of the pathway was modified with the use of site-directed mutagenesis. The results show that during the initial steps of oxygen reduction, electron transfer to the catalytic site (to form the “peroxy” state, Pr) requires charge compensation through the proton pathway, but no proton uptake from the bulk solution. The charge compensation is proposed to involve a movement of the K(I-362) side chain toward the binuclear center. Thus, in contrast to what has been assumed previously, the results indicate that the K-pathway is used during oxygen reduction and that K(I-362) is charged at pH ≈ 7.5. The movement of the Lys is proposed to regulate proton transfer by “shutting off” the protonic connectivity through the K-pathway after initiation of the O2 reduction chemistry. This “shutoff” prevents a short-circuit of the proton-pumping machinery of the enzyme during the subsequent reaction steps.

Neural crest-directed gene transfer demonstrates Wnt1 role in melanocyte expansion and differentiation during mouse development

Wnt1 signaling has been implicated as one factor involved in neural crest-derived melanocyte (NC-M) development. Mice deficient for both Wnt1 and Wnt3a have a marked deficiency in trunk neural crest derivatives including NC-Ms. We have used cell lineage-directed gene targeting of Wnt signaling genes to examine the effects of Wnt signaling in mouse neural crest development. Gene expression was directed to cell lineages by infection with subgroup A avian leukosis virus vectors in lines of transgenic mice that express the retrovirus receptor tv-a. Transgenic mice with tva in either nestin-expressing neural precursor cells (line Ntva) or dopachrome tautomerase (DCT)-expressing melanoblasts (line DCTtva) were analyzed. We overstimulated Wnt signaling in two ways: directed gene transfer of Wnt1 to Ntva+ cells and transfer of β-catenin to DCTtva+ NC-M precursor cells. In both methods, NC-M expansion and differentiation were effected. Significant increases were observed in the number of NC-Ms [melanin+ and tyrosinase-related protein 1 (TYRP1)+ cells], the differentiation of melanin− TYRP1+ cells to melanin+ TYRP1+ NC-Ms, and the intensity of pigmentation per NC-M. These data are consistent with Wnt1 signaling being involved in both expansion and differentiation of migrating NC-Ms in the developing mouse embryo. The use of lineage-directed gene targeting will allow the dissection of signaling molecules involved in NC development and is adaptable to other mammalian developmental systems.

Efficient transfer, integration, and sustained long-term expression of the transgene in adult rat brains injected with a lentiviral vector.

We describe the construction of a safe, replication-defective and efficient lentiviral vector suitable for in vivo gene delivery. The reverse transcription of the vector was found to be a rate-limiting step; therefore, promoting the reaction inside the vector particles before delivery significantly enhanced the efficiency of gene transfer. After injection into the brain of adult rats, sustained long-term expression of the transgene was obtained in the absence of detectable pathology. A high proportion of the neurons in the areas surrounding the injection sites of the vector expressed the transduced beta-galactosidase gene. This pattern was invariant in animals sacrificed several months after a single administration of the vector. Transduction occurs by integration of the vector genome, as it was abolished by a single amino acid substitution in the catalytic site of the integrase protein incorporated in the vector. Development of clinically acceptable derivatives of the lentiviral vector may thus enable the sustained delivery of significant amounts of a therapeutic gene product in a wide variety of somatic tissues.