Comparative molecular field analysis and molecular dynamics studies of α/β hydrolase domain containing 6 (ABHD6) inhibitors

Structure-based virtual screening methods for the identification of novel phytochemical inhibitors targeting furin protease for the management of COVID-19

The coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus, is a highly contagious respiratory disease with widespread societal impact. The symptoms range from cough, fever, and pneumonia to complications affecting various organs, including the heart, kidneys, and nervous system. Despite various ongoing efforts, no effective drug has been developed to stop the spread of the virus. Although various types of medications used to treat bacterial and viral diseases have previously been employed to treat COVID-19 patients, their side effects have also been observed. The way SARS-CoV-2 infects the human body is very specific, as its spike protein plays an important role. The S subunit of virus spike protein cleaved by human proteases, such as furin protein, is an initial and important step for its internalization into a human host. Keeping this context, we attempted to inhibit the furin using phytochemicals that could produce minimal side effects. For this, we screened 408 natural phytochemicals from various plants having antiviral properties, against furin protein, and molecular docking and dynamics simulations were performed. Based on the binding score, the top three compounds (robustaflavone, withanolide, and amentoflavone) were selected for further validation. MM/GBSA energy calculations revealed that withanolide has the lowest binding energy of -57.2 kcal/mol followed by robustaflavone and amentoflavone with a binding energy of -45.2 kcal/mol and -39.68 kcal/mol, respectively. Additionally, ADME analysis showed drug-like properties for these three lead compounds. Hence, these natural compounds robustaflavone, withanolide, and amentoflavone, may have therapeutic potential for the management of SARS-CoV-2 by targeting furin.

Identification of small molecule inhibitors of RAD52 for breast cancer therapy: in silico approach

The breast cancer susceptibility gene 1/2 (BRCA1/2) are the key regulators in maintaining the genomic integrity and mutations in these genes have been associated with development of breast and ovarian cancers. Also, synthetic lethality has been shown in BRCA1/2 deficient cancers, when the RAD52 gene is silenced by shRNA or small molecules aptamers, suggesting a role for RAD52 in the breast cancers pathogenesis. Thus, to find the potential inhibitors of RAD52, a collection of 21,000 compounds from the ChemBridge screening library was screened to conduct molecular docking and molecular dynamics simulation (MD) against RAD52. Further, the results were validated by a density functional theory (DFT) analysis and using post-dynamics free energy calculations. Out of all screened molecules, the docking study revealed five compounds were found to have promising activities against RAD52. Moreover, the catalytic amino acid residues of RAD52 developed stable contacts with compound 8758 and 10593, as anticipated by DFT calculation, MD simulation, and post dynamics MM-GBSA energy calculation. It appears that compound 8758 is the best inhibitor against RAD52 followed by 10593 compared to the other top hits, in terms of the HOMO orbital energy (-1.0966 eV and -1.2136 eV) from DFT and the post dynamics binding free energy calculation (-54.71 and -52.43 Kcal/mol). Furthermore, a drug-like properties of lead molecules (8758 and 10593) were also seen via ADMET analysis. Based on our computational analysis, we hypothesize that a small molecule 8758 and 10593 possess the therapeutic potential in the management for breast cancer patients with a BRCA mutation via targeting RAD52.Communicated by Ramaswamy H. Sarma.

Pharmacoinformatics and Breed-Based De Novo Hybridization Studies to Develop New Neuraminidase Inhibitors as Potential Anti-Influenza Agents

Influenza represents a profoundly transmissible viral ailment primarily afflicting the respiratory system. Neuraminidase inhibitors constitute a class of antiviral therapeutics employed in the management of influenza. These inhibitors impede the liberation of the viral neuraminidase protein, thereby impeding viral dissemination from the infected cell to host cells. As such, neuraminidase has emerged as a pivotal target for mitigating influenza and its associated complications. Here, we apply a de novo hybridization approach based on a breed-centric methodology to elucidate novel neuraminidase inhibitors. The breed technique amalgamates established ligand frameworks with the shared target, neuraminidase, resulting in innovative inhibitor constructs. Molecular docking analysis revealed that the seven synthesized breed molecules (designated Breeds 1-7) formed more robust complexes with the neuraminidase receptor than conventional clinical neuraminidase inhibitors such as zanamivir, oseltamivir, and peramivir. Pharmacokinetic evaluations of the seven breed molecules (Breeds 1-7) demonstrated favorable bioavailability and optimal permeability, all falling within the specified parameters for human application. Molecular dynamics simulations spanning 100 nanoseconds corroborated the stability of these breed molecules within the active site of neuraminidase, shedding light on their structural dynamics. Binding energy assessments, which were conducted through MM-PBSA analysis, substantiated the enduring complexes formed by the seven types of molecules and the neuraminidase receptor. Last, the investigation employed a reaction-based enumeration technique to ascertain the synthetic pathways for the synthesis of the seven breed molecules.

