Sampling of protein dynamics in nanosecond time scale by 15N NMR relaxation and self-diffusion measurements

This paper presents a procedure for detection of intermediate nanosecond internal dynamics in globular proteins. The procedure uses 1H-15N relaxation measurements at several spectrometer frequencies and hydrodynamic calculations based on experimental self-diffusion coefficients. New heteronuclear experiments, using pulse field gradients, are introduced for the measurement of translation diffusion coefficients of 15N labeled proteins. An advanced interpretation of recently published (Luginbühl et al., Biochemistry, 36, 7305-7312 (1997)) backbone amide 15N relaxation data, measured at two spectrometers (400 and 750 MHz for 1H) for N-terminal DNA-binding domain (1-63) of 434 repressor, is presented. Non-applicability of commonly used fast (picosecond) dynamics model (FD) was justified by (i) poor fit of relaxation data by the FD model-free spectral density function both for isotropic and anisotropic models of the overall molecular tumbling; (ii) specific dependence of the overall rotation correlation times calculated from T1/T2 ratio on the spectrometer frequency; (iii) mismatch of the ratio of longitudinal 15N relaxation times T1, measured at different spectrometer frequencies, in comparison with that anticipated for the FD model; (iv) significantly underestimated overall rotation correlation time provided by the FD model (5.50+/-0.15 and 5.80+/-0.15 ns for 750 and 400 MHz spectrometer frequency respectively) in comparison with correlation time obtained from hydrodynamics. On the other hand, all relaxation and hydrodynamics data are in good correspondence with the model of intermediate (nanoseconds) dynamics. Overall rotation correlation time of 7.5+/-0.7 ns was calculated from experimental translation self-diffusion rate using hydrodynamics formalism (Garcia de la Torre, J. and Bloomfield, V.A. Quart. Rev. Biophys., 14, 81-139 (1981)). The statistical analysis of 15N relaxation data along with the hydrodynamic consideration clearly revealed that most of the residues in 434(1-63) repressor are involved in the nanosecond internal dynamics characterized by the the mean order parameters of 0.59+/-0.06 and the correlation times of ca. 5 ns.

Structure and stability of recombinant protein depend on the extra N-terminal methionine residue: S6 permutein from direct and fusion expression systems

Two permuted variants of S6 ribosomal protein were obtained in direct and fusion expression systems, respectively. The product of direct expression contained the extra N-terminal methionine residue. The structural properties and conformational stability of these permuteins were compared using 1-D (1)H-NMR, circular dichroism, intrinsic fluorescence, differential scanning calorimetry and resistance to urea-induced unfolding. A pronounced difference in all the parameters studied has been demonstrated. This means that the structure of recombinant protein can be sensitive to peculiarities of the expression and purification procedures, leading particularly to the presence or absence of the Met at the first position in the target protein sequence.

Two forms of cytotoxin II (cardiotoxin) from Naja naja oxiana in aqueous solution: spatial structures with tightly bound water molecules

1H-NMR spectroscopy data, such as NOE intraprotein and (bound water)/protein contacts, 3J coupling constants and deuterium exchange rates were used to determine the in-solution spatial structure of cytotoxin II from Naja naja oxiana snake venom (CTII). Exploiting information from two 1H-NMR spectral components, shown to be due to cis/trans isomerization of the Val7-Pro8 peptide bond, spatial structures of CTII minor and major forms (1 : 6) were calculated using the torsion angle dynamics algorithm of the DYANA program and then energy refined using the FANTOM program. Each form, major and minor, is represented by 20 resulting conformers, demonstrating mean backbone rmsd values of 0.51 and 0.71 A, respectively. Two forms of CTII preserve the structural skeleton as three large loops, including two beta-sheets with bend regions, and demonstrate structural differences at loop I, where cis/trans isomerization occurs. The CTII side-chain distribution constitutes hydrophilic and hydrophobic belts around the protein, alternating in the trend of the three main loops. Because of the Omega-shaped backbone, formed in participation with two bound water molecules, the tip of loop II bridges the tips of loops I and III. This ensures the continuity of the largest hydrophobic belt, formed with the residues of these tips. Comparison revealed pronounced differences in the spatial organization of the tips of the three main loops between CTII and previous structures of homologous cytotoxins (cardiotoxins) in solution.

