The Art and Science of Molecular Docking.

Molecular docking has become an essential part of a structural biologist's and medicinal chemist's toolkits. Given a chemical compound and the three-dimensional structure of a molecular target-for example, a protein-docking methods fit the compound into the target, predicting the compound's bound structure and binding energy. Docking can be used to discover novel ligands for a target by screening large virtual compound libraries. Docking can also provide a useful starting point for structure-based ligand optimization or for investigating a ligand's mechanism of action. Advances in computational methods, including both physics-based and machine learning approaches, as well as in complementary experimental techniques, are making docking an even more powerful tool. We review how docking works and how it can drive drug discovery and biological research. We also describe its current limitations and ongoing efforts to overcome them.

Multiscale drug screening for cardiac fibrosis identifies MD2 as a therapeutic target.

Cardiac fibrosis impairs cardiac function, but no effective clinical therapies exist. To address this unmet need, we employed a high-throughput screening for antifibrotic compounds using human induced pluripotent stem cell (iPSC)-derived cardiac fibroblasts (CFs). Counter-screening of the initial candidates using iPSC-derived cardiomyocytes and iPSC-derived endothelial cells excluded hits with cardiotoxicity. This screening process identified artesunate as the lead compound. Following profibrotic stimuli, artesunate inhibited proliferation, migration, and contraction in human primary CFs, reduced collagen deposition, and improved contractile function in 3D-engineered heart tissues. Artesunate also attenuated cardiac fibrosis and improved cardiac function in heart failure mouse models. Mechanistically, artesunate targeted myeloid differentiation factor 2 (MD2) and inhibited MD2/Toll-like receptor 4 (TLR4) signaling pathway, alleviating fibrotic gene expression in CFs. Our study leverages multiscale drug screening that integrates a human iPSC platform, tissue engineering, animal models, in silico simulations, and multiomics to identify MD2 as a therapeutic target for cardiac fibrosis.

Positive selection CRISPR screens reveal a druggable pocket in an oligosaccharyltransferase required for inflammatory signaling to NF-κB.

Nuclear factor κB (NF-κB) plays roles in various diseases. Many inflammatory signals, such as circulating lipopolysaccharides (LPSs), activate NF-κB via specific receptors. Using whole-genome CRISPR-Cas9 screens of LPS-treated cells that express an NF-κB-driven suicide gene, we discovered that the LPS receptor Toll-like receptor 4 (TLR4) is specifically dependent on the oligosaccharyltransferase complex OST-A for N-glycosylation and cell-surface localization. The tool compound NGI-1 inhibits OST complexes in vivo, but the underlying molecular mechanism remained unknown. We did a CRISPR base-editor screen for NGI-1-resistant variants of STT3A, the catalytic subunit of OST-A. These variants, in conjunction with cryoelectron microscopy studies, revealed that NGI-1 binds the catalytic site of STT3A, where it traps a molecule of the donor substrate dolichyl-PP-GlcNAc2-Man9-Glc3, suggesting an uncompetitive inhibition mechanism. Our results provide a rationale for and an initial step toward the development of STT3A-specific inhibitors and illustrate the power of contemporaneous base-editor and structural studies to define drug mechanism of action.

Proteome-scale movements and compartment connectivity during the eukaryotic cell cycle.

Cell cycle progression relies on coordinated changes in the composition and subcellular localization of the proteome. By applying two distinct convolutional neural networks on images of millions of live yeast cells, we resolved proteome-level dynamics in both concentration and localization during the cell cycle, with resolution of ∼20 subcellular localization classes. We show that a quarter of the proteome displays cell cycle periodicity, with proteins tending to be controlled either at the level of localization or concentration, but not both. Distinct levels of protein regulation are preferentially utilized for different aspects of the cell cycle, with changes in protein concentration being mostly involved in cell cycle control and changes in protein localization in the biophysical implementation of the cell cycle program. We present a resource for exploring global proteome dynamics during the cell cycle, which will aid in understanding a fundamental biological process at a systems level.

Highly multiplexed bioactivity screening reveals human and microbiota metabolome-GPCRome interactions.

