Structural Insights into Pseudomonas aeruginosa Exotoxin A-Elongation Factor 2 Interactions: A Molecular Dynamics Study.

Exotoxin A (ETA) is an extracellular secreted toxin and a single-chain polypeptide with A and B fragments that is produced by Pseudomonas aeruginosa. It catalyzes the ADP-ribosylation of a post-translationally modified histidine (diphthamide) on eukaryotic elongation factor 2 (eEF2), which results in the inactivation of the latter and the inhibition of protein biosynthesis. Studies show that the imidazole ring of diphthamide plays an important role in the ADP-ribosylation catalyzed by the toxin. In this work, we employ different in silico molecular dynamics (MD) simulation approaches to understand the role of diphthamide versus unmodified histidine in eEF2 on the interaction with ETA. Crystal structures of the eEF2-ETA complexes with three different ligands NAD+, ADP-ribose, and ßTAD were selected and compared in the diphthamide and histidine containing systems. The study shows that NAD+ bound to ETA remains very stable in comparison with other ligands, enabling the transfer of ADP-ribose to the N3 atom of the diphthamide imidazole ring in eEF2 during ribosylation. We also show that unmodified histidine in eEF2 has a negative impact on ETA binding and is not a suitable target for the attachment of ADP-ribose. Analyzing of radius of gyration and COM distances for NAD+, ßTAD, and ADP-ribose complexes revealed that unmodified His affects the structure and destabilizes the complex with all different ligands throughout the MD simulations.

Structural insights in galectin-1-glycan recognition: Relevance of the glycosidic linkage and the N-acetylation pattern of sugar moieties.

Galectins, soluble lectins widely expressed intra- and extracellularly in different cell types, play major roles in deciphering the cellular glycocode. Galectin-1 (Gal-1), a prototype member of this family, presents a carbohydrate recognition domain (CRD) with specific affinity for ß-galactosides such as N-acetyllactosamine (ß-d-Galp-(1 â†’ 4)-d-GlcpNAc), and mediate numerous physiological and pathological processes. In this work, Gal-1 binding affinity for ß-(1 â†’ 6) galactosides, including ß-d-Galp-(1 â†’ 6)-ß-d-GlcpNAc-(1 â†’ 4)-d-GlcpNAc was evaluated, and their performance was compared to that of ß-(1 â†’ 4) and ß-(1 â†’ 3) galactosides. To this end, the trisaccharide ß-d-Galp-(1 â†’ 6)-ß-d-GlcpNAc-(1 â†’ 4)-d-GlcpNAc was enzymatically synthesized, purified and structurally characterized. To evaluate the affinity of Gal-1 for the galactosides, competitive solid phase assays (SPA) and isothermal titration calorimetry (ITC) studies were carried out. The experimental dissociation constants and binding energies obtained were compared to those calculated by molecular docking. These analyses evidenced the critical role of the glycosidic linkage between the terminal galactopyranoside residue and the adjacent monosaccharide, as galactosides bearing ß-(1 â†’ 6) glycosidic linkages showed dissociation constants six- and seven-fold higher than those involving ß-(1 â†’ 4) and ß-(1 â†’ 3) linkages, respectively. Moreover, docking experiments revealed the presence of hydrogen bond interactions between the N-acetyl group of the glucosaminopyranose moiety of the evaluated galactosides and specific amino acid residues of Gal-1, relevant for galectin-glycan affinity. Noticeably, the binding free energies (ΔGbindcalc) derived from the molecular docking were in good agreement with experimental values determined by ITC measurements (ΔGbindexp), evidencing a good correlation between theoretical and experimental approaches, which validates the in silico simulations and constitutes an important tool for the rational design of future optimized ligands.

Engineering of Ancestors as a Tool to Elucidate Structure, Mechanism, and Specificity of Extant Terpene Cyclase.

