PROMALS3D: a tool for multiple protein sequence and structure alignments

Two functionally distinct heme/iron transport systems are virulence determinants of the fish pathogen Flavobacterium psychrophilum

Bacterial pathogens have a critical impact on aquaculture, a sector that accounts for half of the human fish consumption. Flavobacterium psychrophilum (phylum Bacteroidetes) is responsible for bacterial cold-water disease in salmonids worldwide. The molecular factors involved in host invasion, colonization and haemorrhagic septicaemia are mostly unknown. In this study, we identified two new TonB-dependent receptors, HfpR and BfpR, that are required for adaptation to iron conditions encountered during infection and for virulence in rainbow trout. Transcriptional analyses revealed that their expression is tightly controlled and upregulated under specific iron sources and concentrations. Characterization of deletion mutants showed that they act without redundancy: BfpR is required for optimal growth in the presence of high haemoglobin level, while HfpR confers the capacity to acquire nutrient iron from haem or haemoglobin under iron scarcity. The gene hfpY, co-transcribed with hfpR, encodes a protein related to the HmuY family. We demonstrated that HfpY binds haem and contributes significantly to host colonization and disease severity. Overall, these results are consistent with a model in which both BfpR and Hfp systems promote haem uptake and respond to distinct signals to adapt iron acquisition to the different stages of pathogenesis. Our findings give insight into the molecular basis of pathogenicity of a serious pathogen belonging to the understudied family Flavobacteriaceae and point to the newly identified haem receptors as promising targets for antibacterial development.

A novel statovirus identified in fecal samples from wild geladas in the Ethiopian highlands

We present a novel statovirus in geladas (Theropithecus gelada), graminivorous primates endemic to the Ethiopian highlands. Using a high-throughput sequencing approach, we identified contiguous sequences in feces from two adult female geladas in the Simien Mountains National Park, Ethiopia, that share similarities to statoviruses. Our phylogenetic analysis of the whole genome, as well as the RNA-dependent RNA polymerase (RdRp) and capsid protein (CP) amino acid sequences, revealed that the gelada statoviruses cluster with those from other primates (laboratory populations of Macaca nemestrina and Macaca mulatta). As the first report of statovirus in wild primates, this finding contributes to our understanding of the phylogenetic and geographic distribution of statoviruses and their hosts.

A microbial transporter of the dietary antioxidant ergothioneine

Low-molecular-weight (LMW) thiols are small-molecule antioxidants required for the maintenance of intracellular redox homeostasis. However, many host-associated microbes, including the gastric pathogen Helicobacter pylori, unexpectedly lack LMW-thiol biosynthetic pathways. Using reactivity-guided metabolomics, we identified the unusual LMW thiol ergothioneine (EGT) in H. pylori. Dietary EGT accumulates to millimolar levels in human tissues and has been broadly implicated in mitigating disease risk. Although certain microorganisms synthesize EGT, we discovered that H. pylori acquires this LMW thiol from the host environment using a highly selective ATP-binding cassette transporter-EgtUV. EgtUV confers a competitive colonization advantage in vivo and is widely conserved in gastrointestinal microbes. Furthermore, we found that human fecal bacteria metabolize EGT, which may contribute to production of the disease-associated metabolite trimethylamine N-oxide. Collectively, our findings illustrate a previously unappreciated mechanism of microbial redox regulation in the gut and suggest that inter-kingdom competition for dietary EGT may broadly impact human health.

TENT5 cytoplasmic noncanonical poly(A) polymerases regulate the innate immune response in animals

Innate immunity is the first line of host defense against pathogens. Here, through global transcriptome and proteome analyses, we uncover that newly described cytoplasmic poly(A) polymerase TENT-5 (terminal nucleotidyltransferase 5) enhances the expression of secreted innate immunity effector proteins in Caenorhabditis elegans. Direct RNA sequencing revealed that multiple mRNAs with signal peptide-encoding sequences have shorter poly(A) tails in tent-5-deficient worms. Those mRNAs are translated at the endoplasmic reticulum where a fraction of TENT-5 is present, implying that they represent its direct substrates. Loss of tent-5 makes worms more susceptible to bacterial infection. Notably, the role of TENT-5 in innate immunity is evolutionarily conserved. Its orthologs, TENT5A and TENT5C, are expressed in macrophages and induced during their activation. Analysis of macrophages devoid of TENT5A/C revealed their role in the regulation of secreted proteins involved in defense response. In summary, our study reveals cytoplasmic polyadenylation to be a previously unknown component of the posttranscriptional regulation of innate immunity in animals.

Insights into the role of Nup62 and Nup93 in assembling cytoplasmic ring and central transport channel of the nuclear pore complex

The nuclear pore complex (NPC) is a highly modular assembly of 34 distinct nucleoporins (Nups) to form a versatile transport channel between the nucleus and the cytoplasm. Among them, Nup62 is known as an essential component for nuclear transport, Nup93 for proper nuclear envelope assembly. These Nups constitute various NPC subcomplexes such as the central transport channel (CTC), the cytoplasmic ring (CR), and the inner ring (IR). However, how they play their roles in NPC assembly and transport activity is not clear. Here we delineated the interacting regions and conducted biochemical reconstitution and structural characterization of the mammalian CR complex to reveal its intrinsic dynamic behavior and a distinct "4"-shaped architecture resembling the CTC complex. Our in vitro reconstitution data demonstrate that the Nup62 coiled-coil domain is critical to form both Nup62^322-525 •Nup88^517-742 and Nup62^322-525•Nup88^517-742•Nup214^693-926 heterotrimers and both can bind to Nup93^1-150. We therefore propose that Nup93 acts as a "sensor" to bind to Nup62 shared heterotrimers including the Nup62•Nup54 heterotrimer of the CTC, which was not shown previously to be an interacting partner. Altogether, our biochemical study suggests that Nup62 via its coiled-coil domain is central to form compositionally distinct yet structurally similar heterotrimers and Nup93 binds these diverse heterotrimers nonselectively.

Structural Homology-Based Drug Repurposing Approach for Targeting NSP12 SARS-CoV-2

The severe acute respiratory syndrome coronavirus 2, also known as SARS-CoV-2, is the causative agent of the COVID-19 global pandemic. SARS-CoV-2 has a highly conserved non-structural protein 12 (NSP-12) involved in RNA-dependent RNA polymerase (RdRp) activity. For the identification of potential inhibitors for NSP-12, computational approaches such as the identification of homologous proteins that have been previously targeted by FDA-approved antivirals can be employed. Herein, homologous proteins of NSP-12 were retrieved from Protein DataBank (PDB) and the evolutionary conserved sequence and structure similarity of the active site of the RdRp domain of NSP-12 was characterized. The identified homologous structures of NSP-12 belonged to four viral families: Coronaviridae, Flaviviridae, Picornaviridae, and Caliciviridae, and shared evolutionary conserved relationships. The multiple sequences and structural alignment of homologous structures showed highly conserved amino acid residues that were located at the active site of the RdRp domain of NSP-12. The conserved active site of the RdRp domain of NSP-12 was evaluated for binding affinity with the FDA-approved antivirals, i.e., Sofosbuvir and Dasabuvir in a molecular docking study. The molecular docking of Sofosbuvir and Dasabuvir with the active site that contains conserved motifs (motif A-G) of the RdRp domain of NSP-12 revealed significant binding affinity. Furthermore, MD simulation also inferred the potency of Sofosbuvir and Dasabuvir. In conclusion, targeting the active site of the RdRp domain of NSP-12 with Dasabuvir and Sofosbuvir might reduce viral replication and pathogenicity and could be further studied for the treatment of SARS-CoV-2.

The Rhizobium tropici CIAT 899 NodD2 protein promotes symbiosis and extends rhizobial nodulation range by constitutive nodulation factor synthesis

In the symbiotic associations between rhizobia and legumes, the NodD regulators orchestrate the transcription of the specific nodulation genes. This set of genes is involved in the synthesis of nodulation factors, which are responsible for initiating the nodulation process. Rhizobium tropici CIAT 899 is the most successful symbiont of Phaseolus vulgaris and can nodulate a variety of legumes. Among the five NodD regulators present in this rhizobium, only NodD1 and NodD2 seem to have a role in the symbiotic process. However, the individual role of each NodD in the absence of the other proteins has remained elusive. In this work, we show that the CIAT 899 NodD2 does not require activation by inducers to promote the synthesis of nodulation factors. A CIAT 899 strain overexpressing nodD2, but lacking all additional nodD genes, can nodulate three different legumes as efficiently as the wild type. Interestingly, CIAT 899 NodD2-mediated gain of nodulation can be extended to another rhizobial species, since its overproduction in Sinorhizobium fredii HH103 not only increases the number of nitrogen-fixing nodules in two host legumes but also results in nodule development in incompatible legumes. These findings potentially open exciting opportunities to develop rhizobial inoculants and increase legume crop production.

