Carbohydrate-pi interactions: what are they worth?

Protein-carbohydrate interactions play an important role in many biologically important processes. The recognition is mediated by a number of noncovalent interactions, including an interaction between the alpha-face of the carbohydrate and the aromatic side chain of the protein. To elucidate this interaction, it has been studied in the context of a beta-hairpin in aqueous solution, in which the interaction can be investigated in the absence of other cooperative noncovalent interactions. In this beta-hairpin system, both the aromatic side chain and the carbohydrate were varied in an effort to gain greater insight into the driving force and magnitude of the carbohydrate-pi interaction. The magnitude of the interaction was found to vary from -0.5 to -0.8 kcal/mol, depending on the nature of the aromatic ring and the carbohydrate. Replacement of the aromatic ring with an aliphatic group resulted in a decrease in interaction energy to -0.1 kcal/mol, providing evidence for the contribution of CH-pi interactions to the driving force. These findings demonstrate the significance of carbohydrate-pi interactions within biological systems and also their utility as a molecular recognition element in designed systems.     

The recognition of nucleotides with model beta-hairpin receptors: investigation of critical contacts and nucleotide selectivity

[reaction: see text] We have investigated the factors that contribute to binding of ATP by a designed 12-residue beta-hairpin peptide, WKWK, and have determined its selectivity for binding to the naturally occurring nucleotide triphosphates. We have previously shown that WKWK creates an ATP binding pocket on one face of the beta-hairpin consisting of two Trp and two Lys residues. Mutation of the two Lys residues on the binding face of the beta-hairpin resulted in a lower affinity, indicating that each is involved in ATP binding and that each residue contributes approximately -1.5 kcal/mol to the energy of complexation. Replacement of either Trp residue of the ATP binding pocket with Phe or Leu destabilizes the complex formed with ATP by approximately 1 kcal/mol, indicating that both Trp residues participate in interactions with ATP. For binding to the nucleotide triphosphates, the order of binding affinity was shown to follow dTTP > GTP > ATP > CTP, with differences in binding energies spanning as much as 1.6 kcal/mol. NMR analysis demonstrates that both aromatic interactions with the Trp side chains and CH-pi interactions between the ribose protons and the Trp residues may contribute significantly to binding. The results from our model system provide useful thermodynamic information regarding protein-nucleic acid interactions that occur at the surface of a beta-sheet.     

Stereodynamics and edge-to-face CH-pi aromatic interactions in o-phenethyl-substituted biaryls

As part of a search for systems that exhibit intramolecular aromatic edge-to-face interactions, a series of four biaryl compounds containing a phenethyl side chain were prepared. These compounds exist as a slowly interconverting mixture of two atropisomers due to steric hindrance to rotation about the biaryl bond. The more thermodynamically stable isomer exhibits an abnormal shielding of an ortho-proton in solution indicative of an edge-to-face CH-pi interaction with the terminal phenyl ring on the side chain. This tilted-T type of geometry was observed in the X-ray crystal structure of one of the compounds. The edge-to-face conformation in solution is estimated by variable temperature NMR studies to be energetically favored by ca. 1.6 kcal mol(-1) but entropically disfavored by ca. 5.0 cal K(-1) mol(-1).     

A designed beta-hairpin peptide for molecular recognition of ATP in water

A designed 12-residue beta-hairpin peptide with a diagonal tryptophan (Trp) pair was shown to bind ATP in water through a combination of aromatic and electrostatic interactions. The affinity for ATP was 5800 M-1 (DeltaG approximately -5.0 kcal/mol), a remarkable affinity for a short, structured peptide in water, consisting of entirely natural amino acid residues. Proton NMR measurements indicate that the adenine ring of the nucleotide is intercalated between the diagonal tryptophans in the bound state. Delineation of the contributions to ATP binding to the hairpin suggest that aromatic interactions contribute approximately -1.8 kcal/mol, while individual electrostatic interactions involving the ATP phosphates and positively charged side chains of the hairpin contribute approximately -1 kcal/mol each. The designed beta-hairpin receptor presents a novel minimalist system to investigate the energetic contributions to protein-nucleic acid recognition through the surface of a beta-sheet.     

Contribution of aromatic interactions to alpha-helix stability

