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.     

Structure-activity relationships and docking studies of synthetic 2-arylindole derivatives determined with aromatase and quinone reductase 1

In our ongoing effort of discovering anticancer and chemopreventive agents, a series of 2-arylindole derivatives were synthesized and evaluated toward aromatase and quinone reductase 1 (QR1). Biological evaluation revealed that several compounds (e.g., 2d, IC₅₀?=?1.61?μM; 21, IC₅₀?=?3.05?μM; and 27, IC₅₀?=?3.34?μM) showed aromatase inhibitory activity with half maximal inhibitory concentration (IC₅₀) values in the low micromolar concentrations. With regard to the QR1 induction activity, 11 exhibited the highest QR1 induction ratio (IR) with a low concentration to double activity (CD) value (IR?=?8.34, CD?=?2.75?μM), while 7 showed the most potent CD value of 1.12?μM. A dual acting compound 24 showed aromatase inhibition (IC₅₀?=?9.00?μM) as well as QR1 induction (CD?=?5.76?μM) activities. Computational docking studies using CDOCKER (Discovery Studio 3.5) provided insight in regard to the potential binding modes of 2-arylindoles within the aromatase active site. Predominantly, the 2-arylindoles preferred binding with the 2-aryl group toward a small hydrophobic pocket within the active site. The C-5 electron withdrawing group on indole was predicted to have an important role and formed a hydrogen bond with Ser478 (OH). Alternatively, meta-pyridyl analogs may orient with the pyridyl 3′-nitrogen coordinating with the heme group.     

Molecular dynamics simulation of a myosin subfragment-1 docking with an actin filament

Myosins are typical molecular motor proteins, which convert the chemical energy of ATP into mechanical work. The fundamental mechanism of this energy conversion is still unknown. To explain the experimental results observed in molecular motors, Masuda has proposed a theory called the “Driven by Detachment (DbD)” mechanism for the working principle of myosins. Based on this theory, the energy used during the power stroke of the myosins originates from the attractive force between a detached myosin head and an actin filament, and does not directly arise from the energy of ATP. According to this theory, every step in the myosin working process may be reproduced by molecular dynamics (MD) simulations, except for the ATP hydrolysis step. Therefore, MD simulations were conducted to reproduce the docking process of a myosin subfragment-1 (S1) against an actin filament. A myosin S1 directed toward the barbed end of an actin filament was placed at three different positions by shifting it away from the filament axis. After 30 ns of MD simulations, in three cases out of ten trials on average, the myosin made a close contact with two actin monomers by changing the positions and the orientation of both the myosin and the actin as predicted in previous studies. Once the docking was achieved, the distance between the myosin and the actin showed smaller fluctuations, indicating that the docking is stable over time. If the docking was not achieved, the myosin moved randomly around the initial position or moved away from the actin filament. MD simulations thus successfully reproduced the docking of a myosin S1 with an actin filament. By extending the similar MD simulations to the other steps of the myosin working process, the validity of the DbD theory may be computationally demonstrated.     

Acceleration of azo dye decolorization by using quinone reductase activity of azoreductase and quinone redox mediator

This study demonstrated the effective application of intracellular azoreductase in mediated decolorization of azo dyes. Using the quinone reductase activity of overexpressed azoreductase AZR and quinone redox mediators, the decolorization performance of the recombinant strain Escherichia coli YB was significantly enhanced. In the presence of 0.2 mM lawsone, 75% acid red 27 (1 mM) was decolorized by E. coli YB in only 2 h, which was the highest bacterial decolorization rate ever reported. Compared to lawsone, menadione was a less effective redox mediator. Glucose was found to be the best carbon source for mediated decolorization by E. coli YB. The recombinant strain could complete four rounds of mediated decolorization repeatedly in 12 h. In addition, a 10-min pre-incubation of E. coli JM109 and activated sludge with 2-methylhydroquinone resulted in great improvement of mediated decolorization performance, which may be applied in practical treatment.     

Development of a potent 2-oxoamide inhibitor of secreted phospholipase A2 guided by molecular docking calculations and molecular dynamics simulations

