Designing, docking and molecular dynamics simulation studies of novel cloperastine analogues as anti-allergic agents: homology modeling and active site prediction for the human histamine H1 receptor

G-protein signaling of oxytocin receptor as a potential target for cabazitaxel-resistant prostate cancer

Although the treatment armamentarium for patients with metastatic prostate cancer has improved recently, treatment options after progression on cabazitaxel (CBZ) are limited. To identify the mechanisms underlying CBZ resistance and therapeutic targets, we performed single-cell RNA sequencing of circulating tumor cells (CTCs) from patients with CBZ-resistant prostate cancer. Cells were clustered based on gene expression profiles. In silico screening was used to nominate candidate drugs for overcoming CBZ resistance in castration-resistant prostate cancer. CTCs were divided into three to four clusters, reflecting intrapatient tumor heterogeneity in refractory prostate cancer. Pathway analysis revealed that clusters in two cases showed up-regulation of the oxytocin (OXT) receptor-signaling pathway. Spatial gene expression analysis of CBZ-resistant prostate cancer tissues confirmed the heterogeneous expression of OXT-signaling molecules. Cloperastine (CLO) had significant antitumor activity against CBZ-resistant prostate cancer cells. Mass spectrometric phosphoproteome analysis revealed the suppression of OXT signaling specific to CBZ-resistant models. These results support the potential of CLO as a candidate drug for overcoming CBZ-resistant prostate cancer via the inhibition of OXT signaling.

Synthesis, in silico and in vitro antimicrobial efficacy of substituted arylidene-based quinazolin-4(3H)-one motifs

Introduction: Quinazolin-4(3H)-one derivatives have attracted considerable attention in the pharmacological profiling of therapeutic drug targets. The present article reveals the development of arylidene-based quinazolin-4(3H)-one motifs as potential antimicrobial drug candidates. Methods: The synthetic pathway was initiated through thermal cyclization of acetic anhydride on anthranilic acid to produce 2-methyl-4H-3,1-benzoxazan-4-one 1, which (upon condensation with hydrazine hydrate) gave 3-amino-2-methylquinazolin-4(3H)-one 2. The reaction of intermediate 2 at its amino side arm with various benzaldehyde derivatives furnished the final products, in the form of substituted benzylidene-based quinazolin-4(3H)-one motifs 3a-l, and with thiophene-2-carbaldehyde to afford 3 m. The purified targeted products 3a-m were effectively characterized for structural authentication using physicochemical parameters, microanalytical data, and spectroscopic methods, including IR, UV, and 1H- and 13C-NMR, as well as mass spectral data. The substituted arylidene-based quinazolin-4(3H)-one motifs 3a-m were screened for both in silico and in vitro antimicrobial properties against selected bacteria and fungi. The in silico studies carried out consisted of predicted ADMET screening, molecular docking, and molecular dynamics (MD) simulation studies. Furthermore, in vitro experimental validation was performed using the agar diffusion method, and the standard antibacterial and antifungal drugs used were gentamicin and ketoconazole, respectively. Results and discussion: Most of the compounds possessed good binding affinities according to the molecular docking studies, while MD simulation revealed their levels of structural stability in the protein-ligand complexes. 2-methyl-3-((thiophen-2-ylmethylene)amino) quinazolin-4(3H)-one 3 m emerged as both the most active antibacterial agent (with an minimum inhibitory concentration (MIC) value of 1.95 μg/mL) against Staphylococcus aureus and the most active antifungal agent (with an MIC value of 3.90 μg/mL) against Candida albicans, Aspergillus niger, and Rhizopus nigricans.

Molecular docking and in vitro analysis of peptides from Stolephorus indicus with ACE2

