Kinetic analysis of protonation of a specific site on a buffered surface of a macromolecular body

The kinetics of protonation of a specific site on a macromolecular structure (micelle) in buffered solution was studied with the purpose of evaluating the effect of buffer on the observed dynamics. The experimental system consisted of the following elements: Brij 58 micelles serving as homogeneous uncharged macromolecular bodies, bromocresol green, a well-adsorbed proton detector, and 2-naphthol-3,6-disulfonate as a proton emitter in the bulk. Imidazole was the mobile buffer while neutral red, which has a high affinity for the micellar surface, served as the immobile buffer. An intensive laser pulse ejects a proton from the proton emitter, and the subsequent proton-transfer reactions are measured by fast spectrophotometric methods. The dynamics of proton pulse in buffered solution are characterized by a very rapid trapping of the discharged protons by the abundant buffer molecules. This event has a major effect on the kinetic regime of the reaction. During the first 200 ns the proton flux is rate limited by free-proton diffusion. After this period, when the free-proton concentration decayed to the equilibrium level, the relaxation of the system is carried out by the diffusion of buffer. Thus in the buffered biochemical system, at neutral pH, most of proton flux between active sites and bulk is carried out by buffer molecules--not by diffusion of free protons. Surface groups on a high molecular weight body exchange protons among them at a very fast rate. This reaction has a major role on proton transfer from a specific site to the bulk.(ABSTRACT TRUNCATED AT 250 WORDS)     

Proton flux measurements from tissues in buffered solution

Proton movement across plant cell membranes is part of many important physiological processes. The net proton flux to or from tissues can be determined non-invasively by measuring the proton electrochemical potential gradient in the adjacent solution. In buffered solution, some of the protons crossing the tissue boundary diffuse as proto-nated buffer whose flux is not included in the flux calculated from the proton (hydrogen ion) electrochemical gradient. In this theoretical paper, it is shown how experimenters can calculate the protonated buffer flux from the measured proton flux in solution. The ratio of these two components of total proton flux depends on the pH of the solution and on the concentration and pK of the buffer. For a given concentration of a buffer which has a single pK, the flux ratio rises with pH when the solution pH is lower than the buffer pK. The slope is about 2 on a log₁₀ scale. As the pH increases above the pK, the flux ratio levels off to approach its maximum. With mixed buffers, or one having two or more pK values, the flux ratios are additive: each buffer acts independently based on its concentration and its pK value. Unbuffered solutions always have the buffering effects of water itself and also of carbonates due to carbon dioxide dissolved from the atmosphere. In unbuffered solutions at pH 6, the flux carried by water and carbonate is about 1 % of the measured proton flux. This validates measurements of proton flux from tissues, made by a number of workers, in unbuffered solutions below pH 6.     

Kinetic analysis of proton transfer between reactants adsorbed to the same micelle. The effect of proximity on the rate constants

The dense packing of protogenic enzymes on the coupling membrane can furnish a route for a rapid proton flux which may avoid the adjacent bulk phase. In order to evaluate the role of proximity between reactants on the rate constant of proton transfer we generated a model system consisting of 2-naphthol and pH indicator (bromocresol green) both adsorbed on the same micelle of unchanged detergent. Excitation of the 2-naphthol by a short intensive laser pulse lowers its pK with subsequent synchronized proton ejection. The discharged protons are detected by their reaction with the indicator using a fast transient absorption technique. Evidence is produced that under certain conditions all of the observed proton flux represents proton transfer between 2-naphthol and indicator molecules sharing the same micelle. In this model system the entire proton flux proceeds through an aqueous phase fully accessible to phosphate ions. The high proximity of the reactants (the separation can not exceed approximately 6 nm) has a marked effect on the rate constant of the reaction k = 2.0 +/- 0.5 X 10(11) M-1 s-1. In spite of this extremely fast rate of reaction we observe unhindered competition, for the surface discharged proton, between the surface-bound reactants and phosphate ions in the bulk. Thus even in proton transfer between closely packed reactants on an interface, the diffusion of the proton is not limited to the interface. This finding implies that on bioenergetic surface the electrochemical potential of the proton on the surface will equal that of the bulk.     

Proton coupling and the multiscale kinetic mechanism of a peptide transporter

