The optimization principle in phylogenetic analysis tends to give incorrect topologies when the number of nucleotides or amino acids used is small

In the maximum parsimony (MP) and minimum evolution (ME) methods of phylogenetic inference, evolutionary trees are constructed by searching for the topology that shows the minimum number of mutational changes required (M) and the smallest sum of branch lengths (S), respectively, whereas in the maximum likelihood (ML) method the topology showing the highest maximum likelihood (A) of observing a given data set is chosen. However, the theoretical basis of the optimization principle remains unclear. We therefore examined the relationships of M, S, and A for the MP, ME, and ML trees with those for the true tree by using computer simulation. The results show that M and S are generally greater for the true tree than for the MP and ME trees when the number of nucleotides examined (n) is relatively small, whereas A is generally lower for the true tree than for the ML tree. This finding indicates that the optimization principle tends to give incorrect topologies when n is small. To deal with this disturbing property of the optimization principle, we suggest that more attention should be given to testing the statistical reliability of an estimated tree rather than to finding the optimal tree with excessive efforts. When a reliability test is conducted, simplified MP, ME, and ML algorithms such as the neighbor-joining method generally give conclusions about phylogenetic inference very similar to those obtained by the more extensive tree search algorithms.

The hard-soft acid-base principle in enzymatic catalysis: dual reactivity of phosphoenolpyruvate.

In this paper, the chemical reactivity of C3 of phosphoenolpyruvate (PEP) has been analyzed in terms of density functional theory quantified through quantum chemistry calculations. PEP is involved in a number of important enzymatic reactions, in which its C3 atom behaves like a base. In three different enzymatic reactions analyzed here, C3 sometimes behaves like a soft base and sometimes behaves like a hard base in terms of the hard-soft acid-base principle. This dual nature of C3 of PEP was found to be related to the conformational change of the molecule. This leads to a testable hypothesis: that PEP adopts particular conformations in the enzyme-substrate complexes of different PEP-using enzymes, and that the enzymes control the reactivity through controlling the dihedral angle between the carboxylate and the C==C double bond of PEP.

Game theory and reciprocity in some extensive form experimental games

We examine decision making in two-person extensive form game trees using nine treatments that vary matching protocol, payoffs, and payoff information. Our objective is to establish replicable principles of cooperative versus noncooperative behavior that involve the use of signaling, reciprocity, and backward induction strategies, depending on the availability of dominated direct punishing strategies and the probability of repeated interaction with the same partner. Contrary to the predictions of game theory, we find substantial support for cooperation under complete information even in various single-play treatments.

Polysomal ribonuclease 1 exists in a latent form on polysomes prior to estrogen activation of mRNA decay

Estrogen induces a global change in the translation profile of Xenopus hepatocytes, replacing serum protein synthesis with production of the yolk protein precursor vitellogenin. This is accomplished by the coordinate destabilization of serum protein mRNAs and the transcriptional induction and subsequent stabilization of vitellogenin mRNA. Previous work identified an endonuclease activity whose appearance on polysomes correlated with the disappearance of serum protein mRNAs. This enzyme, polysomal ribonuclease 1 (PMR1), is a novel member of the peroxidase gene family. The current study examined the association of PMR1 with its mRNA targets on polysomes and mRNPs. The highest amount of polysome-bound PMR1 was observed prior to estrogen induction of mRNA decay. Its distribution on sucrose density gradients matched the absorbance profile of polysome-bound mRNA, suggesting that PMR1 forms a latent complex with mRNA. Following dissociation with EDTA the 62 kDa PMR1 sedimented with a larger complex of >670 kDa. Estrogen induces a 22-fold increase in unit enzymatic activity of polysome-bound PMR1, and a time-dependent loss of PMR1 from polysomes in a manner that mirrors the disappearance of albumin mRNA. These data suggest that the key step in the extensive estrogen-induced change in mRNA decay in Xenopus liver is activation of a latent mRNA endonuclease associated with its target mRNA.

Site-specific mutations in the mature form of human IL-18 with enhanced biological activity and decreased neutralization by IL-18 binding protein

IL-18 can be considered a proinflammatory cytokine mediating disease as well as an immunostimulatory cytokine that is important for host defense against infection and cancer. The high-affinity, constitutively expressed, and circulating IL-18 binding protein (IL-18BP), which competes with cell surface receptors for IL-18 and neutralizes IL-18 activity, may act as a natural antiinflammatory as well as immunosuppressive molecule. In the present studies, the IL-18 precursor caspase-1 cleavage site was changed to a factor Xa site, and, after expression in Escherichia coli, mature IL-18 was generated by factor Xa cleavage. Mature IL-18 generated by factor Xa cleavage was fully active. Single point mutations in the mature IL-18 peptide were made, and the biological activities of the wild-type (WT) IL-18 were compared with those of the mutants. Mutants E42A and K89A exhibited 2-fold increased activity compared with WT IL-18. A double mutant, E42A plus K89A, exhibited 4-fold greater activity. Unexpectedly, IL-18BP failed to neutralize the double mutant E42A plus K89A compared with WT IL-18. The K89A mutant was intermediate in being neutralized by IL-18BP, whereas neutralization of the E42A mutant was comparable to that in the WT IL-18. The identification of E42 and K89 in the mature IL-18 peptide is consistent with previous modeling studies of IL-18 binding to IL-18BP and explains the unusually high affinity of IL-18BP for IL-18.

