Determination of the base recognition positions of zinc fingers from sequence analysis.

The CC/HH zinc finger is a small independently folded DNA recognition motif found in many eukaryotic proteins, which ligates zinc through two cysteine and two histidine ligands. A database of 1340 zinc fingers from 221 proteins has been constructed and a program for analysis of aligned sequences written. This paper describes sequence analysis aimed at determining the amino acid positions that recognize the DNA bases, by comparing two types of sequence variation. Using the idea that long runs of adjacent zinc fingers have arisen from internal gene duplication, the conservation of each position of the finger within the runs was calculated. The conservation of each position of the finger between homologous proteins from different species was also noted. A correlation of the two types of conservation showed clusters of related amino acids. One cluster of three positions was found to be especially variable within long runs, but highly conserved between corresponding fingers of homologous proteins; these positions are predicted to be the base contact positions. They match the amino acid positions that contact the bases in the co-crystal structure determined by Pavletich and Pabo [Science, 240, 809-817 (1991)]. An adjacent cluster of four positions on the plot may also be associated with DNA binding. This analysis shows that the base recognition positions can be identified even in the absence of a known structure for a zinc finger. These results are applicable to zinc fingers where the structure of the complex is unknown, in particular suggesting that the individual finger--DNA interaction seen in the Zif268--DNA structure has been conserved in many zinc finger--DNA interactions.     

Hydration of DNA bases: analysis of crystallographic data.

We present a systematic analysis of water structure around nucleic acid bases. We have examined 28 crystal structures of oligonucleotides, and have studied the patterns of water around the four bases, guanine, cytosine, adenine, and thymine. The geometries of water positions were calculated up to 4.00 A from base atoms. We have found conformation-dependent differences in both the geometry and extent of hydration of the bases.     

The iminoproton NMR spectrum of yeast tRNA-Phe predicted from crystal coordinates.

The ring current effects on the base paired iminoprotons in yeast tRNA-Phe have been calculated from crystal coordinates. The results in conjunction with independently determined intrinsic positions of the iminoprotons in various base pairs enable us to predict the low field NMR spectrum of yeast tRNA-Phe. It turns out that the calculated NMR spectra are very sensitive to slight changes in structure. Moreover the crystal and solution structure are identical as far as the present methods go.     

Possible configuration of dimers of Gua-Cyt bases. Computation by molecular mechanic models and density functional theory

The energies of interactions between guanine and cytosine in various mutual positions were calculated by the methods of molecular mechanics with refined atom-atom potential functions and the quantum mechanics theory of density functional. Both methods indicate three types of mutual positions of bases in local energy minima. These types correspond to (1) nearly coplanar base positions with intermolecular hydrogen bond formation (base pairing); (2) arrangements of two bases in nearly parallel planes one above another (base stacking); and (3) nearly perpendicular positions of base planes. According to the calculations, the global energy minimum corresponds to the Watson-Crick base pair with three hydrogen bonds. A specific feature of the pair is a transition from many positions of type (2) to positions of type (1) without any energy barrier. This feature is revealed by both methods. Another special feature of this pair is a deviation, for most of mutual base positions, of the amine group atoms from the ring plane, the deviation being more pronounced for Gua. These features are important for understanding the conformational behavior of DNA fragments and the RNA structure.     

Refined crystal structure of Cd, Zn metallothionein at 2.0 A resolution

The crystal structure of Cd5,Zn2-metallothionein from rat liver has been refined at 2.0 A resolution of a R-value of 0.176 for all observed data. The five Cd positions in the asymmetric unit of the crystal create a pseudo-centrosymmetric constellation about a crystallographic 2-fold axis. Consequently, the distribution of anomalous differences is almost ideally centrosymmetric. Therefore, the previously reported metal positions and the protein model derived therefrom are incorrect. Direct methods were applied to the protein amplitudes to locate the Cd positions. The new positions were used to calculate a new electron density map based on the Cd anomalous scattering and partial structure to model the metal clusters and the protein. Phases calculated from this model predict the positions of three sites in a (NH4)2WS4 derivative. Single isomorphous replacement phases calculated with these tungsten sites confirm the positions of the Cd sites from the new direct methods calculations. The refined metallothionein structure has a root-mean-square deviation of 0.016 A from ideality of bonds and normal stereochemistry of phi, phi and chi torsion angles. The metallothionein crystal structure is in agreement with the structures for the alpha and beta domains in solution derived by nuclear magnetic resonance methods. The overall chain folds and all metal to cysteine bonds are the same in the two structure determinations. The handedness of a short helix in the alpha-domain (residues 41 to 45) is the same in both structures. The crystal structure provides information concerning the metal cluster geometry and cysteine solvent accessibility and side-chain stereochemistry. Short cysteine peptide sequences repeated in the structure adopt restricted conformations which favor the formation of amide to sulfur hydrogen bonds. The crystal packing reveals intimate association of molecules about the diagonal 2-fold axes and trapped ions of crystallization (modeled as phosphate and sodium). Variation in the chemical and structural environments of the metal sites is in accord with data for metal exchange reactions in metallothioneins.     

