A divergent protein kinase A regulatory subunit essential for morphogenesis of the human pathogen Leishmania.

Parasitic protozoa of the genus Leishmania cycle between the phagolysosome of mammalian macrophages, where they reside as rounded intracellular amastigotes, and the midgut of female sand flies, which they colonize as elongated extracellular promastigotes. Previous studies indicated that protein kinase A (PKA) plays an important role in the initial steps of promastigote differentiation into amastigotes. Here, we describe a novel regulatory subunit of PKA (which we have named PKAR3) that is unique to Leishmania and most (but not all) other Kinetoplastidae. PKAR3 is localized to subpellicular microtubules (SPMT) in the cell cortex, where it recruits a specific catalytic subunit (PKAC3). Promastigotes of pkar3 or pkac3 null mutants lose their elongated shape and become rounded but remain flagellated. Truncation of an N-terminal formin homology (FH)-like domain of PKAR3 results in its detachment from the SPMT, also leading to rounded promastigotes. Thus, the tethering of PKAC3 via PKAR3 at the cell cortex is essential for maintenance of the elongated shape of promastigotes. This role of PKAR3 is reminiscent of PKARIß and PKARIIß binding to microtubules of mammalian neurons, which is essential for the elongation of dendrites and axons, respectively. Interestingly, PKAR3 binds nucleoside analogs, but not cAMP, with a high affinity similar to the PKAR1 isoform of Trypanosoma. We propose that these early-diverged protists have re-purposed PKA for a novel signaling pathway that spatiotemporally controls microtubule remodeling and cell shape.

Enolase inhibitors as therapeutic leads for Naegleria fowleri infection.

Infections with the pathogenic free-living amoebae Naegleria fowleri can lead to life-threatening illnesses including catastrophic primary amebic meningoencephalitis (PAM). Efficacious treatment options for these infections are lacking and the mortality rate remains >95% in the US. Glycolysis is very important for the infectious trophozoite lifecycle stage and inhibitors of glucose metabolism have been found to be toxic to the pathogen. Recently, human enolase 2 (ENO2) phosphonate inhibitors have been developed as lead agents to treat glioblastoma multiforme (GBM). These compounds, which cure GBM in a rodent model, are well-tolerated in mammals because enolase 1 (ENO1) is the predominant isoform used systemically. Here, we describe findings that demonstrate that these agents are potent inhibitors of N. fowleri ENO ( Nf ENO) and are lethal to amoebae. In particular, (1-hydroxy-2-oxopiperidin-3-yl) phosphonic acid (HEX) was a potent enzyme inhibitor (IC 50 value of 0.14 ± 0.04 µM) that was toxic to trophozoites (EC 50 value of 0.21 ± 0.02 µM) while the reported CC 50 was >300 µM. Molecular docking simulation revealed that HEX binds strongly to the active site of Nf ENO with a binding affinity of -8.6 kcal/mol. Metabolomic studies of parasites treated with HEX revealed a 4.5 to 78-fold accumulation of glycolytic intermediates upstream of Nf ENO. Last, nasal instillation of HEX increased longevity of amoebae-infected rodents. Two days after infection, animals were treated for 10 days with 3 mg/kg HEX, followed by one week of observation. At the conclusion of the experiment, eight of 12 HEX-treated animals remained alive (resulting in an indeterminable median survival time) while one of 12 vehicle-treated rodents remained, yielding a median survival time of 10.9 days. Brains of six of the eight survivors were positive for amoebae, suggesting the agent at the tested dose suppressed, but did not eliminate, infection. These findings suggest that HEX is a promising lead for the treatment of PAM.

Shotgun Kinetic Target-Guided Synthesis Approach Enables the Discovery of Small-Molecule Inhibitors against Pathogenic Free-Living Amoeba Glucokinases.

