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

Functional and Structural Diversity of Acyl-coA Binding Proteins in Oil Crops

Diversities in structure and function of ACBP were discussed in this review. ACBP are important proteins that could transport newly synthesized fatty acid, activated into -coA, from plastid to endoplasmic reticulum, where oil in the form of triacylglycerol occurs. ACBP were detected in various animal and plants species, which indicated their importance in biological function. In fact, involvement of ACBP in important process such as lipid metabolism, regulation of enzyme and gene expression, and in response to plant stresses has been proven in several studies. In this review, findings on ACBP of 11 well-known oil crops were reviewed to comprehend diversity, comparative analyses on ACBP structure were made, and link between structure and function, tissue expression and subcellular location of ACBP were also observed. Incomplete reports in some species were mentioned, which might be encouraging to start or to perform deeper studies. Similar characteristics were found in paralogs ACBP, and orthologs ACBP had different functions, despite the high identity in amino acid sequence. At the end, it is confirmed that ortholog proteins could not necessarily display the same function, even from closely related species.

The first crystal structure of a family 129 glycoside hydrolase from a probiotic bacterium reveals critical residues and metal cofactors

The α-N-acetylgalactosaminidase from the probiotic bacterium Bifidobacterium bifidum (NagBb) belongs to the glycoside hydrolase family 129 and hydrolyzes the glycosidic bond of Tn-antigen (GalNAcα1-Ser/Thr). NagBb is involved in assimilation of O-glycans on mucin glycoproteins by B. bifidum in the human gastrointestinal tract, but its catalytic mechanism has remained elusive because of a lack of sequence homology around putative catalytic residues and of other structural information. Here we report the X-ray crystal structure of NagBb, representing the first GH129 family structure, solved by the single-wavelength anomalous dispersion method based on sulfur atoms of the native protein. We determined ligand-free, GalNAc, and inhibitor complex forms of NagBb and found that Asp-435 and Glu-478 are located in the catalytic domain at appropriate positions for direct nucleophilic attack at the anomeric carbon and proton donation for the glycosidic bond oxygen, respectively. A highly conserved Asp-330 forms a hydrogen bond with the O4 hydroxyl of GalNAc in the -1 subsite, and Trp-398 provides a stacking platform for the GalNAc pyranose ring. Interestingly, a metal ion, presumably Ca²⁺, is involved in the recognition of the GalNAc N-acetyl group. Mutations at Asp-435, Glu-478, Asp-330, and Trp-398 and residues involved in metal coordination (including an all-Ala quadruple mutant) significantly reduced the activity, indicating that these residues and the metal ion play important roles in substrate recognition and catalysis. Interestingly, NagBb exhibited some structural similarities to the GH101 endo-α-N-acetylgalactosaminidases, but several critical differences in substrate recognition and reaction mechanism account for the different activities of these two enzymes.

Molecular docking of carbohydrate ligands to antibodies: structural validation against crystal structures

Cell surface glycoproteins play vital roles in cellular homeostasis and disease. Antibody recognition of glycosylation on different cells and pathogens is critically important for immune surveillance. Conversely, adverse immune reactions resulting from antibody-carbohydrate interactions have been implicated in the development of autoimmune diseases and impact areas such as xenotransplantation and cancer treatment. Understanding the nature of antibody-carbohydrate interactions and the method by which saccharides fit into antibody binding sites is important in understanding the recognition process. In silico techniques offer attractive alternatives to experimental methods (X-ray crystallography and NMR) for the study of antibody-carbohydrate complexes. In particular, molecular docking provides information about protein-ligand interactions in systems that are difficult to study with experimental techniques. Before molecular docking can be used to investigate antibody-carbohydrate complexes, validation of an appropriate docking method is required. In this study, four popular docking programs, Glide, AutoDock, GOLD, and FlexX, were assessed for their ability to accurately dock carbohydrates to antibodies. Comparison of top ranking poses with crystal structures highlighted the strengths and weaknesses of these programs. Rigid docking, in which the protein conformation remains static, and flexible docking, where both the protein and ligand are treated as flexible, were compared. This study has revealed that generally molecular docking of carbohydrates to antibodies has been performed best by Glide.

Computational analyses of the conformational itinerary along the reaction pathway of GH94 cellobiose phosphorylase

