In silico investigation of potential phytoconstituents against ligand- and voltage-gated ion channels as antiepileptic agents.

The most promising anticonvulsant phytocompounds were explored in this work using docking, molecular dynamic (MD) simulation, and Molecular Mechanics-Poisson-Boltzmann Surface Area (MM-PBSA) approaches. A total of 70 phytochemicals were screened against α-amino-3-hydroxyl-5-methyl-4-isoxazole propionic acid (AMPA), N-methyl-d-aspartate (NMDA), voltage-gated sodium ion channels (VGSC), and carbonic anhydrase enzyme II (CA II) receptors, and the docking results were compared to the reference drug phenytoin. Amentoflavone displayed the highest affinity for AMPA and VGSC receptors, with docking scores of - 10.4 and - 10.1 kcal/mol, respectively. Oliganthin H-NMDA and epigallocatechin-3-gallate-CA II complexes showed docking scores of - 10.9 and - 6.9 kcal/mol, respectively. All four complexes depicted a high dock score compared to the phenytoin complex at the binding site of the corresponding proteins. The MD simulation investigated the stabilities and favorable conformation of apoproteins and ligand/reference-bound complexes. The results revealed that proteins AMPA, VGSC, and CA II were more efficiently stabilized by lead phytochemicals than phenytoin binding. Additionally, principal component analysis and MM-PBSA results suggested that these lead phytocompounds have good compactness and strong binding free energy. Further, physicochemical and pharmacokinetic studies revealed that these final lead phytochemicals would be suitable for oral intake, have sufficient intestinal permeability, and have the ability to cross the blood-brain barrier (BBB). Comprehensively, this study predicted amentoflavone as the best lead phytochemical out of the 70 anticonvulsant phytocompounds that can be used to treat epilepsy. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-024-03948-1.

QM/MM reveals the sequence of substrate binding during OPRT action.

Computational investigation of orotate phosphoribosyltransferase (OPRT) action, an enzymatic reaction between phosphoribosyl pyrophosphate (PRPP) and orotic acid (OA) to yield orotidine 5'-monophosphate (OMP), was carried out. Insights into the pathways of the substrate attack step of the reaction were developed under the quantum mechanics/molecular mechanics framework with S. cerevisiae strain as the representative enzyme bearer. Four pathways were proposed for PRPP and OA binding differing in the sequence of PRPP, OA and Mg2+ ion complexation with OPRT. The formation of Mg2+-OPRT complex was accompanied by a small energy change while the largest stabilization was observed for the formation of Mg2+-PRPP complex supporting the experimental observation of Mg2+-PRPP complex as the true substrate for the reaction. Formation of PRPP-OPRT complex was found to be energetically not probable rendering the pathway requiring Mg2+-OA complex not probable. Further, PRPP migration towards the active site was found to be energetically not favoured rendering the pathway involving Mg2+-OA complexation improbable. Migration of OA and Mg2+-PRPP complex towards the active site was found to be energetically probable with a large stabilization of the system when Mg2+-PRPP complex bound to the OA-OPRT complex. This conclusively proved the sequential binding of OA and Mg2+-PRPP complexes during OPRT action.

QM/MM analysis of effect of divalent metal ions on OPRT action.

The role of Mg2+ cofactor in orotate phosphoribosyltransferase (OPRT) catalyzed synthesis of orotidine monophosphate (OMP) from phosphoribosyl pyrophosphate (PRPP) and orotate (OA) in substrate binding and the influence of the identity of the divalent metal ion on the reaction mechanism were addressed in this study using quantum mechanics/molecular mechanics framework. Energetics of migration and binding of different substrate complexes in the active site cavity was established. A quantitative analysis of various processes indicated the reaction pathway to consist of complexation of Mg2+ with PRPP, migration of Mg2+-PRPP and OA towards the active site, binding of OA to OPRT, and binding of Mg2+-PRPP complex to OA-OPRT complex. The mechanism of the reaction was unaltered by the change in the identity of divalent metal ion. Experimentally reported inhibiting character of Co2+ was explained on the basis of large Co2+-PRPP binding and migration energies. Mg2+, Ca2+, Mn2+, Co2+ and Zn2+ ions were screened computationally to assess their inhibiting/activating characteristics. Trends obtained by our computational investigations were in correspondence with experimentally reported trends.