Structural heterogeneity in transmembrane amyloid precursor protein homodimer is a consequence of environmental selection

Machine Learning Derived Collective Variables for the Study of Protein Homodimerization in Membrane

The accurate calculation of equilibrium constants for protein-protein association is of fundamental importance to quantitative biology and remains an outstanding challenge for computational biophysics. Traditionally, equilibrium constants have been computed from one-dimensional free energy surfaces derived from sampling along a single collective variable. Importantly, recent advances in enhanced sampling methodology have facilitated the characterization of multidimensional free energy landscapes, often exposing multiple thermodynamically important minima missed by more restrictive sampling methods. A key to the effectiveness of this multidimensional sampling approach is the identification of collective variables that effectively define the configurational space of dissociated and associated states. Here we present the application of two machine learning methods for the unbiased determination of collective variables for enhanced sampling and analysis of protein-protein association. Our results both validate prior work, based on intuition derived collective variables, and demonstrate the effectiveness of the machine learning methods for the identification of collective variables for association reactions in complex biomolecular systems.

Cholesterol and Lipid Rafts in the Biogenesis of Amyloid-β Protein and Alzheimer's Disease

Cholesterol has been conjectured to be a modulator of the amyloid cascade, the mechanism that produces the amyloid-β (Aβ) peptides implicated in the onset of Alzheimer's disease. We propose that cholesterol impacts the genesis of Aβ not through direct interaction with proteins in the bilayer, but indirectly by inducing the liquid-ordered phase and accompanying liquid-liquid phase separations, which partition proteins in the amyloid cascade to different lipid domains and ultimately to different endocytotic pathways. We explore the full process of Aβ genesis in the context of liquid-ordered phases induced by cholesterol, including protein partitioning into lipid domains, mechanisms of endocytosis experienced by lipid domains and secretases, and pH-controlled activation of amyloid precursor protein secretases in specific endocytotic environments. Outstanding questions on the essential role of cholesterol in the amyloid cascade are identified for future studies.

Cholesterol-dependent amyloid β production: space for multifarious interactions between amyloid precursor protein, secretases, and cholesterol

Amyloid β is considered a key player in the development and progression of Alzheimer's disease (AD). Many studies investigating the effect of statins on lowering cholesterol suggest that there may be a link between cholesterol levels and AD pathology. Since cholesterol is one of the most abundant lipid molecules, especially in brain tissue, it affects most membrane-related processes, including the formation of the most dangerous form of amyloid β, Aβ42. The entire Aβ production system, which includes the amyloid precursor protein (APP), β-secretase, and the complex of γ-secretase, is highly dependent on membrane cholesterol content. Moreover, cholesterol can affect amyloidogenesis in many ways. Cholesterol influences the stability and activity of secretases, but also dictates their partitioning into specific cellular compartments and cholesterol-enriched lipid rafts, where the amyloidogenic machinery is predominantly localized. The most complicated relationships have been found in the interaction between cholesterol and APP, where cholesterol affects not only APP localization but also the precise character of APP dimerization and APP processing by γ-secretase, which is important for the production of Aβ of different lengths. In this review, we describe the intricate web of interdependence between cellular cholesterol levels, cholesterol membrane distribution, and cholesterol-dependent production of Aβ, the major player in AD.

Formation of extramembrane β-strands controls dimerization of transmembrane helices in amyloid precursor protein C99

The 99-residue C-terminal domain of amyloid precursor protein (APP-C99), precursor to amyloid beta (Aβ), is a transmembrane (TM) protein containing intrinsically disordered N- and C-terminal extramembrane domains. Using molecular dynamics (MD) simulations, we show that the structural ensemble of the C99 monomer is best described in terms of thousands of states. The C99 monomer has a propensity to form β-strand in the C-terminal extramembrane domain, which explains the slow spin relaxation times observed in paramagnetic probe NMR experiments. Surprisingly, homodimerization of C99 not only narrows the conformational ensemble from thousands to a few states through the formation of metastable β-strands in extramembrane domains but also stabilizes extramembrane α-helices. The extramembrane domain structure is observed to dramatically impact the homodimerization motif, resulting in the modification of TM domain conformations. Our study provides an atomic-level structural basis for communication between the extramembrane domains of the C99 protein and TM homodimer formation. This finding could serve as a general model for understanding the influence of disordered extramembrane domains on TM protein structure.

