Epistasis as a Determinant of the HIV-1 Protease's Robustness to Mutation

The robustness of phenotypes to mutation is critical to protein evolution; robustness may be an adaptive trait if it promotes evolution. We hypothesised that native proteins subjected to natural selection in vivo should be more robust than proteins generated in vitro in the absence of natural selection. We compared the mutational robustness of two human immunodeficiency virus type 1 (HIV-1) proteases with comparable catalytic efficiencies, one isolated from an infected individual and the second generated in vitro via random mutagenesis. Single mutations in the protease (82 and 60 in the wild-type and mutant backgrounds, respectively) were randomly generated in vitro and the catalytic efficiency of each mutant was determined. No differences were observed between these two protease variants when lethal, neutral, and deleterious mutations were compared (P = 0.8025, chi-squared test). Similarly, average catalytic efficiency (−72.6% and −64.5%, respectively) did not significantly differ between protease mutant libraries (P = 0.3414, Mann Whitney test). Overall, the two parental proteins displayed similar mutational robustness. Importantly, strong and widespread epistatic interactions were observed when the effect of the same mutation was compared in both proteases, suggesting that epistasis can be a key determinant of the robustness displayed by the in vitro generated protease.     

Effects of Topology and Sequence in Protein Folding Linked via Conformational Fluctuations

Experiments have compared the folding of proteins with different amino acid sequences but the same basic structure, or fold. Results indicate that folding is robust to sequence variations for proteins with some nonlocal folds, such as all-β, whereas the folding of more local, all-α proteins typically exhibits a stronger sequence dependence. Here, we use a coarse-grained model to systematically study how variations in sequence perturb the folding energy landscapes of three model sequences with 3α, 4β + α, and β-barrel folds, respectively. These three proteins exhibit folding features in line with experiments, including expected rank order in the cooperativity of the folding transition and stability-dependent shifts in the location of the free-energy barrier to folding. Using a generalized-ensemble simulation approach, we determine the thermodynamics of around 2000 sequence variants representing all possible hydrophobic or polar single- and double-point mutations. From an analysis of the subset of stability-neutral mutations, we find that folding is perturbed in a topology-dependent manner, with the β-barrel protein being the most robust. Our analysis shows, in particular, that the magnitude of mutational perturbations of the transition state is controlled in part by the size or “width” of the underlying conformational ensemble. This result suggests that the mutational robustness of the folding of the β-barrel protein is underpinned by its conformationally restricted transition state ensemble, revealing a link between sequence and topological effects in protein folding.     

Probing protein stability and proteolytic resistance by loop scanning: a comprehensive mutational analysis

Improvement in protein thermostability was often found to be associated with increase in its proteolytic resistance as revealed by comparative studies of homologous proteins from extremophiles or mutational studies. Structural elements of protein responsible for this association are not firmly established although loops are implicated indirectly due to their structural role in protein stability. To get a better insight, a detailed study of protein wide mutants and their influence on stability and proteolytic resistance would be helpful. To generate such a data set, a model protein, Bacillus subtilis lipase was subjected to loop scanning site-saturation mutagenesis on 86 positions spanning all loops including termini. Upon screening of ∼16,000 clones, 17 single mutants with improved thermostability were identified with increment in apparent melting temperature (Tm_app) by 1–6°C resulting in an increase in free energy of unfolding (ΔG_unf) by 0.04–1.16 kcal/mol. Proteolytic resistance of all single mutants upon incubation with nonspecific protease, Subtilisin A, was determined. Upon comparison, post-proteolysis residual activities as well as kinetics of proteolysis of mutants showed excellent correlation with ΔG_unf, (r > 0.9), suggesting that proteolysis was strongly correlated with the global stability of this protein. This significant correlation in this set, with least possible sequence changes (single aa substitution), while covering >60% of protein surface strongly argues for the covariance of these two variables. Compared to studies from extremophiles, with large sequence heterogeneity, the observed correlation in such a narrow sequence space (ΔΔG_unf = 1.57 kcal⁻¹) justifies the robustness of this relation.     

Simple Sequence Repeats and Mucoid Conversion: Biased mucA Mutagenesis in Mismatch Repair-Deficient Pseudomonas aeruginosa

