A semi-grand canonical Monte Carlo simulation model for ion binding to ionizable surfaces: proton binding of carboxylated latex particles as a case study

In this paper, we present a computer simulation study of the ion binding process at an ionizable surface using a semi-grand canonical Monte Carlo method that models the surface as a discrete distribution of charged and neutral functional groups in equilibrium with explicit ions modelled in the context of the primitive model. The parameters of the simulation model were tuned and checked by comparison with experimental titrations of carboxylated latex particles in the presence of different ionic strengths of monovalent ions. The titration of these particles was analysed by calculating the degree of dissociation of the latex functional groups vs. pH curves at different background salt concentrations. As the charge of the titrated surface changes during the simulation, a procedure to keep the electroneutrality of the system is required. Here, two approaches are used with the choice depending on the ion selected to maintain electroneutrality: counterion or coion procedures. We compare and discuss the difference between the procedures. The simulations also provided a microscopic description of the electrostatic double layer (EDL) structure as a function of p H and ionic strength. The results allow us to quantify the effect of the size of the background salt ions and of the surface functional groups on the degree of dissociation. The non-homogeneous structure of the EDL was revealed by plotting the counterion density profiles around charged and neutral surface functional groups.

A semi-grand canonical Monte Carlo simulation model for ion binding to ionizable surfaces: Proton binding of carboxylated latex particles as a case study

In this paper, we present a computer simulation study of the ion binding process at an ionizable surface using a semi-grand canonical Monte Carlo method that models the surface as a discrete distribution of charged and neutral functional groups in equilibrium with explicit ions modelled in the context of the primitive model. The parameters of the simulation model were tuned and checked by comparison with experimental titrations of carboxylated latex particles in the presence of different ionic strengths of monovalent ions. The titration of these particles was analysed by calculating the degree of dissociation of the latex functional groups vs. pH curves at different background salt concentrations. As the charge of the titrated surface changes during the simulation, a procedure to keep the electroneutrality of the system is required. Here, two approaches are used with the choice depending on the ion selected to maintain electroneutrality: counterion or coion procedures. We compare and discuss the difference between the procedures. The simulations also provided a microscopic description of the electrostatic double layer (EDL) structure as a function of pH and ionic strength. The results allow us to quantify the effect of the size of the background salt ions and of the surface functional groups on the degree of dissociation. The non-homogeneous structure of the EDL was revealed by plotting the counterion density profiles around charged and neutral surface functional groups. © 2011 American Institute of Physics.

Propriedades ácido-base e de complexação de ácidos húmico e fúlvico isolados de vermicomposto

Proton binding properties of humic and fulvic acids were studied by potentiometric titration. Carboxylic groups were the predominant ionizable sites in comparison to phenolic and amine groups. Total acidity of fulvic acid was 12 x 10-3 mol g-1, a number significantly higher than that obtained for humic acid (5.2 x 10-3 mol g-1). Copper ion binding was evaluated at pH 4, 5 and 6 by potentiometric titration with an ion selective electrode for Cu(II). Differential stability constants and complexation capacities were systematically higher for humic acid, despite its lower number of ionizable sites in comparison with fulvic acid.

Protonation dynamics of the extracellular and cytoplasmic surface of bacteriorhodopsin in the purple membrane.

The dynamics of proton binding to the extracellular and the cytoplasmic surfaces of the purple membrane were measured by laser-induced proton pulses. Purple membranes, selectively labeled by fluorescein at Lys-129 of bacteriorhodopsin, were pulsed by protons released in the aqueous bulk from excited pyranine (8-hydroxy-1,3,6-pyrenetrisulfonate) and the reaction of protons with the indicators was measured. Kinetic analysis of the data imply that the two faces of the membrane differ in their buffer capacities and in their rates of interaction with bulk protons. The extracellular surface of the purple membrane contains one anionic proton binding site per protein molecule with pK = 5.1. This site is within a Coulomb cage radius (approximately 15 A) from Lys-129. The cytoplasmic surface of the purple membrane bears 4-5 protonable moieties (pK = 5.1) that, due to close proximity, function as a common proton binding site. The reaction of the proton with this cluster is at a very fast rate (3.10(10) M-1.s-1). The proximity between the elements is sufficiently high that even in 100 mM NaCl they still function as a cluster. Extraction of the chromophore retinal from the protein has a marked effect on the carboxylates of the cytoplasmic surface, and two to three of them assume positions that almost bar their reaction with bulk protons. The protonation dynamics determined at the surface of the purple membrane is of relevance both for the vectorial proton transport mechanism of bacteriorhodopsin and for energy coupling, not only in halobacteria, but also in complex chemiosmotic systems such as mitochondrial and thylakoid membranes.

Electrostatic Contributions to the Thermodynamics of Ribonuclease P Protein Folding

