Rationally derived inhibitors of hepatitis C virus (HCV) p7 channel activity reveal prospect for bimodal antiviral therapy.

Since the 1960s, a single class of agent has been licensed targeting virus-encoded ion channels, or 'viroporins', contrasting the success of channel blocking drugs in other areas of medicine. Although resistance arose to these prototypic adamantane inhibitors of the influenza A virus (IAV) M2 proton channel, a growing number of clinically and economically important viruses are now recognised to encode essential viroporins providing potential targets for modern drug discovery. We describe the first rationally designed viroporin inhibitor with a comprehensive structure-activity relationship (SAR). This step-change in understanding not only revealed a second biological function for the p7 viroporin from hepatitis C virus (HCV) during virus entry, but also enabled the synthesis of a labelled tool compound that retained biological activity. Hence, p7 inhibitors (p7i) represent a unique class of HCV antiviral targeting both the spread and establishment of infection, as well as a precedent for future viroporin-targeted drug discovery.

Interaction of antivirals with a heptameric bundle model of the p7 protein of hepatitis C virus.

A series of ligands are known experimentally to affect the infectivity cycle of the hepatitis C virus. The target protein for the ligands is proposed to be p7, a 63 amino acid polytopic channel-forming protein, with possibly two transmembrane domains. Protein p7 is found to assemble into functional oligomers of various sizes, depending on the genotype (GT). Nine ligands are docked to various sites of a computationally derived heptameric bundle of p7 of GT1a. The energy of interaction, here binding energy, is calculated using three different docking programs (Autodock, MOE, LeadIT). Three protein regions are defined to which the ligands are placed, the loop region and the site with the termini as well as the mid-region which is supposed to track poses inside the putative pore. A common feature is that the loop sites and poses either within the pore or at the intermonomer space of the bundle are preferred for all ligands with proposed binding energies smaller than -10 kJ/mol. BIT225, benzamine, amantadine, and NN-DNJ show good overall scoring.

Decoupled side chain and backbone dynamics for proton translocation - M2 of influenza A.

The 97 amino acid bitopic membrane protein M2 of influenza A forms a tetrameric bundle in which two of the monomers are covalently linked via a cysteine bridge. In its tetrameric assembly the protein conducts protons across the viral envelope and within intracellular compartments during the infectivity cycle of the virus. A key residue in the translocation of the protons is His-37 which forms a planar tetrad in the configuration of the bundle accepting and translocating the incoming protons from the N terminal side, exterior of the virus, to the C terminal side, inside the virus. With experimentally available data from NMR spectroscopy of the transmembrane domains of the tetrameric M2 bundle classical MD simulations are conducted with the protein bundle in different protonation stages in respect to His-37. A full correlation analysis (FCA) of the data sets with the His-37 tetrad either in a fully four times unprotonated or protonated state, assumed to mimic high and low pH in vivo, respectively, in both cases reveal asymmetric backbone dynamics. His-37 side chain rotation dynamics is increased at full protonation of the tetrad compared to the dynamics in the fully unprotonated state. The data suggest that proton translocation can be achieved by decoupled side chain or backbone dynamics. Graphical abstract Visualization of the tetrameric bundle of the transmembrane domains of M2 of influenza A after 200 ns of MD simulations (upper left). The four histidine residues 37 are either not protonated as in M20 or fully protonated is in M24+. The asymmetric dynamics of the backbones are shown after a full correlation analysis (FCA) in blue (lower left). The rotational dynamics of the χ2 dihedral angles of the histidines in M20 (upper right) are less than those in M24+ (lower right).

Small molecule ligand docking to genotype specific bundle structures of hepatitis C virus (HCV) p7 protein.

The genome of hepatitis C virus encodes for an essential 63 amino acid polytopic protein p7 of most likely two transmembrane domains (TMDs). The protein is identified to self-assemble thereby rendering lipid membranes permeable to ions. A series of small molecules such as adamantanes, imino sugars and guanidinium compounds are known to interact with p7. A set of 9 of these small molecules is docked against hexameric bundles of genotypes 5a (bundle-5a) and 1b (bundle-1b) using LeadIT. Putative sites for bundle-5a are identified within the pore and at pockets on the outside of the bundle. For bundle-1b preferred sites are found at the site of the loops. Binding energies are in favour of the guanidinium compounds. Rescoring of the identified poses with HYDE reveals a dehydration penalty for the guanidinium compounds, leaving the adamantanes and imino sugar in a better position. Binding energies calculated by HYDE and those by LeadIT indicate that all compounds are moderate binders.

Asymmetric dynamics of ion channel forming proteins - Hepatitis C virus (HCV) p7 bundles.

Protein p7 of hepatitis C virus (HCV) is a short 63 amino acid membrane protein which homo-oligomerises in the lipid membrane to form ion and proton conducting bundles. Two different genotypes (GTs) of p7, 1a and 5a, are used to simulate hexameric bundles of the protein embedded in a fully hydrated lipid bilayer during 400 ns molecular dynamics (MD) simulations. Whilst the bundle of GT 1a is based on a fully computational derived structure, the bundle of GT 5a is based on NMR spectroscopic data. Results of a full correlation analysis (FCA) reveal that albeit structural differences both bundles screen local minima during the simulation. The collective motion of the protein domains is asymmetric. No 'breathing-mode'-like dynamics is observed. The presence of divalent ions, such as Ca-ions affects the dynamics of especially solvent exposed parts of the protein, but leaves the asymmetric domain motion unaffected.

Viral channel forming proteins--How to assemble and depolarize lipid membranes in silico.

Viral channel forming proteins (VCPs) have been discovered in the late 70s and are found in many viruses to date. Usually they are small and have to assemble to form channels which depolarize the lipid membrane of the host cells. Structural information is just about to emerge for just some of them. Thus, computational methods play a pivotal role in generating plausible structures which can be used in the drug development process. In this review the accumulation of structural data is introduced from a historical perspective. Computational performances and their predictive power are reported guided by biological questions such as the assembly, mechanism of function and drug-protein interaction of VCPs. An outlook of how coarse grained simulations can contribute to yet unexplored issues of these proteins is given. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Genotype-specific differences in structural features of hepatitis C virus (HCV) p7 membrane protein.

The 63 amino acid polytopic membrane protein, p7, encoded by hepatitis C virus (HCV) is involved in the modulation of electrochemical gradients across membranes within infected cells. Structural information relating to p7 from multiple genotypes has been generated in silico (e.g. genotype (GT) 1a), as well as obtained from experiments in form of monomeric and hexameric structures (GTs 1b and 5a, respectively). However, sequence diversity and structural differences mean that comparison of their channel gating behaviour has not thus far been simulated. Here, a molecular model of the monomeric GT 1a protein is optimized and assembled into a hexameric bundle for comparison with both the 5a hexamer structure and another hexameric bundle generated using the GT 1b monomer structure. All bundles tend to turn into a compact structure during molecular dynamics (MD) simulations (Gromos96 (ffG45a3)) in hydrated lipid bilayers, as well as when simulated at 'low pH', which may trigger channel opening according to some functional studies. Both GT 1a and 1b channel models are gated via movement of the parallel aligned helices, yet the scenario for the GT 5a protein is more complex, with a short N-terminal helix being involved. However, all bundles display pulsatile dynamics identified by monitoring water dynamics within the pore.