Investigating the role of single point mutations in the human proteome: a computational study

In the post genomic era with the massive production of biological data the understanding of factors affecting protein stability is one of the most important and challenging tasks for highlighting the role of mutations in relation to human maladies. The problem is at the basis of what is referred to as molecular medicine with the underlying idea that pathologies can be detailed at a molecular level. To this purpose scientific efforts focus on characterising mutations that hamper protein functions and by these affect biological processes at the basis of cell physiology. New techniques have been developed with the aim of detailing single nucleotide polymorphisms (SNPs) at large in all the human chromosomes and by this information in specific databases are exponentially increasing. Eventually mutations that can be found at the DNA level, when occurring in transcribed regions may then lead to mutated proteins and this can be a serious medical problem, largely affecting the phenotype. Bioinformatics tools are urgently needed to cope with the flood of genomic data stored in database and in order to analyse the role of SNPs at the protein level. In principle several experimental and theoretical observations are suggesting that protein stability in the solvent-protein space is responsible of the correct protein functioning. Then mutations that are found disease related during DNA analysis are often assumed to perturb protein stability as well. However so far no extensive analysis at the proteome level has investigated whether this is the case. Also computationally methods have been developed to infer whether a mutation is disease related and independently whether it affects protein stability. Therefore whether the perturbation of protein stability is related to what it is routinely referred to as a disease is still a big question mark. In this work we have tried for the first time to explore the relation among mutations at the protein level and their relevance to diseases with a large-scale computational study of the data from different databases. To this aim in the first part of the thesis for each mutation type we have derived two probabilistic indices (for 141 out of 150 possible SNPs): the perturbing index (Pp), which indicates the probability that a given mutation effects protein stability considering all the “in vitro” thermodynamic data available and the disease index (Pd), which indicates the probability of a mutation to be disease related, given all the mutations that have been clinically associated so far. We find with a robust statistics that the two indexes correlate with the exception of all the mutations that are somatic cancer related. By this each mutation of the 150 can be coded by two values that allow a direct comparison with data base information. Furthermore we also implement computational methods that starting from the protein structure is suited to predict the effect of a mutation on protein stability and find that overpasses a set of other predictors performing the same task. The predictor is based on support vector machines and takes as input protein tertiary structures. We show that the predicted data well correlate with the data from the databases. All our efforts therefore add to the SNP annotation process and more importantly found the relationship among protein stability perturbation and the human variome leading to the diseasome.

Single point mutations affect fatty acid block of human myocardial sodium channel α subunit Na+ channels

Suppression of cardiac voltage-gated Na+ currents is probably one of the important factors for the cardioprotective effects of the n-3 polyunsaturated fatty acids (PUFAs) against lethal arrhythmias. The α subunit of the human cardiac Na+ channel (hH1α) and its mutants were expressed in human embryonic kidney (HEK293t) cells. The effects of single amino acid point mutations on fatty acid-induced inhibition of the hH1α Na+ current (INa) were assessed. Eicosapentaenoic acid (EPA, C20:5n-3) significantly reduced INa in HEK293t cells expressing the wild type, Y1767K, and F1760K of hH1α Na+ channels. The inhibition was voltage and concentration-dependent with a significant hyperpolarizing shift of the steady state of INa. In contrast, the mutant N406K was significantly less sensitive to the inhibitory effect of EPA. The values of the shift at 1, 5, and 10 μM EPA were significantly smaller for N406K than for the wild type. Coexpression of the β1 subunit and N406K further decreased the inhibitory effects of EPA on INa in HEK293t cells. In addition, EPA produced a smaller hyperpolarizing shift of the V1/2 of the steady-state inactivation in HEK293t cells coexpressing the β1 subunit and N406K. These results demonstrate that substitution of asparagine with lysine at the site of 406 in the domain-1-segment-6 region (D1-S6) significantly decreased the inhibitory effect of PUFAs on INa, and coexpression with β1 decreased this effect even more. Therefore, asparagine at the 406 site in hH1α may be important for the inhibition by the PUFAs of cardiac voltage-gated Na+ currents, which play a significant role in the antiarrhythmic actions of PUFAs.

Disruption of RB/E2F-1 interaction by single point mutations in E2F-1 enhances S-phase entry and apoptosis.

