Arylated analogues of cypronazole: fungicidal effect and activity on human fibroblasts. Docking analysis and molecular dynamics simulations.

Triadimefon (TDM) and cyproconazole (CPZ) are two triazoles widely used as fungicides. Several azoles were synthesised starting from commercial TDM and CPZ. The compounds were evaluated against phytopathogenic filamentous fungi, including Aspergillus fumigatus (AF), A. niger (AN), A. ustus (AU), A. japonicus (AJ), A. terreus (AT), Fusarium oxysporum and Botrytis cinerea isolated from grapevine in the province of San Juan, Argentina. Three of the synthesised compounds (1-(Biphenyl-4-yloxy)-3,3-dimethyl-1-(1H-1,2,4-triazol-1-yl)butan-2-one, 1; 2-(Biphenyl-4-yl)-3-cyclopropyl-1-(1H-1,2,4-triazol-1-yl)butan-2-ol, 3; 3-Cyclopropyl-2-(4'-fluorobiphenyl-4-yl)-1-(1H-1,2,4-triazol1-yl)butan-2-ol, 4) presented remarkable in vitro fungicidal properties, with better effects than TDM and CPZ on some of the target fungi. Cytotoxicity was assessed using human lung fibroblasts MRC5. Derivative 1, with IC50 values of 389.4 µM, was less toxic towards MRC-5 human lung fibroblasts than commercial TDM (248.5 µM) and CPZ (267.4 µM). Docking analysis and molecular dynamics simulations suggest that the compounds present the same interaction in the binding pocket of the CYP51B enzyme and with the same amino acids as CPZ. The derivatives investigated could be considered broad-spectrum but with some selectivity towards imperfect fungi.

Antiproliferative and antiangiogenic effects of ammonium tetrathiomolybdate in a model of endometriosis.

AIMS: Copper (Cu) is involved in the endometriosis progression. Herein, an experimental endometriosis model was used to evaluate whether its chelation with ammonium tetrathiomolybdate (TM) affects the proliferation and angiogenesis in endometriotic-like lesions and the participation of oxidative stress in these processes. MAIN METHODS: Female C57BL/6 mice were divided into three groups: sham-operated mice, endometriosis-induced mice, and TM-treated endometriosis-induced mice. Each animal in the third group received 0.3 mg of TM/day in their drinking water from the postoperative 15th day. The samples were collected after one month of induced pathology. In peritoneal fluids, Cu and estradiol levels were determined by electrothermal atomic absorption spectrometry and electrochemiluminescence, respectively. Endometriotic-like lesions were processed for the analysis of cell proliferation by PCNA immunohistochemistry, the expression of angiogenic markers by RT-qPCR, the presence of endothelial cells by immunofluorescent staining, and oxidative stress applying spectrophotometric methods. KEY FINDINGS: TM treatment decreased Cu and estradiol levels, which were increased by this pathology. In lesions, TM induced: (a) a decrease in tissue weight and volume, (b) a decrease in PCNA-positive cells, (c) antiangiogenic effects by decreasing the number of blood vessels, the mRNA expression of fibroblast growth factor 2 (Fgf2) and platelet-derived growth factor subunit B (Pdgfb), and the presence of endothelial cells, (d) a decrease in antioxidant activity and an increase in lipid peroxidation. SIGNIFICANCE: TM is a highly effective antiproliferative and antiangiogenic agent, modulating oxidative imbalance in endometriosis. Its anti-endometriotic potential is an attractive feature of TM as a possible non-hormonal treatment.

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.

Design of new quinolin-2-one-pyrimidine hybrids as sphingosine kinases inhibitors.

Sphingosine-1-phosphate is now emerging as an important player in cancer, inflammation, autoimmune, neurological and cardiovascular disorders. Abundance evidence in animal and humans cancer models has shown that SphK1 is linked to cancer. Thus, there is a great interest in the development new SphK1 inhibitors as a potential new treatment for cancer. In a search for new SphK1 inhibitors we selected the well-known SKI-II inhibitor as the starting structure and we synthesized a new inhibitor structurally related to SKI-II with a significant but moderate inhibitory effect. In a second approach, based on our molecular modeling results, we designed new structures based on the structure of PF-543, the most potent known SphK1 inhibitor. Using this approach, we report the design, synthesis and biological evaluation of a new series of compounds with inhibitory activity against both SphK1 and SphK2. These new inhibitors were obtained incorporating new connecting chains between their polar heads and hydrophobic tails. On the other hand, the combined techniques of molecular dynamics simulations and QTAIM calculations provided complete and detailed information about the molecular interactions that stabilize the different complexes of these new inhibitors with the active sites of the SphK1. This information will be useful in the design of new SphK inhibitors.

Theoretical models to predict the inhibitory effect of ligands of sphingosine kinase 1 using QTAIM calculations and hydrogen bond dynamic propensity analysis.

We report here the results of two theoretical models to predict the inhibitory effect of inhibitors of sphingosine kinase 1 that stand on different computational basis. The active site of SphK1 is a complex system and the ligands under the study possess a significant conformational flexibility; therefore for our study we performed extended simulations and proper clusterization process. The two theoretical approaches used here, hydrogen bond dynamics propensity analysis and Quantum Theory of Atoms in Molecules (QTAIM) calculations, exhibit excellent correlations with the experimental data. In the case of the hydrogen bond dynamics propensity analysis, it is remarkable that a rather simple methodology with low computational requirements yields results in excellent accord with experimental data. In turn QTAIM calculations are much more computational demanding and are also more complex and tedious for data analysis than the hydrogen bond dynamic propensity analysis. However, this greater computational effort is justified because the QTAIM study, in addition to giving an excellent correlation with the experimental data, also gives us valuable information about which parts or functional groups of the different ligands are those that should be replaced in order to improve the interactions and thereby to increase the affinity for SphK1. Our results indicate that both approaches can be very useful in order to predict the inhibiting effect of new compounds before they are synthesized.

Tetrahydroisoquinolines functionalized with carbamates as selective ligands of D2 dopamine receptor.

A series of tetrahydroisoquinolines functionalized with carbamates is reported here as highly selective ligands on the dopamine D2 receptor. These compounds were selected by means of a molecular modeling study. The studies were carried out in three stages: first an exploratory study was carried out using combined docking techniques and molecular dynamics simulations. According to these results, the bioassays were performed; these experimental studies corroborated the results obtained by molecular modeling. In the last stage of our study, a QTAIM analysis was performed in order to determine the main molecular interactions that stabilize the different ligand-receptor complexes. Our results show that the adequate use of combined simple techniques is a very useful tool to predict the potential affinity of new ligands at dopamine D1 and D2 receptors. In turn the QTAIM studies show that they are very useful to evaluate in detail the molecular interactions that stabilize the different ligand-receptor complexes; such information is crucial for the design of new ligands.

Molecular design and synthesis of novel peptides from amphibians skin acting as inhibitors of cholinesterase enzymes.

Cholinesterases are a family of enzymes that catalyze the hydrolysis of neurotransmitter acetylcholine. There are two types of cholinesterases, acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), which differ in their distribution in the body. Currently, cholinesterase inhibitors (ChEI) represent the treatment of choice for Alzheimer's disease (AD). In this paper, we report the synthesis and inhibitory effect on both enzymes of four new peptides structurally related to P1-Hp-1971 (amphibian skin peptide found in our previous work. Sequence: TKPTLLGLPLGAGPAAGPGKR-NH2 ). The bioassay data and cytotoxicity test show that some of the compounds possess a significant AChE and BChE inhibition and no toxic effect. The present work demonstrates that diminution of the size of the original peptide could potentially result in new compounds with significant cholinesterase inhibition activity, although it appears that there is an optimal size for the sequence. We also conducted an exhaustive molecular modeling study to better understand the mechanism of action of these compounds by combining docking techniques with molecular dynamics simulations on BChE. This is the first report about amphibian peptides and the second one of natural peptides with ChE inhibitory activity. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.

Alkaloids from Hippeastrum argentinum and Their Cholinesterase-Inhibitory Activities: An in Vitro and in Silico Study.

Two new alkaloids, 4-O-methylnangustine (1) and 7-hydroxyclivonine (2) (montanine and homolycorine types, respectively), and four known alkaloids were isolated from the bulbs of Hippeastrum argentinum, and their cholinesterase-inhibitory activities were evaluated. These compounds were identified using GC-MS, and their structures were defined by physical data analysis. Compound 2 showed weak butyrylcholinesterase (BuChE)-inhibitory activity, with a half-maximal inhibitory concentration (IC50) value of 67.3 ± 0.09 µM. To better understand the experimental results, a molecular modeling study was also performed. The combination of a docking study, molecular dynamics simulations, and quantum theory of atoms in molecules calculations provides new insight into the molecular interactions of compound 2 with BuChE, which were compared to those of galantamine.

Theoretical and Experimental Study of Inclusion Complexes of ß-Cyclodextrins with Chalcone and 2',4'-Dihydroxychalcone.

The inclusion complexes formed by chalcone and 2',4'-dihydroxychalcone with ß-cyclodextrin have been studied combining experimental (phase solubility diagrams, Fourier transform infrared spectroscopy) and molecular modeling (molecular dynamics, quantum mechanics/molecular mechanics calculations) techniques. The formation constants of the complexes were determined at different temperatures, and the thermodynamic parameters of the process were obtained. The inclusion of chalcone in ß-cyclodextrin is an exothermic process, while the inclusion of 2',4'-dihydroxychalcone is endothermic. Free energy profiles, derived from umbrella sampling using molecular dynamics simulations, were constructed to analyze the binding affinity and the complexation reaction at a molecular level. Hybrid QM/MM calculations were also employed to obtain a better description of the energetic and structural aspects of the complexes. The intermolecular interactions that stabilize both inclusion complexes were characterized by means of quantum atoms in molecules theory and reduce density gradient method. The calculated interactions were experimentally observed using FTIR.

Small Peptides Derived from Penetratin as Antibacterial Agents.

The synthesis, in vitro evaluation and conformational study of several small-size peptides acting as antibacterial agents are reported. Among the compounds evaluated, the peptides Arg-Gln-Ile-Lys-Ile-Trp-Arg-Arg-Met-Lys-Trp-Lys-Lys-NH2 , Arg-Gln-Ile-Lys-Ile-Arg-Arg-Met-Lys-Trp-Arg-NH2 , and Arg-Gln-Ile-Trp-Trp-Trp-Trp-Gln-Arg-NH2 exhibited significant antibacterial activity. These were found to be very active antibacterial compounds, considering their small molecular size. In order to better understand the antibacterial activity obtained for these peptides, an exhaustive conformational analysis was performed, using both theoretical calculations and experimental measurements. Molecular dynamics simulations using two different media (water and trifluoroethanol/water) were employed. The results of these theoretical calculations were corroborated by experimental circular dichroism measurements. A brief discussion on the possible mechanism of action of these peptides at molecular level is also presented. Some of the peptides reported here constitute very interesting structures to be used as starting compounds for the design of new small-size peptides possessing antibacterial activity.

New mimetic peptides inhibitors of Αß aggregation. Molecular guidance for rational drug design.

A new series of mimetic peptides possessing a significant Aß aggregation modulating effect was reported here. These compounds were obtained based on a molecular modelling study which allowed us to perform a structural-based virtual selection. Monitoring Aß aggregation by thioflavin T fluorescence and transmission electron microscopy revealed that fibril formation was significantly decreased upon prolonged incubation in presence of the active compounds. Dot blot analysis suggested a decrease of soluble oligomers strongly associated with cognitive decline in Alzheimer's disease. For the molecular dynamics simulations, we used an Aß42 pentameric model where the compounds were docked using a blind docking technique. To analyze the dynamic behaviour of the complexes, extensive molecular dynamics simulations were carried out in explicit water. We also measured parameters or descriptors that allowed us to quantify the effect of these compounds as potential inhibitors of Aß aggregation. Thus, significant alterations in the structure of our Aß42 protofibril model were identified. Among others we observed the destruction of the regular helical twist, the loss of a stabilizing salt bridge and the loss of a stabilizing hydrophobic interaction in the ß1 region. Our results may be helpful in the structural identification and understanding of the minimum structural requirements for these molecules and might provide a guide in the design of new aggregation modulating ligands.

A new series of antibacterial nitrosopyrimidines: synthesis and structure-activity relationship.

New nitrosopyrimidines were synthesized and evaluated as potential antibacterial agents. Different compounds structurally related with 4,6-bis(alkyl or arylamino)-5-nitrosopyrimidines were evaluated. Some of these nitrosopyrimidines displayed significant antibacterial activity against human pathogenic bacteria. Among them compounds 1c, 2a-c, and 9a-c exhibited remarkable activity against methicillin-sensitive and -resistant Staphylococcus aureus, Escherichia coli, Yersinia enterocolitica, and Salmonella enteritidis. A detailed structure-activity relationship study, supported by theoretical calculations, aided us to identify and understand the minimal structural requirements for the antibacterial action of the nitrosopyrimidines reported here. Thus, our results have led us to identify a topographical template that provides a guide for the design of new nitrosopyrimidines with antibacterial effects.

Structural and functional insights into the anti-BACE1 Fab fragment that recognizes the BACE1 exosite.

A molecular modeling study giving structural, functional, and mutagenesis insights into the anti-BACE1 Fab fragment that recognizes the BACE1 exosite is reported. Our results allow extending experimental data resulting from X-ray diffraction experiments in order to examine unknown aspects for the Fab-BACE1 recognition and its binding mode. Thus, the study performed here allows extending the inherently static nature of crystallographic structures in order to gain a deeper understanding of the structural and dynamical basis at the atomic level. The characteristics and strength of the interatomic interactions involved in the immune complex formation are exhaustively analyzed. The results might explain how the anti-BACE1 Fab fragment and other BACE1 exosite binders are capable to produce an allosteric modulation of the BACE1 activity. Our site-directed mutagenesis study indicated that the functional anti-BACE1 paratope, residues Tyr32 (H1), Trp50 (H2), Arg98 (H3), Phe101 (H3), Trp104 (H3) and Tyr94 (L3), strongly dominates the binding energetics with the BACE1 exosite. The mutational studies described in this work might accelerate the development of new BACE1 exosite binders with interesting pharmacological activity.

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.

Synthesis, bioevaluation and structural study of substituted phthalazin-1(2H)-ones acting as antifungal agents.

Twenty-five polysubstituted phthalazinone derivatives were synthesized and tested for their antifungal activity against a panel of pathogenic and clinically important yeasts and filamentous fungi. Among them, the compound 4-(4-chlorobenzyl)-2-methylphthalazin-1(2H)-one (5) exhibited a remarkable antifungal activity against standardised strains of dermatophytes and Cryptococcus neoformans, as well as against some clinical isolates. A physicochemical study performed on compound 5 revealed its conformational and electronic characteristics, providing us with useful data for the future design of novel related antifungal analogues.

Catalytic and molecular properties of rabbit liver carboxylesterase acting on 1,8-cineole derivatives.

Rabbit liver carboxylesterase (rCE) was evaluated as the catalyst for the enantioselective hydrolysis of (+/-)-3-endo-acetyloxy-1 ,8-cineole [(+/-)-4], which yields (1S,3S,4R)-(+)-3-acetyloxy-1,8-cineole [(+)-4] and (1R,3R,4S)-(-)-3-hydroxy-1,8-cineole [(-)-3]. Enantioselective asymmetrization of meso-3,5-diacetoxy-1,8-cineol (5) gives (1S,3S,4R,5R)-(-)-3-acetyloxy-5-hydroxy-1,8-cineole (6), with high enantioselectivity. rCE has been chosen to perform both experiments and molecular modeling simulations. Docking simulations combined with molecular dynamics calculations were used to study rCE-catalyzed enantioselective hydrolysis of cineol derivatives. Both compounds were found to bind with their acetyl groups stabilized by hydrogen bond interactions between their oxygen atoms and Ser221.

Amyloid-ß fibril disruption by C60-molecular guidance for rational drug design.

The WHO has listed Alzheimer's disease among the major neurological disorders with an estimated 35 million people affected worldwide. Amyloid-ß is mostly believed to be the causative factor in Alzheimer's disease and the severity of the disease correlates with the tendency of amyloid-ß to form aggregation patterns-plaques. Lacking effective medication, the identification of any underlying mechanistic principles regarding plaque formation appears to be crucial. Here we carry out computer simulations to study the effect of C60 on structure and stability of an idealised pentameric construct of amyloid-ß units (a model fibril). A binding site on top of the structurally ordered stack of ß-sheets is identified that triggers structural alterations at the turn region of the hook-like ß-sheet assembly. Significant structural alterations are: (i) the destruction of regular helical twist, (ii) the loss of a stabilizing salt bridge and (iii) the loss of a stabilizing hydrophobic interaction close to the turn. Consequently, the main effect of C60 is the induction of sizable destabilization in native fibril structure. These structural insights may serve as a molecular guide for further rational drug design of effective inhibitors targeting fibril formation in Alzheimer's disease.

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.

Tetrahydroisoquinolines acting as dopaminergic ligands. A molecular modeling study using MD simulations and QM calculations.

A molecular modeling study on 16 1-benzyl tetrahydroisoquinolines (BTHIQs) acting as dopaminergic ligands was carried out. By combining molecular dynamics simulations with ab initio and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of BTHIQs interacting with the human dopamine D2 receptor (D2 DR) is reported here, providing a clear picture of the binding interactions of BTHIQs from both structural and energetic viewpoints. Molecular aspects of the binding interactions between BTHIQs and the D2 DR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental pKi values was obtained, predicting the potential dopaminergic effect of non-synthesized BTHIQs.