Search

This study defined chronic recurrent cervical nerve root neurapraxia, the chronic burner syndrome, characterized the clinical findings, and described the responsible pathomechanics. We studied a subset of 55 athletes (mean age, 22 years) for evaluation of recurrent burners. Eleven subjects were professional athletes. The mechanism of injury was extension combined with ipsilateral-lateral deviation in 46 patients (83%). Spurling's sign was positive in 39 patients (70%). Twenty-nine patients (53%) had developmentally narrowed cervical canals, and 48 patients (87%) had evidence of disk disease by magnetic resonance imaging. The disk disease was in the form of a disk bulge, disk protrusion, or a frank disk herniation deforming the cord. Fifty-one patients (93%) had disk disease or narrowing of the intervertebral foramina secondary to degenerative disk disease. Although burners may be the result of a brachial plexus stretch injury in high school and collegiate football players seen with acute symptoms, nerve root compression in the intervertebral foramina secondary to disk disease is a more common cause in collegiate and professional players who have recurrent or chronic burner syndromes. There is a high incidence of cervical canal stenosis in football players with recurrent burner syndrome. The combination of disk disease and cervical spinal canal stenosis may lead to an alteration in normal cervical spine mechanics that may make these athletes more prone to chronic burner syndromes. [ABSTRACT FROM AUTHOR]

We evaluated 55 knees in 51 patients after Elmslie-Trillat-Maquet procedures. The procedure involves medialization of the tibial tubercle on a distal pedicle and elevating the tibial tubercle anteriorly 10 mm with a local bone graft. At a mean followup of 74.2 months (range, 13 to 196), all patients completed postoperative surveys and 38 underwent postoperative examinations. Subjectively, 9 knees (16%) had excellent results, 24 knees (44%) obtained good results, and 13 knees (24%) had fair results for a total of 84% improvement overall. Using Fulkerson's functional knee score, 19 knees (35%) had excellent results, 10 knees (18%) had good results, and 11 knees (20%) had fair results for a total of 73% improvement overall. A total of 24 knees (44%) required later screw removal. The most significant findings of this study include 1) an 84% overall subjective improvement in symptoms; 2) the findings that young patients without evidence of progressive osteoarthrosis and with patella instability as a primary symptom tend to have the most favorable outcome; and 3) 24 knees (44%) required later screw removal. [ABSTRACT FROM AUTHOR]

Background . The antibiotic resistance patterns of young infants with urinary tract infections (UTIs) have evolved over the past 2 decades. Whether current empiric antibiotic regimens are sufficient in this age group is unknown. Methods . A retrospective review of patients aged 0 to 60 days admitted with a UTI discharge diagnosis. Results . Overall susceptibility to empiric antibiotics was 87%. Antibiotic resistance and length of stay were highest among those who were afebrile, those admitted to the intensive care unit, and those with culture diagnosis of enterococcal infection. The sensitivity and specificity of ultrasound as a screening tool for genitourinary anomaly was 70% and 40%, respectively, with a positive predictive value of 31.8%. Conclusions . Empiric antibiotic regimens cover a high percentage of UTIs in infants. However, high rates of resistance and prolonged length of stay in patients with enterococcal infection highlight the need for continued surveillance of such patients in this age group., Competing Interests: Declaration of Conflicting Interests: The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