GH94 cellobiose phosphorylase (CBP) catalyzes the phosphorolysis of cellobiose into alpha-D-glucose 1-phosphate (G1P) and D-glucose with inversion of anomeric configuration. The complex crystal structure of CBP from Cellvibrio gilvus had previously been determined; glycerol, glucose, and phosphate are bound to subsites -1, +1, and the anion binding site, respectively. We performed computational analyses to elucidate the conformational itinerary along the reaction pathway of this enzyme. autodock was used to dock cellobiose with its glycon glucosyl residue in various conformations and with its aglycon glucosyl residue in the low-energy 4C1 conformer. An oxocarbenium ion-like glucose molecule mimicking the transition state was also docked. Based on the clustering analysis, docked energies, and comparison with the crystallographic ligands, we conclude that the reaction proceeds from 1S3 as the pre-transition state conformer (Michaelis complex) via E3 as the transition state candidate to 4C1 as the G1P product conformer. The predicted reaction pathway of the inverting phosphorylase is similar to that proposed for the first-half glycosylation reaction of retaining cellulases, but is different from those for inverting cellulases. NAMD was used to simulate molecular dynamics of the enzyme. The 1S3 pre-transition state conformer is highly stable compared with other conformers, and a conformational change from 4C1 to 1,4B was observed.

Systematic comparison of catalytic mechanisms of hydrolysis and transfer reactions classified in the EzCatDB database

Catalytic mechanisms of 270 enzymes from 131 superfamilies, mainly hydrolases and transferases, were analyzed based on their enzyme structures. A method of systematic comparison and classification of the catalytic reactions was developed. Hydrolysis and transfer reactions closely resemble one another, displaying common mechanisms, single displacement, and double displacement. These displacement mechanisms might be further subclassified according to the type of catalytic factors and nucleophilic substitution involved. Several types of catalytic factors exist: nucleophile, acid, base, stabilizer, modulator, cofactors. Nucleophilic substitution might be categorized as S(N)1/S(N)2 (or dissociative/associative) reactions. The classification indicates that some mechanisms favor particular types of catalytic factors. In hydrolyses of amide bonds and phosphoric ester bonds, mechanisms with single displacement tend to use inorganic cofactors such as zinc and magnesium ions as important catalysts, whereas those with double displacement frequently do not use such cofactors. In contrast, hydrolyses of O-glycoside bond rarely use such cofactors, with one exception. The trypsin-like hydrolytic reaction, which is catalyzed by the classic catalytic triad comprising serine/histidine/aspartate, can be considered as a "super-reaction" because it is observed in at least three nonhomologous enzymes, whereas most reactions are singlets without any nonhomologous enzymes. By dividing complex reactions into several reactions, correlations between active site structures and catalytic functions can be suggested. This classification method is applicable to other reactions such as elimination and isomerization. Furthermore, it will facilitate annotation of enzyme functions from 3D patterns of enzyme active sites. The classification is available at http://mbs.cbrc.jp/EzCatDB/RLCP/index.html.

The fate of β-d-mannopyranose after its formation by endoplasmic reticulum α-(1→2)-mannosidase I catalysis

The automated docking program AutoDock was used to dock all 38 characteristic beta-D-mannopyranose ring conformers into the active site of the yeast endoplasmic reticulum alpha-(1-->2)-mannosidase I, a Family 47 glycoside hydrolase that converts Man9GlcNAc2 to Man8GlcNAc2. The subject of this work is to establish the conformational pathway that allows the cleaved glycon product to leave the enzyme active site and eventually reach the ground-state conformation. Twelve of the 38 conformers optimally dock in the active site where the inhibitors 1-deoxymannonojirimycin and kifunensine are found in enzyme crystal structures. A further 23 optimally dock in a second site on the side of the active-site well, while three dock outside the active-site cavity. It appears, through analysis of the internal energies of different ring conformations, of intermolecular energies between the ligands and enzyme, and of forces exerted on the ligands by the enzyme, that beta-D-mannopyranose follows the path 3E-->1C4-->1H2-->B2,5 before being expelled by the enzyme. The highly conserved second site that strongly binds beta-D-mannopyranose-4C1 may exist to prevent competitive inhibition by the product, and is worthy of further investigation.

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.

