Accelerated evolution as a consequence of transitions to mutualism

Differential rates of nucleotide substitutions among taxa are a common observation in molecular phylogenetic studies, yet links between rates of DNA evolution and traits or behaviors of organisms have proved elusive. Likelihood ratio testing is used here for the first time to evaluate specific hypotheses that account for the induction of shifts in rates of DNA evolution. A molecular phylogenetic investigation of mutualist (lichen-forming fungi and fungi associated with liverworts) and nonmutualist fungi revealed four independent transitions to mutualism. We demonstrate a highly significant association between mutualism and increased rates of nucleotide substitutions in nuclear ribosomal DNA, and we demonstrate that a transition to mutualism preceded the rate acceleration of nuclear ribosomal DNA in these lineages. Our results suggest that the increased rate of evolution after the adoption of a mutualist lifestyle is generalized across the genome of these mutualist fungi.

Somatic mutations in the p53 tumor suppressor gene in rheumatoid arthritis synovium

The factors that regulate the perpetuation and invasiveness of rheumatoid synovitis have been the subject of considerable inquiry, and the possibility that nonimmunologic defects can contribute to the disease has not been rigorously addressed. Using a mismatch detection system, we report that synovial tissue from the joints of severe chronic rheumatoid arthritis patients contain mutant p53 transcripts, which were not found in skin samples from the same patients or in joints of patients with osteoarthritis. Mutant p53 transcripts also were identified in synoviocytes cultured from rheumatoid joints. The predicted amino acid substitutions in p53 were identical or similar to those commonly observed in a variety of tumors and might influence growth and survival of rheumatoid synoviocytes. Thus, mutations in p53 and subsequent selection of the mutant cells may occur in the joints of patients as a consequence of inflammation and contribute to the pathogenesis of the disease.

Identification of the prooxidant site of human ceruloplasmin: A model for oxidative damage by copper bound to protein surfaces

Free transition metal ions oxidize lipids and lipoproteins in vitro; however, recent evidence suggests that free metal ion-independent mechanisms are more likely in vivo. We have shown previously that human ceruloplasmin (Cp), a serum protein containing seven Cu atoms, induces low density lipoprotein oxidation in vitro and that the activity depends on the presence of a single, chelatable Cu atom. We here use biochemical and molecular approaches to determine the site responsible for Cp prooxidant activity. Experiments with the His-specific reagent diethylpyrocarbonate (DEPC) showed that one or more His residues was specifically required. Quantitative [14C]DEPC binding studies indicated the importance of a single His residue because only one was exposed upon removal of the prooxidant Cu. Plasmin digestion of [14C]DEPC-treated Cp (and N-terminal sequence analysis of the fragments) showed that the critical His was in a 17-kDa region containing four His residues in the second major sequence homology domain of Cp. A full length human Cp cDNA was modified by site-directed mutagenesis to give His-to-Ala substitutions at each of the four positions and was transfected into COS-7 cells, and low density lipoprotein oxidation was measured. The prooxidant site was localized to a region containing His426 because CpH426A almost completely lacked prooxidant activity whereas the other mutants expressed normal activity. These observations support the hypothesis that Cu bound at specific sites on protein surfaces can cause oxidative damage to macromolecules in their environment. Cp may serve as a model protein for understanding mechanisms of oxidant damage by copper-containing (or -binding) proteins such as Cu, Zn superoxide dismutase, and amyloid precursor protein.

Inter- and intraspecies polymorphisms in the cholecystokinin-B/gastrin receptor alter drug efficacy