Proton-coupled peptide transporters (POTs) are crucial for the uptake of di- and tripeptides as well as drug and prodrug molecules in prokaryotes and eukaryotic cells. We illustrate from multiscale modeling how transmembrane proton flux couples within a POT protein to drive essential steps of the full functional cycle: 1) protonation of a glutamate on transmembrane helix 7 (TM7) opens the extracellular gate, allowing ligand entry; 2) inward proton flow induces the cytosolic release of ligand by varying the protonation state of a second conserved glutamate on TM10; 3) proton movement between TM7 and TM10 is thermodynamically driven and kinetically permissible via water proton shuttling without the participation of ligand. Our results, for the first time, give direct computational confirmation for the alternating access model of POTs, and point to a quantitative multiscale kinetic picture of the functioning protein mechanism.     

Exceptional characteristics of amino proton exchange in guanosine compounds

Amino 1H NMR line width as a measure of amino proton exchange in guanosine compounds is completely unaffected by the addition of ca. 1 M tris(hydroxymethyl)-aminomethane, imidazole, 2-(N-morpholino)ethanesulfonic acid, glycine, or cacodylate, all shown to be effective buffer catalysts in adenosine and cytidine proton exchange. Line broadening, seen only with phosphate and acetate, is established by intermolecular interactions, as well as by amino to water proton exchange. This absence of buffer catalysis of exchange is accounted for by the relatively small implied effect of G(N-7) protonation on amino acidity, based on similar observations with 7-methylguanosine as a model for endocyclic protonation. The requirement for diffusion-controlled proton transfer in buffer catalysis is achieved by nucleobase protonation in adenine and cytosine, but not in guanine.     

Electrostatic Environment of Proteorhodopsin Affects the pKa of Its Buried Primary Proton Acceptor

The protonation state of embedded charged residues in transmembrane proteins (TMPs) can control the onset of protein function. It is understood that interactions between an embedded charged residue and other charged or polar residues in the moiety would influence its pKa, but how the surrounding environment in which the TMP resides affects the pKa of these residues is unclear. Proteorhodopsin (PR), a light-responsive proton pump from marine bacteria, was used as a model to examine externally accessible factors that tune the pKa of its embedded charged residue, specifically its primary proton acceptor D97. The pKa of D97 was compared between PR reconstituted in liposomes with different net headgroup charges and equilibrated in buffer with different ion concentrations. For PR reconstituted in net positively charged compared to net negatively charged liposomes in low-salt buffer solutions, a drop of the apparent pKa from 7.6 to 5.6 was observed, whereas intrinsic pKa modeled with surface pH calculated from Gouy-Chapman predictions found an opposite trend for the pKa change, suggesting that surface pH does not account for the main changes observed in the apparent pKa. This difference in the pKa of D97 observed from PR reconstituted in oppositely charged liposome environments disappeared when the NaCl concentration was increased to 150 mM. We suggest that protein-intrinsic structural properties must play a role in adjusting the local microenvironment around D97 to affect its pKa, as corroborated with observations of changes in protein side-chain and hydration dynamics around the E-F loop of PR. Understanding the effect of externally controllable factors in tuning the pKa of TMP-embedded charged residues is important for bioengineering and biomedical applications relying on TMP systems, in which the onset of functions can be controlled by the protonation state of embedded residues.     

Time-resolved protonation dynamics of a black lipid membrane monitored by capacitative currents

The laser-induced proton pulse (Gutman, M. (1986) Methods Enzymol. 127, 522-538) was used for transient protonation of one side of a black lipid membrane. The charging of the membrane drives an electric (voltage or current) signal selectively representing the fast proton exchange at the membrane/electrolyte interface. The sensitivity of the electric signal to the presence of buffer indicates that proton transfer is measured, not some dyes or membrane photoelectric artifact. The same event can be visualized in an analogous system consisting of a pH indicator adsorbed to neutral detergent-phospholipid mixed micelles. The time-resolved light absorption transient is equivalent to the electrically determined transient charging of the membrane surface. The sensitivity of the current measurement exceeds the spectrophotometric method by 6-8 orders of magnitudes. As little as 10(-18) mol of H+ reacting with 0.75 mm2 of the membrane surface can be monitored in a time-resolved observation. Both types of observed transients were accurately reconstructed by the numerical solution of coupled, non-linear, differential equations describing the system. The rate constants of the various proton transfer reactions were calculated and found to be of diffusion controlled reactions. There is no evidence for any barrier at the interface which either prevents protons from reaching the membrane, or keeps proton on the interface. The electric measurements can be applied for monitoring proton transfer kinetics of complex biomembrane preparations.     

Evaluation of linked protonation effects in protein binding reactions using isothermal titration calorimetry.

