Ion-dynamics in hepatitis C virus p7 helical transmembrane domains — a molecular dynamics simulation study

Hepatitis C virus sequence divergence preserves p7 viroporin structural and dynamic features

The hepatitis C virus (HCV) viroporin p7 oligomerizes to form ion channels, which are required for the assembly and secretion of infectious viruses. The 63-amino acid p7 monomer has two putative transmembrane domains connected by a cytosolic loop, and has both N- and C- termini exposed to the endoplasmic reticulum (ER) lumen. NMR studies have indicated differences between p7 structures of distantly related HCV genotypes. A critical question is whether these differences arise from the high sequence variation between the different isolates and if so, how the divergent structures can support similar biological functions. Here, we present a side-by-side characterization of p7 derived from genotype 1b (isolate J4) in the detergent 6-cyclohexyl-1-hexylphosphocholine (Cyclofos-6) and p7 derived from genotype 5a (isolate EUH1480) in n-dodecylphosphocholine (DPC). The 5a isolate p7 in conditions previously associated with a disputed oligomeric form exhibits secondary structure, dynamics, and solvent accessibility broadly like those of the monomeric 1b isolate p7. The largest differences occur at the start of the second transmembrane domain, which is destabilized in the 5a isolate. The results show a broad consensus among the p7 variants that have been studied under a range of different conditions and indicate that distantly related HCVs preserve key features of structure and dynamics.

All-atom virus simulations

The constant threat of viral disease can be combated by the development of novel vaccines and therapeutics designed to disrupt key features of virus structure or infection cycle processes. Such development relies on high-resolution characterization of viruses and their dynamical behaviors, which are often challenging to obtain solely by experiment. In response, all-atom molecular dynamics simulations are widely leveraged to study the structural components of viruses, leading to some of the largest simulation endeavors undertaken to date. The present work reviews exemplary all-atom simulation work on viruses, as well as progress toward simulating entire virions.

Perturbations of Native Membrane Protein Structure in Alkyl Phosphocholine Detergents: A Critical Assessment of NMR and Biophysical Studies

Membrane proteins perform a host of vital cellular functions. Deciphering the molecular mechanisms whereby they fulfill these functions requires detailed biophysical and structural investigations. Detergents have proven pivotal to extract the protein from its native surroundings. Yet, they provide a milieu that departs significantly from that of the biological membrane, to the extent that the structure, the dynamics, and the interactions of membrane proteins in detergents may considerably vary, as compared to the native environment. Understanding the impact of detergents on membrane proteins is, therefore, crucial to assess the biological relevance of results obtained in detergents. Here, we review the strengths and weaknesses of alkyl phosphocholines (or foscholines), the most widely used detergent in solution-NMR studies of membrane proteins. While this class of detergents is often successful for membrane protein solubilization, a growing list of examples points to destabilizing and denaturing properties, in particular for α-helical membrane proteins. Our comprehensive analysis stresses the importance of stringent controls when working with this class of detergents and when analyzing the structure and dynamics of membrane proteins in alkyl phosphocholine detergents.

The Unusual Transmembrane Partition of the Hexameric Channel of the Hepatitis C Virus

The p7 protein of the hepatitis C virus (HCV) can oligomerize in membrane to form cation channels. Previous studies showed that the channel assembly in detergent micelles adopts a unique flower-shaped oligomer, but the unusual architecture also presented problems for understanding how this viroporin resides in the membrane. Moreover, the oligomeric state of p7 remains controversial, as both hexamer and heptamer have been proposed. Here we address the above issues using p7 reconstituted in bicelles that mimic a lipid bilayer. We found, using a recently developed oligomer-labeling method, that p7 forms hexamers in the bicelles. Solvent paramagnetic relaxation enhancement analyses showed that the bilayer thickness around the HCV ion channel is substantially smaller than expected, and thus a significant portion of the previously assigned membrane-embedded region is solvent exposed. Our study provides an effective approach for characterizing the transmembrane partition of small ion channels in near lipid bilayer environment.

Weak Selectivity Predicted for Modeled Bundles of Viral Channel-Forming Protein E5 of Human Papillomavirus-16