Peptides from Stolephorus indicus (Anchovies) meat lysate were generated using a Bacillus subtilis cysteine protease. The peptides were generated by enzyme hydrolysis after which the hydrolysate containing peptides were analysed by LC-MS/MS. Computer aided analysis of peptides using CASTp server and GOLD software show four peptides having ACE2 inhibitory activity. Further, peptides 1 (8 amino acids), 2 (8 amino acid), 5 (9 amino acids) and 11 (12 amino acids) showed good docking features for binding to ACE2 enzyme active sites, mainly by hydrogen bonding. Peptide 1 (8 amino acids-octa-peptide) having the highest docking score was tested in vitro for ACE2 binding and showed up to 40 % inhibition of ACE2 activity at a concentration of 10mM. Hence, this octa-peptide has a potential role in applications involving ACE2 inhibition thereby leading to the prevention of binding of spike glycoprotein to ACE2 receptor.

Molecular mechanism of inhibition of COVID-19 main protease by β-adrenoceptor agonists and adenosine deaminase inhibitors using in silico methods

Novel coronavirus (COVID-19) responsible for viral pneumonia which emerged in late 2019 has badly affected the world. No clinically proven drugs are available yet as the targeted therapeutic agents for the treatment of this disease. The viral main protease which helps in replication and transcription inside the host can be an effective drug target. In the present study, we aimed to discover the potential of β-adrenoceptor agonists and adenosine deaminase inhibitors which are used in asthma and cancer/inflammatory disorders, respectively, as repurposing drugs against protease inhibitor by ligand-based and structure-based virtual screening using COVID-19 protease-N3 complex. The AARRR pharmacophore model was used to screen a set of 22,621 molecules to obtain hits, which were subjected to high-throughput virtual screening. Extra precision docking identified four top-scored molecules such as +/--fenoterol, FR236913 and FR230513 with lower binding energy from both categories. Docking identified three major hydrogen bonds with Gly143, Glu166 and Gln189 residues. 100 ns MD simulation was performed for four top-scored molecules to analyze the stability, molecular mechanism and energy requirements. MM/PBSA energy calculation suggested that van der Waals and electrostatic energy components are the main reasons for the stability of complexes. Water-mediated hydrogen bonds between protein-ligand and flexibility of the ligand are found to be responsible for providing extra stability to the complexes. The insights gained from this combinatorial approach can be used to design more potent and bio-available protease inhibitors against novel coronavirus.Communicated by Ramaswamy H. Sarma.

Proteomic analysis of adipose tissue revealing differentially abundant proteins in highly efficient mid-lactating dairy cows

The improvement of nutrient utilization efficiency in dairy cows represents an important task in view of the current rising demand for animal products and sustainable resource usage. In this perspective, the identification of appropriate markers to identify the most efficient animals for dairy production becomes a crucial factor. Residual feed intake (RFI), which represents the difference between predicted and actual intake, is used to define the efficiency of cows. In this study, subcutaneous adipose tissue (AT) was collected from five high efficient (HEF) and five low efficient (LEF) mid-lactation Holstein dairy cows, that represented subgroups of the 20% lowest RFI values (HEF) and highest 20% RFI values (LEF), out of a cohort of 155 cows that were examined for feed efficiency at the individual dairy barn at Volcani Institute, Israel. Adipose samples were examined for proteomic analysis by nano-LC/MS–MS and gene expression by RT-PCR. A total of 101 differential proteins (P ≤ 0.05 and fold change ± 1.5) and two protein networks related to feed efficiency were found between HEF and LEF cows. Among the enriched top canonical pathways, FAT10 signaling, EIF2 signaling, Sirtuin signaling, Acute phase response signaling, Protein ubiquitination and mTOR signaling pathways were related to feed efficiency in AT. Furthermore, abundance of transferrin (TF; FC = 78.35, P = 0.02) enriched pathways, including mTOR signaling, LXR/RXR and FXR/RXR activation was found in AT of HEF cows. Relative mRNA expression of RBM39, which is involved in energy metabolism, was decreased in AT of HEF versus LEF. The relationship found between the AT proteins and/or metabolic pathways and the feed efficiency demonstrates that AT may reflect metabolic adaptations to high efficiency, and suggests that these proteins together with their metabolic mechanisms are suitable candidates as biomarkers to identify efficient cows for dairy production.

