Structure and functional dynamics characterization of the ion channel of the human respiratory syncytial virus (hRSV) small hydrophobic protein (SH) transmembrane domain by combining molecular dynamics with excited normal modes

The human respiratory syncytial virus (hRSV) is the major cause of lower respiratory tract infection in children and elderly people worldwide. Its genome encodes 11 proteins including SH protein, whose functions are not well known. Studies show that SH protein increases RSV virulence degree and permeability to small compounds, suggesting it is involved in the formation of ion channels. The knowledge of SH structure and function is fundamental for a better understanding of its infection mechanism. The aim of this study was to model, characterize, and analyze the structural behavior of SH protein in the phospholipids bilayer environment. Molecular modeling of SH pentameric structure was performed, followed by traditional molecular dynamics (MD) simulations of the protein immersed in the lipid bilayer. Molecular dynamics with excited normal modes (MDeNM) was applied in the resulting system in order to investigate long time scale pore dynamics. MD simulations support that SH protein is stable in its pentameric form. Simulations also showed the presence of water molecules within the bilayer by density distribution, thus confirming that SH protein is a viroporin. This water transport was also observed in MDeNM studies with histidine residues of five chains (His22 and His51), playing a key role in pore permeability. The combination of traditional MD and MDeNM was a very efficient protocol to investigate functional conformational changes of transmembrane proteins that act as molecular channels. This protocol can support future investigations of drug candidates by acting on SH protein to inhibit viral infection.

Open image in new window

Graphical Abstract The ion channel of the human respiratory syncytial virus (hRSV) small hydrophobic protein (SH) transmembrane domain?

     

Modeling the structure of the respiratory syncytial virus small hydrophobic protein by silent-mutation analysis of global searching molecular dynamics

Human respiratory syncytial virus (RSV) encodes a small hydrophobic (SH) protein, whose function in the life cycle of the virus is unknown. Recent channel activity measurements of the protein suggest that like other viroporins, SH may assemble into a homo-oligomeric ion channel. To further our understanding of this potentially important protein, a new strategy was implemented in order to model the transmembrane oligomeric bundle of the protein. Global searching molecular dynamic simulations of SH proteins from eight different viral strains, each at different oligomeric states, as well as different lengths of the putative transmembrane domain, were undertaken. Taken together, a total of 45 different global molecular dynamic simulations pointed to a single pentameric structure for the protein that was found in all of the different variants. The model of the structure obtained is a channel-like homopentamer whose minimal transmembrane pore diameter is 1.46 A.     

The small hydrophobic protein of the human respiratory syncytial virus forms pentameric ion channels

The small hydrophobic (SH) protein is encoded by the human respiratory syncytial virus. Its absence leads to viral attenuation in the context of whole organisms, and it prevents apoptosis in infected cells. Herein, we have examined the structure of SH protein in detergent micelles and in lipid bilayers, by solution NMR and attenuated total reflection-Fourier transform infrared spectroscopy, respectively. We found that SH protein has a single α-helical transmembrane domain and forms homopentamers in several detergents. In detergent micelles, the transmembrane domain is flanked N-terminally by an α-helix that forms a ring around the lumen of the pore and C-terminally by an extended β-turn. SH protein was found in the plasma membrane of transiently expressing HEK 293 cells, which showed pH-dependent (acid-activated) channel activity. Channel activity was abolished in mutants lacking both native His residues, His(22) and His(51), but not when either His was present. Herein, we propose that the pentameric model of SH protein presented is a physiologically relevant conformation, albeit probably not the only one, in which SH contributes to RSV infection and replication. Viroporins are short (～100 amino acids) viral membrane proteins that form oligomers of a defined size, act as proton or ion channels, and in general enhance membrane permeability in the host. However, with some exceptions, their precise biological role of their channel activity is not understood. In general, viroporins resemble poorly specialized proteins but are nevertheless critical for viral fitness. In vivo, viruses lacking viroporins usually exhibit an attenuated or weakened phenotype, altered tropism, and diminished pathological effects. We have chosen to study the SH protein, 64 amino acids long, found in the human respiratory syncytial virus because of the effect of RSV on human health and the lack of adequate antivirals. We show that SH protein forms oligomers that behave as ion channels when activated at low pH. This study adds SH protein to a growing group of viroporins that have been structurally characterized. Although the precise biological role of this pentameric channel is still unknown, this report is nevertheless essential to fill some of the many gaps that exist in the understanding of SH protein function.     

Function of the respiratory syncytial virus small hydrophobic protein

