Molecular Mechanism of ATP Hydrolysis in an ABC Transporter - ACS

Oct 5, 2018 - Hydrolysis of nucleoside triphosphate (NTP) plays a key role for the function of many biomolecular systems. However, the chemistry of th...
0 downloads 0 Views 5MB Size
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.

Research Article Cite This: ACS Cent. Sci. XXXX, XXX, XXX−XXX

http://pubs.acs.org/journal/acscii

Molecular Mechanism of ATP Hydrolysis in an ABC Transporter Marten Prieß,†,§ Hendrik Göddeke,†,§ Gerrit Groenhof,‡ and Lars V. Schäfer*,† †

Theoretical Chemistry, Faculty of Chemistry and Biochemistry, Ruhr-University Bochum, D-44780 Bochum, Germany Department of Chemistry and Nanoscience Center, University of Jyväskylä, P.O. Box 35, FI-40014 Jyväskylä, Finland



Downloaded via 188.72.96.95 on October 6, 2018 at 22:42:54 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

S Supporting Information *

ABSTRACT: Hydrolysis of nucleoside triphosphate (NTP) plays a key role for the function of many biomolecular systems. However, the chemistry of the catalytic reaction in terms of an atomic-level understanding of the structural, dynamic, and free energy changes associated with it often remains unknown. Here, we report the molecular mechanism of adenosine triphosphate (ATP) hydrolysis in the ATP-binding cassette (ABC) transporter BtuCD-F. Free energy profiles obtained from hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations show that the hydrolysis reaction proceeds in a stepwise manner. First, nucleophilic attack of an activated lytic water molecule at the ATP γ-phosphate yields ADP + HPO2− 4 as intermediate product. A conserved glutamate that is located very close to the γ-phosphate transiently accepts a proton and thus acts as catalytic base. In the second step, the proton is transferred back from the catalytic base to the γ-phosphate, yielding ADP + H2PO−4 . These two chemical reaction steps are followed by rearrangements of the hydrogen bond network and the coordination of the Mg2+ ion. The rate constant estimated from the computed free energy barriers is in very good agreement with experiments. The overall free energy change of the reaction is close to zero, suggesting that phosphate bond cleavage itself does not provide a power stroke for conformational changes. Instead, ATP binding is essential for tight dimerization of the nucleotide-binding domains and the transition of the transmembrane domains from inward- to outward-facing, whereas ATP hydrolysis resets the conformational cycle. The mechanism is likely relevant for all ABC transporters and might have implications also for other NTPases, as many residues involved in nucleotide binding and hydrolysis are strictly conserved.



INTRODUCTION Nucleoside triphosphate (NTP) hydrolysis is one of the most fundamental chemical reactions in biology.1,2 The free energy release upon cleavage of NTP into nucleoside diphosphate (NDP) and phosphate (Pi) is used by various proteins such as motors,3,4 signaling proteins,5 and transporters.6 The basis for the function of these proteins is their high catalytic activity: some NTPase enzymes accelerate hydrolysis by up to 9 orders of magnitude over the reaction in aqueous solution. However, despite its importance, a detailed atomic-level understanding of the molecular mechanisms that underlie phosphoryl transfer reactions is still partly lacking.1,2 This discrepancy can at least partly be attributed to experimental challenges linked to achieving the high spatial and temporal resolution required for studying this intricate chemical reaction in biomolecular systems. The hydrolysis reaction involves several key steps, including nucleophilic attack of a lytic water molecule at the γphosphate, cleavage of the scissile phosphate bond, proton transfer (PT) events, and hydrogen bond rearrangements. Theoretical methods such as hybrid quantum mechanical/ molecular mechanical (QM/MM)7 and ab initio molecular dynamics8 approaches can provide the lacking atomic-level insights into the structural and energetic details of chemical © XXXX American Chemical Society

reactions in complex environments. Hence, such methods have been applied to study phosphate hydrolysis in aqueous solution,9−15 guanosine triphosphate (GTP) hydrolysis in signal transduction proteins and translational GTPases (recently reviewed by Carvalho et al.16), and adenosine triphosphate (ATP) hydrolysis in molecular motors such as kinesin,17 myosin (recently reviewed by Kiani and Fischer18), actin,19,20 and F1-ATPase.21−24 These studies focused on the following questions. First, does hydrolysis proceed via an associative, dissociative, or concerted mechanism? These pathways are (i) idealized, and (ii) energetically very close, which has led to controversial discussions in the literature.2,15,16,25,26 Second, what is the nature of the protonaccepting group, and along which pathway(s) does proton transfer occur? A third question concerns the height of the activation energy barrier, and whether the hydrolysis reaction inside the protein is ex- or endergonic. The free energy of ATP hydrolysis is ΔG = −7.3 kcal/mol in aqueous solution under standard conditions, and even more exergonic under physiological conditions.27 However, ΔG of the hydrolysis Received: June 12, 2018

A

DOI: 10.1021/acscentsci.8b00369 ACS Cent. Sci. XXXX, XXX, XXX−XXX

Research Article

ACS Central Science

Figure 1. ATP hydrolysis mechanism in the ABC transporter BtuCD-F. (A−C) Snapshots from QM/MM MD simulations of the reactant state (ES), transition state (TS1), and product state (IS1) of the first reaction step. (D) Potential of mean force (PMF) along the Pγ−OW distance from the simulations with the standard and enlarged QM subsystems (black and gray lines, respectively); the larger QM subsystem additionally includes the side-chains of H191 and K39. (E) Distance between the transferred proton of the attacking water molecule, HW, and the Oϵ atom of E159. (F) Change of the breaking Pγ−Oβ bond distance (upper panel) and distance between the proton-accepting Oϵ atom of E159 and the transferred water proton (lower panel) along the reaction coordinate. The black and gray lines were obtained with the standard and enlarged QM subsystem, respectively. (G) Projection of the ES (black), TS1 (magenta), and IS1 (cyan) trajectories onto the Pγ−Oβ and Oϵ−OW distances.

possible mechanistic scenarios of ATP hydrolysis in ABC transporters.36 However, X-ray crystallography can only provide static snapshots and cannot resolve the positions of the protons. In addition, the energetics of the reaction, in terms of free energy profiles of the individual reaction steps, remains unknown. The link between several human diseases and the malfunction of ABC transporters6 further underlines the significance of understanding their mechanisms of action. To our knowledge, ATP hydrolysis in ABC transporters has been addressed in only three theoretical studies published so far. Figure S1 shows a schematic overview of the different mechanistic scenarios proposed. In the first study, Zhou and co-workers37 used the AM1 semiempirical quantum chemical method within a QM/MM framework to investigate ATP hydrolysis in the Escherichia coli hemolysin toxin transporter HlyB. On the basis of potential energy surface calculations, they concluded that the conserved H-loop histidine residue (H191 in Figure 1, the corresponding residue in HlyB is H622) is directly involved in the hydrolysis. According to Zhou and co-workers, the H-loop histidine first acts as a general acid by donating its NϵH proton to the γ-phosphate, and subsequently accepts a proton from the lytic water molecule.37 This proposed mechanism is surprising, as it proceeds via a negatively charged histidine imidazole ring, which at pH 7 is unfavorable in free energy. Indeed, the

