Backbone dynamics of detergent-solubilized alamethicin from amide hydrogen exchange measurements.

Alamethicin is a 20 amino acid antibiotic peptide produced by the soil fungus Trichoderma viride. The peptide inserts into bacterial membranes and self-associates to form ion channels, but the details of this process are unknown. Residue-specific acid- and base-catalyzed exchange data were obtained for 16 of 18 backbone amides of alamethicin dissolved in sodium dodecyl sulfate micelles using high-resolution 2-dimensional heteronuclear nuclear magnetic resonance spectroscopy. To facilitate interpretation of the exchange data, we synthesized N-acetyl-alpha-aminoisobutyric acid-N'-methyl and N-acetyl-alanine-N'-methyl and measured the pD dependence of their hydrogen-deuterium exchange rates to determine the sequence-dependent inductive and steric effects of the alpha-aminoisobutyric acid residue. Intramolecular H-bonding in alamethicin was monitored through the exchange parameters kmin (minimum exchange rate) which indicate that the backbone is significantly more stable than the backbones of alanine-based helical peptides. Rapid exchange at Gly-11 suggests a highly local conformational flexibility in the middle of the peptide. Interactions with the detergent micelle were revealed by the exchange parameters pDmin (pD of minimum exchange) which suggest that the N-terminus of alamethicin interacts more strongly with the detergent micelle than does the C-terminus. A periodicity in pDmin difference data reveals that one surface of the helix interacts more strongly with the micelle. The surface consists of residues 1, 5, 9, 13, 16, and 20. The opposite face of the helix contains several polar residues (two glutamines and a glycine), suggesting that, on average, this face of the helix is directed toward the solvent. These results serve as a model for the interaction of the peptide with membranes containing anionic lipid. In combination with published molecular dynamics simulations [Gibbs et al. (1997) Biophys. J. 72, 2490-2495], the present results also offer insight into the mechanisms of hydrogen-deuterium exchange in helical peptides.     

Molecular flexibility demonstrated by paramagnetic enhancements of nuclear relaxation. Application to alamethicin: a voltage-gated peptide channel.

A nitroxide spin label attached to the C-terminus of the channel forming peptide alamethicin produces an enhancement of the nuclear spin-lattice relaxation rates of peptide protons as a result of both intermolecular and intramolecular magnetic dipole-dipole interactions. The intermolecular contribution provides evidence that alamethicin monomers collide preferentially in a C-terminal-to-N-terminal configuration in methanol. From the intramolecular paramagnetic enhancement of nuclear spin-lattice relaxation times, effective distances between the unpaired electron on the nitroxide at the C-terminus of alamethicin and protons along the peptide backbone were calculated. These distances are much shorter than distances based on the reported crystal structure of alamethicin, and cannot be accounted for by motion in the bonds that attach the nitroxide to the peptide. In addition, the differences between distances deduced from the nuclear spin relaxation and the distances seen in the crystal structure increase toward the N-terminal end of the peptide. The simplest explanation for these data is that the alamethicin backbone suffers large structural fluctuations that yield shorter effective distances between the C-terminus and positions along the backbone. This finding can be interpreted in terms of a molecular mechanism for the voltage-gating of the alamethicin channel. When the distances between a paramagnetic center and a nucleus fluctuate, paramagnetic enhancements are expected to yield distances that are weighted by r-6, and distances calculated using the Solomon-Bloembergen equations may more nearly represent a distance of closest approach than a time average distance.(ABSTRACT TRUNCATED AT 250 WORDS)     

The role of proline and glycine in determining the backbone flexibility of a channel-forming peptide