The influence of natural and unnatural i, i + 4 aromatic side chain-side chain interactions on alpha-helix stability was determined in Ala-Lys host peptides by circular dichroism (CD). All interactions investigated provided some stability to the helix; however, phenylalanine-phenylalanine (F-F) and phenylalanine-pentafluorophenylalanine (F-f5F) interactions resulted in the greatest enhancement in helicity, doubling the helical content over i, i + 5 control peptides at internal positions. Quantification of these interactions using AGADIR multistate helix-coil algorithm revealed that the F-F and F-f5F interaction energies are equivalent at internal positions in the sequence (deltaGF-F = deltaGF-f5F = -0.27 kcal/mol), despite the differences in their expected geometries. As the strength of a face-to-face stacked phenyl-pentafluorophenyl interaction should surpass an edge-to-face or offset-stacked phenyl-phenyl interaction, we believe this result reflects the inability of the side chains in F-f5F to attain a fully stacked geometry within the context of an alpha-helix. Positioning the interactions at the C-terminus led to much stronger interactions (deltaGF-F = -0.8 kcal/mol; deltaGF-f5F = -0.55 kcal/mol) likely because of favorable chi(1) rotameric preferences for aromatic residues at C-capping regions of alpha-helices, suggesting that aromatic side chain-side chain interactions are an effective alpha-helix C-capping method.     

Aromaticity in metallabenzenes

The electronic structure and bonding situation in 21 metallabenzenes (metal=Os, Ru, Ir, Rh, Pt, and Pd) were investigated at the DFT level (BP86/TZ2P) by using an energy decomposition analysis (EDA) of the interaction energy between various fragments. The aim of the work is to estimate the strength of the pi bonding and the aromatic character of the metallacyclic compounds. Analysis of the electronic structure shows that the metallacyclic moiety has five occupied pi orbitals, two with b1 symmetry and three with a2 symmetry, which describe the pi-bonding interactions. The metallabenzenes are thus 10 pi-electron systems. This holds for 16-electron and for 18-electron complexes. The pi bonding in the metallabenzenes results mainly from the b1 contribution, but the a2 contribution is not negligible. Comparison of the pi-bonding strength in the metallacyclic compounds with acylic reference molecules indicates that metallabenzenes should be considered as aromatic compounds whose extra stabilization due to aromatic conjugation is weaker than in benzene. The calculated aromatic stabilization energies (ASEs) are between 8.7 kcal mol(-1) for 13 and 37.6 kcal mol(-1) for 16 which is nearly as aromatic as benzene (ASE=42.5 kcal mol(-1)). The classical metallabenzene model compounds 1 and 4 exhibit intermediate aromaticity with ASE values of 33.4 and 17.6 kcal mol(-1). The greater stability of the 5d complexes compared with the 4d species appears not to be related to the strength of pi conjugation. From the data reported here there is no apparent trend or pattern which indicates a correlation between aromatic stabilization and particular ligands, metals, coordination numbers or charge. The lower metal-C5H5 binding energy of the 4d complexes correlates rather with weaker sigma-orbital interactions.     

Absolute helical arrangement of sulfonamide in the crystal

[reaction: see text] 1,2-Bis(N-benzenesulfonyl-N-methylamino)benzene (2), which has no fixed asymmetric element, was crystallized from ethyl acetate as chiral crystals belonging to space group P4(1)2(1)2 (No. 92) or P4(3)2(1)2 (No. 96). The array of molecules built by the CH-pi interaction along the c-axis forms an enantiomeric helical superstructure in each individual crystal. The absolute configurations of the chiral crystals of 2 were determined by X-ray crystal structure analysis using the Flack parameter method. The solid-state CD spectra of the chiral crystals in KBr were mirror images. The equilibrium between the two enantiomers in solution is fast during crystallization at ambient temperature, and the energy barrier (DeltaG()) is estimated to be 11.7 +/- 0.3 kcal/mol (233 K).     

Probing phenylalanine/adenine pi-stacking interactions in protein complexes with explicitly correlated and CCSD(T) computations

