PROMALS3D: multiple protein sequence alignment enhanced with evolutionary and three-dimensional structural information

Mutation of the Plasmodium falciparum Flavokinase Confers Resistance to Roseoflavin and 8-Aminoriboflavin

The riboflavin analogues, roseoflavin and 8-aminoriboflavin, inhibit malaria parasite proliferation by targeting riboflavin utilization. To determine their mechanism of action, we generated roseoflavin-resistant parasites by in vitro evolution. Relative to wild-type, these parasites were 4-fold resistant to roseoflavin and cross-resistant to 8-aminoriboflavin. Whole genome sequencing of the resistant parasites revealed a missense mutation leading to an amino acid change (L672H) in the gene coding for a putative flavokinase (PfFK), the enzyme responsible for converting riboflavin into the cofactor flavin mononucleotide (FMN). To confirm that the L672H mutation is responsible for the phenotype, we generated parasites with the missense mutation incorporated into the PfFK gene. The IC50 values for roseoflavin and 8-aminoriboflavin against the roseoflavin-resistant parasites created through in vitro evolution were indistinguishable from those against parasites in which the missense mutation was introduced into the native PfFK. We also generated two parasite lines episomally expressing GFP-tagged versions of either the wild-type or mutant forms of PfFK. We found that PfFK-GFP localizes to the parasite cytosol and that immunopurified PfFK-GFP phosphorylated riboflavin, roseoflavin, and 8-aminoriboflavin. The L672H mutation increased the KM for roseoflavin, explaining the resistance phenotype. Mutant PfFK is no longer capable of phosphorylating 8-aminoriboflavin, but its antiplasmodial activity against resistant parasites can still be antagonized by increasing the extracellular concentration of riboflavin, consistent with it also inhibiting parasite growth through competitive inhibition of PfFK. Our findings, therefore, are consistent with roseoflavin and 8-aminoriboflavin inhibiting parasite proliferation by inhibiting riboflavin phosphorylation and via the generation of toxic flavin cofactor analogues.

Reversible male contraception by targeted inhibition of serine/threonine kinase 33

Men or mice with homozygous serine/threonine kinase 33 (STK33) mutations are sterile owing to defective sperm morphology and motility. To chemically evaluate STK33 for male contraception with STK33-specific inhibitors, we screened our multibillion-compound collection of DNA-encoded chemical libraries, uncovered potent STK33-specific inhibitors, determined the STK33 kinase domain structure bound with a truncated hit CDD-2211, and generated an optimized hit CDD-2807 that demonstrates nanomolar cellular potency (half-maximal inhibitory concentration = 9.2 nanomolar) and favorable metabolic stability. In mice, CDD-2807 exhibited no toxicity, efficiently crossed the blood-testis barrier, did not accumulate in brain, and induced a reversible contraceptive effect that phenocopied genetic STK33 perturbations without altering testis size. Thus, STK33 is a chemically validated, nonhormonal contraceptive target, and CDD-2807 is an effective tool compound.

Impact of Combinatorial Histone Modifications on Acetyllysine Recognition by the ATAD2 and ATAD2B Bromodomains

The ATPase family AAA+ domain containing 2 (ATAD2) protein and its paralog ATAD2B have a C-terminal bromodomain (BRD) that functions as a reader of acetylated lysine residues on histone proteins. Using a structure-function approach, we investigated the ability of the ATAD2/B BRDs to select acetylated lysine among multiple histone post-translational modifications. The ATAD2B BRD can bind acetylated histone ligands that also contain adjacent methylation or phosphorylation marks, while the presence of these modifications significantly weakened the acetyllysine binding activity of the ATAD2 BRD. Our structural studies provide mechanistic insights into how ATAD2/B BRD-binding pocket residues coordinate the acetyllysine group in the context of adjacent post-translational modifications. Furthermore, we investigated how sequence changes in amino acids of the histone ligands impact the recognition of an adjacent acetyllysine residue. Our study highlights how the interplay between multiple combinations of histone modifications influences the reader activity of the ATAD2/B BRDs, resulting in distinct binding modes.

A Practical Guide for the Quality Evaluation of Fluobodies/Chromobodies

BACKGROUND

Fluorescent proteins (FPs) are pivotal reagents for flow cytometry analysis or fluorescent microscopy. A new generation of immunoreagents (fluobodies/chromobodies) has been developed by fusing recombinant nanobodies to FPs.

METHODS

We analyzed the quality of such biomolecules by a combination of gel filtration and SDS-PAGE to identify artefacts due to aggregation or material degradation.

RESULTS

In the SDS-PAGE run, unexpected bands corresponding to separate fluobodies were evidenced and characterized as either degradation products or artefacts that systematically resulted in the presence of specific FPs and some experimental conditions. The elimination of N-terminal methionine from FPs did not impair the appearance of FP fragments, whereas the stability and migration characteristics of some FP constructs were strongly affected by heating in loading buffer, which is a step samples undergo before electrophoretic separation.

CONCLUSIONS

In this work, we provide explanations for some odd results observed during the quality control of fluobodies and summarize practical suggestions for the choice of the most convenient FPs to fuse to antibody fragments.

Streptomyces umbrella toxin particles block hyphal growth of competing species

Streptomyces are a genus of ubiquitous soil bacteria from which the majority of clinically utilized antibiotics derive1. The production of these antibacterial molecules reflects the relentless competition Streptomyces engage in with other bacteria, including other Streptomyces species1,2. Here we show that in addition to small-molecule antibiotics, Streptomyces produce and secrete antibacterial protein complexes that feature a large, degenerate repeat-containing polymorphic toxin protein. A cryo-electron microscopy structure of these particles reveals an extended stalk topped by a ringed crown comprising the toxin repeats scaffolding five lectin-tipped spokes, which led us to name them umbrella particles. Streptomyces coelicolor encodes three umbrella particles with distinct toxin and lectin composition. Notably, supernatant containing these toxins specifically and potently inhibits the growth of select Streptomyces species from among a diverse collection of bacteria screened. For one target, Streptomyces griseus, inhibition relies on a single toxin and that intoxication manifests as rapid cessation of vegetative hyphal growth. Our data show that Streptomyces umbrella particles mediate competition among vegetative mycelia of related species, a function distinct from small-molecule antibiotics, which are produced at the onset of reproductive growth and act broadly3,4. Sequence analyses suggest that this role of umbrella particles extends beyond Streptomyces, as we identified umbrella loci in nearly 1,000 species across Actinobacteria.

Double-stranded RNA sequencing reveals distinct riboviruses associated with thermoacidophilic bacteria from hot springs in Japan

Metatranscriptome sequencing expanded the known diversity of the bacterial RNA virome, suggesting that additional riboviruses infecting bacterial hosts remain to be discovered. Here we employed double-stranded RNA sequencing to recover complete genome sequences of two ribovirus groups from acidic hot springs in Japan. One group, denoted hot spring riboviruses (HsRV), consists of viruses with distinct RNA-directed RNA polymerases (RdRPs) that seem to be intermediates between typical ribovirus RdRPs and viral reverse transcriptases. This group forms a distinct phylum, Artimaviricota, or even kingdom within the realm Riboviria. We identified viruses encoding HsRV-like RdRPs in marine water, river sediments and salt marshes, indicating that this group is widespread beyond extreme ecosystems. The second group, denoted hot spring partiti-like viruses (HsPV), forms a distinct branch within the family Partitiviridae. The genome architectures of HsRV and HsPV and their identification in bacteria-dominated habitats suggest that these viruses infect thermoacidophilic bacteria.

Hiding in plain sight: The discovery of complete genomes of 11 hypothetical spindle-shaped viruses that putatively infect mesophilic ammonia-oxidizing archaea