Structural information is crucial for understanding catalytic mechanisms and to guide enzyme engineering efforts of biocatalysts, such as terpene cyclases. However, low sequence similarity can impede homology modeling, and inherent protein instability presents challenges for structural studies. We hypothesized that X-ray crystallography of engineered thermostable ancestral enzymes can enable access to reliable homology models of extant biocatalysts. We have applied this concept in concert with molecular modeling and enzymatic assays to understand the structure activity relationship of spiroviolene synthase, a class I terpene cyclase, aiming to engineer its specificity. Engineering a surface patch in the reconstructed ancestor afforded a template structure for generation of a high-confidence homology model of the extant enzyme. On the basis of structural considerations, we designed and crystallized ancestral variants with single residue exchanges that exhibited tailored substrate specificity and preserved thermostability. We show how the two single amino acid alterations identified in the ancestral scaffold can be transferred to the extant enzyme, conferring a specificity switch that impacts the extant enzyme's specificity for formation of the diterpene spiroviolene over formation of sesquiterpenes hedycaryol and farnesol by up to 25-fold. This study emphasizes the value of ancestral sequence reconstruction combined with enzyme engineering as a versatile tool in chemical biology.

Enzymatic synthesis of non-natural trisaccharides and galactosides; Insights of their interaction with galectins as a function of their structure.

Galectins are a family of carbohydrate-recognizing proteins that by interacting with specific glycoepitopes can mediate important biological processes, including immune cell homeostasis and activation of tolerogenic circuits. Among the different members of this family, Galectin 1 and 3 have shown pro-tumorigenic effects, being overexpressed in numerous neoplasic diseases, proving to be relevant in tumor immune escape, tumor progression and resistance to drug-induced apoptosis. Thus, generation of specific glycosides that could inhibit their pro-tumorigenic ability by blocking their carbohydrate recognition domain is one of the current major challenges in the field. Considering that galectin-ligand binding strength is closely related to the ligand structure, analysis of this relationship provides valuable information for rational design of high-affinity ligands that could work as effective galectin inhibitors. Taking profit of the ability of glycosidases to catalyze transglycosylation reactions we achieved the enzymatic synthesis of ß-d-Galp-(1 → 6)-ß-d-Galp-(1 → 4)-d-Glcp(2), a mixture of ß-d-Galp-(1 → 6)-ß-d-Glcp-(1 → 4)-d-Glcp(5) and ß-d-Galp-(1 → 3)-ß-d-Glcp-(1 → 4)-d-Glcp(6), and finally benzyl ß-d-galactopyranoside (9), with reaction yields between 16 and 27%. All the galactosides were purified, and characterized using 1H and 13C nuclear magnetic resonance spectroscopy. Docking results performed between the synthesized compounds and human Galectin 1 (hGal-1) and human Galectin 3 (hGal-3) showed that the replacement of a glucose moiety linked to the terminal galactose with a galactose moiety, decreases the affinity for these galectins. Moreover, regarding the interglycosidic bond the most favorable ß-Gal linkage seems to be ß(1 → 4) followed by ß(1 → 3) and ß(1 → 6) for hGal-1, and ß(1 → 4) followed by ß(1 → 6) and ß(1 → 3) for hGal-3. These results were in accordance with the IC50 values obtained with in vitro solid phase inhibition assays. Therefore, docking results obtained in this work proved to be a very good approximation for predicting binding affinity of novel galactosides.

Exploring Polypharmacology in Drug Design.

Nowadays it is widely accepted that one compound can be able to hit several targets at once. This "magic shotgun" approach for drug development properly describes the mechanism of biomolecular recognition. The need to take into account the polypharmacology in structure-based drug design has led to the development of several computational tools. Here we present a computational protocol to identify promising compounds against several biological targets, a protocol known as inverse docking.

QM/MM Studies of Dph5 - A Promiscuous Methyltransferase in the Eukaryotic Biosynthetic Pathway of Diphthamide.

Eukaryotic diphthine synthase, Dph5, is a promiscuous methyltransferase that catalyzes an extraordinary N, O-tetramethylation of 2-(3-carboxy-3-aminopropyl)-l-histidine (ACP) to yield diphthine methyl ester (DTM). These are intermediates in the biosynthesis of the post-translationally modified histidine residue diphthamide (DTA), a unique and essential residue part of the eukaryotic elongation factor 2 (eEF2). Herein, the promiscuity of Saccharomyces cerevisiae Dph5 has been studied with in silico approaches, including homology modeling to provide the structure of Dph5, protein-protein docking and molecular dynamics to construct the Dph5-eEF2 complex, and quantum mechanics/molecular mechanics (QM/MM) calculations to outline a plausible mechanism. The calculations show that the methylation of ACP follows a typical SN2 mechanism, initiating with a complete methylation (trimethylation) at the N-position, followed by the single O-methylation. For each of the three N-methylation reactions, our calculations support a stepwise mechanism, which first involve proton transfer through a bridging water to a conserved aspartate residue D165, followed by a methyl transfer. Once fully methylated, the trimethyl amino group forms a weak electrostatic interaction with D165, which allows the carboxylate group of diphthine to attain the right orientation for the final methylation step to be accomplished.

