crossmark THE JOURNAL OF BIOLOGICAL CHEMISTRY VOL. 290, NO. 39, pp. 23971–23983, September 25, 2015 © 2015 by The American Society for Biochemistry and Molecular Biology, Inc. Published in the U.S.A.

Structural Basis for the ATP-dependent Configuration of Adenylation Active Site in Bacillus subtilis o-Succinylbenzoyl-CoA Synthetase* Received for publication, July 2, 2015, and in revised form, August 8, 2015 Published, JBC Papers in Press, August 14, 2015, DOI 10.1074/jbc.M115.676304

Yaozong Chen, Yueru Sun, Haigang Song, and Zhihong Guo1 From the Department of Chemistry and State Key Laboratory for Molecular Neuroscience, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong SAR, China Background: It is not clear how the adenylation active site is formed in adenylating enzymes. Results: ATP binding induces adenylate-forming conformation, creates a binding site for the carboxylate substrate, and aligns active site residues for catalysis. Conclusion: ATP configures the adenylation active site. Significance: A new structural role is revealed for ATP in forming the active adenylation conformation in the catalysis of adenylating enzymes. o-Succinylbenzoyl-CoA synthetase, or MenE, is an essential adenylate-forming enzyme targeted for development of novel antibiotics in the menaquinone biosynthesis. Using its crystal structures in a ligand-free form or in complex with nucleotides, a conserved pattern is identified in the interaction between ATP and adenylating enzymes, including acyl/aryl-CoA synthetases, adenylation domains of nonribosomal peptide synthetases, and luciferases. It involves tight gripping interactions of the phosphate-binding loop (P-loop) with the ATP triphosphate moiety and an open-closed conformational change to form a compact adenylation active site. In MenE catalysis, this ATP-enzyme interaction creates a new binding site for the carboxylate substrate, allowing revelation of the determinants of substrate specificities and in-line alignment of the two substrates for backside nucleophilic substitution reaction by molecular modeling. In addition, the ATP-enzyme interaction is suggested to play a crucial catalytic role by mutation of the P-loop residues hydrogenbonded to ATP. Moreover, the ATP-enzyme interaction has also clarified the positioning and catalytic role of a conserved lysine residue in stabilization of the transition state. These findings provide new insights into the adenylation half-reaction in the domain alteration catalytic mechanism of the adenylateforming enzymes.

o-Succinylbenzoyl-CoA synthetase, or MenE, is one of the enzymes responsible for biosynthesis of vitamin K in bacteria, cyanobacteria, algae, and plants (1, 2). It converts OSB2 to OSB-