The brain cholecystokinin-B/gastrin receptor (CCK-BR) is a major target for drug development because of its postulated role in modulating anxiety, memory, and the perception of pain. Drug discovery efforts have resulted in the identification of small synthetic molecules that can selectively activate this receptor subtype. These drugs include the peptide-derived compound PD135,158 as well as the nonpeptide benzodiazepine-based ligand, L-740,093 (S enantiomer). We now report that the maximal level of receptor-mediated second messenger signaling that can be achieved by these compounds (drug efficacy) markedly differs among species homologs of the CCK-BR. Further analysis reveals that the observed differences in drug efficacy are in large part explained by single or double aliphatic amino acid substitutions between respective species homologs. This interspecies variability in ligand efficacy introduces the possibility of species differences in receptor-mediated function, an important consideration when selecting animal models for preclinical drug testing. The finding that even single amino acid substitutions can significantly affect drug efficacy prompted us to examine ligand-induced signaling by a known naturally occurring human CCK-BR variant (glutamic acid replaced by lysine in position 288; 288E → K). When examined using the 288E → K receptor, the efficacies of both PD135,158 and L-740,093 (S) were markedly increased compared with values obtained with the wild-type human protein. These observations suggest that functional variability resulting from human receptor polymorphisms may contribute to interindividual differences in drug effects.

Directed evolution of a thermostable esterase

We have used in vitro evolution to probe the relationship between stability and activity in a mesophilic esterase. Previous studies of these properties in homologous enzymes evolved for function at different temperatures have suggested that stability at high temperatures is incompatible with high catalytic activity at low temperatures through mutually exclusive demands on enzyme flexibility. Six generations of random mutagenesis, recombination, and screening stabilized Bacillus subtilis p-nitrobenzyl esterase significantly (>14°C increase in Tm) without compromising its catalytic activity at lower temperatures. Furthermore, analysis of the stabilities and activities of large numbers of random mutants indicates that these properties are not inversely correlated. Although enhanced thermostability does not necessarily come at the cost of activity, the process by which the molecule adapts is important. Mutations that increase thermostability while maintaining low-temperature activity are very rare. Unless both properties are constrained (by natural selection or screening) the evolution of one by the accumulation of single amino acid substitutions typically comes at the cost of the other, regardless of whether the two properties are inversely correlated or not correlated at all.

In vitro evolution of horse heart myoglobin to increase peroxidase activity

