Elementary steps in aqueous proton transfer reactions : a first principles molecular dynamics study

La nature des acides dans un environnement aqueux est primordiale dans de nombreux aspects de la chimie et de la biologie. La caractéristique principale d'un acide est sa capacité à transférer un proton vers une molécule d'eau ou vers n'importe quelle base, mais ce procédé n'est pas aussi simple qu'il y paraît. Il peut au contraire être extrêmement complexe et dépendre de manière cruciale de la solvatation des différents intermédiaires de réaction impliqués. Cette thèse décrit les études computationnelles basées sur des simulations de dynamique moléculaire ab initio qui ont pour but d'obtenir une description à l'échelle moléculaire des divers procédés de transferts de proton entre acide et bases dans un milieu aqueux. Pour cela, nous avons étudié une serie de système, dont l'acide hydrofluorique aqueux, l'acide trifluoroacétique aqueux, et un système modèle constitué d'un phénol et d'une entité carboxylate reliés entre eux par une molécule d'eau en solution aqueuse. Deux états intermédiaires ont été identifiés pour le transfert d'un proton depuis un acide. Ces intermédiaires apparaissent stabilisés par un motif local de solvatation via des ponts H. Leurs signatures spectroscopiques ont été caractérisées au moyen de la spectroscopie infrarouge, en utilisant le formalisme de la dynamique moléculaire ab initio, qui inclut l'effet quantique nucléaire de manière explicite. Cette étude a aussi identifié trois chemins de réaction élémentaire, qui sont responsable pour le transfert d'un proton d'un acide à une base, ainsi que leurs échelles de temps caractéristiques. Les conclusions tirées de ces études sont discutées dans les détails, au niveau moléculaire, avec une emphase sur les comparaisons entre les résultats théoriques et les mesures expérimentales obtenues dans a littérature ou via des collaborateurs.

Temperature- and Pressure-induced Proton Transfer in the 1:1 Adduct Formed between Squaric Acid and 4,4 '-Bipyridine

We have applied a combination of spectroscopic and diffraction methods to study the adduct formed between squaric acid and bypridine, which has been postulated to exhibit proton transfer associated with a single-crystal to single-crystal phase transition at ca. 450 K. A combination of X-ray single-crystal and very-high flux powder neutron diffraction data confirmed that a proton does transfer from the acid to the base in the high-temperature form. Powder X-ray diffraction measurements demonstrated that the transition was reversible but that a significant kinetic energy barrier must be overcome to revert to the original structure. Computational modeling is consistent with these results. Modeling also revealed that, while the proton transfer event would be strongly discouraged in the gas phase, it occurs in the solid state due to the increase in charge state of the molecular ions and their arrangement inside the lattice. The color change is attributed to a narrowing of the squaric acid to bipyridine charge-transfer energy gap. Finally, evidence for the possible existence of two further phases at high pressure is also presented.

Intramolecular proton transfer in [Ph2B-OH2](+). Participation of a phenyl group

Ab initio calculations at the HF/6-31+G* level on [Ph2B-OH2](+) show that in the gas phase the structure with the proton attached to an ipso C is lower in energy than the one with the proton on the oxygen atom by 8.40 kcal mol(-1). The transition states and reaction paths for intramolecular proton transfer in [Ph2B-OH2](+) have also been studied.

Solvent Effects in Chemical Processes. Water-Assisted Proton Transfer Reaction of Pterin in Aqueous Environment

Pterins are members of a family of heterocyclic compounds present in a wide variety of biological systems and may exist in two forms, corresponding to an acid and a basic tautomer. In this work, the proton transfer reaction between these tautomeric forms was investigated in the gas phase and in aqueous solution. In gas phase, the intramolecular mechanism was carried out for die isolated pterin by quantum mechanical second-order Moller-Plesset Perturbation theory (MP2/aug-cc-pVDZ) calculations and it indicates that the acid form is more stable than the basic form by -1.4 kcal/mol with a barrier of 34.2 kcal/mol with respect to the basic form. In aqueous solution, the role of the water molecules in the proton transfer reaction was analyzed in two separated parts, the direct participation of one water molecule in the reaction path, called water-assisted mechanism, and the complementary participation of the aqueous solvation. The water-assisted mechanism was carried out for one pterin-water cluster by quantum mechanical calculations and it indicates that the acid form is still more stable by -3.3 kcal/mol with a drastic reduction of 70% of the barrier, The bulk solution effect on the intramolecular and water-assisted mechanisms was included by free energy perturbation implemented on Monte Carlo simulations. The bulk water effect is found to be substantial and decisive when the reaction path involves the water-assisted mechanism. In this case, the free energy barrier is only 6.7 kcal/mol and the calculated relative Gibbs free energy for the two tautomers is -11.2 kcal/mol. This value is used to calculate the pK(a) value of 8.2 +/- 0.6 that is in excellent agreement with the experimental result of 7.9.

