Secondary Structure Motifs Research Articles

Site-Specific Stabilization and Destabilization of α Helical Peptides upon Phosphorylation and O-GlcNAcylation.

Helices (α-helix) are the most common type of secondary structure motif present in proteins. In this study, we have investigated the structural influence of phosphorylation and O-GlcNAcylation, common intracellular post-translational modifications (PTMs), on the α-helical conformation. The simulation studies were performed on the Baldwin model α-helical peptide sequence (Ac-AKAAAAKAAAAKAA-NH2). The Baldwin sequences were chosen due to the availability of site-specific experimental post-translational data for cross-validation with the simulations. The influence of PTMs was examined across the span of the α-helix, namely, at the N-terminus, position 10 (interior region), and the C-terminus for both serine and threonine residues placed at these positions. Molecular dynamics (MD) simulations revealed that phosphorylation and O-GlcNAcylation at the N-terminus lead to the stabilization of the helical conformation. PTMs in the interior or the C-terminus were found to disrupt helicity, with the disruption being more pronounced for PTMs in the interior region, in accordance with experimental studies. It was found that phosphorylation-derived destabilization was mainly due to the formation of an intraresidue HN-PO32- electrostatic interaction and interactions between the phosphate group and the side chain of adjacent lysine residues (NH3···PO32-). Hydrophobic and steric clashes were the main causes of destabilization in the case of O-GlcNAcylation. The structural disruptions were found to be more pronounced for PTM at the threonine site when compared to the serine site. The salt-bridge-dependent stability of the α-helix was found to be highly position specific, an i → i + 4 interaction stabilizing the helix, with other placements leading to the destabilization of the helix.

RNA Secondary Structures Regulate Adsorption of Fragments onto Flat Substrates.

RNA is a functionally rich molecule with multilevel, hierarchical structures whose role in the adsorption to molecular substrates is only beginning to be elucidated. Here, we introduce a multiscale simulation approach that combines a tractable coarse-grained RNA structural model with an interaction potential of a structureless flat adsorbing substrate. Within this approach, we study the specific role of stem-hairpin and multibranch RNA secondary structure motifs on its adsorption phenomenology. Our findings identify a dual regime of adsorption for short RNA fragments with and without the secondary structure and underline the adsorption efficiency in both cases as a function of the surface interaction strength. The observed behavior results from an interplay between the number of contacts formed at the surface and the conformational entropy of the RNA molecule. The adsorption phenomenology of RNA seems to persist also for much longer RNAs as qualitatively observed by comparing the trends of our simulations with a theoretical approach based on an ideal semiflexible polymer chain.

RNALigands: a database and web server for RNA-ligand interactions.

RNA molecules can fold into complex and stable 3D structures, allowing them to carry out important genetic, structural, and regulatory roles inside the cell. These complex structures often contain 3D pockets made up of secondary structural motifs that can be potentially targeted by small molecule ligands. Indeed, many RNA structures in PDB contain bound small molecules, and high-throughput experimental studies have generated a large number of interacting RNA and ligand pairs. There is considerable interest in developing small molecule lead compounds targeting viral RNAs or those RNAs implicated in neurological diseases or cancer. We hypothesize that RNAs that have similar secondary structural motifs may bind to similar small molecule ligands. Toward this goal, we established a database collecting RNA secondary structural motifs and bound small molecule ligands. We further developed a computational pipeline, which takes as input an RNA sequence, predicts its secondary structure, extracts structural motifs, and searches the database for similar secondary structure motifs and interacting small molecule. We demonstrated the utility of the server by querying α-synuclein mRNA 5' UTR sequence and finding potential matches which were validated as correct. The server is publicly available at http://RNALigands.ccbr.utoronto.ca The source code can also be downloaded at https://github.com/SaisaiSun/RNALigands.

Universal and strain specific structure features of segment 8 genomic RNA of influenza A virus—application of 4-thiouridine photocrosslinking