In Pseudomonas aeruginosa, conversion to the mucoid phenotype marks the onset of an irreversible state of the infection in Cystic Fibrosis (CF) patients. The main pathway for mucoid conversion is mutagenesis of the mucA gene, frequently due to −1 bp deletions in a simple sequence repeat (SSR) of 5 Gs (G⁵-SSR₄₂₆). We have recently observed that this mucA mutation is particularly accentuated in Mismatch Repair System (MRS)-deficient cells grown in vitro. Interestingly, previous reports have shown a high prevalence of hypermutable MRS-deficient strains occurring naturally in CF chronic lung infections. Here, we used mucA as a forward mutation model to systematically evaluate the role of G⁵-SSR₄₂₆ in conversion to mucoidy in a MRS-deficient background, with this being the first analysis combining SSR-dependent localized hypermutability and the acquisition of a particular virulence/persistence trait in P. aeruginosa. In this study, mucA alleles were engineered with different contents of G:C SSRs, and tested for their effect on the mucoid conversion frequency and mucA mutational spectra in a mutS-deficient strain of P. aeruginosa. Importantly, deletion of G⁵-SSR₄₂₆ severely reduced the emergence frequency of mucoid variants, with no preferential site of mutagenesis within mucA. Moreover, although mutagenesis in mucA was not totally removed, this was no longer the main pathway for mucoid conversion, suggesting that G⁵-SSR₄₂₆ biased mutations towards mucA. Mutagenesis in mucA was restored by the addition of a new SSR (C⁶-SSR₄₃₁), and even synergistically increased when G⁵-SSR₄₂₆ and C⁶-SSR₄₃₁ were present simultaneously, with the mucA mutations being restricted to −1 bp deletions within any of both G:C SSRs. These results confirm a critical role for G⁵-SSR₄₂₆ enhancing the mutagenic process of mucA in MRS-deficient cells, and shed light on another mechanism, the SSR- localized hypermutability, contributing to mucoid conversion in P. aeruginosa.     

Effective Inhibition of Human Immunodeficiency Virus 1 Replication by Engineered RNase P Ribozyme

Using an in vitro selection procedure, we have previously isolated RNase P ribozyme variants that efficiently cleave an mRNA sequence in vitro. In this study, a variant was used to target the HIV RNA sequence in the tat region. The variant cleaved the tat RNA sequence in vitro about 20 times more efficiently than the wild type ribozyme. Our results provide the first direct evidence that combined mutations at nucleotide 83 and 340 of RNase P catalytic RNA from Escherichia coli (G₈₃ -> U₈₃ and G₃₄₀ -> A₃₄₀) increase the overall efficiency of the ribozyme in cleaving an HIV RNA sequence. Moreover, the variant is more effective in reducing HIV-1 p24 expression and intracellular viral RNA level in cells than the wild type ribozyme. A reduction of about 90% in viral RNA level and a reduction of 150 fold in viral growth were observed in cells that expressed the variant, while a reduction of less than 10% was observed in cells that either did not express the ribozyme or produced a catalytically inactive ribozyme mutant. Thus, engineered ribozyme variants are effective in inhibiting HIV infection. These results also demonstrate the potential of engineering RNase P ribozymes for anti-HIV application.     

Perspectives on protein evolution from simple exact models

Understanding the evolution of biopolymers is important to rationalise the directed and undirected design of functional molecules. Large scale experiments or detailed computational studies are often impractical. Therefore, simple model systems, such as RNA secondary structure and lattice proteins have been adapted to study general statistical and topological features of genotype (sequence) to phenotype (structure) maps. We review findings from such models that address aspects of thermodynamic and mutational robustness, neutral evolution and recombination of proteins. We compare various modelling approaches, and discuss their generality, parameter dependency and experimental verifications of their predictions. The most striking observation is the universal emergence of neutral nets--sets of phenotypically identical genotypes that are interconnected by series of point mutations. However, fast adaptation by point mutations appears to be problematic for proteins. This may explain why proteins appear to be more specific while RNA is rather versatile. This could even be the reason why RNA had to evolve before proteins. Similar principles of biological organisation are reflected in sequence and structure databases of real proteins. Insights gained from modelling are useful for designing more efficient database organisation and search strategies.     

