Search

A new family of antifibrinolytic drugs has been recently discovered, combining a triazole moiety, an oxadiazolone, and a terminal amine. Two of the molecules of this family have shown activity that is greater than or similar to that of tranexamic acid (TXA), the current antifibrinolytic gold standard, which has been associated with several side effects and whose use is limited in patients with renal impairment. The aim of this work was to thoroughly examine the mechanism of action of the two ideal candidates of the 1,2,3-triazole family and compare them with TXA, to identify an antifibrinolytic alternative active at lower dosages. Specifically, the antifibrinolytic activity of the two compounds ( 1 and 5 ) and TXA was assessed in fibrinolytic isolated systems and in whole blood. Results revealed that despite having an activity pathway comparable to that of TXA, both compounds showed greater activity in blood. These differences could be attributed to a more stable ligand-target binding to the pocket of plasminogen for compounds 1 and 5 , as suggested by molecular dynamic simulations. This work presents further evidence of the antifibrinolytic activity of the two best candidates of the 1,2,3-triazole family and paves the way for incorporating these molecules as new antifibrinolytic therapies.

Peroxidases are essential elements in many biotechnological applications. An especially interesting concept involves split enzymes, where the enzyme is separated into two smaller and inactive proteins that can dimerize into a fully active enzyme. Such split forms were developed for the horseradish peroxidase (HRP) and ascorbate peroxidase (APX) already. Both peroxidases have a high potential for biotechnology applications. In the present study, we performed biophysical comparisons of these two peroxidases and their split analogues. The active site availability is similar for all four structures. The split enzymes are comparable in stability with their native analogues, meaning that they can be used for further biotechnology applications. Also, the tertiary structures of the two peroxidases are similar. However, differences that might help in choosing one system over another for biotechnology applications were noticed. The main difference between the two systems is glycosylation which is not present in the case of APX/sAPEX2, while it has a high impact on the HRP/sHRP stability. Further differences are calcium ions and cysteine bridges that are present only in the case of HRP/sHRP. Finally, computational results identified sAPEX2 as the systems with the smallest structural variations during molecular dynamics simulations showing its dominant stability comparing to other simulated proteins. Taken all together, the sAPEX2 system has a high potential for biotechnological applications due to the lack of glycans and cysteines, as well as due to high stability., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (©2024PublishedbyElsevierB.V.)

Bacterial β-glycosidases are hydrolytic enzymes that depolymerize polysaccharides such as β-cellulose, β-glucans and β-xylans from different sources, offering diverse biomedical and industrial uses. It has been shown that a conformational change of the substrate, from a relaxed 4 C 1 conformation to a distorted 1 S 3 / 1,4 B conformation of the reactive sugar, is necessary for catalysis. However, the molecular determinants that stabilize the substrate's distortion are poorly understood. Here we use quantum mechanics/molecular mechanics (QM/MM)-based molecular dynamics methods to assess the impact of the interaction between the reactive sugar, i. e. the one at subsite -1, and the catalytic nucleophile (a glutamate) on substrate conformation. We show that the hydrogen bond involving the C2 exocyclic group and the nucleophile controls substrate conformation: its presence preserves sugar distortion, whereas its absence (e.g. in an enzyme mutant) knocks it out. We also show that 2-deoxy-2-fluoro derivatives, widely used to trap the reaction intermediates by X-ray crystallography, reproduce the conformation of the hydrolysable substrate at the experimental conditions. These results highlight the importance of the 2-OH⋅⋅⋅nucleophile interaction in substrate recognition and catalysis in endo-glycosidases and can inform mutational campaigns aimed to search for more efficient enzymes., (© 2023 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