Impact of A2T and D23N mutations on C99 homodimer conformations

The proteolytic cleavage of C99 by γ-secretase is the last step in the production of amyloid-β (Aβ) peptides. Previous studies have shown that membrane lipid composition, cholesterol concentration, and mutation in the transmembrane helix modified the structures and fluctuations of C99. In this study, we performed atomistic molecular dynamics simulations of the homodimer of the 55-residue congener of the C-terminal domain of the amyloid protein precursor, C99(1-55), in a POPC-cholesterol lipid bilayer, and we compared the conformational ensemble of WT sequence to those of the A2T and D23N variants. These mutations are particularly interesting as the protective Alzheimer's disease (AD) A2T mutation is known to decrease Aβ production, whereas the early onset AD D23N mutation does not affect Aβ production. We found noticeable differences in the structural ensembles of the three sequences. In particular, A2T varies from both WT and D23N by having long-range effects on the population of the extracellular justamembrane helix, the interface between the G29xxx-G33xxx-G37 motifs and the fluctuations of the transmembrane helical topologies.

Effect of Cholesterol on C99 Dimerization: Revealed by Molecular Dynamics Simulations

C99 is the immediate precursor for amyloid beta (Aβ) and therefore is a central intermediate in the pathway that is believed to result in Alzheimer’s disease (AD). It has been suggested that cholesterol is associated with C99, but the dynamic details of how cholesterol affects C99 assembly and the Aβ formation remain unclear. To investigate this question, we employed coarse-grained and all-atom molecular dynamics simulations to study the effect of cholesterol and membrane composition on C99 dimerization. We found that although the existence of cholesterol delays C99 dimerization, there is no direct competition between C99 dimerization and cholesterol association. In contrast, the existence of cholesterol makes the C99 dimer more stable, which presents a cholesterol binding C99 dimer model. Cholesterol and membrane composition change the dimerization rate and conformation distribution of C99, which will subsequently influence the production of Aβ. Our results provide insights into the potential influence of the physiological environment on the C99 dimerization, which will help us understand Aβ formation and AD’s etiology.

On Computing Equilibrium Binding Constants for Protein-Protein Association in Membranes

Protein association in lipid membranes is fundamental to membrane protein function and of great biomedical relevance. All-atom and coarse-grained models have been extensively used to understand the protein-protein interactions in the membrane and to compute equilibrium association constants. However, slow translational and rotational diffusion of protein in membrane presents challenges to the effective sampling of conformations defining the ensembles of free and bound states contributing to the association equilibrium and the free energy of dimerization. We revisit the homodimerization equilibrium of the TM region of glycophorin A. Conformational sampling is performed using umbrella sampling along previously proposed one-dimensional collective variables and compared with sampling over a two-dimensional collective variable space using the MARTINI v2.2 force field. We demonstrate that the one-dimensional collective variables suffer from restricted sampling of the native homodimer conformations leading to a biased free energy landscape. Conversely, simulations along the two-dimensional collective variable effectively characterize the thermodynamically relevant native and non-native interactions contributing to the association equilibrium. These results demonstrate the challenges associated with accurately characterizing binding equilibria when multiple poses contribute to the bound state ensemble.

Effect of cholesterol on the dimerization of C99-A molecular modeling perspective

