Molecular Modeling of Mycobacterium Tuberculosis DNA Gyrase and its Molecular Docking Study with Gatifloxacin Inhibitors

Determination of antioxidant, DNA protection, enzyme inhibition potential and molecular docking studies of a biomarker ursolic acid in Nepeta species

Ursolic acid (UA), which has many biological properties such as anti-cancer, anti-inflammatory and antioxidant, and regulates some pharmacological processes, has been isolated from the flowers, leaves, berries and fruits of many plant species. In this work, UA was purified from the methanol-chloroform crude extract of Nepeta species (N. aristata, N. baytopii, N. italica, N. trachonitica, N. stenantha) using a silica gel column with chloroform or ethyl acetate solvents via bioactivity-guided isolation. The most active sub-fractions were determined under bioactivities using antioxidant and DNA protection activities and enzyme inhibitions. UA was purified from these fractions and its structure was elucidated by NMR spectroscopy techniques. The highest amount of UA was found in N. stenantha (8.53 mg UA/g), while the lowest amount of UA was found in N. trachonitica (1.92 mg UA/g). The bioactivities of UA were evaluated with antioxidant and DNA protection activities, enzyme inhibitions, kinetics and interactions. The inhibition values (IC50) of α-amylase, α-glucosidase, urease, CA, tyrosinase, lipase, AChE, and BChE were determined between 5.08 and 181.96 µM. In contrast, Ki values of enzyme inhibition kinetics were observed between 0.04 and 0.20 mM. In addition, Ki values of these enzymes for enzyme-UA interactions were calculated as 0.38, 0.86, 0.45, 1.01, 0.23, 0.41, 0.01 and 2.24 µM, respectively. It is supported that UA can be widely used as a good antioxidant against oxidative damage, an effective DNA protector against genetic diseases, and a suitable inhibitor for metabolizing enzymes.Communicated by Ramaswamy H. Sarma.

Cathepsin K inhibitors based on 2-amino-1,3,4-oxadiazole derivatives

Two new series of hitherto unknown dipeptides, containing an electrophilic nitrile or a non-electrophilic 2-amino-1,3,4-oxadiazole moiety were synthesized and evaluated in vitro as Cathepsin K (Cat K) inhibitors. From 14 compounds obtained, the oxadiazole derivatives 10a, 10b, 10e, and 10g acted as enzymatic competitive inhibitors with Ki values between 2.13 and 7.33 µM. Molecular docking calculations were carried out and demonstrated that all inhibitors performed hydrogen bonds with residues from the enzyme active site, such as Asn18. The best inhibitors (10a, 10b, 10g) could also perform these bonds with Cys25, and 10a showed the most stabilizing interaction energy (-134.36 kcal mol^-1) with the active cavity. For the first time, derivatives based in 2-amino-1,3,4-oxadiazole scaffolds were evaluated, and the results suggested that this core displays a remarkable potential as a building block for Cat K inhibitors.

Experimental and theoretical affinity and catalysis studies between halogenated phenols and peroxidases: Understanding the bioremediation potential

Halogenated phenols, such as 2,4-dichlorophenol (2,4-DCP) and 4-bromophenol (4-BP) are pollutants generated by a various industrial sectors like chemical, dye, paper bleaching, pharmaceuticals or in an agriculture as pesticides. The use of Horseradish peroxidase (HRP) in the halogenated phenols treatment has already been mentioned, but it is not well understood how the different phenolic substrates can bind in the peroxidase active site nor how these specific interactions can influence in the bioremediation potential. In this work, different removal efficiencies were obtained for phenolic compounds investigated using HRP as catalyst (93.87 and 59.19% to 4BP and 2,4 DCP, respectively). Thus, to rationalize this result based on the interactions of phenols with active center of HRP, we combine computational and experimental methodologies. The theoretical approaches utilized include density functional theory (DFT) calculations, docking simulation and quantum mechanics/molecular mechanics (QM/MM) technique. Michaelis Menten constant (Km) obtained through experimental methodologies were 2.3 and 0.95 mM to 2,4-DCP and 4-BP, respectively, while the specificity constant (K_cat/K_m) found was 1.44 mM^-1 s^-1 and 0.62 mM^-1 s^-1 for 4-BP and 2,4-DCP, respectively. The experimental parameters appointed to the highest affinity of HRP to 4-BP. According to the molecular docking calculations, both ligands have shown stabilizing intermolecular interaction energies within the HRP active site, however, the 4-BP showed more stabilizing interaction energy (-53.00 kcal mol^-1) than 2,4-dichlorophenol (-49.23 kcal mol^-1). Besides that, oxidative mechanism of 4-BP and 2,4-DCP was investigated by the hybrid QM/MM approach. This study showed that the lowest activation energy values for transition states investigated were obtained for 4-BP. Therefore, by theoretical approach, the compound 4-BP showed the more stabilizing interaction and activation energy values related to the interaction within the enzyme and the oxidative reaction mechanism, respectively, which corroborates with experimental parameters obtained. The combination between experimental and theoretical approaches was essential to understand how the degradation potential of the HRP enzyme depends on the interactions between substrate and the active center cavity of the enzyme.

