Investigation of the intermolecular proton transfer in the supersystems adenine-methanol/ethanol/i-propanol: MP2 and DFT levels study

Twelve H-bonded supersystems constructed between the adenine tautomers and methanol, ethanol, and i-propanol were studied at the B3LYP and MP2 levels of theory using 6-311G(d,p) and 6-311++G(d,p) basis functions. The thermodynamic parameters of the complex formations were calculated in order to estimate the exact stability of the supersystems. It was proven that the calculated energy barriers of the alcohol-assisted proton transfers are about 60% lower than those of the intramolecular proton transfers in adenine found earlier (Gu and Leszczynski in J Phys Chem A 103:2744–2750, 1999). Figure H-bonded complex between i-propanol and adenine     

Theoretical study of the intermolecular H-bonding and intermolecular proton transfer between isocytosine tautomeric forms and <Emphasis Type="Italic">R</Emphasis>,<Emphasis Type="Italic">S</Emphasis>-lactic acid

Eight H-bonded complexes between isocytosine (isoC) tautomeric forms and R/S-lactic acid (LA) have been studied at the B3LYP and HF levels of theory using 6–31+G(d) basis set. The energy barriers of the intermolecular proton transfers were also estimated as the results showed that they are several times lower than those of the intramolecular proton transfers of isoC in the gas phase. Furthermore, the energy barriers of the tautomerizations in which the carboxylic H-atom takes part are several times lower than those in which the LA OH group assists the proton transfer. Figure     

H-bonded complexes between acetylacetone and two molecules of methanol: HF and DFT level study

Five stable H-bonded complexes (supersystems) between acetylacetone and two methanol molecules were investigated at the B3LYP and HF levels of theory using the 6-311G** and 6-11++G** basis sets. The most stable complex was found as the one with the highest relative bonding and interaction energies. All vibrational frequencies resulting from calculations with the 6-311++G** basis set were compared with the recorded IR spectrum of acetylacetone/methanol mixture in a molar ratio 1:2.     

Different types of biological proton transfer reactions studied by quantum chemical methods

Different types of proton transfer occurring in biological systems are described with examples mainly from ribonucleotide reductase (RNR) and cytochrome c oxidase (CcO). Focus is put on situations where electron and proton transfer are rather strongly coupled. In the long range radical transfer in RNR, it is shown that the presence of hydrogen atom transfer (HAT) is the most logical explanation for the experimental observations. In another example from RNR, it is shown that a transition state for concerted motion of both proton and electron can be found even if the donors are separated by a quite long distance. In CcO, the essential proton transfer for the OO bond cleavage, and the most recent modelings of proton translocation are described, indicating a few remaining major problems.     

DFT: a dynamic study of the interaction of ethanol and methanol with platinum

The interaction between alcohol molecules and platinum (Pt) was studied using molecular dynamics (MD; Born-Oppenheimer method). Alcohol molecules like ethanol and methanol present a similar molecular structure, with a methyl group (CH₃) at one end and a fragment of hydroxyl (OH) at the other. This fact generates two orientations that are considered in the interaction with Pt. The MD calculation results for these two orientations indicate a preferential orientation due to energy interactions. A plausible reaction mechanism that takes into account the interaction between Pt and alcohol is presented. The charge transference obtained from the Pt–alcohol interaction was also analyzed. The energy for the two orientations was calculated by indicating the preferential orientation. The methyl and hydroxyl groups are involved in heterolytic breakage of hydrogen bonds, joined to a carbon atom in the former and to an oxygen atom in the latter; however, the methyl group reaction seems to be the most important.     

The interactions of square platinum(II) complexes with guanine and adenine: a quantum-chemical ab initio study of metalated tautomeric forms

The influence of binding of square planar platinum complexes on tautomeric equilibria of the DNA bases guanine and adenine was investigated using the density functional B3LYP method. Neutral trans-dichloro(amine)-, +1 charged chloro(diamine)-, and +2 charged triamine-platinum(II) species were chosen for coordination to bases. Only the N₇ interaction site of the bases was considered. The calculations demonstrate that the neutral platinum adduct does not change the tautomeric equilibria of the bases. Furthermore, N₇ binding of the neutral Pt adduct moderately reduces the probability of protonation of the N₁ position of adenine. Larger effects can be observed for +1 and mainly +2 adducts, but these can be rationalized by electrostatic effects. Since the electrostatic effects are expected to be efficiently compensated for by a charged backbone of DNA and counterions in a polar solvent, no dramatic increase in mispair formation is predicted for Pt(II) adducts, which is in agreement with experiment. The interaction energies between Pt adducts and the nucleobases were also evaluated. These interaction energies range from ca. 210 kJ/mol for neutral adducts, interacting with both bases and their tautomers, up to 500 kJ/mol for the +2 charged adducts, interacting with the keto-guanine tautomer and the anti-imino-adenine tautomer. The surprisingly large interaction energy for the latter structure is due to the strong H-bond between the NH₃ ligand group of the metal adduct and the N₆ nitrogen atom of the base. Received: 6 July 1999 / Accepted: 7 December 1999     

