Structure and inhibition of EV-D68, a virus that causes respiratory illness in children

See allHide authors and affiliations

Science  02 Jan 2015:
Vol. 347, Issue 6217, pp. 71-74
DOI: 10.1126/science.1261962

Targeting EV-D68, a respiratory virus

A recent outbreak of respiratory illness in U.S. children was caused by entorovirus D68 (EV-D68). Enteroviruses also include human pathogens such as human rhinovirus, which causes the common cold, and poliovirus. Most of these viruses are stabilized by a factor that binds in a hydrophobic pocket of the capsid protein VP1, and antiviral compounds can act by displacing this factor. Liu et al. report the crystal structure of EV-D68 and its complex with the antiviral compound peconaril. In EV-D68, the hydrophobic pocket contained a fatty acid that was displaced by peconaril. Peconaril efficiently inhibited EV-D68 infection of cells, making it a possible drug candidate against EV-D68.

Science, this issue p. 71


Enterovirus D68 (EV-D68) is a member of Picornaviridae and is a causative agent of recent outbreaks of respiratory illness in children in the United States. We report here the crystal structures of EV-D68 and its complex with pleconaril, a capsid-binding compound that had been developed as an anti-rhinovirus drug. The hydrophobic drug-binding pocket in viral protein 1 contained density that is consistent with a fatty acid of about 10 carbon atoms. This density could be displaced by pleconaril. We also showed that pleconaril inhibits EV-D68 at a half-maximal effective concentration of 430 nanomolar and might, therefore, be a possible drug candidate to alleviate EV-D68 outbreaks.

Picornaviruses constitute a large family of small icosahedral viruses with a single positive-stranded RNA genome and an external diameter of about 300 Å. The Enterovirus genus includes medically important human pathogens, such as human rhinoviruses (HRVs), polioviruses (PVs) and coxsackieviruses (CVs) (table S1) (1, 2). Many of these enteroviruses (EVs) have been well characterized structurally and functionally (310). However, the species EV-D remains poorly characterized.

An upsurge of EV-D68 cases in the past few years has shown clusters of infections worldwide (11). In August 2014, an outbreak of mild to severe respiratory illnesses occurred among thousands of young children in the United States, of which 1116 cases have been confirmed to be caused by EV-D68. This virus has also been associated with occasional neurological infections (12). Although EV-D68 has emerged as a considerable global public health threat, there is no available vaccine or effective antiviral treatment.

The capsids of EVs consist of 60 copies of each of four different viral proteins: VP1, VP2, VP3, and VP4 (Fig. 1A). Of these, VP1 (about 300 amino acids), VP2 (about 260 amino acids), and VP3 (about 240 amino acids) each have a “jelly roll” fold arranged in the capsid with pseudo T = 3 icosahedral symmetry, where T represents the triangulation number (13). Their organization in the capsid is similar to that of the T = 3 RNA plant viruses, except that the three subunits related by quasi-threefold symmetry have different amino acid sequences in picornaviruses (6, 10). Each roughly 70–amino–acids-long VP4 molecule forms an extended peptide on the internal surface of the capsid shell. The jelly roll fold consists of two antiparallel β sheets, which face each other to form a β barrel with a hydrophobic interior (Fig. 1B).

Fig. 1 The structure of EV-D68.

(A) Diagrammatic representation of the virus. VP1, VP2, and VP3 are colored blue, green, and red, respectively. Each pentagonal icosahedral vertex is surrounded by a canyon colored black. The epitopes for neutralizing antibodies for the homologous HRV14 are labeled as NIm-I, NIm-II, and NIm-III. The NIm-IA and NIm-IB sites are disordered in EV-D68. (B) The VP1 jelly roll in EV-D68 shown as a ribbon diagram. If the β strands along the polypeptide are identified sequentially as A, B, C,….I, then one of the sheets is composed of the antiparallel strands BIDG and the other by the antiparallel strands CHEF. (C) Surface features of EV-D68 compared to HRV16. The canyon of EV-D68 is shallower and narrower than the canyon of HRV16.

EVs have a deep surface depression (“canyon”) circulating around each of the 12 pentameric vertices (Fig. 1A). The canyon was predicted to be the site of receptor binding, because the amino acids outside the canyon that form the external surface of the virus were more exposed and were shown to be accessible to neutralizing antibodies (Fig. 1A) (10, 14). The virus thus could remain faithful to a specific receptor molecule that binds into the canyon while evading the host’s immune system (10, 15). The prediction that the canyon would be the site of binding to cellular receptor molecules was subsequently confirmed for numerous different EVs (16, 17). All of the receptor molecules were found to have an immunoglobulin-like structure.

