Critical assessment of the automated AutoDock as a new docking tool for virtual screening

Unusual Ni⋯Ni interaction in Ni(ii) complexes as potential inhibitors for the development of new anti-SARS-CoV-2 Omicron drugs

Two nickel(ii) coordination complexes [Ni(L)]2(1) and [Ni(L)]n(2) of a tetradentate Schiff base ligand (H2L) derived from 2-hydroxy-1-naphthaldehyde with ethylenediamine were synthesized, designed, and characterized via spectroscopic and single crystal XRD analyses. Both nickel(ii) complexes exhibited unusual Ni⋯Ni interactions and were fully characterized via single-crystal X-ray crystallography. Nickel(ii) complexes [Ni(L)]2(1) and [Ni(L)]n(2) crystallize in monoclinic and triclinic crystal systems with P21/c and P1̄ space groups, respectively, and revealed square planar geometry around each Ni(ii) ion. The structure of both the complexes have established the existence of a new kind of metal system containing nickel(ii)-nickel(ii) interactions with a square planar-like geometry about the nickel(ii) atoms. Both square planar Ni(ii) complexes were often stacked with relatively short Ni⋯Ni distances. The non-bonded Ni-Ni distance (Ni⋯Ni separation) seems to be 3.356 Å and 3.214 Å from the nickel atoms of [Ni(L)]2(1) and [Ni(L)]n(2), respectively. These distances are shorter than the sum of their van der Waals radii (4.80 Å) but longer than the sum of their covalent radii (2.50 Å), indicating that there is a Ni⋯Ni interaction but not a Ni-Ni bond. The discrete molecules are π-stacked and connected via weak intermolecular interactions (C-H⋯O and C-H⋯N). Cyclic voltammetry measurements were obtained for both the complexes, and their pharmacokinetic and chemoinformatics properties were also explored. Detailed structural analysis and non-covalent supramolecular interactions were investigated using single-crystal structure analysis and computational approaches. Both the unique structures show good inhibition performance for the Omicron spike proteins of the SARS CoV-2 virus. To gain insights into potential SARS-CoV-2 Omicron drugs and find inhibitors against the Omicron variants of SARS-CoV-2, we examined the molecular docking of the nickel(ii) complexes [Ni(L)]2(1) and [Ni(L)]n(2) with the SARS-CoV-2 Omicron spike protein (PDB ID: 7WK2 and 7WVO). A strong binding was predicted between Ni(ii) coordination complexes [Ni(L)]2(1) and [Ni(L)]n(2) with the SARS-CoV-2 Omicron variant receptor protein through the negative value of binding affinity. Molecular docking of Nil(ii) complexes [Ni(L)]2(1) and [Ni(L)]n(2) with a DNA duplex (PDB ID: 7D3T) and RNA (PDB ID: 7TDC) binding protein was also studied. Overall, this study suggests that Ni(ii) complexes can be considered as drug candidates against the Omicron variants of SARS-CoV-2.

DockOpt: A Tool for Automatic Optimization of Docking Models

Molecular docking is a widely used technique for leveraging protein structure for ligand discovery, but it remains difficult to utilize due to limitations that have not been adequately addressed. Despite some progress toward automation, docking still requires expert guidance, hindering its adoption by a broader range of investigators. To make docking more accessible, we developed a new utility called DockOpt, which automates the creation, evaluation, and optimization of docking models prior to their deployment in large-scale prospective screens. DockOpt outperforms our previous automated pipeline across all 43 targets in the DUDE-Z benchmark data set, and the generated models for 84% of targets demonstrate sufficient enrichment to warrant their use in prospective screens, with normalized LogAUC values of at least 15%. DockOpt is available as part of the Python package Pydock3 included in the UCSF DOCK 3.8 distribution, which is available for free to academic researchers at https://dock.compbio.ucsf.edu and free for everyone upon registration at https://tldr.docking.org.

Evaluating the inhibitory activity of ferrocenyl Schiff bases derivatives on 5-lipoxygenase: Computational and biological studies

In the search for new 5-LOX inhibitors, two ferrocenyl Schiff base complexes functionalized with catechol ((ƞ5-(E)-C5H4-NCH-3,4-benzodiol)Fe(ƞ5-C5H5) (3a)) and vanillin ((ƞ5-(E)-C5H4-NCH-3-methoxy-4-phenol)Fe(ƞ5-C5H5) (3b)) were obtained. Complexes 3a and 3b were biologically evaluated as 5-LOX inhibitors showed potent inhibition compared to their organic analogs (2a and 2b) and known commercial inhibitors, with IC50 = 0.17 ± 0.05 μM for (3a) and 0.73 ± 0.06 μM for (3b) demonstrated a highly inhibitory and potent effect against 5-LOX due to the incorporation of the ferrocenyl fragment. Molecular dynamic studies showed a preferential orientation of the ferrocenyl fragment toward the non-heme iron of 5-LOX, which, together with electrochemical and in-vitro studies, allowed us to propose a competitive redox deactivation mechanism mediated by water, where Fe(III)-enzyme can be reduced by the ferrocenyl fragment. An Epa/IC50 relationship was observed, and the stability of the Schiff bases was evaluated by SWV in the biological medium, observing that the hydrolysis does not affect the high potency of the complexes, making them interesting alternatives for pharmacological applications.

New nickel(ii) Schiff base complexes as potential tools against SARS-CoV-2 Omicron target proteins: an in silico approach

Herein, we report the in silico design and synthesis of two new nickel(ii) coordination complexes, based on Schiff bases derived from the 2-hydroxy-1-naphthaldehyde moiety.