Glycosyltransferases (GT) catalyse the biosynthesis of complex carbohydrates which are the most abundant group of molecules in nature. They are involved in several key mechanisms such as cell signalling, biofilm formation, host immune system invasion or cell structure and this in both prokaryotic and eukaryotic cells. As a result, research towards complete enzyme mechanisms is valuable to understand and elucidate specific structure-function relationships in this group of molecules. In a next step this knowledge could be used in GT protein engineering, not only for rational drug design but also for multiple biotechnological production processes, such as the biosynthesis of hyaluronan, cellooligosaccharides or chitooligosaccharides. Generation of these poly- and/or oligosaccharides is possible due to a common feature of several of these GTs: processivity. Enzymatic processivity has the ability to hold on to the growing polymer chain and some of these GTs can even control the number of glycosyl transfers. In a first part, recent advances in understanding the mechanism of various processive enzymes are discussed. To this end, an overview is given of possible engineering strategies for the purpose of new industrial and fundamental applications. In the second part of this review, we focused on specific chain length-controlling mechanisms, i.e., key residues or conserved regions, and this for both eukaryotic and prokaryotic enzymes., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Several anthropogenic contaminants have been identified as competing with the thyroid hormone thyroxine (T 4 ) for binding to transport proteins as transthyretin (TTR). This binding can potentially create toxicity mechanisms posing a threat to human health. Many organic UV filters (UVFs) and paraben preservatives (PBs), widely used in personal care products, are chemicals of emerging concern due to their adverse effects as potential thyroid-disrupting compounds. Recently, organic UVFs have been found in paired maternal and fetal samples and PBs have been detected in placenta, which opens the possibility of the involvement of TTR in the transfer of these chemicals across physiological barriers. We aimed to investigate a discrete set of organic UVFs and PBs to identify novel TTR binders. The binding affinities of target UVFs towards TTR were evaluated using in vitro T 4 competitive binding assays. The ligand-TTR affinities were determined by isothermal titration calorimetry (ITC) and compared with known TTR ligands. In parallel, computational studies were used to predict the 3-D structures of the binding modes of these chemicals to TTR. Some organic UVFs, compounds 2,2',4,4'-tetrahydroxybenzophenone (BP2, Kd = 0.43 μM); 2,4-dihydroxybenzophenone (BP1, Kd = 0.60 μM); 4,4'-dihydroxybenzophenone (4DHB, Kd = 0.83 μM), and 4-hydroxybenzophenone (4HB, Kd = 0.93 μM), were found to display a high affinity to TTR, being BP2 the strongest TTR binder (ΔH = -14.93 Kcal/mol). Finally, a correlation was found between the experimental ITC data and the TTR-ligand docking scores obtained by computational studies. The approach integrating in vitro assays and in silico methods constituted a useful tool to find TTR binders among common organic UVFs. Further studies on the involvement of the transporter protein TTR in assisting the transplacental transfer of these chemicals across physiological barriers and the long-term consequences of prenatal exposure to them should be pursued., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Fibrinolysis is a natural process that ensures blood fluidity through the removal of fibrin deposits. However, excessive fibrinolytic activity can lead to complications in different circumstances, such as general surgery or severe trauma. The current antifibrinolytic drugs in the market, aminocaproic acid (EACA) and tranexamic acid (TXA), require high doses repetitively to maintain their therapeutic effect. These high doses are related to a number of side effects such as headaches, nasal symptoms, or gastrointestinal discomfort and severely limit their use in patients with renal impairment. Therefore, the discovery of novel antifibrinolytics with a higher specificity and lower dosage could vastly improve the applicability of these drugs. Herein, we synthesized a total of ten compounds consisting of a combination of three key moieties: an oxadiazolone, a triazole, and a terminal amine. The IC 50 of each compound was calculated in our clot lysis assays, and the best candidate ( 1 ) provided approximately a 2.5-fold improvement over the current gold standard, TXA. Molecular docking and molecular dynamics were used to perform a structure-activity relationship (SAR) analysis with the lysine binding site in the Kringle 1 domain of plasminogen. This analysis revealed that 1,2,3-triazole was crucial for the activity, enhancing the binding affinity through pi-pi stacking and polar interactions with Tyr72. The results presented in this work open the door to further investigate this new family as potential antifibrinolytic drugs.

AbnA is an extracellular GH43 α-L-arabinanase from Geobacillus stearothermophilus, a key bacterial enzyme in the degradation and utilization of arabinan. We present herein its full-length crystal structure, revealing the only ultra-multimodular architecture and the largest structure to be reported so far within the GH43 family. Additionally, the structure of AbnA appears to contain two domains belonging to new uncharacterized carbohydrate-binding module (CBM) families. Three crystallographic conformational states are determined for AbnA, and this conformational flexibility is thoroughly investigated further using the "integrative structure determination" approach, integrating molecular dynamics, metadynamics, normal mode analysis, small angle X-ray scattering, dynamic light scattering, cross-linking, and kinetic experiments to reveal large functional conformational changes for AbnA, involving up to ~100 Å movement in the relative positions of its domains. The integrative structure determination approach demonstrated here may apply also to the conformational study of other ultra-multimodular proteins of diverse functions and structures., (© 2022. The Author(s).)

Bifidobacterium bifidum lacto- N -biosidase (LnbB) is a critical enzyme for the degradation of human milk oligosaccharides in the gut microbiota of breast-fed infants. Guided by recent crystal structures, we unveil its molecular mechanism of catalysis using QM/MM metadynamics. We show that the oligosaccharide substrate follows 1 S 3 / 1,4 B → [ 4 E ] ‡ → 4 C 1 / 4 H 5 and 4 C 1 / 4 H 5 → [ 4 E / 4 H 5 ] ‡ → 1,4 B conformational itineraries for the two successive reaction steps, with reaction free energy barriers in agreement with experiments. The simulations also identify a critical histidine (His263) that switches between two orientations to modulate the p K a of the acid/base residue, facilitating catalysis. The reaction intermediate of LnbB is best depicted as an oxazolinium ion, with a minor population of neutral oxazoline. The present study sheds light on the processing of oligosaccharides of the early life microbiota and will be useful for the engineering of LnbB and similar glycosidases for biocatalysis., Competing Interests: The authors declare no competing financial interest., (© 2022 The Authors. Published by American Chemical Society.)

