Molecular modeling of the structures of human and rat pancreatic cholesterol esterases

Structural models have been generated for rat and human cholesterol esterases by molecular modeling. For rat cholesterol esterase, three separate models were generated according to the following procedure: (1) the cholesterol esterase sequence was aligned with those of three template enzymes: Torpedo californica acetylcholinesterase, Geotrichum candidum lipase and Candida rugosa lipase; (2) the X-ray structure coordinates of the three template enzymes were used to construct cholesterol esterase models by amino acid replacements of matched sequence positions and by making sequence insertions and deletions as required; (3) bad contracts in each of the cholesterol esterase models were relaxed by molecular dynamics and mechanics; (4) the three cholesterol esterase models were merged into one by arithmetic averaging of atomic coordinates; (5) Ramachandran analysis indicated that the model generated from the AChE template possessed the best set of phi/psi angles. Therefore, this model was subjected to molecular dynamics, with harmonic constraints imposed on the C(alpha) coordinates to drive them toward the coordinates of the averaged model. (6) Subsequent relaxation by molecular mechanics produced the final rat cholesterol esterase model. A model for human cholesterol esterase was produced by repeating steps 1-3 above, albeit with the rat cholesterol esterase model as the template. Hydrophobic and electrostatic analyses of the rat and human cholesterol esterase models suggest the structural origins of molecular recognition of hydrophobic substrates and interfaces, of charged interfaces, and of bile salt activators.