Molecular Mechanisms of How Mercury Inhibits Water Permeation through Aquaporin-1: Understanding by Molecular Dynamics Simulation

Aquaporin (AQP) functions as a water-conducting pore. Mercury inhibits the water permeation through AQP. Although site-directed mutagenesis has shown that mercury binds to Cys¹⁸⁹ during the inhibition process, it is not fully understood how this inhibits the water permeation through AQP1. We carried out 40 ns molecular dynamics simulations of bovine AQP1 tetramer with mercury (Hg-AQP1) or without mercury (Free AQP1). In Hg-AQP1, Cys¹⁹¹ (Cys¹⁸⁹ in human AQP1) is converted to Cys-SHg⁺ in each monomer. During each last 10 ns, we observed water permeation events occurred 23 times in Free AQP1 and never in Hg-AQP1. Mercury binding did not influence the whole structure, but did induce a collapse in the orientation of several residues at the ar/R region. In Free AQP1, backbone oxygen atoms of Gly¹⁹⁰, Cys¹⁹¹, and Gly¹⁹² lined, and were oriented to, the surface of the water pore channel. In Hg-AQP1, however, the SHg⁺ of Cys191-SHg⁺ was oriented toward the outside of the water pore. As a result, the backbone oxygen atoms of Gly¹⁹⁰, Cys¹⁹¹, and Gly¹⁹² became disorganized and the ar/R region collapsed, thereby obstructing the permeation of water. We suggest that mercury disrupts the water pore of AQP1 through local conformational changes in the ar/R region.     

The structural pathway for water permeation through sodium-glucose cotransporters

Although water permeation across cell membranes occurs through several types of membrane proteins, the only permeation mechanism resolved at atomic scale is that through aquaporins. Crystallization of the Vibrio parahaemolyticus sodium-galactose transporter (vSGLT) allows investigation of putative water permeation pathways through both vSGLT and the homologous human Na-glucose cotransporter (hSGLT1) using computational methods. Grand canonical Monte Carlo and molecular dynamics simulations were used to stably insert water molecules in both proteins, showing the presence of a water-filled pathway composed of ∼100 water molecules. This provides a structural basis for passive water permeation that is difficult to reconcile with the water cotransport hypothesis. Potential-of-mean-force calculations of water going through the crystal structure of vSGLT shows a single barrier of 7.7 kCal mol⁻¹, in agreement with previously published experimental data for cotransporters of the SGLT family. Electrophysiological and volumetric experiments performed on hSGLT1-expressing Xenopus oocytes showed that the passive permeation pathway exists in different conformational states. In particular, experimental conditions that aim to mimic the conformation of the crystal structure displayed passive water permeability. These results provide groundwork for understanding the structural basis of cotransporter water permeability.     

Comparison of the structure of human recombinant short form stromelysin by multidimensional heteronuclear NMR and X-ray crystallography

Summary Stromelysin-1 is a matrix metalloprotease that has been implicated in a number of degenerative diseases. Here we present the refined NMR solution structure of the catalytic domain of stromelysin-1 complexed with a small inhibitor and compare it to the X-ray crystal structure of the same complex. The structures are similar in global fold and show an unusual bottomless S1' subsite. There are differences, however, in the least well defined regions, Phe⁸³-Ile⁸⁹, His²²⁴-Phe²³² and Pro²⁴⁹-Pro²⁵⁰, reflecting the lack of NOE data and large B-factors. The region His²²⁴-Phe²³² contains residues of the Sl' subsite and, consequently, small differences are observed in this subsite. Hydrogen-bond data show that, in contrast to the crystal structure, the solution structure lacks a hydrogen bond between the amide of Tyr²²³ and the carbonyl of the P3' residue. Analysis of bound water shows two tightly bound water molecules both in the solution and the crystal structure; neither of these waters are in the inhibitor binding site.Abbreviations MMP matrix metalloendoprotease - HMQC heteronuclear multiple quantum coherence - HSQC heteronuclear single quantum coherence - NOE nuclear Overhauser enhancement - NOESY NOE spectroscopy - PFG pulsed field gradient - sfSTR stromelysin-1 (EC 3.4.24.17), truncated at residue 255 The coordinates of the NMR solution structure (file name 2SRT) and the X-ray crystal structure (file name 1SLN) have been deposited in the Brookhaven Protein Data Bank.     

