Identification of potent BRD4-BD1 inhibitors using classical and steered molecular dynamics based free energy analysis.

In the present work a combination of traditional and steered molecular dynamics based techniques were employed to identify potential inhibitors against the human BRD4 protein (BRD4- BD1); an established drug target for multiple illnesses including various malignancies. Quinoline derivatives that were synthesized in-house were tested for their potential as new BRD4-BD1 inhibitors. Initially molecular docking experiments were performed to determine the binding poses of BRD4-BD1 inhibitors. To learn more about the thermodynamics of inhibitor binding to the BRD4-BD1 active site, the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) free energy calculations were conducted afterwards. The findings of the MM-PBSA analysis were further reinforced by performing steered umbrella sampling simulations which revealed crucial details about the binding/unbinding process of the most potent quinoline derivatives at the BRD4-BD1 active site. We report a novel quinoline derivative which can be developed into a fully functional BRD4-BD1 inhibitor after experimental validation. The identified compound (4 g) shows better properties than the standard BRD4-BD1 inhibitors considered in the study. The study also highlights the crucial role of Gln78, Phe79, Trp81, Pro82, Phe83, Gln84, Gln85, Val87, Leu92, Leu94, Tyr97, Met105, Cys136, Asn140, Ile146 and Met149 in inhibitor binding. The study provides a possible lead candidate and key amino acids involved in inhibitor recognition and binding at the active site of BRD4-BD1 protein. The findings might be of significance to medicinal chemists involved in the development of potent BRD4-BD1 inhibitors.

Immuno-informatics analysis predicts B and T cell consensus epitopes for designing peptide vaccine against SARS-CoV-2 with 99.82% global population coverage.

The current global pandemic due to Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has taken a substantial number of lives across the world. Although few vaccines have been rolled-out, a number of vaccine candidates are still under clinical trials at various pharmaceutical companies and laboratories around the world. Considering the intrinsic nature of viruses in mutating and evolving over time, persistent efforts are needed to develop better vaccine candidates. In this study, various immuno-informatics tools and bioinformatics databases were deployed to derive consensus B-cell and T-cell epitope sequences of SARS-CoV-2 spike glycoprotein. This approach has identified four potential epitopes which have the capability to initiate both antibody and cell-mediated immune responses, are non-allergenic and do not trigger autoimmunity. These peptide sequences were also evaluated to show 99.82% of global population coverage based on the genotypic frequencies of HLA binding alleles for both MHC class-I and class-II and are unique for SARS-CoV-2 isolated from human as a host species. Epitope number 2 alone had a global population coverage of 98.2%. Therefore, we further validated binding and interaction of its constituent T-cell epitopes with their corresponding HLA proteins using molecular docking and molecular dynamics simulation experiments, followed by binding free energy calculations with molecular mechanics Poisson-Boltzmann surface area, essential dynamics analysis and free energy landscape analysis. The immuno-informatics pipeline described and the candidate epitopes discovered herein could have significant impact upon efforts to develop globally effective SARS-CoV-2 vaccines.

A lesson for the maestro of the replication fork: Targeting the protein-binding interface of proliferating cell nuclear antigen for anticancer therapy.

The proliferating cell nuclear antigen (PCNA) has emerged as a promising candidate for the development of novel cancer therapeutics. PCNA is a nononcogenic mediator of DNA replication that regulates a diverse range of cellular functions and pathways through a comprehensive list of protein-protein interactions. The hydrophobic binding pocket on PCNA offers an opportunity for the development of inhibitors to target various types of cancers and modulate protein-protein interactions. In the present study, we explored the binding modes and affinity of molecule I1 (standard molecule) with the previously suggested dimer interface pocket and the hydrophobic pocket present on the frontal side of the PCNA monomer. We also identified potential lead molecules from the library of in-house synthesized 3-methylenisoindolin-1-one based molecules to inhibit the protein-protein interactions of PCNA. Our results were based on robust computational methods, including molecular docking, conventional, steered, and umbrella sampling molecular dynamics simulations. Our results suggested that the standard inhibitor I1 interacts with the hydrophobic pocket of PCNA with a higher affinity than the previously suggested binding site. Also, the proposed molecules showed better or comparable binding free energies as calculated by the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach and further validated by enhanced umbrella sampling simulations. In vitro and in vivo methods could test the computationally suggested molecules for advancement in the drug discovery pipeline.

