Flexible scaffold-based cheminformatics approach for polypharmacological drug design.

Effective treatments for complex central nervous system (CNS) disorders require drugs with polypharmacology and multifunctionality, yet designing such drugs remains a challenge. Here, we present a flexible scaffold-based cheminformatics approach (FSCA) for the rational design of polypharmacological drugs. FSCA involves fitting a flexible scaffold to different receptors using different binding poses, as exemplified by IHCH-7179, which adopted a "bending-down" binding pose at 5-HT2AR to act as an antagonist and a "stretching-up" binding pose at 5-HT1AR to function as an agonist. IHCH-7179 demonstrated promising results in alleviating cognitive deficits and psychoactive symptoms in mice by blocking 5-HT2AR for psychoactive symptoms and activating 5-HT1AR to alleviate cognitive deficits. By analyzing aminergic receptor structures, we identified two featured motifs, the "agonist filter" and "conformation shaper," which determine ligand binding pose and predict activity at aminergic receptors. With these motifs, FSCA can be applied to the design of polypharmacological ligands at other receptors.

Structure of the Nanobody-Stabilized Active State of the Kappa Opioid Receptor.

The κ-opioid receptor (KOP) mediates the actions of opioids with hallucinogenic, dysphoric, and analgesic activities. The design of KOP analgesics devoid of hallucinatory and dysphoric effects has been hindered by an incomplete structural and mechanistic understanding of KOP agonist actions. Here, we provide a crystal structure of human KOP in complex with the potent epoxymorphinan opioid agonist MP1104 and an active-state-stabilizing nanobody. Comparisons between inactive- and active-state opioid receptor structures reveal substantial conformational changes in the binding pocket and intracellular and extracellular regions. Extensive structural analysis and experimental validation illuminate key residues that propagate larger-scale structural rearrangements and transducer binding that, collectively, elucidate the structural determinants of KOP pharmacology, function, and biased signaling. These molecular insights promise to accelerate the structure-guided design of safer and more effective κ-opioid receptor therapeutics.

The Psychiatric Cell Map Initiative: A Convergent Systems Biological Approach to Illuminating Key Molecular Pathways in Neuropsychiatric Disorders.

Although gene discovery in neuropsychiatric disorders, including autism spectrum disorder, intellectual disability, epilepsy, schizophrenia, and Tourette disorder, has accelerated, resulting in a large number of molecular clues, it has proven difficult to generate specific hypotheses without the corresponding datasets at the protein complex and functional pathway level. Here, we describe one path forward-an initiative aimed at mapping the physical and genetic interaction networks of these conditions and then using these maps to connect the genomic data to neurobiology and, ultimately, the clinic. These efforts will include a team of geneticists, structural biologists, neurobiologists, systems biologists, and clinicians, leveraging a wide array of experimental approaches and creating a collaborative infrastructure necessary for long-term investigation. This initiative will ultimately intersect with parallel studies that focus on other diseases, as there is a significant overlap with genes implicated in cancer, infectious disease, and congenital heart defects.

Crystal Structure of an LSD-Bound Human Serotonin Receptor.

The prototypical hallucinogen LSD acts via serotonin receptors, and here we describe the crystal structure of LSD in complex with the human serotonin receptor 5-HT2B. The complex reveals conformational rearrangements to accommodate LSD, providing a structural explanation for the conformational selectivity of LSD's key diethylamide moiety. LSD dissociates exceptionally slow from both 5-HT2BR and 5-HT2AR-a major target for its psychoactivity. Molecular dynamics (MD) simulations suggest that LSD's slow binding kinetics may be due to a "lid" formed by extracellular loop 2 (EL2) at the entrance to the binding pocket. A mutation predicted to increase the mobility of this lid greatly accelerates LSD's binding kinetics and selectively dampens LSD-mediated ß-arrestin2 recruitment. This study thus reveals an unexpected binding mode of LSD; illuminates key features of its kinetics, stereochemistry, and signaling; and provides a molecular explanation for LSD's actions at human serotonin receptors. PAPERCLIP.

