Protein backbone flexibility pattern is evolutionarily conserved in the Flaviviridae family: A case of NS3 protease in Flavivirus and Hepacivirus

Piperazine-derived small molecules as potential Flaviviridae NS3 protease inhibitors. In vitro antiviral activity evaluation against Zika and Dengue viruses

Since 2011 Direct Acting antivirals (DAAs) drugs targeting different non-structural (NS) viral proteins (NS3, NS5A or NS5B inhibitors) have been approved for clinical use in HCV therapies. However, currently there are not licensed therapeutics to treat Flavivirus infections and the only licensed DENV vaccine, Dengvaxia, is restricted to patients with preexisting DENV immunity. Similarly to NS5 polymerase, the NS3 catalytic region is evolutionarily conserved among the Flaviviridae family sharing strong structural similarity with other proteases belonging to this family and therefore is an attractive target for the development of pan-flavivirus therapeutics. In this work we present a library of 34 piperazine-derived small molecules as potential Flaviviridae NS3 protease inhibitors. The library was developed through a privileged structures-based design and then biologically screened using a live virus phenotypic assay to determine the half-maximal inhibitor concentration (IC₅₀) of each compound against ZIKV and DENV. Two lead compounds, 42 and 44, with promising broad-spectrum activity against ZIKV (IC₅₀ 6.6 µM and 1.9 µM respectively) and DENV (IC₅₀ 6.7 µM and 1.4 µM respectively) and a good security profile were identified. Besides, molecular docking calculations were performed to provide insights about key interactions with residues in NS3 proteases' active sites.

The role of conserved arginine and proline residues in enterovirus VP1 protein

BACKGROUND

High diversity of VP1 protein among enteroviruses has been a barrier in developing universally effective antiviral drugs. To maintain structure stability during evolution, several residues of VP1 protein of enteroviruses are conserved. Therefore, investigation of highly conserved residues in VP1 protein may provide information for antiviral drug candidates against enteroviruses.

METHODS

To identify highly conserved amino acid sequences of the VP1 in enterovirus genus, the Consurf and CABS-flex 2.0 web software were applied. Through the combination with secondary structure information, we focused on conserved amino acids of VP1 property analysis.

RESULTS

Most conserved residues of VP1 were in the interior and interacted with VP2, VP3 and VP4 capsid proteins. Structure of EV-A71 (PDB code 4AED) showed conserved residues were at hydrophobic pocket and close to the junction between the loop and β-barrel. Interestingly, arginine was the most common conserved residue of VP1. Proline was the second most common conserved residue and was found in the loop and β-barrel intersection areas. VP1 protein flexibility was associated with the secondary structure. Conserved residues of VP1 in β-barrel showed significantly low flexibility.

CONCLUSION

Through large scale sequence analysis, we identified the amino acid distribution and location of conserved residues in VP1. This knowledge can be extrapolated for the Enterovirus genus and may contribute to developing the potential compound as an anti-enteroviral agent.

Targeting the NS2B-NS3 protease of tick-borne encephalitis virus with pan-flaviviral protease inhibitors

Tick-borne encephalitis (TBE) is a severe neurological disorder caused by tick-borne encephalitis virus (TBEV), a member of the Flavivirus genus. Currently, two vaccines are available in Europe against TBEV. However, TBE cases have been rising in Sweden for the past twenty years, and thousands of cases are reported in Europe, emphasizing the need for antiviral treatments against this virus. The NS2B-NS3 protease is essential for flaviviral life cycle and has been studied as a target for the design of inhibitors against several well-known flaviviruses, but not TBEV. In the present study, Compound 86, a known tripeptidic inhibitor of dengue (DENV), West Nile (WNV) and Zika (ZIKV) proteases, was predicted to be active against TBEV protease using a combination of in silico techniques. Further, Compound 86 was found to inhibit recombinant TBEV protease with an IC₅₀ = 0.92 μM in the in vitro enzymatic assay. Additionally, two more peptidic analogues were synthetized and they displayed inhibitory activities against both TBEV and ZIKV proteases. In particular, Compound 104 inhibited ZIKV protease with an IC₅₀ = 0.25 μM. These compounds represent the first reported inhibitors of TBEV protease to date and provides valuable information for the further development of TBEV as well as pan-flavivirus protease inhibitors.

Identification of a C2-symmetric diol based human immunodeficiency virus protease inhibitor targeting Zika virus NS2B-NS3 protease