The genome of a putative Nitrosopumilaceae virus with a hypothetical spindle-shaped particle morphology was identified in the Yangshan Harbour metavirome from the East China Sea through protein similarity comparison and structure analysis. This discovery was accompanied by a set of 10 geographically dispersed close relatives found in the environmental virus datasets from typical locations of ammonia-oxidizing archaeon distribution. Its host prediction was supported by iPHoP prediction and protein sequence similarity. The structure of the predicted major capsid protein, together with the overall N-glycosylation site, the transmembrane helices prediction, the hydrophilicity profile, and the docking simulation of the major capsid proteins, indicate that these viruses resemble spindle-shaped viruses. It suggests a similarly assembled structure and, consequently, a possibly spindle-shaped morphology of these newly discovered archaeal viruses.

Recursive genome engineering decodes the evolutionary origin of an essential thymidylate kinase activity in Pseudomonas putida KT2440

IMPORTANCE

Investigating fundamental aspects of metabolism is vital for advancing our understanding of the diverse biochemical capabilities and biotechnological applications of bacteria. The origin of the essential thymidylate kinase function in the model bacterium Pseudomonas putida KT2440, seemingly interrupted due to the presence of a large genomic island that disrupts the cognate gene, eluded a satisfactory explanation thus far. This is a first-case example of an essential metabolic function, likely acquired by horizontal gene transfer, which "landed" in a locus encoding the same activity. As such, foreign DNA encoding an essential dNMPK could immediately adjust to the recipient host-instead of long-term accommodation and adaptation. Understanding how these functions evolve is a major biological question, and the work presented here is a decisive step toward this direction. Furthermore, identifying essential and accessory genes facilitates removing those deemed irrelevant in industrial settings-yielding genome-reduced cell factories with enhanced properties and genetic stability.

Putative photosensitivity in internal light organs (organs of Pesta) of deep-sea sergestid shrimps

Many marine species can regulate the intensity of bioluminescence from their ventral photophores in order to counterilluminate, a camouflage technique whereby animals closely match the intensity of the downwelling illumination blocked by their bodies, thereby hiding their silhouettes. Recent studies on autogenic cuticular photophores in deep-sea shrimps indicate that the photophores themselves are light sensitive. Here, our results suggest photosensitivity in a second type of autogenic photophore, the internal organs of Pesta, found in deep-sea sergestid shrimps. Experiments were conducted onboard ship on live specimens, exposing the animals to bright light, which resulted in ultrastructural changes that matched those seen in crustacean eyes during the photoreceptor membrane turnover, a process that is crucial for the proper functioning of photosensitive components. In addition, RNA-seq studies demonstrated the expression of visual opsins and phototransduction genes in photophore tissue that are known to play a role in light detection, and electrophysiological measurements indicated that the light organs are responding to light received by the eyes. The long sought after mechanism of counterillumination remains unknown, but evidence of photosensitivity in photophores may indicate a dual functionality of light detection and emission.

A nuclear receptor HR96-related gene underlies large trans-driven differences in detoxification gene expression in a generalist herbivore

The role, magnitude, and molecular nature of trans-driven expression variation underlying the upregulation of detoxification genes in pesticide resistant arthropod populations has remained enigmatic. In this study, we performed expression quantitative trait locus (eQTL) mapping (n = 458) between a pesticide resistant and a susceptible strain of the generalist herbivore and crop pest Tetranychus urticae. We found that a single trans eQTL hotspot controlled large differences in the expression of a subset of genes in different detoxification gene families, as well as other genes associated with host plant use. As established by additional genetic approaches including RNAi gene knockdown, a duplicated gene with a nuclear hormone receptor HR96-related ligand-binding domain was identified as causal for the expression differences between strains. The presence of a large family of HR96-related genes in T. urticae may enable modular control of detoxification and host plant use genes, facilitating this species' known and rapid evolution to diverse pesticides and host plants.

cGLRs are a diverse family of pattern recognition receptors in innate immunity

Cyclic GMP-AMP synthase (cGAS) is an enzyme in human cells that controls an immune response to cytosolic DNA. Upon binding DNA, cGAS synthesizes a nucleotide signal 2'3'-cGAMP that activates STING-dependent downstream immunity. Here, we discover that cGAS-like receptors (cGLRs) constitute a major family of pattern recognition receptors in innate immunity. Building on recent analysis in Drosophila, we identify >3,000 cGLRs present in nearly all metazoan phyla. A forward biochemical screening of 150 animal cGLRs reveals a conserved mechanism of signaling including response to dsDNA and dsRNA ligands and synthesis of isomers of the nucleotide signals cGAMP, c-UMP-AMP, and c-di-AMP. Combining structural biology and in vivo analysis in coral and oyster animals, we explain how synthesis of distinct nucleotide signals enables cells to control discrete cGLR-STING signaling pathways. Our results reveal cGLRs as a widespread family of pattern recognition receptors and establish molecular rules that govern nucleotide signaling in animal immunity.

Distinct groups of RNA viruses associated with thermoacidophilic bacteria

Recent massive metatranscriptome mining substantially expanded the diversity of the bacterial RNA virome, suggesting that additional groups of riboviruses infecting bacterial hosts remain to be discovered. We employed full length double-stranded (ds) RNA sequencing for identification of riboviruses associated with microbial consortia dominated by bacteria and archaea in acidic hot springs in Japan. Whole sequences of two groups of multisegmented riboviruses genomes were obtained. One group, which we denoted hot spring riboviruses (HsRV), consists of unusual viruses with distinct RNA-dependent RNA polymerases (RdRPs) that seem to be intermediates between typical ribovirus RdRPs and viral reverse transcriptases. We also identified viruses encoding HsRV-like RdRPs in moderate aquatic environments, including marine water, river sediments and salt marsh, indicating that this previously overlooked ribovirus group is not restricted to the extreme ecosystem. The HsRV-like viruses are candidates for a distinct phylum or even kingdom within the viral realm Riboviria. The second group, denoted hot spring partiti-like viruses (HsPV), is a distinct branch within the family Partitiviridae. All genome segments in both these groups of viruses display the organization typical of bacterial riboviruses, where multiple open reading frames encoding individual proteins are preceded by ribosome-binding sites. Together with the identification in bacteria-dominated habitats, this genome architecture indicates that riboviruses of these distinct groups infect thermoacidophilic bacterial hosts.

CARD-like domains mediate anti-phage defense in bacterial gasdermin systems

Caspase recruitment domains (CARDs) and pyrin domains are important facilitators of inflammasome activity and pyroptosis. Upon pathogen recognition by NLR proteins, CARDs recruit and activate caspases, which, in turn, activate gasdermin pore forming proteins to and induce pyroptotic cell death. Here we show that CARD-like domains are present in defense systems that protect bacteria against phage. The bacterial CARD is essential for protease-mediated activation of certain bacterial gasdermins, which promote cell death once phage infection is recognized. We further show that multiple anti-phage defense systems utilize CARD-like domains to activate a variety of cell death effectors. We find that these systems are triggered by a conserved immune evasion protein that phages use to overcome the bacterial defense system RexAB, demonstrating that phage proteins inhibiting one defense system can activate another. We also detect a phage protein with a predicted CARD-like structure that can inhibit the CARD-containing bacterial gasdermin system. Our results suggest that CARD domains represent an ancient component of innate immune systems conserved from bacteria to humans, and that CARD-dependent activation of gasdermins is conserved in organisms across the tree of life.

Crystal structure of the 3-ketodihydrosphingosine reductase TSC10 from Cryptococcus neoformans