Picosecond Dynamics of Proton Transfer of a 7-Hydroxyflavylium Salt in Aqueous-Organic Solvent Mixtures

The intermediacy of the geminate base proton pair (A*center dot center dot center dot H(+)) in excited-state proton-transfer (ESPT) reactions (two-step mechanism) has been investigated employing the synthetic flavylium salt 7-hydroxy-4-methyl-flavylium chloride (HMF). In aqueous solution, the ESPT mechanism involves solely the excited acid AH* and base A* forms of HMF as indicated by the fluorescence spectra and double-exponential fluorescence decays (two species, two decay times). However, upon addition of either 1,4-dioxane or 1,2-propylene glycol, the decays become triple-exponential with a term consistent with the presence of the geminate base proton pair A*center dot center dot center dot H(+). The geminate pair becomes detectable because of the increase in the recombination rate constant, k(rec), of (A*center dot center dot center dot H(+)) with increasing the mole fraction of added organic cosolvent. Because the two-step ESPT mechanism splits the intrinsic prototropic reaction rates (deprotonation of AH(+)*, k(d), and recombination, k(rec) of A*center dot center dot center dot H(+)) from the diffusion controlled rates (dissociation, k(diss) and formation, k(diff)[H(+)], of A*center dot center dot center dot H+), the experimental detection of the geminate pair provides a wealth of information on the proton-transfer reaction (k(d) and k(rec)) as well as on proton diffusion/migration (k(diss) and k(diff)).

A theoretical analysis on the intramolecular proton transfer of alpha-alanine in an aqueous medium

Intramolecular proton transfer from oxygen to nitrogen atoms in the alpha-alanine amino acid has been studied by ab initio methods at the HF/6-31G*, HF/6-31 ++ G** and MP2/6-31 ++ G** levels of calculation including the solvent effects by means of self-consistent reaction field theory. An analysis of the results based on the natural bond orbital charges shows that the transition structure presents an imbalance in the sense that the charge shift lags behind the proton transfer and that the bond formation is always in advance with respect to the bond cleavage. All calculation levels show that the barrier height associated with the conformational change on alpha-alanine is larger than the proton transfer process. (C) 1998 Elsevier B.V. B.V. All rights reserved.

Beyond simple proton transfer processes: domino reactions between 4-substituted indoles and nitroethene as a new gateway to ergot alkaloids

The multimodal biology activity of ergot alkaloids is known by humankind since middle ages. Synthetically modified ergot alkaloids are used for the treatment of various medical conditions. Despite the great progress in organic syntheses, the total synthesis of ergot alkaloids remains a great challenge due to the complexity of their polycyclic structure with multiple stereogenic centres. This project has developed a new domino reaction between indoles bearing a Michael acceptor at the 4 position and nitroethene, leading to potential ergot alkaloid precursors in highly enantioenriched form. The reaction was optimised and applied to a large variety of substrate with good results. Even if unfortunately all attempts to further modify the obtained polycyclic structure failed, it was found a reaction able to produce the diastereoisomer of the polycyclic product in excellent yields. The compounds synthetized were characterized by NMR and ESIMS analysis confirming the structure and their enantiomeric excess was determined by chiral stationary phase HPLC. The mechanism of the reaction was evaluated by DFT calculations, showing the formation of a key bicoordinated nitronate intermediate, and fully accounting for the results observed with all substrates. The relative and absolute configuration of the adducts were determined by a combination of NMR, ECD and computational methods.

A coupled proton-transfer and twisting-motion fluorescence probe for lipid bilayers

A new and sensitive molecular probe, 2-(2′-hydroxyphenyl)imidazo[1,2-a]pyridine (HPIP), for monitoring structural changes in lipid bilayers is presented. Migration of HPIP from water into vesicles involves rupture of hydrogen (H) bonds with water and formation of an internal H bond once the probe is inside the vesicle. These structural changes of the dye allow the occurrence of a photoinduced intramolecular proton-transfer reaction and a subsequent twisting/rotational process upon electronic excitation of the probe. The resulting large Stokes-shifted fluorescence band depends on the twisting motion of the zwitterionic phototautomer and is characterized in vesicles of dimyristoyl-phosphatidylcholine and in dipalmitoyl-phosphatidylcholine at the temperature range of interest and in the presence of cholesterol. Because the fluorescence of aqueous HPIP does not interfere in the emission of the probe within the vesicles, HPIP proton-transfer/twisting motion fluorescence directly allows us to monitor and quantify structural changes within bilayers. The static and dynamic fluorescence parameters are sensitive enough to such changes to suggest this photostable dye as a potential molecular probe of the physical properties of lipid bilayers.

