# On the role of VP3-PI3P interaction in birnavirus endosomal membrane targeting

Flavia Zanetti, Ignacio Fernandez, Eduard Baquero, Pablo Guardado-Calvo, Andres Ferrino-Iriarte, Sarah Dubois, Etienne Morel, Victoria Alfonso, Milton Osmar Aguilera, María Celayes, et al.

► To cite this version:

Flavia Zanetti, Ignacio Fernandez, Eduard Baquero, Pablo Guardado-Calvo, Andres Ferrino-Iriarte, et al.. On the role of VP3-PI3P interaction in birnavirus endosomal membrane targeting. eLife, 2025, 13, pp.RP97261. 10.7554/elife.97261. hal-05109598

HAL Id: hal-05109598

https://hal.science/hal-05109598v1

Submitted on 13 Jun 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# On the role of VP3-PI3P interaction in birnavirus endosomal membrane targeting

Flavia A Zanetti¹†, Ignacio Fernandez²†, Eduard Baquero², Pablo Guardado-Calvo², Andres Ferrino-Iriarte², Sarah Dubois³, Etienne Morel³, Victoria Alfonso⁴, Milton Osmar Aguilera⁵, María E Celayes⁵, Luis Mariano Polo⁵, Laila Suhaiman⁶, Vanesa V Galassi⁶,⁷, Maria V Chiarpotti⁶,⁷, Carolina Allende-Ballestero⁸, Javier M Rodriguez⁸, Jose R Castón⁸, Diego Lijavetzky⁹, Oscar Taboga⁴, María I Colombo⁵, Mario Del Pópolo⁶,⁷, Félix A Rey², Laura Ruth Delgui⁵,⁷*

¹Instituto de Ciencia y Tecnología "Dr. Cesar Milstein", Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires, Argentina; ²Institut Pasteur, Université Paris Cité, Structural Virology Unit, Paris, France; ³Université Paris Cité, INSERM UMR-S1151, CNRS UMR-S8253, Institut Necker Enfants Malades, Paris, France; ⁴Instituto de Agrobiotecnología y Biología Molecular (IABIMO), Instituto Nacional de Tecnología Agropecuaria (INTA), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires, Argentina; ⁵Instituto de Histología y Embriología de Mendoza, Universidad Nacional de Cuyo (UNCuyo), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Centro Universitario, Mendoza, Argentina; ⁶Instituto Interdisciplinario de Ciencias Básicas (ICB), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Mendoza, Argentina; ⁷Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo (UNCuyo), Mendoza, Argentina; ⁸Department of Structure of Macromolecules, Centro Nacional de Biotecnología (CNB-CSIC), Madrid, Spain; ⁹Instituto de Biología Agrícola de Mendoza, Universidad Nacional de Cuyo (UNCuyo), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Mendoza, Argentina

## eLife Assessment

Zanetti et al use **convincing** biophysical and cellular assays to investigate the interaction of the birnavirus VP3 protein with the early endosome lipid PI3P. The study provides **valuable** insights and will be of interest to virologists. In future studies, it would be interesting to demonstrate that VP3-PI3P is a specific interaction and not a general interaction with other PIPs.

**Abstract** Birnaviruses are a group of double-stranded RNA (dsRNA) viruses infecting birds, fish, and insects. Early endosomes (EE) constitute the platform for viral replication. Here, we study the mechanism of birnaviral targeting of EE membranes. Using the Infectious Bursal Disease Virus (IBDV) as a model, we validate that the viral protein 3 (VP3) binds to phosphatidylinositol-3-phosphate (PI3P) present in EE membranes. We identify the domain of VP3 involved in PI3P-binding, named P2 and localized in the core of VP3, and establish the critical role of the arginine at position 200 (R200), conserved among all known birnaviruses. Mutating R200 abolishes viral replication. Moreover, we propose a two-stage modular mechanism for VP3 association with EE. Firstly, the carboxy-terminal region of VP3 adsorbs on the membrane, and then the VP3 core reinforces the membrane

*For correspondence:
Idelgui@mendoza-conicet.gob.ar

†These authors contributed
equally to this work

**Competing interest:** See page 21

**Funding:** See page 22

**Preprint posted**

20 February 2024

**Sent for Review**

05 March 2024

**Reviewed preprint posted**

26 April 2024

**Reviewed preprint revised**

05 February 2025

**Version of Record published**

06 March 2025

**Reviewing Editor:** Mauricio
Comas-Garcia, Universidad
Autónoma de San Luis Potosí,
Mexico

© Copyright Zanetti, Fernandez et al. This article is distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use and redistribution provided that the original author and source are credited.
engagement by specifically binding PI3P through its P2 domain, additionally promoting PI3P accumulation.

Introduction

Birnaviruses belong to the *Birnaviridae* family composed of nonenveloped icosahedral dsRNA viruses that infect a wide range of vertebrate and invertebrate hosts. The prototypical and best-characterized family member, the IBDV, is the causative agent of Gumboro disease, a highly contagious immunosup-pressive disease that affects young chickens (*Gallus gallus*) causing significant damage to the poultry industry (*Delmas et al., 2019*).

Birnaviruses are unconventional dsRNA viruses, as the virions lack a 'core,' i.e., an inner protein
layer containing the genome, which instead forms a filamentous ribonucleoprotein complex (RNP)
with the multifunctional protein VP3, the RdRp VP1, and the dsRNA genome segments (Luque et al.,
2009b). IBDV contains a single ~70 nm diameter T=13 I icosahedral capsid surrounding a polyploid
bipartite dsRNA genome (segments A and B, 3.2, and 2.8 kbp, respectively) (Luque et al., 2009a).
The X-ray crystal structure of the capsid revealed that VP2 is the only component. VP2 has two main
domains, a 'shell' domain homologous to the coat proteins of +sRNA such as the T=3 nodaviruses,
and a 'tower' domain, like those present in the T=13 middle-layer protein of dsRNA viruses of the
order Reovirales (Coulibaly et al., 2005). VP3 is a 257-residues long polypeptide that plays multiple
roles during the viral life cycle. In the mature virion, VP3 is a component of the RNPs, where it is asso-
ciated with the dsRNA (Luque et al., 2009b). During viral assembly, VP3 acts as a scaffold protein
via electrostatic interactions between basic residues of pVP2 (the precursor form of VP2), and acidic
residues at the VP3 carboxy-terminal (Ct) region (Saugar et al., 2010). Additionally, VP3 interacts
with VP1 (Lombardo et al., 1999). VP3's multifunctional properties might be mediated by its capacity
to form oligomers (Casañas et al., 2008), or by its intrinsically disordered regions. The X-ray crystal
structure of the central region of the IBDV VP3 (residues 82–222) showed that it consists entirely of
α-helices connected by loops which can be divided into two structural domains linked by a flexible
hinge. The 36 Ct residues of the protein, a highly hydrophilic region rich in charged amino acids and
proline residues, were predicted to be disordered, so it was excluded from the crystallized construct.
The 81 amino-terminal (Nt) region of VP3 was cleaved during the crystallization process and was
absent in the crystal (Casañas et al., 2008). The X-ray crystal structure indicated the presence of two
positively charged domains, termed Patch 1 (P1) and Patch 2 (P2) (Valli et al., 2012). Both patches
are composed of four discontinuous, positively charged residues forming two exposed patches on
the surface of the protein. P1 was demonstrated to be involved in the dsRNA binding-activity of VP3
(Valli et al., 2012).

Previous studies on the subcellular localization of the IBDV replication complex (RC) have shown that IBDV replication components are localized to EE (Delgui et al., 2013). Further analysis revealed that VP3 associates with the cytosolic leaflet of EE membranes via P2, with a critical functional role in the virus life cycle (Gimenez et al., 2018). Additionally, the presence of PI3P in EE membranes was demonstrated to be a key host factor for VP3 association and the consequent establishment of IBDV RCs on EE membranes (Gimenez et al., 2021).

Here, we study the mechanism mediating the birnaviral targeting of EE membranes, focusing
on the VP3-PI3P interaction. Using biophysical approaches in cell-free in vitro systems, we demon-
strate that PI3P constitutes the cellular partner for VP3 association to EE membranes. By intro-
ducing charge-reversal point mutations on VP3 P2, we establish the contribution of each of the four
positively charged residues composing P2, with a key role of the arginine at position 200 (R200), a
residue conserved across all known birnaviruses. Using a reverse genetic system for IBDV, we show
that mutating R200 abolishes viral replication. Finally, by using molecular simulations, we discuss
a potential two-stage modular mechanism for VP3 association with EE. We show that, in the first
step, the Ct flexible tail of VP3 adsorbs to the membrane by electrostatic interactions due to its
positive charge. Then, in the second stage, the VP3 core reinforces the membrane engagement by
recruiting PI3P molecules scattered in the EE membrane to the P2 region, promoting a stable VP3-
PI3P interaction.
# Results

## Biophysical characterization of VP3 binding to PI3P