The second step in the de novo sphingolipid biosynthesis is the reduction of 3-ketodihydrosphingosine by 3-ketodihydrosphingosine reductase (KDSR) to produce dihydrosphingosine (sphinganine). Fungal TSC10 and mammalian KDSR (also named FVT-1) proteins are the enzymes responsible for this process and they belong to the short-chain dehydrogenase/reductase (SDR) superfamily. Albeit that both fungal and mammalian 3-ketodihydrosphingosine reductases were identified more than a decade ago, no structure of these enzymes from any species has been experimentally determined. Here we report the crystal structure of the catalytic domain of TSC10 from Cryptococcus neoformans in complex with NADPH. cnTSC10 adopts a Rossmann fold with a central seven-stranded β-sheet flanked by α-helices on both sides. Several regions are disordered that include the segment connecting the serine and tyrosine residues of the catalytic triad, the so-called 'substrate loop', and the C-terminal region that often participates in homo-tetramerization in other SDRs. In addition, the cofactor NADPH is not fully ordered. These structural features indicate that the catalytic site of cnTSC10 possesses significant flexibility. cnTSC10 is predominantly dimeric in solution while a minor portion of the protein forms homo-tetramer. The crystal structure reveals that the homo-dimer interface involves both hydrophobic and hydrophilic interactions mediated by helices α4 and α5, as well as the loop connecting strand β4 and helix α4. Because residues forming hydrogen bonds and salt bridges in the dimer interface are not conserved between fungal TSC10 and mammalian KDSR proteins, it might be possible to develop inhibitors that selectively target fungal TSC10 dimerization.

A tissue injury sensing and repair pathway distinct from host pathogen defense

Pathogen infection and tissue injury are universal insults that disrupt homeostasis. Innate immunity senses microbial infections and induces cytokines/chemokines to activate resistance mechanisms. Here, we show that, in contrast to most pathogen-induced cytokines, interleukin-24 (IL-24) is predominately induced by barrier epithelial progenitors after tissue injury and is independent of microbiome or adaptive immunity. Moreover, Il24 ablation in mice impedes not only epidermal proliferation and re-epithelialization but also capillary and fibroblast regeneration within the dermal wound bed. Conversely, ectopic IL-24 induction in the homeostatic epidermis triggers global epithelial-mesenchymal tissue repair responses. Mechanistically, Il24 expression depends upon both epithelial IL24-receptor/STAT3 signaling and hypoxia-stabilized HIF1α, which converge following injury to trigger autocrine and paracrine signaling involving IL-24-mediated receptor signaling and metabolic regulation. Thus, parallel to innate immune sensing of pathogens to resolve infections, epithelial stem cells sense injury signals to orchestrate IL-24-mediated tissue repair.

Crystal structure of an L-type lectin domain from archaea

The crystal structures of an L-type lectin domain from Methanocaldococcus jannaschii in apo and mannose-bound forms have been determined. A thorough investigation of L-type lectin domains from several organisms provides insight into the differences in these domains from different kingdoms of life. While the overall fold of the L-type lectin domain is conserved, differences in the lengths of the carbohydrate-binding loops and significant variations in the Mn²⁺ -binding site compared to the Ca²⁺ -binding site are observed. Furthermore, the sequence and phylogenetic analyses suggest that the archaeal L-type lectin domain is evolutionarily closer to the plant legume lectins than to its bacterial or animal counterparts. This is the first report of the biochemical, structural, sequence, and phylogenetic analyses of an L-type lectin domain from archaea and serves to enhance our understanding of the species-specific differences and evolution of L-type lectin domains.

cGLRs are a diverse family of pattern recognition receptors in animal innate immunity

cGAS (cyclic GMP-AMP synthase) is an enzyme in human cells that controls an immune response to cytosolic DNA. Upon binding DNA, cGAS synthesizes a nucleotide signal 2'3'-cGAMP that activates the protein STING and downstream immunity. Here we discover cGAS-like receptors (cGLRs) constitute a major family of pattern recognition receptors in animal innate immunity. Building on recent analysis in Drosophila , we use a bioinformatic approach to identify >3,000 cGLRs present in nearly all metazoan phyla. A forward biochemical screen of 140 animal cGLRs reveals a conserved mechanism of signaling including response to dsDNA and dsRNA ligands and synthesis of alternative nucleotide signals including isomers of cGAMP and cUMP-AMP. Using structural biology, we explain how synthesis of distinct nucleotide signals enables cells to control discrete cGLR-STING signaling pathways. Together our results reveal cGLRs as a widespread family of pattern recognition receptors and establish molecular rules that govern nucleotide signaling in animal immunity.

Archaeal DNA-import apparatus is homologous to bacterial conjugation machinery

Conjugation is a major mechanism of horizontal gene transfer promoting the spread of antibiotic resistance among human pathogens. It involves establishing a junction between a donor and a recipient cell via an extracellular appendage known as the mating pilus. In bacteria, the conjugation machinery is encoded by plasmids or transposons and typically mediates the transfer of cognate mobile genetic elements. Much less is known about conjugation in archaea. Here, we determine atomic structures by cryo-electron microscopy of three conjugative pili, two from hyperthermophilic archaea (Aeropyrum pernix and Pyrobaculum calidifontis) and one encoded by the Ti plasmid of the bacterium Agrobacterium tumefaciens, and show that the archaeal pili are homologous to bacterial mating pili. However, the archaeal conjugation machinery, known as Ced, has been 'domesticated', that is, the genes for the conjugation machinery are encoded on the chromosome rather than on mobile genetic elements, and mediates the transfer of cellular DNA.

Mining metatranscriptomes reveals a vast world of viroid-like circular RNAs

Viroids and viroid-like covalently closed circular (ccc) RNAs are minimal replicators that typically encode no proteins and hijack cellular enzymes for replication. The extent and diversity of viroid-like agents are poorly understood. We developed a computational pipeline to identify viroid-like cccRNAs and applied it to 5,131 metatranscriptomes and 1,344 plant transcriptomes. The search yielded 11,378 viroid-like cccRNAs spanning 4,409 species-level clusters, a 5-fold increase compared to the previously identified viroid-like elements. Within this diverse collection, we discovered numerous putative viroids, satellite RNAs, retrozymes, and ribozy-like viruses. Diverse ribozyme combinations and unusual ribozymes within the cccRNAs were identified. Self-cleaving ribozymes were identified in ambiviruses, some mito-like viruses and capsid-encoding satellite virus-like cccRNAs. The broad presence of viroid-like cccRNAs in diverse transcriptomes and ecosystems implies that their host range is far broader than currently known, and matches to CRISPR spacers suggest that some cccRNAs replicate in prokaryotes.

Two distinct gut microbial pathways contribute to meta-organismal production of phenylacetylglutamine with links to cardiovascular disease

Recent studies show gut microbiota-dependent metabolism of dietary phenylalanine into phenylacetic acid (PAA) is critical in phenylacetylglutamine (PAGln) production, a metabolite linked to atherosclerotic cardiovascular disease (ASCVD). Accordingly, microbial enzymes involved in this transformation are of interest. Using genetic manipulation in selected microbes and monocolonization experiments in gnotobiotic mice, we identify two distinct gut microbial pathways for PAA formation; one is catalyzed by phenylpyruvate:ferredoxin oxidoreductase (PPFOR) and the other by phenylpyruvate decarboxylase (PPDC). PPFOR and PPDC play key roles in gut bacterial PAA production via oxidative and non-oxidative phenylpyruvate decarboxylation, respectively. Metagenomic analyses revealed a significantly higher abundance of both pathways in gut microbiomes of ASCVD patients compared with controls. The present studies show a role for these two divergent microbial catalytic strategies in the meta-organismal production of PAGln. Given the numerous links between PAGln and ASCVD, these findings will assist future efforts to therapeutically target PAGln formation in vivo.

Ribonuclease T2 represents a distinct circularly permutated version of the BECR RNases