Non-canonical Staphylococcus aureus pathogenicity island repression

Mobile genetic elements control their life cycles by the expression of a master repressor, whose function must be disabled to allow the spread of these elements in nature. Here, we describe an unprecedented repression-derepression mechanism involved in the transfer of Staphylococcus aureus pathogenicity islands (SaPIs). Contrary to the classical phage and SaPI repressors, which are dimers, the SaPI1 repressor StlSaPI1 presents a unique tetrameric conformation never seen before. Importantly, not just one but two tetramers are required for SaPI1 repression, which increases the novelty of the system. To derepress SaPI1, the phage-encoded protein Sri binds to and induces a conformational change in the DNA binding domains of StlSaPI1, preventing the binding of the repressor to its cognate StlSaPI1 sites. Finally, our findings demonstrate that this system is not exclusive to SaPI1 but widespread in nature. Overall, our results characterize a novel repression-induction system involved in the transfer of MGE-encoded virulence factors in nature.

Phosphoproteome Dynamics of Streptomyces rimosus during Submerged Growth and Antibiotic Production

Streptomyces rimosus is an industrial streptomycete, best known as a producer of oxytetracycline, one of the most widely used antibiotics. Despite the significant contribution of Streptomyces species to the pharmaceutical industry, most omics analyses have only been conducted on the model organism Streptomyces coelicolor. In recent years, protein phosphorylation on serine, threonine, and tyrosine (Ser, Thr, and Tyr, respectively) has been shown to play a crucial role in the regulation of numerous cellular processes, including metabolic changes leading to antibiotic production and morphological changes. In this study, we performed a comprehensive quantitative (phospho)proteomic analysis during the growth of S. rimosus under conditions of oxytetracycline production and pellet fragmentation. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis combined with phosphopeptide enrichment detected a total of 3,725 proteins, corresponding to 45.6% of the proteome and 417 phosphorylation sites from 230 phosphoproteins. Significant changes in abundance during three distinct growth phases were determined for 494 proteins and 98 phosphorylation sites. Functional analysis revealed changes in phosphorylation events of proteins involved in important cellular processes, including regulatory mechanisms, primary and secondary metabolism, cell division, and stress response. About 80% of the phosphoproteins detected during submerged growth of S. rimosus have not yet been reported in streptomycetes, and 55 phosphoproteins were not reported in any prokaryote studied so far. This enabled the creation of a unique resource that provides novel insights into the dynamics of (phospho)proteins and reveals many potential regulatory events during antibiotic production in liquid culture of an industrially important bacterium. IMPORTANCE Streptomyces rimosus is best known as a primary source of oxytetracycline (OTC). The significant global market value of OTC highlights the need for a better understanding of the regulatory mechanisms that lead to production of this antibiotic. Our study provides, for the first time, a detailed insight into the dynamics of (phospho)proteomic profiles during growth and antibiotic production in liquid culture of S. rimosus. Significant changes in protein synthesis and phosphorylation have been revealed for a number of important cellular proteins during the growth stages that coincide with OTC production and morphological changes of this industrially important bacterium. Most of these proteins have not been detected in previous studies. Therefore, our results significantly expand the insight into phosphorylation events associated with important cellular processes and antibiotic production; they also greatly increase the phosphoproteome of streptomycetes and contribute with newly discovered phosphoproteins to the database of prokaryotic phosphoproteomes. This can consequently lead to the design of novel research directions in elucidation of the complex regulatory network in Streptomyces.

Structural and functional investigation of ABC transporter STE6-2p from Pichia pastoris reveals unexpected interaction with sterol molecules

Adenosine triphosphate (ATP)-binding cassette (ABC) transporters are multidomain transmembrane proteins, which facilitate the transport of various substances across cell membranes using energy derived from ATP hydrolysis. They are important drug targets since they mediate decreased drug susceptibility during pharmacological treatments. For the methylotrophic yeast Pichia pastoris, a model organism that is a widely used host for protein expression, the role and function of its ABC transporters is unexplored. In this work, we investigated the Pichia ABC-B transporter STE6-2p. Functional investigations revealed that STE6-2p is capable of transporting rhodamines in vivo and is active in the presence of verapamil and triazoles in vitro. A phylogenetic analysis displays homology among multidrug resistance (MDR) transporters from pathogenic fungi to human ABC-B transporters. Further, we present high-resolution single-particle electron cryomicroscopy structures of an ABC transporter from P. pastoris in the apo conformation (3.1 Å) and in complex with verapamil and adenylyl imidodiphosphate (AMP-PNP) (3.2 Å). An unknown density between transmembrane helices 4, 5, and 6 in both structures suggests the presence of a sterol-binding site of unknown function.

One fold, many functions-M23 family of peptidoglycan hydrolases

Bacterial cell walls are the guards of cell integrity. They are composed of peptidoglycan that provides rigidity to sustain internal turgor and ensures isolation from the external environment. In addition, they harbor the enzymatic machinery to secure cell wall modulations needed throughout the bacterial lifespan. The main players in this process are peptidoglycan hydrolases, a large group of enzymes with diverse specificities and different mechanisms of action. They are commonly, but not exclusively, found in prokaryotes. Although in most cases, these enzymes share the same molecular function, namely peptidoglycan hydrolysis, they are leveraged to perform a variety of physiological roles. A well-investigated family of peptidoglycan hydrolases is M23 peptidases, which display a very conserved fold, but their spectrum of lytic action is broad and includes both Gram- positive and Gram- negative bacteria. In this review, we summarize the structural, biochemical, and functional studies concerning the M23 family of peptidases based on literature and complement this knowledge by performing large-scale analyses of available protein sequences. This review has led us to gain new insight into the role of surface charge in the activity of this group of enzymes. We present relevant conclusions drawn from the analysis of available structures and indicate the main structural features that play a crucial role in specificity determination and mechanisms of latency. Our work systematizes the knowledge of the M23 family enzymes in the context of their unique antimicrobial potential against drug-resistant pathogens and presents possibilities to modulate and engineer their features to develop perfect antibacterial weapons.

Distinct Mechanisms of Target Search by Endonuclease VIII-like DNA Glycosylases

Proteins that recognize specific DNA sequences or structural elements often find their cognate DNA lesions in a processive mode, in which an enzyme binds DNA non-specifically and then slides along the DNA contour by one-dimensional diffusion. Opposite to the processive mechanism is distributive search, when an enzyme binds, samples and releases DNA without significant lateral movement. Many DNA glycosylases, the repair enzymes that excise damaged bases from DNA, use processive search to find their cognate lesions. Here, using a method based on correlated cleavage of multiply damaged oligonucleotide substrates we investigate the mechanism of lesion search by three structurally related DNA glycosylases—bacterial endonuclease VIII (Nei) and its mammalian homologs NEIL1 and NEIL2. Similarly to another homologous enzyme, bacterial formamidopyrimidine–DNA glycosylase, NEIL1 seems to use a processive mode to locate its targets. However, the processivity of Nei was notably lower, and NEIL2 exhibited almost fully distributive action on all types of substrates. Although one-dimensional diffusion is often regarded as a universal search mechanism, our results indicate that even proteins sharing a common fold may be quite different in the ways they locate their targets in DNA.

Recognition of extracellular DNA by type IV pili promotes biofilm formation by Clostridioides difficile

Clostridioides difficile is a Gram-positive bacillus, which is a frequent cause of gastrointestinal infections triggered by the depletion of the gut microbiome. Because of the frequent recurrence of these infections after antibiotic treatment, mechanisms of C. difficile persistence and recurrence, including biofilm formation, are of increasing interest. Previously, our group and others found that type IV pili, filamentous helical appendages polymerized from protein subunits, promoted microcolony and biofilm formation in C. difficile. In Gram-negative bacteria, the ability of type IV pili to mediate bacterial self-association has been explained through interactions between the pili of adjacent cells, but type IV pili from several Gram-negative species are also required for natural competence through DNA uptake. Here, we report the ability of two C. difficile pilin subunits, PilJ and PilW, to bind to DNA in vitro, as well as the defects in biofilm formation in the pilJ and pilW gene-interruption mutants. Additionally, we have resolved the X-ray crystal structure of PilW, which we use to model possible structural mechanisms for the formation of C. difficile biofilm through interactions between type IV pili and the DNA of the extracellular matrix. Taken together, our results provide further insight into the relationship between type IV pilus function and biofilm formation in C. difficile and, more broadly, suggest that DNA recognition by type IV pili and related structures may have functional importance beyond DNA uptake for natural competence.