Activation and Proton Transport Mechanism in Influenza A M2 Channel

Molecular dynamics trajectories 2 μs in length have been generated for the pH-activated, tetrameric M2 proton channel of the influenza A virus in all protonation states of the pH sensor located at the His³⁷ tetrad. All simulated structures are in very good agreement with high-resolution structures. Changes in the channel caused by progressive protonation of His³⁷ provide insight into the mechanism of proton transport. The channel is closed at both His³⁷ and Trp⁴¹ sites in the singly and doubly protonated states, but it opens at Trp⁴¹ upon further protonation. Anions access the charged His³⁷ and by doing so stabilize the protonated states of the channel. The narrow opening at the His³⁷ site, further blocked by anions, is inconsistent with the water-wire mechanism of proton transport. Instead, conformational interconversions of His³⁷ correlated with hydrogen bonding to water molecules indicate that these residues shuttle protons in high-protonation states. Hydrogen bonds between charged and uncharged histidines are rare. The valve at Val²⁷ remains on average quite narrow in all protonation states but fluctuates sufficiently to support water and proton transport. A proton transport mechanism in which the channel, depending on pH, opens at either the histidine or valine gate is only partially supported by the simulations.     

Activation and Proton Transport Mechanism in Influenza A M2 Channel

Molecular dynamics trajectories 2 μs in length have been generated for the pH-activated, tetrameric M2 proton channel of the influenza A virus in all protonation states of the pH sensor located at the His³⁷ tetrad. All simulated structures are in very good agreement with high-resolution structures. Changes in the channel caused by progressive protonation of His³⁷ provide insight into the mechanism of proton transport. The channel is closed at both His³⁷ and Trp⁴¹ sites in the singly and doubly protonated states, but it opens at Trp⁴¹ upon further protonation. Anions access the charged His³⁷ and by doing so stabilize the protonated states of the channel. The narrow opening at the His³⁷ site, further blocked by anions, is inconsistent with the water-wire mechanism of proton transport. Instead, conformational interconversions of His³⁷ correlated with hydrogen bonding to water molecules indicate that these residues shuttle protons in high-protonation states. Hydrogen bonds between charged and uncharged histidines are rare. The valve at Val²⁷ remains on average quite narrow in all protonation states but fluctuates sufficiently to support water and proton transport. A proton transport mechanism in which the channel, depending on pH, opens at either the histidine or valine gate is only partially supported by the simulations.     

The Gating Mechanism of the Human Aquaporin 5 Revealed by Molecular Dynamics Simulations

Aquaporins are protein channels located across the cell membrane with the role of conducting water or other small sugar alcohol molecules (aquaglyceroporins). The high-resolution X-ray structure of the human aquaporin 5 (HsAQP5) shows that HsAQP5, as all the other known aquaporins, exhibits tetrameric structure. By means of molecular dynamics simulations we analyzed the role of spontaneous fluctuations on the structural behavior of the human AQP5. We found that different conformations within the tetramer lead to a distribution of monomeric channel structures, which can be characterized as open or closed. The switch between the two states of a channel is a tap-like mechanism at the cytoplasmic end which regulates the water passage through the pore. The channel is closed by a translation of the His67 residue inside the pore. Moreover, water permeation rate calculations revealed that the selectivity filter, located at the other end of the channel, regulates the flow rate of water molecules when the channel is open, by locally modifying the orientation of His173. Furthermore, the calculated permeation rates of a fully open channel are in good agreement with the reported experimental value.     

Mode of Binding of the Tuberculosis Prodrug Isoniazid to Heme Peroxidases: BINDING STUDIES AND CRYSTAL STRUCTURE OF BOVINE LACTOPEROXIDASE WITH ISONIAZID AT 2.7 ? RESOLUTION*