Protein E5 is a polytopic 83 amino acid membrane protein with three transmembrane domains (TMDs), encoded by high-risk human papillomavirus-16 (HPV-16). HPV-16 is found to be the causative agent for cervical cancer. Protein E5, among other proteins (e.g., E6, E7), is expressed at an "early" (E) stage when the cell turns malignant. It has been experimentally found that E5 forms hexameric assemblies, which show the characteristics of the class of so-called channel-forming proteins by rendering lipid membranes permeable to ions and small molecules. Protein E5 is used to achieve structural models of the protein in assembled bundles using a force field-based docking approach. Extended molecular dynamics simulations of selected bundles in fully hydrated lipid bilayers suggest the second TMD to be pore-lining, allowing for water columns to exist within the lumen of the pore. Full correlation analysis indicates asymmetric dynamics within the monomers of the bundle. Potential of mean force calculations of a snapshot structure of the putative open pore of the protein bundle propose low selectivity.

Assessing the physiological relevance of alternate architectures of the p7 protein of hepatitis C virus in different environments

The viroporin p7 of the hepatitis C virus forms multimeric channels eligible for ion transport across the endoplasmic reticulum membrane. Currently the subject of many studies and discussion, the molecular assembly of the ion channel and the structural characteristics of the p7 monomer are not yet fully understood. Structural investigation of p7 has been carried out only in detergent environments, making the interpretation of the experimental results somewhat questionable. Here, we analyze by means of molecular dynamics simulations the structure of the p7 monomer as a function of its sequence, initial conformation and environment. We investigate the conductance properties of three models of a hexameric p7 ion channel by examining ion translocation in a pure lipid bilayer. It is noteworthy that although none of the models reflects the experimentally observed trend to conduct preferentially cations, we were able to identify the position and orientation of titratable acidic or basic residues playing a crucial role in ion selectivity and in the overall conductance of the channel. In addition, too compact a packing of the monomers leads to channel collapse rather than formation of a reasonable pore, amenable to ion translocation. The present findings are envisioned to help assess the physiological relevance of p7 ion channel models consisting of multimeric structures obtained in non-native environments.

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.

Investigation of the roles of trace elements during hepatitis C virus infection using protein-protein interactions and a shortest path algorithm

BACKGROUND

Hepatitis is a type of infectious disease that induces inflammation of the liver without pinpointing a particular pathogen or pathogenesis. Type C hepatitis, as a type of hepatitis, has been reported to induce cirrhosis and hepatocellular carcinoma within a very short amount of time. It is a great threat to human health. Some studies have revealed that trace elements are associated with infection with and immune rejection against hepatitis C virus (HCV). However, the mechanism underlying this phenomenon is still unclear.

METHODS

In this study, we aimed to expand our knowledge of this phenomenon by designing a computational method to identify genes that may be related to both HCV and trace element metabolic processes. The searching procedure included three stages. First, a shortest path algorithm was applied to a large network, constructed by protein-protein interactions, to identify potential genes of interest. Second, a permutation test was executed to exclude false discoveries. Finally, some rules based on the betweenness and associations between candidate genes and HCV and trace elements were built to select core genes among the remaining genes.

RESULTS

12 lists of genes, corresponding to 12 types of trace elements, were obtained. These genes are deemed to be associated with HCV infection and trace elements metabolism.

CONCLUSIONS

The analyses indicate that some genes may be related to both HCV and trace element metabolic processes, further confirming the associations between HCV and trace elements. The method was further tested on another set of HCV genes, the results indicate that this method is quite robustness.

GENERAL SIGNIFICANCE

The newly found genes may partially reveal unknown mechanisms between HCV infection and trace element metabolism. This article is part of a Special Issue entitled "System Genetics" Guest Editor: Dr. Yudong Cai and Dr. Tao Huang.

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.

•

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.

The N-terminal Helical Region of the Hepatitis C Virus p7 Ion Channel Protein Is Critical for Infectious Virus Production

The hepatitis C virus (HCV) p7 protein is required for infectious virus production via its role in assembly and ion channel activity. Although NMR structures of p7 have been reported, the location of secondary structural elements and orientation of the p7 transmembrane domains differ among models. Furthermore, the p7 structure-function relationship remains unclear. Here, extensive mutagenesis, coupled with infectious virus production phenotyping and molecular modeling, demonstrates that the N-terminal helical region plays a previously underappreciated yet critical functional role, especially with respect to E2/p7 cleavage efficiency. Interrogation of specific N-terminal helix residues identified as having p7-specific defects and predicted to point toward the channel pore, in a context of independent E2/p7 cleavage, further supports p7 as a structurally plastic, minimalist ion channel. Together, our findings indicate that the p7 N-terminal helical region is critical for E2/p7 processing, protein-protein interactions, ion channel activity, and infectious HCV production.

