Multiscale Modeling of Mesoscale and Interfacial Phenomena

With rapidly emerging technologies that feature interfaces modified at the nanoscale, traditional macroscopic models are pushed to their limits to explain phenomena where molecular processes can play a key role. Often, such problems appear to defy explanation when treated with coarse-grained continuum models alone, yet remain prohibitively expensive from a molecular simulation perspective. A prominent example is surface nanobubbles: nanoscopic gaseous domains typically found on hydrophobic surfaces that have puzzled researchers for over two decades due to their unusually long lifetimes. We show how an entirely macroscopic, non-equilibrium model explains many of their anomalous properties, including their stability and abnormally small gas-side contact angles. From this purely transport perspective, we investigate how factors such as temperature and saturation affect nanobubbles, providing numerous experimentally testable predictions. However, recent work also emphasizes the relevance of molecular-scale phenomena that cannot be described in terms of bulk phases or pristine interfaces. This is true for nanobubbles as well, whose nanoscale heights may require molecular detail to capture the relevant physics, in particular near the bubble three-phase contact line.Therefore, there is a clear need for general ways to link molecular granularity and behavior with large-scale continuum models in the treatment of many interfacial problems. In light of this, we have developed a general set of simulation strategies that couple mesoscale particle-based continuum models to molecular regions simulated through conventional molecular dynamics (MD). In addition, we derived a transport model for binary mixtures that opens the possibility for a wide range of applications in biological and drug delivery problems, and is readily reconciled with our hybrid MD-continuum techniques. Approaches that couple multiple length scales for fluid mixtures are largely absent in the literature, and we provide a novel and general framework for multiscale modeling of systems featuring one or more dissolved species. This makes it possible to retain molecular detail for parts of the problem that require it while using a simple, continuum description for parts where high detail is unnecessary, reducing the number of degrees of freedom (i.e. number of particles) dramatically. This opens the possibility for modeling ion transport in biological processes and biomolecule assembly in ionic solution, as well as electrokinetic phenomena at interfaces such as corrosion. The number of particles in the system is further reduced through an integrated boundary approach, which we apply to colloidal suspensions. In this thesis, we describe this general framework for multiscale modeling single- and multicomponent systems, provide several simple equilibrium and non-equilibrium case studies, and discuss future applications

Discrete charge patterns, Coulomb correlations and interactions in protein solutions

The effective Coulomb interaction between globular proteins is calculated as a function of monovalent salt concentration $c_s$, by explicit Molecular Dynamics simulations of pairs of model proteins in the presence of microscopic co and counterions. For discrete charge patterns of monovalent sites on the surface, the resulting osmotic virial coefficient $B_2$ is found to be a strikingly non-monotonic function of $c_s$. The non-monotonicity follows from a subtle Coulomb correlation effect which is completely missed by conventional non-linear Poisson-Boltzmann theory and explains various experimental findings.Comment: 4 twocolumn pages with 4 figure

The targeted delivery of multicomponent cargos to cancer cells by nanoporous particle-supported lipid bilayers.

Encapsulation of drugs within nanocarriers that selectively target malignant cells promises to mitigate side effects of conventional chemotherapy and to enable delivery of the unique drug combinations needed for personalized medicine. To realize this potential, however, targeted nanocarriers must simultaneously overcome multiple challenges, including specificity, stability and a high capacity for disparate cargos. Here we report porous nanoparticle-supported lipid bilayers (protocells) that synergistically combine properties of liposomes and nanoporous particles. Protocells modified with a targeting peptide that binds to human hepatocellular carcinoma exhibit a 10,000-fold greater affinity for human hepatocellular carcinoma than for hepatocytes, endothelial cells or immune cells. Furthermore, protocells can be loaded with combinations of therapeutic (drugs, small interfering RNA and toxins) and diagnostic (quantum dots) agents and modified to promote endosomal escape and nuclear accumulation of selected cargos. The enormous capacity of the high-surface-area nanoporous core combined with the enhanced targeting efficacy enabled by the fluid supported lipid bilayer enable a single protocell loaded with a drug cocktail to kill a drug-resistant human hepatocellular carcinoma cell, representing a 10(6)-fold improvement over comparable liposomes

Shape-Dependent Thermodynamics and Non-Local Hydrodynamics in a Non-Gibbsian Steady-State of a Drift-Diffusion System

