3-Hydroxybenzoate 6-Hydroxylase from Rhodococcus jostii RHA1 Contains a Phosphatidylinositol Cofactor

Abstract

3-Hydroxybenzoate 6-hydroxylase (3HB6H, EC 1.13.14.26) is a FAD-dependent monooxygenase involved in the catabolism of aromatic compounds in soil microorganisms. 3HB6H is unique among flavoprotein hydroxylases in that it harbors a phospholipid ligand. The purified protein obtained from expressing the gene encoding 3HB6H from Rhodococcus jostii RHA1 in the host Escherichia coli contains a mixture of phosphatidylglycerol and phosphatidylethanolamine, which are the major constituents of E. coli’s cytoplasmic membrane. Here, we purified 3HB6H (RjHB6H) produced in the host R. jostii RHA#2 by employing a newly developed actinomycete expression system. Biochemical and biophysical analysis revealed that Rj3HB6H possesses similar catalytic and structural features as 3HB6H, but now contains phosphatidylinositol, which is a specific constituent of actinomycete membranes. Native mass spectrometry suggests that the lipid cofactor stabilizes monomer-monomer contact. Lipid analysis of 3HB6H from Pseudomonas alcaligenes NCIMB 9867 (Pa3HB6H) produced in E. coli supports the conclusion that 3HB6H enzymes have an intrinsic ability to bind phospholipids with different specificity, reflecting the membrane composition of their bacterial host.

Introduction

Rhodococcus jostii RHA1 is a biotechnologically and environmentally important bacterium from the order Actinomycetales. Together with the genera Nocardia, Corynebacterium and Mycobacterium, Rhodococcus forms a distinct group of bacteria called mycolata (Finnerty, 1992; Brennan and Nikaido, 1995; Chun et al., 1996; Gürtler et al., 2004), characterized by a complex cell envelope (Sutcliffe, 1998; Guerin et al., 2010; De Carvalho et al., 2014) and an impressive catabolic diversity, allowing adaptation to different carbon sources for growth (van der Geize and Dijkhuizen, 2004). In comparison with other mycolata, R. jostii RHA1 is particularly rich in oxygenases (203 putative genes) and ligases (192 putative genes), gained primarily through ancient gene duplications or acquisitions (McLeod et al., 2006; Yam et al., 2010).

We recently reported the crystal structure of R. jostii RHA1 3-hydroxybenzoate 6-hydroxylase (3HB6H), produced as a recombinant protein in Escherichia coli (Montersino et al., 2013). 3HB6H (EC 1.13.14.26) is a NADH and FAD-dependent monooxygenase that catalyzes the para-hydroxylation of 3-hydroxybenzoate to 2,5-dihydroxybenzoate, using a Tyr-His pair for substrate binding and catalysis (Sucharitakul et al., 2015). The crystal structure analysis revealed that 3HB6H has the conserved fold of group A flavoprotein hydroxylases (Montersino et al., 2011; Huijbers et al., 2014), but differs from the other family members in additional binding of phospholipids. The tightly bound phospholipids were identified as a mixture of PG and PE, which are the major constituents of the E. coli cytoplasmic membrane (Pulfer and Murphy, 2003; Oursel et al., 2007). The fatty acyl chains of the phospholipid ligands of 3HB6H protrude into the substrate-binding pockets, whereas the surface-exposed hydrophilic headgroups are involved in enzyme dimerization (Figure 1) (Montersino et al., 2013).

FIGURE 1

To shed more light on the role of these lipid guests, bearing in mind the different lipid compositions of Gram-positive and Gram-negative bacterial membranes (Finnerty, 1992; Sutcliffe, 1998), in the present work we produced Rj3HB6H in a newly developed R. jostii RHA1#2 expression strain and, in addition, 3HB6H from Pseudomonas alcaligenes NCIMB 9867 (Pa3HB6H) in E. coli. Biochemical and biophysical characterization revealed that Rj3HB6H possesses similar catalytic and structural features as 3HB6H, but contains PI as glycerophospholipid ligand. Lipid analysis of Pa3HB6H indicates that lipid binding is an intrinsic property of prokaryotic 3-hydroxybenzoate 6-hydroxylases.

Materials and Methods

Chemicals

Aromatic compounds were purchased from Sigma-Aldrich (St Louis, MO, United States) and Acros Organics (Morris Plains, NJ, United States). Catalase, FAD, FMN, arabinose, antibiotics, Terrific broth (TB) and LB broth (Miller) (LB) were from Sigma-Aldrich (St Louis, MO, United States). Pefabloc SC and DNase I were obtained from Roche Diagnostics GmbH (Mannheim, Germany). Restriction enzymes and Pfu DNA polymerase were from Thermo Fischer Scientific (United States). 4-androstene-3,17-dione was from Merck (Oss, Netherlands). Crystallization kits were purchased from New Hampton (Aliso Viejo, CA, United States). Immobilized metal affinity chromatography columns (His GraviTrap) were from GE Healthcare Bioscience AB (Uppsala, Sweden). All other chemicals were from commercial sources and of the purest grade available.

Bacterial Strains and Primers

All bacterial strains and primers used in this study are listed in Tables 1, 2.

Construction of Rhodococcus Expression Vector Q2+

The E. coli-Rhodococcus shuttle vector pRESQ (van der Geize et al., 2002) was modified by insertion of the RP4 oriT of pK18mobsacB (Schäfer et al., 1994) enabling trans-conjugal transfer of the resulting vector. For this, the oriT-region was amplified from pK18mobsacB by PCR using forward primer oriT-F and reverse primer oriT-R (Table 1). The obtained 549 bp PCR-product was cloned into the SmaI-site of pRESQ, resulting in pQmob. A duplicate region of 424 bp on pQmob was removed by deleting the 760 bp XbaI-BspHI fragment, yielding pQmobΔd. The egfp gene from pIJ8630 (Sun et al., 1999) was amplified by PCR using forward primer egfp-PciI-F, containing a PciI restriction site, and reverse primer egfp-PciI-R, also containing a PciI restriction site (Table 1). The 744 bp PciI-PciI fragment containing the egfp gene was cloned into the PciI-site of pQmobΔd to generate peGFPQ.