We set up a liposome co-flotation assay to assess the binding of VP3 to PI3P in membranes. We expressed and purified His-tagged VP3 full-length (His-VP3 FL) and confirmed its identity by Western blot using specific anti-VP3 and anti-His antibodies. Since we observed four different bands for the purified protein, we confirmed their VP3-identity by mass spectrometry (Figure 1—figure supplement 1). We prepared two populations of liposomes: 'liposomes PI3P(-)' and 'liposomes PI3P(+)'. We confirmed that these preparations were suitable for evaluating the interaction with PI3P by monitoring co-flotation of a recombinant His-2xFYVE domain (Fab1p, YOTB, Vac1p, and EEA1), followed by immunoblot detection of the His-tag. FYVE domains are highly homologous cysteine-rich domains of 70–80 amino acids present in around 30 human proteins that are known to participate in vesicular sorting or endocytosis through direct binding to PI3P from EE (Burd and Emr, 1998; Gaullier et al., 1998; Simonsen et al., 1998; Stenmark et al., 1996; Figure 1A). We observed specific binding of His-VP3 FL to liposomes PI3P(+), evidenced by a significant difference in the ability to bind to both liposome populations (Figure 1B). To test His-VP3 FL PI3P-binding specificity, we performed a co-flotation assay using liposomes where PI3P was replaced by 1,2-dioleoyl-sn-glycero-3-phosphate ('liposomes PA') or [1,2-dioleoyl-sn-glycero-3-phospho-(1'-myo-inositol)] ('liposomes PI') at the same molar ratio. PA contains a POH₄⁻ group also present in the polar head of PI3P but lacks the inositol ring, while PI has the inositol ring but is not phosphorylated. We observed that His-VP3 FL bound to liposomes PI3P(+), but not to liposomes PA or PI, reinforcing the notion that a phosphoinositide is required since neither a single negative charge nor an inositol ring are sufficient to promote VP3 binding to liposomes (Figure 1—figure supplement 2).

Two additional biophysical approaches were implemented. First, we prepared liposomes PI3P(-) and PI3P(+), and incubated them with His-2xFYVE or His-VP3 FL that had been pre-bound to Ni-NTA [nickel (II) nitrilotriacetic acid] gold particles. Subsequently, we inspected the preparations by cryo-electron microscopy (cryo-EM). Gold particles were not observed on liposomes PI3P(-) incubated with either protein (Figure 1C, left panel). We observed gold particles decorating liposomes PI3P(+) in the presence of both, His-2xFYVE and His-VP3 FL (Figure 1C, right panel). Second, we performed bio-layer interferometry (BLI) experiments using purified His-VP3 FL, and liposomes PI3P(-) and PI3P(+) (Shah and Duncan, 2014). In the BLI experiments, the protein is immobilized on the tip of Ni-NTA sensors while the liposomes remain in solution, allowing to monitor the molecular interactions in real-time (Kairys et al., 2019). Purified His-Streptavidin and His-2xFYVE were used as negative and positive controls, respectively. A specific liposome PI3P(+) dose-dependent binding of His-2xFYVE and His-VP3 FL was observed (Figure 1D).

Previous studies indicated that the P2 is responsible for PI3P binding on EE membranes (Gimenez et al., 2021; Gimenez et al., 2018). Yet the Ct region of VP3 is also highly enriched in positively charged residues. We generated a computational model of FL VP3 using AlphaFold2 (Jumper et al., 2021), which predicted a disordered Ct region (pLDTT of 42.6) (Figure 1E). Thus, we hypothesized that both conditions, the positive electrostatic surface charge and the disordered character of the VP3 Ct, might be contributing to the binding to the surface of EE. To test this hypothesis, we made a VP3 construct bearing a Ct truncation (His-VP3 ΔCt) and used it in co-flotation assays and BLI experiments. We did not detect an interaction of this mutant with PI3P(-) nor PI3P(+) liposomes (Figure 1F and G). Taken together, our results validate the role of PI3P in the interaction with VP3 and underscore the significance of the Ct domain in mediating this interaction.

## Role of VP3 P2 in the association of VP3 with the EE membrane

Inspection of the VP3 crystal structure showed that P2 residues (K<sub>157</sub>, R<sub>159</sub>, H<sub>198</sub>, and R<sub>200</sub>) form a positive-charged groove on the protein surface (Figure 2A and B). We prepared VP3 point mutants with reversed charge in a pcDNA vector backbone, giving rise to pcDNA VP3 K<sub>157</sub>D, R<sub>159</sub>D, H<sub>198</sub>D, and R<sub>200</sub>D. The calculated electrostatic potential of P2 for each of these mutants is shown in Figure 2C. Firstly, we used QM7 cells transiently overexpressing VP3 P2 wild-type (WT), P2 (all reversed), or the four single mutation of P2 to quantitatively assess their cellular distribution. We observed a punctuated distribution of VP3 P2 WT, K<sub>157</sub>D, R<sub>159</sub>D, and H<sub>198</sub>D. In contrast, the specific indirect immunofluorescence (IIF) signal of VP3 P2 all reversed and R<sub>200</sub>D showed a diffuse cytosolic distribution (Figure 2D, panels framed in red). Second, we assessed the distribution of the VP3 constructs in QM7
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Biophysical characterization of viral protein 3 (VP3) binding to phosphatidylinositol-3-phosphate (PI3P). (A) Immunoblots of the top (T) and bottom (B) fractions from a liposome PI3P(-) or PI3P(+) OptiprepTM co-floatation assay indicating that His-2xFYVE protein (~35 kDa) specifically binds to liposome PI3P(+). Results are representative of three independent experiments. The bar plot represents the intensity of T/B bands for each liposome preparation. Significant differences (**p<0.01) as determined by one-way ANOVA with Tukey's HSD test.(B) Immunoblots of the three top (T₁, T₂, and
*Figure 1 continued on next page*</figcaption></figure>
Figure 1 continued

$T_3$) and bottom (B) fractions from a liposome PI3P(-) or PI3P(+) OptiprepTM co-floatation assay indicating that His-VP3 FL protein (~32 kDa) specifically binds to liposome PI3P(+). Results are representative of three independent experiments. The bar plot represents the intensity of ($T_1+T_2+T_3$)/B bands for each liposome preparation. Significant differences (*p*<0.05) as determined by one-way ANOVA with Tukey's HSD test. (C) Liposomes were vitrified and only representative cryo-electron microscopy images are shown. Left panel: liposomes PI3P(-) incubated even with the Ni-NTA gold particles reagent alone or pre-incubated with His-2xFYVE of His-VP3 FL. Right panel: liposomes PI3P(+) control (without protein), or incubated with His-2xFYVE- or His-VP3 FL-Ni-NTA gold particles showing gold particles decorating the membrane of the liposomes when His-2xFYVE or His-VP3 FL were present. The bar represents 50 nm. (D) Binding of His-Streptavidin (negative control), His-2xFYVE (positive control), or His-VP3 FL to three different concentrations of liposomes PI3P(-) or PI3P(+). Association and dissociation sensorgrams measured by bio-layer interferometry (BLI), showing the specific interaction of His-2xFYVE and His-VP3 FL with liposomes PI3P(+) in a dose-dependent manner, as indicated. (E) Cartoon representation of AlphaFold2 prediction of VP3 FL. Red region corresponds to the 'core' region of the protein present in the experimental X-ray crystallographic model obtained by Casañas and coworkers PDB: 2R18 (Casañas et al., 2008). Regions not present within the PDB are colored in violet and blue representing the Nt and the Ct of VP3, respectively. pLDTT values lower than 50 are a strong predictor of disorder. (F) Binding of His-VP3 $\Delta$Ct to three different concentrations of liposomes PI3P(-) or PI3P(+). Sensorgrams measured by BLI, showing the absence of binding to either liposomes when the VP3 lacks the Ct region blue in (E). (G) Immunoblots of the three top ($T_1$, $T_2$, and $T_3$) and bottom (B) fractions from a liposome PI3P(-) or PI3P(+) OptiprepTM co-floatation assay of His-VP3 FL protein (~32 kDa) (positive control, left panel) or His-VP3 $\Delta$Ct protein (~28 kDa) (right panel) showing the lack of VP3 $\Delta$Ct binding to both liposomes. Results are representative of three independent experiments. The bar plot represents the intensity of ($T_1+T_2+T_3$)/B bands for each liposome preparation. Significant differences (*p*<0.05; ns p>0.05) as determined by one-way ANOVA with Tukey's HSD test.

The online version of this article includes the following source data and figure supplement(s) for figure 1:

**Source data 1.** Original membranes corresponding to **Figure 1A, B and G**.

**Source data 2.** Individual files corresponding to the original membranes from **Figure 1A, B and G**.

**Figure supplement 1.** His-viral protein 3 (VP3) Full length (FL) purification and characterization.

**Figure supplement 1—source data 1.** Original Coomasie R. Blue-stained polyacrylamide gel and western blot membranes corresponding to **Figure 1—figure supplement 1B and C**.

**Figure supplement 1—source data 2.** Individual files corresponding to the original Coomasie R. Blue-stained polyacrylamide gel and western blot membranes from **Figure 1—figure supplement 1B and C**.

**Figure supplement 2.** His-viral protein 3 (VP3) full length (FL) phosphatidylinositol-3-phosphate (PI3P)-binding specificity.

**Figure supplement 2—source data 1.** Original western blot membranes corresponding to **Figure 1—figure supplement 2**.

**Figure supplement 2—source data 2.** Individual files corresponding to the original western blot membranes from **Figure 1—figure supplement 2**.

cells transiently overexpressing EGFP-Rab5. Rab5 regulates the homotypic tethering and fusion of EEs, and it is considered a genuine marker of EEs (Gorvel *et al.*, 1991). As observed in **Figure 2E**, and the quantitative colocalization analysis (**Figure 2—figure supplement 1**), VP3 P2 WT, $K_{157}D$, $R_{159}D$, and $H_{198}D$ showed a marked colocalization with EGFP-Rab5. In contrast, the VP3 P2 all reversed and $R_{200}D$-derived IIF signal indicated a cytosolic distribution (**Figure 2E**, panels framed in red). It is important to mention that the differences in co-localization shown in **Figure 2—figure supplement 1** are small and not significant, but still show a trend that correlates with observations depicted in **Figure 2D**. Taken together, these results suggest that $R_{200}$ in the P2 patch bears a critical role in the association of VP3 with EE membranes.

## VP3 P2 mediates VP3-PI3P association to EE membranes

To further assess the contribution of P2 residues to PI3P binding, and to strengthen the observations made on cells transiently overexpressing EGFP-Rab5, we set up two additional cell biology assays. One consists of co-transferring cells with EGFP-2xFYVE and VP3 constructs, and detecting them by confocal laser scanning microscopy (CLSM). The second one consists of assessing the ability of purified GST-2xFYVE, which recognizes endogenous PI3P, to co-localize with VP3. We employed QM7 cells transiently co-overexpressing EGFP-2xFYVE and VP3 P2 WT, P2 all reversed, or the four-point mutants to evaluate their colocalization. As shown in **Figure 3A**, we observed a punctuated distribution with co-localization with EGFP-2xFYVE for VP3 P2 WT while for VP3 P2 all reversed the IIF signal indicated a cytosolic distribution of the protein with the loss of colocalization with EGFP-2xFYVE (**Figure 3A**, red framed panels on the upper two rows). Regarding the point mutants, we observed a loss of punctuated distribution and EGFP-2xFYVE co-localization for VP3 FL $R_{200}D$ (**Figure 3A**, red framed panels on the bottom row). Then, we used the GST-2xFYVE purified domain to recognize endogenous PI3P in fixed cells (Nascimbeni *et al.*, 2017). GST-2xFYVE was used as the 'primary antibody,' which was labeled with a fluorescent anti-GST antibody. As an additional control in these
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Viral protein 3 (VP3) P2 involvement in the association of VP3 with the early endosomes (EE) membrane. (A) The cartoon representation of the VP3 dimer PDB 2R18 (Casañas et al., 2008) shows each protomer in different shades of gray. The blue balls depict the residues defining the P2 region. (B) The close-up of P2 region showing residues K<sub>157</sub>, R<sub>159</sub>, H<sub>198</sub>, and R<sub>200</sub>. (C) Electrostatic potential mapped on the surface of the structure of a VP3 dimer in the same orientations as in (A) structures. The close-up shows the impact of P2 residue mutations on the electrostatic potential of the binding site. P2 wild-type (WT) corresponds to the 'UniProt code'. The color-coded electrostatic surface potential of VP3 was drawn using PyMol (blue positive, red negative). (D) QM7 cells transfected with pcDNA VP3 FL (P2 WT), P2 (all reversed), or the four-point mutants (K<sub>157</sub>D, R<sub>159</sub>D, H<sub>198</sub>D, and R<sub>200</sub>D) and immunostained with anti-VP3 showing the distribution of each protein (upper panel). Images were captured using a Confocal Laser Scanning</figcaption></figure>

*Figure 2 continued on next page*
*Figure 2 continued*

Microscopy and then the percentage of cells with punctated fluorescent signal were determined for each protein (lower panel). The red signal shows the VP3 distribution and the blue one shows the nuclei, which were Hoestch-stained. The data were normalized to the P2 WT protein. The box plot represents the percentage of cells with punctuated distribution of VP3. Significant differences (***p<0.001; ns p>0.05.) as determined by one-way ANOVA with Tukey's HSD test. (E) QM7 cells co-transfected with pEGFP-Rab5 and pcDNA VP3 FL (P2 WT), P2 (all reversed), or the four-point mutants ($K_{157}D$, $R_{159}D$, $H_{198}D$, and $R_{200}D$) and immunostained with anti-VP3 showing the distribution of each protein. Representative images, captured using a Confocal Laser Scanning Microscopy are shown where green signal represents Rab5 distribution and the red signal that of VP3. Nuclei were Hoestch-stained and are blue. White bar-scales represent 20 mm. VP3 P2 (all reversed) and $R_{200}D$ depict a cytosolic distribution of the proteins. Quantification of the co-localization of the different VP3 proteins and EGFP-Rab5 is shown in **Figure 2—figure supplement 1**.

The online version of this article includes the following figure supplement(s) for figure 2:

**Figure supplement 1.** VP3-EGFP-Rab5 co-localization quantification.

experiments, we immunostained the EE with antibodies anti-EEA1. As shown in **Figure 3B**, VP3 P2 WT demonstrated a conspicuous PI3P-bearing EEA1 positive endosomes location. In addition, we observed a similar distribution for VP3 $K_{157}D$, $R_{159}D$, and $H_{198}D$ and a drastic change in VP3 FL $R_{200}D$ mutant, which appeared utterly cytosolic, losing its capacity to bind to PI3P-bearing EE membranes (**Figure 3B**, red framed panel on the lowest row). These results further support the role of $R_{200}$ in this interaction.

## VP3 $R_{200}$ is crucial for binding to PI3P

To validate this, we devised a liposome co-flotation assay using a recombinant His-VP3 FL $R_{200}D$ mutant. When assayed in parallel with His-2xFYVE, we observed that His-VP3 FL $R_{200}D$ lost binding to PI3P(+) liposomes (**Figure 4A**). To confirm that the lack of binding was not due to misfolding of the mutant, we compared the circular dichroism spectra of His-VP3 FL and His-VP3 FL $R_{200}D$ proteins, without detecting significant differences (**Figure 4B**). To assess its relevance in IBDV infectivity, we resorted to an IBDV reverse genetics system, based on a modification of the system described by Qi and coworkers (*Qi et al.*, 2007). The full-length sequences of the IBDV RNA segments A and B, flanked by a hammerhead ribozyme at the 5'-end and the hepatitis delta ribozyme at the 3'-end, were expressed under the control of an RNA polymerase II promoter within the plasmids pCAGEN.Hmz.SegA.Hdz (SegA) and pCAGEN.Hmz.SegB.Hdz (SegB). For this specific experiment, we generated a third plasmid, pCAGEN.Hmz.SegA.R$_{200}$.D.Hdz (SegA.R$_{200}$.D), harboring a mutant version of segment A cDNA containing the $R_{200}D$ substitution. Then, QM7 cells were transfected with the plasmids SegA, SegB, or Seg.R$_{200}$.D alone (as controls) or with a mixture of plasmids SegA +SegB (wild-type situation) or SegA.R$_{200}$.D+SegB (mutant situation). At 8 hr post-transfection (p.t.), when the new viruses have been able to assemble out of the two segments of RNA, the cells were recovered and re-plated onto fresh non-transfected cells to reveal the presence of infective viruses. At 72 hr post-plating, the generation of foci forming units (FFUs) was revealed by Coomassie staining. As expected, single-transfections of SegA, SegB, or Seg.R$_{200}$.D did not produce FFUs. As shown in **Figure 4C**, the transfection of SegA +SegB produced detectable FFUs (the three circles in the upper panel) while no FFUs (the three circles in the lower panel) were detected after the transfection of SegA.R$_{200}$.D+SegB (**Figure 4C**). Given the relevance of VP3 $R_{200}$, we assessed its conservation among members of the Birnaviridae family. For this, a multiple sequence alignment of VP3 reference sequences was performed (**Figure 4—figure supplement 1A**). As shown, IBDV VP3 has 32-39% identity to the other selected homologues (**Figure 4—figure supplement 1B**). A closer look at the VP3 P2 revealed the complete conservation of $R_{200}$ (**Figure 4D**). Taken together, these findings show that VP3 $R_{200}$ plays an important role in mediating the VP3-PI3P interaction, highlighting the active participation of EE membranes during the infectious cycle of birnaviruses.

## Electrostatic forces mediate VP3 membrane binding

Computer simulations based on coarse-grain models were performed to assess the binding free-energy and adsorption mechanism of three protein variants: VP3 Ct, VP3 $\Delta$Nt, and VP3 FL. The VP3 Ct region, with a positive electrical charge of +5, was modeled as an intrinsically disordered region according to the AlphaFold2 prediction (**Figure 1E**). VP3 $\Delta$Nt encompasses residues 82-257 and is electroneutral. This variant was chosen primarily because it retains the structural data from the crystallized protein. Similarly to His-VP3 FL, His-VP3 $\Delta$Nt retains the ability to bind to liposomes PI3P(+)
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Viral protein 3 (VP3) P2 involvement in the association between VP3 and early endosomes (EE) phosphatidylinositol-3-phosphate (PI3P). (A) QM7 cells co-transfected with the PI3P biosensor pEGFP-2xFYVE and pcDNA VP3 FL (P2 WT), P2 (all reversed), or the four-point mutants (K<sub>157</sub>D, R<sub>159</sub>D, H<sub>198</sub>D, and R<sub>200</sub>D) and immunostained with anti-VP3 showing the distribution of each protein. Representative images captured using a Confocal Laser Scanning Microscopy are shown where blue signal represents FYVE distribution and the magenta signal that of VP3. Nuclei were Hoestch-stained and are orange (upper panels). White bar-scales represent 20 mm. VP3 P2 (all reversed) and R<sub>200</sub>D depict a cytosolic distribution of the proteins with a significant lower co-localization coefficient. The dot plot in the lower panel depicts the co-localization coefficient for each protein determined as explained in the Materials and methods section. Significant differences (ns p>0.05) as determined by one-way ANOVA with Tukey's HSD test. (B) QM7 cells were transfected with the pcDNA VP3 FL (P2 WT), P2 (all reversed), or the four-point mutants (K<sub>157</sub>D, R<sub>159</sub>D, H<sub>198</sub>D, and R<sub>200</sub>D). The cells were fixed and then the GST-2xFYVE purified peptide and anti-VP3 antibodies were used to recognize endogenous PI3P and VP3, respectively. Additionally, anti-EEA1 antibodies were used to stain the endosomes. GST-2xFYVE was labelled with a fluorescent anti-GST antibody (blue signal), anti-VP3, and anti-EEA1 antibodies with fluorescent secondary antibodies (magenta and green signals, respectively), and Hoestch-stained nuclei in orange. White bar-scales represent 10 mm. VP3 P2 (all reversed) and R<sub>200</sub>D depict a cytosolic distribution of the proteins with a significant lower co-localization coefficient. The dot plot in the lower panel depicts the co-localization coefficient for each protein determined as explained in the Materials and methods section. Significant differences (***p<0.001; ns p>0.05) as determined by one-way ANOVA with Tukey's HSD test.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Biophysical characterization of viral protein 3 (VP3) full length (FL) R~200~D binding to phosphatidylinositol-3-phosphate (PI3P). (A, left panel). Immunoblots of the top (T) and bottom (B) fractions from a liposome PI3P(-) or PI3P(+) OptiprepTM co-floatation assay indicating that His-2xFYVE protein (~35 kDa) specifically binds to liposome PI3P(+). The bar plot represents the intensity of T/B bands for each liposome preparation. Significant differences (ns p>0.05) as determined by one-way ANOVA with Tukey's HSD test. (A, right panel). Immunoblots of the three top (T~1~, T~2~, and T~3~) and bottom (B) fractions from a liposome PI3P(-) or PI3P(+) OptiprepTM co-floatation assay indicating that His-VP3 FL R~200~D protein (~35 kDa) does not bind to liposome PI3P(-) nor PI3P(+). The bar plot represents the intensity of (T~1~+T~2~+T~3~/B bands for each liposome preparation. Significant differences (ns p>0.05) as determined by one-way ANOVA with Tukey's HSD test. (B) Far-UV CD spectra of His-VP3 FL (red line) or His-VP3 FL R~200~D (green line). Spectral acquisitions at 50 nm/min with 0.1 nm steps at 1 s integration time, with a bandwidth of 1 nm were performed four times for the samples as well as for the buffer. The measurements were carried out with constant nitrogen gas flux of 10 ml/min. Acquisitions were averaged and buffer baseline was subtracted with Spectra Manager (JASCO). No smoothing was applied. CDtoolX was used to zero between 255–260 nm and to calibrate the signal amplitude from the fresh CSA signal (Miles and Wallace, 2018). Data are presented as delta epsilon (Δε) per residue (L.mol⁻¹.cm⁻¹.residue⁻¹) calculated using the molar concentration of protein and number of residues. (C) QM7 cells were grown in M24 multi-well plate for 12 hr to approximately 90–95% confluency and then 800 ng of plasmids were transfected [(SegA +SegB) or (SegA.R~200~D+SegB)] in triplicate. At 8 hr post-transfection (p.t.) the supernatants were discarded, and the monolayers were recovered for further plating on M6 multi-well plates containing non transfected QM7 cells. Avicel RC-591 (FMC Biopolymer) was added to the M6 multi-well plates. 72 hr p.i., the monolayers were fixated and stained with Coomassie R250 for revealing the foci forming units. (D) Partial view of amino acid alignment of the VP3 protein for nine reference members of *Birnaviridae* family. Multiple sequence alignment was performed with Clustal OMEGA (v1.2.4) (Sievers et al., 2011) implemented at EMBL's European Bioinformatics Institute (Thakur et al., 2024) (The complete alignment is shown in **Figure 4—figure supplement 1A**). Alignment visualization was done with the R ggmsa</figcaption></figure>

Figure 4 continued on next page
Figure 4 continued

package (Zhou et al., 2022) in with assistance from the RStudio software (RStudio Team, 2020). Amino acids are colored according to their side-chain chemistry. Protein sequence logos annotation is displayed on top of the amino acid alignment. For facilitating the view IBDV VP3 142–210 portion is shown. The black arrow indicates the K<sub>157</sub>, R<sub>159</sub>, H<sub>198</sub>, and R<sub>200</sub>.

The online version of this article includes the following source data and figure supplement(s) for figure 4:

**Source data 1.** Original membranes corresponding to **Figure 4A**.

**Source data 2.** Individual files corresponding to the original membranes from **Figure 4A**.

**Figure supplement 1.** Bioinformatic analysis.

(Figure 5—figure supplement 1). Finally, VP3 FL is the full-length monomeric VP3, using the Alpha-Fold2 predicted structure with an overall electrostatic charge of -3. We employed the Molecular Theory (MT) approach to obtain adsorption free-energy profiles, PMF(z), at different salt concentrations (Chiarpotti et al., 2021; Ramírez et al., 2019). The membrane was modeled as a rigid dielectric slab with 5% of its surface occupied by ternary ionizable groups. **Figure 5** (panels A-C) shows the PMF(z), where z represents the distance between the center of mass of the protein and the center of the membrane, at 50 and 150 mM NaCl. As observed in **Figure 5A**, the positively charged VP3 Ct shows the highest affinity for the membrane due to the deprotonation of the acidic lipids at pH 8, resulting in a net negative charge on the membrane surface (Figure 5—figure supplement 2). However, the electroneutral construct VP3 ΔNt has a lower binding energy than VP3 Ct, and also VP3 FL shows some degree of membrane binding, despite both the protein and the membrane being negatively charged, as evidenced by the shallow minimum in the PMF.

The binding of the three protein constructs was strongly affected by electrostatics, as shown in **Figure 5D**, where the adsorption energy is plotted as a function of the concentration of salt. The effect of salt concentration is most significant for VP3 Ct, but it is also noticeable for VP3 ΔNt and VP3 FL. Also, the interfacial concentration of each construct was higher at lower ionic strengths, as depicted in **Figure 5—figure supplement 3**. The weak binding energy of VP3 FL to the anionic membrane under physiological conditions (~150 mM NaCl, pH 7) was consistent with the experimental observation that VP3 FL did not bind to liposomes PA and PI, which are anionic at the pH of the experiments. However, those same experiments showed that VP3 FL did bind to PI3P(+) liposomes (Figure 1B-C), despite the MT calculations predicting weak binding. Taken together, our results suggest that specific VP3-PI3P interactions enhance the binding of VP3 to PI3P-containing EE membranes, while exhibiting weak and non-specific electrostatic attractions to an otherwise generic anionic membrane. **Figure 5—figure supplement 4** illustrates the distribution of cationic residues of VP3 near the membrane, indicating that the Ct fragment drives the non-specific binding. To identify the electrical charge of different protein regions, **Figure 5—figure supplement 5** shows a map of the electrostatic potential around VP3 FL.

Molecular Dynamics (MD) simulations based on the MARTINI model were employed to further investigate the mechanism by which VP3 Ct and VP3 ΔNt approach the membrane surface. We chose VP3 ΔNt over VP3 FL because VP3 ΔNt has a closer resemblance to the well-established experimental crystal structure. Additionally, the adsorption free energy profiles predicted by MT for both VP3 ΔNt and VP3 FL are notably similar (compare panels B and C of Figure 5). Unlike the MT model, MARTINI considers both electrostatic and van der Waals forces (de Jong et al., 2013). To mimic the liposome PI3P(+) composition used in the co-flotation assays, the lipid bilayer was modeled with a molar ratio of 64:31:5 for the lipids DOPE, DOPC, and PI3P, respectively. DOPE and DOPC are electroneutral, PI3P is anionic. The MD simulations revealed that both VP3 Ct and VP3 ΔNt settled on the membrane surface after 120 ns, as shown in **Figure 6—figure supplement 1**. Figure 6A-D shows a temporal sequence of configurations during the adsorption of VP3 ΔNt. The protein approached the membrane via the Ct fragment and remained adsorbed during the rest of the 500 ns trajectory. VP3 Ct exhibited the same behavior, as demonstrated in **Figure 6—figure supplements 1 and 2**. Notably, VP3 Ct and VP3 ΔNt did not detach from the surface during these MD simulations, suggesting that the adsorption energies calculated by MT were underestimated and represent a lower bound, as van der Waals interactions are absent in the MT model. Also, **Figure 6E** highlights the location of the four positive residues that make up the P2 region, showing their proximity to the membrane surface.
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Adsorption of viral protein 3 (VP3) constructs to phosphatidylinositol-3-phosphate (PI3P)(+) model membranes. (A–C). Adsorption free-energy profiles, PMF(z), for VP3 Ct, VP3 ΔNt, and VP3 FL at 50 and 150 mM of NaCl. (D) Adsorption free energy (ΔG<sub>ad</sub>), computed from the minimum of PMF(z), versus concentration of NaCl. In all cases, the solution pH is 8, and the concentration of protein is 1 µM. The gray areas represent the volume excluded by half the membrane (cis hemilayer, z>0). The membrane surface contains 5% titratable groups, representing PI3P. Each group has three acidic moieties, one with a pKa of 2.5 and the others with 6.5. At 150 mM NaCl and pH 8, more than 90% of the acidic groups are deprotonated (Figure 5—figure supplement 2).</figcaption></figure>

The online version of this article includes the following source data and figure supplement(s) for figure 5:

**Figure supplement 1.** Biophysical characterization of His-viral protein 3 (VP3) DNt binding to phosphatidylinositol-3-phosphate (PI3P).

**Figure supplement 1—source data 1.** Original western blot membranes corresponding to Figure 5—figure supplement 1A.

**Figure supplement 1—source data 2.** Individual files corresponding to the original western blot membranes from Figure 5—figure supplement 1A.

**Figure supplement 2.** Titration curve of the membrane.

**Figure supplement 3.** Interfacial concentration of viral protein 3 (VP3) constructs.

**Figure supplement 4.** Distribution of charged and aromatic residues of viral protein 3 (VP3).

**Figure supplement 5.** Map of electrostatic potential around viral protein 3 (VP3) full length (FL).

## Recruitment of PI3P modulates the membrane binding of VP3

We investigated the impact of VP3 Ct and VP3 ΔNt on the composition of the membrane near the binding point, and on the local curvature of PI3P(+) liposomes. We resorted to our MD simulations and evaluated the in-plane radial distribution function, g(r), between the center of mass of the proteins and the center of mass of PI3P molecules in the (cis) hemilayer, which is in direct contact with the protein. Figure 6F shows that VP3 Ct recruits PI3P within a~2 nm radius, where its local concentration in the cis hemilayer is up to 6 times greater than the average PI3P concentration (Figure 6—figure supplement 3). Also, despite having no net electrical charge, VP3 ΔNt induces a similar concentration enhancement, with PI3P accumulation occurring around the Ct fragment, as evidenced in Figure 6—figure supplement 4. In both cases, the recruitment of PI3P in the cis hemilayer is accompanied by
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Viral protein 3 (VP3) approaching a lipid bilayer, and distortion of the membrane. (A–D) Temporal sequence of configurations illustrating how VP3 ΔNt approaches the negatively charged membrane during a 500 ns molecular dynamics (MD) simulation. The membrane contains 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and phosphatidylinositol-3-phosphate (PI3P) in 64:31:5 molar ratio. The beads making the carboxy-terminal (Ct) fragment are colored blue. (E) Magnification of the protein configuration bound to the membrane, with P2 residues represented in cyan balls. (F) Radial distribution function, g(r), between the center of mass of VP3 ΔNt and the center of PI3P molecules in the cis hemilayer. g(r) greater than 1 implies local enhancement of PI3P concentration. The inset shows an upper view of VP3 ΔNt bound to the membrane, with the area within r=4 nm colored in orange. (G) Cryo-electron microscopy images of cryo-fixated liposomes PI3P(+) control (without protein, left panel), or incubated with His-VP3 FL showing the small pinches or localized thinnings in the bilayer of the liposomes when His-VP3 FL was present (middle panel). The bar represents 50 nm. The right panel represents an enlarged image of the red square on the middle panel. White arrows point to the small pinches or localized thinnings in the bilayer when His-VP3 FL was present.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 6:

Figure 6 continued on next page
Figure 6 continued

**Figure supplement 1.** Distance between membrane and protein constructs computed by molecular dynamics (MD).

**Figure supplement 2.** Initial and final configurations of viral protein 3 (VP3) carboxy-terminal (Ct) on a lipid membrane.

**Figure supplement 3.** Lipids accumulation/depletion around viral protein 3 (VP3) carboxy-terminal (Ct).

**Figure supplement 4.** Lipids accumulation/depletion around the carboxy-terminal (Ct) fragment of viral protein 3 (VP3) DNt.

**Figure supplement 5.** Lipids accumulation/depletion around viral protein 3 (VP3) ΔNt and its mutants.

**Figure supplement 6.** Adsorption free energy of viral protein 3 (VP3) ΔNt and VP3 ΔNt mutants on PI3P(+) model membranes.

a local depletion of PI3P in the trans hemilayer, due to interhemilayer electrostatic repulsions. The effects on PI3P concentration, induced by VP3 ΔNt, result from a synergistic action between the positively charged Ct and the P2 region of the protein. Furthermore, when mutations are introduced to the P2 region - starting with R<sub>200</sub>D and eventually including the entire P2 domain - there is a gradual reduction in its contribution to the recruitment of PI3P (Figure 6—figure supplement 5). Molecular Theory calculations on these mutants also reveal that the binding free energy of both VP3 ΔNt P2 Mut and VP3 ΔNt R<sub>200</sub>D is consistently lower than that of VP3 ΔNt across all salt concentrations explored (Figure 6—figure supplement 6). This indicates a weakening of the interaction with the negatively charged membrane, concomitant with the loss of positively charged amino acids.

Finally, to examine the potential impact of VP3 on the lipid bilayer, we incubated liposomes PI3P(-) and PI3P(+) with purified His-VP3 FL as described previously. The samples were cryo-fixed and cryo-EM was used to observe the resulting structures. As shown in Figure 6G, His-VP3 FL distorts the membrane of liposomes containing PI3P, resulting in small pinches in the bilayer, detectable in 36% of the liposomes PI3P(+), while it is completely absent in PI3P(-) liposomes. Taken together, our findings suggest that the binding of VP3 to EE membranes has the potential to alter the local lipid composition by recruiting PI3P, which could then lead to the localized membrane pinches observed in Figure 6G.

Discussion

Several viruses hijack the PI3P metabolism machinery to complete their infectious cycle. However, only a few cases of viral proteins that bind PI3P have been described. Equine infectious anemia virus (EIAV) is a member of the lentivirus subgroup of retroviruses and replicates in macrophages. Retroviral assembly and release are directed by the structural precursor polyprotein Gag (Wills and Craven, 1991). EIAV Gag localizes to both the cell interior and to the plasma membrane due to its high-affinity interaction with PI3P (Fernandes et al., 2011; Puffer et al., 1998; Tanzi et al., 2003). Moreover, mutation of K<sub>49</sub>, a residue in the PIs-binding pocket of Gag inhibited the release of viral-like particles from the plasma membrane (Fernandes et al., 2011). Another example is the poxviruses, a family of cytoplasmic DNA viruses that also rely on intracellular membranes to develop their envelope, whose morphogenesis requires enzymes from the cellular phosphoinositide metabolic pathway (McNulty et al., 2010). For vaccinia virus (VACV), the prototypic poxvirus, Kolli and coworkers showed that the VACV H7 protein binds to PI3P and PI4P (Kolli et al., 2015). Similarly to our observations in IBDV VP3, the authors uncovered two regions of H7 essential for viral replication. One region is a positive surface patch centered on K<sub>108</sub>, and the other region is the flexible Ct tail (Kolli et al., 2015).

Our biophysical and molecular simulation results suggest a specific but weak interaction of VP3 with PI3P-bearing membranes. In co-floatation assays, we observed that His-VP3 FL reaches the top of the gradient, where liposomes PI3P(+) are located, with a high proportion of protein still found at the bottom of the gradient. This phenomenon was not observed with the His-2xFYVE protein. Additionally, the slopes of the BLI curves show that His-VP3 FL dissociates faster than His-2xFYVE from liposomes PI3P(+). Computer simulations have also revealed a weak attractive interaction between the negatively charged VP3 and membranes that contain PI3P. VP3 bears multiple essential roles during the viral life cycle. In this context, our hypothesis is that an early, specific, but weak interaction with PI3P is required for the formation of RCs associated with EE and represents an advantage for the virus that counts with VP3 to go forward into the replication cycle.

We identified the domain of VP3 involved in PI3P-binding, and established the critical role played by the arginine residue at position 200 of VP3 (R<sub>200</sub>). This arginine is conserved among all known birna-viruses. Mutating R<sub>200</sub> abolishes viral replication. As previously mentioned, VP3 is a multifunctional viral
protein. Therefore, it was initially possible that mutating R<sub>200</sub> might affect other functions beyond PI3P binding and viral replication. However, R<sub>200</sub> is not located in any of the protein regions associated with its other known functions. For instance, the VP3 dimerization domain, located in the second helical domain, involves 81 interprotomeric contacts across three helices, but R<sub>200</sub> does not participate in these contacts (Casañas et al., 2008). Also, VP3 has an oligomerization domain mapped within the 42 C-terminal residues of the polypeptide, i.e., the segment of the protein composed by the residues at positions 216–257 (Maraver et al., 2003); VP3's ability to bind RNA is facilitated by a region of positively-charged amino acids, identified as P1, which includes K<sub>99</sub>, R<sub>102</sub>, K<sub>105</sub>, and K<sub>106</sub> (Valli et al., 2012). Furthermore, our findings indicate that the R<sub>200</sub>D mutant retains a folding pattern similar to the wild-type protein. All these lead us to conclude that the loss of replication capacity of R<sub>200</sub>D viruses results from impaired, or even loss of, VP3-PI3P interaction.

It has recently been reported that the RCs of IBDV exhibit liquid-liquid phase separation (LLPS) (Reddy et al., 2022). In LLPS, one or more biomolecules form a network of mild homo- and heterotypic interactions, resulting in the formation of a distinct liquid phase within a liquid medium (Banani et al., 2017). In light of our results, it is conceivable to hypothesize that viral interaction with the EE membranes acts as a first step in the RCs biogenesis, leading to the growth of biomolecular condensates with LLPS characteristics around EE membranes, as recently proposed (Brodrick and Broadbent, 2023). In support of this model, both the RNA-dependent RNA polymerase VP d the viral dsRNA associate to endosomes and to the LLPS condensates (Delgui et al., 2013; Reddy et al., 2022). Further experiments are needed in order to link the association of VP3 to EE membranes with biomolecular condensate formation.

## Evolutionary links between +sRNA and dsRNA viruses

Birnaviruses integrate the double-strand RNA Baltimore group (Baltimore, 1971), and the newly created Riboviria taxon (a realm), that includes all RNA viruses encoding RdRps (Walker et al., 2021). However, the organization of the birnaviral RdRp (Gorbalenya et al., 2002) and the structure of its capsid (Coulibaly et al., 2005) had already indicated strong links between birnaviruses and +sRNA viruses, placing them in a special branch of the dsRNA viruses that appears to have evolved independently of the others. The membrane association of their replication complexes is a further specificity of birnaviruses that links them more to the +sRNA viruses than to the other dsRNA viruses. Indeed, within the megataonomy merging together Baltimores' and Linnean classifications proposed by Koonin et al., the Birnaviridae family status is uncertain (Koonin et al., 2020).

## Materials and methods

### Cell lines

Sf9 cell line from Spodoptera frugiperda pupal ovarian tissue were purchased from Gibco, Thermo Fisher Scientific, and used for His-VP3 polypeptides production in the Bac-to-Bac Baculovirus Expression System. Quail muscle fibroblasts (QM7; ATCC CRL-1962) were purchased from the American Type Culture Collection, and used for cell biology assays. The absence of mycoplasma was periodically verified by using the MycoAlert Mycoplasma Detection Assays purchased from Lonza.

### Recombinant prokaryotic-expression plasmid

We used an *Escherichia coli*-based prokaryotic system to obtain the His-2xFYVE and GST-2xFYVE fusion proteins. For the first one, the plasmid encoding EGFP-2xFYVE, pEGFP-2xFYVE, generated by Gillooly and collaborators (Gillooly et al., 2000) was kindly provided by Harald A. Stenmark (Norwegian Radium Hospital, Oslo, Norway) to our group and used as a template. The 2xFYVE fingers consisted of duplicated domains of residues 147–223 of mouse hepatocyte growth factor-regulated tyrosine kinase substrate (Hrs) (Komada and Kitamura, 1995) separated by a QGQGS linker. The prokaryotic vector expressing 2xFYVE was obtained following standard PCR methods with the oligonucleotide pair 2xFYVE (Supplementary file 1). The PCR fragment was inserted into the pET32a(+), which incorporates the thioredoxin A (TrxA) gene, a 6xHis tag, and Thrombin site (Ts) within the amino-terminal domain of the 2xFYVE, and the correctness of the insertion was confirmed by sequencing. The complete fusion protein TrxA.His.Ts.2xFYVE (hereafter named His-2xFYVE) is 954 base pairs
(bp)-long and contains 317 residues (~35 kDa). The plasmid encoding GST-2xFYVE was kindly gifted by Harald A. Stenmark (Norwegian Radium Hospital, Oslo, Norway) to the group of Etienne Morel.

## Recombinant eukaryotic-expression plasmids

We used a *Baculovirus*-based eukaryotic system to obtain the His-VP3 FL, His-VP3 FL R₂₀₀D, and His-VP3 ΔCt fusion proteins. The construct of the recombinant baculovirus (rBV) His-VP3 FL was generated from IBDV VP3 FL protein strain Soroa, GenBank AF140705.1, residues 756–1012 (257 residues) from the polyprotein cloned in pcDNA vector. Thus, IBDV VP3 FL protein was subcloned from plasmid pcDNA VP3 FL into the pFast Bac HTb vector (Gibco BRL), which incorporates a 6xHis tag at the N-terminal of the recombinant protein. Briefly, VP3 FL coding sequence was amplified by PCR with the oligonucleotide pairs VP3.Fw and VP3.Rv incorporating *Bam*HI and *Eco*RI restriction sites, respectively (*Supplementary file 1*). The resulting DNA fragment containing the VP3 gene was digested with *Bam*HI and *Eco*RI and ligated to plasmid pFast Bac HTb restricted with the same enzymes, generating the plasmid pFast Bac-his-VP3 FL. This plasmid contained the VP3 gene fused to a histidine tag under the control of the polyhedrin promoter. The complete fusion protein His-VP3 FL is 861 bp-long and contains 286 residues (~32 kDa). The generation of rBV His-VP3 FL R₂₀₀D was performed following the same procedure, but using the pcDNAVP3 FL R₂₀₀D as a template for subcloning. The complete fusion protein His-VP3 FL R₂₀₀D is 861 bp-long and contains 286 residues (~32 kDa). The generation of rBV His-VP3 Δ223–257 (His-VP3 ΔCt) was performed as follows. A 663 nucleotide DNA fragment containing a truncated version of the VP3 sequence encoding a polypeptide lacking the 36 carbox-terminal residues was amplified by PCR from pcDNAVP3 FL with the oligonucleotide pairs VP3.Fw and VP3 Δ223–257.Rv incorporating *Bam*HI and *Eco*RI restriction sites, respectively (*Supplementary file 1*). The DNA fragment was digested with *Bam*HI and *Eco*RI and was ligated to plasmid pFast Bac HTb restricted with the same enzymes, generating the plasmid pFast Bac-his-VP3 Δ223–257, hereafter named His-VP3 ΔCt. All the resulting plasmids were subjected to nucleotide sequencing to assess the correctness of the inserted VP3 sequence, and they were then used to produce the corresponding rBVs by using the Bac-to-Bac system and by following the manufacturer's instructions (Invitrogen). Briefly, the recombinant plasmids were transformed into DH10Bac competent cells to obtain recombinant bacmids. The white colonies on culture plates were picked, confirmed by PCR, and used for Sf9 cell transfection by using Cellfectin II Reagent (Gibco). The transfected cells were cultured at 28 °C for 5 d to produce rBVs. Then, the cells were subjected to Western blot analysis for recombinant protein expression confirmation and the supernatant (P1) was used to amplify the rBVs to a higher titer P2 and scaled up to obtain a higher volume of P3 to be titrated (Plate-Forme Technologique Production et Purification de Protéines Recombinantes, Institut Pasteur) and used for protein production.

## Protein expression and purification

His-VP3 polypeptides were produced in the Bac-to-Bac Baculovirus Expression System. Briefly, Sf9 cells (Gibco, Thermo Fisher Scientific) were infected with recombinant Baculoviruses (rBVs) at a multiplicity of infection (MOI) of 5 PFU/cell. Cells were harvested at 96 hr post-infection, washed once with chilled phosphate-buffered saline (PBS), resuspended in lysis buffer 50 mM Tris-HCl (pH 8.0), 500 mM NaCl, 0.1% nonidet P-40 supplemented with protease inhibitors (Complete Mini, Roche), and maintained on ice for 30 min. Thereafter, extracts were sonicated and centrifuged at 13.000 g for 15 min at 4 °C. Supernatants were collected and subjected to immobilized metal-affinity chromatography (IMAC) purification by using a Ni²⁺ affinity column (HiTrap Fast Flow Crude, GE Healthcare). Resin-bound proteins were eluted with elution buffer [50 mM Tris-HCl (pH 8.0), 500 mM NaCl, 500 mM imidazole]. His-tagged VP3 FL, VP3 ΔCt, or VP3 FL R₂₀₀D -containing fractions were pooled and subjected to size exclusion chromatography using a Superdex 200 column (GE HealthCare). Finally, protein samples were concentrated to a final concentration of ~3 mg/ml by using Centricon filters (Millipore), aliquoted and flash-frozen in liquid nitrogen for -80 °C storage. The His-VP3 FL, His-VP3 ΔCt, and His-VP3 FL R₂₀₀D fusion proteins obtained from Baculoviruses were used for co-flotation, bio-layer interferometry, and cryo-EM assays.

For His-2xFYVE polypeptide production, *E. coli* BL21 DE3 cells containing the plasmid pET32a(+)–2xFYVE were grown at 37 °C until an OD600 of 0.6 was reached. Protein expression was induced overnight at 18 °C using 1 mM IPTG. Cells were harvested via centrifugation, resuspended in buffer [10 mM TRIS-HCl (pH 8.0), 100 mM NaCl] supplemented with protease inhibitors, and sonicated. The
cell lysate was clarified via centrifugation at 18.000 g for 20 min at 4 °C and the resulting superna-
tant was subjected to IMAC purification by using a Ni²⁺ affinity column (HiTrap Fast Flow Crude, GE
Healthcare). Resin-bound proteins were eluted with elution buffer [10 mM TRIS-HCl (pH 8.0), 100 mM
NaCl, 500 mM imidazole]. His-2xFYVE-containing fractions were pooled and subjected to size exclu-
sion chromatography using a Superdex 75 column (GE HealthCare). Finally, protein samples were
concentrated to a final concentration of ~3 mg/ml by using Centricon filters (Millipore), aliquoted and
flash-frozen in liquid nitrogen for -80 °C storage. The His-2xFYVE was used as a positive control in
co-flotation, bio-layer interferometry, and cryo-EM assays.

For GST-2xFYVE polypeptide production, a classical GST fusion protein purification protocol, using pGEX-5X-1 vector (GE Healthcare, 27-4584-01) cloning as described was followed (Nascimbeni et al., 2017). Briefly, E. coli BL21 DE3 strain maxi-cultures were lysed and subjected to glutathione beads binding for GST purification. The glutathione beads were finally treated with benzamidine-coated beads (GE Healthcare Life Science, 17512301) to remove the factor Xa and dialyzed for the final purification in a 75 mM Kac, 30 mM HEPES, pH 7.4, 5 mM MgCl₂ solution. The GST-2xFYVE fusion protein was used as a probe to detect PI3P in immunofluorescence assays.

Purified proteins quantification and quality control assessment

His-VP3 FL, His-VP3 ΔCt, His-VP3 FL R200D, and His-2xFYVE quantification at 280 nm was carried out by recording a full spectrum between 240 and 340 nm. Measurements were done with 60 µL of buffer and sample at 20 °C in a 1 cm quartz cell, reference 105.202-QS.10 (Hellma Analytics, France), using a JASCO V-750 spectrophotometer (JASCO Corporation, Japan). A baseline subtraction at 340 nm was performed with the Spectragryph software to accurately calculate the protein concentration. The assessment of the sample homogeneity was performed by Dynamic Light Scattering (DLS) analysis on a DynaPro Plate Reader III (Wyatt, Santa Barbara, CA, USA) to ensure that the samples did not contain aggregates. A volume of 20 µL of sample was loaded in a 384-well microplate (Corning ref 3540, New-York, USA), with 10 acquisitions of 5 s each at 20 °C, monitored with the DYNAMICS version V7.10.0.21 software (Wyatt, Santa Barbara, CA, USA), three repetitions per measurement. Finally, the protein integrity and purity were assessed by intact mass spectrometry on a Bruker UltrafeXtreme MALDI-TOF/TOF instrument. A volume of 15 µL of protein was passed through a ZipTip C4 and eluted on a MTP 384 ground steel target plate (Bruker-Daltonics, Germany) with 2 µl of a 20 mg/ mα-Cyano-4-hydroxycinnamic acid in 50% acetonitrile, 0.1% trifluoroacetic acid as matrix solution. Data were acquired using Flexcontrol software (Bruker-Daltonics, Germany) and shots were recorded in positive ion linear mode. Mass spectra were externally calibrated in the m/z range of 15–60 kDa with the Protein II (Bruker-Daltonics, Germany) and analyzed with the flexAnalysis software (Bruker).

Liposome preparation

Liposomes were freshly prepared by the freeze-thaw and extrusion method (F Olson, Olson *et al.*, 1979). Briefly, lipids dissolved in chloroform were mixed and placed in glass vials, and the chloroform was evaporated by centrifuging 2 hr at 35 °C in a SpeedVac. The dried lipids were resuspended in a binding buffer composed of 150 mM NaCl, 20 mM Tris-HCl, pH 8 to a final concentration of 5 mM. The preparation was submitted to ten cycles of freezing in liquid nitrogen and thawing in a 37 °C water bath. To generate small unilamellar liposomes, the multilamellar lipids were extruded through polycarbonate membranes (pore size 200 nm, Whatman). The liposomes were used within 1 wk. The lipids POPE (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine), POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine), PI3P [1,2-dioleoyl-sn-glycero-3-phospho'(1'-myo-inosito'-3'-phosphate)], PA (1,2-dioleoyl-sn-glycero-3-phosphate), and PI [1,2-dioleoyl-sn-glycero-3-phospho'(1'-myo-inosit ol)] were purchased from Avanti Polar Lipids. Liposomes PI3P(-) liposomes were prepared as a mixture of POPE and POPC at a molar ratio of 64:36. Liposomes PI3P(+), PA(+), or PI(+) were prepared as a mixture of POPE and POPC, with either PI3P, PA, or PI at a molar ratio of 64:31:5.

