Conformational fluctuations and phases in fused in sarcoma (FUS) low-complexity domain

The well-known phenomenon of phase separation in synthetic polymers and proteins has become a major topic in biophysics because it has been invoked as a mechanism of compartment formation in cells, without the need for membranes. Most of the coacervates (or condensates) are composed of Intrinsically Disordered Proteins (IDPs) or regions that are structureless, often in interaction with RNA and DNA. One of the more intriguing IDPs is the 526-residue RNA-binding protein, Fused in Sarcoma (FUS), whose monomer conformations and condensates exhibit unusual behavior that are sensitive to solution conditions. By focussing principally on the N-terminus low-complexity domain (FUS-LC comprising residues 1-214) and other truncations, we rationalize the findings of solid-state NMR experiments, which show that FUS-LC adopts a non-polymorphic fibril structure (core-1) involving residues 39-95, flanked by fuzzy coats on both the N- and C-terminal ends. An alternate structure (core-2), whose free energy is comparable to core-1, emerges only in the truncated construct (residues 110-214). Both core-1 and core-2 fibrils are stabilized by a Tyrosine ladder as well as hydrophilic interactions. The morphologies (gels, fibrils, and glass-like) adopted by FUS seem to vary greatly, depending on the experimental conditions. The effect of phosphorylation is site-specific. Simulations show that phosphorylation of residues within the fibril has a greater destabilization effect than residues that are outside the fibril region, which accords well with experiments. Many of the peculiarities associated with FUS may also be shared by other IDPs, such as TDP43 and hnRNPA2. We outline a number of problems for which there is no clear molecular explanation.

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. Expected final online publication date for the Annual Review of Biophysics, Volume 53 is May 2024. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The role of structural heterogeneity in the homodimerization of transmembrane proteins

The equilibrium association of transmembrane proteins plays a fundamental role in membrane protein function and cellular signaling. While the study of the equilibrium binding of single pass transmembrane proteins has received significant attention in experiment and simulation, the accurate assessment of equilibrium association constants remains a challenge to experiment and simulation. In experiment, there remain wide variations in association constants derived from experimental studies of the most widely studied transmembrane proteins. In simulation, state-of-the art methods have failed to adequately sample the thermodynamically relevant structures of the dimer state ensembles using coarse-grained models. In addition, all-atom force fields often fail to accurately assess the relative free energies of the dimer and monomer states. Given the importance of this fundamental biophysical process, it is essential to address these shortcomings. In this work, we establish an effective computational protocol for the calculation of equilibrium association constants for transmembrane homodimer formation. A set of transmembrane protein homodimers, used in the parameterization of the MARTINI v3 force field, are simulated using metadynamics, based on three collective variables. The method is found to be accurate and computationally efficient, providing a standard to be used in the future simulation studies using coarse-grained or all-atom models.

Characterizing the transmembrane domains of ADAM10 and BACE1 and the impact of membrane composition

The β-secretase, BACE1, and the α-secretase, ADAM10, are known to competitively cleave amyloid precursor protein (APP) in the amyloid cascades of Alzheimer's disease. Cleavage of APP by BACE1 produces a 99-residue C-terminal peptide (APP-C99) that is subsequently cleaved by γ-secretase to form amyloid-β (Aβ) protein, whereas cleavage of APP by ADAM10 is nonamyloidogenic. It has been speculated that ADAM10/APP and BACE1/APP interactions are regulated by colocalization within and outside of liquid-ordered membrane domains; however, the mechanism of this regulation and the character of the proteins' transmembrane domains are not well understood. In this work, we have developed and characterized minimal congener sequences for the transmembrane domains of ADAM10 and BACE1 using a multiscale modeling approach combining both temperature replica exchange and conventional molecular dynamics simulations based on the coarse-grained Martini2.2 and all-atom CHARMM36 force fields. Our results show that membrane composition impacts the character of the transmembrane domains of BACE1 and ADAM10, adding credence to the speculation that membrane domains are involved in the etiology of Alzheimer's disease.

Synergistic computational and experimental studies of a phosphoglycosyl transferase membrane/ligand ensemble

