Folding of Aquaporin 1: Multiple evidence that helix 3 can shift out of the membrane core

The folding of most integral membrane proteins follows a two‐step process: initially, individual transmembrane helices are inserted into the membrane by the Sec translocon. Thereafter, these helices fold to shape the final conformation of the protein. However, for some proteins, including Aquaporin 1 (AQP1), the folding appears to follow a more complicated path. AQP1 has been reported to first insert as a four‐helical intermediate, where helix 2 and 4 are not inserted into the membrane. In a second step, this intermediate is folded into a six‐helical topology. During this process, the orientation of the third helix is inverted. Here, we propose a mechanism for how this reorientation could be initiated: first, helix 3 slides out from the membrane core resulting in that the preceding loop enters the membrane. The final conformation could then be formed as helix 2, 3, and 4 are inserted into the membrane and the reentrant regions come together. We find support for the first step in this process by showing that the loop preceding helix 3 can insert into the membrane. Further, hydrophobicity curves, experimentally measured insertion efficiencies and MD‐simulations suggest that the barrier between these two hydrophobic regions is relatively low, supporting the idea that helix 3 can slide out of the membrane core, initiating the rearrangement process.     

Topology mapping to characterize cyanobacterial bicarbonate transporters: BicA (SulP/SLC26 family) and SbtA

This mini-review addresses advances in understanding the transmembrane topologies of two unrelated, single-subunit bicarbonate transporters from cyanobacteria, namely BicA and SbtA. BicA is a Na⁺-dependent bicarbonate transporter that belongs to the SulP/SLC26 family that is widespread in both eukaryotes and prokaryotes. Topology mapping of BicA via the phoA/lacZ fusion reporter method identified 12 transmembrane helices with an unresolved hydrophobic region just beyond helix 8. Re-interpreting this data in the light of a recent topology study on rat prestin leads to a consensus topology of 14 transmembrane domains with a 7+7 inverted repeat structure. SbtA is also a Na⁺-dependent bicarbonate transporter, but of considerably higher affinity (K_m 2–5?μM versus >100?μM for BicA). Whilst SbtA is widespread in cyanobacteria and a few bacteria, it appears to be absent from eukaryotes. Topology mapping of SbtA via the phoA/lacZ fusion reporter method identified 10 transmembrane helices. The topology consists of a 5+5 inverted repeat, with the two repeats separated by a large intracellular loop. The unusual location of the N and C-termini outside the cell raises the possibility that SbtA forms a novel fold, not so far identified by structural and topological studies on transport proteins.     

EM-fold: de novo atomic-detail protein structure determination from medium-resolution density maps

Electron density maps of membrane proteins or large macromolecular complexes are frequently only determined at medium resolution between 4?? and 10??, either by cryo-electron microscopy or X-ray crystallography. In these density maps, the general arrangement of secondary structure elements (SSEs) is revealed, whereas their directionality and connectivity remain elusive. We demonstrate that the topology of proteins with up to 250 amino acids can be determined from such density maps when combined with a computational protein folding protocol. Furthermore, we accurately reconstruct atomic detail in loop regions and amino acid side chains not visible in the experimental data. The EM-Fold algorithm assembles the SSEs de novo before atomic detail is added using Rosetta. In a benchmark of 27 proteins, the protocol consistently and reproducibly achieves models with root mean square deviation values <3??.     

Topology and immersion depth of an integral membrane protein by paramagnetic rates from dissolved oxygen

In studies of membrane proteins, knowledge of protein topology can provide useful insight into both structure and function. In this work, we present a solution NMR method for the measurement the tilt angle and average immersion depth of alpha helices in membrane proteins, from analysis of the paramagnetic relaxation rate enhancements arising from dissolved oxygen. No modification to the micelle or protein is necessary, and the topology of both transmembrane and amphipathic helices are readily determined. We apply this method to the measure the topology of a monomeric mutant of phospholamban (AFA-PLN), a 52-residue membrane protein containing both an amphipathic and a transmembrane alpha helix. In dodecylphosphocholine micelles, the amphipathic helix of AFA-PLN was found to have a tilt angle of 87° ± 1° and an average immersion depth of 13.2 ?. The transmembrane helix was found to have an average immersion depth of 5.4 ?, indicating residues 41 and 42 are closest to the micelle centre. The resolution of paramagnetic relaxation rate enhancements from dissolved oxygen compares favourably to those from Ni (II), a hydrophilic paramagnetic species.     

