Single Residue Research Articles

The dithiol mechanism of class I glutaredoxins promotes specificity for glutathione as a reducing agent

Class I glutaredoxins reversibly reduce glutathione- and nonglutathione disulfides with the help of reduced glutathione (GSH) using either a monothiol mechanism or a dithiol mechanism. The monothiol mechanism exclusively involves a single glutathionylated active-site cysteinyl residue, whereas the dithiol mechanism requires the additional formation of an intramolecular disulfide bond between the active-site cysteinyl residue and a resolving cysteinyl residue. While the oxidation of glutaredoxins by glutathione disulfide substrates has been extensively characterized, the enzyme-substrate interactions for the reduction of S-glutathionylated glutaredoxins or intramolecular glutaredoxin disulfides are still poorly characterized. Here we compared the thiol-specificity for the reduction of S-glutathionylated glutaredoxins and the intramolecular glutaredoxin disulfide. We show that S-glutathionylated glutaredoxins rapidly react with a plethora of thiols and that the 2nd glutathione-interaction site of class I glutaredoxins lacks specificity for GSH as a reducing agent. In contrast, the slower reduction of the partially strained intramolecular glutaredoxin disulfide involves specific interactions with both carboxylate groups of GSH at the 1st glutathione-interaction site. Thus, the dithiol mechanism of class I glutaredoxins promotes specificity for GSH as a reducing agent, which might explain the prevalence of dithiol glutaredoxins in pro- and eukaryotes.

Electrically silent mutants unravel the mechanism of binding-gating coupling in Cys-loop receptors.

The transduction of extracellular chemical signals into intracellular events relies on the communication between neighboring domains of membrane receptors. In the particular case of Cys-loop receptor channels, five short stretches of amino acids, one per subunit, link the extracellular and transmembrane domains in such a way that the ion permeability of the latter and the affinity for neurotransmitters of the former become tied to each other. Here, using direct functional approaches, we set out to understand the molecular bases of this crucial interdependence through the characterization of total loss-of-current mutations at the interface between domains. Our results indicate that domain-domain proximity plays a previously unnoticed critical role inasmuch as inserting a single residue in each linker rendered the two domains independent of each other. In marked contrast, loss-of-current mutations that leave the linkers' length unaltered did not compromise the interdomain coupling, but rather, seemed to cause agonist-bound closed receptors to desensitize without appreciably opening.

High fitness paths can connect proteins with low sequence overlap.

The structure and function of a protein are determined by its amino acid sequence. While random mutations change a protein's sequence, evolutionary forces shape its structural fold and biological activity. Studies have shown that neutral networks can connect a local region of sequence space by single residue mutations that preserve viability. However, the larger-scale connectedness of protein morphospace remains poorly understood. Recent advances in artificial intelligence have enabled us to computationally predict a protein's structure and quantify its functional plausibility. Here we build on these tools to develop an algorithm that generates viable paths between distantly related extant protein pairs. The intermediate sequences in these paths differ by single residue changes over subsequent steps - substitutions, insertions and deletions are admissible moves. Their fitness is evaluated using the protein language model ESM2, and maintained as high as possible subject to the constraints of the traversal. We document the qualitative variation across paths generated between progressively divergent protein pairs, some of which do not even acquire the same structural fold. The ease of interpolating between two sequences could be used as a proxy for the likelihood of homology between them.

Site-Specific Histidine Aza-Michael Addition in Proteins Enabled by a Ferritin-Based Metalloenzyme.

Histidine modifications of proteins are broadly based on chemical methods triggering N-substitution reactions such as aza-Michael addition at histidine's moderately nucleophilic imidazole side chain. While recent studies have demonstrated chemoselective, histidine-specific modifications by further exploiting imidazole's electrophilic reactivity to overcome interference from the more nucleophilic lysine and cysteine, achieving site-specific histidine modifications remains a major challenge due to the absence of spatial control over chemical processes. Herein, through X-ray crystallography and cryo-electron microscopy structural studies, we describe the rational design of a nature-inspired, noncanonical amino-acid-incorporated, human ferritin-based metalloenzyme that is capable of introducing site-specific post-translational modifications (PTMs) to histidine in peptides and proteins. Specifically, chemoenzymatic aza-Michael additions on single histidine residues were carried out on eight protein substrates ranging from 10 to 607 amino acids including the insulin peptide hormone. By introducing an insulin-targeting peptide into our metalloenzyme, we further directed modifications to be carried out site-specifically on insulin's B-chain histidine 5. The success of this biocatalysis platform outlines a novel approach in introducing residue- and, moreover, site-specific post-translational modifications to peptides and proteins, which may further enable reactions to be carried out in vivo.