reaction itself might be different in the protein environment. For example, ΔG is close to zero in F1-ATPase.28 None of the above questions have been comprehensively answered for ATP-binding cassette (ABC) transporters, which couple the free energy gain of ATP binding and hydrolysis to the transport of substrate molecules across biological membranes.6,29−32 All ABC transporters share the same architecture, a dimer of two nucleotide-binding domains (NBDs) or ATP-binding cassettes (ABCs) that bind and hydrolyze ATP, and two transmembrane domains (TMDs) that provide the pathway for substrate translocation across the membrane. The NBDs are usually highly conserved among different ABC transporters, suggesting a common ATP hydrolysis mechanism. ABC transporters can be understood as ATP-powered molecular machines. Their working cycle is based on the conversion of the chemical energy stored in ATP into mechanical work, which needs to be transmitted from the NBDs to the TMDs to induce the conformational changes necessary to translocate transport substrates. Whether the actual “power stroke”, in terms of the release of free energy that is necessary to drive the conformational changes, is provided by hydrolysis or binding of ATP remains an open question, the answer to which might depend on the type of ABC transporter.33−35 X-ray crystal structures, including some with transition state analogs, yielded valuable insights into B

DOI: 10.1021/acscentsci.8b00369 ACS Cent. Sci. XXXX, XXX, XXX−XXX

Research Article

ACS Central Science

phosphate cleaves the phosphate bond and yields ADP + HPO2− 4 as intermediate product. Upon formation of the new Pγ−O bond, the attacking water molecule transfers a proton to the highly conserved glutamate E159, which acts as a catalytic base. In the second reaction step, the proton is transferred back − from E159 to HPO2− 4 , leading to ADP + H2PO4 . Finally, in the third step, the hydrogen bond network rearranges such that H2PO−4 is stabilized by hydrogen bonds to ADP and E159; this rearrangement also involves changes of the Mg2+ ion coordination. From our computed free energy profiles, we estimate a reaction rate constant of ca. 7 s−1, in good agreement with the experimental value of ca. 0.5 s−1. The ΔG of the overall reaction is close to zero, ΔG = +1.8 kcal/mol. Hence, cleavage of the phosphate bond itself cannot provide the power stroke for conformational changes of the transporter that are linked to substrate transport. Rather, the free energy of ATP binding is essential for inducing the tight dimerization of the NBDs and the associated transition of the TMDs from an inward- to an outward-facing conformation, whereas ATP hydrolysis is required to reset the conformational cycle.51−53 Furthermore, although not explicitly investigated in this work, we speculate that unbinding of the hydrolysis products Pi and ADP likely plays an important role as well.

activation energy barrier of 22 kcal/mol obtained from the calculations37 is rather high. In addition, the highly conserved Walker B glutamate (E159 in Figure 1), mutations of which impair ATPase activity and which is thought to be the catalytic base,36,38−40 is not directly involved as a general base in this proposed mechanism. The other two studies,41,42 which investigated the maltose transporter MalK, also suggest a direct involvement of the H-loop histidine (H192 in MalK). Huang and Liao employed a QM/MM nudged elastic band method to calculate energy profiles, using the B3LYP density functional for the QM region.41 Assuming a positively charged (doubly protonated) H-loop histidine, they propose that two protons are transferred upon nucleophilic water attack at the γphosphate, one from H192 to the γ-phosphate and the second from the lytic water to the conserved glutamate. These two proton transfers occur at the same point along the chosen reaction coordinate, and the reaction proceeds via a single barrier of 19 kcal/mol. These findings were confirmed by Hsu and co-workers,42 who carried out QM/MM metadynamics simulations using the BLYP density functional; however, the reported reaction barrier of 10.5 kcal/mol is rather low. Furthermore, in this mechanism the protonation states of the H-loop histidine and the conserved glutamate at the end of the reaction differ from the beginning, and resetting the catalytic cycle would thus require additional protonation/deprotonation events. Another peculiarity of this second proposed mechanism is that a doubly protonated H-loop histidine might only be populated at pH < 7. However, at such low pH, ATP could be protonated as well,16 which would have implications for the mechanism. In summary, the mechanism of ATP hydrolysis in ABC transporters remains elusive. Here, we use QM/MM molecular dynamics (MD) simulations to investigate the structural and energetic details of the ATP hydrolysis reaction in the vitamin B12 importer BtuCD-F from E. coli, one of the structurally most thoroughly characterized ABC transporters.43−47 Unlike most previous QM/MM studies, the ATP hydrolysis mechanism was not investigated with minimum energy path calculations, but on the basis of potentials of mean force (PMFs) at room temperature. The MD sampling necessary for obtaining statistically meaningful PMFs makes this approach computationally much more expensive (the accumulated sampling time of our free energy simulations is more than 1 ns, which with the 1 fs time step used to integrate the equations of motion in the MD corresponds to more than 1 million electronic structure calculations), but it takes entropic contributions explicitly into account. Previous studies have shown that entropy can play a significant role for the free energy profile of NTP hydrolysis reactions. 18,48 In the QM/MM MD simulations, the MPW1B95 density functional49 with the 631+G(d,p) basis set was used for the QM part, which consists of the triphosphate chain, Mg2+ ion, and several protein sidechains and water molecules (see the Methods section). MPW1B95 was chosen, because it can very accurately describe the energetics of the hydrolysis reaction, as shown in a benchmark study of phosphodiester hydrolysis that compared 52 density functionals.50 When used together with a basis set that includes diffuse functions on heavy atoms, MPW1B95 energies deviate by less than 2 kcal/mol from approximated CCSD(T)/CBS values.50 Our results show that ATP hydrolysis in the ABC transporter follows a three-step mechanism. In the first step, nucleophilic attack of the lytic water molecule at the γ-



RESULTS AND DISCUSSION Attack of the Lytic Water Molecule. The first step of the hydrolysis reaction is the nucleophilic attack of the lytic water molecule at the γ-phosphate of ATP and breakage of the scissile Pγ−Oβ bond. As a prerequisite, a water molecule has to be present close to the γ-phosphate and adopt a near-linear attack angle (see below). In the X-ray structure,46 the positions of water molecules are not resolved. However, in our force field MD simulations, we found a single water molecule to repeatedly and reversibly adopt such a bridging position between the γ-phosphate and E159, the carboxylate side-chain of which is located very close to Pγ (Figure 1). We find a probability of ca. 15% for such a hydrolysis-competent water configuration (Table S1 and Figure S2 in the Supporting Information). In other words, water molecules frequently exchange between hydrolysis-competent and hydrolysis-noncompetent configurations, but most of the time the transporter is not in a hydrolysis-competent state. We initiated our QM/ MM simulations from a representative hydrolysis-competent snapshot taken from our previous force field MD simulations.54 Figure 1 summarizes the results of our QM/MM MD simulations of the first hydrolysis step. Configurations that are representative for the reactant (ES), transition (TS1), and first intermediate product (IS1) state are shown in Figure 1A−C, and the PMF is shown in Figure 1D. The reaction proceeds from the ES at a Pγ−OW distance of 2.9 Å via a free energy maximum (denoted as TS1) at 2.0 Å to the first product state (IS1), in which the Pγ−OW distance is 1.7 Å; i.e., the formation of the new bond is completed. The products of this first reaction step are ADP + HPO2− 4 . The activation free energy is 14.1 kcal/mol, and IS1 lies 8.8 kcal/mol above the ES. These findings are confirmed by additional PMF simulations with a larger QM subsystem that also includes the side-chains of K39 and H191 (Figure 1D, gray line) and by PMF simulations in which the Pγ−Pβ distance was used as reaction coordinate (Figures S3 and S4 in the Supporting Information). In the ES, the attacking water molecule forms hydrogen bonds with the E159 carboxylate group and the backbone of S163. These hydrogen bonds not only polarize the water C

DOI: 10.1021/acscentsci.8b00369 ACS Cent. Sci. XXXX, XXX, XXX−XXX

