vpu transmembrane peptide structure obtained by site-specific fourier transform infrared dichroism and global molecular dynamics searching

Investigating Potassium Channels in Budding Yeast: A Genetic Sandbox

Like all species, the model eukaryote Saccharomyces cerevisiae, or Bakers' yeast, concentrates potassium in the cytosol as an electrogenic osmolyte and enzyme cofactor. Yeast are capable of robust growth on a wide variety of potassium concentrations, ranging from 10 µM to 2.5 M, due to the presence of a high-affinity potassium uptake system and a battery of cation exchange transporters. Genetic perturbation of either of these systems retards yeast growth on low or high potassium, respectively. However, these potassium-sensitized yeast are a powerful genetic tool, which has been leveraged for diverse studies. Notably, the potassium-sensitive cells can be transformed with plasmids encoding potassium channels from bacteria, plants, or mammals, and subsequent changes in growth rate have been found to correlate with the activity of the introduced potassium channel. Discoveries arising from the use of this assay over the past three decades have increased our understanding of the structure-function relationships of various potassium channels, the mechanisms underlying the regulation of potassium channel function and trafficking, and the chemical basis of potassium channel modulation. In this article, we provide an overview of the major genetic tools used to study potassium channels in S. cerevisiae, a survey of seminal studies utilizing these tools, and a prospective for the future use of this elegant genetic approach.

Modeling Interactions between Multicomponent Vesicles and Antimicrobial Peptide-Inspired Nanoparticles

We examine the interaction between peptide-inspired nanoparticles, or nanopins, and multicomponent vesicles using the dissipative particle dynamics simulation technique. We study the role of nanopin architecture and cholesterol concentration on the binding of the nanopins to the lipid bilayer, their insertion, and postembedding self-organization. We find the insertion to be triggered by enthalpically unfavorable interactions between the hydrophilic solvent and the lipophilic components of the nanopins. The nanopins are observed to form aggregates in solution, insert into the bilayer, and disassemble into the individual nanopins following the insertion process. We examine factors that influence the orientation of the nanopins in the host vesicle. We report the length of the hydrophilic segment of the nanopins to regulate their orientation within the clusters before the embedding process and in the bilayer, after the postinsertion disassembly of the aggregates. The orientation angle distribution for a given nanopin architecture is found to be driven by energy minimization. In addition, higher concentration of cholesterol is observed to constrain the orientation of the nanopins. We also report thermal fluctuations to induce transverse diffusion of nanopins with specific architectures. The incidence of transverse diffusion is observed to decrease with the concentration of cholesterol. Our results can provide guidelines for designing peptide-inspired nanoparticles or macromolecules that can interface with living cells to serve as sensors for applications in medicine, sustainability, and energy.

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.

Folding Simulations of the Transmembrane Helix of Virus Protein U in an Implicit Membrane Model

Vpu is an 81-amino-acid auxiliary membrane protein encoded by human immunodeficiency virus type 1 (HIV-1). One of its roles is to amplify viral release by self-assembling in homo-oligomers to form functional water-filled pores enabling the flux of ions across the membrane. Various NMR and CD studies have shown that the transmembrane domain of Vpu has a helical conformation. With a recently developed implicit membrane model and an efficient Monte Carlo (MC) algorithm using concerted backbone rotations, we simulate the folding of the transmembrane domain of Vpu at atomic resolution. The implicit membrane environment is based on the generalized Born theory and enables very long time scale events, such as folding to be observed using detailed all-atom representation of the protein. Such studies are currently computationally unfeasible with fully explicit lipid bilayer molecular dynamics simulations. The correct helical transmembrane structure of Vpu is predicted from extended conformations and remains stably inserted. Tilt and kink angles agree well with experimental estimates from NMR measurements. The experimentally observed change in tilt angle in membranes of varying hydrophobic width is accurately reproduced. The extensive simulation of a pentamer of the Vpu transmembrane domain in the implicit membrane gives results similar to the ones reported previously for fully explicit bilayer simulations.

Functional roles of HIV-1 Vpu and CD74: Details and implications of the Vpu-CD74 interaction

HIV-1 Vpu has a variety of functions, including CD4 degradation and the downregulation of MHCII. Downregulation of the MHCII occurs through Vpu binding to the cytoplasmic domain of CD74, the chaperone for antigen presentation. The CD74 cytoplasmic domain also plays a vital role in cell signaling through the activation of an NF-κB signal cascade for the maturation, proliferation and survival of B cells as well as by binding the macrophage inhibitory factor. In view of these functions, it follows that the Vpu-CD74 interaction has multiple downstream consequences for the immune system as it not only impairs foreign antigen presentation but may also have an effect on signal transduction cascades. It is thought that Vpu specifically targets intracellular CD74 while other HIV-1 proteins cannot. Therefore, this protein-protein interaction would be a potential drug target in order to reduce viral persistence. We review the functional importance and specific binding site of Vpu and CD74.

