Insight into Early-Stage Unfolding of GPI-Anchored Human Prion Protein

Protein-lipid interactions and protein anchoring modulate the modes of association of the globular domain of the Prion protein and Doppel protein to model membrane patches

The Prion protein is the molecular hallmark of the incurable prion diseases affecting mammals, including humans. The protein-only hypothesis states that the misfolding, accumulation, and deposition of the Prion protein play a critical role in toxicity. The cellular Prion protein (PrPC) anchors to the extracellular leaflet of the plasma membrane and prefers cholesterol- and sphingomyelin-rich membrane domains. Conformational Prion protein conversion into the pathological isoform happens on the cell surface. In vitro and in vivo experiments indicate that Prion protein misfolding, aggregation, and toxicity are sensitive to the lipid composition of plasma membranes and vesicles. A picture of the underlying biophysical driving forces that explain the effect of Prion protein - lipid interactions in physiological conditions is needed to develop a structural model of Prion protein conformational conversion. To this end, we use molecular dynamics simulations that mimic the interactions between the globular domain of PrPC anchored to model membrane patches. In addition, we also simulate the Doppel protein anchored to such membrane patches. The Doppel protein is the closest in the phylogenetic tree to PrPC, localizes in an extracellular milieu similar to that of PrPC, and exhibits a similar topology to PrPC even if the amino acid sequence is only 25% identical. Our simulations show that specific protein-lipid interactions and conformational constraints imposed by GPI anchoring together favor specific binding sites in globular PrPC but not in Doppel. Interestingly, the binding sites we found in PrPC correspond to prion protein loops, which are critical in aggregation and prion disease transmission barrier (β2-α2 loop) and in initial spontaneous misfolding (α2-α3 loop). We also found that the membrane re-arranges locally to accommodate protein residues inserted in the membrane surface as a response to protein binding.

Structural investigation of human cystine/glutamate antiporter system xc− (Sxc−) using homology modeling and molecular dynamics

The cystine/glutamate antiporter system xc− (Sxc−) belongs to the SLC7 family of plasma membrane transporters. It exports intracellular glutamate along the latter’s concentration gradient as a driving force for cellular uptake of cystine. Once imported, cystine is mainly used for the production of glutathione, a tripeptide thiol crucial in maintenance of redox homeostasis and protection of cells against oxidative stress. Overexpression of Sxc− has been found in several cancer cells, where it is thought to counteract the increased oxidative stress. In addition, Sxc− is important in the central nervous system, playing a complex role in regulating glutamatergic neurotransmission and glutamate toxicity. Accordingly, this transporter is considered a potential target for the treatment of cancer as well as neurodegenerative diseases. Till now, no specific inhibitors are available. We herein present four conformations of Sxc− along its transport pathway, obtained using multi-template homology modeling and refined by means of Molecular Dynamics. Comparison with a very recently released cryo-EM structure revealed an excellent agreement with our inward-open conformation. Intriguingly, our models contain a structured N-terminal domain that is unresolved in the experimental structures and is thought to play a gating role in the transport mechanism of other SLC7 family members. In contrast to the inward-open model, there is no direct experimental counterpart for the other three conformations we obtained, although they are in fair agreement with the other stages of the transport mechanism seen in other SLC7 transporters. Therefore, our models open the prospect for targeting alternative Sxc− conformations in structure-based drug design efforts.

The importance of side branches of glycosylphosphatidylinositol anchors: a molecular dynamics perspective

