Genetic manipulation of Patescibacteria provides mechanistic insights into microbial dark matter and the epibiotic lifestyle.

Patescibacteria, also known as the candidate phyla radiation (CPR), are a diverse group of bacteria that constitute a disproportionately large fraction of microbial dark matter. Its few cultivated members, belonging mostly to Saccharibacteria, grow as epibionts on host Actinobacteria. Due to a lack of suitable tools, the genetic basis of this lifestyle and other unique features of Patescibacteira remain unexplored. Here, we show that Saccharibacteria exhibit natural competence, and we exploit this property for their genetic manipulation. Imaging of fluorescent protein-labeled Saccharibacteria provides high spatiotemporal resolution of phenomena accompanying epibiotic growth, and a transposon-insertion sequencing (Tn-seq) genome-wide screen reveals the contribution of enigmatic Saccharibacterial genes to growth on their hosts. Finally, we leverage metagenomic data to provide cutting-edge protein structure-based bioinformatic resources that support the strain Southlakia epibionticum and its corresponding host, Actinomyces israelii, as a model system for unlocking the molecular underpinnings of the epibiotic lifestyle.

A systems-level gene regulatory network model for Plasmodium falciparum.

Many of the gene regulatory processes of Plasmodium falciparum, the deadliest malaria parasite, remain poorly understood. To develop a comprehensive guide for exploring this organism's gene regulatory network, we generated a systems-level model of P. falciparum gene regulation using a well-validated, machine-learning approach for predicting interactions between transcription regulators and their targets. The resulting network accurately predicts expression levels of transcriptionally coherent gene regulatory programs in independent transcriptomic data sets from parasites collected by different research groups in diverse laboratory and field settings. Thus, our results indicate that our gene regulatory model has predictive power and utility as a hypothesis-generating tool for illuminating clinically relevant gene regulatory mechanisms within P. falciparum. Using the set of regulatory programs we identified, we also investigated correlates of artemisinin resistance based on gene expression coherence. We report that resistance is associated with incoherent expression across many regulatory programs, including those controlling genes associated with erythrocyte-host engagement. These results suggest that parasite populations with reduced artemisinin sensitivity are more transcriptionally heterogenous. This pattern is consistent with a model where the parasite utilizes bet-hedging strategies to diversify the population, rendering a subpopulation more able to navigate drug treatment.

Reconsidering the in vivo functions of Clostridial Stickland amino acid fermentations.

Stickland amino acid fermentations occur primarily among species of Clostridia. An ancient form of metabolism, Stickland fermentations use amino acids as electron acceptors in the absence of stronger oxidizing agents and provide metabolic capabilities to support growth when other fermentable substrates, such as carbohydrates, are lacking. The reactions were originally described as paired fermentations of amino acid electron donors, such as the branched-chain amino acids, with recipients that include proline and glycine. We present a redox-focused view of Stickland metabolism following electron flow through metabolically diverse oxidative reactions and the defined-substrate reductase systems, including for proline and glycine, and the role of dual redox pathways for substrates such as leucine and ornithine. Genetic studies and Environment and Gene Regulatory Interaction Network (EGRIN) models for the pathogen Clostridioides difficile have improved our understanding of the regulation and metabolic recruitment of these systems, and their functions in modulating inter-species interactions within host-pathogen-commensal systems and uses in industrial and environmental applications.

The type VI secretion system of Xanthomonas phaseoli pv. manihotis is involved in virulence and in vitro motility.

BACKGROUND: The type VI protein secretion system (T6SS) is important in diverse cellular processes in Gram-negative bacteria, including interactions with other bacteria and with eukaryotic hosts. In this study we analyze the evolution of the T6SS in the genus Xanthomonas and evaluate its importance of the T6SS for virulence and in vitro motility in Xanthomonas phaseoli pv. manihotis (Xpm), the causal agent of bacterial blight in cassava (Manihot esculenta). We delineate the organization of the T6SS gene clusters in Xanthomonas and then characterize proteins of this secretion system in Xpm strain CIO151. RESULTS: We describe the presence of three different clusters in the genus Xanthomonas that vary in their organization and degree of synteny between species. Using a gene knockout strategy, we also found that vgrG and hcp are required for maximal aggressiveness of Xpm on cassava plants while clpV is important for both motility and maximal aggressiveness. CONCLUSION: We characterized the T6SS in 15 different strains in Xanthomonas and our phylogenetic analyses suggest that the T6SS might have been acquired by a very ancient event of horizontal gene transfer and maintained through evolution, hinting at their importance for the adaptation of Xanthomonas to their hosts. Finally, we demonstrated that the T6SS of Xpm is functional, and significantly contributes to virulence and motility. This is the first experimental study that demonstrates the role of the T6SS in the Xpm-cassava interaction and the T6SS organization in the genus Xanthomonas.

