Nucleotides bind to the C-terminus of ClC-5

Mutations in ClC-5 (chloride channel 5), a member of the ClC family of chloride ion channels and antiporters, have been linked to Dent's disease, a renal disease associated with proteinuria. Several of the disease-causing mutations are premature stop mutations which lead to truncation of the C-terminus, pointing to the functional significance of this region. The C-terminus of ClC-5, like that of other eukaryotic ClC proteins, is cytoplasmic and contains a pair of CBS (cystathionine beta-synthase) domains connected by an intervening sequence. The presence of CBS domains implies a regulatory role for nucleotide interaction based on studies of other unrelated proteins bearing these domains [Ignoul and Eggermont (2005) Am. J. Physiol. Cell Physiol. 289, C1369-C1378; Scott, Hawley, Green, Anis, Stewart, Scullion, Norman and Hardie (2004) J. Clin. Invest. 113, 274-284]. However, to date, there has been no direct biochemical or biophysical evidence to support nucleotide interaction with ClC-5. In the present study, we have expressed and purified milligram quantities of the isolated C-terminus of ClC-5 (CIC-5 Ct). CD studies show that the protein is compact, with predominantly alpha-helical structure. We determined, using radiolabelled ATP, that this nucleotide binds the folded protein with low affinity, in the millimolar range, and that this interaction can be competed with 1 muM AMP. CD studies show that binding of these nucleotides causes no significant change in secondary structure, consistent with a model wherein these nucleotides bind to a preformed site. However, both nucleotides induce an increase in thermal stability of ClC-5 Ct, supporting the suggestion that both nucleotides interact with and modify the biophysical properties of this protein.     

Statistical analysis of disease-causing and neutral mutations in human membrane proteins

Mutations in transmembrane proteins (TMPs) have diverse effects on their structure and functions, which may lead to various diseases. In this present study, we have investigated variations in human membrane proteins and found that negatively charged to positively charged/polar and nonpolar to nonpolar changes are dominant in disease-causing and neutral mutations, respectively. Further, we analyzed the top 10 preferred mutations in 14 different disease classes and found that each class has at least two Arg mutations. Moreover, in cardiovascular diseases and congenital disorders of metabolism, Cys mutations occur more frequently in single-pass proteins, whereas Arg and nonpolar residues are more frequently substituted in multi-pass membrane proteins. The immune system diseases are enriched in C → R and C → Y mutations in inside and outside regions. On the other hand, in the membrane region, E → K and R → Q mutations are prevalent. The comparison of mutations in topologically similar regions of globular and membrane proteins showed that Ser and Thr mutations cause deleterious effects in membrane regions, whereas Cys and charged residues, Asp and Arg are prevalent in the buried regions of globular proteins. Our comprehensive analysis of disease-associated mutations in transmembrane proteins will be useful for developing prediction tools.     

Dynamic coupling of residues within proteins as a mechanistic foundation of many enigmatic pathogenic missense variants

Many pathogenic missense mutations are found in protein positions that are neither well-conserved nor fall in any known functional domains. Consequently, we lack any mechanistic underpinning of dysfunction caused by such mutations. We explored the disruption of allosteric dynamic coupling between these positions and the known functional sites as a possible mechanism for pathogenesis. In this study, we present an analysis of 591 pathogenic missense variants in 144 human enzymes that suggests that allosteric dynamic coupling of mutated positions with known active sites is a plausible biophysical mechanism and evidence of their functional importance. We illustrate this mechanism in a case study of β-Glucocerebrosidase (GCase) in which a vast majority of 94 sites harboring Gaucher disease-associated missense variants are located some distance away from the active site. An analysis of the conformational dynamics of GCase suggests that mutations on these distal sites cause changes in the flexibility of active site residues despite their distance, indicating a dynamic communication network throughout the protein. The disruption of the long-distance dynamic coupling caused by missense mutations may provide a plausible general mechanistic explanation for biological dysfunction and disease.     

