Dysregulated innate immune signaling cooperates with RUNX1 mutations to transform an MDS-like disease to AML.

Dysregulated innate immune signaling is linked to preleukemic conditions and myeloid malignancies. However, it is unknown whether sustained innate immune signaling contributes to malignant transformation. Here we show that cell-intrinsic innate immune signaling driven by miR-146a deletion (miR-146aKO), a commonly deleted gene in myelodysplastic syndromes (MDS) and acute myeloid leukemia (AML), cooperates with mutant RUNX1 (RUNX1mut) to initially induce marrow failure and features of MDS. However, miR-146aKO hematopoietic stem and/or progenitor cells (HSPCs) expressing RUNX1mut eventually progress to a fatal AML. miR-146aKO HSPCs exhaust during serial transplantation, while expression of RUNX1mut restored their hematopoietic cell function. Thus, HSPCs exhibiting dysregulated innate immune signaling require a second hit to develop AML. Inhibiting the dysregulated innate immune pathways with a TRAF6-UBE2N inhibitor suppressed leukemic miR-146aKO/RUNX1mut HSPCs, highlighting the necessity of TRAF6-dependent cell-intrinsic innate immune signaling in initiating and maintaining AML. These findings underscore the critical role of dysregulated cell-intrinsic innate immune signaling in driving preleukemic cells toward AML progression.

VExD: a curated resource for human gene expression alterations following viral infection.

Much of the host antiviral response is mediated through changes to host gene expression levels. Likewise, viruses induce changes to host gene expression levels in order to promote the viral life cycle and evade the host immune system. However, there is no resource that specifically collects human gene expression levels pre- and post-virus infection. Further, public gene expression repositories do not contain enough specialized metadata to easily find relevant experiments. Here, we present the Virus Expression Database (VExD), a freely available website and database, that collects human gene expression datasets in response to viral infection. VExD contains ∼8,000 uniformly processed samples obtained from 289 studies examining 51 distinct human viruses. We show that the VExD processing pipeline captures known antiviral responses in the form of interferon-stimulated genes. We further show that the datasets collected in VExD can be used to quickly identify supporting data for experiments performed in human cells or model organisms. VExD is freely available at https://vexd.cchmc.org/.

PU.1-Dependent Enhancer Inhibition Separates Tet2-Deficient Hematopoiesis from Malignant Transformation.

Cytosine hypermethylation in and around DNA-binding sites of master transcription factors, including PU.1, occurs in aging hematopoietic stem cells following acquired loss-of-function mutations of DNA methyl-cytosine dioxygenase ten-eleven translocation-2 (TET2), albeit functional relevance has been unclear. We show that Tet2-deficient mouse hematopoietic stem and progenitor cells undergo malignant transformation upon compromised gene regulation through heterozygous deletion of an upstream regulatory region (UREΔ/WT) of the PU.1 gene. Although compatible with multilineage blood formation at young age, Tet2-deficient PU.1 UREΔ/WT mice develop highly penetrant, transplantable acute myeloid leukemia (AML) during aging. Leukemic stem and progenitor cells show hypermethylation at putative PU.1-binding sites, fail to activate myeloid enhancers, and are hallmarked by a signature of genes with impaired expression shared with human AML. Our study demonstrates that Tet2 and PU.1 jointly suppress leukemogenesis and uncovers a methylation-sensitive PU.1-dependent gene network as a unifying molecular vulnerability associated with AML. SIGNIFICANCE: We identify moderately impaired PU.1 mRNA expression as a biological modality predisposing Tet2-deficient hematopoietic stem and progenitor cells to malignant transformation. Our study furthermore uncovers a methylation-sensitive PU.1 gene network as a common feature of myeloid leukemia potentially allowing for the identification of patients at risk for malignant transformation. See related commentary by Schleicher and Pietras, p. 378. This article is highlighted in the In This Issue feature, p. 369.

Multi-omic characterization of pediatric ARDS via nasal brushings.

