Genomics of Immune Diseases and New Therapies.

Genomic DNA sequencing technologies have been one of the great advances of the 21st century, having decreased in cost by seven orders of magnitude and opening up new fields of investigation throughout research and clinical medicine. Genomics coupled with biochemical investigation has allowed the molecular definition of a growing number of new genetic diseases that reveal new concepts of immune regulation. Also, defining the genetic pathogenesis of these diseases has led to improved diagnosis, prognosis, genetic counseling, and, most importantly, new therapies. We highlight the investigational journey from patient phenotype to treatment using the newly defined XMEN disease, caused by the genetic loss of the MAGT1 magnesium transporter, as an example. This disease illustrates how genomics yields new fundamental immunoregulatory insights as well as how research genomics is integrated into clinical immunology. At the end, we discuss two other recently described diseases, CHAI/LATAIE (CTLA-4 deficiency) and PASLI (PI3K dysregulation), as additional examples of the journey from unknown immunological diseases to new precision medicine treatments using genomics.

A White-Box Machine Learning Approach for Revealing Antibiotic Mechanisms of Action.

Current machine learning techniques enable robust association of biological signals with measured phenotypes, but these approaches are incapable of identifying causal relationships. Here, we develop an integrated "white-box" biochemical screening, network modeling, and machine learning approach for revealing causal mechanisms and apply this approach to understanding antibiotic efficacy. We counter-screen diverse metabolites against bactericidal antibiotics in Escherichia coli and simulate their corresponding metabolic states using a genome-scale metabolic network model. Regression of the measured screening data on model simulations reveals that purine biosynthesis participates in antibiotic lethality, which we validate experimentally. We show that antibiotic-induced adenine limitation increases ATP demand, which elevates central carbon metabolism activity and oxygen consumption, enhancing the killing effects of antibiotics. This work demonstrates how prospective network modeling can couple with machine learning to identify complex causal mechanisms underlying drug efficacy.

Mapping Causal Variants with Single-Nucleotide Resolution Reveals Biochemical Drivers of Phenotypic Change.

Understanding the sequence determinants that give rise to diversity among individuals and species is the central challenge of genetics. However, despite ever greater numbers of sequenced genomes, most genome-wide association studies cannot distinguish causal variants from linked passenger mutations spanning many genes. We report that this inherent challenge can be overcome in model organisms. By pushing the advantages of inbred crossing to its practical limit in Saccharomyces cerevisiae, we improved the statistical resolution of linkage analysis to single nucleotides. This "super-resolution" approach allowed us to map 370 causal variants across 26 quantitative traits. Missense, synonymous, and cis-regulatory mutations collectively gave rise to phenotypic diversity, providing mechanistic insight into the basis of evolutionary divergence. Our data also systematically unmasked complex genetic architectures, revealing that multiple closely linked driver mutations frequently act on the same quantitative trait. Single-nucleotide mapping thus complements traditional deletion and overexpression screening paradigms and opens new frontiers in quantitative genetics.

Reconstituted Postsynaptic Density as a Molecular Platform for Understanding Synapse Formation and Plasticity.

Synapses are semi-membraneless, protein-dense, sub-micron chemical reaction compartments responsible for signal processing in each and every neuron. Proper formation and dynamic responses to stimulations of synapses, both during development and in adult, are fundamental to functions of mammalian brains, although the molecular basis governing formation and modulation of compartmentalized synaptic assemblies is unclear. Here, we used a biochemical reconstitution approach to show that, both in solution and on supported membrane bilayers, multivalent interaction networks formed by major excitatory postsynaptic density (PSD) scaffold proteins led to formation of PSD-like assemblies via phase separation. The reconstituted PSD-like assemblies can cluster receptors, selectively concentrate enzymes, promote actin bundle formation, and expel inhibitory postsynaptic proteins. Additionally, the condensed phase PSD assemblies have features that are distinct from those in homogeneous solutions and fit for synaptic functions. Thus, we have built a molecular platform for understanding how neuronal synapses are formed and dynamically regulated.

Excitable Signal Transduction Networks in Directed Cell Migration.