Understanding the Interaction Modes and Reactivity of Trimedoxime toward MmAChE Inhibited by Nerve Agents: Theoretical and Experimental Aspects

Organophosphorus (OP) compounds are used as both chemical weapons and pesticides. However, these agents are very dangerous and toxic to humans, animals, and the environment. Thus, investigations with reactivators have been deeply developed in order to design new antidotes with better efficiency, as well as a greater spectrum of action in the acetylcholinesterase (AChE) reactivation process. With that in mind, in this work, we investigated the behavior of trimedoxime toward the Mus musculus acetylcholinesterase (MmAChE) inhibited by a range of nerve agents, such as chemical weapons. From experimental assays, reactivation percentages were obtained for the reactivation of different AChE–OP complexes. On the other hand, theoretical calculations were performed to assess the differences in interaction modes and the reactivity of trimedoxime within the AChE active site. Comparing theoretical and experimental data, it is possible to notice that the oxime, in most cases, showed better reactivation percentages at higher concentrations, with the best result for the reactivation of the AChE–VX adduct. From this work, it was revealed that the mechanistic process contributes most to the oxime efficiency than the interaction in the site. In this way, this study is important to better understand the reactivation process through trimedoxime, contributing to the proposal of novel antidotes.

Slight difference in the isomeric oximes K206 and K203 makes huge difference for the reactivation of organophosphorus-inhibited AChE: Theoretical and experimental aspects

Studies with oximes have been extensively developed to design new reactivators with better efficiency, and greater spectrum of action. In this study, we aimed to analyze the influence of the Carbamoyl group position change in two isomeric oximes, K203 and K206, on the reactivation percentage of Mus musculus Acetylcholinesterase (MmAChE), inhibited by different nerve agents. Theoretical calculations were performed to assess the difference for the oxime activity with inhibited AChE-complexes and the factors that govern this difference. Comparing theoretical and experimental data, it is possible to observe that this change between the oximes results in different reactivation percentage for the same nerve agent, due to the different interaction modes and activation energy for the studied systems.

Development of technologies applied to the biodegradation of warfare nerve agents: Theoretical evidence for asymmetric homogeneous catalysis

Organophosphorus compounds have been widely employed to the development of warfare nerve agents and pesticides, resulting in a huge number of people intoxicated annually, being a serious problem of public health. Efforts worldwide have been done in order to design new technologies that are capable of combating or even reversing the poisoning caused by these OP nerve agents. In this line, the bioremediation arises as a promising and efficient alternative for this purpose. As an example of degrading enzymes, there is the organophosphate-degrading (OpdA) enzyme from Agrobacterium radiobacter, which has been quite investigated experimentally due to its high performance in the degradation of neurotoxic nerve agents. This work aims to look into the structural and electronic details that govern the interaction modes of these compounds in the OpdA active site, with the posterior hydrolysis reaction prediction. Our findings have brought about data about the OpdA performance towards different nerve agents, and among them, we may realize that the degradation efficiency strongly depends on the nerve agent structure and its stereochemistry, being in this case the compound Tabun the one more effectively hydrolyzed. By means of the chemical bonds (AIM) and orbitals (FERMO) analysis, it is suggested that the initial reactivity of the OP nerve agents in the OpdA active site does not necessarily dictate the reactivity and interaction modes over the reaction coordinate.

