Missense3D-DB web catalogue: an atom-based analysis and repository of 4M human protein-coding genetic variants

Allele-specific detection of isoniazid metabolism modulating variants of N-acetyltransferase 2 enzyme and their frequencies in the Bangladeshi population

Tuberculosis is one of the oldest diseases that still affects millions of people worldwide and remains a significant public health challenge. The N-acetyltransferase 2 (NAT2) enzyme metabolizes isoniazid (INH), a primary antibiotic in tuberculosis treatment. The single nucleotide polymorphisms (SNPs) of NAT2 affect the metabolism and function of isoniazid. The rs1801280 (T341C) and rs1208 (G803A) variants of NAT2 are associated with INH drug responses. Individuals with the slow-metabolizing rs1801280 variant of the NAT2 enzyme are at a higher risk of INH-induced liver damage and require lower doses or longer treatment regimens. At the same time, individuals with the fast-metabolizing rs1208 variant are at risk of treatment failure due to rapid drug metabolism. Genotyping of the NAT2 variants can help clinicians personalize tuberculosis treatment, optimize drug doses, and thus minimize adverse effects. Under this study, an allele-specific PCR (ASPCR) method was developed for genotyping the NAT2 variants, and the results were validated through targeted sequencing. The allele frequencies at the rs1801280 locus were 0.60 for the T allele and 0.40 for the C. For rs1208, the participants' allele frequencies were 0.27 for the G allele and 0.73 for the A allele. This ASPCR method is quick, affordable, and could be used in routine genotyping to personalize the treatment for tuberculosis patients, leading to more effective and safer treatments. We also used molecular docking to study how the rs1801280 and rs1208 variants affect the interaction between the NAT2 enzyme and drugs. A slight change was visible in the flexibility of the amino acid residues. However, those amino acids were not involved in the ligand binding mechanism.

Identification of a binding pocket of letermovir in the terminase subunit pUL56 of human cytomegalovirus

A key step in replication of human cytomegalovirus (HCMV) is the generation and packaging of unit-length genomes into preformed capsids. The enzymes involved in this process are viral terminases. The HCMV terminase consists of two subunits, the ATPase pUL56 and the nuclease pUL89. A potential third component, pUL51, has been proposed. Letermovir is the first terminase inhibitor available for HCMV prophylaxis to allogenic hematopoietic stem cell recipients. However, mutations in the HCMV terminase subunit pUL56 and, to a lesser extent, in pUL89 or pUL51 lead to resistance. Here we focused on the drug target area in the terminase subunit pUL56. To gain further structural insights into the putative binding site of letermovir, in silico analysis of the structure was performed using Phyre2 and SwissDock. For our analysis, we used three of the most frequent mutations during letermovir treatment, C325F, C325Y and C325W. We demonstrated that all variants have a pronounced cavity reduction, leading to the letermovir binding conformations being "pushed-out" of the binding pocket. This results in a changed distribution of the Gibbs free binding energy. To circumvent the absolute resistance of C325 mutations a further modification of letermovir might solve the problem and leads to optimizing drug targeting capacity.

Comprehensive computational analysis of deleterious nsSNPs in PTEN gene for structural and functional insights

Single nucleotide polymorphisms (SNPs) are pivotal in understanding the genetic basis of complex disorders. Among them, nonsynonymous SNPs (nsSNPs) that alter amino acid sequences can significantly impact protein structure and function. This study focuses on analyzing deleterious nsSNPs in the tumor suppressor gene PTEN (Phosphatase and TENsin Homolog), which plays a central role in regulating the PI3K/Akt signaling pathway and tumorigenesis. Out of 43,855 SNPs in PTEN, 17 deleterious nsSNPs were identified using six computational tools. Protein stability analysis revealed that 15 variants reduce stability, potentially leading to functional impairment. Structural evaluations using HOPE and ConSurf classified mutations into buried structural residues disrupting protein integrity and exposed functional residues affecting molecular interactions. STRING database analysis highlighted PTEN as a central node in an intricate protein network, with deleterious mutations impairing critical interactions with partners such as PIK3CA, AKT1, and TP53. Secondary structure analysis revealed distinct structural deviations, particularly for G129E, which exhibited the most pronounced destabilization. Molecular dynamics simulations confirmed stability variations across mutants, with G129E exhibiting greater instability. This comprehensive analysis enhances understanding of PTEN nsSNP impacts, offering insights for therapeutic interventions and future experimental validation.

Pathogenic single nucleotide polymorphisms in RhoA gene: Insights into structural and functional impacts on RhoA-PLD1 interaction through molecular dynamics simulation

Molecular switches serve as key regulators of biological systems by acting as one of the crucial driving forces in the initiation of signal transduction pathway cascades. The Ras homolog gene family member A (RhoA) is one of the molecular switches that binds with GTP in order to cycle between an active GTP-bound state and an inactive GDP-bound state. Any aberrance in control over this circuit, particularly due to any perturbation in switching, leads to the development of different pathogenicity. Consequently, the single nucleotide polymorphisms (SNPs) within the RhoA gene, especially deleterious genetic variations, are crucial to study to forecast structural alteration and their functional impacts in light of disease onset. In this comprehensive study, we employed a range of computational tools to screen the deleterious SNPs of RhoA from 207 nonsynonymous SNPs (nsSNPs). By utilizing 7 distinct tools for further analysis, 8 common deleterious SNPs were sorted, among them 5 nsSNPs (V9G, G17E, E40K, A61T, F171L) were found to be in the highly conserved regions, with E40K and A61T at G2 and G3 motif of the GTP-binding domain respectively, indicating potential perturbation in GTP/GDP binding ability of the protein. RhoA-GDP complex interacts with the enzyme phospholipase, specifically PLD1, to regulate different cellular activities. PLD1 is also a crucial regulator of thrombosis and cancer. In that line of focus, our initial structural analysis of Y66H, A61T, G17E, I86N, and I151T mutations of RhoA revealed remarkable decreased hydrophobicity from which we further filtered out G17E and I86N which may have potential impact on the RhoA-GDP-PLD1 complex. Intriguingly, the comparative 250 ns (ns) molecular dynamics (MD) simulation of these two mutated complexes revealed overall structural instability and altered interaction patterns. Therefore, further investigation into these deleterious mutations with in vitro and in vivo studies could lead to the identification of potential biomarkers in terms of different pathogenesis and could also be utilized in personalized therapeutic targets in the long run.

LTβR deficiency causes lymph node aplasia and impaired B cell differentiation