A solvent model for simulations of peptides in bilayers. II. Membrane-spanning alpha-helices

We describe application of the implicit solvation model (see the first paper of this series), to Monte Carlo simulations of several peptides in bilayer- and water-mimetic environments, and in vacuum. The membrane-bound peptides chosen were transmembrane segments A and B of bacteriorhodopsin, the hydrophobic segment of surfactant lipoprotein, and magainin2. Their conformations in membrane-like media are known from the experiments. Also, molecular dynamics study of surfactant lipoprotein with different explicit solvents has been reported (Kovacs, H., A. E. Mark, J. Johansson, and W. F. van Gunsteren. 1995. J. Mol. Biol. 247:808-822). The principal goal of this work is to compare the results obtained in the framework of our solvation model with available experimental and computational data. The findings could be summarized as follows: 1) structural and energetic properties of studied molecules strongly depend on the solvent; membrane-mimetic media significantly promote formation of alpha-helices capable of traversing the bilayer, whereas a polar environment destabilizes alpha-helical conformation via reduction of solvent-exposed surface area and packing; 2) the structures calculated in a membrane-like environment agree with the experimental ones; 3) noticeable differences in conformation of surfactant lipoprotein assessed via Monte Carlo simulation with implicit solvent (this work) and molecular dynamics in explicit solvent were observed; 4) in vacuo simulations do not correctly reproduce protein-membrane interactions, and hence should be avoided in modeling membrane proteins.

A solvent model for simulations of peptides in bilayers. I. Membrane-promoting alpha-helix formation

We describe an efficient solvation model for proteins. In this model atomic solvation parameters imitating the hydrocarbon core of a membrane, water, and weak polar solvent (octanol) were developed. An optimal number of solvation parameters was chosen based on analysis of atomic hydrophobicities and fitting experimental free energies of gas-cyclohexane, gas-water, and octanol-water transfer for amino acids. The solvation energy term incorporated into the ECEPP/2 potential energy function was tested in Monte Carlo simulations of a number of small peptides with known energies of bilayer-water and octanol-water transfer. The calculated properties were shown to agree reasonably well with the experimental data. Furthermore, the solvation model was used to assess membrane-promoting alpha-helix formation. To accomplish this, all-atom models of 20-residue homopolypeptides-poly-Leu, poly-Val, poly-Ile, and poly-Gly in initial random coil conformation-were subjected to nonrestrained Monte Carlo conformational search in vacuo and with the solvation terms mimicking the water and hydrophobic parts of the bilayer. All the peptides demonstrated their largest helix-forming tendencies in a nonpolar environment, where the lowest-energy conformers of poly-Leu, Val, Ile revealed 100, 95, and 80% of alpha-helical content, respectively. Energetic and conformational properties of Gly in all environments were shown to be different from those observed for residues with hydrophobic side chains. Applications of the solvation model to simulations of peptides and proteins in the presence of membrane, along with limitations of the approach, are discussed.

NMR spatial structure of alpha-conotoxin ImI reveals a common scaffold in snail and snake toxins recognizing neuronal nicotinic acetylcholine receptors

A 600 MHz NMR study of alpha-conotoxin ImI from Conus imperialis, targeting the alpha7 neuronal nicotinic acetylcholine receptor (nAChR), is presented. ImI backbone spatial structure is well defined basing on the NOEs, spin-spin coupling constants, and amide protons hydrogen-deuterium exchange data: rmsd of the backbone atom coordinates at the 2-12 region is 0.28 A in the 20 best structures. The structure is described as a type I beta-turn (positions 2-5) followed by a distorted helix (positions 5-11). Similar structural patterns can be found in all neuronal-specific alpha-conotoxins. Highly mobile side chains of the Asp-5, Arg-7 and Trp-10 residues form a single site for ImI binding to the alpha7 receptor. When depicted with opposite directions of the polypeptide chains, the ImI helix and the tip of the central loop of long chain snake neurotoxins demonstrate a common scaffold and similar positioning of the functional side chains, both of these structural elements appearing essential for binding to the neuronal nAChRs.