Alamethicin is a helical 20-amino acid voltage-gated channel-forming peptide, which is known to exhibit segmental flexibility in solution along its backbone near alpha-methylalanine (MeA)-10 and Gly-11. In an alpha-helical configuration, MeA at position 10 would normally hydrogen-bond with position 14, but the presence of proline at this position prevents the formation of this interhelical hydrogen bond. To determine whether the presence of proline at position 14 contributes to the flexibility of this helix, two analogs of alamethicin were synthesized, one with proline 14 replaced by alanine and another with both proline 14 and glycine 11 replaced by alanine. The C-termini of these peptides were derivatized with a proxyl nitroxide, and paramagnetic enhancements produced by the nitroxide on the Calpha protons were used to estimate r-6 weighted distances between the nitroxide and the backbone protons. When compared to native alamethicin, the analog lacking proline 14 exhibited similar C-terminal to Calpha proton distances, indicating that substitution of proline alone does not alter the flexibility of this helix; however, the subsequent removal of glycine 11 resulted in a significant increase in the averaged distances between the C- and N-termini. Thus, the G-X-X-P motif found in alamethicin appears to be largely responsible for mediating high-amplitude bending motions that have been observed in the central helical domain of alamethicin in methanol. To determine whether these substitutions alter the channel behavior of alamethicin, the macroscopic and single-channel currents produced by these analogs were compared. Although the substitution of the G-X-X-P motif produces channels with altered characteristics, this motif is not essential to achieve voltage-dependent gating or alamethicin-like behavior.     

TOAC spin labels in the backbone of alamethicin: EPR studies in lipid membranes

Alamethicin is a 19-amino-acid residue hydrophobic peptide that produces voltage-dependent ion channels in membranes. Analogues of the Glu(OMe)(7,18,19) variant of alamethicin F50/5 that are rigidly spin-labeled in the peptide backbone have been synthesized by replacing residue 1, 8, or 16 with 2,2,6,6-tetramethyl-piperidine-1-oxyl-4-amino-4-carboxyl (TOAC), a helicogenic nitroxyl amino acid. Conventional electron paramagnetic resonance spectra are used to determine the insertion and orientation of the TOAC(n) alamethicins in fluid lipid bilayer membranes of dimyristoyl phosphatidylcholine. Isotropic (14)N-hyperfine couplings indicate that TOAC(8) and TOAC(16) are situated in the hydrophobic core of the membrane, whereas the TOAC(1) label resides closer to the membrane surface. Anisotropic hyperfine splittings show that alamethicin is highly ordered in the fluid membranes. Experiments with aligned membranes demonstrate that the principal diffusion axis lies close to the membrane normal, corresponding to a transmembrane orientation. Combination of data from the three spin-labeled positions yields both the dynamic order parameter of the peptide backbone and the intramolecular orientations of the TOAC groups. The latter are compared with x-ray diffraction results from alamethicin crystals. Saturation transfer electron paramagnetic resonance, which is sensitive to microsecond rotational motion, reveals that overall rotation of alamethicin is fast in fluid membranes, with effective correlation times <30 ns. Thus, alamethicin does not form large stable aggregates in fluid membranes, and ionic conductance must arise from transient or voltage-induced associations.     

Backbone dynamics of the channel-forming antibiotic zervamicin IIB studied by 15N NMR relaxation

The backbone dynamics of the channel-forming peptide antibiotic zervamicin IIB (Zrv-IIB) in methanol were studied by 15N nuclear magnetic resonance relaxation measurements at 11.7, 14.1 and 18.8 T magnetic fields. The anisotropic overall rotation of the peptide was characterized based on 15N relaxation data and by hydrodynamic calculations. 'Model-free' analysis of the relaxation data showed that the peptide is fairly rigid on a sub-nanosecond time-scale. The residues from the polar side of Zrv-IIB helix are involved in micro-millisecond time-scale conformational exchange. The conformational exchange observed might indicate intramolecular processes or specific intermolecular interactions of potential relevance to Zrv-IIB ion channel formation.     

Conformational flexibility of angiotensin II. A carbon-13 spin-lattice relaxation study.

Carbon-13 spin-lattice relaxation times (T1) have been determined for the carbon in the octapeptide hormone [5-isoleucine]-angiotensin II in aqueous solution. Two possible models for molecular motion are considered: isotropic overall motion of the hormone with internal motion of some residues and anisotropic overall molecular motion. The data are interpreted in detail using the former model. The alpha carbons of the peptide backbone are all equally restricted in their motion. The correlation time for overall molecular reorientation, calculated from an everage T1 value of 95 msec for the alpha carbons in the peptide backbone, is ca. 5 times 10-10 sec. The carbons in the side chains are more mobile than those in the peptide backbone, with the exception of the side chain of the Tyr residue which does not undergo rapid segmental motion. We propose that [5-isoleucine]-angiotensin II has a restricted backbone conformation and that the alpha carbons of the N- and C-terminal residues are constrained to nearly the same extent as the remaining alpha carbons in the peptide backbone. Chemical shift data indicate that the Pro residue adopts the trans conformation about the His-Pro bond and that the imidazole ring of His has a strong preference for the N-tau -H tautomer.     

