Model-free analysis of protein backbone motion from residual dipolar couplings

On the basis of the measurement of NH residual dipolar couplings (RDCs) in 11 different alignment media, an RDC-based order parameter is derived for each residue in the protein ubiquitin. Dipolar couplings are motionally averaged in the picosecond to millisecond time range and, therefore, reflect motion slower than the inverse overall tumbling correlation time of the protein. It is found that there is considerable motion that is slower than the correlation time and could not be detected with previous NMR methodology. Amplitudes and anisotropies of the motion can be derived from the model-free analysis. The method can be applied provided that at least five sufficiently different alignment media can be found for the biomolecule under investigation.     

Solid-state 19F NMR spectroscopy reveals that Trp41 participates in the gating mechanism of the M2 proton channel of influenza A virus

The integral membrane protein M2 of influenza A virus assembles as a tetrameric bundle to form a proton-conducting channel that is activated by low pH. The side chain of His37 in the transmembrane alpha-helix is known to play an important role in the pH activation of the proton channel. It has also been suggested that Trp41, which is located in an adjacent turn of the helix, forms part of the gating mechanism. Here, a synthetic 25-residue peptide containing the M2 transmembrane domain was labeled with 6F-Trp41 and studied in lipid membranes by solid-state 19F NMR. We monitored the pH-dependent differences in the 19F dipolar couplings and motionally narrowed chemical shift anisotropies of this 6F-Trp41 residue, and we discuss the pH activation mechanism of the H+ channel. At pH 8.0, the structural parameters implicate an inactivated state, while at pH 5.3 the tryptophan conformation represents the activated state. With the aid of COSMOS force field simulations, we have obtained new side-chain torsion angles for Trp41 in the inactivated state (chi1 = -100 degrees +/- 10 degrees , chi2 = +110 degrees +/- 10 degrees ), and we predict a most probable activated state with chi1 = -50 degrees +/- 10 degrees and chi2 = +115 degrees +/- 10 degrees . We have also validated the torsion angles of His37 in the inactivated state as chi1 = -175 degrees +/- 10 degrees and chi2 = -170 degrees +/- 10 degrees .     

Contact model for the prediction of NMR N-H order parameters in globular proteins

An analytical relationship is presented for the estimation of NMR S2 order parameters of N-HN vectors of the protein backbone from high-resolution protein structures. The relationship solely depends on close contacts of the peptide plane to the rest of the protein. Application of the relationship to a number of proteins with high-resolution X-ray and NMR structures yields S2 values that are in good agreement with the ones determined from experimental relaxation data.     

Orientation and conformational preference of leucine-enkephalin at the surface of a hydrated dimyristoylphosphatidylcholine bilayer: NMR and MD simulation

The morphogenic opiate pentapeptide leucine-enkephalin (lenk) in a hydrated dimyristoylphosphatidylcholine (DMPC) bilayer is studied using NMR spectroscopy and molecular dynamics simulation. Contrary to the frequent assumption that the peptide attains a single fixed conformation in the presence of membranes, we find that the lenk molecule is flexible, switching between specific bent conformations. The constraints to the orientation of the aromatic rings that are identified by the NMR experiment are found by the MD simulation to be related to the depth of the peptide in the bilayer. The motion of the N-H vectors of the peptide bonds with respect to the magnetic field direction as observed by MD largely explain the magnitude of the observed residual dipolar coupling (RDC), which are much reduced over the static (15)N-(1)H coupling. The measured RDCs are nevertheless significantly larger than the predicted ones, possibly due the absence of long-time motions in the simulations. The conformational behavior of lenk at the DMPC surface is compared to that in the aqueous solution, both in the neutral and in the zwitterionic forms.     

UV resonance Raman measurements of poly-L-lysine's conformational energy landscapes: dependence on perchlorate concentration and temperature

UV resonance Raman spectroscopy has been used to determine the conformational energy landscape of poly-L-lysine (PLL) in the presence of NaClO4 as a function of temperature. At 1 degree C, in the presence of 0.83 M NaClO4, PLL shows an approximately 86% alpha-helix-like content, which contains alpha-helix and pi-bulge/helix conformations. The high alpha-helix-like content of PLL occurs because of charge screening due to strong ion-pair formation between ClO4- and the lysine side chain -NH3+. As the temperature increases from 1 to 60 degrees C, the alpha-helix and pi-bulge/helix conformations melt into extended conformations (PPII and 2.51-helix). We calculate the Psi Ramachandran angle distribution of the PLL peptide bonds from the UV Raman spectra which allows us to calculate the PLL (un)folding energy landscapes along the Psi reaction coordinate. We observe a basin in the Psi angle conformational space associated with alpha-helix and pi-bulge/helix conformations and another basin for the extended PPII and 2.51-helical conformations.     