Virtual screening of gut microbiome bacteriocins as potential inhibitors of stearoyl-CoA desaturase 1 to regulate adipocyte differentiation and thermogenesis to combat obesity

The gut bacterial strains and their metabolites have been shown to play a significant role in obesity, but the molecular mechanisms underlying this association are largely unresolved. Obesity is a multifactorial problem and is controlled by various mechanisms and pathways to produce and store fat cells. Bacteriocins are secondary metabolites produced by gut bacteria to defend themselves against their competitors. Recently, they have gained great attention due to their role in metabolic disorders, including obesity. Stearoyl-CoA desaturase 1 (SCD1) is a key enzyme involved in the differentiation of adipocytes. The aim of this study is to show the regulation of SCD1 by bacteriocins and thus their importance in obesity control. We screened the human gut bacteriome for the presence of bacteriocins, predicted their structures, and showed their inhibitory role by molecular docking with SCD1. Further, to confirm the docking results, MDS of six top scoring SCD1-bacteriocin complexes were carried out for 100 ns. These six bacteriocins namely, Plantaricin S-beta, Carnolysin, Lactococcin B, Bacteriocin Iic, Plantaricin N, and Thermophilin A, with strong binding affinities, are primarily produced by bacterial strains from the Lactobacillaeacea family. These findings can be the basis of further experiments for enhanced understanding of the underlying mechanisms for obesity control, specifically bacteriocins driven regulation of the SCD1 enzyme. In addition, a consortium of bacterial strains producing these bacteriocins can be developed and used as probiotics for the amelioration of obesity and other metabolic complications.Communicated by Ramaswamy H. Sarma.

Design, Synthesis, and Biological Evaluation of Notopterol Derivatives as Triple Inhibitors of AChE/BACE1/GSK3β for the Treatment of Alzheimer's Disease

The pathogenesis of Alzheimer's disease (AD) is very complex, and there are many hypotheses. Therefore, the development of a multi-target-directed-ligand may be an effective therapeutic strategy. Our previous study showed that notopterol (a natural product from Notopterygium) is a dual BACE1/GSK3β inhibitor. In this study, we designed and synthesized 48 notopterol derivatives with furacoumarin as a scaffold in order to enhance their balanced AChE/BACE1/GSK3β inhibitory activity. Fortunately, 1c showed effective inhibitory activity against AChE (58.7% at 1.0 μM), BACE1 (48.3% at 20 μM), and GSK3β (40.3% at 10 μM). Furthermore, 1c showed good blood-brain barrier penetrability, suitable bioavailability, and oral safety. More importantly, 1c could ameliorate the impaired learning and memory in Aβ-induced AD mice. In conclusion, we reported the triple inhibitor of AChE/BACE1/GSK3β lead compounds based on a furocoumarin scaffold of notopterol for the first time, which provides a potential new strategy for the treatment of AD.

Identification of a promising inhibitor from Illicium verum (star anise) against the main protease of SARS-CoV-2: insights from the computational study

SARS-CoV-2, the causing agent of coronavirus disease (COVID-19), first broke out in Wuhan and rapidly spread worldwide, resulting in a global health emergency. The lack of specific drugs against the coronavirus has made its spread challenging to control. The main protease (M^pro) is a key enzyme of SARS-CoV-2 used as a key target in drug discovery against the coronavirus. Medicines derived from plant phytoconstituents have been widely exploited to treat various diseases. The present study has evaluated the potential of Illicium verum (star anise) phytoconstituents against M^pro by implementing a computational approach. We performed molecular docking and molecular dynamics simulation study with a set of 60 compounds to identify their potential to inhibit the main protease (M^pro) of SARS-CoV-2. DFT study and post dynamics free energy calculations were also performed to strengthen the findings. The identified four compounds by docking study exhibited the highest potential compared to other selected phytoconstituents. Further, density functional theory (DFT) calculation, molecular dynamics simulation and post dynamics MM-GBSA energy calculation predicted Verimol-G as a potential compound, which formed stable interactions through the catalytic dyad residues. The HOMO orbital energy (-0.250038) from DFT and the post dynamics binding free energy calculation (-73.33 Kcal/mol) correlate, suggesting Verimol-G is the best inhibitor compared to the other phytoconstituents. This compound also complies with the ADME properties of drug likeliness. Thus, based on a computational study, we suggest that Verimol G may be developed as a potential inhibitor against the main protease to combat COVID-19.Communicated by Ramaswamy H. Sarma.

