Structural and dynamic effects of cholesterol at preferred sites of interaction with rhodopsin identified from microsecond length molecular dynamics simulations

An unresolved question about GPCR function is the role of membrane components in receptor stability and activation. In particular, cholesterol is known to affect the function of membrane proteins, but the details of its effect on GPCRs are still elusive. Here, we describe how cholesterol modulates the behavior of the TM1‐TM2‐TM7‐helix 8(H8) functional network that comprises the highly conserved NPxxY(x)_5,6F motif, through specific interactions with the receptor. The inferences are based on the analysis of microsecond length molecular dynamics (MD) simulations of rhodopsin in an explicit membrane environment. Three regions on the rhodopsin exhibit the highest cholesterol density throughout the trajectory: the extracellular end of TM7, a location resembling the high‐density sterol area from the electron microscopy data; the intracellular parts of TM1, TM2, and TM4, a region suggested as the cholesterol binding site in the recent X‐ray crystallography data on β₂‐adrenergic GPCR; and the intracellular ends of TM2‐TM3, a location that was categorized as the high cholesterol density area in multiple independent 100 ns MD simulations of the same system. We found that cholesterol primarily affects specific local perturbations of the helical TM domains such as the kinks in TM1, TM2, and TM7. These local distortions, in turn, relate to rigid‐body motions of the TMs in the TM1‐TM2‐TM7‐H8 bundle. The specificity of the effects stems from the nonuniform distribution of cholesterol around the protein. Through correlation analysis we connect local effects of cholesterol on structural perturbations with a regulatory role of cholesterol in the structural rearrangements involved in GPCR function. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Molecular Dynamics Simulations of Homo-oligomeric Bundles Embedded Within a Lipid Bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)₃-NH₂ (LS2) or Ac-(LSSLLSL)₃-NH₂ (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine–serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a “dimer of dimers.” In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.     

NMR determination of protein partitioning into membrane domains with different curvatures and application to the influenza M2 peptide

The M2 protein of the influenza A virus acts both as a drug-sensitive proton channel and mediates virus budding through membrane scission. The segment responsible for causing membrane curvature is an amphipathic helix in the cytoplasmic domain of the protein. Here, we use ³¹P and ¹³C solid-state NMR to examine M2-induced membrane curvature. M2(22–46), which includes only the transmembrane (TM) helix, and M2(21–61), which contains an additional amphipathic helix, are studied. ³¹P chemical shift lineshapes indicate that M2(21–61) causes a high-curvature isotropic phase to both cholesterol-rich virus-mimetic membranes and 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayers, whereas M2(22–46) has minimal effect. The lamellar and isotropic domains have distinct ³¹P isotropic chemical shifts, indicating perturbation of the lipid headgroup conformation by the amphipathic helix. ³¹P- and ¹³C-detected ¹H T₂ relaxation and two-dimensional peptide-lipid correlation spectra show that M2(21–61) preferentially binds to the high-curvature domain. ³¹P linewidths indicate that the isotropic vesicles induced by M2(21–61) are 10–35 nm in diameter, and the virus-mimetic vesicles are smaller than the 1,2-dimyristoyl-sn-glycero-3-phosphocholine vesicles. A strong correlation is found between high membrane curvature and weak drug-binding ability of the TM helix. Thus, the M2 amphipathic helix causes membrane curvature, which in turn perturbs the TM helix conformation, abolishing drug binding. These NMR experiments are applicable to other curvature-inducing membrane proteins such as fusion proteins and antimicrobial peptides.     

Stabilization of the human beta2-adrenergic receptor TM4-TM3-TM5 helix interface by mutagenesis of Glu122(3.41), a critical residue in GPCR structure

G protein-coupled receptor (GPCR) instability represents one of the most profound obstacles in the structural study of GPCRs that bind diffusible ligands. The introduction of targeted mutations at nonconserved residues that lie proximal to helix interfaces has the potential to enhance the fold stability of the receptor helix bundle while maintaining wild-type receptor function. To test this hypothesis, we studied the effect of amino acid substitutions at Glu122^3.41 in the well-studied β₂-adrenergic receptor (β₂AR), which was predicted from sequence conservation to lie at a position equivalent to a tryptophan residue in rhodopsin at the 3,4,5 helix interface among transmembrane (TM) domains 3, 4, and 5. Replacement of Glu122^3.41 with bulky hydrophobic residues, such as tryptophan, tyrosine, and phenylalanine, increases the yield of functionally folded β₂AR by as much as 5-fold. Receptor stability in detergent solution was studied by isothermal denaturation, and it was found that the E122W and E122Y mutations enhanced the β₂AR thermal half-life by 9.3- and 6.7-fold, respectively, at 37 °C. The β₁AR was also stabilized by the introduction of tryptophan at Glu147^3.41, and the effect on protein behavior was similar to the rescue of the unstable wild-type receptor by the antagonist propranolol. Molecular modeling of the E122W and E122Y mutants revealed that the tryptophan ring edge and tyrosine hydroxyl are positioned proximal to the helical break in TM5 introduced by the conserved Pro211^5.50 and may stabilize the helix by interacting favorably with the unpaired carbonyl oxygen of Val206^5.45. Conformational flexibility of TM5 is likely to be a general property of class A GPCRs; therefore, engineering of the TM4-TM3-TM5 interface at the 3.41 position may provide a general strategy for the stabilization of other receptors.     

Molecular Dynamics Simulations of Homo-oligomeric Bundles Embedded Within a Lipid Bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)₃-NH₂ (LS2) or Ac-(LSSLLSL)₃-NH₂ (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine–serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a “dimer of dimers.” In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.     

Simulation of the packing of idealized transmembrane α-helix bundles

The aim of this study is to investigate if the packing motifs of native transmembrane helices can be produced by simulations with simple potentials and to develop a method for the rapid generation of initial candidate models for integral membrane proteins composed of bundles of transmembrane helices. Constituent residues are mapped along the helix axis in order to maintain the amino acid sequence-dependent properties of the helix. Helix packing is optimized according to a semi-empirical potential mainly composed of four components: a bilayer potential, a crossing angle potential, a helix dipole potential and a helix-helix distance potential. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. Preliminary testing of the method has been conducted with idealized seven Ala₂₀ helix bundles. The structures generated show a high degree of compactness. It was observed that both bacteriorhodopsin-like and δ-endotoxin-like structures are generated in seven-helix bundle simulations, within which the composition varies dependent upon the cooling rate. The simulation method has also been employed to explore the packing of N = 4 and N = 12 transmembrane helix bundles. The results suggest that seven and 12 transmembrane helix bundles resembling those observed experimentally (e.g., bacteriorhodopsin, rhodopsin and cytochrome c oxidase subunit I) may be generated by simulations using simple potentials. Received: 16 November 1998 / Revised version: 26 March 1999 / Accepted: 8 April 1999     

Microbial cell disruption for improving lipid recovery using pressurized CO2: Role of CO2 solubility in cell suspension,sugar broth,and spent media

The study of in situ gas explosion to lyse the triglyceride‐rich cells involves the solubilization of gas (e.g., carbon dioxide, CO₂) in lipid‐rich cells under pressure followed by a rapid decompression, which allows the gas inside the cell to rapidly expand and rupture the cell from inside out. The aim of this study was to perform the cell disruption using pressurized CO₂ as well as to determine the solubility of CO₂ in Rhodotorula glutinis cell suspension, sugar broth media, and spent media. Cell disruption of R. glutinis was performed at two pressures of 2,000 and 3,500 kPa, respectively, at 295.2 K, and it was found from both scanning electron microscopy (SEM) and plate count that a substantial amount of R. glutinis was disrupted due to the pressurized CO₂. We also found a considerable portion of lipid present in the aqueous phase after the disruption at P = 3,500 kPa compared to control (no pressure) and P = 2,000 kPa, which implied that more intracellular lipid was released due to the pressurized CO₂. Solubility of CO₂ in R. glutinis cell suspension was found to be higher than the solubility of CO₂ in both sugar broth media and spent media. Experimental solubility was correlated using the extended Henry's law, which showed a good agreement with the experimental data. Enthalpy and entropy of dissolution of CO₂ were found to be ?14.22 kJ mol^?1 and 48.10 kJ mol^?1 K^?1, 9.64 kJ mol^?1 and 32.52 kJ mol^?1 K^?1, and 7.50 kJ mol^?1 and 25.22 kJ mol^?1 K^?1 in R. glutinis, spent media, and sugar broth media, respectively. © 2017 American Institute of Chemical Engineers Biotechnol. Prog., 33:737–748, 2017     

Alamethicin Aggregation in Lipid Membranes

X-ray scattering features induced by aggregates of alamethicin (Alm) were obtained in oriented stacks of model membranes of DOPC(diC18:1PC) and diC22:1PC. The first feature obtained near full hydration was Bragg rod in-plane scattering near 0.11 ?⁻¹ in DOPC and near 0.08 ?⁻¹ in diC22:1PC at a 1:10 Alm:lipid ratio. This feature is interpreted as bundles consisting of n Alm monomers in a barrel-stave configuration surrounding a water pore. Fitting the scattering data to previously published molecular dynamics simulations indicates that the number of peptides per bundle is n = 6 in DOPC and n ≥ 9 in diC22:1PC. The larger bundle size in diC22:1PC is explained by hydrophobic mismatch of Alm with the thicker bilayer. A second diffuse scattering peak located at q _r ≈ 0.7 ?⁻¹ is obtained for both DOPC and diC22:1PC at several peptide concentrations. Theoretical calculations indicate that this peak cannot be caused by the Alm bundle structure. Instead, we interpret it as being due to two-dimensional hexagonally packed clusters in equilibrium with Alm bundles. As the relative humidity was reduced, interactions between Alm in neighboring bilayers produced more peaks with three-dimensional crystallographic character that do not index with the conventional hexagonal space groups.     

The pore-lining regions in cytochrome c oxidases: A computational analysis of caveolin,cholesterol and transmembrane helix contributions to proton movement

Cytochrome c oxidase (CcO) is the terminal enzyme in the electron transfer chain. CcO catalyzes a four electron reduction of O₂ to water at a catalytic site formed by high-spin heme (a₃) and copper atoms (Cu_B). While it is recognized that proton movement is coupled to oxygen reduction, the proton channel(s) have not been well defined. Using computational methods developed to study protein topology, membrane channels and 3D packing arrangements within transmembrane (TM) helix arrays, we find that subunit-1 (COX-1), subunit-2 (COX-2) and subunit-3 (COX-3) contribute 139, 46 and 25 residues, respectively, to channel formation between the mitochondrial matrix and intermembrane space. Nine of 12 TM helices in COX-1, both helices in COX-2 and 5 of the 6 TM helices in COX-3 are pore-lining regions (possible channel formers). Heme a₃ and the Cu_B sites (as well as the Cu_A center of COX-2) are located within the channel that includes TM-6, TM-7, TM-10 and TM-11 of COX-1 and are associated with multiple cholesterol and caveolin-binding (CB) motifs. Sequence analysis identifies five CB motifs within COX-1, two within COX-2 and four within COX-3; each caveolin containing a pore-lining helix C-terminal to a TM helix–turn–helix. Channel formation involves interaction between multiple pore-lining regions within protein subunits and/or dimers. PoreWalker analysis lends support to the D-channel model of proton translocation. Under physiological conditions, caveolins may introduce channel formers juxtaposed to those in COX-1, COX-2 and COX-3, which together with cholesterol may form channel(s) essential for proton translocation through the inner mitochondrial membrane.     

Self-assembly of a simple membrane protein: coarse-grained molecular dynamics simulations of the influenza M2 channel

The transmembrane (TM) domain of the M2 channel protein from influenza A is a homotetrameric bundle of α-helices and provides a model system for computational approaches to self-assembly of membrane proteins. Coarse-grained molecular dynamics (CG-MD) simulations have been used to explore partitioning into a membrane of M2 TM helices during bilayer self-assembly from lipids. CG-MD is also used to explore tetramerization of preinserted M2 TM helices. The M2 helix monomer adopts a membrane spanning orientation in a lipid (DPPC) bilayer. Multiple extended CG-MD simulations (5 × 5 μs) were used to study the tetramerization of inserted M2 helices. The resultant tetramers were evaluated in terms of the most populated conformations and the dynamics of their interconversion. This analysis reveals that the M2 tetramer has 2× rotationally symmetrical packing of the helices. The helices form a left-handed bundle, with a helix tilt angle of ∼16°. The M2 helix bundle generated by CG-MD was converted to an atomistic model. Simulations of this model reveal that the bundle's stability depends on the assumed protonation state of the H37 side chains. These simulations alongside comparison with recent x-ray (3BKD) and NMR (2RLF) structures of the M2 bundle suggest that the model yielded by CG-MD may correspond to a closed state of the channel.     

Structure of a Double Transmembrane Fragment of a G-Protein-Coupled Receptor in Micelles

The structure and dynamic properties of an 80-residue fragment of Ste2p, the G-protein-coupled receptor for α-factor of Saccharomyces cerevisiae, was studied in LPPG micelles with the use of solution NMR spectroscopy. The fragment Ste2p(G31-T110) (TM1-TM2) consisted of 19 residues from the N-terminal domain, the first TM helix (TM1), the first cytoplasmic loop, the second TM helix (TM2), and seven residues from the first extracellular loop. Multidimensional NMR experiments on [¹⁵N], [¹⁵N, ¹³C], [¹⁵N, ¹³C, ²H]-labeled TM1-TM2 and on protein fragments selectively labeled at specific amino acid residues or protonated at selected methyl groups resulted in >95% assignment of backbone and side-chain nuclei. The NMR investigation revealed the secondary structure of specific residues of TM1-TM2. TALOS constraints and NOE connectivities were used to calculate a structure for TM1-TM2 that was highlighted by the presence of three α-helices encompassing residues 39-47, 49-72, and 80-103, with higher flexibility around the internal Arg⁵⁸ site of TM1. RMSD values of individually superimposed helical segments 39-47, 49-72, and 80-103 were 0.25 ± 0.10 Å, 0.40 ± 0.13 Å, and 0.57 ± 0.19 Å, respectively. Several long-range interhelical connectivities supported the folding of TM1-TM2 into a tertiary structure typified by a crossed helix that splays apart toward the extracellular regions and contains considerable flexibility in the G⁵⁶VRSG⁶⁰ region. ¹⁵N-relaxation and hydrogen-deuterium exchange data support a stable fold for the TM parts of TM1-TM2, whereas the solvent-exposed segments are more flexible. The NMR structure is consistent with the results of biochemical experiments that identified the ligand-binding site within this region of the receptor.     

Roles of Subunit NuoK (ND4L) in the Energy-transducing Mechanism of Escherichia coli NDH-1 (NADH:Quinone Oxidoreductase)

The bacterial H⁺-translocating NADH:quinone oxidoreductase (NDH-1) catalyzes electron transfer from NADH to quinone coupled with proton pumping across the cytoplasmic membrane. The NuoK subunit (counterpart of the mitochondrial ND4L subunit) is one of the seven hydrophobic subunits in the membrane domain and bears three transmembrane segments (TM1–3). Two glutamic residues located in the adjacent transmembrane helices of NuoK are important for the energy coupled activity of NDH-1. In particular, mutation of the highly conserved carboxyl residue (_KGlu-36 in TM2) to Ala led to a complete loss of the NDH-1 activities. Mutation of the second conserved carboxyl residue (_KGlu-72 in TM3) moderately reduced the activities. To clarify the contribution of NuoK to the mechanism of proton translocation, we relocated these two conserved residues. When we shifted _KGlu-36 along TM2 to positions 32, 38, 39, and 40, the mutants largely retained energy transducing NDH-1 activities. According to the recent structural information, these positions are located in the vicinity of _KGlu-36, present in the same helix phase, in an immediately before and after helix turn. In an earlier study, a double mutation of two arginine residues located in a short cytoplasmic loop between TM1 and TM2 (loop-1) showed a drastic effect on energy transducing activities. Therefore, the importance of this cytosolic loop of NuoK (_KArg-25, _KArg-26, and _KAsn-27) for the energy transducing activities was extensively studied. The probable roles of subunit NuoK in the energy transducing mechanism of NDH-1 are discussed.     

Anesthetic effects on the structure and dynamics of the second transmembrane domains of nAChR α4β2

Channel functions of the neuronal α4β2 nicotinic acetylcholine receptor (nAChR), one of the most widely expressed subtypes in the brain, can be inhibited by volatile anesthetics. Our Na⁺ flux experiments confirmed that the second transmembrane domains (TM2) of α4 and β2 in 2:3 stoichiometry, (α4)₂(β2)₃, could form pentameric channels, whereas the α4 TM2 alone could not. The structure, topology, and dynamics of the α4 TM2 and (α4)₂(β2)₃ TM2 in magnetically aligned phospholipid bicelles were investigated using solid-state NMR spectroscopy in the absence and presence of halothane and isoflurane, two clinically used volatile anesthetics. ²H NMR demonstrated that anesthetics increased lipid conformational heterogeneity. Such anesthetic effects on lipids became more profound in the presence of transmembrane proteins. PISEMA experiments on the selectively ¹⁵N-labeled α4 TM2 showed that the TM2 formed transmembrane helices with tilt angles of 12° ± 1° and 16° ± 1° relative to the bicelle normal for the α4 and (α4)₂(β2)₃ samples, respectively. Anesthetics changed the tilt angle of the α4 TM2 from 12° ± 1° to 14° ± 1°, but had only a subtle effect on the tilt angle of the (α4)₂(β2)₃ TM2. A small degree of wobbling motion of the helix axis occurred in the (α4)₂(β2)₃ TM2. In addition, a subset of the (α4)₂(β2)₃ TM2 exhibited counterclockwise rotational motion around the helix axis on a time scale slower than 10^- 4 s in the presence of anesthetics. Both helical tilting and rotational motions have been identified computationally as critical elements for ion channel functions. This study suggested that anesthetics could alter these motions to modulate channel functions.     

Assembly of transmembrane helices of simple polytopic membrane proteins from sequence conservation patterns

The transmembrane (TM) domains of most membrane proteins consist of helix bundles. The seemingly simple task of TM helix bundle assembly has turned out to be extremely difficult. This is true even for simple TM helix bundle proteins, i.e., those that have the simple form of compact TM helix bundles. Herein, we present a computational method that is capable of generating native-like structural models for simple TM helix bundle proteins having modest numbers of TM helices based on sequence conservation patterns. Thus, the only requirement for our method is the presence of more than 30 homologous sequences for an accurate extraction of sequence conservation patterns. The prediction method first computes a number of representative well-packed conformations for each pair of contacting TM helices, and then a library of tertiary folds is generated by overlaying overlapping TM helices of the representative conformations. This library is scored using sequence conservation patterns, and a subsequent clustering analysis yields five final models. Assuming that neighboring TM helices in the sequence contact each other (but not that TM helices A and G contact each other), the method produced structural models of Calpha atom root-mean-square deviation (CA RMSD) of 3-5 A from corresponding crystal structures for bacteriorhodopsin, halorhodopsin, sensory rhodopsin II, and rhodopsin. In blind predictions, this type of contact knowledge is not available. Mimicking this, predictions were made for the rotor of the V-type Na(+)-adenosine triphosphatase without such knowledge. The CA RMSD between the best model and its crystal structure is only 3.4 A, and its contact accuracy reaches 55%. Furthermore, the model correctly identifies the binding pocket for sodium ion. These results demonstrate that the method can be readily applied to ab initio structure prediction of simple TM helix bundle proteins having modest numbers of TM helices.     

FoldGPCR: Structure prediction protocol for the transmembrane domain of G protein‐coupled receptors from class A

Building reliable structural models of G protein‐coupled receptors (GPCRs) is a difficult task because of the paucity of suitable templates, low sequence identity, and the wide variety of ligand specificities within the superfamily. Template‐based modeling is known to be the most successful method for protein structure prediction. However, refinement of homology models within 1–3 Å Cα RMSD of the native structure remains a major challenge. Here, we address this problem by developing a novel protocol (foldGPCR) for modeling the transmembrane (TM) region of GPCRs in complex with a ligand, aimed to accurately model the structural divergence between the template and target in the TM helices. The protocol is based on predicted conserved inter‐residue contacts between the template and target, and exploits an all‐atom implicit membrane force field. The placement of the ligand in the binding pocket is guided by biochemical data. The foldGPCR protocol is implemented by a stepwise hierarchical approach, in which the TM helical bundle and the ligand are assembled by simulated annealing trials in the first step, and the receptor‐ligand complex is refined with replica exchange sampling in the second step. The protocol is applied to model the human β₂‐adrenergic receptor (β₂AR) bound to carazolol, using contacts derived from the template structure of bovine rhodopsin. Comparison with the X‐ray crystal structure of the β₂AR shows that our protocol is particularly successful in accurately capturing helix backbone irregularities and helix‐helix packing interactions that distinguish rhodopsin from β₂AR. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Production and initial structural characterization of the TM4TM5 helix‐loop‐helix domain of the translocator protein

Mainly present in the mitochondria, the translocator protein, TSPO, previously known as the peripheral benzodiazepine receptor, is a small essential membrane protein, involved in the translocation of cholesterol across mitochondrial membranes, a rate determining step in steroids biosynthesis. We previously reported the structure of five fragments encompassing the five putative transmembrane helices and showed that each of these fragments constitutes an autonomous folding unit. To further characterize the structural determinants responsible for helix–helix association of this membrane protein, we now investigate the folding of double transmembrane domains in various detergent micelles. Herein, we present the successful biosynthesis of a double transmembrane domain encompassing the last two C‐terminal helices (TM4TM5). For optimal production of this domain in Escherichia coli, the evaluation of various peptide constructs, including TM4TM5 fused to different purification tags or to solubilizing proteins, was necessary. The protocol of production of TM4TM5 with more than 95% purity is reported. This domain was further characterized using circular dichroism and solution state NMR. Far‐UV circular dichroism studies indicate that the secondary structure of TM4TM5 is highly helical when solubilized in various detergent micelles including n‐dodecyl‐β‐d ‐maltoside, n‐octyl‐β‐d ‐glucoside, n‐dodecylphosphocholine, 1,2‐dihexanoyl‐sn‐glycero‐3‐phosphocholine (DHPC), and 1‐palmitoyl‐2‐hydroxy‐sn‐glycero‐3‐phospho‐(1′‐rac‐glycerol). In addition, the solubilization conditions of the domain were optimized for NMR experiments, and preliminary analysis indicates that TM4TM5 adopts a stable tertiary fold within the TM4TM5‐DHPC complex. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Structural Model of the Anion Exchanger 1 (SLC4A1) and Identification of Transmembrane Segments Forming the Transport Site

The anion exchanger 1 (AE1), a member of bicarbonate transporter family SLC4, mediates an electroneutral chloride/bicarbonate exchange in physiological conditions. However, some point mutations in AE1 membrane-spanning domain convert the electroneutral anion exchanger into a Na⁺ and K⁺ conductance or induce a cation leak in a still functional anion exchanger. The molecular determinants that govern ion movement through this transporter are still unknown. The present study was intended to identify the ion translocation pathway within AE1. In the absence of a resolutive three-dimensional structure of AE1 membrane-spanning domain, in silico modeling combined with site-directed mutagenesis experiments was done. A structural model of AE1 membrane-spanning domain is proposed, and this model is based on the structure of a uracil-proton symporter. This model was used to design cysteine-scanning mutagenesis on transmembrane (TM) segments 3 and 5. By measuring AE1 anion exchange activity or cation leak, it is proposed that there is a unique transport site comprising TM3–5 and TM8 that should function as an anion exchanger and a cation leak.     

Design and biophysical characterization of a monomeric four-alpha-helix bundle protein Aα4 with affinity for the volatile anesthetic halothane

A monomeric four-α-helix bundle protein Aα₄ was designed as a step towards investigating the interaction of volatile general anesthetics with their putative membrane protein targets. The alpha helices, connected by glycine loops, have the sequence A, B, B′, A′. The DNA sequence was designed to make the helices with the same amino acid sequences (helix A and A′, B and B′, respectively) as different as possible, while using codons which are favorable for expression in E. coli. The protein was bacterially expressed and purified to homogeneity using reversed-phase HPLC. Protein identity was verified using MALDI–TOF mass spectrometry. Far-UV circular dichroism spectroscopy confirmed the predominantly alpha-helical nature of the protein Aα₄. Guanidinium chloride induced denaturation showed that the monomeric four-α-helix bundle protein Aα₄ is considerably more stable compared to the dimeric di-α-helical protein (Aα₂-L38M)₂. The sigmoidal character of the unfolding reaction is conserved while the sharpness of the transition is increased 1.8-fold. The monomeric four-α-helix bundle protein Aα₄ bound halothane with a dissociation constant (K_d) of 0.93 ± 0.02 mM, as shown by both tryptophan fluorescence quenching and isothermal titration calorimetry. This monomeric four-α-helix bundle protein can now be used as a scaffold to incorporate natural central nervous system membrane protein sequences in order to examine general anesthetic interactions with putative targets in detail.     

Structural characterization of the pore forming protein TatAd of the twin-arginine translocase in membranes by solid-state N-NMR

The transmembrane protein TatA is the pore forming unit of the twin-arginine translocase (Tat), which has the unique ability of transporting folded proteins across the cell membrane. This ATP-independent protein export pathway is a recently discovered alternative to the general secretory (Sec) system of bacteria. To obtain insight in the translocation mechanism, the structure and alignment in the membrane of the well-folded segments 2-45 of TatA_d from Bacillus subtilis was studied here. Using solid-state NMR in bicelles containing anionic lipids, the topology and orientation of TatA_d was determined in an environment mimicking the bacterial membrane. A wheel-like pattern, characteristic for a tilted transmembrane helix, was observed in ¹⁵N chemical shift /¹⁵N-¹H dipolar coupling correlation NMR spectra. Analysis of this PISA wheel revealed a 14-16 residue long N-terminal membrane-spanning helix which is tilted by 17° with respect to the membrane normal. In addition, comparison of uniformly and selectively ¹⁵N-labeled TatA_2-45 samples allowed determination of the helix polarity angle.     

Structural basis for the Trembler-J phenotype of Charcot-Marie-Tooth disease

Mutations in peripheral myelin protein 22 (PMP22) can result in the common peripheral neuropathy Charcot-Marie-Tooth disease (CMTD). The Leu16Pro mutation in PMP22 results in misassembly of the protein, which causes the Trembler-J (TrJ) disease phenotype. Here we elucidate the structural defects present in a partially folded state of TrJ PMP22 that are decisive in promoting CMTD-causing misfolding. In this state, transmembrane helices 2-4 (TM2-4) form a molten globular bundle, while transmembrane helix 1 (TM1) is dissociated from this bundle. The TrJ mutation was seen to profoundly disrupt the TM1 helix, resulting in increased backbone dynamics and changes in the tertiary interactions of TM1 with the PMP22 TM2-4 core in the folded state. Consequently, TM1 undergoes enhanced dissociation from the other transmembrane segments in TrJ PMP22, becoming available for recognition and sequestration by protein-folding quality control, which leads to loss of function and toxic accumulation of aggregates that result in CMTD.