We classified the carboxylic ester hydrolases (CEHs) into families and clans by use of multiple sequence alignments, secondary structure analysis, and tertiary structure superpositions. Our work for the first time has fully established their systematic structural classification. Family members have similar primary, secondary, and tertiary structures, and their active sites and reaction mechanisms are conserved. Families may be gathered into clans by their having similar secondary and tertiary structures, even though primary structures of members of different families are not similar. CEHs were gathered from public databases by use of Basic Local Alignment Search Tool (BLAST) and divided into 91 families, with 36 families being grouped into five clans. Members of one clan have standard α/β-hydrolase folds, while those of other two clans have similar folds but with different sequences of their β-strands. The other two clans have members with six-bladed β-propeller and three-α-helix bundle tertiary structures. Those families not in clans have a large variety of structures or have no members with known structures. At the time of writing, the 91 families contained 321,830 primary structures and 1378 tertiary structures. From these data, we constructed an accessible database: CASTLE (CArboxylic eSTer hydroLasEs, http://www.castle.cbe.iastate.edu)., (© 2016 The Protein Society.)

At present, 69 families of carbohydrate-binding modules (CBMs) have been isolated by statistically significant differences in the amino acid sequences (primary structures) of their members, with most members of different families showing little if any homology. On the other hand, members of the same family have primary and tertiary (three-dimensional) structures that can be computationally aligned, suggesting that they are descended from common protein ancestors. Members of the large majority of CBM families are β-sandwiches. This raises the question of whether members of different families are descended from distant common ancestors, and therefore are members of the same tribe. We have attacked this problem by attempting to computationally superimpose tertiary structure representatives of each of the 53 CBM families that have members with known tertiary structures. When successful, we have aligned locations of secondary structure elements and determined root mean square deviations and percentages of similarity between adjacent amino acid residues in structures from similar families. Further criteria leading to tribal membership are amino acid chain lengths and bound ligands. These considerations have led us to assign 27 families to nine tribes. Eight of the tribes have members with β-sandwich structures, while the ninth is composed of structures with β-trefoils., (© 2014 Wiley Periodicals, Inc.)

Ketoacyl reductases (KRs), hydroxyacyl dehydratases (HDs), and enoyl reductases (ERs) are part of the fatty acid/polyketide synthesis cycle. They are known as acyl dehydrogenases, enoyl hydratases, and hydroxyacyl dehydrogenases, respectively, when catalyzing their reverse reactions. Earlier, we classified these enzymes into four KR, eight HD, and five ER families by statistical criteria. Members of all four KR families and three ER families have Rossmann folds, while five HD family members have HotDog folds. This suggests that those proteins with the same folds in different families may be distantly related, and therefore in clans, even though their amino acid sequences may not be homologous. We have now defined two clans containing three of the four KR families and two of the eight HD families, using manual and statistical tests. One of the ER families is related to the KR clan.

Automated docking allows rapid screening of protein-ligand interactions. A scoring function composed of a force field and linear weights can be used to compute a binding energy from a docked atom configuration. For different force fields or types of molecules, it may be necessary to train a custom scoring function. This chapter describes the data and methods one must consider in developing a custom scoring function for use with AutoDock.

Thioesterases are enzymes that hydrolyze thioester bonds between a carbonyl group and a sulfur atom. They catalyze key steps in fatty acid biosynthesis and metabolism, as well as polyketide biosynthesis. The reaction molecular mechanism of most hotdog-fold acyl-CoA thioesterases remains unknown, but several hypotheses have been put forward in structural and biochemical investigations. The reaction of a human thioesterase (hTHEM2), representing a thioesterase family with a hotdog fold where a coenzyme A moiety is cleaved, was simulated by quantum mechanics/molecular mechanics metadynamics techniques to elucidate atomic and electronic details of its mechanism, its transition-state conformation, and the free energy landscape of the process. A single-displacement acid-base-like mechanism, in which a nucleophilic water molecule is activated by an aspartate residue acting as a base, was found, confirming previous experimental proposals. The results provide unambiguous evidence of the formation of a tetrahedral-like transition state. They also explain the roles of other conserved active-site residues during the reaction, especially that of a nearby histidine/serine pair that protonates the thioester sulfur atom, the participation of which could not be elucidated from mutation analyses alone., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Ketoacyl reductases (KRs), hydroxyacyl dehydratases (HDs) and enoyl reductases (ERs) are part of the fatty acid and polyketide synthesis cycles. Their reverse reactions, catalyzed by acyl dehydrogenases (equivalent to ERs), enoyl hydratases (equivalent to HDs) and hydroxyacyl dehydrogenases (equivalent to KRs), are part of fatty acid degradation by β-oxidation. These enzymes have been classified into families based on similarities in their primary and tertiary structures, and these families and their structures are included in the ThYme (Thioester-active enzYmes) database. Members of each family have strong sequence similarity and have essentially the same tertiary structure, mechanism and catalytic residues.

We gathered primary and tertiary structures of acyl-CoA carboxylases from public databases, and established that members of their biotin carboxylase (BC) and biotin carboxyl carrier protein (BCCP) domains occur in one family each and that members of their carboxyl transferase (CT) domains occur in two families. Protein families have members similar in primary and tertiary structure that probably have descended from the same protein ancestor. The BCCP domains complexed with biotin in acyl and acyl-CoA carboxylases transfer bicarbonate ions from BC domains to CT domains, enabling the latter to carboxylate acyl and acyl-CoA moieties. We separated the BCCP domains into four subfamilies based on more subtle primary structure differences. Members of different BCCP subfamilies often are produced by different types of organisms and are associated with different carboxylases.

All acyl carrier protein primary and tertiary structures were gathered into the ThYme database. They are classified into 16 families by amino acid sequence similarity, with members of the different families having sequences with statistically highly significant differences. These classifications are supported by tertiary structure superposition analysis. Tertiary structures from a number of families are very similar, suggesting that these families may come from a single distant ancestor. Normal vibrational mode analysis was conducted on experimentally determined freestanding structures, showing greater fluctuations at chain termini and loops than in most helices. Their modes overlap more so within families than between different families. The tertiary structures of three acyl carrier protein families that lacked any known structures were predicted as well., (Copyright © 2012 The Protein Society.)

Ketoacyl synthases (KSs) catalyze condensing reactions combining acyl-CoA or acyl-acyl carrier protein (acyl-ACP) with malonyl-CoA to form 3-ketoacyl-CoA or with malonyl-ACP to form 3-ketoacyl-ACP. In each case, the resulting acyl chain is two carbon atoms longer than before, and CO2 and either CoA or ACP are formed. KSs also join other activated molecules in the polyketide synthesis cycle. Our classification of KSs by their primary and tertiary structures instead of by their substrates and the reactions that they catalyze enhances insights into this enzyme group. KSs fall into five families separated by their characteristic primary structures, each having members with the same catalytic residues, mechanisms, and tertiary structures. KS1 members, overwhelmingly named 3-ketoacyl-ACP synthase III or its variants, are produced predominantly by bacteria. Members of KS2 are mainly produced by plants, and they are usually long-chain fatty acid elongases/condensing enzymes and 3-ketoacyl-CoA synthases. KS3, a very large family, is composed of bacterial and eukaryotic 3-ketoacyl-ACP synthases I and II, often found in multidomain fatty acid and polyketide synthases. Most of the chalcone synthases, stilbene synthases, and naringenin-chalcone synthases in KS4 are from eukaryota. KS5 members are all from eukaryota, most are produced by animals, and they are mainly fatty acid elongases. All families except KS3 are split into subfamilies whose members have statistically significant differences in their primary structures. KS1 through KS4 appear to be part of the same clan. KS sequences, tertiary structures, and family classifications are available on the continuously updated ThYme (Thioester-active enzYme) database.

Background: Acyl-acyl carrier protein thioesterases (acyl-ACP TEs) catalyze the hydrolysis of the thioester bond that links the acyl chain to the sulfhydryl group of the phosphopantetheine prosthetic group of ACP. This reaction terminates acyl chain elongation of fatty acid biosynthesis, and in plant seeds it is the biochemical determinant of the fatty acid compositions of storage lipids., Results: To explore acyl-ACP TE diversity and to identify novel acyl ACP-TEs, 31 acyl-ACP TEs from wide-ranging phylogenetic sources were characterized to ascertain their in vivo activities and substrate specificities. These acyl-ACP TEs were chosen by two different approaches: 1) 24 TEs were selected from public databases on the basis of phylogenetic analysis and fatty acid profile knowledge of their source organisms; and 2) seven TEs were molecularly cloned from oil palm (Elaeis guineensis), coconut (Cocos nucifera) and Cuphea viscosissima, organisms that produce medium-chain and short-chain fatty acids in their seeds. The in vivo substrate specificities of the acyl-ACP TEs were determined in E. coli. Based on their specificities, these enzymes were clustered into three classes: 1) Class I acyl-ACP TEs act primarily on 14- and 16-carbon acyl-ACP substrates; 2) Class II acyl-ACP TEs have broad substrate specificities, with major activities toward 8- and 14-carbon acyl-ACP substrates; and 3) Class III acyl-ACP TEs act predominantly on 8-carbon acyl-ACPs. Several novel acyl-ACP TEs act on short-chain and unsaturated acyl-ACP or 3-ketoacyl-ACP substrates, indicating the diversity of enzymatic specificity in this enzyme family., Conclusion: These acyl-ACP TEs can potentially be used to diversify the fatty acid biosynthesis pathway to produce novel fatty acids.

