Structure of Pseudomonas aeruginosa inosine 5′-monophosphate dehydrogenase

The crystal structure of inosine 5′-monophosphate dehydrogenase from P. aeruginosa has been determined to 2.25 Å resolution.


Introduction
Inosine 5 0 -monophosphate dehydrogenase (IMPDH; EC 1.1.1.205) catalyzes the first committed reaction in de novo guanine-nucleotide biosynthesis. This enzyme converts inosine 5 0 -monophosphate (IMP) to xanthosine 5 0 -monophosphate (XMP), which can then be converted to guanosine 5 0 -monophosphate (GMP) by GMP synthetase. Guanine nucleotides are essential for signal transduction and DNA/RNA synthesis and IMPDH therefore plays a pivotal role in the growth and proliferation of both prokaryotic and eukaryotic cells (Hedstrom, 2009). The enzyme has been the subject of thorough kinetic and crystallographic studies. Structures of IMPDH from a variety of organisms are known, including pathogenic bacteria such as Borrelia burgdorferi (McMillan et al., 2000;PDB entry 1eep), Bacillus anthracis (Center for Structural Genomics of Infectious Diseases, unpublished work; PDB entry 3tsb) and Streptococcus pyogenes (Zhang et al., 1999), the protozoan parasites Cryptosporidium parvum ; PDB entry 3ffs) and Tritrichomonas foetus (Whitby et al., 1997;PDB entry 1ak5), and human type II (Colby et al., 1999;PDB entry 1b30).
The enzyme mechanism is thought to start with a catalytic cysteine attacking the C2 position of IMP; an NAD + -dependent dehydrogenation then generates NADH and the covalent product intermediate E-XMP*. The expulsion of NADH promotes the insertion of a mobile flap to occupy the vacated cofactor-binding site. An arginine then acts as a general base catalyst and supports hydrolysis to convert E-XMP* into XMP and to recharge the enzyme (Hedstrom, 2009). IMPDH has been identified as a potential drug target for antimicrobial infections (Hedstrom et al., 2011). Progress has been made in developing highly potent inhibitors of the enzyme from C. parvum, which also display antiparasitic efficacy against Toxoplasma gondii (Gorla et al., 2012). Critically, these researchers have also addressed the issue of selectivity over the human homologue and a selectivity of greater than 1000-fold has been attained for the microbial enzymes.
Our interest is in the opportunistic Gram-negative pathogen Pseudomonas aeruginosa, which is widely recognized as a significant cause of hospital-acquired infections that often affect patients with compromised immune systems (Gaynes & Edwards, 2005). As part of a wide-ranging project to advance early-stage drug discovery against P. aeruginosa (Eadsforth et al., 2012;Moynie et al., 2013), we have adopted a structure-based approach (Hunter, 2009) to drive target assessment. Genetic evidence exists suggesting that the gene encoding IMPDH is essential in several bacteria, including P. aeruginosa (Liberati et al., 2006). We therefore sought to generate a recombinant source of the enzyme and to investigate its structure. Here, we report success in this endeavour and present the crystal structure of apo IMPDH from P. aeruginosa at 2.25 Å resolution. Comparisons with the structures of the enzyme from C. parvum (CpIMPDH) and Homo sapiens suggests that it may be possible to selectively inhibit the bacterial IMPDH over that of the host.