Vpu Protein: The Viroporin Encoded by HIV-1

Viral protein U (Vpu) is a lentiviral viroporin encoded by human immunodeficiency virus type 1 (HIV-1) and some simian immunodeficiency virus (SIV) strains. This small protein of 81 amino acids contains a single transmembrane domain that allows for supramolecular organization via homoligomerization or interaction with other proteins. The topology and trafficking of Vpu through subcellular compartments result in pleiotropic effects in host cells. Notwithstanding the high variability of its amino acid sequence, the functionality of Vpu is well conserved in pandemic virus isolates. This review outlines our current knowledge on the interactions of Vpu with the host cell. The regulation of cellular physiology by Vpu and the validity of this viroporin as a therapeutic target are also discussed.

Viroporins: structure, function and potential as antiviral targets

The channel-forming activity of a family of small, hydrophobic integral membrane proteins termed 'viroporins' is essential to the life cycles of an increasingly diverse range of RNA and DNA viruses, generating significant interest in targeting these proteins for antiviral development. Viroporins vary greatly in terms of their atomic structure and can perform multiple functions during the virus life cycle, including those distinct from their role as oligomeric membrane channels. Recent progress has seen an explosion in both the identification and understanding of many such proteins encoded by highly significant pathogens, yet the prototypic M2 proton channel of influenza A virus remains the only example of a viroporin with provenance as an antiviral drug target. This review attempts to summarize our current understanding of the channel-forming functions for key members of this growing family, including recent progress in structural studies and drug discovery research, as well as novel insights into the life cycles of many viruses revealed by a requirement for viroporin activity. Ultimately, given the successes of drugs targeting ion channels in other areas of medicine, unlocking the therapeutic potential of viroporins represents a valuable goal for many of the most significant viral challenges to human and animal health.

Prediction of the structures of helical membrane proteins based on a minimum unfavorable contacts approach

An understanding of structure-function relationships of membrane proteins continues to be a challenging problem, owing to the difficulty in obtaining their structures experimentally. This study suggests a method for modeling membrane protein structures that can be used to generate a reliable initial conformation prior to the use of other approaches for sampling conformations. It involves optimizing the orientation of hydrophilic residues so as to minimize unfavorable contacts with the hydrophobic tails of the lipid bilayer. Starting with the optimized initial conformation for three different proteins modeled based on this method, two independent approaches have been used for sampling the conformational space of the proteins. Both approaches are able to predict structures reasonably close to experimental structures, indicating that the initial structure enables the sampling of conformations that are close to the native structure. Possible improvements in the method for making it broadly applicable to helical membrane proteins are discussed.

Atomistic detailed mechanism and weak cation-conducting activity of HIV-1 Vpu revealed by free energy calculations

The viral protein U (Vpu) encoded by HIV-1 has been shown to assist in the detachment of virion particles from infected cells. Vpu forms cation-specific ion channels in host cells, and has been proposed as a potential drug target. An understanding of the mechanism of ion transport through Vpu is desirable, but remains limited because of the unavailability of an experimental structure of the channel. Using a structure of the pentameric form of Vpu – modeled and validated based on available experimental data – umbrella sampling molecular dynamics simulations (cumulative simulation time of more than 0.4 µs) were employed to elucidate the energetics and the molecular mechanism of ion transport in Vpu. Free energy profiles corresponding to the permeation of Na⁺ and K⁺ were found to be similar to each other indicating lack of ion selection, consistent with previous experimental studies. The Ser23 residue is shown to enhance ion transport via two mechanisms: creating a weak binding site, and increasing the effective hydrophilic length of the channel, both of which have previously been hypothesized in experiments. A two-dimensional free energy landscape has been computed to model multiple ion permeation, based on which a mechanism for ion conduction is proposed. It is shown that only one ion can pass through the channel at a time. This, along with a stretch of hydrophobic residues in the transmembrane domain of Vpu, explains the slow kinetics of ion conduction. The results are consistent with previous conductance studies that showed Vpu to be a weakly conducting ion channel.

Use of Isotope-Edited FTIR to Derive a Backbone Structure of a Transmembrane Protein