Although directed migration of eukaryotic cells may have evolved to escape nutrient depletion, it has been adopted for an extensive range of physiological events during development and in the adult organism. The subversion of these movements results in disease, such as cancer. Mechanisms of propulsion and sensing are extremely diverse, but most eukaryotic cells move by extending actin-filled protrusions termed macropinosomes, pseudopodia, or lamellipodia or by extension of blebs. In addition to motility, directed migration involves polarity and directional sensing. The hundreds of gene products involved in these processes are organized into networks of parallel and interconnected pathways. Many of these components are activated or inhibited coordinately with stimulation and on each spontaneously extended protrusion. Moreover, these networks display hallmarks of excitability, including all-or-nothing responsiveness and wave propagation. Cellular protrusions result from signal transduction waves that propagate outwardly from an origin and drive cytoskeletal activity. The range of the propagating waves and hence the size of the protrusions can be altered by lowering or raising the threshold for network activation, with larger and wider protrusions favoring gliding or oscillatory behavior over amoeboid migration. Here, we evaluate the variety of models of excitable networks controlling directed migration and outline critical tests. We also discuss the utility of this emerging view in producing cell migration and in integrating the various extrinsic cues that direct migration.

Establishment of dsDNA-dsDNA interactions by the condensin complex.

Condensin is a structural maintenance of chromosomes (SMC) complex family member thought to build mitotic chromosomes by DNA loop extrusion. However, condensin variants unable to extrude loops, yet proficient in chromosome formation, were recently described. Here, we explore how condensin might alternatively build chromosomes. Using bulk biochemical and single-molecule experiments with purified fission yeast condensin, we observe that individual condensins sequentially and topologically entrap two double-stranded DNAs (dsDNAs). Condensin loading transitions through a state requiring DNA bending, as proposed for the related cohesin complex. While cohesin then favors the capture of a second single-stranded DNA (ssDNA), second dsDNA capture emerges as a defining feature of condensin. We provide complementary in vivo evidence for DNA-DNA capture in the form of condensin-dependent chromatin contacts within, as well as between, chromosomes. Our results support a "diffusion capture" model in which condensin acts in mitotic chromosome formation by sequential dsDNA-dsDNA capture.

Mechanistic insight into allosteric activation of human pyruvate carboxylase by acetyl-CoA.

Pyruvate carboxylase (PC) catalyzes the two-step carboxylation of pyruvate to produce oxaloacetate, playing a key role in the maintenance of metabolic homeostasis in cells. Given its involvement in multiple diseases, PC has been regarded as a potential therapeutic target for obesity, diabetes, and cancer. Albeit acetyl-CoA has been recognized as the allosteric regulator of PC for over 60 years, the underlying mechanism of how acetyl-CoA induces PC activation remains enigmatic. Herein, by using time-resolved cryo-electron microscopy, we have captured the snapshots of PC transitional states during its catalytic cycle. These structures and the biochemical studies reveal that acetyl-CoA stabilizes PC in a catalytically competent conformation, which triggers a cascade of events, including ATP hydrolysis and the long-distance communication between the two reactive centers. These findings provide an integrated picture for PC catalysis and unveil the unique allosteric mechanism of acetyl-CoA in an essential biochemical reaction in all kingdoms of life.

Thermophile 90S Pre-ribosome Structures Reveal the Reverse Order of Co-transcriptional 18S rRNA Subdomain Integration.

Eukaryotic ribosome biogenesis involves RNA folding and processing that depend on assembly factors and small nucleolar RNAs (snoRNAs). The 90S (SSU-processome) is the earliest pre-ribosome structurally analyzed, which was suggested to assemble stepwise along the growing pre-rRNA from 5' > 3', but this directionality may not be accurate. Here, by analyzing the structure of a series of 90S assembly intermediates from Chaetomium thermophilum, we discover a reverse order of 18S rRNA subdomain incorporation. Large parts of the 18S rRNA 3' and central domains assemble first into the 90S before the 5' domain is integrated. This final incorporation depends on a contact between a heterotrimer Enp2-Bfr2-Lcp5 recruited to the flexible 5' domain and Kre33, which reconstitutes the Kre33-Enp-Brf2-Lcp5 module on the compacted 90S. Keeping the 5' domain temporarily segregated from the 90S scaffold could provide extra time to complete the multifaceted 5' domain folding, which depends on a distinct set of snoRNAs and processing factors.