Asymmetric biodegradation of the nerve agents Sarin and VX by human dUTPase: chemometrics, molecular docking and hybrid QM/MM calculations

Organophosphorus compounds (OP) nerve agents are among the most toxic chemical substances known. Their toxicity is due to their ability to bind to acetylcholinesterase. Currently, some enzymes, such as phosphotriesterase, human serum paraoxonase 1 and diisopropyl fluorophosphatase, capable of degrading OP, have been characterized. Regarding the importance of bioremediation methods for detoxication of OP, this work aims to study the interaction modes between the human human deoxyuridine triphosphate nucleotidohydrolase (dUTPase) and Sarin and VX, considering their R_p and S_p enantiomers, to evaluate the asymmetric catalysis of those compounds. In previous work, this enzyme has shown good potential to degrade phosphotriesters, and based on this characteristic, we have applied the human dUTPase to the OP degradation. Molecular docking, chemometrics and mixed quantum and molecular mechanics calculations have been employed, showing a good interaction between dUTPase and OP. Two possible reaction mechanisms were tested, and according to our theoretical results, the catalytic degradation of OP by dUTPase can take place via both mechanisms, beyond being stereoselective, that is, dUTPase cleaves one enantiomer preferentially in relation to other. Chemometric techniques provided excellent assistance for performing this theoretical investigation. The dUTPase study shows importance by the fact of it being a human enzyme. Communicated by Ramaswamy H. Sarma.

Theoretical Studies Applied to the Evaluation of the DFPase Bioremediation Potential against Chemical Warfare Agents Intoxication

Organophosphorus compounds (OP) are part of a group of compounds that may be hazardous to health. They are called neurotoxic agents because of their action on the nervous system, inhibiting the acetylcholinesterase (AChE) enzyme and resulting in a cholinergic crisis. Their high toxicity and rapid action lead to irreversible damage to the nervous system, drawing attention to developing new treatment methods. The diisopropyl fluorophosphatase (DFPase) enzyme has been considered as a potent biocatalyst for the hydrolysis of toxic OP and has potential for bioremediation of this kind of intoxication. In order to investigate the degradation process of the nerve agents Tabun, Cyclosarin and Soman through the wild-type DFPase, and taking into account their stereochemistry, theoretical studies were carried out. The intermolecular interaction energy and other parameters obtained from the molecular docking calculations were used to construct a data matrix, which were posteriorly treated by statistical analyzes of chemometrics, using the PCA (Principal Components Analysis) multivariate analysis. The analyzed parameters seem to be quite important for the reaction mechanisms simulation (QM/MM). Our findings showed that the wild-type DFPase enzyme is stereoselective in hydrolysis, showing promising results for the catalytic degradation of the neurotoxic agents under study, with the degradation mechanism performed through two proposed pathways.

A newly developed oxime K203 is the most effective reactivator of tabun-inhibited acetylcholinesterase

BACKGROUND

Based on in vitro and in vivo rat experiments, the newly developed acetylcholinesterase (AChE) reactivator, K203, appears to be much more effective in the treatment of tabun poisonings than currently fielded oximes.

METHODS

To determine if this reactivating efficacy would extend to humans, studies were conducted in vitro using human brain homogenate as the source of AChE. The efficacy of K203 was compared with commercially available oximes; pralidoxime, obidoxime and asoxime (HI-6).

RESULTS

Reactivation studies showed that K203 was the most effective reactivator with a second order kinetic constant (kr) of 2142 min- 1. M- 1, which was 51 times higher than that obtained for obidoxime (kr = 42 min- 1. M- 1). Both pralidoxime and asoxime (HI-6) failed to significantly reactivate tabun-inhibited human AChE.

DISCUSSION

According to these results and previous studies, using K203, it appears that oxime K203 is the most effective reactivator of tabun-inhibited cholinesterase in several species including humans and should be considered as a possible medical countermeasure to tabun exposure.

Identification of biotransformation products of disperse dyes with rat liver microsomes by LC-MS/MS and theoretical studies with DNA: Structure-mutagenicity relationship using Salmonella/microsome assay

