The collection of crystal (X-ray) structure of the enzymes [PDB: 3RX3 (aldose reductase), 3W37 (α-glucosidase), and 1DHK (α-amylase)] were from the RSCB Protein Data Bank (
https://www.rcsb.org/ accessed on 12 December 2020). The UCSF Chimera software V1.14 was used in the preparation of the enzymes in readiness for docking [40 (
link)], PubChem (
https://pubchem.ncbi.nlm.nih.gov/ accessed on 15 December 2020) was used to retrieve the structures of the chromatogram-identified phenolic compounds (sinapic acid, cacticin, hyperoside, 1,3-dicaffeoxyl quinic acid, procyanidin, rutin, epicatechin, isorhamnetin-3-
O-rutinoside, chlorogenic acid, myricetin and luteolin-7-
O-beta-
d-glucoside) and standards (acarbose and ranirestat) and optimization of their three-dimensional structures executed using Avogadro software as previously reported [41 (
link)]. The optimized compounds (ligands) and the enzymes were subsequently subjected to molecular docking.
The docking of the prepared phenolic compounds and standards into binding pockets of the enzymes (α-amylase, α-glucosidase, and aldose reductase) was by Autodock Vina Plugin on Chimera V1.14. Judging by the docking scores, complexes identified to have the best pose for each compound were ranked, selected and further analyzed through 100 ns molecular dynamics simulation (MDS).
The MDS was achieved as recently reported [28 (
link)], using the GPU (force fields) version obtainable in AMBER package, where the description of the system by FF18SB variant of the AMBER force field was carried out [42 (
link)]. With the aid of Restrained Electrostatic Potential (RESP) and the General Amber Force Field (GAFF) methods of the ANTECHAMBER assisted with information on atomic partial charges for the compounds. Hydrogen atoms and Na+ and Cl- counter ions (to neutralize the system) were made possible with Leap module of AMBER 18. The residues were numbered 1–336, 913, and 496, respectively, for aldose reductase, α-glucosidase and α-amylase. The system in each case was then lowered implicitly within an orthorhombic box of TIP3P water molecules such that all atoms were within 8Å of any box edge. MDS total time carried-out were 100 ns. For each simulation, hydrogens atoms were constricted using the SHAKE algorithm. The step size of each simulation was 2 fs, and an SPFP precision model was used. The simulations align with the isobaric-isothermal ensemble (NPT), having randomized seeding, Berendsen barostat maintains 1 bar constant pressure, 2 ps pressure-coupling constant, 300 K temperature and Langevin thermostat with a collision frequency of 1.0 ps [43 (
link)].
Using PTRAJ, the systems were subsequently saved, and each trajectory analyzed every 1 ps, and the RoG, RMSF, and RMSD were analyzed with CPPTRAJ module (AMBER 18 suit).
Molecular Mechanics/GB Surface Area method (MM/GBSA) was adopted to assess the free binding energy while comparison of the systems binding affinity followed afterwards [44 (
link)]. Binding free energy was averaged over 100,000 snapshots extracted from the 100 ns trajectory. The ΔG for each system (enzyme, complex and phenolics) was estimated as earlier reported [45 (
link)].
Sabiu S., Balogun F.O, & Amoo S.O. (2021). Phenolics Profiling of Carpobrotus edulis (L.) N.E.Br. and Insights into Molecular Dynamics of Their Significance in Type 2 Diabetes Therapy and Its Retinopathy Complication. Molecules, 26(16), 4867.