Defending Our Public Biological Databases as a Global Critical Infrastructure.

Progress in modern biology is being driven, in part, by the large amounts of freely available data in public resources such as the International Nucleotide Sequence Database Collaboration (INSDC), the world's primary database of biological sequence (and related) information. INSDC and similar databases have dramatically increased the pace of fundamental biological discovery and enabled a host of innovative therapeutic, diagnostic, and forensic applications. However, as high-value, openly shared resources with a high degree of assumed trust, these repositories share compelling similarities to the early days of the Internet. Consequently, as public biological databases continue to increase in size and importance, we expect that they will face the same threats as undefended cyberspace. There is a unique opportunity, before a significant breach and loss of trust occurs, to ensure they evolve with quality and security as a design philosophy rather than costly "retrofitted" mitigations. This Perspective surveys some potential quality assurance and security weaknesses in existing open genomic and proteomic repositories, describes methods to mitigate the likelihood of both intentional and unintentional errors, and offers recommendations for risk mitigation based on lessons learned from cybersecurity.

Metastable Prepores in Tension-Free Lipid Bilayers.

The formation and closure of aqueous pores in lipid bilayers is a key step in various biophysical processes. Large pores are well described by classical nucleation theory, but the free-energy landscape of small, biologically relevant pores has remained largely unexplored. The existence of small and metastable "prepores" was hypothesized decades ago from electroporation experiments, but resolving metastable prepores from theoretical models remained challenging. Using two complementary methods-atomistic simulations and self-consistent field theory of a minimal lipid model-we determine the parameters for which metastable prepores occur in lipid membranes. Both methods consistently suggest that pore metastability depends on the relative volume ratio between the lipid head group and lipid tails: lipids with a larger head-group volume fraction (or shorter saturated tails) form metastable prepores, whereas lipids with a smaller head-group volume fraction (or longer unsaturated tails) form unstable prepores.

Membrane stress profiles from self-consistent field theory.

Using self-consistent field theory (SCFT), we develop an accurate, local expression for the stress profiles in membranes and soft matter interfaces, in general. The bond stresses are expressed in terms of pre-computed chain propagators, which are used to describe the statistical weight of the molecules and therefore require minimal additional calculations. In addition, we overcome the resolution limit of the molecular bond length by including the Irving and Kirkwood bond assignment and recover a constant normal stress profile across an interface. Using this theory, we find that the membrane lateral stress profile contains repulsive (positive) stresses in the regions of the head and tail groups, and attractive (negative) stresses near the hydrophobic/hydrophilic interface. We also verify that the zeroth and first moments of the stress profile correspond to the thermodynamic tension and product of the bending modulus and the spontaneous curvature, respectively.

Nonequilibrium simulations of model ionomers in an oscillating electric field.

We perform molecular dynamics simulations of a coarse-grained model of ionomer melts in an applied oscillating electric field. The frequency-dependent conductivity and susceptibility are calculated directly from the current density and polarization density, respectively. At high frequencies, we find a peak in the real part of the conductivity due to plasma oscillations of the ions. At lower frequencies, the dynamic response of the ionomers depends on the ionic aggregate morphology in the system, which consists of either percolated or isolated aggregates. We show that the dynamic response of the model ionomers to the applied oscillating field can be understood by comparison with relevant time scales in the systems, obtained from independent calculations.

Amphiphilic triblocks to control assembly of mixed or segregated bilayers and monolayers.

Triblock amphiphilic molecules composed of three distinct segments provide a large parameter space to obtain self-assembled structures beyond what is achievable with conventional amphiphiles. To obtain a molecular understanding of the thermodynamics of self-assembly, we develop a coarse-grained triblock polymer model and apply self-consistent field theory to investigate the packing mechanism into layer structures. By tuning the structural and interaction asymmetry, we are able to obtain bilayers and monolayers, where the latter may additionally be mixed (symmetric) or segregated (asymmetric). Of particular interest for a variety of applications are the asymmetric monolayers, where segregation of end blocks to opposite surfaces is expected to have important implications for the development of functional nanotubes and vesicles with distinct surface chemistries.

Electrostatically tuned self-assembly of branched amphiphilic peptides.

Electrostatics plays an important role in the self-assembly of amphiphilic peptides. To develop a molecular understanding of the role of the electrostatic interactions, we develop a coarse-grained model peptide and apply self-consistent field theory to investigate the peptide assembly into a variety of aggregate nanostructures. We find that the presence and distribution of charged groups on the hydrophilic branches of the peptide can modify the molecular configuration from extended to collapsed. This change in molecular configuration influences the packing into spherical micelles, cylindrical micelles (nanofibers), or planar bilayers. The effects of charge distribution therefore have important implications for the design and utility of functional materials based on peptides.

Nucleation in polymers and soft matter.

Nucleation is a ubiquitous phenomenon in many physical, chemical, and biological processes. In this review, we describe recent progress on the theoretical study of nucleation in polymeric fluids and soft matter, including binary mixtures (polymer blends, polymers in poor solvents, compressible polymer-small molecule mixtures), block copolymer melts, and lipid membranes. We discuss the methodological development for studying nucleation as well as novel insights and new physics obtained in the study of the nucleation behavior in these systems.

