Molecular docking studies and molecular dynamic simulation analysis: To identify novel ATP-competitive inhibition of Glycogen synthase kinase-3β for Alzheimer’s disease

Software used for this computational study

Maestro was used for our computational simulation studies; the graphical interface was Schrödinger. In the present work, utilised different modules that were used in previous studies.: Protein Preparation Wizard, Ligand Preparation with the help of LigPrep, for ADME prediction, Qikprop was used, for Docking studies the Glide module was used, and binding energy prediction Prime was used and MDs studies was done by Desmond.^40-42 Alternatively freely available software are AutoDock 4, ArgusLab, and Gromacs.

Selection and preparation of protein

In the present study PDB I.D: 6Y9R (http://www.rcsb.org/)⁴³ of GSK-3β was selected, which is bound to a co-crystallized ligand with a Resolution of 2.08 Å consisting of only chain A, downloaded from the protein databank site,^41,42 and imported this PDB I. D to the Maestro interface. This is followed by preparing the protein molecule with the help of the protein preparation wizard panel of the Schrödinger suite. Alternatively, the AutoDock 4 which is freely available software can be used. The first step for protein preparation is preprocessing in the workspace the protein structure is selected then Assigned bond orders, using the CCD database, explicit hydrogen to the structure, zero-order bonds to metals, creating disulfide bonds, and optimizing missing side chain atoms by running a Prime job, fill in missing loops by running the Prime job was selected and delete waters that are further than specified distance beyond 5.00 Å from any of the het groups including ions. This is recommended for Glide and virtual screening, but MD applications should keep these water as they will help equilibrate the solvent box faster and generate het states using Epik pH: 7.5±0.0, and clicked for preprocessing. The second step is to review and modify here when a new chain, water, or ht is selected to zoom to fit the selection to the workspace, select waters and hets within 5.0 Å of selected chains, and keep the remaining residues and chains as the default setting and generate a state pH of 7.5±0.0. The third step is to refine the selected sample water orientations in addition to other groups, and the protonation states of residues and ligands were set by using PROPKA pH: 7.5 then it simulated the exact experimental conditions and clicked to automatically optimize hydroxyl, ASN, GLN, HIS states using ProtAssign, removing water that is further than beyond hets 3.0 Å, including ions. This is recommended for Glide and virtual screening, but MD applications should keep this water as it will help equilibrate the solvent box faster. OPLS3e force field for the restrained minimization of proteins was done. Then after the minimization, the glycerol and acetate ions were removed from the minimized protein.

Ligand selection and preparation

IMID 2 scaffolds were used for this study. This scaffold possesses fewer hydrogen bond donors and improved CNS (Central Nervous System) penetration at micromolar concentrations (μM), based on studies by Buonfiglio et al.⁴⁴ To generate the focus library, draw an IMID 2 scaffold structure at Enamine and subjected it to a substructure search. Then, enamine realdb_IMID2_Scaffold_molecules (1400 molecules) was downloaded, followed by ligand preparation using LigPrep in the Schrödinger suite software (Schrödinger 2021-3).⁴⁵ Alternatively, the AutoDock 4 which is freely available software can be used. The 3D coordinates for the compounds were generated using the LigPrep tool, and the Epik module predicted the most probable charge form for the ionization state of compounds at pH: 7.5±0.0, generated tautomers and stereoisomers, determined chiralities from the 3D structure, generated at most 32 per ligand, and finally, the minimization of the drug molecule process was performed using the force field OPLS3e.^46,47

ADME prediction

After ligand preparation, all ligands were subjected to QikProp in the Schrödinger suite software to predict the ADME profile of the compounds. QikProp predicts the drug-like properties of all selected compounds, such as molecular weight, Hydrogen Bond Donors, Hydrogen Bond Acceptor, QPlogS (predicted aqueous solubility), QPlogPo/w (predicted octanol/water partition coefficient), QPlogBB (predicted brain/blood partition coefficient), QPlogPw (predicted water/gas partition coefficient), QPlogPC16 (predicted hexadecane/gas partition coefficient), QPlogPoct (predicted octanol/gas partition coefficient), PSA (Van der Waals surface area of polar nitrogen and oxygen atoms), QPPCaco (predicted apparent Caco-2 cell permeability), and QPPMDCK (predicted apparent MDCK cell permeability).⁴⁸

Molecular docking studies

Here, Schrödinger suite Glide module was used for better prediction of the protein-ligand involved in assessing the fitting of all conformations of compounds at the binding site followed by ranking and modes Schrödinger Release 2021-3. The alternative freely accessible software is AutoDock and Gromacs. The molecular docking study process was used for assessment before molecular docking studies, and the default settings were used for the Receptor Grid generation module. The binding site was recognized using the Receptor Grid generation module, and this binding site was specified as a box of 10 × 10 × 10 ų centered on the centroid of the co-crystallized ligand of GSK-3β, which was taken into consideration for grid generation.⁴⁹ Then subject all prepared compounds to the SP (standard precision) scoring function mode.⁵⁰ Glide SP docking ligand sampling is in flexible mode and selected sample nitrogen inversions, sample ring conformations, and amides only penalize non-planar conformations, adding Epik state penalties to the docking score. In the output file, the pose viewer file that includes the receptor, written out at most one pose per ligand, performed post-docking minimization, the number of poses per ligand to include was five, and the RMSD was computed to input ligand geometries and run the job. Finally, the results were analyzed, separated, and selected based on the VAL135 residue interaction and physicochemical properties reported in the literature. These compounds were then subjected to the prime-MMGBSA module.

MM-GBSA binding energy

MM-GBSA (molecular mechanics-generalized born energies & surface area/accessibility continuum solvation method) was used to calculate the binding free energy and binding analysis of 523 SP docked protein-ligand complexes. MM-GBSA exhibits many energy properties for proteins, ligands, and complex structures. Prime MM-GBSA analysis is based on the solvent model of VSGB 2.0, and the OPLS3e force field is used to calculate the binding affinity of the respective protein-ligand complex.⁵¹

MDs with Desmond

This model was used to precisely predict the interaction of the ligand with 6Y9R (protein), the stability of their binding under physiological conditions, and to analyze their motion at the molecular and atomic levels. The protein-ligand complex was subjected to a system built here, and different solvent models were selected; an orthorhombic box with dimensions of 10 Å × 10 Å × 10 Å was used to determine the structure geometry to have the minimum box volume, and the checked boundary box is shown in the workspace. The next step is to load the previous system builder file and enter 100ns in simulation time, 1000frames were captured, ensemble class is NPT (Normal Pressure Temperature), 300°K temperature, pressure bar is 1, and relax model system before simulation was used for MDs studies.⁵² After MDs, the analysis was done with the help of the RMSD (Root Mean Square Deviation), the RMSD of ligands, RMSD of protein. RMSD is mainly used to determine the stability of the protein-ligand complex, RMSF (Root Mean Square Fluctuations), ligand-protein contacts, and ligand interaction diagram. Alternatively, NAMD molecular simulation along with VMD software can be used which is freely available.

ADME prediction

The ADME or drug-likeness parameters are the most important and helpful parameters for screening compounds after SP docking studies during the drug discovery process.⁵³ Generally, for CNS drugs, the physicochemical properties range from molecular weight 100 Da to 450 Da, hydrogen bond acceptors range from 0 to 5, hydrogen bond donor values range from 0 to 3, and the topological polar surface area value is ≤76 Å2.⁵⁴ The results of this study are listed in Table 2. The present study results are having lesser molecular weight and having similar hydrogen bond donor compared to previous ADME properties of 13 trial compounds with an IMID2 core scaffold molecular weight of 393.48, the hydrogen bond donor was 1. In the14 trial compound of the IMID2 core scaffold, the molecular weight was 451.56, and the hydrogen bond donors were 1. In the 15 trial compounds of the IMID2 core scaffold, the molecular weight was 336.39, and the hydrogen bond donors were 1. In 16 trials of the IMID2 core scaffold, the molecular weight was 394.47, and the number of hydrogen bond donors was 1.⁴⁴

The Z3336252116 with QPlogS value -1.653, QPlogPo/w value was 0.224, QPlogPw value was 14.263, QPlogBB value was -1.042, QPlogPC16 value was 8.612, QPlogPoct values was 15.179, QPPCaco values was 107.434, QPPMDCK values was 226.003. Z3136198649 with QPlogS value -1.249, QPlogPo/w value was 0.180, QPlogPw value was 12.291, QPlogBB value was -0.286, QPlogPC16 value was 8.140, QPlogPoct values was 14.830, QPPCaco values was 155.541, QPPMDCK values was 73.234. PV-005996498401 with QPlogS value -1.416, QPlogPo/w value was 0.223, QPlogPw value was 14.523, QPlogBB values with -1.051, QPlogPC16 value was 8.453, QPlogPoct values was 14.948, QPPCaco values was 159.830, QPPMDCK values was 143.876. PV-0025895352867 with QPlogS value -2.886, QPlogPo/w value was 1.532, QPlogPw value was 11.558, QPlogBB value was -0.639, QPlogPC16 value was 9.016, QPlogPoct values was 14.584, QPPCaco values was 875.662, QPPMDCK values was 428.57. Z3603284828 with QPlogS value -3.257, QPlogPo/w value was 2.117, QPlogPw value was 9.797, QPlogBB value was -0.530, QPlogPC16 value was 8.653, QPlogPoct values was 13.900, QPPCaco values was 1320.626, QPPMDCK values was 668.188. Z3136169169 with QPlogS value -1.258, QPlogPo/w value was 0.181, QPlogPw value was 12.293, QPlogBB value was -0.290, QPlogPC16 value was 8.140, QPlogPoct values was 14.831, QPPCaco values was 155.076, QPPMDCK values was 72.997. Z2497631540 with QPlogS value -3.149, QPlogPo/w value was 1.655, QPlogPw value was 10.508, QPlogBB value was -0.507, QPlogPC16 value was 8.503, QPlogPoct values was 14.075, QPPCaco values was 956.136, QPPMDCK values was 1160.160. Z4468779454 with QPlogS value -2.896, QPlogPo/w value was 1.907, QPlogPw value was 10.321, QPlogBB value was -0.664, QPlogPC16 value was 8.871, QPlogPoct values was 14.100, QPPCaco values was 1311.747, QPPMDCK values was 663.334. Z2708235883 with QPlogS value -2.771, QPlogPo/w value was 1.393, QPlogPw value was 10.536, QPlogBB value was -0.559, QPlogPC16 value was 7.537, QPlogPoct values was 13.626, QPPCaco values was 953.240, QPPMDCK values was 810.911.

Among all there parameter the QPlogBB is a useful parameter to predict or reflect directly the ability of a compound to cross the blood-brain barrier partition coefficient.⁵⁵ The recommended range of QPlogBB should be within -3.0 to 1.2.⁵⁶ If the compound is too polar then the compounds don’t cross the BBB.⁵⁷ The QPlogBB negative values indicate that the compounds are poorly permeable and polar.^58,59

MDs analysis

Compounds were selected based on the screening process described above. Only short-listed compounds were subjected to desmonds for MDs. These studies provide information related to the stability and flexibility of the molecular docking complexes. Performed this study using the MM-GBSA. This study provides information related to protein-ligand interactions, RMSD variation can be assessed, and RMSF fluctuations can be assessed for protein-ligand complexes. The complex was regarded as stable if it fell within the 3 Å range. RMSD generally represents the average difference between the displacements of the atoms at an instant of target structures to the reference structure. As per MD, in comparison to the starting point, how the part of the structure changes or the entire structure over time or to identify large changes in protein structure.

RMSF is a calculation of a particular residue or individual residue fluctuation/flexibility or group of atoms during a simulation, related to the reference structure or average structure of simulation.⁶⁰ These are the compounds.

IMID2 PV-002589535286: During MD simulations, the ligand-protein complex showed hydrogen-bonded interactions. The protein backbone and ligand structure exhibited higher RMSD fluctuations over the first 0-20 ns. The protein backbone and ligand fluctuations stayed within the range of 0.6 Å and 1.3 Å over the last 80 ns. The amino acid residue VAL135 formed a 95% direct hydrogen bonding interaction with the amide carbonyl, the amino acid residue VAL135 formed a 42% direct hydrogen bonding interaction with imidazole in the ring structure, the amino acid residue TYR134 had 33% direct hydrogen bonding interaction, and ARG141 had a 72% PI-cation interaction with 1H-imidazole. The amide carbonyl was exposed to an H₂O molecule through which it interacted with residue ILE62 (49%), as shown in Figure 2; however, it was reported in 13 triads (previous literature) that the amide carbonyl was exposed to water-mediated contacts with residues such as ILE62 and GLN185 for the IMID2 core scaffold. The 16-triad isopropoxy group confirmed the probable water molecule interaction with the residue ASP200. The 17 triads showed an indirect interaction between the ortho substituent and residue ASN186. In the 18 triad, the di-F phenyl group was orientated towards LYS85, thus promoting an electrostatic interaction. However, VAL135 interaction formed direct hydrogen bonding interactions with amide carbonyl in all 13, 16, 17, and 18 trial compounds of previous IMID2 scaffolds.⁴⁴

Figure 2. Pictorial representation of 6Y9R interactions with the IMID2 compound code PV-002589535286 monitored during MDs trajectory (A) Protein-Ligands RMSD; (B) Protein-Ligand Contacts; (C) Ligand-Protein Contacts; (D) Protein-Ligand Contacts described as histogram; (E) Protein RMSF.

Interaction diagram of 6Y9R with IMID2 compound code PV-002589535286 was observed during RMSD, Ligand interaction with amino acid residues of Protein Contact during MDs, protein-ligand contacts (cont.) (On the x-axis trajectory frame, the number is present, and amino acid residues are seen on the y-axis. The amino acid residues were in greater contact with ligands in the trajectory frame that appeared as a dark color shade). The VAL135 interaction fraction approached 1.4, whereas the ASP133 interaction fraction was missing in comparison to the core IMID2 scaffold from the previous literature. Protein root mean square fluctuation (RMSF) is mostly helpful for predicting changes that occur locally along the enzyme chain. These peaks provide information on how much the protein fluctuated during the MDs study.

IMID2 compound code Z3603284828: During the MDs, the ligand-protein complex exhibited hydrogen-bonded interactions. The protein backbone and ligand structure exhibited higher RMSD fluctuations during the first 20 ns. The protein backbone and ligand fluctuations stayed within the range of 0.3 Å and 0.6 Å over the last 70 ns shown in Figure 3. The amino acid residue VAL135 formed 93% direct hydrogen bonding interaction with the amide carbonyl, the amino acid residue VAL135 formed 55% direct hydrogen bonding interaction with imidazole, and the amino acid residue ILE62 formed 35% water-mediated direct hydrogen bonding interaction with 1-hydroxycyclopropyl.

Figure 3. Pictorial representation of 6Y9R interaction with the IMID2 compound code Z3603284828 monitored during MDs trajectory (A) Protein-Ligands RMSD; (B) Protein-Ligand Contacts; (C) Ligand-Protein Contacts; (D) Protein-Ligand Contacts described as histogram; (E) Protein RMSF.

Interaction diagram of 6Y9R with IMID2 compound code Z3603284828 observed during RMSD, Ligand interaction with amino acid residues of Protein Contact during MDs, protein-ligand contacts (cont.) (On the x-axis, a trajectory frame with a number is present, and amino acid residues are seen on the y-axis. The amino acid residues that were in greater contact with ligands in the trajectory frame appeared as a dark color shade). The VAL135 interaction fraction approached 1.4, whereas the ASP133 interaction fraction was missing in comparison to the core IMID2 scaffold from previous literature, and protein RMSF is mostly helpful for predicting the changes that occur locally along the enzyme chain. These peaks provide information on how much the protein fluctuated during the MDs study.

IMID2 compound code Z3136169169: During the MDs, the ligand-protein complex exhibited hydrogen-bonded interactions. The protein backbone and ligand structure exhibited higher RMSD fluctuations during the first 20 ns. The protein backbone and ligand fluctuations stayed within the range of 0.8 Å and 0.3 Å over the last 70 ns shown in Figure 4. The amide carbonyl is exposed to hydrogen bonding and interacts with residue VAL135, and the amino acid residue VAL135 forms 76% direct hydrogen bonding interactions with imidazole.

Figure 4. Pictorial representation of 6Y9R interaction with the IMID2 compound code Z3136169169 monitored during MDs trajectory (A) Protein-Ligand RMSD; (B) Protein-Ligands Contacts; (C) Ligand-Protein Contacts; (D) Protein-Ligand Contacts described as a histogram; (E) Protein RMSF.

Interaction diagram of 6Y9R with IMID2 compound code Z3136169169 observed during RMSD, Ligand interaction with amino acid residues of Protein Contact during MDs, protein-ligand contacts (cont.) (On the x-axis trajectory frame, the number is present, and amino acid residues are seen on the y-axis. The amino acid residues that were in greater contact with ligands in the trajectory frame appeared as a dark color shade). The VAL135 interaction fraction approached 1.6, whereas the ASP133 interaction fraction was missing in comparison to the core IMID2 scaffold from previous literature and protein RMSF is mostly helpful for predicting the changes that occur locally along with the enzyme chain. These peaks provide information on how much the protein fluctuated during the MDs study.

IMID2 compound code Z2497631540: During the MDs, the ligand-protein complex exhibited hydrogen-bonded interactions. The protein Cα and ligand structures exhibited higher RMSD fluctuations over the first 20 ns. The protein backbone and ligand fluctuations stayed within the range of 0.4 Å and 3.7 Å over the last 70 ns shown in Figure 5. The amino acid residue VAL135 formed 84% direct hydrogen bonding interactions with the amide carbonyl, and the amino acid residue VAL135 formed 39% direct hydrogen bonding interactions with imidazole in the ring structure.

Figure 5. Pictorial representation of 6Y9R interaction with the IMID2 compound code Z2497631540 monitored during MDs trajectory (A) Protein-Ligand RMSD; (B) Protein-Ligand Contacts; (C) Ligand-Protein contacts; (D) Protein-Ligand Contacts described as a histogram; (E) Protein RMSF.

Interaction diagram of 6Y9R with IMID2 compound code Z2497631540 observed during RMSD, Ligand interaction with amino acid residues of Protein Contact during MDs, protein-ligand contacts (cont.) (On the x-axis, a trajectory frame with a number is present, and amino acid residues are seen on the y-axis. The amino acid residues that were in greater contact with ligands in the trajectory frame appeared as a dark color shade). The VAL135 interaction fraction approached 1.2, whereas the ASP133 interaction fraction was missing in comparison to the core IMID2 scaffold from previous literature, and protein RMSF is mostly helpful for predicting the changes that occur locally along the enzyme chain. These peaks provide information on how much the protein fluctuated during the MDs study.

IMID2 compound code Z4468779454: During the MDs, the ligand-protein complex exhibited hydrogen-bonded interactions. The protein Cα and ligand structures exhibited higher RMSD fluctuations over the first 20 ns. The protein backbone and ligand fluctuations stayed within the range of 0.5 Å and 0.4 Å over the last 70 ns shown in Figure 6. The amino acid residue VAL135 formed a 96% direct hydrogen bonding interaction with the amide carbonyl, the amino acid residue VAL135 formed a 77% direct hydrogen bonding interaction with imidazole, and the amide carbonyl was exposed to the H₂O molecule through which it interacted with residue ILE62, with 43% and 61% interaction with 1-hydroxypentan-2-yl.

Figure 6. Pictorial representation of 6Y9R interaction with the IMID2 compound code Z4468779454 monitored during MDs trajectory (A) Protein-Ligand RMSD; (B) Protein-Ligand contacts; (C) Ligand-Protein contacts; (D) Protein-Ligand contacts described as a histogram; (E) Protein RMSF.

Interaction diagram of 6Y9R with IMID2 compound code Z4468779454 observed during RMSD, Ligand interaction with amino acid residues of Protein Contact during MDs, protein-ligand contacts (cont.) (On the x-axis trajectory frame, the number is present, and amino acid residues are seen on the y-axis. The amino acid residues were in greater contact with ligands in the trajectory frame that appeared as a dark color shade). The VAL135 interaction fraction approached 1.75, whereas the ASP133 interaction fraction was missing in comparison to the core IMID2 scaffold from previous literature, and protein RMSF is mostly helpful for predicting the changes that occur locally along the enzyme chain. These peaks provide information on how much the protein fluctuated during the MDs study.

IMID2 compound code Z2708235883: During the MDs, the ligand-protein complex exhibited hydrogen-bonded interactions. The protein backbone and ligand structure exhibited higher RMSD fluctuations during the first 20 ns. The protein Cα and ligand fluctuations stayed within the range of 0.5 Å and 0.2 Å over the last 70 ns. The amino acid residue VAL135 formed 93% direct hydrogen bonding interactions with the amide carbonyl, and the amino acid residue VAL135 formed 53% direct hydrogen bonding interactions with imidazole in the ring structure, as shown in Figure 7.

Figure 7. Pictorial representation of 6Y9R interaction with the IMID2 compound code Z2708235883 monitored during MDs trajectory (A) Protein-Ligand RMSD; (B) Protein-Ligand Contacts; (C) Ligand-Protein Contacts; (D) Protein-Ligand Contacts described as a histogram; (E) Protein RMSF.

Interaction diagram of 6Y9R with IMID2 compound code Z2708235883 observed during RMSD, Ligand interaction with amino acid residues of Protein Contact during MDs, protein-ligand contacts (cont.) (On the x-axis trajectory frame, the number is present, and amino acid residues are seen on the y-axis. The amino acid residues that have been in more contact with ligands in the trajectory frame appear to have a dark color shade). The VAL135 interaction fraction approached 1.4, whereas the ASP133 interaction fraction was missing in comparison to the core IMID2 scaffold from previous literature, and protein RMSF is mostly helpful for predicting the changes that occur locally along the enzyme chain. These peaks provide information on how much the protein fluctuated during the MDs study.

IMID2_compound code Z3336252116: During MD simulations, the ligand-protein complex showed hydrogen-bonded interactions. The protein backbone and ligand structure exhibited higher RMSD fluctuations over the first 0-20 ns. The protein backbone and ligand fluctuations stayed within the range of 0.6 Å and 2.2 Å over the last 80 ns. The amino acid residue VAL135 formed 93% direct hydrogen bonding interactions with the amide carbonyl, and the amino acid residue VAL135 formed 48% direct hydrogen bonding interactions with imidazole in the ring structure, which was exposed to the H₂O molecule through which it interacted with residue ILE62 (43%), as shown in Figure 8.

Figure 8. Pictorial representation of 6Y9R interaction with the IMID2 compound code Z3336252116 monitored during MDs trajectory (A) Protein-Ligand RMSD; (B) Protein-Ligand Contacts; (C) Protein-Ligand Contacts described as a histogram; (D) Ligand-Protein Contacts; (E) Protein RMSF.

Interaction diagram of 6Y9R with IMID2_compound code Z3336252116 observed during RMSD, Ligand interaction with amino acid residues of Protein Contact during MDs, protein-ligand contacts (cont.) (On the x-axis trajectory frame, the number is present, and amino acid residues are seen on the y-axis. The amino acid residues that have been in more contact with ligands in the trajectory frame appear to have a dark color shade). The VAL135 interaction fraction approached 1.5, whereas the ASP133 interaction fraction was missing in comparison to the core IMID2 scaffold from the previous literature. Protein RMSF is mostly helpful for predicting changes that occur locally along the enzyme chain. These peaks provide information on how much the protein fluctuated during MDs study.

IMID2 compound code Z3136198649: During MDs, the ligand-protein complex showed hydrogen-bonded interactions. The protein backbone and ligand structure exhibited higher RMSD fluctuations over the first 0-20 ns. The protein backbone and ligand fluctuations stayed within the range of 0.7 Å and 0.5 Å over the last 80 ns. The amino acid residue VAL135 formed 83% direct hydrogen bonding interaction with amide carbonyl, and the amino acid residue VAL135 formed 49% direct hydrogen bonding interaction with imidazole in the ring structure given in Figure 9.

Figure 9. Pictorial representation of 6Y9R interactions with the IMID2 compound code Z3136198649 monitored during MDs trajectory (A) Protein-Ligands RMSD; (B) Protein-Ligand Contacts; (C) Protein-Ligand Contacts described as a histogram; (D) Ligand-Protein Contacts; (E) Protein RMSF.

Interaction diagram of 6Y9R with IMID2 compound code Z3136198649 observed during RMSD, Ligand interaction with amino acid residues of Protein Contact during MDs, protein-ligand contacts (cont.) (On the x-axis trajectory frame, the number is present, and amino acid residues are seen on the y-axis. The amino acid residues that have been in more contact with ligands in the trajectory frame appear to have a dark color shade). The VAL135 interaction fraction approached 1.4, whereas the ASP133 interaction fraction was missing in comparison to the core IMID2 scaffold from the previous literature. Protein RMSF is mostly helpful for predicting changes that occur locally along the enzyme chain. These peaks provide information on how much the protein fluctuated during the MDs study.

IMID2 compound code PV-005996498401: During MDs, the ligand-protein complex showed hydrogen-bonded interactions. The protein backbone and ligand structure exhibited higher RMSD fluctuations over the first 0-20 ns. The protein backbone and ligand fluctuations stayed within the range of 1.0 Å and 0.6 Å over the last 80 ns. The amino acid residue VAL135 formed 90% direct hydrogen bonding interactions with the amide carbonyl, and the amino acid residue VAL135 formed 47% direct hydrogen bonding interactions with imidazole in the ring structure, as shown in Figure 10.

Figure 10. Pictorial representation of 6Y9R interaction with the IMID2 compound code PV-005996498401 monitored during MDs trajectory (A) Protein-Ligand RMSD; (B) Protein-Ligand contacts; (C) Protein-Ligand contacts described as a histogram; (D) Ligand-Protein contacts; (E) Protein RMSF.

Interaction diagram of 6Y9R withIMID2 compound code PV-005996498401 observed during RMSD, Ligand interaction with amino acid residues of Protein Contact during MDs, protein-ligand contacts (cont.) (On the x-axis trajectory frame, the number is present, and amino acid residues are seen on the y-axis. The amino acid residues that have been in more contact with ligands in the trajectory frame appear to have a dark color shade). The VAL135 interaction fraction approached 1.85, whereas the ASP133 interaction fraction was missing in comparison to the core IMID2 scaffold from the previous literature. Protein RMSF is mostly helpful for predicting changes that occur locally along the enzyme chain. These peaks provide information on how much the protein fluctuated during the MDs study.

Ethical considerations

Ethics and written consent were not applicable.