Plastome Engineering of Ribulose-1,5-Bisphosphate Carboxylase/Oxygenase in Tobacco to Form a Sunflower Large Subunit and Tobacco Small Subunit Hybrid1

Targeted gene replacement in plastids was used to explore whether the rbcL gene that codes for the large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase, the key enzyme of photosynthetic CO2 fixation, might be replaced with altered forms of the gene. Tobacco (Nicotiana tabacum) plants were transformed with plastid DNA that contained the rbcL gene from either sunflower (Helianthus annuus) or the cyanobacterium Synechococcus PCC6301, along with a selectable marker. Three stable lines of transformants were regenerated that had altered rbcL genes. Those containing the rbcL gene for cyanobacterial ribulose-1,5-bisphosphate carboxylase/oxygenase produced mRNA but no large subunit protein or enzyme activity. Those tobacco plants expressing the sunflower large subunit synthesized a catalytically active hybrid form of the enzyme composed of sunflower large subunits and tobacco small subunits. A third line expressed a chimeric sunflower/tobacco large subunit arising from homologous recombination within the rbcL gene that had properties similar to the hybrid enzyme. This study demonstrated the feasibility of using a binary system in which different forms of the rbcL gene are constructed in a bacterial host and then introduced into a vector for homologous recombination in transformed chloroplasts to produce an active, chimeric enzyme in vivo.

SoxR, a [2Fe-2S] transcription factor, is active only in its oxidized form.

SoxR protein is known to function both as a sensor and as a transcriptional activator for a superoxide response regulon in Escherichia coli. The activity of SoxR was tested by its ability to enable the transcription of its target gene, soxS, in vitro. The activity of the oxidized form was lost when its [2Fe-2S] clusters were reduced by dithionite under anaerobic conditions, and it was rapidly restored by autooxidation. This result is consistent with the hypothesis that induction of the regulon is effected by the univalent oxidation of the Fe-S centers of SoxR. In vivo, this oxidation may be caused by an alteration of the redox balance of electron chain intermediates that normally maintains soxR in an inactive, reduced state. Oxidized SoxR was about twice as effective as reduced SoxR in protecting the soxS operator from endonucleolytic cleavage. However, this difference could not account for a greater than 50-fold difference in their activities and therefore could not support a model in which oxidation activates SoxR by enabling it to bind to DNA. NADPH, ferredoxin, flavodoxin, or ferredoxin (flavodoxin):NADP+ reductase could not reduce SoxR directly in vitro at a measurable rate. The midpoint potential for SoxR was measured at -283 mV.

In its active form, the GTP-binding protein rab8 interacts with a stress-activated protein kinase.

Rab8 is a small GTP-binding protein that plays a role in vesicular transport from the trans-Golgi network to the basolateral plasma membrane in polarized epithelial cells (MDCK), and to the dendritic surface in hippocampal neurons. As is the case for most other rab proteins, the precise molecular interactions by which rab8 carries out its function remain to be elucidated. Here we report the identification and the complete cDNA-derived amino acid sequence of a murine rab8-interacting protein (rab8ip) that specifically interacts with rab8 in a GTP-dependent manner. Rab8ip displays 93% identity with the GC kinase, a serine/threonine protein kinase recently identified in human lymphoid tissue that is activated in the stress response. Like the GC kinase, rab8ip has protein kinase activity manifested by autophosphorylation and phosphorylation of the classical serine/threonine protein kinase substrates, myelin basic protein and casein. When coexpressed in transfected 293T cells, rab8 and the rab8ip/GC kinase formed a complex that could be recovered by immunoprecipitation with antibodies to rab8. Cell fractionation and immunofluorescence analyses indicate that in MDCK cells endogenous rab8ip is present both in the cytosol and as a peripheral membrane protein concentrated in the Golgi region and basolateral plasma membrane domains, sites where rab8 itself is also located. In light of recent evidence that rab proteins may act by promoting the stabilization of SNARE complexes, the specific GTP-dependent association of rab8 with the rab8ip/GC kinase raises the possibility that rab-regulated protein phosphorylation is important for vesicle targeting or fusion. Moreover, the rab8ip/GC kinase may serve to modulate secretion in response to stress stimuli.

Kinetic equilibrium of forces and molecular events in muscle contraction

