Jee-Ching Wang

Porous Polymer Adsorbent Media Constructed by Molecular Dynamics Modeling and Simulations: the Immobilization of Charged Ligands and Their Effect on Pore Structure and Local Nonelectroneutrality

A molecular dynamics modeling and simulation approach is presented and employed to construct porous dextran polymer ion-exchange adsorbent media. Both the activation step of the surface of the pores of the dextran polymer layer grafted on an agarose surface and the immobilization of charged ligands on the activated surface of the porous dextran polymer layer are considered. For the systems studied in this work, the activation step modifies slightly the pore structure of the base, nonactivated porous dextran polymer, while the immobilization of the ligands on the activated pore surface of the dextran layer changes significantly the pore structure of the activated dextran layer. The density distributions of the counterions and immobilized charged ligands along the direction of net transport in the adsorbent media constructed in this study are found to be nonuniform. The variables that affect the shape and magnitude of the density distributions of the counterions and immobilized charged ligands as well as the total number of charged ligands that can be immobilized on the activated porous dextran layer are identified and presented in this work. Furthermore, the data clearly show that there is local nonelectroneutrality in the porous dextran polymer ion-exchange adsorbent media, and this result has very important practical implications for the operation and performance of separation systems involving ion-exchange adsorbent media (e.g., ion-exchange chromatography systems). Also, the results of this work suggest approaches for (1) controlling the immobilization process of charged ligands and (2) constructing and studying the behavior of chromatographic polymeric monoliths and packed bed columns having a gradient of density of functionalities along the axis of the chromatographic polymeric monolith or packed bed column.

A molecular dynamics study on the transport of a charged biomolecule in a polymeric adsorbent medium and its adsorption onto a charged ligand

The transport of a charged adsorbate biomolecule in a porous polymeric adsorbent medium and its adsorption onto the covalently immobilized ligands have been modeled and investigated using molecular dynamics modeling and simulations as the third part of a novel fundamental methodology developed for studying ion-exchange chromatography based bioseparations. To overcome computational challenges, a novel simulation approach is devised where appropriate atomistic and coarse grain models are employed simultaneously and the transport of the adsorbate is characterized through a number of locations representative of the progress of the transport process. The adsorbate biomolecule for the system studied in this work changes shape, orientation, and lateral position in order to proceed toward the site where adsorption occurs and exhibits decreased mass transport coefficients as it approaches closer to the immobilized ligand. Furthermore, because the ligands are surrounded by counterions carrying the same type of charge as the adsorbate biomolecule, it takes the biomolecule repeated attempts to approach toward a ligand in order to displace the counterions in the proximity of the ligand and to finally become adsorbed. The formed adsorbate-ligand complex interacts with the counterions and polymeric molecules and is found to evolve slowly and continuously from one-site (monovalent) interaction to multisite (multivalent) interactions. Such a transition of the nature of adsorption reduces the overall adsorption capacity of the ligands in the adsorbent medium and results in a type of surface exclusion effect. Also, the adsorption of the biomolecule also presents certain volume exclusion effects by not only directly reducing the pore volume and the availability of the ligands in the adjacent regions, but also causing the polymeric molecules to change to more compact structures that could further shield certain ligands from being accessible to subsequent adsorbate molecules. These findings have significant practical implications to the design and construction of polymeric porous adsorbent media for effective bioseparations and to the synthesis and operation of processes employed in the separation of biomolecules. The modeling and analysis methods presented in this work could also be suitable for the study of biocatalysis where an enzyme is immobilized on the surface of the pores of a porous medium.

Protein adsorption in porous adsorbent particles: A multiscale modeling study on inner radial humps in the concentration profiles of adsorbed protein induced by nonuniform ligand density distributions

Recently published results determined from molecular dynamics (MD) modeling and simulation studies have shown that the spatial distribution of the density of immobilized charged ligands in ion-exchange porous adsorbent particles is most likely nonuniform and the adsorbent particles also exhibit local nonelectroneutrality. In this work, the functional forms of the nonuniform spatial distributions of the density of the immobilized ligands in four different porous adsorbent media that were determined by MD studies were employed in a macroscopic continuum model describing the transport and adsorption of a single protein in the porous particles of the four different adsorbent media. The results clearly show that inner radial humps in the concentration profiles of the adsorbed protein can occur when the spatial distribution of the density of the immobilized ligands in the porous adsorbent particles is nonuniform and also has local maxima or minima along the radial direction in the particle. The results also indicate that the rate at which the equilibrium condition is approached depends significantly on the functional form of the spatial distribution of the density of the immobilized ligands. When adsorption equilibrium has been reached, the concentration profile of the adsorbed protein exhibits the shape of the spatial distribution of the density of the immobilized ligands. The results suggest that the technique of confocal scanning laser microscopy could be used to measure the concentration profile of an adsorbed protein at equilibrium and this measurement could provide the spatial distribution of the density of the immobilized ligands, and such measurements could also be used for quality control of the adsorbent medium. The results in this work have also implications in the modeling, design, analysis, and quality control of systems involving biocatalysis. Furthermore, the results clearly indicate that it is very important to study the dynamic behavior of an adsorption system having a nonuniform spatial distribution in the density of the immobilized charged ligands and where (i) both monovalent and multivalent interactions between the single charged adsorbate and the immobilized charged ligands occur and (ii) the values of the pH and ionic strength are such that the electrophoretic effects are active.

Effects on the dynamic utilization of the adsorptive capacity of chromatographic columns induced by non-uniform ligand density distributions.

