Exploring the conformational space of Vpu from HIV-1: a versatile adaptable protein

The cut-off effect of n-alcohols in lipid rafts: A lipid-dependent phenomenon☆

n-Aliphatic alcohols act as anesthetics only up to a certain chain length, beyond which its biological activity disappears. This is known as the 'cut-off' phenomenon. Although the most accepted explanation is based on action sites in membrane proteins, it is not well understood why alcohols alter their functions. The structural dependence of these protein receptors to lipid domains known as 'lipid rafts', suggests a new approach to tackle the puzzling phenomenon. In this work, by performing molecular dynamic simulations (MDS) to explore the lipid role, we provide relevant molecular details about the membrane-alcohol interaction at the cut-off point regime. Since the high variability of the cut-off points found on protein receptors in neurons may be a consequence of differences in the lipid composition surrounding such proteins, our results could have a clear-cut importance.

TPC2 polymorphisms associated with a hair pigmentation phenotype in humans result in gain of channel function by independent mechanisms

Two-pore channels (TPCs) are endolysosomal cation channels. Two members exist in humans, TPC1 and TPC2. Functional roles associated with the ubiquitously expressed TPCs include VEGF-induced neoangiogenesis, LDL-cholesterol trafficking and degradation, physical endurance under fasting conditions, autophagy regulation, the acrosome reaction in sperm, cancer cell migration, and intracellular trafficking of pathogens such as Ebola virus or bacterial toxins (e.g., cholera toxin). In a genome-wide association study for variants associated with human pigmentation characteristics two coding variants of TPC2, rs35264875 (encoding M484L) and rs3829241 (encoding G734E), have been found to be associated with a shift from brown to blond hair color. In two recent follow-up studies a role for TPC2 in pigmentation has been further confirmed. However, these human polymorphic variants have not been functionally characterized until now. The development of endolysosomal patch-clamp techniques has made it possible to investigate directly ion channel activities and characteristics in isolated endolysosomal organelles. We applied this technique here to scrutinize channel characteristics of the polymorphic TPC2 variants in direct comparison with WT. We found that both polymorphisms lead to a gain of channel function by independent mechanisms. We next conducted a clinical study with more than 100 blond- and brown/black-haired individuals. We performed a genotype/phenotype analysis and subsequently isolated fibroblasts from WT and polymorphic variant carriers for endolysosomal patch-clamp experimentation to confirm key in vitro findings.

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.

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.

Interaction of ions with the luminal sides of wild-type and mutated skeletal muscle ryanodine receptors

Ryanodine receptors (RyRs) are the largest known ion channels, and are of central importance for the release of Ca(2+) from the sarco/endoplasmic reticulum (SR/ER) in a variety of cells. In cardiac and skeletal muscle cells, contraction is triggered by the release of Ca(2+) into the cytoplasm and thus depends crucially on correct RyR function. In this work, in silico mutants of the RyR pore were generated and MD simulations were conducted to examine the impact of the mutations on the Ca(2+) distribution. The Ca(2+) distribution pattern on the luminal side of the RyR was most affected by G4898R, D4899Q, E4900Q, R4913E, and D4917A mutations. MD simulations with our wild-type model and various ion species showed a preference for Ca(2+) over other cations at the luminal pore entrance. This Ca(2+)-accumulating characteristic of the luminal RyR side may be essential to the conductance properties of the channel.

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.

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.

The interaction of local anesthetics with lipid membranes

Molecular Dynamic Simulations are performed to evaluate the interaction of lidocaine, procaine and tetracaine with a lipid membrane. The main interest is to evaluate the structural changes produced by these local anesthetics in the bilayers. Penetration trajectories, interaction energies, entropy changes and an order parameter are calculated to quantify the destabilization of the lipid configurations. We show that such structural parameters give important information to understand how anesthetic agents influence the structure of plasma membranes. Graphic processing units (GPUs) are used in our simulations.

