Interactions of 21 fentanyl derivatives with μ-opioid receptor (μOR) were studied using experimental and theoretical methods. Their binding to μOR was assessed with radioligand competitive binding assay. A uniform set of binding affinity data contains values for two novel and one previously uncharacterized derivative. The data confirms trends known so far and thanks to their uniformity, they facilitate further comparisons. In order to provide structural hypotheses explaining the experimental affinities, the complexes of the studied derivatives with μOR were modeled and subject to molecular dynamics simulations. Five common General Features (GFs) of fentanyls’ binding modes stemmed from these simulations. They include: GF1) the ionic interaction between D147 and the ligands’ piperidine NH⁺ moiety; GF2) the N-chain orientation towards the μOR interior; GF3) the other pole of ligands is directed towards the receptor outlet; GF4) the aromatic anilide ring penetrates the subpocket formed by TM3, TM4, ECL1 and ECL2; GF5) the 4-axial substituent (if present) is directed towards W318. Except for the ionic interaction with D147, the majority of fentanyl-μOR contacts is hydrophobic. Interestingly, it was possible to find nonlinear relationships between the binding affinity and the volume of the N-chain and/or anilide’s aromatic ring. This kind of relationships is consistent with the apolar character of interactions involved in ligand–receptor binding. The affinity reaches the optimum for medium size while it decreases for both large and small substituents. Additionally, a linear correlation between the volumes and the average dihedral angles of W293 and W133 was revealed by the molecular dynamics study. This seems particularly important, as the W293 residue is involved in the activation processes. Further, the Y326 (OH) and D147 (Cγ) distance found in the simulations also depends on the ligands’ size. In contrast, neither RMSF measures nor D114/Y336 hydrations show significant structure-based correlations. They also do not differentiate studied fentanyl derivatives. Eventually, none of 14 popular scoring functions yielded a significant correlation between the predicted and observed affinity data (R < 0.30, n = 28).

Antidepressants target a variety of proteins in the central nervous system (CNS), the most important belonging to the family of G-protein coupled receptors and the family of neurotransmitter transporters. The increasing number of crystallographic structures of these proteins have significantly contributed to the knowledge of their mechanism of action, as well as to the design of new drugs. Several computational approaches such as molecular docking, molecular dynamics, and virtual screening are useful for elucidating the mechanism of drug action and are important for drug design. This review is a survey of molecular targets for antidepressants in the CNS and computer based strategies to discover novel compounds with antidepressant activity.

Systemic lupus erythematosus (SLE) is a chronic and systemic autoimmune disease. SLE is described by production of autoantibodies and causes damage of many organs. T-cells play a crucial role in SLE pathogenesis. T-cells intensify inflammation through a number of processes, which leads to autoimmunization. CCR5 and MECP2 genes are linked with T-cells and pathogenesis of SLE. Polymorphisms in these genes are related with the prognostic factors of risk of disease onset and disease severity. The aim of this study was to estimate the influence of polymorphisms in MECP2 and CCR5 genes on the development and course of systemic lupus erythematosus. We examined 137 SLE patients and 604 healthy controls. We studied polymorphisms for CCR5 gene: rs333 and for MECP2: rs2075596, rs1734787, rs17435, and rs2239464. We genotyped our MECP2 samples and we performed a restriction fragment length polymorphism (RFLP) analysis for CCR5 samples. We showed a risk factor for allele T in rs17435 and for allele A in rs2075596 in MECP2. We noticed that MECP2 rs2075596 G/A, rs1734787 C/A, rs17435 A/T, and rs2239464 G/A polymorphisms are more prevalent in SLE patients than in healthy controls. We believe that above-mentioned MECP2 polymorphisms can be considered as SLE susceptibility factor.

Cancers are the leading cause of deaths worldwide. In 2018, an estimated 18.1 million new cancer cases and 9.6 million cancer-related deaths occurred globally. Several previous studies have shown that the enzyme, leucine aminopeptidase is involved in pathological conditions such as cancer. On the basis of the knowledge that isoquinoline alkaloids have antiproliferative activity and inhibitory activity towards leucine aminopeptidase, the present study was conducted a study which involved database search, virtual screening, and design of new potential leucine aminopeptidase inhibitors with a scaffold based on 3,4-dihydroisoquinoline. These compounds were then filtered through Lipinski’s “rule of five,” and 25 081 of them were then subjected to molecular docking. Next, three-dimensional quantitative structure-activity relationship (3D-QSAR) study was performed for the selected group of compounds with the best binding score results. The developed model, calculated by leave-one-out method, showed acceptable predictive and descriptive capability as represented by standard statistical parameters r² (0.997) and q² (0.717). Further, 35 compounds were identified to have an excellent predictive reliability. Finally, nine selected compounds were evaluated for drug-likeness and different pharmacokinetics parameters such as absorption, distribution, metabolism, excretion, and toxicity. Our methodology suggested that compounds with 3,4-dihydroisoquinoline moiety were potentially active in inhibiting leucine aminopeptidase and could be used for further in-depth in vitro and in vivo studies.