Azo dyes are known as a group of substances with DNA damage potential that depend on the nature and number of azo groups connected to aromatic rings (benzene and naphthalene), chemical properties, e.g. solubility and reactive functional groups, which significantly affect their toxicological and ecological risks. In this paper, we used in vitro models to evaluate the metabolism of selected textile dyes: Disperse Red 73 (DR 73), Disperse Red 78 (DR 78) and Disperse Red 167 (DR 167). To evaluate the mutagenic potential of the textile dyes, the Salmonella mutagenicity assay (Ames test) with strains TA 98 and TA 100 in the presence and absence of the exogenous metabolic system (S9) was used. DR73 was considered the most mutagenic compound, inducing both replacement base pairs (TA 100) and also changing frameshift (TA 98) mutations that are reduced in the presence of the S9 mixture. Furthermore, we used rat liver microsomes in the same experimental conditions of the S9 mixture to metabolize the dyes and the resultant solutions were analyzed using a liquid chromatography coupled to a quadrupole linear ion trap mass spectrometry (LC-MS/MS) to investigate the metabolites formed by the in vitro biotransformation. Based on this experiment, we detected and identified two biotransformation products for each textile dye substrate analyzed. Furthermore, to evaluate the interaction and reactivity of these compounds with DNA, theoretical calculations were also carried out. The results showed that the chemical reaction occurred preferentially at the azo group and the nitro group, indicating that there was a reduction in these groups by the CYP P450 enzymes presented in the rat microsomal medium. Our results clearly demonstrated that the reduction of these dyes by biological systems is a great environmental concern due to increased genotoxicity for the body of living beings, especially for humans.

Peroxoniobium inhibits leukemia cell growth

A new class of polyoxoniobate complex has been synthesized and characterized as a novel anticancer agent for photodynamic therapy. The complex inhibits the growth of chronic myelogenous leukemia cells with an IC₅₀ value of 30 μM, in the dark. However, upon exposure to light (365 nm) there is a fivefold increase in the cytotoxic activity. Light radiation activate the complex with the formation of radical species capable of interacting with DNA according to our experimental and theoretical data.

A new class of polyoxoniobate complex has been synthesized and characterized as a novel anticancer agent for photodynamic therapy.

Can Inhibitors of Snake Venom Phospholipases A₂ Lead to New Insights into Anti-Inflammatory Therapy in Humans? A Theoretical Study

Human phospholipase A₂ (hPLA₂) of the IIA group (HGIIA) catalyzes the hydrolysis of membrane phospholipids, producing arachidonic acid and originating potent inflammatory mediators. Therefore, molecules that can inhibit this enzyme are a source of potential anti-inflammatory drugs, with different action mechanisms of known anti-inflammatory agents. For the study and development of new anti-inflammatory drugs with this action mechanism, snake venom PLA₂ (svPLA₂) can be employed, since the svPLA₂ has high similarity with the human PLA₂ HGIIA. Despite the high similarity between these secretory PLA₂s, it is still not clear if these toxins can really be employed as an experimental model to predict the interactions that occur with the human PLA₂ HGIIA and its inhibitors. Thus, the present study aims to compare and evaluate, by means of theoretical calculations, docking and molecular dynamics simulations, as well as experimental studies, the interactions of human PLA₂ HGIIA and two svPLA₂s,Bothrops toxin II and Crotoxin B (BthTX-II and CB, respectively). Our theoretical findings corroborate experimental data and point out that the human PLA₂ HGIIA and svPLA₂ BthTX-II lead to similar interactions with the studied compounds. From our results, the svPLA₂ BthTX-II can be used as an experimental model for the development of anti-inflammatory drugs for therapy in humans.

Structural Basis and Mechanism of Chiral Benzedrine Molecules Interacting With Third Dopamine Receptor

In order to investigate the chiral benzedrine molecules corresponding to their different characteristics in biochemical systems, we studied their interaction with D3 R using the docking method, molecular dynamic simulation, and quantum chemistry. The obtained results indicate that the active residues for R-benzedrine (RAT) bound with D3 R are Ala132, Asp133, and Tyr55, while Asn57, Asp133, Asp168, Cys172, Gly54, Trp24, and Vall136 act as the active residues for S-benzedrine (SAT). The different active pockets are observed for ART or SAT because they possess different active residues. The binding energies between RAT and SAT with D3 R were determined to be -44.0 kJ.mol(-1) and -71.2 kJ.mol(-1) , respectively. These results demonstrate that SAT within the studied pocket of D3 R has a stronger capability of binding with D3 R, while it is more feasible for RAT to leave from the interior positions of D3 R. In addition, the results suggest that the D3 R protein can recognize chiral benzedrine molecules and influence their different addictive and pharmacological effects in biochemical systems. Chirality 28:674-685, 2016. © 2016 Wiley Periodicals, Inc.

