BRIEF RESEARCH REPORT article

Front. Mol. Biosci., 20 October 2022

Sec. Structural Biology, Biophysics and Evolution

Volume 9 - 2022 | https://doi.org/10.3389/fmolb.2022.967877

The diversity of protein-protein interaction interfaces within T=3 icosahedral viral capsids

  • Simons Centre for the Study of Living Machines, National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bengaluru, India

Abstract

Some non-enveloped virus capsids assemble from multiple copies of a single type of coat-protein (CP). The comparative energetics of the diverse CP-CP interfaces present in such capsids likely govern virus assembly-disassembly mechanisms. The T = 3 icosahedral capsids comprise 180 CP copies arranged about two-, three-, five- and six-fold axes of (quasi-)rotation symmetry. Structurally diverse CPs can assemble into T = 3 capsids. Specifically, the Leviviridae CPs are structurally distinct from the Bromoviridae, Tombusviridae and Tymoviridae CPs which fold into the classic “jelly-roll” fold. However, capsids from across the four families are known to disassemble into dimers. To understand whether the overall symmetry of the capsid or the structural details of the CP determine virus assembly-disassembly mechanisms, we analyze the different CP-CP interfaces that occur in the four virus families. Previous work studied protein homodimer interfaces using interface size (relative to the monomer) and hydrophobicity. Here, we analyze all CP-CP interfaces using these two parameters and find that the dimerization interface (present between two CPs congruent through a two-fold axis of rotation) has a larger relative size in the Leviviridae than in the other viruses. The relative sizes of the other Leviviridae interfaces and all the jelly-roll interfaces are similar. However, the dimerization interfaces across families have slightly higher hydrophobicity, potentially making them stronger than other interfaces. Finally, although the CP-monomers of the jelly-roll viruses are structurally similar, differences in their dimerization interfaces leads to varied dimer flexibility. Overall, differences in CP-structures may induce different modes of swelling and assembly-disassembly in the T = 3 viruses.

Introduction

Non-enveloped viruses enclose their genome in a shell composed exclusively of proteins. In icosahedral viruses, this protein shell, also known as a capsid, has icosahedral symmetry (Louten, 2016). The capsids of some icosahedral viruses are assembled using a single type of protein, called the coat-protein (CP). Icosahedral viruses are further classified based on their triangulation (T) number, which, for simple viruses, is generally equal to the number of unique CP conformers that are present in the unit that is used to tile the icosahedron (). Thus, a T = 3 icosahedral capsid, composed of 180 CPs (Figure 1), has a CP trimer in each tiling unit (also known as the icosahedral asymmetric unit: IAU). Other distinct CP-CP interactions give rise to (quasi-)symmetric dimers, pentamers and hexamers (Figures 1A,B). Henceforth, the interface between the (quasi-)symmetric dimers (CP2) is called the dimerization interface, the interface between any two CPs from an IAU is called a trimerization interface and so on (Figures 1D–G).

FIGURE 1

The process of virus assembly is not completely understood even in viruses whose capsids are made of a single type of CP. Although proteins of the infected-host (HP) may affect virus assembly-disassembly, multiple experimental studies have used capsids expressed from bacterial expression systems, assembled in the absence of HP-CP interactions (Zhao et al., 1995; Sastri et al., 1999; Lokesh et al., 2002; Powell et al., 2012), or have studied disassembly with purified virus particles (Pappachan et al., 2009; Powell et al., 2012; ). In such assembly-disassembly processes, the balance of energetic interactions among CPs and between the CPs and the genome will determine the order of events in assembly. Capsids can also assemble around non-genome cargo (Ren et al., 2006; Sokullu et al., 2019; ) indicating that CP-CP interactions can by themselves direct capsid assembly. In order to understand potential principles that govern such CP-CP interface energetics, we analyzed the strength and nature of all CP-CP interfaces (Figures 1D–G) in diverse T = 3 icosahedral viruses.

The Leviviridae, Bromoviridae, Tombusviridae and Tymoviridae families have positive-sense single-stranded RNA genomes enclosed in T = 3 icosahedral capsids (). However, while the bacteriophage Leviviridae CPs fold to an α+β fold with two α-helices and seven β-strands (Fold d.85 in SCOPe (; ); Figure 2A), the Bromoviridae, Tombusviridae and Tymoviridae are mostly plant viruses whose CPs fold to the jelly-roll fold (SCOPe classification: b.121.4; Figures 2B–D). The jelly-roll fold is composed of two stacked four-stranded antiparallel β-sheets and is present in diverse virus CPs (Rossmann and Johnson, 1989; ; ). Despite the structural differences between the CPs, experimental disassembly data is available for at least one virus from each family which indicates that the viruses disassemble into dimers (; Stockley et al., 2007; Pappachan et al., 2009; Powell et al., 2012; ). The question that arises then is, does the overall capsid symmetry rather than interface energetics drive disassembly processes.

FIGURE 2

In the commonly studied Levivirus bacteriophage MS2, it is known that the capsid assembly unit is not CP but CP2, the homodimer composed of the two symmetry related CPs interacting at the dimerization interface (Ni et al., 1995; Stockley et al., 2007). Our earlier folding simulations (Prakash and Gosavi, 2021) showed that MS2-CP2 is an obligate dimer, i.e., it cannot fold into two independent fully folded monomers. MS2-CP2 also folds in a two-state manner, populating only the unfolded and dimeric folded ensembles without populating either a dimeric or a two-monomer intermediate ensemble. Using folding simulations, a previous study (Levy et al., 2004b) classified diverse homodimers into either two-state folding dimers or three-state folding dimers populating a partially folded ensemble. The authors argued that folding outcomes are probably determined by the nature of the dimeric interface and constructed two interface dependent parameters ad hoc. The first parameter (R), the ratio of the number of inter-monomer contacts to the intra-monomer contacts, is a proxy for the relative size of the interface. The second parameter (H), the average hydrophobicity of the residues in the interface, is a proxy for how buried (and not solvent exposed) the interface residues want to be. It was found that H and R could also classify homodimers into two-state (high H and R values) and three-state folders (lower H and R values).