Detection of homologous relationships among proteins and understanding their mechanisms of diversification are major topics in the fields of protein science, bioinformatics, and phylogenetics. Recent developments in sequence/profile-based and structural similarity-based methods have greatly facilitated the unification and classification of many protein families into superfamilies or folds, yet many proteins remain unclassified in current protein databases. As one of the three earliest identified RNases in biology, ribonuclease T2, also known as RNase I in Escherichia coli, RNase Rh in fungi, or S-RNase in plant, is thought to be an ancient RNase family due to its widespread distribution and distinct structure. In this study, we present evidence that RNase T2 represents a circularly permutated version of the BECR (Barnase-EndoU-Colicin E5/D-RelE) fold RNases. This subtle relationship cannot be detected by traditional methods such as sequence/profile-based comparisons, structure-similarity searches, and circular permutation detections. However, we were able to identify the structural similarity using rational reconstruction of a theoretical RNase T2 ancestor via a reverse circular permutation process, followed by structural modeling using AlphaFold2, and structural comparisons. This relationship is further supported by the fact that RNase T2 and other typical BECR RNases, namely Colicin D, RNase A, and BrnT, share similar catalytic site configurations, all involving an analogous set of conserved residues on the α0 helix and the β4 strand of the BECR fold. This study revealed a hidden root of RNase T2 in bacterial toxin systems and demonstrated that reconstruction and modeling of ancestral topology is an effective strategy to identify remote relationship between proteins.

Identifying Novel Inhibitors for Hepatic Organic Anion Transporting Polypeptides by Machine Learning-Based Virtual Screening

Integration of statistical learning methods with structure-based modeling approaches is a contemporary strategy to identify novel lead compounds in drug discovery. Hepatic organic anion transporting polypeptides (OATP1B1, OATP1B3, and OATP2B1) are classical off-targets, and it is well recognized that their ability to interfere with a wide range of chemically unrelated drugs, environmental chemicals, or food additives can lead to unwanted adverse effects like liver toxicity and drug-drug or drug-food interactions. Therefore, the identification of novel (tool) compounds for hepatic OATPs by virtual screening approaches and subsequent experimental validation is a major asset for elucidating structure-function relationships of (related) transporters: they enhance our understanding about molecular determinants and structural aspects of hepatic OATPs driving ligand binding and selectivity. In the present study, we performed a consensus virtual screening approach by using different types of machine learning models (proteochemometric models, conformal prediction models, and XGBoost models for hepatic OATPs), followed by molecular docking of preselected hits using previously established structural models for hepatic OATPs. Screening the diverse REAL drug-like set (Enamine) shows a comparable hit rate for OATP1B1 (36% actives) and OATP1B3 (32% actives), while the hit rate for OATP2B1 was even higher (66% actives). Percentage inhibition values for 44 selected compounds were determined using dedicated in vitro assays and guided the prioritization of several highly potent novel hepatic OATP inhibitors: six (strong) OATP2B1 inhibitors (IC₅₀ values ranging from 0.04 to 6 μM), three OATP1B1 inhibitors (2.69 to 10 μM), and five OATP1B3 inhibitors (1.53 to 10 μM) were identified. Strikingly, two novel OATP2B1 inhibitors were uncovered (C7 and H5) which show high affinity (IC₅₀ values: 40 nM and 390 nM) comparable to the recently described estrone-based inhibitor (IC₅₀ = 41 nM). A molecularly detailed explanation for the observed differences in ligand binding to the three transporters is given by means of structural comparison of the detected binding sites and docking poses.

Bespoke library docking for 5-HT2A receptor agonists with antidepressant activity

There is considerable interest in screening ultralarge chemical libraries for ligand discovery, both empirically and computationally1-4. Efforts have focused on readily synthesizable molecules, inevitably leaving many chemotypes unexplored. Here we investigate structure-based docking of a bespoke virtual library of tetrahydropyridines-a scaffold that is poorly sampled by a general billion-molecule virtual library but is well suited to many aminergic G-protein-coupled receptors. Using three inputs, each with diverse available derivatives, a one pot C-H alkenylation, electrocyclization and reduction provides the tetrahydropyridine core with up to six sites of derivatization5-7. Docking a virtual library of 75 million tetrahydropyridines against a model of the serotonin 5-HT2A receptor (5-HT2AR) led to the synthesis and testing of 17 initial molecules. Four of these molecules had low-micromolar activities against either the 5-HT2A or the 5-HT2B receptors. Structure-based optimization led to the 5-HT2AR agonists (R)-69 and (R)-70, with half-maximal effective concentration values of 41 nM and 110 nM, respectively, and unusual signalling kinetics that differ from psychedelic 5-HT2AR agonists. Cryo-electron microscopy structural analysis confirmed the predicted binding mode to 5-HT2AR. The favourable physical properties of these new agonists conferred high brain permeability, enabling mouse behavioural assays. Notably, neither had psychedelic activity, in contrast to classic 5-HT2AR agonists, whereas both had potent antidepressant activity in mouse models and had the same efficacy as antidepressants such as fluoxetine at as low as 1/40th of the dose. Prospects for using bespoke virtual libraries to sample pharmacologically relevant chemical space will be considered.

A multidomain connector links the outer membrane and cell wall in phylogenetically deep-branching bacteria

Deinococcus radiodurans is a phylogenetically deep-branching extremophilic bacterium that is remarkably tolerant to numerous environmental stresses, including large doses of ultraviolet (UV) radiation and extreme temperatures. It can even survive in outer space for several years. This endurance of D. radiodurans has been partly ascribed to its atypical cell envelope comprising an inner membrane, a large periplasmic space with a thick peptidoglycan (PG) layer, and an outer membrane (OM) covered by a surface layer (S-layer). Despite intense research, molecular principles governing envelope organization and OM stabilization are unclear in D. radiodurans and related bacteria. Here, we report a electron cryomicroscopy (cryo-EM) structure of the abundant D. radiodurans OM protein SlpA, showing how its C-terminal segment forms homotrimers of 30-stranded β-barrels in the OM, whereas its N-terminal segment forms long, homotrimeric coiled coils linking the OM to the PG layer via S-layer homology (SLH) domains. Furthermore, using protein structure prediction and sequence-based bioinformatic analysis, we show that SlpA-like putative OM-PG connector proteins are widespread in phylogenetically deep-branching Gram-negative bacteria. Finally, combining our atomic structures with fluorescence and electron microscopy of cell envelopes of wild-type and mutant bacterial strains, we report a model for the cell surface of D. radiodurans. Our results will have important implications for understanding the cell surface organization and hyperstability of D. radiodurans and related bacteria and the evolutionary transition between Gram-negative and Gram-positive bacteria.

Global protein dynamics as communication sensors in peptide synthetase domains

Biological activity is governed by the timely redistribution of molecular interactions, and static structural snapshots often appear insufficient to provide the molecular determinants that choreograph communication. This conundrum applies to multidomain enzymatic systems called nonribosomal peptide synthetases (NRPSs), which assemble simple substrates into complex metabolites, where a dynamic domain organization challenges rational design to produce new pharmaceuticals. Using a nuclear magnetic resonance (NMR) atomic-level readout of biochemical transformations, we demonstrate that global structural fluctuations help promote substrate-dependent communication and allosteric responses, and impeding these global dynamics by a point-site mutation hampers allostery and molecular recognition. Our results establish global structural dynamics as sensors of molecular events that can remodel domain interactions, and they provide new perspectives on mechanisms of allostery, protein communication, and NRPS synthesis.

New β-Propellers Are Continuously Amplified From Single Blades in all Major Lineages of the β-Propeller Superfamily

β-Propellers are toroidal folds, in which consecutive supersecondary structure units of four anti-parallel β-strands-called blades-are arranged radially around a central axis. Uniquely among toroidal folds, blades span the full range of sequence symmetry, from near identity to complete divergence, indicating an ongoing process of amplification and differentiation. We have proposed that the major lineages of β-propellers arose through this mechanism and that therefore their last common ancestor was a single blade, not a fully formed β-propeller. Here we show that this process of amplification and differentiation is also widespread within individual lineages, yielding β-propellers with blades of more than 60% pairwise sequence identity in most major β-propeller families. In some cases, the blades are nearly identical, indicating a very recent amplification event, but even in cases where such recently amplified β-propellers have more than 80% overall sequence identity to each other, comparison of their DNA sequence shows that the amplification occurred independently.