Many proteins are anchored to the cell surface of eukaryotes using a unique family of glycolipids called glycosylphosphatidylinositol (GPI) anchors. These glycolipids also exist without a covalently bound protein, in particular on the cell surfaces of protozoan parasites where they are densely populated. GPIs and GPI-anchored proteins participate in multiple cellular processes such as signal transduction, cell adhesion, protein trafficking and pathogenesis of Malaria, Toxoplasmosis, Trypanosomiasis and prion diseases, among others. All GPIs share a common conserved glycan core modified in a cell-dependent manner with additional side glycans or phosphoethanolamine residues. Here, we use atomistic molecular dynamic simulations and perform a systematic study to evaluate the structural properties of GPIs with different side chains inserted in lipid bilayers. Our results show a flop-down orientation of GPIs with respect to the membrane surface and the presentation of the side chain residues to the solvent. This finding agrees well with experiments showing the role of the side residues as active epitopes for recognition of GPIs by macrophages and induction of GPI-glycan-specific immune responses. Protein-GPI interactions were investigated by attaching parasitic GPIs to Green Fluorescent Protein. GPIs are observed to recline on the membrane surface and pull down the attached protein close to the membrane facilitating mutual contacts between protein, GPI and the lipid bilayer. This model is efficient in evaluating the interaction of GPIs and GPI-anchored proteins with membranes and can be extended to study other parasitic GPIs and proteins and develop GPI-based immunoprophylaxis to treat infectious diseases.

Challenges and limitations in the studies of glycoproteins: A computational chemist's perspective

Experimenters face challenges and limitations while analyzing glycoproteins due to their high flexibility, stereochemistry, anisotropic effects, and hydration phenomena. Computational studies complement experiments and have been used in characterization of the structural properties of glycoproteins. However, recent investigations revealed that computational studies face significant challenges as well. Here, we introduce and discuss some of these challenges and weaknesses in the investigations of glycoproteins. We also present requirements of future developments in computational biochemistry and computational biology areas that could be necessary for providing more accurate structural property analyses of glycoproteins using computational tools. Further theoretical strategies that need to be and can be developed are discussed herein.

Molecular foundations of prion strain diversity

Despite being caused by a single protein, prion diseases are strikingly heterogenous. Individual prion variants, known as strains, possess distinct biochemical properties, form aggregates with characteristic morphologies and preferentially seed certain brain regions, causing markedly different disease phenotypes. Strain diversity is determined by protein structure, post-translational modifications and the presence of extracellular matrix components, with single amino acid substitutions or altered protein glycosylation exerting dramatic effects. Here, we review recent advances in the study of prion strains and discuss how a deeper knowledge of the molecular origins of strain heterogeneity is providing a foundation for the development of anti-prion therapeutics.

Prion protein glycans reduce intracerebral fibril formation and spongiosis in prion disease

Posttranslational modifications (PTMs) are common among proteins that aggregate in neurodegenerative disease, yet how PTMs impact the aggregate conformation and disease progression remains unclear. By engineering knockin mice expressing prion protein (PrP) lacking 2 N-linked glycans (Prnp180Q/196Q), we provide evidence that glycans reduce spongiform degeneration and hinder plaque formation in prion disease. Prnp180Q/196Q mice challenged with 2 subfibrillar, non-plaque-forming prion strains instead developed plaques highly enriched in ADAM10-cleaved PrP and heparan sulfate (HS). Intriguingly, a third strain composed of intact, glycophosphatidylinositol-anchored (GPI-anchored) PrP was relatively unchanged, forming diffuse, HS-deficient deposits in both the Prnp180Q/196Q and WT mice, underscoring the pivotal role of the GPI-anchor in driving the aggregate conformation and disease phenotype. Finally, knockin mice expressing triglycosylated PrP (Prnp187N) challenged with a plaque-forming prion strain showed a phenotype reversal, with a striking disease acceleration and switch from plaques to predominantly diffuse, subfibrillar deposits. Our findings suggest that the dominance of subfibrillar aggregates in prion disease is due to the replication of GPI-anchored prions, with fibrillar plaques forming from poorly glycosylated, GPI-anchorless prions that interact with extracellular HS. These studies provide insight into how PTMs impact PrP interactions with polyanionic cofactors, and highlight PTMs as a major force driving the prion disease phenotype.

Characteristics of a Folate Receptor-α Anchored into a Multilipid Bilayer Obtained from Atomistic Molecular Dynamics Simulations

