Quantum-Mechanical Study on the Catalytic Mechanism of Alkaline

Feb 13, 2017 - High-Quality Dataset of Protein-Bound Ligand Conformations and Its Application to Benchmarking Conformer Ensemble Generators. Journal o...
0 downloads 17 Views 2MB Size
Subscriber access provided by McMaster University Library

Article

Quantum Mechanical Study on the Catalytic Mechanism of Alkaline Phosphatases Gabriela L. Borosky J. Chem. Inf. Model., Just Accepted Manuscript • DOI: 10.1021/acs.jcim.6b00755 • Publication Date (Web): 13 Feb 2017 Downloaded from http://pubs.acs.org on February 19, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of Chemical Information and Modeling is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Quantum Mechanical Study on the Catalytic Mechanism of Alkaline Phosphatases

Gabriela L. Borosky*

INFIQC, CONICET and Departamento de Química Teórica y Computacional, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba, Ciudad Universitaria, Córdoba 5000, Argentina. Tel/Fax: +54-351-535-3853. E-mail: [email protected]

Abstract : Alkaline phosphatases (APs) catalyze the hydrolysis and transphosphorylation of phosphate monoesters. The catalytic mechanism was examined by quantum-mechanical calculations, using an active-site model based on the X-ray crystal structure of the human placental AP (PLAP). Free energies of activation and of reaction of the catalytic steps were evaluated for a series of aryl and alkyl phosphate esters, and computational results were compared with experimental values available in the literature. Mechanistical observations previously reported in experimental works were rationalized by the present theoretical study, particularly regarding the difference in the rate-determining step between aryl and alkyl phosphates. The formation rate of the covalent phosphoserine intermediate followed a linear free-energy relationship (LFER) with the pKa of the leaving group. This LFER, that could only be experimentally determined for less reactive alkyl phosphates, was verified by the present calculations to apply for the entire set of aryl and alkyl phosphate substrates.

1 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Introduction

The alkaline phosphatase (AP) superfamily involves a large group of homodimeric metalloenzymes present in almost all living species.1 The amino acid sequence of APs from different organisms is very well conserved, especially residues in the active site and surrounding regions.2,3 As the sequences of mammalian APs fit to the structure of the E. coli bacterial enzyme,2 the catalytic mechanism determined for the latter has been proposed to be similar for eukaryotic APs.4 These enzymes catalyze the hydrolysis and transphosphorylation of a broad range of phosphate monoesters by a mechanism involving formation of a covalent phosphoserine intermediate and release of inorganic phosphate and an alcohol.5 A two-step reaction mechanism (Scheme 1) has been proposed on the basis of kinetic and biochemical data.6-8 The first chemical step implies the generation of the covalent intermediate E-P (k2, Scheme 1), whilst in the second chemical step this phosphoserine is hydrolyzed (k3, Scheme 1). Transphosphorylation to a phosphate acceptor (k5, Scheme 1) takes place in nucleophilic buffers (presence of R´OH).

Scheme 1. General catalytic mechanism of APs.

2 ACS Paragon Plus Environment

Page 2 of 35

Page 3 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

The rate-limiting step is pH dependent for aryl phosphate substrates, the hydrolysis of the covalent intermediate E-P (k3, Scheme 1) being rate-determining at acidic pH < 7.5, whereas the release of phosphate from the noncovalent enzyme-phosphate complex (E.Pi) (k4, Scheme 1) becomes rate-limiting under basic conditions (pH > 7.5).9-11 In contrast, the phosphorylation of the enzyme (k2) has been proposed as the rate-determining step for reaction of alkyl phosphates.12 Human placental AP (PLAP) is one of the four human AP isoenzymes.13 It has been suggested as modulator of fetal growth because increases growth and survival of fetal cells.14,15 Since PLAP is one of the proteins ectopically expressed by tumor cells,16 it has been used as a tumor marker in several clinical reports.17-20 In this way, this enzyme is presumed to play an interesting role in cancer diagnosis and therapy. Molecular docking, molecular dynamics, and quantum mechanical calculations were previously applied to the study of the catalytic mechanism of PLAP.21,22 Activation barriers and free energy values for each reaction step were computed for methylphosphate21 and phenylphosphate22 as substrates. The aim of this work was to achieve a better understanding of the activity of APs by means of additional quantum chemical calculations on the catalytic mechanism of PLAP as a model AP. Continuing our former studies,21,22 in this instance the reaction steps of the mechanism were evaluated for diverse aryl and alkyl phosphate substrates. The existence of a dependence of the rate of catalysis on the pKa of the leaving group was inspected. The present computational results were validated by comparison with experimental assessments informed in the literature.

3 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Computational Methods

General Methodology. Procedures were based on the cluster approach, which has been successfully applied to modeling enzymatic reactions.23-25 In this technique, a limited part of the enzyme is carefully selected to appropriately represent the active site at a quantum mechanical level. The protein residues surrounding the active site that are not explicitly included in the calculations may influence the model in two main modes. First, by causing steric restraints which, if not considered, could generate significant artificial repositioning of the model residues because of the absence of the surrounding amino acids. Therefore, certain key coordinates at the boundary of the system are conserved as in the Xray structures to model the steric interactions. In this way, the positions of the backbone atoms are kept fixed when obtained from high-resolution X-ray structures, while the coordinates of all the other atoms in the model system are relaxed. Secondly, polarization caused by the enzyme surroundings can affect the computed energies; for this reason, polarizable continuum methods are usually employed to take into account electrostatic effects, using a dielectric constant ε = 4 to represent the protein environment.23-25 It has been proved that relative solvation effects decrease very fast as the model size increase, since more groups providing polarization are explicitly added. This combination of a coordinate-locking procedure and continuum solvation effectively accounts for the residues not included in the model, affording precise calculated energies.23-25 The cluster model is the quantum mechanics-only (QM-only) method of the two popular quantum chemical approaches to address enzymatic reaction mechanisms. The other approach is the quantum mechanics/molecular mechanics (QM/MM) model, in which the entire enzyme is