Side chain orientation from methyl 1H-1H residual dipolar couplings measured in highly deuterated proteins

High-level deuteration is a prerequisite for the study of high molecular weight systems using liquid-state NMR. Here, we present new experiments for the measurement of proton-proton dipolar couplings in CH(2)D methyl groups of (13)C labeled, highly deuterated (70-80%) proteins. (1)H-(1)H residual dipolar couplings (RDCs) have been measured in two alignment media for 57 out of 70 possible methyl containing residues in the 167-residue flavodoxin-like domain of the E. coli sulfite reductase. These data yield information on the orientation of the methyl symmetry axis with respect to the molecular alignment frame. The alignment tensor characteristics were obtained very accurately from a set of backbone RDCs measured on the same protein sample. To demonstrate that accurate structural information is obtained from these data, the measured methyl RDCs for Valine residues are analyzed in terms of chi(1) torsion angles and stereospecific assignment of the prochiral methyl groups. On the basis of the previously determined backbone solution structure of this protein, the methyl RDC data proved sufficient to determine the chi(1) torsion angles in seven out of nine valines, assuming a single-rotamer model. Methyl RDCs are complementary to other NMR data, for example, methyl-methyl NOE, to determine side chain conformation in high molecular weight systems.     

Novel,highly polarizable thiophene derivatives for use in nonlinear optical applications

The synthesis of a variety of thiophene-containing materials for use in electro-optic devices requiring highly birefringent materials is reported. The refractive indices were measured with the use of an Abbé refractometer, and from these results the optical anisotropies, polarizabilities and order parameters were determined. The replacement of a phenyl ring by thiophene leads to large enhancements of polarizability. The most significant increases in polarizability anisotropy were observed when the rigid core was a collinear 5,5′-disubstituted-2,2′-dithienyl unit. The changes in the optical properties are discussed in terms of the structural units (thiophene, phenyl, isothiocyanate, nitrile, butylsulfanyl, alkoxy and alkyl moieties) used and their positions with respect to the molecular core.     

Enthalpic and entropic stages in alpha-helical peptide unfolding, from laser T-jump/UV Raman spectroscopy

The alpha-helix is a ubiquitous structural element in proteins, and a number of studies have addressed the mechanism of helix formation and melting in simple peptides. However, fundamental issues remain to be resolved, particularly the temperature (T) dependence of the rate. In this work, we report application of a novel kHz repetition rate solid-state tunable NIR (pump) and deep UV Raman (probe) laser system to study the dynamics of helix unfolding in Ac-GSPEA3KA4KA4-CO-D-Arg-CONH2, a peptide designed for helix stabilization in aqueous solution. Its T-dependent UV resonance Raman (UVRR) spectra, excited at 197 nm for optimal enhancement of amide vibrations, were decomposed into variable contributions from helix and coil spectra. The helix fractions derived from the UVRR spectra and from far UV CD spectra were coincident at low T but deviated increasingly at high T, the UVRR curve giving higher helix content. This difference is consistent with the greater sensitivity of UVRR spectra to local conformation than CD. After a laser-induced T-jump, the UVRR-determined helix fractions defined monoexponential decays, with time-constants of approximately 120 ns, independent of the final T (Tf = 18-61 degrees C), provided the initial T (Ti) was held constant (6 degrees C). However, there was also a prompt loss of helicity, whose amplitude increased with increasing Tf, thereby defining an initial enthalpic phase, distinct from the subsequent entropic phase. These phases are attributed to disruption of H-bonds followed by reorientation of peptide links, as the chain is extended. When Ti was raised in parallel with Tf (10 degrees C T-jumps), the prompt phase merged into an accelerating slow phase, an effect attributable to the shifting distribution of initial helix lengths. Even greater acceleration with rising Ti has been reported in T-jump experiments monitored by IR and fluorescence spectroscopies. This difference is attributable to the longer range character of these probes, whose responses are therefore more strongly weighted toward the H-bond-breaking enthalpic process.     

