High allelic diversity in Arabidopsis NLRs is associated with distinct genomic features.

Plants rely on Nucleotide-binding, Leucine-rich repeat Receptors (NLRs) for pathogen recognition. Highly variable NLRs (hvNLRs) show remarkable intraspecies diversity, while their low-variability paralogs (non-hvNLRs) are conserved between ecotypes. At a population level, hvNLRs provide new pathogen-recognition specificities, but the association between allelic diversity and genomic and epigenomic features has not been established. Our investigation of NLRs in Arabidopsis Col-0 has revealed that hvNLRs show higher expression, less gene body cytosine methylation, and closer proximity to transposable elements than non-hvNLRs. hvNLRs show elevated synonymous and nonsynonymous nucleotide diversity and are in chromatin states associated with an increased probability of mutation. Diversifying selection maintains variability at a subset of codons of hvNLRs, while purifying selection maintains conservation at non-hvNLRs. How these features are established and maintained, and whether they contribute to the observed diversity of hvNLRs is key to understanding the evolution of plant innate immune receptors.

The DUF3715 domain has a conserved role in RNA-directed transposon silencing.

RNA-directed transposon silencing operates in the mammalian soma and germline to safeguard genomic integrity. The piRNA pathway and the HUSH complex identify active transposons through recognition of their nascent transcripts, but mechanistic understanding of how these distinct pathways evolved is lacking. TASOR is an essential component of the HUSH complex. TASOR's DUF3715 domain adopts a pseudo-PARP structure and is required for transposon silencing in a manner independent of complex assembly. TEX15, an essential piRNA pathway factor, also contains the DUF3715 domain. Here, we show that TASOR's and TEX15's DUF3715 domain share extensive structural homology. We found that the DUF3715 domain arose in early eukaryotes and that in vertebrates it is restricted to TEX15, TASOR, and TASORB orthologs. While TASOR-like proteins are found throughout metazoa, TEX15 is vertebrate-specific. The branching of TEX15 and the TASOR-like DUF3715 domain likely occurred in early metazoan evolution. Remarkably, despite this vast evolutionary distance, the DUF3715 domain from divergent TEX15 sequences can functionally substitute the DUF3715 domain of TASOR and mediates transposon silencing. We have thus termed this domain of unknown function as the RNA-directed pseudo-PARP transposon silencing (RDTS) domain. In summary, we show an unexpected functional link between these critical transposon silencing pathways.

Majority of the highly variable NLRs in maize share genomic location and contain additional target-binding domains.

Nucleotide-binding, Leucine Rich Repeat proteins (NLRs) are a major class of immune receptors in plants. NLRs include both conserved and rapidly evolving members, however their evolutionary trajectory in crops remains understudied. Availability of crop pan-genomes enables analysis of the recent events in the evolution of this highly complex gene family within domesticated species. Here, we investigated the NLR complement of 26 nested association mapping (NAM) founder lines of maize. We found that maize has just four main subfamilies containing rapidly evolving highly variable NLR (hvNLR) receptors. Curiously, three of these phylogenetically distinct hvNLR lineages are located in adjacent clusters on chromosome 10. Members of the same hvNLR clade show variable expression and methylation across lines and tissues, consistent with their rapid evolution. By combining sequence diversity analysis and AlphaFold2 computational structure prediction we predicted ligand binding sites in the hvNLRs. We also observed novel insertion domains in the LRR regions of two hvNLR subfamilies that likely contribute to target recogniton. To make this analysis accessible, we created NLRCladeFinder, a Google Colaboratory notebook, that accepts any newly identified NLR sequence, places it in the evolutionary context of the maize pan-NLRome, and provides an updated clade alignment, phylogenetic tree, and sequence diversity information for the gene of interest.

Analysis of intraspecies diversity reveals a subset of highly variable plant immune receptors and predicts their binding sites.

The evolution of recognition specificities by the immune system depends on the generation of receptor diversity and on connecting the binding of new antigens with the initiation of downstream signaling. In plant immunity, the innate Nucleotide-Binding Leucine-Rich Repeat (NLR) receptor family enables antigen binding and immune signaling. In this study, we surveyed the NLR complements of 62 ecotypes of Arabidopsis thaliana and 54 lines of Brachypodium distachyon and identified a limited number of NLR subfamilies that show high allelic diversity. We show that the predicted specificity-determining residues cluster on the surfaces of Leucine-Rich Repeat domains, but the locations of the clusters vary among NLR subfamilies. By comparing NLR phylogeny, allelic diversity, and known functions of the Arabidopsis NLRs, we formulate a hypothesis for the emergence of direct and indirect pathogen-sensing receptors and of the autoimmune NLRs. These findings reveal the recurring patterns of evolution of innate immunity and can inform NLR engineering efforts.

Periphilin self-association underpins epigenetic silencing by the HUSH complex.

Transcription of integrated DNA from viruses or transposable elements is tightly regulated to prevent pathogenesis. The Human Silencing Hub (HUSH), composed of Periphilin, TASOR and MPP8, silences transcriptionally active viral and endogenous transgenes. HUSH recruits effectors that alter the epigenetic landscape and chromatin structure, but how HUSH recognizes target loci and represses their expression remains unclear. We identify the physicochemical properties of Periphilin necessary for HUSH assembly and silencing. A disordered N-terminal domain (NTD) and structured C-terminal domain are essential for silencing. A crystal structure of the Periphilin-TASOR minimal core complex shows Periphilin forms an α-helical homodimer, bound by a single TASOR molecule. The NTD forms insoluble aggregates through an arginine/tyrosine-rich sequence reminiscent of low-complexity regions from self-associating RNA-binding proteins. Residues required for TASOR binding and aggregation were required for HUSH-dependent silencing and genome-wide deposition of repressive mark H3K9me3. The NTD was functionally complemented by low-complexity regions from certain RNA-binding proteins and proteins that form condensates or fibrils. Our work suggests the associative properties of Periphilin promote HUSH aggregation at target loci.

Crystal Structure and Mechanistic Molecular Modeling Studies of Mycobacterium tuberculosis Diterpene Cyclase Rv3377c.

Terpenes make up the largest class of natural products, with extensive chemical and structural diversity. Diterpenes, mostly isolated from plants and rarely prokaryotes, exhibit a variety of important biological activities and valuable applications, including providing antitumor and antibiotic pharmaceuticals. These natural products are constructed by terpene synthases, a class of enzymes that catalyze one of the most complex chemical reactions in biology: converting simple acyclic oligo-isoprenyl diphosphate substrates to complex polycyclic products via carbocation intermediates. Here we obtained the second ever crystal structure of a class II diterpene synthase from bacteria, tuberculosinol pyrophosphate synthase (i.e., Halimadienyl diphosphate synthase, MtHPS, or Rv3377c) from Mycobacterium tuberculosis (Mtb). This enzyme transforms (E,E,E)-geranylgeranyl diphosphate into tuberculosinol pyrophosphate (Halimadienyl diphosphate). Rv3377c is part of the Mtb diterpene pathway along with Rv3378c, which converts tuberculosinol pyrophosphate to 1-tuberculosinyl adenosine (1-TbAd). This pathway was shown to exist only in virulent Mycobacterium species, but not in closely related avirulent species, and was proposed to be involved in phagolysosome maturation arrest. To gain further insight into the reaction pathway and the mechanistically relevant enzyme substrate binding orientation, electronic structure calculation and docking studies of reaction intermediates were carried out. Results reveal a plausible binding mode of the substrate that can provide the information to guide future drug design and anti-infective therapies of this biosynthetic pathway.

Biophysical Characterization of a Disabled Double Mutant of Soybean Lipoxygenase: The "Undoing" of Precise Substrate Positioning Relative to Metal Cofactor and an Identified Dynamical Network.

Soybean lipoxygenase (SLO) has served as a prototype for understanding the molecular origin of enzymatic rate accelerations. The double mutant (DM) L546A/L754A is considered a dramatic outlier, due to the unprecedented size and near temperature-independence of its primary kinetic isotope effect, low catalytic efficiency, and elevated enthalpy of activation. To uncover the physical basis of these features, we herein apply three structural probes: hydrogen-deuterium exchange mass spectrometry, room-temperature X-ray crystallography and EPR spectroscopy on four SLO variants (wild-type (WT) enzyme, DM, and the two parental single mutants, L546A and L754A). DM is found to incorporate features of each parent, with the perturbation at position 546 predominantly influencing thermally activated motions that connect the active site to a protein-solvent interface, while mutation at position 754 disrupts the ligand field and solvation near the cofactor iron. However, the expanded active site in DM leads to more active site water molecules and their associated hydrogen bond network, and the individual features from L546A and L754A alone cannot explain the aggregate kinetic properties for DM. Using recently published QM/MM-derived ground-state SLO-substrate complexes for WT and DM, together with the thorough structural analyses presented herein, we propose that the impairment of DM is the combined result of a repositioning of the reactive carbon of linoleic acid substrate with regard to both the iron cofactor and a catalytically linked dynamic region of protein.

Structural and Genetic Analyses of the Mycobacterium tuberculosis Protein Kinase B Sensor Domain Identify a Potential Ligand-binding Site.

Monitoring the environment with serine/threonine protein kinases is critical for growth and survival of Mycobacterium tuberculosis, a devastating human pathogen. Protein kinase B (PknB) is a transmembrane serine/threonine protein kinase that acts as an essential regulator of mycobacterial growth and division. The PknB extracellular domain (ECD) consists of four repeats homologous to penicillin-binding protein and serine/threonine kinase associated (PASTA) domains, and binds fragments of peptidoglycan. These properties suggest that PknB activity is modulated by ECD binding to peptidoglycan substructures, however, the molecular mechanisms underpinning PknB regulation remain unclear. In this study, we report structural and genetic characterization of the PknB ECD. We determined the crystal structures of overlapping ECD fragments at near atomic resolution, built a model of the full ECD, and discovered a region on the C-terminal PASTA domain that has the properties of a ligand-binding site. Hydrophobic interaction between this surface and a bound molecule of citrate was observed in a crystal structure. Our genetic analyses in M. tuberculosis showed that nonfunctional alleles were produced either by deletion of any of single PASTA domain or by mutation of individual conserved residues lining the putative ligand-binding surface of the C-terminal PASTA repeat. These results define two distinct structural features necessary for PknB signal transduction, a fully extended ECD and a conserved, membrane-distal putative ligand-binding site.

Mycobacterium tuberculosis FtsX extracellular domain activates the peptidoglycan hydrolase, RipC.

Bacterial growth and cell division are coordinated with hydrolysis of the peptidoglycan (PG) layer of the cell wall, but the mechanisms of regulation of extracellular PG hydrolases are not well understood. Here we report the biochemical, structural, and genetic analysis of the Mycobacterium tuberculosis homolog of the transmembrane PG-hydrolase regulator, FtsX. The purified FtsX extracellular domain binds the PG peptidase Rv2190c/RipC N-terminal segment, causing a conformational change that activates the enzyme. Deletion of ftsEX and ripC caused similar phenotypes in Mycobacterium smegmatis, as expected for genes in a single pathway. The crystal structure of the FtsX extracellular domain reveals an unprecedented fold containing two lobes connected by a flexible hinge. Mutations in the hydrophobic cleft between the lobes reduce RipC binding in vitro and inhibit FtsX function in M. smegmatis. These studies suggest how FtsX recognizes RipC and support a model in which a conformational change in FtsX links the cell division apparatus with PG hydrolysis.

Structural and biochemical analyses of Mycobacterium tuberculosis N-acetylmuramyl-L-alanine amidase Rv3717 point to a role in peptidoglycan fragment recycling.

Peptidoglycan hydrolases are key enzymes in bacterial cell wall homeostasis. Understanding the substrate specificity and biochemical activity of peptidoglycan hydrolases in Mycobacterium tuberculosis is of special interest as it can aid in the development of new cell wall targeting therapeutics. In this study, we report biochemical and structural characterization of the mycobacterial N-acetylmuramyl-L-alanine amidase, Rv3717. The crystal structure of Rv3717 in complex with a dipeptide product shows that, compared with previously characterized peptidoglycan amidases, the enzyme contains an extra disulfide-bonded ß-hairpin adjacent to the active site. The structure of two intermediates in assembly reveal that Zn(2+) binding rearranges active site residues, and disulfide formation promotes folding of the ß-hairpin. Although Zn(2+) is required for hydrolysis of muramyl dipeptide, disulfide oxidation is not required for activity on this substrate. The orientation of the product in the active site suggests a role for a conserved glutamate (Glu-200) in catalysis; mutation of this residue abolishes activity. The product binds at the head of a closed tunnel, and the enzyme showed no activity on polymerized peptidoglycan. These results point to a potential role for Rv3717 in peptidoglycan fragment recycling.

High intraspecies allelic diversity in Arabidopsis NLR immune receptors is associated with distinct genomic and epigenomic features.

Plants rely on Nucleotide-binding, Leucine-rich repeat Receptors (NLRs) for pathogen recognition. Highly variable NLRs (hvNLRs) show remarkable intraspecies diversity, while their low variability paralogs (non-hvNLRs) are conserved between ecotypes. At a population level, hvNLRs provide new pathogen recognition specificities, but the association between allelic diversity and genomic and epigenomic features has not been established. Our investigation of NLRs in Arabidopsis Col-0 has revealed that hvNLRs show higher expression, less gene body cytosine methylation, and closer proximity to transposable elements than non-hvNLRs. hvNLRs show elevated synonymous and nonsynonymous nucleotide diversity and are in chromatin states associated with an increased probability of mutation. Diversifying selection maintains variability at a subset of codons of hvNLRs, while purifying selection maintains conservation at non-hvNLRs. How these features are established and maintained, and whether they contribute to the observed diversity of hvNLRs is key to understanding the evolution of plant innate immune receptors.

A Zika virus mutation enhances transmission potential and confers escape from protective dengue virus immunity.

Zika virus (ZIKV) and dengue virus (DENV) are arthropod-borne pathogenic flaviviruses that co-circulate in many countries. To understand some of the pressures that influence ZIKV evolution, we mimic the natural transmission cycle by repeating serial passaging of ZIKV through cultured mosquito cells and either DENV-naive or DENV-immune mice. Compared with wild-type ZIKV, the strains passaged under both conditions exhibit increased pathogenesis in DENV-immune mice. Application of reverse genetics identifies an isoleucine-to-valine mutation (I39V) in the NS2B proteins of both passaged strains that confers enhanced fitness and escape from pre-existing DENV immunity. Introduction of I39V or I39T, a naturally occurring homologous mutation detected in recent ZIKV isolates, increases the replication of wild-type ZIKV in human neuronal precursor cells and laboratory-raised mosquitoes. Our data indicate that ZIKV strains with enhanced transmissibility and pathogenicity can emerge in DENV-naive or -immune settings, and that NS2B-I39 mutants may represent ZIKV variants of interest.

TASOR is a pseudo-PARP that directs HUSH complex assembly and epigenetic transposon control.

The HUSH complex represses retroviruses, transposons and genes to maintain the integrity of vertebrate genomes. HUSH regulates deposition of the epigenetic mark H3K9me3, but how its three core subunits - TASOR, MPP8 and Periphilin - contribute to assembly and targeting of the complex remains unknown. Here, we define the biochemical basis of HUSH assembly and find that its modular architecture resembles the yeast RNA-induced transcriptional silencing complex. TASOR, the central HUSH subunit, associates with RNA processing components. TASOR is required for H3K9me3 deposition over LINE-1 repeats and repetitive exons in transcribed genes. In the context of previous studies, this suggests that an RNA intermediate is important for HUSH activity. We dissect the TASOR and MPP8 domains necessary for transgene repression. Structure-function analyses reveal TASOR bears a catalytically-inactive PARP domain necessary for targeted H3K9me3 deposition. We conclude that TASOR is a multifunctional pseudo-PARP that directs HUSH assembly and epigenetic regulation of repetitive genomic targets.

A mouse-passaged dengue virus strain with reduced affinity for heparan sulfate causes severe disease in mice by establishing increased systemic viral loads.

The four serotypes of dengue virus (DENV1 to DENV4) cause extensive morbidity and mortality. A major obstacle to studying disease pathogenesis and developing therapies has been the lack of a small-animal model. We previously reported isolation of a DENV2 strain, obtained by passaging a clinical isolate between mosquito cells and mice, that caused severe DENV disease in mice and contained multiple mutations, including many in domain II of the envelope (E) protein. Here, we describe a recombinant virus, differing from the non-mouse-passaged virus by two mutations in the E protein, that induces vascular leakage and tumor necrosis factor alpha (TNF-alpha)-mediated lethality, while the non-mouse-passaged virus causes paralysis. This recombinant virus has a weaker affinity for heparan sulfate, resulting in an increased serum half-life, higher systemic viral loads, and high levels of TNF-alpha in the serum of infected mice. These results exemplify the role of the E protein in modulating virion clearance and connect the effect of clearance on the systemic viral loads responsible for severe disease manifestations.

Hydrogen-Deuterium Exchange of Lipoxygenase Uncovers a Relationship between Distal, Solvent Exposed Protein Motions and the Thermal Activation Barrier for Catalytic Proton-Coupled Electron Tunneling.

Defining specific pathways for efficient heat transfer from protein-solvent interfaces to their active sites represents one of the compelling and timely challenges in our quest for a physical description of the origins of enzyme catalysis. Enzymatic hydrogen tunneling reactions constitute excellent systems in which to validate experimental approaches to this important question, given the inherent temperature independence of quantum mechanical wave function overlap. Herein, we present the application of hydrogen-deuterium exchange coupled to mass spectrometry toward the spatial resolution of protein motions that can be related to an enzyme's catalytic parameters. Employing the proton-coupled electron transfer reaction of soybean lipoxygenase as proof of principle, we first corroborate the impact of active site mutations on increased local flexibility and, second, uncover a solvent-exposed loop, 15-34 Å from the reactive ferric center whose temperature-dependent motions are demonstrated to mirror the enthalpic barrier for catalytic C-H bond cleavage. A network that connects this surface loop to the active site is structurally identified and supported by changes in kinetic parameters that result from site-specific mutations.

Mycobacterium tuberculosis RpfE crystal structure reveals a positively charged catalytic cleft.

Resuscitation promoting factor (Rpf) proteins, which hydrolyze the sugar chains in cell-wall peptidoglycan (PG), play key roles in prokaryotic cell elongation, division, and escape from dormancy to vegetative growth. Like other bacteria, Mycobacterium tuberculosis (Mtb) expresses multiple Rpfs, none of which is individually essential. This redundancy has left unclear the distinct functions of the different Rpfs. To explore the distinguishing characteristics of the five Mtb Rpfs, we determined the crystal structure of the RpfE catalytic domain. The protein adopts the characteristic Rpf fold, but the catalytic cleft is narrower compared to Mtb RpfB. Also in contrast to RpfB, in which the substrate-binding surfaces are negatively charged, the corresponding RpfE catalytic pocket and predicted peptide-binding sites are more positively charged at neutral pH. The complete reversal of the electrostatic potential of the substrate-binding site suggests that the different Rpfs function optimally at different pHs or most efficiently hydrolyze different micro-domains of PG. These studies provide insights into the molecular determinants of the evolution of functional specialization in Rpfs.

Subfamily-specific adaptations in the structures of two penicillin-binding proteins from Mycobacterium tuberculosis.

Beta-lactam antibiotics target penicillin-binding proteins including several enzyme classes essential for bacterial cell-wall homeostasis. To better understand the functional and inhibitor-binding specificities of penicillin-binding proteins from the pathogen, Mycobacterium tuberculosis, we carried out structural and phylogenetic analysis of two predicted D,D-carboxypeptidases, Rv2911 and Rv3330. Optimization of Rv2911 for crystallization using directed evolution and the GFP folding reporter method yielded a soluble quadruple mutant. Structures of optimized Rv2911 bound to phenylmethylsulfonyl fluoride and Rv3330 bound to meropenem show that, in contrast to the nonspecific inhibitor, meropenem forms an extended interaction with the enzyme along a conserved surface. Phylogenetic analysis shows that Rv2911 and Rv3330 belong to different clades that emerged in Actinobacteria and are not represented in model organisms such as Escherichia coli and Bacillus subtilis. Clade-specific adaptations allow these enzymes to fulfill distinct physiological roles despite strict conservation of core catalytic residues. The characteristic differences include potential protein-protein interaction surfaces and specificity-determining residues surrounding the catalytic site. Overall, these structural insights lay the groundwork to develop improved beta-lactam therapeutics for tuberculosis.

Structure of the sensor domain of Mycobacterium tuberculosis PknH receptor kinase reveals a conserved binding cleft.

Since their discovery over 20 years ago, eukaryotic-like transmembrane receptor Ser/Thr protein kinases (STPKs) have been shown to play critical roles in the virulence, growth, persistence, and reactivation of many bacteria. Information regarding the signals transmitted by these proteins, however, remains scarce. To enhance understanding of the basis for STPK receptor signaling, we determined the 1.7-Å-resolution crystal structure of the extracellular sensor domain of the Mycobacterium tuberculosis receptor STPK, PknH (Rv1266c). The PknH sensor domain adopts an unanticipated fold containing two intramolecular disulfide bonds and a large hydrophobic and polar cleft. The residues lining the cleft and those surrounding the disulfide bonds are conserved. These results suggest that PknH binds a small-molecule ligand that signals by changing the location or quaternary structure of the kinase domain.