A theoretical development in the evaluation of proton linkage in protein binding reactions by isothermal titration calorimetry (ITC) is presented. For a system in which binding is linked to protonation of an ionizable group on a protein, we show that by performing experiments as a function of pH in buffers with varying ionization enthalpy, one can determine the pK(a)'s of the group responsible for the proton linkage in the free and the liganded states, the protonation enthalpy for this group in these states, as well as the intrinsic energetics for ligand binding (delta H(o), delta S(o), and delta C(p)). Determination of intrinsic energetics in this fashion allows for comparison with energetics calculated empirically from structural information. It is shown that in addition to variation of the ligand binding constant with pH, the observed binding enthalpy and heat capacity change can undergo extreme deviations from their intrinsic values, depending upon pH and buffer conditions.     

Structural transitions in polycytidylic acid: proton buffer capacity data

The pH-dependences of proton buffer capacity of poly(C) were computed on the basis of the literature data. In these curves there were observed four peaks: two narrow and two wide ones. The first narrow peak reflects the process of cooperative formation of double helices, which is induced by protonation of the N3 atom of nucleotide bases. The first wide peak is assigned to noncooperative process of poly(C) double helices protonation at the N3 nitrogen atom. It is proposed that the second wide peak corresponds to noncooperative protonation of the neutral cytosine bases at the oxygen atom. This reaction causes cooperative dissociation of the poly(C) double helices. The second narrow peak reflects the dissociation process.     

General acid-base catalysis in nucleobase amino proton exchange: cytidine

A useful property of DMSO solvent has been exploited to reveal a new catalytic route for cytidine amino proton exchange, relevant to exchange in the macromolecular state, but hidden in aqueous solution. Additional exchange mechanisms in aqueous monomeric cytidine (and adenosine) are obscured by the formation of a fast-exchanging endocyclic-protonated intermediate, which dominates the kinetics. Endocyclic nucleobase protonation could be circumvented in the presence of buffer conjugate acid by the use of DMSO/water solvent, permitting the first unequivocal observation buffer acid-catalyzed exchange from the neutral, unprotonated nucleobase, i.e., general acid catalysis. Because buffer ionization is greatly reduced in DMSO through anion desolvation, nucleobase protonation is suppressed in the presence of buffer acid. Evidence is presented to describe this catalytic route as one involving hydrogen bond formation between the buffer acid and the endocyclic protonation site, C(N-3). Since this same configuration is found in Watson-Crick hydrogen bonding, experiments are presented to demonstrate faster cytidine amino proton exchange with the formation of the G-C base pair in DMSO. The importance of this mechanism in past aqueous monomer studies and in the interpretation of macromolecular (DNA) hydrogen exchange is discussed.     

Interplay between local protein interactions and water bridging of a proton antenna carboxylate cluster

Proteins that bind protons at cell membrane interfaces often expose to the bulk clusters of carboxylate and histidine sidechains that capture protons transiently and, in proton transporters, deliver protons to an internal site. The protonation-coupled dynamics of bulk-exposed carboxylate clusters, also known as proton antennas, is poorly described. An essential open question is how water-mediated bridges between sidechains of the cluster respond to protonation change and facilitate transient proton storage. To address this question, here I studied the protonation-coupled dynamics at the proton-binding antenna of PsbO, a small extrinsinc subunit of the photosystem II complex, with atomistic molecular dynamics simulations and systematic graph-based analyses of dynamic protein and protein-water hydrogen-bond networks. The protonation of specific carboxylate groups is found to impact the dynamics of their local protein-water hydrogen-bond clusters. Regardless of the protonation state considered for PsbO, carboxylate pairs that can sample direct hydrogen bonding, or bridge via short hydrogen-bonded water chains, anchor to nearby basic or polar protein sidechains. As a result, carboxylic sidechains of the hypothesized antenna cluster are part of dynamic hydrogen bond networks that may rearrange rapidly when the protonation changes.     

General Acid-Base Catalysis in Nucleobase Amino Proton Exchange: Cytidine

Abstract

A useful property of DMSO solvent has been exploited to reveal a new catalytic route for cytidine amino proton exchange, relevant to exchange in the macromolecular state, but hidden in aqueous solution. Additional exchange mechanisms in aqueous monomeric cytidine (and adenosine) are obscured by the formation of a fast-exchanging endocyclic-protonated intermediate, which dominates the kinetics. Endocyclic nucleobase protonation could be circumvented in the presence of buffer conjugate acid by the use of DMSO/water solvent, permitting the first unequivocal observation buffer acid-catalyzed exchange from the neutral, unprotonated nucleobase, i.e., general acid catalysis. Because buffer ionization is greatly reduced in DMSO through anion desolvation, nucleobase protonation is supressed m the presence of buffer acid. Evidence is presented to describe this catalytic route as one involving hydrogen bond formation between the buffer acid and the endocyclic protonation site, C(N-3). Since this same configuration is found in Watson-Crick hydrogen bonding, experiments are presented to demonstrate faster cytidine amino proton exchange with the formation of the G-C base pair in DMSO. The importance of this mechanism in past aqueous monomer studies and in the interpretation of macromolecular (DNA) hydrogen exchange is discussed.     