Secondary lymphoid organs (SLOs) provide the confined microenvironment required for stromal cells to interact with immune cells to initiate adaptive immune responses resulting in B cell differentiation. Here, we studied three patients from two families with functional hyposplenism, absence of tonsils, and complete lymph node aplasia, leading to recurrent bacterial and viral infections. We identified biallelic loss-of-function mutations in LTBR, encoding the lymphotoxin beta receptor (LTβR), primarily expressed on stromal cells. Patients with LTβR deficiency had hypogammaglobulinemia, diminished memory B cells, regulatory and follicular T helper cells, and dysregulated expression of several tumor necrosis factor family members. B cell differentiation in an ex vivo coculture system was intact, implying that the observed B cell defects were not intrinsic in nature and instead resulted from LTβR-dependent stromal cell interaction signaling critical for SLO formation. Collectively, we define a human inborn error of immunity caused primarily by a stromal defect affecting the development and function of SLOs.

Discovery of a novel mutation F184S (c.551T>C) in GATA4 gene causing congenital heart disease in a consanguineous Saudi family

Background & aim

Congenital heart disease (CHD) is the most common cause of non-infectious deaths in infants worldwide. However, the molecular mechanisms underlying CHD remain unclear. Approximately 30 % of the causes are believed to be genetic mutations and chromosomal abnormalities. In this study, we aimed to identify the genetic causes of CHD in consanguineous families.

Methods

Fourth-generation pedigrees with CHD were recruited. The main cardiac features of the patient included absence of the right pulmonary artery and a large dilated left pulmonary artery. To determine the underlying genetic cause, whole-exome sequencing was performed and subsequently confirmed using Sanger sequencing and different online databases to study the pathogenesis of the identified gene mutation. An in-silico homology model was created using the Alphafold homology model structure of GATA4 (AF-P43694-F1). The missense3D online program was used to evaluate the structural alterations.

Results

We identified a deleterious mutation c.551T > C (p.Phe184Ser) in GATA4. GATA4 is a highly conserved zinc-finger transcription factor, and its continuous expression is essential for cardiogenesis during embryogenesis. The in-silico model suggested a compromised binding efficiency with other proteins. Several variant interpretation algorithms indicated that the F184S missense variant in GATA4 is damaging, whereas HOPE analysis indicated the functional impairment of DNA binding of transcription factors and zinc-ion binding activities of GATA4.

Conclusion

The variant identified in GATA4 appears to cause recessive CHD in the family. In silico analysis suggested that this variant was damaging and caused multiple structural and functional aberrations. This study may support prenatal screening of the fetus in this family to prevent diseases in new generations.

Unveiling New Genetic Variants Associated with Age at Onset in Alzheimer's Disease and Frontotemporal Lobar Degeneration Due to C9orf72 Repeat Expansions

Alzheimer's disease (AD) and Frontotemporal lobar degeneration (FTLD) represent the most common forms of neurodegenerative dementias with a highly phenotypic variability. Herein, we investigated the role of genetic variants related to the immune system and inflammation as genetic modulators in AD and related dementias. In patients with sporadic AD/FTLD (n = 300) and GRN/C9orf72 mutation carriers (n = 80), we performed a targeted sequencing of 50 genes belonging to the immune system and inflammation, selected based on their high expression in brain regions and low tolerance to genetic variation. The linear regression analyses revealed two genetic variants: (i) the rs1049296 in the transferrin (TF) gene, shown to be significantly associated with age at onset in the sporadic AD group, anticipating the disease onset of 4 years for each SNP allele with respect to the wild-type allele, and (ii) the rs7550295 in the calsyntenin-1 (CLSTN1) gene, which was significantly associated with age at onset in the C9orf72 group, delaying the disease onset of 17 years in patients carrying the SNP allele. In conclusion, our data support the role of genetic variants in iron metabolism (TF) and in the modulation of the calcium signalling/axonal anterograde transport of vesicles (CLSTN1) as genetic modulators in AD and FTLD due to C9orf72 expansions.

ProtVar: mapping and contextualizing human missense variation

Genomic variation can impact normal biological function in complex ways and so understanding variant effects requires a broad range of data to be coherently assimilated. Whilst the volume of human variant data and relevant annotations has increased, the corresponding increase in the breadth of participating fields, standards and versioning mean that moving between genomic, coding, protein and structure positions is increasingly complex. In turn this makes investigating variants in diverse formats and assimilating annotations from different resources challenging. ProtVar addresses these issues to facilitate the contextualization and interpretation of human missense variation with unparalleled flexibility and ease of accessibility for use by the broadest range of researchers. By precalculating all possible variants in the human proteome it offers near instantaneous mapping between all relevant data types. It also combines data and analyses from a plethora of resources to bring together genomic, protein sequence and function annotations as well as structural insights and predictions to better understand the likely effect of missense variation in humans. It is offered as an intuitive web server https://www.ebi.ac.uk/protvar where data can be explored and downloaded, and can be accessed programmatically via an API.

Biallelic variants in COX18 cause a mitochondrial disorder primarily manifesting as peripheral neuropathy

Defects in mitochondrial dynamics are a common cause of Charcot-Marie-Tooth disease (CMT), while primary deficiencies in the mitochondrial respiratory chain (MRC) are rare and atypical for this etiology. This study aims to report COX18 as a novel CMT-causing gene. This gene encodes an assembly factor of mitochondrial Complex IV (CIV) that translocates the C-terminal tail of MTCO2 across the mitochondrial inner membrane. Exome sequencing was performed in four affected individuals. The patients and available family members underwent thorough neurological and electrophysiological assessment. The impact of one of the identified variants on splicing, protein levels, and mitochondrial bioenergetics was investigated in patient-derived lymphoblasts. The functionality of the mutant protein was assessed using a Proteinase K protection assay and immunoblotting. Neuronal relevance of COX18 was assessed in a Drosophila melanogaster knockdown model. Exome sequencing coupled with homozygosity mapping revealed a homozygous splice variant c.435-6A>G in COX18 in two siblings with early-onset progressive axonal sensory-motor peripheral neuropathy. By querying external databases, we identified two additional families with rare deleterious biallelic variants in COX18 . All affected individuals presented with axonal CMT and some patients also exhibited central nervous system symptoms, such as dystonia and spasticity. Functional characterization of the c.435-6A>G variant demonstrated that it leads to the expression of an alternative transcript that lacks exon 2, resulting in a stable but defective COX18 isoform. The mutant protein impairs CIV assembly and activity, leading to a reduction in mitochondrial membrane potential. Downregulation of the COX18 homolog in Drosophila melanogaster displayed signs of neurodegeneration, including locomotor deficit and progressive axonal degeneration of sensory neurons. Our study presents genetic and functional evidence that supports COX18 as a newly identified gene candidate for autosomal recessive axonal CMT with or without central nervous system involvement. These findings emphasize the significance of peripheral neuropathy within the spectrum of primary mitochondrial disorders and the role of mitochondrial CIV in the development of CMT. Our research has important implications for the diagnostic workup of CMT patients.

GTExome: Modeling commonly expressed missense mutations in the human genome