The human body contains thousands of metabolites derived from mammalian cells, the microbiota, food, and medical drugs. Many bioactive metabolites act through the engagement of G-protein-coupled receptors (GPCRs); however, technological limitations constrain current explorations of metabolite-GPCR interactions. Here, we developed a highly multiplexed screening technology called PRESTO-Salsa that enables simultaneous assessment of nearly all conventional GPCRs (>300 receptors) in a single well of a 96-well plate. Using PRESTO-Salsa, we screened 1,041 human-associated metabolites against the GPCRome and uncovered previously unreported endogenous, exogenous, and microbial GPCR agonists. Next, we leveraged PRESTO-Salsa to generate an atlas of microbiome-GPCR interactions across 435 human microbiome strains from multiple body sites, revealing conserved patterns of cross-tissue GPCR engagement and activation of CD97/ADGRE5 by the Porphyromonas gingivalis protease gingipain K. These studies thus establish a highly multiplexed bioactivity screening technology and expose a diverse landscape of human, diet, drug, and microbiota metabolome-GPCRome interactions.

Massively parallel functional dissection of schizophrenia-associated noncoding genetic variants.

Schizophrenia (SCZ) is a highly heritable mental disorder with thousands of associated genetic variants located mostly in the noncoding space of the genome. Translating these associations into insights regarding the underlying pathomechanisms has been challenging because the causal variants, their mechanisms of action, and their target genes remain largely unknown. We implemented a massively parallel variant annotation pipeline (MVAP) to perform SCZ variant-to-function mapping at scale in disease-relevant neural cell types. This approach identified 620 functional variants (1.7%) that operate in a highly developmental context and neuronal-activity-dependent manner. Multimodal integration of epigenomic and CRISPRi screening data enabled us to link these functional variants to target genes, biological processes, and ultimately alterations of neuronal physiology. These results provide a multistage prioritization strategy to map functional single-nucleotide polymorphism (SNP)-to-gene-to-endophenotype relations and offer biological insights into the context-dependent molecular processes modulated by SCZ-associated genetic variation.

Positive selection of somatically mutated clones identifies adaptive pathways in metabolic liver disease.

Somatic mutations in nonmalignant tissues accumulate with age and injury, but whether these mutations are adaptive on the cellular or organismal levels is unclear. To interrogate genes in human metabolic disease, we performed lineage tracing in mice harboring somatic mosaicism subjected to nonalcoholic steatohepatitis (NASH). Proof-of-concept studies with mosaic loss of Mboat7, a membrane lipid acyltransferase, showed that increased steatosis accelerated clonal disappearance. Next, we induced pooled mosaicism in 63 known NASH genes, allowing us to trace mutant clones side by side. This in vivo tracing platform, which we coined MOSAICS, selected for mutations that ameliorate lipotoxicity, including mutant genes identified in human NASH. To prioritize new genes, additional screening of 472 candidates identified 23 somatic perturbations that promoted clonal expansion. In validation studies, liver-wide deletion of Tbx3, Bcl6, or Smyd2 resulted in protection against hepatic steatosis. Selection for clonal fitness in mouse and human livers identifies pathways that regulate metabolic disease.

Multi-chamber cardioids unravel human heart development and cardiac defects.

The number one cause of human fetal death are defects in heart development. Because the human embryonic heart is inaccessible and the impacts of mutations, drugs, and environmental factors on the specialized functions of different heart compartments are not captured by in vitro models, determining the underlying causes is difficult. Here, we established a human cardioid platform that recapitulates the development of all major embryonic heart compartments, including right and left ventricles, atria, outflow tract, and atrioventricular canal. By leveraging 2D and 3D differentiation, we efficiently generated progenitor subsets with distinct first, anterior, and posterior second heart field identities. This advance enabled the reproducible generation of cardioids with compartment-specific in vivo-like gene expression profiles, morphologies, and functions. We used this platform to unravel the ontogeny of signal and contraction propagation between interacting heart chambers and dissect how mutations, teratogens, and drugs cause compartment-specific defects in the developing human heart.

Trellis tree-based analysis reveals stromal regulation of patient-derived organoid drug responses.

