What Makes an Aquaporin a Glycerol Channel? A Comparative Study of AqpZ and GlpF

Structure of urea transporters

Members of the urea transporter (UT) family mediate rapid, selective transport of urea down its concentration gradient. To date, crystal structures of two evolutionarily distant UTs have been solved. These structures reveal a common UT fold involving two structurally homologous domains that encircle a continuous membrane-spanning pore and indicate that UTs transport urea via a channel-like mechanism. Examination of the conserved architecture of the pore, combined with crystal structures of ligand-bound proteins, molecular dynamics simulations, and functional data on permeation and inhibition by a broad range of urea analogs and other small molecules, provides insight into the structural basis of urea permeation and selectivity.

Subangstrom resolution X-ray structure details aquaporin-water interactions

Aquaporins are membrane channels that facilitate the flow of water across biological membranes. Two conserved regions are central for selective function: the dual asparagine-proline-alanine (NPA) aquaporin signature motif and the aromatic and arginine selectivity filter (SF). Here, we present the crystal structure of a yeast aquaporin at 0.88 angstrom resolution. We visualize the H-bond donor interactions of the NPA motif's asparagine residues to passing water molecules; observe a polarized water-water H-bond configuration within the channel; assign the tautomeric states of the SF histidine and arginine residues; and observe four SF water positions too closely spaced to be simultaneously occupied. Strongly correlated movements break the connectivity of SF waters to other water molecules within the channel and prevent proton transport via a Grotthuss mechanism.

The Aquaporin Channel Repertoire of the Tardigrade Milnesium tardigradum

Limno-terrestrial tardigrades are small invertebrates that are subjected to periodic drought of their micro-environment. They have evolved to cope with these unfavorable conditions by anhydrobiosis, an ametabolic state of low cellular water. During drying and rehydration, tardigrades go through drastic changes in cellular water content. By our transcriptome sequencing effort of the limno-terrestrial tardigrade Milnesium tardigradum and by a combination of cloning and targeted sequence assembly, we identified transcripts encoding eleven putative aquaporins. Analysis of these sequences proposed 2 classical aquaporins, 8 aquaglyceroporins and a single potentially intracellular unorthodox aquaporin. Using quantitative real-time PCR we analyzed aquaporin transcript expression in the anhydrobiotic context. We have identified additional unorthodox aquaporins in various insect genomes and have identified a novel common conserved structural feature in these proteins. Analysis of the genomic organization of insect aquaporin genes revealed several conserved gene clusters.

Barriers to superfast water transport in carbon nanotube membranes

Carbon nanotube (CNT) membranes hold the promise of extraordinary fast water transport for applications such as energy efficient filtration and molecular level drug delivery. However, experiments and computations have reported flow rate enhancements over continuum hydrodynamics that contradict each other by orders of magnitude. We perform large scale molecular dynamics simulations emulating for the first time the micrometer thick CNTs membranes used in experiments. We find transport enhancement rates that are length dependent due to entrance and exit losses but asymptote to 2 orders of magnitude over the continuum predictions. These rates are far below those reported experimentally. The results suggest that the reported superfast water transport rates cannot be attributed to interactions of water with pristine CNTs alone.

Phosphorylation of Ser-180 of rat aquaporin-4 shows marginal affect on regulation of water permeability: molecular dynamics study

Water permeation through rat aquaporin-4 (rAQP4), predominantly found in mammalian brain is regulated by phosphorylation of Ser-180. The present study has been carried out to understand the structural mechanism of regulation of water permeability across the channel. Molecular dynamics (MD) simulations have been carried out to investigate the structural changes caused due to phosphorylation of Ser-180 in the tetrameric assembly of rAQP4 along with predicted C-terminal region (255-323). The interactions involving opposite charges are observed between cytoplasmic loops and the C-terminal region during MD simulations. This results in movement of C-terminal region of rAQP4 towards the cytoplasmic mouth of water channel. Despite this movement, there was a gap between C-terminal region and cytoplasmic mouth of the channel through which water molecules were able to gain entry into the channel. The interactions between C-terminus and loop D of neighboring monomers in a tetrameric assembly appear to prevent the complete closure of cytoplasmic mouth of the water channel. Further, the rates of water permeation through phosphorylated and unphosphorylated rAQP4 have also been compared. The simulation studies showed a continuous movement of water in a single file across pore of unphosphorylated as well as phosphorylated rAQP4.

Molecular characterization of branchial aquaporin 1aa and effects of seawater acclimation, emersion or ammonia exposure on its mRNA expression in the gills, gut, kidney and skin of the freshwater climbing perch, Anabas testudineus

We obtained a full cDNA coding sequence of aquaporin 1aa (aqp1aa) from the gills of the freshwater climbing perch, Anabas testudineus, which had the highest expression in the gills and skin, suggesting an important role of Aqp1aa in these organs. Since seawater acclimation had no significant effects on the branchial and intestinal aqp1aa mRNA expression, and since the mRNA expression of aqp1aa in the gut was extremely low, it can be deduced that Aqp1aa, despite being a water channel, did not play a significant osmoregulatory role in A. testudineus. However, terrestrial exposure led to significant increases in the mRNA expression of aqp1aa in the gills and skin of A. testudineus. Since terrestrial exposure would lead to evaporative water loss, these results further support the proposition that Aqp1aa did not function predominantly for the permeation of water through the gills and skin. Rather, increased aqp1aa mRNA expression might be necessary to facilitate increased ammonia excretion during emersion, because A. testudineus is known to utilize amino acids as energy sources for locomotor activity with increased ammonia production on land. Furthermore, ammonia exposure resulted in significant decreases in mRNA expression of aqp1aa in the gills and skin of A. testudineus, presumably to reduce ammonia influx during ammonia loading. This corroborates previous reports on AQP1 being able to facilitate ammonia permeation. However, a molecular characterization of Aqp1aa from A. testudineus revealed that its intrinsic aquapore might not facilitate NH₃ transport. Hence, ammonia probably permeated the central fifth pore of the Aqp1aa tetramer as suggested previously. Taken together, our results indicate that Aqp1aa might have a greater physiological role in ammonia excretion than in osmoregulation in A. testudineus.

