The quantitation of nuclear Overhauser effect methods for total conformational analysis of peptides in solution. Application to gramicidin S.

The [1H:1H] nuclear Overhauser effects (NOE's) and spin-lattice relaxation times (T1's) are reported for the backbone protons of the decapeptide gramicidin S. Several methods for calculating interproton distances from these measurements are presented. Ratios of interproton distances were obtained from [1H:1H] NOE's and from the combination of [1H:1H]NOE'S and T1 values. Actual proton-proton distances were calculated from these ratios either by using the known distance between two geminal protons or distances derived from scalar coupling constants. The interproton distances calculated for gramicidin S are consistent with a II' beta-turn/antiparallel beta-sheet conformation.     

Structure determination of a DNA octamer in solution by NMR spectroscopy. Effect of fast local motions.

NMR structures of biomolecules are primarily based on nuclear Overhauser effects (NOEs) between protons. For the interpretation of NOEs in terms of distances, usually the assumption of a single rotational correlation time corresponding to a rigid molecule approximation is made. Here we investigate the effect of fast internal motions of the interproton vectors in the context of the relaxation matrix approach for structure determination of biomolecules. From molecular dynamics simulations generalized order parameters were calculated for the DNA octamer d(GCGTTCGC).d(CGCAACGC), and these were used in the calculation of NOE intensities. The magnitudes of the order parameters showed some variation for the different types of interproton vectors. The lowest values were observed for the interresidue base H6/H8-H2" proton vectors (S2 = 0.60), while the cytosine H5-H6 interproton vectors were among the most motionally restricted (S2 = 0.92). Inclusion of the motion of the interproton vectors resulted in a much better agreement between theoretically calculated NOE spectra and the experimental spectra measured by 2D NOE spectroscopy. The interproton distances changed only slightly, with a maximum of 10%; nevertheless, the changes were significant and resulted in constraints that were better satisfied. The structure of the DNA octamer was determined by using restrained molecular dynamics simulations with H2O as a solvent, with and without the inclusion of local internal motions. Starting from A- or B-DNA, the structures showed a high local convergence (0.86 A), while the global convergence for the octamer was ca. 2.6 A.     

A 2D NMR study of the internal flexibility of the antifungal peptide stendomycin

A 2-D ¹H NMR study (NOESY, COSY, HOHAHA and ROESY experiments) of the antifungal peptide stendomycin is presented. The variation of the NOESY cross peak intensities is measured as a function of temperature in order to discriminate between constant and fluctuating interproton distances. It is shown that among 71 NOESY cross peaks, only 12 correspond to well defined interproton distances and their correlation time is determined. The other cross peaks cannot be translated accurately in terms of distances owing to internal molecular motions. ¹H}–¹³C nOe measurements confirm the internal mobility of the molecule. Finally a flexibility map of stendomycin can be established. Offprint requests to: J. P. Simorre     

Distance evaluation via heteronuclear SQC-NOESY experiments

Summary A method for the quantitative determination of interproton distances from¹H NOE relayed heteronuclear correlation is presented. Model compounds are investigated. Accurate distances are obtained if all the factors affecting such distance measurements, such as the local mobility and the presence of strong coupling in heteronuclear systems, are properly taken into account.     

Structure, dynamics and mapping of membrane-binding residues of micelle-bound antimicrobial peptides by natural abundance C NMR spectroscopy