To examine the effects of pi-stacking interactions between aromatic amino acid side chains and adenine bearing ligands in crystalline protein structures, 26 toluene/(N9-methyl)adenine model configurations have been constructed from protein/ligand crystal structures. Full geometry optimizations with the MP2 method cause the 26 crystal structures to collapse to six unique structures. The complete basis set (CBS) limit of the CCSD(T) interaction energies has been determined for all 32 structures by combining explicitly correlated MP2-R12 computations with a correction for higher-order correlation effects from CCSD(T) calculations. The CCSD(T) CBS limit interaction energies of the 26 crystal structures range from -3.19 to -6.77 kcal mol (-1) and average -5.01 kcal mol (-1). The CCSD(T) CBS limit interaction energies of the optimized complexes increase by roughly 1.5 kcal mol (-1) on average to -6.54 kcal mol (-1) (ranging from -5.93 to -7.05 kcal mol (-1)). Corrections for higher-order correlation effects are extremely important for both sets of structures and are responsible for the modest increase in the interaction energy after optimization. The MP2 method overbinds the crystal structures by 2.31 kcal mol (-1) on average compared to 4.50 kcal mol (-1) for the optimized structures.     

Influence of water-protein hydrogen bonding on the stability of Trp-cage miniprotein. A comparison between the TIP3P and TIP4P-Ew water models

We report extensive replica exchange molecular dynamics (REMD) simulations on the folding/unfolding equilibrium of Trp-cage miniprotein using the Amber ff99SB all atom forcefield and TIP3P and TIP4P-Ew explicit water solvent models. REMD simulation-lengths in the 500 ns to the microsecond regime per replica are required to adequately sample the folding/unfolding equilibrium. We observe that this equilibrium is significantly affected by the choice of the water model. Compared with experimental data, simulations using the TIP3P solvent describe the stability of the Trp-cage quite realistically, providing a melting point which is just a few Kelvins above the experimental transition temperature of 317 K. The TIP4P-Ew model shifts the equilibrium towards the unfolded state and lowers the free energy of unfolding by about 3 kJ mol(-1) at 280 K, demonstrating the need to fine-tune the protein-forcefield depending on the chosen water model. We report evidence that the main difference between the two water models is mostly due to the different solvation of polar groups of the peptide. The unfolded state of the Trp-cage is stabilized by an increasing number of hydrogen bonds, destabilizing the α-helical part of the molecule and opening the R-D salt bridge. By reweighting the strength of solvent-peptide hydrogen bonds by adding a hydrogen bond square well potential, we can fully recover the effect of the different water models and estimate the shift in population as due to a difference in hydrogen bond-strength of about 0.4 kJ mol(-1) per hydrogen bond.     

Van der Waals interactions in aromatic systems: structure and energetics of dimers and trimers of pyridine.

Full geometry optimizations at the dispersion-corrected DFT-BLYP level of theory were carried out for dimers and trimers of pyridine. The DFT-D interaction energies were checked against results from single-point SCS-MP2/aug-cc-pVTZ calculations. Three stacked structures and a planar H-bonded dimer were found to be very close in energy (interaction energies in the range from -3.4 to -4.0 kcal mol(-1)). Two T-shaped geometries are higher lying, by about 1 kcal mol(-1), which is explained by the more favorable electrostatic interactions in the stacked and H-bonded arrangements. The DFT-D approach has proved to be a reliable and efficient tool to explore the conformational space of aromatic van der Waals complexes and furthermore provides interaction energies with errors of less than 10-20 % of DeltaE. Comparisons with previous results obtained by using only partially optimized model geometries strongly indicate that unconstrained optimizations are mandatory in such weakly bonded low-symmetry systems.     

Controlled aggregation of adenine by sugars: physicochemical studies, molecular modelling simulations of sugar-aromatic CH-pi stacking interactions, and biological significance

CH-Pi stacking interactions between carbohydrates and aromatic compounds play a central role in biomolecular recognition, especially in lectin-sugar and protein-glycolipid systems. In the present study, we have measured the solubility of the sparingly soluble aromatic base adenine in presence of various saccharides as an approach to investigate the interaction between adenine and sugars. Above 82.5 mM, adenine solutions gradually formed a crystalline precipitate which could be quantified by spectrophotometric turbidity measurements. Precipitation of adenine was increased by salts (NaCl and NaF) whereas it was prevented by DMSO, in agreement with the involvement of hydrophobic interactions (pi-pi stacking) in the vertical stacking of adenine molecules. Several monosaccharides and disaccharides were found to increase adenine solubility, with the following order: D-galactose = D-lactose > D-sucrose > D-glucose = D-maltose > D-ribose > D-fructose. Molecular mechanics simulations indicated that the potent cosolvent effect of beta-D-galactopyranose was probably mediated by CH-pi stacking interactions between its apolar surface and the aromatic structure of adenine. The polar OH groups of the sugars interacted with surrounding water molecules, ensuring the solubility of sugar-adenine complexes. In contrast, beta-D-fructofuranose, which has two polar faces, did not stack onto adenine and had a weak cosolvent effect. CH-pi stacking interactions were also demonstrated between 6-methylpurine and the sugar head group of glycolipids (glucosyl-, galactosyl- and lactosylceramide) but not with the charged head group of phosphatidylinositol-4,5-diphosphate. These data indicate that galactose-containing molecules have a high stacking propensity for aromatic compounds such as adenine, due to the specific structure of the galactose cycle.     

