Network Pharmacology, Molecular Docking and Molecular Dynamics to Explore the Potential Immunomodulatory Mechanisms of Deer Antler

Abstract

The use of deer antlers dates back thousands of years in Chinese history. Deer antlers have antitumor, anti-inflammatory, and immunomodulatory properties and can be used in treating neurological diseases. However, only a few studies have reported the immunomodulatory mechanism of deer antler active compounds. Using network pharmacology, molecular docking, and molecular dynamics simulation techniques, we analyzed the underlying mechanism by which deer antlers regulate the immune response. We identified 4 substances and 130 core targets that may play immunomodulatory roles, and the beneficial and non-beneficial effects in the process of immune regulation were analyzed. The targets were enriched in pathways related to cancer, human cytomegalovirus infection, the PI3K-Akt signaling pathway, human T cell leukemia virus 1 infection, and lipids and atherosclerosis. Molecular docking showed that AKT1, MAPK3, and SRC have good binding activity with 17 beta estradiol and estrone. Additionally, the molecular dynamics simulation of the molecular docking result using GROMACS software (version: 2021.2) was performed and we found that the AKT1–estrone complex, 17 beta estradiol–AKT1 complex, estrone–MAPK3 complex, and 17 beta estradiol–MAPK3 complex had relatively good binding stability. Our research sheds light on the immunomodulatory mechanism of deer antlers and provides a theoretical foundation for further exploration of their active compounds.

Elucidating the anti-cancer potential of Cinnamomum tamala essential oil against non-small cell lung cancer: A multifaceted approach involving GC-MS profiling, network pharmacology, and molecular dynamics simulations

“…This interaction was attributed to hydrophobic interaction with ILE48 and LEU173 residues of MAPK3. Previously 17-β estradiol and deoxysappanone B 7,3′-dimethyl ether acetate formed hydrophobic interaction with LEU173 residue of MAPK3 [ 57 ]. Similarly, in the complex involving kaempferol 3-rutinoside-4′-glucoside and MAPK3, there was also evidence of hydrophobic interaction with the ILE48 residue of MAPK3 [ 58 ].…”

Exploring the effect and mechanism of Aloin A against cancer cachexia-induced muscle atrophy via network pharmacology, molecular docking, molecular dynamics and experimental validation

“…The radius of gyration (Rg) reflects the binding tightness and degree of constraint of the system, and a higher Rg value is associated with a better chance of producing flexible ligands. Therefore, the higher the Rg value, the lower the stability [ 30 ]. As shown in Figure 7C , HSP90AA1-AloinA complex has the smallest Rg below 1.75 and is the most stable system compared to the Rg of AKT1-AloinA complex and AKT3-AloinA complex, their values fluctuate around 2.1 nm and 3.1 nm, respectively.…”

3D-QSAR and Molecular Dynamics Study of Isoxazole Derivatives to Identify the Structural Requirements for Farnesoid X Receptor (FXR) Agonists

“…A GAFF force field was employed to parameterize the ligands, while an AMBER ff99SB force field was employed for the FXR structure [ 63 , 64 ]. Each complex was then solvated by a dodecahedral box of TIP3P waters, with the solute 10 Å away from the water box boundary [ 65 ]. After the solvation, each system was neutralized with the addition of counter ions (NaCl), followed by the steepest descent method, in order to minimize the system every 1000 steps.…”