Coevolution at the gene level, as reflected by correlated events of gene loss or gain, can be revealed by phylogenetic profile analysis. The optimal method and metric for comparing phylogenetic profiles, especially in eukaryotic genomes, are not yet established. Here, we describe a procedure suitable for large-scale analysis, which can reveal coevolution based on the assessment of the statistical significance of correlated presence/absence transitions between gene pairs. This metric can identify coevolution in profiles with low overall similarities and is not affected by similarities lacking coevolutionary information. We applied the procedure to a large collection of 60,912 orthologous gene groups (orthogroups) in 1,264 eukaryotic genomes extracted from OrthoDB. We found significant cotransition scores for 7,825 orthogroups associated in 2,401 coevolving modules linking known and unknown genes in protein complexes and biological pathways. To demonstrate the ability of the method to predict hidden gene associations, we validated through experiments the involvement of vertebrate malate synthase-like genes in the conversion of ( S )-ureidoglycolate into glyoxylate and urea, the last step of purine catabolism. This identification explains the presence of glyoxylate cycle genes in metazoa and suggests an anaplerotic role of purine degradation in early eukaryotes.

In cystic fibrosis (CF), the accumulation of viscous lung secretions rich in DNA and actin is a major cause of chronic inflammation and recurrent infections leading to airway obstruction. Mucolytic therapy based on recombinant human DNase1 reduces CF mucus viscosity and promotes airway clearance. However, the marked susceptibility to actin inhibition of this enzyme prompts the research of alternative treatments that could overcome this limitation. Within the human DNase repertoire, DNase1L2 is ideally suited for this purpose because it exhibits metal-dependent endonuclease activity on plasmid DNA in a broad range of pH with acidic optimum and is minimally inhibited by actin. When tested on CF artificial mucus enriched with actin, submicromolar concentrations of DNase1L2 reduces mucus viscosity by 50% in a few seconds. Inspection of superimposed model structures of DNase1 and DNase1L2 highlights differences at the actin-binding interface that justify the increased resistance of DNase1L2 toward actin inhibition. Furthermore, a PEGylated form of the enzyme with preserved enzymatic activity was obtained, showing interesting results in terms of activity. This work represents an effort toward the exploitation of natural DNase variants as promising alternatives to DNase1 for the treatment of CF lung disease.

The cancer epigenome has been studied in cells cultured in two-dimensional (2D) monolayers, but recent studies highlight the impact of the extracellular matrix and the three-dimensional (3D) environment on multiple cellular functions. Here, we report the physical, biochemical, and genomic differences between T47D breast cancer cells cultured in 2D and as 3D spheroids. Cells within 3D spheroids exhibit a rounder nucleus with less accessible, more compacted chromatin, as well as altered expression of ~2000 genes, the majority of which become repressed. Hi-C analysis reveals that cells in 3D are enriched for regions belonging to the B compartment, have decreased chromatin-bound CTCF and increased fusion of topologically associating domains (TADs). Upregulation of the Hippo pathway in 3D spheroids results in the activation of the LATS1 kinase, which promotes phosphorylation and displacement of CTCF from DNA, thereby likely causing the observed TAD fusions. 3D cells show higher chromatin binding of progesterone receptor (PR), leading to an increase in the number of hormone-regulated genes. This effect is in part mediated by LATS1 activation, which favors cytoplasmic retention of YAP and CTCF removal.

The increasing availability of experimental and computational protein structures entices their use for function prediction. Here we develop an automated procedure to identify enzymes involved in metabolic reactions by assessing substrate conformations docked to a library of protein structures. By screening AlphaFold-modeled vitamin B6-dependent enzymes, we find that a metric based on catalytically favorable conformations at the enzyme active site performs best (AUROC Score=0.84) in identifying genes associated with known reactions. Applying this procedure, we identify the mammalian gene encoding hydroxytrimethyllysine aldolase (HTMLA), the second enzyme of carnitine biosynthesis. Upon experimental validation, we find that the top-ranked candidates, serine hydroxymethyl transferase (SHMT) 1 and 2, catalyze the HTMLA reaction. However, a mouse protein absent in humans (threonine aldolase; Tha1) catalyzes the reaction more efficiently. Tha1 did not rank highest based on the AlphaFold model, but its rank improved to second place using the experimental crystal structure we determined at 2.26 Å resolution. Our findings suggest that humans have lost a gene involved in carnitine biosynthesis, with HTMLA activity of SHMT partially compensating for its function.

Excretion of L-xylulose is the hallmark of pentosuria, the fourth of Garrod’s inborn errors of metabolism, yet the molecular basis for L-xylulose formation remains unknown. Here, by projecting coevolutionary data for 511,114 orthogroups across 1,929 eukaryotic genomes onto metabolic maps, we screen for unmapped genes in human metabolism. Among these, we show that the DUF1907 domain of C11orf54 catalyzes formation of L-xylulose by establishing a zinc-coordinated Michaelis complex with β-keto-L-gulonate (BKG). The identification of BKG decarboxylase completes the pentose pathway, in which pentose sugars are produced by decarboxylation of nonphosphorylated hexose precursors. The pathway was present in the unicellular ancestor of animals and is conserved in all deuterostomes, in contrast to the alternative L-ascorbate (vitamin C) biosynthesis pathway. An increased flux toward pentoses may have represented an evolutionary tradeoff, favoring energy metabolism and redox cofactor balance at the expense of ascorbate biosynthesis in organisms, such as humans and other Haplorhini primates, where dietary vitamin C intake prevents scurvy.

Metabolons - transient assemblies of sequential metabolic enzymes - facilitate the reactions of multi-step metabolic pathways, yet, how they mechanistically bolster metabolic flux remains unknown. Here, we investigate the molecular determinants of metabolon formation in coenzyme Q (CoQ) biosynthesis using coarse-grained molecular dynamics simulations and biochemical experiments. We show that the COQ metabolon forms at the critical region of a phase transition, where both metabolon clustering and metabolic flux exhibit coordinated sigmoidal responses to changes in protein-protein interaction strength. These complete metabolons enable substrate channeling between sequential enzymes, leading to a crucial enhancement of CoQ production efficiency. Selectively disrupting protein-protein interactions and randomly shuffling the interaction network demonstrate that protein-proximity rather than fine structure of the metabolon clusters is imperative for substrate channeling. Grounded in both experiment and simulation, these findings provide a framework for understanding the organization and function of metabolons across diverse metabolic pathways.