Significance
Medical visits and missed days of school and work caused by rhinoviruses cost tens of billions of US dollars annually. Currently, there are no antivirals against rhinoviruses, and the available treatments only treat the symptoms. Here, we present the molecular structure of human rhinovirus 14 in complex with its cellular receptor intercellular adhesion molecule 1. The binding of the virus to its receptor initiates the infection. Knowledge of the structure of the human rhinovirus 14intercellular adhesion molecule 1 interface and mechanism of interaction provides the basis for the design of compounds that may block the binding of rhinoviruses to receptors and thus prevent infection.
Most rhinoviruses, which are the leading cause of the common cold, utilize intercellular adhesion molecule-1 (ICAM-1) as a receptor to infect cells. To release their genomes, rhinoviruses convert to activated particles that contain pores in the capsid, lack minor capsid protein VP4, and have an altered genome organization. The binding of rhinoviruses to ICAM-1 promotes virus activation; however, the molecular details of the process remain unknown. Here, we present the structures of virion of rhinovirus 14 and its complex with ICAM-1 determined to resolutions of 2.6 and 2.4 , respectively. The cryo-electron microscopy reconstruction of rhinovirus 14 virions contains the resolved density of octanucleotide segments from the RNA genome that interact with VP2 subunits. We show that the binding of ICAM-1 to rhinovirus 14 is required to prime the virus for activation and genome release at acidic pH. Formation of the rhinovirus 14ICAM-1 complex induces conformational changes to the rhinovirus 14 capsid, including translocation of the C termini of VP4 subunits, which become poised for release through pores that open in the capsids of activated particles. VP4 subunits with altered conformation block the RNAVP2 interactions and expose patches of positively charged residues. The conformational changes to the capsid induce the redistribution of the virus genome by altering the capsidRNA interactions. The restructuring of the rhinovirus 14 capsid and genome prepares the virions for conversion to activated particles. The high-resolution structure of rhinovirus 14 in complex with ICAM-1 explains how the binding of uncoating receptors enables enterovirus genome release.
Human rhinoviruses are the cause of more than half of common colds (1). Medical visits and missed days of school and work cost tens of billions of US dollars annually (2, 3). There is currently no cure for rhinovirus infections, and the available treatments are only symptomatic. Rhinoviruses belong to the family Picornaviridae, genus Enterovirus, and are classified into species A, B, and C (4). Rhinoviruses A and B can belong to either major or minor groups, based on their utilization of intercellular adhesion molecule-1 (ICAM-1) or low-density lipoprotein receptor for cell entry (57). Type C rhinoviruses use CDHR3 as a receptor (8). Rhinovirus 14 belongs to the species rhinovirus B and uses ICAM-1 as a receptor. Receptors recognized by rhinoviruses and other enteroviruses can be divided into two groups based on their function in the infection process (9). Attachment receptors such as DAF, PSGL1, KREMEN1, CDHR3, and sialic acid enable the binding and endocytosis of virus particles into cells (1013). In contrast, uncoating receptors including ICAM-1, CD155, CAR, and SCARB2 enable virus cell entry but also promote genome release from virus particles (5, 1416).
Virions of rhinoviruses are nonenveloped and have icosahedral capsids (17). Genomes of rhinoviruses are 7,000 to 9,000 nucleotide-long single-stranded positive-sense RNA molecules (1, 17). The rhinovirus genome encodes a single polyprotein that is co- and posttranslationally cleaved into functional protein subunits. Capsid proteins VP1, VP3, and VP0, originating from one polyprotein, form a protomer, 60 of which assemble into a pseudo-T = 3 icosahedral capsid. To render the virions mature and infectious, VP0 subunits are cleaved into VP2 and VP4 (18, 19). VP1 subunits form pentamers around fivefold symmetry axes, whereas subunits VP2 and VP3 form heterohexamers centered on threefold symmetry axes. The major capsid proteins VP1 through 3 have a jelly roll -sandwich fold formed by two -sheets, each containing four antiparallel -strands, which are conventionally named B to I (2022). The two -sheets contain the strands BIDG and CHEF, respectively. The C termini of the capsid proteins are located at the virion surface, whereas the N termini mediate interactions between the capsid proteins and the RNA genome on the inner surface of the capsid. VP4 subunits are attached to the inner face of the capsid formed by the major capsid proteins. The surfaces of rhinovirus virions are characterized by circular depressions called canyons, which are centered around fivefold symmetry axes of the capsids (21).
The VP1 subunits of most rhinoviruses, but not those of rhinovirus 14, contain hydrophobic pockets, which are filled by molecules called pocket factors (17, 21, 23, 24). It has been speculated that pocket factors are fatty acids or lipids (25). The pockets are positioned immediately below the canyons. The exposure of rhinoviruses to acidic pH induces expulsion of the pocket factors, which leads to the formation of activated particles and genome release (17, 2632). The activated particles are characterized by capsid expansion, a reduction in interpentamer contacts, the release of VP4 subunits, externalization of N termini of VP1 subunits, and changes in the distribution of RNA genomes (17, 2629, 33, 34). Artificial hydrophobic compounds that bind to VP1 pockets with high affinity inhibit infection by rhinoviruses (35, 36).
ICAM-1 is an endothelial- and leukocyte-associated protein that stabilizes cellcell interactions and facilitates the movement of leukocytes through endothelia (37). ICAM-1 can be divided into an extracellular amino-terminal part composed of five immunoglobulin domains, a single transmembrane helix, and a 29-residuelong carboxyl-terminal cytoplasmic domain. The immunoglobulin domains are characterized by a specific fold that consists of seven to eight -strands, which form two antiparallel -sheets in a sandwich arrangement (3840). The immunoglobulin domains of ICAM-1 are stabilized by disulfide bonds and glycosylation (3841). The connections between the immunoglobulin domains are formed by flexible linkers that enable bending of the extracellular part of ICAM-1. For example, the angle between domains 1 and 2 differs by 8 between molecules in distinct crystal forms (38, 42). As a virus receptor, ICAM-1 enables the virus particles to sequester at the cell surface and mediates their endocytosis (5). The structures of complexes of rhinoviruses 3, 14, and 16, and CVA21 with ICAM-1 have been determined to resolutions of 9 to 28 (4246). It was shown that ICAM-1 molecules bind into the canyons at the rhinovirus surface, approximately between fivefold and twofold symmetry axes (4246). ICAM-1 interacts with residues from all three major capsid proteins. It has been speculated that the binding of ICAM-1 triggers the transition of virions of rhinovirus 14 to activated particles and initiates genome release (45, 47). However, the limited resolution of the previous studies prevented characterization of the corresponding molecular mechanism.
Here, we present the cryo-electron microscopy (cryo-EM) reconstruction of the rhinovirus 14 virion, which contains resolved density of octanucleotide segments of the RNA genome that interact with VP2 subunits. Furthermore, we show that the binding of ICAM-1 to rhinovirus 14 induces changes in its capsid and genome, which are required for subsequent virus activation and genome release at acidic pH.
The cryo-EM reconstruction of rhinovirus 14 virion was determined to a resolution of 2.6 (Fig. 1 A and B and SI Appendix, Fig. S1 and Table S1). The map enabled building of the structure of capsid proteins VP1 through 4 except for residues 1 to 17 and 290 to 293 of VP1, 1 to 6 of VP2, and 1 to 28 of VP4 (Fig. 1B). In addition to the capsid, the cryo-EM map contains resolved density corresponding to octanucleotide segments from the RNA genome (Fig. 1 AC). The quality of the map enabled building of the RNA structure; however, the nucleotide sequence could not be identified. The base of the second nucleotide from the 5 end of the RNA segment is flipped out from the RNA helix and forms a stacking interaction with the sidechain of Trp38 of VP2 (Fig. 1D and SI Appendix, Fig. S2). The residue Trp38 is conserved in the VP2 subunits of numerous picornaviruses, including polioviruses, rhinoviruses 2 and 16, coxsackievirus B3, coxsackievirus A9, and coxsackievirus A21 (SI Appendix, Fig. S3). Virion structures of these viruses contain disk-like densities that stack onto the tryptophane side chains, which were hypothesized to belong to a base of guanine nucleotide (23, 24, 4851). The structure of the RNA segment in the virion of rhinovirus 14 provides evidence that the previous speculations about the densities stacking onto Trp38 side chains were correct. To be consistent with the previous structures, we modeled the second nucleotide of the RNA segment in rhinovirus 14 as guanine (Fig. 1D and SI Appendix, Fig. S2). The stacking interaction between Trp38 and the base of the second nucleotide is the only direct contact between the RNA segment and the capsid (Fig. 1 A, C, and D and SI Appendix, Fig. S2).
Structure of virion of rhinovirus 14 contains resolved density corresponding to octanucleotides from its RNA genome. (A) Surface representation of cryo-EM of reconstruction of virion of rhinovirus 14 with front half of the particle removed to show internal structure. Density corresponding to VP1 is shown in blue, VP2 in green, VP3 in red, VP4 in yellow, and RNA segments in pink. Borders of a selected icosahedral asymmetric unit are indicated with a dashed triangle and positions of selected twofold, threefold, and fivefold symmetry axes are represented by an oval, triangle, and pentagon, respectively. (Scale bar, 5 nm.) (B) Cartoon representation of icosahedral asymmetric unit of rhinovirus 14 viewed from the inside of the capsid. The color coding of individual virus components is the same as in A. Positions of twofold, threefold, and fivefold symmetry axes are represented by an oval, triangle, and pentagon, respectively. (C) Two RNA octanucleotides that interact with each other and with VP2 subunits. Protein and RNA coloring is the same as in A. The cryo-EM density of the RNA segments is shown as a pink semitransparent surface. RNA bases and side chains of Trp38 of VP2 are shown in stick representation, in orange, and indicated with black arrowheads. The position of a twofold symmetry axis is indicated with an oval. (D) Detail of stacking interaction between Gua2 from RNA segment and Trp38 of VP2. The cryo-EM densities of RNA and protein are shown as semitransparent surfaces in pink and gray, respectively. (E) Interaction between N terminus of VP1 and genome. Capsid proteins are shown in cartoon representation with the same coloring as in A. Cryo-EM densities of individual proteins are shown as semitransparent surfaces colored according to the chain they belong to. The density of the RNA genome is shown in gray. The blue arrow indicates the contact between the N terminus of VP1 and the genome. The position of Thr17, the first modeled residue from the N terminus of VP1, is indicated.
Each RNA segment in the virion of rhinovirus 14 is associated with one protomer of capsid proteins VP1 through 4 (Fig. 1 A and B). The RNA is positioned next to a twofold axis, and two of the oligonucleotides interact with each other (Fig. 1 AC). To enable WatsonCrick base pairing between the two segments, the nucleotide at position seven was modeled as uracil and the nucleotide at position eight as adenine (Fig. 1C). Due to the constraints imposed by the interaction of the RNA with the capsid and of the interactions of two segments of the RNA with each other, the sequence of the oligonucleotide was built as 5 UGUUUUUA 3. Nevertheless, other sequences that fulfill the interaction conditions are equally possible.
The reconstruction of virion of rhinovirus 14 provides evidence that the N termini of VP1 subunits interact with the RNA genome (Fig. 1E). A similar function of the N terminus of VP1 was described previously in rhinovirus 2 (26). None of the interactions between the genomes and capsids of enteroviruses identified to date are sequence specific. Nevertheless, even the nonspecific interactions of the N termini of VP1 and Trp38 of VP2 with the RNA may enable the packaging of the enterovirus genome into a particle. Compounds that could prevent the RNAcapsid interaction by, for instance covering the side chain of Trp38, may interfere with the production of infectious virions.
The complex of rhinovirus 14 with ICAM-1 was prepared by mixing the components in phosphate-buffered saline (PBS) of pH 7.4 and incubating them at 34 C for 30 min (Fig. 2). The temperature was chosen to mimic that in the human upper respiratory tract (52). The binding of ICAM-1 did not induce the formation of activated particles or the genome release of rhinovirus 14 (Fig. 2 AC). This is in agreement with previous results showing that the ability of ICAM-1 to trigger genome release depends on temperature, receptor concentration, and buffer composition (47, 53). Experiments showing that ICAM-1 could induce genome release were performed in solutions with nonphysiological concentrations of salts, which may have destabilized the virus particles (42, 46, 47). Since enteroviruses have to deliver their genomes into the cytoplasm of a host cell to initiate infection, it would be detrimental if they released their genomes immediately upon binding to receptors at the cell surface. Our results show that rhinovirus 14 is stable when bound to ICAM-1 under native-like conditions (Fig. 2C), and the induction of genome release requires exposure to acidic pH in endosomes.
Binding of ICAM-1 to rhinovirus 14 is required for efficient genome release at acidic pH. (AD) Electron micrographs (Top), reference-free two-dimensional class averages (Center), and 3D reconstructions (Bottom) of rhinovirus 14 and rhinovirus 14ICAM-1 complex at neutral and acidic pH. The cryo-EM reconstructions are rainbow colored based on the distance of the particle surface from its center. Names above the reconstructions indicate the types of particles. Percentages below the reconstructions indicate the relative abundance of each type of particle in the sample. (A and B) Rhinovirus 14 at pH 7.4 (A) and 6.2 (B). (C and D) Rhinovirus 14ICAM-1 complex at pH 7.4 (C) and 6.2 (D). The top row of D contains two parts of a micrograph to show all types of particles present in the sample. Black arrowheads indicate virions, yellow indicate empty particles, cyan indicate rhinovirus 14ICAM-1 complex, green indicate activated particles, and red indicate open particles. (Scale bars, 30 nm.) (EH) Conformational changes to VP4 are induced by ICAM-1 binding but not by acidic pH. Surface representations of cryo-EM reconstructions showing the inner faces of capsids. The surfaces are color coded according to the capsid proteins with VP1 in blue, VP2 in green, VP3 in red, and VP4 in yellow. Only particles containing VP4 are shown. Positions of icosahedral symmetry axes are indicated with a black pentagon, triangle, and oval for fivefold, threefold, and twofold, respectively. Black arrows highlight C termini of VP4.
The structure of rhinovirus 14 in complex with the soluble ectodomain of ICAM-1 was determined to a resolution of 2.4 using cryo-EM and single-particle reconstruction (Fig. 3 andSI Appendix, Fig. S1 and Table S1). Domain 1 of ICAM-1 bound at the surface of rhinovirus 14 was resolved to a resolution of 2.6 (Fig. 3A and SI Appendix, Fig. S1). Levels of density in the map region corresponding to the domain 1 of ICAM-1 are similar to those in the capsid of the virus, indicating full occupancy of the receptors at the virus surface. Domains 2 and 3 of ICAM-1 are resolved to a resolution of 6 , and domains 4 and 5 are not visible in the cryo-EM reconstruction (Fig. 3A). The low resolution of the region of the map corresponding to domains 2 and 3 and the lack of density for domains 4 and 5 are probably caused by movements of those domains relative to domain 1, which is anchored at the virus surface (Fig. 3A) (38, 42, 43). In agreement with previous studies, domain 1 of ICAM-1 binds into the canyon of rhinovirus 14, approximately in the middle between fivefold and twofold symmetry axes (Fig. 3 B and C) (42, 43, 45). Previous studies of the interactions of rhinoviruses 3, 14, and 16 with ICAM-1 were limited to a resolution of 9.5 or lower (4246). The interpretation of the macromolecular interactions relied on the fitting of high-resolution structures, determined by X-ray crystallography, to the cryo-EM maps of the complex (4246). Therefore, the changes in the structures of the receptor and virus, induced upon their binding, could not be identified. Here, we show that ICAM-1 is wedged 3.4 deeper into the canyon and rotated 7.6 clockwise, when looking along the long axis of domain 1 toward the virus surface, relative to the structure reported previously (Fig. 4 and SI Appendix, Fig. S4A) (42, 43, 45). The interaction interface between ICAM-1 and rhinovirus 14 has a buried surface area of 1,500 2. The core of domain 1 of ICAM-1 is formed by -sheets ABED and GFC (Fig. 4 AC) (11, 54). Residues from the loops BC, DE, and FG and strands B, C, D, E, F, and G of ICAM-1 interact with rhinovirus 14 (SI Appendix, Fig. S4B). The mode of attachment of rhinovirus 14 to ICAM-1 is characteristic for uncoating receptors that bind to enterovirus canyons (5, 9, 14, 15). The uncoating receptors of enteroviruses, including ICAM-1, CD155, and CAR, have elongated molecules formed by domains with an immunoglobulin fold, which enables their insertion into the canyons of enterovirus particles (5, 14, 15).
Structure of rhinovirus 14 in complex with ICAM-1. (A) Surface representation of cryo-EM reconstruction of rhinovirus 14ICAM-1 complex color coded to distinguish individual proteins. Density corresponding to VP1 is shown in blue, VP2 in green, VP3 in red, and ICAM-1 in light magenta. Positions of selected icosahedral symmetry axes are indicated by a pentagon for fivefold, triangle for threefold, and an oval for twofold. The white triangle indicates the border of a selected icosahedral asymmetric unit. The yellow dashed rectangle indicates borders of the area shown in detail in B. (B) Roadmap projection showing residues forming the outer surface of rhinovirus 14 capsid (Left) and residues of domain 1 of ICAM-1 facing toward the virus (Right). Coloring is the same as in A. Residues involved in virusreceptor interactions are shown in bright colors. The polar angles and indicate positions at the capsid surface. (C) Schematic representation of ICAM-1. D1 to D5 indicate extracellular immunoglobulin domains; TM, transmembrane domain; cyt, cytoplasmic domain. Disulfide bridges (S-S) stabilizing the immunoglobulin domains are indicated. Red dashes highlight the binding site for rhinovirus 14. The ectodomain used in this study to determine the rhinovirus 14ICAM-1 interactions included residues 1 to 453.
Conformational changes associated with binding of rhinovirus 14 to ICAM-1. (A) Cartoon representation of icosahedral asymmetric unit of rhinovirus 14 in complex with ICAM-1. The VP1 subunit is shown in blue, VP2 in green, VP3 in red, VP4 in yellow, and domain 1 of ICAM-1 in magenta. The capsid proteins from virion of rhinovirus 14 are superimposed onto those of the rhinovirus 14ICAM-1 complex and are shown in white. The binding of ICAM-1 to rhinovirus 14 induces a 1.7 tilt of VP1 toward VP2 and VP3, which results in a narrowing of the canyon relative to the virion structure. Least-squares planes fitted to VP1 are shown to highlight the rotation of VP1. (B) Conformational change of FG loop of ICAM-1, shown in magenta, is required to prevent clashes with FG loop of VP3, shown in red. The native structure of ICAM-1 clashing with VP2 is shown in gray. (C) Cartoon representation of structure of ICAM-1 bound to rhinovirus 14 with side chains of cysteine residues shown in stick representation with red carbon atoms and yellow sulfur atoms. (D and E) Detail of disulfide bridge between Cys21 and Cys65 (D) and Cys25 and Cys69 (E). Cryo-EM density is shown as a semitransparent blue surface.
The binding of rhinovirus 14 to ICAM-1 is accompanied by the local restructuring of domain 1 of ICAM-1 and surface loops of capsid proteins, as well as by overall changes in the structure of the rhinovirus-14 capsid (Figs. 4 and 5 and SI Appendix, Figs. S5 and S6). VP1 subunits rotate 1.7 toward VP2 and VP3, which results in a contraction of the canyon (Fig. 4A). As a result, the capsid expands by 5 in diameter (Fig. 4A). The binding of ICAM-1 to rhinovirus 14 requires the bending of the FG loop of ICAM-1 8 toward the core of the immunoglobulin domain (Fig. 4B and SI Appendix, Fig. S5). This conformational change is necessary to prevent the clashing of the FG loop of ICAM-1 with residues 178 to 182 from the FG loop of VP3 of rhinovirus 14 (Fig. 4B). The conformational flexibility of the FG loop of ICAM-1 enables enlargement of its interaction interface with the capsid.
Changes of structure of C terminus of VP4 induced by ICAM-1 binding to rhinovirus 14. (A) Surface representation of cryo-EM reconstruction of capsid of rhinovirus 14 in complex with ICAM-1 viewed from inside the virion. Density corresponding to VP1 is shown in pale blue, VP2 in pale green, VP3 in pale red, and VP4 in semitransparent yellow. The structure of VP4 in the rhinovirus 14ICAM-1 complex is shown in cartoon representation in yellow, whereas the structure of VP4 in the virion of rhinovirus 14 is shown in magenta. The positions of selected icosahedral symmetry axes are indicated with a pentagon for fivefold, triangle for threefold, and oval for twofold. Borders of a selected icosahedral asymmetric unit are indicated with a dashed triangle. (B) Capsid structure of an empty particle of rhinovirus 14 containing pores around twofold symmetry axes and between twofold and fivefold symmetry axes through which VP4 may be released from the particle. (CF) Differences in structure of VP4 subunits in virion (C and E) and rhinovirus 14ICAM-1 complex (D and F). Capsid proteins are shown in cartoon representation. VP1 is shown in blue, VP2 in green, VP3 in red, VP4 in yellow, and RNA segments in pink. (C and E) Asn68 from C terminus of VP4 interacts with Asp11 and Arg12 of VP2 in virion of rhinovirus 14. The residues Asp11 and Arg12 are stabilized in position by the underlying loop of VP2 formed by residues 27 to 32 (highlighted in magenta). The side chain of Trp38 (highlighted in orange) forms a stacking interaction with Gua2 that is part of the resolved RNA segment positioned next to a twofold axis. (D and F) Binding of rhinovirus 14 to ICAM-1 induces conformational changes of virus capsid that include movement of residues 27 to 32 of VP2 toward particle center, which prevents interaction of C terminus of VP4 with residues Asp11 and Arg12 of VP2. The C terminus of VP4 acquires a new conformation, which covers the side chain of Trp38 of VP2 and blocks its interaction with RNA.
The structures of domain 1 of ICAM-1 determined to date contain the disulfide bonds Cys21 and Cys65 and Cys25 and Cys69 (3840). Residues Cys25 and Cys69 are located in the vicinity of the virus surface when ICAM-1 binds to rhinovirus 14 (Fig. 4 CE). Cys69 is part of the FG loop, whereas Cys25 is part of the BC loop (Fig. 4 C and E). The density connecting Cys25 and Cys69 of ICAM-1 in the complex with rhinovirus 14 is much weaker than that connecting Cys21 to Cys65 (Fig. 4 C and E). However, the positions of the two cysteines in the cryo-EM density map are consistent with the linkage of their side chains by a disulfide bond (Fig. 4 C and E) (55). Furthermore, mass spectrometry analysis of ICAM-1 molecules from the complex with rhinovirus 14 did not detect any peptides containing free Cys25 and Cys69 (SI Appendix, Fig. S7). However, peptides containing free cysteines were observed after the reduction of the disulfide bonds by dithiothreitol (DTT). This provides evidence that Cys25 and Cys69 of ICAM-1 in complex with rhinovirus 14 are linked by a disulfide bond. The low values of cryo-EM density may be caused by a higher flexibility of this part of ICAM-1, as indicated by the lower resolution than in the core of domain 1 (Fig. 3A). The binding of rhinovirus 14 to ICAM-1 also induces structural changes in the virus proteins. Residues 154 to 162 from the DE loop of VP1 shift 2 toward the core of the subunit (SI Appendix, Fig. S6). This movement helps to accommodate ICAM-1 in the canyon of rhinovirus 14.
The interaction between ICAM-1 and rhinovirus 14 is formed by 36 residues from domain 1 of ICAM-1 and 31, 8, and 13 residues of VP1, VP2, and VP3 of rhinovirus 14, respectively (Fig. 3B). The specificity of the interaction is controlled by a combination of the complementarity of the electrostatic interactions, a network of hydrogen bonds, and the positions of patches of hydrophobic interfaces. There are salt bridges between Lys77 of ICAM-1 and Glu210 from VP1 and Lys39 of ICAM-1 and C-terminal carboxyl group of Glu236 from VP3 (SI Appendix, Fig. S8 A and B). The interface includes a network of 18 hydrogen bonds. Furthermore, the resolution of the cryo-EM reconstruction of the complex was sufficient to enable the identification of water molecules, some of which form hydrogen bonds with both ICAM-1 and rhinovirus 14 and thus mediate interactions between the receptor and virus (SI Appendix, Fig. S8 C and D). For example, the amino group from the side chain of Lys29 of ICAM-1 interacts with two water molecules that form hydrogen bonds with side chains of Thr105 of VP3 and Asn68 of ICAM-1 (SI Appendix, Fig. S8D). It has been shown previously that the mutation of Thr75 of ICAM-1 to Ala reduces the efficiency of binding of rhinovirus 14 by more than 50% (54, 56). No ions were identified at the rhinovirus 14ICAM-1 interface (7, 57).
The interface between rhinovirus 14 and ICAM-1 contains complementary patches of hydrophobic interactions (SI Appendix, Fig. S8 E and F). Previous studies have shown that most mutations of Pro70 from the FG loop of ICAM-1 prevent the binding of rhinovirus 14 (54, 56). We show that Pro70 fits into a small hydrophobic pocket formed by Pro178, Phe86, and Thr180 of VP3 (SI Appendix, Fig. S8F). Fitting Pro70 of ICAM-1 into the hydrophobic cavity in VP3 requires movement and restructuring of the FG loop of ICAM-1 (Fig. 4 C and E and SI Appendix, Fig.S8F). Another residue of ICAM-1 that is critical for the binding of rhinovirus 14 is Leu30 (54, 56). In the complex, the side chain of Leu30 is situated between the side chains of Ile215 and Val217 from VP1, which form a hydrophobic pocket for the leucine side chain (SI Appendix, Fig. S8E). This explains why the mutation of Leu30 to Ser eliminates the binding of ICAM-1 to rhinovirus 14 (54, 56).
The structure of the C terminus of VP4 subunit in the rhinovirus 14ICAM-1 complex differs from that in the native virion (Figs. 2 EH and 5). The two structures of VP4 subunits are similar for residues 29 to 57, with rmsd of C atoms of the corresponding residues of 0.44 . However, residues 58 to 65 of VP4 extend toward a threefold symmetry axis of the capsid in the native virion, whereas the same residues point toward a twofold symmetry axis in the rhinovirus 14ICAM-1 complex (Fig. 5A). The movement of the C terminus of VP4 is induced by conformational changes of the major capsid proteins, which are triggered by ICAM-1 binding to the capsid. In the rhinovirus-14 virion, the C-terminal carboxyl group of Asn68 from VP4 forms a salt bridge with the side chain of Arg12 of VP2 (Fig. 5 C and E). Additionally, the side chain of Asn68 forms two hydrogen bonds with the side chain of Asp11 of VP2 (Fig. 5 C and E). As discussed above, the binding of ICAM-1 to rhinovirus 14 induces a rotation of VP1 toward VP2 and VP3 (Fig. 4A). These movements of capsid proteins bring residues 27 to 33 from the N terminus of VP2 into the space that is occupied by Arg12 of VP2 in the native virion (Fig. 5 D and F). This frees the C terminus of VP4 from the interaction with Arg12 of VP2 and probably enables its translocation toward a twofold axis (Fig. 5 A, D, and F). The restructuring of the C-terminal part of VP4 to point toward a twofold symmetry axis prepares the protein for release through either of the holes that form at and next to the twofold symmetry axes upon particle activation (Fig. 5 A and B) (29, 34).
The binding of ICAM-1 into the canyon of rhinovirus 14 induces the relocation of the C terminus of VP4 toward a twofold symmetry axis of the capsid (Fig. 5A). The movement of the C terminus of VP4 uncovers a patch of positively charged residues at the inner face of the capsid, adjacent to a threefold symmetry axis (Fig. 6 A and B). The positively charged surface attracts genomic RNA, which is represented in the cryo-EM reconstruction as a cylindrical appendage emanating from the spherical genome density filling the center of the virus particle (Fig. 6 CF). This indicates that parts of the RNA genome in various conformations interact nonspecifically with the positively charged regions of the capsid. Furthermore, the N termini of VP2 subunits probably interact with the RNA density positioned on a threefold axis (Fig. 6F). The C terminus of VP4 positioned next to a twofold axis of the capsid covers the side chain of Trp38 of VP2, which in the native virion forms a stacking interaction with a nucleotide from the RNA genome (Figs. 1D and 5 E and F and SI Appendix, Fig. S2). The loss of Trp38RNA contact relaxes the ordering of segments of the RNA genome that interact with the capsid around twofold symmetry axes in native virions (Fig. 6 G and H). A density corresponding to RNA at the twofold axis in the rhinovirus 14ICAM-1 complex does not have resolved features (Fig. 6 G and H). The interactions of the N termini of VP1 subunits with the RNA genome remain preserved even after the binding of rhinovirus 14 to ICAM-1 (Fig. 1E). The binding of rhinovirus 14 to ICAM-1 induces a rearrangement of the RNA genome, which may play a role in particle activation, as discussed below.
Conformational changes of capsid of rhinovirus 14 that are induced by binding to ICAM-1 trigger redistribution of RNA genome in the particle. (A and B) Detail of inner capsid surface around threefold symmetry axis of virion (A) and rhinovirus 14ICAM-1 complex (B). The surfaces are colored according to charge. (C and D) Electron densities of central slices of cryo-EM reconstructions of virion (C) and rhinovirus 14ICAM-1 complex (D) with a thickness of 1 . White represents high density values. Density representing ICAM-1 is highlighted in magenta. Positions of icosahedral symmetry axes are indicated with an oval, triangle, and pentagon for twofold, threefold, and fivefold axes, respectively. Black arrows in D point toward densities on threefold symmetry axes, which are not present in the virion. Red squares indicate positions of details shown at higher magnification in E and F. (Scale bar, 10 nm.) (E and F) Details of cryo-EM density distribution at the inner face of the capsid on a threefold symmetry axis. Capsid proteins are shown in cartoon representation with VP1 in blue, VP2 in green, VP3 in red, and VP4 in yellow. Cryo-EM density is shown as a semitransparent gray surface. Positions of the first resolved residues from the N termini of VP2 subunits are indicated. E and F show sections of particles with a thickness of 20 . (G and H) Comparison of structures of RNA genome interacting with VP2 subunits in virion (G) and rhinovirus 14ICAM-1 complex (H). The virion contains resolved cryo-EM density corresponding to octanucleotides (G). In contrast, there is a featureless density in the rhinovirus 14ICAM-1 complex (H). Capsid proteins are shown in cartoon representation, colored as in E and F.
The binding of ICAM-1 to rhinovirus 14 triggers a cascade of structural changes that prepare the particle for activation and subsequent genome release (Fig. 7 A and B). The rotation of VP1 subunits results in a narrowing of the canyon and transmits the conformational changes to the inside of the capsid (Fig. 7 C and D). C termini of VP4 subunits reposition toward twofold symmetry axes, where they are optimally poised for externalization upon particle activation (Fig. 7 E and F). The conformational change to C termini of VP4 subunits uncovers patches of positively charged residues that attract genomic RNA toward threefold symmetry axes of the capsid (Fig. 7 G and H). The same conformational change prevents the interaction of Trp38 from VP2 with bases from the ordered RNA segments of the genome positioned next to twofold symmetry axes of the capsid (Fig. 7 G and H). Both of these effects result in reorganization of the RNA genome (Fig. 7 G and H). These changes in the capsid and genome structure of rhinovirus 14 induced by ICAM-1 binding are required for efficient genome release at acidic pH (Fig. 2). All in all, 90% of virions of rhinovirus 14 exposed to pH 6.2 remained in their native conformation, whereas the remaining particles were empty (Fig. 2 A, B, E, and F and SI Appendix, Table S1). The structures of rhinovirus 14 in their native conformation and empty capsids at acidic pH were determined to resolutions of 2.8 and 3.9 , respectively (Fig. 2 A and B and SI Appendix, Table S1). The exposure of rhinovirus 14 to acidic pH did not induce structural changes in VP4 subunits (Fig. 2 E and F). In contrast, the exposure of rhinovirus 14ICAM-1 complex to pH 6.2 resulted in activation and genome release from 95% of particles (Fig. 2 C, D, E, and H). The structure of the activated particle was determined to a resolution of 4.0 , empty particle to 3.9 , and open particle to 22 (Fig. 2 C and D and SI Appendix, Table S1). The changes in the capsid and genome of rhinovirus 14, which were induced by ICAM-1 binding, may lower the energy barrier of particle activation so that it can be overcome by random fluctuations in particle structure due to thermal motions termed capsid breathing (58, 59). This provides a putative explanation of how the reduction of capsid dynamics by antiviral compounds, which target VP1 pockets (59), blocks the activation and genome release of enteroviruses.
Overview of structural changes to rhinovirus 14 induced by binding of ICAM-1 that prepare the virus for activation and genome release. (A) Native virion diffuses toward cell membrane (green ribbon) decorated with ICAM-1 molecules (blue sticks with dark blue heads representing domain 1). The virus particle is represented by a central slice with the electron density of VP1 shown in light blue, VP2 in light green, VP3 in light red, VP4 in yellow, the genome in purple, and resolved RNA segments in pink. (B) Rhinovirus 14 is endocytosed by the cell after binding to ICAM-1. (CH) Sequence of structural changes in virion induced by binding to ICAM-1. C, E, and G represent native virion, whereas D, F, and H show rhinovirus 14ICAM-1 complex. (C and D) Binding of ICAM-1 induces rotation of VP1 subunit toward VP2 and VP3. Virus components are colored as in A. (E and F) ICAM-1 binding induces disruption of interactions between C terminus of VP4 and N terminus of VP2. C terminus of VP4 repositions from a threefold symmetry axis (indicated with a triangle) toward a twofold symmetry axis (oval). (G and H) Movements of C termini of VP4 subunits uncover positively charged residues around twofold symmetry axes, which attract negatively charged RNA genome. Furthermore, the C terminus of VP4 in the altered conformation covers Trp38 of VP2 and prevents its specific interaction with structured segments of the RNA genome, which relaxes the structure of RNA adjacent to the twofold symmetry axis.
The sample of complex of rhinovirus 14 with ICAM-1 exposed to acidic pH contained 7% empty particles missing a pentamer of capsid protein protomers (Fig. 2D). Open particles were previously speculated to enable enterovirus genome release (33). The expulsion of pentamers of capsid proteins results in the formation of a large hole in the capsid, which enables the diffusion of the RNA genome from the capsid within a microsecond (33, 60). The rapid release of a genome may be connected to its subsequent transport across the endosome membrane into the cytoplasm (61, 62).
Structural characterization of the rhinovirus 14ICAM-1 complex at atomic resolution provides detailed information about the conformational changes of both the receptor and virus that are required for its binding. Additionally, it provides insight into the structural changes of the virus that enable subsequent particle activation and genome release. In combination, these results provide the basis for the design of compounds that block enterovirus infection.
The extracellular part of ICAM-1 containing domains D1 to D5 was produced using the MutiBac system (Geneva Biotech). The full-length gene of ICAM-1 was a gift from Timothy Springer (Harvard Medical School, Boston, MA) (Addgene plasmid No. 8632; http://addgene.org/8632; RRID: Addgene 8632). The sequence encoding domains D1 to D5, the secretion signal peptide, and the C-terminal 10 His-tag were inserted into the pACEBac1 vector at the restriction site BamHI. The recombinant bacmid with the target sequence was prepared by recombination in DH10EMBacY Escherichia coli cells (Geneva Biotech). The recombinant baculovirus was prepared by transfecting SF9 cells with the recombinant bacmid. A total of 250 mL of the culture of SF9 cells were infected with the recombinant baculovirus and incubated for 96 h at 27 C with 120 rpm shaking. The produced protein was secreted into the medium. Cells and cell debris were pelleted by centrifugation at 20,000 g at 4 C for 15 min. The supernatant was filtered through a 0.2 m filter (Corning) and loaded into a HisTrap column (GE Healthcare) equilibrated in PBS (137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, 1.8 mM KH2PO4, pH 7.4). Most of the impurities were removed by washing with PBS containing 70 mM imidazole. His-tagged ICAM-1 D1 to D5 was eluted using PBS with 500 mM imidazole. The eluted protein was buffer exchanged into PBS using 30 kDa cutoff centrifugal concentrators (Millipore, Merck). The target protein was further purified by size-exclusion chromatography, using a HiLoad 16/600 Superdex 200 pg column (GE Healthcare). Fractions containing ICAM-1 D1 to D5 were pooled and concentrated using centrifugal concentrators (Millipore, Merck) to a final concentration of 3.5 mg/mL.
Rhinovirus-14 strain 1059, obtained from ATCC, was propagated in HeLa-H1 (ATCC CRL195) cells cultivated in Dulbeccos modified Eagles medium enriched with 10% fetal bovine serum. HeLa cells grown to 100% confluence (100 tissue culture dishes, 150 mm diameter) were infected with a multiplicity of infection of 0.1. The infection was allowed to progress for 36 h until a complete cytopathic effect was observed. The media and cells were harvested and centrifuged at 15,000 g at 10 C for 30 min. The resulting pellet was subjected to three freezethaw cycles and resuspended in 5 mL PBS followed by homogenization in a Dounce tissue grinder. Cell debris was separated from the supernatant by centrifugation at 4,000 g for 30 min at 10 C. The virus-containing supernatant was combined with the media from infected cells. The virus particles were precipitated by the addition of PEG-8000 and NaCl to final concentrations of 12.5% (wt/vol) and 500 mM, respectively, and incubation overnight at 10 C. The precipitated virus was pelleted by centrifugation at 15,000 g at 10 C for 30 min. The pellet was resuspended in 20 mL PBS with 5 nm MgCl2. The sample was subjected to DNase (10 g/mL final concentration) and RNase (10 g/mL final concentration) treatment at room temperature for 30 min. Subsequently, trypsin was added to a final concentration of 0.5 g/mL, and the sample was incubated at 37 C for 15 min. EDTA (pH = 9.5) and Nonidet P-40 (Sigma-Aldrich) were added to final concentrations of 15 mM and 1% (vol/vol), respectively. The virus was pelleted through a 30% (wt/vol) sucrose cushion by centrifugation at 210,000 g in an Optima 80 ultracentrifuge (Beckman Coulter) using a Ti50.2 rotor. The pelleted virus particles were resuspended in 2 mL PBS, loaded on the top of a 60% (wt/wt) CsCl solution in PBS, and centrifuged at 160,000 g in an Optima 80 ultracentrifuge using a Beckman Coulter SW41Ti rotor at 10 C for 24 h. Opaque bands containing virus particles were harvested and subjected to buffer exchange in PBS using a Centricon Plus-70 centrifugal filter (Millipore) with a 100 kDa cutoff. The final concentration of rhinovirus 14 was 0.5 mg/mL. Purified virus was kept at 4 C.
The complex of rhinovirus 14 with ICAM-1 was prepared by mixing rhinovirus 14 with ICAM-1 at a molar ratio of 1:100 at pH 7.4 and incubating the mixture for 30 min at 34 C. Rhinovirus 14 and the rhinovirus 14ICAM-1 complex were transferred to acidic pH using DyeEx 2.0 (QIAGEN) spin columns containing PBS with pH 6.2. The samples were applied onto the columns and eluted by 1 min of centrifugation at 1,200 g. The samples were incubated at pH 6.2 at 34 C for 2.5 min, including the centrifugation step.
For vitrification, 3 L virus samples were applied onto a holey carbon-coated copper grid (R2/1, mesh 300, Quantifoil), blotted for 2 s, and vitrified by plunging into liquid ethane using a Vitrobot Mark IV (Thermo Fisher Scientific). Grids for virion reconstruction were prepared by vitrifying a virus sample with a concentration of 0.5 mg/mL. The grids (except for those with rhinovirus 14 at pH 6.2) were then transferred to a Titan Krios electron microscope, operating at 300 kV at cryogenic conditions, equipped with a Falcon III direct electron detector (Thermo Fisher Scientific). The illuminating beam was aligned for parallel illumination in NanoProbe mode. Low-dose imaging was used with a total dose of 84.7 e/ 2. Nominal magnification was set to 75,000, resulting in a calibrated pixel size of 1.063 . The dataset was recorded automatically using EPU software (Thermo Fisher Scientific) in fast acquisition mode, using large image shifts. The samples of rhinovirus 14 and rhinovirus 14ICAM-1 complex were recorded using five acquisitions per hole, nine holes per stage shift. The sample of rhinovirus 14ICAM-1 complex at pH 6.2 was recorded using seven acquisitions per hole, nine holes per stage shift. The exposure time was set to 1 s, and each micrograph was recorded as a movie containing 40 frames. The target defocus range was 0.5 to 2.4 m.
Electron micrographs from the sample of rhinovirus 14 at pH 6.2 were collected using a Talos Arctica electron microscope (Thermo Fisher Scientific), operated at 200 kV under cryogenic conditions, equipped with a Falcon III direct electron detector (Thermo Fisher Scientific). The illuminating beam was aligned for parallel illumination in NanoProbe mode. Low-dose imaging was used with a total dose of 34.1 e/ 2. Nominal magnification was set to 120,000, resulting in a calibrated pixel size of 1.22 . The dataset was recorded automatically using EPU software (Thermo Fisher Scientific) in fast acquisition mode, using five acquisitions per hole, nine holes per stage shift. The exposure time was set to 1 s, and each micrograph was recorded as a movie containing 40 frames. The target defocus range was 0.5 to 2.4 m.
The beam-induced movements within one micrograph were corrected with the software MotionCorr2 using 5 5 patches (63). The motion-corrected micrographs were dose weighted, and defocus values were estimated using the program gCTF (64). Using crYOLO box manager (65), 200 particles were boxed manually and used as a template for ab initio model training. The resulting crYOLO model was used to pick particles. The particles were extracted using Relion3.1 (66) with a box size of 546 px. The particles were binned to a box size of 150 150 px and subjected to reference-free two-dimensional classifications in Relion3.1 (66). Particles from classes exhibiting high-resolution features were used for de novo model calculation with imposed icosahedral symmetry, using stochastic gradient descent as implemented in Relion3.1 (66). The resulting three-dimensional (3D) volume was used as a starting model for autorefinement in Relion3.1. After initial autorefinement, 3D classification in Relion3.1 was performed, omitting the alignment step. Particles belonging to the best class were reextracted and recentered box-size 512 512 px for rhinovirus 14 particles without ICAM-1 and 546 546 px for the rhinovirus 14ICAM-1 complex. Reextracted particles were subjected to another round of autorefinement in Relion3.1. Particles were then sorted into nine optic groups. The optics groups were determined by the position of the image shift used for the acquisition, whereas all the acquisition areas from the same foil hole were considered as one optics group. Therefore, only large (interhole) image shifts were considered as separate optics groups. Magnification correction was performed using Relion3.1, followed by beam-tilt correction, and subsequently by the estimation of third- and fourth-order Zernike polynomials. The aberration-corrected particles were further subjected to per-particle defocus and astigmatism correction and estimation of the CTF envelope function (CTF B-factor fitting). The particles were subjected to autorefinement with imposed icosahedral symmetry. Ewald sphere correction was performed as implemented in relion_reconstruct.py in Relion3.1 (67). The resulting map was used for Bayesian polishing of particles with default parameters. The polished particles were used for 3D autorefinement and CTF refinement followed by another 3D autorefinement. Finally, Ewald sphere correction was performed. The final map was threshold masked, divided by a modulation transfer function, and B-factor sharpened using Relion3.1. Local resolutions were estimated using the program MonoRes implemented in the Scipion software package (68, 69). Map B-factor sharpening based on local resolution estimation from MonoRes was performed using the program LocalDeblur (70). The dataset of rhinovirus 14 with ICAM-1 exposed to pH 6.2 contained empty capsids missing pentamers. These were identified by 3D classification with C5 symmetry, using the complete empty capsid as an initial model. Subsequent 3D autorefinement was performed with C5 symmetry. Neither CTF refinement nor Bayesian polishing were applied to this subset of particles.
The electrostatic potential map from cryo-EM reconstruction was oriented to the standard 222 icosahedral crystallographic orientation. The origin of the map was moved from the 0,0,0, coordinate to the center of the particle using the program mapman (71). The map was normalized and converted to crystallographic space group P23 using the CCP4i software suite (72). The higher-symmetry space group was used to reduce the computational demands of the model refinement. Crystal structures of rhinovirus 14 (Protein Data Bank [PDB]: 4RHV) and domain 1 of ICAM-1 (PDB: 1IC1) were manually fitted into the cryo-EM maps using the program Chimera and refined with the tool Fit in map (73). The cryo-EM structure of an empty particle of rhinovirus 14 in complex with a Fab fragment of antibody (PDB: 5W3O) was used as a starting model for the building of activated and empty particles (74). The fitted models were subjected to multiple rounds of real-space refinement in Phenix (version dev-3765), reciprocal-space refinement in REFMAC5, combined with manual corrections in Coot 0.9 and ISOLDE (7578). Hydrogen atoms were taken into account during the real-space refinement, whereas they were ignored in the reciprocal-space refinement. Waters were added automatically by the program find waters in Coot 0.9 and validated manually. Model validation parameters were calculated using MolProbity server and EMringer as implemented in phenix (79, 80). The RNA octanucleotide sequence in the native virion of rhinovirus 14 was initially built using the program Coot and refined with restraints using the program ISOLDE (78). Structural comparisons were performed in Chimera (73). Hydrogen bonds, salt bridges, and residues involved in the binding interface and buried surface areas were calculated using the program PDBePISA (https://www.ebi.ac.uk/pdbe/pisa/). Roadmaps were produced using the program Rivem (81).
Purified samples of ICAM-1 and the complex of rhinovirus 14 with ICAM-1 were digested with alpha-lytic protease (EC 3.4.21.12, Sigma-Aldrich catalog No. A6362) for 2 h at 37 C with shaking at 700 rpm. Half of the volume of each sample was then reduced using 10 mM DTT (for 45 min at 57 C with shaking at 700 rpm). After adding polyethylene glycol to a final concentration of 0.001%, the peptides were extracted from the vials using 25% formic acid/acetonitrile (1:1 vol/vol mixture) and vacuum concentrated. The peptide mixture was subjected to liquid-chromatography-mass spectrometry (LC-MS)/MS analysis using a RSLCnano system (ThermoFisher Scientific) coupled to an Impact II Qq-Time-Of-Flight mass spectrometer (Bruker). Prior to LC separation, peptides were online concentrated in a trap column (100 m 30 mm) filled with 3.5 m X-Bridge BEH 130 C18 sorbent (Waters). The peptides were separated using an Acclaim Pepmap100 C18 column (3 m particles, 75 m 500 mm; ThermoFisher Scientific) by the following LC gradient program (mobile phase A: 0.1% formic acid in water; mobile phase B: 0.1% formic acid in 80% acetonitrile; 300 nl/min): the gradient elution started at 1% of mobile phase B and increased to 56% over the first 50 min, then increased linearly to 80% of mobile phase B over the next 5 min and remained at this state for the final 10 min. Equilibration of the trapping column and the column was done prior to sample injection into the sample loop. The analytical column outlet was directly connected to a CaptiveSpray nanoBooster ion source (Bruker). The nanoBooster was filled with acetonitrile. MS and MS/MS spectra were acquired in a data-dependent strategy with a 3 s cycle time. The mass range was set to 150 to 2,200 m/z and precursors were selected from 300 to 2,000 m/z. The acquisition speed of MS and MS/MS scans was 2 Hz and 4 to 16 Hz, respectively. The speed of MS/MS spectra acquisition was based on precursor intensity. The preprocessing of the mass spectrometric data including recalibration, compound detection, and charge deconvolution was carried out using DataAnalysis software (version 4.2 SR1; Bruker).
The obtained data were searched with an in-house Mascot search engine (version 2.4.1; Matrixscience) against a custom database involving the ICAM-1 sequence and cRAP entries (downloaded from https://www.thegpm.org/crap/). The database searches were done without enzyme specificity and with oxidation (M) as a variable modification. The mass tolerances for peptides and MS/MS fragments were 10 ppm and 0.1 Da, respectively. Only peptides with a statistically significant peptide score (P < 0.05) were considered, and the obtained MS/MS data were validated manually.
Multiple sequence alignment of capsid proteins of selected viruses from the family Picornaviridae was performed in the Clustal Omega server (82). The multiple sequence alignment was visualized in the software Jalview 2.11.1.3 (83).
The cryo-EM maps and coordinates were deposited under the following accession codes: virion of rhinovirus 14 at neutral pH Electron Microscopy Data Bank EMD-12171 and PDB 7BG6; rhinovirus 14ICAM-1 complex at neutral pH EMD-12172 and PDB 7BG7; rhinovirus 14 in native conformation at acidic pH EMD-12599 and PDB 7NUQ; empty particle of rhinovirus 14 at acidic pH EMD-12597 and PDB 7NUO; rhinovirus 14 in native conformation at acidic pH originating from complex with ICAM-1 EMD-12596 and PDB 7NUN; activated particle originating from complex with ICAM-1 at acidic pH EMD-12594 and PDB 7NUL; empty particle originating from complex with ICAM-1 at acidic pH EMD-12595 and PDB 7NUM; and open particle originating from complex with ICAM-1 at acidic pH EMD-12598.
We gratefully acknowledge the Cryo-Electron Microscopy and Tomography Core Facility of Central European Institute of Technology (CEITEC) supported by Ministry of Education, Youth, and Sports of the Czech Republic (MEYS CR) (Grant LM2018127). This research was carried out under the project CEITEC 2020 (Grant LQ1601), with financial support from the MEYS of the Czech Republic under National Sustainability Program II. This work was supported by the IT4I project (Grant CZ.1.05/1.1.00/02.0070) and funded by the European Regional Development Fund and the national budget of the Czech Republic via the Research and Development for Innovations Operational Programme (RDI-OP) as well as the MEYS via Grant LM2011033. The research leading to these results received funding from Czech Science Foundation Grant GX19-25982X to P.P.
Author contributions: D.H. and P.P. designed research; D.H., L.., A.A., and O.. performed research; D.H., T.F., M.G., O.., Z.Z., and P.P. analyzed data; and D.H., T.F., and P.P. wrote the paper.
The authors declare no competing interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at https://www.pnas.org/lookup/suppl/doi:10.1073/pnas.2024251118/-/DCSupplemental.
- ENCODE: Encyclopedia Of DNA Elements - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- 07.05.2010 - The Human Genome [ Coast To Coast AM ] - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- NOVA scienceNOW : 51 - Public Genomes, Algae Fuel, Mystery of the Gakkel Ridge, Yoky Matsuoka - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- Vincent T. - Genome (Club Remix) - [Preview] - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- Comparing The Human And Chimpanzee Genomes - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- Whole Genome Sequencing and Its Impact on Clinical Care - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- Carlos Bustamante -- "Reconstructing the Great Human Diasporas from Genome Variation Data" - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- 3 Sad Surprises: The Human Genome Project - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- The RFW interviews Genome - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- Science Bulletins: Scientists Peer Inside "Superbug" Genome - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- Genome : Live @ Smu's : June 3 2012 - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- Inoki Genome Federation - Genome 19 - 04 02 2012 - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- THE HUMAN GENOME MUSIC PROJECT - CHROMOSOME 1 - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- Genomic Medicine - Bruce Korf (2012) - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- Human Genome's 'Blockbuster' Potential Undervalued in Bid GSK vs HGSI - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- Announcing the Completion of the First Survey of the Entire Human Genome at the White House - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- DNA analysis Part I. Genomic Sequencing - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- The Genome Question: Moore vs. Jevons with Bud Mishra - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- Genome-Wide Association Studies - Karen Mohlke (2012) - Video [Last Updated On: September 7th, 2012] [Originally Added On: September 7th, 2012]
- New human genome research aids understanding of disease [Last Updated On: September 8th, 2012] [Originally Added On: September 8th, 2012]
- UNC Lineberger scientists lead definition of key lung cancer genome [Last Updated On: September 10th, 2012] [Originally Added On: September 10th, 2012]
- Illumina Announces Expedited Individual Genome Sequencing Service (IGS) [Last Updated On: September 11th, 2012] [Originally Added On: September 11th, 2012]
- Genome research given a boost with opening of bioscience facility [Last Updated On: September 11th, 2012] [Originally Added On: September 11th, 2012]
- Re-Imagining Our Genes: ENCODE Project Reveals Genome as an Information Processing System [Last Updated On: September 11th, 2012] [Originally Added On: September 11th, 2012]
- Illumina unveils upgraded genome sequence service [Last Updated On: September 12th, 2012] [Originally Added On: September 12th, 2012]
- US Personalized Cancer Genome Sequencing Market [Last Updated On: September 18th, 2012] [Originally Added On: September 18th, 2012]
- Yale maps “uncharted” genome regions [Last Updated On: September 18th, 2012] [Originally Added On: September 18th, 2012]
- Research and Markets: US Personalized Cancer Genome Sequencing Market [Last Updated On: September 19th, 2012] [Originally Added On: September 19th, 2012]
- 3Qs: New clues to unlocking the genome [Last Updated On: September 19th, 2012] [Originally Added On: September 19th, 2012]
- Oyster Genome Pries Open Mollusk Evolutionary Shell [Last Updated On: September 20th, 2012] [Originally Added On: September 20th, 2012]
- Bangladeshi scientist decodes genome of deadly fungus [Last Updated On: September 20th, 2012] [Originally Added On: September 20th, 2012]
- Oyster genome uncover the stress adaptation and complexity of shell formation [Last Updated On: September 20th, 2012] [Originally Added On: September 20th, 2012]
- The oyster genome reveals stress adaptation and complexity of shell formation [Last Updated On: September 20th, 2012] [Originally Added On: September 20th, 2012]
- Diseases of aging map to a few 'hotspots' on the human genome [Last Updated On: September 20th, 2012] [Originally Added On: September 20th, 2012]
- GnuBIO Awarded $4.5 Million in Funding from the National Human Genome Research Institute to Develop Lower Cost Genome ... [Last Updated On: September 20th, 2012] [Originally Added On: September 20th, 2012]
- Oyster genome mystery unravelled [Last Updated On: September 20th, 2012] [Originally Added On: September 20th, 2012]
- Devangshu Datta: What's in a genome [Last Updated On: September 20th, 2012] [Originally Added On: September 20th, 2012]
- Pacific Oyster Genome Shows Stress Adaptation And Complexity Of Shell Formation [Last Updated On: September 20th, 2012] [Originally Added On: September 20th, 2012]
- UNC Lineberger scientists lead cancer genome analysis of breast cancer [Last Updated On: September 24th, 2012] [Originally Added On: September 24th, 2012]
- Encoding the human genome [Last Updated On: September 24th, 2012] [Originally Added On: September 24th, 2012]
- Cancer genome analysis of breast cancer: Team identifies genetic causes and similarity to ovarian cancer [Last Updated On: September 24th, 2012] [Originally Added On: September 24th, 2012]
- Fungus genome map paves way for 'Snow White' jute variety [Last Updated On: September 24th, 2012] [Originally Added On: September 24th, 2012]
- New online, open access journal focuses on microbial genome announcements [Last Updated On: September 25th, 2012] [Originally Added On: September 25th, 2012]
- By Simply Sharing, Doctors Could Unlock the Genome's Potential [Last Updated On: September 25th, 2012] [Originally Added On: September 25th, 2012]
- Forget the Cloud—Knome Offers Genome Analysis in a Box [Last Updated On: September 28th, 2012] [Originally Added On: September 28th, 2012]
- BGI@CHOP Joint Genome Center to Offer Clinical Next-Generation Sequencing Services [Last Updated On: September 28th, 2012] [Originally Added On: September 28th, 2012]
- Holy Bat Virus! Genome Hints At Origin Of SARS-Like Virus [Last Updated On: September 29th, 2012] [Originally Added On: September 29th, 2012]
- Community Fundraising Effort Helps Researchers Sequence Parrot Genome [Last Updated On: September 29th, 2012] [Originally Added On: September 29th, 2012]
- UMass Med professors are sleuths of the genome [Last Updated On: September 30th, 2012] [Originally Added On: September 30th, 2012]
- Knome Introduces the knoSYS™100; First Plug-and-Play Human Genome Interpretation System [Last Updated On: September 30th, 2012] [Originally Added On: September 30th, 2012]
- First large scale trial of whole-genome cancer testing for clinical decision-making reported [Last Updated On: October 1st, 2012] [Originally Added On: October 1st, 2012]
- Should You Get Your Genome Mapped? [Last Updated On: October 1st, 2012] [Originally Added On: October 1st, 2012]
- Surprising differences between apples and pears [Last Updated On: October 2nd, 2012] [Originally Added On: October 2nd, 2012]
- 50-Hour Whole Genome Sequencing Provides Rapid Diagnosis for Children With Genetic Disorders [Last Updated On: October 3rd, 2012] [Originally Added On: October 3rd, 2012]
- A map of rice genome variation reveals the origin of cultivated rice [Last Updated On: October 3rd, 2012] [Originally Added On: October 3rd, 2012]
- Genome analysis promises hope for breast cancer patients [Last Updated On: October 3rd, 2012] [Originally Added On: October 3rd, 2012]
- Genome Alberta Welcomes Alberta Minister of Enterprise and Advanced Education, Stephen Khan and Federal Minister of ... [Last Updated On: October 3rd, 2012] [Originally Added On: October 3rd, 2012]
- Fifty-hour whole genome sequencing provides rapid diagnosis for children with genetic disorders [Last Updated On: October 3rd, 2012] [Originally Added On: October 3rd, 2012]
- Will Low-Cost Genome Sequencing Open 'Pandora's Box'? [Last Updated On: October 3rd, 2012] [Originally Added On: October 3rd, 2012]
- Genome testing could help individualize treatments [Last Updated On: October 3rd, 2012] [Originally Added On: October 3rd, 2012]
- Would you get your genome tested? [Last Updated On: October 3rd, 2012] [Originally Added On: October 3rd, 2012]
- The Genome — a Pandora's Box? [Last Updated On: October 4th, 2012] [Originally Added On: October 4th, 2012]
- Fast genome test could help sick newborns [Last Updated On: October 4th, 2012] [Originally Added On: October 4th, 2012]
- In-Depth Genome Analysis Moves Toward The Hospital Bed [Last Updated On: October 5th, 2012] [Originally Added On: October 5th, 2012]
- Your Verdict On Getting A Genome Test? Bring It On [Last Updated On: October 6th, 2012] [Originally Added On: October 6th, 2012]
- Genome-wide study identifies 8 new susceptibility loci for atopic dermatitis [Last Updated On: October 7th, 2012] [Originally Added On: October 7th, 2012]
- Genome-wide study identifies eight new susceptibility loci for atopic dermatitis [Last Updated On: October 7th, 2012] [Originally Added On: October 7th, 2012]
- Genome interpreter vies for place in clinical market [Last Updated On: October 10th, 2012] [Originally Added On: October 10th, 2012]
- The $1,000 Genome: A Bait and Switch? [Last Updated On: October 10th, 2012] [Originally Added On: October 10th, 2012]
- Mount Sinai School of Medicine Offers First-Ever Course with Whole Genome Sequencing [Last Updated On: October 10th, 2012] [Originally Added On: October 10th, 2012]
- First whole genome sequencing of multiple pancreatic cancer patients has been outlined [Last Updated On: October 11th, 2012] [Originally Added On: October 11th, 2012]
- Cheap genome sequences demand new rules on privacy [Last Updated On: October 11th, 2012] [Originally Added On: October 11th, 2012]
- UConn Gets Grant For Genome Research [Last Updated On: October 11th, 2012] [Originally Added On: October 11th, 2012]
- Inconsistent Genome Privacy Laws Need Toughening, Panel Says [Last Updated On: October 12th, 2012] [Originally Added On: October 12th, 2012]
- US panel calls for stronger privacy for genome data [Last Updated On: October 12th, 2012] [Originally Added On: October 12th, 2012]
- Genome Canada Board Appoints New Chair [Last Updated On: October 12th, 2012] [Originally Added On: October 12th, 2012]
- The $1,000 Genome Is Almost Here- Are We Ready? [Last Updated On: October 15th, 2012] [Originally Added On: October 15th, 2012]
- Global genome effort seeks genetic roots of disease [Last Updated On: October 31st, 2012] [Originally Added On: October 31st, 2012]
- Massive encyclopedia helps explain how the human genome works [Last Updated On: October 31st, 2012] [Originally Added On: October 31st, 2012]
- Genome evolution and carbon dioxide dynamics [Last Updated On: October 31st, 2012] [Originally Added On: October 31st, 2012]