Parallel helix bundles and ion channels: molecular modeling via simulated annealing and restrained molecular dynamics.

A parallel bundle of transmembrane (TM) alpha-helices surrounding a central pore is present in several classes of ion channel, including the nicotinic acetylcholine receptor (nAChR). We have modeled bundles of hydrophobic and of amphipathic helices using simulated annealing via restrained molecular dynamics. Bundles of Ala20 helices, with N = 4, 5, or 6 helices/bundle were generated. For all three N values the helices formed left-handed coiled coils, with pitches ranging from 160 A (N = 4) to 240 A (N = 6). Pore radius profiles revealed constrictions at residues 3, 6, 10, 13, and 17. A left-handed coiled coil and a similar pattern of pore constrictions were observed for N = 5 bundles of Leu20. In contrast, N = 5 bundles of Ile20 formed right-handed coiled coils, reflecting loosened packing of helices containing beta-branched side chains. Bundles formed by each of two classes of amphipathic helices were examined: (a) M2a, M2b, and M2c derived from sequences of M2 helices of nAChR; and (b) (LSSLLSL)3, a synthetic channel-forming peptide. Both classes of amphipathic helix formed left-handed coiled coils. For (LSSLLSL)3 the pitch of the coil increased as N increased from 4 to 6. The M2c N = 5 helix bundle is discussed in the context of possible models of the pore domain of nAChR.     

Synthetic model proteins: the relative contribution of leucine residues at the nonequivalent positions of the 3-4 hydrophobic repeat to the stability of the two-stranded alpha-helical coiled-coil.

Our de novo designed coiled-coil model protein consists of two identical 35-residue polypeptide chains arranged in a parallel and in-register alignment via interchain hydrophobic interactions and a disulfide bridge at the position 2 between two helices. To quantitate the relative contribution of leucine residues at the nonequivalent position of the 3-4 hydrophobic repeat to the stability of the two-stranded alpha-helical coiled-coil, a single alanine was systematically substituted for a leucine in each chain at position "a" (9, 16, 23, or 30) or "d" (5, 12, 19, 26, or 33). The formation and stability of the coiled-coils were determined by circular dichroism studies in the absence and presence of guanidine hydrochloride. All the proteins with an alanine substituted at position a have a similar stability ([Gdn.HCl]1/2 ranges from 2.6 to 2.9 M), while all the proteins with an alanine substituted at position d have similar stability ([Gdn.HCl]1/2 ranges from 3.6 to 4.2 M), except for the proteins with an alanine substituted in the C-terminal heptad. The greater decrease in stability observed for a Leu----Ala mutation at position a (the average delta delta Gu value is 3.3 kcal/mol) compared to those where the substitution was effected at position d (the average delta delta Gu value is 2.0 kcal/mol) indicates that an Ala mutation at position a has a greater effect on the side-chain packing and hydrophobic interactions in the coiled-coil than an Ala mutation at position d.(ABSTRACT TRUNCATED AT 250 WORDS)     

The two-stranded alpha-helical coiled-coil is an ideal model for studying protein stability and subunit interactions.

We have designed de novo a two-stranded alpha-helical coiled-coil which consists of two identical 35-residue polypeptide chains arranged in a parallel and in-register alignment. Their structure is stabilized by interchain hydrophobic interactions from hydrophobes at positions "a" and "d" of a repeating heptad sequence. The formation and stability of the coiled-coil is dependent on peptide concentration due to the monomer-dimer equilibrium. In contrast, that coiled-coil containing an inter-helical disulfide bond does not show any concentration dependence in the guanidine hydrochloride denaturation experiments as expected. Replacement of one large hydrophobic Leu residue in each chain with Ala significantly decreases coiled-coil stability in both the reduced and oxidized coiled-coils [decreases in transition midpoint of 1.6M (2.3-0.7) and 2.4M (5.3-2.9), respectively]. A large pH dependence on coiled-coil stability is observed over the pH range 4 to 7 (transition midpoints at pH 4, 5, 5.5, 6 and 7 were 3.8, 3.2, 2.0, 1.2 and 0.7M, respectively). The increasing stability with decreasing pH correlates with the protonation of the Glu acid side-chains and reduction of intrachain repulsions between Glu-Glu side-chains in positions i, i + 3 or i, i + 4 along each alpha-helix of the coiled-coil. In addition, coiled-coil stability increases with increasing ionic strength.     

