Mutation patterns of amino acid tandem repeats in the human proteome

Background: Amino acid tandem repeats are found in nearly one-fifth of human proteins. Abnormal expansion of these regions is associated with several human disorders. To gain further insight into the mutational mechanisms that operate in this type of sequence, we have analyzed a large number of mutation variants derived from human expressed sequence tags (ESTs).Results: We identified 137 polymorphic variants in 115 different amino acid tandem repeats. Of these, 77 contained amino acid substitutions and 60 contained gaps (expansions or contractions of the repeat unit). The analysis showed that at least about 21% of the repeats might be polymorphic in humans. We compared the mutations found in different types of amino acid repeats and in adjacent regions. Overall, repeats showed a five-fold increase in the number of gap mutations compared to adjacent regions, reflecting the action of slippage within the repetitive structures. Gap and substitution mutations were very differently distributed between different amino acid repeat types. Among repeats containing gap variants we identified several disease and candidate disease genes.Conclusion: This is the first report at a genome-wide scale of the types of mutations occurring in the amino acid repeat component of the human proteome. We show that the mutational dynamics of different amino acid repeat types are very diverse. We provide a list of loci with highly variable repeat structures, some of which may be potentially involved in disease.

Ultra-Deep Pyrosequencing Detects Conserved Genomic Sites and Quantifies Linkage of Drug-Resistant Amino Acid Changes in the Hepatitis B Virus Genome

Selection of amino acid substitutions associated with resistance to nucleos(t)ide-analog (NA) therapy in the hepatitis B virus (HBV) reverse transcriptase (RT) and their combination in a single viral genome complicates treatment of chronic HBV infection and may affect the overlapping surface coding region. In this study, the variability of an overlapping polymerase-surface region, critical for NA resistance, is investigated before treatment and under antiviral therapy, with assessment of NA-resistant amino acid changes simultaneously occurring in the same genome (linkage analysis) and their influence on the surface coding region.

A hierarchical Bayesian model for predicting the functional consequences of amino-acid polymorphisms

Genetic polymorphisms in deoxyribonucleic acid coding regions may have a phenotypic effect on the carrier, e.g. by influencing susceptibility to disease. Detection of deleterious mutations via association studies is hampered by the large number of candidate sites; therefore methods are needed to narrow down the search to the most promising sites. For this, a possible approach is to use structural and sequence-based information of the encoded protein to predict whether a mutation at a particular site is likely to disrupt the functionality of the protein itself. We propose a hierarchical Bayesian multivariate adaptive regression spline (BMARS) model for supervised learning in this context and assess its predictive performance by using data from mutagenesis experiments on lac repressor and lysozyme proteins. In these experiments, about 12 amino-acid substitutions were performed at each native amino-acid position and the effect on protein functionality was assessed. The training data thus consist of repeated observations at each position, which the hierarchical framework is needed to account for. The model is trained on the lac repressor data and tested on the lysozyme mutations and vice versa. In particular, we show that the hierarchical BMARS model, by allowing for the clustered nature of the data, yields lower out-of-sample misclassification rates compared with both a BMARS and a frequen-tist MARS model, a support vector machine classifier and an optimally pruned classification tree.

Conservation of hepatitis C virus nonstructural protein 3 amino acid sequence in viral isolates during liver transplantation

As a consequence of selective pressure exerted by the immune response during hepatitis C virus (HCV) infection, a high rate of nucleotide mutations in the viral genome is observed which leads to the emergence of viral escape mutants. The aim of this study was to evaluate the evolution of the amino acid (aa) sequence of the HCV nonstructural protein 3 (NS3) in viral isolates after liver transplantation. Six patients with HCV-induced liver disease undergoing liver transplantation (LT) were followed up for sequence analysis. Hepatitis C recurrence was observed in all patients after LT. The rate of synonymous (dS) nucleotide substitutions was much higher than that of nonsynonymous (dN) ones in the NS3 encoding region. The high values of the dS/dN ratios suggest no sustained adaptive evolution selection pressure and, therefore, absence of specific NS3 viral populations. Clinical genotype assignments were supported by phylogenetic analysis. Serial samples from each patient showed lower mean nucleotide genetic distance when compared with samples of the same HCV genotype and subtype. The NS3 samples studied had an N-terminal aa sequence with several differences as compared with reference ones, mainly in genotype 1b-infected patients. After LT, as compared with the sequences before, a few reverted aa substitutions and several established aa substitutions were observed at the N-terminal of NS3. Sites described to be involved in important functions of NS3, notably those of the catalytic triad and zinc binding, remained unaltered in terms of aa sequence. Rare or frequent aa substitutions occurred indiscriminately in different positions. Several cytotoxic T lymphocyte epitopes described for HCV were present in our 1b samples. Nevertheless, the deduced secondary structure of the NS3 protease showed a few alterations in samples from genotype 3a patients, but none were seen in 1b cases. Our data, obtained from patients under important selective pressure during LT, show that the NS3 protease remains well conserved, mainly in HCV 3a patients. It reinforces its potential use as an antigenic candidate for further studies aiming at the development of a protective immune response.