Path-seq identifies an essential mycolate remodeling program for mycobacterial host adaptation.

The success of Mycobacterium tuberculosis (MTB) stems from its ability to remain hidden from the immune system within macrophages. Here, we report a new technology (Path-seq) to sequence miniscule amounts of MTB transcripts within up to million-fold excess host RNA Using Path-seq and regulatory network analyses, we have discovered a novel transcriptional program for in vivo mycobacterial cell wall remodeling when the pathogen infects alveolar macrophages in mice. We have discovered that MadR transcriptionally modulates two mycolic acid desaturases desA1/desA2 to initially promote cell wall remodeling upon in vitro macrophage infection and, subsequently, reduces mycolate biosynthesis upon entering dormancy. We demonstrate that disrupting MadR program is lethal to diverse mycobacteria making this evolutionarily conserved regulator a prime antitubercular target for both early and late stages of infection.

An experimentally supported model of the Bacillus subtilis global transcriptional regulatory network.

Organisms from all domains of life use gene regulation networks to control cell growth, identity, function, and responses to environmental challenges. Although accurate global regulatory models would provide critical evolutionary and functional insights, they remain incomplete, even for the best studied organisms. Efforts to build comprehensive networks are confounded by challenges including network scale, degree of connectivity, complexity of organism-environment interactions, and difficulty of estimating the activity of regulatory factors. Taking advantage of the large number of known regulatory interactions in Bacillus subtilis and two transcriptomics datasets (including one with 38 separate experiments collected specifically for this study), we use a new combination of network component analysis and model selection to simultaneously estimate transcription factor activities and learn a substantially expanded transcriptional regulatory network for this bacterium. In total, we predict 2,258 novel regulatory interactions and recall 74% of the previously known interactions. We obtained experimental support for 391 (out of 635 evaluated) novel regulatory edges (62% accuracy), thus significantly increasing our understanding of various cell processes, such as spore formation.

3D Hydrogel Culture System Recapitulates Key Tuberculosis Phenotypes and Demonstrates Pyrazinamide Efficacy.

The mortality caused by tuberculosis (TB) infections is a global concern, and there is a need to improve understanding of the disease. Current in vitro infection models to study the disease have limitations such as short investigation durations and divergent transcriptional signatures. This study aims to overcome these limitations by developing a 3D collagen culture system that mimics the biomechanical and extracellular matrix (ECM) of lung microenvironment (collagen fibers, stiffness comparable to in vivo conditions) as the infection primarily manifests in the lungs. The system incorporates Mycobacterium tuberculosis (Mtb) infected human THP-1 or primary monocytes/macrophages. Dual RNA sequencing reveals higher mammalian gene expression similarity with patient samples than 2D macrophage infections. Similarly, bacterial gene expression more accurately recapitulates in vivo gene expression patterns compared to bacteria in 2D infection models. Key phenotypes observed in humans, such as foamy macrophages and mycobacterial cords, are reproduced in the model. This biomaterial system overcomes challenges associated with traditional platforms by modulating immune cells and closely mimicking in vivo infection conditions, including showing efficacy with clinically relevant concentrations of anti-TB drug pyrazinamide, not seen in any other in vitro infection model, making it reliable and readily adoptable for tuberculosis studies and drug screening.

Genome sequence of Xanthomonas fuscans subsp. fuscans strain 4834-R reveals that flagellar motility is not a general feature of xanthomonads.

