Computational Analysis of Solvent Inclusion in Docking Studies of Protein-Glycosaminoglycan Systems

Solvent Model Benchmark for Molecular Dynamics of Glycosaminoglycans

In computational studies of glycosaminoglycans (GAGs), a group of anionic, periodic linear polysaccharides, so far there has been very little discussion about the role of solvent models in the molecular dynamics simulations of these molecules. Predominantly, the TIP3P water model is commonly used as one of the most popular explicit water models in general. However, there are numerous alternative explicit and implicit water models that are neglected in the computational research of GAGs. Since solvent-mediated interactions are particularly important for GAG dynamic and structural properties, it would be of great interest for the GAG community to establish the solvent model that is suited the best in terms of the quality of theoretically obtained GAG parameters and, at the same time, would be reasonably demanding in terms of computational resources required. In this study, heparin (HP) was simulated using five implicit and six explicit solvent models with the aim to find out how different solvent models influence HP's molecular descriptors in the molecular dynamics simulations. Here, we initiate the search for the most appropriate solvent representation for GAG systems and we hope to encourage other groups to contribute to this highly relevant subject.

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.

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.