Thermodynamics of neutral protein evolution

Naturally evolving proteins gradually accumulate mutations while continuing to fold to stable structures. This process of neutral evolution is an important mode of genetic change and forms the basis for the molecular clock. We present a mathematical theory that predicts the number of accumulated mutations, the index of dispersion, and the distribution of stabilities in an evolving protein population from knowledge of the stability effects (delta deltaG values) for single mutations. Our theory quantitatively describes how neutral evolution leads to marginally stable proteins and provides formulas for calculating how fluctuations in stability can overdisperse the molecular clock. It also shows that the structural influences on the rate of sequence evolution observed in earlier simulations can be calculated using just the single-mutation delta deltaG values. We consider both the case when the product of the population size and mutation rate is small and the case when this product is large, and show that in the latter case the proteins evolve excess mutational robustness that is manifested by extra stability and an increase in the rate of sequence evolution. All our theoretical predictions are confirmed by simulations with lattice proteins. Our work provides a mathematical foundation for understanding how protein biophysics shapes the process of evolution.     

Emergence,Retention and Selection: A Trilogy of Origination for Functional De Novo Proteins from Ancestral LncRNAs in Primates

While some human-specific protein-coding genes have been proposed to originate from ancestral lncRNAs, the transition process remains poorly understood. Here we identified 64 hominoid-specific de novo genes and report a mechanism for the origination of functional de novo proteins from ancestral lncRNAs with precise splicing structures and specific tissue expression profiles. Whole-genome sequencing of dozens of rhesus macaque animals revealed that these lncRNAs are generally not more selectively constrained than other lncRNA loci. The existence of these newly-originated de novo proteins is also not beyond anticipation under neutral expectation, as they generally have longer theoretical lifespan than their current age, due to their GC-rich sequence property enabling stable ORFs with lower chance of non-sense mutations. Interestingly, although the emergence and retention of these de novo genes are likely driven by neutral forces, population genetics study in 67 human individuals and 82 macaque animals revealed signatures of purifying selection on these genes specifically in human population, indicating a proportion of these newly-originated proteins are already functional in human. We thus propose a mechanism for creation of functional de novo proteins from ancestral lncRNAs during the primate evolution, which may contribute to human-specific genetic novelties by taking advantage of existed genomic contexts.     

The effect of recombination on the neutral evolution of genetic robustness

Conventional population genetics considers the evolution of a limited number of genotypes corresponding to phenotypes with different fitness. As model phenotypes, in particular RNA secondary structure, have become computationally tractable, however, it has become apparent that the context dependent effect of mutations and the many-to-one nature inherent in these genotype-phenotype maps can have fundamental evolutionary consequences. It has previously been demonstrated that populations of genotypes evolving on the neutral networks corresponding to all genotypes with the same secondary structure only through neutral mutations can evolve mutational robustness [E. van Nimwegen, J.P. Crutchfield, M. Huynen, Neutral evolution of mutational robustness, Proc. Natl. Acad. Sci. USA 96(17), 9716-9720 (1999)], by concentrating the population on regions of high neutrality. Introducing recombination we demonstrate, through numerically calculating the stationary distribution of an infinite population on ensembles of random neutral networks that mutational robustness is significantly enhanced and further that the magnitude of this enhancement is sensitive to details of the neutral network topology. Through the simulation of finite populations of genotypes evolving on random neutral networks and a scaled down microRNA neutral network, we show that even in finite populations recombination will still act to focus the population on regions of locally high neutrality.     

Evolution of mutational robustness

We review recent advances in the understanding of the mutation-selection balance of asexual replicators. For over 30 years, population geneticists thought that an expression derived by Kimura and Maruyama in 1966 fully solved this problem. However, Kimura and Maruyama's result is only correct in the absence of neutral mutations. The inclusion of neutral mutations leads to a wealth of interesting new effects, and, in particular, to a selective pressure to evolve robustness against mutations. We cover recent literature on the population dynamics of asexual replicators on networks of neutral genotypes, on the outcompetition of fast replicators by slower ones with better mutational support, and on the probability of fixation at high mutation rates. We discuss empirical evidence for the evolution of mutational robustness, and speculate on its relevance for higher organisms.     