Electrostatic interactions are of fundamental importance in determining the structure and stability of macromolecules. For example, charge-charge interactions modulate the folding and binding of proteins and influence protein solubility. Electrostatic interactions are highly variable and can be both favorable and unfavorable. The ability to quantify these interactions is challenging but vital to understanding the detailed balance and major roles that they have in different proteins and biological processes. Measuring pKa values of ionizable groups provides a sensitive method for experimentally probing the electrostatic properties of a protein.

pKa values report the free energy of site-specific proton binding and provide a direct means of studying protein folding and pH-dependent stability. Using a combination of NMR, circular dichroism, and fluorescence spectroscopy along with singular value decomposition, we investigated the contributions of electrostatic interactions to the thermodynamic stability and folding of the protein subunit of Bacillus subtilis ribonuclease P, P protein. Taken together, the results suggest that unfavorable electrostatics alone do not account for the fact that P protein is intrinsically unfolded in the absence of ligand because the pKa differences observed between the folded and unfolded state are small. Presumably, multiple factors encoded in the P protein sequence account for its IUP property, which may play an important role in its function.

Thermodynamic evaluation and modeling of proton and water exchange associated with benzamidine and berenil binding to ß-trypsin

Serine-proteases are involved in vital processes in virtually all species. They are important targets for researchers studying the relationships between protein structure and activity, for the rational design of new pharmaceuticals. Trypsin was used as a model to assess a possible differential contribution of hydration water to the binding of two synthetic inhibitors. Thermodynamic parameters for the association of bovine ß-trypsin (homogeneous material, observed 23,294.4 ± 0.2 Da, theoretical 23,292.5 Da) with the inhibitors benzamidine and berenil at pH 8.0, 25ºC and with 25 mM CaCl2, were determined using isothermal titration calorimetry and the osmotic stress method. The association constant for berenil was about 12 times higher compared to the one for benzamidine (binding constants are K = 596,599 ± 25,057 and 49,513 ± 2,732 M-1, respectively; the number of binding sites is the same for both ligands, N = 0.99 ± 0.05). Apparently the driving force responsible for this large difference of affinity is not due to hydrophobic interactions because the variation in heat capacity (DCp), a characteristic signature of these interactions, was similar in both systems tested (-464.7 ± 23.9 and -477.1 ± 86.8 J K-1 mol-1 for berenil and benzamidine, respectively). The results also indicated that the enzyme has a net gain of about 21 water molecules regardless of the inhibitor tested. It was shown that the difference in affinity could be due to a larger number of interactions between berenil and the enzyme based on computational modeling. The data support the view that pharmaceuticals derived from benzamidine that enable hydrogen bond formation outside the catalytic binding pocket of ß-trypsin may result in more effective inhibitors.

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.

Identification of the proton pathway in bacterial reaction centers: Inhibition of proton transfer by binding of Zn2+ or Cd2+

The reaction center (RC) from Rhodobacter sphaeroides converts light into chemical energy through the light induced two-electron, two-proton reduction of a bound quinone molecule QB (the secondary quinone acceptor). A unique pathway for proton transfer to the QB site had so far not been determined. To study the molecular basis for proton transfer, we investigated the effects of exogenous metal ion binding on the kinetics of the proton-assisted electron transfer kAB(2) (QA−•QB−• + H+ → QA(QBH)−, where QA is the primary quinone acceptor). Zn2+ and Cd2+ bound stoichiometrically to the RC (KD ≤ 0.5 μM) and reduced the observed value of kAB(2) 10-fold and 20-fold (pH 8.0), respectively. The bound metal changed the mechanism of the kAB(2) reaction. In native RCs, kAB(2) was previously shown to be rate-limited by electron transfer based on the dependence of kAB(2) on the driving force for electron transfer. Upon addition of Zn2+ or Cd2+, kAB(2) became approximately independent of the electron driving force, implying that the rate of proton transfer was reduced (≥ 102-fold) and has become the rate-limiting step. The lack of an effect of the metal binding on the charge recombination reaction D+•QAQB−• → DQAQB suggests that the binding site is located far (>10 Å) from QB. This hypothesis is confirmed by preliminary x-ray structure analysis. The large change in the rate of proton transfer caused by the stoichiometric binding of the metal ion shows that there is one dominant site of proton entry into the RC from which proton transfer to QB−• occurs.

Calcium is essential for the structural integrity of the cysteine-rich, ligand-binding repeat of the low-density lipoprotein receptor