RATIONALE: While nasal brushing transcriptomics can identify disease subtypes in chronic pulmonary diseases, it is unknown whether this is true in pediatric acute respiratory distress syndrome (PARDS). OBJECTIVES: Determine whether nasal transcriptomics and methylomics can identify clinically meaningful PARDS subgroups that reflect important pathobiological processes. METHODS: Nasal brushings and serum were collected on days 1, 3, 7, and 14 from control and PARDS subjects from two centers. PARDS duration was the primary endpoint. MEASUREMENTS AND MAIN RESULTS: Twenty-four control and 39 PARDS subjects were enrolled. Two nasal methylation patterns were identified. Compared to Methyl Subgroup 1, Subgroup 2 had hypomethylation of inflammatory genes and was enriched for immunocompromised subjects. Four transcriptomic patterns were identified with temporal patterns indicating injury, repair, and regeneration. Over time, both inflammatory (Subgroup B) and cell injury (Subgroup D) patterns transitioned to repair (Subgroup A) and eventually homeostasis (Subgroup C). When control specimens were included, they were largely Subgroup C. In comparison with 17 serum biomarkers, the nasal transcriptome was more predictive of prolonged PARDS. Subjects with initial Transcriptomic Subgroup B or D assignment had median PARDS duration of 8 days compared to 2 in A or C (p = 0.02). For predicting PARDS duration ≥ 3 days, nasal transcriptomics was more sensitive and serum biomarkers more specific. CONCLUSIONS: PARDS nasal transcriptome may reflect distal lung injury, repair, and regeneration. A combined nasal PCR and serum biomarker assay could be useful for predictive and diagnostic enrichment. Trial registration Clinicaltrials.gov NCT03539783 May 29, 2018.

Epigenetic and transcriptional dysregulation in CD4+ T cells in patients with atopic dermatitis.

Atopic dermatitis (AD) is one of the most common skin disorders among children. Disease etiology involves genetic and environmental factors, with 29 independent AD risk loci enriched for risk allele-dependent gene expression in the skin and CD4+ T cell compartments. We investigated the potential epigenetic mechanisms responsible for the genetic susceptibility of CD4+ T cells. To understand the differences in gene regulatory activity in peripheral blood T cells in AD, we measured chromatin accessibility (an assay based on transposase-accessible chromatin sequencing, ATAC-seq), nuclear factor kappa B subunit 1 (NFKB1) binding (chromatin immunoprecipitation with sequencing, ChIP-seq), and gene expression levels (RNA-seq) in stimulated CD4+ T cells from subjects with active moderate-to-severe AD, as well as in age-matched non-allergic controls. Open chromatin regions in stimulated CD4+ T cells were highly enriched for AD genetic risk variants, with almost half of the AD risk loci overlapping AD-dependent ATAC-seq peaks. AD-specific open chromatin regions were strongly enriched for NF-κB DNA-binding motifs. ChIP-seq identified hundreds of NFKB1-occupied genomic loci that were AD- or control-specific. As expected, the AD-specific ChIP-seq peaks were strongly enriched for NF-κB DNA-binding motifs. Surprisingly, control-specific NFKB1 ChIP-seq peaks were not enriched for NFKB1 motifs, but instead contained motifs for other classes of human transcription factors, suggesting a mechanism involving altered indirect NFKB1 binding. Using DNA sequencing data, we identified 63 instances of altered genotype-dependent chromatin accessibility at 36 AD risk variant loci (30% of AD risk loci) that might lead to genotype-dependent gene expression. Based on these findings, we propose that CD4+ T cells respond to stimulation in an AD-specific manner, resulting in disease- and genotype-dependent chromatin accessibility alterations involving NFKB1 binding.

Blocking UBE2N abrogates oncogenic immune signaling in acute myeloid leukemia.