We previously calculated the H and R for MS2-CP2 and these values predicted two-state folding in agreement with simulations (Prakash and Gosavi, 2021). Here, we hypothesize that H and R can also be more generally used to understand how similar CP-CP interaction interfaces are to each other. Accordingly, we calculated H and R for all the interfaces (dimerization, trimerization, pentamerization and hexamerization) of viruses from the four T = 3 families. We find that the dimerization interfaces of the Leviviridae have higher H and R values than all other interfaces indicating that these CP2 dimers are predicted to be two-state folders. The H values of the jelly-roll fold dimerization interfaces are generally marginally higher indicating that the interface residues may be buried more strongly through dimerization possibly enabling disassembly into CP2 dimers.

To further investigate the CP2 dimers, we performed molecular dynamics (MD) simulations of coarse-grained structure-based models () of a representative CP2 from the four virus families. The inter-monomer distance changes the least in MS2-CP2 indicating a rigid interface while this distance varies the most in the Bromoviridae representative. Overall, differences in the structure of the dimerization interface give rise to distinct interface dynamics in the three CP2 jelly-roll fold dimers and this is likely to lead to different modes of assembly and disassembly.

Methods

Obtaining capsid structures for interface analysis

PDB IDs for T = 3 capsids from Leviviridae, Bromoviridae, Tymoviridae and Tombusviridae, were obtained from VIPERdb (http://viperdb.scripps.edu) (Montiel-Garcia et al., 2021) and filtered to eliminate redundant structures. Assembled capsid structures for the selected PDB IDs were then downloaded from RCSB PDB (). From the complete capsid structures, a tri-dimer structure was extracted (Figure 2), which included one triangular IAU with A, B and C chains and one chain each extracted from three adjacent IAUs, such that these chains had a dimerization interface with the chains from the original IAU. This tri-dimer is the minimal unit containing all interactions between adjacent CPs in the T = 3 capsids (Figures 1C–G).

Calculating normalized interface size and hydrophobicity

The H and R parameters (Figure 3) were calculated as done previously for diverse homodimers (; Levy et al., 2004b). When two residues are “close” in structure, they are said to be in contact. Here, interatomic intra-CP and inter-CP contacts were calculated using the CSU software (Sobolev et al., 1999) for every pair of adjacent CPs in the tri-dimer structure. These contacts were then projected onto their corresponding residues. For CPs labeled x and y, the number of intra-CP contacts are Nintra(x) and Nintra(y) respectively, while the inter-CP contacts are Ninter(xy). The ratio of the number of inter-CP to the intra-CP contacts is then

FIGURE 3

Inter-CP interface hydrophobicity was calculated by averaging over the hydrophobicity factor, hi (Pacios, 2001), for every residue i, which participated in an inter-CP contact. Thus,where n = 2×Ninter(xy). It should be noted that if a residue participates in two different inter-CP contacts, then it is counted twice. hi (and H) ranges from 0 to 1, with hi = 1 for the most hydrophobic residue and hi = 0 for the least hydrophobic residue. We calculate the hydrophobicity factors using a scale that was used in the original H parameter calculations for homodimers (Levy et al., 2004b).

Multiple structure alignment of jelly-roll fold CPs

Two ‘C’ CP monomers connected by a two-fold rotation symmetry (CC-CP2) were extracted from the tri-dimer structures of representative viruses (Figure 2): Bromoviridae: Cowpea Chlorotic Mottle Virus: CCMV (PDB ID: 1CWP); Tombusviridae: Tobacco Necrosis Virus: TNV (PDB ID: 1C8N); and Tymoviridae: Turnip Yellow Mosaic Virus: TYMV (PDB ID: 1AUY). For one monomer from each of the CP2s, we performed a multiple structure alignment using the STAMP algorithm (Russell and Barton, 1992) as implemented in the Multiseq extension (Roberts et al., 2006) of the VMD package (). The aligned CPs (Figures 4B–D left chains) were colored based on their alignment from blue through white to red, with well-aligned regions in blue and unaligned regions in red. The second chains (Figures 4B–D right chains) of the dimers were not aligned. These jelly-roll CP2s and the MS2-CP2 (also extracted from the bacteriophage MS2 tri-dimer; PDB ID: 2MS2, Figure 4A) were used to understand dimer dynamics.

FIGURE 4

) and are shown in the aligned position. They are colored from blue through white to red based on how good the alignment is with blue indicating a good alignment and red indicating no alignment. The core β-sheet jelly-roll fold (pale blue to white) aligns well across the three viruses. The change in orientation in the right monomer (orange, blue and green) is due to the distinct dimer interfaces in the three viruses. (E) Histograms of the COM distance calculated from equal time coarse-grained simulations of the dimers. The vertical dashed line (colored like the histogram) shows the COM distance in the native folded structure. MS2-CP2: mean COM distance = 1.67 nm, standard deviation = 0.03 nm. TNV-CP2: mean = 4.42 nm, standard deviation = 0.08 nm. TYMV-CP2: mean = 4.43 nm, standard deviation = 0.06 nm. CCMV-CP2: mean = 3.35 nm, standard deviation = 0.23 nm. A broader histogram such as that seen in CCMV indicates a more flexible dimer interface. For each of these CP2 dimers, the alignment of the native structure with sampled structures whose COM distances are the mean, and mean ± 2×(standard deviation) are shown in Supplementary Figure S4.

Simulating CP2 dynamics

We used a structure-based model coarse-grained to a single Cα-bead per residue (Cα-SBM) () to simulate one CC-CP2 dimer (see previous section) of a representative virus from each family. This Cα-SBM was previously used to simulate the folding of MS2-CP2 and its variants using MD simulations (Prakash and Gosavi, 2021). Details of the Cα-SBM and simulations are described in the supplementary information. The distance between the centers of mass (COM) of the two chains in the dimer was calculated for structures sampled during a simulation. These COM distance values were then binned (bin width 0.01 nm) into a histogram (distribution) which gives the normalized number of structures present in the bin with a given COM distance value.