Identification of acridinedione scaffolds as potential inhibitor of DENV-2 C protein: An in silico strategy to combat dengue.

Dengue is a prominent viral disease transmitted by mosquitoes to humans that affects mainly tropical and subtropical countries worldwide. The global spread of dengue virus (DENV) is mainly occurred by Aedes aegypti and Aedes albopictus mosquitoes. The dengue virus serotypes-2 (DENV-2) is a widely prevalent serotype of DENV, that causes the hemorrhagic fever and bleeding in the mucosa, which can be fatal. In the life cycle of DENV-2, a structural capsid (DENV-2 C) protein forms the nucleocapsid assembly and bind to the viral progeny RNA. For DENV-2 maturation, the nucleocapsid is a vital component. We used virtual ligand screening to filter out the best in-house synthesized acridinedione analogs (DSPD molecules) that could efficiently bind to DENV-2 C protein. The molecular docking and dynamics simulations studies were performed to analyze the effect of DSPD molecules on DENV-2 C protein after binding. Our findings showed that DSPD molecules strongly interacted with DENV-2 C protein, as evident from molecular interactions and several time-dependent molecular dynamics-driven analyses. Moreover, this study was also supported by the thermodynamic binding free energy and steered molecular dynamics simulations. Therefore, we intend to suggest that the DSPD3 molecule could be used as a potential therapeutic molecule against dengue complications as compared to the cocrystallized inhibitor ST-148. However, further studies are required to demonstrate the ability of DSPD3 to induce DENV-2 C tetramer formation.

Inhibition of nonstructural protein 15 of SARS-CoV-2 by golden spice: A computational insight.

The quick widespread of the coronavirus and speedy upsurge in the tally of cases demand the fast development of effective drugs. The uridine-directed endoribonuclease activity of nonstructural protein 15 (Nsp15) of the coronavirus is responsible for the invasion of the host immune system. Therefore, developing potential inhibitors against Nsp15 is a promising strategy. In this concern, the in silico approach can play a significant role, as it is fast and cost-effective in comparison to the trial and error approaches of experimental investigations. In this study, six turmeric derivatives (curcuminoids) were chosen for in silico analysis. The molecular interactions, pharmacokinetics, and drug-likeness of all the curcuminoids were measured. Further, the stability of Nsp15-curcuminoids complexes was appraised by employing molecular dynamics (MD) simulations and MM-PBSA approaches. All the molecules were affirmed to have strong interactions and pharmacokinetic profile. The MD simulations data stated that the Nsp15-curcuminoids complexes were stable during simulations. All the curcuminoids showed stable and high binding affinity, and these curcuminoids could be admitted as potential modulators for Nsp15 inhibition.

Computational targeting of allosteric site of MEK1 by quinoline-based molecules.

MEK1 is an attractive target due to its role in selective extracellular-signal-regulated kinase phosphorylation, which plays a pivotal role in regulating cell proliferation. Another benefit of targeting the MEK protein is its unique hydrophobic pocket that can accommodate highly selective allosteric inhibitors. To date, various MEK1 inhibitors have reached clinical trials against several cancers, but they were discarded due to their severe toxicity and low efficacy. Thus, the development of allosteric inhibitors for MEK1 is the demand of the hour. In this in-silico study, molecular docking, long-term molecular dynamics (5 µs), and molecular mechanics Poisson-Boltzmann surface area analysis were undertaken to address the potential of quinolines as allosteric inhibitors. We selected four reference MEK1 inhibitors for the comparative analysis. The drug-likeness and toxicity of these molecules were also examined based on their ADMET and Toxicity Prediction by Komputer Assisted Technology profiles. The outcome of the analysis revealed that the quinolines (4m, 4o, 4s, and 4n) exhibited better stability and binding affinity while being nontoxic compared to reference inhibitors. We have reached the conclusion that these quinoline molecules could be checked by experimental studies to validate their use as allosteric inhibitors against MEK1.