Bilayer mechanical properties regulate the transmembrane helix mobility and enzymatic state of CD39

CD39 can exist in at least two distinct functional states depending on the presence and intact membrane integration of its two transmembrane helices. In native membranes, the transmembrane helices undergo dynamic rotational motions that are required for enzymatic activity and are regulated by substrate binding. In this study, we show that bilayer mechanical properties regulate conversion between the two enzymatic functional states by modulating transmembrane helix dynamics. Alteration of membrane properties by insertion of cone-shaped or inverse cone-shaped amphiphiles or by cholesterol removal switches CD39 to the same enzymatic state that removal or solubilization of the transmembrane domains does. The same membrane alterations increase the propensity of both transmembrane helices to rotate within the packed structure, resulting in a structure with greater mobility but not an altered primary conformation. Membrane alteration also abolishes the ability of the substrate to stabilize the helices in their primary conformation, indicating a loss of coupling between substrate binding and transmembrane helix dynamics. Removal of either transmembrane helix mimics the effect of membrane alteration on the mobility and substrate sensitivity of the remaining helix, suggesting that the ends of the extracellular domain have intrinsic flexibility. We suggest that a mechanical bilayer property, potentially elasticity, regulates CD39 by altering the balance between the stability and flexibility of its transmembrane helices and, in turn, of its active site.     

Topology of the retinal cone NCKX2 Na/Ca-K exchanger

The Na/Ca-K exchanger (NCKX) is a polytopic membrane protein that plays a critical role in Ca(2+) homeostasis in retinal rod and cone photoreceptors. The NCKX1 isoform is found in rods, while the NCKX2 isoform is found in cones, in retinal ganglion cells, and in various parts of the brain. The topology of the Na/Ca-K exchanger is thought to consist of two large hydrophilic loops and two sets of transmembrane spanning segments (TMs). The first large hydrophilic loop is located extracellularly at the N-terminus; the other is cytoplasmic and separates the two sets of TMs. The TMs consist of either five and five membrane spanning helices or five and six membrane spanning helices, depending upon the predictive algorithm used. Little specific information is yet available on the orientation of the various membrane spanning helices and the localization of the short loops connecting these helices. In this study, we have determined which of the connecting loops are exposed to the extracellular milieu using two different methods: accessibility of substituted cysteine residues and insertion of N-glycosylation sites. The two methods resulted in a consistent NCKX topology in which the two sets of TMs each contain five membrane spanning helices. Our new model places what was previously membrane spanning helix six in the cytoplasm, which places the C-terminus on the extracellular surface. Surprisingly, this NCKX topology model is different from the current NCX topology model with respect to the C-terminal three membrane helices.     

The RTP site shared by the HIV-1 Tat protein and the 11S regulator subunit alpha is crucial for their effects on proteasome function including antigen processing

Nuclease A (NucA) from Anabaena sp. is a non-specific endonuclease able to degrade single and double-stranded DNA and RNA. The endonucleolytic activity is inhibited by the nuclease A inhibitor (NuiA), which binds to NucA with 1:1 stoichiometry and picomolar affinity. In order to better understand the mechanism of inhibition, the solution structure of NuiA was determined by NMR methods. The fold of NuiA is an alpha-beta-alpha sandwich but standard database searches by DALI and TOP revealed no structural homologies. A visual inspection of alpha-beta-alpha folds in the CATH database revealed similarities to the PR-1-like fold (SCOP nomenclature). The similarities include the ordering of secondary structural elements, a single helix on one face of the alpha-beta-alpha sandwich, and three helices on the other face. However, a major difference is in the IV helix, which in the PR-1 fold is short and perpendicular to the I and III helices, but in NuiA is long and parallel to the I and III helices. Additionally, a strand insertion in the beta-sheet makes the NuiA beta-sheet completely antiparallel in organization. The fast time-scale motions of NuiA, characterized by enhanced flexibility of the extended loop between helices III and IV, also show similarities to P14a, which is a PR-1 fold. We propose that the purpose of the PR-1 fold is to form a stable scaffold to present this extended structure for biological interactions with other proteins. This hypothesis is supported by data that show that when NuiA is bound to NucA significant changes in chemical shift occur in the extended loop between helices III and IV.     

