The mechanism of high affinity pentasaccharide binding to antithrombin, insights from Gaussian accelerated molecular dynamics simulations

Glycosaminoglycan microarrays for studying glycosaminoglycan-protein systems

More than 3000 proteins are now known to bind to glycosaminoglycans (GAGs). Yet, GAG-protein systems are rather poorly understood in terms of selectivity of recognition, molecular mechanism of action, and translational promise. High-throughput screening (HTS) technologies are critically needed for studying GAG biology and developing GAG-based therapeutics. Microarrays, developed within the past two decades, have now improved to the point of being the preferred tool in the HTS of biomolecules. GAG microarrays, in which GAG sequences are immobilized on slides, while similar to other microarrays, have their own sets of challenges and considerations. GAG microarrays are rapidly becoming the first choice in studying GAG-protein systems. Here, we review different modalities and applications of GAG microarrays presented to date. We discuss advantages and disadvantages of this technology, explain covalent and non-covalent immobilization strategies using different chemically reactive groups, and present various assay formats for qualitative and quantitative interpretations, including selectivity screening, binding affinity studies, competitive binding studies etc. We also highlight recent advances in implementing this technology, cataloging of data, and project its future promise. Overall, the technology of GAG microarray exhibits enormous potential of evolving into more than a mere screening tool for studying GAG - protein systems.

Distinct antithrombin activation modes for fondaparinux and natural heparins detected using millisecond hydrogen deuterium exchange and collision induced unfolding

Low molecular weight heparin and synthetic mimetics such as fondaparinux show different binding kinetics, protease specificity, and clinical effects. A combination of allosteric and template-mediated bridging mechanisms have been proposed to explain the differences in rate acceleration and specificity. The difficulty in working with heterogeneous heparin species has rendered a crystallographic interpretation of the differences in antithrombin activation between mimetics and natural heparin inaccessible. In this study, we examine the allosteric changes in antithrombin caused by binding fondaparinux, enoxaparin and depolymerized natural heparins using millisecond hydrogen deuterium exchange mass spectrometry (TRESI-HDX MS) and relate these conformational changes to complex stability in the gas phase using collision induced unfolding (CIU). This exploration reveals that in addition to the dynamic changes caused by fondaparinux, long chain heparins reduce structural flexibility proximal to Arg393, the cleavable residue in the reactive centre loop of the protein. These local changes in protein dynamics are associated with an increase in overall complex stability that increases with heparin chain length. Ultimately, these results shed light on the molecular mechanisms underlying differences in activity and specificity between heparin mimetics and natural heparins.

Molecular Mechanisms of the Impaired Heparin Pentasaccharide Interactions in 10 Antithrombin Heparin Binding Site Mutants Revealed by Enhanced Sampling Molecular Dynamics

Antithrombin (AT) is a critical regulator of the coagulation cascade by inhibiting multiple coagulation factors including thrombin and FXa. Binding of heparinoids to this serpin enhances the inhibition considerably. Mutations located in the heparin binding site of AT result in thrombophilia in affected individuals. Our aim was to study 10 antithrombin mutations known to affect their heparin binding in a heparin pentasaccharide bound state using two molecular dynamics (MD) based methods providing enhanced sampling, GaMD and LiGaMD2. The latter provides an additional boost to the ligand and the most important binding site residues. From our GaMD simulations we were able to identify four variants (three affecting amino acid Arg47 and one affecting Lys114) that have a particularly large effect on binding. The additional acceleration provided by LiGaMD2 allowed us to study the consequences of several other mutants including those affecting Arg13 and Arg129. We were able to identify several conformational types by cluster analysis. Analysis of the simulation trajectories revealed the causes of the impaired pentasaccharide binding including pentasaccharide subunit conformational changes and altered allosteric pathways in the AT protein. Our results provide insights into the effects of AT mutations interfering with heparin binding at an atomic level and can facilitate the design or interpretation of in vitro experiments.

Differential Solvent DEEP-STD NMR and MD Simulations Enable the Determinants of the Molecular Recognition of Heparin Oligosaccharides by Antithrombin to Be Disentangled

The interaction of heparin with antithrombin (AT) involves a specific sequence corresponding to the pentasaccharide GlcNAc/NS6S-GlcA-GlcNS3S6S-IdoA2S-GlcNS6S (AGA*IA). Recent studies have revealed that two AGA*IA-containing hexasaccharides, which differ in the sulfation degree of the iduronic acid unit, exhibit similar binding to AT, albeit with different affinities. However, the lack of experimental data concerning the molecular contacts between these ligands and the amino acids within the protein-binding site prevents a detailed description of the complexes. Differential epitope mapping (DEEP)-STD NMR, in combination with MD simulations, enables the experimental observation and comparison of two heparin pentasaccharides interacting with AT, revealing slightly different bound orientations and distinct affinities of both glycans for AT. We demonstrate the effectiveness of the differential solvent DEEP-STD NMR approach in determining the presence of polar residues in the recognition sites of glycosaminoglycan-binding proteins.

