Symmetry‐restrained molecular dynamics simulations improve homology models of potassium channels

Most crystallized homo‐oligomeric ion channels are highly symmetric, which dramatically decreases conformational space and facilitates building homology models (HMs). However, in molecular dynamics (MD) simulations channels deviate from ideal symmetry and accumulate thermal defects, which complicate the refinement of HMs using MD. In this work we evaluate the ability of symmetry constrained MD simulations to improve HMs accuracy, using an approach conceptually similar to Critical Assessment of techniques for protein Structure Prediction (CASP) competition: build HMs of channels with known structure and evaluate the efficiency of proposed methods in improving HMs accuracy (measured as deviation from experimental structure). Results indicate that unrestrained MD does not improve the accuracy of HMs, instantaneous symmetrization improves accuracy but not stability of HMs during subsequent unrestrained MD, while gradually imposing symmetry constraints improves both accuracy (by 5–50%) and stability of HMs. Moreover, accuracy and stability are strongly correlated, making stability a reliable criterion in predicting the accuracy of new HMs. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Association Thermodynamics and Conformational Stability of β-Sheet Amyloid β(17-42) Oligomers: Effects of E22Q (Dutch) Mutation and Charge Neutralization

Amyloid fibrils are associated with many neurodegenerative diseases. It was found that amyloidogenic oligomers, not mature fibrils, are neurotoxic agents related to these diseases. Molecular mechanisms of infectivity, pathways of aggregation, and molecular structure of these oligomers remain elusive. Here, we use all-atom molecular dynamics, molecular mechanics combined with solvation analysis by statistical-mechanical, three-dimensional molecular theory of solvation (also known as 3D-RISM-KH) in a new MM-3D-RISM-KH method to study conformational stability, and association thermodynamics of small wild-type Aβ_17-42 oligomers with different protonation states of Glu²², as well the E22Q (Dutch) mutants. The association free energy of small β-sheet oligomers shows near-linear trend with the dimers being thermodynamically more stable relative to the larger constructs. The linear (within statistical uncertainty) dependence of the association free energy on complex size is a consequence of the unilateral stacking of monomers in the β-sheet oligomers. The charge reduction of the wild-type Aβ_17-42 oligomers upon protonation of the solvent-exposed Glu²² at acidic conditions results in lowering the association free energy compared to the wild-type oligomers at neutral pH and the E22Q mutants. The neutralization of the peptides because of the E22Q mutation only marginally affects the association free energy, with the reduction of the direct electrostatic interactions mostly compensated by the unfavorable electrostatic solvation effects. For the wild-type oligomers at acidic conditions such compensation is not complete, and the electrostatic interactions, along with the gas-phase nonpolar energetic and the overall entropic effects, contribute to the lowering of the association free energy. The differences in the association thermodynamics between the wild-type Aβ_17-42 oligomers at neutral pH and the Dutch mutants, on the one hand, and the Aβ_17-42 oligomers with protonated Glu²², on the other, may be explained by destabilization of the inter- and intrapeptide salt bridges between Asp²³ and Lys²⁸. Peculiarities in the conformational stability and the association thermodynamics for the different models of the Aβ_17-42 oligomers are rationalized based on the analysis of the local physical interactions and the microscopic solvation structure.     

Ionizing radiation and hematopoietic malignancies: Altering the adaptive landscape

Somatic evolution, which underlies tumor progression, is driven by two essential components: (1) diversification of phenotypes through heritable mutations and epigenetic changes and (2) selection for mutant clones which possess higher fitness. Exposure to ionizing radiation (IR) is highly associated with increased risk of carcinogenesis. This link is traditionally attributed to causation of oncogenic mutations through the mutagenic effects of irradiation. On the other hand, potential effects of irradiation on altering fitness and increasing selection for mutant clones are frequently ignored. Recent studies bring the effects of irradiation on fitness and selection into focus, demonstrating that IR exposure results in stable reductions in the fitness of hematopoietic stem and progenitor cell populations. These reductions of fitness are associated with alteration of the adaptive landscape, increasing the selective advantages conferred by certain oncogenic mutations. Therefore, the link between irradiation and carcinogenesis might be more complex than traditionally appreciated: while mutagenic effects of irradiation should increase the probability of occurrence of oncogenic mutations, IR can also work as a tumor promoter, increasing the selective expansion of clones bearing mutations which become advantageous in the irradiation-altered environment, such as activated mutations in Notch1 or disrupting mutations in p53.Key words: Notch, p53, fitness, irradiation, hematopoietic, evolution     

Dissecting muscle and neuronal disorders in a Drosophila model of muscular dystrophy

Perturbation in the Dystroglycan (Dg)-Dystrophin (Dys) complex results in muscular dystrophies and brain abnormalities in human. Here we report that Drosophila is an excellent genetically tractable model to study muscular dystrophies and neuronal abnormalities caused by defects in this complex. Using a fluorescence polarization assay, we show a high conservation in Dg-Dys interaction between human and Drosophila. Genetic and RNAi-induced perturbations of Dg and Dys in Drosophila cause cell polarity and muscular dystrophy phenotypes: decreased mobility, age-dependent muscle degeneration and defective photoreceptor path-finding. Dg and Dys are required in targeting glial cells and neurons for correct neuronal migration. Importantly, we now report that Dg interacts with insulin receptor and Nck/Dock SH2/SH3-adaptor molecule in photoreceptor path-finding. This is the first demonstration of a genetic interaction between Dg and InR.     

