Dissociation Constants of Weak Acids from ab Initio Molecular Dynamics Using Metadynamics: Influence of the Inductive Effect and Hydrogen Bonding on pK(a) Values

The theoretical estimation of the dissociation constant, or pK(a), of weak acids continues to be a challenging field. Here, we show that ab initio CarParrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free-energy profile of the dissociation reaction provide reasonable estimates of the pK(a) value. Water molecules, sufficient to complete the three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. The free-energy profiles exhibit two distinct minima corresponding to the dissociated and neutral states of the acid, and the difference in their values provides the estimate for pK(a). We show for a series of organic acids that CPMD simulations in conjunction with metadynamics can provide reasonable estimates of pK(a) values. The acids investigated were aliphatic carboxylic acids, chlorine-substituted carboxylic acids, cis- and trans-butenedioic acid, and the isomers of hydroxybenzoic acid. These systems were chosen to highlight that the procedure could correctly account for the influence of the inductive effect as well as hydrogen bonding on pK(a) values of weak organic acids. In both situations, the CPMD metadynamics procedure faithfully reproduces the experimentally observed trend and the magnitudes of the pK(a) values.

Kinetic study of solid phase demineralization by weak acids in one-step enzymatic bio-refinery of shrimp cuticles

We describe a one-step bio-refinery process for shrimp composites by-products. Its originality lies in a simple rapid (6 h) biotechnological cuticle fragmentation process that recovers all major compounds (chitins, peptides and minerals in particular calcium). The process consists of a controlled exogenous enzymatic proteolysis in a food-grade acidic medium allowing chitin purification (solid phase), and recovery of peptides and minerals (liquid phase). At a pH of between 3.5 and 4, protease activity is effective, and peptides are preserved. Solid phase demineralization kinetics were followed for phosphoric, hydrochloric, acetic, formic and citric acids with pKa ranging from 2.1 to 4.76. Formic acid met the initial aim of (i) 99 % of demineralization yield and (ii) 95 % deproteinization yield at a pH close to 3.5 and a molar ratio of 1.5. The proposed one-step process is proven to be efficient. To formalize the necessary elements for the future optimization of the process, two models to predict shell demineralization kinetics were studied, one based on simplified physical considerations and a second empirical one. The first model did not accurately describe the kinetics for times exceeding 30 minutes, the empirical one performed adequately.

Dissociation constants of weak organic acids in protic solvents obtained from their first hyperpolarizabilities in solution

The first hyperpolarizabilities (beta) of some weak aromatic organic acids have been measured in protic solvents by the hyper-Rayleigh scattering (HRS) technique at low concentrations. The measured hyperpolarizability (beta(m)) varies between the two extreme limits: the hyperpolarizability of the acid form (beta(HA)) at the lower side and that of the basic form (beta(A-)) at the higher side. The degree of dissociation (alpha) of the acid in a solvent is related to the measured hyperpolarizability, beta(m), by the following relationship: beta(m)(2)=(1-alpha)beta(HA)(2)+alpha beta(A-)(2). The calculated beta's including solvent effects in terms of an Onsager field do not reproduce the experimentally measured hyperpolarizabilities. Other solvent-induced effects like hydrogen bonding and van der Waals interactions seem to influence the first hyperpolarizability and, thus, indirectly the extent of dissociation of these weak acids in these protic solvents.

Phenylboronic acids in crystal engineering: Utility of the energetically unfavorable syn,syn-conformation in co-crystal design

