Looking over Toxin-K+ Channel Interactions. Clues from the Structural and Functional Characterization of alpha-KTx Toxin Tc32, a Kv1.3 Channel Blocker

alpha-KTx toxin Tc32, from the Amazonian scorpion Tityus cambridgei, lacks the dyad motif; including Lys27, characteristic of the family and generally associated with channel blockage. The toxin has been cloned and expressed for the first time. Electrophysiological experiments, by showing that the recombinant form blocks Kv1.3 channels of olfactory bulb periglomerular cells like the natural Tc32 toxin, when tested on the Kv1.3 channel of human T lymphocytes, confirmed it is in an active fold. The nuclear magnetic resonance-derived structure revealed it exhibits an alpha/beta scaffold typical of the members of the alpha-KTx family. TdK2 and TdK3, all belonging to the same alpha-KTx 18 subfamily, share significant sequence identity with Tc32 but diverse selectivity and affinity for Kv1.3 and Kv1.1 channels. To gain insight into the structural features that may justify those differences, we used the recombinant Tc32 nuclear magnetic resonance-derived structure to model the other two toxins, for which no experimental structure is available. Their interaction with Kv1.3 and Kv1.1 has been investigated by means of docking simulations. The results suggest that differences in the electrostatic features of the toxins and channels, in their contact surfaces, and in their total dipole moment orientations govern the affinity and selectivity of toxins. In addition, we found that, regardless of whether the dyad motif is present, it is always a Lys side chain that physically blocks the channels, irrespective of its position in the toxin sequence.

Beta-caryophyllene is a dietary cannabinoid

The psychoactive cannabinoids from Cannabis sativa L. and the arachidonic acid-derived endocannabinoids are nonselective natural ligands for cannabinoid receptor type 1 (CB(1)) and CB(2) receptors. Although the CB(1) receptor is responsible for the psychomodulatory effects, activation of the CB(2) receptor is a potential therapeutic strategy for the treatment of inflammation, pain, atherosclerosis, and osteoporosis. Here, we report that the widespread plant volatile (E)-beta-caryophyllene [(E)-BCP] selectively binds to the CB(2) receptor (K(i) = 155 +/- 4 nM) and that it is a functional CB(2) agonist. Intriguingly, (E)-BCP is a common constituent of the essential oils of numerous spice and food plants and a major component in Cannabis. Molecular docking simulations have identified a putative binding site of (E)-BCP in the CB(2) receptor, showing ligand pi-pi stacking interactions with residues F117 and W258. Upon binding to the CB(2) receptor, (E)-BCP inhibits adenylate cylcase, leads to intracellular calcium transients and weakly activates the mitogen-activated kinases Erk1/2 and p38 in primary human monocytes. (E)-BCP (500 nM) inhibits lipopolysaccharide (LPS)-induced proinflammatory cytokine expression in peripheral blood and attenuates LPS-stimulated Erk1/2 and JNK1/2 phosphorylation in monocytes. Furthermore, peroral (E)-BCP at 5 mg/kg strongly reduces the carrageenan-induced inflammatory response in wild-type mice but not in mice lacking CB(2) receptors, providing evidence that this natural product exerts cannabimimetic effects in vivo. These results identify (E)-BCP as a functional nonpsychoactive CB(2) receptor ligand in foodstuff and as a macrocyclic antiinflammatory cannabinoid in Cannabis.

Unraveling principles of lead discovery: from unfrustrated energy landscapes to novel molecular anchors.

The search for novel leads is a critical step in the drug discovery process. Computational approaches to identify new lead molecules have focused on discovering complete ligands by evaluating the binding affinity of a large number of candidates, a task of considerable complexity. A new computational method is introduced in this work based on the premise that the primary molecular recognition event in the protein binding site may be accomplished by small core fragments that serve as molecular anchors, providing a structurally stable platform that can be subsequently tailored into complete ligands. To fulfill its role, we show that an effective molecular anchor must meet both the thermodynamic requirement of relative energetic stability of a single binding mode and its consistent kinetic accessibility, which may be measured by the structural consensus of multiple docking simulations. From a large number of candidates, this technique is able to identify known core fragments responsible for primary recognition by the FK506 binding protein (FKBP-12), along with a diverse repertoire of novel molecular cores. By contrast, absolute energetic criteria for selecting molecular anchors are found to be promiscuous. A relationship between a minimum frustration principle of binding energy landscapes and receptor-specific molecular anchors in their role as "recognition nuclei" is established, thereby unraveling a mechanism of lead discovery and providing a practical route to receptor-biased computational combinatorial chemistry.

Determination of alpha-conotoxin binding modes on neuronal nicotinic acetylcholine receptors

