Predicting the energy of the water exchange reaction and free energy of solvation for the uranyl ion in aqueous solution

The structures and vibrational frequencies of UO2(H2O)4(2+) and UO2(H2O)5(2+) have been calculated using density functional theory and are in reasonable agreement with experiment. The energies of various reactions were calculated at the density functional theory (DFT) and MP2 levels; the latter provides the best results. Self-consistent reaction field calculations in the PCM and SCIPCM approximations predicted the free energy of the water exchange reaction, UO2(H2O)4(2+) + H2O <--> UO2(H2O)5(2+). The calculated free energies of reaction are very sensitive to the choice of radii (O and H) and isodensity values in the PCM and SCIPCM models, respectively. Results consistent with the experimental HEXS value of -1.19 +/- 0.42 kcal/mol (within 1-3 kcal/mol) are obtained with small cavities. The structures and vibrational frequencies of the clusters with second solvation shell waters: UO2(H2O)4(H2O)8(2+), UO2(H2O)4(H2O)10(2+), UO2(H2O)4(H2O)11(2+), UO2(H2O)5(H2O)7(2+), and UO2(H2O)5(H2O)10(2+), were calculated and are in better agreement with experiment as compared to reactions involving only UO2(H2O)4(2+) and UO2(H2O)5(2+). The MP2 reaction energies for water exchange gave gas-phase results that agreed with experiment in the range -5.5 to +3.3 kcal/mol. The results were improved by inclusion of a standard PCM model with differences of -1.2 to +2.7 kcal/mol. Rearrangement reactions based on an intramolecular isomerization leading to a redistribution of water in the two shells provide good values in comparison to experiment with values of Delta G(exchange) from -2.2 to -0.5 kcal/mol so the inclusion of a second hydration sphere accounts for most solvation effects. Calculation of the free energy of solvation of the uranyl cation yielded an upper bound to the solvation energy of -410 +/- 5 kcal/mol, consistent with the best experimental value of -421 +/- 15 kcal/mol.     

On the determination of monomer dissociation energies of small water clusters from photoionization experiments

Recently, monomer dissociation energies of neutral water clusters were estimated via a thermodynamic cycle that utilized the measured appearance energies of vacuum ultraviolet (VUV) photoionized water clusters and the previously reported dissociation energies of protonated water clusters. The thermodynamic cycle incorrectly assumed that the energy difference between the (H2O)n+ --> (H2O)n-1H+ + OH* asymptotes (the relaxation energy) was zero. We show that these relaxation energies are large and cannot be neglected in the analysis. Thus, the neutral water cluster monomer dissociation energies cannot be directly determined from the measured ionization potentials because they are themselves involved in the appropriate thermodynamic cycle.     

Solvation of monovalent anions in acetonitrile and N,N-dimethylformamide: Parameterization of the IEF-PCM model

The present work reports the parameterization of the polarizable continuum model for predicting the free energies of solvation for monovalent anions in acetonitrile and N,N-dimethylformamide. The parameterization of the model for acetonitrile employed the experimental free energies of solvation for a set of 12 charged solutes, containing H, C, N, O, S, F, Cl, Br, and I atoms. For the N,N-dimethylformamide solutions, experimental solvation free energies for 11 monovalent anions were used. A mean absolute error of 0.7 kcal/mol in the solvation free energies has been achieved for the 12 anions in acetonitrile, whereas the mean absolute error for the 11 anions corresponds to 0.5 kcal/mol in N,N-dimethylformamide. These results indicate that the polarizable continuum model is a suitable methodology for the study of thermodynamic effects in solutions of monovalent anions in both solvents.     

Accurate pK(a) calculations for carboxylic acids using complete basis set and Gaussian-n models combined with CPCM continuum solvation methods