Co-fermentation of dark septate endophytes and Trichoderma viride can effectively degrade cellulose in different Chinese medicinal herbal residues

Chinese medicinal herbal residues (CMHRs) are currently a widely overlooked renewable biomass resource, and developing controllable microbial fermentation technology for CMHRs will help achieve the resource utilization of CMHRs. In this study, the activity of various cellulases, reducing sugar content and cellulose degradation rate were combined to evaluate the solid state fermentation (SSF) potential of dark septate endophytes (DSE) and Trichoderma viride in separate fermentation and co-fermentation of different CMHRs (single Astragalus mongolicus residue and compound residues of Panax notoginseng and Salvia miltiorrhiza). The results showed that the cellulase activity of some DSE strains was higher than that of T. viride. T. viride fermentation alone could increase the content of substrate reducing sugars, while DSE fermentation could reduce the content of substrate reducing sugars. DSE strains such as Macrophomina pseudophaseoline, Paraboeremia selaginellae, Paraphoma chlamydocopiosa and Paraphoma radicina had the potential to co-ferment with T. viride. Optimization of the fermentation process by inoculation of T. viride after 3 d or 6 d of DSE effectively increased the cellulase activity of co-fermented CMHRs. After optimization, the activities of filter paper cellulase, carboxymethyl cellulase, and β-glucosidase in different fermentation combinations reached 2.51–5.61, 2.86–4.72 and 12.85–32.35 U/g, respectively. Further, the co-fermentation of DSE and T. viride effectively degraded cellulose in different CMHRs with a degradation rate of 54.43%–67.38%. Scanning electron microscope confirmed that the co-fermentation of DSE and T. viride could increase the external surface area and porosity of CMHRs, and destroy the structure of CMHRs. The results of this study both provided a basis for the microbial degradation of single and compound CMHR, which is of great significance for solving the treatment of waste such as Chinese herbal residue and further resource utilization.

An Ala/Glu difference in E1 of Cx26 and Cx30 contributes to their differential anionic permeabilities.

Two closely related connexins, Cx26 and Cx30, share widespread expression in the cochlear cellular networks. Gap junction channels formed by these connexins have been shown to have different permeability profiles, with Cx30 showing a strongly reduced preference for anionic tracers. The pore-forming segment of the first extracellular loop, E1, identified by computational studies of the Cx26 crystal structure to form a parahelix and a narrowed region of the pore, differs at a single residue at position 49. Cx26 contains an Ala and Cx30, a charged Glu at this position, and cysteine scanning in hemichannels identified this position to be pore-lining. To assess whether the Ala/Glu difference affects permeability, we modeled and quantified Lucifer Yellow transfer between HeLa cell pairs expressing WT Cx26 and Cx30 and variants that reciprocally substituted Glu and Ala at position 49. Cx26(A49E) and Cx30(E49A) substitutions essentially reversed the Lucifer Yellow permeability profile when accounting for junctional conductance. Moreover, by using a calcein efflux assay in single cells, we observed a similar reduced anionic preference in undocked Cx30 hemichannels and a reversal with reciprocal Ala/Glu substitutions. Thus, our data indicate that Cx26 and Cx30 gap junction channels and undocked hemichannels retain similar permeability characteristics and that a single residue difference in their E1 domains can largely account for their differential permeabilities to anionic tracers. The higher anionic permeability of Cx26 compared with Cx30 suggests that these connexins may serve distinct signaling functions in the cochlea, perhaps reflected in the vastly higher prevalence of Cx26 mutations in human deafness.

A single phosphorylatable amino acid residue is essential forthe recognition of multiple potyviral HCPro effectors by potato Nytbr.

Potato virus Y (PVY, Potyviridae) is among the most important viral pathogens of potato. The potato resistance gene Nytbr confers hypersensitive resistance to the ordinary strain of PVY (PVYO), but not the necrotic strain (PVYN). Here, we unveil that residue 247 of PVY helper component proteinase (HCPro) acts as a central player controlling Nytbr strain-specific activation. We found that substituting the serine at 247 in the HCPro of PVYO (HCProO) with an alanine as in PVYN HCPro (HCProN) disrupts Nytbr recognition. Conversely, an HCProN mutant carrying a serine at position 247 triggers defence. Moreover, we demonstrate that plant defences are induced against HCProO mutants with a phosphomimetic or another phosphorylatable residue at 247, but not with a phosphoablative residue, suggesting that phosphorylation could modulate Nytbr resistance. Extending beyond PVY, we establish that the same response elicited by the PVYO HCPro is also induced by HCPro proteins from other members of the Potyviridae family that have a serine at position 247, but not by those with an alanine. Together, our results provide further insights in the strain-specific PVY resistance in potato and infer a broad-spectrum detection mechanism of plant potyvirus effectors contingent on a single amino acid residue.

