BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michael F. Hagan

Michael F. Hagan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softBiomoleculesBiological Physicscond-mat.stat-mechSubcellular Processescond-mat.mtrl-sciphysics.chem-ph

Catalog footprint

What is connected

39works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

39 published item(s)

preprint2022arXiv

Optimization of Non-Equilibrium Self-Assembly Protocols Using Markov State Models

The promise of self-assembly to enable bottom-up formation of materials with prescribed architectures and functions has driven intensive efforts to uncover rational design principles for maximizing the yield of a target structure. Yet, despite many successful examples of self-assembly, ensuring kinetic accessibility of the target structure remains an unsolved problem in many systems. In particular, long-lived kinetic traps can result in assembly times that vastly exceed experimentally accessible timescales. One proposed solution is to design non-equilibrium assembly protocols in which system parameters change over time to avoid such kinetic traps. Here, we develop a framework to combine Markov state model (MSM) analysis with optimal control theory to compute a time-dependent protocol that maximizes the yield of the target structure at a finite-time. We present an adjoint-based gradient descent method that, in conjunction with MSMs for a system as a function of its control parameters, enables efficiently optimizing the assembly protocol. We also describe an interpolation approach to significantly reduce the number of simulations required to construct the MSMs. We demonstrate our approach on two examples; a simple semi-analytic model for the folding of a polymer of colloidal particles, and a more complex model for capsid assembly. Our results show that optimizing time-dependent protocols can achieve significant improvements in yields of selected structures, including equilibrium free energy minima, long-lived metastable structures, and transient states.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Polymorphic self-assembly of helical tubules is kinetically controlled

In contrast to most self-assembling synthetic materials, which undergo unbounded growth, many biological self-assembly processes are self-limited. That is, the assembled structures have one or more finite dimensions that are much larger than the size scale of the individual monomers. In many such cases, the finite dimension is selected by a preferred curvature of the monomers, which leads to self-closure of the assembly. In this article, we study an example class of self-closing assemblies: cylindrical tubules that assemble from triangular monomers. By combining kinetic Monte Carlo simulations, free energy calculations, and simple theoretical models, we show that a range of programmable size scales can be targeted by controlling the intricate balance between the preferred curvature of the monomers and their interaction strengths. However, their assembly is kinetically controlled - the tubule morphology is essentially fixed shortly after closure, resulting in a distribution of tubule widths that is significantly broader than the equilibrium distribution. We develop a simple kinetic model based on this observation and the underlying free-energy landscape of assembling tubules that quantitatively describes the distributions. Our results are consistent with recent experimental observations of tubule assembly from triangular DNA origami monomers. The modeling framework elucidates design principles for assembling self-limited structures from synthetic components, such as artificial microtubules that have a desired width and chirality.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Tiling a tubule: How increasing complexity improves the yield of self-limited assembly

The ability to design and synthesize ever more complicated colloidal particles opens the possibility of self-assembling a zoo of complex structures, including those with one or more self-limited length scales. An undesirable feature of systems with self-limited length scales is that thermal fluctuations can lead to the assembly of nearby, off-target states. We investigate strategies for limiting off-target assembly by using multiple types of subunits. Using simulations and energetics calculations, we explore this concept by considering the assembly of tubules built from triangular subunits that bind edge to edge. While in principle, a single type of triangle can assemble into tubules with a monodisperse width distribution, in practice, the finite bending rigidity of the binding sites leads to the formation of off-target structures. To increase the assembly specificity, we introduce tiling rules for assembling tubules from multiple species of triangles. We show that the selectivity of the target structure can be dramatically improved by using multiple species of subunits, and provide a prescription for choosing the minimum number of subunit species required for near-perfect yield. Our approach of increasing the system's complexity to reduce the accessibility of neighboring structures should be generalizable to other systems beyond the self-assembly of tubules.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Vesicle shape transformations driven by confined active filaments

