Automated Docking of α-(1,4)- and α-(1,6)-Linked Glucosyl Trisaccharides in the Glucoamylase Active Site

Automated docking to explore subsite binding by glycoside hydrolase family 6 cellobiohydrolases and endoglucanases

Cellooligosaccharides were computationally docked using AutoDock into the active sites of the glycoside hydrolase Family 6 enzymes Hypocrea jecorina (formerly Trichoderma reesei) cellobiohydrolase and Thermobifida fusca endoglucanase. Subsite -2 exerts the greatest intermolecular energy in binding beta-glucosyl residues, with energies progressively decreasing to either side. Cumulative forces imparting processivity exerted by these two enzymes are significantly less than by the equivalent glycoside hydrolase Family 7 enzymes studied previously. Putative subsites -4, -3, +3, and +4 exist in H. jecorina cellobiohydrolase, along with putative subsites -4, -3, and +3 in T. fusca endoglucanase, but they are less important than subsites -2, -1, +1, and +2. In general, binding adds 3-7 kcal/mol to ligand intramolecular energies because of twisting of scissile glycosidic bonds. Distortion of beta-glucosyl residues to the (2)S(O) conformation by binding in subsite -1 adds approximately 7 kcal/mol to substrate intramolecular energies.

Thermodynamic characterization of the adsorption of selected chiral compounds on immobilized amyloglucosidase in liquid chromatography

Immobilized amyloglucosidase was used as a chiral stationary phase (CSP). First, the retention and enantioselectivity of several model chiral amines and acids were investigated. We found that this CSP was unable to separate the enantiomers of acids, though all selected amines could be resolved. The adsorption of (R)- and (S)-propranolol and its influence on column temperature and 2-propanol content in the eluent were then studied in detail, using a three-step methodology. The adsorption was first evaluated using Scatchard plots; thereafter, the adsorption was characterized in detail by calculating the adsorption energy distribution. With this model-independent information, a better judgment could be made of the possible adsorption models selected in the last step, the model fitting to the data. In the case examined, the bi-Langmuir model (containing nonselective and enantioselective sites) describes the system well. The retention of (R)- and (S)-propranolol at low temperatures increases with the content of 2-propanol in the eluent, due to the increased saturation capacity of the enantioselective sites. The retention is an enthalpy-driven process at both types of sites, whereas the enantioseparation is due to differences between the entropy changes of the two enantiomers at the enantioselective sites. The enthalpy of adsorption at the nonselective sites is almost identical at the two concentrations of 2-propanol in the eluent. Enantioselective adsorption, on the other hand, is more exothermic at higher modifier content (20%). Thus, at high temperatures the retention decreases with increasing modifier content, whereas the opposite (unusual) trend is the case at low temperatures.

Docking studies on glycoside hydrolase Family 47 endoplasmic reticulum alpha-(1-->2)-mannosidase I to elucidate the pathway to the substrate transition state

Alpha-(1-->2)-mannosidase I from the endoplasmic reticulum (ERManI), a Family 47 glycoside hydrolase, is a key enzyme in the N-glycan synthesis pathway. Catalytic-domain crystal structures of yeast and human ERMan1s have been determined, the former with a hydrolytic product and the latter without ligands, with the inhibitors 1-deoxymannojirimycin and kifunensine, and with a thiodisaccharide substrate analog. Both inhibitors were bound at the base of the funnel-shaped active site as the unusual 1C4 conformer, while the substrate analog glycon is a 3S1 conformer. In the current study, AutoDock was used to dock alpha-D-mannopyranosyl-(1-->2)-alpha-D-mannopyranose with its glycon in chair (1C4,4C1), half-chair (3H2,3H4,4H3), skew-boat (OS2,3S1,5S1), boat (2,5B,3,OB,B1,4,B2,5), and envelope (3E,4E,E3,E4) conformations into the yeast ERManI active site. Both docked energies and forces on docked ligand atoms were calculated to determine how the ligand distorts to the transition state. From these, we can conclude that (1) both 1C4 and OS2 can be the starting conformers; (2) the most likely binding pathway is 1C4-->3H2-->OS2-->3,OB-->3S1-->3E; (3) the transition state is likely to be close to a 3E conformation.

Modeling protein recognition of carbohydrates

We have a limited understanding of the details of molecular recognition of carbohydrates by proteins, which is critical to a multitude of biological processes. Furthermore, carbohydrate-modifying proteins such as glycosyl hydrolases and phosphorylases are of growing importance as potential drug targets. Interactions between proteins and carbohydrates have complex thermodynamics, and in general the specific positioning of only a few hydroxyl groups determines their binding affinities. A thorough understanding of both carbohydrate and protein structures is thus essential to predict these interactions. An atomic-level view of carbohydrate recognition through structures of carbohydrate-active enzymes complexed with transition-state inhibitors reveals some of the distinctive molecular features unique to protein-carbohydrate complexes. However, the inherent flexibility of carbohydrates and their often water-mediated hydrogen bonding to proteins makes simulation of their complexes difficult. Nonetheless, recent developments such as the parameterization of specific force fields and docking scoring functions have greatly improved our ability to predict protein-carbohydrate interactions. We review protein-carbohydrate complexes having defined molecular requirements for specific carbohydrate recognition by proteins, providing an overview of the different computational techniques available to model them.