A web application, GTExome, is described that quickly identifies, classifies, and models missense mutations in commonly expressed human proteins. GTExome can be used to categorize genomic mutation data with tissue specific expression data from the Genotype-Tissue Expression (GTEx) project. Commonly expressed missense mutations in proteins from a wide range of tissue types can be selected and assessed for modeling suitability. Information about the consequences of each mutation is provided to the user including if disulfide bonds, hydrogen bonds, or salt bridges are broken, buried prolines introduced, buried charges are created or lost, charge is swapped, a buried glycine is replaced, or if the residue that would be removed is a proline in the cis configuration. Also, if the mutation site is in a binding pocket the number of pockets and their volumes are reported. The user can assess this information and then select from available experimental or computationally predicted structures of native proteins to create, visualize, and download a model of the mutated protein using Fast and Accurate Side-chain Protein Repacking (FASPR). For AlphaFold modeled proteins, confidence scores for native proteins are provided. Using this tool, we explored a set of 9,666 common missense mutations from a variety of tissues from GTEx and show that most mutations can be modeled using this tool to facilitate studies of protein-protein and protein-drug interactions. The open-source tool is freely available at https://pharmacogenomics.clas.ucdenver.edu/gtexome/.

The natural history and genotype-phenotype correlations of TMPRSS3 hearing loss: an international, multi-center, cohort analysis

TMPRSS3-related hearing loss presents challenges in correlating genotypic variants with clinical phenotypes due to the small sample sizes of previous studies. We conducted a cross-sectional genomics study coupled with retrospective clinical phenotype analysis on 127 individuals. These individuals were from 16 academic medical centers across 6 countries. Key findings revealed 47 unique TMPRSS3 variants with significant differences in hearing thresholds between those with missense variants versus those with loss-of-function genotypes. The hearing loss progression rate for the DFNB8 subtype was 0.3 dB/year. Post-cochlear implantation, an average word recognition score of 76% was observed. Of the 51 individuals with two missense variants, 10 had DFNB10 with profound hearing loss. These 10 all had at least one of 4 TMPRSS3 variants predicted by computational modeling to be damaging to TMPRSS3 structure and function. To our knowledge, this is the largest study of TMPRSS3 genotype-phenotype correlations. We find significant differences in hearing thresholds, hearing loss progression, and age of presentation, by TMPRSS3 genotype and protein domain affected. Most individuals with TMPRSS3 variants perform well on speech recognition tests after cochlear implant, however increased age at implant is associated with worse outcomes. These findings provide insight for genetic counseling and the on-going design of novel therapeutic approaches.

A rare inherited homozygous missense variant in PLA2G6 influences susceptibility to infantile neuroaxonal dystrophy: a case report

Background

Infantile neuroaxonal dystrophy (INAD) is an ultra-rare early-onset autosomal recessive neurodegenerative disorder due to PLA2G6 variants. The clinical symptoms of INAD patients display considerable diversity, and many PLA2G6 variants are still not thoroughly investigated in relation to their associated clinical presentations.

Case Description

A 16-month-old boy was admitted to our hospital due to regression of acquired motor and speech abilities that had persisted for 4 months. The patient was born to a healthy consanguineous couple after 41 weeks of pregnancy and natural delivery. Before 12 months old, he had normal motor development milestones. On admission, he also showed astasia-abasia, weakness of distal muscles, and diminished patellar tendon reflex. Brain magnetic resonance imaging (MRI) revealed cerebellar atrophy. Auditory brainstem response (ABR) indicated moderately severe hearing loss. With chromosome microarray analysis (CMA), we identified several copy number-neutral regions of runs of homozygosity (ROH) in the patient. Whole-exome sequencing (WES) further revealed that the patient harbored a homozygous missense variant NM_003560.2: c.1778C>T, p.Pro593Leu (rs1451486649) in the PLA2G6 gene. In the patient's asymptomatic parents and brother, the PLA2G6 c.1778C>T variant stayed in heterozygous status as confirmed by Sanger sequencing. The patient was finally diagnosed with INAD.

Conclusions

We report an INAD child with a rare PLA2G6 c.1778C>T homozygous missense variant and associated clinical symptoms. The family-based cosegregation analysis reveals that the PLA2G6 c.1778C>T homozygous variant contributes to the pathogenesis of INAD.

Revisiting Exome Data Identified Missed Splice Site Variant of the Asparagine Synthetase ( ASNS ) Gene

Next-generation sequencing, such as whole-exome sequencing (WES), is increasingly used in the study of Mendelian disorders, yet many are reported as "negative." Inappropriate variant annotation and filtering steps are reasons for missing the molecular diagnosis. Noncoding variants, including splicing mutations, are examples of variants that can be overlooked. Herein, we report a family of four affected newborns, and all presented with severe congenital microcephaly. Initial research WES analysis identified a damaging homozygous variant in NME1 gene as a possible cause of primary microcephaly phenotype in these patients. However, reanalysis of the exome data uncovered a biallelic splice site variant in asparagine synthetase gene which seems to be the possible cause of the phenotype in these patients. This study highlights the importance of revisiting the exome data and the issue of "negative" exome and the afterward approaches to identify and prove new candidate genes.

Biallelic variants in GTF3C5, a regulator of RNA polymerase III-mediated transcription, cause a multisystem developmental disorder

General transcription factor IIIC subunit 5 (GTF3C5) encodes transcription factor IIIC63 (TFIIIC63). It binds to DNA to recruit another transcription factor, TFIIIB, and RNA polymerase III (Pol III) to mediate the transcription of small noncoding RNAs, such as tRNAs. Here, we report four individuals from three families presenting with a multisystem developmental disorder phenotype with biallelic variants in GTF3C5. The overlapping features include growth retardation, developmental delay, intellectual disability, dental anomalies, cerebellar malformations, delayed bone age, skeletal anomalies, and facial dysmorphism. Using lymphoblastoid cell lines (LCLs) from two affected individuals, we observed a reduction in TFIIIC63 protein levels compared to control LCLs. Genome binding of TFIIIC63 protein is also reduced in LCL from one of the affected individuals. Additionally, approximately 40% of Pol III binding regions exhibited reduction in the level of Pol III occupancy in the mutant genome relative to the control, while approximately 54% of target regions showed comparable levels of Pol III occupancy between the two, indicating partial impairment of Pol III occupancy in the mutant genome. Yeasts with subject-specific variants showed temperature sensitivity and impaired growth, supporting the notion that the identified variants have deleterious effects. gtf3c5 mutant zebrafish showed developmental defects, including a smaller body, head, and eyes. Taken together, our data show that GTF3C5 plays an important role in embryonic development, and that biallelic variants in this gene cause a multisystem developmental disorder. Our study adds GTF3C5-related disorder to the growing list of genetic disorders associated with Pol III transcription machinery.