Complete Basis Set and Gaussian-n methods were combined with CPCM continuum solvation methods to calculate pK(a) values for six carboxylic acids. An experimental value of -264.61 kcal/mol for the free energy of solvation of H(+), DeltaG(s)(H(+)), was combined with a value for G(gas)(H(+)) of -6.28 kcal/mol to calculate pK(a) values with Cycle 1. The Complete Basis Set gas-phase methods used to calculate gas-phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed-phase free energies are slightly less accurate than the gas-phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. The use of Cycle 1 and the Complete Basis Set models combined with the CPCM solvation methods yielded pK(a) values accurate to less than half a pK(a) unit.     

Single-ion solvation free energies and the normal hydrogen electrode potential in methanol, acetonitrile, and dimethyl sulfoxide

The division of thermodynamic solvation free energies of electrolytes into contributions from individual ionic constituents is conventionally accomplished by using the single-ion solvation free energy of one reference ion, conventionally the proton, to set the single-ion scales. Thus, the determination of the free energy of solvation of the proton in various solvents is a fundamental issue of central importance in solution chemistry. In the present article, relative solvation free energies of ions and ion-solvent clusters in methanol, acetonitrile, and dimethyl sulfoxide (DMSO) have been determined using a combination of experimental and theoretical gas-phase free energies of formation, solution-phase reduction potentials and acid dissociation constants, and gas-phase clustering free energies. Applying the cluster pair approximation to differences between these relative solvation free energies leads to values of -263.5, -260.2, and -273.3 kcal/mol for the absolute solvation free energy of the proton in methanol, acetonitrile, and DMSO, respectively. The final absolute proton solvation free energies are used to assign absolute values for the normal hydrogen electrode potential and the solvation free energies of other single ions in the solvents mentioned above.     

Treating entropy and conformational changes in implicit solvent simulations of small molecules

Implicit solvent models are increasingly popular for estimating aqueous solvation (hydration) free energies in molecular simulations and other applications. In many cases, parameters for these models are derived to reproduce experimental values for small molecule hydration free energies. Often, these hydration free energies are computed for a single solute conformation, neglecting solute conformational changes upon solvation. Here, we incorporate these effects using alchemical free energy methods. We find significant errors when hydration free energies are estimated using only a single solute conformation, even for relatively small, simple, rigid solutes. For example, we find conformational entropy (TDeltaS) changes of up to 2.3 kcal/mol upon hydration. Interestingly, these changes in conformational entropy correlate poorly (R2 = 0.03) with the number of rotatable bonds. The present study illustrates that implicit solvent modeling can be improved by eliminating the approximation that solutes are rigid.     

Role of solvation in the reduction of the uranyl dication by water: a density functional study

We have studied the solvation of uranyl, UO(2)(2+), and the reduced species UO(OH)(2+) and U(OH)(2)(2+) systematically using three levels of approximation: direct application of a continuum model (M1); explicit quantum-chemical treatment of the first hydration sphere (M2); a combined quantum-chemical/continuum model approach (M3). We have optimized complexes with varying numbers of aquo ligands (n = 4-6) and compared their free energies of solvation. Models M1 and M2 have been found to recover the solvation energy only partially, underestimating it by approximately 100 kcal/mol or more. With our best model M3, the calculated hydration free energy Delta(h)G degrees of UO(2)(2+) is about -420 kcal/mol, which shifts to about -370 kcal/mol when corrected for the expected error of the model. This value agrees well with the experimentally determined interval, -437 kcal/mol < Delta(h)G degrees < -318 kcal/mol. Complexes with 5 and 6 aquo ligands have been found to be about equally favored with models M2 and M3. The same solvation models have been applied to a two-step reduction of UO(2)(2+) by water, previously theoretically studied in the gas phase. Our results show that the solvation contribution to the reaction free energy, about 60 kcal/mol, dominates the endoergicity of the reduction.     

Aqueous solvation free energies of ions and ion-water clusters based on an accurate value for the absolute aqueous solvation free energy of the proton