Cryo-EM reveals a mechanism of USP1 inhibition through a cryptic binding site

Repair of DNA damage is critical to genomic integrity and frequently disrupted in cancers. Ubiquitin-specific protease 1 (USP1), a nucleus-localized deubiquitinase, lies at the interface of multiple DNA repair pathways and is a promising drug target for certain cancers. Although multiple inhibitors of this enzyme, including one in phase 1 clinical trials, have been established, their binding mode is unknown. Here, we use cryo-electron microscopy to study an assembled enzyme-substrate-inhibitor complex of USP1 and the well-established inhibitor, ML323. Achieving 2.5-Å resolution, with and without ML323, we find an unusual binding mode in which the inhibitor disrupts part of the hydrophobic core of USP1. The consequent conformational changes in the secondary structure lead to subtle rearrangements in the active site that underlie the mechanism of inhibition. These structures provide a platform for structure-based drug design targeting USP1.

Mutation in Transmembrane Domain 8 of Human Urate Transporter 1 Disrupts Uric Acid Recognition and Transport

Human urate transporter 1 (hURAT1) is the most pivotal therapeutic target for hyperuricemia. Due to a lack of crystal structure information, the atomic structure of URAT1 is not clearly understood. In this study, a multiple sequence alignment was performed, and K393, a positively charged residue in transmembrane domain (TMD) 8, was observed to be highly conserved in organic anion transporters (OATs). K393 was substituted with a positively, negatively, and neutrally charged amino acid via site-directed mutagenesis and then used to transfect HEK293 cells. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) and enzyme-linked immunosorbent assay (ELISA) analyses indicated that mutants of K393 showed mRNA and protein expression levels similar to those in the WT group. The nonpositively charged mutants K393A, K393D, and K393E eliminated 70-80% of ¹⁴C-uric acid transport capacity, while the K393H mutant showed slight and the K393R mutant showed no reduced transport capacity compared with the WT group. Binding assays indicated that K393A, K393D, and K393E conferred lowered uric acid binding affinity. As indicated by the K _m and V _max values obtained from saturation kinetic experiments, K393A, K393D, and K393E showed increased K _m values, but K393R and K393H showed K _m values similar to those in the WT group. K393 also contributed to a high affinity for benzbromarone (BM) interaction. The inhibitory effects of BM were partly abolished in K393 mutants, with increased IC₅₀ values compared with the WT group. BM also exhibited weaker inhibitory effects on ¹⁴C-uric acid binding in K393R and K393H mutants. In an outward homology model of URAT1, K393 was located in the inner part of the transport tunnel, and further molecular docking analysis indicated that uric acid and BM showed possible hydrogen bonds with K393. Mutants K393R and K393H showed possible interactions with uric acid, and positive charges confer high affinity for uric acid as revealed by their surface electrostatic potential. In conclusion, our data provide evidence that K393 is an important residue for the recognition of uric acid or inhibitors by URAT1.

Mechanism of protein-primed template-independent DNA synthesis by Abi polymerases

Abortive infection (Abi) is a bacterial antiphage defense strategy involving suicide of the infected cell. Some Abi pathways involve polymerases that are related to reverse transcriptases. They are unique in the way they combine the ability to synthesize DNA in a template-independent manner with protein priming. Here, we report crystal and cryo-electron microscopy structures of two Abi polymerases: AbiK and Abi-P2. Both proteins adopt a bilobal structure with an RT-like domain that comprises palm and fingers subdomains and a unique helical domain. AbiK and Abi-P2 adopt a hexameric and trimeric configuration, respectively, which is unprecedented for reverse transcriptases. Biochemical experiments showed that the formation of these oligomers is required for the DNA polymerization activity. The structure of the AbiK–DNA covalent adduct visualized interactions between the 3′ end of DNA and the active site and covalent attachment of the 5′ end of DNA to a tyrosine residue used for protein priming. Our data reveal a structural basis of the mechanism of highly unusual template-independent protein-priming polymerases.

Methyltransferases of Riboviria

Many viruses from the realm Riboviria infecting eukaryotic hosts encode protein domains with sequence similarity to S-adenosylmethionine-dependent methyltransferases. These protein domains are thought to be involved in methylation of the 5′-terminal cap structures in virus mRNAs. Some methyltransferase-like domains of Riboviria are homologous to the widespread cellular FtsJ/RrmJ-like methyltransferases involved in modification of cellular RNAs; other methyltransferases, found in a subset of positive-strand RNA viruses, have been assigned to a separate “Sindbis-like” family; and coronavirus-specific Nsp13/14-like methyltransferases appeared to be different from both those classes. The representative structures of proteins from all three groups belong to a specific variety of the Rossmann fold with a seven-stranded β-sheet, but it was unclear whether this structural similarity extends to the level of conserved sequence signatures. Here I survey methyltransferases in Riboviria and derive a joint sequence alignment model that covers all groups of virus methyltransferases and subsumes the previously defined conserved sequence motifs. Analysis of the spatial structures indicates that two highly conserved residues, a lysine and an aspartate, frequently contact a water molecule, which is located in the enzyme active center next to the methyl group of S-adenosylmethionine cofactor and could play a key role in the catalytic mechanism of the enzyme. Phylogenetic evidence indicates a likely origin of all methyltransferases of Riboviria from cellular RrmJ-like enzymes and their rapid divergence with infrequent horizontal transfer between distantly related viruses.

US-align: universal structure alignments of proteins, nucleic acids, and macromolecular complexes

Structure comparison and alignment are of fundamental importance in structural biology studies. We developed the first universal platform, US-align, to uniformly align monomer and complex structures of different macromolecules-proteins, RNAs and DNAs. The pipeline is built on a uniform TM-score objective function coupled with a heuristic alignment searching algorithm. Large-scale benchmarks demonstrated consistent advantages of US-align over state-of-the-art methods in pairwise and multiple structure alignments of different molecules. Detailed analyses showed that the main advantage of US-align lies in the extensive optimization of the unified objective function powered by efficient heuristic search iterations, which substantially improve the accuracy and speed of the structural alignment process. Meanwhile, the universal protocol fusing different molecular and structural types helps facilitate the heterogeneous oligomer structure comparison and template-based protein-protein and protein-RNA/DNA docking.

Structural analysis of cholesterol binding and sterol selectivity by ABCG5/G8

The ATP-binding cassette (ABC) sterol transporters are responsible for maintaining cholesterol homeostasis in mammals by participating in reverse cholesterol transport (RCT) or transintestinal cholesterol efflux (TICE). The heterodimeric ABCG5/G8 carries out selective sterol excretion, preventing the abnormal accumulation of plant sterols in human bodies, while homodimeric ABCG1 contributes to the biogenesis and metabolism of high-density lipoproteins. A sterol-binding site on ABCG5/G8 was proposed at the interface of the transmembrane domain and the core of lipid bilayers. In this study, we have determined the crystal structure of ABCG5/G8 in a cholesterol-bound state. The structure combined with amino acid sequence analysis shows that in the proximity of the sterol-binding site, a highly conserved phenylalanine array supports functional implications for ABCG cholesterol/sterol transporters. Lastly, in silico docking analysis of cholesterol and stigmasterol (a plant sterol) suggests sterol-binding selectivity on ABCG5/G8, but not ABCG1. Together, our results provide a structural basis for cholesterol binding on ABCG5/G8 and the sterol selectivity by ABCG transporters.

Prokaryotic innate immunity through pattern recognition of conserved viral proteins

Many organisms have evolved specialized immune pattern-recognition receptors, including nucleotide-binding oligomerization domain-like receptors (NLRs) of the STAND superfamily that are ubiquitous in plants, animals, and fungi. Although the roles of NLRs in eukaryotic immunity are well established, it is unknown whether prokaryotes use similar defense mechanisms. Here, we show that antiviral STAND (Avs) homologs in bacteria and archaea detect hallmark viral proteins, triggering Avs tetramerization and the activation of diverse N-terminal effector domains, including DNA endonucleases, to abrogate infection. Cryo-electron microscopy reveals that Avs sensor domains recognize conserved folds, active-site residues, and enzyme ligands, allowing a single Avs receptor to detect a wide variety of viruses. These findings extend the paradigm of pattern recognition of pathogen-specific proteins across all three domains of life.

Discovery of the First Highly Selective Antagonist of the GluK3 Kainate Receptor Subtype