Isoniazid (INH) is an anti-tuberculosis prodrug that is activated by mammalian lactoperoxidase and Mycobacterium tuberculosis catalase peroxidase (MtCP). We report here binding studies, an enzyme assay involving INH, and the crystal structure of the complex of bovine lactoperoxidase (LPO) with INH to illuminate binding properties and INH activation as well as the mode of diffusion and interactions together with a detailed structural and functional comparison with MtCP. The structure determination shows that isoniazid binds to LPO at the substrate binding site on the distal heme side. The substrate binding site is connected to the protein surface through a long hydrophobic channel. The acyl hydrazide moiety of isoniazid interacts with Phe⁴²² O, Gln⁴²³ O^ϵ1, and Phe²⁵⁴ O. In this arrangement, pyridinyl nitrogen forms a hydrogen bond with a water molecule, W-1, which in turn forms three hydrogen bonds with Fe³⁺, His¹⁰⁹ N^ϵ2, and Gln¹⁰⁵ N^ϵ2. The remaining two sides of isoniazid form hydrophobic interactions with the atoms of heme pyrrole ring A, C^β and C^γ atoms of Glu²⁵⁸, and C^γ and C^δ atoms of Arg²⁵⁵. The binding studies indicate that INH binds to LPO with a value of 0.9 × 10⁻⁶ m for the dissociation constant. The nitro blue tetrazolium reduction assay shows that INH is activated by the reaction of LPO-H₂O₂ with INH. This suggests that LPO can be used for INH activation. It also indicates that the conversion of INH into isonicotinoyl radical by LPO may be the cause of INH toxicity.     

Interaction of divalent metal ions with the carboxyl-terminal domain of human voltage-gated proton channel Hv1

The voltage-gated proton channel Hv1 functions as a dimer, in which the intracellular C-terminal domain of the protein is responsible for the dimeric architecture and regulates proton permeability. Although it is well known that divalent metal ions have effect on the proton channel activity, the interaction of divalent metal ions with the channel in detail is not well elucidated. Herein, we investigated the interaction of divalent metal ions with the C-terminal domain of human Hv1 by CD spectra and fluorescence spectroscopy. The divalent metal ions binding induced an obvious conformational change at pH 7 and a pH-sensitive reduction of thermostability in the C-terminal domain. The interactions were further estimated by fluorescence spectroscopy experiments. There are at least two binding sites for divalent metal ions binding to the C-terminal domain of Hv1, either of which is close to His²⁴⁴ or His²⁶⁶ residue. The binding of Zn²⁺ to the two sites both enhanced the fluorescence of the protein at pH 7, whereas the binding of other divalent metal ions to the two sites all resulted fluorescence quenching. The orders of the strength of divalent metal ions binding to the two sites from strong to weak are both Co²⁺, Ca²⁺, Ni²⁺, Mg²⁺, and Mn²⁺. The strength of Ca²⁺, Co²⁺, Mg²⁺, Mn²⁺ and Ni²⁺ binding to the site close to His²⁴⁴ is stronger than that of these divalent metal ions binding to the site close to His²⁶⁶.     

Cu(II)-Binding N-Terminal Sequences of Human Proteins

Proteins anchor copper(II) ions mainly by imidazole from histidine residues located in different positions in the primary protein structures. However, the motifs with histidine in the first three N-terminal positions (His¹, His², and His³) show unique Cu(II)-binding properties, such as availability from the surface of the protein, high flexibility, and high Cu(II) exchangeability with other ligands. It makes such sequences beneficial for the fast exchange of Cu(II) between ligands. Furthermore, sequences with His¹ and His², thus, non-saturating the Cu(II) coordination sphere, are redox-active and may play a role in Cu(II) reduction to Cu(I). All human protein sequences deposited in UniProt Knowledgebase were browsed for those containing His¹, His², or His³. Proteolytically modified sequences (with the removal of a propeptide or Met residue) were taken for the analysis. Finally, the sequences were sorted out according to the subcellular localization of the proteins to match the respective sequences with the probability of interaction with divalent copper.     

Life: Self-Directing with Unlimited Variability on Self-Speeding

Abstract

The crystal structure of the deoxyoctamer d(G-G-^Br U-A-^BrU-A-C-C) was refined to a resolution of 1.7Å using combined diffractometer and synchrotron data. The analysis was carried out independently in two laboratories using different procedures. Although the final results are identical the comparison of the two approaches highlights potential problems in the refinement of oligonucleotides when only limited data are available.