Activated pathways for the directed insertion of patterned nanoparticles into polymer membranes.

We combine the string method with self-consistent field theory to compute the most probable transition pathway, i.e. the minimum free energy path, for the insertion of Janus and protein-like nanoparticles into a polymer membrane bilayer. The method makes no assumptions in the reaction coordinate and overcomes the long timescales challenge associated with simulating rare events. Our study suggests that one approach to building functional polymer­nanoparticle composite membranes with oriented nanoparticles is through electrostatic interactions. In particular, hydrophobic Janus nanoparticles with an asymmetric charge distribution can be made to directionally insert into charged membranes. This process is kinetically driven, and involves overcoming a thermally surmountable activation barrier, which requires favorable interactions between the nanoparticle and the hydrophilic block of the membrane. In contrast, the insertion of protein-like nanoparticles with alternating hydrophilic­hydrophobic­hydrophilic domains into polymer membranes does not occur as a thermally activated event.

Thermodynamic basis for the genome to capsid charge relationship in viral encapsidation.

We establish an appropriate thermodynamic framework for determining the optimal genome length in electrostatically driven viral encapsidation. Importantly, our analysis includes the electrostatic potential due to the Donnan equilibrium, which arises from the semipermeable nature of the viral capsid, i.e., permeable to small mobile ions but impermeable to charged macromolecules. Because most macromolecules in the cellular milieu are negatively charged, the Donnan potential provides an additional driving force for genome encapsidation. In contrast to previous theoretical studies, we find that the optimal genome length is the result of combined effects from the electrostatic interactions of all charged species, the excluded volume and, to a very significant degree, the Donnan potential. In particular, the Donnan potential is essential for obtaining negatively overcharged viruses. The prevalence of overcharged viruses in nature may suggest an evolutionary preference for viruses to increase the amount of genome packaged by utilizing the Donnan potential (through increases in the capsid radius), rather than high charges on the capsid, so that structural stability of the capsid is maintained.

Minimum energy path to membrane pore formation and rupture.

We combine dynamic self-consistent field theory with the string method to calculate the minimum energy path to membrane pore formation and rupture. In the regime where nucleation can occur on experimentally relevant time scales, the structure of the critical nucleus is between a solvophilic stalk and a locally thinned membrane. Classical nucleation theory fails to capture these molecular details and significantly overestimates the free energy barrier. Our results suggest that thermally nucleated rupture may be an important factor for the low rupture strains observed in lipid membranes.

Interactions of a charged nanoparticle with a lipid membrane: implications for gene delivery.

We employ self-consistent field theory to study the thermodynamics of membrane-particle interactions in the context of gene delivery systems, with the aim to guide the design of dendrimers that can overcome the endosomal escape barrier by inserting into membranes and creating pores. We consider the interaction between a model polyamidoamine dendrimer and a membrane under controlled tension as a function of the separation between the dendrimer and the membrane. In all the cases we have studied, the lowest free energy state corresponds to the membrane partially wrapping the dendrimer. However, with moderate tension, we find that a G5 (or larger) generation dendrimer, through thermal fluctuation, can induce the formation of metastable pores. These metastable pores are subsequently shown to significantly lower the critical tension necessary for membrane rupture, thus possibly enhancing the release of the trapped genetic material from the endosome.

A kinetic model for the enzymatic action of cellulase.

We develop a mechanochemical model for the dynamics of cellulase, a two-domain enzyme connected by a peptide linker, as it extracts and hydrolyzes a cellulose polymer from a crystalline substrate. We consider two random walkers, representing the catalytic domain (CD) and the carbohydrate binding module (CBM), whose rates for stepping are biased by the coupling through the linker and the energy required to lift the cellulose polymer from the crystalline surface. Our results show that the linker length and stiffness play a critical role in the cooperative action of the CD and CBM domains and that, for a given linker length, the steady-state hydrolysis shows a maximum at some intermediate linker stiffness. The maximum hydrolysis rate corresponds to a transition of the linker from a compressed to an extended conformation, where the system exhibits maximum fluctuation, as measured by the variance of the separation distance between the two domains and the dispersion around the mean hydrolysis speed. In the range of experimentally known values of the parameters of our model, improving the intrinsic hydrolytic activity of the CD leads to a proportional increase in the overall hydrolysis rate.

Two-dimensional fluorescence resonance energy transfer as a probe for protein folding: a theoretical study.

We describe a two-dimensional (2D), four-color fluorescence resonance energy transfer (FRET) scheme, in which the conformational dynamics of a protein is followed by simultaneously observing the FRET signal from two different donor-acceptor pairs. For a general class of models that assume Markovian conformational dynamics, we relate the properties of the emission correlation functions to the rates of elementary kinetic steps in the model. We further use a toy folding model that treats proteins as chains with breakable cross-links to examine the relationship between the cooperativity of folding and FRET data and to establish what additional information about the folding dynamics can be gleaned from 2D, as opposed to one-dimensional FRET experiments. We finally discuss the potential advantages of the four-color FRET over the three-color FRET technique.