Zika virus (ZIKV) is an emerging mosquito-borne flavivirus and infection by ZIKV Asian lineage is known to cause fetal brain anomalies and Guillain-Barrés syndrome. The WHO declared ZIKV a global public health emergency in 2016. However, currently neither vaccines nor antiviral prophylaxis/treatments are available. In this study, we report the identification of a C2-symmetric diol-based Human immunodeficiency virus type-1 (HIV) protease inhibitor active against ZIKV NS2B-NS3 protease. The compound, referred to as 9b, was identified by in silico screening of a library of 6265 protease inhibitors. Molecular dynamics (MD) simulation studies revealed that compound 9b formed a stable complex with ZIKV protease. Interaction analysis of compound 9b's binding pose from the cluster analysis of MD simulations trajectories predicted that 9b mostly interacted with ZIKV NS3. Although designed as an aspartyl protease inhibitor, compound 9b was found to inhibit ZIKV serine protease in vitro with IC₅₀ = 143.25 ± 5.45 µM, in line with the in silico results. Additionally, linear interaction energy method (LIE) was used to estimate binding affinities of compounds 9b and 86 (a known panflavivirus peptide hybrid with IC₅₀ = 1.64 ± 0.015 µM against ZIKV protease). The LIE method correctly predicted the binding affinity of compound 86 to be lower than that of 9b, proving to be superior to the molecular docking methods in scoring and ranking compounds. Since most of the reported ZIKV protease inhibitors are positively charged peptide-hybrids, with our without electrophilic warheads, compound 9b represents a less polar and more drug-like non-peptide hit compound useful for further optimization.Communicated by Ramaswamy Sarma.

Molecular docking, dynamics, and pharmacology studies on bexarotene as an agonist of ligand-activated transcription factors, retinoid X receptors

Retinoid X receptors (RXRs) belong to the nuclear receptor superfamily, and upon ligand activation, these receptors control gene transcription via either homodimerization with themselves or heterodimerization with the partner-nuclear receptor. The protective effects of RXRs and RXR agonists have been reported in several neurodegenerative diseases, including in the retina. This study was aimed to prioritize compounds from natural and synthetic origin retinoids as potential RXR agonists by molecular docking and molecular dynamic simulation strategies. The docking studies indicated bexarotene as a lead compound that can activate various RXR receptor isoforms (α, β, and γ) and has a strong binding affinity to the receptor protein than retinoic acid, which is known as a natural endogenous RXR agonist. Dynamic simulation studies confirmed that the hydrogen bonding and hydrophobic interactions were highly stable in all the three isoforms of the RXR-bexarotene complex. To further validate the significance of the RXR receptor in neurons, in vitro pharmacological treatment of neuronal SH-SY5Y cells with bexarotene was performed. In vitro data from SH-SY5Y cells confirmed that bexarotene activated RXR-simulated neurite outgrowth significantly. We conclude that bexarotene could be potentially used as an exogenous activator of RXRs and emerge as a good drug target for several neurodegenerative disorders.

TRAF6 Plays a Proviral Role in Tick-Borne Flavivirus Infection through Interaction with the NS3 Protease

Tick-borne flaviviruses (TBFVs) can cause life-threatening encephalitis and hemorrhagic fever. To identify virus-host interactions that may be exploited as therapeutic targets, we analyzed the TBFV polyprotein in silico for antiviral protein-binding motifs. We obtained two putative tumor necrosis factor receptor-associated factor 6 (TRAF6)-binding motifs (TBMs) within the protease domain of the viral nonstructural 3 (NS3) protein. Here, we show that TBFV NS3 interacted with TRAF6 during infection and that TRAF6 supports TBFV replication. The proviral role of TRAF6 was not seen with mosquito-borne flaviviruses, consistent with the lack of conserved TBMs. Mutation of the second TBM within NS3 disrupted TRAF6 binding, coincident with reduced abundance of mature, autocatalytically derived form of the NS3 protease and significant virus attenuation in vitro. Our studies reveal insights into how flaviviruses exploit innate immunity for the purpose of viral replication and identify a potential target for therapeutic design.

•

Langat virus (LGTV) NS3 protease interacts with TRAF6 during infection

•

Tick-borne, unlike mosquito-borne, flaviviruses use TRAF6 for optimal replication

•

E117A mutation of LGTV NS3 reduces TRAF6 binding and mature protease abundance

•

LGTV with a mutated TRAF6-binding motif is attenuated in vitro