As part of the analysis the positions of 84 solvent molecules in the asymmetric unit were established. The DNA molecule is highly solvated, particularly the phosphate-sugar backbone and the functional groups of the bases. The major groove contains, in the central ^BrU-A-^BrU-A region, a ribbon of water molecules forming closed pentagons with shared edges. These water molecules are linked to the base O and N atoms and to the solvent chains connecting the O-1 phosphate oxygen atoms on each strand. The minor groove is also extensively hydrated with a continuous network in the central region and other networks at each end. The pattern of hydration is briefly compared with that observed in the crystal structure of a B-dodecamer.     

Energy distribution and ion selectivity of the bacterial potassium channel

The ion selectivity of the bacterial potassium channel K_CSA is explained upon comparing the energy characteristics of the interaction of cations (Li⁺, Na⁺, K⁺) with atoms of the selectivity filter of the protein pore. Quantum-chemical calculations reveal a deeper potential well for potassium ions, which accounts for preferred K⁺ permeation. It is shown that the conventional methods with AMBER, CHARMM, OPLS force fields in standard parametrization as well as partial re-parametrization give incorrect estimates of ion energy distribution in the channel.     

Characterization of a dipartite iron uptake system from uropathogenic Escherichia coli strain F11

In the uropathogenic Escherichia coli strain F11, in silico genome analysis revealed the dicistronic iron uptake operon fetMP, which is under iron-regulated control mediated by the Fur regulator. The expression of fetMP in a mutant strain lacking known iron uptake systems improved growth under iron depletion and increased cellular iron accumulation. FetM is a member of the iron/lead transporter superfamily and is essential for iron uptake by the Fet system. FetP is a periplasmic protein that enhanced iron uptake by FetM. Recombinant FetP bound Cu(II) and the iron analog Mn(II) at distinct sites. The crystal structure of the FetP dimer reveals a copper site in each FetP subunit that adopts two conformations: CuA with a tetrahedral geometry composed of His⁴⁴, Met⁹⁰, His⁹⁷, and His¹²⁷, and CuB, a second degenerate octahedral geometry with the addition of Glu⁴⁶. The copper ions of each site occupy distinct positions and are separated by ∼1.3 Å. Nearby, a putative additional Cu(I) binding site is proposed as an electron source that may function with CuA/CuB displacement to reduce Fe(III) for transport by FetM. Together, these data indicate that FetMP is an additional iron uptake system composed of a putative iron permease and an iron-scavenging and potentially iron-reducing periplasmic protein.     

The hydration shell of myoglobin

The space in the unit cell of a metmyoglobin crystal not occupied by myoglobin atoms was filled with water using Monte Carlo calculations. Independent calculations with different amounts of water have been performed. Structure factors were calculated using the water coordinates thus obtained and the known coordinates of the myoglobin atoms. A comparison with experimental structure factors showed that both the low and the high resolution regime could be well reproduced with 814 Monte Carlo water molecules per unit cell with a B-value of 50 Ų. The Monte Carlo water molecules yield a smaller standard R-value (0.166) than using a homogeneous electron density for the simulation of the crystal water (R = 0.212). A reciprocal space refinement of the water and the protein coordinates has been performed. Monte Carlo calculations can be used to obtain information for crystallographically invisible parts of the unit cell and yield better coordinates for the visible part in the refinement. Correspondence to: F. Parak     

Probing internal water molecules in proteins using two-dimensional 19F-1H NMR

Summary A simple approach for detecting internal water molecules in proteins in solution is described. This approach combines ¹⁹F-detected heteronuclear Overhauser and exchange spectroscopy (HOESY) with site-specific ¹⁹F substitution. The model system employed was intestinal fatty acid-binding protein complexed with [2-mono-¹⁹F]-palmitate. An intense cross peak was observed between the fluorine and a buried water molecule, as defined in the 1.98 Å crystal structure of the complex. From HOESY spectra, the fluorine-water distance was estimated to be 2.1 Å, in agreement with the crystal structure. This approach may be applicable to macromolecules that are too large for ¹H-detected NMR methods.     

Crystal Structure of the Pig Pancreatic α-Amylase Complexed with Malto-Oligosaccharides

