Folding protein models with a simple hydrophobic energy function: The fundamental importance of monomer inside/outside segregation

The present study explores a “hydrophobic” energy function for folding simulations of the protein lattice model. The contribution of each monomer to conformational energy is the product of its “hydrophobicity” and the number of contacts it makes, i.e., E(h⃗, c⃗) = −Σi=1N cihi = −(h⃗.c⃗) is the negative scalar product between two vectors in N-dimensional cartesian space: h⃗ = (h1, … , hN), which represents monomer hydrophobicities and is sequence-dependent; and c⃗ = (c1, … , cN), which represents the number of contacts made by each monomer and is conformation-dependent. A simple theoretical analysis shows that restrictions are imposed concomitantly on both sequences and native structures if the stability criterion for protein-like behavior is to be satisfied. Given a conformation with vector c⃗, the best sequence is a vector h⃗ on the direction upon which the projection of c⃗ − c̄⃗ is maximal, where c̄⃗ is the diagonal vector with components equal to c̄, the average number of contacts per monomer in the unfolded state. Best native conformations are suggested to be not maximally compact, as assumed in many studies, but the ones with largest variance of contacts among its monomers, i.e., with monomers tending to occupy completely buried or completely exposed positions. This inside/outside segregation is reflected on an apolar/polar distribution on the corresponding sequence. Monte Carlo simulations in two dimensions corroborate this general scheme. Sequences targeted to conformations with large contact variances folded cooperatively with thermodynamics of a two-state transition. Sequences targeted to maximally compact conformations, which have lower contact variance, were either found to have degenerate ground state or to fold with much lower cooperativity.

Implications of a possible clustering of highest-energy cosmic rays

Recently, a possible clustering of a subset of observed ultra-high energy cosmic rays above ≃40 EeV (4 × 1019 eV) in pairs near the supergalactic plane was reported. We show that a confirmation of this effect would provide information on the origin and nature of these events and, in case of charged primaries, imply interesting constraints on the extragalactic magnetic field. Possible implications for the most common models of ultra-high energy cosmic ray production in the literature are discussed.

d-β-Hydroxybutyrate protects neurons in models of Alzheimer's and Parkinson's disease

The heroin analogue 1-methyl-4-phenylpyridinium, MPP+, both in vitro and in vivo, produces death of dopaminergic substantia nigral cells by inhibiting the mitochondrial NADH dehydrogenase multienzyme complex, producing a syndrome indistinguishable from Parkinson's disease. Similarly, a fragment of amyloid protein, Aβ1–42, is lethal to hippocampal cells, producing recent memory deficits characteristic of Alzheimer's disease. Here we show that addition of 4 mM d-β-hydroxybutyrate protected cultured mesencephalic neurons from MPP+ toxicity and hippocampal neurons from Aβ1–42 toxicity. Our previous work in heart showed that ketone bodies, normal metabolites, can correct defects in mitochondrial energy generation. The ability of ketone bodies to protect neurons in culture suggests that defects in mitochondrial energy generation contribute to the pathophysiology of both brain diseases. These findings further suggest that ketone bodies may play a therapeutic role in these most common forms of human neurodegeneration.

The vibrational energy flow transition in organic molecules: Theory meets experiment

Most large dynamical systems are thought to have ergodic dynamics, whereas small systems may not have free interchange of energy between degrees of freedom. This assumption is made in many areas of chemistry and physics, ranging from nuclei to reacting molecules and on to quantum dots. We examine the transition to facile vibrational energy flow in a large set of organic molecules as molecular size is increased. Both analytical and computational results based on local random matrix models describe the transition to unrestricted vibrational energy flow in these molecules. In particular, the models connect the number of states participating in intramolecular energy flow to simple molecular properties such as the molecular size and the distribution of vibrational frequencies. The transition itself is governed by a local anharmonic coupling strength and a local state density. The theoretical results for the transition characteristics compare well with those implied by experimental measurements using IR fluorescence spectroscopy of dilution factors reported by Stewart and McDonald [Stewart, G. M. & McDonald, J. D. (1983) J. Chem. Phys. 78, 3907–3915].

Effects of elevated atmospheric CO2 concentration on leaf dark respiration of Xanthium strumarium in light and in darkness

Leaf dark respiration (R) is an important component of plant carbon balance, but the effects of rising atmospheric CO2 on leaf R during illumination are largely unknown. We studied the effects of elevated CO2 on leaf R in light (RL) and in darkness (RD) in Xanthium strumarium at different developmental stages. Leaf RL was estimated by using the Kok method, whereas leaf RD was measured as the rate of CO2 efflux at zero light. Leaf RL and RD were significantly higher at elevated than at ambient CO2 throughout the growing period. Elevated CO2 increased the ratio of leaf RL to net photosynthesis at saturated light (Amax) when plants were young and also after flowering, but the ratio of leaf RD to Amax was unaffected by CO2 levels. Leaf RN was significantly higher at the beginning but significantly lower at the end of the growing period in elevated CO2-grown plants. The ratio of leaf RL to RD was used to estimate the effect of light on leaf R during the day. We found that light inhibited leaf R at both CO2 concentrations but to a lesser degree for elevated (17–24%) than for ambient (29–35%) CO2-grown plants, presumably because elevated CO2-grown plants had a higher demand for energy and carbon skeletons than ambient CO2-grown plants in light. Our results suggest that using the CO2 efflux rate, determined by shading leaves during the day, as a measure for leaf R is likely to underestimate carbon loss from elevated CO2-grown plants.