Design, synthesis and characterization of novel Ni(II) and Cu(II) complexes as antivirus drug candidates against SARS-CoV-2 and HIV virus

This paper describes the structure-based design, synthesis and anti-virus effect of two new coordination complexes, a Ni(II) complex [Ni(L)2] (1) and a Cu(II) complex [Cu(L)2] (2) of (E)-N-phenyl-2-(thiophen-2-ylmethylene) hydrazine-1-carbothioamide(HL). The synthesized ligand was coordinated to metal ions through the bidentate-N, S donor atoms. The newly synthesized complexes were characterized by various spectroscopic and physiochemical methods, powdered XRD analysis and also X-ray crystallography study. Ni(II) complex [Ni(L)2](1) crystallize in orthorhombic crystal system with the space group Pbca with four molecules in the unit cell (a = 9.857(3) Å, b = 7.749(2) Å, c = 32.292(10) Å, α = 90°, β = 90°, γ = 90°, Z= 4) and reveals a distorted square planar geometry. A Hirshfeld surface and 2D fingerprint plot has been explored in the crystal structure of Ni(II) complex [Ni(L)2] (1). Energy framework computational analysia has also been explored. DFT based calculations have been performed on the Schiff base and its metal complexes to study the structure-property relationship. Furthermore, the molecular docking studies of the ligand and its metal complexes with SARS-CoV-2 virus (PDB ID: 7BZ5) and HIV-1 virus (PDB ID: 6MQA) are also investigated. The molecular docking calculations of the Ni(II) complex [Ni(L)2] (1) and a Cu(II) complex [Cu(L)2] (2) with SARS-CoV-2 virus revealed that the binding affinities at inhibition binding site of receptor protein are 9.7 kcal/mol and -9.3 kcal/mol, respectively. The molecular docking results showed that the binding affinities of Ni(II) complex (1) and Cu(II) complex (2) against SARS-CoV-2 virus were found comparatively higher than the HIV-1 virus (-8.5 kcal/mol and -8.2 kcal/mol, respectively). As potential drug candidates, Swiss-ADME predictions analyses are also studied and the results are compared with Chloroquine (CQ) and Hydroxychloroquine (HCQ) as anti-SARS-CoV-2 drugs.

Dibenzyl (benzo [d] thiazol-2-yl (hydroxy) methyl) phosphonate (DBTMP) showing anti-S. aureus and anti-biofilm properties by elevating activities of serine protease (SspA) and cysteine protease staphopain B (SspB)

Staphylococcus aureus biofilms are the pathogenic factor in the spread of infection and are more pronounced in multidrug-resistant strains of S. aureus, where high expression of proteases is observed. Among various proteases, Serine protease (SspA) and cysteine protease Staphopain B (SspB) are known to play a key role in the biofilm formation and removal of biofilms. In earlier studies, we have reported Dibenzyl (benzo [d] thiazol-2-yl (hydroxy) methyl) phosphonate (DBTMP) exhibits anti-S. aureus and anti-biofilm properties by elevating the expression of the protease. In this study, the effect of DBTMP on the activities of SspA, and SspB of S. aureus was evaluated. The SspA and SspB genes of S. aureus ATCC12600 were sequenced (Genbank accession numbers: MZ456982 and MW574006). In S. aureus active SspA is formed by proteolytic cleavage of immature SspA, to get this mature SspA (mSspA), we have PCR amplified the mSspA sequence from the SspA gene. The mSspA and SspB genes were cloned, expressed, and characterized. The pure recombinant proteins rSspB and rmSspA exhibited a single band in SDS-PAGE with a molecular weight of 40 and 30 KD, respectively. The activities of rmSspA and rSspB are 32.33 and 35.45 Units/mL correspondingly. DBTMP elevated the activities of rmSspA and rSspB by docking with respective enzymes. This compound disrupted the biofilms formed by the multidrug-resistant strains of S. aureus and further prevented biofilm formation. These findings explain that DBTMP possesses anti-S. aureus and anti-biofilm features.

Hinokiflavone Inhibits MDM2 Activity by Targeting the MDM2-MDMX RING Domain

The proto-oncogene MDM2 is frequently amplified in many human cancers and its overexpression is clinically associated with a poor prognosis. The oncogenic activity of MDM2 is demonstrated by its negative regulation of tumor suppressor p53 and the substrate proteins involved in DNA repair, cell cycle control, and apoptosis pathways. Thus, inhibition of MDM2 activity has been pursued as an attractive direction for the development of anti-cancer therapeutics. Virtual screening was performed using the crystal structure of the MDM2-MDMX RING domain dimer against a natural product library and identified a biflavonoid Hinokiflavone as a promising candidate compound targeting MDM2. Hinokiflavone was shown to bind the MDM2-MDMX RING domain and inhibit MDM2-mediated ubiquitination in vitro. Hinokiflavone treatment resulted in the downregulation of MDM2 and MDMX and induction of apoptosis in various cancer cell lines. Hinokiflavone demonstrated p53-dependent and -independent tumor-suppressive activity. This report provides biochemical and cellular evidence demonstrating the anti-cancer effects of Hinokiflavone through targeting the MDM2-MDMX RING domain.

Drug Discovery of Plausible Lead Natural Compounds That Target the Insulin Signaling Pathway: Bioinformatics Approaches

