Improved description of chemical barriers with generalized gradient approximations (GGAs) and meta-GGAs

The performance of exchange and correlation (xc) functionals of the generalized gradient approximation (GGA) type and of the meta-GGA type in the calculation of chemical reactions is related to topological features of the electron density which, in turn, are connected to the orbital structure of chemical bonds within the Kohn-Sham (KS) theory. Seventeen GGA and meta-GGA xc functionals are assessed for 15 hydrogen abstraction reactions and 3 symmetrical S(N)2 reactions. Systems that are problematic for standard GGAs characteristically have enhanced values of the dimensionless gradient argument s(sigma)(2) with local maxima in the bonding region. The origin of this topological feature is the occupation of valence KS orbitals with an antibonding or essentially nonbonding character. The local enhancement of s(sigma)(2) yields too negative exchange-correlation energies with standard GGAs for the transition state of the S(N)2 reaction, which leads to the reduced calculated reaction barriers. The unwarranted localization of the effective xc hole of the standard GGAs, i.e., the nondynamical correlation that is built into them but is spurious in this case, wields its effect by their s(sigma)(2) dependence. Barriers are improved for xc functionals with the exchange functional OPTX as x component, which has a modified dependence on s(sigma)(2). Standard GGAs also underestimate the barriers for the hydrogen abstraction reactions. In this case the barriers are improved by correlation functionals, such as the Laplacian-dependent (LAP3) functional, which has a modified dependence on the Coulomb correlation of the opposite- and like-spin electrons. The best overall performance is established for the combination OLAP3 of OPTX and LAP3.

On the required shape corrections to the local density and generalized gradient approximations to the Kohn-Sham potentials for molecular response calculations of (hyper)polarizabilities and excitation energies

It is well known that shape corrections have to be applied to the local-density (LDA) and generalized gradient (GGA) approximations to the Kohn-Sham exchange-correlation potential in order to obtain reliable response properties in time dependent density functional theory calculations. Here we demonstrate that it is an oversimplified view that these shape corrections concern primarily the asymptotic part of the potential, and that they affect only Rydberg type transitions. The performance is assessed of two shape-corrected Kohn-Sham potentials, the gradient-regulated asymptotic connection procedure applied to the Becke-Perdew potential (BP-GRAC) and the statistical averaging of (model) orbital potentials (SAOP), versus LDA and GGA potentials, in molecular response calculations of the static average polarizability alpha, the Cauchy coefficient S-4, and the static average hyperpolarizability beta. The nature of the distortions of the LDA/GGA potentials is highlighted and it is shown that they introduce many spurious excited states at too low energy which may mix with valence excited states, resulting in wrong excited state compositions. They also lead to wrong oscillator strengths and thus to a wrong spectral structure of properties like the polarizability. LDA, Becke-Lee-Yang-Parr (BLYP), and Becke-Perdew (BP) characteristically underestimate contributions to alpha and S-4 from bound Rydberg-type states and overestimate those from the continuum. Cancellation of the errors in these contributions occasionally produces fortuitously good results. The distortions of the LDA, BLYP, and BP spectra are related to the deficiencies of the LDA/GGA potentials in both the bulk and outer molecular regions. In contrast, both SAOP and BP-GRAC potentials produce high quality polarizabilities for 21 molecules and also reliable Cauchy moments and hyperpolarizabilities for the selected molecules. The analysis for the N-2 molecule shows, that both SAOP and BP-GRAC yield reliable energies omega(i) and oscillator strengths f(i) of individual excitations, so that they reproduce well the spectral structure of alpha and S-4.(C) 2002 American Institute of Physics.

The failure of generalized gradient approximations (GGAs) and meta-GGAs for the two-center three-electron bonds in He-2(+), (H2O)(2)(+), and (NH3)(2)(+)

The radical cations He-2(+) (H2O)(2)(+), and (NH3)(2)(+) with two-center three-electron A-A bonds are investigated at the configuration interaction (CI), accurate Kohn-Sham (KS), generalized gradient approximation (GGA), and meta-GGA levels. Assessment of seven different GGA and six meta-GGA methods shows that the A(2)(+) systems remain a difficult case for density functional theory (DFT). All methods tested consistently overestimate the stability of A(2)(+): the corresponding D-e errors decrease for more diffuse valence densities in the series He-2(+) > (H2O)(2)(+) > (NH3)(2)(+). Upon comparison to the energy terms of the accurate Kohn-Sham solutions, the approximate exchange functionals are found to be responsible for the errors of GGA-type methods, which characteristically overestimate the exchange in A(2)(+). These so-called exchange functionals implicitly use localized holes. Such localized holes do occur if there is left-right correlation, i.e., the exchange functionals then also describe nondynamical correlation. However, in the hemibonded A(2)(+) systems the typical molecular (left-right, nondynamical) correlation of the two-electron pair bond is absent. The nondynamical correlation built into the exchange functionals is then spurious and yields too low energies.

Translation suppression promotes stress granule formation and cell survival in response to cold shock

