Missense mutations and evolutionary conservation of amino acids: evidence that many of the amino acids in factor IX function as "spacer" elements.

We report 31 point mutations in the factor IX gene and explore the relationship between the level of evolutionary conservation of an amino acid and the probability of a mutation causing hemophilia B. From our total sample of 125 hemophiliacs and from those reported by others, we identify 95 independent missense mutations, 94 of which occur at amino acids that are evolutionarily conserved in the available mammalian factor IX sequences. The likelihood of a missense mutation causing hemophilia B depends on whether the residue is also conserved in the factor IX-related proteases: factor VII, factor X, and protein C. Most of the possible missense mutations in generically conserved residues (i.e., those conserved in factor IX and in all the related proteases) should cause disease. In contrast, missense mutations in factor IX-specific residues (i.e., those conserved in human, cow, dog, and mouse factor IX but not in the related proteases) are sixfold less likely to cause disease. Missense mutations at nonconserved residues are 33-fold less likely to cause disease. At least three models are compatible with these observations. A comparison of sequence alignments from four and nine species of factor IX and an examination of the missense mutations occurring at CpG residues suggest a model in which most residues fall on opposite ends of a spectrum. In about 40% of residues, virtually any missense mutation in a minority of the residues will cause disease, while virtually no missense mutations will cause disease in most of the remaining residues. Thus, many of the residues in factor IX are spacers; that is, the main chains are presumably necessary to keep other amino acid interactions in register, but the nature of the side chain is unimportant.     

Mutational and bioinformatics analysis of proline- and glycine-rich motifs in vesicular acetylcholine transporter

The vesicular acetylcholine transporter (VAChT) contains six conserved sequence motifs that are rich in proline and glycine. Because these residues can have special roles in the conformation of polypeptide backbone, the motifs might have special roles in conformational changes during transport. Using published bioinformatics insights, the amino acid sequences of the 12 putative, helical, transmembrane segments of wild-type and mutant VAChTs were analyzed for propensity to form non-alpha-helical conformations and molecular notches. Many instances were found. In particular, high propensity for kinks and notches are robustly predicted for motifs D2, C and C'. Mutations in these motifs either increase or decrease Vmax for transport, but they rarely affect the equilibrium dissociation constants for ACh and the allosteric inhibitor, vesamicol. The near absence of equilibrium effects implies that the mutations do not alter the backbone conformation. In contrast, the Vmax effects demonstrate that the mutations alter the difficulty of a major conformational change in transport. Interestingly, mutation of an alanine to a glycine residue in motif C significantly increases the rates for reorientation across the membrane. These latter rates are deduced from the kinetics model of the transport cycle. This mutation is also predicted to produce a more flexible kink and tighter tandem notches than are present in wild-type. For the full set of mutations, faster reorientation rates correlate with greater predicted propensity for kinks and notches. The results of the study argue that conserved motifs mediate conformational changes in the VAChT backbone during transport.     

Molecular characterization of ura1, a mutant allele for orotidine-5'-monophosphate decarboxylase in Schizophyllum commune

Abstract The basis of the auxotrophic ural phenotype in Schizophyllum commune has been investigated. Two point mutations causing changes in conserved amino acid positions 62 (from lysine to glutamate) and 79 (from leucine to phenylalanine) most likely are the cause for the observed phenotype, whereas the overall gene structure was unchanged. Since reversion rates in this locus are extremely low, a single point mutation could not be expected to be the cause for the mutation. Besides the two point mutations expected to be induced by UV mutagenesis, the two alleles investigated from independently isolated strains differ by approximately 7% in nucleic acid sequence and about 3% in amino acid sequence, indicating a distant relationship between the strains used.     

Protein-DNA interactions: amino acid conservation and the effects of mutations on binding specificity