alpha-Conotoxins, from cone snails, and alpha-neurotoxins, from snakes, are competitive inhibitors of nicotinic acetylcholine receptors (nAChRs) that have overlapping binding sites in the ACh binding pocket. These disulphide-rich peptides are used extensively as tools to localize and pharmacologically characterize specific nAChRs subtypes. Recently, a homology model based on the high-resolution structure of an ACh binding protein (AChBP) allowed the three-fingered alpha-neurotoxins to be docked onto the alpha7 nAChR. To investigate if alpha-conotoxins interact with the nAChR in a similar manner, we built homology models of human alpha7 and alpha3beta2 nAChRs, and performed docking simulations of alpha-conotoxins ImI, PnIB, PnIA and MII using the program GOLD. Docking revealed that alpha-conotoxins have a different mode of interaction compared with alpha-neurotoxins, with surprisingly few nAChR residues in common between their overlapping binding sites. These docking experiments show that Imi and PnIB bind to the ACh binding pocket via a small cavity located above the beta9/beta10 hairpin of the (+)alpha7 nAChR subunit. Interestingly, PnIB, PnIA and MII were found to bind in a similar location on alpha7 or alpha3beta2 receptors mostly through hydrophobic interactions, while ImI bound further from the ACh binding pocket, mostly through electrostatic interactions. These findings, which distinguish alpha-conotoxin and alpha-neurotoxin binding modes, have implications for the rational design of selective nAChR antagonists. Copyright (C) 2004 John Wiley Sons, Ltd.

beta 2 subunit contribution to 4/7 alpha-conotoxin binding to the nicotinic acetylcholine receptor

The structures of acetylcholine-binding protein ( AChBP) and nicotinic acetylcholine receptor ( nAChR) homology models have been used to interpret data from mutagenesis experiments at the nAChR. However, little is known about AChBP-derived structures as predictive tools. Molecular surface analysis of nAChR models has revealed a conserved cleft as the likely binding site for the 4/7 alpha-conotoxins. Here, we used an alpha 3 beta 2 model to identify beta 2 subunit residues in this cleft and investigated their influence on the binding of alpha-conotoxins MII, PnIA, and GID to the alpha 3 beta 2 nAChR by two-electrode voltage clamp analysis. Although a beta 2-L119Q mutation strongly reduced the affinity of all three alpha-conotoxins, beta 2-F117A, beta 2-V109A, and beta 2-V109G mutations selectively enhanced the binding of MII and GID. An increased activity of alpha-conotoxins GID and MII was also observed when the beta 2-F117A mutant was combined with the alpha 4 instead of the alpha 3 subunit. Investigation of A10L-PnIA indicated that high affinity binding to beta 2-F117A, beta 2-V109A, and beta 2-V109G mutants was conferred by amino acids with a long side chain in position 10 (PnIA numbering). Docking simulations of 4/7 alpha-conotoxin binding to the alpha 3 beta 2 model supported a direct interaction between mutated nAChR residues and alpha-conotoxin residues 6, 7, and 10. Taken together, these data provide evidence that the beta subunit contributes to alpha-conotoxin binding and selectivity and demonstrate that a small cleft leading to the agonist binding site is targeted by alpha-conotoxins to block the nAChR.

Toxin insights into nicotinic acetylcholine receptors

Venomous species have evolved cocktails of bioactive peptides to facilitate prey capture. Given their often exquisite potency and target selectivity, venom peptides provide unique biochemical tools for probing the function of membrane proteins at the molecular level. in the field of the nicotinic acetylcholine receptors (nAChRs), the subtype specific snake alpha-neurotoxins and cone snail alpha-conotoxins have been widely used to probe receptor structure and function in native tissues and recombinant systems. However, only recently has it been possible to generate an accurate molecular view of these nAChR-toxin interactions. Crystal structures of AChBP, a homologue of the nAChR ligand binding domain, have now been solved in complex with alpha-cobratoxin, alpha-conotoxin PnIA and alpha-conotoxin Iml. The orientation of all three toxins in the ACh binding site confirms many of the predictions obtained from mutagenesis and docking simulations on homology models of mammalian nAChR. The precise understanding of the molecular determinants of these complexes is expected to contribute to the development of more selective nAChR modulators. In this commentary, we review the structural data on nAChR-toxin interactions and discuss their implications for the design of novel ligands acting at the nAChR. (c) 2006 Elsevier Inc. All rights reserved.

Conformational dynamics associated with ligand binding to vertebrate hexa-coordinate hemoglobins