Pathogenic free-living amoebae (pFLA) can cause life-threatening central nervous system (CNS) infections and warrant the investigation of new chemical agents to combat the rise of infection from these pathogens. Naegleria fowleri glucokinase (NfGlck), a key metabolic enzyme involved in generating glucose-6-phosphate, was previously identified as a potential target due to its limited sequence similarity with human Glck (HsGlck). Herein, we used our previously demonstrated multifragment kinetic target-guided synthesis (KTGS) screening strategy to identify inhibitors against pFLA glucokinases. Unlike the majority of previous KTGS reports, our current study implements a "shotgun" approach, where fragments were not biased by predetermined binding potentials. The study resulted in the identification of 12 inhibitors against 3 pFLA glucokinase enzymes─NfGlck, Balamuthia mandrillaris Glck (BmGlck), and Acanthamoeba castellanii Glck (AcGlck). This work demonstrates the utility of KTGS to identify small-molecule binders for biological targets where resolved X-ray crystal structures are not readily accessible.

Characterization of Glucokinases from Pathogenic Free-Living Amoebae.

Infection with pathogenic free-living amoebae, including Naegleria fowleri, Acanthamoeba spp., and Balamuthia mandrillaris, can lead to life-threatening illnesses, primarily because of catastrophic central nervous system involvement. Efficacious treatment options for these infections are lacking, and the mortality rate due to infection is high. Previously, we evaluated the N. fowleri glucokinase (NfGlck) as a potential target for therapeutic intervention, as glucose metabolism is critical for in vitro viability. Here, we extended these studies to the glucokinases from two other pathogenic free-living amoebae, including Acanthamoeba castellanii (AcGlck) and B. mandrillaris (BmGlck). While these enzymes are similar (49.3% identical at the amino acid level), they have distinct kinetic properties that distinguish them from each other. For ATP, AcGlck and BmGlck have apparent Km values of 472.5 and 41.0 µM, while Homo sapiens Glck (HsGlck) has a value of 310 µM. Both parasite enzymes also have a higher apparent affinity for glucose than the human counterpart, with apparent Km values of 45.9 µM (AcGlck) and 124 µM (BmGlck) compared to ~8 mM for HsGlck. Additionally, AcGlck and BmGlck differ from each other and other Glcks in their sensitivity to small molecule inhibitors, suggesting that inhibitors with pan-amoebic activity could be challenging to generate.

The transcriptome of Balamuthia mandrillaris trophozoites for structure-guided drug design.

Balamuthia mandrillaris, a pathogenic free-living amoeba, causes cutaneous skin lesions as well as granulomatous amoebic encephalitis, a 'brain-eating' disease. As with the other known pathogenic free-living amoebas (Naegleria fowleri and Acanthamoeba species), drug discovery efforts to combat Balamuthia infections of the central nervous system are sparse; few targets have been validated or characterized at the molecular level, and little is known about the biochemical pathways necessary for parasite survival. Current treatments of encephalitis due to B. mandrillaris lack efficacy, leading to case fatality rates above 90%. Using our recently published methodology to discover potential drugs against pathogenic amoebas, we screened a collection of 85 compounds with known antiparasitic activity and identified 59 compounds that impacted the growth of Balamuthia trophozoites at concentrations below 220 µM. Since there is no fully annotated genome or proteome of B. mandrillaris, we sequenced and assembled its transcriptome from a high-throughput RNA-sequencing (RNA-Seq) experiment and located the coding sequences of the genes potentially targeted by the growth inhibitors from our compound screens. We determined the sequence of 17 of these target genes and obtained expression clones for 15 that we validated by direct sequencing. These will be used in the future in combination with the identified hits in structure guided drug discovery campaigns to develop new approaches for the treatment of Balamuthia infections.

In silico detection of SARS-CoV-2 specific B-cell epitopes and validation in ELISA for serological diagnosis of COVID-19.

Rapid generation of diagnostics is paramount to understand epidemiology and to control the spread of emerging infectious diseases such as COVID-19. Computational methods to predict serodiagnostic epitopes that are specific for the pathogen could help accelerate the development of new diagnostics. A systematic survey of 27 SARS-CoV-2 proteins was conducted to assess whether existing B-cell epitope prediction methods, combined with comprehensive mining of sequence databases and structural data, could predict whether a particular protein would be suitable for serodiagnosis. Nine of the predictions were validated with recombinant SARS-CoV-2 proteins in the ELISA format using plasma and sera from patients with SARS-CoV-2 infection, and a further 11 predictions were compared to the recent literature. Results appeared to be in agreement with 12 of the predictions, in disagreement with 3, while a further 5 were deemed inconclusive. We showed that two of our top five candidates, the N-terminal fragment of the nucleoprotein and the receptor-binding domain of the spike protein, have the highest sensitivity and specificity and signal-to-noise ratio for detecting COVID-19 sera/plasma by ELISA. Mixing the two antigens together for coating ELISA plates led to a sensitivity of 94% (N = 80 samples from persons with RT-PCR confirmed SARS-CoV-2 infection), and a specificity of 97.2% (N = 106 control samples).