The retinoblastoma protein (RB) has been proposed to function as a negative regulator of cell proliferation by complexing with cellular proteins such as the transcription factor E2F. To study the biological consequences of the RB/E2F-1 interaction, point mutants of E2F-1 which fail to bind to RB were isolated by using the yeast two-hybrid system. Sequence analysis revealed that within the minimal 18-amino acid peptide of E2F-1 required for RB binding, five residues, Tyr (position 411), Glu (419), and Asp-Leu-Phe (423-425), are critical. These amino acids are conserved among the known E2F family members. While mutation of any of these five amino acids abolished binding to RB, all mutants retained their full transactivation potential. Expression of mutated E2F-1, when compared with that of wild-type, significantly accelerated entry into S phase and subsequent apoptosis. These results provide direct genetic evidence for the biological significance of the RB/E2F interaction and strongly suggest that the interplay between RB and E2F is critical for proper cell cycle progression.

Detection of point mutations by non-isotopic single strand conformation polymorphism

We have developed a procedure for nonradioactive single strand conformation polymorphism analysis and applied it to the detection of point mutations in the human tumor suppressor gene p53. The protocol does not require any particular facilities or equipment, such as radioactive handling, large gel units for sequencing, or a semiautomated electrophoresis system. This technique consists of amplification of DNA fragments by PCR with specific oligonucleotide primers, denaturation, and electrophoresis on small neutral polyacrylamide gels, followed by silver staining. The sensitivity of this procedure is comparable to other described techniques and the method is easy to perform and applicable to a variety of tissue specimens.

A single point mutation in the pore region of the epithelial Na+ channel changes ion selectivity by modifying molecular sieving.

The epithelial Na+ channel (ENaC) belongs to a new class of channel proteins called the ENaC/DEG superfamily involved in epithelial Na+ transport, mechanotransduction, and neurotransmission. The role of ENaC in Na+ homeostasis and in the control of blood pressure has been demonstrated recently by the identification of mutations in ENaC beta and gamma subunits causing hypertension. The function of ENaC in Na+ reabsorption depends critically on its ability to discriminate between Na+ and other ions like K+ or Ca2+. ENaC is virtually impermeant to K+ ions, and the molecular basis for its high ionic selectivity is largely unknown. We have identified a conserved Ser residue in the second transmembrane domain of the ENaC alpha subunit (alphaS589), which when mutated allows larger ions such as K+, Rb+, Cs+, and divalent cations to pass through the channel. The relative ion permeability of each of the alphaS589 mutants is related inversely to the ionic radius of the permeant ion, indicating that alphaS589 mutations increase the molecular cutoff of the channel by modifying the pore geometry at the selectivity filter. Proper geometry of the pore is required to tightly accommodate Na+ and Li+ ions and to exclude larger cations. We provide evidence that ENaC discriminates between cations mainly on the basis of their size and the energy of dehydration.

Clonal expansion of early to mid-life mitochondrial DNA point mutations drives mitochondrial dysfunction during human ageing

Age-related decline in the integrity of mitochondria is an important contributor to the human ageing process. In a number of ageing stem cell populations, this decline in mitochondrial function is due to clonal expansion of individual mitochondrial DNA (mtDNA) point mutations within single cells. However the dynamics of this process and when these mtDNA mutations occur initially are poorly understood. Using human colorectal epithelium as an exemplar tissue with a well-defined stem cell population, we analysed samples from 207 healthy participants aged 17-78 years using a combination of techniques (Random Mutation Capture, Next Generation Sequencing and mitochondrial enzyme histochemistry), and show that: 1) non-pathogenic mtDNA mutations are present from early embryogenesis or may be transmitted through the germline, whereas pathogenic mtDNA mutations are detected in the somatic cells, providing evidence for purifying selection in humans, 2) pathogenic mtDNA mutations are present from early adulthood (<20 years of age), at both low levels and as clonal expansions, 3) low level mtDNA mutation frequency does not change significantly with age, suggesting that mtDNA mutation rate does not increase significantly with age, and 4) clonally expanded mtDNA mutations increase dramatically with age. These data confirm that clonal expansion of mtDNA mutations, some of which are generated very early in life, is the major driving force behind the mitochondrial dysfunction associated with ageing of the human colorectal epithelium.

A single point mutation in the pore region of the epithelial Na+ channel changes ion selectivity by modifying molecular sieving