Nonequilibration of membrane-associated protons with the internal aqueous space in dark-maintained chloroplast thylakoids

Isolated spinach thylakoids retain a slowly equilibrating pool of protons in the dark which are predominantly bound to buffering groups, probably amines, with low pKa values. We have measured the effects of permeant buffers, salts, sucrose, and uncouplers on the retention of the proton pool. Acetic anhydride, which reacts with neutral primary amine groups, was used to determine the protonation state of the amine buffering groups. It was previously shown by Bakeret al. that the extent of inhibition of photosystem II water-oxidizing capacity by acetic anhydride and the increase in derivatization by the anhydride are proportional to, and dependent on, the deprotonated state of the amine buffering pool. Therefore, acetic anhydride inhibition of water oxidation activity may be used as a measure of the protonation state of the amine buffering pool. By this method it is inferred that protons, in a metastable state, were retained by membranes suspended in high pH buffer for several hours in the dark. When both the internal and external aqueous phases were equilibrated with pH 8.8 buffer, the proton pool was released only upon addition of a protonophore. The osmotic strength of the suspension buffer affected uncoupler-induced proton release while ionic strength had little influence. The acetic anhydride-sensitive buffering group(s) of the water-oxidizing apparatus had an apparent pKa of 7.8. We conclude that an array of protein buffering groups reside either within the membrane matrix, or in proteins at the membrane surface, not in equilibrium with the bulk aqueous phases, and is responsible for the retention of the proton pool in dark maintained chloroplasts.     

Kinetics and mechanism of the acid transition of the active site in plastocyanin

Exchange on the microsecond time scale between the protonated and deprotonated forms of His92 in the copper site of reduced plastocyanin from the cyanobacteria Anabaena variabilis was monitored using 15N NMR relaxation measurements. On the basis of the dependence of the kinetics on pH and phosphate buffer concentration, we propose a two-step model for the protonation of the copper site in agreement with previous crystallographic studies. It is shown that the proton transfer is the rate-limiting step in the reaction at low buffer concentrations, whereas at high buffer concentrations, another step becomes rate-limiting. We suggest that the latter step is a concerted dissociation of His92 from the Cu(I) ion and a 180 degrees rotation of the imidazole ring, which precede the protonation. The first-order rate constant for the dissociation of His92 from the Cu(I) ion is estimated to be 2.4 x 10(4) s(-1). Also, a cooperative effect of the protonation of the remote His61 on the protonation of His92 and the redox properties of the protein was investigated by substituting His61 with asparagine. The mutation causes a modest change in both the pKa value of His92 and the redox potential of the protein.     

Gaugement of the inner space of the apomyoglobin's heme binding site by a single free diffusing proton. II. Interaction with a bulk proton.

The reaction mechanism and the dynamic aspects of protonation of a defined moiety located inside an aqueous cavity in a protein were monitored by time resolved spectroscopy using the pyranine apomyoglobin complex as a model (Shimoni, Tsfadia, Nachliel, and Gutman, 1993, Biophys. J. 64:472-479). The reaction was synchronized by a short laser pulse and the reprotonation of the ground state pyranine anion (phi O-) was monitored, in the microsecond time scale, by its transient absorption at 457 nm. The observed signal was reconstructed by a numeric solution of nonlinear, coupled differential equations which account for the direct reaction of phi O- with bulk proton and by proton transfer from the nearby amino acids: His 64, Asp 44, Asp 60, and Glu 59. A unique combination of rate constant was obtained which quantitates the contribution of each pathway to the overall relaxation process. In the first phase of the dynamics phi O- abstracts a proton from the nearby protonated histidine. The bulk proton interacts preferentially with the cluster of three carboxylates and immediately shuttled to the deprotonated histidine. The high proximity of the reactive groups and the strong electrostatic forces operating inside the heme binding cavity render the rate of proton transfer in the site ultrafast.     

Catalyzed decomposition of urea. Molecular dynamics simulations of the binding of urea to urease