In silico study of human aquaporin AQP11 and AQP12 channels

AQP11 and AQP12 are the most distantly related paralogs of the aquaporin family in human. They share indeed a low sequence similarity with other aquaporins and exhibit a modified N-terminal NPA signature motif. Furthermore, they have an anomalous subcellular localization. The AQP11 and AQP12 biological role remains to be fully clarified and their ability to allow transport of water is still debated. We have built accurate 3D-models for AQP11 and AQP12 and comprehensively compared their sequence and structure to other known aquaporins. In order to investigate whether they appear compatible or not with water permeability, we especially focused on the amino acid composition and electrostatics of their channels, keeping the structure of the low-water efficiency AQP0 as a reference system. Our analysis points out a possible alternative ar/R site and shows that these aquaporins feature unique residues at key pore-lining positions that make the shape, composition and electrostatics of their channel peculiar. Such residues can represent pivotal hints to study and explain the AQP11 and AQP12 biological and molecular function.

Multiscale molecular dynamics simulations of membrane proteins

The time and length scales accessible by biomolecular simulations continue to increase. This is in part due to improvements in algorithms and computing performance, but is also the result of the emergence of coarse-grained (CG) potentials, which complement and extend the information obtainable from fully detailed models. CG methods have already proven successful for a range of applications that benefit from the ability to rapidly simulate spontaneous self-assembly within a lipid membrane environment, including the insertion and/or oligomerization of a range of "toy models," transmembrane peptides, and single- and multi-domain proteins. While these simplified approaches sacrifice atomistic level detail, it is now straightforward to "reverse map" from CG to atomistic descriptions, providing a strategy to assemble membrane proteins within a lipid environment, prior to all-atom simulation. Moreover, recent developments have been made in "dual resolution" techniques, allowing different molecules in the system to be modeled with atomistic or CG resolution simultaneously.

Intestinal receptor targeting for peptide delivery: an expert's personal perspective on reasons for failure and new opportunities

The technology has been available more than 25 years that would enable the oral delivery of vaccines, proteins and peptides, thus avoiding the need for injection. To this day, injection is still the mode of delivery, yet not the main mode of choice. This review focuses on several of the potential modes for oral delivery of peptides, proteins and vaccines. Additionally, the review will provide the reader with an insight into the problems and potential solutions for several of these modes of oral delivery of peptides and proteins.

Structure and permeation mechanism of a mammalian urea transporter

As an adaptation to infrequent access to water, terrestrial mammals produce urine that is hyperosmotic to plasma. To prevent osmotic diuresis by the large quantity of urea generated by protein catabolism, the kidney epithelia contain facilitative urea transporters (UTs) that allow rapid equilibration between the urinary space and the hyperosmotic interstitium. Here we report the first X-ray crystal structure of a mammalian UT, UT-B, at a resolution of 2.36 Å. UT-B is a homotrimer and each protomer contains a urea conduction pore with a narrow selectivity filter. Structural analyses and molecular dynamics simulations showed that the selectivity filter has two urea binding sites separated by an approximately 5.0 kcal/mol energy barrier. Functional studies showed that the rate of urea conduction in UT-B is increased by hypoosmotic stress, and that the site of osmoregulation coincides with the location of the energy barrier.

Insights into the mechanisms of the selectivity filter of Escherichia coli aquaporin Z

Aquaporin Z (AQPZ) is a tetrameric protein that forms water channels in the cell membrane of Escherichia coli. The histidine residue (residue 174) in the selectivity filter (SF) region plays an important role in the transport of water across the membrane. In this work, we perform equilibrium molecular dynamics (MD) simulations to illustrate the gating mechanism of the SF and the influences of residue 174 in two different protonation states: Hsd174 with the proton at Nδ, and Hse174 with the proton at Nε. We calculate the pore radii in the SF region versus the simulation time. We perform steered MD to compute the free-energy profile, i.e., the potential of mean force (PMF) of a water molecule through the SF region. We conduct a quantum mechanics calculation of the binding energy of one water molecule with the residues in the SF region. The hydrogen bonds formed between the side chain of Hsd174 and the side chain of residue 189 (Arg189) play important roles in the selectivity mechanism of AQPZ. The radii of the pores, the hydrogen-bond analysis, and the free energies show that it is easier for water molecules to permeate through the SF region of AQPZ with residue 174 in the Hse state than in the Hsd state.

Molecular wire of urea in carbon nanotube: a molecular dynamics study

