Intermediate motions and dipolar couplings as studied by Lee-Goldburg cross-polarization NMR: Hartmann-Hahn matching profiles

In this article, we evaluate the use of simple Lee-Goldburg cross-polarization (LG-CP) NMR experiments for obtaining quantitative information of molecular motion in the intermediate regime. In particular, we introduce the measurement of Hartmann-Hahn matching profiles for the assessment of heteronuclear dipolar couplings as well as dynamics as a reliable and robust alternative to the more common analysis of build-up curves. We have carried out dynamic spin dynamics simulations in order to test the method's sensitivity to intermediate motion and address its limitations concerning possible experimental imperfections. We further demonstrate the successful use of simple theoretical concepts, most prominently Anderson-Weiss (AW) theory, to analyze the data. We further propose an alternative way to estimate activation energies of molecular motions, based upon the acquisition of only two LG-CP spectra per temperature at different temperatures. As experimental tests, molecular jumps in imidazole methyl sulfonate, trimethylsulfoxonium iodide, and bisphenol A polycarbonate were investigated with the new method.

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.

Fast Structure-Based Assignment of 15N HSQC Spectra of Selectively 15N-Labeled Paramagnetic Proteins

A novel strategy for fast NMR resonance assignment of N-15 HSQC spectra of proteins is presented. It requires the structure coordinates of the protein, a paramagnetic center, and one or more residue-selectively N-15-labeled samples. Comparison of sensitive undecoupled N-15 HSQC spectra recorded of paramagnetic and diamagnetic samples yields data for every cross-peak on pseudocontact shift, paramagnetic relaxation enhancement, cross-correlation between Curie-spin and dipole-dipole relaxation, and residual dipolar coupling. Comparison of these four different paramagnetic quantities with predictions from the three-dimensional structure simultaneously yields the resonance assignment and the anisotropy of the susceptibility tensor of the paramagnetic center. The method is demonstrated with the 30 kDa complex between the N-terminal domain of the epsilon subunit and the theta subunit of Escherichia Coll DNA polymerase III. The program PLATYPUS was developed to perform the assignment, provide a measure of reliability of the assignment, and determine the susceptibility tensor anisotropy.

Application of NMR RDC’s, relaxation and self-diffusion for the study of dynamic processes of small molecules in solution

Spin-lattice Relaxation, self-Diffusion coefficients and Residual Dipolar Couplings (RDC’s) are the basis of well established Nuclear Magnetic Resonance techniques for the physicochemical study of small molecules (typically organic compounds and natural products with MW < 1000 Da), as they proved to be a powerful and complementary source of information about structural dynamic processes in solution. The work developed in this thesis consists in the application of the earlier-mentioned NMR techniques to explore, analyze and systematize patterns of the molecular dynamic behavior of selected small molecules in particular experimental conditions. Two systems were chosen to investigate molecular dynamic behavior by these techniques: the dynamics of ion-pair formation and ion interaction in ionic liquids (IL) and the dynamics of molecular reorientation when molecules are placed in oriented phases (alignment media). The application of NMR spin-lattice relaxation and self-diffusion measurements was applied to study the rotational and translational molecular dynamics of the IL: 1-butyl-3-methylimidazolium tetrafluoroborate [BMIM][BF4]. The study of the cation-anion dynamics in neat and IL-water mixtures was systematically investigated by a combination of multinuclear NMR relaxation techniques with diffusion data (using by H1, C13 and F19 NMR spectroscopy). Spin-lattice relaxation time (T1), self-diffusion coefficients and nuclear Overhauser effect experiments were combined to determine the conditions that favor the formation of long lived [BMIM][BF4] ion-pairs in water. For this purpose and using the self-diffusion coefficients of cation and anion as a probe, different IL-water compositions were screened (from neat IL to infinite dilution) to find the conditions where both cation and anion present equal diffusion coefficients (8% water fraction at 25 ºC). This condition as well as the neat IL and the infinite dilution were then further studied by 13C NMR relaxation in order to determine correlation times (c) for the molecular reorientational motion using a mathematical iterative procedure and experimental data obtained in a temperature range between 273 and 353 K. The behavior of self-diffusion and relaxation data obtained in our experiments point at the combining parameters of molar fraction 8 % and temperature 298 K as the most favorable condition for the formation of long lived ion-pairs. When molecules are subjected to soft anisotropic motion by being placed in some special media, Residual Dipolar Couplings (RDCs), can be measured, because of the partial alignment induced by this media. RDCs are emerging as a powerful routine tool employed in conformational analysis, as it complements and even outperforms the approaches based on the classical NMR NOE or J3 couplings. In this work, three different alignment media have been characterized and evaluated in terms of integrity using 2H and 1H 1D-NMR spectroscopy, namely the stretched and compressed gel PMMA, and the lyotropic liquid crystals CpCl/n-hexanol/brine and cromolyn/water. The influence that different media and degrees of alignment have on the dynamic properties of several molecules was explored. Different sized sugars were used and their self-diffusion was determined as well as conformation features using RDCs. The results obtained indicate that no influence is felt by the small molecules diffusion and conformational features studied within the alignment degree range studied, which was the 3, 5 and 6 % CpCl/n-hexanol/brine for diffusion, and 5 and 7.5 % CpCl/n-hexanol/brine for conformation. It was also possible to determine that the small molecules diffusion verified in the alignment media presented close values to the ones observed in water, reinforcing the idea of no conditioning of molecular properties in such media.