Chromatin-Associated Protein Complexes Link DNA Base J and Transcription Termination in Leishmania.

Unlike most other eukaryotes, Leishmania and other trypanosomatid protozoa have largely eschewed transcriptional control of gene expression, relying instead on posttranscriptional regulation of mRNAs derived from polycistronic transcription units (PTUs). In these parasites, a novel modified nucleotide base (ß-d-glucopyranosyloxymethyluracil) known as J plays a critical role in ensuring that transcription termination occurs only at the end of each PTU, rather than at the polyadenylation sites of individual genes. To further understand the biology of J-associated processes, we used tandem affinity purification (TAP) tagging and mass spectrometry to reveal proteins that interact with the glucosyltransferase performing the final step in J synthesis. These studies identified four proteins reminiscent of subunits in the PTW/PP1 complex that controls transcription termination in higher eukaryotes. Moreover, bioinformatic analyses identified the DNA-binding subunit of Leishmania PTW/PP1 as a novel J-binding protein (JBP3), which is also part of another complex containing proteins with domains suggestive of a role in chromatin modification/remodeling. Additionally, JBP3 associates (albeit transiently and/or indirectly) with the trypanosomatid equivalent of the PAF1 complex involved in the regulation of transcription in other eukaryotes. The downregulation of JBP3 expression levels in Leishmania resulted in a substantial increase in transcriptional readthrough at the 3' end of most PTUs. We propose that JBP3 recruits one or more of these complexes to the J-containing regions at the end of PTUs, where they halt the progression of the RNA polymerase. This decoupling of transcription termination from the splicing of individual genes enables the parasites' unique reliance on polycistronic transcription and posttranscriptional regulation of gene expression.IMPORTANCELeishmania parasites cause a variety of serious human diseases, with no effective vaccine and emerging resistance to current drug therapy. We have previously shown that a novel DNA base called J is critical for transcription termination at the ends of the polycistronic gene clusters that are a hallmark of Leishmania and related trypanosomatids. Here, we describe a new J-binding protein (JBP3) associated with three different protein complexes that are reminiscent of those involved in the control of transcription in other eukaryotes. However, the parasite complexes have been reprogrammed to regulate transcription and gene expression in trypanosomatids differently than in the mammalian hosts, providing new opportunities to develop novel chemotherapeutic treatments against these important pathogens.

Multiplexed Functional Assessment of Genetic Variants in CARD11.

Genetic testing has increased the number of variants identified in disease genes, but the diagnostic utility is limited by lack of understanding variant function. CARD11 encodes an adaptor protein that expresses dominant-negative and gain-of-function variants associated with distinct immunodeficiencies. Here, we used a "cloning-free" saturation genome editing approach in a diploid cell line to simultaneously score 2,542 variants for decreased or increased function in the region of CARD11 associated with immunodeficiency. We also described an exon-skipping mechanism for CARD11 dominant-negative activity. The classification of reported clinical variants was sensitive (94.6%) and specific (88.9%), which rendered the data immediately useful for interpretation of seven coding and splicing variants implicated in immunodeficiency found in our clinic. This approach is generalizable for variant interpretation in many other clinically actionable genes, in any relevant cell type.

A doublecortin-domain protein of Toxoplasma and its orthologues bind to and modify the structure and organization of tubulin polymers.