Respiratory syncytial virus (RSV), a member of the Paramyxoviridae family, encodes a small hydrophobic (SH) protein of unknown function. Parainfluenza virus 5 (PIV5), a prototypical paramyxovirus, also encodes an SH protein, which inhibits tumor necrosis factor alpha (TNF-α) signaling. In this study, recombinant PIV5 viruses without their own SH but containing RSV SH (from RSV strain A2 or B1) in its place (PIV5ΔSH-RSV SH) and RSV lacking its own SH (RSVΔSH) were generated and analyzed. The results indicate that the SH protein of RSV has a function similar to that of PIV5 SH and that it can inhibit TNF-α signaling.     

Direct visualization of the small hydrophobic protein of human respiratory syncytial virus reveals the structural basis for membrane permeability

Human respiratory syncytial virus (HRSV) is the leading cause of lower respiratory tract disease in infants. The HRSV small hydrophobic (SH) protein plays an important role in HRSV pathogenesis, although its mode of action is unclear. Analysis of the ability of SH protein to induce membrane permeability and form homo-oligomers suggests it acts as a viroporin. For the first time, we directly observed functional SH protein using electron microscopy, which revealed SH forms multimeric ring-like objects with a prominent central stained region. Based on current and existing functional data, we propose this region represents the channel that mediates membrane permeability.

Structured summary

MINT-7890792, MINT-7890805: SH (uniprotkb:P04852) and SH (uniprotkb:P04852) bind (MI:0407) by chromatography technology (MI:0091)MINT-7890784, MINT-7890776: SH (uniprotkb:P04852) and SH (uniprotkb:P04852) bind (MI:0407) by electron microscopy (MI:0040)     

The passage of an ion through a synthetic ion channel

The simulated system consisted of a fatty acid bilayer membrane dividing two electrolyte layers each containing ions, and a channel composed of linked 15-crown-5 ether rings. The Na⁺ and F^?  ions in the aqueous electrolyte layers were too large to enter the channel, but the Li⁺ ions entered and were transported. Conditions that optimised the passive, electric-field-induced transport of Li⁺ ions through the channel were investigated. It was calculated and rationalised that the higher the numerical value of the electrostatic charge on the oxygen atoms of the crown ether rings, the more easily does the channel convey the Li⁺ ions.     

Cystic fibrosis transmembrane conductance regulator (CFTR): Making an ion channel out of an active transporter structure

Cystic fibrosis is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR). CFTR is a member of the ATP-binding cassette (ABC) family of membrane transport proteins, most members of which function as ATP-dependent pumps. CFTR is unique among human ABC proteins in functioning not as a pump, but as an ion channel. Recent structural data has indicated that CFTR shares broadly similar overall architecture and ATP-dependent conformational changes as other ABC proteins. Functional investigations suggest that CFTR has a unique open portal connecting the cytoplasm to the transmembrane channel pore, that allows for a continuous pathway for Cl^? ions to cross the membrane in one conformation. This lateral portal may be what allows CFTR to function as an ion channel rather than as a pump, suggesting a plausible mechanism by which channel function may have evolved in CFTR.     

Simulations of ion current in realistic models of ion channels: the KcsA potassium channel

Realistic studies of ion current in biologic channels present a major challenge for computer simulation approaches. All-atom molecular dynamics simulations involve serious time limitations that prevent their use in direct evaluation of ion current in channels with significant barriers. The alternative use of Brownian dynamics (BD) simulations can provide the current for simplified macroscopic models. However, the time needed for accurate calculations of electrostatic energies can make BD simulations of ion current expensive. The present work develops an approach that overcomes some of the above challenges and allows one to simulate ion currents in models of biologic channels. Our method provides a fast and reliable estimate of the energetics of the system by combining semimacroscopic calculations of the self-energy of each ion and an implicit treatment of the interactions between the ions, as well as the interactions between the ions and the protein-ionizable groups. This treatment involves the use of the semimacroscopic version of the protein dipole Langevin dipole (PDLD/S) model in its linear response approximation (LRA) implementation, which reduces the uncertainties about the value of the protein "dielectric constant." The resulting free energy surface is used to generate the forces for on-the-fly BD simulations of the corresponding ion currents. Our model is examined in a preliminary simulation of the ion current in the KcsA potassium channel. The complete free energy profile for a single ion transport reflects reasonable energetics and captures the effect of the protein-ionized groups. This calculated profile indicates that we are dealing with the channel in its closed state. Reducing the barrier at the gate region allows us to simulate the ion current in a reasonable computational time. Several limiting cases are examined, including those that reproduce the observed current, and the nature of the productive trajectories is considered. The ability to simulate the current in realistic models of ion channels should provide a powerful tool for studies of the biologic function of such systems, including the analysis of the effect of mutations, pH, and electric potentials.     