The growing smooth talk in the field of natural compounds is due to the ancient and current interest in herbal medicine and their potentially positive effects on health. Dozens of antidiabetic natural compounds were reported and tested in vivo, in silico, and in vitro. The role of these natural compounds, their actions on the insulin signaling pathway, and the stimulation of the glucose transporter-4 (GLUT4) insulin-responsive translocation to the plasma membrane (PM) are all crucial in the treatment of diabetes and insulin resistance. In this review, we collected and summarized a group of available in vivo and in vitro studies which targeted isolated phytochemicals with possible antidiabetic activity. Moreover, the in silico docking of natural compounds with some of the insulin signaling cascade key proteins is also summarized based on the current literature. In this review, hundreds of recent studies on pure natural compounds that alleviate type II diabetes mellitus (type II DM) were revised. We focused on natural compounds that could potentially regulate blood glucose and stimulate GLUT4 translocation through the phosphoinositide 3-kinase (PI3K)/protein kinase B (Akt) pathway. On attempt to point out potential new natural antidiabetic compounds, this review also focuses on natural ingredients that were shown to interact with proteins in the insulin signaling pathway in silico, regardless of their in vitro/in vivo antidiabetic activity. We invite interested researchers to test these compounds as potential novel type II DM drugs and explore their therapeutic mechanisms.

Target-Based In Silico Screening for Phytoactive Compounds Targeting SARS-CoV-2

Coronavirus disease 2019 (COVID-19), resulting from infection by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), can cause severe and fatal pneumonia along with other life-threatening complications. The COVID-19 pandemic has taken a heavy toll on the healthcare system globally and has hit the economy hard in all affected countries. As a result, there is an unmet medical need for both the prevention and treatment of COVID-19 infection. Several herbal remedies have claimed to show promising clinical results, but the mechanisms of action are not clear. We set out to identify the anti-viral natural products of these herbal remedies that presumably inhibit the life cycle of SARS-CoV-2. Particularly we chose four key SARS-CoV-2 viral enzymes as targets: Papain-like protease, Main protease, RNA dependent RNA polymerase, and 2'-O-ribose methyltransferase, which were subjected to an unbiased in silico screening against a small molecule library of 33,765 compounds originating from herbs and medicinal plants. The small molecules were then ranked based on their free energy of fitting into the "druggable" pockets on the surface of each target protein. We have analyzed the best "fit" molecules and annotated them according to their plant sources and pharmacokinetic properties. Here we present a list of potential anti-viral ingredients of herbal remedies targeting SARS-CoV-2 and explore the potential mechanisms of action of these compounds as a framework for further development of chemoprophylaxis agents against COVID-19.

Modelling and targeting mitochondrial protein tyrosine phosphatase 1: a computational approach

The present research scintillates on the homology modelling of rat mitochondrial protein tyrosine phosphatase 1 (PTPMT1) and targeting its activity using flavonoids through a computational docking approach. PTPMT1 is a dual-specificity phosphatase responsible for protein phosphorylation and plays a vital role in the metabolism of cardiolipin biosynthesis, insulin regulation, etc. The inhibition of PTPMT1 has also shown enhanced insulin levels. The three-dimensional structure of the protein is not yet known. The homology modelling was performed using SWISS-MODEL and Geno3D webservers to compare the efficiencies. The PROCHECK for protein modelled using SWISS-MODEL showed 91.6% of amino acids in the most favoured region, 0.7% residues in the disallowed region that was found to be significant compared to the model built using Geno3D. 210 common flavonoids were docked in the modelled protein using the AutoDock 4.2.6 along with a control drug alexidine dihydrochloride. Our results show promising candidates that bind protein tyrosine phosphatase 1, including, prunin (- 8.66 kcal/mol); oroxindin (- 8.56 kcal/mol); luteolin 7-rutinoside (- 8.47 kcal/mol); 3(2H)-isoflavenes (- 8.36 kcal/mol); nicotiflorin (- 8.29 kcal/mol), ranked top in the docking experiments. We predicted the pharmacokinetic and Lipinski properties of the top ten compounds with the lowest binding energies. To further validate the stability of the modelled protein and docked complexes molecular dynamics simulations were performed using Desmond, Schrodinger for 150 ns in conjunction with MM-GBSA. Thus, flavonoids could act as potential inhibitors of PTPMT1, and further, in-vitro and in-vivo studies are essential to complete the drug development process.

Structure-based design and synthesis of copper(ii) complexes as antivirus drug candidates targeting SARS CoV-2 and HIV

This paper describes the structure-based design and synthesis of two novel square-planar trans-N2O2 Cu(ii) complexes [Cu(L1)2] (1) and [Cu(L2)2] (2) of 2-((Z)-(4-methoxyphenylimino)methyl)-4,6-dichlorophenol (L1H) and 2-((Z)-(2,4-dibromophenylimino)methyl)-4-bromophenol (L2H) as potential inhibitors against the main protease of the SARS-CoV-2 and HIV viruses.

The effect of ligand affinity to the contact dynamics of the ligand binding domain of thyroid hormone receptor - retinoid X receptor

Ligand-based allostery has been gaining attention for its importance in protein regulation and implication in drug design. One of the interesting cases of protein allostery is the thyroid hormone receptor - retinoid x receptor (TR:RXR), which regulates the gene expression of important physiological processes, such as development and metabolism. It is regulated by the TR native ligand triiodothyronine (T3), which displays anticooperative behavior to the RXR ligand 9-cis retinoic acid (9C). In contrast to this anticooperative behavior, 9C has been shown to increase the activity of TR:RXR. Here we probed the influence of the affinity and the interactions of the TR ligand to the allostery of the TR:RXR through contact dynamics and residue networks. The TR ligand analogs were designed to have higher (G2) and lower (N1) binding energies than T3 when docked to the TR:RXR(9C) complex. The aqueous TR(N1/T3/G2):RXR(9C) complexes were subjected to 30 ns all-atom simulations using theNAMD. The program CAMERRA was used to capture the subtle perturbations of TR:RXR by mapping the residue contact dynamics. Various parts of the TR ligands; including the hydrophilic head, the iodine substituents, and the ligand tail; have been probed for their significance in ligand affinity. The results on the T3 and G2 complexes suggest that ligand affinity can be utilized as a predictor for anticooperative systems on which ligand is more likely to dissociate or remain bound. All 3 complexes also display distinct contact networks for cross-dimer signalling and ligand communication. Understanding ligand-based allostery could potentially unveil secrets of ligand-regulated protein dynamics, a foundation for the design of better and more efficient allosteric drugs.