Role of aromatic side chains in the folding and thermodynamic stability of integral membrane proteins

Aromatic residues are frequently found in helical and beta-barrel integral membrane proteins enriched at the membrane-water interface. Although the importance of these residues in membrane protein folding has been rationalized by thermodynamic partition measurements using peptide model systems, their contribution to the stability of bona fide membrane proteins has never been demonstrated. Here, we have investigated the contribution of interfacial aromatic residues to the thermodynamic stability of the beta-barrel outer membrane protein OmpA from Escherichia coli in lipid bilayers by performing extensive mutagenesis and equilibrium folding experiments. Isolated interfacial tryptophanes contribute -2.0 kcal/mol, isolated interfacial tyrosines contribute -2.6 kcal/mol, and isolated interfacial phenylalanines contribute -1.0 kcal/mol to the stability of this protein. These values agree well with the prediction from the Wimley-White interfacial hydrophobicity scale, except for tyrosine residues, which contribute more than has been expected from the peptide models. Double mutant cycle analysis reveals that interactions between aromatic side chains become significant when their centroids are separated by less than 6 A but are nearly insignificant above 7 A. Aromatic-aromatic side chain interactions are on the order of -1.0 to -1.4 kcal/mol and do not appear to depend on the type of aromatic residue. These results suggest that the clustering of aromatic side chains at membrane interfaces provides an additional heretofore not yet recognized driving force for the folding and stability of integral membrane proteins.     

Electronic control of amide cis-trans isomerism via the aromatic-prolyl interaction

The cis-trans isomerization of prolyl amide bonds results in large structural and functional changes in proteins and is a rate-determining step in protein folding. We describe a novel electronic strategy to control cis-trans isomerization, based on the demonstration that interactions between aromatic residues and proline are tunable by aromatic electronics. A series of peptides of sequence TXPN, X = Trp, pyridylalanine, pentafluorophenylalanine, or 4-Z-phenylalanine derivatives (Z = electron-donating, electron-withdrawing, or electron-neutral substituents), was synthesized and Ktrans/cis analyzed by NMR. Electron-rich aromatic residues stabilized cis amide bond formation, while electron-poor aromatics relatively favored trans amide bond formation. A Hammett correlation between aromatic electronics and cis-trans isomerization was observed. These results indicate that the interaction between aromatic residues and proline, which is observed to stabilize cis amide bonds and is also a general stabilizing interaction ubiquitous in proteins and protein-protein complexes, is not stabilized exclusively by a classical hydrophobic effect. To a large extent, the aromatic-prolyl interaction is driven and controllable by an electronic effect between the aromatic ring pi-electrons and the proline ring, consistent with a C-H-pi interaction as the key stabilizing force. The aromatic-prolyl interaction is electronically tunable by 0.9 kcal/mol and is enthalpic in nature. In addition, by combining aromatic ring electronics and stereoelectronic effects using 4-fluoroprolines, we demonstrate broad tuning (2.0 kcal/mol) of cis-trans isomerism in tetrapeptides. We demonstrate a simple tetrapeptide, TWflpN, that exhibits 60% cis amide bond and adopts a type VIa1 beta-turn conformation.     

The fragment molecular orbital method for geometry optimizations of polypeptides and proteins

