On Ohta's hypothesis: Most amino acid substitutions are deleterious

Ohta's hypothesis that most amino acid substitutions are deleterious grew out of a class of population-genetics models called shift models. Recently, shift models have been shown to be biologically unreasonable and have been replaced by a more plausible house-of-cards model. In this paper, the simplest form of the house-of-cards models is shown to be incompatible with most of the major features of protein evolution. Moreover, this model is shown to not be a model of exclusively deleterious-allele evolution, but rather to be a model with an equal mix of deleterious and advantageous substitutions.     

Chromosome Rearrangements Recovered following Transformation of Neurospora Crassa

Analyses of evolution and maintenance of quantitative genetic variation depend on the mutation models assumed. Currently two polygenic mutation models have been used in theoretical analyses. One is the random walk mutation model and the other is the house-of-cards mutation model. Although in the short term the two models give similar results for the evolution of neutral genetic variation within and between populations, the predictions of the changes of the variation are qualitatively different in the long term. In this paper a more general mutation model, called the regression mutation model, is proposed to bridge the gap of the two models. The model regards the regression coefficient, γ, of the effect of an allele after mutation on the effect of the allele before mutation as a parameter. When γ = 1 or 0, the model becomes the random walk model or the house-of-cards model, respectively. The additive genetic variances within and between populations are formulated for this mutation model, and some insights are gained by looking at the changes of the genetic variances as γ changes. The effects of γ on the statistical test of selection for quantitative characters during macroevolution are also discussed. The results suggest that the random walk mutation model should not be interpreted as a null hypothesis of neutrality for testing against alternative hypotheses of selection during macroevolution because it can potentially allocate too much variation for the change of population means under neutrality.     

Rational Design of Antibiotic Treatment Plans: A Treatment Strategy for Managing Evolution and Reversing Resistance

The development of reliable methods for restoring susceptibility after antibiotic resistance arises has proven elusive. A greater understanding of the relationship between antibiotic administration and the evolution of resistance is key to overcoming this challenge. Here we present a data-driven mathematical approach for developing antibiotic treatment plans that can reverse the evolution of antibiotic resistance determinants. We have generated adaptive landscapes for 16 genotypes of the TEM β-lactamase that vary from the wild type genotype “TEM-1” through all combinations of four amino acid substitutions. We determined the growth rate of each genotype when treated with each of 15 β-lactam antibiotics. By using growth rates as a measure of fitness, we computed the probability of each amino acid substitution in each β-lactam treatment using two different models named the Correlated Probability Model (CPM) and the Equal Probability Model (EPM). We then performed an exhaustive search through the 15 treatments for substitution paths leading from each of the 16 genotypes back to the wild type TEM-1. We identified optimized treatment paths that returned the highest probabilities of selecting for reversions of amino acid substitutions and returning TEM to the wild type state. For the CPM model, the optimized probabilities ranged between 0.6 and 1.0. For the EPM model, the optimized probabilities ranged between 0.38 and 1.0. For cyclical CPM treatment plans in which the starting and ending genotype was the wild type, the probabilities were between 0.62 and 0.7. Overall this study shows that there is promise for reversing the evolution of resistance through antibiotic treatment plans.     

A Mutation–Selection Model of Protein Evolution under Persistent Positive Selection

We use first principles of population genetics to model the evolution of proteins under persistent positive selection (PPS). PPS may occur when organisms are subjected to persistent environmental change, during adaptive radiations, or in host–pathogen interactions. Our mutation–selection model indicates protein evolution under PPS is an irreversible Markov process, and thus proteins under PPS show a strongly asymmetrical distribution of selection coefficients among amino acid substitutions. Our model shows the criteria ω>1 (where ω is the ratio of nonsynonymous over synonymous codon substitution rates) to detect positive selection is conservative and indeed arbitrary, because in real proteins many mutations are highly deleterious and are removed by selection even at positively selected sites. We use a penalized-likelihood implementation of the PPS model to successfully detect PPS in plant RuBisCO and influenza HA proteins. By directly estimating selection coefficients at protein sites, our inference procedure bypasses the need for using ω as a surrogate measure of selection and improves our ability to detect molecular adaptation in proteins.     