The R. jostii strain RHA1 genomic region consisting of gene ro00440, its promoter region and the prmA promoter (PprmA) (here referred to as region reg1-PprmA; GenBank accession number CP000431: nt 521345 - nt 523358) was amplified from genomic DNA of R. jostii RHA1 by PCR using forward primer reg1-BglII-F, containing a BglII restriction site, and reverse primer prmA-NdeI-R, containing an NdeI restriction site (Table 1). The 2014 bp BglII-NdeI reg1-PprmA fragment was cloned into the BglII-NdeI sites of peGFPQ, yielding prMOeGFPQ1.

The R. jostii RHA1 gene ro00452 and its promoter region (here referred to as region reg2; CP000431: nt 534363 – nt 536227) were amplified by PCR using forward primer reg2-PciI-F, containing a PciI restriction site and reverse primer reg2-PciI-R, also containing a PciI restriction site (Table 1). The 1880 bp PciI-PciI reg2 fragment was cloned into the PciI-site of prMOeGFPQ1, resulting in prMOeGFPQ2.

For construction of expression vector Q2+, the egfp gene of prMOeGFPQ2 was replaced with a multiple cloning site (MCS). For this, part of the MCS of pBluescript KS was amplified by PCR using forward primer MCS-NdeI-F, containing an NdeI restriction site and reverse primer MCS-PciI-R, containing a PciI restriction site (Table 1). The 125 bp NdeI-PciI MCS fragment was cloned into the NdeI–PciI site of prMOeGFPQ2, replacing the NdeI–PciI region containing egfp, resulting in the expression vector Q2+.

Cloning and Production of 3HB6H in R. jostii RHA1#2

The 1321 bp NdeI-XmnI fragment of pBAD-3HB6H-His₆ (Montersino and van Berkel, 2012) containing the 3HB6H gene including the His₆-tag, was cloned into the NdeI-HindIII (Klenow-fragment treated) site of expression vector Q2+ to generate Q2+-3HB6H-His₆.

Rhodococcus cells were electroporated as described previously (van der Geize et al., 2000). Prior to electroporation, plasmid DNA was desalted by dialyzing 10 μL plasmid DNA for 30 min on a Millipore “V” Series filter disk (0.025 μm) floating on MiliQ water.

Cultures of R. jostii RHA1#2 were grown in LB broth supplemented with 50 μg⋅mL^-1 kanamycin at 30°C at 200 rpm. R. jostii RHA1#2 cells harboring Q2+3HB6H-His₆ were grown overnight in 3 mL LB broth, diluted 1:300 in 300 mL LB broth in a 3 L Erlenmeyer flask and grown for 20–24 h. Cultures were induced by adding 2 mM 4-androstene-3,17-dione dissolved in acetone. Growth was continued for 48 h after induction. Cells were harvested by centrifugation at 4°C and pellets were washed once with ice-cold 20 mM potassium phosphate, pH 7.2, containing 300 mM NaCl. After centrifugation at 4°C, cells were stored at -20°C.

Cloning and Production of 3HB6H from Pseudomonas alcaligenes NCIMB 9867

The xlnD gene sequence encoding for Pa3HB6H (UniProt: Q9F131) was synthesized and subcloned in a pBAD vector by GeneArt (Invitrogen, Carlsbad, CA, United States). The resulting construct (pBAD-Pa3HB6H-His₆) was verified by automated sequencing of both strands and electroporated into E. coli TOP10 cells for recombinant expression.

For enzyme production, E. coli TOP10 cells, harboring the pBAD-Pa3HB6H-His₆ plasmid, were grown in TB medium at 37°C supplemented with 100 μg⋅mL^-1 ampicillin until an optical density (OD_600nm) of 0.8 was reached. Expression was induced by the addition of 0.02% (w/v) arabinose and incubation was continued for 40 h at 17°C. Cells were harvested by centrifugation at 4°C and stored at -20°C.

Enzyme Purification

Rj3HB6H was purified to apparent homogeneity using an Åkta Explorer chromatography system (GE-Healthcare). R. jostii RHA1#2 cells containing the recombinant protein were suspended in ice-cold 20 mM potassium phosphate, pH 7.2, containing 300 mM NaCl, 1 mM Pefabloc SC, 1 mg DNAse and 100 μM MgCl₂, and subsequently passed twice through a precooled French Pressure cell (SLM Aminco, SLM Instruments, Urbana, IL, United States) at 16,000 psi. The resulting homogenate was centrifuged at 25,000 × g for 45 min at 4°C to remove cell debris, and the supernatant was applied onto a Ni-NTA agarose column (16 mm × 50 mm) equilibrated with 20 mM potassium phosphate, pH 7.2, containing 300 mM NaCl. After washing with five volumes of equilibration buffer, the enzyme was eluted with 300 mM imidazole in equilibration buffer. The resulting Rj3HB6H fraction was supplemented with 100 μM FAD and loaded onto a Source Q-15 anion exchange column (16 mm × 90 mm), pre-equilibrated with 50 mM Bis-Tris, 0.1 mM EDTA, pH 7.2. After washing with two volumes of equilibrating buffer, the enzyme was eluted with a linear gradient of 0–1 M NaCl in the same buffer. Active fractions were pooled, concentrated to 10 mg⋅mL^-1 using ultrafiltration (Amicon 30 kDa cutoff filter), and applied onto a Superdex S-200 (26 mm × 600 mm) column running in 50 mM potassium phosphate, 150 mM NaCl, pH 7.2. Active fractions were concentrated to 10 mg⋅mL^-1 using ultrafiltration (Amicon 30 kDa cutoff filter) and dialyzed at 4°C against 50 mM Bis-Tris, pH 7.2. The final Rj3HB6H preparation showed a single band after SDS-PAGE. The specific activity of the purified enzyme was 21 U mg^-1 using the standard activity assay (Table 3A).