Recognizing misfolded and distorted protein structures by the assumption-based similarity score

A new similarity score (sigma-score) is proposed which is able to find the correct protein structure among the very close alternatives and to distinguish between correct and deliberately misfolded structures. This score is based on the general principle 'similar likes similar', and it favors hydrophobic and hydrophilic contacts, and disfavors hydrophobic-to-hydrophilic contacts in proteins. The values of sigma-scores calculated for the high-resolution protein structures from the representative set are compared with those of alternatives: (i) very close alternatives which are only slightly distorted by conformational energy minimization in vacuo; (ii) alternatives with subsequently growing distortions, generated by molecular dynamics simulations in vacuo; (iii) structures derived by molecular dynamics simulation in solvent at 300 K; (iv) deliberately misfolded protein models. In nearly all tested cases the similarity score can successfully distinguish between experimental structure and its alternatives, even if the root mean square displacement of all heavy atoms is less than 1 A. The confidence interval of the similarity score was estimated using the high-resolution X-ray structures of domain pairs related by non-crystallographic symmetry. The similarity score can be used for the evaluation of the general quality of the protein models, choosing the correct structures among the very close alternatives, characterization of models simulating folding/unfolding, etc.

Molecular modeling of HIV-1 coreceptor CCR5 and exploring of conformational space of its extracellular domain in molecular dynamics simulation

The chemokine receptor CCR5 functions as a major fusion coreceptor for macrophage-tropic human immunodeficiency virus entry into cell. Here we report a three-dimensional model of CCR5 built using molecular modeling approach. Because the virus binds to extracellular domain of the receptor, special attention was given to conformational flexibility, hydrogen bonding, and environmental polarity properties of this protein part. Such data were obtained in the result of molecular dynamics study of the extracellular domain. It was shown that during the simulation the extracellular segments form a compact globular domain with numerous long-range hydrogen bonds between them. First loop of the receptor stays quite rigid while N-terminal region and loops 2, 3 are rather flexible. A number of amino acid residues disposed in unfavourable environment and, therefore, potentially involved in binding of CCR5 to viral glycoproteins and chemokines, was delineated. Comparison of the results with available experimental data permits a proposal that such residues in loop-1 and N-terminal part of the receptor are important for HIV-1 entry, while those in loops 2 and 3 participate in ligand binding. Perspectives of rational alteration of virus-binding activity of CCR5 are discussed.

Three-dimensional structure of binase in solution

We present the spatial structure of binase, a small extracellular ribonuclease, derived from 1H-NMR* data in aqueous solution. The total of 20 structures were obtained via torsion angle dynamics using DYANA program with experimental NOE and hydrogen bond distance constraints and phi and chi1 dihedral angle constraints. The final structures were energy minimised with ECEPP/2 potential in FANTOM program. Binase consists of three alpha-helices in N-terminal part (residues 6-16, 26-32 and 41-44), five-stranded antiparallel beta-sheet in C-terminal part (residues 51-55, 70-75, 86-90, 94-99 and 104-108) and two-stranded parallel beta-sheet (residues 22-24 and 49-51). Three loops (residues 36-39, 56-67 and 76-83), which play significant role in biological functioning of binase, are flexible in solution. The differences between binase and barnase spatial structures in solution explain the differences in thermostability of binase, barnase and their hybrids.

Spatial structure of the M3 transmembrane segment of the nicotinic acetylcholine receptor alpha subunit

