Capturing Subdiffusive Solute Dynamics and Predicting Selectivity in Nanoscale Pores with Time Series Modeling.

Abstract

Fitting mathematical models with a direct connection to experimental observables to the outputs of molecular simulations can be a powerful tool for extracting important physical information from them. In this study, we present two new approaches that use stochastic time series modeling to predict long-time-scale behavior and macroscopic properties from molecular simulation, which can be generalized to other molecular systems where complex diffusion occurs. In our previous work, we studied long molecular dynamics (MD) simulation trajectories of a cross-linked HII phase lyotropic liquid crystal (LLC) membrane, where we observed subdiffusive solute transport behavior characterized by intermittent hops separated by periods of entrapment. In this work, we use our models to parameterize the behavior of the same systems, so we can generate characteristic trajectory realizations that can be used to predict solute mean-squared displacements (MSDs), solute flux, and solute selectivity in macroscopic length pores. First, using anomalous diffusion theory, we show how solute dynamics can be modeled as a fractional diffusion process subordinate to a continuous time random walk. From the MD simulations, we parameterize the distribution of dwell times, hop lengths between dwells, and correlation between hops. We explore two variations of the anomalous diffusion modeling approach. The first variation applies a single set of parameters to the solute displacements and the second applies two sets of parameters based on the solute's radial distance from the closest pore center. Next, we present an approach that generalizes Markov state models, treating the configurational states of the system as a Markov process where each state has distinct transport properties. For each state and transition between states, we parameterize the distribution and temporal correlation structure of positional fluctuations as a means of characterization and to allow us to predict solute MSDs. We show that both stochastic models reasonably reproduce the MSDs calculated from MD simulations. However, qualitative differences between MD and Markov state-dependent model-generated trajectories may in some cases limit their usefulness. With these parameterized stochastic models, we demonstrate how one can estimate the flux of a solute across a macroscopic length pore and, based on these quantities, the membrane's selectivity toward each solute. This work therefore helps to connect microscopic, chemically dependent solute motions that do not follow simple diffusive behavior with long-time-scale behavior, in an approach generalizable to many types of molecular systems with complex dynamics.