Shape-dependent thermodynamics and non-local hydrodynamics are argued to occur in dissipative steady states of driven diffusive systems. These predictions are confirmed by numerical simulations. Unlike power-law correlations, these phenomena cannot be explained by a hypothesis of ``criticality''. Instead, they require the effective Hamiltonian of the system to contain very long-range potentials, making the invariant probability measures formally ``non-Gibbsian''.Comment: 4 pages, Latex Version 2.09, 1 Postscript figur

Non-monotonic variation with salt concentration of the second virial coefficient in protein solutions

The osmotic virial coefficient $B_2$ of globular protein solutions is calculated as a function of added salt concentration at fixed pH by computer simulations of the ``primitive model''. The salt and counter-ions as well as a discrete charge pattern on the protein surface are explicitly incorporated. For parameters roughly corresponding to lysozyme, we find that $B_2$ first decreases with added salt concentration up to a threshold concentration, then increases to a maximum, and then decreases again upon further raising the ionic strength. Our studies demonstrate that the existence of a discrete charge pattern on the protein surface profoundly influences the effective interactions and that non-linear Poisson Boltzmann and Derjaguin-Landau-Verwey-Overbeek (DLVO) theory fail for large ionic strength. The observed non-monotonicity of $B_2$ is compared to experiments. Implications for protein crystallization are discussed.Comment: 43 pages, including 17 figure

Adsorption of mono- and multivalent cat- and anions on DNA molecules

Adsorption of monovalent and multivalent cat- and anions on a deoxyribose nucleic acid (DNA) molecule from a salt solution is investigated by computer simulation. The ions are modelled as charged hard spheres, the DNA molecule as a point charge pattern following the double-helical phosphate strands. The geometrical shape of the DNA molecules is modelled on different levels ranging from a simple cylindrical shape to structured models which include the major and minor grooves between the phosphate strands. The densities of the ions adsorbed on the phosphate strands, in the major and in the minor grooves are calculated. First, we find that the adsorption pattern on the DNA surface depends strongly on its geometrical shape: counterions adsorb preferentially along the phosphate strands for a cylindrical model shape, but in the minor groove for a geometrically structured model. Second, we find that an addition of monovalent salt ions results in an increase of the charge density in the minor groove while the total charge density of ions adsorbed in the major groove stays unchanged. The adsorbed ion densities are highly structured along the minor groove while they are almost smeared along the major groove. Furthermore, for a fixed amount of added salt, the major groove cationic charge is independent on the counterion valency. For increasing salt concentration the major groove is neutralized while the total charge adsorbed in the minor groove is constant. DNA overcharging is detected for multivalent salt. Simulations for a larger ion radii, which mimic the effect of the ion hydration, indicate an increased adsorbtion of cations in the major groove.Comment: 34 pages with 14 figure

Calibrative approaches to protein solubility modeling of a mutant series using physicochemical descriptors

A set of physicochemical properties describing a protein of known structure is employed for a calibrative approach to protein solubility. Common hydrodynamic and electrophoretic properties routinely measured in the bio-analytical laboratory such as zeta potential, dipole moment, the second osmotic virial coefficient are first estimated in silico as a function a pH and solution ionic strength starting with the protein crystal structure. The utility of these descriptors in understanding the solubility of a series of ribonuclease Sa mutants is investigated. A simple two parameter model was trained using solubility data of the wild type protein measured at a restricted number of solution pHs. Solubility estimates of the mutants demonstrate that zeta potential and dipole moment may be used to rationalize solubility trends over a wide pH range. Additionally a calibrative model based on the protein’s second osmotic virial coefficient, B22 was developed. A modified DVLO type potential along with a simplified representation of the protein allowed for efficient computation of the second viral coefficient. The standard error of prediction for both models was on the order of 0.3 log S units. These results are very encouraging and demonstrate that these models may be trained with a small number of samples and employed extrapolatively for estimating mutant solubilities

An Algorithm for Emulsion Stability Simulations: Account of Flocculation, Coalescence, Surfactant Adsorption and the Process of Ostwald Ripening

The first algorithm for Emulsion Stability Simulations (ESS) was presented at the V Conferencia Iberoamericana sobre Equilibrio de Fases y Diseño de Procesos [Luis, J.; García-Sucre, M.; Urbina-Villalba, G. Brownian Dynamics Simulation of Emulsion Stability In: Equifase 99. Libro de Actas, 1st Ed., Tojo J., Arce, A., Eds.; Solucion’s: Vigo, Spain, 1999; Volume 2, pp. 364–369]. The former version of the program consisted on a minor modification of the Brownian Dynamics algorithm to account for the coalescence of drops. The present version of the program contains elaborate routines for time-dependent surfactant adsorption, average diffusion constants, and Ostwald ripening