A whole-slide foundation model for digital pathology from real-world data.

Digital pathology poses unique computational challenges, as a standard gigapixel slide may comprise tens of thousands of image tiles1-3. Prior models have often resorted to subsampling a small portion of tiles for each slide, thus missing the important slide-level context4. Here we present Prov-GigaPath, a whole-slide pathology foundation model pretrained on 1.3 billion 256 × 256 pathology image tiles in 171,189 whole slides from Providence, a large US health network comprising 28 cancer centres. The slides originated from more than 30,000 patients covering 31 major tissue types. To pretrain Prov-GigaPath, we propose GigaPath, a novel vision transformer architecture for pretraining gigapixel pathology slides. To scale GigaPath for slide-level learning with tens of thousands of image tiles, GigaPath adapts the newly developed LongNet5 method to digital pathology. To evaluate Prov-GigaPath, we construct a digital pathology benchmark comprising 9 cancer subtyping tasks and 17 pathomics tasks, using both Providence and TCGA data6. With large-scale pretraining and ultra-large-context modelling, Prov-GigaPath attains state-of-the-art performance on 25 out of 26 tasks, with significant improvement over the second-best method on 18 tasks. We further demonstrate the potential of Prov-GigaPath on vision-language pretraining for pathology7,8 by incorporating the pathology reports. In sum, Prov-GigaPath is an open-weight foundation model that achieves state-of-the-art performance on various digital pathology tasks, demonstrating the importance of real-world data and whole-slide modelling.

The rise of baobab trees in Madagascar.

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

Recognition of methamphetamine and other amines by trace amine receptor TAAR1.

Trace amine-associated receptor 1 (TAAR1), the founding member of a nine-member family of trace amine receptors, is responsible for recognizing a range of biogenic amines in the brain, including the endogenous ß-phenylethylamine (ß-PEA)1 as well as methamphetamine2, an abused substance that has posed a severe threat to human health and society3. Given its unique physiological role in the brain, TAAR1 is also an emerging target for a range of neurological disorders including schizophrenia, depression and drug addiction2,4,5. Here we report structures of human TAAR1-G-protein complexes bound to methamphetamine and ß-PEA as well as complexes bound to RO5256390, a TAAR1-selective agonist, and SEP-363856, a clinical-stage dual agonist for TAAR1 and serotonin receptor 5-HT1AR (refs. 6,7). Together with systematic mutagenesis and functional studies, the structures reveal the molecular basis of methamphetamine recognition and underlying mechanisms of ligand selectivity and polypharmacology between TAAR1 and other monoamine receptors. We identify a lid-like extracellular loop 2 helix/loop structure and a hydrogen-bonding network in the ligand-binding pockets, which may contribute to the ligand recognition in TAAR1. These findings shed light on the ligand recognition mode and activation mechanism for TAAR1 and should guide the development of next-generation therapeutics for drug addiction and various neurological disorders.

De novo protein design with a denoising diffusion network independent of pretrained structure prediction models.

The recent success of RFdiffusion, a method for protein structure design with a denoising diffusion probabilistic model, has relied on fine-tuning the RoseTTAFold structure prediction network for protein backbone denoising. Here, we introduce SCUBA-diffusion (SCUBA-D), a protein backbone denoising diffusion probabilistic model freshly trained by considering co-diffusion of sequence representation to enhance model regularization and adversarial losses to minimize data-out-of-distribution errors. While matching the performance of the pretrained RoseTTAFold-based RFdiffusion in generating experimentally realizable protein structures, SCUBA-D readily generates protein structures with not-yet-observed overall folds that are different from those predictable with RoseTTAFold. The accuracy of SCUBA-D was confirmed by the X-ray structures of 16 designed proteins and a protein complex, and by experiments validating designed heme-binding proteins and Ras-binding proteins. Our work shows that deep generative models of images or texts can be fruitfully extended to complex physical objects like protein structures by addressing outstanding issues such as the data-out-of-distribution errors.