Dysregulation of innate immune signaling pathways is implicated in various hematologic malignancies. However, these pathways have not been systematically examined in acute myeloid leukemia (AML). We report that AML hematopoietic stem and progenitor cells (HSPCs) exhibit a high frequency of dysregulated innate immune-related and inflammatory pathways, referred to as oncogenic immune signaling states. Through gene expression analyses and functional studies in human AML cell lines and patient-derived samples, we found that the ubiquitin-conjugating enzyme UBE2N is required for leukemic cell function in vitro and in vivo by maintaining oncogenic immune signaling states. It is known that the enzyme function of UBE2N can be inhibited by interfering with thioester formation between ubiquitin and the active site. We performed in silico structure-based and cellular-based screens and identified two related small-molecule inhibitors UC-764864/65 that targeted UBE2N at its active site. Using these small-molecule inhibitors as chemical probes, we further revealed the therapeutic efficacy of interfering with UBE2N function. This resulted in the blocking of ubiquitination of innate immune- and inflammatory-related substrates in human AML cell lines. Inhibition of UBE2N function disrupted oncogenic immune signaling by promoting cell death of leukemic HSPCs while sparing normal HSPCs in vitro. Moreover, baseline oncogenic immune signaling states in leukemic cells derived from discrete subsets of patients with AML exhibited a selective dependency on UBE2N function in vitro and in vivo. Our study reveals that interfering with UBE2N abrogates leukemic HSPC function and underscores the dependency of AML cells on UBE2N-dependent oncogenic immune signaling states.

Runx1 shapes the chromatin landscape via a cascade of direct and indirect targets.

Runt-related transcription factor 1 (Runx1) can act as both an activator and a repressor. Here we show that CRISPR-mediated deletion of Runx1 in mouse metanephric mesenchyme-derived mK4 cells results in large-scale genome-wide changes to chromatin accessibility and gene expression. Open chromatin regions near down-regulated loci enriched for Runx sites in mK4 cells lose chromatin accessibility in Runx1 knockout cells, despite remaining Runx2-bound. Unexpectedly, regions near upregulated genes are depleted of Runx sites and are instead enriched for Zeb transcription factor binding sites. Re-expressing Zeb2 in Runx1 knockout cells restores suppression, and CRISPR mediated deletion of Zeb1 and Zeb2 phenocopies the gained expression and chromatin accessibility changes seen in Runx1KO due in part to subsequent activation of factors like Grhl2. These data confirm that Runx1 activity is uniquely needed to maintain open chromatin at many loci, and demonstrate that Zeb proteins are required and sufficient to maintain Runx1-dependent genome-scale repression.

Global discovery of lupus genetic risk variant allelic enhancer activity.

Genome-wide association studies of Systemic Lupus Erythematosus (SLE) nominate 3073 genetic variants at 91 risk loci. To systematically screen these variants for allelic transcriptional enhancer activity, we construct a massively parallel reporter assay (MPRA) library comprising 12,396 DNA oligonucleotides containing the genomic context around every allele of each SLE variant. Transfection into the Epstein-Barr virus-transformed B cell line GM12878 reveals 482 variants with enhancer activity, with 51 variants showing genotype-dependent (allelic) enhancer activity at 27 risk loci. Comparison of MPRA results in GM12878 and Jurkat T cell lines highlights shared and unique allelic transcriptional regulatory mechanisms at SLE risk loci. In-depth analysis of allelic transcription factor (TF) binding at and around allelic variants identifies one class of TFs whose DNA-binding motif tends to be directly altered by the risk variant and a second class of TFs that bind allelically without direct alteration of their motif by the variant. Collectively, our approach provides a blueprint for the discovery of allelic gene regulation at risk loci for any disease and offers insight into the transcriptional regulatory mechanisms underlying SLE.

Comparative analysis of human microglial models for studies of HIV replication and pathogenesis.

BACKGROUND: HIV associated neurocognitive disorders cause significant morbidity and mortality despite the advent of highly active antiretroviral therapy. A deeper understanding of fundamental mechanisms underlying HIV infection and pathogenesis in the central nervous system is warranted. Microglia are resident myeloid cells of the brain that are readily infected by HIV and may constitute a CNS reservoir. We evaluated two microglial model cell lines (C20, HMC3) and two sources of primary cell-derived microglia (monocyte-derived microglia [MMG] and induced pluripotent stem cell-derived microglia [iPSC-MG]) as potential model systems for studying HIV-microglia interactions. RESULTS: All four microglial model cells expressed typical myeloid markers with the exception of low or absent CD45 and CD11b expression by C20 and HMC3, and all four expressed the microglia-specific markers P2RY12 and TMEM119. Marked differences were observed upon gene expression profiling, however, indicating that MMG and iPSC-MG cluster closely together with primary human microglial cells, while C20 and HMC3 were similar to each other but very different from primary microglia. Expression of HIV-relevant genes also revealed important differences, with iPSC-MG and MMG expressing relevant genes at levels more closely resembling primary microglia. iPSC-MG and MMG were readily infected with R5-tropic HIV, while C20 and HMC3 lack CD4 and require pseudotyping for infection. Despite many similarities, HIV replication dynamics and HIV-1 particle capture by Siglec-1 differed markedly between the MMG and iPSC-MG. CONCLUSIONS: MMG and iPSC-MG appear to be viable microglial models that are susceptible to HIV infection and bear more similarities to authentic microglia than two transformed microglia cell lines. The observed differences in HIV replication and particle capture between MMG and iPSC-MG warrant further study.