Protein expression and purification
The full-length gene (PA3770) for P. aeruginosa IMPDH (PaIMPDH; amino acids 1-489; UniProt ID Q9HXM5) was synthesized and codon-optimized for expression in Escherichia coli (GenScript). It was cloned into a modified pET15b (Novagen) cloning vector which encodes an N-terminal His tag followed by a tobacco etch virus (TEV) protease cleavage site and the protein of interest (Rao et al., 2011). The integrity of the construct was verified by DNA sequencing (DNA Sequencing Unit, University of Dundee). Recombinant PaIMPDH was produced in E. coli BL21 (DE3) pLysS cells (Stratagene) at 310 K in LB medium supplemented with 100 mg ml À1 ampicillin. The cells were grown to an OD 600 of $0.7 before induction with isopropyl -d-1-thiogalactopyranoside at a final concentration of 1 mM and growth at 293 K for a further 24 h. The cells were harvested by centrifugation (3500g at 277 K for 40 min). The cell pellet was resuspended in buffer (25 mM Tris-HCl pH 8.5, 500 mM NaCl, 20 mM imidazole) supplemented with an EDTA-free protease-inhibitor cocktail tablet (Calbiochem. The cells were lysed by passage through a continuous-flow cell disruptor (Constant Systems) at 207 MPa and the cell debris was removed following centrifugation (40 000g at 277 K for 20 min). PaIMPDH was purified using affinity chromatography on a 5 ml HisTrap HP column (GE Healthcare) pre-charged with Ni 2+ and equilibrated in buffer A (25 mM Tris-HCl pH 8.5, 500 mM NaCl, 20 mM imidazole). A linear concentration gradient of imidazole was applied using buffer B (25 mM Tris-HCl pH 8.5, 500 mM NaCl, 500 mM imidazole) to elute the protein, which was then dialyzed against buffer C (25 mM Tris-HCl pH 7.5, 250 mM NaCl) at 277 K. Fractions were analyzed using SDS-PAGE and those containing PaIMPDH were pooled. The protein was further purified, with the His tag still present, by sizeexclusion chromatography using a Superdex 75 26/60 column (GE Healthcare) equilibrated with buffer C on an Ä KTApurifier (GE Healthcare). This column had previously been calibrated with molecular-weight standards: blue dextran (>2000 kDa), thyroglobulin (669 kDa), ferritin (440 kDa), aldolase (158 kDa), conalbumin (75 kDa), ovalbumin (43 kDa), carbonic anhydrase (29.5 kDa), ribonuclease A (13.7 kDa) and aprotinin (6.5 kDa) (GE Healthcare; data not shown). The protein eluted as a single species of approximate mass 200 kDa. Fractions containing the protein were pooled and concentrated to 22 mg ml À1 using Amicon Ultra devices (Millipore) for subsequent use. A final yield of 25 mg per litre of bacterial culture was obtained. The purity of the protein was confirmed by SDS-PAGE and mass spectrometry (Fingerprint Proteomics Facility, University of Dundee). A theoretical extinction coefficient of 22 450 M À1 cm À1 at 280 nm was used to estimate the protein concentration (ProtParam; Gasteiger et al., 2005); the theoretical mass of one subunit was estimated as 54 kDa, with a calculated isoelectric point of 6.8. The purified protein was stored at 277 K in buffer C for subsequent use.

Crystallization, data collection and structure determination
A variety of commercial screens were employed to identify appropriate starting conditions for crystal growth. Diffraction-quality crystals of PaIMPDH were obtained at room temperature by sittingdrop vapour diffusion with a reservoir solution consisting of 0.2 M sodium acetate pH 4.5, 40% 2-methyl-2,4-pentanediol (MPD). Prismatic crystals appeared after a month from a mixture of equal volumes of protein solution (9 mg ml À1 in 25 mM Tris-HCl pH 7.5, 250 mM NaCl from the stock protein solution) and reservoir solution. The concentration of protein in the drop is therefore about 4.5 mg ml À1 . The crystals grew to a maximum dimension of approximately 0.15 mm. The high concentration of MPD in the reservoir solution proved to be a suitable cryoprotectant and data were collected at 100 K using a Rigaku MicroMax-007 rotating-anode X-ray generator (Cu K, = 1.5418 Å ) coupled to an R-AXIS IV ++ image-plate detector. The data were indexed with XDS (Kabsch, 2010) and scaled using SCALA (Evans, 2006) from the CCP4 program suite .   (hkl) is the intensity of the ith measurement of reflection hkl and hI(hkl)i is the mean value of I i (hkl) for all i measurements. ‡ R work = P hkl jF obs j À jF calc j = P hkl jF obs j, where F obs is the observed structure factor and F calc is the calculated structure factor. § R free is the same as R cryst except that it was calculated with a subset (5%) of data that were excluded from the refinement calculations. } Diffraction precision index (Cruickshank, 1999). † † Engh & Huber (1991).
(TLS) was applied (Painter & Merritt, 2006). The refinement proceeded with the incorporation of water molecules following conservative criteria (Leonard & Hunter, 1993), a chloride and a number of side chains with dual conformers. The refinement was terminated when there were no significant changes in the R work and R free values and when inspection of the difference map suggested that no further corrections or additions were required. The stereochemistry and quality of the model were validated using MolProbity (Chen et al., 2010). Data-collection and structure-refinement statistics are shown in Table 1.

