Transmembrane peptide NB of influenza B: a simulation, structure, and conductance study

Natural History of Influenza B Virus-Current Knowledge on Treatment, Resistance and Therapeutic Options

Influenza B virus (IBV) significantly impacts the health and the economy of the global population. WHO global health estimates project 1 billion flu cases annually, with 3 to 5 million resulting in severe disease and 0.3 to 0.5 million influenza-related deaths worldwide. Influenza B virus epidemics result in significant economic losses due to healthcare expenses, reduced workforce productivity, and strain on healthcare systems. Influenza B virus epidemics, such as the 1987-1988 Yamagata lineage outbreak and the 2001-2002 Victoria lineage outbreak, had a significant global impact. IBV's fast mutation and replication rates facilitate rapid adaptation to the environment, enabling the evasion of existing immunity and the development of resistance to virus-targeting treatments. This leads to annual outbreaks and necessitates the development of new vaccination formulations. This review aims to elucidate IBV's evolutionary genomic organization and life cycle and provide an overview of anti-IBV drugs, resistance, treatment options, and prospects for IBV biology, emphasizing challenges in preventing and treating IBV infection.

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.

•

Early references about the discovery of viral channel forming proteins.

•

Latest structural information about the class of proteins.

•

Identification of structural motifs, assembly mechanism of function and drug action using computational methods.

•

Outlook for the use of coarse grained techniques to address assembly and integration into cellular processes.

NB protein does not affect influenza B virus replication in vitro and is not required for replication in or transmission between ferrets

The influenza B virus encodes a unique protein, NB, a membrane protein whose function in the replication cycle is not, as yet, understood. We engineered a recombinant influenza B virus lacking NB expression, with no concomitant difference in expression or activity of viral neuraminidase (NA) protein, an important caveat since NA is encoded on the same segment and initiated from a start codon just 4 nt downstream of NB. Replication of the virus lacking NB was not different to wild-type virus with full-length NB in clonal immortalized or complex primary cell cultures. In the mouse model, virus lacking NB induced slightly lower IFN-α levels in infected lungs, but this did not affect virus titres or weight loss. In ferrets infected with a mixture of viruses that did or did not express NB, there was no fitness advantage for the virus that retained NB. Moreover, virus lacking NB protein was transmitted following respiratory droplet exposure of sentinel animals. These data suggest no role for NB in supporting replication or transmission in vivo in this animal model. The role of NB and the nature of selection to retain it in all natural influenza B viruses remain unclear.

Knowns and unknowns of influenza B viruses

Influenza B viruses (IBVs) circulate annually along with influenza A (IAV) strains during seasonal epidemics. IBV can dominate influenza seasons and cause severe disease, particularly in children and adolescents. Research has revealed interesting aspects of IBV and highlighted the importance of these viruses in clinical settings. Yet, many important questions remain unanswered. In this review, the clinical relevance of IBV is emphasized, unique features in epidemiology, host range and virology are highlighted and gaps in knowledge pinpointed. Multiple aspects of IBV epidemiology, evolution, virology and immunology are discussed. Future research into IBV is needed to understand how we can prevent severe disease in high-risk groups, especially children and elderly.

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.

Genetic characterisation of influenza B viruses detected in Singapore, 2004 to 2009

BACKGROUND

Influenza B viruses are classified into two main lineages: Yamagata-like and Victoria-like, which differ antigenically and phylogenetically. To understand the evolution of influenza B viruses in South East Asia as well as to determine the vaccine efficacy, we genetically characterised gene segments 4, 6 and 8 from non-tissue culture adapted influenza B viruses detected in Singapore from 2004 to 2009.

METHODS

vRNA were extracted from the nasopharyngeal swabs or nasal washes of SAF servicemen displaying febrile and respiratory symptoms, and subjected to PCR assay to test for the presence of influenza B virus. The PCR-positive specimens were next subjected to sequencing of the full gene segments 4 (HA), 6 (NA/NB) and 8 (NS1/NEP). The nucleotide sequences were aligned together with that of other specimens isolated from South East Asia as well as the vaccine strains. Phylogenetic trees of each gene segment were constructed and the amino acid alignments were analysed.

