CryptoCEN: A Co-Expression Network for Cryptococcus neoformans reveals novel proteins involved in DNA damage repair

Elucidating gene function is a major goal in biology, especially among non-model organisms. However, doing so is complicated by the fact that molecular conservation does not always mirror functional conservation, and that complex relationships among genes are responsible for encoding pathways and higher-order biological processes. Co-expression, a promising approach for predicting gene function, relies on the general principal that genes with similar expression patterns across multiple conditions will likely be involved in the same biological process. For Cryptococcus neoformans, a prevalent human fungal pathogen greatly diverged from model yeasts, approximately 60% of the predicted genes in the genome lack functional annotations. Here, we leveraged a large amount of publicly available transcriptomic data to generate a C. neoformans Co-Expression Network (CryptoCEN), successfully recapitulating known protein networks, predicting gene function, and enabling insights into the principles influencing co-expression. With 100% predictive accuracy, we used CryptoCEN to identify 13 new DNA damage response genes, underscoring the utility of guilt-by-association for determining gene function. Overall, co-expression is a powerful tool for uncovering gene function, and decreases the experimental tests needed to identify functions for currently under-annotated genes.

Multiple ParA/MinD ATPases coordinate the positioning of disparate cargos in a bacterial cell

In eukaryotes, linear motor proteins govern intracellular transport and organization. In bacteria, where linear motors involved in spatial regulation are absent, the ParA/MinD family of ATPases organize an array of genetic- and protein-based cellular cargos. The positioning of these cargos has been independently investigated to varying degrees in several bacterial species. However, it remains unclear how multiple ParA/MinD ATPases can coordinate the positioning of diverse cargos in the same cell. Here, we find that over a third of sequenced bacterial genomes encode multiple ParA/MinD ATPases. We identify an organism (Halothiobacillus neapolitanus) with seven ParA/MinD ATPases, demonstrate that five of these are each dedicated to the spatial regulation of a single cellular cargo, and define potential specificity determinants for each system. Furthermore, we show how these positioning reactions can influence each other, stressing the importance of understanding how organelle trafficking, chromosome segregation, and cell division are coordinated in bacterial cells. Together, our data show how multiple ParA/MinD ATPases coexist and function to position a diverse set of fundamental cargos in the same bacterial cell.

Candida albicans selection for human commensalism results in substantial within-host diversity without decreasing fitness for invasive disease

Candida albicans is a frequent colonizer of human mucosal surfaces as well as an opportunistic pathogen. C. albicans is remarkably versatile in its ability to colonize diverse host sites with differences in oxygen and nutrient availability, pH, immune responses, and resident microbes, among other cues. It is unclear how the genetic background of a commensal colonizing population can influence the shift to pathogenicity. Therefore, we examined 910 commensal isolates from 35 healthy donors to identify host niche-specific adaptations. We demonstrate that healthy people are reservoirs for genotypically and phenotypically diverse C. albicans strains. Using limited diversity exploitation, we identified a single nucleotide change in the uncharacterized ZMS1 transcription factor that was sufficient to drive hyper invasion into agar. We found that SC5314 was significantly different from the majority of both commensal and bloodstream isolates in its ability to induce host cell death. However, our commensal strains retained the capacity to cause disease in the Galleria model of systemic infection, including outcompeting the SC5314 reference strain during systemic competition assays. This study provides a global view of commensal strain variation and within-host strain diversity of C. albicans and suggests that selection for commensalism in humans does not result in a fitness cost for invasive disease.

The life and times of endogenous opioid peptides: Updated understanding of synthesis, spatiotemporal dynamics, and the clinical impact in alcohol use disorder

The opioid G-protein coupled receptors (GPCRs) strongly modulate many of the central nervous system structures that contribute to neurological and psychiatric disorders including pain, major depressive disorder, and substance use disorders. To better treat these and related diseases, it is essential to understand the signaling of their endogenous ligands. In this review, we focus on what is known and unknown about the regulation of the over two dozen endogenous peptides with high affinity for one or more of the opioid receptors. We briefly describe which peptides are produced, with a particular focus on the recently proposed possible synthesis pathways for the endomorphins. Next, we describe examples of endogenous opioid peptide expression organization in several neural circuits and how they appear to be released from specific neural compartments that vary across brain regions. We discuss current knowledge regarding the strength of neural activity required to drive endogenous opioid peptide release, clues about how far peptides diffuse from release sites, and their extracellular lifetime after release. Finally, as a translational example, we discuss the mechanisms of action of naltrexone (NTX), which is used clinically to treat alcohol use disorder. NTX is a synthetic morphine analog that non-specifically antagonizes the action of most endogenous opioid peptides developed in the 1960s and FDA approved in the 1980s. We review recent studies clarifying the precise endogenous activity that NTX prevents. Together, the works described here highlight the challenges and opportunities the complex opioid system presents as a therapeutic target.

