PyLDM - An open source package for lifetime density analysis of time-resolved spectroscopic data

Ultrafast spectroscopy offers temporal resolution for probing processes in the femto- and picosecond regimes. This has allowed for investigation of energy and charge transfer in numerous photoactive compounds and complexes. However, analysis of the resultant data can be complicated, particularly in more complex biological systems, such as photosystems. Historically, the dual approach of global analysis and target modelling has been used to elucidate kinetic descriptions of the system, and the identity of transient species respectively. With regards to the former, the technique of lifetime density analysis (LDA) offers an appealing alternative. While global analysis approximates the data to the sum of a small number of exponential decays, typically on the order of 2-4, LDA uses a semi-continuous distribution of 100 lifetimes. This allows for the elucidation of lifetime distributions, which may be expected from investigation of complex systems with many chromophores, as opposed to averages. Furthermore, the inherent assumption of linear combinations of decays in global analysis means the technique is unable to describe dynamic motion, a process which is resolvable with LDA. The technique was introduced to the field of photosynthesis over a decade ago by the Holzwarth group. The analysis has been demonstrated to be an important tool to evaluate complex dynamics such as photosynthetic energy transfer, and complements traditional global and target analysis techniques. Although theory has been well described, no open source code has so far been available to perform lifetime density analysis. Therefore, we introduce a python (2.7) based package, PyLDM, to address this need. We furthermore provide a direct comparison of the capabilities of LDA with those of the more familiar global analysis, as well as providing a number of statistical techniques for dealing with the regularization of noisy data.

[1]  Aleksander Siemiarczuk,et al.  Comparison of the maximum entropy and exponential series methods for the recovery of distributions of lifetimes from fluorescence lifetime data , 1990 .

[2]  R. Tibshirani,et al.  Sparsity and smoothness via the fused lasso , 2005 .

[3]  L. Eldén Algorithms for the regularization of ill-conditioned least squares problems , 1977 .

[4]  Gabriel F. Dorlhiac,et al.  Femtosecond Visible Transient Absorption Spectroscopy of Chlorophyll f-Containing Photosystem I , 2017, Biophysical journal.

[5]  W. Lubitz,et al.  Ultrafast transient absorption studies on photosystem I reaction centers from Chlamydomonas reinhardtii. 2: mutations near the P700 reaction center chlorophylls provide new insight into the nature of the primary electron donor. , 2006, Biophysical journal.

[6]  W. Lubitz,et al.  Ultrafast transient absorption studies on Photosystem I reaction centers from Chlamydomonas reinhardtii. 1. A new interpretation of the energy trapping and early electron transfer steps in Photosystem I. , 2003, Biophysical journal.

[7]  T. Ulrych,et al.  Time series modeling and maximum entropy , 1976 .

[8]  A. Livesey,et al.  Maximum entropy analysis of quasielastic light scattering from colloidal dispersions , 1986 .

[9]  Shohei Tateishi,et al.  Tuning parameter selection in sparse regression modeling , 2013, Comput. Stat. Data Anal..

[10]  J. V. van Thor,et al.  Modelling Multi-Pulse Population Dynamics from Ultrafast Spectroscopy , 2011, PloS one.

[11]  M. G. Müller,et al.  Kinetics and mechanism of electron transfer in intact photosystem II and in the isolated reaction center: pheophytin is the primary electron acceptor. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[12]  A. J. Mello,et al.  Evanescent wave-induced fluorescence study of Rhodamine 101 at dielectric interfaces , 1996 .

[13]  G. Zucchelli,et al.  Energy transfer processes in the isolated core antenna complexes CP43 and CP47 of photosystem II. , 2010, Biochimica et biophysica acta.

[14]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[15]  H. Zou,et al.  Regularization and variable selection via the elastic net , 2005 .

[16]  Arthur E. Hoerl,et al.  Ridge Regression: Biased Estimation for Nonorthogonal Problems , 2000, Technometrics.

[17]  Eric Jones,et al.  SciPy: Open Source Scientific Tools for Python , 2001 .

[18]  R. Tibshirani,et al.  Regression shrinkage and selection via the lasso: a retrospective , 2011 .