ABHD6 Controls Amphetamine-Stimulated Hyperlocomotion: Involvement of CB1 Receptors

Introduction: Activation of cannabinoid 1 receptors (CB₁Rs) by endocannabinoids (eCBs) is controlled by both eCB production and eCB inactivation. Accordingly, inhibition of eCB hydrolyzing enzymes, monoacylglycerol lipase (MAGL) and α/β-hydrolase domain containing 6 (ABHD6), enhances eCB accumulation and CB₁R activation. It is known that inhibition of MAGL regulates select CB₁R-dependent behaviors in mice, including locomotor behaviors and their modulation by psychostimulants, but much less is known about the effect of inhibiting ABHD6 activity on such behaviors. Methods: We report a new mouse line that carries a genetic deletion of Abhd6 and evaluated its effect on spontaneous locomotion measured in a home cage monitoring system, motor coordination measured on a Rotarod, and amphetamine-stimulated hyperlocomotion and amphetamine sensitization (AS) measured in an open-field chamber. Results: ABHD6 knockout (KO) mice reached adulthood without exhibiting overt behavioral impairment, and we measured only mild reduction in spontaneous locomotion and motor coordination in adult ABHD6 KO mice compared to wild-type (WT) mice. Significantly, amphetamine-stimulated hyperlocomotion was enhanced by twofold in ABHD6 KO mice compared to WT mice and yet ABHD6 KO mice expressed AS to the same extent as WT mice. A twofold increase in amphetamine-stimulated hyperlocomotion was also measured in ABHD6 heterozygote mice and in WT mice treated with the ABHD6 inhibitor KT-182. It is known that amphetamine-stimulated hyperlocomotion is not affected by the CB₁R antagonist, SR141617, and we discovered that the enhanced amphetamine-stimulated hyperlocomotion resulting from ABHD6 inhibition is blocked by SR141617. Conclusions: Our study suggests that ABHD6 controls amphetamine-stimulated hyperlocomotion by a mechanistic switch to a CB₁R-dependent mechanism.

In silico identification of promising inhibitor against RNA-dependent RNA polymerase target of SARS-CoV-2

The severe acute respiratory syndrome is a viral respiratory disease recognised as COVID-19, caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Formerly, no precise remedies are available, and many studies regarding COVID-19 prevention and treatment are under development. Several targets for the design of drugs are identified, and studies are in headway to explore the potential target. RNA-dependent RNA polymerase (RdRp) protein identified as a promising target against SARS-CoV-2 infection for the drug design due to its significant role in viral replication. The present study focuses on identifying the binding effect of previously known RdRp inhibitors with RdRp of SARS-CoV-2 using molecular docking and molecular dynamics simulation approaches. Molecular docking and binding free energy calculations against RdRp enzyme identified suramin as a potential compound that showed the highest docking score of -7.83 Kcal/mole and binding energy of -80.83 Kcal/mole as a comparison to other compounds. Further, molecular dynamics simulation studies were moreover showed the stable binding behaviour of suramin docked complex in the protein active site. Thus, the study concludes that suramin might be helpful as a potential inhibitor against RNA-dependent RNA polymerase of SRAS-CoV-2. However, further investigation is needed to assess the possible effect of inhibitors on RdRp through in vitro and in vivo experiments.

ABHD6 Inhibition Rescues a Sex-Dependent Deficit in Motor Coordination in The HdhQ200/200 Mouse Model of Huntington's Disease

Huntington's Disease is associated with motor behavior deficits that are lessened by few therapeutic options. This preliminary study tested if pharmacological inhibition of α/β-hydrolase domain containing 6 (ABHD6), a multifunctional enzyme expressed in the striatum, rescues behavioral deficits in HdhQ200/200 mice. Previous work has shown that this model exhibits a reduction in spontaneous locomotion and motor coordination at 8 and 10 months of age, with a more severe phenotype in female mice. Semi-quantitative immunohistochemistry analysis indicated no change in striatal ABHD6 expression at 8 months of age, but a 40% reduction by 10 months in female HdhQ200/200 mice compared to female wild-type (WT) littermates. At 8 months of age, acute ABHD6 inhibition rescued motor coordination deficits in female HdhQ200/200 mice without affecting WT performance. ABHD6 inhibition did not impact spontaneous locomotion, grip strength, or overall weight in either group, showing that effects were specific to motor coordination. At 10 months of age, semi-chronic ABHD6 inhibition by osmotic pump delivery also rescued motor coordination deficits in female HdhQ200/200 mice without affecting female WT littermates. Our preliminary study suggests that ABHD6 inhibition improves motor performance in female HdhQ200/200 mice.