RNA structure in the influenza A virus (IAV) has been the focus of several studies that have shown connections between conserved secondary structure motifs and their biological function in the virus replication cycle. Questions have arisen on how to best recognize and understand the pandemic properties of IAV strains from an RNA perspective, but determination of the RNA secondary structure has been challenging. Herein, we used chemical mapping to determine the secondary structure of segment 8 viral RNA (vRNA) of the pandemic A/California/04/2009 (H1N1) strain of IAV. Additionally, this long, naturally occurring RNA served as a model to evaluate RNA mapping with 4-thiouridine (4sU) crosslinking. We explored 4-thiouridine as a probe of nucleotides in close proximity, through its incorporation into newly transcribed RNA and subsequent photoactivation. RNA secondary structural features both universal to type A strains and unique to the A/California/04/2009 (H1N1) strain were recognized. 4sU mapping confirmed and facilitated RNA structure prediction, according to several rules: 4sU photocross-linking forms efficiently in the double-stranded region of RNA with some flexibility, in the ends of helices, and across bulges and loops when their structural mobility is permitted. This method highlighted three-dimensional properties of segment 8 vRNA secondary structure motifs and allowed to propose several long-range three-dimensional interactions. 4sU mapping combined with chemical mapping and bioinformatic analysis could be used to enhance the RNA structure determination as well as recognition of target regions for antisense strategies or viral RNA detection.

Mechanistic investigations in sonoenzymatic synthesis of n-butyl levulinate

n-Butyl levulinate (n-BL) has numerous applications in biofuel and green solvents. In this paper, we have reported mechanistic accounts of sonoenzymatic synthesis of n-BL using Novozym 435. Esterification parameters were optimized with statistical design of experiments. Application of 35 kHz sonication boosted n-BL yield from 70.9% to 92.22%. Mechanistic investigation using Ping-Pong Bi-Bi kinetics model and molecular docking simulations revealed interesting influence of sonication on esterification reaction. Exposure to sonication induced significant changes in the secondary structure of enzyme, as revealed in ATR-FTIR analysis. Sonication increased α-helix content of enzyme that increased enzyme activity by opening of flapping lid and widening the catalytic cavity. Sonication also caused unfolding of secondary structural motifs with rise in random coil content that favored formation of enzyme-ligand complexes. These effects enhanced reaction velocity and substrate affinity with reduction in inhibition and unfavorable dissociation of intermediate complexes, which ultimately enhanced n-BL yield.

E2EDNA: Simulation Protocol for DNA Aptamers with Ligands.

We present E2EDNA, a simulation protocol and accompanying code for the molecular biophysics and materials science communities. This protocol is both easy to use and sufficiently efficient to simulate single-stranded (ss)DNA and small analyte systems that are important to cellular processes and nanotechnologies such as DNA aptamer-based sensors. Existing computational tools used for aptamer design focus on cost-effective secondary structure prediction and motif analysis in the large data sets produced by SELEX experiments. As a rule, they do not offer flexibility with respect to the choice of the theoretical engine or direct access to the simulation platform. Practical aptamer optimization often requires higher accuracy predictions for only a small subset of sequences suggested, e.g., by SELEX experiments, but in the absence of a streamlined procedure, this task is extremely time and expertise intensive. We address this gap by introducing E2EDNA, a computational framework that accepts a DNA sequence in the FASTA format and the structures of the desired ligands and performs approximate folding followed by a refining step, analyte complexation, and molecular dynamics sampling at the desired level of accuracy. As a case study, we simulate a DNA-UTP (uridine triphosphate) complex in water using the state-of-the-art AMOEBA polarizable force field. The code is available at https://github.com/InfluenceFunctional/E2EDNA.

Insights into Structure and Aggregation Behavior of Elastin-like Polypeptide Coacervates: All-Atom Molecular Dynamics Simulations.

The stimuli-responsive character of elastin-like polypeptides (ELP) has led to their use in a wide range of applications. The temperature-triggered aggregation, or LCST behavior, of ELPs is a complex and multistep phenomenon, which proposed to include the structural transitions, loss of hydrophobic hydration, expulsion of water molecules and physical association of chains. Thus, the origin and detailed mechanism of LCST in ELPs is difficult to elucidate. Here, to gain insights into structure and dynamics of coacervates, we performed all-atom molecular dynamics simulations of 27 90-mer ELPs in explicit water at 350 K. Two sequences, poly(VGPVG)18 and poly(VPGVG)18, were examined due to their experimentally observed differences in thermal hysteresis albeit identical overall composition but different arrangement of amino acids. The simulation results indicate that surface hydrophobicity of poly(VGPVG) aggregate is less than that of the poly(VPGVG) aggregate, and there are marked changes in torsion angles and the propensities of secondary structural motifs during the aggregation process. Moreover, there are significant differences between structure of a single polypeptide in water and structure within the aggregate. Overall, the aggregation process is driven by the formation of peptide-peptide interactions whereas the average hydration of peptides remains almost the same between dissolved and aggregated states. Even though the aggregation is driven by the hydrophobic interactions, ELP coacervate has no hydrophobic core and contains many water molecules. Overall, our findings provide an insight into the sequence-dependent structure of coacervates and molecular behavior of individual peptides during aggregation.

