Development of Computational Approaches with a Fragment-Based Drug Design Strategy: In Silico Hsp90 Inhibitors Discovery.

A semi-exhaustive approach and a heuristic search algorithm use a fragment-based drug design (FBDD) strategy for designing new inhibitors in an in silico process. A deconstruction reconstruction process uses a set of known Hsp90 ligands for generating new ones. The deconstruction process consists of cutting off a known ligand in fragments. The reconstruction process consists of coupling fragments to develop a new set of ligands. For evaluating the approaches, we compare the binding energy of the new ligands with the known ligands.

The catalytic effects of a thiazolium salt in the oxa-Diels-Alder reaction between benzaldehyde and Danishefsky's diene: a molecular electron density theory study.

The oxa-Diels-Alder (ODA) reaction of benzaldehyde with Danishefsky's diene in the presence of a [thiazolium][Cl] salt, as a model of an ionic liquid, has been studied within Molecular Electron Density Theory (MEDT) at the M06-2X/6-311G(d,p) computational level. The formation of two hydrogen bonds (HBs) between the thiazolium cation and the carbonyl oxygen of benzaldehyde modifies neither the electrophilic character of benzaldehyde nor its electronic structure substantially but accelerates the reaction considerably. This ODA reaction presents an activation energy of 4.5 kcal mol-1; the formation of the only observed dihydropyranone is strongly exothermic by -28.8 kcal mol-1. The presence of the [thiazolium][Cl] salt decreases the Gibbs free energy of activation of the ODA reaction between benzaldehyde and Danishefsky's diene by 5.9 kcal mol-1. This ODA reaction presents total para regioselectivity and high endo stereoselectivity. This ODA reaction takes place through a highly asynchronous polar transition state structure (TS) associated with a non-concerted two-stage one-step mechanism. ELF analysis of para/endo TSs associated with the ODA reactions in the absence and presence of the [thiazolium][Cl] salt shows that the formation of the HBs at the TSs does not modify their electronic structure substantially. This MEDT study makes it possible to conclude that the acceleration found in the ODA reaction of benzaldehyde with Danishefsky's diene in ILs is a consequence of an increase of the global electron density transfer at TS3-pn, resulting from HB formation, and the greater strength of the HBs at the polar TS3-pn compared to that at the benzaldehyde : [thiazolium][Cl] complex, and that the strength in the HB formed is more relevant that than an increase of the electrophilic character of the interaction between reagent.

New benzimidazolequinones as trypanosomicidal agents.

Herein, the design and synthesis of new 2-phenyl(pyridinyl)benzimidazolequinones and their 5-phenoxy derivatives as potential anti-Trypanosoma cruzi agents are described. The compounds were evaluated in vitro against the epimastigotes and trypomastigote forms of Trypanosoma cruzi. The replacing of a benzene moiety in the naphthoquinone system by an imidazole enhanced the trypanosomicidal activity against Trypanosoma cruzi. Three of the tested compounds (11a-c) showed potent trypanosomicidal activity and compound 11a, with IC50 of 0.65 µM on the trypomastigote form of T. cruzi, proved to be 15 times more active than nifurtimox. Additionally, molecular docking studies indicate that the quinone derivatives 11a-c could have a multitarget profile interacting preferentially with trypanothione reductase and Old Yellow Enzyme.

Identification of the naphthoquinone derivative inhibitors binding site in heat shock protein 90: an induced-fit docking, molecular dynamics and 3D-QSAR study.

The combination of molecular modeling methods to identify the putative binding site of inhibitors constitutes an important tool in drug discovery. In this work, we used these analyses to understand the potent inhibitory effect of naphthoquinone derivatives on heat shock protein 90 (Hsp90), one of the proteins involved in many types of cancer. Molecular docking results indicated that some favorable interactions of key amino acid residues at the binding site of Hsp90 with these quinones would be responsible for the inhibition of Hsp90 activity. Molecular docking and molecular dynamics simulation were carried out to further understand the binding modes and the interactions between the protein and these inhibitors. The main residues of the internal cavity were Val136, Phe138, Tyr139, Val150, Trp162 and Val186. The high concordance between the docking results and 3D-QSAR contour maps gives us helpful information about the environment of the binding site. Our results provide the bases for a rational modification of new molecules based in quinone scaffold, in order to design more potent Hsp90 inhibitors, which would exhibit highly potent antitumor activity.Communicated by Ramaswamy H. Sarma.

