Molecular dynamics insights into human aquaporin 2 water channel

Functional regulation of aquaporin dynamics by lipid bilayer composition

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics are examined. We demonstrate that SoPIP2;1's structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

Functional Regulation of Aquaporin Dynamics by Lipid Bilayer Composition

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics were examined. We demonstrate that SoPIP2;1s structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

Non-equilibrium molecular dynamics study of human aquaporin-2 in the static external electric fields

In this paper, non-equilibrium MD simulations (NEMD) of human aquaporin-2 (AQP2) in the presence of an external static electric field have been performed along + z and - z directions of the pore axis. The impacts of the electric field direction on the gating mechanism corresponding to the selectivity filter (SF) region of AQP2 have been studied. Besides, the effects of applied external electric field on the PMF profile of water molecules translocation, water permeability, and molecules dipole orientation are investigated. Our results showed that when the external electric field is implemented along the + z direction of the channels, the selectivity filter region is kept in the wide conformation for the majority of the time. Therefore, a remarkable increase in the overall water permeability can be seen compared to the case without any external electric field. This is in contrast to the effects of - z-directed electric field on the conformations of the selectivity filter, which induces mostly narrow conformations in this constriction region. A substantial higher energy barrier emerged in the middle of the AQP2's pores under the effect of -z-directed electric field in comparison with the zero and + z-directed electric field strengths, which is mainly ascribed to the deviation from bipolar dipole orientation within the AQP2's pores.Communicated by Ramaswamy H. Sarma.

Deciphering the role of aquaporin 1 in the adaptation of the stinging catfish Heteropneustes fossilis to environmental hypertonicity by molecular dynamics simulation studies

A thorough investigation of the water permeability of H. fossilis aquaporin 1 (hfAQP1) in a hypertonic environment can provide a useful insight into the understanding of the underlying molecular mechanism of its high tolerance to salinity. Here, we constructed a 3 D homology model of hfAQP1 by taking Bos taurus AQP1, AQP0, and human AQP2 as templates using I-TASSER. The model obtained has similar structural organizations with mammalian AQP1s in all aspects. We investigated the water permeability of the modeled hfAQP1 in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane under neutral and 100 mM hypersalinity by subjecting each system to a 100 ns molecular dynamics simulation. Our results show that hypersalinity hinders water permeation across the membrane through the hfAQP1 channel. A change in the intermolecular distance between key residues of the ar/R selectivity filter along with charge redistribution resulted in the accommodation of only 2-6 water molecules inside the channel at once under hypersaline conditions. We investigated the mRNA expression pattern of hfaqp1 in osmoregulatory organs of H. fossilis in response to 100 mM hypertonicity by using qPCR analysis. The transcript was downregulated in kidney and GI tract, but upregulated in the Gills. Thus, the catfish survive in a hypertonic environment by reducing the transport of water in its cellular systems and downregulating the expression of the hfaqp1 gene. The results observed in our study can shed more light on the functionality of AQP1 in catfishes under salinity stress and aid in future researches on solving more gating mechanisms involved in its regulation.Communicated by Ramaswamy H. Sarma.

Molecular dynamics study of water transport through AQP5-R188C mutant causing palmoplantar keratoderma (PPK) using the gating mechanism concept

It is widely known that any disruption to the water regulation in aquaporins (AQPs) leads to numerous important diseases. However, studies of dynamics and energetics of disease-causing mutations in the aquaporins on the molecular level are still limited. In the present work, the effects of a skin disease-causing mutant, R188C, on the structure of AQP5 and water transport mechanism within this mutated aquaporin are investigated using the concept of gating mechanism. Our results have revealed that the R188C mutation causes a remarkable increase in the pore radius inside the selectivity filter (SF) region facilitating the passage of water molecules. This observation is supported by plotting the free energy profiles of water molecules transport and calculating permeability values through AQP5-R188C, such that the energy barrier in the SF region of the pores was substantially reduced by this mutation, and therefore, the translocation of water molecules was improved. The total averaged osmotic permeability for R188C has been computed as about 11-fold of the wild-type permeability. However, a comparison between the osmotic permeability values related to the open conformation of CE revealed that this coefficient for AQP5-R188C is about 6.5 times larger than that of wt-AQP5, which can be a more accurate value according to the gating mechanism associated with the constriction region of the aquaporin.