Spectroscopy-Based Modelling of the 3D Structure of the β Subunit of the High Affinity IgE Receptor

Abstract

The high affinity IgE receptor, possesses a tetrameric structure. The 243 residue β subunit is a polytopic protein with four hydrophobic membrane-spanning segments, whereas the individual α and γ subunits are bitopic proteins each containing one transmembrane domain in their monomeric form. In the proposed topographical model (Blank et al., 1989), the four trans-membrane α helices of the β subunit are connected by three loop sequences.

To study the individual subunits and intact receptor, this membrane protein was divided into domains such as its loop peptides, cytoplasmic peptides and transmembrane helices according to Blank et al., 1989. The 3D structure of the synthesized loop peptides and cytoplasmic peptides were calculated; CD and/or NMR data were used as appropriate to generate the resultant structures which were then used as data basis for the higher level calculations.

The four individual transmembrane helices of the β subunit were characterised, first of all, by mapping the relative lipophilicity of their surfaces using lipophilic probes. A second procedure, docking of the individual helices in pairs, was used to predict helix–helix interactions.

The data on the relative lipophilicity of the surfaces as well as the surfaces that favoured helix–helix interactions were used in combination with the spectroscopy-based structures of the loops and cytoplasmic domains to calculate via molecular dynamics, the helix arrangement and 3D structure of the β subunit of the high affinity IgE receptor. In the final analysis, the molecular simulations yielded two structures of the β subunit, which should form a basis for the modelling of the whole high affinity IgE receptor.     

Fluorescence resonance energy transfer shows a close helix-helix distance in the transmembrane M13 procoat protein

During the membrane insertion process the major coat protein of bacteriophage M13 assumes a conformation in which two transmembrane helices corresponding to the leader sequence and the anchor region in the mature part of the protein coming into close contact with each other. Previous studies on the molecular mechanism of membrane insertion of M13 procoat protein have shown that this interaction between the two helices might drive the actual translocation process. We investigated the intramolecular distance between the two helices of the transmembrane procoat protein by measuring fluorescence resonance energy transfer (FRET) between the donor (Tyr) placed in one helix and the acceptor (Trp) placed in the other helix. Various mutant procoat proteins with differently positioned donor-acceptor pairs were generated, purified, and reconstituted into artificial lipid bilayers. The results obtained from the FRET measurements, combined with molecular modeling, show that the transmembrane helices are in close contact on the order of 1-1.5 nm. The present approach might be of general interest for determining the topology and the folding of membrane proteins.     

Deuterium/hydrogen exchange factors measured by solution nuclear magnetic resonance spectroscopy as indicators of the structure and topology of membrane proteins

Deuterium/hydrogen exchange factors (chi) were measured for the backbone amide sites of the membrane-bound forms of the 50-residue fd coat protein and the 23-residue magainin2 peptide in lipid micelles by solution nuclear magnetic resonance spectroscopy. By combining kinetic and thermodynamic effects, deuterium/hydrogen exchange factors overcome the principal limitations encountered in the measurements of kinetic protection factors and thermodynamic fractionation factors for membrane proteins. The magnitudes of the exchange factors can be correlated with the structure and topology of membrane-associated polypeptides. In fd coat protein, residues in the transmembrane helix have exchange factors that are substantially smaller than those in the amphipathic surface helix or the loop connecting the two helices. For the amphipathic helical peptide, magainin2, the exchange factors of residues exposed to the solvent are appreciably larger than those that face the hydrocarbon portion of membrane bilayers. These examples demonstrate that deuterium/hydrogen exchange factors can be measured by solution NMR spectroscopy and used to identify residues in transmembrane helices as well as to determine the polarity of amphipathic helices in membrane proteins.     