Efficient Fizeau drag from Dirac electrons in monolayer graphene.

Fizeau demonstrated in 1850 that the speed of light can be modified when it is propagating in moving media1. However, such control of the light speed has not been achieved efficiently with a fast-moving electron media by passing an electrical current. Because the strong electromagnetic coupling between the electron and light leads to the collective excitation of plasmon polaritons, it is hypothesized that Fizeau drag in electron flow systems manifests as a plasmonic Doppler effect. Experimental observation of the plasmonic Doppler effect in electronic systems has been challenge because the plasmon propagation speed is much faster than the electron drift velocity in conventional noble metals. Here we report direct observation of Fizeau drag of plasmon polaritons in strongly biased monolayer graphene by exploiting the high electron mobility and the slow plasmon propagation of massless Dirac electrons. The large bias current in graphene creates a fast-drifting Dirac electron medium hosting the plasmon polariton. This results in non-reciprocal plasmon propagation, where plasmons moving with the drifting electron media propagate at an enhanced speed. We measure the Doppler-shifted plasmon wavelength using cryogenic near-field infrared nanoscopy, which directly images the plasmon polariton mode in the biased graphene at low temperature. We observe a plasmon wavelength difference of up to 3.6 per cent between a plasmon moving with and a plasmon moving against the drifting electron media. Our findings on the plasmonic Doppler effect provide opportunities for electrical control of non-reciprocal surface plasmon polaritons in non-equilibrium systems.

An additional proofreader contributes to DNA replication fidelity in mycobacteria.

The removal of mis-incorporated nucleotides by proofreading activity ensures DNA replication fidelity. Whereas the ε-exonuclease DnaQ is a well-established proofreader in the model organism Escherichia coli, it has been shown that proofreading in a majority of bacteria relies on the polymerase and histidinol phosphatase (PHP) domain of replicative polymerase, despite the presence of a DnaQ homolog that is structurally and functionally distinct from E. coli DnaQ. However, the biological functions of this type of noncanonical DnaQ remain unclear. Here, we provide independent evidence that noncanonical DnaQ functions as an additional proofreader for mycobacteria. Using the mutation accumulation assay in combination with whole-genome sequencing, we showed that depletion of DnaQ in Mycolicibacterium smegmatis leads to an increased mutation rate, resulting in AT-biased mutagenesis and increased insertions/deletions in the homopolymer tract. Our results showed that mycobacterial DnaQ binds to the ß clamp and functions synergistically with the PHP domain proofreader to correct replication errors. Furthermore, the loss of dnaQ results in replication fork dysfunction, leading to attenuated growth and increased mutagenesis on subinhibitory fluoroquinolones potentially due to increased vulnerability to fork collapse. By analyzing the sequence polymorphism of dnaQ in clinical isolates of Mycobacterium tuberculosis (Mtb), we demonstrated that a naturally evolved DnaQ variant prevalent in Mtb lineage 4.3 may enable hypermutability and is associated with drug resistance. These results establish a coproofreading model and suggest a division of labor between DnaQ and PHP domain proofreader. This study also provides real-world evidence that a mutator-driven evolutionary pathway may exist during the adaptation of Mtb.

Aging and comprehensive molecular profiling in acute myeloid leukemia.

Acute myeloid leukemia (AML) is an aging-related and heterogeneous hematopoietic malignancy. In this study, a total of 1,474 newly diagnosed AML patients with RNA sequencing data were enrolled, and targeted or whole exome sequencing data were obtained in 94% cases. The correlation of aging-related factors including age and clonal hematopoiesis (CH), gender, and genomic/transcriptomic profiles (gene fusions, genetic mutations, and gene expression networks or pathways) was systematically analyzed. Overall, AML patients aged 60 y and older showed an apparently dismal prognosis. Alongside age, the frequency of gene fusions defined in the World Health Organization classification decreased, while the positive rate of gene mutations, especially CH-related ones, increased. Additionally, the number of genetic mutations was higher in gene fusion-negative (GF-) patients than those with GF. Based on the status of CH- and myelodysplastic syndromes (MDS)-related mutations, three mutant subgroups were identified among the GF- AML cohort, namely, CH-AML, CH-MDS-AML, and other GF- AML. Notably, CH-MDS-AML demonstrated a predominance of elderly and male cases, cytopenia, and significantly adverse clinical outcomes. Besides, gene expression networks including HOXA/B, platelet factors, and inflammatory responses were most striking features associated with aging and poor prognosis in AML. Our work has thus unraveled the intricate regulatory circuitry of interactions among different age, gender, and molecular groups of AML.

