Theoretical study of HOCl-catalyzed keto-enol tautomerization of ß-cyclopentanedione in an explicit water environment.

The mechanism of acid-catalyzed keto-enol tautomerization of ß-cyclopentanedione (CPD) in solution is studied computationally. Reaction profiles are first calculated for a limited solvation environment using ab initio and density functional methods. Barrier heights for systems including up to one hydration shell of explicit water molecules depend strongly on the number of waters involved in proton transfer and to a lesser but significant extent on the number of waters forming hydrogen bonds with waters in the proton-transfer chain (each such water reduces the barrier by 4.4 kcal/mol on average). Barriers of 8-13 kcal/mol were obtained when a full or nearly full hydration shell was present, consistent with calculations for nonacid-catalyzed keto-enol tautomerization of related molecules. The presence of HOCl reduced the barrier by 4.5 kcal/mol in relation to the gas phase, consistent with the well-known principle that keto-enol tautomerization can be acid- or base-catalyzed. The reaction was also modeled beginning with snapshots of reactant conformations taken from a 300 K molecular dynamics simulation of CPD, HOCl, and 324 explicit waters. Reaction profiles were calculated at a QM/MM level with waters in the first hydration shell either fixed or energy-minimized at each step along the reaction coordinate. A substantial variation in barrier height was observed in both cases, depending primarily on electrostatic interactions (hydrogen bonding) with first-hydration-shell waters and, to a lesser extent, on electrostatic interactions with more distant waters and geometric distortion effects. For the lowest barriers, the extent of barrier reduction by waters involved in proton transfer is consistent with the limited solvation results, but further barrier reduction due to hydrogen bonding to waters involved in proton transfer is not observed. It is postulated that this is because highly flexible structures such as extensive hydrogen bonding networks optimal for reaction are entropically disfavored and so may not contribute significantly to the observed reaction rate.

Enantiospecificity of chloroperoxidase-catalyzed epoxidation: biased molecular dynamics study of a cis-ß-methylstyrene/chloroperoxidase-compound I complex.

Molecular dynamics simulations of an explicitly solvated cis-ß-methylstyrene/chloroperoxidase-Compound I complex are performed to determine the cause of the high enantiospecificity of epoxidation. From the simulations, a two-dimensional free energy potential is calculated to distinguish binding potential wells from which reaction to 1S2R and 1R2S epoxide products may occur. Convergence of the free energy potential is accelerated with an adaptive biasing potential. Analysis of binding is followed by analysis of 1S2R and 1R2S reaction precursor structures in which the substrate, having left the binding wells, places its reactive double bond in steric proximity to the oxyferryl heme center. Structural analysis of binding and reaction precursor conformations is presented. We find that 1), a distortion of Glu(183) is important for CPO-catalyzed epoxidation as was postulated previously based on experimental results; 2), the free energy of binding does not provide significant differentiation between structures leading to the respective epoxide enantiomers; and 3), CPO's enantiospecificity toward cis-ß-methylstyrene is likely to be caused by a specific group of residues which form a hydrophobic core surrounding the oxyferryl heme center.

Correlation times and adiabatic barriers for methyl rotation in SNase.

The relation of rotational correlation times to adiabatic rotational barriers for alanine methyl groups in staphylococcal nuclease (SNase) is investigated. The hypothesis that methyl rotational barriers may be useful probes of local packing in proteins is supported by an analysis of ten X-ray crystal structures of SNase mutants. The barrier heights are consistent across a set of ten structures of a native SNase and mutants containing single-point mutations or single or double insertions, most in a ternary SNase complex. The barriers for different methyls have a range of 7.5 kcal/mol, which at 300 K would correspond to a five-order-of-magnitude range in correlation time. It is demonstrated that adiabatic rotational barriers can fluctuate significantly during an MD simulation of hydrated SNase, but that a Boltzmann weighted average is predictive of rotational correlation times determined from correlation functions. Even if a given methyl is on average quite sterically hindered, infrequently sampled low-barrier conformations may dominate the Boltzmann distribution. This result is consistent with the observed uniformity of NMR correlation times for (13)C-labeled methyls. The methyl barriers in simulation fluctuate on multiple time scales, which can make the precise relationship between methyl rotational correlation time and methyl rotation barriers complicated. The implications of these issues for the interpretation of correlation times determined from NMR and simulation are discussed.

Methyl dynamics in crystalline amino acids: MD and NMR.

Correlation times for rotation of deuterated methyls in crystalline leucine, valine, and cyclo-L-alanyl-L-alanine are calculated with molecular dynamics and compared with NMR data. The simulations distinguish between methyls having different steric environments in the crystal, yielding correlation times differing by a factor of up to 30 for methyls within a given crystal. MD and NMR correlation times agree to within a factor of 2. However, averaging over nonequivalent methyls can yield correlation functions that, although actually multiexponential, are well fit by single exponentials. This may have significance for interpreting NMR data; previous NMR data did not distinguish between the methyls in these crystals. Adiabatic rotational barriers calculated with the X-ray structure differ from effective barriers during simulation by up to +/-1 kcal/mol; the difference indicates that dynamical effects have a significant role in determining rotational correlation times.