Effects of Environmental Heat Load on Endocannabinoid System Components in Adipose Tissue of High Yielding Dairy Cows

Simple Summary

We hypothesized that environmental heat load (HL) may affect the endocannabinoid system (ECS), a central regulator of metabolism and the stress response, in adipose tissue (AT), plasma and milk of dairy cows. In AT of summer vs. winter calving cows, gene expression of ECS components was decreased, but this was not translated to differences in protein abundance or in levels of endocannabinoids. In late-lactation cows that were not cooled vs. cooled, AT protein abundance of the heat sensitive, and ECS receptor, transient-receptor-potential-cation-channel-subfamily-V-member-1 (TRPV1) tended to be lower, and milk levels of 2-arachidonoylglycerol (2-AG) tended to increase in cows that were not cooled; but other ECS components were not different between groups. This suggests that HL is associated with limited alterations in the ECS of AT in dairy cows, either directly or via reduced feed intake.

Abstract

Environmental heat load (HL) adversely affects the performance of dairy cows. The endocannabinoid system (ECS) regulates metabolism and the stress response, thus we hypothesized that HL may affect the ECS of dairy cows. Our objective was to determine the levels of endocannabinoids (eCBs) and gene and protein expressions of the ECS components in adipose tissue (AT) and plasma of early postpartum (PP) and late-lactation cows. In addition, we examined eCBs in milk, and studied the interaction of eCBs with bovine cannabinoids receptors CB1 and CB2. In the first experiment, plasma and AT were sampled from cows calving during summer (S, n = 9) or winter (W, n = 9). Dry matter intake (DMI) and energy balance (EB) were lower in S vs. W, and relative gene expressions of transient-receptor-potential-cation-channel-subfamily-V-member-1 (TRPV1), the cannabinoid receptors CNR1 (CB1) and CNR2 (CB2), and monoglyceride lipase (MGLL) were decreased in AT of S compared to W. Protein abundance of peroxisome proliferator-activated-receptor-alpha (PPAR-α) was decreased, while tumor-necrosis factor-α (TNF-α) was increased in AT of S vs. W. Other components of the ECS were not different between S and W calving cows. To study whether the degree of HL may affect the ECS, we performed a second experiment with 24 late-lactation cows that were either cooled (CL) or not cooled (heat-stressed; HS) during summer. DMI was lower in HS vs. CL, AT protein abundance of PPAR-α was lower, and TRPV1 tended to be lower in HS vs. CL, but other components of the ECS were not different between groups. Milk levels of 2-arachidonoylglycerol (2-AG) tended to increase in HS vs. CL. Additionally, modeling of the bovine cannabinoid receptors demonstrated their binding to anandamide and 2-AG. Environmental HL, possibly via lower intake, is associated with limited alterations in ECS components in AT of dairy cows.

Antiviral peptides against the main protease of SARS-CoV-2: A molecular docking and dynamics study

The recent coronavirus outbreak has changed the world's economy and health sectors due to the high mortality and transmission rates. Because the development of new effective vaccines or treatments against the virus can take time, an urgent need exists for the rapid development and design of new drug candidates to combat this pathogen. Here, we obtained antiviral peptides obtained from the data repository of antimicrobial peptides (DRAMP) and screened their predicted tertiary structures for the ability to inhibit the main protease of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) using multiple combinatorial docking programs, including PatchDock, FireDock, and ClusPro. The four best peptides, DRAMP00877, DRAMP02333, DRAMP02669, and DRAMP03804, had binding energies of -1125.3, -1084.5, -1005.2, and -924.2 Kcal/mol, respectively, as determined using ClusPro, and binding energies of -55.37, -50.96, -49.25, -54.81 Kcal/mol, respectively, as determined using FireDock, which were better binding energy values than observed for other peptide molecules. These peptides were found to bind with the active cavity of the SARS-CoV-2 main protease; at Glu166, Cys145, Asn142, Phe140, and Met165, in addition to the substrate-binding sites, Domain 2 and Domain 3, whereas fewer interactions were observed with Domain 1. The docking studies were further confirmed by a molecular dynamics simulation study, in which several descriptors, including the root-mean-square difference (RMSD), root-mean-square fluctuation (RMSF), solvent-accessible surface area (SASA), radius of gyration (Rg), and hydrogen bond formation, confirmed the stable nature of the peptide-main protease complexes. Toxicity and allergenicity studies confirmed the non-allergenic nature of the peptides. This present study suggests that these identified antiviral peptide molecules might inhibit the main protease of SARS-CoV-2, although further wet-lab experiments remain necessary to verify these findings.