BACKGROUND: Xanthomonads are plant-associated bacteria responsible for diseases on economically important crops. Xanthomonas fuscans subsp. fuscans (Xff) is one of the causal agents of common bacterial blight of bean. In this study, the complete genome sequence of strain Xff 4834-R was determined and compared to other Xanthomonas genome sequences. RESULTS: Comparative genomics analyses revealed core characteristics shared between Xff 4834-R and other xanthomonads including chemotaxis elements, two-component systems, TonB-dependent transporters, secretion systems (from T1SS to T6SS) and multiple effectors. For instance a repertoire of 29 Type 3 Effectors (T3Es) with two Transcription Activator-Like Effectors was predicted. Mobile elements were associated with major modifications in the genome structure and gene content in comparison to other Xanthomonas genomes. Notably, a deletion of 33 kbp affects flagellum biosynthesis in Xff 4834-R. The presence of a complete flagellar cluster was assessed in a collection of more than 300 strains representing different species and pathovars of Xanthomonas. Five percent of the tested strains presented a deletion in the flagellar cluster and were non-motile. Moreover, half of the Xff strains isolated from the same epidemic than 4834-R was non-motile and this ratio was conserved in the strains colonizing the next bean seed generations. CONCLUSIONS: This work describes the first genome of a Xanthomonas strain pathogenic on bean and reports the existence of non-motile xanthomonads belonging to different species and pathovars. Isolation of such Xff variants from a natural epidemic may suggest that flagellar motility is not a key function for in planta fitness.

Disrupting the ArcA Regulatory Network Amplifies the Fitness Cost of Tetracycline Resistance in Escherichia coli.

There is an urgent need for strategies to discover secondary drugs to prevent or disrupt antimicrobial resistance (AMR), which is causing >700,000 deaths annually. Here, we demonstrate that tetracycline-resistant (TetR) Escherichia coli undergoes global transcriptional and metabolic remodeling, including downregulation of tricarboxylic acid cycle and disruption of redox homeostasis, to support consumption of the proton motive force for tetracycline efflux. Using a pooled genome-wide library of single-gene deletion strains, at least 308 genes, including four transcriptional regulators identified by our network analysis, were confirmed as essential for restoring the fitness of TetR E. coli during treatment with tetracycline. Targeted knockout of ArcA, identified by network analysis as a master regulator of this new compensatory physiological state, significantly compromised fitness of TetR E. coli during tetracycline treatment. A drug, sertraline, which generated a similar metabolome profile as the arcA knockout strain, also resensitized TetR E. coli to tetracycline. We discovered that the potentiating effect of sertraline was eliminated upon knocking out arcA, demonstrating that the mechanism of potential synergy was through action of sertraline on the tetracycline-induced ArcA network in the TetR strain. Our findings demonstrate that therapies that target mechanistic drivers of compensatory physiological states could resensitize AMR pathogens to lost antibiotics. IMPORTANCE Antimicrobial resistance (AMR) is projected to be the cause of >10 million deaths annually by 2050. While efforts to find new potent antibiotics are effective, they are expensive and outpaced by the rate at which new resistant strains emerge. There is desperate need for a rational approach to accelerate the discovery of drugs and drug combinations that effectively clear AMR pathogens and even prevent the emergence of new resistant strains. Using tetracycline-resistant (TetR) Escherichia coli, we demonstrate that gaining resistance is accompanied by loss of fitness, which is restored by compensatory physiological changes. We demonstrate that transcriptional regulators of the compensatory physiologic state are promising drug targets because their disruption increases the susceptibility of TetR E. coli to tetracycline. Thus, we describe a generalizable systems biology approach to identify new vulnerabilities within AMR strains to rationally accelerate the discovery of therapeutics that extend the life span of existing antibiotics.

MtrA modulates Mycobacterium tuberculosis cell division in host microenvironments to mediate intrinsic resistance and drug tolerance.

The success of Mycobacterium tuberculosis (Mtb) is largely attributed to its ability to physiologically adapt and withstand diverse localized stresses within host microenvironments. Here, we present a data-driven model (EGRIN 2.0) that captures the dynamic interplay of environmental cues and genome-encoded regulatory programs in Mtb. Analysis of EGRIN 2.0 shows how modulation of the MtrAB two-component signaling system tunes Mtb growth in response to related host microenvironmental cues. Disruption of MtrAB by tunable CRISPR interference confirms that the signaling system regulates multiple peptidoglycan hydrolases, among other targets, that are important for cell division. Further, MtrA decreases the effectiveness of antibiotics by mechanisms of both intrinsic resistance and drug tolerance. Together, the model-enabled dissection of complex MtrA regulation highlights its importance as a drug target and illustrates how EGRIN 2.0 facilitates discovery and mechanistic characterization of Mtb adaptation to specific host microenvironments within the host.

Genetic manipulation of candidate phyla radiation bacteria provides functional insights into microbial dark matter.