Liposome co-flotation assay

For liposome co-flotation assays, the following proteins were used: His-VP3 FL, His-VP3 ΔCt, His-VP3 FL R200D produced in the Bac-to-Bac Baculovirus Expression System (Invitrogen), and His-2xFYVE produced using the *Escherichia coli* prokaryotic system. Then, 3 mM liposomes were incubated with 1 µM protein in binding buffer composed of 150 mM NaCl, 20 mM Tris-HCl, pH 8 at 4 °C overnight
(on), in a final volume of 100 µl. Liposomes and liposome-protein complexes were separated from unbound protein by ultracentrifugation of the reaction mixture in 5 mL Optiprep (Axis-Shield PoC AS) density gradients for 1 hr at 40,000×g (SW55Ti Beckman Coulter rotor). The liposome-protein mixture was adjusted to 34% Optiprep concentration in a volume of 300 µl, and loaded onto the bottom of the tube that contained 4.5 ml of 20% Optiprep in 150 mM NaCl, 20 mM Tris-HCl, pH 8, overlaid with 200 µl of the buffer. The three top and the bottom fractions, containing 400 µl of the gradient were used for SDS-PAGE gel analysis and immunoblotting against His-tag and anti-VP3 rabbit polyclonal serum.

## Cryo-electron microscopy

