Unraveling the Mechanism of Action of Myricetin in the Inhibition of hUba1∼Ubiquitin Thioester Bond Formation via In Silico Molecular Modeling Techniques.

Ubiquitination is a crucial type of protein modification which helps to control substrate degradation and maintain cell homeostasis. Recent studies suggest that ubiquitination and deubiquitination are involved in regulating metabolic reprogramming in cancer cells and maintaining cancer stem cells. Uba1, a crucial protein in the ubiquitination cascade, can be targeted to develop effective inhibitors for cancer treatment. In previous work, we showed that myricetin (Myr) acts as a potential human Uba1 (hUba1) inhibitor. In this study, we have utilized computational modeling techniques to attempt to illustrate the mechanism of action of Myr. Through extra-precision docking, we confirmed that Myr binds to the adenosine triphosphate (ATP)-binding site of hUba1 (referred to as hotspot 1) with the highest binding affinity. The dynamics of this interaction revealed that hUba1 undergoes a conformational shift from open to closed upon binding of Myr. Myr also migrates outward to interact with the crossover loop simultaneously as the rotational shift of the ubiquitin fold domain (UFD) takes place, thereby blocking access to the ubiquitin binding interface of hUba1 and the crossover loop. The outward migration also explains the reversible nature of Myr binding to hUba1 in previous experiments. We hypothesize that Myr acts as an inhibitor of Uba1∼Ub thioester bond formation by causing a large domain shift toward a closed conformation. Few other analogues of Myr containing the same flavone skeleton showed promising docking scores against hUba1 and could be considered for further validation. We propose that Myr and some of its analogues reported in this study may be promising candidates for developing effective Uba1 inhibitors for cancer treatment.

Mode of inhibitory binding of epigallocatechin gallate to the ubiquitin-activating enzyme Uba1 via accelerated molecular dynamics.

The green tea polyphenol (-)-epigallocatechin-3-gallate (EGCG) and some of its analogs potently inhibit the ubiquitin-activating enzyme Uba1. In an effort to understand the possible molecular basis of inhibitory activity of EGCG, we conducted a molecular docking and molecular dynamics simulation study. We found that EGCG and its two selected analogs, (-)-epicatechin-3-gallate (ECG) and (-)-epigallocatechin (EGC), bind favorably at two likely hot spots for small-molecule ligand binding on human Uba1. The compounds bind with energetics that mirror their experimental potency for inhibition of Uba1∼ubiquitin thioester formation. The binding of EGCG, ECG, and EGC at one of the hot spots, in particular, recapitulates the rank order of potency determined experimentally and suggests a possible mechanism for inhibition. A hinge-like conformational change of the second catalytic cysteine domain and the opposing ubiquitin-fold domain observed during accelerated molecular dynamics simulations of the EGCG-bound Uba1 complex that results in disruption of the ubiquitin-binding interfaces could explain the compounds' inhibitory activity. These results shed light on the possible molecular mechanism of EGCG and related catechins in the inhibition of Uba1.

Analysis of mutations of defensin protein using accelerated molecular dynamics simulations.

Plant defensins possess diverse biological functions that include antifungal and antibacterial activities and α-amylase and trypsin inhibitory properties. Two mutations, G9R and V39R, were confirmed to increase the antifungal activity of Raphanus sativus antifungal protein 2 (RsAFP2). Accelerated Molecular Dynamics (aMD) were carried out to examine the conformational changes present in these RsAFP2 mutants, and its two closest homologs compared to the wild-type protein. Specifically, the root mean square fluctuation values for the eight cysteine amino acids involved in the four disulfide bonds were low in the V39R mutant compared to the wild-type. Additionally, analysis of the free energy change revealed that G9R and V39R mutations exert a neutral and stabilizing effect on RsAFP2 conformation, and this is supported by the observed lower total energy of mutants compared to the wild-type, suggesting that enhanced stability of the mutants. However, MD simulations to a longer time scale would aid in capturing more conformational state of the wild-type and mutants defensin protein. Furthermore, the aMD simulations on fungal mimic membranes with RsAFP2 and its mutants and homologs showed that the mutant proteins caused higher deformation and water diffusion than the native RsAFP2, especially the V39R mutant. The mutant variants seem to interact by specifically targeting the POPC and POPI lipids amongst others. This work highlights the stabilizing effect of mutations at the 9th and 39th positions of RsAFP2 and their increased membrane deformation activity.

Natural polyphenolic inhibitors against the antiapoptotic BCL-2.

