Hao-Jen Hsu

PEDV ORF3 encodes an ion channel protein and regulates virus production

Several studies suggest that the open reading frame 3 (ORF3) gene of porcine epidemic diarrhea virus (PEDV) is related to viral infectivity and pathogenicity, but its function remains unknown. Here, we propose a structure model of the ORF3 protein consisting of four TM domains and forming a tetrameric assembly. ORF3 protein can be detected in PEDV‐infected cells and it functions as an ion channel in both Xenopus laevis oocytes and yeast. Mutation analysis showed that Tyr170 in TM4 is important for potassium channel activity. Furthermore, viral production is reduced in infected Vero cells when ORF3 gene is silenced by siRNA. Interestingly, the ORF3 gene from an attenuated PEDV encodes a truncated protein with 49 nucleotide deletions, which lacks the ion channel activity.

Nucleolin is important for Epstein–Barr virus nuclear antigen 1-mediated episome binding, maintenance, and transcription

Significance Epstein–Barr Virus (EBV) genome persistence and expression in human cells can result in human lymphocyte or epithelial cell malignancies. EBV genome persistence and gene expression depends on the EBV encoded Epstein–Barr virus nuclear antigen 1 (EBNA1) protein binding to cognate sites in EBV episomes. Search for a cell protein necessary for EBNA1 interaction with EBV episomes identified a low-level nucleolin (NCL) association with EBNA1. Despite low-level stable association, NCL shRNA knockdowns revealed wild-type NCL, but not an NCL ATP binding site point mutant to be critical for EBNA1 binding to EBV episomes, transcription, and genome maintenance. These experiments identify NCL and NCL K429 as potential targets for inhibition of EBNA1’s role in EBV oncogenesis. Epstein–Barr virus (EBV) nuclear antigen 1 (EBNA1) is essential for EBV episome maintenance, replication, and transcription. These effects are mediated by EBNA1 binding to cognate oriP DNA, which comprise 20 imperfect copies of a 30-bp dyad symmetry enhancer and an origin for DNA replication. To identify cell proteins essential for these EBNA1 functions, EBNA1 associated cell proteins were immune precipitated and analyzed by liquid chromatography-tandem mass spectrometry. Nucleolin (NCL) was identified to be EBNA1 associated. EBNA1s N-terminal 100 aa and NCLs RNA-binding domains were critical for EBNA1/NCL interaction. Lentivirus shRNA-mediated NCL depletion substantially reduced EBNA1 recruitment to oriP DNA, EBNA1-dependent transcription of an EBV oriP luciferase reporter, and EBV genome maintenance in lymphoblastoid cell lines. NCL RNA-binding domain K429 was critical for ATP and EBNA1 binding. NCL overexpression increased EBNA1 binding to oriP and transcription, whereas NCL K429A was deficient. Moreover, NCL silencing impaired lymphoblastoid cell line growth. These experiments reveal a surprisingly critical role for NCL K429 in EBNA1 episome maintenance and transcription, which may be a target for therapeutic intervention.

Viral channel forming proteins - modeling the target.

Abstract 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.

In Silico Analysis Reveals Sequential Interactions and Protein Conformational Changes during the Binding of Chemokine CXCL-8 to Its Receptor CXCR1

Chemokine CXCL-8 plays a central role in human immune response by binding to and activate its cognate receptor CXCR1, a member of the G-protein coupled receptor (GPCR) family. The full-length structure of CXCR1 is modeled by combining the structures of previous NMR experiments with those from homology modeling. Molecular docking is performed to search favorable binding sites of monomeric and dimeric CXCL-8 with CXCR1 and a mutated form of it. The receptor-ligand complex is embedded into a lipid bilayer and used in multi ns molecular dynamics (MD) simulations. A multi-steps binding mode is proposed: (i) the N-loop of CXCL-8 initially binds to the N-terminal domain of receptor CXCR1 driven predominantly by electrostatic interactions; (ii) hydrophobic interactions allow the N-terminal Glu-Leu-Arg (ELR) motif of CXCL-8 to move closer to the extracellular loops of CXCR1; (iii) electrostatic interactions finally dominate the interaction between the N-terminal ELR motif of CXCL-8 and the EC-loops of CXCR1. Mutation of CXCR1 abrogates this mode of binding. The detailed binding process may help to facilitate the discovery of agonists and antagonists for rational drug design.

ORF8a of SARS-CoV forms an ion channel: Experiments and molecular dynamics simulations

Abstract ORF8a protein is 39 residues long and contains a single transmembrane domain. The protein is synthesized using solid phase peptide synthesis and reconstituted into artificial lipid bilayers that forms cation-selective ion channels with a main conductance level of 8.9±0.8pS at elevated temperature (38.5°C). Computational modeling studies including multi nanosecond molecular dynamics simulations in a hydrated POPC lipid bilayer are done with a 22 amino acid transmembrane helix to predict a putative homooligomeric helical bundle model. A structural model of a pentameric bundle is proposed with cysteines, serines and threonines facing the pore.