Solving structures of membrane proteins has always been a formidable challenge, yet even upon success, the results are normally obtained in a mimetic environment that can be substantially different from a biological membrane. Herein, we use noninvasive isotope-edited FTIR spectroscopy to derive a structural model for the SARS coronavirus E protein transmembrane domain in lipid bilayers. Molecular-dynamics-based structural refinement, incorporating the IR-derived orientational restraints points to the formation of a helical hairpin structure. Disulfide cross-linking and X-ray reflectivity depth profiling provide independent support of the results. The unusually short helical hairpin structure of the protein might explain its ability to deform bilayers and is reminiscent of other peptides with membrane disrupting functionalities. Taken together, we show that isotope-edited FTIR is a powerful tool to analyze small membrane proteins in their native environment, enabling us to relate the unusual structure of the SARS E protein to its function.

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.

Molecular dynamics simulations reveal the HIV-1 Vpu transmembrane protein to form stable pentamers

The human immunodeficiency virus type I (HIV-1) Vpu protein is 81 residues long and has two cytoplasmic and one transmembrane (TM) helical domains. The TM domain oligomerizes to form a monovalent cation selective ion channel and facilitates viral release from host cells. Exactly how many TM domains oligomerize to form the pore is still not understood, with experimental studies indicating the existence of a variety of oligomerization states. In this study, molecular dynamics (MD) simulations were performed to investigate the propensity of the Vpu TM domain to exist in tetrameric, pentameric, and hexameric forms. Starting with an idealized α-helical representation of the TM domain, a thorough search for the possible orientations of the monomer units within each oligomeric form was carried out using replica-exchange MD simulations in an implicit membrane environment. Extensive simulations in a fully hydrated lipid bilayer environment on representative structures obtained from the above approach showed the pentamer to be the most stable oligomeric state, with interhelical van der Waals interactions being critical for stability of the pentamer. Atomic details of the factors responsible for stable pentamer structures are presented. The structural features of the pentamer models are consistent with existing experimental information on the ion channel activity, existence of a kink around the Ile17, and the location of tetherin binding residues. Ser23 is proposed to play an important role in ion channel activity of Vpu and possibly in virus propagation.

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.

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.

Environment Polarity in Proteins Mapped Noninvasively by FTIR Spectroscopy

The polarity pattern of a macromolecule is of utmost importance to its structure and function. For example, one of the main driving forces for protein folding is the burial of hydrophobic residues. Yet polarity remains a difficult property to measure experimentally, due in part to its non-uniformity in the protein interior. Herein, we show that FTIR linewidth analysis of noninvasive 1-(13)C=(18)O labels can be used to obtain a reliable measure of the local polarity, even in a highly multi-phasic system, such as a membrane protein. We show that in the Influenza M2 H(+) channel, residues that line the pore are located in an environment that is as polar as fully solvated residues, while residues that face the lipid acyl chains are located in an apolar environment. Taken together, FTIR linewidth analysis is a powerful, yet chemically non-perturbing approach to examine one of the most important properties in proteins - polarity.

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.

Viral channel forming proteins - modeling the target

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

Oligomerization state and supramolecular structure of the HIV-1 Vpu protein transmembrane segment in phospholipid bilayers

HIV-1 Vpu is an 81-residue protein with a single N-terminal transmembrane (TM) helical segment that is involved in the release of new virions from host cell membranes. Vpu and its TM segment form ion channels in phospholipid bilayers, presumably by oligomerization of TM helices into a pore-like structure. We describe measurements that provide new constraints on the oligomerization state and supramolecular structure of residues 1-40 of Vpu (Vpu(1-40)), including analytical ultracentrifugation measurements to investigate oligomerization in detergent micelles, photo-induced crosslinking experiments to investigate oligomerization in bilayers, and solid-state nuclear magnetic resonance measurements to obtain constraints on intermolecular contacts between and orientations of TM helices in bilayers. From these data, we develop molecular models for Vpu TM oligomers. The data indicate that a variety of oligomers coexist in phospholipid bilayers, so that a unique supramolecular structure can not be defined. Nonetheless, since oligomers of various sizes have similar intermolecular contacts and orientations, molecular models developed from our data are most likely representative of Vpu TM oligomers that exist in host cell membranes.

HIV-1 Vpu targets cell surface markers CD4 and BST-2 through distinct mechanisms

Vpu is a small integral membrane protein encoded by HIV-1 and some SIV isolates. The protein is known to induce degradation of the viral receptor molecule CD4 and to enhance the release of newly formed virions from the cell surface. Vpu accomplishes these two functions through two distinct mechanisms. In the case of CD4, Vpu acts as a molecular adaptor to connect CD4 to an E3 ubiquitin ligase complex resulting in CD4 degradation by cellular proteasomes. This requires signals located in Vpu's cytoplasmic domain. Enhancement of virus release on the other hand involves the neutralization of a cellular host factor, BST-2 (also known as CD317, HM1.24, or tetherin) and requires Vpu's TM domain. The current review discusses recent advances on the role of Vpu in controlling degradation of CD4 and in regulating virus release.

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.