Pleiotropy of autism-associated chromatin regulators.

Gene ontology analyses of high-confidence autism spectrum disorder (ASD) risk genes highlight chromatin regulation and synaptic function as major contributors to pathobiology. Our recent functional work in vivo has additionally implicated tubulin biology and cellular proliferation. As many chromatin regulators, including the ASD risk genes ADNP and CHD3, are known to directly regulate both tubulins and histones, we studied the five chromatin regulators most strongly associated with ASD (ADNP, CHD8, CHD2, POGZ and KMT5B) specifically with respect to tubulin biology. We observe that all five localize to microtubules of the mitotic spindle in vitro in human cells and in vivo in Xenopus. Investigation of CHD2 provides evidence that mutations present in individuals with ASD cause a range of microtubule-related phenotypes, including disrupted localization of the protein at mitotic spindles, cell cycle stalling, DNA damage and cell death. Lastly, we observe that ASD genetic risk is significantly enriched among tubulin-associated proteins, suggesting broader relevance. Together, these results provide additional evidence that the role of tubulin biology and cellular proliferation in ASD warrants further investigation and highlight the pitfalls of relying solely on annotated gene functions in the search for pathological mechanisms.

Comparative analyses of gene networks mediating cancer metastatic potentials across lineage types.

Studies have identified genes and molecular pathways regulating cancer metastasis. However, it remains largely unknown whether metastatic potentials of cancer cells from different lineage types are driven by the same or different gene networks. Here, we aim to address this question through integrative analyses of 493 human cancer cells' transcriptomic profiles and their metastatic potentials in vivo. Using an unsupervised approach and considering both gene coexpression and protein-protein interaction networks, we identify different gene networks associated with various biological pathways (i.e. inflammation, cell cycle, and RNA translation), the expression of which are correlated with metastatic potentials across subsets of lineage types. By developing a regularized random forest regression model, we show that the combination of the gene module features expressed in the native cancer cells can predict their metastatic potentials with an overall Pearson correlation coefficient of 0.90. By analyzing transcriptomic profile data from cancer patients, we show that these networks are conserved in vivo and contribute to cancer aggressiveness. The intrinsic expression levels of these networks are correlated with drug sensitivity. Altogether, our study provides novel comparative insights into cancer cells' intrinsic gene networks mediating metastatic potentials across different lineage types, and our results can potentially be useful for designing personalized treatments for metastatic cancers.

A new paradigm for applying deep learning to protein-ligand interaction prediction.

Protein-ligand interaction prediction presents a significant challenge in drug design. Numerous machine learning and deep learning (DL) models have been developed to accurately identify docking poses of ligands and active compounds against specific targets. However, current models often suffer from inadequate accuracy or lack practical physical significance in their scoring systems. In this research paper, we introduce IGModel, a novel approach that utilizes the geometric information of protein-ligand complexes as input for predicting the root mean square deviation of docking poses and the binding strength (pKd, the negative value of the logarithm of binding affinity) within the same prediction framework. This ensures that the output scores carry intuitive meaning. We extensively evaluate the performance of IGModel on various docking power test sets, including the CASF-2016 benchmark, PDBbind-CrossDocked-Core and DISCO set, consistently achieving state-of-the-art accuracies. Furthermore, we assess IGModel's generalizability and robustness by evaluating it on unbiased test sets and sets containing target structures generated by AlphaFold2. The exceptional performance of IGModel on these sets demonstrates its efficacy. Additionally, we visualize the latent space of protein-ligand interactions encoded by IGModel and conduct interpretability analysis, providing valuable insights. This study presents a novel framework for DL-based prediction of protein-ligand interactions, contributing to the advancement of this field. The IGModel is available at GitHub repository https://github.com/zchwang/IGModel.

