Structural basis of pesticide detection by enzymatic biosensing: a molecular docking and MD simulation study

Designing of rapid, facile, selective, and cost-effective biosensor technology is a growing area for the detection of various classes of pesticides. The biosensor with these features can be achieved only through the various bio-components using different transducers. This study, therefore, focuses on the usage of molecular docking, specificity tendencies, and capabilities of proteins for the detection of pesticides. Accordingly, the four transducers, acetylcholinesterase (ACH), cytochromes P450 (CYP), glutathione S-transferase (GST), and protein kinase C (PKC) were selected based on their applications including neurotransmitter, metabolism, detoxification enzyme, and protein phosphorylation. Then after molecular docking of the pesticides, fenobucarb, dichlorodiphenyltrichloroethane (DDT), and parathion onto each enzyme, the conformational behavior of the most stable complexes was further analyzed using 50 ns Molecular Dynamics (MD) simulations carried out under explicit water conditions. In the case of protein kinase C (PKC) and cytochrome P450 3A4 enzyme (CYP), the fenobucarb complex showed the most suitable combination of free energy of binding and inhibition constant ?4.42 kcal/mol (573.73 μM) and ?5.1 kcal/mol (183.49 μM), respectively. Parathion dominated for acetylcholinesterase (ACH) with ?4.57 kcal/mol (448.09 μM) and lastly dichlorodiphenyltrichloroethane for glutathione S-transferase (GST), ?5.43 kcal/mol (103.88 μM). The RMSD variations were critical for understanding the impact of pesticides as they distinctively influence the energetic attributes of the proteins. Overall, the outcomes from the extensive analysis provide an insight into the structural features of the proteins studied, thereby highlighting their potential use as a substrate in biorecognition sensing of pesticide compounds.     

Dynamics of fluoroquinolones induced resistance in DNA gyrase of Mycobacterium tuberculosis

DNA gyrase is a validated target of fluoroquinolones which are key components of multidrug resistance tuberculosis (TB) treatment. Most frequent occurring mutations associated with high level of resistance to fluoroquinolone in clinical isolates of TB patients are A90V, D94G, and A90V–D94G (double mutant [DM]), present in the larger subunit of DNA Gyrase. In order to explicate the molecular mechanism of drug resistance corresponding to these mutations, molecular dynamics (MD) and mechanics approach was applied. Structure-based molecular docking of complex comprised of DNA bound with Gyrase A (large subunit) and Gyrase C (small subunit) with moxifloxacin (MFX) revealed high binding affinity to wild type with considerably high Glide XP docking score of ?7.88 kcal/mol. MFX affinity decreases toward single mutants and was minimum toward the DM with a docking score of ?3.82 kcal/mol. Docking studies were also performed against 8-Methyl-moxifloxacin which exhibited higher binding affinity against wild and mutants DNA gyrase when compared to MFX. Molecular Mechanics/Generalized Born Surface Area method predicted the binding free energy of the wild, A90V, D94G, and DM complexes to be ?55.81, ?25.87, ?20.45, and ?12.29 kcal/mol, respectively. These complexes were further subjected to 30 ns long MD simulations to examine significant interactions and conformational flexibilities in terms of root mean square deviation, root mean square fluctuation, and strength of hydrogen bond formed. This comparative drug interaction analysis provides systematic insights into the mechanism behind drug resistance and also paves way toward identifying potent lead compounds that could combat drug resistance of DNA gyrase due to mutations.     

In silico identification of mimicking molecule(s) triggering von Willebrand factor in human: a molecular drug target for regulating coagulation pathway

Abstract

Blood coagulation is a complex and dynamic process wherein the body activates its emergency mechanism to stop bleeding and wound healing via the interactions of prothrombotic and antithrombotic agents. von Willebrand factor (VWF) is a complex glycoprotein and initial component of the hemostasis pathway which serves a multipurpose role in blood coagulation process. There are reports of various plants that contain several bioactive compounds possessing properties of inducing blood coagulation directly or indirectly. In the present study, efforts have been made to identify bioactive compounds that may play a significant role in regulation of the coagulation cascade by accelerating VWF and thus enhance the hemostasis process. An antidiuretic peptide drug, Desmopressin, works on VWF and releases them in circulation. Forty-seven compounds from different plant sources were screened through molecular docking, out of which two compounds, Emodin and Peruvianoside II, showed more binding affinity than the reference drug Desmopressin. Emodin and Peruvianoside II showed binding energies ?7.2 and ?7.0?kcal/mol, respectively, when docked with VWF, whereas Desmopressin displayed less binding energy (?6.9?kcal/mol). Emodin belongs to anthraquinone from Rumex hastasus and Peruvianoside II belongs to flavanone glycosides from Thevetia peruviana. The mimicking potential of top identified molecules with respect to the drug was confirmed through simulation analysis. Besides, the molecular dynamics simulation (MDS) study (for 20?ns) showed that the Peruvianoside II protein complex was energetically more stable than Emodin protein complex. Based on the results, Peruvianoside II possesses great potential and thus may be considered for development of drugs for hemostasis.     