The study of bacteria has yielded fundamental insights into cellular biology and physiology, biotechnological advances and many therapeutics. Yet due to a lack of suitable tools, the significant portion of bacterial diversity held within the candidate phyla radiation (CPR) remains inaccessible to such pursuits. Here we show that CPR bacteria belonging to the phylum Saccharibacteria exhibit natural competence. We exploit this property to develop methods for their genetic manipulation, including the insertion of heterologous sequences and the construction of targeted gene deletions. Imaging of fluorescent protein-labeled Saccharibacteria provides high spatiotemporal resolution of phenomena accompanying epibiotic growth and a transposon insertion sequencing genome-wide screen reveals the contribution of enigmatic Saccharibacterial genes to growth on their Actinobacteria hosts. Finally, we leverage metagenomic data to provide cutting-edge protein structure-based bioinformatic resources that support the strain Southlakia epibionticum and its corresponding host, Actinomyces israelii , as a model system for unlocking the molecular underpinnings of the epibiotic lifestyle.

Genomes-based phylogeny of the genus Xanthomonas.

BACKGROUND: The genus Xanthomonas comprises several plant pathogenic bacteria affecting a wide range of hosts. Despite the economic, industrial and biological importance of Xanthomonas, the classification and phylogenetic relationships within the genus are still under active debate. Some of the relationships between pathovars and species have not been thoroughly clarified, with old pathovars becoming new species. A change in the genus name has been recently suggested for Xanthomonas albilineans, an early branching species currently located in this genus, but a thorough phylogenomic reconstruction would aid in solving these and other discrepancies in this genus. RESULTS: Here we report the results of the genome-wide analysis of DNA sequences from 989 orthologous groups from 17 Xanthomonas spp. genomes available to date, representing all major lineages within the genus. The phylogenetic and computational analyses used in this study have been automated in a Perl package designated Unus, which provides a framework for phylogenomic analyses which can be applied to other datasets at the genomic level. Unus can also be easily incorporated into other phylogenomic pipelines. CONCLUSIONS: Our phylogeny agrees with previous phylogenetic topologies on the genus, but revealed that the genomes of Xanthomonas citri and Xanthomonas fuscans belong to the same species, and that of Xanthomonas albilineans is basal to the joint clade of Xanthomonas and Xylella fastidiosa. Genome reduction was identified in the species Xanthomonas vasicola in addition to the previously identified reduction in Xanthomonas albilineans. Lateral gene transfer was also observed in two gene clusters.

Mathematical models to study the biology of pathogens and the infectious diseases they cause.

Mathematical models have many applications in infectious diseases: epidemiologists use them to forecast outbreaks and design containment strategies; systems biologists use them to study complex processes sustaining pathogens, from the metabolic networks empowering microbial cells to ecological networks in the microbiome that protects its host. Here, we (1) review important models relevant to infectious diseases, (2) draw parallels among models ranging widely in scale. We end by discussing a minimal set of information for a model to promote its use by others and to enable predictions that help us better fight pathogens and the diseases they cause.

Mycobacterium abscessus biofilms produce an extracellular matrix and have a distinct mycolic acid profile.

A non-tuberculous mycobacterium, Mycobacterium abscessus is an emerging opportunistic pathogen associated with difficult to treat pulmonary infections, particularly in patients suffering from cystic fibrosis. It is capable of forming biofilms in vitro that result in an increase of already high levels of antibiotic resistance in this bacterium. Evidence that M. abscessus forms biofilm-like microcolonies in patient lungs and on medical devices further implicated the need to investigate this biofilm in detail. Therefore, in this study we characterized the M. abscessus pellicular biofilm, formed on a liquid-air interface, by studying its molecular composition, and its transcriptional profile in comparison to planktonic cells. Using scanning electron micrographs and fluorescence microscopy, we showed that M. abscessus biofilms produce an extracellular matrix composed of lipids, proteins, carbohydrates and extracellular DNA. Transcriptomic analysis of biofilms revealed an upregulation of pathways involved in the glyoxylate shunt, redox metabolism and mycolic acid biosynthesis. Genes involved in elongation and desaturation of mycolic acids were highly upregulated in biofilms and, mirroring those findings, biochemical analysis of mycolates revealed molecular changes and an increase in mycolic acid chain length. Together these results give us an insight into the complex structure of M. abscessus biofilms, the understanding of which may be adapted for clinical use in treatment of biofilm infections, including strategies for dispersing the extracellular matrix, allowing antibiotics to gain access to bacteria within the biofilm.

