Stochastic dynamics of virus capsid formation: direct versus hierarchical self-assembly
© Baschek et al.; licensee BioMed Central Ltd. 2012
Received: 5 July 2012
Accepted: 28 November 2012
Published: 17 December 2012
Skip to main content
© Baschek et al.; licensee BioMed Central Ltd. 2012
Received: 5 July 2012
Accepted: 28 November 2012
Published: 17 December 2012
In order to replicate within their cellular host, many viruses have developed self-assembly strategies for their capsids which are sufficiently robust as to be reconstituted in vitro. Mathematical models for virus self-assembly usually assume that the bonds leading to cluster formation have constant reactivity over the time course of assembly (direct assembly). In some cases, however, binding sites between the capsomers have been reported to be activated during the self-assembly process (hierarchical assembly).
In order to study possible advantages of such hierarchical schemes for icosahedral virus capsid assembly, we use Brownian dynamics simulations of a patchy particle model that allows us to switch binding sites on and off during assembly. For T1 viruses, we implement a hierarchical assembly scheme where inter-capsomer bonds become active only if a complete pentamer has been assembled. We find direct assembly to be favorable for reversible bonds allowing for repeated structural reorganizations, while hierarchical assembly is favorable for strong bonds with small dissociation rate, as this situation is less prone to kinetic trapping. However, at the same time it is more vulnerable to monomer starvation during the final phase. Increasing the number of initial monomers does have only a weak effect on these general features. The differences between the two assembly schemes become more pronounced for more complex virus geometries, as shown here for T3 viruses, which assemble through homogeneous pentamers and heterogeneous hexamers in the hierarchical scheme. In order to complement the simulations for this more complicated case, we introduce a master equation approach that agrees well with the simulation results.
Our analysis shows for which molecular parameters hierarchical assembly schemes can outperform direct ones and suggests that viruses with high bond stability might prefer hierarchical assembly schemes. These insights increase our physical understanding of an essential biological process, with many interesting potential applications in medicine and materials science.
The structure and dynamics of viruses are a fascinating research subject not only from a biological, but also from a physical perspective . In particular, they are a very instructive model system to study self-assembly of large protein complexes with a relatively clear biological function. As viruses do not show metabolic activity of their own, they need to infect host organisms in order to replicate. One key step during the replication process is the formation of the protein shell containing the viral genome. For many viruses, the capsid formation is sufficiently autonomous that it occurs even in vitro. This robustness of the process guarantees successful replication within the dynamic and heterogeneous environment of a living cell. Although virus shell formation is considered as a paradigm for the self-assembly of protein complexes , its underlying principles are far from being fully understood. Progress in our understanding of virus assembly would increase our knowledge of a process of large biological and medical relevance as well as help to advance new self-assembly strategies in materials science applications .
A large variety of mathematical models and simulation approaches has been developed to gain insight into the dynamics of capsid formation from a theoretical perspective. In these approaches the characteristics of protein association and dissociation processes were analyzed depending on parameters like interaction strength, subunit geometry or temperature. The employed techniques range from large-scale Molecular Dynamics (MD) simulations with only a modest amount of coarse-graining of the atomic details  through various schemes of coarse-grained MD [6–11] to patchy particle simulations with interaction potentials [12, 13]. A thermodynamic framework for assembly of icosahedral viruses has been established by Zlotnick and coworkers [14–20]. In general, these studies have revealed that the formation of complete virus capsids requires intermediate bond stability. If interaction strength is too high (or, equivalently, temperature too low), the system becomes kinetically trapped in intermediates which cannot rearrange anymore due to the strong binding. If interaction strength is too low (or, equivalently, temperature too high), the target structure is not sufficiently stable. Another mechanism which can prevent complete capsid formation is the occurrence of misfits, leading to structural polymorphism as often studied with MD-schemes allowing for cluster distortions [21–23].
Due to the large number of single building blocks assembling during virus formation (the simplest icosahedral capsid, T1, has already 60 protein components), there is a multitude of topologically possible assembly pathways. Similar to protein folding, the dominance of few key structures is believed to limit the number of pathways and to speed up the process . In this respect it has been observed that some viruses have developed mechanisms to orchestrate self-assembly by regulating the reactivity of their binding sites [24, 25]. This switching establishes a hierarchy in the formation of transient intermediates during the assembly process. In a number of experiments, partly supported by theoretical calculations, it has been shown that intermediates of pentameric and hexameric symmetry are of special importance for the assembly process of icosahedral viruses [3, 26–32]. Early observations of in vitro assembly of phages and small viruses revealed pentamer sub-structures to play a key role [26, 27]. Experiments on Brome Mosaic Virus , Cowpea Chlorotic Mottle Virus , Human Papillomavirus  and Simian Virus 40 in vivo and in vitro explicitly report capsid assembly from pentameric capsomers. A model for the assembly of Cowpea Chlorotic Mottle Virus suggests that its protein shell assembles from pentamers as well as from trimers of dimers (hexamers) .
Despite the described variety of computational approaches used for virus assembly, to our knowledge the effect of a state-dependent activation of binding sites during the assembly process (hierarchical assembly) has not been explored yet from the theoretical point of view. Although some models consider assembly from pentameric and hexameric clusters, these subunits at the same time represent the smallest entities of the system and their formation from single proteins is not included [10, 21, 33]. Here we investigate the effect of a binding hierarchy on the assembly of icosahedral viruses by comparison of hierarchical and non-hierarchical (direct) assembly from single monomers. We use Brownian Dynamics simulations with reaction patches which have previously been used to study transport-limited protein reactions [34, 35]. Our approach assumes well-defined capsid geometries (in the spirit of local rules) and does not require the use of interaction potentials. This makes our simulations relatively fast, but does not allow us to study structural polymorphism. One particular strength of our approach is that it implements the correct mobility matrix for each possible geometry of the assembling clusters . Another advantage which is exploited here is that one can easily implement hierarchical assembly by an event-driven switching of patch reactivity.
This paper is organized as follows. We first give an overview of the simulation framework and the implemented geometries. Then we present our results for direct and hierarchical assembly of T1 virus capsids. The analysis of T1-assembly is completed with a comparison of the two assembly mechanisms and a discussion of the effect of an increased number of initial monomers. We then explain our results for direct and hierarchical assembly of the more complex T3 virus. They are followed by a detailed analysis of the formation of individual capsomers, which includes a master equation approach. The paper closes with concluding remarks and an outlook to potential future applications of our approach.
The capsid geometries follow the well-established Caspar-Klug scheme where the quasi-equivalent positions in the scaffold of an icosahedral capsid are represented by different types of monomers . The structural complexity is described by the triangulation number T derived from the capsid geometry. T is restricted to certain integer values (T=1,3,4,7,9,...) and denotes the number of protein types which are needed to form a full icosahedral shell. The total number of monomers per capsid is n f =60 T. The icosahedron vertices represent points around which the proteins cluster into close-packed arrays. The proteins grouped around the twelve vertices (which represent axes of fivefold symmetry) form pentamers, while the triangular faces of the icosahedron are covered with hexamers. Every scaffold consists of 12 pentameric and 10·(T−1)hexameric capsomers. In our description we restrict the effect of growing complexity (T > 1) to the hexamers. Thus, every hexamer contains (T−1) different proteins, so that the number of hexameric subunits as well as the number of individual components of each hexamer increase with T. We want to point out that this scheme for dividing virus geometries into ringlike subunits represents only one out of several possible realizations. Following previous approaches to the characterization of icosahedral geometries , we use a set of local rules to define the bond angles between the individual particles. In this way, the resulting structure is encoded in the bond properties of the elementary subunits. Due to the high symmetry of the viral capsid, only a small number of different bonds is sufficient to define a unique target geometry. We note that the exact definition of bond properties impedes the formation of aberrant cluster structures. Therefore the approach used here does not allow us to study structural polymorphism.
Direct assembly is defined as the formation of a capsid from monomers whose bond properties remain unchanged throughout the whole simulation. Thus every reaction patch is active at all times. In contrast to this unconstrained assembly mechanism, hierarchical assembly is decomposed into multiple steps of switching of patch reactivity depending on the configuration of the particles. Since pentameric and hexameric rings have been identified as key subunits in the assembly process, we implement hierarchical assembly as switching of reactivity upon formation of these structures. Initially only the two patches leading to assembly of pentameric or hexameric ring structures are active on every single protein (intra-capsomer bonds). These patches are locked once the ring has closed, so that the formation of these subunits is irreversible. Simultaneously to the locking of intra-capsomer patches the binding sites which connect the pentamers and hexamers with each other (inter-capsomer bonds) are activated so that in a second step, formation of the capsid proceeds via association of the capsomer rings. This collective switching in binding properties should not be confused with the conformational switching of individual subunits which has been used before to study structural polymorphism [21–23].
In this section we investigate the dynamics of virus capsid formation for direct and hierarchical assembly and compare them in order to identify their generic differences. To characterize the assembly performance, the yield (i.e. the relative number of successful trajectories within a given simulation time) and the first passage times (FPTs) of selected intermediates are recorded for different model parameters. We systematically compare both assembly mechanisms in a parameter space ranging from k a =3.0 ns−1 to 9.0 ns−1 and from k d =1.5·10−3 ns−1 to 1.95·10−2 ns−1. For the simulations of assembly of T1 capsids we use an initial monomer concentration of c=4.5 mM (60 particles in a cubic box with side length L=28nm). Investigation of T3 is carried out at an initial concentration of c=1.7 mM (n f =180, L=55nm). To classify different assembly regimes we distinguish between three different phases: the early, intermediate and final phases which we define to be delineated by the emergence of cluster sizes 1/3 n f , 2/3 n f and n f , respectively.
The main difference between the two parameter sets used in Figure 2 is that the second (unfavorable) case leads to more stable intermediates (higher k a , lower k d ). In Figure 2c we systematically investigate the effects of the bond parameters on direct assembly by comparing yield and assembly speed for different combinations of k a and k d . The upper left corner of the parameter plots represent strong bonds (high k a , low k d ), while weak bonds are found in the lower right corner (low k a , high k d ). The left plot shows the relative yield averaged over an ensemble of 40 trajectories. Direct assembly of T1 shows a large region of high yield for dissociation rate values above a threshold of around k d =10.5·10−3 ns−1. Below this value almost no successful assembly is observed. This is due to the limited possibilities of the intermediates to reorganize, which results in the occurrence of kinetically trapped structures. For low dissociation rates we also observe a dependency of the yield on the choice of k a . In this region, lowering of the bond breaking rate k d can at least partly be compensated by lowering of the association rate k a .
In the right plot of Figure 2c, we show the relative assembly speed as a function of model parameters. The assembly speed is defined as the inverse of the completion time of the capsid, v = 1/FPT(n f ). Because this quantity can be obtained only for successful assemblies, here we average only over completed trajectories. In contrast to the relative yield, we see a clear dependence of the assembly speed on the association rate k a across the whole parameter space. Fastest assembly is observed for relatively low values of k a . The observation that relatively high association rates lead to slower assembly can be explained by the increasing tendency to form more than one large cluster in the early and intermediate phases. Thus, even for high dissociation rates, the necessary rearrangement of the clusters slows down the assembly process considerably. We also record a relatively high assembly speed at low k d values where only low yield is observed. Since the relative speed values are obtained by averaging over successful trajectories only, these results show that, if a full capsid is formed, it is completed within a short time.
To conclude, we see that the success of assembly in terms of yield is mostly determined by the choice of the dissociation rate k d . For low values of k d the system becomes kinetically trapped, while large values of k d allow for the reorganization of the clusters. The relative assembly speed of successful trajectories is strongly influenced by the choice of k a . Here we identify an optimum at k a =5.0 ns−1, with speed being worse both at larger and smaller values. In agreement with previous studies, we observe that most efficient assembly (i.e. high yield combined with fast capsid completion) occurs at intermediate bond stability and that bond reversibility is an important requirement for successful capsid formation [9, 40, 41].
The assembly dynamics shown in Figure 3b for an unfavorable parameter combination does not lead to complete assembly within the given simulation time. In contrast to the favorable case (Figure 3a), the association rate is lower and the dissociation rate is higher, which results in a reduced overall bond stability. This is found to strongly hinder the formation of the late pentamers. Although slower, the overall course of the assembly process is not substantially different from the successful case in Figure 3a. The main difference between the two parameter combinations becomes clear by looking at the completion times of the pentamers which are shown in the inset of Figure 3b. We see that the pentamer FPTs of both cases follow the same shape during the early and intermediate phases, but that for low bond stability the completion times in the late phase are delayed. This delay grows with ongoing assembly, so that the final pentamer does not close within the simulation time. Here the negative effect of low monomer concentration on capsomer assembly, which was discussed earlier, is amplified by the low bond stability.
To quantify the effects of different combinations of k a and k d on hierarchical assembly of T1, we again perform a systematic investigation of the bond parameter space as shown in Figure 3c. Considering the relative yield of complete capsids (Figure 3c, left image), we observe that only a narrow range of parameters leads to a considerable fraction of successful trajectories. High yield is only observed at high bond stabilities (high k a , low k d ) in the upper left corner of the parameter plot. To take into account the critical role of the formation of the last pentamer in our simulations, we also show the yield of almost finished capsids (n=55) at the end of the simulation time (Figure 3c, right image). The region where we observe almost finished capsid is considerably expanded and a large fraction of trajectories reaches n=55 in the upper left corner of parameter space. The yield decreases along the diagonal from high towards low bond stability values (lower right corner). It becomes clear that the unfavorable parameter combinations do not show kinetically trapped states as they occur in direct assembly, and that most trajectories are close to capsid completion. The high yield of almost finished capsids and the lack of trapped trajectories suggests that the bond hierarchy promotes successful capsid completion, but is vulnerable to monomer starvation.
For the first emergence of intermediates of half the capsid size (FPT(30), Figure 4a), direct assembly is faster throughout the whole parameter space. This is related to the earlier observation of an extended initial phase of hierarchical assembly when compared to direct assembly (see Figure 3). It can be explained by the fact that the monomers in direct assembly exhibit three active binding sites and thus easily form clusters of considerable size. Since hierarchically assembling monomers are designed to form flat pentamer rings, only the two patches forming intra-capsomer bonds are active until full capsomers are formed. Thus the number of fruitful encounters is reduced remarkably, which leads to the observed slow-down of the initial phase.
Looking at the FPTs for the two-third assembled capsid (FPT(40), Figure 4b), we see a large region in the upper left part of the parameter space (high k a , low k d ) where hierarchical assembly is now able to overtake direct assembly. This can be attributed to two effects. Firstly, hierarchical assembly speeds up once a pool of capsomers is available. Secondly, direct assembly is slowed down at high bond stabilities. Since the combination of fast formation of large, stable clusters in the early phase (due to high k a ) and slow dissociation of small clusters leads to a small number of free monomers, the dominant cluster grows only slowly. In the region of lower bond stability, direct assembly remains faster. Here the increased ability of un- and rebinding of single proteins allows for fast rearrangement, leading to a sufficiently large supply of free monomers so that the dominant cluster can easily grow beyond n=40. Simultaneously the pentamer rings in the hierarchical setup form slower than at high bond stabilities.
The difference between the two assembly mechanisms becomes even more evident when looking at FPT(50) (Figure 4c). At low k d values, direct assembly experiences kinetic trapping. As a consequence, hierarchical assembly is superior for almost all small k d values, also at points where the question of dominance remained undecided for FPT(40). The parameter region of weak bonds where direct assembly is faster than hierarchical one is observed to extend during the step from FPT(40) to FPT(50) (lower right corner of parameter space). Under these conditions the effect of beginning monomer starvation delays the pentamer completion of hierarchical assembly.
For the assembly speed of the complete virus capsid (FPT(60), Figure 4d), direct assembly dominates again across almost the whole parameter space. Only at very high bond stabilities hierarchical assembly shows lower or comparable FPT values. This is not surprising taking into account the results for the overall yield of hierarchical assembly (Figure 3c) and underlines the large impact of monomer starvation on hierarchical assembly.
We conclude that hierarchical assembly is not always better than direct assembly. Direct assembly performs better both in the initial and final phases. During the intermediate phase, however, hierarchical assembly is more successful, because it does not suffer from stable bonds preventing structural rearrangements. Due to the limited number of possible interactions, hierarchical assembly is unlikely to get trapped in sub-pentameric units. In general, for hierarchical assembly parameter combinations resulting in high bond stability are favorable. At these values we observe kinetic trapping of most of the directly assembling systems. In addition, the symmetry of the pentamers themselves and the low complexity of their interactions prevent them from getting trapped in large clusters incompatible with the final capsid. For T1, this favors the step-wise build-up of the target structure.
Figures 5d-5f show the yield of the first capsid within simulation time for an initial number of N=60, N=80 and N=120 monomers, respectively. Here again red indicates a higher yield of hierarchical assembly while blue indicates a higher yield of direct assembly. Gray marks parameter pairs with the same yield. In contrast to the FPT(50), we can clearly see that increasing the initial number of monomers results in a largely expanded parameter space in which hierarchical assembly is favorable. This shows that monomer starvation affects the final phase of hierarchical assembly in particular as it has been inferred in the previous section. In fact hierarchical assembly performs well throughout the whole parameter space and shows high yield for intermediate and weak bonds. At very high bond strength we even observe some trapping for hierarchical assembly. However, direct assembly still strongly suffers from kinetic trapping so that the parameter space corresponding to high bond strength remains clearly dominated by hierarchical assembly. We also note that increasing the initial number of monomers from N=80 to N=120 does not lead to a further promotion of hierarchical assembly, presumably because now two capsids form in parallel, each drawing monomers in a similar manner as before for N=60.
To conclude, we find that our main results from the previous section remain valid for an increased number of initial monomers. Hierarchical assembly is favorable at high bond strength due to the decreased possibility of trapping while direct assembly is favorable at low bond strength allowing for fast reorganization of large clusters. In general, we expect that our results also carry over to even larger systems.
In order to further investigate the effects that slow down hierarchical assembly, we now analyze the dynamics of hexamer and pentamer formation both with computer simulations and a master equation approach. To compare the FPTs for pentamer and hexamer formation, we scale these values with the number of monomers per capsomer ring. This linear scaling is based on the assumption that the mean time for a net addition of monomers to small ring-forming clusters is independent of the cluster size. This simplification in particular neglects the increased number of decay paths of hexamers compared to pentamers. However, the assumption seems justified for the present case of high bond stabilities (high k a , low k d ), at least for the early and intermediate phase of assembly.
Numerical evaluation with the initial condition ν 1(t=0)=Ngives the time evolution of all cluster size populations ν n . By fitting the set of equations to the course of all ν n (t) from the complementary simulation (n f =6, N=60), we obtain parameter combinations a, b which reproduce the observed assembly dynamics. Under the constraint that the dissociation rate b per bond is constant for identical and T3-like hexamers (since the simulations apply the same k d ), we find the following parameters: a id =2.4·10−6 ns−1 (identical hexamers), a T3=1.2·10−6 ns−1 (T3-like hexamers) and b=9·10−5 ns−1. In general, all ν n (t)are reproduced well. This suggests that the assumption of a constant association rate a per bond, independent of the sizes of the encountering clusters, is a reasonable approximation for the formation of small rings. The results for ν 1(t) and ν 6(t) are displayed in Figure 9b together with the simulation data points for both types of hexamers. The early phase is the region which exhibits the largest discrepancies between data and ME results, while the final phase of assembly shows a high level of consistency. This can be explained by the fact that the rate equation framework does not include any spatial constraints and is thus not able to reproduce the same sort of lag time before the first protein reactions as was observed in the simulations, where the randomly distributed particles react only after diffusional mixing leads to the first encounter events. This is also the reason why the difference between the cases of T3-like and identical hexamers becomes visible in the simulation data only after a certain time, while the ME results differ from the very first iteration step (see Figure 9b). Since the rate equations do not contain a diffusional component, the coefficients a and b cannot be directly related to the simulation parameters k a and k d . While k a determines the rate of transition from encounter to a bound state, a as well includes the formation of diffusional patch overlap. We estimate a=k a /(N A ·V), where V is the simulation box volume. Using the initial concentration c=N/(N A ·V), we find the expression . a is thus, as expected, proportional to the initial monomer concentration in the simulation box. Applying this relation to the fit parameters using the effective initial concentration of the protein types, we estimate the overall association rate values to be and . The fact that the association rate a id for identical hexamers is about twice the value found for a T3confirms that the difference in assembly dynamics for identical and T3-like hexamers has its origin in a reduced association rate, caused by reduced encounter of matching protein types. Our observations suggest that the association rate decreases linearly with increasing number of bond partners in the system and thus the number of different protein types needed to form a capsomer ring. When comparing our values for the diffusional encounter rate to data from experiments, we see that we overestimate the association rate. In general, the association rate for bimolecular binding reactions is experimentally found to lie between 4·106 and 107 s−1 M−1. Absence of long-ranged forces, as it is the case for our simulation framework, is predicted to push the rates below 106 s−1 M−1. The reason for our relatively high estimates for the encounter rate could be the treatment of dissociation as a stochastic event without immediate relocation of the partners. In the present implementation, two patches stay in an encounter after dissociation and their movement is subject to the cluster mobility. We assume this to cause an overestimation of rebinding frequencies which results in an increased association constant. Whereas the association rate constant can be related to other results, there is no such argument for the value of b.
Understanding the biophysical principles underlying the self-assembly of virus capsids is of fundamental importance for biology and medicine, and might also promote novel applications in materials science. Here we have presented a Brownian dynamics study of the assembly of icosahedral virus capsids. Using a patchy particle model without potentials, our simulations are relatively fast and therefore we are able to obtain good statistics with relatively modest computing times. One special strength of our approach is the rigorous treatment of translational and rotational diffusion, with the motility matrices for any cluster shape calculated on the fly. Our approach is particularly suited to focus on the effect of a bonding hierarchy on the performance of the assembly process. The hierarchy was established by an event-driven switching of bond characteristics upon the formation of capsomer rings, which have earlier been identified as key intermediate structures of the assembly pathway of some icosahedral viruses [3, 26–28, 30–32]. We first conducted a detailed comparison of direct versus hierarchical assembly for T1 viruses. To elucidate the effects of an increased complexity of the capsid geometry on the formation of the capsomer rings, we then performed a detailed analysis of capsomer assembly for T3 viruses, including a master equation approach complementing the computer simulations.
Our results for direct assembly of T1 virus capsids show that capsid completion is only successful if the bonds are weak enough to allow for a sufficient number of unbinding and reorganization events. Otherwise kinetically trapped clusters appear. These findings are in good agreement with the results of previous approaches [6, 9, 40, 41]. In marked contrast, hierarchical assembly performs better for high bond stabilities, as the imposed hierarchy reduces kinetic trapping. However, hierarchical assembly is more vulnerable to monomer starvation in the final phase. This effect has previously been observed in other approaches for direct assembly [6, 43], but it is even more severe for hierarchical assembly specifically studied here. Comparison of direct and hierarchical assembly reveals that hierarchical assembly, although slower in the early phases, is able to outperform direct assembly at high and intermediate bond strength. This is due to the fact that capsids assembling from highly symmetric capsomers do not require fundamental reorganizations to achieve large cluster size, as it is the case in direct assembly.
The analysis of T3 virus assembly shows that the effects apparent for T1 viruses become amplified by the increased complexity of the capsid geometry. In general, the assembly process of T3 viruses is slower due to the size of the capsid and the increased complexity of the protein interactions. Starting with exactly 180 monomers the parameter space for successful direct assembly is narrowed and we do not observe any complete capsids in hierarchical assembly within the used simulation time. Investigation of the course of assembly of the two mechanisms reveals that they both perform best in distinct regions of the parameter space. Increasing the initial number of monomers we find that hierarchical assembly performs better while direct assembly remains widely unaffected. However, we still observe that both mechanisms are favorable in distinct regions of the parameter space. To analyze the effect of geometric complexity on capsomer formation during hierarchical assembly, we perform a closer analysis of assembly of different capsomer types. The results show a significant slow-down of capsomer formation with increasing structural complexity, which explains why we do not observe any full T3 virus capsids in the hierarchical setup within the given simulation time. These findings suggest a further slow-down for the assembly dynamics of more complex capsids such as T4 and T7 for hierarchical assembly.
Computer simulations of virus assembly are usually carried out with a fixed number of initial monomers and therefore necessarily lead to monomer starvation in the final phase. In vivo, this constraint should be less relevant than in our simulations. Once a cell is infected by a virus, one expects to see a constant production rate for viral proteins, and therefore monomer starvation should be less of an issue. It would be interesting to test if in such a situation, hierarchical assembly becomes even more favorable than found here. In computer simulations, this could be done by continuously adding new monomers and removing completed capsids. We leave this to future studies as it would entail to introduce at least two more model parameters, namely the rates for monomer injection and capsid removal, as well as explicit rules on the spatial positioning of the new monomers. A similar study could be done experimentally for viruses which self-assemble in vitro, although here too there might be technical problems to implement such procedures. For in vivo systems, such studies would depend very much on the details of the virus assembly of interest, in particular on the spatial coordination in regard to the different cellular compartments.
Our simulation framework has great potential for further investigation of assembly of icosahedral viruses, i.e. capsids with higher T-number (T4, T7,...). Although larger simulation times become necessary, they are potentially much smaller than the ones required for less coarse-grained approaches, including patchy particle models with interaction potentials or coarse-grained molecular dynamics simulations. A particular strength of our approach is the possibility to switch patch reactivity during the assembly process. This suggests to investigate even more complex ways to build virus capsids. Our approach could also be applied to other interesting cases of protein assembly, for example to the actin cytoskeleton, for which different regulatory proteins lead to changes in local reactivity.
We conclude that it might be beneficial for icosahedral viruses to assemble hierarchically, since it prevents kinetic trapping and allows for faster formation of larger structures. Our results suggest that hierarchical assembly performs better than direct assembly for high and intermediate bond stability, while direct assembly is favorable for weak bonds allowing for fast reorganization. For complex viruses, our study suggests that the problem of monomer starvation and critical concentrations has to be addressed for each type of monomer separately, thus making the process more vulnerable for fluctuations in the supply chain and imposing limits to the overall degree of complexity. The partitioning of parameter space into favorable regions for direct versus hierarchical schemes becomes even stronger for more complex capsid geometries and suggests ways to design optimal assembly schemes for different molecular species.
We thank Jakob Schluttig for providing the initial computer code and for helpful discussions. HCRK was supported by fellowships from the Landesgraduiertenstiftung Baden-Wuerttemberg and from the Cusanuswerk. USS is member of the CellNetworks cluster of excellence at Heidelberg.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.