Hydrophobic modulation of heme properties in heme protein maquettes

We have investigated the properties of the two hemes bound to histidine in the H10 positions of the uniquely structured apo form of the heme binding four-helix bundle protein maquette [H10H24-L6I,L13F](2), here called [I(6)F(13)H(24)](2) for the amino acids at positions 6 (I), 13 (F) and 24 (H), respectively. The primary structure of each alpha-helix, alpha-SH, in [I(6)F(13)H(24)](2) is Ac-CGGGEI(6)WKL.H(10)EEF(13)LKK.FEELLKL.H(24)EERLKK.L-CONH(2). In our nomenclature, [I(6)F(13)H(24)] represents the disulfide-bridged di-alpha-helical homodimer of this sequence, i.e., (alpha-SS-alpha), and [I(6)F(13)H(24)](2) represents the dimeric four helix bundle composed of two di-alpha-helical subunits, i.e., (alpha-SS-alpha)(2). We replaced the histidines at positions H24 in [I(6)F(13)H(24)](2) with hydrophobic amino acids incompetent for heme ligation. These maquette variants, [I(6)F(13)I(24)](2), [I(6)F(13)A(24)](2), and [I(6)F(13)F(24)](2), are distinguished from the tetraheme binding parent peptide, [I(6)F(13)H(24)](2), by a reduction in the heme:four-helix bundle stoichiometry from 4:1 to 2:1. Iterative redesign has identified phenylalanine as the optimal amino acid replacement for H24 in the context of apo state conformational specificity. Furthermore, the novel second generation diheme [I(6)F(13)F(24)](2) maquette was related to the first generation diheme [H10A24](2) prototype, [L(6)L(13)A(24)](2) in the present nomenclature, via a sequential path in sequence space to evaluate the effects of conservative hydrophobic amino acid changes on heme properties. Each of the disulfide-linked dipeptides studied was highly helical (>77% as determined from circular dichroism spectroscopy), self-associates in solution to form a dimer (as determined by size exclusion chromatography), is thermodynamically stable (-DeltaG(H)2(O) >18 kcal/mol), and possesses conformational specificity that NMR data indicate can vary from multistructured to single structured. Each peptide binds one heme with a dissociation constant, K(d1) value, tighter than 65 nM forming a series of monoheme maquettes. Addition of a second equivalent of heme results in heme binding with a K(d2) in the range of 35-800 nM forming the diheme maquette state. Single conservative amino acid changes between peptide sequences are responsible for up to 10-fold changes in K(d) values. The equilibrium reduction midpoint potential (E(m7.5)) determined in the monoheme state ranges from -156 to -210 mV vs SHE and in the diheme state ranges from -144 to -288 mV. An observed heme-heme electrostatic interaction (>70 mV) in the diheme state indicates a syn global topology of the di-alpha-helical monomers. The heme affinity and electrochemistry of the three H24 variants studied identify the tight binding sites (K(d1) and K(d2) values <200 nM) having the lower reduction midpoint potentials (E(m7.5) values of -155 and -260 mV) with the H10 bound hemes in the parent tetraheme state of [H10H24-L6I,L13F](2), here called [I(6)F(13)H(24)](2). The results of this study illustrate that conservative hydrophobic amino acid changes near the heme binding site can modulate the E(m) by up to +/-50 mV and the K(d) by an order of magnitude. Furthermore, the effects of multiple single amino acid changes on E(m) and K(d) do not appear to be additive.     

Defining the minimum size of a hydrophobic cluster in two-stranded alpha-helical coiled-coils: effects on protein stability

