Ratiometric fluorescent scaffold giving discrete response towards iodide ion: a combined experimental and DFT study

A novel ratiometric fluorescent receptor (Z)‐2‐(4‐[diethylamino]‐2‐hydroxybenzylideneamino) pyridine‐3‐carbaldehyde (3) bearing one phenolic OH and one aldehyde group as recognition sites was synthesized and characterized. The anion recognition behaviour of receptor 3 was evaluated by various spectroscopic (UV‐visible, fluorescence and ¹H nuclear magnetic resonance) methods and was validated by computational studies. The receptor showed fast response time, excellent selectivity and reproducibility towards iodide ion detection among the other surveyed anions, with a binding constant of 6.12 × 10⁴ M^?1 and a detection limit of 0.24 μM, thus confirming its potential applicability as a fluorescent sensor for iodide Copyright © 2014 John Wiley & Sons, Ltd.     

The Role of Hydrated Divalent Metal Ions in the Bridging of Two Anionic Groups. An ab initio Quantum Chemical and Molecular Mechanics Study of Dimethyl Phosphate and Formate Bridged by Calcium and Magnesium Ions

Abstract

Ab initio quantum chemical (Gaussian82) and molecular mechanics (AMBER2.0) computational techniques are employed to investigate the interaction of twoanions (formate and dimethyl phosphate) and a central divalent metal cation (magnesium or calcium). These systems are models for the essential GDP binding unit of the G-proteins (e.g., EF-Tu or the ras oncogene proteins) and for protein/phospholipid interactions, both of which are mediated by divalent metal cations. Various levels of hydration are utilized to examine coordination of differences between magnesium and calcium ions. Two different orientations of formate and dimethyl phosphate in direct ion contact with a magnesium ion and two waters of hydration were energy minimized with both quantum and molecular mechanics techniques. The structures and energy differences between the two orientations determined by either of the computational techniques are similar. Magnesium ion has a strong propensity to assume six coordination whereas calcium ion preferentially assumes a coordination greater than six. Likewise, water molecules attached to magnesium ion are held more rigidly than those to calcium ion, thus calcium ion is more accommodating in the exchange of water for negative ligands.     

Experimental determination of the vertical alignment between the second and third transmembrane segments of muscle nicotinic acetylcholine receptors

Nicotinic acetylcholine receptors (nAChR) are members of the Cys‐loop ligand‐gated ion channel superfamily. Muscle nAChR are heteropentamers that assemble from two α, and one each of β, γ, and δ subunits. Each subunit is composed of three domains, extracellular, transmembrane and intracellular. The transmembrane domain consists of four α‐helical segments (M1–M4). Pioneering structural information was obtained using electronmicroscopy of Torpedo nAChR. The recently solved X‐ray structure of the first eukaryotic Cys‐loop receptor, a truncated (intracellular domain missing) glutamate‐gated chloride channel α (GluClα) showed the same overall architecture. However, a significant difference with regard to the vertical alignment between the channel‐lining segment M2 and segment M3 was observed. Here, we used functional studies utilizing disulfide trapping experiments in muscle nAChR to determine the spatial orientation between M2 and M3. Our results are in agreement with the vertical alignment as obtained when using the GluClα structure as a template to homology model muscle nAChR, however, they cannot be reconciled with the current Torpedo nAChR model. The vertical M2–M3 alignments as observed in X‐ray structures of prokaryotic Gloeobacter violaceus ligand‐gated ion channel and GluClα are in agreement. Our results further confirm that this alignment in Cys‐loop receptors is conserved between prokaryotes and eukaryotes.     

Modeling the human PTC bitter-taste receptor interactions with bitter tastants

