1
|
Grazioli G, Tao A, Bhatia I, Regan P. Genetic Algorithm for Automated Parameterization of Network Hamiltonian Models of Amyloid Fibril Formation. J Phys Chem B 2024; 128:1854-1865. [PMID: 38359362 PMCID: PMC10910512 DOI: 10.1021/acs.jpcb.3c07322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2023] [Revised: 01/07/2024] [Accepted: 02/05/2024] [Indexed: 02/17/2024]
Abstract
The time scales of long-time atomistic molecular dynamics simulations are typically reported in microseconds, while the time scales for experiments studying the kinetics of amyloid fibril formation are typically reported in minutes or hours. This time scale deficit of roughly 9 orders of magnitude presents a major challenge in the design of computer simulation methods for studying protein aggregation events. Coarse-grained molecular simulations offer a computationally tractable path forward for exploring the molecular mechanism driving the formation of these structures, which are implicated in diseases such as Alzheimer's, Parkinson's, and type-II diabetes. Network Hamiltonian models of aggregation are centered around a Hamiltonian function that returns the total energy of a system of aggregating proteins, given the graph structure of the system as an input. In the graph, or network, representation of the system, each protein molecule is represented as a node, and noncovalent bonds between proteins are represented as edges. The parameter, i.e., a set of coefficients that determine the degree to which each topological degree of freedom is favored or disfavored, must be determined for each network Hamiltonian model, and is a well-known technical challenge. The methodology is first demonstrated by beginning with an initial set of randomly parametrized models of low fibril fraction (<5% fibrillar), and evolving to subsequent generations of models, ultimately leading to high fibril fraction models (>70% fibrillar). The methodology is also demonstrated by applying it to optimizing previously published network Hamiltonian models for the 5 key amyloid fibril topologies that have been reported in the Protein Data Bank (PDB). The models generated by the AI produced fibril fractions that surpass previously published fibril fractions in 3 of 5 cases, including the most naturally abundant amyloid fibril topology, the 1,2 2-ribbon, which features a steric zipper. The authors also aim to encourage more widespread use of the network Hamiltonian methodology for fitting a wide variety of self-assembling systems by releasing a free open-source implementation of the genetic algorithm introduced here.
Collapse
Affiliation(s)
- Gianmarc Grazioli
- Department of Chemistry, San
José State University, San Jose, California 95192, United States
| | - Andy Tao
- Department of Chemistry, San
José State University, San Jose, California 95192, United States
| | - Inika Bhatia
- Department of Chemistry, San
José State University, San Jose, California 95192, United States
| | - Patrick Regan
- Department of Chemistry, San
José State University, San Jose, California 95192, United States
| |
Collapse
|
2
|
Tolstova AP, Makarov AA, Adzhubei AA. Structure Comparison of Beta Amyloid Peptide Aβ 1-42 Isoforms. Molecular Dynamics Modeling. J Chem Inf Model 2024; 64:918-932. [PMID: 38241093 DOI: 10.1021/acs.jcim.3c01624] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2024]
Abstract
Beta amyloid peptide Aβ 1-42 (Aβ42) has a unique dual role in the human organism, as both the peptide with an important physiological function and one of the most toxic biological compounds provoking Alzheimer's disease (AD). There are several known Aβ42 isoforms that we discuss here that are highly neurotoxic and lead to the early onset of AD. Aβ42 is an intrinsically disordered protein with no experimentally solved structure under physiological conditions. The objective of this research was to establish the appropriate molecular dynamics (MD) methodology and model a uniform set of structures for the Aβ42 isoforms that form the core of this study. For that purpose, force field selection and verification including convergence testing for MD simulations was made. Replica exchange MD and conventional MD modeling of several Aβ42 and Aβ16 isoforms that have neurotoxic and amyloidogenic effects impacting the severity of Alzheimer's disease were carried out with the optimal force field and solvent parameters. A standardized ensemble of structures for the Aβ42 and Aβ16 isoforms covering 30-50% of the conformational ensembles extracted from the free energy minima was calculated from MD trajectories. The resulting data set of modeled structures includes Aβ42 wild type, isoD7, pS8, D7H, and H6R-Aβ42 and Aβ16 wild type, isoD7, pS8, D7H, and H6R-Aβ16. The representative structures are given in the Supporting Information; they are open for public access. In the study, we also evaluated the differences between the structures of Aβ42 isoforms and speculate on their possible relevance to the known functions. Utilizing several representative structures for a single disordered protein for docking, with their subsequent averaging by conformations, would markedly increase the reliability of docking results.
Collapse
Affiliation(s)
- Anna P Tolstova
- Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, Moscow 119991, Russia
| | - Alexander A Makarov
- Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, Moscow 119991, Russia
| | - Alexei A Adzhubei
- Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, Moscow 119991, Russia
- Washington University School of Medicine and Health Sciences, Washington 20052, D.C., United States
| |
Collapse
|
3
|
Macke AC, Stump JE, Kelly MS, Rowley J, Herath V, Mullen S, Dima RI. Searching for Structure: Characterizing the Protein Conformational Landscape with Clustering-Based Algorithms. J Chem Inf Model 2024; 64:470-482. [PMID: 38173388 DOI: 10.1021/acs.jcim.3c01511] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
The identification and characterization of the main conformations from a protein population are a challenging and inherently high-dimensional problem. Here, we evaluate the performance of the Secondary sTructural Ensembles with machine LeArning (StELa) double-clustering method, which clusters protein structures based on the relationship between the φ and ψ dihedral angles in a protein backbone and the secondary structure of the protein, thus focusing on the local properties of protein structures. The classification of states as vectors composed of the clusters' indices arising naturally from the Ramachandran plot is followed by the hierarchical clustering of the vectors to allow for the identification of the main features of the corresponding free energy landscape (FEL). We compare the performance of StELa with the established root-mean-squared-deviation (RMSD)-based clustering algorithm, which focuses on global properties of protein structures and with Combinatorial Averaged Transient Structure (CATS), the combinatorial averaged transient structure clustering method based on distributions of the φ and ψ dihedral angle coordinates. Using ensembles of conformations from molecular dynamics simulations of intrinsically disordered proteins (IDPs) of various lengths (tau protein fragments) or short fragments from a globular protein, we show that StELa is the clustering method that identifies many of the minima and relevant energy states around the minima from the corresponding FELs. In contrast, the RMSD-based algorithm yields a large number of clusters that usually cover most of the FEL, thus being unable to distinguish between states, while CATS does not sample well the FELs for long IDPs and fragments from globular proteins.
Collapse
Affiliation(s)
- Amanda C Macke
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| | - Jacob E Stump
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| | - Maria S Kelly
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| | - Jamie Rowley
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| | - Vageesha Herath
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| | - Sarah Mullen
- Department of Chemistry, The College of Wooster, Wooster, Ohio 44691, United States
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Ruxandra I Dima
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| |
Collapse
|
4
|
Lalmansingh JM, Keeley AT, Ruff KM, Pappu RV, Holehouse AS. SOURSOP: A Python Package for the Analysis of Simulations of Intrinsically Disordered Proteins. J Chem Theory Comput 2023; 19:5609-5620. [PMID: 37463458 PMCID: PMC11188088 DOI: 10.1021/acs.jctc.3c00190] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2023]
Abstract
Conformational heterogeneity is a defining hallmark of intrinsically disordered proteins and protein regions (IDRs). The functions of IDRs and the emergent cellular phenotypes they control are associated with sequence-specific conformational ensembles. Simulations of conformational ensembles that are based on atomistic and coarse-grained models are routinely used to uncover the sequence-specific interactions that may contribute to IDR functions. These simulations are performed either independently or in conjunction with data from experiments. Functionally relevant features of IDRs can span a range of length scales. Extracting these features requires analysis routines that quantify a range of properties. Here, we describe a new analysis suite simulation analysis of unfolded regions of proteins (SOURSOP), an object-oriented and open-source toolkit designed for the analysis of simulated conformational ensembles of IDRs. SOURSOP implements several analysis routines motivated by principles in polymer physics, offering a unique collection of simple-to-use functions to characterize IDR ensembles. As an extendable framework, SOURSOP supports the development and implementation of new analysis routines that can be easily packaged and shared.
Collapse
Affiliation(s)
- Jared M. Lalmansingh
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130, USA
- Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO 63130, USA
| | - Alex T. Keeley
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130, USA
- Department of Chemistry, University of Illinois Urbana-Champaign, Urbana-Champaign, IL 61801, USA
| | - Kiersten M. Ruff
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130, USA
- Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO 63130, USA
| | - Rohit V. Pappu
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130, USA
- Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO 63130, USA
| | - Alex S. Holehouse
- Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO 63130, USA
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO, 63110, USA
| |
Collapse
|
5
|
Zhang P, Yang W. Toward a general neural network force field for protein simulations: Refining the intramolecular interaction in protein. J Chem Phys 2023; 159:024118. [PMID: 37431910 PMCID: PMC10481389 DOI: 10.1063/5.0142280] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2023] [Accepted: 06/22/2023] [Indexed: 07/12/2023] Open
Abstract
Molecular dynamics (MD) is an extremely powerful, highly effective, and widely used approach to understanding the nature of chemical processes in atomic details for proteins. The accuracy of results from MD simulations is highly dependent on force fields. Currently, molecular mechanical (MM) force fields are mainly utilized in MD simulations because of their low computational cost. Quantum mechanical (QM) calculation has high accuracy, but it is exceedingly time consuming for protein simulations. Machine learning (ML) provides the capability for generating accurate potential at the QM level without increasing much computational effort for specific systems that can be studied at the QM level. However, the construction of general machine learned force fields, needed for broad applications and large and complex systems, is still challenging. Here, general and transferable neural network (NN) force fields based on CHARMM force fields, named CHARMM-NN, are constructed for proteins by training NN models on 27 fragments partitioned from the residue-based systematic molecular fragmentation (rSMF) method. The NN for each fragment is based on atom types and uses new input features that are similar to MM inputs, including bonds, angles, dihedrals, and non-bonded terms, which enhance the compatibility of CHARMM-NN to MM MD and enable the implementation of CHARMM-NN force fields in different MD programs. While the main part of the energy of the protein is based on rSMF and NN, the nonbonded interactions between the fragments and with water are taken from the CHARMM force field through mechanical embedding. The validations of the method for dipeptides on geometric data, relative potential energies, and structural reorganization energies demonstrate that the CHARMM-NN local minima on the potential energy surface are very accurate approximations to QM, showing the success of CHARMM-NN for bonded interactions. However, the MD simulations on peptides and proteins indicate that more accurate methods to represent protein-water interactions in fragments and non-bonded interactions between fragments should be considered in the future improvement of CHARMM-NN, which can increase the accuracy of approximation beyond the current mechanical embedding QM/MM level.
Collapse
Affiliation(s)
- Pan Zhang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
6
|
Zhang Y, Liu X, Chen J. Re-Balancing Replica Exchange with Solute Tempering for Sampling Dynamic Protein Conformations. J Chem Theory Comput 2023; 19:1602-1614. [PMID: 36791464 PMCID: PMC10795075 DOI: 10.1021/acs.jctc.2c01139] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
Abstract
Replica exchange with solute tempering (REST) is a highly effective variant of replica exchange for enhanced sampling in explicit solvent simulations of biomolecules. By scaling the Hamiltonian for a selected "solute" region of the system, REST effectively applies tempering only to the degrees of freedom of interest but not the rest of the system ("solvent"), allowing fewer replicas for covering the same temperature range. A key consideration of REST is how the solute-solvent interactions are scaled together with the solute-solute interactions. Here, we critically evaluate the performance of the latest REST2 protocol for sampling large-scale conformation fluctuations of intrinsically disordered proteins (IDPs). The results show that REST2 promotes artificial protein conformational collapse at high effective temperatures, which seems to be a designed feature originally to promote the sampling of reversible folding of small proteins. The collapse is particularly severe with larger IDPs, leading to replica segregation in the effective temperature space and hindering effective sampling of large-scale conformational changes. We propose that the scaling of the solute-solvent interactions can be treated as free parameters in REST, which can be tuned to control the solute conformational properties (e.g., chain expansion) at different effective temperatures and achieve more effective sampling. To this end, we derive a new REST3 protocol, where the strengths of the solute-solvent van der Waals interactions are recalibrated to reproduce the levels of protein chain expansion at high effective temperatures. The efficiency of REST3 is examined using two IDPs with nontrivial local and long-range structural features, including the p53 N-terminal domain and the kinase inducible transactivation domain of transcription factor CREB. The results suggest that REST3 leads to a much more efficient temperature random walk and improved sampling efficiency, which also further reduces the number of replicas required. Nonetheless, our analysis also reveals significant challenges of relying on tempering alone for sampling large-scale conformational fluctuations of disordered proteins. It is likely that more efficient sampling protocols will require incorporating more sophisticated Hamiltonian replica exchange schemes in addition to tempering.
Collapse
Affiliation(s)
- Yumeng Zhang
- Department of Chemistry, University of Massachusetts, Amherst, MA 01003, USA
| | - Xiaorong Liu
- Corresponding Authors: (XL), (JC), Phone: (413) 545-3386 (JC)
| | - Jianhan Chen
- Department of Chemistry, University of Massachusetts, Amherst, MA 01003, USA
| |
Collapse
|
7
|
Lalmansingh JM, Keeley AT, Ruff KM, Pappu RV, Holehouse AS. SOURSOP: A Python package for the analysis of simulations of intrinsically disordered proteins. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.02.16.528879. [PMID: 36824878 PMCID: PMC9949127 DOI: 10.1101/2023.02.16.528879] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/19/2023]
Abstract
Conformational heterogeneity is a defining hallmark of intrinsically disordered proteins and protein regions (IDRs). The functions of IDRs and the emergent cellular phenotypes they control are associated with sequence-specific conformational ensembles. Simulations of conformational ensembles that are based on atomistic and coarse-grained models are routinely used to uncover the sequence-specific interactions that may contribute to IDR functions. These simulations are performed either independently or in conjunction with data from experiments. Functionally relevant features of IDRs can span a range of length scales. Extracting these features requires analysis routines that quantify a range of properties. Here, we describe a new analysis suite SOURSOP, an object-oriented and open-source toolkit designed for the analysis of simulated conformational ensembles of IDRs. SOURSOP implements several analysis routines motivated by principles in polymer physics, offering a unique collection of simple-to-use functions to characterize IDR ensembles. As an extendable framework, SOURSOP supports the development and implementation of new analysis routines that can be easily packaged and shared.
Collapse
|
8
|
Depraz Depland A, Stroganova I, Wootton CA, Rijs AM. Developments in Trapped Ion Mobility Mass Spectrometry to Probe the Early Stages of Peptide Aggregation. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2023; 34:193-204. [PMID: 36633834 PMCID: PMC9896548 DOI: 10.1021/jasms.2c00253] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/07/2022] [Revised: 11/15/2022] [Accepted: 12/16/2022] [Indexed: 06/17/2023]
Abstract
Ion mobility mass spectrometry (IM-MS) has proven to be an excellent method to characterize the structure of amyloidogenic protein and peptide aggregates, which are formed in coincidence with the development of neurodegenerative diseases. However, it remains a challenge to obtain detailed structural information on all conformational intermediates, originating from the early onset of those pathologies, due to their complex and heterogeneous environment. One way to enhance the insights and the identification of these early stage oligomers is by employing high resolution ion mobility mass spectrometry experiments. This would allow us to enhance the mobility separation and MS characterization. Trapped ion mobility spectrometry (TIMS) is an ion mobility technique known for its inherently high resolution and has successfully been applied to the analysis of protein conformations among others. To obtain conformational information on fragile peptide aggregates, the instrumental parameters of the TIMS-Quadrupole-Time-of-Flight mass spectrometer (TIMS-qToF-MS) have to be optimized to allow the study of intact aggregates and ensure their transmission toward the detector. Here, we investigate the suitability and application of TIMS to probe the aggregation process, targeting the well-characterized M307-N319 peptide segment of the TDP-43 protein, which is involved in the development of amyotrophic lateral sclerosis. By studying the influence of key parameters over the full mass spectrometer, such as source temperature, applied voltages or RFs among others, we demonstrate that by using an optimized instrumental method TIMS can be used to probe peptide aggregation.
Collapse
Affiliation(s)
- Agathe Depraz Depland
- Division
of Bioanalytical Chemistry, Amsterdam Institute of Molecular and Life
Sciences, Vrije Universiteit Amsterdam, De Boelelaan 1105, 1081 HV Amsterdam, The Netherlands
| | - Iuliia Stroganova
- Division
of Bioanalytical Chemistry, Amsterdam Institute of Molecular and Life
Sciences, Vrije Universiteit Amsterdam, De Boelelaan 1105, 1081 HV Amsterdam, The Netherlands
| | | | - Anouk M. Rijs
- Division
of Bioanalytical Chemistry, Amsterdam Institute of Molecular and Life
Sciences, Vrije Universiteit Amsterdam, De Boelelaan 1105, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
9
|
Zhang Y, Liu X, Chen J. Toward Accurate Coarse-Grained Simulations of Disordered Proteins and Their Dynamic Interactions. J Chem Inf Model 2022; 62:4523-4536. [PMID: 36083825 PMCID: PMC9910785 DOI: 10.1021/acs.jcim.2c00974] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Intrinsically disordered proteins (IDPs) play crucial roles in cellular regulatory networks and are now recognized to often remain highly dynamic even in specific interactions and assemblies. Accurate description of these dynamic interactions is extremely challenging using atomistic simulations because of the prohibitive computational cost. Efficient coarse-grained approaches could offer an effective solution to overcome this bottleneck if they could provide an accurate description of key local and global properties of IDPs in both unbound and bound states. The recently developed hybrid-resolution (HyRes) protein model has been shown to be capable of providing a semiquantitative description of the secondary structure propensities of IDPs. Here, we show that greatly improved description of global structures and transient interactions can be achieved by introducing a solvent-accessible surface area-based implicit solvent term followed by reoptimization of effective interaction strengths. The new model, termed HyRes II, can semiquantitatively reproduce a wide range of local and global structural properties of a set of IDPs of various lengths and complexities. It can also distinguish the level of compaction between folded proteins and IDPs. In particular, applied to the disordered N-terminal transactivation domain (TAD) of tumor suppressor p53, HyRes II is able to recapitulate various nontrivial structural properties compared to experimental results, some of them to a level of accuracy that is almost comparable to results from atomistic explicit solvent simulations. Furthermore, we demonstrate that HyRes II can be used to simulate the dynamic interactions of TAD with the DNA-binding domain of p53, generating structural ensembles that are highly consistent with existing NMR data. We anticipate that HyRes II will provide an efficient and relatively reliable tool toward accurate coarse-grained simulations of dynamic protein interactions.
Collapse
Affiliation(s)
- Yumeng Zhang
- Department of Chemistry, University of Massachusetts, Amherst, MA 01003, USA
| | - Xiaorong Liu
- Department of Chemistry, University of Massachusetts, Amherst, MA 01003, USA
| | - Jianhan Chen
- Department of Chemistry, University of Massachusetts, Amherst, MA 01003, USA
| |
Collapse
|
10
|
Latham AP, Zhang B. Unifying coarse-grained force fields for folded and disordered proteins. Curr Opin Struct Biol 2022; 72:63-70. [PMID: 34536913 PMCID: PMC9057422 DOI: 10.1016/j.sbi.2021.08.006] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2021] [Revised: 08/08/2021] [Accepted: 08/17/2021] [Indexed: 12/22/2022]
Abstract
Liquid-liquid phase separation drives the formation of biological condensates that play essential roles in transcriptional regulation and signal sensing. Computational modeling could provide high-resolution structural characterizations of these condensates and help uncover physicochemical interactions that dictate their stability. However, many protein molecules involved in phase separation often contain multiple ordered domains connected with flexible, structureless linkers. Simulating such proteins necessitates force fields with consistent accuracy for both folded and disordered proteins. We provide a critical review of existing coarse-grained force fields for disordered proteins and highlight the challenges in their application to folded proteins. After discussing existing algorithms for force field parameterization, we propose an optimization strategy that should lead to computer models with improved transferability across protein types.
Collapse
Affiliation(s)
- Andrew P Latham
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA.
| |
Collapse
|
11
|
Statistical potentials from the Gaussian scaling behaviour of chain fragments buried within protein globules. PLoS One 2022; 17:e0254969. [PMID: 35085247 PMCID: PMC8794220 DOI: 10.1371/journal.pone.0254969] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2021] [Accepted: 10/28/2021] [Indexed: 11/19/2022] Open
Abstract
Knowledge-based approaches use the statistics collected from protein data-bank structures to estimate effective interaction potentials between amino acid pairs. Empirical relations are typically employed that are based on the crucial choice of a reference state associated to the null interaction case. Despite their significant effectiveness, the physical interpretation of knowledge-based potentials has been repeatedly questioned, with no consensus on the choice of the reference state. Here we use the fact that the Flory theorem, originally derived for chains in a dense polymer melt, holds also for chain fragments within the core of globular proteins, if the average over buried fragments collected from different non-redundant native structures is considered. After verifying that the ensuing Gaussian statistics, a hallmark of effectively non-interacting polymer chains, holds for a wide range of fragment lengths, although with significant deviations at short spatial scales, we use it to define a ‘bona fide’ reference state. Notably, despite the latter does depend on fragment length, deviations from it do not. This allows to estimate an effective interaction potential which is not biased by the presence of correlations due to the connectivity of the protein chain. We show how different sequence-independent effective statistical potentials can be derived using this approach by coarse-graining the protein representation at varying levels. The possibility of defining sequence-dependent potentials is explored.
Collapse
|
12
|
Nguyen PH, Derreumaux P. Computer Simulations Aimed at Exploring Protein Aggregation and Dissociation. Methods Mol Biol 2022; 2340:175-196. [PMID: 35167075 DOI: 10.1007/978-1-0716-1546-1_9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.
Collapse
Affiliation(s)
- Phuong H Nguyen
- Laboratoire de Biochimie Théorique, UPR 9080, CNRS, Université de Paris, Paris, France
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, UPR 9080, CNRS, Université de Paris, Paris, France.
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris, France.
| |
Collapse
|
13
|
Nguyen PH, Tufféry P, Derreumaux P. Dynamics of Amyloid Formation from Simplified Representation to Atomistic Simulations. Methods Mol Biol 2022; 2405:95-113. [PMID: 35298810 DOI: 10.1007/978-1-0716-1855-4_5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Amyloid fibril formation is an intrinsic property of short peptides, non-disease proteins, and proteins associated with neurodegenerative diseases. Aggregates of the Aβ and tau proteins, the α-synuclein protein, and the prion protein are observed in the brain of Alzheimer's, Parkinson's, and prion disease patients, respectively. Due to the transient short-range and long-range interactions of all species and their high aggregation propensities, the conformational ensemble of these devastating proteins, the exception being for the monomeric prion protein, remains elusive by standard structural biology methods in bulk solution and in lipid membranes. To overcome these limitations, an increasing number of simulations using different sampling methods and protein models have been performed. In this chapter, we first review our main contributions to the field of amyloid protein simulations aimed at understanding the early aggregation steps of short linear amyloid peptides, the conformational ensemble of the Aβ40/42 dimers in bulk solution, and the stability of Aβ aggregates in lipid membrane models. Then we focus on our studies on the interactions of amyloid peptides/inhibitors to prevent aggregation, and long amyloid sequences, including new results on a monomeric tau construct.
Collapse
Affiliation(s)
- Phuong Hoang Nguyen
- Laboratoire de Biochimie Théorique, CNRS, Université de Paris, UPR 9080, Paris, France
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | - Pierre Tufféry
- Université de Paris, BFA, UMR 8251, CNRS, ERL U1133, Inserm, RPBS, Paris, France
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, CNRS, Université de Paris, UPR 9080, Paris, France.
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris, France.
| |
Collapse
|
14
|
Kognole AA, Lee J, Park SJ, Jo S, Chatterjee P, Lemkul JA, Huang J, MacKerell AD, Im W. CHARMM-GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field. J Comput Chem 2021; 43:359-375. [PMID: 34874077 DOI: 10.1002/jcc.26795] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2021] [Revised: 11/10/2021] [Accepted: 11/25/2021] [Indexed: 12/18/2022]
Abstract
Explicit treatment of electronic polarizability in empirical force fields (FFs) represents an extension over a traditional additive or pairwise FF and provides a more realistic model of the variations in electronic structure in condensed phase, macromolecular simulations. To facilitate utilization of the polarizable FF based on the classical Drude oscillator model, Drude Prepper has been developed in CHARMM-GUI. Drude Prepper ingests additive CHARMM protein structures file (PSF) and pre-equilibrated coordinates in CHARMM, PDB, or NAMD format, from which the molecular components of the system are identified. These include all residues and patches connecting those residues along with water, ions, and other solute molecules. This information is then used to construct the Drude FF-based PSF using molecular generation capabilities in CHARMM, followed by minimization and equilibration. In addition, inputs are generated for molecular dynamics (MD) simulations using CHARMM, GROMACS, NAMD, and OpenMM. Validation of the Drude Prepper protocol and inputs is performed through conversion and MD simulations of various heterogeneous systems that include proteins, nucleic acids, lipids, polysaccharides, and atomic ions using the aforementioned simulation packages. Stable simulations are obtained in all studied systems, including 5 μs simulation of ubiquitin, verifying the integrity of the generated Drude PSFs. In addition, the ability of the Drude FF to model variations in electronic structure is shown through dipole moment analysis in selected systems. The capabilities and availability of Drude Prepper in CHARMM-GUI is anticipated to greatly facilitate the application of the Drude FF to a range of condensed phase, macromolecular systems.
Collapse
Affiliation(s)
- Abhishek A Kognole
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, USA
| | - Jumin Lee
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania, USA
| | - Sang-Jun Park
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania, USA
| | - Sunhwan Jo
- Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois, USA
| | - Payal Chatterjee
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, USA
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Tech, Blacksburg, Virginia, USA
| | - Jing Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, Zhejiang, Hangzhou, China
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, USA
| | - Wonpil Im
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania, USA
| |
Collapse
|
15
|
Choi B, Kim NH, Jin GY, Kim YS, Kim YH, Eom K. Sequence-dependent aggregation-prone conformations of islet amyloid polypeptide. Phys Chem Chem Phys 2021; 23:22532-22542. [PMID: 34590645 DOI: 10.1039/d1cp01061a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Amyloid proteins, which aggregate to form highly ordered structures, play a crucial role in various disease pathologies. Despite many previous studies on amyloid fibrils, which are an end product of protein aggregation, the structural characteristics of amyloid proteins in the early stage of aggregation and their related aggregation mechanism still remain elusive. The role of the amino acid sequence in the aggregation-prone structures of amyloid proteins at such a stage is not understood. Here, we have studied the sequence-dependent structural characteristics of islet amyloid polypeptide based on atomistic simulations and spectroscopic experiments. We show that the amino acid sequence determines non-bonded interactions that play a leading role in the formation of aggregation-prone conformations. Specifically, a single point mutation critically changes the population of aggregation-prone conformations, resulting in a change of the aggregation mechanism. Our simulation results were supported by experimental results suggesting that mutation affects the kinetics of aggregation and the structural characteristics of amyloid aggregates. Our study provides an insight into the role of sequence-dependent aggregation-prone conformations in the underlying mechanisms of amyloid aggregation.
Collapse
Affiliation(s)
- Bumjoon Choi
- Biomechanics Laboratory, College of Sport Science, Sungkyunkwan University (SKKU), Suwon 16419, Republic of Korea.
| | - Nam Hyeong Kim
- SKKU Advanced Institute of Nano Technology (SAINT), Sungkyunkwan University (SKKU), Suwon 16419, Republic of Korea.
| | - Geun Young Jin
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Yung Sam Kim
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Yong Ho Kim
- SKKU Advanced Institute of Nano Technology (SAINT), Sungkyunkwan University (SKKU), Suwon 16419, Republic of Korea. .,Department of Nano Engineering, Sungkyunkwan University (SKKU), Suwon 16419, Republic of Korea.,Center for Neuroscience Imaging Research, Institute for Basic Science (IBS), Suwon 16419, Republic of Korea
| | - Kilho Eom
- Biomechanics Laboratory, College of Sport Science, Sungkyunkwan University (SKKU), Suwon 16419, Republic of Korea.
| |
Collapse
|
16
|
Gong X, Zhang Y, Chen J. Advanced Sampling Methods for Multiscale Simulation of Disordered Proteins and Dynamic Interactions. Biomolecules 2021; 11:1416. [PMID: 34680048 PMCID: PMC8533332 DOI: 10.3390/biom11101416] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Revised: 09/22/2021] [Accepted: 09/24/2021] [Indexed: 11/16/2022] Open
Abstract
Intrinsically disordered proteins (IDPs) are highly prevalent and play important roles in biology and human diseases. It is now also recognized that many IDPs remain dynamic even in specific complexes and functional assemblies. Computer simulations are essential for deriving a molecular description of the disordered protein ensembles and dynamic interactions for a mechanistic understanding of IDPs in biology, diseases, and therapeutics. Here, we provide an in-depth review of recent advances in the multi-scale simulation of disordered protein states, with a particular emphasis on the development and application of advanced sampling techniques for studying IDPs. These techniques are critical for adequate sampling of the manifold functionally relevant conformational spaces of IDPs. Together with dramatically improved protein force fields, these advanced simulation approaches have achieved substantial success and demonstrated significant promise towards the quantitative and predictive modeling of IDPs and their dynamic interactions. We will also discuss important challenges remaining in the atomistic simulation of larger systems and how various coarse-grained approaches may help to bridge the remaining gaps in the accessible time- and length-scales of IDP simulations.
Collapse
Affiliation(s)
- Xiping Gong
- Department of Chemistry, University of Massachusetts Amherst, Amherst, MA 01003, USA; (X.G.); (Y.Z.)
| | - Yumeng Zhang
- Department of Chemistry, University of Massachusetts Amherst, Amherst, MA 01003, USA; (X.G.); (Y.Z.)
| | - Jianhan Chen
- Department of Chemistry, University of Massachusetts Amherst, Amherst, MA 01003, USA; (X.G.); (Y.Z.)
- Department of Biochemistry and Molecular Biology, University of Massachusetts Amherst, Amherst, MA 01003, USA
| |
Collapse
|
17
|
Alston JJ, Soranno A, Holehouse AS. Integrating single-molecule spectroscopy and simulations for the study of intrinsically disordered proteins. Methods 2021; 193:116-135. [PMID: 33831596 PMCID: PMC8713295 DOI: 10.1016/j.ymeth.2021.03.018] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2020] [Revised: 03/25/2021] [Accepted: 03/31/2021] [Indexed: 12/21/2022] Open
Abstract
Over the last two decades, intrinsically disordered proteins and protein regions (IDRs) have emerged from a niche corner of biophysics to be recognized as essential drivers of cellular function. Various techniques have provided fundamental insight into the function and dysfunction of IDRs. Among these techniques, single-molecule fluorescence spectroscopy and molecular simulations have played a major role in shaping our modern understanding of the sequence-encoded conformational behavior of disordered proteins. While both techniques are frequently used in isolation, when combined they offer synergistic and complementary information that can help uncover complex molecular details. Here we offer an overview of single-molecule fluorescence spectroscopy and molecular simulations in the context of studying disordered proteins. We discuss the various means in which simulations and single-molecule spectroscopy can be integrated, and consider a number of studies in which this integration has uncovered biological and biophysical mechanisms.
Collapse
Affiliation(s)
- Jhullian J Alston
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis 63110, MO, USA; Center for Science and Engineering of Living Systems (CSELS), Washington University in St. Louis, St. Louis 63130, MO, USA
| | - Andrea Soranno
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis 63110, MO, USA; Center for Science and Engineering of Living Systems (CSELS), Washington University in St. Louis, St. Louis 63130, MO, USA.
| | - Alex S Holehouse
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis 63110, MO, USA; Center for Science and Engineering of Living Systems (CSELS), Washington University in St. Louis, St. Louis 63130, MO, USA.
| |
Collapse
|
18
|
Song J, Li J, Chan HS. Small-Angle X-ray Scattering Signatures of Conformational Heterogeneity and Homogeneity of Disordered Protein Ensembles. J Phys Chem B 2021; 125:6451-6478. [PMID: 34115515 DOI: 10.1021/acs.jpcb.1c02453] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
An accurate account of disordered protein conformations is of central importance to deciphering the physicochemical basis of biological functions of intrinsically disordered proteins and the folding-unfolding energetics of globular proteins. Physically, disordered ensembles of nonhomopolymeric polypeptides are expected to be heterogeneous, i.e., they should differ from those homogeneous ensembles of homopolymers that harbor an essentially unique relationship between average values of end-to-end distance REE and radius of gyration Rg. It was posited recently, however, that small-angle X-ray scattering (SAXS) data on conformational dimensions of disordered proteins can be rationalized almost exclusively by homopolymer ensembles. Assessing this perspective, chain-model simulations are used to evaluate the discriminatory power of SAXS-determined molecular form factors (MFFs) with regard to homogeneous versus heterogeneous ensembles. The general approach adopted here is not bound by any assumption about ensemble encodability, in that the postulated heterogeneous ensembles we evaluated are not restricted to those entailed by simple interaction schemes. Our analysis of MFFs for certain heterogeneous ensembles with more narrowly distributed REE and Rg indicates that while they deviate from MFFs of homogeneous ensembles, the differences can be rather small. Remarkably, some heterogeneous ensembles with asphericity and REE drastically different from those of homogeneous ensembles can nonetheless exhibit practically identical MFFs, demonstrating that SAXS MFFs do not afford unique characterizations of basic properties of conformational ensembles in general. In other words, the ensemble to MFF mapping is practically many-to-one and likely nonsmooth. Heteropolymeric variations of the REE-Rg relationship were further showcased using an analytical perturbation theory developed here for flexible heteropolymers. Ramifications of our findings for interpretation of experimental data are discussed.
Collapse
Affiliation(s)
- Jianhui Song
- School of Polymer Science and Engineering, Qingdao University of Science and Technology, 53 Zhengzhou Road, Qingdao 266042, China
| | - Jichen Li
- School of Polymer Science and Engineering, Qingdao University of Science and Technology, 53 Zhengzhou Road, Qingdao 266042, China
| | - Hue Sun Chan
- Department of Biochemistry, University of Toronto Faculty of Medicine, Toronto, Ontario M5S 1A8, Canada
| |
Collapse
|
19
|
Computational methods for exploring protein conformations. Biochem Soc Trans 2021; 48:1707-1724. [PMID: 32756904 PMCID: PMC7458412 DOI: 10.1042/bst20200193] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2020] [Revised: 07/07/2020] [Accepted: 07/09/2020] [Indexed: 12/13/2022]
Abstract
Proteins are dynamic molecules that can transition between a potentially wide range of structures comprising their conformational ensemble. The nature of these conformations and their relative probabilities are described by a high-dimensional free energy landscape. While computer simulation techniques such as molecular dynamics simulations allow characterisation of the metastable conformational states and the transitions between them, and thus free energy landscapes, to be characterised, the barriers between states can be high, precluding efficient sampling without substantial computational resources. Over the past decades, a dizzying array of methods have emerged for enhancing conformational sampling, and for projecting the free energy landscape onto a reduced set of dimensions that allow conformational states to be distinguished, known as collective variables (CVs), along which sampling may be directed. Here, a brief description of what biomolecular simulation entails is followed by a more detailed exposition of the nature of CVs and methods for determining these, and, lastly, an overview of the myriad different approaches for enhancing conformational sampling, most of which rely upon CVs, including new advances in both CV determination and conformational sampling due to machine learning.
Collapse
|
20
|
Regy RM, Thompson J, Kim YC, Mittal J. Improved coarse-grained model for studying sequence dependent phase separation of disordered proteins. Protein Sci 2021; 30:1371-1379. [PMID: 33934416 DOI: 10.1002/pro.4094] [Citation(s) in RCA: 65] [Impact Index Per Article: 21.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2021] [Revised: 04/27/2021] [Accepted: 04/30/2021] [Indexed: 12/28/2022]
Abstract
We present improvements to the hydropathy scale (HPS) coarse-grained (CG) model for simulating sequence-specific behavior of intrinsically disordered proteins (IDPs), including their liquid-liquid phase separation (LLPS). The previous model based on an atomistic hydropathy scale by Kapcha and Rossky (KR scale) is not able to capture some well-known LLPS trends such as reduced phase separation propensity upon mutations (R-to-K and Y-to-F). Here, we propose to use the Urry hydropathy scale instead, which was derived from the inverse temperature transitions in a model polypeptide with guest residues X. We introduce two free parameters to shift (Δ) and scale (µ) the overall interaction strengths for the new model (HPS-Urry) and use the experimental radius of gyration for a diverse group of IDPs to find their optimal values. Interestingly, many possible (Δ, µ) combinations can be used for typical IDPs, but the phase behavior of a low-complexity (LC) sequence FUS is only well described by one of these models, which highlights the need for a careful validation strategy based on multiple proteins. The CG HPS-Urry model should enable accurate simulations of protein LLPS and provide a microscopically detailed view of molecular interactions.
Collapse
Affiliation(s)
- Roshan Mammen Regy
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania, USA
| | - Jacob Thompson
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania, USA
| | - Young C Kim
- Center for Materials Physics and Technology, Naval Research Laboratory, Washington, District of Columbia, USA
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania, USA
| |
Collapse
|
21
|
Shea JE, Best RB, Mittal J. Physics-based computational and theoretical approaches to intrinsically disordered proteins. Curr Opin Struct Biol 2021; 67:219-225. [PMID: 33545530 PMCID: PMC8150118 DOI: 10.1016/j.sbi.2020.12.012] [Citation(s) in RCA: 72] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2020] [Revised: 12/24/2020] [Accepted: 12/28/2020] [Indexed: 02/06/2023]
Abstract
Intrinsically disordered proteins (IDPs) are an important class of proteins that do not fold to a well-defined three-dimensional shape but rather adopt an ensemble of inter-converting conformations. This feature makes their experimental characterization challenging and invites a theoretical and computational approach to complement experimental studies. In this review, we highlight the recent progress in developing new computational and theoretical approaches to study the structure and dynamics of monomeric and order higher assemblies of IDPs, with a particular emphasis on their phase separation into protein-rich condensates.
Collapse
Affiliation(s)
- Joan-Emma Shea
- Department of Chemistry & Biochemistry, University of California, Santa Barbara, CA 93106, United States; Department of Physics, University of California, Santa Barbara, CA 93106, United States.
| | - Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, United States
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, 111 Research Drive, Bethlehem, PA 18015, United States.
| |
Collapse
|
22
|
Mu J, Liu H, Zhang J, Luo R, Chen HF. Recent Force Field Strategies for Intrinsically Disordered Proteins. J Chem Inf Model 2021; 61:1037-1047. [PMID: 33591749 PMCID: PMC8256680 DOI: 10.1021/acs.jcim.0c01175] [Citation(s) in RCA: 42] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Intrinsically disordered proteins (IDPs) are widely distributed across eukaryotic cells, playing important roles in molecular recognition, molecular assembly, post-translational modification, and other biological processes. IDPs are also associated with many diseases such as cancers, cardiovascular diseases, and neurodegenerative diseases. Due to their structural flexibility, conventional experimental methods cannot reliably capture their heterogeneous structures. Molecular dynamics simulation becomes an important complementary tool to quantify IDP structures. This review covers recent force field strategies proposed for more accurate molecular dynamics simulations of IDPs. The strategies include adjusting dihedral parameters, adding grid-based energy correction map (CMAP) parameters, refining protein-water interactions, and others. Different force fields were found to perform well on specific observables of specific IDPs but also are limited in reproducing all available experimental observables consistently for all tested IDPs. We conclude the review with perspective areas for improvements for future force fields for IDPs.
Collapse
Affiliation(s)
- Junxi Mu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Hao Liu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Jian Zhang
- Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, School of Medicine, Shanghai Jiao Tong University, Shanghai 20025, China
| | - Ray Luo
- Departments of Molecular Biology and Biochemistry, Chemical and Molecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine, California 92697-3900, United States
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| |
Collapse
|
23
|
Liu X, Gong X, Chen J. Accelerating atomistic simulations of proteins using multiscale enhanced sampling with independent tempering. J Comput Chem 2021; 42:358-364. [PMID: 33301208 DOI: 10.1002/jcc.26461] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2020] [Revised: 10/07/2020] [Accepted: 11/22/2020] [Indexed: 02/06/2023]
Abstract
Efficient sampling of the conformational space is essential for quantitative simulations of proteins. The multiscale enhanced sampling (MSES) method accelerates atomistic sampling by coupling it to a coarse-grained (CG) simulation. Bias from coupling to the CG model is removed using Hamiltonian replica exchange, such that one could benefit simultaneously from the high accuracy of atomistic models and fast dynamics of CG ones. Here, we extend MSES to allow independent control of the effective temperatures of atomistic and CG simulations, by directly scaling the atomistic and CG Hamiltonians. The new algorithm, named MSES with independent tempering (MSES-IT), supports more sophisticated Hamiltonian and temperature replica exchange protocols to further improve the sampling efficiency. Using a small but nontrivial β-hairpin, we show that setting the effective temperature of CG model in all conditions to its melting temperature maximizes structural transition rates at the CG level and promotes more efficient replica exchange and diffusion in the condition space. As the result, MSES-IT drive faster reversible transitions at the atomic level and leads to significant improvement in generating converged conformational ensembles compared to the original MSES scheme.
Collapse
Affiliation(s)
- Xiaorong Liu
- Department of Chemistry, University of Massachusetts Amherst, Amherst, Massachusetts, USA
| | - Xiping Gong
- Department of Chemistry, University of Massachusetts Amherst, Amherst, Massachusetts, USA
| | - Jianhan Chen
- Department of Chemistry, University of Massachusetts Amherst, Amherst, Massachusetts, USA.,Department of Biochemistry and Molecular Biology, University of Massachusetts Amherst, Amherst, Massachusetts, USA
| |
Collapse
|
24
|
Wang W. Recent advances in atomic molecular dynamics simulation of intrinsically disordered proteins. Phys Chem Chem Phys 2021; 23:777-784. [PMID: 33355572 DOI: 10.1039/d0cp05818a] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Intrinsically disordered proteins (IDPs) play important roles in cellular functions. The inherent structural heterogeneity of IDPs makes the high-resolution experimental characterization of IDPs extremely difficult. Molecular dynamics (MD) simulation could provide the atomic-level description of the structural and dynamic properties of IDPs. This perspective reviews the recent progress in atomic MD simulation studies of IDPs, including the development of force fields and sampling methods, as well as applications in IDP-involved protein-protein interactions. The employment of large-scale simulations and advanced sampling techniques allows more accurate estimation of the thermodynamics and kinetics of IDP-mediated protein interactions, and the holistic landscape of the binding process of IDPs is emerging.
Collapse
Affiliation(s)
- Wenning Wang
- Department of Chemistry, Multiscale Research Institute of Complex Systems and Institute of Biomedical Sciences, Fudan University, Shanghai 200438, China.
| |
Collapse
|
25
|
Konovalov A, Symons BCB, Popelier PLA. On the many-body nature of intramolecular forces in FFLUX and its implications. J Comput Chem 2020; 42:107-116. [PMID: 33107993 DOI: 10.1002/jcc.26438] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2020] [Revised: 09/28/2020] [Accepted: 09/28/2020] [Indexed: 12/24/2022]
Abstract
FFLUX is a biomolecular force field under construction, based on Quantum Chemical Topology (QCT) and machine learning (kriging), with a minimalistic and physically motivated design. A detailed analysis of the forces within the kriging models as treated in FFLUX is presented, taking as a test example a liquid water model. The energies of topological atoms are modeled as 3Natoms -6 dimensional potential energy surfaces, using atomic local frames to represent the internal degrees of freedom. As a result, the forces within the kriging models in FFLUX are inherently N-body in nature where N refers to Natoms . This provides a fuller picture that is closer to a true quantum mechanical representation of interactions between atoms. The presented computational example quantitatively showcases the non-negligible (as much as 9%) three-body nature of bonded forces and angular forces in a water molecule. We discuss the practical impact on the pressure calculation with N-body forces and periodic boundary conditions (PBC) in molecular dynamics, as opposed to classical force fields with two-body forces. The equivalence between the PBC-related correction terms in the general virial equation is shown mathematically.
Collapse
Affiliation(s)
- Anton Konovalov
- Manchester Institute of Biotechnology (MIB), Manchester, United Kingdom.,Department of Chemistry, University of Manchester, Manchester, United Kingdom
| | - Benjamin C B Symons
- Manchester Institute of Biotechnology (MIB), Manchester, United Kingdom.,Department of Chemistry, University of Manchester, Manchester, United Kingdom
| | - Paul L A Popelier
- Manchester Institute of Biotechnology (MIB), Manchester, United Kingdom.,Department of Chemistry, University of Manchester, Manchester, United Kingdom
| |
Collapse
|
26
|
Tang WS, Fawzi NL, Mittal J. Refining All-Atom Protein Force Fields for Polar-Rich, Prion-like, Low-Complexity Intrinsically Disordered Proteins. J Phys Chem B 2020; 124:9505-9512. [PMID: 33078950 DOI: 10.1021/acs.jpcb.0c07545] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Significant efforts in the past decade have given us highly accurate all-atom protein force fields for molecular dynamics (MD) simulations of folded and disordered proteins. These simulations, complemented with experimental data, provide new insights into molecular interactions that underlie the physical properties of proteins, especially for intrinsically disordered proteins (IDPs) for which defining the heterogeneous structural ensemble is hugely challenging by experiments alone. Consequently, the accuracy of these protein force fields is of utmost importance to ensure reliable simulated conformational data. Here, we first assess the accuracy of current state-of-the-art force fields for IDPs (ff99SBws and ff03ws) applied to disordered proteins of low amino acid sequence complexity that can undergo liquid-liquid phase separation. On the basis of a detailed comparison of NMR chemical shifts between simulation and experiment on several IDPs, we find that regions surrounding specific polar residues result in simulated ensembles with exaggerated helicity when compared to experiment. To resolve this discrepancy, we introduce residue-specific modifications to the backbone torsion potential of three residues (Ser, Thr, and Gln) in the ff99SBws force field. The modified force field, ff99SBws-STQ, provides a more accurate representation of helical structure propensity in these LC domains without compromising faithful representation of helicity in a region with distinct sequence composition. Our refinement strategy also suggests a path forward for integrating experimental data in the assessment of residue-specific deficiencies in the current physics-based force fields and improves these force fields further for their broader applicability.
Collapse
Affiliation(s)
- Wai Shing Tang
- Department of Physics, Brown University, Providence, Rhode Island 02912, United States
| | - Nicolas L Fawzi
- Department of Molecular Pharmacology, Physiology and Biotechnology, Brown University, Providence, Rhode Island 02912, United States
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|
27
|
Lou H, Cukier RI. A maximum entropy principle approach to a joint probability model for sequences with known neighbor and next neighbor pair probabilities. Chem Phys 2020. [DOI: 10.1016/j.chemphys.2020.110872] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
28
|
Optimization of Molecular Dynamics Simulations of c-MYC 1-88-An Intrinsically Disordered System. Life (Basel) 2020; 10:life10070109. [PMID: 32664335 PMCID: PMC7400636 DOI: 10.3390/life10070109] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2020] [Revised: 07/06/2020] [Accepted: 07/07/2020] [Indexed: 11/28/2022] Open
Abstract
Many of the proteins involved in key cellular regulatory events contain extensive intrinsically disordered regions that are not readily amenable to conventional structure/function dissection. The oncoprotein c-MYC plays a key role in controlling cell proliferation and apoptosis and more than 70% of the primary sequence is disordered. Computational approaches that shed light on the range of secondary and tertiary structural conformations therefore provide the only realistic chance to study such proteins. Here, we describe the results of several tests of force fields and water models employed in molecular dynamics simulations for the N-terminal 88 amino acids of c-MYC. Comparisons of the simulation data with experimental secondary structure assignments obtained by NMR establish a particular implicit solvation approach as highly congruent. The results provide insights into the structural dynamics of c-MYC1-88, which will be useful for guiding future experimental approaches. The protocols for trajectory analysis described here will be applicable for the analysis of a variety of computational simulations of intrinsically disordered proteins.
Collapse
|
29
|
Chen J, Liu X, Chen J. Targeting Intrinsically Disordered Proteins through Dynamic Interactions. Biomolecules 2020; 10:E743. [PMID: 32403216 PMCID: PMC7277182 DOI: 10.3390/biom10050743] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2020] [Revised: 05/04/2020] [Accepted: 05/09/2020] [Indexed: 12/18/2022] Open
Abstract
Intrinsically disordered proteins (IDPs) are over-represented in major disease pathways and have attracted significant interest in understanding if and how they may be targeted using small molecules for therapeutic purposes. While most existing studies have focused on extending the traditional structure-centric drug design strategies and emphasized exploring pre-existing structure features of IDPs for specific binding, several examples have also emerged to suggest that small molecules could achieve specificity in binding IDPs and affect their function through dynamic and transient interactions. These dynamic interactions can modulate the disordered conformational ensemble and often lead to modest compaction to shield functionally important interaction sites. Much work remains to be done on further elucidation of the molecular basis of the dynamic small molecule-IDP interaction and determining how it can be exploited for targeting IDPs in practice. These efforts will rely critically on an integrated experimental and computational framework for disordered protein ensemble characterization. In particular, exciting advances have been made in recent years in enhanced sampling techniques, Graphic Processing Unit (GPU)-computing, and protein force field optimization, which have now allowed rigorous physics-based atomistic simulations to generate reliable structure ensembles for nontrivial IDPs of modest sizes. Such de novo atomistic simulations will play crucial roles in exploring the exciting opportunity of targeting IDPs through dynamic interactions.
Collapse
Affiliation(s)
- Jianlin Chen
- Department of Hematology, Taizhou Central Hospital (Taizhou University Hospital), Taizhou 318000, Zhejiang, China;
| | - Xiaorong Liu
- Department of Chemistry, University of Massachusetts, Amherst, MA 01003, USA;
| | - Jianhan Chen
- Department of Chemistry, University of Massachusetts, Amherst, MA 01003, USA;
- Department of Biochemistry and Molecular Biology, University of Massachusetts, Amherst, MA 01003, USA
| |
Collapse
|
30
|
Soranno A. Physical basis of the disorder-order transition. Arch Biochem Biophys 2020; 685:108305. [DOI: 10.1016/j.abb.2020.108305] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2019] [Revised: 02/10/2020] [Accepted: 02/14/2020] [Indexed: 12/29/2022]
|
31
|
Nguyen PH, Sterpone F, Derreumaux P. Aggregation of disease-related peptides. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:435-460. [PMID: 32145950 DOI: 10.1016/bs.pmbts.2019.12.002] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Protein misfolding and aggregation of amyloid proteins is the fundamental cause of more than 20 diseases. Molecular mechanisms of the self-assembly and the formation of the toxic aggregates are still elusive. Computer simulations have been intensively used to study the aggregation of amyloid peptides of various amino acid lengths related to neurodegenerative diseases. We review atomistic and coarse-grained simulations of short amyloid peptides aimed at determining their transient oligomeric structures and the early and late aggregation steps.
Collapse
Affiliation(s)
- Phuong H Nguyen
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | - Fabio Sterpone
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | - Philippe Derreumaux
- Laboratory of Theoretical Chemistry, Ton Duc Thang University, Ho Chi Minh City, Vietnam; Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam.
| |
Collapse
|
32
|
Abstract
The oligomerization of Aβ16-22 peptide, which is the hydrophobic core region of full-length Aβ1-42, causes Alzheimer's disease (AD). This progressive neurodegenerative disease affects over 44 million people worldwide. However, very few synthesized drug molecules are available to inhibit the aggregation of Aβ. Recently, experimental studies have shown that the biological ATP molecule prevents Aβ fibrillation at the millimolar scale; however, the significance of ATP molecules on Aβ fibrillation and the mechanism behind it remain elusive. We have carried out a total of 7.5 μs extensive all-atom molecular dynamics and 8.82 μs of umbrella sampling in explicit water using AMBER14SB, AMBER99SB-ILDN, and AMBER-FB15 force fields for Aβ16-22 peptide, to investigate the role of ATP on the disruption of Aβ16-22 prefibrils. From various analyses, such as secondary structure analysis, residue-wise contact map, SASA, and interaction energies, we have observed that, in the presence of ATP, the aggregation of Aβ16-22 peptide is very unfavorable. Moreover, the biological molecule ATP interacts with the Aβ16-22 peptide via hydrogen bonding, π-π stacking, and NH-π interactions which, ultimately, prevent the aggregation of Aβ16-22 peptide. Hence, we assume that the deficiency of ATP may cause Alzheimer's disease (AD).
Collapse
Affiliation(s)
- Saikat Pal
- Department of Chemistry , Indian Institute of Technology , Guwahati , Assam 781039 , India
| | - Sandip Paul
- Department of Chemistry , Indian Institute of Technology , Guwahati , Assam 781039 , India
| |
Collapse
|
33
|
Yang J, Gao M, Xiong J, Su Z, Huang Y. Features of molecular recognition of intrinsically disordered proteins via coupled folding and binding. Protein Sci 2019; 28:1952-1965. [PMID: 31441158 PMCID: PMC6798136 DOI: 10.1002/pro.3718] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2019] [Revised: 08/16/2019] [Accepted: 08/20/2019] [Indexed: 12/12/2022]
Abstract
The sequence-structure-function paradigm of proteins has been revolutionized by the discovery of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs). In contrast to traditional ordered proteins, IDPs/IDRs are unstructured under physiological conditions. The absence of well-defined three-dimensional structures in the free state of IDPs/IDRs is fundamental to their function. Folding upon binding is an important mode of molecular recognition for IDPs/IDRs. While great efforts have been devoted to investigating the complex structures and binding kinetics and affinities, our knowledge on the binding mechanisms of IDPs/IDRs remains very limited. Here, we review recent advances on the binding mechanisms of IDPs/IDRs. The structures and kinetic parameters of IDPs/IDRs can vary greatly, and the binding mechanisms can be highly dependent on the structural properties of IDPs/IDRs. IDPs/IDRs can employ various combinations of conformational selection and induced fit in a binding process, which can be templated by the target and/or encoded by the IDP/IDR. Further studies should provide deeper insights into the molecular recognition of IDPs/IDRs and enable the rational design of IDP/IDR binding mechanisms in the future.
Collapse
Affiliation(s)
- Jing Yang
- Department of Biological Engineering and Key Laboratory of Industrial Fermentation (Ministry of Education)Hubei University of TechnologyWuhanHubeiChina
- Institute of Biomedical and Pharmaceutical SciencesHubei University of TechnologyWuhanHubeiChina
| | - Meng Gao
- Department of Biological Engineering and Key Laboratory of Industrial Fermentation (Ministry of Education)Hubei University of TechnologyWuhanHubeiChina
- Institute of Biomedical and Pharmaceutical SciencesHubei University of TechnologyWuhanHubeiChina
| | - Junwen Xiong
- Department of Biological Engineering and Key Laboratory of Industrial Fermentation (Ministry of Education)Hubei University of TechnologyWuhanHubeiChina
- Institute of Biomedical and Pharmaceutical SciencesHubei University of TechnologyWuhanHubeiChina
| | - Zhengding Su
- Department of Biological Engineering and Key Laboratory of Industrial Fermentation (Ministry of Education)Hubei University of TechnologyWuhanHubeiChina
- Institute of Biomedical and Pharmaceutical SciencesHubei University of TechnologyWuhanHubeiChina
| | - Yongqi Huang
- Department of Biological Engineering and Key Laboratory of Industrial Fermentation (Ministry of Education)Hubei University of TechnologyWuhanHubeiChina
- Institute of Biomedical and Pharmaceutical SciencesHubei University of TechnologyWuhanHubeiChina
| |
Collapse
|
34
|
Lohia R, Salari R, Brannigan G. Sequence specificity despite intrinsic disorder: How a disease-associated Val/Met polymorphism rearranges tertiary interactions in a long disordered protein. PLoS Comput Biol 2019; 15:e1007390. [PMID: 31626641 PMCID: PMC6821141 DOI: 10.1371/journal.pcbi.1007390] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2019] [Revised: 10/30/2019] [Accepted: 09/10/2019] [Indexed: 11/24/2022] Open
Abstract
The role of electrostatic interactions and mutations that change charge states in intrinsically disordered proteins (IDPs) is well-established, but many disease-associated mutations in IDPs are charge-neutral. The Val66Met single nucleotide polymorphism (SNP) in precursor brain-derived neurotrophic factor (BDNF) is one of the earliest SNPs to be associated with neuropsychiatric disorders, and the underlying molecular mechanism is unknown. Here we report on over 250 μs of fully-atomistic, explicit solvent, temperature replica-exchange molecular dynamics (MD) simulations of the 91 residue BDNF prodomain, for both the V66 and M66 sequence. The simulations were able to correctly reproduce the location of both local and non-local secondary structure changes due to the Val66Met mutation, when compared with NMR spectroscopy. We find that the change in local structure is mediated via entropic and sequence specific effects. We developed a hierarchical sequence-based framework for analysis and conceptualization, which first identifies “blobs” of 4-15 residues representing local globular regions or linkers. We use this framework within a novel test for enrichment of higher-order (tertiary) structure in disordered proteins; the size and shape of each blob is extracted from MD simulation of the real protein (RP), and used to parameterize a self-avoiding heterogenous polymer (SAHP). The SAHP version of the BDNF prodomain suggested a protein segmented into three regions, with a central long, highly disordered polyampholyte linker separating two globular regions. This effective segmentation was also observed in full simulations of the RP, but the Val66Met substitution significantly increased interactions across the linker, as well as the number of participating residues. The Val66Met substitution replaces β-bridging between V66 and V94 (on either side of the linker) with specific side-chain interactions between M66 and M95. The protein backbone in the vicinity of M95 is then free to form β-bridges with residues 31-41 near the N-terminus, which condenses the protein. A significant role for Met/Met interactions is consistent with previously-observed non-local effects of the Val66Met SNP, as well as established interactions between the Met66 sequence and a Met-rich receptor that initiates neuronal growth cone retraction. Intrinsically disordered proteins are proteins that have no well-defined structure in at least one functional form. Mutations in one amino acid may still affect their function significantly, especially in subtle ways with cumulative adverse effects on health. Here we report on molecular dynamics simulations of a protein that is critical for neuronal health throughout adulthood (brain-derived neurotrophic factor). We investigate the effects of a mutation carried by 30% of human population, which has been widely studied for its association with aging-related and stress-related disorders, reduced volume of the hippocampus, and variations in episodic memory. We identify a molecular mechanism in which the mutation may change the global conformations of the protein and its ability to bind to receptors.
Collapse
Affiliation(s)
- Ruchi Lohia
- Center for Computational and Integrative Biology, Rutgers University, Camden, New Jersey, United States of America
| | - Reza Salari
- Center for Computational and Integrative Biology, Rutgers University, Camden, New Jersey, United States of America
| | - Grace Brannigan
- Center for Computational and Integrative Biology, Rutgers University, Camden, New Jersey, United States of America
- Department of Physics, Rutgers University, Camden, New Jersey, United States of America
- * E-mail:
| |
Collapse
|
35
|
Feng CJ, Dhayalan B, Tokmakoff A. Refinement of Peptide Conformational Ensembles by 2D IR Spectroscopy: Application to Ala‒Ala‒Ala. Biophys J 2019; 114:2820-2832. [PMID: 29925019 DOI: 10.1016/j.bpj.2018.05.003] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2018] [Revised: 04/27/2018] [Accepted: 05/02/2018] [Indexed: 10/28/2022] Open
Abstract
Characterizing ensembles of intrinsically disordered proteins is experimentally challenging because of the ill-conditioned nature of ensemble determination with limited data and the intrinsic fast dynamics of the conformational ensemble. Amide I two-dimensional infrared (2D IR) spectroscopy has picosecond time resolution to freeze structural ensembles as needed for probing disordered-protein ensembles and conformational dynamics. Also, developments in amide I computational spectroscopy now allow a quantitative and direct prediction of amide I spectra based on conformational distributions drawn from molecular dynamics simulations, providing a route to ensemble refinement against experimental spectra. We performed a Bayesian ensemble refinement method on Ala-Ala-Ala against isotope-edited Fourier-transform infrared spectroscopy and 2D IR spectroscopy and tested potential factors affecting the quality of ensemble refinements. We found that isotope-edited 2D IR spectroscopy provides a stringent constraint on Ala-Ala-Ala conformations and returns consistent conformational ensembles with the dominant ppII conformer across varying prior distributions from many molecular dynamics force fields and water models. The dominant factor influencing ensemble refinements is the systematic frequency uncertainty from spectroscopic maps. However, the uncertainty of conformer populations can be significantly reduced by incorporating 2D IR spectra in addition to traditional Fourier-transform infrared spectra. Bayesian ensemble refinement against isotope-edited 2D IR spectroscopy thus provides a route to probe equilibrium-complex protein ensembles and potentially nonequilibrium conformational dynamics.
Collapse
Affiliation(s)
- Chi-Jui Feng
- Department of Chemistry, James Franck Institute and Institute for Biophysical Dynamics, University of Chicago, Chicago, Illinois
| | - Balamurugan Dhayalan
- Department of Chemistry, James Franck Institute and Institute for Biophysical Dynamics, University of Chicago, Chicago, Illinois
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute and Institute for Biophysical Dynamics, University of Chicago, Chicago, Illinois.
| |
Collapse
|
36
|
Levine ZA, Teranishi K, Okada AK, Langen R, Shea JE. The Mitochondrial Peptide Humanin Targets but Does Not Denature Amyloid Oligomers in Type II Diabetes. J Am Chem Soc 2019; 141:14168-14179. [DOI: 10.1021/jacs.9b04995] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Zachary A. Levine
- Department of Pathology, Yale School of Medicine, New Haven, Connecticut 06520, United States
- Department of Molecular Biophysics & Biochemistry, Yale University, New Haven, Connecticut 06520, United States
| | | | | | | | | |
Collapse
|
37
|
Lee M, Yoon J, Shin S. Computational Study on Structure and Aggregation Pathway of Aβ42 Amyloid Protofibril. J Phys Chem B 2019; 123:7859-7868. [DOI: 10.1021/acs.jpcb.9b07195] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- MinJun Lee
- Department of Chemistry, Seoul National University, Seoul 08826, Korea
| | - Jeseong Yoon
- Department of Chemistry, Seoul National University, Seoul 08826, Korea
| | - Seokmin Shin
- Department of Chemistry, Seoul National University, Seoul 08826, Korea
| |
Collapse
|
38
|
Thirumalai D, Samanta HS, Maity H, Reddy G. Universal Nature of Collapsibility in the Context of Protein Folding and Evolution. Trends Biochem Sci 2019; 44:675-687. [DOI: 10.1016/j.tibs.2019.04.003] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2018] [Revised: 04/03/2019] [Accepted: 04/04/2019] [Indexed: 12/21/2022]
|
39
|
Kasahara K, Terazawa H, Takahashi T, Higo J. Studies on Molecular Dynamics of Intrinsically Disordered Proteins and Their Fuzzy Complexes: A Mini-Review. Comput Struct Biotechnol J 2019; 17:712-720. [PMID: 31303975 PMCID: PMC6603302 DOI: 10.1016/j.csbj.2019.06.009] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2019] [Revised: 05/29/2019] [Accepted: 06/11/2019] [Indexed: 11/19/2022] Open
Abstract
The molecular dynamics (MD) method is a promising approach toward elucidating the molecular mechanisms of intrinsically disordered regions (IDRs) of proteins and their fuzzy complexes. This mini-review introduces recent studies that apply MD simulations to investigate the molecular recognition of IDRs. Firstly, methodological issues by which MD simulations treat IDRs, such as developing force fields, treating periodic boundary conditions, and enhanced sampling approaches, are discussed. Then, several examples of the applications of MD to investigate molecular interactions of IDRs in terms of the two kinds of complex formations; coupled-folding and binding and fuzzy complex. MD simulations provide insight into the molecular mechanisms of these binding processes by sampling conformational ensembles of flexible IDRs. In particular, we focused on all-atom explicit-solvent MD simulations except for studies of higher-order assembly of IDRs. Recent advances in MD methods, and computational power make it possible to dissect the molecular details of realistic molecular systems involving the dynamic behavior of IDRs.
Collapse
Affiliation(s)
- Kota Kasahara
- College of Life Sciences, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
- Corresponding author.
| | - Hiroki Terazawa
- Graduate School of Life Sciences, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
| | - Takuya Takahashi
- College of Life Sciences, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
| | - Junichi Higo
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima-minamimachi, Chuo-ku, Kobe 650-0047, Japan
| |
Collapse
|
40
|
Zerze GH, Zheng W, Best RB, Mittal J. Evolution of All-Atom Protein Force Fields to Improve Local and Global Properties. J Phys Chem Lett 2019; 10:2227-2234. [PMID: 30990694 PMCID: PMC7507668 DOI: 10.1021/acs.jpclett.9b00850] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
Experimental studies on intrinsically disordered and unfolded proteins have shown that in isolation they typically have low populations of secondary structure and exhibit distance scalings suggesting that they are at near-theta-solvent conditions. Until recently, however, all-atom force fields failed to reproduce these fundamental properties of intrinsically disordered proteins (IDPs). Recent improvements by refining against ensemble-averaged experimental observables for polypeptides in aqueous solution have addressed deficiencies including secondary structure bias, global conformational properties, and thermodynamic parameters of biophysical reactions such as folding and collapse. To date, studies utilizing these improved all-atom force fields have mostly been limited to a small set of unfolded or disordered proteins. Here, we present data generated for a diverse library of unfolded or disordered proteins using three progressively improved generations of Amber03 force fields, and we explore how global and local properties are affected by each successive change in the force field. We find that the most recent force field refinements significantly improve the agreement of the global properties such as radii of gyration and end-to-end distances with experimental estimates. However, these global properties are largely independent of the local secondary structure propensity. This result stresses the need to validate force fields with reference to a combination of experimental data providing information about both local and global structure formation.
Collapse
Affiliation(s)
- Gül H Zerze
- Department of Chemical and Biomolecular Engineering , Lehigh University , Bethlehem , Pennsylvania 18015 , United States
| | - Wenwei Zheng
- College of Integrative Sciences and Arts , Arizona State University , Mesa , Arizona 85212 , United States
| | - Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering , Lehigh University , Bethlehem , Pennsylvania 18015 , United States
| |
Collapse
|
41
|
Baul U, Chakraborty D, Mugnai ML, Straub JE, Thirumalai D. Sequence Effects on Size, Shape, and Structural Heterogeneity in Intrinsically Disordered Proteins. J Phys Chem B 2019; 123:3462-3474. [PMID: 30913885 PMCID: PMC6920032 DOI: 10.1021/acs.jpcb.9b02575] [Citation(s) in RCA: 103] [Impact Index Per Article: 20.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
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.
Collapse
Affiliation(s)
- Upayan Baul
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Debayan Chakraborty
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Mauro L. Mugnai
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - John E. Straub
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - D. Thirumalai
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
42
|
Blinov N, Wishart DS, Kovalenko A. Solvent Composition Effects on the Structural Properties of the Aβ42 Monomer from the 3D-RISM-KH Molecular Theory of Solvation. J Phys Chem B 2019; 123:2491-2506. [PMID: 30811210 DOI: 10.1021/acs.jpcb.9b00480] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Structural characterization of amyloid (A)β peptides implicated in Alzheimer's disease is a challenging problem due to their intrinsically disordered nature and their high propensity for aggregation. Only limited information is currently available from experiments on conformational properties and aggregation pathways of the peptides in cellular environments. In silico modeling complements experimental information, providing atomistic insight into structure and dynamics of different Aβ species. All-atom explicit solvent molecular dynamics (MD) simulations with a properly selected force field can deliver reliable structural and dynamic information. In the case of intrinsically disordered Aβ peptides, enhanced sampling simulations beyond the nanosecond time scale are required to obtain statistically meaningful results even for simple solvent conditions. To overcome the challenges of conformational sampling in crowded cellular environments, alternative approaches have to be used, including postprocessing of MD data. In this study, we employ the statistical-mechanical, three-dimensional reference interaction site model with the Kovalenko-Hirata closure integral equation molecular theory of solvation to describe solvent composition effects on the conformational equilibrium in a structural ensemble of the Aβ42 (covering residues 1-42) monomer based on a statistical reweighting technique. The methodology enables a computationally efficient prediction on how different factors in the cellular environment, such as solvent composition, nonpolar solvation, and macromolecular crowding, affect the structural properties of the monomer. Similarities have been identified between changes in the structural ensemble caused by nonpolar solvation and crowded environments modeled by ionic solution with large negative ions. In particular, both solvent conditions reduce the random coil content and enhance the helical structure content of the monomer. In contrast to the previous studies, which reported increased α-helical content of peptides in crowded environments, this work attributes these structural features to the difference in solvent exposure of hydrophilic residues of the monomer for different secondary structure elements, rather than to (entropic) excluded volume effects.
Collapse
Affiliation(s)
- Nikolay Blinov
- Department of Mechanical Engineering , Edmonton , Alberta T6G 1H9 , Canada.,Nanotechnology Research Centre , Edmonton , Alberta T6G 2M9 , Canada
| | - David S Wishart
- Departments of Computing Science and Biological Sciences , University of Alberta , Edmonton , Alberta T6G 2E8 , Canada
| | - Andriy Kovalenko
- Department of Mechanical Engineering , Edmonton , Alberta T6G 1H9 , Canada.,Nanotechnology Research Centre , Edmonton , Alberta T6G 2M9 , Canada
| |
Collapse
|
43
|
Extreme Fuzziness: Direct Interactions between Two IDPs. Biomolecules 2019; 9:biom9030081. [PMID: 30813629 PMCID: PMC6468500 DOI: 10.3390/biom9030081] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2018] [Revised: 02/10/2019] [Accepted: 02/18/2019] [Indexed: 01/06/2023] Open
Abstract
Protein interactions involving intrinsically disordered proteins (IDPs) greatly extend the range of binding mechanisms available to proteins. In interactions employing coupled folding and binding, IDPs undergo disorder-to-order transitions to form a complex with a well-defined structure. In many other cases, IDPs retain structural plasticity in the final complexes, which have been defined as the fuzzy complexes. While a large number of fuzzy complexes have been characterized with variety of fuzzy patterns, many of the interactions are between an IDP and a structured protein. Thus, whether two IDPs can interact directly to form a fuzzy complex without disorder-to-order transition remains an open question. Recently, two studies of interactions between IDPs (4.1G-CTD/NuMA and H1/ProTα) have found a definite answer to this question. Detailed characterizations combined with nuclear magnetic resonance (NMR), single-molecule Förster resonance energy transfer (smFRET) and molecular dynamics (MD) simulation demonstrate that direct interactions between these two pairs of IDPs do form fuzzy complexes while retaining the conformational dynamics of the isolated proteins, which we name as the extremely fuzzy complexes. Extreme fuzziness completes the full spectrum of protein-protein interaction modes, suggesting that a more generalized model beyond existing binding mechanisms is required. Previous models of protein interaction could be applicable to some aspects of the extremely fuzzy interactions, but in more general sense, the distinction between native and nonnative contacts, which was used to understand protein folding and binding, becomes obscure. Exploring the phenomenon of extreme fuzziness may shed new light on molecular recognition and drug design.
Collapse
|
44
|
Falahati H, Haji-Akbari A. Thermodynamically driven assemblies and liquid-liquid phase separations in biology. SOFT MATTER 2019; 15:1135-1154. [PMID: 30672955 DOI: 10.1039/c8sm02285b] [Citation(s) in RCA: 52] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The sustenance of life depends on the high degree of organization that prevails through different levels of living organisms, from subcellular structures such as biomolecular complexes and organelles to tissues and organs. The physical origin of such organization is not fully understood, and even though it is clear that cells and organisms cannot maintain their integrity without consuming energy, there is growing evidence that individual assembly processes can be thermodynamically driven and occur spontaneously due to changes in thermodynamic variables such as intermolecular interactions and concentration. Understanding the phase separation in vivo requires a multidisciplinary approach, integrating the theory and physics of phase separation with experimental and computational techniques. This paper aims at providing a brief overview of the physics of phase separation and its biological implications, with a particular focus on the assembly of membraneless organelles. We discuss the underlying physical principles of phase separation from its thermodynamics to its kinetics. We also overview the wide range of methods utilized for experimental verification and characterization of phase separation of membraneless organelles, as well as the utility of molecular simulations rooted in thermodynamics and statistical physics in understanding the governing principles of thermodynamically driven biological self-assembly processes.
Collapse
Affiliation(s)
- Hanieh Falahati
- Department of Neuroscience, Yale School of Medicine, New Haven, CT 06510, USA.
| | | |
Collapse
|
45
|
Moradi S, Nowroozi A, Shahlaei M. Shedding light on the structural properties of lipid bilayers using molecular dynamics simulation: a review study. RSC Adv 2019; 9:4644-4658. [PMID: 35520151 PMCID: PMC9060685 DOI: 10.1039/c8ra08441f] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2018] [Revised: 03/12/2019] [Accepted: 01/16/2019] [Indexed: 12/20/2022] Open
Abstract
This review gives an overview about the some of the most important possible analyzes, technical challenges, and existing protocols that can be performed on the biological membrane by the molecular dynamics simulation.
Collapse
Affiliation(s)
- Sajad Moradi
- Nano Drug Delivery Research Center
- Kermanshah University of Medical Sciences
- Kermanshah
- Iran
| | - Amin Nowroozi
- Pharmaceutical Sciences Research Center
- Faculty of Pharmacy
- Kermanshah University of Medical Sciences
- Kermanshah
- Iran
| | - Mohsen Shahlaei
- Medical Biology Research Center
- Kermanshah University of Medical Sciences
- Kermanshah
- Iran
| |
Collapse
|
46
|
Bakels S, Meijer E, Greuell M, Porskamp SBA, Rouwhorst G, Mahé J, Gaigeot MP, Rijs AM. Interactions of aggregating peptides probed by IR-UV action spectroscopy. Faraday Discuss 2019; 217:322-341. [DOI: 10.1039/c8fd00208h] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The interplay between intramolecular and formed inter-sheet hydrogen bonds and the effect of dispersion interactions on the formation of peptide dimers is studied using IR-UV action spectroscopy.
Collapse
Affiliation(s)
- Sjors Bakels
- Radboud University
- Institute for Molecules and Materials
- FELIX Laboratory
- 6525 ED Nijmegen
- The Netherlands
| | - Eline M. Meijer
- Radboud University
- Institute for Molecules and Materials
- FELIX Laboratory
- 6525 ED Nijmegen
- The Netherlands
| | - Mart Greuell
- Radboud University
- Institute for Molecules and Materials
- FELIX Laboratory
- 6525 ED Nijmegen
- The Netherlands
| | - Sebastiaan B. A. Porskamp
- Radboud University
- Institute for Molecules and Materials
- FELIX Laboratory
- 6525 ED Nijmegen
- The Netherlands
| | - George Rouwhorst
- Radboud University
- Institute for Molecules and Materials
- FELIX Laboratory
- 6525 ED Nijmegen
- The Netherlands
| | - Jerôme Mahé
- LAMBE CNRS UMR8587
- Université d’Evry val d’Essonne
- 91025 Evry
- France
| | | | - Anouk M. Rijs
- Radboud University
- Institute for Molecules and Materials
- FELIX Laboratory
- 6525 ED Nijmegen
- The Netherlands
| |
Collapse
|
47
|
Das S, Amin AN, Lin YH, Chan HS. Coarse-grained residue-based models of disordered protein condensates: utility and limitations of simple charge pattern parameters. Phys Chem Chem Phys 2018; 20:28558-28574. [PMID: 30397688 DOI: 10.1039/c8cp05095c] [Citation(s) in RCA: 76] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Biomolecular condensates undergirded by phase separations of proteins and nucleic acids serve crucial biological functions. To gain physical insights into their genetic basis, we study how liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) depends on their sequence charge patterns using a continuum Langevin chain model wherein each amino acid residue is represented by a single bead. Charge patterns are characterized by the "blockiness" measure κ and the "sequence charge decoration" (SCD) parameter. Consistent with random phase approximation (RPA) theory and lattice simulations, LLPS propensity as characterized by critical temperature Tcr* increases with increasingly negative SCD for a set of sequences showing a positive correlation between κ and -SCD. Relative to RPA, the simulated sequence-dependent variation in Tcr* is often-though not always-smaller, whereas the simulated critical volume fractions are higher. However, for a set of sequences exhibiting an anti-correlation between κ and -SCD, the simulated Tcr*'s are quite insensitive to either parameter. Additionally, we find that blocky sequences that allow for strong electrostatic repulsion can lead to coexistence curves with upward concavity as stipulated by RPA, but the LLPS propensity of a strictly alternating charge sequence was likely overestimated by RPA and lattice models because interchain stabilization of this sequence requires spatial alignments that are difficult to achieve in real space. These results help delineate the utility and limitations of the charge pattern parameters and of RPA, pointing to further efforts necessary for rationalizing the newly observed subtleties.
Collapse
Affiliation(s)
- Suman Das
- Department of Biochemistry, University of Toronto, Medical Sciences Building - 5th Fl., 1 King's College Circle, Toronto, Ontario M5S 1A8, Canada.
| | | | | | | |
Collapse
|
48
|
Samanta HS, Chakraborty D, Thirumalai D. Charge fluctuation effects on the shape of flexible polyampholytes with applications to intrinsically disordered proteins. J Chem Phys 2018; 149:163323. [PMID: 30384718 DOI: 10.1063/1.5035428] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
Random polyampholytes (PAs) contain positively and negatively charged monomers that are distributed randomly along the polymer chain. The interaction between charges is assumed to be given by the Debye-Huckel potential. We show that the size of the PA is determined by an interplay between electrostatic interactions, giving rise to the polyelectrolyte effect due to net charge per monomer (σ) and an effective attractive PA interaction due to charge fluctuations, δσ. The interplay between these terms gives rise to non-monotonic dependence of the radius of gyration, R g , on the inverse Debye length, κ, when PA effects are important ( δ σ σ > 1 ). In the opposite limit, R g decreases monotonically with increasing κ. Simulations of PA chains, using a charged bead-spring model, further corroborate our theoretical predictions. The simulations unambiguously show that conformational heterogeneity manifests itself among sequences that have identical PA parameters. A clear implication is that the phases of PA sequences, and by inference intrinsically disordered proteins (IDPs), cannot be determined using only the bare PA parameters (σ and δσ). The theory is used to calculate the changes in R g on N, the number of residues for a set of IDPs. For a certain class of IDPs, with N between 24 and 441, the size grows as R g ∼ N 0.6, which agrees with data from small angle X-ray scattering experiments.
Collapse
Affiliation(s)
- Himadri S Samanta
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| | - Debayan Chakraborty
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| | - D Thirumalai
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| |
Collapse
|
49
|
Joseph JA, Wales DJ. Intrinsically Disordered Landscapes for Human CD4 Receptor Peptide. J Phys Chem B 2018; 122:11906-11921. [DOI: 10.1021/acs.jpcb.8b08371] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Jerelle A. Joseph
- Department of Chemistry, University of Cambridge, Lenfield Road, Cambridge CB2 1EW, United Kingdom
| | - David J. Wales
- Department of Chemistry, University of Cambridge, Lenfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
50
|
Conformational preferences and phase behavior of intrinsically disordered low complexity sequences: insights from multiscale simulations. Curr Opin Struct Biol 2018; 56:1-10. [PMID: 30439585 DOI: 10.1016/j.sbi.2018.10.003] [Citation(s) in RCA: 64] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2018] [Revised: 10/17/2018] [Accepted: 10/19/2018] [Indexed: 11/22/2022]
Abstract
While many proteins and protein regions utilize a complex repertoire of amino acids to achieve their biological function, a subset of protein sequences are enriched in a reduced set of amino acids. These so-called low complexity (LC) sequences, specifically intrinsically disordered variants of LC sequences, have been the focus of recent investigations owing to their roles in a range of biological functions, specifically phase separation. Computational studies of LC sequences have provided rich insights into their behavior both as individual proteins in dilute solutions and as the drivers and modulators of higher-order assemblies. Here, we review how simulations performed across distinct resolutions have provided different types of insights into the biological role of LC sequences.
Collapse
|