Homology model of the human tRNA splicing ligase RtcB.

RtcB is an essential human tRNA ligase required for ligating the 2',3'-cyclic phosphate and 5'-hydroxyl termini of cleaved tRNA halves during tRNA splicing and XBP1 fragments during endoplasmic reticulum stress. Activation of XBP1 has been implicated in various human tumors including breast cancer. Here we present, for the first time, a homology model of human RtcB (hRtcB) in complex with manganese and covalently bound GMP built from the Pyrococcus horikoshii RtcB (bRtcB) crystal structure, PDB ID 4DWQA. The structure is analyzed in terms of stereochemical quality, folding reliability, secondary structure similarity with bRtcB, druggability of the active site binding pocket and its metal-binding microenvironment. In comparison with bRtcB, loss of a manganese-coordinating water and movement of Asn226 (Asn202 in 4DWQA) to form metal-ligand coordination, demonstrates the uniqueness of the hRtcB model. Rotation of GMP leads to the formation of an additional metal-ligand coordination (Mn-O). Umbrella sampling simulations of Mn binding in wild type and the catalytically inactive C122A mutant reveal a clear reduction of Mn binding ability in the mutant, thus explaining the loss of activity therein. Our results furthermore clearly show that the GTP binding site of the enzyme is a well-defined pocket that can be utilized as target site for in silico drug discovery.

Structural insights into human microsomal epoxide hydrolase by combined homology modeling, molecular dynamics simulations, and molecular docking calculations.

A new homology model of human microsomal epoxide hydrolase was derived based on multiple templates. The model obtained was fully evaluated, including MD simulations and ensemble-based docking, showing that the quality of the structure is better than that of only previously known model. Particularly, a catalytic triad was clearly identified, in agreement with the experimental information available. Analysis of intermediates in the enzymatic mechanism led to the identification of key residues for substrate binding, stereoselectivity, and intermediate stabilization during the reaction. In particular, we have confirmed the role of the oxyanion hole and the conserved motif (HGXP) in epoxide hydrolases, in excellent agreement with known experimental and computational data on similar systems. The model obtained is the first one that fully agrees with all the experimental observations on the system. Proteins 2017; 85:720-730. © 2016 Wiley Periodicals, Inc.

Ligand Selectivity between the ADP-Ribosylating Toxins: An Inverse-Docking Study for Multitarget Drug Discovery.

Bacterial adenosine 5'-diphosphate-ribosylating toxins are encoded by several human pathogens, such as Pseudomonas aeruginosa (exotoxin A (ETA)), Corynebacterium diphtheriae (diphtheria toxin (DT)), and Vibrio cholerae (cholix toxin (CT)). The toxins modify eukaryotic elongation factor 2, an essential human enzyme in protein synthesis, thereby causing cell death. Targeting external virulence factors, such as the above toxins, is a promising alternative for developing new antibiotics, while at the same time avoiding drug resistance. This study aims to establish a reliable computational methodology to find a "silver bullet" able to target all three toxins. Herein, we have undertaken a detailed analysis of the active sites of ETA, DT, and CT, followed by the determination of the most appropriate selection of the size of the docking sphere. Thereafter, we tested two different approaches for normalizing the docking scores and used these to verify the best target (toxin) for each ligand. The results indicate that the methodology is suitable for identifying selective as well as multitoxin inhibitors, further validating the robustness of inverse docking for target-fishing experiments.

Structures, Properties, and Dynamics of Intermediates in eEF2-Diphthamide Biosynthesis.