Phenylboronic acids can exist, in principle, in three different conformers (syn,syn; syn,anti and anti,anti) with distinct energy profiles. In their native state, these compounds prefer the energetically favored syn, anti-conformation. In molecular complexes, however, the functionality exhibits conformational diversity. In this paper we report a series of co-crystals, with N-donor compounds, prepared by a design strategy involving the synthons based on the syn, syn-conformation of the boronic acid functionality. For this purpose, we employed compounds with the 1,2-diazo fragment (alprazolam, 1H-tetrazole, acetazolamide and benzotriazole), 1,10-phenanthroline and 2,2'-bipyridine for the co-crystallization experiments. However, our study shows that the mere presence of the 1,2-diazo fragment in the coformer does not guarantee the successful formation of co-crystals with a syn, syn-conformation of the boronic acid. [GRAPHICS] The -B(OH)(2) fragment makes unsymmetrical O-H center dot center dot center dot N heterosynthons with alprazolam (ALP) and 1,10-phenanthroline (PHEN). In the co-crystals of phenylboronic acids with 1H-tetrazole (TETR) and 2,2'-bipyridine (BPY), the symmetrical boronic acid dimer is the major synthon. In the BPY complex, boronic acid forms linear chains and the pyridine compound interacts with the lateral OH of boronic acid dimers that acts as a connector, thus forming a ladder structure. In the TETR complex, each heterocycle interacts with three boronic acids. While two boronic acids interact using the phenolic group, the third molecule generates O-H center dot center dot center dot N hydrogen bonds using the extra OH group, of -B(OH)(2) fragment, left after the dimer formation. Thus, although molecules were selected retrosynthetically with the 1,2-diazo fragment or with nearby hetero-atoms to induce co-crystal formation using the syn,syn-orientation of the -B(OH)(2) functionality, co-crystal formation is in fact selective and is probably driven by energy factors. Acetazolamide (ACET) contains self-complementary functional groups and hence creates stable homosynthons. Phenylboronic acids being weak competitors fail to perturb the homosynthons and hence the components crystallize separately. Therefore, besides the availability of possible hydrogen bond acceptors in the required position and orientation, the ability of the phenyl-boronic acid to perturb the existing interactions is also a prerequisite to form co-crystals. This is illustrated in the table below. In the case of ALP, PHEN and BPY, the native structures are stabilized by weak interactions and may be influenced by the boronic acid fragment. Thus phenylboronic acids can attain co-crystals with those compounds, wherein the cyclic O-H center dot center dot center dot N hydrogen bonds are stronger than the individual homo-interactions. This can lower the lattice energy of the molecular complex as compared with the individual crystals. [GRAPHICS] Phenylboronic acids show some selectivity in the formation of co-crystals with N-heterocycles. The differences in solubility of the components fall short to provide a possible reason for the selective formation of co-crystals only with certain compounds. These compounds, being weak acids, do not follow the Delta pK(a) analysis and hence fail to provide any conclusive observation. Theoretical results show that of the three conformers possible, the syn,anti conformer is the most stable. The relative stabilities of the three conformers syn,anti,syn,syn and anti,anti are 0.0, 2.18 and 3.14 kcal/mol, respectively. The theoretical calculations corroborate the fact that only energetically favorable synthons can induce the formation of heterosynthons, as in ALP and PHEN complexes. From a theoretical and structural analysis it is seen that phenylboronic acids will form interactions with those molecules wherein the heterocyclic and acidic fragments can interrupt the homosynthons. However, the energy profile is shallow and can be perturbed easily by the presence of competing functional groups (such as OH and COOH) in the vicinity. [GRAPHICS] .

Estimating successive pK(a) values of polyprotic acids from ab initio molecular dynamics using metadynamics: the dissociation of phthalic acid and its isomers

Estimation of the dissociation constant, or pK(a), of weak acids continues to be a central goal in theoretical chemistry. Here we show that ab initio Car-Parrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free energy profile of the dissociation reaction can provide reasonable estimates of the successive pK(a) values of polyprotic acids. We use the distance-dependent coordination number of the protons bound to the hydroxyl oxygen of the carboxylic group as the collective variable to explore the free energy profile of the dissociation process. Water molecules, sufficient to complete three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. Two distinct minima corresponding to the dissociated and un-dissociated states of the acid are observed and the difference in their free energy values provides the estimate for pK(a), the acid dissociation constant. We show that the method predicts the pK(a) value of benzoic acid in good agreement with experiment and then show using phthalic acid (benzene dicarboxylic acid) as a test system that both the first and second pK(a) values as well, as the subtle difference in their values for different isomers can be predicted in reasonable agreement with experimental data.

Characterization of groundwater chemistry under the influence of lithologic and anthropogenic factors along a climatic gradient in Upper Cauvery basin, South India