Molecular dynamics guided study of salt bridge length dependence in both fluorinated and non-fluorinated parallel dimeric coiled-coils

The alpha-helical coiled-coil is one of the most common oligomerization motifs found in both native and engineered proteins. To better understand the stability and dynamics of the coiled-coil motifs, including those modified by fluorination, several fluorinated and nonfluorinated parallel dimeric coiled-coil protein structures were designed and modeled. We also attempt to investigate how changing the length and geometry of the important stabilizing salt bridges influences the coiled-coil protein structure. Molecular dynamics (MD) and free energy simulations with AMBER used a particle mesh Ewald treatment of the electrostatics in explicit TIP3P solvent with balanced force field treatments. Preliminary studies with legacy force fields (ff94, ff96, and ff99) show a profound instability of the coiled-coil structures in short MD simulation. Significantly, better behavior is evident with the more balanced ff99SB and ff03 protein force fields. Overall, the results suggest that the coiled-coil structures can readily accommodate the larger acidic arginine or S-2,7-diaminoheptanedoic acid mutants in the salt bridge, whereas substitution of the smaller L-ornithine residue leads to rapid disruption of the coiled-coil structure on the MD simulation time scale. This structural distortion of the secondary structure allows both the formation of large hydration pockets proximal to the charged groups and within the hydrophobic core. Moreover, the increased structural fluctuations and movement lead to a decrease in the water occupancy lifetimes in the hydration pockets. In contrast, analysis of the hydration in the stable dimeric coiled-coils shows high occupancy water sites along the backbone residues with no water occupancy in the hydrophobic core, although transitory water interactions with the salt bridge residues are evident. The simulations of the fluorinated coiled-coils suggest that in some cases fluorination electrostatically stabilizes the intermolecular coiled-coil salt bridges. Structural analyses also reveal different side chain rotamer preferences for leucine when compared with 5,5,5,5',5',5'-hexafluoroleucine mutants. These observed differences in the side chain rotamer populations suggest differential changes in the side chain conformational entropy upon coiled-coil formation when the protein is fluorinated. The free energy of hydration of the isolated 5,5,5,5',5',5'-hexafluoroleucine amino acid is calculated to be 1.1 kcal/mol less stable than leucine; this hydrophobic penalty in the monomer may provide a driving force for coiled-coil dimer formation. Estimation of the ellipticity at 222 nm from a series of snapshots from the MD simulations with DicroCalc shows distinct increases in the ellipticity when the coiled-coil is fluorinated, which suggests that the helicity in the folded coiled-coils is greater when fluorinated.

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.

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

The dynamic behavior of monomeric Vpu(1-32) from HIV-1 in different lipid environments has been studied. The peptide shows highly flexible behavior during the simulations and easily adapts to changing lipid environments as it experiences when travelling through the Golgi apparatus. Protein-lipid interactions do not show any significant correlation towards lipid type or thickness based on multiple 10 ns simulations. The averaged structure of a series of 16 independent simulations suggest kink around Ser-24, which compensates the polarity of its side chain by forming hydrogen bonds with the carbonyl backbone of adjacent amino acids towards the N-terminus.

Application of solid-state NMR restraint potentials in membrane protein modeling

We have developed a set of orientational restraint potentials for solid-state NMR observables including (15)N chemical shift and (15)N-(1)H dipolar coupling. Torsion angle molecular dynamics simulations with available experimental (15)N chemical shift and (15)N-(1)H dipolar coupling as target values have been performed to determine orientational information of four membrane proteins and to model the structures of some of these systems in oligomer states. The results suggest that incorporation of the orientational restraint potentials into molecular dynamics provides an efficient means to the determination of structures that optimally satisfy the experimental observables without an extensive geometrical search.

Sequence-specific conformational flexibility of SNARE transmembrane helices probed by hydrogen/deuterium exchange

SNARE proteins mediate fusion of intracellular eukaryotic membranes and their alpha-helical transmembrane domains are known to contribute to lipid bilayer mixing. Synthetic transmembrane domain peptides were previously shown to mimic the function of SNARE proteins in that they trigger liposome fusion in a sequence-specific fashion. Here, we performed a detailed investigation of the conformational dynamics of the transmembrane helices of the presynaptic SNAREs synaptobrevin II and syntaxin 1a. To this end, we recorded deuterium/hydrogen-exchange kinetics in isotropic solution as well as in the membrane-embedded state. In solution, the exchange kinetics of each peptide can be described by three different classes of amide deuteriums that exchange with different rate constants. These are likely to originate from exchange at different domains of the helices. Interestingly, the rate constants of each class vary with the TMD sequence. Thus, the exchange rate is position-specific and sequence-specific. Further, the rate constants correlate with the previously determined membrane fusogenicities. In membranes, exchange is retarded and a significant proportion of amide hydrogens are protected from exchange. We conclude that the conformational dynamics of SNARE TMD helices is mechanistically linked to their ability to drive lipid mixing.