Phosphorylation of Piezo1 at a single residue, serine-1612, regulates its mechanosensitivity and invivo mechanotransduction function.

Piezo1 is a mechanically activated cation channel that converts mechanical force into diverse physiological processes. Owing to its large protein size of more than 2,500 amino acids and complex 38-transmembrane helix topology, how Piezo1 is post-translationally modified for regulating its invivo mechanotransduction functions remains largely unexplored. Here, we show that PKA activation potentiates the mechanosensitivity and slows the inactivation kinetics of mouse Piezo1 and identify the major phosphorylation site, serine-1612 (S1612), that also responds to PKC activation and shear stress. Mutating S1612 abolishes PKA and PKC regulation of Piezo1 activities. Primary endothelial cells derived from the Piezo1-S1612A knockin mice lost PKA- and PKC-dependent phosphorylation and functional potentiation of Piezo1. The mutant mice show activity-dependent elevation of blood pressure and compromised exercise endurance, resembling endothelial-specific Piezo1 knockout mice. Taken together, we identify the major PKA and PKC phosphorylation site in Piezo1 and demonstrate its contribution to Piezo1-mediated physiological functions.

Cohesive Living Bacterial Films with Tunable Mechanical Properties from Cell Surface Protein Display.

Engineered living materials (ELMs) constitute a novel class of functional materials that contain living organisms. The mechanical properties of many such systems are dominated by the polymeric matrices used to encapsulate the cellular components of the material, making it hard to tune the mechanical behavior through genetic manipulation. To address this issue, we have developed living materials in which mechanical properties are controlled by the cell-surface display of engineered proteins. Here, we show that engineered Esherichia coli cells outfitted with surface-displayed elastin-like proteins (ELPs, designated E6) grow into soft, cohesive bacterial films with biaxial moduli around 14 kPa. When subjected to bulge-testing, such films yielded at strains of approximately 10%. Introduction of a single cysteine residue near the exposed N-terminus of the ELP (to afford a protein designated CE6) increases the film modulus 3-fold to 44 kPa and eliminates the yielding behavior. When subjected to oscillatory stress, films prepared from E. coli strains bearing CE6 exhibit modest hysteresis and full strain recovery; in E6 films much more significant hysteresis and substantial plastic deformation are observed. CE6 films heal autonomously after damage, with the biaxial modulus fully restored after a few hours. This work establishes an approach to living materials with genetically programmable mechanical properties and a capacity for self-healing. Such materials may find application in biomanufacturing, biosensing, and bioremediation.

Short Fragmented Peptides from Pardachirus Marmoratus Exhibit Stronger Anticancer Activities in In Silico Residue Replacement and Analyses.

Cancer is a worldwide issue. It has been observed that conventional therapies face many problems, such as side effects and drug resistance. Recent research reportedly used marine-derived products to treat various diseases and explored their potential in treating cancers. This study aims to discover short-length anticancer peptides derived from pardaxin 6 through an in silico approach. Fragmented peptides ranging from 5 to 15 amino acids were derived from the pardaxin 6 parental peptide. These peptides were further replaced with one residue and, along with the original fragmented peptides, were predicted for their SVM scores and physicochemical properties. The top 5 derivative peptides were further examined for their toxicity, hemolytic probability, peptide structures, docking models, and energy scores using various web servers. The trend of in silico analysis outputs across 5 to 15 amino acid fragments was further analyzed. Results showed that when the amino acids were increased, SVM scores of the original fragmented peptides were also increased. Designed peptides had increased SVM scores, which was aligned with previous studies where the single residue replacement transformed the non-anticancer peptide into an anticancer agent. Moreover, in vitro studies validated that the designed peptides retained or enhanced anticancer effects against different cancer cell lines. Interestingly, a decreasing trend was observed in those fragmented derivative peptides. Single residue replacement in fragmented pardaxin 6 was found to produce stronger anticancer agents through in silico predictions. Through bioinformatics tools, fragmented peptides improved the efficiency of marine-derived drugs with higher efficacy and lower hemolytic effects in treating cancers.