Structural effect of the H992D/H418D mutation of angiotensin-converting enzyme in the Indian population: implications for health and disease

The Non synonymous SNPs (nsSNPs) of the renin-angiotensin-system (RAS) pathway, unique to the Indian population were investigated in view of its importance as an endocrine system. nsSNPs of the RAS pathway genes were mined from the IndiGenome database. Damaging nsSNPs were predicted using SIFT, PredictSNP, SNP and GO, Snap2 and Protein Variation Effect Analyzer. Loss of function was predicted based on protein stability change using I mutant, PremPS and CONSURF. The structural impact of the nsSNPs was predicted using HOPE and Missense3d followed by modeling, refinement, and energy minimization. Molecular Dynamics studies were carried out using Gromacsv2021.1. 23 Indian nsSNPs of the RAS pathway genes were selected for structural analysis and 8 were predicted to be damaging. Further sequence analysis showed that HEMGH zinc binding motif changes to HEMGD in somatic ACE-C domain (sACE-C) H992D and Testis ACE (tACE) H418D resulted in loss of zinc coordination, which is essential for enzymatic activity in this metalloprotease. There was a loss of internal interactions around the zinc coordination residues in the protein structural network. This was also confirmed by Principal Component Analysis, Free Energy Landscape and residue contact maps. Both mutations lead to broadening of the AngI binding cavity. The H992D mutation in sACE-C is likely to be favorable for cardiovascular health, but may lead to renal abnormalities with secondary impact on the heart. H418D in tACE is potentially associated with male infertility.

Missense3D-TM: Predicting the Effect of Missense Variants in Helical Transmembrane Protein Regions Using 3D Protein Structures

Variant effect predictors assess if a substitution is pathogenic or benign. Most predictors, including those that are structure-based, are designed for globular proteins in aqueous environments and do not consider that the variant residue is located within the membrane. We report Missense3D-TM that provides a structure-based assessment of the impact of a missense variant located within a membrane. On a dataset of 2,078 pathogenic and 1,060 benign variants, spanning 711 proteins from 706 structures, Missense3D-TM achieved an accuracy of 66%, Mathews correlation coefficient of 0.37, sensitivity of 58% and specificity of 81%. Missense3D-TM performed similarly to mCSM-membrane: accuracy 66% vs 61% (p = 0.02) on an unbalanced test set and 70% vs 67% (p = 0.20) on a balanced test set. The Missense3D-TM website provides an analysis of the structural effects of the variant along with its predicted position within the membrane. The web server is available at http://missense3d.bc.ic.ac.uk/.

Protein engineering in the computational age: An open source framework for exploring mutational landscapes in silico

The field of protein engineering has seen tremendous expansion in the last decade, with researchers developing novel proteins with specialised functionalities for a range of uses, from drug discovery to industrial biotechnology. The emergence of computational tools and high-throughput screening technology has substantially sped up the process of protein engineering. However, much of the expertise required to engage in such projects is still concentrated in the hands of a few specialised individuals, including computational biologists and structural biochemists. The international Genetically Engineered Machine (iGEM) competition represents a platform for undergraduate students to innovate in synthetic biology. Yet, due to their complexity, arduous protein engineering projects are hindered by the resources available and strict timelines of the competition. The authors highlight how the 2022 iGEM Team, 'Sporadicate', set out to develop InFinity 1.0, a computational framework for increased accessibility to effective protein engineering, hoping to increase awareness and accessibility to novel in silico tools.

GTExome: Modeling commonly expressed missense mutations in the human genome

A web application, GTExome, is described that quickly identifies, classifies, and models missense mutations in commonly expressed human proteins. GTExome can be used to categorize genomic mutation data with tissue specific expression data from the Genotype-Tissue Expression (GTEx) project. Commonly expressed missense mutations in proteins from a wide range of tissue types can be selected and assessed for modeling suitability. Information about the consequences of each mutation is provided to the user including if disulfide bonds, hydrogen bonds, or salt bridges are broken, buried prolines introduced, buried charges are created or lost, charge is swapped, a buried glycine is replaced, or if the residue that would be removed is a proline in the cis configuration. Also, if the mutation site is in a binding pocket the number of pockets and their volumes are reported. The user can assess this information and then select from available experimental or computationally predicted structures of native proteins to create, visualize, and download a model of the mutated protein using Fast and Accurate Side-chain Protein Repacking (FASPR). For AlphaFold modeled proteins, confidence scores for native proteins are provided. Using this tool, we explored a set of 9,666 common missense mutations from a variety of tissues from GTEx and show that most mutations can be modeled using this tool to facilitate studies of protein-protein and protein-drug interactions. The open-source tool is freely available at https://pharmacogenomics.clas.ucdenver.edu/gtexome/.

Epithelial IFNγ signalling and compartmentalized antigen presentation orchestrate gut immunity

All nucleated cells express major histocompatibility complex I and interferon-γ (IFNγ) receptor1, but an epithelial cell-specific function of IFNγ signalling or antigen presentation by means of major histocompatibility complex I has not been explored. We show here that on sensing IFNγ, colonic epithelial cells productively present pathogen and self-derived antigens to cognate intra-epithelial T cells, which are critically located at the epithelial barrier. Antigen presentation by the epithelial cells confers extracellular ATPase expression in cognate intra-epithelial T cells, which limits the accumulation of extracellular adenosine triphosphate and consequent activation of the NLRP3 inflammasome in tissue macrophages. By contrast, antigen presentation by the tissue macrophages alongside inflammasome-associated interleukin-1α and interleukin-1β production promotes a pathogenic transformation of CD4+ T cells into granulocyte-macrophage colony-stimulating-factor (GM-CSF)-producing T cells in vivo, which promotes colitis and colorectal cancer. Taken together, our study unravels critical checkpoints requiring IFNγ sensing and antigen presentation by epithelial cells that control the development of pathogenic CD4+ T cell responses in vivo.

DICER1 platform domain missense variants inhibit miRNA biogenesis and lead to tumor susceptibility

The endoribonuclease DICER1 plays an essential role in the microRNA (miRNA) biogenesis pathway, cleaving precursor miRNA (pre-miRNA) stem-loops to generate mature single-stranded miRNAs. Germline pathogenic variants (GPVs) in DICER1 result in DICER1 tumor predisposition syndrome (DTPS), a mainly childhood-onset tumor susceptibility disorder. Most DTPS-causing GPVs are nonsense or frameshifting, with tumor development requiring a second somatic missense hit that impairs the DICER1 RNase IIIb domain. Interestingly, germline DICER1 missense variants that cluster in the DICER1 Platform domain have been identified in some persons affected by tumors that also associate with DTPS. Here, we demonstrate that four of these Platform domain variants prevent DICER1 from producing mature miRNAs and as a result impair miRNA-mediated gene silencing. Importantly, we show that in contrast to canonical somatic missense variants that alter DICER1 cleavage activity, DICER1 proteins harboring these Platform variants fail to bind to pre-miRNA stem-loops. Taken together, this work sheds light upon a unique subset of GPVs causing DTPS and provides new insights into how alterations in the DICER1 Platform domain can impact miRNA biogenesis.