The alpha-helical coiled-coil motif is characterized by a heptad repeat pattern (abcdefg)(n) in which residues a and d form the hydrophobic core. Long coiled-coils (e.g., tropomyosin, 284 residues per polypeptide chain) typically do not have a continuous hydrophobic core of stabilizing residues, but rather one that consists of alternating clusters of stabilizing and destabilizing residues. We have arbitrarily defined a cluster as a minimum of three consecutive stabilizing or destabilizing residues in the hydrophobic core. We report here on a series of two-stranded, disulfide-bridged parallel alpha-helical coiled-coils that contain a central cassette of three consecutive hydrophobic core positions (d, a, and d) with a destabilizing cluster of three consecutive Ala residues in the hydrophobic core on each side of the cassette. The effect of adding one to three stabilizing hydrophobes in these positions (Leu or Ile; denoted as [see text]) was investigated. Alanine residues (denoted as [see text]) are used to represent destabilizing residues. The peptide with three Ala residues in the d a d cassette positions ([see text]) was among the least stable coiled-coil (T(m) = 39.3 degrees C and Urea(1/2) = 1.9 M). Surprisingly, the addition of one stabilizing hydrophobe (Leu) to the cassette or two stabilizing hydrophobes (Leu), still interspersed by an Ala in the cassette ([see text]), also did not lead to any gain in stability. However, peptides with two adjacent hydrophobes in the cassette ([see text])([see text]) did show a gain in stability of 0.9 kcal/mole over the peptide with two interspersed hydrophobes ([see text]). Because the latter three peptides have the same inherent hydrophobicity, the juxtaposition of stabilizing hydrophobes leads to a synergistic effect, and thus a clustering effect. The addition of a third stabilizing hydrophobe to the cassette ([see text]) resulted in a further synergistic gain in stability of 1.7 kcal/mole (T(m) = 54.1 degrees C and Urea(1/2) = 3.3M). Therefore, the role of hydrophobicity in the hydrophobic core of coiled-coils is extremely context dependent and clustering is an important aspect of protein folding and stability.     

Transmembrane helix packing of ErbB/Neu receptor in membrane environment: a molecular dynamics study

Dimerization or oligomerization of the ErbB/Neu receptors are necessary but not sufficient for initiation of receptor signaling. The two intracellular domains must be properly oriented for the juxtaposition of the kinase domains allowing trans-phosphorylation. This suggests that the transmembrane (TM) domain acts as a guide for defining the proper orientation of the intracellular domains. Two structural models, with the two helices either in left-handed or in right-handed coiling have been proposed as the TM domain structure of the active receptor. Because experimental data do not distinguish clearly helix-helix packing, molecular dynamics (MD) simulations are used to investigate the energetic factors that drive Neu TM-TM interactions of the wild and the oncogenic receptor (Val664/Glu mutation) in DMPC or in POPC environments. MD results indicate that helix-lipid interactions in the bilayer core are extremely similar in the two environments and raise the role of the juxtamembrane residues in helix insertion and helix-helix packing. The TM domain shows a greater propensity to adopt a left-handed structure in DMPC, with helices in optimal position for strong inter-helical Hbonds induced by the Glu mutation. In POPC, the right-handed structure is preferentially formed with the participation of water in inter-helical Hbonds. The two structural arrangements of the Neu(TM) helices both with GG4 residue motif in close contact at the interface are permissible in the membrane environment. According to the hypothesis of a monomer-dimer equilibrium of the proteins it is likely that the bilayer imposes structural constraints that favor dimerization-competent structure responsible of the proper topology necessary for receptor activation.     

Design of amphiphilic protein maquettes: controlling assembly, membrane insertion, and cofactor interactions

We have designed polypeptides combining selected lipophilic (LP) and hydrophilic (HP) sequences that assemble into amphiphilic (AP) alpha-helical bundles to reproduce key structure characteristics and functional elements of natural membrane proteins. The principal AP maquette (AP1) developed here joins 14 residues of a heme binding sequence from a structured diheme-four-alpha-helical bundle (HP1), with 24 residues of a membrane-spanning LP domain from the natural four-alpha-helical M2 channel of the influenza virus, through a flexible linking sequence (GGNG) to make a 42 amino acid peptide. The individual AP1 helices (without connecting loops) assemble in detergent into four-alpha-helical bundles as observed by analytical ultracentrifugation. The helices are oriented parallel as indicated by interactions typical of adjacent hemes. AP1 orients vectorially at nonpolar-polar interfaces and readily incorporates into phospholipid vesicles with >97% efficiency, although most probably without vectorial bias. Mono- and diheme-AP1 in membranes enhance functional elements well established in related HP analogues. These include strong redox charge coupling of heme with interior glutamates and internal electric field effects eliciting a remarkable 160 mV splitting of the redox potentials of adjacent hemes that leads to differential heme binding affinities. The AP maquette variants, AP2 and AP3, removed heme-ligating histidines from the HP domain and included heme-ligating histidines in LP domains by selecting the b(H) heme binding sequence from the membrane-spanning d-helix of respiratory cytochrome bc(1). These represent the first examples of AP maquettes with heme and bacteriochlorophyll binding sites located within the LP domains.     

