1.6 Å Crystal structure of a PA2721 protein from <i>pseudomonas aeruginosa</i>—A potential drug‐resistance protein

Drug-resistant microorganisms are an important health issue, as many human pathogens have gained resistance to a wide range of commercially available pharmaceuticals. A number of targets from pathogenic organisms have been selected for structural studies at the Midwest Structural Biology Center with the objective to broaden our knowledge about these organisms, acquire information on the mechanisms that trigger their resistances to known antibiotics, and potentially provide a foundation for drug discovery. Pseudomonas aeruginosa is an opportunistic human pathogen that has a natural resistance to many antibiotics and disinfectants.1 It causes infections throughout the body, including the urinary tract, respiratory system, the skin, soft tissue, bones, and joints. Herein, we report the crystal structure of PA2721, a 17.6 kDa uncharacterized protein from P. aeruginosa PAO1. The structure of PA2721 has been determined and refined to 1.60 Å by the single wavelength anomalous dispersion (SAD) method using a selenomethionine (SeMet) derivative of PA2721. Data collection, refinement, and model validation statistics are summarized in Table I. The final model includes four protein molecules in the asymmetric unit that are organized into two pairs of homodimers [Fig. 1(a)]. The dimension of the dimer is approximately 71 × 43 × 30 Å. Each monomer is composed of eight β-strands (arranged in two β-sheets β1–β4 and β5–β8) and five α-helices, and is organized into two domains. The N-terminal domain is made up of a β1↑α1↓β2↑β3↓β4↑ motif, whereas the C-terminal domain has a β5↓α2↑β6↓β7↑β8↓ motif plus an extra three α-helices at the C-terminus (α3–α5) [Fig. 1(b)]. The repeat of the βαβββ motif within the monomer could indicate the possibility that this structure evolved by gene duplication and fusion events, as it was suggested by Bergdoll et al.2 for enzymes containing a similar fold. A: Ribbon diagrams of the dimer of PA2721viewed in a side orientation (top) and viewed down the noncrystallographic twofold relating the subunits (bottom). B: Diagram showing the secondary structure elements and primary sequence of PA2721. The strands that belong to each β-sheet are labeled A,B; β-hairpins are indicated as red loops, β-bulges are marked by β, γ-turns are marked by γ.20 Although there is no significant sequence similarity between the domains, 42 Cα pairs superimposed with a root-mean-square deviation (RMSD) of 1.67 Å. Crystallographic packing analysis, as well as the solvent accessible surface area calculation using AreaIMol from the CCP4 program suite, showed extensive interaction between tightly associated monomers with ∼1,900 Å2 of buried surface of the monomer.3 Therefore, it is likely that the homodimer is the biologically relevant oligomeric form. The dimer interface consists of an antiparallel interaction between the β4-strand of the N-terminal domain from one subunit and the β5-strand of the C-terminal domain of the opposite subunit, which leads to the formation of two curved 8-stranded β-sheets [Fig. 1(a)]. The fold recognition programs 3DHit and DALI server revealed that PA2721 bore significant similarity to members of the βαβββ superfamily such as the fosfomycin resistance proteins FosA and FosX and the metalloenzyme glyoxalase I (Fig. 2).4, 5 The closest matches are the structures of fosfomycin resistance protein [Protein Data Bank (PDB) code: 1LQP; Z-score 5.6, sequence identity 20.6%, RMSD 2.35 Å], and Ni(II)-bound glyoxalase I (PDB code: 1F9Z; Z-score 5.3, sequence identity 18%, RMSD 2.28 Å).6 Structural analysis of the members of the βαβββ superfamily indicates that substrate binding sites of these enzymes are located at the subunit–subunit interface, within the inside cleft formed by the curved 8-stranded β-sheets.2 Despite the low sequence similarity among members of the family and the high diversity of reactions they catalyze, there is good conservation of the active center residues and their location within the enzymes. This conservation is particularly strong among members that require metal ions for catalytic activity. For example, the residues involved in Fe+2 coordination in extradiol dioxygenase (E.C.1.13.11.39) (H146, H210, and E260) superimposed within 0.7–1.4 Å (Cα atoms) of the Zn+2 coordinating residues in glyoxalase I (Q33, E99, H126, E172) (E.C.4.4.1.5).2, 7 However, PA2721 is unlikely to bind metal ions because superimposition of PA2721 with these enzymes using the SSM server showed that residues H35, G146, F76, E117 occupy equivalent positions.8 These residues are unlikely to provide a good environment for metal ion coordination. Ribbon diagram of a superposition of a monomer of PA2721 (red) with a monomer of the fosfomycin resistance protein (FosA) (PDB code: 1lqk) (blue) and glyoxylase I (PDB code: 1f9z) which have been identified as the closest matching folds. The images were generated with the program PyMol.21 To determine possible functionally significant motifs/active centers in PA2721, we applied a nest analysis as implemented in Profunc and active center pattern recognition searches using the program pvSoar. Profunc identified two functionally significant structural motifs formed by the following residues: S13, A14, A15, and W112, G113, I114.9, 10 The comparison of active center pockets, surface patterns, and conservation of structurally important residues between PA2721 and structures in PDB database with pvSoar revealed a pocket on the PA2721 surface that is similar to the active center pocket of phosphatidylinositol-specific phospholipase C from Bacillus cereus (PDB code: 5PTD).11 The PA2721 pocket is positioned at the domain–domain interface, within the monomer and is lined by Y11, P68, R71, H72 in addition to the residues identified by Profunc. The active center pocket of the phospholipase C is positioned at the C-terminal end of the β-barrel with residues H32, R69, H82, K115, R163, W178, Y200 involved in interaction with the inhibitor inositol.12 Calculation of solvent-accessible areas for the active center revealed that in PA2721 the pocket has a solvent accessible area 204.6 A2 and an accessible volume 178.1 A3 in comparison to a solvent accessible area 197.721 A2 and an accessible volume 111.070 A3 for the active site pocket of phosphatidylinositol-specific phospholipase C. Therefore, both ProFunc and pvSoar detected similar regions in the protein as functionally important; however, the chemical nature of the potential ligand as well as the specific biochemical function of PA2721 is not clear. In summary, the PA2721 protein has been identified as a new member of βαβββ family with structural similarity to fosfomycin resistance protein and glyoxalase I, and has a potential active site similar to phosphatidylinositol-specific phospholipase C, suggesting that it binds a phosphatidylinositol-like compound. Further biochemical and biophysical studies are needed to gain detailed understanding of this protein's function. The ORF of PA2721 was amplified by polymerase chain reaction from Escherichia coli genomic DNA (ATCC). The gene was cloned into the NdeI and BamHI sites of a modified pET15b cloning vector (Novagen) in which the TEV protease cleavage site replaced the thrombin cleavage site and a double stop codon was introduced downstream from the BamHI site. This construct provides an N-terminal His6-tag separated from the gene by a TEV protease recognition site (ENLYFQ↓G). The fusion protein was overexpressed in E. coli BL21-Gold [DE3 (Stratagene)] harboring an extra plasmid encoding three rare tRNAs (AGG and AGA for Arg, ATA for Ile). The cells were grown in LB at 37°C to an OD600 of ∼0.6 and protein expression induced with 0.4 mM isopropyl-β-D-thiogalactopyranosid. After induction, the cells were incubated with shaking overnight at 15°C. The harvested cells were resuspended in binding buffer (500 mM NaCl, 5% glycerol, 50 mM HEPES pH 7.5, 5 mM imidazole), flash-frozen in liquid N2, and stored at −70°C. The thawed cells were lysed by sonication after the addition of 0.5% NP-40 and 1 mM each of PMSF and benzamidine. The lysate was clarified by centrifugation (30 min at 27,000g) and passed through a DE52 column preequilibrated in binding buffer. The flow-through fraction was then applied to a metal chelate affinity-column charged with Ni2+. The His6-tagged protein was eluted from the column in elution buffer (500 mM NaCl, 5% glycerol, 50 mM HEPES pH 7.5, 500 mM imidazole) and the tag then cleaved from the protein by treatment with recombinant His-tagged TEV protease. The cleaved protein was then resolved from the cleaved His6-tag and the His6-tagged protease by flowing the mixture through a second Ni2+-column. The PA2721 protein was dialyzed in 10 mM HEPES pH 7.5, 500 mM NaCl, and concentrated using a BioMax concentrator (Millipore). Before crystallization, any particulate matter was removed from the sample by passing through a 0.2-μm Ultrafree-MC centrifugal filter (Millipore). For the preparation of SeMet labeled protein, PA2721 was expressed in the methionine auxotroph strain B834(DE3) of E. coli (Novagen) and purified under the same conditions as the native protein in supplemented M9 media. The reducing reagent β-mercaptoethanol (5 mM) was added to all purification buffers. The protein was crystallized by the vapor diffusion method in hanging drops by mixing 2 μL of the protein solution with 2 μL of 0.2 M di-ammonium hydrogen citrate and 20% polyethylene glycol 3350, and equilibrated at 20°C over 500 μL of this solution. PA2721 crystallizes in the triclinic space group P1 with a cell dimension a = 44.55 Å, b = 55.02 Å, c = 67.07 Å, α = 70.71°, β = 89.89°, γ = 85.61°. Before data collection, crystals were soaked for ∼20 s in a cryoprotectant consisting of a 25% ethylene glycol in the crystal mother liquor, and then flash-frozen in liquid nitrogen. Diffraction data were collected at 100 K at the 19ID beamline of the Structural Biology Center at the Advanced Photon Source, Argonne National Laboratory. SAD data to 1.6 Å [12.6603 keV (0.9794 Å)] were collected from a SeMet-substituted protein crystal as described earlier by Walsh et al.13 One crystal was used for data collection at 100 K and a total oscillation range of 360° was obtained using inverse beam geometry. Data were integrated, reduced, and scaled with the HKL2000 suite.14 Data statistics are summarized in Table I. The structure was determined by SAD phasing using autoSHARP and refined against the peak data with REFMAC5.2, 15 The initial model was autotraced with 1.6 Å data using the program ARP/wARP.16 The partial model obtained from ARP/wARP was subsequently manually extended, and rebuilt with O and COOT.17, 18 The electron density map contoured at 1.0 σ is well connected for most of the main chain, except for the surface-exposed residues 38–48 of molecules C and D, which have been omitted from the model. The stereochemistry of the structure was checked with PROCHECK.19 Final refinement statistics are presented in Table I. Atomic coordinates have been deposited into the PDB as 1U69. The authors thank Andrew Binkowski for the use of pvSoar server and all members of the Structural Biology Center at Argonne National Laboratory for their help in conducting experiments and Lindy Keller for help in preparation of the manuscript.

Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.