Protein-Protein Interface Topology as a Predictor of Secondary Structure and Molecular Function Using Convolutional Deep Learning.

To power the specific recognition and binding of protein partners into functional complexes, a wealth of information about the structure and function of the partners is necessarily encoded into the global shape of protein-protein interfaces and their local topological features. To identify whether this is the case, this study uses convolutional deep learning methods (typically leveraged for 2D image recognition) on 3D voxel representations of protein-protein interfaces colored by burial depth. A novel two-stage network fed with voxelizations of each interface at two distinct resolutions achieves balance between performance and computational cost. From the shape of the interfaces, the network tries to predict the presence of secondary structure motifs at the interface and the molecular function of the corresponding complex. Secondary structure and certain classes of function are found to be very well predicted, validating the hypothesis that interface shape is a conveyor of higher-level information. Interface patterns triggering the recognition of specific classes are also identified and described.

BRIO: a web server for RNA sequence and structure motif scan.

The interaction between RNA and RNA-binding proteins (RBPs) has a key role in the regulation of gene expression, in RNA stability, and in many other biological processes. RBPs accomplish these functions by binding target RNA molecules through specific sequence and structure motifs. The identification of these binding motifs is therefore fundamental to improve our knowledge of the cellular processes and how they are regulated. Here, we present BRIO (BEAM RNA Interaction mOtifs), a new web server designed for the identification of sequence and structure RNA-binding motifs in one or more RNA molecules of interest. BRIO enables the user to scan over 2508 sequence motifs and 2296 secondary structure motifs identified in Homo sapiens and Mus musculus, in three different types of experiments (PAR-CLIP, eCLIP, HITS). The motifs are associated with the binding of 186 RBPs and 69 protein domains. The web server is freely available at http://brio.bio.uniroma2.it.

Conserved Structural Motifs of Two Distant IAV Subtypes in Genomic Segment 5 RNA.

The functionality of RNA is fully dependent on its structure. For the influenza A virus (IAV), there are confirmed structural motifs mediating processes which are important for the viral replication cycle, including genome assembly and viral packaging. Although the RNA of strains originating from distant IAV subtypes might fold differently, some structural motifs are conserved, and thus, are functionally important. Nowadays, NGS-based structure modeling is a source of new in vivo data helping to understand RNA biology. However, for accurate modeling of in vivo RNA structures, these high-throughput methods should be supported with other analyses facilitating data interpretation. In vitro RNA structural models complement such approaches and offer RNA structures based on experimental data obtained in a simplified environment, which are needed for proper optimization and analysis. Herein, we present the secondary structure of the influenza A virus segment 5 vRNA of A/California/04/2009 (H1N1) strain, based on experimental data from DMS chemical mapping and SHAPE using NMIA, supported by base-pairing probability calculations and bioinformatic analyses. A comparison of the available vRNA5 structures among distant IAV strains revealed that a number of motifs present in the A/California/04/2009 (H1N1) vRNA5 model are highly conserved despite sequence differences, located within previously identified packaging signals, and the formation of which in in virio conditions has been confirmed. These results support functional roles of the RNA secondary structure motifs, which may serve as candidates for universal RNA-targeting inhibitory methods.

Searching for Aquamelt Behavior among Silklike Biomimetics during Fibrillation under Flow.

