Search

The Ca

The Ca 2+ -activated SK4 K + channel is gated by Ca 2+ –calmodulin (CaM) and is expressed in immune cells, brain, and heart. A cryoelectron microscopy (cryo-EM) structure of the human SK4 K + channel recently revealed four CaM molecules per channel tetramer, where the apo CaM C-lobe and the holo CaM N -lobe interact with the proximal carboxyl terminus and the linker S4–S5, respectively, to gate the channel. Here, we show that phosphatidylinositol 4-5 bisphosphate (PIP2) potently activates SK4 channels by docking to the boundary of the CaM-binding domain. An allosteric blocker, BA6b9, was designed to act to the CaM–PIP2-binding domain, a previously untargeted region of SK4 channels, at the interface of the proximal carboxyl terminus and the linker S4–S5. Site-directed mutagenesis, molecular docking, and patch-clamp electrophysiology indicate that BA6b9 inhibits SK4 channels by interacting with two specific residues, Arg191 and His192 in the linker S4–S5, not conserved in SK1–SK3 subunits, thereby conferring selectivity and preventing the Ca 2+ –CaM N -lobe from properly interacting with the channel linker region. Immunohistochemistry of the SK4 channel protein in rat hearts showed a widespread expression in the sarcolemma of atrial myocytes, with a sarcomeric striated Z-band pattern, and a weaker occurrence in the ventricle but a marked incidence at the intercalated discs. BA6b9 significantly prolonged atrial and atrioventricular effective refractory periods in rat isolated hearts and reduced atrial fibrillation induction ex vivo. Our work suggests that inhibition of SK4 K + channels by targeting drugs to the CaM–PIP2-binding domain provides a promising anti-arrhythmic therapy.

Human γD-crystallin protein is abundant in the lens and is essential for preserving lens transparency. With age the protein may lose its native structure resulting in the formation of cataract. We recently reported an aggregative peptide

The formation of amyloid-like structures by metabolites is associated with several inborn errors of metabolism (IEMs). These structures display most of the biological, chemical and physical properties of protein amyloids. However, the molecular interactions underlying the assembly remain elusive, and so far, no modulating therapeutic agents are available for clinical use. Chemical chaperones are known to inhibit protein and peptide amyloid formation and stabilize misfolded enzymes. Here, we provide an in-depth characterization of the inhibitory effect of osmolytes and hydrophobic chemical chaperones on metabolite assemblies, thus extending their functional repertoire. We applied a combined in vivo-in vitro-in silico approach and show their ability to inhibit metabolite amyloid-induced toxicity and reduce cellular amyloid content in yeast. We further used various biophysical techniques demonstrating direct inhibition of adenine self-assembly and alteration of fibril morphology by chemical chaperones. Using a scaffold-based approach, we analyzed the physiochemical properties of various dimethyl sulfoxide derivatives and their role in inhibiting metabolite self-assembly. Lastly, we employed whole-atom molecular dynamics simulations to elucidate the role of hydrogen bonds in osmolyte inhibition. Our results imply a dual mode of action of chemical chaperones as IEMs therapeutics, that could be implemented in the rational design of novel lead-like molecules.

International audience; Neurofibrillary tangles of the Tau protein and plaques of the amyloid β peptide are hallmarks of Alzheimer's disease (AD), which is characterized by the conversion of monomeric proteins/peptides into misfolded β-sheet rich fibrils. Halting the fibrillation process and disrupting the existing aggregates are key challenges for AD drug development. Previously, we performed in vitro high-throughput screening for the identification of potent inhibitors of Tau aggregation using a proxy model, a highly aggregation-prone hexapeptide fragment 306VQIVYK311 (termed PHF6) derived from Tau. Here we have characterized a hit molecule from that screen as a modulator of Tau aggregation using in vitro, in silico, and in vivo techniques. This molecule, an anthraquinone derivative named Purpurin, inhibited ~ 50% of PHF6 fibrillization in vitro at equimolar concentration and disassembled pre-formed PHF6 fibrils. In silico studies showed that Purpurin interacted with key residues of PHF6, which are responsible for maintaining its β-sheets conformation. Isothermal titration calorimetry and surface plasmon resonance experiments with PHF6 and full-length Tau (FL-Tau), respectively, indicated that Purpurin interacted with PHF6 predominantly via hydrophobic contacts and displayed a dose-dependent complexation with FL-Tau. Purpurin was non-toxic when fed to Drosophila and it significantly ameliorated the AD-related neurotoxic symptoms of transgenic flies expressing WT-FL human Tau (hTau) plausibly by inhibiting Tau accumulation and reducing Tau phosphorylation. Purpurin also reduced hTau accumulation in cell culture overexpressing hTau. Importantly, Purpurin efficiently crossed an in vitro human blood-brain barrier model. Our findings suggest that Purpurin could be a potential lead molecule for AD therapeutics.