α/β-Hydrolase Domain (ABHD) Inhibitors as New Potential Therapeutic Options against Lipid-Related Diseases

Much of the experimental evidence in the literature has linked altered lipid metabolism to severe diseases such as cancer, obesity, cardiovascular pathologies, diabetes, and neurodegenerative diseases. Therefore, targeting key effectors of the dysregulated lipid metabolism may represent an effective strategy to counteract these pathological conditions. In this context, α/β-hydrolase domain (ABHD) enzymes represent an important and diversified family of proteins, which are involved in the complex environment of lipid signaling, metabolism, and regulation. Moreover, some members of the ABHD family play an important role in the endocannabinoid system, being designated to terminate the signaling of the key endocannabinoid regulator 2-arachidonoylglycerol. This Perspective summarizes the research progress in the development of ABHD inhibitors and modulators: design strategies, structure-activity relationships, action mechanisms, and biological studies of the main ABHD ligands will be highlighted.

Docking-Based 3D-QSAR Studies for 1,3,4-oxadiazol-2-one Derivatives as FAAH Inhibitors

This work aimed to construct 3D-QSAR CoMFA and CoMSIA models for a series of 31 FAAH inhibitors, containing the 1,3,4-oxadiazol-2-one moiety. The obtained models were characterized by good statistical parameters: CoMFA Q2 = 0.61, R2 = 0.98; CoMSIA Q2 = 0.64, R2 = 0.93. The CoMFA model field contributions were 54.1% and 45.9% for steric and electrostatic fields, respectively. In the CoMSIA model, electrostatic, steric, hydrogen bond donor, and hydrogen acceptor properties were equal to 34.6%, 23.9%, 23.4%, and 18.0%, respectively. These models were validated by applying the leave-one-out technique, the seven-element test set (CoMFA r2test-set = 0.91; CoMSIA r2test-set = 0.91), a progressive scrambling test, and external validation criteria developed by Golbraikh and Tropsha (CoMFA r20 = 0.98, k = 0.95; CoMSIA r20 = 0.98, k = 0.89). As the statistical significance of the obtained model was confirmed, the results of the CoMFA and CoMSIA field calculation were mapped onto the enzyme binding site. It gave us the opportunity to discuss the structure-activity relationship based on the ligand-enzyme interactions. In particular, examination of the electrostatic properties of the established CoMFA model revealed fields that correspond to the regions where electropositive substituents are not desired, e.g., in the neighborhood of the 1,3,4-oxadiazol-2-one moiety. This highlights the importance of heterocycle, a highly electronegative moiety in this area of each ligand. Examination of hydrogen bond donor and acceptor properties contour maps revealed several spots where the implementation of another hydrogen-bond-donating moiety will positively impact molecules' binding affinity, e.g., in the neighborhood of the 1,3,4-oxadiazol-2-one ring. On the other hand, there is a large isopleth that refers to the favorable H-bond properties close to the terminal phenoxy group of a ligand, which means that, generally speaking, H-bond acceptors are desired in this area.

Binding of SARS-CoV Covalent Non-Covalent Inhibitors to the SARS-CoV-2 Papain-Like Protease and Ovarian Tumor Domain Deubiquitinases

The urgent need for novel and effective drugs against the SARS-CoV-2 coronavirus pandemic has stimulated research worldwide. The Papain-like protease (PLpro), which is essential for viral replication, shares a similar active site structural architecture to other cysteine proteases. Here, we have used representatives of the Ovarian Tumor Domain deubiquitinase family OTUB1 and OTUB2 along with the PLpro of SARS-CoV-2 to validate and rationalize the binding of inhibitors from previous SARS-CoV candidate compounds. By forming a new chemical bond with the cysteine residue of the catalytic triad, covalent inhibitors irreversibly suppress the protein’s activity. Modeling covalent inhibitor binding requires detailed knowledge about the compounds’ reactivities and binding. Molecular Dynamics refinement simulations of top poses reveal detailed ligand-protein interactions and show their stability over time. The recently discovered selective OTUB2 covalent inhibitors were used to establish and validate the computational protocol. Structural parameters and ligand dynamics are in excellent agreement with the ligand-bound OTUB2 crystal structures. For SARS-CoV-2 PLpro, recent covalent peptidomimetic inhibitors were simulated and reveal that the ligand-protein interaction is very dynamic. The covalent and non-covalent docking plus subsequent MD refinement of known SARS-CoV inhibitors into DUBs and the SARS-CoV-2 PLpro point out a possible approach to target the PLpro cysteine protease from SARS-CoV-2. The results show that such an approach gives insight into ligand-protein interactions, their dynamic character, and indicates a path for selective ligand design.