We perform molecular dynamics simulations of narrow single-walled carbon nanotubes (SWNTs) in aqueous urea to investigate the structure and dynamical behavior of urea molecules inside the SWNT. Even at low urea concentrations (e.g., 0.5 M), we have observed spontaneous and continuous filling of SWNT with a one-dimensional urea wire (leaving very few water molecules inside the SWNT). The urea wire is structurally ordered, both translationally and orientationally, with a contiguous hydrogen-bonded network and concerted urea's dipole orientations. Interestingly, despite the symmetric nature of the whole system, the potential energy profile of urea along the SWNT is asymmetric, arising from the ordering of asymmetric urea partial charge distribution (or dipole moment) in confined environment. Furthermore, we study the kinetics of confined urea and find that the permeation of urea molecules through the SWNT decreases significantly (by a factor of ∼20) compared to that of water molecules, due to the stronger dispersion interaction of urea with SWNT than water, and a maximum in urea permeation happens around a concentration of 5 M. These findings might shed some light on the better understanding of unique properties of molecular wires (particularly the wires formed by polar organic small molecules) confined within both artificial and biological nanochannels, and are expected to have practical applications such as the electronic devices for signal transduction and multiplication at the nanoscale.

Exploring substrate diffusion in channels using biased molecular dynamics simulations

Substrate transport and diffusion through membrane-bound channels are processes that can span a range of time scales, with only the fastest ones being amenable to most atomic-scale equilibrium molecular dynamics (MD) simulations. However, the application of forces within a simulation can greatly accelerate diffusion processes, revealing important structural and energetic features of the channel. Here, we demonstrate the use of two methods for applying biases to a substrate in a simulation, using the ammonia/ammonium transporter AmtB as an example. The first method, steered MD, applies a constant force or velocity constraint to the substrate, permitting the exploration of potential substrate pathways and the barriers encountered, although typically far outside of equilibrium. On the other hand, the second method, adaptive biasing forces, is quasi-equilibrium, permitting the derivation of a potential of mean force, which characterizes the free energy of the substrate during transport.

Conformational dynamics of a membrane transport protein probed by H/D exchange and covalent labeling: the glycerol facilitator

Glycerol facilitator (GF) is a tetrameric membrane protein responsible for the selective permeation of glycerol and water. Each of the four GF subunits forms a transmembrane channel. Every subunit consists of six helices that completely span the lipid bilayer, as well as two half-helices (TM7 and TM3). X-ray crystallography has revealed that the selectivity of GF is due to its unique amphipathic channel interior. To explore the structural dynamics of GF, we employ hydrogen/deuterium exchange (HDX) and oxidative labeling with mass spectrometry (MS). HDX-MS reveals that transmembrane helices are generally more protected than extramembrane segments, consistent with data previously obtained for other membrane proteins. Interestingly, TM7 does not follow this trend. Instead, this half-helix undergoes rapid deuteration, indicative of a highly dynamic local structure. The oxidative labeling behavior of most GF residues is consistent with the static crystal structure. However, the side chains of C134 and M237 undergo labeling although they should be inaccessible according to the X-ray structure. In agreement with our HDX-MS data, this observation attests to the fact that TM7 is only marginally stable. We propose that the highly mobile nature of TM7 aids in the efficient diffusion of guest molecules through the channel ("molecular lubrication"). In the absence of such dynamics, host-guest molecular recognition would favor semipermanent binding of molecules inside the channel, thereby impeding transport. The current work highlights the complementary nature of HDX, covalent labeling, and X-ray crystallography for the characterization of membrane proteins.

The aromatic/arginine selectivity filter of NIP aquaporins plays a critical role in substrate selectivity for silicon, boron, and arsenic

Nodulin-26-like intrinsic proteins (NIPs) of the aquaporin family are involved in the transport of diverse solutes, but the mechanisms controlling the selectivity of transport substrates are poorly understood. The purpose of this study was to investigate how the aromatic/arginine (ar/R) selectivity filter influences the substrate selectivity of two NIP aquaporins; the silicic acid (Si) transporter OsLsi1 (OsNIP2;1) from rice and the boric acid (B) transporter AtNIP5;1 from Arabidopsis; both proteins are also permeable to arsenite. Native and site-directed mutagenized variants of the two genes were expressed in Xenopus oocytes and the transport activities for Si, B, arsenite, and water were assayed. Substitution of the amino acid at the ar/R second helix (H2) position of OsLsi1 did not affect the transport activities for Si, B, and arsenite, but that at the H5 position resulted in a total loss of Si and B transport activities and a partial loss of arsenite transport activity. Conversely, changes of the AtNIP5;1 ar/R selectivity filter and the NPA motifs to the OsLsi1 type did not result in a gain of Si transport activity. B transport activity was partially lost in the H5 mutant but unaffected in the H2 mutant of AtNIP5;1. In contrast, both the single and double mutations at the H2 and/or H5 positions of AtNIP5;1 did not affect arsenite transport activity. The results reveal that the residue at the H5 position of the ar/R filter of both OsLsi1 and AtNIP5;1 plays a key role in the permeability to Si and B, but there is a relatively low selectivity for arsenite.

Arabidopsis thaliana NIP7;1: an anther-specific boric acid transporter of the aquaporin superfamily regulated by an unusual tyrosine in helix 2 of the transport pore

