Granular packings with sliding, rolling, and twisting friction.

Intuition tells us that a rolling or spinning sphere will eventually stop due to the presence of friction and other dissipative interactions. The resistance to rolling and spinning or twisting torque that stops a sphere also changes the microstructure of a granular packing of frictional spheres by increasing the number of constraints on the degrees of freedom of motion. We perform discrete element modeling simulations to construct sphere packings implementing a range of frictional constraints under a pressure-controlled protocol. Mechanically stable packings are achievable at volume fractions and average coordination numbers as low as 0.53 and 2.5, respectively, when the particles experience high resistance to sliding, rolling, and twisting. Only when the particle model includes rolling and twisting friction were experimental volume fractions reproduced.

Heterogeneous partition of cellular blood-borne nanoparticles through microvascular bifurcations.

Blood flowing through microvascular bifurcations has been an active research topic for many decades, while the partitioning pattern of nanoscale solutes in the blood remains relatively unexplored. Here we demonstrate a multiscale computational framework for direct numerical simulation of the nanoparticle (NP) partitioning through physiologically relevant vascular bifurcations in the presence of red blood cells (RBCs). The computational framework is established by embedding a particulate suspension inflow-outflow boundary condition into a multiscale blood flow solver. The computational framework is verified by recovering a tubular blood flow without a bifurcation and validated against the experimental measurement of an intravital bifurcation flow. The classic Zweifach-Fung (ZF) effect is shown to be well captured by the method. Moreover, we observe that NPs exhibit a ZF-like heterogeneous partition in response to the heterogeneous partition of the RBC phase. The NP partitioning prioritizes the high-flow-rate daughter branch except for extreme (large or small) suspension flow partition ratios under which the complete phase separation tends to occur. By analyzing the flow field and the particle trajectories, we show that the ZF-like heterogeneity in the NP partition can be explained by the RBC-entrainment effect caused by the deviation of the flow separatrix preceded by the tank treading of RBCs near the bifurcation junction. The recovery of homogeneity in the NP partition under extreme flow partition ratios is due to the plasma skimming of NPs in the cell-free layer. These findings, based on the multiscale computational framework, provide biophysical insights to the heterogeneous distribution of NPs in microvascular beds that are observed pathophysiologically.

Controlling Binder Adhesion to Impact Electrode Mesostructures and Transport.

The complex three-phase composition of lithium-ion battery electrodes, containing an ion-conducting pore phase, a nanoporous electron-conducting carbon binder domain (CBD) phase, and an active material (AM) phase, provides several avenues of mesostructural engineering to enhance battery performance. We demonstrate a promising strategy for engineering electrode mesostructures by controlling the strength of adhesion between the AM and CBD phases. Using high-fidelity, physics-based colloidal and granular dynamics simulations, we predict that this strategy can provide significant control over electrochemical transport-relevant properties such as ionic conductivity, electronic conductivity, and available AM-electrolyte interface area. Importantly, the proposed strategy could be experimentally realized through surface functionalization of the AM and CBD phases and would be compatible with traditional electrode manufacturing methods.

Random walks on jammed networks: Spectral properties.

Using random walk analyses we explore diffusive transport on networks obtained from contacts between isotropically compressed, monodisperse, frictionless sphere packings generated over a range of pressures in the vicinity of the jamming transition p→0. For conductive particles in an insulating medium, conduction is determined by the particle contact network with nodes representing particle centers and edges contacts between particles. The transition rate is not homogeneous, but is distributed inhomogeneously due to the randomness of packing and concomitant disorder of the contact network, e.g., the distribution of the coordination number. A narrow escape time scale is used to write a Markov process for random walks on the particle contact network. This stochastic process is analyzed in terms of spectral density of the random, sparse, Euclidean and real, symmetric, positive, semidefinite transition rate matrix. Results show network structures derived from jammed particles have properties similar to ordered, euclidean lattices but also some unique properties that distinguish them from other structures that are in some sense more homogeneous. In particular, the distribution of eigenvalues of the transition rate matrix follow a power law with spectral dimension 3. However, quantitative details of the statistics of the eigenvectors show subtle differences with homogeneous lattices and allow us to distinguish between topological and geometric sources of disorder in the network.

Diffusion in Jammed Particle Packs.

