Characterizing alpha helical properties of Ebola viral proteins as potential targets for inhibition of alpha-helix mediated protein-protein interactions

Computational analysis of perturbations in the post-fusion Dengue virus envelope protein highlights known epitopes and conserved residues in the Zika virus

The dramatic transformation of the Zika virus (ZIKV) from a relatively unknown virus to a pathogen generating global-wide panic has exposed the dearth of detailed knowledge about this virus. Decades of research in the related Dengue virus (DENV), finally culminating in a vaccine registered for use in endemic regions (CYD-TDV) in three countries, provides key insights in developing strategies for tackling ZIKV, which has caused global panic to microcephaly and Guillain-Barre Syndrome. Dengue virus (DENV), a member of the family Flaviviridae, the causal agent of the self-limiting Dengue fever and the potentially fatal hemorrhagic fever/dengue shock syndrome, has been a scourge in tropical countries for many centuries. The recently solved structure of mature ZIKV (PDB ID:5IRE) has provided key insights into the structure of the envelope (E) and membrane (M) proteins, the primary target of neutralizing antibodies. The previously established MEPP methodology compares two conformations of the same protein and identifies residues with significant spatial and electrostatic perturbations. In the current work, MEPP analyzed the pre-and post-fusion DENV type 2 envelope (E) protein, and identified several known epitopes (His317, Tyr299, Glu26, Arg188, etc.) (MEPPitope). These residues are overwhelmingly conserved in ZIKV and all DENV serotypes, and also enumerates residue pairs that undergo significant polarity reversal. Characterization of α-helices in E-proteins show that α1 is not conserved in the sequence space of ZIKV and DENV. Furthermore, perturbation of α1 in the post-fusion DENV structure includes a known epitope Asp215, a residue absent in the pre-fusion α1. A cationic β-sheet in the GAG-binding domain that is stereochemically equivalent in ZIKV and all DENV serotypes is also highlighted due to a residue pair (Arg286-Arg288) that has a significant electrostatic polarity reversal upon fusion. Finally, two highly conserved residues (Thr32 and Thr40), with little emphasis in existing literature, are found to have significant electrostatic perturbation. Thus, a combination of different computational methods enable the rapid and rational detection of critical residues as epitopes in the search for an elusive therapy or vaccine that neutralizes multiple members of the Flaviviridae family. These secondary structures are conserved in the related Dengue virus (DENV), and possibly rationalize isolation techniques particle adsorption on magnetic beads coated with anionic polymers and anionic antiviral agents (viprolaxikine) for DENV. These amphipathic α-helices could enable design of molecules for inhibiting α-helix mediated protein-protein interactions. Finally, comparison of these secondary structures in proteins from related families illuminate subtle changes in the proteins that might render them ineffective to previously successful drugs and vaccines, which are difficult to identify by a simple sequence or structural alignment. Finally, conflicting results about residues that are involved in neutralizing a DENV-E protein by the potent antibody 5J7 (PDB ID:3J6U) are reported.

Computational analysis of perturbations in the post-fusion Dengue virus envelope protein highlights known epitopes and conserved residues in the Zika virus

The dramatic transformation of the Zika virus (ZIKV) from a relatively unknown virus to a pathogen generating global-wide panic has exposed the dearth of detailed knowledge about this virus. Decades of research in the related Dengue virus (DENV), finally culminating in a vaccine registered for use in endemic regions (CYD-TDV) in three countries, provides key insights in developing strategies for tackling ZIKV, which has caused global panic to microcephaly and Guillain-Barre Syndrome. Dengue virus (DENV), a member of the family Flaviviridae, the causal agent of the self-limiting Dengue fever and the potentially fatal hemorrhagic fever/dengue shock syndrome, has been a scourge in tropical countries for many centuries. The recently solved structure of mature ZIKV (PDB ID:5IRE) has provided key insights into the structure of the envelope (E) and membrane (M) proteins, the primary target of neutralizing antibodies. The previously established MEPP methodology compares two conformations of the same protein and identifies residues with significant spatial and electrostatic perturbations. In the current work, MEPP analyzed the pre-and post-fusion DENV type 2 envelope (E) protein, and identified several known epitopes (His317, Tyr299, Glu26, Arg188, etc.) (MEPPitope). These residues are overwhelmingly conserved in ZIKV and all DENV serotypes, and also enumerates residue pairs that undergo significant polarity reversal. Characterization of α-helices in E-proteins show that α1 is not conserved in the sequence space of ZIKV and DENV. Furthermore, perturbation of α1 in the post-fusion DENV structure includes a known epitope Asp215, a residue absent in the pre-fusion α1. A cationic β-sheet in the GAG-binding domain that is stereochemically equivalent in ZIKV and all DENV serotypes is also highlighted due to a residue pair (Arg286-Arg288) that has a significant electrostatic polarity reversal upon fusion. Finally, two highly conserved residues (Thr32 and Thr40), with little emphasis in existing literature, are found to have significant electrostatic perturbation. Thus, a combination of different computational methods enable the rapid and rational detection of critical residues as epitopes in the search for an elusive therapy or vaccine that neutralizes multiple members of the Flaviviridae family. These secondary structures are conserved in the related Dengue virus (DENV), and possibly rationalize isolation techniques particle adsorption on magnetic beads coated with anionic polymers and anionic antiviral agents (viprolaxikine) for DENV. These amphipathic α-helices could enable design of molecules for inhibiting α-helix mediated protein-protein interactions. Finally, comparison of these secondary structures in proteins from related families illuminate subtle changes in the proteins that might render them ineffective to previously successful drugs and vaccines, which are difficult to identify by a simple sequence or structural alignment. Finally, conflicting results about residues that are involved in neutralizing a DENV-E protein by the potent antibody 5J7 (PDB ID:3J6U) are reported.