Worldwide bacterial resistance to traditional antibiotics has drawn much research attention to naturally occurring antimicrobial peptides (AMPs) owing to their potential as alternative antimicrobials. Structural studies of AMPs are essential for an in-depth understanding of their activity, mechanism of action, and in guiding peptide design. Two-dimensional solution proton NMR spectroscopy has been the major tool. In this article, we describe the applications of natural abundance ¹³C NMR spectroscopy that provides complementary information to 2D ¹H NMR. The correlation of ¹³Cα secondary shifts with both 3D structure and heteronuclear ¹⁵N NOE values indicates that natural abundance carbon chemical shifts are useful probes for backbone structure and dynamics of membrane peptides. Using human LL-37-derived peptides (GF-17, KR-12, and RI-10), as well as amphibian antimicrobial and anticancer peptide aurein 1.2 and its analog LLAA, as models, we show that the cross peak intensity plots of 2D ¹H-¹³Cα HSQC spectra versus residue number present a wave-like pattern (HSQC wave) where key hydrophobic residues of micelle-bound peptides are located in the troughs with weaker intensities, probably due to fast exchange between the free and bound forms. In all the cases, the identification of aromatic phenylalanines as a key membrane-binding residue is consistent with previous intermolecular Phe-lipid NOE observations. Furthermore, mutation of one of the key hydrophobic residues of KR-12 to Ala significantly reduced the antibacterial activity of the peptide mutants. These results illustrate that natural abundance heteronuclear-correlated NMR spectroscopy can be utilized to probe backbone structure and dynamics, and perhaps to map key membrane-binding residues of peptides in complex with micelles. ¹H-¹³Cα HSQC wave, along with other NMR waves such as dipolar wave and chemical shift wave, offers novel insights into peptide-membrane interactions from different angles.     

Internal mobility in a double-stranded B DNA hexamer and undecamer. A time-dependent proton-proton nuclear Overhauser enhancement study

The internal mobility of the deoxyribose H2'-H2" and base C(H5)-C(H6) and T(CH3)-T(H6) vectors has been investigated by means of time-dependent nuclear Overhauser enhancement (NOE) measurements in a B DNA hexamer and undecamer. Cross-relaxation rates between these proton pairs are determined from the initial slopes of the time development of the NOEs, and, as the interproton distances between these proton pairs are fixed, apparent correlation times for the 3 interproton vectors are calculated from the cross-relaxation rate data. It is shown that there is little residue to residue variation in the cross-relaxation rates of the interproton vectors within each oligonucleotide, that the mean apparent correlation times of the C(H5)-C(H6) and T(CH3)-T(H6) vectors are approximately equal and significantly greater than that of the H2'-H2" vectors, and that the data for the H2'-H2" vectors of both oligonucleotides and the C(H5)-C(H6) and T(CH3)-T(H6) vectors of the undecamer cannot be accounted for by isotropic tumbling alone. The data are analysed in terms of a two motion model with isotropic tumbling and a single internal motion. The relaxation time of the internal motion at 23 degrees C is less than or equal to 1 ns for the H2'-H2" vectors of both oligonucleotides and less than or equal to 3 ns for the C(H5)-C(H6) and T(CH3)-T(H6) vectors of the undecamer. In the case of the H2'-H2" vectors, however, the amplitude of the internal motion is found to be too large to be compatible with the known stereochemistry of DNA. This finding can only be explained by invoking additional degrees of internal freedom with a larger number of internal motions of small amplitude of the type deduced from the analysis of crystallographic thermal factors [(1984) J. Mol. Biol. 173, 361-388].     

Prospects for resonance assignments in multidimensional solid-state NMR spectra of uniformly labeled proteins

Summary The feasibility of assigning the backbone ¹⁵N and ¹³C NMR chemical shifts in multidimensional magic angle spinning NMR spectra of uniformly isotopically labeled proteins and peptides in unoriented solid samples is assessed by means of numerical simulations. The goal of these simulations is to examine how the upper limit on the size of a peptide for which unique assignments can be made depends on the spectral resolution, i.e., the NMR line widths. Sets of simulated three-dimensional chemical shift correlation spectra for artificial peptides of varying length are constructed from published liquid-state NMR chemical shift data for ubiquitin, a well-characterized soluble protein. Resonance assignments consistent with these spectra to within the assumed spectral resolution are found by a numerical search algorithm. The dependence of the number of consistent assignments on the assumed spectral resolution and on the length of the peptide is reported. If only three-dimensional chemical shift correlation data for backbone ¹⁵N and ¹³C nuclei are used, and no residue-specific chemical shift information, information from amino acid side-chain signals, and proton chemical shift information are available, a spectral resolution of 1 ppm or less is generally required for a unique assignment of backbone chemical shifts for a peptide of 30 amino acid residues.     