Complex glycans serve essential functions in all living systems. Many of these intricate and byzantine biomolecules are assembled employing biosynthetic pathways wherein the constituent enzymes are membrane-associated. A signature feature of the stepwise assembly processes is the essentiality of unusual linear long-chain polyprenol phosphate-linked substrates of specific isoprene unit geometry, such as undecaprenol phosphate (UndP) in bacteria. How these enzymes and substrates interact within a lipid bilayer needs further investigation. Here, we focus on a small enzyme, PglC from Campylobacter, structurally characterized for the first time in 2018 as a detergent-solubilized construct. PglC is a monotopic phosphoglycosyl transferase that embodies the functional core structure of the entire enzyme superfamily and catalyzes the first membrane-committed step in a glycoprotein assembly pathway. The size of the enzyme is significant as it enables high-level computation and relatively facile, for a membrane protein, experimental analysis. Our ensemble computational and experimental results provided a high-level view of the membrane-embedded PglC/UndP complex. The findings suggested that it is advantageous for the polyprenol phosphate to adopt a conformation in the same leaflet where the monotopic membrane protein resides as opposed to additionally disrupting the opposing leaflet of the bilayer. Further, the analysis showed that electrostatic steering acts as a major driving force contributing to the recognition and binding of both UndP and the soluble nucleotide sugar substrate. Iterative computational and experimental mutagenesis support a specific interaction of UndP with phosphoglycosyl transferase cationic residues and suggest a role for critical conformational transitions in substrate binding and specificity.

Conformational Fluctuations and Phases in Fused in Sarcoma (FUS) Low-Complexity Domain

The well known phenomenon of phase separation in synthetic polymers and proteins has become a major topic in biophysics because it has been invoked as a mechanism of compartment formation in cells, without the need for membranes. Most of the coacervates (or condensates) are composed of Intrinsically Disordered Proteins (IDPs) or regions that are structureless, often in interaction with RNA and DNA. One of the more intriguing IDPs is the 526-residue RNA binding protein, Fused In Sarcoma (FUS), whose monomer conformations and condensates exhibit unusual behavior that are sensitive to solution conditions. By focussing principally on the N-terminus low complexity domain (FUS-LC comprising residues 1-214) and other truncations, we rationalize the findings of solid state NMR experiments, which show that FUS-LC adopts a non-polymorphic fibril (core-1) involving residues 39-95, flanked by fuzzy coats on both the N- and C- terminal ends. An alternate structure (core-2), whose free energy is comparable to core-1, emerges only in the truncated construct (residues 110-214). Both core-1 and core-2 fibrils are stabilized by a Tyrosine ladder as well as hydrophilic interactions. The morphologies (gels, fibrils, and glass-like behavior) adopted by FUS seem to vary greatly, depending on the experimental conditions. The effect of phosphorylation is site specific and affects the stability of the fibril depending on the sites that are phosphorylated. Many of the peculiarities associated with FUS may also be shared by other IDPs, such as TDP43 and hnRNPA2. We outline a number of problems for which there is no clear molecular understanding.

Energy landscapes of Aβ monomers are sculpted in accordance with Ostwald's rule of stages

The transition from a disordered to an assembly-competent monomeric state (N*) in amyloidogenic sequences is a crucial event in the aggregation cascade. Using a well-calibrated model for intrinsically disordered proteins (IDPs), we show that the N* states, which bear considerable resemblance to the polymorphic fibril structures found in experiments, not only appear as excitations in the free energy landscapes of Aβ40 and Aβ42, but also initiate the aggregation cascade. For Aβ42, the transitions to the different N* states are in accord with Ostwald's rule of stages, with the least stable structures forming ahead of thermodynamically favored ones. The Aβ40 and Aβ42 monomer landscapes exhibit different extents of local frustration, which we show have profound implications in dictating subsequent self-assembly. Using kinetic transition networks, we illustrate that the most favored dimerization routes proceed via N* states. We argue that Ostwald's rule also holds for the aggregation of fused in sarcoma and polyglutamine proteins.

Efficient calculation of the free energy for protein partitioning using restraining potentials

An approach for the efficient simulation of phase-separated lipid bilayers, for use in the calculation of equilibrium free energies of partitioning between lipid domains, is proposed. The methodology exploits restraint potentials and rectangular aspect ratios that enforce lipid phase separation, allowing for the simulation of smaller systems that approximately reproduce bulk behavior. The utility of this approach is demonstrated through the calculation of potentials of mean force for the translation of a transmembrane protein between lipid domains. The impact of the imposed restraints on lipid tail ordering and lipid packing are explored, providing insight into how restraints can best be employed to compute accurate free-energy surfaces. This approach should be useful in the accurate calculation of equilibrium partition coefficients for transmembrane protein partitioning in heterogeneous membranes, providing insight into the thermodynamic driving forces that control this fundamental biophysical phenomenon.