Two forms of Ruminococcus flavefaciens FD-1 endoglucanase B, a member of glycoside hydrolase family 44, one with only a catalytic domain and the other with a catalytic domain and a carbohydrate binding domain (CBM), were produced. Both forms hydrolyzed cellotetraose, cellopentaose, cellohexaose, carboxymethylcellulose (CMC), birchwood and larchwood xylan, xyloglucan, lichenan, and Avicel but not cellobiose, cellotriose, mannan, or pullulan. Addition of the CBM increased catalytic efficiencies on both CMC and birchwood xylan but not on xyloglucan, and it decreased rates of cellopentaose and cellohexaose hydrolysis. Catalytic efficiencies were much higher on xyloglucan than on other polysaccharides. Hydrolysis rates increased with increasing cellooligosaccharide chain length. Cellotetraose hydrolysis yielded only cellotriose and glucose. Hydrolysis of cellopentaose gave large amounts of cellotetraose and glucose, somewhat more of the former than of the latter, and much smaller amounts of cellobiose and cellotriose. Cellohexaose hydrolysis yielded much more cellotetraose than cellobiose and small amounts of glucose and cellotriose, along with a low and transient amount of cellopentaose., (Copyright © 2010 Elsevier Inc. All rights reserved.)

The ThYme (Thioester-active enzYme; http://www.enzyme.cbirc.iastate.edu) database has been constructed to bring together amino acid sequences and 3D (tertiary) structures of all the enzymes constituting the fatty acid synthesis and polyketide synthesis cycles. These enzymes are active on thioester-containing substrates, specifically those that are parts of the acyl-CoA synthase, acyl-CoA carboxylase, acyl transferase, ketoacyl synthase, ketoacyl reductase, hydroxyacyl dehydratase, enoyl reductase and thioesterase enzyme groups. These groups have been classified into families, members of which are similar in sequences, tertiary structures and catalytic mechanisms, implying common protein ancestry. ThYme is continually updated as sequences and tertiary structures become available.

Glycoside hydrolases cleave the glycosidic linkage between two carbohydrate moieties. They are among the most efficient enzymes currently known. β-Xylosidases from glycoside hydrolase family 43 hydrolyze the nonreducing ends of xylooligomers using an inverting mechanism. Although the general mechanism and catalytic amino acid residues of β-xylosidases are known, the nature of the reaction's transition state and the conformations adopted by the glycon xylopyranosyl ring along the reaction pathway are still elusive. In this work, the xylobiose hydrolysis reaction catalyzed by XynB3, a β-xylosidase produced by Geobacillus stearothermophilus T-6, was explicitly modeled using first-principles quantum mechanics/molecular mechanics Car-Parrinello metadynamics. We present the reaction's free energy surface and its previously undetermined reaction pathway. The simulations also show that the glycon xylopyranosyl ring proceeds through a (2,5)B-type transition state with significant oxacarbenium ion character.

Thioesterases (TEs) are classified into EC 3.1.2.1 through EC 3.1.2.27 based on their activities on different substrates, with many remaining unclassified (EC 3.1.2.-). Analysis of primary and tertiary structures of known TEs casts a new light on this enzyme group. We used strong primary sequence conservation based on experimentally proved proteins as the main criterion, followed by verification with tertiary structure superpositions, mechanisms, and catalytic residue positions, to accurately define TE families. At present, TEs fall into 23 families almost completely unrelated to each other by primary structure. It is assumed that all members of the same family have essentially the same tertiary structure; however, TEs in different families can have markedly different folds and mechanisms. Conversely, the latter sometimes have very similar tertiary structures and catalytic mechanisms despite being only slightly or not at all related by primary structure, indicating that they have common distant ancestors and can be grouped into clans. At present, four clans encompass 12 TE families. The new constantly updated ThYme (Thioester-active enzYmes) database contains TE primary and tertiary structures, classified into families and clans that are different from those currently found in the literature or in other databases. We review all types of TEs, including those cleaving CoA, ACP, glutathione, and other protein molecules, and we discuss their structures, functions, and mechanisms.

Golgi alpha-mannosidase II (GMII), a member of glycoside hydrolase family 38, cleaves two mannosyl residues from GlcNAcMan(5)GlcNAc(2) as part of the N-linked glycosylation pathway. To elucidate the molecular and electronic details of the reaction mechanism, in particular the conformation of the substrate at the transition state, we performed quantum mechanics/molecular mechanics metadynamics simulations of the glycosylation reaction catalyzed by GMII. The calculated free energy of activation for mannosyl glycosylation (23 kcal/mol) agrees very well with experiments, as does the conformation of the glycon mannosyl ring in the product of the glycosylation reaction (the covalent intermediate). In addition, we provide insight into the electronic aspects of the molecular mechanism that were not previously available. We show that the substrate adopts an (O)S(2)/B(2,5) conformation in the GMII Michaelis complex and that the nucleophilic attack occurs before complete departure of the leaving group, consistent with a D(N)A(N) reaction mechanism. The transition state has a clear oxacarbenium ion (OCI) character, with the glycosylation reaction following an (O)S(2)/B(2,5) --> B(2,5) [TS] --> (1)S(5) itinerary, agreeing with an earlier proposal based on comparing alpha- and beta-mannanases. The simulations also demonstrate that an active-site Zn ion helps to lengthen the O2'-H(O2') bond when the substrate acquires OCI character, relieving the electron deficiency of the OCI-like species. Our results can be used to explain the potency of recently formulated GMII anticancer inhibitors, and they are potentially relevant in deriving new inhibitors.