The epithelial Na+ channel (ENaC) belongs to a new class of channel proteins called the ENaC/DEG superfamily involved in epithelial Na+ transport, mechanotransduction, and neurotransmission. The role of ENaC in Na+ homeostasis and in the control of blood pressure has been demonstrated recently by the identification of mutations in ENaC β and γ subunits causing hypertension. The function of ENaC in Na+ reabsorption depends critically on its ability to discriminate between Na+ and other ions like K+ or Ca2+. ENaC is virtually impermeant to K+ ions, and the molecular basis for its high ionic selectivity is largely unknown. We have identified a conserved Ser residue in the second transmembrane domain of the ENaC α subunit (αS589), which when mutated allows larger ions such as K+, Rb+, Cs+, and divalent cations to pass through the channel. The relative ion permeability of each of the αS589 mutants is related inversely to the ionic radius of the permeant ion, indicating that αS589 mutations increase the molecular cutoff of the channel by modifying the pore geometry at the selectivity filter. Proper geometry of the pore is required to tightly accommodate Na+ and Li+ ions and to exclude larger cations. We provide evidence that ENaC discriminates between cations mainly on the basis of their size and the energy of dehydration.

Low frequency vibrational modes in proteins: Changes induced by point-mutations in the protein-cofactor matrix of bacterial reaction centers

As a step toward understanding their functional role, the low frequency vibrational motions (<300 cm−1) that are coupled to optical excitation of the primary donor bacteriochlorophyll cofactors in the reaction center from Rhodobacter sphaeroides were investigated. The pattern of hydrogen-bonding interaction between these bacteriochlorophylls and the surrounding protein was altered in several ways by mutation of single amino acids. The spectrum of low frequency vibrational modes identified by femtosecond coherence spectroscopy varied strongly between the different reaction center complexes, including between different mutants where the pattern of hydrogen bonds was the same. It is argued that these variations are primarily due to changes in the nature of the individual modes, rather than to changes in the charge distribution in the electronic states involved in the optical excitation. Pronounced effects of point mutations on the low frequency vibrational modes active in a protein-cofactor system have not been reported previously. The changes in frequency observed indicate a strong involvement of the protein in these nuclear motions and demonstrate that the protein matrix can increase or decrease the fluctuations of the cofactor along specific directions.

Visualization of oligonucleotide probes and point mutations in interphase nuclei and DNA fibers using rolling circle DNA amplification

Rolling circle amplification (RCA) is a surface-anchored DNA replication reaction that can be exploited to visualize single molecular recognition events. Here we report the use of RCA to visualize target DNA sequences as small as 50 nts in peripheral blood lymphocytes or in stretched DNA fibers. Three unique target sequences within the cystic fibrosis transmembrane conductance regulator gene could be detected simultaneously in interphase nuclei, and could be ordered in a linear map in stretched DNA. Allele-discriminating oligonucleotide probes in conjunction with RCA also were used to discriminate wild-type and mutant alleles in the cystic fibrosis transmembrane conductance regulator, p53, BRCA-1, and Gorlin syndrome genes in the nuclei of cultured cells or in DNA fibers. These observations demonstrate that signal amplification by RCA can be coupled to nucleic acid hybridization and multicolor fluorescence imaging to detect single nucleotide changes in DNA within a cytological context or in single DNA molecules. This provides a means for direct physical haplotyping and the analysis of somatic mutations on a cell-by-cell basis.

Novel point mutations in the ERG11 gene in clinical isolates of azole resistant Candida species

The azoles are the class of medications most commonly used to fight infections caused by Candida sp. Typically, resistance can be attributed to mutations in ERG11 gene (CYP51) which encodes the cytochrome P450 14α-demethylase, the primary target for the activity of azoles. The objective of this study was to identify mutations in the coding region of theERG11 gene in clinical isolates of Candidaspecies known to be resistant to azoles. We identified three new synonymous mutations in the ERG11 gene in the isolates of Candida glabrata (C108G, C423T and A1581G) and two new nonsynonymous mutations in the isolates of Candida krusei - A497C (Y166S) and G1570A (G524R). The functional consequence of these nonsynonymous mutations was predicted using evolutionary conservation scores. The G524R mutation did not have effect on 14α-demethylase functionality, while the Y166S mutation was found to affect the enzyme. This observation suggests a possible link between the mutation and dose-dependent sensitivity to voriconazole in the clinical isolate of C. krusei. Although the presence of the Y166S in phenotype of reduced azole sensitivity observed in isolate C. kruseidemands investigation, it might contribute to the search of new therapeutic agents against resistant Candida isolates.

Point Mutations in the Monocarboxylate Transporter SLC16A12 Lead to Juvenile and Age-Related Cataract

