Predicted effects of missense mutations on native-state stability account for phenotypic outcome in phenylketonuria, a paradigm of misfolding diseases

Significance of utilizing in silico structural analysis and phenotypic data to characterize phenylalanine hydroxylase variants: A PAH landscape

Phenylketonuria (PKU) is a genetic disorder caused by variations in the phenylalanine hydroxylase (PAH) gene. Among the 3369 reported PAH variants, 33.7% are missense alterations. Unfortunately, 30% of these missense variants are classified as variants of unknown significance (VUS), posing challenges for genetic risk assessment. In our study, we focused on analyzing 836 missense PAH variants following the American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) guidelines specified by ClinGen PAH Variant Curation Expert Panel (VCEP) criteria. We utilized and compared variant annotator tools like Franklin and Varsome, conducted 3D structural analysis of PAH, and examined active and regulatory site hotspots. In addition, we assessed potential splicing effect of apparent missense variants. By evaluating phenotype data from 22962 PKU patients, our aim was to reassess the pathogenicity of missense variants. Our comprehensive approach successfully reclassified 309 VUSs out of 836 missense variants as likely pathogenic or pathogenic (37%), upgraded 370 likely pathogenic variants to pathogenic, and reclassified one previously considered likely benign variant as likely pathogenic. Phenotypic information was available for 636 missense variants, with 441 undergoing 3D structural analysis and active site hotspot identification for 180 variants. After our analysis, only 6% of missense variants were classified as VUSs, and three of them (c.23A>C/p.Asn8Thr, c.59_60delinsCC/p.Gln20Pro, and c.278A >T/p.Asn93Ile) may be influenced by abnormal splicing. Moreover, a pathogenic variant (c.168G>T/p.Glu56Asp) was identified to have a risk exceeding 98% for modifications of the consensus splice site, with high scores indicating a donor loss of 0.94. The integration of ACMG/AMP guidelines with in silico structural analysis and phenotypic data significantly reduced the number of missense VUSs, providing a strong basis for genetic counseling and emphasizing the importance of metabolic phenotype information in variant curation. This study also sheds light on the current landscape of PAH variants.

Delays in Newborn Screening for Phenylketonuria from Birth to Diagnosis and Factors Affecting This

This study aims to evaluate the process of neonatal phenylketonuria (PKU) screening from birth to admission to the pediatric metabolism polyclinic, determining delays in the screening program and the factors influencing them. This study was conducted during 2021-2023. Blood collection dates, results, and probable parameters causing delays in the screening program were recorded. This study included 118 infants. Admission time to the polyclinic was (mean ± SD) 25.2 ± 12.6 days (min-max: 3.4-78.9 days). Admission time was significantly high for refugees, those whose parents were consanguineous, and those who had more heel-prick blood samples taken (p < 0.001, p = 0.005, and p < 0.001, respectively). The first heel-prick blood phenylalanine (phe) level was significantly negatively correlated with the admission time (p < 0.001). Patients' admission time whose first blood phe level < 240 μmol/L was statistically significantly higher than in those with ≥240 μmol/L (p < 0.001). We determined that there were delays in PKU screening from birth to admission to the polyclinic. Being a refugee, the presence of consanguineous marriages, the increase in the number of heel-prick tests, and blood phe levels at a range of 120-240 μmol/L were the factors that played a role in this delay. Taking steps to reduce the impact of these parameters can prevent delays in newborn PKU screening and increase the success of the screening program.

Genetic Landscape and Clinical Features of Hyperphenylalaninemia in North Ossetia-Alania: High Frequency of P281L and P211T Genetic Variants in the PAH Gene

This study, conducted in the Republic of North Ossetia-Alania (RNOA), aimed to explore the genetic landscape of hyperphenylalaninemia (HPA) and phenylketonuria (PKU) in the Ossetian population using data from newborn screening (NBS). Through comprehensive molecular genetic analysis of 29 patients with HPA from diverse ethnic backgrounds, two major genetic variants in the PAH gene, P281L and P211T, were identified, constituting 50% of all detected pathogenic alleles in Ossetian patients. Remarkably, these variants exhibited an exceptionally high frequency in the Ossetian population, surpassing global prevalence rates. This study unveiled a notable prevalence of mild forms of HPA (78%), underscoring the importance of genetic counseling for carriers of pathogenic variants in the PAH gene. Moreover, the findings emphasized the necessity for ongoing monitoring of patients with mild forms, as they may lack significant symptoms for diagnosis, potentially impacting offspring. Overall, this research offers valuable insights into the genetic landscape of HPA and PKU in the Ossetian population.

A novel homozygous missense variant in TBC1D31 in a consanguineous family with congenital anomalies of the kidney and urinary tract (CAKUT)

Congenital anomalies of the kidney and urinary tract (CAKUT) is the leading cause of chronic kidney disease in the first three decades of life. Until now, more than 180 monogenic causes of isolated or syndromic CAKUT have been described. In addition, copy number variants (CNV) have also been implicated, however, all of these causative factors only explain a small fraction of patients with CAKUT, suggesting that additional yet-to-be-discovered novel genes are present. Herein, we report three siblings (two of them are monozygotic twin) of a consanguineous family with CAKUT. Whole-exome sequencing identified a homozygous variant in TBC1D31. Three dimensional protein modeling as well as molecular dynamics simulations predicted it as pathogenic. We therefore showed for the first time an association between a homozygous TBC1D31 variant with CAKUT in humans, expanding its genetic spectrum.

Comparing Supervised Learning and Rigorous Approach for Predicting Protein Stability upon Point Mutations in Difficult Targets

Accurate prediction of protein stability upon a point mutation has important applications in drug discovery and personalized medicine. It remains a challenging issue in computational biology. Existing computational prediction methods, which range from mechanistic to supervised learning approaches, have experienced limited progress over the last few decades. This stagnation is largely due to their heavy reliance on both the quantity and quality of the training data. This is evident in recent state-of-the-art methods that continue to yield substantial errors on two challenging blind test sets: frataxin and p53, with average root-mean-square errors exceeding 3 and 1.5 kcal/mol, respectively, which is still above the theoretical 1 kcal/mol prediction barrier. Rigorous approaches, on the other hand, offer greater potential for accuracy without relying on training data but are computationally demanding and require both wild-type and mutant structure information. Although they showed high accuracy for conserving mutations, their performance is still limited for charge-changing mutation cases. This might be due to the lack of an available mutant structure, often represented by a simplified capped peptide. The recent advances in protein structure prediction methods now make it possible to obtain structures comparable to experimental ones, including complete mutant structure information. In this work, we compare the performance of supervised learning-based methods and rigorous approaches for predicting protein stability on point mutations in difficult targets: frataxin and p53. The rigorous alchemical method significantly surpasses state-of-the-art techniques in terms of both the root-mean-squared error and Pearson correlation coefficient in these two challenging blind test sets. Additionally, we propose an improved alchemical method that employs the pmx double-system/single-box approach to accurately predict the folding free energy change upon both conserving and charge-changing mutations. The enhanced protocol can accurately predict both types of mutations, thereby outperforming existing state-of-the-art methods in overall performance.

Deubiquitinase USP19 enhances phenylalanine hydroxylase protein stability and its enzymatic activity