Landscape of mSWI/SNF chromatin remodeling complex perturbations in neurodevelopmental disorders

DNA sequencing-based studies of neurodevelopmental disorders (NDDs) have identified a wide range of genetic determinants. However, a comprehensive analysis of these data, in aggregate, has not to date been performed. Here, we find that genes encoding the mammalian SWI/SNF (mSWI/SNF or BAF) family of ATP-dependent chromatin remodeling protein complexes harbor the greatest number of de novo missense and protein-truncating variants among nuclear protein complexes. Non-truncating NDD-associated protein variants predominantly disrupt the cBAF subcomplex and cluster in four key structural regions associated with high disease severity, including mSWI/SNF-nucleosome interfaces, the ATPase-core ARID-armadillo repeat (ARM) module insertion site, the Arp module and DNA-binding domains. Although over 70% of the residues perturbed in NDDs overlap with those mutated in cancer, ~60% of amino acid changes are NDD-specific. These findings provide a foundation to functionally group variants and link complex aberrancies to phenotypic severity, serving as a resource for the chromatin, clinical genetics and neurodevelopment communities.

Missense3D-PPI: A Web Resource to Predict the Impact of Missense Variants at Protein Interfaces Using 3D Structural Data

In 2019, we released Missense3D which identifies stereochemical features that are disrupted by a missense variant, such as introducing a buried charge. Missense3D analyses the effect of a missense variant on a single structure and thus may fail to identify as damaging surface variants disrupting a protein interface i.e., a protein-protein interaction (PPI) site. Here we present Missense3D-PPI designed to predict missense variants at PPI interfaces. Our development dataset comprised of 1,279 missense variants (pathogenic n = 733, benign n = 546) in 434 proteins and 545 experimental structures of PPI complexes. Benchmarking of Missense3D-PPI was performed after dividing the dataset in training (320 benign and 320 pathogenic variants) and testing (226 benign and 413 pathogenic). Structural features affecting PPI, such as disruption of interchain bonds and introduction of unbalanced charged interface residues, were analysed to assess the impact of the variant at PPI. The performance of Missense3D-PPI was superior to that of Missense3D: sensitivity 44 % versus 8% and accuracy 58% versus 40%, p = 4.23 × 10-16. However, the specificity of Missense3D-PPI was lower compared to Missense3D (84% versus 98%). On our dataset, Missense3D-PPI's accuracy was superior to BeAtMuSiC (p = 3.4 × 10-5), mCSM-PPI2 (p = 1.5 × 10-12) and MutaBind2 (p = 0.0025). Missense3D-PPI represents a valuable tool for predicting the structural effect of missense variants on biological protein networks and is available at the Missense3D web portal (http://missense3d.bc.ic.ac.uk).

A genotype-to-phenotype approach suggests under-reporting of single nucleotide variants in nephrocystin-1 (NPHP1) related disease (UK 100,000 Genomes Project)

Autosomal recessive whole gene deletions of nephrocystin-1 (NPHP1) result in abnormal structure and function of the primary cilia. These deletions can result in a tubulointerstitial kidney disease known as nephronophthisis and retinal (Senior-Løken syndrome) and neurological (Joubert syndrome) diseases. Nephronophthisis is a common cause of end-stage kidney disease (ESKD) in children and up to 1% of adult onset ESKD. Single nucleotide variants (SNVs) and small insertions and deletions (Indels) have been less well characterised. We used a gene pathogenicity scoring system (GenePy) and a genotype-to-phenotype approach on individuals recruited to the UK Genomics England (GEL) 100,000 Genomes Project (100kGP) (n = 78,050). This approach identified all participants with NPHP1-related diseases reported by NHS Genomics Medical Centres and an additional eight participants. Extreme NPHP1 gene scores, often underpinned by clear recessive inheritance, were observed in patients from diverse recruitment categories, including cancer, suggesting the possibility of a more widespread disease than previously appreciated. In total, ten participants had homozygous CNV deletions with eight homozygous or compound heterozygous with SNVs. Our data also reveals strong in-silico evidence that approximately 44% of NPHP1 related disease may be due to SNVs with AlphaFold structural modelling evidence for a significant impact on protein structure. This study suggests historical under-reporting of SNVS in NPHP1 related diseases compared with CNVs.

Protein structure-based evaluation of missense variants: Resources, challenges and future directions

We provide an overview of the methods that can be used for protein structure-based evaluation of missense variants. The algorithms can be broadly divided into those that calculate the difference in free energy (ΔΔG) between the wild type and variant structures and those that use structural features to predict the damaging effect of a variant without providing a ΔΔG. A wide range of machine learning approaches have been employed to develop those algorithms. We also discuss challenges and opportunities for variant interpretation in view of the recent breakthrough in three-dimensional structural modelling using deep learning.

Orthogonal analysis of variants in APOE gene using in-silico approaches reveals novel disrupting variants

Introduction: Alzheimer's disease (AD) is one of the most prominent medical conditions in the world. Understanding the genetic component of the disease can greatly advance our knowledge regarding its progression, treatment and prognosis. Single amino-acid variants (SAVs) in the APOE gene have been widely investigated as a risk factor for AD Studies, including genome-wide association studies, meta-analysis based studies, and in-vivo animal studies, were carried out to investigate the functional importance and pathogenesis potential of APOE SAVs. However, given the high cost of such large-scale or experimental studies, there are only a handful of variants being reported that have definite explanations. The recent development of in-silico analytical approaches, especially large-scale deep learning models, has opened new opportunities for us to probe the structural and functional importance of APOE variants extensively. Method: In this study, we are taking an ensemble approach that simultaneously uses large-scale protein sequence-based models, including Evolutionary Scale Model and AlphaFold, together with a few in-silico functional prediction web services to investigate the known and possibly disease-causing SAVs in APOE and evaluate their likelihood of being functional and structurally disruptive. Results: As a result, using an ensemble approach with little to no prior field-specific knowledge, we reported 5 SAVs in APOE gene to be potentially disruptive, one of which (C112R) was classificed by previous studies as a key risk factor for AD. Discussion: Our study provided a novel framework to analyze and prioritize the functional and structural importance of SAVs for future experimental and functional validation.

Rare variant analyses across multiethnic cohorts identify novel genes for refractive error

