Mode matches and their locations in the hydrophobic free energy sequences of peptide ligands and their receptor eigenfunctions

Patterns in sequences of amino acid hydrophobic free energies predict secondary structures in proteins. In protein folding, matches in hydrophobic free energy statistical wavelengths appear to contribute to selective aggregation of secondary structures in “hydrophobic zippers.” In a similar setting, the use of Fourier analysis to characterize the dominant statistical wavelengths of peptide ligands’ and receptor proteins’ hydrophobic modes to predict such matches has been limited by the aliasing and end effects of short peptide lengths, as well as the broad-band, mode multiplicity of many of their frequency (power) spectra. In addition, the sequence locations of the matching modes are lost in this transformation. We make new use of three techniques to address these difficulties: (i) eigenfunction construction from the linear decomposition of the lagged covariance matrices of the ligands and receptors as hydrophobic free energy sequences; (ii) maximum entropy, complex poles power spectra, which select the dominant modes of the hydrophobic free energy sequences or their eigenfunctions; and (iii) discrete, best bases, trigonometric wavelet transformations, which confirm the dominant spectral frequencies of the eigenfunctions and locate them as (absolute valued) moduli in the peptide or receptor sequence. The leading eigenfunction of the covariance matrix of a transmembrane receptor sequence locates the same transmembrane segments seen in n-block-averaged hydropathy plots while leaving the remaining hydrophobic modes unsmoothed and available for further analyses as secondary eigenfunctions. In these receptor eigenfunctions, we find a set of statistical wavelength matches between peptide ligands and their G-protein and tyrosine kinase coupled receptors, ranging across examples from 13.10 amino acids in acid fibroblast growth factor to 2.18 residues in corticotropin releasing factor. We find that the wavelet-located receptor modes in the extracellular loops are compatible with studies of receptor chimeric exchanges and point mutations. A nonbinding corticotropin-releasing factor receptor mutant is shown to have lost the signatory mode common to the normal receptor and its ligand. Hydrophobic free energy eigenfunctions and their transformations offer new quantitative physical homologies in database searches for peptide-receptor matches.

Multidrug resistance mediated by a bacterial homolog of the human multidrug transporter MDR1.

Resistance of Lactococcus lactis to cytotoxic compounds shares features with the multidrug resistance phenotype of mammalian tumor cells. Here, we report the gene cloning and functional characterization in Escherichia coli of LmrA, a lactococcal structural and functional homolog of the human multidrug resistance P-glycoprotein MDR1. LmrA is a 590-aa polypeptide that has a putative topology of six alpha-helical transmembrane segments in the N-terminal hydrophobic domain, followed by a hydrophilic domain containing the ATP-binding site. LmrA is similar to each of the two halves of MDR1 and may function as a homodimer. The sequence conservation between LmrA and MDR1 includes particular regions in the transmembrane domains and connecting loops, which, in MDR1 and the MDR1 homologs in other mammalian species, have been implicated as determinants of drug recognition and binding. LmrA and MDR1 extrude a similar spectrum of amphiphilic cationic compounds, and the activity of both systems is reversed by reserpine and verapamil. As LmrA can be functionally expressed in E. coli, it offers a useful prokaryotic model for future studies on the molecular mechanism of MDR1-like multidrug transporters.

Interaction of the protein import and folding machineries of the chloroplast.