In biomolecular systems, the mechanical transfer of free energy occurs with both high efficiency and high speed. It is shown here that such a transfer can be achieved only if the participating free-energy-storing elements exhibit opposing relationships between their content of free energy and the force they exert in the transfer direction. A kinetic equilibrium of forces (KEF) results, in which the transfer of free energy is mediated essentially by thermal molecular motion. On the basis of present evidence, KEF is used as a guiding principle in developing a mechanical model of the crossbridge cycle in muscle contraction. The model allows the basic features of molecular events to be visualized in terms of plausible structures. Real understanding of the process will require identification of the elements that perform the functions described here. Besides chemomechanical energy transduction, KEF may have a role in other biomolecular processes in which free energy is transferred mechanically over large distances.

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

Unraveling principles of lead discovery: from unfrustrated energy landscapes to novel molecular anchors.

The search for novel leads is a critical step in the drug discovery process. Computational approaches to identify new lead molecules have focused on discovering complete ligands by evaluating the binding affinity of a large number of candidates, a task of considerable complexity. A new computational method is introduced in this work based on the premise that the primary molecular recognition event in the protein binding site may be accomplished by small core fragments that serve as molecular anchors, providing a structurally stable platform that can be subsequently tailored into complete ligands. To fulfill its role, we show that an effective molecular anchor must meet both the thermodynamic requirement of relative energetic stability of a single binding mode and its consistent kinetic accessibility, which may be measured by the structural consensus of multiple docking simulations. From a large number of candidates, this technique is able to identify known core fragments responsible for primary recognition by the FK506 binding protein (FKBP-12), along with a diverse repertoire of novel molecular cores. By contrast, absolute energetic criteria for selecting molecular anchors are found to be promiscuous. A relationship between a minimum frustration principle of binding energy landscapes and receptor-specific molecular anchors in their role as "recognition nuclei" is established, thereby unraveling a mechanism of lead discovery and providing a practical route to receptor-biased computational combinatorial chemistry.

DNA strand breakage, activation of poly (ADP-ribose) synthetase, and cellular energy depletion are involved in the cytotoxicity of macrophages and smooth muscle cells exposed to peroxynitrite.

The free radicals nitric oxide and superoxide anion react to form peroxynitrite (ONOO-), a highly toxic oxidant species. In vivo formation of ONOO- has been demonstrated in shock and inflammation. Herein we provide evidence that cytotoxicity in cells exposed to ONOO- is mediated by DNA strand breakage and the subsequent activation of the DNA repair enzyme poly(ADP ribose) synthetase (PARS). Exposure to ONOO- (100 microM to 1 mM) inhibited mitochondrial respiration in cultured J774 macrophages and in rat aortic smooth muscle cells. The loss of cellular respiration was rapid, peaking 1-3 h after ONOO- exposure, and reversible, with recovery after a period of 6-24 h. The inhibition of mitochondrial respiration was paralleled by a dose-dependent increase in DNA strand breakage, reaching its maximum at 20-30 min after exposure to ONOO-. We observed a dose-dependent increase in the activity of PARS in cells exposed to ONOO-. Inhibitors of PARS such as 3-aminobenzamide (1 mM) prevented the inhibition of cellular respiration in cells exposed to ONOO-. Activation of PARS by ONOO--mediated DNA strand breakage resulted in a significant decrease in intracellular energy stores, as reflected by a decline of intracellular NAD+ and ATP content. 3-Aminobenzamide prevented the loss of NAD+ and ATP in cells exposed to ONOO-. In contrast, impairment of cellular respiration by the addition of the nitric oxide donors S-nitroso-N-acetyl-DL-penicillamine or diethyltriamine nitric oxide complex, was not associated with the development of DNA strand breaks, in concentrations up to 1 mM, and was largely refractory to PARS inhibition. Our results suggest that DNA damage and activation of PARS, an energy-consuming futile repair cycle, play a central role in ONOO--mediated cellular injury.

An alternative pathway for rod signals in the rodent retina: Rod photoreceptors, cone bipolar cells, and the localization of glutamate receptors

In the mammalian retina, extensive processing of spatiotemporal and chromatic information occurs. One key principle in signal transfer through the retina is parallel processing. Two of these parallel pathways are the ON- and OFF-channels transmitting light and dark signals. This dual system is created in the outer plexiform layer, the first relay station in retinal signal transfer. Photoreceptors release glutamate onto ON- and OFF-type bipolar cells, which are functionally distinguished by their postsynaptic expression of different types of glutamate receptors, namely ionotropic and metabotropic glutamate receptors. In the current concept, rod photoreceptors connect only to rod bipolar cells (ON-type) and cone photoreceptors connect only to cone bipolar cells (ON- and OFF-type). We have studied the distribution of (RS)-α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) glutamate receptor subunits at the synapses in the outer plexiform layer of the rodent retina by immunoelectron microscopy and serial section reconstruction. We report a non-classical synaptic contact and an alternative pathway for rod signals in the retina. Rod photoreceptors made synaptic contact with putative OFF-cone bipolar cells that expressed the AMPA glutamate receptor subunits GluR1 and GluR2 on their dendrites. Thus, in the retina of mouse and rat, an alternative pathway for rod signals exists, where rod photoreceptors bypass the rod bipolar cell and directly excite OFF-cone bipolar cells through an ionotropic sign-conserving AMPA glutamate receptor.

Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and 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.