A comparative modeling and molecular docking study on Mycobacterium tuberculosis targets involved in peptidoglycan biosynthesis

An alarming rise of multidrug-resistant Mycobacterium tuberculosis strains and the continuous high global morbidity of tuberculosis have reinvigorated the need to identify novel targets to combat the disease. The enzymes that catalyze the biosynthesis of peptidoglycan in M. tuberculosis are essential and noteworthy therapeutic targets. In this study, the biochemical function and homology modeling of MurI, MurG, MraY, DapE, DapA, Alr, and Ddl enzymes of the CDC1551 M. tuberculosis strain involved in the biosynthesis of peptidoglycan cell wall are reported. Generation of the 3D structures was achieved with Modeller 9.13. To assess the structural quality of the obtained homology modeled targets, the models were validated using PROCHECK, PDBsum, QMEAN, and ERRAT scores. Molecular dynamics simulations were performed to calculate root mean square deviation (RMSD) and radius of gyration (Rg) of MurI and MurG target proteins and their corresponding templates. For further model validation, RMSD and Rg for selected targets/templates were investigated to compare the close proximity of their dynamic behavior in terms of protein stability and average distances. To identify the potential binding mode required for molecular docking, binding site information of all modeled targets was obtained using two prediction algorithms. A docking study was performed for MurI to determine the potential mode of interaction between the inhibitor and the active site residues. This study presents the first accounts of the 3D structural information for the selected M. tuberculosis targets involved in peptidoglycan biosynthesis.

Computer-Aided Drug Discovery Approaches against the Tropical Infectious Diseases Malaria, Tuberculosis, Trypanosomiasis, and Leishmaniasis

Despite the tremendous improvement in overall global health heralded by the adoption of the Millennium Declaration in the year 2000, tropical infections remain a major health problem in the developing world. Recent estimates indicate that the major tropical infectious diseases, namely, malaria, tuberculosis, trypanosomiasis, and leishmaniasis, account for more than 2.2 million deaths and a loss of approximately 85 million disability-adjusted life years annually. The crucial role of chemotherapy in curtailing the deleterious health and economic impacts of these infections has invigorated the search for new drugs against tropical infectious diseases. The research efforts have involved increased application of computational technologies in mainstream drug discovery programs at the hit identification, hit-to-lead, and lead optimization stages. This review highlights various computer-aided drug discovery approaches that have been utilized in efforts to identify novel antimalarial, antitubercular, antitrypanosomal, and antileishmanial agents. The focus is largely on developments over the past 5 years (2010-2014).

Theoretical insight into the structural mechanism for the binding of vinblastine with tubulin

Vinblastine (VLB) is one of vinca alkaloids with high cytotoxicity toward cancer cells approved for clinical use. However, because of drug resistance, toxicity, and other side effects caused from the use of VLB, new vinca alkaloids with higher cytotoxicity toward cancer cells and other good qualities need to develop. One strategy is to further study and better understand the essence why VLB possesses the high cytotoxicity toward cancer cells. In present work, by using molecular simulation, molecular docking, density functional calculation, and the crystal structure of α,β-tubulin complex, we find two modes labeled in catharanthine moiety (CM) and vindoline moiety (VM) modes of VLB bound with the interface of α,β-tubulin to probe the essence why VLB has the high cytotoxicity toward cancer cells. In the CM mode, nine key residues B-Ser178, B-Asp179, B-Glu183, B-Tyr210, B-Asp226, C-Lys326, C-Asp327, C-Lys336, and C-Lys352 from the α,β-tubulin complex are determined as the active sites for the interaction of VLB with α,β-tubulin. Some of them such as B-Ser178, B-Glu183, B-Tyr210, B-Asp226, C-Lys326, C-Asp327, and C-Lys336 are newly identified as the active sites in present work. The affinity between VLB and the active pocket within the interface of α,β-tubulin is -60.8 kJ mol(-1) in the CM mode. In the VM mode, that is a new mode established in present paper, nine similar key residues B-Lys176, B-Ser178, B-Asp179, B-Glu183, B-Tyr210, B-Asp226, C-Lys326, C-Asp327, and C-Lys336 from the α,β-tubulin complex are found as the active sites for the interaction with VLB. The difference is from one key residue C-Lys352 in the CM mode changed to the key residue B-Lys176 in the VM mode. The affinity between VLB and the active pocket within the interface of α,β-tubulin is -96.3 kJ mol(-1) in the VM mode. Based on the results obtained in present work, and because VLB looks like two faces, composed of CM and VM both to have similar polar active groups, to interact with the active sites, we suggest double-faces sticking mechanism for the binding of VLB to the interface of α,β-tubulin. The double-faces sticking mechanism can be used to qualitatively explain high cytotoxicity toward cancer cells of vinca alkaloids including vinblastine, vincristine, vindestine, and vinorelbine approved for clinical use and vinflunine still in a phase III clinical trial. Furthermore, this mechanism will be applied to develop novel vinca alkaloids with much higher cytotoxicity toward cancer cells.