Hydrogeological and climatic effect on chemical behavior of groundwater along a climatic gradient is studied along a river basin. `Semi-arid' (500-800 mm of mean annual rainfall), `sub-humid' (800-1,200 mm/year) and `humid' (1,200-1,500 mm/year) are the climatic zones chosen along the granito-gneissic plains of Kabini basin in South India for the present analysis. Data on groundwater chemistry is initially checked for its quality using NICB ratio (<+/- 5 %), EC versus TZ+ (similar to 0.85 correlation), EC versus TDS and EC versus TH analysis. Groundwater in the three climatic zones is `hard' to `very hard' in terms of Ca-Mg hardness. Polluted wells are identified (> 40 % of pollution) and eliminated for the characterization. Piper's diagram with mean concentrations indicates the evolution of CaNaHCO3 (semi-arid) from CaHCO3 (humid zone) along the climatic gradient. Carbonates dominate other anions and strong acids exceeded weak acids in the region. Mule Hole SEW, an experimental watershed in sub-humid zone, is characterized initially using hydrogeochemistry and is observed to be a replica of entire sub-humid zone (with 25 wells). Extension of the studies for the entire basin (120 wells) showed a chemical gradient along the climatic gradient with sub-humid zone bridging semi-arid and humid zones. Ca/Na molar ratio varies by more than 100 times from semi-arid to humid zones. Semi-arid zone is more silicaceous than sub-humid while humid zone is more carbonaceous (Ca/Cl similar to 14). Along the climatic gradient, groundwater is undersaturated (humid), saturated (sub-humid) and slightly supersaturated (semi-arid) with calcite and dolomite. Concentration-depth profiles are in support of the geological stratification i.e., not approximate to 18 m of saprolite and similar to 25 m of fracture rock with parent gneiss beneath. All the wells are classified into four groups based on groundwater fluctuations and further into `deep' and `shallow' based on the depth to groundwater. Higher the fluctuations, larger is its impact on groundwater chemistry. Actual seasonal patterns are identified using `recharge-discharge' concept based on rainfall intensity instead of traditional monsoon-non-monsoon concept. Non-pumped wells have low Na/Cl and Ca/Cl ratios in recharge period than in discharge period (Dilution). Few other wells, which are subjected to pumping, still exhibit dilution chemistry though water level fluctuations are high due to annual recharge. Other wells which do not receive sufficient rainfall and are constantly pumped showed high concentrations in recharge period rather than in discharge period (Anti-dilution). In summary, recharge-discharge concept demarcates the pumped wells from natural deep wells thus, characterizing the basin.

Determination of arsenic species by capillary zone electrophoresis with large-volume field-amplified stacking injection

Determination of arsenic species by large-volume field amplified stacking injection-capillary zone electrophoresis (LV-FASI-CZE) is reported in this paper. Whole column injection was employed. The optimum buffer pH for the separation of weak acids was discussed. It was found that the optimum buffer to analyze the stacked arsenate (As(V)), monomethylarsonate (MMA), and dimethylarsinate (DMA) was 25 mm phosphate at pH 6.5. However, the optimum buffer to analyze the concentrated arsenite (As(III)) was 20 mm phosphate - 10 mm borate at pH 9.28. The limits of detection of the method developed were 0.026 mg/L for As(III), 0.023 mg/L for As(V), 0.043 mg/L for MMA, and 0.018 mg/L for DMA. An enrichment factor of 34-100 for several arsenic species was obtained. In the end, this method was applied to determine the arsenic concentration in the environmental reference materials to show the usefulness of the method developed.

Problems in the use of benzoic acid for estimating the internal pH of yeasts

Some yeasts have the peculiar ability to grow in the presence of weak acids at rather low pH. These conditions are predominant in preserved foods and beverages such as fruit concentrates, juices, wine, where these yeasts are responsible for spoilage. The main preservatives currently utilized by food industries are sorbic, propionic, benzoic acids and SO2. It is usually assumed that weak acids diffuse through the plasma membrane in the undissociated form. In the cytoplasm, where the pH is higher, dissociation occurs resulting in accumulation of the lipid-insoluble anion and internal acidification. This is probably a very general mechanism of preventing microbial growth in foods.

Electrogenic proton/L-glutamate symport in isolated asparagus sprengeri mesophyll cell