Neuroglobin (Ngb) and cytoglobin (Cygb) are two new additions to the globin family, exhibiting heme iron hexa-coordination, a disulfide bond and large internal cavities. These proteins are implicated in cytoprotection under hypoxic-ischemic conditions, but the molecular basis of their cytoprotective function is unclear. Herein, a photothermal and spectroscopic study of the interactions of diatomic ligands with Ngb, Cygb, myoglobin and hemoglobin is presented. The impact of the disulfide bond in Ngb and Cygb and role of conserved residues in Ngb His64, Val68, Cys55, Cys120 and Tyr44 on conformational dynamics associated with ligand binding/dissociation were investigated. Transient absorption and photoacoustic calorimetry studies indicate that CO photo-dissociation from Ngb leads to a volume expansion (13.4±0.9 mL mol-1), whereas a smaller volume change was determined for Ngb with reduced Cys (ΔV=4.6±0.3 mL mol-1). Furthermore, Val68 side chain regulates ligand migration between the distal pocket and internal hydrophobic cavities since Val68Phe geminate quantum yield is ∼2.7 times larger than that of WT Ngb. His64Gln and Tyr44Phe mutations alter the thermodynamic parameters associated with CO photo-release indicating that electrostatic/hydrogen binding network that includes heme propionate groups, Lys 67, His64, and Tyr 44 in Ngb modulates the energetics of CO photo-dissociation. In Cygb, CO escape from the protein matrix is fast (< 40 ns) with a ΔH of 18±2 kcal mol-1 in Cygbred, whereas disulfide bridge formation promotes a biphasic ligand escape associated with an overall enthalpy change of 9±4 kcal mol-1. Therefore, the disulfide bond modulates conformational dynamics in Ngb and Cygb. I propose that in Cygb with reduced Cys the photo-dissociated ligand escapes through the hydrophobic tunnel as occurs in Ngb, whereas the CO preferentially migrates through the His64 gate in Cygbox. To characterize Cygb surface 1,8-ANS interactions with Cygb were investigated employing fluorescence spectroscopy, ITC and docking simulations. Two 1,8-ANS binding sites were identified. One binding site is located close to the extended N-terminus of Cygb and was also identified as a binding site for oleate. Furthermore, guanidinium hydrochloride-induced unfolding studies of Cygb reveal that the disulfide bond does not impact Cygb stability, whereas binding of cyanide slightly increases the protein stability.

Développement et applications d’un outil bio-informatique pour la détection de similarités de champs d’interaction moléculaire

Résumé : Les méthodes de détection de similarités de sites de liaison servent entre autres à la prédiction de fonction et à la prédiction de cibles croisées. Ces méthodes peuvent aider à prévenir les effets secondaires, suggérer le repositionnement de médicament existants, identifier des cibles polypharmacologiques et des remplacements bio-isostériques. La plupart des méthodes utilisent des représentations basées sur les atomes, même si les champs d’interaction moléculaire (MIFs) représentent plus directement ce qui cherche à être identifié. Nous avons développé une méthode bio-informatique, IsoMif, qui détecte les similarités de MIF entre différents sites de liaisons et qui ne nécessite aucun alignement de séquence ou de structure. Sa performance a été comparée à d’autres méthodes avec des bancs d’essais, ce qui n’a jamais été fait pour une méthode basée sur les MIFs. IsoMif performe mieux en moyenne et est plus robuste. Nous avons noté des limites intrinsèques à la méthodologie et d’autres qui proviennent de la nature. L’impact de choix de conception sur la performance est discuté. Nous avons développé une interface en ligne qui permet la détection de similarités entre une protéine et différents ensembles de MIFs précalculés ou à des MIFs choisis par l’utilisateur. Des sessions PyMOL peuvent être téléchargées afin de visualiser les similarités identifiées pour différentes interactions intermoléculaires. Nous avons appliqué IsoMif pour identifier des cibles croisées potentielles de drogues lors d’une analyse à large échelle (5,6 millions de comparaisons). Des simulations d’arrimage moléculaire ont également été effectuées pour les prédictions significatives. L’objectif est de générer des hypothèses de repositionnement et de mécanismes d’effets secondaires observés. Plusieurs exemples sont présentés à cet égard.

Design and synthesis of novel quinolones directed to the Plasmodium falciparum bc1 protein complex

Tese de Doutoramento, Química, Especialização em Química Orgânica, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2016

