The mechanism of glycerol conduction in aquaglyceroporins

The Aquaporin gene family of the yellow fever mosquito, Aedes aegypti

Background

The mosquito, Aedes aegypti, is the principal vector of the Dengue and yellow fever viruses. During feeding, an adult female can take up more than its own body weight in vertebrate blood. After a blood meal females excrete large amounts of urine through their excretion system, the Malpighian tubules (MT). Diuresis starts within seconds after the mosquito starts feeding. Aquaporins (AQPs) are a family of membrane transporters that regulate the flow of water, glycerol and other small molecules across cellular membranes in both prokaryotic and eukaryotic cells. Our aim was to identify aquaporins that function as water channels, mediating transcellular water transport in MTs of adult female Ae. aegypti.

Methodology/Principal Findings

Using a bioinformatics approach we screened genome databases and identified six putative AQPs in the genome of Ae. aegypti. Phylogenetic analysis showed that five of the six Ae. aegypti AQPs have high similarity to classical water-transporting AQPs of vertebrates. Using microarray, reverse transcription and real time PCR analysis we found that all six AQPs are expressed in distinct patterns in mosquito tissues/body parts. AaAQP1, 4, and 5 are strongly expressed in the adult female MT. RNAi-mediated knockdown of the MT-expressed mosquito AQPs resulted in significantly reduced diuresis.

Conclusions/Significance

Our results support the notion that AQP1, 4, and 5 function as water transporters in the MTs of adult female Ae. aegypti mosquitoes. Our results demonstrate the importance of these AQPs for mosquito diuresis after blood ingestion and highlight their potential as targets for the development of novel vector control strategies.

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.

Free-energy landscape of glycerol permeation through aquaglyceroporin GlpF determined from steered molecular dynamics simulations

The free-energy landscape of glycerol permeation through the aquaglyceroporin GlpF has been estimated in the literature by the nonequilibrium method of steered molecular dynamics (SMD) simulations and by the equilibrium method of adaptive biasing force (ABF) simulations. However, the ABF results qualitatively disagree with the SMD results that were based on the Jarzynski equality (JE) relating the equilibrium free-energy difference to the nonequilibrium work of the irreversible pulling experiments. In this paper, I present a new SMD study of the glycerol permeation through GlpF to explore the free-energy profile of glycerol along the permeation channel. Instead of the JE in terms of thermodynamic work, I use the fluctuation-dissipation theorem (FDT) of Brownian dynamics (BD), in terms of mechanical work, for extracting the free-energy difference from the nonequilibrium work of irreversible pulling experiments. The results of this new SMD-BD-FDT study are in agreement with the experimental data and with the ABF results.

A repulsive electrostatic mechanism for protein export through the type III secretion apparatus

Many Gram-negative bacteria initiate infections by injecting effector proteins into host cells through the type III secretion apparatus, which is comprised of a basal body, a needle, and a tip. The needle channel is formed by the assembly of a single needle protein. To explore the export mechanisms of MxiH needle protein through the needle of Shigella flexneri, an essential step during needle assembly, we have performed steered molecular dynamics simulations in implicit solvent. The trajectories reveal a screwlike rotation motion during the export of nativelike helix-turn-helix conformations. Interestingly, the channel interior with excessive electronegative potential creates an energy barrier for MxiH to enter the channel, whereas the same may facilitate the ejection of the effectors into host cells. Structurally known basal regions and ATPase underneath the basal region also have electronegative interiors. Effector proteins also have considerable electronegative potential patches on their surfaces. From these observations, we propose a repulsive electrostatic mechanism for protein translocation through the type III secretion apparatus. Based on this mechanism, the ATPase activity and/or proton motive force could be used to energize the protein translocation through these nanomachines. A similar mechanism may be applicable to macromolecular channels in other secretion systems or viruses through which proteins or nucleic acids are transported.

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.

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.

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.

Diffusion of glycerol through Escherichia coli aquaglyceroporin GlpF

The glycerol uptake facilitator, GlpF, a major intrinsic protein found in Escherichia coli, selectively conducts water and glycerol across the inner membrane. The free energy landscape characterizing the assisted transport of glycerol by this homotetrameric aquaglyceroporin has been explored by means of equilibrium molecular dynamics over a timescale spanning 0.12 micros. To overcome the free energy barriers of the conduction pathway, an adaptive biasing force is applied to the glycerol molecule confined in each of the four channels. The results illuminate the critical role played by intramolecular relaxation on the diffusion properties of the permeant. These free energy calculations reveal that glycerol tumbles and isomerizes on a timescale comparable to that spanned by its adaptive-biasing-force-assisted conduction in GlpF. As a result, reorientation and conformational equilibrium of glycerol in GlpF constitute a bottleneck in the molecular simulations of the permeation event. A profile characterizing the position-dependent diffusion of the permeant has been determined, allowing reaction rate theory to be applied for investigating conduction kinetics based on the measured free energy landscape.

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.

