In Silico Determination of Gas Permeabilities by Non-Equilibrium Molecular Dynamics: CO2 and He through PIM-1.
Ontology highlight
ABSTRACT:
SUBMITTER: Frentrup H
PROVIDER: S-EPMC4384093 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description: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).
| S-EPMC5493307 | biostudies-literature
Project description:Island biotas emerge from the interplay between colonisation, speciation and extinction and are often the scene of spectacular adaptive radiations. A common assumption is that insular diversity is at a dynamic equilibrium, but for remote islands, such as Hawaii or Galápagos, this idea remains untested. Here, we reconstruct the temporal accumulation of terrestrial bird species of the Galápagos using a novel phylogenetic method that estimates rates of biota assembly for an entire community. We show that species richness on the archipelago is in an ascending phase and does not tend towards equilibrium. The majority of the avifauna diversifies at a slow rate, without detectable ecological limits. However, Darwin's finches form an exception: they rapidly reach a carrying capacity and subsequently follow a coalescent-like diversification process. Together, these results suggest that avian diversity of remote islands is rising, and challenge the mutual exclusivity of the non-equilibrium and equilibrium ecological paradigms.
| S-EPMC4745040 | biostudies-literature
Project description:Grand canonical Monte Carlo (GCMC) simulations are widely used with equilibrium molecular dynamics (EMD) to predict gas adsorption and diffusion in single-crystals of metal-organic frameworks (MOFs). Adsorption and diffusion data obtained from these simulations are then combined to predict gas permeabilities and selectivities of MOF membranes. This GCMC + EMD approach is highly useful to screen a large number of MOFs for a target membrane-based gas separation process. External field non-equilibrium molecular dynamics (NEMD) simulations, on the other hand, can directly compute gas permeation by providing an accurate representation of MOF membranes but they are computationally demanding and require long simulation times. In this work, we performed NEMD simulations to investigate H2/CH4 separation performances of MOF membranes. Both single-component and binary mixture permeabilities of H2 and CH4 were computed using the NEMD approach and results were compared with the predictions of the GCMC + EMD approach and experimental measurements reported in the literature. Our results showed that there is a good agreement between NEMD simulations and experiments for the permeability and selectivity of MOF membranes. NEMD simulations provided the direct observation of the mass transfer resistances on the pore mouth of MOF membranes, which is neglected in the GCMC + EMD approach. Our results suggested that once the very large numbers of MOF materials were screened using the GCMC + EMD approach, more detailed NEMD calculations can be performed for the best membrane candidates to unlock the actual gas transport mechanism before the experimental fabrication of MOF membranes.
| S-EPMC6395021 | biostudies-literature
Project description:To support photosynthetic CO2 fixation by Rubisco, the chloroplast must be fed with inorganic carbon in the form of CO2 or bicarbonate. However, the mechanisms allowing the rapid passage of this gas and this charged molecule through the bounding membranes of the chloroplast envelope are not yet completely elucidated. We describe here a method allowing us to measure the permeability of these two molecules through the chloroplast envelope using a membrane inlet mass spectrometer and 18O-labelled inorganic carbon. We established that the internal stromal carbonic anhydrase activity is not limiting for this technique, and precisely measured the chloroplast surface area and permeability values for CO2 and bicarbonate. This was performed on chloroplasts from several plant species, with values ranging from 2.3 × 10-4 m s-1 to 8 × 10-4 m s-1 permeability for CO2 and 1 × 10-8 m s-1 for bicarbonate. We were able to apply our method to chloroplasts from an Arabidopsis aquaporin mutant, and this showed that CO2 permeability was reduced 50% in the mutant compared with the wild-type reference.
| S-EPMC5853524 | biostudies-literature
Project description:Time-lapse microscopy imaging provides direct access to the dynamics of soft and living systems. At mesoscopic scales, such microscopy experiments reveal intrinsic thermal and non-equilibrium fluctuations. These fluctuations, together with measurement noise, pose a challenge for the dynamical analysis of these Brownian movies. Traditionally, methods to analyze such experimental data rely on tracking embedded or endogenous probes. However, it is in general unclear, especially in complex many-body systems, which degrees of freedom are the most informative about their non-equilibrium nature. Here, we introduce an alternative, tracking-free approach that overcomes these difficulties via an unsupervised analysis of the Brownian movie. We develop a dimensional reduction scheme selecting a basis of modes based on dissipation. Subsequently, we learn the non-equilibrium dynamics, thereby estimating the entropy production rate and time-resolved force maps. After benchmarking our method against a minimal model, we illustrate its broader applicability with an example inspired by active biopolymer gels.
| S-EPMC7585442 | biostudies-literature
Project description:A popular approach to analyze the dynamics of high-dimensional many-body systems, such as macromolecules, is to project the trajectories onto a space of slowly varying collective variables, where subsequent analyses are made, such as clustering or estimation of free energy profiles or Markov state models. However, existing "dynamical" dimension reduction methods, such as the time-lagged independent component analysis (TICA), are only valid if the dynamics obeys detailed balance (microscopic reversibility) and typically require long, equilibrated simulation trajectories. Here, we develop a dimension reduction method for non-equilibrium dynamics based on the recently developed Variational Approach for Markov Processes (VAMP) by Wu and Noé. VAMP is illustrated by obtaining a low-dimensional description of a single file ion diffusion model and by identifying long-lived states from molecular dynamics simulations of the KcsA channel protein in an external electrochemical potential. This analysis provides detailed insights into the coupling of conformational dynamics, the configuration of the selectivity filter, and the conductance of the channel. We recommend VAMP as a replacement for the less general TICA method.
| S-EPMC6486394 | biostudies-literature
Project description:CO2 enhanced oil recovery (CO2-EOR) has become significantly crucial to the petroleum industry, in particular, CO2 miscible flooding can greatly improve the efficiency of EOR. Minimum miscibility pressure (MMP) is a vital factor affecting CO2 flooding, which determines the yield and economic benefit of oil recovery. Therefore, it is important to predict this property for a successful field development plan. In this study, a novel model based on molecular dynamics to determine MMP was developed. The model characterized a miscible state by calculating the ratio of CO2 and crude oil atoms that pass through the initial interface. The whole process was not affected by other external objective factors. We compared our model with several famous empirical correlations, and obtained satisfactory results-the relative errors were 8.53% and 13.71% for the two equations derived from our model. Furthermore, we found the MMPs predicted by different reference materials (i.e., CO2/crude oil) were approximately linear (R2 = 0.955). We also confirmed the linear relationship between MMP and reservoir temperature (TR). The correlation coefficient was about 0.15 MPa/K in the present study.
| S-EPMC8401628 | biostudies-literature
Project description:BackgroundWe consider a focal adhesion to be made up of molecular complexes, each consisting of a ligand, an integrin molecule, and associated plaque proteins. Free energy changes drive the binding and unbinding of these complexes and thereby controls the focal adhesion's dynamic modes of growth, treadmilling and resorption.Principal findingsWe have identified a competition among four thermodynamic driving forces for focal adhesion dynamics: (i) the work done during the addition of a single molecular complex of a certain size, (ii) the chemical free energy change associated with the addition of a molecular complex, (iii) the elastic free energy change associated with deformation of focal adhesions and the cell membrane, and (iv) the work done on a molecular conformational change. We have developed a theoretical treatment of focal adhesion dynamics as a nonlinear rate process governed by a classical kinetic model. We also express the rates as being driven by out-of-equilibrium thermodynamic driving forces, and modulated by kinetics. The mechanisms governed by the above four effects allow focal adhesions to exhibit a rich variety of behavior without the need to introduce special constitutive assumptions for their response. For the reaction-limited case growth, treadmilling and resorption are all predicted by a very simple chemo-mechanical model. Treadmilling requires symmetry breaking between the ends of the focal adhesion, and is achieved by driving force (i) above. In contrast, depending on its numerical value (ii) causes symmetric growth, resorption or is neutral, (iii) causes symmetric resorption, and (iv) causes symmetric growth. These findings hold for a range of conditions: temporally-constant force or stress, and for spatially-uniform and non-uniform stress distribution over the FA. The symmetric growth mode dominates for temporally-constant stress, with a reduced treadmilling regime.SignificanceIn addition to explaining focal adhesion dynamics, this treatment can be coupled with models of cytoskeleton dynamics and contribute to the understanding of cell motility.
| S-EPMC2923603 | biostudies-literature
Project description:X-ray photon correlation spectroscopy (XPCS) is a powerful tool in the investigation of dynamics covering a broad time and length scale. It has been widely used to probe dynamics for systems in both equilibrium and non-equilibrium states; in particular, for systems undergoing a phase transition where the structural growth kinetics and the microscopic dynamics are strongly intertwined. The resulting time-dependent dynamic behavior can be described using the two-time correlation function (TTC), which, however, often contains more interesting features than the component along the diagonal, and cannot be easily interpreted via the classical simulation methods. Here, a reverse engineering (RE) approach is proposed based on particle-based heuristic simulations. This approach is applied to an XPCS measurement on a protein solution undergoing a liquid-liquid phase separation. It is demonstrated that the rich features of experimental TTCs can be well connected with the key control parameters including size distribution, concentration, viscosity and mobility of domains. The dynamic information obtained from this RE analysis goes beyond the existing theory. The RE approach established in this work is applicable for other processes such as film growth, coarsening or evolving systems.
| S-EPMC9252156 | biostudies-literature
Project description:Why are some traits and trait combinations exceptionally common across the tree of life, whereas others are vanishingly rare? The distribution of trait diversity across a clade at any time depends on the ancestral state of the clade, the rate at which new phenotypes evolve, the differences in speciation and extinction rates across lineages, and whether an equilibrium has been reached. Here we examine the role of transition rates, differential diversification (speciation minus extinction) and non-equilibrium dynamics on the evolutionary history of angiosperms, a clade well known for the abundance of some trait combinations and the rarity of others. Our analysis reveals that three character states (corolla present, bilateral symmetry, reduced stamen number) act synergistically as a key innovation, doubling diversification rates for lineages in which this combination occurs. However, this combination is currently less common than predicted at equilibrium because the individual characters evolve infrequently. Simulations suggest that angiosperms will remain far from the equilibrium frequencies of character states well into the future. Such non-equilibrium dynamics may be common when major innovations evolve rarely, allowing lineages with ancestral forms to persist, and even outnumber those with diversification-enhancing states, for tens of millions of years.
| S-EPMC4874697 | biostudies-literature