Dynamical Models of Chemical Exchange in Nuclear Magnetic Resonance Spectroscopy

Chemical exchange line-broadening is an important phenomenon in nuclear magnetic resonance (NMR) spectroscopy, in which a nuclear spin experiences more than one magnetic environment as a result of chemical or conformational changes of a molecule. The dynamic process of chemical exchange strongly affects the sensitivity and resolution of NMR experiments, and increasingly provides a powerful probe of the inter-conversion between chemical and conformational states of proteins, nucleic acids, and other biological macromolecules. A simple and often used theoretical description of chemical exchange in NMR spectroscopy is based on an idealized two-state jump model (the random-phase or telegraph signal). However, chemical exchange can also be represented as a barrier-crossing event that can be modeled using chemical reaction rate theory. The time scale of crossing is determined by the barrier height, the temperature, and the dissipation modeled as collisional or frictional damping. This tutorial explores the connection between the NMR theory of chemical exchange line-broadening and strong-collision models for chemical kinetics in statistical mechanics. Theoretical modeling and numerical simulation are used to map the rate of barrier-crossing dynamics of a particle on a potential energy surface to the chemical exchange relaxation rate constant. By developing explicit models for the exchange dynamics, the tutorial aims to elucidate the underlying dynamical processes that give rise to the rich phenomenology of chemical exchange observed in NMR spectroscopy. Software for generating and analyzing the numerical simulations is provided in the form of Python and Fortran source codes.

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.

Sequence Determines the Switch in the Fibril Forming Regions in the Low-Complexity FUS Protein and Its Variants

Residues spanning distinct regions of the low-complexity domain of the RNA-binding protein, Fused in Sarcoma (FUS-LC), form fibril structures with different core morphologies. Solid-state NMR experiments show that the 214-residue FUS-LC forms a fibril with an S-bend (core-1, residues 39-95), while the rest of the protein is disordered. In contrast, the fibrils of the C-terminal variant (FUS-LC-C; residues 111-214) have a U-bend topology (core-2, residues 112-150). Absence of the U-bend in FUS-LC implies that the two fibril cores do not coexist. Computer simulations show that these perplexing findings could be understood in terms of the population of sparsely populated fibril-like excited states in the monomer. The propensity to form core-1 is higher compared to core-2. We predict that core-2 forms only in truncated variants that do not contain the core-1 sequence. At the monomer level, sequence-dependent enthalpic effects determine the relative stabilities of the core-1 and core-2 topologies.

Finite-Size Effects and Optimal System Sizes in Simulations of Surfactant Micelle Self-Assembly

The spontaneous formation of micelles in aqueous solutions is governed by the amphipathic nature of surfactants and is practically interesting due to the regular use of micelles as membrane mimics, for the characterization of protein structure, and for drug design and delivery. We performed a systematic characterization of the finite-size effect observed in single-component dodecylphosphocholine (DPC) micelles with the coarse-grained MARTINI model. Of multiple coarse-grained solvent models investigated using large system sizes, the nonpolarizable solvent model was found to most accurately reproduce SANS spectra of 100 mM DPC in aqueous solution. We systematically investigated the finite-size effect at constant 100 mM concentration in 23 systems of sizes 40-150 DPC, confirming the finite-size effect to manifest as an oscillation in the mean micelle aggregation number about the thermodynamic aggregation number as the system size increases. The oscillations in aggregation number mostly diminish once the system supports the formation of three micelles. Similar oscillations were observed in the estimated critical micelle concentration with a mean value of 1.10 mM, which is in agreement with experiment to 0.1 mM. The accuracy of using a multiscale simulation approach to avoid finite-size effects in the micelle size distribution and SANS spectra using MARTINI and CHARMM36 was explored using multiple long time scale 500 DPC coarse-grained simulations, which were back-mapped to CHARMM36 all-atom systems. It was found that the MARTINI model generally occupies more volume than the all-atom model, leading to the formation of micelles that are of a reasonable radius of gyration but are smaller in aggregation number. The systematic characterization of the finite-size effect and exploration of multiscale modeling presented in this work provide guidance for the accurate modeling of micelles in simulations.

Addressing the Excessive Aggregation of Membrane Proteins in the MARTINI Model