Helix bending in alamethicin: molecular dynamics simulations and amide hydrogen exchange in methanol

Molecular dynamics simulations of alamethicin in methanol were carried out with either a regular alpha-helical conformation or the x-ray crystal structure as starting structures. The structures rapidly converged to a well-defined hydrogen-bonding pattern with mixed alpha-helical and 3(10)-helical hydrogen bonds, consistent with NMR structural characterization, and did not unfold throughout the 1-ns simulation, despite some sizable backbone fluctuations involving reversible breaking of helical hydrogen bonds. Bending of the helical structure around residues Aib10-Aib13 was associated with reversible flips of the peptide bonds involving G11 (Aib10-G11 or G11-L12 peptide bonds), yielding discrete structural states in which the Aib10 carbonyl or (rarely) the G11 carbonyl was oriented away from the peptide helix. These peptide bond reversals could be accommodated without greatly perturbing the adjacent helical structure, and intramolecular hydrogen bonding was generally maintained in bent states through the formation of new (non-alpha or 3[10]) hydrogen bonds with good geometries: G11 NH-V9 CO (inverse gamma turn), Aib13 NH-Aib8 CO (pi-helix) and, rarely, L12 NH- Q7 NH (pi-helix). These observations may reconcile potentially conflicting NMR structural information for alamethicin in methanol, in which evidence for conformational flexibility in the peptide sequence before P14 (G11-Aib13) contrasts with the stability of backbone amide NH groups to exchange with solvent. Similar reversible reorientation of the Thr11-Gly12 peptide bond of melittin is also observed in dynamics simulations in methanol (R. B. Sessions, N. Gibbs, and C. E. Dempsey, submitted). This phenomenon may have some role in the orientation of the peptide carbonyl in solvating the channel lumen in membrane ion channel states of these peptides.     

Backbone dynamics of alamethicin bound to lipid membranes: spin-echo electron paramagnetic resonance of TOAC-spin labels

Alamethicin F50/5 is a hydrophobic peptide that is devoid of charged residues and that induces voltage-dependent ion channels in lipid membranes. The peptide backbone is likely to be involved in the ion conduction pathway. Electron spin-echo spectroscopy of alamethicin F50/5 analogs in which a selected Aib residue (at position n = 1, 8, or 16) is replaced by the TOAC amino-acid spin label was used to study torsional dynamics of the peptide backbone in association with phosphatidylcholine bilayer membranes. Rapid librational motions of limited angular amplitude were observed at each of the three TOAC sites by recording echo-detected spectra as a function of echo delay time, 2τ. Simulation of the time-resolved spectra, combined with conventional EPR measurements of the librational amplitude, shows that torsional fluctuations of the peptide backbone take place on the subnanosecond to nanosecond timescale, with little temperature dependence. Associated fluctuations in polar fields from the peptide could facilitate ion permeation.     

Functional analysis of myosin mutations that cause familial hypertrophic cardiomyopathy.

Two approaches employing nuclear magnetic resonance (NMR) were used to investigate the transmembrane migration rate of the C-terminal end of native alamethicin and a more hydrophobic analog called L1. Native alamethicin exhibits a very slow transmembrane migration rate when bound to phosphatidylcholine vesicles, which is no greater than 1 x 10(-4) min(-1). This rate is much slower than expected, based on the hydrophobic partition energies of the amino acid side chains and the backbone of the exposed C-terminal end of alamethicin. The alamethicin analog L1 exhibits crossing rates that are at least 1000 times faster than that of native alamethicin. A comparison of the equilibrium positions of these two peptides shows that L1 sits approximately 3-4 A deeper in the membrane than does native alamethicin (Barranger-Mathys and Cafiso. 1996. Biochemistry. 35:489). The slow rate of alamethicin crossing can be explained if the peptide helix is irregular at its C-terminus and hydrogen bonded to solvent or lipid. We postulate that L1 does not experience as large a barrier to transport because its C-terminus is already buried within the membrane interface. This difference is most easily explained by conformational differences between L1 and alamethicin rather than differences in hydrophobicity. The results obtained here demonstrate that side-chain hydrophobicity alone cannot account for the energy barriers to peptide and protein transport across membranes.     

