Advanced Molecular Dynamics Approaches to Model a Tertiary Complex APRIL/TACI with Long Glycosaminoglycans

Molecular architecture and platelet-activating properties of small immune complexes assembled on heparin and platelet factor 4

Heparin-induced thrombocytopenia (HIT) is an adverse reaction to heparin leading to a reduction in circulating platelets with an increased risk of thrombosis. It is precipitated by polymerized immune complexes consisting of pathogenic antibodies that recognize a small chemokine platelet factor 4 (PF4) bound to heparin. Characterization of these immune complexes is extremely challenging due to the enormous structural heterogeneity of such macromolecular assemblies and their constituents. Native mass spectrometry demonstrates that up to three PF4 tetramers can be assembled on a heparin chain, consistent with the molecular modeling studies showing facile polyanion wrapping along the polycationic belt on the PF4 surface. Although these assemblies can accommodate a maximum of only two antibodies, the resulting immune complexes are capable of platelet activation despite their modest size. Taken together, these studies provide further insight into molecular mechanisms of HIT and other immune disorders where anti-PF4 antibodies play a central role.

Molecular architecture and platelet-activating properties of small immune complexes assembled on intact heparin and their possible involvement in heparin-induced thrombocytopenia

Heparin-induced thrombocytopenia (HIT) is an adverse reaction to heparin leading to a reduction in circulating platelets with an increased risk of thrombosis. It is precipitated by polymerized immune complexes consisting of pathogenic antibodies that recognize a small chemokine platelet factor 4 (PF4) bound to heparin, which trigger platelet activation and a hypercoagulable state. Characterization of these immune complexes is extremely challenging due to the enormous structural heterogeneity of such macromolecular assemblies and their constituents (especially heparin). We use native mass spectrometry to characterize small immune complexes formed by PF4, heparin and monoclonal HIT-specific antibodies. Up to three PF4 tetramers can be assembled on a heparin chain, consistent with the results of molecular modeling studies showing facile polyanion wrapping along the polycationic belt on the PF4 surface. Although these assemblies can accommodate a maximum of only two antibodies, the resulting immune complexes are capable of platelet activation despite their modest size. Taken together, these studies provide further insight into molecular mechanisms of HIT and other immune disorders where anti-PF4 antibodies play a central role.

Repulsive Scaling Replica Exchange Molecular Dynamics in Modeling Protein-Glycosaminoglycan Complexes

Glycosaminoglycans are long linear periodic anionic polysaccharides consisting of disaccharide units exhibiting different sulfation patterns forming a highly heterogeneous group of molecules. Due to their flexibility, length, high charge, and periodicity, they are challenging for computational approaches. Despite their biological significance in terms of the important role in various diseases (e.g., Alzheimer, cancer, SARS-CoV-2) and proper cell functioning (e.g., proliferation, maturation), there is a lack of effective molecular docking tools designed specifically for glycosaminoglycans due to their challenging physical-chemical nature. In this chapter we present protocols for the Repulsive Scaling Replica Exchange Molecular Dynamics (RS-REMD) methods to dock glycosaminoglycans with both implicit and explicit solvent models implemented. This novel molecular dynamics-based replica exchange technique should help to elevate our current knowledge on the complexes and interactions between glycosaminoglycans and their protein receptors.

TNFRSF13B in B cell responses to organ transplantation

Antibodies directed against organ transplants are thought to pose the most vexing hurdle to enduring function and survival of the transplants, particularly organ xenotransplants, and accordingly basic and clinical investigation has focused on elucidating the specificity and pathogenicity of graft-specific antibodies. While much has been learned about these matters, far less is known about the B cells producing graft-specific antibodies and why these antibodies appear to injure some grafts but not others. With the goal of addressing those questions, we have investigated the properties of tumor necrosis factor receptor super family-13B (TNFRSF13B), which regulates various aspects of B cell responses. A full understanding of the functions of TNFRSF13B however is hindered by extreme polymorphism and by diversity of interactions of the protein. Nevertheless, TNFRSF13B variants have been found to exert distinct impact on natural and elicited antibody responses and host defense and mutations of TNFRSF13B have been found to influence the propensity for development of antibody-mediated rejection of organ transplants. Because B cell responses potentially limit application of xenotransplantation, understanding how TNFRSF13B diversity and TNFRSF13B variants govern immunity in xenotransplantation could inspire development of novel therapeutics that could in turn accelerate clinical implementation of xenotransplantation.

Explicit solvent repulsive scaling replica exchange molecular dynamics (RS-REMD) in molecular modeling of protein-glycosaminoglycan complexes

Glycosaminoglcyans (GAGs), linear anionic periodic polysaccharides, are crucial for many biologically relevant functions in the extracellular matrix. By interacting with proteins GAGs mediate processes such as cancer development, cell proliferation and the onset of neurodegenerative diseases. Despite this eminent importance of GAGs, they still represent a limited focus for the computational community in comparison to other classes of biomolecules. Therefore, there is a lack of modeling tools designed specifically for docking GAGs. One has to rely on existing docking software developed mostly for small drug molecules substantially differing from GAGs in their basic physico-chemical properties. In this study, we present an updated protocol for docking GAGs based on the Repulsive Scaling Replica Exchange Molecular Dynamics (RS-REMD) that includes explicit solvent description. The use of this water model improved docking performance both in terms of its accuracy and speed. This method represents a significant computational progress in GAG-related research.