This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.105.076737v1 90/9/3029    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Hemberg, M. Articles by Barahona, M. Search for Related Content PubMed PubMed Citation Articles by Hemberg, M. Articles by Barahona, M.

Stochastic Kinetics of Viral Capsid Assembly Based on Detailed Protein Structures

Martin Hemberg ^*, Sophia N. Yaliraki  and Mauricio Barahona ^*

^* Department of Bioengineering and Department of Chemistry, Imperial College London, London, United Kingdom

Correspondence: Address reprint requests to M. Barahona, Dept. of Bioengineering, Imperial College London, South Kensington Campus, London, SW7 2AZ United Kingdom. E-mail: m.barahona{at}imperial.ac.uk.

We present a generic computational framework for the simulation of viral capsid assembly which is quantitative and specific. Starting from PDB files containing atomic coordinates, the algorithm builds a coarse-grained description of protein oligomers based on graph rigidity. These reduced protein descriptions are used in an extended Gillespie algorithm to investigate the stochastic kinetics of the assembly process. The association rates are obtained from a diffusive Smoluchowski equation for rapid coagulation, modified to account for water shielding and protein structure. The dissociation rates are derived by interpreting the splitting of oligomers as a process of graph partitioning akin to the escape from a multidimensional well. This modular framework is quantitative yet computationally tractable, with a small number of physically motivated parameters. The methodology is illustrated using two different viruses which are shown to follow quantitatively different assembly pathways. We also show how in this model the quasi-stationary kinetics of assembly can be described as a Markovian cascading process, in which only a few intermediates and a small proportion of pathways are present. The observed pathways and intermediates can be related a posteriori to structural and energetic properties of the capsid oligomers.

This article has been cited by other articles:

				H. Gohlke and M. F. Thorpe A Natural Coarse Graining for Simulating Large Biomolecular Motion Biophys. J., September 15, 2006; 91(6): 2115 - 2120. [Abstract] [Full Text] [PDF]

				M. F. Hagan and D. Chandler Dynamic Pathways for Viral Capsid Assembly Biophys. J., July 1, 2006; 91(1): 42 - 54. [Abstract] [Full Text] [PDF]