C99, a naturally occurring peptide, is a precursor of the amyloid β-peptide (Aβ) and plays an important role in the so-called amyloidogenic pathway of degradation of amyloid precursor protein. While the effect of C99's dimerization is not clearly determined, it has been hypothesized that the dimerization protects C99 from being cleaved further. Cholesterol (CHOL) is known to interact with C99 and its presence in high concentrations has been linked to an increase in the production of Aβ; however, to what extent this is correlated, and how, has not yet been determined. In this study, we systematically examine the effect of increasing cholesterol concentration on the homodimerization propensity of C99, combining unbiased atomistic molecular dynamics simulations with biased simulations using a coarse grained resolution. Through the use of umbrella sampling, we show how the presence of high levels of CHOL destabilizes the interaction between two C99 monomers. The interaction pattern between the two C99s has shifted several residues, from the N-terminal end of the transmembrane region toward the corresponding C-terminal in the presence of CHOL. The umbrella sampling shows that the presence of high levels of CHOL led to a decrease of the disassociation energy by approximately 3 kJ/mol. In conclusion, this suggests that increasing CHOL destabilizes the interaction between the two C99 monomers, which may possibly cause an increase in the production of Aβ₄₂.

Synthesis, Characterization, and Simulation of Four-Armed Megamolecules

This paper describes the synthesis, characterization, and modeling of a series of molecules having four protein domains attached to a central core. The molecules were assembled with the "megamolecule" strategy, wherein enzymes react with their covalent inhibitors that are substituted on a linker. Three linkers were synthesized, where each had four oligo(ethylene glycol)-based arms terminated in a para-nitrophenyl phosphonate group that is a covalent inhibitor for cutinase. This enzyme is a serine hydrolase and reacts efficiently with the phosphonate to give a new ester linkage at the Ser-120 residue in the active site of the enzyme. Negative-stain transmission electron microscopy (TEM) images confirmed the architecture of the four-armed megamolecules. These cutinase tetramers were also characterized by X-ray crystallography, which confirmed the active-site serine-phosphonate linkage by electron-density maps. Molecular dynamics simulations of the tetracutinase megamolecules using three different force field setups were performed and compared with the TEM observations. Using the Amberff99SB-disp + pH7 force field, the two-dimensional projection distances of the megamolecules were found to agree with the measured dimensions from TEM. The study described here, which combines high-resolution characterization with molecular dynamics simulations, will lead to a comprehensive understanding of the molecular structures and dynamics for this new class of molecules.

The C99 domain of the amyloid precursor protein resides in the disordered membrane phase

Processing of the amyloid precursor protein (APP) via the amyloidogenic pathway is associated with the etiology of Alzheimer's disease. The cleavage of APP by β-secretase to generate the transmembrane 99-residue C-terminal fragment (C99) and subsequent processing of C99 by γ-secretase to yield amyloid-β (Aβ) peptides are essential steps in this pathway. Biochemical evidence suggests that amyloidogenic processing of C99 occurs in cholesterol- and sphingolipid-enriched liquid-ordered phase membrane rafts. However, direct evidence that C99 preferentially associates with these rafts has remained elusive. Here, we tested this by quantifying the affinity of C99-GFP for raft domains in cell-derived giant plasma membrane vesicles (GPMVs). We found that C99 was essentially excluded from ordered domains in vesicles from HeLa cells, undifferentiated SH-SY5Y cells, or SH-SY5Y-derived neurons; instead, ∼90% of C99 partitioned into disordered domains. The strong association of C99 with disordered domains occurred independently of its cholesterol-binding activity or homodimerization, or of the presence of the familial Alzheimer disease Arctic mutation (APP E693G). Finally, through biochemical studies we confirmed previous results, which showed that C99 is processed in the plasma membrane by α-secretase, in addition to the well-known γ-secretase. These findings suggest that C99 itself lacks an intrinsic affinity for raft domains, implying that either i) amyloidogenic processing of the protein occurs in disordered regions of the membrane, ii) processing involves a marginal subpopulation of C99 found in rafts, or iii) as-yet-unidentified protein-protein interactions with C99 in living cells drive this protein into membrane rafts to promote its cleavage therein.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Bicelles Rich in both Sphingolipids and Cholesterol and Their Use in Studies of Membrane Proteins