Discovery and in silico evaluation of aminoarylbenzosuberene molecules as novel checkpoint kinase 1 inhibitor determinants.

Checkpoint kinase 1 (CHK1) is an essential kinase with a critical function in cell cycle arrest. Several potent inhibitors targeting CHK1 have been published, but most of them have failed in clinical trials. Acknowledging the emerging consequence of CHK1 inhibitors in medication of cancer, there is a demand for widening the chemical range of CHK1 inhibitors. In this research, we considered a set of in-house plant based semi-synthetic aminoarylbenzosuberene molecules as potential CHK1 inhibitors. Based on a combined computational research that consolidates molecular docking and binding free energy computations we recognized the crucial determinants for their receptor binding. The drug likeness of these molecules were also scrutinized based on their toxicity and bioavailibilty profile. The computational strategy indicates that the Bch10 could be regarded as a potential CHK1 inhibitor in comparison with top five co-crystallize molecules. Bch10 signifies a promising outlet for the development of potent inhibitors for CHK1.

Mechanochemical Approach towards Multi-Functionalized 1,2,3-Triazoles and Anti-Seizure Drug Rufinamide Analogs Using Copper Beads.

Highly regiospecific, copper-salt-free and neat conditions have been demonstrated for the 1,3-dipolar azide-alkyne cycloaddition (AAC) reactions under mechanochemical conditions. A group of structurally challenging alkynes and heterocyclic derivatives was efficiently implemented to achieve highly functionalized 1,4-disubstituted-1,2,3-triazoles in good to excellent yield by using the Cu beads without generation of unwanted byproducts. Furthermore, the high-speed ball milling (HSBM) strategy has also been extended to the synthesis of the commercially available pharmaceutical agent, Rufinamide, an antiepileptic drug (AED) and its analogues. The same strategy was also applied for the synthesis of the Cl-derivative of Rufinamide. Analysis of the single crystal XRD data of the triazole was also performed for the final structural confirmation. The Cu beads are easily recoverable from the reaction mixture and used for the further reactions without any special treatment.

Benzosuberene-sulfone analogues synthesis from Cedrus deodara oil and their therapeutic evaluation by computational analysis to treat type 2 diabetes.

Benzosuberene-sulfone (BSS) analogues have been semi-synthesized following green approaches from himachalenes, which has been extracted from essential oil of Cedrus deodara. In this process, benzosuberene in presence of different aryl or alkyl sodium sulfinates, I2 and potassium persulfate (K2S2O8) in acetonitrile-water solvent conditions gave BSS-analogues at room temperature. Under this reaction, a facile endocyclic ß-H elimination has been noticed for BSS-analogues synthesis instead of vinyl sulfones and the reason may be due to its specific structure and electronic environment. The BSS-compounds were obtained with moderate to excellent yields under mild conditions. All the compounds were computationally subjected to drug likeliness and toxicity prediction studies. Further, the synthesized molecules were evaluated under in-silico studies for their binding affinity towards the native Peroxisome Proliferator-Activated Receptor Gamma (PPARG), and two PPARG mutants (R357A and V290M). Both the mutant forms of PPARG are deficient in eliciting a response to treatment with full and partial agonists. Our computational studies suggested that the molecule 3q performed better than the standard drug (Rosiglitazone) in all three protein structures. This implies that our suggested molecule could act as a more potent antagonist to native PPARG and could also be developed to treat type-2 diabetes patients with R357A and V290M mutations, which didn't elicit any response to currently available drugs in the market.

Plant-based analogues identified as potential inhibitor against tobacco mosaic virus: A biosimulation approach.