A human liver organoid screening platform for DILI risk prediction

BACKGROUND & AIMS

Drug-induced liver injury (DILI), both intrinsic and idiosyncratic, causes frequent morbidity, mortality, clinical trial failures and post-approval withdrawal. This suggests an unmet need for improved in vitro models for DILI risk prediction that can account for diverse host genetics and other clinical factors. In this study, we evaluated the utility of human liver organoids (HLOs) for high-throughput DILI risk prediction and in an organ-on-chip system.

METHODS

HLOs were derived from three separate iPSC lines and benchmarked on two platforms for their ability to model in vitro liver function and identify hepatotoxic compounds using biochemical assays for albumin, ALT, AST, microscopy-based morphological profiling, and single-cell transcriptomics: i) HLOs dispersed in 384-well-formatted plates and exposed to a library of compounds; ii) HLOs adapted to a liver-on-chip system.

RESULTS

Dispersed HLOs derived from the three iPSC lines had similar DILI predictive capacity as intact HLOs in a high-throughput screening format, allowing for measurable IC₅₀ values of compound cytotoxicity. Distinct morphological differences were observed in cells treated with drugs exerting differing mechanisms of toxicity. On-chip HLOs significantly increased albumin production, CYP450 expression, and ALT/AST release when treated with known hepatoxic drugs compared to dispersed HLOs and primary human hepatocytes. On-chip HLOs were able to predict the synergistic hepatotoxicity of tenofovir-inarigivir and displayed steatosis and mitochondrial perturbation, via phenotypic and transcriptomic analysis, on exposure to fialuridine and acetaminophen, respectively.

CONCLUSIONS

The high-throughput and liver-on-chip systems exhibit enhanced in vivo-like functions and demonstrate the potential utility of these platforms for DILI risk assessment. Tenofovir-inarigivr-associated hepatotoxicity was observed and correlates with the clinical manifestation of DILI observed in patients.

IMPACT AND IMPLICATIONS

Idiosyncratic (spontaneous, patient-specific) drug-induced liver injury (DILI) is difficult to study due to the lack of liver models that function as human liver tissue and are adaptable for large-scale drug screening. Human liver organoids grown from patient stem cells respond to known DILI-causing drugs in both a high-throughput and on a physiological "chip" culture system. These platforms show promise for researchers in their use as predictive models for novel drugs before entering clinical trials and as a potential in vitro diagnostic tool. Our findings support further development of patient-derived liver organoid lines and their use in the context of DILI research.

Development of an automated screen for Kv7.2 potassium channels and discovery of a new agonist chemotype

To identify pore domain ligands on Kv7.2 potassium ion channels, we compared wild-type (WT) and W236L mutant Kv7.2 channels in a series of assays with previously validated and novel agonist chemotypes. Positive controls were retigabine, flupirtine, and RL-81; i.e. Kv7.2 channel activators that significantly shift voltage-dependent activation to more negative potentials (ΔV50) at 5 µM. We identified 6 new compounds that exhibited differential enhancing activity between WT and W236L mutant channels. Whole cell patch-clamp electrophysiology studies were conducted to identify Kv7.2. Kv7.2/3, Kv7.4, and Kv7.5 selectivity. Our results validate the SyncroPatch platform and establish new structure activity relationships (SAR). Specifically, in addition to selective Kv7.2, Kv7.2/3, Kv7.4. and Kv7.5 agonists, we identified a novel chemotype, ZK-21, a 4-aminotetrahydroquinoline that is distinct from any of the previously described Kv7 channel modifiers. Using flexible receptor docking, ZK-21 was predicted to be stabilized by W236 and bind perpendicular to retigabine, burying the benzyl carbamate group into a tunnel reaching the core of the pore domain.

Drug-induced phospholipidosis confounds drug repurposing for SARS-CoV-2

Repurposing drugs as treatments for COVID-19, the disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has drawn much attention. Beginning with sigma receptor ligands and expanding to other drugs from screening in the field, we became concerned that phospholipidosis was a shared mechanism underlying the antiviral activity of many repurposed drugs. For all of the 23 cationic amphiphilic drugs we tested, including hydroxychloroquine, azithromycin, amiodarone, and four others already in clinical trials, phospholipidosis was monotonically correlated with antiviral efficacy. Conversely, drugs active against the same targets that did not induce phospholipidosis were not antiviral. Phospholipidosis depends on the physicochemical properties of drugs and does not reflect specific target-based activities-rather, it may be considered a toxic confound in early drug discovery. Early detection of phospholipidosis could eliminate these artifacts, enabling a focus on molecules with therapeutic potential.