We investigate the conservation of amino acid residue sequences in 21 DNA-binding protein families and study the effects that mutations have on DNA-sequence recognition. The observations are best understood by assigning each protein family to one of three classes: (i) non-specific, where binding is independent of DNA sequence; (ii) highly specific, where binding is specific and all members of the family target the same DNA sequence; and (iii) multi-specific, where binding is also specific, but individual family members target different DNA sequences. Overall, protein residues in contact with the DNA are better conserved than the rest of the protein surface, but there is a complex underlying trend of conservation for individual residue positions. Amino acid residues that interact with the DNA backbone are well conserved across all protein families and provide a core of stabilising contacts for homologous protein-DNA complexes. In contrast, amino acid residues that interact with DNA bases have variable levels of conservation depending on the family classification. In non-specific families, base-contacting residues are well conserved and interactions are always found in the minor groove where there is little discrimination between base types. In highly specific families, base-contacting residues are highly conserved and allow member proteins to recognise the same target sequence. In multi-specific families, base-contacting residues undergo frequent mutations and enable different proteins to recognise distinct target sequences. Finally, we report that interactions with bases in the target sequence often follow (though not always) a universal code of amino acid-base recognition and the effects of amino acid mutations can be most easily understood for these interactions.     

Identification of biochemically neutral positions in liver pyruvate kinase

Understanding how each residue position contributes to protein function has been a long-standing goal in protein science. Substitution studies have historically focused on conserved protein positions. However, substitutions of nonconserved positions can also modify function. Indeed, we recently identified nonconserved positions that have large substitution effects in human liver pyruvate kinase (hLPYK), including altered allosteric coupling. To facilitate a comparison of which characteristics determine when a nonconserved position does vs does not contribute to function, the goal of the current work was to identify neutral positions in hLPYK. However, existing hLPYK data showed that three features commonly associated with neutral positions—high sequence entropy, high surface exposure, and alanine scanning—lacked the sensitivity needed to guide experimental studies. We used multiple evolutionary patterns identified in a sequence alignment of the PYK family to identify which positions were least patterned, reasoning that these were most likely to be neutral. Nine positions were tested with a total of 117 amino acid substitutions. Although exploring all potential functions is not feasible for any protein, five parameters associated with substrate/effector affinities and allosteric coupling were measured for hLPYK variants. For each position, the aggregate functional outcomes of all variants were used to quantify a “neutrality” score. Three positions showed perfect neutral scores for all five parameters. Furthermore, the nine positions showed larger neutral scores than 17 positions located near allosteric binding sites. Thus, our strategy successfully enriched the dataset for positions with neutral and modest substitutions.     

Interspecies relationships among ADP-ribosylation factors (ARFs): Evidence of evolutionary pressure to maintain individual identities

ADP-ribosylation factors (ARFs) are ～20-kDa guanine nucleotide-binding proteins that are allosteric activators of the NAD:arginine ADP-ribosyltransferase activity of cholera toxin and appear to play a role in intracellular vesicular trafficking. Although the physiological roles of these proteins have not been defined, it has been presumed that each has a specific intracellular function. To obtain genetic evidence that each ARF is under evolutionary pressure to maintain its structure, and presumably function, rat ARF cDNA clones were isolated and their nucleotide and deduced amino acid sequences were compared to those of other mammalian ARFs. Deduced amino acid sequences for rat ARFs 1, 2, 3, 5 and 6 were identical to those of the known cognate human and bovine ARFs; rat ARF4 was 96% identical to human ARF4. Nucleotide sequences of both the untranslated as well as the coding regions were highly conserved. These results indicate that the ARF proteins are, as a family, extraordinarily well conserved across mammalian species. The unusually high degree of conservation of the untranslated regions is consistent with these regions having important regulatory roles and that individual ARFs contain structurally unique elements required for specific functions.     

Characterisation of XlCdc1, a Xenopus homologue of the small (PolD2) subunit of DNA polymerase delta; identification of ten conserved regions I-X based on protein sequence comparisons across ten eukaryotic species

DNA polymerase delta (Pol delta), which plays keys roles in DNA replication, repair and recombination in eukaryotic cells, comprises at least two essential subunits - a large catalytic subunit (PolD1) possessing both DNA polymerase and 3'-5' exonuclease activities, and a smaller subunit (PolD2) whose function is not yet clear. Here we describe the cloning and sequencing of a Xenopus cDNA encoding a homologue of the PolD2 subunit. This protein (designated XlCdc1) is 69% identical to the human PolD2 protein and 34% identical to fission yeast Cdc1. Alignment of PolD2 protein sequences across ten eukaryotic species identifies 36 invariant amino-acid positions. These 36 residues are located within ten conserved regions (designated I-X) likely to have key functional roles. Consistent with this, the mutations in six previously identified yeast mutant PolD2 proteins map within conserved regions III, VI, VII and VIII. Several of the invariant amino acids are also conserved across the archaeal DNA polymerase II DP1 protein family.     