Changes in holopelagic Sargassum spp. biomass composition across an unusual year.

The year 2021 marked a decade of holopelagic sargassum (morphotypes Sargassum natans I and VIII, and Sargassum fluitans III) stranding on the Caribbean and West African coasts. Beaching of millions of tons of sargassum negatively impacts coastal ecosystems, economies, and human health. Additionally, the La Soufrière volcano erupted in St. Vincent in April 2021, at the start of the sargassum season. We investigated potential monthly variations in morphotype abundance and biomass composition of sargassum harvested in Jamaica and assessed the influence of processing methods (shade-drying vs. frozen samples) and of volcanic ash exposure on biochemical and elemental components. S. fluitans III was the most abundant morphotype across the year. Limited monthly variations were observed for key brown algal components (phlorotannins, fucoxanthin, and alginate). Shade-drying did not significantly alter the contents of proteins but affected levels of phlorotannins, fucoxanthin, mannitol, and alginate. Simulation of sargassum and volcanic ash drift combined with age statistics suggested that sargassum potentially shared the surface layer with ash for ~50 d, approximately 100 d before stranding in Jamaica. Integrated elemental analysis of volcanic ash, ambient seawater, and sargassum biomass showed that algae harvested from August had accumulated P, Al, Fe, Mn, Zn, and Ni, probably from the ash, and contained less As. This ash fingerprint confirmed the geographical origin and drift timescale of sargassum. Since environmental conditions and processing methods influence biomass composition, efforts should continue to improve understanding, forecasting, monitoring, and valorizing sargassum, particularly as strandings of sargassum show no sign of abating.

Temperature compensation through kinetic regulation in biochemical oscillators.

Nearly all circadian clocks maintain a period that is insensitive to temperature changes, a phenomenon known as temperature compensation (TC). Yet, it is unclear whether there is any common feature among different systems that exhibit TC. From a general timescale invariance, we show that TC relies on the existence of certain period-lengthening reactions wherein the period of the system increases strongly with the rates in these reactions. By studying several generic oscillator models, we show that this counterintuitive dependence is nonetheless a common feature of oscillators in the nonlinear (far-from-onset) regime where the oscillation can be separated into fast and slow phases. The increase of the period with the period-lengthening reaction rates occurs when the amplitude of the slow phase in the oscillation increases with these rates while the progression speed in the slow phase is controlled by other rates of the system. The positive dependence of the period on the period-lengthening rates balances its inverse dependence on other kinetic rates in the system, which gives rise to robust TC in a wide range of parameters. We demonstrate the existence of such period-lengthening reactions and their relevance for TC in all four model systems we considered. Theoretical results for a model of the Kai system are supported by experimental data. A study of the energy dissipation also shows that better TC performance requires higher energy consumption. Our study unveils a general mechanism by which a biochemical oscillator achieves TC by operating in parameter regimes far from the onset where period-lengthening reactions exist.

The oscillation of mitotic kinase governs cell cycle latches in mammalian cells.

The mammalian cell cycle alternates between two phases - S-G2-M with high levels of A- and B-type cyclins (CycA and CycB, respectively) bound to cyclin-dependent kinases (CDKs), and G1 with persistent degradation of CycA and CycB by an activated anaphase promoting complex/cyclosome (APC/C) bound to Cdh1 (also known as FZR1 in mammals; denoted APC/C:Cdh1). Because CDKs phosphorylate and inactivate Cdh1, these two phases are mutually exclusive. This 'toggle switch' is flipped from G1 to S by cyclin-E bound to a CDK (CycE:CDK), which is not degraded by APC/C:Cdh1, and from M to G1 by Cdc20-bound APC/C (APC/C:Cdc20), which is not inactivated by CycA:CDK or CycB:CDK. After flipping the switch, cyclin E is degraded and APC/C:Cdc20 is inactivated. Combining mathematical modelling with single-cell timelapse imaging, we show that dysregulation of CycB:CDK disrupts strict alternation of the G1-S and M-G1 switches. Inhibition of CycB:CDK results in Cdc20-independent Cdh1 'endocycles', and sustained activity of CycB:CDK drives Cdh1-independent Cdc20 endocycles. Our model provides a mechanistic explanation for how whole-genome doubling can arise, a common event in tumorigenesis that can drive tumour evolution.