Structure-based screening and molecular dynamics simulations offer novel natural compounds as potential inhibitors of Mycobacterium tuberculosis isocitrate lyase

Mycobacterium tuberculosis is the etiological agent of tuberculosis in humans and is responsible for more than two million deaths annually. M. tuberculosis isocitrate lyase (MtbICL) catalyzes the first step in the glyoxylate cycle, plays a pivotal role in the persistence of M. tuberculosis, which acts as a potential target for an anti-tubercular drug. To identify the potential anti-tuberculosis compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,748) against the MtbICL structure. The ligands were docked against MtbICL in three sequential docking modes that resulted in 340 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 27 compounds were found to fit well with re-docking studies. After refinement by molecular docking and drug-likeness analyses, three potential inhibitors (ZINC1306071, ZINC2111081, and ZINC2134917) were identified. These three ligands and the reference compounds were further subjected to molecular dynamics simulation and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the MtbICL and bound complexes. The binding free energy analyses were calculated to validate and capture the intermolecular interactions. The results suggested that the three compounds had a negative binding energy with ?96.462, ?143.549, and ?122.526 kJ mol^?1 for compounds with IDs ZINC1306071, ZINC2111081, and ZINC2134917, respectively. These lead compounds displayed substantial pharmacological and structural properties to be drug candidates. We concluded that ZINC2111081 has a great potential to inhibit MtbICL and would add to the drug discovery process against tuberculosis.     

Comparative docking of dual conformations in human fatty acid synthase thioesterase domain reveals potential binding cavity for virtual screening of ligands

Human fatty acid synthase (hFASN), a homo dimeric lipogenic enzyme with seven catalytic domains, is an important clinical target in cancer, metabolic syndrome and infections. Here, molecular modelling and docking methods were implemented to examine the inter-molecular interactions of thioesterase (TE) domain in hFASN with its physiological substrate, and to identify potential chemical inhibitors. TE catalyses the hydrolysis of thioester bond between palmitate and the 4’ phosphopantetheine of acyl carrier protein, releasing 16-carbon palmitate. The crystal structure of hFASN TE in two inhibitory conformations (A and B) were geometry-optimized and used for molecular docking with palmitate, orlistat (a known FASN inhibitor) and virtual screening against compounds from National Cancer Institute (NCI) database. Relatively, low binding affinity was observed during the complex formation of palmitate with A (?.164 kcal/mol) and B (?.332 kcal/mol) forms of TE, when compared with orlistat-docked TE (A form: ?5.872 kcal/mol and B form: ?5.484 kcal/mol), clearly indicating that the native inhibited conformation (crystal structure) was unfavourable for substrate binding. We used these orlistat dual binding modes as positive controls for prioritizing the ligands during virtual screening. From 2, 31,617 molecules in the NCI database, 916 high-scoring compounds (hit ligands) were obtained for A-form and 4582 for B-form of the TE-domain, which were then ranked according to glide docking score, XP H bond score, absorption, distribution, metabolism and excretion and binding free energy (Prime/MM-GBSA). Consequently, two top scoring ligands (NSC: 319661 and NSC: 153166) emerged as promising drug candidates that may be tested in FASN-over-expressing diseases.     

Designing of benzothiazole derivatives as promising EGFR tyrosine kinase inhibitors: a pharmacoinformatics study

Abstract

