Molecular dynamics simulation and bioinformatics study on yeast aquaporin Aqy1 from Pichia pastoris

In the present study, an equilibrated system for the Aqy1 tetramer was developed, and molecular biophysics modeling showed that the Aqy1 channel was blocked by Tyr-31 in the N-terminus, which was also supported by the free energy profiles. However, bioinformatics analysis of the amino acid sequence of Aqy1 indicated this Tyr-31 is not conserved across all fungi. Analysis of the equilibrated structure showed that the central pore along the four-fold axis of the tetramers is formed with hydrophobic amino acid residues. In particular, Phe-90, Trp-198, and Phe-202 form the narrowest part of the pore. Therefore, water molecules are not expected to translocate through the central pore, a hypothesis that we confirmed by molecular dynamics simulations. Each monomer of the Aqy1 tetramers forms a channel whose walls consist mostly of hydrophilic residues, transporting through the selectivity filter containing Arg-227, His-212, Phe-92, and Ala-221, and the two conserved Asn-Pro-Ala (NPA) motifs containing asparagines 224 and 112. In summary, not all fungal aquaporins share the same gating mechanism by a tyrosine residue in the N-terminus, and the structural analysis in the present study should aid our understanding of aquaporin structure and its functional implications.     

Water transport in AQP0 aquaporin: molecular dynamics studies

A double lipid bilayer structure containing opposing tetramers of AQP0 aquaporin, in contact through extracellular face loop regions, was recently modeled using an intermediate-resolution map obtained by electron crystallographic methods. The pores of these water channels were found to be critically narrow in three regions and subsequently interpreted to be those of a closed state of the channel. The subsequent determination of a high-resolution AQP0 tetramer structure by X-ray crystallographic methods yielded a pore model featuring two of the three constrictions as noted in the EM work and water molecules within the channel pore. The extracellular-side constriction region of this AQP0 structure was significantly larger than that of the EM-based model and similar to that of the highly water permeable AQP1. The X-ray-based study of AQP0 however could not ascertain if the water molecules found in the pore were the result of water entering from one or both ends of the channel, nor whether water could freely pass through all constriction points. Additionally, this X-ray-based structure could not provide an answer to the question of whether the double lipid bilayer configuration of AQP0 could functionally maintain a water impermeable state of the channel. To address these questions we conducted molecular dynamics simulations to compare the time-dependent behavior of the AQP0 and AQP1 channels within lipid bilayers. The simulations demonstrate that AQP0, in single or double lipid bilayers, is not closed to water transport and that thermal motions of critical side-chains are sufficient to facilitate the movement of water past any of its constriction regions. These motional requirements do however lead to significant free energy barriers and help explain physiological observations that found water permeability in AQP0 to be substantially lower than in the AQP1 pore.     

Exploring gas permeability of cellular membranes and membrane channels with molecular dynamics

Aquaporins are a family of membrane proteins specialized in rapid water conduction across biological membranes. Whether these channels also conduct gas molecules and the physiological significance of this potential function have not been well understood. Here we report 140 ns of molecular dynamics simulations of membrane-embedded AQP1 and of a pure POPE bilayer addressing these questions. The permeability of AQP1 to two types of gas molecules, O2 and CO2, was investigated using two complementary methods, namely, explicit gas diffusion simulation and implicit ligand sampling. The simulations show that the central (tetrameric) pore of AQP1 can be readily used by either gas molecule to permeate the channel. The two approaches produced similar free energy profiles associated with gas permeation through the central pore: a -0.4 to -1.7 kcal/mol energy well in the middle, and a 3.6-4.6 kcal/mol energy barrier in the periplasmic vestibule. The barrier appears to be mainly due to a dense cluster of water molecules anchored in the periplasmic mouth of the central pore by four aspartate residues. Water pores show a very low permeability to O2, but may contribute to the overall permeation of CO2 due to its more hydrophilic nature. Although the central pore of AQP1 is found to be gas permeable, the pure POPE bilayer provides a much larger cross-sectional area, thus exhibiting a much lower free energy barrier for CO2 and O2 permeation. As such, gas conduction through AQP1 may only be physiologically relevant either in membranes of low gas permeability, or in cells where a major fraction of the cellular membrane is occupied by AQPs.     

