Essential and virulence-related protein interactions of pathogens revealed through deep learning.

Identification of bacterial protein-protein interactions and predicting the structures of the complexes could aid in the understanding of pathogenicity mechanisms and developing treatments for infectious diseases. Here, we developed a deep learning-based pipeline that leverages residue-residue coevolution and protein structure prediction to systematically identify and structurally characterize protein-protein interactions at the proteome-wide scale. Using this pipeline, we searched through 78 million pairs of proteins across 19 human bacterial pathogens and identified 1923 confidently predicted complexes involving essential genes and 256 involving virulence factors. Many of these complexes were not previously known; we experimentally tested 12 such predictions, and half of them were validated. The predicted interactions span core metabolic and virulence pathways ranging from post-transcriptional modification to acid neutralization to outer membrane machinery and should contribute to our understanding of the biology of these important pathogens and the design of drugs to combat them.

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than fivefold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here, we sought to define the roles of PAX1 and newly identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); p=7.07E-11, OR = 1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wild-type. By genetic targeting we found that wild-type Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, the latter suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2 or tamoxifen treatment significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a PAX1-COL11a1-MMP3 signaling axis in spinal chondrocytes.

Adolescent idiopathic scoliosis (AIS) is a twisting deformity of the spine that occurs during periods of rapid growth in children worldwide. Children with severe cases of AIS require surgery to stop it from getting worse, presenting a significant financial burden to health systems and families. Although AIS is known to cluster in families, its genetic causes and its inheritance pattern have remained elusive. Additionally, AIS is known to be more prevalent in females, a bias that has not been explained. Advances in techniques to study the genetics underlying diseases have revealed that certain variations that increase the risk of AIS affect cartilage and connective tissue. In humans, one such variation is near a gene called Pax1, and it is female-specific. The extracellular matrix is a network of proteins and other molecules in the space between cells that help connect tissues together, and it is particularly important in cartilage and other connective tissues. One of the main components of the extracellular matrix is collagen. Yu, Kanshour, Ushiki et al. hypothesized that changes in the extracellular matrix could affect the cartilage and connective tissues of the spine, leading to AIS. To show this, the scientists screened over 100,000 individuals and found that AIS is associated with variants in two genes coding for extracellular matrix proteins. One of these variants was found in a gene called Col11a1, which codes for one of the proteins that makes up collagen. To understand the relationship between Pax1 and Col11a1, Yu, Kanshour, Ushiki et al. genetically modified mice so that they would lack the Pax1 gene. In these mice, the activation of Col11a1 was reduced in the mouse spine. They also found that the form of Col11a1 associated with AIS could not suppress the activation of a gene called Mmp3 in mouse cartilage cells as effectively as unmutated Col11a1. Going one step further, the researchers found that lowering the levels of an estrogen receptor altered the activation patterns of Pax1, Col11a1, and Mmp3 in mouse cartilage cells. These findings suggest a possible mechanism for AIS, particularly in females. The findings of Yu, Kanshour, Ushiki et al. highlight that cartilage cells in the spine are particularly relevant in AIS. The results also point to specific molecules within the extracellular matrix as important for maintaining proper alignment in the spine when children are growing rapidly. This information may guide future therapies aimed at maintaining healthy spinal cells in adolescent children, particularly girls.

Computational analysis of protein-protein interactions of cancer drivers in renal cell carcinoma.

Renal cell carcinoma (RCC) is the most common type of kidney cancer with rising cases in recent years. Extensive research has identified various cancer driver proteins associated with different subtypes of RCC. Most RCC drivers are encoded by tumor suppressor genes and exhibit enrichment in functional categories such as protein degradation, chromatin remodeling, and transcription. To further our understanding of RCC, we utilized powerful deep-learning methods based on AlphaFold to predict protein-protein interactions (PPIs) involving RCC drivers. We predicted high-confidence complexes formed by various RCC drivers, including TCEB1, KMT2C/D and KDM6A of the COMPASS-related complexes, TSC1 of the MTOR pathway, and TRRAP. These predictions provide valuable structural insights into the interaction interfaces, some of which are promising targets for cancer drug design, such as the NRF2-MAFK interface. Cancer somatic missense mutations from large datasets of genome sequencing of RCCs were mapped to the interfaces of predicted and experimental structures of PPIs involving RCC drivers, and their effects on the binding affinity were evaluated. We observed more than 100 cancer somatic mutations affecting the binding affinity of complexes formed by key RCC drivers such as VHL and TCEB1. These findings emphasize the importance of these mutations in RCC pathogenesis and potentially offer new avenues for targeted therapies.