Structural Evidence of the Similarity of Sb(OH)3 and As(OH)3 with Glycerol: Implications for Their Uptake

Recent experimental results suggest that As(III) and Sb(III) transport in prokaryotes and eukaryotes might be facilitated by aquaglyceroporins. GlpF, the glycerol facilitator in Escherichia coli was the first to be identified as a trivalent metalloid transporter. Quantum calculations have been performed to study the possible existence of common structural properties between As(OH) 3 and Sb(OH) 3 and glycerol. Because the mechanism of substrate migration is primarily related to the successive formation of hydrogen bonds between the substrate and the hydrophilic part of the channel wall, this study has focused on the structural, thermodynamic, and electrostatic comparison of the main As(III) and Sb(III) compounds present in aqueous solution at physiological pH values, As(OH) 3 and Sb(OH) 3, with the glycerol conformation closest to the structures of these As- and Sb-containing compounds. This particular glycerol conformation has then been compared to three known experimental conformations of glycerol present in the protein channel. Besides their stoichiometry and electroneutral condition, As(OH) 3 and Sb(OH) 3 show very strong similarities to both each other and the studied conformation of the glycerol molecule: Namely, they show a similar charge distribution and a slightly smaller volume than glycerol. Their smaller size can be an additional advantage for the transit through the narrowest region of the GlpF channel. However, the metalloid hydroxyl groups lack the flexibility of glycerol, which probably helps this molecule to adapt its conformation to the topology of the GlpF channel.

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.

Molecular mechanisms of conduction and selectivity in aquaporin water channels

Aquaporins (AQP) are a family of membrane channels primarily responsible for conducting water across cellular membranes. The availability of a large body of high resolution structural data along with numerous atomic-scale simulation studies have resulted in an unprecedented level of understanding of the mechanism of function and selectivity in AQP. In this article, after summarizing major highlights of structure-functional studies of AQP, we will report on some of our recent large-scale molecular dynamics simulations investigating the mechanisms of permeation of various substances through pure lipid bilayers and through multiple pathways provided by tetrameric structures of different AQP. Comparison of the results obtained for structurally highly homologous, but functionally distinct, AQP allowed us to identify novel mechanisms of gating and selectivity of these channels and to design mutants with experimentally verified, altered properties. When applicable, special attention will be given to specific aromatic amino acids and their involvement in various functional aspects of AQP.

Molecular dynamics simulations of GlpF in a micelle vs in a bilayer: conformational dynamics of a membrane protein as a function of environment

Octyl glucoside (OG) is a detergent widely employed in structural and functional studies of membrane proteins. To better understand the nature of protein-OG interactions, molecular dynamics simulations (duration 10 ns) have been used to explore an alpha-helical membrane protein, GlpF, in OG micelles and in DMPC bilayers. Greater conformational drift of the extramembraneous protein loops, from the initial X-ray structure, is seen for the GlpF-OG simulations than for the GlpF-DMPC simulation. The mobility of the transmembrane alpha-helices is approximately 1.3x higher in the GlpF-OG than the GlpF-DMPC simulations. The detergent is seen to form an irregular torus around the protein. The presence of the protein leads to a small perturbation in the behavior of the alkyl chains in the OG micelle, namely an approximately 15% increase in the trans-gauche(-)-gauche(+) transition time. Aromatic side chains (Trp, Tyr) and basic side chains (Arg, Lys) play an important role in both protein-detergent (OG) and protein-lipid (DMPC) interactions.

A combined structural dynamics approach identifies a putative switch in factor VIIa employed by tissue factor to initiate blood coagulation

Coagulation factor VIIa (FVIIa) requires tissue factor (TF) to attain full catalytic competency and to initiate blood coagulation. In this study, the mechanism by which TF allosterically activates FVIIa is investigated by a structural dynamics approach that combines molecular dynamics (MD) simulations and hydrogen/deuterium exchange (HX) mass spectrometry on free and TF-bound FVIIa. The differences in conformational dynamics from MD simulations are shown to be confined to regions of FVIIa observed to undergo structural stabilization as judged by HX experiments, especially implicating activation loop 3 (residues 365-374{216-225}) of the so-called activation domain and the 170-loop (residues 313-322{170A-175}) succeeding the TF-binding helix. The latter finding is corroborated by experiments demonstrating rapid deglycosylation of Asn322 in free FVIIa by PNGase F but almost complete protection in the presence of TF or an active-site inhibitor. Based on MD simulations, a key switch of the TF-induced structural changes is identified as the interacting pair Leu305{163} and Phe374{225} in FVIIa, whose mutual conformations are guided by the presence of TF and observed to be closely linked to the structural stability of activation loop 3. Altogether, our findings strongly support an allosteric activation mechanism initiated by the stabilization of the Leu305{163}/Phe374{225} pair, which, in turn, stabilizes activation loop 3 and the S(1) and S(3) substrate pockets, the activation pocket, and N-terminal insertion.