Experimental and Theoretical Approaches in the Study of Phenanthroline-Tetrahydroquinolines for Alzheimer's Disease.

The imino-Diels-Alder reaction is one of the most common strategies in organic chemistry and is an important tool for providing a broad spectrum of biologically active heterocyclic systems. A combined theoretical and experimental study of the imino-Diels-Alder reaction is described. The new phenanthroline-tetrahydroquinolines were evaluated as cholinesterase inhibitors. Their cytotoxicity in human neuroblastoma SH-SY5Y cells was also evaluated. The theoretical results suggest that compounds formation in stages can be explained by endo cycloadducts under the established reaction conditions, thereby confirming experimental results obtained for percentage yield. These results allowed us to establish that pyridine substituent remarkably influences activation energy and reaction yield, as well as in acetylcholinesterase (AChE) activity. Among these derivatives, compounds with 4-pyridyl and 4-nitrophenyl showed favorable AChE activity and proved to be non-cytotoxic.

State of the art of Smo antagonists for cancer therapy: advances in the target receptor and new ligand structures.

Since the Hedgehog signaling pathway has been associated with cancer, it has emerged as a therapeutic target for cancer therapy. The main target among the key Hedgehog proteins is the GPCR-like Smo receptor. Therefore, some Smo antagonists that have entered clinical trials, including the US FDA-approved drugs vismodegib and sonidegib, to treat basal cell carcinoma and medulloblastoma. However, early resistance of these drugs has spawned the need to understand the molecular bases of this phenomena. We therefore reviewed details about Smo receptor structures and the best Smo antagonist chemical structures. In addition, we discussed strategies that should be considered to develop new, safer generations of Smo antagonists that avoid current clinical limitations.

Ion-Ion Repulsions and Charge-Shielding Effects Dominate the Permeation Mechanism through the OmpF Porin Channel.

OmpF is a wide channel bacterial porin frequently employed to study selective ionic translocation. The cationic preference of this porin is mainly determined by electrostatic forces between the translocated ion and the protein and the formation of ion pairs (e.g., K+···Cl-) being previously pointed as the main cause to favor the cationic transport through the constriction zone. Hybrid quantum mechanics/molecular mechanics-molecular dynamics simulations, which have provided polarization-containing potentials of mean force profiles for different permeation scenarios, reveal significant new insights related with the ion translocation mechanism. Results show that the permeation is dominated by electrostatic interactions, which in turn affect ion-protein interactions at the constriction zone. However, it is observed that ion flow is favored by ion-ion repulsions and, in a lesser extent, by charge-shielding effects, instead of the previously pointed ionic pair formation.

Structure-Activity Relationships Based on 3D-QSAR CoMFA/CoMSIA and Design of Aryloxypropanol-Amine Agonists with Selectivity for the Human ß3-Adrenergic Receptor and Anti-Obesity and Anti-Diabetic Profiles.

The wide tissue distribution of the adrenergic ß3 receptor makes it a potential target for the treatment of multiple pathologies such as diabetes, obesity, depression, overactive bladder (OAB), and cancer. Currently, there is only one drug on the market, mirabegron, approved for the treatment of OAB. In the present study, we have carried out an extensive structure-activity relationship analysis of a series of 41 aryloxypropanolamine compounds based on three-dimensional quantitative structure-activity relationship (3D-QSAR) techniques. This is the first combined comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) study in a series of selective aryloxypropanolamines displaying anti-diabetes and anti-obesity pharmacological profiles. The best CoMFA and CoMSIA models presented values of r²ncv = 0.993 and 0.984 and values of r²test = 0.865 and 0.918, respectively. The results obtained were subjected to extensive external validation (q², r², r²m, etc.) and a final series of compounds was designed and their biological activity was predicted (best pEC50 = 8.561).

Combined molecular modelling and 3D-QSAR study for understanding the inhibition of NQO1 by heterocyclic quinone derivatives.