Protein-liposome complexes were prepared as described for co-flotation assays. Then, 4 µl of the samples were spotted on glow-discharged lacey grids (S166-3, EMS) and cryo-fixed by plunge freezing at -180 °C in liquid ethane using a Leica EMGP (Leica, Austria). Grids were observed with a Tecnai F20 electron microscope (Thermo Fisher Scientific). The Tecnai F20 was operated at 200 kV and images were acquired under low-dose conditions using the software EPU (Thermo Fisher Scientific) and a direct detector Falcon II (Thermo Fisher Scientific). For immunogold staining, the purified protein were first incubated with 5 µL of Ni-NTA [nickel (II) nitrilotriacetic acid] gold particles 0.5 mM in a five times gold particles/purified protein proportion in binding buffer to a final volume of 10 µl. This reagent comprises 5 nm gold particles covered with multiple Ni-NTA functionalities incorporated into the ligands on the surface of gold particles. Each $Ni^{2+}$ coordinates with one NTA and two histidines from the His-tagged recombinant protein to form a stable complex with extremely low dissociation constants (Reddy et al., 2005). Then, 4 µl of this mix were incubated with 6 µl of 3 mM liposome preparation at 4 °C. Finally, 4 µl of each complex was subjected to cryo-fixation and cryo-EM analyzed as mentioned.