Solvent structure in crystals of trypsin determined by X-ray and neutron diffraction.

The solvent structure in orthorhombic crystals of bovine trypsin has been independently determined by X-ray diffraction to 1.35 A resolution and by neutron diffraction to 2.1 A resolution. A consensus model of the water molecule positions was obtained using oxygen positions identified in the electron density map determined by X-ray diffraction, which were verified by comparison to D2O-H2O difference neutron scattering density. Six of 184 water molecules in the X-ray structure, all with B-factors greater than 50 A2, were found to be spurious after comparison with neutron results. Roughly two-thirds of the water of hydration expected from thermodynamic data for proteins was localized by neutron diffraction; approximately one-half of the water of hydration was located by X-ray diffraction. Polar regions of the protein are well hydrated, and significant D2O-H2O difference density is seen for a small number of water molecules in a second shell of hydration. Hydrogen bond lengths and angles calculated from unconstrained refinement of water positions are distributed about values typically seen in small molecule structures. Solvent models found in seven other bovine trypsin and trypsinogen and rat trypsin structures determined by X-ray diffraction were compared. Internal water molecules are well conserved in all trypsin structures including anionic rat trypsin, which is 65% homologous to bovine trypsin. Of the 22 conserved waters in trypsin, 19 were also found in trypsinogen, suggesting that they are located in regions of the apoprotein that are structurally conserved in the transition to the mature protein. Seven waters were displaced upon activation of trypsinogen. Water structure at crystal contacts is not generally conserved in different crystal forms. Three groups of integral structural water molecules are highly conserved in all solvent structures, including a spline of water molecules inserted between two beta-strands, which may resemble an intermediate in the formation of beta sheets during the folding of a protein.     

X-ray diffraction studies on chromatophore membrane from photosynthetic bacteria. I. Diffraction pattern of the photoreaction unit isolated from Rhodospirillum rubrum chromatophore and some characteristics of the structure

The X-ray diffraction pattern from chromatophore membranes of a photosynthetic bacterium, Rhodospirillum rubrum, indicates that a highly organized protein assembly exists in the membrane. The X-ray scatterer was solubilized from chromatophores by a mixture of cholate and deoxycholate. The basic component was identified as the photoreaction unit, which consists of light-harvesting bacteriochlorophyll proteins and a reaction center. The radial autocorrelation function, calculated directly from the X-ray intensity dats, made it possible to deduce certain structural features of the X-ray scatterer. 1. The maximum dimension of the X-ray scatterer is estimated to be 110-130 A. 2. The arrangement of the units in the chromatophore membrane is random. 3. Protein molecules in the unit form a rigid structure, being arranged mutually in fixed positions to give a distinct X-ray diffraction pattern. 4. The most probable structure is one which has rotational symmetry.     

Intramolecular interactions in pancreatic ribonucleases.