Quantitative prediction of conditional vulnerabilities in regulatory and metabolic networks using PRIME.

The ability of Mycobacterium tuberculosis (Mtb) to adopt heterogeneous physiological states underlies its success in evading the immune system and tolerating antibiotic killing. Drug tolerant phenotypes are a major reason why the tuberculosis (TB) mortality rate is so high, with over 1.8 million deaths annually. To develop new TB therapeutics that better treat the infection (faster and more completely), a systems-level approach is needed to reveal the complexity of network-based adaptations of Mtb. Here, we report a new predictive model called PRIME (Phenotype of Regulatory influences Integrated with Metabolism and Environment) to uncover environment-specific vulnerabilities within the regulatory and metabolic networks of Mtb. Through extensive performance evaluations using genome-wide fitness screens, we demonstrate that PRIME makes mechanistically accurate predictions of context-specific vulnerabilities within the integrated regulatory and metabolic networks of Mtb, accurately rank-ordering targets for potentiating treatment with frontline drugs.

Predictive regulatory and metabolic network models for systems analysis of Clostridioides difficile.

We present predictive models for comprehensive systems analysis of Clostridioides difficile, the etiology of pseudomembranous colitis. By leveraging 151 published transcriptomes, we generated an EGRIN model that organizes 90% of C. difficile genes into a transcriptional regulatory network of 297 co-regulated modules, implicating genes in sporulation, carbohydrate transport, and metabolism. By advancing a metabolic model through addition and curation of metabolic reactions including nutrient uptake, we discovered 14 amino acids, diverse carbohydrates, and 10 metabolic genes as essential for C. difficile growth in the intestinal environment. Finally, we developed a PRIME model to uncover how EGRIN-inferred combinatorial gene regulation by transcription factors, such as CcpA and CodY, modulates essential metabolic processes to enable C. difficile growth relative to commensal colonization. The C. difficile interactive web portal provides access to these model resources to support collaborative systems-level studies of context-specific virulence mechanisms in C. difficile.

In vivo commensal control of Clostridioides difficile virulence.

Leveraging systems biology approaches, we illustrate how metabolically distinct species of Clostridia protect against or worsen Clostridioides difficile infection in mice by modulating the pathogen's colonization, growth, and virulence to impact host survival. Gnotobiotic mice colonized with the amino acid fermenter Paraclostridium bifermentans survive infection with reduced disease severity, while mice colonized with the butyrate-producer, Clostridium sardiniense, succumb more rapidly. Systematic in vivo analyses revealed how each commensal alters the gut-nutrient environment to modulate the pathogen's metabolism, gene regulatory networks, and toxin production. Oral administration of P. bifermentans rescues conventional, clindamycin-treated mice from lethal C. difficile infection in a manner similar to that of monocolonized animals, thereby supporting the therapeutic potential of this commensal species. Our findings lay the foundation for mechanistically informed therapies to counter C. difficile disease using systems biology approaches to define host-commensal-pathogen interactions in vivo.

PerSort Facilitates Characterization and Elimination of Persister Subpopulation in Mycobacteria.

Mycobacterium tuberculosis (MTB) generates phenotypic diversity to persist and survive the harsh conditions encountered during infection. MTB avoids immune effectors and antibacterial killing by entering into distinct physiological states. The surviving cells, persisters, are a major barrier to the timely and relapse-free treatment of tuberculosis (TB). We present for the first time, PerSort, a method to isolate and characterize persisters in the absence of antibiotic or other pressure. We demonstrate the value of PerSort to isolate translationally dormant cells that preexisted in small numbers within Mycobacterium species cultures growing under optimal conditions but that dramatically increased in proportion under stress conditions. The translationally dormant subpopulation exhibited multidrug tolerance and regrowth properties consistent with those of persister cells. Furthermore, PerSort enabled single-cell transcriptional profiling that provided evidence that the translationally dormant persisters were generated through a variety of mechanisms, including vapC30, mazF, and relA/spoT overexpression. Finally, we demonstrate that notwithstanding the varied mechanisms by which the persister cells were generated, they converge on a similar low-oxygen metabolic state that was reversed through activation of respiration to rapidly eliminate persisters fostered under host-relevant stress conditions. We conclude that PerSort provides a new tool to study MTB persisters, enabling targeted strategies to improve and shorten the treatment of TB.IMPORTANCE Mycobacterium tuberculosis (MTB) persists and survives antibiotic treatments by generating phenotypically heterogeneous drug-tolerant subpopulations. The surviving cells, persisters, are a major barrier to the relapse-free treatment of tuberculosis (TB), which is already killing >1.8 million people every year and becoming deadlier with the emergence of multidrug-resistant strains. This study describes PerSort, a cell sorting method to isolate and characterize, without antibiotic treatment, translationally dormant persisters that preexist in small numbers within Mycobacterium cultures. Characterization of this subpopulation has discovered multiple mechanisms by which mycobacterial persisters emerge and unveiled the physiological basis for their dormant and multidrug-tolerant physiological state. This analysis has discovered that activating oxygen respiratory physiology using l-cysteine eliminates preexisting persister subpopulations, potentiating rapid antibiotic killing of mycobacteria under host-relevant stress. PerSort serves as a new tool to study MTB persisters for enabling targeted strategies to improve and shorten the treatment of TB.