Using random walk simulations we explore diffusive transport through monodisperse sphere packings over a range of packing fractions ϕ in the vicinity of the jamming transition at ϕ(c). Various diffusion properties are computed over several orders of magnitude in both time and packing pressure. Two well-separated regimes of normal "Fickian" diffusion, where the mean squared displacement is linear in time, are observed. The first corresponds to diffusion inside individual spheres, while the latter is the long-time bulk diffusion. The intermediate anomalous diffusion regime and the long-time value of the diffusion coefficient are both shown to be controlled by particle contacts, which in turn depend on proximity to ϕ(c). The time required to recover normal diffusion t* scales as (ϕ-ϕ(c))(-0.5) and the long-time diffusivity D(∞)∼(ϕ-ϕ(c))0.5, or D(∞)∼1/t*. It is shown that the distribution of mean first passage times associated with the escape of random walkers between neighboring particles controls both t* and D(∞) in the limit ϕ→ϕ(c).

Ligand structure and mechanical properties of single-nanoparticle-thick membranes.

The high mechanical stiffness of single-nanoparticle-thick membranes is believed to result from the local structure of ligand coatings that mediate interactions between nanoparticles. These ligand structures are not directly observable experimentally. We use molecular dynamics simulations to observe variations in ligand structure and simultaneously measure variations in membrane mechanical properties. We have shown previously that ligand end group has a large impact on ligand structure and membrane mechanical properties. Here we introduce and apply quantitative molecular structure measures to these membranes and extend analysis to multiple nanoparticle core sizes and ligand lengths. Simulations of nanoparticle membranes with a nanoparticle core diameter of 4 or 6 nm, a ligand length of 11 or 17 methylenes, and either carboxyl (COOH) or methyl (CH(3)) ligand end groups are presented. In carboxyl-terminated ligand systems, structure and interactions are dominated by an end-to-end orientation of ligands. In methyl-terminated ligand systems large ordered ligand structures form, but nanoparticle interactions are dominated by disordered, partially interdigitated ligands. Core size and ligand length also affect both ligand arrangement within the membrane and the membrane's macroscopic mechanical response, but are secondary to the role of the ligand end group. Moreover, the particular end group (COOH or CH(3)) alters the nature of how ligand length, in turn, affects the membrane properties. The effect of core size does not depend on the ligand end group, with larger cores always leading to stiffer membranes. Asymmetry in the stress and ligand density is observed in membranes during preparation at a water-vapor interface, with the stress asymmetry persisting in all membranes after drying.

Investigation of Changes in Tetracycline Repressor Binding upon Mutations in the Tetracycline Operator.

The tetracycline operon is an important gene network component, commonly used in synthetic biology applications because of its switch-like character. At the heart of this system is the highly specific interaction of the tet repressor protein (TetR) with its cognate DNA sequence (tetO). TetR binding on tetO practically stops expression of genes downstream of tetO by excluding RNA polymerase from binding the promoter and initiating transcription. Mutating the tetO sequence alters the strength of TetR-tetO binding and thus provides a tool to synthetic biologists to manipulate gene expression levels. We employ molecular dynamics (MD) simulations coupled with the free energy perturbation method to investigate the binding affinity of TetR to different tetO mutants. We also carry out in vivo tests in Escherichia coli for a series of promoters based on these mutants. We obtain reasonable agreement between experimental green fluorescent protein (GFP) repression levels and binding free energy differences computed from molecular simulations. In all cases, the wild-type tetO sequence yields the strongest TetR binding, which is observed both experimentally, in terms of GFP levels, and in simulation, in terms of free energy changes. Two of the four tetO mutants we tested yield relatively strong binding, whereas the other two mutants tend to be significantly weaker. The clustering and relative ranking of this subset of tetO mutants is generally consistent between our own experimental data, previous experiments with different systems and the free energy changes computed from our simulations. Overall, this work offers insights into an important synthetic biological system and demonstrates the potential, as well as limitations of molecular simulations to quantitatively explain biologically relevant behavior.

Effects of functional groups and ionization on the structure of alkanethiol-coated gold nanoparticles.

We report classical atomistic molecular dynamics simulations of alkanethiol-coated gold nanoparticles solvated in water and decane, as well as at water/vapor interfaces. The structure of the coatings is analyzed as a function of various functional end groups, including amine and carboxyl groups in various ionization states. We study both neutral and charged end groups for two different chain lengths (9 and 17 carbons). For the charged end groups, we simulated both mono- and divalent counterions. For the longer alkanes, we find significant local bundling of chains on the nanoparticle surface, which results in highly asymmetric coatings. In general, the charged end groups attenuate this effect by enhancing the water solubility of the nanoparticles. On the basis of the coating structures and density profiles, we can qualitatively infer the overall solubility of the nanoparticles. This asymmetry in the alkanethiol coatings is likely to have a significant effect on aggregation behavior. Our simulations elucidate the mechanism by which modulating the end group charge state can be used to control coating structure and therefore nanoparticle solubility and aggregation behavior.