In active matter systems, deformable boundaries provide a mechanism to organize internal active stresses and perform work on the external environment. To study a minimal model of such a system, we perform particle-based simulations of an elastic vesicle containing a collection of polar active filaments. The interplay between the active stress organization due to interparticle interactions and that due to the deformability of the confinement leads to a variety of filament spatiotemporal organizations that have not been observed in bulk systems or under rigid confinement, including highly-aligned rings and caps. In turn, these filament assemblies drive dramatic and tunable transformations of the vesicle shape and its dynamics. We present simple scaling models that reveal the mechanisms underlying these emergent behaviors and yield design principles for engineering active materials with targeted shape dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Statistical properties of a tangentially driven active filament

Active polymers play a central role in many biological systems, from bacterial flagella to cellular cytoskeletons. Minimal models of semiflexible active filaments have been used to study a variety of interesting phenomena in active systems, such as defect dynamics in active nematics, clustering and laning in motility assays, and conformational properties of chromatin in eukaryotic cells. In this paper, we map a semiflexible polymer to an exactly solvable active Rouse chain, which enables us to analytically compute configurational and dynamical properties of active polymers with arbitrary rigidity.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Topological structure and dynamics of three dimensional active nematics

Point-like motile topological defects control the universal dynamics of diverse two-dimensional active nematics ranging from shaken granular rods to cellular monolayers. A comparable understanding in higher dimensions has yet to emerge. We report the creation of three-dimensional active nematics by dispersing extensile microtubule bundles in a passive colloidal liquid crystal. Light-sheet microscopy reveals the millimeter-scale structure of active nematics with a single bundle resolution and the temporal evolution of the associated nematic director field. The dominant excitations of three-dimensional active nematics are extended charge-neutral disclination loops that undergo complex dynamics and recombination events. These studies introduce a new class of non-equilibrium systems whose turbulent-like dynamics arises from the interplay between internally generated active stresses, the chaotic flows and the topological structure of the constituent defects.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Active Particles on Curved Surfaces

Recent studies have highlighted the sensitivity of active matter to boundaries and their geometries. Here we develop a general theory for the dynamics and statistics of active particles on curved surfaces and illustrate it on two examples. We first show that active particles moving on a surface with no ability to probe its curvature only exhibit steady-state inhomogeneities in the presence of orientational order. We then consider a strongly confined 3D ideal active gas and compute its steady-state density distribution in a box of arbitrary convex shape.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Allosteric control in icosahedral capsid assembly

During the lifecycle of a virus, viral proteins and other components self-assemble to form a symmetric protein shell called a capsid. This assembly process is subject to multiple competing constraints, including the need to form a thermostable shell while avoiding kinetic traps. It has been proposed that viral assembly satisfies these constraints through allosteric regulation, including the interconversion of capsid proteins among conformations with different propensities for assembly. In this article we use computational and theoretical modeling to explore how such allostery affects the assembly of icosahedral shells. We simulate assembly under a wide range of protein concentrations, protein binding affinities, and two different mechanisms of allosteric control. We find that, above a threshold strength of allosteric control, assembly becomes robust over a broad range of subunit binding affinities and concentrations, allowing the formation of highly thermostable capsids. Our results suggest that allostery can significantly shift the range of protein binding affinities that lead to successful assembly, and thus should be accounted for in models that are used to estimate interaction parameters from experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamics and density distribution of strongly confined noninteracting nonaligning self-propelled particles in a nonconvex boundary

We study the dynamics of non-aligning, non-interacting self-propelled particles confined to a box in two dimensions. In the strong confinement limit, when the persistence length of the active particles is much larger than the size of the box, particles stay on the boundary and align with the local boundary normal. It is then possible to derive the steady-state density on the boundary for arbitrary box shapes. In non-convex boxes, the non-uniqueness of the boundary normal results in hysteretic dynamics and the density is non-local, i.e. it depends on the global geometry of the box. These findings establish a general connection between the geometry of a confining box and the behavior of an ideal active gas it confines, thus providing a powerful tool to understand and design such confinements.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Emergent Self-organization in Active Materials

Biological systems exhibit large-scale self-organized dynamics and structures which enable organisms to perform the functions of life. The field of active matter strives to develop and understand microscopically-driven nonequilibrium materials, with emergent properties comparable to those of living systems. This review will describe two recently developed classes of active matter systems, in which simple building blocks --- self-propelled colloidal particles or extensile rod-like particles --- self-organize to form macroscopic structures with features not possible in equilibrium systems. We summarize the recent experimental and theoretical progress on each of these systems, and we present simple descriptions of the physics underlying their emergent behaviors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Equilibrium mappings in polar-isotropic confined active particles