Inference of Bacterial Small RNA Regulatory Networks and Integration with Transcription Factor-Driven Regulatory Networks.

Small noncoding RNAs (sRNAs) are key regulators of bacterial gene expression. Through complementary base pairing, sRNAs affect mRNA stability and translation efficiency. Here, we describe a network inference approach designed to identify sRNA-mediated regulation of transcript levels. We use existing transcriptional data sets and prior knowledge to infer sRNA regulons using our network inference tool, the Inferelator This approach produces genome-wide gene regulatory networks that include contributions by both transcription factors and sRNAs. We show the benefits of estimating and incorporating sRNA activities into network inference pipelines using available experimental data. We also demonstrate how these estimated sRNA regulatory activities can be mined to identify the experimental conditions where sRNAs are most active. We uncover 45 novel experimentally supported sRNA-mRNA interactions in Escherichia coli, outperforming previous network-based efforts. Additionally, our pipeline complements sequence-based sRNA-mRNA interaction prediction methods by adding a data-driven filtering step. Finally, we show the general applicability of our approach by identifying 24 novel, experimentally supported, sRNA-mRNA interactions in Pseudomonas aeruginosa, Staphylococcus aureus, and Bacillus subtilis Overall, our strategy generates novel insights into the functional context of sRNA regulation in multiple bacterial species.IMPORTANCE Individual bacterial genomes can have dozens of small noncoding RNAs with largely unexplored regulatory functions. Although bacterial sRNAs influence a wide range of biological processes, including antibiotic resistance and pathogenicity, our current understanding of sRNA-mediated regulation is far from complete. Most of the available information is restricted to a few well-studied bacterial species; and even in those species, only partial sets of sRNA targets have been characterized in detail. To close this information gap, we developed a computational strategy that takes advantage of available transcriptional data and knowledge about validated and putative sRNA-mRNA interactions for inferring expanded sRNA regulons. Our approach facilitates the identification of experimentally supported novel interactions while filtering out false-positive results. Due to its data-driven nature, our method prioritizes biologically relevant interactions among lists of candidate sRNA-target pairs predicted in silico from sequence analysis or derived from sRNA-mRNA binding experiments.

Antibody-recruiting protein-catalyzed capture agents to combat antibiotic-resistant bacteria.

Antibiotic resistant infections are projected to cause over 10 million deaths by 2050, yet the development of new antibiotics has slowed. This points to an urgent need for methodologies for the rapid development of antibiotics against emerging drug resistant pathogens. We report on a generalizable combined computational and synthetic approach, called antibody-recruiting protein-catalyzed capture agents (AR-PCCs), to address this challenge. We applied the combinatorial protein catalyzed capture agent (PCC) technology to identify macrocyclic peptide ligands against highly conserved surface protein epitopes of carbapenem-resistant Klebsiella pneumoniae, an opportunistic Gram-negative pathogen with drug resistant strains. Multi-omic data combined with bioinformatic analyses identified epitopes of the highly expressed MrkA surface protein of K. pneumoniae for targeting in PCC screens. The top-performing ligand exhibited high-affinity (EC50 ∼50 nM) to full-length MrkA, and selectively bound to MrkA-expressing K. pneumoniae, but not to other pathogenic bacterial species. AR-PCCs that bear a hapten moiety promoted antibody recruitment to K. pneumoniae, leading to enhanced phagocytosis and phagocytic killing by macrophages. The rapid development of this highly targeted antibiotic implies that the integrated computational and synthetic toolkit described here can be used for the accelerated production of antibiotics against drug resistant bacteria.