Pa3HB6H was purified to apparent homogeneity, applying essentially the same procedure as described above for Rj3HB6H. The final Pa3HB6H preparation showed a single band after SDS-PAGE. The specific activity of the purified enzyme was 34 U⋅mg^-1 using the standard activity assay (Table 3B).

Purified enzymes were flash frozen in 1 mL aliquots in liquid nitrogen and stored at -80°C. Before use, thawed enzyme samples were incubated with 0.1 mM FAD and excess FAD was removed using a gel filtration column (10 mm × 100 mm) containing Bio-Gel P-6DG.

Biochemical Characterization

Molar absorption coefficients of protein-bound FAD were determined from absorption spectra of Rj3HB6H and Pa3HB6H recorded in the presence and absence of 0.1% (w/v) SDS, assuming a molar absorption coefficient for free FAD of 11.3 mM^-1⋅cm^-1 at 450 nm. The enzyme concentration of Rj3HB6H was determined by measuring the absorbance at 453 nm using a molar absorption coefficient for protein-bound FAD of 10.3 mM^-1⋅cm^-1. The enzyme concentration of Pa3HB6H was determined by measuring the absorbance at 450 nm using a molar absorption coefficient for protein-bound FAD of 11.0 mM^-1⋅cm^-1. Rj3HB6H and Pa3HB6H activity was determined at 25°C by measuring NADH consumption at 360 nm (Montersino and van Berkel, 2012). The standard assay mixture contained 50 mM Tris-SO₄, pH 8.0, 200 μM 4-hydroxybenzoate and 250 μM NADH. Steady-state kinetic parameters were determined from measurements at 25°C in 50 mM Tris-SO₄, pH 8.0. Hydroxylation efficiencies were determined by oxygen consumption experiments, essentially as described before (Montersino and van Berkel, 2012).

Crystallization and Structure Determination

Crystals of Rj3HB6H for structure determination were obtained by the sitting drop vapor diffusion method at 20°C by mixing equal volumes (2 μL) of protein and reservoir solutions. Protein solutions consisted of 30 mg⋅mL^-1 enzyme in 1 mM FAD, 2 mM 3-hydroxybenzoate, and 50 mM Bis-Tris, pH 7.2, whereas precipitant solutions consisted of 30% PEG 4000, 0.2 M lithium sulfate, and 0.1 M Tris-HCl, pH 8.5. Yellow crystals grew in 1 day.

X-ray diffraction data were collected at Grenoble and processed with the CCP4 package (Winn et al., 2011). The Rj3HB6H structure was solved by molecular replacement using the structure of a monomer of 3HB6H (pdb entry: 4BJZ) as search model. Crystallographic computing and model analysis were performed with COOT (Emsley et al., 2010), PHENIX (Adams et al., 2010) and the CCP4 package (Potterton et al., 2004). Pictures were generated with Pymol (Schrodinger, 2015) and CCP4 (Potterton et al., 2004). Data collection parameters and refinement statistics are presented in Table 4.

The atomic coordinates and structure factors of Rj3HB6H (code 5HYM) have been deposited in the Protein Data Bank¹.

Lipid Identification and Native ESI-MS Experiments

Extraction and identification of protein-bound lipids from Rj3HB6H and Pa3HB6H was performed as described for 3HB6H (Montersino et al., 2013). For nanoflow ESI-MS analysis under native conditions, enzyme samples were prepared in 50 mM ammonium acetate, pH 6.8. For analysis under denaturing conditions, enzyme samples were diluted either in 50% acetonitrile with 0.2% formic acid or in 5% formic acid. Native MS analysis was performed using a LC-T nanoflow ESI orthogonal TOF mass spectrometer (Micromass, Manchester, United Kingdom) in positive ion mode with a capillary voltage of 1.3 kV. The cone voltage was varied between 90 and 150 V and source pressure was set to 6.9 mbar to enhance transmission of large ions. Lipid identification was performed using a Quattro Ultima nanoflow triple quadrupole mass spectrometer (Micromass, Manchester, United Kingdom) in negative ion mode, with a capillary voltage of 1.3 kV and a cone voltage of 150 V. For MS/MS analysis, argon was supplied in the collision cell (2.0 × 10^-3 bar). Collision energy was adjusted to gain optimal fragmentation. Both mass spectrometers were equipped with a Z-spray nano-electrospray ionization source. Measurements were performed by using gold-coated needles, made from borosilicate glass capillaries (Kwik-Fill; World precision Instruments, Sarasota) on a P-97 puller (from Sutter Instruments, Novato, CA, United States). Needles were coated with a gold layer using an Edwards Scancoat six Pirani 501 sputter coater (Edwards laboratories, Milpitas, CA, United States). All TOF spectra were mass calibrated by using an aqueous solution of cesium iodide (25 mg⋅mL^-1).

Sequence Comparison

Protein sequences were retrieved using protein resources from the National Centre for Biotechnology Information² and UniProt Database³. Multiple sequence alignments were made using CLUSTALW (Thompson et al., 1994). Phylogenetic plots were made using FigTree⁴.

Results

Biochemical Properties of Rj3HB6H