Despite their fundamentally non-equilibrium nature, the individual and collective behavior of active systems with polar propulsion and isotropic interactions (polar-isotropic active systems) are remarkably well captured by equilibrium mapping techniques. Here we examine two signatures of equilibrium systems -- the existence of a local free energy function and the independence of the coarse- grained behavior on the details of the microscopic dynamics -- in polar-isotropic active particles confined by hard walls of arbitrary geometry at the one-particle level. We find that boundaries that possess concave regions make the density profile strongly dynamics-dependent and give it a nonlocal dependence on the geometry of the confining box. This in turn constrains the scope of equilibrium mapping techniques in polar-isotropic active systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Faceted particles formed by the frustrated packing of anisotropic colloids on curved surfaces

We use computer simulations and simple theoretical models to analyze the morphologies that result when rod-like particles end-attach onto a curved surface, creating a finite-thickness monolayer aligned with the surface normal. This geometry leads to two forms of frustration, one associated with the incompatibility of hexagonal order on surfaces with Gaussian curvature, and the second reflecting the deformation of a layer with finite thickness on a surface with non-zero mean curvature. We show that the latter effect leads to a faceting mechanism. Above threshold values of the inter-particle attraction strength and surface mean curvature, the adsorbed layer undergoes a transition from orientational disorder to an ordered state that is demarcated by reproducible patterns of line defects. The number of facets is controlled by the competition between line defect energy and intra-facet strain. Tuning control parameters thus leads to a rich variety of morphologies, including icosahedral particles and irregular polyhedra. In addition to suggesting a new strategy for the synthesis of aspherical particles with tunable symmetries, our results may shed light on recent experiments in which rod-like HIV GAG proteins assemble around nanoscale particles.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Interaction of chiral rafts in self-assembled colloidal membranes

Colloidal membranes are monolayer assemblies of rodlike particles that capture the long-wavelength properties of lipid bilayer membranes on the colloidal scale. Recent experiments on colloidal membranes formed by chiral rodlike viruses showed that introducing a second species of virus with different length and opposite chirality leads to the formation of rafts --- micron-sized domains of one virus species floating in a background of the other viruses [Sharma et al., Nature 513, 77 (2014)]. In this article we study the interaction of such rafts using liquid crystal elasticity theory. By numerically minimizing the director elastic free energy, we predict the tilt angle profile for both a single raft and two rafts in a background membrane, and the interaction between two rafts as a function of their separation. We find that the chiral penetration depth in the background membrane sets the scale for the range of the interaction. We compare our results with the experimental data and find good agreement for the strength and range of the interaction. Unlike the experiments, however, we do not observe a complete collapse of the data when rescaled by the tilt angle at the raft edge.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Recent advances in coarse-grained modeling of virus assembly

In many virus families, tens to thousands of proteins assemble spontaneously into a capsid (protein shell) while packaging the genomic nucleic acid. This review summarizes recent advances in computational modeling of these dynamical processes. We present an overview of recent technological and algorithmic developments, which are enabling simulations to describe the large ranges of length-and time-scales relevant to assembly, under conditions more closely matched to experiments than in earlier work. We then describe two examples in which computational modeling has recently provided an important complement to experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Allosteric Switching Mechanism in Bacteriophage MS2

In this article we use all-atom simulations to elucidate the mechanisms underlying conformational switching and allostery within the coat protein of the bacteriophage MS2. Assembly of most icosahedral virus capsids requires that the capsid protein adopt different conformations at precise locations within the capsid. It has been shown that a 19 nucleotide stem loop (TR) from the MS2 genome acts as an allosteric effector, guiding conformational switching of the coat protein during capsid assembly. Since the principal conformational changes occur far from the TR binding site, it is important to understand the molecular mechanism underlying this allosteric communication. To this end, we use all-atom simulations with explicit water combined with a path sampling technique to sample the MS2 coat protein conformational transition, in the presence and absence of TR-binding. The calculations find that TR binding strongly alters the transition free energy profile, leading to a switch in the favored conformation. We discuss changes in molecular interactions responsible for this shift. We then identify networks of amino acids with correlated motions to reveal the mechanism by which effects of TR binding span the protein. The analysis predicts amino acids whose substitution by mutagenesis could alter populations of the conformational substates or their transition rates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Measuring cohesion between macromolecular filaments, one pair at a time: Depletion-induced microtubule bundling