RESULTS

A majority of the Singaporean specimens analysed in this study, from 2004-2009, had gene segment 4 from the Victoria-like lineage and gene segment 6 from Yamagata-like lineage. Some of these specimens had both gene segments from the Yamagata lineage and this resulted in several vaccine mismatches. Gene segment 8 from majority of these specimens clustered separately from both the Yamagata and Victoria strains. The HA protein of most of the Singaporean specimens isolated post 2000 contained a glycosylation site at position 211, which was not dominant prior to 2000. No amino acid substitution conferring drug-resistance was found in either the HA or NA proteins.

CONCLUSIONS

The presence of both lineages co-circulating post 2000, suggests that a trivalent vaccine is not enough to confer immunity to the general public, strongly endorsing the inclusion of both lineages in the vaccine. Several amino acid substitutions were observed, prompting in depth functional analyses.

Characterization of the nucleocytoplasmic shuttle of the matrix protein of influenza B virus

Influenza B virus is an enveloped negative-strand RNA virus that contributes considerably to annual influenza illnesses in human. The matrix protein of influenza B virus (BM1) acts as a cytoplasmic-nuclear shuttling protein during the early and late stages of infection. The mechanism of this intracellular transport of BM1 was revealed through the identification of two leucine-rich CRM1-dependent nuclear export signals (NESs) (3 to 14 amino acids [aa] and 124 to 133 aa), one bipartite nuclear localization signal (NLS) (76 to 94 aa), and two phosphorylation sites (80T and 84S) in BM1. The biological function of the NLS and NES regions were determined through the observation of the intracellular distribution of enhanced green fluorescent protein (EGFP)-tagged signal peptides, and wild-type, NES-mutant, and NLS-mutant EGFP-BM1. Furthermore, the NLS phosphorylation sites 80T and 84S, were found to be required for the nuclear accumulation of EGFP-NLS and for the efficient binding of EGFP-BM1 to human importin-α1. Moreover, all of these regions/sites were required for the generation of viable influenza B virus in a 12-plasmid virus rescue system.

IMPORTANCE

This study expands our understanding of the life cycle of influenza B virus by defining the dynamic mechanism of the nucleocytoplasmic shuttle of BM1 and could provide a scientific basis for the development of antiviral medication.

Viral channel forming proteins - modeling the target

The cellular and subcellular membranes encounter an important playground for the activity of membrane proteins encoded by viruses. Viral membrane proteins, similar to their host companions, can be integral or attached to the membrane. They are involved in directing the cellular and viral reproduction, the fusion and budding processes. This review focuses especially on those integral viral membrane proteins which form channels or pores, the classification to be so, modeling by in silico methods and potential drug candidates. The sequence of an isolate of Vpu from HIV-1 is aligned with host ion channels and a toxin. The focus is on the alignment of the transmembrane domains. The results of the alignment are mapped onto the 3D structures of the respective channels and toxin. The results of the mapping support the idea of a 'channel-pore dualism' for Vpu.

Viral proteins function as ion channels

Viral ion channels are short membrane proteins with 50–120 amino acids and play an important role either in regulating virus replication, such as virus entry, assembly and release or modulating the electrochemical balance in the subcellular compartments of host cells. This review summarizes the recent advances in viral encoded ion channel proteins (or viroporins), including PBCV-1 KcV, influenza M2, HIV-1 Vpu, HCV p7, picornavirus 2B, and coronavirus E and 3a. We focus on their function and mechanisms, and also discuss viral ion channel protein serving as a potential drug target.

Model generation of viral channel forming 2B protein bundles from polio and coxsackie viruses