Periodicity in residual dipolar couplings and nucleic acid structures

The periodicity in nucleic acid duplex structures is shown to be correlated to the periodicity in residual dipolar couplings (RDCs) in the form of an "RDC wave". This "RDC wave" is characteristic of the alignment of the duplex in the magnetic field, and hence fitting of the data allows the duplex global orientation (, Phi) to be extracted. Further, because the "RDC wave" is fit as a data set of a corresponding secondary structure element, the degeneracy problem is greatly reduced. Consequently, with the global orientation (, Phi) determined, local bond vector conformations are defined. The fit is demonstrated in the examples of the imino RDCs of the negative regulator of splicing RNA fragment (NRS23) and for the C1'H1' RDCs of the Dickerson dodecamer.     

Chirality and helix stability of polyglutamic acid enantiomers

In this work the chirality and the relative thermal stability of ordered micellar aggregates of poly-L- and poly-D-glutamic acids with the cationic surfactant C14TAB is examined. The complexed mesophases poly-L-Glu/C14TAB and poly-d-Glu/C14TAB were characterized by circular dichroism (CD) in the temperature range 10-70 degrees C for their chirality and thermal stability as well as by X-ray diffraction (XRD) for the micellar ordered structure. Low angle XRD analysis showed that both micellar aggregates poly-L-Glu/C14TAB and poly-D-Glu/C14TAB are hexagonally packed in a MCM-41 fashion with an intermicellar distance identical and equal to 3.55+/-0.10 nm. The CD spectra indicated that both complexes poly-L-Glu/C14TAB and poly-D-Glu/C14TAB possess a mainly alpha-helix structure and are exact mirror images to each other. The same mirror images and a mainly alpha-helix configuration were also observed by CD for the free poly-l- and poly-d-glutamic acids at room temperature. As the temperature increases from 10 up to 70 degrees C the alpha-helix of the poly-l-glutamic acid is gradually transformed to a mixture containing increased amounts of the 3(10)-helix while the alpha-helix structure of the poly-d-glutamic acid is constantly degenerated. In contrast the alpha-helices of the corresponding complexes poly-L-Glu/C14TAB and poly-d-Glu/C14TAB are degenerated upon heating without appreciable increase of the 3(10)-helices as an intermediate configuration. This difference in helix conservation is attributed to increase protection of the l-enantiomers, compared to d-enantiomers, which might be related to the survival of l-aminoacids in the living world.     

RDC‐enhanced structure calculation of a β‐heptapeptide in methanol

Residual dipolar couplings (RDCs) are a rich source of structural information that goes beyond the range covered by the nuclear Overhauser effect or scalar coupling constants. They can only be measured in partially oriented samples. RDC studies of peptides in organic solvents have so far been focused on samples in chloroform or DMSO. Here, we show that stretched poly(vinyl acetate) can be used for the partial alignment of a linear β‐peptide with proteinogenic side chains in methanol. ¹D_CH, ¹D_NH, and ²D_HH RDCs were collected with this sample and included as restraints in a simulated annealing calculation. Incorporation of RDCs in the structure calculation process improves the long‐range definition in the backbone of the resulting 3₁₄‐helix and uncovers side‐chain mobility. Experimental side‐chain RDCs of the central leucine and valine residues are in good agreement with predicted values from a local three‐state model. Copyright © 2016 John Wiley & Sons, Ltd.     

Rotor modulations and recoupling strategies in 13C solid-state magic-angle-spinning NMR spectroscopy: probing molecular orientation and dynamics.

Recoupling strategies for anisotropic interactions enable the investigation of molecular structure, order and dynamics in a sensitive and site-specific fashion by solid-state NMR spectroscopy. Whereas magic-angle spinning (MAS) efficiently averages anisotropic interactions and enhances spectral resolution, recoupling pulse sequences selectively restore certain parts of rotor-modulated dipole-dipole couplings or chemical shift anisotropies (CSA). More specifically, it is possible to recouple either the omegaR- or the 2omegaR-modulated terms of an interaction Hamiltonian, which exhibit different orientation dependencies and, in this way, provide a means of distinguishing whether the observed NMR spectra are affected by molecular motion or by molecular orientation. Sideband patterns generated by reconversion rotor encoding allow for a precise and selective determination of coupling constants and anisotropies, which contain site-specific information on structure, orientation and/or dynamics of individual molecular segments. Corresponding recoupling schemes are presented in a common context, and the possibilities of exploiting these effects for the determination of order parameters of oriented materials, such as oriented polymer chains or extruded fibres of a discotic mesogen, are discussed. The obtained orientational order parameters are compared to results from two-dimensional wide angle X-ray scattering (WAXS).     