The structural X-ray map of a pig pancreatic α-amylase crystal soaked (and flash-frozen) with a maltopentaose substrate showed a pattern of electron density corresponding to the binding of oligosaccharides at the active site and at three surface binding sites. The electron density region observed at the active site, filling subsites ?3 through ?1, was interpreted in terms of the process of enzyme-catalyzed hydrolysis undergone by maltopentaose. Because the expected conformational changes in the “flexible loop” that constitutes the surface edge of the active site were not observed, the movement of the loop may depend on aglycone site being filled. The crystal structure was refined at 2.01 å resolution to an R factor of 17.0% (R _free factor of 19.8%). The final model consists of 3910 protein atoms, one calcium ion, two chloride ions, 103 oligosaccharide atoms, 761 atoms of water molecules, and 23 ethylene glycol atoms.     

Display and interpretation of solvent electron density distributions in insulin crystals

In macromolecular crystallography, three-dimensional contour surfaces are useful for interactive computer graphics displays of the protein electron density but are less effective for presenting static images of large volumes of solvent density. A raster-based computer graphics program which displays depth-cued projections of continuous density distributions has been developed to analyze the distribution of solvent atoms in macromolecular crystals. Maps of the water distribution in the cubic insulin crystal show some well-ordered waters, which are bound to surrounding protein atoms by multiple hydrogen bonds, and an ill-defined solvent structure at a greater distance from the protein surface. Molecular dynamics calculations were used to assist in the interpretation of the time-varying solvent structure within two enclosed cavities in the crystal. Two water molecules that ligate a sodium ion were almost immobile during the simulation but the majority of water molecules were found to move rapidly between the density maxima identified from the crystallographic refinement.     

Comparative study of the energetics of ion permeation in Kv1.2 and KcsA potassium channels

Biological ion channels rely on a multi-ion transport mechanism for fast yet selective permeation of ions. The crystal structure of the KcsA potassium channel provided the first microscopic picture of this process. A similar mechanism is assumed to operate in all potassium channels, but the validity of this assumption has not been well investigated. Here, we examine the energetics of ion permeation in Shaker Kv1.2 and KcsA channels, which exemplify the six-transmembrane voltage-gated and two-transmembrane inward-rectifier channels. We study the feasibility of binding a third ion to the filter and the concerted motion of ions in the channel by constructing the potential of mean force for K⁺ ions in various configurations. For both channels, we find that a pair of K⁺ ions can move almost freely within the filter, but a relatively large free-energy barrier hinders the K⁺ ion from stepping outside the filter. We discuss the effect of the CMAP dihedral energy correction that was recently incorporated into the CHARMM force field on ion permeation dynamics.     

Crystal structure of HugZ, a novel heme oxygenase from Helicobacter pylori

The crystal structure of a heme oxygenase (HO) HugZ from Helicobacter pylori complexed with heme has been solved and refined at 1.8 Å resolution. HugZ is part of the iron acquisition mechanism of H. pylori, a major pathogen of human gastroenteric diseases. It is required for the adaptive colonization of H. pylori in hosts. Here, we report that HugZ is distinct from all other characterized HOs. It exists as a dimer in solution and in crystals, and the dimer adopts a split-barrel fold that is often found in FMN-binding proteins but has not been observed in hemoproteins. The heme is located at the intermonomer interface and is bound by both monomers. The heme iron is coordinated by the side chain of His²⁴⁵ and an azide molecule when it is present in crystallization conditions. Experiments show that Arg¹⁶⁶, which is involved in azide binding, is essential for HugZ enzymatic activity, whereas His²⁴⁵, surprisingly, is not, implying that HugZ has an enzymatic mechanism distinct from other HOs. The placement of the azide corroborates the observed γ-meso specificity for the heme degradation reaction, in contrast to most known HOs that have α-meso specificity. We demonstrate through sequence and structural comparisons that HugZ belongs to a new heme-binding protein family with a split-barrel fold. Members of this family are widespread in pathogenic bacteria and may play important roles in the iron acquisition of these bacteria.     

A theory for the proton transport of the influenza virus M2 protein: extensive test against conductance data

A theory for calculating the proton flux through the influenza virus M2 channel is tested here against an extensive set of conductance data. The flux is determined by the rate constants for binding to the His³⁷ tetrad from the two sides of the membrane and the corresponding unbinding rate constants. The rate constants are calculated by explicitly treating the structure and dynamics of the protein. Important features revealed by previous studies, such as a gating role for Val²⁷ at the entrance to the channel pore, and channel activation by viral exterior pH, are incorporated in this theory. This study demonstrates that the conductance function of the M2 proton channel can now be quantitatively rationalized by the structure and dynamics of the protein.