Current status of membrane protein structure classification

For over 2 decades, continuous efforts to organize the jungle of available protein structures have been underway. Although a number of discrepancies between different classification approaches for soluble proteins have been reported, the classification of membrane proteins has so far not been comparatively studied because of the limited amount of available structural data. Here, we present an analysis of α‐helical membrane protein classification in the SCOP and CATH databases. In the current set of 63 α‐helical membrane protein chains having between 1 and 13 transmembrane helices, we observed a number of differently classified proteins both regarding their domain and fold assignment. The majority of all discrepancies affect single transmembrane helix, two helix hairpin, and four helix bundle domains, while domains with more than five helices are mostly classified consistently between SCOP and CATH. It thus appears that the structural constraints imposed by the lipid bilayer complicate the classification of membrane proteins with only few membrane‐spanning regions. This problem seems to be specific for membrane proteins as soluble four helix bundles, not restrained by the membrane, are more consistently classified by SCOP and CATH. Our findings indicate that the structural space of small membrane helix bundles is highly continuous such that even minor differences in individual classification procedures may lead to a significantly different classification. Membrane proteins with few helices and limited structural diversity only seem to be reasonably classifiable if the definition of a fold is adapted to include more fine‐grained structural features such as helix–helix interactions and reentrant regions. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Mapping the structure and function of the E1 and E2 glycoproteins in alphaviruses

The 9 A resolution cryo-electron microscopy map of Sindbis virus presented here provides structural information on the polypeptide topology of the E2 protein, on the interactions between the E1 and E2 glycoproteins in the formation of a heterodimer, on the difference in conformation of the two types of trimeric spikes, on the interaction between the transmembrane helices of the E1 and E2 proteins, and on the conformational changes that occur when fusing with a host cell. The positions of various markers on the E2 protein established the approximate topology of the E2 structure. The largest conformational differences between the icosahedral surface spikes at icosahedral 3-fold and quasi-3-fold positions are associated with the monomers closest to the 5-fold axes. The long E2 monomers, containing the cell receptor recognition motif at their extremities, are shown to rotate by about 180 degrees and to move away from the center of the spikes during fusion.     

A Continuum theory for the prediction of lateral and rotational positioning of α-helices in membrane proteins: Bacteriorhodopsin

We have developed a new method for the prediction of the lateral and the rotational positioning of transmembrane helices, based upon the present status of knowledge about the dominant interaction of the tertiary structure formation. The basic assumption about the interaction is that the interhelix binding is due to the polar interactions and that very short extramembrane loop segments restrict the relative position of the helices. Another assumption is made for the simplification of the prediction that a helix may be regarded as a continuum rod having polar interaction fields around it. The polar interaction field is calculated by a probe helix method, using a copolymer of serine and alanine as probe helices. The lateral position of helices is determined by the strength of the interhelix binding estimated from the polar interaction field together with the length of linking loop segments. The rotational positioning is determined by the polar interaction field, assuming the optimum lateral configuration. The structural change due to the binding of a prosthetic group is calculated, fixing the rotational freedom of a helix that is connected to the prosthetic group. Applying this method to bacteriorhodopsin, the optimum lateral and rotational positioning of transmembrane helices that are very similar to the experimental configuration was obtained. This method was implemented by a software system, which was developed for this work, and automatic calculation became possible for membrane proteins comprised of several transmembrane helices. © 1995 Wiley-Liss, Inc.     

Ribosome-associated factor Y adopts a fold resembling a double-stranded RNA binding domain scaffold.

