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 α-maltose and three for different forms of methyl α-acarviosinide to provide energies for twisting of α-(1,4) glycosidic bonds. Maps were also made for β-thiocellobiose and for three β-cellobiose conformers having different glycon ring shapes to model distortions of β-(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.

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.

Conformations and inversion pathways leading to racemization of all the tautomers of gossypol, gossypolone, anhydrogossypol, and a diethylamine Schiff's base of gossypol were investigated with MM3(2000). All forms have hindered rotation because of clashes between the methyl carbon atom and oxygen-containing moieties ortho to the bond linking the two naphthalene rings. Inversion energies generally agree with available experimental data. Gossypol preferentially inverts in its dihemiacetal tautomeric form through the cis pathway (where similar groups clash). Gossypolone inverts more easily than gossypol, and preferentially through the trans pathway (where dissimilar groups clash) when one of its outer rings has an enol-keto group and the other has an aldehyde group. Anhydrogossypol racemizes through the cis pathway. The bridge bond and the ortho exo-cyclic bonds in all the structures bend from planarity, and the inner naphthalene rings pucker to accommodate the inversion. For gossypol, the transition is achieved through greater bending of the exo-cyclic bonds (up to 12°) and less distortion of the inner benzyl rings (q≤0.34 Å), (up to 12.7°) . For gossypolone the transition occurs with greater distortion of the inner benzyl rings (q≤0.63 Å) and less out-of-plane bending (up to 8.4°). By isolating individual clashes, their contribution to the overall barrier can be analyzed, as shown for the dialdehyde tautomer of gossypol.

The first crystal structures of a two-domain, prokaryotic glucoamylase were determined to high resolution from the clostridial species Thermoanaerobacterium thermosaccharolyticum with and without acarbose. The N-terminal domain has 18 antiparallel strands arranged in β-sheets of a super-β-sandwich. The C-terminal domain is an (α/α)₆ barrel, lacking the peripheral subdomain of eukaryotic glucoamylases. Interdomain contacts are common to all prokaryotic Family GH15 proteins. Domains similar to those of prokaryotic glucoamylases in maltose phosphorylases (Family GH65) and glycoaminoglycan lyases (Family PL8) suggest evolution from a common ancestor. Eukaryotic glucoamylases may have evolved from prokaryotic glucoamylases by the substitution of the N-terminal domain with the peripheral subdomain and by the addition of a starch-binding domain.

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 β-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 ²S_O/^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 D_N*A_N-like mechanism and a β-²S_O → ^2,5B [TS] → α-⁵S₁ 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.

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.

Surface tensions, critical micelle concentrations (CMCs), contact angles on hydrophobic polyethylene, and foaming characteristics of phosphatidic acids, phosphatidylcholines, phosphatidylethanolamines, and phosphatidylglycerols were measured to determine their suitability as substitutes for traditional surfactants. These phospholipids have fatty acid chains of 5 to 12 carbon atoms, a range over which they are soluble at room temperature. Their surface tensions decrease with increasing concentrations until their CMCs are reached, above which their plateau surface tensions are as low as 21 mN/m, indicating excellent surface activities. In general, plateau surface tensions decrease with increasing chain length within each phospholipid type. The classical relationship for In CMC vs. chain length is followed with slopes typical of anionic surfactants for phosphatidic acids and phosphatidylglycerols and resembling zwitterionic surfactants for phosphatidylcholines and phosphatidylethanolamines, consistent with the charge on the hydrophilic group. The wetting capabilities of aqueous solutions on polyethylene are good and foam heights and stabilities are high, the latter two properties being comparable to traditional anionic (sodium dodecylsulfate) and nonionic (octylphenol polyethoxylate) surfactants. Some anomalies are observed regarding the effect of chain length on wetting and foaming, probably due to the depletion effect. Many phospholipids slowly degrade in aqueous solution. We conclude that short-chain phospholipids exhibit excellent surfactant properties and may be useful in many applications.

It has been difficult to identify the proton donor and nucleophilic assistant/base of endoplasmic reticulum α-(1→2)-mannosidase I, a member of glycoside hydrolase Family 47, which cleaves the glycosidic bond between two α-(1→2)-linked mannosyl residues by the inverting mechanism, trimming Man₉GlcNAc₂ to Man₈GlcNAc₂ 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 α-d-³S₁-mannopyranosyl-(1→2)-α-d-⁴C₁-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).

The automated docking program AutoDock was used to dock all 38 characteristic β-d-mannopyranose ring conformers into the active site of the yeast endoplasmic reticulum α-(1→2)-mannosidase I, a Family 47 glycoside hydrolase that converts Man₉GlcNAc₂ to Man₈GlcNAc₂. 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 β-d-mannopyranose follows the path ³E→¹C₄→¹H₂→B_2,5 before being expelled by the enzyme. The highly conserved second site that strongly binds β-d-mannopyranose-⁴C₁ may exist to prevent competitive inhibition by the product, and is worthy of further investigation.

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 Arg³⁴³ in hSP-D. To test whether Arg³⁴³contributes to Glc > GlcNAc recognition specificity, we constructed a computational model of Arg³⁴³→Val (R343V) mutant hSP-D based on homology with MBP-C. Automated docking of α-Me-Glc and α-Me-GlcNAc into wild-type hSP-D and the R343V mutant of hSP-D suggests that Arg³⁴³ 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 Arg³⁴³ 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.