## Bio-layer interferometry (BLI) assays

His-2xFYVE, His-VP3 FL, or His-VP3 $\Delta$Ct were immobilized on Octet NTA biosensor (Sartorius) incubated in an Octet$^{Red}$ 384 system (ForteBio). The functionalized sensors were then washed for 120 s in binding buffer composed of 150 mM NaCl, 20 mM Tris-HCl, pH 8 - BSA (1 mg/mL) before dipping for 600 s into solutions of liposomes diluted in binding buffer - BSA to a final concentration of 5.0, 1.5, or 0.5 mM. Dissociation of the liposome/protein complexes was then monitored in binding buffer - BSA for 300 s. All the experiments were performed in duplicate on two different sensors to account for potential experimental artifacts due to inter-sensor variability. In all cases, potential non-specific interactions were monitored using a sensor coated with the fusion His-Streptavidin as a negative control. The curves were processed using BiaEvaluation software (Biacore).

## Recombinant mammalian expression plasmids

Plasmids encoding EGFP-Rab5 was kindly provided by Philip D. Stahl (Washington University, St. Louis, MO) to our group. Plasmids pcDNA VP3 FL and pcDNA VP3 P2 (P2 all reversed) were constructed by Valli and collaborators in a previous study (Valli et al., 2012) and kindly provided by José F. Rodríguez (CSIC, Madrid, Spain) to our group. The plasmids encoding point mutant VP3 proteins ($K_{157}$D, $R_{159}$D, $H_{198}$D, and $R_{200}$D) were obtained by site-directed mutagenesis. For this purpose, the plasmid pcDNA VP3 FL was used as template DNA in polymerase chain reactions (PCR) together with oligonucleotide primer pairs containing changed nucleotides in order to generate the point mutations (*Supplementary file 1*). PCRs were carried out using Pfu DNA polymerase (PB-L, Argentina) and the following conditions: initial DNA denaturation at 95 °C for 5 min followed by 30 cycles of 95 °C for 1 min, 65 °C for 1 min and 72 °C for 6 min and a final extension at 72 °C for 10 min. The size of the amplified products was corroborated by agarose gel electrophoresis. In order to degrade the plasmid mold, the PCR mixtures were incubated with *Dpn*I endonuclease (Thermo Fisher Scientific, USA) at 37 °C for 2 hr followed by heating at 80 °C for 20 min. Aliquots of those reactions were transformed into competent *E. coli* DH10B cells (Invitrogen, USA). Colonies carrying the mutated plasmids were selected in LB plates containing ampicillin (100 µg/mL). Plasmids were purified using the EasyPure Plasmid MiniPrep Kit (TransGen Biotech Co., Ltd, China), and the introduction of mutations was confirmed by nucleotide sequencing (Macrogen service, South Korea). Alignments and prediction of amino acid sequences were performed using the Aligner and Translator tools from the JustBio website (https://justbio.com/).
Transient transfections and indirect immunofluorescence assay

QM7 cells were grown on coverslips in an M24 multi-well plate for 12 hr to approximately 70% conflu-
ence and then plasmids were transfected employing Lipofectamine 3000 (number L300015; Thermo
Fisher, Argentina) following the manufacturer's recommendations. After 24 hr p.i., the cells were fixed
with 4% paraformaldehyde (PFA) for 15 min at room temperature (RT), washed with PBS (pH 7.4), and
permeabilized with 0.05% saponin in PBS containing 0.2% BSA for 20 min at RT. Then the cells were
incubated with anti-VP3 primary antibodies overnight (ON) at 4 °C in a humidity chamber, and after
extensive washing, cells were incubated with secondary antibodies conjugated with Alexa Fluor Cy3
for 1 hr 30 min, followed by extensive washes in PBS. The coverslips were mounted in Mowiol plus
Hoechst and analyzed by Confocal Laser Scanning Microscopy (CLSM). Images were captured using
an Olympus FluoView TM FV1000 confocal microscope (Olympus, Argentina) with FV10-ASW (version
01.07.00.16) software and processed using the ImageJ software (Schindelin et al., 2012). For quan-
tification of VP3 cellular distribution (Figure 2D), the percentage of cells with vesicular distribution of
the red signal was calculated out of approximately 30 cells per construct and experiment. For quan-
tification of the co-localization of VP3 and either EGFP-Rab5 or EGFP-2xFYVE, the Manders M<sub>2</sub> coef-
ficient was calculated out of approximately 30 cells per construct and experiment (Figure 2—figure
supplement 1, Figure 3A, respectively). The M2 coefficient, which reflects co-localization of signals,
is defined as the ratio of the total intensities of magenta image pixels for which the intensity in the
blue channel is above zero to the total intensity in the magenta channel. JACoP plugin was utilized
to determine M2 (Bolte and Cordelières, 2006). For VP3 puncta co-distributing with EEA1 and GST-
2xFYVE the number of puncta co-distributing for the three signals was manually determined out of
approximately 40 cells per construct and experiment per 200 µm<sup>2</sup> (Figure 3B).

Reverse genetics

For the generation of recombinant IBDV, we used a modification of the reverse genetics system
described by Qi and coworkers (Qi et al., 2007), in which the full-length sequences of the IBDV
segments, flanked by a hammerhead ribozyme at the 5'-end and a hepatitis delta ribozyme at the
3'-end, are expressed under the control of an RNA polymerase II promoter. For its construction, the
complete segments of IBDV Soroa were amplified from viral dsRNA using a sequence-independent
amplification method (Potgieter et al., 2009) and cloned in plasmid pJET1.2 (Thermo Fisher Scien-
tific). The hammerhead ribozymes, specific for each segment, are based in the HamRz-R, described
by Yanay and cols. (Yanai et al., 2006), with 10 nucleotides of sequence complementary to each IBDV
segment stabilizing Stem I. At the 3'-end, we used the sequence of the antigenomic hepatitis delta
ribozyme HDVagrz SC (Ghanem et al., 2012). DNA fragments containing the ribozyme sequences
were constructed by PCR-based gene synthesis with oligonucleotides designed by the Assembly PCR
Oligo Maker software (Rydzanicz et al., 2005). Each PCR-amplified IBDV full-length segment were
cloned, flanked by the corresponding hammerhead ribozyme and HDVagrz SC, between the EcoRI
and NotI sites of plasmid pCAGEN a gift from Connie Cepko; Addgene plasmid # 11160; (Matsuda
and Cepko, 2004 by Gibson assembly Gibson et al., 2009). The correctness of the sequences in
both plasmids, pCAGEN.Hmz.SegA.Hdz (SegA) and pCAGEN.Hmz.SegB.Hdz (SegB), was assessed
by sequencing (Eurofins Scientific). The plasmid pCAGEN.Hmz.SegA.R<sub>200</sub>D.Hdz (SegA.R<sub>200</sub>D) was
purchased to GenScript and the mutation was verified by sequencing of the complete plasmid (Euro-
fins Scientific).

Recombinant viruses foci forming units assay

QM7 cells were grown in M24 multi-well plates for 12 hr to approximately 90-95% confluency and then 800 ng (400 ng of plasmid with each segment) were transfected by using Lipofectamine 2000 (number 11668027; Thermo Fisher, Spain) following the manufacturer's recommendations. When performing single segment transfections (SegA, SegB, and SegA.R200D), 400 ng of the plasmid pCAGIG a gift from Connie Cepko; Addgene plasmid #11159; *Matsuda and Cepko, 2004* was added to 400 ng of those containing the viral segments to reach a total of 800 ng of total DNA per transfection. pCAGIG is a mammalian expression vector that shares the backbone with the pCAGEN vector used for the expression of the viral segments. The following transfections were performed in triplicate: SegA, SegB, SegA.R200D, SegA + SegB, or SegA.R200D+SegB. At 8 hr post-transfection (p.t.) the supernatants were discarded and the monolayers were recovered for further plating on M6 multi-well plates
containing non transfected QM7 cells. For monolayers transfected with SegA, SegB, SegA.R200D, and SegA.R200D+SegB half of the transfected monolayer was plated on fresh cells. For SegA +SegB, which leads to the production of wild-type virus, several dilutions were tested, always in triplicate. To limit virus diffusion and facilitating the foci forming Avicel RC-591 (FMC Biopolymer) was added to the M6 multi-well plates. 72 hr p.i., the monolayers were 16% formaldehyde-fixated and stained with Coomassie R250 (BioRad).