Escherichia coli protein Y (pY) binds to the small ribosomal subunit and stabilizes ribosomes against dissociation when bacteria experience environmental stress. pY inhibits translation in vitro, most probably by interfering with the binding of the aminoacyl-tRNA to the ribosomal A site. Such a translational arrest may mediate overall adaptation of cells to environmental conditions. We have determined the 3D solution structure of a 112-residue pY and have studied its backbone dynamic by NMR spectroscopy. The structure has a betaalphabetabetabetaalpha topology and represents a compact two-layered sandwich of two nearly parallel alpha helices packed against the same side of a four-stranded beta sheet. The 23 C-terminal residues of the protein are disordered. Long-range angular constraints provided by residual dipolar coupling data proved critical for precisely defining the position of helix 1. Our data establish that the C-terminal region of helix 1 and the loop linking this helix with strand beta2 show significant conformational exchange in the ms- micro s time scale, which may have relevance to the interaction of pY with ribosomal subunits. Distribution of the conserved residues on the protein surface highlights a positively charged region towards the C-terminal segments of both alpha helices, which most probably constitutes an RNA binding site. The observed betaalphabetabetabetaalpha topology of pY resembles the alphabetabetabetaalpha topology of double-stranded RNA-binding domains, despite limited sequence similarity. It appears probable that functional properties of pY are not identical to those of dsRBDs, as the postulated RNA-binding site in pY does not coincide with the RNA-binding surface of the dsRBDs.     

Vacuolar type H(+)-ATPase genes: presence of four genes including pseudogenes for the 16-kDa proteolipid subunit in the human genome.

Genes for the human vacuolar type H(+)-ATPase proteolipid (16-kDa) subunit were cloned and their nucleotide sequences were determined. Comparison of the deduced sequences indicated that at least four genes including pseudogenes are present in the human genome. One of them corresponded to that for the 16-kDa subunit expressed in HeLa cells. The coding sequence was separated by two introns. The second intron was located in the DNA segment giving a loop between the second and third transmembrane helices, supporting the idea that the 16-kDa subunit was evolved by gene duplication. The primary sequence determined from the second clone had a termination codon behind the third transmembrane helix. Possible translation products from the other two clones had no putative acidic residues essential for proton transport function of the 16-kDa subunit. Thus, it is interesting to know whether these genes are transcribed, since they may have unique cellular functions.     

Transmembrane topology of the Acr3 family arsenite transporter from Bacillus subtilis

The transmembrane topology of the Acr3 family arsenite transporter Acr3 from Bacillus subtilis was analysed experimentally using translational fusions with alkaline phosphatase and green fluorescent protein and in silico by topology modelling. Initial topology prediction resulted in two models with 9 and 10 TM helices respectively. 32 fusion constructs were made between truncated forms of acr3 and the reporter genes at 17 different sites throughout the acr3 sequence to discriminate between these models. Nine strong reporter protein signals provided information about the majority of the locations of the cytoplasmic and extracellular loops of Acr3 and showed that both the N- and the C-termini are located in the cytoplasm. Two ambiguous data points indicated the possibility of an alternative 8 helix topology. This possibility was investigated using another 10 fusion variants, but no experimental support for the 8 TM topology was obtained. We therefore conclude that Acr3 has 10 transmembrane helices. Overall, the loops which connect the membrane spanning segments are short, with cytoplasmic loops being somewhat longer than the extracellular loops. The study provides the first ever experimentally derived structural information on a protein of the Acr3 family which constitutes one of the largest classes of arsenite transporters.     

Simulation studies on bacteriorhodopsin bundle of transmembrane α segments

Bacteriorhodopsin (BR) is a membrane protein which pumps protons through the plasma membrane. Seven transmembrane BR helical segments are subjected to simulation studies in order to investigate the packing process of transmembrane helices. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Helix packing is optimized according to a semi-empirical potential mainly composed of six components: a bilayer potential, a crossing angle potential, a helix dipole potential, a helix-helix distance potential, a helix orientation potential and a helix-helix distance restraint potential (a loop potential). Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. The structures from the simulations are compared with the experimentally determined structures in terms of geometry. The structures generated show similar shapes to the experimentally suggested structure even without the helix-helix distance restraint potential. However, the relative locations of individual helices were reproduced only when the helix-helix distance restraint potential was used with restraint conditions. Our results suggest that transmembrane helix bundles resembling those observed experimentally may be generated by simulations using simple potentials. Received: 19 April 1999 / Revised version: 6 September 1999 / Accepted: 17 September 1999     