Probing the structure and function of polymerase θ helicase-like domain

DNA Polymerase θ is the key actuator of the recently identified double-strand break repair pathway, theta-mediated end joining (TMEJ). It is the only known polymerase to have a 3-domain architecture containing an independently functional family A DNA polymerase tethered by a long central region to an N-terminal helicase-like domain (HLD). Full-length polymerase θ and the isolated HLD hydrolyze ATP in the presence of DNA, but no processive DNA duplex unwinding has been observed. Based on sequence and structure conservation, the HLD is classified as a member of helicase superfamily II and, more specifically, the Ski2-like family. The specific subdomain composition and organization most closely resemble that of archaeal DNA repair helicases Hel308 and Hjm. The underlying structural basis as to why the HLD is not able to processively unwind duplex DNA, despite its similarity to bona fide helicases, remains elusive. Activities of the HLD include ATP hydrolysis, protein displacement, and annealing of complementary DNA. These observations have led to speculation about the role of the HLD within the context of double-strand break repair via TMEJ, such as removal of single-stranded DNA binding proteins like RPA and RAD51 and microhomology alignment. This review summarizes the structural classification and organization of the polymerase θ HLD and its homologs and explores emerging data on its biochemical activities. We conclude with a simple, speculative model for the HLD's role in TMEJ.

Resurrection of an ancient inflammatory locus reveals switch to caspase-1 specificity on a caspase-4 scaffold

Pyroptosis is a mechanism of inflammatory cell death mediated by the activation of the prolytic protein gasdermin D by caspase-1, caspase-4, and caspase-5 in human, and caspase-1 and caspase-11 in mouse. In addition, caspase-1 amplifies inflammation by proteolytic activation of cytokine interleukin-1β (IL-1β). Modern mammals of the order Carnivora lack the caspase-1 catalytic domain but express an unusual version of caspase-4 that can activate both gasdermin D and IL-1β. Seeking to understand the evolutionary origin of this caspase, we utilized the large amount of data available in public databases to perform ancestral sequence reconstruction of an inflammatory caspase of a Carnivora ancestor. We expressed the catalytic domain of this putative ancestor in Escherichia coli, purified it, and compared its substrate specificity on synthetic and protein substrates to extant caspases. We demonstrated that it activates gasdermin D but has reduced ability to activate IL-1β. Our reconstruction suggests that caspase-1 was lost in a Carnivora ancestor, perhaps upon a selective pressure for which the generation of biologically active IL-1β by caspase-1 was detrimental. We speculate that later, a Carnivora encountered selective pressures that required the production of IL-1β, and caspase-4 subsequently gained this activity. This hypothesis would explain why extant Carnivora possess an inflammatory caspase with caspase-1 catalytic function placed on a caspase-4 scaffold.

Spindle-shaped archaeal viruses evolved from rod-shaped ancestors to package a larger genome

Spindle- or lemon-shaped viruses infect archaea in diverse environments. Due to the highly pleomorphic nature of these virions, which can be found with cylindrical tails emanating from the spindle-shaped body, structural studies of these capsids have been challenging. We have determined the atomic structure of the capsid of Sulfolobus monocaudavirus 1, a virus that infects hosts living in nearly boiling acid. A highly hydrophobic protein, likely integrated into the host membrane before the virions assemble, forms 7 strands that slide past each other in both the tails and the spindle body. We observe the discrete steps that occur as the tail tubes expand, and these are due to highly conserved quasiequivalent interactions with neighboring subunits maintained despite significant diameter changes. Our results show how helical assemblies can vary their diameters, becoming nearly spherical to package a larger genome and suggest how all spindle-shaped viruses have evolved from archaeal rod-like viruses.

Structure-Based Design of a Chemical Probe Set for the 5-HT5A Serotonin Receptor

The 5-HT_5A receptor (5-HT_5AR), for which no selective agonists and a few antagonists exist, remains the least understood serotonin receptor. A single commercial antagonist, SB-699551, has been widely used to investigate the 5-HT_5AR function in neurological disorders, including pain, but this molecule has substantial liabilities as a chemical probe. Accordingly, we sought to develop an internally controlled probe set. Docking over 6 million molecules against a 5-HT_5AR homology model identified 5 mid-μM ligands, one of which was optimized to UCSF678, a 42 nM arrestin-biased partial agonist at the 5-HT_5AR with a more restricted off-target profile and decreased assay liabilities versus SB-699551. Site-directed mutagenesis supported the docked pose of UCSF678. Surprisingly, analogs of UCSF678 that lost the 5-HT_5AR activity revealed that 5-HT_5AR engagement is nonessential for alleviating pain, contrary to studies with less-selective ligands. UCSF678 and analogs constitute a selective probe set with which to study the function of the 5-HT_5AR.

Molecular Insights Into The Interaction Between Human Nicotinamide Phosphoribosyltransferase and Toll-Like Receptor 4

The secreted form of the enzyme nicotinamide phosphoribosyltransferase (NAMPT), which catalyzes a key reaction in intracellular NAD biosynthesis, acts as a damage-associated molecular pattern triggering Toll-like receptor 4 (TLR4)-mediated inflammatory responses. However, the precise mechanism of interaction is unclear. Using an integrated approach combining bioinformatics and functional and structural analyses, we investigated the interaction between NAMPT and TLR4 at the molecular level. Starting from previous evidence that the bacterial ortholog of NAMPT cannot elicit the inflammatory response, despite a high degree of structural conservation, two positively charged areas unique to the human enzyme (the α1-α2 and β1-β2 loops) were identified as likely candidates for TLR4 binding. However, alanine substitution of the positively charged residues within these loops did not affect either the oligomeric state or the catalytic efficiency of the enzyme. The kinetics of the binding of wildtype and mutated NAMPT to biosensor-tethered TLR4 was analyzed. We found that mutations in the α1-α2 loop strongly decreased the association rate, increasing the K_D value from 18 nM, as determined for the wildtype, to 1.3 μM. In addition, mutations in the β1-β2 loop or its deletion increased the dissociation rate, yielding K_D values of 0.63 and 0.22 μM, respectively. Mutations also impaired the ability of NAMPT to trigger the NF-κB inflammatory signaling pathway in human cultured macrophages. Finally, the involvement of the two loops in receptor binding was supported by NAMPT-TLR4 docking simulations. This study paves the way for future development of compounds that selectively target eNAMPT/TLR4 signaling in inflammatory disorders.

Slipknot or Crystallographic Error: A Computational Analysis of the Plasmodium falciparum DHFR Structural Folds

The presence of protein structures with atypical folds in the Protein Data Bank (PDB) is rare and may result from naturally occurring knots or crystallographic errors. Proper characterisation of such folds is imperative to understanding the basis of naturally existing knots and correcting crystallographic errors. If left uncorrected, such errors can frustrate downstream experiments that depend on the structures containing them. An atypical fold has been identified in P. falciparum dihydrofolate reductase (PfDHFR) between residues 20–51 (loop 1) and residues 191–205 (loop 2). This enzyme is key to drug discovery efforts in the parasite, necessitating a thorough characterisation of these folds. Using multiple sequence alignments (MSA), a unique insert was identified in loop 1 that exacerbates the appearance of the atypical fold-giving it a slipknot-like topology. However, PfDHFR has not been deposited in the knotted proteins database, and processing its structure failed to identify any knots within its folds. The application of protein homology modelling and molecular dynamics simulations on the DHFR domain of P. falciparum and those of two other organisms (E. coli and M. tuberculosis) that were used as molecular replacement templates in solving the PfDHFR structure revealed plausible unentangled or open conformations of these loops. These results will serve as guides for crystallographic experiments to provide further insights into the atypical folds identified.