Plant nodulin-26 intrinsic proteins (NIPs) are members of the aquaporin superfamily that serve as multifunctional transporters of uncharged metabolites. In Arabidopsis thaliana, a specific NIP pore subclass, known as the NIP II proteins, is represented by AtNIP5;1 and AtNIP6;1, which encode channel proteins expressed in roots and leaf nodes, respectively, that participate in the transport of the critical cell wall nutrient boric acid. Modeling of the protein encoded by the AtNIP7;1 gene shows that it is a third member of the NIP II pore subclass in Arabidopsis. However, unlike AtNIP5;1 and AtNIP6;1 proteins, which form constitutive boric acid channels, AtNIP7;1 forms a channel with an extremely low intrinsic boric acid transport activity. Molecular modeling and molecular dynamics simulations of AtNIP7;1 suggest that a conserved tyrosine residue (Tyr81) located in transmembrane helix 2 adjacent to the aromatic arginine (ar/R) pore selectivity region stabilizes a closed pore conformation through interaction with the canonical Arg220 in ar/R region. Substitution of Tyr81 with a Cys residue, characteristic of established NIP boric acid channels, results in opening of the AtNIP7;1 pore that acquires a robust, transport activity for boric acid as well as other NIP II test solutes (glycerol and urea). Substitution of a Phe for Tyr81 also opens the channel, supporting the prediction from MD simulations that hydrogen bond interaction between the Tyr81 phenol group and the ar/R Arg may contribute to the stabilization of a closed pore state. Expression analyses show that AtNIP7;1 is selectively expressed in developing anther tissues of young floral buds of A. thaliana, principally in developing pollen grains of stage 9-11 anthers. Because boric acid is both an essential nutrient as well as a toxic compound at high concentrations, it is proposed that Tyr81 modulates transport and may provide an additional level of regulation for this transporter in male gametophyte development.

Evagination of cells controls bio-silica formation and maturation during spicule formation in sponges

The enzymatic-silicatein mediated formation of the skeletal elements, the spicules of siliceous sponges starts intracellularly and is completed extracellularly. With Suberites domuncula we show that the axial growth of the spicules proceeds in three phases: (I) formation of an axial canal; (II) evagination of a cell process into the axial canal, and (III) assembly of the axial filament composed of silicatein. During these phases the core part of the spicule is synthesized. Silicatein and its substrate silicate are stored in silicasomes, found both inside and outside of the cellular extension within the axial canal, as well as all around the spicule. The membranes of the silicasomes are interspersed by pores of ≈2 nm that are likely associated with aquaporin channels which are implicated in the hardening of the initial bio-silica products formed by silicatein. We can summarize the sequence of events that govern spicule formation as follows: differential genetic readout (of silicatein) → fractal association of the silicateins → evagination of cells by hydro-mechanical forces into the axial canal → and finally processive bio-silica polycondensation around the axial canal. We termed this process, occurring sequentially or in parallel, bio-inorganic self-organization.

Probing osmolyte participation in the unfolding transition state of a protein

Understanding the molecular mechanisms of osmolyte protection in protein stability has proved to be challenging. In particular, little is known about the role of osmolytes in the structure of the unfolding transition state of a protein, the main determinant of its dynamics. We have developed an experimental protocol to directly probe the transition state of a protein in a range of osmolyte environments. We use an atomic force microscope in force-clamp mode to apply mechanical forces to the protein I27 and obtain force-dependent rate constants of protein unfolding. We measure the distance to the unfolding transition state, Δx(u), along a 1D reaction coordinate imposed by mechanical force. We find that for the small osmolytes, ethylene glycol, propylene glycol, and glycerol, Δx(u) scales with the size of the molecule, whereas for larger osmolytes, sorbitol and sucrose, Δx(u) remains the same as that measured in water. These results are in agreement with steered molecular dynamics simulations that show that small osmolytes act as solvent bridges in the unfolding transition state structure, whereas only water molecules act as solvent bridges in large osmolyte environments. These results demonstrate that novel force protocols combined with solvent substitution can directly probe angstrom changes in unfolding transition state structure. This approach creates new opportunities to gain molecular level understanding of the action of osmolytes in biomolecular processes.

Voltage-regulated water flux through aquaporin channels in silico

Aquaporins (AQPs) facilitate the passive flux of water across biological membranes in response to an osmotic pressure. A number of AQPs, for instance in plants and yeast, have been proposed to be regulated by phosphorylation, cation concentration, pH change, or membrane-mediated mechanical stress. Here we report an extensive set of molecular dynamics simulations of AQP1 and AQP4 subject to large membrane potentials in the range of ±1.5 V, suggesting that AQPs may in addition be regulated by an electrostatic potential. As the regulatory mechanism we identified the relative population of two different states of the conserved arginine in the aromatic/arginine constriction region. A positive membrane potential was found to stabilize the arginine in an up-state, which allows rapid water flux, whereas a negative potential favors a down-state, which reduces the single-channel water permeability.

Functional characterization and hyperosmotic regulation of aquaporin in Synechocystis sp. PCC 6803

The genome of the cyanobacterium Synechocystis sp. PCC 6803 (hereafter, Synechocystis) contains an aqpZ gene (slr2057) which encodes an aquaporin (SsAqpZ), a membrane channel protein that might play a role in osmotic water transport and therefore the growth of Synechocystis. Structural characterization of SsAqpZ by protein sequence analysis and homology modelling revealed that it was more similar to bacterial aquaporin Z than the glycerol facilitator. To understand the functional role of SsAqpZ, the aqpZ knockout (KO) and myc-tagged aqpZ knockin (KI) Synechocystis were constructed. Water channel activity assays showed that SsAqpZ facilitated water transportation. SsAqpZ-mediated changes in cell volume were observed in wild-type (WT) and KI Synechocystis. Expression of SsAqpZ in KI Synechocystis was induced by extracellular hyperosmolarity. In the absence of hyperosmolarity, WT, KO and KI Synechocystis showed the same pattern of growth and no morphological or phenotypical perturbations. Under hyperosmotic condition, while the WT and also KI cells maintained a similar growth rate throughout the entire exponential phase, KO cells grew significantly slower. These results indicate that SsAqpZ has water channel activity and is involved in the adaptation and maintenance of growth of Synechocystis in a hyperosmotic environment.