4 ACS Paragon Plus Environment

Page 4 of 35

Page 5 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

considered with a small core described quantum mechanically, while the rest is treated by molecular mechanics.26,27 In this work theoretical calculations were performed with a model of 219 quantum atoms, a size which is expected to yield accurate results by the cluster method.23-25 The starting geometry was built from an X-ray structure with an atomic resolution of 1.8 Å (PDB entry 1EW2).28 Former molecular dynamics simulations had verified the stability of the crystal structure, particularly within the active site region.21 The coordinates of the backbone atoms were kept fixed to avoid distortion of the active site during geometry optimization of the side chains of all residues, as this procedure had been previously confirmed to be appropriate.21,22 The influence of the environment was considered by using a polarizable implicit continuum solvation model with a dielectric constant ε = 4.0 to represent the rest of the amino acids surrounding the active site; this value is generally considered to be a good representation of the protein environment.24,25

Enzyme model. The three-dimensional structure of PLAP was obtained from the Protein Data Bank (PDB code 1EW2).28 The computational model was composed by the catalytic metal triad (two Zn2+ ions (M1 and M2) and one Mg2+ ion (M3)) with their ligands Asp316, His320, His432, His358, Asp357, Asp42, Glu311, Ser155, and the nucleophilic Ser92, as well as the significant residues Arg166 and Glu429 (hydrophilic pocket). Valence at cut peptide bonds were completed with hydrogen atoms. Six water molecules were conserved: three to complete Mg2+ coordination (one as a hydroxide ion), and three coordinated to Glu429. The phosphate ion in 1EW2 was utilized as a template for building the dianionic monoesters (charge -2e), since experimental observations had determined that these substrates bind and react as dianions.12 The side chains of arginine, aspartate, and 5 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

glutamate residues were considered as ionized, and histidines were defined as neutral, singly protonated on Hδ. The complete model system had a net charge of -1, resulting from the sum of seven positive charges (two Zn2+ cations, one Mg2+, and one arginine) and eight negative charges (three aspartates, two glutamates, one hydroxide anion, and the phosphate monoester dianion substrate).

Quantum-mechanical methods. ONIOM29 calculations employing two quantummechanical layers (QM:QM) were carried out with the Gaussian 09 suite of programs.30 Density Functional Theory (DFT) optimizations with the B3LYP31-33 functional were performed for the high layer, consisting of three metal cations, the phosphate monoester dianions, five water molecules, one hydroxide ion, carboxylate groups of glutamates and aspartates, -CH2OH groups of serines, and the -C(NH2)2 group of arginine (around 40 heavy atoms and 25 hydrogens). The 6-31+G* basis was employed for C, O, N, P, Mg and H atoms, and the pseudopotential Lanl2DZ was applied for Zn atoms. For the low layer (73 heavy atoms and 78 hydrogens) energy minimizations with the semi-empirical method PM3MM34 were performed. The ONIOM(B3LYP:PM3MM) methodology has been shown to be appropriate in our prior studies.21,22 The electrostatic effect of the environment was considered via polarized continuum model (IEFPCM)35-38 optimizations, using a dielectric constant ε = 4.0 to simulate the influence of the residues surrounding the active site. The positions of the backbone atoms (involved in peptide bonds) were held fixed to preserve the active site structure, whereas the coordinates of the rest of the atoms were fully optimized. Minima and transition states (TSs) on the potential energy surfaces were characterized by harmonic vibrational frequency calculations, which also afforded the corresponding zeropoint vibrational energies, thermal corrections, enthalpies, entropies, and free energies. 6 ACS Paragon Plus Environment

Page 6 of 35

Page 7 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

ONIOM(B3LYP/6-311+G(2d,p):PM3MM)-IEFPCM

single-point

and

frequency

computations were performed for all the optimized stationary points in order to obtain more accurate energy and free energy values.

Results and Discussion

The reaction steps of the catalytic mechanism of APs were computationally evaluated within the active site of PLAP, which includes the catalytic Ser92, the metal triplet (two Zn2+ and one Mg2+), Arg166, Glu429, and some other relevant proximate amino acids (Figure 1). The hydrophilic pocket formed by Arg166 and Glu429 is presumed to stabilize the hydrophilic moiety of the phosphate monoester ligand. Residue Glu429 also stabilizes the water molecules that bridge the gap to the phosphate group of the phosphoseryl intermediate.39 One of these waters is significantly conserved because it is the nucleophile involved in the hydrolysis of the covalent phosphoserine.

7 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. Active site of PLAP with relevant adjacent residues. Fixed atoms are indicated with asterisks. 8 ACS Paragon Plus Environment

Page 8 of 35

Page 9 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

In this work the mechanism of catalysis in the active site of PLAP was modelled for the dianions of several primary phosphate esters: two aryl phosphates, and nine alkyl phosphates (Figure 2). These substrates were selected from a previous experimental study where the corresponding kcat/KM values for the E. coli AP-catalyzed hydrolysis were reported.12

Figure 2.

Dianions of phosphate esters selected as substrates.

9 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

Page 10 of 35