Patient-derived organoids (PDOs) can model personalized therapy responses; however, current screening technologies cannot reveal drug response mechanisms or how tumor microenvironment cells alter therapeutic performance. To address this, we developed a highly multiplexed mass cytometry platform to measure post-translational modification (PTM) signaling, DNA damage, cell-cycle activity, and apoptosis in >2,500 colorectal cancer (CRC) PDOs and cancer-associated fibroblasts (CAFs) in response to clinical therapies at single-cell resolution. To compare patient- and microenvironment-specific drug responses in thousands of single-cell datasets, we developed "Trellis"-a highly scalable, tree-based treatment effect analysis method. Trellis single-cell screening revealed that on-target cell-cycle blockage and DNA-damage drug effects are common, even in chemorefractory PDOs. However, drug-induced apoptosis is rarer, patient-specific, and aligns with cancer cell PTM signaling. We find that CAFs can regulate PDO plasticity-shifting proliferative colonic stem cells (proCSCs) to slow-cycling revival colonic stem cells (revCSCs) to protect cancer cells from chemotherapy.

Comparative landscape of genetic dependencies in human and chimpanzee stem cells.

Comparative studies of great apes provide a window into our evolutionary past, but the extent and identity of cellular differences that emerged during hominin evolution remain largely unexplored. We established a comparative loss-of-function approach to evaluate whether human cells exhibit distinct genetic dependencies. By performing genome-wide CRISPR interference screens in human and chimpanzee pluripotent stem cells, we identified 75 genes with species-specific effects on cellular proliferation. These genes comprised coherent processes, including cell-cycle progression and lysosomal signaling, which we determined to be human-derived by comparison with orangutan cells. Human-specific robustness to CDK2 and CCNE1 depletion persisted in neural progenitor cells and cerebral organoids, supporting the G1-phase length hypothesis as a potential evolutionary mechanism in human brain expansion. Our findings demonstrate that evolutionary changes in human cells reshaped the landscape of essential genes and establish a platform for systematically uncovering latent cellular and molecular differences between species.

A transcription factor atlas of directed differentiation.

Transcription factors (TFs) regulate gene programs, thereby controlling diverse cellular processes and cell states. To comprehensively understand TFs and the programs they control, we created a barcoded library of all annotated human TF splice isoforms (>3,500) and applied it to build a TF Atlas charting expression profiles of human embryonic stem cells (hESCs) overexpressing each TF at single-cell resolution. We mapped TF-induced expression profiles to reference cell types and validated candidate TFs for generation of diverse cell types, spanning all three germ layers and trophoblasts. Targeted screens with subsets of the library allowed us to create a tailored cellular disease model and integrate mRNA expression and chromatin accessibility data to identify downstream regulators. Finally, we characterized the effects of combinatorial TF overexpression by developing and validating a strategy for predicting combinations of TFs that produce target expression profiles matching reference cell types to accelerate cellular engineering efforts.

TScan-II: A genome-scale platform for the de novo identification of CD4+ T cell epitopes.

CD4+ T cells play fundamental roles in orchestrating immune responses and tissue homeostasis. However, our inability to associate peptide human leukocyte antigen class-II (HLA-II) complexes with their cognate T cell receptors (TCRs) in an unbiased manner has hampered our understanding of CD4+ T cell function and role in pathologies. Here, we introduce TScan-II, a highly sensitive genome-scale CD4+ antigen discovery platform. This platform seamlessly integrates the endogenous HLA-II antigen-processing machinery in synthetic antigen-presenting cells and TCR signaling in T cells, enabling the simultaneous screening of multiple HLAs and TCRs. Leveraging genome-scale human, virome, and epitope mutagenesis libraries, TScan-II facilitates de novo antigen discovery and deep exploration of TCR specificity. We demonstrate TScan-II's potential for basic and translational research by identifying a non-canonical antigen for a cancer-reactive CD4+ T cell clone. Additionally, we identified two antigens for clonally expanded CD4+ T cells in Sjögren's disease, which bind distinct HLAs and are expressed in HLA-II-positive ductal cells within affected salivary glands.