Random mutagenesis and screening for enzymatic activity has been used to engineer horse heart myoglobin to enhance its intrinsic peroxidase activity. A chemically synthesized gene encoding horse heart myoglobin was subjected to successive cycles of PCR random mutagenesis. The mutated myoglobin gene was expressed in Escherichia coli LE392, and the variants were screened for peroxidase activity with a plate assay. Four cycles of mutagenesis and screening produced a series of single, double, triple, and quadruple variants with enhanced peroxidase activity. Steady-state kinetics analysis demonstrated that the quadruple variant T39I/K45D/F46L/I107F exhibits peroxidase activity significantly greater than that of the wild-type protein with k1 (for H2O2 oxidation of metmyoglobin) of 1.34 × 104 M−1 s−1 (≈25-fold that of wild-type myoglobin) and k3 [for reducing the substrate (2, 2′-azino-di-(3-ethyl)benzthiazoline-6-sulfonic acid] of 1.4 × 106 M−1 s−1 (1.6-fold that of wild-type myoglobin). Thermal stability of these variants as measured with circular dichroism spectroscopy demonstrated that the Tm of the quadruple variant is decreased only slightly compared with wild-type (74.1°C vs. 76.5°C). The rate constants for binding of dioxygen exhibited by the quadruple variant are identical to the those observed for wild-type myoglobin (kon, 22.2 × 10−6 M−1 s−1 vs. 22.3 × 10−6 M−1 s−1; koff, 24.3 s−1 vs. 24.2 s−1; KO2, 0.91 × 10−6 M−1 vs. 0.92 × 10−6 M−1). The affinity of the quadruple variant for CO is increased slightly (kon, 0.90 × 10−6 M−1s−1 vs. 0.51 × 10−6 M−1s−1; koff, 5.08 s−1 vs. 3.51 s−1; KCO, 1.77 × 10−7 M−1 vs. 1.45 × 10−7 M−1). All four substitutions are in the heme pocket and within 5 Å of the heme group.

Identification of active sites in amidase: Evolutionary relationship between amide bond- and peptide bond-cleaving enzymes

Mainly based on various inhibitor studies previously performed, amidases came to be regarded as sulfhydryl enzymes. Not completely satisfied with this generally accepted interpretation, we performed a series of site-directed mutagenesis studies on one particular amidase of Rhodococcus rhodochrous J1 that was involved in its nitrile metabolism. For these experiments, the recombinant amidase was produced as the inclusion body in Escherichia coli to greatly facilitate its recovery and subsequent purification. With regard to the presumptive active site residue Cys203, a Cys203 → Ala mutant enzyme still retained 11.5% of the original specific activity. In sharp contrast, substitutions in certain other positions in the neighborhood of Cys203 had a far more dramatic effect on the amidase. Glutamic acid substitution of Asp191 reduced the specific activity of the mutant enzyme to 1.33% of the wild-type activity. Furthermore, Asp191 → Asn substitution as well as Ser195 → Ala substitution completely abolished the specific activity. It would thus appear that, among various conserved residues residing within the so-called signature sequence common to all amidases, the real active site residues are Asp191 and Ser195 rather than Cys203. Inasmuch as an amide bond (CO-NH2) in the amide substrate is not too far structurally removed from a peptide bond (CO-NH-), the signature sequences of various amidases were compared with the active site sequences of various types of proteases. It was found that aspartic acid and serine residues corresponding to Asp191 and Ser195 of the Rhodococcus amidase are present within the active site sequences of aspartic proteinases, thus suggesting the evolutionary relationship between the two.

DHR3, an ecdysone-inducible early-late gene encoding a Drosophila nuclear receptor, is required for embryogenesis

Response to the steroid hormone ecdysone in Drosophila is controlled by genetic regulatory hierarchies that include eight members of the nuclear receptor protein family. The DHR3 gene, located within the 46F early-late ecdysone-inducible chromosome puff, encodes an orphan nuclear receptor that recently has been shown to exert both positive and negative regulatory effects in the ecdysone-induced genetic hierarchies at metamorphosis. We used a reverse genetics approach to identify 11 DHR3 mutants from a pool of lethal mutations in the 46F region on the second chromosome. Two DHR3 mutations result in amino acid substitutions within the conserved DNA binding domain. Analysis of DHR3 mutants reveals that DHR3 function is required to complete embryogenesis. All DHR3 alleles examined result in nervous system defects in the embryo.

Prion protein NMR structure and species barrier for prion diseases

The structural basis of species specificity of transmissible spongiform encephalopathies, such as bovine spongiform encephalopathy or “mad cow disease” and Creutzfeldt–Jakob disease in humans, has been investigated using the refined NMR structure of the C-terminal domain of the mouse prion protein with residues 121–231. A database search for mammalian prion proteins yielded 23 different sequences for the fragment 124–226, which display a high degree of sequence identity and show relevant amino acid substitutions in only 18 of the 103 positions. Except for a unique isolated negative surface charge in the bovine protein, the amino acid differences are clustered in three distinct regions of the three-dimensional structure of the cellular form of the prion protein. Two of these regions represent potential species-dependent surface recognition sites for protein–protein interactions, which have independently been implicated from in vitro and in vivo studies of prion protein transformation. The third region consists of a cluster of interior hydrophobic side chains that may affect prion protein transformation at later stages, after initial conformational changes in the cellular protein.

A single amino acid substitution converts a carboxylesterase to an organophosphorus hydrolase and confers insecticide resistance on a blowfly

Resistance to organophosphorus (OP) insecticides is associated with decreased carboxylesterase activity in several insect species. It has been proposed that the resistance may be the result of a mutation in a carboxylesterase that simultaneously reduces its carboxylesterase activity and confers an OP hydrolase activity (the “mutant ali-esterase hypothesis”). In the sheep blowfly, Lucilia cuprina, the association is due to a change in a specific esterase isozyme, E3, which, in resistant flies, has a null phenotype on gels stained using standard carboxylesterase substrates. Here we show that an OP-resistant allele of the gene that encodes E3 differs at five amino acid replacement sites from a previously described OP-susceptible allele. Knowledge of the structure of a related enzyme (acetylcholinesterase) suggests that one of these substitutions (Gly137 → Asp) lies within the active site of the enzyme. The occurrence of this substitution is completely correlated with resistance across 15 isogenic strains. In vitro expression of two natural and two synthetic chimeric alleles shows that the Asp137 substitution alone is responsible for both the loss of E3’s carboxylesterase activity and the acquisition of a novel OP hydrolase activity. Modeling of Asp137 in the homologous position in acetylcholinesterase suggests that Asp137 may act as a base to orientate a water molecule in the appropriate position for hydrolysis of the phosphorylated enzyme intermediate.

“Switch I” mutant forms of the bacterial enhancer-binding protein NtrC that perturb the response to DNA

NtrC (nitrogen regulatory protein C) is a bacterial enhancer-binding protein of 469 residues that activates transcription by σ54-holoenzyme. A region of its transcriptional activation (central) domain that is highly conserved among homologous activators of σ54-holoenzyme—residues 206–220—is essential for interaction with this RNA polymerase: it is required for contact with the polymerase and/or for coupling the energy from ATP hydrolysis to a change in the conformation of the polymerase that allows it to form transcriptionally productive open complexes. Several mutant NtrC proteins with amino acid substitutions in this region, including NtrCA216V and NtrCG219K, have normal ATPase activity but fail in transcriptional activation. We now report that other mutant forms carrying amino acid substitutions at these same positions, NtrCA216C and NtrCG219C, are capable of activating transcription when they are not bound to a DNA template (non-DNA-binding derivatives with an altered helix–turn–helix DNA-binding motif at the C terminus of the protein) but are unable to do so when they are bound to a DNA template, whether or not it carries a specific enhancer. Enhancer DNA remains a positive allosteric effector of ATP hydrolysis, as it is for wild-type NtrC but, surprisingly, appears to have become a negative allosteric effector for some aspect of interaction with σ54-holoenzyme. The conserved region in which these amino acid substitutions occur (206–220) is equivalent to the Switch I region of a large group of purine nucleotide-binding proteins. Interesting analogies can be drawn between the Switch I region of NtrC and that of p21ras.

Horizontal gene transfer and mutation: Ngrol genes in the genome of Nicotiana glauca

Ngrol genes (NgrolB, NgrolC, NgORF13, and NgORF14) that are similar in sequence to genes in the left transferred DNA (TL-DNA) of Agrobacterium rhizogenes have been found in the genome of untransformed plants of Nicotiana glauca. It has been suggested that a bacterial infection resulted in transformation of Ngrol genes early in the evolution of the genus Nicotiana. Although the corresponding four rol genes in TL-DNA provoked hairy-root syndrome in plants, present-day N. glauca and plants transformed with Ngrol genes did not exhibit this phenotype. Sequenced complementation analysis revealed that the NgrolB gene did not induce adventitious roots because it contained two point mutations. Single-base site-directed mutagenesis at these two positions restored the capacity for root induction to the NgrolB gene. When the NgrolB, with these two base substitutions, was positioned under the control of the cauliflower mosaic virus 35S promoter (P35S), transgenic tobacco plants exhibited morphological abnormalities that were not observed in P35s-RirolB plants. In contrast, the activity of the NgrolC gene may have been conserved after an ancient infection by bacteria. Discussed is the effect of the horizontal gene transfer of the Ngrol genes and mutations in the NgrolB gene on the phenotype of ancient plants during the evolution of N. glauca.

Different population dynamics of human T cell lymphotropic virus type II in intravenous drug users compared with endemically infected tribes

The phylogeny of human T cell lymphotropic virus type II (HTLV-II) was investigated by using strains isolated from Amerindian and Pygmy tribes, in which the virus is maintained primarily through mother-to-child transmission via breast-feeding, and strains from intravenous drug users (IDUs), in which spread is mainly blood-borne via needle sharing. Molecular clock analysis showed that HTLV-II has two different evolutionary rates with the molecular clock for the virus in IDUs ticking 150–350 times faster than the one in endemically infected tribes: 2.7 × 10−4 compared with 1.71/7.31 × 10−7 nucleotide substitutions per site per year in the long terminal repeat region. This dramatic acceleration of the evolutionary rate seems to be related with the mode of transmission. Mathematical models showed the correlation of these two molecular clocks with an endemic spread of HTLV-II in infected tribes compared with the epidemic spread in IDUs. We also noted a sharp increase in the population size of the virus among IDUs during the last decades probably caused by the worldwide increase in intravenous drug use.

Allelic variation of a dehydrin gene cosegregates with chilling tolerance during seedling emergence

Dehydrins (DHNs, LEA D-11) are plant proteins present during environmental stresses associated with dehydration or low temperatures and during seed maturation. Functions of DHNs have not yet been defined. Earlier, we hypothesized that a ≈35-kDa DHN and membrane properties that reduce electrolyte leakage from seeds confer chilling tolerance during seedling emergence of cowpea (Vigna unguiculata L. Walp.) in an additive and independent manner. Evidence for this hypothesis was not rigorous because it was based on correlations of presence/absence of the DHN and slow electrolyte leakage with chilling tolerance in closely related cowpea lines that have some other genetic differences. Here, we provide more compelling genetic evidence for involvement of the DHN in chilling tolerance of cowpea. We developed near-isogenic lines by backcrossing. We isolated and determined the sequence of a cDNA corresponding to the ≈35-kDa DHN and used gene-specific oligonucleotides derived from it to test the genetic linkage between the DHN presence/absence trait and the DHN structural gene. We tested for association between the DHN presence/absence trait and both low-temperature seed emergence and electrolyte leakage. We show that allelic differences in the Dhn structural gene map to the same position as the DHN protein presence/absence trait and that the presence of the ≈35-kDa DHN is indeed associated with chilling tolerance during seedling emergence, independent of electrolyte leakage effects. Two types of allelic variation in the Dhn gene were identified in the protein-coding region, deletion of one Φ-segment from the DHN-negative lines and two single amino acid substitutions.

A zinc finger directory for high-affinity DNA recognition

We have used two monovalent phage display libraries containing variants of the Zif268 DNA-binding domain to obtain families of zinc fingers that bind to alterations in the last 4 bp of the DNA sequence of the Zif268 consensus operator, GCG TGGGCG. Affinity selection was performed by altering the Zif268 operator three base pairs at a time, and simultaneously selecting for sets of 16 related DNA sequences. In this way, only four experiments were required to select for all possible 64 combinations of DNA triplet sequences. The results show that (i) for high-affinity DNA binding in the range observed for the Zif268 wild-type complex (Kd = 0.5–5 nM), finger 1 specifically requires the arginine at the carboxy terminus of its recognition helix that forms a bidentate hydrogen-bond with the guanine base (G) in the crystal structure of Zif268 complexed to its DNA operator sequence GCG TGG GCG; (ii) when the guanine base (G) is replaced by A, C, or T, a lower-affinity family (Kd ⩾ 50 nM) can be detected that shows an overall tendency to bind G-rich DNA; (iii) the residues at position 2 on the finger 2 recognition helix do not appear to interact strongly with the complementary 5′ base in the finger 1 binding site; and (iv) unexpected substitutions at the amino terminus of finger 1 can occasionally result in specificity for the 3′ base in the finger 1 binding site. A DNA recognition directory was constructed for high-affinity zinc fingers that recognize all three bases in a DNA triplet for seven sequences of the type GNN. Similar approaches may be applied to other zinc fingers to broaden the scope of the directory.