In silico prediction of human carboxylesterase-1 (hCES1) metabolism combining docking analyses and MD simulations.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted in developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas little has been done to predict the hydrolytic activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES1. The study involves both docking analyses of known substrates to develop predictive models, and molecular dynamics (MD) simulations to reveal the in situ behavior of substrates and products, with particular attention being paid to the influence of their ionization state. The results emphasize some crucial properties of the hCES1 catalytic cavity, confirming that as a trend with several exceptions, hCES1 prefers substrates with relatively smaller and somewhat polar alkyl/aryl groups and larger hydrophobic acyl moieties. The docking results underline the usefulness of the hydrophobic interaction score proposed here, which allows a robust prediction of hCES1 catalysis, while the MD simulations show the different behavior of substrates and products in the enzyme cavity, suggesting in particular that basic substrates interact with the enzyme in their unprotonated form.

Histidine hydrogen bonding in MHC at pH 5 and pH 7 modeled by molecular docking and molecular dynamics simulations

Hydrogen bonds play important roles in maintaining the structure of proteins and in the formation of most biomolecular protein-ligand complexes. All amino acids can act as hydrogen bond donors and acceptors. Among amino acids, Histidine is unique, as it can exist in neutral or positively charged forms within the physiological pH range of 5.0 to 7.0. Histidine can thus interact with other aromatic residues as well as forming hydrogen bonds with polar and charged residues. The ability of His to exchange a proton lies at the heart of many important functional biomolecular interactions, including immunological ones. By using molecular docking and molecular dynamics simulation, we examine the influence of His protonation/deprotonation on peptide binding affinity to MHC class II proteins from locus HLA-DP. Peptide-MHC interaction underlies the adaptive cellular immune response, upon which the next generation of commercially-important vaccines will depend. Consistent with experiment, we find that peptides containing protonated His residues bind better to HLA-DP proteins than those with unprotonated His. Enhanced binding at pH 5.0 is due, in part, to additional hydrogen bonds formed between peptide His+ and DP proteins. In acidic endosomes, protein His79β is predominantly protonated. As a result, the peptide binding cleft narrows in the vicinity of His79β, which stabilizes the peptide - HLA-DP protein complex. © 2014 Bentham Science Publishers.

Data of the molecular dynamics simulations of mutations in the human connexin46 docking interface

The structure of hCx26 derived from the X-ray analysis was used to generate a homology model for hCx46. Interacting connexin molecules were used as starting model for the molecular dynamics (MD) simulation using NAMD and allowed us to predict the dynamic behavior of hCx46wt and the cataract related mutant hCx46N188T as well as two artificial mutants hCx46N188Q and hCx46N188D. Within the 50 ns simulation time the docked complex composed of the mutants dissociate while hCx46wt remains stable. The data indicates that one hCx46 molecule forms 5-7 hydrogen bonds (HBs) with the counterpart connexin of the opposing connexon. These HBs appear essential for a stable docking of the connexons as shown by the simulation of an entire gap junction channel and were lost for all the tested mutants. The data described here are related to the research article entitled "The cataract related mutation N188T in human connexin46 (hCx46) revealed a critical role for residue N188 in the docking process of gap junction channels" (Schadzek et al., 2015) [1].

Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

PACKMOL: A Package for Building Initial Configurations for Molecular Dynamics Simulations

Adequate initial configurations for molecular dynamics simulations consist of arrangements of molecules distributed in space in such a way to approximately represent the system`s overall structure. In order that the simulations are not disrupted by large van der Waals repulsive interactions, atoms from different molecules Must keep safe pairwise distances. Obtaining Such a molecular arrangement can be considered it packing problem: Each type molecule must satisfy spatial constraints related to the geometry of the system, and the distance between atoms of different molecules Must be greater than some specified tolerance. We have developed a code able to pack millions of atoms. grouped in arbitrarily complex molecules, inside a variety of three-dimensional regions. The regions may be intersections of spheres, ellipses, cylinders, planes, or boxes. The user must provide only the structure of one molecule of each type and the geometrical constraints that each type of molecule must satisfy. Building complex mixtures, interfaces, solvating biomolecules in water, other solvents, or mixtures of solvents, is straight forward. In addition. different atoms belonging to the same molecule may also be restricted to different spatial regions, in Such a way that more ordered molecular arrangements call be built, as micelles. lipid double-layers, etc. The packing time for state-of-the-art molecular dynamics systems varies front a few seconds to a few Minutes in a personal Computer. The input files are simple and Currently compatible with PDB, Tinker, Molden, or Moldy coordinate files. The package is distributed as free software and call be downloaded front http://www.ime.unicamp.br/similar to martinez/packmol/. (C) 2009 Wiley Periodicals. Inc. J Comput Chem 30: 2157-2164, 2009

Molecular Docking and Molecular Dynamic Studies of Semi-Synthetic Piperidine Alkaloids as Acetylcholinesterase Inhibitors

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)