Antioxidant Resveratrol Increases Lipolytic and Reduces Lipogenic Gene Expression under In Vitro Heat Stress Conditions in Dedifferentiated Adipocyte-Derived Progeny Cells from Dairy Cows

Heat stress (HS) induces oxidative stress by increasing reactive oxygen species (ROS), and the polyphenol resveratrol (RSV) has been shown to have antioxidant properties by reducing ROS. Hence, we aimed to examine the effects of RSV, HS and their interaction on bovine adipocytes. We generated bovine dedifferentiated adipocyte-derived progeny (DFAT) cells from subcutaneous adipose tissue and examined the effects of RSV (100 µM), heat conditions: isothermal (ISO-37 °C), short heat (SH-41.2 °C for 1 h) and long HS (LH-41.2 °C for 16 h), and their interaction on gene expression in DFAT-cells. In medium of DFAT-cells treated with RSV, malondialdehyde levels were reduced and oxygen-radical absorbance-capacity levels were increased compared to control. Treating DFAT-cells with RSV increased the relative mRNA expression of stress-induced-phosphoprotein-1 (STIP1) and the expression of hormone-sensitive-lipase (LIPE) and perilipin-1 (PLIN1), whereas it reduced the expressions of fatty-acid-synthase (FASN) and of pro-inflammatory chemotactic-C-C-motif-chemokine-ligand-2 (CCL2) also under HS. Moreover, reduced protein abundance of FASN was found in RSV-treated DFAT-cells compared to controls. Molecular docking of RSV with FASN confirmed its possible binding to FASN active site. This work demonstrates that RSV has an antioxidant effect on bovine DFAT cells and may induce adipose lipolysis and reduce lipogenesis also under in vitro HS conditions.

Identification of Salicylic Acid Mechanism against Leaf Blight Disease in Oryza sativa by SR-FTIR Microspectroscopic and Docking Studies

The present study was to investigate the application and mechanism of salicylic acid (SA) as SA-Ricemate for the control of leaf blight disease using a Synchrotron Radiation-based Fourier-Transform Infra-Red (SR-FTIR) microspectroscopy and docking studies. After treating rice plants cv. KDML 105 with SA-Ricemate, the leaves were inoculated with Xanthomonas oryzae pv. oryzae, the causal agent of leaf blight, and disease severity were assessed. The leaves were also used to detect changes in endogenous SA content. The results indicated that SA-Ricemate, as an activated compound, reduced disease severity by 60% at three weeks post-inoculation and increased endogenous content by 50%. The SR-FTIR analysis of changes in the mesophyll of leaves (treated and untreated) showed that the groups of lipids, pectins, and proteins amide I and amide II occurred at higher values, and polysaccharides were shown at lower values in treated compared to untreated. Besides, docking studies were used to model a three-dimensional structure for Pathogenesis-related (PR1b) protein and further identify its interaction with SA. The results showed that ASP28, ARG31, LEU32, GLN97, and ALA93 are important residues that have strong hydrogen bonds with SA. The docking results showed that SA has a good interaction, confirming its role in expression.

Molecular Modeling of Histamine Receptors-Recent Advances in Drug Discovery