Hepatitis C virus (HCV) infection can lead to significant liver disease and, without a vaccine, continues to pose a significant public health threat. The viral p7 protein is a multifunctional protein that is required for infectious virus production via its role in orchestrating virion assembly and its activity as an ion channel. However, while there is accumulating structural information on p7, there is no consensus on which conformation(s) exist during a natural infection or how structural elements relate to p7 functions. By comparing two prominent, yet highly divergent models of p7, we identified one region of structural similarity–the N-terminal helical region. While mutagenesis screening of other regions of the protein are in keeping with p7 conformational flexibility, mutations within the N-terminal helical region had a significant impact on infectious virus production, due in part to a loss of efficient E2/p7 cleavage. We further postulated the precise functional impact of mutations throughout p7 by homology modeling and demonstrated tolerance for diverse amino acid substitutions for specific N-terminal helix residues with putative ion channel defects. Together, these data not only support p7 as a structurally plastic, minimalistic ion channel, but also provide extensive insight into the p7 structure-function relationship and highlight the importance of the N-terminal helical region in E2/p7 processing, protein-protein interactions, ion channel activity, and infectious HCV production.

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.

Structure based computational assessment of channel properties of assembled ORF-8a from SARS-CoV

ORF 8a is a short 39 amino acid bitopic membrane protein encoded by severe acute respiratory syndrome causing corona virus (SARS‐CoV). It has been identified to increase permeability of the lipid membrane for cations. Permeability is suggested to occur due to the assembly of helical bundles. Computational models of a pentameric assembly of 8a peptides are generated using the first 22 amino acids, which include the transmembrane domain. Low energy structures reveal a hydrophilic pore mantled by residues Thr‐8, and −18, Ser‐11, Cys‐13, and Arg‐22. Potential of mean force (PMF) profiles for mono (Na⁺, K⁺, Cl⁻) and divalent (Ca²⁺) ions along the pore are calculated. The data support experimental findings of a weak cation selectivity of the channel. Calculations on 8a are compared to data derived for a pentameric bundle consisting of the M2 helices of the bacterial pentameric ligand gated ion channel GLIC (3EHZ). PMF curves of both, bundles 8a and M2, show sigmoidal shaped profiles. In comparison to the data for the M2‐GLIC model, data of the 8a bundle show lower amplitude of the PMF values between maximum and minimum and less discrimination amongst ions. Proteins 2015; 83:300–308. © 2014 Wiley Periodicals, Inc.

Insights into the binding modes of CC chemokine receptor 4 (CCR4) inhibitors: a combined approach involving homology modelling, docking, and molecular dynamics simulation studies

CC chemokine receptor 4 (CCR4), a G protein-coupled receptor (GPCR), plays a vital role in the progression of asthma, T-cell lymphoma, inflammation, and Alzheimer's disease. To date, the structure of CCR4 has not been determined. Therefore, the nature of the interactions between inhibitors and CCR4 is not well known. In this study, we used CCR5 as a template to model the structure of CCR4. Docking studies were performed for four naphthalene-sulphonamide derivatives and crucial ligand-protein interactions were analysed. Molecular dynamics (MD) simulations of these complexes (100 ns each) were carried out to gain insights into the interactions between ligands and CCR4. MD simulations revealed that the residues identified by the docking were displaced and new residues were inserted near the ligands. Results of a principal component analysis (PCA) suggested that CCR4 unfolds at the extracellular site surrounding the ligands. Our simulations identified crucial residues involved in CCR4 antagonism, which were supported by previous mutational studies. Additionally, we identified Ser3.29, Leu3.33, Ser5.39, Phe6.47, Ile7.35, Thr7.38, Thr7.40, and Ala7.42 as residues that play crucial roles in CCR4 antagonism. Mutational studies will help elucidate the significance of these residues in CCR4 antagonism. An understanding of ligand-CCR4 interactions might aid in the design of novel CCR4 inhibitors.

Docking assay of small molecule antivirals to p7 of HCV

Protein p7 of HCV is a 63 amino acid channel forming membrane protein essential for the progression of viral infection. With this momentousness, p7 emerges as an important target for antiviral therapy. A series of small molecule drugs, such as amantadine, rimantadine, amiloride, hexamethylene amiloride, NN-DNJ and BIT225 have been found to affect the channel activity. These compounds are docked against monomeric and hexameric structures of p7 taken at various time steps from a molecular dynamics simulation of the protein embedded in a hydrated lipid bilayer. The energetics of binding identifies the guanidine based ligands as the most potent ligands. The adamantanes and NN-DNJ show weaker binding energies. The lowest energy poses are those at the site of the loop region for the monomer and hexamer. For the latter, the poses show a tendency of the ligand to face the lumen of the pore. The mode of binding is that of a balance between hydrophobic interactions and hydrogen bond formation with backbone atoms of the protein.