We present the results of molecular dynamics simulations on the urea/urease system. The starting structure was prepared from the 2.0 A crystal structure of Benini et al. [(1999) Struct. Folding Des. 7, 205-216] of DAP-inhibited urease (PDB code ), and the trimeric structure (2479 residues) resulted in 180K atoms after solvation by water. The force field parameters were derived using the bonded model approach described by Hoops et al. [(1991) J. Am. Chem. Soc. 113, 8262-8270]. Three different systems were analyzed, each one modeling a different protonation pattern for the His320 and His219 residues. In each case, the three monomers of urease have been analyzed separately. The time-averaged structures observed in the three monomers suggest that urease could follow two different competitive mechanisms. A "protein-assisted proton transfer" mechanism points to Asp221 as crucial for catalysis. An "Asp-mediated proton transfer" involves the transfer of a proton from the bridging OH to an NH2 moiety of urea, assisted by Asp360 in the active site. The impact of the simulation results on our understanding of urease catalysis is discussed in detail.     

Unraveling the mechanism of proton translocation in the extracellular half‐channel of bacteriorhodopsin

Bacteriorhodopsin, a light activated protein that creates a proton gradient in halobacteria, has long served as a simple model of proton pumps. Within bacteriorhodopsin, several key sites undergo protonation changes during the photocycle, moving protons from the higher pH cytoplasm to the lower pH extracellular side. The mechanism underlying the long‐range proton translocation between the central (the retinal Schiff base SB216, D85, and D212) and exit clusters (E194 and E204) remains elusive. To obtain a dynamic view of the key factors controlling proton translocation, a systematic study using molecular dynamics simulation was performed for eight bacteriorhodopsin models varying in retinal isomer and protonation states of the SB216, D85, D212, and E204. The side‐chain orientation of R82 is determined primarily by the protonation states of the residues in the EC. The side‐chain reorientation of R82 modulates the hydrogen‐bond network and consequently possible pathways of proton transfer. Quantum mechanical intrinsic reaction coordinate calculations of proton‐transfer in the methyl guanidinium‐hydronium‐hydroxide model system show that proton transfer via a guanidinium group requires an initial geometry permitting proton donation and acceptance by the same amine. In all the bacteriorhodopsin models, R82 can form proton wires with both the CC and the EC connected by the same amine. Alternatively, rare proton wires for proton transfer from the CC to the EC without involving R82 were found in an O′ state where the proton on D85 is transferred to D212. Proteins 2016; 84:639–654. © 2016 Wiley Periodicals, Inc.     

Proton motive function of the terminal antiporter-like subunit in respiratory complex I

In the aerobic respiratory chains of many organisms, complex I functions as the first electron input. By reducing ubiquinone (Q) to ubiquinol, it catalyzes the translocation of protons across the membrane as far as ~200 Å from the site of redox reactions. Despite significant amount of structural and biochemical data, the details of redox coupled proton pumping in complex I are poorly understood. In particular, the proton transfer pathways are extremely difficult to characterize with the current structural and biochemical techniques. Here, we applied multiscale computational approaches to identify the proton transfer paths in the terminal antiporter-like subunit of complex I. Data from combined classical and quantum chemical simulations reveal for the first time structural elements that are exclusive to the subunit, and enables the enzyme to achieve coupling between the spatially separated Q redox reactions and proton pumping. By studying long time scale protonation and hydration dependent conformational dynamics of key amino acid residues, we provide novel insights into the proton pumping mechanism of complex I.     

Subsecond proton-hole propagation in bacteriorhodopsin

The dynamics of proton transfer between the surface of purple membrane and the aqueous bulk have recently been investigated by the Laser Induced Proton Pulse Method. Following a Delta-function release of protons to the bulk, the system was seen to regain its state of equilibrium within a few hundreds of microseconds. These measurements set the time frame for the relaxation of any state of acid-base disequilibrium between the bacteriorhodopsin's surface and the bulk. It was also deduced that the released protons react with the various proton binding within less than 10 micro s. In the present study, we monitored the photocycle and the proton-cycle of photo-excited bacteriorhodopsin, in the absence of added buffer, and calculated the proton balance between the Schiff base and the bulk phase in a time-resolved mode. It was noticed that the late phase of the M decay (beyond 1 ms) is characterized by a slow (subsecond) relaxation of disequilibrium, where the Schiff base is already reprotonated but the pyranine still retains protons. Thus, it appears that the protonation of D96 is a slow rate-limiting process that generates a "proton hole" in the cytoplasmic section of the protein. The velocity of the hole propagation is modulated by the ionic strength of the solution and by selective replacements of charged residues on the interhelical loops of the protein, at domains that seems to be remote from the intraprotein proton conduction trajectory.