We employed the first principles computational method MembStruk and homology modeling techniques to predict the 3D structures of the human phenylthiocarbamide (PTC) taste receptor. This protein is a seven-transmembrane-domain G protein-coupled receptor that exists in two main forms worldwide, designated taster and nontaster, which differ from each other at three amino-acid positions. 3D models were generated with and without structural similarity comparison to bovine rhodopsin. We used computational tools (HierDock and ScanBindSite) to generate models of the receptor bound to PTC ligand to estimate binding sites and binding energies. In these models, PTC binds at a site distant from the variant amino acids, and PTC binding energy was equivalent for both the taster and nontaster forms of the protein. These models suggest that the inability of humans to taste PTC is due to a failure of G protein activation rather than decreased binding affinity of the receptor for PTC. Amino-acid substitutions in the sixth and seventh transmembrane domains of the nontaster form of the protein may produce increased steric hindrance between these two α-helices and reduce the motion of the sixth helix required for G protein activation.     

A BEST example of channel structure annotation by molecular simulation

An increasing number of ion channel structures are being determined. This generates a need for computational tools to enable functional annotation of channel structures. However, several studies of ion channel and model pores have indicated that the physical dimensions of a pore are not always a reliable indicator of its conductive status. This is due to the unusual behavior of water within nano-confined spaces, resulting in a phenomenon referred to as “hydrophobic gating”. We have recently demonstrated how simulating the behavior of water within an ion channel pore can be used to predict its conductive status. In this addendum to our study, we apply this method to compare the recently solved structure of a mutant of the bestrophin chloride channel BEST1 with that of the wild-type channel. Our results support the hypothesis of a hydrophobic gate within the narrow neck of BEST1. This provides further validation that this simulation approach provides the basis for an accurate and computationally efficient tool for the functional annotation of ion channel structures.     

Crystallographic and Spectroscopic Studies of d(CCGGTACCGG)

The decanucleotide sequence d(CCGGTACCGG) crystallizes as a four-way junction at low cobalt ion concentrations (i.e., 1 mM). When the cobalt concentration in the crystallization solution is increased to 5 mM, the sequence crystallizes as resolved B-DNA duplexes. Gel retardation studies of the decamer show both a faint slow-moving band and a much thicker fast-moving band at low cobalt ion concentrations, and only the intense fast-moving band at higher ion concentration. Circular dichroism (CD) spectroscopy of the decamer indicates a structural transition as the cobalt ion concentration in the solution is increased, probably from B-type to A-type DNA. These studies revealed that the oligomer sequence has several conformations and structures accessible to it, in a manner dependent on sequence, ion concentration, and DNA concentration.

[Supplementary materials are available for this article. Go to the publisher's online edition of Nucleosides, Nucleotides & Nucleic Acids for the following free supplemental resources(s): Supplementary Figures 1, 2, and 3.]     

A Nanosecond Molecular Dynamics Study of Antiparallel d(G)7 Quadruplex Structures: Effect of the Coordinated Cations

Abstract

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G₇) quadruplex structures with different coordinated ions, namely Na⁺ and K⁺ ion, water and Na⁺ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of ～ 1.5 Å from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K⁺ ions are quite immobile. MD studies at 400K indicate that K⁺ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K⁺ ion coordinated structure is more favourable than that for the Na⁺ ion coordinated antiparallel quadruplex structure.     

Structure of the zinc‐bound amino‐terminal domain of the NMDA receptor NR2B subunit

N‐methyl‐D ‐aspartate (NMDA) receptors belong to the family of ionotropic glutamate receptors (iGluRs) that mediate the majority of fast excitatory synaptic transmission in the mammalian brain. One of the hallmarks for the function of NMDA receptors is that their ion channel activity is allosterically regulated by binding of modulator compounds to the extracellular amino‐terminal domain (ATD) distinct from the L ‐glutamate‐binding domain. The molecular basis for the ATD‐mediated allosteric regulation has been enigmatic because of a complete lack of structural information on NMDA receptor ATDs. Here, we report the crystal structures of ATD from the NR2B NMDA receptor subunit in the zinc‐free and zinc‐bound states. The structures reveal the overall clamshell‐like architecture distinct from the non‐NMDA receptor ATDs and molecular determinants for the zinc‐binding site, ion‐binding sites, and the architecture of the putative phenylethanolamine‐binding site.     

Crystal Structural Framework of Lithium Super‐Ionic Conductors