Research Article

ACS Central Science molecule and thus activate it for proton abstraction and nucleophilic attack, but also hold it in line with the γphosphate (OW−Pγ−Oβ angle of about 170°). A similar geometry of the hydrolysis-competent state was proposed previously on the basis of force field MD simulations of the ABC exporter Sav1866.55 Upon the attack of the lytic water at Pγ, the Pγ−Oβ bond breaks, and a proton is transferred from the water molecule to E159, which acts as catalytic base. This proton transfer (PT) is not imposed by the chosen RC, which is the Pγ−OW distance. Transfer of the water proton to the Oϵ atom of E159 occurs spontaneously (i.e., without any bias) at Pγ−OW distances ≤2.0 Å, i.e., at the TS and beyond (Figure 1E). Figure 1F shows that the concerted cleavage of the Pγ−Oβ phosphate bond and PT to E159 occur in a largely synchronous fashion in the Pγ−OW distance range 1.8−2.1 Å. For larger distances, the Pγ−Oβ bond is elongated but not yet completely broken. These two modes are thus tightly coupled, underlining the importance of PT for phosphate bond cleavage. In Figure 1G, the ES, TS1, and IS1 trajectories are projected onto the configuration space spanned by the Pγ−Oβ distance and the distance between the E159 Oϵ atom and the water oxygen. At the TS, a broad range of Pγ−Oβ distances is sampled (magenta curve), such that the corresponding distance distribution almost overlaps with those of ES and IS1, which are much more narrow. This demonstrates that, at neighboring points along the chosen RC, there is overlap of the sampled distributions along these degrees of freedom, which are not part of the chosen RC but are of crucial importance for the reaction. Furthermore, Figure 1G shows that the Oϵ−OW distance is tightened at TS1, indicating that the transition state is stabilized by a strengthening of the hydrogen bond between E159 and the attacking water molecule. This observation agrees with the interpretation of X-ray crystal structures of the ABC transporter MalK trapped with transition state analogs.36 In addition to the attack angle, also the orientation of the lytic water molecule with respect to the γ-phosphate is important. To enable P−O bond formation, the water oxygen needs to point toward Pγ. We observed the reorientation of the water dipole in our simulations (Figure 2). Upon an increase in the Pγ−OW distance from 2.9 to 3.3 Å, the free energy increases by about 1.6 kcal/mol (Figure 1D). This increase is due to the rearrangement of hydrogen bonds of the attacking water molecule. The hydrogen bond to the backbone of S163, which is present in the hydrolysis-competent geometry at the ES, is broken and replaced by a hydrogen bond with an Oγ atom of ATP (Figure 2A). These rearrangements occur at Pγ− OW distances beyond 3.3 Å; see Figure 2B for a representative trajectory. At smaller distances, the water molecule preferentially adopts the hydrolysis-competent orientation. This was confirmed by additional control simulations in which the PMF was calculated backward from 3.3 to 2.9 Å; no hysteresis was found. The hydrogen bond between the lytic water molecule and the carbonyl group of the S163 backbone, which is present in the hydrolysis-competent state, is broken in the noncompetent state. As a consequence, the carbonyl group of S163 contacts a Hα atom of G130 (Figure 2A), suggesting that this interaction, which is likely very weak, has to be broken to stabilize the attacking water molecule in the hydrolysis-competent orientation. Deprotonation of the Catalytic Base. After the first reaction step, E159 is protonated and hydrogen-bonded to the

Figure 2. Reorientation of the attacking water molecule. (A) Superposition of representative structures of the hydrolysis-competent (thick) and hydrolysis-noncompetent (transparent) states. (B) At a Pγ−OW distance of 3.3 Å, the water molecule reorients and rearranges its hydrogen bond network. The hydrogen bond between the catalytic glutamate and the water molecule is maintained.

OH group of HPO2− 4 (Figure 1C). This IS1 intermediate is further stabilized by hydrogen bonds of HPO2− 4 with S163, H191, and K39. However, the reference pKa values of the carboxylic acid and phosphoric acid in water (ca. 4.5 and 7.2, respectively) suggest deprotonated E159 and H2PO−4 as likely final reaction products. Back-transfer of the proton from E159 to the γ-phosphate requires the rearrangement of hydrogen bonds, because the protonated OH group of the γ-phosphate cannot accept a second proton. We thus considered two Oγ atoms as possible proton acceptors: the O1γ atom that coordinates the Mg2+ ion and the O3γ atom that is hydrogen-bonded to H191 (Figure 3); the O2γ atom is pointing away from E159 and is at a too large distance for PT. Figure 3 shows that PT from E159 to O1γ proceeds via a second barrier of 6.5 kcal/mol (i.e., TS2 is 15.3 kcal/mol above the initial ES) and leads to another metastable intermediate state (IS2) that is 11.8 kcal/mol above the initial ES. In IS2, the formed H2PO−4 anion is stabilized by hydrogen bonds with H191 and E159, and the coordinate bond between the protonated O1γ atom and Mg2+ is elongated (Figure 3E), priming it for dissociation (see below). The finding that the O1γ atom that coordinates Mg2+ accepts the proton from E159 might seem counterintuitive at first. However, recent NMR experiments on a GTPase showed that the O atom that coordinates the Mg2+ ion has the largest electron density,56 suggesting that it is the strongest base in the PT step. The alternative PT to the O3γ atom, which is hydrogen-bonded to H191, involves a higher activation free energy barrier and is thus unlikely because of a larger distance over which the proton needs to be transferred, and because of structural strain (Figure S5 in the Supporting Information). The described PT resets the catalytic cycle in terms of the protonation states of the involved amino acid side-chains. However, the free energy of IS2 is still rather high, +11.8 kcal/ D

DOI: 10.1021/acscentsci.8b00369 ACS Cent. Sci. XXXX, XXX, XXX−XXX

Research Article

ACS Central Science

− Figure 3. Proton transfer from E159 to HPO2− 4 and formation of H2PO4 + ADP. (A) Reactant, (B) transition, and (C) product states. (D) PMF along the distance between the E159 proton and O1γ as RC. (E) Distance of the Mg2+−O1γ coordinate bond during simulations of the reactant (IS1, black line) and product (IS2, cyan line) states.

Figure 4. Final rearrangement of hydrogen bond network and Mg2+ coordination. (A) Reactant, (B) transition, and (C) product states. (D) PMF along the distance between O2γ and Mg2+ as RC. (E) Distance between the proton on H2PO−4 and ADP during simulation of the product state.

The final product, H2PO−4 + ADP, is +1.8 kcal/mol higher in free energy than the initial reactant state; i.e., ΔG is close to zero. This suggests that (i) the reaction is reversible up to this point, and (ii) phosphate bond cleavage does not provide a power stroke, in the sense that the free energy change that is directly associated with it cannot trigger the (partial) opening of the NBD dimer. However, one should bear in mind that, in the final product state reached in our simulations, H2PO−4 and ADP are still bound to the transporter and coordinate the Mg2+ ion. Unbinding of these hydrolysis products involves reorganization of the Mg2+ coordination and individual solvation of the separated molecules. This is expected to be associated with a further lowering of the free energy. Our results thus suggest that release and solvation of ADP + Pi are essential for the exergonic nature of the overall reaction, a hypothesis that will be tested in future work. Our previous