Expression of the 3HB6H gene from R. jostii RHA1#2 yielded about 7 mg of purified Rj3HB6H protein from 10 g wet cells (Table 3A). Rj3HB6H displayed the same absorption spectrum as 3HB6H, with maxima at 274, 383, and 453 nm and a shoulder at 480 nm (Montersino and van Berkel, 2012). A molar absorption coefficient of protein-bound flavin, 𝜀₄₅₃ = 10.3 mM^-1 cm^-1, was used for both proteins.

Determination of steady-state kinetic parameters revealed that Rj3HB6H behaves similarly as 3HB6H using 3-hydroxybenzoate as variable substrate and fixed NADH concentration (k_cat = (20 ± 1) s^-1; K_M = (35 ± 3) μM; k_cat/K_M = (5.7 ± 0.8) × 10⁵ s^-1⋅M^-1) and with variable concentration of NADH (preferred coenzyme) and fixed 3-hydroxybenzoate concentration (k_cat = (20 ± 1) s^-1; K_M = (68 ± 5) μM; k_cat/K_M = (3.0 ± 0.4) × 10⁵ s^-1⋅M^-1). Rj3HB6H displays a very low NADH oxidase activity (<1 U⋅mg^-1). Uncoupling of hydroxylation of 3-hydroxybenzoate occurs to a minor extent (less than 10%), while 2,5-dihydroxybenzoate is a strong non-substrate effector (k_cat = (6 ± 0.8) s^-1; K_M = (150 ± 30) μM; k_cat/K_M = (4.0 ± 1.3) × 10⁴ s^-1⋅M^-1), efficiently stimulating the rate of flavin reduction by NADH (Montersino and van Berkel, 2012; Sucharitakul et al., 2012, 2013; Ni et al., 2016).

Structural Characterization

Rj3HB6H crystals grew in similar conditions as found for 3HB6H, and are isomorphous to those of 3HB6H, where lithium sulfate was present instead of sodium acetate. The three-dimensional structure of Rj3HB6H was solved at 2.3 Å resolution by molecular replacement (Table 4). The isoalloxazine moiety of FAD was refined with full occupancy in the in conformation. Similar to the crystallographic analysis of 3HB6H, no substrate could be detected in the active site of the enzyme, despite presence of excess 3-hydroxybenzoate in the crystallization drop. The protein crystallizes as a dimer, just as 3HB6H (Montersino et al., 2013), and contains a phospholipid molecule in each subunit. The electron density of the phospholipid in the crystal structure was refined as two acyl chains, one of twelve and one of seventeen carbon units. Superimposition of the Rj3HB6H and 3HB6H models (root mean square deviation = 0.22 Å) shows minor deviations (Figure 2A). The phospholipid is located in a tunnel, which runs from the dimer interface to the active site (Figure 2B), and interacts with the opposite monomer. The phosphate group resides at the protein surface near Arg350 and Lys385, and the electron density of the headgroup is consistent with the presence of a cyclohexanehexol moiety (Figure 3A).

FIGURE 2

FIGURE 3

Identification of Protein-Bound Lipid Molecules

Assignment of protein-bound phospholipids was achieved by ESI-MS analysis of the low molecular weight components extracted from denatured Rj3HB6H. The mass spectrum in negative mode (Figure 3B) displayed three main peaks with m/z values of 823, 837, and 851. From the MS pattern it was evident that Rj3HB6H binds a phospholipid with a bigger headgroup compared to that of the lipid found in 3HB6H.