Engineered Nanobodies Bind Bismuth, Indium and Gallium for Applications in Theranostics.

Targeted theranostics heavily rely on metal isotopes conjugated to antibodies. Single-domain antibodies, known as nanobodies, are much smaller in size without compromising specificity and affinity. The conventional way of conjugating metals to nanobodies involves non-specific modification of amino acid residues with bifunctional chelating agents. We demonstrate that mutagenesis of a single residue in a nanobody creates a triple cysteine motif that selectively binds bismuth which is, for example, used in targeted alpha therapy. Two mutations create a quadruple cysteine mutant specific for gallium and indium used in positron emission tomography and single-photon emission computed tomography, respectively. Labelling is quantitative within a few minutes. The metal nanobodies maintain structural integrity and stability over weeks, resist competition from endogenous metal binders like glutathione, and retain functionality.

Computational Tools for Hydrogen-Deuterium Exchange Mass Spectrometry Data Analysis.

Hydrogen-deuterium exchange (HDX) has become a pivotal method for investigating the structural and dynamic properties of proteins. The versatility and sensitivity of mass spectrometry (MS) made the technique the ideal companion for HDX, and today HDX-MS is addressing a growing number of applications in both academic research and industrial settings. The prolific generation of experimental data has spurred the concurrent development of numerous computational tools, designed to automate parts of the workflow while employing different strategies to achieve common objectives. Various computational methods are available to perform automated peptide searches and identification; different statistical tests have been implemented to quantify differences in the exchange pattern between two or more experimental conditions; alternative strategies have been developed to deconvolve and analyze peptides showing multimodal behavior; and different algorithms have been proposed to computationally increase the resolution of HDX-MS data, with the ultimate aim to provide information at the level of the single residue. This review delves into a comprehensive examination of the merits and drawbacks associated with the diverse strategies implemented by software tools for the analysis of HDX-MS data.

Engineered Nanobodies Bind Bismuth, Indium and Gallium for Applications in Theranostics

Targeted theranostics heavily rely on metal isotopes conjugated to antibodies. Single‐domain antibodies, known as nanobodies, are much smaller in size without compromising specificity and affinity. The conventional way of conjugating metals to nanobodies involves non‐specific modification of amino acid residues with bifunctional chelating agents. We demonstrate that mutagenesis of a single residue in a nanobody creates a triple cysteine motif that selectively binds bismuth which is, for example, used in targeted alpha therapy. Two mutations create a quadruple cysteine mutant specific for gallium and indium used in positron emission tomography and single‐photon emission computed tomography, respectively. Labelling is quantitative within a few minutes. The metal nanobodies maintain structural integrity and stability over weeks, resist competition from endogenous metal binders like glutathione, and retain functionality.

Diversification of sphingolipid synthase activities in kinetoplastid protozoa

Phosphosphingolipids (PSL) are essential components of eukaryotic membranes. The major PSL in fungi and protists is inositol phosphorylceramide (IPC), while sphingomyelin (SM), and to a lesser extent ethanolamine phosphorylceramide (EPC) predominate in mammals. Most kinetoplastid protozoa have a syntenic locus that encodes a single sphingolipid synthase (SLS) gene. Uniquely, among the kinetoplastids, the salivarian (African) trypanosomes have expanded this locus from a single gene in Trypanosoma vivax (TvSLS) to four genes in T. brucei (TbSLS1-4). We have previously shown that one of these is an IPC synthase, while the others are SM/EPC synthases, and that specificity is controlled by a single signature residue (IPC, serine; SM/EPC, phenylalanine). This residue is serine in T. cruzi and Leishmania major SLSs, both of which are demonstrated IPC synthases. However, T. vivax has a tyrosine at this residue raising the issue of specificity. Using a liposome-supplemented in vitro translation system we now show that T. vivax SLS is an SM/EPC synthase, and that the basal kinetoplastid Bodo saltans SLS is an IPC synthase (serine). We use these data, and a multiple alignment of available sequences, to discuss the evolution of kinetoplastid SLSs and their unique expansion in T. brucei and related salivarian trypanosomes.

PreMLS: The undersampling technique based on ClusterCentroids to predict multiple lysine sites.