Surface binding of alamethicin stabilizes its helical structure: molecular dynamics simulations.

Alamethicin is an amphipathic alpha-helical peptide that forms ion channels. An early event in channel formation is believed to be the binding of alamethicin to the surface of a lipid bilayer. Molecular dynamics simulations are used to compare the structural and dynamic properties of alamethicin in water and alamethicin bound to the surface of a phosphatidylcholine bilayer. The bilayer surface simulation corresponded to a loosely bound alamethicin molecule that interacted with lipid headgroups but did not penetrate the hydrophobic core of the bilayer. Both simulations started with the peptide molecule in an alpha-helical conformation and lasted 2 ns. In water, the helix started to unfold after approximately 300 ps and by the end of the simulation only the N-terminal region of the peptide remained alpha-helical and the molecule had collapsed into a more compact form. At the surface of the bilayer, loss of helicity was restricted to the C-terminal third of the molecule and the rod-shaped structure of the peptide was retained. In the surface simulation about 10% of the peptide/water H-bonds were replaced by peptide/lipid H-bonds. These simulations suggest that some degree of stabilization of an amphipathic alpha-helix occurs at a bilayer surface even without interactions between hydrophobic side chains and the acyl chain core of the bilayer.     

Calculations suggest a pathway for the transverse diffusion of a hydrophobic peptide across a lipid bilayer

Alamethicin is a hydrophobic antibiotic peptide 20 amino acids in length. It is predominantly helical and partitions into lipid bilayers mostly in transmembrane orientations. The rate of the peptide transverse diffusion (flip-flop) in palmitoyl-oleyl-phosphatidylcholine vesicles has been measured recently and the results suggest that it involves an energy barrier, presumably due to the free energy of transfer of the peptide termini across the bilayer. We used continuum-solvent model calculations, the known x-ray crystal structure of alamethicin and a simplified representation of the lipid bilayer as a slab of low dielectric constant to calculate the flip-flop rate. We assumed that the lipids adjust rapidly to each configuration of alamethicin in the bilayer because their motions are significantly faster than the average peptide flip-flop time. Thus, we considered the process as a sequence of discrete peptide-membrane configurations, representing critical steps in the diffusion, and estimated the transmembrane flip-flop rate from the calculated free energy of the system in each configuration. Our calculations indicate that the simplest possible pathway, i.e., the rotation of the helix around the bilayer midplane, involving the simultaneous burial of the two termini in the membrane, is energetically unfavorable. The most plausible alternative is a two-step process, comprised of a rotation of alamethicin around its C-terminus residue from the initial transmembrane orientation to a surface orientation, followed by a rotation around the N-terminus residue from the surface to the final reversed transmembrane orientation. This process involves the burial of one terminus at a time and is much more likely than the rotation of the helix around the bilayer midplane. Our calculations give flip-flop rates of approximately 10(-7)/s for this pathway, in accord with the measured value of 1.7 x 10(-6)/s.     

Interactions of the channel forming peptide alamethicin with artificial and natural membranes

Alamethicin and related α-aminoisobutyric acid peptides form transmembrane channels across lipid bilayers. This article briefly reviews studies on the effect of alamethicin on lipid phase transitions in lipid bilayers and on mitochondrial oxidative phosphorylation. Fluorescence polarization studies, employing 1,6-diphenyl-1,3,5-hexatriene as a probe, suggest that alamethicin fluidizes lipid bilayers below the phase transition t-emperature, but has little effect above the gel-liquid crystal transition point. Alamethicin is shown to function as an uncoupler of oxidative phosphorylation in rat liver mitochondria. The influence of alamethicin on mitochondrial respiration is modulated by the phosphate ion concentration in the medium. Classical uncoupling activity is evident at low phosphate levels while inhibitory effects set in at higher phosphate concentrations. Time-dependent changes in respiration rates following peptide addition are rationalized in terms of alamethicin interactions with mitochondrial membrane components.     