Transmembrane protein models based on high-throughput molecular dynamics simulations with experimental constraints

Elucidating the structure of transmembrane proteins domains with high-resolution methods is a difficult and sometimes impossible task. Here, we explain the method of combining a limited amount of experimental data with automated high-throughput molecular dynamics (MD) simulations of alpha-helical transmembrane bundles in an explicit lipid bilayer/water environment. The procedure uses a systematic conformational search of the helix rotation with experimentally constrained MDs simulations. The experimentally determined helix tilt and rotational angle of a labeled residue with site-specific infrared dichroism allows us to select a unique high-resolution model from a number of possible energy minima encountered in the systematic conformational search.

Phospholemman Transmembrane Structure Reveals Potential Interactions with Na+/K+-ATPase

Phospholemman (PLM) is a 72-residue bitopic cardiac transmembrane protein, which acts as a modulator of the Na(+)/K(+)-ATPase and the Na(+)/Ca(2+) exchanger and possibly forms taurine channels in nonheart tissue. This work presents a high resolution structural model obtained from a combination of site-specific infrared spectroscopy and experimentally constrained high throughput molecular dynamics (MD) simulations. Altogether, 37 experimental constraints, including nine long range orientational constraints, have been used during MD simulations in an explicit lipid bilayer/water system. The resulting tetrameric alpha-helical bundle has an average helix tilt of 7.3 degrees and a crossing angle close to 0 degrees . It does not reveal a hydrophilic pore, but instead strong interactions between various residues occlude any pore. The helix-helix packing is unusual, with Gly(19) and Gly(20) pointing to the outside of the helical bundle, facilitating potential interaction with other transmembrane proteins, thus providing a structural basis for the modulatory effect of PLM on the Na(+)/K(+)-ATPase. A two-stage model of interaction between PLM and the Na(+)/K(+)-ATPase is discussed involving PLM-ATPase interaction and subsequent formation of an unstable PLM trimer, which readily interacts with surrounding ATPase molecules. Further unconstrained MD simulations identified other packing models of PLM, one of which could potentially undergo a conformational transition to an open pore.

Coarse-grained molecular dynamics simulations of membrane proteins and peptides

Molecular dynamics (MD) simulations provide a valuable approach to the dynamics, structure, and stability of membrane-protein systems. Coarse-grained (CG) models, in which small groups of atoms are treated as single particles, enable extended (>100 ns) timescales to be addressed. In this study, we explore how CG-MD methods that have been developed for detergents and lipids may be extended to membrane proteins. In particular, CG-MD simulations of a number of membrane peptides and proteins are used to characterize their interactions with lipid bilayers. CG-MD is used to simulate the insertion of synthetic model membrane peptides (WALPs and LS3) into a lipid (PC) bilayer. WALP peptides insert in a transmembrane orientation, whilst the LS3 peptide adopts an interfacial location, both in agreement with experimental biophysical data. This approach is extended to a transmembrane fragment of the Vpu protein from HIV-1, and to the coat protein from fd phage. Again, simulated protein/membrane interactions are in good agreement with solid state NMR data for these proteins. CG-MD has also been applied to an M3-M4 fragment from the CFTR protein. Simulations of CFTR M3-M4 in a detergent micelle reveal formation of an alpha-helical hairpin, consistent with a variety of biophysical data. In an I231D mutant, the M3-M4 hairpin is additionally stabilized via an inter-helix Q207/D231 interaction. Finally, CG-MD simulations are extended to a more complex membrane protein, the bacterial sugar transporter LacY. Comparison of a 200 ns CG-MD simulation of LacY in a DPPC bilayer with a 50 ns atomistic simulation of the same protein in a DMPC bilayer shows that the two methods yield comparable predictions of lipid-protein interactions. Taken together, these results demonstrate the utility of CG-MD simulations for studies of membrane/protein interactions.

Structure, dynamics and topology of membrane polypeptides by oriented 2H solid-state NMR spectroscopy

