Molecular modeling of the enantioselectivity in lipase-catalyzed transesterification reactions.

Two strategies based on the use of subsets for calculating the enantioselectivity in lipase-catalyzed transesterifications using the CHARMM force field were investigated. Molecular dynamics was used in our search for low energy conformations. Molecular mechanics was used for refining these low energy conformations. A tetrahedral intermediate with a rigid central part was used for mimicking the transition state. The energy differences between the transition states of the diastereomeric enzyme-substrate complexes were calculated. The way of defining the subsets was based on two fundamentally different strategies. The first strategy used predefined parts of the enzyme and the substrate as subsets. The second approach formed energy-based subsets, varying in size with the substrates studied. The selection of residues to be included in these energy-based subsets was based on the energy of the interaction between the specific residue or water molecule and the transition state. The reaction studied was the kinetic resolution of secondary alcohols in transesterifications using the Candida antarctica lipase B as chiral biocatalyst. The secondary alcohols used in the study were 2-butanol, 3-methyl-2-butanol, and 3,3-dimethyl-2-butanol.