Thermochemical cycles that involve pKa, gas-phase acidities, aqueous solvation free energies of neutral species, and gas-phase clustering free energies have been used with the cluster pair approximation to determine the absolute aqueous solvation free energy of the proton. The best value obtained in this work is in good agreement with the value reported by Tissandier et al. (Tissandier, M. D.; Cowen, K. A.; Feng, W. Y.; Gundlach, E.; Cohen, M. J.; Earhart, A. D.; Coe, J. V. J. Phys. Chem. A 1998, 102, 7787), who applied the cluster pair approximation to a less diverse and smaller data set of ions. We agree with previous workers who advocated the value of -265.9 kcal/mol for the absolute aqueous solvation free energy of the proton. Considering the uncertainties associated with the experimental gas-phase free energies of ions that are required to use the cluster pair approximation as well as analyses of various subsets of data, we estimate an uncertainty for the absolute aqueous solvation free energy of the proton of no less than 2 kcal/mol. Using a value of -265.9 kcal/mol for the absolute aqueous solvation free energy of the proton, we expand and update our previous compilation of absolute aqueous solvation free energies; this new data set contains conventional and absolute aqueous solvation free energies for 121 unclustered ions (not including the proton) and 147 conventional and absolute aqueous solvation free energies for 51 clustered ions containing from 1 to 6 water molecules. When tested against the same set of ions that was recently used to develop the SM6 continuum solvation model, SM6 retains its previously determined high accuracy; indeed, in most cases the mean unsigned error improves when it is tested against the more accurate reference data.     

Predicting binding affinities of host-guest systems in the SAMPL3 blind challenge: the performance of relative free energy calculations

Relative free energy calculations based on molecular dynamics simulations are combined with available experimental binding free energies to predict unknown binding affinities of acyclic Cucurbituril complexes in the blind SAMPL3 competition. The predictions yield root mean square errors between 2.6 and 3.2 kcal/mol for seven host-guest systems. Those deviations are comparable to results for solvation free energies of small organic molecules. However, the standard deviations found in our simulations range from 0.4 to 2.4 kcal/mol, which indicates the need for better sampling. Three different approaches are compared. Bennett's Acceptance Ratio Method and thermodynamic integration based on the trapezoidal rule with 12 λ-points exhibit a root mean square error of 2.6 kcal/mol, while thermodynamic integration with Simpson's rule and 11 λ-points leads to a root mean square error of 3.2 kcal/mol. In terms of absolute median errors, Bennett's Acceptance Ratio Method performs better than thermodynamic integration with the trapezoidal rule (1.7 vs. 2.9 kcal/mol). Simulations of the deprotonated forms of the guest molecules exhibit a poorer correspondence to experimental results with a root mean square error of 5.2 kcal/mol. In addition, a decrease of the buffer concentration by approximately 20 mM in the simulations raises the root mean square error to 3.8 kcal/mol.     

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.     

Accuracy of free energies of hydration for organic molecules from 6-31g*-derived partial charges

Absolute free energies of hydration have been computed for 13 diverse organic molecules using partial charges derived from ab initio 6-31G* wave functions. Both Mulliken charges and charges fit to the electrostatic potential surface (EPS) were considered in conjunction with OPLS Lennard–Jones parameters for the organic molecules and the TIP4P model of water. Monte Carlo simulations with statistical perturbation theory yielded relative free energies of hydration. These were converted to absolute quantities through perturbations to reference molecules for which absolute free energies of hydration had been obtained previously in TIP4P water. The average errors in the computed absolute free energies of hydration are 1.1 kcal/mol for the 6-31G* EPS charges and 4.0 kcal/mol for the Mulliken charges. For the EPS charges, the largest individual errors are under 2 kcal/mol except for acetamide, in which case the error is 3.7 kcal/mol. The hydrogen bonding between the organic solutes and water has also been characterized. © John Wiley & Sons, Inc.     

Experimentation with different thermodynamic cycles used for pKa calculations on carboxylic acids using complete basis set and Gaussian‐n models combined with CPCM continuum solvation methods