Kainate receptors belong to the family of glutamate receptors ion channels, which are responsible for the majority of rapid excitatory synaptic transmission in the central nervous system. The therapeutic potential of kainate receptors is still poorly understood, which is also due to the lack of potent and subunit-selective pharmacological tools. In search of selective ligands for the GluK3 kainate receptor subtype, a series of quinoxaline-2,3-dione analogues was synthesized and pharmacologically characterized at selected recombinant ionotropic glutamate receptors. Among them, compound 28 was found to be a competitive GluK3 antagonist with submicromolar affinity and unprecedented high binding selectivity, showing a 400-fold preference for GluK3 over other homomeric receptors GluK1, GluK2, GluK5 and GluA2. Furthermore, in functional assays performed for selected metabotropic glutamate receptor subtypes, 28 did not show agonist or antagonist activity. The molecular determinants underlying the observed affinity profile of 28 were analyzed using molecular docking and molecular dynamics simulations performed for individual GluK1 and GluK3 ligand-binding domains.

Structural and dynamic mechanisms of GABAA receptor modulators with opposing activities

γ-Aminobutyric acid type A (GABA_A) receptors are pentameric ligand-gated ion channels abundant in the central nervous system and are prolific drug targets for treating anxiety, sleep disorders and epilepsy. Diverse small molecules exert a spectrum of effects on γ-aminobutyric acid type A (GABA_A) receptors by acting at the classical benzodiazepine site. They can potentiate the response to GABA, attenuate channel activity, or counteract modulation by other ligands. Structural mechanisms underlying the actions of these drugs are not fully understood. Here we present two high-resolution structures of GABA_A receptors in complex with zolpidem, a positive allosteric modulator and heavily prescribed hypnotic, and DMCM, a negative allosteric modulator with convulsant and anxiogenic properties. These two drugs share the extracellular benzodiazepine site at the α/γ subunit interface and two transmembrane sites at β/α interfaces. Structural analyses reveal a basis for the subtype selectivity of zolpidem that underlies its clinical success. Molecular dynamics simulations provide insight into how DMCM switches from a negative to a positive modulator as a function of binding site occupancy. Together, these findings expand our understanding of how GABA_A receptor allosteric modulators acting through a common site can have diverging activities.

GABA_A receptors are important targets for anxiety, sedation and anesthesia. Here, the authors present structures bound by zolpidem (Ambien), the most prescribed hypnotic in the US, and DMCM, a negative modulator, providing insights into receptor modulation.

Atypical integrative element with strand-biased circularization activity assists interspecies antimicrobial resistance gene transfer from Vibrio alfacsensis

The exchange of antimicrobial resistance (AMR) genes between aquaculture and terrestrial microbial populations has emerged as a serious public health concern. However, the nature of the mobile genetic elements in marine bacteria is poorly documented. To gain insight into the genetic mechanisms underlying AMR gene transfer from marine bacteria, we mated a multidrug-resistant Vibrio alfacsensis strain with an Escherichia coli strain, and then determined the complete genome sequences of the donor and the transconjugant strains. Sequence analysis revealed a conjugative multidrug resistance plasmid in the donor strain, which was integrated into the chromosome of the recipient. The plasmid backbone in the transconjugant chromosome was flanked by two copies of a 7.1 kb unclassifiable integrative element harboring a β-lactamase gene. The 7.1 kb element and the previously reported element Tn6283 share four coding sequences, two of which encode the catalytic R-H-R-Y motif of tyrosine recombinases. Polymerase chain reaction and sequencing experiments revealed that these elements generate a circular copy of one specific strand without leaving an empty site on the donor molecule, in contrast to the movement of integron gene cassettes or ICE/IMEs discovered to date. These elements are termed SEs (strand-biased circularizing integrative elements): SE-6945 (the 7.1 kb element) and SE-6283 (Tn6283). The copy number and location of SE-6945 in the chromosome affected the antibiotic resistance levels of the transconjugants. SEs were identified in the genomes of other Vibrio species. Overall, these results suggest that SEs are involved in the spread of AMR genes among marine bacteria.

Structure, mechanism, and phylogeny of LysM-chitinase conjugates specifically found in fern plants

A unique GH18 chitinase containing two N-terminal lysin motifs (PrLysM1 and PrLysM2) was first found in fern, Pteris ryukyuensis (Onaga and Taira, Glycobiology, 18, 414-423, 2008). This type of LysM-chitinase conjugates is not usually found in plants but in fungi. Here, we produced a similar GH18 chitinase with one N-terminal LysM module (EaLysM) from the fern, Equisetum arvense (EaChiA, Inamine et al., Biosci. Biotechnol. Biochem., 79, 1296-1304, 2015), using an Escherichia coli expression system and characterized for its structure and mechanism of action. The crystal structure of EaLysM exhibited an almost identical fold (βααβ) to that of PrLysM2. From isothermal titration calorimetry and nuclear magnetic resonance, the binding mode and affinities of EaLysM for chitooligosaccharides (GlcNAc)_n (3, 4, 5, and 6) were found to be comparable to those of PrLysM2. The LysM module in EaChiA is likely to bind (GlcNAc)_n almost independently through CH-π stacking of a Tyr residue with the pyranose ring. The (GlcNAc)_n-binding mode of LysMs in the LysM-chitinase conjugates from fern plants appears to differ from that of plant LysMs acting in chitin- or Nod-signal perception, in which multiple LysMs cooperatively act on (GlcNAc)_n. Phylogenetic analysis suggested that LysM-GH18 conjugates of fern plants formed a monophyletic group and had been separated earlier than forming the clade of fungal chitinases with LysMs.

Redefining the architecture of ferlin proteins: Insights into multi-domain protein structure and function

Ferlins are complex, multi-domain proteins, involved in membrane trafficking, membrane repair, and exocytosis. The large size of ferlin proteins and the lack of consensus regarding domain boundaries have slowed progress in understanding molecular-level details of ferlin protein structure and function. However, in silico protein folding techniques have significantly enhanced our understanding of the complex ferlin family domain structure. We used RoseTTAFold to assemble full-length models for the six human ferlin proteins (dysferlin, myoferlin, otoferlin, Fer1L4, Fer1L5, and Fer1L6). Our full-length ferlin models were used to obtain objective domain boundaries, and these boundaries were supported by AlphaFold2 predictions. Despite the differences in amino acid sequence between the ferlin proteins, the domain ranges and distinct subdomains in the ferlin domains are remarkably consistent. Further, the RoseTTAFold/AlphaFold2 in silico boundary predictions allowed us to describe and characterize a previously unknown C2 domain, ubiquitous in all human ferlins, which we refer to as C2-FerA. At present, the ferlin domain-domain interactions implied by the full-length in silico models are predicted to have a low accuracy; however, the use of RoseTTAFold and AlphaFold2 as a domain finder has proven to be a powerful research tool for understanding ferlin structure.

In silico binding affinity prediction for metabotropic glutamate receptors using both endpoint free energy methods and a machine learning-based scoring function

Metabotropic glutamate receptors (mGluRs) play an important role in regulating glutamate signal pathways, which are involved in neuropathy and periphery homeostasis. mGluR4, which belongs to Group III mGluRs, is most widely distributed in the periphery among all the mGluRs. It has been proved that the regulation of this receptor is involved in diabetes, colorectal carcinoma and many other diseases. However, the application of structure-based drug design to identify small molecules to regulate the mGluR4 receptor is limited due to the absence of a resolved mGluR4 protein structure. In this work, we first built a homology model of mGluR4 based on a crystal structure of mGluR8, and then conducted hierarchical virtual screening (HVS) to identify possible active ligands for mGluR4. The HVS protocol consists of three hierarchical filters including Glide docking, molecular dynamic (MD) simulation and binding free energy calculation. We successfully prioritized active ligands of mGluR4 from a set of screening compounds using HVS. The predicted active ligands based on binding affinities can almost cover all the experiment-determined active ligands, with only one ligand missed. The correlation between the measured and predicted binding affinities is significantly improved for the MM-PB/GBSA-WSAS methods compared to the Glide docking method. More importantly, we have identified hotspots for ligand binding, and we found that SER157 and GLY158 tend to contribute to the selectivity of mGluR4 ligands, while ALA154 and ALA155 could account for the ligand selectivity to mGluR8. We also recognized other 5 key residues that are critical for ligand potency. The difference of the binding profiles between mGluR4 and mGluR8 can guide us to develop more potent and selective modulators. Moreover, we evaluated the performance of IPSF, a novel type of scoring function trained by a machine learning algorithm on residue-ligand interaction profiles, in guiding drug lead optimization. The cross-validation root-mean-square errors (RMSEs) are much smaller than those by the endpoint methods, and the correlation coefficients are comparable to the best endpoint methods for both mGluRs. Thus, machine learning-based IPSF can be applied to guide lead optimization, albeit the total number of actives/inactives are not big, a typical scenario in drug discovery projects.

