Transmembrane fragment structures of amyloid precursor protein depend on membrane surface curvature

Exploring Free Energy Landscapes for Protein Partitioning into Membrane Domains in All-Atom and Coarse-Grained Simulations

It is known that membrane environment can impact the structure and function of integral membrane proteins. As such, elucidation of the thermodynamic driving forces governing protein partitioning between membrane domains of varying lipid composition is a fundamental topic in membrane biophysics. Molecular dynamics simulations provide valuable tools for quantitatively characterizing the free energy landscapes governing protein partitioning at the molecular level. In this study, we propose an efficient simulation methodology for the calculation of free energies for the partitioning of transmembrane proteins between liquid-disorder (Ld) and liquid-ordered (Lo) domains in all-atom (AA) phase-separated lipid bilayers. The computed potential of mean force defining the equilibrium partition coefficients is compared for AA and coarse-grained systems. Energy decomposition is used to identify differences in the underlying thermodynamics. Our findings highlight the importance of employing AA models to accurately estimate relevant free energy changes during protein translation between membrane domains.

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.

Nanodisc Reconstitution and Characterization of Amyloid-β Precursor Protein C99

Amyloid precursor protein (APP) plays a pivotal role in the pathology of Alzheimer's disease (AD). Since the fragmentation of the membrane-bound APP that results in the production of amyloid-β peptides is the starting point for amyloid toxicity in AD, it is important to investigate the structure and dynamics of APP in a near-native lipid-bilayer environment. However, the reconstitution of APP into a stable and suitable membrane-mimicking lipid environment is a challenging task. In this study, the 99-residue C-terminal domain of APP is successfully reconstituted into polymer nanodiscs and characterized using size-exclusion chromatography, mass spectrometry, solution NMR, and magic-angle spinning solid-state NMR. In addition, the feasibility of using lipid-solubilizing polymers for isolating and characterizing APP in the native Escherichia. coli membrane environment is demonstrated.

Nanodisc reconstitution and characterization of amyloid-β precursor protein C99

Amyloid precursor protein (APP) plays a pivotal role in the pathology of Alzheimer's disease. Since the fragmentation of the membrane-bound APP that results in the production of amyloid-beta peptides is the starting point for amyloid toxicity in AD, it is important to investigate the structure and dynamics of APP in a near-native lipid-bilayer environment. However, the reconstitution of APP into a stable/suitable membrane-mimicking lipid environment is a challenging task. In this study, the 99-residue C-terminal domain of APP is successfully reconstituted into polymer nanodiscs and characterized using size-exclusion chromatography, mass spectrometry, solution NMR, and magic-angle spinning solid-state NMR. In addition, the feasibility of using lipid-solubilizing polymers for isolating and characterizing APP in native E. coli membrane environment is demonstrated.

Interaction of Substrates with γ-Secretase at the Level of Individual Transmembrane Helices-A Methodological Approach

Intramembrane proteases, such as γ secretase, typically recruit multiple substrates from an excess of single-span membrane proteins. It is currently unclear to which extent substrate recognition depends on specific interactions of their transmembrane domains (TMDs) with TMDs of a protease. Here, we investigated a large number of potential pairwise interactions between TMDs of γ secretase and a diverse set of its substrates using two different configurations of BLaTM, a genetic reporter system. Our results reveal significant interactions between TMD2 of presenilin, the enzymatic subunit of γ secretase, and the TMD of the amyloid precursor protein, as well as of several other substrates. Presenilin TMD2 is a prime candidate for substrate recruitment, as has been shown from previous studies. In addition, the amyloid precursor protein TMD enters interactions with presenilin TMD 4 as well as with the TMD of nicastrin. Interestingly, the Gly-rich interfaces between the amyloid precursor protein TMD and presenilin TMDs 2 and 4 are highly similar to its homodimerization interface. In terms of methodology, the economics of the newly developed library-based method could prove to be a useful feature in related future work for identifying heterotypic TMD-TMD interactions within other biological contexts.

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.

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.

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.

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.