How does a voltage sensor interact with a lipid bilayer? Simulations of a potassium channel domain

The nature of voltage sensing by voltage-activated ion channels is a key problem in membrane protein structural biology. The way in which the voltage-sensor (VS) domain interacts with its membrane environment remains unclear. In particular, the known structures of Kv channels do not readily explain how a positively charged S4 helix is able to stably span a lipid bilayer. Extended (2 × 50 ns) molecular dynamics simulations of the high-resolution structure of the isolated VS domain from the archaebacterial potassium channel KvAP, embedded in zwitterionic and in anionic lipid bilayers, have been used to explore VS/lipid interactions at atomic resolution. The simulations reveal penetration of water into the center of the VS and bilayer. Furthermore, there is significant local deformation of the lipid bilayer by interactions between lipid phosphate groups and arginine side chains of S4. As a consequence of this, the electrostatic field is “focused” across the center of the bilayer.

Sugar transport across lactose permease probed by steered molecular dynamics

Escherichia coli lactose permease (LacY) transports sugar across the inner membrane of the bacterium using the proton motive force to accumulate sugar in the cytosol. We have probed lactose conduction across LacY using steered molecular dynamics, permitting us to follow molecular and energetic details of lactose interaction with the lumen of LacY during its permeation. Lactose induces a widening of the narrowest parts of the channel during permeation, the widening being largest within the periplasmic half-channel. During permeation, the water-filled lumen of LacY only partially hydrates lactose, forcing it to interact with channel lining residues. Lactose forms a multitude of direct sugar-channel hydrogen bonds, predominantly with residues of the flexible N-domain, which is known to contribute a major part of LacY's affinity for lactose. In the periplasmic half-channel lactose predominantly interacts with hydrophobic channel lining residues, whereas in the cytoplasmic half-channel key protein-substrate interactions are mediated by ionic residues. A major energy barrier against transport is found within a tight segment of the periplasmic half-channel where sugar hydration is minimal and protein-sugar interaction maximal. Upon unbinding from the binding pocket, lactose undergoes a rotation to permeate either half-channel with its long axis aligned parallel to the channel axis. The results hint at the possibility of a transport mechanism, in which lactose permeates LacY through a narrow periplasmic half-channel and a wide cytoplasmic half-channel, the opening of which is controlled by changes in protonation states of key protein side groups.

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.

Mechanism of gating and ion conductivity of a possible tetrameric pore in aquaporin-1

While substrate permeation through monomeric pores of aquaporins is well characterized, little is known about the possible tetrameric pore. AQP1 has been suggested to function as an ion channel upon cGMP activation, although this idea has been controversial. Taking a theoretical and experimental approach, we demonstrate that the current might arise through the tetrameric pore and propose a plausible mechanism for conduction and gating. In response to simulated ion permeation, immediate hydration of the putative central pore was facilitated by moderate conformational changes of pore-lining residues. cGMP is found to interact with an unusually arginine-rich, cytoplasmic loop (loop D) facilitating its outward motion, which is hypothesized to trigger the opening of a cytoplasmic gate. Physiological analyses of wild-type AQP1 and a designed mutant in which two arginines of the gating loop are replaced by alanine provide experimental support for identifying a key component of the proposed mechanism.

Oleic acid is elevated in cell membranes during rapid cold-hardening and pupal diapause in the flesh fly, Sarcophaga crassipalpis

The integrity of cellular membranes is critical to the survival of insects at low temperatures, thus an advantage is conferred to insects that can adjust their composition of membrane fatty acids (FAs). Such changes contribute to homeoviscous adaption, a process that allows cellular membranes to maintain a liquid-crystalline state at temperatures that are potentially low enough to cause the membrane to enter the gel state and thereby lose its ability to maintain homeostasis. Flesh flies (Sarcophaga crassipalpis) were subjected to two experimental conditions that elicit low temperature tolerance: rapid cold-hardening and diapause. FAs were isolated and analyzed using gas chromatography-mass spectrometry. FAs changed in response to both rapid cold-hardening and diapause. In response to rapid cold-hardening (8 h at 4 degrees C), the proportion of oleic acid (18:1n-9) in pharate adults increased from 30% to 47% of the total FA pool. The proportion of almost every other FA was reduced. By entering diapause, pupae experienced an even greater increase in oleic acid proportion, to 58% of the total FA pool. Oleic acid not only promotes membrane fluidity at low temperature but also allows the cell membrane to maintain a liquid crystalline state if temperatures increase.