The structure, energetics, and aromaticity of c.a. 100 constitutional isomers and tautomers of pyrido[m,n]diazepines (m = 1, 2; n = 2, 3, 4, 5; m ≠ n) were studied at the B3LYP/cc-pVTZ level. The pyrido[1,3]diazepines appear the most, while pyrido[2,4]diazepines are the least stable (ca. 26 kcal/mol). In the pyrido[1,n]diazepine group (n = 2–5), the [1,5] isomers are higher in energy by ca. 4.5 kcal/mol and the [1,4] ones by ca. 7 kcal/mol, and the pyrido[1,2]diazepines are the least stable (ca. 20 kcal/mol). All the most stable pyrido[1,n]diazepines have N-atoms near the ring’s junction bond but on opposite sites. The most stable [2,n]-forms are also those with the pyridine ring N6-atom near the junction bond. Surprisingly, for the [1,2]-, [1,3]-, and [1,4]-isomer condensation types of pyridine and diazepine rings, the same N9 > N7 > N6 > N8 stability pattern obeys. The stability remains similar in a water medium simulated with the Polarizable Continuum Model of the solvent and is conserved when calculated using the CAM-B3LYP or BHandHlyp functionals. The ring’s aromaticity in the pyridine[m,n]diazepines was established based on the integral INICS index resulting from the NICSzz-scan curves’ integration. The integral INICS index is physically justified through its relation to the ringcurrent as demonstrated by Berger, R.J.F., et al. Phys. Chem. Chem. Phys. 2022, 24, 624. The six-membered pyrido rings have negative INICS_ZZ indices and can be aromatic only if they are not protonated at the N-atom. All protonated pyrido and seven-membered rings exhibit meaningful positive INICS_ZZ values and can be assigned as antiaromatic. However, some non-protonated pyrido rings also have substantial positive INICS_ZZ indices and are antiaromatic. A weak linear correlation (R² = 0.72) between the INICS_ZZ values of the pyridine I(6) and diazepine I(7) rings exists and is a consequence of the communication between the π-electron systems of the two rings. The juxtaposition of the INICS descriptor of the six- and seven-membered rings and diverse electron density parameters at the Ring Critical Points (RCP) revealed good correlations only with the Electrostatic Potentials from the electrons and nuclei (ESPe and ESPn). The relationships with other RCP parameters like electron density and its Laplacian, total energy, and the Hamiltonian form of kinetic energy density were split into two parts: one nearly constant for the six-membered rings and one linearly correlating for the seven-membered rings. Thus, most of the electron density parameters at the RCP of the six-membered rings of pyridodiazepines practically do not change with the diazepine type and the labile proton position. In contrast, those of the seven-membered rings display aromaticity changes in the antiaromatic diazepine with its ring structural modifications.

The cytotoxic activity of several serotonin transporter (SERT) inhibitors and subtype of serotonin receptor 1A (5-HT1A receptor) ligands have been examined in androgen-insensitive human PC-3 prostate and neuroblastoma SH-SY5Y cancer cells. Almost all of the studied compounds (except 5-HT1A receptor agonist (2R)-(+)-8-Hydroxy-2-(di-n-propylamino)tetralin hydrobromide (8-OH-DPAT)) exhibited absolute cytotoxic activity against the examined cancer cells. The compound 4-Fluoro-N-[2-[4-(7-methoxy-1-naphthalenyl)-1-piperazinyl]ethyl]benzamide hydrochloride (S14506) that showed highest activity against neuroblastoma tumors was the 5-HT1A receptor agonist (although not alike other 5-HT1A receptor agonists). On the other hand, the compound 6-nitro-2-(4-undecylpiperazin-1-yl)quinoline hydrochloride (AZ07) that had the highest activity against PC-3 prostate cancer cells was a compound exhibiting antagonistic activity against the 5-HT1A receptor. Thus, compounds of oncotoxic properties S14506 and AZ07 should be evaluated further for their potential use in the prevention and treatment of cancer. Most of the 15 compounds tested exhibited either agonistic or antagonistic activity for both the cyclic adenosine monophosphate (cAMP) and extracellular signal-regulated kinase 1 and 2 (ERK1/2) pathways in human embryonic kidney 293 (HEK293) cells that overexpress the 5HT1AR gene. However, compounds paroxetine, N-Ac-paroxetine and 2-[4-(cyclobutylmethyl)piperazin-1-yl]-6-nitroquinoline hydrochloride (AB22) simultaneously exhibited antagonistic activity on the cAMP pathway and agonistic activity on the ERK1/2 pathway. Fluoxetine relative to compound AZ07 had almost three times lower cytotoxic activity against PC-3 prostate cancer cells. However, the proapoptotic activity of fluoxetine compared to compound AZ07 is almost two times higher which would suggest that the cytotoxic activity of both compounds may be dependent on different cell death mechanisms. Compound S14506 was found to be an antagonist of the serine-threonine protein kinase B (Akt) pathway. Prosurvival Akt activity may be reversed by Akt antagonists. Therefore, the antagonistic activity of S14506 on the Akt pathway may evoke caspase-3 expression and cytotoxicity. It appears that one should not expect a straightforward relationship between the activation of particular serotonergic pathways by selective serotonin reuptake inhibitors (SSRIs) and 5-HT1A receptor ligands and their cytotoxic or cytoprotective activity. Additionally, nuclear transcription factor κB (NF-κB), which may be involved in 5-HT-dependent biochemical pathways by coordinating different subunits in the formation of a dimer, may regulate the transcription of different transduction pathways. Therefore, it can be suggested that the mechanism of the cytotoxic activity of certain compounds (serotonergic against nonserotonergic) may depend on the compound and cancer type being examined. Docking studies showed that S14506, buspirone and spiperone bind in similar ways in the 5-HT1A receptor model and interacted with similar 5-HT1A receptor residues. S14506 and spiperone were found to be located closer to both phenylalanines in TM6 than buspirone, thus exhibiting more antagonist binding modes.