The apoptotic mechanism is regulated by the BCL-2 family of proteins, such as BCL-2 or Bcl-xL, which block apoptosis while Bad, Bak, Bax, Bid, Bim or Hrk induce apoptosis. The overexpression of BCL-2 was found to be related to the progression of cancer and also providing resistance towards chemotherapeutic treatments. In the present study, we found that all polyphenols (apigenin, fisetin, galangin and luteolin) bind to the hydrophobic groove of BCL-2 and the interaction is stable throughout MD simulation run. Luteolin was found to bind with highest negative binding energy and thus, claimed highest potency towards BCL-2 inhibition followed by fisetin. The hydrophobic interactions were found to be critical for stable complex formation as revealed by the vdW energy and ligplot analysis. Finally, on the basis of data obtained during the study, it can be concluded that these polyphenols have the potential to be used as lead molecules for BCL-2 inhibition.

Hydrophobic Interactions Are a Key to MDM2 Inhibition by Polyphenols as Revealed by Molecular Dynamics Simulations and MM/PBSA Free Energy Calculations.

p53, a tumor suppressor protein, has been proven to regulate the cell cycle, apoptosis, and DNA repair to prevent malignant transformation. MDM2 regulates activity of p53 and inhibits its binding to DNA. In the present study, we elucidated the MDM2 inhibition potential of polyphenols (Apigenin, Fisetin, Galangin and Luteolin) by MD simulation and MM/PBSA free energy calculations. All polyphenols bind to hydrophobic groove of MDM2 and the binding was found to be stable throughout MD simulation. Luteolin showed the highest negative binding free energy value of -173.80 kJ/mol followed by Fisetin with value of -172.25 kJ/mol. It was found by free energy calculations, that hydrophobic interactions (vdW energy) have major contribution in binding free energy.

Fragment based group QSAR and molecular dynamics mechanistic studies on arylthioindole derivatives targeting the α-ß interfacial site of human tubulin.

BACKGROUND: A number of microtubule disassembly blocking agents and inhibitors of tubulin polymerization have been elements of great interest in anti-cancer therapy, some of them even entering into the clinical trials. One such class of tubulin assembly inhibitors is of arylthioindole derivatives which results in effective microtubule disorganization responsible for cell apoptosis by interacting with the colchicine binding site of the ß-unit of tubulin close to the interface with the α unit. We modelled the human tubulin ß unit (chain D) protein and performed docking studies to elucidate the detailed binding mode of actions associated with their inhibition. The activity enhancing structural aspects were evaluated using a fragment-based Group QSAR (G-QSAR) model and was validated statistically to determine its robustness. A combinatorial library was generated keeping the arylthioindole moiety as the template and their activities were predicted. RESULTS: The G-QSAR model obtained was statistically significant with r2 value of 0.85, cross validated correlation coefficient q2 value of 0.71 and pred_r2 (r2 value for test set) value of 0.89. A high F test value of 65.76 suggests robustness of the model. Screening of the combinatorial library on the basis of predicted activity values yielded two compounds HPI (predicted pIC50 = 6.042) and MSI (predicted pIC50 = 6.001) whose interactions with the D chain of modelled human tubulin protein were evaluated in detail. A toxicity evaluation resulted in MSI being less toxic in comparison to HPI. CONCLUSIONS: The study provides an insight into the crucial structural requirements and the necessary chemical substitutions required for the arylthioindole moiety to exhibit enhanced inhibitory activity against human tubulin. The two reported compounds HPI and MSI showed promising anti cancer activities and thus can be considered as potent leads against cancer. The toxicity evaluation of these compounds suggests that MSI is a promising therapeutic candidate. This study provided another stepping stone in the direction of evaluating tubulin inhibition and microtubule disassembly degeneration as viable targets for development of novel therapeutics against cancer.

Group-based QSAR and molecular dynamics mechanistic analysis revealing the mode of action of novel piperidinone derived protein-protein inhibitors of p53-MDM2.

Tumour suppressor p53 is known to play a central role in prevention of tumour development, DNA repair, senescence and apoptosis which is in normal cells maintained by negative feedback regulator MDM2 (Murine Double Minute 2). In case of dysfunctioning of this regulatory loop, tumour development starts thus resulting in cancerous condition. Inhibition of p53-MDM2 binding would result in activation of the tumour suppressor. In this study, a novel robust fragment-based QSAR model has been developed for piperidinone derived compounds experimentally known to inhibit p53-MDM2 interaction. The QSAR model developed showed satisfactory statistical parameters for the experimentally reported dataset (r(2)=0.9415, q(2)=0.8958, pred_r(2)=0.8894 and F-test=112.7314), thus judging the robustness of the model. Low standard error values (r(2)_se=0.3003, q(2)_se=0.4009 and pred_r(2)_se=0.3315) confirmed the accuracy of the developed model. The regression equation obtained constituted three descriptors (R2-DeltaEpsilonA, R1-RotatableBondCount and R2-SssOCount), two of which had positive contribution while third showed negative correlation. Based on the developed QSAR model, a combinatorial library was generated and activities of the compounds were predicted. These compounds were docked with MDM2 and two top scoring compounds with binding affinities of -10.13 and -9.80kcal/mol were selected. The binding modes of actions of these complexes were analyzed using molecular dynamics simulations. Analysis of the developed fragment-based QSAR model revealed that addition of unsaturated electronegative groups at R2 site and groups with more rotatable bonds at R1 improved the inhibitory activity of these potent lead compounds. The detailed analysis carried out in this study provides a considerable basis for the design and development of novel piperidinone-based lead molecules against cancer and also provides mechanistic insights into their mode of actions.

