Free energy modeling computational studies

The enzyme mechanism, however, remains elusive. Quantum mechanical models generally disfavor C6-protonation, but 02, 04, and C5-protonation mechanisms remain possibilities. Free energy computations also appear to indicate that C5-protonation is a feasible mechanism, as is direct decarboxylation without preprotonation O-protonation mechanisms have yet to be explored with these methods. Controversy remains, however, as to the roles of ground state destabilization, transition state stabilization, and dynamic effects. Because free energy models do take into account the entire enzyme active site, a comprehensive study of the relative energetics of pre-protonation and concerted protonation-decarboxylation at 02, 04, and C5 should be undertaken with such methods. In addition, quantum mechanical isotope effects are also likely to figure prominently in the ultimate identification of the operative ODCase mechanism. 

We have previously calculated conformational free energy differences for a well-suited model system, the catalytic subunit of cAMP-dependent protein kinase (cAPK), which is the best characterized member of the protein kinase family. It has been crystallized in three different conformations and our main focus was on how ligand binding shifts the equilibrium among these ([Helms and McCammon 1997]). As an example using state-of-the-art computational techniques, we summarize the main conclusions of this study and discuss a variety of methods that may be used to extend this study into the dynamic regime of protein domain motion. 

Solvents exert their influence on organic reactions through a complicated mixture of all possible types of noncovalent interactions. Chemists have tried to unravel this entanglement and, ideally, want to assess the relative importance of all interactions separately. In a typical approach, a property of a reaction (e.g. its rate or selectivity) is measured in a laige number of different solvents. All these solvents have unique characteristics, quantified by their physical properties (i.e. refractive index, dielectric constant) or empirical parameters (e.g. ET(30)-value, AN). Linear correlations between a reaction property and one or more of these solvent properties (Linear Free Energy Relationships - LFER) reveal which noncovalent interactions are of major importance. The major drawback of this approach lies in the fact that the solvent parameters are often not independent. Alternatively, theoretical models and computer simulations can provide valuable information. Both methods have been applied successfully in studies of the solvent effects on Diels-Alder reactions. 

The present chapter thus provides an overview of the current status of continuum models of solvation. We review available continuum models and computational techniques implementing such models for both electrostatic and non-electrostatic components of the free energy of solvation. We then consider a number of case studies, with particular focus on the prediction of heterocyclic tautomeric equilibria. In the discussion of the latter we center attention on the subtleties of actual chemical systems and some of the dangers of applying continuum models uncritically. We hope the reader will emerge with a balanced appreciation of the power and limitations of these methods. 

---