Results and discussion

The Leviviridae CPs are obligate dimers

As stated in the introduction, two parameters, R, the ratio of interface to monomer contacts, and H, the average interface hydrophobicity, were previously shown to be able to classify a set of diverse homodimers as being two-state or three-state folding in simulations (Levy et al., 2004b). Since R and H are physically intuitive interface properties with R being a proxy for interface size and H for its “stickiness”, we use these parameters more generally to understand the strength and nature of CP-CP interfaces.

We first calculated the H and R values of the set of homodimers used previously (Levy et al., 2004b) and drew an arbitrary line in the H-R space such that homodimers which lie on one side of the line are two-state folders while those on the other side are three state folders (Figure 3). However, it should be underlined that this division is arbitrary and there are insufficient data points to demarcate the boundary exactly. So, homodimers that lie near the boundary-line, on either side, could either be two-state with a small folding barrier or three-state with a weakly populated intermediate. Since we intend to use H and R more generally here, we assume that homodimers that lie in the two-state region have stronger interfaces than homodimers that lie near the boundary which in turn have stronger interfaces than homodimers that lie in the three-state folding region.

Our previous folding simulations of the Levivirus MS2-CP2 agreed with the H and R values and showed that MS2-CP2 was a two-state dimer indicating that MS2-CP2 has an extensive interface which couples the folding of the two CPs (Prakash and Gosavi, 2021). Other Leviviridae CP2 also have high H and R values (Figure 3A) indicating that CP2 is the capsid assembly unit in the entire family. In contrast, all other Leviviridae CP-CP interfaces (the trimerization, pentamerization and hexamerization interfaces) have lower R values. Further, trimerization interfaces are the least hydrophobic (Figure 3A).

An examination of MS2-CP2 shows a two-layer interface (Figure 4A). The inner layer is formed by the anti-parallel β-β hydrogen bonding interactions between the edge β-strands (named βG in MS2) of the two extended β-sheets, one from each CP. The outer layer comprises the C-terminal α-helices of each CP stacked over the β-sheet of the partner CP. The CP helices are also sandwiched between the partner CP’s N-terminal β-hairpin and C-terminal α-helices, creating an interdigitated topology (Figure 4A). This two-layer structure gives rise to an extensive interface area with R varying between 0.95–1.50. Moreover, 50–70% residues participating in interface contacts are non-polar, giving rise to large H values of 0.35–0.45 (Supplementary Table S1). Previously (Prakash and Gosavi, 2021), we had computationally designed two MS2-CP2 variants, each having reduced interactions in one of the interface layers. One variant was a pseudo-circular permutant (CiP) CP2 which had the C-terminal helices from each monomer excised and linked to the N-terminus of the partner monomer (Supplementary Figure S1A). This preserved all dimer contacts but converted inter-CP β-sheet-helix contacts into intra-CP contacts. The H and R values of CiP-CP2 lie close to the line dividing the two-state and three-state folders (Supplementary Figure S1B). In agreement, CiP-CP2 is a two-state folder with a barrier much smaller than that of wild type (WT) MS2-CP2. Further unlike WT-MS2-CP, CiP-CP monomers can fold completely (Prakash and Gosavi, 2021). A similar effect (two-state folding with a low barrier and high H values with lower R values) was seen in the other variant of MS2-CP2 (termed MS2-ΔβG:βG-CP2) in which the interactions between the two antiparallel βG edge strands were deleted.

In summary, for each of the CP-CP interfaces (dimerization, trimerization, pentamerization and hexamerization), the H and R values are similar across Leviviridae. Additionally, the Leviviridae have intertwined dimerization interfaces which lead to higher R values and two-state folding of MS2-CP2. We next analyze the CP-CP interfaces of the other three families and compare these with those of MS2-CP2.

The CP2 interfaces of the jelly-roll fold viruses are similar to their other interfaces

The sizes of the dimerization interfaces (CP2; Figures 2B–D) of the CPs from the three jelly-roll virus families are similar to their corresponding pentamerization and hexamerization interfaces with R ranging from 0.07 to 0.50 (Supplementary Tables S2–S4). However, except in Bromoviridae (Figure 3B, Supplementary Table S2), the dimerization interface is generally marginally more hydrophobic than the pentamerization and hexamerization interfaces (Figures 3C,D, Supplementary Tables S3, S4). The trimerization interface in Bromoviridae and Tombusviridae is the least hydrophobic (Figures 3B,C) while in Tymoviridae, the trimerization interface is closer in hydrophobicity to the pentamerization and hexamerization interfaces (Figure 3D, Supplementary Table S4). Overall the H and R values of all the CP-CP interfaces of the jelly-roll viruses are also similar to the non-CP2 (non-dimerization) interfaces of the Leviviridae.

What then could be the reason for CP2 either forming first or persisting after capsid disassembly in disassembly/reassembly experiments with jelly-roll fold viruses (; Zhao et al., 1995; Sastri et al., 1999; Zlotnick et al., 2000; Tang et al., 2006; Powell et al., 2012) as well as in experiments involving assembly defective mutants (Tang et al., 2006; Pappachan et al., 2009)? We find that although the H and R values for all the CP-CP interfaces of the jelly-roll fold families lie in the three-state folding region, the dimerization interfaces lie close to the arbitrary line (dashed-line in Figures 3A–D) separating this region from the two-state folders. Thus, the interface of CP2 is likely to be stronger than interfaces that lie in the three-state folding region and is likely to remain intact in conditions where the other interfaces break, leading to early formation of and late dissociation of CP2 dimers in assembly/disassembly experiments. Consequently, we focus on the CP2 interfaces and examine them in detail.

Despite having structurally-similar monomers, the dimerization interfaces of the jelly-roll fold viruses are distinct

