Abstract

Pseudomonas aeruginosa is the most common pathogen that persists in the cystic fibrosis lungs. Bacteria such as P. aeruginosa secrete siderophores (iron-chelating molecules) and the host limits bacterial growth by producing neutrophil-gelatinase-associated lipocalin (NGAL) that specifically scavenges bacterial siderophores, therefore preventing bacteria from establishing infection. P. aeruginosa produces a major siderophore known as pyoverdine, found to be important for bacterial virulence and biofilm development. We report that pyoverdine did not bind to NGAL, as measured by tryptophan fluorescence quenching, while enterobactin bound to NGAL effectively causing a strong response. The experimental data indicate that pyoverdine evades NGAL recognition. We then employed a molecular modeling approach to simulate the binding of pyoverdine to human NGAL using NGAL’s published crystal structures. The docking of pyoverdine to NGAL predicted nine different docking positions; however, neither apo- nor ferric forms of pyoverdine docked into the ligand-binding site in the calyx of NGAL where siderophores are known to bind. The molecular modeling results offer structural support that pyoverdine does not bind to NGAL, confirming the results obtained in the tryptophan quenching assay. The data suggest that pyoverdine is a stealth siderophore that evades NGAL recognition allowing P. aeruginosa to establish chronic infections in CF lungs.

1. Introduction

Cystic fibrosis (CF) is an inherited disease resulting in the formation of thick, sticky mucus in the lungs and digestive pathways [1]. CF lungs are breeding grounds for a variety of bacterial pathogens. Pseudomonas aeruginosa is the most common bacterial isolate that persists in CF lungs by evading host defenses and clearance mechanisms, contributing to declining lung function [1, 2]. In order to grow and thrive, bacteria need iron as a cofactor for many metabolic enzymes. However, in the case of infection, iron is not freely available in the host and is predominantly found in iron-binding proteins such as transferrin, lactoferrin, and ferritin [3]. Thus, due to the low bioavailability of free iron, bacteria produce iron-chelating small molecules known as siderophores that are secreted in apo-form and imported as ferric complexes [1, 4]. P. aeruginosa produces two siderophores, pyoverdine (Figures 1(a), and 1(b)) and pyochelin, that are found to be important for biofilm development and bacterial virulence [4].

As a host defense strategy, mammals have evolved strategies to limit bacterial growth by producing antibacterial iron-depleting innate immune defense proteins, such as neutrophil-gelatinase-associated lipocalin (NGAL), that specifically scavenge bacterial ferric and apo-siderophores, therefore preventing bacteria from establishing infection [5, 6]. NGAL is a 25 kDa protein also known as siderocalin, lipocalin 2, uterocalin, and 24p3, that exhibits antibacterial activity [6]. NGAL expression is upregulated following bacterial infection in an NFκB-TLR4-dependent manner [7]. NGAL-deficient mice are more susceptible to bacterial infections [8] and NGAL is required for pulmonary host defense against Klebsiella pneumoniae [7]. Lipocalins such as NGAL share a common tertiary structure, the so-called “lipocalin fold,” characterized by an eight-strand antiparallel 𝛽 -barrel with a calyx that contains the ligand-binding site. NGAL is expressed and secreted upon stimulation by many cells and also acts as a growth and differentiation factor in various tissues, mainly in the kidney [9]. One study has shown that NGAL levels are highly elevated in CF serum and bronchoalveolar lavage fluid (BALF) [10]. In spite of the abundance of NGAL in serum and BALF of CF patients, P. aeruginosa infections persist for many years.