Identification of bioactive molecule from Withania somnifera (Ashwagandha) as SARS-CoV-2 main protease inhibitor

SARS-CoV-2 is the causative agent of COVID-19 and has been declared as pandemic disease by World Health Organization. Lack of targeted therapeutics and vaccines for COVID-2019 have triggered the scientific community to develop new vaccines or drugs against this novel virus. Many synthetic compounds and antimalarial drugs are undergoing clinical trials. The traditional medical practitioners widely use Indian medicinal plant Withania somnifera (Ashwagandha) natural constituents, called withanolides for curing various diseases. The main protease (M^pro) of SARS-CoV-2 plays a vital role in disease propagation by processing the polyproteins which are required for its replication. Hence, it denotes a significant target for drug discovery. In the present study, we evaluate the potential of 40 natural chemical constituents of Ashwagandha to explore a possible inhibitor against main protease of SARS-CoV-2 by adopting the computational approach. The docking study revealed that four constituents of Ashwagandha; Withanoside II (-11.30 Kcal/mol), Withanoside IV (-11.02 Kcal/mol), Withanoside V (-8.96 Kcal/mol) and Sitoindoside IX (-8.37 Kcal/mol) exhibited the highest docking energy among the selected natural constituents. Further, MD simulation study of 100 ns predicts Withanoside V possess strong binding affinity and hydrogen-bonding interactions with the protein active site and indicates its stability in the active site. The binding free energy score also correlates with the highest score of -87.01 ± 5.01 Kcal/mol as compared to other selected compounds. In conclusion, our study suggests that Withanoside V in Ashwagandha may be serve as a potential inhibitor against M^pro of SARS-CoV-2 to combat COVID-19 and may have an antiviral effect on nCoV.Communicated by Ramaswamy H. Sarma.

The Universal 3D QSAR Model for Dopamine D2 Receptor Antagonists

In order to search for novel antipsychotics acting through the D₂ receptor, it is necessary to know the structure–activity relationships for dopamine D₂ receptor antagonists. In this context, we constructed the universal three-dimensional quantitative structure–activity relationship (3D- QSAR) model for competitive dopamine D₂ receptor antagonists. We took 176 compounds from chemically different groups characterized by the half maximal inhibitory concentration (IC₅₀)from the CHEMBL database and docked them to the X-ray structure of the human D₂ receptor in the inactive state. Selected docking poses were applied for Comparative Molecular Field Analysis (CoMFA) alignment. The obtained CoMFA model is characterized by a cross-validated coefficient Q² of 0.76 with an optimal component of 5, R² of 0.92, and an F value of 338.9. The steric and electrostatic field contributions are 67.4% and 32.6%, respectively. The statistics obtained prove that the CoMFA model is significant. Next, the IC₅₀ of the 16 compounds from the test set was predicted with R² of 0.95. Finally, a progressive scrambling test was carried out for additional validation. The CoMFA fields were mapped onto the dopamine D₂ receptor binding site, which enabled a discussion of the structure–activity relationship based on ligand–receptor interactions. In particular, it was found that one of the desired steric interactions covers the area of a putative common allosteric pocket suggested for some other G protein-coupled receptors (GPCRs), which would suggest that some of the known dopamine receptor antagonists are bitopic in their essence. The CoMFA model can be applied to predict the potential activity of novel dopamine D₂ receptor antagonists.

Computational Approaches for Elucidating Protein-Protein Interactions in Cation Channel Signaling

BACKGROUND

The lipid bilayer of the plasma membrane is impermeable to ions, yet changes in the flux of ions across the cell membrane are critical regulatory events in cells. Because of their regulatory roles in a range of physiological processes, such as electrical signaling in muscles and neurons, to name a few, these proteins are one of the most important drug targets.

OBJECTIVE

This review mainly focused on the computational approaches for elucidating proteinprotein interactions in cation channel signaling.

DISCUSSION

Due to continuously advanced facilities and technologies in computer sciences, the physical contacts of macromolecules of channel structures have been virtually visualized. Indeed, techniques like protein-protein docking, homology modeling, and molecular dynamics simulation are valuable tools for predicting the protein complex and refining channels with unreleased structures. Undoubtedly, these approaches will greatly expand the cation channel signaling research, thereby speeding up structure-based drug design and discovery.

CONCLUSION

We introduced a series of valuable computational tools for elucidating protein-protein interactions in cation channel signaling, including molecular graphics, protein-protein docking, homology modeling, and molecular dynamics simulation.

ABHD6: Its Place in Endocannabinoid Signaling and Beyond