Knowledge of the structure, dynamics and interactions of polypeptides when associated with phospholipid bilayers is key to understanding the functional mechanisms of channels, antibiotics, signal- or translocation peptides. Solid-state NMR spectroscopy on samples uniaxially aligned relative to the magnetic field direction offers means to determine the alignment of polypeptide bonds and domains relative to the bilayer normal. Using this approach the (15)N chemical shift of amide bonds provides a direct indicator of the approximate helical tilt, whereas the (2)H solid-state NMR spectra acquired from peptides labelled with 3,3,3-(2)H(3)-alanines contain valuable complimentary information for a more accurate analysis of tilt and rotation pitch angles. The deuterium NMR line shapes are highly sensitive to small variations in the alignment of the C(alpha)-C(beta) bond relative to the magnetic field direction and, therefore, also the orientational distribution of helices relative to the membrane normal. When the oriented membrane samples are investigated with their normal perpendicular to the magnetic field direction, the rate of rotational diffusion can be determined in a semi-quantitative manner and thereby the aggregation state of the peptides can be analysed. Here the deuterium NMR approach is first introduced showing results from model amphipathic helices. Thereafter investigations of the viral channel peptides Vpu(1-27) and Influenza A M2(22-46) are shown. Whereas the (15)N chemical shift data confirm the transmembrane helix alignments of these hydrophobic sequences, the deuterium spectra indicate considerable mosaic spread in the helix orientations. At least two peptide populations with differing rotational correlation times are apparent in the deuterium spectra of the viral channels suggesting an equilibrium between monomeric peptides and oligomeric channel configurations under conditions where solid-state NMR structural studies of these peptides have previously been performed.

Evidence for the formation of a heptameric ion channel complex by the hepatitis C virus p7 protein in vitro

The p7 protein of hepatitis C virus functions as an ion channel both in vitro and in cell-based assays and is inhibited by amantadine, long alkyl chain imino-sugar derivatives, and amiloride compounds. Future drug design will be greatly aided by information on the stoichiometry and high resolution structure of p7 ion channel complexes. Here, we have refined a bacterial expression system for p7 based on a glutathione S-transferase fusion methodology that circumvents the inherent problems of hydrophobic protein purification and the limitations of chemical synthesis. Rotational averaging and harmonic analysis of transmission electron micrographs of glutathione S-transferase-FLAG-p7 fusion proteins in liposomes revealed a heptameric stoichiometry. The oligomerization of p7 protein was then confirmed by SDS-PAGE and mass spectrometry analysis of pure, concentrated FLAG-p7. The same protein was also confirmed to function as an ion channel in suspended lipid bilayers and was inhibited by amantadine. These data validate this system as a means of generating high resolution structural information on the p7 ion channel complex.

Viral ion channel proteins in model membranes: a comparative study by X-ray reflectivity

We have investigated the effect of the transmembrane domain of three viral ion channel proteins on the lipid bilayer structure by X-ray reflectivity and scattering from oriented planar bilayers. The proteins show a similar effect on the lipid bilayer structural parameters: an increase in the lipid bilayer hydrophobic core, a decrease in the amplitude of the vertical density profile and a systematic change in the ordering of the acyl chains as a function of protein-to-lipid ratio. These results are discussed in a comparative view.

Isotope-edited IR spectroscopy for the study of membrane proteins

Fourier transform infrared (FTIR) spectroscopy has long been a powerful tool for structural analysis of membrane proteins. However, because of difficulties in resolving contributions from individual residues, most of the derived measurements tend to yield average properties for the system under study. Isotope editing, through its ability to resolve individual vibrations, establishes FTIR as a method that is capable of yielding accurate structural data on individual sites in a protein.

The transmembrane domain of the oncogenic mutant ErbB-2 receptor: a structure obtained from site-specific infrared dichroism and molecular dynamics

ErbB-2 is a member of the family of epidermal growth factor receptors, which shows an oncogenic mutation in the rat gene neu, Val664Glu in the transmembrane domain that causes permanent dimerisation and subsequently leads to uncontrollable cell division and tumour formation. We have obtained the alpha-helical structure of the mutant transmembrane domain dimer experimentally with site-specific infrared dichroism (SSID) based on six transmembrane peptides with 13C18O carbonyl group-labelled residues. The derived orientational data indicate a local helix tilt ranging from 28(+/-6) degrees to 22(+/-4) degrees. Altogether using orientational constraints from SSID and experimental alpha-helical constraints while performing a systematic conformational search including molecular dynamics simulation in a lipid bilayer, we have obtained a unique experimentally defined atomic structure. The resulting structure consists of a right handed alpha-helical bundle with the residues Ile659, Val663, Leu667, Ile671, Val674 and Leu679 in the dimerisation interface. The right-handed bundle is in contrast to the left-handed structures obtained in previous modelling efforts. In order to facilitate tight helical packing, the spacious Glu664 residues do not interact directly but with water molecules that enter the bilayer.

Three-dimensional structure of the transmembrane domain of Vpu from HIV-1 in aligned phospholipid bicelles