The dynamic behavior of the breakthrough curves of a single adsorbate obtained from columns employing adsorbent media which differ from one another only on the spatial distribution of the immobilized ligands in the porous particles is examined. The spatial distributions of the immobilized ligands considered in this study are uniform and non-uniform, but the total number of immobilized ligands in the particles has the same value whether the spatial distribution is uniform or non-uniform. The results clearly show that the columns employing adsorbent particles in which the spatial distribution of the immobilized ligands is non-uniform and such that the concentration of the immobilized ligands increases monotonically from the center of the particle to the outer particle surface, exhibit (i) larger breakthrough times, (ii) steeper breakthrough curves, and (iii) higher dynamic utilization of the adsorptive capacity of the column as the superficial velocity of the flowing fluid stream in the column increases (throughput increase) than the columns using adsorbent particles in which the spatial distribution of the immobilized ligands is uniform. The importance of employing in the columns adsorbent media whose spatial ligand density distributions satisfy the mathematical property of monotonically increasing ligand concentration with increasing from the particle center radial position, will be significantly enhanced when (i) the size of the particle radius is increased, and (ii) continuous counter-current and periodic counter-current (simulated moving beds) operations are employed.

Modeling the construction of polymeric adsorbent media: Effects of counter-ions on ligand immobilization and pore structure

Molecular dynamics modeling and simulations are employed to study the effects of counter-ions on the dynamic spatial density distribution and total loading of immobilized ligands as well as on the pore structure of the resultant ion exchange chromatography adsorbent media. The results show that the porous adsorbent media formed by polymeric chain molecules involve transport mechanisms and steric resistances which cause the charged ligands and counter-ions not to follow stoichiometric distributions so that (i) a gradient in the local nonelectroneutrality occurs, (ii) non-uniform spatial density distributions of immobilized ligands and counter-ions are formed, and (iii) clouds of counter-ions outside the porous structure could be formed. The magnitude of these counter-ion effects depends on several characteristics associated with the size, structure, and valence of the counter-ions. Small spherical counter-ions with large valence encounter the least resistance to enter a porous structure and their effects result in the formation of small gradients in the local nonelectroneutrality, higher ligand loadings, and more uniform spatial density distributions of immobilized ligands, while the formation of exterior counter-ion clouds by these types of counter-ions is minimized. Counter-ions with lower valence charges, significantly larger sizes, and elongated shapes, encounter substantially greater steric resistances in entering a porous structure and lead to the formation of larger gradients in the local nonelectroneutrality, lower ligand loadings, and less uniform spatial density distributions of immobilized ligands, as well as substantial in size exterior counter-ion clouds. The effects of lower counter-ion valence on pore structure, local nonelectroneutrality, spatial ligand density distribution, and exterior counter-ion cloud formation are further enhanced by the increased size and structure of the counter-ion. Thus, the design, construction, and functionality of polymeric porous adsorbent media will significantly depend, for a given desirable ligand to be immobilized and represent the adsorption active sites, on the type of counter-ion that is used during the ligand immobilization process. Therefore, the molecular dynamics modeling and simulation approach presented in this work could contribute positively by representing an engineering science methodology to the design and construction of polymeric porous adsorbent media which could provide high intraparticle mass transfer and adsorption rates for the adsorbate biomolecules of interest which are desired to be separated by an adsorption process.

Molecular modeling of polymeric adsorbent media: The effects of counter‐ions on ligand immobilization and pore structure

Molecular dynamics modeling and simulations are employed to study the immobilization of ligands on the surface of the pores of a base porous polymeric matrix. The results show the significant effects that the counter-ions have on the spatial distribution of the density of immobilized ligands as well as on the pore size and pore connectivity distributions of the porous adsorbent medium being constructed. The results for the systems studied in this work indicate that by using doubly charged counter-ions whose numbers during ligand immobilization are half to those of singly charged counter-ions, the ligand immobilization process proceeds faster and the magnitude of local nonelectroneutrality becomes smaller. More importantly, the pore structures of the adsorbent media resulting from the system using doubly charged counter-ions have porous structures that are characterized by more mid-sized pores and higher pore connectivity than the porous adsorbent structures generated by the system employing singly charged counter-ions and, furthermore, the density distribution of the immobilized ligands in the porous structures where doubly charged counter-ions are employed tends to be more uniform laterally and the ligands are surrounded by fewer counter-ions. These characteristics affected by the use of doubly charged counter-ions could provide important advantages with respect to the transport and adsorption of adsorbate biomolecules of interest. Furthermore, the results of this work indicate that the type of counter-ions being used in the ligand immobilization process could represent an additional control variable for affecting the ligand density distribution as well as the pore size and pore connectivity distributions of the porous structure of the adsorbent medium being constructed.

Practical benefits for the separation process of drying that can be realized from novel multiscale modeling procedures that utilize the scientific information and results obtained from molecular dynamics modeling and simulation studies

ABSTRACT Novel multiscale modeling procedures are constructed and presented that use the scientific information and results determined from microscopic molecular dynamics (MD) modeling and simulation studies to calculate local effective values for the parameters that characterize the heat and mass transfer mechanisms of dynamic macroscopic continuum models (Euler physics of continua) that are used in practice to describe and predict the dynamic behavior of large scale in time and space (e.g., industrial scale), separation (e.g., drying; adsorption), and chemical and biochemical reaction engineering (e.g., chemical catalysis; biocatalysis; immobilized cell bioreactor systems) processes involving porous media whose pore structure is formed either by a solid rigid matter or by a solid soft matter. Furthermore, the results determined from MD modeling and simulation studies with regard to the energies of interaction between the molecules of the different species of the porous media during the time evolution (time varying) of the drying process can be used to design a time optimally controlled heat input system that could appropriately and accurately supply at any time during drying the amount of heat necessary to provide a desired drying rate with respect to both free and bound water and to satisfy the constraints that safeguard the quality properties of the product.