Qualitative computational bioanalytics: assembly of viral channel-forming peptides around mono and divalent ions

A fine-grained docking protocol was used to generate a bundle-like structure of the bitopic membrane protein Vpu from HIV-1. Vpu is a type I membrane protein with 81 amino acids. It is proposed that Vpu forms ion- and substrate-conducting bundles, which are located at the plasma membrane in the infected cell. The Vpu1-32 peptide that includes the transmembrane domain (TMD) is assembled into homo-pentameric bundles around prepositioned Na, K, Ca or Cl ions. For bundles with the lowest energy, the TMDs generate a hydrophobic pore. Bundles in which Ser-24 faces the pore have higher energy. The tilt of the helices in the lowest energy bundles is larger than bundles with serines facing the pore. Left-handed bundles are lowest in energy where the ions are located at the serines.

Correlation of biological activity with computationally derived structural features from transmembrane hetero-dimers of HIV-1 Vpu with host factors

Vpu is an 81 amino acid type I integral membrane protein encoded by human immunodeficiency virus type 1 (HIV-1). It is identified to support viral release by potentially forming ion and substrate conducting channels and by modulating the function of host factors. The focus is on the interaction of the transmembrane domains of Vpu with those of host factors using a combination of molecular dynamics simulations and docking approach. Binding poses and adopted tilt angles of the dimers are analyzed and correlated with experimentally derived activity data from literature. Vpu activity is driven by dimerization with the host protein via its alanine rim Ala-8/11/15/19. Tight binding is shown by an almost parallel alignment of the helices in the dimers. Less parallel alignment is proposed to correlate with lower activity. This article is part of a Special Issue entitled: Viral Membrane Proteins - Channels for Cellular Networking.

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 TMD1_1-32 and a region around Gly-46 to Trp-48 for TMD2_36-58. Mutations of the tyrosine residues in TMD2_36-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 NN-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.

Membrane Packing Problems: A short Review on computational Membrane Modeling Methods and Tools

The use of model membranes is currently part of the daily workflow for many biochemical and biophysical disciplines. These membranes are used to analyze the behavior of small substances, to simulate transport processes, to study the structure of macromolecules or for illustrative purposes. But, how can these membrane structures be generated? This mini review discusses a number of ways to obtain these structures. First, the problem will be formulated as the Membrane Packing Problem. It will be shown that the theoretical problem of placing proteins and lipids onto a membrane area differ significantly. Thus, two sub-problems will be defined and discussed. Then, different - partly historical - membrane modeling methods will be introduced. And finally, membrane modeling tools will be evaluated which are able to semi-automatically generate these model membranes and thus, drastically accelerate and simplify the membrane generation process. The mini review concludes with advice about which tool is appropriate for which application case.

The influence of non polar and polar molecules in mouse motile cells membranes and pure lipid bilayers

We report an experimental study of mouse sperm motility that shows chief aspects characteristic of neurons: the anesthetic (produced by tetracaine) and excitatory (produced by either caffeine or calcium) effects and their antagonic action. While tetracaine inhibits sperm motility and caffeine has an excitatory action, the combination of these two substances balance the effects, producing a motility quite similar to that of control cells. We also study the effects of these agents (anesthetic and excitatory) on the melting points of pure lipid liposomes constituted by 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and dipalmitoyl phosphatidic acid (DPPA). Tetracaine induces a large fluidization of the membrane, shifting the liposomes melting transition temperature to much lower values. The effect of caffeine is null, but its addition to tetracaine-doped liposomes greatly screen the fluidization effect. A high calcium concentration stiffens pure lipid membranes and strongly reduces the effect of tetracaine. Molecular Dynamics Simulations are performed to further understand our experimental findings at the molecular level. We find a strong correlation between the effect of antagonic molecules that could explain how the mechanical properties suitable for normal cell functioning are affected and recovered.

