Relation between protein stability, evolution and structure, as probed by carboxylic acid mutations

Native proteins are marginally stable. Low thermodynamic stability may actually be advantageous, although the accumulation of neutral, destabilizing mutations may have also contributed to it. In any case, once marginal stability has been reached, it appears plausible that mutations at non-constrained positions become fixed in the course of evolution (due to random drift) with frequencies that roughly reflect the mutation effects on stability ("pseudo-equilibrium hypothesis"). We have found that all glutamate-->aspartate mutations in wild-type Escherichia coli thioredoxin are destabilizing, as well as most of the aspartate-->glutamate mutations. Furthermore, the effect of these mutations on thioredoxin thermodynamic stability shows a robust correlation with the frequencies of occurrence of the involved residues in several-hundred sequence alignments derived from a BLAST search. These results provide direct and quantitative experimental evidence for the pseudo-equilibrium hypothesis and should have general consequences for the interpretation of mutation effects on protein stability, as they suggest that residue environments in proteins may be optimized for stabilizing interactions to a remarkable degree of specificity. We also provide evidence that such stabilizing interactions may be detected in sequence alignments, and briefly discuss the implications of this possibility for the derivation of structural information (on native and denatured states) from comparative sequence analyses.     

Genic mutation rates in mammals: local similarity,chromosomal heterogeneity,and X-versus-autosome disparity

The reduction of mutation rates on the mammalian X chromosome relative to autosomes is most often explained in the literature as evidence of male-driven evolution. This hypothesis attributes lowered mutation rates on the X chromosome to the fact that this chromosome spends less time in the germline of males than in the germline of females. In contrast to this majority view, two articles argued that the patterns of mutation rates across chromosomes are inconsistent with male-driven evolution. One article reported a 40% reduction in synonymous substitution rates (Ks) for X-linked genes relative to autosomes in the mouse-rat lineage. The authors argued that this reduction is too dramatic to be explained by male-driven evolution and concluded that selection has systematically reduced mutation rate on the X chromosome to a level optimal for this male-hemizygous chromosome. More recently, a second article found that chromosomal mutation rates in both the human-mouse and mouse-rat lineages were so heterogeneous that the X chromosome was not an outlier. Here again, the authors argued that this is at odds with male-driven evolution and suggested that selection has modulated chromosomal mutation rates to locally optimal levels, thus extending the argument of the first mentioned article to include autosomes. Here, we reexamine these conclusions using mouse-rat and human-mouse coding-region data. We find a more modest reduction of Ks on the X chromosome, but our results contradict the finding that the X chromosome is not distinct from autosomes. Multiple statistical tests show that Ks rates on the X chromosome differ systematically from the autosomes in both lineages. We conclude that the moderate reduction of mutation rate on the X chromosome of both lineages is consistent with male-driven evolution; however, the large variance in mutation rates across chromosomes suggests that mutation rates are affected by additional factors besides male-driven evolution. Investigation of mutation rates by synteny reveals that synteny blocks, rather than entire chromosomes, might represent the unit of mutation rate variation.     

Molecular mechanisms of functional rescue mediated by P53 tumor suppressor mutations

We have utilized both molecular dynamics simulations and solution biophysical measurements to investigate the rescue mechanism of mutation N235K, which plays a key role in the recently identified global suppressor motif of K235/Y239/R240 in the human p53 DNA-binding domain (DBD). Previous genetic analysis indicates that N235K alone rescues five out of six destabilized cancer mutants. However, the solution biophysical measurement shows that N235K generates only a slight increase to the stability of DBD, implying a rescue mechanism that is not a simple additive contribution to thermodynamic stability. Our molecular simulations show that the N235K substitution generates two non-native salt bridges with residues D186 and E198. We find that the nonnative salt bridges, D186-K235 and E198-K235, and a native salt bridge, E171-R249, are mutually exclusive, thus resulting in only a marginal increase in stability as compared to the wild type protein. When a destabilized V157F is paired with N235K, the native salt bridge E171-R249 is retained. In this context, the non-native salt bridges, D186-K235 and E198-K235, produce a net increase in stability as compared to V157F alone. A similar rescue mechanism may explain how N235K stabilize other highly unstable β-sandwich cancer mutants.     

Accepted mutations in a gene family: Evolutionary diversification of duplicated DNA

