Abstract:The aquaglyceroporin GlpF is a member of the aquaporin family. It selectively conducts small molecules, such as glycerol, across the cell membrane under a concentration gradient of the substrate. Atomistic molecular dynamics (MD) simulation would provide great insight into the substrate transport mechanism of GlpF and membrane channels alike. Ideally, non-equilibrium simulations under various concentration gradients of glycerol are desired to emulate the transportation in cells, but this kind of simulation is difficult due to a complicated system setup and high computational cost. Here, we present a new strategy to extract non-equilibrium kinetic information from equilibrium MD simulation. We first performed long-time (totally 22.5 μs) multi-copy equilibrium MD simulations of glycerol conduction through GlpF. Tens of times the spontaneous permeation of glycerol through GlpF was observed, allowing us to elucidate the detailed mechanism of the stereoselectivity for glycerol. Then we employed Markov state model (MSM) analysis of the MD trajectories to identify the intermediate states during glycerol transport and calculate the inter-state transition rate constants. Based on the results of MSM analysis, we built the kinetic models of glycerol transport and calculated the glycerol fluxes under various concentration gradients by solving the master equations. The results agree well with the experimental measurement at a certain glycerol concentration, and provide holistic information on the glycerol conduction capacity of GlpF. Our work demonstrates that long-time atomistic MD simulations can now bridge the microscopic dynamics and the kinetic description of substance transport through membrane channels, hopefully facilitating the engineering of new selective channels for various molecules.

Evidence for a Permanent Single-Event Upset Mechanism

Anomalous behaviour of transport properties in a supercooled water–glycerol mixture

Anomalous behaviour of transport properties in a supercooled water-glycerol mixture

Molecular velocity auto-correlations in glycerol/water mixtures studied by NMR MGSE method

Density scaling and isodynes in glycerol-water mixtures

A molecular dynamics study of cryoprotective agent – Water–sodium chloride ternary solutions

Molecular Dynamics Simulation Studies of Cryoprotective Agent Solutions: the Relation Between Melting Temperature and the Ratio of Hydrogen Bonding Acceptor to Donor Number

Orientational dynamics in supercooled glycerol computed from MD simulations: self and cross contributions

Terahertz dynamics in the glycerol-water system

Glycerol transport through the aquaglyceroporin GlpF: bridging dynamics and kinetics with atomic simulation.

Self-Diffusion of Glycerol in $γ$-Alumina Nanopores: Understanding the Effect of Pore Saturation on the Dynamics of Confined Polyalcohols

Slow rheological mode in glycerol and glycerol-water mixtures

Effect of temperature and glycerol on the hydrogen-bond dynamics of water

Molecular packing, hydrogen bonding, and fast dynamics in lysozyme/trehalose/glycerol and trehalose/glycerol glasses at low hydration

Supercritical density fluctuations and structural heterogeneity in supercooled water-glycerol microdroplets

Interfacial properties of binary azeotropic mixtures of simple fluids: Molecular dynamics simulation and density gradient theory.

Microscopic structure and interaction analysis for supercritical carbon dioxide-ethanol mixtures: a Monte Carlo simulation study.

Temperature dependent structure and dynamics of isopropanol-water liquid mixtures at low alcohol content

Liquid dynamics in partially crystalline glycerol

Auditory Evoked Potentials in the Assessment of Central Nervous System Effects of Antimigraine Drugs