The three-dimensional structure of a synthetic peptide corresponding to the putative transmembrane segment M3 (amino acid residues 277-301) of the alpha subunit of the nicotinic acetylcholine receptor from Torpedo californica has been studied by means of two-dimensional 1H-NMR spectroscopy in a chloroform/methanol (1:1) mixture containing 0.1 M LiClO4. Complete resonance assignment has been performed using double-quantum-filtered COSY (DQF-COSY), TOCSY and NOESY spectra. The spatial structure has been calculated using the Diana program on the basis of integrated intensities of NOESY spectra. HN-C(alpha)H and HC(alpha)-C(beta)H spin-spin coupling constants. Residues 279-297 of M3 form a right-handed helix (root mean square deviation is 0.032 nm for backbone atoms and 0.088 nm for all heavy atoms). The conformations of the 17 side chains have been unambiguously determined. The obtained structure is in accord with the photolabeling pattern of the membrane nicotinic acetylcholine receptor (nAChR) which suggests alpha-helical structure of M3 in the labeled portion [Blanton, M. P. & Cohen, J. B. (1994) Biochemistry 33, 2859-2872].

Effect of Coupling between Rotational and Translational Brownian Motions on NMR Spin Relaxation: Consideration Using Green Function of Rigid Body Diffusion

Using the Green function of arbitrary rigid Brownian diffusion (Goldstein, Biopolymers 33, 409-436, 1993), it was analytically shown that coupling between translation and rotation diffusion degrees of freedom does not affect the correlation functions relevant to the NMR intramolecular relaxation. It follows that spectral densities usually used for the anisotropic rotation diffusion (Woessner, J. Chem. Phys. 37, 647-654, 1962) can be regarded as exact in respect to the rotation-translation coupling for the spin system connected with a rigid body. Copyright 1998 Academic Press.

Two distinct structures of alpha-conotoxin GI in aqueous solution

The detailed analysis of conformational space of alpha-conotoxin GI in aqueous solution has been performed on the basis of two-dimensional NMR spectroscopy data using multiconformational approach. As the result, two topologically distinct interconvertible sets of GI conformations (populations of 78% and 22%) have been found. A common feature of the two sets is the Asn4-Cys7 beta-turn. The Gly8 to Tyrll region has a structure of right-handed helical turn in the major set and two sequential bends in the minor one. N-terminus and C-terminus also have different orientations, anti-parallel in the major conformational set and parallel in the minor one. An average pairwise rmsd for backbone heavy atoms is 0.56 A in the major set, 0.23 A in the minor, and 1.85 A between the structures of the two sets. The X-ray structure of GI [Guddat, L. W., Martin, J. A., Shan, L., Edmundson, A. B. & Gray, W. R. (1996) Biochemistry 35, 11329 - 11335] has the same folding pattern as the major NMR set, the average backbone rmsd between the two structures being 0.77 A.

Labeling of Torpedo californica nicotinic acetylcholine receptor subunits by cobratoxin derivatives with photoactivatable groups of different chemical nature at Lys23

Different photoactivatable derivatives of toxin 3 (CTX) Naja naja siamensis were obtained after CTX reaction with N-hydroxysuccinimide esters of p-azidobenzoic, p-azidotetraflourobenzoic, p-benzoylbenzoic and p-[3-(trifluoromethyl)-3H-diazirin-3-yl]benzoic acids. The ion-exchange HPLC profiles for the reaction products were very similar in four cases, with one predominant peak corresponding to the derivative containing the label at Lys23. After [125I]iodination, CTX photoactivatable derivatives were cross-linked to the nicotinic acetylcholine receptor from Torpedo californica under optimized conditions. The highest cross-linking yield (up to 16% of the bound toxin) was observed for azidobenzoyl-Lys23-CTX. Different receptor subunits were found to be labelled depending on the nature of the photoactivatable group: the azido derivatives labelled the gamma and delta subunits, benzoylbenzoyl derivative labelled the alpha and delta subunits, while p-[3-(trifluoromethyl)-3H-diazirin-3-yl]benzoyl derivative reacted with alpha, gamma and delta subunits. The cross-linking experiments in the presence of varying concentrations of (+)-tubocurarine demonstrated that the Lys23-attached diazirinyl group contacts the delta and alpha subunits in one ligand-binding site, whereas at the other site, for another CTX molecule, the contacts of the Lys23-diazirinyl are with gamma and alpha subunits. This means that the central loop in the two CTX molecules binds at the alpha/gamma and alpha/delta interfaces. Calculation of the sterically possible displacement of diazirinyl nitrogen, basing on the known X-ray structure of CTX, showed that this value does not exceed 13 A. The results obtained favor the disposition of the ligand-binding sites at the subunit interfaces, with the distance between alpha and delta, or alpha and gamma subunits at these sites being not more than 13 A.