Complete basis set and Gaussian‐n methods were combined with Barone and Cossi's implementation of the polarizable conductor model (CPCM) continuum solvation methods to calculate pK_a values for six carboxylic acids. Four different thermodynamic cycles were considered in this work. An experimental value of ?264.61 kcal/mol for the free energy of solvation of H⁺, ΔG_s(H⁺), was combined with a value for G_gas(H⁺) of ?6.28 kcal/mol, to calculate pK_a values with cycle 1. The complete basis set gas‐phase methods used to calculate gas‐phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed‐phase free energies are slightly less accurate than the gas‐phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. Thermodynamic cycles that include an explicit water in the cycle are not accurate when the free energy of solvation of a water molecule is used, but appear to become accurate when the experimental free energy of vaporization of water is used. This apparent improvement is an artifact of the standard state used in the calculation. Geometry relaxation in solution does not improve the results when using these later cycles. The use of cycle 1 and the complete basis set models combined with the CPCM solvation methods yielded pK_a values accurate to less than half a pK_a unit. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Electronic structure and solvation of copper and silver ions: a theoretical picture of a model aqueous redox reaction

Electronic states and solvation of Cu and Ag aqua ions are investigated by comparing the Cu(2+) + e(-)--> Cu(+) and Ag(2+) + e(-)--> Ag(+) redox reactions using density functional-based computational methods. The coordination number of aqueous Cu(2+) is found to fluctuate between 5 and 6 and reduces to 2 for Cu(+), which forms a tightly bound linear dihydrate. Reduction of Ag(2+) changes the coordination number from 5 to 4. The energetics of the oxidation reactions is analyzed by comparing vertical ionization potentials, relaxation energies, and vertical electron affinities. The model is validated by a computation of the free energy of the full redox reaction Ag(2+) + Cu(+) --> Ag(+) + Cu(2+). Investigation of the one-electron states shows that the redox active frontier orbitals are confined to the energy gap between occupied and empty states of the pure solvent and localized on the metal ion hydration complex. The effect of solvent fluctuations on the electronic states is highlighted in a computation of the UV absorption spectrum of Cu(+) and Ag(+).     

Calculation of the solvation free energy of neutral and ionic molecules in diverse solvents

The solvation free energy density (SFED) model was modified to extend its applicability and predictability. The parametrization process was performed with a large, diverse set of solvation free energies that included highly polar and ionic molecules. The mean absolute error for 1200 solvation free energies of the 379 neutral molecules in 9 organic solvents and water was 0.40 kcal/mol, and for 90 hydration free energies of ions was 1.7 kcal/mol. Overall, the calculated solvation free energies of a wide range of solute functional groups in diverse solvents were consistent with experimental data.     

Free energy of solvation of simple ions: molecular-dynamics study of solvation of Cl- and Na+ in the ice/water interface

Molecular-dynamics simulations of Cl(-) and Na(+) ions are performed to calculate ionic solvation free energies in both bulk simple point-charge/extended water and ice 1 h at several different temperatures, and at the basal ice 1 h/water interface. For the interface we calculate the free energy of "transfer" of the ions across the ice/water interface. For the ions in bulk water in the NPT ensemble at 298 K and 1 atm, results are found to be in good agreement with experiments, and with other simulation results. Simulations performed in the NVT ensemble are shown to give equivalent solvation free energies, and this ensemble is used for the interfacial simulations. Solvation free energies of Cl(-) and Na(+) ions in ice at 150 K are found to be approximately 30 and approximately 20 kcal mol(-1), respectively, less favorable than for water at room temperature. Near the melting point of the model the solvation of the ions in water is the same (within statistical error) as that measured at room temperature, and in the ice is equivalent and approximately 10 kcal mol(-1) less favorable than the liquid. The free energy of transfer for each ion across ice/water interface is calculated and is in good agreement with the bulk observations for the Cl(-) ion. However, for the model of Na(+) the long-range electrostatic contribution to the free energy was more negative in the ice than the liquid, in contrast with the results observed in the bulk calculations.     

Calculation of the water-cyclohexane transfer free energies of neutral amino acid side-chain analogs using the OPLS all-atom force field