Altered Hinge Conformations in APP Transmembrane Helix Mutants May Affect Enzyme-Substrate Interactions of γ-Secretase

Cleavage of substrates by γ-secretase is an inherently slow process where substrate-enzyme affinities cannot be broken down into specific sequence requirements in contrast to soluble proteases. Nevertheless, despite its apparent sequence tolerance single point mutations in amyloid precursor protein can severely affect cleavage efficiencies and change product line preferences. We have determined by NMR spectroscopy the structures of the transmembrane domain of amyloid precursor protein in TFE/water and compared it to that of four mutants: two FAD mutants, V44M and I45T, and the two diglycine hinge mutants, G38L and G38P. In accordance with previous publications, the transmembrane domain is composed of two helical segments connected by the diglycine hinge. Mutations alter kink angles and structural flexibility. Furthermore, to our surprise, we observe different, but specific mutual orientations of N- and C-terminal helical segments in the four mutants compared to the wildtype. We speculate that the observed orientations for G38L, G38P, V44M, and I45T lead to unfavorable interactions with γ-secretase exosites during substrate movement to the enzyme's active site in presenilin and/or for the accommodation into the substrate-binding cavity of presenilin.

Modulating Hinge Flexibility in the APP Transmembrane Domain Alters γ-Secretase Cleavage

Intramembrane cleavage of the β-amyloid precursor protein C99 substrate by γ-secretase is implicated in Alzheimer's disease pathogenesis. Biophysical data have suggested that the N-terminal part of the C99 transmembrane domain (TMD) is separated from the C-terminal cleavage domain by a di-glycine hinge. Because the flexibility of this hinge might be critical for γ-secretase cleavage, we mutated one of the glycine residues, G38, to a helix-stabilizing leucine and to a helix-distorting proline. Both mutants impaired γ-secretase cleavage and also altered its cleavage specificity. Circular dichroism, NMR, and backbone amide hydrogen/deuterium exchange measurements as well as molecular dynamics simulations showed that the mutations distinctly altered the intrinsic structural and dynamical properties of the substrate TMD. Although helix destabilization and/or unfolding was not observed at the initial ε-cleavage sites of C99, subtle changes in hinge flexibility were identified that substantially affected helix bending and twisting motions in the entire TMD. These resulted in altered orientation of the distal cleavage domain relative to the N-terminal TMD part. Our data suggest that both enhancing and reducing local helix flexibility of the di-glycine hinge may decrease the occurrence of enzyme-substrate complex conformations required for normal catalysis and that hinge mobility can thus be conducive for productive substrate-enzyme interactions.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

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.

Molecular dynamics of C99-bound γ-secretase reveal two binding modes with distinct compactness, stability, and active-site retention: implications for Aβ production

The membrane protease γ-secretase cleaves the C99 fragment of the amyloid precursor protein, thus producing the Aβ peptides central to Alzheimer's disease. Cryo-electron microscopy has provided the topology but misses the membrane and loop parts that contribute to substrate binding. We report here an essentially complete atomic model of C99 within wild-type γ-secretase that respects all the experimental constraints and additionally describes loop, helix, and C99 substrate dynamics in a realistic all-atom membrane. Our model represents the matured auto-cleaved state required for catalysis. From two independent 500-ns molecular dynamic simulations, we identify two conformation states of C99 in equilibrium, a compact and a loose state. Our simulations provide a basis for C99 processing and Aβ formation and explain the production of longer and shorter Aβ, as the compact state retains C99 for longer and thus probably trims to shorter Aβ peptides. We expect pathogenic presenilin mutations to stabilize the loose over the compact state. The simulations detail the role of the Lys53-Lys54-Lys55 anchor for C99 binding, a loss of helicity of bound C99, and positioning of Thr48 and Leu49 leading to alternative trimming pathways on opposite sides of the C99 helix in three amino acid steps. The C99 binding topology resembles that of C83-bound γ-secretase without membrane but lacks a presenilin 1-C99 β-sheet, which could be induced by C83's stronger binding. The loose state should be selectively disfavored by γ-secretase modulators to increase C99 trimming and reduce the formation of longer Aβ, a strategy that is currently much explored but has lacked a structural basis.