1 2 3 For the aryl phosphate substrates 1 and 2, the catalytic mechanism illustrated in 4 5 Scheme 2 was evaluated inside the active site, starting from the Michaelis complex. The 6 7 8 initial coordinates for this structure were determined by previous molecular docking 9 10 calculations.22 In this complex, the oxygen of the ester leaving group was coordinated to 11 12 Zn1, one nonbridging oxygen was coordinated to Zn2, and the other two nonbridging 13 14 15 oxygen atoms presented hydrogen bond contacts with the guanidinium group of Arg166. 16 17 After energy minimization, the water molecule originally included as a hydroxide ion was 18 19 protonated by the hydroxyl group of catalytic Ser92, which became ionized. This 20 21 22 spontaneous proton transfer, also observed in our previous works,21,22 had been proposed in 23 24 the literature.4,12,40 25 26 27 28 29 30 2+ Zn1 2+ 2+ 31 Zn1 O Zn1 O O 32 O-Ser92 Step 1 33 Step 2 P P O-Ser92 ArO OP O-Ser92 34 O O ArO H O O 2 35 OOO 36 O O ArO C Michaelis complex 37 Covalent intermediate C C 38 Glu429 O Glu429 Glu429 O 39 O Step 3 40 41 42 43 2+ 2+ 2+ 44 Zn1 Zn1 Zn1 O O 45 O Step 4 Step 5 O-Ser92 46 O-Ser92 P O-Ser92 P P O 47 HO O HO OHO O48 OOH OH OH 49 O C 50 C C 51 Glu429 Glu429 O 52 Glu429 O O 53 54 55 Scheme 2. Mechanism for the hydrolysis of aryl phosphate dianions. 56 57 58 59 10 60 ACS Paragon Plus Environment

Page 11 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

After formation of the covalent phosphoserine intermediate (Scheme 2, Step 1), the aromatic oxide leaving group was stabilized by coordination to Zn1. This anion was displaced from metal coordination by one water molecule in Step 2, and a Zn1-coordinated hydroxide ion was formed in Step 3 by proton transfer to Glu429. In Step 4, the nucleophilic hydroxide attacked the phosphoserine, releasing a hydrogen phosphate dianion and regenerating the nucleophilic serine oxyanion. A proton transfer returned Glu429 to its initial ionized form, generating dihydrogen phosphate ion as the final product in Step 5. The structures of selected stationary points are shown in Figure 3, and calculated free energies are presented in Table 1.

11 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. Model stationary points for p-NO2-phenyl phosphate as substrate (bond distances in Å). (a) Michaelis complex. (b) TS for Ser92 nucleophilic attack. (c) Covalent phosphoserine intermediate, water complexed to Zn1. (d) TS for phosphoserine hydrolysis. 12 ACS Paragon Plus Environment

Page 12 of 35

Page 13 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Table 1. Reaction and activation free energies computed for the hydrolysis of phosphate ester dianions (kcal/mol).a Reaction steps are shown in Schemes 2 and 3. Substrate (pKa,b kcat/KM,c ∆G≠obsd) 1 (7.1, 3.3 x 107, 7.2) 2 (10.0, 2.4 x 107, 7.4) 3 7

(12.4, 2.0 x 10 , 7.5) 4 (13.6, 2.8 x 107, 7.3) 5 (14.0, 7.5 x 106, 8.1) 6 (14.2, 2.3 x 106, 8.8)

Step 1

Step 2

Step 3e

∆G≠

∆Gr

∆Gr

∆Gr

∆G≠

∆Gr

-9.7

4.8

-6.1

10.8

8.3

15.9

-12.8

-1.9

6.4

-8.8

10.8

8.3

15.9

-12.8

5.9

14.3

-9.3

1.2

-10.5

11.3

(3.5)

(-4.5)

(15.9)

7.6

13.5

-7.9

-0.4

-10.5

11.3

(0.4)

(-4.5)

(15.9)

9.3

15.2

-10.0

0.5

-10.5

11.3

(0.8)

(-4.5)

(15.9)

6.7

13.8

-4.5

0.2

-10.5

11.3

(-2.1)

(-4.5)

(15.9)

13.1

16.3

-6.6

-2.6

-10.5

11.3

(-6.4)

(-4.5)

(15.9)

14.4

17.4

-10.1

-0.6

-10.5

11.3

(-4.2)

(-4.5)

(15.9)

13.9

18.8

-10.7

-1.7

-10.5

11.3

(-3.1)

(-4.5)

(15.9)

14.6

18.0

-11.0

-0.1

-10.5

11.3

(-3.5)

(-4.5)

(15.9)

15.4

18.9

-12.2

-0.1

-10.5

11.3

(-3.1)

(-4.5)

(15.9)

5f

3.5 x 10 , 9.9) 8 (15.9, 1.5 x 105, 10.4) 9 5

(16.1, 1.3 x 10 , 10.5) 10 (16.1, 8.6 x 104, 10.7) 11 4

(16.1, 4.2 x 10 , 11.1) a b

Step 5

∆Gr

