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].

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.