A Neural Network Approach to the Rapid Computation of Rotational
Correlation Times from Slow Motional ESR Spectra
Gary V. Martinez and Glenn L. Millhauser
1
Department of Chemistry and Biochemistry, University of California, Santa Cruz, Santa Cruz, California 95064
Received August 18, 1997; revised April 7, 1998
We explore the use of feed forward artificial neural networks for
determining rotational correlation times from slow motional ni-
troxide electron spin resonance spectra. This approach is rapid
and potentially eliminates the need for traditional iterative fitting
procedures. Two networks are examined:the radial basis network
and the multilayerperceptron. Although the radial basis network
trains rapidly and performs well on simulated spectra, it is less
satisfactory when applied to experimental spectra. In contrast, the
multilayer perceptron trains slowly but is excellent at extracting
correlation times from experimental spectra. In addition, the mul-
tilayerperceptron operates well in the presence of noise as long as
the signal-to-noise ratio is greater than approximately 200/1.
These findings suggest neural networks offer a promising ap-
proach for rapidly extracting correlation times without the need
for iterative simulations. © 1998 Academic Press
Key Words: neural networks; ESR; slow motional spectra; radial
basis network; multilayer perceptron; rotational correlation time.
INTRODUCTION
Nitroxide spin labels, in conjunction with electron spin
resonance (ESR), serve as probes for exploring dynamics in
nucleic acids, peptides and proteins (1–3). The ESR lineshape
is often used to determine the rotational correlation time (
R
) of
the nitroxide, thereby revealing motion at the label site. When
R
is less than approximately 1 ns, ESR spectra are character-
ized by three motionally narrowed hyperfine lines, and
straightforward lineshape measurement gives accurate values
for the correlation time. However, when
R
is between 1 ns and
100 ns—the so-called slow motional regime—lineshape anal-
ysis is substantially more complicated (4). The correlation
times of most macromolecules fall within this regime and
consequently enormous effort has been directed toward ex-
tracting dynamic information from slow motional spectra.
Leading efforts in this field have come from Freed and
co-workers (5) (see also Chapter 3 in Ref. 1). Throughout the
1970s and 1980s they developed and refined slow motional
simulation techniques based on the stochastic Liouville equa-
tion. Despite their great success, determining dynamic param-
eters from slow motional spectra remains a challenge. The
essential problem lies with the iterative approach one must
employ when simulating spectra. Given an experimental spec-
trum, one must first guess at a value for the correlation time (as
well as other parameters, including motional anisotropy, local
ordering, and sample heterogeneity) and then perform a sim-
ulation. The result is compared to the experimental spectrum.
If the agreement is not satisfactory, the input values are ad-
justed and another simulation is performed. This process is
repeated until good overlap between the experimental and
simulated spectra is achieved.
Recently, Budil et al. developed a nonlinear least squares
approach to automate the procedure of fitting slow motional
spectra (6). They applied a “model trust region” modification
of the Levenberg–Marquardt algorithm and showed that com-
plicated spectra could be successfully simulated using a com-
puter workstation. Nevertheless, their method still relies on
iteration.
The difference in effort required for a noniterative approach,
such as that used for analyzing fast motional spectra, and the
iterative approach required for slow motional spectra is strik-
ing. It would be quite helpful if a method could be developed
that used a noniterative approach for analyzing slow motional
spectra. Toward this goal we explore the use of artificial neural
networks (7). Neural networks have emerged as remarkable
tools for pattern recognition in scientific applications. They
have been applied with good success to spectroscopic problems
in nuclear magnetic resonance (8 –10), circular dichroism (11–
13), and infrared spectroscopy (14).
Artificial neural networks were inspired by research into the
interplay that takes place among networks of real biological
neurons. In an artificial neural network, the synaptic connec-
tions between neurons are represented by numerical weights,
which measure the strength of a connection, and a transfer
function that emulates the firing of the neuron. Training a
network involves establishing a set of numerical weights that
successfully connect a training input with a desired output.
Once trained, an artificial neural network can be an effective
tool for recognizing and extracting key features from previ-
ously unseen input.
In spectroscopic applications, neural networks have been
1
To whom correspondence should be addressed. Fax: (408) 459-2935.
E-mail: glennm@hydrogen.ucsc.edu
JOURNAL OF MAGNETIC RESONANCE 134, 124 –130 (1998)
ARTICLE NO. MN981496
124
1090-7807/98 $25.00
Copyright © 1998 by Academic Press
All rights of reproduction in any form reserved.