Thorough computational description of the properties of membrane-anchored protein receptors, which are important for example in the context of active targeting drug delivery, may be achieved by models representing as close as possible the immediate environment of these macromolecules. An all-atom bilayer, including 35 different lipid types asymmetrically distributed among the two monolayers, is suggested as a model neoplastic cell membrane. One molecule of folate receptor-α (FRα) is anchored into its outer leaflet, and the behavior of the system is explored by atomistic molecular dynamics simulations. The total number of atoms in the model is ∼185 000. Three 1-μs-long simulations are carried out, where physiological conditions (310 K and 1 bar) are maintained with three different pressure scaling schemes. To evaluate the structure and the phase state of the membrane, the density profiles of the system, the average area per lipid, and the deuterium order parameter of the lipid tails are calculated. The bilayer is in liquid ordered state, and the specific arrangement varies between the three trajectories. The changes in the structure of FRα are investigated and are found time- and ensemble-dependent. The volume of the ligand binding pocket fluctuates with time, but this variation remains independent of the more global structural alterations. The latter are mostly "waving" motions of the protein, which periodically approaches and retreats from the membrane. The semi-isotropic pressure scaling perturbs the receptor most significantly, while the isotropic algorithm induces rather slow changes. Maintaining constant nonzero surface tension leads to behavior closest to the experimentally observed one.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

CHARMM-GUI Membrane Builder for Complex Biological Membrane Simulations with Glycolipids and Lipoglycans

Glycolipids (such as glycoglycerolipids, glycosphingolipids, and glycosylphosphatidylinositol) and lipoglycans (such as lipopolysaccharides (LPS), lipooligosaccharides (LOS), mycobacterial lipoarabinomannan, and mycoplasma lipoglycans) are typically found on the surface of cell membranes and play crucial roles in various cellular functions. Characterizing their structure and dynamics at the molecular level is essential to understand their biological roles, but systematic generation of glycolipid and lipoglycan structures is challenging because of great variations in lipid structures and glycan sequences (i.e., carbohydrate types and their linkages). To facilitate the generation of all-atom glycolipid/LPS/LOS structures, we have developed Glycolipid Modeler and LPS Modeler in CHARMM-GUI ( http://www.charmm-gui.org ), a web-based interface that simplifies building of complex biological simulation systems. In addition, we have incorporated these modules into Membrane Builder so that users can readily build a complex symmetric or asymmetric biological membrane system with various glycolipids and LPS/LOS. These tools are expected to be useful in innovative and novel glycolipid/LPS/LOS modeling and simulation research by easing tedious and intricate steps in modeling complex biological systems and shall provide insight into structures, dynamics, and underlying mechanisms of complex glycolipid-/LPS-/LOS-containing biological membrane systems.

Glycosylation Significantly Inhibits the Aggregation of Human Prion Protein and Decreases Its Cytotoxicity

Prion diseases are primarily caused by the misfolding of prion proteins in humans, cattle, sheep, and cervid species. The effects of glycosylation on prion protein (PrP) structure and function have not been thoroughly elucidated to date. In this study, we attempt to elucidate the effects of glycosylation on the aggregation and toxicity of human PrP. As revealed by immunocytochemical staining, wild-type PrP and its monoglycosylated mutants N181D, N197D, and T199N/N181D/N197D are primarily attached to the plasma membrane. In contrast, PrP F198S, a pathological mutant with an altered residue within the glycosylation site, and an unglycosylated PrP mutant, N181D/N197D, primarily exist in the cytoplasm. In the pathological mutant V180I, there is an equal mix of membranous and cytoplasmic PrP, indicating that N-linked glycosylation deficiency impairs the correct localization of human PrP at the plasma membrane. As shown by immunoblotting and flow cytometry, human PrP located in the cytoplasm displays considerably greater PK resistance and aggregation ability and is associated with considerably higher cellular ROS levels than PrP located on the plasma membrane. Furthermore, glycosylation deficiency enhances human PrP cytotoxicity induced by MG132 or the toxic prion peptide PrP 106-126. Therefore, we propose that glycosylation acts as a necessary cofactor in determining PrP localization on the plasma membrane and that it significantly inhibits the aggregation of human PrP and decreases its cytotoxicity.

Role of the cell membrane interface in modulating production and uptake of Alzheimer's beta amyloid protein

The beta amyloid protein (Aβ) plays a central role in Alzheimer's disease (AD) pathogenesis and its interaction with cell membranes in known to promote mutually disruptive structural perturbations that contribute to amyloid deposition and neurodegeneration in the brain. In addition to protein aggregation at the membrane interface and disruption of membrane integrity, growing reports demonstrate an important role for the membrane in modulating Aβ production and uptake into cells. The aim of this review is to highlight and summarize recent literature that have contributed insight into the implications of altered membrane composition on amyloid precursor protein (APP) proteolysis, production of Aβ, its internalization in to cells via permeabilization and receptor mediated uptake. Here, we also review the various membrane model systems and experimental tools used for probing Aβ-membrane interactions to investigate the key mechanistic aspects underlying the accumulation and toxicity of Aβ in AD.

CHARMM-GUI 10 years for biomolecular modeling and simulation

CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM-GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1) PDB Reader & Manipulator, Glycan Reader, and Ligand Reader & Modeler for reading and modifying molecules; (2) Quick MD Simulator, Membrane Builder, Nanodisc Builder, HMMM Builder, Monolayer Builder, Micelle Builder, and Hex Phase Builder for building all-atom simulation systems in various environments; (3) PACE CG Builder and Martini Maker for building coarse-grained simulation systems; (4) DEER Facilitator and MDFF/xMDFF Utilizer for experimentally guided simulations; (5) Implicit Solvent Modeler, PBEQ-Solver, and GCMC/BD Ion Simulator for implicit solvent related calculations; (6) Ligand Binder for ligand solvation and binding free energy simulations; and (7) Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM-GUI, such as Glycolipid Modeler for generation of various glycolipid structures, and LPS Modeler for generation of lipopolysaccharide structures from various Gram-negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM-GUI development project. © 2016 Wiley Periodicals, Inc.