Pathogenic bacteria evolved under the pressure of iron limitation to survive and evade host immune systems by producing structurally modified siderophores that do not bind to NGAL [11]. Pathogenic Gram-negative bacteria harbor the iroA gene cluster that mediates bacterial evasion of NGAL [12, 13]. For example, Salmonella typhimurium, a Gram-negative bacterium, produces a glycosylated enterobactin known as salmochelin [11, 14]. No homologs of the iroA gene cluster were found in the published Pseudomonas genomes, although P. aeruginosa produces three types of pyoverdines [1]. Other pathogenic bacteria also produce structurally distinct siderophores that evade NGAL recognition. For example, Bacillus anthracis, a Gram-positive bacterium, produces a siderophore known as petrobactin that does not bind to NGAL [15]; Klebsiella pneumonia, a Gram-negative bacterium, produces a siderophore known as yersiniabactin [16] that also does not bind to NGAL and thereby evades the host’s antibacterial iron-depleting defense. Further, Vibrio cholerae produces vibriobactin, a catecholate siderophore that evades NGAL recognition by unique iron coordination in the ferric vibriobactin form [17].

It is not known if pyoverdine, the major siderophore produced by P. aeruginosa, evades this innate defense mechanism and whether this contributes to inability of CF patients to clear this bacterium from their lungs. Therefore, the aim of this work was to determine if pyoverdine binds to the antibacterial protein NGAL. We hypothesized that NGAL does not bind pyoverdine, allowing P. aeruginosa to establish infections and persist in CF lungs. Hence, this iron-depleting innate host defense mechanism may be dysfunctional and impaired in CF. In this paper we employed a fluorescence-based assay to characterize the binding of pyoverdine to human NGAL and molecular modeling of NGAL pyoverdine complexes to predict possible protein-ligand interactions. Here we report that neither apo- nor ferric pyoverdine bind to NGAL, results that are further supported by computational predictions indicating lack of interaction of pyoverdine with NGAL at its conventional ligand-binding site. We conclude that pyoverdine is a stealth siderophore that evades NGAL recognition, thereby allowing P. aeruginosa to establish infections.

2. Materials and Methods

2.1. Reagents

Recombinant human NGAL (rhNGAL, catalog no. 1757-LC) produced in a mammalian system rather than a bacterial source was purchased from R&D Systems (Minneapolis, MN, USA). Highly purified and structurally defined pyoverdines were obtained from Sigma-Aldrich (St Louis, MO, USA) in apo-form, that is, aferric pyoverdine (catalog no. P8124), and iron-bound ferric pyoverdine (catalog no. P8374). Enterobactin (catalog no. E3910), iron standard (catalog no. 16596), 2,3-dihydroxybenzoic acid (2,3-DHBA, catalog no. 0926 KB), and sterile bovine serum albumin (BSA, endotoxin free tissue culture grade, CAS 9048468) were purchased from Sigma-Aldrich.

2.2. Tryptophan Fluorescence Quenching

Fluorescence measurements were made using a SpectraMax microtiter plate reader by Molecular Devices. Human NGAL contains two tryptophan residues: W31 located at the base of the hydrophobic calyx and W79 located inside the calyx (i.e., the ligand-binding site). If pyoverdine binds to NGAL it will cause quenching of the intrinsic tryptophan fluorescence. Therefore, we employed an established tryptophan fluorescence-quenching assay [13, 15] to determine the binding of pyoverdine to NGAL. Since 2,3-DHBA is the core functional moiety in pyoverdine that binds iron, the fluorescence assay was first optimized by testing ferric iron and 2,3-DHBA complexes diluted in 0.1 M Tris, pH 8.0, as the ligand for potentially quenching intrinsic tryptophan fluorescence in rhNGAL and in sterile BSA. Enterobactin (ENT), the major siderophore from E. coli and other enteric bacteria, was used as a positive control. The enterobactin stock solution was dissolved in DMSO, whereas apo-pyoverdine (PVD) and ferric pyoverdine (PVD-Fe) stocks were dissolved in sterile water; then, serial dilutions of the siderophores at 100, 50, 25, and 12.5  𝜇 M were made in an assay buffer (50 mM Tris, 10 mM CaCl2, 150 mM NaCl pH 7.5). Fifty 𝜇 L of each dilution was then transferred into triplicate wells (F16 Black Maxisorp Plates from Nunc). Fifty 𝜇 L of rhNGAL was added into each well at a final concentration of 2  𝜇 M and incubated for 30 min at room temperature. Recombinant hNGAL alone was incubated with assay buffer without any siderophore and used to measure intrinsic tryptophan fluorescence. The plate was then read at excitation and emission wavelengths of 280 and 340 nm, respectively, in the endpoint mode using a microplate reader with fluorescence detection. Results are presented as the percent of tryptophan fluorescence remaining compared to rhNGAL alone.