Purpose: Previously we reported on a premature termination mutation in SLC16A12 that leads to dominant juvenile cataract and renal glucosuria. To assess the mutation rate and genotype-phenotype correlations of SLC16A12 in juvenile or age-related forms of cataract, we performed a mutation screen in cataract patients. Methods: Clinical data of approximately 660 patients were collected, genomic DNA was isolated and analyzed. Exons 3 to 8 including flanking intron sequences of SLC16A12 were PCR amplified and DNA sequence was determined. Selected mutations were tested by cell culture assays, in silico analysis and RT-PCR. Results: We found sequence alterations at a rate of approximately 1/75 patients. None of them was found in 360 control alleles. Alterations affect splice site and regulatory region but most mutations caused an amino acid substitution. The majority of the coding region mutations maps to trans-membrane domains. One mutation located to the 5'UTR. It affects translational efficiency of SLC16A12. In addition, we identified a cataract-predisposing SNP in the non-coding region that causes allele-specific splicing of the 5'UTR region. Conclusions: Altered translational efficiency of the solute carrier SLC16A12 and its allele-specific splicing strongly support a model of challenged homeostasis to cause various forms of cataract. In addition, the pathogenic property of the here reported sequence alterations is supported by the lack of known sequence variations within the coding region of SLC16A12. Due to the relatively high mutation rate, we suggest to include SLC16A12 in diagnostic cataract screening. Generally, our data recommend the assessment of regulatory sequences for diagnostic purposes.

Gain of function mutations in CgPDR1 of Candida glabrata not only mediate antifungal resistance but also enhance virulence.

CgPdr1p is a Candida glabrata Zn(2)-Cys(6) transcription factor involved in the regulation of the ABC-transporter genes CgCDR1, CgCDR2, and CgSNQ2, which are mediators of azole resistance. Single-point mutations in CgPDR1 are known to increase the expression of at least CgCDR1 and CgCDR2 and thus to contribute to azole resistance of clinical isolates. In this study, we investigated the incidence of CgPDR1 mutations in a large collection of clinical isolates and tested their relevance, not only to azole resistance in vitro and in vivo, but also to virulence. The comparison of CgPDR1 alleles from azole-susceptible and azole-resistant matched isolates enabled the identification of 57 amino acid substitutions, each positioned in distinct CgPDR1 alleles. These substitutions, which could be grouped into three different "hot spots," were gain of function (GOF) mutations since they conferred hyperactivity to CgPdr1p revealed by constitutive high expression of ABC-transporter genes. Interestingly, the major transporters involved in azole resistance (CgCDR1, CgCDR2, and CgSNQ2) were not always coordinately expressed in presence of specific CgPDR1 GOF mutations, thus suggesting that these are rather trans-acting elements (GOF in CgPDR1) than cis-acting elements (promoters) that lead to azole resistance by upregulating specific combinations of ABC-transporter genes. Moreover, C. glabrata isolates complemented with CgPDR1 hyperactive alleles were not only more virulent in mice than those with wild type alleles, but they also gained fitness in the same animal model. The presence of CgPDR1 hyperactive alleles also contributed to fluconazole treatment failure in the mouse model. In conclusion, this study shows for the first time that CgPDR1 mutations are not only responsible for in vitro/in vivo azole resistance but that they can also confer a selective advantage under host conditions.

Single Point of Entry Long-Term Living Resource System Team Report

Pursuant to House File 451 the Single Point of Entry Long-Term Living Resources System Team, involving several state agencies as well as interested associations, submitted a report to the legislature on recommendations to establish a single point of entry system.

Probing the role of point mutations in the cyp51A gene from Aspergillus fumigatus in the model yeast Saccharomyces cerevisiae.

Azole-resistant strains of Aspergillus fumigatus have been detected and the underlying molecular mechanisms of resistance characterized. Point mutations in the cyp51A gene have been proved to be related to azole resistance in A. fumigatus clinical strains and with different resistance profiles depending on the amino acid change (G54E, G54V, G54R, G54W, M220V, M220K, M220T, M220I). The aim of this work was to express A. fumigatus cyp51A genes in the yeast Saccharomyces cerevisiae in order to better assess the contribution of each independent amino acid substitution to resistance. A tetracycline regulatable system allowing repression of the endogenous essential ERG11 gene was used. The expression of Aspergillus cyp51A alleles could efficiently restore the absence of ERG11 in S. cerevisiae. In general, S. cerevisiae clones expressing. A. fumigatus cyp51A alleles from azole-resistant isolates showed higher MICs to all azoles tested than those expressing alleles from susceptible isolates. The azole susceptibility profiles obtained in S. cerevisiae upon expression of specific cyp51A alleles recapitulated susceptibility profiles observed from their A. fumigatus origins. In conclusion this work supports the concept that characteristics of specific A. fumigatus cyp51A alleles could be investigated in the heterologous host S. cerevisiae.