The hemoprotein myoglobin is a model system to study protein dynamics. We used time-resolved serial femtosecond crystallography at an x-ray free-electron laser to resolve the ultrafast structural changes taking place in the carbonmonoxy myoglobin complex upon photolysis of the Fe-CO bond. Structural changes appear throughout the protein within 500 fs with the C-, F- and H-helices moving away from the heme and the E- and A-helices moving toward it. These collective movements are predicted by quantum mechanics/molecular mechanics simulations. Together with the observed oscillations of residues contacting the heme, the calculations support predictions that an immediate collective response of the protein takes place upon ligand dissociation due to coupling of vibrational modes of the heme to global modes of the protein

Using ab initio metadynamics we have computed the conformational free energy landscape of beta-D-glucopyranose as a function of the puckering coordinates. We show that the correspondence between the free energy and the Stoddard's pseudorotational itinerary for the system is rather poor. The number of free energy minima (9) is smaller than the number of ideal structures (13). Moreover, only six minima correspond to a canonical conformation. The structural features, the electronic properties, and the relative stability of the predicted conformers permit the rationalization of the occurrence of distorted sugar conformations in all the available X-ray structures of beta-glucoside hydrolase Michaelis complexes. We show that these enzymes recognize the most stable distorted conformers of the isolated substrate and at the same time the ones better prepared for catalysis in terms of bond elongation/shrinking and charge distribution. This suggests that the factors governing the distortions present in these complexes are largely dictated by the intrinsic properties of a single glucose unit.

AKR1B10 is a human aldo-keto reductase (AKR) found to be elevated in several cancer types and in precancerous lesions. In vitro, AKR1B10 exhibits a much higher retinaldehyde reductase activity than any other human AKR, including AKR1B1 (aldose reductase). We here demonstrate that AKR1B10 also acts as a retinaldehyde reductase in vivo. This activity may be relevant in controlling the first step of retinoic acid synthesis. Up-regulation of AKR1B10, resulting in retinoic acid depletion, may lead to cellular proliferation. Both in vitro and in vivo activities of AKR1B10 were inhibited by tolrestat, an AKR1B1 inhibitor developed for diabetes treatment. The crystal structure of the ternary complex AKR1B10 -NADP ؉ -tolrestat was determined at 1.25-Å resolution. Molecular dynamics models of AKR1B10 and AKR1B1 with retinaldehyde isomers and site-directed mutagenesis show that subtle differences at the entrance of the retinoid-binding site, especially at position 125, are determinant for the all-transretinaldehyde specificity of AKR1B10. Substitutions in the retinaldehyde cyclohexene ring also influence the specificity. These structural features should facilitate the design of specific inhibitors, with potential use in cancer and diabetes treatments.aldo-keto reductases ͉ aldose reductase ͉ protein structure ͉ retinoic acid ͉ tolrestat

Carbohydrate-active enzymes such as glycoside hydrolases (GHs) and glycosyltransferases (GTs) are of growing importance as drug targets. The development of efficient competitive inhibitors and chaperones to treat diseases related to these enzymes requires a detailed knowledge of their mechanisms of action. In recent years, sophisticated first-principles modeling approaches have significantly advanced in our understanding of the catalytic mechanisms of GHs and GTs, not only the molecular details of chemical reactions but also the significant implications that just the conformational dynamics of a sugar ring can have on these mechanisms. Here we provide an overview of the progress that has been made in the past decade, combining molecular dynamics simulations with density functional theory to solve these sweet mysteries of nature.

A quantum leap: By means of quantum mechanics/molecular mechanics metadynamics simulations, a front‐face SNi‐type reaction for glycosyl transfer with retention of the anomeric configuration is shown to be feasible. A short‐lived oxocarbenium‐like species (see picture; O red, P gold, N blue, C black) is identified and provides the complete itinerary of this long sought after molecular mechanism.

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 2 S O / 2,5 B 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 boattype conformation of the glucosyl unit in subsite −1. The overall catalytic mechanism follows a D N *A N -like mechanism and a β-2 S O → 2,5 B [TS] → α-5 S 1 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. Disciplines Biochemical and Biomolecular Engineering | Biological Engineering | Chemical Engineering CommentsPosted with permission from The Journal of Physical Chemistry B, 113, no. 20 (2009) ReceiVed: December 29, 2008; ReVised Manuscript ReceiVed: March 18, 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 -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 distortedB 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 inf...

Inorganic surfaces and nanoparticles can accelerate or inhibit the fibrillation process of proteins and peptides, including the biomedically relevant amyloid β peptide. However, the microscopic mechanisms that determine such an effect are still poorly understood. By means of large-scale, state-of-the-art enhanced sampling molecular dynamics simulations, here we identify an interaction mechanism between the segments 16-22 of the amyloid β peptide, known to be fibrillogenic by itself, and the Au(111) surface in water that leads to the suppression of fiber-like conformations from the peptide conformational ensemble. Moreover, thanks to advanced simulation analysis techniques, we characterize the conformational selection vs. induced fit nature of the gold effect. Our results disclose an inhibition mechanism that is rooted in the details of the microscopic peptide-surface interaction rather than in general phenomena such as peptide sequestration from the solution.

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.