Increased H-Bond Stability Relates to Altered ε-Cleavage Efficiency and Aβ Levels in the I45T Familial Alzheimer's Disease Mutant of APP

Cleavage of the amyloid precursor protein's (APP) transmembrane domain (TMD) by γ-secretase is a crucial step in the aetiology of Alzheimer's Disease (AD). Mutations in the APP TMD alter cleavage and lead to familial forms of AD (FAD). The majority of FAD mutations shift the preference of initial cleavage from ε49 to ε48, thus raising the AD-related Aβ42/Aβ40 ratio. The I45T mutation is among the few FAD mutations that do not alter ε-site preference, while it dramatically reduces the efficiency of ε-cleavage. Here, we investigate the impact of the I45T mutation on the backbone dynamics of the substrate TMD. Amide exchange experiments and molecular dynamics simulations in solvent and a lipid bilayer reveal an increased stability of amide hydrogen bonds at the ζ- and γ-cleavage sites. Stiffening of the H-bond network is caused by an additional H-bond between the T45 side chain and the TMD backbone, which alters dynamics within the cleavage domain. In particular, the increased H-bond stability inhibits an upward movement of the ε-sites in the I45T mutant. Thus, an altered presentation of ε-sites to the active site of γ-secretase as a consequence of restricted local flexibility provides a rationale for reduced ε-cleavage efficiency of the I45T mutant.

Structural Modeling of γ-Secretase Aβ n Complex Formation and Substrate Processing

The intramembrane aspartyl protease γ-secretase (GSEC) cleaves single-span transmembrane helices including the C-terminal fragment of the amyloid precursor protein (APP). This substrate is initially cleaved at the ϵ-site followed by successive processing (trimming) events mostly in steps of three amino acids. GSEC is responsible for the formation of N-terminal APP amyloid-β (Aβ) peptides of different length (e.g., Aβ₄₂) that can form aggregates involved in Alzheimer's disease pathogenesis. The molecular mechanism of GSEC-APP substrate recognition is key for understanding how different peptide products are formed and could help in designing APP-selective modulators. Based on the known structure of apo GSEC and the APP-C99 fragment we have generated putative structural models of the initial binding in three different possible modes using extensive molecular dynamics (MD) simulations. The binding mode with the substrate helix located in a cleft between the transmembrane helices 2 and 3 of the presenilin subunit was identified as a most likely binding mode. Based on this arrangement, the processing steps were investigated using restraint MD simulations to pull the scissile bond (for each processing step) into a transition like (cleavable) state. This allowed us to analyze in detail the motions and energetic contributions of participating residues. The structural model agrees qualitatively well with the influence of many mutations in GSEC and C99. It also explains the effects of inhibitors, cross-linking, as well as spectroscopic data on GSEC substrate binding and can serve as working model for the future planning of structural and biochemical studies.

Amyloid-β(29-42) Dimeric Conformations in Membranes Rich in Omega-3 and Omega-6 Polyunsaturated Fatty Acids

The omega-3 and omega-6 polyunsaturated fatty acids are two important components of cell membranes in human brains. When incorporated into phospholipids, omega-3 slows the progression of Alzheimer's disease (AD), whereas omega-6 is linked to increased risk of AD. Little is known on the amyloid-β (Aβ) conformations in membranes rich in omega-3 and omega-6 phospholipids. Herein, the structural properties of the Aβ_29-42 dimer embedded in both fatty acid membranes were comparatively studied to a 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine (POPC) bilayer using all-atom molecular dynamics (MD) simulations. Starting from α-helix, both omega-6 and omega-3 membranes promote new orientations and conformations of the dimer, in agreement with the observed dependence of Aβ production upon addition of these two fatty acids. This conformational result is corroborated by atomistic MD simulations of the dimer of the 99 amino acid C-terminal fragment of amyloid precursor protein spanning the residues 15-55. Starting from β-sheet, omega-6 membrane promotes helical and disordered structures of Aβ_29-42 dimer, whereas omega-3 membrane preserves the β-sheet structures differing however from those observed in POPC. Remarkably, the mixture of the two fatty acids and POPC depicts another conformational ensemble of the Aβ_29-42 dimer. This finding demonstrates that variation in the abundance of the molecular phospholipids, which changes with age, modulates membrane-embedded Aβ oligomerization.

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.