The p7 protein of hepatitis C virus forms structurally plastic, minimalist ion channels

Hepatitis C virus (HCV) p7 is a membrane-associated oligomeric protein harboring ion channel activity. It is essential for effective assembly and release of infectious HCV particles and an attractive target for antiviral intervention. Yet, the self-assembly and molecular mechanism of p7 ion channelling are currently only partially understood. Using molecular dynamics simulations (aggregate time 1.2 µs), we show that p7 can form stable oligomers of four to seven subunits, with a bias towards six or seven subunits, and suggest that p7 self-assembles in a sequential manner, with tetrameric and pentameric complexes forming as intermediate states leading to the final hexameric or heptameric assembly. We describe a model of a hexameric p7 complex, which forms a transiently-open channel capable of conducting ions in simulation. We investigate the ability of the hexameric model to flexibly rearrange to adapt to the local lipid environment, and demonstrate how this model can be reconciled with low-resolution electron microscopy data. In the light of these results, a view of p7 oligomerization is proposed, wherein hexameric and heptameric complexes may coexist, forming minimalist, yet robust functional ion channels. In the absence of a high-resolution p7 structure, the models presented in this paper can prove valuable as a substitute structure in future studies of p7 function, or in the search for p7-inhibiting drugs.

Protein intrinsic disorder as a flexible armor and a weapon of HIV-1

Many proteins and protein regions are disordered in their native, biologically active states. These proteins/regions are abundant in different organisms and carry out important biological functions that complement the functional repertoire of ordered proteins. Viruses, with their highly compact genomes, small proteomes, and high adaptability for fast change in their biological and physical environment utilize many of the advantages of intrinsic disorder. In fact, viral proteins are generally rich in intrinsic disorder, and intrinsically disordered regions are commonly used by viruses to invade the host organisms, to hijack various host systems, and to help viruses in accommodation to their hostile habitats and to manage their economic usage of genetic material. In this review, we focus on the structural peculiarities of HIV-1 proteins, on the abundance of intrinsic disorder in viral proteins, and on the role of intrinsic disorder in their functions.

Molecular dynamics simulations and conductance studies of the interaction of VP1 N-terminus from Polio virus and gp41 fusion peptide from HIV-1 with lipid membranes

The icosahedral Polio virus capsid consists of 60 copies of each of the coat proteins VP1, VP2, VP3 and myristolyated VP4 (myrVP4). Catalyzed by the host cell receptor the Polio virus enters the host cell via externalization of myrVP4 and the N terminal part of VP1. There are several assumptions about the individual role of both of the proteins in the mechanism of membrane attachment and genome injection. We use the first 32 N terminal amino acids of VP1 and applied molecular dynamics simulations to assess its mechanism of function when attached and inserted into hydrated lipid membranes (POPC). Helical models are placed in various positions in regard to the lipid membrane to start with. As a comparison, the first 33 amino acids of the fusion peptide of gp41 of HIV-1 are simulated under identical conditions. Computational data support the idea that VP1 is not penetrating into the membrane to form a pore; it rather lays on the membrane surface and only perturbs the membrane. Furthermore, this idea is strengthened by channel recordings of both peptides showing irregular openings.

Mechanism of function of viral channel proteins and implications for drug development

Viral channel-forming proteins comprise a class of viral proteins which, similar to their host companions, are made to alter electrochemical or substrate gradients across lipid membranes. These proteins are active during all stages of the cellular life cycle of viruses. An increasing number of proteins are identified as channel proteins, but the precise role in the viral life cycle is yet unknown for the majority of them. This review presents an overview about these proteins with an emphasis on those with available structural information. A concept is introduced which aligns the transmembrane domains of viral channel proteins with those of host channels and toxins to give insights into the mechanism of function of the viral proteins from potential sequence identities. A summary of to date investigations on drugs targeting these proteins is given and discussed in respect of their mode of action in vivo.