A combination of three-dimensional quantitative structure-activity relationship (3D-QSAR), and molecular modelling methods were used to understand the potent inhibitory NAD(P)H:quinone oxidoreductase 1 (NQO1) activity of a set of 52 heterocyclic quinones. Molecular docking results indicated that some favourable interactions of key amino acid residues at the binding site of NQO1 with these quinones would be responsible for an improvement of the NQO1 activity of these compounds. The main interactions involved are hydrogen bond of the amino group of residue Tyr128, π-stacking interactions with Phe106 and Phe178, and electrostatic interactions with flavin adenine dinucleotide (FADH) cofactor. Three models were prepared by 3D-QSAR analysis. The models derived from Model I and Model III, shown leave-one-out cross-validation correlation coefficients (q2LOO ) of .75 and .73 as well as conventional correlation coefficients (R2 ) of .93 and .95, respectively. In addition, the external predictive abilities of these models were evaluated using a test set, producing the predicted correlation coefficients (r2pred ) of .76 and .74, respectively. The good concordance between the docking results and 3D-QSAR contour maps provides helpful information about a rational modification of new molecules based in quinone scaffold, in order to design more potent NQO1 inhibitors, which would exhibit highly potent antitumor activity.

2D-QSAR and 3D-QSAR/CoMSIA Studies on a Series of (R)-2-((2-(1H-Indol-2-yl)ethyl)amino)-1-Phenylethan-1-ol with Human ß3-Adrenergic Activity.

The ß3 adrenergic receptor is raising as an important drug target for the treatment of pathologies such as diabetes, obesity, depression, and cardiac diseases among others. Several attempts to obtain selective and high affinity ligands have been made. Currently, Mirabegron is the only available drug on the market that targets this receptor approved for the treatment of overactive bladder. However, the FDA (Food and Drug Administration) in USA and the MHRA (Medicines and Healthcare products Regulatory Agency) in UK have made reports of potentially life-threatening side effects associated with the administration of Mirabegron, casting doubts on the continuity of this compound. Therefore, it is of utmost importance to gather information for the rational design and synthesis of new ß3 adrenergic ligands. Herein, we present the first combined 2D-QSAR (two-dimensional Quantitative Structure-Activity Relationship) and 3D-QSAR/CoMSIA (three-dimensional Quantitative Structure-Activity Relationship/Comparative Molecular Similarity Index Analysis) study on a series of potent ß3 adrenergic agonists of indole-alkylamine structure. We found a series of changes that can be made in the steric, hydrogen-bond donor and acceptor, lipophilicity and molar refractivity properties of the compounds to generate new promising molecules. Finally, based on our analysis, a summary and a regiospecific description of the requirements for improving ß3 adrenergic activity is given.

Structural analysis and molecular docking of trypanocidal aryloxy-quinones in trypanothione and glutathione reductases: a comparison with biochemical data.

A set of aryloxy-quinones, previously synthesized and evaluated against Trypanosoma cruzi epimastigotes cultures, were found more potent and selective than nifurtimox. One of the possible mechanisms of the trypanocidal activity of these quinones could be inhibition of trypanothione reductase (TR). Considering that glutathione reductase (GR) is the equivalent of TR in humans, biochemical, kinetic, and molecular docking studies in TR and GR were envisaged and compared with the trypanocidal and cytotoxic data of a set of aryloxy-quinones. Biochemical assays indicated that three naphthoquinones (Nq-h, Nq-g, and Nq-d) selectively inhibit TR and the TR kinetic analyses indicated that Nq-h inhibit TR in a noncompetitive mechanism. Molecular dockings were performed in TR and GR in the following three putative binding sites: the catalytic site, the dimer interface, and the nicotinamide adenine dinucleotide phosphate-binding site. In TR and GR, the aryloxy-quinones were found to exhibit high affinity for a site near it cognate-binding site in a place in which the noncompetitive kinetics could be justified. Taking as examples the three compounds with TR specificity (TRS) (Nq-h, Nq-g, and Nq-d), the presence of a network of contacts with the quinonic ring sustained by the triad of Lys62, Met400', Ser464' residues, seems to contribute hardly to the TRS. Compound Nq-b, a naphthoquinone with nitrophenoxy substituent, proved to be the best scaffold for the design of trypanocidal compounds with low toxicity. However, the compound displayed only a poor and non-selective effect toward TR indicating that TR inhibition is not the main reason for the antiparasitic activity of the aryloxy-quinones.