Identification of amino acid residues in the alpha, beta, and gamma subunits of the epithelial sodium channel (ENaC) involved in amiloride block and ion permeation.

The amiloride-sensitive epithelial Na channel (ENaC) is a heteromultimeric channel made of three alpha beta gamma subunits. The structures involved in the ion permeation pathway have only been partially identified, and the respective contributions of each subunit in the formation of the conduction pore has not yet been established. Using a site-directed mutagenesis approach, we have identified in a short segment preceding the second membrane-spanning domain (the pre-M2 segment) amino acid residues involved in ion permeation and critical for channel block by amiloride. Cys substitutions of Gly residues in beta and gamma subunits at position beta G525 and gamma G537 increased the apparent inhibitory constant (Ki) for amiloride by > 1,000-fold and decreased channel unitary current without affecting ion selectivity. The corresponding mutation S583 to C in the alpha subunit increased amiloride Ki by 20-fold, without changing channel conducting properties. Coexpression of these mutated alpha beta gamma subunits resulted in a non-conducting channel expressed at the cell surface. Finally, these Cys substitutions increased channel affinity for block by external Zn2+ ions, in particular the alpha S583C mutant showing a Ki for Zn2+ of 29 microM. Mutations of residues alpha W582L, or beta G522D also increased amiloride Ki, the later mutation generating a Ca2+ blocking site located 15% within the membrane electric field. These experiments provide strong evidence that alpha beta gamma ENaCs are pore-forming subunits involved in ion permeation through the channel. The pre-M2 segment of alpha beta gamma subunits may form a pore loop structure at the extracellular face of the channel, where amiloride binds within the channel lumen. We propose that amiloride interacts with Na+ ions at an external Na+ binding site preventing ion permeation through the channel pore.

Amino acids substitutions in σ1 and μ1 outer capsid proteins of a Vero cell-adapted mammalian orthoreovirus are required for optimal virus binding and disassembly

In a recent study, the serotype 3 Dearing strain of mammalian orthoreovirus was adapted to Vero cells; cells that exhibit a limited ability to support the early steps of reovirus uncoating and are unable to produce interferon as an antiviral response upon infection. The Vero cell-adapted virus (VeroAV) exhibits amino acids substitutions in both the σ1 and μ1 outer capsid proteins but no changes in the σ3 protein. Accordingly, the virus was shown not to behave as a classical uncoating mutant. In the present study, an increased ability of the virus to bind at the Vero cell surface was observed and is likely associated with an increased ability to bind onto cell-surface sialic acid residues. In addition, the kinetics of μ1 disassembly from the virions appears to be altered. The plasmid-based reverse genetics approach confirmed the importance of σ1 amino acids substitutions in VeroAV's ability to efficiently infect Vero cells, although μ1 co-adaptation appears necessary to optimize viral infection. This approach of combining in vitro selection of reoviruses with reverse genetics to identify pertinent amino acids substitutions appears promising in the context of eventual reovirus modification to increase its potential as an oncolytic virus.

Mutational analyses of human eIF5A-1-identification of amino acid residues critical for eIF5A activity and hypusine modification