Benzothiazole derivatives represent an important class of therapeutic chemical agents and are widely used for interesting biological activities and therapeutic functions including anticancer, antitumor and antimicrobial. In this study, we have performed similarity/substructure-based search of eMolecule database to find out promising benzothiazole derivatives as EGFR tyrosine kinase inhibitors. Several screening criteria that included molecular docking, pharmacokinetics and synthetic accessibility were used on initially derived about 7000 molecules consisting of benzothiazole as major component. Finally, four molecules were found to be promising EGFR tyrosine kinase inhibitors. The best docked pose of each molecule was considered for binding interactions followed by molecular dynamics (MD) and binding energy calculation. Molecular docking clearly showed the final proposed derivatives potential to form a number of binding interactions. MD simulation trajectories undoubtedly indicated that the EGFR protein becomes stable when proposed derivatives bind to the receptor cavity. Strong binding affinity was found for all molecules toward the EGFR which was substantiated by the binding energy calculation using the MM-PBSA approach. Therefore, proposed benzothiazole derivatives may be promising EGFR tyrosine kinase inhibitors for potential application as cancer therapy.

Communicated by Ramaswamy H. Sarma     

Molecular docking studies of benzimidazopyrimidine and coumarin substituted benzimidazopyrimidine derivatives: As potential human Aurora A kinase inhibitors

Protein kinases are important drug targets in human cancers, inflammation and metabolic diseases. Docking studies was performed for all the benzimidazopyrimidine and coumarin substituted benzimidazopyridimine derivatives with human Aurora A kinase target (3FDN) employing flexible ligand docking approach by using AutoDock 4.2. All the compounds were found to have minimum binding energy ranging from -6.26 to -9.29 kJ/mol. Among the molecules tested for docking study, 10-(6-Bromo-2-oxo- 2H-chromen-4-ylmethyl)-2-isopropyl-10H-benzo[4,5]imidazo[1,2-a]pyrimidin-4-one (2k) showed minimum binding energy (-9.29 kJ/mol) with ligand efficiency of -0.31. All the ligands were docked deeply within the binding pocket region of 3FDN showing hydrogen bonds with Ala 213 and Asn 261. The docking study results showed that these derivatives are excellent inhibitor of human Aurora A kinase target; and also all these docked compounds have good inhibition constant, vdW + Hbond + desolv energy with best RMSD value.     

Molecular docking of 1H-pyrazole derivatives to receptor tyrosine kinase and protein kinase for screening potential inhibitors

Tyrosine kinase receptor and protein kinases drawn much attention for the scientific fraternity in drug discovery due to its important role in different cancer, cardiovascular diseases and other hyper-proliferative disorders. Docking studies of pyrazole derivatives with tyrosine kinase and different serine/threonine protein kinases were employed by using flexible ligand docking approach of AutoDock 4.2. Among the molecules tested for docking study, 2-(4-chlorophenyl)-5-(3-(4-chlorophenyl)-5-methyl-1- phenyl-1H-pyrazol-4-yl)-1,3,4-thiadiazole (1b), 2-(4-methoxyphenyl)-5-(3-(4-methoxyphenyl)-5-methyl-1-phenyl-1H-pyrazol-4-yl)- 1,3,4-thiadiazole (1d) and 2-(4-chlorophenyl)-5-(3-(4-chlorophenyl)-5-methyl-1-phenyl-1H-pyrazol-4-yl)-1,3,4-thiadiazole (2b) revealed minimum binding energy of -10.09, -8.57 and -10.35 kJ/mol with VEGFR-2 (2QU5), Aurora A (2W1G) and CDK2 (2VTO) protein targets, respectively. These proteins are representatives of plausible models of interactions with different anticancer agents. All the ligands were docked deeply within the binding pocket region of all the three proteins, showing reasonable hydrogen bonds. The docking study results showed that these pyrazole derivatives are potential inhibitor of all the three protein targets; and also all these docked compounds have good inhibition constant, vdW + Hbond + desolv energy with best RMSD value.     

An explorative study on Staphylococcus aureus MurE inhibitor: induced fit docking,binding free energy calculation,and molecular dynamics

Staphylococcus aureus MurE enzyme catalyzes the addition of l-lysine as third residue of the peptidoglycan peptide moiety. Due to the high substrate specificity and its ubiquitous nature among bacteria, MurE enzyme is considered as one of the potential target for the development of new therapeutic agents. In the present work, induced fit docking (IFD), binding free energy calculation, and molecular dynamics (MD) simulation were carried out to elucidate the inhibition potential of 2-thioxothiazolidin-4-one based inhibitor 1 against S. aureus MurE enzyme. The inhibitor 1 formed majority of hydrogen bonds with the central domain residues Asn151, Thr152, Ser180, Arg187, and Lys219. Binding free-energy calculation by MM-GBSA approach showed that van der Waals (ΔG_vdW, ?57.30?kcal/mol) and electrostatic solvation (ΔG_solv, ?36.86?kcal/mol) energy terms are major contributors for the inhibitor binding. Further, 30-ns MD simulation was performed to validate the stability of ligand–protein complex and also to get structural insight into mode of binding. Based on the IFD and MD simulation results, we designed four new compounds D1–D4 with promising binding affinity for the S. aureus MurE enzyme. The designed compounds were subjected to the extra-precision docking and binding free energy was calculated for complexes. Further, a 30-ns MD simulation was performed for D1/4C13 complex.     