Crystal structure of recombinant human interleukin-4.

The crystal structure of recombinant human interleukin-4 (rhuIL-4) was initially determined at 3.5-A resolution by multiple isomorphous replacement techniques and subsequently refined to a resolution of 2.35 A by simulated annealing. The final crystallographic R-factor, based on all data in the range 6.0-2.35 A (7470 reflections), is 0.232. Bond lengths and bond angles in the molecule have root mean square deviations from ideal values of 0.016 A and 2.4 degrees, respectively. The overall structure is highly compact and globular with a predominantly hydrophobic core. The main structural feature of rhuIL-4 is a four alpha-helix bundle, which composes approximately 58% of the structure. The helices are arranged in a left-handed antiparallel bundle with two overhand connections. Within these connections is a two-stranded antiparallel beta-sheet. Both the tertiary and secondary structures of rhuIL-4 are similar to those of human granulocyte-macrophage colony-stimulating factor. Critical regions for receptor binding are proposed.     

Four-alpha-helix bundle with designed anesthetic binding pockets. Part I: structural and dynamical analyses

The four-α-helix bundle mimics the transmembrane domain of the Cys-loop receptor family believed to be the protein target for general anesthetics. Using high resolution NMR, we solved the structure (Protein Data Bank ID: 2I7U) of a prototypical dimeric four-α-helix bundle, (Aα₂-L1M/L38M)₂, with designed specific binding pockets for volatile anesthetics. Two monomers of the helix-turn-helix motif form an antiparallel dimer as originally designed, but the high-resolution structure exhibits an asymmetric quaternary arrangement of the four helices. The two helices from the N-terminus to the linker (helices 1 and 1′) are associated with each other in the dimer by the side-chain ring stacking of F12 and W15 along the long hydrophobic core and by a nearly perfect stretch of hydrophobic interactions between the complementary pairs of L4, L11, L18, and L25, all of which are located at the heptad e position along the helix-helix dimer interface. In comparison, the axes of the two helices from the linker to the C-terminus (helices 2 and 2′) are wider apart from each other, creating a lateral access pathway around K47 from the aqueous phase to the center of the designed hydrophobic core. The site of the L38M mutation, which was previously shown to increase the halothane binding affinity by ∼3.5-fold, is not part of the hydrophobic core presumably involved in the anesthetic binding but shows an elevated transverse relaxation (R₂) rate. Qualitative analysis of the protein dynamics by reduced spectral density mapping revealed exchange contributions to the relaxation at many residues in the helices. This observation was confirmed by the quantitative analysis using the Modelfree approach and by the NMR relaxation dispersion measurements. The NMR structures and Autodock analysis suggest that the pocket with the most favorable amphipathic property for anesthetic binding is located between the W15 side chains at the center of the dimeric hydrophobic core, with the possibility of two additional minor binding sites between the F12 and F52 ring stacks of each monomer. The high-resolution structure of the designed anesthetic-binding protein offers unprecedented atomistic details about possible sites for anesthetic-protein interactions that are essential to the understanding of molecular mechanisms of general anesthesia.     

Comparison of helix interactions in membrane and soluble alpha-bundle proteins

Helix-helix interactions are important for the folding, stability, and function of membrane proteins. Here, two independent and complementary methods are used to investigate the nature and distribution of amino acids that mediate helix-helix interactions in membrane and soluble alpha-bundle proteins. The first method characterizes the packing density of individual amino acids in helical proteins based on the van der Waals surface area occluded by surrounding atoms. We have recently used this method to show that transmembrane helices pack more tightly, on average, than helices in soluble proteins. These studies are extended here to characterize the packing of interfacial and noninterfacial amino acids and the packing of amino acids in the interfaces of helices that have either right- or left-handed crossing angles, and either parallel or antiparallel orientations. We show that the most abundant tightly packed interfacial residues in membrane proteins are Gly, Ala, and Ser, and that helices with left-handed crossing angles are more tightly packed on average than helices with right-handed crossing angles. The second method used to characterize helix-helix interactions involves the use of helix contact plots. We find that helices in membrane proteins exhibit a broader distribution of interhelical contacts than helices in soluble proteins. Both helical membrane and soluble proteins make use of a general motif for helix interactions that relies mainly on four residues (Leu, Ala, Ile, Val) to mediate helix interactions in a fashion characteristic of left-handed helical coiled coils. However, a second motif for mediating helix interactions is revealed by the high occurrence and high average packing values of small and polar residues (Ala, Gly, Ser, Thr) in the helix interfaces of membrane proteins. Finally, we show that there is a strong linear correlation between the occurrence of residues in helix-helix interfaces and their packing values, and discuss these results with respect to membrane protein structure prediction and membrane protein stability.     