Dissecting conformational changes in APP's transmembrane domain linked to ε-efficiency in familial Alzheimer's disease

The mechanism by which familial Alzheimer's disease (FAD) mutations within the transmembrane domain (TMD) of the Amyloid Precursor Protein (APP) affect ε-endoproteolysis is only poorly understood. Thereby, mutations in the cleavage domain reduce ε-efficiency of γ-secretase cleavage and some even shift entry into production lines. Since cleavage occurs within the TMD, a relationship between processing and TMD structure and dynamics seems obvious. Using molecular dynamic simulations, we dissect the dynamic features of wild-type and seven FAD-mutants into local and global components. Mutations consistently enhance hydrogen-bond fluctuations upstream of the ε-cleavage sites but maintain strong helicity there. Dynamic perturbation-response scanning reveals that FAD-mutants target backbone motions utilized in the bound state. Those motions, obscured by large-scale motions in the pre-bound state, provide (i) a dynamic mechanism underlying the proposed coupling between binding and ε-cleavage, (ii) key sites consistent with experimentally determined docking sites, and (iii) the distinction between mutants and wild-type.

Radial distribution function of semiflexible oligomers with stretching flexibility

The radial distribution of the end-to-end distance R_ee is crucial for quantifying the global size and flexibility of a linear polymer. For semiflexible polymers, several analytical formulas have been derived for the radial distribution of R_ee ignoring the stretching flexibility. However, for semiflexible oligomers, such as DNA or RNA, the stretching flexibility can be rather pronounced and can significantly affect the radial distribution of R_ee. In this study, we obtained an extended formula that includes the stretch modulus to describe the distribution of R_ee for semiflexible oligomers on the basis of previous formulas for semiflexible polymers without stretching flexibility. The extended formula was validated by extensive Monte Carlo simulations over wide ranges of the stretch modulus and persistence length, as well as all-atom molecular dynamics simulations of short DNAs and RNAs. Additionally, our analyses showed that the effect of stretching flexibility on the distribution of R_ee becomes negligible for DNAs longer than ∼130 base pairs and RNAs longer than ∼240 base pairs.

Impact of membrane curvature on amyloid aggregation

The misfolding, amyloid aggregation, and fibril formation of intrinsically disordered proteins/peptides (or amyloid proteins) have been shown to cause a number of disorders. The underlying mechanisms of amyloid fibrillation and structural properties of amyloidogenic precursors, intermediates, and amyloid fibrils have been elucidated in detail; however, in-depth examinations on physiologically relevant contributing factors that induce amyloidogenesis and lead to cell death remain challenging. A large number of studies have attempted to characterize the roles of biomembranes on protein aggregation and membrane-mediated cell death by designing various membrane components, such as gangliosides, cholesterol, and other lipid compositions, and by using various membrane mimetics, including liposomes, bicelles, and different types of lipid-nanodiscs. We herein review the dynamic effects of membrane curvature on amyloid generation and the inhibition of amyloidogenic proteins and peptides, and also discuss how amyloid formation affects membrane curvature and integrity, which are key for understanding relationships with cell death. Small unilamellar vesicles with high curvature and large unilamellar vesicles with low curvature have been demonstrated to exhibit different capabilities to induce the nucleation, amyloid formation, and inhibition of amyloid-β peptides and α-synuclein. Polymorphic amyloidogenesis in small unilamellar vesicles was revealed and may be viewed as one of the generic properties of interprotein interaction-dominated amyloid formation. Several mechanical models and phase diagrams are comprehensively shown to better explain experimental findings. The negative membrane curvature-mediated mechanisms responsible for the toxicity of pancreatic β cells by the amyloid aggregation of human islet amyloid polypeptide (IAPP) and binding of the precursors of the semen-derived enhancer of viral infection (SEVI) are also described. The curvature-dependent binding modes of several types of islet amyloid polypeptides with high-resolution NMR structures are also discussed.

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.