The kinetics of a weakly electron-coupled proton transfer in azurin

The thermodynamics and kinetics of the weakly electron-coupled proton transfer of Pseudomonas aeruginosa azurin have been studied and quantitatively determined using fast scan cyclic voltammetry (CV), in which the protein is adsorbed on two types of electrode (pyrolytic graphite ‘edge’ (PGE) or 1-decanethiol modified gold). Electron transfer is coupled to a slow protonation of His35, which is not a ligand to the copper ion but is located approximately 8 Å away. Protonation of His35 produces small changes in the reduction potential of the copper site, which are time-resolved within the scan rate range 0.01-100 V s⁻¹.     

DFT and MP2 study of low barrier proton transfer in hydrazide schiff base tautomers via water bridges and in the gas

MP2 and DFT studies were performed on the tautomers of N′-ethylideneacetohydrazide in different environments including gas phase, continuum solvent and microhydrated environment. The ground electronic state structures of the tautomers were optimized at the MP2 and B3LYP levels of theory using 6-311++G(d,p), separately. The optimized geometries of the transition states of different tautomerism processes, which occur through the proton transfer (PT) reaction, were determined using the QST3 approach at the same levels of theory. The same stability order as T1Z〉 T1E〉 T2ZE〉 T2ZZ〉 T2EE〉 T2EZ〉 T6〉 T4E was found for the tautomers in the gas phase, continuum solvent and microhydrated environment for both B3LYP and MP2 levels of theory. In addition, the variations of the Gibbs free energies of tautomeric processes, the activation Gibbs free energies of the forward and reverse tautomeric processes with the polarity of the solvent (in continuum solvent model) and the number of water molecules (in microhydrated environment) were investigated. It was found that the reverse tautomeric process is more favorable in all considered environments thermodynamically and kinetically. In addition, it was shown that the rate constants of the reverse and forward considered tautomeric processes decrease with the solvent polarity in the continuum solvent model and the process becomes more difficult than the gas phase. The opposite trend is seen in the microhydrated environment.     

Formation of hydrogen peroxide by VUV-photolysis of water and aqueous solutions with methanol

The hydrogen peroxide production upon vacuum ultraviolet (VUV) irradiation of water is reviewed, because published results from the last 10 years lead to conflicting mechanistic interpretations. This work confirms that in pure water, hydrogen peroxide is only produced in the presence of molecular oxygen. Mechanistic schemes explain these findings and confirm earlier statements that recombination of hydroxyl radicals is kinetically disfavoured. In agreement with other recent publications, this work confirms that enhanced hydrogen peroxide production takes place upon VUV irradiation of aqueous solutions of organic compounds. For these investigations, methanol was chosen as an organic model compound. During photolyses, hydrogen peroxide, dissolved molecular oxygen, pH-value of the reaction system, methanol and its products of oxidative degradation were analyzed, and kinetic studies were undertaken to explain the evolution of the concentrations of these components.     

Proton pumping in cytochrome c oxidase: Energetic requirements and the role of two proton channels

Cytochrome c oxidase is a superfamily of membrane bound enzymes catalyzing the exergonic reduction of molecular oxygen to water, producing an electrochemical gradient across the membrane. The gradient is formed both by the electrogenic chemistry, taking electrons and protons from opposite sides of the membrane, and by proton pumping across the entire membrane. In the most efficient subfamily, the A-family of oxidases, one proton is pumped in each reduction step, which is surprising considering the fact that two of the reduction steps most likely are only weakly exergonic. Based on a combination of quantum chemical calculations and experimental information, it is here shown that from both a thermodynamic and a kinetic point of view, it should be possible to pump one proton per electron also with such an uneven distribution of the free energy release over the reduction steps, at least up to half the maximum gradient. A previously suggested pumping mechanism is developed further to suggest a reason for the use of two proton transfer channels in the A-family. Since the rate of proton transfer to the binuclear center through the D-channel is redox dependent, it might become too slow for the steps with low exergonicity. Therefore, a second channel, the K-channel, where the rate is redox-independent is needed. A redox-dependent leakage possibility is also suggested, which might be important for efficient energy conservation at a high gradient. A mechanism for the variation in proton pumping stoichiometry over the different subfamilies of cytochrome oxidase is also suggested. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.     

Mechanism of proton transfer through the KC proton pathway in the Vibrio cholerae cbb3 terminal oxidase