mol, suggesting that further rearrangements are necessary to yield the final product. Indeed, Figure 4 shows that rotation of H2PO−4 , such that the O2γ atom is replacing the protonated O1γ atom in the Mg2+ coordination shell, has a low activation free energy barrier of 3.2 kcal/mol and yields a final PS that is favorable in free energy, only +1.8 kcal/mol above the initial ES. We assign the low activation barrier to the weakening of the Mg2+−O1γ bond due to the previous PT step. In the final product configuration, the hydrogen bond network rearranges such that H2PO−4 forms hydrogen bonds with ADP, E159, and H191 (Figure 4C). A proton is shared between H2PO−4 and ADP, hopping back and forth between the two moieties in the PS simulations (Figure 4E). The populations of H2PO−4 + ADP and HPO2− 4 + H-ADP are about 75% and 25%, respectively, suggesting a free energy difference of about 1 kBT and a very low barrier between these two states. E

DOI: 10.1021/acscentsci.8b00369 ACS Cent. Sci. XXXX, XXX, XXX−XXX

Research Article

ACS Central Science work54 showed that H191 acts not only as a linchpin36 to hold the ATP molecule in place, but also as a gatekeeper that regulates water access to the nucleotide-binding pocket. After breaking of the hydrogen bond to the γ-phosphate, H191 can adopt an outward-pointing conformation that is similar to the one observed in X-ray crystal structures of nucleotide-free BtuCD.43,44 In this conformation, the occupancy of the binding pocket with water molecules is significantly enhanced.54 This increased hydration of the binding pocket could be one of the initial steps on the way to the release of the reaction products, as water molecules can replace hydrogen bonds of ADP and Pi with the protein. Two-Water Mechanism is Unfavorable. Our results show that, in BtuCD, a single water molecule can hydrolyze ATP, thereby transiently transferring a proton to the catalytic glutamate in a direct manner, i.e., without the help of other water molecules or amino acid side-chains. This is possible because the carboxylate group of the glutamate is located very close to the γ-phosphate. This structural feature is found in all ABC transporters and discriminates them from other NTPases, where PT along longer water wires has been suggested.16,18,19,24,57,58 To investigate this possibility in BtuCD, we analyzed our force field MD simulations for the formation of a two-water bridge between the γ-phosphate of ATP and E159. Such a conformation was indeed observed (Figure S2 and Table S2 in the Supporting Information); its probability of ca. 14% is similar to that of the reactive conformation with one bridging water molecule (15%, see above). In the two-water bridged configuration, the E159 side-chain is oriented away from the γ-phosphate and at a larger distance from it (the mean distance between the Cδ atom of E159 and Pγ is 5.7 and 7.3 Å for the one- and two-water bridged conformations, respectively). Starting from this two-water bridged ES, we calculated the PMF for nucleophilic water attack on Pγ, using the Pγ−OW1 distance as RC and including the second water molecule in the QM subsystem. However, unlike for the one-water mechanism described above, we did not observe PT from the lytic water molecule (WAT1) to the second, assisting water (WAT2), suggesting a barrier along this degree of freedom that cannot be overcome in the short QM/ MM simulations. Therefore, we calculated a two-dimensional PMF that additionally includes the intramolecular OW−HW distance of WAT1. Figure S6 (Supporting Information) shows that PT from WAT1 via WAT2 onto E159 occurs at a Pγ−OW1 distance of 1.8 Å. The free energy barrier associated with this process is 24.9 kcal/mol and thus much higher than for the one-water mechanism. The intermediate formed by this reaction, HPO2− 4 + ADP, is formally the same as previously but lies higher in free energy (22.8 kcal/mol). We assign this difference to a different hydrogen bond network (compare Figure 3A to Figure S6). The two-water mechanism involves a significantly higher free energy barrier than the one-water mechanism (24.9 versus 15.3 kcal/mol) and is thus considered unlikely. This finding can explain the reduced ATPase activity of ATP-binding cassettes in which the catalytic glutamate is replaced by an aspartate.59−61 This is an intriguing and up to now unexplained finding, because the E-to-D mutation does not remove the carboxylate moiety that is required as a catalytic base but merely positions it slightly further away from the γ-phosphate. Heterodimeric ABC transporters, which feature an ATPase active consensus site and an inactive degenerate site, have an aspartate residue in the latter. However, the degenerate site in

ABC heterodimers bears additional noncanonical residues in the ABC signature motif, which might also contribute to the strongly reduced ATPase activity. The above results suggest that the shorter aspartate sidechain requires a two-water bridge for nucleophilic attack of the lytic water molecule and proton shuttling to the carboxylate group. We speculate that, in this case, hydrolysis has to proceed via a mechanism that involves a higher free energy barrier and thus a slower rate. We carried out additional force field MD simulations of the E159D mutant of BtuCD-F. Figure S7 and Table S3 in the Supporting Information show that, indeed, configurations with two bridging water molecules between D159 and Pγ are frequently observed. However, the orientation of the aspartate carboxylate group and the precise arrangement of the two water molecules differ from the wildtype (E159), which might alter the actual activation energy barrier. Further studies of this topic will be the subject of future work. Comparison to Experiments. The proposed mechanism is in agreement with, and can provide atomistic explanations for, available experimental data. The reaction rate estimated from the computed free energy profiles can be compared to experimental data for wild-type BtuCD-F from Locher and coworkers,46,62 who reported a turnover number of about 1 ATP per second, or 0.5 s−1 per NBD monomer. Under the assumption that ATP hydrolysis is the rate-determining step,63,64 this apparent experimental rate constant can be compared to our computations. We estimate a rate constant of ca. 7 s−1 from the highest computed activation free energy barrier of 15.3 kcal/mol, as obtained for the second reaction step (Figure 3), using the Eyring equation k = pwkBTh−1 exp[−ΔG‡/(RT)], where pw = 15% is the probability to actually find a hydrolysis-competent water molecule close to Pγ (see above). Within the statistical uncertainty in ΔG‡ of about ±0.8 kcal/mol (Figure 3D), the computational rate constants are between 2 and 27 s−1. Given that the rate constant is exponentially sensitive to the barrier height, the agreement between our simulations and experiment is very encouraging, although maybe even somewhat fortuitous. Additional strong support for our proposed mechanism with the Walker B glutamate as catalytic base comes from sitedirected mutagenesis studies of various ATP-binding cassettes.39,40,59,60,65 Furthermore, several studies report on the pH dependence of ATPase activity. For example, Zaitseva and co-workers63 found that the ATPase activities of the wild-type ABC dimer module of HlyB and the E-to-Q mutant, in which the Walker B glutamate is replaced by a glutamine, have similar pH dependencies with an optimum activity at around pH 7. At first sight, this finding seems to speak against the involvement of the glutamate as catalytic base. However, in a different ABC module (the TAP1 dimer), an aspartate-to-glutamate mutation in the catalytic dyad does alter the pH profile and increases ATPase activity,61,66 underlining the important role of the Walker B glutamate. Another study,67 which reports the pH dependence of the ATPase activity of the ABC transporter Pdr5, showed that mutation of the catalytic glutamate abolished ATPase activity, whereas mutation of the switch histidine had no strong effect, in agreement with our mechanism. Interestingly, Pdr5 has an optimal ATPase activity at pH 9,67 clearly speaking against an essential functional role of a doubly protonated switch histidine. However, the detailed mechanistic interpretation of pH-dependent ATPase essays is complicated. For example, the NBD dimerization equilibrium F

DOI: 10.1021/acscentsci.8b00369 ACS Cent. Sci. XXXX, XXX, XXX−XXX

Research Article

ACS Central Science