The roles of the two proton input channels in cytochrome c oxidase from Rhodobacter sphaeroides probed by the effects of site-directed mutations on time-resolved electrogenic intraprotein proton transfer

The crystal structures of cytochrome c oxidase from both bovine and Paracoccus denitrificans reveal two putative proton input channels that connect the heme-copper center, where dioxygen is reduced, to the internal aqueous phase. In this work we have examined the role of these two channels, looking at the effects of site-directed mutations of residues observed in each of the channels of the cytochrome c oxidase from Rhodobacter sphaeroides. A photoelectric technique was used to monitor the time-resolved electrogenic proton transfer steps associated with the photo-induced reduction of the ferryl-oxo form of heme a3 (Fe4+ = O2−) to the oxidized form (Fe3+OH−). This redox step requires the delivery of a “chemical” H+ to protonate the reduced oxygen atom and is also coupled to proton pumping. It is found that mutations in the K channel (K362M and T359A) have virtually no effect on the ferryl-oxo-to-oxidized (F-to-Ox) transition, although steady-state turnover is severely limited. In contrast, electrogenic proton transfer at this step is strongly suppressed by mutations in the D channel. The results strongly suggest that the functional roles of the two channels are not the separate delivery of chemical or pumped protons, as proposed recently [Iwata, S., Ostermeier, C., Ludwig, B. & Michel, H. (1995) Nature (London) 376, 660–669]. The D channel is likely to be involved in the uptake of both “chemical” and “pumped” protons in the F-to-Ox transition, whereas the K channel is probably idle at this partial reaction and is likely to be used for loading the enzyme with protons at some earlier steps of the catalytic cycle. This conclusion agrees with different redox states of heme a3 in the K362M and E286Q mutants under aerobic steady-state turnover conditions.

Femtosecond cluster studies of the solvated 7-azaindole excited state double-proton transfer

Presented here are femtosecond pump-probe studies on the water-solvated 7-azaindole dimer, a model DNA base pair. In particular, studies are presented that further elucidate the nature of the reactive and nonreactive dimers and also provide new insights establishing that the excited state double-proton transfer in the dimer occurs in a stepwise rather than a concerted manner. A major question addressed is whether the incorporation of a water molecule with the dimer results in the formation of species that are unable to undergo excited state double-proton transfer, as suggested by a recent study reported in the literature [Nakajima, A., Hirano, M., Hasumi, R., Kaya, K., Watanabe, H., Carter, C. C., Williamson, J. M. & Miller, T. (1997) J. Phys. Chem. 101, 392–398]. In contrast to this earlier work, our present findings reveal that both reactive and nonreactive dimers can coexist in the molecular beam under the same experimental conditions and definitively show that the clustering of water does not induce the formation of the nonreactive dimer. Rather, when present with a species already determined to be a nonreactive dimer, the addition of water can actually facilitate the occurrence of the proton transfer reaction. Furthermore, on attaining a critical hydration number, the data for the nonreactive dimer suggest a solvation-induced conformational structure change leading to proton transfer on the photoexcited half of the 7-azaindole dimer.

Proton transfer from the bulk to the bound ubiquinone QB of the reaction center in chromatophores of Rhodobacter sphaeroides: Retarded conveyance by neutral water

The mechanism of proton transfer from the bulk into the membrane protein interior was studied. The light-induced reduction of a bound ubiquinone molecule QB by the photosynthetic reaction center is accompanied by proton trapping. We used kinetic spectroscopy to measure (i) the electron transfer to QB (at 450 nm), (ii) the electrogenic proton delivery from the surface to the QB site (by electrochromic carotenoid response at 524 nm), and (iii) the disappearance of protons from the bulk solution (by pH indicators). The electron transfer to QB− and the proton-related electrogenesis proceeded with the same time constant of ≈100 μs (at pH 6.2), whereas the alkalinization in the bulk was distinctly delayed (τ ≈ 400 μs). We investigated the latter reaction as a function of the pH indicator concentration, the added pH buffers, and the temperature. The results led us to the following conclusions: (i) proton transfer from the surface-located acidic groups into the QB site followed the reduction of QB without measurable delay; (ii) the reprotonation of these surface groups by pH indicators and hydronium ions was impeded, supposedly, because of their slow diffusion in the surface water layer; and (iii) as a result, the protons were slowly donated by neutral water to refill the proton vacancies at the surface. It is conceivable that the same mechanism accounts for the delayed relaxation of the surface pH changes into the bulk observed previously with bacteriorhodopsin membranes and thylakoids. Concerning the coupling between proton pumps in bioenergetic membranes, our results imply a tendency for the transient confinement of protons at the membrane surface.