In this paper, we elucidate a generic mechanism behind strain-induced phase transition in aqueous solutions of silk-inspired biomimetics by atomistic molecular dynamics simulations. We show the results of modeling of homopeptides polyglycine Gly30 and polyalanine Ala30 and a heteropeptide (Gly-Ala-Gly-Ala-Gly-Ser)5, i.e., the simplest and yet relevant sequences that could mimic the behavior of natural silk under stress conditions. First, we analyze hydrophobicities of the sequences by calculating the Gibbs free energy of hydration and inspecting the interchain hydrogen bonding and hydration by water. Second, the force-extension profiles are scanned and compared with the results for poly(ethylene oxide), the synthetic polymer for which the aquamelt behavior has been proved recently. Additionally, the conformational transitions of oligopeptides from coiled to extended states are characterized by a generalized order parameter and by the dependence of the solvent-accessible surface area of the chains on applied stretching. Fibrillation itself is surveyed using both the two-dimensional interchain pair correlation function and the SAXS/WAXS patterns for the aggregates formed under stress. These are compared with experimental data found in the literature on fibril structure of silk composite materials doped with oligoalanine peptides. Our results show that tensile stress introduced into aqueous oligopeptide solutions facilitates interchain interactions. The oligopeptides display both a greater resistance to extension as compared to poly(ethylene oxide) and a reduced ability for hydrogen bonding of the stretched chains between oligomers and with water. Fiber formation is proved for all simulated objects, but the most structured one is made of a heteropeptide (Gly-Ala-Gly-Ala-Gly-Ser)5: For this sequence, we obtain the highest degree of the secondary structure motifs in the fiber. We conclude that this is the most promising candidate among considered sequences to find the aquamelt behavior in further experimental studies.

DeepCys: Structure-based multiple cysteine function prediction method trained on deep neural network: Case study on domains of unknown functions belonging to COX2 domains.