Revisiting Robustness and Evolvability: Evolution in Weighted Genotype Spaces

Robustness and evolvability are highly intertwined properties of biological systems. The relationship between these properties determines how biological systems are able to withstand mutations and show variation in response to them. Computational studies have explored the relationship between these two properties using neutral networks of RNA sequences (genotype) and their secondary structures (phenotype) as a model system. However, these studies have assumed every mutation to a sequence to be equally likely; the differences in the likelihood of the occurrence of various mutations, and the consequence of probabilistic nature of the mutations in such a system have previously been ignored. Associating probabilities to mutations essentially results in the weighting of genotype space. We here perform a comparative analysis of weighted and unweighted neutral networks of RNA sequences, and subsequently explore the relationship between robustness and evolvability. We show that assuming an equal likelihood for all mutations (as in an unweighted network), underestimates robustness and overestimates evolvability of a system. In spite of discarding this assumption, we observe that a negative correlation between sequence (genotype) robustness and sequence evolvability persists, and also that structure (phenotype) robustness promotes structure evolvability, as observed in earlier studies using unweighted networks. We also study the effects of base composition bias on robustness and evolvability. Particularly, we explore the association between robustness and evolvability in a sequence space that is AU-rich – sequences with an AU content of 80% or higher, compared to a normal (unbiased) sequence space. We find that evolvability of both sequences and structures in an AU-rich space is lesser compared to the normal space, and robustness higher. We also observe that AU-rich populations evolving on neutral networks of phenotypes, can access less phenotypic variation compared to normal populations evolving on neutral networks.     

Neutralism and selectionism: a network-based reconciliation

Neutralism and selectionism are extremes of an explanatory spectrum for understanding patterns of molecular evolution and the emergence of evolutionary innovation. Although recent genome-scale data from protein-coding genes argue against neutralism, molecular engineering and protein evolution data argue that neutral mutations and mutational robustness are important for evolutionary innovation. Here I propose a reconciliation in which neutral mutations prepare the ground for later evolutionary adaptation. Key to this perspective is an explicit understanding of molecular phenotypes that has only become accessible in recent years.     

Structures of complexes of 5S RNA with ribosomal proteins L5, L18 and L25 from Escherichia coli: identification of kethoxal-reactive sites on the 5S RNA.

The topography of Escherichia coli 5S RNA has been examined in the presence of ribosomal proteins L5, L18 and L25 and their different combinations, by comparing the kethoxal modification characteristics of the various RNA-protein complexes with those of the free A-conformer of 5S RNA (Noller &; Garrett, 1979, accompanying paper).Two of the four most reactive guanines, G₁₃ and G₄₁, are unaffected by the protein, in accord with the finding that these are the only two guanines that are accessible in the 50S subunit (Noller &; Herr, 1974). The other two very reactive guanines, G₂₄ and G₆₉, are strongly protected by protein L18, either in the presence or absence of proteins L5 and L25. Protein binding studies with kethoxal-modified 5S RNA provide evidence that one or both of these two guanines are directly involved in the protein-RNA interactions, and this conclusion is supported by the occurrence of guanines in these two positions in all the other sequenced prokaryotic 5S RNAs.The group of less reactive guanines, G₁₆, G₂₃, G₄₄, G₈₆ and G₁₀₇, are protected to some extent by each of the proteins L5, L18 and L25; the strongest effect is with L18. We suggest that this is attributable to a small increase in the conformational homogeneity of the 5S RNA and that L18, in particular, induces some tightening of the RNA structure.Only one guanine, G₆₉, is rendered more accessible by the proteins. This effect is produced by protein L25, which is known to cause some destructuring of the 5S RNA (Bear et al., 1977). There was no other evidence for any destructuring of the 5S RNA. In particular, the sequence 72 to 83, which is complementary to a sequence in 23S RNA (Herr &; Noller, 1975), is not modified. However, in contrast to an earlier report (Erdmann et al., 1973), the conserved sequence G₄₄-A-A-C, which has been implicated in tRNA binding, was not rendered more accessible by the proteins.     