Predicting Mendelian Disease-Causing Non-Synonymous Single Nucleotide Variants in Exome Sequencing Studies

Exome sequencing is becoming a standard tool for mapping Mendelian disease-causing (or pathogenic) non-synonymous single nucleotide variants (nsSNVs). Minor allele frequency (MAF) filtering approach and functional prediction methods are commonly used to identify candidate pathogenic mutations in these studies. Combining multiple functional prediction methods may increase accuracy in prediction. Here, we propose to use a logit model to combine multiple prediction methods and compute an unbiased probability of a rare variant being pathogenic. Also, for the first time we assess the predictive power of seven prediction methods (including SIFT, PolyPhen2, CONDEL, and logit) in predicting pathogenic nsSNVs from other rare variants, which reflects the situation after MAF filtering is done in exome-sequencing studies. We found that a logit model combining all or some original prediction methods outperforms other methods examined, but is unable to discriminate between autosomal dominant and autosomal recessive disease mutations. Finally, based on the predictions of the logit model, we estimate that an individual has around 5% of rare nsSNVs that are pathogenic and carries ∼22 pathogenic derived alleles at least, which if made homozygous by consanguineous marriages may lead to recessive diseases.     

Genome-wide analysis of human disease alleles reveals that their locations are correlated in paralogous proteins

The millions of mutations and polymorphisms that occur in human populations are potential predictors of disease, of our reactions to drugs, of predisposition to microbial infections, and of age-related conditions such as impaired brain and cardiovascular functions. However, predicting the phenotypic consequences and eventual clinical significance of a sequence variant is not an easy task. Computational approaches have found perturbation of conserved amino acids to be a useful criterion for identifying variants likely to have phenotypic consequences. To our knowledge, however, no study to date has explored the potential of variants that occur at homologous positions within paralogous human proteins as a means of identifying polymorphisms with likely phenotypic consequences. In order to investigate the potential of this approach, we have assembled a unique collection of known disease-causing variants from OMIM and the Human Genome Mutation Database (HGMD) and used them to identify and characterize pairs of sequence variants that occur at homologous positions within paralogous human proteins. Our analyses demonstrate that the locations of variants are correlated in paralogous proteins. Moreover, if one member of a variant-pair is disease-causing, its partner is likely to be disease-causing as well. Thus, information about variant-pairs can be used to identify potentially disease-causing variants, extend existing procedures for polymorphism prioritization, and provide a suite of candidates for further diagnostic and therapeutic purposes.     

Relating the Disease Mutation Spectrum to the Evolution of the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR)

Cystic fibrosis (CF) is the most common genetic disease among Caucasians, and accordingly the cystic fibrosis transmembrane conductance regulator (CFTR) protein has perhaps the best characterized disease mutation spectrum with more than 1,500 causative mutations having been identified. In this study, we took advantage of that wealth of mutational information in an effort to relate site-specific evolutionary parameters with the propensity and severity of CFTR disease-causing mutations. To do this, we devised a scoring scheme for known CFTR disease-causing mutations based on the Grantham amino acid chemical difference matrix. CFTR site-specific evolutionary constraint values were then computed for seven different evolutionary metrics across a range of increasing evolutionary depths. The CFTR mutational scores and the various site-specific evolutionary constraint values were compared in order to evaluate which evolutionary measures best reflect the disease-causing mutation spectrum. Site-specific evolutionary constraint values from the widely used comparative method PolyPhen2 show the best correlation with the CFTR mutation score spectrum, whereas more straightforward conservation based measures (ConSurf and ScoreCons) show the greatest ability to predict individual CFTR disease-causing mutations. While far greater than could be expected by chance alone, the fraction of the variability in mutation scores explained by the PolyPhen2 metric (3.6%), along with the best set of paired sensitivity (58%) and specificity (60%) values for the prediction of disease-causing residues, were marginal. These data indicate that evolutionary constraint levels are informative but far from determinant with respect to disease-causing mutations in CFTR. Nevertheless, this work shows that, when combined with additional lines of evidence, information on site-specific evolutionary conservation can and should be used to guide site-directed mutagenesis experiments by more narrowly defining the set of target residues, resulting in a potential savings of both time and money.     