Validating Small Molecule Chemical Probes for Biological Discovery.

Small molecule chemical probes are valuable tools for interrogating protein biological functions and relevance as a therapeutic target. Rigorous validation of chemical probe parameters such as cellular potency and selectivity is critical to unequivocally linking biological and phenotypic data resulting from treatment with a chemical probe to the function of a specific target protein. A variety of modern technologies are available to evaluate cellular potency and selectivity, target engagement, and functional response biomarkers of chemical probe compounds. Here, we review these technologies and the rationales behind using them for the characterization and validation of chemical probes. In addition, large-scale phenotypic characterization of chemical probes through chemical genetic screening is increasingly leading to a wealth of information on the cellular pharmacology and disease involvement of potential therapeutic targets. Extensive compound validation approaches and integration of phenotypic information will lay foundations for further use of chemical probes in biological discovery.

The phenotypic landscape of essential human genes.

Understanding the basis for cellular growth, proliferation, and function requires determining the roles of essential genes in diverse cellular processes, including visualizing their contributions to cellular organization and morphology. Here, we combined pooled CRISPR-Cas9-based functional screening of 5,072 fitness-conferring genes in human HeLa cells with microscopy-based imaging of DNA, the DNA damage response, actin, and microtubules. Analysis of >31 million individual cells identified measurable phenotypes for >90% of gene knockouts, implicating gene targets in specific cellular processes. Clustering of phenotypic similarities based on hundreds of quantitative parameters further revealed co-functional genes across diverse cellular activities, providing predictions for gene functions and associations. By conducting pooled live-cell screening of ∼450,000 cell division events for 239 genes, we additionally identified diverse genes with functional contributions to chromosome segregation. Our work establishes a resource detailing the consequences of disrupting core cellular processes that represents the functional landscape of essential human genes.

Cannabinoid receptor 1 antagonist genistein attenuates marijuana-induced vascular inflammation.

Epidemiological studies reveal that marijuana increases the risk of cardiovascular disease (CVD); however, little is known about the mechanism. Δ9-tetrahydrocannabinol (Δ9-THC), the psychoactive component of marijuana, binds to cannabinoid receptor 1 (CB1/CNR1) in the vasculature and is implicated in CVD. A UK Biobank analysis found that cannabis was an risk factor for CVD. We found that marijuana smoking activated inflammatory cytokines implicated in CVD. In silico virtual screening identified genistein, a soybean isoflavone, as a putative CB1 antagonist. Human-induced pluripotent stem cell-derived endothelial cells were used to model Δ9-THC-induced inflammation and oxidative stress via NF-κB signaling. Knockdown of the CB1 receptor with siRNA, CRISPR interference, and genistein attenuated the effects of Δ9-THC. In mice, genistein blocked Δ9-THC-induced endothelial dysfunction in wire myograph, reduced atherosclerotic plaque, and had minimal penetration of the central nervous system. Genistein is a CB1 antagonist that attenuates Δ9-THC-induced atherosclerosis.

Sustained deep-tissue voltage recording using a fast indicator evolved for two-photon microscopy.

Genetically encoded voltage indicators are emerging tools for monitoring voltage dynamics with cell-type specificity. However, current indicators enable a narrow range of applications due to poor performance under two-photon microscopy, a method of choice for deep-tissue recording. To improve indicators, we developed a multiparameter high-throughput platform to optimize voltage indicators for two-photon microscopy. Using this system, we identified JEDI-2P, an indicator that is faster, brighter, and more sensitive and photostable than its predecessors. We demonstrate that JEDI-2P can report light-evoked responses in axonal termini of Drosophila interneurons and the dendrites and somata of amacrine cells of isolated mouse retina. JEDI-2P can also optically record the voltage dynamics of individual cortical neurons in awake behaving mice for more than 30 min using both resonant-scanning and ULoVE random-access microscopy. Finally, ULoVE recording of JEDI-2P can robustly detect spikes at depths exceeding 400 µm and report voltage correlations in pairs of neurons.

In vivo CRISPR screens identify the E3 ligase Cop1 as a modulator of macrophage infiltration and cancer immunotherapy target.