As technologically important materials for solid‐state batteries, Li super‐ionic conductors are a class of materials exhibiting exceptionally high ionic conductivity at room temperature. These materials have unique crystal structural frameworks hosting a highly conductive Li sublattice. However, it is not understood why certain crystal structures of the super‐ionic conductors lead to high conductivity in the Li sublattice. In this study, using topological analysis and ab initio molecular dynamics simulations, the crystal structures of all Li‐conducting oxides and sulfides are studied systematically and the key features pertaining to fast‐ion conduction are quantified. In particular, a unique feature of enlarged Li sites caused by large local spaces in the crystal structural framework is identified, promoting fast conduction in the Li‐ion sublattice. Based on these quantified features, the high‐throughput screening identifies many new structures as fast Li‐ion conductors, which are further confirmed by ab initio molecular dynamics simulations. This study provides new insights and a systematic quantitative understanding of the crystal structural frameworks of fast ion‐conductor materials and motivates future experimental and computational studies on new fast‐ion conductors.     

Prediction of structures of zinc‐binding proteins through explicit modeling of metal coordination geometry

Metal ions play an essential role in stabilizing protein structures and contributing to protein function. Ions such as zinc have well‐defined coordination geometries, but it has not been easy to take advantage of this knowledge in protein structure prediction efforts. Here, we present a computational method to predict structures of zinc‐binding proteins given knowledge of the positions of zinc‐coordinating residues in the amino acid sequence. The method takes advantage of the “atom‐tree” representation of molecular systems and modular architecture of the Rosetta3 software suite to incorporate explicit metal ion coordination geometry into previously developed de novo prediction and loop modeling protocols. Zinc cofactors are tethered to their interacting residues based on coordination geometries observed in natural zinc‐binding proteins. The incorporation of explicit zinc atoms and their coordination geometry in both de novo structure prediction and loop modeling significantly improves sampling near the native conformation. The method can be readily extended to predict protein structures bound to other metal and/or small chemical cofactors with well‐defined coordination or ligation geometry.     

K+ and Na+ Conduction in Selective and Nonselective Ion Channels Via Molecular Dynamics Simulations

Generations of scientists have been captivated by ion channels and how they control the workings of the cell by admitting ions from one side of the cell membrane to the other. Elucidating the molecular determinants of ion conduction and selectivity are two of the most fundamental issues in the field of biophysics. Combined with ongoing progress in structural studies, modeling and simulation have been an integral part of the development of the field. As of this writing, the relentless growth in computational power, the development of new algorithms to tackle the so-called rare events, improved force-field parameters, and the concomitant increasing availability of membrane protein structures, allow simulations to contribute even further, providing more-complete models of ion conduction and selectivity in ion channels. In this report, we give an overview of the recent progress made by simulation studies on the understanding of ion permeation in selective and nonselective ion channels.     

K and Na Conduction in Selective and Nonselective Ion Channels Via Molecular Dynamics Simulations

Generations of scientists have been captivated by ion channels and how they control the workings of the cell by admitting ions from one side of the cell membrane to the other. Elucidating the molecular determinants of ion conduction and selectivity are two of the most fundamental issues in the field of biophysics. Combined with ongoing progress in structural studies, modeling and simulation have been an integral part of the development of the field. As of this writing, the relentless growth in computational power, the development of new algorithms to tackle the so-called rare events, improved force-field parameters, and the concomitant increasing availability of membrane protein structures, allow simulations to contribute even further, providing more-complete models of ion conduction and selectivity in ion channels. In this report, we give an overview of the recent progress made by simulation studies on the understanding of ion permeation in selective and nonselective ion channels.     

Analysis and modeling of heat-labile enterotoxins of Escherichia coli suggests a novel space with insights into receptor preference

Features of heat-labile enterotoxins of Escherichia coli which make them fit to use as novel receptors for antidiarrheals are not completely explored. Data-set of 14 different serovars of enterotoxigenic Escherichia coli producing heat-labile toxins were taken from NCBI Genbank database and used in the study. Sequence analysis showed mutations in different subunits and also at their interface residues. As these toxins lack crystallography structures, homology modeling using Modeller 9.11 led to the structural approximation for the E. coli producing heat-labile toxins. Interaction of modeled toxin subunits with proanthocyanidin, an antidiarrheal showed several strong hydrogen bonding interactions at the cost of minimized energy. The hits were subsequently characterized by molecular dynamics simulation studies to monitor their binding stabilities. This study looks into novel space where the ligand can choose the receptor preference not as a whole but as an individual subunit. Mutation at interface residues and interaction among subunits along with the binding of ligand to individual subunits would help to design a non-toxic labile toxin and also to improve the therapeutics.     