Refractive error, measured here as mean spherical equivalent (SER), is a complex eye condition caused by both genetic and environmental factors. Individuals with strong positive or negative values of SER require spectacles or other approaches for vision correction. Common genetic risk factors have been identified by genome-wide association studies (GWAS), but a great part of the refractive error heritability is still missing. Some of this heritability may be explained by rare variants (minor allele frequency [MAF] ≤ 0.01.). We performed multiple gene-based association tests of mean Spherical Equivalent with rare variants in exome array data from the Consortium for Refractive Error and Myopia (CREAM). The dataset consisted of over 27,000 total subjects from five cohorts of Indo-European and Eastern Asian ethnicity. We identified 129 unique genes associated with refractive error, many of which were replicated in multiple cohorts. Our best novel candidates included the retina expressed PDCD6IP, the circadian rhythm gene PER3, and P4HTM, which affects eye morphology. Future work will include functional studies and validation. Identification of genes contributing to refractive error and future understanding of their function may lead to better treatment and prevention of refractive errors, which themselves are important risk factors for various blinding conditions.

Homozygous GRHPR C.494G>A mutation is deleterious that causes early onset of nephrolithiasis in West Bengal, India

Pediatric nephrolithiasis (NL) or Kidney stone disease (KSD) is an untethered topic in Asian population. In Western countries, the annual incidence of paediatric NL is around 6-10%. Here, we present data from West Bengal, India, on lower age (LA, 0-20 years) NL and its prevalence for the first time. To discover the mutations associated with KSD, twenty-four (18 + 6) rare LA-NL patients were selected for Whole Exome Sequencing (WES) and Sanger sequencing, respectively. It was found that GRHPR c. 494G>A mutation (MZ826703) is predominant in our study cohort. This specific homozygous mutation is functionally studied for the first time directly from human peripheral mononuclear cell (PBMC) samples. Using expression study with biochemical activity and computational analysis we assumed that the mutation is pathogenic with loss of function. Moreover, three genes, AGXT, HOGA1 and GRHPR with Novel variants known to cause hyperoxaluria were found frequently in the study cohort. Our study analyses the genes and variations that cause LA-NL, as well as the molecular function of the GRHPR mutation, which may serve as a clinical marker in the population of West Bengal, Eastern India.

Predicting and Understanding the Pathology of Single Nucleotide Variants in Human COQ Genes

Coenzyme Q (CoQ) is a vital lipid that functions as an electron carrier in the mitochondrial electron transport chain and as a membrane-soluble antioxidant. Deficiencies in CoQ lead to metabolic diseases with a wide range of clinical manifestations. There are currently few treatments that can slow or stop disease progression. Primary CoQ₁₀ deficiency can arise from mutations in any of the COQ genes responsible for CoQ biosynthesis. While many mutations in these genes have been identified, the clinical significance of most of them remains unclear. Here we analyzed the structural and functional impact of 429 human missense single nucleotide variants (SNVs) that give rise to amino acid substitutions in the conserved and functional regions of human genes encoding a high molecular weight complex known as the CoQ synthome (or Complex Q), consisting of the COQ3-COQ7 and COQ9 gene products. Using structures of COQ polypeptides, close homologs, and AlphaFold models, we identified 115 SNVs that are potentially pathogenic. Further biochemical characterizations in model organisms such as Saccharomyces cerevisiae are required to validate the pathogenicity of the identified SNVs. Collectively, our results will provide a resource for clinicians during patient diagnosis and guide therapeutic efforts toward combating primary CoQ₁₀ deficiency.

PDBe and PDBe-KB: Providing high-quality, up-to-date and integrated resources of macromolecular structures to support basic and applied research and education

The archiving and dissemination of protein and nucleic acid structures as well as their structural, functional and biophysical annotations is an essential task that enables the broader scientific community to conduct impactful research in multiple fields of the life sciences. The Protein Data Bank in Europe (PDBe; pdbe.org) team develops and maintains several databases and web services to address this fundamental need. From data archiving as a member of the Worldwide PDB consortium (wwPDB; wwpdb.org), to the PDBe Knowledge Base (PDBe-KB; pdbekb.org), we provide data, data-access mechanisms, and visualizations that facilitate basic and applied research and education across the life sciences. Here, we provide an overview of the structural data and annotations that we integrate and make freely available. We describe the web services and data visualization tools we offer, and provide information on how to effectively use or even further develop them. Finally, we discuss the direction of our data services, and how we aim to tackle new challenges that arise from the recent, unprecedented advances in the field of structure determination and protein structure modeling.

Interpreting the spectrum of gamma-secretase complex missense variation in the context of hidradenitis suppurativa—An in-silico study

Hidradenitis suppurativa (HS) is a disease of the pilosebaceous unit characterized by recurrent nodules, abscesses and draining tunnels with a predilection to intertriginous skin. The pathophysiology of HS is complex. However, it is known that inflammation and hyperkeratinization at the hair follicle play crucial roles in disease manifestation. Genetic and environmental factors are considered the main drivers of these two pathophysiological processes. Despite a considerable proportion of patients having a positive family history of disease, only a minority of patients suffering from HS have been found to harbor monogenic variants which segregate to affected kindreds. Most of these variants are in the ɣ secretase complex (GSC) protein-coding genes. In this manuscript, we set out to characterize the burden of missense pathogenic variants in healthy reference population using large scale genomic dataset thereby providing a standard for comparing genomic variation in GSC protein-coding genes in the HS patient cohort.

A comprehensive in silico exploration of the impacts of missense variants on two different conformations of human pirin protein

BACKGROUND

Pirin, a member of the cupin superfamily, is an iron-binding non-heme protein. It acts as a coregulator of several transcription factors, especially the members of NFκB transcription factor family. Based on the redox state of its iron cofactor, it can assume two different conformations and thereby act as a redox sensor inside the nucleus. Previous studies suggested that pirin may be associated with cancer, inflammatory diseases as well as COVID-19 severities. Hence, it is important to explore the pathogenicity of its missense variants. In this study, we used a number of in silico tools to investigate the effects of missense variants of pirin on its structure, stability, metal cofactor binding affinity and interactions with partner proteins. In addition, we used protein dynamics simulation to elucidate the effects of selected variants on its dynamics. Furthermore, we calculated the frequencies of haplotypes containing pirin missense variants across five major super-populations (African, Admixed American, East Asian, European and South Asian).

RESULTS

Among a total of 153 missense variants of pirin, 45 were uniformly predicted to be pathogenic. Of these, seven variants can be considered for further experimental studies. Variants R59P and L116P were predicted to significantly destabilize and damage pirin structure, substantially reduce its affinity to its binding partners and alter pirin residue fluctuation profile via changing the flexibility of several key residues. Additionally, variants R59Q, F78V, G98D, V151D and L220P were found to impact pirin structure and function in multiple ways. As no haplotype was identified to be harboring more than one missense variant, further interrogation of the individual effects of these seven missense variants is highly recommended.

CONCLUSIONS