The endocannabinoid (eCB) signaling system modulates neurotransmission and inflammation, among other physiological functions. Its newest member, α/β-hydrolase domain-containing 6 (ABHD6), has emerged as a promising therapeutic target to treat several devastating diseases, including epilepsy. Here, we review the molecular mechanisms that mediate and control eCB signaling and, within it, the specific role of ABHD6. We also discuss how ABHD6 controls the abundance of additional lipids and the trafficking of ionotropic receptors to plasma membranes. We finish with several unexplored questions regarding this novel enzyme. Our current understanding of the molecular mechanism and biological function of ABHD6 provides a strong foundation for the development of small-molecule therapeutics to treat devastating diseases.

Structural properties and role of the endocannabinoid lipases ABHD6 and ABHD12 in lipid signalling and disease

The endocannabinoid (eCB) system is an important part of both the human central nervous system (CNS) and peripheral tissues. It is involved in the regulation of various physiological and neuronal processes and has been associated with various diseases. The eCB system is a complex network composed of receptor molecules, their cannabinoid ligands, and enzymes regulating the synthesis, release, uptake, and degradation of the signalling molecules. Although the eCB system and the molecular processes of eCB signalling have been studied extensively over the past decades, the involved molecules and underlying signalling mechanisms have not been described in full detail. An example pose the two poorly characterised eCB-degrading enzymes α/β-hydrolase domain protein six (ABHD6) and ABHD12, which have been shown to hydrolyse 2-arachidonoyl glycerol-the main eCB in the CNS. We review the current knowledge about the eCB system and the role of ABHD6 and ABHD12 within this important signalling system and associated diseases. Homology modelling and multiple sequence alignments highlight the structural features of the studied enzymes and their similarities, as well as the structural basis of disease-related ABHD12 mutations. However, homologies within the ABHD family are very low, and even the closest homologues have widely varying substrate preferences. Detailed experimental analyses at the molecular level will be necessary to understand these important enzymes in full detail.

Identification of isobavachalcone as a potential drug for rice blast disease caused by the fungus Magnaporthe grisea

The rice blast disease caused by the fungus Magnaporthe grisea is one of the most devastating rice diseases, but there is no effective fungicide toward chitinase which is a key enzyme of M. grisea. In this study, we observed that distortion and cell-wall damage of M. grisea hyphae were significantly under the scanning electron micrograph after a 24-h treatment with 10 mg/L isobavachalcone (IBC) extracted from Psoralea corylifolia L. To further explore the effect of IBC on the cell wall of M. grisea, we examined changes in enzymes associated with cell wall degradation by enzyme activity experiments, treated liquid culture mycelia with 10 mg/L IBC for 1 h. Results displayed that chitinase was obviously more active than control group. To illustrate the interactions between IBC and chitinase, the studies of homology modeling and molecular docking were carried out successively. The results revealed that IBC had hydrogen bonds with residues ASP267 and ARG276 of chitinase. Besides, these nonpolar residues TYR270, PRO271, VAL272, LEU310, PRO311, TYR316, and LEU317 were able to form strong hydrophobic interactions. Binding energies of the chitinase-IBC complexes were calculated by MM-GBSA showed that the ΔG_bind score of molecular dynamics had lower binding energy and more stable than docking complexes. All above, IBC owns significant agonistic activity in chitinase and would be a potent fungicide to inhibit the growth of M. grisea. We hope the above information provides an important insight for understanding the interactions between IBC and chitinase, which may be useful in the discovery of a novel potent agonist. Communicated by Ramaswamy H. Sarma.

In silico exploration of aryl sulfonamide analogs as voltage-gated sodium channel 1.7 inhibitors by using 3D-QSAR, molecular docking study, and molecular dynamics simulations

It has been demonstrated by human genetics that the voltage-gated sodium channel Nav1.7 is currently a promising target for the treatment of pain. In this research, we performed molecular simulation works on a series of classic aryl sulfonamide Nav1.7 inhibitors using three-dimensional quantitative structure-activity relationships (3D-QSAR), molecular docking and molecular dynamics (MD) simulations for the first time to explore the correlation between their structures and activities. The results of the relevant statistical parameters of comparative molecular field analyses (CoMFA) and comparative molecular similarity indices analyses (CoMSIA) had been verified to be reasonable, and the deep relationship between the structures and activities of these inhibitors was obtained by analyzing the contour maps. The generated 3D-QSAR model showed a good predictive ability and provided valuable clues for the rational modification of molecules. The interactions between compounds and proteins were modeled by molecular docking studies. Finally, accuracy of the docking results and stability of the complexes were verified by 100 ns MD simulations. Detailed information on the key residues at the binding site and the types of interactions they participate in involved was obtained. The van der Waals energy contributed the most in the molecular binding process according to the calculation of binding free energy. All research results provided a good basis for further research on novel and effective Nav1.7 inhibitors.