TIFAB Regulates USP15-Mediated p53 Signaling during Stressed and Malignant Hematopoiesis.

TRAF-interacting protein with a forkhead-associated domain B (TIFAB) is implicated in myeloid malignancies with deletion of chromosome 5q. Employing a combination of proteomic and genetic approaches, we find that TIFAB regulates ubiquitin-specific peptidase 15 (USP15) ubiquitin hydrolase activity. Expression of TIFAB in hematopoietic stem/progenitor cells (HSPCs) permits USP15 signaling to substrates, including MDM2 and KEAP1, and mitigates p53 expression. Consequently, TIFAB-deficient HSPCs exhibit compromised USP15 signaling and are sensitized to hematopoietic stress by derepression of p53. In MLL-AF9 leukemia, deletion of TIFAB increases p53 signaling and correspondingly decreases leukemic cell function and development of leukemia. Restoring USP15 expression partially rescues the function of TIFAB-deficient MLL-AF9 cells. Conversely, elevated TIFAB represses p53, increases leukemic progenitor function, and correlates with MLL gene expression programs in leukemia patients. Our studies uncover a function of TIFAB as an effector of USP15 activity and rheostat of p53 signaling in stressed and malignant HSPCs.

Genomic programming of IRF4-expressing human Langerhans cells.

Langerhans cells (LC) can prime tolerogenic as well as immunogenic responses in skin, but the genomic states and transcription factors (TF) regulating these context-specific responses are unclear. Bulk and single-cell transcriptional profiling demonstrates that human migratory LCs are robustly programmed for MHC-I and MHC-II antigen presentation. Chromatin analysis reveals enrichment of ETS-IRF and AP1-IRF composite regulatory elements in antigen-presentation genes, coinciding with expression of the TFs, PU.1, IRF4 and BATF3 but not IRF8. Migration of LCs from the epidermis is accompanied by upregulation of IRF4, antigen processing components and co-stimulatory molecules. TNF stimulation augments LC cross-presentation while attenuating IRF4 expression. CRISPR-mediated editing reveals IRF4 to positively regulate the LC activation programme, but repress NF2EL2 and NF-kB pathway genes that promote responsiveness to oxidative stress and inflammatory cytokines. Thus, IRF4-dependent genomic programming of human migratory LCs appears to enable LC maturation while attenuating excessive inflammatory and immunogenic responses in the epidermis.

CCL2 mobilizes ALIX to facilitate Gag-p6 mediated HIV-1 virion release.

Cellular ESCRT machinery plays pivotal role in HIV-1 budding and release. Extracellular stimuli that modulate HIV-1 egress are currently unknown. We found that CCL2 induced by HIV-1 clade B (HIV-1B) infection of macrophages enhanced virus production, while CCL2 immuno-depletion reversed this effect. Additionally, HIV-1 clade C (HIV-1C) was refractory to CCL2 levels. We show that CCL2-mediated increase in virus production requires Gag late motif LYPX present in HIV-1B, but absent in HIV-1C, and ALIX protein that recruits ESCRT III complex. CCL2 immuno-depletion sequestered ALIX to F-actin structures, while CCL2 addition mobilized it to cytoplasm facilitating Gag-ALIX binding. The LYPX motif improves virus replication and its absence renders the virus less fit. Interestingly, novel variants of HIV-1C with PYRE/PYKE tetrapeptide insertions in Gag-p6 conferred ALIX binding, CCL2-responsiveness and enhanced virus replication. These results, for the first time, indicate that CCL2 mediates ALIX mobilization from F-actin and enhances HIV-1 release and fitness.