Explorer: efficient DNA coding by De Bruijn graph toward arbitrary local and global biochemical constraints.

With the exponential growth of digital data, there is a pressing need for innovative storage media and techniques. DNA molecules, due to their stability, storage capacity, and density, offer a promising solution for information storage. However, DNA storage also faces numerous challenges, such as complex biochemical constraints and encoding efficiency. This paper presents Explorer, a high-efficiency DNA coding algorithm based on the De Bruijn graph, which leverages its capability to characterize local sequences. Explorer enables coding under various biochemical constraints, such as homopolymers, GC content, and undesired motifs. This paper also introduces Codeformer, a fast decoding algorithm based on the transformer architecture, to further enhance decoding efficiency. Numerical experiments indicate that, compared with other advanced algorithms, Explorer not only achieves stable encoding and decoding under various biochemical constraints but also increases the encoding efficiency and bit rate by ¿10%. Additionally, Codeformer demonstrates the ability to efficiently decode large quantities of DNA sequences. Under different parameter settings, its decoding efficiency exceeds that of traditional algorithms by more than two-fold. When Codeformer is combined with Reed-Solomon code, its decoding accuracy exceeds 99%, making it a good choice for high-speed decoding applications. These advancements are expected to contribute to the development of DNA-based data storage systems and the broader exploration of DNA as a novel information storage medium.

Defining the Transcriptional Ecosystem.

Fine tuning of the transcriptional program requires the competing action of multiple protein complexes in a well-organized environment. Genome folding creates proximity between genes, leading to accumulation of regulatory factors and formation of local microenvironments. Many roles of this complex organization controlling gene transcription remain to be explored. In this Perspective, we are proposing the existence of a transcriptional ecosystem equilibrium: a mechanism balancing transcriptional regulation between connected genes during environmental disturbances. This model is derived from chromosome architecture studies assigning genes to specific DNA structures and evidence establishing that the transcription machinery and coregulators create dynamic phase separation droplets surrounding active genes. Defining connected genes as ecosystems rather than individuals will cement that transcriptional regulation is a biochemical equilibrium and force a reassessment of direct and indirect responses to environmental disturbances.

Mechanistic analysis of a novel membrane-interacting variable loop in the pleckstrin-homology domain critical for dynamin function.

Classical dynamins are best understood for their ability to generate vesicles by membrane fission. During clathrin-mediated endocytosis (CME), dynamin is recruited to the membrane through multivalent protein and lipid interactions between its proline-rich domain (PRD) with SRC Homology 3 (SH3) domains in endocytic proteins and its pleckstrin-homology domain (PHD) with membrane lipids. Variable loops (VL) in the PHD bind lipids and partially insert into the membrane thereby anchoring the PHD to the membrane. Recent molecular dynamics (MD) simulations reveal a novel VL4 that interacts with the membrane. Importantly, a missense mutation that reduces VL4 hydrophobicity is linked to an autosomal dominant form of Charcot-Marie-Tooth (CMT) neuropathy. We analyzed the orientation and function of the VL4 to mechanistically link data from simulations with the CMT neuropathy. Structural modeling of PHDs in the cryo-electron microscopy (cryo-EM) cryoEM map of the membrane-bound dynamin polymer confirms VL4 as a membrane-interacting loop. In assays that rely solely on lipid-based membrane recruitment, VL4 mutants with reduced hydrophobicity showed an acute membrane curvature-dependent binding and a catalytic defect in fission. Remarkably, in assays that mimic a physiological multivalent lipid- and protein-based recruitment, VL4 mutants were completely defective in fission across a range of membrane curvatures. Importantly, expression of these mutants in cells inhibited CME, consistent with the autosomal dominant phenotype associated with the CMT neuropathy. Together, our results emphasize the significance of finely tuned lipid and protein interactions for efficient dynamin function.

AMP-dependent phosphite dehydrogenase, a phosphorylating enzyme in dissimilatory phosphite oxidation.

