Tabulation as a high-resolution alternative to coarse-graining protein interactions: Initial application to virus capsid subunits

Markovian Weighted Ensemble Milestoning (M-WEM): Long-Time Kinetics from Short Trajectories

We introduce a rare-event sampling scheme, named Markovian Weighted Ensemble Milestoning (M-WEM), which inlays a weighted ensemble framework within a Markovian milestoning theory to efficiently calculate thermodynamic and kinetic properties of long-time-scale biomolecular processes from short atomistic molecular dynamics simulations. M-WEM is tested on the Müller-Brown potential model, the conformational switching in alanine dipeptide, and the millisecond time-scale protein-ligand unbinding in a trypsin-benzamidine complex. Not only can M-WEM predict the kinetics of these processes with quantitative accuracy but it also allows for a scheme to reconstruct a multidimensional free-energy landscape along additional degrees of freedom, which are not part of the milestoning progress coordinate. For the ligand-receptor system, the experimental residence time, association and dissociation kinetics, and binding free energy could be reproduced using M-WEM within a simulation time of a few hundreds of nanoseconds, which is a fraction of the computational cost of other currently available methods, and close to 4 orders of magnitude less than the experimental residence time. Due to the high accuracy and low computational cost, the M-WEM approach can find potential applications in kinetics and free-energy-based computational drug design.

Should Virus Capsids Assemble Perfectly? Theory and Observation of Defects

Although published structural models of viral capsids generally exhibit a high degree of regularity or symmetry, structural defects might be expected because of the fluctuating environment in which capsids assemble and the requirement of some capsids for disassembly before genome delivery. Defective structures are observed in computer simulations, and are evident in single-particle cryoelectron microscopy studies. Here, we quantify the conditions under which defects might be expected, using a statistical mechanics model allowing for ideal, defective, and vacant sites. The model displays a threshold in affinity parameters below which there is an appreciable population of defective capsids. Even when defective sites are not allowed, there is generally some population of vacancies. Analysis of single particles in cryoelectron microscopy micrographs yields a confirmatory ≳15% of defective particles. Our findings suggest structural heterogeneity in virus capsids may be under-appreciated, and also points to a nontraditional strategy for assembly inhibition.

Quantitative computational models of molecular self-assembly in systems biology

Molecular self-assembly is the dominant form of chemical reaction in living systems, yet efforts at systems biology modeling are only beginning to appreciate the need for and challenges to accurate quantitative modeling of self-assembly. Self-assembly reactions are essential to nearly every important process in cell and molecular biology and handling them is thus a necessary step in building comprehensive models of complex cellular systems. They present exceptional challenges, however, to standard methods for simulating complex systems. While the general systems biology world is just beginning to deal with these challenges, there is an extensive literature dealing with them for more specialized self-assembly modeling. This review will examine the challenges of self-assembly modeling, nascent efforts to deal with these challenges in the systems modeling community, and some of the solutions offered in prior work on self-assembly specifically. The review concludes with some consideration of the likely role of self-assembly in the future of complex biological system models more generally.

Weighted Ensemble Simulation: Review of Methodology, Applications, and Software

The weighted ensemble (WE) methodology orchestrates quasi-independent parallel simulations run with intermittent communication that can enhance sampling of rare events such as protein conformational changes, folding, and binding. The WE strategy can achieve superlinear scaling-the unbiased estimation of key observables such as rate constants and equilibrium state populations to greater precision than would be possible with ordinary parallel simulation. WE software can be used to control any dynamics engine, such as standard molecular dynamics and cell-modeling packages. This article reviews the theoretical basis of WE and goes on to describe successful applications to a number of complex biological processes-protein conformational transitions, (un)binding, and assembly processes, as well as cell-scale processes in systems biology. We furthermore discuss the challenges that need to be overcome in the next phase of WE methodological development. Overall, the combined advances in WE methodology and software have enabled the simulation of long-timescale processes that would otherwise not be practical on typical computing resources using standard simulation.

Coarse-graining strategy for molecular pair interactions: A reaction coordinate study for two- and three-dimensional systems

We investigate and provide optimal sets of reaction coordinates for mixed pairs of molecules displaying polar, uniaxial, or spherical symmetry in two and three dimensions. These coordinates are non-redundant, i.e., they implicitly involve the molecules' symmetries. By tabulating pair interactions in these coordinates, resulting tables are thus minimal in length and require a minimal memory space. The intended fields of application are computer simulations of large ensembles of molecules or colloids with rather complex interactions in a fluid or liquid crystalline phase at low densities. Using effective interactions directly in the form of tables can help bridging the time and length scales without introducing errors stemming from any modeling procedure. Finally, we outline an exemplary computational methodology for gaining an effective pair potential in these coordinates, based on the Boltzmann inversion principle, by providing a step-by-step recipe.