BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Y. C. Zhou

Y. C. Zhou contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAQuantitative Methodscond-mat.softBiomoleculesmath.APphysics.comp-phBiological PhysicsNumerical Analysis

Trust snapshot

Quick read

Trust 21 - Emerging
10works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2020arXiv

A Computational Model of Protein Induced Membrane Morphology with Geodesic Curvature Driven Protein-Membrane Interface

Continuum or hybrid modeling of bilayer membrane morphological dynamics induced by embedded proteins necessitates the identification of protein-membrane interfaces and coupling of deformations of two surfaces. In this article we developed (i) a minimal total geodesic curvature model to describe these interfaces, and (ii) a numerical one-one mapping between two surface through a conformal mapping of each surface to the common middle annulus. Our work provides the first computational tractable approach for determining the interfaces between bilayer and embedded proteins. The one-one mapping allows a convenient coupling of the morphology of two surfaces. We integrated these two new developments into the energetic model of protein-membrane interactions, and developed the full set of numerical methods for the coupled system. Numerical examples are presented to demonstrate (1) the efficiency and robustness of our methods in locating the curves with minimal total geodesic curvature on highly complicated protein surfaces, (2) the usefulness of these interfaces as interior boundaries for membrane deformation, and (3) the rich morphology of bilayer surfaces for different protein-membrane interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Enriched Gradient Recovery for Interface Solutions of the Poisson-Boltzmann Equation

Accurate calculation of electrostatic potential and gradient on the molecular surface is highly desirable for the continuum and hybrid modeling of large scale deformation of biomolecules in solvent. In this article a new numerical method is proposed to calculate these quantities on the dielectric interface from the numerical solutions of the Poisson-Boltzmann equation. Our method reconstructs a potential field locally in the least square sense on the polynomial basis enriched with Green's functions, the latter characterize the Coulomb potential induced by charges near the position of reconstruction. This enrichment resembles the decomposition of electrostatic potential into singular Coulomb component and the regular reaction field in the Generalized Born methods. Numerical experiments demonstrate that the enrichment recovery produces drastically more accurate and stable potential gradients on molecular surfaces compared to classical recovery techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Curvature-driven molecular flows on membrane surfaces

Morphological change of bilayer membrane in vivo is not a spontaneous procedure but modulated by various types of proteins in general. Most of these modulations are associated with the localization of related proteins in the crowded lipid environment in bilayer membrane. This work presents an mathematical model for the localization of multiple species of diffusion molecules on membrane surfaces. We start with the energetic description of the distributions of molecules on curved membrane surface, by assembling the bending energy of bilayer membrane and the entropic energy of diffusive molecules. We introduce the spontaneous curvature of molecules in membrane, and define the spontaneous curvature of bilayer membrane as a function of the molecule concentrations on membrane surfaces. This connection gives rise to a drift-diffusion equation to govern the gradient flows of the surface molecule concentrations. We recast the energetic formulation and the related governing equations in the Eulerian framework by using a phase field function that defines the membrane morphology. Computational simulations with the proposed mathematical model and related numerical techniques predict the molecular localization on membrane surfaces at locations with preferred mean curvature.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Interface Formulation and High Order Numerical Solutions of PDEs with Low Regularity

Linear elastic fracture mechanics admit analytic solutions that have low regularity at crack tips. Current numerical methods for partial differential equations (PDEs) of this type suffer from the constraint of such low regularity, and fail to deliver optimal high order rate of convergence. We approach the problem by (i) choosing an artificial interface to enclose the center of the low regularity; and (ii) representing the solution in the interior of artificial interface as unknown linear combination of known modes of low regular solutions. This gives rise to an interface formulation of the original PDE, and the linear combination are represented the interface conditions. By enforcing the smooth component of numerical solution in the interior domain to be approximately zero, a least square problem is obtained for the unknown coefficients. The solution of this least square problem will provide approximate interface conditions for the numerical solution of the PDE in the exterior domain. The potential of our interface formulation is favorably demonstrated by numerical experiments on 1-D and 2-D Poisson equations with low regular solutions. High order numerical solutions of unknown coefficients and PDEs are obtained. This proves the potential of the proposed interface formulation as the theoretical basis for solving linear elastic fracture mechanics problems. We indicate the relations between our interface formulation and domain decomposition methods as well as a regularization strategy for the Poisson-Boltzmann equation with singular charge density.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Variational Methods for Biomolecular Modeling