Applying the chain formation model to magnetic properties of aggregated ferrofluids

The magnetization properties of aggregated ferrofluids are calculated by combining the chain formation model developed by Zubarev with the modified mean-field theory. Using moderate assumptions for the inter- and intrachain interactions we obtain expressions for the magnetization and initial susceptibility. When comparing the results of our theory to molecular dynamics simulations of the same model we find that at large dipolar couplings (lambda>3) the chain formation model appears to give better predictions than other analytical approaches. This supports the idea that chain formation is an important structural ingredient of strongly interacting dipolar particles.

Molecular dynamics study on the equilibrium magnetization properties and structure of ferrofluids

We investigate in detail the initial susceptibility, magnetization curves, and microstructure of ferrofluids in various concentration and particle dipole moment ranges by means of molecular dynamics simulations. We use the Ewald summation for the long-range dipolar interactions, take explicitly into account the translational and rotational degrees of freedom, coupled to a Langevin thermostat. When the dipolar interaction energy is comparable with the thermal energy, the simulation results on the magnetization properties agree with the theoretical predictions very well. For stronger dipolar couplings, however, we find systematic deviations from the theoretical curves. We analyze in detail the observed microstructure of the fluids under different conditions. The formation of clusters is found to enhance the magnetization at weak fields and thus leads to a larger initial susceptibility. The influence of the particle aggregation is isolated by studying ferro-solids, which consist of magnetic dipoles frozen in at random locations but which are free to rotate. Due to the artificial suppression of clusters in ferrosolids the observed susceptibility is considerably lowered when compared to ferrofluids.

A new member of the ribbon-helix-helix transcription factor superfamily from the plant pathogen Xanthomonas axonopodis pv. citri

XACb0070 is an uncharacterized protein coded by the two large plasmids isolated from Xanthomonas axonopodis pv. cirri, the agent of citrus canker and responsible for important economical losses in citrus world production. XACb0070 presents sequence homology only with other hypothetical proteins belonging to plant pathogens, none of which have their structure determined. The NMR-derived solution structure reveals this protein is a homodimer in which each monomer presents two domains with different structural and dynamic properties: a folded N-terminal domain with beta alpha alpha topology which mediates dimerization and a long disordered C-terminal tail. The folded domain shows high structural similarity to the ribbon-helix-helix transcriptional repressors, a family of DNA-binding proteins of conserved 3D fold but low sequence homology: indeed XACb0070 binds DNA. Primary sequence and fold comparison of XACb0070 with other proteins of the ribbon-helix-helix family together with examination of the genes in the vicinity of xacb0070 suggest the protein might be the component of a toxin-antitoxin system. (C) 2010 Elsevier Inc. All rights reserved.

Recoupled separated-local-field experiments and applications to study intermediate-regime molecular motions

A specific separated-local-field NMR experiment, dubbed Dipolar-Chemical-Shift Correlation (DIPSHIFT) is frequently used to study molecular motions by probing reorientations through the changes in XH dipolar coupling and T-2. In systems where the coupling is weak or the reorientation angle is small, a recoupled variant of the DIPSHIFT experiment is applied, where the effective dipolar coupling is amplified by a REDOR-like pi-pulse train. However, a previously described constant-time variant of this experiment is not sensitive to the motion-induced T-2 effect, which precludes the observation of motions over a large range of rates ranging from hundreds of Hz to around a MHz. We present a DIPSHIFT implementation which amplifies the dipolar couplings and is still sensitive to T-2 effects. Spin dynamics simulations, analytical calculations and experiments demonstrate the sensitivity of the technique to molecular motions, and suggest the best experimental conditions to avoid imperfections. Furthermore, an in-depth theoretical analysis of the interplay of REDOR-like recoupling and proton decoupling based on Average-Hamiltonian Theory was performed, which allowed explaining the origin of many artifacts found in literature data. (C) 2012 Elsevier Inc. All rights reserved.