Protein conformation and disease: pathological consequences of analogous mutations in homologous proteins

The antibody light chain variable domain (V(L))(1) and myelin protein zero (MPZ) are representatives of the functionally diverse immunoglobulin superfamily. The V(L) is a subunit of the antigen-binding component of antibodies, while MPZ is the major membrane-linked constituent of the myelin sheaths that coat peripheral nerves. Despite limited amino acid sequence homology, the conformations of the core structures of the two proteins are largely superimposable. Amino acid variations in V(L) account for various conformational disease outcomes, including amyloidosis. However, the specific amino acid changes in V(L) that are responsible for disease have been obscured by multiple concurrent primary structure alterations. Recently, certain demyelination disorders have been linked to point mutations and single amino acid polymorphisms in MPZ. We demonstrate here that some pathogenic variations in MPZ correspond to changes suspected of determining amyloidosis in V(L). This unanticipated observation suggests that studies of the biophysical origin of conformational disease in one member of a superfamily of homologous proteins may have implications throughout the superfamily. In some cases, findings may account for overt disease; in other cases, due to the natural repertoire of inherited polymorphisms, variations in a representative protein may predict subclinical impairment of homologous proteins.     

CELF proteins regulate CFTR pre-mRNA splicing: essential role of the divergent domain of ETR-3

Cystic fibrosis is a prominent genetic disease caused by mutations of the cystic fibrosis transmembrane conductance regulator (CFTR) gene. Among the many disease-causing alterations are pre-mRNA splicing defects that can hamper mandatory exon inclusion. CFTR exon 9 splicing depends in part on a polymorphic UG(m)U(n) sequence at the end of intron 8, which can be bound by TDP-43, leading to partial exon 9 skipping. CELF proteins, like CUG-BP1 and ETR-3, can also bind UG repeats and regulate splicing. We show here that ETR-3, but not CUG-BP1, strongly stimulates exon 9 skipping, although both proteins bind efficiently to the same RNA motif as TDP-43 and with higher affinity. We further show that the skipping of this exon may be due to the functional antagonism between U2AF⁶⁵ and ETR-3 binding onto the polymorphic U or UG stretch, respectively. Importantly, we demonstrate that the divergent domain of ETR-3 is critical for CFTR exon 9 skipping, as shown by deletion and domain-swapping experiments. We propose a model whereby several RNA-binding events account for the complex regulation of CFTR exon 9 inclusion, with strikingly distinct activities of ETR-3 and CUG-BP1, related to the structure of their divergent domain.     

On the linkage of thermodynamics and pathogenicity

This review outlines the effect of disease-causing mutations on proteins’ thermodynamics. Two major thermodynamics quantities, which are essential for structural integrity, the folding and binding free energy changes caused by missense mutations, are considered. It is emphasized that disease effects in case of complex diseases may originate from several mutations over several genes, while monogenic diseases are caused by mutation is a single gene. Nevertheless, in both cases it is shown that pathogenic mutations cause larger perturbations of the above-mentioned thermodynamics quantities as compared with the benign mutations. Recent works demonstrating the effect of pathogenic mutations on the above-mentioned thermodynamics quantities, as well as on structural dynamics and allosteric pathways, are reviewed.     

Missense mutations in transmembrane domains of proteins: phenotypic propensity of polar residues for human disease

