skip to main content

Title: Hybrid PDE‐kMC modeling approach to simulate multivalent lectin‐glycan binding process

Glycans are the major components of the cellular membranes and mediate many cellular processes via their interactions with lectins. A kinetic Monte Carlo (kMC) model was proposed previously to incorporate the key features of glycan‐lectin interactions such as multivalency and glycan diffusion, and its accuracy has been validated by experiments. However, computational cost of the kMC model is its major bottleneck. In this study, a hybrid model combining a partial differential equation (PDE) with the kMC model is proposed to greatly reduce the computational cost while preserving the accuracy. Specifically, glycan diffusion is simulated by the PDE for improving computational efficiency since the glycan diffusion execution through the kMC is computationally expensive. The hybrid PDE‐kMC model is employed to simulate the binding dynamics between cholera toxin subunit B and gangliosides on cellular membranes. The accuracy and efficiency of the proposed model was demonstrated by comparing with the sole kMC model.

more » « less
Award ID(s):
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
AIChE Journal
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model’s predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.

    more » « less
  2. Gilbert, Jack (Ed.)
    ABSTRACT The cyanobacterium Trichodesmium is an important contributor of new nitrogen (N) to the surface ocean, but its strategies for protecting the nitrogenase enzyme from inhibition by oxygen (O 2 ) remain poorly understood. We present a dynamic physiological model to evaluate hypothesized conditions that would allow Trichodesmium to carry out its two conflicting metabolic processes of N 2 fixation and photosynthesis. First, the model indicates that managing cellular O 2 to permit N 2 fixation requires high rates of respiratory O 2 consumption. The energetic cost amounts to ∼80% of daily C fixation, comparable to the observed diminution of the growth rate of Trichodesmium relative to other phytoplankton. Second, by forming a trichome of connected cells, Trichodesmium can segregate N 2 fixation from photosynthesis. The transfer of stored C to N-fixing cells fuels the respiratory O 2 consumption that protects nitrogenase, while the reciprocal transfer of newly fixed N to C-fixing cells supports cellular growth. Third, despite Trichodesmium lacking the structural barrier found in heterocystous species, the model predicts low diffusivity of cell membranes, a function that may be explained by the presence of Gram-negative membrane, production of extracellular polysaccharide substances (EPS), and “buffer cells” that intervene between N 2 -fixing and photosynthetic cells. Our results suggest that all three factors—respiratory protection, trichome formation, and diffusion barriers—represent essential strategies that, despite their energetic costs, facilitate the growth of Trichodesmium in the oligotrophic aerobic ocean and permit it to be a major source of new reactive nitrogen. IMPORTANCE Trichodesmium is a major nitrogen-fixing cyanobacterium and exerts a significant influence on the oceanic nitrogen cycle. It is also a widely used model organism in laboratory studies. Since the nitrogen-fixing enzyme nitrogenase is extremely sensitive to oxygen, how these surface-dwelling plankton manage the two conflicting processes of nitrogen fixation and photosynthesis has been a long-standing question. In this study, we developed a simple model of metabolic fluxes of Trichodesmium capturing observed daily cycles of photosynthesis, nitrogen fixation, and boundary layer oxygen concentrations. The model suggests that forming a chain of cells for spatially segregating nitrogen fixation and photosynthesis is essential but not sufficient. It also requires a barrier against oxygen diffusion and high rates of oxygen scavenging by respiration. Finally, the model indicates an extremely short life span of oxygen-enabling cells to instantly create a low-oxygen environment upon deactivation of photosynthesis. 
    more » « less
  3. Summary

    Quasicontinuum (QC) is a partitioned‐domain multiscale method where the full‐atomistic description is retained in highly deformable subdomains while the energy in the remaining region is approximated using a continuum constitutive relation called the Cauchy‐Born (CB) rule. While the CB rule can capture the full nonlinearity of the strain energy, its computational cost is one or two orders of magnitude higher than the linear elasticity (LE) constitutive relation. In this paper, we propose an extension to the QC method, referred to as “hybrid QC,” where both the CB and LE relations are combined to increase the efficiency of the original QC method. This extension is summed up to a two‐step process. First, the continuum region is further divided into subregions with CB being restricted to subregions adjoining the atomistic region while LE is adopted in the remaining subregions. Second, a corrective scheme is applied to the LE subregion to recover the higher‐order accuracy. The hybrid QC method is validated through two examples: (i) a crystal containing a dipole of Lomer dislocations and (ii) a nanoindentation system. The simulation results show that the proposed hybrid QC method is numerically more efficient than the original QC method while maintaining virtually the same accuracy.

    more » « less
  4. Lateral stiffness of structural components, such as reinforced concrete (RC) columns, plays an important role in resisting the lateral earthquake loads. The lateral stiffness relates the lateral force to the lateral deformation, having a critical effect on the accuracy of the lateral seismic response predictions. The classical methods (e.g. fiber beam–column model) to estimate the lateral stiffness require calculations from section, element, and structural levels, which is time-consuming. Moreover, the shear deformation and bond-slip effect may also need to be included to more accurately calculate the lateral stiffness, which further increases the modeling difficulties and the computational cost. To reduce the computational time and enhance the accuracy of the predictions, this article proposes a novel data-driven method to predict the laterally seismic response based on the estimated lateral stiffness. The proposed method integrates the machine learning (ML) approach with the hysteretic model, where ML is used to compute the parameters that govern the nonlinear properties of the lateral response of target structural components directly from a training set composed of experimental data (i.e. data-driven procedure) and the hysteretic model is used to directly output the lateral stiffness based on the computed parameters and then to perform the seismic analysis. We apply the proposed method to predict the lateral seismic response of various types of RC columns subjected to cyclic loading and ground motions. We present the detailed model formulation for the application, including the developments of a modified hysteretic model, a hybrid optimization algorithm, and two data-driven seismic response solvers. The results predicted by the proposed method are compared with those obtained by classical methods with the experimental data serving as the ground truth, showing that the proposed method significantly outperforms the classical methods in both generalized prediction capabilities and computational efficiency.

    more » « less
  5. This article addresses the operational optimization of industrial steam systems under device efficiency uncertainty using a data‐driven adaptive robust optimization approach. A semiempirical model of steam turbine is first developed based on process mechanism and operational data. Uncertain parameters of the proposed steam turbine model are further derived from the historical process data. A robust kernel density estimation method is then used to construct the uncertainty sets for modeling these uncertain parameters. The data‐driven uncertainty sets are incorporated into a two‐stage adaptive robust mixed‐integer linear programming (MILP) framework for operational optimization of steam systems to minimize the total operating cost. Integer variables are introduced to model the on/off decisions of the steam turbines and electrical motors, which are the major energy consumers of the steam system. By applying the affine decision rule, the proposed multilevel optimization model is transformed into its robust counterpart, which is a single‐level MILP problem. The proposed framework is applied to the steam system of a real‐world ethylene plant to demonstrate its applicability. © 2018 American Institute of Chemical EngineersAIChE J, 65: e16500 2019

    more » « less