Silybins are stereospecific regulators of the 20S proteasome

A reduced proteasome activity tiles excessive amyloid growth during the progress of protein conformational diseases (PCDs). Hence, the development of safe and effective proteasome enhancers represents an attractive target for the therapeutic treatment of these chronic disorders. Here we analyze two natural diastereoisomers belonging to the family of flavonolignans, Sil A and Sil B, by evaluating their capacity to increase proteasome activity. Enzyme assays carried out on yeast 20S (y20S) proteasome and in parallel on a permanently "open gate" mutant (α3ΔN) evidenced that Sil B is a more efficient 20S activator than Sil A. Conversely, in the case of human 20S proteasome (h20S) a higher affinity and more efficient activation is observed for Sil A. Driven by experimental data, computational studies further demonstrated that the taxifolin group of both diastereoisomers plays a crucial role in their anchoring to the α5/α6 groove of the outer α-ring. However, due to the different stereochemistry at C-7" and C-8" of ring D, only Sil A was able to reproduce the interactions responsible for h20S proteasome activation induced by their cognate regulatory particles. The provided silybins/h20S interaction models allowed us to rationalize their different ability to activate the peptidase activities of h20S and y20S. Our results provide structural details concerning the important role played by stereospecific interactions in driving Sil A and Sil B binding to the 20S proteasome and may support future rational design of proteasome enhancers.

Comparing the folding landscapes of evolutionarily divergent procaspase-3

All caspases evolved from a common ancestor and subsequently developed into two general classes, inflammatory or apoptotic caspases. The caspase-hemoglobinase fold has been conserved throughout nearly one billion years of evolution and is utilized for both the monomeric and dimeric subfamilies of apoptotic caspases, called initiator and effector caspases, respectively. We compared the folding and assembly of procaspase-3b from zebrafish to that of human effector procaspases in order to examine the conservation of the folding landscape. Urea-induced equilibrium folding/unfolding of procaspase-3b showed a minimum three-state folding pathway, where the native dimer isomerizes to a partially folded dimeric intermediate, which then unfolds. A partially folded monomeric intermediate observed in the folding landscape of human procaspase-3 is not well-populated in zebrafish procaspase-3b. By comparing effector caspases from different species, we show that the effector procaspase dimer undergoes a pH-dependent conformational change, and that the conformational species in the folding landscape exhibit similar free energies. Together, the data show that the landscape for the caspase-hemoglobinase fold is conserved, yet it provides flexibility for species-specific stabilization or destabilization of folding intermediates resulting in changes in stability. The common pH-dependent conformational change in the native dimer, which yields an enzymatically inactive species, may provide an additional, albeit reversible, mechanism for controlling caspase activity in the cell.

Modularity of the hydrophobic core and evolution of functional diversity in fold A glycosyltransferases

Hydrophobic cores are fundamental structural properties of proteins typically associated with protein folding and stability; however, how the hydrophobic core shapes protein evolution and function is poorly understood. Here, we investigated the role of conserved hydrophobic cores in fold-A glycosyltransferases (GT-As), a large superfamily of enzymes that catalyze formation of glycosidic linkages between diverse donor and acceptor substrates through distinct catalytic mechanisms (inverting versus retaining). Using hidden Markov models and protein structural alignments, we identify similarities in the phosphate-binding cassette (PBC) of GT-As and unrelated nucleotide-binding proteins, such as UDP-sugar pyrophosphorylases. We demonstrate that GT-As have diverged from other nucleotide-binding proteins through structural elaboration of the PBC and its unique hydrophobic tethering to the F-helix, which harbors the catalytic base (xED-Asp). While the hydrophobic tethering is conserved across diverse GT-A fold enzymes, some families, such as B3GNT2, display variations in tethering interactions and core packing. We evaluated the structural and functional impact of these core variations through experimental mutational analysis and molecular dynamics simulations and find that some of the core mutations (T336I in B3GNT2) increase catalytic efficiency by modulating the conformational occupancy of the catalytic base between “D-in” and acceptor-accessible “D-out” conformation. Taken together, our studies support a model of evolution in which the GT-A core evolved progressively through elaboration upon an ancient PBC found in diverse nucleotide-binding proteins, and malleability of this core provided the structural framework for evolving new catalytic and substrate-binding functions in extant GT-A fold enzymes.

Structures and mechanism of the plant PIN-FORMED auxin transporter

Auxins are hormones that have central roles and control nearly all aspects of growth and development in plants^1–3. The proteins in the PIN-FORMED (PIN) family (also known as the auxin efflux carrier family) are key participants in this process and control auxin export from the cytosol to the extracellular space^4–9. Owing to a lack of structural and biochemical data, the molecular mechanism of PIN-mediated auxin transport is not understood. Here we present biophysical analysis together with three structures of Arabidopsis thaliana PIN8: two outward-facing conformations with and without auxin, and one inward-facing conformation bound to the herbicide naphthylphthalamic acid. The structure forms a homodimer, with each monomer divided into a transport and scaffold domain with a clearly defined auxin binding site. Next to the binding site, a proline–proline crossover is a pivot point for structural changes associated with transport, which we show to be independent of proton and ion gradients and probably driven by the negative charge of the auxin. The structures and biochemical data reveal an elevator-type transport mechanism reminiscent of bile acid/sodium symporters, bicarbonate/sodium symporters and sodium/proton antiporters. Our results provide a comprehensive molecular model for auxin recognition and transport by PINs, link and expand on a well-known conceptual framework for transport, and explain a central mechanism of polar auxin transport, a core feature of plant physiology, growth and development.

Structural and biophysical analysis of the Arabidopsis thaliana auxin transporter PIN8 reveal that PIN transporters export auxin using an elevator mechanism.

A new paradigm of multiheme cytochrome evolution by grafting and pruning protein modules

Multiheme cytochromes play key roles in diverse biogeochemical cycles, but understanding the origin and evolution of these proteins is a challenge due to their ancient origin and complex structure. Up until now, the evolution of multiheme cytochromes composed by multiple redox modules in a single polypeptide chain was proposed to occur by gene fusion events. In this context, the pentaheme nitrite reductase NrfA and the tetraheme cytochrome c₅₅₄ were previously proposed to be at the origin of the extant octa- and nonaheme cytochrome c involved in metabolic pathways that contribute to the nitrogen, sulfur, and iron biogeochemical cycles by a gene fusion event. Here, we combine structural and character-based phylogenetic analysis with an unbiased root placement method to refine the evolutionary relationships between these multiheme cytochromes. The evidence show that NrfA and cytochrome c₅₅₄ belong to different clades, which suggests that these two multiheme cytochromes are products of truncation of ancestral octaheme cytochromes related to extant octaheme nitrite reductase and MccA, respectively. From our phylogenetic analysis, the last common ancestor is predicted to be an octaheme cytochrome with nitrite reduction ability. Evolution from this octaheme framework led to the great diversity of extant multiheme cytochromes analyzed here by pruning and grafting of protein modules and hemes. By shedding light into the evolution of multiheme cytochromes that intervene in different biogeochemical cycles, this work contributes to our understanding about the interplay between biology and geochemistry across large time scales in the history of Earth.

Structural variations between small alarmone hydrolase dimers support different modes of regulation of the stringent response

The bacterial stringent response involves wide-ranging metabolic reprogramming aimed at increasing long-term survivability during stress conditions. One of the hallmarks of the stringent response is the production of a set of modified nucleotides, known as alarmones, which affect a multitude of cellular pathways in diverse ways. Production and degradation of these molecules depend on the activity of enzymes from the RelA/SpoT homologous family, which come in both bifunctional (containing domains to both synthesize and hydrolyze alarmones) and monofunctional (consisting of only synthetase or hydrolase domain) variants, of which the structure, activity, and regulation of the bifunctional RelA/SpoT homologs have been studied most intensely. Despite playing an important role in guanosine nucleotide homeostasis in particular, mechanisms of regulation of the small alarmone hydrolases (SAHs) are still rather unclear. Here, we present crystal structures of SAH enzymes from Corynebacterium glutamicum (RelH_Cg) and Leptospira levettii (RelH_Ll) and show that while being highly similar, structural differences in substrate access and dimer conformations might be important for regulating their activity. We propose that a varied dimer form is a general property of the SAH family, based on current structural information as well as prediction models for this class of enzymes. Finally, subtle structural variations between monofunctional and bifunctional enzymes point to how these different classes of enzymes are regulated.