A variety of small ~350-dalton hydrophobic molecules that are inhibitors of EV infection were shown to bind into the hydrophobic pocket in the center of the VP1 jelly roll (18) (Fig. 1B). They stabilize the virus (19) by filling the pocket with a well-fitting hydrophobic molecule, thereby inhibiting uncoating of the virus and the release of the genome into the infected cell. These compounds also prevent viral attachment to cells by altering the surface features of the canyon floor, where the virus attaches to a cellular receptor (17, 19). The floor of the canyon is formed primarily by the GH loop of VP1, the connecting residues between β-strands G and H (18).

Most infectious EVs contain a small molecule or “pocket factor,” probably a fatty acid in the VP1 binding pocket (7, 18, 20). Like the capsid-binding compounds, the pocket factor presumably stabilizes the virus by filling the VP1 hydrophobic pocket. Thus, the virus is stabilized while being transmitted to a new host. However, when a receptor molecule binds to the floor of the canyon, it depresses the floor (that also forms the roof of the VP1 binding pocket), which then squeezes the binding pocket, probably expelling the pocket factor (17). Thus, the attachment of the virus to a cell surface initiates uncoating, causing the release of the viral genome into the cell’s cytosol.

In most EVs that have been investigated, these capsid-binding antiviral compounds displace the pocket factor (e,g. polioviruses, coxsackieviruses A or B, and many HRVs). The rhinoviruses HRV14 (18) and HRV3 (21) do not contain a pocket factor, and the GH loop is displaced relative to its position in other EVs that contain a pocket factor, reducing the volume of the VP1 pocket. Either neither HRV14 nor HRV3 binds a pocket factor in vivo or it was lost during the purification procedure. Much effort was made between 1985 and 2000 to design a compound that fits into the VP1 pocket and would inhibit the maximum number of rhinovirus serotypes (22). The final optimal structure was pleconaril (fig. S1), which not only had good efficacy but was also stable enough to maintain good bioavailability in clinical tests (23). However, pleconaril was not licensed, primarily because it put women using birth control drugs at risk of conception.

The EV-D68 prototype strain Fermon CA62-1 was propagated in human rhabdomyosarcoma cells at 33°C, which had previously been shown to be the optimal growth temperature for EV-D68, indicating that EV-D68 behaves much like members of the RV-A and RV-B species that are responsible for common colds (24). The virus was purified and crystalized as described in the supplementary materials and methods. The cubic-looking crystals had a diameter between 0.1 and 0.2 mm. After soaking in glycerol, they were flash-frozen in liquid nitrogen. X-ray diffraction data were collected at sector 14 of the Advanced Photon Source. The data extended to 2.0 Å resolution and were processed with HKL2000 (25). The crystal symmetry was I222, with two particles per unit cell, implying that the particles were located on a 222 symmetry position. A rotation function (26) differentiated between the two orthogonal possible orientations of the icosahedron. Initial phases were calculated based on the structure of HRV2 (27) (Protein Data Bank accession number 1FPN) after the pocket factor was removed. The phases were then extended in small steps to 2.0 Å resolution, using 15-fold averaging and solvent flattening. A model of the structure was built using the program Coot (28) and refined with the program CNS (29). The final value of Rwork was 27.5%. Rfree is essentially the same as Rwork in the presence of high noncrystallographic symmetry redundancy as is the case here (table S2).

Comparison of the amino acid sequences of EV-D68 with those of other EVs shows that VP3 has a short C-terminal α helix not present in other EVs (fig. S2). The EV-D68 electron density map showed that this helix decorates the north side of the canyon in the neighboring, fivefold-related, icosahedral asymmetric unit. As a result, the canyon is narrower than in other EVs and might therefore not be large enough to accommodate immunoglobulinlike receptors (Fig. 1C).

The BC and DE loops of VP1 are structurally the most variable among known picornaviruses. The EV-D68 VP1 has two disordered regions corresponding to residues 80 to 86 (EV-D68 numbering) in the BC loop and 129 to 136 in the DE loop, both of which are near fivefold axes. These regions harbor the neutralizing immunogen sites NIm-IA and NIm-IB on HRV14, respectively (10) (Fig. 1A). Thus the flexible immunogenic regions around the fivefold axes might be an alternative mechanism for evading host humoral immune responses.