Mechanistic insights into mode of actions of novel oligopeptidase B inhibitors for combating leishmaniasis.

Leishmaniasis is an endemic disease caused by infection with one of several different species of protozoan parasite Leishmania. Oligopeptidase B (OPB) is a serine peptidase which plays a vital role in survival of the Leishmania parasite in the host (human) macrophage and help in attaining complete virulence. Inhibition of this peptidase would check the parasite growth inside the host organism and would thus control its infection. Lack of efficient and cheap drugs has led to an urgent need for development of new anti-leishmanial drugs and this study is a step forward in this direction. Using a structure-based approach we virtually screened a large naturally-occurring compound library against OPB and subjected two top scoring compounds with high binding affinity to molecular dynamics simulations which showed a stable RMSD trajectory. The first compound COP (Glide score: -13.183) was found stable for 15 ns at RMSD of 2.5 Å while the second compound TOA (Glide score: -10.308) was stable for 8 ns at RMSD of 1.5 Å. The screened compounds interacted with some crucial residues of OPB such as COP interacted with Ser577 and His697 (part of the catalytic triad), Tyr499 (responsible for substrate stability), Arg576 (conserved in protozoan family) and Arg664 (plays a role in stabilization of the bound inhibitor). TOA also interacted with Glu669 (conserved in protozoan family) in addition to the residues interacted with COA. These interactions are crucial for OPB inhibition. This study identified naturally-occurring compound leads against OPB with good binding affinity and low toxicity to human cells.

Mechanistic insights into mode of action of novel natural cathepsin L inhibitors.

BACKGROUND: Development of a cancerous cell takes place when it ceases to respond to growth-inhibiting signals and multiplies uncontrollably and can detach and move to other parts of the body; the process called as metastasis. A particular set of cysteine proteases are very active during cancer metastasis, Cathepsins being one of them. They are involved in tumor growth and malignancy and have also been reported to be overexpressed in tumor cell lines. In the present study, a combinatorial approach comprising three-dimensional quantitative structure-activity relationship (3D QSAR), ligand-based pharmacophore modelling and search followed by cathepsin L structure-based high throughput screening was carried out using an initial set of 28 congeneric thiosemicarbazone derivatives as cathepsin L inhibitors. A 3D QSAR was derived using the alignment of a common thiosemicarbazone substructure. Essential structural features responsible for biological activity were taken into account for development of a pharmacophore model based on 29 congeneric thiosemicarbazone derivatives. This model was used to carry out an exhaustive search on a large dataset of natural compounds. A further cathepsin L structure-based screen identified two top scoring compounds as potent anti-cancer leads. RESULTS: The generated 3D QSAR model showed statistically significant results with an r(2) value of 0.8267, cross-validated correlation coefficient q(2) of 0.7232, and a pred_r(2) (r(2) value for test set) of 0.7460. Apart from these, a high F test value of 30.2078 suggested low probability of the model's failure. The pharmacophoric hypothesis chosen for searching the natural compound libraries was identified as DDHRR, where two Ds denote 2 hydrogen donors, H represents a hydrophobic group and two Rs represent aromatic rings, all of which are essential for the biological activity. We report two potential drug leads ZINC08764437 (NFP) and ZINC03846634 (APQ) obtained after a combined approach of pharmacophore-based search and structure-based virtual screen. These two compounds displayed extra precision docking scores of -7.972908 and -7.575686 respectively suggesting considerable binding affinity for cathepsin L. High activity values of 5.72 and 5.75 predicted using the 3D QSAR model further substantiated the inhibitory potential of these identified leads. CONCLUSION: The present study attempts to correlate the structural features of thiosemicarbazone group with their biological activity by development of a robust 3D QSAR model. Being statistically valid, this model provides near accurate values of the activities predicted for the congeneric set on which it is based. These predicted activities are good for the test set compounds making it indeed a statistically sound 3D QSAR model. The identified pharmacophore model DDHRR.8 comprised of all the essential features required to interact with the catalytic triad of cathepsin L. A search for natural compounds based on this pharmacophore followed by docking studies further screened out two top scoring candidates: NFP and AFQ. The high binding affinity and presence of essential structural features in these two compounds make them ideal for consideration as natural anti-tumoral agents. Activity prediction using 3D QSAR model further validated their potential as worthy drug candidates against cathepsin L for treatment of cancer.