Bacterial gasdermins reveal an ancient mechanism of cell death

[Figure: see text].

Catalytically active holo Homo sapiens adenosine deaminase I adopts a closed conformation

Homo sapiens adenosine deaminase 1 (HsADA1; UniProt P00813) is an immunologically relevant enzyme with roles in T-cell activation and modulation of adenosine metabolism and signaling. Patients with genetic deficiency in HsADA1 suffer from severe combined immunodeficiency, and HsADA1 is a therapeutic target in hairy cell leukemias. Historically, insights into the catalytic mechanism and the structural attributes of HsADA1 have been derived from studies of its homologs from Bos taurus (BtADA) and Mus musculus (MmADA). Here, the structure of holo HsADA1 is presented, as well as biochemical characterization that confirms its high activity and shows that it is active across a broad pH range. Structurally, holo HsADA1 adopts a closed conformation distinct from the open conformation of holo BtADA. Comparison of holo HsADA1 and MmADA reveals that MmADA also adopts a closed conformation. These findings challenge previous assumptions gleaned from BtADA regarding the conformation of HsADA1 that may be relevant to its immunological interactions, particularly its ability to bind adenosine receptors. From a broader perspective, the structural analysis of HsADA1 presents a cautionary tale for reliance on homologs to make structural inferences relevant to applications such as protein engineering or drug development.

Genomic discovery and structural dissection of a novel type of polymorphic toxin system in gram-positive bacteria

Bacteria have developed several molecular conflict systems to facilitate kin recognition and non-kin competition to gain advantages in the acquisition of growth niches and of limited resources. One such example is a large class of so-called polymorphic toxin systems (PTSs), which comprise a variety of the toxin proteins secreted via T2SS, T5SS, T6SS, T7SS and many others. These systems are highly divergent in terms of sequence/structure, domain architecture, toxin-immunity association, and organization of the toxin loci, which makes it difficult to identify and characterize novel systems using traditional experimental and bioinformatic strategies. In recent years, we have been developing and utilizing unique genome-mining strategies and pipelines, based on the organizational principles of both domain architectures and genomic loci of PTSs, for an effective and comprehensive discovery of novel PTSs, dissection of their components, and prediction of their structures and functions. In this study, we present our systematic discovery of a new type of PTS (S8-PTS) in several gram-positive bacteria. We show that the S8-PTS contains three components: a peptidase of the S8 family (subtilases), a polymorphic toxin, and an immunity protein. We delineated the typical organization of these polymorphic toxins, in which a N-terminal signal peptide is followed by a potential receptor binding domain, BetaH, and one of 16 toxin domains. We classified each toxin domain by the distinct superfamily to which it belongs, identifying nine BECR ribonucleases, one Restriction Endonuclease, one HNH nuclease, two novel toxin domains homologous to the VOC enzymes, one toxin domain with the Frataxin-like fold, and several other unique toxin families such as Ntox33 and HicA. Accordingly, we identified 20 immunity families and classified them into different classes of folds. Further, we show that the S8-PTS-associated peptidases are analogous to many other processing peptidases found in T5SS, T7SS, T9SS, and many proprotein-processing peptidases, indicating that they function to release the toxin domains during secretion. The S8-PTSs are mostly found in animal and plant-associated bacteria, including many pathogens. We propose S8-PTSs will facilitate the competition of these bacteria with other microbes or contribute to the pathogen-host interactions.

Spontaneous Mutation in the Movement Protein of Citrus Leprosis Virus C2, in a Heterologous Virus Infection Context, Increases Cell-to-Cell Transport and Generates Fitness Advantage

Previous results using a movement defective alfalfa mosaic virus (AMV) vector revealed that citrus leprosis virus C (CiLV-C) movement protein (MP) generates a more efficient local movement, but not more systemic transport, than citrus leprosis virus C2 (CiLV-C2) MP, MPs belonging to two important viruses for the citrus industry. Here, competition experiment assays in transgenic tobacco plants (P12) between transcripts of AMV constructs expressing the cilevirus MPs, followed by several biological passages, showed the prevalence of the AMV construct carrying the CiLV-C2 MP. The analysis of AMV RNA 3 progeny recovered from P12 plant at the second viral passage revealed the presence of a mix of progeny encompassing the CiLV-C2 MP wild type (MP^WT) and two variants carrying serines instead phenylalanines at positions 72 (MP^S72F) or 259 (MP^S259F), respectively. We evaluated the effects of each modified residue in virus replication, and cell-to-cell and long-distance movements. Results indicated that phenylalanine at position 259 favors viral cell-to-cell transport with an improvement in viral fitness, but has no effect on viral replication, whereas mutation at position 72 (MP^S72F) has a penalty in the viral fitness. Our findings indicate that the prevalence of a viral population may be correlated with its greater efficiency in cell-to-cell and systemic movements.

Cytochrome P450-catalyzed biosynthesis of furanoditerpenoids in the bioenergy crop switchgrass (Panicum virgatum L.)

Specialized diterpenoid metabolites are important mediators of plant-environment interactions in monocot crops. To understand metabolite functions in plant environmental adaptation that ultimately can enable crop improvement strategies, a deeper knowledge of the underlying species-specific biosynthetic pathways is required. Here, we report the genomics-enabled discovery of five cytochrome P450 monooxygenases (CYP71Z25-CYP71Z29) that form previously unknown furanoditerpenoids in the monocot bioenergy crop Panicum virgatum (switchgrass). Combinatorial pathway reconstruction showed that CYP71Z25-CYP71Z29 catalyze furan ring addition directly to primary diterpene alcohol intermediates derived from distinct class II diterpene synthase products. Transcriptional co-expression patterns and the presence of select diterpenoids in switchgrass roots support the occurrence of P450-derived furanoditerpenoids in planta. Integrating molecular dynamics, structural analysis and targeted mutagenesis identified active site determinants that contribute to the distinct catalytic specificities underlying the broad substrate promiscuity of CYP71Z25-CYP71Z29 for native and non-native diterpenoids.

A novel Actinidia cytorhabdovirus characterized using genomic and viral protein interaction features

A novel cytorhabdovirus, tentatively named Actinidia virus D (AcVD), was identified from kiwifruit (Actinidia chinensis) in China using high-throughput sequencing technology. The genome of AcVD consists of 13,589 nucleotides and is organized into seven open reading frames (ORFs) in its antisense strand, coding for proteins in the order N-P-P3-M-G-P6-L. The ORFs were flanked by a 3' leader sequence and a 5' trailer sequence and are separated by conserved intergenic junctions. The genome sequence of AcVD was 44.6%-51.5% identical to those of reported cytorhabdoviruses. The proteins encoded by AcVD shared the highest sequence identities, ranging from 27.3% (P6) to 44.5% (L), with the respective proteins encoded by reported cytorhabdoviruses. Phylogenetic analysis revealed that AcVD clustered together with the cytorhabdovirus Wuhan insect virus 4. The subcellular locations of the viral proteins N, P, P3, M, G, and P6 in epidermal cells of Nicotiana benthamiana leaves were determined. The M protein of AcVD uniquely formed filament structures and was associated with microtubules. Bimolecular fluorescence complementation assays showed that three proteins, N, P, and M, self-interact, protein N plays a role in the formation of cytoplasm viroplasm, and protein M recruits N, P, P3, and G to microtubules. In addition, numerous paired proteins interact in the nucleus. This study presents the first evidence of a cytorhabdovirus infecting kiwifruit plants and full location and interaction maps to gain insight into viral protein functions.

Differences in Gluco and Galacto Substrate-Binding Interactions in a Dual 6Pβ-Glucosidase/6Pβ-Galactosidase Glycoside Hydrolase 1 Enzyme from Bacillus licheniformis