The translated protein undergoes a specific modification process, which involves the formation of covalent bonds on lysine residues and the attachment of small chemical moieties. The protein's fundamental physicochemical properties undergo a significant alteration. The change significantly alters the proteins' 3D structure and activity, enabling them to modulate key physiological processes. The modulation encompasses inhibiting cancer cell growth, delaying ovarian aging, regulating metabolic diseases, and ameliorating depression. Consequently, the identification and comprehension of post-translational lysine modifications hold substantial value in the realms of biological research and drug development. Post-translational modifications (PTMs) at lysine (K) sites are among the most common protein modifications. However, research on K-PTMs has been largely centered on identifying individual modification types, with a relative scarcity of balanced data analysis techniques. In this study, a classification system is developed for the prediction of concurrent multiple modifications at a single lysine residue. Initially, a well-established multi-label position-specific triad amino acid propensity algorithm is utilized for feature encoding. Subsequently, PreMLS: a novel ClusterCentroids undersampling algorithm based on MiniBatchKmeans was introduced to eliminate redundant or similar major class samples, thereby mitigating the issue of class imbalance. A convolutional neural network architecture was specifically constructed for the analysis of biological sequences to predict multiple lysine modification sites. The model, evaluated through five-fold cross-validation and independent testing, was found to significantly outperform existing models such as iMul-kSite and predML-Site. The results presented here aid in prioritizing potential lysine modification sites, facilitating subsequent biological assays and advancing pharmaceutical research. To enhance accessibility, an open-access predictive script has been crafted for the multi-label predictive model developed in this study.

The 2024 IPD-MHC database update: a comprehensive resource for major histocompatibility complex studies.

The IPD-MHC Database project (http://www.ebi.ac.uk/ipd/mhc/) serves as a comprehensive and expertly curated repository for major histocompatibility complex (MHC) sequences from non-human species, providing the necessary infrastructure and tools to study the function and evolution of this highly polymorphic genomic region. In its latest version, the IPD-MHC database has expanded both in content and in the tools for data visualization and comparison. The database now hosts over 18 000 MHC alleles from 125 species, organized into eleven taxonomic groups, all manually curated and named by the Comparative MHC Nomenclature Committee. A cetacean section has recently been included, offering researchers valuable data to study the immune system of whales, dolphins, and porpoises, as well establishing the official nomenclature platform for the Cetacea Leukocyte Antigens (CeLA). In response to user demand and reflecting broader trends in bioinformatics and immunogenetics, IPD-MHC now includes the predicted tertiary structure of over 8000 alleles and allows comparison and visualisation of allele variation within and between species at single residue resolution. These latest developments maintain the critically important link between official nomenclature of curated alleles and the ability to analyse this complex polymorphism using the most up to date methods within a single repository.

Computational Approach to Identifying New Chemical Entities as Elastase Inhibitors with Potential Antiaging Effects.

In the aging process, skin morphology might be affected by wrinkle formation due to the loss of elasticity and resilience of connective tissues linked to the cleavage of elastin by the enzymatic activity of elastase. Little information is available about the structural requirements to efficiently inhibit elastase 1 (EC 3.4.21.36) expressed in skin keratinocytes. In this study, a structure-based approach led to the identification to the pharmacophoric hypotheses that described the main structural requirements for binding to porcine pancreatic elastase as a valuable tool for the development of skin therapeutic agents due to its similarity with human elastase 1. The obtained models were subsequently refined through the application of computational alanine-scanning mutagenesis to evaluate the effect of single residues on the binding affinity and protein stability; in turn, molecular dynamic simulations were carried out; these procedures led to a simplified model bearing few essential features, enabling a reliable collection of chemical features for their interactions with elastase. Then, a virtual screening campaign on the in-house library of synthetic compounds led to the identification of a nonpeptide-based inhibitor (IC50 = 60.4 µM) belonging to the class of N-substituted-1H-benzimidazol-2-yl]thio]acetamides, which might be further exploited to obtain more efficient ligands of elastase for therapeutic applications.

Mechanistic Basis for a Single Amino Acid Residue Mutation Causing Human DNA Ligase 1 Deficiency, A Rare Pediatric Disease