Energy transfer to a proton-transfer fluorescence probe: Tryptophan to a flavonol in human serum albumin

A protein fluorescence probe system, coupling excited-state intermolecular Förster energy transfer and intramolecular proton transfer (PT), is presented. As an energy donor for this system, we used tryptophan, which transfers its excitation energy to 3-hydroxyflavone (3-HF) as a flavonol prototype, an acceptor exhibiting excited-state intramolecular PT. We demonstrate such a coupling in human serum albumin–3-HF complexes, excited via the single intrinsic tryptophan (Trp-214). Besides the PT tautomer fluorescence (λmax = 526 nm), these protein–probe complexes exhibit a 3-HF anion emission (λmax = 500 nm). Analysis of spectroscopic data leads to the conclusion that two binding sites are involved in the human serum albumin–3-HF interaction. The 3-HF molecule bound in the higher affinity binding site, located in the IIIA subdomain, has the association constant (k1) of 7.2 × 105 M−1 and predominantly exists as an anion. The lower affinity site (k2 = 2.5 × 105 M−1), situated in the IIA subdomain, is occupied by the neutral form of 3-HF (normal tautomer). Since Trp-214 is situated in the immediate vicinity of the 3-HF normal tautomer bound in the IIA subdomain, the intermolecular energy transfer for this donor/acceptor pair has a 100% efficiency and is followed by the PT tautomer fluorescence. Intermolecular energy transfer from the Trp-214 to the 3-HF anion bound in the IIIA subdomain is less efficient and has the rate of 1.61 × 108 s−1, thus giving for the donor/acceptor distance a value of 25.5 Å.

Large kinetic isotope effects in enzymatic proton transfer and the role of substrate oscillations

We propose an interpretation of the experimental findings of Klinman and coworkers [Cha, Y., Murray, C. J. & Klinman, J. P. (1989) Science 243, 1325–1330; Grant, K. L. & Klinman, J. P. (1989) Biochemistry 28, 6597–6605; and Bahnson, B. J. & Klinman, J. P. (1995) Methods Enzymol. 249, 373–397], who showed that proton transfer reactions that are catalyzed by bovine serum amine oxidase proceed through tunneling. We show that two different tunneling models are consistent with the experiments. In the first model, the proton tunnels from the ground state. The temperature dependence of the kinetic isotope effect is caused by a thermally excited substrate mode that modulates the barrier, as has been suggested by Borgis and Hynes [Borgis, D. & Hynes, J. T. (1991) J. Chem. Phys. 94, 3619–3628]. In the second model, there is both over-the-barrier transfer and tunneling from excited states. Finally, we propose two experiments that can distinguish between the possible mechanisms.

Determination of the binding sites of the proton transfer inhibitors Cd2+ and Zn2+ in bacterial reaction centers

The reaction center (RC) from Rhodobacter sphaeroides couples light-driven electron transfer to protonation of a bound quinone acceptor molecule, QB, within the RC. The binding of Cd2+ or Zn2+ has been previously shown to inhibit the rate of reduction and protonation of QB. We report here on the metal binding site, determined by x-ray diffraction at 2.5-Å resolution, obtained from RC crystals that were soaked in the presence of the metal. The structures were refined to R factors of 23% and 24% for the Cd2+ and Zn2+ complexes, respectively. Both metals bind to the same location, coordinating to Asp-H124, His-H126, and His-H128. The rate of electron transfer from QA− to QB was measured in the Cd2+-soaked crystal and found to be the same as in solution in the presence of Cd2+. In addition to the changes in the kinetics, a structural effect of Cd2+ on Glu-H173 was observed. This residue was well resolved in the x-ray structure—i.e., ordered—with Cd2+ bound to the RC, in contrast to its disordered state in the absence of Cd2+, which suggests that the mobility of Glu-H173 plays an important role in the rate of reduction of QB. The position of the Cd2+ and Zn2+ localizes the proton entry into the RC near Asp-H124, His-H126, and His-H128. Based on the location of the metal, likely pathways of proton transfer from the aqueous surface to QB⨪ are proposed.