Molecular dynamics simulation studies of the conformation and lateral mobility of a charged adsorbate biomolecule: implications for estimating the critical value of the radius of a pore in porous media.

The conformations, the values of the lateral transport coefficient of a charged biomolecule (desmopressin) in the adsorbed layer and in the liquid layers above the adsorbed layer, the potential energies of the interaction between the biomolecules located in different liquid layers with the charged solid surface and with the biomolecules in the adsorbed layer, the potential energies of the interaction between water molecules in the hydration layers surrounding the conformations of the biomolecules in different layers, as well as the structure and number of hydration layers between the different conformations of desmopressin, were determined by molecular dynamics simulation studies. The results show that the lateral mobility of the adsorbed desmopressin is approximately equal to zero and the value of the lateral transport coefficient of the biomolecule in the liquid layers located above the adsorbed layer increases as the distance of the liquid layer from the charged solid surface increases. But the values of the lateral transport coefficient of the biomolecule in the liquid layers above the adsorbed layer are lower in magnitude than the value of the transport coefficient of desmopressin along the direction normal to the charged solid surface in the liquid phase located above the vacant charged sites of the solid surface, and these differences in the values of the transport coefficients have important implications with respect to the replenishment of the biomolecules in the inner parts of a channel (pore), the overall rate of adsorption, and the form of the constitutive equations that would have to be used in macroscopic models to describe the mechanisms of mass transfer and adsorption in the pores of adsorbent media. Furthermore, a novel method is presented in this work that utilizes the information about the sizes of the conformations of the biomolecule in the adsorbed layer and in the liquid layers above the adsorbed layer along the direction that is normal to the charged solid surface, as well as the number and size of the hydration layers along the same direction, and could be used to estimate the value of the lower bound of the linear characteristic dimension of a pore (i.e., pore radius) in porous adsorbent media (e.g., porous adsorbent particles; skeletons of porous monoliths) in order to realize effective transport and overall adsorption rate.

Construction by molecular dynamics modeling and simulations of the porous structures formed by dextran polymer chains attached on the surface of the pores of a base matrix: characterization of porous structures.

Significant increases in the separation of bioactive molecules by using ion-exchange chromatography are realized by utilizing porous adsorbent particles in which the affinity group/ligand is linked to the base matrix of the porous particle via a polymeric extender. To study and understand the behavior of such systems, the M3B model is modified and used in molecular dynamics (MD) simulation studies to construct porous dextran layers on the surface of a base matrix, where the dextran polymer chains and the surface are covered by water. Two different porous polymer layers having 25 and 40 monomers per main polymer chain of dextran, respectively, are constructed, and their three-dimensional (3D) porous structures are characterized with respect to porosity, pore size distribution, and number of conducting pathways along the direction of net transport. It is found that the more desirable practical implications with respect to structural properties exhibited by the porous polymer layer having 40 monomers per main polymer chain, are mainly due to the higher flexibility of the polymer chains of this system, especially in the upper region of the porous structure. The characterization and analysis of the porous structures have suggested a useful definition for the physical meaning and implications of the pore connectivity of a real porous medium that is significantly different than the artificial physical meaning associated with the pore connectivity parameter employed in pore network models and whose physical limitations are discussed; furthermore, the methodology developed for the characterization of the three-dimensional structures of real porous media could be used to analyze the experimental data obtained from high-resolution noninvasive three-dimensional methods like high-resolution optical microscopy. The MD modeling and simulations methodology presented here could be used, considering that the type and size of affinity group/ligand as well as the size of the biomolecule to be adsorbed onto the affinity group/ligand are known, to construct different porous dextran layers by varying the length of the polymeric chain of dextran, the number of attachment points to the base matrix, the degree of side branching, and the number of main polymeric chains immobilized per unit surface area of base matrix. After the characterization of the porous structures of the different porous dextran layers is performed, then only a few promising structures would be selected for studying the immobilization of adsorption sites on the pore surfaces and the subsequent adsorption of the bioactive molecules onto the immobilized affinity groups/ligands.

Molecular dynamics simulation studies of the transport and adsorption of a charged macromolecule onto a charged adsorbent solid surface immersed in an electrolytic solution.

Molecular dynamics simulations were performed in order to study the transport and adsorption of a charged macromolecule (desmopressin) onto a charged solid surface in an electrolytic solution. The strong Coulombic interaction from the charged solid surface represents the major force for accelerating, orienting, entrapping in the electrical double layer, and adsorbing the macromolecule onto the charged solid surface. The macromolecule is flattened as it approaches the charged surface, giving rise to a stronger surface exclusion effect that shields surface sites. When adsorbed, the macromolecule is restrained by a surface interaction more than one hundred times stronger than the thermal energy, of which 99.8% results from the strong dominant Coulombic interaction, and trapped by a hydration layer adjacent to the surface. This leads to zero lateral displacement of the adsorbed macromolecule and indicates that surface diffusion is a physically implausible mechanism in similar systems. Explicit solvent is required for realistic representation of the macromolecular structure and the surface interaction energy. The adsorbed macromolecule also decreased the electrostatic potential gradient perpendicular to the charged solid surface and introduced additional electrostatic potential gradients laterally. The results obtained from the molecular dynamics simulations confirm the importance of electrophoretic migration and support the physical mechanisms used in a macroscopic continuum model that predicts an overshoot in the concentration of a charged macromolecule in the adsorbed phase under certain conditions of pH and ionic strength.

Analysis and parametric sensitivity of the behavior of overshoots in the concentration of a charged adsorbate in the adsorbed phase of charged adsorbent particles: practical implications for separations of charged solutes.

In this work, an analysis of the parametric sensitivity of the overshoot in the concentration of the adsorbate in the adsorbed phase, which occurs under certain conditions during an ion-exchange adsorption process, is presented and used to suggest practical implications of the concentration overshoot phenomenon on operational policies and configurations of chromatographic columns and finite bath adsorption systems. The results presented in this work demonstrate and explain how the development of an overshoot in the concentration of the adsorbate in the adsorbed phase could be enhanced or suppressed by (i) varying the diffusion coefficient, D3, of the adsorbate relative to the diffusion coefficients, D1 and D2, of the cations and anions, respectively, of the background/buffer electrolyte, (ii) altering the initial surface charge density, delta0, of the charged adsorbent particles, (iii) varying the Debye length, lambda, and (iv) changing the initial concentration, Cd3(0), of the adsorbate in the bulk liquid of the finite bath. The influence of the pH and ionic strength, Iinfinity, of the liquid solution on the development of an overshoot in the concentration of the adsorbate in the adsorbed phase is also presented and discussed through the relationships of these parameters to delta0 and lambda, respectively. Furthermore, a detailed explanation of the effects of each parameter on the interplay between the diffusive and electrophoretic molar fluxes, as well as on the structure and functioning of the electrical double layer, which are responsible for the concentration overshoot phenomenon, is presented.