A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus.

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.

Transcription factors operate across disease loci, with EBNA2 implicated in autoimmunity.

Explaining the genetics of many diseases is challenging because most associations localize to incompletely characterized regulatory regions. Using new computational methods, we show that transcription factors (TFs) occupy multiple loci associated with individual complex genetic disorders. Application to 213 phenotypes and 1,544 TF binding datasets identified 2,264 relationships between hundreds of TFs and 94 phenotypes, including androgen receptor in prostate cancer and GATA3 in breast cancer. Strikingly, nearly half of systemic lupus erythematosus risk loci are occupied by the Epstein-Barr virus EBNA2 protein and many coclustering human TFs, showing gene-environment interaction. Similar EBNA2-anchored associations exist in multiple sclerosis, rheumatoid arthritis, inflammatory bowel disease, type 1 diabetes, juvenile idiopathic arthritis and celiac disease. Instances of allele-dependent DNA binding and downstream effects on gene expression at plausibly causal variants support genetic mechanisms dependent on EBNA2. Our results nominate mechanisms that operate across risk loci within disease phenotypes, suggesting new models for disease origins.

Control of mitochondrial function and cell growth by the atypical cadherin Fat1.

Mitochondrial products such as ATP, reactive oxygen species, and aspartate are key regulators of cellular metabolism and growth. Abnormal mitochondrial function compromises integrated growth-related processes such as development and tissue repair, as well as homeostatic mechanisms that counteract ageing and neurodegeneration, cardiovascular disease, and cancer. Physiologic mechanisms that control mitochondrial activity in such settings remain incompletely understood. Here we show that the atypical Fat1 cadherin acts as a molecular 'brake' on mitochondrial respiration that regulates vascular smooth muscle cell (SMC) proliferation after arterial injury. Fragments of Fat1 accumulate in SMC mitochondria, and the Fat1 intracellular domain interacts with multiple mitochondrial proteins, including critical factors associated with the inner mitochondrial membrane. SMCs lacking Fat1 (Fat1KO) grow faster, consume more oxygen for ATP production, and contain more aspartate. Notably, expression in Fat1KO cells of a modified Fat1 intracellular domain that localizes exclusively to mitochondria largely normalizes oxygen consumption, and the growth advantage of these cells can be suppressed by inhibition of mitochondrial respiration, which suggest that a Fat1-mediated growth control mechanism is intrinsic to mitochondria. Consistent with this idea, Fat1 species associate with multiple respiratory complexes, and Fat1 deletion both increases the activity of complexes I and II and promotes the formation of complex-I-containing supercomplexes. In vivo, Fat1 is expressed in injured human and mouse arteries, and inactivation of SMC Fat1 in mice potentiates the response to vascular damage, with markedly increased medial hyperplasia and neointimal growth, and evidence of higher SMC mitochondrial respiration. These studies suggest that Fat1 controls mitochondrial activity to restrain cell growth during the reparative, proliferative state induced by vascular injury. Given recent reports linking Fat1 to cancer, abnormal kidney and muscle development, and neuropsychiatric disease, this Fat1 function may have importance in other settings of altered cell growth and metabolism.

Iterative sorting of apical and basolateral cargo in Madin-Darby canine kidney cells.

For several decades, the trans-Golgi network (TGN) was considered the most distal stop and hence the ultimate protein-sorting station for distinct apical and basolateral transport carriers that reach their respective surface domains in the direct trafficking pathway. However, recent reports of apical and basolateral cargoes traversing post-Golgi compartments accessible to endocytic ligands before their arrival at the cell surface and the post-TGN breakup of large pleomorphic membrane fragments that exit the Golgi region toward the surface raised the possibility that compartments distal to the TGN mediate or contribute to biosynthetic sorting. Here we describe the development of a novel assay that quantitatively distinguishes different cargo pairs by their degree of colocalization at the TGN and by the evolution of colocalization during their TGN-to-surface transport. Keys to the high resolution of our approach are 1) conversion of perinuclear organelle clustering into a two-dimensional microsomal spread and 2) identification of TGN and post-TGN cargo without the need for a TGN marker that universally cosegregates with all cargo. Using our assay, we provide the first evidence that apical NTRp75 and basolateral VSVG in Madin-Darby canine kidney cells still undergo progressive sorting after they exit the TGN toward the cell surface.

