Leading coordinate analysis of reaction pathways in proton chain transfer: Application to a two-proton transfer model for the green fluorescent protein

The ‘leading coordinate’ approach to computing an approximate reaction pathway, with subsequent determination of the true minimum energy profile, is applied to a two-proton chain transfer model based on the chromophore and its surrounding moieties within the green fluorescent protein (GFP). Using an ab initio quantum chemical method, a number of different relaxed energy profiles are found for several plausible guesses at leading coordinates. The results obtained for different trial leading coordinates are rationalized through the calculation of a two-dimensional relaxed potential energy surface (PES) for the system. Analysis of the 2-D relaxed PES reveals that two of the trial pathways are entirely spurious, while two others contain useful information and can be used to furnish starting points for successful saddle-point searches. Implications for selection of trial leading coordinates in this class of proton chain transfer reactions are discussed, and a simple diagnostic function is proposed for revealing whether or not a relaxed pathway based on a trial leading coordinate is likely to furnish useful information.

4-(2-Hydroxyethyl) anilinium 3,5-dinitrobenzoate

In the title compound, C8H12NO+ C7H3N2O6-, the anilinium and hydroxyl protons of the cation result in N-H...O, N-H..(O,O) and O-H...O hydrogen-bonding interactions with carboxylate O atom acceptors, forming a two-dimensional network structure. An intermolecular C-H...O interaction is also present.

4-Chloroanilinium 2-carboxy-4,5-dichlorobenzoate

The structure of the 1:1 proton-transfer compound of 4-chloroaniline with 4,5-dichlorophthalic acid (DCPA), viz. C6H7ClN+ C8H3Cl2O4-, has been determined at 130 K. The non-planar hydrogen phthalate anions and the 4-chloroanilinium cations form two-dimensional O-H...O and N-H...O hydrogen-bonded substructures which have no peripheral extension. Between the sheets there are weak \p--\p associations between alternating cation--anion aromatic ring systems [shortest centroid separation, 3.735(4)A].

1,10-phenanthrolin-1-ium 2-carboxy-4,5-dichlorobenzoate

In the structure of the 1:1 proton-transfer compound of 1,10-phenanthroline with 4,5-dichlorophthalic acid, C12H9N2+ C8H3Cl2O4-, determined at 130 K, the 1,10-phenanthroline cation and the hydrogen 4,5-dichlorophthalate anion associate through a single N-H...O(carboxyl) hydrogen bond giving discrete units which have no extension except through a number of weak cation C-H...O(anion) associations and weak cation--anion aromatic ring pi-pi interactions [minimum centroid separation, 3.6815(12)A]. The anions are essentially planar [maximum deviation 0.214(1)A (a carboxyl O)] with the syn-related H atom of the carboxyl group forming a short intramolecular O-H...O(carboxyl) hydrogen bond.

Isopropylaminium 2-carboxy-4, 5-dichlorobenzoate

In the structure of the 1:1 proton-transfer compound of isopropylamine with 4,5-dichlorophthalic acid, C3H10N+·C8H3Cl2O4-, the three cation H-atom donors associate with three separate carboxyl O-atom anion acceptors, giving conjoint cyclic R44(12), R44(16) hydrogen-bonding cation-anion interactions in a one-dimensional ribbon structure. In the anions, the carboxyl groups lie slightly out of the plane of the benzene ring [maximum deviations = 0.439 (1) for a carboxylic acid O atom and 0.433 (1) Å for a carboxylate O atom]. However, the syn-related proton of the carboxylic acid group forms the common short intramolecular O-HOcarboxyl hydrogen bond.

4-(4-Nitrobenzyl)pyridinium 5-nitrosalicylate

In the title salt, C12H11N2O2+·C7H4NO5-, the cations and anions interact through asymmetric cyclic pyridinium-carboxylate N-HO,O' hydrogen-bonding associations [graph set R12(4)], giving discrete heterodimers having weak cation-anion - aromatic ring interactions [minimum ring centroid separation = 3.7116 (9) Å]