The MARTINI model is a widely used coarse-grained force field popular for its capacity to represent a diverse array of complex biomolecules. However, efforts to simulate increasingly realistic models of membranes, involving complex lipid mixtures and multiple proteins, suggest that membrane protein aggregates are overstabilized by the MARTINI v2.2 force field. In this study, we address this shortcoming of the MARTINI model. We determined the free energy of dimerization of four transmembrane protein systems using the nonpolarizable MARTINI model. Comparison with experimental FRET-based estimates of the dimerization free energy was used to quantify the significant overstabilization of each protein homodimer studied. To improve the agreement between simulation and experiment, a single uniform scaling factor, α, was used to enhance the protein-lipid Lennard-Jones interaction. A value of α = 1.04-1.045 was found to provide the best fit to the dimerization free energies for the proteins studied while maintaining the specificity of contacts at the dimer interface. To further validate the modified force field, we performed a multiprotein simulation using both MARTINI v2.2 and the reparameterized MARTINI model. While the original MARTINI model predicts oligomerization of protein into a single aggregate, the reparameterized MARTINI model maintains a dynamic equilibrium between monomers and dimers as predicted by experimental studies. The proposed reparameterization is an alternative to the standard MARTINI model for use in simulations of realistic models of a biological membrane containing diverse lipids and proteins.

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.

Impact of Cholesterol Concentration and Lipid Phase on Structure and Fluctuation of Amyloid Precursor Protein

Elevated levels of cellular cholesterol have been identified as one factor contributing to the onset of Alzheimer's disease (AD). Specific interaction between cholesterol and the amyloid precursor protein (APP), investigated via NMR experiments and computational studies, has been proposed to play a critical role in the processing of APP by secretases and the biogenesis of amyloid-β (Aβ) protein. We present all-atom molecular dynamics simulations of the 40-residue congener of the C-terminal domain of APP, C9916-55 (C99), in cholesterol-enriched DMPC lipid bilayers. We investigated the effect of cholesterol concentration on the conformational ensemble of wild-type C99 and C99-cholesterol associations at the low pH of endosomal environments, at which residues E22 and D23 are neutral. C99 was also characterized in liquid ordered domains for Dutch (E22Q) and Iowa (D23N) Familial AD mutants at low pH and for the wild-type sequence using protonation states characteristic of neutral pH. Our results reproduce the equilibrium constant of past NMR characterizations of the C99-cholesterol interaction but are not consistent with the C99-cholesterol binding hypothesis. We find that the lifetimes of both DMPC and cholesterol complexed with C99 display a power-law distribution of residence lifetimes. Longer-lived C99-DMPC and C99-cholesterol complexes are primarily stabilized by salt bridges and hydrogen bonds of lysine amines to phosphate and hydroxyl groups. Nevertheless, specific interfaces for C99-cholesterol association which are not present for DMPC can be identified. Changes to C99-cholesterol interfaces are found to depend on C99 tilt angle and orientation of the juxtamembrane domain of C99 containing residues E22 and D23. These observations support a more nuanced view of the C99-cholesterol interaction than has previously been suggested. We propose that cholesterol modulates the conformation and activity of C99 and other small transmembrane proteins indirectly through induction of the liquid ordered phase and directly through hydrogen bonding. This suggests a critical role for membrane heterogeneity introduced by cholesterol in modulating the structural ensemble of C99 and the production of Aβ.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

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.

Structural Basis for Lipid Binding and Function by an Evolutionarily Conserved Protein, Serum Amyloid A

Serum amyloid A (SAA) is a plasma protein that transports lipids during inflammation. To explore SAA solution conformations and lipid-binding mechanism, we used hydrogen-deuterium exchange mass spectrometry, lipoprotein reconstitution, amino acid sequence analysis, and molecular dynamics simulations. Solution conformations of lipid-bound and lipid-free mSAA1 at pH~7.4 agreed in details with the crystal structures but also showed important differences. The results revealed that amphipathic α-helices h1 and h3 comprise a lipid-binding site that is partially pre-formed in solution, is stabilized upon binding lipids, and shows lipid-induced folding of h3. This site sequesters apolar ligands via a concave hydrophobic surface in SAA oligomers. The largely disordered/dynamic C-terminal region is conjectured to mediate the promiscuous binding of other ligands. The h1-h2 linker region is predicted to form an unexpected β-hairpin that may represent an early amyloidogenic intermediate. The results help establish structural underpinnings for understanding SAA interactions with its key functional ligands, its evolutional conservation, and its transition to amyloid.

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.