Multiple ligand-binding modes in bacterial R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR), a bacterial plasmid-encoded enzyme associated with resistance to the drug trimethoprim, shows neither sequence nor structural homology with the chromosomal DHFR. It presents a highly symmetrical toroidal structure, where four identical monomers contribute to the unique central active-site pore. Two reactants (dihydrofolate, DHF), two cofactors (NADPH) or one of each (R67*DHF*NADPH) can be found simultaneously within the active site, the last one being the reactive ternary complex. As the positioning of the ligands has proven elusive to empirical determination, we addressed the problem from a theoretical perspective. Several potential structures of the ternary complex were generated using the docking programs AutoDock and FlexX. The variability among the final poses, many of which conformed to experimental data, prompted us to perform a comparative scoring analysis and molecular dynamics simulations to assess the stability of the complexes. Analysis of ligand-ligand and ligand-protein interactions along the 4 ns trajectories of eight different structures allowed us to identify important inter-ligand contacts and key protein residues. Our results, combined with published empirical data, clearly suggest that multipe binding modes of the ligands are possible within R67 DHFR. While the pterin ring of DHF and the nicotinamide ring of NADPH assume a stacked endo-conformation at the centre of the pore, probably assisted by V66, Q67 and I68, the tails of the molecules extend towards opposite ends of the cavity, adopting multiple configurations in a solvent rich-environment where hydrogen-bond interactions with K32 and Y69 may play important roles.

Force calculations in automated docking: Enzyme-substrate interactions in Fusarium oxysporum Cel7B

AutoDock is a small-molecule docking program that uses an energy function to score docked ligands. Here AutoDock's grid-based method for energy evaluation was exploited to evaluate the force exerted by Fusarium oxysporum Cel7B on the atoms of docked cellooligosaccharides and a thiooligosaccharide substrate analog. Coupled with the interaction energies evaluated for each docked ligand, these forces give insight into the dynamics of the ligand in the active site, and help to elucidate the relative importance of specific enzyme-substrate interactions in stabilizing the substrate transition-state conformation. The processive force on the docked substrate in the F. oxysporum Cel7B active site is less than half of that on the docked substrate in the Hypocrea jecorina Cel7A active site. Hydrogen bonding interactions of the enzyme with the C2 hydroxyl group of the glucosyl residue in subsite -2 and with the C3 hydroxyl group of the glucosyl residue in subsite +1 are the most significant in stabilizing the distorted14B transition-state conformation of the glucosyl residue in subsite -1. The force calculations also help to elucidate the mechanism that prevents the active site from fouling.

Hypocrea jecorina (Trichoderma reesei) Cel7A as a molecular machine: A docking study

Hypocrea jecorina (formerly Trichoderma reesei) Cel7A has a catalytic domain (CD) and a cellulose-binding domain (CBD) separated by a highly glycosylated linker. Very little is known of how the 2 domains interact to degrade crystalline cellulose. Based on the interaction energies and forces on cello-oligosaccharides computationally docked to the CD and CBD, we propose a molecular machine model, where the CBD wedges itself under a free chain end on the crystalline cellulose surface and feeds it to the CD active site tunnel. Enzyme-substrate interactions produce the forces required to pull cellulose chains from the surface and also to help the enzyme move on the cellulose chain for processive hydrolysis. The energy to generate these forces is ultimately derived from the chemical energy of glycosidic bond breakage.

Mechanism of Class 1 (Glycosylhydrolase Family 47) α-Mannosidases Involved in N-Glycan Processing and Endoplasmic Reticulum Quality Control

Quality control in the endoplasmic reticulum (ER) determines the fate of newly synthesized glycoproteins toward either correct folding or disposal by ER-associated degradation. Initiation of the disposal process involves selective trimming of N-glycans attached to misfolded glycoproteins by ER alpha-mannosidase I and subsequent recognition by the ER degradation-enhancing alpha-mannosidase-like protein family of lectins, both members of glycosylhydrolase family 47. The unusual inverting hydrolytic mechanism catalyzed by members of this family is investigated here by a combination of kinetic and binding analyses of wild type and mutant forms of human ER alpha-mannosidase I as well as by structural analysis of a co-complex with an uncleaved thiodisaccharide substrate analog. These data reveal the roles of potential catalytic acid and base residues and the identification of a novel (3)S(1) sugar conformation for the bound substrate analog. The co-crystal structure described here, in combination with the (1)C(4) conformation of a previously identified co-complex with the glycone mimic, 1-deoxymannojirimycin, indicates that glycoside bond cleavage proceeds through a least motion conformational twist of a properly predisposed substrate in the -1 subsite. A novel (3)H(4) conformation is proposed as the exploded transition state.

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.