2B is a 99 amino acid membrane protein encoded by enteroviruses such as polio and coxsackie viruses with two transmembrane domains. The protein is found to make membranes of infected cells permeable. Using a computational approach which positions the models and assesses stability by molecular dynamics (MD) simulations a putative tetrameric bundle model of 2B is generated. The bundles show a pore lining motif of three lysines followed by a serine. The bundle is discussed in terms of different possible orientations of the helices in the membrane and the consequences this has on the in vivo activity of 2B.

Viral channel-forming proteins

Channel-forming proteins are found in a number of viral genomes. In some cases, their role in the viral life cycle is well understood, in some cases it needs still to be elucidated. A common theme is that their mode of action involves a change of electrochemical or proton gradient across the lipid membrane which modulates the viral or cellular activity. Blocking these proteins can be a suitable therapeutic strategy as for some viruses this may be "lethal." Besides the many biological relevant questions still to be answered, there are also many open questions concerning the biophysical side as well as structural information and the mechanism of function on a molecular level. The immanent biophysical issues are addressed and the work in the field is summarized.

Cytopathic mechanisms of HIV-1

The human immunodeficiency virus type 1 (HIV-1) has been intensely investigated since its discovery in 1983 as the cause of acquired immune deficiency syndrome (AIDS). With relatively few proteins made by the virus, it is able to accomplish many tasks, with each protein serving multiple functions. The Envelope glycoprotein, composed of the two noncovalently linked subunits, SU (surface glycoprotein) and TM (transmembrane glycoprotein) is largely responsible for host cell recognition and entry respectively. While the roles of the N-terminal residues of TM is well established as a fusion pore and anchor for Env into cell membranes, the role of the C-terminus of the protein is not well understood and is fiercely debated. This review gathers information on TM in an attempt to shed some light on the functional regions of this protein.

The Alphavirus 6K Protein Activates Endogenous Ionic Conductances when Expressed in Xenopus Oocytes

The Alphavirus Sindbis 6K protein is involved in several functions. It contributes to the processing and membrane insertion of E1 and PE2 viral envelope glycoproteins and to virus budding. It also permeabilizes Escherichia coli and mammalian cells. These viroporin-like properties have been proposed to help virus budding by modifying membrane permeabilities. We expressed Sindbis virus 6K cRNA in Xenopus oocytes to further characterize the effect of 6K on membrane conductances and permeabilization. Although no intrinsic channel properties were seen, cell shrinkage was observed within 24 h. Voltage-clamp experiments showed that 6K upregulated endogenous currents: a hyperpolarization-activated inward current (I (in)) and a calcium-dependent chloride current (I (Cl)). 6K was located at both the plasma and the endoplasmic reticulum membranes. The plasma membrane current upregulation likely results from disruption of the calcium homeostasis of the cell at the endoplasmic reticulum level. Indeed, 6K cRNA expression induced reticular calcium store depletion and capacitative calcium entry activation. By experimental modifications of the incubation medium, we showed that downstream of these events cell shrinkage resulted from a 6K -induced KCl efflux (I (Cl) upregulation leads to chloride efflux, which itself electrically drives potassium efflux), which was responsible for an osmotic water efflux. Our data confirm that 6K specifically triggers a sequential cascade of events that leads to cytoplasmic calcium elevation and cell permeabilization, which likely play a role in the Sindbis virus life cycle.

Transmembrane helices of membrane proteins may flex to satisfy hydrophobic mismatch