quantum mechanically (see below), the force field parameters of Meagher and co-workers76 were used. The protonation states of titratable groups were assigned with PROPKA77,78 in the presence of ATP. For the X-ray crystal structure, these computations yielded a pKa of 5.9 for the H-loop histidine; using 500 snapshots from previous 500 ns force field MD simulations of the ATP-bound state54 yielded an average pKa of 5.2. H191 was thus modeled as singly protonated at the Nϵ atom (the cytoplasmic pH of E. coli is 7.2−7.879). In additional control simulations with doubly protonated H191, the imidazole ring rotated and formed a hydrogen bond between the (protonated) Nδ−H and E159 in addition to the one between N ϵ −H and the ATP γ-phosphate. Such a conformation is not found in X-ray structures.36 For the QM/MM simulations, we used the standard QM/ MM functionality of GROMACS (version 4.5.7)80 in combination with a script that interfaces to GAUSSIAN09.81 While the QM/MM routines have been an integral part of the GROMACS MD program since version 3.3, the interface script is available in the Supporting Information or for download at wwwuser.gwdg.de/~ggroenh/qmmm.html. The MPW1B95 density functional49 was used together with a D3 dispersion correction82 and the 6-31+G(d,p) basis set. To validate this quantum chemical method, we recalculated the reaction profile for dimethylphosphate hydrolysis, as described in the benchmark study of Ribeiro and co-workers.50 The results (Table S4 in the Supporting Information) show that MPW1B95/631+G(d,p) energies closely agree with CCSD(T)/CBS values, with mean unsigned errors below 2 kcal/mol. In our QM/MM simulations of the ABC transporter, the QM subsystem consisted of the triphosphate chain of ATP; the Mg2+ ion and its two coordinating water molecules; the side-chains of S40, Q80, and E159; and one (or two) attacking water molecule(s), yielding 49 (52) QM atoms in total (Figure 1; an example GAUSSIAN09 input file is contained in Supporting Information). Additional control simulations were carried out with a larger QM region that also contained the H191 and K39 side-chains (66 QM atoms). In terms of the water molecules close to the γ-phosphate, the chosen QM subsystem was complete, as there were no additional water molecules present in this ATP-binding site.54 Electrostatic embedding of the QM subsystem into the MM point-charge surrounding83 was used, and hydrogen link atoms were used to saturate the QM subsystem at cuts through covalent bonds across the QM/MM boundary (only C−C single bonds were cut). The partial charges of the MM atoms of the CHn group next to the hydrogen link atom were set to zero, and the resulting (small) charge difference was distributed over the neighboring MM atoms such that the total charge was unchanged. Prior to MD simulation, all systems were energy-minimized (60 steepest descent steps). In the QM/MM MD simulations, the equations of motion were integrated with 1 fs time steps. An SCF convergence criterion of 10−8 Hartree was applied. For the standard QM subsystem (49 atoms), a single MD step took 6 min on a 2 × 10 core Xeon E5-2640 2.4 GHz node. NpT ensembles were simulated by coupling to a temperature bath at 300 K with a velocity rescaling84 thermostat (τT = 0.1 ps). For constant pressure, semi-isotropic coupling was applied by separately coupling the lateral (x, y) and normal (z) directions of the periodic simulation cell to a pressure bath at 1 bar using a Berendsen barostat (τp = 2.0 ps) and compressibility 4.5 × 10−5 bar−1. All nonbonded interactions were truncated at a large cutoff of 3.5 nm, and the nonbonded pair-list was

is also pH-dependent, and this can affect the measured ATPase activities.66



CONCLUSIONS The QM/MM free energy simulations presented in this work provide key insights into the catalytic strategy of ABC transporters. The ATP hydrolysis reaction proceeds in three steps. First, a polarized lytic water molecule attacks the ATP γphosphate, whichaccompanied by proton transfer to the catalytic glutamateleads to cleavage of the scissile Pγ−Oβ bond and formation of ADP + HPO2− 4 . In the second step, the proton is transferred back from the catalytic glutamate to the γphosphate, yielding ADP + H2PO−4 and resetting the catalytic cycle in terms of the protonation states of the involved amino acids. Third, these two reaction steps are followed by hydrogen bond rearrangements, which also involve transient changes in the coordination of the Mg2+ ion. The rate constant estimated from the computed activation free energy barrier agrees with experiments. The described mechanism explains the catalytic activity by providing atomic-level details of the structure, dynamics, and energetics underlying the ATP hydrolysis reaction. The nucleotide-binding domains are highly conserved in sequence and structure among all ABC transporters, 6 and the mechanism is thus very likely relevant not only for BtuCD (or other type II importers), but also for a broad range of ABC transporters. However, although many residues involved in nucleotide binding and/or hydrolysis are conserved even across different NTPases,5,6,18,68,69 the general mechanistic picture seems to be somewhat more diverse. For example, kinesin, myosin, and F1-ATPase share a similar ATP hydrolysis mechanism, in which a second assisting water molecule is involved in shuttling the proton from the lytic water molecule to the catalytic base.70 In GTPases, diverse mechanisms have been proposed,16 including substrate-assisted catalysis mechanisms in which the substrate itself acts as a base. These mechanisms differ from the one proposed here for ATP hydrolysis in ABC transporters. One specific feature that distinguishes ABC transporters from the previous systems is the close distance of the catalytic glutamate to the ATP γphosphate, and longer proton relay pathways are thus not required. We conclude that some, but not all, aspects of the hydrolysis reaction mechanism described here might be relevant also for other proteins, thus contributing to our understanding of the chemomechanical energy conversion in NTP-driven molecular machines.



METHODS The QM/MM MD simulations were started from preequilibrated snapshots taken from our recent force field MD study,54 and we hence refer to this previous study for the details of the simulation setup. In brief, the complete BtuCD-F complex was simulated in a fully solvated palmitoyloleoylphosphatidylcholine (POPC) lipid bilayer. The simulations were initiated from a nucleotide-bound X-ray crystal structure (PDB entry 4FI346) after replacing 5′-(β,γ-imido)triphosphate (AMP-PNP) by ATP and changing glutamine at position 159 to glutamate and cysteine at position 162 to asparagine to create the wild-type protein. The Amber ff99SB-ILDN force field71,72 was used, in combination with the Berger lipid parameters73,74 and the TIP4P-Ew75 water model. For the ATP molecule in the second binding site that was not treated G

DOI: 10.1021/acscentsci.8b00369 ACS Cent. Sci. XXXX, XXX, XXX−XXX

ACS Central Science



updated at every time step. The bond lengths in the MM subsystem were constrained using LINCS,85 and SETTLE86 was used to constrain the internal degrees of freedom of the MM water molecules. No bond-length constraints were applied to the bonds in the QM subsystem. For the PMFs, three sets of free energy simulations were carried out. In the first set (reaction step 1), we chose the distance between the oxygen atom of the lytic water molecule, OW, and the Pγ atom of ATP as reaction coordinate (RC). To test the influence of this choice, we carried out additional PMF simulations with the Pγ−Pβ distance as RC, with very similar results (Figures S3 and S4 in the Supporting Information). In the second set (reaction step 2), the distance between the proton on E159 and an O atom of the γ-phosphate was chosen as RC. In the third set (reaction step 3), the distance between an Oγ atom and the Mg2+ ion was chosen as RC. In addition, for the alternative reaction mechanism involving two water molecules, a two-dimensional PMF was calculated. In all cases, the RC was treated with a constraint, and the force acting on the constraint was recorded in the course of the MD simulation (by using the pull code of GROMACS). Subsequent structures along the RC to initialize the MD simulations were generated by incrementing the constraint distance (a spacing of 0.1 Å was used unless otherwise noted; additional points were introduced where necessary), followed by 60 steps steepest descent energy minimization. Then, 4 ps of QM/MM MD sampling (2.5 ps for the larger QM subsystem) at 300 K was carried out at each value of the RC. Finally, the PMF can be obtained by integrating over the mean force, PMF = ∫ ⟨fc⟩r dr − kBT ln ⟨z(r)−1/2⟩, where r is the variable along the path, ⟨fc⟩r is the ensemble average of the constraint force at a particular point of the RC, and z(r) is the Fixman determinant of the coordinate transformation used,87,88 which in principle has to be taken into account. However, Schlitter and co-workers 13 have shown for phosphoester hydrolysis in water that the Fixman correction due to constraining two degrees of freedom, one of which was the Pγ−OW distance, is less than 1 kcal/mol. In our case, the correction to the PMF is even smaller because (usually) only one degree of freedom was constrained. As this correction is smaller than the statistical error, it was neglected. For the twodimensional PMF, the forces acting on the two constrained coordinates were recorded, and the PMF in the subspace of these two coordinates was obtained as described previously.13 The first 2 ps (first 1.25 ps for the larger QM subsystem) at each distance were discarded from the analysis of fc (Figure S8 in the Supporting Information). Statistical errors were estimated using a block averaging procedure,89 as implemented in the gmx analyze tool of GROMACS.



