CD39, NTPDase 1, is attached to the plasma membrane by two transmembrane domains. Why?

Since the identification of CD39 and other members of the e-NTPDase (ecto-nucleoside triphosphate diphosphohydrolase) family as the primary enzymes responsible for cell surface nucleotide hydrolysis, one of their most intriguing features has been their unusual topology. The active site lies in the large extracellular region, but instead of being anchored in the membrane by a single transmembrane domain or lipid link like other ectoenzymes, CD39 has two transmembrane domains, one at each end. In this review we discuss evidence that the structure and dynamics of the transmembrane helices are intricately connected to enzymatic function. Removal of either or both transmembrane domains or disruption of their native state by detergent solubilization reduces activity by 90%, indicating that native function requires both transmembrane domains to be present and in the membrane. Enzymatic and mutational analysis of the native and truncated forms has shown that the active site can exist in distinct functional states characterized by different total activities, substrate specificities, hydrolysis mechanisms, and intermediate ADP release during ATP hydrolysis, depending on the state of the transmembrane domains. Disulfide crosslinking of cysteines introduced within the transmembrane helices revealed that they interact within and between molecules, in particular near the extracellular domain, and that activity depends on their organization. Both helices exhibit a high degree of rotational mobility, and the ability to undergo dynamic motions is required for activity and regulated by substrate binding. Recent reports suggest that membrane composition can regulate NTPDase activity. We propose that mechanical bilayer properties, potentially elasticity, might regulate CD39 by altering the balance between stability and mobility of its transmembrane domains.     

Sequence-dependent backbone dynamics of a viral fusogen transmembrane helix

The transmembrane domains of membrane fusogenic proteins are known to contribute to lipid bilayer mixing as indicated by mutational studies and functional reconstitution of peptide mimics. Here, we demonstrate that mutations of a GxxxG motif or of Ile residues, that were previously shown to compromise the fusogenicity of the Vesicular Stomatitis virus G-protein transmembrane helix, reduce its backbone dynamics as determined by deuterium/hydrogen-exchange kinetics. Thus, the backbone dynamics of these helices may be linked to their fusogenicity which is consistent with the known over-representation of Gly and Ile in viral fusogen transmembrane helices. The transmembrane domains of membrane fusogenic proteins are known to contribute to lipid bilayer mixing. Our present results demonstrate that mutations of certain residues, that were previously shown to compromise the fusogenicity of the Vesicular Stomatitis virus G-protein transmembrane helix, reduce its backbone dynamics. Thus, the data suggest a relationship between sequence, backbone dynamics, and fusogenicity of transmembrane segments of viral fusogenic proteins.     

Thermostability of membrane protein helix-helix interaction elucidated by statistical analysis

A prerequisite for the survival of (micro)organisms at high temperatures is an adaptation of protein stability to extreme environmental conditions. In contrast to soluble proteins, where many factors have already been identified, the mechanisms by which the thermostability of membrane proteins is enhanced are almost unknown. The hydrophobic membrane environment constrains possible stabilizing factors for transmembrane domains, so that a difference might be expected between soluble and membrane proteins. Here we present sequence analysis of predicted transmembrane helices of the genomes from eight thermophilic and 12 mesophilic organisms. A comparison of the amino acid compositions indicates that more polar residues can be found in the transmembrane helices of thermophilic organisms. Particularly, the amino acids aspartic acid and glutamic acid replace the corresponding amides. Cysteine residues are found to be significantly decreased by about 70% in thermophilic membrane domains suggesting a non-specific function of most cysteine residues in transmembrane domains of mesophilic organisms. By a pair-motif analysis of the two sets of transmembrane helices, we found that the small residues glycine and serine contribute more to transmembrane helix-helix interactions in thermophilic organisms. This may result in a tighter packing of the helices allowing more hydrogen bond formation.     

GrinchGEF--a novel Rho-specific guanine nucleotide exchange factor

RhoGTPases, which are activated by specific guanine nucleotide exchange factors (GEFs), play pivotal roles in several cellular functions. We identified a new RhoGEF (GrinchGEF) containing the typical Dbl homology domain, a putative WD40-like domain, and two predicted transmembrane helices. In contrast to most other RhoGEFs, it exhibits no sequence similarities to known pleckstrin homology domains. GrinchGEF mRNA was highly abundant in skeletal muscle and pancreas. Despite the predicted transmembrane domains, subcellular localization studies revealed a cytosolic distribution. In vitro, GrinchGEF induced the GDP/GTP exchange at RhoA, but not at Rac1 or Cdc42. In intact cells, GrinchGEF induced specifically Rho activation and enhanced RhoA-C-specific downstream effects.     