Previous experiments on the cystic fibrosis transmembrane conductance regulator suggested that non-native polar residues within membrane domains can compromise protein structure/function. However, depending on context, replacement of a native residue by a non-native residue can result either in genetic disease or in benign effects (e.g., polymorphisms). Knowledge of missense mutations that frequently cause protein malfunction and subsequent disease can accordingly reveal information as to the impact of these residues in local protein environments. We exploited this concept by performing a statistical comparison of disease-causing mutations in protein membrane-spanning domains versus soluble domains. Using the Human Gene Mutation Database of 240 proteins (including 80 membrane proteins) associated with human disease, we compared the relative phenotypic propensity to cause disease of the 20 naturally occurring amino acids when removed from-or inserted into-native protein sequences. We found that in transmembrane domains (TMDs), mutations involving polar residues, and ionizable residues in particular (notably arginine), are more often associated with protein malfunction than soluble proteins. To further test the hypothesis that interhelical cross-links formed by membrane-embedded polar residues stabilize TMDs, we compared the occurrence of such residues in the TMDs of mesophilic and thermophilic prokaryotes. Results showed a significantly higher proportion of ionizable residues in thermophilic organisms, reinforcing the notion that membrane-embedded electrostatic interactions play critical roles in TMD stability.     

Generalized arterial calcification of infancy and pseudoxanthoma elasticum can be caused by mutations in either ENPP1 or ABCC6

Spontaneous pathologic arterial calcifications in childhood can occur in generalized arterial calcification of infancy (GACI) or in pseudoxanthoma elasticum (PXE). GACI is associated with biallelic mutations in ENPP1 in the majority of cases, whereas mutations in ABCC6 are known to cause PXE. However, the genetic basis in subsets of both disease phenotypes remains elusive. We hypothesized that GACI and PXE are in a closely related spectrum of disease. We used a standardized questionnaire to retrospectively evaluate the phenotype of 92 probands with a clinical history of GACI. We obtained the ENPP1 genotype by conventional sequencing. In those patients with less than two disease-causing ENPP1 mutations, we sequenced ABCC6. We observed that three GACI patients who carried biallelic ENPP1 mutations developed typical signs of PXE between 5 and 8 years of age; these signs included angioid streaks and pseudoxanthomatous skin lesions. In 28 patients, no disease-causing ENPP1 mutation was found. In 14 of these patients, we detected pathogenic ABCC6 mutations (biallelic mutations in eight patients, monoallelic mutations in six patients). Thus, ABCC6 mutations account for a significant subset of GACI patients, and ENPP1 mutations can also be associated with PXE lesions in school-aged children. Based on the considerable overlap of genotype and phenotype of GACI and PXE, both entities appear to reflect two ends of a clinical spectrum of ectopic calcification and other organ pathologies, rather than two distinct disorders. ABCC6 and ENPP1 mutations might lead to alterations of the same physiological pathways in tissues beyond the artery.     

Collagens, the basic proteins of the human body

Collagens are structural elements of many tissues in the human body. The family of collagens can be divided into fibrillar and non-fibrillar collagens. The criterion of the classification is the structure of these proteins. Mutations in the genes encoding collagens cause a variety of human diseases that include osteogenesis imperfecta, some forms of osteoporosis, chondrodysplasias, some types of Ehlers-Danlos syndrome, arterial and intracranial aneurysms, epidermolysis bullosa and the renal disease known as Alport syndrome. The detection of mutations is important both scientifically and clinically. Defining the molecular defects underlying a disorder helps in the understanding of not only the properties of the mutated protein but also the function of the normal protein. Even though many mutations in the genes encoding collagens have been described, the pathogenic consequences of some of the mutations are not fully understood. The important rationale for mutation detection is the clinical use of molecular diagnostics in genetic counselling and differential diagnosis.     

Characterization of intrinsically disordered regions in proteins informed by human genetic diversity