Replacement of C305 in Heart/Muscle-Type Isozyme of Human Carnitine Palmitoyltransferase I with Aspartic Acid and Other Amino Acids

Liver- and heart/muscle-type isozymes of human carnitine palmitoyltransferase I (L- and M-CPTI, respectively) show a certain similarity in their amino acid sequences, and mutation studies on the conserved amino acids between these two isozymes often show essentially the same effects on their enzymatic properties. Earlier mutation studies on C305 in human M-CPTI and its counterpart residue, C304, in human L-CPTI showed distinct effects of the mutations, especially in the aspect of enzyme stability; however, simple comparison of these effects on the conserved Cys residue between L- and M-CPTI was difficult, because these studies were carried out using different expression systems and distinct amino acids as replacements. In the present study, we carried out mutation studies on the C305 in human M-CPTI using COS cells for the expression system. Our results showed that C305 was replaceable with aspartic acid but that substitution with other amino acids caused both loss of function and reduced expression.     

Correlated Mutations: A Hallmark of Phenotypic Amino Acid Substitutions

Point mutations resulting in the substitution of a single amino acid can cause severe functional consequences, but can also be completely harmless. Understanding what determines the phenotypical impact is important both for planning targeted mutation experiments in the laboratory and for analyzing naturally occurring mutations found in patients. Common wisdom suggests using the extent of evolutionary conservation of a residue or a sequence motif as an indicator of its functional importance and thus vulnerability in case of mutation. In this work, we put forward the hypothesis that in addition to conservation, co-evolution of residues in a protein influences the likelihood of a residue to be functionally important and thus associated with disease. While the basic idea of a relation between co-evolution and functional sites has been explored before, we have conducted the first systematic and comprehensive analysis of point mutations causing disease in humans with respect to correlated mutations. We included 14,211 distinct positions with known disease-causing point mutations in 1,153 human proteins in our analysis. Our data show that (1) correlated positions are significantly more likely to be disease-associated than expected by chance, and that (2) this signal cannot be explained by conservation patterns of individual sequence positions. Although correlated residues have primarily been used to predict contact sites, our data are in agreement with previous observations that (3) many such correlations do not relate to physical contacts between amino acid residues. Access to our analysis results are provided at http://webclu.bio.wzw.tum.de/~pagel/supplements/correlated-positions/.     

Human germline mutation in the factor IX gene

The molecular epidemiology of factor IX germline mutations in patients with hemophilia B has been studied in detail because it is an advantageous model for analyzing recent germline mutations in humans. It is estimated that mutations have been defined in the majority of nucleotides that are the target for mutation. The likelihood that a factor IX missense mutation will cause disease correlates with the degree of evolutionary conservation of the amino acid. Mutation rates per base-pair have been estimated after careful consideration and correction for biases, predicting about 76 de novo mutations per generation per individual resulting in 0.3 deleterious changes. The male-to-female sex ratio of mutation varies with the type of mutation. There is evidence for a maternal age effect and an excess of non-CpG G:C to A:T transitions. The factor IX mutation pattern is similar among geographically, racially and ethnically diverse human populations. The data support primarily endogenous mechanisms of germline mutation in the factor IX gene. Mutations at splice junctions are compatible with simple rules for predicting disease causing mutations.     

High conservation of sequences involved in cystic fibrosis mutations in five mammalian species