2.3. Molecular Modeling
2.3.1. Protein File Selection

Among 18 published crystal structures of NGAL, we selected two human NGAL structures solved by Strong and coworkers [18, 19] to serve as the basis for the molecular modeling study files 1L6M and 1QQS from the Protein Data Bank (PDB) [20]. For simplicity, we refer to the 1L6M structure as the NGAL “trimer” and the 1QQS structure as the NGAL “monomer” in this paper. Both structures were solved to 2.4  Å resolution, the highest resolution solved to date for an NGAL structure bound to enterobactin, and in the same space group (P41212). The two structures are virtually identical in terms of conformation even though differences were observed in crystal packing in spite of matching space groups. Three monomers were observed in the 1L6M structure while the 1QQS structure was a homodimer with only one monomer. Moreover, NGAL has been found in a 25 kD monomeric and 45 kD dimeric biological forms, as well as a 135 kD heterodimer with MMP9 [2123]. Due to the numerous structural similarities, minor distinctions in crystal packing, and the availability of different quaternary forms of NGAL, the 1L6M and 1QQS structures were selected as interesting models for predicting different types of interactions that pyoverdine might be able to make with NGAL and for proposing structure-activity relationships.

Another important structural feature relevant to this study is the presence and location of tryptophan residues in the NGAL structures. Biologically, human NGAL has two tryptophan residues W31 and W79. However, coordinates for the NGAL monomer (1QQS) contained a W79A where the side chain of W79 was disordered in the structure and modeled as an alanine. Otherwise, all other tryptophans present in the crystal structures were used for molecular modeling.

2.3.2. Protein File Preparation

Coordinates for the NGAL trimer were extracted from PDB file 1L6M. All heteroatoms including water molecules and definitions for symmetry were excluded from the file. Coordinates for the monomeric form of NGAL were extracted from PDB file 1QQS. All water atoms were excluded from the file. Covalent attachment of the sugars N-acetyl-D-glucosamine and α-D-mannose and the lipid decanoic acid to the monomer protein were retained. Polar hydrogens were added to each protein and each set of modified PDB coordinates was converted to PDBQT format using AutoDock Tools [24]. PDBQT format includes information about partial atomic charges and atom types in addition to the otherwise conventional PDB format.

2.3.3. Ligand File Preparation

Coordinates for ferric pyoverdine (PVD-Fe) were extracted from PDB file 2W78, which describes a complex between iron-bound PVD-ATCC13535 and its transporter protein from P. aeruginosa. The specific form of pyoverdine used has a chemical formula of C17H17N3O6. All other heteroatoms including the sulfate counterion were excluded. The iron atom was excluded from the working ferric pyoverdine file to generate the analogous coordinates for pyoverdine without iron bound (i.e., apo-pyoverdine). Coordinates for the ferric enterobactin control were extracted from PDB file 3CMP. Other heteroatoms present in the 3CMP file including the sodium and sulfate counterions were excluded. Connectivity definitions between atoms for each ligand file were included. Modified PDB files were converted to PDBQT format using AutoDock Tools.

2.3.4. Docking

AutoDock Vina [24] was used to dock the ligand into the receptor protein. The protein was kept rigid in all experiments. Prior to the docking experiment, a grid box was selected in AutoDock Tools to encompass all atoms within either the modified NGAL monomer or trimer structure. The search space for AutoDock Vina was defined to be the dimensions of the 3D grid box. Output files were visualized using Python Molecular Viewer (PMV 1.5.4; [27]).

3. Results

3.1. Pyoverdine Does Not Bind to NGAL

NGAL is an antibacterial innate defense protein that limits bacterial growth by binding to and sequestering siderophores [6]. We hypothesized that NGAL does not bind pyoverdine, allowing P. aeruginosa to persist in CF lungs. To test this hypothesis, we measured the direct binding of pyoverdine (Figure 1) to rhNGAL using the tryptophan fluorescence-quenching method [13]. Human NGAL contains two tryptophan residues: W31, which is located at the base of the hydrophobic calyx, and W79, located at the ligand-binding site. If pyoverdine binds to NGAL it would be expected to cause the quenching of the intrinsic fluorescence of tryptophan.