Aquaporin gating

An acceleration in the rate at which new aquaporin structures are determined means that structural models are now available for mammalian AQP0, AQP1, AQP2 and AQP4, bacterial GlpF, AqpM and AQPZ, and the plant SoPIP2;1. With an apparent consensus emerging concerning the mechanism of selective water transport and proton extrusion, emphasis has shifted towards the issues of substrate selectivity and the mechanisms of aquaporin regulation. In particular, recently determined structures of plant SoPIP2;1, sheep and bovine AQP0, and Escherichia coli AQPZ provide new insights into the underlying structural mechanisms by which water transport rates are regulated in diverse organisms. From these results, two distinct pictures of 'capping' and 'pinching' have emerged to describe aquaporin gating.

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.

Single-channel water permeabilities of Escherichia coli aquaporins AqpZ and GlpF

From equilibrium molecular dynamics simulations we have determined single-channel water permeabilities for Escherichia coli aquaporin Z (AqpZ) and aquaglyceroporin GlpF with the channels embedded in lipid bilayers. GlpF's osmotic water permeability constant pf exceeds by 2-3 times that of AqpZ and the diffusive permeability constant (pd) of GlpF is found to exceed that of AqpZ 2-9-fold. Achieving complete water selectivity in AqpZ consequently implies lower transport rates overall relative to the less selective, wider channel of GlpF. For AqpZ, the ratio pf/pd congruent with 12 is close to the average number of water molecules in the channel lumen, whereas for GlpF, pf/pd congruent with 4. This implies that single-file structure of the luminal water is more pronounced for AqpZ, the narrower channel of the two. Electrostatics profiles across the pore lumens reveal that AqpZ significantly reinforces water-channel interactions, and weaker water-water interactions in turn suppress water-water correlations relative to GlpF. Consequently, suppressed water-water correlations across the narrow selectivity filter become a key structural determinant for water permeation causing luminal water to permeate slower across AqpZ.

Modelling of proteins in membranes

This review describes some recent theories and simulations of mesoscopic and microscopic models of lipid membranes with embedded or attached proteins. We summarize results supporting our understanding of phenomena for which the activities of proteins in membranes are expected to be significantly affected by the lipid environment. Theoretical predictions are pointed out, and compared to experimental findings, if available. Among others, the following phenomena are discussed: interactions of interfacially adsorbed peptides, pore-forming amphipathic peptides, adsorption of charged proteins onto oppositely charged lipid membranes, lipid-induced tilting of proteins embedded in lipid bilayers, protein-induced bilayer deformations, protein insertion and assembly, and lipid-controlled functioning of membrane proteins.

Aquaporins in the hepatobiliary system

The review focuses on the potential physiological and pathophysiological roles of aquaporins (AQPs), a family of water channel proteins, in the hepatobiliary system. Among 13 aquaporins (AQP0-AQP12) cloned in mammals, seven AQPs have been identified in the liver and biliary tree. Accumulating evidence suggests that AQPs are likely involved in canalicular and ductal bile secretion, gluconeogenesis and microbial infection and may have other novel roles that affect liver function.

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.

Structural mechanism of plant aquaporin gating

Plants counteract fluctuations in water supply by regulating all aquaporins in the cell plasma membrane. Channel closure results either from the dephosphorylation of two conserved serine residues under conditions of drought stress, or from the protonation of a conserved histidine residue following a drop in cytoplasmic pH due to anoxia during flooding. Here we report the X-ray structure of the spinach plasma membrane aquaporin SoPIP2;1 in its closed conformation at 2.1 A resolution and in its open conformation at 3.9 A resolution, and molecular dynamics simulations of the initial events governing gating. In the closed conformation loop D caps the channel from the cytoplasm and thereby occludes the pore. In the open conformation loop D is displaced up to 16 A and this movement opens a hydrophobic gate blocking the channel entrance from the cytoplasm. These results reveal a molecular gating mechanism which appears conserved throughout all plant plasma membrane aquaporins.

Scalable molecular dynamics with NAMD

NAMD is a parallel molecular dynamics code designed for high-performance simulation of large biomolecular systems. NAMD scales to hundreds of processors on high-end parallel platforms, as well as tens of processors on low-cost commodity clusters, and also runs on individual desktop and laptop computers. NAMD works with AMBER and CHARMM potential functions, parameters, and file formats. This article, directed to novices as well as experts, first introduces concepts and methods used in the NAMD program, describing the classical molecular dynamics force field, equations of motion, and integration methods along with the efficient electrostatics evaluation algorithms employed and temperature and pressure controls used. Features for steering the simulation across barriers and for calculating both alchemical and conformational free energy differences are presented. The motivations for and a roadmap to the internal design of NAMD, implemented in C++ and based on Charm++ parallel objects, are outlined. The factors affecting the serial and parallel performance of a simulation are discussed. Finally, typical NAMD use is illustrated with representative applications to a small, a medium, and a large biomolecular system, highlighting particular features of NAMD, for example, the Tcl scripting language. The article also provides a list of the key features of NAMD and discusses the benefits of combining NAMD with the molecular graphics/sequence analysis software VMD and the grid computing/collaboratory software BioCoRE. NAMD is distributed free of charge with source code at www.ks.uiuc.edu.