Residual-Chemical-Shift-Anisotropy-Based Enantiodifferentiation in Lyotropic Liquid Crystalline Phases Based on Helically Chiral Polyacetylenes

Anisotropic NMR spectroscopy, revealing residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs) has emerged as a powerful tool to determine the configurations of synthetic and complex natural compounds. The deduction of the absolute in addition to the relative configuration is one of the primary goals in the field. Therefore, the investigation of the enantiodiscriminating capabilities of chiral alignment media becomes essential. While RDCs and RCSAs are now used for the determination of the relative configuration routinely, RCSAs have not been measured in chiral alignment media such as chiral liquid crystals. Herein, we present this application by measuring RCSAs for chiral analytes such as indanol and isopinocampheol in the lyotropic liquid crystalline phase of an L-valine derived helically chiral polyacetylenes. We have also demonstrated that a single 1D ¹³C−{¹H} NMR spectrum suffices to get the RCSAs circumventing the necessity to acquire two spectra at two alignment conditions.     

Limited variations in 15N CSA magnitudes and orientations in ubiquitin are revealed by joint analysis of longitudinal and transverse NMR relaxation

The site-specific magnitudes and orientations of the chemical shift tensors have been estimated for 70 backbone (15)N-nuclei in human ubiquitin from the field dependence of dynamic independent ratios between relaxation rates, both longitudinal and transverse, measured at 9.4, 11.7, 14.1, and 18.8 T. The results were jointly analyzed with previously published relaxation data [Fushman; Tjandra; Cowburn. J.Am. Chem. Soc. 1998, 120, 10947-10952] [Kover; Batta. J. Mag. Reson. 2001, 150, 137-146]. The effective magnitudes of the anisotropies distribute around 169 ppm with a variability of 5 ppm. The orientation factors, reflecting the orientation of the CSA relative to the NH bond, distribute around -0.80 with a variability of 0.04, which corresponds to an angle between the symmetry axis of an assumed axially symmetric shielding tensor and the NH bond of 21.4 degrees, and a variability of 2.3 degrees. Correlations with the isotropic (15)N-chemical shifts are observed. Variations in the shielding anisotropies add uncertainty to the obtained order parameters proportional to the square of the magnetic field, when data are analyzed using an assumed invariant CSA tensor for all sites. Around 3% additional uncertainty in the order parameters for 800 MHz data is expected. The optimal TROSY field for amide nitrogen TROSY is estimated, with only marginal variations due to site-to-site variations. Variations in the shielding tensors add uncertainty to the exchange terms calculated from cross-correlation rates. An approach for estimating the exchange terms is suggested, where the uncertainty due to CSA-variations is minimized.     

Hydrodynamic and conformational properties of cellulose valerate molecules in dilute solution

Using the methods of molecular hydrodynamics, structural-conformational studies have been performed for a number of cellulose valerate and acetovalerate samples in the molecular-mass range M = (58.1–464.3) × 10³ and with a mean degree of substitution of 182.4 with respect to valeric acid isomers. The conformations of valerate-substituted cellulose molecules are characterized by an increased local packing density of monomer units. The molecular conformations are quantitatively described in terms of the helix formed by the succession of vectors connecting glycoside oxygens of a chain. The molecular-hydrodynamic and conformational characteristics of cellulose valerate are compared with the corresponding characteristics of cellulose myristate.     

15N T2' relaxation times of bacteriorhodopsin transmembrane amide nitrogens