The fragment molecular orbital method (FMO) has been used with a large number of wave functions for single-point calculations, and its high accuracy in comparison to ab initio methods has been well established. We have developed the analytic derivative of the electrostatic interaction between far separated fragments and performed a number of restricted Hartree-Fock (RHF) geometry optimizations using FMO and ab initio methods. In particular, the alpha-helix, beta-turn, and extended conformers of a 10-residue polyalanine were studied and the good FMO accuracy was established (the rms deviations for the former two forms were about 0.2 A and for the latter structure about 0.001 A). Met-enkephalin dimer was used as a model for the polypeptide binding and computed at the 3-21G and 6-31G* levels with a similar accuracy achieved; the error in the binding energy predictions (FMO vs ab initio) was 1-3 kcal/mol. Chignolin (PDB: 1uao) and an agonist polypeptide of the erythropoietin receptor protein (emp1) were optimized at the 3-21(+)G level, with the rms deviation from ab initio of about 0.2 A, or 0.5 degrees in terms of bond angles. The effect of solvation on the structure optimization was studied in chignolin and the Trp-cage miniprotein construct (PDB:1l2y), by describing water with TIP3P. The computed structures in gas phase and solution are compared to each other and experiment.     

Cooperativity, partially bound states, and enthalpy-entropy compensation

Efforts to develop a quantitative understanding of molecular recognition rely on the additivity of individual intermolecular interactions, and cooperativity represents one of the major potential stumbling blocks. A chemical double-mutant cycle has been used to experimentally measure cooperativity between functional group interactions within a complex framework. The interaction between two aromatic groups varies by 0.2 +/- 0.4 kJ mol(-1) in synthetic H-bonded complexes that differ by 8-13 kJ mol(-1) in overall stability. In these systems, the free energies associated with individual intermolecular interactions can therefore be reliably treated in an additive fashion. The results suggest that alternative explanations should be considered for cooperative phenomena observed in other systems, and a rationale based on the population of partially bound states in flexible molecules is proposed to account for the enthalpic chelate effect and enthalpy-entropy compensation.     

Modelling of carbohydrate–aromatic interactions: <Emphasis Type="Italic">ab initio </Emphasis>energetics and force field performance

Summary Aromatic amino acid residues are often present in carbohydrate-binding sites of proteins. These binding sites are characterized by a placement of a carbohydrate moiety in a stacking orientation to an aromatic ring. This arrangement is an example of CH/π interactions. Ab initio interaction energies for 20 carbohydrate–aromatic complexes taken from 6 selected ultra-high resolution X-ray structures of glycosidases and carbohydrate-binding proteins were calculated. All interaction energies of a pyranose moiety with a side chain of an aromatic residue were calculated as attractive with interaction energy ranging from −2.8 to −12.3 kcal/mol as calculated at the MP2/6-311+G(d) level. Strong attractive interactions were observed for a wide range of orientations of carbohydrate and aromatic ring as present in selected X-ray structures. The most attractive interaction was associated with apparent combination of CH/π interactions and classical H-bonds. The failure of Hartree–Fock method (interaction energies from +1.0 to −6.9 kcal/mol) can be explained by a dispersion nature of a majority of the studied complexes. We also present a comparison of interaction energies calculated at the MP2 level with those calculated using molecular mechanics force fields (OPLS, GROMOS, CSFF/CHARMM, CHEAT/CHARMM, Glycam/AMBER, MM2 and MM3). For a majority of force fields there was a strong correlation with MP2 values. RMSD between MP2 and force field values were 1.0 for CSFF/CHARMM, 1.2 for Glycam/AMBER, 1.2 for GROMOS, 1.3 for MM3, 1.4 for MM2, 1.5 for OPLS and to 2.3 for CHEAT/CHARMM (in kcal/mol). These results show that molecular mechanics approximates interaction energies very well and support an application of molecular mechanics methods in the area of glycochemistry and glycobiology.     

A quantum chemistry study of binding carotenoids in the bacterial light-harvesting complexes

Carotenoids play the dual function of light harvesting and photoprotection in photosynthetic organisms. Despite their functional importance, the molecular basis for binding of carotenoids in the photosynthetic proteins is poorly understood. We have discovered that all carotenoids are surrounded either by aromatic residues or by chlorophylls in all known crystal structures of the photosynthetic pigment-protein complexes. The intermolecular pi-pi stacking interactions between carotenoids and the surrounding aromatic residues in the light-harvesting complex II (LH-II) of Rhodospirillum molischianum were analyzed by high level ab initio electronic structure calculations. Intermolecular interaction energies were calculated with the second-order M?ller-Plesset perturbation method (MP2) using the modified 6-31G*(0.25) basis set with diffuse d-polarization by Hobza and co-workers. The MP2/6-31G*(0.25) calculations yield a total stabilization energy of -15.66 kcal/mol between the carotenoid molecule and the four surrounding aromatic residues (alpha-Trp-23, beta-Phe-20, beta-Phe-24, beta-Phe-27). It is thus concluded that pi-pi stacking interactions between carotenoids and the aromatic residues play an essential role in binding carotenoids in the LH-II complex of Rhodospirillum molischianum. The physical nature of the pi-pi stacking interactions was further analyzed, and the dispersion interactions were found to be the dominant intermolecular attraction force. There is also a substantial electrostatic contribution to the overall intermolecular stabilization energy.     