* This work was supported by Grants GRF601413 and N_HKUST621/13 from the Research Grants Council and Grant SBI14SC05 from the University Grants Council of the Government of the Hong Kong Special Administrative Region. The authors declare that they have no conflicts of interest with the contents of this article. The atomic coordinates and structure factors (codes 5BUQ, 5BUR, and 5BUS) have been deposited in the Protein Data Bank (http://wwpdb.org/). 1 To whom correspondence should be addressed. Tel.: 852-23587352; Fax: 852-23581594; E-mail: [email protected]. 2 The abbreviations used are: OSB, o-succinylbenzoyl; r.m.s.d., root mean square deviation; ACS, acyl/aryl-CoA synthetase; NRPS, nonribosomal peptide synthetase; PDB, Protein Data Bank; AMPCPP, ␣,␤-methyleneadeno-

SEPTEMBER 25, 2015 • VOLUME 290 • NUMBER 39

CoA and is essential in a number of bacteria that rely on vitamin K2, or menaquinone, for electron transportation in their respiratory chains, including microbial pathogens such as Haemophilus influenzae and Staphylococcus aureus (3–5). This enzyme is absent in mammals and thus is an attractive target for development of new antibacterial agents. In recent years, a number of potent inhibitors of MenE have been synthesized, and some of them have been shown to exhibit high antitubercular activities (6 – 8). MenE is homologous to acyl/aryl-CoA synthetases (ACS) and catalyzes a similar two-step reaction. In the first adenylation step, the OSB substrate reacts with Mg2⫹-ATP to form the adenylate intermediate, OSB-AMP, which subsequently reacts with the third substrate, coenzyme A, to form OSB-CoA in the second thioesterification step (Fig. 1). The same adenylation step is found in the reactions catalyzed by the adenylation domains of nonribosomal peptide synthetases (NRPS) and firefly luciferases, which together with acyl/aryl-CoA synthetases form the ANL family of adenylating enzymes (9). Over the last several decades, several members of this enzyme family, including MenE, have been kinetically analyzed to adopt an ordered Bi Uni Uni Bi ping-pong mechanism, which involves ATP as the first binding substrate and a large conformational change during the catalysis (6, 10 –15). Early structural and related biochemical studies of adenylating enzymes in the ANL family were focused on luciferases and the adenylation domains of NRPS, which were found to be highly conserved in the three-dimensional structure consisting of a large N-terminal domain (N-domain) containing ⬃400 – 550 residues and a smaller C-terminal domain (C-domain) containing ⬃130 residues (16 –18). Ten conserved motifs, namely A1–A10, were identified to play important structural and catalytic roles in the enzymatic reactions, of which the A3 motif is the most conserved among the ANL family (19). This conserved A3 motif overlaps with “motif I,” which defines a consensus

sine 5⬘-triphosphate; 4CBL, 4-chlorobenzoyl:coenzyme A ligase; LC-FACS, long chain fatty acyl-CoA synthetase.

JOURNAL OF BIOLOGICAL CHEMISTRY

23971

ATP Configures Adenylate-forming Conformation

FIGURE 1. Reaction catalyzed by OSB-CoA synthetase (MenE). A, MenE-catalyzed two-step reaction. B, ordered Bi Uni Uni Bi ping-pong catalytic mechanism.

sequence of (T/S)(S/G)G(T/S)(T/S/E)GX(P/Q/M/T)K (20) and is homologous to the conserved phosphate-binding loop (P-loop) involved in binding of triphosphate by ATP- and GTPbinding kinases (21). However, the active site formed at the interface of the N- and C-domains in the early crystal structures is only suitable for the first-step formation of the adenylate intermediate from the carboxylate and ATP substrates. In 2003, Gulick and co-workers (22) used a nonhydrolyzable adenylate, adenosine 5⬘-propyl phosphate, and coenzyme A to trap the bacterial acetyl-CoA synthetase in a different conformation in which the C-domain is rotated around a hinge residue by 140° to form a new active site suitable for the secondstep thioesterification. On this basis, a domain alteration mechanism was proposed for all adenylating enzymes in the ANL family involving two distinct active conformations, one for adenylate intermediate formation and the other for thioester formation (9). These two conformations form two distinct active sites with the same set of N-domain residues but a different set of residues on two opposite sides of the small C-domain, which rotates a large angle after adenylate formation to achieve conformational change and active site re-configuration. This domain alteration mechanism is strongly supported by the finding that mutation of residues on a different side of the C-domain specifically affects either the first or the second step of the reaction catalyzed by a number of adenylating enzymes such as 4-chlorobenzoyl:coenzyme A ligase (4CBL) (15), bacterial acetyl-CoA synthetase (23), luciferase (24), and the standalone adenylation domain EntE of the enterobactin synthetase (25). At present, both the adenylate-forming and thioesterforming conformations have been structurally characterized for several ANL family members, including DltA from Bacillus subtilis (26, 27), 4CBL (28, 29), human medium-chain acyl-CoA synthetase 2A (ACSM2A) (30), and firefly luciferase (16, 31). An interesting question in the domain alteration mechanism is how ATP configures the adenylation active site as the first binding substrate. At present, there is limited information about this process probably due to the difficulty to avoid autohydrolysis in co-crystallizing ATP with the enzymes. In the few available complex crystal structures (30, 32, 33), ATP or its analog is bound to the ANL enzymes in different modes and only the P-loop of ACSM2A is found to form gripping interactions with the triphosphate of the ligand as expected. In this study, we use MenE to further explore the binding and configuration of the adenylation active site by ATP. We have solved the crystal structures of B. subtilis MenE (BsMenE) in a new

23972 JOURNAL OF BIOLOGICAL CHEMISTRY

ligand-free form and in complex with ATP or AMP. These structures have led to a general ATP binding mode and an open-closed conformational change in ATP configuration of the adenylation active site.

Experimental Procedures Chemicals—The MenE substrate OSB was prepared chemoenzymatically from chorismic acid that was isolated from an engineered Escherichia coli strain KA12 (34). It was synthesized in one pot using recombinant EntC (35), MenC (36, 37), MenD (36, 37), and MenH (38 – 40) and purified by reverse-phase HPLC. The isolated product was lyophilized and re-dissolved in 10 mM Tris (pH 8.0) buffer for determination of its concentration by the UV absorbance at 250 nm (the extinction coefficient ⑀ ⫽ 4900 M⫺1 cm⫺1). Other chemicals used in this study, including ATP, AMP, CoA, NaHCO3, isopropyl ␤-D-1-thiogalactopyranoside, polyethylene glycol (PEG), (⫾)-2-methyl-2,4pentadiol, buffers, and other salts, were purchased from Sigma. Expression and Purification of BsMenE and Its Mutants— The menE gene was amplified from the genomic DNA of B. subtilis strain 168 using 5⬘-ATT TGC GGC CGC ATG CTG ACA GAA CAG CCC AA-3⬘ and 5⬘-CCG CTC GAG TCA TAG CAG TTC TCC TTT ACG-3⬘ as primers and was inserted into pET28a between NotI and BamHI sites. The subcloned gene was sequenced by Beijing Genomics Institute (BGI, Shenzhen, China) and found to be identical to the sequence deposited in GenBankTM (accession number NP_390957). The protein with a C-terminal hexahistidine tag was overexpressed in E. coli strain BL21 (DE3) in Luria broth (LB) containing 25 ␮g/ml kanamycin for 8 h at 37 °C after induction by 0.2 mM isopropyl ␤-D-1-thiogalactopyranoside and purified to homogeneity by Ni2⫹-chelating column chromatography and size exclusion chromatography (HiPrep 16/60 Sephacryl, GE Healthcare). The purified protein was concentrated, flash-frozen in liquid nitrogen, and stored at ⫺20 °C in 10 mM Tris-HCl (pH 8.0) buffer containing 2 mM 2-mercaptoethanol and 10% glycerol. The protein concentration was determined by PierceTM Coomassie Plus (Bradford) assay kit (Thermo Scientific). The expression constructs of BsMenE mutants R382A, R382K, T152A, G154P, T155A/T156A, T155A, T156A, and G157P were generated by QuikChangeTM site-directed mutagenesis kit (Stratagene) using the BsMenE-expressing plasmid in pET28a as template. All mutants were sequenced by BGI to confirm that only the desired mutations had been incorVOLUME 290 • NUMBER 39 • SEPTEMBER 25, 2015

ATP Configures Adenylate-forming Conformation porated into the menE gene. The expression and purification of these mutants followed the same protocol described above for the wild-type BsMenE. Activity Assay and Steady-state Kinetics—The BsMenE activity was determined by a coupling assay using E. coli MenB (6), which was obtained as an untagged recombinant protein as reported previously (41, 42). In the assay, the BsMenE product OSB-CoA was transformed by excess MenB to 1,4-dihydroxyl2-naphthoic acid-CoA for measurement at 392 nm (molar extinction coefficient ⑀ ⫽ 4000 M⫺1cm⫺1) (43). For the single substrate kinetics, a 140-␮l standard assay mixture contained 50 mM Tris-HCl (pH 7.5), 20 mM NaCl, 5 mM MgCl2, 80 ␮M MenB, 20 mM NaHCO3 and a varied concentration of the substrate OSB, ATP, or CoA with the other two substrates fixed at a saturation concentration (OSB at 1.0 mM, ATP and CoA at 2.0 mM). The reactions were initiated by addition of BsMenE to a final concentration of 250 nM after a 3-min incubation. The real time absorbance increase at 392 nm was measured for at least 1 min, and the reaction rate was calculated from the slope of the initial linear region of the reaction progress curve. All activity assays were carried out in triplicate at 24 °C. Subsequently, the data were fit to Michaelis-Menten Equation 1, v ⫽ V max关S兴/共KM ⫹ 关S兴兲

(Eq. 1)

where ␯ is the reaction rate (in ␮M/min); Vmax is the maximum reaction rate (in ␮M/min); [S] is the substrate concentration (in ␮M), and KM is the Michaelis-Menten constant (in ␮M). Crystallization—Crystallization conditions were screened using Index Screen and Crystal Screen I and II from Hampton Research or WizardTM Classic 3 and 4 from Rigaku. Crystallization was carried out with the hanging drop vapor diffusion method at 16 °C. A rod-like crystal (150 ⫻ 20 ⫻ 30 ␮m) of ligand-free BsMenE was obtained after a 10-day growth in a 1:1 mixture of a solution of BsMenE at 10 mg/ml and a reservoir solution at pH 6.15 containing 8.1% (w/v) PEG8000, 0.09 M cacodylate/HCl, 0.144 M calcium acetate, 12.6% glycerol, 0.049 M sodium phosphate monobasic, and 0.091 M dibasic potassium phosphate. A trapezoid-like crystal (200 ⫻ 100 ⫻ 70 ␮m) of BsMenE䡠ATP-Mg2⫹ complex was grown within a week in a 1:1 mixture of a solution of 10 mg/ml BsMenE preincubated with 3 mM ATP and 5 mM MgCl2 and a reservoir solution at pH 7.7 containing 7.2% (v/v) ethylene glycol, 9.9% (w/v) PEG8000, 0.09 M HEPES, 3.5% (v/v) 2-methyl-2,4-pentanediol, and 0.01 M acetate. Single hexagonal crystals of the BsMenE-AMP complex with the longest dimension up to 50 ␮m were grown in the 1:1 mixture of a solution of 10 mg/ml BsMenE preincubated with 3 mM AMP and a reservoir solution at pH 7.5 containing 0.10 M HEPES, 10% (w/v) PEG6000, 5% (v/v) (⫾)-2-methyl-2,4-pentadiol. These crystals were soaked in cryoprotectant containing 20% (v/v) glycerol in the mother liquor and then flash-frozen in liquid nitrogen before data collection. Data Collection, Processing, and Structural Determination— X-ray diffraction data were collected for the obtained BsMenE single crystals at the beamline BL17U of the Shanghai Synchrotron Radiation Facility with an ADSC Quantum 315R chargecoupled device detector. Diffraction images were indexed, integrated, and scaled using HKL2000 and iMOSFLM in the CCP4 SEPTEMBER 25, 2015 • VOLUME 290 • NUMBER 39

suite (44). TRUNCATE of the CCP4 package was used to convert the intensities to structure factors. The ligand-free BsMenE structure was solved by molecular replacement with Phaser (45) in the CCP4 suite using the backbone of the N-domain (residue 1–394) in chain A of S. aureus MenE (PDB code 3IPL) as the search model. The initial electron density map showed that two BsMenE polypeptide chains were correctly located in the asymmetric unit in the space group P1211 for the apo-form. After refinement by PHENIX (47), the C-domain(383– 486) of apo-MenE was extended via manual model building using Coot (48) followed by refinement with REFMAC (49) and PHENIX. When the refinement was close to completion, TLS anisotropic refinement (50) was performed for both the N- and C-domains. For the binary complex with AMP and ternary complex with ATP and Mg2⫹, the diffraction datasets were processed similarly, and chain A of the full-length apo-BsMenE was used as the search model for structural determination by molecular replacement. Both structures were determined to have a space group of P43 and refined using the same strategy as for the apo-protein. Restraints of ATP and AMP were generated and refined using eLBOW (51). The overall quality of the three structural models was assessed by PDB validation, Procheck (52), and MolProbity (53). Statistics for the data collection and structural refinement are summarized in Table 1. Structural Analysis and Sequence Alignment—BsMenE homologs were identified using a PSI-BLAST search and aligned with ClustalW 2.0 (54). Multiple sequence alignments were displayed with ESPript 3.0 (55). PyMOL version 1.3 was used in the structural analysis and generation of all graphics (56). PISA was used to analyze the protein interfaces (57), and CASTp was used to identify the cavities in the protein structures and to calculate their volume (58). Docking Simulation—All docking simulations were run with AutoDock Vina 1.1.2 using the ATP-BsMenE complex as searching macromolecule (59). The grid maps with 1.000 Å spacing were generated by AutoDock Tools (ADT) with the center set at the coordinates 54.113, 157.401, and 8.502 in the x, y, and z directions, respectively, and the grid box was set at a size of 14, 6, and 6 for the three dimensions. Hydrogen bonds and van der Waals interactions were modeled with a Lennard-Jones parameter of 12–10 and 12– 6, respectively. The electrostatic grid maps were calculated using the distance-dependent dielectric permittivity of Mehler and Solmajer (60). The docking simulation was run independently 10 times, and the resulting conformations within a root mean square deviation (r.m.s.d.) of 2.5 Å were clustered and analyzed.

Results Overall Structure—The crystal structures of ligand-free, AMP-bound, and ATP-bound BsMenE were determined and refined to a resolution of 1.98, 2.60, and 2.82 Å, respectively (Table 1). The structural models have good overall stereochemistry with most residues in the most favorable or allowed regions of the Ramachandran diagram. Ser-293 is the only Ramachandran outlier in every subunit of all three structures, which also takes the same unfavorable conformation in other BsMenE structures determined previously (61). Curiously, the JOURNAL OF BIOLOGICAL CHEMISTRY

23973

ATP Configures Adenylate-forming Conformation TABLE 1 Data collection and refinement statistics Values for the highest resolution shell are given in parentheses. apo-BsMenE PDB codes Data collection Space Group Unit cell dimensions a, b, c (Å) ␣, ␤, ␥ (°) Redundancy (%) Completeness (%) Reflections (unique) I/␴I Rmerge Refinement Resolution range (Å) No. of non-hydrogen Atoms Protein Water Ligands/ions Average B factor (Å2) Protein Water Ligand Rwork/Rfree r.m.s.d. for ideal values in Bond length (Å) Bond angle (degree) Ramachandran statisticsa a

ATP-BsMenE

AMP-BsMenE

5BUQ

5BUR

5BUS

P1211

P43

P43

82.8, 76.9, 93.9 90.0, 114.9, 90.0 2.3 98.83 (99.15) 141,919 (73,620) 22.1 (9.0) 0.048 (0.204)

121.8, 121.8, 97.8 90, 90, 90 6.2 99.74 (100.00) 606,644 (49,305) 30.0 (3.7) 0.091 (0.692)

121.8, 121.8, 98.0 90, 90, 90 3.1 99.46 (99.34) 2,042,714 (43,794) 6.6 (4.9) 0.136 (0.864)

33.27–1.98 8104 7064 1035 5 21.00 19.00 34.40 49.00 0.1605/0.1997

36.39–2.82 6269 6146 44 79 70.90 70.80 63.50 84.50 0.2160/0.2658

38.52–2.60 6500 6434 62 46 72.70 72.90 58.20 109.20 0.1904/0.2335

0.008 1.03 97.00/2.78/0.22%

0.011 1.35 93.00/6.64/0.36%

0.009 1.21 95.00/4.53/0.47%

Ramachandran statistics indicate the fraction of residues in the most favored, allowed, and disallowed regions of the Ramachandran diagram.

corresponding residue of BsMenE Ser-293 is also a Ramachandran outlier in the resolved crystal structures of other adenylate-forming enzymes (22). In the three new BsMenE structures, two substantially different polypeptide chains without the histidine tag and the first two residues at the N terminus are present in the asymmetric units. These two subunits form an interface with extensive hydrogen bonding and ionic interactions that are judged to be biologically relevant by PISA (Fig. 2A) (57). This dimeric assembly in the crystal structures is consistent with the homodimeric quaternary structure determined in solution by size exclusion column chromatography. Both BsMenE subunits in the ligand-free structure (apoBsMenE) contain a large N-domain (residues 3–380) and a small C-domain (residues 387– 486), whereas only one fulllength subunit (chain A) is found in the nucleotide-bound structures. The other subunit contains only the large N-domain with many disordered loops in both ATP-BsMenE and AMPBsMenE structures. A significantly higher number of residues are disordered in the linker region between the two domains and loops in the small C-domain in the ligand-free subunits than in the full-length nucleotide-bound subunits of the other two structures. Besides the nucleotide ligand, an ethylene glycol molecule from the crystallization buffer is found in a cavity next to ATP in the ATP-BsMenE structure. Interestingly, the C-domain takes a very different orientation relative to the N-domain to form a “closed” conformation in chain A and an “open” conformation in chain B in the ligand-free structure. By comparison, both nucleotide-bound BsMenE structures take a closed conformation similar to chain A of the apo-structure (Fig. 2B). Open bsMenE Conformations—The open form of the ligandfree BsMenE (chain B of the apo-structure structure, Fig. 2) is a new addition to the several open conformations found in both subunits of another apo-BsMenE structure in a different space group and the sulfate-bound BsMenE structure, which are dif-

23974 JOURNAL OF BIOLOGICAL CHEMISTRY

ferent in the distance between the small C-domain and the large N-domain (61). Although it is not possible to compare this new open apo-BsMenE structure with the previously determined structures due to unavailability of the latter in the Protein Data Bank, the multiplicity of the open conformations clearly shows the intrinsic property of the enzyme to be present in these structural forms in the absence of a ligand. Noticeably, a similar open conformation is also present for the apo-form of another MenE orthologue from S. aureus subsp. aureus Mu¨ 50 (PDB code 3IPL), demonstrating that the open conformation is not unique to BsMenE but may be a general structural form for all MenE orthologues. Actually, open conformations similar to those of BsMenE are also present for other adenylate-forming enzymes (16, 34). A common feature of these conformations is that the small C-domain rotates away from the N-domain at a large angle in comparison with the closed adenylate-forming conformation. However, the rotation angle varies in a wide range and is different for each enzyme. The abundance of these open conformational states clearly indicates the high conformational dynamics of the small C-domain and strongly suggests that it is an intrinsic property of the adenylate-forming enzymes to populate various open conformations with a variable distance between its two domains when substrates or other ligands are absent. In their resting state, these enzymes may also populate closed conformations as suggested by the closed conformations observed in the absence of a ligand for MenE (chain A of apo-BsMenE, Fig. 2B) and other adenylate-forming enzymes such as DhbE (18), 4CBL (29), and DltA (62). These resting conformations form the structural basis for a potential substrate recognition mechanism in which the active closed adenylate-forming conformation is either induced or stabilized from the resting open and closed conformations by binding of the carboxylate or ATP in the first-step adenylation reaction. VOLUME 290 • NUMBER 39 • SEPTEMBER 25, 2015

ATP Configures Adenylate-forming Conformation

FIGURE 2. Crystal structure of BsMenE. A, dimeric assembly of ligand-free BsMenE in its asymmetric unit. Chain A (green) is in a closed conformation, and chain B (brown) is in an open conformation. B, superposition of different BsMenE conformations. The open (brown) and closed (green) forms of ligand-free BsMenE are aligned with AMP-BsMenE (cyan) and ATP-BsMenE (gray) according to their N-domains (in C␣ lines) to demonstrate the different orientation of the C-domains (schematic). ATP of the ATP-BsMenE structure is presented as sticks to indicate the position of the interdomain active site. C, ATP bound to a tunnel with two ends open to the bulk solvent in the ATP-BsMenE structure. The complex is represented in surface and colored according to domains and conserved motifs as labeled. The arrow denotes a probable exit route for the pyrophosphate product after the adenylation reaction. ATP is shown as sticks with the coordinated Mg2⫹ ion denoted as a green sphere.

Closed BsMenE Structures—All three closed BsMenE subunits in the solved crystal structures are closely similar to the adenylate-forming conformations identified for a number of ACS, adenylation domains of NRPS, and firefly luciferases with their interdomain cavities corresponding to the identified active sites of the adenylation step. However, they exhibit a significant difference in the exact orientation of the C-domain (Fig. 2B), which rotates from the large N-domain in a slightly different angle around the linker region between the two domains. This difference is more pronounced between the ligand-free BsMenE and the nucleotide-bound structures with a calculated interdomain cavity of 6880 Å3 for the former and 6521 Å3 for the AMP-bound complex using CASTp (58), although it is much smaller between the two nucleotide-bound structures with the smallest interdomain cavity formed in the ATP-BsMenE structure. Notably, there is another closed BsMenE subunit in the previously determined crystal structure in complex with the OSB-AMP intermediate (61), which is not available for comparison. From the fact that the small C-domain takes two very different orientations in this OSBAMP-bound BsMenE structure as in the ligand-free BsMenE structure, it is likely that the interdomain cavity in the OSBAMP-bound closed conformation is similar to that in the ligand-free or AMP-bound BsMenE structure and larger than that in the ATP-BsMenE structure. The most compact ATP-BsMenE structure contains a ligand-binding tunnel with both ends extending to the protein surface (Fig. 2C), which consists of two segments crossed at an almost perpendicular angle that binds the nucleoside and SEPTEMBER 25, 2015 • VOLUME 290 • NUMBER 39

triphosphate moieties of the ATP ligand, respectively. In comparison, the binding tunnel is leaky in the AMP-BsMenE structure with an additional opening to the bulk solvent, while it is hardly formed in the ligand-free closed conformation due to the different C-domain orientation. The compact ATP-binding site is a strong indicator that the adenylation active site is properly formed only in the ATP-BsMenE structure, demonstrating the crucial role of the substrate in configuring the active site. In contrast, the closed apo-BsMenE and AMP-BsMenE structures are merely precursor conformations of this active conformation. In the ATP-bound structure, the buffer molecule ethylene glycol occupies a new pocket formed after ATP binding, which should be the binding site of the second substrate, OSB. Nucleotide-binding Site—The nucleotide-binding site of BsMenE is composed of many conserved amino acid residues mostly from the large N-domain, including the highly conserved nine-residue P-loop (152TSGTTGKPK160) in the A3 motif (amino acids 149 –160) (Fig. 3). As shown in Fig. 4A, the adenine ring of the ligand in the ATP-BsMenE structure makes direct hydrophobic contact with the side chain benzene ring of the conserved Tyr-286 in the A5 motif oriented at an almost perpendicular angle on one side and the backbone of the tripeptide Gly-263–Gly-264 –Gly-265 on the opposite side, although its 6-amino group forms a hydrogen bond with the backbone carbonyl group of Ser-285. The ATP ribose takes a 3⬘-endo pucker envelope and is stabilized by bifurcated hydrogen bonds between the side chain carboxylate of the invariant Asp-367 in the A7 motif (amino acids 365–367) and its 2⬘- and 3⬘-hydroxyl groups. The ␣-phosphate of ATP is bound by JOURNAL OF BIOLOGICAL CHEMISTRY

23975

ATP Configures Adenylate-forming Conformation

FIGURE 3. Alignment of sequences associated with substrate binding, structural changes, and transition state stabilization. Representative members with their PDB identification numbers in the ANL family and representative MenE orthologues are aligned separately with ClustalW 2.0, and only the relevant motifs are presented using ESPript 3.0. The BsMenE_5BUQ in the upper alignment is identical to the MenE orthologue from B. subtilis in the lower alignment.

hydrogen bonds with side chains of Thr-289 and Ser-153 at the N terminus of the P-loop, and likely also Lys-471 of the A10 motif (amino acids 465– 474) from the C-domain with its ⑀-NH3⫹ at a distance of 3.6 and 3.8 Å to one of the acidic oxygen atoms and the bridging oxygen between ␣- and ␤-phosphates, respectively. Interestingly, as shown in Fig. 4B, the ␤,␥-pyrophosphate of the ATP ligand is fixed to the position by a large number of polar interactions, including hydrogen-bonding interactions with five consecutive residues (152TSGTT156) of the P-loop, a salt bridge to the ␥-phosphate with the invariant Lys-160 at the P-loop C terminus, direct chelating interactions of ␤- and ␥-phosphates with the divalent Mg2⫹ ion anchored by the side chain carboxylate of the strictly conserved Glu-290 through mediation of two water molecules, and bidendate hydrogen bonds to the ␤-phosphate with the positive side chain of the invariant Arg-382 from the linker A8 motif. These tight wrapping interactions for the triphosphate moiety of ATP are reminiscent of the nucleotide-binding interaction of the P-loop in ATP- and GTP-dependent kinases and are likely necessary for precisely positioning the ␣-phosphate group for in-line nucleophilic reaction with the carboxylate substrate. AMP takes an orientation similar to ATP when bound to the nucleotide-binding site, sharing all the interactions for the adenine, the ribose, and the 5⬘-phosphate moieties but losing all the strong interactions for the ␤,␥-pyrophosphates of the latter (Fig. 4C). Consequently, the P-loop is significantly displaced from its position in the ATP-BsMenE structure, and the AMP phosphate group is dislocated by 0.5 Å in the AMP-BsMenE structure. The movement of the AMP phosphate allows it to form a 2.6 Å hydrogen bond with the imidazole ring of His-196 that takes different rotameric conformations in other closed structures. This histidine residue is equivalent to the A4 histidine residue conserved among long chain ACS, aryl-CoA synthetases, and luciferases, which corresponds to a phenylalanine residue in the adenylation domains of NRPS and a tryptophan

23976 JOURNAL OF BIOLOGICAL CHEMISTRY

residue in shorter chain ACS, and has been proposed to be involved in both the adenylation and thioesterification reactions. In the OSB-AMP-bound BsMenE structure, this histidine residue takes another slightly different rotamer conformation by forming a hydrogen bond with the oxygen atom of the carbonyl group linked to the phosphate group (61), suggesting that it may be involved in binding of the carboxylate substrate, OSB, and stabilizing the transition state of the adenylation reaction. Besides the different P-loop orientation and the associated different interactions, AMP-BsMenE is different from the ATP-bound structure in the absence of the ligand-mediated domain interaction involving Arg-382, which is apparently correlated to its incomplete closure of the interdomain active site cavity. In the ATP-BsMenE structure, the Arg-382 side chain strongly interacts with the ATP ␤-phosphate group via the bidentate hydrogen bonds with strong attractive electrostatic interaction and thus likely provides the driving force for the complete active site closure. This charged residue likely exerts this structural effect by affecting the backbone torsional angles of the downstream hinge residue Ser-384 to move the C-domain closer to the N-domain. Consistent with this proposed role, this arginine residue is strictly conserved among all ANL family members (Fig. 3). Comparison of ATP-bound Adenylate-forming Enzymes— Several crystal structures with bound ATP or its nonhydrolyzable analogues have been determined for human medium chain acyl-CoA synthase (ACSM2A) (30), D-alanyl carrier protein ligase DltA from Bacillus cereus (32), and long chain fatty acylCoA synthetase (LC-FACS) from Thermus thermophilus (33). The LC-FACS structure is unique in taking a thioester-forming conformation with the ligand bound to a site far from the P-loop, whereas all other structures are present in the adenylate-forming conformation with the ligand bound close to the P-loop. Among these structures, the ATP-bound BsMenE is most similar to ACSM2A bound to ATP or AMPCPP with a VOLUME 290 • NUMBER 39 • SEPTEMBER 25, 2015

ATP Configures Adenylate-forming Conformation

FIGURE 4. Stereo diagram of the BsMenE nucleotide-binding site. A, ATP in its binding site. ATP is represented as sticks (carbon atoms in green) and its mFo ⫺ DFc omit map is contoured at 2.5␴ in blue mesh, and the binding site residues are represented as sticks (carbon atoms in gray). B, binding interactions for the triphosphate moiety of ATP. This is magnified from a portion of the ATP-binding site in A and viewed in stereo from a slightly different perspective. C, comparison of AMP and ATP binding to BsMenE. The AMP-BsMenE (cyan) and ATP-BsMenE (gray) complexes are superimposed to reveal different P-loop orientation, different side chain orientation of the A8 Arg-382 and A10 Lys-471, and different protein-ligand interactions at the ␣-phosphate. The proteinligand interactions for the adenosine moiety of the ligands are similar in the complexes as shown in A but omitted in C for simplicity. In all three panels, the yellow broken lines denote hydrogen bonds with a distance ⱕ3.6 Å and the green spheres denote Mg2⫹; red dots denote bound water molecules.

similar C-domain orientation relative to the N-domain, despite a limited 19.6% sequence identity between the two proteins. In the ATP-ACSM2A structure (PDB code 3C5E) (30), most of the ligand-enzyme interactions are very similar to those in the ATP-BsMenE structure and involve the same number and the same type of interactions to the same atoms of the ligand (Fig. 5). The only obvious difference is that the invariant A10 Lys577 forms a hydrogen bond with the ␤-phosphate in ACSM2A SEPTEMBER 25, 2015 • VOLUME 290 • NUMBER 39

rather than a hydrogen bond to the ␣-phosphate of ATP by the corresponding A10 residue in BsMenE. Both the ATP-bound BsMenE and ACSM2A are significantly different from the ATP-DltA structure in the C-domain orientation by as much as ⬃48°, although they are all in the adenylate-forming conformation (32). Their nucleotide-binding sites also exhibit significant differences, although they share similar ligand-enzyme interactions for the adenine and ribose moieties JOURNAL OF BIOLOGICAL CHEMISTRY

23977

ATP Configures Adenylate-forming Conformation

FIGURE 5. Structural comparison of BsMenE and ACSM2A in ATP binding. The binding site residues from BsMenE (in gray sticks) are labeled, whereas their counterparts in ACSM2A (in yellow sticks) are presented without numbering. Hydrogen bonding interactions involving the catalytic A10 lysine (Lys-471 in BsMenE) and the A8 arginine (Arg-382 in BsMenE) are denoted by yellow broken lines.

and for the invariant A8 arginine (Arg-397 in DltA) and the Mg2⫹-chelating aspartate (Asp-298 in DltA). Specifically, there are few hydrogen bonds between the triphosphate and the P-loop in the ATP-DltA structure, in contrast to the tight wrapping of the triphosphate by the P-loop in BsMenE and ACSM2A. Actually, the P-loop is located at a distance from the triphosphate in DltA and a large part of it is disordered without contact with ATP. The reason for these different ATP interactions in DltA is unclear, and adverse crystal packing interactions, as implicated previously (32), are one of the possible explanations. Mutation of Residues in the P-loop—The P-loop is highly conserved among all members of the ANL family with a consensus sequence of (T/S)(S/G)G(T/S)(T/S/E)GX(P/Q/M/T)K (residues Thr-152 to Lys-160 in BsMenE) (20). Previous functional characterization of the conserved residues mainly focused on the penultimate proline residue and the strictly conserved glycine and lysine residues. It was found that introduction of hydrophobic isoleucine or negatively charged aspartate to the positions of glycine residues causes more than a 4 orders of magnitude decrease in catalytic efficiency in either 4CBL from Pseudomonas sp. strain CBS3 (20) or adenylation domains in NRPS (19). In contrast, changes at the invariant lysine to methionine or other polar residues are tolerated with a moderate decrease in catalytic activity (20, 63, 64), whereas mutations at the proline residue cause enzyme-dependent effects that range from no significant activity change for a proline-to-valine change in an NRPS adenylation domain (65) to complete inactivation for a proline-to-alanine mutation in 4CBL (20). From the ATP-enzyme interactions in BsMenE and ACSM2A (30), it is easy to understand that the two glycine residues in the P-loop have to be conserved to avoid steric repulsion with the bound ATP, whereas the proline and lysine residues are important in forming the ␤-turn in the P-loop hairpin to tightly wrap around the triphosphate of the ligand, for which the former directly facilitates the directional change of the polypeptide chain, and the latter stabilizes the turn by forming hydrogen-bonding interactions with the ligand as well as other P-loop residues. To better understand the catalytic contribution of the intense enzyme-ATP interaction (Fig. 4, A and B), the conserved serine/threonine residues in the P-loop interacting with the ␤,␥pyrophosphate of the substrate were mutated into alanine individually or in combination. For comparison, the invariant A8

23978 JOURNAL OF BIOLOGICAL CHEMISTRY

Arg-382 hydrogen-bonded to the ␤-phosphate was mutated to alanine or lysine. In addition, proline was also introduced into the positions of the conserved glycine residues. The resulting mutants, namely T152A, T155A, T156A, T155A/T156A, R382A, R382K, G154P, and G157P, were successfully obtained in high purity with the same stability and circular dichroism spectra as the wild-type enzyme (data not shown), indicating the absence of major structural change caused by the mutations. The single-substrate kinetic constants kcat, KM, and kcat/KM of the mutants and the wild-type enzyme are listed in Table 2. It is interesting to note that the kinetic constants are moderately affected for the proline mutants of the conserved glycine residues, leading to a 20 –95-fold decrease for ATP and less than a 15-fold decrease for other substrates in kcat/KM in comparison with the wild-type enzyme. These moderate changes indicate that the introduced proline causes a much smaller steric repulsion in comparison with isoleucine in 4CBL (20), and it has only a mild effect on the P-loop backbone orientation, which likely depends more on the multiple hydrogen bonds between ATP and the P-loop residues. In contrast, the threonine-to-alanine mutations cause a much larger decrease in catalytic activity. The resulting mutants exhibit a general kcat decrease of about 2 orders of magnitude for all substrates and a significant substrate-dependent KM increase, which is 100 –200-fold for ATP, 2–7-fold for OSB, and 1–6.8-fold for CoA-SH. Overall, the decrease in the catalytic efficiency (kcat/KM) is 4800 –24,000-fold for ATP, 180 –1700-fold for OSB, and 62–190-fold for CoA-SH in single mutants and significantly higher for the double mutant T155A/T156A with a 81,000-fold decrease for ATP and a 1700fold decrease for the other two substrates. In comparison, the R382A and R382K mutations cause similar kcat and kcat/KM decreases for both ATP and OSB with an extent comparable with that of the T156A mutant (Table 2). However, these two mutants cause a higher KM increase for the coenzyme A substrate than T156A, suggesting that Arg-382 plays a critical role in both the adenylation and the thioesterification step. In contrast, all the P-loop mutants mainly affect the first-step reaction. It is interesting to note that kcat is significantly decreased by about 2 orders of magnitude when one or two hydrogen bonds are removed by the mutations from the numerous binding interactions for the ␤,␥-pyrophosphate of the ATP substrate VOLUME 290 • NUMBER 39 • SEPTEMBER 25, 2015

ATP Configures Adenylate-forming Conformation TABLE 2 The single substrate kinetic constants of wild-type and mutant BsMenE Protein

Substrate

kcat min

WT T152A G154P T155A T156A T155A/T156A G157P R382A R382K

OSB ATP CoA OSB ATP CoA OSB ATP CoA OSB ATP CoA OSB ATP CoA OSB ATP CoA OSB ATP CoA OSB ATP CoA OSB ATP CoA

kcat(WT)/kcat

KM

⫺1

(7.2 ⫾ 0.3) ⫻ 102 (6.3 ⫾ 0.2) ⫻ 102 (7.3 ⫾ 0.8) ⫻ 102 8.2 ⫾ 0.7 8.1 ⫾ 1.2 12.4 ⫾ 1.5 (1.3 ⫾ 0.1) ⫻ 102 99 ⫾ 14 (1.3 ⫾ 0.1) ⫻ 102 7.5 ⫾ 0.9 4.6 ⫾ 2.5 6.2 ⫾ 0.6 12.8 ⫾ 0.3 13.1 ⫾ 0.6 8.7 ⫾ 0.1 2.7 ⫾ 0.2 1.1 ⫾ 2.0 3.0 ⫾ 0.6 (2.8 ⫾ 0.2) ⫻ 102 (2.4 ⫾ 0.3) ⫻ 102 (2.6 ⫾ 0.2) ⫻ 102 9.0 ⫾ 0.2 13.1 ⫾ 2.6 17 ⫾ 1 13.7 ⫾ 0.8 9.8 ⫾ 0.5 9.7 ⫾ 4.8

␮M

1 1 1 88 78 59 5.5 6.3 5.6 96 137 118 56 48 84 267 573 243 2.6 27 2.8 80 48 43 53 64 75

44 ⫾ 19 24 ⫾ 4 (2.4 ⫾ 0.8) ⫻ 102 (8.8 ⫾ 1.1) ⫻ 102 (3.5 ⫾ 1.2) ⫻ 103 (2.6 ⫾ 0.7) ⫻ 102 (1.2 ⫾ 0.2) ⫻ 102 (3.4 ⫾ 1.0) ⫻ 102 (2. 7 ⫾ 1.0) ⫻ 102 87 ⫾ 40 (4.2 ⫾ 3.9) ⫻ 103 (4.1 ⫾ 0.8) ⫻ 102 (1.4 ⫾ 0.2) ⫻ 102 (2.3 ⫾ 0.2) ⫻ 103 (2.1 ⫾ 0.3) ⫻ 102 (2.8 ⫾ 0.5) ⫻ 102 (3.3 ⫾ 1.1) ⫻ 103 (1.7 ⫾ 0.9) ⫻ 103 80 ⫾ 11 (1.9 ⫾ 0.8) ⫻ 102 (2.5 ⫾ 1.2) ⫻ 102 57 ⫾ 10 (5.1 ⫾ 1.8) ⫻ 102 (2.7 ⫾ 0.8) ⫻ 103 66 ⫾ 19 (2.8 ⫾ 0.2) ⫻ 102 (3.0 ⫾ 2.3) ⫻ 103

(Fig. 4B), demonstrating that the conserved residues not only are involved in ATP binding but also affect subsequent adenylation reaction in the catalytic process. These residues may contribute to the chemical reaction by keeping the ATP substrate in optimal position and orientation for the reaction and by stabilizing the transition state through their interactions with the leaving ␤-phosphate group. In addition, these residues together with other active site residues and the metal cofactor may also be involved in activating the scissile phosphorus– oxygen bond by straining it via the observed tight binding of the ATP substrate (Fig. 4A). Modeling of the OSB Substrate into Its Binding Site—The OSB substrate was successfully modeled by AutoDock Vina into the binding site adjacent to the bound ATP in its complex with BsMenE after removing the buffer molecule, ethylene glycol. As shown in Fig. 6, the final model was chosen from five structures resulting from the docking simulations on the basis of its higher similarity to the OSB-AMP-bound BsMenE structure (61) in the interactions between the OSB moiety and the enzyme. As in the adenylate-bound BsMenE structure, OSB is bound between a loop containing the tripeptide 261LLG263 with the hydrophobic leucine side chains interacting with the benzene ring of the ligand and the loop containing the A4 motif (196HISG199). All these binding site residues are specifically conserved among the MenE family members but not among different adenylating enzymes (Fig. 3). Other residues lining the substrate-binding site include Ser-293, Gln-294, Ser-237, Gly287, Lys-471, and the bound ATP. Despite their similarities, the docked structure and the adenylate-bound structure exhibit marked differences. Notably, OSB in the docked structure takes a conformation with its benzene ring twisted at an angle to both the 2-carboxylate and the carbonyl function of the succinyl group, whereas the OSB moiety in the adenylate-bound structure takes a different conSEPTEMBER 25, 2015 • VOLUME 290 • NUMBER 39

KM/KM(WT)

kcat/KM min

1 1 1 20 146 1.1 2.7 14 1.1 2.0 175 1.7 3.2 96 0.9 6.4 137 7.1 1.8 7.9 1.0 1.3 21 11 1.5 12 12

⫺1

kcat(WT) KM/kcat KM(WT)

䡠M⫺1

(1.6 ⫾ 0.6) ⫻ 107 (2.7 ⫾ 0.4) ⫻ 107 (3.0 ⫾ 0.6) ⫻ 106 (9.4 ⫾ 0.6) ⫻ 103 (2.3 ⫾ 0.4) ⫻ 103 (4.8 ⫾ 0.8) ⫻ 104 (1.1 ⫾ 0.1) ⫻ 106 (2.9 ⫾ 0.6) ⫻ 105 (4.8 ⫾ 1.7) ⫻ 105 (8.6 ⫾ 2.5) ⫻ 104 (1.1 ⫾ 0.6) ⫻ 103 (1.5 ⫾ 0.2) ⫻ 104 (9.2 ⫾ 1.3) ⫻ 104 (5.6 ⫾ 0.3) ⫻ 103 (4.1 ⫾ 0.7) ⫻ 104 (9.6 ⫾ 0.8) ⫻ 103 (3.3 ⫾ 1.4) ⫻ 102 (1.8 ⫾ 1.4) ⫻ 103 (3.5 ⫾ 0.4) ⫻ 106 (1.3 ⫾ 0.3) ⫻ 106 (1.0 ⫾ 0.4) ⫻ 106 (1.6 ⫾ 0.2) ⫻ 105 (2.6 ⫾ 1.4) ⫻ 104 (6.3 ⫾ 1.2) ⫻ 103 (2.1 ⫾ 0.4) ⫻ 105 (3.5 ⫾ 0.4) ⫻ 104 (3.2 ⫾ 0.9) ⫻ 103

1 1 1 1.7 ⫻ 103 1.2 ⫻ 104 62 15 92 6.2 1.9 ⫻ 102 2.4 ⫻ 104 2.0 ⫻ 102 1.8 ⫻ 102 4.8 ⫻ 103 70 1.7 ⫻ 103 8.1 ⫻ 104 1.7 ⫻ 103 4.6 20 3.0 1.0 ⫻ 102 1.0 ⫻ 103 4.7 ⫻ 102 79 7.7 ⫻ 102 9.3 ⫻ 102

formation with its benzene ring perpendicular to the 2-carboxylate and coplanar with the carbonyl function of the succinyl group. In addition, the conformation of the ethylene moiety of the ligand is gauche in the former and anti in the latter. Other differences include the absence of hydrogen bonding of the His196 side chain imidazole and the Ser-293 side chain hydroxyl to the ligand in the modeled structure, which is present in the adenylate-bound structure due to rotameric flipping of the side chain functional groups and is believed to occur in the binding of OSB by the ATP-bound BsMenE. In the modeled complex, one oxygen atom of the succinyl carboxylate is close to the ␣-phosphate of ATP with a short oxygen-to-phosphorus distance of 3.5 Å and is opposite the pyrophosphate group of ATP. This positioning of OSB allows backside nucleophilic attack on the ATP ␣-phosphate by its carboxylate from the other side of the pyrophosphate leaving group in an almost linear manner, which is consistent with the in-line reaction mechanism for the adenylation half-reaction proposed for enzymes in the ANL family. Interestingly, the ⑀amino group of the invariant A10 Lys-471, which has been known to interact with the ␣-phosphate of ATP (Fig. 3B), is found to be within a hydrogen-bonding distance from the other oxygen atom on the attacking carboxylate group in OSB. Fully consistent with its proposed catalytic role, the simultaneous interaction of this invariant residue with both reacting functional groups gives it the best position to stabilize the transition state.

Discussion Adenylating enzymes in the ANL family catalyze diverse chemical transformations in cellular metabolism with two distinct active sites configured by domain alteration, which is consistent with an ordered Bi Uni Uni Bi ping-pong kinetic mechanism reported for various enzymes in the ACS subfamily such JOURNAL OF BIOLOGICAL CHEMISTRY

23979

ATP Configures Adenylate-forming Conformation

FIGURE 6. ATP and the docked OSB in their binding sites in BsMenE. ATP and OSB are represented in green and yellow sticks for their carbon atoms, respectively. Residues contacting OSB are presented in gray sticks for their carbon atoms, and those interacting with ATP are omitted. The yellow and blue broken lines denote hydrogen bonds and the shortest distance (3.5 Å) between the OSB succinyl carboxylate and ␣-phosphate of ATP, respectively.

as MenE (6), 4CBL (15), and ACS with various chain length specificities (11–14). Binding of ATP is the critical first step in the complex catalytic mechanism of these enzymes but is not fully understood. Reported structures of three adenylating enzymes complexed with ATP or its analogues have presented three different pictures for the important ATP-enzyme interaction. In this work, we have successfully determined the structure of ATP-bound BsMenE with its P-loop tightly wrapping around the triphosphate as expected for the same structural motifs in ATP- or GTP-dependent kinases (21). This enzyme is found to be highly similar to human ACSM2A in binding ATP with all but one of the numerous ligand-enzyme interactions conserved between the evolutionarily distal proteins. Remarkably, all the amino acid residues involved in the ATP-enzyme interaction are highly or even strictly conserved, suggesting an identical ATP-binding mode and the same functional and structural effects for ATP binding in all ANL family members. Thus, together with the previous findings for the interaction of ATP with human ACSM2A (30), this ATP-BsMenE complex structure provides a highly conserved structural model for ATP interaction with ANL enzymes to form the active adenylateforming conformation. Currently, it is not clear why the nucleotide substrate is bound to a site far from the P-loop in LCFACS (33). An important component of this conserved ATP-enzyme interaction is the ligand-dependent open-closed conformational change that emerges from current and numerous previous studies. For BsMenE, co-existence of both open and closed conformations in the ligand-free enzyme (Fig. 2), together with the open structures previously determined (61), is a clear indication that the enzyme populates both the open and closed conformations in its resting state. Considering the tight gripping interactions for the ATP triphosphate that makes a 90° turn in the binding tunnel, the compact ATP-BsMenE complex is much more likely to form from an open conformation rather than a closed conformation similar to the active adenylateforming conformation. On the basis of the fact that open con-

23980 JOURNAL OF BIOLOGICAL CHEMISTRY

formations similar to those of BsMenE have been found in the absence of a ligand for other adenylating enzymes in the ANL family (16, 34), such an open-closed conformational change is proposed to generally occur as the first step to form the adenylate-forming active site in catalysis of all ANL enzymes. Actually, a similar open-closed conformational change has been proposed for the formation of the adenylation active site in the catalysis of DltA, although no open conformations have been found for the ligand-free enzyme (26). The slight C-domain orientation difference between AMPand ATP-bound complexes (Fig. 2B) indicates that there exist structural elements safeguarding the formation of the compact adenylation active site in the open-closed conformational change. Previously, the strictly conserved A10 lysine residue was proposed to play such a role on the basis of its hydrogen bonding to the ATP ␤-phosphate in ACSM2A (30). However, the A10 Lys-471 in BsMenE makes the same hydrogen-bonding contact with the ␣-phosphate of the ligands as in other adenylating enzymes (17, 18) and is unable to account for the slight structural difference between the two nucleotide-bound complexes. In contrast, the invariant A8 Arg-382 forms bidentate hydrogen bonding with the ATP ␤-phosphate in the ATPBsMenE complex but makes no contact with the ligand in the AMP-BsMenE complex. This is the only residue that makes differential interactions in the two complexes, which may impact the C-domain orientation by affecting the torsional angles of the downstream hinge residue (Ser-384 in BsMenE) in the linker region between the N- and C-domains. This invariant A8 residue is therefore believed to be crucial for proper formation of the adenylation active site, which is consistent with its pronounced effect on the enzyme activity (Table 2). Other residues at the domain interface are involved in the same interactions in the two nucleotide-bound complexes. Another important part of the ATP configuration of the adenylation active site is to create a new binding site for the carboxylate substrate. This new binding pocket is successfully defined in BsMenE by molecular docking with the OSB subVOLUME 290 • NUMBER 39 • SEPTEMBER 25, 2015

ATP Configures Adenylate-forming Conformation strate and comparison with the OSB-AMP-bound BsMenE structure (Fig. 6), which is composed of amino acid residues specifically conserved among the MenE family (Fig. 3). Importantly, the two substrates are perfectly positioned for in-line backside nucleophilic substitution reaction with a short oxygen-to-phosphorus distance between the reacting functional groups (Fig. 6). This allows us a glimpse of how the enzyme sequentially binds the two substrates and aligns them just before the reaction happens. It also provides a general model for alignment of the two substrates for the first half-adenylation reaction catalyzed by other ANL enzymes. In addition to creating a new substrate binding, the ATP binding also positions conserved amino acid residues for stabilization of the transition state. Noticeably, the ATP-bound BsMenE shows for the first time hydrogen bonding interaction of the invariant A10 Lys-471 with the ␣-phosphoryl group of the ATP substrate, strongly supporting the proposed stabilization of the developing negative charge in forming the transition state of the first step adenylation (32). The fact that this strictly conserved lysine interacts with the ␤-phosphate of ATP in DltA and ACSM2A shows that its side chain is highly dynamic and flexible. The crucial catalytic role of the invariant A10 lysine residue is strongly supported by the complete activity loss in the adenylation half-reaction without affecting the thioesterification half-reaction in several ANL enzymes when it is mutated (8, 13, 20, 66, 67). It is also supported by the hydrogen bonding interaction of this residue with the ␣-phosphoryl group in complexes of the adenylate-AMP intermediate or analogues with several ANL enzymes, including BsMenE (17, 18, 61, 68) and the use of acetylation of this residue as a regulation control of acetyl-CoA synthetase in mammals (46). Noticeably, the hydrogen bond between this invariant A10 lysine residue and the OSB carboxylate group found in the modeled ATP-Mg2⫹䡠 OSB䡠-BsMenE complex should align the substrates for the reaction and put the side chain in the proper position to stabilize the transition state, although it may reduce the nucleophilicity of the attacking carboxylate. With the transition state stabilization by the invariant A10 lysine, the reaction between the two aligned substrates as visualized in the modeled ATP-OSBBsMenE complex (Fig. 6), should proceed smoothly to form the adenylate intermediate and pyrophosphate. Noticeably, a structure relaxation mechanism is intrinsically installed in the compact ATP-configured adenylate-forming conformation for its easy disassembly, which is required after adenylation to allow transition to the thioester-forming conformation in the domain alteration mechanism. When the adenylation is complete in the adenylate-forming conformation, the pyrophosphate product is released easily from one end of the ATP-binding tunnel, which is very close to the protein surface (Fig. 2C), to conform to the ordered Bi Uni Uni Bi ping-pong mechanism. This product release eliminates the strong interaction of the A8 Arg-382 to the ␤-phosphate of ATP, which controls the correct orientation of the C-domain in the conserved ATP-dependent configuration of the adenylation active site. Consequently, the interdomain interaction is significantly weakened, similar to that observed for the loosened AMPBsMenE structure (Fig. 2B), to reach a state that allows another conformational change to accept the coenzyme A substrate and SEPTEMBER 25, 2015 • VOLUME 290 • NUMBER 39

form a new active site for the second-step thioesterification. Regaining the C-domain flexibility after adenylation is demonstrated by the two different C-domain orientations found in the OSB-AMP-BsMenE structure (61). In summary, besides being the first substrate of the ANL enzymes, ATP plays a pivotal structural role in configuring the active adenylate-forming conformation to impact on the adenylation half-reaction in the ordered substrate binding, the in-line alignment of the two substrates for backside nucleophilic substitution, and stabilization of the transition state. This is achieved by its intense interaction with the P-loop and other structural elements, which are highly conserved in the ANL family. For this conserved mode of ATP binding in configuration of the adenylation active site, it is important to note that the strongest interaction affecting the C-domain orientation occurs at the ␤-phosphate of the bound ATP. This mode of C-domain control allows easy disassembly of the adenylateforming conformation to release the C-domain for a large angle rotation in configuration of the second active site for thioesterification, once the first half-reaction releases its pyrophosphate product to eliminate the domain-affecting interaction. Consequently, the ATP configuration of the adenylate-forming conformation is transient and is intimately coupled to progression of the adenylation half-reaction. Revelation of this structural role of ATP has provided novel mechanistic insights into the adenylation half-reaction and the transition to the thioesterification half-reaction to allow us a better understanding of the domain alteration mechanism of the adenylating enzymes in the ANL family. Author Contributions—Z. G. conceived and coordinated the study. Z. G. and Y. C. wrote the paper. Y. C. designed, performed, and analyzed experiments in all figures. Y. S. and H. S. participated in protein purification, protein crystallization, data collection, and structural determination. All authors reviewed the results and approved the final version of the manuscript. Acknowledgments—We are grateful for access to beamline BL17U1 at Shanghai Synchrotron Radiation Facility (SSRF), and we thank the beamline staff for technical support. References 1. Nowicka, B., and Kruk, J. (2010) Occurrence, biosynthesis and function of isoprenoid quinones. Biochim. Biophys. Acta 1797, 1587–1605 2. Gross, J., Meurer, J., and Bhattacharya, D. (2008) Evidence of a chimeric genome in the cyanobacterial ancestor of plastids. BMC Evol. Biol. 8, 117 3. Akerley, B. J., Rubin, E. J., Novick, V. L., Amaya, K., Judson, N., and Mekalanos, J. J. (2002) A genome-scale analysis for identification of genes required for growth or survival of Haemophilus influenzae. Proc. Natl. Acad. Sci. U.S.A. 99, 966 –971 4. Forsyth, R. A., Haselbeck, R. J., Ohlsen, K. L., Yamamoto, R. T., Xu, H., Trawick, J. D., Wall, D., Wang, L., Brown-Driver, V., Froelich, J. M., C, K. G., King, P., McCarthy, M., Malone, C., Misiner, B., et al. (2002) A genome-wide strategy for the identification of essential genes in Staphylococcus aureus. Mol. Microbiol. 43, 1387–1400 5. Kobayashi, K., Ehrlich S. D., Albertini, A., Amati, G., Andersen, K. K., Arnaud, M., Asai, K., Ashikaga, S., Aymerich, S., Bessieres, P., Boland, F., Brignell, S. C., Bron, S., Bunai, K., Chapuis, J., et al. (2003) Essential Bacillus subtilis genes. Proc. Natl. Acad. Sci. U.S.A. 100, 4678 – 4683 6. Tian, Y., Suk, D.-H., Cai, F., Crich, D., and Mesecar, A. D. (2008) Bacillus anthracis o-succinylbenzoyl-CoA synthetase: reaction kinetics and a

JOURNAL OF BIOLOGICAL CHEMISTRY

23981

ATP Configures Adenylate-forming Conformation

7.

8.

9.

10. 11.

12.

13.

14.

15.

16.

17.

18.

19.

20.

21.

22.

23.

24.

25.

26.

novel inhibitor mimicking its reaction intermediate. Biochemistry 47, 12434 –12447 Lu, X., Zhang, H., Tonge, P. J., and Tan, D. S. (2008) Mechanism-based inhibitors of MenE, an acyl-CoA synthetase involved in bacterial menaquinone biosynthesis. Bioorg. Med. Chem. Lett. 18, 5963–5966 Lu, X., Zhou, R., Sharma, I., Li, X., Kumar, G., Swaminathan, S., Tonge, P. J., and Tan, D. S. (2012) Stable analogues of OSB-AMP: potent inhibitors of MenE, the o-succinylbenzoate-CoA synthetase from bacterial menaquinone biosynthesis. ChemBioChem 13, 129 –136 Gulick, A. M. (2009) Conformational dynamics in the acyl-CoA synthetases, adenylation domains of nonribosomal peptide synthetases, and firefly luciferase. ACS Chem. Biol. 4, 811– 827 Berg, P. (1956) Acyl adenylates: an enzymatic mechanism of acetate activation. J. Biol. Chem. 222, 991–1013 Kim, Y. S., and Kang, S. W. (1994) Steady-state kinetics of malonyl-CoA synthetase from Bradyrhizobium japonicum and evidence for malonylAMP formation in the reaction. Biochem. J. 297, 327–333 Vessey, D. A., and Kelley, M. (2001) Characterization of the reaction mechanism for the XL-I form of bovine liver xenobiotic/medium-chain fatty acid:CoA ligase. Biochem. J. 357, 283–288 Horswill, A. R., and Escalante-Semerena, J. C. (2002) Characterization of the propionyl-CoA synthetase (PrpE) enzyme of Salmonella enterica: residue Lys592 is required for propionyl-AMP synthesis. Biochemistry 41, 2379 –2387 Li, H., Melton, E. M., Quackenbush, S., DiRusso, C. C., and Black, P. N. (2007) Mechanistic studies of the long chain acyl-CoA synthetase Faa1p from Saccharomyces cerevisiae. Biochim. Biophys. Acta 1771, 1246 –1253 Wu, R., Cao, J., Lu, X., Reger, A. S., Gulick, A. M., and Dunaway-Mariano, D. (2008) Mechanism of 4-chlorobenzoate:coenzyme A ligase catalysis. Biochemistry 47, 8026 – 8039 Conti, E., Franks, N. P., and Brick, P. (1996) Crystal structure of firefly luciferase throws light on a superfamily of adenylate-forming enzymes. Structure 4, 287–298 Conti, E., Stachelhaus, T., Marahiel, M. A., and Brick, P. (1997) Structural basis for the activation of phenylalanine in the nonribosomal biosynthesis of gramicidin S. EMBO J. 16, 4174 – 4183 May, J. J., Kessler, N., Marahiel, M. A., and Stubbs, M. T. (2002) Crystal structure of DhbE, an archetype for aryl acid activating domains of modular nonribosomal peptide synthetases. Proc. Natl. Acad. Sci. U.S.A. 99, 12120 –12125 Marahiel, M. A., Stachelhaus, T., and Mootz, H. D. (1997) Modular peptide synthetases involved in nonribosomal peptide synthesis. Chem. Rev. 97, 2651–2674 Chang, K.-H., Xiang, H., and Dunaway-Mariano, D. (1997) Acyl-adenylate motif of the acyl-adenylate/thioester-forming enzyme superfamily: a sitedirected mutagenesis study with the Pseudomonas sp. strain CBS3 4-chlorobenzoate:coenzyme A ligase. Biochemistry 36, 15650 –15659 Saraste, M., Sibbald, P. R., and Wittinghofer, A. (1990) The P-loop–a common motif in ATP- and GTP-binding proteins. Trends Biochem. Sci. 15, 430 – 434 Gulick, A. M., Starai, V. J., Horswill, A. R., Homick, K. M., and EscalanteSemerena, J. C. (2003) The 1.75 Å crystal structure of acetyl-CoA synthetase bound to adenosine-5⬘-propylphosphate and coenzyme A. Biochemistry 42, 2866 –2873 Reger, A. S., Carney, J. M., and Gulick, A. M. (2007) Biochemical and crystallographic analysis of substrate binding and conformational changes in acetyl-CoA synthetase. Biochemistry 46, 6536 – 6546 Branchini, B. R., Southworth, T. L., Murtiashaw, M. H., Wilkinson, S. R., Khattak, N. F., Rosenberg, J. C., and Zimmer, M. (2005) Mutagenesis evidence that the partial reactions of firefly bioluminescence are catalyzed by different conformations of the luciferase C-terminal domain. Biochemistry 44, 1385–1393 Drake, E. J., Nicolai, D. A., and Gulick, A. M. (2006) Structure of the EntB multidomain nonribosomal peptide synthetase and functional analysis of its interaction with the EntE adenylation domain. Chem. Biol. 13, 409 – 419 Yonus, H., Neumann, P., Zimmermann, S., May, J. J., Marahiel, M. A., and Stubbs, M. T. (2008) Crystal structure of DltA. Implications for the reac-

23982 JOURNAL OF BIOLOGICAL CHEMISTRY

27.

28.

29.

30.

31.

32.

33.

34.

35.

36.

37.

38.

39.

40.

41.

42.

43.

44.

tion mechanism of nonribosomal peptide synthetase adenylation domains. J. Biol. Chem. 283, 32484 –32491 Du, L., He, Y., and Luo, Y. (2008) Crystal structure and enantiomer selection by D-alanyl carrier protein ligase DltA from Bacillus cereus. Biochemistry 47, 11473–11480 Reger, A. S., Wu, R., Dunaway-Mariano, D., and Gulick, A. M. (2008) Structural characterization of a 140° domain movement in the two-step reaction catalyzed by 4-chlorobenzoate:CoA ligase. Biochemistry 47, 8016 – 8025 Gulick, A. M., Lu, X., and Dunaway-Mariano, D. (2004) Crystal structure of 4-chlorobenzoate:CoA ligase/synthetase in the unliganded and aryl substrate-bound states. Biochemistry 43, 8670 – 8679 Kochan, G., Pilka, E. S., von Delft, F., Oppermann, U., and Yue, W. W. (2009) Structural snapshots for the conformation-dependent catalysis by human medium-chain acyl-coenzyme A synthetase ACSM2A. J. Mol. Biol. 388, 997–1008 Sundlov, J. A., Fontaine, D. M., Southworth, T. L., Branchini, B. R., and Gulick, A. M. (2012) Crystal structure of firefly luciferase in a second catalytic conformation supports a domain alternation mechanism. Biochemistry 51, 6493– 6495 Osman, K. T., Du, L., He, Y., and Luo, Y. (2009) Crystal structure of Bacillus cereus D-alanyl carrier protein ligase (DltA) in complex with ATP. J. Mol. Biol. 388, 345–355 Hisanaga, Y., Ago, H., Nakagawa, N., Hamada, K., Ida, K., Yamamoto, M., Hori, T., Arii, Y., Sugahara, M., Kuramitsu, S., Yokoyama, S., and Miyano, M. (2004) Structural basis of the substrate specific two-step catalysis of long chain fatty acyl-CoA synthetase dimer. J. Biol. Chem. 279, 31717–31726 Grisostomi, G., Kast, P., Pulido, R., Huynh, J., and Hilvert, D. (1997) Efficient in vivo synthesis and rapid purification of chorismic acid using an engineered Escherichia coli strain. Bioorg. Chem. 25, 297–305 Jiang, M., and Guo, Z. (2007) Effects of macromolecular crowding on the intrinsic catalytic efficiency and structure of enterobactin-specific isochorismate synthase. J. Am. Chem. Soc. 129, 730 –731 Jiang, M., Cao, Y., Guo, Z.-F., Chen, M., Chen, X., and Guo, Z. (2007) Menaquinone biosynthesis in Escherichia coli: identification of 2-succinyl-5-enolpyruvyl-6-hydroxyl-3-cyclohexene-1-carboxylate (SEPHCHC) as a novel intermediate and re-evaluation of MenD activity. Biochemistry 46, 10979 –10989 Jiang, M., Chen, M., Cao, Y., Yang, Y., Sze, K. H., Chen, X., and Guo, Z. (2007) Determination of the stereochemistry of 2-succinyl-5-enolpyruvyl6-hydroxy-3-cyclohexene-1-carboxylic acid, a key intermediate in menaquinone biosynthesis. Org. Lett. 9, 4765– 4767 Jiang, M., Chen X., Guo Z.-F., Cao Y., Chen M., and Guo, Z. (2008) Identification and characterization of (1R,6R)-2-succinyl-6-hydroxy-2,4-cyclohexadiene-1-carboxylate synthase in the menaquinone biosynthesis of. Escherichia coli. Biochemistry 47, 3426 –3434 Jiang, M., Chen, X., Wu, X. H., Chen, M., Wu, Y. D., and Guo, Z. (2009) Catalytic mechanism of SHCHC synthase in the menaquinone biosynthesis of Escherichia coli: identification and mutational analysis of the active site residues. Biochemistry 48, 6921– 6931 Sun, Y., Yin, S., Feng, Y., Li, J., Zhou, J., Liu, C., Zhu, G., and Guo, Z. (2014) Molecular basis of the general base catalysis of an ␣/␤-hydrolase catalytic triad. J. Biol. Chem. 289, 15867–15879 Jiang, M., Chen, M., Guo, Z. F., and Guo, Z. (2010) A bicarbonate cofactor modulates 1,4-dihydroxy-2-naphthoyl coenzyme A synthase in menaquinone biosynthesis of Escherichia coli. J. Biol. Chem. 285, 30159 –30169 Sun, Y., Song, H., Li, J., Jiang, M., Li, Y., Zhou, J., and Guo, Z. (2012) Active site binding and catalytic role of bicarbonate in 1,4-dihydroxy-2-naphthoyl coenzyme A synthases from vitamin K biosynthetic pathways. Biochemistry 51, 4580 – 4589 Chen, M., Jiang, M., Sun, Y., Guo, Z. F., and Guo, Z. (2011) Stabilization of the second oxyanion intermediate by 1,4-dihydroxy-2-naphthoyl coenzyme A synthase of the menaquinone pathway: spectroscopic evidence of the involvement of a conserved aspartic acid. Biochemistry 50, 5893–5904 Winn, M. D., Ballard, C. C., Cowtan, K. D., Dodson, E. J., Emsley, P., Evans, P. R., Keegan, R. M., Krissinel, E. B., Leslie, A. G., McCoy, A., McNicholas, S. J., Murshudov, G. N., Pannu, N. S., Potterton, E. A., Powell, H. R., et al.

VOLUME 290 • NUMBER 39 • SEPTEMBER 25, 2015

ATP Configures Adenylate-forming Conformation

45.

46.

47.

48. 49.

50.

51.

52.

53.

54.

55.

(2011) Overview of the CCP4 suite and current developments. Acta Crystallogr. D Biol. Crystallogr. 67, 235–242 McCoy, A. J., Grosse-Kunstleve, R. W., Adams, P. D., Winn, M. D., Storoni, L. C., and Read, R. J. (2007) Phaser crystallographic software. J. Appl. Crystallogr. 40, 658 – 674 Hallows, W. C., Lee, S., and Denu, J. M. (2006) Sirtuins deacetylate and activate mammalian acetyl-CoA synthetases. Proc. Natl. Acad. Sci. U.S.A. 103, 10230 –10235 Adams, P. D., Afonine, P. V., Bunko´czi, G., Chen, V. B., Davis, I. W., Echols, N., Headd, J. J., Hung, L. W., Kapral, G. J., Grosse-Kunstleve, R. W., McCoy, A. J., Moriarty, N. W., Oeffner, R., Read, R. J., Richardson, D. C., et al. (2010) PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr. D Biol. Crystallogr. 66, 213–221 Emsley, P., Lohkamp, B., Scott, W. G., and Cowtan, K. (2010) Features and development of Coot. Acta Crystallogr. D Biol. Crystallogr. 66, 486 –501 Vagin, A. A., Steiner, R. A., Lebedev, A. A., Potterton, L., McNicholas, S., Long, F., and Murshudov, G. N. (2004) REFMAC5 dictionary: organization of prior chemical knowledge and guidelines for its use. Acta Crystallogr. D Biol. Crystallogr. 60, 2184 –2195 Winn, M. D., Isupov, M. N., and Murshudov, G. N. (2001) Use of TLS parameters to model anisotropic displacements in macromolecular refinement. Acta Crystallogr. D Biol. Crystallogr. 57, 122–133 Moriarty, N. W., Grosse-Kunstleve, R. W., and Adams, P. D. (2009) Electronic ligand builder and optimization workbench (eLBOW): a tool for ligand coordinate and restraint generation. Acta Crystallogr. D Biol. Crystallogr. 65, 1074 –1080 Laskowski, R. A., Macarthur, M. W., Moss, D. S., and Thornton, J. M. (1993) PROCHECK: a program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 26, 283–291 Chen, V. B., Arendall, W. B., 3rd, Headd, J. J., Keedy, D. A., Immormino, R. M., Kapral, G. J., Murray, L. W., Richardson, J. S., and Richardson, D. C. (2010) MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallogr. D Biol. Crystallogr. 66, 12–21 Larkin, M. A., Blackshields, G., Brown, N. P., Chenna, R., McGettigan, P. A., McWilliam, H., Valentin, F., Wallace, I. M., Wilm, A., Lopez, R., Thompson, J. D., Gibson, T. J., and Higgins, D. G. (2007) ClustalW and ClustalX version 2.0. Bioinformatics 23, 2947–2948 Robert, X., and Gouet, P. (2014) Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res. 42, W320 –W324

SEPTEMBER 25, 2015 • VOLUME 290 • NUMBER 39

56. DeLano, W. L. (2002) The PyMOL Molecular Graphics System, DeLano Scientific, San Carlos, CA 57. Krissinel, E., and Henrick, K. (2007) Inference of macromolecular assemblies from crystalline state. J. Mol. Biol. 372, 774 –797 58. Dundas, J., Ouyang, Z., Tseng, J., Binkowski, A., Turpaz, Y., and Liang, J. (2006) CASTp: computed atlas of surface topography of proteins with structural and topographical mapping of functionally annotated residues. Nucleic Acids Res. 34, W116 –W118 59. Trott, O., and Olson, A. J. (2010) AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 31, 455– 461 60. Mehler, E. L., and Solmajer, T. (1991) Electrostatic effects in proteins: comparison of dielectric and charge models. Protein Eng. 4, 903–910 61. Tian, Y. (2008) The Reaction Kinetics, Crystal Structures and Novel Inhibitors of Bacterial OSB-CoA Synthetase. Ph.D. thesis, University of Illinois at Chicago 62. Du, L., and Luo, Y. (2014) Thiolation-enhanced substrate recognition by D-alanyl carrier protein ligase DltA from Bacillus cereus. F1000Research 3, 106 63. Gocht, M., and Marahiel, M. A. (1994) Analysis of core sequences in the D-Phe activating domain of the multifunctional peptide synthetase TycA by site-directed mutagenesis. J. Bacteriol. 176, 2654 –2662 64. Hamoen, L. W., Eshuis, H., Jongbloed, J., Venema, G., and van Sinderen, D. (1995) A small gene, designated comS, located within the coding region of the fourth amino acid-activation domain of srfA, is required for competence development in Bacillus subtilis. Mol. Microbiol. 15, 55– 63 65. Saito, M., Hori, K., Kurotsu, T., Kanda, M., and Saito, Y. (1995) Three conserved glycine residues in valine activation of gramicidin S synthetase 2 from Bacillus brevis. J. Biochem. 117, 276 –282 66. Stuible, H., Bu¨ttner, D., Ehlting, J., Hahlbrock, K., and Kombrink, E. (2000) Mutational analysis of 4-coumarate: CoA ligase identifies functionally important amino acids and verifies its close relationship to other adenylateforming enzymes. FEBS Lett. 467, 117–122 67. Branchini, B. R., Murtiashaw, M. H., Magyar, R. A., and Anderson, S. M. (2000) The role of lysine 529, a conserved residue of the acyl-adenylateforming enzyme superfamily, in firefly luciferase. Biochemistry 39, 5433–5440 68. Jogl, G., and Tong, L. (2004) Crystal structure of yeast acetyl-coenzyme A synthetase in complex with AMP. Biochemistry 43, 1425–1431

JOURNAL OF BIOLOGICAL CHEMISTRY

23983