Phenylalanine hydroxylase (PAH) is the key enzyme in phenylalanine metabolism, deficiency of which is associated with the most common metabolic phenotype of phenylketonuria (PKU) and hyperphenylalaninemia (HPA). A bulk of PKU disease-associated missense mutations in the PAH gene have been studied, and the consequence of each PAH variant vary immensely. Prior research established that PKU-associated variants possess defects in protein folding with reduced cellular stability leading to rapid degradation. However, recent evidence revealed that PAH tetramers exist as a mixture of resting state and activated state whose transition depends upon the phenylalanine concentration and certain PAH variants that fail to modulate the structural equilibrium are associated with PKU disease. Collectively, these findings framed our understanding of the complex genotype-phenotype correlation in PKU. In the current study, we substantiate a link between PAH protein stability and its degradation by the ubiquitin-mediated proteasomal degradation system. Here, we provide an evidence that PAH protein undergoes ubiquitination and proteasomal degradation, which can be reversed by deubiquitinating enzymes (DUBs). We identified USP19 as a novel DUB that regulates PAH protein stability. We found that ectopic expression of USP19 increased PAH protein level, whereas depletion of USP19 promoted PAH protein degradation. Our study indicates that USP19 interacts with PAH and prevents polyubiquitination of PAH subsequently extending the half-life of PAH protein. Finally, the increase in the level of PAH protein by the deubiquitinating activity of USP19 resulted in enhanced metabolic function of PAH. In summary, our study identifies the role of USP19 in regulating PAH protein stability and promotes its metabolic activity. Graphical highlights 1. E3 ligase Cdh1 promotes PAH protein degradation leading to insufficient cellular amount of PAH causing PKU. 2. A balance between E3 ligase and DUB is important to regulate the proteostasis of PAH. 3. USP19 deubiquitinates and stabilizes PAH further protecting it from rapid degradation. 4. USP19 increases the enzymatic activity of PAH, thus maintaining normal Phe levels.

Lynch syndrome, molecular mechanisms and variant classification

Patients with the heritable cancer disease, Lynch syndrome, carry germline variants in the MLH1, MSH2, MSH6 and PMS2 genes, encoding the central components of the DNA mismatch repair system. Loss-of-function variants disrupt the DNA mismatch repair system and give rise to a detrimental increase in the cellular mutational burden and cancer development. The treatment prospects for Lynch syndrome rely heavily on early diagnosis; however, accurate diagnosis is inextricably linked to correct clinical interpretation of individual variants. Protein variant classification traditionally relies on cumulative information from occurrence in patients, as well as experimental testing of the individual variants. The complexity of variant classification is due to (1) that variants of unknown significance are rare in the population and phenotypic information on the specific variants is missing, and (2) that individual variant testing is challenging, costly and slow. Here, we summarise recent developments in high-throughput technologies and computational prediction tools for the assessment of variants of unknown significance in Lynch syndrome. These approaches may vastly increase the number of interpretable variants and could also provide important mechanistic insights into the disease. These insights may in turn pave the road towards developing personalised treatment approaches for Lynch syndrome.

Coevolving stability and activity of LsCR by a single point mutation and constructing neat substrate bioreaction system

Carbonyl reductase (CR)-catalyzed bioreduction in the organic phase and the neat substrate reaction system is a lasting challenge, placing higher requirements on the performance of enzymes. Protein engineering is an effective method to enhance the properties of enzymes for industrial applications. In the present work, a single point mutation E145A on our previously constructed CR mutant LsCR_M3 , coevolved thermostability, and activity. Compared with LsCR_M3 , the catalytic efficiency k_cat /K_M of LsCR_M3 -E145A (LsCR_M4 ) was increased from 6.6 to 21.9 s^-1  mM^-1 . Moreover, E145A prolonged the half-life t_1/2 at 40°C from 4.1 to 117 h, T m ${T}_{m}$ was increased by 5°C, T 50 30 ${T}_{50}^{30}$ was increased by 14.6°C, and T_opt was increased by 15°C. Only 1 g/L of lyophilized Escherichia coli cells expressing LsCR_M4 completely reduced up to 600 g/L 2-chloro-1-(3,4-difluorophenyl)ethanone (CFPO) within 13 h at 45°C, yielding the corresponding (1S)-2-chloro-1-(3,4-difluorophenyl)ethanol ((S)-CFPL) in 99.5% ee_P , with a space-time yield of 1.0 kg/L d, the substrate to catalyst ratios (S/C) of 600 g/g. Compared with LsCR_M3 , the substrate loading was increased by 50%, with the S/C increased by 14 times. Compared with LsCR_WT , the substrate loading was increased by 6.5 times. In contrast, LsCR_M4 completely converted 600 g/L CFPO within 12 h in the neat substrate bioreaction system.

Interpreting the molecular mechanisms of disease variants in human transmembrane proteins

Next-generation sequencing of human genomes reveals millions of missense variants, some of which may lead to loss of protein function and ultimately disease. Here, we investigate missense variants in membrane proteins-key drivers in cell signaling and recognition. We find enrichment of pathogenic variants in the transmembrane region across 19,000 functionally classified variants in human membrane proteins. To accurately predict variant consequences, one fundamentally needs to understand the underlying molecular processes. A key mechanism underlying pathogenicity in missense variants of soluble proteins has been shown to be loss of stability. Membrane proteins, however, are widely understudied. Here, we interpret variant effects on a larger scale by performing structure-based estimations of changes in thermodynamic stability using a membrane-specific energy function and analyses of sequence conservation during evolution of 15 transmembrane proteins. We find evidence for loss of stability being the cause of pathogenicity in more than half of the pathogenic variants, indicating that this is a driving factor also in membrane-protein-associated diseases. Our findings show how computational tools aid in gaining mechanistic insights into variant consequences for membrane proteins. To enable broader analyses of disease-related and population variants, we include variant mappings for the entire human proteome.

Phosphorylation of Thr9 Affects the Folding Landscape of the N-Terminal Segment of Human AGT Enhancing Protein Aggregation of Disease-Causing Mutants

The mutations G170R and I244T are the most common disease cause in primary hyperoxaluria type I (PH1). These mutations cause the misfolding of the AGT protein in the minor allele AGT-LM that contains the P11L polymorphism, which may affect the folding of the N-terminal segment (NTT-AGT). The NTT-AGT is phosphorylated at T9, although the role of this event in PH1 is unknown. In this work, phosphorylation of T9 was mimicked by introducing the T9E mutation in the NTT-AGT peptide and the full-length protein. The NTT-AGT conformational landscape was studied by circular dichroism, NMR, and statistical mechanical methods. Functional and stability effects on the full-length AGT protein were characterized by spectroscopic methods. The T9E and P11L mutations together reshaped the conformational landscape of the isolated NTT-AGT peptide by stabilizing ordered conformations. In the context of the full-length AGT protein, the T9E mutation had no effect on the overall AGT function or conformation, but enhanced aggregation of the minor allele (LM) protein and synergized with the mutations G170R and I244T. Our findings indicate that phosphorylation of T9 may affect the conformation of the NTT-AGT and synergize with PH1-causing mutations to promote aggregation in a genotype-specific manner. Phosphorylation should be considered a novel regulatory mechanism in PH1 pathogenesis.

TMAO to the rescue of pathogenic protein variants

Trimethylamine N-oxide (TMAO) is a chemical chaperone found in various organisms including humans. Various studies unveiled that it is an excellent protein-stabilizing agent, and induces folding of unstructured proteins. It is also well established that it can counteract the deleterious effects of urea, salt, and hydrostatic pressure on macromolecular integrity. There is also existence of large body of data regarding its ability to restore functional deficiency of various mutant proteins or pathogenic variants by correcting misfolding defects and inhibiting the formation of high-order toxic protein oligomers. Since an important class of human disease called "protein conformational disorders" is due to protein misfolding and/or formation of high-order oligomers, TMAO stands as a promising molecule for the therapeutic intervention of such diseases. The present review has been designed to gather a comprehensive knowledge of the TMAO's effect on the functional restoration of various mutants, identify its shortcomings and explore its potentiality as a lead molecule. Future prospects have also been suitably incorporated.

Disease-linked mutations cause exposure of a protein quality control degron

More than half of disease-causing missense variants are thought to lead to protein degradation, but the molecular mechanism of how these variants are recognized by the cell remains enigmatic. Degrons are stretches of amino acids that help mediate recognition by E3 ligases and thus confer protein degradation via the ubiquitin-proteasome system. While degrons that mediate controlled degradation of, for example, signaling components and cell-cycle regulators are well described, so-called protein-quality-control degrons that mediate the degradation of destabilized proteins are poorly understood. Here, we show that disease-linked dihydrofolate reductase (DHFR) missense variants are structurally destabilized and chaperone-dependent proteasome targets. We find two regions in DHFR that act as degrons, and the proteasomal turnover of one of these was dependent on the molecular chaperone Hsp70. Structural analyses by nuclear magnetic resonance (NMR) and hydrogen/deuterium exchange revealed that this degron is buried in wild-type DHFR but becomes transiently exposed in the disease-linked missense variants.