The Xanthophyll Cycle Modulates the Kinetics of Nonphotochemical Energy Dissipation in Isolated Light-Harvesting Complexes, Intact Chloroplasts, and Leaves of Spinach1

We analyzed the kinetics of nonphotochemical quenching of chlorophyll fluorescence (qN) in spinach (Spinacia oleracea) leaves, chloroplasts, and purified light-harvesting complexes. The characteristic biphasic pattern of fluorescence quenching in dark-adapted leaves, which was removed by preillumination, was evidence of light activation of qN, a process correlated with the de-epoxidation state of the xanthophyll cycle carotenoids. Chloroplasts isolated from dark-adapted and light-activated leaves confirmed the nature of light activation: faster and greater quenching at a subsaturating transthylakoid pH gradient. The light-harvesting chlorophyll a/b-binding complexes of photosystem II were isolated from dark-adapted and light-activated leaves. When isolated from light-activated leaves, these complexes showed an increase in the rate of quenching in vitro compared with samples prepared from dark-adapted leaves. In all cases, the quenching kinetics were fitted to a single component hyperbolic function. For leaves, chloroplasts, and light-harvesting complexes, the presence of zeaxanthin was associated with an increased rate constant for the induction of quenching. We discuss the significance of these observations in terms of the mechanism and control of qN.

Cosmic microwave background theory

A long-standing goal of theorists has been to constrain cosmological parameters that define the structure formation theory from cosmic microwave background (CMB) anisotropy experiments and large-scale structure (LSS) observations. The status and future promise of this enterprise is described. Current band-powers in ℓ-space are consistent with a ΔT flat in frequency and broadly follow inflation-based expectations. That the levels are ∼(10−5)2 provides strong support for the gravitational instability theory, while the Far Infrared Absolute Spectrophotometer (FIRAS) constraints on energy injection rule out cosmic explosions as a dominant source of LSS. Band-powers at ℓ ≳ 100 suggest that the universe could not have re-ionized too early. To get the LSS of Cosmic Background Explorer (COBE)-normalized fluctuations right provides encouraging support that the initial fluctuation spectrum was not far off the scale invariant form that inflation models prefer: e.g., for tilted Λ cold dark matter sequences of fixed 13-Gyr age (with the Hubble constant H0 marginalized), ns = 1.17 ± 0.3 for Differential Microwave Radiometer (DMR) only; 1.15 ± 0.08 for DMR plus the SK95 experiment; 1.00 ± 0.04 for DMR plus all smaller angle experiments; 1.00 ± 0.05 when LSS constraints are included as well. The CMB alone currently gives weak constraints on Λ and moderate constraints on Ωtot, but theoretical forecasts of future long duration balloon and satellite experiments are shown which predict percent-level accuracy among a large fraction of the 10+ parameters characterizing the cosmic structure formation theory, at least if it is an inflation variant.

Turbulent diffusion phase transition is due to singular energy spectrum.

The phase transition for turbulent diffusion, reported by Avellaneda and Majda [Avellaneda, M. & Majda, A. J. (1994) Philos. Trans. R. Soc. London A 346, 205-233, and several earlier papers], is traced to a modeling assumption in which the energy spectrum of the turbulent fluid is singularly dependent on the viscosity in the inertial range. Phenomenological models of turbulence and intermittency, by contrast, require that the energy spectrum be independent of the viscosity in the inertial range. When the energy spectrum is assumed to be consistent with the phenomenological models, there is no phase transition for turbulent diffusion.

The CuA center of cytochrome-c oxidase: electronic structure and spectra of models compared to the properties of CuA domains.

The electronic structure and spectrum of several models of the binuclear metal site in soluble CuA domains of cytochrome-c oxidase have been calculated by the use of an extended version of the complete neglect of differential overlap/spectroscopic method. The experimental spectra have two strong transitions of nearly equal intensity around 500 nm and a near-IR transition close to 800 nm. The model that best reproduces these features consists of a dimer of two blue (type 1) copper centers, in which each Cu atom replaces the missing imidazole on the other Cu atom. Thus, both Cu atoms have one cysteine sulfur atom and one imidazole nitrogen atom as ligands, and there are no bridging ligands but a direct Cu-Cu bond. According to the calculations, the two strong bands in the visible region originate from exciton coupling of the dipoles of the two copper monomers, and the near-IR band is a charge-transfer transition between the two Cu atoms. The known amino acid sequence has been used to construct a molecular model of the CuA site by the use of a template and energy minimization. In this model, the two ligand cysteine residues are in one turn of an alpha-helix, whereas one ligand histidine is in a loop following this helix and the other one is in a beta-strand.