The quenching of the intrinsic tryptophan fluorescence of rhNGAL in a dose-dependent manner was first optimized using iron complexed to 2,3-dihydroxybenzoic acid (Fe-DHBA) (Figure 2(a)). As expected, increasing concentrations of Fe-DHBA led to increased quenching, suggesting binding of the siderophore core iron-binding structure to NGAL. As another positive control, enterobactin induced a response similar to that of Fe-DHBA. Ferric enterobactin-NGAL complexes showed a fluorescence decrease as the siderophore concentration increased, indicating fluorescence quenching and binding of the ligand to NGAL (Figure 2(b)). In contrast to enterobactin, neither apo-pyoverdine nor ferric pyoverdine reduced the intrinsic tryptophan fluorescence of rhNGAL significantly (Figure 2(b)). Tryptophan fluorescence was only diminished by less than 10% of the NGAL baseline fluorescence even when a high dose of pyoverdine was used, suggesting that pyoverdine does not bind to NGAL even when complexed with iron which is known to increase the affinity of siderophore binding to NGAL [12]. The data suggest that NGAL does not bind pyoverdine. Thus, NGAL does not exert an antibacterial function against P. aeruginosa and, therefore, pyoverdine is a stealth siderophore.

3.2. Pyoverdine-NGAL Interactions Predicted from Molecular Modeling
3.2.1. Apo-Pyoverdine Docking to NGAL

Docking of apo-pyoverdine (Figure 1(a)) to the NGAL trimer predicted nine different docking sites or modes (Figure 3(a)). Six of the nine modes (nos. 1, 2, 4, 5, 7, 8, Figure 3) clustered at an interface between subunits within the NGAL trimer. The remaining binding modes (nos. 3, 6, 9, Figure 3) were at the periphery of a subunit near the NGAL calyx without being fully buried within the hydrophobic binding pocket. When apo-pyoverdine docked at the periphery of a subunit, the ligand was also near a tryptophan residue. Docking modes 6 and 9 for apo-pyoverdine were in proximity to both W31 and W79 as best seen in the 180° rotational view (Figure 3, right). In docking mode 3, apo-pyoverdine is only in proximity to W79. All nine of the docking modes for the apo-pyoverdine ligand localized to the outer surface of the NGAL monomer (data not shown). In eight of the nine docking modes, apo-pyoverdine docked near the opening of the calyx but not within the ligand-binding pocket.

3.2.2. Ferric Pyoverdine Docking to NGAL

As with the apo-pyoverdine structures, AutoDock Vina predicted nine positions in which ferric pyoverdine could dock into the NGAL trimer; however, none of the ferric pyoverdine molecules were predicted to bind near or inside NGAL’s ligand binding pocket (Figure 4(a)). Each ferric pyoverdine docked between the NGAL subunits. Further, AutoDock Vina predicted nine positions in which ferric pyoverdine could dock onto the NGAL monomer structure. In each case, however, the position of the docked ligand was at the surface of the protein and not near or inside the calyx (Figure 4(b)). This computational data strongly suggests that ferric pyoverdine does not bind to NGAL, which is consistent with the experimental results.

3.2.3. Enterobactin Docking to NGAL