The Interaction of Factor Xa and IXa with Non-Activated Antithrombin in Michaelis Complex: Insights from Enhanced-Sampling Molecular Dynamics Simulations

The interaction between coagulation factors Xa and IXa and the activated state of their inhibitor, antithrombin (AT),have been investigated using X-ray diffraction studies. However, only mutagenesis data are available for non-activated AT. Our aim was to propose a model based on docking and advanced-sampling molecular dynamics simulations that can reveal the conformational behavior of the systems when AT is not binding a pentasaccharide. We built the initial structure for non-activated AT-FXa and AT-FIXa complexes using HADDOCK 2.4. The conformational behavior was studied using Gaussian accelerated molecular dynamics simulations. In addition to the docked complexes, two systems based on the X-ray structures were also simulated, with and without the ligand. The simulations revealed large variability in conformation for both factors. In the docking-based complex of AT-FIXa, conformations with stable Arg150-AT interactions can exist for longer time periods but the system also has a higher tendency for reaching states with very limited interaction with the "exosite" of AT. By comparing simulations with or without the pentasaccharide, we were able to gain insights into the effects of conformational activation on the Michaelis complexes. RMSF analysis and correlation calculations for the alpha-carbon atoms revealed important details of the allosteric mechanisms. Our simulations provide atomistic models for better understanding the conformational activation mechanism of AT against its target factors.

Identification and characterization of two SERPINC1 mutations causing congenital antithrombin deficiency

BACKGROUND

Antithrombin (AT) is the main physiological anticoagulant involved in hemostasis. Hereditary AT deficiency is a rare autosomal dominant thrombotic disease mainly caused by mutations in SERPINC1, which was usually manifested as venous thrombosis and pulmonary embolism. In this study, we analyzed the clinical characteristics and screened for mutant genes in two pedigrees with hereditary AT deficiency, and the functional effects of the pathogenic mutations were evaluated.

METHODS

Candidate gene variants were analyzed by next-generation sequencing to screen pathogenic mutations in probands, followed by segregation analysis in families by Sanger sequencing. Mutant and wild-type plasmids were constructed and transfected into HEK293T cells to observe protein expression and cellular localization of SERPINC1. The structure and function of the mutations were analyzed by bioinformatic analyses.

RESULTS

The proband of pedigree A with AT deficiency carried a heterozygous frameshift mutation c.1377delC (p.Asn460Thrfs*20) in SERPINC1 (NM000488.3), a 1377C base deletion in exon 7 resulting in a backward shift of the open reading frame, with termination after translation of 20 residues, and a different residue sequence translated after the frameshift. Bioinformatics analysis suggests that the missing amino acid sequence caused by the frameshift mutation might disrupt the disulfide bond between Cys279 and Cys462 and affect the structural function of the protein. This newly discovered variant is not currently included in the ClinVar and HGMD databases. p.Arg229* resulted in a premature stop codon in exon 4, and bioinformatics analysis suggests that the truncated protein structure lost its domain of interaction with factor IX (Ala414 site) after the deletion of nonsense mutations. However, considering the AT truncation protein resulting from the p.Arg229* variant loss a great proportion of the molecule, we speculate the variant may affect two functional domains HBS and RCL and lack of the corresponding function. The thrombophilia and decreased-AT-activity phenotypes of the two pedigrees were separated from their genetic variants. After lentiviral plasmid transfection into HEK293T cells, the expression level of AT protein decreased in the constructed c.1377delC mutant cells compared to that in the wild-type, which was not only reduced in c.685C > T mutant cells but also showed a significant band at 35 kDa, suggesting a truncated protein. Immunofluorescence localization showed no significant differences in protein localization before and after the mutation.

CONCLUSIONS

The p.Asn460Thrfs*20 and p.Arg229* variants of SERPINC1 were responsible for the two hereditary AT deficiency pedigrees, which led to AT deficiency by different mechanisms. The p.Asn460Thrfs*20 variant is reported for the first time.

Gaussian accelerated molecular dynamics (GaMD): principles and applications

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

Conformational Analysis of Heparin-Analogue Pentasaccharides by Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulations

Elucidation and improvement of the blood coagulant properties of heparin are the focus of intense research. In this study, we performed conformational analysis using nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations on the heparin pentasaccharide analogue idraparinux, its disulfonatomethyl analogue, which features a slightly improved blood coagulation property, and a trisulfonatomethyl analogue, in which the activity has been totally abolished. As the ring conformation of the G subunit has been suggested as a major determinant of the biological properties, we analyzed the sugar ring conformations and dynamics of the interglycosidic linkages. We found that the conformation of the G ring is dominated by the ²S_O skewed boat next to the ¹C₄ chair in all three derivatives. Both the thermodynamics and the kinetics of the conformational states were found to be highly similar in the three derivatives. Molecular kinetic analysis showed that the ²S_O skewed boat state of the G ring is equally favorable in the three analogues, resulting in similar ²S_O populations. Also, the transition kinetics from the ¹C₄ chair to the ²S_O skewed boat was found to be comparable in the derivatives, which indicates a similar energy barrier between the two states of the G subunit. We also identified a slower conformational transition between the dominant ⁴C₁ chair and the boat conformations on the E subunit. Both G and E ring flips are also accompanied by changes along the interglycosidic linkages, which take place highly synchronously with the ring flips. These findings indicate that conformational plasticity of the G ring and the dominance of the ²S_O skewed boat populations do not necessarily warrant the biological activity of the derivatives and hence the impact of other factors also needs to be considered.