The eukaryotic translation initiation factor 5A (eIF5A) is the only protein that contains hypusine [N-epsilon-(4-amino-2-hydroxybutyl)lysine], which is required for its activity. Hypusine is formed by post-translational modification of one specific lysine (Lys50 for human eIF5A) by deoxyhypusine synthase and deoxyhypusine hydroxylase. To investigate the features of eIF5A required for its activity, we generated 49 mutations in human eIF5A-1, with a single amino acid substitution at the highly conserved residues or with N-terminal or C-terminal truncations, and tested mutant proteins in complementing the growth of a Saccharomyces cerevisiae eIF5A null strain. Growth-supporting activity was abolished in only a few mutant eIF5As (K47D, G49A, K50A, K50D, K50I, K50R, G52A and K55A), with substitutions at or near the hypusine modification site or with truncation of 21 amino acids from either the N-terminus or C-terminus. The inactivity of the Lys50 substitution proteins is obviously due to lack of deoxyhypusine modification. In contrast, K47D and G49A were effective substrates for deoxyhypusine synthase, yet failed to support growth, suggesting critical roles of Lys47 and Gly49 in eIF5A activity, possibly in its interaction with effector(s). By use of a UBHY-R strain harboring genetically engineered unstable eIF5A, we present evidence for the primary function of eIF5A in protein synthesis. When selected eIF5A mutant proteins were tested for their activity in protein synthesis, a close correlation was observed between their ability to enhance protein synthesis and growth, lending further support for a central role of eIF5A in translation.

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^

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

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

Anti-peptide aptamers recognize amino acid sequence and bind a protein epitope.

In vitro selection of nucleic acid binding species (aptamers) is superficially similar to the immune response. Both processes produce biopolymers that can recognize targets with high affinity and specificity. While antibodies are known to recognize the sequence and conformation of protein surface features (epitopes), very little is known about the precise interactions between aptamers and their epitopes. Therefore, aptamers that could recognize a particular epitope, a peptide fragment of human immunodeficiency virus type I Rev, were selected from a random sequence RNA pool. Several of the selected RNAs could bind the free peptide more tightly than a natural RNA ligand, the Rev-binding element. In accord with the hypothesis that protein and nucleic acid binding cusps are functionally similar, interactions between aptamers and the peptide target could be disrupted by sequence substitutions. Moreover, the aptamers appeared to be able to bind peptides with different solution conformations, implying an induced fit mechanism for binding. Just as anti-peptide antibodies can sometimes recognize the corresponding epitope when presented in a protein, the anti-peptide aptamers were found to specifically bind to Rev.

Inducible nitric oxide synthase: identification of amino acid residues essential for dimerization and binding of tetrahydrobiopterin.

Nitric oxide synthases (NOSs) require tetrahydrobiopterin (BH4) for dimerization and NO production. Mutation analysis of mouse inducible NOS (iNOS; NOS2) identified Gly-450 and Ala-453 as critical for NO production, dimer formation, and BH4 binding. Substitutions at five neighboring positions were tolerated, and normal binding of heme, calmodulin, and NADPH militated against major distortions affecting the NH2-terminal portion, midzone, or COOH terminus of the inactive mutants. Direct involvement of residues 450 and 453 in the binding of BH4 is supported by the striking homology of residues 448-480 to a region extensively shared by the three BH4-utilizing aromatic amino acid hydroxylases and is consistent with the conservation of these residues among all 10 reported NOS sequences, including mammalian NOSs 1, 2, and 3, as well as avian and insect NOSs. Altered binding of BH4 and/or L-arginine may explain how the addition of a single methyl group to the side chain of residue 450 or the addition of three methylenes to residue 453 can each abolish an enzymatic activity that reflects the concerted function of 1143 other residues.

Structure-function analysis of amino acid 74 of human RAMP1 and RAMP3 and its role in peptide interactions with adrenomedullin and calcitonin gene-related peptide receptors

The receptors for calcitonin gene-related peptide (CGRP) and adrenomedullin (AM) are complexes of the calcitonin receptor-like receptor (CLR) and receptor activity-modifying proteins (RAMP). The CGRP receptor is a CLR/RAMP1 pairing whereas CLR/RAMP2 and CLR/RAMP3 constitute two subtypes of AM receptor: AM(1) and AM(2), respectively. Previous studies identified Glu74 in RAMP3 to be important for AM binding and potency. To further understand the importance of this residue and its equivalent in RAMP1 (Trp74) we substituted the native amino acids with several others. In RAMP3, these were Trp, Phe, Tyr, Ala, Ser, Thr, Arg and Asn; in RAMP1, Glu, Phe, Tyr, Ala and Asn substitutions were made. The mutant RAMPs were co-expressed with CLR in Cos7 cells; receptor function in response to AM, AM(2)/intermedin and CGRP was measured in a cAMP assay and cell surface expression was determined by ELISA. Phe reduced AM potency in RAMP3 but had no effect in RAMP1. In contrast, Tyr had no effect in RAMP3 but enhanced AM potency in RAMP1. Most other substitutions had a small effect on AM potency in both receptors whereas there was little impact on CGRP or AM(2) potency. Overall, these data suggest that the geometry and charge of the residue at position 74 contribute to how AM interacts with the AM(2) and CGRP receptors and confirms the role of this position in dictating differential AM pharmacology at the AM(2) and CGRP receptors.