Enterobactin, the major siderophore in E. coli and many other pathogenic enteric bacteria, binds to NGAL with very high affinity [18, 28]. The binding and molecular docking of ferric enterobactin to NGAL is well established and extensively investigated [18]. In order to validate the docking results obtained for pyoverdine, ferric enterobactin (Figure 1(d)) was used as a positive control for computational docking to the NGAL trimer (Figure 5) and monomer (data not shown). AutoDock Vina predicted 9 positions in which ferric enterobactin could dock into the trimeric NGAL protein crystal structure PDB 1L6M. Seven different docking modes (nos. 1, 2, 3, 4, 5, 6, and 9) indicate that ferric enterobactin could dock inside the NGAL calyx as been reported in the crystal structure of ferric enterobactin-bound NGAL [18]. The remaining two docking modes (nos. 7 and 8) predicted ferric enterobactin docking away from the calyx and between the NGAL subunits (Figure 5). Overall, AutoDock Vina accurately simulated the binding of ferric enterobactin to the NGAL trimer at the ligand binding site (i.e., inside the calyx of NGAL). Further, ferric pyoverdine was superimposed on ferric enterobactin docked into NGAL calyx with a transparent surface presentation. The results reflect again the inability of ferric pyoverdine to dock into the calyx (Figure 6).

3.3. Predicted Binding Affinities

AutoDock Vina provides a computed binding affinity for each docking mode predicted. Binding affinity data are summarized in Table 1. The more negative the value, the tighter the predicted binding. For the apo-pyoverdine interaction with the NGAL trimer, predicted binding affinities were very similar for the nine modes, ranging from −6.4 to −6.9 kcal/mol. Predicted binding affinities for apo-pyoverdine with the NGAL monomer (−5.8 to −6.5 kcal/mol) were slightly lower in range than those for the analogous trimer complex (Table 1). Data for ferric pyoverdine interactions gave slightly tighter predicted binding than that for apo-pyoverdine interactions: predicted binding energies of −8.2 to −9.2 kcal/mol for the trimer and −7.6 to −8.9 kcal/mol for the monomer. Consistent with biological trends, iron-bound siderophores were predicted to have stronger binding interactions with NGAL monomer or trimer than their corresponding apo-forms. In all permutations, docking mode 1 represented the highest affinity conformation of the ligand interacting with the NGAL. In every case, docking mode 1 was computed to have a higher binding energy in interaction with the trimer over the monomer, although even this mode for pyoverdine gave low predicted affinity. In contrast, the tightest predicted binding affinities were computed for the ferric enterobactin interactions with the NGAL trimer (−9.7 to −11.8 kcal/mol). While all of the predicted binding affinities are moderate in scale, none of the pyoverdine values are comparable to the range predicted for ferric enterobactin binding with NGAL.

4. Discussion

P. aeruginosa is an opportunistic pathogen that causes severe infection not only in CF patients but also in critically ill patients in intensive care units where it can cause ventilator-associated infections. It is a major nosocomial pathogen associated with significant morbidity and mortality (40–60%) in immunocompromised patients and in certain hospital units like burn and palliative care [29]. P. aeruginosa possesses several virulence factors, one of which is pyoverdine, the major siderophore that sequesters iron from the host to facilitate bacterial growth [30, 31]. Pyoverdine is a virulence factor required for establishing infection and biofilm formation [32]. In response to infection, the host secretes NGAL that sequesters bacterial siderophores thereby preventing bacteria from establishing infection [6].

The role of NGAL in innate immunity is well established and although NGAL-deficient mice develop normally, they are more susceptible to bacterial infections. NGAL is required for pulmonary host defense against Klebsiella pneumonaie [7], indicating that NGAL is an important mediator in innate host defense [8]. However, not all bacterial siderophores are recognized by NGAL and some pathogenic bacteria modify their siderophore structure to evade NGAL recognition [1113, 16, 17]. It has been shown that salmochelin, a major siderophore in pathogenic Salmonella, is a glucosylated form of enterobactin [11, 12, 33]. The glycosylation renders the molecule bulky such that it would not fit within the NGAL calyx, thus evading innate immune defense. Pyoverdine, a catecholate-based molecule, is not a glycosylated siderophore. P. aeruginosa does not contain the iroA gene cluster responsible for the glycosylation of catecholate-based siderophores that do not bind to NGAL [1, 3]. Compared to other dihydroxybenzoic-acid-(DHBA-) based catecholates that bind to NGAL such as enterobactin, pyoverdine is a structurally bulky molecule due to its dihydroxyquinoline chromophore [1]. We postulated that pyoverdine is a stealth siderophore that does not bind to NGAL. Here we show that pyoverdine evades NGAL recognition and therefore likely facilitates iron acquisition by P. aeruginosa which supports establishing infection.