In presence of non-adsorbing polymers, colloidal particles experience ubiquitous attractive interactions induced by depletion forces. Here, we measure the depletion interaction between a pair of microtubule filaments using a method that combines single filament imaging with optical trapping. By quantifying the dependence of filament cohesion on both polymer concentration and solution ionic strength, we demonstrate that the minimal model of depletion, based on the Asakura-Oosawa theory, fails to quantitatively describe the experimental data. By measuring the cohesion strength in two- and three- filament bundles we verify pairwise additivity of depletion interactions for the specific experimental conditions. The described experimental technique can be used to measure pairwise interactions between various biological or synthetic filaments and complements information extracted from bulk osmotic stress experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Orientational order of motile defects in active nematics

The study of liquid crystals at equilibrium has led to fundamental insights into the nature of ordered materials, as well as to practical applications such as display technologies. Active nematics are a fundamentally different class of liquid crystals, driven away from equilibrium by the autonomous motion of their constituent rod-like particles. This internally generated activity powers the continuous creation and annihilation of topological defects, which leads to complex streaming flows whose chaotic dynamics appear to destroy long-range order. Here, we study these dynamics in experimental and computational realizations of active nematics. By tracking thousands of defects over centimetre-scale distances in microtubule-based active nematics, we identify a non-equilibrium phase characterized by system-spanning orientational order of defects. This emergent order persists over hours despite defect lifetimes of only seconds. Similar dynamical structures are observed in coarse-grained simulations, suggesting that defect-ordered phases are a generic feature of active nematics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Probing a self-assembled fd virus membrane with a microtubule

The self-assembly of highly anisotropic colloidal particles leads to a rich variety of morphologies, whose properties are just beginning to be understood. This article uses computer simulations to probe a particle-scale perturbation of a commonly studied colloidal assembly, a monolayer membrane composed of rodlike fd viruses in the presence of a polymer depletant. Motivated by experiments currently in progress, we simulate the interaction between a microtubule and a monolayer membrane as the microtubule "pokes" and penetrates the membrane face-on. Both the viruses and the microtubule are modeled as hard spherocylinders of the same diameter, while the depletant is modeled using ghost spheres. We find that the force exerted on the microtubule by the membrane is zero either when the microtubule is completely outside the membrane or when it has fully penetrated the membrane. The microtubule is initially repelled by the membrane as it begins to penetrate but experiences an attractive force as it penetrates further. We assess the roles played by translational and rotational fluctuations of the viruses and the osmotic pressure of the polymer depletant. We find that rotational fluctuations play a more important role than the translational ones. The dependence on the osmotic pressure of the depletant of the width and height of the repulsive barrier and the depth of the attractive potential well is consistent with the assumed depletion-induced attractive interaction between the microtubule and viruses. We discuss the relevance of these studies to the experimental investigations. \end{abstract}

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamics of Self-Propelled Particles Under Strong Confinement

We develop a statistical theory for the dynamics of non-aligning, non-interacting self-propelled particles confined in a convex box in two dimensions. We find that when the size of the box is small compared to the persistence length of a particle's trajectory (strong confinement), the steady-state density is zero in the bulk and proportional to the local curvature on the boundary. Conversely, the theory may be used to construct the box shape that yields any desired density distribution on the boundary. When the curvature variations are small, we also predict the distribution of orientations at the boundary and the exponential decay of pressure as a function of box size recently observed in 3D simulations in a spherical box.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hierarchical organization of chiral rafts in colloidal membranes