Research Article

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Phone: +49 234 3221582. Fax: +49 234 3214045. ORCID

Lars V. Schäfer: 0000-0002-8498-3061 Author Contributions §

M.P. and H.G. contributed equally to this work.

Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We thank Rachel Glaves, Dominik Marx, Nisanth Nair, Patrick Nuernberger, Lutz Schmitt, and Ravi Tripathi for insightful discussions, and António Ribeiro and Pedro Fernandes for sharing G03 output files from ref 50. The CSC−IT Center for Science, Finland, and Steinbuch Centre for Computing (SCC), Karlsruhe/Germany, are acknowledged for providing computational resources. This work was funded by the Deutsche Forschungsgemeinschaft (DFG) through an Emmy Noether grant to L.V.S. (SCHA 1574/3-1) and Cluster of Excellence RESOLV (EXC 1069). G.G. acknowledges the Academy of Finland for support (Grant 304455).



REFERENCES

(1) Lassila, J. K.; Zalatan, J. G.; Herschlag, D. Biological phosphoryltransfer reactions: understanding mechanism and catalysis. Annu. Rev. Biochem. 2011, 80, 669−702. (2) Kamerlin, S. C. L.; Sharma, P. K.; Prasad, R. B.; Warshel, A. Why nature really chose phosphate. Q. Rev. Biophys. 2013, 46, 1−132. (3) Howard, J. Molecular motors: structural adaptations to cellular functions. Nature 1997, 389, 561−567. (4) Kinosita, K., Jr; Adachi, K.; Itoh, H. Rotation of F1-ATPase: how an ATP-driven molecular machine may work. Annu. Rev. Biophys. Biomol. Struct. 2004, 33, 245−268. (5) Sprang, S. R. G protein mechanisms: insights from structural analysis. Annu. Rev. Biochem. 1997, 66, 639−678. (6) Davidson, A. L.; Dassa, E.; Orelle, C.; Chen, J. Structure, function, and evolution of bacterial ATP-binding cassette systems. Microbiol. Mol. Biol. Rev. 2008, 72, 317−364. (7) Senn, H. M.; Thiel, W. QM/MM methods for biomolecular systems. Angew. Chem., Int. Ed. 2009, 48, 1198−1229. (8) Marx, D.; Hutter, J. Ab initio Molecular Dynamics: Basic Theory and Advanced Methods; Cambridge University Press: New York, 2009. (9) Akola, J.; Jones, R. O. ATP hydrolysis in water − a density functional study. J. Phys. Chem. B 2003, 107, 11774−11783. (10) Klähn, M.; Rosta, E.; Warshel, A. On the mechanism of hydrolysis of phosphate monoesters dianions in solutions and proteins. J. Am. Chem. Soc. 2006, 128, 15310−15323. (11) Rosta, E.; Kamerlin, S. C. L.; Warshel, A. On the interpretation of the observed linear free energy relationship in phosphate hydrolysis: a thorough computational study of phosphate diester hydrolysis in solution. Biochemistry 2008, 47, 3725−3735. (12) Glaves, R.; Mathias, G.; Marx, D. Mechanistic insights into the hydrolysis of a nucleoside triphosphate model in neutral and acidic solution. J. Am. Chem. Soc. 2012, 134, 6995−7000. (13) Li, W.; Rudack, T.; Gerwert, K.; Gräter, F.; Schlitter, J. Exploring the multidimensional free energy surface of phosphoester hydrolysis with constrained QM/MM dynamics. J. Chem. Theory Comput. 2012, 8, 3596−3604. (14) Harrison, C. B.; Schulten, K. Quantum and classical dynamics simulations of ATP hydrolysis in solution. J. Chem. Theory Comput. 2012, 8, 2328−2335.

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acscentsci.8b00369. Additional data and figures including schematic, force field MD simulations, comparison of PMFs and reaction mechanisms, 2D-PMF of first step of ATP hydrolysis, energies for dimethylphosphate hydrolysis, and equilibration during QM/MM simulation (PDF) Example GAUSSIAN09 input file (ZIP) QM/MM bash script (ZIP) H

DOI: 10.1021/acscentsci.8b00369 ACS Cent. Sci. XXXX, XXX, XXX−XXX

Research Article

