Simulations of NMR pulse sequences during equilibrium and non-equilibrium chemical exchange

The McConnell equations combine the differential equations for a simple two-state chemical exchange process with the Bloch differential equations for a classical description of the behavior of nuclear spins in a magnetic field. This equation system provides a useful starting point for the analysis of slow, intermediate and fast chemical exchange studied using a variety of NMR experiments. The McConnell equations are in the mathematical form of an inhomogeneous system of first-order differential equations. Here we rewrite the McConnell equations in a homogeneous form in order to facilitate fast and simple numerical calculation of the solution to the equation system. The McConnell equations can only treat equilibrium chemical exchange. We therefore also present a homogeneous equation system that can handle both equilibrium and non-equilibrium chemical processes correctly, as long as the kinetics is of first-order. Finally, the same method of rewriting the inhomogeneous form of the McConnell equations into a homogeneous form is applied to a quantum mechanical treatment of a spin system in chemical exchange. In order to illustrate the homogeneous McConnell equations, we have simulated pulse sequences useful for measuring exchange rates in slow, intermediate and fast chemical exchange processes. A stopped-flow NMR experiment was simulated using the equations for non-equilibrium chemical exchange. The quantum mechanical treatment was tested by the simulation of a sensitivity enhanced 15N-HSQC with pulsed field gradients during slow chemical exchange and by the simulation of the transfer efficiency of a two-dimensional heteronuclear cross-polarization based experiment as a function of both chemical shift difference and exchange rate constants.