4-Methylpyridinium 2-carboxy-4,5-dichlorobenzoate

In the structure of the 1:1 proton-transfer compound of 4-methylpyridine (\g-picoline) with 4,5-dichlorophthalic acid, C6H8N+ C8H3Cl2O4- . H2O, determined at 200 K, the 4,5-dichlorophthalate anions are bridged by the water molecule through O--H...O~carboxyl~ hydrogen bonds, giving zig-zag chains which extend along the c axial direction of the unit cell. The 4-methylpyridine cations are linked to the chains through single N--H...O~water~ hydrogen bonds and occupy the voids within the chains in the one-dimensional structure. The anions have the common 'planar' conformation with the short intramolecular O--H...O(carboxyl) hydrogen bond.

4-Aminopyridinium cis-2-carboxycyclohexane-1-carboxylate

In the structure of the title compound, C5H7N2+ C8H11O4-, the cis-monoanions associate through short carboxylic acid-carboxyl O-H...O hydrogen bonds [graph set C(7)], forming zigzag chains which extend along c and are inter-linked through pyridinium and amine N-H...O(carboxyl) hydrogen bonds giving a three-dimensional network structure.

Ethane-1,2-diaminium 4,4'-sulfonyldibenzoate

In the title compound 2(C2H10N2+) C14H8O6S2-, both the ethylenediaminium cations and the 4,4'-sulfonyldibenzoate dianions have crystallographic two-fold rotational symmetry and are interlinked by aminium N-H...O(carboxyl) hydrogen-bonding associations giving sheets which extend along (101) and are further linked down b forming a three-dimensional structure.

4-(4-Nitrobenzyl)pyridinium 3-carboxy-4-hydroxybenzenesulfonate

In the title salt, C12H11N2O2+·C7H5O6S-, the dihedral angle between the benzene and pyridine rings in the 4-(4-nitrobenzyl)pyridinium cation is 82.7 (2)°. Within the anion there is an intramolecular hydroxy-O-HO(carboxylic acid) bond. In the crystal, the cation forms a single N+-HOsulfonate hydrogen bond with the anion. These cation-anion pairs interact through duplex anion carboxylic acid O-HOsulfonate hydrogen bonds, giving a centrosymmetric cyclic association [graph set R22(16)]. The crystals studied were non-merohedrally twinned.

Proton translocation coupled to electron transfer reactions in terminal oxidases

Terminal oxidases are the final proteins of the respiratory chain in eukaryotes and some bacteria. They catalyze most of the biological oxygen consumption on Earth done by aerobic organisms. During the catalytic reaction terminal oxidases reduce dioxygen to water and use the energy released in this process to maintain the electrochemical proton gradient by functioning as a redox-driven proton pump. This membrane gradient of protons is extremely important for cells as it is used for many cellular processes, such as transportation of substrates and ATP synthesis. Even though the structures of several terminal oxidases are known, they are not sufficient in themselves to explain the molecular mechanism of proton pumping. In this work we have applied a complex approach using a variety of different techniques to address the properties and the mechanism of proton translocation by the terminal oxidases. The combination of direct measurements of pH changes during catalytic turnover, time-resolved potentiometric electrometry and optical spectroscopy, made it possible to obtain valuable information about various aspects of oxidase functioning. We compared oxygen binding properties of terminal oxidases from the distinct heme-copper (CcO) and cytochrome bd families and found that cytochrome bd has a high affinity for oxygen, which is 3 orders of magnitude higher than that of CcO. Interestingly, the difference between CcO and cytochrome bd is not only in higher affinity of the latter to oxygen, but also in the way that each of these enzymes traps oxygen during catalysis. CcO traps oxygen kinetically - the molecule of bound dioxygen is rapidly reduced before it can dissociate. Alternatively, cytochrome bd employs an alternative mechanism of oxygen trapping - part of the redox energy is invested into tight oxygen binding, and the price paid for this is the lack of proton pumping. A single cycle of oxygen reduction to water is characterized by translocation of four protons across the membrane. Our results make it possible to assign the pumping steps to discrete transitions of the catalytic cycle and indicate that during in vivo turnover of the oxidase these four protons are transferred, one at a time, during the P→F, F→OH, Oh→Eh, and Eh→R transitions. At the same time, each individual proton translocation step in the catalytic cycle is not just a single reaction catalyzed by CcO, but rather a complicated sequence of interdependent electron and proton transfers. We assume that each single proton translocation cycle of CcO is assured by internal proton transfer from the conserved Glu-278 to an as yet unidentified pump site above the hemes. Delivery of a proton to the pump site serves as a driving reaction that forces the proton translocation cycle to continue.

