Rigorous study of molecular dynamics of a single dsDNA confined in a nanochannel: Introduction of a critical mobility behaviour

A novel molecular dynamics study of CO2 permeation through aquaporin-5

Aquaporins (AQPs) are protein channels which facilitate rapid water permeation across cell membrane. The AQPs are very vital for biological organs, as their malfunction causes severe diseases in human body. A particular family of AQPs, that is AQP5, has a significant role in lung fluid transport due to submucosal glands structure. However, it has not been yet well understood whether these protein channels can conduct gas molecules. Here, Molecular Dynamics (MD) simulations are used to investigate the CO₂ permeability and diffusion in AQP5 during a 40-nanosecond period. For the first time, equilibrium and Steered MD (SMD) are used to simulate self and force-induced diffusion of CO₂ molecules across AQP5 and POPE lipid bilayer. According to PMFs profile associated to CO₂ permeation, the hydrophobic central pore provides a more suitable pathway for gas molecules compared to other AQP5 channels. Although CO₂ molecules can also permeate across AQP5 water channels, the rate of CO₂ permeation through four channels of the AQP5 monomers is much lower than the central pore. The rate of CO₂ permeation through four AQP5 water channels is even lower than CO₂ diffusion through POPE lipid membrane. The results reported in this investigation demonstrate that MD simulations of human AQP5 provide valuable insights into the gas permeation mechanism for both the equilibrium self-diffusion, and quasi-equilibrium condition.

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.

Monte Carlo simulation algorithm for B-DNA

Understanding the structure-function relationship of biomolecules containing DNA has motivated experiments aimed at determining molecular structure using methods such as small-angle X-ray and neutron scattering (SAXS and SANS). SAXS and SANS are useful for determining macromolecular shape in solution, a process which benefits by using atomistic models that reproduce the scattering data. The variety of algorithms available for creating and modifying model DNA structures lack the ability to rapidly modify all-atom models to generate structure ensembles. This article describes a Monte Carlo algorithm for simulating DNA, not with the goal of predicting an equilibrium structure, but rather to generate an ensemble of plausible structures which can be filtered using experimental results to identify a sub-ensemble of conformations that reproduce the solution scattering of DNA macromolecules. The algorithm generates an ensemble of atomic structures through an iterative cycle in which B-DNA is represented using a wormlike bead-rod model, new configurations are generated by sampling bend and twist moves, then atomic detail is recovered by back mapping from the final coarse-grained configuration. Using this algorithm on commodity computing hardware, one can rapidly generate an ensemble of atomic level models, each model representing a physically realistic configuration that could be further studied using molecular dynamics. © 2016 Wiley Periodicals, Inc.

Numerical investigation of molecular nano-array in potential-energy profile for a single dsDNA

A Rigorous numerical investigation on dsDNA translocation in quasi-2-dimensional nano-array filter is performed using Molecular Dynamics (MD) method. Various dsDNA molecules with different sizes are chosen in order to model Ogston sieving in a nano-array filter. The radius of gyration of dsDNA molecule is less than the characteristic length of the shallow region in nano-array. The dsDNA molecule is assumed to be in the 0.05M NaCl electrolyte. MD shows acceptable results for potential-energy profile for nano-array filter. According to the MD outcomes, the dsDNA electrophoretic mobility decreases almost linearly with dsDNA size and show the same trend as Ogston sieving for gel electrophoresis. In addition, different shapes for nano-array filter are studied for a unique dsDNA molecule. It is concluded that steeping the nano-array wall can cause the retardation of dsDNA translocation and decreases dsDNA electrophoretic mobility.