Hydrogen bond stabilities in membrane-reconstituted alamethicin from amide-resolved hydrogen-exchange measurements.

Amide-resolved hydrogen-deuterium exchange-rate constants were measured for backbone amides of alamethicin reconstituted in dioleoylphosphatidylcholine vesicles by an exchange-trapping method combined with high-resolution nuclear magnetic resonance spectroscopy. In vesicles containing alamethicin at molar ratios between 1:20 and 1:100 relative to lipid, the exchange-rate constants increased with increasing volume of the D20 buffer in which the vesicles were suspended, indicating that exchange under these conditions is dominated by partitioning of the peptide into the aqueous phase. This was supported by observation of a linear relationship between the exchange-rate constants for amides in membrane-reconstituted alamethicin and those for amides in alamethicin dissolved directly into D2O buffer. Significant protection of amides from exchange with D2O buffer in membrane-reconstituted alamethicin is interpreted in terms of stabilization by helical hydrogen bonding. Under conditions in which amide exchange occurred by partitioning of the peptide into solution, only lower limits for hydrogen-bond stabilities in the membrane were determined; all the potentially hydrogen-bonded amides of alamethicin are at least 1000-fold exchange protected in the membrane-bound state. When partitioning of alamethicin into the aqueous phase was suppressed by hydration of reconstituted vesicles in a limiting volume of water [D2O:dioleoylphosphatidylcholine:alamethicin; 220:1:0.05; (M:M:M)], the exchange-protection factors exhibited helical periodicity with highly exchange-protected, and less well-protected, amides on the nonpolar and polar helix faces, respectively. The exchange data indicate that, under the conditions studied, alamethicin adopts a stable helical structure in DOPC bilayers in which all the potentially hydrogen-bonded amides are stabilized by helical hydrogen bonds. The protection factors define the orientation of the peptide helix with respect to an aqueous phase, which is either the bulk solution or water within parallel or antiparallel transmembrane arrays of reconstituted alamethicin.     

Repeat motions and backbone flexibility in designed proteins with different numbers of identical consensus tetratricopeptide repeats

The tetratricopeptide repeat (TPR) is a 34-residue helix-turn-helix motif that occurs as three or more tandem repeats in a wide variety of proteins. We have determined the repeat motions and backbone fluctuations of proteins containing two or three consensus TPR repeats (CTPR2 and CPTR3, respectively) using 15N NMR relaxation measurements. Rotational diffusion tensors calculated from these data for each repeat within each TPR protein indicate that there is a high degree of motional correlation between different repeats in the same protein. This is consistent with the prevailing view that repeat proteins, such as CTPR2 and CTPR3, behave as single cooperatively folded domains. The internal motions of backbone NH groups were determined using the Lipari-Szabo model-free formalism. For most residues, there was a clear separation between the influence of internal motion and the influence of global rotational tumbling on the observed magnetic relaxation. The local internal motions are highly restricted in most of the helical elements, with slightly greater flexibility in the linker elements. Comparisons between CTPR2 and CTPR3 indicate that an addition of a TPR repeat to the C-terminus (before the solvation helix) of CTPR2 slightly reduces the flexibility of the preceding helix.     

Collisions between helical peptides in membranes monitored using electron paramagnetic resonance: evidence that alamethicin is monomeric in the absence of a membrane potential.

Alamethicin is a 20-amino-acid peptide that produces a voltage-dependent conductance in membranes. We investigated the state of aggregation of alamethicin in egg phosphatidylcholine and dioleoylphosphatidylcholine membranes by examining the EPR spectra obtained from an active analog of this peptide that is spin-labeled at its C-terminus. The dependence of both the linewidth and signal intensity as a function of peptide concentration exhibit exchange broadening as the peptide concentration is increased; however, the exchange rates are linear with peptide concentration as is expected for the simple diffusion of monomers. In addition, the spin-exchange rates obtained from the linebroadening are consistent with collisional rates that are predicted from free Brownian diffusion. The results provide strong evidence that in the absence of a membrane potential, alamethicin is largely monomeric in these membranes.     