We report the molecular cloning of import intermediate associated protein (IAP) 100, a 100-kDa protein of the chloroplast protein import machinery of peas. IAP100 contains two potential alpha-helical transmembrane segments and also behaves like an integral membrane protein. It was localized to the inner chloroplast envelope membrane. Immunoprecipitation experiments using monospecific anti-IAP100 antibodies and a nonionic detergent-generated chloroplast lysate gave the following results. (i) The four integral membrane proteins of the outer chloroplast import machinery were not coprecipitated with IAP100 indicating that the inner and outer membrane import machineries are not coupled in isolated chloroplasts. (ii) the major protein that coprecipitated with IAP100 was identified as stromal chaperonin 60 (cpn60); the association of IAP100 and cpn60 was specific and was abolished when immunoprecipitation was carried out in the presence of ATP. (iii) In a lysate from chloroplasts that had been preincubated for various lengths of time in an import reaction with radiolabeled precursor (pS) of the small subunit of Rubisco, we detected coimmunoprecipitation of IAP100, cpn60, and the imported mature form (S) of precursor. Relative to the time course of import, coprecipitation of S first increased and then decreased, consistent with a transient association of the newly imported S with the chaperonin bound to IAP100. These data suggest that IAP100 serves in recruiting chaperonin for folding of newly imported proteins.

Mutational analysis and molecular modeling of the nonapeptide hormone binding domains of the [Arg8]vasotocin receptor.

To identify determinants that form nonapeptide hormone binding domains of the white sucker Catostomus commersoni [Arg8]vasotocin receptor, chimeric constructs encoding parts of the vasotocin receptor and parts of the isotocin receptor have been analyzed by [(3,5-3H)Tyr2, Arg8]vasotocin binding to membranes of human embryonic kidney cells previously transfected with the different cDNA constructs and by functional expression studies in Xenopus laevis oocytes injected with mutant cRNAs. The results indicate that the N terminus and a region spanning the second extracellular loop and its flanking transmembrane segments, which contains a number of amino acid residues that are conserved throughout the nonapeptide receptor family, contribute to the affinity of the receptor for its ligand. Nonapeptide selectivity, however, is mainly defined by transmembrane region VI and the third extracellular loop. These results are complemented by a molecular model of the vasotocin receptor obtained by aligning its sequence with those of other G-protein coupled receptors as well as that of bacteriorhodopsin. The model indicates that amino acid residues of transmembrane regions II-VII that are located close to the extracellular surface also contribute to the binding of vasotocin.

Signal transduction in the archaeon Halobacterium salinarium is processed through three subfamilies of 13 soluble and membrane-bound transducer proteins.

Eubacterial transducers are transmembrane, methyl-accepting proteins central to chemotaxis systems and share common structural features. We identified a large family of transducer proteins in the Archaeon Halobacterium salinarium using a site-specific multiple antigenic peptide antibody raised against 23 amino acids, representing the highest homology region of eubacterial transducers. This immunological observation was confirmed by isolating 13 methyl-accepting taxis genes using a 27-mer oligonucleotide probe, corresponding to conserved regions between the eubacterial and first halobacterial phototaxis transducer gene htrI. On the basis of the comparison of the predicted structural domains of these transducers, we propose that at least three distinct subfamilies of transducers exist in the Archaeon H. salinarium: (i) a eubacterial chemotaxis transducer type with two hydrophobic membrane-spanning segments connecting sizable domains in the periplasm and cytoplasm; (ii) a cytoplasmic domain and two or more hydrophobic transmembrane segments without periplasmic domains; and (iii) a cytoplasmic domain without hydrophobic transmembrane segments. We fractionated the halobacterial cell lysate into soluble and membrane fractions and localized different halobacterial methyl-accepting taxis proteins in both fractions.

Multidrug resistance proteins QacA and QacB from Staphylococcus aureus: membrane topology and identification of residues involved in substrate specificity.

The closely related multidrug efflux pumps QacA and QacB, from the bacterial pathogen Staphylococcus aureus, both confer resistance to various toxic organic cations but differ in that QacB mediates lower levels of resistance to divalent cations. Cloning and nucleotide sequencing of the qacB gene revealed that qacB differs from qacA by only seven nucleotide substitutions. Random hydroxylamine mutagenesis of qacB was undertaken, selecting for variants that conferred increased resistance to divalent cations. Both QacA and the QacB mutants capable of conferring resistance to divalent cations contain an acidic residue at either amino acid 322 or 323, whereas QacB contains uncharged residues in these positions. Site-directed mutagenesis of qacA confirmed the importance of an acidic residue within this region of QacA in conferring resistance to divalent cations. Membrane topological analysis using alkaline phosphatase and beta-galactosidase fusions indicated that the QacA protein contains 14 transmembrane segments. Thus, QacA represents the first membrane transport protein shown to contain 14 transmembrane segments, and confirms that the major facilitator superfamily contains a family of proteins with 14 transmembrane segments.