Pirin is involved in the transcriptional regulation of several genes and can play an important role in inflammatory responses. The variants predicted to be pathogenic in this study may thus contribute to a better understanding of the underlying molecular mechanisms of various inflammatory diseases. Future studies should be focused on clarifying if any of these variants can be used as disease biomarkers.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1186/s42269-022-00917-7.

Assessing the clinical utility of protein structural analysis in genomic variant classification: experiences from a diagnostic laboratory

BACKGROUND

The widespread clinical application of genome-wide sequencing has resulted in many new diagnoses for rare genetic conditions, but testing regularly identifies variants of uncertain significance (VUS). The remarkable rise in the amount of genomic data has been paralleled by a rise in the number of protein structures that are now publicly available, which may have clinical utility for the interpretation of missense and in-frame insertions or deletions.

METHODS

Within a UK National Health Service genomic medicine diagnostic laboratory, we investigated the number of VUS over a 5-year period that were evaluated using protein structural analysis and how often this analysis aided variant classification.

RESULTS

We found 99 novel missense and in-frame variants across 67 genes that were initially classified as VUS by our diagnostic laboratory using standard variant classification guidelines and for which further analysis of protein structure was requested. Evidence from protein structural analysis was used in the re-assessment of 64 variants, of which 47 were subsequently reclassified as pathogenic or likely pathogenic and 17 remained as VUS. We identified several case studies where protein structural analysis aided variant interpretation by predicting disease mechanisms that were consistent with the observed phenotypes, including loss-of-function through thermodynamic destabilisation or disruption of ligand binding, and gain-of-function through de-repression or escape from proteasomal degradation.

CONCLUSIONS

We have shown that using in silico protein structural analysis can aid classification of VUS and give insights into the mechanisms of pathogenicity. Based on our experience, we propose a generic evidence-based workflow for incorporating protein structural information into diagnostic practice to facilitate variant classification.

Pathogenicity of missense variants affecting the collagen IV α5 carboxy non-collagenous domain in X-linked Alport syndrome

X-linked Alport syndrome is a genetic kidney disease caused by pathogenic COL4A5 variants, but little is known of the consequences of missense variants affecting the NC1 domain of the corresponding collagen IV α5 chain. This study examined these variants in a normal (gnomAD) and other databases (LOVD, Clin Var and 100,000 Genomes Project) to determine their pathogenicity and clinical significance. Males with Cys substitutions in the collagen IV α5 NC1 domain reported in LOVD (n = 25) were examined for typical Alport features, including age at kidney failure. All NC1 variants in LOVD (n = 86) were then assessed for structural damage using an online computational tool, Missense3D. Variants in the ClinVar, gnomAD and 100,000 Genomes Project databases were also examined for structural effects. Predicted damage associated with NC1 substitutions was then correlated with the level of conservation of the affected residues. Cys substitutions in males were associated with the typical features of X-linked Alport syndrome, with a median age at kidney failure of 31 years. NC1 substitutions predicted to cause structural damage were overrepresented in LOVD (p < 0.001), and those affecting Cys residues or 'buried' Gly residues were more common than expected (both p < 0.001). Most NC1 substitutions in gnomAD (88%) were predicted to be structurally-neutral. Substitutions affecting conserved residues resulted in more structural damage than those affecting non-conserved residues (p < 0.001). Many pathogenic missense variants affecting the collagen IV α5 NC1 domain have their effect through molecular structural damage and 3D modelling is a useful tool in their assessment.

Complex effects on CaV2.1 channel gating caused by a CACNA1A variant associated with a severe neurodevelopmental disorder

P/Q-type Ca²⁺ currents mediated by Ca_V2.1 channels are essential for active neurotransmitter release at neuromuscular junctions and many central synapses. Mutations in CACNA1A, the gene encoding the principal Ca_V2.1 α_1A subunit, cause a broad spectrum of neurological disorders. Typically, gain-of-function (GOF) mutations are associated with migraine and epilepsy while loss-of-function (LOF) mutations are causative for episodic and congenital ataxias. However, a cluster of severe Ca_V2.1 channelopathies have overlapping presentations which suggests that channel dysfunction in these disorders cannot always be defined bimodally as GOF or LOF. In particular, the R1667P mutation causes focal seizures, generalized hypotonia, dysarthria, congenital ataxia and, in one case, cerebral edema leading ultimately to death. Here, we demonstrate that the R1667P mutation causes both channel GOF (hyperpolarizing voltage-dependence of activation, slowed deactivation) and LOF (slowed activation kinetics) when expressed heterologously in tsA-201 cells. We also observed a substantial reduction in Ca²⁺ current density in this heterologous system. These changes in channel gating and availability/expression manifested in diminished Ca²⁺ flux during action potential-like stimuli. However, the integrated Ca²⁺ fluxes were no different when normalized to tail current amplitude measured upon repolarization from the reversal potential. In summary, our findings indicate a complex functional effect of R1667P and support the idea that pathological missense mutations in Ca_V2.1 may not represent exclusively GOF or LOF.

Cochlear nerve deficiency in SOX11-related Coffin-Siris syndrome

The phenotypic spectrum of SOX11-related Coffin-Siris syndrome (CSS) is expanding with reports of new associations. SOX11 is implicated in neurogenesis and inner ear development. Cochlear nerve deficiency, absence or hypoplasia, is commonly associated with cochlear canal stenosis or with CHARGE syndrome, a monogenic condition that affects inner ear development. SOX11 is a transcription factor essential for neuronal identity, highly correlated with the expression of CHD7, which regulates SOX11. We present two unrelated probands, each with novel de novo SOX11 likely pathogenic variants and phenotypic manifestations of CSS including global developmental delay, growth deficiency, and hypoplastic nails. They have unilateral sensorineural hearing loss due to cochlear nerve deficiency confirmed on MRI. SOX11 is implicated in sensory neuron survival and maturation. It is highly expressed in the developing inner ear. Homozygous ablation of SOX11 in a mouse model resulted in a reduction in sensory neuron survival and decreased axonal growth. A heterozygous knockout mice model had hearing impairment with grossly normal inner ear structures like the two probands reported. We propose cochlear nerve deficiency as a new phenotypic feature of SOX11-related CSS. Magnetic resonance imaging is useful in delineating the cochlear nerve deficiency and other CSS-related brain malformations.

pyFoldX: enabling biomolecular analysis and engineering along structural ensembles

SUMMARY

Recent years have seen an increase in the number of structures available, not only for new proteins but also for the same protein crystallized with different molecules and proteins. While protein design software has proven to be successful in designing and modifying proteins, they can also be overly sensitive to small conformational differences between structures of the same protein. To cope with this, we introduce here pyFoldX, a python library that allows the integrative analysis of structures of the same protein using FoldX, an established forcefield and modelling software. The library offers new functionalities for handling different structures of the same protein, an improved molecular parametrization module and an easy integration with the data analysis ecosystem of the python programming language.