Experimental evolution of viruses: Microviridae as a model system

φX174 was developed as a model system for experimental studies of evolution because of its small genome size and ease of cultivation. It has been used extensively to address statistical questions about the dynamics of adaptive evolution. Molecular changes seen during experimental evolution of φX174 under a variety of conditions were compiled from 10 experiments comprising 58 lineages, where whole genomes were sequenced. A total of 667 substitutions was seen. Parallel evolution was rampant, with over 50 per cent of substitutions occurring at sites with three or more events. Comparisons of experimentally evolved sites to variation seen among wild phage suggest that at least some of the adaptive mechanisms seen in the laboratory are relevant to adaptation in nature. Elucidation of these mechanisms is aided by the availability of capsid and pro-capsid structures for φX174 and builds on years of genetic studies of the phage life history.     

The Amount of DNA Polymorphism Maintained in a Finite Population When the Neutral Mutation Rate Varies among Sites

The expectations of the average number of nucleotide differences per site (π), the proportion of segregating site (s), the minimum number of mutations per site (s*) and some other quantities were derived under the finite site models with and without rate variation among sites, where the finite site models include Jukes and Cantor's model, the equal-input model and Kimura's model. As a model of rate variation, the gamma distribution was used. The results indicate that if distribution parameter α is small, the effect of rate variation on these quantities are substantial, so that the estimates of θ based on the infinite site model are substantially underestimated, where θ = 4Nv, N is the effective population size and v is the mutation rate per site per generation. New methods for estimating θ are also presented, which are based on the finite site models with and without rate variation. Using these methods, underestimation can be corrected.     

Alternative Splice Variants in TIM Barrel Proteins from Human Genome Correlate with the Structural and Evolutionary Modularity of this Versatile Protein Fold

After the surprisingly low number of genes identified in the human genome, alternative splicing emerged as a major mechanism to generate protein diversity in higher eukaryotes. However, it is still not known if its prevalence along the genome evolution has contributed to the overall functional protein diversity or if it simply reflects splicing noise. The (βα)₈ barrel or TIM barrel is one of the most frequent, versatile, and ancient fold encountered among enzymes. Here, we analyze the structural modifications present in TIM barrel proteins from the human genome product of alternative splicing events. We found that 87% of all splicing events involved deletions; most of these events resulted in protein fragments that corresponded to the (βα)₂, (βα)₄, (βα)₅, (βα)₆, and (βα)₇ subdomains of TIM barrels. Because approximately 7% of all the splicing events involved internal β-strand substitutions, we decided, based on the genomic data, to design β-strand and α-helix substitutions in a well-studied TIM barrel enzyme. The biochemical characterization of one of the chimeric variants suggests that some of the splice variants in the human genome with β-strand substitutions may be evolving novel functions via either the oligomeric state or substrate specificity. We provide results of how the splice variants represent subdomains that correlate with the independently folding and evolving structural units previously reported. This work is the first to observe a link between the structural features of the barrel and a recurrent genetic mechanism. Our results suggest that it is reasonable to expect that a sizeable fraction of splice variants found in the human genome represent structurally viable functional proteins. Our data provide additional support for the hypothesis of the origin of the TIM barrel fold through the assembly of smaller subdomains. We suggest a model of how nature explores new proteins through alternative splicing as a mechanism to diversify the proteins encoded in the human genome.     

Cysteine Substitutions Define Etomidate Binding and Gating Linkages in the α-M1 Domain of γ-Aminobutyric Acid Type A (GABAA) Receptors

