RESUMO
The accurate calculation of equilibrium constants for protein-protein association is of fundamental importance to quantitative biology and remains an outstanding challenge for computational biophysics. Traditionally, equilibrium constants have been computed from one-dimensional free energy surfaces derived from sampling along a single collective variable. Importantly, recent advances in enhanced sampling methodology have facilitated the characterization of multidimensional free energy landscapes, often exposing multiple thermodynamically important minima missed by more restrictive sampling methods. A key to the effectiveness of this multidimensional sampling approach is the identification of collective variables that effectively define the configurational space of dissociated and associated states. Here we present the application of two machine learning methods for the unbiased determination of collective variables for enhanced sampling and analysis of protein-protein association. Our results both validate prior work, based on intuition derived collective variables, and demonstrate the effectiveness of the machine learning methods for the identification of collective variables for association reactions in complex biomolecular systems.
Assuntos
Aprendizado de Máquina , Termodinâmica , Multimerização Proteica , Simulação de Dinâmica Molecular , Proteínas de Membrana/química , Proteínas de Membrana/metabolismoRESUMO
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.
Assuntos
Doença de Alzheimer , Peptídeos beta-Amiloides , Colesterol , Microdomínios da Membrana , Doença de Alzheimer/metabolismo , Colesterol/metabolismo , Colesterol/química , Humanos , Peptídeos beta-Amiloides/metabolismo , Peptídeos beta-Amiloides/química , Microdomínios da Membrana/metabolismo , Animais , Endocitose/fisiologiaRESUMO
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.
Assuntos
Proteínas Intrinsicamente Desordenadas , Sarcoma , Humanos , Domínios Proteicos , Espectroscopia de Ressonância Magnética , Conformação Proteica , Fosforilação , Proteínas Intrinsicamente Desordenadas/química , Proteína FUS de Ligação a RNA/genética , Proteína FUS de Ligação a RNA/química , Proteína FUS de Ligação a RNA/metabolismoRESUMO
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.
Assuntos
Proteínas de Membrana , Polímeros , Simulação por ComputadorRESUMO
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.
Assuntos
Doença de Alzheimer , Secretases da Proteína Precursora do Amiloide , Humanos , Secretases da Proteína Precursora do Amiloide/metabolismo , Doença de Alzheimer/metabolismo , Proteínas de Membrana/metabolismo , Ácido Aspártico Endopeptidases/metabolismo , Precursor de Proteína beta-Amiloide/metabolismo , Proteína ADAM10/metabolismo , Peptídeos beta-Amiloides/metabolismoRESUMO
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.
Assuntos
Membrana Celular , Poliprenois , Transferases , Ligantes , Proteínas de Membrana , Fosfatos , Poliprenois/metabolismo , Transferases/química , Fosfatos de Poli-Isoprenil/química , Membrana Celular/química , Bactérias/química , Bactérias/citologiaRESUMO
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.
Assuntos
Peptídeos beta-Amiloides , Fragmentos de Peptídeos , Peptídeos beta-Amiloides/metabolismo , Fragmentos de Peptídeos/metabolismo , EntropiaRESUMO
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.
RESUMO
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.
Assuntos
Bicamadas Lipídicas , Proteínas de Membrana , Bicamadas Lipídicas/metabolismo , Termodinâmica , Simulação por Computador , Proteínas de Membrana/metabolismo , Membranas/metabolismoRESUMO
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.
Assuntos
Peptídeos beta-Amiloides , Precursor de Proteína beta-Amiloide , Precursor de Proteína beta-Amiloide/metabolismo , Dimerização , Peptídeos beta-Amiloides/metabolismo , Conformação Proteica em Folha beta , Domínios Proteicos , Secretases da Proteína Precursora do Amiloide/metabolismoRESUMO
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.
Assuntos
Proteínas de Membrana , Simulação de Dinâmica Molecular , Entropia , Proteínas de Membrana/química , Conformação Molecular , TermodinâmicaRESUMO
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.
RESUMO
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.
Assuntos
Amiloide/metabolismo , Proteína FUS de Ligação a RNA/metabolismo , Sequência de Aminoácidos , Amiloide/química , Humanos , Simulação de Dinâmica Molecular , Mutagênese , Ressonância Magnética Nuclear Biomolecular , Conformação Proteica , Domínios Proteicos , Proteína FUS de Ligação a RNA/química , Proteína FUS de Ligação a RNA/genéticaRESUMO
Lipid-coated noble metal nanoparticles (L-NPs) combine the biomimetic surface properties of a self-assembled lipid membrane with the plasmonic properties of a nanoparticle (NP) core. In this work, we investigate derivatives of cholesterol, which can be found in high concentrations in biological membranes, and other terpenoids, as tunable, synthetic platforms to functionalize L-NPs. Side chains of different length and polarity, with a terminal alkyne group as Raman label, are introduced into cholesterol and betulin frameworks. The synthesized tags are shown to coexist in two conformations in the lipid layer of the L-NPs, identified as "head-out" and "head-in" orientations, whose relative ratio is determined by their interactions with the lipid-water hydrogen-bonding network. The orientational dimorphism of the tags introduces orthogonal functionalities into the NP surface for selective targeting and plasmon-enhanced Raman sensing, which is utilized for the identification and Raman imaging of epidermal growth factor receptor-overexpressing cancer cells.
Assuntos
Lipídeos/química , Lipossomos/química , Nanopartículas Metálicas/química , Nanopartículas/química , Química Click , Bicamadas Lipídicas/química , Simulação de Dinâmica MolecularRESUMO
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.
Assuntos
Micelas , Tensoativos , Simulação por Computador , Solventes , TermodinâmicaRESUMO
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.
Assuntos
Glicoforinas/química , Proteínas de Membrana/química , Simulação de Dinâmica Molecular , Transferência Ressonante de Energia de Fluorescência , Bicamadas Lipídicas/química , Agregados Proteicos , TermodinâmicaRESUMO
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.