Genome interpretation using in silico predictors of variant impact

Estimating the effects of variants found in disease driver genes opens the door to personalized therapeutic opportunities. Clinical associations and laboratory experiments can only characterize a tiny fraction of all the available variants, leaving the majority as variants of unknown significance (VUS). In silico methods bridge this gap by providing instant estimates on a large scale, most often based on the numerous genetic differences between species. Despite concerns that these methods may lack reliability in individual subjects, their numerous practical applications over cohorts suggest they are already helpful and have a role to play in genome interpretation when used at the proper scale and context. In this review, we aim to gain insights into the training and validation of these variant effect predicting methods and illustrate representative types of experimental and clinical applications. Objective performance assessments using various datasets that are not yet published indicate the strengths and limitations of each method. These show that cautious use of in silico variant impact predictors is essential for addressing genome interpretation challenges.

Predicting and interpreting large-scale mutagenesis data using analyses of protein stability and conservation

Understanding and predicting the functional consequences of single amino acid changes is central in many areas of protein science. Here, we collect and analyze experimental measurements of effects of >150,000 variants in 29 proteins. We use biophysical calculations to predict changes in stability for each variant and assess them in light of sequence conservation. We find that the sequence analyses give more accurate prediction of variant effects than predictions of stability and that about half of the variants that show loss of function do so due to stability effects. We construct a machine learning model to predict variant effects from protein structure and sequence alignments and show how the two sources of information support one another and enable mechanistic interpretations. Together, our results show how one can leverage large-scale experimental assessments of variant effects to gain deeper and general insights into the mechanisms that cause loss of function.

Characterization of phenylalanine hydroxylase gene variants and analysis of genotype-phenotype correlation in patients with phenylalanine hydroxylase deficiency from Fujian Province, Southeastern China

BACKGROUND

Phenylalanine hydroxylase deficiency (PAHD) is the most prevalent inherited disorder of amino acid metabolism in China. Its complex phenotype includes many variants and genotypes among different populations.

METHODS AND RESULTS

In this study, we analyzed the phenylalanine hydroxylase gene (PAH) variants in a cohort of 93 PAHD patients from Fujian Province. We also assessed genotype and phenotype correlation in patients with PAHD. A total of 44 different pathogenic variants were identified, including five novel variants. The three most prevalent variants among all patents were c.158G > A, p.(Arg53His) (18.03%), c.721C > T, p.(Arg241Cys) (14.75%), and c.728G > A, p.(Arg243Gln) (7.65%). The frequency of the c.158G > A, p.(Arg53His) variant was highest in patients with mild hyperphenylalaninemia, whereas the frequency of the c.1197A > T, p.(Val399 =) and c.331C > T, p.(Arg111Ter) variants was highest in patients with classic phenylketonuria. The most abundant genotypes observed in PAHD patients were c.[158G > A];[728G > A], c.[158G > A];[442-1G > A], and c.[158G > A];[721C > T]. Comparing allelic phenotype to genotypic phenotype values yielded fairly accurate predictions of phenotype, with an overall consistency rate was 85.71% for PAHD patients.

CONCLUSIONS

Our study identified a PAH variant spectrum in PAHD patients from Fujian Province, Southeastern China. Quantitative correlation analysis between genotype and phenotype severity is helpful for genetic counseling and management.

Unravelling the Complex Denaturant and Thermal-Induced Unfolding Equilibria of Human Phenylalanine Hydroxylase

Human phenylalanine hydroxylase (PAH) is a metabolic enzyme involved in the catabolism of L-Phe in liver. Loss of conformational stability and decreased enzymatic activity in PAH variants result in the autosomal recessive disorder phenylketonuria (PKU), characterized by developmental and psychological problems if not treated early. One current therapeutic approach to treat PKU is based on pharmacological chaperones (PCs), small molecules that can displace the folding equilibrium of unstable PAH variants toward the native state, thereby rescuing the physiological function of the enzyme. Understanding the PAH folding equilibrium is essential to develop new PCs for different forms of the disease. We investigate here the urea and the thermal-induced denaturation of full-length PAH and of a truncated form lacking the regulatory and the tetramerization domains. For either protein construction, two distinct transitions are seen in chemical denaturation followed by fluorescence emission, indicating the accumulation of equilibrium unfolding intermediates where the catalytic domains are partly unfolded and dissociated from each other. According to analytical centrifugation, the chemical denaturation intermediates of either construction are not well-defined species but highly polydisperse ensembles of protein aggregates. On the other hand, each protein construction similarly shows two transitions in thermal denaturation measured by fluorescence or differential scanning calorimetry, also indicating the accumulation of equilibrium unfolding intermediates. The similar temperatures of mid denaturation of the two constructions, together with their apparent lack of response to protein concentration, indicate the catalytic domains are unfolded in the full-length PAH thermal intermediate, where they remain associated. That the catalytic domain unfolds in the first thermal transition is relevant for the choice of PCs identified in high throughput screening of chemical libraries using differential scanning fluorimetry.

Mapping the degradation pathway of a disease-linked aspartoacylase variant

Canavan disease is a severe progressive neurodegenerative disorder that is characterized by swelling and spongy degeneration of brain white matter. The disease is genetically linked to polymorphisms in the aspartoacylase (ASPA) gene, including the substitution C152W. ASPA C152W is associated with greatly reduced protein levels in cells, yet biophysical experiments suggest a wild-type like thermal stability. Here, we use ASPA C152W as a model to investigate the degradation pathway of a disease-causing protein variant. When we expressed ASPA C152W in Saccharomyces cerevisiae, we found a decreased steady state compared to wild-type ASPA as a result of increased proteasomal degradation. However, molecular dynamics simulations of ASPA C152W did not substantially deviate from wild-type ASPA, indicating that the native state is structurally preserved. Instead, we suggest that the C152W substitution interferes with the de novo folding pathway resulting in increased proteasomal degradation before reaching its stable conformation. Systematic mapping of the protein quality control components acting on misfolded and aggregation-prone species of C152W, revealed that the degradation is highly dependent on the molecular chaperone Hsp70, its co-chaperone Hsp110 as well as several quality control E3 ubiquitin-protein ligases, including Ubr1. In addition, the disaggregase Hsp104 facilitated refolding of aggregated ASPA C152W, while Cdc48 mediated degradation of insoluble ASPA protein. In human cells, ASPA C152W displayed increased proteasomal turnover that was similarly dependent on Hsp70 and Hsp110. Our findings underscore the use of yeast to determine the protein quality control components involved in the degradation of human pathogenic variants in order to identify potential therapeutic targets.

Galactosemia: Towards Pharmacological Chaperones

Galactosemia is a rare inherited metabolic disease resulting from mutations in the four genes which encode enzymes involved in the metabolism of galactose. The current therapy, the removal of galactose from the diet, is inadequate. Consequently, many patients suffer lifelong physical and cognitive disability. The phenotype varies from almost asymptomatic to life-threatening disability. The fundamental biochemical cause of the disease is a decrease in enzymatic activity due to failure of the affected protein to fold and/or function correctly. Many novel therapies have been proposed for the treatment of galactosemia. Often, these are designed to treat the symptoms and not the fundamental cause. Pharmacological chaperones (PC) (small molecules which correct the folding of misfolded proteins) represent an exciting potential therapy for galactosemia. In theory, they would restore enzyme function, thus preventing downstream pathological consequences. In practice, no PCs have been identified for potential application in galactosemia. Here, we review the biochemical basis of the disease, identify opportunities for the application of PCs and describe how these might be discovered. We will conclude by considering some of the clinical issues which will affect the future use of PCs in the treatment of galactosemia.

The effect of improved dietary control on cognitive and psychiatric functioning in adults with phenylketonuria: the ReDAPT study