Sequence context and modified hydrophobic moment plots help identify 'horizontal' surface helices in transmembrane protein structure prediction

Transmembrane proteins make up at least one-fifth of the genome of most organisms and are critical components of key pathways for cell survival and interactions with the environment. The function of helices found at the membrane surface in transmembrane proteins has not been greatly explored, but it is likely that they play an ancillary role to membrane spanning helices and are analogous to the surface active helices of peripheral membrane proteins, being involved in: lipid association, membrane perturbation, transmembrane signal transduction and regulation, and transmembrane helical bundle formation. Due to the difficulties in obtaining high-resolution structural data for this class of proteins, structure-from-sequence predictive methods continue to be developed as a means to obtain structural models for these largely intractable systems. A simple but effective variant of the hydrophobic moment analysis of amino acid sequences is described here as part of a protocol for distinguishing helical sequences that are parallel to or 'horizontal' at the membrane bilayer/aqueous phase interface from helices that are membrane-embedded or located in extra-membranous domains. This protocol when tested on transmembrane spanning protein amino acid sequences not used in its development, was found to be 84-91% accurate when the results were compared to the partition locations in the corresponding structures determined by X-ray crystallography, and 72% accurate in determining which helices lie horizontal or near horizontal at the lipid interface.     

Structural basis of ligand activation in a cyclic nucleotide regulated potassium channel

Here we describe the initial functional characterization of a cyclic nucleotide regulated ion channel from the bacterium Mesorhizobium loti and present two structures of its cyclic nucleotide binding domain, with and without cAMP. The domains are organized as dimers with the interface formed by the linker regions that connect the nucleotide binding pocket to the pore domain. Together, structural and functional data suggest the domains form two dimers on the cytoplasmic face of the channel. We propose a model for gating in which ligand binding alters the structural relationship within a dimer, directly affecting the position of the adjacent transmembrane helices.     

Major structural determinants of transmembrane proteins identified by principal component analysis

We identify amino acid characteristics important in determining the secondary structures of transmembrane proteins, and compare them with characteristics important for cytoplasmic proteins. Using information derived from multiple sequence alignments, we perform a principal component analysis (PCA) to identify the directions in the 20-dimensional amino acid frequency space that comprise the most variance within each protein secondary structure. These vectors represent the important position-specific properties of the amino acids for coils, turns, beta sheets, and alpha helices. As expected, the most important axis for most of the datasets was hydrophobicity. Additional axes, distinct from hydrophobicity, are surprising, especially in the case of transmembrane alpha helices, where the effects of aromaticity and beta-branching are the next two most significant characteristics. The axis representing beta-branching also has equal importance in cytoplasmic and transmembrane helices, a finding that contrasts with some experimental results in membrane-like environments. In a further analysis, we examine trends for some of the PCA axes over averaged transmembrane alpha helices, and find interesting results for aromaticity.     

Phenylalanine promotes interaction of transmembrane domains via GxxxG motifs

Interactions of transmembrane helices play a crucial role in the folding and oligomerisation of integral membrane proteins. In order to uncover novel sequence motifs mediating these interactions, we randomised one face of a transmembrane helix with a set of non-polar or moderately polar amino acids. Those sequences capable of self-interaction upon integration into bacterial inner membranes were selected by means of the ToxR/POSSYCCAT system. A comparison between low/medium-affinity and high-affinity sequences reveals that high-affinity sequences are strongly enriched in phenylalanine residues that are frequently observed at the − 3 position of GxxxG motifs, thus yielding FxxGxxxG motifs. Mutation of Phe or GxxxG in selected sequences significantly reduces self-interaction of the transmembrane domains without affecting their efficiency of membrane integration. Conversely, grafting FxxGxxxG onto unrelated transmembrane domains strongly enhances their interaction. Further, we find that FxxGxxxG is significantly over-represented in transmembrane domains of bitopic membrane proteins. The same motif contributes to self-interaction of the vesicular stomatitis virus G protein transmembrane domain. We conclude that Phe stabilises membrane-spanning GxxxG motifs. This is one example of how the role of certain side-chains in helix-helix interfaces is modulated by sequence context.     

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

SNARE proteins mediate fusion of intracellular eukaryotic membranes and their α-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.     

Membrane assembly of simple helix homo-oligomers studied via molecular dynamics simulations