Bacterial glycoside hydrolase 1 (GH1) enzymes with 6-phospho-β-galactosidase and 6-phospho-β-glucosidase activities have the important task of releasing phosphorylated and nonphosphorylated monosaccharides into the cytoplasm. Curiously, dual 6-phospho-β-galactosidase/6-phospho-β-glucosidase (dual-phospho) enzymes have broad specificity and are able to hydrolyze galacto- and gluco-derived substrates. This study investigates the structure and substrate specificity of a GH family 1 enzyme from Bacillus licheniformis, hereafter known as BlBglC. The enzyme structure has been solved, and sequence analysis, molecular dynamics simulations, and binding free energy calculations offered evidence of dual-phospho activity. Both test ligands p-nitrophenyl-β-d-galactoside-6-phosphate (PNP6Pgal) and p-nitrophenyl-β-d-glucoside-6-phosphate (PNP6Pglc) demonstrated strong binding to BlBglC although the pose and interactions of the PNP6Pglc triplicates were slightly more consistent. Interestingly, known specificity-inducing residues, Gln23 and Trp433, bind strongly to the ligand O3 hydroxyl group in the PNP6Pgal-BlBglC complex and to the ligand O4 hydroxyl group in the PNP6Pglc-BlBglC complex. Additionally, the BlBglC-His124 residue is a major contributor of hydrogen bonds to the PNP6Pgal O3 hydroxyl group but does not form any hydrogen bonds with PNP6Pglc. On the other hand, BlBglC residues Tyr173, Tyr301, Gln302, and Thr321 form hydrogen bonds with PNP6Pglc but not PNP6Pgal. These findings provide important details of the broad specificity of dual-phospho activity GH1 enzymes.

Proteomic Tools for the Analysis of Cytoskeleton Proteins

Proteomic analyses have become an essential part of the toolkit of the molecular biologist, given the widespread availability of genomic data and open source or freely accessible bioinformatics software. Tools are available for detecting homologous sequences, recognizing functional domains, and modeling the three-dimensional structure for any given protein sequence, as well as for predicting interactions with other proteins or macromolecules. Although a wealth of structural and functional information is available for many cytoskeletal proteins, with representatives spanning all of the major subfamilies, the majority of cytoskeletal proteins remain partially or totally uncharacterized. Moreover, bioinformatics tools provide a means for studying the effects of synthetic mutations or naturally occurring variants of these cytoskeletal proteins. This chapter discusses various freely available proteomic analysis tools, with a focus on in silico prediction of protein structure and function. The selected tools are notable for providing an easily accessible interface for the novice while retaining advanced functionality for more experienced computational biologists.

Recombinant antibodies recognize conformation-dependent epitopes of the leucine zipper of misfolding-prone myocilin

Recombinant antibodies with well-characterized epitopes and known conformational specificities are critical reagents to support robust interpretation and reproducibility of immunoassays across biomedical research. For myocilin, a protein prone to misfolding that is associated with glaucoma and an emerging player in other human diseases, currently available antibodies are unable to differentiate among the numerous disease-associated protein states. This fundamentally constrains efforts to understand the connection between myocilin structure, function, and disease. To address this concern, we used protein engineering methods to develop new recombinant antibodies that detect the N-terminal leucine zipper structural domain of myocilin and that are cross-reactive for human and mouse myocilin. After harvesting spleens from immunized mice and in vitro library panning, we identified two antibodies, 2A4 and 1G12. 2A4 specifically recognizes a folded epitope while 1G12 recognizes a range of conformations. We matured antibody 2A4 for improved biophysical properties, resulting in variant 2H2. In a human IgG1 format, 2A4, 1G12, and 2H2 immunoprecipitate full-length folded myocilin present in the spent media of human trabecular meshwork (TM) cells, and 2H2 can visualize myocilin in fixed human TM cells using fluorescence microscopy. These new antibodies should find broad application in glaucoma and other research across multiple species platforms.

A novel heteromeric pantothenate kinase complex in apicomplexan parasites

Coenzyme A is synthesised from pantothenate via five enzyme-mediated steps. The first step is catalysed by pantothenate kinase (PanK). All PanKs characterised to date form homodimers. Many organisms express multiple PanKs. In some cases, these PanKs are not functionally redundant, and some appear to be non-functional. Here, we investigate the PanKs in two pathogenic apicomplexan parasites, Plasmodium falciparum and Toxoplasma gondii. Each of these organisms express two PanK homologues (PanK1 and PanK2). We demonstrate that PfPanK1 and PfPanK2 associate, forming a single, functional PanK complex that includes the multi-functional protein, Pf14-3-3I. Similarly, we demonstrate that TgPanK1 and TgPanK2 form a single complex that possesses PanK activity. Both TgPanK1 and TgPanK2 are essential for T. gondii proliferation, specifically due to their PanK activity. Our study constitutes the first examples of heteromeric PanK complexes in nature and provides an explanation for the presence of multiple PanKs within certain organisms.

ProteoVision: web server for advanced visualization of ribosomal proteins

ProteoVision is a web server designed to explore protein structure and evolution through simultaneous visualization of multiple sequence alignments, topology diagrams and 3D structures. Starting with a multiple sequence alignment, ProteoVision computes conservation scores and a variety of physicochemical properties and simultaneously maps and visualizes alignments and other data on multiple levels of representation. The web server calculates and displays frequencies of amino acids. ProteoVision is optimized for ribosomal proteins but is applicable to analysis of any protein. ProteoVision handles internally generated and user uploaded alignments and connects them with a selected structure, found in the PDB or uploaded by the user. It can generate de novo topology diagrams from three-dimensional structures. All displayed data is interactive and can be saved in various formats as publication quality images or external datasets or PyMol Scripts. ProteoVision enables detailed study of protein fragments defined by Evolutionary Classification of protein Domains (ECOD) classification. ProteoVision is available at http://proteovision.chemistry.gatech.edu/.

Early Metazoan Origin and Multiple Losses of a Novel Clade of RIM Presynaptic Calcium Channel Scaffolding Protein Homologs

The precise localization of Ca_V2 voltage-gated calcium channels at the synapse active zone requires various interacting proteins, of which, Rab3-interacting molecule or RIM is considered particularly important. In vertebrates, RIM interacts with Ca_V2 channels in vitro via a PDZ domain that binds to the extreme C-termini of the channels at acidic ligand motifs of D/E-D/E/H-WC-_COOH, and knockout of RIM in vertebrates and invertebrates disrupts Ca_V2 channel synaptic localization and synapse function. Here, we describe a previously uncharacterized clade of RIM proteins bearing domain architectures homologous to those of known RIM homologs, but with some notable differences including key amino acids associated with PDZ domain ligand specificity. This novel RIM emerged near the stem lineage of metazoans and underwent extensive losses, but is retained in select animals including the early-diverging placozoan Trichoplax adhaerens, and molluscs. RNA expression and localization studies in Trichoplax and the mollusc snail Lymnaea stagnalis indicate differential regional/tissue type expression, but overlapping expression in single isolated neurons from Lymnaea. Ctenophores, the most early-diverging animals with synapses, are unique among animals with nervous systems in that they lack the canonical RIM, bearing only the newly identified homolog. Through phylogenetic analysis, we find that Ca_V2 channel D/E-D/E/H-WC-_COOH like PDZ ligand motifs were present in the common ancestor of cnidarians and bilaterians, and delineate some deeply conserved C-terminal structures that distinguish Ca_V1 from Ca_V2 channels, and Ca_V1/Ca_V2 from Ca_V3 channels.

A Mechanically Transmitted DNA Mycovirus Is Targeted by the Defence Machinery of Its Host, Botrytis cinerea