Both Epistasis and Diversifying Selection Drive the Structural Evolution of the Ebola Virus Glycoprotein Mucin-Like Domain

Throughout the last 3 decades, Ebola virus (EBOV) outbreaks have been confined to isolated areas within Central Africa; however, the 2014 variant reached unprecedented transmission and mortality rates. While the outbreak was still under way, it was reported that the variant leading up to this outbreak evolved faster than previous EBOV variants, but evidence for diversifying selection was undetermined. Here, we test this selection hypothesis and show that while previous EBOV outbreaks were preceded by bursts of diversification, evidence for site-specific diversifying selection during the emergence of the 2014 EBOV clade is weak. However, we show strong evidence supporting an interplay between selection and correlated evolution (epistasis), particularly in the mucin-like domain (MLD) of the EBOV glycoprotein. By reconstructing ancestral structures of the MLD, we further propose a structural mechanism explaining how the substitutions that accumulated between 1918 and 1969 distorted the MLD, while more recent epistatic substitutions restored part of the structure, with the most recent substitution being adaptive. We suggest that it is this complex interplay between weak selection, epistasis, and structural constraints that has shaped the evolution of the 2014 EBOV variant.

IMPORTANCE

The role that selection plays in the emergence of viral epidemics remains debated, particularly in the context of the 2014 EBOV outbreak. Most critically, should such evidence exist, it is generally unclear how this relates to function and increased virulence. Here, we show that the viral lineage leading up to the 2014 outbreak underwent a complex interplay between selection and correlated evolution (epistasis) in a protein region that is critical for immune evasion. We then reconstructed the three-dimensional structure of this domain and showed that the initial mutations in this lineage deformed the structure, while subsequent mutations restored part of the structure. Along this mutational path, the first and last mutations were adaptive, while the intervening ones were epistatic. Altogether, we provide a mechanistic model that explains how selection and epistasis acted on the structural constraints that materialized during the 2014 EBOV outbreak.

The PDB database is a rich source of alpha-helical anti-microbial peptides to combat disease causing pathogens

The therapeutic potential of α-helical anti-microbial peptides (AH-AMP) to combat pathogens is fast gaining prominence. Based on recently published open access software for characterizing α-helical peptides (PAGAL), we elucidate a search methodology (SCALPEL) that leverages the massive structural data pre-existing in the PDB database to obtain AH-AMPs belonging to the host proteome. We provide in vitro validation of SCALPEL on plant pathogens ( Xylella fastidiosa, Xanthomonas arboricola and Liberibacter crescens) by identifying AH-AMPs that mirror the function and properties of cecropin B, a well-studied AH-AMP. The identified peptides include a linear AH-AMP present within the existing structure of phosphoenolpyruvate carboxylase (PPC20), and an AH-AMP mimicing the properties of the two α-helices of cecropin B from chitinase (CHITI25). The minimum inhibitory concentration of these peptides are comparable to that of cecropin B, while anionic peptides used as control failed to show any inhibitory effect on these pathogens. Substitute therapies in place of conventional chemotherapies using membrane permeabilizing peptides like these might also prove effective to target cancer cells. The use of native structures from the same organism largely ensures that administration of such peptides will be better tolerated and not elicit an adverse immune response. We suggest a similar approach to target Ebola epitopes, enumerated using PAGAL recently, by selecting suitable peptides from the human proteome, especially in wake of recent reports of cationic amphiphiles inhibiting virus entry and infection.