Structure, function and dynamics of many biomolecular systems can be characterized by the energetic variational principle and the corresponding systems of partial differential equations (PDEs). This principle allows us to focus on the identification of essential energetic components, the optimal parametrization of energies, and the efficient computational implementation of energy variation or minimization. Given the fact that complex biomolecular systems are structurally non-uniform and their interactions occur through contact interfaces, their free energies are associated with various interfaces as well, such as solute-solvent interface, molecular binding interface, lipid domain interface, and membrane surfaces. This fact motivates the inclusion of interface geometry, particular its curvatures, to the parametrization of free energies. Applications of such interface geometry based energetic variational principles are illustrated through three concrete topics: the multiscale modeling of biomolecular electrostatics and solvation that includes the curvature energy of the molecular surface, the formation of microdomains on lipid membrane due to the geometric and molecular mechanics at the lipid interface, and the mean curvature driven protein localization on membrane surfaces. By further implicitly representing the interface using a phase field function over the entire domain, one can simulate the dynamics of the interface and the corresponding energy variation by evolving the phase field function, achieving significant reduction of the number of degrees of freedom and computational complexity. Strategies for improving the efficiency of computational implementations and for extending applications to coarse-graining or multiscale molecular simulations are outlined.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fast Algorithm for Simulating Lipid Vesicle Deformation I: Spherical Harmonic Approximation

Lipid vesicles appear ubiquitously in biological systems. Understanding how the mechanical and intermolecular interations deform vesicle membrane is a fundamental question in biophysics. In this article we developed a fast algorithm to compute the surface configurations of lipid vesicles by introducing the surface harmonic functions to approximate the surfaces. This parameterization of the surfaces allows an analytical computation of the membrane curvature energy and its gradient for the efficient minimization of the curvature energy using a nonlinear conjugate gradient method. Our approach drastically reduces the degrees of freedom for approximating the membrane surfaces compared to the previously developed finite element and finite difference methods. Vesicle deformations with a reduced volume larger than 0.65 can be well approximated by using as small as 49 surface harmonic functions. The method thus has a great potential to reduce the computational expense of tracking multiple vesicles which deform for their interaction with external fields.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electrostatic forces on charged surfaces of bilayer lipid membranes

Simulating protein-membrane interactions is an important and dynamic area of research. A proper definition of electrostatic forces on membrane surfaces is necessary for developing electromechanical models of protein-membrane interactions. Here we modeled the bilayer membrane as a continuum with general continuous distributions of lipids charges on membrane surfaces. A new electrostatic potential energy functional was then defined for this solvated protein-membrane system. We investigated the geometrical transformation properties of the membrane surfaces under a smooth velocity field. These properties allows us to apply the Hadamard-Zolesio structure theorem, and the electrostatic forces on membrane surfaces can be computed as the shape derivative of the electrostatic energy functional.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Genetic Exponentially Fitted Method for Solving Multi-dimensional Drift-diffusion Equations

A general approach was proposed in this article to develop high-order exponentially fitted basis functions for finite element approximations of multi-dimensional drift-diffusion equations for modeling biomolecular electrodiffusion processes. Such methods are highly desirable for achieving numerical stability and efficiency. We found that by utilizing the one-one correspondence between continuous piecewise polynomial space of degree $k+1$ and the divergence-free vector space of degree $k$, one can construct high-order 2-D exponentially fitted basis functions that are strictly interpolative at a selected node set but are discontinuous on edges in general, spanning nonconforming finite element spaces. First order convergence was proved for the methods constructed from divergence-free Raviart-Thomas space $RT_0^0$ at two different node sets

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Electrodiffusion of Lipids on Membrane Surfaces

Random lateral translocation of lipids and proteins is a universal process on membrane surfaces. Local aggregation or organization of lipids and proteins can be induced when this lateral random diffusion is mediated by the electrostatic interactions and membrane curvature. Though the lateral diffusion rates of lipids on membrane of various compositions are measured and the electrostatic free energies of predetermined protein-membrane-lipid systems can be computed, the process of the aggregation and the evolution to the electrostatically favorable states remain undetermined. Here we propose an electrodiffusion model, based on the variational principle of free energy functional, for the self-consistent lateral drift-diffusion of multiple species of charged lipids on membrane surfaces. Finite sizes of lipids are modeled to enforce the geometrical constraint of the lipid concentration on membrane surfaces. A surface finite element method is developed to appropriate the Laplace-Beltrami operators in the partial differential equations (PDEs) of the model. Our model properly describes the saturation of lipids on membrane surface, and correctly predicts that the MARCKS peptide can consistently sequester three multivalent phosphatidylinositol 4,5-bisphosphate (\pip2) lipids through its basic amino acid residues, even there is a large fraction of monovalent phosphatidylserine (PS) in the membrane. Solutions of the PDEs also show the salt-dependence of the lipid sequestration.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Continuum Electromechanical Modeling of Protein-Membrane Interaction

A continuum electromechanical model is proposed to describe the membrane curvature induced by electrostatic interactions in a solvated protein-membrane system. The model couples the macroscopic strain energy of membrane and the electrostatic solvation energy of the system, and equilibrium membrane deformation is obtained by minimizing the electro-elastic energy functional with respect to the dielectric interface. The model is illustrated with the systems with increasing geometry complexity and captures the sensitivity of membrane curvature to the permanent and mobile charge distributions.

0 citations0 reviewsNo code linkedOpen access