Systematic probing of an atomic charge set of sialic acid disaccharides for the rational molecular modeling of avian influenza virus based on molecular dynamics simulations

Density functional theory based probe of the affinity interaction of saccharide ligands with extra-cellular sialic acid residues

Changes in glycosylation pattern leads to malignant transformations among the cells. In combination with upregulated actions of sialyltransferases, it ultimately leads to differential expression of sialic acid (SA) at cell surface. Given its negative charge and localization to extracellular domain, SA has been exploited for the development of targeted theranostics using approaches, such as, cationization and appending recognition saccharides on carrier surface. In this study, we have performed quantum mechanical calculations based on density functional theory (DFT) to study the interaction of saccharides with extracellular SA. Gradient-corrected DFT with the three parameter function (B3) was utilized for the calculation of Lee-Yang-Parr (LYP) correlation function. Atomic charge, vibrational frequencies and energy of the optimized structures were calculated through B3LYP. Our calculations demonstrate a stronger galactose-sialic acid interaction at tumour-relevant low pH and hyperthermic condition. These results support the application of pH responsive delivery vehicles and targeted hyperthermic chemotherapy for eradicating solid tumour deposits. These studies, conducted a priori, can guide the formulation scientists over appropriate choice of ligands and their applications in the design of 'smart' theranostic tools.

Enhanced conformational sampling of carbohydrates by Hamiltonian replica-exchange simulation

Knowledge of the structure and conformational flexibility of carbohydrates in an aqueous solvent is important to improving our understanding of how carbohydrates function in biological systems. In this study, we extend a variant of the Hamiltonian replica-exchange molecular dynamics (MD) simulation to improve the conformational sampling of saccharides in an explicit solvent. During the simulations, a biasing potential along the glycosidic-dihedral linkage between the saccharide monomer units in an oligomer is applied at various levels along the replica runs to enable effective transitions between various conformations. One reference replica runs under the control of the original force field. The method was tested on disaccharide structures and further validated on biologically relevant blood group B, Lewis X and Lewis A trisaccharides. The biasing potential-based replica-exchange molecular dynamics (BP-REMD) method provided a significantly improved sampling of relevant conformational states compared with standard continuous MD simulations, with modest computational costs. Thus, the proposed BP-REMD approach adds a new dimension to existing carbohydrate conformational sampling approaches by enhancing conformational sampling in the presence of solvent molecules explicitly at relatively low computational cost.