Heme redox potential control in de novo designed four-alpha-helix bundle proteins

The effects of various mechanisms of metalloporphyrin reduction potential modulation were investigated experimentally using a robust, well-characterized heme protein maquette, synthetic protein scaffold H10A24 [?CH(3)()CONH-CGGGELWKL.HEELLKK.FEELLKL.AEERLKK. L-CONH(2)()?(2)](2). Removal of the iron porphyrin macrocycle from the high dielectric aqueous environment and sequestration within the hydrophobic core of the H10A24 maquette raises the equilibrium reduction midpoint potential by 36-138 mV depending on the hydrophobicity of the metalloporphyrin structure. By incorporating various natural and synthetic metalloporphyrins into a single protein scaffold, we demonstrate a 300-mV range in reduction potential modulation due to the electron-donating/withdrawing character of the peripheral macrocycle substituents. Solution pH is used to modulate the metalloporphyrin reduction potential by 160 mV, regardless of the macrocycle architecture, by controlling the protonation state of the glutamate involved in partial charge compensation of the ferric heme. Attempts to control the reduction potential by inserting charged amino acids into the hydrophobic core at close proximity to the metalloporphyrin lead to varied success, with H10A24-L13E lowering the E(m8.5) by 40 mV, H10A24-E11Q raising it by 50 mV, and H10A24-L13R remaining surprisingly unaltered. Modifying the charge of the adjacent metalloporphyrin, +1 for iron(III) protoporphyrin IX or neutral for zinc(II) protoporphyrin IX resulted in a loss of 70 mV [Fe(III)PPIX](+) - [Fe(III)PPIX](+) interaction observed in maquettes. Using these factors in combination, we illustrate a 435-mV variation of the metalloporphyrin reduction midpoint potential in a simple heme maquette relative to the about 800-mV range observed for natural cytochromes. Comparison between the reduction potentials of the heme maquettes and other de novo designed heme proteins reveals global trends in the E(m) values of synthetic cytochromes.     

The structure of ColE1 rop in solution

Summary The structure of the ColE1 repressor of primer (rop) protein in solution was determined from the proton nuclear magnetic resonance data by a combined use of distance geometry and restrained molecular dynamics calculations. A set of structures was determined with low internal energy and virtually no violations of the experimental distance restraints. Rop forms homodimers: Two helical hairpins are arranged as an antiparallel four helix bundle with a left-handed rope-like twist of the helix axes and with left-handed bundle topology. The very compact packing of the side chains in the helix interfaces of the rop coiled-coil structure may well account for its high stability. Overall, the solution structure is highly similar to the recently determined X-ray structure (Banner, D.W., Kokkinidis, M. and Tsernoglou, D. (1987)J. Mol. Biol.,196, 657–675), although there are minor differences in regions where packing forces appear to influence the crystal structure.Abbreviations rop repressor of primer - NMR nuclear magnetic resonance - NOE nuclear Overhauser enhancement - NOESY NOE spectroscopy - RAN Set Structures generated from random choice of the dihedrai angles - HEL Set Structures generated from random choice of the dihedral angles restricted to ranges allowed for helices - MD molecular dynamics - EM energy minimization - RMSD root-mean-square deviation of atomic positions     

Clustering of large hydrophobes in the hydrophobic core of two-stranded alpha-helical coiled-coils controls protein folding and stability