The assembly of simple transmembrane helix homo-oligomers is studied by combining a generalized Born implicit membrane model with replica exchange molecular dynamics simulations to sample the conformational space of various oligomerization states and the native oligomeric conformation. Our approach is applied to predict the structures of transmembrane helices of three proteins--glycophorin A, the M2 proton channel, and phospholamban--using only peptide sequence and the native oligomerization state information. In every case, the methodology reproduces native conformations that are in good agreement with available experimental structural data. Thus, our method should be useful in the prediction of native structures of transmembrane domains of other peptides. When we ignore the experimental constraint on the native oligomerization state and attempt de novo prediction of the structure and oligomerization state based only on sequence and simple energetic considerations, we identify the pentamer as the most stable oligomer for phospholamban. However, for the glycophorin A and the M2 proton channels, we tend to predict higher oligomers as more stable. Our studies demonstrate that reliable predictions of the structure of transmembrane helical oligomers can be achieved when the observed oligomerization state is imposed as a constraint, but that further efforts are needed for the de novo prediction of both structure and oligomeric state.     

Alpha-helical coiled-coil oligomerization domains are almost ubiquitous in the collagen superfamily

Alpha-helical coiled-coils are widely occurring protein oligomerization motifs. Here we show that most members of the collagen superfamily contain short, repeating heptad sequences typical of coiled coils. Such sequences are found at the N-terminal ends of the C-propeptide domains in all fibrillar procollagens. When fused C-terminal to a reporter molecule containing a collagen-like sequence that does not spontaneously trimerize, the C-propeptide heptad repeats induced trimerization. C-terminal heptad repeats were also found in the oligomerization domains of the multiplexins (collagens XV and XVIII). N-terminal heptad repeats are known to drive trimerization in transmembrane collagens, whereas fibril-associated collagens with interrupted triple helices, as well as collagens VII, XIII, XXIII, and XXV, were found to contain heptad repeats between collagen domains. Finally, heptad repeats were found in the von Willebrand factor A domains known to be involved in trimerization of collagen VI, as well as in collagen VII. These observations suggest that coiled-coil oligomerization domains are widely used in the assembly of collagens and collagen-like proteins.     

Identification of a new gene (rat TM6P1) encoding a fasting-inducible, integral membrane protein with six transmembrane domains

We describe the isolation of a new gene that encodes a membrane-integrated protein with six transmembrane domains, termed TM6P1. A 403-bp expressed sequence tag was isolated from fasted rat liver subtracted cDNA library, and its full-length cDNA is 1482 bp long. It contains an open reading frame of 816 bp and is predicted to encode a 271-amino acid protein with a deduced mass of 29520 Da. A sequence homology search failed to show significant correspondence to any known protein in the databank. TM6P1 has six highly hydrophobic domains that are predicted to be transmembrane helices. Consistent with this prediction, the TM6P1-EGFP fusion protein was shown to localize to the plasma membrane. TM6P1 mRNA is widely expressed in rat tissues, with placenta and liver being the most abundant sites. Fasting increased TM6P1 mRNA nearly two-fold in liver. Taken together, our data suggest that TM6P1 is a unique new membrane integral protein that might have a function important during fasting-induced catabolism.     

Nucleotide bias causes a genomewide bias in the amino acid composition of proteins

We analyzed the nucleotide contents of several completely sequenced genomes, and we show that nucleotide bias can have a dramatic effect on the amino acid composition of the encoded proteins. By surveying the genes in 21 completely sequenced eubacterial and archaeal genomes, along with the entire Saccharomyces cerevisiae genome and two Plasmodium falciparum chromosomes, we show that biased DNA encodes biased proteins on a genomewide scale. The predicted bias affects virtually all genes within the genome, and it could be clearly seen even when we limited the analysis to sets of homologous gene sequences. Parallel patterns of compositional bias were found within the archaea and the eubacteria. We also found a positive correlation between the degree of amino acid bias and the magnitude of protein sequence divergence. We conclude that mutational bias can have a major effect on the molecular evolution of proteins. These results could have important implications for the interpretation of protein-based molecular phylogenies and for the inference of functional protein adaptation from comparative sequence data.     

Nuclear magnetic resonance studies of the hammerhead ribozyme domain. Secondary structure formation and magnesium ion dependence