Valproic acid-induced changes of 4D nuclear morphology in astrocyte cells

Histone deacetylase inhibitors, such as valproic acid (VPA), have important clinical therapeutic and cellular reprogramming applications. They induce chromatin reorganization that is associated with altered cellular morphology. However, there is a lack of comprehensive characterization of VPA-induced changes of nuclear size and shape. Here, we quantify 3D nuclear morphology of primary human astrocyte cells treated with VPA over time (hence, 4D). We compared volumetric and surface-based representations and identified seven features that jointly discriminate between normal and treated cells with 85% accuracy on day 7. From day 3, treated nuclei were more elongated and flattened and then continued to morphologically diverge from controls over time, becoming larger and more irregular. On day 7, most of the size and shape descriptors demonstrated significant differences between treated and untreated cells, including a 24% increase in volume and 6% reduction in extent (shape regularity) for treated nuclei. Overall, we show that 4D morphometry can capture how chromatin reorganization modulates the size and shape of the nucleus over time. These nuclear structural alterations may serve as a biomarker for histone (de-)acetylation events and provide insights into mechanisms of astrocytes-to-neurons reprogramming.

Morphological Cell Profiling of SARS-CoV-2 Infection Identifies Drug Repurposing Candidates for COVID-19

The global spread of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and the associated disease COVID-19, requires therapeutic interventions that can be rapidly identified and translated to clinical care. Traditional drug discovery methods have a >90% failure rate and can take 10–15 years from target identification to clinical use. In contrast, drug repurposing can significantly accelerate translation. We developed a quantitative high-throughput screen to identify efficacious agents against SARS-CoV-2. From a library of 1,425 FDA-approved compounds and clinical candidates, we identified 17 dose-responsive compounds with in vitro antiviral efficacy in human liver Huh7 cells and confirmed antiviral efficacy in human colon carcinoma Caco-2, human prostate adenocarcinoma LNCaP, and in a physiologic relevant model of alveolar epithelial type 2 cells (iAEC2s). Additionally, we found that inhibitors of the Ras/Raf/MEK/ERK signaling pathway exacerbate SARS-CoV-2 infection in vitro. Notably, we discovered that lactoferrin, a glycoprotein classically found in secretory fluids, including mammalian milk, inhibits SARS-CoV-2 infection in the nanomolar range in all cell models with multiple modes of action, including blockage of virus attachment to cellular heparan sulfate and enhancement of interferon responses. Given its safety profile, lactoferrin is a readily translatable therapeutic option for the management of COVID-19.

Comparative host-coronavirus protein interaction networks reveal pan-viral disease mechanisms

The COVID-19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a grave threat to public health and the global economy. SARS-CoV-2 is closely related to the more lethal but less transmissible coronaviruses SARS-CoV-1 and Middle East respiratory syndrome coronavirus (MERS-CoV). Here, we have carried out comparative viral-human protein-protein interaction and viral protein localization analyses for all three viruses. Subsequent functional genetic screening identified host factors that functionally impinge on coronavirus proliferation, including Tom70, a mitochondrial chaperone protein that interacts with both SARS-CoV-1 and SARS-CoV-2 ORF9b, an interaction we structurally characterized using cryo-electron microscopy. Combining genetically validated host factors with both COVID-19 patient genetic data and medical billing records identified molecular mechanisms and potential drug treatments that merit further molecular and clinical study.

A SARS-CoV-2 protein interaction map reveals targets for drug repurposing

A newly described coronavirus named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is the causative agent of coronavirus disease 2019 (COVID-19), has infected over 2.3 million people, led to the death of more than 160,000 individuals and caused worldwide social and economic disruption1,2. There are no antiviral drugs with proven clinical efficacy for the treatment of COVID-19, nor are there any vaccines that prevent infection with SARS-CoV-2, and efforts to develop drugs and vaccines are hampered by the limited knowledge of the molecular details of how SARS-CoV-2 infects cells. Here we cloned, tagged and expressed 26 of the 29 SARS-CoV-2 proteins in human cells and identified the human proteins that physically associated with each of the SARS-CoV-2 proteins using affinity-purification mass spectrometry, identifying 332 high-confidence protein-protein interactions between SARS-CoV-2 and human proteins. Among these, we identify 66 druggable human proteins or host factors targeted by 69 compounds (of which, 29 drugs are approved by the US Food and Drug Administration, 12 are in clinical trials and 28 are preclinical compounds). We screened a subset of these in multiple viral assays and found two sets of pharmacological agents that displayed antiviral activity: inhibitors of mRNA translation and predicted regulators of the sigma-1 and sigma-2 receptors. Further studies of these host-factor-targeting agents, including their combination with drugs that directly target viral enzymes, could lead to a therapeutic regimen to treat COVID-19.