Seven cysteine-rich repeats form the ligand-binding region of the low-density lipoprotein (LDL) receptor. Each of these repeats is assumed to bind a calcium ion, which is needed for association of the receptor with its ligands, LDL and beta-VLDL. The effects of metal ions on the folding of the reduced N-terminal cysteine-rich repeat have been examined by using reverse-phase high-performance liquid chromatography to follow the formation of fully oxidized isomers with different disulfide connectivities. in the absence of calcium many of the 15 possible isomers formed on oxidation, whereas in its presence the predominant product at equilibrium had the native disulfide bond connectivities. Other metals were far less effective at directing disulfide bond formation: Mn2+ partly mimicked the action of Ca2+, but Ba2+, Sr2+, and Mg2+ had little effect. This metal-ion specificity was also observed in two-dimensional H-1 NMR spectral studies: only Ca2+ induced the native three-dimensional fold. The two paramagnetic ions, Gd3+ and Mn2+, and Cd2+ did not promote adoption of a well-defined structure, and the two paramagnetic ions did not displace calcium ions. The location of calcium ion binding sites in the repeat was also explored by NMR spectroscopy. The absence of chemical shift changes for the side chain proton resonances of Asp26, Asp36, and Glu37 from pH 3.9 to 6.8 in the presence of calcium ions and their proximal location in the NMR structures implicated these side chains as calcium ligands. Deuterium exchange NMR experiments also revealed a network of hydrogen bonds that stabilizes the putative calcium-binding loop.

Microscopic-macroscopic approach for binding energies with the Wigner-Kirkwood method. II. Deformed nuclei

The binding energies of deformed even-even nuclei have been analyzed within the framework of a recently proposed microscopic-macroscopic model. We have used the semiclassical Wigner-Kirkwood ̄h expansion up to fourth order, instead of the usual Strutinsky averaging scheme, to compute the shell corrections in a deformed Woods-Saxon potential including the spin-orbit contribution. For a large set of 561 even-even nuclei with Z 8 and N 8, we find an rms deviation from the experiment of 610 keV in binding energies, comparable to the one found for the same set of nuclei using the finite range droplet model of Moller and Nix (656 keV). As applications of our model, we explore its predictive power near the proton and neutron drip lines as well as in the superheavy mass region. Next, we systematically explore the fourth-order Wigner-Kirkwood corrections to the smooth part of the energy. It is found that the ratio of the fourth-order to the second-order corrections behaves in a very regular manner as a function of the asymmetry parameter I=(N−Z)/A. This allows us to absorb the fourth-order corrections into the second-order contributions to the binding energy, which enables us us to simplify and speed up the calculation of deformed nuclei.

Proton translocation by cytochrome c oxidase reconstituted into proteoliposomes

Cytochrome c oxidase .inserted into proteoliposomes translocates protons with a stoichiometry of approx-, imately 0.4-0.6 H+/e- in the presence of valinomycin plus pottasium. The existance .ofsuchproton translocation is .supportedby experiments with lauryl maltoside which abolished the pulses but~~d not inhibit cyt. c binding .or oxidase turnover. Pulses with K3FeCN6 did not induce acidification further supporting vectorial proton transport by cyt ..aa3 . Upon lowering the ionic strength and pulsing with ferrocytochrome c, H+/eratios increased. This increase is attributed to scaler proton release consequent upon cyt.c-phospholipid binding. Oxygen pulses at low ionic strength however did not exhibit this large scaler increase in H+/e- ratios.A-small increase was observed upon .02 pul'sing at·low ionic strengt.h. This increase was KeN and, ,pcep sensitive and thus possibly due to a redox linked scaler deprotonation. Increases in the H+/e- ratio also occurred ifp~lses ,were performed in the presence of nonactin rather.than valinomycin. The fluorescent pH indicator pyranine was internally trapped inaa3 conta~ning "proteoliposomes. Internal alkalinization, as mon,itored by pyranine fluorescence leads to a of approx.imately 0.35 units, which is proportional to electron flux. This internal alkalinization was also DCCD sensitive, being inhibited by approximately 50%. This 50% inhibition of internal alkalinization supports the existance of vectorial proton transport.

Evaluation of the binding ability of a novel dioxatetraazamacrocyclic receptor that contains two phenanthroline units: selective uptake of carboxylate anions