A detailed analysis of the composition and properties of hydrophobic nuclei and microclusters in pancreatic ribonuclease A (RNase A) has been carried out. Distance calculations for all noncovalently bonded atoms revealed that the average number of nonpolar contacts between a side chain of an amino acid and its neighbors is substantially larger if it involves hydrophobic residues rather than nonhydrophobic ones. However, the difference decreased when the number of contacts per nonpolar group and/or atom were calculated. Three main nuclei and five microclusters were identified, and their quantitative parameters were calculated. These nuclei include hydrophobic residues with a substantial number of nonpolar contacts with the environment (Phe 8, Phe 120, Phe 46, Tyr 25, Tyr 97, Ile 107, Leu 35, Ile 81, Val 54, Val 108, Met 29, Met 30). Hydrophobic nuclei of RNase A differ in shape and in composition, in the number of intranuclear contacts and of associated residues, as well as in their internal mobility. All eight cysteine residues are involved in nonpolar interactions with amino acid residues of hydrophobic nuclei. Active site amino acid residues of RNase A form a noncovalent contact network comprised of themselves, as well as of many conserved residues from hydrophobic nuclei. Sequence alignment with some other members of the RNase A family of proteins shows remarkable similarity in positions and in conservation of the main nonpolar residues, comprising cores of two (out of three) hydrophobic nuclei. A correlation was shown to exist between the average density of contacts for side-chain atoms and the number of amino acids to be found in the appropriate positions in the sequences of related mammalian ribonucleases. However, there are certain amino acid positions in the third, smaller nucleus, which are highly variable within the family. Taking into account that this nucleus is composed of residues belonging to different elements of the secondary structure, it is likely that the mutual orientation of these elements can be somehow different for these proteins.     

Bivariate flow karyotyping of human chromosomes: evaluation of variation in Hoechst 33258 fluorescence, chromomycin A3 fluorescence, and relative chromosomal DNA content.

The total variation of chromosome peak positions, in bivariate distributions of Hoechst 33258 and chromomycin A3 fluorescence of 19 healthy individuals, was compared with the experimental variation, determined from 23 bivariate distributions of chromosomes prepared separately from a single cell lineage. The experimental variation in Hoechst and chromomycin fluorescence and the relative chromosomal DNA content were determined from experiments performed over several days. The additional variance contributed by time was the same as the daily variance. The accuracy by which the relative chromosomal DNA content can be calculated from bivariate peak positions was investigated. A least squares method was used to fit the distributions of relative DNA content, obtained, respectively, from mono- and bivariate flow analyses of chromosomes from the same cell lineage. In general the DNA contents match quite well, but for a few chromosomes a difference was found, statistically discernible at the 5% level. The average relative chromosomal DNA content of the chromosomes from the 19 normal individuals, calculated from bivariate peak positions, showed a linear relation with the estimates published by other investigators.     

Structure of a fluid dioleoylphosphatidylcholine bilayer determined by joint refinement of x-ray and neutron diffraction data. III. Complete structure.

We present in this paper the complete structure of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) in the L alpha phase (66% RH, 23 degrees C) obtained by the joint refinement of neutron and x-ray lamellar diffraction data. The structural details obtained have previously required a large number of neutron diffraction experiments, using numerous specifically-deuterated phospholipid isomorphs (Büldt et al., 1978. Nature (Lond.). 271:182-184). The joint-refinement approach minimizes specific deuteration by utilizing independent neutron and x-ray data sets. The method yields a quasimolecular structure consisting of a series of multiatomic fragments that are each represented by one or several Gaussian distributions whose positions and widths can be determined to within 0.06 to 0.52 A exclusive of the methylene region. The image of DOPC at 66% RH (5.36 +/- 0.08 waters per lipid) is consistent with many aspects of bilayer structure previously determined by structural and spectroscopic studies. The most striking feature of the structure is the large amount of transbilayer thermal motion suggested by the widths and overlaps of the Gaussian envelopes of the quasimolecular fragments. We discuss the "dynamic bilayer thickness" which describes the minimum effective thickness of the hydrocarbon permeability barrier in terms of the thermal motion of the water. A gradient of thermal motion exists that increases in either direction away from the glycerol backbone which is the most constrained portion of the bilayer. The steric interactions between headgroups of apposed bilayers, expected at the hydration level of our experiments, are clearly revealed. A useful consequence of the quasimolecular structure is that average boundaries within bilayers calculated using composition and volumetric data and ad hoc assumptions can be related to the positions of the principal structural groups. Several measures of "bilayer thickness" in common use can be identified as the positions of the cholines for Luzzati's d1 (Luzzati and Husson. 1962. J. Cell Biol. 12:207-219) and the glycerols for Small's dL (Small. 1967. J. Lipid Res. 8:551-556). We do not know if these relations will be true at other hydrations or for other lipids. Of particular interest is the fact that the position of the carbonyl groups marks the average hydrocarbon/headgroup boundary. It must be emphasized, however, that this region of the bilayer must be generally characterized as one of tumultuous chemical heterogeneity because of the thermal motion of the bilayer.     