A gene encoding a glycoside hydrolase family 44 (GH44) protein from Clostridium acetobutylicum ATCC 824 was synthesized and transformed into Escherichia coli. The previously uncharacterized protein was expressed with a C-terminal His tag and purified by nickel-nitrilotriacetic acid affinity chromatography. Crystallization and X-ray diffraction to a 2.2-A resolution revealed a triose phosphate isomerase (TIM) barrel-like structure with additional Greek key and beta-sandwich folds, similar to other GH44 crystal structures. The enzyme hydrolyzes cellotetraose and larger cellooligosaccharides, yielding an unbalanced product distribution, including some glucose. It attacks carboxymethylcellulose and xylan at approximately the same rates. Its activity on carboxymethylcellulose is much higher than that of the isolated C. acetobutylicum cellulosome. It also extensively converts lichenan to oligosaccharides of intermediate size and attacks Avicel to a limited extent. The enzyme has an optimal temperature in a 10-min assay of 55 degrees C and an optimal pH of 5.0.

Patterns of scissile bond twisting have been found in crystal structures of glycoside hydrolases (GHs) that are complexed with substrates and inhibitors. To estimate the increased potential energy in the substrates that results from this twisting, we have plotted torsion angles for the scissile bonds on hybrid Quantum Mechanics::Molecular Mechanics energy surfaces. Eight such maps were constructed, including one for alpha-maltose and three for different forms of methyl alpha-acarviosinide to provide energies for twisting of alpha-(1,4) glycosidic bonds. Maps were also made for beta-thiocellobiose and for three beta-cellobiose conformers having different glycon ring shapes to model distortions of beta-(1,4) glycosidic bonds. Different GH families twist scissile glycosidic bonds differently, increasing their potential energies from 0.5 to 9.5 kcal/mol. In general, the direction of twisting of the glycosidic bond away from the conformation of lowest intramolecular energy correlates with the position (syn or anti) of the proton donor with respect to the glycon's ring oxygen atom. This correlation suggests that glycosidic bond distortion is important for the optimal orientation of one of the glycosidic oxygen lone pairs toward the enzyme's proton donor.