This work employs Adult Polyglucosan Body Disease (APBD) models to explore the efficacy and mechanism of action of 144DG11, a new polyglucosan-reducing lead compound discovered by a high-throughput screen (HTS). APBD is an adult onset glycogen storage disorder (GSD) manifesting as a debilitating progressive axonopathic leukodystrophy. APBD is caused by glycogen branching enzyme (GBE) deficiency leading to poorly branched and insoluble glycogen inclusions, which precipitate as neuropathogenic polyglucosans (PG). 144DG11 led to prolonged survival and improved motor parameters in a GBE knockin (Gbeys/ys) APBD mouse model. Histopathologically, 144DG11 reduced PG and glycogen levels in brain, liver, heart, and peripheral nerve. Indirect calorimetry experiments revealed that 144DG11 increases carbohydrate burn at the expense of fat burn, suggesting metabolic mobilization of pathogenic PG. These results were also reflected at the cellular level by increased glycolytic, mitochondrial and total ATP production. Mechanistically, we show that the molecular target of 144DG11 is the lysosomal membrane protein LAMP1, whose interaction with the compound, similar to LAMP1 knockdown, enhanced autolysosomal degradation of glycogen and lysosomal acidification. Enhanced mitochondrial activity and lysosomal modifications were also the most pronounced effects of 144DG11 in APBD patient fibroblasts as discovered by image-based multiparametric phenotyping analysis and corroborated by proteomics. In summary, this work presents a broad mechanistic and target-based characterization of 144DG11 in in vivo and cell models of the prototypical GSD APBD. This investigation warrants development of 144DG11 into a safe and efficacious GSD therapy.One Sentence SummaryA new compound, demonstrated to ameliorate APBD in vivo and ex vivo by autophagic catabolism of glycogen, may potentially become a universal drug for glycogen storage disorders.

The Target of rapamycin (TOR) protein kinase forms part of TOR complex 1 (TORC1) and TOR complex 2 (TORC2), two multi-subunit protein complexes that regulate growth, proliferation, survival and developmental processes by phosphorylation and activation of AGC-family kinases. In the fission yeast, Schizosaccharomyces pombe, TORC2 and its target, the AGC kinase Gad8 (an orthologue of human AKT or SGK1) are required for viability under stress conditions and for developmental processes in response to starvation cues. In this study, we describe the isolation of gad8 mutant alleles that bypass the requirement for TORC2 and reveal a separation of function of TORC2 and Gad8 under stress conditions. In particular, osmotic and nutritional stress responses appear to form a separate branch from genotoxic stress responses downstream of TORC2-Gad8. Interestingly, TORC2-independent mutations map into the regulatory PIF pocket of Gad8, a highly conserved motif in AGC kinases that regulates substrate binding in PDK1 (phosphoinositide dependent kinase-1) and kinase activity in several AGC kinases. Gad8 activation is thought to require a two-step mechanism, in which phosphorylation by TORC2 allows further phosphorylation and activation by Ksg1 (an orthologue of PDK1). We focus on the Gad8-K263C mutation and demonstrate that it renders the Gad8 kinase activity independent of TORC2 in vitro and independent of the phosphorylation sites of TORC2 in vivo. Molecular dynamics simulations of Gad8-K263C revealed abnormal high flexibility at T387, the phosphorylation site for Ksg1, suggesting a mechanism for the TORC2-independent Gad8 activity. Significantly, the K263 residue is highly conserved in the family of AGC-kinases, which may suggest a general way of keeping their activity in check when acting downstream of TOR complexes., Author summary Protein kinases catalyze the transfer of phosphate from high-energy, phosphate-donating molecules, such as ATP, to their substrates. This process is pivotal for regulation of almost any aspect of cellular biology. Many human diseases are associated with aberrant functions of protein kinases due to mutations. Accordingly, there is a growing number of kinase inhibitors that have been approved for clinical use. A better understanding of how protein kinases become active and how their activity is relayed to regulate their cellular functions is much needed for rational design of kinase inhibitors and for their optimal use in the clinic. The AGC-family of protein kinases play key roles in regulating cellular growth, proliferation and survival. In human cells, as well as in the fission yeast, our cellular model system, a subgroup of the AGC kinases is activated by the TOR protein kinases. Here we report the isolation of mutations in the AGC kinase Gad8 (AKT or SGK1 in human) that bypass the requirement for activation by TOR. Analyses of how these mutations affect cellular growth revealed separate branches of stress response mechanisms downstream of Gad8, while computer simulation methods suggested a molecular mechanism that keeps the activity of Gad8 in check.