Phtx-ii A Basic Myotoxic Phospholipase A2 From Porthidium Hyoprora Snake Venom, Pharmacological Characterization And Amino Acid Sequence By Mass Spectrometry.

A monomeric basic PLA2 (PhTX-II) of 14149.08 Da molecular weight was purified to homogeneity from Porthidium hyoprora venom. Amino acid sequence by in tandem mass spectrometry revealed that PhTX-II belongs to Asp49 PLA2 enzyme class and displays conserved domains as the catalytic network, Ca2+-binding loop and the hydrophobic channel of access to the catalytic site, reflected in the high catalytic activity displayed by the enzyme. Moreover, PhTX-II PLA2 showed an allosteric behavior and its enzymatic activity was dependent on Ca2+. Examination of PhTX-II PLA2 by CD spectroscopy indicated a high content of alpha-helical structures, similar to the known structure of secreted phospholipase IIA group suggesting a similar folding. PhTX-II PLA2 causes neuromuscular blockade in avian neuromuscular preparations with a significant direct action on skeletal muscle function, as well as, induced local edema and myotoxicity, in mice. The treatment of PhTX-II by BPB resulted in complete loss of their catalytic activity that was accompanied by loss of their edematogenic effect. On the other hand, enzymatic activity of PhTX-II contributes to this neuromuscular blockade and local myotoxicity is dependent not only on enzymatic activity. These results show that PhTX-II is a myotoxic Asp49 PLA2 that contributes with toxic actions caused by P. hyoprora venom.

Energy and amino acid content in phase 1 nursery diet: piglet performance and body chemical composition

It was evaluated the effects of metabolizable energy (ME) and digestible lysine (dLYS) densities on performance and body composition of weaned piglets. The study used 114 piglets weaned at 7.4 ± 0.80 kg, out of which 108 were allotted in the nursery and 6 were slaughtered on the weaning day to determine comparative data of body chemical composition. Six nutrients densities were stipulated from a previous study based on the highest nitrogen retention, maintaining the following ME:LYS relationship in the experimental diets: 3,390:1.291; 3,450:1.409; 3,650:1.411; 3,780:1.461; 3,940:1.507; and 4,109 kcal/kg ME:1.564% dLYS. The experimental diets were offered for 13 days when the piglets reached 12.986 ± 1.449 kg of body weight. The probable residual effects of nutritional density on the subsequent performance of the piglets were evaluated. At the end of initial phase 1, six piglets from each density were slaughtered to determine their chemical composition in body fractions and empty body. There was no significant influence of nutritional levels on the performance of the piglets at the end of the evaluation. The results of food conversion and body composition confirm the level indicated in the previous study, 4 g dLYS/Mcal of ME. The increase of energy and lysine densities confirms the need for a correct relationship among both of them to assure better performance of the piglets at the beginning of the growing phase.

A rational protocol for the successful crystallization of l-amino-acid oxidase from Bothrops atrox

Despite the valuable contributions of robotics and high-throughput approaches to protein crystallization, the role of an experienced crystallographer in the evaluation and rationalization of a crystallization process is still crucial to obtaining crystals suitable for X-ray diffraction measurements. In this work, the difficult task of crystallizing the flavoenzyme l-amino-acid oxidase purified from Bothrops atrox snake venom was overcome by the development of a protocol that first required the identification of a non-amorphous precipitate as a promising crystallization condition followed by the implementation of a methodology that combined crystallization in the presence of oil and seeding techniques. Crystals were obtained and a complete data set was collected to 2.3 A resolution. The crystals belonged to space group P2(1), with unit-cell parameters a = 73.64, b = 123.92, c = 105.08 A, beta = 96.03 degrees. There were four protein subunits in the asymmetric unit, which gave a Matthews coefficient V (M) of 2.12 A3 Da-1, corresponding to 42% solvent content. The structure has been solved by molecular-replacement techniques.