Structure determination of proteins and peptides in solution: simulation, chirality and NMR studies

The study of protein fold is a central problem in life science, leading in the last years to several attempts for improving our knowledge of the protein structures. In this thesis this challenging problem is tackled by means of molecular dynamics, chirality and NMR studies. In the last decades, many algorithms were designed for the protein secondary structure assignment, which reveals the local protein shape adopted by segments of amino acids. In this regard, the use of local chirality for the protein secondary structure assignment was demonstreted, trying to correlate as well the propensity of a given amino acid for a particular secondary structure. The protein fold can be studied also by Nuclear Magnetic Resonance (NMR) investigations, finding the average structure adopted from a protein. In this context, the effect of Residual Dipolar Couplings (RDCs) in the structure refinement was shown, revealing a strong improvement of structure resolution. A wide extent of this thesis is devoted to the study of avian prion protein. Prion protein is the main responsible of a vast class of neurodegenerative diseases, known as Bovine Spongiform Encephalopathy (BSE), present in mammals, but not in avian species and it is caused from the conversion of cellular prion protein to the pathogenic misfolded isoform, accumulating in the brain in form of amiloyd plaques. In particular, the N-terminal region, namely the initial part of the protein, is quite different between mammal and avian species but both of them contain multimeric sequences called Repeats, octameric in mammals and hexameric in avians. However, such repeat regions show differences in the contained amino acids, in particular only avian hexarepeats contain tyrosine residues. The chirality analysis of avian prion protein configurations obtained from molecular dynamics reveals a high stiffness of the avian protein, which tends to preserve its regular secondary structure. This is due to the presence of prolines, histidines and especially tyrosines, which form a hydrogen bond network in the hexarepeat region, only possible in the avian protein, and thus probably hampering the aggregation.

Heteronuclear recoupling methods in solid-state NMR

This thesis is focused on the development of heteronuclear correlation methods in solid-state NMR spectroscopy, where the spatial dependence of the dipolar coupling is exploited to obtain structural and dynamical information in solids. Quantitative results on dipolar coupling constants are extracted by means of spinning sideband analysis in the indirect dimension of the two-dimensional experiments. The principles of sideband analysis were established and are currently widely used in the group of Prof. Spiess for the special case of homonuclear 1H double-quantum spectroscopy. The generalization of these principles to the heteronuclear case is presented, with special emphasis on naturally abundant 13C-1H systems. For proton spectroscopy in the solid state, line-narrowing is of particular importance, and is here achieved by very-fast sample rotation at the magic angle (MAS), with frequencies up to 35 kHz. Therefore, the heteronuclear dipolar couplings are suppressed and have to be recoupled in order to achieve an efficient excitation of the observed multiple-quantum modes. Heteronuclear recoupling is most straightforwardly accomplished by performing the known REDOR experiment, where pi-pulses are applied every half rotor period. This experiment was modified by the insertion of an additional spectroscopic dimension, such that heteronuclear multiple-quantum experiments can be carried out, which, as shown experimentally and theoretically, closely resemble homonuclear double-quantum experiments. Variants are presented which are well-suited for the recording of high-resolution 13C-1H shift correlation and spinning-sideband spectra, by means of which spatial proximities and quantitative dipolar coupling constants, respectively, of heteronuclear spin pairs can be determined. Spectral editing of 13C spectra is shown to be feasible with these techniques. Moreover, order phenomena and dynamics in columnar mesophases with 13C in natural abundance were investigated. Two further modifications of the REDOR concept allow the correlation of 13C with quadrupolar nuclei, such as 2H. The spectroscopic handling of these nuclei is challenging in that they cover large frequency ranges, and with the new experiments it is shown how the excitation problem can be tackled or circumvented altogether, respectively. As an example, one of the techniques is used for the identification of a yet unknown motional process of the H-bonded protons in the crystalline parts of poly(vinyl alcohol).

A study of nitronyl and imino nitroxide radicals attached to heterocyclic cores