Prediction of DNA binding motifs from 3D models of transcription factors; identifying TLX3 regulated genes.

Proper cell functioning depends on the precise spatio-temporal expression of its genetic material. Gene expression is controlled to a great extent by sequence-specific transcription factors (TFs). Our current knowledge on where and how TFs bind and associate to regulate gene expression is incomplete. A structure-based computational algorithm (TF2DNA) is developed to identify binding specificities of TFs. The method constructs homology models of TFs bound to DNA and assesses the relative binding affinity for all possible DNA sequences using a knowledge-based potential, after optimization in a molecular mechanics force field. TF2DNA predictions were benchmarked against experimentally determined binding motifs. Success rates range from 45% to 81% and primarily depend on the sequence identity of aligned target sequences and template structures, TF2DNA was used to predict 1321 motifs for 1825 putative human TF proteins, facilitating the reconstruction of most of the human gene regulatory network. As an illustration, the predicted DNA binding site for the poorly characterized T-cell leukemia homeobox 3 (TLX3) TF was confirmed with gel shift assay experiments. TLX3 motif searches in human promoter regions identified a group of genes enriched in functions relating to hematopoiesis, tissue morphology, endocrine system and connective tissue development and function.

The underlying molecular and network level mechanisms in the evolution of robustness in gene regulatory networks.

Gene regulatory networks show robustness to perturbations. Previous works identified robustness as an emergent property of gene network evolution but the underlying molecular mechanisms are poorly understood. We used a multi-tier modeling approach that integrates molecular sequence and structure information with network architecture and population dynamics. Structural models of transcription factor-DNA complexes are used to estimate relative binding specificities. In this model, mutations in the DNA cause changes on two levels: (a) at the sequence level in individual binding sites (modulating binding specificity), and (b) at the network level (creating and destroying binding sites). We used this model to dissect the underlying mechanisms responsible for the evolution of robustness in gene regulatory networks. Results suggest that in sparse architectures (represented by short promoters), a mixture of local-sequence and network-architecture level changes are exploited. At the local-sequence level, robustness evolves by decreasing the probabilities of both the destruction of existent and generation of new binding sites. Meanwhile, in highly interconnected architectures (represented by long promoters), robustness evolves almost entirely via network level changes, deleting and creating binding sites that modify the network architecture.

Identification of linear epitopes in Bacillus anthracis protective antigen bound by neutralizing antibodies.

Protective antigen (PA), the binding subunit of anthrax toxin, is the major component in the current anthrax vaccine, but the fine antigenic structure of PA is not well defined. To identify linear neutralizing epitopes of PA, 145 overlapping peptides covering the entire sequence of the protein were synthesized. Six monoclonal antibodies (mAbs) and antisera from mice specific for PA were tested for their reactivity to the peptides by enzyme-linked immunosorbent assays. Three major linear immunodominant B-cell epitopes were mapped to residues Leu(156) to Ser(170), Val(196) to Ile(210), and Ser(312) to Asn(326) of the PA protein. Two mAbs with toxin-neutralizing activity recognized two different epitopes in close proximity to the furin cleavage site in domain 1. The three-dimensional complex structure of PA and its neutralizing mAbs 7.5G and 19D9 were modeled using the molecular docking method providing models for the interacting epitope and paratope residues. For both mAbs, LeTx neutralization was associated with interference with furin cleavage, but they differed in effectiveness depending on whether they bound on the N- or C-terminal aspect of the cleaved products. The two peptides containing these epitopes that include amino acids Leu(156)-Ser(170) and Val(196)-Ile(210) were immunogenic and elicited neutralizing antibody responses to PA. These results identify the first linear neutralizing epitopes of PA and show that peptides containing epitope sequences can elicit neutralizing antibody responses, a finding that could be exploited for vaccine design.