The recent developments of fast reliable docking, virtual screening and other algorithms gave rise to discovery of many novel ligands of histamine receptors that could be used for treatment of allergic inflammatory disorders, central nervous system pathologies, pain, cancer and obesity. Furthermore, the pharmacological profiles of ligands clearly indicate that these receptors may be considered as targets not only for selective but also for multi-target drugs that could be used for treatment of complex disorders such as Alzheimer's disease. Therefore, analysis of protein-ligand recognition in the binding site of histamine receptors and also other molecular targets has become a valuable tool in drug design toolkit. This review covers the period 2014-2020 in the field of theoretical investigations of histamine receptors mostly based on molecular modeling as well as the experimental characterization of novel ligands of these receptors.

Beneficial effects of minocycline on the ovary of polycystic ovary syndrome mouse model: Molecular docking analysis and evaluation of TNF-α, TNFR2, TLR-4 gene expression

Polycystic ovary syndrome (PCOS) is the most common cause of ovulatory infertility. Inflammation may be involved in the pathogenesis and development of PCOS. We investigated the anti-inflammatory effect of minocycline on TNF-α, TNFR2, and TLR4 expression levels and the key features of PCOS in a mouse model. Molecular docking was performed by Molecular Operating Environment software. PCOS was induced by estradiol valerate injection (EV) (2 mg/kg/day) in 40 mice. After 28 days, the mice were divided into five groups, including control, PCOS, minocycline control, minocycline PCOS model (50 mg/kg), and letrozole PCOS (0.5 mg/kg). The Levels of FSH, LH, E2, and testosterone were determined by ELISA. H&E staining was used for histological analysis in the ovarian tissues. Docking scores were -10.35, -10.57, and -12.45 kcal/mol for TNFα, TLR-4, and TNFR2, respectively. The expression levels of TNF-α, TNFR2, and TLR4 were detected by Real-Time PCR. PCOS models exhibited acyclicity, a significant increase in E2 levels (P < 0.01), and no difference in FSH, LH, and testosterone. The expression levels of TNF-α, TNFR2, and TLR-4 significantly increased in PCOS (2.70, 7.90, and 14.83-fold, respectively). EV treatment significantly increased graafian follicles (P < 0.001) and decreased corpus luteum (CL) (P < 0.01). Minocycline treatment in PCOS led to a significant decrease in E2 (P < 0.01) and graafian follicles (P < 0.001) and a significant increase in the CL numbers (P < 0.05). Our findings showed the positive effects of minocycline on estradiol level, CL and graafian follicles counts, suggesting that minocycline might inhibit these proteins and improve ovulation in our mouse model of PCOS.

Structure prediction, molecular simulations of RmlD from Mycobacterium tuberculosis, and interaction studies of Rhodanine derivatives for anti-tuberculosis activity

Tuberculosis is the most dangerous disease causing maximum deaths than any other, caused by single infectious agent. Due to multidrug resistant of Mycobacterium tuberculosis strains, there is a need of new drugs and drug targets. In this work, we have selected RmlD (α-dTDP-6-deoxy-lyxo-4-hexulose reductase) in the dTDP Rhamnose pathway as drug target to control tuberculosis using Rhodanine analogues. In order to study interaction of RmlD with Rhodanine analogues, a three-dimensional model based on crystal structures such as 1VLO from Clostridium, 1KBZ from Salmonella typhimurium, and 2GGS from Sulfolobus was generated using Modeller 9v7. The modeled structure reliability has been checked using programs such as Procheck, What if, Prosa, Verify 3D, and Errat. In an attempt to find new inhibitors for RmlD enzyme, docking studies were done with a series of Rhodanine and its analogues. Detailed analysis of enzyme-inhibitor interactions identified specific key residues, SER5, VAL9, ILE51, HIS54, and GLY55 which were important in forming hydrogen bonds in binding affinity. Homology modeling and docking studies on RmlD model provided valuable insight information for designing better inhibitors as novel anti-tuberculosis drugs by rational method.

Histamine H1 receptor antagonists selectively kill cisplatin-resistant human cancer cells