Sequence Effects on Size, Shape, and Structural Heterogeneity in Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) lack well-defined three-dimensional structures, thus challenging the archetypal notion of structure-function relationships. Determining the ensemble of conformations that IDPs explore under physiological conditions is the first step toward understanding their diverse cellular functions. Here, we quantitatively characterize the structural features of IDPs as a function of sequence and length using coarse-grained simulations. For diverse IDP sequences, with the number of residues ( NT) ranging from 20 to 441, our simulations not only reproduce the radii of gyration ( Rg) obtained from experiments, but also predict the full scattering intensity profiles in excellent agreement with small-angle X-ray scattering experiments. The Rg values are well-described by the standard Flory scaling law, Rg = Rg0 NTν, with ν ≈ 0.588, making it tempting to assert that IDPs behave as polymers in a good solvent. However, clustering analysis reveals that the menagerie of structures explored by IDPs is diverse, with the extent of heterogeneity being highly sequence-dependent, even though ensemble-averaged properties, such as the dependence of Rg on chain length, may suggest synthetic polymer-like behavior in a good solvent. For example, we show that for the highly charged Prothymosin-α, a substantial fraction of conformations is highly compact. Even if the sequence compositions are similar, as is the case for α-Synuclein and a truncated construct from the Tau protein, there are substantial differences in the conformational heterogeneity. Taken together, these observations imply that metrics based on net charge or related quantities alone cannot be used to anticipate the phases of IDPs, either in isolation or in complex with partner IDPs or RNA. Our work sets the stage for probing the interactions of IDPs with each other, with folded protein domains, or with partner RNAs, which are critical for describing the structures of stress granules and biomolecular condensates with important cellular functions.

Aerosol-OT Surfactant Forms Stable Reverse Micelles in Apolar Solvent in the Absence of Water

Normal micelle aggregates of amphiphilic surfactant in aqueous solvents are formed by a process of entropically driven self-assembly. The self-assembly of reverse micelles from amphiphilic surfactant in a nonpolar solvent in the presence of water is considered to be an enthalpically driven process. Although the formation of normal and reverse surfactant micelles has been well characterized in theory and experiment, the nature of dry micelle formation, from amphiphilic surfactant in a nonpolar solvent in the absence of water, is poorly understood. In this study, a theory of dry reverse micelle formation is developed. Variation in free energy during micelle assembly is derived for the specific case of aerosol-OT surfactant in isooctane solvent using atomistic molecular dynamics simulation analyzed using the energy representation method. The existence and thermodynamic stability of dry reverse micelles of limited size are confirmed. The abrupt occurrence of monodisperse aggregates is a clear signature of a critical micelle concentration, commonly observed in the formation of normal surfactant micelles. The morphology of large dry micelles provides insight into the nature of the thermodynamic driving forces stabilizing the formation of the surfactant aggregates. Overall, this study provides detailed insight into the structure and stability of dry reverse micelles assembly in a nonpolar solvent.

Design, synthesis, and biomedical applications of synthetic sulphated polysaccharides

This review summarizes the synthetic methods to sulphated polysaccharides, describes their compositional and structural diversity in regards to activity, and showcases their biomedical applications.

Enhanced sampling method in molecular simulations using genetic algorithm for biomolecular systems

We propose a molecular simulation method using genetic algorithm (GA) for biomolecular systems to obtain ensemble averages efficiently. In this method, we incorporate the genetic crossover, which is one of the operations of GA, to any simulation method such as conventional molecular dynamics (MD), Monte Carlo, and other simulation methods. The genetic crossover proposes candidate conformations by exchanging parts of conformations of a target molecule between a pair of conformations during the simulation. If the candidate conformations are accepted, the simulation resumes from the accepted ones. While conventional simulations are based on local update of conformations, the genetic crossover introduces global update of conformations. As an example of the present approach, we incorporated genetic crossover to MD simulations. We tested the validity of the method by calculating ensemble averages and the sampling efficiency by using two kinds of peptides, ALA3 and (AAQAA)₃ . The results show that for ALA3 system, the distribution probabilities of backbone dihedral angles are in good agreement with those of the conventional MD and replica-exchange MD simulations. In the case of (AAQAA)₃ system, our method showed lower structural correlation of α-helix structures than the other two methods and more flexibility in the backbone ψ angles than the conventional MD simulation. These results suggest that our method gives more efficient conformational sampling than conventional simulation methods based on local update of conformations. © 2018 Wiley Periodicals, Inc.

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.