Inhibition of group IIA secreted phospholipase A₂ (GIIA sPLA₂) has been an important objective for medicinal chemists. We have previously shown that inhibitors incorporating the 2-oxoamide functionality may inhibit human and mouse GIIA sPLA₂s. Herein, the development of new potent inhibitors by molecular docking calculations using the structure of the known inhibitor 7 as scaffold, are described. Synthesis and biological evaluation of the new compounds revealed that the long chain 2-oxoamide based on (S)-valine GK241 led to improved activity (IC₅₀ = 143 nM and 68 nM against human and mouse GIIA sPLA₂, respectively). In addition, molecular dynamics simulations were employed to shed light on GK241 potent and selective inhibitory activity.     

Effects of hydroxysafflor yellow A on ALDH1: Inhibition kinetics and molecular dynamics simulation

Hydroxysafflor yellow A (HSYA) is a potent natural antioxidant that displays important neuroprotective activity. Inhibition of acetaldehyde dehydrogenase 1 (ALDH1) has attracted the attention of researchers due to its overexpression in several types of cancers. We studied the effects of HSYA on ALDH1 by evaluating the inhibitory kinetics based on its antioxidant properties and performing computational simulation integrating methods. HSYA reversibly inhibited human recombinant ALDH1 via non-competitive inhibition (K_i = 0.267 ± 0.024 mM). We also investigated the tertiary structural changes via measuring intrinsic and ANS-binding fluorescence. The results indicated that the inactivation induced by HSYA was associated with structural changes. To obtain further information, we simulated the 3D structure of ALDH1 and conducted computational docking simulations as well as molecular dynamics simulations. The results indicated that 4 rings of HSYA interact with several residues near the ALDH1 active site. Our study provides insight into the inhibition of ALDH1 accompanied by structural changes. Based on its ALDH1-inhibiting effect and its potential as a natural antioxidant, HSYA is a potential agent for treating ALDH1-associated cancers.     

Structural study of biologically significant ligands with major birch pollen allergen Betv1 by docking and molecular dynamics simulation

The major birch pollen allergen, Betv1 of Betula verrucosa is the main causative agent of birch pollen allergy in humans. Betv1 is capable of binding several physiological ligands including fatty acids, flavones, cytokinins and sterols. Until now, no structural information from crystallography or NMR is available regarding binding mode of any of these ligands into the binding pocket of Betv1. In the present study thirteen ligands have been successfully docked into the hydrophobic cavity of Betv1 and binding free energies of the complexes have been calculated using AutoDock 3.0.5. A linear relationship with correlation coefficient (R2) of 0.6 is obtained between ΔG(b)s values plotted against their corresponding IC50 values. The complex formed between Betv1 and the best docking pose for each ligand has been optimized by molecular dynamics simulation. Here, we describe the ligand binding of Betv1, which provides insight into the biological function of this protein. This knowledge is required for structural alteration or inhibition of some of these ligands in order to modify the allergenic properties of this protein.     

Molecular docking and molecular dynamics simulation approaches for identifying new lead compounds as potential AChE inhibitors

To provide hints for the design of novel acetylcholinesterase (AChE) inhibitors with higher potency and specificity, the binding modes of the (RS, S)-17b and (RS, R)-17b enantiomers on AChE were chosen to investigate by molecular docking and molecular dynamics simulation. The results show that the binding modes of (RS, S)-17b and (RS, R)-17b are clearly different from each other. In particular, the (RS, S)-17b and (RS, R)-17b enantiomers tend to be planar and bend conformations to interact with AChE, respectively. Furthermore, based on the binding mode on AChE and structure modification of (RS, S)-17b, two novel inhibitors (1 and 2) with higher inhibitory activity were designed. Our design strategy suggests that the number of N and O atoms should be increased, the 5, 6-dimethoxy should be transformed into ring and the indanone moiety should be ring-opening, which would result in generating potent and selective AChE inhibitors.     

Reduced catalytic activity of human CYP2C9 natural alleles for gliclazide: molecular dynamics simulation and docking studies