Metal ligand aromatic cation-pi interactions in metalloproteins: ligands coordinated to metal interact with aromatic residues

Cation-pi interactions between aromatic residues and cationic amino groups in side chains and have been recognized as noncovalent bonding interactions relevant for molecular recognition and for stabilization and definition of the native structure of proteins. We propose a novel type of cation-pi interaction in metalloproteins; namely interaction between ligands coordinated to a metal cation--which gain positive charge from the metal--and aromatic groups in amino acid side chains. Investigation of crystal structures of metalloproteins in the Protein Data Bank (PDB) has revealed that there exist quite a number of metalloproteins in which aromatic rings of phenylalanine, tyrosine, and tryptophan are situated close to a metal center interacting with coordinated ligands. Among these ligands are amino acids such as asparagine, aspartate, glutamate, histidine, and threonine, but also water and substrates like ethanol. These interactions play a role in the stability and conformation of metalloproteins, and in some cases may also be directly involved in the mechanism of enzymatic reactions, which occur at the metal center. For the enzyme superoxide dismutase, we used quantum chemical computation to calculate that Trp163 has an interaction energy of 10.09 kcal mol(-1) with the ligands coordinated to iron.     

Polar group enhanced gas-phase acidities of carboxylic acids: an investigation of intramolecular electrostatic interaction

We studied the effects of polar groups on the gas-phase acidities of carboxylic acids experimentally and computationally. In this connection, the gas-phase acidities (DeltaH(acid), the enthalpy of deprotonation, and DeltaG(acid), the deprotonation free energy) of borane-complexed methylaminoacetic acid ((CH(3))2N(BH(3))CH(2)CO(2)H) and methylthioacetic acid (CH(3)S(BH(3))CH(2)CO(2)H) were measured using the kinetic method in a flowing afterglow-triple quadrupole mass spectrometer. The values of DeltaH(acid) and DeltaG(acid) of (CH(3))2N(BH(3))CH(2)CO(2)H were determined to be 328.8 +/- 1.9 and 322.1 +/- 1.9 kcal/mol, and those of CH(3)S(BH(3))CH(2)CO(2)H were determined to be 325.8 +/- 1.9 and 319.2 +/- 1.9 kcal/mol, respectively. The theoretical enthalpies of deprotonation of (CH(3))2N(BH(3))CH(2)CO(2)H (329.2 kcal/mol) and CH(3)S(BH(3))CH(2)CO(2)H (325.5 kcal/mol) were calculated at the B3LYP/6-31+G(d) level of theory. The calculated enthalpies of deprotonation of N-oxide-acetic acid (CH(3)NOCH(2)CO(2)H, 329.4 kcal/mol) and S-oxide-acetic acid (CH(3)SOCH(2)CO(2)H, 328.6 kcal/mol) are comparable to the experimental results for borane-complexed methylamino- and methylthioacetic acids. The enthalpy of deprotonation of sulfone-acetic acid (CH(3)SO2CH(2)CO(2)H, 326.1 kcal/mol) is about 2 kcal/mol lower than the S-oxide-acetic acid. The calculated enthalpy of deprotonation of sulfoniumacetic acid, (CH(3))2S+CH(2)CO(2)H, is 243.0 kcal/mol. Compared to the corresponding reference molecules, CH(3)NHCH(2)CO(2)H and CH(3)SCH(2)CO(2)H, the dipolar group and the monopolar group substituted carboxylic acids are stronger acids by 11-14 and 97 kcal/mol, respectively. We correlated the changes of the acidity upon a polar group substitution to the electrostatic free energy within the carboxylate anion. The acidity enhancements in polar group substituted carboxylic acids are the results of the favorable electrostatic interactions between the polar group and the developing charge at the carboxyl group.