Alternative splicing encodes functional intracellular CD59 isoforms that mediate insulin secretion and are down-regulated in diabetic islets

This project describes the existence of previously unknown non–GPI-anchored CD59 isoforms required for insulin secretion, named CD59–IRIS-1 and CD59–IRIS-2, and finds reduced expression of CD59-IRIS isoforms in human diabetic islets, showing a link between dysregulation of IRIS isoforms and defects in insulin secretion in diabetic patients. These data open a path for future studies into CD59-IRIS expression and function in additional cell types capable of regulated secretion. Identification of additional specific CD59-IRIS binding partners within the cell could provide therapeutic targets for enhancement of insulin secretion in T2D.

Human pancreatic islets highly express CD59, which is a glycosylphosphatidylinositol (GPI)-anchored cell-surface protein and is required for insulin secretion. How cell-surface CD59 could interact with intracellular exocytotic machinery has so far not been described. We now demonstrate the existence of CD59 splice variants in human pancreatic islets, which have unique C-terminal domains replacing the GPI-anchoring signal sequence. These isoforms are found in the cytosol of β-cells, interact with SNARE proteins VAMP2 and SNAP25, colocalize with insulin granules, and rescue insulin secretion in CD59-knockout (KO) cells. We therefore named these isoforms IRIS-1 and IRIS-2 (Isoforms Rescuing Insulin Secretion 1 and 2). Antibodies raised against each isoform revealed that expression of both IRIS-1 and IRIS-2 is significantly lower in islets isolated from human type 2 diabetes (T2D) patients, as compared to healthy controls. Further, glucotoxicity induced in primary, healthy human islets led to a significant decrease of IRIS-1 expression, suggesting that hyperglycemia (raised glucose levels) and subsequent decreased IRIS-1 expression may contribute to relative insulin deficiency in T2D patients. Similar isoforms were also identified in the mouse CD59B gene, and targeted CRISPR/Cas9-mediated knockout showed that these intracellular isoforms, but not canonical CD59B, are involved in insulin secretion from mouse β-cells. Mouse IRIS-2 is also down-regulated in diabetic db/db mouse islets. These findings establish the endogenous existence of previously undescribed non–GPI-anchored intracellular isoforms of human CD59 and mouse CD59B, which are required for normal insulin secretion.

Architecture of the linker-scaffold in the nuclear pore

INTRODUCTION In eukaryotic cells, the selective bidirectional transport of macromolecules between the nucleus and cytoplasm occurs through the nuclear pore complex (NPC). Embedded in nuclear envelope pores, the ~110-MDa human NPC is an ~1200-Å-wide and ~750-Å-tall assembly of ~1000 proteins, collectively termed nucleoporins. Because of the NPC's eightfold rotational symmetry along the nucleocytoplasmic axis, each of the ~34 different nucleoporins occurs in multiples of eight. Architecturally, the NPC's symmetric core is composed of an inner ring encircling the central transport channel and two outer rings anchored on both sides of the nuclear envelope. Because of its central role in the flow of genetic information from DNA to RNA to protein, the NPC is commonly targeted in viral infections and its nucleoporin constituents are associated with a plethora of diseases. RATIONALE Although the arrangement of most scaffold nucleoporins in the NPC's symmetric core was determined by quantitative docking of crystal structures into cryo-electron tomographic (cryo-ET) maps of intact NPCs, the topology and molecular details of their cohesion by multivalent linker nucleoporins have remained elusive. Recently, in situ cryo-ET reconstructions of NPCs from various species have indicated that the NPC's inner ring is capable of reversible constriction and dilation in response to variations in nuclear envelope membrane tension, thereby modulating the diameter of the central transport channel by ~200 Å. We combined biochemical reconstitution, high-resolution crystal and single-particle cryo-electron microscopy (cryo-EM) structure determination, docking into cryo-ET maps, and physiological validation to elucidate the molecular architecture of the linker-scaffold interaction network that not only is essential for the NPC's integrity but also confers the plasticity and robustness necessary to allow and withstand such large-scale conformational changes. RESULTS By biochemically mapping scaffold-binding regions of all fungal and human linker nucleoporins and determining crystal and single-particle cryo-EM structures of linker-scaffold complexes, we completed the characterization of the biochemically tractable linker-scaffold network and established its evolutionary conservation, despite considerable sequence divergence. We determined a series of crystal and single-particle cryo-EM structures of the intact Nup188 and Nup192 scaffold hubs bound to their Nic96, Nup145N, and Nup53 linker nucleoporin binding regions, revealing that both proteins form distinct question mark-shaped keystones of two evolutionarily conserved hetero‑octameric inner ring complexes. Linkers bind to scaffold surface pockets through short defined motifs, with flanking regions commonly forming additional disperse interactions that reinforce the binding. Using a structure‑guided functional analysis in Saccharomyces cerevisiae, we confirmed the robustness of linker‑scaffold interactions and established the physiological relevance of our biochemical and structural findings. The near-atomic composite structures resulting from quantitative docking of experimental structures into human and S. cerevisiae cryo-ET maps of constricted and dilated NPCs structurally disambiguated the positioning of the Nup188 and Nup192 hubs in the intact fungal and human NPC and revealed the topology of the linker-scaffold network. The linker-scaffold gives rise to eight relatively rigid inner ring spokes that are flexibly interconnected to allow for the formation of lateral channels. Unexpectedly, we uncovered that linker‑scaffold interactions play an opposing role in the outer rings by forming tight cross-link staples between the eight nuclear and cytoplasmic outer ring spokes, thereby limiting the dilatory movements to the inner ring. CONCLUSION We have substantially advanced the structural and biochemical characterization of the symmetric core of the S. cerevisiae and human NPCs and determined near-atomic composite structures. The composite structures uncover the molecular mechanism by which the evolutionarily conserved linker‑scaffold establishes the NPC's integrity while simultaneously allowing for the observed plasticity of the central transport channel. The composite structures are roadmaps for the mechanistic dissection of NPC assembly and disassembly, the etiology of NPC‑associated diseases, the role of NPC dilation in nucleocytoplasmic transport of soluble and integral membrane protein cargos, and the anchoring of asymmetric nucleoporins. [Figure: see text].

Cloning the barley nec3 disease lesion mimic mutant using complementation by sequencing

Disease lesion mimic (DLM) or necrotic mutants display necrotic lesions in the absence of pathogen infections. They can show improved resistance to some pathogens and their molecular dissection can contribute to revealing components of plant defense pathways. Although forward-genetics strategies to find genes causal to mutant phenotypes are available in crops, these strategies require the production of experimental cross populations, mutagenesis, or gene editing and are time- and resource-consuming or may have to deal with regulated plant materials. In this study, we described a collection of 34 DLM mutants in barley (Hordeum vulgare L.) and applied a novel method called complementation by sequencing (CBS), which enables the identification of the gene responsible for a mutant phenotype given the availability of two or more chemically mutagenized individuals showing the same phenotype. Complementation by sequencing relies on the feasibility to obtain all induced mutations present in chemical mutants and on the low probability that different individuals share the same mutated genes. By CBS, we identified a cytochrome P450 CYP71P1 gene as responsible for orange blotch DLM mutants, including the historical barley nec3 locus. By comparative phylogenetic analysis we showed that CYP71P1 gene family emerged early in angiosperm evolution but has been recurrently lost in some lineages including Arabidopsis thaliana (L.) Heynh. Complementation by sequencing is a straightforward cost-effective approach to clone genes controlling phenotypes in a chemically mutagenized collection. The TILLMore (TM) collection will be instrumental for understanding the molecular basis of DLM phenotypes and to contribute knowledge about mechanisms of host-pathogen interaction.

Structural evolution of the ancient enzyme, dissimilatory sulfite reductase