Computational analysis of regulatory regions in human protein kinases.

Eukaryotic proteins often feature modular domain structures comprising globular domains that are connected by linker regions and intrinsically disordered regions that may contain important functional motifs. The intramolecular interactions of globular domains and nonglobular regions can play critical roles in different aspects of protein function. However, studying these interactions and their regulatory roles can be challenging due to the flexibility of nonglobular regions, the long insertions separating interacting modules, and the transient nature of some interactions. Obtaining the experimental structures of multiple domains and functional regions is more difficult than determining the structures of individual globular domains. High-quality structural models generated by AlphaFold offer a unique opportunity to study intramolecular interactions in eukaryotic proteins. In this study, we systematically explored intramolecular interactions between human protein kinase domains (KDs) and potential regulatory regions, including globular domains, N- and C-terminal tails, long insertions, and distal nonglobular regions. Our analysis identified intramolecular interactions between human KDs and 35 different types of globular domains, exhibiting a variety of interaction modes that could contribute to orthosteric or allosteric regulation of kinase activity. We also identified prevalent interactions between human KDs and their flanking regions (N- and C-terminal tails). These interactions exhibit group-specific characteristics and can vary within each specific kinase group. Although long-range interactions between KDs and nonglobular regions are relatively rare, structural details of these interactions offer new insights into the regulation mechanisms of several kinases, such as HASPIN, MAPK7, MAPK15, and SIK1B.

Pathogenic mutation hotspots in protein kinase domain structure.

Control of eukaryotic cellular function is heavily reliant on the phosphorylation of proteins at specific amino acid residues, such as serine, threonine, tyrosine, and histidine. Protein kinases that are responsible for this process comprise one of the largest families of evolutionarily related proteins. Dysregulation of protein kinase signaling pathways is a frequent cause of a large variety of human diseases including cancer, autoimmune, neurodegenerative, and cardiovascular disorders. In this study, we mapped all pathogenic mutations in 497 human protein kinase domains from the ClinVar database to the reference structure of Aurora kinase A (AURKA) and grouped them by the relevance to the disease type. Our study revealed that the majority of mutation hotspots associated with cancer are situated within the catalytic and activation loops of the kinase domain, whereas non-cancer-related hotspots tend to be located outside of these regions. Additionally, we identified a hotspot at residue R371 of the AURKA structure that has the highest number of exclusively non-cancer-related pathogenic mutations (21) and has not been previously discussed.

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than five-fold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here we sought to define the roles of PAX1 and newly-identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); P=7.07e-11, OR=1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc (IVD)-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wildtype. By genetic targeting we found that wildtype Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, this suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2, or tamoxifen treatment, significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a Pax1-Col11a1-Mmp3 signaling axis in spinal chondrocytes.

Recent advances in predicting and modeling protein-protein interactions.

Protein-protein interactions (PPIs) drive biological processes, and disruption of PPIs can cause disease. With recent breakthroughs in structure prediction and a deluge of genomic sequence data, computational methods to predict PPIs and model spatial structures of protein complexes are now approaching the accuracy of experimental approaches for permanent interactions and show promise for elucidating transient interactions. As we describe here, the key to this success is rich evolutionary information deciphered from thousands of homologous sequences that coevolve in interacting partners. This covariation signal, revealed by sophisticated statistical and machine learning (ML) algorithms, predicts physiological interactions. Accurate artificial intelligence (AI)-based modeling of protein structures promises to provide accurate 3D models of PPIs at a proteome-wide scale.

AFTM: a database of transmembrane regions in the human proteome predicted by AlphaFold.