The three-dimensional backbone structure of the transmembrane domain of Vpu from HIV-1 was determined by solid-state NMR spectroscopy in two magnetically-aligned phospholipid bilayer environments (bicelles) that differed in their hydrophobic thickness. Isotopically labeled samples of Vpu(2-30+), a 36-residue polypeptide containing residues 2-30 from the N-terminus of Vpu, were incorporated into large (q = 3.2 or 3.0) phospholipid bicelles composed of long-chain ether-linked lipids (14-O-PC or 16-O-PC) and short-chain lipids (6-O-PC). The protein-containing bicelles are aligned in the static magnetic field of the NMR spectrometer. Wheel-like patterns of resonances characteristic of tilted transmembrane helices were observed in two-dimensional (1)H/(15)N PISEMA spectra of uniformly (15)N-labeled Vpu(2-30+) obtained on bicelle samples with their bilayer normals aligned perpendicular or parallel to the direction of the magnetic field. The NMR experiments were performed at a (1)H resonance frequency of 900 MHz, and this resulted in improved data compared to lower-resonance frequencies. Analysis of the polarity-index slant-angle wheels and dipolar waves demonstrates the presence of a transmembrane alpha-helix spanning residues 8-25 in both 14-O-PC and 16-O-PC bicelles, which is consistent with results obtained previously in micelles by solution NMR and mechanically aligned lipid bilayers by solid-state NMR. The three-dimensional backbone structures were obtained by structural fitting to the orientation-dependent (15)N chemical shift and (1)H-(15)N dipolar coupling frequencies. Tilt angles of 30 degrees and 21 degrees are observed in 14-O-PC and 16-O-PC bicelles, respectively, which are consistent with the values previously determined for the same polypeptide in mechanically-aligned DMPC and DOPC bilayers. The difference in tilt angle in C14 and C16 bilayer environments is also consistent with previous results indicating that the transmembrane helix of Vpu responds to hydrophobic mismatch by changing its tilt angle. The kink found in the middle of the helix in the longer-chain C18 bilayers aligned on glass plates was not found in either of these shorter-chain (C14 or C16) bilayers.

Model of a putative pore: the pentameric alpha-helical bundle of SARS coronavirus E protein in lipid bilayers

The coronavirus responsible for the severe acute respiratory syndrome contains a small envelope protein, E, with putative involvement in host apoptosis and virus morphogenesis. To perform these functions, it has been suggested that protein E can form a membrane destabilizing transmembrane (TM) hairpin, or homooligomerize to form a TM pore. Indeed, in a recent study we reported that the α-helical putative transmembrane domain of E protein (ETM) forms several SDS-resistant TM interactions: a dimer, a trimer, and two pentameric forms. Further, these interactions were found to be evolutionarily conserved. Herein, we have studied multiple isotopically labeled ETM peptides reconstituted in model lipid bilayers, using the orientational parameters derived from infrared dichroic data. We show that the topology of ETM is consistent with a regular TM α-helix. Further, the orientational parameters obtained unequivocally correspond to a homopentameric model, by comparison with previous predictions. We have independently confirmed that the full polypeptide of E protein can also aggregate as pentamers after expression in Escherichia coli. This interaction must be stabilized, at least partially, at the TM domain. The model we report for this pentameric α-helical bundle may explain some of the permabilizing properties of protein E, and should be the basis of mutagenesis efforts in future functional studies.

Full Length Vpu from HIV-1: Combining Molecular Dynamics Simulations with NMR Spectroscopy

Based on structures made available by solution NMR, molecular models of the protein Vpu from HIV-1 were built and refined by 6 ns MD simulations in a fully hydrated lipid bilayer. Vpu is an 81 amino acid type I integral membrane protein encoded by the human immunodeficiency virus type-1 (HIV-1) and closely related simian immunodeficiency viruses (SIVs). Its role is to amplify viral release. Upon phosphorylation, the cytoplasmic domain adopts a more compact shape with helices 2 and 3 becoming almost parallel to each other. A loss of helicity for several residues belonging to the helices adjacent to both ends of the loop region containing serines 53 and 57 is observed. A fourth helix, present in one of the NMR-based structures of the cytoplasmic domain and located near the C-terminus, is lost upon phosphorylation.

Energetics of the native and non-native states of the glycophorin transmembrane helix dimer