Multiple pore lining residues modulate water permeability of GlpF

The water permeability of aquaporins (AQPs) varies by more than an order of magnitude even though the pore structure, geometry, as well as the channel lining residues are highly conserved. However, channel gating by pH, divalent ions or phosphorylation was only shown for a minority of AQPs. Structural and in silico indications of water flux modulation by flexible side chains of channel lining residues have not been experimentally confirmed yet. Hence, the aquaporin “open state” is still considered to be a continuously open pore with water molecules permeating in a single‐file fashion. Using protein mutations outside the selectivity filter in the aqua(glycerol)facilitator GlpF of Escherichia coli we, to the best of our knowledge, for the first time, modulate the position of the highly conserved Arg in the selectivity filter. This in turn enhances or reduces the unitary water permeability of GlpF as shown in silico by molecular dynamics (MD) simulations and in vitro with purified and reconstituted GlpF. This finding suggests that AQP water permeability can indeed be regulated by lipid bilayer asymmetry and the transmembrane potential. Strikingly, our long‐term MD simulations reveal that not only the conserved Arg in the selectivity filter, but the position and dynamics of multiple other pore lining residues modulate water passage through GlpF. This finding is expected to trigger a wealth of future investigations on permeability and regulation of AQPs among others with the aim to tune water permeability for biotechnological applications.     

Comparative simulations of aquaporin family: AQP1, AQPZ, AQP0 and GlpF

Molecular dynamics simulations were performed for four members of the aquaporin family (AQP1, AQPZ, AQP0, and GlpF) in the explicit membrane environment. The single-channel water permeability, pf, was evaluated to be GlpF approximately AQPZ > AQP1 > AQP0, while their relative pore sizes were GlpF > AQP1 > AQPZ > AQP0. This relation between pf and pore size indicates that water permeability was determined not only by the channel radius, but also another competing factor. Analysis of water dynamics revealed that this factor was the single-file nature of water transport.     

Selectivity and conductance among the glycerol and water conducting aquaporin family of channels

The atomic structures of a transmembrane water plus glycerol conducting channel (GlpF), and now of aquaporin Z (AqpZ) from the same species, Escherichia coli, bring the total to three atomic resolution structures in the aquaporin (AQP) family. Members of the AQP family each assemble as tetramers of four channels. Common helical axes support a wider channel in the glycerol plus water channel paradigm, GlpF. Water molecules form a single hydrogen bonded file throughout the 28 A long channel in AqpZ. The basis for absolute exclusion of proton or hydronium ion conductance through the line of water is explored using simulations.     

Simulation of the N-terminus of HIV-1 glycoprotein 41000 fusion peptide in micelles.

In this paper, the N-terminus of glycoprotein-41, the HIV-1 fusion peptide, was studied by molecular dynamics simulations in an explicit sodium dodecyl sulfate micelle. The simulation provides a detailed picture of the equilibrium structure and peptide stability as it interacts with the micelle. The equilibrium location of the peptide shows the peptide at the surface of the micelle with hydrophobic residues interacting with the micelle's core. At equilibrium, the peptide adopts an alpha-helical structure from residues 5-16 and a type-1 beta-turn from 17-20 with the other residues exhibiting more flexible conformations. The primary hydrophobic interactions with the micelle are from the leucine and phenylalanine residues (Leu-7, Phe-8, Leu-9, Phe-11, Leu-12) while the alanine and glycine residues (Ala-1, Gly-3, Gly-5, Ala-6, Gly-10, Gly-13, Ala-14, Ala-15, Gly-16, Gly-10, Ala-21) interact favorably with water molecules. The results suggest that Phe-8, part of the highly conserved FLG motif of the fusion peptide, plays a key role in the interaction of the peptide with membranes. Our simulations corroborate experimental investigations of the fusion peptide in SDS micelles, providing a high-resolution picture that explains the experimental findings.     