Bacillus subtilis PdaC ( Bs PdaC) is a membrane-bound, multidomain peptidoglycan N- deacetylase acting on N -acetylmuramic acid (MurNAc) residues and conferring lysozyme resistance to modified cell wall peptidoglycans. Bs PdaC contains a C-terminal family 4 carbohydrate esterase (CE4) catalytic domain, but unlike other MurNAc deacetylases, Bs PdaC also has GlcNAc deacetylase activity on chitooligosaccharides (COSs), characteristic of chitin deacetylases. To uncover the molecular basis of this dual activity, here we determined the X-ray structure of the Bs PdaC CE4 domain at 1.54 Å resolution and analyzed its mode of action on COS substrates. We found that the minimal substrate is GlcNAc 3 and that activity increases with the degree of glycan polymerization. COS deacetylation kinetics revealed that Bs PdaC operates by a multiple-chain mechanism starting at the internal GlcNAc units and leading to deacetylation of all but the reducing-end GlcNAc residues. Interestingly, Bs PdaC shares higher sequence similarity with the peptidoglycan GlcNAc deacetylase Sp PgdaA than with other MurNAc deacetylases. Therefore, we used ligand-docking simulations to analyze the dual GlcNAc- and MurNAc-binding specificities of Bs PdaC and compared them with those of Sp PgdA and Bs PdaA, representing peptidoglycan deacetylases highly specific for GlcNAc or MurNAc residues, respectively. Bs PdaC retains the conserved Asp-His-His metal-binding triad characteristic of CE4 enzymes acting on GlcNAc residues, differing from MurNAc deacetylases that lack the metal-coordinating Asp residue. Bs PdaC contains short loops similar to those in Sp PgdA, resulting in an open binding cleft that can accommodate polymeric substrates. We propose that PdaC is the first member of a new subclass of peptidoglycan MurNAc deacetylases., (© 2019 Grifoll-Romero et al.)

Because of the lack of cell wall, Micoplasma species require a fine control of membrane fluidity and integrity. mg517 is an essential gene of Mycoplasma genitalium responsible for the biosynthesis of membrane glycoglycerolipids. It encodes for a unique glycosyltransferase (MG517) with processive activity, transferring activated glycosyl donors to either nude diacylglycerol or already glycosylated diacylglycerol. This dual activity, asserted to different enzymes in other species, is sensitive to and regulated by the presence of anionic lipid vesicles in vitro. We present here a computational model of the C-terminus domain of MG517 that complements a previous structural model of the N-terminus domain. By means of sequence analysis, molecular dynamics and metadynamics simulations, we have identified a short α-helix at the apical C-terminus of MG517 with clear amphipathic character. Binding to a membrane model is thermodynamically favored which suggests that this structural element guides the adhesion of MG517 to the cell membrane. We have experimentally verified that truncation of part of this helix causes a substantial reduction of glycoglycerolipids synthesis. The model proposes that MG517 recognizes and binds the diacylglycerol substrate embedded in the membrane by means of this α-helix at the C-terminus together with a previously identified binding pocket at the N-terminus.

The enzymatic hydrolysis of chitin, one of the most abundant carbohydrates in nature, is achieved by chitinases, enzymes of increasing importance in biomedicine and industry. Unlike most retaining glycosidases, family GH18 chitinases follow a substrate-assisted mechanism in which the 2-acetamido group of one N-acetylglucosamine monomer, rather than a basic residue of the enzyme, reacts with the sugar anomeric carbon, forming an intermediate that has been described as an oxazolinium ion. Based on QM/MM metadynamics simulations on chitinase B from Serratia marcescens, we show that the reaction intermediate of GH18 chitinases features instead a neutral oxazoline in a 4 C 1 / 4 H 5 conformation, with an oxazolinium ion being formed on the pathway towards the reaction products. The role of a well-defined hydrogen-bond network that operates around the N-acetyl group, orchestrating catalysis by protonation events, is discussed., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Depolymerization and de- N -acetylation of chitin by chitinases and deacetylases generates a series of derivatives including chitosans and chitooligosaccharides (COS), which are involved in molecular recognition events such as modulation of cell signaling and morphogenesis, immune responses, and host-pathogen interactions. Chitosans and COS are also attractive scaffolds for the development of bionanomaterials for drug/gene delivery and tissue engineering applications. Most of the biological activities associated with COS seem to be largely dependent not only on the degree of polymerization but also on the acetylation pattern, which defines the charge density and distribution of GlcNAc and GlcNH₂ moieties in chitosans and COS. Chitin de- N -acetylases (CDAs) catalyze the hydrolysis of the acetamido group in GlcNAc residues of chitin, chitosan, and COS. The deacetylation patterns are diverse, some CDAs being specific for single positions, others showing multiple attack, processivity or random actions. This review summarizes the current knowledge on substrate specificity of bacterial and fungal CDAs, focusing on the structural and molecular aspects of their modes of action. Understanding the structural determinants of specificity will not only contribute to unravelling structure-function relationships, but also to use and engineer CDAs as biocatalysts for the production of tailor-made chitosans and COS for a growing number of applications., Competing Interests: The authors declare no conflict of interest.

