Background {#Sec1}
==========

*Tospovirus* is a thrips-borne plant-infecting genus within the family *Bunyaviridae*, which also includes four animal-infecting genera: *Hanta/Nairo/Phlebo-* and *Orthobunyavirus* \[[@CR1]\]. GRSV (*Groundnut ringspot virus)* is an emerging tospovirus, that has caused severe diseases in distinct vegetable crops in South America and is phylogenetically close to the tospovirus type-species TSWV (*Tomato spotted wilt virus*) \[[@CR2]\]. Like all tospoviruses, GRSV contain a trisegmented negative single-stranded RNA (ssRNA) genome that encodes the viral RNA-dependent RNA polymerase (RdRp), two glycoproteins (Gn/Gc), the movement protein (NSm), the RNA silencing suppressor protein (NSs) and the nucleoprotein (N) \[[@CR3]\]. N is a multifunctional protein involved in RNA protection, particle assembly, intracellular movement and might play a role in transcription/replication regulation \[[@CR4]--[@CR14]\]. Multiple copies of the N protein form oligomers that interact with the viral RNAs to build ribonucleoprotein complexes (RNPs) that are proposed to be transported via plasmodesmata and are functional templates for RNA replication and transcription \[[@CR6], [@CR15], [@CR16]\].

Multiple studies have attempted to identify candidate N protein regions involved in RNA binding and protein multimerization for TSWV using yeast two-hybrid systems (Y2HS) and site-directed mutagenesis \[[@CR4], [@CR6], [@CR17], [@CR18]\], but the tospovirus RNPs remains largely uncharacterized at the molecular level and the lack of structural information prevents detailed insight into these interactions. The lack of a reverse genetics system, which is available for other bunyaviruses, has hampered tospovirus research. The N protein crystal structures of related RNA virus families (*Arena*/*Orthomyxo*/*Bunyaviridae*) have been elucidated \[[@CR8], [@CR19]--[@CR26]\] and despite different size and distinct N-folding structures, there are common features and architectural principles by which these proteins form N-N multimers and N-RNA complexes \[[@CR27]\]. Therefore, these available structures were used to predict a three-dimensional model for GRSV N (the most important and prevalent tospovirus in Brazil) using homology modeling.

Results and discussion {#Sec2}
======================

Three-dimensional model of GRSV N and oligomerization {#Sec3}
-----------------------------------------------------

The GRSV N and LACV N have similar protein fold with the predicted GRSV N monomer forming thirteen helical segments and two small beta-sheets (Figs. [1](#Fig1){ref-type="fig"}, 2a and e-f). The protein has a globular core domain (26--223 aa) containing a deep positively charged groove with the two chain terminals forming an N-terminus arm (1--25 aa) and a C-terminus arm (224--258 aa) (Fig. [2a-b](#Fig2){ref-type="fig"} and Fig. [3](#Fig3){ref-type="fig"}). The N- and C-arms extend outwards from the globular core domain and interacts with the globular core domain of neighboring monomers to mediate the multimerization, supporting the "head-to-tail" model proposed by \[[@CR18]\]. Amino acids S2-V12 of the N-arm interact with the Q61-N82 of the core domain of one neighboring monomer (Fig. [2c and e](#Fig2){ref-type="fig"}) while K227-K249 of the C-arm interact with the K173--K198 of the core domain the other neighboring monomer (Fig. [2d and f](#Fig2){ref-type="fig"}). Specific residue-residue interactions have been listed in Table [1](#Tab1){ref-type="table"} for the two independent interfaces. According to PISA, the intermolecular interactions were mainly hydrogen bonds, but van der Waals and hydrophobic interactions also contribute to hold the monomers together (data not shown). This interaction model is further corroborated by the available mutational studies on TSWV N \[[@CR4], [@CR17], [@CR18]\].Fig. 1*Groundnut ringspot virus* (GRSV) and *La crosse virus* (LACV) Nucleoproteins sequence alignment. Key residues for GRSV N and LACV N oligomerization and for ssRNA binding are colored as indicated by the colored bars. The secondary structure of LACV N is shown above, and every 10 residues are indicated with a dot (.). Strictly conserved residues are highlighted in red with white letter and highly conserved residues are displayed by red letters. GRSV N GenBank accession number is AF251271 and LACV N UniProt accession code is P04873 Fig. 2Monomeric and tetrameric structure of the *Groundnut ringspot virus* (GRSV) nucleoprotein (N). **a** Cartoon representation of monomeric GRSV N with rainbow coloring from N- (blue) to C-terminus (red). **b** Electrostatic surface of the GRSV N with a positively charged groove in complex with RNA shown as yellow (carbons) and red (oxygens) sticks. Positive and negative charges are blue and red, respectively. **c** N-terminus interaction surface representation of four GRSV N monomers A, B, C, D shown in color pink, yellow, cyan and green, respectively. **d** C-terminus interaction surface representation of the GRSV N tetrameric ring. The RNA is shown in black sticks deeply bound inside the tetrameric ring. **e** Cartoon representation with the N-arm oligomerization interface showing interacting residues. The N-terminal arm is in pink and the globular region is in green. The intermolecular hydrogen bonds are shown as yellow dotted lines. **f** 180° rotation of Fig. 2e, C-arm oligomerization interface showing interacting residues. The C-terminal arm is in cyan and the globular region is in green Fig. 3Sequence alignment of representative tospoviruses Nucleoproteins (N). The secondary structure of *Groundnut ringspot virus* (GRSV) is shown above and of *La crosse virus* (LACV) is shown at the bottom. Key residues for GRSV N ssRNA binding are marked with yellow triangles. GRSV N- and C-arms are marked with blue and green boxes respectively, with key residues for oligomerization highlighted. Strictly conserved residues are highlighted in red with white letter and highly conserved residues with red letter. I: *Tospovirus* American clade I; II: *Tospovirus* American clade II; III: *Tospovirus* Eurasian clade; IV: *Orthobunyavirus*. The sequence codes are supplied at the Additional file [2](#MOESM2){ref-type="media"}: Table S1 Table 1Pairs of interacting residues for GRSV N-N oligomerizationN-arm^a^N-arm binding site^b^C-arm^c^C-arm binding site^b^S2^d^S83A226K183V4N82S229K183T7Q61D233T186T7S62Y235K183K8S62N238K198N10T73Y243N185V12T73V246K175V12G75V248K173K249Y174^a^N-arm amino acids residues of GRSV N^b^Interacting amino acids residues of GRSV N globular core domain^c^C-arm amino acids residues of GRSV N^d^Amino acids residues position in the GRSV N sequence

Actually, the first assay to map functional domains of TSWV N, performing Y2HS and random serial deletions, showed that both the N- (1--39 aa) and C-terminals (233--248 aa) were important for N-N interaction \[[@CR18]\], in clear agreement with the structural results presented here. Furthermore, \[[@CR17]\] identified three crucial intermonomer binding regions: 42--56, 132--152 and 222--248 which have a clear correspondence with the predicted interaction residues of GRSV N located at N- and C-arms, or buried in the core of the model (Fig. [3](#Fig3){ref-type="fig"}). Moreover, amino acids residues located at the regions K103-A119 and L132-V135 are solvent accessible and therefore are able to interact with NSm, glycoproteins, viral polymerase or host proteins \[[@CR6], [@CR7]\]. Recently, studies have been performed attempting to identify N-NSm interactions \[[@CR28], [@CR29]\] which results are in perfect congruence with the GRSV N protein model. In both cases, the model proposed here represents an efficient tool to assist in planning experiments with mutations and deletion in the N protein.

In addition, the obtained model for N protein was submitted to molecular dynamics simulations in order to both refine the structure in aqueous solvent \[[@CR30], [@CR31]\] and access the protein conformational ensemble, further exploring its structural and functional roles. During the simulation time, the globular core domain did not reveal any loss of secondary structure, increase of radius of gyration or persistent increments on RMSD values, which supports the model quality. It is worthy to mention that RMSF calculations indicate the N-terminal arm (1--25 aa) as a very flexible region (Fig. [4c](#Fig4){ref-type="fig"}).Fig. 4Molecular Dynamics of monomeric Nucleoprotein (N) of *Groundnut ringspot virus* (GRSV). Root Mean Square Deviation (RMSD) calculations for different set of atoms in both presence **a** and absence **b** of RNA. **c** Root Mean Square Fluctuations (RMSF) calculations for the entire N protein in both presence (*red*) and absence (*black*) of RNA. **d** Plot of α-helix content as function of time in both presence (*red*) and absence (*black*) of RNA

RNA interaction {#Sec4}
---------------

According to the GRSV N protein model, the RNA is primarily bound at the central RNA-binding groove (Fig. [2b](#Fig2){ref-type="fig"}), and the key residues for this interaction (K3, K5, Q17, K58, R60, Q61, R94, R95, K183, Y184, K187, K192 and K227) are mainly located in this positively charged groove. This positively charged groove is only possible because residues F37, F56, F72, F74, I79, M91, F93 and L96 form a hydrophobic core, which is indispensable to stabilize the protein folding and to correctly orient the RNA interacting residues towards the groove. Importantly, these residues are highly conserved among all tospoviruses (Fig. [3](#Fig3){ref-type="fig"}). Note that the N-terminal arm is also involved in RNA binding and shielding RNA from the solvent (Fig. [2c-d](#Fig2){ref-type="fig"}). Residues F23, L54, F56, L57 and F93 were observed to modulate the RNA nucleobases dynamics during the performed simulation, while the N-terminal arm seems to play a stabilization role during MD simulations of GRSV N protein (Fig. [4a and b](#Fig4){ref-type="fig"}). In addition, the content of alpha-helices in GRSV N protein bound to RNA increased 25 % during the simulation in comparison to the free monomer (Fig. [4d](#Fig4){ref-type="fig"}), suggesting that, in the simulated timescale, the monomeric state does not present a lack of conformational stability in detriment of oligomeric states, as observed experimentally for other viruses \[[@CR32], [@CR33]\].

Recently, the residues R60, R94, and R95 were confirmed to interact with RNA \[[@CR33]\], which also supports our results. RNA is strongly bent at each N-N interface and is largely solvent-inaccessible in the tetramer (Fig. [2d](#Fig2){ref-type="fig"}). The dimensions of the groove can accommodate ssRNA and PISA analysis showed that the majority of residue-nucleotide interactions occur with the ribose and the phosphate moieties, suggesting a non-sequence-specific RNA interaction. Indeed, Richmond et al. \[[@CR4]\] carried out mutagenesis and gel shift assay studies to identify N regions important for ssRNA binding and demonstrated that the N-RNA complex is highly stable and non-sequence-specific, further supporting these results.

Conclusions {#Sec5}
===========

Taken together, these data will help designing further and more accurate mutational and functional studies of tospovirus N proteins. In addition, the proposed model may shed light on the mechanisms of RNP shaping and could allow the identification of essential amino acid residues as potential targets for tospovirus control strategies.

Methods {#Sec6}
=======

*In silico* homology modeling and model optimization {#Sec7}
----------------------------------------------------

A template for modeling the GRSV N protein was searched in expasy SWISS-MODEL server \[[@CR34]\] using the amino acid sequence of GRSV N as a reference. Template crystal structures of *Orthobunyavirus* genus were chosen due to their genetic relationship. The LACV (*La Crosse virus-Orthobunyavirus*) N tetrameric crystal structure in complex with ssRNA (PDB ID 4BHH) was selected as the template \[[@CR20]\], aligned with GRSV N using T-Coffee server \[[@CR35]\] and the resulting alignment was manually improved using BioEdit \[[@CR36]\]. Aligned sequences were used with MODELLERv9.10 \[[@CR37]\] to develop high quality tetrameric models along with or without RNA.

Optimization of the models was achieved using energy minimization protocols available at Yasara \[[@CR38]\] and Chiron \[[@CR39]\] servers. Quality of the 3D models were evaluated with ERRAT (version 2.0) \[[@CR40]\] and MOL probity \[[@CR41]\]. Ramachandran plots for the models were assessed and Ramachandran outlier residues were fixed with COOT \[[@CR42]\] and energy minimization. The highest quality model with 90.1 % residues in favored region and 8.4 % in allowed region while 1.5 % outlier at Ramachandran plot was selected after visual inspection (see Additional file [1](#MOESM1){ref-type="media"}: Figure S1). The model was subjected to the PISA program \[[@CR43]\] for interface analysis at EBI-EMBL server and the retrieved PISA data was analyzed for binding patterns using PyMOL \[[@CR44]\].

Molecular dynamics {#Sec8}
------------------

Molecular dynamics techniques were applied using GROMACS suite \[[@CR45]\] in order to evaluate the stability and consistency of the obtained N protein monomeric model and investigate GRSV N protein-RNA interactions over time. Therefore, N protein model was simulated in the presence and absence of the modeled RNA, in two analytical systems. Amber99SB-ILDN force field \[[@CR46]\] was used to generate proper topologies. The models were placed at the center of a dodecahedral box and solvated with TIP3P water model \[[@CR47]\]. Counterions were used to neutralize the net charge of the system, and 0.15 M of NaCl was added to the box in order to simulate cellular ionic environment.

After a minimization protocol using steepest descent and conjugate gradient to eliminate possible clashes and bad contacts, NVT ensemble with restraint forces of 1000 kJ/mol was carried for 4 ns at 300 K. Moreover, five subsequent equilibration steps in NPT ensemble were carried out at 1 bar with restraint forces of 800 kJ/mol on heavy atoms, 600 Kcal/(mol x nm) and 400 kJ/mol on mainchain, 200 kJ/mol on backbone and 100 kJ/mol on alpha-carbons, totalizing 13 ns. Finally, production runs with no restraints were carried for 50 ns using an integration step of 2 fs and LINCS algorithm \[[@CR48]\]. Also, Particle Mesh Ewald method \[[@CR49]\] was applied for Coulombic and Lennard-Jones interactions longer than 1 nm.

Additional files {#Sec9}
================

Additional file 1:Ramachandran plot analysis of predicted structure of *Groundnut ringspot virus* (GRSV) N protein. The regions covered by light blue lines show most favored regions, while the regions covered by dark blue lines show allowed regions. Other regions of the plot show the disallowed region. The pink dots show the outliers (PNG 87 kb) Additional file 2:The Genbank acession numbers of the viruses used at this work (TABLEDOCX 19 kb)

Declarations {#FPar1}
============

This article has been published as part of BMC Bioinformatics Volume 17 Supplement 18, 2016. Proceedings of X-meeting 2015: 11th International Conference of the AB3C + Brazilian Symposium on Bioinformatics: bioinformatics. The full contents of the supplement are available online <https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume-17-supplement-18>.

Funding {#FPar2}
=======

Publication of this paper has been funded by CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico), CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior), and FAPDF (Fundação de Apoio à Pesquisa do Distrito Federal), Brazil.

Availability of data and materials {#FPar3}
==================================

The datasets supporting the conclusions of this article are included within the article and its supplementary files.

Authors' contributions {#FPar4}
======================

Conceived and designed the experiments: JARGB MF RNL. Performed the homology modeling: MF RNL. Performed the Molecular Dynamics: HV MDP. Analyzed the data: FLM HV JARGB MDP MF RNL ROR. Wrote the paper: FLM HV JARGB MDP MF RNL ROR. All authors read and approved the final manuscript.

Competing interests {#FPar5}
===================

The authors declare that they have no competing interests.

Consent for publication {#FPar6}
=======================

Not applicable.

Ethics approval and consent to participate {#FPar7}
==========================================

Not applicable.