Multiscale Treatment for the Molecular Mechanism of a Diels-Alder Reaction in Solution: A QM/MM-MD Study.

Thermodynamics and the solvent role in the acceleration of the Diels-Alder reaction between cyclopentadiene (CPD) and methyl vinyl ketone (MVK) have been revisited. In this work we use an ab initio hybrid QM/MM-MD scheme combined with multiple steered molecular dynamics to extract the free energy pofile in water and methanol using the bidirectional Minh-Adib estimator. We obtain 18.7 kcal mol-1 and 20.8 kcal mol-1 free energy barrier for the reaction in water and methanol, respectively. This methodology reproduces experimental values with an absolute error of about 0.8 kcal mol-1. The experimental difference between the activation free-energy barriers of water and methanol is also reproduced with an absolute error of about 0.1 kcal mol-1. We explore the charge transfer evolution along reaction coordinates to characterize the electronic behavior for this reaction. It is shown that the solvent molecules around the reaction system produce a global polarization along the reaction coordinate which is consistent with the solvent polarity. The results highlight the role of hydrogen bonding formed in the transition state to stabilize the system charge reorganization in the reaction process.

Synthesis and biological characterization of new aryloxyindole-4,9-diones as potent trypanosomicidal agents.

A new indole-4,9-dione and their phenoxy derivatives were synthesized and evaluated in vitro against the epimastigote form of Trypanosoma cruzi, Y strain. All of these novel compounds were found to be extremely potent and selective that the standard drug nifurtimox. Interestingly, phenoxyindole-4,9-dione 9d displayed excellent nanomolar inhibitory activity, IC50=20 nM, and high selectivity index, SI=625. In silico studies using MOE program were performed to generate a preliminary pharmacophore model.

Theoretical study on the molecular mechanism of the [5 + 2] vs. [4 + 2] cyclization mediated by Lewis acid in the quinone system.

The thermal and Lewis acid (LA) catalyzed cyclizations of quinone 1 involved in the synthesis of Colombiasin A and Elipsaterosin B have been theoretically studied using DFT methods at the B3LYP/6-311G(d,p) computational level. B3LYP calculations suggest that the formal endo [4 + 2] cycloadduct allowing the synthesis of Colombiasin A is preferentially formed under thermal conditions, while in the presence of the BF3 LA catalyst the formal [5 + 2] cycloadduct is seen, allowing the synthesis of Elipsaterosin B. The BF3 LA catalyst not only accelerates the nucleophilic attack on the C2 carbon of the quinone framework through a more polar C-C bond formation, but also provokes a different electron density rearrangement along the nucleophilic attack favoring the subsequent C-C bond formation at the C4 carbon with the formation of the formal [5 + 2] cycloadduct. ELF bonding analysis along these cyclizations indicates that the C-C single bond formation takes place in the range of 1.91-2.1 Å by C-to-C coupling of two pseudoradical centers. Along the formation of the first C2-C9 single bond, these pseudoradical centers appear at one of the most electrophilic and at one of the most nucleophilic centers of quinone 1, C2 and C9 carbons, respectively. Analysis of the Parr functions suggests that although the most favorable electrophilic/nucleophilic interaction is that involving the C2 carbon of quinone and the C12 carbon of the butadiene framework, the intramolecular nature of the cyclization prevents the corresponding reactive channel.

On the catalytic effect of water in the intramolecular Diels­Alder reaction of quinone systems: a theoretical study.

The mechanism of the intramolecular Diels–Alder (IMDA) reaction of benzoquinone 1, in the absence and in the presence of three water molecules, 1w, has been studied by means of density functional theory (DFT) methods, using the M05-2X and B3LYP functionals for exploration of the potential energy surface (PES). The energy and geometrical results obtained are complemented with a population analysis using the NBO method, and an analysis based on the global, local and group electrophilicity and nucleophilicity indices. Both implicit and explicit solvation emphasize the increase of the polarity of the reaction and the reduction of activation free energies associated with the transition states (TSs) of this IMDA process. These results are reinforced by the analysis of the reactivity indices derived from the conceptual DFT, which show that the increase of the electrophilicity of the quinone framework by the hydrogen-bond formation correctly explains the high polar character of this intramolecular process. Large polarization at the TSs promoted by hydrogen-bonds and implicit solvation by water together with a high electrophilicity-nucleophilicity difference consistently explains the catalytic effects of water molecules.