Effects of Moisture on Glass Transition and Microstructure of Glycerol-Plasticized Soy Protein

The glass transition behavior of the glycerol-plasticized soy protein sheets (SL series) at various relative humidity (RH) was investigated by using differential scanning calorimetry with the aluminum pan and O-ring-sealed stainless steel capsule, and the microstructure of these sheets was detected on small-angle X-ray scattering. The results revealed that there were three glass transitions (Tg1, Tg2 and Tg3), corresponding to glycerol-rich, protein-rich and protein-water domains, in the protein-glycerol-water ternary system. The Tg1 values of the SL-series sheets at 75% RH decreased from -49.3 to -83.8 degrees C with an increase of glycerol content from 10 to 50 wt.-%, whereas Tg2 and Tg3 were almost invariable at about 60 degrees C and 3 degrees C, respectively. In addition, the Tg1, Tg2 and Tg3 values of the SL-25 containing 25 wt.-% glycerol at 0%, 35%, 58%, 75% and 98% RH were in the range of -12.7 - -83.2 degrees C, 65.8 - 53.1 degrees C and 3.5 - 1.9 degrees C, respectively. The result from small-angle X-ray scattering indicated that the radii of gyration (Rg) of protein-rich domain were in the range of 60-63 nm; this suggested the existence of protein macromolecules as aggregates in the stable protein-rich and protein-water domains. With an increase of RH, the tensile strength and Tg values of the SL-series sheets decreased, but the elongation at break increased. In view of the results above, the moisture in ambient environment significantly influenced the Tg values and microstructures of the glycerol-plasticized soy protein sheets, leading to the changes of the mechanical and thermal properties.

Hydroxide and proton migration in aquaporins

Hypothetical hydroxide and proton migration along the linear water chain in Aquaporin GlpF from Escherichia coli are studied by ab initio Car-Parrinello molecular dynamics simulations. It is found that the protein stabilizes a bipolar single file of water. The single file features a contiguous set of water-water hydrogen bonds in which polarization of the water molecules vary with position along the channel axis. Deprotonation of the water chain promotes the reorientation of water molecules while the hydroxide ion rapidly migrates by sequentially accepting protons from the neighboring water molecules. The hydroxide ion is not attracted by a conserved, channel-lining arginine residue, but is immobilized at two centrally located, conserved Asparagine-Proline-Alanine motifs where fourfold coordination stabilizes the ion. Hydroxide transition from the channel vestibules into the channel lumen is strongly influenced by electrostatic coupling to two conserved oppositely aligned macrodipoles. This suggests that the macrodipole's negative poles play a role in preventing hydroxide ions from entering into the channel's inner vestibules. Water protonation within the lumen facilitates water reorientation and subsequent proton expelling occurs. In the periplasmic half-channel, expelling occurs via the Grotthuss mechanism. Protonation within the cytoplasmic half-channel implies wire-breakage at the Asn-Pro-Ala motifs. The proton is here diffusively rejected as (H(5)O(2))(+).

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

The recent availability of high-resolution structures of two structurally highly homologous, but functionally distinct aquaporins from the same species, namely Escherichia coli AqpZ, a pure water channel, and GlpF, a glycerol channel, presents a unique opportunity to understand the mechanism of substrate selectivity in these channels. Comparison of the free energy profile of glycerol conduction through AqpZ and GlpF reveals a much larger barrier in AqpZ (22.8 kcal/mol) than in GlpF (7.3 kcal/mol). In either channel, the highest barrier is located at the selectivity filter. Analysis of substrate-protein interactions suggests that steric restriction of AqpZ is the main contribution to this large barrier. Another important difference is the presence of a deep energy well at the periplasmic vestibule of GlpF, which was not found in AqpZ. The latter difference can be attributed to the more pronounced structural asymmetry of GlpF, which may play a role in attracting glycerol.

Simulation studies of protein-induced bilayer deformations, and lipid-induced protein tilting, on a mesoscopic model for lipid bilayers with embedded proteins

