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.

Use of dimensionless residence time to study variations in breakthrough behaviour in expanded beds formed from varied particle size distributions.

This work demonstrates an experimental method for studying breakthrough behaviour in expanded beds. The behaviour of beds made with differently sized particles were studied at varying flowrates. The use of a dimensionless residence time measurement allowed a more valid comparison of breakthrough characteristics in expanded bed operation by compensating for the changes in bed volume that occur during expansion. We demonstrate that bed breakthrough behaviour can be compared directly even when the beds contain different-sized particles and hence have different expanded volumes. By utilising this concept we demonstrate that, in the case of the Alcohol Dehydrogenase (ADH) / STREAMLINE Phenyl system used here, there was little or no variation in ADH breakthrough behaviour between beds of differently sized particles operating at flowrates above 100 cm/h. This suggests that the higher specific surface area and hence binding capacity of smaller particles is negated in this case due to mass transfer limitations and the increase in system void volume even at normal operating flowrates of 200-300 cm/h.

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.

Modeling and analysis of the dynamic behavior of mechanisms that result in the development of inner radial humps in the concentration of a single adsorbate in the adsorbed phase of porous adsorbent particles observed in confocal scanning laser microscopy experiments: diffusional mass transfer and adsorption in the presence of an electrical double layer.

A theoretical model for adsorption of a single charged adsorbate that accounts for the presence of an electrical double layer in the pores of adsorbent particles is constructed and solved. The dynamic behavior of the mechanisms of the model can result in the development of inner radial humps (concentration rings) in the concentration of a single charged analyte (adsorbate) in the adsorbed phase of porous adsorbent particles. The results of the present work demonstrate the implication of the concept regarding the effect of the presence of an electrical double layer in the pores of adsorbent particles and the induced interactions between the electrostatic potential distribution and the mechanisms of mass transport of the species by diffusion, electrophoretic migration, and adsorption. Furthermore, the mechanisms of the model could explain qualitatively the development of the concentration ring (hump) observed in confocal scanning laser microscopy experiments.

A method for the decrease of phenolic content in commercial canola meal using an enzyme preparation secreted by the white-rot fungus Trametes versicolor.

An enzymatic process for upgrading the quality of canola meal (CM) by decreasing its phenolic content was investigated. The new method was based on the addition of the enzyme preparation from white-rot fungus Trametes versicolor to the meal-buffer slurry. A 98% decrease in the concentration of SAE was observed after 1 h of the treatment. The following process variables were considered for optimizing the process: pH, temperature, enzyme, meal, and oxygen concentrations. It was found that: (1) the natural buffering capacity of CM resulted in a negligible effect of the pH of the buffer, which was used as the continuous phase in the process, on the extent of decrease in sinapic acid esters (SAE); (2) the system was saturated with the enzyme when its concentration was 4 nkat/mL of the continuous phase; and (3) the optimum temperature was 50 degrees C. The process could be carried out even at higher temperatures due to the protective action of CM, which resulted in an increase in the thermal stability of the enzyme. The particle size influenced the extraction of the SAE from the meal, indicating that, at lower SAE concentrations, the process became diffusion limited. This result, together with those showing no effect of the intensity of agitation, indicated that the enzymatic process can be characterized by high Biot numbers. During the enzymatic process, the molar concentration of available oxygen can become a limiting factor when it is more than four times lower than the molar concentration of phenolics in the treated meal. The new enzymatic method was compared with other methods reported in the literature for the decrease in the phenolic content of rapeseed meals. It was found that, among the methods tested, the enzymatic treatment was the most effective, followed by the lime treatment. The enzymatic process did not reduce the quality of the protein isolates prepared from the CM. After the addition of a simple acetone-washing step, the isolate from the enzymatically treated meal had even better properties. Copyright 1999 John Wiley & Sons, Inc.

Transformation of 3,5-dimethoxy,4-hydroxy cinnamic acid by polyphenol oxidase from the fungus Trametes versicolor: product elucidation studies.

Sinapic acid (SA), 3,5-dimethoxy,4-hydroxy cinnamic acid, was incubated with a crude polyphenol oxidase from the fungus Trametes versicolor. Some products of this transformation were isolated and their structures identified using mass spectrometry, nuclear magnetic resonance and Fourier transform infrared spectroscopy, and X-ray crystallography. It was found that the enzymatic oxidation of SA includes two distinct phases. In the initial phase SA is enzymatically transformed to r-1H-2c,6c-bis-(4'-hydroxy-3', 5'-dimethoxyphenyl)-3,7-dioxabicyclo-[3,3,0]-octane-4,8-dione, dehydrodisinapic acid dilactone. The mechanism of this reaction may involve coupling of two phenoxy radicals by the beta-beta mode and subsequent intramolecular nucleophilic attack. In the second phase dehydrodisinapic acid dilactone is transformed by polyphenol oxidase into several intermediate products, including 4-(4-(3, 5-dimethoxy-4-oxo-2,5-cyclohexadienyliden)-1, 4-dihydroxy-(E)-2-butenylidene)-2,6-dimethoxy-2, 5-cyclohexadien-1-one. The final product of the overall transformation of SA is 2,6-dimethoxy-p-benzoquinone. The obtained results were used to propose a part of the transformation pathway for the enzymatic oxidation of SA by polyphenol oxidase.

Modeling the enzymatic transformation of 3,5-dimethoxy,4-hydroxy cinnamic acid by polyphenoloxidase from the white-rot fungus Trametes versicolor.

A mechanism for transforming sinapic acid by a polyphenoloxidase from Trametes versicolor was investigated using changes in sinapic acid and oxygen concentrations during the reaction. The experiments were performed in a closed system without supplemental oxygen. The effects of temperature and initial oxygen concentration on the reaction rates were examined. To compare the obtained results with those from spectrophotometric studies, some runs were performed using an open system with supplemental oxygen. Sinapic acid transformation can be described by the Theorell-Chance Bi-Bi or Ordered Bi-Bi mechanisms. This reacting system consisted also of additional enzymatic reactions between the products of sinapic acid transformation and oxygen. A mathematical model was developed using four ordinary differential equations that represent the Theorell-Chance Bi-Bi mechanism with three alternate substrates. Model parameters (i.e., rate constants) were determined using the data collected at three different temperatures. On the basis of the transition state theory, relationships between these constants and temperature were established. It is shown that, in the open system, the observed change in the enzyme activity at higher temperatures was caused by two opposing phenomena: an Arrhenius effect which increased the rate, and a solubility effect which reduced the rate due to a lower oxygen concentration. This finding allows us to recommend better conditions for spectrophotometric methods, the assay most commonly used to evaluate this and similar enzymes.

Comparison of 3 methods for the determination of sinapic acid ester content in enzymatically treated canola meals.

The enzymatic reduction of sinapic acid ester content in canola meal using polyphenol oxidase from the fungus T. versicolor was investigated. To determine the effectiveness of this new process, the results obtained using two spectrophotometric methods and an HPLC analytical method for assaying sinapic acid ester content in the treated and untreated meals were compared. It was found that all the methods gave practically the same results when the samples from untreated canola meals were analysed. However, both of the spectrophotometric methods overestimated the sinapic acid ester content in the enzymatically treated meal by 7%-20%, as compared to the results obtained using HPLC. It was found that the sensitivity limits for the spectrophotometric methods used for the determination of sinapic acid ester content in enzymatically treated canola meals were 2.67 g and 1.47 g phenolics/kg meal for the direct and chemical spectrophotometric methods respectively. A correlation between the results obtained using the spectrophotometric and HPLC methods is given. The enzymatic treatment resulted in a negligible amount of phenolics in the treated meal.