Cancer therapy is often hampered by the disease's development of resistance to anticancer drugs. We previously showed that the autonomously upregulated product of fibroblast growth factor 13 gene (FGF13; also known as FGF homologous factor 2 (FHF2)) is responsible for the cisplatin resistance of HeLa cisR cells and that it is likely responsible for the poor prognosis of cervical cancer patients treated with cisplatin. Here we show that cloperastine and two other histamine H1 receptor antagonists selectively kill HeLa cisR cells at concentrations that little affect parental HeLa S cells. The sensitivity of HeLa cisR cells to cloperastine was abolished by knocking down FGF13 expression. Cisplatin-resistant A549 cisR cells were similarly susceptible to cloperastine. H2, H3, and H4 receptor antagonists showed less or no cytotoxicity toward HeLa cisR or A549 cisR cells. These results indicate that histamine H1 receptor antagonists selectively kill cisplatin-resistant human cancer cells and suggest that this effect is exerted through a molecular mechanism involving autocrine histamine activity and high-level expression of FGF13. We think this represents a potential opportunity to utilize H1 receptor antagonists in combination with anticancer agents to treat cancers in which emergent drug-resistance is preventing effective treatment.

Exploration of stilbenoid trimers as potential inhibitors of sirtuin1 enzyme using a molecular docking and molecular dynamics simulation approach

A combination of molecular docking and molecular dynamics simulation (250 ns) has been carried out to study the interaction of stilbenoid trimer compounds with the SIRT1 enzyme as the target protein. SIRT1 expression regulates cellular stress responses that lead to the development of cancer. Redocking showed a good native ligand pose with an RMSD value of 1.40 Å at the receptor active site's coordinates. The molecular docking score uses a grid score functional (kcal mol⁻¹), which shows results of 1NS: 79.56, TS1: −26.83, TS2: −87.77, and TS3: −83.67. The TS2 and TS3 candidates were chosen for further analysis because they had a lower grid score than the native ligand (1NS). Furthermore, prediction of binding free energy (kcal mol⁻¹) using the Quantum Mechanics/generalized Born Surface Area (QM/MM-GBSA) method shows the results of 1NS: −31.52 ± 0.39, TS2: −58.99 ± 0.34, and TS3: −43.38 ± 0.35. These results indicate that the TS2 and TS3 compounds have good potential as inhibitors of the SIRT1 enzyme. Additionally, the amino acid residues were responsible for the inhibition mechanism through hydrogen bond interactions at the molecular level, including ASP22, PHE91, PRO11, ILE165, ASP166, and VAL230. The observations made in this study provide theoretical information for exploring the stilbenoid trimers as anticancer agents by targeting the SIRT1 enzyme.

A combination of molecular docking and molecular dynamics simulation (250 ns) has been carried out to study the interaction of stilbenoid trimer compounds with the SIRT1 enzyme as the target protein.

Molecular Modeling of µ Opioid Receptor Ligands with Various Functional Properties: PZM21, SR-17018, Morphine, and Fentanyl-Simulated Interaction Patterns Confronted with Experimental Data

Molecular modeling approaches are an indispensable part of the drug design process. They not only support the process of searching for new ligands of a given receptor, but they also play an important role in explaining particular activity pathways of a compound. In this study, a comprehensive molecular modeling protocol was developed to explain the observed activity profiles of selected µ opioid receptor agents: two G protein-biased µ opioid receptor agonists(PZM21 and SR-17018), unbiased morphine, and the β-arrestin-2-biased agonist,fentanyl. The study involved docking and molecular dynamics simulations carried out for three crystal structures of the target at a microsecond scale, followed by the statistical analysis of ligand-protein contacts. The interaction frequency between the modeled compounds and the subsequent residues of a protein during the simulation was also correlated with the output of in vitro and in vivo tests, resulting in the set of amino acids with the highest Pearson correlation coefficient values. Such indicated positions may serve as a guide for designing new G protein-biased ligands of the µ opioid receptor.