Biological membranes are complex and highly cooperative structures. To relate biomembrane structure to their biological function it is often necessary to consider simpler systems. Lipid bilayers composed of one or two lipid species, and with embedded proteins, provide a model system for biological membranes. Here we present a mesoscopic model for lipid bilayers with embedded proteins, which we have studied with the help of the dissipative particle dynamics simulation technique. Because hydrophobic matching is believed to be one of the main physical mechanisms regulating lipid-protein interactions in membranes, we considered proteins of different hydrophobic length (as well as different sizes). We studied the cooperative behavior of the lipid-protein system at mesoscopic time- and lengthscales. In particular, we correlated in a systematic way the protein-induced bilayer perturbation, and the lipid-induced protein tilt, with the hydrophobic mismatch (positive and negative) between the protein hydrophobic length and the pure lipid bilayer hydrophobic thickness. The protein-induced bilayer perturbation was quantified in terms of a coherence length, xi(P), of the lipid bilayer hydrophobic thickness profile around the protein. The dependence on temperature of xi(P), and the protein tilt-angle, were studied above the main-transition temperature of the pure system, i.e., in the fluid phase. We found that xi(P) depends on mismatch, i.e., the higher the mismatch is, the longer xi(P) becomes, at least for positive values of mismatch; a dependence on the protein size appears as well. In the case of large model proteins experiencing extreme mismatch conditions, in the region next to the so-called lipid annulus, there appears an undershooting (or overshooting) region where the bilayer hydrophobic thickness is locally lower (or higher) than in the unperturbed bilayer, depending on whether the protein hydrophobic length is longer (or shorter) than the pure lipid bilayer hydrophobic thickness. Proteins may tilt when embedded in a too-thin bilayer. Our simulation data suggest that, when the embedded protein has a small size, the main mechanism to compensate for a large hydrophobic mismatch is the tilt, whereas large proteins react to negative mismatch by causing an increase of the hydrophobic thickness of the nearby bilayer. Furthermore, for the case of small, peptidelike proteins, we found the same type of functional dependence of the protein tilt-angle on mismatch, as was recently detected by fluorescence spectroscopy measurements.

Filter flexibility and distortion in a bacterial inward rectifier K+ channel: simulation studies of KirBac1.1

The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (<0.1 nm) filter flexibility during the concerted motion of ions and water molecules within the filter, such local changes in conformation occurring on an approximately 1-ns timescale. In the absence of K(+) ions, the KirBac selectivity filter undergoes more substantial distortions. These resemble those seen in comparable simulations of other channels (e.g., KcsA and KcsA-based homology models) and are likely to lead to functional closure of the channel. This suggests filter distortions may provide a mechanism of K-channel gating in addition to changes in the hydrophobic gate formed at the intracellular crossing point of the M2 helices. The simulation data also provide evidence for interactions of the "slide" (pre-M1) helix of KirBac with phospholipid headgroups.

Computational studies of membrane channels

The determination of the structure of several members of the K+ channel and aquaporin family represents a unique opportunity to explain the mechanism of these biomolecular systems. With their ability to go beyond static structures, molecular dynamics simulations offer a unique route for relating functional properties to membrane channel structure. The recent progress in this area is reviewed.

Lipids do influence protein function-the hydrophobic matching hypothesis revisited

A topical review of the current state of lipid-protein interactions is given with focus on the physical interactions between lipids and integral proteins in lipid-bilayer membranes. The concepts of hydrophobic matching and curvature stress are revisited in light of recent data obtained from experimental and theoretical studies which demonstrate that not only do integral proteins perturb the lipids, but the physical state of the lipids does also actively influence protein function. The case of the trans-membrane water-channel protein aquaporin GlpF from E. coli imbedded in lipid-bilayer membranes is discussed in some detail. Numerical data obtained from Molecular Dynamics simulations show on the one side that the lipid bilayer adapts to the channel by a hydrophobic matching condition which reflects the propensity of the lipid molecules for forming curved structures. On the other side, it is demonstrated that the transport function of the channel is modulated by the matching condition and/or the curvature stress in a lipid-specific manner.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

Fluctuation-driven molecular transport through an asymmetric membrane channel

Channel proteins that selectively conduct molecules across cell membranes often exhibit an asymmetric structure. By means of a stochastic model, we argue that channel asymmetry in the presence of nonequilibrium fluctuations, fueled by the cell's metabolism as observed recently, can dramatically influence the transport through such channels by a ratchetlike mechanism. For an aquaglyceroporin that conducts water and glycerol, we show that a previously determined asymmetric glycerol potential leads to enhanced inward transport of glycerol, but for unfavorably high glycerol concentrations also to enhanced outward transport that protects a cell against poisoning.

Molecular Interaction Model for the C1B Domain of Protein Kinase C-γ in the Complex with Its Activator Phorbol-12-myristate-13-acetate in Water Solution and Lipid Bilayer

