Proton Transfer DOI: 10.1002/anie.200603583 Dynamic Protonation Equilibrium of Solvated Acetic Acid** WeiGu,TomasoFrigato,TjerkP.Straatsma,andVolkhardHelms* The biological functions of many proteins are crucially coupled to protonation equilibria, for instance, in the enzymatic reactions of serine proteases [1–3] and carbonic anhydrase, [4] and in integral membrane proton pumps such as bacteriorhodopsin, [5] cytochrome c oxidase (COX), [6,7] and F O F 1 -ATP synthase. [8] Proton-transfer (PT) reactions are also crucial in many other areas of chemistry, such as membrane permeation in hydrogen fuel cells and in polymers. [9] Despite their tremendous importance, many aspects of PTreactions in biomolecules remain poorly understood. As the direct observation of PT reactions by experimental techniques faces fundamental and/or technical difficulties, it is highly desirable to use computational methods as a complement. In the past decades, various efficient computational methods have been developed to calculate the pK a values of amino acid side chains as well as to perform simulations of proteins at constant pH, [10–17] using, for example, fractional charges [10] and implicit solvent models. [16,17] However, these methods do not model explicit proton-exchange reactions between the titratable sites and the surrounding aqueous solution or the exchange between different titratable sites; this makes it difficult to identify PT pathways and to characterize the mechanisms of PT reactions. This is the area where dynamic simulations of explicit proton-transfer reactions come into play. Tuckerman etal. studied the shared proton in hydrogen bonds [18] and a hydrated excess proton in water [19] using the Car–Parrinello molecular dynamics (CPMD) method. [20,21] Lobaugh and Voth investigated proton transport in water by simulating an excess proton in a box of water molecules [22] within the centroid molecular dynamics [23] framework and by using a multistate empirical valence bond (MS-EVB) model for proton transfer. [24–26] A recent study also presented the dynamic simulation of pK a values for amino acid side-chain analogues [27] using the MS-EVB model and the umbrella sampling technique. [28,29] The deviation between their com- puted value and the experimental pK a was 1–2 pK a units. CPMD combined with metadynamics and transition-path sampling was employed to compute free-energy profiles for the deprotonation of acetic acid in water. [30] Several further applications showed the importance and success of studying PT in protein systems by theoretical approaches, for instance, PT in bacteriorhodopsin, [31] PT in gramicidin A, [32,33] PT along a water chain in the D-pathway of COX, [34] and proton translocation in carbonic anhydrase. [35] A simulation model of intermediate accuracy, the Q-HOP MD method, was introduced by our research group in 2001 to study dynamic proton transport between general titratable sites in biomolecular systems. [36–39] In the Q-HOP scheme, the PT probabilities for each proton-donor and -acceptor pair are calculated using a semiempirical approach during the MD simulation (see the Supporting Information for details). Depending on whether PT occurs (by comparing the PT probability to a random number), the topology of the system is modified or kept unchanged before the next step of MD simulation. The transfer probabilities depend on the momen- tary donor–acceptor distance (R DA ) and the energy difference between the minima at the donor and acceptor (E 12 ). This method has been applied successfully to study the proton shuttle in green-fluorescent protein [40] and to understand the mechanism of proton blockage in aquaporin. [41] Herein we present the application of Q-HOP-MD to study the explicit protonation equilibrium of solvated acetic acid on a time scale of tens of nanoseconds &German: 50 nm & and at a reasonable pH (pH 1). The pK a of acetic acid is calculated based on the relative populations of protonated and deprotonated states observed during a 50-ns-long Q- HOP MD simulation. By analyzing the unbiased MD simulation, we can also identify the proton-hopping mecha- nism and the driving force of the activated processes of proton transfer. This study thus serves as a proof of principle of the method, and it provides detailed mechanistic insight into atomistic protonation equilibria on two separate time scales, femtoseconds and nanoseconds. During the Q-HOP MD simulation, two types of proto- nation equilibria were observed. Figure 1a,b show the position of the “free” proton and the distance between the hydronium ion and the deprotonated acetic acid (when the proton stays on hydronium ion) during the first 10 ns. (The results for the full length of the simulation are presented in Figure 5 of the Supporting Information.) Two different situations can be distinguished. The first type of protonation equilibrium, “proton swapping”, involves only the acetic acid and a nearby H 2 O/H 3 O + molecule, which forms a hydrogen Dateiname: Z603583e Pagina: 1 Pfad: l:/daten/verlage/vch/ach/hefte/pool/ Seite: 1 te von 6 Status Neusatz Umfang (Seiten): 6 Datum: 8 KW., 19. Februar 2007 (Montag) Zeit: 14:01:52 Uhr [*] W. Gu,Dr. T. Frigato, Prof.Dr. V. Helms Zentrum für Bioinformatik Universität des Saarlandes 66041 Saarbrücken (Germany) Fax: (+ 49) 681-302-64180 E-mail: volkhard.helms@bioinformatik.uni-saarland.de Homepage: http://gepard.bioinformatik.uni-saarland.de Dr. T. P. Straatsma Computational Sciences and Mathematics Division Pacific Northwest National Laboratory Richland, WA 99352 (USA) [**] NWchem Version 4.7, as developed and distributed by the Pacific Northwest National Laboratory (PNNL). P.O. Box 999, Richland, WA 99352 (USA), and funded by the US Department of Energy, was used for the calculation. We thank Elena Herzog from MPI Frankfurt for computation of the cis conformation of acetic acid. We thank Dr. Michael Hutter and Dr. TihamØr Geyer for valuable discussion. We thank the EMSL Grand Challenge Project (project gc3551) from PNNL for providing the computing resources. Supporting information for this article is available on the WWW under http://www.angewandte.org or from the author. Angewandte Chemie 1 Angew. Chem. Int. Ed. 2007, 46,1–6 # 2007 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim These are not the final page numbers! Ü Ü 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 11 11 12 12 13 13 14 14 15 15 16 16 17 17 18 18 19 19 20 20 21 21 22 22 23 23 24 24 25 25 26 26 27 27 28 28 29 29 30 30 31 31 32 32 33 33 34 34 35 35 36 36 37 37 38 38 39 39 40 40 41 41 42 42 43 43 44 44 45 45 46 46 47 47 48 48 49 49 50 50 51 51 52 52 53 53 54 54 55 55 56 56 57 57 58 58 59 59