Chitin deacetylases (CDAs) act on chitin polymers and low molecular weight oligomers producing chitosans and chitosan oligosaccharides. Structurally-defined, partially deacetylated chitooligosaccharides produced by enzymatic methods are of current interest as bioactive molecules for a variety of applications. Among Pochonia chlamydosporia (Pc) annotated CDAs, gene pc_2566 was predicted to encode for an extracellular CE4 deacetylase with two CBM18 chitin binding modules. Chitosan formation during nematode egg infection by this nematophagous fungus suggests a role for their CDAs in pathogenicity. The P. chlamydosporia CDA catalytic domain (PcCDA) was expressed in E. coli BL21, recovered from inclusion bodies, and purified by affinity chromatography. It displays deacetylase activity on chitooligosaccharides with a degree of polymerization (DP) larger than 3, generating mono- and di-deacetylated products with a pattern different from those of closely related fungal CDAs. This is the first report of a CDA from a nematophagous fungus. On a DP5 substrate, PcCDA gave a single mono-deacetylated product in the penultimate position from the non-reducing end (ADAAA) which was then transformed into a di-deacetylated product (ADDAA). This novel deacetylation pattern expands our toolbox of specific CDAs for biotechnological applications, and will provide further insights into the determinants of substrate specificity in this family of enzymes.

Carbohydrate esterases family 4 (CE4 enzymes) includes chitin and peptidoglycan deacetylases, acetylxylan esterases, and poly- N -acetylglucosamine deacetylases that act on structural polysaccharides, altering their physicochemical properties, and participating in diverse biological functions. Chitin and peptidoglycan deacetylases are not only involved in cell wall morphogenesis and remodeling in fungi and bacteria, but they are also used by pathogenic microorganisms to evade host defense mechanisms. Likewise, biofilm formation in bacteria requires partial deacetylation of extracellular polysaccharides mediated by poly- N -acetylglucosamine deacetylases. Such biological functions make these enzymes attractive targets for drug design against pathogenic fungi and bacteria. On the other side, acetylxylan esterases deacetylate plant cell wall complex xylans to make them accessible to hydrolases, making them attractive biocatalysts for biomass utilization. CE4 family members are metal-dependent hydrolases. They are highly specific for their particular substrates, and show diverse modes of action, exhibiting either processive, multiple attack, or patterned deacetylation mechanisms. However, the determinants of substrate specificity remain poorly understood. Here, we review the current knowledge on the structure, activity, and specificity of CE4 enzymes, focusing on chitin deacetylases and related enzymes active on N -acetylglucosamine-containing oligo and polysaccharides., Competing Interests: The authors declare no conflict of interest. The founding sponsors had no role in the design of the study, in the writing of the manuscript, and in the decision to publish the results.

Fast genome sequencing offers invaluable opportunities for building updated and improved models of protein sequence evolution. We here show that Single Nucleotide Polymorphisms (SNPs) can be used to build a model capable of predicting the probability of substitution between amino acids in variants of the same protein in different species. The model is based on a substitution matrix inferred from the frequency of codon interchanges observed in a suitably selected subset of human SNPs, and predicts the substitution probabilities observed in alignments between Homo sapiens and related species at 85-100% of sequence identity better than any other approach we are aware of. The model gradually loses its predictive power at lower sequence identity. Our results suggest that SNPs can be employed, together with multiple sequence alignment data, to model protein sequence evolution. The SNP-based substitution matrix developed in this work can be exploited to better align protein sequences of related organisms, to refine the estimate of the evolutionary distance between protein variants from related species in phylogenetic trees and, in perspective, might become a useful tool for population analysis., (Copyright © 2017 by the Genetics Society of America.)

Glycosyltransferases (GTs) play a central role in nature. They catalyze the transfer of a sugar moiety to a broad range of acceptor substrates. GTs are highly selective enzymes, allowing the recognition of subtle structural differences in the sequences and stereochemistry of their sugar and acceptor substrates. We report here a series of structural snapshots of the reaction center of the retaining glucosyl-3-phosphoglycerate synthase (GpgS). During this sequence of events, we visualize how the enzyme guides the substrates into the reaction center where the glycosyl transfer reaction takes place, and unveil the mechanism of product release, involving multiple conformational changes not only in the substrates/products but also in the enzyme. The structural data are further complemented by metadynamics free-energy calculations, revealing how the equilibrium of loop conformations is modulated along these itineraries. The information reported here represent an important contribution for the understanding of GT enzymes at the molecular level., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Glycosyltransferases (GTs) comprise a prominent family of enzymes that play critical roles in a variety of cellular processes, including cell signaling, cell development, and host-pathogen interactions. Glycosyl transfer can proceed with either inversion or retention of the anomeric configuration with respect to the reaction substrates and products. The elucidation of the catalytic mechanism of retaining GTs remains a major challenge. A native ternary complex of a GT in a productive mode for catalysis is reported, that of the retaining glucosyl-3-phosphoglycerate synthase GpgS from M. tuberculosis in the presence of the sugar donor UDP-Glc, the acceptor substrate phosphoglycerate, and the divalent cation cofactor. Through a combination of structural, chemical, enzymatic, molecular dynamics, and quantum-mechanics/molecular-mechanics (QM/MM) calculations, the catalytic mechanism was unraveled, thereby providing a strong experimental support for a front-side substrate-assisted SN i-type reaction., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Hexosaminidases are involved in important biological processes catalyzing the hydrolysis of N-acetyl-hexosaminyl residues in glycosaminoglycans and glycoconjugates. The GH20 enzymes present diverse domain organizations for which we propose two minimal model architectures: Model A containing at least a non-catalytic GH20b domain and the catalytic one (GH20) always accompanied with an extra α-helix (GH20b-GH20-α), and Model B with only the catalytic GH20 domain. The large Bifidobacterium bifidum lacto-N-biosidase was used as a model protein to evaluate the minimal functional unit due to its interest and structural complexity. By expressing different truncated forms of this enzyme, we show that Model A architectures cannot be reduced to Model B. In particular, there are two structural requirements general to GH20 enzymes with Model A architecture. First, the non-catalytic domain GH20b at the N-terminus of the catalytic GH20 domain is required for expression and seems to stabilize it. Second, the substrate-binding cavity at the GH20 domain always involves a remote element provided by a long loop from the catalytic domain itself or, when this loop is short, by an element from another domain of the multidomain structure or from the dimeric partner. Particularly, the lacto-N-biosidase requires GH20b and the lectin-like domain at the N- and C-termini of the catalytic GH20 domain to be fully soluble and functional. The lectin domain provides this remote element to the active site. We demonstrate restoration of activity of the inactive GH20b-GH20-α construct (model A architecture) by a complementation assay with the lectin-like domain. The engineering of minimal functional units of multidomain GH20 enzymes must consider these structural requirements.

Chitin and chitosan oligomers have diverse biological activities with potentially valuable applications in fields like medicine, cosmetics, or agriculture. These properties may depend not only on the degrees of polymerization and acetylation, but also on a specific pattern of acetylation (PA) that cannot be controlled when the oligomers are produced by chemical hydrolysis. To determine the influence of the PA on the biological activities, defined chitosan oligomers in sufficient amounts are needed. Chitosan oligomers with specific PA can be produced by enzymatic deacetylation of chitin oligomers, but the diversity is limited by the low number of chitin deacetylases available. We have produced specific chitosan oligomers which are deacetylated at the first two units starting from the non-reducing end by the combined use of two different chitin deacetylases, namely NodB from Rhizobium sp. GRH2 that deacetylates the first unit and COD from Vibrio cholerae that deacetylates the second unit starting from the non-reducing end. Both chitin deacetylases accept the product of each other resulting in production of chitosan oligomers with a novel and defined PA. When extended to further chitin deacetylases, this approach has the potential to yield a large range of novel chitosan oligomers with a fully defined architecture.

Prions are responsible for a heterogeneous group of fatal neurodegenerative diseases, involving post-translational modifications of the cellular prion protein. Epidemiological studies on Creutzfeldt-Jakob disease, a prototype prion disorder, show a majority of cases being sporadic, while the remaining occurrences are either genetic or iatrogenic. The molecular mechanisms by which PrP(C) is converted into its pathological isoform have not yet been established. While point mutations and seeds trigger the protein to cross the energy barriers, thus causing genetic and infectious transmissible spongiform encephalopathies, respectively, the mechanism responsible for sporadic forms remains unclear. Since prion diseases are protein-misfolding disorders, we investigated prion protein folding and stability as functions of different milieus. Using spectroscopic techniques and atomistic simulations, we dissected the contribution of major structural determinants, also defining the energy landscape of prion protein. In particular, we elucidated (i) the essential role of the octapeptide region in prion protein folding and stability, (ii) the presence of a very enthalpically stable intermediate in prion-susceptible species, and (iii) the role of the disulfide bridge in prion protein folding., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

