Structure and Dynamics of Coil-like Molecules by Residual Dipolar Couplings

Nuclear magnetic resonance (NMR) spectroscopy provides us with many means to study biological macromolecules in solution. Proteins in particular are the most intriguing targets for NMR studies. Protein functions are usually ascribed to specific three-dimensional structures but more recently tails, long loops and non-structural polypeptides have also been shown to be biologically active. Examples include prions, -synuclein, amylin and the NEF HIV-protein. However, conformational preferences in coil-like molecules are difficult to study by traditional methods. Residual dipolar couplings (RDCs) have opened up new opportunities; however their analysis is not trivial. Here we show how to interpret RDCs from these weakly structured molecules. The most notable residual dipolar couplings arise from steric obstruction effects. In dilute liquid crystalline media as well as in anisotropic gels polypeptides encounter nematogens. The shape of a polypeptide conformation limits the encounter with the nematogen. The most elongated conformations may come closest whereas the most compact remain furthest away. As a result there is slightly more room in the solution for the extended than for the compact conformations. This conformation-dependent concentration effect leads to a bias in the measured data. The measured values are not arithmetic averages but essentially weighted averages over conformations. The overall effect can be calculated for random flight chains and simulated for more realistic molecular models. Earlier there was an implicit thought that weakly structured or non-structural molecules would not yield to any observable residual dipolar couplings. However, in the pioneering study by Shortle and Ackerman RDCs were clearly observed. We repeated the study for urea-denatured protein at high temperature and also observed indisputably RDCs. This was very convincing to us but we could not possibly accept the proposed reason for the non-zero RDCs, namely that there would be some residual structure left in the protein that to our understanding was fully denatured. We proceeded to gain understanding via simulations and elementary experiments. In measurements we used simple homopolymers with only two labelled residues and we simulated the data to learn more about the origin of RDCs. We realized that RDCs depend on the position of the residue as well as on the length of the polypeptide. Investigations resulted in a theoretical model for RDCs from coil-like molecules. Later we extended the studies by molecular dynamics. Somewhat surprisingly the effects are small for non-structured molecules whereas the bias may be large for a small compact protein. All in all the work gave clear and unambiguous results on how to interpret RDCs as structural and dynamic parameters of weakly structured proteins.

Dynamic nature of proteins - interpretation of residual dipolar couplings

Protein conformations and dynamics can be studied by nuclear magnetic resonance spectroscopy using dilute liquid crystalline samples. This work clarifies the interpretation of residual dipolar coupling data yielded by the experiments. It was discovered that unfolded proteins without any additional structure beyond that of a mere polypeptide chain exhibit residual dipolar couplings. Also, it was found that molecular dynamics induce fluctuations in the molecular alignment and doing so affect residual dipolar couplings. The finding clarified the origins of low order parameter values observed earlier. The work required the development of new analytical and computational methods for the prediction of intrinsic residual dipolar coupling profiles for unfolded proteins. The presented characteristic chain model is able to reproduce the general trend of experimental residual dipolar couplings for denatured proteins. The details of experimental residual dipolar coupling profiles are beyond the analytical model, but improvements are proposed to achieve greater accuracy. A computational method for rapid prediction of unfolded protein residual dipolar couplings was also developed. Protein dynamics were shown to modulate the effective molecular alignment in a dilute liquid crystalline medium. The effects were investigated from experimental and molecular dynamics generated conformational ensembles of folded proteins. It was noted that dynamics induced alignment is significant especially for the interpretation of molecular dynamics in small, globular proteins. A method of correction was presented. Residual dipolar couplings offer an attractive possibility for the direct observation of protein conformational preferences and dynamics. The presented models and methods of analysis provide significant advances in the interpretation of residual dipolar coupling data from proteins.

2(n)-SEMA-a robust solid state nuclear magnetic resonance experiment for measuring heteronuclear dipolar couplings in static oriented systems using effective transverse spin-lock

Separated local field (SLF) spectroscopy is a powerful technique to measure heteronuclear dipolar couplings. The method provides site-specific dipolar couplings for oriented samples such as membrane proteins oriented in lipid bilayers and liquid crystals. A majority of the SLF techniques utilize the well-known Polarization Inversion Spin Exchange at Magic Angle (PISEMA) pulse scheme which employs spin exchange at the magic angle under Hartmann-Hahn match. Though PISEMA provides a relatively large scaling factor for the heteronuclear dipolar coupling and a better resolution along the dipolar dimension, it has a few shortcomings. One of the major problems with PISEMA is that the sequence is very much sensitive to proton carrier offset and the measured dipolar coupling changes dramatically with the change in the carrier frequency. The study presented here focuses on modified PISEMA sequences which are relatively insensitive to proton offsets over a large range. In the proposed sequences, the proton magnetization is cycled through two quadrants while the effective field is cycled through either two or four quadrants. The modified sequences have been named as 2(n)-SEMA where n represents the number of quadrants the effective field is cycled through. Experiments carried out on a liquid crystal and a single crystal of a model peptide demonstrate the usefulness of the modified sequences. A systematic study under various offsets and Hartmann-Hahn mismatch conditions has been carried out and the performance is compared with PISEMA under similar conditions.