Benzosuberene compounds with a pyrrolone group adhered to it are compounds extracted from the oils of Cedrus deodara plant, that bear inhibitory capabilities. Tobacco mosaic virus is known to affect crop production every year. The currently known inhibitors against TMV have a weak inhibition effect and also tend to be toxic towards non-target living organisms as well as the environment. Thus, the requirement of non-toxic potent inhibitors is the need of the hour, which led us to test our benzosuberene molecules on the binding site of TMV and check their affinity as well as stability. The non-toxic nature of these molecules has already been experimentally established. Through in-silico analysis involving docking and simulation experiments, we compared the interaction pattern of these ligand molecules with the already present inhibitors. Our investigation proved that the reported ligands (ligands 3, 7, 9, and 17 obtained -177.103, -228.632, -184.134, and - 188.075 kJ/mol binding energies, respectively) interacted with the binding site of TMV much efficiently than the known inhibitors (Ribavirin and Zhao et al. 2020 obtained 121.561 and - 221.393 kJ/mol binding energies, respectively). Moreover, they acquired a stable conformation inside the binding pocket, where a higher number of binding site residues contributed towards interaction. Thus, their structural framework can be optimized for the exploration of their antiviral properties to develop potent botanical viricides against plant virus infection.

Structural changes induced by substitution of amino acid 129 in the coat protein of Cucumber mosaic virus.

Cucumber mosaic virus infection leads to mosaic symptoms on a broad range of crop plants. Mutation at positions 129 in the coat protein of virus causes alterations in the severity of symptoms caused by the viral infection. In our investigation, we performed long term molecular dynamics simulations to elucidate the effect of different amino acid substitutes (infectious and non-infectious) at position 129 in the coat protein of Cucumber mosaic virus using various structural parameters. We found that the contagious mutants displayed more flexibility at loops ßE-αEF (129-136) and ßF-ßG loop (155-163) as compared to the non-infectious and native structures. This specific study at the atomic level yields innovative ideas for designing new therapeutic agents against the pathogen, which would further pave the path for researchers to control this devastating plant virus.

Elimination of bitter-off taste of stevioside through structure modification and computational interventions.

Stevioside is a natural non-caloric sweetener obtained from Stevia rebaudiana Bertoni plant. The major challenge in the commercialization of stevioside as a natural sweetener is its bitter-off taste. In this study, we prepared molecular models of potential taste receptors of stevioside, both sweet and bitter. With appropriate modifications on the stevioside backbone, we performed molecular docking of prepared ligands with both sweet and bitter taste receptors. Based on binding energy, we found that one of the potential substituents, Kamiya-8, shows a good affinity towards sweet taste and a weak affinity for bitter receptors. Further, we selected Kamiya-8 for molecular dynamics simulations to improve the prediction of binding energy and to check the binding strength of Kamiya-8 with taste receptors. Moreover, we also performed MM-PBSA calculation for calculating the end state free energies of molecules in solvent and found that Kamiya-8 gives a 2-fold effect as it interacts with sweet receptors (T1R2, T1R3) with lowest binding energy conformation (-285.265 kcal/mol, -571.481 kcal/mol). Secondly, it gives high binding energy (-273.319 kcal/mol, -355.500 kcal/mol) with bitter taste receptors (T2R4, T2R14) as compared to stevioside. Based on this study, we found that Kamiya-8 can be the potential substituent that can improve the palatability of stevioside.

Iodine(iii) promoted ring-rearrangement reaction of 1-arylamino-2-oxocyclopentane-1-carbonitriles to synthesize N-aryl-δ-valerolactams.

The first example of N-aryl-δ-valerolactam synthesis via an intramolecular ring rearrangement reaction of 1-arylamino-2-oxocyclopentane-1-carbonitrile promoted by phenyliodine bis(trifluoroacetate) (PIFA) was reported. We show that this unprecedented regio-selective ring-rearrangement reaction driven by hypervalent iodine (PIFA) involves C5-H elimination, C1-C2 bond opening, and C1-N bond rearrangement steps and restraining the leaving tendency of the CN group. The structure of the lactam was further confirmed by single crystal X-ray diffraction (XRD) analysis. The present protocol showed a diverse array of functional group tolerance under the reaction conditions and offered good to excellent yields of lactams.