Two glycoside hydrolase (GH) families were analyzed to detect the presence of functional divergence using the program DIVERGE. These two families, GH7 and GH16, each contain members related by amino acid sequence similarity, retaining hydrolytic mechanisms, and catalytic residue identity. GH7 and GH16 comprise GH Clan B, with a shared beta-jelly roll topology and mechanism. GH7 contains fungal cellobiohydrolases and endoglucanases and is divided into five main subfamilies, four of the former and one of the latter. Cluster comparisons between three of the cellobiohydrolase subfamilies and the endoglucanase subfamily identified specific amino acid residues that play a role in the functional divergence between the two enzyme types. GH16 contains subfamilies of bacterial agarases, xyloglucosyl transferases, 1,3-beta-D-glucanases, lichenases, and other enzymes with various substrate specificities and product profiles. Four cluster comparisons between these four main subfamilies again have identified amino acid residues involved in functional divergence between the subfamilies. (c) 2009 Wiley Periodicals, Inc. Biopolymers 91: 478-495, 2009.

A detailed understanding of the catalytic strategy of cellulases is key to finding alternative ways to hydrolyze cellulose to mono-, di-, and oligosaccharides. Endoglucanases from glycoside hydrolase family 8 (GH8) catalyze the hydrolysis of beta-1,4-glycosidic bonds in cellulose by an inverting mechanism believed to involve a oxacarbenium ion-like transition state (TS) with a boat-type conformation of the glucosyl unit in subsite -1. In this work, hydrolysis by Clostridium thermocellum endo-1,4-glucanase A was computationally simulated with quantum mechanics/molecular mechanics metadynamics based on density functional theory. Our calculations show that the glucosyl residue in subsite -1 in the Michaelis complex is in a distorted 2SO/2,5B ring conformation, agreeing well with its crystal structure. In addition, our simulations capture the cationic oxacarbenium ion-like character of the TS with a partially formed double bond between the ring oxygen and C5' carbon atoms. They also provide previously unknown structural information of important states along the reaction pathway. The simulations clearly show for the first time in GH8 members that the TS features a boat-type conformation of the glucosyl unit in subsite -1. The overall catalytic mechanism follows a DN*AN-like mechanism and a beta-2SO-->2,5B [TS]-->alpha-5S1 conformational itinerary along the reaction coordinate, consistent with the anti-periplanar lone pair hypothesis. Because of the structural similarities and sequence homology among all GH8 members, our results can be extended to all GH8 cellulases, xylanases, and other endoglucanases. In addition, we provide evidence supporting the role of Asp278 as the catalytic proton acceptor (general base) for GH8a subfamily members.