Dynamic and energetic mechanisms for the distinct permeation rate in AQP1 and AQP0

Despite sharing overall sequence and structural similarities, water channel aquaporin 0 (AQP0) transports water more slowly than other aquaporins. Using molecular dynamics simulations of AQP0 and AQP1, we find that there is a sudden decrease in the distribution profile of water density along the pore of AQP0 in the region of residue Tyr23, which significantly disrupts the single file water chain by forming hydrogen bond with permeating water molecules. Comparisons of free-energy and interaction-energy profiles for water conduction between AQP0 and AQP1 indicate that this interruption of the water chain causes a huge energy barrier opposing water translocation through AQP0. We further show that a mutation of Tyr23 to phenylalanine leads to a 2- to 4-fold enhancement in water permeability of AQP0, from (0.5 ± 0.2) × 10^− 14 cm³s^− 1 to (1.9 ± 0.6) × 10^− 14 cm³s^− 1. Therefore, Tyr23 is a dominate factor leading to the low water permeability in AQP0.     

Influence of the Cytoplasmic Domains of Aquaporin-4 on Water Conduction and Array Formation

Phosphorylation of Ser180 in cytoplasmic loop D has been shown to reduce the water permeability of aquaporin (AQP) 4, the predominant water channel in the brain. However, when the structure of the S180D mutant (AQP4M23S180D), which was generated to mimic phosphorylated Ser180, was determined to 2.8 Å resolution using electron diffraction patterns, it showed no significant differences from the structure of the wild-type channel. High-resolution density maps usually do not resolve protein regions that are only partially ordered, but these can sometimes be seen in lower-resolution density maps calculated from electron micrographs. We therefore used images of two-dimensional crystals and determined the structure of AQP4M23S180D at 10 Å resolution. The features of the 10-Å density map are consistent with those of the previously determined atomic model; in particular, there were no indications of any obstruction near the cytoplasmic pore entrance. In addition, water conductance measurements, both in vitro and in vivo, show the same water permeability for wild-type and mutant AQP4M23, suggesting that the S180D mutation neither reduces water conduction through a conformational change nor reduces water conduction by interacting with a protein that would obstruct the cytoplasmic channel entrance. Finally, the 10-Å map shows a cytoplasmic density in between four adjacent tetramers that most likely represents the association of four N termini. This finding supports the critical role of the N terminus of AQP4 in the stabilization of orthogonal arrays, as well as their interference through lipid modification of cysteine residues in the longer N-terminal isoform.     

Water transport in aquaporins: osmotic permeability matrix analysis of molecular dynamics simulations

Single-channel osmotic water permeability (p(f)) is a key quantity for investigating the transport capability of the water channel protein, aquaporin. However, the direct connection between the single scalar quantity p(f) and the channel structure remains unclear. In this study, based on molecular dynamics simulations, we propose a p(f)-matrix method, in which p(f) is decomposed into contributions from each local region of the channel. Diagonal elements of the p(f) matrix are equivalent to the local permeability at each region of the channel, and off-diagonal elements represent correlated motions of water molecules in different regions. Averaging both diagonal and off-diagonal elements of the p(f) matrix recovers p(f) for the entire channel; this implies that correlated motions between distantly-separated water molecules, as well as adjacent water molecules, influence the osmotic permeability. The p(f) matrices from molecular dynamics simulations of five aquaporins (AQP0, AQP1, AQP4, AqpZ, and GlpF) indicated that the reduction in the water correlation across the Asn-Pro-Ala region, and the small local permeability around the ar/R region, characterize the transport efficiency of water. These structural determinants in water permeation were confirmed in molecular dynamics simulations of three mutants of AqpZ, which mimic AQP1.     

A new gating site in human aquaporin-4: Insights from molecular dynamics simulations

