Construction and assessment of reaction models of class I EPSP synthase: molecular docking and density functional theoretical calculations

Virtual screening of natural products against 5-enolpyruvylshikimate-3-phosphate synthase using the Anagreen herbicide-like natural compound library

The shikimate pathway enzyme 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) catalyzes a reaction involved in the production of amino acids essential for plant growth and survival. EPSPS is the main target of glyphosate, a broad-spectrum herbicide that acts as a competitive inhibitor concerning phosphoenolpyruvate (PEP), which is the natural substrate of EPSPS. In the present study, we introduce a natural compound library, named Anagreen, which is a compendium of herbicide-like compounds obtained from different natural product databases. Herein, we combined the structure- and ligand-based virtual screening strategies to explore Anagreen against EPSPS using the structure of glyphosate complexed with a T102I/P106S mutant of EPSPS from Eleusine indica (EiEPSPS) as a starting point. First, ligand-based pharmacophore screening was performed to select compounds with a similar pharmacophore to glyphosate. Then, structure-based pharmacophore modeling was applied to build a model which represents the molecular features of glyphosate. Then, consensus docking was performed to rank the best poses of the natural compounds against the PEP binding site, and then molecular dynamics simulations were performed to analyze the stability of EPSPS complexed with the selected ligands. Finally, we have investigated the binding affinity of the complexes using free energy calculations. The selected hit compound, namely AG332841, showed a stable conformation and binding affinity to the EPSPS structure and showed no structural similarity to the already known weed EPSPS inhibitors. Our computational study aims to clarify the inhibition of the mutant EiEPSPS, which is resistant to glyphosate, and identify new potential herbicides from natural products.

In Silico Pesticide Discovery for New Anti-Tobacco Mosaic Virus Agents: Reactivity, Molecular Docking, and Molecular Dynamics Simulations

Considerable data are available regarding the molecular genetics of the tobacco mosaic virus. The disease caused by the tobacco mosaic virus is still out of control due to the lack of an efficient functional antagonist chemical molecule. Extensive research was carried out to try to find effective new anti-tobacco mosaic virus agents, however no study could find an effective agent which could completely inhibit the disease caused by the virus. In recent years, molecular docking, combined with molecular dynamics, which is considered to be one of the most important methods of drug discovery and design, were used to evaluate the type of binding between the ligand and its protein enzyme. The aim of the current work was to assess the in silico anti-tobacco mosaic virus activity for a selection of 41 new and 2 reference standard compounds. These compounds were chosen to examine their reactivity and binding efficiency with the tobacco mosaic virus coat protein (PDB ID: 2OM3). A comparison was made between the activity of the selected compounds and that for ningnanmycin and ribavirin, which are common inhibitors of plant viruses. The simulation results obtained from the molecular docking and molecular dynamics showed that two compounds of the antofine analogues could bind with the tobacco mosaic virus coat protein receptor better than ningnanmycin and ribavirin.

Molecular docking, molecular dynamics simulations and reactivity, studies on approved drugs library targeting ACE2 and SARS-CoV-2 binding with ACE2

The recent new contagion coronavirus 2019 (COVID-19) disease is a new generation of severe acute respiratory syndrome coronavirus-2 SARS-CoV-2 which infected millions confirmed cases and hundreds of thousands death cases around the world so far. Molecular docking combined with molecular dynamics is one of the most important tools of drug discovery and drug design, which it used to examine the type of binding between the ligand and its protein enzyme. Global reactivity has important properties, which enable chemists to understand the chemical reactivity and kinetic stability of compounds. In this study, molecular docking and reactivity were applied for eighteen drugs, which are similar in structure to chloroquine and hydroxychloroquine, the potential inhibitors to angiotensin-converting enzyme (ACE2). Those drugs were selected from DrugBank. The reactivity, molecular docking and molecular dynamics were performed for two receptors ACE2 and [SARS-CoV-2/ACE2] complex receptor in two active sites to find a ligand, which may inhibit COVID-19. The results obtained from this study showed that Ramipril, Delapril and Lisinopril could bind with ACE2 receptor and [SARS-CoV-2/ACE2] complex better than chloroquine and hydroxychloroquine. This new understanding should help to improve predictions of the impact of such alternatives on COVID-19.Communicated by Ramaswamy H. Sarma.