Conformational independence of N- and C-domains in ribosomal protein L7/L12 and in the complex with protein L10

Isolated N- (1-37) and C-terminal (47-120) fragments of L7 protein, and pentameric (L7)4L10 complex were studied by NMR spectroscopy in solution. The results indicate that the dimer state of the 1-37 fragment with a helical hairpin conformation is identical to the N-terminal structure of the intact L7 dimer. The C-terminal domain of the L7 protein does not participate in (L7)4L10 complex formation. The overall motions of the L7 C-domains are essentially independent both in the L7 dimer and in the (L7)4L10 complex. Conformational motions on a millisecond time scale are detected in the (L7)4L10 complex. The possible relevance of these motions to the biological function of L7/L12 is discussed.

A new method to characterize hydrophobic organization of proteins: application to rational protein engineering of barnase

We present a new algorithm for characterization of protein spatial structure basing on the molecular hydrophobicity potential approach. The method is illustrated by the analysis of three-dimensional structure of barnase and barnase-barstar complex. Current approach enables identification of amino acid residues situated in unfavorable environment (these residues may be "active" for binding), and to map quantitatively hydrophobic, hydrophilic and unfavorable hydrophobic-hydrophilic intra- and inter-molecular contacts involving backbone and side-chain segments of amino acid residues. Calculation of individual contributions of amino acid residues to such contacts permits identification of structurally-important residues. The contact plots obtained with molecular hydrophobicity potential calculations, provide easy rules to choose sites for mutations, which can increase a strength of intra- or inter-molecular hydrophobic interactions. The unfavorable hydrophobic-hydrophilic contact can be mutated to favorable hydrophobic, and already existing weak hydrophobic contact can be strengthen by increasing hydrophobicity of residues in contact. Basing on the analysis of the contact plots, we suggest several mutations of barnase which are supposed to increase intramolecular hydrophobic interactions, and thus might lead to increased stability of the protein. Part of these mutations was studied previously experimentally, and indeed stabilized barnase. The other of predicted mutations were not studied experimentally yet. Several new mutations of barnase and barstar are also proposed to enhance the hydrophobic interactions on their binding interface.

Model-free approach beyond the borders of its applicability

Model calculations presented in this article show that commonly used methodology of 15N relaxation data analysis completely fails in detecting nanosecond time scale motions if the major part of the molecule is involved in these motions. New criteria are introduced for the detection of such cases, based on the dependence of the apparent overall correlation time, derived from the T1/T2 ratio, on the spectrometer frequency. Correctly estimating the overall rotation correlation time tauR was shown to play the key role in model-free data analysis. It is found, however, that in cases of slow internal motions with characteristic times of more than 3-4 ns, the effective tauR provided by the T1/T2 ratio for individual amide nitrogens can be used for the characterization of the fast picosecond internal dynamics.

Atomic solvation parameters for proteins in a membrane environment. Application to transmembrane alpha-helices

Several sets of atomic solvation parameters imitating: (i) nonpolar environment of hydrocarbon core of a membrane, (ii) aqueous solution, and (iii) weakly-polar solvents have been developed. The parameters have been incorporated into the ECEPP/2 and CHARMM force fields and employed in non-restrained Monte Carlo and molecular dynamics simulations of membrane-spanning alpha-helical peptides (segment A of bacteriorhodopsin, melittin). Through these simulations, the structure and energetics of the helices have been examined as a function of the solvation term in the potential energy function. For the peptides under study, the set (i) of atomic solvation parameters reveals good retention of the alpha-helical conformation. By contrast, the simulations in vacuum or with the parameters imitating a polar solvent (sets (ii) or (iii)) show fast helix destabilization and tight packing of the structure accompanied by significant decreasing of the surface area accessible to solvent. Increased helical propensity for amino acid residues, population of side-chain rotamers as well as hydrogen-bonding pattern in nonpolar membrane-like environment agree well with available experimental and computational data. The problems related to further applications of the membrane-mimicking sets of atomic solvation parameters to simulations of membrane proteins and peptides are addressed.