Aquaporin-4 (AQP4) is the predominant water channel in different organs and tissues. An alteration of its physiological functioning is responsible for several disorders of water regulation and, thus, is considered an attractive target with a promising therapeutic and diagnostic potential. Molecular dynamics (MD) simulations performed on the AQP4 tetramer embedded in a bilayer of lipid molecules allowed us to analyze the role of spontaneous fluctuations occurring inside the pore. Following the approach by Hashido et al. [Hashido M, Kidera A, Ikeguchi M (2007) Biophys J 93: 373–385], our analysis on 200 ns trajectory discloses three domains inside the pore as key elements for water permeation. Herein, we describe the gating mechanism associated with the well-known selectivity filter on the extracellular side of the pore and the crucial regulation ensured by the NPA motifs (asparagine, proline, alanine). Notably, on the cytoplasmic side, we find a putative gate formed by two residues, namely, a cysteine belonging to the loop D (C178) and a histidine from loop B (H95). We observed that the spontaneous reorientation of the imidazole ring of H95 acts as a molecular switch enabling H-bond interaction with C178. The occurrence of such local interaction seems to be responsible for the narrowing of the pore and thus of a remarkable decrease in water flux rate. Our results are in agreement with recent experimental observations and may represent a promising starting point to pave the way for the discovery of chemical modulators of AQP4 water permeability.     

1,3‐propanediol binds deep inside the channel to inhibit water permeation through aquaporins

Aquaporins and aquaglyceroporins (AQPs) are membrane channel proteins responsible for transport of water and for transport of glycerol in addition to water across the cell membrane, respectively. They are expressed throughout the human body and also in other forms of life. Inhibitors of human AQPs have been sought for therapeutic treatment for various medical conditions including hypertension, refractory edema, neurotoxic brain edema, and so forth. Conducting all‐atom molecular dynamics simulations, we computed the binding affinity of acetazolamide to human AQP4 that agrees closely with in vitro experiments. Using this validated computational method, we found that 1,3‐propanediol (PDO) binds deep inside the AQP4 channel to inhibit that particular aquaporin efficaciously. Furthermore, we used the same method to compute the affinities of PDO binding to four other AQPs and one aquaglyceroporin whose atomic coordinates are available from the protein data bank (PDB). For bovine AQP1, human AQP2, AQP4, AQP5, and Plasmodium falciparum PfAQP whose structures were resolved with high resolution, we obtained definitive predictions on the PDO dissociation constant. For human AQP1 whose PDB coordinates are less accurate, we estimated the dissociation constant with a rather large error bar. Taking into account the fact that PDO is generally recognized as safe by the US FDA, we predict that PDO can be an effective diuretic which directly modulates water flow through the protein channels. It should be free from the serious side effects associated with other diuretics that change the hydro‐homeostasis indirectly by altering the osmotic gradients.     

The influences of protonation state of histidine on aromatic/arginine region of aquaporin-1 protein

Aquaporin-1 (AQP1) is widely distributed in the epithelial tissue for water absorption and secretion. The histidine (His182) residue in the aromatic/arginine (ar/R) constriction region plays an important role in transporting water through the membrane. In this study, we have performed a total of 46 ns equilibrium molecular dynamics (MD) simulations, and obtained the influence of His182 in two protonation states (Hsd is the proton at N_δ and Hse is the proton at N_?) on the ar/R region. Water permeation rate shows that it is easier for water molecules to permeate the ar/R region of the AQP1 with residue in the Hsd state than in the Hse state. The minimum radii of the pore in the ar/R region were calculated during the last 10 ns MD simulation. We have analysed the correlation among the state of the pore (open or close), the minimum radius of the ar/R region and the dihedral angles < C_β-C_γ-C_δ-N_? of Arg197. The results show that the minimum radius can be used to mark the state of the pore.     

Ion exclusion mechanism in aquaporin at an atomistic level

