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.

Alpha-helical, but not beta-sheet, propensity of proline is determined by peptide environment.

Proline is established as a potent breaker of both alpha-helical and beta-sheet structures in soluble (globular) proteins. Thus, the frequent occurrence of the Pro residue in the putative transmembrane helices of integral membrane proteins, particularly transport proteins, presents a structural dilemma. We propose that this phenomenon results from the fact that the structural propensity of a given amino acid may be altered to conform to changes imposed by molecular environment. To test this hypothesis on proline, we synthesized model peptides of generic sequence H2N-(Ser-LyS)2-Ala- Leu-Z-Ala-Leu-Z-Trp-Ala-Leu-Z-(Lys-Ser)3-OH (Z = Ala and/or Pro). Peptide conformations were analyzed by circular dichroism spectroscopy in aqueous buffer, SDS, lysophosphatidylglycerol micelles, and organic solvents (methanol, trifluoroethanol, and 2-propanol). The helical propensity of Pro was found to be greatly enhanced in the membrane-mimetic environments of both lipid micelles and organic solvents. Proline was found to stabilize the alpha-helical conformation relative to Ala at elevated temperatures in 2-propanol, an observation that argues against the doctrine that Pro is the most potent alpha-helix breaker as established in aqueous media. Parallel studies in deoxycholate micelles of the temperature-induced conformational transitions of the single-spanning membrane bacteriophage IKe major coat protein, in which the Pro-containing wild type was compared with Pro30 --> Ala mutant, Pro was found to protect the helix, but disrupt the beta-sheet structure as effectively as it does to model peptides in water. The intrinsic capacity of Pro to disrupt beta-sheets was further reflected in a survey of porins where Pro was found to be selectively excluded from the core of membrane-spanning beta-sheet barrels. The overall data provide a rationale for predicting and understanding the structural consequences when Pro occurs in the context of a membrane.

Thermodynamic parameters for the helix-coil transition of oligopeptides: molecular dynamics simulation with the peptide growth method.

The helix-coil transition equilibrium of polypeptides in aqueous solution was studied by molecular dynamics simulation. The peptide growth simulation method was introduced to generate dynamic models of polypeptide chains in a statistical (random) coil or an alpha-helical conformation. The key element of this method is to build up a polypeptide chain during the course of a molecular transformation simulation, successively adding whole amino acid residues to the chain in a predefined conformation state (e.g., alpha-helical or statistical coil). Thus, oligopeptides of the same length and composition, but having different conformations, can be incrementally grown from a common precursor, and their relative conformational free energies can be calculated as the difference between the free energies for growing the individual peptides. This affords a straightforward calculation of the Zimm-Bragg sigma and s parameters for helix initiation and helix growth. The calculated sigma and s parameters for the polyalanine alpha-helix are in reasonable agreement with the experimental measurements. The peptide growth simulation method is an effective way to study quantitatively the thermodynamics of local protein folding.

Interaction of batrachotoxin with the local anesthetic receptor site in transmembrane segment IVS6 of the voltage-gated sodium channel

The voltage-gated sodium channel is the site of action of more than six classes of neurotoxins and drugs that alter its function by interaction with distinct, allosterically coupled receptor sites. Batrachotoxin (BTX) is a steroidal alkaloid that binds to neurotoxin receptor site 2 and causes persistent activation. BTX binding is inhibited allosterically by local anesthetics. We have investigated the interaction of BTX with amino acid residues I1760, F1764, and Y1771, which form part of local anesthetic receptor site in transmembrane segment IVS6 of type IIA sodium channels. Alanine substitution for F1764 (mutant F1764A) reduces tritiated BTX-A-20-α-benzoate binding affinity, causing a 60-fold increase in Kd. Alanine substitution for I1760, which is adjacent to F1764 in the predicted IVS6 transmembrane alpha helix, causes only a 4-fold increase in Kd. In contrast, mutant Y1771A shows no change in BTX binding affinity. For wild-type and mutant Y1771A, BTX shifted the voltage for half-maximal activation ≈40 mV in the hyperpolarizing direction and increased the percentage of noninactivating sodium current to ≈60%. In contrast, these BTX effects were eliminated completely for the F1764A mutant and were reduced substantially for mutant I1760A. Our data suggest that the BTX receptor site shares overlapping but nonidentical molecular determinants with the local anesthetic receptor site in transmembrane segment IVS6 as well as having unique molecular determinants in transmembrane segment IS6, as demonstrated in previous work. Evidently, BTX conforms to a domain–interface allosteric model of ligand binding and action, as previously proposed for calcium agonist and antagonist drugs acting on l-type calcium channels.

Redesigning secondary structure to invert coenzyme specificity in isopropylmalate dehydrogenase.