Since the three jelly-roll fold virus families have structurally similar CPs, we performed a multiple structure alignment on one monomer from the CP2 of one representative virus from each of the three families (see Methods; Figures 4B–D). The core β-sheet region forming the jelly-roll fold is fairly well-aligned in all three virus CPs, but each CP has additional secondary structure elements and a completely different dimerization interface (Figures 4B–D). The TNV and TYMV CP2s have an interface involving side-by-side interactions between the edge β-strands of the two CPs (Figures 4C,D), reminiscent of the antiparallel β-sheet inner layer of MS2-CP2 (Figure 4A). In addition to these β-β interactions, the TNV-CP2 interface also contains small α-helices from each monomer stacked over its own β-sheet that are also packed against each other (Figure 4C). This makes the interface similar to that of MS2-CiP-CP2 (Supplementary Figure S1A) and as expected, their H and R values are also similar (Supplementary Figure S1B).

The CCMV-CP2 interface is distinct from the TYMV and TNV CP2 interfaces and lacks β-β interactions. Each monomer has a C-terminal tail which reaches out and interacts with the β-sheet of the core jelly-roll fold of the partner monomer (Figure 4B). This is similar to the variant, MS2-ΔβG:βG-CP2, in which the interchain β-β contacts were deleted (Prakash and Gosavi, 2021). The role of the swapped C-terminal helices in MS2-ΔβG:βG-CP2 is played by the swapped C-terminal tails in CCMV. The lack of β-β contacts and the length of the C-terminal tails creates a space or a cavity between the two CPs of the CCMV-CP2 (Supplementary Figure S2B) and this gap should make the CCMV-CP2 interface (and that of other Bromoviridae) more flexible.

To test this hypothesis, we simulated the CP2s from MS2, CCMV, TNV and TYMV using Cα-SBMs and calculated the distributions of the distance between the centers of mass of the two monomers (COM distance) from the simulation trajectory (see Methods). As expected, the two-layer MS2-CP2 interface has the most sharply peaked distribution (Figure 4E, Supplementary Figure S3) with the least standard deviation from the mean COM distance. The stiffness of this interface is pictorially depicted in Supplementary Figure S4A. TNV-CP2 and TYMV-CP2 have distributions which are more similar to each other (Figure 4E, Supplementary Figures S4C,D) than they are to CCMV-CP2. CCMV-CP2 (Figure 4E) has a broad COM distance distribution with the highest standard deviation. This flexibility in the CCMV-CP2 dimerization interface (Supplementary Figure S4B) has been shown experimentally and computationally (Tama and Brooks, 2002; ) to allow capsids with varying sizes. Capsids from Bromoviridae have also been shown to swell at neutral pH = 7.0 and at low divalent metal ion (Ca2+ or Mg2+) concentrations (Speir et al., 1995). The modelled pseudoatomic structure of the swollen capsid (Tama and Brooks, 2002) retains the pentamerization and hexamerization interfaces but a change in conformation at the dimerization interface expands the pore at the trimerization axis (QC3, Figure 2). Thus, Bromoviridae CP2 flexibility may be responsible for the ability of the capsid to swell. An additional role of the capsid flexibly could be to allow it to package genomes of slightly varying sizes. Unlike Tymoviridae and Tombusviridae, viruses belonging to Bromoviridae have their genome divided in four parts that are packaged separately in three different capsids (tripartite genomes) (; ). It was also shown that the C-terminal tail which forms the dimerization interface in Bromoviridae may be important for cell-to-cell movement of the viral capsids (Okinaka et al., 2001). In fact, Bromoviridae reassembly studies have shown aggregates other than regular T = 3 or T = 1 capsids (; ; Tang et al., 2006).

Capsid asssembly and disassembly for understanding infections and in protein engineering

Due to the size and strength of the Leviviridae dimerization interfaces, they are expected to be stronger than the trimerization, hexamerization and pentamerization interfaces, implying that during disassembly, CP2-CP2 interactions will break before interactions within a CP2 dimer (intra-CP2), as was shown in previous acid disassembly experiments (Stockley et al., 2007). In contrast, because of the relative chemical similarity between the dimerization, pentamerization and hexamerization interfaces in jelly-roll viruses, these interactions are likely to break at similar denaturant conditions, while the trimerization interfaces may break in milder conditions (Szoverfi and Fejer, 2022). The swelling of the CCMV capsid by increasing the size of the pore at the trimerization axis (Speir et al., 1995; Tama and Brooks, 2002) is in accordance with this observation. A previous computational study (Polles et al., 2013) determining the rigid mechanical units in virus capsids, showed that while CP2 was the mechanically stable unit in MS2, pentamers and hexamers were stable in CCMV. Although the pentamers and hexamers are mechanically rigid units, the swapped C-terminal tails in CCMV-CP2 link hexamers to both adjacent hexamers and pentamers potentially buffering the breaking apart of viruses during disassembly.

The difference in the protein interfaces may also enable different infection mechanisms. Leviviridae infection is believed to involve disassociation of the maturation protein (which replaces one CC-CP2 in the WT virus) along with the genomic RNA from the capsid (Meng et al., 2019). In contrast, for Bromoviridae, infection is hypothesized to follow co-translational disassembly, where the genomic RNA is exposed for translation due to swelling of the capsid at the neutral pH of the host cell (Roenhorst et al., 1989; Speir et al., 1995; Tama and Brooks, 2002; ). Interestingly, similar swelling is also observed in both Tombusviridae (; Llauró et al., 2015) and Tymoviridae (Virudachalam et al., 1985; Sastri et al., 1999), which do not have flexible CP2 interfaces. Thus, although it is assumed that they too co-translationally disassemble like Bromoviridae, the rates of swelling and infection could be very different.

Overall, even coarse-grained parameters and models which do not contain atomic details of the interfaces can be used to both understand and propose hypotheses in assembly-disassembly mechanisms especially when such mechanisms are studied in vitro. However, it is important to note that the present coarse-grained approaches do not include metal ions, virus RNA or host proteins and RNA. For instance, ordered genomic RNA such as is present in Leviviridae (; Rolfsson et al., 2016; ) can also modulate the mechanism of viral assembly by coding specific sites for RNA-protein interactions called packaging signals (; Stockley et al., 2016; Martín Garrido et al., 2020). Packaging signals are also present in Bromoviridae (; ; ), Tombusviridae (Qu and Morris, 1997) and Tymoviridae (; ) although it is known that the RNA in Bromoviridae is much less ordered (). We do not consider the effect of these signals here.