Enantiodiscrimination and extraction of short and long range homo- and hetero-nuclear residual dipolar couplings by a spin selective correlation experiment

A two dimensional correlation experiment for the measurement of short and long range homo- and hetero- nuclear residual dipolar couplings (RDCs) from the broad and featureless proton NMR spectra including C-13 satellites is proposed. The method employs a single natural abundant C-13 spin as a spy nucleus to probe all the coupled protons and permits the determination of RDCs of negligible strengths. The technique has been demonstrated for the study of organic chiral molecules aligned in chiral liquid crystal, where additional challenge is to unravel the overlapped spectrum of enantiomers. The significant advantage of the method is demonstrated in better chiral discrimination using homonuclear RDCs as additional parameters. (C) 2010 Elsevier B.V. All rights reserved.

Spin-Selective Correlation Experiment for Measurement of Long-Range J Couplings and for Assignment of (R/S) Enantiomers from the Residual Dipolar Couplings and DFT

We report the C-HETSERF experiment for determination of long- and short-range homo- and heteronuclear scalar couplings ((n)J(HH) and (n)J(XH), n >= 1) of organic molecules with a low sensitivity dilute heteronucleus in natural abundance. The method finds significant advantage in measurement of relative signs of long-range heteronuclear total couplings in chiral organic liquid crystal. The advantage of the method is demonstrated for the measurement of residual dipolar couplings (RDCs) in enantiomers oriented in the chiral liquid crystal with a focus to unambiguously assign R/S designation in a 2D spectrum. The alignment tensor calculated from the experimental RDCs and with the computed structures of enantiomers obtained by DFT calculations provides the size of the back-calculated RDCs. Smaller root-mean-square deviations (rmsd) between experimental and calculated RDCs indicate better agreement with the input structure and its correct designation of the stereogenic center.

Developments in quantum information processing by nuclear magnetic resonance: Use of quadrupolar and dipolar couplings

Use of dipolar and quadrupolar couplings for quantum information processing (QIP) by nuclear magnetic resonance (NMR) is described. In these cases, instead of the individual spins being qubits, the 2(n) energy levels of the spin-system can be treated as an n-qubit system. It is demonstrated that QIP in such systems can be carried out using transition-selective pulses, in (CHCN)-C-3, (CH3CN)-C-13, Li-7 (I = 3/2) and Cs-133 (I = 7/2), oriented in liquid crystals yielding 2 and 3 qubit systems. Creation of pseudopure states, implementation of logic gates and arithmetic operations (half-adder and subtractor) have been carried out in these systems using transition-selective pulses.

C-13-H-1 dipolar couplings for probing rod-like hydrogen bonded mesogens

Hydrogen bonding is the most important non-covalent interaction utilised in building supramolecular assemblies and is preferred often as a means of construction of molecular, oligomeric as well as polymeric materials that show liquid crystalline properties. In this work, a pyridine based nematogenic acceptor has been synthesized and mixed with non-mesogenic 4-methoxy benzoic acid to get a hydrogen bonded mesogen. The existence of hydrogen bonding between the pyridyl unit and the carboxylic acid was established using FT-IR spectroscopy from the observation of characteristic stretching vibrations of unionized type at 2425 and 1927 cm(-1). The mesogenic acceptor and the complex have been investigated using C-13 NMR in solution, solid and liquid crystalline states. Together with the 2D separated local field NMR experiments, the studies confirm the molecular structure in the mesophase and yield the local orientational order parameters. It is observed that the insertion of 4-methoxy benzoic acid not only enhances the mesophase stability but also induces a smectic phase due to an increase in the core length of the hydrogen bonded mesogen.

Molecular topology of three ring nematogens from C-13-H-1 dipolar couplings

Molecules in their liquid crystalline phase undergo rotational motion about the long axis of the molecule and the shape adopted by the rotating molecule plays an important role in influencing the mesophase morphology. In this context, obtaining the topology and the relative orientation of the different sub-units are important steps. For studying the liquid crystalline phase, C-13 NMR spectroscopy is a convenient method and for certain specifically designed nematogens, 2-dimensional separated local field (2D-SLF) NMR spectroscopy provides a particularly simple and straightforward means of arriving at the molecular topology. We demonstrate this approach on two three ring based nematogens designed with a phenyl or a thiophene ring at one of the termini. From the C-13-H-1 dipolar couplings of the terminal carbon obtained using the 2D-SLF NMR technique, the order parameter of the local symmetry axis of the terminal phenyl ring as well as of the long molecular axis could be easily estimated. For the thiophene nematogen, the lack of symmetry of the thiophene moiety necessitates some additional computational steps. The results indicate that the thiophene unit has its local ordering axis oriented away from the long molecular axis by a small angle, consistent with a bent structure expected in view of the thiophene geometry. The experiment also demonstrates the ability of 2D-SLF NMR to provide high resolution spectra by separation of several overlapped resonances in terms of their C-13-H-1 dipolar couplings. The results are consistent with a rod-like topology of the core of the investigated mesogens. The investigation demonstrates the potential of 2D-SLF NMR C-13 spectroscopy for obtaining atomistic level information and its utility for topological studies of different mesogens.