Transmembrane proteins (TMPs), with diverse cellular functions, are difficult targets for structural determination. Predictions of TMPs and the locations of transmembrane segments using computational methods could be unreliable due to the potential for false positives and false negatives and show inconsistencies across different programs. Recent advances in protein structure prediction methods have made it possible to identify TMPs and their membrane-spanning regions using high-quality structural models. We developed the AlphaFold Transmembrane proteins (AFTM) database of candidate human TMPs by identifying transmembrane regions in AlphaFold structural models of human proteins and their domains using the positioning of proteins in membranes, version 3 program, followed by automatic corrections inspired by manual analysis of the results. We compared our results to annotations from the UniProt database and the Human Transmembrane Proteome (HTP) database. While AFTM did not identify transmembrane regions in some single-pass TMPs, it identified more transmembrane regions for multipass TMPs than UniProt and HTP. AFTM also showed more consistent results with experimental structures, as benchmarked against the Protein Data Bank Transmembrane proteins (PDBTM) database. In addition, some proteins previously annotated as TMPs were suggested to be non-TMPs by AFTM. We report the results of AFTM together with those of UniProt, HTP, TmAlphaFold, PDBTM and Membranome in the online AFTM database compiled as a comprehensive resource of candidate human TMPs with structural models. Database URL http://conglab.swmed.edu/AFTM.

Classification of domains in predicted structures of the human proteome.