Cell signaling and other biological activities of chitooligosaccharides (COSs) seem to be dependent not only on the degree of polymerization, but markedly on the specific de-N-acetylation pattern. Chitin de-N-acetylases (CDAs) catalyze the hydrolysis of the acetamido group in GlcNAc residues of chitin, chitosan, and COS. A major challenge is to understand how CDAs specifically define the distribution of GlcNAc and GlcNH2 moieties in the oligomeric chain. We report the crystal structure of the Vibrio cholerae CDA in four relevant states of its catalytic cycle. The two enzyme complexes with chitobiose and chitotriose represent the first 3D structures of a CDA with its natural substrates in a productive mode for catalysis, thereby unraveling an induced-fit mechanism with a significant conformational change of a loop closing the active site. We propose that the deacetylation pattern exhibited by different CDAs is governed by critical loops that shape and differentially block accessible subsites in the binding cleft of CE4 enzymes., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Glycosynthases from more than 16 glycosidase families have been developed for the efficient synthesis of oligosaccharides and glycoconjugates. β-1,3-1,4-Glucan oligo- and polysaccharides with defined sequences can be quantitatively achieved with the glycosynthases derived from Bacillus licheniformis β-1,3-1,4-glucanase. The screening of a nucleophile saturation library of this enzyme yielded the unexpected E134D mutant which has high glycosynthase efficiency (25% higher kcat than the best glycosynthase to date, E134S) but also retains some hydrolase activity (2% relative to the wild-type enzyme). Here, we report the biochemical and structural analyses of this mutant compared to E134S and wild-type enzymes. E134D shows a pH profile of general base catalysis for the glycosynthase activity, with a kinetic pKa (on kcat/KM) assigned to Glu138 of 5.8, whereas the same residue acts as a general acid in the hydrolase activity with the same pKa value. The pKa of Glu138 in the wt enzyme was 7.0, a high value due to the presence of the catalytic nucleophile Glu134 which destabilizes the conjugate base of Glu138. Thus, the pKa of Glu138 drops 1.1 pH units in the mutant relative to the wild-type enzyme meaning that the larger distance between carboxylates in positions 138 and 134 (5.6Å for wt, 7.0Å for E134D) and/or a new hydrogen bonding interaction with a third Asp residue (Asp136) in the mutant reduces the effect of the negatively charged Asp134. In consequence, the pKa of Glu138 has a similar pKa value in the E134D mutant than in the other glycosynthase mutants having a neutral residue in position 134. The behavior of the E134D mutant shows that shortening the side chain of the nucleophile, despite maintaining a carboxylate group, confers glycosynthase activity. Therefore E134D is a transitional hydrolase to glycosynthase mutation., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

The identification of causes underlying intellectual disability (ID) is one of the most demanding challenges for clinical Geneticists and Researchers. Despite molecular diagnostics improvements, the vast majority of patients still remain without genetic diagnosis. Here, we report the results obtained using Whole Exome and Target Sequencing on nine patients affected by isolated ID without pathological copy number variations, which were accurately selected from an initial cohort of 236 patients. Three patterns of inheritance were used to search for: (1) de novo, (2) X-linked, and (3) autosomal recessive variants. In three of the nine proband-parent trios analyzed, we identified and validated two de novo and one X-linked potentially causative mutations located in three ID-related genes. We proposed three genes as ID candidate, carrying one de novo and three X-linked mutations. Overall, this systematic proband-parent trio approach using next generation sequencing could explain a consistent percentage of patients with isolated ID, thus increasing our knowledge on the molecular bases of this disease and opening new perspectives for a better diagnosis, counseling, and treatment., (© 2013 Wiley Periodicals, Inc.)

Glycoglycerolipids are structural components of mycoplasma membranes with a fundamental role in membrane properties and stability. Their biosynthesis is mediated by glycosyltransferases (GT) that catalyze the transfer of glycosyl units from a sugar nucleotide donor to diacylglycerol. The essential function of glycolipid synthases in mycoplasma viability, and the absence of glycoglycerolipids in animal host cells make these GT enzymes a target for drug discovery by designing specific inhibitors. However, rational drug design has been hampered by the lack of structural information for any mycoplasma GT. Most of the annotated GTs in pathogenic mycoplasmas belong to family GT2. We had previously shown that MG517 in Mycoplasma genitalium is a GT-A family GT2 membrane-associated glycolipid synthase. We present here a series of structural models of MG517 obtained by homology modeling following a multiple-template approach. The models have been validated by mutational analysis and refined by long scale molecular dynamics simulations. Based on the models, key structure-function relationships have been identified: The N-terminal GT domain has a GT-A topology that includes a non-conserved variable region involved in acceptor substrate binding. Glu193 is proposed as the catalytic base in the GT mechanism, and Asp40, Tyr126, Tyr169, Ile170 and Tyr218 define the substrates binding site. Mutation Y169F increases the enzyme activity and significantly alters the processivity (or sequential transferase activity) of the enzyme. This is the first structural model of a GT-A glycoglycerolipid synthase and provides preliminary insights into structure and function relationships in this family of enzymes.

By extended atomistic simulations in explicit solvent and bias-exchange metadynamics, we study the aggregation process of 18 chains of the C-terminal segment of amyloid-β, an intrinsically disordered protein involved in Alzheimer's disease and prone to form fibrils. Starting from a disordered aggregate, we are able to observe the formation of an ordered nucleus rich in beta sheets. The rate limiting step in the nucleation pathway involves crossing a barrier of approximately 40 kcal/mol and is associated with the formation of a very specific interdigitation of the side chains belonging to different sheets. This structural pattern is different from the one observed experimentally in a microcrystal of the same system, indicating that the structure of a "nascent" fibril may differ from the one of an "extended" fibril.

Retaining glycoside hydrolases (GHs), key enzymes in the metabolism of polysaccharides and glycoconjugates and common biocatalysts used in chemoenzymatic oligosaccharide synthesis, operate via a double-displacement mechanism with the formation of a glycosyl-enzyme intermediate. However, the degree of oxocarbenium ion character of the reaction transition state and the precise conformational itinerary of the substrate during the reaction, pivotal in the design of efficient inhibitors, remain elusive for many GHs. By means of QM/MM metadynamics, we unravel the catalytic itinerary of 1,3-1,4-β-glucanase, one of the most active GHs, belonging to family 16. We show that, in the Michaelis complex, the enzyme environment restricts the conformational motion of the substrate to stabilize a (1,4)B/(1)S(3) conformation of the saccharide ring at the -1 subsite, confirming that this distortion preactivates the substrate for catalysis. The metadynamics simulation of the enzymatic reaction captures the complete conformational itinerary of the substrate during the glycosylation reaction ((1,4)B/(1)S(3) -(4)E/(4)H(3) - (4)C(1)) and shows that the transition state is not the point of maximum charge development at the anomeric carbon. The overall catalytic mechanism is of dissociative type, and proton transfer to the glycosidic oxygen is a late event, clarifying previous kinetic studies of this enzyme., (© 2011 American Chemical Society)

Metadynamics is emerging as a useful free energy method in physics, chemistry and biology. Recently, it has been applied also to investigate ligand binding to biomolecules of pharmacological interest. Here, after introducing the basic idea of the method, we review applications to challenging targets for pharmaceutical intervention. We show that this methodology, especially when combined with a variety of other computational approaches such as molecular docking and/or molecular dynamics simulation, may be useful to predict structure and energetics of ligand/target complexes even when the targets lack a deep binding cavity, such as DNA and proteins undergoing fibrillation in neurodegenerative diseases. Furthermore, the method allows investigating the routes of molecular recognition and the associated binding energy profiles, providing a molecular interpretation to experimental data.

The mechanism of glycosidic bond cleavage by glycosidases involves substrate ring distortions in the Michaelis complex that favor catalysis. Retaining β-mannosidases bind the substrate in a (1)S(5) conformation, and recent experiments have proposed an unusual substrate conformational pathway ((1)S(5) → B(2,5) → (O)S(2)) for the hydrolysis reaction. By means of Car-Parrinello metadynamics simulations, we have obtained the conformational free-energy surface (FES) of a β-d-mannopyranose molecule associated with the ideal Stoddart conformational diagram. We have found that (1)S(5) is among the most stable conformers and simultaneously is the most preactivated conformation in terms of elongation/shortening of the C1-O1/C1-O5 bonds, C1-O1 orientation, and charge development at the anomeric carbon. Analysis of the computed FES gives support to the proposed (1)S(5) → B(2,5) → (O)S(2) catalytic itinerary, showing that the degree of preactivation of the substrate in glycoside hydrolases (GHs) is related to the properties of an isolated sugar ring. We introduce a simple preactivation index integrating several structural, electronic, and energetic properties that can be used to predict the conformation of the substrate in the Michaelis complex of any GH.

The key step in the catalytic cycle of methionine synthase (MetH) is the transfer of a methyl group from the methylcobalamin (MeCbl) cofactor to homocysteine (Hcy). This mechanism has been traditionally viewed as an S(N)2-type reaction, but a different mechanism based on one-electron reduction of the cofactor (reductive cleavage) has been recently proposed. In this work, we analyze whether this mechanism is plausible from a theoretical point of view. By means of a combination of gas-phase as well as hybrid QM/MM calculations, we show that cleavage of the Co-C bond in a MeCbl···Hcy complex (Hcy = methylthiolate substrate (Me-S(-)), a structural mimic of deprotonated homocysteine) proceeds via a [Co(III)(corrin(*-))]-Me···*S-Me diradical configuration, involving electron transfer (ET) from a π*(corrin)-type state to a σ*(Co-C) one, and the methyl transfer displays an energy barrier ≤8.5 kcal/mol. This value is comparable to the one previously computed for the alternative S(N)2 reaction pathway (10.5 kcal/mol). However, the ET-based reductive cleavage pathway does not impose specific geometrical and distance constraints with respect to substrate and cofactor, as does the S(N)2 pathway. This might be advantageous from the enzymatic point of view because in that case, a methyl group can be transferred efficiently at longer distances.

Humans' bitter taste perception is mediated by the hTAS2R subfamily of the G protein-coupled membrane receptors (GPCRs). Structural information on these receptors is currently limited. Here we identify residues involved in the binding of phenylthiocarbamide (PTC) and in receptor activation in one of the most widely studied hTAS2Rs (hTAS2R38) by means of structural bioinformatics and molecular docking. The predictions are validated by site-directed mutagenesis experiments that involve specific residues located in the putative binding site and trans-membrane (TM) helices 6 and 7 putatively involved in receptor activation. Based on our measurements, we suggest that (i) residue N103 participates actively in PTC binding, in line with previous computational studies. (ii) W99, M100 and S259 contribute to define the size and shape of the binding cavity. (iii) W99 and M100, along with F255 and V296, play a key role for receptor activation, providing insights on bitter taste receptor activation not emerging from the previously reported computational models.