Cysteine (Cys) is the most reactive amino acid participating in a wide range of biological functions. In-silico predictions complement the experiments to meet the need of functional characterization. Multiple Cys function prediction algorithm is scarce, in contrast to specific function prediction algorithms. Here we present a deep neural network-based multiple Cys function prediction, available on web-server (DeepCys) (https://deepcys.herokuapp.com/). DeepCys model was trained and tested on two independent datasets curated from protein crystal structures. This prediction method requires three inputs, namely, PDB identifier (ID), chain ID and residue ID for a given Cys and outputs the probabilities of four cysteine functions, namely, disulphide, metal-binding, thioether and sulphenylation and predicts the most probable Cys function. The algorithm exploits the local and global protein properties, like, sequence and secondary structure motifs, buried fractions, microenvironments and protein/enzyme class. DeepCys outperformed most of the multiple and specific Cys function algorithms. This method can predict maximum number of cysteine functions. Moreover, for the first time, explicitly predicts thioether function. This tool was used to elucidate the cysteine functions on domains of unknown functions belonging to cytochrome C oxidase subunit-II like transmembrane domains. Apart from the web-server, a standalone program is also available on GitHub (https://github.com/vam-sin/deepcys).

Structural and biophysical correlation of anti-NANP antibodies with in vivo protection against P. falciparum

The most advanced P. falciparum circumsporozoite protein-based malaria vaccine, RTS,S/AS01 (RTS,S), confers partial protection but with antibody titers that wane relatively rapidly, highlighting the need to elicit more potent and durable antibody responses. Here, we elucidate crystal structures, binding affinities and kinetics, and in vivo protection of eight anti-NANP antibodies derived from an RTS,S phase 2a trial and encoded by three different heavy-chain germline genes. The structures reinforce the importance of homotypic Fab-Fab interactions in protective antibodies and the overwhelmingly dominant preference for a germline-encoded aromatic residue for recognition of the NANP motif. In this study, antibody apparent affinity correlates best with protection in an in vivo mouse model, with the more potent antibodies also recognizing epitopes with repeating secondary structural motifs of type I β- and Asn pseudo 310 turns; such insights can be incorporated into design of more effective immunogens and antibodies for passive immunization.

Connecting Structure with Kinetics of Single-Stranded DNA Fluctuations that Provide Potential Binding Sites for Replication Proteins

The local conformations of DNA at and near replication forks and primer-template (p/t) junctions likely regulate the non-sequence-specific binding of replication proteins. Thermal fluctuations can drive the exploration of local DNA conformations, thereby providing opportunities for proteins to stably bind to transiently formed secondary structure motifs. In this work, we apply microsecond-resolved single-molecule Förster resonance energy transfer (sm FRET) experiments to study the equilibrium distribution and conformational transitions of the end-to-end distances of an oligo-deoxythymidine chain (dT)15 that serves as the single-stranded (ss) region of a Cy3-Cy5 labeled p/t DNA construct. By comparing our experimental data to the results of a kinetic network model, we determined that the system exists as three distinct macrostates with different end-to-end extensions, and that transitions between these macrostates occur over time scales ranging from ∼ 10 microseconds to ∼ 100 milliseconds. In addition, we measured the relative stabilities and transition barrier heights that characterize the free energy landscape of the (dT)15 p/t DNA system. The magnitudes of the activation barriers between macrostates of varying chain extension likely depend on the rates of rearrangement of the various populations of base stacked segments within the single-stranded portion of the p/t DNA construct. In ongoing work, we are investigating the physical nature of the transition barriers between macrostates in terms of local base-stacking interactions, using a combination of Monte Carlo simulation and spectroscopy of fluorescent base analogue substituted DNA constructs. Certain conformations within these macrostates are likely to represent targets for the binding of proteins, such as the single-stranded DNA binding protein (gp32) of bacteriophage T4, to ssDNA templates in the assembly and function of the replication complex.

Target-binding behavior of IDPs via pre-structured motifs.

Pre-Structured Motifs (PreSMos) are transient secondary structures observed in many intrinsically disordered proteins (IDPs) and serve as protein target-binding hot spots. The prefix "pre" highlights that PreSMos exist a priori in the target-unbound state of IDPs as the active pockets of globular proteins pre-exist before target binding. Therefore, a PreSMo is an "active site" of an IDP; it is not a spatial pocket, but rather a secondary structural motif. The classical and perhaps the most effective approach to understand the function of a protein has been to determine and investigate its structure. Ironically or by definition IDPs do not possess structure (here structure refers to tertiary structure only). Are IDPs then entirely structureless? The PreSMos provide us with an atomic-resolution answer to this question. For target binding, IDPs do not rely on the spatial pockets afforded by tertiary or higher structures. Instead, they utilize the PreSMos possessing particular conformations that highly presage the target-bound conformations. PreSMos are recognized or captured by targets via conformational selection (CS) before their conformations eventually become stabilized via structural induction into more ordered bound structures. Using PreSMos, a number of, if not all, IDPs can bind targets following a sequential pathway of CS followed by an induced fit (IF). This chapter presents several important PreSMos implicated in cancers, neurodegenerative diseases, and other diseases along with discussions on their conformational details that mediate target binding, a structural rationale for unstructured proteins.

Development of fructose-1,6-bisphosphate aldolase enzyme peptide mimics as biocatalysts in direct asymmetric aldol reactions.

This study describes the design and synthesis of mimetic peptides modelled on the catalytic active site of the fructose-1,6-bisphosphate aldolase (FBPA) enzyme. The synthesized peptides consisting of the turn motifs and catalytic site amino acids of FBPA enzyme were evaluated for catalytic activity in direct asymmetric aldol reactions of ketones and aldehydes. The influence of substrate scope, catalyst loading and solvents including water, on the reaction were also investigated. Nuclear magnetic resonance (NMR) and circular dichroism (CD) were used to determine the secondary structure of the peptides to provide an understanding of the structure–activity relationship. The peptides showed catalytic activity and the aldol products were obtained in low yields (up to 44%), but excellent enantioselectivity (up to 93%) and moderate diastereoselectivity (65 : 35).

Polymorphism of methyl 4-amino-3-phenylisothiazole-5-carboxylate: an experimental and theoretical study

Being a close analogue of amflutizole, methyl 4-amino-3-phenylisothiazole-5-carboxylate (C11H10N2O2S) was assumed to be capable of forming polymorphic structures. Noncentrosymmetric and centrosymmetric polymorphs have been obtained by crystallization from a series of more volatile solvents and from denser tetrachloromethane, respectively. Identical conformations of the molecule are found in both structures. The two polymorphs differ mainly in the intermolecular interactions formed by the amino group and in the type of stacking interactions between the π-systems. The most effective method for revealing packing motifs in structures with intermolecular interactions of different types (hydrogen bonding, stacking, dispersion, etc.) is to study the pairwise interaction energies using quantum chemical calculations. Molecules form a column as the primary basic structural motif due to stacking interactions in both polymorphic structures under study. The character of a column (straight or zigzag) is determined by the orientations of the stacked molecules (in a `head-to-head' or `head-to-tail' manner). Columns bound by intermolecular N-H...O and N-H...N hydrogen bonds form a double column as the main structural motif in the noncentrosymmetric structure. Double columns in the noncentrosymmetric structure and columns in the centrosymmetric structure interact strongly within the ab crystallographic plane, forming a layer as a secondary basic structural motif. The noncentrosymmetric structure has a lower density and a lower (by 0.59 kJ mol-1) lattice energy, calculated using periodic calculations, compared to the centrosymmetric structure.

Protein secondary structure motifs: A kinematic construction.

The kinematic geometry of protein backbone structures, constrained by either single or multiple hydrogen bonds (H-bonds), possibly in a periodic array, is discussed. These structures include regular secondary structure elements α-helices and β-sheets but also include other short H-bond stabilized irregular structural elements like β-turns. The work here shows that the variations observed in such structures have simple geometrical correlations consistent with constrained motion kinematics. A new classification of the ideal helices is given, in terms of the parameter α, the angle at a Cα atom to its two neighboring Cα 's along the helix, and shown how it can be generalized to include nonideal helices. Specifically, we derive an analytical expression of the backbone dihedrals, (ϕ, ψ), in terms of the parameter α subject to the constraint that the peptide planes are parallel to the helical axis. Helices constructed in this way exhibit near-vertical alignment of the C = O and N - H units and are the canonical objects of this study. These expressions are easily modifiable to include perturbations of parameters relevant to nonplanar peptide units and noncanonical angles. The addition of a second parameter, ε0 , inclination of successive peptide planes along a helix with respect to the helical axis leads to a generalization of the previous expression and provides an efficient parametrization of such structures in terms of coordinates consistent with H-bond parameters. An analogs parametrization of β-turns, using inverse kinematic methods, is also given. Besides offering a unifying viewpoint, our results may find useful applications to protein and peptide design.

Incorporating a Thiophosphate Modification into a Common RNA Tetraloop Motif Causes an Unanticipated Stability Boost

GNRA (N = A, C, G, or U; R = A or G) tetraloops are common RNA secondary structural motifs and feature a phosphate stacked atop a nucleobase. The rRNA sarcin/ricin loop (SRL) is capped by GApGA, and the phosphate p stacks on G. We recently found that regiospecific incorporation of a single dithiophosphate (PS2) but not a monothiophosphate (PSO) instead of phosphate in the backbone of RNA aptamers dramatically increases the binding affinity for their targets. In the RNA:thrombin complex, the key contribution to the 1000-fold tighter binding stems from an edge-on contact between PS2 and a phenylalanine ring. Here we investigated the consequences of replacing the SRL phosphate engaged in a face-on interaction with guanine with either PS2 or PSO for stability. We found that PS2···G and Rp-PSO···G contacts stabilize modified SRLs compared to the parent loop to unexpected levels: up to 6.3 °C in melting temperature Tm and -4.7 kcal/mol in ΔΔG°. Crystal structures demonstrate that the vertical distance to guanine for the closest sulfur is just 0.05 Å longer on average compared to that of oxygen despite the larger van der Waals radius of the former (1.80 Å for S vs 1.52 Å for O). The higher stability is enthalpy-based, and the negative charge as assessed by a neutral methylphosphonate modification plays only a minor role. Quantum mechanical/molecular mechanical calculations are supportive of favorable dispersion attraction interactions by sulfur making the dominant contribution. A stacking interaction between phosphate and guanine (SRL) or uracil (U-turn) is also found in newly classified RNA tetraloop families besides GNRA.

Organization of the Influenza A Virus Genomic RNA in the Viral Replication Cycle-Structure, Interactions, and Implications for the Emergence of New Strains.

The influenza A virus is a human pathogen causing respiratory infections. The ability of this virus to trigger seasonal epidemics and sporadic pandemics is a result of its high genetic variability, leading to the ineffectiveness of vaccinations and current therapies. The source of this variability is the accumulation of mutations in viral genes and reassortment enabled by its segmented genome. The latter process can induce major changes and the production of new strains with pandemic potential. However, not all genetic combinations are tolerated and lead to the assembly of complete infectious virions. Reports have shown that viral RNA segments co-segregate in particular circumstances. This tendency is a consequence of the complex and selective genome packaging process, which takes place in the final stages of the viral replication cycle. It has been shown that genome packaging is governed by RNA–RNA interactions. Intersegment contacts create a network, characterized by the presence of common and strain-specific interaction sites. Recent studies have revealed certain RNA regions, and conserved secondary structure motifs within them, which may play functional roles in virion assembly. Growing knowledge on RNA structure and interactions facilitates our understanding of the appearance of new genome variants, and may allow for the prediction of potential reassortment outcomes and the emergence of new strains in the future.