Phloretin and phloridzin improve insulin sensitivity and enhance glucose uptake by subverting PPARγ/Cdk5 interaction in differentiated adipocytes.

Activators of peroxisome proliferator-activated receptor-γ (PPARγ agonists) are therapeutically promising candidates against insulin resistance and hyperglycemia. Synthetic PPARγ agonists are known to effectively enhance insulin sensitivity, but these are also associated with adverse side-effects and rising cost of treatment. Therefore, natural PPARγ targeting ligands are desirable alternatives for the management of insulin resistance associated with type 2 diabetes. Phloretin (PT) and Phloridzin (PZ) are predominant apple phenolics, which are recognized for their various pharmacological functions. The present study assessed the potential of PT and PZ in enhancing insulin sensitivity and glucose uptake by inhibiting Cdk5 activation and corresponding PPARγ phosphorylation in differentiated 3T3L1 cells. In silico docking and subsequent validation using 3T3L1 cells revealed that PT and PZ not only block the ser273 site of PPARγ but also inhibit the activation of Cdk5 itself, thereby, indicating their potent PPARγ regulatory attributes. Corroborating this, application of PT and PZ significantly enhanced the accumulation of cellular triglycerides as well as expression of insulin-sensitizing genes in adipocytes ultimately resulting in improved glucose uptake. Taken together, the present study reports that PT and PZ inhibit Cdk5 activation, which could be directly influencing the apparent PPARγ inhibition at ser273, ultimately resulting in improved insulin sensitivity and glucose uptake.

Styryl-cinnamate hybrid inhibits glioma by alleviating translation, bioenergetics and other key cellular responses leading to apoptosis.

Gliomas are lethal and aggressive form of brain tumors with resistance to conventional radiation and cytotoxic chemotherapies; inviting continuous efforts for drug discovery and drug delivery. Interestingly, small molecule hybrids are one such pharmacophore that continues to capture interest owing to their pluripotent medicinal effects. Accordingly, we earlier reported synthesis of potent Styryl-cinnamate hybrids (analogues of Salvianolic acid F) along with its plausible mode of action (MOA). We explored iTRAQ-LC/MS-MS technique to deduce differentially expressed landscape of native & phospho-proteins in treated glioma cells. Based on this, Protein-Protein Interactome (PPI) was looked into by employing computational tools and further validated in vitro. We hereby report that the Styryl-cinnamate hybrid, an analogue of natural Salvianolic acid F, alters key regulatory proteins involved in translation, cytoskeleton development, bioenergetics, DNA repair, angiogenesis and ubiquitination. Cell cycle analysis dictates arrest at G0/G1 stage along with reduced levels of cyclin D; involved in G1 progression. We discovered that Styryl-cinnamate hybrid targets glioma by intrinsically triggering metabolite-mediated stress. Various oncological circuits alleviated by the potential drug candidate strongly supports the role of such pharmacophores as anticancer drugs. Although, further analysis of SC hybrid in treating xenografts or solid tumors is yet to be explored but their candidature has gained huge impetus through this study. This study equips us better in understanding the shift in proteomic landscape after treating glioma cells with SC hybrid. It also allows us to elicit molecular targets of this potential drug before progressing to preclinical studies.

Gain of native conformation of Aurora A S155R mutant by small molecules.

Aurora A is a mitotic serine/threonine kinase protein that is a proposed target of the first-line anticancer drug design. It has been found to be overexpressed in many human cancer cells, including hematological, breast, and colorectal. Here, we focus on a particular somatic mutant S155R of Aurora kinase A protein, whose activity decreases because of loss of interaction with a TPX2 protein that results in ectopic expression of the Aurora kinase A protein, which contributes chromosome instability, centrosome amplification, and oncogenic transformation. The primary target of this study is to select a drug molecule whose binding results in gaining S155R mutant interaction with TPX2. The computational methodology applied in this study involves mapping of hotspots (for uncompetitive binding), virtual screening, protein-ligand docking, postdocking optimization, and protein-protein docking approach. In this study, we screen and validate ZINC968264, which acts as a potential molecule that can improve the loss of function occurred because of mutation (S155R) in Aurora A. Our approaches pave a suitable path to design a potential drug against physiological condition manifested because of S155R mutant in Aurora A.