We calculated the free energy of solvation of the neutral analogs of 18 amino acid side-chains (not including glycine and proline) using the OPLS all-atom force field in TIP4P water, SPC water, and cyclohexane by molecular dynamics simulation and thermodynamic integration. The average unsigned errors in the free energies of solvation in TIP4P, SPC, and cyclohexane are 4.4, 4.9, and 2.1 kJ/mol respectively. Most of the calculated hydration free energies are not favorable enough compared to experiment. The largest errors are found for tryptophan, histidine, glutamic acid, and glutamine. The average unsigned errors in the free energy of transfer from TIP4P to cyclohexane and from SPC to cyclohexane are 4.0 and 4.1 kJ/mol, respectively. The largest errors, of more than 7.5 kJ/mol, are found for histidine, glutamine, and glutamatic acid.     

Theoretical study of [Ni (H2O)n]2+(H2O)m (n ≤ 6, m ≤ 18)

The [Ni-(H(2)O)(n)](2+)(H(2)O)(m) (n ≤ 6, m ≤ 18) complexes were studied by means of first-principles all-electron calculations performed with the BPW91 gradient corrected functional and the 6-311+G(d,p) basis sets for the H, O, and Ni atoms. Triplet states were found as low-lying states for each (n, m) combination. The estimated Ni(2+)-(H(2)O)(n) binding energies (112.8-57.4 kcal/mol for the first layer and 52.0-23.0 kcal/mol for the second one) decreases and the Ni(2+)-OH(2) bond lengths lengthen as n + m increases. With six H(2)O moieties the Ni(2+) ion furnishes its first coordination sphere of octahedral geometry. Further water addition renders the formation of the second layer. The effect of Ni(2+) on the (H(2)O)(n)···(H(2)O)(m) hydrogen bond formation for several "n" and "m" combinations was studied, revealing an enhancement of this kind of bonding, which is of key importance for the stabilization and growth of the clusters. For some n + m isomers the second layer appears before the first octahedral layer is fully formed. For example, the square planar Ni(2+)-(H(2)O)(4) core originates two-dimensional 4 + 2 and 4 + 4 isomers, where each outer water molecule accepts two H-bonds, lying 2.0 kcal/mol above the 6 and 6 + 2 ground states. The clusters were also studied by IR spectra; the OH stretching vibrational frequencies allowed us to identify the outer solvation shells by the presence of red-shifted hydrogen bond regions.     

On the viability of small endohedral hydrocarbon cage complexes: X@C4H4, X@C8H8, X@C8H14, X@C10H16, X@C12H12, AND X@C16H16

Small hydrocarbon complexes (X@cage) incorporating cage-centered endohedral atoms and ions (X = H(+), H, He, Ne, Ar, Li(0,+), Be(0,+,2+), Na(0,+), Mg(0,+,2+)) have been studied at the B3LYP/6-31G(d) hybrid HF/DFT level of theory. No tetrahedrane (C(4)H(4), T(d)()) endohedral complexes are minima, not even with the very small hydrogen atom or beryllium dication. Cubane (C(8)H(8), O(h)()) and bicyclo[2.2.2]octane (C(8)H(14), D(3)(h)()) minima are limited to encapsulating species smaller than Ne and Na(+). Despite its intermediate size, adamantane (C(10)H(16), T(d)()) can enclose a wide variety of endohedral atoms and ions including H, He, Ne, Li(0,+), Be(0,+,2+), Na(0,+), and Mg(2+). In contrast, the truncated tetrahedrane (C(12)H(12), T(d)()) encapsulates fewer species, while the D(4)(d)() symmetric C(16)H(16) hydrocarbon cage (see Table of Contents graphic) encapsulates all but the larger Be, Mg, and Mg(+) species. The host cages have more compact geometries when metal atoms, rather than cations, are inside. This is due to electron donation from the endohedral metals into C-C bonding and C-H antibonding cage molecular orbitals. The relative stabilities of endohedral minima are evaluated by comparing their energies (E(endo)) to the sum of their isolated components (E(inc) = E(endo) - E(cage) - E(x)) and to their exohedral isomer energies (E(isom) = E(endo) - E(exo)). Although exohedral binding is preferred to endohedral encapsulation without exception (i.e., E(isom) is always exothermic), Be(2+)@C(10)H(16) (T(d)(); -235.5 kcal/mol), Li(+)@C(12)H(12) (T(d)(); 50.2 kcal/mol), Be(2+)@C(12)H(12) (T(d)(); -181.2 kcal/mol), Mg(2+)@C(12)H(12) (T(d)(); -45.0 kcal/mol), Li(+)@C(16)H(16) (D(4)(d)(); 13.3 kcal/mol), Be(+)@C(16)H(16) (C(4)(v)(); 31.8 kcal/mol), Be(2+)@C(16)H(16) (D(4)(d)(); -239.2 kcal/mol), and Mg(2+)@C(16)H(16) (D(4)(d)(); -37.7 kcal/mol) are relatively stable as compared to experimentally known He@C(20)H(20) (I(h)()), which has an E(inc) = 37.9 kcal/mol and E(isom) = -35.4 kcal/mol. Overall, endohedral cage complexes with low parent cage strain energies, large cage internal cavity volumes, and a small, highly charged guest species are the most viable synthetic targets.     