Membrane-wrapped nanoparticles probe divergent roles of GM3 and phosphatidylserine in lipid-mediated viral entry pathways

Gold nanoparticles (NPs) wrapped in a membrane can be utilized as artificial virus nanoparticles (AVNs) that combine the large nonblinking or bleaching optical cross-section of the NP core with the biological surface properties and functionalities provided by a self-assembled lipid membrane. We used these hybrid nanomaterials to test the roles of monosialodihexosylganglioside (GM3) and phosphatidylserine (PS) for a lipid-mediated targeting of virus-containing compartments (VCCs) in macrophages. GM3-presenting AVNs bind to CD169 (Siglec-1)-expressing macrophages, but inclusion of PS in the GM3-containing AVN membrane decreases binding. Molecular dynamics simulations of the AVN membrane and experimental binding studies of CD169 to GM3-presenting AVNs reveal Na+-mediated interactions between GM3 and PS as a potential cause of the antagonistic action on binding by the two negatively charged lipids. GM3-functionalized AVNs with no or low PS content localize to tetherin+, CD9+ VCC in a membrane composition-depending fashion, but increasing amounts of PS in the AVN membrane redirect the NP to lysosomal compartments. Interestingly, this compartmentalization is highly GM3 specific. Even AVNs presenting the related monosialotetrahexosylganglioside (GM1) fail to achieve an accumulation in VCC. AVN localization to VCC was observed for AVN with gold NP core but not for liposomes, suggesting that NP sequestration into VCC has additional requirements beyond ligand (GM3)-receptor (CD169) recognition that are related to the physical properties of the NP core. Our results confirm AVN as a scalable platform for elucidating the mechanisms of lipid-mediated viral entry pathways and for selective intracellular targeting.

Molecular Insights into Human Hereditary Apolipoprotein A-I Amyloidosis Caused by the Glu34Lys Mutation

Hereditary apolipoprotein A-I (apoA-I) amyloidosis is a life-threatening incurable genetic disorder whose molecular underpinnings are unclear. In this disease, variant apoA-I, the major structural and functional protein of high-density lipoprotein, is released in a free form, undergoes an α-helix to intermolecular cross-β-sheet conversion along with a proteolytic cleavage, and is deposited as amyloid fibrils in various organs, which can cause organ damage and death. Glu34Lys is the only known charge inversion mutation in apoA-I that causes human amyloidosis. To elucidate the structural underpinnings of the amyloidogenic behavior of Glu34Lys apoA-I, we generated its recombinant globular N-terminal domain (residues 1-184) and compared the conformation and dynamics of its lipid-free form with those of two other naturally occurring apoA-I variants, Phe71Tyr (amyloidogenic) and Leu159Arg (non-amyloidogenic). All variants showed reduced structural stability and altered aromatic residue packing. The greatest decrease in stability was observed in the non-amyloidogenic variant, suggesting that amyloid formation is driven by local structural perturbations at sensitive sites. Molecular dynamics simulations revealed local helical unfolding and suggested that transient opening of the Trp72 side chain induced mutation-dependent structural perturbations in a sensitive region, including the major amyloid hot spot residues Leu14-Leu22. We posit that a shift from the "closed" to the "open" orientation of the Trp72 side chain modulates structural protection of amyloid hot spots, suggesting a previously unknown early step in the protein misfolding pathway.

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.

Influence of membrane lipid composition on the structure and activity of γ-secretase

Insight into the effect of bilayer lipid composition on the GS structural ensemble and function using multiscale computational modeling approaches.

Critical size dependence of domain formation observed in coarse-grained simulations of bilayers composed of ternary lipid mixtures

Model cellular membranes are known to form micro- and macroscale lipid domains dependent on molecular composition. The formation of macroscopic lipid domains by lipid mixtures has been the subject of many simulation investigations. We present a critical study of system size impact on lipid domain phase separation into liquid-ordered and liquid-disordered macroscale domains in ternary lipid mixtures. In the popular di-C16:0 PC:di-C18:2 PC:cholesterol at 35:35:30 ratio mixture, we find systems with a minimum of 1480 lipids to be necessary for the formation of macroscopic phase separated domains and systems of 10 000 lipids to achieve structurally converged conformations similar to the thermodynamic limit. To understand these results and predict the behavior of any mixture forming two phases, we develop and investigate an analytical Flory-Huggins model which is recursively validated using simulation and experimental data. We find that micro- and macroscale domains can coexist in ternary mixtures. Additionally, we analyze the distributions of specific lipid-lipid interactions in each phase, characterizing domain structures proposed based on past experimental studies. These findings offer guidance in selecting appropriate system sizes for the study of phase separations and provide new insights into the nature of domain structure for a popular ternary lipid mixture.

