Antinociceptive effect of cyclic and linear diterpenoids as new atypical agonists of κ-opioid receptors obtained from four species of the Baccharis genus, and vehiculated in nanometric niosomes

New natural analgesic compounds that act in KORs are important alternatives for potential therapeutical use in medicine. In this work, we report and compare here the antinociceptive activity displayed by cyclic and linear diterpenes, obtained from the genus Baccharis. The antinociceptive activities determined were relatively strong, in comparison whit morphine. The antinociceptive mechanism of action was made through naloxone administration (a non-selective antagonist of opioid receptors). The more active compounds were vehiculized successfully in niosomes at nanometric scale. The observed antinociceptive activity for Bartemidiolide oxide (BARTO), obtain from Baccharis artemisioides, was greater than Flabeloic acid dimer (DACD), the first compound isolated from Baccharis flabellata that was reported possessing antinociceptive effects. We also conducted docking calculations and molecular dynamics simulations, which suggested that the newly identified diterpenes might share the molecular action mechanism reported for Salvinorin A (SalA). Molecular simulations have allowed us to appreciate some subtle differences between molecular interactions of these ligands stabilizing their respective complexes; such information might be useful for designing and searching for new inhibitors of KORs.

Evaluating the conformational space of the active site of D2 dopamine receptor. Scope and limitations of the standard docking methods

We report here for the first time the potential energy surfaces (PES) of phenyletilamine (PEA) and meta-tyramine (m-OH-PEA) at the D₂ dopamine receptor (D2DR) binding site. PESs not only allow us to observe all the critical points of the surface (minimums, maximums, and transition states), but also to note the ease or difficulty that each local minima have for their conformational inter-conversions and therefore know the conformational flexibility that these ligands have in their active sites. Taking advantage of possessing this valuable information, we analyze how accurate a standard docking study is in these cases. Our results indicate that although we have to be careful in how to carry out this type of study and to consider performing some extra-simulations, docking calculations can be satisfactory. In order to analyze in detail the different molecular interactions that are stabilizing the different ligand-receptor (L-R) complexes, we carried out quantum theory of atoms in molecules (QTAIM) computations and NMR shielding calculations. Although some of these techniques are a bit tedious and require more computational time, our results demonstrate the importance of performing computational simulations using different types of combined techniques (docking/MD/hybrid QM-MM/QTAIM and NMR shielding calculations) in order to obtain more accurate results. Our results allow us to understand in details the molecular interactions stabilizing and destabilizing the different L-R complexes reported here. Thus, the different activities observed for dopamine (DA), m-OH-PEA, and PEA can be clearly explained at molecular level.

Quinazoline-tethered hydrazone: A versatile scaffold toward dual anti-TB and EGFR inhibition activities in NSCLC

Globally, lung cancer and tuberculosis are considered to be very serious and complex diseases. Evidence suggests that chronic infection with tuberculosis (TB) can often lead to lung tumors; therefore, developing drugs that target both diseases is of great clinical significance. In our study, we designed and synthesized a suite of 14 new quinazolinones (5a-n) and performed biological investigations of these compounds in Mycobacterium tuberculosis (MTB) and cancer cell lines. In addition, we conducted a molecular modeling study to determine the mechanism of action of these compounds at the molecular level. Compounds that showed anticancer activity in the preliminary screening were further evaluated in three cancer cell lines (A549, Calu-3, and BT-474 cells) and characterized in an epidermal growth factor receptor (EGFR) binding assay. Cytotoxicity in noncancerous lung fibroblast cells was also evaluated to obtain safety data. Our theoretical and experimental studies indicated that our compounds showed a mechanism of action similar to that of erlotinib by inhibiting the EGFR tyrosine kinase. In turn, the antituberculosis activity of these compounds would be produced by the inhibition of enoyl-ACP-reductase. From our findings, we were able to identify two potential lead compounds (5i and 5l) with dual activity and elevated safety toward noncancerous lung fibroblast cells. In addition, our data identified three compounds with excellent anti-TB activities (compounds 5i, 5l, and 5n).

Molecular modeling study of dihydrofolate reductase inhibitors. Molecular dynamics simulations, quantum mechanical calculations, and experimental corroboration

A molecular modeling study on dihydrofolate reductase (DHFR) inhibitors was carried out. By combining molecular dynamics simulations with semiempirical (PM6), ab initio, and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of DHFR inhibitors interacting with the human enzyme is reported here, providing a clear picture of the binding interactions of these ligands from both structural and energetic viewpoints. A reduced model for the binding pocket was used. This approach allows us to perform more accurate quantum mechanical calculations as well as to obtain a detailed electronic analysis using the quantum theory of atoms in molecules (QTAIM) technique. Thus, molecular aspects of the binding interactions between inhibitors and the DHFR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental IC₅₀ values was obtained, predicting with an acceptable qualitative accuracy the potential inhibitor effect of nonsynthesized compounds. Such correlation was experimentally corroborated synthesizing and testing two new inhibitors reported in this paper.

Searching the "biologically relevant"conformation of dopamine: a computational approach

We report here an exhaustive and complete conformational study on the conformational potential energy hypersurface (PEHS) of dopamine (DA) interacting with the dopamine D2 receptor (D2-DR). A reduced 3D model for the binding pocket of the human D2-DR was constructed on the basis of the theoretical model structure of bacteriorhodopsin. In our reduced model system, only 13 amino acids were included to perform the quantum mechanics calculations. To obtain the different complexes of DA/D2-DR, we combined semiempirical (PM6), DFT (B3LYP/6-31G(d)), and QTAIM calculations. The molecular flexibility of DA interacting with the D2-DR was evaluated from potential energy surfaces and potential energy curves. A comparative study between the molecular flexibility of DA in the gas phase and at D2-DR was carried out. In addition, several molecular dynamics simulations were carried out to evaluate the molecular flexibility of the different complexes obtained. Our results allow us to postulate the complexes of type A as the "biologically relevant conformations" of DA. In addition, the theoretical calculations reported here suggested that a mechanistic stepwise process takes place for DA in which the protonated nitrogen group (in any conformation) acts as the anchoring portion, and this process is followed by a rapid rearrangement of the conformation allowing the interaction of the catecholic OH groups.

Ab initio and DFT study of the conformational energy hypersurface of cyclic Gly-Gly-Gly

The multidimensional conformational potential energy hypersurface (PEHS) of cyclic Gly-Gly-Gly (1,4,7-triazonane-2,5,8-trione) was comprehensively investigated at the Hartree-Fock (RHF/6-31G(d)) level of theory. The equilibrium structures, their relative stability, and the transition state (TS) structures involved in the conformational interconversion pathways were analyzed. aug-cc-pVTZ//B3LYP/6-311++G** single point calculations predict a trans-cis-cis conformation as the energetically preferred form for this compound. However, all of the levels of theory employed here predicted that two forms, a trans-cis-cis and a cis-cis-cis (crown), of conformers contribute significantly to the equilibrium mixture at room temperature. The conformational interconversion between the global minimum and the symmetric cis-cis-cis crown form requires 12.49 kcal/mol at the RHF 6-31G(d) level of theory, whereas the conformational interconversion between the cis-cis-cis crown and cis-cis-cis boat form requires 18.70 kcal/mol. An exploratory topological analysis of the PEHS was also carried out. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of these cyclic tripeptides, describing the conformations as well as the conformational interconversion processes in these hypersurfaces.