The eukaryotic translation Elongation Factor 2 (eEF2) is an essential enzyme in protein synthesis. eEF2 contains a unique modification of a histidine (His699 in yeast; HIS) into diphthamide (DTA), obtained via 3-amino-3-carboxypropyl (ACP) and diphthine (DTI) intermediates in the biosynthetic pathway. This essential and unique modification is also vulnerable, in that it can be efficiently targeted by NAD(+)-dependent ADP-ribosylase toxins, such as diphtheria toxin (DT). However, none of the intermediates in the biosynthesis path is equally vulnerable against the toxins. This study aims to address the different susceptibility of DTA and its precursors against bacterial toxins. We have herein undertaken a detailed in silico study of the structural features and dynamic motion of different His699 intermediates along the diphthamide synthesis pathway (HIS, ACP, DTI, DTA). The study points out that DTA forms a strong hydrogen bond with an asparagine which might explain the ADP-ribosylation mechanism caused by the diphtheria toxin (DT). Finally, in silico mutagenesis studies were performed on the DTA modified protein, in order to hamper the formation of such a hydrogen bond. The results indicate that the mutant structure might in fact be less susceptible to attack by DT and thereby behave similarly to DTI in this respect.

Improved Homology Model of the Human all-trans Retinoic Acid Metabolizing Enzyme CYP26A1.

A new CYP26A1 homology model was built based on the crystal structure of cyanobacterial CYP120A1. The model quality was examined for stereochemical accuracy, folding reliability, and absolute quality using a variety of different bioinformatics tools. Furthermore, the docking capabilities of the model were assessed by docking of the natural substrate all-trans-retinoic acid (atRA), and a group of known azole- and tetralone-based CYP26A1 inhibitors. The preferred binding pose of atRA suggests the (4S)-OH-atRA metabolite production, in agreement with recently available experimental data. The distances between the ligands and the heme group iron of the enzyme are in agreement with corresponding distances obtained for substrates and azole inhibitors for other cytochrome systems. The calculated theoretical binding energies agree with recently reported experimental data and show that the model is capable of discriminating between natural substrate, strong inhibitors (R116010 and R115866), and weak inhibitors (liarozole, fluconazole, tetralone derivatives).

Improved homology model of cyclohexanone monooxygenase from Acinetobacter calcoaceticus based on multiple templates.

A new homology model of cyclohexanone monooxygenase (CHMO) from Acinetobacter calcoaceticus is derived based on multiple templates, and in particular the crystal structure of CHMO from Rhodococcus sp. The derived model was fully evaluated, showing that the quality of the new structure was improved over previous models. Critically, the nicotinamide cofactor is included in the model for the first time. Analysis of several molecular dynamics snapshots of intermediates in the enzymatic mechanism led to a description of key residues for cofactor binding and intermediate stabilization during the reaction, in particular Arg327 and the well known conserved motif (FxGxxxHxxxW) in Baeyer-Villiger monooxygenases, in excellent agreement with known experimental and computational data.

Homology models of human all-trans retinoic acid metabolizing enzymes CYP26B1 and CYP26B1 spliced variant.

Homology models of CYP26B1 (cytochrome P450RAI2) and CYP26B1 spliced variant were derived using the crystal structure of cyanobacterial CYP120A1 as template for the model building. The quality of the homology models generated were carefully evaluated, and the natural substrate all-trans-retinoic acid (atRA), several tetralone-derived retinoic acid metabolizing blocking agents (RAMBAs), and a well-known potent inhibitor of CYP26B1 (R115866) were docked into the homology model of full-length cytochrome P450 26B1. The results show that in the model of the full-length CYP26B1, the protein is capable of distinguishing between the natural substrate (atRA), R115866, and the tetralone derivatives. The spliced variant of CYP26B1 model displays a reduced affinity for atRA compared to the full-length enzyme, in accordance with recently described experimental information.

Cloning and functional studies of a splice variant of CYP26B1 expressed in vascular cells.