Do new N-substituted 3-amino-4-phenyl-5-oxo-pyrazolinecarboxamide derivatives exhibit antitubercular potential?

As a continuation of previous tests concerning new N-substituted 3-amino-4-phenyl-5-oxo-pyrazolinecarboxamide derivatives (R3, R4 and R8) of notable antibacterial activity, their antitubercular potential against different mycobacterial strains was estimated. Tests performed on virulent (reference and clinical) strains of Mycobacterium bovis and Mycobacterium tuberculosis revealed the highest therapeutic potential of R8 derivative: MIC within the range 7.8-15.6 μg/ml and TI (therapeutic index) within the range 46.5-93. Moreover, the synergistic interaction was found between R3, R4 and R8 derivatives and rifampicin, one of the front-line antitubercular drugs. R8/rifampicin mixture in concentrations effective in inhibition of Mycobacterium tuberculosis strain was non-cytotoxic against GMK cells, displaying cell viability approximately 88-97% when compared to control. Molecular docking study enabled to conclude that enoyl acyl carrier protein reductase (InhA) can be considered as a potential molecular target of tested pyrazole derivatives. Although further modifications of chemical structure of the investigated pyrazole derivatives is required, in order to increase their antitubercular efficacy and therapeutic safety, these compounds, in particular R8 compound, can be promising for the treatment of human and bovine tuberculosis.

Computational investigation of TGF-β receptor inhibitors for treatment of idiopathic pulmonary fibrosis: Field-based QSAR model and molecular dynamics simulation

The discovery of drugs relevant to transforming growth factor β (TGF-β) receptor inhibitors have been considered as a considerable challenge during therapy idiopathic pulmonary fibrosis diseases. For the first time, herein we illustrate a field-based quantitative structure-activity relationship (QSAR) model and molecular dynamics (MD) simulations for novel 7-substituted-pyrazolo [4, 3-b] pyridine derivatives with biological activity for the TGF-β receptor, with an attempt of elucidating the 3D structural features that are essential for the activity. Results demonstrate that the field-based model (Q² = 0.548, R²_training = 0.840, R²_test = 0.750) are acceptable with good predictive capabilities. In addition, MD studies were also carried out on the training set with the aim of exploring their binding modes in the active pocket of TGF-β receptor, resulting in some of the crucial structural fragments which are responsible for inhibitory activity. Therefore, we summarized the following features required for TGF-β receptor inhibition: electronegative in region1, bulky groups in region2 and smaller groups in region3. Based on the model and related information, we hope the above information provides an important insight for understanding the interactions of the inhibitors and TGF-β receptor, which may be useful in discovering novel potent inhibitors.

Molecular docking, 3D-QSAR and structural optimization on imidazo-pyridine derivatives dually targeting AT1 and PPARg

Telmisartan, a bifunctional agent of blood pressure lowering and glycemia reduction, was previously reported to antagonize angiotensin II type 1 (AT1) receptor and partially activate peroxisome proliferator-activated receptor γ (PPARγ) simultaneously. Through the modification to telmisartan, researchers designed and obtained imidazo-\pyridine derivatives with the IC₅₀s of 0.49∼94.1 nM against AT1 and EC₅₀s of 20∼3640 nM towards PPARγ partial activation. For minutely inquiring the interaction modes with the relevant receptor and analyzing the structure-activity relationships, molecular docking and 3D-QSAR (Quantitative structure-activity relationships) analysis of these imidazo-\pyridines on dual targets were conducted in this work. Docking approaches of these derivatives with both receptors provided explicit interaction behaviors and excellent matching degree with the binding pockets. The best CoMFA (Comparative Molecular Field Analysis) models exhibited predictive results of q²=0.553, r²=0.954, SEE=0.127, r²_pred=0.779 for AT1 and q²=0.503, r²=1.00, SEE=0.019, r²_pred=0.604 for PPARγ, respectively. The contour maps from the optimal model showed detailed information of structural features (steric and electrostatic fields) towards the biological activity. Combining the bioisosterism with the valuable information from above studies, we designed six molecules with better predicted activities towards AT1 and PPARγ partial activation. Overall, these results could be useful for designing potential dual AT1 antagonists and partial PPARγ agonists.

Comparative molecular field analysis and molecular dynamics studies of the dopamine D2 receptor antagonists without a protonatable nitrogen atom