In silico investigations of possible routes of assembly of ORF 3a from SARS-CoV

ORF 3a of human severe acute respiratory syndrome corona virus (SARS-CoV) has been identified as a 274 amino acid membrane protein. When expressed in Xenopus oocytes the protein forms channels. Based on bioinformatics approaches the topology has been identified to include three transmembrane domains (TMDs). Since structural models from experiments are still lacking, computational methods can be challenged to generate such models. In this study, a ‘sequential approach’ for the assembly is proposed in which the individual TMDs are assembled one by one. This protocol is compared with a concerted protocol in which all TMDs are assembled simultaneously. The role of the loops between the TMDs during assembly of the monomers into a bundle is investigated. Molecular dynamics simulations for 20 ns are performed as a short equilibration to assess the bundle stability in a lipid environment. The results suggest that bundles are likely with the second TMD facing the putative pore. All the putative bundles show water molecules trapped within the lumen of the pore with only occasional events of complete crossing.

CELLmicrocosmos 2.2 MembraneEditor: a modular interactive shape-based software approach to solve heterogeneous membrane packing problems

New perspectives have been developed to understand the processes of modeling heterogeneous membranes. These are crucial steps prior to applying advanced techniques like molecular dynamic simulations of whole membrane systems. Lipid, protein, and membrane packing problems are addressed based on biochemical properties in combination with computational optimization techniques. The CELLmicrocosmos 2.2 MembraneEditor (CmME) is introduced as an appropriate framework to handle such problems by offering diverse algorithmic approaches. Its algorithm plug-in-interface enables modelers to generate problem-specific algorithms. Good solutions concerning runtime and lipid density are realized by focusing on the outer shapes of the PDB-based molecules. Application cases are presented like the publication-based modeling of inner and outer mitochondrial membrane-fragments, semiautomatic incorporation of proteins, and the assembly of rafts. Concerning geometrical aspects of the lipids, the achieved results are consistent with experimental observations related to lipid densities and distributions. Finally, two membranes simulated with GROMACS are analyzed and compared: the first is generated with conventional scripting techniques, the second with the CmME Distributor algorithm. The examples prove that CmME is a valuable and versatile tool for a broad set of applications in analysis and visualization of biomembranes.

BST-2 is rapidly down-regulated from the cell surface by the HIV-1 protein Vpu: evidence for a post-ER mechanism of Vpu-action

Recent evidence suggests that transmembrane domain (TMD) interactions are essential for HIV-1 Vpu-mediated antagonism of the restriction factor BST-2/tetherin. We made Vpu TMD mutants to study the mechanism of BST-2 antagonism. Vpu-I17A, -A18F, -W22L, and -S23L co-localized with BST-2 within endosomal membranes while effectively enhancing virion release and down-regulating surface BST-2. However, Vpu-A18H was confined to an endoplasmic reticulum (ER)-like distribution, resulting in impaired down-regulation of BST-2 and reduced virion release. Brefeldin A confined wild type Vpu to the ER, resulting in a similarly impaired phenotype, as did the addition of a C-terminal ER-retention signal to Vpu. We determined the half-life of cell-surface BST-2 to be ~8 hours, whereas Vpu mediated an ~80% reduction of surface BST-2 within 6 hours, suggesting that TMD interactions between Vpu and BST-2 occur within post-ER membranes to directly and rapidly remove BST-2 from the cell surface and relieve restricted virion release.

Helical integrity and microsolvation of transmembrane domains from Flaviviridae envelope glycoproteins