A novel mechanism for membrane modulation of transmembrane protein structure, and consequently function, is suggested in which mismatch between the hydrophobic surface of the protein and the hydrophobic interior of the lipid bilayer induces a flexing or bending of a transmembrane segment of the protein. Studies on model hydrophobic transmembrane peptides predict that helices tilt to submerge the hydrophobic surface within the lipid bilayer to satisfy the hydrophobic effect if the helix length exceeds the bilayer width. The hydrophobic surface of transmembrane helix 1 (TM1) of lactose permease, LacY, is accessible to the bilayer, and too long to be accommodated in the hydrophobic portion of a typical lipid bilayer if oriented perpendicular to the membrane surface. Hence, nuclear magnetic resonance (NMR) data and molecular dynamics simulations show that TM1 from LacY may flex as well as tilt to satisfy the hydrophobic mismatch with the bilayer. In an analogous study of the hydrophobic mismatch of TM7 of bovine rhodopsin, similar flexing of the transmembrane segment near the conserved NPxxY sequence is observed. As a control, NMR data on TM5 of lacY, which is much shorter than TM1, show that TM5 is likely to tilt, but not flex, consistent with the close match between the extent of hydrophobic surface of the peptide and the hydrophobic thickness of the bilayer. These data suggest mechanisms by which the lipid bilayer in which the protein is embedded modulates conformation, and thus function, of integral membrane proteins through interactions with the hydrophobic transmembrane helices.

Peptides as transmembrane segments: Decrypting the determinants for helix–helix interactions in membrane proteins

Although the structural analysis of membrane proteins is advancing, an understanding of the basic principles that underlie their folding and assembly remains limited because of the high insolubility intrinsic to these molecules and concomitant challenges in obtaining crystals. Fortunately, from an experimental standpoint, membrane protein folding can be approximated as the rigid-body docking of pre-formed alpha-helical transmembrane segments one with another to form the final functional protein structure. Peptides derived from the sequences of native alpha-helical transmembrane segments and those that mimic their properties are therefore valuable in the experimental evaluation of protein folding within the membrane. Here we present an overview of the progress made in our laboratory and elsewhere in using peptide models toward defining the sequence requirements and forces stabilizing membrane protein folds.

Planar bilayer studies reveal multiple conductance states for synthetic anion transporters

Compounds of the general type R(1)(2)NCOCH(2)OCH(2)CO-(Gly)(3)-Pro-(Gly)(3)-OCH(2)Ph insert in phospholipid bilayers and conduct ions. Different levels of activity were observed when R(1) was either decyl or octadecyl, as judged either by Cl(-) release, detected by ion selective electrodes, or carboxyfluorescein dequenching, detected by fluorescence. Either method reports average behavior for all ionophores over all liposomes. These methods also show that at least two ionophores are involved in the formation of each pore. Planar bilayer experiments reported here confirm pore formation by these compounds but identify more than one conductance state for each. The pseudo-dimer, in which two molecules of the type shown above are covalently linked, shows only two conductance states, of which one is dominant. This state has been characterized by use of a current-voltage plot.

Computational analysis of mutations in the transmembrane region of Vpu from HIV-1

Vpu is an 81 amino acid integral membrane protein encoded by HIV-1. Its alpha-helical transmembrane (TM) domain (residues approximately 6-28) enhances virion release by oligomerizing into bundles and forming ion-conducting channels across the plasma membrane. Its cytoplasmic domain (residues approximately 29-81) is also alpha-helical and binds to the transmembrane protein CD4, inducing its degradation. Mutations within the TM domain have been found to abrogate enhanced particle release from the infected cell (Tiganos et al. Virology (1998) 251 96-107). A series of computational models of monomeric, pentameric and hexameric Vpu(1-31) mutants have been constructed, embedded in fully hydrated lipid bilayers and subjected to a 3 ns molecular dynamics (MD) simulation. None of the mutations has any destabilizing effect on the secondary and tertiary structure. One of the mutants, in which the position of a tryptophan residue within the TM domain is altered, is known not to induce CD4 degradation; an extended kinked model of this mutant has been generated (Vpu(1-52)IVW-k) and during subsequent MD simulations, the bend between the TM and a part of the cytoplasmic domain is found to unwind and a complex salt bridge involving Lys-37 is formed.

An amino-acid substitution in the influenza-B NB protein affects ion-channel gating