The de novo design and biophysical characterization of two 60-residue peptides that dimerize to fold as parallel coiled-coils with different hydrophobic core clustering is described. Our goal was to investigate whether designing coiled-coils with identical hydrophobicity but with different hydrophobic clustering of non-polar core residues (each contained 6 Leu, 3 Ile, and 7 Ala residues in the hydrophobic core) would affect helical content and protein stability. The disulfide-bridged P3 and P2 differed dramatically in alpha-helical structure in benign conditions. P3 with three hydrophobic clusters was 98% alpha-helical, whereas P2 was only 65% alpha-helical. The stability profiles of these two analogs were compared, and the enthalpy and heat capacity changes upon denaturation were determined by measuring the temperature dependence by circular dichroism spectroscopy and confirmed by differential scanning calorimetry. The results showed that P3 assembled into a stable alpha-helical two-stranded coiled-coil and exhibited a native protein-like cooperative two-state transition in thermal melting, chemical denaturation, and calorimetry experiments. Although both peptides have identical inherent hydrophobicity (the hydrophobic burial of identical non-polar residues in equivalent heptad coiled-coil positions), we found that the context dependence of an additional hydrophobic cluster dramatically increased stability of P3 (Delta Tm approximately equal to 18 degrees C and Delta[urea](1/2) approximately equal to 1.5 M) as compared with P2. These results suggested that hydrophobic clustering significantly stabilized the coiled-coil structure and may explain how long fibrous proteins like tropomyosin maintain chain integrity while accommodating polar or charged residues in regions of the protein hydrophobic core.     

Loopless Rop: structure and dynamics of an engineered homotetrameric variant of the repressor of primer protein

The repressor of primer (Rop) protein has become a steady source of surprises concerning the relationship between the sequences and the structures of several of its mutants and variants. Here we add another piece to the puzzle of Rop by showing that an engineered deletion mutant of the protein (corresponding to a deletion of residues 30-34 of the wild-type protein and designed to restore the heptad periodicity at the turn region) results in a complete reorganization of the bundle which is converted from a homodimer to a homotetramer. In contrast (and as previously shown), a two-residue insertion, which also restores the heptad periodicity, is essentially identical with wild-type Rop. The new deletion mutant structure is a canonical, left-handed, all-antiparallel bundle with a completely different hydrophobic core and distinct surface properties. The structure agrees and qualitatively explains the results from functional, thermodynamic, and kinetic studies which indicated that this deletion mutant is a biologically inactive hyperstable homotetramer. Additional insight into the stability and dynamics of the mutant structure has been obtained from extensive molecular dynamics simulations in explicit water and with full treatment of electrostatics.     

Fusion core structure of the severe acute respiratory syndrome coronavirus (SARS-CoV): in search of potent SARS-CoV entry inhibitors

Severe acute respiratory coronavirus (SARS-CoV) spike (S) glycoprotein fusion core consists of a six-helix bundle with the three C-terminal heptad repeat (HR2) helices packed against a central coiled-coil of the other three N-terminal heptad repeat (HR1) helices. Each of the three peripheral HR2 helices shows prominent contacts with the hydrophobic surface of the central HR1 coiled-coil. The concerted protein-protein interactions among the HR helices are responsible for the fusion event that leads to the release of the SARS-CoV nucleocapsid into the target host-cell. In this investigation, we applied recombinant protein and synthetic peptide-based biophysical assays to characterize the biological activities of the HR helices. In a parallel experiment, we employed a HIV-luc/SARS pseudotyped virus entry inhibition assay to screen for potent inhibitory activities on HR peptides derived from the SARS-CoV S protein HR regions and a series of other small-molecule drugs. Three HR peptides and five small-molecule drugs were identified as potential inhibitors. ADS-J1, which has been used to interfere with the fusogenesis of HIV-1 onto CD4+ cells, demonstrated the highest HIV-luc/SARS pseudotyped virus-entry inhibition activity among the other small-molecule drugs. Molecular modeling analysis suggested that ADS-J1 may bind to the deep pocket of the hydrophobic groove on the surface of the central coiled-coil of SARS-CoV S HR protein and prevent the entrance of the SARS-CoV into the host cells.     

Transmembrane Helix Packing of ErbB/Neu Receptor in Membrane Environment: A Molecular Dynamics Study

Abstract

Dimerization or oligomerization of the ErbB/Neu receptors are necessary but not sufficient for initiation of receptor signaling. The two intracellular domains must be properly oriented for the juxtaposition of the kinase domains allowing trans-phosphorylation. This suggests that the transmembrane (TM) domain acts as a guide for defining the proper orientation of the intracellular domains.

