The interaction of six low-molecular tissue-clearing agents (1,2 and 1,3-propanediol, ethylene glycol, glycerol, xylitol, sorbitol) with the collagen mimetic peptide (GPH)3 was studied by applying the methods of classical molecular dynamics (GROMACS), molecular docking (AutoDock Vina) and quantum chemistry (PM6 and B3LYP). The spatial configurations of intermolecular complexes were determined and interaction energies calculated. The dependence of the volume occupied by the collagen peptide on the clearing agent concentration in an aqueous solution was calculated. This dependence is not linear, and has a maximum for almost all the agents in the study. The correlations between the optical clearing potential and intermolecular interactions parameters, such as the time of an agent being in a hydrogen-bonded state, and the relative probability of formation of double hydrogen bonds and interaction energies, were determined. Using the correlations determined, we predicted the numeric value of the optical clearing potential of dextrose molecules in rat skin, which correlates with experimental data. A molecular mechanism of tissue optical clearing within the post-diffusion stage is suggested. Graphical abstract The molecular modeling of the interaction between clearing agents and collagen.
Keywords: Hydrogen binding; Modeling; Molecular dynamics; Tissue molecular structure; Tissue optical clearing.