How the distinctive lipid composition of mammalian plasma membranes impacts membrane protein structure is largely unexplored, partly because of the dearth of isotropic model membrane systems that contain abundant sphingolipids and cholesterol. This gap is addressed by showing that sphingomyelin and cholesterol-rich (SCOR) lipid mixtures with phosphatidylcholine can be cosolubilized by n-dodecyl-β-melibioside to form bicelles. Small-angle X-ray and neutron scattering, as well as cryo-electron microscopy, demonstrate that these assemblies are stable over a wide range of conditions and exhibit the bilayered-disc morphology of ideal bicelles even at low lipid-to-detergent mole ratios. SCOR bicelles are shown to be compatible with a wide array of experimental techniques, as applied to the transmembrane human amyloid precursor C99 protein in this medium. These studies reveal an equilibrium between low-order oligomer structures that differ significantly from previous experimental structures of C99, providing an example of how ordered membranes alter membrane protein structure.

Regulation of the alternative β-secretase meprin β by ADAM-mediated shedding

Alzheimer's Disease (AD) is the sixth-leading cause of death in industrialized countries. Neurotoxic amyloid-β (Aβ) plaques are one of the pathological hallmarks in AD patient brains. Aβ accumulates in the brain upon sequential, proteolytic processing of the amyloid precursor protein (APP) by β- and γ-secretases. However, so far disease-modifying drugs targeting β- and γ-secretase pathways seeking a decrease in the production of toxic Aβ peptides have failed in clinics. It has been demonstrated that the metalloproteinase meprin β acts as an alternative β-secretase, capable of generating truncated Aβ2-x peptides that have been described to be increased in AD patients. This indicates an important β-site cleaving enzyme 1 (BACE-1)-independent contribution of the metalloprotease meprin β within the amyloidogenic pathway and may lead to novel drug targeting avenues. However, meprin β itself is embedded in a complex regulatory network. Remarkably, the anti-amyloidogenic α-secretase a disintegrin and metalloproteinase domain-containing protein 10 (ADAM10) is a direct competitor for APP at the cell surface, but also a sheddase of inactive pro-meprin β. Overall, we highlight the current cellular, molecular and structural understanding of meprin β as alternative β-secretase within the complex protease web, regulating APP processing in health and disease.

Structural Investigation of Human Prolactin Receptor Transmembrane Domain Homodimerization in a Membrane Environment through Multiscale Simulations

It is well established that prolactin (PRL) and its receptor (PRLR) are associated with hundreds of biological functions. They have been postulated to be linked to breast and prostate cancers, and PRLR signaling has attracted considerable medical and pharmaceutical interest in the development of compounds targeting PRLR. Dimerization of the receptor through its transmembrane (TM) domain is a key step for understanding its signaling and related issues. Our multiscale simulation results revealed that its TM domain can form dimers in a membrane environment with distinct states stabilized by different residue motifs. On the basis of the simulated data, an activation mechanism of PRL with the importance of two symmetrical tryptophan residues was proposed in detail to determine the conformational change of its receptor, which is essential for signal transduction. The better knowledge of PRLR structure and its protein-protein interaction can considerably contribute to a further understanding of PRLR signaling action and thereby help to develop some new PRLR signaling-based strategies for PRL-related diseases.

Exploring the impact of proteins on the line tension of a phase-separating ternary lipid mixture

The separation of lipid mixtures into thermodynamically stable phase-separated domains is dependent on lipid composition, temperature, and system size. Using molecular dynamics simulations, the line tension between thermodynamically stable lipid domains formed from ternary mixtures of di-C16:0 PC:di-C18:2 PC:cholesterol at 40:40:20 mol. % ratio was investigated via two theoretical approaches. The line tension was found to be 3.1 ± 0.2 pN by capillary wave theory and 4.7 ± 3.7 pN by pressure tensor anisotropy approaches for coarse-grained models based on the Martini force field. Using an all-atom model of the lipid membrane based on the CHARMM36 force field, the line tension was found to be 3.6 ± 0.9 pN using capillary wave theory and 1.8 ± 2.2 pN using pressure anisotropy approaches. The discrepancy between estimates of the line tension based on capillary wave theory and pressure tensor anisotropy methods is discussed. Inclusion of protein in Martini membrane lipid mixtures was found to reduce the line tension by 25%-35% as calculated by the capillary wave theory approach. To further understand and predict the behavior of proteins in phase-separated membranes, we have formulated an analytical Flory-Huggins model and parameterized it against the simulation results. Taken together these results suggest a general role for proteins in reducing the thermodynamic cost associated with domain formation in lipid mixtures and quantifies the thermodynamic driving force promoting the association of proteins to domain interfaces.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