Molecular details of spontaneous insertion and interaction of HCV non-structure 3 protease protein domain with PIP2-containing membrane

Hepatitis C virus (HCV), known as the leading cause of liver cirrhosis, viral hepatitis, and hepatocellular carcinoma, has been affecting more than 150 million people globally. The HCV non-structure 3 (NS3) protease protein domain plays a key role in HCV replication and pathogenesis; and is currently a primary target for HCV antiviral therapy. Through unbiased molecular dynamics simulations which take advantage of the novel highly mobile membrane mimetic model, we constructed the membrane-bound state of the protein domain at the atomic level. Our results indicated that protease domain of HCV NS3 protein can spontaneously bind and penetrate to an endoplasmic reticulum complex membrane containing phosphatidylinositol 4,5-bisphosphate (PIP2). An amphipathic helix α₀ and loop S1 show their anchoring role to keep the protein on the membrane surface. Proper orientation of the protein domain at membrane surface was identified through measuring tilt angles of two specific vectors, wherein residue R161 plays a crucial role in its final orientation. Remarkably, PIP2 molecules were observed to bind to three main sites of the protease domain via specific electrostatic contacts and hydrogen bonds. PIP2-interaction determines the protein orientation at the membrane while both hydrophobic interplay and PIP2-interaction can stabilize the NS3 - membrane complex. Simulated results provide us with a detailed characterization of insertion, orientation and PIP2-interaction of NS3 protease domain at membrane environment, thus enhancing our understanding of structural functions and mechanism for the association of HCV non-structure 3 protein with respect to ER membranes.

Structural Characterization of the Amyloid Precursor Protein Transmembrane Domain and Its γ-Cleavage Site

Alzheimer's disease is the most common form of dementia that affects about 50 million of sufferers worldwide. A major role for the initiation and progression of Alzheimer's disease has been associated with the amyloid β-peptide (Aβ), which is a protease cleavage product of the amyloid precursor protein. The amyloid precursor protein is an integral membrane protein with a single transmembrane domain. Here, we assessed the structural integrity of the transmembrane domain within oriented phosphatidylcholine lipid bilayers and determined the tilt angle distribution and dynamics of various subdomains using solid-state NMR and attenuated total reflectance Fourier transform infrared spectroscopies. Although the overall secondary structure of the transmembrane domain is α-helical, pronounced conformational and topological heterogeneities were observed for the γ- and, to a lesser extent, the ζ-cleavage site, with pronounced implications for the production of Aβ and related peptides, the development of the disease, and pharmaceutical innovation.

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.

Membrane-Accelerated Amyloid-β Aggregation and Formation of Cross-β Sheets

Amyloid- β aggregates play a causative role in Alzheimer's disease. These aggregates are a product of the physical environment provided by the basic neuronal membrane, composed of a lipid bilayer. The intrinsic properties of the lipid bilayer allow amyloid- β peptides to nucleate and form well-ordered cross- β sheets within the membrane. Here, we correlate the aggregation of the hydrophobic fragment of the amyloid- β protein, A β 25 - 35 , with the hydrophobicity, fluidity, and charge density of a lipid bilayer. We summarize recent biophysical studies of model membranes and relate these to the process of aggregation in physiological systems.

Initial Substrate Binding of γ-Secretase: The Role of Substrate Flexibility