Molecular Docking and Dynamic Simulation to Identify α7nAChR Binding Affinity of Flavonoids for the Treatment of Alzheimer's Disease

Background : Alpha-7-nicotinic acetylcholine receptor (α7nAChR), a ligand-gated ion channel is one of the important parts of the cholinergic pathway in the brain and has a remarkable role in Alzheimer's disease (AD). It has been documented that the modulation of α7nAChR with the help of phytoconstituent can be helpful in the treatment of AD. Method : The binding efficacy of fifty flavonoids was evaluated for human α7nAChR using molecular docking. The best two flavonoids shortlisted from docking analysis were then subjected to molecular dynamic simulations for 100 ns to analyze conformational binding stability with the target protein. Further, the druggability of the selected flavonoids was checked using in silico ADMET studies. Result : The top two flavonoids selected based on binding affinity toward the binding site of α7nAChR from molecular docking were amentoflavone (–9.1 kcal/mol) and gallocatechin (–8.8 kcal/mol). The molecular dynamics simulation revealed that amentoflavone and gallocatechin have a stable state during overall simulation time, lesser root mean deviation (RMSD) and root mean square fluctuation (RMSF), and complex of both compounds with protein is stable until 100 ns. Conclusion : The two flavonoids amentoflavone and gallocatechin are potential lead molecules that could be utilized as effective agonists of α7nAChR to combat Alzheimer's disease. Future in vitro and in vivo analyses are required to confirm their effectiveness.     

Study on the binding mode of a pyrrolotriazin derivative with JAK2 by docking and MD simulation

The pyrrolotriazin derivative 2-(4-(4-((7-(3-(N-methylmethylsulfonamido)phenyl)pyrrolo [2,1-f][1,2,4]triazin-2-yl)amino)phenyl)piperidin-1-yl)acetamide (PPA) is a potential Janus kinase 2 (JAK2) inhibitor. The binding mode between PPA and JAK2 was investigated by using a combined method of docking, molecular dynamics (MD) simulation and binding free-energy calculation. The docking calculations preliminarily indicated that there were two possible binding modes 1 and 2; MD simulations and binding free-energy calculations identified that binding mode 1 was more stable and favourable, with the lower MM-PBSA binding free energy of ?34.00?±?0.17?kcal/mol. Moreover, some valuable binding information is revealed as follows: the inhibitor PPA is suitably located at the ATP-binding site of JAK2 and the hydrophobic interaction plays an essential role. PPA not only interacts with residues Leu855, Val863, Ala880, Tyr931, Leu932 and Leu983 via hydrophobic interaction but also interacts with Ser936 and Asp994 by hydrogen bonds. These two factors are advantageous for PPA to strongly bind to JAK2. These results help to understand the action mechanisms and designing new compounds with a higher affinity to JAK2.     

Computational selection of flavonoid compounds as inhibitors against SARS-CoV-2 main protease,RNA-dependent RNA polymerase and spike proteins: A molecular docking study

