Computational Approaches to Modeling Viral Structure and Assembly

Biomolecular engineering of virus-like particles aided by computational chemistry methods

Multi-scale investigation of VLP self-assembly aided by computational methods is facilitating the design, redesign, and modification of functionalized VLPs.

Perspective: Coarse-grained models for biomolecular systems

By focusing on essential features, while averaging over less important details, coarse-grained (CG) models provide significant computational and conceptual advantages with respect to more detailed models. Consequently, despite dramatic advances in computational methodologies and resources, CG models enjoy surging popularity and are becoming increasingly equal partners to atomically detailed models. This perspective surveys the rapidly developing landscape of CG models for biomolecular systems. In particular, this review seeks to provide a balanced, coherent, and unified presentation of several distinct approaches for developing CG models, including top-down, network-based, native-centric, knowledge-based, and bottom-up modeling strategies. The review summarizes their basic philosophies, theoretical foundations, typical applications, and recent developments. Additionally, the review identifies fundamental inter-relationships among the diverse approaches and discusses outstanding challenges in the field. When carefully applied and assessed, current CG models provide highly efficient means for investigating the biological consequences of basic physicochemical principles. Moreover, rigorous bottom-up approaches hold great promise for further improving the accuracy and scope of CG models for biomolecular systems.

Comment on the letter by A. Ben-Shaul: "entropy, energy, and bending of DNA in viral capsids"

The conformational entropic penalty associated with packaging double-stranded DNA into viral capsids remains an issue of contention. So far, models based on a continuum approximation for DNA have either left the question unexamined, or they have assumed that the entropic penalty is negligible, following an early analysis by Riemer and Bloomfield. In contrast, molecular-dynamics (MD) simulations using bead-and-spring models consistently show a large penalty. A recent letter from Ben-Shaul attempts to reconcile the differences. While the letter makes some valid points, the issue of how to include conformational entropy in the continuum models remains unresolved. In this Comment, I show that the free energy decomposition from continuum models could be brought into line with the decomposition from the MD simulations with two adjustments. First, the entropy from Flory-Huggins theory should be replaced by the estimate of the entropic penalty given in Ben-Shaul's letter, which corresponds closely to that from the MD simulations. Second, the DNA-DNA repulsions are well described by the empirical relationship given by the Cal Tech group, but the strength of these should be reduced by about half, using parameters based on the Rau-Parsegian experiments, rather than treating them as "fitting parameters (tuned) to fit the data from (single molecule pulling) experiments."

Effects of pulling forces, osmotic pressure, condensing agents and viscosity on the thermodynamics and kinetics of DNA ejection from bacteriophages to bacterial cells: a computational study

In this work, we report on simulations of double-stranded DNA (dsDNA) ejection from bacteriophage φ29 into a bacterial cell. The ejection was studied with a coarse-grained model, in which viral dsDNA was represented by beads on a torsion-less string. The bacteriophage's capsid and the bacterial cell were defined by sets of spherical constraints. To account for the effects of the viscous medium inside the bacterial cell, the simulations were carried out using a Langevin dynamics protocol. Our simplest simulations (involving constant viscosity and no external biasing forces) produced results compatible with the push-pull model of DNA ejection, with an ejection rate significantly higher in the first part of ejection than in the latter parts. Additionally, we performed more complicated simulations, in which we included additional factors such as external forces, osmotic pressure, condensing agents and ejection-dependent viscosity. The effects of these factors (independently and in combination) on the thermodynamics and kinetics of DNA ejection were studied. We found that, in general, the dependence of ejection forces and ejection rates on the amount of DNA ejected becomes more complex if the ejection is modeled with a broader, more realistic set of parameters and influences (such as variation in the solvent's viscosity and the application of an external force). However, certain combinations of factors and numerical parameters led to the opposition of some ejection-driving and ejection-inhibiting influences, ultimately causing an apparent simplification of the ejection profiles.

The entropic cost of polymer confinement

The confinement of a polymer into a small space is thermodynamically unfavorable because of the reduction in the number of conformational states. The entropic penalty affects a variety of biological processes, and it plays an important role in polymer transport properties and in microfluidic devices. We determine the entropic penalty for the confinement of elastic polymer of persistence length P in the long-chain limit. We examine three geometries: (1) parallel planes separated by a distance d (a slit); (2) a circular tube of diameter d; and (3) a sphere of diameter d. We first consider infinitely thin (ideal) chains. As d/P drops from 100 to 0.01, TΔS rises from ∼5 × 10(-4) kT to ∼30 kT per persistence length for cases (1) and (2), with the entropic penalty for case (2) being consistently about twice that for case (1). TΔS is ∼5 kT per persistence length for confinement to a sphere when d = P, about twice the value predicted by mean field theory. For all three geometries, in the limit d/P ≫ 1, the asymptotic behavior of ΔS vs d is consistent with the d(-2) behavior predicted by theory. In the limit d/P ≪ 1, the scaling of ΔS for slits and tubes is also consistent with earlier predictions (d(-2/3)). Finally, we treat volume exclusion effects, examining chains of diameter D > 0. Confinement to a narrow slit or tube (d/P ≪ 1) has the same entropic penalty as that for an ideal chain in a slit or tube with d' = d - D; in the weak confinement regime (d/P ≫ 1), the entropic penalties are significantly larger than those for infinitely thin chains. When a chain of finite diameter is forced into a sphere or other closed cavity, the entropic confinement penalty rises without limit because there are no configurations available to the chain once its volume exceeds that of the cavity.

Coarse-graining of multiprotein assemblies

Multiscale models are important tools to elucidate how small changes in local subunit conformations may propagate to affect the properties of macromolecular complexes. We review recent advances in coarse-graining methods for poly-protein assemblies, systems that are composed of many copies of relatively few components, with a particular focus on viral capsids and cytoskeletal filaments. These methods are grouped into two broad categories-mapping methods, which use information from one scale of representation to parameterize a lower resolution model, and bridging methods, which repeatedly connect different scales during simulation-and we provide examples of both classes at different levels of complexity. Collectively, these models illustrate the numerous approaches to information transfer between scales and demonstrate that the complexity required of the model depends in general on the nature of the information sought.