[19]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[20]  Yongxin Pan,et al.  Physics of the Earth and Planetary Interiors , 2015 .

[21]  A. Istratov,et al.  Exponential analysis in physical phenomena , 1999 .

[22]  Lothar Reichel,et al.  Inverse problems for regularization matrices , 2012, Numerical Algorithms.

[23]  Igor Jurisica,et al.  Optimized application of penalized regression methods to diverse genomic data , 2011, Bioinform..

[24]  John T. M. Kennis,et al.  Ultrafast transient absorption spectroscopy: principles and application to photosynthetic systems , 2009, Photosynthesis Research.

[25]  Alfred R. Holzwarth,et al.  Data Analysis of Time-Resolved Measurements , 1996 .

[26]  Pascal Pernot,et al.  Comprehensive data analysis of femtosecond transient absorption spectra: A review , 2012 .

[27]  Travis E. Oliphant,et al.  Python for Scientific Computing , 2007, Computing in Science & Engineering.

[28]  Bin Dai,et al.  Orthogonalizing EM: A Design-Based Least Squares Algorithm , 2016, Technometrics.

[29]  Rienk van Grondelle,et al.  Global and target analysis of time-resolved spectra. , 2004, Biochimica et biophysica acta.

[30]  C. Slavov,et al.  Implementation and evaluation of data analysis strategies for time-resolved optical spectroscopy. , 2015, Analytical chemistry.

[31]  Gaël Varoquaux,et al.  The NumPy Array: A Structure for Efficient Numerical Computation , 2011, Computing in Science & Engineering.

[32]  HansenPer Christian The truncated SVD as a method for regularization , 1987 .

[33]  Rajiv Luthra,et al.  Independent initiation of primary electron transfer in the two branches of the photosystem I reaction center , 2010, Proceedings of the National Academy of Sciences.

[34]  Peter R. Johnston,et al.  Selecting the corner in the L-curve approach to Tikhonov regularization , 2000, IEEE Transactions on Biomedical Engineering.

[35]  Heinz W. Engl,et al.  Distribution of event times in time-resolved fluorescence: the exponential series approach—algorithm, regularization, analysis , 1991 .

[36]  R. Tibshirani Regression Shrinkage and Selection via the Lasso , 1996 .

[37]  H. Kandori,et al.  Transformation of Time-Resolved Spectra to Lifetime-Resolved Spectra by Maximum Entropy Inversion of the Laplace Transform , 2006, Applied spectroscopy.

[38]  P. Hansen The discrete picard condition for discrete ill-posed problems , 1990 .

[39]  Thomas Renger,et al.  Light harvesting in photosystem II core complexes is limited by the transfer to the trap: can the core complex turn into a photoprotective mode? , 2008, Journal of the American Chemical Society.

[40]  A. V. Schaik,et al.  L1 regularization method in electrical impedance tomography by using the L1-curve (Pareto frontier curve) , 2012 .

[41]  M. G. Müller,et al.  Carotenoid-to-chlorophyll energy transfer in recombinant major light-harvesting complex (LHCII) of higher plants. I. Femtosecond transient absorption measurements. , 2001, Biophysical journal.

[42]  A. Visser,et al.  Time-resolved fluorescence studies of dityrosine in the outer layer of intact yeast ascospores. , 1994, Biophysical journal.

[43]  Gene H. Golub,et al.  The differentiation of pseudo-inverses and non-linear least squares problems whose variables separate , 1972, Milestones in Matrix Computation.

[44]  R. Tibshirani,et al.  The solution path of the generalized lasso , 2010, 1005.1971.

[45]  Vincent A Voelz,et al.  Calculation of rate spectra from noisy time series data , 2012, Proteins.

[46]  H. Voss,et al.  Solving regularized total least squares problems based on eigenproblems , 2010 .

[47]  Michael Wahl,et al.  Time-Correlated Single Photon Counting , 2009 .

[48]  G. E. Thomas,et al.  Measurement of the Time Dependence of Scintillation Intensity by a Delayed‐Coincidence Method , 1961 .