Proton nuclear magnetic resonance (n.m.r.) experiments were used to probe base-pair formation in several hammerhead RNA enzyme (ribozyme) domains. The hammerhead domains consist of a 34 nucleotide ribozyme bound to a complementary 13 nucleotide non-cleavable DNA substrate. Three hammerhead domains were studied that differ in the sequence and stability of one of the helices involved in recognition of the substrate by the ribozyme. The n.m.r. data show a 1:1 stoichiometry for the ribozyme-substrate complexes. The imino proton resonances in the hammerhead complexes were assigned by two-dimensional nuclear Overhauser effect experiments. These data confirm the presence of two of the three helical regions in the hammerhead domain, predicted from phylogenetic data; and are also consistent with the formation of the third helix. Since a divalent cation is required for efficient catalytic activity of the hammerhead domain, the magnesium ion dependence of the n.m.r. spectra was studied for two of the hammerhead complexes. One of the complexes showed very large spectral changes upon addition of magnesium ions. However, the complex that has the most C.G base-pairs in one of the recognition helices shows essentially no spectral (and therefore presumably structural) changes upon addition of magnesium. These data are consistent with a model where the magnesium binding site already exists in the magnesium-free complex, suggesting that the magnesium ion serves primarily a catalytic, and not a structural, role under the conditions used here.     

Three-dimensional structures of the mammalian multidrug resistance P-glycoprotein demonstrate major conformational changes in the transmembrane domains upon nucleotide binding

P-glycoprotein is an ATP-binding cassette transporter that is associated with multidrug resistance and the failure of chemotherapy in human patients. We have previously shown, based on two-dimensional projection maps, that P-glycoprotein undergoes conformational changes upon binding of nucleotide to the intracellular nucleotide binding domains. Here we present the three-dimensional structures of P-glycoprotein in the presence and absence of nucleotide, at a resolution limit of approximately 2 nm, determined by electron crystallography of negatively stained crystals. The data reveal a major reorganization of the transmembrane domains throughout the entire depth of the membrane upon binding of nucleotide. In the absence of nucleotide, the two transmembrane domains form a single barrel 5-6 nm in diameter and about 5 nm deep with a central pore that is open to the extracellular surface and spans much of the membrane depth. Upon binding nucleotide, the transmembrane domains reorganize into three compact domains that are each 2-3 nm in diameter and 5-6 nm deep. This reorganization opens the central pore along its length in a manner that could allow access of hydrophobic drugs (transport substrates) directly from the lipid bilayer to the central pore of the transporter.     

How important are transmembrane helices of bitopic membrane proteins?

The topology of a bitopic membrane protein consists of a single transmembrane helix connecting two extra-membranous domains. As opposed to helices from polytopic proteins, the transmembrane helices of bitopic proteins were initially considered as merely hydrophobic anchors, while more recent studies have begun to shed light on their role in the protein's function. Herein the overall importance of transmembrane helices from bitopic membrane proteins was analyzed using a relative conservation analysis. Interestingly, the transmembrane domains of bitopic proteins are on average, significantly more conserved than the remainder of the protein, even when taking into account their smaller amino acid repertoire. Analysis of highly conserved transmembrane domains did not reveal any unifying consensus, pointing to a great diversity in their conservation patterns. However, Fourier power spectrum analysis was able to show that regardless of the conservation motif, in most sequences a significant conservation moment was observed, in that one side of the helix was conserved while the other was not. Taken together, it may be possible to conclude that a significant proportion of transmembrane helices from bitopic membrane proteins participate in specific interactions, in a variety of modes in the plane of the lipid bilayer.     

How important are transmembrane helices of bitopic membrane proteins?

The topology of a bitopic membrane protein consists of a single transmembrane helix connecting two extra-membranous domains. As opposed to helices from polytopic proteins, the transmembrane helices of bitopic proteins were initially considered as merely hydrophobic anchors, while more recent studies have begun to shed light on their role in the protein's function. Herein the overall importance of transmembrane helices from bitopic membrane proteins was analyzed using a relative conservation analysis. Interestingly, the transmembrane domains of bitopic proteins are on average, significantly more conserved than the remainder of the protein, even when taking into account their smaller amino acid repertoire. Analysis of highly conserved transmembrane domains did not reveal any unifying consensus, pointing to a great diversity in their conservation patterns. However, Fourier power spectrum analysis was able to show that regardless of the conservation motif, in most sequences a significant conservation moment was observed, in that one side of the helix was conserved while the other was not. Taken together, it may be possible to conclude that a significant proportion of transmembrane helices from bitopic membrane proteins participate in specific interactions, in a variety of modes in the plane of the lipid bilayer.     

Structural consequences of ATP hydrolysis on the ABC transporter NBD dimer: molecular dynamics studies of HlyB