Background

Phenylketonuria (PKU) is an autosomal recessive inherited disorder characterised by a deficiency in phenylalanine hydroxylase. Untreated, PKU is associated with a wide range of cognitive and psychiatric sequelae. Contemporary management guidelines recommend lifetime dietary control of phenylalanine (Phe) levels, however many individuals who discontinue dietary control subsequently suffer symptoms of anxiety, depression and disturbances to cognition. We undertook a prospective cohort study of patients with early-treated phenylketonuria who had ceased dietary control to test the hypothesis that resumption of dietary control of PKU is associated with improvements in measures of psychiatric morbidity and cognitive functioning.

Methods

We re-initiated dietary control for early-treated patients with PKU and monitored cognitive and psychiatric outcomes over a twelve-month period. Assessments included objective cognitive function (measured by cognitive proficiency index (CPI)), anxiety and depression scales. General linear mixed model (GLMM) analyses were performed to assess change in psychometric variables from baseline over twelve months after resumption of dietary control.

Results

A total of nine patients were recruited. Mean age was 33 years (SD = 8.75), five were female. Mean time off dietary control was 19.1 years (SD = 11.3), and mean baseline phenylalanine (Phe) levels were 1108 µmol/L (SD = 293). GLMM analysis demonstrated a positive relationship between CPI and time on diet (b = 0.56 [95% CI = 0.17, 0.95]). Age, time off diet, Phe levels and depression scores were not associated with cognitive function. There was a negative relationship between time on diet and anxiety (b = − 0.88 95% CI = [− 1.26, − 0.50]) and depression ratings (b = − 0.61, 95% CI = [− 0.95, − 0.26]).

Conclusions

This study demonstrated improvements in cognitive function, anxiety, and depression ratings associated with resumption of dietary control of PKU. Raw Phe levels were not strongly associated with psychiatric or cognitive scores in this cohort. These findings support the importance of lifelong treatment for PKU in improving the cognitive and psychiatric sequelae of the disease.

Understanding protein structural changes for oncogenic missense variants

Understanding and predicting the changes of protein structure and function upon mutation and their relationship to human health is a critical element to translate the genomic revolution into actionable interventions. Therefore, it is pertinent to explore how mutations result in structural changes leading to pathogenic proteins, but due to the protein structural knowledge gap, experimental approaches are lacking. Protein structure prediction methods, such as I-TASSER, have made it possible to predict the structure of a given amino acid sequence, thus opening a new way to explore protein structure changes upon mutations when experimental information is not available. Using known mutations from the Catalogue of Somatic Mutation in Cancer (COSMIC) and ClinVar databases, we compare predicted structure-derived properties from wild type (WT) and mutated proteins and find differences between the local and global 3D protein structures of the WT and the mutants. The studies in this relatively small sample reveal that the structural changes are quite diverse.

Naturally-Occurring Rare Mutations Cause Mild to Catastrophic Effects in the Multifunctional and Cancer-Associated NQO1 Protein

The functional and pathological implications of the enormous genetic diversity of the human genome are mostly unknown, primarily due to our unability to predict pathogenicity in a high-throughput manner. In this work, we characterized the phenotypic consequences of eight naturally-occurring missense variants on the multifunctional and disease-associated NQO1 protein using biophysical and structural analyses on several protein traits. Mutations found in both exome-sequencing initiatives and in cancer cell lines cause mild to catastrophic effects on NQO1 stability and function. Importantly, some mutations perturb functional features located structurally far from the mutated site. These effects are well rationalized by considering the nature of the mutation, its location in protein structure and the local stability of its environment. Using a set of 22 experimentally characterized mutations in NQO1, we generated experimental scores for pathogenicity that correlate reasonably well with bioinformatic scores derived from a set of commonly used algorithms, although the latter fail to semiquantitatively predict the phenotypic alterations caused by a significant fraction of mutations individually. These results provide insight into the propagation of mutational effects on multifunctional proteins, the implementation of in silico approaches for establishing genotype-phenotype correlations and the molecular determinants underlying loss-of-function in genetic diseases.

Towards Accurate Genotype-Phenotype Correlations in the CYP2D6 Gene

Establishing accurate and large-scale genotype-phenotype correlations and predictions of individual response to pharmacological treatments are two of the holy grails of Personalized Medicine. These tasks are challenging and require an integrated knowledge of the complex processes that regulate gene expression and, ultimately, protein functionality in vivo, the effects of mutations/polymorphisms and the different sources of interindividual phenotypic variability. A remarkable example of our advances in these challenging tasks is the highly polymorphic CYP2D6 gene, which encodes a cytochrome P450 enzyme involved in the metabolization of many of the most marketed drugs (including SARS-Cov-2 therapies such as hydroxychloroquine). Since the introduction of simple activity scores (AS) over 10 years ago, its ability to establish genotype-phenotype correlations on the drug metabolizing capacity of this enzyme in human population has provided lessons that will help to improve this type of score for this, and likely many other human genes and proteins. Multidisciplinary research emerges as the best approach to incorporate additional concepts to refine and improve such functional/activity scores for the CYP2D6 gene, as well as for many other human genes associated with simple and complex genetic diseases.

Investigating the structural impacts of a novel missense variant identified with whole exome sequencing in an Egyptian patient with propionic acidemia

Propionic Acidemia (PA) is an inborn error of metabolism caused by variants in the PCCA or PCCB genes, leading to mitochondrial accumulation of propionyl-CoA and its by-products. Here, we report a 2 year-old Egyptian boy with PA who was born to consanguineous parents. Biochemical analysis was performed using tandem mass spectrometry (MS/MS) on the patient's dried blood spots (DBS) followed by urine examination of amino acids using gas chromatography/mass spectrometry (GC/MS). Molecular genetic analysis was carried out using whole-exome sequencing (WES). The PCCA gene sequencing revealed a novel homozygous missense variant affecting the locus (chr13:100962160) of exon 16 of the PCCA gene, resulting in the substitution of the amino acid arginine with proline at site 476 (p.Arg476Pro). Computational analysis revealed that the novel variant might be pathogenic and attributed to decrease the stability and also has an effect on the biotin carboxylase c-terminal domain of the propionyl carboxylase enzyme. The physicochemical properties analysis using NCBI amino acid explorer study revealed restrictions in the side chain and loss of hydrogen bonds due to the variant. On the structural level, the loss of beta-sheet was observed due to the variant proline, which has further led to the loss of surrounding interactions. This loss of beta-sheet and the surrounding interactions might serve the purpose of the structural stability changes. The current study demonstrates that a combination of whole-exome sequencing (WES) and computational analysis are potent tools for validation of diagnosis and classification of disease-causing variants.

Protein Degradation and the Pathologic Basis of Phenylketonuria and Hereditary Tyrosinemia

A delicate intracellular balance among protein synthesis, folding, and degradation is essential to maintaining protein homeostasis or proteostasis, and it is challenged by genetic and environmental factors. Molecular chaperones and the ubiquitin proteasome system (UPS) play a vital role in proteostasis for normal cellular function. As part of protein quality control, molecular chaperones recognize misfolded proteins and assist in their refolding. Proteins that are beyond repair or refolding undergo degradation, which is largely mediated by the UPS. The importance of protein quality control is becoming ever clearer, but it can also be a disease-causing mechanism. Diseases such as phenylketonuria (PKU) and hereditary tyrosinemia-I (HT1) are caused due to mutations in PAH and FAH gene, resulting in reduced protein stability, misfolding, accelerated degradation, and deficiency in functional proteins. Misfolded or partially unfolded proteins do not necessarily lose their functional activity completely. Thus, partially functional proteins can be rescued from degradation by molecular chaperones and deubiquitinating enzymes (DUBs). Deubiquitination is an important mechanism of the UPS that can reverse the degradation of a substrate protein by covalently removing its attached ubiquitin molecule. In this review, we discuss the importance of molecular chaperones and DUBs in reducing the severity of PKU and HT1 by stabilizing and rescuing mutant proteins.

Common mechanism of thermostability in small α- and β-proteins studied by molecular dynamics