γ-Secretase cleaves transmembrane domains (TMD) of amyloid precursor protein (APP), producing pathologically relevant amyloid-β proteins. Initial substrate binding represents a key step of the γ-secretase cleavage whose mechanism remains elusive. Through long time scale coarse-grained and atomic simulations, we have found that the APP TMD can bind to the catalytic subunit presenilin 1 (PS1) on an extended surface covering PS1's TMD2/6/9 and PAL motif that are all known to be essential for enzymatic activity. This initial substrate binding could lead to reduction in the vertical gap between APP's ε-cleavage sites and γ-secretase's active center, enhanced flexibility and hydration levels around the ε-sites, and the presentation of these sites to the enzyme. There are heterogeneous substrate binding poses in which the substrate is found to bind to either the N- or C-terminal parts of PS1, or both. Moreover, we also find that the stability of the binding poses can be modulated by the flexibility of substrate TMD. Especially, the APP substrate, when deprived of bending fluctuation, does not bind to TMD9 at PS1's C-terminus. Our simulations have revealed further that another substrate of γ-secretase, namely, notch receptors, though bearing a rigid TMD, can still bind to PS1 TMD9, but by a different mechanism, suggesting that the influence of substrate flexibility is context-dependent. Together, these findings shed light on the mechanism of initial substrate docking of γ-secretase and the role of substrate flexibility in this process.

Simulations of disordered proteins and systems with conformational heterogeneity

Intrinsically disordered proteins (IDPs) and protein regions can facilitate a wide variety of complex physiological processes such as binding, signaling, and formation of membraneless organelles. They can however also play pathological roles by aggregating into cytotoxic oligomers and fibrils. Characterizing the structure and function of disordered proteins is an onerous task, primarily because these proteins adopt transient structures, which are difficult to capture in experiments. Simulations have emerged as a powerful tool for interpreting and augmenting experimental measurements of IDPs. In this review we focus on computer simulations of disordered protein structures, functions, assemblies, and emerging questions that, taken together, give an overview of the field as it exists today.

Observation of helix associations for insertion of a retinal molecule and distortions of helix structures in bacteriorhodopsin

We applied a newly proposed prediction method for membrane protein structures to bacteriorhodopsin that has distorted transmembrane helices in the native structure. This method uses an implicit membrane model, which restricts sampling space during folding in a membrane region, and includes helix bending. Replica-exchange simulations were performed with seven transmembrane helices only without a retinal molecule. Obtained structures were classified into clusters of similar structures, which correspond to local-minimum free energy states. The two lowest free energy states corresponded to a native-like structure with the correct empty space for retinal and a structure with this empty space filled with a helix. Previous experiments of bacteriorhodopsin suggested that association of transmembrane helices enables them to make a room for insertion of a retinal. Our results are consistent with these results. Moreover, distortions of helices in the native-like structures were successfully reproduced. In the distortions, whereas the locations of kinks for all helices were similar to those of Protein Data Bank's data, the amount of bends was more similar for helices away from the retinal than for those close to the retinal in the native structure. This suggests a hypothesis that the amino-acid sequence specifies the location of kinks in transmembrane helices and that the amount of distortions depends on the interactions with the surrounding molecules such as neighboring helices, lipids, and retinal.

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.

Conformational Changes in the Epidermal Growth Factor Receptor: Role of the Transmembrane Domain Investigated by Coarse-Grained MetaDynamics Free Energy Calculations

The epidermal growth factor receptor (EGFR) is a dimeric membrane protein that regulates key aspects of cellular function. Activation of the EGFR is linked to changes in the conformation of the transmembrane (TM) domain, brought about by changes in interactions of the TM helices of the membrane lipid bilayer. Using an advanced computational approach that combines Coarse-Grained molecular dynamics and well-tempered MetaDynamics (CG-MetaD), we characterize the large-scale motions of the TM helices, simulating multiple association and dissociation events between the helices in membrane, thus leading to a free energy landscape of the dimerization process. The lowest energy state of the TM domain is a right-handed dimer structure in which the TM helices interact through the N-terminal small-X₃-small sequence motif. In addition to this state, which is thought to correspond to the active form of the receptor, we have identified further low-energy states that allow us to integrate with a high level of detail a range of previous experimental observations. These conformations may lead to the active state via two possible activation pathways, which involve pivoting and rotational motions of the helices, respectively. Molecular dynamics also reveals correlation between the conformational changes of the TM domains and of the intracellular juxtamembrane domains, paving the way for a comprehensive understanding of EGFR signaling at the cell 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.