Mechanically encoded cellular shapes for synthesis of anisotropic mesoporous particles.

The asymmetry that pervades molecular mechanisms of living systems increasingly informs the aims of synthetic chemistry, particularly in the development of catalysts, particles, nanomaterials, and their assemblies. For particle synthesis, overcoming viscous forces to produce complex, nonspherical shapes is particularly challenging; a problem that is continuously solved in nature when observing dynamic biological entities such as cells. Here we bridge these dynamics to synthetic chemistry and show that the intrinsic asymmetric shapes of erythrocytes can be directed, captured, and translated into composites and inorganic particles using a process of nanoscale silica-bioreplication. We show that crucial aspects in particle design such as particle-particle interactions, pore size, and macromolecular accessibility can be tuned using cellular responses. The durability of resultant particles provides opportunities for shape-preserving transformations into metallic, semiconductive, and ferromagnetic particles and assemblies. The ability to use cellular responses as "structure directing agents" offers an unprecedented toolset to design colloidal-scale materials.

Hydrogen-bonded aggregates in precise acid copolymers.

We perform atomistic molecular dynamics simulations of melts of four precise acid copolymers, two poly(ethylene-co-acrylic acid) (PEAA) copolymers, and two poly(ethylene-co-sulfonic acid) (PESA) copolymers. The acid groups are spaced by either 9 or 21 carbons along the polymer backbones. Hydrogen bonding causes the acid groups to form aggregates. These aggregates give rise to a low wavevector peak in the structure factors, in agreement with X-ray scattering data for the PEAA materials. The structure factors for the PESA copolymers are very similar to those for the PEAA copolymers, indicating a similar distance between aggregates which depends on the spacer length but not on the nature of the acid group. The PEAA copolymers are found to form more dimers and other small aggregates than do the PESA copolymers, while the PESA copolymers have both more free acid groups and more large aggregates.

High strength, molecularly thin nanoparticle membranes.

The unique strength observed in molecular thin films consisting of assemblies of nanoparticles encoded with short organic chains opens an intriguing new realm of controllable materials. Here the fundamental mechanisms underlying this unique mechanical strength are probed by molecular dynamics simulations. Using nanoparticles encoded with short hydrocarbon chains, we show that the mechanical response and failure of single nanoparticle thick membranes depend on subtle changes of the coating. Extremely high moduli were observed in agreement with experiment. We calculate Young's modulus for the membrane system based on properties of the individual components and find that ligand end-group interactions explain the observed changes in mechanical properties. Specifically, for dodecanethiol chains on 6 nm diameter gold cores, Young's modulus is 2.5 GPa for CH_{3} terminated chains and increases by 50% when end groups are replaced by COOH.

Atomistic Simulations Predict a Surprising Variety of Morphologies in Precise Ionomers.

The nature of ionic aggregates in ionomers remains an important open question, particularly considering its significance to their unique electrical and mechanical properties. We have carried out fully atomistic molecular dynamics simulations of melts of lithium-neutralized precise ionomers that reveal the structural features of ionic aggregates in unprecedented detail. In particular, we observe a rich variety of aggregate morphologies depending on neutralization level and ionic content, including string-like and percolated aggregates. The traditional assumption of spherical ionic aggregates with liquid-like ordering that is typically used to interpret experimental scattering data is too simplistic; a more rich and complex set of structures exist that also fit the scattering data.

Multiscale models of the antimicrobial peptide protegrin-1 on gram-negative bacteria membranes.

Antimicrobial peptides (AMPs) are naturally-occurring molecules that exhibit strong antibiotic properties against numerous infectious bacterial strains. Because of their unique mechanism of action, they have been touted as a potential source for novel antibiotic drugs. We present a summary of computational investigations in our lab aimed at understanding this unique mechanism of action, in particular the development of models that provide a quantitative connection between molecular-level biophysical phenomena and relevant biological effects. Our work is focused on protegrins, a potent class of AMPs that attack bacteria by associating with the bacterial membrane and forming transmembrane pores that facilitate the unrestricted transport of ions. Using fully atomistic molecular dynamics simulations, we have computed the thermodynamics of peptide-membrane association and insertion, as well as peptide aggregation. We also present a multi-scale analysis of the ion transport properties of protegrin pores, ranging from atomistic molecular dynamics simulations to mesoscale continuum models of single-pore electrodiffusion to models of transient ion transport from bacterial cells. Overall, this work provides a quantitative mechanistic description of the mechanism of action of protegrin antimicrobial peptides across multiple length and time scales.