Rational engineering of enzymes involves introducing key amino acids guided by a knowledge of protein structure to effect a desirable change in function. To date, all successful attempts to change specificity have been limited to substituting individual amino acids within a protein fold. However, the infant field of protein engineering will only reach maturity when changes in function can be generated by rationally engineering secondary structures. Guided by x-ray crystal structures and molecular modeling, site-directed mutagenesis has been used to systematically invert the coenzyme specificity of Thermus thermophilus isopropylmalate dehydrogenase from a 100-fold preference for NAD to a 1000-fold preference for NADP. The engineered mutant, which is twice as active as wild type, contains four amino acid substitutions and an alpha-helix and loop that replaces the original beta-turn. These results demonstrate that rational engineering of secondary structures to produce enzymes with novel properties is feasible.

Protein-protein interactions in the rigor actomyosin complex.

Since it has not been possible to crystallize the actomyosin complex, the x-ray structures of the individual proteins together with data obtained by fiber diffraction and electron microscopy have been used to build detailed models of filamentous actin (f-actin) and the actomyosin rigor complex. In the f-actin model, a single monomer uses 10 surface loops and two alpha-helices to make sometimes complicated interactions with its four neighbors. In the myosin molecule, both the essential and regulatory light chains show considerable structural homology to calmodulin. General principles are evident in their mode of attachment to the target alpha-helix of the myosin heavy chain. The essential light chain also makes contacts with other parts of the heavy chain and with the regulatory light chain. The actomyosin rigor interface is extensive, involving interaction of a single myosin head with regions on two adjacent actin monomers. A number of hydrophobic residues on the apposing faces of actin and myosin contribute to the main binding site. This site is flanked on three sides by charged myosin surface loops that form predominantly ionic interactions with adjacent regions of actin. Hydrogen bonding is likely to play a significant role in actin-actin and actin-myosin interactions since many of the contacts involve loops. The model building approach used with actomyosin is applicable to other multicomponent assemblies of biological interest and is a powerful method for revealing molecular interactions and providing insights into the mode of action of the assemblies.

Three distinct structural environments of a transmembrane domain in the inwardly rectifying potassium channel ROMK1 defined by perturbation.

To probe the protein environment of an ion channel, we have perturbed the structure of a transmembrane domain by substituting side chains with those of two different sizes by using site-specific mutagenesis. We have used Trp and Ala as a high- and a low-impact perturbation probe, respectively, to replace each of 18 consecutive residues within the putative second transmembrane segment, M2, of an inwardly rectifying potassium channel, ROMK1. Our rationale is that a change in the channel function as a consequence of these mutations at a particular position will reflect the structural environment of the altered side chain. Each position can then be assigned to one of three classes of environments, as grated by different levels of perturbation: very tolerant (channel functions with both Trp and Ala substitutions), tolerant (function preserved with Ala but not with Trp substitution), and intolerant (either Ala or Trp substitution destroys function). We identify the very tolerant environment as being lipid-facing, tolerant as protein-interior-facing, and intolerant as pore-facing. We observe a strikingly ordered pattern of perturbation of all three environmental classes. This result indicates that M2 is a straight alpha-helix.

Optimization of rates of protein folding: the nucleation-condensation mechanism and its implications.

Small, single-module proteins that fold in a single cooperative step may be paradigms for understanding early events in protein-folding pathways generally. Recent experimental studies of the 64-residue chymotrypsin inhibitor 2 (CI2) support a nucleation mechanism for folding, as do some computer stimulations. CI2 has a nucleation site that develops only in the transition state for folding. The nucleus is composed of a set of adjacent residues (an alpha-helix), stabilized by long-range interactions that are formed as the rest of the protein collapses around it. A simple analysis of the optimization of the rate of protein folding predicts that rates are highest when the denatured state has little residual structure under physiological conditions and no intermediates accumulate. This implies that any potential nucleation site that is composed mainly of adjacent residues should be just weakly populated in the denatured state and become structured only in a high-energy intermediate or transition state when it is stabilized by interactions elsewhere in the protein. Hierarchical mechanisms of folding in which stable elements of structure accrete are unfavorable. The nucleation-condensation mechanism of CI2 fulfills the criteria for fast folding. On the other hand, stable intermediates do form in the folding of more complex proteins, and this may be an unavoidable consequence of increasing size and nucleation at more than one site.

Lyophilization-induced reversible changes in the secondary structure of proteins.

Changes in the secondary structure of some dozen different proteins upon lyophilization of their aqueous solutions have been investigated by means of Fourier-transform infrared spectroscopy in the amide III band region. Dehydration markedly (but reversibly) alters the secondary structure of all the proteins studied, as revealed by both the quantitative analysis of the second derivative spectra and the Gaussian curve fitting of the original infrared spectra. Lyophilization substantially increases the beta-sheet content and lowers the alpha-helix content of all proteins. In all but one case, proteins become more ordered upon lyophilization.

Specificity of dimer formation in tropomyosins: influence of alternatively spliced exons on homodimer and heterodimer assembly.

