Molecular Structure Refinement Based on Residual Dipolar Couplings: A Comparison of the Molecular Rotational-Sampling Method with the Alignment-Tensor Approach

In NMR experiments, residual dipolar couplings (RDCs) in a molecule can be measured by averaging the dipolar couplings (DCs) over the rotational motion of a molecule in an environment that induces a slight anisotropic orientation distribution of the molecule. Since the shape of the anisotropic distribution cannot be measured, it is standard practice to use a particular orientation distribution of the molecule with respect to the magnetic field, in the form of a so-called alignment tensor (AT), to calculate RDC-values for the molecule. Since the same alignment tensor is commonly used to calculate the different RDCs of a molecule, this approach rests on the assumption that the rotational motion of the molecule is decoupled from its internal motions and that the molecule is rigid. The validity of these two assumptions is investigated for a small, simple molecule, using a relatively rigid atomic interaction function or force field and a more flexible one. By simulating the molecule using an orientation-biasing force an anisotropic rotational distribution can be generated, for which RDCs can be obtained. Using these RDCs as target RDCs when applying one of the two approaches of structure refinement based on RDCs, it can be investigated how well the target RDCs are approximated in the RDC restraining and whether the corresponding nonuniform orientation distribution is reproduced. For the relatively rigid version of the molecule, the AT approach reproduces the target RDC-values, although the nonuniform orientation distribution of the angle theta(ab,H) between the vector r(ab)(->) connecting two atoms a and b in the molecule and the vector representing the direction of the magnetic field H--> as generated in the orientation-biasing simulation cannot be reproduced in the AT RDC-restraining simulation. For the relatively flexible version of the molecule, the AT approach fails to reproduce both the target RDC values and the nonuniform orientation distribution. For biomolecules with flexible parts, the application of the AT approach is thus not recommended. Instead, a method based on sampling of the rotational and internal degrees of freedom of the molecule should be applied in molecular structure determination or refinement based on measured RDCs.