Three-dimensional structure of the ion-chanel forming peptide trichorzianin TA VII bound to sodium dodecyl sulfate micelles

Trichorzianin TA VII, Ac⁰ U¹ A² A³ U⁴ J⁵ Q⁶ U⁷ U⁸ U⁹ S¹⁰ L¹¹ U¹² P¹³ V¹⁴ U¹⁵ I¹⁶ Q¹⁷ Q¹⁸ Fol¹⁹, is a nonadecapeptide member of the peptaibol antibiotics biosynthesized by Trichoderma soil fungi, which is characterized by a high proportion of the α,α-dialkylated amino acids, α-aminoisobutyric acid (Aib, U) and isovaline (Iva, J), an acetylated N-terminus and a C-terminal phenylalaninol (Pheol, Fol). The main interest in such peptides stems from their ability to interact with phospholipid bilayers and form voltage-dependent transmembrane channels in planar lipid bilayers. In order to provide insights into the lipid-peptide interaction promoting the voltage gating, the conformational study of TA VII in the presence of perdeuterated sodium dodecyl sulfate (SDS-d25) micelles has been carried out. ¹H sequential assignments have been performed with the use of two-dimensional homo- and -heteronuclear nmr techniques including double quantum filtered correlated spectroscopy, homonuclear Hartmann-Hahn, nuclear Overhauser effect spectroscopy, ¹H-¹³C heteronuclear single quantum correlation, and heteronuclear multiple bond correlation. Conformational parameters, such as ³ J _NHCαH coupling constants, temperature coefficients of amide protons (Δδ/Δ T _NH) and quantitative nuclear Overhauser enhancement data, lead to detailed structural information. Ninety-eight three-dimensional structures consistent with the nmr data were generated from 231 interproton distances and six Φ dihedral angle restraints, using restrained molecular dynamics and energy minimization calculations. The average rms deviation between the 98 refined structures and the energy-minimized average structure is 0.59 Å for the backbone atoms. The structure of trichorzianin TA VII associated with SDS micelles, as determined by these methods, is characterized by two right-handed helical segments involving residues 1–8 and 11–19, linked by a β-turn that leads to an angle about 90°–100° between the two helix axes; residues 18 and 19 at the end of the C-terminal helix exhibit multiple conformations. © 1998 John Wiley & Sons, Inc. Biopoly 46: 75–88, 1998     

Proton Nuclear Magnetic Resonance Studies on Huwentoxin-I from the Venom of the Spider Selenocosmia huwena: 2. Three-Dimensional Structure in Solution

The three-dimensional structure in aqueous solution of native huwentoxin-I, a neurotoxin from the venom of the spider Selenocosmia huwena, has been determined from two-dimensional ¹H NMR data recorded at 500 and 600 MHz. Structural constraints consisting of interproton distances inferred from NOEs and dihedral angles from spin–spin coupling constants were used as input for distance geometry calculation with the program XPLOR 3.1. The best 10 structures have NOE violations <0.3 Å, dihedral violations <2°, and pairwise root-mean-square differences of 1.08 (±0.20) Å over backbone atoms (N, C^α;, C). The molecule adopts a compact structure consisting of a small triple-stranded antiparallel β-sheet and five β-turns. A small hydrophobic patch consisting of Phe 6, Trp 28, and Trp 31 is located on one side of the molecule. All six lysine residues are distributed on the molecular surface. The three disulfidc bridges are buried within the molecule. The structure contains an “inhibitor cystine knot motif” which is adopted by several other small proteins, such as ω-conotoxin, agatoxin IVA, and gurmarin.     

1H, 13C and 15N NMR assignments and secondary structure of the paramagnetic form of rat cytochrome b 5