Exploring transmembrane diffusion pathways with molecular dynamics

Transmembrane exchange of materials is a fundamental process in biology. Molecular dynamics provides a powerful method to investigate in great detail various aspects of the phenomenon, particularly the permeation of small uncharged molecules, which continues to pose a challenge to experimental studies. We will discuss some of the recent simulation studies investigating the role of lipid-mediated and protein-mediated mechanisms in permeation of water and gas molecules across the membrane.

Structural context shapes the aquaporin selectivity filter

Aquaporins are transmembrane channels that facilitate the permeation of water and small, uncharged amphipathic molecules across cellular membranes. One distinct aquaporin subfamily contains pure water channels, whereas a second subfamily contains channels that conduct small alditols such as glycerol, in addition to water. Distinction between these substrates is central to aquaporin function, though the contributions of protein structural motifs required for selectivity are not yet fully characterized. To address this question, we sequentially engineered three signature amino acids of the glycerol-conducting subfamily into the Escherichia coli water channel aquaporin Z (AqpZ). Functional analysis of these mutant channels showed a decrease in water permeability but not the expected increase in glycerol conduction. Using X-ray crystallography, we determined the atomic resolution structures of the mutant channels. The structures revealed a channel surprisingly similar in size to the wild-type AqpZ pore. Comparison with measured rates of transport showed that, as the size of the selectivity filter region of the channel approaches that of water, channel hydrophilicity dominated water conduction energetics. In contrast, the major determinant of selectivity for larger amphipathic molecules such as glycerol was channel cross-section size. Finally, we find that, although the selectivity filter region is indeed central to substrate transport, other structural elements that do not directly interact with the substrates, such as the loop connecting helices M6 and M7, and the C loop between helices C4 and C5, play an essential role in facilitating selectivity.

Arsenic transport in prokaryotes and eukaryotic microbes

Aquaporins (AQPs) and aquaglyceroporins facilitate transport of a broad spectrum of substrates such as water, glycerol and other small uncharged solutes. More recently, AQPs ave also been shown to facilitate diffusion of metalloids such as arsenic (As) and antimony (Sb). At neutral pH, the trivalent forms of these metalloids are structurally similar to glycerol and hence they can enter cells through AQPs. As- and Sb-containing compounds are toxic to cells, yet both metalloids are used as chemotherapeutic agents for treating acute promyelocytic leukemia and diseases caused by protozoan parasites. In this chapter, we will review the role of AQPs and other proteins in metalloid transport in prokaryotes and eukaryotic microbes.

Protein dynamics and ligand migration interplay as studied by computer simulation

Since proteins are dynamic systems in living organisms, the employment of methodologies contemplating this crucial characteristic results fundamental to allow revealing several aspects of their function. In this work, we present results obtained using classical mechanical atomistic simulation tools applied to understand the connection between protein dynamics and ligand migration. Firstly, we will present a review of the different sampling schemes used in the last years to obtain both ligand migration pathways and the thermodynamic information associated with the process. Secondly, we will focus on representative examples in which the schemes previously presented are employed, concerning the following: i) ligand migration, tunnels, and cavities in myoglobin and neuroglobin; ii) ligand migration in truncated hemoglobin members; iii) NO escape and conformational changes in nitrophorins; iv) ligand selectivity in catalase and hydrogenase; and v) larger ligand migration: the P450 and haloalkane dehalogenase cases. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Involvement of silicon influx transporter OsNIP2;1 in selenite uptake in rice

Rice (Oryza sativa) as a staple food, provides a major source of dietary selenium (Se) for humans, which essentially requires Se, however, the molecular mechanism for Se uptake is still poorly understood. Herein, we show evidence that the uptake of selenite, a main bioavailable form of Se in paddy soils, is mediated by a silicon (Si) influx transporter Lsi1 (OsNIP2;1) in rice. Defect of OsNIP2;1 resulted in a significant decrease in the Se concentration of the shoots and xylem sap when selenite was given. However, there was no difference in the Se concentration between the wild-type rice and mutant of OsNIP2;1 when selenate was supplied. A short-term uptake experiment showed that selenite uptake greatly increased with decreasing pH in the external solution. Si as silicic acid did not inhibit the Se uptake from selenite in both rice and yeast (Saccharomyces cerevisiae) at low pHs. Expression of OsNIP2;1 in yeast enhanced the selenite uptake at pH 3.5 and 5.5 but not at pH 7.5. On the other hand, defect of Si efflux transporter Lsi2 did not affect the uptake of Se either from selenite or selenate. Taken together, our results indicate that Si influx transporter OsNIP2;1 is permeable to selenite.

Deciphering ionic current signatures of DNA transport through a nanopore