The heme?copper oxidases (HCuOs) are terminal components of the respiratory chain, catalyzing oxygen reduction coupled to the generation of a proton motive force. The C-family HCuOs, found in many pathogenic bacteria under low oxygen tension, utilize a single proton uptake pathway to deliver protons both for O₂ reduction and for proton pumping. This pathway, called the K^C-pathway, starts at Glu-49^P in the accessory subunit CcoP, and connects into the catalytic subunit CcoN via the polar residues Tyr-(Y)-227, Asn (N)-293, Ser (S)-244, Tyr (Y)-321 and internal water molecules, and continues to the active site. However, although the residues are known to be functionally important, little is known about the mechanism and dynamics of proton transfer in the K^C-pathway. Here, we studied variants of Y227, N293 and Y321. Our results show that in the N293L variant, proton-coupled electron transfer is slowed during single-turnover oxygen reduction, and moreover it shows a pH dependence that is not observed in wildtype. This suggests that there is a shift in the pK_a of an internal proton donor into an experimentally accessible range, from >10 in wildtype to ~8.8 in N293L. Furthermore, we show that there are distinct roles for the conserved Y321 and Y227. In Y321F, proton uptake from bulk solution is greatly impaired, whereas Y227F shows wildtype-like rates and retains ~50% turnover activity. These tyrosines have evolutionary counterparts in the K-pathway of B-family HCuOs, but they do not have the same roles, indicating diversity in the proton transfer dynamics in the HCuO superfamily.     

A cluster of carboxylic groups in PsbO protein is involved in proton transfer from the water oxidizing complex of Photosystem II

The hypothesis presented here for proton transfer away from the water oxidation complex of Photosystem II (PSII) is supported by biochemical experiments on the isolated PsbO protein in solution, theoretical analyses of better understood proton transfer systems like bacteriorhodopsin and cytochrome oxidase, and the recently published 3D structure of PS II (Pdb entry 1S5L). We propose that a cluster of conserved glutamic and aspartic acid residues in the PsbO protein acts as a buffering network providing efficient acceptors of protons derived from substrate water molecules. The charge delocalization of the cluster ensures readiness to promptly accept the protons liberated from substrate water. Therefore protons generated at the catalytic centre of PSII need not be released into the thylakoid lumen as generally thought. The cluster is the beginning of a localized, fast proton transfer conduit on the lumenal side of the thylakoid membrane. Proton-dependent conformational changes of PsbO may play a role in the regulation of both supply of substrate water to the water oxidizing complex and the resultant proton transfer.     

Tuning the properties of M(diimine)(dithiolate) complexes - The role of the metal and solvent effect. A combined experimental, DFT and TDDFT study

A combined experimental and DFT/TDDFT study on a series of seven M(diimine)(dithiolate) complexes of Group VIII metals has been performed. This analysis focuses mainly on three aspects: (a) the role of the metal, (b) the connection of NLO properties, through the intrinsic hyperpolarizability, with the slope of solvatochromic plots and (c) the effect of solvation on the valence orbitals of the complexes. Besides, the molecular and the electronic structure as well as the bonding scheme of these complexes are also explored. The three aforementioned aspects were not satisfactorily clarified in the literature up to now. Moreover, as it is proved, they possess a central role in the experimental design of the reactions in which these complexes take place or in the design of the synthesis of compounds with pre-required properties.     

Role of protein motions on proton transfer pathways in human carbonic anhydrase II

We report here a theoretical study on the formation of long-range proton transfer pathways in proteins due to side chain conformational fluctuations of amino acid residues and reorganization of interior hydration positions. The proton transfer pathways in such systems may be modeled as fluctuating hydrogen-bonded networks with both short- and long-lived connections between the networked nodes, the latter being formed by polar protein atoms and water molecules. It is known that these fluctuations may extend over several decades of time ranging from a few femtoseconds to a few milliseconds. We have shown in this article how the use of a variety of theoretical methods may be utilized to detect a generic set of pathways and assess the feasibility of forming one or more transient connections. We demonstrate the application of these methods to the enzyme human carbonic anhydrase II and its mutants. Our results reveal several alternative pathways in addition to the one mediated by His-64. We also probe at length the mechanism of key conformational fluctuations contributing to the formation of the detected pathways.     

High resolution crystal structure of Paracoccus denitrificans cytochrome c oxidase: New insights into the active site and the proton transfer pathways