## Statistical analysis

One-way ANOVA followed by a Tukey's HSD (honestly significant difference) test were performed using AnalystSoft Inc, StatPlus:mac - Version v8 (AnalystSoft Inc, 2025). Plotting was done using the ggplot2 package Wickham, 2016 in R, with assistance from the RStudio software (RStudio Team, 2020).

## Bioinformatics analysis

Multiple sequence alignment and the percent identity matrix were performed with Clustal OMEGA (v1.2.4) (Sievers et al., 2011) implemented at EMBL's European Bioinformatics Institute ('EMBL's European Bioinformatics Institute,' n.d.). Accession numbers of the sequences: Infectious bursal disease virus (AF_140705.1); Blotched snakehead virus (YP_052872.1); Victorian trout aquabirnavirus (YP_009255397.1); Tellina virus 2 (YP_010084301.1); Yellowtail ascites virus (NP_690805.1); Infectious pancreatic necrosis virus (NP_047196.1); Tasmanian aquabirnavirus (YP_009177608.1); Lates calcarifer birnavirus (YP_010086267.1); and Tellina virus 1 (YP_009509080.1). Alignment visualization was done with the R ggmsa package (Zhou et al., 2022) with the assistance from the RStudio software (RStudio Team, 2020). Amino acids are colored according to their side-chain chemistry. Protein sequence logos annotation is displayed on top of the amino acid alignment.

The model of VP3 FL from IBDV [strain Soroa, GenBank AF140705.1, residues 756-1012 (257 residues) from the polyprotein] was generated using a local copy of AlphaFold-2 Jumper et al., 2021, installed using the open-source code available at https://github.com/deepmind/alphafold (Zidek and Tomlinson, 2021; accessed on 1 September 2021). Runs were performed on a Centos 7 workstation with 64 CPUs and 4 GeForce RTX 2080 Ti GPUs, using the casp14 preset and including all PDB templates present in the database. The program produces a per-residue confidence metric termed pLDTT on a scale from 0 to 100. pLDTT values higher than 70 reflect reliable models with correct backbone predictions, and those with values higher than 90 correspond to models with both reliable backbone and side-chain orientation predictions. pLDTT values lower than 50 are a strong predictor of disorder, and regions with such values are either unstructured under physiological conditions or only structured as part of a complex.

## Reagents and antibodies

Western blot and confocal laser scanning microscopy (CLSM) analyses were carried out using rabbit anti-VP3 specific sera followed by horseradish peroxidase (HRP)-conjugated anti-rabbit secondary antibodies (number A0545), purchased from Sigma-Aldrich (Buenos Aires, Argentina). For early endosomes detection by CLSM, we used anti-EEA1 mouse monoclonal antibody (1:200 dilution; BD Transduction Laboratories, 610456), and for anti-GST detection, we used Alexa Fluor 649-labeled goat anti-GST (Rockland-Inc, 600-143-200). Fluorescently labeled secondary antibodies, all diluted 1:200 from Invitrogen, included Alexa Fluor 546 donkey anti-mouse (A10036), and Alexa Fluor 546 donkey anti-rabbit (A10040).

## Far-UV circular dichroism (CD) of VP3 proteins

His-VP3 FL and His-VP3 FL R200D at a concentration of 2 mg/ml in 50 mM Tris-HCl pH 8, 500 mM NaCl were analyzed by far-UV CD (180-260 nm). CD spectra were obtained with an AVIV CD spectropolarimeter model 215 using a 0.02 cm path length cell at RT. Five successive scans were averaged and the background spectrum of the sample buffer, acquired under identical conditions, was subtracted. The resulting corrected CD intensities were then converted to Δε per residue. Secondary structure contents were estimated from the far-UV CD spectra using the CDSSTR routine (Johnson, 1999) of the DICHROWEB server (Whitmore and Wallace, 2008; Whitmore and Wallace, 2004) run on the SP175 reference dataset (Lees et al., 2006), containing 72 proteins representing a large panel
of secondary structures. Similar results were obtained on different datasets (Sreerama and Woody, 2000) or by using the CONTIN/LL routine (Provencher and Glöckner, 1981).

## Coarse grain (CG) models of VP3 constructs

Three main protein constructs derived from VP3 were considered: 1. VP3 Ct. A peptide comprising the last 36 residues of VP3 (223–257). VP3 Ct has an overall electrical charge of +5, and is absent from the available crystal structure of VP3. Consequently, it was assumed to be a disordered domain and was modeled as unstructured. 2. VP3 ΔNt. A shorter version of VP3 that lacks the first 81 residues leading to the Nt. It consists of the crystallized structure of the protein core (residues 82–222, available from the Protein Data Bank PDB ID 2R18), concatenated with the VP3 Ct peptide. As the first nine residues of the protein core are also missing from resolved structure, they were modeled as unstructured, similar to the VP3 Ct fragment. VP3 ΔNt was selected primarily because it contains the structural information of the crystallized protein, which provides valuable insights into the protein's overall structure and function. In addition, the construct includes the Ct peptide (VP3 Ct), which is crucial for membrane binding. Although the Nt 81 residues are absent in this construct, their exclusion was found to be inconsequential in terms of evaluating the protein's binding energy and positioning on the membrane surface (compare, for example, VP3 ΔNt and VP3 FL in Figure 5). 3. VP3 FL. The full-length VP3 protein, which was obtained as an AlphaFold2 prediction (Jumper et al., 2021). The structure of the protein predicted by AlphaFold2 not only provided insight into the overall architecture of VP3, but it also reinforced our model of unstructured VP3 Ct. It is worth noting that the Nt 81 residues, which are absent in VP3 ΔNt, are predicted to be in a modular conformation, constituting a globular domain separated by a loop from the globular domain of the core region (crystallized residues 82–222). The structural information provided by the full-length VP3 model is, therefore, useful in understanding the interaction with lipid membranes. In addition to VP3 Ct, VP3 ΔNt, and VP3 FL, two mutants of VP3 ΔNt were simulated: VP3 ΔNt R200D, in which R200 was replaced by aspartate (D), and VP3 ΔNt P2 Mut, where the whole P2 region (R200, H198, R159, and K157) was replaced by D. All proteins were modeled with the MARTINI force-field (de Jong et al., 2013) [whenever needed, the secondary structure of the non-disordered protein domains was restrained by means of the ElNeDyn elastic network model (elnedyn22)] (Periole et al., 2009; Poma et al., 2017). In this case, harmonic potentials with a force constant of 1000 kJ mol⁻¹ nm² were applied between backbone particles separated by less than 1.4 nm. In the construction of VP3 ΔNt, elastic network restrictions were applied to residues 91–220, which was the fragment resolved from by crystallography. However, in VP3 FL, the elastic network was applied to all globular domains, while leaving the linkers and the last 36 residues at the Ct unrestricted.

## Molecular theory calculations

The molecular theory (MT) used to assess protein adsorption is detailed in references (Chiarpotti et al., 2021; Ramírez et al., 2019). Briefly, this mean-field method yields the distribution of macro-molecules in an anisotropic system, through the minimization of a model free-energy functional that accounts for electrostatic and steric interactions, the configurational and translational entropy of the macromolecules, the translational entropy of the free species (protons, hydroxyls, and salt ions) and the acid-base equilibrium of all ionizable molecules in the system. In the present case, the membrane was modeled as a dielectric slab of thickness h, with a dielectric constant ε<sub>M</sub>, which is in contact with a solution of dielectric constant ε<sub>S</sub>. Far from the membrane (z→∞), the solution is in contact with an infinite bulk that contains the proteins (1 µM concentration), all the free species, and a fixed pH = 8. The rigid membrane contained 15% of ionizable groups on its surface, in order to mimic the 5% of ternary ionizable PI3P. One-third of these groups were assigned a pK<sub>a</sub> of 2.5, corresponding to the phosphate group of the glycerol moiety, while the pK<sub>a</sub> for the rest was set to 6.5, corresponding to the phosphoether groups of inositol. MT requires the input of a large set of protein conformations, which are subsequently reweighted by the theory. These conformations were generated by Molecular Dynamics simulations of the proteins in water. VP3 Ct, VP3 ΔNt, VP3 FL, VP3 ΔNt R200D, and VP3 ΔNt P2 Mut, were modeled using MARTINI-2.2 parameters (Marrink et al., 2007). The proteins were hydrated in MARTINI polarizable water. As stated before, the proteins globular segments (residues 92–222 in VP3 FL, VP3 ΔNt and its mutants, and residues 7–79 additionally in VP3 FL) were restricted through the application of the elastic network ElNeDyn (Marrink et al., 2007) with a spring constant
of 1000 kJ.mol⁻¹ applied to backbone beads (BB) within a cutoff distance of 1.4 nm. Each protein was simulated for 100 ns, after equilibration, in the NPT ensemble using GROMACS 2021.3. A cluster analysis was performed in order to obtain the most representative conformations of the simulated ensemble. This method groups conformations by structural similarity, using the root mean squared deviation (RMSD) of the coordinates calculated for all conformer pairs along the trajectory, clustering structures with RMSD ≤0.3 nm. This procedure leads to 50 conformers for each protein construct. Each conformer was subsequently rotated 30 times randomly to produce a set of 1500 configurations for the MT calculations.

## Molecular dynamics simulations of VP3 on a lipid bilayer

MARTINI-2.2 was also used to model the lipid molecules in our MD simulations. The bilayer was assembled using an insane utility (Wassenaar et al., 2015), while hydration with polarizable water molecules and ions was performed with the solvate utility of the GROMACS package. The bilayer had a symmetric lipid composition of 30% 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), 65% of 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE), and 5% 1,2-dioleoyl-sn-glycero-3-phosphatidyli nositol-3-phosphate (PI3P), leading to a lipid composition of 111:231:18 (DOPC, DOPE, and PI3P) in each leaflet. The protein was placed at a distance of 10 nm from the membrane surface, to avoid a possible initial binding bias. The system was then equilibrated and finally simulated for 500 ns. Both equilibration and production were in the NPT ensemble. A semi-isotropic Parinello-Rahman barostat (Parrinello and Rahman, 1981) with a 12 ps time constant, and lateral and normal compressibilities of 3×10⁻⁴ bar⁻¹ were used to keep the pressure at 1.0 bar. The v-rescale thermostat, with a time constant of 1 ps, was used to keep the temperature at 325 K. In all cases, the integration timestep was 20 fs. Long range electrostatics were calculated using the particle mesh Ewald (PME) method (Darden et al., 1993) with a real-space cut-off of 1.2 nm.

## Acknowledgements

This project was partially supported by: the National Agency for Scientific, and Technological Promotion, Ministry of Science, Technology and Innovation, through grants PICT 2016-0528 and 2019-01324 to LRD; PICT 2015-2210 to FAZ; PICT 2019-01889 to LMP, and PICT 2020-3795 to MGDP; the National Scientific and Technical Research Council through grants PIP 2015-2017 11220150100114CO and 2021-2023 11220200103139CO to LRD; 2021-2023 11220200103195CO to LMP and 11220200103223CO to MGDP; and by the National University of Cuyo through grants 2013-2015 M006, 2016-2018 M029, 2019-2021 M071 and 2022-2024 M012 to LRD. For stays performed by LRD in the laboratory of Structural Virology at the Institut Pasteur in Paris, we are grateful for the support of National University of Cuyo through the Teaching Staff Mobility Program in 2018-2019, the Ministry of Education through the Program scholarships for training abroad in Science and Technology (BEC.AR 2018), and the International Union of Biochemistry and Molecular Biology (IUBMB) through a '2020 Mid-Career Research Fellowship.' The work in Madrid was supported by grants from the Spanish Ministry of Science and Innovation (PID2020-113287RB-100) and the Comunidad Autónoma de Madrid (P2018/NMT-4389) to JRC. At the IHEM, we sincerely appreciate Elisa Bocanegra, Norberto Domizio, and Jorge Ibañez for valuable technical assistance in CLSM handling. At the Institut Pasteur in Paris, we are grateful to Gérard Pehau-Arnaudet and the staff at the Ultrastuctural BioImaging core facility for image acquisition and analysis; to Patrick England and Sébastien Brûlé and the staff of the Molecular Biophysics platform for purified protein analyses, CD and BLI assays, to Stéphane Petres and members at the Production and Purification of Recombinant Proteins Technological Platform, and to the staff of the Proteomics platform.

## Additional information

### Competing interests

Laura Ruth Delgui: Reviewing editor, eLife. The other authors declare that no competing interests exist.
## Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>Ministry of Science, Technology and Innovation</td><td>PICT 2016-0528</td><td>Laura Ruth Delgui</td></tr><tr><td>Ministry of Science, Technology and Innovation</td><td>PICT 2019-01324</td><td>Laura Ruth Delgui</td></tr><tr><td>Ministry of Science, Technology and Innovation</td><td>PICT 2015-2210</td><td>Flavia A Zanetti</td></tr><tr><td>Ministry of Science, Technology and Innovation</td><td>PICT 2019-01889</td><td>Luis Mariano Polo</td></tr><tr><td>Ministry of Science, Technology and Innovation</td><td>PICT 2020-3795</td><td>Mario Del Pópolo</td></tr><tr><td>National Scientific and Technical Research Council</td><td>PIP 2015-2017<br>11220150100114CO</td><td>Laura Ruth Delgui</td></tr><tr><td>National Scientific and Technical Research Council</td><td>PIP 2021-2023<br>11220200103195CO</td><td>Laura Ruth Delgui</td></tr><tr><td>National Scientific and Technical Research Council</td><td>PIP 11220200103223CO</td><td>Mario Del Pópolo</td></tr><tr><td>National Scientific and Technical Research Council</td><td>2021-2023<br>11220200103139CO</td><td>Laura Ruth Delgui</td></tr><tr><td>National University of Cuyo</td><td>2013-2015 M006</td><td>Laura Ruth Delgui</td></tr><tr><td>National University of Cuyo</td><td>2016-2018 M029</td><td>Laura Ruth Delgui</td></tr><tr><td>National University of Cuyo</td><td>2019-2021 M071</td><td>Laura Ruth Delgui</td></tr><tr><td>National University of Cuyo</td><td>2022-2024 M012</td><td>Laura Ruth Delgui</td></tr><tr><td>National University of Cuyo</td><td>Teaching Staff Mobility Program in 2018-2019</td><td>Laura Ruth Delgui</td></tr><tr><td>Ministry of Education</td><td>BEC.AR 2018</td><td>Laura Ruth Delgui</td></tr><tr><td>International Union of Biochemistry and Molecular Biology</td><td>2020 Mid-Career Research Fellowship</td><td>Laura Ruth Delgui</td></tr><tr><td>Spanish Ministry of Science and Innovation</td><td>PID2020-113287RB-I00</td><td>Jose R Castón</td></tr><tr><td>Comunidad Autónoma de Madrid</td><td>P2018/NMT-4389</td><td>Jose R Castón</td></tr></tbody></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