Protein thermostability is important to evolution, diseases, and industrial applications. Proteins use diverse molecular strategies to achieve stability at high temperature, yet reducing the entropy of unfolding seems required. We investigated five small α-proteins and five β-proteins with known, distinct structures and thermostability (T_m ) using multi-seed molecular dynamics simulations at 300, 350, and 400 K. The proteins displayed diverse changes in hydrogen bonding, solvent exposure, and secondary structure with no simple relationship to T_m . Our dynamics were in good agreement with experimental B-factors at 300 K and insensitive to force-field choice. Despite the very distinct structures, the native-state (300 + 350 K) free-energy landscapes (FELs) were significantly broader for the two most thermostable proteins and smallest for the three least stable proteins in both the α- and β-group and with both force fields studied independently (tailed t-test, 95% confidence level). Our results suggest that entropic ensembles stabilize proteins at high temperature due to reduced entropy of unfolding, viz., ΔG = ΔH - TΔS. Supporting this mechanism, the most thermostable proteins were also the least kinetically stable, consistent with broader FELs, typified by villin headpiece and confirmed by specific comparison to a mesophilic ortholog of Thermus thermophilus apo-pyrophosphate phosphohydrolase. We propose that molecular strategies of protein thermostabilization, although diverse, tend to converge toward highest possible entropy in the native state consistent with the functional requirements. We speculate that this tendency may explain why many proteins are not optimally structured and why molten-globule states resemble native proteins so much.

Computing the Pathogenicity of Wilson's Disease ATP7B Mutations: Implications for Disease Prevalence

Genetic variations in the gene encoding the copper-transport protein ATP7B are the primary cause of Wilson's disease. Controversially, clinical prevalence seems much smaller than the prevalence estimated by genetic screening tools, causing fear that many people are undiagnosed, although early diagnosis and treatment is essential. To address this issue, we benchmarked 16 state-of-the-art computational disease-prediction methods against established data of missense ATP7B mutations. Our results show that the quality of the methods varies widely. We show the importance of optimizing the threshold of the methods used to distinguish pathogenic from nonpathogenic mutations against data of clinically confirmed pathogenic and nonpathogenic mutations. We find that most methods use thresholds that predict too many ATP7B mutations to be pathogenic. Thus, our findings explain the current controversy on Wilson's disease prevalence because meta-analysis and text search methods include many computational estimates that lead to higher disease prevalence than clinically observed. As proteins and diseases differ widely, a one-size-fits-all threshold cannot distinguish pathogenic and nonpathogenic mutations efficiently, as shown here. We also show that amino acid changes with small evolutionary substitution probability, mainly due to amino acid volume, are more associated with the disease, implying a pathological effect on the conformational state of the protein, which could affect copper transport or adenosine triphosphate recognition and hydrolysis. These findings may be a first step toward a more quantitative genotype-phenotype relationship of Wilson's disease.

Visual impairment and progressive phthisis bulbi caused by recessive pathogenic variant in MARK3

Developmental eye defects often severely reduce vision. Despite extensive efforts, for a substantial fraction of these cases the molecular causes are unknown. Recessive eye disorders are frequent in consanguineous populations and such large families with multiple affected individuals provide an opportunity to identify recessive causative genes. We studied a Pakistani consanguineous family with three affected individuals with congenital vision loss and progressive eye degeneration. The family was analyzed by exome sequencing of one affected individual and genotyping of all family members. We have identified a non-synonymous homozygous variant (NM_001128918.2: c.1708C > G: p.Arg570Gly) in the MARK3 gene as the likely cause of the phenotype. Given that MARK3 is highly conserved in flies (I: 55%; S: 67%) we knocked down the MARK3 homologue, par-1, in the eye during development. This leads to a significant reduction in eye size, a severe loss of photoreceptors and loss of vision based on electroretinogram (ERG) recordings. Expression of the par-1 p.Arg792Gly mutation (equivalent to the MARK3 variant found in patients) in developing fly eyes also induces loss of eye tissue and reduces the ERG signals. The data in flies and human indicate that the MARK3 variant corresponds to a loss of function. We conclude that the identified mutation in MARK3 establishes a new gene-disease link, since it likely causes structural abnormalities during eye development and visual impairment in humans, and that the function of MARK3/par-1 is evolutionarily conserved in eye development.

Computational and cellular studies reveal structural destabilization and degradation of MLH1 variants in Lynch syndrome

Defective mismatch repair leads to increased mutation rates, and germline loss-of-function variants in the repair component MLH1 cause the hereditary cancer predisposition disorder known as Lynch syndrome. Early diagnosis is important, but complicated by many variants being of unknown significance. Here we show that a majority of the disease-linked MLH1 variants we studied are present at reduced cellular levels. We show that destabilized MLH1 variants are targeted for chaperone-assisted proteasomal degradation, resulting also in degradation of co-factors PMS1 and PMS2. In silico saturation mutagenesis and computational predictions of thermodynamic stability of MLH1 missense variants revealed a correlation between structural destabilization, reduced steady-state levels and loss-of-function. Thus, we suggest that loss of stability and cellular degradation is an important mechanism underlying many MLH1 variants in Lynch syndrome. Combined with analyses of conservation, the thermodynamic stability predictions separate disease-linked from benign MLH1 variants, and therefore hold potential for Lynch syndrome diagnostics.

Biophysical and Mechanistic Models for Disease-Causing Protein Variants

The rapid decrease in DNA sequencing cost is revolutionizing medicine and science. In medicine, genome sequencing has revealed millions of missense variants that change protein sequences, yet we only understand the molecular and phenotypic consequences of a small fraction. Within protein science, high-throughput deep mutational scanning experiments enable us to probe thousands of variants in a single, multiplexed experiment. We review efforts that bring together these topics via experimental and computational approaches to determine the consequences of missense variants in proteins. We focus on the role of changes in protein stability as a driver for disease, and how experiments, biophysical models, and computation are providing a framework for understanding and predicting how changes in protein sequence affect cellular protein stability.

Structural and Functional Impact of Seven Missense Variants of Phenylalanine Hydroxylase

The molecular genetics of well-characterized inherited diseases, such as phenylketonuria (PKU) and hyperphenylalaninemia (HPA) predominantly caused by mutations in the phenylalanine hydroxylase (PAH) gene, is often complicated by the identification of many novel variants, often with no obvious impact on the associated disorder. To date, more than 1100 PAH variants have been identified of which a substantial portion have unknown clinical significance. In this work, we study the functionality of seven yet uncharacterized PAH missense variants p.Asn167Tyr, p.Thr200Asn, p.Asp229Gly, p.Gly239Ala, p.Phe263Ser, p.Ala342Pro, and p.Ile406Met first identified in the Czech PKU/HPA patients. From all tested variants, three of them, namely p.Asn167Tyr, p.Thr200Asn, and p.Ile406Met, exerted residual enzymatic activity in vitro similar to wild type (WT) PAH, however, when expressed in HepG2 cells, their protein level reached a maximum of 72.1% ± 4.9%, 11.2% ± 4.2%, and 36.6% ± 7.3% compared to WT PAH, respectively. Remaining variants were null with no enzyme activity and decreased protein levels in HepG2 cells. The chaperone-like effect of applied BH4 precursor increased protein level significantly for p.Asn167Tyr, p.Asp229Gly, p.Ala342Pro, and p.Ile406Met. Taken together, our results of functional characterization in combination with in silico prediction suggest that while p.Asn167Tyr, p.Thr200Asn, and p.Ile406Met PAH variants have a mild impact on the protein, p.Asp229Gly, p.Gly239Ala, p.Phe263Ser, and p.Ala342Pro severely affect protein structure and function.

Structure of full-length human phenylalanine hydroxylase in complex with tetrahydrobiopterin