Liquid-liquid phase separation is ubiquitous in suspensions of nanoparticles, proteins and colloids. With a few notable exceptions, surface-tension-minimizing liquid droplets in bulk suspensions continuously coalesce, increasing in size without bound until achieving macroscale phase separation. In comparison, the phase behavior of colloids, nanoparticles or proteins confined to interfaces, surfaces or membranes is significantly more complex. Inclusions distort the local interface structure leading to interactions that are fundamentally different from the well-studied interactions mediated by isotropic solvents. Here, we investigate liquid-liquid phase separation in monolayer membranes composed of dissimilar chiral colloidal rods. We demonstrate that colloidal rafts are a ubiquitous feature of binary colloidal membranes. We measure the raft free energy landscape by visualizing its assembly kinetics. Subsequently, we quantify repulsive raft-raft interactions and relate them to directly imaged raft-induced membrane distortions, demonstrating that particle chirality plays a key role in this microphase separation. At high densities, rafts assemble into cluster crystals which constantly exchange monomeric rods with the background reservoir to maintain a self-limited size. Lastly, we demonstrate that rafts can form bonds to assemble into higher-order supra-structures. Our work demonstrates that membrane-mediated liquid-liquid phase separation can be fundamentally different from the well-characterized behavior of bulk liquids. It outlines a robust membrane-based pathway for assembly of monodisperse liquid clusters which is complementary to existing methods which take place in bulk suspensions. Finally, it reveals that chiral inclusions in membranes acquire long-ranged repulsive interactions, which might play a role in stabilizing assemblages of finite size.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Using Markov State Models to Study Self-Assembly

Markov state models (MSMs) have been demonstrated to be a powerful method for computationally studying intramolecular processes such as protein folding and macromolecular conformational changes. In this article, we present a new approach to construct MSMs that is applicable to modeling a broad class of multi-molecular assembly reactions. Distinct structures formed during assembly are distinguished by their undirected graphs, which are defined by strong subunit interactions. Spatial inhomogeneities of free subunits are accounted for using a recently developed Gaussian-based signature. Simplifications to this state identification are also investigated. The feasibility of this approach is demonstrated on two different coarse-grained models for virus self-assembly. We find good agreement between the dynamics predicted by the MSMs and long, unbiased simulations, and that the MSMs can reduce overall simulation time by orders of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Virus Assembly on a Membrane is Facilitated by Membrane Microdomains

For many viruses assembly and budding occur simultaneously during virion formation. Understanding the mechanisms underlying this process could promote biomedical efforts to block viral propagation and enable use of capsids in nanomaterials applications. To this end, we have performed molecular dynamics simulations on a coarse-grained model that describes virus assembly on a fluctuating lipid membrane. Our simulations show that the membrane can promote association of adsorbed subunits through dimensional reduction, but also can introduce barriers that inhibit complete assembly. We find several mechanisms, including one not anticipated by equilibrium theories, by which membrane microdomains, such as lipid rafts, can enhance assembly by reducing these barriers. We show how this predicted mechanism can be experimentally tested. Furthermore, the simulations demonstrate that assembly and budding depend crucially on the system dynamics via multiple timescales related to membrane deformation, protein diffusion, association, and adsorption onto the membrane.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modeling Viral Capsid Assembly

I present a review of the theoretical and computational methodologies that have been used to model the assembly of viral capsids. I discuss the capabilities and limitations of approaches ranging from equilibrium continuum theories to molecular dynamics simulations, and I give an overview of some of the important conclusions about virus assembly that have resulted from these modeling efforts. Topics include the assembly of empty viral shells, assembly around single-stranded nucleic acids to form viral particles, and assembly around synthetic polymers or charged nanoparticles for nanotechnology or biomedical applications. I present some examples in which modeling efforts have promoted experimental breakthroughs, as well as directions in which the connection between modeling and experiment can be strengthened.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reentrant Phase Behavior in Active Colloids with Attraction

Motivated by recent experiments, we study a system of self-propelled colloids that experience short-range attractive interactions and are confined to a surface. Using simulations we find that the phase behavior for such a system is reentrant as a function of activity: phase-separated states exist in both the low- and high-activity regimes, with a homogeneous active fluid in between. To understand the physical origins of reentrance, we develop a kinetic model for the system's steady-state dynamics whose solution captures the main features of the phase behavior. We also describe the varied kinetics of phase separation, which range from the familiar nucleation and growth of clusters to the complex coarsening of active particle gels.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Viral genome structures are optimal for capsid assembly