Synthesis, structure and antibacterial activity of potent DNA gyrase inhibitors: N'-benzoyl-3-(4-bromophenyl)-1H-pyrazole-5-carbohydrazide derivatives

A total of 19 novel (3a-3s) N'-benzoyl-3-(4-bromophenyl)-1H-pyrazole-5-carbohydrazide analogs were designed, synthesized, and evaluated for biological activities as potential DNA gyrase inhibitors. The results showed that compound 3k can strongly inhibit Staphylococcus aureus DNA gyrase and Bacillus subtilis DNA gyrase (with IC50 of 0.15 µg/mL and 0.25 µg/mL, respectively). Structure-activity relationships were also discussed base on the biological and docking simulation results.

The effect of fucoidan on tyrosinase: computational molecular dynamics integrating inhibition kinetics

Fucoidan is a complex sulfated polysaccharide extracted from brown seaweed and has a wide variety of biological activities. In this study, we investigated the inhibitory effect of fucoidan on tyrosinase via a combination of inhibition kinetics and computational simulations. Fucoidan reversibly inhibited tyrosinase in a mixed-type manner. Time-interval kinetics showed that the inhibition was processed as first order with biphasic processes. For further insight, we simulated dockings with various sizes of molecular models (monomer to decamer) of fucoidan and showed that the best binding energy change results were obtained from the pentamer (-1.89 kcal/mol) and the hexamer (-1.97 kcal/mol) models of AutoDock Vina. The molecular dynamics simulation confirmed the binding mechanisms between tyrosinase and fucoidan and suggested that fucoidan mostly interacts with several residues including copper ions located in the active site. Our study suggests that fucoidan might be a potential natural antipigment agent.

Traditional Chinese medicine as dual guardians against hypertension and cancer?