Glycoside hydrolase family 1 consists of beta-glucosidases, beta-galactosidases, 6-phospho-beta-galactosidases, myrosinases, and other enzymes having similar primary and tertiary structures but diverse specificities. Among these enzymes, beta-glucosidases hydrolyze cellobiose to glucose, and therefore they are key players in any cellulose to glucose process. All family members attack beta-glycosidic bonds between a pyranosyl glycon and an aglycon, but most have little specificity for the aglycon or for the bond configuration. Furthermore, glycon specificity is not absolute. Sixteen family members (six beta-glucosidases, two cyanogenic beta-glucosidases, one 6-phospho-beta-galactosidase, two myrosinases, and five beta-glycosidases) have known tertiary structures. We have used automated docking to computationally bind disaccharides with allopyranosyl, galactopyranosyl, glucopyranosyl, mannopyranosyl, 6-phosphogalactopyranosyl, and 6-phosphoglucopyranosyl glycons, all linked by beta-(1,2), beta-(1,3), beta-(1,4), and beta-(1,6)-glycosidic bonds to beta-glucopyranoside aglycons, along with beta-(1,1-thio)-allopyranosyl, -galactopyranosyl, -glucopyranosyl, and -mannopyranosyl) beta-glucopyranosides, into all of these structures to investigate the structural determinants of their enzyme specificities. The following are the eight active-site residues: Glu191, Thr194, Phe205, Asn285, Arg336, Asn376, Trp378, and Trp465 (Zea mays beta-glucosidase numbering), that control a significant amount of glycon specificity.