MED: a new non-supervised gene prediction algorithm for bacterial and archaeal genomes

Background  

Despite a remarkable success in the computational prediction of genes in Bacteria and Archaea, a lack of comprehensive understanding of prokaryotic gene structures prevents from further elucidation of differences among genomes. It continues to be interesting to develop new ab initio algorithms which not only accurately predict genes, but also facilitate comparative studies of prokaryotic genomes.     

SANA: an algorithm for sequential and non-sequential protein structure alignment

Protein structure alignment algorithms play an important role in the studies of protein structure and function. In this paper, a novel approach for structure alignment is presented. Specifically, core regions in two protein structures are first aligned by identifying connected components in a network of neighboring geometrically compatible aligned fragment pairs. The initial alignments then are refined through a multi-objective optimization method. The algorithm can produce both sequential and non-sequential alignments. We show the superior performance of the proposed algorithm by the computational experiments on several benchmark datasets and the comparisons with the well-known structure alignment algorithms such as DALI, CE and MATT. The proposed method can obtain accurate and biologically significant alignment results for the case with occurrence of internal repeats or indels, identify the circular permutations, and reveal conserved functional sites. A ranking criterion of our algorithm for fold similarity is presented and found to be comparable or superior to the Z-score of CE in most cases from the numerical experiments. The software and supplementary data of computational results are available at .     

Correlated Migration Invokes Higher Na+‐Ion Conductivity in NaSICON‐Type Solid Electrolytes

Na super ion conductor (NaSICON), Na_1+nZr₂Si_nP_3–nO₁₂ is considered one of the most promising solid electrolytes; however, the underlying mechanism governing ion transport is still not fully understood. Here, the existence of a previously unreported Na5 site in monoclinic Na₃Zr₂Si₂PO₁₂ is unveiled. It is revealed that Na⁺‐ions tend to migrate in a correlated mechanism, as suggested by a much lower energy barrier compared to the single‐ion migration barrier. Furthermore, computational work uncovers the origin of the improved conductivity in the NaSICON structure, that is, the enhanced correlated migration induced by increasing the Na⁺‐ion concentration. Systematic impedance studies on doped NaSICON materials bolster this finding. Significant improvements in both the bulk and total ion conductivity (e.g., σ_bulk = 4.0 mS cm^?1, σ_total = 2.4 mS cm^?1 at 25 °C) are achieved by increasing the Na content from 3.0 to 3.30–3.55 mol formula unit^?1. These improvements stem from the enhanced correlated migration invoked by the increased Coulombic repulsions when more Na⁺‐ions populate the structure rather than solely from the increased mobile ion carrier concentration. The studies also verify a strategy to enhance ion conductivity, namely, pushing the cations into high energy sites to therefore lower the energy barrier for cation migration.     

The Influence of Monovalent Cation Size on the Stability of RNA Tertiary Structures

Many RNA tertiary structures are stable in the presence of monovalent ions alone. To evaluate the degree to which ions at or near the surfaces of such RNAs contribute to stability, the salt-dependent stability of a variety of RNA structures was measured with each of the five group I cations. The stability of hairpin secondary structures and a pseudoknot tertiary structure are insensitive to the ion identity, but the tertiary structures of two other RNAs, an adenine riboswitch and a kissing loop complex, become more stable by 2-3 kcal/mol as ion size decreases. This "default" trend is attributed to the ability of smaller ions to approach the RNA surface more closely. The degree of cation accumulation around the kissing loop complex was also inversely proportional to ion radius, perhaps because of the presence of sterically restricted pockets that can be accessed only by smaller ions. An RNA containing the tetraloop-receptor motif shows a strong (up to ∼ 3 kcal/mol) preference for Na⁺ or K⁺ over other group I ions, consistent with the chelation of K⁺ by this motif in some crystal structures. This RNA reverts to the default dependence on ion size when a base forming part of the chelation site is mutated. Lastly, an RNA aptamer for cobinamide, which was originally selected in the presence of high concentrations of LiCl, binds ligand more strongly in the presence of Li⁺ than other monovalent ions.On the basis of these trends in RNA stability with group I ion size, it is argued that two features of RNA tertiary structures may promote strong interactions with ions at or near the RNA surface: negative charge densities that are higher than that in secondary structures, and the occasional presence of chelation sites, which are electronegative pockets that selectively bind ions of an optimum size.     