Tropomyosins consist of nearly 100% alpha-helix and assemble into parallel and in-register coiled-coil dimers. In vitro it has been established that nonmuscle as well as native muscle tropomyosins can form homodimers. However, a mixture of muscle alpha and beta tropomyosin subunits results in the formation of the thermodynamically more stable alpha/beta heterodimer. Although the assembly preference of the muscle tropomyosin heterodimer can be understood thermodynamically, the presence of multiple tropomyosin isoforms expressed in nonmuscle cells points toward a more complex principle for determining dimer formation. We have investigated the dimerization of rat tropomyosins in living cells by the use of epitope tagging with a 16-aa sequence of the influenza hemagglutinin. Employing transfection and immunoprecipitation techniques, we have analyzed the dimers formed by muscle and nonmuscle tropomyosins in rat fibroblasts. We demonstrate that the information for homo- versus heterodimerization is contained within the tropomyosin molecule itself and that the information for the selectivity is conferred by the alternatively spliced exons. These results have important implications for models of the regulation of cytoskeletal dynamics.

The CuA center of cytochrome-c oxidase: electronic structure and spectra of models compared to the properties of CuA domains.

The electronic structure and spectrum of several models of the binuclear metal site in soluble CuA domains of cytochrome-c oxidase have been calculated by the use of an extended version of the complete neglect of differential overlap/spectroscopic method. The experimental spectra have two strong transitions of nearly equal intensity around 500 nm and a near-IR transition close to 800 nm. The model that best reproduces these features consists of a dimer of two blue (type 1) copper centers, in which each Cu atom replaces the missing imidazole on the other Cu atom. Thus, both Cu atoms have one cysteine sulfur atom and one imidazole nitrogen atom as ligands, and there are no bridging ligands but a direct Cu-Cu bond. According to the calculations, the two strong bands in the visible region originate from exciton coupling of the dipoles of the two copper monomers, and the near-IR band is a charge-transfer transition between the two Cu atoms. The known amino acid sequence has been used to construct a molecular model of the CuA site by the use of a template and energy minimization. In this model, the two ligand cysteine residues are in one turn of an alpha-helix, whereas one ligand histidine is in a loop following this helix and the other one is in a beta-strand.

Residues in the pathway through a membrane transporter.

The structure of solute transporters is understood largely from analysis of their amino acid sequences, and more direct information is greatly needed. Here we report work that applies cysteine scanning mutagenesis to describe structure-function relations in UhpT, a bacterial membrane transporter. By using an impermeant SH-reactive agent to probe single-cysteine variants, we show that UhpT transmembrane segment 7 spans the membrane as an alpha-helix and that the central portion of this helix is exposed to both membrane surfaces, forming part of the translocation pathway through this transporter.

Structural variety of arginine-rich RNA-binding peptides.

Arginine-rich domains are used by a variety of RNA-binding proteins to recognize specific RNA hairpins. It has been shown previously that a 17-aa arginine-rich peptide from the human immunodeficiency virus Rev protein binds specifically to its RNA site when the peptide is in an alpha-helical conformation. Here we show that related peptides from splicing factors, viral coat proteins, and bacteriophage antiterminators (the N proteins) also have propensities to form alpha-helices and that the N peptides require helical conformations to bind to their cognate RNAs. In contrast, introducing proline mutations into the arginine-rich domain of the human immunodeficiency virus Tat protein abolishes its potential to form an alpha-helix but does not affect RNA-binding affinity in vitro or in vivo. Based on results from several peptide-RNA model systems, we suggest that helical peptides may be used to recognize RNA structures having particularly wide major grooves, such as those found near loops or large bulges, and that nonhelical or extended peptides may be used to recognize less accessible grooves.

Specific mutations in the estrogen receptor change the properties of antiestrogens to full agonists.

The estrogen receptor (ER) stimulates transcription of target genes by means of its two transcriptional activation domains, AF-1 in the N-terminal part of the receptor and AF-2 in its ligand-binding domain. AF-2 activity is dependent upon a putative amphipathic alpha-helix between residues 538 and 552 in the mouse ER. Point mutagenesis of conserved hydrophobic residues within this region reduces estrogen-dependent transcriptional activation without affecting hormone and DNA binding significantly. Here we show that these mutations dramatically alter the pharmacology of estrogen antagonists. Both tamoxifen and ICI 164,384 behave as strong agonists in HeLa cells expressing the ER mutants. In contrast to the wild-type ER, the mutant receptors maintain nuclear localization and DNA-binding activity after ICI 164,384 treatment. Structural alterations in AF-2 caused by gene mutations such as those described herein or by estrogen-independent signaling pathways may account for the insensitivity of some breast cancers to tamoxifen treatment.

Folding of a nascent polypeptide chain in vitro: cooperative formation of structure in a protein module.

We have prepared a family of peptide fragments of the 64-residue chymotrypsin inhibitor 2, corresponding to its progressive elongation from the N terminus. The growing polypeptide chain has little tendency to form stable structure until it is largely synthesized, and what structures are formed are nonnative and lack, in particular, the native secondary structural elements of alpha-helix and beta-sheet. These elements then develop as sufficient tertiary interactions are made in the nearly full-length chain. The growth of structure in the small module is highly cooperative and does not result from the hierarchical accretion of substructures.