## Author contributions

Flavia A Zanetti, Pablo Guardado-Calvo, Andres Ferrino-Iriarte, Sarah Dubois, Etienne Morel, Victoria Alfonso, Milton Osmar Aguilera, María E Celayes, Luis Mariano Polo, Laila Suhaiman, Vanesa V Galassi, Maria V Chiarpotti, Carolina Allende-Ballesteros, Javier M Rodriguez, Oscar Taboga, Investigation; Ignacio Fernandez, Conceptualization, Formal analysis; Eduard Baquero, Conceptualization; Jose R Castón, Resources; Diego Lijavetzky, Data curation, Investigation; María I Colombo, Resources, Writing – review and editing; Mario Del Pópolo, Conceptualization, Data curation, Writing – original draft, Writing – review and editing; Félix A Rey, Data curation, Writing – review and editing; Laura Ruth Delgui, Conceptualization, Data curation, Formal analysis, Supervision, Funding acquisition, Investigation, Methodology, Writing – original draft, Writing – review and editing

## Author ORCIDs

Sarah Dubois https://orcid.org/0009-0000-7504-7135

Etienne Morel https://orcid.org/0000-0002-4763-4954

Javier M Rodriguez https://orcid.org/0000-0003-0146-9903

Jose R Castón https://orcid.org/0000-0003-2350-9048
Mario Del Pópolo https://orcid.org/0000-0002-1435-2424

Félix A Rey https://orcid.org/0000-0002-9953-7988

Laura Ruth Delgui https://orcid.org/0000-0002-3647-3593

## Peer review material

Reviewer #1 (Public review): https://doi.org/10.7554/eLife.97261.3.sa1

Reviewer #3 (Public review): https://doi.org/10.7554/eLife.97261.3.sa2

Author response https://doi.org/10.7554/eLife.97261.3.sa3

## Additional files

### Supplementary files

Supplementary file 1. Primers. All the primers used for this work. * For the first four-point mutants, the nucleotides that allow to introduce amino acid changes in the VP3 protein are indicated in italic bold letters. The underlined nucleotides indicate restriction sites.

MDAR checklist

### Data availability

All data generated or analysed during this study are included in the manuscript and supporting files.

The following previously published datasets were used:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset URL</td><td>Database and Identifier</td></tr></thead><tbody><tr><td>Da Costa B, Soignier S, Chevalier C, Henry C, Thory C, Huet JC, Delmas B</td><td>2018</td><td>Blotched snakehead virus is a new aquatic birnavirus that is slightly more related to avibirnavirus than to aquabirnavirus</td><td>https://www.ncbi.nlm.nih.gov/protein/YP_052872.1</td><td>NCBI Protein, YP_052872.1</td></tr><tr><td>Mohr PG, Moody NJ, Williams LM, Hoad J, Crane MS</td><td>2018</td><td>Molecular characterization of Tasmanian aquabirnaviruses from 1998 to 2013</td><td>https://www.ncbi.nlm.nih.gov/protein/YP_009255397.1</td><td>NCBI Protein, YP_009255397.1</td></tr><tr><td>Blake S, Ma JY, Caporale DA, Jairath S, Nicholson BL</td><td>2021</td><td>Phylogenetic relationships of aquatic birnaviruses based on deduced amino acid sequences of genome segment A cDNA</td><td>https://www.ncbi.nlm.nih.gov/protein/YP_010084301.1</td><td>NCBI Protein, YP_010084301.1</td></tr><tr><td>Suzuki S, Kimura M, Kusuda R</td><td>2018</td><td>The complete nucleotide sequence of the polyprotein and VP5 gene of a marine birnavirus</td><td>https://www.ncbi.nlm.nih.gov/protein/YP_690805.1</td><td>NCBI Protein, NP_690805.1</td></tr><tr><td>Duncan R, Dobos P</td><td>2018</td><td>The nucleotide sequence of infectious pancreatic necrosis virus (IPNV) dsRNA segment A reveals one large ORF encoding a precursor polyprotein</td><td>https://www.ncbi.nlm.nih.gov/protein/YP_047196.1</td><td>NCBI Protein, NP_047196.1</td></tr><tr><td>Mohr PG, Moody NJ, Williams LM, Hoad J, Crane MS</td><td>2018</td><td>Molecular characterization of Tasmanian aquabirnaviruses from 1998 to 2013</td><td>https://www.ncbi.nlm.nih.gov/protein/YP_009177608.1</td><td>NCBI Protein, YP_009177608.1</td></tr><tr><td>Nobiron I, Galloux M, Henry C, Torhy C, Boudinot P, Lejal N, Da Costa B, Delmas B</td><td>2018</td><td>Genome and polypeptides characterization of Tellina virus 1 reveals a fifth genetic cluster in the Birnaviridae family</td><td>https://www.ncbi.nlm.nih.gov/protein/YP_009509080.1</td><td>NCBI Protein, YP_009509080.1</td></tr><tr><td>Chen J, Toh X, Ong J, Wang Y, Teo X, Lee B, Wong P, Chee S, Wee D, Wang A, Ng Y, Tan BM, Khor D, Chong S</td><td>2022</td><td>Detection and characterisation of a novel marine Birnavirus isolated from Asian Seabass in Singapore</td><td>https://www.ncbi.nlm.nih.gov/protein/YP_010086267.1</td><td>NCBI Protein, YP_010086267.1</td></tr></tbody></table>
## References

* AnalystSoft Inc. 2025. StatPlus:mac. v8. StatPlus.

* Baltimore D. 1971. Expression of animal virus genomes. *Bacteriological Reviews* **35**:235-241. DOI: https://doi.org/10.1128/br.35.3.235-241.1971, PMID: 4329869

* Banani SF, Lee HO, Hyman AA, Rosen MK. 2017. Biomolecular condensates: organizers of cellular biochemistry. *Nature Reviews. Molecular Cell Biology* **18**:285-298. DOI: https://doi.org/10.1038/nrm.2017.7, PMID: 28225081

* Bolte S, Cordelières FP. 2006. A guided tour into subcellular colocalization analysis in light microscopy. *Journal of Microscopy* **224**:213-232. DOI: https://doi.org/10.1111/j.1365-2818.2006.01706.x, PMID: 17210054

* Brodrick AJ, Broadbent AJ. 2023. The formation and function of Birnaviridae virus factories. *International Journal of Molecular Sciences* **24**:8471. DOI: https://doi.org/10.3390/ijms24108471, PMID: 37239817

* Burd CG, Emr SD. 1998. Phosphatidylinositol(3)-phosphate signaling mediated by specific binding to RING FYVE domains. *Molecular Cell* **2**:157-162. DOI: https://doi.org/10.1016/s1097-2765(00)80125-2, PMID: 9702203

* Casañas A, Navarro A, Ferrer-Orta C, González D, Rodríguez JF, Verdaguer N. 2008. Structural insights into the multifunctional protein VP3 of birnaviruses. *Structure* **16**:29-37. DOI: https://doi.org/10.1016/j.str.2007.10.023, PMID: 18184581

* Chiarpotti MV, Longo GS, Del Pópolo MG. 2021. Nanoparticles modified with cell penetrating peptides: assessing adsorption on membranes containing acidic lipids. *Colloids and Surfaces. B, Biointerfaces* **197**:111373. DOI: https://doi.org/10.1016/j.col surfb.2020.111373, PMID: 33045543

* Coulibaly F, Chevalier C, Gutsche I, Pous J, Navaza J, Bressanelli S, Delmas B, Rey FA. 2005. The birnavirus crystal structure reveals structural relationships among icosahedral viruses. *Cell* **120**:761-772. DOI: https://doi.org/10.1016/j.cell.2005.01.009, PMID: 15797378

* Darden T, York D, Pedersen L. 1993. Particle mesh Ewald: An N-log(N) method for Ewald sums in large systems. *The Journal of Chemical Physics* **98**:10089-10092. DOI: https://doi.org/10.1063/1.464397

* de Jong DH, Singh G, Bennett WFD, Arnarez C, Wassenaar TA, Schäfer LV, Periole X, Tieleman DP, Marrink SJ. 2013. Improved parameters for the martini coarse-grained protein force field. *Journal of Chemical Theory and Computation* **9**:687-697. DOI: https://doi.org/10.1021/ct300646g, PMID: 26589065

* Delgui LR, Rodríguez JF, Colombo MI. 2013. The endosomal pathway and the Golgi complex are involved in the infectious bursal disease virus life cycle. *Journal of Virology* **87**:8993-9007. DOI: https://doi.org/10.1128/JVI.03152-12, PMID: 23741000

* Delmas B, Attoui H, Ghosh S, Malik YS, Mundt E, Vakharia VN. 2019. ICTV virus taxonomy profile: birnaviridae. *The Journal of General Virology* **100**:5-6. DOI: https://doi.org/10.1099/jgv.0.001185, PMID: 30484762

* Dolinsky TJ, Nielsen JE, McCammon JA, Baker NA. 2004. PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations. *Nucleic Acids Research* **32**:W665-W667. DOI: https://doi.org/10.1093/nar/gkh381, PMID: 15215472

* Fernandes F, Chen K, Ehrlich LS, Jin J, Chen MH, Medina GN, Symons M, Montelaro R, Donaldson J, Tjandra N, Carter CA. 2011. Phosphoinositides direct equine infectious anemia virus gag trafficking and release. *Traffic* **12**:438-451. DOI: https://doi.org/10.1111/j.1600-0854.2010.01153.x, PMID: 21176037

* Gaullier JM, Simonsen A, D'Arrigo A, Bremnes B, Stenmark H, Aasland R. 1998. FYVE fingers bind PtdIns(3)P. *Nature* **394**:432-433. DOI: https://doi.org/10.1038/28767, PMID: 9697764

* Ghanem A, Kern A, Conzelmann KK. 2012. Significantly improved rescue of rabies virus from cDNA plasmids. *European Journal of Cell Biology* **91**:10-16. DOI: https://doi.org/10.1016/j.eijcb.2011.01.008, PMID: 21397981

* Gibson DG, Young L, Chuang R-Y, Venter JC, Hutchison CA III, Smith HO. 2009. Enzymatic assembly of DNA molecules up to several hundred kilobases. *Nature Methods* **6**:343-345. DOI: https://doi.org/10.1038/nmeth.1318

* Gillooly DJ, Morrow IC, Lindsay M, Gould R, Bryant NJ, Gaullier JM, Parton RG, Stenmark H. 2000. Localization of phosphatidylinositol 3-phosphate in yeast and mammalian cells. *The EMBO Journal* **19**:4577-4588. DOI: https://doi.org/10.1093/emboj/19.17.4577, PMID: 10970851

* Gimenez MC, Zanetti FA, Terebiznik MR, Colombo MI, Delgui LR. 2018. Infectious bursal disease virus hijacks endosomal membranes as the scaffolding structure for viral replication. *Journal of Virology* **92**:e01964-17. DOI: https://doi.org/10.1128/JVI.01964-17, PMID: 29540593

* Gimenez MC, Issa M, Sheth J, Colombo MI, Terebiznik MR, Delgui LR. 2021. Phosphatidylinositol 3-phosphate mediates the establishment of infectious bursal disease virus replication complexes in association with early endosomes. *Journal of Virology* **95**:e02313-20. DOI: https://doi.org/10.1128/JVI.02313-20

* Gorbalenya AE, Pringle FM, Zeddam JL, Luke BT, Cameron CE, Kalmakoff J, Hanzlik TN, Gordon KHJ, Ward VK. 2002. The palm subdomain-based active site is internally permuted in viral RNA-dependent RNA polymerases of an ancient lineage. *Journal of Molecular Biology* **324**:47-62. DOI: https://doi.org/10.1016/s0022-2836(02)01033-1, PMID: 12421558

* Gorvel JP, Chavrier P, Zerial M, Gruenberg J. 1991. rab5 controls early endosome fusion in vitro. *Cell* **64**:915-925. DOI: https://doi.org/10.1016/0092-8674(91)90316-q, PMID: 1900457

* Johnson WC. 1999. Analyzing protein circular dichroism spectra for accurate secondary structures. *Proteins* **35**:307-312. DOI: https://doi.org/10.1002/(SICI)1097-0134(19990515)35:3, PMID: 10328265

* Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, Tunyasuvunakool K, Bates R, Židek A, Potapenko A, Bridgland A, Meyer C, Kohl SAA, Ballard AJ, Cowie A, Romera-Paredes B, Nikolov S, Jain R, Adler J, Back T, et al. 2021. Highly accurate protein structure prediction with AlphaFold. *Nature* **596**:583-589. DOI: https://doi.org/10.1038/s41586-021-03819-2, PMID: 34265844
Kairys V, Baranauskiene L, Kazlauskiene M, Matulis D, Kazlauskas E. 2019. Binding affinity in drug design: experimental and computational techniques. *Expert Opinion on Drug Discovery* **14**:755–768. DOI: https://doi.org/10.1080/17460441.2019.1623202, PMID: 31146609

Kolli S, Meng X, Wu X, Shengjuler D, Cameron CE, Xiang Y, Deng J. 2015. Structure-function analysis of vaccinia virus h7 protein reveals a novel phosphoinositide binding fold essential for poxvirus replication. *Journal of Virology* **89**:2209–2219. DOI: https://doi.org/10.1128/JVI.03073-14/ASSET/F1575F5D-784E-4FD0-9DDB-A2B374935CCA/ASSETS/GRAPHIC/ZJV9990900510007.JPEG

