Structural glycobiology of heparin dynamics on the exosite 2 of coagulation cascade proteases: Implications for glycosaminoglycans antithrombotic activity

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.

Polyelectrolytes for Enzyme Immobilization and the Regulation of Their Properties

In this review, we considered aspects related to the application of polyelectrolytes, primarily synthetic polyanions and polycations, to immobilize enzymes and regulate their properties. We mainly focused on the description of works in which polyelectrolytes were used to create complex and unusual systems (self-regulated enzyme-polyelectrolyte complexes, artificial chaperones, polyelectrolyte brushes, layer-by-layer immobilization and others). These works represent the field of "smart polymers", whilst the trivial use of charged polymers as carriers for adsorption or covalent immobilization of proteins is beyond the scope of this short review. In addition, we have included a section on the molecular modeling of interactions between proteins and polyelectrolytes, as modeling the binding of proteins with a strictly defined, and already known, spatial structure, to disordered polymeric molecules has its own unique characteristics.

Evidence of a putative glycosaminoglycan binding site on the glycosylated SARS-CoV-2 spike protein N-terminal domain

SARS-CoV-2 has rapidly spread throughout the world's population since its initial discovery in 2019. The virus infects cells via a glycosylated spike protein located on its surface. The protein primarily binds to the angiotensin-converting enzyme-2 (ACE2) receptor, using glycosaminoglycans (GAGs) as co-receptors. Here, we performed bioinformatics and molecular dynamics simulations of the spike protein to investigate the existence of additional GAG binding sites on the receptor-binding domain (RBD), separate from previously reported heparin-binding sites. A putative GAG binding site in the N-terminal domain (NTD) of the protein was identified, encompassing residues 245-246. We hypothesized that GAGs of a sufficient length might bridge the gap between this site and the PRRARS furin cleavage site, including the mutation S247R. Docking studies using GlycoTorch Vina and subsequent MD simulations of the spike trimer in the presence of dodecasaccharides of the GAGs heparin and heparan sulfate supported this possibility. The heparan sulfate chain bridged the gap, binding the furin cleavage site and S247R. In contrast, the heparin chain bound the furin cleavage site and surrounding glycosylation structures, but not S247R. These findings identify a site in the spike protein that favors heparan sulfate binding that may be particularly pertinent for a better understanding of the recent UK and South African strains. This will also assist in future targeted therapy programs that could include repurposing clinical heparan sulfate mimetics.

Integration of network models and evolutionary analysis into high-throughput modeling of protein dynamics and allosteric regulation: theory, tools and applications

Proteins are dynamical entities that undergo a plethora of conformational changes, accomplishing their biological functions. Molecular dynamics simulation and normal mode analysis methods have become the gold standard for studying protein dynamics, analyzing molecular mechanism and allosteric regulation of biological systems. The enormous amount of the ensemble-based experimental and computational data on protein structure and dynamics has presented a major challenge for the high-throughput modeling of protein regulation and molecular mechanisms. In parallel, bioinformatics and systems biology approaches including genomic analysis, coevolution and network-based modeling have provided an array of powerful tools that complemented and enriched biophysical insights by enabling high-throughput analysis of biological data and dissection of global molecular signatures underlying mechanisms of protein function and interactions in the cellular environment. These developments have provided a powerful interdisciplinary framework for quantifying the relationships between protein dynamics and allosteric regulation, allowing for high-throughput modeling and engineering of molecular mechanisms. Here, we review fundamental advances in protein dynamics, network theory and coevolutionary analysis that have provided foundation for rapidly growing computational tools for modeling of allosteric regulation. We discuss recent developments in these interdisciplinary areas bridging computational biophysics and network biology, focusing on promising applications in allosteric regulations, including the investigation of allosteric communication pathways, protein–DNA/RNA interactions and disease mutations in genomic medicine. We conclude by formulating and discussing future directions and potential challenges facing quantitative computational investigations of allosteric regulatory mechanisms in protein systems.

Protein Interaction with Charged Macromolecules: From Model Polymers to Unfolded Proteins and Post-Translational Modifications

Interaction of proteins with charged macromolecules is involved in many processes in cells. Firstly, there are many naturally occurred charged polymers such as DNA and RNA, polyphosphates, sulfated glycosaminoglycans, etc., as well as pronouncedly charged proteins such as histones or actin. Electrostatic interactions are also important for "generic" proteins, which are not generally considered as polyanions or polycations. Finally, protein behavior can be altered due to post-translational modifications such as phosphorylation, sulfation, and glycation, which change a local charge of the protein region. Herein we review molecular modeling for the investigation of such interactions, from model polyanions and polycations to unfolded proteins. We will show that electrostatic interactions are ubiquitous, and molecular dynamics simulations provide an outstanding opportunity to look inside binding and reveal the contribution of electrostatic interactions. Since a molecular dynamics simulation is only a model, we will comprehensively consider its relationship with the experimental data.

Mass Spectrometry Reveals a Multifaceted Role of Glycosaminoglycan Chains in Factor Xa Inactivation by Antithrombin