The structure of the two-subunit cytochrome c oxidase from Paracoccus denitrificans has been refined using X-ray cryodata to 2.25 Å resolution in order to gain further insights into its mechanism of action. The refined structural model shows a number of new features including many additional solvent and detergent molecules. The electron density bridging the heme a₃ iron and Cu_B of the active site is fitted best by a peroxo-group or a chloride ion. Two waters or OH⁻ groups do not fit, one water (or OH⁻) does not provide sufficient electron density. The analysis of crystals of cytochrome c oxidase isolated in the presence of bromide instead of chloride appears to exclude chloride as the bridging ligand. In the D-pathway a hydrogen bonded chain of six water molecules connects Asn131 and Glu278, but the access for protons to this water chain is blocked by Asn113, Asn131 and Asn199. The K-pathway contains two firmly bound water molecules, an additional water chain seems to form its entrance. Above the hemes a cluster of 13 water molecules is observed which potentially form multiple exit pathways for pumped protons. The hydrogen bond pattern excludes that the Cu_B ligand His326 is present in the imidazolate form.     

Glycocardiolipin modulates the surface interaction of the proton pumped by bacteriorhodopsin in purple membrane preparations

Glycocardiolipin is an archaeal analogue of mitochondrial cardiolipin, having an extraordinary affinity for bacteriorhodopsin, the photoactivated proton pump in the purple membrane of Halobacterium salinarum. Here purple membranes have been isolated by osmotic shock from either cells or envelopes of Hbt. salinarum. We show that purple membranes isolated from envelopes have a lower content of glycocardiolipin than standard purple membranes isolated from cells. The properties of bacteriorhodopsin in the two different purple membrane preparations are compared; although some differences in the absorption spectrum and the kinetic of the dark adaptation process are present, the reduction of native membrane glycocardiolipin content does not significantly affect the photocycle (M-intermediate rise and decay) as well as proton pumping of bacteriorhodopsin. However, interaction of the pumped proton with the membrane surface and its equilibration with the aqueous bulk phase are altered.     

How inter-subunit contacts in the membrane domain of complex I affect proton transfer energetics

The respiratory complex I is a redox-driven proton pump that employs the free energy released from quinone reduction to pump protons across its complete ca. 200?Å wide membrane domain. Despite recently resolved structures and molecular simulations, the exact mechanism for the proton transport process remains unclear. Here we combine large-scale molecular simulations with quantum chemical density functional theory (DFT) models to study how contacts between neighboring antiporter-like subunits in the membrane domain of complex I affect the proton transfer energetics. Our combined results suggest that opening of conserved Lys/Glu ion pairs within each antiporter-like subunit modulates the barrier for the lateral proton transfer reactions. Our work provides a mechanistic suggestion for key coupling effects in the long-range force propagation process of complex I.     

Reaction field effects on proton transfer in the active site of actinidin

The feasibility of the inclusion of reaction field effects in accurate ab initio self-consistent field-molecular orbital calculations was studied in the case of proton transfer in the active site of actinidin. The effects of the polarizability of the environment were included, using the direct reaction field model, which treats the environment as a set of interacting polarizable atoms. Up to 1000 of these atoms could be treated but about 300 were sufficient. The full geometry of the active site and the environment was taken into account. The stabilization of the ion pair was calculated to be 3.5 kcal. but this value may be 10 kcal depending on the geometry used. The effect of the static field from the long alpha-helix present in the enzyme was also studied. Dispersion effects are shown to be unimportant. The orientational polarizability of side chains and water molecules was not included.     

Interspecific hydrogen transfer during methanol degradation by Sporomusa acidovorans and hydrogenophilic anaerobes

In the presence of active hydrogenophilic sulfate-reducing bacteria, the homoacetogenic bacterium Sporomusa acidovorans did not produce acetate during methanol degradation. H₂S and presumably CO₂ were the only end products. Since the sulfate-reducer did not degrade methnol or acetate, the sulfidogenesis from methanol was related to a complete interspecific hydrogen transfer between both species.In coculture with hydrogenophilic methanogenic bacteria (Methanobacterium formicicum, Methanospirillum hungatei), the interspecific hydrogen transfer with S. acidovorans was incomplete. Beside CH₄ and presumably CO₂, acetate was produced. The results suggested that H₂-production and H₂-consumption were involved during anaerobic methanol degradation by S. acidovorans and the hydrogenophilic anaerobes play an important role during methanol degradation by homoacetogenic bacteria in anoxic environments.     

Toward an expanded oxygen atom transfer reactivity scale: Computational investigation of the energetics of oxo transfer reaction couples

The computational prediction of gas phase enthalpy (neutral substrates) and aqueous free energy (anion substrates) changes has been evaluated for the oxygen atom transfer reaction X + 1/2O₂ → XO. Several density functionals (SVWN, BP86, B3LYP) at double- and triple-ζ levels were surveyed, along with one composite ab initio method (G3(MP2)). Results are presented for extensive main group test sets for which experimental thermochemistry is available. In addition, several minimal reaction couples of the type [M^IVOL₂]/[M^VIO₂L₂] (M = Mo, W) have been examined. Overall, the results suggest a computational approach to the energetics of oxo transfer is feasible, potentially affording an expanded oxo transfer reactivity scale.