Understanding how virus capsids assemble around their nucleic acid (NA) genomes could promote efforts to block viral propagation or to reengineer capsids for gene therapy applications. We develop a coarse-grained model of capsid proteins and NAs with which we investigate assembly dynamics and thermodynamics. In contrast to recent theoretical models, we find that capsids spontaneously `overcharge'; i.e., the negative charge of the NA exceeds the positive charge on capsid. When applied to specific viruses, the optimal NA lengths closely correspond to the natural genome lengths. Calculations based on linear polyelectrolytes rather than base-paired NAs underpredict the optimal length, demonstrating the importance of NA structure to capsid assembly. These results suggest that electrostatics, excluded volume, and NA tertiary structure are sufficient to predict assembly thermodynamics and that the ability of viruses to selectively encapsidate their genomic NAs can be explained, at least in part, on a thermodynamic basis.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Diffusion-limited rates on low-dimensional manifolds with extreme aspect ratios

We consider a single-species diffusion-limited annihilation reaction with reactants confined to a two-dimensional surface with one arbitrarily large dimension and the other comparable in size to interparticle distances. This situation could describe reactants which undergo both longitudinal and transverse diffusion on long filamentous molecules (such as microtubules), or molecules that undergo truly one-dimensional translational diffusion (e.g. a transcription factor on DNA) but simultaneously exhibit diffusive behavior in a second dimension corresponding to a rotational or conformational degree of freedom. We combine simple analytical arguments and Monte Carlo simulations to show that the reaction rate law exhibits a crossover from one-dimensional to two-dimensional diffusion as a function of particle concentration and the size of the smaller dimension. In the case of a reversible binding reaction, the diffusion-limited reaction rate is given by the Smoluchowski expression, but the crossover is revealed in the statistics of particle collision histories. The results can also be applied to a particle-antiparticle annihilation reaction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Flagellar dynamics of a connected chain of active, Brownian particles

Eukaryotic flagella are active structures with a complex architecture of microtubules, motor proteins and elastic links. They are capable of whiplike motions driven by motors sliding along filaments that are themselves constrained at an end. Here, we show that active, self-propelled particles that are connected together to form a single chain that is anchored at one end can produce the graceful beating motions of flagella. We use a combination of numerical simulations, scaling analysis and mean field continuum elastic theory to demarcate the phase diagram for this type of oscillatory motion as a function of the filament length, passive elasticity, propulsion force and longitudinal persistence of propulsion directions. Depending on the nature of the anchoring, we show that filament either undergoes flagella-like beating or assumes a steadily rotating coiled conformation. Our system is simpler than its biological inspiration, and thus could be experimentally realized using a variety of self-propelled particles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mechanisms of budding of nanoscale particles through lipid bilayers

We examine the budding of a nanoscale particle through a lipid bilayer using molecular dynamics simulations, free energy calculations, and an elastic theory, with the aim of determining the extent to which equilibrium elasticity theory can describe the factors that control the mechanism and efficiency of budding. The particle is a smooth sphere which experiences attractive interactions to the lipid head groups. Depending on the parameters, we observe four classes of dynamical trajectories: particle adhesion to the membrane, stalled partially wrapped states, budding followed by scission, and membrane rupture. In most regions of parameter space we find that the elastic theory agrees nearly quantitatively with the simulated phase behavior as a function of adhesion strength, membrane bending rigidity, and particle radius. However, at parameter values near the transition between particle adhesion and budding, we observe long-lived partially wrapped states which are not captured by existing elastic theories. These states could constrain the accessible system parameters for those enveloped viruses or drug delivery vehicles which rely on exo- or endocytosis for membrane transport.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulations of HIV capsid protein dimerization reveal the effect of chemistry and topography on the mechanism of hydrophobic protein association

Recent work has shown that the hydrophobic protein surfaces in aqueous solution sit near a drying transition. The tendency for these surfaces to expel water from their vicinity leads to self assembly of macromolecular complexes. In this article we show with a realistic model for a biologically pertinent system how this phenomenon appears at the molecular level. We focus on the association of the C-terminal domain (CA-C) of the human immunodeficiency virus (HIV) capsid protein. By combining all-atom simulations with specialized sampling techniques we measure the water density distribution during the approach of two CA-C proteins as a function of separation and amino acid sequence in the interfacial region. The simulations demonstrate that CA-C protein-protein interactions sit at the edge of a dewetting transition and that this mesoscopic manifestation of the underlying liquid-vapor phase transition can be readily manipulated by biology or protein engineering to significantly affect association behavior. While the wild type protein remains wet until contact, we identify a set of in silico mutations, in which three hydrophilic amino acids are replaced with nonpolar residues, that leads to dewetting prior to association. The existence of dewetting depends on the size and relative locations of substituted residues separated by nm length scales, indicating long range cooperativity and a sensitivity to surface topography. These observations identify important details which are missing from descriptions of protein association based on buried hydrophobic surface area.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamical Self-regulation in Self-propelled Particle Flows

We study a continuum model of overdamped self-propelled particles with an aligning interaction in two dimensions. By combining analytical and numerical work, we map out the phase diagram for generic parameters. We find that the system self-organizes into two robust structures in different regions of parameter space: solitary waves of ordered swarms moving through a low density disordered background, and stationary asters. The self-regulating nature of the flow yields phase separation, ubiquitous in this class of systems, and controls the formation of solitary waves. Self-propulsion and the associated active convection play a crucial role in aster formation. A new result of our work is a phase diagram that displays these different regimes in a unified manner.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mechanisms of kinetic trapping in self-assembly and phase transformation

In self-assembly processes, kinetic trapping effects often hinder the formation of thermodynamically stable ordered states. In a model of viral capsid assembly and in the phase transformation of a lattice gas, we show how simulations in a self-assembling steady state can be used to identify two distinct mechanisms of kinetic trapping. We argue that one of these mechanisms can be adequately captured by kinetic rate equations, while the other involves a breakdown of theories that rely on cluster size as a reaction coordinate. We discuss how these observations might be useful in designing and optimising self-assembly reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Self-assembly of 2D membranes from mixtures of hard rods and depleting polymers

We elucidate the molecular forces leading to assembly of two dimensional membrane-like structures composed of a one rod-length thick monolayer of aligned rods from an immiscible suspension of hard rods and depleting polymers. We perform simulations which predict that monolayer membranes are thermodynamically stable above a critical rod aspect ratio and below a critical depletion interaction length scale. Outside of these conditions alternative structures such as stacked smectic columns or nematic droplets are thermodynamically stable. These predictions are confirmed using an experimental model system of virus rod-like molecules and non-adsorbing polymer. Our work demonstrates that collective molecular protrusion fluctuations alone are sufficient to stabilize membranes composed of homogenous rods with simple excluded volume interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sitting at the edge: How biomolecules use hydrophobicity to tune their interactions and function

Water near hydrophobic surfaces is like that at a liquid-vapor interface, where fluctuations in water density are substantially enhanced compared to that in bulk water. Here we use molecular simulations with specialized sampling techniques to show that water density fluctuations are similarly enhanced, even near hydrophobic surfaces of complex biomolecules, situating them at the edge of a dewetting transition. Consequently, water near these surfaces is sensitive to subtle changes in surface conformation, topology, and chemistry, any of which can tip the balance towards or away from the wet state, and thus significantly alter biomolecular interactions and function. Our work also resolves the long-standing puzzle of why some biological surfaces dewet and other seemingly similar surfaces do not.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spontaneous Segregation of Self-Propelled Particles with Different Motilities

We study mixtures of self-propelled and passive rod-like particles in two dimensions using Brownian dynamics simulations. The simulations demonstrate that the two species spontaneously segregate to generate a rich array of dynamical domain structures whose properties depend on the propulsion velocity, density, and composition. In addition to presenting phase diagrams as a function of the system parameters, we investigate the mechanisms driving segregation. We show that the difference in collision frequencies between self-propelled and passive rods provides a driving force for segregation, which is amplified by the tendency of the self-propelled rods to swarm or cluster. Finally, both self-propelled and passive rods exhibit giant number fluctuations for sufficient propulsion velocities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Theoretical calculation of the phase behavior of colloidal membranes

We formulate a density functional theory that describes the phase behavior of hard rods and depleting polymers, as realized in recent experiments on suspensions of \emph{fd} virus and non-adsorbing polymer. The theory predicts the relative stability of nematic droplets, stacked smectic columns, and a recently discovered phase of isolated monolayers of rods, or colloidal membranes. We find that a minimum rod aspect ratio is required for stability of colloidal membranes and that collective protrusion undulations are the dominant effect that stabilizes this phase. The theoretical predictions are shown to be qualitatively consistent with experimental and computational results.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Encapsulation of a polymer by an icosahedral virus

The coat proteins of many viruses spontaneously form icosahedral capsids around nucleic acids or other polymers. Elucidating the role of the packaged polymer in capsid formation could promote biomedical efforts to block viral replication and enable use of capsids in nanomaterials applications. To this end, we perform Brownian dynamics on a coarse-grained model that describes the dynamics of icosahedral capsid assembly around a flexible polymer. We identify several mechanisms by which the polymer plays an active role in its encapsulation, including cooperative polymer-protein motions. These mechanisms are related to experimentally controllable parameters such as polymer length, protein concentration, and solution conditions. Furthermore, the simulations demonstrate that assembly mechanisms are correlated to encapsulation efficiency, and we present a phase diagram that predicts assembly outcomes as a function of experimental parameters. We anticipate that our simulation results will provide a framework for designing in vitro assembly experiments on single-stranded RNA virus capsids.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mechanisms of viral capsid assembly around a polymer

Capsids of many viruses assemble around nucleic acids or other polymers. Understanding how the properties of the packaged polymer affect the assembly process could promote biomedical efforts to prevent viral assembly or nanomaterials applications that exploit assembly. To this end, we simulate on a lattice the dynamical assembly of closed, hollow shells composed of several hundred to 1000 subunits, around a flexible polymer. We find that assembly is most efficient at an optimum polymer length that scales with the surface area of the capsid; significantly longer than optimal polymers often lead to partial-capsids with unpackaged polymer `tails' or a competition between multiple partial-capsids attached to a single polymer. These predictions can be tested with bulk experiments in which capsid proteins assemble around homopolymeric RNA or synthetic polyelectrolytes. We also find that the polymer can increase the net rate of subunit accretion to a growing capsid both by stabilizing the addition of new subunits and by enhancing the incoming flux of subunits; the effects of these processes may be distinguishable with experiments that monitor the assembly of individual capsids.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Self-limited self-assembly of chiral filaments

