The value of the proton hydration free energy, ΔGhyd(H+), has been quoted in the literature to be from −252.6 to −262.5 kcal/mol. In this article, we present a theoretical model for calculating the hydration free energy of ions in aqueous solvent and use this model to calculate the proton hydration free energy, ΔGhyd(H+), in an effort to resolve the uncertainty concerning its exact value. In the model we define ΔGhyd(H+) as the free energy change associated with the following process: ΔG[H+(gas)+H2nOn(aq)→H+(H2nOn)(aq)], where the solvent is represented by a neutral n-water cluster embedded in a dielectric continuum and the solvated proton is represented by a protonated n-water cluster also in the continuum. All solvated species are treated as quantum mechanical solutes coupled to a dielectric continuum using a self consistent reaction field cycle. We investigated the behavior of ΔGhyd(H+) as the number of explicit waters of hydration is increased from n=1 to n=6. As n increases from 1 to 3, the hydration free energy decreases dramatically. However, for n=4–6 the hydration free energy maintains a relatively constant value of −262.23 kcal/mol. These results indicate that the first hydration shell of the proton is composed of at least four water molecules. The constant value of the hydration free energy for n⩾4 strongly suggests that the proton hydration free energy is at the far lower end of the range of values that have been proposed in the literature.

We calculate free energy changes of ionization reactions in aqueous solvent using a self-consistent reaction field method. In the calculations all species are treated as quantum mechanical solutes coupled to a solvent dielectric continuum. We show for a series of substituted imidazole compounds that both absolute and relative pK a values for the deprotonation of nitrogen on the imidazole ring can be obtained with an average absolute deviation of 0.8 units from experiment. This degree of accuracy is possible only if the solutes are treated at the correlated level using either G2 type or density functional theory. Inconsistencies in published experimental free energies of hydration that might undermine the reliability of the calculated absolute pK a values are discussed.

The mechanism of the hydrolysis reaction of guanosine triphosphate (GTP) by the protein complex Ras-GAP (p21(ras) - p120(GAP)) has been modeled by the quantum mechanical-molecular mechanical (QM/MM) and ab initio quantum calculations. Initial geometry configurations have been prompted by atomic coordinates of a structural analog (PDBID:1WQ1). It is shown that the minimum energy reaction path is consistent with an assumption of two-step chemical transformations. At the first stage, a unified motion of Arg789 of GAP, Gln61, Thr35 of Ras, and the lytic water molecule results in a substantial spatial separation of the gamma-phosphate group of GTP from the rest of the molecule (GDP). This phase of hydrolysis process proceeds through the low-barrier transition state TS1. At the second stage, Gln61 abstracts and releases protons within the subsystem including Gln61, the lytic water molecule and the gamma-phosphate group of GTP through the corresponding transition state TS2. Direct quantum calculations show that, in this particular environment, the reaction GTP + H(2)O --> GDP + H(2)PO(4) (-) can proceed with reasonable activation barriers of less than 15 kcal/mol at every stage. This conclusion leads to a better understanding of the anticatalytic effect of cancer-causing mutations of Ras, which has been debated in recent years.

The energy of dimerization of two N-methylacetamide (NMA) molecules in vacuum is calculated using density functional theory. Natural orbital analysis suggests that the dimerization energy of -6.6 kcal/mol is predominantly due to the (NsH‚‚‚OdC) donor-acceptor interaction. The gas phase to water hydration free energies and the free energies of transfer from the aqueous phase to liquid alkane of hydrogen bonded, (NsH‚‚‚OdC), and nonbonded, (NsH,OdC), groups are calculated using a continuum solvent model. On the basis of these calculations, we estimate the free energy of forming an amide hydrogen bond in the context of the NMA dimer in water and in liquid alkane as ∼-1 and ∼-5 kcal/mol, respectively. The relevance of these calculations to processes such as protein folding and membrane insertion of proteins is discussed.

We have used an 'activated' molecular dynamics approach to simulate flap opening in HIV-1 protease. An initial impulse for flap opening was provided by applying harmonic constraints to non-flap residues. After an initial 'melting' phase, the two beta-hairpin structures that constitute the flaps opened to a 25 A gap within 200 ps of simulation. Analysis of backbone torsion angles suggests that flap opening is related to conformational changes at Lys 45, Met 46, Gly 52 and Phe 53. In contrast, similar molecular dynamics simulations on the M46I mutant, which is associated with drug resistance, indicates that this mutation stabilizes the flaps in a closed conformation.