The PDB database is a rich source of alpha-helical anti-microbial peptides to combat disease causing pathogens

The therapeutic potential of α-helical anti-microbial peptides (AH-AMP) to combat pathogens is fast gaining prominence. Based on recently published open access software for characterizing α-helical peptides (PAGAL), we elucidate a search methodology (SCALPEL) that leverages the massive structural data pre-existing in the PDB database to obtain AH-AMPs belonging to the host proteome. We provide in vitro validation of SCALPEL on plant pathogens ( Xylella fastidiosa, Xanthomonas arboricola and Liberibacter crescens) by identifying AH-AMPs that mirror the function and properties of cecropin B, a well-studied AH-AMP. The identified peptides include a linear AH-AMP present within the existing structure of phosphoenolpyruvate carboxylase (PPC20), and an AH-AMP mimicing the properties of the two α-helices of cecropin B from chitinase (CHITI25). The minimum inhibitory concentration of these peptides are comparable to that of cecropin B, while anionic peptides used as control failed to show any inhibitory effect on these pathogens. Substitute therapies in place of conventional chemotherapies using membrane permeabilizing peptides like these might also prove effective to target cancer cells. The use of native structures from the same organism could possibly ensure that administration of such peptides will be better tolerated and not elicit an adverse immune response. We suggest a similar approach to target Ebola epitopes, enumerated using PAGAL recently, by selecting suitable peptides from the human proteome, especially in wake of recent reports of cationic amphiphiles inhibiting virus entry and infection.

A common feature pharmacophore for FDA-approved drugs inhibiting the Ebola virus

We are currently faced with a global infectious disease crisis which has been anticipated for decades. While many promising biotherapeutics are being tested, the search for a small molecule has yet to deliver an approved drug or therapeutic for the Ebola or similar filoviruses that cause haemorrhagic fever. Two recent high throughput screens published in 2013 did however identify several hits that progressed to animal studies that are FDA approved drugs used for other indications. The current computational analysis uses these molecules from two different structural classes to construct a common features pharmacophore. This ligand-based pharmacophore implicates a possible common target or mechanism that could be further explored. A recent structure based design project yielded nine co-crystal structures of pyrrolidinone inhibitors bound to the viral protein 35 (VP35). When receptor-ligand pharmacophores based on the analogs of these molecules and the protein structures were constructed, the molecular features partially overlapped with the common features of solely ligand-based pharmacophore models based on FDA approved drugs. These previously identified FDA approved drugs with activity against Ebola were therefore docked into this protein. The antimalarials chloroquine and amodiaquine docked favorably in VP35. We propose that these drugs identified to date as inhibitors of the Ebola virus may be targeting VP35. These computational models may provide preliminary insights into the molecular features that are responsible for their activity against Ebola virus in vitro and in vivo and we propose that this hypothesis could be readily tested.

A common feature pharmacophore for FDA-approved drugs inhibiting the Ebola virus

We are currently faced with a global infectious disease crisis which has been anticipated for decades. While many promising biotherapeutics are being tested, the search for a small molecule has yet to deliver an approved drug or therapeutic for the Ebola or similar filoviruses that cause haemorrhagic fever. Two recent high throughput screens published in 2013 did however identify several hits that progressed to animal studies that are FDA approved drugs used for other indications. The current computational analysis uses these molecules from two different structural classes to construct a common features pharmacophore. This ligand-based pharmacophore implicates a possible common target or mechanism that could be further explored. A recent structure based design project yielded nine co-crystal structures of pyrrolidinone inhibitors bound to the viral protein 35 (VP35). When receptor-ligand pharmacophores based on the analogs of these molecules and the protein structures were constructed, the molecular features partially overlapped with the common features of solely ligand-based pharmacophore models based on FDA approved drugs. These previously identified FDA approved drugs with activity against Ebola were therefore docked into this protein. The antimalarials chloroquine and amodiaquine docked favorably in VP35. We propose that these drugs identified to date as inhibitors of the Ebola virus may be targeting VP35. These computational models may provide preliminary insights into the molecular features that are responsible for their activity against Ebola virus in vitro and in vivo and we propose that this hypothesis could be readily tested.