ACS Central Science (15) Prasad, B. R.; Plotnikov, N. V.; Warshel, A. Addressing open questions about phosphate hydrolysis pathways by careful free energy mapping. J. Phys. Chem. B 2013, 117, 153−163. (16) Carvalho, A. T. P.; Szeler, K.; Vavitsas, K.; Åqvist, J.; Kamerlin, S. C. L. Modeling the mechanisms of biological GTP hydrolysis. Arch. Biochem. Biophys. 2015, 582, 80−90. (17) McGrath, M. J.; Kuo, I.-F. W.; Hayashi, S.; Takada, S. Adenosine triphosphate hydrolysis mechanism in kinesin studied by combined quantum-mechanical/molecularmechanical metadynamics simulations. J. Am. Chem. Soc. 2013, 135, 8908−8919. (18) Kiani, F. A.; Fischer, S. Advances in quantum simulations of ATPase catalysis in the myosin motor. Curr. Opin. Struct. Biol. 2015, 31, 115−123. (19) McCullagh, M.; Saunders, M. G.; Voth, G. A. Unraveling the mystery of ATP hydrolysis in actin filaments. J. Am. Chem. Soc. 2014, 136, 13053−13058. (20) Sun, R.; Sode, O.; Dama, J. F.; Voth, G. A. Simulating protein mediated hydrolysis of ATP and other nucleoside triphosphates by combining QM/MM molecular dynamics with advances in metadynamics. J. Chem. Theory Comput. 2017, 13, 2332−2341. (21) Strajbl, M.; Shurki, A.; Warshel, A. Converting conformational changes to electrostatic energy in molecular motors: the energetics of ATP synthase. Proc. Natl. Acad. Sci. U. S. A. 2003, 100, 14834−14839. (22) Dittrich, M.; Hayashi, S.; Schulten, K. On the mechanism of ATP hydrolysis in F1- ATPase. Biophys. J. 2003, 85, 2253−2266. (23) Beke-Somfai, T.; Lincoln, P.; Nordén, B. Double-lock ratchet mechanism revealing the role of αSER-344 in F0F1 ATP synthase. Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 4828−4833. (24) Hayashi, S.; Ueno, H.; Shaikh, A. R.; Umemura, M.; Kamiya, M.; Ito, Y.; Ikeguchi, M.; Komoriya, Y.; Iino, R.; Noji, H. Molecular mechanism of ATP hydrolysis in F1-ATPase revealed by molecular simulations and single-molecule observations. J. Am. Chem. Soc. 2012, 134, 8447−8454. (25) Yang, Y.; Cui, Q. The hydrolysis activity of adenosine triphosphate in myosin: a theoretical analysis of anomeric effects and the nature of the transition state. J. Phys. Chem. A 2009, 113, 12439−12446. (26) Duarte, F.; Åqvist, J.; Williams, N. H.; Kamerlin, S. C. L. Resolving apparent conflicts between theoretical and experimental models of phosphate monoester hydrolysis. J. Am. Chem. Soc. 2015, 137, 1081−1093. (27) Rosing, J.; Slater, E. C. The value of ΔG0 for the hydrolysis of ATP. Biochim. Biophys. Acta, Bioenerg. 1971, 267, 275−290. (28) Watanabe, R.; Iino, R.; Noji, H. Phosphate release in F1ATPase catalytic cycle follows ADP release. Nat. Chem. Biol. 2010, 6, 814−820. (29) Rees, D. C.; Johnson, E.; Lewinson, O. ABC transporters: the power to change. Nat. Rev. Mol. Cell Biol. 2009, 10, 218−227. (30) Jones, P. M.; George, A. M. ABC transporters: a riddle wrapped in a mystery inside an enigma. Trends Biochem. Sci. 2009, 34, 520− 531. (31) ter Beek, J.; Guskov, A.; Slotboom, D. J. Structural diversity of ABC transporters. J. Gen. Physiol. 2014, 143, 419−435. (32) Locher, K. P. Mechanistic diversity in ATP-binding cassette (ABC) transporters. Nat. Struct. Mol. Biol. 2016, 23, 487−493. (33) Dong, J.; Yang, G.; McHaourab, H. S. Structural basis of energy transduction in the transport cycle of MsbA. Science 2005, 308, 1023− 1028. (34) Zou, P.; Mchaourab, H. S. Alternating access of the putative substrate-binding chamber in the ABC transporter MsbA. J. Mol. Biol. 2009, 393, 574−585. (35) Mishra, S.; Verhalen, B.; Stein, R. A.; Wen, P.-C.; Tajkhorshid, E.; Mchaourab, H. S. Conformational dynamics of the nucleotide binding domains and the power stroke of a heterodimeric ABC transporter. eLife 2014, 3, e02740. (36) Oldham, M. L.; Chen, J. Snapshots of the maltose transporter during ATP hydrolysis. Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 15152−15156.

(37) Zhou, Y.; Ojeda-May, P.; Pu, J. H-loop histidine catalyzes ATP hydrolysis in the E. coli ABC-transporter HlyB. Phys. Chem. Chem. Phys. 2013, 15, 15811−15815. (38) Geourjon, C.; Orelle, C.; Steinfels, E.; Blanchet, C.; Deléage, G.; Pietro, A. D.; Jault, J.-M. A common mechanism for ATP hydrolysis in ABC transporter and helicase superfamilies. Trends Biochem. Sci. 2001, 26, 539−544. (39) Moody, J. E.; Millen, L.; Binns, D.; Hunt, J. F.; Thomas, P. J. Cooperative, ATPdependent association of the nucleotide binding cassettes during the catalytic cycle of ATP-binding cassette transporters. J. Biol. Chem. 2002, 277, 21111−21114. (40) Orelle, C.; Dalmas, O.; Gros, P.; Di Pietro, A.; Jault, J. M. The conserved glutamate residue adjacent to the Walker-B motif is the catalytic base for ATP hydrolysis in the ATP-binding cassette transporter BmrA. J. Biol. Chem. 2003, 278, 47002−47008. (41) Huang, W.; Liao, J.-L. Catalytic mechanism of the maltose transporter hydrolyzing ATP. Biochemistry 2016, 55, 224−231. (42) Hsu, W.-L.; Furuta, T.; Sakurai, M. ATP hydrolysis mechanism in a maltose transporter explored by QM/MM metadynamics simulation. J. Phys. Chem. B 2016, 120, 11102−11112. (43) Locher, K. P.; Lee, A. T.; Rees, D. C. The E. coli BtuCD structure: a framework for ABC transporter architecture and mechanism. Science 2002, 296, 1091−1098. (44) Hvorup, R. N.; Goetz, B. A.; Niederer, M.; Hollenstein, K.; Perozo, E.; Locher, K. P. Asymmetry in the structure of the ABC transporter-binding protein complex BtuCDBtuF. Science 2007, 317, 1387−1390. (45) Korkhov, V. M.; Mireku, S. A.; Hvorup, R. N.; Locher, K. P. Asymmetric states of vitamin B12 transporter BtuCD are not discriminated by its cognate substrate binding protein BtuF. FEBS Lett. 2012, 586, 972−976. (46) Korkhov, V. M.; Mireku, S. A.; Locher, K. P. Structure of AMPPNP-bound vitamin B12 transporter BtuCD-F. Nature 2012, 490, 367−372. (47) Korkhov, V. M.; Mireku, S. A.; Veprintsev, D. B.; Locher, K. P. Structure of AMP-PNP−bound BtuCD and mechanism of ATPpowered vitamin B12 transport by BtuCD−F. Nat. Struct. Mol. Biol. 2014, 21, 1097−1099. (48) Åqvist, J.; Kamerlin, S. C. L. Exceptionally large entropy contributions enable the high rates of GTP hydrolysis on the ribosome. Sci. Rep. 2015, 5, 15817. (49) Zhao, Y.; Truhlar, D. G. Hybrid meta density functional theory methods for thermochemistry, thermochemical kinetics, and noncovalent interactions: the MPW1B95 and MPWB1K models and comparative assessments for hydrogen bonding and van der Waals interactions. J. Phys. Chem. A 2004, 108, 6908−6918. (50) Ribeiro, A. J. M.; Ramos, M. J.; Fernandes, P. A. Benchmarking of DFT functionals for the hydrolysis of phosphodiester bonds. J. Chem. Theory Comput. 2010, 6, 2281−2292. (51) Timachi, M. H.; Hutter, C. A.; Hohl, M.; Assafa, T.; Böhm, S.; Mittal, A.; Seeger, M. A.; Bordignon, E. Exploring conformational equilibria of a heterodimeric ABC transporter. eLife 2017, 6, e20236. (52) Göddeke, H.; Timachi, M. H.; Hutter, C. A. J.; Galazzo, L.; Seeger, M. A.; Karttunen, M.; Bordignon, E.; Schäfer, L. V. Atomistic mechanism of large-scale conformational transition in a heterodimeric ABC exporter. J. Am. Chem. Soc. 2018, 140, 4543−4551. (53) Johnson, Z. L.; Chen, J. ATP binding enables substrate release from Multidrug Resistance Protein 1. Cell 2018, 172, 81−89. (54) Prieß, M.; Schäfer, L. V. Release of entropic spring reveals conformational coupling mechanism in the ABC transporter BtuCDF. Biophys. J. 2016, 110, 2407−2418. (55) Jones, P. M.; George, A. M. Role of the D-loops in allosteric control of ATP hydrolysis in an ABC transporter. J. Phys. Chem. A 2012, 116, 3004−3013. (56) Jin, Y.; Molt, R. W.; Waltho, J. P.; Richards, N. G. J.; Blackburn, G. M. 19F NMR and DFT analysis reveal structural and electronic transition state features for RhoAcatalyzed GTP hydrolysis. Angew. Chem., Int. Ed. 2016, 55, 3318−3322. I