The electron density map of EV-D68 showed density in the VP1 pocket (Fig. 2A). The height of this pocket factor density was about three standard deviations above the mean of the noncrystallographic symmetry averaged map, as compared with about five standard deviations of most of the main-chain density. As in other EVs that have a pocket factor, the conformation of the GH loop of VP1, which defines the interface between the VP1 pocket and the floor of the canyon, is pushed into the canyon relative to the empty pocket in HRV14 (10, 18) or HRV3 (21).

Fig. 2

The pocket factor. (A) One icosahedral asymmetric unit of the EV-D68 structure, showing the Cα atom backbone for VP1, VP2, and VP3 in marine blue, green and red, respectively. The pocket factor electron density outline is shown in gray. (B) Enlargement of the pocket factor density with a fitted putative C10 fatty acid. Shown also are the amino acids lining the pocket. (C) Comparison of the putative pocket factor structures in six known EV structures. Shown also is the residue I184 in EV-D68 that alters the direction of the shorter pocket factors in some of the viruses as compared to others with longer pocket factors and a smaller residue at the equivalent position of 184. Atoms at the head of the pocket factors and of the VP1 residue side chains are colored red (oxygen), dark blue (nitrogen), and yellow (sulfur).

The length of the pocket factor density in the EV-D68 map corresponded to a fatty acid with an aliphatic chain of about 10 carbon atoms (Fig. 2B). Similarly, 12–carbon-atom-long pocket factors were observed for HRV16 (9) and HRV2 (27). In contrast, the well-formed pocket factors in poliovirus 1 (6, 7), coxsackievirus B3 (5), and EV-A71 (3) corresponded to longer fatty acids with 18, 16, and 18 carbon atoms, respectively. Furthermore, the orientation of the pocket factor tails in EV-D68, HRV16, and HRV2 is slightly different than that of the long pocket factors of EV-A71, poliovirus 1, and coxsackievirus B3 (Fig. 2C). The orientation is controlled in part by VP1 residue 184, which is Leu or Ile for the viruses with small, short pocket factors and Val for viruses with longer pocket factors (fig. S2). The larger Leu or Ile residue pushes the pocket factor sideways. The similarity between the properties of the pocket factor in EV-D68 and that in the HRVs might partially explain why these viruses are less stable.

The anti–EV-D68 activity of two capsid-binding compounds, pirodavir and BTA-188 (fig. S1), that had significant anti-rhinovirus activity was compared with pleconaril, using plaque reduction assays in HeLa cells (Fig. 3, A and B). The half-maximal effective concentration (EC50) value of these two compounds was found to be comparable to previous results using cytopathic effect inhibition assays against EV-D68 (30). However, pleconaril was found to be more potent against EV-D68 than pirodavir and BTA-188 (Fig. 3A). The inhibitory effect of pleconaril is similar against EV-D68, HRV16, and HRV14, but better than against EV-A71. It is therefore noteworthy that pleconaril was an effective inhibitor in extensive clinical tests for the treatment of common colds (23, 31). Furthermore, fluorescence-based thermal stability assays indicated that when EV-D68 was incubated with pleconaril at either 10 or 50 μg/ml, 4°C higher temperatures were required to release the RNA genome than when no pleconaril was present (Fig. 3C and materials and methods in the supplementary materials). Thus, pleconaril stabilizes EV-D68 capsids, preventing the virus from uncoating during viral entry.

Fig. 3

Anti-EV activity of three capsid-binding compounds. (A) List of EC50 values (30, 32, 33). Each experiment was performed at least three times. (B) Plot of the percent of plaque reduction of pleconaril as a function of the log of its concentration. Error bars indicate standard deviations. (C) Release of EV-D68 genome upon increase in temperature as monitored by Sybr green II. The experiments were done in triplicate. For each data point, the ratio between mean fluorescence intensity and standard deviation is at least 7.6. Shown are curves fitted with a sigmoidal function for the native virus and for the virus after incubation with pleconaril.