Summary We report and compare the DNA sequences of 14 silkmoth (Antheraea polyphemus) chorion genes, derived from either cDNA or chromosomal DNA clones. Seven of these genes are members of the A multigene family, and seven are members of the B family. Where available, the previously reported (Jones and Kafatos 1980) intronic and extragenic flanking DNA sequences are also considered. Closely related sequences are compared, revealing the types of spontaneous mutations that were fixed during paralogous evolution. Segmental mutations (i.e. mutations other than substitutions) are nearly always interpretable as small duplications or deletions. related to small direct repeats. Segmental mutations are strongly constrained in the coding regions, although they do occur. Nucleotide substitutions also appear to be under selective constraints: relatively few substitutions leading to amino acid replacements are accepted, silent substitutions leading to some codons (especially purine-terminated ones) are disfavored, and different compositional biases are maintained in different parts of the sequences. Other sequence differences can be interpreted as indicative of neutral drift, including most differences in non-coding regions and most T/C transitions in third-base positions. In the non-coding regions, which are thought to be only loosely constrained by selection, transitions are observed more frequently than might be expected: they account for 52% of all substitutions, and they appear to be favored two to threefold over transversions when allowance is made for the skewed base composition of these regions.     

Naturally occurring mutations at the acetylcholine receptor binding site independently alter ACh binding and channel gating

By defining functional defects in a congenital myasthenic syndrome (CMS), we show that two mutant residues, located in a binding site region of the acetylcholine receptor (AChR) epsilon subunit, exert opposite effects on ACh binding and suppress channel gating. Single channel kinetic analysis reveals that the first mutation, epsilon N182Y, increases ACh affinity for receptors in the resting closed state, which promotes sequential occupancy of the binding sites and discloses rate constants for ACh occupancy of the nonmutant alphadelta site. Studies of the analogous mutation in the delta subunit, deltaN187Y, disclose rate constants for ACh occupancy of the nonmutant alpha epsilon site. The second CMS mutation, epsilon D175N, reduces ACh affinity for receptors in the resting closed state; occupancy of the mutant site still promotes gating because a large difference in affinity is maintained between closed and open states. epsilon D175N impairs overall gating, however, through an effect independent of ACh occupancy. When mapped on a structural model of the AChR binding site, epsilon N182Y localizes to the interface with the alpha subunit, and epsilon D175 to the entrance of the ACh binding cavity. Both epsilon N182Y and epsilon D175 show state specificity in affecting closed relative to desensitized state affinities, suggesting that the protein chain harboring epsilon N182 and epsilon D175 rearranges in the course of receptor desensitization. The overall results show that key residues at the ACh binding site differentially stabilize the agonist bound to closed, open and desensitized states, and provide a set point for gating of the channel.     

The distribution of fitness effects of new beneficial mutations in Pseudomonas fluorescens

Theoretical studies of adaptation emphasize the importance of understanding the distribution of fitness effects (DFE) of new mutations. We report the isolation of 100 adaptive mutants—without the biasing influence of natural selection—from an ancestral genotype whose fitness in the niche occupied by the derived type is extremely low. The fitness of each derived genotype was determined relative to a single reference type and the fitness effects found to conform to a normal distribution. When fitness was measured in a different environment, the rank order changed, but not the shape of the distribution. We argue that, even with detailed knowledge of the genetic architecture underpinning the adaptive types (as is the case here), the DFEs remain unpredictable, and we discuss the possibility that general explanations for the shape of the DFE might not be possible in the absence of organism-specific biological details.     

Common conformational effects of p53 mutations

The tumor suppressor gene p53 has been identified as the most frequent target of genetic alterations in human cancers. Most of these mutations occur in highly conserved regions in the DNA-binding core domain of the p53 protein, suggesting that the amino acid residues in these regions are critical for maintaining normal p53 structure and function. We previously used molecular dynamics calculations to demonstrate that several amino acid substitutions in these regions that are induced by environmental carcinogens and found in human tumors produce certain common conformational changes in the mutant proteins that differ substantially from the wild-type structure. In order to determine whether these conformational changes are consistent for other p53 mutants, we have now used molecular dynamics to determine the structure of the DNA-binding core domain of seven other environmentally induced, cancer-related p53 mutants, namely His 175, Asp 245, Asn 245, Trp 248, Met 249, Ser 278, and Lys 286. The results indicate that all of these mutants differ substantially from the wild-type structure in certain discrete regions and that some of these conformational changes are similar for these mutants as well as those determined previously. The changes are also consistent with experimental evidence for alterations in structure in p53 mutants determined by epitope detectability using monoclonal antibodies directed against these regions of predicted conformational change.     