It has been difficult to identify the proton donor and nucleophilic assistant/base of endoplasmic reticulum alpha-(1-->2)-mannosidase I, a member of glycoside hydrolase Family 47, which cleaves the glycosidic bond between two alpha-(1-->2)-linked mannosyl residues by the inverting mechanism, trimming Man(9)GlcNAc(2) to Man(8)GlcNAc(2) isomer B. Part of the difficulty is caused by the enzyme's use of a water molecule to transmit the proton that attacks the glycosidic oxygen atom. We earlier used automated docking to conclusively determine that Glu435 in the yeast enzyme (Glu599 in the corresponding human enzyme) is the nucleophilic assistant. The commonly accepted proton donor has been Glu330 in the human enzyme (Glu132 in the yeast enzyme). However, for theoretical reasons this conclusion is untenable. Theory, automated docking of alpha-d-(3)S(1)-mannopyranosyl-(1-->2)-alpha-d-(4)C(1)-mannopyranose and water molecules associated with candidate proton donors, and estimation of dissociation constants of the latter have shown that the true proton donor is Asp463 in the human enzyme (Asp275 in the yeast enzyme).

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.

Thermodynamic information can be inferred from static atomic configurations. To model the thermodynamics of carbohydrate binding to proteins accurately, a large binding data set has been assembled from the literature. The data set contains information from 262 unique protein-carbohydrate crystal structures for which experimental binding information is known. Hydrogen atoms were added to the structures and training conformations were generated with the automated docking program AutoDock 3.06, resulting in a training set of 225,920 all-atom conformations. In all, 288 formulations of the AutoDock 3.0 free energy model were trained against the data set, testing each of four alternate methods of computing the van der Waals, solvation, and hydrogen-bonding energetic components. The van der Waals parameters from AutoDock 1 produced the lowest errors, and an entropic model derived from statistical mechanics produced the only models with five physically and statistically significant coefficients. Eight models predict the Gibbs free energy of binding with an error of less than 40% of the error of any similar models previously published., (Copyright 2007 Wiley Periodicals, Inc.)

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.

Bovine rumen fluid was fermented anaerobically with 25 mM R-propylene glycol, S-propylene glycol, or glycerol added. After 24 h, all of the propylene glycol enantiomers and approximately 80% of the glycerol were metabolized. Acetate, propionate, butyrate, valerate, and caproate concentrations, in decreasing order, all increased with incubation time. Addition of any of the three substrates somewhat decreased acetate formation, while addition of either propylene glycol increased propionate formation but decreased that of butyrate. R- and S-propylene glycol did not differ significantly in either their rates of disappearance or the products formed when they were added to the fermentation medium. Fermentations of rumen fluid containing propylene glycol emitted the sulfur-containing gases 1-propanethiol, 1-(methylthio)propane, methylthiirane, 2,4-dimethylthiophene, 1-(methylthio)-1-propanethiol, dipropyl disulfide, 1-(propylthio)-1-propanethiol, dipropyl trisulfide, 3,5-diethyl-1,2,4-trithiolane, 2-ethyl-1,3-dithiane, and 2,4,6-triethyl-1,3,5-trithiane. Metabolic pathways that yield each of these gases are proposed. The sulfur-containing gases produced during propylene glycol fermentation in the rumen may contribute to the toxic effects seen in cattle when high doses are administered for therapeutic purposes.

We describe a new method of describing the pucker of an N-member monocyclic ring using N - 3 parameters. To accomplish this, three ring atoms define a reference plane, and the remainder of the ring is decomposed into triangular flaps. The angle of incidence for each flap upon the reference plane is then measured. The combination of these angles is characteristic of the ring's pucker. This puckering coordinate system is compared to existing reduced parameter systems to describe rings using a cyclohexane molecule. We show that this method has the same descriptive power of previous systems while offering advantages in molecular simulations.

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.

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.

Different programs and methods were employed to superimpose protein structures, using members of four very different protein families as test subjects, and the results of these efforts were compared. Algorithms based on human identification of key amino acid residues on which to base the superpositions were nearly always more successful than programs that used automated techniques to identify key residues. Among those programs automatically identifying key residues, MASS could not superimpose all members of some families, but was very efficient with other families. MODELLER, MultiProt, and STAMP had varying levels of success. A genetic algorithm program written for this project did not improve superpositions when results from neighbor-joining and pseudostar algorithms were used as its starting cases, but it always improved superpositions obained by MODELLER and STAMP. A program entitled PyMSS is presented that includes three superposition algorithms featuring human interaction., ((c) 2006 Wiley-Liss, Inc.)