Medium' alkaliniiation occurred -lipon the addition of L-Glu to mechanically isolated Asparagus sprenger-i mesophyll cells suspended in 1 mM CaS04. Alkalinization resulted from the coupled entry of H+ and L-Glu anion into the cells. This H+ IL-Glu symport did not stimulate K+ efflux. K+ efflux has been observed during H~ lamino acid symport in other systems. The stimulation of K+ efflux by proton coupled symport is regarded as an indicator of a plasma membrane depolarizing electrogenic symport process. H+ IL-Glu symport in Asparagus sprengerimesophyl1 cells was investigated to determine whether or not the process was electrogenic. The rate of uptake of 0.25 11M 3H-MTPP+ ( Methyltriphenylphosphonium, methyl-3H ) is a probe for monitoring changes in the membrane potential. 3HMTPP+ uptake was reduced by K+ or CCCP, agents known to depolarize the membrane potential. Uptake of 3H-MTPP+ was also inhibited by L-Glu but not by D-Glu. Conversely, 10 mM external MTPP+ inhibited the uptake of 14C-U-LGlu. Simultaneous measurements of the rates of 14C-U-L-Glu uptake and L-Glu dependent H+ influx showed that the molar stoichiometry of H+ IL-Glu symport was 2 to 1. K+ or Na+ stimulated H+ efflux was completely inhibited by DCCD, DES, oligomycin and antimycin reagents which inhibit ATP driven H+ efflux. The H+ efflux \Vas also stimulate.d by the weak acids, butyric acid and acetic acid, which are known fo-aCidify the cytoplasm. This weak acid stimulated H+ efflux was also completely inhibited by oligomycin. It was calculated that net L-Glu dependent H+ influx increased by 100% in the presence of oligomycin and that despite net medium alkalinization H+ IL-Glu symport stimulates ATP dependent H+ efflux. 11 The data presented in this study indicate that H+ IL-Glu symport is electrogenic. The data also show that ATP dependent Ht efflux rather than K+ efflux is the- process compensating for thi~ electrogenic H+ IL-Glu symport.

The production of 4-adminobutyrate in response to treatments reducing cytosolic pH and the regulation of intracellular pH

GABA (4-aminobutyrate) is synthesized through the decarboxylation of LGlu- (L-Glu-+ H+ ---> GABA + C02), and compared to many free amino acids is present in high concentrations in plant cells. GABA levels rise rapidly and dramatically in response to varied stress conditions including anaerobiosis. Recent papers suggest that GABA production and associated H+ consumption are parts of a metabolic pH-stat mechanism which ameliorates the intracellular pH decline associated with anaerobiosis or other treatments. To test this hypothesis GABA production and efflux have been measured in isolated Asparagus sprengeri cells in response to three treatments which potentially cause intracellular acidification. Acid loads were imposed using 60 min of (i) anaerobiosis, (ii) H+/LGlu- cotransport, and (iii) treatment with permeant weak acids (butyric, acetic and propionic). Both intra- and extracellular GABA concentrations increased more than 100% after anaerobiosis, almost 1000% after H+/L-Glu- cotransport (light or dark) and almost 5000/0 after addition of 5 mM butyric acid at pH 5.0. HPLC analysis of amino acids indicates that as GABA concentrations increased in response to butyric acid addition, glutamate concentrations decreased. Time-course studies demonstrated that added butyric acid stimulates GABA production by 2800/0 within 15 seconds. A fluorescent determination of cytosolic pH indicates that addition of butyric or other weak acids resulted in a rapid reduction in cytosolic pH of 0.6 pH units. The half time for the response to butyric acid addition is 2.1 seconds, indicating that the decline in cytosolic pH is rapid enough to account for the rapid stimulation of GABA production. The acid load in response to butyric acid addition was assayed by measurements of 14C-butyric acid uptake. Calculations indicate that GABA production accounted for 45% of the imposed acid load. The biological significance of GABA efflux is not yet understood. The results support the original hypothesis suggesting a role for GABA production in cellular pH regulation.

Development of an Analytical Methodology for Quantification of Strong Acid in Diesel Oil

Many factors can affect the quality of diesel oil, in particular the degradation processes that are directly related to some organosulfur compounds. During the degradation process, these compounds are oxidized into their corresponding sulfonic acids, generating a strong acid content during the process. p-Toluene sulfonic acid analysis was performed using the linear sweep voltammetry technique with a platinum ultramicroelectrode in aqueous solution containing 3 mol L(-1) potassium chloride. An extraction step was introduced prior to the voltammetric detection in order to avoid the adsorption of organic molecules, which inhibit the electrochemical response. The extraction step promoted the transference of sulfonic acid from the diesel oil to an aqueous phase. The method was accurate and reproducible, with detection and quantification limits of 5 ppm and 15 ppm, respectively. Recovery of sulfonic acid was around 90%.