Summary Modern multidimensional double- and triple-resonance NMR methods have been applied to assign the backbone and side-chain ¹³C resonances for both equilibrium conformers of the paramagnetic form of rat liver microsomal cytochrome b ₅. The assignment of backbone ¹³C resonances was used to confirm previous ¹H and ¹⁵N resonance assignments [Guiles, R.D. et al. (1993) Biochemistry, 32, 8329–8340]. On the basis of short- and medium-range NOEs and backbone ¹³C chemical shifts, the solution secondary structure of rat cytochrome b ₅ has been determined. The striking similarity of backbone ¹³C resonances for both equilibrium forms strongly suggests that the secondary structures of the two isomers are virtually identical. It has been found that the ¹³C chemical shifts of both backbone and side-chain atoms are relatively insensitive to paramagnetic effects. The reliability of such methods in anisotropic paramagnetic systems, where large pseudocontact shifts can be observed, is evaluated through calculations of the magnitude of such shifts.Abbreviations DANTE delays alternating with nutation for tailored excitation - DEAE diethylaminoethyl - DQF-COSY 2D double-quantum-filtered correlation spectroscopy - EDTA ethylenediaminetetraacetic acid - HCCH-TOCSY 3D proton-correlated carbon TOCSY experiment - HMQC 2D heteronuclear multiple-quantum correlation spectroscopy - HNCA 3D triple-resonance experiment correlating amide protons, amide nitrogens and alpha carbons - HNCO 3D triple-resonance experiment correlating amide protons, amide nitrogens and carbonyl carbons - HNCOCA 3D triple-resonance experiment correlating amide protons, amide nitrogens and alpha carbons via carbonyl carbons - HOHAHA 2D homonuclear Hartmann-Hahn spectroscopy - HOHAHA-HMQC 3D HOHAHA relayed HMQC - HSQC 2D heteronuclear single-quantum correlation spectroscopy - IPTG isopropyl thiogalactoside - NOESY 2D nuclear Overhauser enhancement spectroscopy - NOESY-HSQC 3D NOESY relayed HSQC - TOCSY 2D total correlation spectroscopy - TPPI time-proportional phase incrementation - TSP trimethyl silyl propionate     

1H-nmr studies of polyoxyethylene-bound homo-oligo-L-methionines

The use of ¹H-nmr spectroscopy is demonstrated to be a useful analytical method to characterize the structure of synthetic peptides attached to soluble, macromolecular polyoxyethylene (POE) supports in the liquid-phase method (LPM) of peptide synthesis. We report an extensive 360-MHz ¹H-nmr study of POE-bound homo-oligo-L -methionine peptides. A combination of high field and selective saturation or Redfield pulse methods allows resolution of individual backbone NH and α-CH resonances of dilute peptides in the presence of strong resonances from macromolecular POE and/or protonated solvents. The nmr spectra for the POE-bound peptides in CDCl₃ are qualitatively similar to those of the low-molecular-weight Boc-L -Met_n-OMe peptide esters. This corroborates other observations that POE has little effect on peptide stucture. The backbone α-CH region of peptides is overlapped by signals from the terminal oxyethylene group of POE, but the peptide side-chain and low-field backbone NH resonances are well resolved. In trifluoroethanol the Boc-(L -Met)_n-NH-POE heptamer and octamer adopt the right-handed α-helical structure, and the present nmr studies provide evidence for two strong intramolecular hydrogen bonds to stabilize the helices. In water, the N-deblocked derivatives, (L -Met)_n-NH-POE oligomers adopt β-sheet structure and manifest well-resolved nonequivalent NH resonances with 6–7 Hz ³J_NH-CH coupling constants.     

Molecular mobility of polypeptides containing proline as determined by 13C magnetic resonance

