Dynamics of Growth and Form in Prebiotic Vesicles

The growth, form, and division of prebiotic vesicles, membraneous bags of fluid of varying components and shapes is hypothesized to have served as the substrate for the origin of life. The dynamics of these out-of-equilibrium structures is controlled by physicochemical processes that include the intercalation of amphiphiles into the membrane, fluid flow across the membrane, and elastic deformations of the membrane. To understand prebiotic vesicular forms and their dynamics, we construct a minimal model that couples membrane growth, deformation, and fluid permeation, ultimately couched in terms of two dimensionless parameters that characterize the relative rate of membrane growth and the membrane permeability. Numerical simulations show that our model captures the morphological diversity seen in extant precursor mimics of cellular life, and might provide simple guidelines for the synthesis of these complex shapes from simple ingredients.

Organ size control via hydraulically gated oscillations

Hollow vesicular tissues of various sizes and shapes arise in biological organs such as ears, guts, hearts, brains and even entire organisms. Regulating their size and shape is crucial for their function. Although chemical signaling has been thought to play a role in the regulation of cellular processes that feed into larger scales, it is increasingly recognized that mechanical forces are involved in the modulation of size and shape at larger length scales. Motivated by a variety of examples of tissue cyst formation and size control that show simultaneous growth and size oscillations, we create a minimal theoretical framework for the growth and dynamics of a soft, fluid-permeable, spherical shell. We show that these shells can relieve internal pressure by bursting intermittently, shrinking and re-growing, providing a simple mechanism by which hydraulically gated oscillations can regulate size. To test our theory, we develop an in vitro experimental set-up to monitor the growth and oscillations of a hollow tissue spheroid growing freely or when confined. A simple generalization of our theory to account for irreversible deformations allows us to explain the time scales and the amplitudes of oscillations in terms of the geometry and mechanical properties of the tissue shells. Taken together, our theory and experimental observations show how soft hydraulics can regulate the size of growing tissue shells.

SiteOut: An Online Tool to Design Binding Site-Free DNA Sequences

DNA-binding proteins control many fundamental biological processes such as transcription, recombination and replication. A major goal is to decipher the role that DNA sequence plays in orchestrating the binding and activity of such regulatory proteins. To address this goal, it is useful to rationally design DNA sequences with desired numbers, affinities and arrangements of protein binding sites. However, removing binding sites from DNA is computationally non-trivial since one risks creating new sites in the process of deleting or moving others. Here we present an online binding site removal tool, SiteOut, that enables users to design arbitrary DNA sequences that entirely lack binding sites for factors of interest. SiteOut can also be used to delete sites from a specific sequence, or to introduce site-free spacers between functional sequences without creating new sites at the junctions. In combination with commercial DNA synthesis services, SiteOut provides a powerful and flexible platform for synthetic projects that interrogate regulatory DNA. Here we describe the algorithm and illustrate the ways in which SiteOut can be used; it is publicly available at https://depace.med.harvard.edu/siteout/.

Simulations show that virus assembly and budding are 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 it also introduces thermodynamic and kinetic effects that can inhibit complete assembly. We find several mechanisms by which membrane microdomains, such as lipid rafts, reduce these effects, and thus, enhance assembly. We show how these predicted mechanisms 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.

A tunable coarse-grained model for ligand-receptor interaction

Cell-surface receptors are the most common target for therapeutic drugs. The design and optimization of next generation synthetic drugs require a detailed understanding of the interaction with their corresponding receptors. Mathematical approximations to study ligand-receptor systems based on reaction kinetics strongly simplify the spatial constraints of the interaction, while full atomistic ligand-receptor models do not allow for a statistical many-particle analysis, due to their high computational requirements. Here we present a generic coarse-grained model for ligand-receptor systems that accounts for the essential spatial characteristics of the interaction, while allowing statistical analysis. The model captures the main features of ligand-receptor kinetics, such as diffusion dependence of affinity and dissociation rates. Our model is used to characterize chimeric compounds, designed to take advantage of the receptor over-expression phenotype of certain diseases to selectively target unhealthy cells. Molecular dynamics simulations of chimeric ligands are used to study how selectivity can be optimized based on receptor abundance, ligand-receptor affinity and length of the linker between both ligand subunits. Overall, this coarse-grained model is a useful approximation in the study of systems with complex ligand-receptor interactions or spatial constraints.

The current importance of cell surface receptors as primary targets for drug treatment explains the increasing interest in a mathematical and quantitative description of the process of ligand-receptor interaction. Recently, a new generation of synthetic chimeric ligands has been developed to selectively target unhealthy cells, without harming healthy tissue. To understand these and other types of complex ligand-receptor systems, conventional chemical interaction models often rely on simplifications and assumptions about the spatial characteristics of the system, while full atomistic molecular dynamics simulations are too computationally demanding to perform many particle statistical analysis. In this paper, we describe a novel approach to model the interaction between ligands and receptors based on a coarse grained approximation that includes explicitly both spatial and kinetic details of the interaction, while allowing many-particle simulations and therefore, statistical analysis at biologically relevant time scales. The model is used to study the binding properties of generic chimeric ligands to understand how cell specificity is achieved, its dependence on receptor concentration and the influence of the distance between subunits of the chimera. Overall, this approach proves optimal to study other ligand-receptor systems with complex spatial regulation, such as receptor clustering, multimerization and multivalent asymmetric ligand binding.

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.