Transient hydrophobic exposure in the molecular dynamics of Abeta peptide at low water concentration

Abeta is a disordered peptide central to Alzheimer's Disease. Aggregation of Abeta has been widely explored, but its molecular crowding less so. The synaptic cleft where Abeta locates only holds 60-70 water molecules along its width. We subjected Abeta40 to 100 different simulations with variable water cell size. We show that even for this disordered aggregation-prone peptide, many properties are not cell-size dependent, i.e. a small cell is easily justified. The radius of gyration, intra-peptide, and peptide-water hydrogen bonds are well-sampled by short (50 ns) time scales at any cell size. Abeta is mainly disordered with 0-30% alpha helix but undergoes consistent alpha-beta transitions up to 14% strand in 5-10% of the simulations regardless of cell size. The similar prevalence in long and short simulations indicate small diffusion barriers for structural transitions in contrast to folded globular proteins, which we suggest is a defining hallmark of intrinsically disordered proteins. Importantly, the hydrophobic surface increases significantly in small cells (confidence level 95%, two-tailed t-test), as does the variation in exposure and backbone conformations (>40% and >27% increased standard deviations). Whereas hydrophilic exposure dominates hydrophobic exposure in large cells, this tendency breaks down at low water concentration. We interpret these findings as a concentration-dependent hydrophobic effect, with the small water layer unable to keep the protein unexposed, an effect mainly caused by the layered water-water interactions, not by the peptide dynamics. The exposure correlates with radius of gyration (R2 0.35-0.50) and could be important in crowded environments, e.g. the synaptic cleft

The Influence of Selection for Protein Stability on dN/dS Estimations