Rational Computational Design of Fourth-Generation EGFR Inhibitors to Combat Drug-Resistant Non-Small Cell Lung Cancer

Although the inhibitors of singly mutated epidermal growth factor receptor (EGFR) kinase are effective for the treatment of non-small cell lung cancer (NSCLC), their clinical efficacy has been limited due to the emergence of various double and triple EGFR mutants with drug resistance. It has thus become urgent to identify potent and selective inhibitors of triple mutant EGFRs resistant to first-, second-, and third-generation EGFR inhibitors. Herein, we report the discovery of potent and highly selective inhibitors of EGFR exon 19 p.E746_A750del/EGFR exon 20 p.T790M/EGFR exon 20 p.C797S (d746-750/T790M/C797S) mutant, which were derived via two-track virtual screening and de novo design. This two-track approach was performed so as to maximize and minimize the inhibitory activity against the triple mutant and the wild type, respectively. Extensive chemical modifications of the initial hit compounds led to the identification of several low-nanomolar inhibitors of the d746-750/T790M/C797S mutant. Among them, two compounds exhibited more than 10⁴-fold selectivity in the inhibition of EGFR^{d746-750/T790M/C797S} over the wild type. The formations of a hydrogen bond with the mutated residue Ser797 and the van der Waals contact with the mutated residue Met790 were found to be a common feature in the interactions between EGFR^{d746-750/T790M/C797S} and the fourth-generation inhibitors. Such an exceptionally high selectivity could also be attributed to the formation of the hydrophobic contact with a Gly loop residue or the hydrogen bond with Asp855 in the activation loop. The discovery of the potent and selective EGFR^{d746-750/T790M/C797S} inhibitors were actually made possible by virtue of the modified protein-ligand binding free energy function involving a new hydration free energy term with enhanced accuracy. The fourth-generation EGFR inhibitors found in this work are anticipated to serve as a new starting point for the discovery of anti-NSCLC medicines to overcome the problematic drug resistance.

Atomic Charge Calculator II: web-based tool for the calculation of partial atomic charges

Partial atomic charges serve as a simple model for the electrostatic distribution of a molecule that drives its interactions with its surroundings. Since partial atomic charges are frequently used in computational chemistry, chemoinformatics and bioinformatics, many computational approaches for calculating them have been introduced. The most applicable are fast and reasonably accurate empirical charge calculation approaches. Here, we introduce Atomic Charge Calculator II (ACC II), a web application that enables the calculation of partial atomic charges via all the main empirical approaches and for all types of molecules. ACC II implements 17 empirical charge calculation methods, including the highly cited (QEq, EEM), the recently published (EQeq, EQeq+C), and the old but still often used (PEOE). ACC II enables the fast calculation of charges even for large macromolecular structures. The web server also offers charge visualization, courtesy of the powerful LiteMol viewer. The calculation setup of ACC II is very straightforward and enables the quick calculation of high-quality partial charges. The application is available at https://acc2.ncbr.muni.cz.

Citronellal perception and transmission by Anopheles gambiae s.s. (Diptera: Culicidae) females

Anopheles gambiae s.s. is a key vector of Plasmodium parasites. Repellents, which may be a promising alternative to pesticides used to control malaria mosquitoes. Although citronellal is a known mosquito repellent, its repellency characteristics are largely unknown. Determining the specific odorant-binding proteins (OBPs) and odorant receptors (ORs) that detect and transfer the citronellal molecule in A. gambiae s.s. will help to define the mode of action of this compound. In this research, we assessed the repellent activity of citronellal in A. gambiae s.s. using a Y-tube olfactory meter, screened candidate citronellal-binding OBPs and ORs using reverse molecular docking, clarified the binding properties of predicted proteins for citronellal using fluorescence competition binding assay. Results showed that citronellal had a dosage effect on repelling A. gambiae s.s.. The 50% repellent rate was determined to be 4.02 nmol. Results of simulated molecular docking showed that the only proteins that bound tightly with citronellal were AgamOBP4 and AgamORC7. Fluorescence competitive binding assays confirmed the simulations. This research determined that citronellal was captured by AgamOBP4 and transmitted to AgamORC7 in A. gambiae s.s.. Our study will be beneficial in the further understanding the repellent mechanism of citronellal against A. gambiae s.s..

Interleukin-26 activates macrophages and facilitates killing of Mycobacterium tuberculosis

Tuberculosis-causing Mycobacterium tuberculosis (Mtb) is transmitted via airborne droplets followed by a primary infection of macrophages and dendritic cells. During the activation of host defence mechanisms also neutrophils and T helper 1 (T_H1) and T_H17 cells are recruited to the site of infection. The T_H17 cell-derived interleukin (IL)-17 in turn induces the cathelicidin LL37 which shows direct antimycobacterial effects. Here, we investigated the role of IL-26, a T_H1- and T_H17-associated cytokine that exhibits antimicrobial activity. We found that both IL-26 mRNA and protein are strongly increased in tuberculous lymph nodes. Furthermore, IL-26 is able to directly kill Mtb and decrease the infection rate in macrophages. Binding of IL-26 to lipoarabinomannan might be one important mechanism in extracellular killing of Mtb. Macrophages and dendritic cells respond to IL-26 with secretion of tumor necrosis factor (TNF)-α and chemokines such as CCL20, CXCL2 and CXCL8. In dendritic cells but not in macrophages cytokine induction by IL-26 is partly mediated via Toll like receptor (TLR) 2. Taken together, IL-26 strengthens the defense against Mtb in two ways: firstly, directly due to its antimycobacterial properties and secondly indirectly by activating innate immune mechanisms.