Evidence for a Slow-Turnover Form of the Ca2+-Independent Phosphoenolpyruvate Carboxylase Kinase in the Aleurone-Endosperm Tissue of Germinating Barley Seeds1

Phosphoenolpyruvate carboxylase (PEPC) activity was detected in aleurone-endosperm extracts of barley (Hordeum vulgare) seeds during germination, and specific anti-sorghum (Sorghum bicolor) C4 PEPC polyclonal antibodies immunodecorated constitutive 103-kD and inducible 108-kD PEPC polypeptides in western analysis. The 103- and 108-kD polypeptides were radiolabeled in situ after imbibition for up to 1.5 d in 32P-labeled inorganic phosphate. In vitro phosphorylation by a Ca2+-independent PEPC protein kinase (PK) in crude extracts enhanced the enzyme's velocity and decreased its sensitivity to l-malate at suboptimal pH and [PEP]. Isolated aleurone cell protoplasts contained both phosphorylated PEPC and a Ca2+-independent PEPC-PK that was partially purified by affinity chromatography on blue dextran-agarose. This PK activity was present in dry seeds, and PEPC phosphorylation in situ during imbibition was not affected by the cytosolic protein-synthesis inhibitor cycloheximide, by weak acids, or by various pharmacological reagents that had proven to be effective blockers of the light signal transduction chain and PEPC phosphorylation in C4 mesophyll protoplasts. These collective data support the hypothesis that this Ca2+-independent PEPC-PK was formed during maturation of barley seeds and that its presumed underlying signaling elements were no longer operative during germination.

A simple light-driven transmembrane proton pump.

Light-induced lipophilic porphyrin/aqueous acceptor charge separation across a single lipid-water interface can pump protons across the lipid bilayer when the hydrophobic weak acids, carbonylcyanide m-chlorophenylhydrazone and its p-trifluoromethoxyphenyl analogue, are present. These compounds act as proton carriers across lipid bilayers. In their symmetric presence across the bilayer, the positive currents and voltages produced by the photogeneration of porphyrin cations are replaced by larger negative currents and voltages. The maximum negative current and voltage occur at the pH of maximum dark conductance. The reversed larger current and voltage show a positive ionic charge transport in the same direction as the electron transfer. This transport can form an ion concentration gradient. The movement of protons is verified by an unusual D2O isotope effect that increases the negative ionic current by 2- to 3-fold. These effects suggest that an interfacial pK shift of the weak acid caused by the local electric field of photoformed porphyrin cations/acceptor anions functions as the driving force. The estimated pumping efficiency is 10-30%. Time-resolved results show that proton pumping across the bilayer occurs on the millisecond time scale, similar to that of biological pumps. This light-driven proteinless pump offers a simple model for a prebiological energy transducer.

Determination of the transiently lowered pKa of the retinal Schiff base during the photocycle of bacteriorhodopsin.

Reprotonation of the transiently deprotonated retinal Schiff base in the bacteriorhodopsin photocycle is greatly slowed when the proton donor Asp-96 is removed with site-specific mutagenesis, but its rate is restored upon adding azide or other weak acids such as formate and cyanate. As expected, between pH 3 and 7 the rate of Schiff base protonation in the photocycle of the D96N mutant correlates with the concentrations of the acid forms of these agents. Dissection of the rates in the biexponential reprotonation kinetics of the Schiff base between pH 7 and 9 yielded calculated rate constants for the protonation equilibrium. Their dependencies on pH and azide or cyanate concentrations are consistent with both earlier suggested mechanisms: (i) azide and other weak acids may function as proton carriers in the protonation equilibrium of the Schiff base, or (ii) the binding of their anionic forms may catalyze proton conduction to and from the Schiff base. The measured rate constants allow the calculation of the pKa of the Schiff base during its reprotonation in the photocycle of D96N. It is 8.2-8.3, a value much below the pKa determined earlier in unphotolyzed bacteriorhodopsin.