Three-dimensional structure of toxin OSK1 from Orthochirus scrobiculosus scorpion venom

A 600 MHz 1H NMR study of toxin OSK1, blocker of small-conductance Ca2+-activated K+ channels, is presented. The unambiguous sequential assignment of all the protons of the toxin was obtained using TOCSY, DQF-COSY, and NOESY experiments at pH 3.0 (10, 30, and 45 degrees C) in aqueous solution. 3J(N alpha), 3J(alphabeta) vicinal spin coupling constants were determined in high-resolution spectra. The cross-peak volumes in NOESY spectra and the coupling constants were used to define the local structure of the protein by the program HABAS and to generate torsion angle and interproton distance constraints for the program DIANA. Hydrogen-deuterium exchange rates of amide protons showed possible locations of hydrogen bonds. The hydrogen bond acceptors and disulfide bridges between residues 8-28, 14-33, and 18-35 were determined when analyzing distance distribution in preliminary DIANA structures. All constraints were used to obtain a set of 30 structures by DIANA. The resulting rms deviations over 30 structures are 1.30 A for the heavy atoms and 0.42 A for the backbone heavy atoms. The structures were refined by constrained energy minimization using the SYBYL program. Their analysis indicated the existence of the alpha-helix (residues 10-21) slightly distorted at the Cys14 residue, two main strands of the antiparallel beta-sheet (24-29, 32-38), and the extended fragment (2-6). The motif is stabilized by the disulfide bridges in the way, common to all known scorpion toxins. Using the fine spatial toxin structure, alignment of the homologues, mutagenesis analysis, and comparison of scorpion toxin family functions, we delineate some differences significant for the toxin specificity.

1H-15N backbone resonance assignments of bacteriorhodopsin

Des-(232-248)-bacteriorhodopsin was solubilized in a membrane mimicking environment of methanol-chloroform (1:1) containing 0.1 M 2HCO2N2H4 and 1H-15N backbone resonance assignment was obtained using 2D HMQC, 3D NOESY-HMQC, 3D TOCSY-HMQC and 3D HMQC-NOESY-HMQC NMR experiments. 87 cross-peaks out of 117 present in the HMQC spectrum were assigned to particular residues in 1-73 and 195-231 parts of the protein. For these residues also signals of C alpha H and C beta H protons were assigned.

Conformation of surface exposed N-terminus part of bacteriorhodopsin studied by transferred NOE technique

Interaction of the monoclonal antibody A5 raised against native bacteriorhodopsin (BR) with the synthetic peptide pGlu1-Ala-Gln-Ile-Thr-Gly-Arg7-NH2, corresponding to the amino acid sequence 1-7 was studied by transferred nuclear Overhauser effect (TRNOE) spectroscopy. The denaturing reagents and the specially designed pulse sequences which eliminate broad signals from the TRNOE spectra were used to favour evaluation of the TRNOE peaks. On the basis of the data obtained, the conformation of peptide bound with A5 was calculated. A model of the mutual arrangement of bacteriorhodopsin N-terminus and the first transmembrane alpha-helical segment 8-32 was proposed.

Topology of the secondary structure elements of ribosomal protein L7/L12 from E. coli in solution