Potential therapeutic target identification in the novel 2019 coronavirus: insight from homology modeling and blind docking study

Background

The 2019-nCoV which is regarded as a novel coronavirus is a positive-sense single-stranded RNA virus. It is infectious to humans and is the cause of the ongoing coronavirus outbreak which has elicited an emergency in public health and a call for immediate international concern has been linked to it. The coronavirus main proteinase which is also known as the 3C-like protease (3CLpro) is a very important protein in all coronaviruses for the role it plays in the replication of the virus and the proteolytic processing of the viral polyproteins. The resultant cytotoxic effect which is a product of consistent viral replication and proteolytic processing of polyproteins can be greatly reduced through the inhibition of the viral main proteinase activities. This makes the 3C-like protease of the coronavirus a potential and promising target for therapeutic agents against the viral infection.

Results

This study describes the detailed computational process by which the 2019-nCoV main proteinase coding sequence was mapped out from the viral full genome, translated and the resultant amino acid sequence used in modeling the protein 3D structure. Comparative physiochemical studies were carried out on the resultant target protein and its template while selected HIV protease inhibitors were docked against the protein binding sites which contained no co-crystallized ligand.

Conclusion

In line with results from this study which has shown great consistency with other scientific findings on coronaviruses, we recommend the administration of the selected HIV protease inhibitors as first-line therapeutic agents for the treatment of the current coronavirus epidemic.

Binding Affinity Prediction by Pairwise Function Based on Neural Network

We present a new approach to estimate the binding affinity from given three-dimensional poses of protein-ligand complexes. In this scheme, every protein-ligand atom pair makes an additive free-energy contribution. The sum of these pairwise contributions then gives the total binding free energy or the logarithm of the dissociation constant. The pairwise contribution is calculated by a function implemented via a neural network that takes the properties of the two atoms and their distance as input. The pairwise function is trained using a portion of the PDBbind 2018 data set. The model achieves good accuracy for affinity predictions when evaluated with PDBbind 2018 and with the CASF-2016 benchmark, comparing favorably to many scoring functions such as that of AutoDock Vina. The framework here may be extended to incorporate other factors to further improve its accuracy and power.

Phenanthridine derivatives as promising new anticancer agents: synthesis, biological evaluation and binding studies

Aim: Phenanthridines are an essential class of nitrogenous heterocycles with extensive applications in medicinal chemistry. The development of efficient and eco-friendly methods for the preparation of chirally pure dihydropyrrolo[1,2-f]phenanthridines (5a-h), and their in vitro evaluation and modeling studies as potential anticancer, antioxidant and DNA cleavage agents is reported. Methodology & results: Compounds 5a-h were prepared through a facile one-pot synthesis and characterized by infrared, high resolution mass spectrometry, ¹H and ¹³C nuclear magnetic resonance. The molecules were subjected to virtual screening and docking analysis against selected human molecular targets. Compound 5g displayed good binding properties as well as significant anticancer and DNA cleavage activity. Conclusion: Compound 5g has been identified as a potential lead candidate for further testing against additional cancer cell lines and animal models in future.

Conjugation of a gold(iii) complex with vitamin B1 and chlorambucil derivatives: anticancer evaluation and mechanistic insights

A gold(iii) complex containing vitamin B1 and chlorambucil derivatives was investigated for mechanistic insights in colon and breast cancer treatment.

Three-Dimensional Analysis of Binding Sites for Predicting Binding Affinities in Drug Design

Understanding the interaction between drug molecules and proteins is one of the main challenges in drug design. Several tools have been developed recently to decrease the complexity of the process. Artificial intelligence and machine learning methods offer promising results in predicting the binding affinities. It becomes possible to do accurate predictions by using the known protein-ligand interactions. In this study, the electrostatic potential values extracted from 3-dimensional grid cubes of the drug-protein binding sites are used for predicting binding affinities of related complexes. A new algorithm with a dynamic feature selection method was implemented, which is derived from Compressed Images For Affinity Prediction (CIFAP) study, to predict binding affinities of Checkpoint Kinase 1 and Caspase 3 inhibitors.

Molecular docking analysis of apigenin and quercetin from ethylacetate fraction of Adansonia digitata with malaria-associated calcium transport protein: An in silico approach

Background

The investigation and knowledge of calcium handling mechanisms in the plasmodium has been considered as a potential biological target against malaria.

Objective

This study deals with the evaluation of inhibitory activity of secondary metabolites of ethylacetate partitioned-fraction of Adansonia digitata stem bark extract on malaria-associated protein using in silico docking studies.

Materials and methods

Molecular docking and virtual screening was performed to understand the mechanism of ligand binding and to identify potent calcium transporter inhibitors. The stem bark extracts of A. digitata contains rich sources of phytochemicals. The secondary metabolites were determined by HPLC-DAD and HRGC-MS analysis. The major chemical constituent present in the ethylacetate partitioned-fraction of A. digitata stem bark extract were examined for their antiplasmodial activity and were also involved in docking study.

Results

The secondary metabolites, quercetin and apigenin inhibited the formation of β-hematin. The results showed that all the selected compounds in the A. digitata showed binding energy ranging between -6.5 kcal/mol and -7.1 kcal/mol. Among the two chemical constituents, apigenin has the highest docking score along with the highest number of hydrogen bonds formed when compared to quercetin. Analysis of the results suggests that apigenin and quercetin could act as an anti-malaria agent.