Cellular prion protein targets amyloid-β fibril ends via its C-terminal domain to prevent elongation

Oligomeric forms of the amyloid-β (Aβ) peptide are thought to represent the primary synaptotoxic species underlying the neurodegenerative changes seen in Alzheimer's disease. It has been proposed that the cellular prion protein (PrP^C) functions as a cell-surface receptor, which binds to Aβ oligomers and transduces their toxic effects. However, the molecular details of the PrP^C-Aβ interaction remain uncertain. Here, we investigated the effect of PrP^C on polymerization of Aβ under rigorously controlled conditions in which Aβ converts from a monomeric to a fibrillar state via a series of kinetically defined steps. We demonstrated that PrP^C specifically inhibited elongation of Aβ fibrils, most likely by binding to the ends of growing fibrils. Surprisingly, this inhibitory effect required the globular C-terminal domain of PrP^C, which has not been previously implicated in interactions with Aβ. Our results suggest that PrP^C recognizes structural features common to both Aβ oligomers and fibril ends and that this interaction could contribute to the neurotoxic effect of Aβ aggregates. Additionally, our results identify the C terminus of PrP^C as a new and potentially more druggable molecular target for treating Alzheimer's disease.

Characterizing the structural ensemble of γ-secretase using a multiscale molecular dynamics approach

γ-Secretase is an intramembrane-cleaving aspartyl protease that plays an essential role in the processing of a variety of integral membrane proteins. Its role in the ultimate cleavage step in the processing of amyloid precursor protein to form amyloid-β (Aβ) peptide makes it an important therapeutic target in Alzheimer's disease research. Significant recent advances have been made in structural studies of this critical membrane protein complex. However, details of the mechanism of activation of the enzyme complex remain unclear. Using a multiscale computational modeling approach, combining multiple coarse-grained microsecond dynamic trajectories with all-atom models, the structure and two conformational states of the γ-secretase complex were evaluated. The transition between enzymatic state 1 and state 2 is shown to critically depend on the protonation states of the key catalytic residues Asp257 and Asp385 in the active site domain. The active site formation, related to our γ-secretase state 2, is observed to involve a concerted movement of four transmembrane helices from the catalytic subunit, resulting in the required localization of the catalytic residues. Global analysis of the structural ensemble of the enzyme complex was used to identify collective fluctuations important to the mechanism of substrate recognition and demonstrate that the corresponding fluctuations observed were uncorrelated with structural changes associated with enzyme activation. Overall, this computational study provides essential insight into the role of structure and dynamics in the activation and function of γ-secretase.

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.

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.

On the use of mass scaling for stable and efficient simulated tempering with molecular dynamics

Simulated tempering (ST) is a generalized-ensemble algorithm that employs trajectories exploring a range of temperatures to effectively sample rugged energy landscapes. When implemented using the molecular dynamics method, ST can require the use of short time steps for ensuring the stability of trajectories at high temperatures. To address this shortcoming, a mass-scaling ST (MSST) method is presented in which the particle mass is scaled in proportion to the temperature. Mass scaling in the MSST method leads to velocity distributions that are independent of temperature and eliminates the need for velocity scaling after the accepted temperature updates that are required in conventional ST simulations. The homogeneity in time scales with changing temperature improves the stability of simulations and allows for the use of longer time steps at high temperatures. As a result, the MSST is found to be more efficient than the standard ST method, particularly for cases in which a large temperature range is employed. © 2016 Wiley Periodicals, Inc.

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.

Exploring the role of hydration and confinement in the aggregation of amyloidogenic peptides Aβ(16-22) and Sup35(7-13) in AOT reverse micelles