In mammalian cells, DNA ligase 1 (LIG1) functions as the primary DNA ligase in both genomic replication and single-strand break repair. Several reported mutations in human LIG1, including R305Q, R641L, and R771W, cause LIG1 syndrome, a primary immunodeficiency. While the R641L and R771W mutations, respectively located in the nucleotidyl transferase and oligonucleotide binding domains, have been biochemically characterized and shown to reduce catalytic efficiency, the recently reported R305Q mutation within the DNA binding domain (DBD) remains mechanistically unexplored. The R641L and R771W mutations are known to decrease the catalytic activity of LIG1 by affecting both interdomain interactions and DNA binding during catalysis, without significantly impacting overall DNA affinity. To elucidate the molecular basis of the LIG1 syndrome-causing R305Q mutation, we purified this single-residue mutant protein and investigated its secondary structure, protein stability, DNA binding affinity, and catalytic efficiency. Our findings reveal that the R305Q mutation significantly impairs the function of LIG1 by disrupting the DBD-DNA interactions, leading to a 7–21-fold lower DNA binding affinity and a 33–300-fold reduced catalytic efficiency of LIG1. Additionally, the R305Q mutation slightly decreases LIG1’s protein stability by 2 to 3.6 °C, on par with the effect observed previously with either the R641L or R771W mutant. Collectively, our results uncover a new mechanism whereby the R305Q mutation impairs LIG1-catalyzed nicked DNA ligation, resulting in LIG1 syndrome, and highlight the crucial roles of the DBD-DNA interactions in tight DNA binding and efficient LIG1 catalysis.

A new species of alien land flatworm in the Southern United States.

Specimens of a flat and dark brown land planarian were found in a plant nursery in North Carolina, USA in 2020. On the basis of examination of photographs of the live specimens only, the specimens were considered as belonging to Obama nungara, a species originally from South America, which has now invaded a large part of Europe. Unexpectedly, a molecular analysis revealed that the specimens did not belong to this species, neither to the genus Obama. We then undertook its histological study, which finally confirmed that the species is a member of the genus Amaga: the species is herein described as a new species, Amaga pseudobama n. sp. The species has been found in three locations in North Carolina and some infested plants were from Georgia. We reinvestigated specimens collected in Florida in 2015 and found that they also belong to this species. Citizen science observations suggest its presence in other states. Therefore, it is likely that A. pseudobama has already invaded a part of south-east USA and that the invasion took place more than ten years ago. The complete 14,909 bp long mitochondrial genome was obtained. The mitogenome is colinear with those of other Geoplanidae and it was possible to find and annotate a tRNA-Thr, which has been reported missing in several geoplanids. Amaga pseudobama shares with other Geoplaninae the presence of alternative start codons in three protein-coding genes of its mitogenome. The availability of this new genome helped us to improve our annotations of the ND3 gene, for which an ATT start codon is now suggested. Also, the sequence of the ATP6 gene raised questions concerning the use of genetic code 9 to translate the protein-coding genes of Geoplanidae, as the whole translated protein would not contain a single methionine residue when using this code. Two maximum likelihood phylogenies were obtained from genomic data. The first one was based on concatenated alignments of the partial 28S, Elongation Factor 1-alpha (EF1) and cox1 genes. The second was obtained from a concatenated alignment of the mitochondrial proteins. Both strictly discriminate A. pseudobama from O. nungara and instead associate it with Amaga expatria. We note that the nine species currently accepted within Amaga can be separated into two groups, one with extrabulbar prostatic apparatus, including the type species A. amagensis, and one with intrabulbar prostatic apparatus, including the new species A. pseudobama. This suggests that species of the latter group should be separated from Amaga and constitute a new genus. This finding again illustrates the possible emergence of new invasive species in regions naturally devoid of large land planarians, such as North America. Amaga pseudobama thus deserves to be monitored in the USA, although its superficial resemblance to O. nungara and Geoplana arkalabamensis will complicate the use of photographs obtained from citizen science. Our molecular information provides tools for this monitoring.

Functional characterization of all CDKN2A missense variants and comparison to in silico models of pathogenicity.

Interpretation of variants identified during genetic testing is a significant clinical challenge. In this study, we developed a high-throughput CDKN2A functional assay and characterized all possible CDKN2A missense variants. We found that 17.7% of all missense variants were functionally deleterious. We also used our functional classifications to assess the performance of in silico models that predict the effect of variants, including recently reported models based on machine learning. Notably, we found that all in silico models performed similarly when compared to our functional classifications with accuracies of 39.5-85.4%. Furthermore, while we found that functionally deleterious variants were enriched within ankyrin repeats, we did not identify any residues where all missense variants were functionally deleterious. Our functional classifications are a resource to aid the interpretation of CDKN2A variants and have important implications for the application of variant interpretation guidelines, particularly the use of in silico models for clinical variant interpretation.