Nramp defines a family of membrane proteins.

Nramp (natural resistance-associated macrophage protein) is a newly identified family of integral membrane proteins whose biochemical function is unknown. We report on the identification of Nramp homologs from the fly Drosophila melanogaster, the plant Oryza sativa, and the yeast Saccharomyces cerevisiae. Optimal alignment of protein sequences required insertion of very few gaps and revealed remarkable sequence identity of 28% (yeast), 40% (plant), and 55% (fly) with the mammalian proteins (46%, 58%, and 73% similarity), as well as a common predicted transmembrane topology. This family is defined by a highly conserved hydrophobic core encoding 10 transmembrane segments. Other features of this hydrophobic core include several invariant charged residues, helical periodicity of sequence conservation suggesting conserved and nonconserved faces for several transmembrane helices, a consensus transport signature on the intracytoplasmic face of the membrane, and structural determinants previously described in ion channels. These characteristics suggest that the Nramp polypeptides form part of a group of transporters or channels that act on as yet unidentified substrates.

Structural conservation of ion conduction pathways in K channels and glutamate receptors.

Single channel recordings demonstrate that ion channels switch stochastically between an open and a closed pore conformation. In search of a structural explanation for this universal open/close behavior, we have uncovered a striking degree of amino acid homology across the pore-forming regions of voltage-gated K channels and glutamate receptors. This suggested that the pores of these otherwise unrelated classes of channels could be structurally conserved. Strong experimental evidence supports a hairpin structure for the pore-forming region of K channels. Consequently, we hypothesized the existence of a similar structure for the pore of glutamate receptors. In ligand-gated channels, the pore is formed by M2, the second of four putative transmembrane segments. A hairpin structure for M2 would affect the subsequent membrane topology, inverting the proposed orientation of the next segments, M3. We have tested this idea for the NR1 subunit of the N-methyl-D-aspartate receptor. Mutations that affected the glycosylation pattern of the NR1 subunit localize both extremes of the M3-M4 linker to the extracellular space. Whole cell currents and apparent agonist affinities were not affected by these mutations. Therefore it can be assumed that they represent the native transmembrane topology. The extracellular assignment of the M3-M4 linker challenged the current topology model by inverting M3. Taken together, the amino acid homology and the new topology suggest that the pore-forming M2 segment of glutamate receptors does not transverse the membrane but, rather, forms a hairpin structure, similar to that found in K channels.

The LAR/PTP delta/PTP sigma subfamily of transmembrane protein-tyrosine-phosphatases: multiple human LAR, PTP delta, and PTP sigma isoforms are expressed in a tissue-specific manner and associate with the LAR-interacting protein LIP.1.

The transmembrane protein-tyrosine-phosphatases (PTPases) LAR, PTP delta, and PTP sigma each contain two intracellular PTPase domains and an extracellular region consisting of Ig-like and fibronectin type III-like domains. We describe the cloning and characterization of human PTP sigma (HPTP sigma) and compare the structure, alternative splicing, tissue distribution, and PTPase activity of LAR, HPTP delta, and HPTP sigma, as well their ability to associate with the intracellular coiled-coil LAR-interacting protein LIP.1. Overall, these three PTPases are structurally very similar, sharing 64% amino acid identity. Multiple isoforms of LAR, HPTP delta, and HPTP sigma appear to be generated by tissue-specific alternative splicing of up to four mini-exon segments that encode peptides of 4-16 aa located in both the extracellular and intracellular regions. Alternative usage of these peptides varies depending on the tissue mRNA analyzed. Short isoforms of both HPTP sigma and HPTP delta were also detected that contain only four of the eight fibronectin type III-like domains. Northern blot analysis indicates that LAR and HPTP sigma are broadly distributed whereas HPTP delta expression is largely restricted to brain, as is the short HPTP sigma isoform containing only four fibronectin type III-like domains. LAR, HPTP delta, and HPTP sigma exhibit similar in vitro PTPase activities and all three interact with LIP.1, which has been postulated to recruit LAR to focal adhesions. Thus, these closely related PTPases may perform similar functions in various tissues.