DOI: 10.1021/acscentsci.8b00369 ACS Cent. Sci. XXXX, XXX, XXX−XXX

Research Article

ACS Central Science

(76) Meagher, K. L.; Redman, L. T.; Carlson, H. A. Development of polyphosphate parameters for use with the AMBER force field. J. Comput. Chem. 2003, 24, 1016−1025. (77) Søndergaard, C. R.; Olsson, M. H.; Rostkowski, M.; Jensen, J. H. Improved treatment of ligands and coupling effects in empirical calculation and rationalization of pKa values. J. Chem. Theory Comput. 2011, 7, 2284−2295. (78) Olsson, M. H.; Søndergaard, C. R.; Rostkowski, M.; Jensen, J. H. PROPKA3: consistent treatment of internal and surface residues in empirical pKa predictions. J. Chem. Theory Comput. 2011, 7, 525− 537. (79) Wilks, J. C.; Slonczewski, J. L. pH of the cytoplasm and Periplasm of Escherichia coli: rapid measurement by green fluorescent protein fluorimetry. J. Bacteriol. 2007, 189, 5601−5607. (80) Hess, B.; Kutzner, C.; van der Spoel, D.; Lindahl, E. GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 2008, 4, 435− 447. (81) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; 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.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, O.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09, revision D.01; Gaussian, Inc.: Wallingford, CT, 2009. (82) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, S. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. (83) Field, M. J.; Bash, P. A.; Karplus, M. A combined quantummechanical and molecularmechanical potential for molecular dynamics simulations. J. Comput. Chem. 1990, 11, 700−733. (84) Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 014101. (85) Hess, B. P-LINCS: a parallel linear constraint solver for molecular simulation. J. Chem. Theory Comput. 2008, 4, 116−122. (86) Miyamoto, S.; Kollman, P. A. SETTLE: an analytical version of the SHAKE and RATTLE algorithm for rigid water models. J. Comput. Chem. 1992, 13, 952−962. (87) Fixman, M. Classical statistical mechanics of constraints: a theorem and application to polymers. Proc. Natl. Acad. Sci. U. S. A. 1974, 71, 3050−3053. (88) Schlitter, J.; Klähn, M. A new concise expression for the free energy of a reaction coordinate. J. Chem. Phys. 2003, 118, 2057−2060. (89) Hess, B. Determining the shear viscosity of model liquids from molecular dynamics simulations. J. Chem. Phys. 2002, 116, 209−217.

(57) Tripathi, R.; Glaves, R.; Marx, D. The GTPase hGBP1 converts GTP to GMP in two steps via proton shuttle mechanisms. Chem. Sci. 2017, 8, 371−380. (58) Vithani, N.; Batra, S.; Prakash, B.; Nair, N. N. Elucidating the GTP hydrolysis mechanism in FeoB: a hydrophobic amino-acid substituted GTPase. ACS Catal. 2017, 7, 902−906. (59) Payen, L. F.; Gao, M.; Westlake, C. J.; Cole, S. P. C.; Deeley, R. G. Role of carboxylate residues adjacent to the conserved core Walker B motifs in the catalytic cycle of Multidrug Resistance Protein 1 (ABCC1). J. Biol. Chem. 2003, 278, 38537−38547. (60) Tombline, G.; Bartholomew, L. A.; Tyndall, G. A.; Gimi, K.; Urbatsch, I. L.; Senior, A. E. Properties of P-glycoprotein with mutations in the catalytic carboxylate glutamate residues. J. Biol. Chem. 2004, 279, 46518−46526. (61) Procko, E.; Ferrin-O’Connell, I.; Ng, S.-L.; Gaudet, R. Distinct structural and functional properties of the ATPase sites in an asymmetric ABC transporter. Mol. Cell 2006, 24, 51−62. (62) Borths, E. L.; Poolman, B.; Hvorup, R. N.; Locher, K. P.; Rees, D. C. In vitro functional characterization of BtuCD-F, the Escherichia coli ABC transporter for vitamin B12 uptake. Biochemistry 2005, 44, 16301−16309. (63) Zaitseva, J.; Jenewein, S.; Jumpertz, T.; Holland, I. B.; Schmitt, L. H662 is the linchpin of ATP hydrolysis in the nucleotide-binding domain of the ABC transporter HlyB. EMBO J. 2005, 24, 1901−1910. (64) Syberg, F.; Suveyzdis, Y.; Kötting, C.; Gerwert, K.; Hofmann, E. Time-resolved Fourier transform infrared spectroscopy of the nucleotide-binding domain from the ATP-binding cassette transporter MsbA: ATP hydrolysis is the rate-limiting step in the catalytic cycle. J. Biol. Chem. 2012, 287, 23923−23931. (65) Janas, E.; Hofacker, M.; Chen, M.; Gompf, S.; van der Does, C.; Tampé, R. The ATP hydrolysis cycle of the nucleotide-binding domain of the mitochondrial ATP-binding cassette transporter Mdl1p. J. Biol. Chem. 2003, 278, 26862−26869. (66) Ernst, R.; Koch, J.; Horn, C.; Tampé, R.; Schmitt, L. Engineering ATPase activity in the isolated ABC cassette of human TAP1. J. Biol. Chem. 2006, 281, 27471−27480. (67) Ernst, R.; Kueppers, P.; Klein, C. M.; Schwarzmueller, T.; Kuchler, K.; Schmitt, L. A mutation of the H-loop selectively affects rhodamine transport by the yeast multidrug ABC transporter Pdr5. Proc. Natl. Acad. Sci. U. S. A. 2008, 105, 5069−5074. (68) Leipe, D. D.; Wolf, Y. I.; Koonin, E. V.; Aravind, L. Classification and evolution of P-loop GTPases and related ATPases. J. Mol. Biol. 2002, 317, 41−72. (69) Vorobiev, S.; Strokopytov, B.; Drubin, D. G.; Frieden, C.; Ono, S.; Condeelis, J.; Rubenstein, P. A.; Almo, S. C. The structure of nonvertebrate actin: implications for the ATP hydrolytic mechanism. Proc. Natl. Acad. Sci. U. S. A. 2003, 100, 5760−5765. (70) Kiani, F. A.; Fischer, S. Comparing the catalytic strategy of ATP hydrolysis in biomolecular motors. Phys. Chem. Chem. Phys. 2016, 18, 20219−20233. (71) Hornak, V.; Abel, R.; Okur, A.; Strockbine, B.; Roitberg, A.; Simmerling, C. Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins: Struct., Funct., Genet. 2006, 65, 712−725. (72) Lindorff-Larsen, K.; Piana, S.; Palmo, K.; Maragakis, P.; Klepeis, J. L.; Dror, R. O.; Shaw, D. E. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins: Struct., Funct., Genet. 2010, 78, 1950−1958. (73) Berger, O.; Edholm, O.; Jähnig, F. Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature. Biophys. J. 1997, 72, 2002−2013. (74) Cordomí, A.; Caltabiano, G.; Pardo, L. Membrane protein simulations using AMBER force field and Berger lipid parameters. J. Chem. Theory Comput. 2012, 8, 948−958. (75) Horn, H. W.; Swope, W. C.; Pitera, J. W.; Madura, J. D.; Dick, T. J.; Hura, G. L.; Head-Gordon, T. Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. J. Chem. Phys. 2004, 120, 9665−9678. J

DOI: 10.1021/acscentsci.8b00369 ACS Cent. Sci. XXXX, XXX, XXX−XXX