Eukaryotic circular single-stranded DNA (ssDNA) viruses were known only to infect plants and vertebrates until the discovery of the isolated DNA mycovirus from the fungus Sclerotinia sclerotiorum. Similar viral sequences were reported from several other sources and classified in ten genera within the Genomoviridae family. The current study reports two circular ssDNA mycoviruses isolated from the phytopathogen Botrytis cinerea, and their assignment to a newly created genus tentatively named Gemydayirivirus. The mycoviruses, tentatively named botrytis gemydayirivirus 1 (BGDaV1) and BGDaV2, are 1701 and 1693 nt long and encode three and two open reading frames (ORFs), respectively. Of the predicted ORFs, only ORF I, which codes for a replication initiation protein (Rep), shared identity with other proteins in GenBank. BGDaV1 is infective as cell-free purified particles and confers hypovirulence on its natural host. Investigation revealed that BGDaV1 is a target for RNA silencing and genomic DNA methylation, keeping the virus at very low titre. The discovery of BGDaV1 expands our knowledge of the diversity of genomoviruses and their interaction with fungal hosts.

b2bTools: online predictions for protein biophysical features and their conservation

We provide integrated protein sequence-based predictions via https://bio2byte.be/b2btools/. The aim of our predictions is to identify the biophysical behaviour or features of proteins that are not readily captured by structural biology and/or molecular dynamics approaches. Upload of a FASTA file or text input of a sequence provides integrated predictions from DynaMine backbone and side-chain dynamics, conformational propensities, and derived EFoldMine early folding, DisoMine disorder, and Agmata β-sheet aggregation. These predictions, several of which were previously not available online, capture 'emergent' properties of proteins, i.e. the inherent biophysical propensities encoded in their sequence, rather than context-dependent behaviour (e.g. final folded state). In addition, upload of a multiple sequence alignment (MSA) in a variety of formats enables exploration of the biophysical variation observed in homologous proteins. The associated plots indicate the biophysical limits of functionally relevant protein behaviour, with unusual residues flagged by a Gaussian mixture model analysis. The prediction results are available as JSON or CSV files and directly accessible via an API. Online visualisation is available as interactive plots, with brief explanations and tutorial pages included. The server and API employ an email-free token-based system that can be used to anonymously access previously generated results.

Data-Driven Ensemble Docking to Map Molecular Interactions of Steroid Analogs with Hepatic Organic Anion Transporting Polypeptides

Hepatic organic anion transporting polypeptides—OATP1B1, OATP1B3, and OATP2B1—are expressed at the basolateral membrane of hepatocytes, being responsible for the uptake of a wide range of natural substrates and structurally unrelated pharmaceuticals. Impaired function of hepatic OATPs has been linked to clinically relevant drug–drug interactions leading to altered pharmacokinetics of administered drugs. Therefore, understanding the commonalities and differences across the three transporters represents useful knowledge to guide the drug discovery process at an early stage. Unfortunately, such efforts remain challenging because of the lack of experimentally resolved protein structures for any member of the OATP family. In this study, we established a rigorous computational protocol to generate and validate structural models for hepatic OATPs. The multistep procedure is based on the systematic exploration of available protein structures with shared protein folding using normal-mode analysis, the calculation of multiple template backbones from elastic network models, the utilization of multiple template conformations to generate OATP structural models with various degrees of conformational flexibility, and the prioritization of models on the basis of enrichment docking. We employed the resulting OATP models of OATP1B1, OATP1B3, and OATP2B1 to elucidate binding modes of steroid analogs in the three transporters. Steroid conjugates have been recognized as endogenous substrates of these transporters. Thus, investigating this data set delivers insights into mechanisms of substrate recognition. In silico predictions were complemented with in vitro studies measuring the bioactivity of a compound set on OATP expressing cell lines. Important structural determinants conferring shared and distinct binding patterns of steroid analogs in the three transporters have been identified. Overall, this comparative study provides novel insights into hepatic OATP-ligand interactions and selectivity. Furthermore, the integrative computational workflow for structure-based modeling can be leveraged for other pharmaceutical targets of interest.

Distribution of Merlin in eukaryotes and first report of DNA transposons in kinetoplastid protists

DNA transposons are defined as repeated DNA sequences that can move within the host genome through the action of transposases. The transposon superfamily Merlin was originally found mainly in animal genomes. Here, we describe a global distribution of the Merlin in animals, fungi, plants and protists, reporting for the first time their presence in Rhodophyceae, Metamonada, Discoba and Alveolata. We identified a great variety of potentially active Merlin families, some containing highly imperfect terminal inverted repeats and internal tandem repeats. Merlin-related sequences with no evidence of mobilization capacity were also observed and may be products of domestication. The evolutionary trees support that Merlin is likely an ancient superfamily, with early events of diversification and secondary losses, although repeated re-invasions probably occurred in some groups, which would explain its diversity and discontinuous distribution. We cannot rule out the possibility that the Merlin superfamily is the product of multiple horizontal transfers of related prokaryotic insertion sequences. Moreover, this is the first account of a DNA transposon in kinetoplastid flagellates, with conserved Merlin transposase identified in Bodo saltans and Perkinsela sp., whereas it is absent in trypanosomatids. Based on the level of conservation of the transposase and overlaps of putative open reading frames with Merlin, we propose that in protists it may serve as a raw material for gene emergence.

Structural characterization of a PCP-R didomain from an archaeal nonribosomal peptide synthetase reveals novel interdomain interactions

Nonribosomal peptide synthetases (NRPSs) are multimodular enzymes that produce a wide range of bioactive peptides, such as siderophores, toxins, and antibacterial and insecticidal agents. NRPSs are dynamic proteins characterized by extensive interdomain communications as a consequence of their assembly-line mode of synthesis. Hence, crystal structures of multidomain fragments of NRPSs have aided in elucidating crucial interdomain interactions that occur during different steps of the NRPS catalytic cycle. One crucial yet unexplored interaction is that between the reductase (R) domain and the peptide carrier protein (PCP) domain. R domains are members of the short-chain dehydrogenase/reductase family and function as termination domains that catalyze the reductive release of the final peptide product from the terminal PCP domain of the NRPS. Here, we report the crystal structure of an archaeal NRPS PCP-R didomain construct. This is the first NRPS R domain structure to be determined together with the upstream PCP domain and is also the first structure of an archaeal NRPS to be reported. The structure reveals that a novel helix-turn-helix motif, found in NRPS R domains but not in other short-chain dehydrogenase/reductase family members, plays a major role in the interface between the PCP and R domains. The information derived from the described PCP-R interface will aid in gaining further mechanistic insights into the peptide termination reaction catalyzed by the R domain and may have implications in engineering NRPSs to synthesize novel peptide products.

cGAS-like receptors sense RNA and control 3'2'-cGAMP signalling in Drosophila

Cyclic GMP-AMP synthase (cGAS) is a cytosolic DNA sensor that produces the second messenger cG[2'-5']pA[3'-5']p (2'3'-cGAMP) and controls activation of innate immunity in mammalian cells^1-5. Animal genomes typically encode multiple proteins with predicted homology to cGAS^6-10, but the function of these uncharacterized enzymes is unknown. Here we show that cGAS-like receptors (cGLRs) are innate immune sensors that are capable of recognizing divergent molecular patterns and catalysing synthesis of distinct nucleotide second messenger signals. Crystal structures of human and insect cGLRs reveal a nucleotidyltransferase signalling core shared with cGAS and a diversified primary ligand-binding surface modified with notable insertions and deletions. We demonstrate that surface remodelling of cGLRs enables altered ligand specificity and used a forward biochemical screen to identify cGLR1 as a double-stranded RNA sensor in the model organism Drosophila melanogaster. We show that RNA recognition activates Drosophila cGLR1 to synthesize the novel product cG[3'-5']pA[2'-5']p (3'2'-cGAMP). A crystal structure of Drosophila stimulator of interferon genes (dSTING) in complex with 3'2'-cGAMP explains selective isomer recognition, and 3'2'-cGAMP induces an enhanced antiviral state in vivo that protects from viral infection. Similar to radiation of Toll-like receptors in pathogen immunity, our results establish cGLRs as a diverse family of metazoan pattern recognition receptors.