Solution structure of the sixth transmembrane helix of the G-protein-coupled receptor, rhodopsin

Low resolution electron density maps have revealed the general orientation of the transmembrane helices of rhodopsin. However, high resolution structural information for the transmembrane domain of the G-protein-coupled receptor, rhodopsin, is as yet unavailable. In this study, a high resolution solution structure is reported for a 15 residue portion of the sixth transmembrane helix of rhodopsin (rhovih) as a free peptide. Helix 6 is one of the transmembrane helices of rhodopsin that contains a proline (amino acid residue 267) and the influence of this proline on the structure of this transmembrane domain was unknown. The structure obtained shows an alpha-helix through most of the sequence. The proline apparently induces only a modest distortion in the helix. Previously, the structure of the intradiskal loop connected to helix 6 was solved. The sequence of this loop contained five residues in common (residues 268-272) with the peptide reported here from the rhovih. The five residues in common between these two structures were superimposed to connect these two structures. The superposition showed a root mean square deviation of 0.2 A. Thus, this five residue sequence formed the same structure in both peptides, indicating that the structure of this region is governed primarily by short range interactions.     

Topography of helices 5-7 in membrane-inserted diphtheria toxin T domain: identification and insertion boundaries of two hydrophobic sequences that do not form a stable transmembrane hairpin

The T domain of diphtheria toxin undergoes a low pH-induced conformational change that allows it to penetrate cell membranes. T domain hydrophobic helices 8 and 9 can adopt two conformations, one close to the membrane surface (P state) and a second in which they apparently form a transmembrane hairpin (TM state). We have now studied T domain helices 5-7, a second cluster of hydrophobic helices, using Cys-scanning mutagenesis. After fluorescently labeling a series of Cys residues, penetration into a non-polar environment, accessibility to externally added antibodies, and relative depth in the bilayer were monitored. It was found that helices 5-7 insert shallowly in the P state and deeply in the TM state. Thus, the conformational changes in helices 5-7 are both similar and somehow linked to those in helices 8 and 9. The boundaries of deeply inserting sequences were also identified. One deeply inserted segment was found to span residues 270 to 290, which overlaps helix 5, and a second spanned residues 300 to 320, which includes most of helix 6 and all of helix 7. This indicates that helices 6 and 7 form a continuous hydrophobic segment despite their separation by a Pro-containing kink. Additionally, it is found that in the TM state some residues in the hydrophilic loop between helices 5 and 6 become more highly exposed than they are in the P state. Their exposure to external solution in the TM state indicates that helices 5-7 do not form a stable transmembrane hairpin. However, helix 5 and/or helices 6 plus 7 could form transmembrane structures that are in equilibrium with non-transmembrane states, or be kinetically prevented from forming a transmembrane structure. How helices 5-7 might influence the mechanism by which the T domain aids translocation of the diphtheria toxin A chain across membranes is discussed.     

Novel inter- and intrasubunit contacts between transport-relevant residues of the homodimeric mitochondrial phosphate transport protein

Ser158 is located near the middle of the matrix loop connecting transmembrane helices C and D of the mitochondrial phosphate transport protein (PTP). The mutant Ser158Thr PTP is transport-inactive. His32 is located near the middle of transmembrane helix A and Thr79 is located 5 residues away from transmembrane helix B and its N-terminal (matrix end). Single site mutant PTPs that have either residue replaced with Ala are transport-inactive. Based on the high resolution structure of a subunit of the bovine ADP/ATP translocase, on sequence similarities between members of the mitochondrial transport protein family, and on the PTP subunit/subunit contact site between transmembrane A helices, it is now suggested that the Ser158 site is at the PTP subunit/subunit contact site. This contact site is essential for keeping the transport cycles catalyzed by the two PTP subunits 180 degrees out of phase. The data also suggest that His32 and Thr79 of the same subunit interact and couple the phosphate and the proton transport paths.