Several mutations have been identified in the first nucleocide binding fold (NBF) of the cystic fibrosis (CF) transmembrane conductance regulator (CFTR) gene. We have analyzed the DNA sequences of exons 10 and 11 in five different mammalian species, marmoset, mouse, cow, pig, and sheep; the amino acid conservation studied for nine disease mutations; and two “benign” mutations. For exon 10,87% homology at the DNA level and 93.5% at the amino acid level were found for these species. For exon 11, the lowest homology (70%), as found in mouse and the highest in marmoset (93%), whereas the amino acid sequence conservation ranged from 82.5 to 100%. All codons involved in CF mutations are highly conserved throughout evolution.     

Dissection of the effector-binding site and complementation studies of Escherichia coli phosphofructokinase using site-directed mutagenesis

A systematic study by site-directed mutagenesis has been conducted on the effector site of phosphofructokinase from Escherichia coli to delineate the role of side chains in binding the allosteric activator, GDP, and inhibitor, PEP, and to search for key residues in the allosteric transtion. Target residues were identified from the crystal structure of the enzyme-nucleoside diphosphate complex. It is found that both activator and inhibitor bind to the same set of amino acid side chains. Deletion of positively charged groups (Arg21, Arg25, Arg54, Arg154, and Lys213 mutated to alanine) weakens binding of both effectors by 2-3 kcal/mol, consistent with the disruption of charged hydrogen bonds. Residue Glu187, which is known from the crystal structure to bind the coordinated Mg2+ ion of GDP, is found to have a unique behavior on mutation and appears to be crucial in triggering the allosteric transition. All other residues mutated simply weaken binding of both PEP and GDP in a parallel manner. However, mutation of Glu----Ala187 reverses the roles of GDP and PEP, causing GDP to become an allosteric inhibitor and PEP an activator. Mutation of Glu----Gln187 has only a small effect on the binding of PEP, and both PEP and GDP are inhibitors. Studies are described in which mutations in different subunits of a tetrameric complex complement each other. The effector site is composed of residues from two subunits. In particular, Arg21 and Lys213 in each site are from different subunits. Mutations of either one of these residues abolishes activation by GDP of the homotetramer.(ABSTRACT TRUNCATED AT 250 WORDS)     

Mutations in ribosomal protein L3 and 23S ribosomal RNA at the peptidyl transferase centre are associated with reduced susceptibility to tiamulin in Brachyspira spp. isolates

The pleuromutilin antibiotic tiamulin binds to the ribosomal peptidyl transferase centre. Three groups of Brachyspira spp. isolates with reduced tiamulin susceptibility were analysed to define resistance mechanisms to the drug. Mutations were identified in genes encoding ribosomal protein L3 and 23S rRNA at positions proximal to the peptidyl transferase centre. In two groups of laboratory-selected mutants, mutations were found at nucleotide positions 2032, 2055, 2447, 2499, 2504 and 2572 of 23S rRNA (Escherichia coli numbering) and at amino acid positions 148 and 149 of ribosomal protein L3 (Brachyspira pilosicoli numbering). In a third group of clinical B. hyodysenteriae isolates, only a single mutation at amino acid 148 of ribosomal protein L3 was detected. Chemical footprinting experiments show a reduced binding of tiamulin to ribosomal subunits from mutants with decreased susceptibility to the drug. This reduction in drug binding is likely the resistance mechanism for these strains. Hence, the identified mutations located near the tiamulin binding site are predicted to be responsible for the resistance phenotype. The positions of the mutated residues relative to the bound drug advocate a model where the mutations affect tiamulin binding indirectly through perturbation of nucleotide U2504.     

Amino acid conservation and clinical severity of human glucose-6-phosphate dehydrogenase mutations

More than a hundred naturally occurring mutations of human glucose-6-phosphate dehydrogenase (G6PD) have been identified at the amino acid level. The abundance of distinct mutation sites and their clinical manifestations make this enzyme ideal for structure-function analysis studies. We present here a sequence and structure combined analysis by which the severity of clinical symptoms resulting from point mutations of this enzyme is correlated with quantified degrees of amino acid conservation within 23 G6PD sequences from different organisms. Our analysis verifies, on a quantitative basis, a widely held notion that clinically severer mutations of G6PD usually occur at conserved amino acids. However, marked exceptions to this general trend exist which are most notably revealed by a number of mutations associated with chronic nonspherocytic hemolytic anemia (class I variants). When mapped onto a homology-derived structural model of human G6PD, these class I mutational sites of low amino acid conservation appear to localize in two spatially distinct clusters, both of which are populated with mutations consisting mainly of clinically severer variants (i.e. class I and class II). These results of computer-assisted analyses contribute to a further understanding of the structure-function relationships of human G6PD deficiency.     