The effects of site-directed mutations in NB, a protein encoded by the influenza B virus that has been shown to form cation-selective ion channels at pH 6.0, were studied on ion channel characteristics in artificial lipid bilayers. It was thought that the residues in the hydrophobic region of NB we selected for mutation might be involved in the transport of cations across the channel and that changes in these residues might affect channel properties such as gating and ion-selectivity. Serine residues at positions 20 and 28, threonine at position 24 and cysteine at position 26 were replaced by alanine. We found that the mutation S20A gave channels that did not gate and that remained open most of the time. Proton permeability of NB channels, as detected by fluorescence quenching, was also altered by the mutation S20A: channels were no longer proton-permeable. The other mutations, S28A, T24A and C26A, did not have any detectable effect on the activity or proton permeability of channels formed by NB. The results indicate that serine 20 may have an important role in normal function of NB channels.

Vpu from HIV-1 on an atomic scale: experiments and computer simulations

Vpu is an 81 amino acid protein encoded by HIV-1. Its role is to amplify viral release by two mechanisms: (i) docking to CD4 with the consequence of targeting CD4 for ubiquitine-mediated degradation, and (ii) formation of ion channels to enhance particle release. The intensive research on its in vivo function, combined with structural investigations, makes this viral membrane protein one of the better characterised membrane proteins. The wealth of structural information enables the use of computational methods to elucidate the mechanisms of function on an atomic scale. The discovery of Vpu and the development of structural models in a chronological order is summarised and first efforts on investigating the mechanics are outlined.

The hepatitis C virus p7 protein forms an ion channel that is inhibited by long-alkyl-chain iminosugar derivatives

We show that hepatitis C virus (HCV) p7 protein forms ion channels in black lipid membranes. HCV p7 ion channels are inhibited by long-alkyl-chain iminosugar derivatives, which have antiviral activity against the HCV surrogate bovine viral diarrhea virus. HCV p7 presents a potential target for antiviral therapy.

The NB protein of influenza B virus is not necessary for virus replication in vitro

The NB protein of influenza B virus is thought to function as an ion channel and therefore would be expected to have an essential function in viral replication. Because direct evidence for its absolute requirement in the viral life cycle is lacking, we generated NB knockout viruses by reverse genetics and tested their growth properties both in vitro and in vivo. Mutants not expressing NB replicated as efficiently as the wild-type virus in cell culture, whereas in mice they showed restricted growth compared with findings for the wild-type virus. Thus, the NB protein is not essential for influenza B virus replication in cell culture but promotes efficient growth in mice.

Interactions of the designed antimicrobial peptide MB21 and truncated dermaseptin S3 with lipid bilayers: molecular-dynamics simulations

Molecular-dynamics simulations covering 30 ns of both a natural and a synthetic antimicrobial peptide in the presence of a zwitterionic lipid bilayer were performed. In both simulations, copies of the peptides were placed in an alpha-helical conformation on either side of the bilayer about 10 A (1 A=0.1 nm) from the interface, with either the hydrophobic or the positively charged face of the helix directed toward the bilayer surface. The degree of peptide-lipid interaction was dependent on the starting configuration: surface binding and subsequent penetration of the bilayer was observed for the hydrophobically oriented peptides, while the charge-oriented peptides demonstrated at most partial surface binding. Aromatic residues near the N-termini of the peptides appear to play an important role in driving peptide-lipid interactions. A correlation between the extent of peptide-lipid interactions and helical stability was observed in the simulations. Insertion of the peptides into the bilayer caused a dramatic increase in the lateral area per lipid and decrease in the bilayer thickness, resulting in substantial disordering of the lipid chains. Results from the simulations are consistent with early stages of proposed mechanisms for the lytic activity of antimicrobial peptides. In addition to these 'free' simulations, 25 ns simulations were carried out with the peptides constrained at three different distances relative to the bilayer interface. The constraint forces are in agreement with the extent of peptide-bilayer insertion observed in the free simulations.

Gauging of the PhoE channel by a single freely diffusing proton