The molecular conformations and dynamics of poly(L -prolyl), poly(hydroxyl-L -prolyl), poly(L -prolyl-glycyl), poly(hydroxyl-L -prolyl), and poly(glycyl-glycyl-L -prolyl-glycyl), in aqueous solution, have been studied using ¹³C pulse Fourier transform nmr spectroscopy. From a measurement of the intensities of major and minor resonances in the spectra of the copolypeptides, it was determined that 15–20% of the glycyl-prolyl and glycyl-hydroxyprolyl peptide bonds are cis. Effective rotational correlation times (τ_eff), obtained from measurements of spin-lattice relaxation times (T₁) of individual backbone and side-chain carbons, demonstrated that backbone reorientation is approximately isotropic for the five polypeptides and is characterized by correlation times of ca. 0.3–0.6 nanoseconds as a result of rapid segmental motion. In a given polypeptide glycyl and pyrrolidine residues were found to have the same backbone correlation times, but backbone carbon τ_eff values did decrease as the glycyl content of the peptides increased. A semi-quantitative analysis of C^β, C^γ, and C^δ correlation times suggests that rapid ring motion in both prolyl and hydroxyprolyl involves primarily C^γ and C^β, with the prolyl ring being more mobile than the hydroxyprolyl ring.     

Internal mobility of cyclic RGD hexapeptides studied by 13C NMR relaxation and the model-free approach

Summary The internal mobility of three isomeric cyclic RGD hexapeptides designed to contain two -turns in defined positions, cyclo(Arg-Gly-Asp-Gly-d-Pro-Pro) (I), cyclo(Arg-Gly-Asp-d-Pro-Gly-Pro) (II) and cyclo(Arg-Gly-Asp-d-Pro-Pro-Gly) (III), have been studied by ¹³C NMR longitudinal and transverse relaxation experiments and measurements of steady-state heteronuclear {¹H}-¹³C NOE enhancement with ¹³C at natural abundance. The data were interpreted according to the model-free formalism of Lipari and Szabo, which is usually applied to data from macromolecules or larger sized peptides with overall rotational correlation times exceeding 1 ns, to yield information about internal motions on the 10–100 ps time scale. The applicability of the model-free analysis with acceptable uncertainties to these small peptides, with overall rotational correlation times slightly below 0.3 ns, was demonstrated for this specific instance. Chemical exchange contributions to T₂ from slower motions were also identified in the process. According to the order parameters obtained for its backbone -carbon atoms, II has the most rigid backbone conformation on the 10–100 ps time scale, and I the most flexible. This result coincides with the results of earlier NMR-constrained conformational searches, which indicated greatest uncertainty in the structure of I and least in II.     

Sequence-specific assignment of the backbone 1H-and 31P-NMR lines in a short DNA duplex with homo- and heteronuclear correlated spectroscopy

We report the assignment of the backbone ¹H- and ³¹P-nmr lines in the synthetic hexadeoxyribonucleotide pentaphosphate duplex d(GCATGC)₂, using double quantum filtered ¹H-¹H correlation spectroscopy, ¹H observed ¹H-³¹P heteronuclear correlation spectroscopy, and ³¹P relayed ¹H-¹H correlation spectroscopy. The strategy used enables one to make sequence-specific resonance assignments without reference to a known or assumed conformation of the DNA fragment.     

Fast evaluation of protein dynamics from deficient <Superscript>15</Superscript>N relaxation data

Simple and convenient method of protein dynamics evaluation from the insufficient experimental ¹⁵N relaxation data is presented basing on the ratios, products, and differences of longitudinal and transverse ¹⁵N relaxation rates obtained at a single magnetic field. Firstly, the proposed approach allows evaluating overall tumbling correlation time (nanosecond time scale). Next, local parameters of the model-free approach characterizing local mobility of backbone amide N–H vectors on two different time scales, S² and R _ex , can be elucidated. The generalized order parameter, S², describes motions on the time scale faster than the overall tumbling correlation time (pico- to nanoseconds), while the chemical exchange term, R _ex , identifies processes slower than the overall tumbling correlation time (micro- to milliseconds). Advantages and disadvantages of different methods of data handling are thoroughly discussed.     