Pharmacological and clinical application of heparin progress: An essential drug for modern medicine

Heparin is the earliest and most widely used anticoagulant and antithrombotic drug that is still used in a variety of clinical indications. Since it was discovered in 1916, after more than a century of repeated exploration, heparin has not been replaced by other drugs, but a great progress has been made in its basic research and clinical application. Besides anticoagulant and antithrombotic effects, heparin also has antitumor, anti-inflammatory, antiviral, and other pharmacological activities. It is widely used clinically in cardiovascular and cerebrovascular diseases, lung diseases, kidney diseases, cancer, etc., as the first anticoagulant medicine in COVID-19 exerts anticoagulant, anti-inflammatory and antiviral effects. At the same time, however, it also leads to a lot of adverse reactions, such as bleeding, thrombocytopenia, elevated transaminase, allergic reactions, and others. This article comprehensively reviews the modern research progress of heparin compounds; discusses the structure, preparation, and adverse reactions of heparin; emphasizes the pharmacological activity and clinical application of heparin; reveals the possible mechanism of the therapeutic effect of heparin in related clinical applications; provides evidence support for the clinical application of heparin; and hints on the significance of exploring the wider application fields of heparin.

Investigation of the Differences in Antithrombin to Heparin Binding among Antithrombin Budapest 3, Basel, and Padua Mutations by Biochemical and In Silico Methods

Antithrombin (AT) is a serine protease inhibitor, its activity is highly accelerated by heparin. Mutations at the heparin-binding region lead to functional defect, type II heparin-binding site (IIHBS) AT deficiency. The aim of this study was to investigate and compare the molecular background of AT Budapest 3 (p.Leu131Phe, ATBp3), AT Basel (p.Pro73Leu), and AT Padua (p.Arg79His) mutations. Advanced in silico methods and heparin-binding studies of recombinant AT proteins using surface plasmon resonance method were used. Crossed immunoelectrophoresis and Differential Scanning Fluorimetry (NanoDSF) were performed in plasma samples. Heparin affinity of AT Padua was the lowest (KD = 1.08 × 10^-6 M) and had the most severe consequences affecting the allosteric pathways of activation, moreover significant destabilizing effects on AT were also observed. KD values for AT Basel, ATBp3 and wild-type AT were 7.64 × 10^-7 M, 2.15 × 10^-8 M and 6.4 × 10^-10 M, respectively. Heparin-binding of AT Basel was slower, however once the complex was formed the mutation had only minor effect on the secondary and tertiary structures. Allosteric activation of ATBp3 was altered, moreover decreased thermostability in ATBp3 homozygous plasma and increased fluctuations in multiple regions of ATBp3 were observed by in silico methods suggesting the presence of a quantitative component in the pathogenicity of this mutation due to molecular instability.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

Conformational flexibility correlates with glucose tolerance for point mutations in β-glucosidases - a computational study

β-glucosidases (EC 3.2.1.21) have been described as essential to second-generation biofuel production. They act in the last step of the lignocellulosic saccharification, cleaving the β - 1,4 glycosidic bonds in cellobiose to produce two molecules of glucose. However, β-glucosidases have been described as strongly inhibited by glucose, causing an increment of cellobiose concentration. Also, cellobiose is an inhibitor of other enzymes used in this process, such as exoglucanases and endoglucanases. Hence, the engineering of thermostable and glucose-tolerant β-glucosidases has been targeted by many studies. In this study, we performed high sampling accelerated molecular dynamics for a wild glucose-tolerant GH1 β-glucosidase (Bgl1A), a wild non-tolerant (Bgl1B), and a set of glucose-tolerant Bgl1B's mutants: V302F, N301Q/V302F, F172I, V227M, G246S, T299S, and H228T. Our results suggest that point mutations promissory to induce glucose tolerance trend to enhance the mobility of the flexible loops around the active site. Mutations affected B and C loops regions, and an αβ-hairpin motif between them. Conformational clusters and free energy landscape profiles suggest that the mobility acquired by mutants allows a higher closure of the substrate channel. This closure is compatible with a higher impedance for glucose entrance and stimulus of its withdrawal. Based on mutants' structural analyses, we inferred that both the direct stereochemical effect on the glucose path and the changes in the mobility affect glucose tolerance. We hope these results be useful for the rational design of glucose-tolerant and industrially promising enzymes.Communicated by Ramaswamy H. Sarma.