Empirical valence bond method EVB

For this reason, there has been much work on empirical potentials suitable for use on a wide range of systems. These take a sensible functional form with parameters fitted to reproduce available data. Many different potentials, known as molecular mechanics (MM) potentials, have been developed for ground-state organic and biochemical systems [58-60], They have the advantages of simplicity, and are transferable between systems, but do suffer firom inaccuracies and rigidity—no reactions are possible. Schemes have been developed to correct for these deficiencies. The empirical valence bond (EVB) method of Warshel [61,62], and the molecular mechanics-valence bond (MMVB) of Bemardi et al. [63,64] try to extend MM to include excited-state effects and reactions. The MMVB Hamiltonian is parameterized against CASSCF calculations, and is thus particularly suited to photochemistry. 

The approach presented above is referred to as the empirical valence bond (EVB) method (Ref. 6). This approach exploits the simple physical picture of the VB model which allows for a convenient representation of the diagonal matrix elements by classical force fields and convenient incorporation of realistic solvent models in the solute Hamiltonian. A key point about the EVB method is its unique calibration using well-defined experimental information. That is, after evaluating the free-energy surface with the initial parameter a , we can use conveniently the fact that the free energy of the proton transfer reaction is given by... 

The empirical valence bond (EVB) method of Warshel [19] has features of both the structurally and thermodynamically coupled QM/MM method. In the EVB method the different states of the process studied are described in terms of relevant covalent and ionic resonance structures. The potential energy surface of the QM system is calibrated to reproduce the known experimental... 

Valence bond (VB) theories or empirical valence bond (EVB) methods have been developed in order to solve this problem with bond potential functions that (i) allow the change of the valence bond network over time and (ii) are simple enough to be used efficiently in an otherwise classical MD simulation code. In an EVB scheme, the chemical bond in a dissociating molecule is described as the superposition of two states a less-polar bonded state and an ionic dissociated state. One of the descriptions is given by Walbran and Kornyshev in modeling of the water dissociation process.4,5 As... 

Warshel and coworkers have employed the empirical valence-bond (EVB) method [49] to simulate FERs for PT [50] and other reactions [51]. The PT step between two water molecules in the mechanism of the reaction catalysed by carbonic anhy-drase was described as an effective two-state problem involving reactant-like (H0H)(0H2) and product-like (HO )(HOH2+) VB structures [50a], Diabatic energy curves for these two VB structures were calibrated to reproduce the experimental free energy change for autodissociation in water, and the mixing of the... 

During our search for reliable methods for studies of enzymatic reactions it became apparent that, in studies of chemical reactions, it is more physical to calibrate surfaces that reflect bond properties (that is, valence bond-based, VB, surfaces) than to calibrate surfaces that reflect atomic properties (for example, MO-based surfaces). Furthermore, it appears to be very advantageous to force the potential surfaces to reproduce the experimental results of the broken fragments at infinite separation in solution. This can be effectively accomplished with the VB picture. The resulting empirical valence bond (EVB) method has been discussed extensively elsewhere [3, 24], but its main features will be outlined below, because it provides the most direct microscopic connection to PT processes. 

Enzyme catalyzed reactions can also be modeled using the empirical valence bond (EVB) method [21]. EVB methods are central in modeling enzyme reactions and catalysis. They are reviewed well elsewhere [22,12] and so are not the focus of this review. 

Warshel and co-workers > ° >i studied the dynamics of this Sn2 reaction using an interesting and different approach the empirical valence bond (EVB) method, which has been described in detail in a recent book by Warshel. The fundamental idea behind the application of the EVB method to this Sn2 reaction is that the reaction can be treated as a two state system, where the reactants and products are each taken to be separate quantum mechanical states with Hamiltonians and Hi- These states can be coupled together by an empirical coupling Hamiltonian so that when the two state Hamiltonian is diagonalized, the correct features of the ground state surface on which the reaction occurs are obtained. H, and H2 are taken by Warshel and coworkers to have analytic forms based on the gas phase parameters. 

One example of the use of linear response theory has been that of Hwang et al. in their studies of an reaction in solution. > o In their work, based on the empirical valence bond (EVB) method discussed earlier, they defined their reaction coordinate Q as the electrostatic contribution to the energy gap between the two valence bond states that are coupled together to create the potential energy surface on which the reaction occurs. Thus, the solvent coordinate is zero at the point where both valence states are solvated equivalently (i.e., at the transition state). Hwang et al. studied the time dependence of this coordinate through both molecular dynamics simulations and through a linear response treatment ... 

Warshel et improved and extended their QM/MM approach to the so-called empirical valence bond (EVB) method in which (in the last versions) the empirical parameters for an enzyme-catalyzed reaction are calibrated using ab initio calculations for this reaction in solutions. Using the EVB method, several catalytic mechanisms of enzyme activity were considered . 

These aspects are discussed in section 5.2, where we review different methodological approaches to describing reactions in enzymes as well as in solution. The emphasis is given to the empirical valence bond (EVB) method which is presented in some detail. In section 5.3, we examine the catalytic... 

More than 30 years ago Warshel proposed, on the basis of semiempirical simulations, an isomerization mechanism that could explain how this process can occur in the restricted space of the Rh binding pocket (Warshel 1976). Since two adjacent double bonds were found to isomerize simultaneously the mechanism reveal a so-called bicycle pedal motion. Due to the concerted rotation of two double bonds in opposite directions the overall conformational change is minimized and hence this mechanism was found to be space-saving. The empirical valence bond (EVB) method (Warshel and Levitt 1976) was used to compute the excited state potential energy surface of the chromophore during a trajectory calculation where the steric effects of the protein matrix were modeled by specific restraints on the retinal atoms. Since then, Warshel and his coworkers have improved the model employing better structural data and new computational developments (Warshel and Barboy 1982 Warshel and Chu 2001 Warshel et al. 1991). The main refinement of the bicycle pedal mechanism was that the simultaneous rotation of the adjacent double bonds is aborted at a twist of 40° and leads to the isomerization of only one bond (Warshel and Barboy 1982). 

First, we show how the FEG method realizes the FE changes through the process of stmctural optimization in solution at a finite temperature. As an example, we take the glycine zwitterion (ZW) in aqueous solution [6] which is the most simple but interesting amino acid. To describe the glycine ZW potential energy surface in aqueous solution, we employed the empirical valence bond (EVB) method which was prepared so as to reproduce a set of energies and forces calculated at the... 

---