The requirement of sterol glucoside for pexophagy in yeast is dependent on the species and nature of peroxisome inducers

Sterol glucosyltransferase, Ugt51/Atg26, is essential for both micropexophagy and macropexophagy of methanol-induced peroxisomes in Pichia pastoris. However, the role of this protein in pexophagy in other yeast remained unclear. We show that oleate- and amine-induced peroxisomes in Yarrowia lipolytica are degraded by Atg26-independent macropexophagy. Surprisingly, Atg26 was also not essential for macropexophagy of oleate- and amine-induced peroxisomes in P. pastoris, suggesting that the function of sterol glucoside (SG) in pexophagy is both species and peroxisome inducer specific. However, the rates of degradation of oleate- and amine-induced peroxisomes in P. pastoris were reduced in the absence of SG, indicating that P. pastoris specifically uses sterol conversion by Atg26 to enhance selective degradation of peroxisomes. However, methanol-induced peroxisomes apparently have lost the redundant ability to be degraded without SG. We also show that the P. pastoris Vac8 armadillo repeat protein is not essential for macropexophagy of methanol-, oleate-, or amine-induced peroxisomes, which makes PpVac8 the first known protein required for the micropexophagy, but not for the macropexophagy, machinery. The uniqueness of Atg26 and Vac8 functions under different pexophagy conditions demonstrates that not only pexophagy inducers, such as glucose or ethanol, but also the inducers of peroxisomes, such as methanol, oleate, or primary amines, determine the requirements for subsequent pexophagy in yeast.     

Molecular Biophysics of Orai Store-Operated Ca2+ Channels

Upon endoplasmic reticulum Ca²⁺ store depletion, Orai channels in the plasma membrane are activated directly by endoplasmic reticulum-resident STIM proteins to generate the Ca²⁺-selective, Ca²⁺ release-activated Ca²⁺ (CRAC) current. After the molecular identification of Orai, a plethora of functional and biochemical studies sought to compare Orai homologs, determine their stoichiometry, identify structural domains responsible for the biophysical fingerprint of the CRAC current, identify the physiological functions, and investigate Orai homologs as potential therapeutic targets. Subsequently, the solved crystal structure of Drosophila Orai (dOrai) substantiated many findings from structure-function studies, but also revealed an unexpected hexameric structure. In this review, we explore Orai channels as elucidated by functional and biochemical studies, analyze the dOrai crystal structure and its implications for Orai channel function, and present newly available information from molecular dynamics simulations that shed light on Orai channel gating and permeation.     

CCZ1, MON1 and YPT7 genes are involved in pexophagy, the Cvt pathway and non-specific macroautophagy in the methylotrophic yeast Pichia pastoris

Orthologues of Saccharomyces cerevisiae CCZ1, MON1 and YPT7 genes in the methylotrophic yeast, Pichia pastoris, have been identified. These genes encode proteins, which act as a complex, being involved in degradation of oleate-induced peroxisomes, Cvt (cytoplasm to vacuole targeting) pathway and non-specific macroautophagy in S. cerevisiae. CCZ1, MON1 and YPT7 gene orthologues are essential for multiple delivery pathways in P. pastoris. Strains with deletion of either of these genes displayed complete deficiency in pexophagy, non-specific macroautophagy and the biosynthetic Cvt pathway. The data suggest that CCZ1, MON1 and YPT7 genes are involved in degradation of both small oleate-induced and large methanol-induced peroxisomes. The data suggest conservative functions of CCZ1, MON1 and YPT7 genes among yeast species.     

Structural models of TREK channels and their gating mechanism

Mechanosensitive TREK channels belong to the family of K2P channels, a family of widely distributed, well modulated channels that uniquely have two similar or identical subunits, each with two TM1-P-TM2 motifs. Our goal is to build viable structural models of TREK channels, as representatives of K2P channels family. The structures available to be used as templates belong to the 2TM channels superfamily. These have low sequence similarity and different structural features: four symmetrically arranged subunits, each having one TM1-P-TM2 motif. Our model building strategy used two subunits of the template (KcsA) to build one subunit of the target (TREK-1). Our models of the Closed channel were adjusted to differ substantially from those of the template, e.g., TM2 of the 2nd repeat is near the axis of the pore whereas TM2 of the 1st repeat is far from the axis. Segments linking the two repeats and immediately following the last TM segment were modeled ab initio as α-helices based on helical periodicities of hydrophobic and hydrophilic residues, highly conserved and poorly conserved residues, and statistically related positions from multiple sequence alignments. The models were further refined by two-fold symmetry-constrained MD simulations using a protocol we developed previously. We also built models of the Open state and suggest a possible tension-activated gating mechanism characterized by helical motion with two-fold symmetry. Our models are consistent with deletion/truncation mutagenesis and thermodynamic analysis of gating described in the accompanying paper.