Although the mechanism of proton exclusion in aquaporin is investigated by many researchers, the detailed molecular mechanism for ion exclusion in aquaporin is still not completely understood. In the present work, a detailed mechanism for ion exclusion in aquaporin-1 (AQP1) at an atomistic level is investigated by calculating the free energy for transport of ions in AQP1 using an atomistic molecular dynamics simulation. For this purpose, sodium and chloride ions are chosen as representatives for nonprotonic ions. The simulation shows that the free energy barrier showing its maximum is located at the NPA region for sodium ion while it is located at both the front and the rear for chloride ion and that the barrier height is 18 and 9 kcal/mol, respectively, indicating that the ions are not able to pass through aquaporin. Analysis of the pair interaction energy between the permeating ion and its environment reveals that sodium ion is excluded by the positive charge generated by two alpha-helical macro-dipoles, while chloride ion is expelled by carbonyl oxygen atoms protruding from pore-making residues before it reaches the NPA motif. It is also found that the number of water molecules hydrating the ions is reduced as the ions enter the pore, implying that the energetic cost for detaching water molecules from a permeating ion also contributes to the free energy barriers of ion transport in AQP1.     

Ammonium recruitment and ammonia transport by E. coli ammonia channel AmtB

To investigate substrate recruitment and transport across the Escherichia coli Ammonia transporter B (AmtB) protein, we performed molecular dynamics simulations of the AmtB trimer. We have identified residues important in recruitment of ammonium and intraluminal binding sites selective of ammonium, which provide a means of cation selectivity. Our results indicate that A162 guides translocation of an extraluminal ammonium into the pore lumen. We propose a mechanism for transporting the intraluminally recruited proton back to periplasm. Our mechanism conforms to net transport of ammonia and can explain why ammonia conduction is lost upon mutation of the conserved residue D160. We unify previous suggestions of D160 having either a structural or an ammonium binding function. Finally, our simulations show that the channel lumen is hydrated from the cytoplasmic side via the formation of single file water, while the F107/F215 stack at the inner-most part of the periplasmic vestibule constitutes a hydrophobic filter preventing AmtB from conducting water.     

Water channel structures analysed by electron crystallography

Background

The mechanisms underlying water transport through aquaporin (AQP) have been debated for two decades. The water permeation phenomenon of AQP seems inexplicable because the Grotthuss mechanism does not allow for simultaneous fast water permeability and inhibition of proton transfer through the hydrogen bonds of water molecules.

Scope of review

The AQP1 structure determined by electron crystallography provided the first insights into the proton exclusion mechanism despite fast water permeation. Although several studies have provided clues about the mechanism based on the AQP structure, each proposed mechanism remains incomplete. The present review is focused on AQP function and structure solved by electron crystallography in an attempt to fill the gaps between the findings in the absence and presence of lipids.

Major conclusions

Many AQP structures can be superimposed regardless of the determination method. The AQP fold is preserved even under conditions lacking lipids, but the water arrangement in the channel pore differs. The differences might be explained by dipole moments formed by the two short helices in the lipid bilayer. In addition, structure analyses of double-layered two-dimensional crystals of AQP suggest an array formation and cell adhesive function.

General significance

Electron crystallography findings not only have contributed to resolve some of the water permeation mechanisms, but have also elucidated the multiple functions of AQPs in the membrane. The roles of AQPs in the brain remain obscure, but their multiple activities might be important in the regulation of brain and other biological functions. This article is part of a Special Issue entitled Aquaporins.     

Nitric oxide conduction by the brain aquaporin AQP4