Paired natural cysteine mutation mapping: aid to constraining models of protein tertiary structure.

This paper discusses the benefit of mapping paired cysteine mutation patterns as a guide to identifying the positions of protein disulfide bonds. This information can facilitate the computer modeling of protein tertiary structure. First, a simple, paired natural-cysteine-mutation map is presented that identifies the positions of putative disulfide bonds in protein families. The method is based on the observation that if, during the process of evolution, a disulfide-bonded cysteine residue is not conserved, then it is likely that its counterpart will also be mutated. For each target protein, protein databases were searched for the primary amino acid sequences of all known members of distinct protein families. Primary sequence alignment was carried out using PileUp algorithms in the GCG package. To search for correlated mutations, we listed only the positions where cysteine residues were highly conserved and emphasized the mutated residues. In proteins of known three-dimensional structure, a striking pattern of paired cysteine mutations correlated with the positions of known disulfide bridges. For proteins of unknown architecture, the mutation maps showed several positions where disulfide bridging might occur.     

Comparison of class A and D G protein-coupled receptors: common features in structure and activation

All G protein-coupled receptors (GPCRs) share a common seven TM helix architecture and the ability to activate heterotrimeric G proteins. Nevertheless, these receptors have widely divergent sequences with no significant homology. We present a detailed structure-function comparison of the very divergent Class A and D receptors to address whether there is a common activation mechanism across the GPCR superfamily. The Class A and D receptors are represented by the vertebrate visual pigment rhodopsin and the yeast alpha-factor pheromone receptor Ste2, respectively. Conserved amino acids within each specific receptor class and amino acids where mutation alters receptor function were located in the structures of rhodopsin and Ste2 to assess whether there are functionally equivalent positions or regions within these receptors. We find several general similarities that are quite striking. First, strongly polar amino acids mediate helix interactions. Their mutation generally leads to loss of function or constitutive activity. Second, small and weakly polar amino acids facilitate tight helix packing. Third, proline is essential at similar positions in transmembrane helices 6 and 7 of both receptors. Mapping the specific location of the conserved amino acids and sites of constitutively active mutations identified conserved microdomains on transmembrane helices H3, H6, and H7, suggesting that there are underlying similarities in the mechanism of the widely divergent Class A and Class D receptors.     

Large scale mtDNA sequencing reveals sequence and functional conservation as major determinants of homoplasmic mtDNA variant distribution

The mitochondrial DNA (mtDNA) is highly variable, containing large numbers of pathogenic mutations and neutral polymorphisms. The spectrum of homoplasmic mtDNA variation was characterized in 730 subjects and compared with known pathogenic sites. The frequency and distribution of variants in protein coding genes were inversely correlated with conservation at the amino acid level. Analysis of tRNA secondary structures indicated a preference of variants for the loops and some acceptor stem positions. This comprehensive overview of mtDNA variants distinguishes between regions and positions which are likely not critical, mainly conserved regions with pathogenic mutations and essential regions containing no mutations at all.     

Genetic analysis of the human immunodeficiency virus type 1 integrase protein.

Single-amino-acid changes in a highly conserved central region of the human immunodeficiency virus type 1 (HIV-1) integrase protein were analyzed for their effects on viral protein synthesis, virion morphogenesis, and viral replication. Alteration of two amino acids that are invariant among retroviral integrases, D116 and E152 of HIV-1, as well as a mutation of the highly conserved amino acid S147 blocked viral replication in two CD4+ human T-cell lines. Mutations of four other highly conserved amino acids in the region had no detectable effect on viral replication, whereas mutations at two positions, N117 and Y143, resulted in viruses with a delayed-replication phenotype. Defects in virion precursor polypeptide processing, virion morphology, or viral DNA synthesis were observed for all of the replication-defective mutants, indicating that changes in integrase can have pleiotropic effects on viral replication.