Helix-Switch Enables C99 Dimer Transition between the Multiple Conformations

C99 is the immediate precursor of amyloid-β (Aβ) and therefore is a central intermediate in the pathway that is believed to result in Alzheimer's disease (AD). Recent studies have shown that C99 dimerization changes the Aβ ratio, but the mechanism remains unclear. Previous studies of the C99 dimer have produced controversial structure models. To address these questions, we investigated C99 dimerization using molecular dynamics (MD) simulations. A helix-switch model was revealed in the formation and transition of the C99 dimer, and six types of conformations were identified. The different conformations show differential exposures of γ-cleavage sites and insertion depths in the bilayer, which may modulate γ-cleavage of C99 and lead to different Aβ levels. Our results redefine C99 dimerization, provide a framework to mediate the current controversial results, and give insights into the understanding of the relationship between C99 dimerization and Aβ formation.

Regimes of Complex Lipid Bilayer Phases Induced by Cholesterol Concentration in MD Simulation

Cholesterol is essential to the formation of phase-separated lipid domains in membranes. Lipid domains can exist in different thermodynamic phases depending on the molecular composition and play significant roles in determining structure and function of membrane proteins. We investigate the role of cholesterol in the structure and dynamics of ternary lipid mixtures displaying phase separation using molecular dynamics simulations, employing a physiologically relevant span of cholesterol concentration. We find that cholesterol can induce formation of three regimes of phase behavior: 1) miscible liquid-disordered bulk, 2) phase-separated, domain-registered coexistence of liquid-disordered and liquid-ordered domains, and 3) phase-separated, domain-antiregistered coexistence of liquid-disordered and newly identified nanoscopic gel domains composed of cholesterol threads we name "cholesterolic gel" domains. These findings are validated and discussed in the context of current experimental knowledge, models of cholesterol spatial distributions, and models of ternary lipid-mixture phase separation.

Bexarotene Binds to the Amyloid Precursor Protein Transmembrane Domain, Alters Its α-Helical Conformation, and Inhibits γ-Secretase Nonselectively in Liposomes

Bexarotene is a pleiotropic molecule that has been proposed as an amyloid-β (Aβ)-lowering drug for the treatment of Alzheimer's disease (AD). It acts by upregulation of an apolipoprotein E (apoE)-mediated Aβ clearance mechanism. However, whether bexarotene induces removal of Aβ plaques in mouse models of AD has been controversial. Here, we show by NMR and CD spectroscopy that bexarotene directly interacts with and stabilizes the transmembrane domain α-helix of the amyloid precursor protein (APP) in a region where cholesterol binds. This effect is not mediated by changes in membrane lipid packing, as bexarotene does not share with cholesterol the property of inducing phospholipid condensation. Bexarotene inhibited the intramembrane cleavage by γ-secretase of the APP C-terminal fragment C99 to release Aβ in cell-free assays of the reconstituted enzyme in liposomes, but not in cells, and only at very high micromolar concentrations. Surprisingly, in vitro, bexarotene also inhibited the cleavage of Notch1, another major γ-secretase substrate, demonstrating that its inhibition of γ-secretase is not substrate specific and not mediated by acting via the cholesterol binding site of C99. Our data suggest that bexarotene is a pleiotropic molecule that interfere with Aβ metabolism through multiple mechanisms.

Structure of APP-C991-99 and implications for role of extra-membrane domains in function and oligomerization