Amongst sulfonylureas, gliclazide is one of the mostly prescribed drugs to diabetic patients and is metabolized extensively by P450 CYP2C9. Among 24-CYP2C9 alleles, the *2/*2 and *3/*3 genotypes showed significantly lower gliclazide clearances with reductions of 25 and 57%, respectively. However, the reason for the change in drug-metabolizing activity induced by these natural alleles is unknown. In the present study, we used molecular dynamics simulation and autodocking studies to provide models for gliclazide-bound complexes of CYP2C9*2, *3 and *2/*3 mutants, which give insight into CYP2C9-gliclazide interactions and explain the reduced enzymatic activity seen in these variants. Our data shows that the size of the substrate-access entry site is significantly reduced in mutants, which limits the access of gliclazide to heme and the active site. The distance from the substrate oxidation site and heme is >5 Å in *3 and *2/*3. Therefore, the addition of an active oxygen molecule by heme-Fe is hindered. The absence of F100, F114 and F476 in the interacting amino acid pocket in *3 reduces catalytic efficiency toward gliclazide. In *1, gliclazide is stabilized by the formation of two hydrogen bonds with R108 while it is absent in mutants. Further in *3 and *2/*3, the key heme-stabilizing residue, R97 stabilization is greatly reduced. Therefore, the decreased catalytic activity of these variants can be explained from the reduced access of the gliclazide to heme, and the interaction between heme and substrate is affected due to their instability in the active site.     

A fast virtual screening approach to identify structurally diverse inhibitors of trypanothione reductase

Trypanothione reductase (TryR) is one of the favorite targets for those designing drugs for the treatment of Chagas disease. We present the application of a fast virtual screening approach for designing hit compounds active against TryR. Our protocol combines information derived from structurally known inhibitors and from the TryR receptor structure. Five structurally diverse hit compounds active against TryR and holding promise for the treatment of Chagas disease are reported.     

Design,synthesis, molecular docking and biological evaluation of thiophen-2-iminothiazolidine derivatives for use against Trypanosoma cruzi

In this study, we designed and synthesized a series of thiophen-2-iminothiazolidine derivatives from thiophen-2-thioureic with good anti-Trypanosoma cruzi activity. Several of the final compounds displayed remarkable trypanocidal activity. The ability of the new compounds to inhibit the activity of the enzyme cruzain, the major cysteine protease of T. cruzi, was also explored. The compounds 3b, 4b, 8b and 8c were the most active derivatives against amastigote form, with significant IC₅₀ values between 9.7 and 6.03 μM. The 8c derivative showed the highest potency against cruzain (IC₅₀ = 2.4 μM). Molecular docking study showed that this compound can interact with subsites S1 and S2 simultaneously, and the negative values for the theoretical energy binding (E_b = −7.39 kcal·mol⁻¹) indicates interaction (via dipole–dipole) between the hybridized sulfur sp³ atom at the thiazolidine ring and Gly66. Finally, the results suggest that the thiophen-2-iminothiazolidines synthesized are important lead compounds for the continuing battle against Chagas disease.     

Thermal expansion of ZnSe by molecular dynamics simulation

Thermal expansion of bulk ZnSe crystal, grown by physical vapour transport technique, is experimentally measured for various temperatures and molecular dynamics simulation using a three body Tersoff potential is also developed for the same. The resulting expansivity compares well with our experimental results and other reported values.     

Interaction between dihydropyridines and phospholipid bilayers: a molecular dynamics simulation

Interaction of the calcium-channel antagonist dihydropyridines (DHPs), lacidipine and nifedipine, with a phospholipid bilayer was studied using 600 ps molecular dynamic simulations. We have constructed a double layer membrane model composed of 42 dimirystoyl-phosphatidylcholine molecules. The DHP molecules locate at about 7 Å from the centre of the membrane, inducing an asymmetry in the bilayer. While lacidipine did not induce significant local perturbations as judged by the gauche-trans isomerisation rate, nifedipine significantly decreased this rate, probably by producing a local rigidity of the membrane in the vicinity of the DHP.     

Docking and molecular dynamics simulation of the Azurin-Cytochrome c551 electron transfer complex

We coupled protein-protein docking procedure with molecular dynamics (MD) simulation to investigate the electron transfer (ET) complex Azurin-Cytochrome c551 whose transient character makes difficult a direct experimental investigation. The ensemble of complexes generated by the docking algorithm are filtered according to both the distance between the metal ions in the redox centres of the two proteins and to the involvement of suitable residues at the interface. The resulting best complex (BC) is characterized by a distance of 1.59 nm and involves Val23 and Ile59 of Cytochrome c551. The ET properties have been evaluated in the framework of the Pathways model and compared with experimental data. A 60 ns long MD simulation, carried on at full hydration, evidenced that the two protein molecules retain their mutual spatial positions upon forming the complex. An analysis of the ET properties of the complex, monitored at regular time intervals, has revealed that several different ET paths are possible, with the occasional intervening of water molecules. Furthermore, the temporal evolution of the geometric distance between the two redox centres is characterized by very fast fluctuations around an average value of 1.6 nm, with periodic jumps at 2 nm with a frequency of about 70 MHz. Such a behaviour is discussed in connection with a nonlinear dynamics of protein systems and its possible implications in the ET process are explored.     

Structure and dynamics of the antifungal molecules Syringotoxin-B and Syringopeptin-25A from molecular dynamics simulation

The phytopathogen Pseudomonas syringae pv. syringae produces toxic cyclic lipodepsipeptides (CLPs): nona-peptides and syringopeptins. All CLPs inhibit the growth of many fungal species, including human pathogens, although different fungi display different degrees of sensitivity. The best studied CLPs are Syringomycin-E (SR-E), Syringotoxin-B (ST-B) and Syringopeptin-25A (SP-25A). Their biological activity is affected by membrane composition and their structural differences. We previously (Matyus et al. in Eur Biophys J 35:459–467, 2006) reported the molecular features and structural preferences of SR-E in water and octane environments. Here we investigate in atomic detail the molecular features of the two other main CLP components, ST-B and SP-25A, in water and octane by 200 ns molecular dynamics simulations (MD), using distance restraints derived from NMR NOE data (Ballio et al. in Eur J Biochem 234:747–758, 1995). We have obtained three-dimensional models of ST-B and SP-25A CLPs in different environments. These models can now be used as a basis to investigate the interactions of ST-B and SP-25A with lipid membranes an important further step towards a better understanding of the antifungal and antibacterial activity of these peptides. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Combined docking and molecular dynamics simulations to enlighten the capacity of Pseudomonas cepacia and Candida antarctica lipases to catalyze quercetin acetylation

A combined docking and molecular dynamics protocol was applied to investigate quercetin binding modes within the catalytic cavity of Candida antarctica lipase B (CALB) and Pseudomonas cepacia lipase (PCL), aiming to explain the difference of specificity of these enzymes in acetylation reaction. For both lipases, docking of quercetin yielded two families of conformers with either the quercetin A or B-ring pointing towards the catalytic residues. Molecular dynamics (MD) calculations were subsequently performed on several complexes of each family. MD trajectories were analyzed focusing on the orientation of the acyl donor bound to the catalytic serine towards the oxyanion hole residues and the proximity of quercetin hydroxyl groups to the catalytic residues. Results showed that with CALB, the acetate was not correctly positioned within the oxyanion hole whatever the orientation of quercetin, suggesting that no product could be obtained. With PCL, the acetate remained within the oxyanion hole during all MD trajectories. Depending on quercetin orientation, either the 7-OH group or the 3, 5, 3′, 4′-OH groups came alternatively near the catalytic residues, suggesting that all of them could be acylated. The capacity of models to explain the regioselectivity of the reaction was discussed. Key residues and interactions involved in quercetin binding modes were identified and related to the reaction feasibility.     

Metronidazole aryloxy,carboxy and azole derivatives: Synthesis,anti-tumor activity,QSAR, molecular docking and dynamics studies

A series of novel metronidazole aryloxy, carboxy and azole derivatives has been synthesized and their cytotoxic activities on three cancer cell lines were evaluated by MTT assay. Compounds 4m, 4l and 4d showed the most potent cytotoxic activity (IC₅₀s?less than?100?µg/mL). Apoptosis was also detected for these compounds by flow cytometry. Docking studies were performed in order to propose the probable target protein. In the next step, molecular dynamics simulation was carried out on the proposed target protein, focal adhesion kinase (FAK, PDB code: 2ETM), bound to compound 4m. As, 4m showed a potent cytotoxic activity and an acceptable apoptotic effect, it can be a potential anticancer candidate that may work through inhibition of FAK.