Visualizing complexes of phospholipids with Streptomyces phospholipase D by automated docking

The automated docking program AutoDock was used to dock nine phosphatidic acids (PAs), six phosphatidylcholines, five phosphatidylethanolamines, four phosphatidylglycerols, one phosphatidylinositol and two phosphatidylserines, which have two identical saturated fatty acid residues with an even numbers of carbon atoms, onto the active site of Streptomyces sp. PMF phospholipase D (PLD). Two PAs with one double bond on the fatty acid chain linked to the C2 of the glycerol residue were also docked. In general, binding energies become progressively more negative as fatty acid residues become longer. When these residues are of sufficient length, one is coiled against a hydrophobic cliff in a well that also holds the glycerol and phosphate residues and the head group, while the other generally is bound by a hydrophobic surface outside the well. Phosphatidylcholines have the only head group that is firmly bound by the active site, giving a possible structural explanation for the low selectivity of Streptomyces PLD for other phospholipid substrates.

Understanding protein structure-function relationships in Family 47 alpha-1,2-mannosidases through computational docking of ligands

Family 47 alpha-1,2-mannosidases are crucial enzymes involved in N-glycan maturation in the endoplasmic reticula and Golgi apparati of eukaryotic cells. High-resolution crystal structures of the human and yeast endoplasmic reticulum alpha-1,2-mannosidases have been recently determined, the former complexed with the inhibitors 1-deoxymannojirimycin and kifunensine, both of which bind in its active site in the unusual 1C4 conformation. However, unambiguous identification of the catalytic proton donor and nucleophile involved in glycoside bond hydrolysis was not possible from this structural information. In this work, alpha-D-galactose, alpha-D-glucose, and alpha-D-mannose were computationally docked in the active site in the energetically stable 4C1 conformation as well as in the 1C4 conformation to compare their interaction energetics. From these docked structures, a model for substrate and conformer selectivity based on the dimensions of the active site was proposed. Alpha-D-galactopyranosyl-(1-->2)-alpha-D-mannopyranose, alpha-D-glucopyranosyl-(1-->2)-alpha-D-mannopyranose, and alpha-D-mannopyranosyl-(1-->2)-alpha-D-mannopyranose were also docked into the active site with their nonreducing-end residues in the 1C4 and E4 (representing the transition state) conformations. Based on the docked structure of alpha-D-mannopyranosyl-E4-(1-->2)-alpha-D-mannopyranose, the catalytic acid and base are Glu132 and Glu435, respectively.

Automated docking of alpha-(1-->4)- and alpha-(1-->6)-linked glucosyl trisaccharides and maltopentaose into the soybean beta-amylase active site

The Lamarckian genetic algorithm of AutoDock 3.0 was used to dock alpha-maltotriose, methyl alpha-panoside, methyl alpha-isopanoside, methyl alpha-isomaltotrioside, methyl alpha-(6(1)-alpha-glucopyranosyl)-maltoside, and alpha-maltopentaose into the closed and, except for alpha-maltopentaose, into the open conformation of the soybean beta-amylase active site. In the closed conformation, the hinged flap at the mouth of the active site closes over the substrate. The nonreducing end of alpha-maltotriose docks preferentially to subsites -2 or +1, the latter yielding nonproductive binding. Some ligands dock into less optimal conformations with the nonreducing end at subsite -1. The reducing-end glucosyl residue of nonproductively-bound alpha-maltotriose is close to residue Gln194, which likely contributes to binding to subsite +3. In the open conformation, the substrate hydrogen-bonds with several residues of the open flap. When the flap closes, the substrate productively docks if the nonreducing end is near subsites -2 or -1. Trisaccharides with alpha-(1-->6) bonds do not successfully dock except for methyl alpha-isopanoside, whose first and second glucosyl rings dock exceptionally well into subsites -2 and -1. The alpha-(1-->6) bond between the second and third glucosyl units causes the latter to be improperly positioned into subsite +1; the fact that isopanose is not a substrate of beta-amylase indicates that binding to this subsite is critical for hydrolysis.

Automated docking of maltose, 2-deoxymaltose, and maltotetraose into the soybean beta-amylase active site

In this study, products and substrates were docked into the active site of beta-amylase using the simulated annealing algorithm AutoDock. Lowest-energy conformers reproduced known crystallographic atom positions within 0.4 to 0.8 A rmsd. Docking studies were carried out with both open and closed configurations of the beta-amylase mobile flap, a loop comprising residues 96 to 103. Ligands with two rings docked within the cleft near the active site when the flap was open, but those with four rings did not. The flap must be closed for alpha-maltotetraose to adopt a conformation allowing it to dock near the crystallographically determined subsites. The closed flap is necessary for productive but not for nonproductive binding, and therefore it plays a essential role in catalysis. The gain in total binding energy upon closing of the flap for alpha-maltose docked to subsites -2, -1 and +1, +2 is about 22 kcal/mol, indicating more favorable interactions are possible with the flap closed. Larger intermolecular interaction energies are observed for two alpha-maltose molecules docked to subsites -2, -1 and +1, +2 than for one alpha-maltotetraose molecule docked from subsites -2 to +2, suggesting that it is only upon cleavage of the alpha-1,4 linkage that optimal closed-flap binding can occur with the crytallographically determined enzyme structure.