Finally, we discuss possible applications of understanding CP-CP interactions for capsid engineering. MS2 (Leviviridae) capsids have been used to design vessels for delivering therapeutic molecules (Wu et al., 1995; ; Wu et al., 2005; ; Yan et al., 2014; Pumpens et al., 2016) and for antigen presentation in vaccines (Mastico et al., 1993; Wu et al., 1995; Peabody, 1997; Tumban et al., 2015; ; Pumpens et al., 2016). Similar studies have also been performed with T = 3 jelly-roll capsids (Ren et al., 2006; ; Lockney et al., 2011; Yildiz et al., 2012; Zeng et al., 2013; Sokullu et al., 2019; ). Understanding protein-protein interactions is an important step in improvising such systems, as was shown in a recent study on MS2 (), where insertions were made to tweak the dimerization interface, allowing the formation of larger capsid structures.

Conclusions

The capsids of the Leviviridae, Bromoviridae, Tymoviridae and Tombusviridae families have T = 3 icosahedral structures made up, almost completely, of a single type of coat-protein (CP). However, the structures of the Leviviridae CPs are distinct from those of the CPs of the other three virus families which fold into the jelly-roll fold. We examined the average hydrophobicity (H) and relative size (R) of all the CP-CP interfaces (dimerization, trimerization, pentamerization and hexamerization) and found that the dimerization interface of the Leviviridae CPs is distinctly different from all other interfaces with high H and R values indicating that these CPs are likely to be present as CP2 dimers. On an average, H values of the dimerization interfaces of the other three families (Bromoviridae, Tymoviridae and Tombusviridae) are also higher than that of the other interfaces potentially making this interface stronger. This could be the reason for capsid disassembly into CP2 in in vitro experiments. Although the CP monomers from these three virus families are structurally similar, differences in dimer interface structure make the dynamics of the dimerization interfaces distinct across the three families. Thus, even in viruses with the same overall capsid symmetry, a diversity in CP structure and interactions may lead to diverse assembly-disassembly mechanisms.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

DLP and SG designed the study. DLP performed data analysis and simulations. DLP and SG wrote the manuscript and approved the submitted version.

Funding

This work was funded by the Department of Atomic Energy, Government of India through the Tata Institute of Fundamental Research from Project Identification No. RTI 4006. SG is supported by a MATRICS grant from the Government of India-SERB (MTR/2019/001555). DLP is funded through an SPM fellowship from the Council of Scientific and Industrial Research, India. We also acknowledge support of the Simons Foundation (Grant No. 287975).

Acknowledgments

We thank Professor M. R. N. Murthy and Professor H. S. Savithri (IISc Bangalore, India) for detailed discussions and sharing their valuable insights on T = 3 icosahedral capsids.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.967877/full#supplementary-material