Stabile organische Radikale mit zusätzlichen Funktionalitäten wie Donor/Akzepotor Eigenschaften und Ligandeneignung für Übergangsmetallkomplexierung repräsentieren eine synthetische Herausforderung beim Streben nach der Konstruktion hochdimensionaler heterospin Strukturen. In diesem Hinblick wurden acht neue Hochspinbiradikal-Moleküle zusammen mit ihren Monoradikal- Pendants in dieser Arbeit hergestellt. Die Wahl der Liganden als organische Distanzhalter der Radikaleinheiten wurde auf stickstoffhaltige Heterozyklen (Pyridin und Pyrazol) gelenkt. Diese wurden weiterhin mit den stabilen Spinträgern Nitronylnitroxid- (NN) und Iminonitroxidfragmenten (IN) dekoriert. Ihre Synthese beinhaltete mehrstufige Umsetzungen (Brominierung, Iodierung, N- und Carbaldehyd Schutzgruppen, Stille-Kupplung, Grignard Reaktion, etc.) um die Mono- und Dicarbaldehyd-heterocyclenderivate als Schlüsselvorläufer der Radikaleinheiten zu gewinnen. Die Carbaldehyd-Zwischenstufen wurden Kondensationsreaktionen mit 2,3-Dimethyl-2,3-bis(hydroxylamino)-butan unterworfen (üblicherweise in Dioxan unter Argon für ~ 7 Tage), gefolgt von der Oxidation der Bis-hydroxylimidazolidin-Vorläufer unter Phasentransferkatalyse (NaIO4/H2O). Die Radikalmoleküle wurden mit verschiedenen spektroskopischen Methoden untersucht (FT/IR, UV/Vis/ EPR etc.) und ihre Einkristalle mit Röntgenstrahlbeugung gemessen. Die UV/VIS- Lösungsspektren zeigten in einem breiten Bereich verschiedener Lösungsmittelpolaritäten keine spezifische Wechselwirkung zwischen Lösungsmittel und Radikaleinheit, während ihre Stabilitäten in protischen Lösunsgmitteln wie MeOH stark abnahmen. Als Pulver konnten sie jedoch im Kühlschrank an der Luft für eine Jahr gelagert werden, ohne sich zu zersetzen. Die spektroskopischen Fingerabdrücke der Radikale wurden eindeutig identifiziert and erschienen stark abhängig vom Typ des pi-Ringsystems an das die Spinträger gekoppelt wurden. Basierend auf diesen Informationen wurde ein schnelles Protokoll etabliert, das eine direkte Zuordnung der Art der Radikale und ihrer Anzahl ermöglicht, sowie ihre Reinheit und Verunreinigungen zu definieren. In Lösung bestätigte die Analyse der EPR Spektren der Biradikale die starke Austauschwechselwirkung J zwischen den Radikalfragmenten über die Kopplungseinheiten (J >> an, an ist die Stickstoffhyperfeinkopplungskonstante). Dies wurde weiter unterstützt durch die Beobachtungen in gefrorener Lösung über die Nullfeldaufspaltungen und verbotenen Halbfeldübergänge (Δms = 2). Die Temperaturabhängigkeiten der Δms = 2 - EPR Signale wurden bis herunter auf 4 K gemessen und das exakte Vorzeichen und die Größe von J ermittelt. Diese Arbeit unterstreicht die Möglichkeit über synthetische Chemie eine Feineinstellung der „through bond“ Austauschwechselwirkung zwischen verwandten pi- und sigma- konjugierten Heterozyklen zu erreichen, in denen der S = 1 Grundzustand angenommen wird. Zusätzlich zeigten diese Resultate, dass die Übertragung der Spinpolarisation durch verschiedene Koppler sehr effektiv war.

Struktur-Dynamik-Beziehungen in amorphen festen Lösungen