Etomidate is a potent general anesthetic that acts as an allosteric co-agonist at GABA_A receptors. Photoreactive etomidate derivatives labeled αMet-236 in transmembrane domain M1, which structural models locate in the β+/α- subunit interface. Other nearby residues may also contribute to etomidate binding and/or transduction through rearrangement of the site. In human α1β2γ2L GABA_A receptors, we applied the substituted cysteine accessibility method to α1-M1 domain residues extending from α1Gln-229 to α1Gln-242. We used electrophysiology to characterize each mutant''s sensitivity to GABA and etomidate. We also measured rates of sulfhydryl modification by p-chloromercuribenzenesulfonate (pCMBS) with and without GABA and tested if etomidate blocks modification of pCMBS-accessible cysteines. Cys substitutions in the outer α1-M1 domain impaired GABA activation and variably affected etomidate sensitivity. In seven of eight residues where pCMBS modification was evident, rates of modification were accelerated by GABA co-application, indicating that channel activation increases water and/or pCMBS access. Etomidate reduced the rate of modification for cysteine substitutions at α1Met-236, α1Leu-232 and α1Thr-237. We infer that these residues, predicted to face β2-M3 or M2 domains, contribute to etomidate binding. Thus, etomidate interacts with a short segment of the outer α1-M1 helix within a subdomain that undergoes significant structural rearrangement during channel gating. Our results are consistent with in silico docking calculations in a homology model that orient the long axis of etomidate approximately orthogonal to the transmembrane axis.     

Accessibility of ENaC extracellular domain central core residues

The epithelial Na⁺ channel (ENaC)/degenerin family has a similar extracellular architecture, where specific regulatory factors interact and alter channel gating behavior. The extracellular palm domain serves as a key link to the channel pore. In this study, we used cysteine-scanning mutagenesis to assess the functional effects of Cys-modifying reagents on palm domain β10 strand residues in mouse ENaC. Of the 13 ENaC α subunit mutants with Cys substitutions examined, only mutants at sites in the proximal region of β10 exhibited changes in channel activity in response to methanethiosulfonate reagents. Additionally, Cys substitutions at three proximal sites of β and γ subunit β10 strands also rendered mutant channels methanethiosulfonate-responsive. Moreover, multiple Cys mutants were activated by low concentrations of thiophilic Cd²⁺. Using the Na⁺ self-inhibition response to assess ENaC gating behavior, we identified four α, two β, and two γ subunit β10 strand mutations that changed the Na⁺ self-inhibition response. Our results suggest that the proximal regions of β10 strands in all three subunits are accessible to small aqueous compounds and Cd²⁺ and have a role in modulating ENaC gating. These results are consistent with a structural model of mouse ENaC that predicts the presence of aqueous tunnels adjacent to the proximal part of β10 and with previously resolved structures of a related family member where palm domain structural transitions were observed with channels in an open or closed state.     

Similarity in Shape Dictates Signature Intrinsic Dynamics Despite No Functional Conservation in TIM Barrel Enzymes

The conservation of the intrinsic dynamics of proteins emerges as we attempt to understand the relationship between sequence, structure and functional conservation. We characterise the conservation of such dynamics in a case where the structure is conserved but function differs greatly. The triosephosphate isomerase barrel fold (TBF), renowned for its 8 β-strand-α-helix repeats that close to form a barrel, is one of the most diverse and abundant folds found in known protein structures. Proteins with this fold have diverse enzymatic functions spanning five of six Enzyme Commission classes, and we have picked five different superfamily candidates for our analysis using elastic network models. We find that the overall shape is a large determinant in the similarity of the intrinsic dynamics, regardless of function. In particular, the β-barrel core is highly rigid, while the α-helices that flank the β-strands have greater relative mobility, allowing for the many possibilities for placement of catalytic residues. We find that these elements correlate with each other via the loops that link them, as opposed to being directly correlated. We are also able to analyse the types of motions encoded by the normal mode vectors of the α-helices. We suggest that the global conservation of the intrinsic dynamics in the TBF contributes greatly to its success as an enzymatic scaffold both through evolution and enzyme design.     

Crystal Structure of Allophycocyanin from Marine Cyanobacterium Phormidium sp. A09DM