Two structural models, with the two helices either in left-handed or in right-handed coiling have been proposed as the TM domain structure of the active receptor. Because experimental data do not distinguish clearly helix-helix packing, molecular dynamics (MD) simulations are used to investigate the energetic factors that drive Neu TM-TM interactions of the wild and the oncogenic receptor (Val664/Glu mutation) in DMPC or in POPC environments. MD results indicate that helix-lipid interactions in the bilayer core are extremely similar in the two environments and raise the role of the juxtamembrane residues in helix insertion and helix-helix packing. The TM domain shows a greater propensity to adopt a left-handed structure in DMPC, with helices in optimal position for strong inter-helical Hbonds induced by the Glu mutation. In POPC, the right-handed structure is preferentially formed with the participation of water in inter-helical Hbonds. The two structural arrangements of the NeuTM helices both with GG4 residue motif in close contact at the interface are permissible in the membrane environment. According to the hypothesis of a monomer-dimer equilibrium of the proteins it is likely that the bilayer imposes structural constraints that favor dimerization- competent structure responsible of the proper topology necessary for receptor activation.     

Protein destabilization by electrostatic repulsions in the two-stranded alpha-helical coiled-coil/leucine zipper.

The destabilizing effect of electrostatic repulsions on protein stability has been studied by using synthetic two-stranded alpha-helical coiled-coils as a model system. The native coiled-coil consists of two identical 35-residue polypeptide chains with a heptad repeat QgVaGbAcLdQeKf and a Cys residue at position 2 to allow formation of an interchain disulfide bridge. This peptide, designed to contain no intrahelical or interhelical electrostatic interactions, forms a stable coiled-coil structure at 20 degrees C in benign medium (50 mM KCl, 25 mM PO4, pH 7) with a [urea]1/2 value of 6.1 M. Four mutant coiled-coils were designed to contain one or two Glu substitutions for Gln per polypeptide chain. The resulting coiled-coils contained potential i to i' + 5 Glu-Glu interchain repulsions (denoted as peptide E2(15,20)), i to i' + 2 Glu-Glu interchain repulsions (denoted E2(20,22)), or no interchain ionic interactions (denoted E2(13,22) and E1(20)). The stabilities of the coiled-coils were determined by measuring the ellipticities at 222 nm as a function of urea or guanidine hydrochloride concentration at 20 degrees C in the presence and absence of an interchain disulfide bridge. At pH 7, in the presence of urea, the stabilities of E2(13,22) and E2(20,22) were identical suggesting that the potential i to i' + 2 interchain Glu-Glu repulsion in the E2(20,22) coiled-coil does not occur. In contrast, the mutant E2(15,20) is substantially less stable than E2(13,22) or E2(15,20) by 0.9 kcal/mol due to the presence of two i to i' + 5 interchain Glu-Glu repulsions, which destabilize the coiled-coil by 0.45 kcal/mol each. At pH 3 the coiled-coils were found to increase in stability as the number of Glu substitutions were increased. This, combined with reversed-phase HPLC results at pH 7 and pH 2, supports the conclusion that the protonated Glu side chains present at low pH are significantly more hydrophobic than Gln side chains which are in turn more hydrophobic than the ionized Glu side chains present at neutral pH. The protonated Glu residues increase the hydrophobicity of the coiled-coil interface leading to higher coiled-coil stability. The guanidine hydrochloride results at pH 7 show similar stabilities between the native and mutant coiled-coils indicating that guanidine hydrochloride masks electrostatic repulsions due to its ionic nature and that Glu and Gln in the e and g positions of the heptad repeat have very similar effects on coiled-coil stability in the presence of GdnHCl.     

Molecular dynamics simulation of the M2 helices within the nicotinic acetylcholine receptor transmembrane domain: structure and collective motions

Multiple nanosecond duration molecular dynamics simulations were performed on the transmembrane region of the Torpedo nicotinic acetylcholine receptor embedded within a bilayer mimetic octane slab. The M2 helices and M2-M3 loop regions were free to move, whereas the outer (M1, M3, M4) helix bundle was backbone restrained. The M2 helices largely retain their hydrogen-bonding pattern throughout the simulation, with some distortions in the helical end and loop regions. All of the M2 helices exhibit bending motions, with the hinge point in the vicinity of the central hydrophobic gate region (corresponding to residues alphaL251 and alphaV255). The bending motions of the M2 helices lead to a degree of dynamic narrowing of the pore in the region of the proposed hydrophobic gate. Calculations of Born energy profiles for various structures along the simulation trajectory suggest that the conformations of the M2 bundle sampled correspond to a closed conformation of the channel. Principal components analyses of each of the M2 helices, and of the five-helix M2 bundle, reveal concerted motions that may be relevant to channel function. Normal mode analyses using the anisotropic network model reveal collective motions similar to those identified by principal components analyses.