Search

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

The 99-residue C-terminal domain of amyloid precursor protein (APP-C99), precursor to amyloid beta (A

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.

Cholesterol is a ubiquitous component of mammalian cell membranes and affects membrane protein function. Although cholesterol-mediated formation of ordered membrane domains has been extensively studied, molecular-level structural information about cholesterol self-association has been absent. Here, we combine solid-state nuclear magnetic resonance (NMR) spectroscopy with all-atom molecular dynamics simulations to determine the oligomeric structure of cholesterol in phospholipid bilayers. Two-dimensional 13C-13C correlation spectra of differentially labeled cholesterol indicate that cholesterol self-associates in a face-to-face fashion at membrane concentrations from 17 to 44 mol %. 2D 13C and 19F spin-counting experiments allowed us to measure the average oligomeric number of these cholesterol clusters. At low cholesterol concentrations of ∼20%, the average cluster size is centered on dimers. At a high cholesterol concentration of 44%, which is representative of virus lipid envelopes and liquid-ordered domains of cell membranes, both dimers and tetramers are observed. The cholesterol dimers are found in both phase-separated membranes that contain sphingomyelin and in disordered and miscible membranes that are free of sphingomyelin. Molecular dynamics simulations support these experimental observations and moreover provide the lifetimes, stabilities, distributions, and structures of these nanoscopic cholesterol clusters. Taken together, these NMR and MD data strongly suggest that dimers are the basic structural unit of cholesterol in phospholipid bilayers. The direct observation of cholesterol dimers and tetramers provides a revised framework for studying cholesterol interactions with membrane proteins to regulate protein functions and for understanding the pathogenic role of cholesterol in diseases.

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 co-solubilized 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.

Normal micelle aggregates of amphiphilic surfactant in aqueous solvent are formed by a process of entropically driven self-assembly. The self-assembly of reverse micelles from amphiphilic surfactant in non-polar solvent in the presence of water is considered to be an enthalpically driven process. While 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 non-polar 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 AOT 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 occurence of monodisperse aggregates is a clear signature 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 non-polar solvent., Comment: 34 pages, 6 figure. In Ver2, GAP joined and re-plot figures, and we changed the title

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 non-polarizable 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 to 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, mostly diminishing once the system supports formation of three micelles. The accuracy of employing 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 timescale 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 provides guidance for the accurate modeling of micelles in simulations., 19 pages, 5 figures, 2 tables, Corresponding Author: John E. Straub

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.

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.

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, I) miscible liquid disordered bulk, II) phase separated, domain registered coexistence of liquid disordered and liquid ordered and domains, and III) phase separated, domain-anti-registered 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.

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 G37G38 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.

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 increasingly decorrelated to each other. Mean and variance of the transmembrane and G(37)G(38) 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. The N-terminus, which 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.

The p53-MDM2 complex is both a major target for cancer drug development and a valuable model system for computational predictions of protein-ligand binding. To investigate the accuracy of molecular simulations of MDM2 and its complex with p53, we performed a number of long (200 ns to 1 µs) explicit-solvent simulations using a range of force fields. We systematically compared nine popular force fields (AMBER ff03, ff12sb, ff14sb, ff99sb, ff99sb-ildn, ff99sb-ildn-nmr, ff99sb-ildn-phi, CHARMM22*, and CHARMM36) against experimental chemical shift data, and found similarly accurate results, with microsecond simulations achieving better agreement compared to 200-ns trajectories. Although the experimentally determined apo structure has a closed binding cleft, simulations in all force fields suggest the apo state of MDM2 is highly flexible, and able to sample holo-like conformations, consistent with a conformational selection model. Initial structuring of the MDM2 lid region, known to competitively bind the binding cleft, is also observed in long simulations. Taken together, these results show molecular simulations can accurately sample conformations relevant for ligand binding. We expect this study to inform future computational work on folding and binding of MDM2 ligands. Proteins 2015; 83:1665–1676. © 2015 Wiley Periodicals, Inc.

Under normal cellular conditions, the tumor suppressor protein p53 is kept at low levels in part due to ubiquitination by MDM2, a process initiated by binding of MDM2 to the intrinsically disordered transactivation domain (TAD) of p53. Many experimental and simulation studies suggest that disordered domains such as p53 TAD bind their targets nonspecifically before folding to a tightly associated conformation, but the microscopic details are unclear. Toward a detailed prediction of binding mechanisms, pathways, and rates, we have performed large-scale unbiased all-atom simulations of p53-MDM2 binding. Markov state models (MSMs) constructed from the trajectory data predict p53 TAD binding pathways and on-rates in good agreement with experiment. The MSM reveals that two key bound intermediates, each with a nonnative arrangement of hydrophobic residues in the MDM2 binding cleft, control the overall on-rate. Using microscopic rate information from the MSM, we parameterize a simple four-state kinetic model to 1) determine that induced-fit pathways dominate the binding flux over a large range of concentrations, and 2) predict how modulation of residual p53 helicity affects binding, in good agreement with experiment. These results suggest new ways in which microscopic models of peptide binding, coupled with simple few-state binding flux models, can be used to understand biological function in physiological contexts.

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.