Using an implicit membrane model (IMM1), we examine whether the structure of the transmembrane domain of Glycophorin A (GpA) could be predicted based on energetic considerations alone. The energetics of native GpA shows that van der Waals interactions make the largest contribution to stability. Although specific electrostatic interactions are stabilizing, the overall electrostatic contribution is close to zero. The GXXXG motif contributes significantly to stability, but residues outside this motif contribute almost twice as much. To generate non-native states a global conformational search was done on two segments of GpA: an 18-residue peptide (GpA74-91) that is embedded in the membrane and a 29-residue peptide (GpA70-98) that has additional polar residues flanking the transmembrane region. Simulated annealing was done on a large number of conformations generated from parallel, antiparallel, left- and right-handed starting structures by rotating each helix at 20 degrees intervals around its helical axis. Several crossing points along the helix dimer were considered. For 18-residue parallel topology, an ensemble of native-like structures was found at the lowest effective energy region; the effective energy is lowest for a right-handed structure with an RMSD of 1.0 A from the solid-state NMR structure with correct orientation of the helices. For the 29-residue peptide, the effective energies of several left-handed structures were lower than that of the native, right-handed structure. This could be due to deficiencies in modeling the interactions between charged sidechains and/or omission of the sidechain entropy contribution to the free energy. For 18-residue antiparallel topology, both IMM1 and a Generalized Born model give effective energies that are lower than that of the native structure. In contrast, the Poisson-Boltzmann solvation model gives lower effective energy for the parallel topology, largely because the electrostatic solvation energy is more favorable for the parallel structure. IMM1 seems to underestimate the solvation free energy advantage when the CO and NH dipoles just outside the membrane are parallel. This highlights the importance of electrostatic interactions even when these are not obvious by looking at the structures.

Protein-protein interactions: modeling the hepatitis C virus ion channel p7

The p7 protein is a small ion-channel-forming membrane polypeptide encoded by the hepatitis C virus which consists of two transmembrane alpha-helices, TM1 and TM2, and can be blocked by long-alkyl-chain iminosugar derivatives. The length of TM1 and TM2 was estimated by employing different secondary structure prediction algorithms and is proposed to span from Ala-10 to Leu-32 for TM1 and from Trp-36 to Pro-58 for TM2. A configurational search protocol based on simulated annealing combined with short restrained molecular dynamics simulations is used in addition to protein-protein docking to investigate the packing of TM1/TM2. Full p7 oligomeric bundles were generated, and in the most plausible models serines and threonines are facing the hydrophilic pore. In these models, His-17 would be a pore-facing residue, suggesting that p7 may be sensitive to pH in respect to its function.

Structure and Dynamics of the HIV-1 Vpu Transmembrane Domain Revealed by Solid-State NMR with Magic-Angle Spinning†

We report solid-state nuclear magnetic resonance (NMR) measurements on the peptide Vpu(1-40), comprising residues 1-40 of the 81-residue type 1 integral membrane protein Vpu encoded by the HIV-1 genome. On the basis of a combination of 13C and 15N NMR chemical shifts under magic-angle spinning (MAS), effects of local mobility on NMR signal intensities, site-specific MAS NMR line widths, and NMR-detected hydrogen-deuterium exchange, we develop a model for the structure and dynamics of the Vpu(1-40) monomer in phospholipid bilayer membranes. Our data are largely consistent with earlier structural studies of Vpu peptides by Opella and co-workers, in which solution NMR and solid-state NMR without MAS were used, but our data provide new information about local variations in the degree of mobility and structural order. In addition, our data indicate that the transmembrane alpha-helix of Vpu(1-40) extends beyond the hydrophobic core of the bilayer. We find no evidence for heterogeneity in the conformation and intermolecular contacts of the transmembrane alpha-helix, with the exception of two distinct chemical shifts observed for the C alpha and C beta atoms of A18 that may reflect distinct modes of helix-helix interaction. These results have possible implications for the supramolecular structure of Vpu oligomers that form cation-selective ion channels.

Evaluating tilt angles of membrane-associated helices: comparison of computational and NMR techniques

A computational method to calculate the orientation of membrane-associated alpha-helices with respect to a lipid bilayer has been developed. It is based on a previously derived implicit membrane representation, which was parameterized using the structures of 46 alpha-helical membrane proteins. The method is validated by comparison with an independent data set of six transmembrane and nine antimicrobial peptides of known structure and orientation. The minimum energy orientations of the transmembrane helices were found to be in good agreement with tilt and rotation angles known from solid-state NMR experiments. Analysis of the free-energy landscape found two types of minima for transmembrane peptides: i), Surface-bound configurations with the helix long axis parallel to the membrane, and ii), inserted configurations with the helix spanning the membrane in a perpendicular orientation. In all cases the inserted configuration also contained the global energy minimum. Repeating the calculations with a set of solution NMR structures showed that the membrane model correctly distinguishes native transmembrane from nonnative conformers. All antimicrobial peptides investigated were found to orient parallel and bind to the membrane surface, in agreement with experimental data. In all cases insertion into the membrane entailed a significant free-energy penalty. An analysis of the contributions of the individual residue types confirmed that hydrophobic residues are the main driving force behind membrane protein insertion, whereas polar, charged, and aromatic residues were found to be important for the correct orientation of the helix inside the membrane.