The assembly of filamentous bundles with controlled diameters is common in biological systems and desirable for the development of nanomaterials. We discuss dynamical simulations and free energy calculations on patchy spheres with chiral pair interactions that spontaneously assemble into filamentous bundles. The chirality frustrates long-range crystal order by introducing twist between interacting subunits. For some ranges of system parameters this constraint leads to bundles with a finite diameter as the equilibrium state, and in other cases frustration is relieved by the formation of defects. While some self-limited structures can be modeled as twisted filaments arranged with local hexagonal symmetry, other structures are surprising in their complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Understanding the Concentration Dependence of Viral Capsid Assembly Kinetics - the Origin of the Lag Time and Identifying the Critical Nucleus Size

The kinetics for the assembly of viral proteins into a population of capsids can be measured in vitro with size exclusion chromatography or dynamic light scattering, but extracting mechanistic information from these studies is challenging. For example, it is not straightforward to determine the critical nucleus size or the elongation time (the time required for a nucleated partial capsid to grow completion). We show that, for two theoretical models of capsid assembly, the critical nucleus size can be determined from the concentration dependence of the assembly reaction half-life and the elongation time is revealed by the length of the lag phase. Furthermore, we find that the system becomes kinetically trapped when nucleation becomes fast compared to elongation. Implications of this constraint for determining elongation mechanisms from experimental assembly data are discussed.

0 citations0 reviewsNo code linkedOpen access