Probing the transmembrane topology of cyclic nucleotide-gated ion channels with a gene fusion approach.

Cyclic nucleotide-gated (CNG) cation channels contain two short sequence motifs--a residual voltage-sensor (S4) and a pore-forming (P) segment--that are reminiscent of similar segments in voltage-activated Shaker-type K+ channels. It has been tacitly assumed that CNG channels and this K+ channel subfamily share a common overall topology, characterized by a hydrophobic domain comprising six membrane-spanning segments. We have systematically investigated the topology of CNG channels from bovine rod photoreceptor and Drosophila melanogaster by a gene fusion approach using the bacterial reporter enzymes alkaline phosphatase and beta-galactosidase, which are active only in the periplasm and only in the cytoplasm, respectively. Enzymatic activity was determined after expression of fusion constructs in Escherichia coli. CNG channels were found to have six membrane-spanning segments, suggesting that CNG and Shaker-type K+ channels, albeit distant relatives within a gene superfamily of ion channels, share a common topology.

On the distribution of amino acid residues in transmembrane alpha-helix bundles.

The periodic distribution of residues in the sequence of 469 putative transmembrane alpha-helices from eukaryotic plasma membrane polytopic proteins has been analyzed with correlation matrices. The method does not involve any a priori assumption about the secondary structure of the segments or about the physicochemical properties of individual amino acid residues. Maximal correlation is observed at 3.6 residues per period, characteristic of alpha-helices. A scale extracted from the data describes the propensity of the various residues to lie on the same or on opposite helix faces. The most polar face of transmembrane helices, presumably that buried in the protein core, shows a strong enrichment in aromatic residues, while residues likely to face the fatty acyl chains of lipids are largely aliphatic.

Molecular basis of fast inactivation in voltage and Ca2+-activated K+ channels: A transmembrane β-subunit homolog

Voltage-dependent and calcium-sensitive K+ (MaxiK) channels are key regulators of neuronal excitability, secretion, and vascular tone because of their ability to sense transmembrane voltage and intracellular Ca2+. In most tissues, their stimulation results in a noninactivating hyperpolarizing K+ current that reduces excitability. In addition to noninactivating MaxiK currents, an inactivating MaxiK channel phenotype is found in cells like chromaffin cells and hippocampal neurons. The molecular determinants underlying inactivating MaxiK channels remain unknown. Herein, we report a transmembrane β subunit (β2) that yields inactivating MaxiK currents on coexpression with the pore-forming α subunit of MaxiK channels. Intracellular application of trypsin as well as deletion of 19 N-terminal amino acids of the β2 subunit abolished inactivation of the α subunit. Conversely, fusion of these N-terminal amino acids to the noninactivating smooth muscle β1 subunit leads to an inactivating phenotype of MaxiK channels. Furthermore, addition of a synthetic N-terminal peptide of the β2 subunit causes inactivation of the MaxiK channel α subunit by occluding its K+-conducting pore resembling the inactivation caused by the “ball” peptide in voltage-dependent K+ channels. Thus, the inactivating phenotype of MaxiK channels in native tissues can result from the association with different β subunits.

A mutation in the transmembrane/luminal domain of the ryanodine receptor is associated with abnormal Ca2+ release channel function and severe central core disease