CH-RES-TOCSY: Enantiomers spectral resolution and measurement of heteronuclear residual dipolar couplings

A new 2D NMR technique cited as CH-RES-TOCSY, for complete unraveling the spectra of enantiomers and for the measurement of structurally important C-HRDC sis reported. The spectral overlap and complexity of peaks were reduced by the blend of selective excitation and homo-decoupling. Differential values of C-H RDCs of enantiomers (R and S) are exploited to separate the enantiomeric peaks. The complete unraveling of the spectra of both the enantiomers is achieved by incorporating a TOCSY mixing blockprior to signal acquisition. The additional application of the method is demonstrated for the assignment of symmetric isomers. (C) 2015 Elsevier B.V. All rights reserved.

Measurement of Large Dipolar Couplings of a Liquid Crystal with Terminal Phenyl Rings and Estimation of the Order Parameters

NMR spectroscopy is a powerful means of studying liquid-crystalline systems at atomic resolutions. Of the many parameters that can provide information on the dynamics and order of the systems, H-1-C-13 dipolar couplings are an important means of obtaining such information. Depending on the details of the molecular structure and the magnitude of the order parameters, the dipolar couplings can vary over a wide range of values. Thus the method employed to estimate the dipolar couplings should be capable of estimating both large and small dipolar couplings at the same time. For this purpose, we consider here a two-dimensional NMR experiment that works similar to the insensitive nuclei enhanced by polarization transfer (INEPT) experiment in solution. With the incorporation of a modification proposed earlier for experiments with low radio frequency power, the scheme is observed to enable a wide range of dipolar couplings to be estimated at the same time. We utilized this approach to obtain dipolar couplings in a liquid crystal with phenyl rings attached to either end of the molecule, and estimated its local order parameters.

13C NMR of a single molecule magnet: analysis of pseudocontact shifts and residual dipolar couplings

Paramagnetic triple decker complexes of lanthanides are promising Single Molecule Magnets (SMMs), with many potential uses. Some of them show preferable relaxation behavior, which enables the recording of well resolved NMR spectra. These axially symmetric complexes are also strongly magnetically anisotropic, and this property can be described with the axial component of the magnetic susceptibility tensor, χa. For triple decker complexes with phthalocyanine based ligands, the Fermi˗contact contribution is small. Hence, together with the axial symmetry, the experimental chemical shifts in 1H and 13C NMR spectra can be modeled easily by considering pseudocontact and orbital shifts alone. This results in the determination of the χa value, which is also responsible for molecular alignment and consequently for the observation of residual dipolar couplings (RDCs). A detailed analysis of the experimental 1H-13C and 1H-1H couplings revealed that contributions from RDCs (positive and negative) and from dynamic frequency shifts (negative for all observed couplings) have to be considered. Whilst the pseudocontact shifts depend on the average positions of 1H and 13C nuclei relative to the lanthanide ions, the RDCs are related to the mobility of nuclei they correspond to. This phenomenon allows for the measurement of the internal mobility of the various groups in the SMMs.

On Provable Algorithms for Determination of Continuous Protein Interdomain Motions from Residual Dipolar Couplings

Dynamics of biomolecules over various spatial and time scales are essential for biological functions such as molecular recognition, catalysis and signaling. However, reconstruction of biomolecular dynamics from experimental observables requires the determination of a conformational probability distribution. Unfortunately, these distributions cannot be fully constrained by the limited information from experiments, making the problem an ill-posed one in the terminology of Hadamard. The ill-posed nature of the problem comes from the fact that it has no unique solution. Multiple or even an infinite number of solutions may exist. To avoid the ill-posed nature, the problem needs to be regularized by making assumptions, which inevitably introduce biases into the result.

Here, I present two continuous probability density function approaches to solve an important inverse problem called the RDC trigonometric moment problem. By focusing on interdomain orientations we reduced the problem to determination of a distribution on the 3D rotational space from residual dipolar couplings (RDCs). We derived an analytical equation that relates alignment tensors of adjacent domains, which serves as the foundation of the two methods. In the first approach, the ill-posed nature of the problem was avoided by introducing a continuous distribution model, which enjoys a smoothness assumption. To find the optimal solution for the distribution, we also designed an efficient branch-and-bound algorithm that exploits the mathematical structure of the analytical solutions. The algorithm is guaranteed to find the distribution that best satisfies the analytical relationship. We observed good performance of the method when tested under various levels of experimental noise and when applied to two protein systems. The second approach avoids the use of any model by employing maximum entropy principles. This 'model-free' approach delivers the least biased result which presents our state of knowledge. In this approach, the solution is an exponential function of Lagrange multipliers. To determine the multipliers, a convex objective function is constructed. Consequently, the maximum entropy solution can be found easily by gradient descent methods. Both algorithms can be applied to biomolecular RDC data in general, including data from RNA and DNA molecules.