This study utilizes the comprehensive traditional Chinese medicine database TCM Database@Taiwan ( http://tcm.cmu.edu.tw/ ) in conjunction with structure-based and ligand-based drug design to identify multi-function Src inhibitors. The three potential TCM candidates identified as having suitable docking conformations and bioactivity profiles were Angeliferulate, (3R)-2'-hydroxy-3',4'-dimethoxyisoflavan-7-O-beta-D-glucoside (HMID), and 3-[2',6-dihydroxy-5'-(2-propenyl)[1,1'-biphenyl]3-yl]-(E)-2-propenoic acid (3PA). Molecular dynamics simulation demonstrated that the TCM candidates have more stable interactions with the cleft and in complex with Src kinase compared to Saracatinib. Angeliferulate and HMID, both originated from Angelica sinensis, not only interact with Lys298 and amino acids from different loops in the cleft, but also with Asp407 located on the activation loop. These interactions are important to reduce the opening of the activation loop due to phosphorylation, hence stabilize the Src kinase cleft structure and inhibit activation. The TCM candidates also exhibited high affinity to other cancer-related target proteins (EGFR, HER2, and HSP90). Our observations suggest that the TCM candidates might have multi-targeting effects in hypertension and cancer.

Inhibition of mushroom tyrosinase by a newly synthesized ligand: inhibition kinetics and computational simulations

Alterations in the synthesis of melanin contribute to a number of diseases; therefore, the design of new tyrosinase inhibitors is very important. Mushroom tyrosinase (MT) is a metalloenzyme, which plays an important role in melanin biosynthesis. In this study, the inhibitory effect of a novel designed compound, i.e. 2-((1Z)-(2-(2,4-dinitrophenyl)hydrazin-1-ylidene)methyl)phenol, as a specific ligand which can bind to the copper ion of MT, has been assessed. The ligand was found to competitively inhibit both the cresolase and catecholase activities of MT, with small inhibition constants of 2.8 and 2.6 μM, respectively. Intrinsic fluorescence studies were performed to gain more information on the binding constants. Docking results indicated that the ligand binds to copper ions in the active site of MT via the OH group of the ligand. The ligand makes four hydrogen bonds with aspartic acid and one hydrogen bond with the histidine residue in the active site. Molecular dynamics results show that ligand binds to the MT via both electrostatic and hydrophobic interactions with its different parts.

Biological evaluation and interaction of a newly designed anti-cancer Pd(II) complex and human serum albumin

The pharmacokinetics and pharmacodynamics of any drug will depend, largely, on the interaction that has with human serum albumin (HSA), the most abundant plasma protein. The interaction between newly synthesized Pd(II) complexe, 2,2'-bipyridin Butylglycinato Pd(II) nitrate, an anti-tumor component, with HSA was studied at different temperatures by fluorescence, far UV circular dichroism (CD), UV-visible spectrophotometry and theoretical approaches. The Pd(II) complex has a strong ability to quench the intrinsic fluorescence of HSA through a dynamic quenching procedure. The binding parameters and thermodynamic parameters, including δH°, δS° and δG° were calculated by fluorescence quenching method, indicated that hydrophobic forces play a major role in the interaction of Pd(II) complex with HSA. Based on Autodock, FRET (fluorescence resonance energy transfer) and fluorescence quenching data, it may be concluded that one of the binding sites in the complex of HSA is near the only one Trp of HSA (Trp214) in sub domain IIA of the protein. Far-UV-CD results indicated that Pd(II)-complex induced increase in the α-helical content of the protein. The anti-tumor property of the synthesized Pd(II) complex was studied by testing it on human tumor cell line K562. The 50% cytotoxic concentration (Cc₅₀) of complex was determined using MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide) assay. Also, fluorescence staining with DAPI (4,6-diamidino-2-phenylindole) revealed some typical nuclear changes that are characteristic of apoptosis which is induced at Cc₅₀ concentration of Pd(II) complex in K562 cell line after 24 h incubation. Our results suggest that Pd(II) complex is a promising anti-proliferative agent and should execute its biological effects by inducing apoptosis.

Molecular dynamics and docking studies on cardiac troponin C

Cardiac troponin C (cTnC) is the Ca²⁺ dependent switch for contraction in heart muscle making it a potential target for drug research in the therapy of heart failure. Calcium binding on Troponin C (TnC) triggers a series of conformational changes exposing a hydrophobic pocket in the N-domain of TnC (cNTnC), which leads to force generation. Mutations and acidic pH have been related to altering the sensitivity of TnC affecting the efficiency of the heart. Bepridil, identified as a calcium sensitizer to TnC, has been experimentally found to bind to the N-domain pocket of TnC but with negative cooperativity. Screening and de novo design were carried out using LUDI and AUTOLUDI programs in this work to identify and design potential ligands that can bind to the hydrophobic pocket of TnC. Two docking centers and multiple searching radii including 5 Å, 5.5 Å, 6 Å, 6.5 Å, 7.0 Å and 7.5 Å were used in LUDI to screen the ZINC database. Based on the LUDI docking results, 8 molecules were identified from the database with good potential to bind into the binding pocket and they were used as template molecules to generate a series of new molecules by AUTOLUDI design. Out of all the newly-designed molecules, 14 new ligands were recognized to be potential ligands that can bind and fit well into the binding pocket. These molecules can be used as starting molecules to develop TnC ligands. The binding stability and binding affinity of these molecules to the protein was further analyzed by molecular dynamics simulations. The results show that the binding energies, interactions and complex stabilities of 6 ligands are comparable to or better than bepridil.

Biphasic association of T7 RNA polymerase and a nucleotide analogue, cibacron blue as a model to understand the role of initiating nucleotide in the mechanism of enzyme action

T7 RNA polymerase (T7 RNAP) is an enzyme that utilizes ribonucleotides to synthesize the nascent RNA chain in a template-dependent manner. Here we have studied the interaction of T7 RNAP with cibacron blue, an anthraquinone monochlorotriazine dye, its effect on the function of the enzyme and the probable mode of binding of the dye. We have used difference absorption spectroscopy and isothermal titration calorimetry to show that the dye binds T7 RNAP in a biphasic manner. The first phase of the binding is characterized by inactivation of the enzyme. The second binding site overlaps with the common substrate-binding site of the enzyme. We have carried out docking experiment to map the binding site of the dye in the promoter bound protein. Competitive displacement of the dye from the high affinity site by labeled GTP and isothermal titration calorimetry of high affinity GTP bound enzyme with the dye suggests a strong correlation between the high affinity dye binding and the high affinity GTP binding in T7 RNAP reported earlier from our laboratory.

In search of allosteric modulators of a7-nAChR by solvent density guided virtual screening

Nicotinic acetylcholine receptors (nAChR) are pentameric ligand gated ion channels whose activity can be modulated by endogenous neurotransmitters as well as by synthetic ligands that bind the same or distinct sites from the natural ligand. The subtype of α7 nAChR has been considered as a potenial therapeutic target for Alzheimer's disease, schizophrenia and other neurological and psychiatric disorders. Here we have developed a homology model of α7 nAChR based on two high resolution crystal structures with Brookhaven Protein Data Bank (PDB) codes 2QC1 and 2WN9 for threading on one monomer and then for building a pentamer, respectively. A number of small molecule binding sites are identified using Pocket Finder (J. An, M. Tortov, and R. Abagyan, Molecular & Cellular Proteomics, 4.6, 752-761 (2005)) of Internal Coordinate Mechanics (ICM). Remarkably, these computer-identified sites match perfectly with ordered solvent densities found in the high-resolution crystal structure of α1 nAChR, suggesting that the surface cavities in the α7 nAChR model are likely binding sites of small molecules. A high throughput virtual screening by flexible ligand docking of 5008 small molecule compounds was performed at three potential allosteric modulator (AM) binding sites of α7 nAChR using Molsoft ICM software (R. Abagyan, M. Tortov and D. Kuznetsov, J Comput Chem 15, 488-506, (1994)). Some experimentally verified allosteric modulators of α7 like CCMI comp-6, LY 7082101, 5-HI, TQS, PNU-120596, genistein, and NS-1738 ranked among top 100 compounds, while the rest of the compounds in the list could guide further search for new allosteric modulators.

The key residues of active sites on the catalytic fragment for paclitaxel interacting with poly (ADP-ribose) polymerase

Poly(ADP-ribose) polymerase (PARP) is regarded as a target protein for paclitaxel (PTX) to bind. An important issue is to identify the key residues as active sites for PTX interacting with PARP, which will help to understand the potential drug activity of PTX against cancer cells. Using docking method and MD simulation, we have constructed a refined structure of PTX docked on the catalytic function domain of PARP (PDB code: 1A26). The residues Glu327(988), Tyr246(907), Lys242(903), His165(826), Asp105(766), Gln102(763) and Gln98(759) in PARP are identified as potential sites involved in interaction with PTX according to binding energy (E(b)) between PTX and single residue calculated with B3LYP/6-31G(d,p). These residues form an active binding pocket located on the surface of the catalytic fragment, possibly interacting with the required groups of PTX leading to its activity against cancer cells. It is noted that most of the active sites make conatct with the "southern hemisphere" of PTX except for one residue, Tyr246(907), which interacts with the "northern hemisphere" of PTX. The conformation of PTX in complex with the catalytic fragment is observed as being T-shaped, similar to that complexed with β-tubulin. The total Eb of -269.9 kJ/mol represents the potent interaction between PTX and the catalytic fragment, implying that PTX can readily bind to the active pocket. The tight association of PTX with the catalytic fragment would inhibit PARP activation, suggesting a potential application of PTX as an effective antineoplastic agent.