Phenylalanine hydroxylase (PAH) is a key enzyme in the catabolism of phenylalanine, and mutations in this enzyme cause phenylketonuria (PKU), a genetic disorder that leads to brain damage and mental retardation if untreated. Some patients benefit from supplementation with a synthetic formulation of the cofactor tetrahydrobiopterin (BH4) that partly acts as a pharmacological chaperone. Here we present structures of full-length human PAH (hPAH) both unbound and complexed with BH4 in the precatalytic state. Crystal structures, solved at 3.18-Å resolution, show the interactions between the cofactor and PAH, explaining the negative regulation exerted by BH4 BH4 forms several H-bonds with the N-terminal autoregulatory tail but is far from the catalytic FeII Upon BH4 binding a polar and salt-bridge interaction network links the three PAH domains, explaining the stability conferred by BH4 Importantly, BH4 binding modulates the interaction between subunits, providing information about PAH allostery. Moreover, we also show that the cryo-EM structure of hPAH in absence of BH4 reveals a highly dynamic conformation for the tetramers. Structural analyses of the hPAH:BH4 subunits revealed that the substrate-induced movement of Tyr138 into the active site could be coupled to the displacement of BH4 from the precatalytic toward the active conformation, a molecular mechanism that was supported by site-directed mutagenesis and targeted molecular dynamics simulations. Finally, comparison of the rat and human PAH structures show that hPAH is more dynamic, which is related to amino acid substitutions that enhance the flexibility of hPAH and may increase the susceptibility to PKU-associated mutations.

Mutation spectrum of PAH gene in phenylketonuria patients in Northwest China: identification of twenty novel variants

This study was performed to analyze the mutational spectrum of the phenylalanine hydroxylase (PAH) gene in phenylketonuria (PKU) patients in Northwest China, to identify mutational hot spots, and to determine the correlation between variants and clinical phenotypes of PKU. A large cohort of 475 PKU families in Northwest China was enrolled to analyze PAH gene variants using Sanger sequencing, Multiplex ligation-dependent probe amplification (MLPA), and gap-PCR. Bioinformatics software was used to predict the pathogenicity of novel variants and analyze the correlations between PAH gene variants and phenotypes of PKU patients. A total of 895 variants were detected in the 950 alleles of 475 patients with PKU (detection rate: 94.21%), 20 of which were novel variants. Other 108, previously known variants, were also identified, with the three most frequent variants being p.Arg243Gln (14.00%), c.611A > G (5.58%), and p.Tyr356* (4.95%). Seven different large deletion/duplication variants were identified by the MLPA method, including the large deletion c.-4163_-406del3758 with high frequency. A correlation analysis between patient phenotype and gene variant frequency showed that p.Arg53His and p.Gln419Arg were correlated with mild hyperphenylalaninemia (MHP). In conclusion, the mutational spectrum underlying PKU in Northwest China was established for the first time. Functional analysis of 20 novel PAH gene variants enriched the PAH gene mutational spectrum. Correlation analysis between variants frequencies in compound heterozygous patients and phenotype severity is helpful for phenotypic prediction.

In vitro residual activities in 20 variants of phenylalanine hydroxylase and genotype-phenotype correlation in phenylketonuria patients

Phenylketonuria (PKU), caused by phenylalanine hydroxylase (PAH) gene variants, is a common autosomal inherited metabolic disease. So far, 1111 PAH variants have been revealed. The residual activity of the PAH variants is the key determinant of the metabolic phenotype and BH₄ responsiveness in PKU patients. In this study, the spectrum of PAH variants in 1083 Chinese PKU patients was analyzed. Then 20 variants (p.L52F, p.R86P, p.L128P, p.L142P, p.D163N, p.C203G, p.E214G, p.F260L, p.M276T, p.L311R, p.P314A, p.L364F, p.Q375H, p.F382I, p.A395S, p.V412D, p.E108*, p.C203*, p.C284* and p.E353*) were expressed in COS-7 cells. The residual activities and protein expression levels were detected by isotope-dilution liquid chromatography-electrospray ionization tandem mass spectrometry (LC-ESI-MS/MS) and Western blotting, respectively. We compared the results of the phenotypic prediction based on APV and PAH activity respectively, and further explored the relationship between residual activity and phenotype in PKU patients. We reported 9 newly discovered PAH variants for the first time, thereby expanding the spectrum of PAH variants. Among the 20 variants in our assay, 8 variants showed mild impaired residual activities (48-92%) and approximately normal protein expression levels compared to the wild-type PAH. In contrast, 9 variants showed severely impaired residual activities (0-34%) and reduced protein expression. However, three variants (p.L52F, p.F260L and p.P314A) showed impaired residual activities (5%, 32% and 29%), although the proteins were well expressed. We assigned APV scores for 14 variants, in which the results of the phenotypic prediction were consistent for 12/14 (86%) variants based on APV and residual activity respectively, and the residual activity correctly predicted 17/22 (77%) of the patients. Our study helped to further understand the genotype-phenotype correlation in PKU patients.

NQO1: A target for the treatment of cancer and neurological diseases, and a model to understand loss of function disease mechanisms

NAD(P)H quinone oxidoreductase 1 (NQO1) is a multi-functional protein that catalyses the reduction of quinones (and other molecules), thus playing roles in xenobiotic detoxification and redox balance, and also has roles in stabilising apoptosis regulators such as p53. The structure and enzymology of NQO1 is well-characterised, showing a substituted enzyme mechanism in which NAD(P)H binds first and reduces an FAD cofactor in the active site, assisted by a charge relay system involving Tyr-155 and His-161. Protein dynamics play important role in physio-pathological aspects of this protein. NQO1 is a good target to treat cancer due to its overexpression in cancer cells. A polymorphic form of NQO1 (p.P187S) is associated with increased cancer risk and certain neurological disorders (such as multiple sclerosis and Alzheimer´s disease), possibly due to its roles in the antioxidant defence. p.P187S has greatly reduced FAD affinity and stability, due to destabilization of the flavin binding site and the C-terminal domain, which leading to reduced activity and enhanced degradation. Suppressor mutations partially restore the activity of p.P187S by local stabilization of these regions, and showing long-range allosteric communication within the protein. Consequently, the correction of NQO1 misfolding by pharmacological chaperones is a viable strategy, which may be useful to treat cancer and some neurological conditions, targeting structural spots linked to specific disease-mechanisms. Thus, NQO1 emerges as a good model to investigate loss of function mechanisms in genetic diseases as well as to improve strategies to discriminate between neutral and pathogenic variants in genome-wide sequencing studies.

Genotype-phenotype correlations and BH4 estimated responsiveness in patients with phenylketonuria from Rio de Janeiro, Southeast Brazil

BACKGROUND

Genetic heterogeneity and compound heterozygosis give rise to a continuous spectrum of phenylalanine hydroxylase deficiency and metabolic phenotypes in phenylketonuria (PKU). The most used parameters for evaluating phenotype in PKU are pretreatment phenylalanine (Phe) levels, tolerance for dietary Phe, and Phe overloading test. Phenotype can vary from a "classic" (severe) form to mild hyperphenylalaninemia, which does not require dietary treatment. A subset of patients is responsive to treatment by the cofactor tetrahydrobiopterin (BH4 ). Genotypes of PKU patients from Rio de Janeiro, Brazil, were compared to predicted and observed phenotypes. Genotype-based estimations of responsiveness to BH4 were also conducted.

METHODS

Phenotype was defined by pretreatment Phe levels. A standard prediction system based on arbitrary assigned values was employed to measure genotype-phenotype concordance. Patients were also estimated as BH4 -responders according to the responsiveness previously reported for their mutations and genotypes.

RESULTS

A 48.3% concordance rate between genotype-predicted and observed phenotypes was found. When the predicted phenotypes included those reported at the BIOPKU database, the concordance rate reached 77%. A total of 18 genotypes from 30 patients (29.4%) were estimated as of potential or probable BH4 responsiveness. Inconsistencies were observed in genotypic combinations including the common "moderate" mutations p.R261Q, p.V388M, and p.I65T and the mild mutations p.L48S, p.R68S, and p.L249F.

CONCLUSION

