Search

In order to elucidate how phosphatidylserine (PS6) interacts with AQP5 in a cell membrane, we developed a hybrid steered molecular dynamics (hSMD) method that involved: (1) Simultaneously steering two centers of mass of two selected segments of the ligand, and (2) equilibrating the ligand-protein complex with and without biasing the system. Validating hSMD, we first studied vascular endothelial growth factor receptor 1 (VEGFR1) in complex with N-(4-Chlorophenyl)-2-((pyridin-4-ylmethyl)amino)benzamide (8ST), for which the binding energy is known fromin vitroexperiments. In this study, our computed binding energy well agreed with the experimental value. Knowing the accuracy of this hSMD method, we applied it to the AQP5-lipid-bilayer system to answer an outstanding question relevant to AQP5’s physiological function: Will the PS6, a lipid having a single long hydrocarbon tail that was found in the central pore of the AQP5 tetramer crystal, actually bind to and inhibit AQP5’s central pore under near-physiological conditions, namely, when AQP5 tetramer is embedded in a lipid bilayer? We found,in silico, using the CHARMM 36 force field, that binding PS6 to AQP5 was a factor of 3 million weaker than “binding” it in the lipid bilayer. This suggests that AQP5’s central pore will not be inhibited by PS6 or a similar lipid in a physiological environment. [ABSTRACT FROM PUBLISHER]

Abstract For its fundamental relevance, transport of water and glycerol across the erythrocyte membrane has long been investigated before and after the discovery of aquaporins (AQPs), the membrane proteins responsible for water and glycerol transport. AQP1 is abundantly expressed in the human erythrocyte for maintaining its hydrohomeostasis where AQP3 is also expressed (at a level ~30-folds lower than AQP1) facilitating glycerol transport. This research is focused on two of the remaining questions: How permeable is AQP3 to water? What is the glycerol-AQP3 affinity under near-physiological conditions? Through atomistic modelling and large-scale simulations, we found that AQP3 is two to three times more permeable to water than AQP1 and that the glycerol-AQP3 affinity is approximately 500/M. Using these computed values along with the data from the latest literature on AQP1 and on erythrocyte proteomics, we estimated the water and glycerol transport rates across the membrane of an entire erythrocyte. We used these rates to predict the time courses of erythrocyte swelling-shrinking in response to inward and outward osmotic gradients. Experimentally, we monitored the time course of human erythrocytes when subject to an osmotic or glycerol gradient with light scattering in a stopped-flow spectrometer. We observed close agreement between the experimentally measured and the computationally predicted time courses of erythrocytes, which corroborated our computational conclusions on the AQP3 water-permeability and the glycerol-AQP3 affinity. Graphical abstract Unlabelled Image Highlights • Human Aquaporin 3 is functionally similar to its bacterium and parasite analogues. • Aquaporin 3, a glycerol channel, is more permeable than Aquaporin 1. • Aquaporin 3 has an affinity of 500/M for glycerol. [ABSTRACT FROM AUTHOR]

Measuring or computing the single-channel permeability of aquaporins/aquaglyceroporins (AQPs) has long been a challenge. The measured values scatter over an order of magnitude but the corresponding Arrhenius activation energies converge in the current literature. Osmotic flux through an AQP was simulated as water current forced through the channel by kilobar hydraulic pressure or theoretically approximated as single-file diffusion. In this paper, we report large scale simulations of osmotic current under sub M gradient through three AQPs (water channels AQP4 and AQP5 and glycerol-water channel GlpF) using the mature particle mesh Ewald technique (PME) for which the established force fields have been optimized with known accuracy. These simulations were implemented with hybrid periodic boundary conditions devised to avoid the artifactitious mixing across the membrane in a regular PME simulation. The computed single-channel permeabilities at 5 °C and 25 °C are in agreement with recently refined experiments on GlpF. The Arrhenius activation energies extracted from our simulations for all the three AQPs agree with the in vitro measurements. The single-file diffusion approximations from our large-scale simulations are consistent with the current literature on smaller systems. From these unambiguous agreements among the in vitro and in silico studies, we observe the quantitative accuracy of the all-atom force fields of the current literature for water-channel biology. We also observe that AQP4, that is particularly rich in the central nervous system, is more efficient in water conduction and more temperature-sensitive than other water-only channels (excluding glycerol channels that also conduct water when not inhibited by glycerol). [ABSTRACT FROM AUTHOR]