A SARS-CoV-2 protein interaction map reveals targets for drug repurposing

A newly described coronavirus named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is the causative agent of coronavirus disease 2019 (COVID-19), has infected over 2.3 million people, led to the death of more than 160,000 individuals and caused worldwide social and economic disruption^1,2. There are no antiviral drugs with proven clinical efficacy for the treatment of COVID-19, nor are there any vaccines that prevent infection with SARS-CoV-2, and efforts to develop drugs and vaccines are hampered by the limited knowledge of the molecular details of how SARS-CoV-2 infects cells. Here we cloned, tagged and expressed 26 of the 29 SARS-CoV-2 proteins in human cells and identified the human proteins that physically associated with each of the SARS-CoV-2 proteins using affinity-purification mass spectrometry, identifying 332 high-confidence protein-protein interactions between SARS-CoV-2 and human proteins. Among these, we identify 66 druggable human proteins or host factors targeted by 69 compounds (of which, 29 drugs are approved by the US Food and Drug Administration, 12 are in clinical trials and 28 are preclinical compounds). We screened a subset of these in multiple viral assays and found two sets of pharmacological agents that displayed antiviral activity: inhibitors of mRNA translation and predicted regulators of the sigma-1 and sigma-2 receptors. Further studies of these host-factor-targeting agents, including their combination with drugs that directly target viral enzymes, could lead to a therapeutic regimen to treat COVID-19.

A SARS-CoV-2-Human Protein-Protein Interaction Map Reveals Drug Targets and Potential Drug-Repurposing

An outbreak of the novel coronavirus SARS-CoV-2, the causative agent of COVID-19 respiratory disease, has infected over 290,000 people since the end of 2019, killed over 12,000, and caused worldwide social and economic disruption 1,2 . There are currently no antiviral drugs with proven efficacy nor are there vaccines for its prevention. Unfortunately, the scientific community has little knowledge of the molecular details of SARS-CoV-2 infection. To illuminate this, we cloned, tagged and expressed 26 of the 29 viral proteins in human cells and identified the human proteins physically associated with each using affinity-purification mass spectrometry (AP-MS), which identified 332 high confidence SARS-CoV-2-human protein-protein interactions (PPIs). Among these, we identify 66 druggable human proteins or host factors targeted by 69 existing FDA-approved drugs, drugs in clinical trials and/or preclinical compounds, that we are currently evaluating for efficacy in live SARS-CoV-2 infection assays. The identification of host dependency factors mediating virus infection may provide key insights into effective molecular targets for developing broadly acting antiviral therapeutics against SARS-CoV-2 and other deadly coronavirus strains.

Ultra-large library docking for discovering new chemotypes

Despite intense interest in expanding chemical space, libraries containing hundreds-of-millions to billions of diverse molecules have remained inaccessible. Here we investigate structure-based docking of 170 million make-on-demand compounds from 130 well-characterized reactions. The resulting library is diverse, representing over 10.7 million scaffolds that are otherwise unavailable. For each compound in the library, docking against AmpC β-lactamase (AmpC) and the D4 dopamine receptor were simulated. From the top-ranking molecules, 44 and 549 compounds were synthesized and tested for interactions with AmpC and the D4 dopamine receptor, respectively. We found a phenolate inhibitor of AmpC, which revealed a group of inhibitors without known precedent. This molecule was optimized to 77 nM, which places it among the most potent non-covalent AmpC inhibitors known. Crystal structures of this and other AmpC inhibitors confirmed the docking predictions. Against the D4 dopamine receptor, hit rates fell almost monotonically with docking score, and a hit-rate versus score curve predicted that the library contained 453,000 ligands for the D4 dopamine receptor. Of 81 new chemotypes discovered, 30 showed submicromolar activity, including a 180-pM subtype-selective agonist of the D4 dopamine receptor.

Prediction of enzymatic pathways by integrative pathway mapping