Review: amino acid biosynthesis as a target for herbicide development

There are three amino acid biosynthesis pathways that are targeted by current herbicides, namely those leading to the production of aromatic amino acids, branched chain amino acids and glutamine. However, their efficacy is diminishing as a result of the increasing number of resistant weeds. Indeed, resistance to most classes of herbicides is on the rise, posing a significant threat to the utility of current herbicides to sustain effective weed management. This review provides an overview of potential herbicide targets within amino acid biosynthesis that remain unexploited commercially, and recent inhibitor discovery efforts. Despite contemporary approaches to herbicide discovery, such as chemical repurposing and the use of omics technologies, there have been no new products introduced to the market that inhibit amino acid biosynthesis over the past three decades. This highlights the chasm that exists between identifying a potent inhibitor and introducing a commercial herbicide. The unpredictability of a mode of action at the systemic level, as well as poor physicochemical properties, often contribute to a lack of progression beyond the target inhibition stage. Nevertheless, it will be important to overcome these obstacles for the development of new herbicides to protect our agricultural industry and ensure food security for an increasing world population. © 2020 Society of Chemical Industry.

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.

Investigation of the target-site resistance of EPSP synthase mutants P106T and T102I/P106S against glyphosate

The shikimate pathway enzyme 5-enolpyruvyl shikimate-3-phosphate synthase (EPSPS) catalyzes the reaction involved in the production of amino acids essential for plant growth and survival. Thus, EPSPS is the main target of various herbicides, including glyphosate, a broad-spectrum herbicide that acts as a competitive inhibitor of phosphoenolpyruvate (PEP), which is the natural substrate of EPSPS. However, punctual mutations in the EPSPS gene have led to glyphosate resistance in some plants. Here, we investigated the mechanism of EPSPS resistance to glyphosate in mutants of two weed species, Conyza sumatrensis (mutant, P106T) and Eleusine indica (mutant, T102I/P106S), both of which have an economic impact on industrial crops. Molecular dynamics (MD) simulations and binding free energy calculations revealed the influence of the mutations on the affinity of glyphosate in the PEP-binding site. The amino acid residues of the EPSPS protein in both species involved in glyphosate resistance were elucidated as well as other residues that could be useful for protein engineering. In addition, during MD simulations, we identified conformational changes in glyphosate when complexed with resistant EPSPS, related to loss of herbicide activity and binding affinity. Our computational findings are consistent with previous experimental results and clarify the inhibitory activity of glyphosate as well as the structural target-site resistance of EPSPS against glyphosate.

Single or double EPSP synthase mutations lead glyphosate to undergo conformational changes that limit its inhibitory action.

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.

Resistance determination of the ACCase-inhibiting herbicide of clodinafop propargyl in Avena ludoviciana (Durieu), and study of their interaction using molecular docking and simulation

Structural mutations providing herbicide resistance may cause a modification of the three dimensional structure of a protein which will lead to a decrease in the herbicide efficacy. Wild oat (Avena ludoviciana Durieu.) is an increasingly disruptive weed in areas of intensive cereal production, thus the aim of this research was to identify mutations conferring resistance to ACCase-inhibitor herbicides at greenhouse, laboratory and in silico scales. Among the selected biotypes, No. 3 in the position 1781 (Ile1781-Leu) and No. 14 in the position 2041 (Ile2041-Asn), showed resistance to ACCase-inhibitor. The above mutations were confirmed using the specific primers and PCR-based methods. Analysis of molecular docking indicated that residues of Trp1948 and Pro2001 are important in the binding site and showed remarkable variation in the mutation types. Using molecular dynamic simulation analysis, we demonstrated that mutation types changed the conformation of the enzyme. These changes resulted in compressed conformation in the active site, which limited the availability of binding herbicide-enzyme. In present, no crystallography molecular structure and modeling reported on the ACCase of plants and this study investigated interactions of clodinafop propargyl and ACCase CT domain in A. ludoviciana by modeling, docking and simulations for the first time. Totally, bioinformatics analysis as well as PCR-based method confirmed that herbicide resistance conferred by nucleotide mutations in the gene sequence.