Assembling viral channel forming proteins: Vpu from HIV-1.

Different routes of assembly are probed for the transmembrane domain (TMD) of the bitopic membrane protein Vpu from HIV-1. Vpu is responsible for the amplification of viral release from the host cell. The mode of action includes (i) heteroassembly with host factors and (ii) the formation of homo-oligomers, which are able to conduct ions across the lipid membrane. Two different routes of assembling short sequences of the N terminus, including the TMD of Vpu, Vpu1-32, and Vpu8-26, are presented by using a combination of classical molecular dynamics (MD) simulations combined with a docking approach. The rim of alanines (Ala-8, -11, -15, and -19) resembles an interlocking motif for the sequential assembly into a dimer and trimer. Simultaneous assembly results in oligomeric bundles (trimers to pentamers) with either tryptophans (Trp-23) or purely hydrophobic residues facing the center. Bundles, with serines facing the pore (Ser-24), are energetically not the lowest structures. For pentameric bundles with Ser-24 facing the pore, no water column develops during a short 25 ns MD simulation.

Membrane undulation induced by NS4A of Dengue virus: a molecular dynamics simulation study

Nonstructural protein 4A (NS4A) of Dengue virus (DENV) is a membrane protein involved in rearrangements of the endoplasmic reticulum membrane that are required for formation of replication vesicles. NS4A is composed most likely of three membrane domains. The N- and C-terminal domains are supposed to traverse the lipid membrane whereas the central one is thought to reside on the membrane surface, thus forming a u-shaped protein. All three membrane domains are proposed to be helical by secondary structure prediction programs. After performing multi nanosecond molecular dynamics (MD) simulations at various temperatures (300, 310, and 315.15 K) with each of the individual domains, they are used in a docking approach to define putative association motifs of the transmembrane domains (TMDs). Two structures of the u-shaped protein are generated by separating two assembled TMDs linking them with the membrane-attached domain. Lipid undulation is monitored with the structures embedded in a fully hydrated lipid bilayer applying multiple 200 ns MD simulations at 310 K. An intact structure of the protein supports membrane undulation. The strong unwinding of the helices in the domain-linking section of one of the structures lowers its capability to induce membrane curvature. Unwinding of the link region is due to interactions of two tryptophan residues, Trp-96 and 104. These results provide first insights into the membrane-altering properties of DENV NS4A.

Computational modeling of the p7 monomer from HCV and its interaction with small molecule drugs.

Hepatitis C virus p7 protein is a 63 amino acid polytopic protein with two transmembrane domains (TMDs) and one of the prime targets for anti HCV drug development. A bio-inspired modeling pathway is used to generate plausible computational models of the two TMDs forming the monomeric protein model. A flexible region between Leu-13 and Gly-15 is identified for TMD11-32 and a region around Gly-46 to Trp-48 for TMD236-58. Mutations of the tyrosine residues in TMD236-58 into phenylalanine and serine are simulated to identify their role in shaping TMD2. Lowest energy structures of the two TMDs connected with the loop residues are used for a posing study in which small molecule drugs BIT225, amantadine, rimantadine and N N-DNJ, are identified to bind to the loop region. BIT225 is identified to interact with the backbone of the functionally important residues Arg-35 and Trp-36.

Peptides derived from CXCL8 based on in silico analysis inhibit CXCL8 interactions with its receptor CXCR1.

Chemokine CXCL8 is crucial for regulation of inflammatory and immune responses via activating its cognate receptor CXCR1. In this study, molecular docking and binding free energy calculations were combined to predict the initial binding event of CXCL8 to CXCR1 for peptide drug design. The simulations reveal that in the initial binding, the N-loop of CXCL8 interacts with the N-terminus of CXCR1, which is dominated by electrostatic interactions. The derived peptides from the binding region of CXCL8 are synthesized for further confirmation. Surface plasmon resonance analyses indicate that the CXCL8 derived peptide with 14 residues is able to bind to the receptor CXCR1 derived peptide with equilibrium KD of 252 μM while the peptide encompassing a CXCL8 K15A mutation hardly binds to CXCR1 derived peptide (KD = 1553 μM). The cell experiments show that the designed peptide inhibits CXCL8-induced and LPS-activated monocytes adhesion and transmigration. However, when the peptides were mutated on two lysine residues (K15 and K20), the inhibition effects were greatly reduced indicating these two amino acids are key residues for the initial binding of CXCL8 to CXCR1. This study demonstrates that in silico prediction based functional peptide design can be effective for developing anti-inflammation drugs.

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 (Ca2+) 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.