Cells respond to different types of stress by inhibition of protein synthesis and subsequent assembly of stress granules (SGs), cytoplasmic aggregates that contain stalled translation preinitiation complexes. Global translation is regulated through the translation initiation factor eukaryotic initiation factor 2a (eIF2a) and the mTOR pathway. Here we identify cold shock as a novel trigger of SG assembly in yeast and mammals. Whereas cold shock-induced SGs take hours to form, they dissolve within minutes when cells are returned to optimal growth temperatures. Cold shock causes eIF2a phosphorylation through the kinase PERK in mammalian cells, yet this pathway is not alone responsible for translation arrest and SG formation. In addition, cold shock leads to reduced mitochondrial function, energy depletion, concomitant activation of AMP-activated protein kinase (AMPK), and inhibition of mTOR signaling. Compound C, a pharmacological inhibitor of AMPK, prevents the formation of SGs and strongly reduces cellular survival in a translation-dependent manner. Our results demonstrate that cells actively suppress protein synthesis by parallel pathways, which induce SG formation and ensure cellular survival during hypothermia.

Dynamic oscillation of translation and stress granule formation mark the cellular response to virus infection

Virus infection-induced global protein synthesis suppression is linked to assembly of stress granules (SGs), cytosolic aggregates of stalled translation preinitiation complexes. To study long-term stress responses, we developed an imaging approach for extended observation and analysis of SG dynamics during persistent hepatitis C virus (HCV) infection. In combination with type 1 interferon, HCV infection induces highly dynamic assembly/disassembly of cytoplasmic SGs, concomitant with phases of active and stalled translation, delayed cell division, and prolonged cell survival. Double-stranded RNA (dsRNA), independent of viral replication, is sufficient to trigger these oscillations. Translation initiation factor eIF2a phosphorylation by protein kinase R mediates SG formation and translation arrest. This is antagonized by the upregulation of GADD34, the regulatory subunit of protein phosphatase 1 dephosphorylating eIF2a. Stress response oscillation is a general mechanism to prevent long-lasting translation repression and a conserved host cell reaction to multiple RNA viruses, which HCV may exploit to establish persistence.

Modeling continuous self-report measures of perceived emotion using generalized additive mixed models

Emotion research has long been dominated by the “standard method” of displaying posed or acted static images of facial expressions of emotion. While this method has been useful it is unable to investigate the dynamic nature of emotion expression. Although continuous self-report traces have enabled the measurement of dynamic expressions of emotion, a consensus has not been reached on the correct statistical techniques that permit inferences to be made with such measures. We propose Generalized Additive Models and Generalized Additive Mixed Models as techniques that can account for the dynamic nature of such continuous measures. These models allow us to hold constant shared components of responses that are due to perceived emotion across time, while enabling inference concerning linear differences between groups. The mixed model GAMM approach is preferred as it can account for autocorrelation in time series data and allows emotion decoding participants to be modelled as random effects. To increase confidence in linear differences we assess the methods that address interactions between categorical variables and dynamic changes over time. In addition we provide comments on the use of Generalized Additive Models to assess the effect size of shared perceived emotion and discuss sample sizes. Finally we address additional uses, the inference of feature detection, continuous variable interactions, and measurement of ambiguity.

Age-dependent increase in ortho-tyrosine and methionine sulfoxide in human skin collagen is not accelerated in diabetes. Evidence against a generalized increase in oxidative stress in diabetes

The glycoxidation products Nepsilon-(carboxymethyl)lysine and pentosidine increase in skin collagen with age and at an accelerated rate in diabetes. Their age-adjusted concentrations in skin collagen are correlated with the severity of diabetic complications. To determine the relative roles of increased glycation and/or oxidation in the accelerated formation of glycoxidation products in diabetes, we measured levels of amino acid oxidation products, distinct from glycoxidative modifications of amino acids, as independent indicators of oxidative stress and damage to collagen in aging and diabetes. We show that ortho-tyrosine and methionine sulfoxide are formed in concert with Nepsilon-(carboxymethyl)lysine and pentosidine during glycoxidation of collagen in vitro, and that they also increase with age in human skin collagen. The age-adjusted levels of these oxidized amino acids in collagen was the same in diabetic and nondiabetic subjects, arguing that diabetes per se does not cause an increase in oxidative stress or damage to extracellular matrix proteins. These results provide evidence for an age-dependent increase in oxidative damage to collagen and support previous conclusions that the increase in glycoxidation products in skin collagen in diabetes can be explained by the increase in glycemia alone, without invoking a generalized, diabetes-dependent increase in oxidative stress.

A generalized functional response for predators that switch between multiple prey species

We develop a theory for the food intake of a predator that can switch between multiple prey species. The theory addresses empirical observations of prey switching and is based on the behavioural assumption that a predator tends to continue feeding on prey that are similar to the prey it has consumed last, in terms of, e.g., their morphology, defences, location, habitat choice, or behaviour. From a predator's dietary history and the assumed similarity relationship among prey species, we derive a general closed-form multi-species functional response for describing predators switching between multiple prey species. Our theory includes the Holling type II functional response as a special case and makes consistent predictions when populations of equivalent prey are aggregated or split. An analysis of the derived functional response enables us to highlight the following five main findings. (1) Prey switching leads to an approximate power-law relationship between ratios of prey abundance and prey intake, consistent with experimental data. (2) In agreement with empirical observations, the theory predicts an upper limit of 2 for the exponent of such power laws. (3) Our theory predicts deviations from power-law switching at very low and very high prey-abundance ratios. (4) The theory can predict the diet composition of a predator feeding on multiple prey species from diet observations for predators feeding only on pairs of prey species. (5) Predators foraging on more prey species will show less pronounced prey switching than predators foraging on fewer prey species, thus providing a natural explanation for the known difficulties of observing prey switching in the field. (C) 2013 Elsevier Ltd. All rights reserved.