The functions of most proteins are yet to be determined. The function of an enzyme is often defined by its interacting partners, including its substrate and product, and its role in larger metabolic networks. Here, we describe a computational method that predicts the functions of orphan enzymes by organizing them into a linear metabolic pathway. Given candidate enzyme and metabolite pathway members, this aim is achieved by finding those pathways that satisfy structural and network restraints implied by varied input information, including that from virtual screening, chemoinformatics, genomic context analysis, and ligand -binding experiments. We demonstrate this integrative pathway mapping method by predicting the L-gulonate catabolic pathway in Haemophilus influenzae Rd KW20. The prediction was subsequently validated experimentally by enzymology, crystallography, and metabolomics. Integrative pathway mapping by satisfaction of structural and network restraints is extensible to molecular networks in general and thus formally bridges the gap between structural biology and systems biology.

The Rosetta All-Atom Energy Function for Macromolecular Modeling and Design

Over the past decade, the Rosetta biomolecular modeling suite has informed diverse biological questions and engineering challenges ranging from interpretation of low-resolution structural data to design of nanomaterials, protein therapeutics, and vaccines. Central to Rosetta's success is the energy function: a model parametrized from small-molecule and X-ray crystal structure data used to approximate the energy associated with each biomolecule conformation. This paper describes the mathematical models and physical concepts that underlie the latest Rosetta energy function, called the Rosetta Energy Function 2015 (REF15). Applying these concepts, we explain how to use Rosetta energies to identify and analyze the features of biomolecular models. Finally, we discuss the latest advances in the energy function that extend its capabilities from soluble proteins to also include membrane proteins, peptides containing noncanonical amino acids, small molecules, carbohydrates, nucleic acids, and other macromolecules.

A Web Resource for Standardized Benchmark Datasets, Metrics, and Rosetta Protocols for Macromolecular Modeling and Design

The development and validation of computational macromolecular modeling and design methods depend on suitable benchmark datasets and informative metrics for comparing protocols. In addition, if a method is intended to be adopted broadly in diverse biological applications, there needs to be information on appropriate parameters for each protocol, as well as metrics describing the expected accuracy compared to experimental data. In certain disciplines, there exist established benchmarks and public resources where experts in a particular methodology are encouraged to supply their most efficient implementation of each particular benchmark. We aim to provide such a resource for protocols in macromolecular modeling and design. We present a freely accessible web resource (https://kortemmelab.ucsf.edu/benchmarks) to guide the development of protocols for protein modeling and design. The site provides benchmark datasets and metrics to compare the performance of a variety of modeling protocols using different computational sampling methods and energy functions, providing a “best practice” set of parameters for each method. Each benchmark has an associated downloadable benchmark capture archive containing the input files, analysis scripts, and tutorials for running the benchmark. The captures may be run with any suitable modeling method; we supply command lines for running the benchmarks using the Rosetta software suite. We have compiled initial benchmarks for the resource spanning three key areas: prediction of energetic effects of mutations, protein design, and protein structure prediction, each with associated state-of-the-art modeling protocols. With the help of the wider macromolecular modeling community, we hope to expand the variety of benchmarks included on the website and continue to evaluate new iterations of current methods as they become available.

Role of electrostatic repulsion in controlling pH-dependent conformational changes of viral fusion proteins

Viral fusion proteins undergo dramatic conformational transitions during membrane fusion. For viruses that enter through the endosome, these conformational rearrangements are typically pH sensitive. Here, we provide a comprehensive review of the molecular interactions that govern pH-dependent rearrangements and introduce a paradigm for electrostatic residue pairings that regulate progress through the viral fusion coordinate. Analysis of structural data demonstrates a significant role for side-chain protonation in triggering conformational change. To characterize this behavior, we identify two distinct residue pairings, which we define as Histidine-Cation (HisCat) and Anion-Anion (AniAni) interactions. These side-chain pairings destabilize a particular conformation via electrostatic repulsion through side-chain protonation. Furthermore, two energetic control mechanisms, thermodynamic and kinetic, regulate these structural transitions. This review expands on the current literature by identification of these residue clusters, discussion of data demonstrating their function, and speculation of how these residue pairings contribute to the energetic controls.

Scientific benchmarks for guiding macromolecular energy function improvement

Accurate energy functions are critical to macromolecular modeling and design. We describe new tools for identifying inaccuracies in energy functions and guiding their improvement, and illustrate the application of these tools to the improvement of the Rosetta energy function. The feature analysis tool identifies discrepancies between structures deposited in the PDB and low-energy structures generated by Rosetta; these likely arise from inaccuracies in the energy function. The optE tool optimizes the weights on the different components of the energy function by maximizing the recapitulation of a wide range of experimental observations. We use the tools to examine three proposed modifications to the Rosetta energy function: improving the unfolded state energy model (reference energies), using bicubic spline interpolation to generate knowledge-based torisonal potentials, and incorporating the recently developed Dunbrack 2010 rotamer library (Shapovalov & Dunbrack, 2011).