Involvement of aquaporins in gas conduction across the membrane and the physiological significance of this process have attracted marked attention from both experimental and theoretical studies. Previous work demonstrated that AQP1 is permeable to both CO₂ and O₂. Here we employ various simulation techniques to examine the permeability of the brain aquaporin AQP4 to NO and O₂ and to describe energetics and pathways associated with these phenomena. The energy barrier to NO and O₂ permeation through AQP4 central pore is found to be only ～3 kcal mol^?1. The results suggest that the central pore of AQP4, similar to that of AQP1, can indeed conduct gas molecules. Interestingly, despite a longer and narrower central pore, AQP4 appears to provide an energetically more favorable permeation pathway for gas molecules than AQP1, mainly due to the different orientation of its charged residues near the pore entrance. Although the low barrier against gas permeation through AQP4 indicates that it can participate in gas conduction across the cellular membrane, physiological relevance of the phenomenon remains to be established experimentally, particularly since pure lipid bilayers appear to present a more favorable pathway for gas conduction across the membrane. With an energy well of ?1.8 kcal mol^?1, the central pore of AQP4 may also act as a reservoir for NO molecules to accumulate in the membrane. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Energetics of Transport through the Nuclear Pore Complex

Molecular transport across the nuclear envelope in eukaryotic cells is solely controlled by the nuclear pore complex (NPC). The NPC provides two types of nucleocytoplasmic transport: passive diffusion of small molecules and active chaperon-mediated translocation of large molecules. It has been shown that the interaction between intrinsically disordered proteins that line the central channel of the NPC and the transporting cargoes is the determining factor, but the exact mechanism of transport is yet unknown. Here, we use coarse-grained molecular dynamics simulations to quantify the energy barrier that has to be overcome for molecules to pass through the NPC. We focus on two aspects of transport. First, the passive transport of model cargo molecules with different sizes is studied and the size selectivity feature of the NPC is investigated. Our results show that the transport probability of cargoes is significantly reduced when they are larger than ∼5 nm in diameter. Secondly, we show that incorporating hydrophobic binding spots on the surface of the cargo effectively decreases the energy barrier of the pore. Finally, a simple transport model is proposed which characterizes the energy barrier of the NPC as a function of diameter and hydrophobicity of the transporting particles.     

短杆菌肽A-DMPC通道内离子输运的分子动力学模拟

用最近提出的构建膜体系初始构象的有效方法 ,构建了在DMPC脂膜环境下短杆菌肽A通道模型 (GA -DMPC)。通过对Na 、Ca2 、Cl-三种不同离子在GA -DMPC通道内不同位置的分子动力学模拟 ,研究离子在通道内输运过程中与通道及通道内水分子的相互作用 ,从分子动力学的角度阐明离子在通道内的输运机制。主要计算结果表明 :(1)离子在通道内的输运使GA的构象发生变化 ,GA的柔性是离子在通道内通透的重要因素 ;(2)Cl- 离子可扩大通道半径 ,Na 离子和Ca2 离子则减小通道半径。Cl-离子不能在GA通道内通透 ;(3)离子的出现使通道内水分子的偶极方向发生变化。上述结果均与实验相符。     

Formation of aquaporin-4 arrays is inhibited by palmitoylation of N-terminal cysteine residues

Tetramers of the mammalian water channel aquaporin-4 (AQP4) assemble into square arrays and mediate bidirectional water transport across the blood-brain interface. The aqp4 gene expresses two splicing isoforms. Only the shorter AQP4M23 isoform assembles into square arrays, while the longer AQP4M1 isoform interferes with array formation, presumably due to the additional 22 N-terminal residues. To understand why the N-terminus of AQP4M1 interferes with array formation, we constructed a series of N-terminal deletion mutants and examined their ability to form square arrays in Chinese hamster ovary (CHO) cells using SDS-digested freeze fracture replica labeling. Mutants with deletions of less than seventeen N-terminal residues did not form square arrays and showed dispersed immunogold labels against AQP4 molecules, whereas more deletions led to the formation of square arrays labeled with immunogolds. Furthermore, mutagenic substitution of the two cysteine residues at the position 13 and 17 in the N-terminus of AQP4M1 also resulted in the square array formation. Biochemical analysis and metabolic labeling of transfected CHO cells revealed that the two N-terminal cysteines of AQP4M1 are palmitoylated. These results suggest that palmitoylation of the N-terminal cysteines is the reason for the inability of AQP4M1 to form square arrays.