Structural studies on the 30 S ribosomal subunit from Escherichia coli

Small-angle X-ray scattering curves computed for various 30 S subunit structures have been compared with the experimental scattering curve. The curve from the 30 S subunit is best approximated by that calculated for a 1:3.6:3.6 ellipsoidal structure. The rather prolate ellipsoidal model suggested by recent electron microscope studies gives a scattering curve considerably different from the 30 S curve, suggesting that the electron microscope model is not that found in solution. Analysis of the more extended portions of the experimental scattering curve suggests some internal structure. A model is proposed that contains RNA and protein in positions such that the calculated scattering curve shows more extensive, yet similar internal structure. Resultant constraints on the structure of the 30 S subunit in solution are given.     

The urea amidolyase (DUR1,2) gene of Saccharomyces cerevisiae.

The DNA sequence of the urea amidolyase (DUR1,2) gene from S. cerevisiae has been determined. The polypeptide structure deduced from the DNA sequence contains 1,835 amino acid residues and possesses a calculated weight of 201,665 daltons which favorably correlates with that predicted from compositional analysis of purified protein (1,881 amino acid residues and a molecular weight of 203,900). The C-terminal 57 residues of the polypeptide exhibit significant homology with similarly situated sequences found in five other biotin carboxylases whose primary structures have been determined or deduced from protein and DNA sequence data, respectively. Major S1 nuclease protection fragments derived from DUR1,2 RNA-DNA hybrids exhibit apparent termini at positions -140 and -141 upstream of the coding region. The termini of minor protection fragments also occur at eleven other positions as well.     

Solution structure of B. subtilis acyl carrier protein

BACKGROUND: Acyl carrier protein (ACP) is a fundamental component of fatty acid biosynthesis in which the fatty acid chain is elongated by the fatty acid synthetase system while attached to the 4'-phosphopantetheine prosthetic group (4'-PP) of ACP. Activation of ACP is mediated by holo-acyl carrier protein synthase (ACPS) when ACPS transfers the 4'-PP moiety from coenzyme A (CoA) to Ser36 of apo-ACP. Both ACP and ACPS have been identified as essential for E. coli viability and potential targets for development of antibiotics. RESULTS: The solution structure of B. subtilis ACP (9 kDa) has been determined using two-dimensional and three-dimensional heteronuclear NMR spectroscopy. A total of 22 structures were calculated by means of hybrid distance geometry-simulated annealing using a total of 1,050 experimental NMR restraints. The atomic rmsd about the mean coordinate positions for the 22 structures is 0.45 +/- 0.08 A for the backbone atoms and 0.93 +/- 0.07 A for all atoms. The overall ACP structure consists of a four alpha-helical bundle in which 4'-PP is attached to the conserved Ser36 that is located in alpha helix II. CONCLUSIONS: Structural data were collected for both the apo and holo forms of ACP that suggest that the two forms of ACP are essentially identical. Comparison of the published structures for E. coli ACP and actinorhodin polyketide synthase acyl carrier protein (act apo-ACP) from Streptomyces coelicolor A3(2) with B. subtilis ACP indicates similar secondary structure elements but an extremely large rmsd between the three ACP structures (>4.3 A). The structural difference between B. subtilis ACP and both E. coli and act apo-ACP is not attributed to an inherent difference in the proteins, but is probably a result of a limitation in the methodology available for the analysis for E. coli and act apo-ACP. Comparison of the structure of free ACP with the bound form of ACP in the ACP-ACPS complex reveals a displacement of helix II in the vicinity of Ser36. The induced perturbation of ACP by ACPS positions Ser36 proximal to coenzyme A and aligns the dipole of helix II to initiate transfer of 4'-PP to ACP.     

Aspartic protease inhibitors. An integrated approach for the design andsynthesis of diaminodiol-based peptidomimetics.