An extensive simulation study of lipid bilayer properties with different head groups, acyl chain lengths, and chain saturations

Previous MD simulations of six phosphocholine (PC) lipid bilayers demonstrated the accuracy of the CHARMM36 force field (C36FF) for PC bilayer simulation at varied temperatures (BBA-Biomembranes, 1838 (2014): 2520-2529). In this work, we further examine the accuracy of C36FF over a wide temperature range for a broader range of lipid types such as various head groups (phosphatidic acid (PA), PC, phosphoethanolamine (PE), phosphoglycerol (PG), and phosphoserine (PS)), and tails (saturated, mono-, mixed- and poly-unsaturated acyl chains with varied chain lengths). The structural properties (surface area per lipid (SA/lip), overall bilayer thickness, hydrophobic thickness, headgroup-to-headgroup thickness, deuterium order parameter (S_CD), and spin-lattice relaxation time (T₁)) obtained from simulations agree well with nearly all available experimental data. Our analyses indicate that PS lipids have the most inter-lipid hydrogen bonds, while PG lipids have the most intra-lipid hydrogen bonds, which play the main role in their low SA/lip in PS lipids and low thicknesses in PG lipids, respectively. PS, PE, and PA lipids have the largest contact clusters with on average 5-8 lipids per cluster, while PC and PG have clusters of 4 lipids based on a cutoff distance of 6.5Å. PS lipids have much slower lipid wobble (i.e., higher correlation time) than other head groups at a given temperature as the hydrogen bonded network significantly reduces a lipid's mobility, and the rate of lipid wobble increases dramatically as temperature increases. These in-depth analyses facilitate further understanding of lipid bilayers at the atomic level.

Challenges in structural approaches to cell modeling