Extension of the platform of applicability of the SM5.42R universal solvation model

We present eight new parameterizations of the SM5.42R solvation model: in particular we present parameterizations for HF/MIDI!, HF/6-31G^*, HF/6-31+G^*, HF/cc-pVDZ, AM1, PM3, BPW91/MIDI!, and B3LYP/MIDI!. Two of the new cases are parameterized using the reaction-field operator presented previously, and six of the new cases are parameterized with a simplified reaction-field operator; results obtained by the two methods are compared for selected examples. For a training set of 2135 data for 275 neutral solutes containing H, C, N, O, F, S, P, Cl, Br, and I in 91 solvents (water and 90 nonaqueous solvents), seven of the eight new parameterizations give mean unsigned errors in the range 0.43–0.46 kcal/mol, and the eighth – for a basis set containing diffuse functions – gives a mean unsigned error of 0.53 kcal/mol. The mean unsigned error for 49 ionic solutes (containing the same elements) in water is 3.5–3.9 kcal/mol for the Hartree–Fock, Becke–Perdew–Wang-1991 and Becke three-parameter Lee–Yang–Parr cases and 4.1 and 4.0 kcal/mol for parameterized model 3 and Austin model 1, respectively. The methods are tested for sensitivity of solvation free energies to geometry and for predicting partition coefficients of carbonates, which were not included in the training set. Received: 24 November 1998 / Accepted: 31 December 1998 / Published online: 7 June 1999     

Is the peptide bond formation activated by Cu(2+) interactions? Insights from density functional calculations

The catalytic role that Cu(2+) cations play in the peptide bond formation has been addressed by means of density functional calculations. First, the Cu(2+)-(glycine)2 --> Cu(2+)-(glycylglycine) + H2O reaction was investigated since mass spectrometry low collision activated dissociation (CAD) spectra of Cu(2+)-(glycine)2 led to the elimination of a water molecule, which suggested that an intracomplex peptide bond formation might have occurred. Results show that this intracomplex condensation is associated to a very high free energy barrier (97 kcal mol(-1)) and reaction free energy (66 kcal mol(-1)) because of the loss of metal coordination during the reaction. Second, on the basis of the salt-induced peptide formation theory, the condensation reaction between two glycines was studied in aqueous solution using discrete water molecules and the conductor polarized continuum model (CPCM) continuous method. It is found that the synergy between the interaction of glycines with Cu(2+) and the presence of water molecules acting as proton-transfer helpers significantly lower the activation barrier (from 55 kcal/mol for the uncatalyzed system to 20 kcal/mol for the Cu(2+) solvated system) which largely favors the formation of the peptide bond.