Microsecond simulation of human aquaporin 2 reveals structural determinants of water permeability and selectivity

Many kinases for controlling the water channel aquaporin-2

Aquaporin-2 (AQP2) is a member of the aquaporin water channel family. In the kidney, AQP2 is expressed in collecting duct principal cells where it facilitates water reabsorption in response to antidiuretic hormone (arginine vasopressin, AVP). AVP induces the redistribution of AQP2 from intracellular vesicles and its incorporation into the plasma membrane. The plasma membrane insertion of AQP2 represents the crucial step in AVP-mediated water reabsorption. Dysregulation of the system preventing the AQP2 plasma membrane insertion causes diabetes insipidus (DI), a disease characterised by an impaired urine concentrating ability and polydipsia. There is no satisfactory treatment of DI available. This review discusses kinases that control the localisation of AQP2 and points out potential kinase-directed targets for the treatment of DI.

Aquaporins Display a Diversity in their Substrates

Aquaporins constitute a family of transmembrane proteins that function to transport water and other small solutes across the cell membrane. Aquaporins family members are found in diverse life forms. Aquaporins share the common structural fold consisting of six transmembrane alpha helices with a central water-transporting channel. Four such monomers assemble together to form tetramers as their biological unit. Initially, aquaporins were discovered as water-transporting channels, but several studies supported their involvement in mediating the facilitated diffusion of different solutes. The so-called water channel is able to transport a variety of substrates ranging from a neutral molecule to a charged molecule or a small molecule to a bulky molecule or even a gas molecule. This article gives an overview of a diverse range of substrates conducted by aquaporin family members. Prime focus is on human aquaporins where aquaporins show a wide tissue distribution and substrate specificity leading to various physiological functions. This review also highlights the structural mechanisms leading to the transport of water and glycerol. More research is needed to understand how one common fold enables the aquaporins to transport an array of solutes.

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.

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.

Molecular characterization and structural dynamics of Aquaporin1 from walking catfish in lipid bilayers

Aquaporin's (AQPs) are the major superfamily of small integral membrane proteins that facilitates transportation of water, urea, ammonia, glycerol and ions across biological cell membranes. Despite of recent advancements made in understanding the biology of Aquaporin's, only few isoforms of aquaporin 1 (AQP1) some of the teleost fish species have been characterized at molecular scale. In this study, we made an attempt to elucidate the molecular mechanism of water transportation in AQP1 from walking catfish (Clarias batrachus), a model species capable of breathing in air and inhabits in challenging environments. Using state-of-the-art computational modelling and all-atoms molecular dynamics simulation, we explored the structural dynamics of full-length aquaporin 1 from walking catfish (CbAQP1) in lipid mimetic bilayers. Unlike AQP1 of human and bovine, structural ensembles of CbAQP1 from MD revealed discrete positioning of pore lining residues at the intracellular end. Snapshots from MD simulation displayed differential dynamics of aromatic/arginine (ar/R) filter and extracellular loop C bridging transmembrane (TM) helix H3 and H4. Distinct conformation of large extracellular loops, loop bridging TM2 domain and HB helix along with positioning of selectivity filter lining residues controls the permeability of water across the bilayer. Moreover, the identified unique and conserved lipid binding sites with 100% lipid occupancy signifies lipid mediated structural dynamics of CbAQP1. All-together, this is the first ever report on structural-dynamics of aquaporin 1 in walking catfish which will be useful to understand the molecular basis of transportation of water and other small molecules under varying degree of hyperosmotic environment.

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.

Urea-aromatic interactions in biology

Noncovalent interactions are key determinants in both chemical and biological processes. Among such processes, the hydrophobic interactions play an eminent role in folding of proteins, nucleic acids, formation of membranes, protein-ligand recognition, etc.. Though this interaction is mediated through the aqueous solvent, the stability of the above biomolecules can be highly sensitive to any small external perturbations, such as temperature, pressure, pH, or even cosolvent additives, like, urea-a highly soluble small organic molecule utilized by various living organisms to regulate osmotic pressure. A plethora of detailed studies exist covering both experimental and theoretical regimes, to understand how urea modulates the stability of biological macromolecules. While experimentalists have been primarily focusing on the thermodynamic and kinetic aspects, theoretical modeling predominantly involves mechanistic information at the molecular level, calculating atomistic details applying the force field approach to the high level electronic details using the quantum mechanical methods. The review focuses mainly on examples with biological relevance, such as (1) urea-assisted protein unfolding, (2) urea-assisted RNA unfolding, (3) urea lesion interaction within damaged DNA, (4) urea conduction through membrane proteins, and (5) protein-ligand interactions those explicitly address the vitality of hydrophobic interactions involving exclusively the urea-aromatic moiety.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

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.