Stepwise manipulation of DNA specificity in Flp recombinase: progressively adapting Flp to individual and combinatorial mutations in its target site

The Flp protein from Saccharomyces cerevisiae is one of the site-specific tyrosine family recombinases that are used widely in genomic engineering. As a first step towards mediating directed DNA rearrangements at non-native Flp recombination targets (mFRTs), we have evolved three separate groups of Flp variants that preferentially act on mFRTs containing substitutions at the first, seventh or both positions of the Flp-binding elements. The variants that recombine the double-mutant mFRT contain a subset of the mutations present in those that are active on the single-mutant mFRTs, plus additional mutations. Specificity for and discrimination between target sites, effected primarily by amino acid residues that contact DNA, can be modulated by those that do not interact with DNA or with a DNA-contacting residue. The degree of modulation can range from relaxed DNA specificity to almost completely altered specificity. Our results suggest that combined DNA shuffling and mutagenesis of libraries of Flp variants active on distinct mFRTs can yield variants that can recombine mFRTs containing combinations of the individual mutations.     

Functional significance of conserved amino acid residues.

A systemic study of single amino acid substitutions in bacteriophage T4 lysozyme permitted a test of the concept that conserved amino acid residues are more functionally important than nonconserved residues. Substitutions of amino acid residues that are conserved among five bacteriophage-encoded lysozymes were found to lead more frequently to loss of function than substitutions of nonconserved residues. Of 163 residues tested, only 74 (45%) are sensitive to at least one substitution; however, all 14 residues that are fully conserved are sensitive to substitutions.     

Genetic analysis of mutations that alter cell fates in maize leaves: Dominant Liguleless mutations

The three major components of the maize leaf are the blade, the sheath, and at their junction, the ligular region. Each exhibits specific cell types and organization. Four dominant Liguleless (Lg) mutations (Lg3-O, Lg4-O, Lg*347, and Lg*9167) in at least three different genes cause a similar morphological phenotype in leaves, although each mutation affects a distinct domain of the blade. Mutant leaves display regions of altered cell fate in the blade, occompanied by elimination of ligule and auricle at their wild-type positions and development of ligule and auricle in the blade at the borders of the altered regions. The affected blade cells are transformed into sheath-like cells, as determined by morphological and genetic tests. Lg4-O expressivity is highly dependent on genetic background. For example, two different backgrounds may specify converse patterns of phenotypic expression. Lg4-O expressivity is also affected by the heterochronic mutation Teopod2 (Tp2). Gene dosage experiments indicate that Lg4-O is a neomorph. Interactions between recessive lg mutations (which eliminate ligular structures) and the dominant Lg mutations suggest that the lg⁺ genes act after the Lg mutations. Lg3-O and Lg4-O act semidominantly, and interact with each other and with other mutations in the Knotted1 (Kn1)-like family (a family in which dominant mutant alleles cause blade to sheath transformation phenotypes). These interactions suggest that the above Kn1-like mutations may function similarly in the leaf. We discuss the similarities between the Lg mutations and the other mutations of the Kn1-like family, which led us to postulate that lg3 and lg4 are members of a growing family of kn1-like (knox) homeobox genes that are identified by dominant mutant alleles causing leaf transformation phenotypes. We also propose that certain key characteristics of this family of dominant neomorphic mutations are important for generating meaningful morphological changes during evolution. © 1996 Wiley-Liss, Inc.     

Effects of frontotemporal dementia FTDP-17 mutations on heparin-induced assembly of tau filaments

Missense mutations and intronic mutations in the gene for microtubule-associated protein tau cause frontotemporal dementia and Parkinsonism linked to chromosome 17 (FTDP-17). Most missense mutations have as likely primary effect a reduced ability of tau to interact with microtubules. We report here an additional effect of several missense mutations, namely the stimulation of heparin-induced filament assembly of recombinant tau, despite the absence of any change in structure indicated by circular dichroism. These findings indicate that missense mutations in tau lead to frontotemporal dementia through potentially multiple mechanisms.     

A three‐dimensional model of mammalian tyrosinase active site accounting for loss of function mutations