Progress and opportunities of foundation models in bioinformatics.

Bioinformatics has undergone a paradigm shift in artificial intelligence (AI), particularly through foundation models (FMs), which address longstanding challenges in bioinformatics such as limited annotated data and data noise. These AI techniques have demonstrated remarkable efficacy across various downstream validation tasks, effectively representing diverse biological entities and heralding a new era in computational biology. The primary goal of this survey is to conduct a general investigation and summary of FMs in bioinformatics, tracing their evolutionary trajectory, current research landscape, and methodological frameworks. Our primary focus is on elucidating the application of FMs to specific biological problems, offering insights to guide the research community in choosing appropriate FMs for tasks like sequence analysis, structure prediction, and function annotation. Each section delves into the intricacies of the targeted challenges, contrasting the architectures and advancements of FMs with conventional methods and showcasing their utility across different biological domains. Further, this review scrutinizes the hurdles and constraints encountered by FMs in biology, including issues of data noise, model interpretability, and potential biases. This analysis provides a theoretical groundwork for understanding the circumstances under which certain FMs may exhibit suboptimal performance. Lastly, we outline prospective pathways and methodologies for the future development of FMs in biological research, facilitating ongoing innovation in the field. This comprehensive examination not only serves as an academic reference but also as a roadmap for forthcoming explorations and applications of FMs in biology.

Serotonin 2A Receptor (5-HT2AR) Agonists: Psychedelics and Non-Hallucinogenic Analogues as Emerging Antidepressants.

Psychedelics make up a group of psychoactive compounds that induce hallucinogenic effects by activating the serotonin 2A receptor (5-HT2AR). Clinical trials have demonstrated the traditional psychedelic substances like psilocybin as a class of rapid-acting and long-lasting antidepressants. However, there is a pressing need for rationally designed 5-HT2AR agonists that possess optimal pharmacological profiles in order to fully reveal the therapeutic potential of these agonists and identify safer drug candidates devoid of hallucinogenic effects. This Perspective provides an overview of the structure-activity relationships of existing 5-HT2AR agonists based on their chemical classifications and discusses recent advancements in understanding their molecular pharmacology at a structural level. The encouraging clinical outcomes of psychedelics in depression treatment have sparked drug discovery endeavors aimed at developing novel 5-HT2AR agonists with improved subtype selectivity and signaling bias properties, which could serve as safer and potentially nonhallucinogenic antidepressants. These efforts can be significantly expedited through the utilization of structure-based methods and functional selectivity-directed screening.

Nucleic Acid Sensing by STING Induces an IFN-like Antiviral Response in a Marine Invertebrate.

The cytosolic detection of pathogen-derived nucleic acids has evolved as an essential strategy for host innate immune defense in mammals. One crucial component in this process is the stimulator of IFN genes (STING), which acts as a vital signaling adaptor, connecting the cytosolic detection of DNA by cyclic GMP-AMP (cGAMP) synthase (cGAS) to the downstream type I IFN signaling pathway. However, this process remains elusive in invertebrates. In this study, we present evidence demonstrating that STING, an ortholog found in a marine invertebrate (shrimp) called Litopenaeus vannamei, can directly detect DNA and initiate an IFN-like antiviral response. Unlike its homologs in other eukaryotic organisms, which exclusively function as sensors for cyclic dinucleotides, shrimp STING has the ability to bind to both double-stranded DNA and cyclic dinucleotides, including 2'3'-cGAMP. In vivo, shrimp STING can directly sense DNA nucleic acids from an infected virus, accelerate IFN regulatory factor dimerization and nuclear translocation, induce the expression of an IFN functional analog protein (Vago4), and finally establish an antiviral state. Taken together, our findings unveil a novel double-stranded DNA-STING-IKKε-IRF-Vago antiviral axis in an arthropod, providing valuable insights into the functional origins of DNA-sensing pathways in evolution.