High Membrane Curvature Enhances Binding, Conformational Changes, and Fibrillation of Amyloid-β on Lipid Bilayer Surfaces

Aggregation of the amyloid-β (Aβ) protein and the formation of toxic aggregates are the possible pathogenic pathways in Alzheimer's disease. Accumulating evidence suggests that lipid membranes play key roles in protein aggregation, although the intermolecular forces that drive the interactions between Aβ-(1-40) and the membranes vary in different membrane systems. Here, we observed that a high positive curvature of lipid vesicles with diameters of ∼30 nm enhanced the association of Aβ with anionic phosphatidylglycerol membranes in the liquid-crystalline phase and with zwitterionic phosphatidylcholine membranes in the gel phase. The binding modes of Aβ to these membranes differ in terms of the location of the protein on the membrane and of the protein secondary structure. The fibrillation of Aβ was accelerated in the presence of the vesicles and at high protein-to-lipid ratios. Under these conditions, the protein accumulated on the surfaces, as demonstrated by a high (10(7) M(-1)) binding constant. Our findings suggest that packing defects on membranes with high curvatures, such as the intraluminal vesicles in multivesicular bodies and the exosomes, might result in the accumulation of toxic protein aggregates.

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.

Understanding intramembrane proteolysis: from protein dynamics to reaction kinetics

Intramembrane proteolysis - cleavage of proteins within the plane of a membrane - is a widespread phenomenon that can contribute to the functional activation of substrates and is involved in several diseases. Although different families of intramembrane proteases have been discovered and characterized, we currently do not know how these enzymes discriminate between substrates and non-substrates, how site-specific cleavage is achieved, or which factors determine the rate of proteolysis. Focusing on γ-secretase and rhomboid proteases, we argue that answers to these questions may emerge from connecting experimental readouts, such as reaction kinetics and the determination of cleavage sites, to the structures and the conformational dynamics of substrates and enzymes.

Modulation of Alzheimer's Aβ protofilament-membrane interactions by lipid headgroups

The molecular pathogenesis of Alzheimer's disease (AD) is complex and sparsely understood. The relationship between AD's amyloid β (Aβ) peptides and neuronal membranes is central to Aβ's cytotoxicity and is directly modulated by the composition of the lipid headgroups. Molecular studies of the insertion of model Aβ40 protofilaments in lipid bilayers revealed strong interactions that affect the structural integrity of both the membranes and the ordered amyloid aggregates. In particular, electrostatics plays a crucial role in the interaction between Aβ protofilaments and palmytoil-oleoyl-phosphatidylethanolamine (POPE) lipids, a common component of neuronal plasma membranes. Here, we use all-atom molecular dynamics and steered molecular dynamics simulations to systematically compare the effects that POPE and palmytoil-oleoyl-phosphatidylcholine (POPC) headgroups have on the Aβ-lipid interactions. We find that Aβ protofilaments exhibit weaker electrostatic interactions with POPC headgroups and establish significantly shorter-lived contacts with the POPC bilayer. This illustrates the crucial yet complex role of electrostatic and hydrogen bonding interactions in modulating the anchoring and insertion of Aβ peptides into lipid bilayers. Our study reveals the atomistic details behind the barrier created by the lipid headgroup region in impeding solution-aggregated fibrillar oligomers to spontaneously insert into POPC bilayers, in contrast to the POPE case. While the biological reality is notoriously more complex (e.g., including other factors such as cholesterol), our results evidence a simple experimentally and computationally testable case for probing the factors that control the insertion of Aβ oligomeric aggregates in neuronal cell membranes--a process central to their neurotoxicity.

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.