Within just a decade from the pioneering work demonstrating the utility of nanopores for molecular sensing, nanopores have emerged as versatile systems for single-molecule manipulation and analysis. In a typical setup, a gradient of the electrostatic potential captures charged solutes from the solution and forces them to move through a single nanopore, across an otherwise impermeable membrane. The ionic current blockades resulting from the presence of a solute in a nanopore can reveal the type of the solute, for example, the nucleotide makeup of a DNA strand. Despite great success, the microscopic mechanisms underlying the functionality of such stochastic sensors remain largely unknown, as it is not currently possible to characterize the microscopic conformations of single biomolecules directly in a nanopore and thereby unequivocally establish the causal relationship between the observables and the microscopic events. Such a relationship can be determined using molecular dynamics-a computational method that can accurately predict the time evolution of a molecular system starting from a given microscopic state. This article describes recent applications of this method to the process of DNA transport through biological and synthetic nanopores.

Electrostatics of aquaporin and aquaglyceroporin channels correlates with their transport selectivity

Aquaporins are homotetrameric channel proteins, which allow the diffusion of water and small solutes across biological membranes. According to their transport function, aquaporins can be divided into "orthodox aquaporins", which allow the flux of water molecules only, and "aquaglyceroporins", which facilitate the diffusion of glycerol and other small solutes in addition to water. The contribution of individual residues in the pore to the selectivity of orthodox aquaporins and aquaglyceroporins is not yet fully understood. To gain insights into aquaporin selectivity, we focused on the sequence variation and electrostatics of their channels. The continuum Poisson-Boltzmann electrostatic potential along the channel was calculated and compared for ten three-dimensional-structures which are representatives of different aquaporin subfamilies, and a panel of functionally characterized mutants, for which high-accuracy three-dimensional-models could be derived. Interestingly, specific electrostatic profiles associated with the main selectivity to water or glycerol could be identified. In particular: (i) orthodox aquaporins showed a distinctive electrostatic potential maximum at the periplasmic side of the channel around the aromatic/Arg (ar/R) constriction site; (ii) aquaporin-0 (AQP0), a mammalian aquaporin with considerably low water permeability, had an additional deep minimum at the cytoplasmic side; (iii) aquaglyceroporins showed a rather flat potential all along the channel; and (iv) the bifunctional protozoan PfAQP had an unusual all negative profile. Evaluation of electrostatics of the mutants, along with a thorough sequence analysis of the aquaporin pore-lining residues, illuminated the contribution of specific residues to the electrostatics of the channels and possibly to their selectivity.

A single glutamate residue controls the oligomerization, function, and stability of the aquaglyceroporin GlpF

Like many other alpha-helical membrane proteins, the monomeric Escherichia coli aquaglyceroporin GlpF associates within cellular membranes and forms higher-order oligomeric structures. A potential impact of the oligomeric state on the protein function remains enigmatic. We have analyzed the role of residues W42 and E43 in the oligomerization of the E. coli GlpF protein in vitro and in vivo. In contrast to W42, the polar glutamate residue at position 43 appears to be critical for oligomerization. While other polar residues can substitute for the function of E43, replacement of E43 with alanine results in a greatly reduced GlpF oligomerization propensity. The reduced interaction propensity of GlpF E43A correlates with an impaired in vivo function as well as a decreased in vivo stability. Therefore, E43 is critical for the proper oligomerization of GlpF, and protein oligomerization appears to be crucial for the channel function as well as for the in vivo stability of the protein.

Genome-wide analysis of major intrinsic proteins in the tree plant Populus trichocarpa: characterization of XIP subfamily of aquaporins from evolutionary perspective

BACKGROUND

Members of major intrinsic proteins (MIPs) include water-conducting aquaporins and glycerol-transporting aquaglyceroporins. MIPs play important role in plant-water relations. The model plants Arabidopsis thaliana, rice and maize contain more than 30 MIPs and based on phylogenetic analysis they can be divided into at least four subfamilies. Populus trichocarpa is a model tree species and provides an opportunity to investigate several tree-specific traits. In this study, we have investigated Populus MIPs (PtMIPs) and compared them with their counterparts in Arabidopsis, rice and maize.

RESULTS

Fifty five full-length MIPs have been identified in Populus genome. Phylogenetic analysis reveals that Populus has a fifth uncharacterized subfamily (XIPs). Three-dimensional models of all 55 PtMIPs were constructed using homology modeling technique. Aromatic/arginine (ar/R) selectivity filters, characteristics of loops responsible for solute selectivity (loop C) and gating (loop D) and group conservation of small and weakly polar interfacial residues have been analyzed. Majority of the non-XIP PtMIPs are similar to those in Arabidopsis, rice and maize. Additional XIPs were identified from database search and 35 XIP sequences from dicots, fungi, moss and protozoa were analyzed. Ar/R selectivity filters of dicots XIPs are more hydrophobic compared to fungi and moss XIPs and hence they are likely to transport hydrophobic solutes. Loop C is longer in one of the subgroups of dicot XIPs and most probably has a significant role in solute selectivity. Loop D in dicot XIPs has higher number of basic residues. Intron loss is observed on two occasions: once between two subfamilies of eudicots and monocot and in the second instance, when dicot and moss XIPs diverged from fungi. Expression analysis of Populus MIPs indicates that Populus XIPs don't show any tissue-specific transcript abundance.

CONCLUSION

Due to whole genome duplication, Populus has the largest number of MIPs identified in any single species. Non-XIP MIPs are similar in all four plant species considered in this study. Small and weakly polar residues at the helix-helix interface are group conserved presumably to maintain the hourglass fold of MIP channels. Substitutions in ar/R selectivity filter, insertion/deletion in loop C, increasing basic nature of loop D and loss of introns are some of the events occurred during the evolution of dicot XIPs.