Despite remarkable clinical efficacy of immune checkpoint blockade (ICB) in cancer treatment, ICB benefits for triple-negative breast cancer (TNBC) remain limited. Through pooled in vivo CRISPR knockout (KO) screens in syngeneic TNBC mouse models, we found that deletion of the E3 ubiquitin ligase Cop1 in cancer cells decreases secretion of macrophage-associated chemokines, reduces tumor macrophage infiltration, enhances anti-tumor immunity, and strengthens ICB response. Transcriptomics, epigenomics, and proteomics analyses revealed that Cop1 functions through proteasomal degradation of the C/ebpδ protein. The Cop1 substrate Trib2 functions as a scaffold linking Cop1 and C/ebpδ, which leads to polyubiquitination of C/ebpδ. In addition, deletion of the E3 ubiquitin ligase Cop1 in cancer cells stabilizes C/ebpδ to suppress expression of macrophage chemoattractant genes. Our integrated approach implicates Cop1 as a target for improving cancer immunotherapy efficacy in TNBC by regulating chemokine secretion and macrophage infiltration in the tumor microenvironment.

In vivo CRISPR screening reveals nutrient signaling processes underpinning CD8+ T cell fate decisions.

How early events in effector T cell (TEFF) subsets tune memory T cell (TMEM) responses remains incompletely understood. Here, we systematically investigated metabolic factors in fate determination of TEFF and TMEM cells using in vivo pooled CRISPR screening, focusing on negative regulators of TMEM responses. We found that amino acid transporters Slc7a1 and Slc38a2 dampened the magnitude of TMEM differentiation, in part through modulating mTORC1 signaling. By integrating genetic and systems approaches, we identified cellular and metabolic heterogeneity among TEFF cells, with terminal effector differentiation associated with establishment of metabolic quiescence and exit from the cell cycle. Importantly, Pofut1 (protein-O-fucosyltransferase-1) linked GDP-fucose availability to downstream Notch-Rbpj signaling, and perturbation of this nutrient signaling axis blocked terminal effector differentiation but drove context-dependent TEFF proliferation and TMEM development. Our study establishes that nutrient uptake and signaling are key determinants of T cell fate and shape the quantity and quality of TMEM responses.

Visualizing and Modulating Mitophagy for Therapeutic Studies of Neurodegeneration.

Dysfunctional mitochondria accumulate in many human diseases. Accordingly, mitophagy, which removes these mitochondria through lysosomal degradation, is attracting broad attention. Due to uncertainties in the operational principles of conventional mitophagy probes, however, the specificity and quantitativeness of their readouts are disputable. Thorough investigation of the behaviors and fates of fluorescent proteins inside and outside lysosomes enabled us to develop an indicator for mitophagy, mito-SRAI. Through strict control of its mitochondrial targeting, we were able to monitor mitophagy in fixed biological samples more reproducibly than before. Large-scale image-based high-throughput screening led to the discovery of a hit compound that induces selective mitophagy of damaged mitochondria. In a mouse model of Parkinsons disease, we found that dopaminergic neurons selectively failed to execute mitophagy that promoted their survival within lesions. These results show that mito-SRAI is an essential tool for quantitative studies of mitochondrial quality control.

High-Throughput Discovery and Characterization of Human Transcriptional Effectors.

Thousands of proteins localize to the nucleus; however, it remains unclear which contain transcriptional effectors. Here, we develop HT-recruit, a pooled assay where protein libraries are recruited to a reporter, and their transcriptional effects are measured by sequencing. Using this approach, we measure gene silencing and activation for thousands of domains. We find a relationship between repressor function and evolutionary age for the KRAB domains, discover that Homeodomain repressor strength is collinear with Hox genetic organization, and identify activities for several domains of unknown function. Deep mutational scanning of the CRISPRi KRAB maps the co-repressor binding surface and identifies substitutions that improve stability/silencing. By tiling 238 proteins, we find repressors as short as ten amino acids. Finally, we report new activator domains, including a divergent KRAB. These results provide a resource of 600 human proteins containing effectors and demonstrate a scalable strategy for assigning functions to protein domains.