Bond Fukui functions as descriptor of the electron density reorganization in π conjugated systems.

The bond Fukui function is introduced and tested as a new reactivity index capable of predicting the evolution of bond breaking and formation processes during an organic reaction involving π conjugated systems. As an illustration, we examine many cases where substituted ethylenes and dienes may respond to different reagents to yield cycloaddition, Michael addition, and other reactions at double bonds.

Quantitative characterization of group electrophilicity and nucleophilicity for intramolecular Diels-Alder reactions.

In a previous work (L. R. Domingo, M. J. Aurell, P. Perez and R. Contreras, Tetrahedron 2002, 58, 4417) we proposed that the difference in global electrophilicity index be taken as a measure of the polarity at the transition state in intermolecular Diels-Alder reactions. We herein extend this model to deal with intramolecular Diels-Alder (IMDA) processes. The transferability of the empirical reactivity rules established for the intermolecular DA reactions to the IMDA reactions is discussed. The analysis based on group electrophilicity and nucleophilicity in general fails because having two different reactivity patterns within the same molecule hampers a clean classification of electrophilicity and nucleophilicity of the interacting fragments. We introduce dual philicity indexes E1 and E2 that solve this problem by separating a series of 30 IMDA reactions into two families, namely the diene to dienophile electron flow (DDpF) and the dienophile to diene electron flow (DpDF) processes. The new indexes correctly describe the charge transfer at the transition state and the reaction mechanism expected for the title reactions.

Microbial transformation of marine halogenated sesquiterpenes.

The sesquiterpene pacifenol is one of the main constituents of the red alga Laurencia claviformis. Earlier work on the semisynthetic derivatives of pacifenol afforded a series of halogenated sesquiterpenes. The aim of the present work was to obtain new hydroxylated derivatives of halogenated sesquiterpenes by means of microbial transformation using Aspergillus niger, Gibberella fujikuroi and Mucor plumbeus. The best results were obtained with M. plumbeus. The microbiological transformation by M. plumbeus of pacifenol, and two semisynthetic derivatives, is described. The structures of the new compounds obtained were determined by spectroscopic means.

An unusual halogenated meroditerpenoid from Stypopodium flabelliforme: studies by NMR spectroscopic and computational methods.

Meroditerpenoids, 2-[2'(E)-3',7',11',15'-tetramethylhexadec-2-en-1'-yl]-6-methyl-1,4-benzohydroquinone diacetate and 4'-chlorostypotriol triacetate, along with eight known compounds isolated from the dichloromethane extract of the brown alga Stypopodium flabelliforme after peracetylation are reported. One of them, 2-(1-oxo-hexadecyl)-1,3,5-trihydroxybenzene, is described for the first time within this genus. Structural elucidation was carried out on the basis of spectroscopic data and theoretical studies using GIAO/DFT analysis at B3LYP/6-31G(d) and mPW1PW91/6-31G(d) levels of theory for 4'-chlorostypotriol. This isomer is the first metabolite from the Stypopodium genus possessing one halogen atom.

Antioxidant properties and free radical-scavenging reactivity of a family of hydroxynaphthalenones and dihydroxyanthracenones.

This study was undertaken to investigate the free radical-scavenging and antioxidant activities of various structurally related hydroquinones including hydroxynaphthalenones and dihydroxyanthracenones. Electron spin resonance spectroscopy and spin trapping techniques were used to evaluate the ability of hydroquinones to scavenge hydroxyl, diphenylpicrylhydrazyl, and galvinoxyl radicals. In addition, the oxygen radical absorbing capacity assay using fluorescein (ORAC-FL) was used to obtain the relative antioxidant capacity of these radicals. The rate constants of the first H atom abstraction by 2,2-diphenyl-2-picrylhydrazyl (k(2)), were obtained under pseudo-first-order conditions. The free radical-scavenging activities and k(2) values discriminate well between hydroxynaphthalenones and dihydroxyanthracenones, showing that the latter have better antioxidant properties. The aforementioned experimental data agree with quantum-chemical results demonstrating the relevance of intramolecular H bonding to radical-scavenging activities.