The 99 amino acid C-terminal fragment of Amyloid Precursor Protein APP-C99 (C99) is cleaved by γ-secretase to form Aβ peptide, which plays a critical role in the etiology of Alzheimer's Disease (AD). The structure of C99 consists of a single transmembrane domain flanked by intra and intercellular domains. While the structure of the transmembrane domain has been well characterized, little is known about the structure of the flanking domains and their role in C99 processing by γ-secretase. To gain insight into the structure of full-length C99, REMD simulations were performed for monomeric C99 in model membranes of varying thickness. We find equilibrium ensembles of C99 from simulation agree with experimentally-inferred residue insertion depths and protein backbone chemical shifts. In thin membranes, the transmembrane domain structure is correlated with extra-membrane structural states and the extra-membrane domain structural states become less correlated to each other. Mean and variance of the transmembrane and G₃₇G₃₈ hinge angles are found to increase with thinning membrane. The N-terminus of C99 forms β-strands that may seed aggregation of Aβ on the membrane surface, promoting amyloid formation. In thicker membranes the N-terminus forms α-helices that interact with the nicastrin domain of γ-secretase. The C-terminus of C99 becomes more α-helical as the membrane thickens, forming structures that may be suitable for binding by cytoplasmic proteins, while C-terminal residues essential to cytotoxic function become α-helical as the membrane thins. The heterogeneous but discrete extra-membrane domain states analyzed here open the path to new investigations of the role of C99 structure and membrane in amyloidogenesis. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy.

Substrate processing in intramembrane proteolysis by γ-secretase - the role of protein dynamics

Intramembrane proteases comprise a number of different membrane proteins with different types of catalytic sites. Their common denominator is cleavage within the plane of the membrane, which usually results in peptide bond scission within the transmembrane helices of their substrates. Despite recent progress in the determination of high-resolution structures, as illustrated here for the γ-secretase complex and its substrate C99, it is still unknown how these enzymes function and how they distinguish between substrates and non-substrates. In principle, substrate/non-substrate discrimination could occur at the level of substrate binding and/or cleavage. Focusing on the γ-secretase/C99 pair, we will discuss recent observations suggesting that global motions within a substrate transmembrane helix may be much more important for defining a substrate than local unraveling at cleavage sites.

Dimerization of the transmembrane domain of amyloid precursor protein is determined by residues around the γ-secretase cleavage sites

One of the hallmarks of Alzheimer's disease is the formation of extracellular amyloid plaques that consist mainly of abnormally aggregated forms of amyloid β (Aβ) peptides. These peptides are generated by γ-secretase-catalyzed cleavage of a dimeric membrane-bound C-terminal fragment (C99) of the amyloid precursor protein. Although C99 homodimerization has been linked to Aβ production and changes in the aggregation-determining Aβ42/Aβ40 ratio, the motif through which C99 dimerizes has remained controversial. Here, we have used two independent assays to gain insight into C99 homodimerization in the context of the membrane of live cells: bioluminescence resonance energy transfer and Tango membrane protein-protein interaction assays, which were further confirmed by traditional pull-down assays. Our results indicate a four-amino acid region within the C99 transmembrane helix (⁴³TVIV⁴⁶) as well as its local secondary structure as critical determinants for homodimerization. These four amino acids are also a hot spot of familial Alzheimer's disease-linked mutations that both decrease C99 homodimerization and γ-secretase cleavage and alter the initial cleavage site to increase the Aβ42/40 ratio.

Exploring the structure and stability of cholesterol dimer formation in multicomponent lipid bilayers

For 40 years, the existence and possible functional importance of cholesterol dimer formation has been discussed. Due to challenges associated with structural studies of membrane lipids, there has as yet been no direct experimental verification of the existence and relevance of the cholesterol dimer. Building on recent advances in lipid force fields for molecular simulation, in this work the structure and stability of the cholesterol dimer is characterized in POPC bilayers in absence and presence of sphingomyelin. The cholesterol dimer structural ensemble is found to consist of sub-states that reflect, but also differ from, previously proposed dimer structures. While face-to-face dimer structures predominate, no evidence is found for the existence of tail-to-tail dimers in POPC lipid bilayers. Near stoichiometric complex formation of cholesterol with sphingomyelin is found to effect cholesterol dimer structure without impacting population. Comparison with NMR-derived order parameters provide validation for the simulation model employed and conclusions drawn related to the structure and stability of cholesterol dimers in multicomponent lipid bilayers. © 2016 Wiley Periodicals, Inc.