Isolated phycobilisome (PBS) sub-assemblies have been widely subjected to X-ray crystallography analysis to obtain greater insights into the structure-function relationship of this light harvesting complex. Allophycocyanin (APC) is the phycobiliprotein always found in the PBS core complex. Phycocyanobilin (PCB) chromophores, covalently bound to conserved Cys residues of α- and β- subunits of APC, are responsible for solar energy absorption from phycocyanin and for transfer to photosynthetic apparatus. In the known APC structures, heterodimers of α- and β- subunits (known as αβ monomers) assemble as trimer or hexamer. We here for the first time report the crystal structure of APC isolated from a marine cyanobacterium (Phormidium sp. A09DM). The crystal structure has been refined against all the observed data to the resolution of 2.51 Å to R_work (R_free) of 0.158 (0.229) with good stereochemistry of the atomic model. The Phormidium protein exists as a trimer of αβ monomers in solution and in crystal lattice. The overall tertiary structures of α- and β- subunits, and trimeric quaternary fold of the Phormidium protein resemble the other known APC structures. Also, configuration and conformation of the two covalently bound PCB chromophores in the marine APC are same as those observed in fresh water cyanobacteria and marine red algae. More hydrophobic residues, however, constitute the environment of the chromophore bound to α-subunit of the Phormidium protein, owing mainly to amino acid substitutions in the marine protein.     

Signatures of neutral evolution in exponentially growing tumors: A theoretical perspective

Recent work of Sottoriva, Graham, and collaborators have led to the controversial claim that exponentially growing tumors have a site frequency spectrum that follows the 1/f law consistent with neutral evolution. This conclusion has been criticized based on data quality issues, statistical considerations, and simulation results. Here, we use rigorous mathematical arguments to investigate the site frequency spectrum in the two-type model of clonal evolution. If the fitnesses of the two types are λ₀ < λ₁, then the site frequency spectrum is c/f^α where α = λ₀/λ₁. This is due to the advantageous mutations that produce the founders of the type 1 population. Mutations within the growing type 0 and type 1 populations follow the 1/f law. Our results show that, in contrast to published criticisms, neutral evolution in an exponentially growing tumor can be distinguished from the two-type model using the site frequency spectrum.     

Assisted Evolution Enables HIV-1 to Overcome a High TRIM5α-Imposed Genetic Barrier to Rhesus Macaque Tropism

Diversification of antiretroviral factors during host evolution has erected formidable barriers to cross-species retrovirus transmission. This phenomenon likely protects humans from infection by many modern retroviruses, but it has also impaired the development of primate models of HIV-1 infection. Indeed, rhesus macaques are resistant to HIV-1, in part due to restriction imposed by the TRIM5α protein (rhTRIM5α). Initially, we attempted to derive rhTRIM5α-resistant HIV-1 strains using two strategies. First, HIV-1 was passaged in engineered human cells expressing rhTRIM5α. Second, a library of randomly mutagenized capsid protein (CA) sequences was screened for mutations that reduced rhTRIM5α sensitivity. Both approaches identified several individual mutations in CA that reduced rhTRIM5α sensitivity. However, neither approach yielded mutants that were fully resistant, perhaps because the locations of the mutations suggested that TRIM5α recognizes multiple determinants on the capsid surface. Moreover, even though additive effects of various CA mutations on HIV-1 resistance to rhTRIM5α were observed, combinations that gave full resistance were highly detrimental to fitness. Therefore, we employed an ‘assisted evolution’ approach in which individual CA mutations that reduced rhTRIM5α sensitivity without fitness penalties were randomly assorted in a library of viral clones containing synthetic CA sequences. Subsequent passage of the viral library in rhTRIM5α-expressing cells resulted in the selection of individual viral species that were fully fit and resistant to rhTRIM5α. These viruses encoded combinations of five mutations in CA that conferred complete or near complete resistance to the disruptive effects of rhTRIM5α on incoming viral cores, by abolishing recognition of the viral capsid. Importantly, HIV-1 variants encoding these CA substitutions and SIV_mac239 Vif replicated efficiently in primary rhesus macaque lymphocytes. These findings demonstrate that rhTRIM5α is difficult to but not impossible to evade, and doing so should facilitate the development of primate models of HIV-1 infection.     

Evidence for Pervasive Adaptive Protein Evolution in Wild Mice