Detailed molecular models of the free C1B domain of protein kinase C-gamma (PKC-gamma) and the C1B domain with its activator phorbol-12-myristate-13-acetate (PMA) in water solution and in the presence of dipalmitoylphosphatidylcholine (DPPC) bilayer are presented. Molecular dynamics of the free C1B domain reveals hydrogen bonds, which are critical for the forming of the diacylglycerols/phorbol esters binding site, and indicates the important role of Gln27 for the geometry of this site. According to the model, PMA interacts with the C1B domain by hydrophobic interactions with Pro11 and Tyr22 and by three persistent hydrogen bonds between the C3 carbonyl group of PMA and Gly23 and between the C20 hydroxyl group of PMA and the Leu21 and Thr12 residues of the C1B domain. The C9 hydroxyl group of PMA does not interact with the C1B domain, but it is involved in the interaction with the DPPC bilayer. Two preferential orientations of the C1B-PMA complex toward the DPPC bilayer resulted from our molecular modeling study.

Water and ion permeation in bAQP1 and GlpF channels: a kinetic Monte Carlo study

The kinetic Monte Carlo reaction-path-following technique is applied to determine the lowest-energy water pathway and the coordinating amino acids in bAQP1 and GlpF channels, both treated as rigid. In bAQP1, water molecules pass through the pore between the asparagine-proline-alanine (NPA) and selectivity filter (SF) sites one at a time. The water chain is interrupted at the SF where one water forms three stable hydrogen bonds with protein atoms. In this SF, water's conformation depends on the protonation locus of H182. In GlpF, two water molecules bond simultaneously to the NPA asparagines and pass through the SF in zigzag fashion. No water single-file forms in rigid GlpF. To accommodate a single file of waters requires narrowing the GlpF pore. Our results reveal that in both proteins a proposed bipolar water arrangement is thermally disrupted in the NPA region, especially in the cytoplasmic part of the pore. The equilibrium hydrogen-bonded chain is occasionally interrupted in the hydrophobic zones adjacent to the NPA motifs. The permeation of alkali cations through bAQP1 and GlpF is barred due to a large free-energy barrier in the NPA region as well as a large energy barrier blocking entry from the cytoplasm. Permeation of halides is prevented due to two large energy barriers in the cytoplasmic and periplasmic pores as well as a large free-energy barrier barring entry from the periplasm. Our results, based on modeling charge permeation, support an electrostatic rather than orientational basis for proton exclusion. Binding within the aquaporin pore cannot compensate sufficiently for dehydration of the protonic charge; there is also an electrostatic barrier in the NPA region blocking proton transport. The highly ordered single file of waters, which is drastically interrupted at the SF of bAQP1, may also contribute to proton block.

Electrostatic tuning of permeation and selectivity in aquaporin water channels

Water permeation and electrostatic interactions between water and channel are investigated in the Escherichia coli glycerol uptake facilitator GlpF, a member of the aquaporin water channel family, by molecular dynamics simulations. A tetrameric model of the channel embedded in a 16:0/18:1c9-palmitoyloleylphosphatidylethanolamine membrane was used for the simulations. During the simulations, water molecules pass through the channel in single file. The movement of the single file water molecules through the channel is concerted, and we show that it can be described by a continuous-time random-walk model. The integrity of the single file remains intact during the permeation, indicating that a disrupted water chain is unlikely to be the mechanism of proton exclusion in aquaporins. Specific hydrogen bonds between permeating water and protein at the channel center (at two conserved Asp-Pro-Ala "NPA" motifs), together with the protein electrostatic fields enforce a bipolar water configuration inside the channel with dipole inversion at the NPA motifs. At the NPA motifs water-protein electrostatic interactions facilitate this inversion. Furthermore, water-water electrostatic interactions are in all regions inside the channel stronger than water-protein interactions, except near a conserved, positively charged Arg residue. We find that variations of the protein electrostatic field through the channel, owing to preserved structural features, completely explain the bipolar orientation of water. This orientation persists despite water translocation in single file and blocks proton transport. Furthermore, we find that for permeation of a cation, ion-protein electrostatic interactions are more unfavorable at the conserved NPA motifs than at the conserved Arg, suggesting that the major barrier against proton transport in aquaporins is faced at the NPA motifs.

Glycerol conductance and physical asymmetry of the Escherichia coli glycerol facilitator GlpF

The aquaglyceroporin GlpF is a transmembrane channel of Escherichia coli that facilitates the uptake of glycerol by the cell. Its high glycerol uptake rate is crucial for the cell to survive in very low glycerol concentrations. Although GlpF allows both influx and outflux of glycerol, its structure, similar to the structure of maltoporin, exhibits a significant degree of asymmetry. The potential of mean force characterizing glycerol in the channel shows a corresponding asymmetry with an attractive vestibule only at the periplasmic side. In this study, we analyze the potential of mean force, showing that a simplified six-step model captures the kinetics and yields a glycerol conduction rate that agrees well with observation. The vestibule improves the conduction rate by 40% and 75% at 10- micro M and 10-mM periplasmic glycerol concentrations, respectively. In addition, neither the conduction rate nor the conduction probability for a single glycerol (efficiency) depends on the orientation of GlpF. GlpF appears to conduct equally well in both directions under physiological conditions.