Genetic variation of blood proteins within and between local populations of grivet monkey (Cercopithecus aethiops aethiops) in central Ethiopia

A population genetic study by blood protein electrophoresis revealed that populations of wild grivet monkeys in central Ethiopia show a comparatively low level of variability and less differentiation among local populations over broad geographical areas. This is evaluated by comparison with other wild primate studies using the same electrophoretic technique. A total of 196 blood samples, collected from 10 local populations comprising 11 troops distributed along approximately 600 km of the Awash River, were examined for 33 genetic loci. The low level of variability was indicated by the proportion of polymorphic loci (P_poly), which was on average 11% with an average heterozygosity (H) of 3%. A tendency for lower genetic differentiation among local populations was shown by theG _ST value of 0.09, an averageF _ST of 0.08, andNei's genetic distance; ranging from 0.002 to 0.023. Considering paleoclimatological studies of the area and ecological traits of this species, these findings can be explained as a consequence of a comparatively recent and repeated series of rapid habitat expansions following severe climatic conditions.     

Nucleotide excision repair of 2-acetylaminofluorene- and 2-aminofluorene-(C8)-guanine adducts: molecular dynamics simulations elucidate how lesion structure and base sequence context impact repair efficiencies

Nucleotide excision repair (NER) efficiencies of DNA lesions can vary by orders of magnitude, for reasons that remain unclear. An example is the pair of N-(2′-deoxyguanosin-8-yl)-2-aminofluorene (dG-C8-AF) and N-(2′-deoxyguanosin-8-yl)-2-acetylaminofluorene (dG-C8-AAF) adducts that differ by a single acetyl group. The NER efficiencies in human HeLa cell extracts of these lesions are significantly different when placed at G₁, G₂ or G₃ in the duplex sequence (5′-CTCG₁G₂CG₃CCATC-3′) containing the NarI mutational hot spot. Furthermore, the dG-C8-AAF adduct is a better substrate of NER than dG-C8-AF in all three NarI sequence contexts. The conformations of each of these adducts were investigated by Molecular dynamics (MD) simulation methods. In the base-displaced conformational family, the greater repair susceptibility of dG-C8-AAF in all sequences stems from steric hindrance effects of the acetyl group which significantly diminish the adduct-base stabilizing van der Waals stacking interactions relative to the dG-C8-AF case. Base sequence context effects for each adduct are caused by differences in helix untwisting and minor groove opening that are derived from the differences in stacking patterns. Overall, the greater NER efficiencies are correlated with greater extents of base sequence-dependent local untwisting and minor groove opening together with weaker stacking interactions.     

Cloning and analysis of the cDNA encoding the rod G-protein transducin α, β1 and γ1 subunits from the canine retina

The canine (Canis familiaris) retinal rod transducin (G_T) α, β1 and γ1 subunits were sequenced. Cloning of the cDNAs was accomplished by polymerase chain reaction (PCR) using degenerate and wild type retinal cDNA libraries as templates. The deduced amino acid sequences were highly similar to rod transducins from other species: G_Tα differed by 5 amino acids from the corresponding human sequence, whereas β1 and γ1 were identical to human sequences. The coding sequence of rod transducin was evaluated as a possible cause for the recessively inherited retinal rod-cone degeneration: there were no nucleotide differences between the wild type and retinal degenerate strains.     

The Effects of Demography and Long-Term Selection on the Accuracy of Genomic Prediction with Sequence Data

The use of dense SNPs to predict the genetic value of an individual for a complex trait is often referred to as “genomic selection” in livestock and crops, but is also relevant to human genetics to predict, for example, complex genetic disease risk. The accuracy of prediction depends on the strength of linkage disequilibrium (LD) between SNPs and causal mutations. If sequence data were used instead of dense SNPs, accuracy should increase because causal mutations are present, but demographic history and long-term negative selection also influence accuracy. We therefore evaluated genomic prediction, using simulated sequence in two contrasting populations: one reducing from an ancestrally large effective population size (N_e) to a small one, with high LD common in domestic livestock, while the second had a large constant-sized N_e with low LD similar to that in some human or outbred plant populations. There were two scenarios in each population; causal variants were either neutral or under long-term negative selection. For large N_e, sequence data led to a 22% increase in accuracy relative to ∼600K SNP chip data with a Bayesian analysis and a more modest advantage with a BLUP analysis. This advantage increased when causal variants were influenced by negative selection, and accuracy persisted when 10 generations separated reference and validation populations. However, in the reducing N_e population, there was little advantage for sequence even with negative selection. This study demonstrates the joint influence of demography and selection on accuracy of prediction and improves our understanding of how best to exploit sequence for genomic prediction.     