Komada M, Kitamura N. 1995. Growth factor-induced tyrosine phosphorylation of hrs, a novel 115-kilodalton protein with a structurally conserved putative zinc finger domain. *Molecular and Cellular Biology* **15**:6213–6221. DOI: https://doi.org/10.1128/MCB.15.11.6213

Koonin EV, Dolja VV, Krupovic M, Varsani A, Wolf YI, Yutin N, Zerbini FM, Kuhn JH. 2020. Global organization and proposed megataonomy of the virus world. *Microbiology and Molecular Biology Reviews* **84**:e00061-19. DOI: https://doi.org/10.1128/MMBR.00061-19, PMID: 32132243

Lees JG, Miles AJ, Wien F, Wallace BA. 2006. A reference database for circular dichroism spectroscopy covering fold and secondary structure space. *Bioinformatics* **22**:1955–1962. DOI: https://doi.org/10.1093/bioinformatics/btl327, PMID: 16787970

Lombardo E, Maraver A, Castón JR, Rivera J, Fernández-Arias A, Serrano A, Carrascosa JL, Rodríguez JF. 1999. VP1, the putative RNA-dependent RNA polymerase of infectious bursal disease virus, forms complexes with the capsid protein VP3, leading to efficient encapsidation into virus-like particles. *Journal of Virology* **73**:6973–6983. DOI: https://doi.org/10.1128/JVI.73.8.6973-6983.1999

Luque D, Rivas G, Alfonso C, Carrascosa JL, Rodríguez JF, Castón JR. 2009a. Infectious bursal disease virus is an icosahedral polyploid dsRNA virus. *PNAS* **106**:2148–2152. DOI: https://doi.org/10.1073/pnas.0808498106, PMID: 19164552

Luque D, Saugar I, Rejas MT, Carrascosa JL, Rodríguez JF, Castón JR. 2009b. Infectious Bursal disease virus: ribonucleoprotein complexes of a double-stranded RNA virus. *Journal of Molecular Biology* **386**:891–901. DOI: https://doi.org/10.1016/j.jmb.2008.11.029, PMID: 19063900

Maraver A, Oña A, Abaitua F, González D, Clemente R, Ruiz-Díaz JA, Castón JR, Pazos F, Rodríguez JF. 2003. The oligomerization domain of VP3, the scaffolding protein of infectious bursal disease virus, plays a critical role in capsid assembly. *Journal of Virology* **77**:6438–6449. DOI: https://doi.org/10.1128/JVI.77.11.6438-6449.2003

Marrink SJ, Risselada HJ, Yefimov S, Tieleman DP, de Vries AH. 2007. The MARTINI force field: coarse grained model for biomolecular simulations. *The Journal of Physical Chemistry. B* **111**:7812–7824. DOI: https://doi.org/10.1021/jp071097f, PMID: 17569554

Matsuda T, Cepko CL. 2004. Electroporation and RNA interference in the rodent retina in vivo and in vitro. *PNAS* **101**:16–22. DOI: https://doi.org/10.1073/pnas.2235688100, PMID: 14603031

McNulty S, Bornmann W, Schriewer J, Werner C, Smith SK, Olson VA, Damon IK, Buller RM, Heuser J, Kalman D. 2010. Multiple phosphatidylinositol 3-kinases regulate vaccinia virus morphogenesis. *PLOS ONE* **5**:e10884. DOI: https://doi.org/10.1371/journal.pone.0010884, PMID: 20526370

Miles AJ, Wallace BA. 2018. CDtoolX, a downloadable software package for processing and analyses of circular dichroism spectroscopic data. *Protein Science* **27**:1717–1722. DOI: https://doi.org/10.1002/pro.3474, PMID: 30168221

Nascimbeni AC, Codogno P, Morel E. 2017. Local detection of PtdIns3P at autophagosome biogenesis membrane platforms. *Autophagy* **13**:1602–1612. DOI: https://doi.org/10.1080/15548627.2017.1341465, PMID: 28813193

Olson F, Hunt CA, Szoka FC, Vail WJ, Papahadjopoulos D. 1979. Preparation of liposomes of defined size distribution by extrusion through polycarbonate membranes. *Biochimica et Biophysica Acta* **557**:9–23. DOI: https://doi.org/10.1016/0005-2736(79)90085-3, PMID: 95096

Parrinello M, Rahman A. 1981. Polymorphic transitions in single crystals: a new molecular dynamics method. *Journal of Applied Physics* **52**:7182–7190. DOI: https://doi.org/10.1063/1.328693

Periole X, Cavalli M, Marrink SJ, Ceruso MA. 2009. Combining an elastic network with a coarse-grained molecular force field: structure, dynamics, and intermolecular recognition. *Journal of Chemical Theory and Computation* **5**:2531–2543. DOI: https://doi.org/10.1021/CT9002114/SUPPL_FILE/CT9002114_SI_001.PDF

Poma AB, Cieplak M, Theodorakis PE. 2017. Combining the MARTINI and structure-based coarse-grained approaches for the molecular dynamics studies of conformational transitions in proteins. *Journal of Chemical Theory and Computation* **13**:1366–1374. DOI: https://doi.org/10.1021/ACS.JCTC.6B00986/ASSET/IMAGES/LARGE/CT-2016-00986K_0005.JPEG

Potgieter AC, Page NA, Liebenberg J, Wright IM, Landt O, van Dijk AA. 2009. Improved strategies for sequence-independent amplification and sequencing of viral double-stranded RNA genomes. *The Journal of General Virology* **90**:1423–1432. DOI: https://doi.org/10.1099/vir.0.009381-0, PMID: 19264638

Provencher SW, Glöckner J. 1981. Estimation of globular protein secondary structure from circular dichroism. *Biochemistry* **20**:33–37. DOI: https://doi.org/10.1021/bi00504a006, PMID: 7470476

Puffer BA, Watkins SC, Montelaro RC. 1998. Equine infectious anemia virus Gag polyprotein late domain specifically recruits cellular AP-2 adapter protein complexes during virion assembly. *Journal of Virology* **72**:10218–10221. DOI: https://doi.org/10.1128/JVI.72.12.10218-10221.1998, PMID: 9811764

Qi X, Gao Y, Gao H, Deng X, Bu Z, Wang X, Fu C, Wang X. 2007. An improved method for infectious bursal disease virus rescue using RNA polymerase II system. *Journal of Virological Methods* **142**:81–88. DOI: https://doi.org/10.1016/j.jviromet.2007.01.021, PMID: 17383019
Ramírez PG, Del Pópolo MG, Vila JA, Szleifer I, Longo GS. 2019. Adsorption and insertion of polyarginine peptides into membrane pores: The trade-off between electrostatics, acid-base chemistry and pore formation energy. *Journal of Colloid and Interface Science* **552**:701–711. DOI: https://doi.org/10.1016/j.jcis.2019.05.087, PMID: 31176053

Reddy V, Lymar E, Hu M, Hainfeld JF. 2005. 5 nm gold-Ni-NTA binds his tags. *Microscopy and Microanalysis* **11**:1118–1119. DOI: https://doi.org/10.1017/S1431927605507712

Reddy V, Campbell EA, Wells J, Simpson J, Nazki S, Hawes PC, Broadbent AJ. 2022. Birnaviridae virus factories show features of liquid-liquid phase separation and are distinct from paracrystalline arrays of virions observed by electron microscopy. *Journal of Virology* **96**:e0202421. DOI: https://doi.org/10.1128/jvi.02024-21, PMID: 35138130

RStudio Team. 2020. RStudio: integrated development environment for R. RStudio. https://www.scirp.org/reference/referencespapers?referenceid=3542223

Rydzanicz R, Zhao XS, Johnson PE. 2005. Assembly PCR oligo maker: a tool for designing oligodeoxynucleotides for constructing long DNA molecules for RNA production. *Nucleic Acids Research* **33**:W521-5. DOI: https://doi.org/10.1093/nar/gki380, PMID: 15980526

Saugar I, Irigoyen N, Luque D, Carrascosa JL, Rodríguez JF, Castón JR. 2010. Electrostatic interactions between capsid and scaffolding proteins mediate the structural polymorphism of a double-stranded RNA virus. *The Journal of Biological Chemistry* **285**:3643–3650. DOI: https://doi.org/10.1074/jbc.M109.075994, PMID: 19933276

Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, Preibisch S, Rueden C, Saalfeld S, Schmid B, Tinevez JY, White DJ, Hartenstein V, Eliceiri K, Tomancak P, Cardona A. 2012. Fiji: an open-source platform for biological-image analysis. *Nature Methods* **9**:676–682. DOI: https://doi.org/10.1038/nmeth.2019, PMID: 22743772

Shah NB, Duncan TM. 2014. Bio-layer interferometry for measuring kinetics of protein-protein interactions and allosteric ligand effects. *Journal of Visualized Experiments* **01**:e51383. DOI: https://doi.org/10.3791/51383, PMID: 24638157

Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Söding J, Thompson JD, Higgins DG. 2011. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. *Molecular Systems Biology* **7**:539. DOI: https://doi.org/10.1038/msb.2011.75, PMID: 21988835

Simonsen A, Lippé R, Christoforidis S, Gaullier JM, Brech A, Callaghan J, Toh BH, Murphy C, Zeral M, Stenmark H. 1998. EEA1 links PI(3)K function to Rab5 regulation of endosome fusion. *Nature* **394**:494–498. DOI: https://doi.org/10.1038/28879, PMID: 9697774

Sreerama N, Woody RW. 2000. Estimation of protein secondary structure from circular dichroism spectra: comparison of CONTIN, SELCON, and CDSSTR methods with an expanded reference set. *Analytical Biochemistry* **287**:252–260. DOI: https://doi.org/10.1006/abio.2000.4880, PMID: 11112271

Stenmark H, Aasland R, Toh BH, D'Arrigo A. 1996. Endosomal localization of the autoantigen EEA1 is mediated by a zinc-binding FYVE finger. *The Journal of Biological Chemistry* **271**:24048–24054. DOI: https://doi.org/10.1074/jbc.271.39.24048, PMID: 8798641

Tanzi GO, Piefer AJ, Bates P. 2003. Equine infectious anemia virus utilizes host vesicular protein sorting machinery during particle release. *Journal of Virology* **77**:8440–8447. DOI: https://doi.org/10.1128/jvi.77.15.8440-8447.2003, PMID: 12857913

Thakur M, Buniello A, Brooksbank C, Gurwitz KT, Hall M, Hartley M, Hulcoop DG, Leach AR, Marques D, Martin M, Mithani A, McDonagh EM, Mutasa-Gottgens E, Ochoa D, Perez-Riverol Y, Stephenson J, Varadi M, Velankar S, Vizcaino JA, Witham R, et al. 2024. EMBL's European Bioinformatics Institute (EMBL-EBI) in 2023. *Nucleic Acids Research* **52**:D10–D17. DOI: https://doi.org/10.1093/nar/gkad1088, PMID: 38015445

Valli A, Busnadiego I, Maliogka V, Ferrero D, Castón JR, Rodríguez JF, García JA. 2012. The VP3 factor from viruses of birnaviridae family suppresses rna silencing by binding both long and small rna duplexes. *PLOS ONE* **7**:e45957. DOI: https://doi.org/10.1371/journal.pone.0045957

Walker PJ, Siddell SG, Lefkowitz EJ, Mushegian AR, Adriaenssens EM, Alfenas-Zerbini P, Davison AJ, Dempsey DM, Dutilh BE, García ML, Harrach B, Harrison RL, Hendrickson RC, Junglen S, Knowles NJ, Krupovic M, Kuhn JH, Lambert AJ, Łobocka M, Nibert ML, et al. 2021. Changes to virus taxonomy and to the international code of virus classification and nomenclature ratified by the international committee on taxonomy of viruses (2021). *Archives of Virology* **166**:2633–2648. DOI: https://doi.org/10.1007/s00705-021-05156-1, PMID: 34231026

Wassenaar TA, Ingólfsson HI, Böckmann RA, Tieleman DP, Marrink SJ. 2015. Computational lipidomics with insane: a versatile tool for generating custom membranes for molecular simulations. *Journal of Chemical Theory and Computation* **11**:2144–2155. DOI: https://doi.org/10.1021/acs.jctc.5b00209, PMID: 26574417

Whitmore L, Wallace BA. 2004. DICHROWEB, an online server for protein secondary structure analyses from circular dichroism spectroscopic data. *Nucleic Acids Research* **32**:e371. DOI: https://doi.org/10.1093/nar/gkh371, PMID: 15215473

Whitmore L, Wallace BA. 2008. Protein secondary structure analyses from circular dichroism spectroscopy: methods and reference databases. *Biopolymers* **89**:392–400. DOI: https://doi.org/10.1002/bip.20853, PMID: 17896349

Wickham H. 2016. Ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag. DOI: https://doi.org/10.1007/978-3-319-24277-4
Wills JW, Craven RC. 1991. Form, function, and use of retroviral gag proteins. *AIDS* **5**:639-654. DOI: https://doi.org/10.1097/00002030-199106000-00002, PMID: 1883539

Yanai H, Hayashi Y, Watanabe Y, Ohtaki N, Kobayashi T, Nozaki Y, Ikuta K, Tomonaga K. 2006. Development of a novel Borna disease virus reverse genetics system using RNA polymerase II promoter and SV40 nuclear import signal. *Microbes and Infection* **8**:1522-1529. DOI: https://doi.org/10.1016/j.micinf.2006.01.010, PMID: 16697679

Zhou L, Feng T, Xu S, Gao F, Lam TT, Wang Q, Wu T, Huang H, Zhan L, Li L, Guan Y, Dai Z, Yu G. 2022. ggmsa: a visual exploration tool for multiple sequence alignment and associated data. *Briefings in Bioinformatics* **23**:bbac222. DOI: https://doi.org/10.1093/bib/bbac222, PMID: 35671504

Zidek A, Tomlinson H. 2021. Alphafold. v2.0.0. GitHub. https://github.com/google-deepmind/alphafold