BACKGROUND: TgDCX is a doublecortin-domain protein associated with the conoid fibers, a set of strongly curved non-tubular tubulin-polymers in Toxoplasma. TgDCX deletion impairs conoid structure and parasite invasion. TgDCX contains two tubulin-binding domains: a partial P25α and the DCX/doublecortin domain. Orthologues are found in apicomplexans and their free-living relatives Chromera and Vitrella. RESULTS: We report that isolated TgDCX-containing conoid fibers retain their pronounced curvature, but loss of TgDCX destabilizes the fibers. We crystallized and determined the 3D-structure of the DCX-domain, which is similar to those of human doublecortin and well-conserved among TgDCX orthologues. However, the orthologues vary widely in targeting to the conoid in Toxoplasma and in modulating microtubule organization in Xenopus cells. Several orthologues bind to microtubules in Xenopus cells, but only TgDCX generates short, strongly curved microtubule arcs. EM analysis shows microtubules decorated with TgDCX bundled into rafts, often bordered on one edge by a "C"-shaped incomplete tube. A Chromera orthologue closely mimics TgDCX targeting in Toxoplasma and binds to microtubules in Xenopus cells, but does not generate arcs or "C"-shaped tubes, and fails to rescue the defects of the TgDCX-knockout parasite. CONCLUSIONS: These observations suggest that species-specific features of TgDCX enable it to generate strongly curved tubulin-polymers to support efficient host-cell invasion.

Structures of glyceraldehyde 3-phosphate dehydrogenase in Neisseria gonorrhoeae and Chlamydia trachomatis.

Neisseria gonorrhoeae (Ng) and Chlamydia trachomatis (Ct) are the most commonly reported sexually transmitted bacteria worldwide and usually present as co-infections. Increasing resistance of Ng to currently recommended dual therapy of azithromycin and ceftriaxone presents therapeutic challenges for syndromic management of Ng-Ct co-infections. Development of a safe, effective, and inexpensive dual therapy for Ng-Ct co-infections is an effective strategy for the global control and prevention of these two most prevalent bacterial sexually transmitted infections. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is a validated drug target with two approved drugs for indications other than antibacterials. Nonetheless, any new drugs targeting GAPDH in Ng and Ct must be specific inhibitors of bacterial GAPDH that do not inhibit human GAPDH, and structural information of Ng and Ct GAPDH will aid in finding such selective inhibitors. Here, we report the X-ray crystal structures of Ng and Ct GAPDH. Analysis of the structures demonstrates significant differences in amino acid residues in the active sites of human GAPDH from those of the two bacterial enzymes suggesting design of compounds to selectively inhibit Ng and Ct is possible. We also describe an efficient in vitro assay of recombinant GAPDH enzyme activity amenable to high-throughput drug screening to aid in identifying inhibitory compounds and begin to address selectivity.

Toward a structome of Acinetobacter baumannii drug targets.

Acinetobacter baumannii is well known for causing hospital-associated infections due in part to its intrinsic antibiotic resistance as well as its ability to remain viable on surfaces and resist cleaning agents. In a previous publication, A. baumannii strain AB5075 was studied by transposon mutagenesis and 438 essential gene candidates for growth on rich-medium were identified. The Seattle Structural Genomics Center for Infectious Disease entered 342 of these candidate essential genes into our pipeline for structure determination, in which 306 were successfully cloned into expression vectors, 192 were detectably expressed, 165 screened as soluble, 121 were purified, 52 crystalized, 30 provided diffraction data, and 29 structures were deposited in the Protein Data Bank. Here, we report these structures, compare them with human orthologs where applicable, and discuss their potential as drug targets for antibiotic development against A. baumannii.

Structural diversity in the Mycobacteria DUF3349 superfamily.

A protein superfamily with a "Domain of Unknown Function,", DUF3349 (PF11829), is present predominately in Mycobacterium and Rhodococcus bacterial species suggesting that these proteins may have a biological function unique to these bacteria. We previously reported the inaugural structure of a DUF3349 superfamily member, Mycobacterium tuberculosis Rv0543c. Here, we report the structures determined for three additional DUF3349 proteins: Mycobacterium smegmatis MSMEG_1063 and MSMEG_1066 and Mycobacterium abscessus MAB_3403c. Like Rv0543c, the NMR solution structure of MSMEG_1063 revealed a monomeric five α-helix bundle with a similar overall topology. Conversely, the crystal structure of MSMEG_1066 revealed a five α-helix protein with a strikingly different topology and a tetrameric quaternary structure that was confirmed by size exclusion chromatography. The NMR solution structure of a fourth member of the DUF3349 superfamily, MAB_3403c, with 18 residues missing at the N-terminus, revealed a monomeric α-helical protein with a folding topology similar to the three C-terminal helices in the protomer of the MSMEG_1066 tetramer. These structures, together with a GREMLIN-based bioinformatics analysis of the DUF3349 primary amino acid sequences, suggest two subfamilies within the DUF3349 family. The division of the DUF3349 into two distinct subfamilies would have been lost if structure solution had stopped with the first structure in the DUF3349 family, highlighting the insights generated by solving multiple structures within a protein superfamily. Future studies will determine if the structural diversity at the tertiary and quaternary levels in the DUF3349 protein superfamily have functional roles in Mycobacteria and Rhodococcus species with potential implications for structure-based drug discovery.

The domain architecture of the protozoan protein J-DNA-binding protein 1 suggests synergy between base J DNA binding and thymidine hydroxylase activity.

J-DNA-binding protein 1 (JBP1) contributes to the biosynthesis and maintenance of base J (ß-d-glucosyl-hydroxymethyluracil), an epigenetic modification of thymidine (T) confined to pathogenic protozoa such as Trypanosoma and Leishmania JBP1 has two known functional domains: an N-terminal T hydroxylase (TH) homologous to the 5-methylcytosine hydroxylase domain in TET proteins and a J-DNA-binding domain (JDBD) that resides in the middle of JBP1. Here, we show that removing JDBD from JBP1 results in a soluble protein (Δ-JDBD) with the N- and C-terminal regions tightly associated together in a well-ordered structure. We found that this Δ-JDBD domain retains TH activity in vitro but displays a 15-fold lower apparent rate of hydroxylation compared with JBP1. Small-angle X-ray scattering (SAXS) experiments on JBP1 and JDBD in the presence or absence of J-DNA and on Δ-JDBD enabled us to generate low-resolution three-dimensional models. We conclude that Δ-JDBD, and not the N-terminal region of JBP1 alone, is a distinct folding unit. Our SAXS-based model supports the notion that binding of JDBD specifically to J-DNA can facilitate T hydroxylation 12-14 bp downstream on the complementary strand of the J-recognition site. We postulate that insertion of the JDBD module into the Δ-JDBD scaffold during evolution provided a mechanism that synergized J recognition and T hydroxylation, ensuring inheritance of base J in specific sequence patterns following DNA replication in kinetoplastid parasites.

Enzymatic and Structural Characterization of the Naegleria fowleri Glucokinase.

Infection with the free-living amoeba Naegleria fowleri leads to life-threatening primary amoebic meningoencephalitis. Efficacious treatment options for these infections are limited, and the mortality rate is very high (∼98%). Parasite metabolism may provide suitable targets for therapeutic design. Like most other organisms, glucose metabolism is critical for parasite viability, being required for growth in culture. The first enzyme required for glucose metabolism is typically a hexokinase (HK), which transfers a phosphate from ATP to glucose. The products of this enzyme are required for both glycolysis and the pentose phosphate pathway. However, the N. fowleri genome lacks an obvious HK homolog and instead harbors a glucokinase (Glck). The N. fowleri Glck (NfGlck) shares limited (25%) amino acid identity with the mammalian host enzyme (Homo sapiens Glck), suggesting that parasite-specific inhibitors with anti-amoeba activity can be generated. Following heterologous expression, NfGlck was found to have a limited hexose substrate range, with the greatest activity observed with glucose. The enzyme had apparent Km values of 42.5 ± 7.3 µM and 141.6 ± 9.9 µM for glucose and ATP, respectively. The NfGlck structure was determined and refined to 2.2-Å resolution, revealing that the enzyme shares greatest structural similarity with the Trypanosoma cruzi Glck. These similarities include binding modes and binding environments for substrates. To identify inhibitors of NfGlck, we screened a small collection of inhibitors of glucose-phosphorylating enzymes and identified several small molecules with 50% inhibitory concentration values of <1 µM that may prove useful as hit chemotypes for further leads and therapeutic development against N. fowleri.

Mycobacterium tuberculosis Rv3651 is a triple sensor-domain protein.

The genome of the human pathogen Mycobacterium tuberculosis (Mtb) encodes ∼4,400 proteins, but one third of them have unknown functions. We solved the crystal structure of Rv3651, a hypothetical protein with no discernible similarity to proteins with known function. Rv3651 has a three-domain architecture that combines one cGMP-specific phosphodiesterases, adenylyl cyclases and FhlA (GAF) domain and two Per-ARNT-Sim (PAS) domains. GAF and PAS domains are sensor domains that are typically linked to signaling effector molecules. Unlike these sensor-effector proteins, Rv3651 is an unusual sensor domain-only protein with highly divergent sequence. The structure suggests that Rv3651 integrates multiple different signals and serves as a scaffold to facilitate signal transfer.

Calmodulin-like proteins localized to the conoid regulate motility and cell invasion by Toxoplasma gondii.

Toxoplasma gondii contains an expanded number of calmodulin (CaM)-like proteins whose functions are poorly understood. Using a combination of CRISPR/Cas9-mediated gene editing and a plant-like auxin-induced degron (AID) system, we examined the roles of three apically localized CaMs. CaM1 and CaM2 were individually dispensable, but loss of both resulted in a synthetic lethal phenotype. CaM3 was refractory to deletion, suggesting it is essential. Consistent with this prediction auxin-induced degradation of CaM3 blocked growth. Phenotypic analysis revealed that all three CaMs contribute to parasite motility, invasion, and egress from host cells, and that they act downstream of microneme and rhoptry secretion. Super-resolution microscopy localized all three CaMs to the conoid where they overlap with myosin H (MyoH), a motor protein that is required for invasion. Biotinylation using BirA fusions with the CaMs labeled a number of apical proteins including MyoH and its light chain MLC7, suggesting they may interact. Consistent with this hypothesis, disruption of MyoH led to degradation of CaM3, or redistribution of CaM1 and CaM2. Collectively, our findings suggest these CaMs may interact with MyoH to control motility and cell invasion.

Structure of acid deoxyribonuclease.

Deoxyribonuclease II (DNase II) is also known as acid deoxyribonuclease because it has optimal activity at the low pH environment of lysosomes where it is typically found in higher eukaryotes. Interestingly, DNase II has also been identified in a few genera of bacteria and is believed to have arisen via horizontal transfer. Here, we demonstrate that recombinant Burkholderia thailandensis DNase II is highly active at low pH in the absence of divalent metal ions, similar to eukaryotic DNase II. The crystal structure of B. thailandensis DNase II shows a dimeric quaternary structure which appears capable of binding double-stranded DNA. Each monomer of B. thailandensis DNase II exhibits a similar overall fold as phospholipase D (PLD), phosphatidylserine synthase (PSS) and tyrosyl-DNA phosphodiesterase (TDP), and conserved catalytic residues imply a similar mechanism. The structural and biochemical data presented here provide insights into the atomic structure and catalytic mechanism of DNase II.

Mycofactocin-associated mycobacterial dehydrogenases with non-exchangeable NAD cofactors.

During human infection, Mycobacterium tuberculosis (Mtb) survives the normally bacteriocidal phagosome of macrophages. Mtb and related species may be able to combat this harsh acidic environment which contains reactive oxygen species due to the mycobacterial genomes encoding a large number of dehydrogenases. Typically, dehydrogenase cofactor binding sites are open to solvent, which allows NAD/NADH exchange to support multiple turnover. Interestingly, mycobacterial short chain dehydrogenases/reductases (SDRs) within family TIGR03971 contain an insertion at the NAD binding site. Here we present crystal structures of 9 mycobacterial SDRs in which the insertion buries the NAD cofactor except for a small portion of the nicotinamide ring. Line broadening and STD-NMR experiments did not show NAD or NADH exchange on the NMR timescale. STD-NMR demonstrated binding of the potential substrate carveol, the potential product carvone, the inhibitor tricyclazol, and an external redox partner 2,6-dichloroindophenol (DCIP). Therefore, these SDRs appear to contain a non-exchangeable NAD cofactor and may rely on an external redox partner, rather than cofactor exchange, for multiple turnover. Incidentally, these genes always appear in conjunction with the mftA gene, which encodes the short peptide MftA, and with other genes proposed to convert MftA into the external redox partner mycofactocin.