Dynamics of methyl groups in proteins as studied by proton-detected 13C NMR spectroscopy. Application to the leucine residues of staphylococcal nuclease.

This paper describes the application of recently developed nuclear magnetic resonance (NMR) pulse sequences to obtain information about the internal dynamics of isotopically enriched hydrophobic side chains in proteins. The two-dimensional spectra provided by the pulse sequences enable one to make accurate measurements of nuclear Overhauser effects (NOE) and longitudinal (T1) and transverse (T2) relaxation times of enriched methyl carbons in proteins. Herein, these techniques are used to investigate the internal dynamics of the 11 leucine side chains of staphylococcal nuclease (SNase), a small enzyme having Mr = 16.8K, in the absence and presence of ligands thymidine 3',5'-bisphosphate (pdTp) and Ca2+. We report the synthesis of [5,5'-13C2]leucine, the preparation of SNase containing the labeled leucine, the sequential assignment of the leucine methyl carbons and protons in the liganded and unliganded proteins, and the measurement of the 13C T1, T2, and NOE values for the SNase leucine methyl carbons. Analysis of the relaxation parameters using the formalism of Lipari and Szabo shows that the internal motions of the leucine methyl carbons are characterized by effective correlation times tau f (5-80 ps) and tau s (less than 2 ns). The fast motion is identified with the rapid rotation of the methyl group about the C gamma-C delta bond axis, while the slow motion is associated with reorientation of the C gamma-C delta bond axis itself. The mean squared order parameters associated with the latter motion, Ss2, lie in the range 0.34-0.92. The values of Ss2 correlate reasonably well with the temperature factors of the leucine methyl carbons obtained from the crystal structures, but some are smaller than anticipated on the basis of the fact that nearly all leucine methyl carbons are buried and have temperature factors no larger than that of the leucine backbone atoms. Five leucine residues in liganded SNase and eight in unliganded SNase have values of Ss2 less than 0.71. These order parameters correspond to large amplitude motions (angular excursions of 27-67 degrees) of the C gamma-C delta bond axis. These results indicate that, in solution, the internal motions of the leucine side chains of SNase are significantly larger than suggested by the X-ray structures or by qualitative analysis of NOESY spectra. Comparison of Ss2 values obtained from liganded and unliganded SNase reveals a strong correlation between delta Ss2 and distance between the leucine methyl carbon and the ligands.(ABSTRACT TRUNCATED AT 400 WORDS)     

13C Nuclear magnetic relaxation study of segmental dynamics of hyaluronan in aqueous solutions

(13)C spin-lattice relaxation times (T(1)) and nuclear Overhauser enhancements (NOE) were measured as a function of temperature and magnetic field strength for the hetero-polysaccharide hyaluronan in water solutions. The relaxation data of the endocyclic ring carbons were successfully interpreted in terms of chain segmental motions by using the bimodal time-correlation function of Dejean de la Batie, Laupretre and Monnerie. On the basis of the calculated correlation times for segmental motion and amplitudes of librational motions of the C-H vectors at the various carbon sites of the HA repeating unit, we concluded that intramolecular hydrogen bonding of the secondary structure of HA plays a major role in the conformational flexibility of this carbohydrate molecule. The internal rotation of the free hydroxymethyl groups about the exocyclic C-5-C-6 bonds superimposed on segmental motion has been described as a diffusion process of restricted amplitude. The rate and amplitude of the internal rotation indicate that the hydroxymethyl groups are not involved in intramolecular hydrogen bonding. Finally, the motional parameters describing the local dynamics of the HA chain were correlated with the secondary structure of HA in aqueous solutions.     

Thermal motions of surface alpha-helices in the D-galactose chemosensory receptor. Detection by disulfide trapping.

The D-galactose chemosensory receptor of Escherichia coli is a .32 kDa globular protein possessing two distinct structural domains, each organized in an alpha/beta folding motif. Helices I and X lie at adjacent approximately parallel positions on the surface of the N-terminal domain, near the hinge region. In order to analyze the relative thermal motions of these two helices, the present study utilizes a generalizable disulfide trapping approach: first, site-directed mutagenesis is used to place a pair of cysteine residues at locations of interest on the protein surface, then disulfide bond formation is used to trap intramolecular cysteine-cysteine collisions resulting from thermal motions. Specifically, four engineered di-cysteine receptors have been constructed, each possessing one cysteine at position 26 on helix I, and a second cysteine at varying positions on helix X. A fifth control receptor possesses one cysteine at position 26, and a second on the opposite surface of the molecule. These surface cysteine substitutions have little or no effect on the measurable receptor parameters as judged by ligand binding equilibria and kinetics, protein stability, and 19F nuclear magnetic resonance, indicating that the engineered receptors are useful probes of native backbone dynamics. Spatial and kinetic features of backbone motions have been investigated by measuring intramolecular disulfide formation rates for cysteine pairs in the fully liganded receptor. The resulting rates decrease monotonically with increasing distance between cysteines in the crystal structure, while no disulfide formation is observed for the control pair unless the molecule is unfolded. The minimum translational amplitudes of the observed backbone motions range from 4.5 to 15.2 A, and the minimum rotational amplitudes are as large as 35 degrees. For each motion the rate of intramolecular sulfhydryl-sulfhydryl collision has been estimated from the measured rate of disulfide formation: the 4.5 and 15.2 A translations yield approximately 10(4) and approximately 10 collisions s-1 molecule-1, respectively. These collision rates, which are faster than ligand dissociation, likely underestimate the actual motional frequencies since only an undetermined fraction of the total motions yield collisions. The simplest plausible trajectory capable of producing such collisions is a rate-limiting translation of one or both helices along their long axes, coupled with minor helix rotations. When sugar is removed from the receptor, a substantial increase in backbone dynamics is observed, indicating the presence of new long-range backbone trajectories. Overall, the results suggest that internal motions in proteins may have larger amplitudes than previously observed.     

Voltage-dependent insertion of alamethicin at phospholipid/water and octane/water interfaces

Understanding the binding and insertion of peptides in lipid bilayers is a prerequisite for understanding phenomena such as antimicrobial activity and membrane-protein folding. We describe molecular dynamics simulations of the antimicrobial peptide alamethicin in lipid/water and octane/water environments, taking into account an external electric field to mimic the membrane potential. At cis-positive potentials, alamethicin does not insert into a phospholipid bilayer in 10 ns of simulation, due to the slow dynamics of the peptide and lipids. However, in octane N-terminal insertion occurs at field strengths from 0.33 V/nm and higher, in simulations of up to 100 ns duration. Insertion of alamethicin occurs in two steps, corresponding to desolvation of the Gln7 side chain, and the backbone of Aib10 and Gly11. The proline induced helix kink angle does not change significantly during insertion. Polyalanine and alamethicin form stable helices both when inserted in octane and at the water/octane interface, where they partition in the same location. In water, both polyalanine and alamethicin partially unfold in multiple simulations. We present a detailed analysis of the insertion of alamethicin into the octane slab and the influence of the external field on the peptide structure. Our findings give new insight into the mechanism of channel formation by alamethicin and the structure and dynamics of membrane-associated helices.     

13C multiplet nuclear magnetic resonance relaxation-derived ring puckering and backbone dynamics in proline-containing glycine-based peptides.

13CH2-multiplet nuclear magnetic resonance relaxation studies on proline (P)-containing glycine (G)-based peptides, GP, PG, GPG, PGG, and GPGG, provided numerous dipolar auto- and cross-correlation times for various motional model analyses of backbone and proline-ring bond rotations. Molecular dynamics simulations and bond rotation energy profiles were calculated to assess which motions could contribute most to observed relaxation phenomena. Results indicate that proline restricts backbone psi 1, psi 2, and phi 2 motions by 50% relative to those found for a polyglycine control peptide. psi 1 rotations are more restricted in the trans-proline isomer state than in the cis form. A two-state jump model best approximates proline ring puckering which in water could occur either by the C gamma endo-exo or by the C2 interconversion mechanism. The temperature dependence (5 degrees to 75 degrees C) of C beta, and C gamma, and C delta angular changes is rather flat, suggesting a near zero enthalpic contribution to the ring puckering process. In lower dielectric solvents, dimethylsulfoxide and methanol, which may mimic the hydrophobic environment within a protein, the endo-exo mechanism is preferred.