The relative contributions of neutral and adaptive substitutions to molecular evolution has been one of the most controversial issues in evolutionary biology for more than 40 years. The analysis of within-species nucleotide polymorphism and between-species divergence data supports a widespread role for adaptive protein evolution in certain taxa. For example, estimates of the proportion of adaptive amino acid substitutions (α) are 50% or more in enteric bacteria and Drosophila. In contrast, recent estimates of α for hominids have been at most 13%. Here, we estimate α for protein sequences of murid rodents based on nucleotide polymorphism data from multiple genes in a population of the house mouse subspecies Mus musculus castaneus, which inhabits the ancestral range of the Mus species complex and nucleotide divergence between M. m. castaneus and M. famulus or the rat. We estimate that 57% of amino acid substitutions in murids have been driven by positive selection. Hominids, therefore, are exceptional in having low apparent levels of adaptive protein evolution. The high frequency of adaptive amino acid substitutions in wild mice is consistent with their large effective population size, leading to effective natural selection at the molecular level. Effective natural selection also manifests itself as a paucity of effectively neutral nonsynonymous mutations in M. m. castaneus compared to humans.     

Distinctive Patterns of Evolution of the δ-Globin Gene (HBD) in Primates

In most vertebrates, hemoglobin (Hb) is a heterotetramer composed of two dissimilar globin chains, which change during development according to the patterns of expression of α- and β-globin family members. In placental mammals, the β-globin cluster includes three early-expressed genes, ε(HBE)-γ(HBG)-ψβ(HBBP1), and the late expressed genes, δ (HBD) and β (HBB). While HBB encodes the major adult β-globin chain, HBD is weakly expressed or totally silent. Paradoxically, in human populations HBD shows high levels of conservation typical of genes under strong evolutionary constraints, possibly due to a regulatory role in the fetal-to-adult switch unique of Anthropoid primates. In this study, we have performed a comprehensive phylogenetic and comparative analysis of the two adult β-like globin genes in a set of diverse mammalian taxa, focusing on the evolution and functional divergence of HBD in primates. Our analysis revealed that anthropoids are an exception to a general pattern of concerted evolution in placental mammals, showing a high level of sequence conservation at HBD, less frequent and shorter gene conversion events. Moreover, this lineage is unique in the retention of a functional GATA-1 motif, known to be involved in the control of the developmental expression of the β-like globin genes. We further show that not only the mode but also the rate of evolution of the δ-globin gene in higher primates are strictly associated with the fetal/adult β-cluster developmental switch. To gain further insight into the possible functional constraints that have been shaping the evolutionary history of HBD in primates, we calculated dN/dS (ω) ratios under alternative models of gene evolution. Although our results indicate that HBD might have experienced different selective pressures throughout primate evolution, as shown by different ω values between apes and Old World Monkeys + New World Monkeys (0.06 versus 0.43, respectively), these estimates corroborated a constrained evolution for HBD in Anthropoid lineages, which is unlikely to be related to protein function. Collectively, these findings suggest that sequence change at the δ-globin gene has been under strong selective constraints over 65 Myr of primate evolution, likely due to a regulatory role in ontogenic switches of gene expression.     

Acid Hydrolysis and Molecular Density of Phytoglycogen and Liver Glycogen Helps Understand the Bonding in Glycogen α (Composite) Particles

Phytoglycogen (from certain mutant plants) and animal glycogen are highly branched glucose polymers with similarities in structural features and molecular size range. Both appear to form composite α particles from smaller β particles. The molecular size distribution of liver glycogen is bimodal, with distinct α and β components, while that of phytoglycogen is monomodal. This study aims to enhance our understanding of the nature of the link between liver-glycogen β particles resulting in the formation of large α particles. It examines the time evolution of the size distribution of these molecules during acid hydrolysis, and the size dependence of the molecular density of both glucans. The monomodal distribution of phytoglycogen decreases uniformly in time with hydrolysis, while with glycogen, the large particles degrade significantly more quickly. The size dependence of the molecular density shows qualitatively different shapes for these two types of molecules. The data, combined with a quantitative model for the evolution of the distribution during degradation, suggest that the bonding between β into α particles is different between phytoglycogen and liver glycogen, with the formation of a glycosidic linkage for phytoglycogen and a covalent or strong non-covalent linkage, most probably involving a protein, for glycogen as most likely. This finding is of importance for diabetes, where α-particle structure is impaired.     