Feste Lösungen homogen dispergierter Wirkstoffmoleküle in amorphen Polymermatrizen sind wichtige Materialien in vielen pharmazeutischen Anwendungen, bei denen eine kontrollierte Abgabe wasserunlöslicher Wirkstoffe in wässrige Systeme eine Rolle spielt. Die intermolekulare Bindungs-stärke zwischen Polymer- und Wirkstoffmolekülgruppen bestimmt die Stabilität der festen Lösung und steuert somit die biologische Aktivität der Wirkstoffmoleküle. In festen Lösungen, die aus acryl-säurehaltigen Copolymeren (Protonendonoren) und basischen Wirkstoffmolekülen (Protonenakzepto-ren) hergestellt werden, sind intermolekulare Wasserstoffbrücken zwischen den Systemkomponenten Triebkraft für die Bildung einer stabilen homogenen Dispersion und für die Entstehung struktureller Merkmale zwischen den Molekülgruppen der Systemkomponenten. Zudem ist die Bindungsstärke der Wasserstoffbrücken im Hinblick auf die kontrollierte Abgabe der Wirkstoffe von Bedeutung. Da dynamische chemische Gleichgewichte bei der Bildung der Wasserstoffbrücken eine wichtige Rolle spielen müssen neben strukturellen Parametern auch dynamische Faktoren beleuchtet werden. Ziel dieser Arbeit ist neben der Ermittlung von intermolekularen Bindungsstärken vor allem die Identifika-tion struktureller Verhältnisse zwischen den Systemkomponenten auf molekularer Ebene. Die Be-stimmung der Abhängigkeit dieser Parameter von der Struktur der verwendeten Polymere und einer Vielzahl weiterer Einflüsse wie z.B. Feuchtigkeit, Lagerdauer oder Wirkstoffkonzentration soll ein kontrolliertes Design fester Lösungen mit definierten anwendungsspezifischen Eigenschaften ermögli-chen. Temperaturabhängige 1H-Festkörper-MAS-NMR (Magic Angle Spinning Nuclear Magnetic Resonance) Experimente an festen Lösungen mit unterschiedlichen Copolymer-Zusammensetzungen weisen die Existenz dynamischer chemischer Gleichgewichte in den komplexen Wasserstoffbrücken-netzwerken nach. Veränderungen in der chemischen Verschiebung und in der Linienform der Reso-nanzlinien acider Protonen erlauben einen tiefen Einblick in die Architektur dieser Netzwerke und legen die Bindungsverhältnisse unter Berücksichtigung der Polymerchemie und der Mobilität der Systemkomponenten dar, wobei die Befunde mithilfe quantenchemischer Rechnungen untermauert werden können. Die Gegenwart acider Protonen ermöglicht einen einfachen 1H-2H-Austausch, wor-aufhin mithilfe rotorsynchronisierter temperaturabhängiger 2H-MAS-NMR Experimente die Wasser-stoffbrückenbindungsstärke bestimmt werden kann. Mit 1H-1H-Korrelationsexperimenten (Doppelquantenspektroskopie) stehen Methoden für die Bestimmung homonuklearer dipolarer 1H-1H-Kopplungen zur Verfügung, die strukturelle Aussagen aufgrund von bevorzugten räumlichen Kontak-ten bestimmter Molekülgruppen ermöglichen. Weiterhin können diese Experimente verwendet werden, um Wasserstoffbrücken zwischen Polymergruppen von Polymer-Wirkstoff-Wasserstoffbrücken zu unterscheiden, wodurch eine quantitative Beschreibung des Bindungsnetzwerks und der Konkurrenz-prozesse zwischen den einzelnen wasserstoffverbrückten Spezies ermöglicht wird. Eine Kristallisation der Wirkstoffmoleküle ist in vielen Anwendungen unerwünscht, da sie die biologische Verfügbarkeit des Wirkstoffs reduzieren. Mit 1H-Festkörper-MAS-NMR Experimenten können kristalline von amorph dispergierten Wirkstoffmolekülen unterschieden werden, wodurch eine Quantifizierung der Destabilisierungsprozesse ermöglicht wird, die durch Exposition der festen Lösungen mit Wasserdampf ausgelöst werden können. Die Zeit- und Konzentrationsabhängigkeit der Wasseraufnahme kann mit NMR-Experimenten verfolgt werden, wobei unterschiedlich mobile Was-serspezies an unterschiedlichen Bindungsorten identifiziert werden können, was zum molekularen Verständnis der Destabilisierungsprozesse beiträgt. Zusätzlich wird die Mobilität der Wirkstoffmole-küle bestimmt, die sich – wie auch die Wirkstoffkonzentration - als wichtige Größe in der Beschrei-bung der Destabilisierung erweist. Aufbauend auf den Beobachtungen wird ein Zusammenhang zwischen der Copolymerzusammensetzung und einer kritischen Wirkstoffkonzentration hergestellt, der für die Anwendungen amorpher fester Lösungen in biologischen Systemen von großer Bedeutung ist.

New methods of structure refinement for macromolecular structure determination by NMR

Recent advances in multidimensional NMR methodology have permitted solution structures of proteins in excess of 250 residues to be solved. In this paper, we discuss several methods of structure refinement that promise to increase the accuracy of macromolecular structures determined by NMR. These methods include the use of a conformational database potential and direct refinement against three-bond coupling constants, secondary 13C shifts, 1H shifts, T1/T2 ratios, and residual dipolar couplings. The latter two measurements provide long range restraints that are not accessible by other solution NMR parameters.