Fragmentation analysis and comparison of data to reference lipid MS spectra led to a match of the obtained mass peaks with those of PI, having aliphatic chains containing 15 to 19 carbons. Peak 1 (m/z 823) is assigned to PI (15:0/18:0), peak 2 (m/z 837) is assigned to PI (16:0/18:0) (Sharp et al., 2007; Morita et al., 2011) and peak 3 ((m/z 851) is assigned to PI (16:0/19:0) with alternate acylate form (tuberculostearic acid) (Drage et al., 2010).

Typical fragmentation of PI was visible in the MS/MS spectra by signature peaks with m/z values of 153, 223, 241, and 297, representing glycerol phosphate water (m/z 153) and inositol headgroup fragments (Lang and Philp, 1998; Pulfer and Murphy, 2003; Oursel et al., 2007) (data not shown). Minor peaks at approximately 2 m/z values lower than the identified peaks represent the same PI, containing one unsaturated bond.

Protein Oligomeric Composition

To gain further insight into the enzyme–lipid interaction, we determined the oligomeric protein composition of 3HB6H and Rj3HB6H using native MS (Leney and Heck, 2017). As a first step, we determined the experimental masses of the denatured proteins. The measured values (46,766 ± 4 Da for 3HB6H and 46,761 ± 2 Da for Rj3HB6H) agree with the mass deduced from the primary sequence, lacking the N-terminal methionine.

Native MS of 3HB6H showed eight charge states corresponding to five different protein forms (Figure 4A and Table 5). The charge state distribution ions +12, +13, and +14 represent the monomeric apoprotein (average mass 46,835 ± 5 Da; red stars), the monomeric holoprotein (average mass 47603 ± 6 Da; green stars), and the monomeric holoprotein containing additionally one PG/PE molecule (average mass 48,312 ± 7 Da; blue stars). The charge state distribution ions +18, +19, +20, +21, and +22 predominantly represent the dimeric holoprotein with either one or two PG/PE molecules bound (average mass 95,868 ± 16 Da; orange stars, and 96,643 ± 14 Da; purple stars, respectively). Tandem MS experiments revealed that one, two, three, or four ligands can be expelled from 3HB6H. The assignment of bound ligands was made on the basis of total mass increase and comparison with the mass of the native apoprotein (Table 5).

FIGURE 4

Native MS of Rj3HB6H also showed a range of charge state distributions (Figure 4B and Table 5). The charge state distribution ions +12, +13, and +14 represent the monomeric apoprotein (average mass 46,829 ± 1 Da; red star), the monomeric holoprotein (average mass 47613 ± 1 Da; green stars), and the monomeric holoprotein containing one PI molecule (average mass 48,458 ± 2 Da; blue star). The latter species differs from the related 3HB6H species (Figure 4A, blue stars) because it has a bigger lipid headgroup. The major species in the native mass spectrum of Rj3HB6H corresponds to the holo-Rj3HB6H dimer with PI bound to both subunits (average mass 96,938 ± 32 Da; Figure 4B, yellow stars). Only by magnification it is possible to detect a minor peak representing the holo-Rj3HB6H dimer with one PI bound (average mass 96,128 ± 5 Da; Figure 4B, orange stars). A cartoon of the different subunit compositions of 3HB6H is presented in Figure 4D.

Conservation of Lipid Binding Site

To analyze whether the lipid-binding site of 3HB6H is conserved among species, we explored the natural diversity of 3HB6H enzymes. 3HB6H activity has been reported for Gram-positive and Gram-negative bacteria and for yeasts. Besides from the R. jostii prototype, the enzymes from Klebsiella pneumonia M5a1 (Suárez et al., 1995; Liu et al., 2005), Pseudomonas alcaligenes NCIMB 9867 (Gao et al., 2005), Polaromonas naphthalenivorans CJ2 (Park et al., 2007), Corynebacterium glutamicum ATCC 12032 (Yang et al., 2010), Rhodococcus sp. NCIMB 12038 (Liu et al., 2011) and Candida parapsilosis (Holesova et al., 2011) have been characterized to some extent.

From the structural data of the R. jostii 3HB6H enzyme and the multiple sequence alignment presented in Figure 5 it can be inferred that residues directly involved in lipid binding in Rj3HB6H are not always conserved in the orthologs; among the bacterial enzymes studied, most sequence divergence occurs in 3HB6H from P. alcaligenes NCIMB 9867 (Pa3HB6H). This prompted us to study the lipid binding properties of the Pseudomonas 3HB6H enzyme.

FIGURE 5

Expression of the Pa3HB6H gene in E. coli TOP10 cells yielded about 10 mg of enzyme from a 1 L batch culture. Purified Pa3HB6H had a specific activity of 34 U⋅mg^-1 (Table 3B) and migrated in SDS-PAGE as a single band with an apparent subunit mass of 47 kDa (not shown). ESI-MS established that native Pa3HB6H is a dimer, and not a trimer as suggested earlier (Gao et al., 2005), and that the enzyme indeed contains lipids (Figure 4C and Table 5). The mass spectrum of extracted lipids showed peaks with m/z values characteristic of PG and PE with aliphatic chains ranging from 14 to 19 carbons, similar to the previously identified lipids in 3HB6H from R. jostii RHA1 produced in E. coli (Montersino et al., 2013).

Discussion

3HB6H is a flavoenzyme that catalyzes the para-hydroxylation of 3-hydroxybenzoate to gentisate, a key step in the catabolism of lignin-derived aromatic compounds in the soil (Pérez-Pantoja et al., 2010). Up to now, 3HB6H is the only flavoprotein monooxygenase that has been found to bind phospholipids (Montersino et al., 2013). Structural analysis showed that the hydrophobic tails of the phospholipids deeply penetrate into the substrate-binding domains, whereas the hydrophilic parts are exposed on the protein surface, connecting the dimerization domains (Figure 1). Attempts to obtain native lipid-free protein were not successful, indicating that the phospholipids are important to attain a properly folded protein (Montersino et al., 2013).

3HB6H binds a mixture of PG and PE, the major constituents of the E. coli inner membrane (Montersino et al., 2013). By expressing its gene in R. jostii RHA1#2, we aimed at unraveling the lipid binding abilities of 3HB6H in the original host. Although E. coli gives considerable higher yields (Montersino and van Berkel, 2012), significant quantities of soluble His-tagged Rj3HB6H were obtained. The difference in enzyme yield could be linked to the type of induction and promoter strength used in the R. jostii RHA1#2 strain, which is based on the propane monooxygenase operon (Sharp et al., 2007). Nevertheless, our results show that the newly developed R. jostii RHA1#2 strain opens new prospects for actinomycetes as host cells for production of recombinant proteins (Nakashima et al., 2005).

Rj3HB6H displayed similar catalytic and structural properties as 3HB6H, and the mode of lipid binding was highly conserved (Figure 2). Gratifyingly, the crystallographic data and mass spectrometry analysis provided clear evidence that Rj3HB6H contains PI as natural glycerophospholipid cofactor (Figure 3). The crystal structure showed that the inositol headgroups of the phospholipids are located at the protein surface, and that the sn-2 acyl moieties are in contact with helix 11 of the other subunit (Figure 1). Based on MS/MS analysis, we identified the bound phospholipids as a mixture of PIs with carbon chains between 15 and 19 carbons. One of the extracted lipids was identified as tuberculostearic acid, an alternative acylated form of palmitate present in the membranes of Rhodococcus and Mycobacterium (Drage et al., 2010).

Rj3HB6H is a dimer both in solution and in crystal form, but native MS showed a ratio of monomer to dimer of about 1:3 (Figure 4B). Release of only one PI from the dimer resulted in monomerization in the gas phase. A similar observation was made with 3HB6H, but with this enzyme more dimers containing only one bound lipid (PG or PE) were detected (Figure 4A). 3HB6H dimers containing two phospholipids seem to be more stable in the gas phase than dimers containing one phospholipid. This strongly supports that lipid binding near the dimer interface stabilizes monomer contacts.

Native MS-analysis showed that Pa3HB6H is a homodimer and not a trimer as postulated earlier (Gao et al., 2005). The dimeric nature is in agreement with the structural properties of 3HB6H and Rj3HB6H. MS-analysis also revealed that recombinant Pa3HB6H binds the same type of phospholipids as 3HB6H. This supports that lipid binding is an intrinsic property of 3HB6Hs. As a main result, it appears that the 3HB6H family uses phospholipids as a common tool to increase their dimerization strength. Phospholipid binding is independent of the type of lipid headgroup, but relies on the presence of hydrophobic tunnels running from the protein surface to the active site.

Like PG/PE in E. coli (Oursel et al., 2007), PI is the major lipid membrane component in Rhodococcus (Nigou et al., 2003). This may explain why PG/PE are found as lipid ligands of 3HB6H and Pa3HB6H, while PI is found in Rj3HB6H. PI is the precursor for lipoarabinomannan and PI-mannoside synthesis. Glycolipid synthesis and reorganization of membrane composition allow Rhodococcus to adapt to environmental changes (Lang and Philp, 1998; Sharp et al., 2007; Guerin et al., 2010; Morita et al., 2011; De Carvalho et al., 2014). Binding of PI may localize 3HB6H at the cytoplasmic membrane, via inositol recognition of other proteins or specific phospholipid patching on the inner side of the membrane (Morita et al., 2011). At those specific spots, uptake of aromatic compounds from the environment may be coupled more efficiently to their catabolism.

Taking together, phospholipids do not have a direct catalytic role in 3HB6H, but are involved in stabilizing the dimer contact and, possibly, substrate orientation (Montersino et al., 2013). At this stage, we cannot exclude that bound phospholipids have some other function, for instance in directing the cytoplasmic membrane localization or in guiding/protecting molecules from entering the active site. In addition, the R. jostii RHA1#2 expression strain described here represents a useful alternative for the production of (whole-cell) biocatalysts.

References

• 1

  AdamsP. D.AfonineP. V.BunkócziG.ChenV. B.DavisI. W.EcholsN.et al (2010). PHENIX: a comprehensive Python-based system for macromolecular structure solution.Acta Crystallogr. D Biol.66(Pt 2)213–221. 10.1107/s0907444909052925

  • CrossRef
  • Google Scholar

• 2

  BrennanP. J.NikaidoH. (1995). The envelope of mycobacteria.Annu. Rev. Biochem.6429–63. 10.1146/Annurev.Bi.64.070195.000333

  • CrossRef
  • Google Scholar

• 3

  ChunJ.KangS. O.HahY. C.GoodfellowM. (1996). Phylogeny of mycolic acid-containing actinomycetes.J. Ind. Microbiol.17205–213. 10.1007/Bf01574694

  • CrossRef
  • Google Scholar

• 4

  De CarvalhoC. C.CostaS. S.FernandesP.CoutoI.ViveirosM. (2014). Membrane transport systems and the biodegradation potential and pathogenicity of genus Rhodococcus.Front. Physiol.5:133. 10.3389/fphys.2014.00133

  • CrossRef
  • Google Scholar

• 5

  DrageM. G.TsaiH. C.PecoraN. D.ChengT. Y.AridaA. R.ShuklaS.et al (2010). Mycobacterium tuberculosis lipoprotein LprG (Rv1411c) binds triacylated glycolipid agonists of Toll-like receptor 2.Nat. Struct. Mol. Biol.171088–1095. 10.1038/nsmb.1869

  • CrossRef
  • Google Scholar

• 6

  EmsleyP.LohkampB.ScottW. G.CowtanK. (2010). Features and development of Coot.Acta Crystallogr. D Biol.66486–501. 10.1107/S0907444910007493

  • CrossRef
  • Google Scholar

• 7

  EppinkM.SchreuderH.van BerkelW. (1997). Identification of a novel conserved sequence motif in flavoprotein hydroxylases with a putative dual function in FAD/NAD(P)H binding.Protein Sci.62454–2458. 10.1002/pro.5560061119

  • CrossRef
  • Google Scholar

• 8

  FinnertyW. R. (1992). The biology and genetics of the genus Rhodococcus.Annu. Rev. Microbiol.46193–218. 10.1146/Annurev.Micro.46.1.193

  • CrossRef
  • Google Scholar

• 9

  GaoX. L.TanC. L.YeoC. C.PohC. L. (2005). Molecular and biochemical characterization of the xlnD-encoded 3-hydroxybenzoate6-hydroxylase involved in the degradation of 2,5-xylenol via the gentisate pathway in Pseudomonas alcaligenes NCIMB 9867.J. Bacteriol.1877696–7702. 10.1128/JB.187.22.7696-7702.2005

  • CrossRef
  • Google Scholar

• 10

  GuerinM. E.KordulákováJ.AlzariP. M.BrennanP. J.JacksonM. (2010). Molecular basis of phosphatidyl-myo-inositol mannoside biosynthesis and regulation in mycobacteria.J. Biol. Chem.28533577–33583. 10.1074/jbc.R110.168328

  • CrossRef
  • Google Scholar

• 11

  GürtlerV.MayallB. C.SeviourR. (2004). Can whole genome analysis refine the taxonomy of the genus Rhodococcus?FEMS Microbiol. Rev.28377–403. 10.1016/j.femsre.2004.01.001

  • CrossRef
  • Google Scholar

• 12

  HolesovaZ.JakubkovaM.ZavadiakovaI.ZemanI.TomaskaL.NosekJ. (2011). Gentisate and 3-oxoadipate pathways in the yeast Candida parapsilosis: identification and functional analysis of the genes coding for 3-hydroxybenzoate 6-hydroxylase and 4-hydroxybenzoate 1-hydroxylase.Microbiology1572152–2163. 10.1099/mic.0.048215-0

  • CrossRef
  • Google Scholar

• 13

  HuijbersM. M. E.MontersinoS.WestphalA. H.TischlerD.van BerkelW. J. H. (2014). Flavin dependent monooxygenases.Arch. Biochem. Biophys.5442–17. 10.1016/j.abb.2013.12.005

  • CrossRef
  • Google Scholar

• 14

  LangS.PhilpJ. C. (1998). Surface-active lipids in rhodococci.Antonie van Leeuwenhoek7459–70. 10.1023/A:1001799711799

  • CrossRef
  • Google Scholar

• 15

  LeneyA. C.HeckA. J. R. (2017). Native mass spectrometry: what is in the name?J. Am. Soc. Mass Spectrom.285–13. 10.1007/s13361-016-1545-3

  • CrossRef
  • Google Scholar

• 16

  LiuD. Q.LiuH.GaoX. L.LeakD. J.ZhouN. Y. (2005). Arg(169) is essential for catalytic activity of 3-hydroxybenzoate 6-hydroxylase from Klebsiella pneumoniae M5a1.Microbiol. Res.16053–59. 10.1016/j.micres.2004.09.003

  • CrossRef
  • Google Scholar

• 17

  LiuT. T.XuY.LiuH.LuoS.YinY. J.LiuS. J.et al (2011). Functional characterization of a gene cluster involved in gentisate catabolism in Rhodococcus sp. strain NCIMB 12038.Appl. Microbiol. Biotechnol.90671–678. 10.1007/s00253-010-3033-1

  • CrossRef
  • Google Scholar

• 18

  McLeodM. P.WarrenR. L.HsiaoW. W. L.ArakiN.MyhreM.FernandesC.et al (2006). The complete genome of Rhodococcus sp RHA1 provides insights into a catabolic powerhouse.Proc. Natl. Acad. Sci. U.S.A.10315582–15587. 10.1073/pnas.0607048103

  • CrossRef
  • Google Scholar

• 19

  MontersinoS. (2012). Structural and Biochemical Characterization of 3-Hydroxybenzoate 6-Hydroxylase.Doctoral thesis, Wageningen UniversityWageningen.

  • Google Scholar

• 20

  MontersinoS.OrruR.BarendregtA.WestphalA. H.van DuijnE.MatteviA.et al (2013). Crystal structure of 3-hydroxybenzoate 6-hydroxylase uncovers lipid-assisted flavoprotein strategy for regioselective aromatic hydroxylation.J. Biol. Chem.28826235–26245. 10.1074/jbc.M113.479303

  • CrossRef
  • Google Scholar

• 21

  MontersinoS.TischlerD.GassnerG. T.van BerkelW. J. H. (2011). Catalytic and structural features of flavoprotein hydroxylases and epoxidases.Adv. Synth. Catal.3532301–2319. 10.1002/adsc.201100384

  • CrossRef
  • Google Scholar

• 22

  MontersinoS.van BerkelW. J. H. (2012). Functional annotation and characterization of 3-hydroxybenzoate 6-hydroxylase from Rhodococcus jostii RHAl.Biochim. Biophys. Acta1824433–442. 10.1016/j.bbapap.2011.12.003

  • CrossRef
  • Google Scholar

• 23

  MoritaY. S.FukudaT.SenaC. B. C.Yamaryo-BotteY.McConvilleM. J.KinoshitaT. (2011). Inositol lipid metabolism in mycobacteria: biosynthesis and regulatory mechanisms.Biochim. Biophys. Acta1810630–641. 10.1016/j.bbagen.2011.03.017

  • CrossRef
  • Google Scholar

• 24

  NakashimaN.MitaniY.TamuraT. (2005). Actinomycetes as host cells for production of recombinant proteins.Microb. Cell Fact.4:7. 10.1186/1475-2859-4-7

  • CrossRef
  • Google Scholar

• 25

  NiY.Fernandez-FueyoE.BaraibarA. G.UllrichR.HofrichterM.YanaseH.et al (2016). Peroxygenase-catalyzed oxyfunctionalization reactions promoted by the complete oxidation of methanol.Angew. Chem. Int. Ed. Engl.55798–801. 10.1002/anie.201507881

  • CrossRef
  • Google Scholar

• 26

  NigouJ.GilleronM.PuzoG. (2003). Lipoarabinomannans: from structure to biosynthesis.Biochimie85153–166. 10.1016/S0300-9084(03)00048-8

  • CrossRef
  • Google Scholar

• 27

  OurselD.Loutelier-BourhisC.OrangeN.ChevalierS.NorrisV.LangeC. M. (2007). Lipid composition of membranes of Escherichia coli by liquid chromatography/tandem mass spectrometry using negative electrospray ionization.Rapid Commun. Mass Spectrom.211721–1728. 10.1002/rcm.3013

  • CrossRef
  • Google Scholar

• 28

  ParkM.JeonY.JangH. H.RoH.-S.ParkW.MadsenE. L.et al (2007). Molecular and biochemical characterization of 3-hydroxybenzoate 6-hydroxylase from Polaromonas naphthalenivorans CJ2.Appl. Environ. Microbiol.735146–5152. 10.1128/AEM.00782-07

  • CrossRef
  • Google Scholar

• 29

  Pérez-PantojaD.GonzálezB.PieperD. H. (2010). “Aerobic Degradation of aromatic hydrocarbons,” inHandbook of Hydrocarbon and Lipid Microbiologyed.TimmisK. N. (Berlin: Springer) 799–837. 10.1007/978-3-540-77587-4_60

  • CrossRef
  • Google Scholar

• 30

  PottertonL.McNicholasS.KrissinelE.GruberJ.CowtanK.EmsleyP.et al (2004). Developments in the CCP4 molecular-graphics project.Acta Crystallogr. D Biol.602288–2294. 10.1107/S0907444904023716

  • CrossRef
  • Google Scholar

• 31

  PulferM.MurphyR. C. (2003). Electrospray mass spectrometry of phospholipids.Mass Spectrom. Rev.22332–364. 10.1002/mas.10061

  • CrossRef
  • Google Scholar

• 32

  RobertX.GouetP. (2014). Deciphering key features in protein structures with the new ENDscript server.Nucleic Acids Res.42W320–W324. 10.1093/nar/gku316

  • CrossRef
  • Google Scholar

• 33

  SchäferA.TauchA.JägerW.KalinowskiJ.ThierbachG.PühlerA. (1994). Small mobilizable multipurpose cloning vectors derived from the Escherichia coli plasmids pK18 and pK19: selection of defined deletions in the chromosome of Corynebacterium glutamicum.Gene14569–73. 10.1016/0378-1119(94)90324-7

  • CrossRef
  • Google Scholar

• 34

  SchrodingerL. L. C. (2015). ”The PyMOL Molecular Graphics System, Version 1.8”.

  • Google Scholar

• 35

  SharpJ. O.SalesC. M.LeBlancJ. C.LiuJ.WoodT. K.EltisL. D.et al (2007). An inducible propane monooxygenase is responsible for N-nitrosodimethylamine degradation by Rhodococcus sp strain RHA1.Appl. Environ. Microbiol.736930–6938. 10.1128/AEM.01697-07

  • CrossRef
  • Google Scholar

• 36

  SuárezM.FerrerE.GarridopertierraA.MartínM. (1995). Purification and characterization of the 3-hydroxybenzoate 6-hydroxylase from Klebsiella pneumoniae.FEMS Microbiol. Lett.126283–290. 10.1111/J.1574-6968.1995.Tb07431.X

  • CrossRef
  • Google Scholar

• 37

  SucharitakulJ.MedhanavynD.PakotipraphaD.van BerkelW. J.ChaiyenP. (2015). Tyr217 and His213 are important for substrate binding and hydroxylation of 3-hydroxybenzoate 6-hydroxylase from Rhodococcus jostii RHA1.FEBS J.283860–881. 10.1111/febs.13636

  • CrossRef
  • Google Scholar

• 38

  SucharitakulJ.TongsookC.PakotipraphaD.van BerkelW. J. H.ChaiyenP. (2013). The reaction kinetics of 3-hydroxybenzoate 6-hydroxylase from Rhodococcus jostii RHA1 provide an understanding of the para-hydroxylation enzyme catalytic cycle.J. Biol. Chem.28835210–35221. 10.1074/jbc.M113.515205

  • CrossRef
  • Google Scholar

• 39

  SucharitakulJ.WongnateT.MontersinoS.van BerkelW. J. H.ChaiyenP. (2012). Reduction kinetics of 3-hydroxybenzoate 6-hydroxylase from Rhodococcus jostii RHA1.Biochemistry514309–4321. 10.1021/bi201823c

  • CrossRef
  • Google Scholar

• 40

  SunJ. H.KelemenG. H.Fernández-AbalosJ. M.BibbM. J. (1999). Green fluorescent protein as a reporter for spatial and temporal gene expression in Streptomyces coelicolor A3(2).Microbiology1452221–2227. 10.1099/00221287-145-9-2221

  • CrossRef
  • Google Scholar

• 41

  SutcliffeI. C. (1998). Cell envelope composition and organisation in the genus Rhodococcus.Antonie van Leeuwenhoek7449–58. 10.1023/A:1001747726820

  • CrossRef
  • Google Scholar

• 42

  ThompsonJ. D.HigginsD. G.GibsonT. J. (1994). Clustal-W - improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice.Nucleic Acids Res.224673–4680. 10.1093/Nar/22.22.4673

  • CrossRef
  • Google Scholar

• 43

  van der GeizeR.DijkhuizenL. (2004). Harnessing the catabolic diversity of Rhodococci for environmental and biotechnological applications.Curr. Opin. Microbiol.7255–261. 10.1016/j.mib.2004.04.001

  • CrossRef
  • Google Scholar

• 44

  van der GeizeR.HesselsG. I.van GerwenR.van der MeijdenR.DijkhuizenL. (2002). Molecular and functional characterization of kshA and kshB, encoding two components of 3-ketosteroid 1-alpha-hydroxylase, a class IA monooxygenase, in Rhodococcus erythropolis strain SQ1.Mol. Microbiol.451007–1018. 10.1046/j.1365-2958.2002.03069.x

  • CrossRef
  • Google Scholar

• 45

  van der GeizeR.HesselsG. I.van GerwenR.VrijbloedJ. W.van der MeijdenP.DijkhuizenL. (2000). Targeted disruption of the kstD gene encoding a 3-ketosteroid delta(1)-dehydrogenase isoenzyme of Rhodococcus erythropolis strain SQ1.Appl. Environ. Microbiol.662029–2036. 10.1128/Aem.66.5.2029-2036.2000

  • CrossRef
  • Google Scholar

• 46

  WinnM. D.BallardC. C.CowtanK. D.DodsonE. J.EmsleyP.EvansP. R.et al (2011). Overview of the CCP4 suite and current developments.Acta Crystallogr. D Biol.67235–242. 10.1107/S0907444910045749

  • CrossRef
  • Google Scholar

• 47

  YamK. C.GeizeR.EltisL. D. (2010). “Catabolism of aromatic compounds and steroids by Rhodococcus,” inBiology of Rhodococcused.AlvarezM. H. (Berlin: Springer) 133–169.

  • Google Scholar

• 48

  YangY. F.ZhangJ. J.WangS. H.ZhouN. Y. (2010). Purification and characterization of the ncgl2923-encoded 3-hydroxybenzoate6-hydroxylase from Corynebacterium glutamicum.J. Basic Microbiol.50599–604. 10.1002/jobm.201000053

  • CrossRef
  • Google Scholar