Structural and biochemical differences between the Notch and the amyloid precursor protein transmembrane domains

γ-Secretase cleavage of the Notch receptor transmembrane domain is a critical signaling event for various cellular processes. Efforts to develop inhibitors of γ-secretase cleavage of the amyloid-β precursor C99 protein as potential Alzheimer's disease therapeutics have been confounded by toxicity resulting from the inhibition of normal cleavage of Notch. We present biochemical and structural data for the combined transmembrane and juxtamembrane Notch domains (Notch-TMD) that illuminate Notch signaling and that can be compared and contrasted with the corresponding traits of C99. The Notch-TMD and C99 have very different conformations, adapt differently to changes in model membrane hydrophobic span, and exhibit different cholesterol-binding properties. These differences may be exploited in the design of agents that inhibit cleavage of C99 while allowing Notch cleavage.

Specific Binding of Cholesterol to C99 Domain of Amyloid Precursor Protein Depends Critically on Charge State of Protein

Recent NMR chemical shift measurements of the 99 residue C-terminal fragment of amyloid precursor protein (APP-C99) in the presence of cholesterol provide evidence of binary complex formation between C99 and cholesterol in membrane mimetic environments. It has also been observed that the production of Aβ protein is enhanced under conditions of high cholesterol concentration. In this study, we investigated the impact of the charge state of C99 on the structure and stability of the C99-cholesterol complex. We observed that the binding of C99 to cholesterol depends critically on the charge state of Glu 693 (E22) and Asp 694 (D23). Evaluation of the pKa values of the Asp and Glu side chains suggests that these residues may be predominantly neutral in existing experimental observations of a stable C99-cholesterol complex at lower pH (characteristic of the endosomal environment), while binding is destabilized near neutral pH (characteristic of the cytoplasm). These observations suggest that specific binding of cholesterol to C99 is a sensitive function of the pH encountered in vivo, with key E22 and D23 residues serving as a "pH switch" controlling C99-cholesterol binding.

Effect of the Synaptic Plasma Membrane on the Stability of the Amyloid Precursor Protein Homodimer

The proteolytic cleavage of the transmembrane (TM) domain of the amyloid precursor protein (APP) releases amyloid-β (Aβ) peptides, which accumulation in the brain tissue is an early indicator of Alzheimer's disease. We used multiscale molecular dynamics simulations to investigate the stability of APP-TM dimer in realistic models of the synaptic plasma membrane (SPM). Between the two possible dimerization motifs proposed by NMR and EPR, namely G709XXXA713 and G700XXXG704XXXG708, our study revealed that the dimer promoted by the G709XXXA713 motif is not stable in the SPM due to the competition with highly unsaturated lipids that constitute the SPM. Under the same conditions, the dimer promoted by the G700XXXG704XXXG708 motif is instead the most stable species and likely the most biologically relevant. Independently of the dimerization state, both these motifs can be involved in the recruitment of cholesterol molecules.

Impact of membrane lipid composition on the structure and stability of the transmembrane domain of amyloid precursor protein