Dissimilatory sulfite reductase is an ancient enzyme that has linked the global sulfur and carbon biogeochemical cycles since at least 3.47 Gya. While much has been learned about the phylogenetic distribution and diversity of DsrAB across environmental gradients, far less is known about the structural changes that occurred to maintain DsrAB function as the enzyme accompanied diversification of sulfate/sulfite reducing organisms (SRO) into new environments. Analyses of available crystal structures of DsrAB from Archaeoglobus fulgidus and Desulfovibrio vulgaris, representing early and late evolving lineages, respectively, show that certain features of DsrAB are structurally conserved, including active siro-heme binding motifs. Whether such structural features are conserved among DsrAB recovered from varied environments, including hot spring environments that host representatives of the earliest evolving SRO lineage (e.g., MV2-Eury), is not known. To begin to overcome these gaps in our understanding of the evolution of DsrAB, structural models from MV2.Eury were generated and evolutionary sequence co-variance analyses were conducted on a curated DsrAB database. Phylogenetically diverse DsrAB harbor many conserved functional residues including those that ligate active siro-heme(s). However, evolutionary co-variance analysis of monomeric DsrAB subunits revealed several False Positive Evolutionary Couplings (FPEC) that correspond to residues that have co-evolved despite being too spatially distant in the monomeric structure to allow for direct contact. One set of FPECs corresponds to residues that form a structural path between the two active siro-heme moieties across the interface between heterodimers, suggesting the potential for allostery or electron transfer within the enzyme complex. Other FPECs correspond to structural loops and gaps that may have been selected to stabilize enzyme function in different environments. These structural bioinformatics results suggest that DsrAB has maintained allosteric communication pathways between subunits as SRO diversified into new environments. The observations outlined here provide a framework for future biochemical and structural analyses of DsrAB to examine potential allosteric control of this enzyme.

Type II BMP and activin receptors BMPR2 and ACVR2A share a conserved mode of growth factor recognition

BMPR2 is a type II Transforming Growth Factor (TGF)-β family receptor that is fundamentally associated with pulmonary arterial hypertension (PAH) in humans. BMPR2 shares functional similarities with the type II activin receptors ACVR2A and ACVR2B, as it interacts with an overlapping group of TGF-β family growth factors (GFs). However, how BMPR2 recognizes GFs remains poorly understood. Here, we solved crystal structures of BMPR2 in complex with the GF activin B and of ACVR2A in complex with the related GF activin A. We show that both BMPR2 and ACVR2A bind GFs with nearly identical geometry using a conserved hydrophobic hot spot, while differences in contacting residues are predominantly found in loop areas. Upon further exploration of the GF-binding spectrum of the two receptors, we found that although many GFs bind both receptors, the high-affinity BMPR2 GFs comprise BMP15, BMP10, and Nodal, whereas those of ACVR2A are activin A, activin B, and GDF11. Lastly, we evaluated GF-binding domain BMPR2 variants found in human PAH patients. We demonstrate that mutations within the GF-binding interface resulted in loss of GF binding, while mutations in loop areas allowed BMPR2 to retain the ability to bind cognate GFs with high affinity. In conclusion, the in vitro activities of BMPR2 variants and the crystal structures reported here indicate biochemically relevant complexes that explain how some GF-binding domain variants can lead to PAH.

Structural basis of transposon end recognition explains central features of Tn7 transposition systems

Tn7 is a bacterial transposon with relatives containing element-encoded CRISPR-Cas systems mediating RNA-guided transposon insertion. Here, we present the 2.7 Å cryoelectron microscopy structure of prototypic Tn7 transposase TnsB interacting with the transposon end DNA. When TnsB interacts across repeating binding sites, it adopts a beads-on-a-string architecture, where the DNA-binding and catalytic domains are arranged in a tiled and intertwined fashion. The DNA-binding domains form few base-specific contacts leading to a binding preference that requires multiple weakly conserved sites at the appropriate spacing to achieve DNA sequence specificity. TnsB binding imparts differences in the global structure of the protein-bound DNA ends dictated by the spacing or overlap of binding sites explaining functional differences in the left and right ends of the element. We propose a model of the strand-transfer complex in which the terminal TnsB molecule is rearranged so that its catalytic domain is in a position conducive to transposition.

The W-Acidic Motif of Histidine Kinase WalK Is Required for Signaling and Transcriptional Regulation in Streptococcus mutans

In Streptococcus mutans, we find that the histidine kinase WalK possesses the longest C-terminal tail (CTT) among all 14 TCSs, and this tail plays a key role in the interaction of WalK with its response regulator WalR. We demonstrate that the intrinsically disordered CTT is characterized by a conserved tryptophan residue surrounded by acidic amino acids. Mutation in the tryptophan not only disrupts the stable interaction, but also impairs the efficient phosphotransferase and phosphatase activities of WalRK. In addition, the tryptophan is important for WalK to compete with DNA containing a WalR binding motif for the WalR interaction. We further show that the tryptophan is important for in vivo transcriptional regulation and bacterial biofilm formation by S. mutans. Moreover, Staphylococcus aureus WalK also has a characteristic CTT, albeit relatively shorter, with a conserved W-acidic motif, that is required for the WalRK interaction in vitro. Together, these data reveal that the W-acidic motif of WalK is indispensable for its interaction with WalR, thereby playing a key role in the WalRK-dependent signal transduction, transcriptional regulation and biofilm formation.

Molecular Mingling: Multimodal Predictions of Ligand Promiscuity in Pentameric Ligand-Gated Ion Channels

Background: Human pentameric ligand-gated ion channels (pLGICs) comprise nicotinic acetylcholine receptors (nAChRs), 5-hydroxytryptamine type 3 receptors (5-HT3Rs), zinc-activated channels (ZAC), γ-aminobutyric acid type A receptors (GABAARs) and glycine receptors (GlyRs). They are recognized therapeutic targets of some of the most prescribed drugs like general anesthetics, anxiolytics, smoking cessation aids, antiemetics and many more. Currently, approximately 100 experimental structures of pLGICs with ligands bound exist in the protein data bank (PDB). These atomic-level 3D structures enable the generation of a comprehensive binding site inventory for the superfamily and the in silico prediction of binding site properties. Methods: A panel of high throughput in silico methods including pharmacophore screening, conformation analysis and descriptor calculation was applied to a selection of allosteric binding sites for which in vitro screens are lacking. Variant abundance near binding site forming regions and computational docking complement the approach. Results: The structural data reflects known and novel binding sites, some of which may be unique to individual receptors, while others are broadly conserved. The membrane spanning domain, comprising four highly conserved segments, contains ligand interaction sites for which in vitro assays suitable for high throughput screenings are critically lacking. This is also the case for structurally more variable novel sites in the extracellular domain. Our computational results suggest that the phytocannabinoid Δ9-tetrahydrocannabinol (Δ9-THC) can utilize multiple pockets which are likely to exist on most superfamily members. Conclusion: With this study, we explore the potential for polypharmacology among pLGICs. Our data suggest that ligands can display two forms of promiscuity to an extent greater than what has been realized: 1) Ligands can interact with homologous sites in many members of the superfamily, which bears toxicological relevance. 2) Multiple pockets in distinct localizations of individual receptor subtypes share common ligands, which counteracts efforts to develop selective agents. Moreover, conformational states need to be considered for in silico drug screening, as certain binding sites display considerable flexibility. In total, this work contributes to a better understanding of polypharmacology across pLGICs and provides a basis for improved structure guided in silico drug development and drug derisking.

Identification and Functional Characterization of Peptides With Antimicrobial Activity From the Syphilis Spirochete, Treponema pallidum

The etiological agent of syphilis, Treponema pallidum ssp. pallidum, is a highly invasive “stealth” pathogen that can evade the host immune response and persist within the host for decades. This obligate human pathogen is adept at establishing infection and surviving at sites within the host that have a multitude of competing microbes, sometimes including pathogens. One survival strategy employed by bacteria found at polymicrobial sites is elimination of competing microorganisms by production of antimicrobial peptides (AMPs). Antimicrobial peptides are low molecular weight proteins (miniproteins) that function directly via inhibition and killing of microbes and/or indirectly via modulation of the host immune response, which can facilitate immune evasion. In the current study, we used bioinformatics to show that approximately 7% of the T. pallidum proteome is comprised of miniproteins of 150 amino acids or less with unknown functions. To investigate the possibility that AMP production is an unrecognized defense strategy used by T. pallidum during infection, we developed a bioinformatics pipeline to analyze the complement of T. pallidum miniproteins of unknown function for the identification of potential AMPs. This analysis identified 45 T. pallidum AMP candidates; of these, Tp0451a and Tp0749 were subjected to further bioinformatic analyses to identify AMP critical core regions (AMPCCRs). Four potential AMPCCRs from the two predicted AMPs were identified and peptides corresponding to these AMPCCRs were experimentally confirmed to exhibit bacteriostatic and bactericidal activity against a panel of biologically relevant Gram-positive and Gram-negative bacteria. Immunomodulation assays performed under inflammatory conditions demonstrated that one of the AMPCCRs was also capable of differentially regulating expression of two pro-inflammatory chemokines [monocyte chemoattractant protein-1 (MCP-1) and interleukin-8 (IL-8)]. These findings demonstrate proof-of-concept for our developed AMP identification pipeline and are consistent with the novel concept that T. pallidum expresses AMPs to defend against competing microbes and modulate the host immune response.