γD-crystallin is among the most abundant γ-crystallins in the human eye lens which are essential for preserving its transparency. Aging, and environmental changes, cause crystallins to lose their native soluble structure and aggregate, resulting in the formation of cataract. Current treatment of cataract is surgical removal which is costly. Pharmaceutical therapeutics of cataract is an unmet need. We report a screen for small molecules capable of inhibiting aggregation of human γD-crystallin. Using a highly amyloidogenic hexapeptide model 41GCWMLY46 derived from the full-length protein, we screened a library of 68 anthraquinone molecules using ThT fluorescence assay. A leading hit, the cochineal Carmine, effectively reduced aggregation of the model GDC6 peptide in dose dependent manner. Similar effect was observed toward aggregation of the full-length γD-crystallin. Transmission electron microscopy, intrinsic Tryptophan fluorescence and ANS fluorescence assays corroborated these results. Insights obtained from molecular docking suggested that Carmine interaction with monomeric GDC6 involved hydrogen bonding with Ace group, Cys, Met residues and hydrophobic contact with Trp residue. Carmine was non-toxic toward retinal cells in culture. It also reduced ex vivo the turbidity of human extracted cataract material. Collectively, our results indicate that Carmine could be used for developing new therapeutics to treat cataract.

One of the pathways of the unfolded protein response, initiated by PKR-like endoplasmic reticulum kinase (PERK), is key to neuronal homeostasis in neurodegenerative diseases. PERK pathway activation is usually accomplished by inhibiting eIF2α-P dephosphorylation, after its phosphorylation by PERK. Less tried is an approach involving direct PERK activation without compromising long-term recovery of eIF2α function by dephosphorylation. Here we show major improvement in cellular (STHdhQ111/111) and mouse (R6/2) Huntington's disease (HD) models using a potent small molecule PERK activator that we developed, MK-28. MK-28 showed PERK selectivity in vitro on a 391-kinase panel and rescued cells (but not PERK−/− cells) from ER stress-induced apoptosis. Cells were also rescued by the commercial PERK activator CCT020312 but MK-28 was significantly more potent. Computational docking suggested MK-28 interaction with the PERK activation loop. MK-28 exhibited remarkable pharmacokinetic properties and high BBB penetration in mice. Transient subcutaneous delivery of MK-28 significantly improved motor and executive functions and delayed death onset in R6/2 mice, showing no toxicity. Therefore, PERK activation can treat a most aggressive HD model, suggesting a possible approach for HD therapy and worth exploring for other neurodegenerative disorders.

Although aminoglycoside antibiotics are effective against Gram-negative infections, these drugs often cause irreversible hearing damage. Binding to the decoding site of the eukaryotic ribosomes appears to result in ototoxicity, but there is evidence that other effects are involved. Here, we show how chemical modifications of apramycin and geneticin, considered among the least and most toxic aminoglycosides, respectively, reduce auditory cell damage. Using molecular dynamics simulations, we studied how modified aminoglycosides influence the essential freedom of movement of the decoding site of the ribosome, the region targeted by aminoglycosides. By determining the ratio of a protein translated in mitochondria to that of a protein translated in the cytoplasm, we showed that aminoglycosides can paradoxically elevate rather than reduce protein levels. We showed that certain aminoglycosides induce rapid plasma membrane permeabilization and that this nonribosomal effect can also be reduced through chemical modifications. The results presented suggest a new paradigm for the development of safer aminoglycoside antibiotics.