Dimerization of the EphA1 receptor tyrosine kinase transmembrane domain: Insights into the mechanism of receptor activation

EphA1 is a receptor tyrosine kinase (RTK) that plays a key role in developmental processes, including guidance of the migration of axons and cells in the nervous system. EphA1, in common with other RTKs, contains an N-terminal extracellular domain, a single transmembrane (TM) α-helix, and a C-terminal intracellular kinase domain. The TM helix forms a dimer, as seen in recent NMR studies. We have modeled the EphA1 TM dimer using a multiscale approach combining coarse-grain (CG) and atomistic molecular dynamics (MD) simulations. The one-dimensional potential of mean force (PMF) for this system, based on interhelix separation, has been calculated using CG MD simulations. This provides a view of the free energy landscape for helix–helix interactions of the TM dimer in a lipid bilayer. The resulting PMF profiles suggest two states, consistent with a rotation-coupled activation mechanism. The more stable state corresponds to a right-handed helix dimer interacting via an N-terminal glycine zipper motif, consistent with a recent NMR structure (2K1K). A second metastable state corresponds to a structure in which the glycine zipper motif is not involved. Analysis of unrestrained CG MD simulations based on representative models from the PMF calculations or on the NMR structure reveals possible pathways of interconversion between these two states, involving helix rotations about their long axes. This suggests that the interaction of TM helices in EphA1 dimers may be intrinsically dynamic. This provides a potential mechanism for signaling whereby extracellular events drive a shift in the repopulation of the underlying TM helix dimer energy landscape.

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

The 99 amino acid C-terminal fragment of amyloid precursor protein (C99), consisting of a single transmembrane (TM) helix, is known to form homodimers. Homodimers can be processed by γ-secretase to produce amyloid-β (Aβ) protein, which is implicated in Alzheimer's disease (AD). While knowledge of the structure of C99 homodimers is of great importance, experimental NMR studies and simulations have produced varying structural models, including right-handed and left-handed coiled-coils. In order to investigate the structure of this critical protein complex, simulations of the C99(15-55) homodimer in POPC membrane bilayer and DPC surfactant micelle environments were performed using a multiscale approach that blends atomistic and coarse-grained models. The C99(15-55) homodimer adopts a dominant right-handed coiled-coil topology consisting of three characteristic structural states in a bilayer, only one of which is dominant in the micelle. Our structural study, which provides a self-consistent framework for understanding a number of experiments, shows that the energy landscape of the C99 homodimer supports a variety of slowly interconverting structural states. The relative importance of any given state can be modulated through environmental selection realized by altering the membrane or micelle characteristics.

Impact of bilayer lipid composition on the structure and topology of the transmembrane amyloid precursor C99 protein

C99 (also known as β-CTF) is the 99 residue transmembrane C-terminal domain (residues 672–770) of the amyloid precursor protein and is the immediate precursor of the amyloid-β (Aβ) polypeptides. To test the dependence of the C99 structure on the composition of the host model membranes, NMR studies of C99 were conducted both in anionic lyso-myristoylphosphatidylglycerol (LMPG) micelles and in a series of five zwitterionic bicelle compositions involving phosphatidylcholine and sphingomyelin in which the acyl chain lengths of these lipid components varied from 14 to 24 carbons. Some of these mixtures are reported for the first time in this work and should be of broad utility in membrane protein research. The site-specific backbone ¹⁵N and ¹H chemical shifts for C99 in LMPG and in all five bicelle mixtures were seen to be remarkably similar, indicating little dependence of the backbone structure of C99 on the composition of the host model membrane. However, the length of the transmembrane span was seen to vary in a manner that alters the positioning of the γ-secretase cleavage sites with respect to the center of the bilayer. This observation may contribute to the known dependency of the Aβ42-to-Aβ40 production ratio on both membrane thickness and the length of the C99 transmembrane domain.