The importance of reactant positioning in enzyme catalysis: A hybrid quantum mechanics/molecular mechanics study of a haloalkane dehalogenase

Hybrid quantum mechanics/molecular mechanics calculations using Austin Model 1 system-specific parameters were performed to study the SN2 displacement reaction of chloride from 1,2-dichloroethane (DCE) by nucleophilic attack of the carboxylate of acetate in the gas phase and by Asp-124 in the active site of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10. The activation barrier for nucleophilic attack of acetate on DCE depends greatly on the reactants having a geometry resembling that in the enzyme or an optimized gas-phase structure. It was found in the gas-phase calculations that the activation barrier is 9 kcal/mol lower when dihedral constraints are used to restrict the carboxylate nucleophile geometry to that in the enzyme relative to the geometries for the reactants without dihedral constraints. The calculated quantum mechanics/molecular mechanics activation barriers for the enzymatic reaction are 16.2 and 19.4 kcal/mol when the geometry of the reactants is in a near attack conformer from molecular dynamics and in a conformer similar to the crystal structure (DCE is gauche), respectively. This haloalkane dehalogenase lowers the activation barrier for dehalogenation of DCE by 2–4 kcal/mol relative to the single point energies of the enzyme's quantum mechanics atoms in the gas phase. SN2 displacements of this sort in water are infinitely slower than in the gas phase. The modest lowering of the activation barrier by the enzyme relative to the reaction in the gas phase is consistent with mutation experiments.

Complex flexibility of the transforming growth factor beta superfamily.

The transforming growth factors beta (TGF-beta s) are important modulators of growth and differentiation. They are intermolecular disulfide-bonded homodimeric molecules. The monomer fold has a conserved cystine knot and lacks a hydrophobic core. The biological specificity of a given member of the family is believed to be determined by the conformational flexibility of the variable loop regions of the monomer. The monomer subunit assembly in the dimer is stabilized mainly by hydrophobic contacts and a few hydrogen bonds. Since these interactions are nondirectional, we examined subunit assemblies of TGF-beta by using conformational analysis. The different subunit assemblies in TGF-beta 2 dimer were characterized in terms of the intersubunit disulfide torsion. Our analyses show that the subunit assemblies fall into two states: the crystallographically observed gauche+conformation and the previously not reported gauche--conformation, both having almost identical interaction energies. Furthermore, there is significant flexibility in the subunit assembly within the gauche+ and the gauche- states of the disulfide bond. The monomer subunit assembly is independent of the variations about the loop regions. The variations in the loop regions, coupled with flexibility in the monomer assembly, lead to a complex flexibility in the dimer of the TGF-beta superfamily. For the TGF-beta superfamily, the cystine knot acts as a scaffold and complex flexibility provides for biological selectivity. Complex flexibility might provide an explanation for the diverse range of biological activities that these important molecules display.