Under normal cellular conditions, the tumor suppressor protein p53 is kept at a low levels in part due to ubiquitination by MDM2, a process initiated by binding of MDM2 to the intrinsically disordered transactivation domain (TAD) of p53. Although many experimental and simulation studies suggest that disordered domains such as p53 TAD bind their targets nonspecifically before folding to a tightly-associated conformation, the molecular details are unclear. Toward a detailed prediction of binding mechanism, pathways and rates, we have performed large-scale unbiased all-atom simulations of p53-MDM2 binding. Markov State Models (MSMs) constructed from the trajectory data predict p53 TAD peptide binding pathways and on-rates in good agreement with experiment. The MSM reveals that two key bound intermediates, each with a non-native arrangement of hydrophobic residues in the MDM2 binding cleft, control the overall on-rate. Using microscopic rate information from the MSM, we parameterize a simple four-state kinetic model to (1) determine that induced-fit pathways dominate the binding flux over a large range of concentrations, and (2) predict how modulation of residual p53 helicity affects binding, in good agreement with experiment. These results suggest new ways in which microscopic models of bound-state ensembles can be used to understand biological function on a macroscopic scale.AUTHOR SUMMARYMany cell signaling pathways involve protein-protein interactions in which an intrinsically disordered peptide folds upon binding its target. Determining the molecular mechanisms that control these binding rates is important for understanding how such systems are regulated. In this paper, we show how extensive all-atom simulations combined with kinetic network models provide a detailed mechanistic understanding of how tumor suppressor protein p53 binds to MDM2, an important target of new cancer therapeutics. A simple four-state model parameterized from the simulations shows a binding-then-folding mechanism, and recapitulates experiments in which residual helicity boosts binding. This work goes beyond previous simulations of small-molecule binding, to achieve pathways and binding rates for a large peptide, in good agreement with experiment.

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.

MDM2 is a negative regulator of p53 activity and an important target for cancer therapeutics. The N-terminal lid region of MDM2 modulates interactions with p53 via competition for its binding cleft, exchanging slowly between docked and undocked conformations in the absence of p53. To better understand these dynamics, we constructed Markov State Models (MSMs) from large collections of unbiased simulation trajectories of apo-MDM2, and find strong evidence for diffuse, yet two-state folding and binding of the N-terminal region to the p53 receptor site. The MSM also identifies holo-like receptor conformations highly suitable for computational docking, despite initiating trajectories from closed-cleft receptor structures unsuitable for docking. Fixed-anchor docking studies using a test set of high-affinity small molecules and peptides show simulated receptor ensembles achieve docking successes comparable to cross-docking studies using crystal structures of receptors bound by alternative ligands. For p53, the best-scoring receptor structures have the N-terminal region lid region bound in a helical conformation mimicking the bound structure of p53, suggesting lid region association induces receptor conformations suitable for binding. These results suggest that MD + MSM approaches can sample binding-competent receptor conformations suitable for computational peptidomimetic design, and that inclusion of disordered regions may be essential to capturing the correct receptor dynamics.

Reverse micelles (RMs) are recognized as a paradigm of molecular self-assembly and used in a variety of applications, such as chemical synthesis and molecular structure refinement. Nevertheless, many fundamental properties including their equilibrium size distribution, internal structure, and mechanism of self-assembly remain poorly understood. To provide an enhanced microscopic understanding of the assembly process and resulting structural distribution, we perform multiple nonequilibrium molecular dynamics simulations of dioctyl sulfosuccinate sodium salt (AOT) RM assembly, quantifying RM size, water core structure, and dynamics. Rapid assembly of smaller RM from a random mixture is observed to establish a constant AOT water loading within a nanosecond consistent with a diffusion-adsorption mechanism validated through the Monte-Carlo simulation of a model system. The structure of RM water cores and RM molecular volume during RM assembly is characterized during the AOT assembly process. A moment-closure equation is developed from a novel master equation model to elucidate the elementary events underlying the AOT self-assembly process. The resulting kinetic model is used to explore the role of monomer addition and dissociation, RM association and dissociation, and RM collision-induced exchange, all dependent on average RM size, which provides fundamental insight regarding the mechanisms and time scales for AOT RM self-assembly. The nascent dynamics that rapidly establish water loading, intermediate time scales of RM fusion, and longer time scale dynamics of inter-RM exchange essential in establishing the equilibrium condition are quantified through these kinetic models. Overall, this work provides insight into AOT RM self-assembly and provides a general theoretical framework for the analysis of the molecular self-assembly dynamics and mechanism.

The p53-MDM2 complex is both a major target for cancer drug development and a valuable model system for computational predictions of protein-ligand binding. To investigate the accuracy of molecular simulations of MDM2 and its complex with p53, we performed a number of long (200 ns to 1 µs) explicit-solvent simulations using a range of force fields. We systematically compared nine popular force fields (AMBER ff03, ff12sb, ff14sb, ff99sb, ff99sb-ildn, ff99sb-ildn-nmr, ff99sb-ildn-phi, CHARMM22*, and CHARMM36) against experimental chemical shift data, and found similarly accurate results, with microsecond simulations achieving better agreement compared to 200-ns trajectories. Although the experimentally determined apo structure has a closed binding cleft, simulations in all force fields suggest the apo state of MDM2 is highly flexible, and able to sample holo-like conformations, consistent with a conformational selection model. Initial structuring of the MDM2 lid region, known to competitively bind the binding cleft, is also observed in long simulations. Taken together, these results show molecular simulations can accurately sample conformations relevant for ligand binding. We expect this study to inform future computational work on folding and binding of MDM2 ligands.

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.