References

  • 1

    AdolphK. W.ButlerP. J. G. (1974). Studies on the assembly of a spherical plant virus. I. States of aggregation of the isolated protein. J. Mol. Biol.88, 327341. 10.1016/0022-2836(74)90485-9

  • 2

    AramayoR.MérigouxC.LarquetE.BronP.PérezJ.DumasC.et al (2005). Divalent ion-dependent swelling of tomato bushy stunt virus: A multi-approach study. Biochim. Biophys. Acta1724, 345354. 10.1016/j.bbagen.2005.05.020

  • 3

    BancroftJ. B.HillsG. J.MarkhamR. (1967). A study of the self-assembly process in a small spherical virus formation of organized structures from protein subunits in vitro. Virology31, 354379. 10.1016/0042-6822(67)90180-8

  • 4

    BarnhillH. N.ReutherR.FergusonP. L.DreherT.WangQ. (2007). Turnip yellow mosaic virus as a chemoaddressable bionanoparticle. Bioconjug. Chem.18, 852859. 10.1021/bc060391s

  • 5

    BerenC.CuiY.ChakravartyA.YangX.RaoA. L. N.KnoblerC. M.et al (2020). Genome organization and interaction with coat‐protein in a multipartite RNA virus. Proc. Natl. Acad. Sci. U. S. A.117, 1067310680. 10.1073/pnas.1915078117

  • 6

    BermanH. M.WestbrookJ.FengZ.GillilandG.BhatT. N.WeissigH.et al (2000). The protein data bank. Nucleic Acids Res.28, 235242. 10.1093/nar/28.1.235

  • 7

    BielaA. P.NaskalskaA.FatehiF.TwarockR.HeddleJ. G. (2022). Programmable polymorphism of a virus-like particle. Commun. Mat.3, 7. 10.1038/s43246-022-00229-3

  • 8

    BinkH. H. J.SchirawskiJ.HaenniA.-L.PleijC. W. A. (2003). The 5′-proximal hairpin of turnip yellow mosaic virus RNA: Its role in translation and encapsidation. J. Virol.77, 74527458. 10.1128/jvi.77.13.7452-7458.2003

  • 9

    BondK. M.LykteyN. A.TsvetkovaI. B.DragneaB.JarroldM. F. (2020). Disassembly intermediates of the brome mosaic virus identified by charge detection mass spectrometry. J. Phys. Chem. B124, 21242131. 10.1021/acs.jpcb.0c00008

  • 10

    BrownW. L.MasticoR. A.WuM.HealK. G.AdamsC. J.MurrayJ. B.et al (2002). RNA bacteriophage capsid-mediated drug delivery and epitope presentation. Intervirology45, 371380. 10.1159/000067930

  • 11

    CasparD. L.KlugA. (1962). Physical principles in the construction of regular viruses. Cold Spring Harb. Symp. Quant. Biol.27, 124. 10.1101/SQB.1962.027.001.005

  • 12

    ChakrabortyD.HoriN.ThirumalaiD. (2018). Sequence-dependent three interaction site model for single- and double-stranded DNA. J. Chem. Theory Comput.14, 37633779. 10.1021/acs.jctc.8b00091

  • 13

    ChakravartyA.RaoA. L. (2022). “Bromoviridae: A family of plant viruses with tripartite genomes,” in eLS. (John Wiley & Sons, Ltd (Ed.)). 10.1002/9780470015902.a0029278

  • 14

    ChandoniaJ. M.GuanL.LinS.YuC.FoxN. K.BrennerS. E. (2022). SCOPe: Improvements to the structural classification of proteins-extended database to facilitate variant interpretation and machine learning. Nucleic Acids Res.50, D553D559. 10.1093/nar/gkab1054

  • 15

    ChangJ. Y.GorzelnikK. v.ThongcholJ.ZhangJ. (2022). Structural assembly of qβ virion and its diverse forms of virus-like particles. Viruses14, 225. 10.3390/v14020225

  • 16

    ChaturvediS.RaoA. L. N. (2018). Molecular and biological factors regulating the genome packaging in single-strand positive-sense tripartite RNA plant viruses. Curr. Opin. Virol.33, 113119. 10.1016/j.coviro.2018.07.019

  • 17

    ChengS.BrooksC. L. (2013). Viral coat‐protein are segregated in structural fold space. PLoS Comput. Biol.9, e1002905. 10.1371/journal.pcbi.1002905

  • 18

    ChoiY. G.DreherT. W.RaoA. L. N. (2002). tRNA elements mediate the assembly of an icosahedral RNA virus. Proc. Natl. Acad. Sci. U. S. A.99, 655660. 10.1073/pnas.022618199

  • 19

    ChoiY. G.RaoA. L. N. (2003). Packaging of brome mosaic virus RNA3 is mediated through a bipartite signal. J. Virol.77, 97509757. 10.1128/jvi.77.18.9750-9757.2003

  • 20

    ClementiC.NymeyerH.OnuchicJ. N. (2000). Topological and energetic factors: What determines the structural details of the transition state ensemble and “en-route” intermediates for protein folding? An investigation for small globular proteins. J. Mol. Biol.298, 937953. 10.1006/jmbi.2000.3693

  • 21

    Durán-MezaA. L.Escamilla-RuizM. I.Segovia-GonzálezX. F.Villagrana-EscareñoM. v.Vega-AcostaJ. R.Ruiz-GarciaJ. (2020). Encapsidation of different plasmonic gold nanoparticles by the CCMV CP. Molecules25, E2628. 10.3390/molecules25112628

  • 22

    ElradO. M.HaganM. F. (2008). Mechanisms of size control and polymorphism in viral capsid assembly. Nano Lett.8, 38503857. 10.1021/nl802269a

  • 23

    ElSawyK. M.CavesL. S. D.TwarockR. (2010). The impact of viral RNA on the association rates of coat‐protein assembly: Bacteriophage MS2 as a case study. J. Mol. Biol.400, 935947. 10.1016/j.jmb.2010.05.037

  • 24

    FoxN. K.BrennerS. E.ChandoniaJ. M. (2014). SCOPe: Structural Classification of Proteins - extended, integrating SCOP and ASTRAL data and classification of new structures. Nucleic Acids Res.42, D304D309. 10.1093/nar/gkt1240

  • 25

    FrietzeK. M.PeabodyD. S.ChackerianB. (2016). Engineering virus-like particles as vaccine platforms. Curr. Opin. Virol.18, 4449. 10.1016/j.coviro.2016.03.001

  • 26

    GalawayF. A.StockleyP. G. (2013). MS2 viruslike particles: A robust, semisynthetic targeted drug delivery platform. Mol. Pharm.10, 5968. 10.1021/mp3003368

  • 27

    HumphreyW.DalkeA.SchultenK. (1996). Vmd: Visual molecular dynamics. J. Mol. Graph.14, 3338. 10.1016/0263-7855(96)00018-5

  • 28

    JohnsonJ. E.ChiuW. (2000). Structures of virus and virus-like particles. Curr. Opin. Struct. Biol.10, 229235. 10.1016/S0959-440X(00)00073-7

  • 29

    KingA. M. Q.LefkowitzE.AdamsM. J.CarstensE. B. (2012). “Virus taxonomy,” in Reference Module in life sciences (Amsterdam, Netherlands: Elsevier). 10.1016/b978-0-12-809633-8.21231-4

  • 30

    KoningR. I.Gomez-BlancoJ.AkopjanaI.VargasJ.KazaksA.TarsK.et al (2016). Asymmetric cryo-EM reconstruction of phage MS2 reveals genome structure in situ. Nat. Commun.7, 12524. 10.1038/ncomms12524

  • 31

    LarsonS. B.LucasR. W.GreenwoodA.McPhersonA. (2005). The RNA of turnip yellow mosaic virus exhibits icosahedral order. Virology334, 245254. 10.1016/j.virol.2005.01.036

  • 32

    LevyY.PapoianG. A.OnuchicJ. N.WolynesP. G. (2004a). Energy landscape analysis of protein dimers. Isr. J. Chem.44, 281297. 10.1560/XGCB-WUHN-JRJC-BV0Y

  • 33

    LevyY.WolynesP. G.OnuchicJ. N. (2004b). Protein topology determines binding mechanism. Proc. Natl. Acad. Sci. U. S. A.101, 511516. 10.1073/pnas.2534828100

  • 34

    LlauróA.CoppariE.ImperatoriF.BizzarriA. R.CastónJ. R.SantiL.et al (2015). Calcium ions modulate the mechanics of tomato bushy stunt virus. Biophys. J.109, 390397. 10.1016/j.bpj.2015.05.039

  • 35

    LockneyD. M.GuentherR. N.LooL.OvertonW.AntonelliR.ClarkJ.et al (2011). The red clover necrotic mosaic virus capsid as a multifunctional cell targeting plant viral nanoparticle. Bioconjug. Chem.22, 6773. 10.1021/bc100361z

  • 36

    LokeshG. L.GowriT. D. S.SatheshkumarP. S.MurthyM. R. N.SavithriH. S. (2002). A molecular switch in the coat‐protein controls the particle polymorphism in an icosahedral virus. Virology292, 211223. 10.1006/viro.2001.1242

  • 37

    LoutenJ. (2016). “Virus structure and classification,” in Essential human virology (Amsterdam, Netherlands: Elsevier Science). 10.1016/b978-0-12-800947-5.00002-8

  • 38

    Martín GarridoN.CroneM. A.RamlaulK.SimpsonP. A.FreemontP. S.AylettC. H. S. (2020). Bacteriophage MS2 displays unreported capsid variability assembling T = 4 and mixed capsids. Mol. Microbiol.113, 143152. 10.1111/mmi.14406

  • 39

    MasticoR. A.TalbotS. J.StockleyP. G. (1993). Multiple presentation of foreign peptides on the surface of an RNA-free spherical bacteriophage capsid. J. Gen. Virol.74, 541548. 10.1099/0022-1317-74-4-541

  • 40

    MengR.JiangM.CuiZ.ChangJ. Y.YangK.JakanaJ.et al (2019). Structural basis for the adsorption of a single-stranded RNA bacteriophage. Nat. Commun.10, 3130. 10.1038/s41467-019-11126-8

  • 41

    Montiel-GarciaD.Santoyo-RiveraN.HoP.Carrillo-TrippM.BrooksC. L.JohnsonJ. E.et al (2021). VIPERdb v3.0: A structure-based data analytics platform for viral capsids. Nucleic Acids Res.49, D809D816. 10.1093/nar/gkaa1096

  • 42

    NiC.-Z.SyedR.KodandapaniR.WickershamJ.PeabodyD. S.ElyK. R. (1995). Crystal structure of the MS2 coat‐protein dimer: Implications for RNA binding and virus assembly. Structure3, 255263. 10.1016/S0969-2126(01)00156-3

  • 43

    OkinakaY.MiseK.SuzukiE.OkunoT.FurusawaI. (2001). The C terminus of brome mosaic virus coat‐protein controls viral cell-to-cell and long-distance movement. J. Virol.75, 53855390. 10.1128/jvi.75.11.5385-5390.2001

  • 44

    PaciosL. F. (2001). Distinct molecular surfaces and hydrophobicity of amino acid residues in proteins. J. Chem. Inf. Comput. Sci.41, 14271435. 10.1021/ci010369n

  • 45

    PappachanA.ChinnathambiS.SatheshkumarP. S.SavithriH. S.MurthyM. R. N. (2009). A single point mutation disrupts the capsid assembly in Sesbania Mosaic Virus resulting in a stable isolated dimer. Virology392, 215221. 10.1016/j.virol.2009.06.047

  • 46

    PeabodyD. S. (1997). Subunit fusion confers tolerance to peptide insertions in a virus coat‐protein. Arch. Biochem. Biophys.347, 8592. 10.1006/abbi.1997.0312

  • 47

    PettersenE. F.GoddardT. D.HuangC. C.CouchG. S.GreenblattD. M.MengE. C.et al (2004). UCSF Chimera—A visualization system for exploratory research and analysis. J. Comput. Chem.25, 16051612. 10.1002/jcc.20084

  • 48

    PollesG.IndelicatoG.PotestioR.CermelliP.TwarockR.MichelettiC. (2013). Mechanical and assembly units of viral capsids identified via quasi-rigid domain decomposition. PLoS Comput. Biol.9, e1003331. 10.1371/journal.pcbi.1003331

  • 49

    PowellJ. D.BarbarE.DreherT. W. (2012). Turnip yellow mosaic virus forms infectious particles without the native beta-annulus structure and flexible coat‐protein N-terminus. Virology422, 165173. 10.1016/j.virol.2011.10.019

  • 50

    PrakashD. L.GosaviS. (2021). Understanding the folding mediated assembly of the bacteriophage MS2 coat‐protein dimers. J. Phys. Chem. B125, 87228732. 10.1021/acs.jpcb.1c03928

  • 51

    PumpensP.RenhofaR.DishlersA.KozlovskaT.OseV.PushkoP.et al (2016). The true story and advantages of RNA phage capsids as nanotools. Intervirology59, 74110. 10.1159/000449503

  • 52

    QuF.MorrisT. J. (1997). Encapsidation of turnip crinkle virus is defined by a specific packaging signal and RNA size. J. Virol.71, 14281435. 10.1128/jvi.71.2.1428-1435.1997

  • 53

    RenY.WongS. M.LimL. Y. (2006). In vitro-reassembled plant virus-like particles for loading of polyacids. J. Gen. Virol.87, 27492754. 10.1099/vir.0.81944-0

  • 54

    RobertsE.EargleJ.WrightD.Luthey-SchultenZ. (2006). MultiSeq: Unifying sequence and structure data for evolutionary analysis. BMC Bioinforma.7, 382. 10.1186/1471-2105-7-382

  • 55

    RoenhorstJ. W.VerduinB. J. M.GoldbachR. W. (1989). Virus-ribosome complexes from cell-free translation systems supplemented with cowpea chlorotic mottle virus particles. Virology168, 138146. 10.1016/0042-6822(89)90412-1

  • 56

    RolfssonÓ.MiddletonS.ManfieldI. W.WhiteS. J.FanB.VaughanR.et al (2016). Direct evidence for packaging signal-mediated assembly of bacteriophage MS2. J. Mol. Biol.428, 431448. 10.1016/j.jmb.2015.11.014

  • 57

    RossmannM. G.JohnsonJ. E. (1989). Icosahedral RNA virus structure. Annu. Rev. Biochem.58, 533573. 10.1146/ANNUREV.BI.58.070189.002533

  • 58

    RussellR. B.BartonG. J. (1992). Multiple protein sequence alignment from tertiary structure comparison: Assignment of global and residue confidence levels. Proteins14, 309323. 10.1002/prot.340140216

  • 59

    SastriM.ReddyD. S.KrishnaS. S.MurthyM. R. N.SavithriH. S. (1999). Identification of a discrete intermediate in the assembly/disassembly of physalis mottle tymovirus through mutational analysis. J. Mol. Biol.289, 905918. 10.1006/jmbi.1999.2786

  • 60

    SobolevV.SorokineA.PriluskyJ.AbolaE. E.EdelmanM. (1999). Automated analysis of interatomic contacts in proteins. Bioinformatics15, 327332. 10.1093/bioinformatics/15.4.327

  • 61

    SokulluE.AbyanehH. S.GauthierM. A. (2019). Plant/bacterial virus-based drug discovery, drug delivery, and therapeutics. Pharmaceutics11, E211. 10.3390/pharmaceutics11050211

  • 62

    SpeirJ. A.MunshiS.WangG.BakerT. S.JohnsonJ. E. (1995). Structures of the native and swollen forms of cowpea chlorotic mottle virus determined by X-ray crystallography and cryo-electron microscopy. Structure3, 6378. 10.1016/S0969-2126(01)00135-6

  • 63

    StockleyP. G.RolfssonO.ThompsonG. S.BasnakG.FranceseS.StonehouseN. J.et al (2007). A simple, RNA-mediated allosteric switch controls the pathway to formation of a T = 3 viral capsid. J. Mol. Biol.369, 541552. 10.1016/j.jmb.2007.03.020

  • 64

    StockleyP. G.WhiteS. J.DykemanE.ManfieldI.RolfssonO.PatelN.et al (2016). Bacteriophage MS2 genomic RNA encodes an assembly instruction manual for its capsid. Bacteriophage6, e1157666. 10.1080/21597081.2016.1157666

  • 65

    SzoverfiJ.FejerS. N. (2022). Dynamic stability of salt stable cowpea chlorotic mottle virus coat‐protein dimers and pentamers of dimers. Sci. Rep.12 (1), 14251. 10.1038/s41598-022-18019-9

  • 66

    TamaF.BrooksC. L. (2002). The mechanism and pathway of pH induced swelling in cowpea chlorotic mottle virus. J. Mol. Biol.318, 733747. 10.1016/S0022-2836(02)00135-3

  • 67

    TangJ.JohnsonJ. M.DrydenK. A.YoungM. J.ZlotnickA.JohnsonJ. E. (2006). The role of subunit hinges and molecular “switches” in the control of viral capsid polymorphism. J. Struct. Biol.154, 5967. 10.1016/j.jsb.2005.10.013

  • 68

    TumbanE.MuttilP.EscobarC. A. A.PeabodyJ.WafulaD.PeabodyD. S.et al (2015). Preclinical refinements of a broadly protective VLP-based HPV vaccine targeting the minor coat‐protein, L2. Vaccine33, 33463353. 10.1016/j.vaccine.2015.05.016

  • 69

    VirudachalamR.LowP. S.ArgosP.MarkleyJ. L. (1985). Turnip yellow mosaic virus and its capsid have thermal stabilities with opposite ph dependence: Studies by differential scanning calorimetry and 31P nuclear magnetic resonance spectroscopy. Virology146, 213220. 10.1016/0042-6822(85)90005-4

  • 70

    WuM.BrownW. L.StockleyP. G. (1995). Cell-specific delivery of bacteriophage-encapsidated ricin a chain. Bioconjug. Chem.6, 587595. 10.1021/bc00035a013

  • 71

    WuM.SherwinT.BrownW. L.StockleyP. G. (2005). Delivery of antisense oligonucleotides to leukemia cells by RNA bacteriophage capsids. Nanomedicine1, 6776. 10.1016/j.nano.2004.11.011

  • 72

    YanR.HallamA.StockleyP. G.BoyesJ. (2014). Oncogene dependency and the potential of targeted RNAi-based anti-cancer therapy. Biochem. J.461, 113. 10.1042/BJ20140173

  • 73

    YildizI.TsvetkovaI.WenA. M.ShuklaS.MasarapuM. H.DragneaB.et al (2012). Engineering of Brome mosaic virus for biomedical applications. RSC Adv.2, 36703677. 10.1039/c2ra01376b

  • 74

    ZengQ.WenH.WenQ.ChenX.WangY.XuanW.et al (2013). Cucumber mosaic virus as drug delivery vehicle for doxorubicin. Biomaterials34, 46324642. 10.1016/j.biomaterials.2013.03.017

  • 75

    ZhaoX.FoxJ. M.OlsonN. H.BakerT. S.YoungM. J. (1995). In vitro assembly of cowpea chlorotic mottle virus from coat‐protein expressed in Escherichia coli and in vitro-transcribed viral cDNA. Virology207, 486494. 10.1006/viro.1995.1108

  • 76

    ZlotnickA.AldrichR.JohnsonJ. M.CeresP.YoungM. J. (2000). Mechanism of capsid assembly for an icosahedral plant virus. Virology277, 450456. 10.1006/viro.2000.0619

Summary

Keywords

T=3 icosahedral viruses, coat-proteins, protein-protein interactions, dimerization interfaces, Leviviridae, Bromoviridae, Tombusviridae, Tymoviridae

Citation

Prakash DL and Gosavi S (2022) The diversity of protein-protein interaction interfaces within T=3 icosahedral viral capsids. Front. Mol. Biosci. 9:967877. doi: 10.3389/fmolb.2022.967877

Received

13 June 2022

Accepted

08 September 2022

Published

20 October 2022

Volume

9 - 2022

Edited by

Michael Blaber, Florida State University, United States

Reviewed by

Ajeet Sharma, Indian Institute of Technology Jammu, India

Sanjay Kumar Dey, University of Delhi, India

Updates

Copyright

*Correspondence: Shachi Gosavi,

This article was submitted to Structural Biology, a section of the journal Frontiers in Molecular Biosciences

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics