Unravelling the Antifibrinolytic Mechanism of Action of the 1,2,3-Triazole Derivatives.

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.

Control of Substrate Conformation by Hydrogen Bonding in a Retaining ß-Endoglycosidase.

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 C1 conformation to a distorted 1 S3 /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.

Binding of common organic UV-filters to the thyroid hormone transport protein transthyretin using in vitro and in silico studies: Potential implications in health.

Several anthropogenic contaminants have been identified as competing with the thyroid hormone thyroxine (T4) 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 T4 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.

1,2,3-Triazole Derivatives as Novel Antifibrinolytic Drugs.

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 IC50 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.

Structure-function relationships underlying the dual N-acetylmuramic and N-acetylglucosamine specificities of the bacterial peptidoglycan deacetylase PdaC.

Bacillus subtilis PdaC (BsPdaC) is a membrane-bound, multidomain peptidoglycan N-deacetylase acting on N-acetylmuramic acid (MurNAc) residues and conferring lysozyme resistance to modified cell wall peptidoglycans. BsPdaC contains a C-terminal family 4 carbohydrate esterase (CE4) catalytic domain, but unlike other MurNAc deacetylases, BsPdaC 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 BsPdaC CE4 domain at 1.54 Å resolution and analyzed its mode of action on COS substrates. We found that the minimal substrate is GlcNAc3 and that activity increases with the degree of glycan polymerization. COS deacetylation kinetics revealed that BsPdaC 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, BsPdaC shares higher sequence similarity with the peptidoglycan GlcNAc deacetylase SpPgdaA than with other MurNAc deacetylases. Therefore, we used ligand-docking simulations to analyze the dual GlcNAc- and MurNAc-binding specificities of BsPdaC and compared them with those of SpPgdA and BsPdaA, representing peptidoglycan deacetylases highly specific for GlcNAc or MurNAc residues, respectively. BsPdaC 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. BsPdaC contains short loops similar to those in SpPgdA, 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.

Oxazoline or Oxazolinium Ion? The Protonation State and Conformation of the Reaction Intermediate of Chitinase Enzymes Revisited.

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 C1 /4 H5 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.

Substrate Recognition and Specificity of Chitin Deacetylases and Related Family 4 Carbohydrate Esterases.

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.

A Native Ternary Complex Trapped in a Crystal Reveals the Catalytic Mechanism of a Retaining Glycosyltransferase.

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.

Next generation sequencing in nonsyndromic intellectual disability: from a negative molecular karyotype to a possible causative mutation detection.

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.

Structural basis of chitin oligosaccharide deacetylation.

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.

Comparison of two peroxidases with high potential for biotechnology applications - HRP vs. APEX2.

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.

Nucleation process of a fibril precursor in the C-terminal segment of amyloid-ß.

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.

Controlled processivity in glycosyltransferases: A way to expand the enzymatic toolbox.

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.

Multidimensional view of amyloid fibril nucleation in atomistic detail.

Starting from a disordered aggregate, we have simulated the formation of ordered amyloid-like beta structures in a system formed by 18 polyvaline chains in explicit solvent, employing molecular dynamics accelerated by bias-exchange metadynamics. We exploited 8 different collective variables to compute the free energy of hundreds of putative aggregate structures, with variable content of parallel and antiparallel ß-sheets and different packing among the sheets. This allowed characterizing in detail a possible nucleation pathway for the formation of amyloid fibrils: first the system forms a relatively large ordered nucleus of antiparallel ß-sheets, and then a few parallel sheets start appearing. The relevant nucleation process culminates at this point: when a sufficient number of parallel sheets is formed, the free energy starts to decrease toward a new minimum in which this structure is predominant. The complex nucleation pathway we found cannot be described within classical nucleation theory, namely employing a unique simple reaction coordinate like the total content of ß-sheets.

Enzymatic Hydrolysis of Human Milk Oligosaccharides. The Molecular Mechanism of Bifidobacterium Bifidum Lacto-N-biosidase.

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 pK 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.

Integrative structure determination reveals functional global flexibility for an ultra-multimodular arabinanase.

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.

Catalytic itinerary in 1,3-1,4-ß-glucanase unraveled by QM/MM metadynamics. Charge is not yet fully developed at the oxocarbenium ion-like transition state.

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.

Molecular motions in drug design: the coming age of the metadynamics method.

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 conformational free-energy landscape of ß-D-mannopyranose: evidence for a (1)S(5) → B(2,5) → (O)S(2) catalytic itinerary in ß-mannosidases.

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 molecular mechanism of the catalase reaction.

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.