Binding and Incision Activities of UvrABC Excinuclease on Slipped DNA Intermediates that Generate Frameshift Mutations

Previousin vivostudies involving sequence 5′-CCCG₁G₂G₃-3′ (SmaI site) have demonstrated that adducts ofN-2-acetylaminofluorene (AAF) to any of the three guanine residues of theSmaI sequence induce, with different efficiencies, two classes of −1 frameshift events, namely −G and −C mutations, referred to as targeted and semitargeted mutations, respectively. It has been proposed that both events occur during replication as a consequence of slippage events involving slipped mutagenic intermediates (SMIs). In order to evaluate the potential role of the UvrABC excinuclease in frameshift mutagenesis, we have studied the interaction of this enzyme with DNA molecules mimicking SMIsin vitro.In all of our constructions, when present, the AAF adduct was located on the third guanine residue of theSmaI site (5′-CCCG₁G₂G₃-3′). This strand was referred to as the top strand, the complementary strand being the bottom strand. Double-stranded heteroduplexes mimicking the targeted and semitargeted SMIs contained a deletion of a C and a G within theSmaI sequence in the bottom strand and were designated ΔC/3 and ΔG/3 when modified with the AAF on the third guanine residue in the top strand or ΔC/O and ΔG/O when unmodified. The modified homoduplex was designatedSmaI/3.ΔC/O and ΔG/O were weakly recognized by UvrA₂B, but not incised. All three AAF-modified substrates were recognized with similar efficiency and much more efficiently than unmodified heteroduplexes. With AAF-monomodified substrates, dissociation of UvrA₂from the UvrA₂B- DNA complex occurred more readily in heteroduplexes than in the homoduplex.SmaI/3 and ΔC/3 were incised with equal efficiency, while ΔG/3 was less incised. The position of the AAF lesion dictated the position of the incised phosphodiester bonds, suggesting that the presence of a bulge can modulate the yield but not the incision pattern of AAF-modified substrates. The finding that UvrABC excinuclease acts on substrates that mimic SMIs suggests that the nucleotide excision repair pathway may help in fixing frameshift mutations before the following round of replication.     

The effect of scale-free topology on the robustness and evolvability of genetic regulatory networks

We investigate how scale-free (SF) and Erd?s-Rényi (ER) topologies affect the interplay between evolvability and robustness of model gene regulatory networks with Boolean threshold dynamics. In agreement with Oikonomou and Cluzel (2006) we find that networks with SF_in topologies, that is SF topology for incoming nodes and ER topology for outgoing nodes, are significantly more evolvable towards specific oscillatory targets than networks with ER topology for both incoming and outgoing nodes. Similar results are found for networks with SF_both and SF_out topologies. The functionality of the SF_out topology, which most closely resembles the structure of biological gene networks (Babu et al., 2004), is compared to the ER topology in further detail through an extension to multiple target outputs, with either an oscillatory or a non-oscillatory nature. For multiple oscillatory targets of the same length, the differences between SF_out and ER networks are enhanced, but for non-oscillatory targets both types of networks show fairly similar evolvability. We find that SF networks generate oscillations much more easily than ER networks do, and this may explain why SF networks are more evolvable than ER networks are for oscillatory phenotypes. In spite of their greater evolvability, we find that networks with SF_out topologies are also more robust to mutations (mutational robustness) than ER networks. Furthermore, the SF_out topologies are more robust to changes in initial conditions (environmental robustness). For both topologies, we find that once a population of networks has reached the target state, further neutral evolution can lead to an increase in both the mutational robustness and the environmental robustness to changes in initial conditions.