Biotransformation of disperse dyes using nitroreductase immobilized on magnetic particles modified with tosyl group: Identification of products by LC-MS-MS and theoretical studies conducted with DNA

The present work evaluates the action of nitroreductase enzyme immobilized on Tosylactivated magnetic particles (MP-Tosyl) on three disperse dyes which contain nitro and azo groups. The dyes included Disperse Red 73 (DR 73), Disperse Red 78 (DR 78), and Disperse Red 167 (DR 167). The use of a magnet enabled the rapid and easy removal of the immobilized enzyme after biotransformation; this facilitated the identification of the products generated using high-performance liquid chromatography with diode array detector (HPLC-DAD) and mass spectrometry (LC-MS/MS). The main products formed by the in vitro biotransformation were identified as the product of nitro group reduction to the correspondent amine groups, which were denoted as follows: 50% of 2-(2-(4-((2-cyanoethyl)(ethyl)amino)phenyl)hydrazinyl)-5-nitrobenzonitrile, 98% of 3-((4-((4-amino-2-chlorophenyl) diazenyl)phenyl) (ethyl)amino)propanenitrile and 99% of (3-acetamido-4 - ((4-amino-2-chlorophenyl) diazenyl) phenyl) azanediyl) bis (ethane-2,1-diyl) for DR 73, DR 78 and DR 167, respectively. Based on the docking studies, the dyes investigated were found to be biotransformed by nitroreductase enzyme due to their favorable interaction with the active site of the enzyme. Theoretical results show that DR73 dye exhibits a relatively lower rate of degradation; this is attributed to the cyanide substituent which affects the electron density of the azo group. The docking studies also indicate that all the dyes presented significant reactivity towards DNA. However, Disperse Red 73 was found to exhibit a substantially higher reactivity compared to the other dyes; this implies that the dye possesses a relatively higher mutagenic power. The docking results also show that DR 73, DR 78 and DR 167 may be harmful to both humans and the environment, since the mutagenicity of nitro compounds is associated with the products formed during the reduction of nitro groups. These products can interact with biomolecules, including DNA, causing toxic and mutagenic effects.

On the degradation pathway of glyphosate and glycine

The degradation in water of the most widespread herbicide, glyphosate, is still under debate. Experimental disagreements on this process exist and there are only a few theoretical studies to support any conclusions. Moreover, the relationship between glyphosate and glycine is underestimated. Besides the structural similarity, glycine is a product of glyphosate degradation; hence, their studies are complementary. In this study, two mechanisms for the decomposition of the glyphosate molecule and glycine molecule in water are proposed. These mechanisms were explored by using quantum mechanical calculations. A combined microsolvation/PCM approach was employed to find and characterize their transition states, by which the reaction pathways were determined via the IRC method. The results have shown that the degradation processes might occur via a C-C bond cleavage, through a concerted mechanism, whereby the proton transfers and the CO2 detachments occur simultaneously. The second mechanism had two consecutive steps, a decarboxylation followed by the proton transfers. The water molecules served as a conduit for the proton transfers, away from the amine group (or the phosphonate, glyphosate case). Their function was to assist the reactions in a water-mediated decarboxylation. In these particular cases, the free energy of activation was 42.68 and 42.28 kcal mol-1 for the glycine structure and the glyphosate structure, respectively. These results agreed with the photodegradation and thermodegradation of glyphosate, as well as with the spontaneous decarboxylation of glycine. A concerted mechanism might be expected to yield C-P and C-N bond cleavages in the glyphosate molecule.

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.

Mechanistic studies of new oximes reactivators of human butyryl cholinesterase inhibited by cyclosarin and sarin