Dynamics and energetics of water transport through aquaporin mutants causing nephrogenic diabetes insipidus (NDI): A molecular dynamics study

Human aquaporin-2 (AQP2) is the principal water channel in the human kidney. Any alteration of its physiological functioning may lead to the water imbalance and consequently diseases in humans, especially nephrogenic diabetes insipidus (NDI). Although many of the mutations associated with NDI are experimentally discovered and examined, a molecular level characterization of the structure and transport mechanism is still missing. In this paper, the gating effects of selectivity filter (SF) as wide/narrow states on the mechanism and dynamics of water permeation within the wild-type AQP2 and two NDI causing mutants as AQP2-V168M and AQP2-G64R are studied for the first time. The analysis of the 200 ns trajectory shows that the SF region in AQP2 is not only a selectivity filter, as previously reported but also it performs as a gating site depending on the side-chain conformation of His172. The assignment of the wide/narrow states of SF is supported by computing the free energy and permeability through the AQP2. Moreover, by exploring the effects of V168M and G64R mutants on the AQP2 structure during 200 ns trajectories, remarkable increases of energy barriers are observed in the middle and cytoplasmic side of the pore, respectively. Interestingly, it is found that due to the variable conformations of the SF region as wide/narrow, the effect of the NDI causing mutants on the average water permeability can be revealed with notably better accuracy by finding the wide states in the wild-type and mutated types of AQP2 and comparing the osmotic permeabilities for this state.Communicated by Ramaswamy H. Sarma.

Investigating the Mechanisms of AquaporinZ Reconstitution through Polymeric Vesicle Composition for a Biomimetic Membrane

Aquaporin-Z (AqpZ) are water channel proteins with excellent water permeability and solute rejection properties. AqpZ can be reconstituted into vesicles utilizing cell-like bilayer membranes assembled from amphiphilic block copolymers, for the preparation of high-performance biomimetic membranes. However, only a few copolymers have been found suitable to act as the membrane matrix for protein reconstitution. Hence, this work analyzes the mechanism of protein reconstitution based on a composition-reconstitution relationship. The vesicle formation and AqpZ reconstitution processes in various amphiphilic block copolymers were investigated in terms of size, morphology, stability, polymeric bilayer membrane rigidity, and thermal behavior. Overall, this study contributes to the understanding of the composition-reconstitution relationship of biomimetic membranes based on AqpZ-reconstituted polymeric vesicles.

Advances in aquaporin-2 trafficking mechanisms and their implications for treatment of water balance disorders

In mammals, conservation of body water is critical for survival and is dependent on the kidneys' ability to minimize water loss in the urine during periods of water deprivation. The collecting duct water channel aquaporin-2 (AQP2) plays an essential role in this homeostatic response by facilitating water reabsorption along osmotic gradients. The ability to increase the levels of AQP2 in the apical plasma membrane following an increase in plasma osmolality is a rate-limiting step in water reabsorption, a process that is tightly regulated by the antidiuretic hormone arginine vasopressin (AVP). In this review, the focus is on the role of the carboxyl-terminus of AQP2 as a key regulatory point for AQP2 trafficking. We provide an overview of AQP2 structure, disease-causing mutations in the AQP2 carboxyl-terminus, the role of posttranslational modifications such as phosphorylation and ubiquitylation in the tail domain, and their implications for balanced trafficking of AQP2. Finally, we discuss how various modifications of the AQP2 tail facilitate selective protein-protein interactions that modulate the AQP2 trafficking mechanism.