Aspartic proteases play key roles in a variety of pathologies, including acquired immunodeficiency syndrome. Peptidomimetic inhibitors can act as drugs to combat these pathologies. We have developed an integrated methodology for preparing human immunodeficiency virus (HIV)-1 aspartic protease diaminodiol inhibitors, based on a computational method that predicts the potential inhibitory activity of the designed structures in terms of calculated enzyme-inhibitor complexation energies. This is combined with a versatile synthetic strategy that couples a high degree of stereochemical control in the central diaminodiol module with complete flexibility in the choice of side chains in the core and in flanking residues. A series of 23 tetrameric, pentameric and hexameric inhibitors, with a wide range of calculated relative complexation energies (-47.2 to +117 kJ.mol-1) and predicted hydrophobicities (logPo/w = 1.8-8.4) was thus assembled from readily available amino acids and carboxylic acids. The IC50 values for these compounds ranged from 3.2 nM to 90 microM, allowing study of correlations between structure and activity, and individuation of factors other than calculated complexation energies that determine the inhibition potency. Multivariable regression analysis revealed the importance of side-chain bulkiness and rigidity at the P2, P2' positions, suggesting possible improvements for the prediction process used to select candidate structures.     

NMR structure of the Euplotes raikovi pheromone Er-23 and identification of its five disulfide bonds.

The NMR solution structure of the 51 residue pheromone Er-23 from the ciliated protozoan Euplotes raikovi (Er) was calculated with the torsion angle dynamics program DYANA from 582 nuclear Overhauser enhancement (NOE) upper limit distance constraints, 46 dihedral angle constraints and 30 disulfide bond constraints. The disulfide bridges had not been assigned by chemical methods, and initially were assigned tentatively on the basis of inspection of the positioning of the Cys sulfhydryl groups in a bundle of 20 conformers that was calculated without disulfide bond constraints. The assignment of disulfide bridges was then validated by structure calculations that assessed the compatibility of plausible alternative Cys-Cys disulfide combinations with the input of NOE upper distance constraints and dihedral angle constraints. For a group of 20 conformers used to characterize the solution structure, the average pairwise root-mean-square distances from the mean coordinates calculated for the backbone heavy atoms N, C(alpha) and C' of resideus 1-51 is 0.38 A. The molecular architecture consists of a three-dimensional arrangement of five helices comprised of residues 2-8, 14-17, 26-29, 34-36 and 38-47, with five disulfide bridges in the positions 3-24, 6-16, 13-47, 27-40, and 35-51, which has so far not been represented in the Protein Data Bank. Er-23 is unique among presently known Er-pheromones with respect to size, sequence, the number of disulfide bonds and the three-dimensional structure, thus providing a new structural basis for rationalizing the physiological functions of this protein family.     

The molecular structure of a dimer composed of the variable portions of the Bence-Jones protein REI refined at 2.0-A resolution.

The structure of the variable portions of a K-type Bence-Jones protein REI forming a dimer has been determined by X-ray diffraction to a resolution of 2.0 A. The structure has been refined using a constrained crystallographic refinement procedure. The final R value is 0.24 for 15000 significantly measured reflections; the estimated standard deviation of atomic positions is 0.09 A. A more objective assessment of the error in the atomic positions is possible by comparing the two independently refined monomers. The mean deviation of main-chain atoms of the two chains in internal segments in 0.22 A, of main-chain dihedral angles 6.3 degrees for these segments. The unrefined molecular structure of the VREI dimer has been published (Epp, O., Colman, P., Fehlhammer, H., Bode, W., Schiffer, M., Huber, R., and Palm, W. (1974), Eur. J. Biochem. 45, 513). Now a detailed analysis is presented in terms of hydrogen bonds and conformational angles. Secondary structural elements (antiparallel beta structure, reverse turns) are defined. A more precise atomic arrangement of the amino acid residues forming the contact region and the hapten binding site is given as well as the localization of solvent molecules. Two cis-prolines (Pro-8 and Pro-95) were detected. The intrachain disulfide bridge (Cys-23-Cys-88) occurs statistically in two alternative conformations. The structure suggests reasons for strong conservation of several amino acid residues. The knowledge of the refined molecular structure enables crystal structure analyses of related molecules to be made by Patterson search techniques. The calculated phases based on the refined structure are much improved compared to isomorphous phases. Therefore the effects of hapten binding on the molecular structure can be analyzed by the difference Fourier technique with more reliability. Hapten binding studies have been started.     

SuperStar: a knowledge-based approach for identifying interaction sites in proteins.