The structure of EV-D68 was also determined to 2.3 Å resolution when co-crystalized with pleconaril (fig. S3). The crystallographic procedure was the same as for the native structure determination (materials and methods in the supplementary materials). The electron density of the ligand inside the hydrophobic pocket was of the same height but much longer than that of the pocket factor in the native EV-D68 structure (Fig. 4A) and could be easily fitted with the structure of pleconaril (Fig. 4B). This demonstrated that the native pocket factor was replaced by pleconaril. However, the densities of the ligand and of the pocket factor are slightly lower than that of the polypeptide main chain. In part this is to be expected, because the ligand or fatty acid pocket factor is composed primarily of only carbon atoms as opposed to the heavier combination of carbon and nitrogen plus oxygen atoms of the main chain. In part, the slightly lower density of the ligand might be due to some variations of conformation within the VP1 binding pocket.

Fig. 4 Structure of pleconaril bound into the VP1 pocket of EV-D68.

(A) Pocket factor density (gray) compared to the pleconaril density (magenta). (B) Pleconaril (green) fitted to density in the structure of the complex. (C) Conformational change of the VP1 GH loop as a consequence of the presence of pleconaril. The native and complex structures are shown in marine blue and baby blue, respectively. Oxygen, nitrogen, sulfur, and fluorine atoms are shown in red, dark blue, yellow, and light green, respectively.

Part of the VP1 GH loop (residues 212 to 215) that forms the entrance to the VP1 pocket had become less ordered in the EV-D68-pleconaril complex. The Cα atom of residue 211 had moved 1.2 Å toward the inside of the pocket relative to that of the native structure, possibly blocking the entrance to the pocket once pleconaril had entered (Fig. 4C). Thus, the dynamics of the GH loop might be a consideration for future structure-based design of EV-D68 capsid-binding inhibitors.

A comparison of the EV-D68–pleconaril, HRV14-pleconaril, and HRV16-pleconaril structures showed a similar binding mode for pleconaril in the VP1 pocket of these three viruses (fig. S4 and table S3). This may explain why pleconaril is similarly effective against these three EVs. To investigate why pleconaril is more effective against EV-D68 than pirodavir or BTA-188, we performed in silico docking experiments. The presence of a trifluoromethyl-substituted oxadiazole moiety in pleconaril, rather than a more hydrophilic group in either pirodavir (ethyl carboxylate group) or BTA-188 (O-ethyloxime group) at structurally equivalent positions (fig. S5), probably contributes to more favorable interactions of pleconaril with the hydrophobic residues deep inside the VP1 pocket of EV-D68.

Our results show that the structure of EV-D68 has considerable similarities to those of the well-studied HRVs for which pleconaril was specifically designed. We also show that pleconaril replaces the pocket factor and is a potent inhibitor of EV-D68, with an EC50 value of 430 nM. The size and location of the pocket factor lodged in the VP1 pocket are similar to those found in other HRVs and different from those of the pocket factors found in poliovirus 1 and EV-A71. This correlates with the observation that pleconaril is far more active when the natural pocket factor is short, as in the HRVs and in EV-D68. Furthermore, sequence alignment of 188 EV-D68 strains found between 1962 and 2013 indicates that residues in VP1 that interact with pleconaril, as identified from the complex structure, are completely conserved, with one exception. Therefore, pleconaril is likely to inhibit not only the prototype strain examined here but also many other strains. In view of the previous extensive clinical trials that have established its safety, pleconaril would be a possible drug candidate to alleviate EV-D68 outbreaks.

Supplementary Materials

Materials and Methods

Figs. S1 to S5

Tables S1 to S3

References (3449)

References and Notes

  1. Acknowledgments: We thank M. Steven Oberste of the Centers for Disease Control and Prevention for supplying us with the prototype strain of EV-D68; M. A. McKinlay of the Task Force for Global Health for helpful discussion and suggestions; S. Kelly for help with the manuscript preparation; and V. Srajer, R. Henning, and the other staff of the Advanced Photon Source BioCARS beamline 14 for help with x-ray diffraction data collection. Use of BioCARS sector 14 was supported by the National Institutes of Health, National Center for Research Resources (NIH/NCRR) grant RR007707. Use of the Advanced Photon Source was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under contract DE-AC02-006CH11357. This study was supported by NIH grant award AI11219 to M.G.R. Coordinates for native EV-D68 and EV-D68–pleconaril structures were deposited with the Protein Data Bank with accession numbers 4WM8 and 4WM7, respectively.
View Abstract

Stay Connected to Science

Navigate This Article