Cleavage of the amyloid precursor protein (APP) by γ-secretase is a crucial first step in the evolution of Alzheimer's disease. To discover the cleavage mechanism, it is urgent to predict the structures of APP monomers and dimers in varying membrane environments. We determined the structures of the C9923-55 monomer and homodimer as a function of membrane lipid composition using a multiscale simulation approach that blends atomistic and coarse-grained models. We demonstrate that the C9923-55 homodimer structures form a heterogeneous ensemble with multiple conformational states, each stabilized by characteristic interpeptide interactions. The relative probabilities of each conformational state are sensitive to the membrane environment, leading to substantial variation in homodimer peptide structure as a function of membrane lipid composition or the presence of an anionic lipid environment. In contrast, the helicity of the transmembrane domain of monomeric C991-55 is relatively insensitive to the membrane lipid composition, in agreement with experimental observations. The dimer structures of human EphA2 receptor depend on the lipid environment, which we show is linked to the location of the structural motifs in the dimer interface, thereby establishing that both sequence and membrane composition modulate the complete energy landscape of membrane-bound proteins. As a by-product of our work, we explain the discrepancy in structures predicted for C99 congener homodimers in membrane and micelle environments. Our study provides insight into the observed dependence of C99 protein cleavage by γ-secretase, critical to the formation of amyloid-β protein, on membrane thickness and lipid composition.

Growth Factor Identity Is Encoded by Discrete Coiled-Coil Rotamers in the EGFR Juxtamembrane Region

Binding of transforming growth factor α (TGF-α) to the epidermal growth factor receptor (EGFR) extracellular domain is encoded through the formation of a unique antiparallel coiled coil within the juxtamembrane segment. This new coiled coil is an "inside-out" version of the coiled coil formed in the presence of epidermal growth factor (EGF). A third, intermediary coiled-coil interface is formed in the juxtamembrane region when EGFR is stimulated with betacellulin. The seven growth factors that activate EGFR in mammalian systems (EGF, TGF-α, epigen, epiregulin, betacellulin, heparin-binding EGF, and amphiregulin) fall into distinct categories in which the structure of the coiled coil induced within the juxtamembrane region correlates with cell state. The observation that coiled-coil state tracks with the downstream signaling profiles for each ligand provides evidence for growth factor functional selectivity by EGFR. Encoding growth factor identity in alternative coiled-coil rotamers provides a simple and elegant method for communicating chemical information across the plasma membrane.

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.

Extension of a protein docking algorithm to membranes and applications to amyloid precursor protein dimerization

Novel adjustments are introduced to the docking algorithm, DOCK/PIERR, for the purpose of predicting structures of transmembrane protein complexes. Incorporating knowledge about the membrane environment is shown to significantly improve docking accuracy. The extended version of DOCK/PIERR is shown to perform comparably to other leading docking packages. This membrane version of DOCK/PIERR is applied to the prediction of coiled-coil homodimer structures of the transmembrane region of the C-terminal peptide of amyloid precursor protein (C99). Results from MD simulation of the C99 homodimer in POPC bilayer and docking are compared. Docking results are found to capture key aspects of the homodimer ensemble, including the existence of three topologically distinct conformers. Furthermore, the extended version of DOCK/PIERR is successful in capturing the effects of solvation in membrane and micelle. Specifically, DOCK/PIERR reproduces essential differences in the homodimer ensembles simulated in POPC bilayer and DPC micelle, where configurational entropy and surface curvature effects bias the handedness and topology of the homodimer ensemble.

Specific Binding of Cholesterol to the Amyloid Precursor Protein: Structure of the Complex and Driving Forces Characterized in Molecular Detail

C99 is the C-terminal membrane-bound fragment of the amyloid precursor protein that is cleaved by γ-secretase to release Aβ peptides, the hallmark of Alzheimer's disease (AD). Specific interactions of C99 with cholesterol have been proposed to underlie the recognized role of cholesterol in promoting amyloidogenesis. By using molecular dynamics simulations, we studied cholesterol binding to C99 in a lipid bilayer. We determined the free-energy profile of binding and analyzed the structure of C99/cholesterol complexes in two low-energy binding modes. We also examined the complexation driving forces and found, unexpectedly, that the interactions between the GxxxG dimerization motif and the cholesterol ring system are not sufficient for binding and that further stabilization mediated by the C99 N-terminal domain is essential. Taken together, our results strongly support the view that C99 specifically binds cholesterol in the cell membrane; the detailed information on the structure and energetics of the complex may assist in the design of new anti-AD drugs.