All proteomes contain both proteins and polypeptide segments that don’t form a defined three-dimensional structure yet are biologically active—called intrinsically disordered proteins and regions (IDPs and IDRs). Most of these IDPs/IDRs lack useful functional annotation limiting our understanding of their importance for organism fitness. Here we characterized IDRs using protein sequence annotations of functional sites and regions available in the UniProt knowledgebase (“UniProt features”: active site, ligand-binding pocket, regions mediating protein-protein interactions, etc.). By measuring the statistical enrichment of twenty-five UniProt features in 981 IDRs of 561 human proteins, we identified eight features that are commonly located in IDRs. We then collected the genetic variant data from the general population and patient-based databases and evaluated the prevalence of population and pathogenic variations in IDPs/IDRs. We observed that some IDRs tolerate 2 to 12-times more single amino acid-substituting missense mutations than synonymous changes in the general population. However, we also found that 37% of all germline pathogenic mutations are located in disordered regions of 96 proteins. Based on the observed-to-expected frequency of mutations, we categorized 34 IDRs in 20 proteins (DDX3X, KIT, RB1, etc.) as intolerant to mutation. Finally, using statistical analysis and a machine learning approach, we demonstrate that mutation-intolerant IDRs carry a distinct signature of functional features. Our study presents a novel approach to assign functional importance to IDRs by leveraging the wealth of available genetic data, which will aid in a deeper understating of the role of IDRs in biological processes and disease mechanisms.     

Compensated Pathogenic Deviations: Analysis of Structural Effects

Pathogenic deviations (PDs) in humans are disease-causing missense mutations. However, in some cases, these disease-associated residues occur as the wild-type residues in functionally equivalent proteins in other species and these cases are termed ‘compensated pathogenic deviations’ (CPDs). The lack of pathogenicity in a non-human protein is presumed to be explained in most cases by the presence of compensatory mutations, most commonly within the same protein. Identification of structural features of CPDs and detection of specific compensatory events will help us to understand traversal along fitness landscape valleys in protein evolution.We divided mutations listed in the OMIM (Online Mendelian Inheritance in Man) database into PD and CPD data sets and performed two independent analyses: (i) We searched for potential compensatory mutations spatially close to the CPDs and, (ii) using our SAAPdb database, we examined likely structural effects to try to explain why mutations are pathogenic, comparing PDs and CPDs. Our data sets were obtained from a set of 245 human proteins of known structure and contained a total of 2328 mutations of which 453 (from 85 structures) were seen to be compensated in at least one functionally equivalent protein in another (non-human) species.Structural analysis results confirm previous findings that CPDs are, on average, ‘milder’ in their likely structural effects than uncompensated PDs and tend to be on the protein surface. We also showed that the residues surrounding the CPD residue in the folded protein are more often mutated than the residues surrounding an uncompensated mutation, supporting the hypothesis that compensation is largely a result of structurally local mutations.     

Ryanodine receptors

RyRs are large homotetrameric proteins that are approximately 4/5 cytoplasmic and approximately 1/5 transmembrane and luminal in mass. Mutations in RyRs produce human disease and many of these disease-causing mutations are in the cytoplasmic domains. To elucidate the mechanisms of a disease and to develop interventions, it is crucial to determine how the alterations in the cytoplasmic domains communicate with the transmembrane pore of this channel. One of the major activators of all three RyR isoforms is Ca2+ and some of the disease-causing mutations are thought to alter the sensitivity of the channels to Ca2+ activation. This review examines the current state of structural understanding of the RyR channel activation.     

Expression and characterization of 14 GLB1 mutant alleles found in GM1-gangliosidosis and Morquio B patients