Understanding the relative contributions of various evolutionary processes—purifying selection, neutral drift, and adaptation—is fundamental to evolutionary biology. A common metric to distinguish these processes is the ratio of nonsynonymous to synonymous substitutions (i.e., dN/dS) interpreted from the neutral theory as a null model. However, from biophysical considerations, mutations have non-negligible effects on the biophysical properties of proteins such as folding stability. In this work, we investigated how stability affects the rate of protein evolution in phylogenetic trees by using simulations that combine explicit protein sequences with associated stability changes. We first simulated myoglobin evolution in phylogenetic trees with a biophysically realistic approach that accounts for 3D structural information and estimates of changes in stability upon mutation. We then compared evolutionary rates inferred directly from simulation to those estimated using maximum-likelihood (ML) methods. We found that the dN/dS estimated by ML methods (ωML) is highly predictive of the per gene dN/dS inferred from the simulated phylogenetic trees. This agreement is strong in the regime of high stability where protein evolution is neutral. At low folding stabilities and under mutation-selection balance, we observe deviations from neutrality (per gene dN/dS > 1 and dN/dS 1. Altogether, we show how protein biophysics affects the dN/dS estimations and its subsequent interpretation. These results are important for improving the current approaches for detecting positive selection

Positively selected sites in cetacean myoglobins contribute to protein stability.

Since divergence ∼50 Ma ago from their terrestrial ancestors, cetaceans underwent a series of adaptations such as a ∼10–20 fold increase in myoglobin (Mb) concentration in skeletal muscle, critical for increasing oxygen storage capacity and prolonging dive time. Whereas the $O_2$-binding affinity of Mbs is not significantly different among mammals (with typical oxygenation constants of ∼0.8–1.2 $µM^{−1}$), folding stabilities of cetacean Mbs are ∼2–4 kcal/mol higher than for terrestrial Mbs. Using ancestral sequence reconstruction, maximum likelihood and Bayesian tests to describe the evolution of cetacean Mbs, and experimentally calibrated computation of stability effects of mutations, we observe accelerated evolution in cetaceans and identify seven positively selected sites in Mb. Overall, these sites contribute to Mb stabilization with a conditional probability of 0.8. We observe a correlation between Mb folding stability and protein abundance, suggesting that a selection pressure for stability acts proportionally to higher expression. We also identify a major divergence event leading to the common ancestor of whales, during which major stabilization occurred. Most of the positively selected sites that occur later act against other destabilizing mutations to maintain stability across the clade, except for the shallow divers, where late stability relaxation occurs, probably due to the shorter aerobic dive limits of these species. The three main positively selected sites 66, 5, and 35 undergo changes that favor hydrophobic folding, structural integrity, and intra-helical hydrogen bonds.Chemistry and Chemical Biolog

Bacterial porin disrupts mitochondrial membrane potential and sensitizes host cells to apoptosis

The bacterial PorB porin, an ATP-binding beta-barrel protein of pathogenic Neisseria gonorrhoeae, triggers host cell apoptosis by an unknown mechanism. PorB is targeted to and imported by host cell mitochondria, causing the breakdown of the mitochondrial membrane potential (delta psi m). Here, we show that PorB induces the condensation of the mitochondrial matrix and the loss of cristae structures, sensitizing cells to the induction of apoptosis via signaling pathways activated by BH3-only proteins. PorB is imported into mitochondria through the general translocase TOM but, unexpectedly, is not recognized by the SAM sorting machinery, usually required for the assembly of beta-barrel proteins in the mitochondrial outer membrane. PorB integrates into the mitochondrial inner membrane, leading to the breakdown of delta psi m. The PorB channel is regulated by nucleotides and an isogenic PorB mutant defective in ATP-binding failed to induce delta psi m loss and apoptosis, demonstrating that dissipation of delta psi m is a requirement for cell death caused by neisserial infection

Survival of the cheapest: How proteome cost minimization drives evolution

Darwin's theory of evolution emphasized that positive selection of functional proficiency provides the fitness that ultimately determines the structure of life, a view that has dominated biochemical thinking of enzymes as perfectly optimized for their specific functions. The 20th-century modern synthesis, structural biology, and the central dogma explained the machinery of evolution, and nearly neutral theory explained how selection competes with random fixation dynamics that produce molecular clocks essential e.g. for dating evolutionary histories. However, the quantitative proteomics revealed that fitness effects not related to functional proficiency play much larger roles on long evolutionary time scales than previously thought, with particular evidence that some universal biophysical selection pressures act via protein expression levels. This paper first summarizes recent progress in the 21st century towards recovering this universal selection pressure. Then, the paper argues that proteome cost minimization is the dominant, underlying "non-function" selection pressure controlling most of the evolution of already functionally adapted living systems. A theory of proteome cost minimization is described and argued to have consequences for understanding evolutionary trade-offs, aging, cancer, and neurodegenerative protein-misfolding diseases

Bim and Bmf synergize to induce apoptosis in Neisseria gonorrhoeae infection

Abstract: Bcl-2 family proteins including the pro-apoptotic BH3-only proteins are central regulators of apoptotic cell death. Here we show by a focused siRNA miniscreen that the synergistic action of the BH3-only proteins Bim and Bmf is required for apoptosis induced by infection with Neisseria gonorrhoeae (Ngo). While Bim and Bmf were associated with the cytoskeleton of healthy cells, they both were released upon Ngo infection. Loss of Bim and Bmf from the cytoskeleton fraction required the activation of Jun-N-terminal kinase-1 (JNK-1), which in turn depended on Rac-1. Depletion and inhibition of Rac-1, JNK-1, Bim, or Bmf prevented the activation of Bak and Bax and the subsequent activation of caspases. Apoptosis could be reconstituted in Bim-depleted and Bmf-depleted cells by additional silencing of antiapoptotic Mcl-1 and Bcl-XL, respectively. Our data indicate a synergistic role for both cytoskeletal-associated BH3-only proteins, Bim, and Bmf, in an apoptotic pathway leading to the clearance of Ngo-infected cells. Author Summary: A variety of physiological death signals, as well as pathological insults, trigger apoptosis, a genetically programmed form of cell death. Pathogens often induce host cell apoptosis to establish a successful infection. Neisseria gonorrhoeae (Ngo), the etiological agent of the sexually transmitted disease gonorrhoea, is a highly adapted obligate human-specific pathogen and has been shown to induce apoptosis in infected cells. Here we unveil the molecular mechanisms leading to apoptosis of infected cells. We show that Ngo-mediated apoptosis requires a special subset of proapoptotic proteins from the group of BH3-only proteins. BH3-only proteins act as stress sensors to translate toxic environmental signals to the initiation of apoptosis. In a siRNA-based miniscreen, we found Bim and Bmf, BH3-only proteins associated with the cytoskeleton, necessary to induce host cell apoptosis upon infection. Bim and Bmf inactivated different inhibitors of apoptosis and thereby induced cell death in response to infection. Our data unveil a novel pathway of infection-induced apoptosis that enhances our understanding of the mechanism by which BH3-only proteins control apoptotic cell death

Exploring Spirituality in Teaching Within a Christian School Context Through Collaborative Action Research

This article reports on a collaborative action research project conducted in New Zealand, during 2012, exploring spirituality in teaching within a Christian school context. The experienced primary school teacher participant chose to take action around the issue of personal fear and insecurity which were believed to be hindering professional growth and relationships. Through self-directed inquiry, critical reflective journaling, Bible study, fellowship and prayer with trusted friends, the teacher experienced a renewed sense of peace and freedom in Christ. This personal transformation was believed to be influential on subsequent professional practice, assisting the teacher to become more relational, responsive and compassionate. The findings provide a rich description of the participant’s spirituality, the lived reality of a person’s spiritual life. This report will be of interest to teachers, teacher-leaders and teacher-educators who desire to explore Christian spirituality through practitioner-led inquiry

Allele frequencies of variants in Ultra Conserved Elements identify selective pressure on transcription factor binding

Ultra-conserved genes or elements (UCGs/UCEs) in the human genome are extreme examples of conservation. We characterized natural variations in 2884 UCEs and UCGs in two distinct populations ; Singaporean Chinese (n=280) and Italian (n=501) by using a pooled sample, targeted capture, sequencing approach. We identify, with high confidence, in these regions the abundance of rare SNVs (MAF<0.5%) of which 75% is not present in dbSNP137. UCEs association studies for complex human traits can use this information to model expected background variation and thus necessary power for association studies. By combining our data with 1000 Genome Project data, we show in three independent datasets that prevalent UCE variants (MAF>5%) are more often found in relatively less-conserved nucleotides within UCEs, compared to rare variants. Moreover, prevalent variants are less likely to overlap transcription factor binding site. Using SNPfold we found no significant influence of RNA secondary structure on UCE conservation. All together, these results suggest UCEs are not under selective pressure as a stretch of DNA but are under differential evolutionary pressure on the single nucleotide level

Using Electronegativity and Hardness to Test Density Functional Universality

Density functional theory (DFT) is used in thousands of papers each year, yet lack of universality reduces DFT's predictive capacity, and functionals may produce energy-density imbalances. The absolute electronegativity (\chi) and hardness (\eta) directly reflect the energy-density relationship via the chemical potential dE/dN and we thus hypothesized that they probe universality. We studied \chi and \eta for atoms Z = 1-36 using 50 diverse functionals covering all major classes. Very few functionals describe both \chi and \eta well. \eta benefits from error cancelation whereas \chi is marred by error propagation from IP and EA; thus almost all standard GGA and hybrid functionals display a plateau in the MAE at 0.2-0.3 eV for \eta. In contrast, variable performance for \chi indicates problems in describing the chemical potential by DFT. The accuracy and precision of a functional is far from linearly related, yet for a universal functional we expect linearity. Popular functionals such as B3LYP, PBE, and revPBE, perform poorly for both properties. Density sensitivity calculations indicate large density-derived errors as occupation of degenerate p- and d-orbitals causes "non-universality" and large dependency on exact exchange. Thus, we argue that performance for \chi for the same systems is a hallmark of universality by probing dE/dN. With this metric, B98, B97-1, PW6B95D3, APFD are the most "universal" tested functionals. B98 and B97-1 are accurate for very diverse metal-ligand bonds, supporting that a balanced description of dE/dN and dE2/dN2, via \chi and \eta, is probably a first simple probe of universality