An outbreak of Coronavirus disease 2019 (COVID-19) caused by SARS-CoV-2 has been recognized as a global health concern. Since, no specific antiviral drug is proven effective for treatment against COVID-19, identification of new therapeutics is an urgent need. In this study, flavonoid compounds were analyzed for its inhibitory potential against important protein targets of SARS-CoV-2 using computational approaches. Virtual docking was performed for screening of flavonoid compounds retrieved from PubChem against the main protease of SARS-CoV-2 using COVID-19 docking server. The cut off of dock score was set to >?9 kcal/mol and screened compounds were individually docked against main protease, RNA-dependent RNA polymerase, and spike proteins using AutoDock 4.1 software. Finally, lead flavonoid compounds were subjected to ADMET analysis. A total of 458 flavonoid compounds were virtually screened against main protease target and 36 compounds were selected based on the interaction energy value >?9 kcal/mol. Furthermore, these compounds were individually docked against protein targets and top 10 lead compounds were identified. Among the lead compounds, agathisflavone showed highest binding energy value of ?8.4 kcal/mol against main protease, Albireodelphin showed highest dock score of ?9.8 kcal/mol and ?11.2 kcal/mol against RdRp, and spike proteins, respectively. Based on the high dock score and ADMET properties, top 5 lead molecules such as Albireodelphin, Apigenin 7-(6″-malonylglucoside), Cyanidin-3-(p-coumaroyl)-rutinoside-5-glucoside, Delphinidin 3-O-beta-D-glucoside 5-O-(6-coumaroyl-beta-D-glucoside) and (-)-Maackiain-3-O-glucosyl-6″-O-malonate were identified as potent inhibitors against main protease, RdRp, and spike protein targets of SARS-CoV-2. These all compounds are having non-carcinogenic and non-mutagenic properties. This study finding suggests that the screened compounds include Albireodelphin, Apigenin 7-(6″-malonylglucoside), Cyanidin-3-(p-coumaroyl)-rutinoside-5-glucoside, Delphinidin 3-O-beta-D-glucoside 5-O-(6-coumaroyl-beta-D-glucoside) and (-)-Maackiain-3-O-glucosyl-6″-O-malonate could be the potent inhibitors of SARS-CoV-2 targets.     

Abrogation of AuroraA-TPX2 by novel natural inhibitors: molecular dynamics-based mechanistic analysis

Abstract

Introduction: Cancer is characterized by uncontrolled cell growth and genetic instabilities. The human Aurora-A kinase protein plays a crucial role in spindle assembly during mitosis and is activated by another candidate oncogene, targeting protein for Xklp2 (TPX2). It has been proposed that dissociation of Aurora A–TPX2 complex leads to disruption of mitotic spindle apparatus, thereby preventing cell division and further tumor growth. Materials and methods: A large natural compound library was docked against the active site of Aurora A–TPX2 complex. The protein–ligand complexes were subjected to molecular dynamics simulation to ascertain their binding stability. The drug properties of the compounds were analyzed to observe their drug-like properties. Results: The virtual screening of natural compound library yielded two high scoring compounds, the first compound CTOM [ZINC ID: 38143674] (Glide score: –9.49) was stable for 17 ns while the second TTOM (Glide score: ?9.07) was stable for 15 ns. While CTOM interacted with His280, Thr288 of Aurora A and Tyr34, Lys38 of TPX2, TTOM interacted with Arg285 and Arg286 in addition to the residues involved with CTOM. Conclusions: We report two natural compounds as potential drugs leads for the disruption of this complex. These ligands show a preferable docking score and have many drugs like properties within in the range of 95% of known drugs. The study provides evidence that CTOM and TTOM can efficiently inhibit the TPX2-mediated activation of Aurora A. Thus, it paves way for an elaborate investigation and establishes the importance of computational approaches as time- and cost-effective techniques.     

High throughput screening,docking, and molecular dynamics studies to identify potential inhibitors of human calcium/calmodulin-dependent protein kinase IV

Calcium/calmodulin-dependent protein kinase IV (CAMKIV) is associated with many diseases including cancer and neurodegenerative disorders and thus being considered as a potential drug target. Here, we have employed the knowledge of three-dimensional structure of CAMKIV to identify new inhibitors for possible therapeutic intervention. We have employed virtual high throughput screening of 12,500 natural compounds of Zinc database to screen the best possible inhibitors of CAMKIV. Subsequently, 40 compounds which showed significant docking scores (?11.6 to ?10.0?kcal/mol) were selected and further filtered through Lipinski rule and drug likeness parameter to get best inhibitors of CAMKIV. Docking results are indicating that ligands are binding to the hydrophobic cavity of the kinase domain of CAMKIV and forming a significant number of non-covalent interactions. Four compounds, ZINC02098378, ZINC12866674, ZINC04293413, and ZINC13403020, showing excellent binding affinity and drug likeness were subjected to molecular dynamics simulation to evaluate their mechanism of interaction and stability of protein-ligand complex. Our observations clearly suggesting that these selected ligands may be further employed for therapeutic intervention to address CAMKIV associated diseases.

Communicated by Ramaswamy H. Sarma     

In silico docking and molecular dynamics simulation of 3-dehydroquinate synthase (DHQS) from <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis>