Conclusion

Molecular docking analysis could lead to further development of potent calcium transporter inhibitors for the prevention and treatment of malaria and related conditions.

Synthesis and biological evaluation of acylthiourea against DUSP1 inhibition

Structure based virtual screening attempts to discover DUSP1 inhibitors have yielded a scaffold featuring benzoxazole and acylthiourea pharmacophore. A series of its analogues were synthesized to explore structure activity relationship (SAR) of DUSP1 inhibition.

Experimental and theoretical investigations of novel oxidovanadium(IV) juglone complex: DNA/HSA interaction and cytotoxic activity

A new oxidovanadium(IV) complex VO(L)(Jug) (HL = 5-methoxy-1,3-bis (1-methyl-1H-benzo[d]imidazol-2-yl)benzene, Jug = juglone) was synthesized and characterized. Interactions of the V(IV) complex with calf thymus DNA (CT DNA) and human serum albumin were studied using different techniques such as UV-vis and fluorescence emission spectroscopy. The experimental results were confirmed by the molecular docking study. The oxidovanadium(IV) complex can efficiently cleave pUC19 DNA in the presence of Hydrogen peroxide. Also, the in vitro cytotoxicity properties of the oxidovanadium(IV) complex was evaluated against MCF-7, HPG-2 and HT-29 cancer cell lines and HEK293 non-malignant fibroblasts were evaluated and compared with free ligands, VOSO₄ and cisplatin as reference drugs.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulation analysis of conessine against multi drug resistant Serratia marcescens

Ornithine decarboxylase (ODC) is an immediate precursor of polyamine biosynthesis in Serratia marcescens and a potential target for inhibition of its growth. We predicted the 3D structural conformation of ODC enzyme and validated it using MDS in our previous study. In this current study, the potential inhibitors of ODC were obtained by virtual screening of potential inhibitors from ZINC database and studied in depth for their different binding pose. Among the ten virtually screened inhibitors, Conessine exhibited the best binding with ODC and its inhibition property was studied further by MDS studies. The natural compound conessine is isolated from plant Holarrhena antidysenterica and it is studied against ODC of Serratia marcenses for its inhibitory potentials. This revealed unforeseen twisted position in root mean square fluctuation (RMSF) and ODC modelled conformation that influenced ligand binding. Both predicted model and ligand bound model were compared and found to be stable with Root Mean Square Deviation (RMSD) of approximately 7 nm and 0.25 nm to that of crystallographic structure over simulation time of 55 ns and 70 ns respectively. This work paves the way for future development of new drugs against nosocomial diseases caused by Serratia marcescens.

Shining Light on Molecular Mechanism for Odor-selectivity of CNT-immobilized Olfactory Receptor

Olfactory receptor (OR)-based bioelectronic nose is a new type of bio-affinity sensor applied for detecting numerous odorant molecules. In order to elucidate the effect of the adsorption of nanomaterial carriers on the receptor structure and its selectivity to odors, we used a systematic computation-scheme to study two OR models immobilized onto carbon nanotube. Our result indicates that there is a multistep OR-adsorption process driven by hydrophobic interaction. Many allosteric communication pathways exist between the absorbed residues and the pocket ones, leading to a significant shrinkage of the pocket. Consequently, the size-selectivity of the receptor to the odors is changed to some extent. But, the odor size and its hydrophobicity, rather than specific functional groups of the odor, still play a determinant role in binding OR, at least for the 132 odors under study. Regardless of the limitation for the odor size in initial recognition, the different-size odors could induce significant changes in the pocket conformation so that it could better match the pocket space, indicating the importance of the ligand-fit binding. Due to the CNT-induced shrinkage of the pocket, the CNT immobilization could increase the binding affinity through enhancing van der Waals interaction, in particular for the large odors.

Role of proton balance in formation of self-assembled chitosan nanoparticles

Researchers have explored the ability of chitosan to form nanoparticles, to suit varying applications, ranging from wound-healing to gene delivery. Ionic gelation is a widely used method for formulating chitosan nanoparticles, where self-assembly plays a crucial role. This self-assembly is initially promoted by hydrophilic-hydrophobic parity amongst individual chitosan residues, along with electrostatic and Van der Waals interactions with the cross-linker. However, until now the intrinsic ability of chitosan to self-assemble is not widely studied; hence, we investigate the self-assembly of chitosan, based on proton balance between its protonated and deprotonated residues, to promote facile nanoparticle synthesis. This is one of the first reports that highlights subtle but critical influence of proton balance in the chitosan polymer on the formation of chitosan nanoparticles.

New HSP27 inhibitors efficiently suppress drug resistance development in cancer cells

Drug resistance is an important open problem in cancer treatment. In recent years, the heat shock protein HSP27 (HSPB1) was identified as a key player driving resistance development. HSP27 is overexpressed in many cancer types and influences cellular processes such as apoptosis, DNA repair, recombination, and formation of metastases. As a result cancer cells are able to suppress apoptosis and develop resistance to cytostatic drugs. To identify HSP27 inhibitors we follow a novel computational drug repositioning approach. We exploit a similarity between a predicted HSP27 binding site to a viral thymidine kinase to generate lead inhibitors for HSP27. Six of these leads were verified experimentally. They bind HSP27 and down-regulate its chaperone activity. Most importantly, all six compounds inhibit development of drug resistance in cellular assays. One of the leads – chlorpromazine – is an antipsychotic, which has a positive effect on survival time in human breast cancer. In summary, we make two important contributions: First, we put forward six novel leads, which inhibit HSP27 and tackle drug resistance. Second, we demonstrate the power of computational drug repositioning.

