Exploring the interaction of N-(benzothiazol-2-yl)pyrrolamide DNA gyrase inhibitors with the GyrB ATP-binding site lipophilic floor: A medicinal chemistry and QTAIM study.

N-(Benzothiazole-2-yl)pyrrolamide DNA gyrase inhibitors with benzyl or phenethyl substituents attached to position 3 of the benzothiazole ring or to the carboxamide nitrogen atom were prepared and studied for their inhibition of Escherichia coli DNA gyrase by supercoiling assay. Compared to inhibitors bearing the substituents at position 4 of the benzothiazole ring, the inhibition was attenuated by moving the substituent to position 3 and further to the carboxamide nitrogen atom. A co-crystal structure of (Z)-3-benzyl-2-((4,5-dibromo-1H-pyrrole-2-carbonyl)imino)-2,3-dihydrobenzo[d]-thiazole-6-carboxylic acid (I) in complex with E. coli GyrB24 (ATPase subdomain) was solved, revealing the binding mode of this type of inhibitor to the ATP-binding pocket of the E. coli GyrB subunit. The key binding interactions were identified and their contribution to binding was rationalised by quantum theory of atoms in molecules (QTAIM) analysis. Our study shows that the benzyl or phenethyl substituents bound to the benzothiazole core interact with the lipophilic floor of the active site, which consists mainly of residues Gly101, Gly102, Lys103 and Ser108. Compounds with substituents at position 3 of the benzothiazole core were up to two orders of magnitude more effective than compounds with substituents at the carboxamide nitrogen. In addition, the 6-oxalylamino compounds were more potent inhibitors of E. coli DNA gyrase than the corresponding 6-acetamido analogues.

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 D2 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).

Combined MD/QTAIM techniques to evaluate ligand-receptor interactions. Scope and limitations.

In medicinal chemistry, it is extremely important to evaluate, as accurately as possible, the molecular interactions involved in the formation of different ligand-receptor (L-R) complexes. Evaluating the different molecular interactions by quantum mechanics calculations is not a simple task, since formation of an L-R complex is a dynamic process. In this case, the use of combined techniques of molecular dynamics (MD) and quantum calculations is one the best possible approaches. In this work we report a comparative study using combined MD and QTAIM (Quantum Theory of Atoms In Molecules) calculations for five biological systems with different levels of structural complexity. We have studied Acetylcholinesterase (AChE), D2 Dopamine Receptor (D2DR), beta Secretase (BACE1), Dihydrofolate Reductase (DHFR) and Sphingosine Kinase 1 (SphK1). In these molecular targets, we have analyzed different ligands with diverse structural characteristics. The inhibitory activities of most of them have been previously measured in our laboratory. Our results indicate that QTAIM calculations can be extremely useful for in silico studies. It is possible to obtain very accurate information about the strength of the molecular interactions that stabilize the formation of the different L-R complexes. Better correlations can be obtained between theoretical and experimental data by using QTAIM calculations, allowing us to discriminate among ligands with similar affinities. QTAIM analysis gives fairly accurate information for weak interactions which are not well described by MD simulations. QTAIM study also allowed us to evaluate and determine which parts of the ligand need to be modified in order to increase its interactions with the molecular target. In this study we have discussed the importance of combined MD/QTAIM calculations for this type of simulations, showing their scopes and limitations.

Conformational and electronic study of dopamine interacting with the D2 dopamine receptor.

We report an exhaustive conformational and electronic study on dopamine (DA) interacting with the D2 dopamine receptor (D2 DR). For the first time, the complete surface of the conformational potential energy of the complex DA/D2 DR is reported. Such a surface was obtained through the use of QM/MM calculations. A detailed study of the molecular interactions that stabilize and destabilize the different molecular complexes was carried out using two techniques: Quantum Theory of Atoms in Molecules computations and nuclear magnetic shielding constants calculations. A comparative study of the behavior of DA in the gas phase, aqueous solution, and in the active site of D2 DR has allowed us to evaluate the degree of deformation suffered by the ligand and, therefore, analyze how rustic are the lock-key model and the induced fit theory in this case. Our results allow us to propose one of the conformations obtained as the "biologically relevant" conformation of DA when it is interacting with the D2 DR.

Antinociceptive effect of neo-clerodane diterpenes obtained from Baccharis flabellata.

We report here for the first time antinociceptive effects of extracts from Baccharis flabellata. Two extracts in this analysis, one obtained in summer and the other during winter time. Our results indicate that both extract show strong antinociceptive effects, being the extracts obtained during the summer significantly more active. Our results suggest that this activity is mainly due to the presence of the diene-acid clerodane ent-15,16-epoxy-19-hydroxy-1,3,13(16),14-clerodatetraen-18-oic acid (DAC) and its dimer called DACD. Employing naloxone as an antagonist of opioid receptors, we demonstrated that both compounds act on opioid receptors, being the antinociceptive effect of DACD stronger than DAC. Thus, the antinociceptive activity of DACD was almost two times stronger than DAC (44.8 over 24.6 s in the hot-plate test) after one hour of treatments. In order to better understand the mechanism of action at molecular level of these compounds, we conducted a molecular modeling study analyzing the molecular interactions of DAC and DACD complexes with the κ-ORs. Our results suggest interactions for both DAC and DACD with Gln115, Val118, Tyr119, Asn122 and Tyr313 stabilizing their complexes; however, these interactions are significantly stronger for DACD with respect to DAC. This finding could explain why DACD have a higher affinity for the κ-ORs. These results are in agreement with the obtained antinociceptive effect. In addition, our results indicate that these neoclerodanes would have a mechanism of action similar to that of salvinorin A; such information can be very useful for the design of new inhibitors of κ-ORs.

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 IC50 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.