AVAILABILITY AND IMPLEMENTATION

pyFoldX rely on the FoldX software for energy calculations and modelling, which can be downloaded upon registration in http://foldxsuite.crg.eu/ and its licence is free of charge for academics. The pyFoldX library is open-source. Full details on installation, tutorials covering the library functionality and the scripts used to generate the data and figures presented in this paper are available at https://github.com/leandroradusky/pyFoldX.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

PDBe-KB consortium, Varadi M, Anyango S, Armstrong D, Berrisford J, Choudhary P, Deshpande M, Nadzirin N, Nair SS, Pravda L, Tanweer A, Al-Lazikani B, Andreini C, Barton GJ, Bednar D, Berka K, Blundell T, Brock KP, Carazo JM, Damborsky J, David A, Dey S, Dunbrack R, Recio JF, Fraternali F, Gibson T, Helmer-Citterich M, Hoksza D, Hopf T, Jakubec D, Kannan N, Krivak R, Kumar M, Levy ED, London N, Macias JR, Srivatsan MM, Marks DS, Martens L, McGowan SA, McGreig JE, Modi V, Parra RG, Pepe G, Piovesan D, Prilusky J, Putignano V, Radusky LG, Ramasamy P, Rausch AO, Reuter N, Rodriguez LA, Rollins NJ, Rosato A, Rubach P, Serrano L, Singh G, Skoda P, Sorzano COS, Stourac J, Sulkowska JI, Svobodova R, Tichshenko N, Tosatto SCE, Vranken W, Wass MN, Xue D, Zaidman D, Thornton J, Sternberg M, Orengo C, Velankar S. PDBe-KB: collaboratively defining the biological context of structural data

The Protein Data Bank in Europe - Knowledge Base (PDBe-KB, https://pdbe-kb.org) is an open collaboration between world-leading specialist data resources contributing functional and biophysical annotations derived from or relevant to the Protein Data Bank (PDB). The goal of PDBe-KB is to place macromolecular structure data in their biological context by developing standardised data exchange formats and integrating functional annotations from the contributing partner resources into a knowledge graph that can provide valuable biological insights. Since we described PDBe-KB in 2019, there have been significant improvements in the variety of available annotation data sets and user functionality. Here, we provide an overview of the consortium, highlighting the addition of annotations such as predicted covalent binders, phosphorylation sites, effects of mutations on the protein structure and energetic local frustration. In addition, we describe a library of reusable web-based visualisation components and introduce new features such as a bulk download data service and a novel superposition service that generates clusters of superposed protein chains weekly for the whole PDB archive.

Metabolovigilance: Associating Drug Metabolites with Adverse Drug Reactions

The Metabolovigilance database (https://pharmacogenomics.clas.ucdenver.edu/pharmacogenomics/side-effect/) is a single repository of information on over 15,920 pharmaceuticals and the compounds expected to result from metabolism of these drugs. Metabolovigilance functions as both a web server, providing data directly to users and as a web application, applying user inputs to create logic statements that curate the data presented or downloaded. Using this tool, it is easy to collect information on drugs, their side effects, and the metabolites associated with specific side effects. Information on these compounds can be sorted based on physical properties of the drugs and their metabolites. All of this information can be viewed, sorted, and downloaded for use in other applications. This open-access tool will facilitate molecular studies on the causes of adverse drug reactions and is well suited to integrate with genomic data furthering the goals of personalized medicine.

Structural model of human PORCN illuminates disease-associated variants and drug-binding sites

Wnt signaling is essential for normal development and is a therapeutic target in cancer. The enzyme PORCN, or porcupine, is a membrane-bound O-acyltransferase (MBOAT) that is required for the post-translational modification of all Wnts, adding an essential mono-unsaturated palmitoleic acid to a serine on the tip of Wnt hairpin 2. Inherited mutations in PORCN cause focal dermal hypoplasia, and therapeutic inhibition of PORCN slows the growth of Wnt-dependent cancers. Based on homology to mammalian MBOAT proteins, we developed and validated a structural model of human PORCN. The model accommodates palmitoleoyl-CoA and Wnt hairpin 2 in two tunnels in the conserved catalytic core, shedding light on the catalytic mechanism. The model predicts how previously uncharacterized human variants of uncertain significance can alter PORCN function. Drugs including ETC-159, IWP-L6 and LGK-974 dock in the PORCN catalytic site, providing insights into PORCN pharmacologic inhibition. This structural model enhances our mechanistic understanding of PORCN substrate recognition and catalysis, as well as the inhibition of its enzymatic activity, and can facilitate the development of improved inhibitors and the understanding of disease-relevant PORCN mutants. This article has an associated First Person interview with the joint first authors of the paper.

Comprehensive Study of Human FBXW7 Deleterious nsSNP's Functional Inference and Susceptibility to Gynaecological Cancer

Cancer is one of the world's major causes of mortality, and it plays a most important role in the world's declining life expectancy. F-box and WD-40 domain protein 7 (FBXW7), a typical participant of the F-box family of proteins, has been considered as an anti-tumor protein and one of the maximum deregulated ubiquitin-proteasome system proteins in uterine carcinosarcoma, endometrial clear cell carcinoma and cervical carcinoma with the greatest prevalence of alterations. FBXW7 variants with known clinical significance, as well as non-synonymous single nucleotide polymorphisms (nsSNPs) in the F-Box and WD40 domains, were evaluated using functionality prediction web resources. Upon analysing the seventy-three deleterious nsSNP's impact on protein stability and function, we identified that forty-one nsSNPs of WD40 domain and three of F-Box domain imply decreased stability of the FBXW7 structure. Next to TP53 and PTEN, FBXW7 was reported with the highest percentage of arginine substitution among mutations related to cancer. The current research concentrated on two arginine residue locations (Arg465, Arg505) within the WD40-repeat domain, which is vital for substrate binding. Computational analysis revealed significant deviation in stability and structural configuration of mutants R505L, R465H, R465P, R505G, R505C, R465C, R505S and R505L structures. Protein-protein interaction network of FBXW7 populated with promising hub proteins NOTCH1, c-Myc, CCNE1, STYX, KLG5, SREB1, NFKB2, SKP1 and CUL1; thus, alteration in the FBXW7 leads to aberration in their signalling pathways as well as their substrate binding ability makes this protein as attractive target for personalized therapeutic intervention.

The AlphaFold Database of Protein Structures: A Biologist's Guide

AlphaFold, the deep learning algorithm developed by DeepMind, recently released the three-dimensional models of the whole human proteome to the scientific community. Here we discuss the advantages, limitations and the still unsolved challenges of the AlphaFold models from the perspective of a biologist, who may not be an expert in structural biology.