A structure for the trimeric MHC class II-associated invariant chain transmembrane domain

The major histocompatibility complex (MHC)-associated invariant chain (Ii) contains a single transmembrane domain that forms trimers. Ii is involved in the assembly of the MHC and antigen presentation, and is thus central to the function of the immune system. Here, we show by attenuated total reflectance, Fourier transform infrared (ATR-FTIR) spectroscopy that the transmembrane domain is alpha-helical and we provide a structural model of the transmembrane domain obtained by a combination of site-specific infrared dichroism and molecular dynamics (MD) simulations. This work resolves the backbone structure of a transmembrane peptide by multiple (13)C=(18)O labelling at ten different residues. A second purely computational approach, based on MD simulations of Ii transmembrane homologous sequences, yields a similar structure that is consistent with our experimental results. The structure presented forms a left-handed coiled coil with an average helix tilt of 13(+/-6) degrees; the residue Gln47 implicated in trimer formation forms strong interhelical contacts, Thr50 points to the inside of the trimeric coil and forms a network of hydrogen bonds.     

机械敏感性离子通道蛋白Piezo2的研究进展

2010年科学家在小鼠神经瘤母细胞中筛选鉴定出Piezo1和Piezo2蛋白.Piezo2是一个由机械刺激直接激活且可将机械刺激转换为电信号进而形成机械敏感性电流通道的蛋白质.Piezo2自发现以来一直受到广泛的关注,在触觉、本体感觉、痛觉、肿瘤癌症等多种生理病理过程中发挥重要作用.本文在前期研究的基础上阐述了机械敏感...     

Terminal ion pairs stabilize the second beta-hairpin of the B1 domain of protein G

The effects of terminal ion pairs on the stability of a beta-hairpin peptide corresponding to the C-terminal residues of the B1 domain of protein G were determined using thermal unfolding as monitored by nuclear magnetic resonance and circular dichroism spectroscopy. Molecular dynamics (MD) simulations were also performed to examine the effect of ion pairs on the structures. Eight peptides were studied including the wild type (G41) and the N-terminal modified sequences that had the first residue deleted (E42), replaced with a Lys (K41), or extended by an additional Gly (G40). Acetylated variants were made to examine the effect of removing the positive N-terminal charge on beta-hairpin stability. The rank in stability determined experimentally is K41 > E42 approximately G41 approximately G40 > Ac-K41 > Ac-E42 approximately Ac-G41 > Ac-G40. The Tm of the K41 peptide is 12 degrees C higher than G41, while the Tm values for the acetylated peptides are less than their unacetylated forms by more than 15 degrees C. NOE cross-peaks between side-chain methylene groups at the N- and C-termini and larger CalphaH shifts compared to random values are seen for K41. The addition of 20% methanol increases the stability in K41 and G41. The MD studies complement these results by showing that the charged N-terminus is important to stability. The type of ion pair observed varies with peptide, and when formed the simulations show that the ion pair can prevent fraying of the beta-strands through electrostatic and hydrophobic contacts. Therefore, introducing favorable electrostatic interactions at the N- and C-termini can substantially enhance beta-hairpin stability and help define the structure.     

Small single transmembrane domain (STMD) proteins organize the hydrophobic subunits of large membrane protein complexes

The large membrane protein complexes of mitochondrial oxidative phosphorylation are composed of central subunits that are essential for their bioenergetic core function and accessory subunits that may assist in regulation, assembly or stabilization. Although sequence conservation is low, a significant proportion of the accessory subunits is characterized by a common single transmembrane (STMD) topology. The STMD signature is also found in subunits of other membrane protein complexes. We hypothesize that the general function of STMD subunits is to organize the hydrophobic subunits of large membrane protein complexes in specialized environments like the inner mitochondrial membrane.     

Structure and functional characterization of the periplasmic N-terminal polypeptide domain of the sugar-specific ion channel protein (ScrY porin)