The shikimate pathway is as an attractive target because it is present in bacteria, algae, fungi, and plants but does not occur in mammals. In Mycobacterium tuberculosis (MTB), the shikimate pathway is integral to the biosynthesis of naphthoquinones, menaquinones, and mycobactin. In these study, novel inhibitors of 3-dehydroquinate synthase (DHQS), an enzyme that catalyzes the second step of the shikimate pathway in MTB, were determined. 12,165 compounds were selected from two public databases through virtual screening and molecular docking analysis using PyRx 8.0 and Autodock 4.2, respectively. A total of 18 compounds with the best binding energies (?13.23 to ?8.22 kcal/mol) were then selected and screened for absorption, distribution, metabolism, excretion, and toxicity (ADMET) analysis, and nine of those compounds were found to satisfy all of the ADME and toxicity criteria. Among those nine, the three compounds—ZINC633887?(binding energy =??10.29 kcal/mol), ZINC08983432?(?9.34 kcal/mol), and PubChem73393?(?8.61 kcal/mol)—with the best binding energies were further selected for molecular dynamics (MD) simulation analysis. The results of the 50-ns MD simulations showed that the two compounds ZINC633887 and PubChem73393 formed stable complexes with DHQS and that the structures of those two ligands remained largely unchanged at the ligand-binding site during the simulations. These two compounds identified through docking and MD simulation are potential candidates for the treatment of TB, and should undergo validation in vivo and in vitro.     

Structure‐guided screening of chemical database to identify NS3‐NS2B inhibitors for effective therapeutic application in dengue infection

Dengue infection is the most common arthropod‐borne disease caused by dengue viruses, predominantly affecting millions of human beings annually. To find out promising chemical entities for therapeutic application in Dengue, in the current research, a multi‐step virtual screening effort was conceived to screen out the entire “screening library” of the Asinex database. Initially, through “Lipinski rule of five” filtration criterion almost 0.6 million compounds were collected and docked with NS3‐NS2B protein. Thereby, the chemical space was reduced to about 3500 compounds through the analysis of binding affinity obtained from molecular docking study in AutoDock Vina. Further, the “Virtual Screening Workflow” (VSW) utility of Schrödinger suite was used, which follows a stepwise multiple docking programs such as ‐ high‐throughput virtual screening (HTVS), standard precision (SP), and extra precision (XP) docking, and in postprocessing analysis the MM‐GBSA based free binding energy calculation. Finally, five potent molecules were proposed as potential inhibitors for the dengue NS3‐NS2B protein based on the investigation of molecular interactions map and protein‐ligand fingerprint analyses. Different pharmacokinetics and drug‐likeness parameters were also checked, which favour the potentiality of selected molecules for being drug‐like candidates. The molecular dynamics (MD) simulation analyses of protein‐ligand complexes were explained that NS3‐NS2B bound with proposed molecules quite stable in dynamic states as observed from the root means square deviation (RMSD) and root means square fluctuation (RMSF) parameters. The binding free energy was calculated using MM‐GBSA method from the MD simulation trajectories revealed that all proposed molecules possess such a strong binding affinity towards the dengue NS3‐NS2B protein. Therefore, proposed molecules may be potential chemical components for effective inhibition of dengue NS3‐NS2B protein subjected to experimental validation.     

Anti-methanogenic effect of rhubarb (Rheum spp.) – An in silico docking studies on methyl-coenzyme M reductase (MCR)

The present study explored anti-methanogenic properties of rhubarb compounds using in silico analysis on methyl-coenzyme M reductase (MCR) for identifying its anti-methanogen mechanism. To identify pharmacokinetics of 35 compounds from rhubarb, molecular docking and ADME analysis were performed against MCR using AutoDockVina, FAFDrugs3 and PROTOX programs. Docking results successfully indicated three possible candidate compounds 9,10-anthracenedione, 1,8-dihydroxy-3-methyl (?6.92 kcal/mol); phthalic acid isobutyl octadecyl ester (?5.26 kcal/mol); and diisooctyl phthalate (?5.61 kcal/mol) showed minimum binding energy (kcal/mol) with the target protein MCR which catalyze the biosynthesis of rumen methane. In conclusion, the identified compounds showed the most docking fitness score against the target methyl-coenzyme M reductase and the decrease in ruminal methane emission by rhubarb might be a result of these compounds by inhibition of methanogenesis.