Lithium diffusion in olivine phosphates has been widely studied both experimentally and theoretically. However, nuclear quantum effects (NQEs) of the Li ions have not been accounted for in theoretical studies thus far. In the current work, we compared Li and Na diffusion in Li0.25FePO4 and Na0.25FePO4 by computing density functional theory based classical diffusion barriers in conjunction with NQEs for the Li and Na ions. The NQEs are computed using a novel three-dimensional wave function method based on a path integral formulation. The calculations of both the potential and free energy diffusion barriers suggest that Li diffusion is faster than Na diffusion, in agreement with recent experiments. The NQEs for lithium ions in Li0.25FePO4 are higher than those for sodium ions in Na0.25FePO4. Although the contribution of NQEs to the computed Li and Na ion diffusion rates is rather small, the quantum behavior of the Li ions is unusual. Indeed, we observe a reduction in the computed diffusion rate for Li ions ...

A practical approach to treat nuclear quantum mechanical effects in simulations of condensed phases, such as enzymes, is via Feynman path integral (PI) formulations. Typically, the standard primitive approximation (PA) is employed in enzymatic PI simulations. Nonetheless, these PI simulations are computationally demanding due to the large number of beads required to obtain converged results. The efficiency of PI simulations may be greatly improved if higher-order factorizations of the density matrix operator are employed. Herein, we compare the results of model calculations obtained employing the standard PA, the improved operator of Takahashi and Imada (TI), and a gradient-based forward corrector algorithm due to Chin (CH). The quantum transmission coefficient is computed for the Eckart potential while the partition functions and rate constant are computed for the H 2 + H hydrogen transfer reaction. These potentials are simple models for chemical reactions. The study of the different factorization methods reveals that in most cases the higher-order action converges faster than the PA and TI approaches, at a moderate computational cost.

Background Alzheimer's disease (AD) is the most common neurodegenerative disorder which is characterized by the deposits of intra-cellular tau protein and extra-cellular amyloid-β (Aβ) peptides in the human brain. Understanding the mechanism of protein aggregation and finding compounds that are capable of inhibiting its aggregation is considered to be highly important for disease therapy. Methods We used an in vitro High-Throughput Screening for the identification of potent inhibitors of tau aggregation using a proxy model; a highly aggregation-prone hexapeptide fragment 306VQIVYK311 derived from tau. Using ThS fluorescence assay we screened a library of 2401 FDA approved, bio-active and natural compounds in attempt to find molecules which can efficiently modulate tau aggregation. Results Among the screened compounds, palmatine chloride (PC) alkaloid was able to dramatically reduce the aggregation propensity of PHF6 at sub-molar concentrations. PC was also able to disassemble preformed aggregates of PHF6 and reduce the amyloid content in a dose-dependent manner. Insights obtained from MD simulation showed that PC interacted with the key residues of PHF6 responsible for β-sheet formation, which could likely be the mechanism of inhibition and disassembly. Furthermore, PC could effectively inhibit the aggregation of full-length tau and disassemble preformed aggregates. Conclusions We found that PC possesses “dual functionality” towards PHF6 and full-length tau, i.e. inhibit their aggregation and disassemble pre-formed fibrils. General significance The “dual functionality” of PC is valuable as a disease modifying strategy for AD, and other tauopathies, by inhibiting their progress and reducing the effect of fibrils already present in the brain.

Enzymes are remarkably efficient catalysts evolved to perform well-defined and highly specific chemical transformations. Studying the nature of enzymatic rate enhancements is highly important from several aspects, including the rational design of synthetic catalysts and transition-state inhibitors. Herein, we describe recent progress in our group in the development of multiscale simulation methods and their application to several enzyme systems. In particular, we describe the use of combined quantum mechanics/molecular mechanics (QM/MM) methods in classical and quantum simulations. The development of various novel path-integral methods is reviewed. These methods are tailor-made for enzyme systems, where only a few degrees of freedom involved in the chemistry need to be quantized. The application of the hybrid QM/MM quantum-classical simulation approach to three case studies is presented. The first case involves proton transfer in nitroalkane oxidase, where the enzyme employs tunneling as a catalytic fine-tuning tool. The second case presented involves orotidine 5′-monophosphate decarboxylase, where multidimensional free energy simulations together with kinetic isotope effects are combined in the study of the reaction mechanism. Finally, we discuss the monoterpene cyclase bornyl diphosphate synthase, where non-statistical dynamics is a key component in enzyme function.