Molecular basis of pH and Ca2+ regulation of aquaporin water permeability

Aquaporins facilitate the diffusion of water across cell membranes. We previously showed that acid pH or low Ca²⁺ increase the water permeability of bovine AQP0 expressed in Xenopus oocytes. We now show that external histidines in loops A and C mediate the pH dependence. Furthermore, the position of histidines in different members of the aquaporin family can “tune” the pH sensitivity toward alkaline or acid pH ranges. In bovine AQP0, replacement of His40 in loop A by Cys, while keeping His122 in loop C, shifted the pH sensitivity from acid to alkaline. In the killifish AQP0 homologue, MIPfun, with His at position 39 in loop A, alkaline rather than acid pH increased water permeability. Moving His39 to His40 in MIPfun, to mimic bovine AQP0 loop A, shifted the pH sensitivity back to the acid range. pH regulation was also found in two other members of the aquaporin family. Alkaline pH increased the water permeability of AQP4 that contains His at position 129 in loop C. Acid and alkaline pH sensitivity was induced in AQP1 by adding histidines 48 (in loop A) and 130 (in loop C). We conclude that external histidines in loops A and C that span the outer vestibule contribute to pH sensitivity. In addition, we show that when AQP0 (bovine or killifish) and a crippled calmodulin mutant were coexpressed, Ca²⁺ sensitivity was lost but pH sensitivity was maintained. These results demonstrate that Ca²⁺ and pH modulation are separable and arise from processes on opposite sides of the membrane.

Homology modelling and molecular dynamics simulations: comparative studies of human aquaporin-1

The structures of the mammalian water transport protein Aqp1 and of its bacterial homologue GlpF enables us to test whether homology models can be used to explore relationships between structure, dynamics and function in mammalian transport proteins. Molecular dynamics simulations (totalling almost 40 ns) were performed starting from: the X-ray structure of Aqp1; a homology model of Aqp1 based on the GlpF structure; and intermediate resolution structures of Aqp1 derived from electron microscopy. Comparisons of protein RMSDs vs. time suggest that the homology models are of comparable conformational stability to the X-ray structure, whereas the intermediate resolution structures exhibit significant conformation drift. For simulations based on the X-ray structure and on homology models, the flexibility profile vs. residue number correlates well with the crystallographic B-values for each residue. In the simulations based on intermediate resolution structures, mobility of the highly conserved NPA loops is substantially higher than in the simulations based on the X-ray structure or the homology models. Pore radius profiles remained relatively constant in the X-ray and homology model simulations but showed substantial fluctuations (reflecting the higher NPA loop mobility) in the intermediate resolution simulations. The orientation of the dipoles of water molecules within the pore is of key importance in maintaining low proton permeability through Aqp1. This property seems to be quite robust to the starting model used in the simulation. These simulations suggest that homology models based on bacterial homologues may be used to derive functionally relevant information on the structural dynamics of mammalian transport proteins.

Understanding nature's design for a nanosyringe

Synthetic and natural peptide assemblies can possess transport or conductance activity across biomembranes through the formation of nanopores. The fundamental mechanisms of membrane insertion necessary for antimicrobial or synthetic pore formation are poorly understood. We observe a lipid-assisted mechanism for passive insertion into a model membrane from molecular dynamics simulations. The assembly used in the study, a generic nanotube functionalized with hydrophilic termini, is assisted in crossing the membrane core by transleaflet lipid flips. Lipid tails occlude a purely hydrophobic nanotube. The observed insertion mechanism requirements for hydrophobic-hydrophilic matching have implications for the design of synthetic channels and antibiotics.

Molecular Basis of Proton Blockage in Aquaporins

Water transport channels in membrane proteins of the aquaporin superfamily are impermeable to ions, including H+ and OH-. We examine the molecular basis for the blockage of proton translocation through the single-file water chain in the pore of a bacterial aquaporin, GlpF. We compute the reversible thermodynamic work for the two complementary steps of the Grotthuss "hop-and-turn" relay mechanism: consecutive transfers of H+ along the hydrogen-bonded chain (hop) and conformational reorganization of the chain (turn). In the absence of H+, the strong preference for the bipolar orientation of water around the two Asn-Pro-Ala (NPA) motifs lining the pore over both unidirectional polarization states of the chain precludes the reorganization of the hydrogen-bonded network. Inversely, translocation of an excess proton in either direction is opposed by a free-energy barrier centered at the NPA region. Both hop and turn steps of proton translocation are opposed by the electrostatic field of the channel.