BACKGROUND: All-trans retinoic acid (atRA) plays an essential role in the regulation of gene expression, cell growth and differentiation and is also important for normal cardiovascular development but may in turn be involved in cardiovascular diseases, i.e. atherosclerosis and restenosis. The cellular atRA levels are under strict control involving several cytochromes P450 isoforms (CYPs). CYP26 may be the most important regulator of atRA catabolism in vascular cells. The present study describes the molecular cloning, characterization and function of atRA-induced expression of a spliced variant of the CYP26B1 gene. METHODOLOGY/PRINCIPAL FINDINGS: The coding region of the spliced CYP26B1 lacking exon 2 was amplified from cDNA synthesized from atRA-treated human aortic smooth muscle cells and sequenced. Both the spliced variant and full length CYP26B1 was found to be expressed in cultured human endothelial and smooth muscle cells, and in normal and atherosclerotic vessel. atRA induced both variants of CYP26B1 in cultured vascular cells. Furthermore, the levels of spliced mRNA transcript were 4.5 times higher in the atherosclerotic lesion compared to normal arteries and the expression in the lesions was increased 20-fold upon atRA treatment. The spliced CYP26B1 still has the capability to degrade atRA, but at an initial rate one-third that of the corresponding full length enzyme. Transfection of COS-1 and THP-1 cells with the CYP26B1 spliced variant indicated either an increase or a decrease in the catabolism of atRA, probably depending on the expression of other atRA catabolizing enzymes in the cells. CONCLUSIONS/SIGNIFICANCE: Vascular cells express the spliced variant of CYP26B1 lacking exon 2 and it is also increased in atherosclerotic lesions. The spliced variant displays a slower and reduced degradation of atRA as compared to the full-length enzyme. Further studies are needed, however, to clarify the substrate specificity and role of the CYP26B1 splice variant in health and disease.

COMPARE analysis of the toxicity of an iminoquinone derivative of the imidazo[5,4-f]benzimidazoles with NAD(P)H:quinone oxidoreductase 1 (NQO1) activity and computational docking of quinones as NQO1 substrates.

Synthesis and cytotoxicity of imidazo[5,4-f]benzimidazolequinones and iminoquinone derivatives is described, enabling structure-activity relationships to be obtained. The most promising compound (an iminoquinone derivative) has undergone National Cancer Institute (NCI) 60 cell line (single and five dose) screening, and using the NCI COMPARE program, has shown correlation to NQO1 activity and to other NQO1 substrates. Common structural features suggest that the iminoquinone moiety is significant with regard to NQO1 specificity. Computational docking into the active site of NQO1 was performed, and the first comprehensive mitomycin C (MMC)-NQO1 docking study is presented. Small distances for hydride reduction and high binding affinities are characteristic of MMC and of iminoquinones showing correlations with NQO1 via COMPARE analysis. Docking also indicated that the presence of a substituent capable of hydrogen bonding to the His194 residue is important in influencing the orientation of the substrate in the NQO1 active site, leading to more efficient reduction.

Properties and behaviour of tetracyclic allopsoralen derivatives inside a DPPC lipid bilayer model.

Allopsoralens are angular psoralen derivatives presenting advantages over the parent compound because of monofunctional DNA-photobinding and consequent lower toxicity. Allopsoralen molecules with three different substituents and different protonation states were studied using the molecular dynamics technique. The location of these molecules when inside the lipid bilayer is of major importance because their photochemical properties can change with the environment. Also, the ability of psoralens to form photoadducts with unsaturated phospholipids depends on the preference of the molecules to locate themselves closer to the bilayer middle were the double bond functionality can be found. Herein we show that the allopsoralens tend to accumulate inside the lipid bilayer closer to the water interface when protonated or closer to the interface middle otherwise. Allopsoralens containing one amine terminated carbon chain tend to have different rotational and orientational behaviour and an orientation preference close to the ones shown by the lipids. The size and chemical nature of the substituent also affect the molecular mobility and capacity to interact with water molecules and the nitrogen or phosphorus atoms of the lipids.

Mechanism of the organocatalyzed decarboxylative Knoevenagel-Doebner reaction. A theoretical study.

We have investigated important intermediates and key transition states of the organocatalyzed Knoevenagel condensation using density functional theory and two different basis sets (6-31 G(d,p) and 6-311++G(2df,2pd)), both in gas phase and simulating the bulk solvent (pyridine) using the PCM method. Calculated structures for reactants, intermediates, and key transition states suggest that the secondary amine catalyst is essential, both for activating the aldehyde for nucleophilic attack, and in the possible decarboxylation pathways. The calculated results are shown to agree with available experimental information. On the basis of the results obtained, the studied mechanism may be important in the understanding of vinylphenol production during malting and brewing of wheat and barley grains.