Central core disease is a rare, nonprogressive myopathy that is characterized by hypotonia and proximal muscle weakness. In a large Mexican kindred with an unusually severe and highly penetrant form of the disorder, DNA sequencing identified an I4898T mutation in the C-terminal transmembrane/luminal region of the RyR1 protein that constitutes the skeletal muscle ryanodine receptor. All previously reported RYR1 mutations are located either in the cytoplasmic N terminus or in a central cytoplasmic region of the 5,038-aa protein. The I4898T mutation was introduced into a rabbit RYR1 cDNA and expressed in HEK-293 cells. The response of the mutant RyR1 Ca2+ channel to the agonists halothane and caffeine in a Ca2+ photometry assay was completely abolished. Coexpression of normal and mutant RYR1 cDNAs in a 1:1 ratio, however, produced RyR1 channels with normal halothane and caffeine sensitivities, but maximal levels of Ca2+ release were reduced by 67%. [3H]Ryanodine binding indicated that the heterozygous channel is activated by Ca2+ concentrations 4-fold lower than normal. Single-cell analysis of cotransfected cells showed a significantly increased resting cytoplasmic Ca2+ level and a significantly reduced luminal Ca2+ level. These data are indicative of a leaky channel, possibly caused by a reduction in the Ca2+ concentration required for channel activation. Comparison with two other coexpressed mutant/normal channels suggests that the I4898T mutation produces one of the most abnormal RyR1 channels yet investigated, and this level of abnormality is reflected in the severe and penetrant phenotype of affected central core disease individuals.

Membrane-bound state of the colicin E1 channel domain as an extended two-dimensional helical array

Atomic level structures have been determined for the soluble forms of several colicins and toxins, but the structural changes that occur after membrane binding have not been well characterized. Changes occurring in the transition from the soluble to membrane-bound state of the C-terminal 190-residue channel polypeptide of colicin E1 (P190) bound to anionic membranes are described. In the membrane-bound state, the α-helical content increases from 60–64% to 80–90%, with a concomitant increase in the average length of the helical segments from 12 to 16 or 17 residues, close to the length required to span the membrane bilayer in the open channel state. The average distance between helical segments is increased and interhelix interactions are weakened, as shown by a major loss of tertiary structure interactions, decreased efficiency of fluorescence resonance energy transfer from an energy donor on helix V of P190 to an acceptor on helix IX, and decreased resonance energy transfer at higher temperatures, not observed in soluble P190, implying freedom of motion of helical segments. Weaker interactions are also shown by a calorimetric thermal transition of low cooperativity, and the extended nature of the helical array is shown by a 3- to 4-fold increase in the average area subtended per molecule to 4,200 Å2 on the membrane surface. The latter, with analysis of the heat capacity changes, implies the absence of a developed hydrophobic core in the membrane-bound P190. The membrane interfacial layer thus serves to promote formation of a highly helical extended two-dimensional flexible net. The properties of the membrane-bound state of the colicin channel domain (i.e., hydrophobic anchor, lengthened and loosely coupled α-helices, and close association with the membrane interfacial layer) are plausible structural features for the state that is a prerequisite for voltage gating, formation of transmembrane helices, and channel opening.

Additional evidence for an eight-transmembrane-domain topology for Caenorhabditis elegans and human presenilins

Presenilins have been implicated in the genesis of Alzheimer’s disease and in facilitating LIN-12/Notch activity during development. All presenilins have multiple hydrophobic regions that could theoretically span a membrane, and a description of the membrane topology is a crucial step toward deducing the mechanism of presenilin function. Previously, we proposed an eight-transmembrane-domain model for presenilin, based on studies of the Caenorhabditis elegans SEL-12 presenilin. Here, we describe experiments that support the view that two of the hydrophobic regions of SEL-12 function as the seventh and eighth transmembrane domains. Furthermore, we have shown that human presenilin 1 behaves like SEL-12 presenilin when analyzed by our methods. Our results provide additional experimental support for the eight-transmembrane-domain model of presenilin topology.