The high discordance rate between genotype-predicted and observed metabolic phenotypes in this study seems to be due partially to the high frequency of the so-called "moderate" common mutations, p.R261Q, p.V388M, and p.I65T, which are reported to be associated to erratic or more severe than expected metabolic phenotypes. Although our results of BH4 estimated responsiveness must be regarded as tentative, it should be emphasized that genotyping and genotype-phenotype association studies are important in selecting patients to be offered a BH4 overload test, especially in low-resource settings like Brazil.

EU-OPENSCREEN: A Novel Collaborative Approach to Facilitate Chemical Biology

Compound screening in biological assays and subsequent optimization of hits is indispensable for the development of new molecular research tools and drug candidates. To facilitate such discoveries, the European Research Infrastructure EU-OPENSCREEN was founded recently with the support of its member countries and the European Commission. Its distributed character harnesses complementary knowledge, expertise, and instrumentation in the discipline of chemical biology from 20 European partners, and its open working model ensures that academia and industry can readily access EU-OPENSCREEN's compound collection, equipment, and generated data. To demonstrate the power of this collaborative approach, this perspective article highlights recent projects from EU-OPENSCREEN partner institutions. These studies yielded (1) 2-aminoquinazolin-4(3 H)-ones as potential lead structures for new antimalarial drugs, (2) a novel lipodepsipeptide specifically inducing apoptosis in cells deficient for the pVHL tumor suppressor, (3) small-molecule-based ROCK inhibitors that induce definitive endoderm formation and can potentially be used for regenerative medicine, (4) potential pharmacological chaperones for inborn errors of metabolism and a familiar form of acute myeloid leukemia (AML), and (5) novel tankyrase inhibitors that entered a lead-to-candidate program. Collectively, these findings highlight the benefits of small-molecule screening, the plethora of assay designs, and the close connection between screening and medicinal chemistry within EU-OPENSCREEN.

Phenylalanine hydroxylase variants interact with the co-chaperone DNAJC12

DNAJC12, a type III member of the HSP40/DNAJ family, has been identified as the specific co-chaperone of phenylalanine hydroxylase (PAH) and the other aromatic amino acid hydroxylases. DNAJ proteins work together with molecular chaperones of the HSP70 family to assist in proper folding and maintenance of intracellular stability of their clients. Autosomal recessive mutations in DNAJC12 were found to reduce PAH levels, leading to hyperphenylalaninemia (HPA) in patients without mutations in PAH. In this work, we investigated the interaction of normal wild-type DNAJC12 with mutant PAH in cells expressing several PAH variants associated with HPA in humans, as well as in the Enu^1/1 mouse model, homozygous for the V106A-Pah variant, which leads to severe protein instability, accelerated PAH degradation and mild HPA. We found that mutant PAH exhibits increased ubiquitination, instability, and aggregation compared with normal PAH. In mouse liver lysates, we showed that DNAJC12 interacts with monoubiquitin-tagged PAH. This form represented a major fraction of PAH in the Enu^1/1 but was also present in liver of wild-type PAH mice. Our results support a role of DNAJC12 in the processing of misfolded ubiquitinated PAH by the ubiquitin-dependent proteasome/autophagy systems and add to the evidence that the DNAJ proteins are important players both for proper folding and degradation of their clients.

Toward mechanistic models for genotype-phenotype correlations in phenylketonuria using protein stability calculations

Phenylketonuria (PKU) is a genetic disorder caused by variants in the gene encoding phenylalanine hydroxylase (PAH), resulting in accumulation of phenylalanine to neurotoxic levels. Here, we analyzed the cellular stability, localization, and interaction with wild-type PAH of 20 selected PKU-linked PAH protein missense variants. Several were present at reduced levels in human cells, and the levels increased in the presence of a proteasome inhibitor, indicating that proteins are proteasome targets. We found that all the tested PAH variants retained their ability to associate with wild-type PAH, and none formed aggregates, suggesting that they are only mildly destabilized in structure. In all cases, PAH variants were stabilized by the cofactor tetrahydrobiopterin (BH₄ ), a molecule known to alleviate symptoms in certain PKU patients. Biophysical calculations on all possible single-site missense variants using the full-length structure of PAH revealed a strong correlation between the predicted protein stability and the observed stability in cells. This observation rationalizes previously observed correlations between predicted loss of protein destabilization and disease severity, a correlation that we also observed using new calculations. We thus propose that many disease-linked PAH variants are structurally destabilized, which in turn leads to proteasomal degradation and insufficient amounts of cellular PAH protein.

Psychiatric and Cognitive Aspects of Phenylketonuria: The Limitations of Diet and Promise of New Treatments

Phenylketonuria (PKU) is a recessive disorder of phenylalanine metabolism due to mutations in the gene for phenylalanine hydroxylase (PAH). Reduced PAH activity results in significant hyperphenylalaninemia, which leads to alterations in cerebral myelin and protein synthesis, as well as reduced levels of serotonin, dopamine, and noradrenaline in the brain. When untreated, brain development is grossly disrupted and significant intellectual impairment and behavioral disturbance occur. The advent of neonatal heel prick screening has allowed for diagnosis at birth, and the institution of a phenylalanine restricted diet. Dietary treatment, particularly when maintained across neurodevelopment and well into adulthood, has resulted in markedly improved outcomes at a cognitive and psychiatric level for individuals with PKU. However, few individuals can maintain full dietary control lifelong, and even with good control, an elevated risk remains of-in particular-mood, anxiety, and attentional disorders across the lifespan. Increasingly, dietary recommendations focus on maintaining continuous dietary treatment lifelong to optimize psychiatric and cognitive outcomes, although the effect of long-term protein restricted diets on brain function remains unknown. While psychiatric illness is very common in adult PKU populations, very little data exist to guide clinicians on optimal treatment. The advent of new treatments that do not require restrictive dietary management, such as the enzyme therapy Pegvaliase, holds the promise of allowing patients a relatively normal diet alongside optimized mental health and cognitive functioning.

Structural and functional insights on the roles of molecular chaperones in the mistargeting and aggregation phenotypes associated with primary hyperoxaluria type I

To carry out their biological function in cells, proteins must be folded and targeted to the appropriate subcellular location. These processes are controlled by a vast collection of interacting proteins collectively known as the protein homeostasis network, in which molecular chaperones play a prominent role. Protein homeostasis can be impaired by inherited mutations leading to genetic diseases. In this chapter, we focus on a particular disease, primary hyperoxaluria type 1 (PH1), in which disease-associated mutations exacerbate protein aggregation in the cell and mistarget the peroxisomal alanine:glyoxylate aminotransferase (AGT) protein to mitochondria, in part due to native state destabilization and enhanced interaction with Hsp60, 70 and 90 chaperone systems. After a general introduction of molecular chaperones and PH1, we review our current knowledge on the structural and energetic features of PH1-causing mutants that lead to these particular pathogenic mechanisms. From this perspective, and in the context of the key role of molecular chaperones in PH1 pathogenesis, we present and discuss current and future perspectives for pharmacological treatments for this disease.

Pharmacoperones as Novel Therapeutics for Diverse Protein Conformational Diseases

After synthesis, proteins are folded into their native conformations aided by molecular chaperones. Dysfunction in folding caused by genetic mutations in numerous genes causes protein conformational diseases. Membrane proteins are more prone to misfolding due to their more intricate folding than soluble proteins. Misfolded proteins are detected by the cellular quality control systems, especially in the endoplasmic reticulum, and proteins may be retained there for eventual degradation by the ubiquitin-proteasome system or through autophagy. Some misfolded proteins aggregate, leading to pathologies in numerous neurological diseases. In vitro, modulating mutant protein folding by altering molecular chaperone expression can ameliorate some misfolding. Some small molecules known as chemical chaperones also correct mutant protein misfolding in vitro and in vivo. However, due to their lack of specificity, their potential as therapeutics is limited. Another class of compounds, known as pharmacological chaperones (pharmacoperones), binds with high specificity to misfolded proteins, either as enzyme substrates or receptor ligands, leading to decreased folding energy barriers and correction of the misfolding. Because many of the misfolded proteins are misrouted but do not have defects in function per se, pharmacoperones have promising potential in advancing to the clinic as therapeutics, since correcting routing may ameliorate the underlying mechanism of disease. This review will comprehensively summarize this exciting area of research, surveying the literature from in vitro studies in cell lines to transgenic animal models and clinical trials in several protein misfolding diseases.