7 (15.5, 1.2 x 106, 9.2

Step 4e

At the ONIOM(B3LYP/6-311+G(2d,p):PM3MM//B3LYP/6-31+G(d):PM3MM)-IEFPCM level. pKa of the alcohol leaving group from reference 41. c From reference 12 (M-1s-1). d Calculated 13 ACS Paragon Plus Environment

-

-

-

-

-

-

-

-

-

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

from kcat/KM values. e Values in parenthesis were obtained by using a neutral active site model, similar as for aryl phosphates (see text). f From reference 42.

The mechanism in Scheme 3 was considered for the alkyl phosphate substrates. With these alkyl oxide leaving groups, the difference in pKa values allowed the proton transfer from the Zn1-coordinated water to generate an alkyl alcohol and hydroxide ion in Step 3 (pKa H2O = 15.7, pKa MeOH = 15.5, pKa BuOH = 16.1; on the other hand, pKa PhOH = 10.0 (Scheme 2)).41 Therefore, hydrogen phosphate dianion was the final product for these substrates. Computed free energies are displayed in Table 1, and the structures of model stationary points are illustrated in Figure 4.

Scheme 3. Mechanism for the hydrolysis of alkyl phosphate dianions.

14 ACS Paragon Plus Environment

Page 14 of 35

Page 15 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Figure 4. Model stationary points for ethyl phosphate as substrate (bond distances in Å). (a) Michaelis complex. (b) TS for Ser92 nucleophilic attack. (c) Covalent phosphoserine intermediate, water complexed to Zn1. (d) TS for phosphoserine hydrolysis. 15 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The pathway was computed from the respective Michaelis complexes, whose initial coordinates were derived from previous molecular docking calculations.21 The potential energy surfaces calculated for all the phosphate esters in this study shared some common features. The Michaelis noncovalent complex converted into the covalent phosphoserine through a small rotation of the side chain of Ser92. The TS for the first chemical step, involving generation of the covalent phosphoserine (Step 1 in Schemes 2 and 3), presented a trigonal bipyramidal configuration corresponding to an in-line displacement reaction, as proposed by Holtz et al.,43 with bond forming/breaking oxygens in an opposite axial disposition. The three nonbridging oxygens of the transferred phosphoryl group bisected the axial plane and formed stabilizing interactions with Arg166, Zn2, and one Mg3-bound water. In all cases, departure of the leaving group was assisted by coordination with Zn1. The nucleophilic attack by a hydroxide anion coordinated to Zn1 on the phosphoseryl intermediate (second chemical step, Step 4 in Schemes 2 and 3) also exhibited a trigonal bipyramidal TS corresponding to another in-line displacement. The stabilizing interactions of the phosphoryl group in this transition structure were similar to those observed in the first TS. In general, the positions of the model residues experienced only minor variations, with negligible modifications in the atomic coordinates during the course of the catalytic mechanism. This observation suggests that the active site structure presents an optimal configuration to assist the hydrolysis of phosphate monoesters. Keynote theoretical investigations on the mechanism of APs have been recently reviewed.44,45 Comparison with experimental assessments allows validation of computational results. Reproduction of experimental observations is a proper test of the accuracy of theoretical procedures, and supports their interpretations and predictions. However, 16 ACS Paragon Plus Environment

Page 16 of 35

Page 17 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

comparison of theoretical and experimental kinetic data is not straightforward. Whereas computed barriers correspond to the free energy difference between the TS and the Michaelis complex (∆Gcat≠) and are related to kcat, experimental reported values are usually kcat/KM, which provide the free energy barrier as measured from the ground state of the free enzyme and substrate in solution to the TS for the first irreversible reaction step (∆Gobs≠). Both free energy barriers are related through the binding free energy of the substrate (∆Gobs≠ = ∆Gcat≠ - ∆GBinding). Thus, ∆Gobs≠ derived from the experimental kcat/KM value is a lower bound for the calculated free energy of activation of a chemical step (∆Gcat≠). Table 1 presents the free energies of activation afforded by the present computations, along with the corresponding values from experimental studies reported in the literature12 (∆G≠obs values were determined from kcat/KM by transition-state theory46). Calculated activation barriers were higher than those derived from kcat/KM values, supporting the reliability of the computational methodology and level of theory employed in this work. Although hydrolysis of the phosphoserine is the same reaction for both type of substrates (Step 4 in Schemes 2 and 3), the activation barrier for this step was 4.6 kcal/mol lower with alkyl phosphates than for aryl phosphates. This hydrolysis reaction was exothermic for alkyl phosphates, but endothermic with the aromatic substrates (Table 1). Even though the model system was the same for both type of substrates, the total charge was -1e after departure of an alkyl alcohol in Scheme 3, while the system became neutral after removal of an aryl oxide in Scheme 2. It has been suggested that activation energies for phosphoryl transfer reactions in metalloenzymes are very sensitive to the charge balance around the active center, with neutral models being more suitable for reproducing experimental data.47 An activation free energy of 16.2 kcal/mol was determined for the hydrolysis of the covalent intermediate by measuring the rate constant for breakdown of a 17 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

radioactive E-P complex (k3 in Scheme 1).48 It is worth of noting the very good agreement between this value and the free energy barrier of 15.9 kcal/mol calculated for aryl phosphate substrates (Table 1). Phosphorylation of the enzyme by alkyl phosphates should have a barrier higher than 16.2 kcal/mol, since experimental evidences have indicated k2 as rate-limiting.12 The present computational results yielded barriers of 16.3-18.9 kcal/mol for methyl, ethyl, propyl, butyl and iso-butyl phosphates, in accordance with these observations. Free energy profiles are displayed in Figure 5. The two chemical steps of the catalytic mechanism are Step 1 (Schemes 2 and 3, corresponding to k2 in Scheme 1), and Step 4 (Schemes 2 and 3, equivalent to k3 in Scheme 1). Activation free energies in Table 1 indicate that the hydrolysis of the covalent intermediate (k3) is the rate-limiting step for aryl phosphate monoesters. On the other hand, for alkyl phosphate substrates the computed values point to the nucleophilic attack of the serine alkoxide to phosphorus (k2) as the ratedetermining step. While for aryl phosphates the rate-determining step at pH > 7.5 corresponds to product release (k4 in Scheme 1), and to the hydrolysis of the phosphoserine (k3) at pH < 7.5,9-11 phosphorylation of the enzyme (k2) has been hinted as rate-limiting for alkyl phosphates.12 Hence, the present computational results are in accordance with the mechanistic interpretations arisen from experimental observations.

18 ACS Paragon Plus Environment

Page 18 of 35

Page 19 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Figure 5. Free energy profiles for the catalytic mechanism. (a) Free energy barriers for Step 4 calculated according to Scheme 3 for alkyl phosphates. (b) Free energy barriers for Step 4 computed according to Scheme 2 for all substrates (neutral model). 19 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The more reactive alkyl phosphate substrates 3-6 represent an intermediate case. According to calculations for Scheme 3, the higher barrier corresponds to the serine attack to phosphorus in Step 1 (k2 in Scheme 1, Figure 5a). However, if a neutral active site model is considered in the second chemical step (as for aryl phosphates in Scheme 2), the TS for hydrolysis of the phosphoserine is the highest in energy (Step 4, k3 in Scheme 1, Figure 5b). Taking into account that trifluoroethyl (3) and propargyl (4) phosphates have values of kcat/KM that are similar to those for p-nitrophenyl and phenyl phosphates, these most reactive alkyl phosphates were suspected to be limited by the same process as aryl phosphates; nevertheless, test experiments for 3 did not evidence a change in ratedetermining step as compared to ethyl phosphate (8).12 By inspection of the potential energy surfaces calculated for both type of substrates (Table 1, Figure 5), the following factors appears to determine the change in rate-limiting step between alkyl and aryl phosphates. The basicity of the corresponding leaving groups seems to be the major aspect influencing the barrier height of the first chemical step (k2), as departure of a more stable aromatic anion is thermodynamically favored. Regarding the second chemical step (k3), the presence of a more basic alkyl oxide anion would help the proton transfer from the Zn1-coordinated water molecule to generate the nucleophilic hydroxide (Step 3 in Scheme 3), assisting the hydrolysis of the phosphoserine intermediate; on the other hand, formation of the hydroxide ion would be less facile in the presence of an aryl substrate (Step 3 in Scheme 2). Therefore, the second energy barrier becomes the highest for aryl phosphates, with the enzymatic rate-determining step being highly influenced by the pH of the media (at pH < 7.5 k3 is rate-limiting; at pH > 7.5 k4 is ratelimiting).9-11 A pKa of 8.0 has been suggested for the zinc-coordinated water.12

20 ACS Paragon Plus Environment

Page 20 of 35

Page 21 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Linear free-energy relationships (LFERs). When applied to enzymatic reactions, LFERs or Brønsted correlations can provide information about the nature of the TS if the chemical step is rate-limiting and substituent-specific effects are negligible.49 For a series of rate constants with different leaving groups, the Brønsted parameter βlg, which is the slope of the plot of log(k) versus pKa of the leaving group, describes the change in effective charge in going from reactants to the TS. The AP-catalyzed hydrolysis of aryl Ophosphorothioates follows a steep leaving group dependence (βlg = -0.77),50 whereas for alkyl phosphates βlg = -0.85,12 both values being consistent with a dissociative TS. Based on the above experimental observations, log(k2) values (derived from the calculated activation free energies by TS theory,46 Table 2) were plotted against the corresponding pKas in order to gauge the existence of a LFER. A linear correlation with a coefficient R2 = 0.933 and a slope of -1.15 was determined (Figure 6), which confirms the strong dependence in the formation rate of the covalent phosphoserine intermediate with leaving group stabilization. It is worth of mentioning that, as the AP-catalyzed reactions of aryl phosphates are not limited by this chemical step, the relationship of kcat/KM with the leaving group pKa could not be experimentally observed with aromatic substrates.12 As the pKa depends on the electrostatic properties of the local environment,51 the pKa values of the leaving groups within a protein should differ from the values in solution reported in reference 41. However, computational approaches to calculate pKa values in enzymes are very demanding and converge slowly (see, as examples, references 52 and 53). Nevertheless, a very good correlation was obtained with the solution values employed in the present work.

21 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 35

Table 2. Computed values corresponding to Step 1 (Schemes 2 and 3, kcal/mol).a

Bond lengths (Å) Substrate (pKa)b

∆Gr

∆G≠

Log(k)

(Mulliken charge density at O) P-OSer92

P-Olg

1 (7.1)

-9.7

4.8

9.3

2.381 (-0.572)

2.049 (-0.472)

2 (10.0)

-1.9

6.4

8.1

2.213 (-0.581)

2.126 (-0.468)

3 (12.4)

5.9

14.3

2.3

2.159 (-0.588)

2.154 (-0.509)

4 (13.6)

7.6

13.5

2.9

2.083 (-0.593)

2.198 (-0.518)

5 (14.0)

9.3

15.2

1.6

2.092 (-0.594)

2.181 (-0.524)

6 (14.2)

6.7

13.8

2.7

2.068 (-0.594)

2.206 (-0.517)

7 (15.5)

13.1

16.3

0.8

2.011 (-0.598)

2.275 (-0.507)

8 (15.9)

14.4

17.4

0.04

2.028 (-0.597)

2.257 (-0.517)

9 (16.1)

13.9

18.8

-1.0

2.024 (-0.591)

2.225 (-0.521)

10 (16.1)

14.6

18.0

-0.4

2.026 (-0.592)

2.222 (-0.518)

11 (16.1)

15.4

18.9

-1.0

2.026 (-0.597)

2.257 (-0.516)

a

At the ONIOM(B3LYP/6-311+G(2d,p):PM3MM//B3LYP/6-31+G(d):PM3MM)-

IEFPCM level. b pKa of the alcohol leaving group from reference 41.

22 ACS Paragon Plus Environment

Page 23 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

log(k) 10 1

2

8 6 4 4

6

3

2 5

7 8

0

10 9 11

-2 6

8

10

12

14

16

pKa

Figure 6. Leaving group dependence for Step 1 of the AP-catalyzed hydrolysis of organic phosphates.

The TSs for serine attack to trifluoroethyl (3), fluoroethyl (6), and cyanoethyl phosphate (5) presented a stabilizing hydrogen bond interaction between one of the water molecules and one fluorine or nitrogen atom of the substrate. As these compounds were among the most reactive alkyl phosphates examined, this specific interaction between the active site and these particular leaving groups was taken into account to check any perturbation effect on the observed LFER. It should be noted that consistent deviations for 3 had been detected in experimental LFERs, and exclusion of this phosphate from the analyzed data improved the linear correlation.12 Thus, exclusion of 3 from correlation

23 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

yielded a coefficient R2 = 0.956 for the linear fit in Figure 6, whereas further excluding 3, 5 and 6 afforded a value of R2 = 0.963. The steep dependence for the AP-catalyzed hydrolysis of organic phosphates with the electron withdrawing properties of the leaving group suggests substantial buildup of negative charge on this group in the TS. This observation is consistent with a dissociative TS, with significant bond breaking and charge accumulation on the leaving group. Table 2 summarizes the bond distances and Mulliken charge densities in the TSs for Step 1 characterized for the eleven substrates considered. Computations indicate that a lower pKa of the leaving group gives rise to a more exothermic reaction with a lower activation barrier, and an earlier TS. That is, the general trend show that bonding between the P atom and the O of the nucleophilic Ser92 (P-OSer92), as well as bond cleavage to the leaving group (P-Olg) are both less advanced for a lower pKa value.

Concluding Remarks

The catalytic mechanism of hydrolysis in the active site of PLAP as a model AP was inspected by quantum mechanical computations for eleven phosphate esters as substrates. Stationary points on the corresponding potential energy surfaces (intermediates and TSs along the reaction pathway) were optimized and characterized. The computational free energy barriers were in good accordance with experimental studies for the APcatalyzed hydrolysis of these substrates.12 Calculated activation barriers resulted higher than those derived from kcat/KM assessments, and consistent with those derived from kcat values. The correspondence between the activation free energy of 15.9 kcal/mol calculated 24 ACS Paragon Plus Environment

Page 24 of 35

Page 25 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

for the hydrolysis of the covalent intermediate (k3 in Scheme 1) and the experimental value of 16.2 kcal/mol determined by measuring the rate constant for breakdown of a radioactive E-P complex48 was especially remarkable. Calculated results pointed to the hydrolysis of the covalent intermediate (k3) as the rate-limiting step for aryl phosphate monoesters, while the nucleophilic attack of the serine alkoxide to phosphorus (k2) is the rate-determining step proposed for alkyl phosphate substrates. These observations agree with mechanistic interpretations derived from previous experimental studies, which suggest these different rate-limiting steps for aryl and alkyl phosphates.9-12 The change in the rate-determining step between alkyl and aryl phosphates was ascribed to the basicity of the corresponding leaving groups. Regarding the first chemical step (k2), departure of a more stable aromatic alkoxide is preferred, affording a lower barrier. The presence of a more basic alkyl oxide anion favors the proton transfer from a Zn1-coordinated water to form the nucleophilic hydroxide for the second chemical step (k3, Step 3 in Scheme 3) and assists the hydrolysis of the phosphoserine intermediate; on the other hand, generation of the hydroxide ion was more energetically costly for the aryl phosphate substrates. A lower pKa of the leaving group increased the exothermicity of the covalent phosphoserine formation reaction, causing a lower activation barrier with an earlier TS. A linear correlation with a coefficient R2 = 0.933 and a slope of -1.15 was observed, indicating a strong dependence in the formation rate of the phosphoserine intermediate with the stability of the leaving group. Because the AP-catalyzed reactions of aryl phosphates are not limited by this chemical step, the relationship of kcat/KM with the leaving group pKa could only be experimentally determined for alkyl phosphate substrates.12 Interestingly, the

25 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

present computational results confirm the existence of this LFER for the complete set of aryl and alkyl phosphate esters.

Acknowledgements

The author gratefully acknowledges financial support from Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET) and the Secretaría de Ciencia y Tecnología de la Universidad Nacional de Córdoba (Secyt-UNC). Access to computational resources at Mendieta and Cristina clusters from CCAD-UNC, which is part of SNCADMinCyT, Argentina, is also acknowledged.

Supporting Information: Cartesian coordinates for relevant optimized stationary points. This material is available free of charge via the Internet at http://pubs.acs.org.

26 ACS Paragon Plus Environment

Page 26 of 35

Page 27 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

References

(1) McComb, R. B.; Bowers, G. N., Jr.; Posen, S. Measurement of Alkaline Phosphatase Activity. In Alkaline Phosphatases; Plenum Press: New York, 1979; pp 986989. (2) Kim, E. E.; Wyckoff, H. W. Structure of Alkaline Phosphatases. Clin. Chim. Acta 1990, 186, 175-178. (3) Murphy, J. E.; Tibbitts, T. T.; Kantrowitz, E. R. Mutations at Positions 153 and 328 in Escherichia Coli Alkaline Phosphatase Provide Insight Towards the Structure and Function of Mammalian and Yeast Alkaline Phosphatases. J. Mol. Biol. 1995, 253, 604617. (4) Kim, E. E.; Wyckoff, H. W. Reaction Mechanism of Alkaline Phosphatase Based on Crystal Structures : Two-metal Ion Catalysis. J. Mol. Biol. 1991, 218, 449-464. (5) Schwartz, J. H.; Lipmann, F. Phosphate Incorporation into Alkaline Phosphatase of E. coli. Proc. Natl. Acad. Sci. USA 1961, 47, 1996-2005. (6) Holtz, K. M.; Kantrowitz, E. R. The Mechanism of the Alkaline Phosphatase Reaction: Insights from NMR, Crystallography and Site-specific Mutagenesis. FEBS Lett. 1999, 462, 7-11. (7) Reid, T. W.; Wilson, I. B. Escherichia Coli Alkaline Phosphatase. In The Enzymes; Boyer, P. D., Ed.; Academic Press: New York, 1971; Vol. 4, pp 373-415. (8) Coleman, J. E. Structure and Mechanism of Alkaline Phosphatase. Annu. Rev. Biophys. Biomol. Struct. 1992, 21, 441-483.

27 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(9) Bloch, W.; Gorby, M. S. Catalytic Mechanism of Escherichia Coli Alkaline Phosphatase: Resolution of Three Variants of the Acyl-enzyme Mechanism. Biochemistry 1980, 19, 5008-5018. (10) Gettins, P.; Coleman, J. E. 31P Nuclear Magnetic Resonance of Phosphoenzyme Intermediates of Alkaline Phosphatase. J. Biol. Chem. 1983, 258, 408-416. (11) Hull, W. E.; Halford, S. E.; Gutfreund, H.; Sykes, B. D. Phosphorus-31 Nuclear Magnetic Resonance Study of Alkaline Phosphatase: The Role of Inorganic Phosphate in Limiting the Enzyme Turnover Rate at Alkaline pH. Biochemistry 1976, 15, 1547-1561. (12) O´Brien, P. J.; Herschlag, D. Alkaline Phosphatase Revisited: Hydrolysis of Alkyl Phosphates. Biochemistry 2002, 41, 3207-3225. (13) Millán, J. L. Mammalian Alkaline Phosphatases: From Biology to Applications in Medicine and Biotechnology; Wiley-VCH Verlag GmbH & Co.: Weinheim, Germany, 2006; pp 1-322. (14) She, Q.-B.; Mukherjee, J. J.; Huang, J.; Crilly, K. S.; Kiss, Z. Growth Factorlike Effects of Placental Alkaline Phosphatase in Human Fetus and Mouse Embryo Fibroblasts. FEBS Lett. 2000, 469, 163-167. (15) She, Q.-B.; Mukherjee, J. J.; Chung, T.; Kiss, Z. Placental Alkaline Phosphatase, Insulin, and Adenine Nucleotides or Adenosine Synergistically Promote Long-term Survival of Serum-starved Mouse Embryo and Human Fetus Fibroblasts. Cell. Signal. 2000, 12, 659-665 (16) Fishman, W. H.; Inglis, N. R.; Stolbach, L. L.; Krant, M. J. A Serum Alkaline Phosphatase Isoenzyme of Human Neoplastic Cell Origin. Cancer Res. 1968, 36, 935-938.

28 ACS Paragon Plus Environment

Page 28 of 35

Page 29 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

(17) Dempo, K.; Elliot, K. A. C.; Desmond, W.; Fishman, W. H. Demonstration of Gamma-glutamyl Transferase, Alkaline Phosphatase, CEA and HCG in Human Lung Cancer. Oncodev. Biol. Med. 1981, 2, 21-37. (18) Nouwen, E. J.; Pollet, D. E.; Schelstraete, J. B.; Eerdekens, M. W.; Hänsch, C.; Van de Voorde, A.; De Broe, M. E. Human Placental Alkaline Phosphatase in Benign and Malignant Ovarian Neoplasia. Cancer Res. 1985, 45, 892-902. (19) Nouwen, E. J.; Pollet, D. E. ; Eerdekens, M. W. ; Hendrix, P. G. ; Briers, T. W. ; De Broe, M. E. Immunohistochemical Localization of Placental Alkaline Phosphatase, Carcinoembryonic Antigen, and Cancer Antigen 125 in Normal and Neoplastic Human Lung. Cancer Res. 1986, 46, 866-876. (20) Wick, M. R.; Swanson, P. E.; Manivel, J. C. Placental-like Alkaline Phosphatase Reactivity in Human Tumors: An Immunohistochemical Study of 520 Cases. Hum. Pathol. 1987, 18, 946-954. (21) Borosky, G. L.; Lin, S. Computational Modeling of the Catalytic Mechanism of Human Placental Alkaline Phosphatase (PLAP). J. Chem. Inf. Model. 2011, 51, 2538-2548. (22) Borosky, G. L. Catalytic Activity of Human Placental Alkaline Phosphatase (PLAP): Insights from a Computational Study. J. Phys. Chem. B 2014, 118, 14302-14313. (23) Blomberg, M. R. A.; Borowski, T.; Himo, F.; Liao, R.-Z. Siegbahn, P. E. M. Quantum Chemical Studies of Mechanisms for Metalloenzymes. Chem. Rev. 2014, 114, 3601-3658. (24) Siegbahn, P. E. M.; Himo, F. The Quantum Chemical Cluster Approach for Modeling Enzyme Reactions. WIREs Comp. Mol. Sci. 2011, 1, 323-336 and references therein.

29 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(25)

Siegbahn, P. E. M.; Himo, F. Recent Developments of the Quantum

Chemical Cluster Approach for Modeling Enzyme Reactions. J. Biol. Inorg. Chem. 2009, 14, 643-451 and references therein. (26)

Warshel, A.; Levitt, M. Theoretical Studies of Enzymic Reactions:

Dielectric, Electrostatic, and Steric Stabilization of the Carbonium Ion in the Reaction of Lysozyme. J. Mol. Biol. 1976, 103, 227-249. (27)

Warshel, A.; Sharma, P. K.; Kato, M.; Xiang, Y.; Liu, H.; Olsson, M. H. M.

Electrostatic Basis for Enzyme Catalysis. Chem. Rev. 2006, 106, 3210-3225 and references therein. (28) Le Du, M. H.; Stigbrand, T.; Taussig, M. J.; Ménez, A.; Stura, E. A. Crystal Structure of Alkaline Phosphatase from Human Placenta at 1.8 Å Resolution. J. Biol. Chem. 2001, 276, 9158-9165. (29) Dapprich, S.; Komáromi, I.; Byun, K. S.; Morokuma, K.; Frisch, M. J. A New ONIOM Implementation in Gaussian98. Part I. The Calculation of Energies, Gradients, Vibrational Frequencies and Electric Field Derivatives. J. Mol. Struct. (Theochem) 1999, 461-462, 1–21. (30) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery, Jr., J. A.; Vreven, T.; Kudin, K. N.; Burant, J. C.; Millam, J. M.; Iyengar, S. S.; Tomasi, J.; Barone, V.; Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.; Petersson, G. A.; Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Klene, M.; Li, X.; Knox, J. E.; Hratchian, H. P.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Ayala, P. Y.; Morokuma, K.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Zakrzewski, V. 30 ACS Paragon Plus Environment

Page 30 of 35

Page 31 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

G.; Dapprich, S.; Daniels, A. D.; Strain, M. C.; Farkas, O.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Ortiz, J. V.; Cui, Q.; Baboul, A. G.; Clifford, S.; Cioslowski, J.; Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Gonzalez, C.; Pople, J. A. Gaussian 09, Revision A.01; Gaussian, Inc.: Wallingford, CT, 2009. (31) Becke, A. D. Density-functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648-5652. (32) Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti Correlationenergy Formula into a Functional of the Electron Density. Physical Review B 1988, 37, 785-789. (33) Miehlich, B.; Savin, A.; Stoll, H.; Preuss, H. Results Obtained with the Correlation Energy Density Functionals of Becke and Lee, Yang and Parr. Chem. Phys. Lett. 1989, 157, 200-206. (34) Stewart, J. J. P. Optimization of Parameters for Semiempirical Methods. I. Method. J. Comp. Chem. 1989; 10, 209-220. (35) Cancès, E.; Mennucci, B.; Tomasi, J. A New Integral Equation Formalism for the Polarizable Continuum Model: Theoretical Background and Applications to Isotropic and Anisotropic Dielectrics. J. Chem. Phys. 1997, 107, 3032-3041. (36) Mennucci, B.; Tomasi, J. Continuum Solvation Models: A New Approach to the Problem of Solute's Charge Distribution and Cavity Boundaries. J. Chem. Phys. 1997, 106, 5151-5158. (37) Mennucci, B.; Cancès, E.; Tomasi, J. Evaluation of Solvent Effects in Isotropic and Anisotropic Dielectrics and in Ionic Solutions with a Unified Integral Equation 31 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Method: Theoretical Bases, Computational Implementation, and Numerical Applications. J. Phys. Chem. B 1997, 101, 10506-10517. (38) Tomasi, J.; Mennucci, B.; Cancès, E. The IEF Version of the PCM Solvation Method: An Overview of a New Method Addressed to Study Molecular Solutes at the QM Ab Initio Level. J. Mol. Struct. (Theochem) 1999, 464, 211-226. (39) Llinas, P.; Stura, E. A.; Ménez, A.; Kiss, Z.; Stigbrand, T.; Millán, J. L.; Le Du, M. H. Structural Studies of Human Placental Alkaline Phosphatase in Complex with Functional Ligands. J. Mol. Biol. 2005, 350, 441-451. (40) Stec, B.; Holtz, K. M.; Kantrowitz, E. R. A revised Mechanism for the Alkaline Phosphatase Reaction Involving Three Metal Ions. J. Mol. Biol. 2000, 299, 1303-1311. (41) Jencks, W. P.; Regenstein, J. Ionization Constants of Acids and Bases. In Handbook of Biochemistry and Molecular Biology; Fasman, G. D., Ed.; CRC Press: Cleveland, OH, 1976; pp 305-351. (42) Sunden, F.; Peck, A.; Salzman, J.; Ressl, S.; Herschlag, D. Extensive SiteDirected Mutagenesis Reveals Interconnected Functional Units in the Alkaline Phosphatase Active Site. eLife 2015, 4, e06181. (43) Holtz, K. M.; Stec, B.; Kantrowitz, E. R. A Model of the Transition State in the Alkaline Phosphatase Reaction. J. Biol. Chem. 1999, 274, 8351-8354. (44) Duarte, F.; Amrein, B. A.; Kamerlin, S. C. L. Modeling Catalytic Promiscuity in the Alkaline Phosphatase Superfamily. Phys. Chem. Chem. Phys. 2013, 15, 11160-11177 and references therein. (45) Pabis, A.; Kamerlin, S. C. L. Promiscuity and Electrostatic Flexibility in the Alkaline Phosphatase Superfamily. Curr. Opin. Struct. Biol. 2016, 37, 14-21 and references therein. 32 ACS Paragon Plus Environment

Page 32 of 35

Page 33 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

(46) The classical transition-state theory expresses the rate constant for a reaction as k = (kBT/h)exp(-∆G/RT), in which kB is the Boltzmann constant, R is the gas constant, T is the absolute temperature, h is the Planck constant, and ∆G is the free energy of activation. (47) Ribeiro, A. J. M.; Alberto, M. E.; Ramos, M. J.; Fernandes, P. A.; Russo, N. The Catalytic Mechanism of Protein Phosphatase 5 Established by DFT Calculations. Chem. Eur. J. 2013, 19, 14081-14089 and references therein. (48) O´Brien, P. J.; Lassila, J. K.; Fenn, T. D.; Zalatan, J. G.; Herschlag, D. Arginine Coordination in Enzymatic Phosphoryl Transfer: Evaluation of the Effect of Arg166 Mutations in Escherichia Coli Alkaline Phosphatase. Biochemistry 2008, 47, 76637672. (49) Kirsch, J. F. In Advances in Linear Free Energy Relationships; Chapman, N. B., Shorter, J., Eds.; Plenum, New York, 1972, pp 369-400. (50) Hollfelder, F.; Herschlag, D. The Nature of the Transition State for EnzymeCatalyzed Phosphoryl Transfer. Hydrolysis of O-aryl Phosphorothioates by Alkaline Phosphatase. Biochemistry 1995, 34, 12255-12264. (51) Sham, Y.Y.; Chu, Z. T.; Warshel, A. Consistent Calculations of pKa’s of Ionizable Residues in Proteins: Semi-microscopic and Microscopic Approaches. J. Phys. Chem. B 1997, 101, 4458-4472. (52) Repič, M.; Purg, M.; Vianello, R.; Mavri, J. Examining Electrostatic Preorganization in Monoamine Oxidases A and B by Structural Comparison and pKa Calculations. J. Phys. Chem. B 2014, 118, 4326-4332.

33 ACS Paragon Plus Environment

Journal of Chemical Information and Modeling

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(53) Wu, X.; Lee, J.; Brooks, B. R. Origin of pKa Shifts of Internal Lysine Residues in SNase Studied Via Equal-Molar VMMS Simulations in Explicit Water. J. Phys. Chem. B 2017, in press, DOI: 10.1021/acs.jpcb.6b08249.

34 ACS Paragon Plus Environment

Page 34 of 35

Page 35 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

For Table of Contents Use Only

Quantum Mechanical Study on the Catalytic Mechanism of Alkaline Phosphatases

Gabriela L. Borosky*

35 ACS Paragon Plus Environment