Molecular Dynamics Simulation of the Effect of Angle Variation on Water Permeability through Hourglass-Shaped Nanopores.
Ontology highlight
ABSTRACT:
SUBMITTER: Tang D
PROVIDER: S-EPMC5458919 | biostudies-other | 2015 Oct
REPOSITORIES: biostudies-other
Similar Datasets
Project description:This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.
| S-EPMC7517714 | biostudies-literature
Project description:The integration of local heat sources with solid-state nanopores offers new means for controlling the transmembrane transport of charged biomacromolecules. In the case of electrophoretic transport of DNA, recent experimental studies revealed unexpected temperature dependences of the DNA capture rate, the DNA translocation velocity, and the ionic current blockades produced by the presence of DNA in the nanopore. Here, we report the results of all-atom molecular dynamics simulations that elucidated the effect of temperature on the key microscopic processes governing electric field-driven transport of DNA through nanopores. Mimicking the experimental setup, we simulated the capture and subsequent translocation of short DNA duplexes through a locally heated nanopore at several temperatures and electrolyte conditions. The temperature dependence of ion mobility at the DNA surface was found to cause the dependence of the relative conductance blockades on temperature. To the first order, the effective force on DNA in the nanopore was found to be independent of temperature, despite a considerable reduction of solution viscosity. The temperature dependence of the solution viscosity was found to make DNA translocations faster for a uniformly heated system but not in the case of local heating that does not affect viscosity of solution surrounding the untranslocated part of the molecule. Increasing solution temperature was also found to reduce the lifetime of bonds formed between cations and DNA. Using a flow suppression algorithm, we were able to separate the effects of electro-osmotic flow and direct ion binding, finding the reduced durations of DNA-ion bonds to increase, albeit weakly, the effective force experienced by DNA in an electric field. Unexpectedly, our simulations revealed a considerable temperature dependence of solvent velocity at the DNA surface-slip velocity, an effect that can alter hydrodynamic coupling between the motion of DNA and the surrounding fluid.
| S-EPMC4880514 | biostudies-literature
Project description:The alpha-hemolysin (AHL) nanochannel is a non-selective channel that allows for uncontrolled transport of small molecules across membranes leading to cell death. Although it is a bacterial toxin, it has promising applications, ranging from drug delivery systems to nano-sensing devices. This study focuses on the transport of water molecules through an AHL nanochannel using molecular dynamics (MD) simulations. Our results show that AHL can quickly transport water across membranes. The first-passage time approach was used to estimate the diffusion coefficient and the mean exit time. To study the energetics of transport, the potential of mean force (PMF) of a water molecule along the AHL nanochannel was calculated. The results show that the energy barriers of water permeation across a nanopore are always positive along the channel and the values are close to thermal energy (kBT). These findings suggest that the observed quick permeation of water is due to small energy barriers and a hydrophobic inner channel surface resulting in smaller friction. We speculate that these physical mechanisms are important in how AHL causes cell death.
| S-EPMC4713412 | biostudies-literature
Project description:Engineered protein nanopores, such as those based on α-hemolysin from Staphylococcus aureus have shown great promise as components of next-generation DNA sequencing devices. However, before such protein nanopores can be used to their full potential, the conformational dynamics and translocation pathway of the DNA within them must be characterized at the individual molecule level. Here, we employ atomistic molecular dynamics simulations of single-stranded DNA movement through a model α-hemolysin pore under an applied electric field. The simulations enable characterization of the conformations adopted by single-stranded DNA, and allow exploration of how the conformations may impact on translocation within the wild-type model pore and a number of mutants. Our results show that specific interactions between the protein nanopore and the DNA can have a significant impact on the DNA conformation often leading to localized coiling, which in turn, can alter the order in which the DNA bases exit the nanopore. Thus, our simulations show that strategies to control the conformation of DNA within a protein nanopore would be a distinct advantage for the purposes of DNA sequencing.
| S-EPMC3433622 | biostudies-literature
Project description:Aquaporin (AQP) functions as a water-conducting pore. Mercury inhibits the water permeation through AQP. Although site-directed mutagenesis has shown that mercury binds to Cys189 during the inhibition process, it is not fully understood how this inhibits the water permeation through AQP1. We carried out 40 ns molecular dynamics simulations of bovine AQP1 tetramer with mercury (Hg-AQP1) or without mercury (Free AQP1). In Hg-AQP1, Cys191 (Cys189 in human AQP1) is converted to Cys-SHg+ in each monomer. During each last 10 ns, we observed water permeation events occurred 23 times in Free AQP1 and never in Hg-AQP1. Mercury binding did not influence the whole structure, but did induce a collapse in the orientation of several residues at the ar/R region. In Free AQP1, backbone oxygen atoms of Gly190, Cys191, and Gly192 lined, and were oriented to, the surface of the water pore channel. In Hg-AQP1, however, the SHg+ of Cys191-SHg+ was oriented toward the outside of the water pore. As a result, the backbone oxygen atoms of Gly190, Cys191, and Gly192 became disorganized and the ar/R region collapsed, thereby obstructing the permeation of water. We suggest that mercury disrupts the water pore of AQP1 through local conformational changes in the ar/R region.
| S-EPMC2856166 | biostudies-literature
Project description:The anisotropic diffusive behavior of nanoparticles with complex shapes attracts great interest due to its potential applications in many fields ranging from bionics to aeronautic industry. Although molecular dynamics (MD) simulations are used widely to investigate nanoparticle diffusion properties, universal methods to describe the diffusion process comprehensively are still lacking. Here, we address this problem by introducing diffusion tensor as it can describe translational and rotational diffusion in three dimensions both individually and their coupling. We take carbon triple sphere suspended in argon fluid as our model system. The consistency of our results and velocity autocorrelation function(VAF) method validates our simulations. The coupling between translational and rotational diffusion is observed directly from analyzing diffusion tensor, and quantified by coupling diffusion coefficient. Our simulation reveals non-trivial effect of some factors in diffusion at nanoscale, which was not considered in previous theories. In addition to introducing an effective method to calculate the diffusion tensor in MD simulations, our work also provides insights for understanding the diffusion process of arbitrary-shaped particles in nanoengineering.
| S-EPMC6908649 | biostudies-literature
Project description:The porous structure and mass transport characteristics of disordered silicate porous media were investigated via a geometry based analysis of water confined in the pores. Disordered silicate porous media were constructed to mimic the dissolution behavior of an alkali aluminoborosilicate glass, i.e., soluble Na and B were removed from the bulk glass, and then water molecules and Na were introduced into the pores to provide a complex porous structure filled with water. This modelling approach revealed large surface areas of disordered porous media. In addition, a number of isolated water molecules were observed in the pores, despite accessible porous connectivity. As the fraction of mobile water was approximately 1%, the main water dynamics corresponded to vibrational motion in a confined space. This significantly reduced water mobility was due to strong hydrogen-bonding water-surface interactions resulting from the large surface area. This original approach provides a method for predicting the porous structure and water transport characteristics of disordered silicate porous media.
| S-EPMC5830603 | biostudies-literature
Project description: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.
| S-EPMC1896254 | biostudies-literature
Project description:Ion transports through ion channels, biological nanopores, are essential for life: Living cells generate electrical signals by utilizing ion permeation through channels. The measured current-voltage (i-V) relations through most ion channels are sublinear, however, its physical meaning is still elusive. Here we calculated the i-V curves through anion-doped carbon nanotubes, a model of an ion channel, using molecular dynamics simulation. It was found the i-V curve reflects the physical origin of the rate-determining step: the i-V curve is sublinear when the permeation is entropy bottlenecked, while it is superlinear in the case of the energy bottlenecked permeation. Based on this finding, we discuss the relation between the molecular mechanism of ion permeation through the biological K(+) channels and the shape of the i-V curves through them. This work also provides a clue for a novel design of nanopores that show current rectification.
| S-EPMC4863172 | biostudies-other
Project description:Recently, molecular dynamics (MD) simulations have been utilized to investigate the barrier properties of human skin stratum corneum (SC) lipid bilayers. Different MD methods and force fields have been utilized, with predicted permeabilities varying by few orders of magnitude. In this work, we compare constrained MD simulations with restrained MD simulations to obtain the potential of the mean force and the diffusion coefficient profile for the case of a water molecule permeating across an SC lipid bilayer. Corresponding permeabilities of the simulated lipid bilayer are calculated via the inhomogeneous solubility diffusion model. Results show that both methods perform similarly, but restrained MD simulations have proven to be the more robust approach for predicting the potential of the mean force profile. Critical to both methods are the sampling of the whole trans-bilayer axis and the following symmetrization process. Re-analysis of the previously reported free energy profiles showed that some of the discrepancies in the reported permeability values is due to misquotation of units, while some are due to the inaccurately obtained potential of the mean force. By using the existing microscopic geometrical models via the intercellular lipid pathway, the permeation through the whole SC is predicted from the MD simulation results, and the predicted barrier properties have been compared to experimental data from the literature with good agreement.
| S-EPMC8717400 | biostudies-literature