Knowledge of how intermolecular interactions of amyloid-forming proteins cause protein aggregation and how those interactions are affected by sequence and solution conditions is essential to our understanding of the onset of many degenerative diseases. Of particular interest is the aggregation of the amyloid-β (Aβ) peptide, linked to Alzheimer's disease, and the aggregation of the Sup35 yeast prion peptide, which resembles the mammalian prion protein linked to spongiform encephalopathies. To facilitate the study of these important peptides, experimentalists have identified small peptide congeners of the full-length proteins that exhibit amyloidogenic behavior, including the KLVFFAE sub-sequence, Aβ16-22, and the GNNQQNY subsequence, Sup357-13. In this study, molecular dynamics simulations were used to examine these peptide fragments encapsulated in reverse micelles (RMs) in order to identify the fundamental principles that govern how sequence and solution environment influence peptide aggregation. Aβ16-22 and Sup357-13 are observed to organize into anti-parallel and parallel β-sheet arrangements. Confinement in the sodium bis(2-ethylhexyl) sulfosuccinate (AOT) reverse micelles is shown to stabilize extended peptide conformations and enhance peptide aggregation. Substantial fluctuations in the reverse micelle shape are observed, in agreement with earlier studies. Shape fluctuations are found to facilitate peptide solvation through interactions between the peptide and AOT surfactant, including direct interaction between non-polar peptide residues and the aliphatic surfactant tails. Computed amide I IR spectra are compared with experimental spectra and found to reflect changes in the peptide structures induced by confinement in the RM environment. Furthermore, examination of the rotational anisotropy decay of water in the RM demonstrates that the water dynamics are sensitive to the presence of peptide as well as the peptide sequence. Overall, our results demonstrate that the RM is a complex confining environment where substantial direct interaction between the surfactant and peptides plays an important role in determining the resulting ensemble of peptide conformations. By extension the results suggest that similarly complex sequence-dependent interactions may determine conformational ensembles of amyloid-forming peptides in a cellular environment.

Role of Charge and Solvation in the Structure and Dynamics of Alanine-Rich Peptide AKA2 in AOT Reverse Micelles

The propensity of peptides to form α-helices has been intensely studied using theory, computation, and experiment. Important model peptides for the study of the coil-to-helix transition have been alanine-lysine (AKA) peptides in which the lysine residues are placed on opposite sides of the helix avoiding charge repulsion while enhancing solubility. In this study, the effects of capped versus zwitterionic peptide termini on the secondary structure of alanine-rich peptides in reverse micelles are explored. The reverse micelles are found to undergo substantial shape fluctuations, a property observed in previous studies of AOT reverse micelles in the absence of solvated peptide. The peptides are observed to interact with water, as well as the AOT surfactant, including interactions between the nonpolar residues and the aliphatic surfactant tails. Computation of IR spectra for the amide I band of the peptide allows for direct comparison with experimental spectra. The results demonstrate that capped AKA2 peptides form more stable α helices than zwitterionic AKA2 peptides in reverse micelles. The rotational anisotropy decay of water is found to be distinctly different in the presence or absence of peptide within the reverse micelle, suggesting that the introduction of peptide significantly alters the number of free waters within the reverse micelle nanopool. However, neither the nature of the peptide termini (capped or charged) nor the degree of peptide helicity is found to significantly alter the balance of interactions between the peptides and the environment. Observed changes in the degree of helicity in AKA2 peptides in bulk solution and in reverse micelle environments result from changes in peptide confinement and hydration as well as direct nonpolar and polar interactions with the water-surfactant interface.

Propensity to form amyloid fibrils is encoded as excitations in the free energy landscape of monomeric proteins

Protein aggregation, linked to many of diseases, is initiated when monomers access rogue conformations that are poised to form amyloid fibrils. We show, using simulations of src SH3 domain, that mechanical force enhances the population of the aggregation-prone (N(⁎)) states, which are rarely populated under force free native conditions but are encoded in the spectrum of native fluctuations. The folding phase diagrams of SH3 as a function of denaturant concentration ([C]), mechanical force (f), and temperature exhibit an apparent two-state behavior, without revealing the presence of the elusive N(⁎) states. Interestingly, the phase boundaries separating the folded and unfolded states at all [C] and f fall on a master curve, which can be quantitatively described using an analogy to superconductors in a magnetic field. The free energy profiles as a function of the molecular extension (R), which are accessible in pulling experiments, (R), reveal the presence of a native-like N(⁎) with a disordered solvent-exposed amino-terminal β-strand. The structure of the N(⁎) state is identical with that found in Fyn SH3 by NMR dispersion experiments. We show that the timescale for fibril formation can be estimated from the population of the N(⁎) state, determined by the free energy gap separating the native structure and the N(⁎) state, a finding that can be used to assess fibril forming tendencies of proteins. The structures of the N(⁎) state are used to show that oligomer formation and likely route to fibrils occur by a domain-swap mechanism in SH3 domain.

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.