Tyrosinases are the first and rate‐limiting enzymes in the synthesis of melanin pigments responsible for colouring hair, skin and eyes. Mutation of tyrosinases often decreases melanin production resulting in albinism, but the effects are not always understood at the molecular level. Homology modelling of mouse tyrosinase based on recently published crystal structures of non‐mammalian tyrosinases provides an active site model accounting for loss‐of‐function mutations. According to the model, the copper‐binding histidines are located in a helix bundle comprising four densely packed helices. A loop containing residues M374, S375 and V377 connects the CuA and CuB centres, with the peptide oxygens of M374 and V377 serving as hydrogen acceptors for the NH‐groups of the imidazole rings of the copper‐binding His367 and His180. Therefore, this loop is essential for the stability of the active site architecture. A double substitution ³⁷⁴MS³⁷⁵→³⁷⁴GG³⁷⁵ or a single M374G mutation lead to a local perturbation of the protein matrix at the active site affecting the orientation of the H367 side chain, that may be unable to bind CuB reliably, resulting in loss of activity. The model also accounts for loss of function in two naturally occurring albino mutations, S380P and V393F. The hydroxyl group in S380 contributes to the correct orientation of M374, and the substitution of V393 for a bulkier phenylalanine sterically impedes correct side chain packing at the active site. Therefore, our model explains the mechanistic necessity for conservation of not only active site histidines but also adjacent amino acids in tyrosinase.     

Effects of active site mutations in haemoglobin I from Lucina pectinata: a molecular dynamic study

Haemoglobin I from Lucina pectinata is a monomeric protein consisting of 142 amino acids. Its active site contains a peculiar arrangement of phenylalanine residues (PheB10, PheCD1 and PheE11) and a distal Gln at position E7. Active site mutations at positions B10, E7 and E11 were performed in deoxy haemoglobin I (HbI), followed by 10 ns molecular dynamic simulations. The results showed that the mutations induced changes in domains far from the active site producing more flexible structures than the native HbI. Distance analyses revealed that the heme pocket amino acids at positions E7 and B10 are extremely sensitive to any heme pocket residue mutation. The high flexibility observed by the E7 position suggests an important role in the ligand binding kinetics in ferrous HbI, while both positions play a major role in the ligand stabilisation processes. Furthermore, our results showed that E11Phe plays a pivotal role in protein stability.     

Towards the identification of the allosteric Phe-binding site in phenylalanine hydroxylase

The enzyme phenylalanine hydroxylase (PAH) is defective in the inherited disorder phenylketonuria. PAH, a tetrameric enzyme, is highly regulated and displays positive cooperativity for its substrate, Phe. Whether Phe binds to an allosteric site is a matter of debate, despite several studies worldwide. To address this issue, we generated a dimeric model for Phe–PAH interactions, by performing molecular docking combined with molecular dynamics simulations on human and rat wild-type sequences and also on a human G46S mutant. Our results suggest that the allosteric Phe-binding site lies at the dimeric interface between the regulatory and the catalytic domains of two adjacent subunits. The structural and dynamical features of the site were characterized in depth and described. Interestingly, our findings provide evidence for lower allosteric Phe-binding ability of the G46S mutant than the human wild-type enzyme. This also explains the disease-causing nature of this mutant.     

Positive natural selection has driven the evolution of the Hsp70s in Diguetia spiders

Hsp70s are a ubiquitous family of highly conserved proteins. Hsp70s are chaperones and have important roles in both protein folding and thermotolerance. It has been widely assumed that Hsp70 sequence evolution is governed by the strong functional constraints imposed by its crucial cellular functions. In this study of cytosolic heat-inducible Hsp70s from three spider families, we have found clear evidence of positive natural selection altering Hsp70s in desert-dwelling and heat-loving Diguetidae spiders. These spiders are a small family restricted to deserts. They display heat-tolerant behaviours not seen in their closest relatives, the Pholcidae and Plectreuridae.     

Scaling of Folding Properties in Go Models of Proteins

Insights about scaling of folding properties of proteins are obtained bystudying folding in heteropolymers described by Go-like Hamiltonians. Bothlattice and continuum space models are considered. In the latter case, themonomer-monomer interactions correspond to the Lennard-Jones potential.Several statistical ensembles of the two- and three-dimensional targetnative conformations are considered. Among them are maximally compactconformations which are confined to a lattice and those which are obtainedeither through quenching or annealing of homopolymers to their compactlocal energy minima. Characteristic folding times are found to grow aspower laws with the system size. The corresponding exponents are notuniversal. The size related deterioration of foldability is found to beconsistent with the scaling behavior of the characteristic temperatures:asymptotically, the folding temperature becomes much lower than thetemperature at which glassy kinetics become important. The helicalconformations are found to have the lowest overall scaling exponent andthe best foldability among the classes of conformations studied. Thescaling properties of the Go-like models of the protein conformationsstored in the Protein Data Bank suggest that proteins are not optimizedkinetically.