An Introduction to the 'Special Volume Spectroscopy and Chemometrics in R'

This special volume collates ten issues under the rubric "Spectroscopy and Chemometrics in R". In so doing, it provides an overview of the breadth, depth and state of the art of R-based software projects for spectroscopy and chemometrics applications. Just as the authors have contributed to R their documentation and source code, so has R contributed to the quality, standardization and dissemination of their software, as this volume attests. We hope that the volume is inspiring to both computational statisticians interested in applications of their methodologies and to spectroscopists or chemometricians in need of solutions to their data analysis problems.

Global and target analysis of time-resolved spectra

AbstractIn biological/bioenergetics research the response of a complex system to an externally applied perturbation is often studied. Spectroscopic measurements at multiple wavelengths are used to monitor the kinetics. These time-resolved spectra are considered as an example of multiway data. In this paper, the methodology for global and target analysis of time-resolved spectra is reviewed. To fully extract the information from the overwhelming amount of data, a model-based analysis is mandatory. This analysis is based upon assumptions regarding the measurement process and upon a physicochemical model for the complex system. This model is composed of building blocks representing scientific knowledge and assumptions. Building blocks are the instrument response function (IRF), the components of the system connected in a kinetic scheme, and anisotropy properties of the components. The combination of a model for the kinetics and for the spectra of the components results in a more powerful spectrotemporal model. The model parameters, like rate constants and spectra, can be estimated from the data, thus providing a concise description of the complex system dynamics. This spectrotemporal modeling approach is illustrated with an elaborate case study of the ultrafast dynamics of the photoactive yellow protein

Mechanisms of drought-induced dissipation of excitation energy in sun- and shade-adapted drought-tolerant mosses studied by fluorescence yield change and global and target analysis of fluorescence decay kinetics

Some mosses stay green and survive long even under desiccation. Dissipation mechanisms of excess excitation energy were studied in two drought-tolerant moss species adapted to contrasting niches: shade-adapted Rhytidiadelphus squarrosus and sun-adapted Rhytidium rugosum in the same family. (1) Under wet conditions, a light-induced nonphotochemical quenching (NPQ) mechanism decreased the yield of photosystem II (PSII) fluorescence in both species. The NPQ extent saturated at a lower illumination intensity in R. squarrosus, suggesting a larger PSII antenna size. (2) Desiccation reduced the fluorescence intensities giving significantly lower F 0 levels and shortened the overall fluorescence lifetimes in both R. squarrosus and R. rugosum, at room temperature. (3) At 77 K, desiccation strongly reduced the PSII fluorescence intensity. This reduction was smaller in R. squarrosus than in R. rugosum. (4) Global and target analysis indicated two different mechanisms of energy dissipation in PSII under desiccation: the energy dissipation to a desiccation-formed strong fluorescence quencher in the PSII core in sun-adapted R. rugosum (type-A quenching) and (5) the moderate energy dissipation in the light-harvesting complex/PSII in shade-adapted R. squarrosus (type-B quenching). The two mechanisms are consistent with the different ecological niches of the two mosses

Vibronic Wavepackets and Energy Transfer in Cryptophyte Light-Harvesting Complexes

Determining the key features of high-efficiency photosynthetic energy transfer remains an ongoing task. Recently, there has been evidence for the role of vibronic coherence in linking donor and acceptor states to redistribute oscillator strength for enhanced energy transfer. To gain further insights into the interplay between vibronic wavepackets and energy-transfer dynamics, we systematically compare four structurally related phycobiliproteins from cryptophyte algae by broad-band pump-probe spectroscopy and extend a parametric model based on global analysis to include vibrational wavepacket characterization. The four phycobiliproteins isolated from cryptophyte algae are two "open" structures and two "closed" structures. The closed structures exhibit strong exciton coupling in the central dimer. The dominant energy-transfer pathway occurs on the subpicosecond timescale across the largest energy gap in each of the proteins, from central to peripheral chromophores. All proteins exhibit a strong 1585 cm-1 coherent oscillation whose relative amplitude, a measure of vibronic intensity borrowing from resonance between donor and acceptor states, scales with both energy-transfer rates and damping rates. Central exciton splitting may aid in bringing the vibronically linked donor and acceptor states into better resonance resulting in the observed doubled rate in the closed structures. Several excited-state vibrational wavepackets persist on timescales relevant to energy transfer, highlighting the importance of further investigation of the interplay between electronic coupling and nuclear degrees of freedom in studies on high-efficiency photosynthesis

FluxSimulator: An R Package to Simulate Isotopomer Distributions in Metabolic Networks

The representation of biochemical knowledge in terms of fluxes (transformation rates) in a metabolic network is often a crucial step in the development of new drugs and efficient bioreactors. Mass spectroscopy (MS) and nuclear magnetic resonance spectroscopy (NMRS) in combination with ^13C labeled substrates are experimental techniques resulting in data that may be used to quantify fluxes in the metabolic network underlying a process. The massive amount of data generated by spectroscopic experiments increasingly requires software which models the dynamics of the underlying biological system. In this work we present an approach to handle isotopomer distributions in metabolic networks using an object-oriented programming approach, implemented using S4 classes in R. The developed package is called FluxSimulator and provides a user friendly interface to specify the topological information of the metabolic network as well as carbon atom transitions in plain text files. The package automatically derives the mathematical representation of the formulated network, and assembles a set of ordinary differential equations (ODEs) describing the change of each isotopomer pool over time. These ODEs are subsequently solved numerically. In a case study FluxSimulator was applied to an example network. Our results indicate that the package is able to reproduce exact changes in isotopomer compositions of the metabolite pools over time at given flux rates.

Exploring the mechanism(s) of energy dissipation in the light harvesting complex of the photosynthetic algae Cyclotella meneghiniana

AbstractPhotosynthetic organisms have developed vital strategies which allow them to switch from a light-harvesting to an energy dissipative state at the level of the antenna system in order to survive the detrimental effects of excess light illumination. These mechanisms are particularly relevant in diatoms, which grow in highly fluctuating light environments and thus require fast and strong response to changing light conditions. We performed transient absorption spectroscopy on FCPa, the main light-harvesting antenna from the diatom Cyclotella meneghiniana, in the unquenched and quenched state. Our results show that in quenched FCPa two quenching channels are active and are characterized by differing rate constants and distinct spectroscopic signatures. One channel is associated with a faster quenching rate (16ns−1) and virtually no difference in spectral shape compared to the bulk unquenched chlorophylls, while a second channel is associated with a slower quenching rate (2.7ns−1) and exhibits an increased population of red-emitting states. We discuss the origin of the two processes in the context of the models proposed for the regulation of photosynthetic light-harvesting. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: Keys to Produce Clean Energy

Pathways for Energy Transfer in the Core Light-Harvesting Complexes CP43 and CP47 of Photosystem II

AbstractThe pigment-protein complexes CP43 and CP47 transfer excitation energy from the peripheral antenna of photosystem II toward the photochemical reaction center. We measured the excitation dynamics of the chlorophylls in isolated CP43 and CP47 complexes at 77K by time-resolved absorbance-difference and fluorescence spectroscopy. The spectral relaxation appeared to occur with rates of 0.2–0.4ps and 2–3ps in both complexes, whereas an additional relaxation of 17ps was observed only in CP47. Using the 3.8-Å crystal structure of the photosystem II core complex from Synechococcus elongatus (A. Zouni, H.-T. Witt, J. Kern, P. Fromme, N. Krauss, W. Saenger, and P. Orth, 2001, Nature, 409:739–743), excitation energy transfer kinetics were calculated and a Monte Carlo simulation of the absorption spectra was performed. In both complexes, the rate of 0.2–0.4ps can be ascribed to excitation energy transfer within a layer of chlorophylls near the stromal side of the membrane, and the slower 2–3-ps process to excitation energy transfer to the calculated lowest excitonic state. We conclude that excitation energy transfer within CP43 and CP47 is fast and does not contribute significantly to the well-known slow trapping of excitation energy in photosystem II

Fluorescence Lifetime Imaging Microscopy (FLIM) Data Analysis with TIMP

Fluorescence Lifetime Imaging Microscopy (FLIM) allows fluorescence lifetime images of biological objects to be collected at 250 nm spatial resolution and at (sub-)nanosecond temporal resolution. Often n_comp kinetic processes underlie the observed fluorescence at all locations, but the intensity of the fluorescence associated with each process varies per-location, i.e., per-pixel imaged. Then the statistical challenge is global analysis of the image: use of the fluorescence decay in time at all locations to estimate the n_comp lifetimes associated with the kinetic processes, as well as the amplitude of each kinetic process at each location. Given that typical FLIM images represent on the order of 10^2 timepoints and 10^3 locations, meeting this challenge is computationally intensive. Here the utility of the TIMP package for R to solve parameter estimation problems arising in FLIM image analysis is demonstrated. Case studies on simulated and real data evidence the applicability of the partitioned variable projection algorithm implemented in TIMP to the problem domain, and showcase options included in the package for the visual validation of models for FLIM data.