Computational modeling is essential for structural characterization of biomolecular mechanisms across the broad spectrum of scales. Adequate understanding of biomolecular mechanisms inherently involves our ability to model them. Structural modeling of individual biomolecules and their interactions has been rapidly progressing. However, in terms of the broader picture, the focus is shifting toward larger systems, up to the level of a cell. Such modeling involves a more dynamic and realistic representation of the interactomes in vivo, in a crowded cellular environment, as well as membranes and membrane proteins, and other cellular components. Structural modeling of a cell complements computational approaches to cellular mechanisms based on differential equations, graph models, and other techniques to model biological networks, imaging data, etc. Structural modeling along with other computational and experimental approaches will provide a fundamental understanding of life at the molecular level and lead to important applications to biology and medicine. A cross section of diverse approaches presented in this review illustrates the developing shift from the structural modeling of individual molecules to that of cell biology. Studies in several related areas are covered: biological networks; automated construction of three-dimensional cell models using experimental data; modeling of protein complexes; prediction of non-specific and transient protein interactions; thermodynamic and kinetic effects of crowding; cellular membrane modeling; and modeling of chromosomes. The review presents an expert opinion on the current state-of-the-art in these various aspects of structural modeling in cellular biology, and the prospects of future developments in this emerging field.

Energy Decomposition Analysis Based on Absolutely Localized Molecular Orbitals for Large-Scale Density Functional Theory Calculations in Drug Design

We report the development and implementation of an energy decomposition analysis (EDA) scheme in the ONETEP linear-scaling electronic structure package. Our approach is hybrid as it combines the localized molecular orbital EDA (Su, P.; Li, H. J. Chem. Phys., 2009, 131, 014102) and the absolutely localized molecular orbital EDA (Khaliullin, R. Z.; et al. J. Phys. Chem. A, 2007, 111, 8753-8765) to partition the intermolecular interaction energy into chemically distinct components (electrostatic, exchange, correlation, Pauli repulsion, polarization, and charge transfer). Limitations shared in EDA approaches such as the issue of basis set dependence in polarization and charge transfer are discussed, and a remedy to this problem is proposed that exploits the strictly localized property of the ONETEP orbitals. Our method is validated on a range of complexes with interactions relevant to drug design. We demonstrate the capabilities for large-scale calculations with our approach on complexes of thrombin with an inhibitor comprised of up to 4975 atoms. Given the capability of ONETEP for large-scale calculations, such as on entire proteins, we expect that our EDA scheme can be applied in a large range of biomolecular problems, especially in the context of drug design.

Conformational dynamics of a hydrophobic prion fragment (113-127) in different pH and osmolyte solutions

Prion diseases are characterized by a conformational change in prion protein from its native state into beta-sheet rich aggregates that are neurotoxic. The central domain that contain a highly conserved hydrophobic region of the protein play an important role in the toxicity. The conformation of the proteins is largely influenced by various solvent environments. Here we report results of study of hydrophobic prion fragment peptide PrP(113-127) under different pH and osmolytes solution conditions. The secondary structure and the folding of PrP(113-127) was determined using circular dichroism and fluorescence spectroscopic methods. The results indicate that PrP(113-127) adopts a random coil conformation in aqueous buffer at neutral pH and that converted into beta sheet on aging. Even though the initial random coil conformation was similar in different pH conditions, the acidic as well as basic pH conditions delays the conformational transition to beta sheet. FRET results indicate that the distance between N and C-terminal regions increased on aging due to unfolding by self-assembly of the peptide into an organized beta sheet structure. Presence of osmolytes, prevented or decelerated the aggregation process of PrP(113-127) peptide.

Molecular dynamics simulations of biological membranes and membrane proteins using enhanced conformational sampling algorithms

This paper reviews various enhanced conformational sampling methods and explicit/implicit solvent/membrane models, as well as their recent applications to the exploration of the structure and dynamics of membranes and membrane proteins. Molecular dynamics simulations have become an essential tool to investigate biological problems, and their success relies on proper molecular models together with efficient conformational sampling methods. The implicit representation of solvent/membrane environments is reasonable approximation to the explicit all-atom models, considering the balance between computational cost and simulation accuracy. Implicit models can be easily combined with replica-exchange molecular dynamics methods to explore a wider conformational space of a protein. Other molecular models and enhanced conformational sampling methods are also briefly discussed. As application examples, we introduce recent simulation studies of glycophorin A, phospholamban, amyloid precursor protein, and mixed lipid bilayers and discuss the accuracy and efficiency of each simulation model and method. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.