RDBSB: a database for catalytic bioparts with experimental evidence.

Catalytic bioparts are fundamental to the design, construction and optimization of biological systems for specific metabolic pathways. However, the functional characterization information of these bioparts is frequently dispersed across multiple databases and literature sources, posing significant challenges to the effective design and optimization of specific chassis or cell factories. We developed the Registry and Database of Bioparts for Synthetic Biology (RDBSB), a comprehensive resource encompassing 83 193 curated catalytic bioparts with experimental evidences. RDBSB offers their detailed qualitative and quantitative catalytic information, including critical parameters such as activities, substrates, optimal pH and temperature, and chassis specificity. The platform features an interactive search engine, visualization tools and analysis utilities such as biopart finder, structure prediction and pathway design tools. Additionally, RDBSB promotes community engagement through a catalytic bioparts submission system to facilitate rapid data sharing and utilization. To date, RDBSB has supported the contribution of >1000 catalytic bioparts. We anticipate that the database will significantly enhance the resources available for pathway design in synthetic biology and serve essential tools for researchers. RDBSB is freely available at https://www.biosino.org/rdbsb/.

Bioorthogonal labeling and profiling of N6-isopentenyladenosine (i6A) modified RNA.

Chemical modifications in RNAs play crucial roles in diversifying their structures and regulating numerous biochemical processes. Since the 1990s, several hydrophobic prenyl-modifications have been discovered in various RNAs. Prenyl groups serve as precursors for terpenes and many other biological molecules. The processes of prenylation in different macromolecules have been extensively studied. We introduce here a novel chemical biology toolkit that not only labels i6A, a prenyl-modified RNA residue, by leveraging the unique reactivity of the prenyl group, but also provides a general strategy to incorporate fluorescence functionalities into RNAs for molecular tracking purposes. Our findings revealed that iodine-mediated cyclization reactions of the prenyl group occur rapidly, transforming i6A from a hydrogen-bond acceptor to a donor. Based on this reactivity, we developed an Iodine-Mediated Cyclization and Reverse Transcription (IMCRT) tRNA-seq method, which can profile all nine endogenous tRNAs containing i6A residues in Saccharomyces cerevisiae with single-base resolution. Furthermore, under stress conditions, we observed a decline in i6A levels in budding yeast, accompanied by significant decrease of mutation rate at A37 position. Thus, the IMCRT tRNA-seq method not only permits semi-quantification of i6A levels in tRNAs but also holds potential for transcriptome-wide detection and analysis of various RNA species containing i6A modifications.

Closed-loop recyclable membranes enabled by covalent adaptable networks for water purification.

In the state-of-the-art membrane industry, membranes have linear life cycles and are commonly disposed of by landfill or incineration, sacrificing their sustainability. To date, little or no thought is given in the design phase to the end-of-life management of membranes. For the first time, we have innovated high-performance sustainable membranes, which can be closed-loop recycled after long-term usage for water purification. By synergizing membrane technology and dynamic covalent chemistry, covalent adaptable networks (CANs) with thermally reversible Diels-Alder (DA) adducts were synthesized and employed to fabricate integrally skinned asymmetric membranes via the nonsolvent-induced phase separation technique. Due to the stable and reversible features of CAN, the closed-loop recyclable membranes exhibit excellent mechanical properties and thermal and chemical stabilities as well as separation performance, which are comparable to or even higher than the state-of-the-art nonrecyclable membranes. Moreover, the used membranes can be closed-loop recycled with consistent properties and separation performance by depolymerization to remove contaminants, followed by refabrication into new membranes through the dissociation and reformation of DA adducts. This study may fill in the gaps in closed-loop recycling of membranes and inspire the advancement of sustainable membranes for a green membrane industry.