Pyoverdine is uniquely produced by Pseudomonas bacteria. Using the tryptophan fluorescence-quenching method, we found that pyoverdine did not bind to NGAL, neither in the apo-form nor when complexed with iron. In contrast, enterobactin bound to NGAL causing strong tryptophan quenching. Using molecular docking we showed that pyoverdine did not dock into the calyx of NGAL while identical approaches indicate that enterobactin docks to NGAL with high affinity. We conclude that pyoverdine is a virulence factor and a stealth siderophore which allows P. aeruginosa to establish infection.

Lipocalins act as carrier proteins for a variety of ligands including small-molecule drugs [34, 35]. Other members of the lipocalin family have also been shown to act as siderocalins, that is, binding and sequestering bacterial siderophores [36, 37]. Tear lipocalin, also called LCN1, has been shown to sequester enterobactin but not pyoverdine which lends strong support to our results [38, 39]. Unlike NGAL that preferably binds catecholate-based siderophores [28], tear lipocalin can bind wide range of bacterial and fungal siderophores including catecholate-based and hydroxamate-based siderophores [39]. In contrast to NGAL and tear lipocalin, galline Ex-FABP, a lipocalin expressed in chicken serum and egg white, has been reported to sequester both glucosylated enterobactin and pyoverdine. It is also a sensor for lysophosphatidic acid and is reported to have dual ligand specificities [40]. The recently published crystal structure of galline Ex-FABP explained the dual ligand specificity due its expanded multichambered cavity [40]. NGAL is reported to be a rigid and quite stable lipocalin with a smaller calyx that does not undergo conformational changes upon ligand binding [18, 19] which explains why pyoverdine binds to Ex-FABP but not to NGAL. Similar to NGAL, galline Ex-FABP was shown to exert strong bacteriostatic activity against E. coli and B. subtilis but not against Pseudomonas aeruginosa due to low affinity of interaction between Ex-FABP and ferric pyoverdine [40]. We also performed similar bacterial growth inhibition assay with NGAL and P. aeruginosa and found only 10% reduction in bacterial growth (data not shown); therefore we concluded that NGAL is not bacteriostatic for P. aeruginosa.

P. aeruginosa produces another siderophore called pyochelin which is a small chatecholate-based molecule with iron-binding affinity [32]. Recently, Strong and coworkers crystallized NGAL bound to pyochelin (PDB 3U03, deposited at the RCSB protein data bank) and the crystal structure clearly shows that pyochelin docks inside the calyx of NGAL but without binding avidly. However, pyochelin is not essential for establishing infection since a P. aeruginosa mutant that lacks pyochelin successfully established infection and formed a biofilm [4]. Chronic infections caused by P. aeruginosa as seen in CF are associated with biofilms that protect the bacteria against the host killing mechanisms as well as increase resistance to antibiotics [4, 30]. In established biofilms of CF airways, P. aeruginosa thrives and downregulates pyoverdine since iron is abundant in the airway surface liquid due to chronic inflammation [41, 42]. Therefore, pyoverdine is important for establishing infections but not for well-established infections like that in biofilms.

In summary, we show that pyoverdine, the major siderophore in P. aeruginosa, is a stealth siderophore that evades NGAL recognition.

Authors’ Contribution

S. M. Zughaier, M. E. Peek, and N. A. McCarty conceived and designed the experiments. S. M. Zughaier, M. E. Peek, and A. Bhatnagar performed the experiments. S. M. Zughaier, M. E. Peek, and N. A. McCarty contributed the reagents/materials/analysis tool. S. M. Zughaier and M. E. Peek analyzed the data. S. M. Zughaier and M. E. Peek wrote the paper. M. Zughaier, M. E. Peek, and N. A. McCarty critically reviewed the paper.

Conflict of Interests

The authors have no conflict of interests.

Acknowledgments

The authors are very thankful to Dr. Roland Strong (Division of Basic Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA) for critically reviewing the paper. This work has been supported by an Emory-Egleston Children’s Research Center Grant to S. M. Zughaier.