Butyryl cholinesterase (BChE) has been seen as a key enzyme in the search for new strategies in the treatment of poisoning by organophosphates (OPs), since human BChE (HssBChE), complexed with the appropriate oxime, can be a suitable scavenger and deactivator for OPs in the blood stream. However, the efficacy of HssBChE is limited by its strict stoichiometric scavenging, slow reactivation, and propensity for aging. The improvement of the reactivation rate by new and more efficient oximes could contribute to mitigate this problem and increase the HssBChE efficiency as scavenger. Several oximes have been synthesized and tested with this goal, some with promising results, but the mechanistic aspects of the reactivation reaction are not fully understood yet. In order to better investigate this mechanism, docking and mixed quantum and molecular mechanics combined with principal components analysis were performed here to evaluate the capacity of reactivation and determine the preferred route for the reactivation reaction of two new oximes on HssBChE inhibited by the neurotoxic agents cyclosarin and sarin. Plots of potential energies were calculated and all the transition states of the reactional mechanism were determined. Our results showed a good correlation with experimental data and pointed to the most efficient oxime with both OPs. The protocol used could be a suitable tool for a preliminary evaluation of the HssBChE reactivation rates by new oximes.

Peculiarities of copper binding to alpha-synuclein

Heavy metals have been implicated as the causative agents for the pathogenesis of the most prevalent neurodegenerative disease. Various mechanisms have been proposed to explain the toxic effects of metals ranging from metal-induced oxidation of protein to metal-induced changes in the protein conformation. Aggregation of a-synuclein is implicated in Parkinson's disease (PD), and various metals, including copper, constitute a prominent group of alpha-synuclein aggregation enhancers. In this study, we have systematically characterized the a-synuclein-Cu21 binding sites and analyzed the possible role of metal binding in a-synuclein fibrillation using a set of biophysical techniques, such as electron paramagnetic resonance (EPR), electron spin-echo envelope modulation (ESEEM), circular dichroism (CD), and size exclusion chromatography (SEC). Our analyses indicated that a-synuclein possesses at least two binding sites for Cu21. We have been able to locate one of the binding sites in the N-terminal region. Furthermore, based on the EPR studies of model peptides and Beta-synuclein, we concluded that the suspected His residue did not appear to participate in strong Cu21 binding.

Molecular modeling and in vitro reactivation study between the oxime BI-6 and acetylcholinesterase inhibited by different nerve agents

Nerve agents are organophosphates acting as potent inhibitors of acetylcholinesterase (AChE), the enzyme responsible for the hydrolysis of acetylcholine and, consequently, the termination of the transmission of nerve impulses. The inhibition of AChE by an organophosphate can be reversed by a nucleophilic agent able to dephosphorylate a serine residue in the active site of AChE. In this sense, the oximes are compounds capable of removing the nerve agent and reactivate the enzyme. Here, we have applied a methodology involving theoretical docking and Quantum Mechanics/Molecular Mechanics, using the softwares Molegro(®) and Spartan(®), to evaluate the kinetic constants of reactivation and the interactions of the oxime BI-6 with AChE inhibited by different organophosphorus compounds in comparison to in vitro data. Results confirm that this method is suitable for the prediction of kinetic and thermodynamic parameters of oximes, which may be useful in the design and selection of new and more effective oximes.

Construction and assessment of reaction models of Class I EPSP synthase. Part II: investigation of the EPSP ketal

Although the proposed mechanisms are reasonable, there are still many questions about the 5-enolpyruvyl shikimate-3-phosphate (EPSP) synthase mechanism that are difficult to answer by experimental means alone. EPSP synthase is a key enzyme in the shikimic acid pathway, which is found only in plants and some micro-organisms and is also molecular target of glyphosate, active component of one of the top-selling herbicides. In the study of reaction mechanism of EPSP synthase, in addition to inorganic phosphate and EPSP products, after long time at equilibrium, it was shown that a side product is formed, the EPSP ketal. In this line, studies using density functional theory (DFT) techniques were performed to investigate the reaction mechanism of formation of EPSP and the corresponding ketal. Our findings indicate some key amino acid residues in the EPSP synthase mechanism and a possible route for the formation of the EPSP ketal.

Mechanism of the activation step of the aminoacylation reaction: a significant difference between class I and class II synthetases