Factor Xa (fXa) inhibition by antithrombin (AT) enabled by heparin or heparan sulfate is critical for controlling blood coagulation. AT activation by heparin has been investigated extensively, while interaction of heparin with trapped AT/fXa intermediates has received relatively little attention. We use native electrospray ionization mass spectrometry to study the role of heparin chains of varying length [hexa-, octa-, deca-, and eicosasaccharides (dp6, dp8, dp10, and dp20, respectively)] in AT/fXa complex assembly. Despite being critical promoters of AT/Xa binding, shorter heparin chains are excluded from the final products (trapped intermediates). However, replacement of short heparin segments with dp20 gives rise to a prominent ionic signal of ternary complexes. These species are also observed when the trapped intermediate is initially prepared in the presence of a short oligoheparin (dp6), followed by addition of a longer heparin chain (dp20), indicating that binding of heparin to AT/fXa complexes takes place after the inhibition event. The importance of the heparin chain length for its ability to associate with the trapped intermediate suggests that the binding likely occurs in a bidentate fashion (where two distinct segments of oligoheparin make contacts with the protein components, while the part of the chain separating these two segments is extended into solution to minimize electrostatic repulsion). This model is corroborated by both molecular dynamics simulations with an explicit solvent and ion mobility measurements in the gas phase. The observed post-inhibition binding of heparin to the trapped AT/fXa intermediates hints at the likely role played by heparan sulfate in their catabolism.

Multiscale modeling of glycosaminoglycan structure and dynamics: current methods and challenges

Glycosaminoglycans are long unbranched and complex polysaccharides that are an essential component of mammalian extracellular matrices. Characterization of their molecular structure, dynamics and interactions are essential to understand important biological phenomena in health and disease, and will lead to novel therapeutics and medical devices. However, this has proven to be a challenge experimentally and theoretical techniques are needed to develop new hypotheses, and interpret experiments. This review aims to examine the current theoretical (rather than experimental) methods used by researchers to investigate glycosaminoglycan structure, dynamics and interactions, from the monosaccharide to the macromolecular scale. It will consider techniques such as quantum mechanics, molecular mechanics, molecular dynamics, coarse graining and docking.

Effects of Glycosylation on the Enzymatic Activity and Mechanisms of Proteases

Posttranslational modifications are an important feature of most proteases in higher organisms, such as the conversion of inactive zymogens into active proteases. To date, little information is available on the role of glycosylation and functional implications for secreted proteases. Besides a stabilizing effect and protection against proteolysis, several proteases show a significant influence of glycosylation on the catalytic activity. Glycans can alter the substrate recognition, the specificity and binding affinity, as well as the turnover rates. However, there is currently no known general pattern, since glycosylation can have both stimulating and inhibiting effects on activity. Thus, a comparative analysis of individual cases with sufficient enzyme kinetic and structural data is a first approach to describe mechanistic principles that govern the effects of glycosylation on the function of proteases. The understanding of glycan functions becomes highly significant in proteomic and glycomic studies, which demonstrated that cancer-associated proteases, such as kallikrein-related peptidase 3, exhibit strongly altered glycosylation patterns in pathological cases. Such findings can contribute to a variety of future biomedical applications.

Heparan sulfate proteoglycans mediate factor XIIa binding to the cell surface

Hageman factor (FXIIa) initiates the intrinsic coagulation pathway and triggers the kallikrein-kinin and the complement systems. In addition, it functions as a growth factor by expressing promitogenic activities toward several cell types. FXIIa binds to the cell surface via a number of structurally unrelated surface receptors; however, the underlying mechanisms are not yet fully understood. Here, we demonstrate that FXIIa utilizes cell membrane-bound glycosaminoglycans to interact with the cell surface of human lung fibroblasts (HLF). The combination of enzymatic, inhibitory, and overexpression approaches identified a heparan sulfate (HS) component of proteoglycans as an important determinant of the FXIIa binding capacity of HLF. Moreover, cell-free assays and competition experiments revealed preferential binding of FXIIa to HS and heparin over dextran sulfate, dermatan sulfate, and chondroitin sulfate A and C. Finally, we demonstrate that fibroblasts isolated from the lungs of the patients suffering from idiopathic pulmonary fibrosis (IPF) exhibit enhanced FXIIa binding capacity. Increased sulfation of HS resulting from elevated HS 6-O-sulfotransferase-1 expression in IPF HLF accounted, in part, for this phenomenon. Application of RNA interference technology and inhibitors of intracellular sulfation revealed the cooperative action of cell surface-associated HS and urokinase-type plasminogen activator receptor in the accumulation of FXIIa on the cell surface of IPF HLF. Moreover, FXIIa stimulated IPF HLF migration, which was abrogated by pretreatment of cells with heparinase I. Collectively, our study uncovers a novel role of HS-type glycosaminoglycans in a local accumulation of FXIIa on the cell membrane. The enhanced association of FXIIa with IPF HLF suggests its contribution to fibrogenesis.

Surface modification and endothelialization of biomaterials as potential scaffolds for vascular tissue engineering applications

This review highlights the recent developments of surface modification and endothelialization of biomaterials in vascular tissue engineering applications.