Topology of the secondary structure elements of ribosomal protein L7/L12 has been studied. The sequential assignment was obtained for main and side chain resonances. This allows the overall secondary structure to be described. The results of high resolution NMR studies show that dimer of the ribosomal protein L7/L12 from Escherichia coli has a parallel (head-to-head) orientation of subunits, and N-terminal domain (NTD, residues Ser1-Ser33) has no contacts with the C-terminal domain (CTD, residues Lys51-Lys120). The NMR data for CTD are in line with crystallographic structure. The flexible interdomain (hinge) region (residues Ala34-Glu50) has an unordered structure, the Pro44 forming both cis and trans peptide bonds. Due to the conformational exchange the intensities of the peaks from the NTD are low. The conformation of the NTD, which is responsible for the formation of the L7/L12 dimer, is alpha-helical hairpin. the NTD dimer forms an antiparallel four-alpha-helix bundle.

The secondary structure of phospholamban: a two-dimensional NMR study

Phospholamban (PLN) is an intrinsic membrane protein of 52 amino acids which regulates the Ca2+ pump of the sarcoplasmic reticulum of heart, slow-twitch and smooth muscle (SR): it is normally assumed to exist in the membrane as a homopentamer. A monomeric analogue of phospholamban PLN(C41F), in which Cys41 was replaced by a Phe, was synthesized and its conformation studied by 1H NMR spectroscopy in a 1:1 mixture of chloroform/methanol. Most of the resonances in the 1H NMR spectra were assigned. The work has shown that the C-terminal hydrophobic portion forms a very stable alpha-helix. The hydrophilic N-terminal part adopts an alpha-helix configuration which is much less stable except for the stretch containing the phosphorylation sites.

Amino acid residue: is it structural or functional?

A new approach is suggested for delineating the structural and functional amino acid residues in proteins with known three-dimensional structure, basing on the involvement of residues in intramolecular hydrophobic and hydrophilic interactions and additional information about the conservativity of the residues. The approach is applied to the families of homologous neurotoxins and cardiotoxins. The results obtained concerning the role of amino acid residues in both families of toxins accord well with the similarity of their fold, but different mechanisms of action. Current approach can be used for detailed characterization of protein spatial structures, as well as for rational protein engineering.

Backbone dynamics of (1-71)- and (1-36)bacterioopsin studied by two-dimensional (1)H- (15)N NMR spectroscopy

The backbone dynamics of uniformly (15)N-labelled fragments (residues 1-71 and 1-36) of bacterioopsin, solubilized in two media (methanol-chloroform (1:1), 0.1 M (2)HCO(2)NH(4), or SDS micelles) have been investigated using 2D proton-detected heteronuclear (1)H-(15)N NMR spectroscopy at two spectrometer frequencies, 600 and 400 MHz. Contributions of the conformational exchange to the transverse relaxation rates of individual nitrogens were elucidated using a set of different rates of the CPMG spin-lock pulse train and were essentially suppressed by the high-frequency CPMG spin-lock. We found that most of the backbone amide groups of (1-71)bacterioopsin in SDS micelles are involved in the conformational exchange process over a rate range of 10(3) to 10(4) s(-1). This conformational exchange is supposed to be due to an interaction between two α-helixes of (1-71)bacterioopsin, since the hydrolysis of the peptide bond in the loop region results in the disappearance of exchange line broadening. (15)N relaxation rates and (1)H-(15)N NOE values were interpreted using the model-free approach of Lipari and Szabo [Lipari, G. and Szabo, A. (1982) J. Am. Chem. Soc., 104, 4546-4559]. In addition to overall rotation of the molecule, the backbone N-H vectors of the peptides are involved in two types of internal motions: fast, on a time scale <20 ps, and intermediate, on a time scale close to 1 ns. The intermediate dynamics in the α-helical stretches was mostly attributed to bending motions. A decrease in the order parameter of intermediate motions was also observed for residues next to Pro(50), indicating an anisotropy of the overall rotational diffusion of the molecule. Distinctly mobile regions are identified by a large decrease in the order parameter of intermediate motions and correspond to the N- and C-termini, and to a loop connecting the α-helixes of (1-71)bacterioopsin. The internal dynamics of the α-helixes on the millisecond and nanosecond time scales should be taken into account in the development of a model of the functioning bacteriorhodopsin.