Insight into factors directing high production of eukaryotic membrane proteins; production of 13 human AQPs in Pichia pastoris

Membrane proteins are key players in all living cells. To achieve a better understanding of membrane protein function, significant amounts of purified protein are required for functional and structural analyses. Overproduction of eukaryotic membrane proteins, in particular, is thus an essential yet non-trivial task. Hence, improved understanding of factors which direct a high production of eukaryotic membrane proteins is desirable. In this study we have compared the overproduction of all human aquaporins in the eukaryotic host Pichia pastoris. We report quantitated production levels of each homologue and the extent of their membrane localization. Our results show that the protein production levels vary substantially, even between highly homologous aquaporins. A correlation between the extents of membrane insertion with protein function also emerged, with a higher extent of membrane insertion for pure water transporters compared to aquaporin family members with other substrate specificity. Nevertheless, the nucleic acid sequence of the second codon appears to play an important role in overproduction. Constructs containing guanine at the first position of this codon (being part of the mammalian Kozak sequence) are generally produced at a higher level, which is confirmed for hAQP8. In addition, mimicking the yeast consensus sequence (ATGTCT) apparently has a negative influence on the production level, as shown for hAQP1. Moreover, by mutational analysis we show that the yield of hAQP4 can be heavily improved by directing the protein folding pathway as well as stabilizing the aquaporin tetramer.

Water channel proteins (later called aquaporins) and relatives: past, present, and future

Water channels or water channel proteins (WCPs) are transmembrane proteins that have a specific three-dimensional structure with a pore that can be permeated by water molecules. WCPs are large families (over 450 members) that are present in all kingdoms of life. The first WCP was discovered in the human red blood cell (RBC) membrane in 1980s. In 1990s other WCPs were discovered in plants, microorganisms, various animals, and humans; and it became obvious that the WCPs belong to the superfamily of major intrinsic proteins (MIPs, over 800 members). WCPs include three subfamilies: (a) aquaporins (AQPs), which are water specific (or selective water channels); (b) aquaglyceroporins (and glycerol facilitators), which are permeable to water and/or other small molecules; and (c) "superaquaporins" or subcellular AQPs. WCPs (and MIPs) have several structural characteristics which were better understood after the atomic structure of some MIPs was deciphered. The structure-function relationships of MIPs expressed in microorganisms (bacteria, archaea, yeast, and protozoa), plants, and some multicellular animal species [nematodes, insects, fishes, amphibians, mammals (and humans)] are described. A synthetic overview on the WCPs from RBCs from various species is provided. The physiological roles of WCPs in kidney, gastrointestinal system, respiratory apparatus, central nervous system, eye, adipose tissue, skin are described, and some implications of WCPs in various diseases are briefly presented. References of detailed reviews on each topic are given. This is the first review providing in a condensed form an overview of the whole WCP field that became in the last 20 years a very hot area of research in biochemistry and molecular cell biology, with wide and increasing implications.

Dynamics and energetics of permeation through aquaporins. What do we learn from molecular dynamics simulations?

Aquaporins (AQPs) are a family of integral membrane proteins, which facilitate the rapid and yet highly selective flux of water and other small solutes across biological membranes. Molecular dynamics (MD) simulations contributed substantially to the understanding of the molecular mechanisms that underlie this remarkable efficiency and selectivity of aquaporin channels. This chapter reviews the current state of MD simulations of aquaporins and related aquaglyceroporins as well as the insights these simulations have provided. The mechanism of water permeation through AQPs and methods to determine channel permeabilities from simulations are described. Protons are strictly excluded from AQPs by a large electrostatic barrier and not by an interruption of the Grotthuss mechanism inside the pore. Both the protein's electric field and desolvation effects contribute to this barrier. Permeation of apolar gas molecules such as CO(2) through AQPs is accompanied by a large energetic barrier and thus can only be expected in membranes with a low intrinsic gas permeability. Additionally, the insights from simulations into the mechanism of glycerol permeation through the glycerol facilitator GlpF from E. coli are summarized. Finally, MD simulations are discussed that revealed that the aro-matic/arginine constriction region is generally the filter for uncharged solutes, and that AQP selectivity is controlled by a hydrophobic effect and steric restraints.

Dynamic control of slow water transport by aquaporin 0: implications for hydration and junction stability in the eye lens

Aquaporin 0 (AQP0), the most abundant membrane protein in mammalian lens fiber cells, not only serves as the primary water channel in this tissue but also appears to mediate the formation of thin junctions between fiber cells. AQP0 is remarkably less water permeable than other aquaporins, but the structural basis and biological significance of this low permeability remain uncertain, as does the permeability of the protein in a reported junctional form. To address these issues, we performed molecular dynamics (MD) simulations of water transport through membrane-embedded AQP0 in both its (octameric) junctional and (tetrameric) nonjunctional forms. From our simulations, we measured an osmotic permeability for the nonjunctional form that agrees with experiment and found that the distinct dynamics of the conserved, lumen-protruding side chains of Tyr-23 and Tyr-149 modulate water passage, accounting for the slow permeation. The junctional and nonjunctional forms conducted water equivalently, in contrast to a previous suggestion based on static crystal structures that water conduction is lost on junction formation. Our analysis suggests that the low water permeability of AQP0 may help maintain the mechanical stability of the junction. We hypothesize that the structural features leading to low permeability may have evolved in part to allow AQP0 to form junctions that both conduct water and contribute to the organizational structure of the fiber cell tissue and microcirculation within it, as required to maintain transparency of the lens.