Recent advances in protein structure prediction have generated accurate structures of previously uncharacterized human proteins. Identifying domains in these predicted structures and classifying them into an evolutionary hierarchy can reveal biological insights. Here, we describe the detection and classification of domains from the human proteome. Our classification indicates that only 62% of residues are located in globular domains. We further classify these globular domains and observe that the majority (65%) can be classified among known folds by sequence, with a smaller fraction (33%) requiring structural data to refine the domain boundaries and/or to support their homology. A relatively small number (966 domains) cannot be confidently assigned using our automatic pipelines, thus demanding manual inspection. We classify 47,576 domains, of which only 23% have been included in experimental structures. A portion (6.3%) of these classified globular domains lack sequence-based annotation in InterPro. A quarter (23%) have not been structurally modeled by homology, and they contain 2,540 known disease-causing single amino acid variations whose pathogenesis can now be inferred using AF models. A comparison of classified domains from a series of model organisms revealed expansions of several immune response-related domains in humans and a depletion of olfactory receptors. Finally, we use this classification to expand well-known protein families of biological significance. These classifications are presented on the ECOD website (http://prodata.swmed.edu/ecod/index_human.php).

Impact of Asp/Glu-ADP-ribosylation on protein-protein interaction and protein function.

PARylation plays critical role in regulating multiple cellular processes such as DNA damage response and repair, transcription, RNA processing, and stress response. More than 300 human proteins have been found to be modified by PARylation on acidic residues, that is, Asp (D) and Glu (E). We used the deep-learning tool AlphaFold to predict protein-protein interactions (PPIs) and their interfaces for these proteins based on coevolution signals from joint multiple sequence alignments (MSAs). AlphaFold predicted 260 confident PPIs involving PARylated proteins, and about one quarter of these PPIs have D/E-PARylation sites in their predicted PPI interfaces. AlphaFold predictions offer novel insights into the mechanisms of PARylation regulations by providing structural details of the PPI interfaces. D/E-PARylation sites have a preference to occur in coil regions and disordered regions, and PPI interfaces containing D/E-PARylation sites tend to occur between short linear sequence motifs in disordered regions and globular domains. The hub protein PCNA is predicted to interact with more than 20 proteins via the common PIP box motif and the structurally variable flanking regions. D/E-PARylation sites were found in the interfaces of key components of the RNA transcription and export complex, the SF3a spliceosome complex, and H/ACA and C/D small nucleolar ribonucleoprotein complexes, suggesting that systematic PARylation have a profound effect in regulating multiple RNA-related processes such as RNA nuclear export, splicing, and modification. Finally, PARylation of SUMO2 could modulate its interaction with CHAF1A, thereby representing a potential mechanism for the cross-talk between PARylation and SUMOylation in regulation of chromatin remodeling.

Computed cancer interactome explains the effects of somatic mutations in cancers.

Protein-protein interactions (PPIs) are involved in almost all essential cellular processes. Perturbation of PPI networks plays critical roles in tumorigenesis, cancer progression, and metastasis. While numerous high-throughput experiments have produced a vast amount of data for PPIs, these data sets suffer from high false positive rates and exhibit a high degree of discrepancy. Coevolution of amino acid positions between protein pairs has proven to be useful in identifying interacting proteins and providing structural details of the interaction interfaces with the help of deep learning methods like AlphaFold (AF). In this study, we applied AF to investigate the cancer protein-protein interactome. We predicted 1,798 PPIs for cancer driver proteins involved in diverse cellular processes such as transcription regulation, signal transduction, DNA repair, and cell cycle. We modeled the spatial structures for the predicted binary protein complexes, 1,087 of which lacked previous 3D structure information. Our predictions offer novel structural insight into many cancer-related processes such as the MAP kinase cascade and Fanconi anemia pathway. We further investigated the cancer mutation landscape by mapping somatic missense mutations (SMMs) in cancer to the predicted PPI interfaces and performing enrichment and depletion analyses. Interfaces enriched or depleted with SMMs exhibit different preferences for functional categories. Interfaces enriched in mutations tend to function in pathways that are deregulated in cancers and they may help explain the molecular mechanisms of cancers in patients; interfaces lacking mutations appear to be essential for the survival of cancer cells and thus may be future targets for PPI modulating drugs.

Evolution of a chordate-specific mechanism for myoblast fusion.

Vertebrate myoblast fusion allows for multinucleated muscle fibers to compound the size and strength of mononucleated cells, but the evolution of this important process is unknown. We investigated the evolutionary origins and function of membrane-coalescing agents Myomaker and Myomixer in various groups of chordates. Here, we report that Myomaker likely arose through gene duplication in the last common ancestor of tunicates and vertebrates, while Myomixer appears to have evolved de novo in early vertebrates. Functional tests revealed a complex evolutionary history of myoblast fusion. A prevertebrate phase of muscle multinucleation driven by Myomaker was followed by the later emergence of Myomixer that enables the highly efficient fusion system of vertebrates. Evolutionary comparisons between vertebrate and nonvertebrate Myomaker revealed key structural and mechanistic insights into myoblast fusion. Thus, our findings suggest an evolutionary model of chordate fusogens and illustrate how new genes shape the emergence of novel morphogenetic traits and mechanisms.

Human mitochondrial protein complexes revealed by large-scale coevolution analysis and deep learning-based structure modeling.

MOTIVATION: Recent development of deep-learning methods has led to a breakthrough in the prediction accuracy of 3D protein structures. Extending these methods to protein pairs is expected to allow large-scale detection of protein-protein interactions (PPIs) and modeling protein complexes at the proteome level. RESULTS: We applied RoseTTAFold and AlphaFold, two of the latest deep-learning methods for structure predictions, to analyze coevolution of human proteins residing in mitochondria, an organelle of vital importance in many cellular processes including energy production, metabolism, cell death and antiviral response. Variations in mitochondrial proteins have been linked to a plethora of human diseases and genetic conditions. RoseTTAFold, with high computational speed, was used to predict the coevolution of about 95% of mitochondrial protein pairs. Top-ranked pairs were further subject to modeling of the complex structures by AlphaFold, which also produced contact probability with high precision and in many cases consistent with RoseTTAFold. Most top-ranked pairs with high contact probability were supported by known PPIs and/or similarities to experimental structural complexes. For high-scoring pairs without experimental complex structures, our coevolution analyses and structural models shed light on the details of their interfaces, including CHCHD4-AIFM1, MTERF3-TRUB2, FMC1-ATPAF2 and ECSIT-NDUFAF1. We also identified novel PPIs (PYURF-NDUFAF5, LYRM1-MTRF1L and COA8-COX10) for several proteins without experimentally characterized interaction partners, leading to predictions of their molecular functions and the biological processes they are involved in. AVAILABILITY AND IMPLEMENTATION: Data of mitochondrial proteins and their interactions are available at: http://conglab.swmed.edu/mitochondria. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Evolutionary origin and sequence signatures of the heterodimeric ABCG5/ABCG8 transporter.

ATP-binding cassette (ABC) systems, characterized by ABC-type nucleotide-binding domains (NBDs), play crucial roles in various aspects of human physiology. Human ABCG5 and ABCG8 form a heterodimeric transporter that functions in the efflux of sterols. We used sequence similarity search, multiple sequence alignment, phylogenetic analysis, and structure comparison to study the evolutionary origin and sequence signatures of ABCG5 and ABCG8. Orthologs of ABCG5 and ABCG8, supported by phylogenetic analysis and signature residues, were identified in bilaterian animals, Filasterea, Fungi, and Amoebozoa. Such a phylogenetic distribution suggests that ABCG5 and ABCG8 could have originated in the last common ancestor of Amorphea (the unikonts), the eukaryotic group including Amoebozoa and Opisthokonta. ABCG5 and ABCG8 were missing in genomes of various lineages such as snakes, jawless vertebrates, non-vertebrate chordates, echinoderms, and basal metazoan groups. Amino-acid changes in key positions in ABCG8 Walker A motif and/or ABCG5 C-loop were observed in most tetrapod organisms, likely resulted in the loss of ATPase activity at one nucleotide-binding site. ABCG5 and ABCG8 in Ecdysozoa (such as insects) exhibit elevated evolutionary rates and accumulate various changes in their NBD functional motifs. Alignment inspection revealed several residue positions that show different amino-acid usages in ABCG5/ABCG8 compared to other ABCG subfamily proteins. These residues were mapped to the structural cores of transmembrane domains (TMDs), the NBD-TMD interface, and the interface between TMDs. They serve as sequence signatures to differentiate ABCG5/ABCG8 from other ABCG subfamily proteins, and some of them may contribute to substrate specificity of the ABCG5/ABCG8 transporter.

DisEnrich: database of enriched regions in human dark proteome.

MOTIVATION: Intrinsically disordered proteins (IDPs) are involved in numerous processes crucial for living organisms. Bias in amino acid composition of these proteins determines their unique biophysical and functional features. Distinct intrinsically disordered regions (IDRs) with compositional bias play different important roles in various biological processes. IDRs enriched in particular amino acids in human proteome have not been described consistently. RESULTS: We developed DisEnrich-the database of human proteome IDRs that are significantly enriched in particular amino acids. Each human protein is described using Gene Ontology (GO) function terms, disorder prediction for the full-length sequence using three methods, enriched IDR composition and ranks of human proteins with similar enriched IDRs. Distribution analysis of enriched IDRs among broad functional categories revealed significant overrepresentation of R- and Y-enriched IDRs in metabolic and enzymatic activities and F-enriched IDRs in transport. About 75% of functional categories contain IDPs with IDRs significantly enriched in hydrophobic residues that are important for protein-protein interactions. AVAILABILITY AND IMPLEMENTATION: The database is available at http://prodata.swmed.edu/DisEnrichDB/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics Advances online.

Computed structures of core eukaryotic protein complexes.

Protein-protein interactions play critical roles in biology, but the structures of many eukaryotic protein complexes are unknown, and there are likely many interactions not yet identified. We take advantage of advances in proteome-wide amino acid coevolution analysis and deep-learning­based structure modeling to systematically identify and build accurate models of core eukaryotic protein complexes within the Saccharomyces cerevisiae proteome. We use a combination of RoseTTAFold and AlphaFold to screen through paired multiple sequence alignments for 8.3 million pairs of yeast proteins, identify 1505 likely to interact, and build structure models for 106 previously unidentified assemblies and 806 that have not been structurally characterized. These complexes, which have as many as five subunits, play roles in almost all key processes in eukaryotic cells and provide broad insights into biological function.

TMEM120A is a coenzyme A-binding membrane protein with structural similarities to ELOVL fatty acid elongase.

TMEM120A, also named as TACAN, is a novel membrane protein highly conserved in vertebrates and was recently proposed to be a mechanosensitive channel involved in sensing mechanical pain. Here we present the single-particle cryogenic electron microscopy (cryo-EM) structure of human TMEM120A, which forms a tightly packed dimer with extensive interactions mediated by the N-terminal coiled coil domain (CCD), the C-terminal transmembrane domain (TMD), and the re-entrant loop between the two domains. The TMD of each TMEM120A subunit contains six transmembrane helices (TMs) and has no clear structural feature of a channel protein. Instead, the six TMs form an α-barrel with a deep pocket where a coenzyme A (CoA) molecule is bound. Intriguingly, some structural features of TMEM120A resemble those of elongase for very long-chain fatty acids (ELOVL) despite the low sequence homology between them, pointing to the possibility that TMEM120A may function as an enzyme for fatty acid metabolism, rather than a mechanosensitive channel.

Topology evaluation of models for difficult targets in the 14th round of the critical assessment of protein structure prediction (CASP14).

This report describes the tertiary structure prediction assessment of difficult modeling targets in the 14th round of the Critical Assessment of Structure Prediction (CASP14). We implemented an official ranking scheme that used the same scores as the previous CASP topology-based assessment, but combined these scores with one that emphasized physically realistic models. The top performing AlphaFold2 group outperformed the rest of the prediction community on all but two of the difficult targets considered in this assessment. They provided high quality models for most of the targets (86% over GDT_TS 70), including larger targets above 150 residues, and they correctly predicted the topology of almost all the rest. AlphaFold2 performance was followed by two manual Baker methods, a Feig method that refined Zhang-server models, two notable automated Zhang server methods (QUARK and Zhang-server), and a Zhang manual group. Despite the remarkable progress in protein structure prediction of difficult targets, both the prediction community and AlphaFold2, to a lesser extent, faced challenges with flexible regions and obligate oligomeric assemblies. The official ranking of top-performing methods was supported by performance generated PCA and heatmap clusters that gave insight into target difficulties and the most successful state-of-the-art structure prediction methodologies.

Completeness and Consistency in Structural Domain Classifications.

Domain classifications are a useful resource for computational analysis of the protein structure, but elements of their composition are often opaque to potential users. We perform a comparative analysis of our classification ECOD against the SCOPe, SCOP2, and CATH domain classifications with respect to their constituent domain boundaries and hierarchal organization. The coverage of these domain classifications with respect to ECOD and to the PDB was assessed by structure and by sequence. We also conducted domain pair analysis to determine broad differences in hierarchy between domains shared by ECOD and other classifications. Finally, we present domains from the major facilitator superfamily (MFS) of transporter proteins and provide evidence that supports their split into domains and for multiple conformations within these families. We find that the ECOD and CATH provide the most extensive structural coverage of the PDB. ECOD and SCOPe have the most consistent domain boundary conditions, whereas CATH and SCOP2 both differ significantly.

Metabolic impact of pathogenic variants in the mitochondrial glutamyl-tRNA synthetase EARS2.

Glutamyl-tRNA synthetase 2 (encoded by EARS2) is a mitochondrial aminoacyl-tRNA synthetase required to translate the 13 subunits of the electron transport chain encoded by the mitochondrial DNA. Pathogenic EARS2 variants cause combined oxidative phosphorylation deficiency, subtype 12 (COXPD12), an autosomal recessive disorder involving lactic acidosis, intellectual disability, and other features of mitochondrial compromise. Patients with EARS2 deficiency present with variable phenotypes ranging from neonatal lethality to a mitigated disease with clinical improvement in early childhood. Here, we report a neonate homozygous for a rare pathogenic variant in EARS2 (c.949G>T; p.G317C). Metabolomics in primary fibroblasts from this patient revealed expected abnormalities in TCA cycle metabolites, as well as numerous changes in purine, pyrimidine, and fatty acid metabolism. To examine genotype-phenotype correlations in COXPD12, we compared the metabolic impact of reconstituting these fibroblasts with wild-type EARS2 versus four additional EARS2 variants from COXPD12 patients with varying clinical severity. Metabolomics identified a group of signature metabolites, mostly from the TCA cycle and amino acid metabolism, that discriminate between EARS2 variants causing relatively mild and severe COXPD12. Taken together, these findings indicate that metabolomics in patient-derived fibroblasts may help establish genotype-phenotype correlations in EARS2 deficiency and likely other mitochondrial disorders.