Identification of 34 novel mutations in propionic acidemia: Functional characterization of missense variants and phenotype associations

Propionic acidemia (PA) is caused by mutations in the PCCA and PCCB genes, encoding α and β subunits, respectively, of the mitochondrial enzyme propionyl-CoA carboxylase (PCC). Up to date, >200 pathogenic mutations have been identified, mostly missense defects. Genetic analysis in PA patients referred to the laboratory for the past 15 years identified 20 novel variants in the PCCA gene and 14 in the PCCB gene. 21 missense variants were predicted as probably disease-causing by different bioinformatics algorithms. Structural analysis in the available 3D model of the PCC enzyme indicated potential instability for most of them. Functional analysis in a eukaryotic system confirmed the pathogenic effect for the missense variants and for one amino acid deletion, as they all exhibited reduced or null PCC activity and protein levels compared to wild-type constructs. PCCB variants p.E168del, p.Q58P and p.I460T resulted in medium-high protein levels and no activity. Variants p.R230C and p.C712S in PCCA, and p.G188A, p.R272W and p.H534R in PCCB retained both partial PCC activity and medium-high protein levels. Available patients-derived fibroblasts carriers of some of these mutations were grown at 28 °C or 37 °C and a slight increase in PCC activity or protein could be detected in some cases at the folding-permissive conditions. Examination of available clinical data showed correlation of the results of the functional analysis with disease severity for most mutations, with some notable exceptions, confirming the notion that the final phenotypic outcome in PA is not easily predicted.

Insight into the specificity and severity of pathogenic mechanisms associated with missense mutations through experimental and structural perturbation analyses

Most pathogenic missense mutations cause specific molecular phenotypes through protein destabilization. However, how protein destabilization is manifested as a given molecular phenotype is not well understood. We develop here a structural and energetic approach to describe mutational effects on specific traits such as function, regulation, stability, subcellular targeting or aggregation propensity. This approach is tested using large-scale experimental and structural perturbation analyses in over thirty mutations in three different proteins (cancer-associated NQO1, transthyretin related with amyloidosis and AGT linked to primary hyperoxaluria type I) and comprising five very common pathogenic mechanisms (loss-of-function and gain-of-toxic function aggregation, enzyme inactivation, protein mistargeting and accelerated degradation). Our results revealed that the magnitude of destabilizing effects and, particularly, their propagation through the structure to promote disease-associated conformational states largely determine the severity and molecular mechanisms of disease-associated missense mutations. Modulation of the structural perturbation at a mutated site is also shown to cause switches between different molecular phenotypes. When very common disease-associated missense mutations were investigated, we also found that they were not among the most deleterious possible missense mutations at those sites, and required additional contributions from codon bias and effects of CpG sites to explain their high frequency in patients. Our work sheds light on the molecular basis of pathogenic mechanisms and genotype–phenotype relationships, with implications for discriminating between pathogenic and neutral changes within human genome variability from whole genome sequencing studies.

Biophysical and functional perturbation analyses at cancer-associated P187 and K240 sites of the multifunctional NADP(H):quinone oxidoreductase 1

Once whole-genome sequencing has reached the clinical practice, a main challenge ahead is the high-throughput and accurate prediction of the pathogenicity of genetic variants. However, current prediction tools do not consider explicitly a well-known property of disease-causing mutations: their ability to affect multiple functional sites distant in the protein structure. Here we carried out an extensive biophysical characterization of fourteen mutant variants at two cancer-associated sites of the enzyme NQO1, a paradigm of multi-functional protein. We showed that the magnitude of destabilizing effects, their molecular origins (structural vs. dynamic) and their efficient propagation through the protein structure gradually led to functional perturbations at different sites. Modulation of these structural perturbations also led to switches between molecular phenotypes. Our work supports that experimental and computational perturbation analyses would improve our understanding of the molecular basis of many loss-of-function genetic diseases as well as our ability to accurately predict the pathogenicity of genetic variants in a high-throughput fashion.

Tetrahydrobiopterin treatment reduces brain L-Phe but only partially improves serotonin in hyperphenylalaninemic ENU1/2 mice

Hyperphenylalaninemia (HPA) caused by hepatic phenylalanine hydroxylase (PAH) deficiency has severe consequences on brain monoamine neurotransmitter metabolism. We have studied monoamine neurotransmitter status and the effect of tetrahydrobiopterin (BH4) treatment in Pahenu1/enu2 (ENU1/2) mice, a model of partial PAH deficiency. These mice exhibit elevated blood L-phenylalanine (L-Phe) concentrations similar to that of mild hyperphenylalaninemia (HPA), but brain levels of L-Phe are still ~5-fold elevated compared to wild-type. We found that brain L-tyrosine, L-tryptophan, BH4 cofactor and catecholamine concentrations, and brain tyrosine hydroxylase (TH) activity were normal in these mice but that brain serotonin, 5-hydroxyindolacetic acid (5HIAA) and 3-methoxy-4-hydroxyphenylglycol (MHPG) content, and brain TH protein, as well as tryptophan hydroxylase type 2 (TPH2) protein levels and activity were reduced in comparison to wild-type mice. Parenteral L-Phe loading conditions did not lead to significant changes in brain neurometabolite concentrations. Remarkably, enteral BH4 treatment, which normalized brain L-Phe levels in ENU1/2 mice, lead to only partial recovery of brain serotonin and 5HIAA concentrations. Furthermore, indirect evidence indicated that the GTP cyclohydrolase I (GTPCH) feedback regulatory protein (GFRP) complex may be a sensor for brain L-Phe elevation to ameliorate the toxic effects of HPA. We conclude that BH4 treatment of HPA toward systemic L-Phe lowering reverses elevated brain L-Phe content but the recovery of TPH2 protein and activity as well as serotonin levels is suboptimal, indicating that patients with mild HPA and mood problems (depression or anxiety) treated with the current diet may benefit from supplementation with BH4 and 5-OH-tryptophan.

Melanoma cases demonstrate increased carrier frequency of phenylketonuria/hyperphenylalanemia mutations

Identifying novel melanoma genetic risk factors informs screening and prevention efforts. Mutations in the phenylalanine hydroxylase gene (the causative gene in phenylketonuria) lead to reduced pigmentation in untreated phenylketonuria patients, and reduced pigmentation is associated with greater melanoma risk. Therefore, we sought to characterize the relationship between phenylketonuria carrier status and melanoma risk. Using National Newborn Screening Reports, we determined the United States phenylketonuria/hyperphenylalanemia carrier frequency in Caucasians to be 1.76%. We examined three publically available melanoma datasets for germline mutations in the phenylalanine hydroxylase gene associated with classic phenylketonuria and/or hyperphenylalanemia. Mutations were identified in 29/814 melanoma patients, with a carrier frequency of 3.56%. There was a twofold enrichment (p-value = 3.4 × 10^-5 ) compared to the Caucasian frequency of hyperphenylalanemia/phenylketonuria carriers. These data demonstrate a novel association between phenylalanine hydroxylase carrier status and melanoma risk. Further, functional investigation is warranted to determine the link between phenylalanine hydroxylase mutations and melanomagenesis.

FoldX as Protein Engineering Tool: Better Than Random Based Approaches?

Improving protein stability is an important goal for basic research as well as for clinical and industrial applications but no commonly accepted and widely used strategy for efficient engineering is known. Beside random approaches like error prone PCR or physical techniques to stabilize proteins, e.g. by immobilization, in silico approaches are gaining more attention to apply target-oriented mutagenesis. In this review different algorithms for the prediction of beneficial mutation sites to enhance protein stability are summarized and the advantages and disadvantages of FoldX are highlighted. The question whether the prediction of mutation sites by the algorithm FoldX is more accurate than random based approaches is addressed.