In the present study we combined a continuum approximation with a detailed mapping of the electrostatic potential inside an ionic channel to define the most probable trajectory for proton propagation through the channel (propagation along a structure-supported trajectory (PSST)). The conversion of the three-dimensional diffusion space into propagation along a one-dimensional pathway permits reconstruction of an ion motion by a short calculation (a few seconds on a state-of-the-art workstation) rather than a laborious, time-consuming random walk simulations. The experimental system selected for testing the accuracy of this concept was the reversible dissociation of a proton from a single pyranine molecule (8-hydroxypyrene-1,2,3-trisulfonate) bound by electrostatic forces inside the PhoE ionic channel of the Escherichia coli outer membrane. The crystal structure coordinates were used for calculation of the intra-cavity electrostatic potential, and the reconstruction of the observed fluorescence decay curve was carried out using the dielectric constant of the intra-cavity space as an adjustable parameter. The fitting of past experimental observations (Shimoni, E., Y. Tsfadia, E. Nachliel, and M. Gutman. 1993. Biophys. J. 64:472-479) was carried out by a modified version of the Agmon geminate recombination program (Krissinel, E. B., and N. Agmon. 1996. J. Comp. Chem. 17:1085-1098), where the gradient of the electrostatic potential and the entropic terms were calculated by the PSST program. The best-fitted reconstruction of the observed dynamics was attained when the water in the cavity was assigned epsilon </= 55, corroborating the theoretical estimation of Sansom (Breed, J. R., I. D. Kerr, and M. S. P. Sansom. 1996. Biophys. J. 70:1643-1661). The dielectric constant calculated for reversed micelles of comparable size (Cohen, B., D. Huppert, K. M. Solntsev, Y. Tsfadia, E. Nachliel, and M. Gutman. 2002. JACS. 124:7539-7547) allows us to set a margin of epsilon = 50 +/- 5.

Viral ion channels: structure and function

Viral ion channels are short auxiliary membrane proteins with a length of ca. 100 amino acids. They are found in enveloped viruses from influenza A, influenza B and influenza C (Orthomyxoviridae), and the human immunodeficiency virus type 1 (HIV-1, Retroviridae). The channels are called M2 (influenza A), NB (influenza B), CM2 (influenza C) and Vpu (HIV-1). Recently, in Paramecium bursaria chlorella virus (PBCV-1, Phycodnaviridae), a K+ selective ion channel has been discovered. The viral channels form homo oligomers to allow an ion flux and represent miniaturised systems. Proton conductivity of M2 is established; NB, Vpu and the potassium channel from PBC-1 conduct ions; for CM2 ion conductivity is still under proof. This review summarises the current knowledge of these short viral membrane proteins. Their discovery is outlined and experimental evidence for their structure and function is discussed. Studies using computational methods are presented as well as investigations of drug-protein interactions.

Molecular dynamics simulations of wild-type and mutant forms of the Mycobacterium tuberculosis MscL channel

The crystal structure of the Mycobacterium tuberculosis homolog of the bacterial mechanosensitive channel of large conductance (Tb-MscL) provides a unique opportunity to consider mechanosensitive signal transduction at the atomic level. Molecular dynamics simulations of the Tb-MscL channel embedded in an explicit lipid bilayer and of its C-terminal helical bundle alone in aqueous solvent were performed. C-terminal calculations imply that although the helix bundle structure is relatively unstable at physiological pH, it may have been stabilized under low pH conditions such as those used in the crystallization of the channel. Specific mutations to the C-terminal region, which cause a similar conservation of the crystal structure conformation, have also been identified. Full channel simulations were performed for the wild-type channel and two experimentally characterized gain-of-function mutants, V21A and Q51E. The wild-type Tb-MscL trajectory gives insight into regions of relative structural stability and instability in the channel structure. Channel mutations led to observable changes in the trajectories, such as an alteration of intersubunit interactions in the Q51E mutant. In addition, interesting patterns of protein-lipid interactions, such as hydrogen bonding, arose in the simulations. These and other observations from the simulations are relevant to previous and ongoing experimental studies focusing on characterization of the channel.