Side-chain dynamics are critical for water permeation through aquaporin-1

Molecular dynamics simulations of aquaporin-1 embedded in a solvated lipid bilayer were carried out to investigate the mechanism of water permeation. The 2.2 A resolution crystal structure of the bovine protein was used for five independent trajectories. During the equilibration and preparatory steps in which the protein was held fixed, water molecules inside the water channel adopted the same positions as observed in the crystal structure but they did not pass through the channel, suggesting that the dynamic motion of the protein is critical for water permeation. When the protein atoms were allowed to move, the side chains of the two asparagines in the two conserved Asn-Pro-Ala motifs near the center of the channel formed hydrogen bonds with water and helped water molecules move through the channel by actively aligning them for transport. The main-chain oxygen atoms, which were exposed to the pore surface in the crystal structure, also contributed to water transfer. Besides the constriction region observed in the crystal structure (Arg(197), Phe(58), His(182), and Cys(191)), we found that His(76) and Val(155) act as a valve by dynamically blocking water permeation and helping control flow.

Electrostatic funneling of substrate in mitochondrial inner membrane carriers

Exchange of ATP and ADP across mitochondrial membrane replenishes the cytoplasm with newly synthesized ATP and provides the mitochondria with the substrate ADP for oxidative phosphorylation. The sole means of this exchange is the mitochondrial ADP/ATP carrier (AAC), a membrane protein that is suggested to cycle between two conformationally distinct states, cytosolic-open (c-state) and matrix-open (m-state), thereby shuttling nucleotides across the inner mitochondrial membrane. However, the c-state is the only structurally resolved state, and the binding site of ADP remains elusive. Here, we present approximately 0.3 mus of all-atom MD simulations of the c-state revealing rapid, spontaneous binding of ADP to deeply positioned binding sites within the AAC lumen. To our knowledge, a complete ligand-binding event has heretofore not been described in full atomic detail in unbiased simulations. The identified ADP-bound state and additional simulations shed light on key structural elements and the initial steps involved in conversion to the m-state. Electrostatic analysis of trajectories reveals the presence of an unusually strong positive electrostatic potential in the lumen of AAC that appears to be the main driving force for the observed spontaneous binding of ADP. We provide evidence that the positive electrostatic potential is likely a common attribute among the entire family of mitochondrial carriers. In addition to playing a key role in substrate recruitment and translocation, the electropositivity of mitochondrial carriers might also be critical for their binding to the negatively charged environment of the inner mitochondrial membrane.

Free energy calculations applied to membrane proteins

Selected applications of free energy calculations to the realm of membrane proteins are reviewed. The theoretical underpinnings of these calculations are described, focusing on free energy perturbation and the use of thermodynamic integration to determine free energy changes along well-delineated order parameters. Current strategies for improving the reliability of free energy calculations, while making them somewhat more affordable are outlined. Application of the free energy methodology to understand the structure and function of membrane proteins is illustrated in three concrete examples: The binding of an agonist ligand to a G protein-coupled receptor, the assisted transport of a small permeant through a membrane channel, and the recognition and association of transmembrane alpha-helical domains.

Crystal structure of the aquaglyceroporin PfAQP from the malarial parasite Plasmodium falciparum

The 2.05-Å resolution structure of the aquaglyceroporin from the malarial parasite Plasmodium falciparum (PfAQP), a protein important in the parasite’s life cycle, has been solved. The structure provides key evidence for the basis of water versus glycerol selectivity in aquaporin family members. Unlike its closest homolog of known structure, GlpF, the channel conducts both glycerol and water at high rates, framing the question of what determines high water conductance in aquaporin channels. The universally conserved arginine in the selectivity filter is constrained by only two hydrogen bonds in GlpF, whereas there are three in all water-selective aquaporins and in PfAQP. The decreased cost of dehydrating the triply-satisfied arginine cation may provide the basis for high water conductance. The two Asn-Pro-Ala (NPA) regions of PfAQP, which bear rare substitutions to Asn-Leu-Ala (NLA) and Asn-Pro-Ser (NPS), participate in preserving the orientation of the selectivity filter asparagines in the center of the channel.

Solvent molecules bridge the mechanical unfolding transition state of a protein

We demonstrate a combination of single molecule force spectroscopy and solvent substitution that captures the presence of solvent molecules in the transition state structure. We measure the effect of solvent substitution on the rate of unfolding of the I27 titin module, placed under a constant stretching force. From the force dependency of the unfolding rate, we determine Deltax(u), the distance to the transition state. Unfolding the I27 protein in water gives a Deltax(u) = 2.5 A, a distance that compares well to the size of a water molecule. Although the height of the activation energy barrier to unfolding is greatly increased in both glycerol and deuterium oxide solutions, Deltax(u) depends on the size of the solvent molecules. Upon replacement of water by increasing amounts of the larger glycerol molecules, Deltax(u) increases rapidly and plateaus at its maximum value of 4.4 A. In contrast, replacement of water by the similarly sized deuterium oxide does not change the value of Deltax(u). From these results we estimate that six to eight water molecules form part of the unfolding transition state structure of the I27 protein, and that the presence of just one glycerol molecule in the transition state is enough to lengthen Deltax(u). Our results show that solvent composition is important for the mechanical function of proteins. Furthermore, given that solvent composition is actively regulated in vivo, it may represent an important modulatory pathway for the regulation of tissue elasticity and other important functions in cellular mechanics.