No-slip boundary conditions and forced flow in multiparticle collision dynamics.

Multiparticle collision dynamics (MPCD) is a particle-based fluid simulation technique that is becoming increasingly popular for mesoscale fluid modeling. However, some confusion and conflicting results persist in literature regarding several important methodological details, in particular the enforcement of the no-slip condition and thermostatting in forced flow. These issues persist in simple flows past stationary boundaries, which we exclusively focus on here. We discuss the parametrization of MPCD fluids and its consequences for fluid-solid boundaries in great detail, and show that the method of virtual particles proposed by Lamura et al. and adopted by many others is required only for parameter choices that lead to viscosities dominated by collisional contributions. We test several implementations of the virtual particle method and discuss how to completely eliminate slip at stationary boundaries. We also show that stochastic boundary reflection rules are inherently problematic for forced flow and suggest a possible remedy. Finally, we discuss the most robust way to achieve forced flow and evaluate several thermostatting methods in the process. All discussion is limited to solid objects that do not move as a result of collisions with MPCD particles (i.e., walls). However, the results can be extended to solutes that experience forces and torques due to interactions with MPCD particles (e.g., colloids). The detailed analysis presented for this simple case provides the level of rigor and accuracy to the MPCD method required for the study of more complex systems.

Computational studies of protegrin antimicrobial peptides: a review.

Antimicrobial peptides (AMPs) are small, naturally occurring peptides that exhibit strong antibacterial properties generally believed to be a result of selective bacterial membrane disruption. As a result, there has been significant interest in the development of therapeutic antibiotics based on AMPs; however, the poor understanding of the fundamental mechanism of action of these peptides has largely hampered such efforts. We present a summary of computational and theoretical investigations of protegrin, a particularly potent peptide that is both an excellent model for the mechanism of action of AMPs and a promising therapeutic candidate. Experimental investigations have shed light on many of the key steps in the action of protegrin: protegrin monomers are known to dimerize in various lipid environments; protegrin peptides interact strongly with lipid bilayer membranes, particularly anionic lipids; protegrins have been shown to form pores in lipid bilayers, which results in uncontrolled ion transport and may be a key factor in bacterial death. In this work, we present a comprehensive review of the computational and theoretical studies that have complemented and extended the information obtained from experimental work with protegrins, as well as a brief survey of the experimental biophysical studies that are most pertinent to such computational work. We show that a consistent, mechanistic description of the bactericidal mechanism of action of protegrins is emerging, and briefly outline areas where the current understanding is deficient. We hope that the research reviewed herein offers compelling evidence of the benefits of computational investigations of protegrins and other AMPs, as well as providing a useful guide to future work in this area.

Poisson-Nernst-Planck models of nonequilibrium ion electrodiffusion through a protegrin transmembrane pore.

Protegrin peptides are potent antimicrobial agents believed to act against a variety of pathogens by forming nonselective transmembrane pores in the bacterial cell membrane. We have employed 3D Poisson-Nernst-Planck (PNP) calculations to determine the steady-state ion conduction characteristics of such pores at applied voltages in the range of -100 to +100 mV in 0.1 M KCl bath solutions. We have tested a variety of pore structures extracted from molecular dynamics (MD) simulations based on an experimentally proposed octomeric pore structure. The computed single-channel conductance values were in the range of 290-680 pS. Better agreement with the experimental range of 40-360 pS was obtained using structures from the last 40 ns of the MD simulation, where conductance values range from 280 to 430 pS. We observed no significant variation of the conductance with applied voltage in any of the structures that we tested, suggesting that the voltage dependence observed experimentally is a result of voltage-dependent channel formation rather than an inherent feature of the open pore structure. We have found the pore to be highly selective for anions, with anionic to cationic current ratios (I(Cl-)/I(K+)) on the order of 10(3). This is consistent with the highly cationic nature of the pore but surprisingly in disagreement with the experimental finding of only slight anionic selectivity. We have additionally tested the sensitivity of our PNP model to several parameters and found the ion diffusion coefficients to have a significant influence on conductance characteristics. The best agreement with experimental data was obtained using a diffusion coefficient for each ion set to 10% of the bulk literature value everywhere inside the channel, a scaling used by several other studies employing PNP calculations. Overall, this work presents a useful link between previous work focused on the structure of protegrin pores and experimental efforts aimed at investigating their conductance characteristics.