Crankshaft motions of the polypeptide backbone in molecular dynamics simulations of human type-α transforming growth factor

Summary Order parameters for the backbone N–H and C–H bond vectors have been calculated from a 150 ps molecular dynamics (MD) simulation of human type- transforming growth factor in H₂O solvent. Two kinds of crankshaft motions of the polypeptide backbone are observed in this MD trajectory. The first involves small-amplitude rocking of the rigid peptide bond due to correlated changes in the backbone dihedral angles _i–1 and _i. These high-frequency librational crankshaft motions are correlated with systematically smaller values of motional order parameters for backbone N–H bond vectors compared to C–H bond vectors. In addition, infrequent crankshaft flips of the peptide bond from one local minimum to another are observed for several amino acid residues. These MD simulations demonstrate that comparisons of N–H and C–H order parameters provide a useful approach for identifying crank-shaft librational motions in proteins.     

Detecting antimicrobial peptides by exploring the mutual information of their sequences

Abstract

The rise of antibiotic resistance in pathogenic bacteria is a growing concern for every part of the world. The present study shows the prediction efficiency of mutual information for the classification of antimicrobial peptides. The proven role of antimicrobial peptides (AMPs) to fight against multidrug-resistant pathogens and AMP’s low toxic properties laid the foundation of computational methods to play their role in detecting AMPs from non-AMPs. Mutual information vectors (MIV) were created for AMP/non-AMP sequences and then fed to different machine learning classifiers out of which a random forest (RF) classifier showed best results for predicting AMPs. Random forest classifiers were evaluated on benchmark datasets by 10-fold cross-validation. The proposed MIV-RF method showed better prediction accuracy, MCC (Matthews correlation coefficient), and AUC-ROC (Area Under The Curve-Receiver Operating Characteristics) than available methods for detecting AMPs.

Communicated by Ramaswamy H. Sarma     

Proton and phosphorus magnetic relaxation studies on the dynamic structure of single-stranded polyriboadenylic acid

Proton and phosphorus-31 nuclear spin-lattice relaxation times (T₁) have been measured with the Fourier-transform method at 100 and 40.5 MHz, respectively, on single-stranded polyriboadenylic acid (poly(A)) in a neutral D₂O solution in the temperature range of 14–82°C. T₁ minimum is observed around 35–45°C for H(8), H(1′), and phosphorus resonances. Rotational correlation times have been deduced from the T₁ data, which indicate that the sugar–phosphate backbone as well as the base–sugar segment is undergoing rapid internal motion of 10^?8–10^?10 sec range. The molecular motion of the sugar–phosphate backbone as deduced from the phosphorus relaxation is well-characterized by a single activation enthalpy of 8.1 kal/mole for the whole temperature range of 14–82°C. Activation enthalpies of similar magnitude have been obtained for the motion of the adenine–ribose moiety from H(8) and H(1′) relaxation. The relative magnitude of T₁ for H(8) and H(1′) infers that the poly(A) nucleotide exists on the average as anti in the single-stranded form. The phosphorus T₁ value is consistent with a conformation such that both C(4′)–C(5′) and C(4′)–C(3′) bonds are nearly trans to their connected O–P bonds.     

Protein dynamics by solid-state nuclear magnetic resonance spectroscopy: Peptide backbone of the coat protein in fd bacteriophage

¹³C and ¹⁵N chemical shift anisotropy and ¹⁵N¹H dipolar powder patterns from backbone sites of the coat protein in fd bacteriophage are not averaged by motion. This means that the polypeptide backbone of the protein has no large amplitude motions rapid compared to 10⁴ Hz. Relaxation studies on the ¹³C_α and ¹⁵N amide resonances indicate the presence of motions on the 10⁹ Hz timescale. These results are reconciled with a model where an otherwise rigid backbone undergoes small amplitude, rapid motions.