Oxidation of phosphite (HPO32-) to phosphate (HPO42-) releases electrons at a very low redox potential (E0'= -690 mV) which renders phosphite an excellent electron donor for microbial energy metabolism. To date, two pure cultures of strictly anaerobic bacteria have been isolated that run their energy metabolism on the basis of phosphite oxidation, the Gram-negative Desulfotignum phosphitoxidans (DSM 13687) and the Gram-positive Phosphitispora fastidiosa (DSM 112739). Here, we describe the key enzyme for dissimilatory phosphite oxidation in these bacteria. The enzyme catalyzed phosphite oxidation in the presence of adenosine monophosphate (AMP) to form adenosine diphosphate (ADP), with concomitant reduction of oxidized nicotinamide adenine dinucleotide (NAD+) to reduced nicotinamide adenine dinucleotide (NADH). The enzyme of P. fastidiosa was heterologously expressed in Escherichia coli. It has a molecular mass of 35.2 kDa and a high affinity for phosphite and NAD+. Its activity was enhanced more than 100-fold by addition of ADP-consuming adenylate kinase (myokinase) to a maximal activity between 30 and 80 mU x mg protein-1. A similar NAD-dependent enzyme oxidizing phosphite to phosphate with concomitant phosphorylation of AMP to ADP is found in D. phosphitoxidans, but this enzyme could not be heterologously expressed. Based on sequence analysis, these phosphite-oxidizing enzymes are related to nucleotide-diphosphate-sugar epimerases and indeed represent AMP-dependent phosphite dehydrogenases (ApdA). A reaction mechanism is proposed for this unusual type of substrate-level phosphorylation reaction.

Tissue interplay during morphogenesis.

The process by which biological systems such as cells, tissues and organisms acquire shape has been named as morphogenesis and it is central to a plethora of biological contexts including embryo development, wound healing, or even cancer. Morphogenesis relies in both self-organising properties of the system and in environmental inputs (biochemical and biophysical). The classical view of morphogenesis is based on the study of external biochemical molecules, such as morphogens. However, recent studies are establishing that the mechanical environment is also used by cells to communicate within tissues, suggesting that this mechanical crosstalk is essential to synchronise morphogenetic transitions and self-organisation. In this article we discuss how tissue interaction drive robust morphogenesis, starting from a classical biochemical view, to finalise with more recent advances on how the biophysical properties of a tissue feedback with their surroundings to allow form acquisition. We also comment on how in silico models aid to integrate and predict changes in cell and tissue behaviour. Finally, considering recent advances from the developmental biomechanics field showing that mechanical inputs work as cues that promote morphogenesis, we invite to revisit the concept of morphogen.

Analysis of Poly(ethylene terephthalate) degradation kinetics of evolved IsPETase variants using a surface crowding model.

Poly(ethylene terephthalate) (PET) is a major plastic polymer utilized in the single-use and textile industries. The discovery of PET-degrading enzymes (PETases) has led to an increased interest in the biological recycling of PET in addition to mechanical recycling. IsPETase from Ideonella sakaiensis is a candidate catalyst, but little is understood about its structure-function relationships with regards to PET degradation. To understand the effects of mutations on IsPETase productivity, we develop a directed evolution assay to identify mutations beneficial to PET film degradation at 30 °C. IsPETase also displays enzyme concentration-dependent inhibition effects, and surface crowding has been proposed as a causal phenomenon. Based on total internal reflectance fluorescence microscopy and adsorption experiments, IsPETase is likely experiencing crowded conditions on PET films. Molecular dynamics simulations of IsPETase variants reveal a decrease in active site flexibility in free enzymes and reduced probability of productive active site formation in substrate-bound enzymes under crowding. Hence, we develop a surface crowding model to analyze the biochemical effects of three hit mutations (T116P, S238N, S290P) that enhanced ambient temperature activity and/or thermostability. We find that T116P decreases susceptibility to crowding, resulting in higher PET degradation product accumulation despite no change in intrinsic catalytic rate. In conclusion, we show that a macromolecular crowding-based biochemical model can be used to analyze the effects of mutations on properties of PETases and that crowding behavior is a major property to be targeted for enzyme engineering for improved PET degradation.

Reconstitution of the alternative pathway of the complement system enables rapid delineation of the mechanism of action of novel inhibitors.

The complement system plays a critical role in the innate immune response, acting as a first line of defense against invading pathogens. However, dysregulation of the complement system is implicated in the pathogenesis of numerous diseases, ranging from Alzheimer's to age-related macular degeneration and rare blood disorders. As such, complement inhibitors have enormous potential to alleviate disease burden. While a few complement inhibitors are in clinical use, there is still a significant unmet medical need for the discovery and development of novel inhibitors to treat patients suffering from disorders of the complement system. A key hurdle in the development of complement inhibitors has been the determination of their mechanism of action. Progression along the complement cascade involves the formation of numerous multimeric protein complexes, creating the potential for inhibitors to act at multiple nodes in the pathway. This is especially true for molecules that target the central component C3 and its fragment C3b, which serve a dual role as a substrate for the C3 convertases and as a scaffolding protein in both the C3 and C5 convertases. Here, we report a step-by-step in vitro reconstitution of the complement alternative pathway using bio-layer interferometry. By physically uncoupling each step in the pathway, we were able to determine the kinetic signature of inhibitors that act at single steps in the pathway and delineate the full mechanism of action of known and novel C3 inhibitors. The method could have utility in drug discovery and further elucidating the biochemistry of the complement system.

Mapping of specialized metabolite terms onto a plant phylogeny using text mining and large language models.

Plants produce a staggering array of chemicals that are the basis for organismal function and important human nutrients and medicines. However, it is poorly defined how these compounds evolved and are distributed across the plant kingdom, hindering a systematic view and understanding of plant chemical diversity. Recent advances in plant genome/transcriptome sequencing have provided a well-defined molecular phylogeny of plants, on which the presence of diverse natural products can be mapped to systematically determine their phylogenetic distribution. Here, we built a proof-of-concept workflow where previously reported diverse tyrosine-derived plant natural products were mapped onto the plant tree of life. Plant chemical-species associations were mined from literature, filtered, evaluated through manual inspection of over 2500 scientific articles, and mapped onto the plant phylogeny. The resulting "phylochemical" map confirmed several highly lineage-specific compound class distributions, such as betalain pigments and Amaryllidaceae alkaloids. The map also highlighted several lineages enriched in dopamine-derived compounds, including the orders Caryophyllales, Liliales, and Fabales. Additionally, the application of large language models, using our manually curated data as a ground truth set, showed that post-mining processing can largely be automated with a low false-positive rate, critical for generating a reliable phylochemical map. Although a high false-negative rate remains a challenge, our study demonstrates that combining text mining with language model-based processing can generate broader phylochemical maps, which will serve as a valuable community resource to uncover key evolutionary events that underlie plant chemical diversity and enable system-level views of nature's millions of years of chemical experimentation.

Phenylalanine ammonia-lyases and 4-coumaric acid coenzyme A ligases in Chara braunii, Marchantia polymorpha, and Physcomitrium patens as extant model organisms for plant terrestrialization.

The conquest of land posed severe problems to plants which they had to cope with by adapting biosynthetic capacities. Adaptations to respond to UV irradiation, water loss, pathogen and herbivore defense, and the earth's pull were essential. Chemical compounds alleviating these problems can be synthesized by the phenylpropanoid pathway, the core of which are three enzymes: phenylalanine ammonia-lyase (PAL), cinnamic acid 4-hydroxylase, and 4-coumaric acid coenzyme A-ligase (4CL). The genomes of model organisms, Chara braunii as aquatic alga and the two bryophytes Physcomitrium patens and Marchantia polymorpha, were searched for sequences encoding PAL and 4CL and selected sequences heterologously expressed in Escherichia coli for biochemical characterization. Several possible isoforms were identified for both enzymes in Marchantia polymorpha and Physcomitrium patens, while only one or two isoforms could be retrieved for Chara braunii. Active forms of both enzymes were found in all three organisms, although the catalytic efficiencies varied in a wide range. l-Phenylalanine was accepted as best substrate by all PAL-like enzymes, despite annotations in some cases suggesting different activities. The substrate spectrum of 4CLs was more diverse, but caffeic and/or 4-coumaric acids generally were the best-accepted substrates. Our investigations show that PAL and 4CL, important enzymes for the formation of phenolic compounds, are present and active in extant charophytes and bryophytes as model organisms for the conquest of land.