The dopaminergic hypothesis of schizophrenia is the main concept explaining the direct reasons of schizophrenia and the effectiveness of current antipsychotics. All antipsychotics present on the market are potent dopamine D₂ receptor antagonists or partial agonists. In this work we investigate a series of dopamine D₂ receptor antagonists which do not fulfill the criteria of the classical pharmacophore model as they do not possess a protonatable nitrogen atom necessary to interact with the conserved Asp(3.32). Such compounds are interesting, inter alia, due to possible better pharmacokinetic profile when compared to basic, ionizable molecules. By means of homology modeling, molecular docking and molecular dynamics we determined that the compounds investigated interact with Asp(3.32) via their amide nitrogen atom. It was found that the studied compounds stabilize the receptor inactive conformation through the effect on the ionic lock, which is typical for GPCR antagonists. We constructed a CoMFA model for the studied compounds with the following statistics: R² = 0.95, Q² = 0.63. The quality of the CoMFA model was confirmed by high value of R² of the test set, equal 0.96. The CoMFA model indicated two regions where bulky substituents are favored and two regions where bulky substituents are not beneficial. Two red contour regions near carbonyl groups were identified meaning that negative charge would be favored here. Furthermore, the S-oxide group is connected with blue contour region meaning that positive charge is favored in this position. These findings may be applied for further optimization of the studied compound series.

On the virtues of automated quantitative structure-activity relationship: the new kid on the block

Quantitative structure-activity relationship (QSAR) has proved to be an invaluable tool in medicinal chemistry. Data availability at unprecedented levels through various databases have collaborated to a resurgence in the interest for QSAR. In this context, rapid generation of quality predictive models is highly desirable for hit identification and lead optimization. We showcase the application of an automated QSAR approach, which randomly selects multiple training/test sets and utilizes machine-learning algorithms to generate predictive models. Results demonstrate that AutoQSAR produces models of improved or similar quality to those generated by practitioners in the field but in just a fraction of the time. Despite the potential of the concept to the benefit of the community, the AutoQSAR opportunity has been largely undervalued.

Design, sythesis and evaluation of a series of 3- or 4-alkoxy substituted phenoxy derivatives as PPARs agonists

Peroxisome proliferators-activated receptors (PPARα, γ and δ) are potentially effective targets for Type 2 diabetes mellitus therapy. The severe effects of known glitazones and the successfully approved agents (saroglitazar and lobeglitazone) motivated us to study novelly potent PPARs drugs with improved safety profile. In this work, we received 15 carboxylic acids based on the combination principle to integrate the polar head of bezafibrate with the hydrophobic tail of pioglitazone. Another 12 tetrazoles based on the bioisosterism principle were obtained accordingly. Furthermore, in vitro PPARs transactivation assays on these 3- or 4-alkoxy substituted phenoxy derivatives afforded six compounds. Interactions and binding stability from the docking analysis and 20 ns molecular dynamic simulations confirmed the representative compounds to be suitable and plausible for PPARs pockets. The above-mentioned results demonstrated that the compounds may be used as reference for further optimization for enhanced PPARs activities and wide safety range.

Identification of dual ligands targeting angiotensin II type 1 receptor and peroxisome proliferator-activated receptor-γ by core hopping of telmisartan

It has been reported previously that some angiotensin II receptor blockers not only antagonize angiotensin II type 1 receptor (AT₁R), but also exert stimulation in peroxisome proliferator-activated receptor γ (PPARγ) partial activation, among which telmisartan displays the best. Telmisartan has been tested as a bifunctional ligand with antihypertensive and hypoglycemic activity. Aiming at more potent leads with selective AT₁R antagonism and PPARγ partial agonism, the three parts of telmisartan including the distal benzimidazole ring, the biphenyl moiety, and the carboxylic acid group experienced modification by core hopping method in our study. The central benzimidazole ring, however, remained intact considering its great affinity toward AT₁R and PPARγ. We utilized computational techniques for the sake of details on the binding interactions and conformational stability. Standard precision docking analysis and absorption, distribution, metabolism, excretion, and toxicity prediction received 10 molecules with higher Glide scores, similar interactions, and improved pharmacokinetic profiles compared to telmisartan. Comp#91 with highest scores for AT₁R (-11.92 kcal/mol) and PPARγ (-13.88 kcal/mol) exhibited excellent binding modes and pharmacokinetic parameters. Molecular dynamics trajectories on best docking pose of comp#91 confirmed the docking results and verified the conformational stability with both receptors throughout the course of 20-ns simulations. Thus, comp#91 could be identified as a promising lead in the development of dual AT₁R antagonist and PPARγ partial agonist against hypertension and type 2 diabetes.

Patent Highlight

A snapshot of noteworthy recent developments in the patent literature of relevance to pharmaceutical and medical research and development.