Structural dynamics of an ionotropic glutamate receptor

Ionotropic glutamate receptors (iGluRs) are postsynaptic ion channels involved in excitatory neurotransmission. iGluRs play important roles in development and in forms of synaptic plasticity that underlie higher order processes such as learning and memory. Neurobiological and biochemical studies have long characterized iGluRs in detail. However, the structural basis for the function of iGluRs has not yet been investigated, because there is insufficient information about their three-dimensional structures. In 1998, a crystal structure called S1S2 lobes was first solved for the extracellular bilobed ligand-binding domain of the GluR2 subunit. Since then, the crystal structures for the S1S2 lobes both in the apo and in various liganded states have been reported, and recent biophysical studies have further elucidated the dynamic aspects of the structure of the S1S2 lobes. In this review, the dynamic structures of the S1S2 lobes and their ligands are summarized, and the importance of their structural flexibility and fluctuation is discussed in light of the mechanisms of ligand recognition, activation, and desensitization of the receptor.     

Synthesis and Computational Exploration of Morpholine Bearing Halogenated Sulfonamides as Potential Tyrosinase Inhibitors

In the presented work, a new series of three different 4-((3,5-dichloro-2-[(2/4-halobenzyl)oxy]phenyl)sulfonyl)morpholines was synthesized and the structure of these compounds were corroborated by ¹H-NMR & ¹³C-NMR studies. The in vitro results established all the three compounds as potent tyrosinase inhibitors relative to the standard. The Kinetics mechanism plots established that compound 8 inhibited the enzyme non-competitively. The inhibition constants K_i calculated from Dixon plots for this compound was 0.0025 μM. Additionally, computational techniques were used to explore electronic structures of synthesized compounds. Fully optimized geometries were further docked with tyrosinase enzyme for inhibition studies. Reasonably good binding/interaction energies and intermolecular interactions were obtained. Finally, drug likeness was also predicted using the rule of five (RO5) and Chemical absorption, distribution, metabolism, excretion, and toxicity (ADMET) characteristics. It is anticipated that current experimental and computational investigations will evoke the scientific interest of the research community for the above-entitled compounds.     

Analysis of secondary structure and self‐assembly of amelogenin by variable temperature circular dichroism and isothermal titration calorimetry

Amelogenin is a proline‐rich enamel matrix protein known to play an important role in the oriented growth of enamel crystals. Amelogenin self‐assembles to form nanospheres and higher order structures mediated by hydrophobic interactions. This study aims to obtain a better insight into the relationship between primary–secondary structure and self‐assembly of amelogenin by applying computational and biophysical methods. Variable temperature circular dichroism studies indicated that under physiological pH recombinant full‐length porcine amelogenin contains unordered structures in equilibrium with polyproline type II (PPII) structure, the latter being more populated at lower temperatures. Increasing the concentration of rP172 resulted in the promotion of folding to an ordered β‐structured assembly. Isothermal titration calorimetry dilution studies revealed that at all temperatures, self‐assembly is entropically driven due to the hydrophobic effect and the molar heat of assembly (ΔH_A) decreases with temperature. Using a computational approach, a profile of domains in the amino acid sequence that have a high propensity to assemble and to have PPII structures has been identified. We conclude that the assembly properties of amelogenin are due to complementarity between the hydrophobic and PPII helix prone regions. Proteins 2009. © 2009 Wiley‐Liss, Inc.