Charged and polar amino acids in the transmembrane domains of integral membrane proteins can be crucial for protein function and also promote helix-helix association or protein oligomerization. Yet, our current understanding is still limited on how these hydrophilic amino acids are efficiently translocated from the Sec61/SecY translocon into the cell membrane during the biogenesis of membrane proteins. In hepatitis C virus, the putative transmembrane segments of envelope glycoproteins E1 and E2 were suggested to heterodimerize via a Lys-Asp ion-pair in the host endoplasmic reticulum. Therefore in this work, we carried out molecular dynamic simulations in explicit lipid bilayer and solvent environment to explore the stability of all possible bridging ion-pairs using the model of H-segment helix dimers. We observed that, frequently, several water molecules penetrated from the interface into the membrane core to stabilize the charged and polar pairs. The hydration time and amount of water molecules in the membrane core depended on the position of the charged residues as well as on the type of ion-pairs. Similar microsolvation events were observed in simulations of the putative E1-E2 transmembrane helix dimers. Simulations of helix monomers from other members of the Flaviviridae family suggest that these systems show similar behaviors. Thus this study illustrates the important contribution of water microsolvation to overcome the unfavorable energetic cost of burying charged and polar amino acids in membrane lipid bilayers. Also, it emphasizes the novel role of bridging charged or polar interactions stabilized by water molecules in the hydrophobic lipid bilayer core that has an important biological function for helix dimerization in several envelope glycoproteins from the family of Flaviviridae viruses.

Membrane raft association of the Vpu protein of human immunodeficiency virus type 1 correlates with enhanced virus release

The Vpu protein of human immunodeficiency virus type 1 (HIV-1) is known to enhance virion release from certain cell types. To accomplish this function, Vpu interacts with the restriction factor known as bone marrow stromal cell antigen 2 (BST-2)/tetherin. In this study, we analyzed whether the Vpu protein is associated with microdomains known as lipid or membrane rafts. Our results indicate that Vpu partially partitions into detergent-resistant membrane (DRM) fractions when expressed alone or in the context of simian-human immunodeficiency virus (SHIV) infection. The ability to be partitioned into rafts was observed with both subtype B and C Vpu proteins. The use of cholesterol lowering lovastatin/M-β-cyclodextrin and co-patching experiments confirmed that Vpu can be detected in cholesterol rich regions of membranes. Finally, we present data showing that raft association-defective transmembrane mutants of Vpu have impaired enhanced virus release function, but still maintain the ability to down-regulate CD4.

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

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.

Ligand-protein docking studies of potential HIV-1 drug compounds using the algorithm FlexX

Four compounds are docked to a pentameric bundle representing the transmembrane part of the Vpu protein from HIV-1. Employing the docking algorithm FlexX, their free energy of binding is estimated leading to the conclusion that potential drug candidates need to form H-bonds either with neighbouring or with n + 2 helices at the site of the serines within the bundle.

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.

Interaction and conformational dynamics of membrane-spanning protein helices

Within 1 or 2 decades, the reputation of membrane-spanning alpha-helices has changed dramatically. Once mostly regarded as dull membrane anchors, transmembrane domains are now recognized as major instigators of protein-protein interaction. These interactions may be of exquisite specificity in mediating assembly of stable membrane protein complexes from cognate subunits. Further, they can be reversible and regulatable by external factors to allow for dynamic changes of protein conformation in biological function. Finally, these helices are increasingly regarded as dynamic domains. These domains can move relative to each other in different functional protein conformations. In addition, small-scale backbone fluctuations may affect their function and their impact on surrounding lipid shells. Elucidating the ways by which these intricate structural features are encoded by the amino acid sequences will be a fascinating subject of research for years to come.

Structural implications of mutations assessed by molecular dynamics: Vpu1-32 from HIV-1

Structural pore models are generated for Vpu(1-32)WT from HIV-1 as well as for three mutants W23L, S24L and R31V. A computational methodology is employed which samples the whole conformational space of the pentameric assemblies of Vpu. The analysis of the related energy landscape reveals a small set of reasonable pore models, which are thoroughly investigated regarding their structural properties as well as their putative stability under native-like conditions. The models are also discussed in respect of earlier experimental findings about their channel activities. The study proposes functional pores reflecting the experimentally found conductance states of Vpu and its mutants.

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.