The structure of the sucrose-specific porin (ScrY) from Salmonella typhimurium has been elucidated by X-ray crystallography to consist of 18 antiparallel beta-strands, associated as a trimer complex similar to ion-transport channels. However, the 71-amino-acid-residue N-terminal periplasmic domain was not determined from the crystal structure due to the absence of sufficient electron density. The N-terminal polypeptide contains a coiled-coil structural motif and has been assumed to play a role in the sugar binding of ScrY porin. In this study the proteolytic stability and a specific proteolytic truncation site at the N-terminal domain were identified by the complete primary structure characterization of ScrY porin, using MALDI mass spectrometry and post-source-decay fragmentation. The secondary structure and supramolecular association of the coiled-coil N-terminal domain were determined by chemical synthesis of the complete N-terminal polypeptide and several partial sequences and their spectroscopic, biophysical, and mass spectrometric characterization. Circular dichroism spectra revealed predominant alpha-helical conformation for the putative coiled-coil domain comprising residues 4-46. Specific association to both dimer and trimer complexes was identified by electrospray ionization mass spectra and was ascertained by dynamic light scattering and electrophoresis data. The role of the N-terminal domain in sugar binding was examined by comparative TR-NOE-NMR spectroscopy of the complete ScrY porin and a recombinant mutant, ScrY(delta1-62), lacking the N-terminal polypeptide. The TR-NOE-NMR data showed a strong influence of ScrY porin on the sugar-binding affinity and suggested a possible function of the periplasmic N terminus for supramolecular stabilization and low-affinity sugar binding.     

Expression and purification of human respiratory syncytial virus recombinant fusion protein

The Human Respiratory Syncytial Virus (HRSV) fusion protein (F) was expressed in Escherichia coli BL21A using the pET28a vector at 37 °C. The protein was purified from the soluble fraction using affinity resin. The structural quality of the recombinant fusion protein and the estimation of its secondary structure were obtained by circular dichroism. Structural models of the fusion protein presented 46% of the helices in agreement with the spectra by circular dichroism analysis. There are only few studies that succeeded in expressing the HRSV fusion protein in bacteria. This is a report on human fusion protein expression in E. coli and structure analysis, representing a step forward in the development of fusion protein F inhibitors and the production of antibodies.     

Determinants of anion permeation in the second transmembrane domain of the mouse bestrophin-2 chloride channel

Bestrophins have been proposed to constitute a new family of Cl channels that are activated by cytosolic Ca. We showed previously that mutation of serine-79 to cysteine in mouse bestrophin-2 (mBest2) altered the relative permeability and conductance to SCN. In this paper, we have overexpressed various mutant constructs of mBest2 in HEK-293 cells to explore the contributions to anion selectivity of serine-79 and other amino acids (V78, F80, G83, F84, V86, and T87) located in the putative second transmembrane domain (TMD2). Residues selected for mutagenesis were distributed throughout TMD2, but mutations at all positions changed the selectivity. The effects on selectivity were rather modest. Replacement of residues 78, 79, 80, 83, 84, 86, or 87 with cysteine had similar effects: the permeability of the channel to SCN relative to Cl (PSCN/PCl) was decreased three- to fourfold and the relative SCN conductance (GSCN/GCl) was increased five- to tenfold. Side chains at positions 78 and 80 appeared to be situated close to the permeant anion, because the electrostatic charge at these positions affected permeation in specific ways. The effects of charged sulfhydryl-reactive MTS reagents were the opposite in the V78C and F80C mutants and the effects were partially mimicked by substitution of F80 with charged amino acids. In S79T, switching from Cl to SCN caused slow changes in GSCN/GCl (tau = 16.6 s), suggesting that SCN binding to the channel altered channel gating as well as conductance. The data in this paper and other data support a model in which TMD2 plays an important role in forming the bestrophin pore. We suggest that the major determinant in anion permeation involves partitioning of the permeant anion into an aqueous pore whose structural features are rather flexible. Furthermore, anion permeation and gating may be linked.     

The strong dimerization of the transmembrane domain of the fibroblast growth factor receptor (FGFR) is modulated by C‐terminal juxtamembrane residues

The fibroblast growth factor receptor 3 (FGFR3) is a member of the FGFR subfamily of the receptor tyrosine kinases (RTKs) involved in signaling across the plasma membrane. Generally, ligand binding leads to receptor dimerization and activation. Dimerization involves the transmembrane (TM) domain, where mutations can lead to constitutive activation in certain cancer types and also in skeletal malformations. Thus, it has been postulated that FGFR homodimerization must be inherently weak to allow regulation, a feature reminiscent of α and β integrin TM interactions. However, we show herein that in FGFR3‐TM, four C‐terminal residues, CRLR, have a profound destabilizing effect in an otherwise strongly dimerizing TM peptide. In the absence of these four residues, the dimerizing propensity of FGFR3‐TM is comparable to glycophorin, as shown using various detergents. In addition, the expected enhanced dimerization induced by the mutation associated to the Crouzon syndrome A391E, was observed only when these four C‐terminal residues were present. In the absence of these four residues, A391E was dimer‐destabilizing. Finally, using site specific infrared dichroism and convergence with evolutionary conservation data, we have determined the backbone model of the FGFR3‐TM homodimer in model lipid bilayers. This model is consistent with, and correlates with the effects of, most known pathological mutations found in FGFR‐TM.     