A bipartite chromatophore transit peptide and N-terminal protein processing in the Paulinella chromatophore

The amoeba Paulinella chromatophora contains photosynthetic organelles, termed chromatophores, which evolved independently from plastids in plants and algae. At least one-third of the chromatophore proteome consists of nucleus-encoded (NE) proteins that are imported across the chromatophore double envelope membranes. Chromatophore-targeted proteins exceeding 250 amino acids (aa) carry a conserved N-terminal extension presumably involved in protein targeting, termed the chromatophore transit peptide (crTP). Short imported proteins do not carry discernable targeting signals. To explore whether the import of proteins is accompanied by their N-terminal processing, here we identified N-termini of 208 chromatophore-localized proteins by a mass spectrometry-based approach. Our study revealed extensive N-terminal acetylation and proteolytic processing in both NE and chromatophore-encoded (CE) fractions of the chromatophore proteome. Mature N-termini of 37 crTP-carrying proteins were identified, of which 30 were cleaved in a common processing region. Surprisingly, only the N-terminal ∼50 aa (part 1) become cleaved upon import. This part contains a conserved adaptor protein-1 complex-binding motif known to mediate protein sorting at the trans-Golgi network followed by a predicted transmembrane helix, implying that part 1 anchors the protein co-translationally in the endoplasmic reticulum and mediates trafficking to the chromatophore via the Golgi. The C-terminal part 2 contains conserved secondary structural elements, remains attached to the mature proteins, and might mediate translocation across the chromatophore inner membrane. Short imported proteins remain largely unprocessed. Finally, this work illuminates N-terminal processing of proteins encoded in an evolutionary-early-stage organelle and suggests host-derived posttranslationally acting factors involved in regulation of the CE chromatophore proteome.

Domain Expansion and Functional Diversification in Vertebrate Reproductive Proteins

The rapid evolution of fertilization proteins has generated remarkable diversity in molecular structure and function. Glycoproteins of vertebrate egg coats contain multiple zona pellucida (ZP)-N domains (1-6 copies) that facilitate multiple reproductive functions, including species-specific sperm recognition. In this report, we integrate phylogenetics and machine learning to investigate how ZP-N domains diversify in structure and function. The most C-terminal ZP-N domain of each paralog is associated with another domain type (ZP-C), which together form a "ZP module." All modular ZP-N domains are phylogenetically distinct from nonmodular or free ZP-N domains. Machine learning-based classification identifies eight residues that form a stabilizing network in modular ZP-N domains that is absent in free domains. Positive selection is identified in some free ZP-N domains. Our findings support that strong purifying selection has conserved an essential structural core in modular ZP-N domains, with the relaxation of this structural constraint allowing free N-terminal domains to functionally diversify.

Comparison of PsbQ and Psb27 in photosystem II provides insight into their roles

Photosystem II (PSII) catalyzes the oxidation of water at its active site that harbors a high-valent inorganic Mn₄CaO_x cluster called the oxygen-evolving complex (OEC). Extrinsic subunits generally serve to protect the OEC from reductants and stabilize the structure, but diversity in the extrinsic subunits exists between phototrophs. Recent cryo-electron microscopy experiments have provided new molecular structures of PSII with varied extrinsic subunits. We focus on the extrinsic subunit PsbQ, that binds to the mature PSII complex, and on Psb27, an extrinsic subunit involved in PSII biogenesis. PsbQ and Psb27 share a similar binding site and have a four-helix bundle tertiary structure, suggesting they are related. Here, we use sequence alignments, structural analyses, and binding simulations to compare PsbQ and Psb27 from different organisms. We find no evidence that PsbQ and Psb27 are related despite their similar structures and binding sites. Evolutionary divergence within PsbQ homologs from different lineages is high, probably due to their interactions with other extrinsic subunits that themselves exhibit vast diversity between lineages. This may result in functional variation as exemplified by large differences in their calculated binding energies. Psb27 homologs generally exhibit less divergence, which may be due to stronger evolutionary selection for certain residues that maintain its function during PSII biogenesis and this is consistent with their more similar calculated binding energies between organisms. Previous experimental inconsistencies, low confidence binding simulations, and recent structural data suggest that Psb27 is likely to exhibit flexibility that may be an important characteristic of its activity. The analysis provides insight into the functions and evolution of PsbQ and Psb27, and an unusual example of proteins with similar tertiary structures and binding sites that probably serve different roles.

Molecular Organisation of Tick-Borne Encephalitis Virus

Tick-borne encephalitis virus (TBEV) is a pathogenic, enveloped, positive-stranded RNA virus in the family Flaviviridae. Structural studies of flavivirus virions have primarily focused on mosquito-borne species, with only one cryo-electron microscopy (cryo-EM) structure of a tick-borne species published. Here, we present a 3.3 Å cryo-EM structure of the TBEV virion of the Kuutsalo-14 isolate, confirming the overall organisation of the virus. We observe conformational switching of the peripheral and transmembrane helices of M protein, which can explain the quasi-equivalent packing of the viral proteins and highlights their importance in stabilising membrane protein arrangement in the virion. The residues responsible for M protein interactions are highly conserved in TBEV but not in the structurally studied Hypr strain, nor in mosquito-borne flaviviruses. These interactions may compensate for the lower number of hydrogen bonds between E proteins in TBEV compared to the mosquito-borne flaviviruses. The structure reveals two lipids bound in the E protein which are important for virus assembly. The lipid pockets are comparable to those recently described in mosquito-borne Zika, Spondweni, Dengue, and Usutu viruses. Our results thus advance the understanding of tick-borne flavivirus architecture and virion-stabilising interactions.

Recurrent Duplication and Diversification of Acrosomal Fertilization Proteins in Abalone

Reproductive proteins mediating fertilization commonly exhibit rapid sequence diversification driven by positive selection. This pattern has been observed among nearly all taxonomic groups, including mammals, invertebrates, and plants, and is remarkable given the essential nature of the molecular interactions mediating fertilization. Gene duplication is another important mechanism that facilitates the generation of molecular novelty through functional divergence. Following duplication, paralogs may partition ancestral gene function (subfunctionalization) or acquire new roles (neofunctionalization). However, the contributions of duplication followed by sequence diversification to the molecular diversity of gamete recognition genes has been understudied in many models of fertilization. The marine gastropod mollusk abalone is a classic model for fertilization. Its two acrosomal proteins (lysin and sp18) are ancient gene duplicates with unique gamete recognition functions. Through detailed genomic and bioinformatic analyses we show how duplication events followed by sequence diversification has played an ongoing role in the evolution of abalone acrosomal proteins. The common ancestor of abalone had four members of its acrosomal protein family in a tandem gene array that repeatedly experienced positive selection. We find that both sp18 paralogs contain positively selected sites located in different regions of the paralogs, suggestive of functional divergence where selection acted upon distinct binding interfaces in each paralog. Further, a more recent species-specific duplication of both lysin and sp18 in the European abalone H. tuberculata is described. Despite clade-specific acrosomal protein paralogs, there are no concomitant duplications of egg coat proteins in H. tuberculata, indicating that duplication of egg proteins per se is not responsible for retention of duplicated acrosomal proteins. We hypothesize that, in a manner analogous to host/pathogen evolution, sperm proteins are selected for increased diversity through extensive sequence divergence and recurrent duplication driven by conflict mechanisms.

The biosynthetic origin of ribofuranose in bacterial polysaccharides

Bacterial surface polysaccharides are assembled by glycosyltransferase enzymes that typically use sugar nucleotide or polyprenyl-monophosphosugar activated donors. Characterized representatives exist for many monosaccharides but neither the donor nor the corresponding glycosyltransferases have been definitively identified for ribofuranose residues found in some polysaccharides. Klebsiella pneumoniae O-antigen polysaccharides provided prototypes to identify dual-domain ribofuranosyltransferase proteins catalyzing a two-step reaction sequence. Phosphoribosyl-5-phospho-D-ribosyl-α-1-diphosphate serves as the donor for a glycan acceptor-specific phosphoribosyl transferase (gPRT), and a more promiscuous phosphoribosyl-phosphatase (PRP) then removes the residual 5'-phosphate. The 2.5-Å resolution crystal structure of a dual-domain ribofuranosyltransferase ortholog from Thermobacillus composti revealed a PRP domain that conserves many features of the phosphatase members of the haloacid dehalogenase family, and a gPRT domain that diverges substantially from all previously characterized phosphoribosyl transferases. The gPRT represents a new glycosyltransferase fold conserved in the most abundant ribofuranosyltransferase family.