Discussion:
R2 from R1rho & R1
Justin Lecher
2014-03-26 11:19:16 UTC
Permalink
Hi all,

could someone please give me a short hint how I can calculate R2 values
from R1 and R1rho using relax?

Thanks,
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
Troels Emtekær Linnet
2014-03-26 11:47:39 UTC
Permalink
Hi Justin.

If one have R1, and R1rho, R2 (R1rho') is fitted in:
http://wiki.nmr-relax.com/DPL94
http://wiki.nmr-relax.com/TAP03
http://wiki.nmr-relax.com/TP02
http://wiki.nmr-relax.com/NS_R1rho_2-site

R1+(R2) can be made by:
http://wiki.nmr-relax.com/Tutorial_for_R1/R2_Relaxation_curve-fitting_analysis_on_varian_recorded_as_fid_interleaved
(A tutorial I made for our students in our lab. We have Varian at our
place).

Best
Troels
Post by Justin Lecher
Hi all,
could someone please give me a short hint how I can calculate R2 values
from R1 and R1rho using relax?
Thanks,
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
_______________________________________________
relax (http://www.nmr-relax.com)
This is the relax-users mailing list
To unsubscribe from this list, get a password
reminder, or change your subscription options,
visit the list information page at
https://mail.gna.org/listinfo/relax-users
Justin Lecher
2014-03-26 11:58:08 UTC
Permalink
Post by Troels Emtekær Linnet
Hi Justin.
http://wiki.nmr-relax.com/DPL94
http://wiki.nmr-relax.com/TAP03
http://wiki.nmr-relax.com/TP02
http://wiki.nmr-relax.com/NS_R1rho_2-site
http://wiki.nmr-relax.com/Tutorial_for_R1/R2_Relaxation_curve-fitting_analysis_on_varian_recorded_as_fid_interleaved
(A tutorial I made for our students in our lab. We have Varian at our
place).
Best
Troels
Hi Troels,

perfect, thanks! I found the wiki pages, but I wasn't sure which one is
the best to choose. Any general suggestion or does it strictly depend on
the type of model which applies best to my system?

thanks
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
Edward d'Auvergne
2014-03-26 12:39:31 UTC
Permalink
Hi,

The key here is what is meant by R2. There are many different
definitions. In relax, the parameter name R2 is defined as:

- In the standard relaxation equations, spin-spin relaxation
component. I.e. the part influenced by the spectral density function
J(w).

- Again in the standard relaxation equations, R2 is used for R2*,
where R2* = R2 + Rex. R2* is also defined differently if you look at
the exponential decay curves or if you looking at peak widths (for the
later there are additional factors broadening the peaks which add to
the R2* equation).

- For relaxation dispersion, R2 is currently defined as R1rho', i.e.
just the parts influenced by the spectral density function. Here
R1rho' is the on-resonance part of R1rho, excluding exchange.

R2 = R1rho' + Rex could also be defined, which is probably what you
are after. The value of R1rho' + Rex has been named many different
things by the field and there seems to be no consensus. I have
discussed this with Troels at
http://thread.gmane.org/gmane.science.nmr.relax.devel/5119/focus=5207.
This is a long thread with many discussions about implementing this
as an automatically calculated parameter - which in the end would be a
great feature.

So which R2 value are you after?

Regards,

Edward
Post by Justin Lecher
Post by Troels Emtekær Linnet
Hi Justin.
http://wiki.nmr-relax.com/DPL94
http://wiki.nmr-relax.com/TAP03
http://wiki.nmr-relax.com/TP02
http://wiki.nmr-relax.com/NS_R1rho_2-site
http://wiki.nmr-relax.com/Tutorial_for_R1/R2_Relaxation_curve-fitting_analysis_on_varian_recorded_as_fid_interleaved
(A tutorial I made for our students in our lab. We have Varian at our
place).
Best
Troels
Hi Troels,
perfect, thanks! I found the wiki pages, but I wasn't sure which one is
the best to choose. Any general suggestion or does it strictly depend on
the type of model which applies best to my system?
thanks
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
_______________________________________________
relax (http://www.nmr-relax.com)
This is the relax-users mailing list
To unsubscribe from this list, get a password
reminder, or change your subscription options,
visit the list information page at
https://mail.gna.org/listinfo/relax-users
Justin Lecher
2014-03-26 13:22:46 UTC
Permalink
Post by Edward d'Auvergne
Hi,
The key here is what is meant by R2. There are many different
- In the standard relaxation equations, spin-spin relaxation
component. I.e. the part influenced by the spectral density function
J(w).
- Again in the standard relaxation equations, R2 is used for R2*,
where R2* = R2 + Rex. R2* is also defined differently if you look at
the exponential decay curves or if you looking at peak widths (for the
later there are additional factors broadening the peaks which add to
the R2* equation).
- For relaxation dispersion, R2 is currently defined as R1rho', i.e.
just the parts influenced by the spectral density function. Here
R1rho' is the on-resonance part of R1rho, excluding exchange.
R2 = R1rho' + Rex could also be defined, which is probably what you
are after. The value of R1rho' + Rex has been named many different
things by the field and there seems to be no consensus. I have
discussed this with Troels at
http://thread.gmane.org/gmane.science.nmr.relax.devel/5119/focus=5207.
This is a long thread with many discussions about implementing this
as an automatically calculated parameter - which in the end would be a
great feature.
So which R2 value are you after?
Regards,
Edward
Hi Edward,

Principally I am looking for pure R2 without Rex contribution.

Currently we are using some home made scripts which basically are using
equation 20 in Palmer, Massi, 2006 (PMID 16683750).


So for the technical side, how do I proceed in relax? Do I need to fix
the decay for R1 and R1rho and then calculate R2 or just fit R1 and use
my raw intensities from the R1rho measurements for the R2 calculation?

Thanks,
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
Edward d'Auvergne
2014-03-26 13:42:38 UTC
Permalink
Hi,

Equation 20 in that paper (http://dx.doi.org/10.1021/cr0404287) is
actually the sum of the spin-spin relaxation and Rex! This is exactly
what Troels has been working with recently. See
http://wiki.nmr-relax.com/Matplotlib_DPL94_R1rho_R2eff. In the
relaxation dispersion analysis in relax we should really have two
parameters automatically generated for the user, R2 and Rex (defined
as R2 = R1rho' + Rex). These could then be used to generate
dispersion plots of R2 or Rex verses Omega_eff. This might take a
while though, in the mean time have a look at Troels' script on that
wiki page.

Regards,

Edward
Post by Justin Lecher
Post by Edward d'Auvergne
Hi,
The key here is what is meant by R2. There are many different
- In the standard relaxation equations, spin-spin relaxation
component. I.e. the part influenced by the spectral density function
J(w).
- Again in the standard relaxation equations, R2 is used for R2*,
where R2* = R2 + Rex. R2* is also defined differently if you look at
the exponential decay curves or if you looking at peak widths (for the
later there are additional factors broadening the peaks which add to
the R2* equation).
- For relaxation dispersion, R2 is currently defined as R1rho', i.e.
just the parts influenced by the spectral density function. Here
R1rho' is the on-resonance part of R1rho, excluding exchange.
R2 = R1rho' + Rex could also be defined, which is probably what you
are after. The value of R1rho' + Rex has been named many different
things by the field and there seems to be no consensus. I have
discussed this with Troels at
http://thread.gmane.org/gmane.science.nmr.relax.devel/5119/focus=5207.
This is a long thread with many discussions about implementing this
as an automatically calculated parameter - which in the end would be a
great feature.
So which R2 value are you after?
Regards,
Edward
Hi Edward,
Principally I am looking for pure R2 without Rex contribution.
Currently we are using some home made scripts which basically are using
equation 20 in Palmer, Massi, 2006 (PMID 16683750).
So for the technical side, how do I proceed in relax? Do I need to fix
the decay for R1 and R1rho and then calculate R2 or just fit R1 and use
my raw intensities from the R1rho measurements for the R2 calculation?
Thanks,
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
Troels Emtekær Linnet
2014-03-26 20:45:58 UTC
Permalink
Hi Edward.

I would advise not to save Rex as a value, since it is correlated with the
effective field.

Ah ja, that was the problem!
Theta and the effective field are correlated (obvious).

So one can only interpolate, by determining to ramp either spin lock field
or offset. I wonder how that looks as a 3D graph?

Best
Troels
Post by Edward d'Auvergne
Hi,
Equation 20 in that paper (http://dx.doi.org/10.1021/cr0404287) is
actually the sum of the spin-spin relaxation and Rex! This is exactly
what Troels has been working with recently. See
http://wiki.nmr-relax.com/Matplotlib_DPL94_R1rho_R2eff. In the
relaxation dispersion analysis in relax we should really have two
parameters automatically generated for the user, R2 and Rex (defined
as R2 = R1rho' + Rex). These could then be used to generate
dispersion plots of R2 or Rex verses Omega_eff. This might take a
while though, in the mean time have a look at Troels' script on that
wiki page.
Regards,
Edward
Post by Justin Lecher
Post by Edward d'Auvergne
Hi,
The key here is what is meant by R2. There are many different
- In the standard relaxation equations, spin-spin relaxation
component. I.e. the part influenced by the spectral density function
J(w).
- Again in the standard relaxation equations, R2 is used for R2*,
where R2* = R2 + Rex. R2* is also defined differently if you look at
the exponential decay curves or if you looking at peak widths (for the
later there are additional factors broadening the peaks which add to
the R2* equation).
- For relaxation dispersion, R2 is currently defined as R1rho', i.e.
just the parts influenced by the spectral density function. Here
R1rho' is the on-resonance part of R1rho, excluding exchange.
R2 = R1rho' + Rex could also be defined, which is probably what you
are after. The value of R1rho' + Rex has been named many different
things by the field and there seems to be no consensus. I have
discussed this with Troels at
http://thread.gmane.org/gmane.science.nmr.relax.devel/5119/focus=5207.
This is a long thread with many discussions about implementing this
as an automatically calculated parameter - which in the end would be a
great feature.
So which R2 value are you after?
Regards,
Edward
Hi Edward,
Principally I am looking for pure R2 without Rex contribution.
Currently we are using some home made scripts which basically are using
equation 20 in Palmer, Massi, 2006 (PMID 16683750).
So for the technical side, how do I proceed in relax? Do I need to fix
the decay for R1 and R1rho and then calculate R2 or just fit R1 and use
my raw intensities from the R1rho measurements for the R2 calculation?
Thanks,
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
Edward d'Auvergne
2014-03-27 07:51:34 UTC
Permalink
Hi Troels,

Storing such secondary dispersion parameters (Rex, R1rho'+Rex, etc.)
is not a good idea as they can be very difficult to keep up to date.
They are dependent on the base parameters of the model as well as the
experimental details (offsets, field strengths, etc.), so any change
in these values would require the secondary parameters to be updated -
but this approach is often a big source of future bugs. It is better
to generate them when asked by a user function, as 'theta' and 'w_eff'
currently are.

A 3D graph as a surface would be quite interesting! I have seen such
plots before, but I can't remember off the top of my head in which
papers that was. For this we would need to add some basic
infrastructure to relax for plotting with matplotlib. Grace does not
do 3D surface graphs. Adding this infrastructure is rather trivial as
the code for the Grace user functions can simply be copied and
modified. The hardest part is simply that I don't have experience
with matplotlib, so I don't know how to command it to produce high
quality plots. Maybe I'll post something on the relax-devel mailing
list about this.

Regards,

Edward
Post by Troels Emtekær Linnet
Hi Edward.
I would advise not to save Rex as a value, since it is correlated with the
effective field.
Ah ja, that was the problem!
Theta and the effective field are correlated (obvious).
So one can only interpolate, by determining to ramp either spin lock field
or offset. I wonder how that looks as a 3D graph?
Best
Troels
Post by Edward d'Auvergne
Hi,
Equation 20 in that paper (http://dx.doi.org/10.1021/cr0404287) is
actually the sum of the spin-spin relaxation and Rex! This is exactly
what Troels has been working with recently. See
http://wiki.nmr-relax.com/Matplotlib_DPL94_R1rho_R2eff. In the
relaxation dispersion analysis in relax we should really have two
parameters automatically generated for the user, R2 and Rex (defined
as R2 = R1rho' + Rex). These could then be used to generate
dispersion plots of R2 or Rex verses Omega_eff. This might take a
while though, in the mean time have a look at Troels' script on that
wiki page.
Regards,
Edward
Post by Justin Lecher
Post by Edward d'Auvergne
Hi,
The key here is what is meant by R2. There are many different
- In the standard relaxation equations, spin-spin relaxation
component. I.e. the part influenced by the spectral density function
J(w).
- Again in the standard relaxation equations, R2 is used for R2*,
where R2* = R2 + Rex. R2* is also defined differently if you look at
the exponential decay curves or if you looking at peak widths (for the
later there are additional factors broadening the peaks which add to
the R2* equation).
- For relaxation dispersion, R2 is currently defined as R1rho', i.e.
just the parts influenced by the spectral density function. Here
R1rho' is the on-resonance part of R1rho, excluding exchange.
R2 = R1rho' + Rex could also be defined, which is probably what you
are after. The value of R1rho' + Rex has been named many different
things by the field and there seems to be no consensus. I have
discussed this with Troels at
http://thread.gmane.org/gmane.science.nmr.relax.devel/5119/focus=5207.
This is a long thread with many discussions about implementing this
as an automatically calculated parameter - which in the end would be a
great feature.
So which R2 value are you after?
Regards,
Edward
Hi Edward,
Principally I am looking for pure R2 without Rex contribution.
Currently we are using some home made scripts which basically are using
equation 20 in Palmer, Massi, 2006 (PMID 16683750).
So for the technical side, how do I proceed in relax? Do I need to fix
the decay for R1 and R1rho and then calculate R2 or just fit R1 and use
my raw intensities from the R1rho measurements for the R2 calculation?
Thanks,
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
Edward d'Auvergne
2014-03-27 09:01:20 UTC
Permalink
Hi Troels,

Please see below. It might be best to shift these discussions to the
relax-devel mailing list as well.
Yep, I think some "porcelain functions" that will work on
relax base "plumbing" functions would be the best solution?
I'm not quite sure what you mean here, but it sounds good ;) In relax
there are a lot of API interfaces between the different parts of the
code. By keeping a clean separation between the different parts of
relax, adding interpolation to one part requires zero changes to other
parts of the program. I guess this is along the same lines.
That could return info to screen, and return values in dictionarys, where
the dictionary keys are also stored.
Actually, an extension of "calc_rotating_frame_params" in
"specific_analyses.relax_disp.disp_data"
which should take interpolated offset data as a possibility.
Extending calc_rotating_frame_params() is not necessary and is more of
a complication. A single interpolation function where you pass in the
offset array with arguments for how much interpolation and how much
extension to add, and then that function returns the interpolated
array. I.e. exactly as is performed in the
relax_disp.plot_disp_curves user function right now. This is after
assembling the data, so we have jumped through one API already and the
dictionaries and keys no longer exist at this stage. This makes life
a lot easier :)
For 3D, I would prefer that we use gnuplot instead ?
In relax there is currently Grace and OpenDX plotting abilities.
Adding matplotlib and gnuplot would be complementary. For matplotlib
we would have to detect if this Python package is installed - which is
rather trivial. But for gnuplot, this can be available at all times.
The most difficult part for gnuplot graphing support would be to find
a relax developer who knows the gnuplot commands for creating graphs!

This can also be saved for the relax-devel mailing list, but a plan
for the future would be to collect all these user functions into one
'plotting' user function class, for example:

plotting.dx_execute
plotting.dx_map
plotting.gnuplot_view
plotting.gnuplot_write
plotting.grace_view
plotting.grace_write
plotting.matplotlib_view
plotting.matplotlib_write

This is simply to make it easier for the user to find and choose.
I really really start liking the idea of the grace files.
Here you in "one" file carry around both the data and definitions to the
plot.
It is rather old software, but it is incredibly powerful. And it
produces the highest quality, professional-looking 2D plots for
papers.
For matplotlib, you would need to acquire that the user has relax and
matplotlib installed at the same time.
That is not so good, and would be a problem for the compiled sources.
I am writing something at the moment for the relax-devel mailing list.
But this is actually trivial. If matplotlib is not present, we just
give a RelaxWarning to the user when they try to create a matplotlib
graph (this includes in the auto-analyses which are simply very large
relax scripts).
Unless you export python files, with the data inside, and then the plot
commands to matplotlib.
This would need to be discussed on the relax-devel mailing list. As I
don't know matplotlib, I don't know the best way to do this. Do you
create Python files with data so the user can change the plot, or do
everything in relax and output to SVG or some other format, as is
currently done in NESSY? We can save this discussion for the
relax-devel mailing list.
Hm hm.
What is best?
Expect the user to know gnuplot , and have that installed?
Or that the user have matplotlib ?
The best here would be to give everything to the user. Then you will
hear less from the users ;) But this can be a gradual process. We
could even one day add a new chapter to the relax manual about data
plotting, just to help the user understand what the best tools are and
how to create and visualise these plots.
I know that matplotlib on our RHEL system is only 0.99, and is not moving.
Gnuplot is old... And proven to work.
Reach for the future or backwards-compability?
relax is flexible enough to do everything! For the record:

- gnuplot is from 1986 but is still developed.
- Grace is from 1991 but is also still developed (just).
- OpenDX is from 1991 but is abandoned. Though it is insanely
powerful for multi-dimensional data plotting!
- matplotlib is from 2002 and is still active.

For the future though, you can never know when something will be
abandoned or when abandoned open source software will come back in a
more powerful incarnation. So instead of choosing what would be best
and forcing that on a user, a relax developer who prefers one plotting
software over another can simply add support for that software. And a
relax user with some scripting knowledge can easily become a relax
developer and add it themselves :)

Regards,

Edward
Troels Emtekær Linnet
2014-03-27 08:25:19 UTC
Permalink
Hi Edward.

Yep, I think some "porcelain functions" that will work on
relax base "plumbing" functions would be the best solution?

That could return info to screen, and return values in dictionarys, where
the dictionary keys are also stored.
Actually, an extension of "calc_rotating_frame_params"
in "specific_analyses.relax_disp.disp_data"
which should take interpolated offset data as a possibility.

For 3D, I would prefer that we use gnuplot instead ?

I really really start liking the idea of the grace files.
Here you in "one" file carry around both the data and definitions to the
plot.

For matplotlib, you would need to acquire that the user has relax and
matplotlib installed at the same time.
That is not so good, and would be a problem for the compiled sources.

Unless you export python files, with the data inside, and then the plot
commands to matplotlib.

Hm hm.

What is best?

Expect the user to know gnuplot , and have that installed?
Or that the user have matplotlib ?

I know that matplotlib on our RHEL system is only 0.99, and is not moving.
Gnuplot is old... And proven to work.

Reach for the future or backwards-compability?

Best
Troels



Troels Emtekær Linnet
Post by Justin Lecher
Hi Troels,
Storing such secondary dispersion parameters (Rex, R1rho'+Rex, etc.)
is not a good idea as they can be very difficult to keep up to date.
They are dependent on the base parameters of the model as well as the
experimental details (offsets, field strengths, etc.), so any change
in these values would require the secondary parameters to be updated -
but this approach is often a big source of future bugs. It is better
to generate them when asked by a user function, as 'theta' and 'w_eff'
currently are.
A 3D graph as a surface would be quite interesting! I have seen such
plots before, but I can't remember off the top of my head in which
papers that was. For this we would need to add some basic
infrastructure to relax for plotting with matplotlib. Grace does not
do 3D surface graphs. Adding this infrastructure is rather trivial as
the code for the Grace user functions can simply be copied and
modified. The hardest part is simply that I don't have experience
with matplotlib, so I don't know how to command it to produce high
quality plots. Maybe I'll post something on the relax-devel mailing
list about this.
Regards,
Edward
Post by Troels Emtekær Linnet
Hi Edward.
I would advise not to save Rex as a value, since it is correlated with
the
Post by Troels Emtekær Linnet
effective field.
Ah ja, that was the problem!
Theta and the effective field are correlated (obvious).
So one can only interpolate, by determining to ramp either spin lock
field
Post by Troels Emtekær Linnet
or offset. I wonder how that looks as a 3D graph?
Best
Troels
Post by Edward d'Auvergne
Hi,
Equation 20 in that paper (http://dx.doi.org/10.1021/cr0404287) is
actually the sum of the spin-spin relaxation and Rex! This is exactly
what Troels has been working with recently. See
http://wiki.nmr-relax.com/Matplotlib_DPL94_R1rho_R2eff. In the
relaxation dispersion analysis in relax we should really have two
parameters automatically generated for the user, R2 and Rex (defined
as R2 = R1rho' + Rex). These could then be used to generate
dispersion plots of R2 or Rex verses Omega_eff. This might take a
while though, in the mean time have a look at Troels' script on that
wiki page.
Regards,
Edward
Post by Justin Lecher
Post by Edward d'Auvergne
Hi,
The key here is what is meant by R2. There are many different
- In the standard relaxation equations, spin-spin relaxation
component. I.e. the part influenced by the spectral density function
J(w).
- Again in the standard relaxation equations, R2 is used for R2*,
where R2* = R2 + Rex. R2* is also defined differently if you look at
the exponential decay curves or if you looking at peak widths (for
the
Post by Troels Emtekær Linnet
Post by Edward d'Auvergne
Post by Justin Lecher
Post by Edward d'Auvergne
later there are additional factors broadening the peaks which add to
the R2* equation).
- For relaxation dispersion, R2 is currently defined as R1rho', i.e.
just the parts influenced by the spectral density function. Here
R1rho' is the on-resonance part of R1rho, excluding exchange.
R2 = R1rho' + Rex could also be defined, which is probably what you
are after. The value of R1rho' + Rex has been named many different
things by the field and there seems to be no consensus. I have
discussed this with Troels at
http://thread.gmane.org/gmane.science.nmr.relax.devel/5119/focus=5207.
Post by Troels Emtekær Linnet
Post by Edward d'Auvergne
Post by Justin Lecher
Post by Edward d'Auvergne
This is a long thread with many discussions about implementing this
as an automatically calculated parameter - which in the end would be
a
Post by Troels Emtekær Linnet
Post by Edward d'Auvergne
Post by Justin Lecher
Post by Edward d'Auvergne
great feature.
So which R2 value are you after?
Regards,
Edward
Hi Edward,
Principally I am looking for pure R2 without Rex contribution.
Currently we are using some home made scripts which basically are
using
Post by Troels Emtekær Linnet
Post by Edward d'Auvergne
Post by Justin Lecher
equation 20 in Palmer, Massi, 2006 (PMID 16683750).
So for the technical side, how do I proceed in relax? Do I need to fix
the decay for R1 and R1rho and then calculate R2 or just fit R1 and
use
Post by Troels Emtekær Linnet
Post by Edward d'Auvergne
Post by Justin Lecher
my raw intensities from the R1rho measurements for the R2 calculation?
Thanks,
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
Troels Emtekær Linnet
2014-03-26 20:35:37 UTC
Permalink
Dear Justin.

If you have R1 values measured from an experiment.

And you have R1rho intensities measured as variations of spin-lock time,
spin-lock offset and spin-lock field/power, you can perform R1rho analysis
in relax.

If we just think of model DPL94, you would in relax get fitted R20(or
R1rhoprime) and Rex (determined by fitted phi_ex, kex and calculated
effective field).

Then you would like to interpolate the graph.

One could now interpolate R1rho values, by just varying theta, since we all
other parameters.

The current implementation in relax is by ramping the spin-lock field, but
for the graphs you refer to, it should instead ramp the offset.
I faced that exact same problem a week ago, and made this little page to
show how the graph change according to that:
http://wiki.nmr-relax.com/DPL94_math

So with the parameters in hand from relax, you could in a spreadsheet just
vary theta and make the graphs.

But I aimed for the auto generation of these plots as grace files, and was
stuck until realizing what needed to be modified in relax, since just
making the graphs looked wrong.

For an implementation in relax, we need something like:
Function to ramp field values
Function to ramp offset values
Make function that return calculated R1rho data accept interpolated data
Make back calculation accept as well
Expand the grace write function

So, it takes a little before we are there. :) But now there is a plan.

Best
Troels
Post by Justin Lecher
Post by Edward d'Auvergne
Hi,
The key here is what is meant by R2. There are many different
- In the standard relaxation equations, spin-spin relaxation
component. I.e. the part influenced by the spectral density function
J(w).
- Again in the standard relaxation equations, R2 is used for R2*,
where R2* = R2 + Rex. R2* is also defined differently if you look at
the exponential decay curves or if you looking at peak widths (for the
later there are additional factors broadening the peaks which add to
the R2* equation).
- For relaxation dispersion, R2 is currently defined as R1rho', i.e.
just the parts influenced by the spectral density function. Here
R1rho' is the on-resonance part of R1rho, excluding exchange.
R2 = R1rho' + Rex could also be defined, which is probably what you
are after. The value of R1rho' + Rex has been named many different
things by the field and there seems to be no consensus. I have
discussed this with Troels at
http://thread.gmane.org/gmane.science.nmr.relax.devel/5119/focus=5207.
This is a long thread with many discussions about implementing this
as an automatically calculated parameter - which in the end would be a
great feature.
So which R2 value are you after?
Regards,
Edward
Hi Edward,
Principally I am looking for pure R2 without Rex contribution.
Currently we are using some home made scripts which basically are using
equation 20 in Palmer, Massi, 2006 (PMID 16683750).
So for the technical side, how do I proceed in relax? Do I need to fix
the decay for R1 and R1rho and then calculate R2 or just fit R1 and use
my raw intensities from the R1rho measurements for the R2 calculation?
Thanks,
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
Justin Lecher
2014-03-27 07:06:03 UTC
Permalink
Post by Troels Emtekær Linnet
Dear Justin.
If you have R1 values measured from an experiment.
And you have R1rho intensities measured as variations of spin-lock time,
spin-lock offset and spin-lock field/power, you can perform R1rho
analysis in relax.
If we just think of model DPL94, you would in relax get fitted R20(or
R1rhoprime) and Rex (determined by fitted phi_ex, kex and calculated
effective field).
Then you would like to interpolate the graph.
One could now interpolate R1rho values, by just varying theta, since we
all other parameters.
The current implementation in relax is by ramping the spin-lock field,
but for the graphs you refer to, it should instead ramp the offset.
I faced that exact same problem a week ago, and made this little page to
http://wiki.nmr-relax.com/DPL94_math
So with the parameters in hand from relax, you could in a spreadsheet
just vary theta and make the graphs.
But I aimed for the auto generation of these plots as grace files, and
was stuck until realizing what needed to be modified in relax, since
just making the graphs looked wrong.
Function to ramp field values
Function to ramp offset values
Make function that return calculated R1rho data accept interpolated data
Make back calculation accept as well
Expand the grace write function
So, it takes a little before we are there. :) But now there is a plan.
Best
Troels
Thanks a lot, Troels.

To things I learned, first I really need to learn more on the whole
theory and second I will try to do the full R1rho analysis and see what
the models bring to me. I will also have a look into the Matplotlib
script and see how it helps.

Thanks,
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
Edward d'Auvergne
2014-03-27 08:25:28 UTC
Permalink
Hi Justin,

For the models you use, you need to be careful and selective. Almost
all dispersion models are implemented in relax, but that does not mean
that they should all be used
(http://www.nmr-relax.com/manual/Dispersion_GUI_mode_choosing_models_optimise.html).
You can think of relax as a Swiss army knife - just because scissors
are there, that does not mean you should unscrew screws with it ;) Be
sure to read the reference paper for any model you wish to use - they
will not all be suitable for your data or your system. Easy links to
these references can be found on the relax wiki
(http://wiki.nmr-relax.com/Category:Relaxation_dispersion).

Regards,

Edward
Post by Justin Lecher
Post by Troels Emtekær Linnet
Dear Justin.
If you have R1 values measured from an experiment.
And you have R1rho intensities measured as variations of spin-lock time,
spin-lock offset and spin-lock field/power, you can perform R1rho
analysis in relax.
If we just think of model DPL94, you would in relax get fitted R20(or
R1rhoprime) and Rex (determined by fitted phi_ex, kex and calculated
effective field).
Then you would like to interpolate the graph.
One could now interpolate R1rho values, by just varying theta, since we
all other parameters.
The current implementation in relax is by ramping the spin-lock field,
but for the graphs you refer to, it should instead ramp the offset.
I faced that exact same problem a week ago, and made this little page to
http://wiki.nmr-relax.com/DPL94_math
So with the parameters in hand from relax, you could in a spreadsheet
just vary theta and make the graphs.
But I aimed for the auto generation of these plots as grace files, and
was stuck until realizing what needed to be modified in relax, since
just making the graphs looked wrong.
Function to ramp field values
Function to ramp offset values
Make function that return calculated R1rho data accept interpolated data
Make back calculation accept as well
Expand the grace write function
So, it takes a little before we are there. :) But now there is a plan.
Best
Troels
Thanks a lot, Troels.
To things I learned, first I really need to learn more on the whole
theory and second I will try to do the full R1rho analysis and see what
the models bring to me. I will also have a look into the Matplotlib
script and see how it helps.
Thanks,
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
Edward d'Auvergne
2014-03-27 08:19:39 UTC
Permalink
Hi Troels,
Post by Troels Emtekær Linnet
But I aimed for the auto generation of these plots as grace files, and was
stuck until realizing what needed to be modified in relax, since just making
the graphs looked wrong.
It's a pity that this doesn't work yet because the implementation is
so close. All of the planning has been completed and you have written
most of the code for the required infrastructure.
Post by Troels Emtekær Linnet
Function to ramp field values
This spin-lock field strength interpolation (and extension) code
exists, but just needs to be shifted into its own function to simplify
things. Well, there is a different interpolation used for the
analytic and numeric CPMG models for the nu_CPMG values, so additional
functions would be required (this could be put into one or two
functions, that is not important).
Post by Troels Emtekær Linnet
Function to ramp offset values
After shifting the field strength interpolation into its own function,
this can be duplicated and changed from the spin-lock field strength
into the spin-lock offset.
Post by Troels Emtekær Linnet
Make function that return calculated R1rho data accept interpolated data
This already exists. The current code of the back_calc_r2eff()
function only needs a single line of code for the calculation of the
offsets shifted into the earlier function where the interpolation
happens, and then the data passed back into this function with or
without interpolation, and then everything will just work.
Post by Troels Emtekær Linnet
Make back calculation accept as well
There is nothing to do here, once the single line of code has been
shifted. If the interpolated data structures are created and passed
into it, this will just work. This code will not even know the
difference.
Post by Troels Emtekær Linnet
Expand the grace write function
This also requires zero work, everything is already in place and the
code should not be changed :) If you mean the
relax_disp.plot_disp_curves user function, then this just needs to
call the interpolation functions above to change the list of offsets.
If the matplotlib plotting infrastructure was added to relax, it would
be then easy to have this user function create a 3D surface graph of
the R1rho dispersion curves.

Cheers,

Edward
Troels Emtekær Linnet
2014-03-26 13:09:58 UTC
Permalink
Well, you can chose to try all models.

Then relax will do a model selection, due to AIC criterion.
Relax will make model.out file in the "final" directory, and perform monte
carlo error analysis.

Then one can then consider doing a global fitting/clustering analysis for
some residues, choosing one of the models.
Reading in the old results directory from the previous run. That will
average the values for the global parameter kex.
This will probably modify the R2 values, since they are now fitted under
restriction of the global parameters.

If have not yet completed a tutorial for setting up R1rho.
http://wiki.nmr-relax.com/Tutorial_for_Relaxation_dispersion_analysis_r1rho_fixed_time_recorded_on_varian_as_sequential_spectra

But you can see the system test in:
test_suite/system_tests/relax_disp.py under def test_r1rho_kjaergaard(self):

Best
Troels
Post by Troels Emtekær Linnet
Post by Troels Emtekær Linnet
Hi Justin.
http://wiki.nmr-relax.com/DPL94
http://wiki.nmr-relax.com/TAP03
http://wiki.nmr-relax.com/TP02
http://wiki.nmr-relax.com/NS_R1rho_2-site
http://wiki.nmr-relax.com/Tutorial_for_R1/R2_Relaxation_curve-fitting_analysis_on_varian_recorded_as_fid_interleaved
Post by Troels Emtekær Linnet
(A tutorial I made for our students in our lab. We have Varian at our
place).
Best
Troels
Hi Troels,
perfect, thanks! I found the wiki pages, but I wasn't sure which one is
the best to choose. Any general suggestion or does it strictly depend on
the type of model which applies best to my system?
thanks
Justin
--
Justin Lecher
Institute of Complex Systems
ICS-6 Structural Biochemistry
Research Centre Juelich
52425 Juelich, Germany
phone: +49 2461 61 2117
Continue reading on narkive:
Loading...