A novel molecular dynamics study of CO2 permeation through aquaporin-5

Aquaporins (AQPs) are protein channels which facilitate rapid water permeation across cell membrane. The AQPs are very vital for biological organs, as their malfunction causes severe diseases in human body. A particular family of AQPs, that is AQP5, has a significant role in lung fluid transport due to submucosal glands structure. However, it has not been yet well understood whether these protein channels can conduct gas molecules. Here, Molecular Dynamics (MD) simulations are used to investigate the CO₂ permeability and diffusion in AQP5 during a 40-nanosecond period. For the first time, equilibrium and Steered MD (SMD) are used to simulate self and force-induced diffusion of CO₂ molecules across AQP5 and POPE lipid bilayer. According to PMFs profile associated to CO₂ permeation, the hydrophobic central pore provides a more suitable pathway for gas molecules compared to other AQP5 channels. Although CO₂ molecules can also permeate across AQP5 water channels, the rate of CO₂ permeation through four channels of the AQP5 monomers is much lower than the central pore. The rate of CO₂ permeation through four AQP5 water channels is even lower than CO₂ diffusion through POPE lipid membrane. The results reported in this investigation demonstrate that MD simulations of human AQP5 provide valuable insights into the gas permeation mechanism for both the equilibrium self-diffusion, and quasi-equilibrium condition.

Forced diffusion of water molecules through aquaporin-5 biomembrane; a molecular dynamics study

Aquaporins (AQPs) are protein channels located across the cell membrane which conduct the water permeation through the cell membrane. Different types of AQPs exist in human organs and play vital roles, as the malfunction of such protein membranes can lead to life-threatening conditions. A specific type of AQP, identified as AQP5, is particularly essential to the generation of saliva, tears and pulmonary secretions. We have adopted Molecular Dynamics (MD) simulation to analyze the water permeation and diffusion in AQP5 structure in a 0.5 microsecond simulation time window. The MD numerical simulation shows the water permeability of the human AQP5 is in the nominal range for other members of human aquaporins family. In addition, we have considered the effect of the osmotic water diffusion and the diffusion occurred by pressure gradient on the protein membrane. The water permeability grows monotonically as the applied pressure on the solvent increases. Furthermore, the forced diffusion increases the minimum radius of Selectivity Filter (SF) region of region AQP5 up to 20% and consequently the permeability coefficients enhance enormously compared to osmotic self-diffusion in AQP5 tetramer. Finally, it is revealed that the MD simulation of human AQP5 provides useful insights into the mechanisms of water regulation through alveolar cells under the different physical conditions; osmotic self-diffusion and forced diffusion condition.

Molecular Dynamics Simulation Study of Polyamide Membrane Structures and RO/FO Water Permeation Properties

Polyamide (PA) membranes possess properties that allow for selective water permeation and salt rejection, and these are widely used for reverse osmotic (RO) desalination of sea water to produce drinking water. In order to design high-performance RO membranes with high levels of water permeability and salt rejection, an understanding of microscopic PA membrane structures is indispensable, and this includes water transport and ion rejection mechanisms on a molecular scale. In this study, two types of virtual PA membranes with different structures and densities were constructed on a computer, and water molecular transport properties through PA membranes were examined on a molecular level via direct reverse/forward osmosis (RO/FO) filtration molecular dynamics (MD) simulations. A quasi-non-equilibrium MD simulation technique that uses applied (RO mode) or osmotic (FO mode) pressure differences of several MPa was conducted to estimate water permeability through PA membranes. A simple NVT (Number, Volume, and Temperature constant ensemble)-RO MD simulation method was presented and verified. The simulations of RO and FO water permeability for a dense PA membrane model without a support layer agreed with the experimental value in the RO mode. This PA membrane completely rejected Na⁺ and Cl⁻ ions during a simulation time of several nano-seconds. The naturally dense PA structure showed excellent ion rejection. The effect that the void size of PA structure exerted on water permeability was also examined.