The novel dioxatetraaza macrocycle [26]phen(2)N(4)O(2), which incorporates two phenanthroline units, has been synthesized, and its acid-base behavior has been evaluated by potentiometric and H-1 NMR methods. Six protonation constants were determined, and the protonation sequence was established by NMR. The location of the fifth proton on the phen nitrogen was confirmed by X-ray determinations of the crystal structures of the receptor as bromide and chloride salts. The two compounds have the general molecular formula {(H-5[26]phen(2)N(4)O(2))X-n(H2O)(5-n)}X(n-1)(.)mH(2)O, where X = Cl, n = 3, and m = 6 or X = Br, n = 4, and m = 5.5. In the solid state, the (H-5[26]phen(2)N(4)O(2))(5+) cation adopts a "horseshoe" topology with sufficient room to encapsulate three or four halogen anions through the several N-(HX)-X-... hydrogen-bonding interactions. Two supermolecules {(H-5[26]phen(2)N(4)O(2))X-n(H2O)(5-n)}((5-n)+) form an interpenetrating dimeric species, which was also found by ESI mass spectrum. Binding studies of the protonated macrocycle with aliphatic (ox(2-), mal(2-), suc(2-), cit(3-), cta(3-)) and aromatic (bzc(-), naphc(-), anthc(-), pyrc(-), ph(2-), iph(2-), tph(2-), btc(3-)) anions were determined in water by potentiometric methods. These studies were complemented by H-1 NMR titrations in D2O of the receptor with selected anions. The H-i[26]phen(2)N(4)O(2)(i+) receptor can selectively uptake highly charged or extended aromatic carboxylate anions, such as btc(3-) and pyrc(-), in the pH ranges of 4.0-8.5 and < 4.0, respectively, from aqueous solution that contain the remaining anions as pollutants or contaminants. To obtain further insight into these structural and experimental findings, molecular dynamics (MD) simulations were carried out in water solution.

Highly proton-ordered water structures on oxygen precovered Ru{0001}

We present helium scattering measurements of a water ad-layer grown on a O(2 1)/Ru(0001) surface. The adsorbed water layer results in a well ordered helium diffraction pattern with systematic extinctions of diffraction spots due to glide line symmetries. The data reflects a well-defined surface structure that maintains proton order even at surprisingly high temperatures of 140 K. The diffraction data we measure is consistent with a structure recently derived from STM measurements performed at 6 K. Comparison with recent DFT calculation is in partial agreement, suggesting that these calculations might be underestimating the contribution of relative water molecule orientations to the binding energy.

The Chameleon-Like Nature of Zwitterionic Micelles: Effect of Cation Binding

Addition of salts, especially perchlorates, to zwitterionic micelles of SB3-14, C(14)H(29)NMe(2)(+)(CH(2))(3)SO(3)(-), induces anionic character and uptake of H(3)O(+) by SB3-14 micelles. Thus, the addition of alkali metal perchlorates accelerates the acid hydrolysis of 2-(p-heptoxypheny1)-1,3-dioxolane, HPD, in the presence of SB3-14 micelles, which depends on the local proton concentration at the micelle surface. The addition of metal chlorides to solutions of such perchlorate-modified SB3-14 micelles decreases both the negative zeta potential of the micelles and the observed rate constant for acid hydrolysis of HPD. The effect of the monovalent cations Li(+), Na(+), and K(+) is smaller than that of the divalent cations Be(2+), Mg(2+), and Ca(2+), and much smaller than that of the trivalent cations Al(3+), La(3+), and Er(3+). The major factor responsible for this cation valence dependence of these effects is shown to be electrostatic in nature, reflecting the strong dependence of the micellar surface potential on the cation valence. The fact that the salt effects are not identical after correction for the electrostatic effects indicates that additional secondary nonelectrostatic effects may contribute as well.

Three strategically placed hydrogen-bonding residues convert a proton pump into a sensory receptor.

In haloarchaea, light-driven ion transporters have been modified by evolution to produce sensory receptors that relay light signals to transducer proteins controlling motility behavior. The proton pump bacteriorhodopsin and the phototaxis receptor sensory rhodopsin II (SRII) differ by 74% of their residues, with nearly all conserved residues within the photoreactive retinal-binding pocket in the membrane-embedded center of the proteins. Here, we show that three residues in bacteriorhodopsin replaced by the corresponding residues in SRII enable bacteriorhodopsin to efficiently relay the retinal photoisomerization signal to the SRII integral membrane transducer (HtrII) and induce robust phototaxis responses. A single replacement (Ala-215-Thr), bridging the retinal and the membrane-embedded surface, confers weak phototaxis signaling activity, and the additional two (surface substitutions Pro-200-Thr and Val-210-Tyr), expected to align bacteriorhodopsin and HtrII in similar juxtaposition as SRII and HtrII, greatly enhance the signaling. In SRII, the three residues form a chain of hydrogen bonds from the retinal's photoisomerized C(13)=C(14) double bond to residues in the membrane-embedded alpha-helices of HtrII. The results suggest a chemical mechanism for signaling that entails initial storage of energy of photoisomerization in SRII's hydrogen bond between Tyr-174, which is in contact with the retinal, and Thr-204, which borders residues on the SRII surface in contact with HtrII, followed by transfer of this chemical energy to drive structural transitions in the transducer helices. The results demonstrate that evolution accomplished an elegant but simple conversion: The essential differences between transport and signaling proteins in the rhodopsin family are far less than previously imagined.