Multifunctional iminochromene-2H-carboxamide derivatives containing different aminomethylene triazole with BACE1 inhibitory, neuroprotective and metal chelating properties targeting Alzheimer's disease

Alzheimer's disease (AD) is a neurodegenerative disorder known for the presence of amyloid beta plaques resulting from the sequential action of β-secretase and γ-secretase on amyloid precursor protein. We developed and synthesized, through click reactions, a new family of iminochromene carboxamides containing different aminomethylene triazole. The BACE1 inhibition, neuroprotective capacity and metal chelation of these derivatives make them ideal candidates against AD. Most of the synthesized compounds were shown to have potent BACE1 inhibitory activity in a FRET assay, with an IC₅₀ value of 2.2 μM for the most potent compound. Moreover, molecular modeling evaluation of these BACE1 inhibitors demonstrates the vital role of the amine and amide linkers through hydrogen bond interactions with key amino acids in the BACE1 active site. Our in vitro neuroprotective evaluations in PC12 neuronal cells of Aβ-induced neuroprotection demonstrated promising activity for most of the compounds as neuroprotective agents. Based on our findings, we propose that introduction of a phthalimide substitute on the triazole ring shown to be interesting multifunctional lead compound worthy of further study.

Two-track virtual screening approach to identify both competitive and allosteric inhibitors of human small C-terminal domain phosphatase 1

Despite a wealth of persuasive evidence for the involvement of human small C-terminal domain phosphatase 1 (Scp1) in the impairment of neuronal differentiation and in Huntington's disease, small-molecule inhibitors of Scp1 have been rarely reported so far. This study aims to the discovery of both competitive and allosteric Scp1 inhibitors through the two-track virtual screening procedure. By virtue of the improvement of the scoring function by implementing a new molecular solvation energy term and by reoptimizing the atomic charges for the active-site Mg²⁺ ion cluster, we have been able to identify three allosteric and five competitive Scp1 inhibitors with low-micromolar inhibitory activity. Consistent with the results of kinetic studies on the inhibitory mechanisms, the allosteric inhibitors appear to be accommodated in the peripheral binding pocket through the hydrophobic interactions with the nonpolar residues whereas the competitive ones bind tightly in the active site with a direct coordination to the central Mg²⁺ ion. Some structural modifications to improve the biochemical potency of the newly identified inhibitors are proposed based on the binding modes estimated with docking simulations.

An experimental and theoretical study on the interaction of DNA and BSA with novel Ni2+, Cu2+ and VO2+ complexes derived from vanillin bidentate Schiff base ligand

In this investigation, the structure of bidentate N,N-Schiff base ligand of vanillin, (E)-4-(((2-amino-5-nitrophenyl)imino)methyl)-2-methoxyphenol (HL) was determined by single crystal X-ray diffraction. The interaction of new [CuL₂], [NiL₂] and [VOL₂] complexes with DNA and BSA was explored through UV-Vis and fluorescence spectroscopy. The electronic spectra changes displayed an isosbestic point for the complexes upon titration with DNA. The K_b values for the complexes [CuL₂], [NiL₂] and [VOL₂] were 2.4×10⁵, 1.9×10⁵ and 4.2×10⁴, respectively. [CuL₂] complex was bound more toughly than [NiL₂] and [VOL₂] complexes. These complexes had a significant interaction with Bovine Serum Albumin (BSA) and the results demonstrated that the quenching mechanism was a static procedure. Also, the complexes interacted with BSA by more than one binding site (n>1). Finally, the theoretical studies were performed using the docking method to calculate the binding constants and recognize the binding site of the DNA and BSA with the complexes. The ligand and complexes including Ni²⁺, Cu²⁺ and VO²⁺ ions were colonized by fungal growth.

Molecular mechanism of carbon nanotube to activate Subtilisin Carlsberg in polar and non-polar organic media

In the work, we mainly used molecular dynamics (MD) simulation and protein structure network (PSN) to study subtilisin Carlsberg (SC) immobilized onto carbon nanotube (CNT) in water, acetonitrile and heptane solvents, in order to explore activation mechanism of enzymes in non-aqueous media. The result indicates that the affinity of SC with CNT follows the decreasing order of water > acetonitrile > heptane. The overall structure of SC and the catalytic triad display strong robustness to the change of environments, responsible for the activity retaining. However, the distances between two β-strands of substrate-binding pocket are significantly expanded by the immobilization in the increasing order of water < acetonitrile < heptane, contributing to the highest substrate-binding energy in heptane media. PSN analysis further reveals that the immobilization enhances structural communication paths to the substrate-binding pocket, leading to its larger change than the free-enzymes. Interestingly, the increase in the number of the pathways upon immobilization is not dependent on the absorbed extent but the desorbed one, indicating significant role of shifting process of experimental operations in influencing the functional region. In addition, some conserved and important hot-residues in the paths are identified, providing molecular information for functional modification.

NEEMP: software for validation, accurate calculation and fast parameterization of EEM charges

Background

The concept of partial atomic charges was first applied in physical and organic chemistry and was later also adopted in computational chemistry, bioinformatics and chemoinformatics. The electronegativity equalization method (EEM) is the most frequently used approach for calculating partial atomic charges. EEM is fast and its accuracy is comparable to the quantum mechanical charge calculation method for which it was parameterized. Several EEM parameter sets for various types of molecules and QM charge calculation approaches have been published and new ones are still needed and produced. Methodologies for EEM parameterization have been described in a few articles, but a software tool for EEM parameterization and EEM parameter sets validation has not been available until now.