The hydrolysis reaction of guanosine triphosphate (GTP) by p21(ras) (Ras) has been modeled by using the ab initio type quantum mechanical-molecular mechanical simulations. Initial geometry configurations have been prompted by atomic coordinates of the crystal structure (PDBID: 1QRA) corresponding to the prehydrolysis state of Ras in complex with GTP. Multiple searches of minimum energy geometry configurations consistent with the hydrogen bond networks have been performed, resulting in a series of stationary points on the potential energy surface for reaction intermediates and transition states. It is shown that the minimum energy reaction path is consistent with an assumption of a two-step mechanism of GTP hydrolysis. At the first stage, a unified action of the nearest residues of Ras and the nearest water molecules results in a substantial spatial separation of the gamma-phosphate group of GTP from the rest of the molecule (GDP). This phase of hydrolysis process proceeds through the low barrier (16.7 kcal/mol) transition state TS1. At the second stage, the inorganic phosphate is formed in consequence of proton transfers mediated by two water molecules and assisted by the Gln61 residue from Ras. The highest transition state at this segment, TS3, is estimated to have an energy 7.5 kcal/mol above the enzyme-substrate complex. The results of simulations are compared to the previous findings for the GTP hydrolysis in the Ras-GAP (p21(ras)-p120(GAP)) protein complex. Conclusions of the modeling lead to a better understanding of the anticatalytic effect of cancer causing mutation of Gln61 from Ras, which has been debated in recent years.

The (salen)Mn(III)-catalyzed epoxidation reaction mechanism has been investigated using density functional theory (DFT). There is considerable interest in and controversy over the mechanism of this reaction. The results of experimental studies have offered some support for three different reaction mechanisms: concerted, stepwise radical, and metallooxetane mediated. In this paper, a theoretical examination of the reaction suggests a novel mechanism that describes the reaction as a multichannel process combining both concerted and stepwise radical pathways. The competing channels have different spin states: the singlet, the triplet, and the quintet. The singlet reaction pathway corresponds to a concerted mechanism and leads exclusively to a cis epoxide product. In contrast, the triplet and quintet reactions follow a stepwise mechanism and lead to a product mixture of cis and trans epoxides. We show that the experimentally observed dependence of isomer product ratios on electronic effects connected with the substitution of the catalyst ligands is due to changing the relative position and, hence, the relative activities of the channels with different cis-trans yields. Because the results and conclusions of the present work dramatically differ from the results and conclusion of the recent DFT theoretical investigation (Linde, C.; Akermark, B; Norrby, P.-O.; Svensson, M. J. Am. Chem. Soc. 1999, 121, 5083.), we studied possible sources for the deep contradictions between the two works. The choice of the DFT functional and a model has been shown to be crucial for accurate results. Using high level ab initio calculations (coupled cluster-CCSD(T)), we show that the computational procedure employed in this study generates significantly more reliable numerical results. It is also shown that the smaller cationic model without a chlorine ligand that was used by Linde et al. is too oversimplified with respect to our larger neutral model. For this reason, using the cationic model led to a qualitatively wrong quintet reaction profile that played a key role in theoretical postulates in the earlier work.

The hydration free energies relative to that of the proton are calculated for a representative set of monatomic ions Z±. These include cationic forms of the alkali earth elements Li, Na, and K, and anionic forms of the halogens F, Cl, and Br. In the current model the relative ion hydration free energy is defined as Δ[ΔGhyd(Z±)]=G(Z±[H2O]n(aq))−G(H+[H2O]n(aq))−G(Z±(gas))−G(H+(gas)), where the solvated ions are represented by ion–water clusters coupled to a dielectric continuum using a self-consistent reaction field cycle. An investigation of the behavior of Δ[ΔGhyd(Z±)] as the number of explicit waters of hydration is increased reveals convergence by n=4. This convergence indicates that the free energy change for the addition of water to a solvated proton–water complex is the same as the free energy change associated with the addition of water to a solvated Z±–water complex. This is true as long as there are four explicitly solvating waters associated with the ion. This convergence is independent of the type of monatomic ion studied and it occurs before the first hydration shell of the ions (typically ⩾6) is satisfied. Structural analysis of the ion–water clusters reveals that the waters within the cluster are more likely to form hydrogen bonds with themselves when clustering around anions than when clustering around cations. This suggests that for small ion–water clusters, anions are more likely to be externally solvated than cations.