A robust method for estimating cross-relaxation rates from simultaneous fits to build-up and decay curves

This paper introduces a novel computational method for estimating relaxation rates among pairs of spin orders. This method simultaneously estimates all the auto- and cross-relaxation rates from the same measurements, and avoids the ill-conditioning problems associated with multiexponential fits. The method models the relaxation dynamics by a system of linear differential equations, and assumes that measurements of the spin orders have been made at an equally spaced sequence of time points. It computes a nonlinear least-squares fit of the exponential of the rate matrix at the shortest time point to these measurements. Preliminary estimates of the exponential matrix and initial spin orders from which to start the computations are obtained by solving simpler linear-least-squares problems. The performance of the method on simulated 2 x 2 test problems indicates that when measurements at eight or more equally spaced times spanning the maximum and inflection points of the build-up curves are available, the relative errors in the rates are usually less than the relative errors in the measurements. The method is further demonstrated by applying it to the problem of determining the cross correlation-induced cross-relaxation rates between the in-phase and antiphase coherence of the amide groups in the 15N-labeled protein oxidized flavodoxin. Finally, the possibility of extending the method to other kinds of relaxation measurements and larger spin systems is discussed.