The Rickettsia type IV secretion system: unrealized complexity mired by gene family expansion.

Many prokaryotes utilize type IV secretion systems (T4SSs) to translocate substrates (e.g. nucleoprotein, DNA, protein) across the cell envelope, and/or to elaborate surface structures (i.e. pili or adhesins). Among eight distinct T4SS classes, P-T4SSs are typified by the Agrobacterium tumefaciens vir T4SS, which is comprised of 12 scaffold components (VirB1-VirB11, VirD4). While most P-T4SSs include all 12 Vir proteins, some differ from the vir archetype by either containing additional scaffold components not analogous to Vir proteins or lacking one or more of the Vir proteins. In a special case, the Rickettsiales vir homolog (rvh) P-T4SS comprises unprecedented gene family expansion. rvh contains three families of gene duplications (rvhB9, rvhB8, rvhB4): RvhB9,8,4-I are conserved relative to equivalents in other P-T4SSs, while RvhB9,8,4-II have evolved atypical features that deviate substantially from other homologs. Furthermore, rvh contains five VirB6-like genes (rvhB6a-e), which are tandemly arrayed and contain large N- and C-terminal extensions. Our work herein focuses on the complexity underpinned by rvh gene family expansion. Furthermore, we describe an RvhB10 insertion, which occurs in a region that forms the T4SS pore. The significance of these curious properties to rvh structure and function is evaluated, shedding light on a highly complex T4SS.

Structural Insight into How Bacteria Prevent Interference between Multiple Divergent Type IV Secretion Systems.

UNLABELLED: Prokaryotes use type IV secretion systems (T4SSs) to translocate substrates (e.g., nucleoprotein, DNA, and protein) and/or elaborate surface structures (i.e., pili or adhesins). Bacterial genomes may encode multiple T4SSs, e.g., there are three functionally divergent T4SSs in some Bartonella species (vir, vbh, and trw). In a unique case, most rickettsial species encode a T4SS (rvh) enriched with gene duplication. Within single genomes, the evolutionary and functional implications of cross-system interchangeability of analogous T4SS protein components remains poorly understood. To lend insight into cross-system interchangeability, we analyzed the VirB8 family of T4SS channel proteins. Crystal structures of three VirB8 and two TrwG Bartonella proteins revealed highly conserved C-terminal periplasmic domain folds and dimerization interfaces, despite tremendous sequence divergence. This implies remarkable structural constraints for VirB8 components in the assembly of a functional T4SS. VirB8/TrwG heterodimers, determined via bacterial two-hybrid assays and molecular modeling, indicate that differential expression of trw and vir systems is the likely barrier to VirB8-TrwG interchangeability. We also determined the crystal structure of Rickettsia typhi RvhB8-II and modeled its coexpressed divergent paralog RvhB8-I. Remarkably, while RvhB8-I dimerizes and is structurally similar to other VirB8 proteins, the RvhB8-II dimer interface deviates substantially from other VirB8 structures, potentially preventing RvhB8-I/RvhB8-II heterodimerization. For the rvh T4SS, the evolution of divergent VirB8 paralogs implies a functional diversification that is unknown in other T4SSs. Collectively, our data identify two different constraints (spatiotemporal for Bartonella trw and vir T4SSs and structural for rvh T4SSs) that mediate the functionality of multiple divergent T4SSs within a single bacterium. IMPORTANCE: Assembly of multiprotein complexes at the right time and at the right cellular location is a fundamentally important task for any organism. In this respect, bacteria that express multiple analogous type IV secretion systems (T4SSs), each composed of around 12 different components, face an overwhelming complexity. Our work here presents the first structural investigation on factors regulating the maintenance of multiple T4SSs within a single bacterium. The structural data imply that the T4SS-expressing bacteria rely on two strategies to prevent cross-system interchangeability: (i) tight temporal regulation of expression or (ii) rapid diversification of the T4SS components. T4SSs are ideal drug targets provided that no analogous counterparts are known from eukaryotes. Drugs targeting the barriers to cross-system interchangeability (i.e., regulators) could dysregulate the structural and functional independence of discrete systems, potentially creating interference that prevents their efficient coordination throughout bacterial infection.