We have developed a series of upper undergraduate/graduate lecture and laboratory courses on biotechnological topics to supplement existing biochemical engineering, bioseparations, and biomedical engineering lecture courses. The laboratory courses are based on problem-based learning techniques, featuring two- and three-person teams, journaling, and performance rubrics for guidance and assessment. Participants initially have found them to be difficult, since they had little experience with problem-based learning. To increase enrollment, we are combining the laboratory courses into 2-credit groupings and allowing students to substitute one of them for the second of our 2-credit chemical engineering unit operations laboratory courses.

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., ((c) 2005 Wiley-Liss, Inc.)

Multiple sequence alignment separates members of glycoside hydrolase Family 6 into eight subfamilies: one of mainly actinobacterial endoglucanases (EGs), one of ascomycotal EGs, one of chytridiomycotal EGs and cellobiohydrolases (CBHs), one of actinobacterial and proteobacterial CBHs, one of chytridiomycotal CBHs, two of ascomycotal CBHs, and one of basidiomycotal CBHs. Each also has some proteins of unknown function. Multiple sequence alignment also extends to all of Family 6 the observation that lengths of loops that form the active-site tunnel in CBHs vary among subfamilies, and along with loop conformations, determine enzyme function., (Copyright 2005 Wiley Periodicals, Inc)

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., (Copyright 2005 Wiley-Liss, Inc.)

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., (Copyright 2005 Wiley-Liss, Inc.)

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., (Copyright 2004 Wiley-Liss, Inc.)

Surfactant protein D (SP-D), one of the members of the collectin family of C-type lectins, is an important component of pulmonary innate immunity. SP-D binds carbohydrates in a calcium-dependent manner, but the mechanisms governing its ligand recognition specificity are not well understood. SP-D binds glucose (Glc) stronger than N-acetylglucosamine (GlcNAc). Structural superimposition of hSP-D with mannose- binding protein C (MBP-C) complexed with GlcNAc reveals steric clashes between the ligand and the side chain of Arg343 in hSP-D. To test whether Arg343 contributes to Glc > GlcNAc recognition specificity, we constructed a computational model of Arg343-->Val (R343V) mutant hSP-D based on homology with MBP-C. Automated docking of alpha-Me-Glc and alpha-Me-GlcNAc into wild-type hSP-D and the R343V mutant of hSP-D suggests that Arg343 is critical in determining ligand-binding specificity by sterically prohibiting one binding orientation. To empirically test the docking predictions, an R343V mutant recombinant hSP-D was constructed. Inhibition analysis shows that the R343V mutant binds both Glc and GlcNAc with higher affinity than the wild-type protein and that the R343V mutant binds Glc and GlcNAc equally well. These data demonstrate that Arg343 is critical for hSP-D recognition specificity and plays a key role in defining ligand specificity differences between MBP and SP-D. Additionally, our results suggest that the number of binding orientations contributes to monosaccharide binding affinity.

We present an automated docking protocol specifically optimized to predict the structure and affinity of a protein-carbohydrate complex. A scoring function was developed based on a training set of 30 protein-carbohydrate complexes of known structure and affinity. Combinations of several models for hydrogen bonding, torsional entropy loss, and solvation were tested for their ability to fit the training set data, and the best model was used with AutoDock. The electrostatic empirical coefficient is larger than in a previously obtained model using a training set comprised of various types of protein-ligand complexes, indicating that electrostatic interactions play a more important role in determining the affinity between a carbohydrate and a protein. The differences in the relative weighting of the empirical coefficients in the model yields predicted free energies for the training set with a standard error of 1.403 kcal/mol. The new scoring function was tested on 17 Aspergillus niger glucoamylase inhibitors for which binding energies had been determined experimentally. Free energies of complex formation were predicted with a residual standard error of 1.101 kcal/mol. The new scoring function therefore provides a robust method for predicting free energies of formation and optimal conformations of carbohydrate-protein complexes., (Copyright 2003 Wiley Periodicals, Inc.)