Prion diseases are fatal neurodegenerative and infectious disorders for which effective pharmacological tools are not yet available. This unmet challenge and the recently proposed interplay between prion diseases and Alzheimer's have led to a more urgent demand for new antiprion agents. Herein, we report the identification of a novel bifunctional diketopiperazine (DKP) derivative 1 d, which exhibits activity in the low micromolar range against prion replication in ScGT1 cells, while showing low cytotoxicity. Supported by properly addressed molecular modeling studies, we hypothesized that a planar conformation is the major determinant for activity in this class of compounds. Moreover, studies aimed at assessing the mechanism-of-action at the molecular level showed that 1 d might interact directly with recombinant prion protein (recPrP) to prevent its conversion to the pathogenic misfolded prion protein (PrP(Sc))-like form. This investigation suggests that DKP based antiprion compounds can serve as a promising lead scaffold in developing new drugs to combat prion diseases.

Molecular docking of ligands targeting proteins undergoing fibrillization in neurodegenerative diseases is difficult because of the lack of deep binding sites. Here we extend standard docking methods with free energy simulations in explicit solvent to address this issue in the context of the prion protein surface. We focus on a specific ligand (2-pyrrolidin-1-yl-N-[4-[4-(2-pyrrolidin-1-yl-acetylamino)-benzyl]-phenyl]-acetamide), which binds to the structured part of the protein as shown by NMR (Kuwata, K. et al. Proc Natl Acad Sci U.S.A. 2007, 104, 11921-11926). The calculated free energy of dissociation (7.8 ± 0.9 kcal/mol) is in good agreement with the value derived by the experimental dissociation constant (Kd = 3.9 μM, corresponding to ΔG(0) = -7.5 kcal/mol). Several binding poses are predicted, including the one reported previously. Our prediction is fully consistent with the presence of multiple binding sites, emerging from NMR measurements. Our molecular simulation-based approach emerges, therefore, as a useful tool to predict poses and affinities of ligand binding to protein surfaces.

Catalases are ubiquitous enzymes that prevent cell oxidative damage by degrading hydrogen peroxide to water and oxygen (2H(2)O(2) --> 2 H(2)O + O(2)) with high efficiency. The enzyme is first oxidized to a high-valent iron intermediate, known as Compound I (Cpd I) which, in contrast to other hydroperoxidases, is reduced back to the resting state by further reacting with H(2)O(2). By means of hybrid QM/MM Car-Parrinello metadynamics simulations, we have investigated the mechanism of the reduction of Compound I by H(2)O(2) in Helicobacter pylori catalase (HPC) and Penicillium vitale catalase (PVC). We found that the Cpd I-H(2)O(2) complex evolves to a Cpd II-like species through the transfer of a hydrogen atom from the peroxide to the oxoferryl unit. To complete the reaction, two mechanisms may be operative: a His-mediated (Fita-Rossmann) mechanism, which involves the distal His as an acid-base catalyst mediating the transfer of a proton (associated with an electron transfer), and a direct mechanism, in which a hydrogen atom transfer occurs. Independently of the mechanism, the reaction proceeds by two one-electron transfers rather than one two-electron transfer, as has long been the lore. The calculations provide a detailed view of the atomic and electronic reorganizations during the reaction, and highlight the key role of the distal residues to assist the reaction. Additional calculations on the in silico HPC His56Gly mutant and gas-phase models provide clues to understand the requirements for the reaction to proceed with low barriers.

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.

The structure and dynamics of the enzyme-substrate complex of Bacillus 1,3-1,4-beta-glucanase, one of the most active glycoside hydrolases, is investigated by means of Car-Parrinello molecular dynamics simulations (CPMD) combined with force field molecular dynamics (QM/MM CPMD). It is found that the substrate sugar ring located at the -1 subsite adopts a distorted 1S3 skew-boat conformation upon binding to the enzyme. With respect to the undistorted 4C1 chair conformation, the 1S3 skew-boat conformation is characterized by: (a) an increase of charge at the anomeric carbon (C1), (b) an increase of the distance between C1 and the leaving group, and (c) a decrease of the intraring O5-C1 distance. Therefore, our results clearly show that the distorted conformation resembles both structurally and electronically the transition state of the reaction in which the substrate acquires oxocarbenium ion character, and the glycosidic bond is partially broken. Together with analysis of the substrate conformational dynamics, it is concluded that the main determinants of substrate distortion have a structural origin. To fit into the binding pocket, it is necessary that the aglycon leaving group is oriented toward the beta region, and the skew-boat conformation naturally fulfills this premise. Only when the aglycon is removed from the calculation the substrate recovers the all-chair conformation, in agreement with the recent determination of the enzyme product structure. The QM/MM protocol developed here is able to predict the conformational distortion of substrate binding in glycoside hydrolases because it accounts for polarization and charge reorganization at the -1 sugar ring. It thus provides a powerful tool to model E.S complexes for which experimental information is not yet available.

The properties of the Co-C bond in methylcobalamin (MeCbl) are analyzed by means of first-principles molecular dynamics. The optimized structure is in very good agreement with experiments, reproducing the bent-up deformation of the corrin ring as well as the metal-ligand bond distances. The analysis of the binding energies, bond orders, and vibrational stretching frequencies shows that the axial base slightly weakens the Co-C bond (by 4%), while the alkyl ligand substantially reinforces the Co-axial base bond (by 90%). These findings support several experiments and provide insight into the conversion between the base-on and base-off forms of the MeCbl cofactor.