The C-terminal alphaO helix of human Ogg1 is essential for 8-oxoguanine DNA glycosylase activity: the mitochondrial beta-Ogg1 lacks this domain and does not have glycosylase activity

The human Ogg1 glycosylase is responsible for repairing 8-oxo-7,8-dihydroguanine (8-oxoG) in both nuclear and mitochondrial DNA. Two distinct Ogg1 isoforms are present; α-Ogg1, which mainly localizes to the nucleus and β-Ogg1, which localizes only to mitochondria. We recently showed that mitochondria from ρ⁰ cells, which lack mitochondrial DNA, have similar 8-oxoG DNA glycosylase activity to that of wild-type cells. Here, we show that β-Ogg1 protein levels are ~80% reduced in ρ⁰ cells, suggesting β-Ogg1 is not responsible for 8-oxoG incision in mitochondria. Thus, we characterized the biochemical properties of recombinant β-Ogg1. Surprisingly, recombinant β-Ogg1 did not show any significant 8-oxoG DNA glycosylase activity in vitro. Since β-Ogg1 lacks the C-terminal αO helix present in α-Ogg1, we generated mutant proteins with various amino acid substitutions in this domain. Of the seven amino acid positions substituted (317–323), we identified Val-317 as a novel critical residue for 8-oxoG binding and incision. Our results suggest that the αO helix is absolutely necessary for 8-oxoG DNA glycosylase activity, and thus its absence may explain why β-Ogg1 does not catalyze 8-oxoG incision in vitro. Western blot analysis revealed the presence of significant amounts of α-Ogg1 in human mitochondria. Together with previous localization studies in vivo, this suggests that α-Ogg1 protein may provide the 8-oxoG DNA glycosylase activity for the repair of these lesions in human mitochondrial DNA. β-Ogg1 may play a novel role in human mitochondria.     

Low-fidelity DNA synthesis by human DNA polymerase theta

Human DNA polymerase theta (pol θ or POLQ) is a proofreading-deficient family A enzyme implicated in translesion synthesis (TLS) and perhaps in somatic hypermutation (SHM) of immunoglobulin genes. These proposed functions and kinetic studies imply that pol θ may synthesize DNA with low fidelity. Here, we show that when copying undamaged DNA, pol θ generates single base errors at rates 10- to more than 100-fold higher than for other family A members. Pol θ adds single nucleotides to homopolymeric runs at particularly high rates, exceeding 1% in certain sequence contexts, and generates single base substitutions at an average rate of 2.4 × 10⁻³, comparable to inaccurate family Y human pol κ (5.8 × 10⁻³) also implicated in TLS. Like pol κ, pol θ is processive, implying that it may be tightly regulated to avoid deleterious mutagenesis. Pol θ also generates certain base substitutions at high rates within sequence contexts similar to those inferred to be copied by pol θ during SHM of immunoglobulin genes in mice. Thus, pol θ is an exception among family A polymerases, and its low fidelity is consistent with its proposed roles in TLS and SHM.     

Formation and Function of the Rbl2p–β-Tubulin Complex

The yeast protein Rbl2p suppresses the deleterious effects of excess β-tubulin as efficiently as does α-tubulin. Both in vivo and in vitro, Rbl2p forms a complex with β-tubulin that does not contain α-tubulin, thus defining a second pool of β-tubulin in the cell. Formation of the complex depends upon the conformation of β-tubulin. Newly synthesized β-tubulin can bind to Rbl2p before it binds to α-tubulin. Rbl2p can also bind β-tubulin from the α/β-tubulin heterodimer, apparently by competing with α-tubulin. The Rbl2p–β-tubulin complex has a half-life of ~2.5 h and is less stable than the α/β-tubulin heterodimer. The results of our experiments explain both how excess Rbl2p can rescue cells overexpressing β-tubulin and how it can be deleterious in a wild-type background. They also suggest that the Rbl2p–β-tubulin complex is part of a cellular mechanism for regulating the levels and dimerization of tubulin chains.