In the present work we report, for the first time, a novel difference in the molecular mechanism of the activation step of aminoacylation reaction between the class I and class II aminoacyl tRNA synthetases (aaRSs). The observed difference is in the mode of nucleophilic attack by the oxygen atom of the carboxylic group of the substrate amino acid (AA) to the αP atom of adenosine triphosphate (ATP). The syn oxygen atom of the carboxylic group attacks the α-phosphorous atom (αP) of ATP in all class I aaRSs (except TrpRS) investigated, while the anti oxygen atom attacks in the case of class II aaRSs. The class I aaRSs investigated are GluRS, GlnRS, TyrRS, TrpRS, LeuRS, ValRS, IleRS, CysRS, and MetRS and class II aaRSs investigated are HisRS, LysRS, ProRS, AspRS, AsnRS, AlaRS, GlyRS, PheRS, and ThrRS. The variation of the electron density at bond critical points as a function of the conformation of the attacking oxygen atom measured by the dihedral angle ψ (C(α)-C') conclusively proves this. The result shows that the strength of the interaction of syn oxygen and αP is stronger than the interaction with the anti oxygen for class I aaRSs. This indicates that the syn oxygen is the most probable candidate for the nucleophilic attack in class I aaRSs. The result is further supported by the computation of the variation of the nonbonded interaction energies between αP atom and anti oxygen as well as syn oxygen in class I and II aaRSs, respectively. The difference in mechanism is explained based on the analysis of the electrostatic potential of the AA and ATP which shows that the relative arrangement of the ATP with respect to the AA is opposite in class I and class II aaRSs, which is correlated with the organization of the active site in respective aaRSs. A comparative study of the reaction mechanisms of the activation step in a class I aaRS (Glutaminyl tRNA synthetase) and in a class II aaRS (Histidyl tRNA synthetase) is carried out by the transition state analysis. The atoms in molecule analysis of the interaction between active site residues or ions and substrates are carried out in the reactant state and the transition state. The result shows that the observed novel difference in the mechanism is correlated with the organizations of the active sites of the respective aaRSs. The result has implication in understanding the experimentally observed different modes of tRNA binding in the two classes of aaRSs.

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.

Potent inhibitor design against H1N1 swine influenza: structure-based and molecular dynamics analysis for M2 inhibitors from traditional Chinese medicine database

The rapid spread of influenza virus subtype H1N1 poses a great threat to million lives worldwide. To search for new anti-influenza compounds, we performed molecular docking and molecular dynamics simulation to identify potential traditional Chinese medicine (TCM) constituents that could block influenza M2 channel activity. Quinic acid, genipin, syringic acid, cucurbitine, fagarine, and methyl isoferulate all have extremely well docking results as compared to control amantadine. Further de novo drug design suggests that derivatives of genipin and methyl isoferulate could have enhanced binding affinity towards M2 channel. Selected molecular dynamics simulations of M2-derivative complexes show stable hydrogen bond interactions between the derivatives and M2 residues, Ser10 and Ala9. To our best knowledge, this is the first study on the anti-viral activity of the above listed TCM compounds.

Anti-infectious bronchitis virus (IBV) activity of 1,8-cineole: effect on nucleocapsid (N) protein

In the present study, anti-IBV (infectious bronchitis virus) activity of 1,8-cineole was studied by MTT assay, as well as docking and molecular dynamic (MD) simulations. The CC50 of 1,8-cineole was above 10 mM. And the maximum noncytotoxic concentration (TD0) of 1,8-cineole was determined to be 3.90 ± 0.22 mM, which was much higher than that of ribavirin (0.78 ± 0.15 mM). 1,8-cineole could inhibit IBV with an IC(50) of 0.61 mM. MTT assay showed that the inhibition of IBV by 1, 8-cineole appears to occur moderately before entering the cell but much strongly after penetration of the virus into the cell. In silico simulations indicated that the binding site of 1,8-cineole was located at the N terminus of phosphorylated nucleocapsid (N) protein, with interaction energy equaling -40.33 kcal mol(-1). The residues TyrA92, ProA134, PheA137, AspA138 and TyrA140 had important roles during the binding process and are fully or partially conserved in various IBV strains. Based on spatial and energetic criteria, 1,8-cineole interfered with the binding between RNA and IBV N-protein. Results presented here may suggest that 1,8-cineole possesses anti-IBV properties, and therefore is a potential source of anti-IBV ingredients for the pharmaceutical industry.