Results

We provide the software tool NEEMP (http://ncbr.muni.cz/NEEMP), which offers three main functionalities: EEM parameterization [via linear regression (LR) and differential evolution with local minimization (DE-MIN)]; EEM parameter set validation (i.e., validation of coverage and quality) and EEM charge calculation. NEEMP functionality is shown using a parameterization and a validation case study. The parameterization case study demonstrated that LR is an appropriate approach for smaller and homogeneous datasets and DE-MIN is a suitable solution for larger and heterogeneous datasets. The validation case study showed that EEM parameter set coverage and quality can still be problematic. Therefore, it makes sense to verify the coverage and quality of EEM parameter sets before their use, and NEEMP is an appropriate tool for such verification. Moreover, it seems from both case studies that new EEM parameterizations need to be performed and new EEM parameter sets obtained with high quality and coverage for key structural databases.

Conclusion

We provide the software tool NEEMP, which is to the best of our knowledge the only available software package that enables EEM parameterization and EEM parameter set validation. Additionally, its DE-MIN parameterization method is an innovative approach, developed by ourselves and first published in this work. In addition, we also prepared four high-quality EEM parameter sets tailored to ligand molecules.Graphical abstract

Identification of small molecules that inhibit the histone chaperone Asf1 and its chromatin function

The eukaryotic genome is packed into chromatin, which is important for the genomic integrity and gene regulation. Chromatin structures are maintained through assembly and disassembly of nucleosomes catalyzed by histone chaperones. Asf1 (anti-silencing function 1) is a highly conserved histone chaperone that mediates histone transfer on/off DNA and promotes histone H3 lysine 56 acetylation at globular core domain of histone H3. To elucidate the role of Asf1 in the modulation of chromatin structure, we screened and identified small molecules that inhibit Asf1 and H3K56 acetylation without affecting other histone modifications. These pyrimidine-2,4,6-trione derivative molecules inhibited the nucleosome assembly mediated by Asf1 in vitro, and reduced the H3K56 acetylation in HeLa cells. Furthermore, production of HSV viral particles was reduced by these compounds. As Asf1 is implicated in genome integrity, cell proliferation, and cancer, current Asf1 inhibitor molecules may offer an opportunity for the therapeutic development for treatment of diseases. [BMB Reports 2015; 48(12): 685-690]

The Ascaris suum nicotinic receptor, ACR-16, as a drug target: Four novel negative allosteric modulators from virtual screening

Soil-transmitted helminth infections in humans and livestock cause significant debility, reduced productivity and economic losses globally. There are a limited number of effective anthelmintic drugs available for treating helminths infections, and their frequent use has led to the development of resistance in many parasite species. There is an urgent need for novel therapeutic drugs for treating these parasites. We have chosen the ACR-16 nicotinic acetylcholine receptor of Ascaris suum (Asu-ACR-16), as a drug target and have developed three-dimensional models of this transmembrane protein receptor to facilitate the search for new bioactive compounds. Using the human α7 nAChR chimeras and Torpedo marmorata nAChR for homology modeling, we defined orthosteric and allosteric binding sites on the Asu-ACR-16 receptor for virtual screening. We identified four ligands that bind to sites on Asu-ACR-16 and tested their activity using electrophysiological recording from Asu-ACR-16 receptors expressed in Xenopus oocytes. The four ligands were acetylcholine inhibitors (SB-277011-A, IC₅₀, 3.12 ± 1.29 μM; (+)-butaclamol Cl, IC₅₀, 9.85 ± 2.37 μM; fmoc-1, IC₅₀, 10.00 ± 1.38 μM; fmoc-2, IC₅₀, 16.67 ± 1.95 μM) that behaved like negative allosteric modulators. Our work illustrates a structure-based in silico screening method for seeking anthelmintic hits, which can then be tested electrophysiologically for further characterization.

•

Three-dimensional structural models of the Ascaris nicotinic (Asu-ACR-16) receptor made by homology modeling.

•

High affinity ligands selected by in silico screening.

•

Four ligands validated by electrophysiological studies as negative allosteric modulators.

High-quality and universal empirical atomic charges for chemoinformatics applications

BACKGROUND

Partial atomic charges describe the distribution of electron density in a molecule and therefore provide clues to the chemical behaviour of molecules. Recently, these charges have become popular in chemoinformatics, as they are informative descriptors that can be utilised in pharmacophore design, virtual screening, similarity searches etc. Especially conformationally-dependent charges perform very successfully. In particular, their fast and accurate calculation via the Electronegativity Equalization Method (EEM) seems very promising for chemoinformatics applications. Unfortunately, published EEM parameter sets include only parameters for basic atom types and they often miss parameters for halogens, phosphorus, sulphur, triple bonded carbon etc. Therefore their applicability for drug-like molecules is limited.

RESULTS

We have prepared six EEM parameter sets which enable the user to calculate EEM charges in a quality comparable to quantum mechanics (QM) charges based on the most common charge calculation schemes (i.e., MPA, NPA and AIM) and a robust QM approach (HF/6-311G, B3LYP/6-311G). The calculated EEM parameters exhibited very good quality on a training set ([Formula: see text]) and also on a test set ([Formula: see text]). They are applicable for at least 95 % of molecules in key drug databases (DrugBank, ChEMBL, Pubchem and ZINC) compared to less than 60 % of the molecules from these databases for which currently used EEM parameters are applicable.

CONCLUSIONS

We developed EEM parameters enabling the fast calculation of high-quality partial atomic charges for almost all drug-like molecules. In parallel, we provide a software solution for their easy computation (http://ncbr.muni.cz/eem_parameters). It enables the direct application of EEM in chemoinformatics.