Stochastic modeling of flexible biomolecules applied to NMR relaxation. I. Internal dynamics of cyclodextrins: γ-cyclodextrin as a case study

Stochastic Modelling of 13C NMR Spin Relaxation Experiments in Oligosaccharides

A framework for the stochastic description of relaxation processes in flexible macromolecules including dissipative effects has been recently introduced, starting from an atomistic view, describing the joint relaxation of internal coordinates and global degrees of freedom, and depending on parameters recoverable from classic force fields (energetics) and medium modelling at the continuum level (friction tensors). The new approach provides a rational context for the interpretation of magnetic resonance relaxation experiments. In its simplest formulation, the semi-flexible Brownian (SFB) model has been until now shown to reproduce correctly correlation functions and spectral densities related to orientational properties obtained by direct molecular dynamics simulations of peptides. Here, for the first time, we applied directly the SFB approach to the practical evaluation of high-quality ¹³C nuclear magnetic resonance relaxation parameters, T1 and T2, and the heteronuclear NOE of several oligosaccharides, which were previously interpreted on the basis of refined ad hoc modelling. The calculated NMR relaxation parameters were in agreement with the experimental data, showing that this general approach can be applied to diverse classes of molecular systems, with the minimal usage of adjustable parameters.

Stochastic modeling of macromolecules in solution. I. Relaxation processes

A framework for the stochastic description of relaxation processes in flexible macromolecules, including dissipative effects, is introduced from an atomistic point of view. Projection-operator techniques are employed to obtain multidimensional Fokker-Planck operators governing the relaxation of internal coordinates and global degrees of freedom and depending upon parameters fully recoverable from classic force fields (energetics) and continuum models (friction tensors). A hierarchy of approaches of different complexity is proposed in this unified context, aimed primarily at the interpretation of magnetic resonance relaxation experiments. In particular, a model based on a harmonic internal Hamiltonian is discussed as a test case.

Flexibility at a glycosidic linkage revealed by molecular dynamics, stochastic modeling, and (13)C NMR spin relaxation: conformational preferences of α-L-Rhap-α-(1 → 2)-α-L-Rhap-OMe in water and dimethyl sulfoxide solutions

The monosaccharide L-rhamnose is common in bacterial polysaccharides and the disaccharide α-L-Rhap-α-(1 → 2)-α-L-Rhap-OMe represents a structural model for a part of Shigella flexneri O-antigen polysaccharides. Utilization of [1'-(13)C]-site-specific labeling in the anomeric position at the glycosidic linkage between the two sugar residues facilitated the determination of transglycosidic NMR (3)JCH and (3)JCC coupling constants. Based on these spin-spin couplings the major state and the conformational distribution could be determined with respect to the ψ torsion angle, which changed between water and dimethyl sulfoxide (DMSO) as solvents, a finding mirrored by molecular dynamics (MD) simulations with explicit solvent molecules. The (13)C NMR spin relaxation parameters T1, T2, and heteronuclear NOE of the probe were measured for the disaccharide in DMSO-d6 at two magnetic field strengths, with standard deviations ≤1%. The combination of MD simulation and a stochastic description based on the diffusive chain model resulted in excellent agreement between calculated and experimentally observed (13)C relaxation parameters, with an average error of <2%. The coupling between the global reorientation of the molecule and the local motion of the spin probe is deemed essential if reproduction of NMR relaxation parameters should succeed, since decoupling of the two modes of motion results in significantly worse agreement. Calculation of (13)C relaxation parameters based on the correlation functions obtained directly from the MD simulation of the solute molecule in DMSO as solvent showed satisfactory agreement with errors on the order of 10% or less.

Toward a General Yet Effective Computational Approach for Diffusive Problems: Variable Diffusion Tensor and DVR Solution of the Smoluchowski Equation along a General One-Dimensional Coordinate

A generalization to arbitrary large amplitude motions of a recent approach to the evaluation of diffusion tensors [ J. Comput. Chem. , 2009 , 30 , 2 - 13 ] is presented and implemented in a widely available package for electronic structure computations. A fully black-box tool is obtained, which, starting from the generation of geometric structures along different kinds of paths, proceeds toward the evaluation of an effective diffusion tensor and to the solution of one-dimensional Smoluchowski equations by a robust numerical approach rooted in the discrete variable representation (DVR). Application to a number of case studies shows that the results issuing from our approach are identical to those delivered by previous software (in particular DiTe) for rigid scans along a dihedral angle, but can be improved by employing relaxed scans (i.e., constrained geometry optimizations) or even more general large amplitude paths. The theoretical and numerical background is robust and general enough to allow quite straightforward extensions in several directions (e.g., inclusion of reactive paths, solution of Fokker-Planck or stochastic Liouville equations, multidimensional problems, free-energy rather than electronic-energy driven processes).

Reducing acquisition times in multidimensional NMR with a time-optimized Fourier encoding algorithm

Speeding up the acquisition of multidimensional nuclear magnetic resonance (NMR) spectra is an important topic in contemporary NMR, with central roles in high-throughput investigations and analyses of marginally stable samples. A variety of fast NMR techniques have been developed, including methods based on non-uniform sampling and Hadamard encoding, that overcome the long sampling times inherent to schemes based on fast-Fourier-transform (FFT) methods. Here, we explore the potential of an alternative fast acquisition method that leverages a priori knowledge, to tailor polychromatic pulses and customized time delays for an efficient Fourier encoding of the indirect domain of an NMR experiment. By porting the encoding of the indirect-domain to the excitation process, this strategy avoids potential artifacts associated with non-uniform sampling schemes and uses a minimum number of scans equal to the number of resonances present in the indirect dimension. An added convenience is afforded by the fact that a usual 2D FFT can be used to process the generated data. Acquisitions of 2D heteronuclear correlation NMR spectra on quinine and on the anti-inflammatory drug isobutyl propionic phenolic acid illustrate the new method's performance. This method can be readily automated to deal with complex samples such as those occurring in metabolomics, in in-cell as well as in in vivo NMR applications, where speed and temporal stability are often primary concerns.

Stochastic modeling of flexible biomolecules applied to NMR relaxation. 2. Interpretation of complex dynamics in linear oligosaccharides

A computational stochastic approach is applied to the description of flexible molecules. By combining (i) molecular dynamics simulations, (ii) hydrodynamics approaches, and (iii) a multidimensional diffusive description for internal and global dynamics, it is possible to build an efficient integrated approach to the interpretation of relaxation processes in flexible systems. In particular, the model is applied to the interpretation of nuclear magnetic relaxation measurements of linear oligosaccharides, namely a mannose-containing trisaccharide and the pentasaccharide LNF-1. Experimental data are reproduced with sufficient accuracy without free model parameters.