An empirical method for identifying interaction sites in proteins is described and validated. The method is based entirely on experimental information about non-bonded interactions occurring in small-molecule crystal structures. These data are used in the form of scatterplots that show the experimentally observed distribution of one functional group (the "contact group" or "probe") around another. A template molecule (e.g. a protein binding site) is broken down into structure fragments and the scatterplots, showing the distribution of a chosen probe around these structure fragments, are superimposed on the corresponding parts of the template. The scatterplots are then translated into a three-dimensional map that shows the propensity of the probe at different positions around the template molecule. The method is illustrated for l -arabinose-binding protein, complexed with l -arabinose and with d -fucose, and for dihydrofolate reductase complexed with methotrexate. The method is validated on 122 X-ray structures of protein-ligand complexes. For all the binding sites of these proteins, propensity maps are generated for four different probes: a charged NH+3nitrogen, a carbonyl oxygen, a hydroxyl oxygen and a methyl carbon atom. Next, the maps are compared with the experimentally observed positions of ligand atoms of these types. For 74% of these ligand atoms (84% of the solvent-inaccessible ones) the calculated propensity of the matching probe at the experimental positions is higher than expected by chance. For 68% of the atoms (82% of the solvent-inaccessible ones) the propensity of the matching probe is higher than that of the other three probes. These results indicate that the approach generally gives good predictions for protein-ligand interactions. The potential applications of the propensity maps range from an aid in manual docking and structure-based drug design to their use in pharmacophore development.     

The amino acid sequence of human insulin-like growth factor I and its structural homology with proinsulin

The complete amino acid sequence of human insulin-like growth factor I (IGF-I), a polypeptide isolated from serum, has been determined. IGF-I is a single chain polypeptide of 70 amino acid residues cross-linked by three disulfide bridges. The calculated molecular weight is 7649. IGF-I displays obvious homology to proinsulin: positions 1 to 29 are homologous to insulin B chain and positions 42 to 62 to insulin A chain. A shortened "connecting" peptide with 12 residues (positions 30 to 41) compared to 30 to 35 in proinsulins shows no homology to proinsulin C peptide. An octapeptide sequence at the COOH-terminal end is also a feature not found in proinsulins. The number of differences in amino acid positions between IGF-I and insulins suggests that duplication of the gene of the common ancestor of proinsulin and IGF occurred before the time of appearance of the vertebrates. Of the 19 residues known to be invariant in all insulins so far sequenced, only glutamine A5 and asparagine A21 are replaced in IGF-I by glutamic acid and alanine, respectively. The fact that all half-cystine and glycine residues and most nonpolar core residues of the insulin monomer are conserved is compatible with a three-dimensional structure of IGF-I similar to that of insulin.     

Effect of anisotropy and anharmonicity on protein crystallographic refinement. An evaluation by molecular dynamics

Molecular dynamics simulations are employed to determine the errors introduced by anharmonicity and anisotropy in the structure and temperature factors obtained for proteins by refinement of X-ray diffraction data. Simulations (25 ps and 300 ps) of metmyoglobin are used to generate time-averaged diffraction data at 1.5 A resolution. The crystallographic restrained-parameter least-squares refinement program PROLSQ is used to refine models against these simulated data. The resulting atomic positions and isotropic temperature factors are compared with the average structure and fluctuations calculated directly from the simulations. It is found that significant errors in the atomic positions and fluctuations are introduced by the refinement, and that the errors increase with the magnitude of the atomic fluctuations. Of particular interest is the fact that the refinement generally underestimates the atomic motions. Moreover, while the actual fluctuations go up to a mean-square value of about 5 A2, the X-ray results never go above approximately 2 A2. This systematic deviation in the motional parameters appears to be due to the use of a single-site isotropic model for the atomic fluctuations. Many atoms have multiple peaks in their probability distribution functions. For some atoms, the multiple peaks are seen in difference electron density maps and it is possible to include these in the refinement as disordered residues. However, for most atoms the refinement fits only one peak and neglects the rest, leading to the observed errors in position and temperature factor. The use of strict stereochemical restraints is inconsistent with the average dynamical structure; nevertheless, refinement with tight restraints results in structures that are comparable to those obtained with loose restraints and better than those obtained with no restraints. The results support the use of tight stereochemical restraints, but indicate that restraints on the variation of temperature factors are too restrictive.     

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)