The role of vascular endothelial protein tyrosine phosphatase on nitric oxide synthase function in diabetes: from molecular biology to the clinic

Endothelial nitric oxide synthase (eNOS) and receptor-type vascular endothelial protein tyrosine phosphatase (VE-PTP) are one of the majors signaling pathways related to endothelial health in diabetes. Several reports have shown that the inhibition of VE-PTP can lead the nitric oxide production, although repeated studies showed that VE-PTP regulated the eNOS exclusive at Ser1177 in indirect-manner. A recent, exciting paper (Siragusa et al. in Cardiovasc Res, 2020. https://doi.org/10.1093/cvr/cvaa213), showing that VE-PTP regulates eNOS in a direct-manner, dephosphorylating eNOS at Tyr81 and indirect at Ser1177 and the effects of a VE-PTP inhibitor, AKB-9778, in the blood pressure from diabetic patients.     

Structural insights into the binding mode of flavonols with the active site of matrix metalloproteinase-9 through molecular docking and molecular dynamic simulations studies

Cartilage degradation in rheumatoid arthritis is mediated principally by the collagenases and gelatinases. Gelatinase B (also called matrix metalloproteinase 9 – MMP-9), is a valid target molecule which is known to participate in cartilage degradation as well as angiogenesis associated with the disease and inhibition of its activity shall prevent cartilage damage and angiogenesis. The focus of this study is to investigate the possibilities of MMP-9 inhibition by flavonol class of bioflavonoids by studying their crucial binding interactions at the active site of MMP 9 using molecular docking (Glide XP and QPLD) and further improvisation by post-docking MM-GBSA and molecular dynamic (MD) simulations. The results show that flavonols can convincingly bind to active site of MMP-9 as demonstrated by their stable interactions at the S1′ specificity pocket and favourable binding energies. Gossypin has emerged as a promising candidate with a docking score of ?14.618 kcal/mol, binding energy of ?79.97 kcal/mol and a stable MD pattern over 15 ns. In addition, interaction mechanisms with respect to catalytic site zinc are also discussed. Further, the drug-like characters of the ligands were also analysed using ADME analysis.     

Ugi Bis-Amide Derivatives as Tyrosinase Inhibitor; Synthesis,Biology Assessment,and in Silico Analysis

Herein, a straightforward synthetic strategy mediated by Ugi reaction was developed to synthesize novel series of compounds as tyrosinase inhibitors. The structures of all compounds were confirmed by FT-IR, ¹H-NMR, ¹³C-NMR, and CHNOS techniques. The tyrosinase inhibitory activities of all synthesized derivatives 5a – m were determined against mushroom tyrosinase and it was found that derivative 5c possesses the best inhibition with an IC₅₀ value of 69.53±0.042 μM compared to the rest of the synthesized derivatives. Structure–activity relationships (SARs) showed that the presence of 4-MeO or 4-NO₂ at the R² position plays a key role in tyrosinase inhibitory activities. The enzyme kinetics studies showed that compound 5c is an noncompetitive inhibitor. For in silico study, the allosteric site detection was first applied to find the appropriate binding site and then molecular docking and molecular dynamic studies were performed to reveal the position and interactions of 5c as the most potent inhibitor within the tyrosinase active site. The results showed that 5c bind well with the proposed binding site and formed a stable complex with the target protein.     

Stoichiometric noncovalent interaction in molecular imprinting

In this review article the function of the binding site monomers in the molecular imprinting procedure is discussed. Especially, new developments towards stoichiometric noncovalent interactions are highlighted. In stoichiometric noncovalent interactions template and binding site monomer in an 1:1 molar ratio are nearly completely bound to each other. This is only possible if the association constants are considerably high (K _ass < 900 M^–1). Using this type of interaction in molecular imprinting no excess of binding sites is necessary and binding sites are only located inside the imprinted cavity. Since all cavities can be reloaded these polymers show high capacity (e.g., for preparative application) and are especially suited for the synthesis of catalytically active imprinted polymers. Discussed are binding site interactions based on amidines (and guanidines), multiple hydrogen bonding, charge-transfer interactions, and host–guest inclusion. The systematic investigation of the underlying binding reaction is described in detail. With low-molecular weight model substances the thermodynamics of the association can be conveniently investigated, e.g., by NMR spectroscopy.     

The study on interactions between levofloxacin and model proteins by using multi-spectroscopic and molecular docking methods

The interactions of levofloxacin (LEV) with lysozyme (LYZ), trypsin and bovine hemoglobin (BHb) were investigated, respectively, by using multi-spectral techniques and molecular docking in vitro. Fluorescence studies showed that LEV quenched LYZ/trypsin fluorescence in a combined quenching ways and BHb fluorescence in a static quenching with binding constants of .14, .51 and .20 × 10⁵ L mol^?1 at 298 K, respectively. The thermodynamic parameters demonstrated that hydrophobic forces, hydrogen bonds, and van der Waals forces played the major role in the binding process. The binding distances between LEV and the inner tryptophan residues of LYZ, trypsin, and BHb were calculated to be 4.04, 3.38, and 4.52 nm, respectively. Furthermore, the results of circular dichroism spectra (CD), UV–vis, and three-dimensional fluorescence spectra indicated that the secondary structures of LYZ, trypsin, and BHb were partially changed by LEV with the α-helix percentage of LYZ-LEV system increased while that of BHb-LEV system was decreased, the β-sheet percentage of trypsin-LEV system increased from 41.3 to 42.9%. UV–vis spectral results showed that the binding interactions could cause conformational and some micro-environmental changes of LYZ, trypsin, and BHb. The results of molecular docking revealed that in LYZ and trypsin systems, LEV bound to the active sites residues GLU 35 and ASP 52 of LYZ and trypsin at the active site SER 195, and in BHb system, LEV was located in the central cavity, which was consistent with the results of synchronous fluorescence experiment. Besides, LEV made the activity of LYZ decrease while the activity of trypsin increased.     

Hydrogen bonding-assisted interaction between amitriptyline hydrochloride and hemoglobin: spectroscopic and molecular dynamics studies

Herein, we have explored the interaction between amitriptyline hydrochloride (AMT) and hemoglobin (Hb), using steady-state and time-resolved fluorescence spectroscopy, UV–visible spectroscopy, and circular dichroism spectroscopy, in combination with molecular docking and molecular dynamic (MD) simulation methods. The steady-state fluorescence reveals the static quenching mechanism in the interaction system, which was further confirmed by UV–visible and time-resolved fluorescence spectroscopy. The binding constant, number of binding sites, and thermodynamic parameters viz. ΔG, ΔH, ΔS are also considered; result confirms that the binding of the AMT with Hb is a spontaneous process, involving hydrogen bonding and van der Waals interactions with a single binding site, as also confirmed by molecular docking study. Synchronous fluorescence, CD data, and MD simulation results contribute toward understanding the effect of AMT on Hb to interpret the conformational change in Hb upon binding in aqueous solution.     

PDZ binding to the BAR domain of PICK1 is elucidated by coarse-grained molecular dynamics

A key regulator of α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid (AMPA) receptor traffic, PICK1 is known to interact with over 40 other proteins, including receptors, transporters and ionic channels, and to be active mostly as a homodimer. The current lack of a complete PICK1 structure determined at atomic resolution hinders the elucidation of its functional mechanisms. Here, we identify interactions between the component PDZ and BAR domains of PICK1 by calculating possible binding sites for the PDZ domain of PICK1 (PICK1-PDZ) to the homology-modeled, crescent-shaped dimer of the PICK1-BAR domain using multiplexed replica-exchange molecular dynamics (MREMD) and canonical molecular dynamics simulations with the coarse-grained UNRES force field. The MREMD results show that the preferred binding site for the single PDZ domain is the concave cavity of the BAR dimer. A second possible binding site is near the N-terminus of the BAR domain that is linked directly to the PDZ domain. Subsequent short canonical molecular dynamics simulations used to determine how the PICK1-PDZ domain moves to the preferred binding site on the BAR domain of PICK1 revealed that initial hydrophobic interactions drive the progress of the simulated binding. Thus, the concave face of the BAR dimer accommodates the PDZ domain first by weak hydrophobic interactions and then the PDZ domain slides to the center of the concave face, where more favorable hydrophobic interactions take over.     

Insights into the interaction of methotrexate and human serum albumin: A spectroscopic and molecular modeling approach

In this study, fluorescence spectroscopy and molecular modeling approaches were employed to investigate the binding of methotrexate to human serum albumin (HSA) under physiological conditions. From the mechanism, it was demonstrated that fluorescence quenching of HSA by methotrexate results from the formation of a methotrexate/HSA complex. Binding parameters calculated using the Stern–Volmer method and the Scatchard method showed that methotrexate binds to HSA with binding affinities in the order 10⁴ L·mol^?1. Thermodynamic parameter studies revealed that the binding reaction is spontaneous, and that hydrogen bonds and van der Waals interactions play a major role in the reaction. Site marker competitive displacement experiments and a molecular modeling approach demonstrated that methotrexate binds with appropriate affinity to site I (subdomain IIA) of HSA. Furthermore, we discuss some factors that influence methotrexate binding to HSA.     

Identification of the hot spot residues for pyridine derivative inhibitor CCT251455 and ATP substrate binding on monopolar spindle 1 (MPS1) kinase by molecular dynamic simulation

Protein kinase monopolar spindle 1 plays an important role in spindle assembly checkpoint at the onset of mitosis. Over expression of MPS1 correlated with a wide range of human tumors makes it an attractive target for finding an effective and specific inhibitor. In this work, we performed molecular dynamics simulations of protein MPS1 itself as well as protein bound systems with the inhibitor and natural substrate based on crystal structures. The reported orally bioavailable 1 h-pyrrolo [3,2-c] pyridine inhibitors of MPS1 maintained stable binding in the catalytic site, while natural substrate ATP could not stay. Comparative study of stability and flexibility of three systems reveals position shifting of β-sheet region within the catalytic site, which indicates inhibition mechanism was through stabilizing the β-sheet region. Binding free energies calculated with MM-GB/PBSA method shows different binding affinity for inhibitor and ATP. Finally, interactions between protein and inhibitor during molecular dynamic simulations were measured and counted. Residue Gly605 and Leu654 were suggested as important hot spots for stable binding of inhibitor by molecular dynamic simulation. Our results reveal an important position shifting within catalytic site for non-inhibited proteins. Together with hot spots found by molecular dynamic simulation, the results provide important information of inhibition mechanism and will be referenced for designing novel inhibitors.     

Elucidating the interaction of clofazimine with bovine liver catalase; a comprehensive spectroscopic and molecular docking approach

Nowadays, understanding of interface between protein and drugs has become an active research area of interest. These types of interactions provide structural guidelines in drug design with greater clinical efficacy. Thus, structural changes in catalase induced by clofazimine were monitored by various biophysical techniques including UV‐visible spectrometer, fluorescence spectroscopy, circular dichroism, and dynamic light scattering techniques. Increase in absorption spectra (UV‐visible spectrum) confers the complex formation between drug and protein. Fluorescence quenching with a binding constants of 2.47 × 10⁴ M⁻¹ revealed that clofazimine binds with protein. Using fluorescence resonance energy transfer, the distance (r ) between the protein (donor) and drug (acceptor) was found to be 2.89 nm. Negative Gibbs free energy change (ΔG °) revealed that binding process is spontaneous. In addition, an increase in α‐helicity was observed by far‐UV circular dichroism spectra by adding clofazimine to protein. Dynamic light scattering results indicate that topology of bovine liver catalase was slightly altered in the presence of clofazimine. Hydrophobic interactions are the main forces between clofazimine and catalase interaction as depicted by molecular docking studies. Apart from hydrophobic interactions, some hydrogen bonding was also observed during docking method. The results obtained from the present study may establish abundant in optimizing the properties of ligand‐protein mixtures relevant for numerous formulations.     

Insights into the interaction of ulipristal acetate and human serum albumin using multi-spectroscopic methods,molecular docking,and dynamic simulation

Interaction between ulipristal acetate (UPA) and human serum albumin (HSA) was investigated in simulated physiological environment using multi-spectroscopic and computational methods. Fluorescence experiments showed that the quenching mechanism was static quenching, which was confirmed by the time-resolved fluorescence. Binding constants (K_a) were found to be 1?×?10⁵ L mol^?1, and fluorescence data showed one binding site. Thermodynamic constants suggested the binding process was mainly controlled by electrostatic interactions. Results from the competition experiments indicated that UPA bound to site I of HSA. Fourier transform infrared spectra, circular dichroism spectra, synchronous fluorescence spectra, and 3D fluorescence indicated that UPA can induce conformation change in the HSA. The content of α-helix and β-sheet increased, while β-turn decreased. Hydrophobicity around the tryptophan residues declined, whereas its polarity increased. Molecular docking results were consistent with the experimental results. Results suggested that UPA located at the hydrophobic cavity site I of HSA, and hydrophobic force played the key role in the binding process. Moreover, molecular dynamics simulation was performed to determine the stability of free HSA and HSA-UPA system. Results indicated that UPA can stabilize HSA to a certain degree and enhance the flexibility of residues around site I.

Communicated by Ramaswamy H. Sarma     

Probing the acting mode and advantages of RMC-4550 as an Src-homology 2 domain-containing protein tyrosine phosphatase (SHP2) inhibitor at molecular level through molecular docking and molecular dynamics

Abstract

The over-activation of Ras/mitogen-activated protein kinase (MAPK) signaling pathway associated with a variety of cancers is usually related with abnormal activation of Src-homology 2 domain-containing protein tyrosine phosphatase (SHP2). For this purpose, SHP2 has attracted extensive interest as a potential target for cancer treatment. RMC-4550, as a newly developed selective inhibitor of SHP2, possesses an overwhelming advantage over the previous generation inhibitor SHP099 in terms of in vitro activity. However, the binding mode of SHP2 with RMC-4550 and the reason for the high efficiency of RMC-4550 as SHP2 inhibitor at molecular level are still unclear. Therefore, in this study, the binding mode of RMC-4550 with SHP2 and the superiorities of RMC-4550 as inhibitor at binding affinity and dynamic interactive behavior with SHP2 were probed by molecular docking and molecular dynamics (MD) simulations. By comparing the results of molecular docking, it was found that SHP2 formed more tight interaction with RMC-4550 than that with SHP099. Subsequently, a series of post-dynamic analyses on three simulation trajectories (SHP2^WT, SHP2^SHP099 and SHP2^RMC-4550) were performed and found that the SHP2 protein bound with RMC-4550 maintained a firmer interaction between N-Src-homology 2 (N-SH2) and PTP domain throughout the MD simulation, leading to a more stable protein conformation. The finding here provides new clues for the design of SHP2 inhibitor against the over-activation of Ras/MAPK pathway.

Communicated by Ramaswamy H. Sarma     

Evaluating the suitability of RNA intervention mechanism exerted by some flavonoid molecules against dengue virus MTase RNA capping site: a molecular docking,molecular dynamics simulation,and binding free energy study

Abstract

Dengue virus (DENV) is one of the most dangerous mosquito-borne human pathogens known to the mankind. Currently, no vaccines or standard therapy is avaliable to treate DENV infection. This makes the drug development against DENV more significant and challenging. The MTase domain of DENV RNA RdRp NS5 is a promising drug target, because this domain hosts the RNA capping process of DENV RNA to escape from human immune system. In the present study, we have analysed the RNA intervention mechanism exerted by flavoniod molecules against NS5 MTase RNA capping site by using molecular docking, molecular dynamics simulation and the binding free energy calculations. The results from the docking analysis confirmed that the RNA intervention mecanism is exerted by the quercetagetin (QGN) molecule with all necessary intermolecular interactions and high binding affinity. Notably, QGN forms strong hydrogen bonding interactions with Asn18, Leu20 and Ser150 residues and π???π stacking interaction with Phe25 residue. The apo and QGN bound NS5 MTase and QGN-NS5 MTase complex were used for MD simulation. The results of MD simulation reveal that the RMSD and RMSF values of QGN-MTase complex have increased on comparing the apo protein due to the effect of ligand binding. The binding free energy calulation includes prediction of total binding free energy of ligand-protein complex and per-residue free energy decomposition. The QGN binding to NS5 MTase affects it’s native motion, this result is found from Principal component analysis.

Communicated by Ramaswamy H. Sarma     

Underlying molecular interaction of bovine serum albumin and linezolid: a biophysical outlook

Linezolid, one of the reserve antibiotic of oxazolidinone class has wide range of antimicrobial activity. Here we have conducted a fundamental study concerning the dynamics of its interaction with bovine serum albumin (BSA), and the post binding modification of the later by employing different spectroscopic (absorption, fluorescence and circular dichroism (CD) spectroscopy) and molecular docking tools. Gradual quenching of the tryptophan (Trp) fluorescence upon addition of linezolid to BSA confirms their interaction. Analysis of fluorescence quenching at different temperature indicates that the interaction is made by static complex formation and the BSA has one binding site for the drug. The negative Gibbs energy change (ΔG⁰), and positive values of enthalpy change (ΔH⁰) and entropy change (ΔS⁰) strongly suggest that it is an entropy driven spontaneous and endothermic reaction. The reaction involves hydrophobic pocket of the protein, which is further stabilized by hydrogen bonding and electrostatic interactions as evidenced from 8-anilino-1-napthalene sulfonic acid, sucrose and NaCl binding studies. These findings also support the molecular docking study using AutoDock 4.2. The influence of this interaction on the secondary structure of the protein is negligible as evidenced by CD spectroscopy. So, from these findings, we conclude that linezolid interacts with BSA in 1:1 ratio through hydrophobic, hydrogen bonding and ionic interactions, and this may not affect the secondary structure of the protein.     

Understanding the conformational flexibility and electrostatic properties of curcumin in the active site of rhAChE via molecular docking,molecular dynamics,and charge density analysis

Acetylcholinesterase (AChE) is an important enzyme responsible for Alzheimer’s disease, as per report, keto-enol form of curcumin inhibits this enzyme. The present study aims to understand the binding mechanism of keto-enol curcumin with the recombinant human Acetylcholinesterase (rhAChE) from its conformational flexibility, intermolecular interactions, charge density distribution, and the electrostatic properties at the active site of rhAChE. To accomplish this, a molecular docking analysis of curcumin with the rhAChE was performed, which gives the structure and conformation of curcumin in the active site of rhAChE. Further, the charge density distribution and the electrostatic properties of curcumin molecule (lifted from the active site of rhAChE) were determined from the high level density functional theory (DFT) calculations coupled with the charge density analysis. On the other hand, the curcumin molecule was optimized (gas phase) using DFT method and further, the structure and charge density analysis were also carried out. On comparing the conformation, charge density distribution and the electrostatic potential of the active site form of curcumin with the corresponding gas phase form reveals that the above said properties are significantly altered when curcumin is present in the active site of rhAChE. The conformational stability and the interaction of curcumin in the active site are also studied using molecular dynamics simulation, which shows a large variation in the conformational geometry of curcumin as well as the intermolecular interactions.     

Insight into the structural requirements of thiophene-3-carbonitriles-based MurF inhibitors by 3D-QSAR,molecular docking and molecular dynamics study

The discovery of clinically relevant inhibitors against MurF enzyme has proven to be a challenging task. In order to get further insight into the structural features required for the MurF inhibitory activity, we performed pharmacophore and atom-based three-dimensional quantitative structure–activity relationship studies for novel thiophene-3-carbonitriles based MurF inhibitors. The five-feature pharmacophore model was generated using 48 inhibitors having IC₅₀ values ranging from 0.18 to 663?μm. The best-fitted model showed a higher coefficient of determination (R²?=?0.978), cross-validation coefficient (Q²?=?0.8835) and Pearson coefficient (0.9406) at four component partial least-squares factor. The model was validated with external data set and enrichment study. The effectiveness of the docking protocol was validated by docking the co-crystallized ligand into the catalytic pocket of MurF enzyme. Further, binding free energy calculated by the molecular mechanics generalized Born surface area approach showed that van der Waals and non-polar solvation energy terms are the main contributors to ligand binding in the active site of MurF enzyme. A 10-ns molecular dynamic simulation was performed to confirm the stability of the 3ZM6-ligand complex. Four new molecules are also designed as potent MurF inhibitors. These results provide insights regarding the development of novel MurF inhibitors with better binding affinity.     

Molecular modeling study for conformational changes of Sirtuin 2 due to substrate and inhibitor binding

Sirtuin is a member of NAD⁺-dependent deacetylase family. The structural details of Sirtuin 2 (SIRT2) complex will be very useful to discover the drug which might have beneficial effects on various diseases like cancer, diabetes, etc. Unfortunately, SIRT2 complex structure is not available yet, hence molecular docking was carried out to dock the substrate (NAD⁺ and acetylated lysine) and inhibitor (sirtinol) in the NAD⁺ binding site. The suitable binding orientation of substrate and inhibitor in the SIRT2 active site was selected and subjected to 5?ns molecular dynamics simulations to adjust the binding orientation of inhibitor and substrate as well as to identify the conformational changes in the active site. The result provides an insight about 3D SIRT2 structural details as well as the importance of F96 in deacetylation function. In addition, our simulations revealed the displacement of F96 upon substrate and inhibitor binding, inducing an extended conformation of loop3 and changing its interactions with the rest of SIRT2. We believe that our study could be helpful to gain a structural insight of SIRT2 and to design the receptor-based inhibitors.     

Probing the behavior of bovine serum albumin upon binding to atenolol: insights from spectroscopic and molecular docking approaches

Molecular interaction of atenolol, a selective β₁ receptor antagonist with the major carrier protein, bovine serum albumin (BSA), was investigated under imitated physiological conditions (pH 7.4) by means of fluorescence spectroscopy, UV absorption spectroscopy, Fourier transform infrared spectroscopy (FT-IR), and molecular modeling studies. The steady-state fluorescence spectra manifested that static type, due to formation of the atenolol-BSA complex, was the dominant mechanism for fluorescence quenching. The characteristic information about the binding interaction of atenolol with BSA in terms of binding constant (K_b) were determined by the UV–vis absorption titration, and were found to be in the order of 10³ M^?1 at different temperatures, indicating the existence of a weak binding in this system. Thermodynamic analysis revealed that the binding process was primarily mediated by van der Waals force and hydrogen bonds due to the negative sign for enthalpy change (ΔH⁰), entropy change (ΔS⁰). The molecular docking results elucidated that atenolol preferred binding on the site II of BSA according to the findings observed in competitive binding experiments. Moreover, via alterations in synchronous fluorescence, three-dimensional fluorescence and FT-IR spectral properties, it was concluded that atenolol could arouse slight configurational and micro-environmental changes of BSA.     

Insight into the key active sites of afChiA1 based on molecular dynamics simulations and free energy calculations

In this study, the binding of the enzyme chitinase A1 (afChiA1) from the plant-type Aspergillus fumigatus with four potent inhibitors, allosamidin (ASM), acetazolamide (AZM), 8-chloro-theophylline (CTP) and kinetin (KIT) is investigated by molecular docking, molecular dynamics simulation and binding free energy calculation. The results reveal that the electrostatic interactions play an important role in the stabilisation of the binding of afChiA1 with inhibitors. Based on the binding energy of afChiA1-ligands, the key residues (Gln37 and Trp312) in the active binding pocket of the complex systems are confirmed by molecular mechanics/Poisson–Boltzmann surface area method, and the active inhibitors, ASM and AZM, both could form strong interaction with Gln37 and Trp312, and the non-active ligands, CTP and KIT, could not interact with these two residues, which is consistent with the result of experimental report. Then, it is identified that Gln37 and Trp312 should be one of the important active site residues of afChiA1.     

Binding mechanism of spinosine and venenatine molecules with p300 HAT enzyme: Molecular screening,molecular dynamics and free-energy analysis

The chromatin modification is regulated by the histone acetyltransferase (HAT) and histone deacetyltransferase (HDAC) enzymes; abnormal function of these enzymes leads to several malignant diseases. The inhibition of these enzymes using natural ligand molecules is an emerging technique to cure these diseases. The in vitro analysis of natural molecules, venenatine, spinosine, palmatine and taxodione are giving the best inhibition rate against p300 HAT enzyme. However, the detailed understanding of binding and the stability of these molecules with p300 HAT is not yet known. The aim of the present study is focused to determine the binding strength of the molecules from molecular dynamics simulation analysis. The docking analysis confirms that, the venenatine (−6.97 kcal/mol - conformer 8), spinosine (−6.52 kcal/mol conformer −10), palmatine (−5.72 kcal/mol conformer-3) and taxodione (−4.99 kcal/mol conformer-4) molecules form strong hydrogen bonding interactions with the key amino acid residues (Arg1410, Thr1411 and Trp1466) present in the active site of p300. In the molecular dynamics (MD) simulation, the spinosine retain these key interactions with the active site amino acid residues (Arg1410, Thr1411, and Trp1466) than venenatine and are stable throughout the simulation. The RMSD value of spinosine (0.5 to 1.3 Å) and venenatine (0.3 to 1.3 Å) are almost equal during the MD simulation. However, during the MD simulation, the intermolecular interaction between venenatine and the active site amino acid residues (Arg1410, Thr1411, and Trp1466) decreased on comparing with the spinosine-p300 interaction. The binding free energy of the spinosine (−15.30 kcal/mol) is relatively higher than the venenatine (−11.8 kcal/mol); this increment is attributed to the strong hydrogen bonding interactions of spinosine molecule with the active site amino acid residues of p300.