Direct simulation of proton-coupled electron transfer reaction dynamics and mechanisms

Proton-coupled electron transfer (PCET) reactions are ubiquitous throughout chemistry and biology. However, challenges arise in both the the experimental and theoretical investigation of PCET reactions; the rare-event nature of the reactions and the coupling between quantum mechanical electron- and proton-transfer with the slower classical dynamics of the surrounding environment necessitates the development of robust simulation methodology. In the following dissertation, novel path-integral based methods are developed and employed for the direct simulation of the reaction dynamics and mechanisms of condensed-phase PCET.

Proton uptake controls electron transfer in cytochrome c oxidase

In cytochrome c oxidase, a requirement for proton pumping is a tight coupling between electron and proton transfer, which could be accomplished if internal electron-transfer rates were controlled by uptake of protons. During reaction of the fully reduced enzyme with oxygen, concomitant with the “peroxy” to “oxoferryl” transition, internal transfer of the fourth electron from CuA to heme a has the same rate as proton uptake from the bulk solution (8,000 s−1). The question was therefore raised whether the proton uptake controls electron transfer or vice versa. To resolve this question, we have studied a site-specific mutant of the Rhodobacter sphaeroides enzyme in which methionine 263 (SU II), a CuA ligand, was replaced by leucine, which resulted in an increased redox potential of CuA. During reaction of the reduced mutant enzyme with O2, a proton was taken up at the same rate as in the wild-type enzyme (8,000 s−1), whereas electron transfer from CuA to heme a was impaired. Together with results from studies of the EQ(I-286) mutant enzyme, in which both proton uptake and electron transfer from CuA to heme a were blocked, the results from this study show that the CuA → heme a electron transfer is controlled by the proton uptake and not vice versa. This mechanism prevents further electron transfer to heme a3–CuB before a proton is taken up, which assures a tight coupling of electron transfer to proton pumping.

Ethane-1,2-diaminium 4,5-dichlorophthalate

In the structure of the title compound, C2H10N22+·C8H2Cl2O42-, the dications and dianions form hydrogen-bonded ribbon substructures which enclose conjoint cyclic R21(7), R12(7) and R42(8) associations and extend down the c-axis direction. These ribbons inter-associate down b, giving a two-dimensional sheet structure. In the dianions, one of the carboxylate groups is essentially coplanar with the benzene ring, while the other is normal to it [C-C-C-O torsion angles = 177.67 (12) and 81.94 (17)°, respectively].

Bis[2-(2,4-dinitrobenzyl)pyridinium] biphenyl-4,4'-disulfonate

In the structure of the title compound, the salt 2(C12H10N3O4+) (C12H8O6S2)2- . 3H2O, determined at 173 K, the biphenyl-4,4'-disulfonate dianions lie across crystallographic inversion centres with the sulfonate groups interacting head-to-head through centrosymmetric cyclic bis(water)-bridged hydrogen-bonding associations [graph set R4/4(11)], forming chain structures. The 2-(2,4-dinitrobenzyl)pyridinium cations are linked to these chains through N+-H...O(water) hydrogen bonds and a two-dimensional network structure is formed through water bridges between sulfonate and 2-nitro O atoms, while the structure also has weak cation--anion pi-pi aromatic ring interactions [minimum ring centroid separation 3.8441(13)A].