ABC transporters are a large and important family of membrane proteins involved in substrate transport across the membrane. The transported substrates are quite diverse, ranging from monatomic ions to large biomolecules. Consequently, some ABC transporters are involved in biomedically relevant situations, from genetic diseases to multidrug resistance. The most conserved domains in ABC transporters are the nucleotide binding domains (NBDs), which form a dimer responsible for the binding and hydrolysis of ATP, concomitantly with substrate translocation. To elucidate how ATP hydrolysis structurally affects the NBD dimer, and consequently the transporter, we performed a molecular dynamics study on the NBD dimer of the HlyB ABC exporter. We have observed a change in the contact surface between the monomers after hydrolysis, even though we have not seen dimer opening in any of the five 100 ns simulations. We have also identified specific regions that respond to ATP hydrolysis, in particular the X-loop motif of ABC exporters, which has been shown to be in contact with the coupling helices of the transmembrane domains (TMDs). We propose that this motif is an important part of the NBD-TMD communication in ABC exporters. Through nonequilibrium analysis, we have also identified gradual conformational changes within a short time scale after ATP hydrolysis.     

Calculation of the three-dimensional structure of Saccharomyces cerevisiae cytochrome b inserted in a lipid matrix

Cytochrome b is an integral membrane protein, which forms the core of the ubiquinol-cytochrome c oxidoreductase (cytochrome bc1) complex. A computer-aided three-dimensional modeling procedure was carried out in four steps. First, the candidate hydrophobic helices were searched for throughout the protein primary sequence by a computer procedure based upon the method of Eisenberg; second, a secondary helical structure was imposed to the transmembrane peptides; third, the helical segments at a lipid-water interface were oriented, and finally the possible interactions between helices with similar properties were investigated. This procedure enabled the identification of nine hydrophobic segments, of which eight are membrane-spanning helices while one has amphipathic properties. Three hydrophilic receptor-binding domains were also identified. Based upon their hydrophobicity profiles, the transmembrane helices could be associated in pairs inside the lipid bilayer. In our folding model proposed for cytochrome b, all mutation sites are not only located on the same side of the membrane but are also in close proximity in the three-dimensional structure. Inhibitor resistance mutational sites which were recently characterized (di Rago, J.-P., and Colson, A.-M. (1988) J. Biol. Chem. 263, 12564-12570) have been located on this model. Moreover, the receptor-binding domains and the mutation sites are close neighbors in the three-dimensional spatial representation.     

Two ancient bacterial endosymbionts have coevolved with the planthoppers (Insecta: Hemiptera: Fulgoroidea)

Background

Pseudoscorpions are chelicerates and have historically been viewed as being most closely related to solifuges, harvestmen, and scorpions. No mitochondrial genomes of pseudoscorpions have been published, but the mitochondrial genomes of some lineages of Chelicerata possess unusual features, including short rRNA genes and tRNA genes that lack sequence to encode arms of the canonical cloverleaf-shaped tRNA. Additionally, some chelicerates possess an atypical guanine-thymine nucleotide bias on the major coding strand of their mitochondrial genomes.

Results

We sequenced the mitochondrial genomes of two divergent taxa from the chelicerate order Pseudoscorpiones. We find that these genomes possess unusually short tRNA genes that do not encode cloverleaf-shaped tRNA structures. Indeed, in one genome, all 22 tRNA genes lack sequence to encode canonical cloverleaf structures. We also find that the large ribosomal RNA genes are substantially shorter than those of most arthropods. We inferred secondary structures of the LSU rRNAs from both pseudoscorpions, and find that they have lost multiple helices. Based on comparisons with the crystal structure of the bacterial ribosome, two of these helices were likely contact points with tRNA T-arms or D-arms as they pass through the ribosome during protein synthesis. The mitochondrial gene arrangements of both pseudoscorpions differ from the ancestral chelicerate gene arrangement. One genome is rearranged with respect to the location of protein-coding genes, the small rRNA gene, and at least 8 tRNA genes. The other genome contains 6 tRNA genes in novel locations. Most chelicerates with rearranged mitochondrial genes show a genome-wide reversal of the CA nucleotide bias typical for arthropods on their major coding strand, and instead possess a GT bias. Yet despite their extensive rearrangement, these pseudoscorpion mitochondrial genomes possess a CA bias on the major coding strand. Phylogenetic analyses of all 13 mitochondrial protein-coding gene sequences consistently yield trees that place pseudoscorpions as sister to acariform mites.

Conclusion

The well-supported phylogenetic placement of pseudoscorpions as sister to Acariformes differs from some previous analyses based on morphology. However, these two lineages share multiple molecular evolutionary traits, including substantial mitochondrial genome rearrangements, extensive nucleotide substitution, and loss of helices in their inferred tRNA and rRNA structures.