Mutational Analysis on Membrane Associated Transporter Protein (MATP) and Their Structural Consequences in Oculocutaeous Albinism Type 4 (OCA4)-A Molecular Dynamics Approach.

Oculocutaneous albinism type IV (OCA4) is an autosomal recessive inherited disorder which is characterized by reduced biosynthesis of melanin pigmentation in skin, hair, and eyes and caused by the genetic mutations in the membrane-associated transporter protein (MATP) encoded by SLC45A2 gene. The MATP protein consists of 530 amino acids which contains 12 putative transmembrane domains and plays an important role in pigmentation and probably functions as a membrane transporter in melanosomes. We scrutinized the most OCA4 disease-associated mutation and their structural consequences on SLC45A2 gene. To understand the atomic arrangement in 3D space, the native and mutant structures were modeled. Further the structural behavior of native and mutant MATP protein was investigated by molecular dynamics simulation (MDS) approach in explicit lipid and water background. We found Y317C as the most deleterious and disease-associated SNP on SLC45A2 gene. In MDS, mutations in MATP protein showed loss of stability and became more flexible, which alter its structural conformation and function. This phenomenon has indicated a significant role in inducing OCA4. Our study explored the understanding of molecular mechanism of MATP protein upon mutation at atomic level and further helps in the field of pharmacogenomics to develop a personalized medicine for OCA4 disorder. J. Cell. Biochem. 117: 2608-2619, 2016. © 2016 Wiley Periodicals, Inc.

Impact of point mutation P29S in RAC1 on tumorigenesis.

A point mutation (P29S) in the RAS-related C3 botulinum toxin substrate 1 (RAC1) was considered to be a trigger for melanoma, a form of skin cancer with highest mortality rate. In this study, we have investigated the pathogenic role of P29S based on the conformational behavior of RAC1 protein toward guanosine triphosphate (GTP). Molecular interaction, molecular dynamics trajectory analysis (RMSD, RMSF, Rg, SASA, DSSP, and PCA), and shape analysis of binding pocket were performed to analyze the interaction energy and the dynamic behavior of native and mutant RAC1 at the atomic level. Due to this mutation, the RAC1 switch I region acquired more flexibility and, to compensate it, the switch II region becomes rigid in their conformational space, as a result of which the interaction energy of the protein for GTP increased. The overall results strongly implied that the changes in atomic conformation of the switch I and II regions in mutant RAC1 protein were a significant reason for its malignant transformation and tumorigenesis. We raised the opportunity for researchers to design possible therapeutic molecule by considering our findings.

Computational investigation of molecular mechanism and neuropathological implications in Huntington disease.

Huntington's disorder (HD), caused by mutations of the IT-15 gene, is an autosomal genetic disease that causes the breakdown of the nerve cells in the brain. The IT-15 gene encodes the huntingtin (Htt) protein. Htt, along with its interacting partners, are involved in maintaining proper communication among neurons. Our work is based on the interaction behavior between Htt (in three polyglutamine (polyQ) states that is Htt 0Q, 17Q and 36Q) and SH3GL3 interacting protein by using computational methods. We used the HADDOCK docking platform to find out the extent of interaction between Htt polyQ models and SH3GL3. The Htt36Q (mutated) showed higher interaction than Htt17Q (native) with SH3GL3. Molecular dynamics simulation was performed to uncover the structural fluctuations of polyQ models and their complexes. RMSD, Rg, SASA, and total interaction energy graph showed significant results, where as mutant Htt showed higher fluctuations and flexibility than native Htt. The increase in the length of polyQ was found to affect the stability, flexibility, and compactness of the protein and its complex. Our research provided a propitious approach to understand the consequence of polyglutamination in Htt and its relation with HD.