Molecular dynamics simulation and experimental validation of the effect of pH on protein desorption in hydrophobic charge induction chromatography

The ligands in hydrophobic charge induction chromatography (HCIC) are hydrophobic and ionisable. Thus, the pH is crucial for the separation performance in HCIC, especially for elution. However, it is difficult to obtain the microscopic information in HCIC through experimental means. In this work, molecular dynamics simulations are performed to examine the effect of pH on elution and protein conformational transition in HCIC, using a 46-bead β-barrel coarse-grained model protein and an HCIC adsorbent pore model constructed in an earlier work. Corresponding experiments are carried out for the validation of simulation results, using lysozyme and MEP Hypercel. Both the activities and fluorescence of lysozyme are examined to evaluate the conformational transition. The simulations indicate that the elution efficiency of protein increases with decreasing pH value in a non-linear manner. This is qualitatively consistent with the experimental results. MD simulations indicate that protein unfolding occurs in elution at all pH values. However, the experimental data show that the activity and conformation of lysozyme is independent of pH of the elution buffer. The microscopic information from simulation shows that protein unfolding is mainly observed on the adsorbent surface, but it cannot be detected in the experiments that only probe the proteins in the bulk liquid.     

Conformation and dynamics of the three-helix bundle UBA domain of p62 from experiment and simulation

The ubiquitin associated domain of p62 is a small three-helix bundle of approximately 50 residues that mediates the recognition of polyubiquitin chains and ubiquitylated substrates. The solution structure of a 52 residue construct containing this domain has been characterized using heteronuclear nuclear magnetic resonance (NMR) methods. The resulting ensemble of NMR-derived structures was used in molecular dynamics (MD) simulations to investigate the equilibrium conformation and dynamics of this domain. NOE and (15)N relaxation data have been used to validate the structural ensemble produced by the MD simulations and show a good correlation for residues in regions of secondary structure. A similar approach was taken using an ensemble of structures from the MD simulations to calculate electronic circular dichroism (CD) and IR spectra from first principles with an encouraging correlation with the experimental CD and IR data.     

Chaperone-like activity of the N-terminal region of a human small heat shock protein and chaperone-functionalized nanoparticles

Small heat shock proteins (sHsps) are molecular chaperones employed to interact with a diverse range of substrates as the first line of defense against cellular protein aggregation. The N-terminal region (NTR) is implicated in defining features of sHsps; notably in their ability to form dynamic and polydisperse oligomers, and chaperone activity. The physiological relevance of oligomerization and chemical-scale mode(s) of chaperone function remain undefined. We present novel chemical tools to investigate chaperone activity and substrate specificity of human HspB1 (B1NTR), through isolation of B1NTR and development of peptide-conjugated gold nanoparticles (AuNPs). We demonstrate that B1NTR exhibits chaperone capacity for some substrates, determined by anti-aggregation assays and size-exclusion chromatography. The importance of protein dynamics and multivalency on chaperone capacity was investigated using B1NTR-conjugated AuNPs, which exhibit concentration-dependent chaperone activity for some substrates. Our results implicate sHsp NTRs in chaperone activity, and demonstrate the therapeutic potential of sHsp-AuNPs in rescuing aberrant protein aggregation.     

The transmembrane domain of the respiratory syncytial virus F protein is an orientation-independent apical plasma membrane sorting sequence

The processes that facilitate transport of integral membrane proteins though the secretory pathway and subsequently target them to particular cellular membranes are relevant to almost every field of biology. These transport processes involve integration of proteins into the membrane of the endoplasmic reticulum (ER), passage from the ER to the Golgi, and post-Golgi trafficking. The respiratory syncytial virus (RSV) fusion (F) protein is a type I integral membrane protein that is uniformly distributed on the surface of infected nonpolarized cells and localizes to the apical plasma membrane of polarized epithelial cells. We expressed wild-type or altered RSV F proteins to gain a better understanding of secretory transport and plasma membrane targeting of type I membrane proteins in polarized and nonpolarized epithelial cells. Our findings reveal a novel, orientation-independent apical plasma membrane targeting function for the transmembrane domain of the RSV F protein in polarized epithelial cells. This work provides a basis for a more complete understanding of the role of the transmembrane domain and cytoplasmic tail of viral type I integral membrane proteins in secretory transport and plasma membrane targeting in polarized and nonpolarized cells.