The inclusion of nuclear quantum effects such as zero-point energy and tunneling is of great importance in studying condensed phase chemical reactions involving the transfer of protons, hydrogen atoms, and hydride ions. In the current work, we derive an efficient quantum simulation approach for the computation of the momentum distribution in condensed phase chemical reactions. The method is based on a quantum-classical approach wherein quantum and classical simulations are performed separately. The classical simulations use standard sampling techniques, whereas the quantum simulations employ an open polymer chain path integral formulation which is computed using an efficient Monte Carlo staging algorithm. The approach is validated by applying it to a one-dimensional harmonic oscillator and symmetric double-well potential. Subsequently, the method is applied to the dihydrofolate reductase (DHFR) catalyzed reduction of 7,8-dihydrofolate by nicotinamide adenine dinucleotide phosphate hydride (NADPH) to yield S-5,6,7,8-tetrahydrofolate and NADP(+). The key chemical step in the catalytic cycle of DHFR involves a stereospecific hydride transfer. In order to estimate the amount of quantum delocalization, we compute the position and momentum distributions for the transferring hydride ion in the reactant state (RS) and transition state (TS) using a recently developed hybrid semiempirical quantum mechanics-molecular mechanics potential energy surface. Additionally, we examine the effect of compression of the donor-acceptor distance (DAD) in the TS on the momentum distribution. The present results suggest differential quantum delocalization in the RS and TS, as well as reduced tunneling upon DAD compression.

A practical approach to treat nuclear quantum mechanical (QM) effects in simulations of condensed phases, such as enzymes, is via Feynman path integral (PI) formulations. Typically, the standard primitive approximation (PA) is employed in enzymatic PI simulations. Nonetheless, these PI simulations are computationally demanding due to the large number of discretizations, or beads, required to obtain converged results. The efficiency of PI simulations may be greatly improved if higher order factorizations of the density matrix operator are employed. Herein, we compare the results of model calculations obtained employing the standard PA, the improved operator of Takahashi and Imada (TI), and several gradient-based forward corrector algorithms due to Chin (CH). The quantum partition function is computed for the harmonic oscillator, Morse, symmetric, and asymmetric double well potentials. These potentials are simple models for nuclear quantum effects, such as zero-point energy and tunneling. It is shown that a unique set of CH parameters may be employed for a variety of systems. Additionally, the nuclear QM effects of a water molecule, treated with density functional theory, are computed. Finally, we derive a practical perturbation expression for efficient computation of isotope effects in chemical systems using the staging algorithm. This new isotope effect approach is tested in conjunction with the PA, TI, and CH methods to compute the equilibrium isotope effect in the Schiff base-oxyanion keto-enol tautomerism in the cofactor pyridoxal-5'-phosphate in the enzyme alanine racemase. The study of the different factorization methods reveals that the higher-order actions converge substantially faster than the PA approach, at a moderate computational cost.

Formate dehydrogenase (FDH) catalyzes the oxidation of formic acid to carbon dioxide using nicotinamide adenine dinucleotide (NAD(+)) as a cofactor. In the current work we present extensive benchmark calculations for several model reactions in the gas phase that are relevant to the FDH catalyzed hydride transfer. To this end we employ G4MP2 and CBS-QB3 ab initio calculations as well as density functional theory methods. Using these results we develop a specific reaction parameter (SRP) Hamiltonian based on the semiempirical AM1 method. The SRP semiempirical Hamiltonian is subsequently used in hybrid quantum mechanics/molecular mechanics simulations of the FDH catalyzed reaction in Pseudomonas sp. 101 (PseFDH). The classical potential of mean force (PMF) is computed as a function of structural progress coordinates during the course of the hydride transfer reaction: The antisymmetric reactive stretch, the donor-acceptor distance, and an orbital rehybridization coordinate. The quantum PMF is computed using a centroid Feynman path-integral (PI) approach. Subsequently, kinetic isotope effects are computed using a mass-perturbation based PI method. Finally, the antisymmetric stretch vibrational frequency is computed for an azide ion in FDH and in aqueous solution.

Semiclassical integral representations, analogous to initial value expressions for the propagator, are presented for the Clebsch-Gordan angular momentum coupling coefficients. Two forms (L and R types) of the approximation are presented. For each form, new non-Gaussian expressions, which are specifically adapted to the nature of angular momentum variables, are proposed in place of the familiar Gaussian coherent state functions. With these non-Gaussian kernels, it is found that the present treatments are capable of accuracy similar to that obtained from a uniform Airy approximation. Although the present semiclassical approximations involve only real-valued angle variables, associated with sets of angular momenta that are related by ordinary, real, classical transformations, the treatments produce accurate results not only for classically allowed choices of quantum numbers but also for very strongly classically forbidden values.