Structure determination
A highly efficient recombinant protein-expression system has been prepared and protocols for purification and crystallization have been established. The crystal structure of PaIMPDH was determined at 2.25 Å resolution. The asymmetric unit consists of a single polypeptide chain comprising residues 1-91, 205-370 and 426-467, with an estimated solvent content of 60% and a V M (Matthews, 1968) of 2.94 Å 3 Da À1 . The initial model consisted of 261 residues, with a correlation coefficient of 0.75 and R work and R free values of 37.2% and 41.8%, respectively. Subsequent model building and refinement extended this to 293 residues, with a correlation coefficient of 0.96 and improved R work and R free values of 14.8% and 19.0%, respec-    tively. Analysis of the Ramachandran plot revealed that 99.3% of the residues were in the allowed region, with only two outliers: Gly226 and Gln458. These two residues are located on flexible loops linking 8 to 9 and 13 to 12, respectively, and the associated electron density is relatively poor. A greater degree of disorder was evident at several positions, where it was not even possible to interpret the electron density. Consequently, 196 residues are absent from the final model. A further comment on this will be made below.

Overall structure
The crystal structures of IMPDHs from several species have previously been determined and an excellent review has been published (Hedstrom, 2009). In brief, IMPDH forms homotetramers. The subunit is constructed from two segments: a catalytic domain which displays the (/) 8 -barrel fold and a subdomain containing two cystathionine -synthase (CBS)/Bateman domains (Bateman, 1997;Hedstrom, 2009). The function of the CBS domains in this enzyme system is a source of debate (Hedstrom et al., 2011). It has been determined that the deletion of these subdomains has no effect on the enzymatic activity (Nimmesgern et al., 1999). This subdomain protrudes out from the catalytic domain and is usually disordered in IMPDH structures (Hedstrom, 2009).

The active site
The active site of PaIMPDH is located at the C-terminal end of the -barrel, as is common for this (/) 8 -barrel fold (Hall et al., 2002), and at the edge there is a flexible loop (residues 297-315; Fig. 2b). This loop adopts different conformations at different stages of the catalytic cycle and is often disordered in crystal structures (Hedstrom, 2009). Recent studies have shown that the active-site loop controls the transition between the dehydrogenase and hydrolase activities (Josephine et al., 2010). This loop is fully ordered in PaIMPDH, with an average B factor of 33.6 Å 2 , which is comparable to the average for the whole protein (33.7 Å 2 ), and adopts an 'open' conformation in contrast to the 'closed' conformation observed in other IMPDH structures (Fig. 2b).

Potential for selective inhibition of PaIMPDH
A sequence alignment of PaIMPDH with selected homologues shows that the key residues involved in catalysis are conserved across species (Fig. 3). However, variation in the residues interacting with the adenosine and pyrophosphate portions of the NAD + cofactor have been observed (Fig. 3) and are believed to be good targets for achieving selective IMPDH inhibition (Hedstrom et al., 2011). Furthermore, PaIMPDH possesses a structural motif that indicates that this enzyme belongs to a group of IMPDHs that are susceptible to inhibitors developed to target CpIMPDH Sequence alignment of PaIMPDH with homologues from T. foetus (GenBank ABI11203.1), H. sapiens (NCBI Reference Sequence NP_000875.2), B. anthracis (UniProt Q81W29) and S. pyogenes (UniProt P0C0H6). Residues involved in interactions with IMP, including the catalytic cysteine, are coloured green. Conserved residues interacting with NAD + are coloured blue and nonconserved residues interacting with the adenosine and pyrophosphate portions are coloured magenta and orange, respectively. Two conserved residues on the mobile flap are coloured red. Residues involved in selectively binding the inhibitor C64 are highlighted in yellow. Residues missing in the PaIMPDH structure have a dashed line beneath them. The sequence alignment was generated and annotated using ClustalW (Larkin et al., 2007) and ALINE (Bond & Schü ttelkopf, 2009), respectively. residues that interact with the bromoaniline moiety of C64: Ala165 and Tyr358 (Gollapalli et al., 2010;Fig. 4). The alanine and tyrosine residues create a hydrophobic surface that forms van der Waals interactions with the bromoaniline group. These residues are conserved in bacterial homologues, including PaIMPDH (Ala249 and Tyr446; Figs. 3 and 4), but differ in the human enzyme (Ser276 and Asp470; Fig. 3). A further seven residues in CpIMPDH are implicated in assisting C64 binding and these are Ser22, Pro26, Thr221, Ser354, Gly357, Met308 and Glu329 (data not shown). In PaIMPDH these residues correspond to Gly19, Pro25, Thr306, Ala442, Gly445, Met393 and Glu417 (Fig. 3). This high degree of conservation extends to other bacterial homologues but not to the human enzyme (Fig. 3). This reinforces the idea that selectivity for the bacterial enzymes over human IMPDH is possible. With an efficient recombinant expression system for PaIMPDH and established crystallization conditions, the scene is set for a structure-based approach, noting that the reagents and data generated in the study of protozoan homologues might now be exploited to fast-track the discovery of new IMPDH inhibitors relevant to antibacterial drug discovery.