15N T(2)' relaxation times of bacteriorhodopsin (BR) amide nitrogens were determined in the temperature range from 40 to -60 degrees C using a Hahn echo pulse sequence and proton decoupling during the echo and detection times. Using oriented membrane samples, with their bilayer normal parallel to the external magnetic field, the (15)N amide nitrogens belonging to the transmembrane helices could be selected for the analysis. The experiments were performed on purple membrane fragments (in which BR is organized in a 2D crystalline network) and on BR reconstituted into dimyristoylphosphatidylcholine at a 1:150 molar ratio (in which BR is in a freely diffusing monomeric state at 40 degrees C and in an aggregated state at 4 degrees C and below). The results are discussed in terms of helix dynamics, mosaic spread and resolution of the (15)N spectra for the various samples and experimental conditions.     

Sequence-dependent bending of DNA induced by cisplatin: NMR structures of an A.T-rich 14-mer duplex

The NMR solution structure of the A.T rich DNA 14-mer duplex d(ATACATGGTACATA).d(TATGTACCATGTAT) is reported. This is compared with the NMR structure of the same duplex intrastrand cross-linked at the d(G*pG*) site by cis-(Pt(NH3)2?2+, derived from the anticancer drug cisplatin. The unmodified duplex has B-DNA geometry, but there is a large positive base-pair roll (roll angle 24 +/- 2 degrees) at the T9-A10 step on the 3' side of the central GG site. Platination of the DNA duplex causes the adjacent guanine bases to roll toward one another (roll angle 44 +/- 4 degrees), leading to an overall helix bend of 52 +/- 9 degrees. The platinum atom is displaced from the planes of the coordinated G7* and G8* by 0.8 A and 0.3 A, respectively. The minor groove opposite the platinum lesion is widened and flattened, with geometric parameters similar to those of A-form DNA. The unwinding of the helix at the platination site is 26 degrees. Platination causes the DNA duplex to bend toward the 3'-end (with respect to the G*G* strand), in contrast to G C-rich structures reported previously, which bend toward the 5'-end. This difference can be attributed to the predisposition of the A.T rich duplex toward bending in this region. Protein recognition of bent platinated G*G* lesions may therefore exhibit a strong dependence on the local DNA structure.     

Cross-linking-induced permanently perpendicular helix orientation in surface-grafted polyglutamate films

Using a chemical cross-linking procedure, surface-grafted polyglutamate films with a permanently perpendicular helix orientation were prepared. A surface-grafted alpha-helical polyglutamate film containing polymerizable side groups was synthesized by ring-opening terpolymerization of 50 molar% gamma-methyl-L-glutamate N-carboxyanhydride (NCA), 30% gamma-stearyl-L-glutamate NCA and 20% gamma-4-vinylbenzyl-L-glutamate NCA initiated from a silicon substrate functionalized with primary amino groups. The average tilt angle of the end-grafted helices in this film is approximately 66 degrees , indicating a nearly parallel helix orientation with respect to the substrate surface. After swelling of the grafted terpolyglutamate film in stearyl methacrylate and subsequent radical cross-linking, the average helix tilt angle decreases to about 11 degrees, indicating an almost perpendicular helix orientation. The film thickness increases accordingly from 151 A before to approximately 390 A after cross-linking. Extensive solvent treatment of the cross-linked film shows that the perpendicular helix orientation is permanent.     

Simultaneous assignment of all diastereotopic protons in strychnine using RDCs: PELG as alignment medium for organic molecules

The concept of using residual dipolar couplings (RDCs) for the structure determination of organic molecules is applied to the simultaneous assignment of all diastereotopic protons in strychnine. To use this important NMR parameter the molecule has to be aligned in the magnetic field. Here we present a new alignment medium for organic substrates. The optimization of the alignment properties of mixtures of poly-gamma-ethyl-L-glutamate (PELG) and CDCl(3) are described and the alignment properties of PELG at different concentrations are evaluated. A comparison of PELG with poly-gamma-benzyl-L-glutamate (PBLG) shows considerable differences in the magnitude of alignment for strychnine in the two alignment media. PELG induces a lower degree of order and makes the measurement of residual dipolar couplings (RDCs) in strychnine possible. All one-bond C-H RDCs of strychnine in PELG were determined by using 2D heteronuclear single quantum coherence (HSQC) spectroscopy. The strategy for the extraction of RDCs for methylene groups is described in detail. The RDCs and order parameters are used to assign pairs of diastereotopic protons. This methodology can distinguish not only one pair of diastereotopic protons but it can be used to assign all pairs of diastereotopic protons simultaneously. Two different calculation approaches to achieve this task are described in detail.