GM1-gangliosidosis and Morquio B disease are lysosomal storage disorders caused by beta-galactosidase deficiency attributable to mutations in the GLB1 gene. On reaching the endosomal-lysosomal compartment, the beta-galactosidase protein associates with the protective protein/cathepsin A (PPCA) and neuraminidase proteins to form the lysosomal multienzyme complex (LMC). The correct interaction of these proteins in the complex is essential for their activity. More than 100 mutations have been described in GM1-gangliosidosis and Morquio B patients, but few have been further characterized. We expressed 12 mutations suspected to be pathogenic, one known polymorphic change (p.S532G), and a variant described as either a pathogenic or a polymorphic change (p.R521C). Ten of them had not been expressed before. The expression analysis confirmed the pathogenicity of the 12 mutations, whereas the relatively high activity of p.S532G is consistent with its definition as a polymorphism. The results for p.R521C suggest that this change is a low-penetrant disease-causing allele. Furthermore, the effect of these beta-galactosidase changes on the LMC was also studied by coimmunoprecipitations and Western blotting. The alteration of neuraminidase and PPCA patterns in several of the Western blotting analyses performed on patient protein extracts indicated that the LMC is affected in at least some GM1-gangliosidosis and Morquio B patients.     

Mutant prion proteins are partially retained in the endoplasmic reticulum

Familial prion diseases are linked to point and insertional mutations in the prion protein (PrP) gene that are presumed to favor conversion of the cellular isoform of PrP to the infectious isoform. In this report, we have investigated the subcellular localization of PrP molecules carrying pathogenic mutations using immunofluorescence staining, immunogold labeling, and PrP-green fluorescent protein chimeras. To facilitate visualization of the mutant proteins, we have utilized a novel Sindbis viral replicon engineered to produce high protein levels without cytopathology. We demonstrate that several different pathogenic mutations have a common effect on the trafficking of PrP, impairing delivery of the molecules to the cell surface and causing a portion of them to accumulate in the endoplasmic reticulum. These observations suggest that protein quality control in the endoplasmic reticulum may play an important role in prion diseases, as it does in some other inherited human disorders. Our experiments also show that chimeric PrP molecules with the sequence of green fluorescent protein inserted adjacent to the glycolipidation site are post-translationally modified and localized normally, thus documenting the utility of these constructs in cell biological studies of PrP.     

Desmosomal protein structure and function and the impact of disease-causing mutations

In this graphical review we focus on the structural characteristics of desmosomal proteins, their interactions with each other and with the intermediate filament cytoskeleton. The wealth of structural information that is now available allows predictions to be made about the pathogenic effect of disease-causing mutations. We have selected representative examples of missense mutations that are buried, semi-buried or surface exposed, and demonstrate how such variants could affect the structural fold of desmosomal proteins that are expressed in the heart. We explain how such alterations could compromise desmosomal adhesion, resulting in life threatening diseases including arrhythmogenic right ventricular cardiomyopathy.     

In silico investigation of the ATP7B gene: insights from functional prediction of non-synonymous substitution to protein structure

ATP7B is a copper-transporting ATPase that plays a key role in the regulation of copper homeostasis. Mutations in the ATP7B gene are causative for Wilson’s disease, and recent reports have suggested that genetic variants are associated with susceptibility to Alzheimer’s disease. Unfortunately, it is difficult to profile experimentally novel genetic variants in the ATP7B gene, because the human protein X-ray structure is not yet entirely understood. In order to investigate ATP7B non-synonymous substitutions, we used an in silico amino acid sequence-based approach. Specifically, we analyzed 337 ATP7B non-synonymous substitutions, which included Wilson’s disease-causing mutations (DVs) and non Wilson’s disease-causing variants (NDVs), with an algorithm that estimated a combined probability (cP_del) of an amino acidic change to be deleterious for the protein function. This approach appeared to reliably indentify the probability of DVs and NDVs to be deleterious and to profile still unknown gene variants. Specifically, after analyzing ATP7B protein domains with the cP_del method, we found results in line with the predicted–modeled domains and some new suggestions. In conclusion, a functional survey of amino acid changes in the ATP7B protein is provided herein, and we suggest that this bioinformatic method can furnish information about novel ATP7B mutations. Furthermore, the same approach can be applied to other uncharacterized proteins.