1
|
Borges-Araújo L, Pereira GP, Valério M, Souza PCT. Assessing the Martini 3 protein model: A review of its path and potential. BIOCHIMICA ET BIOPHYSICA ACTA. PROTEINS AND PROTEOMICS 2024; 1872:141014. [PMID: 38670324 DOI: 10.1016/j.bbapap.2024.141014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Revised: 03/13/2024] [Accepted: 04/09/2024] [Indexed: 04/28/2024]
Abstract
Coarse-grained (CG) protein models have become indispensable tools for studying many biological protein details, from conformational dynamics to the organization of protein macro-complexes, and even the interaction of proteins with other molecules. The Martini force field is one of the most widely used CG models for bio-molecular simulations, partly because of the enormous success of its protein model. With the recent release of a new and improved version of the Martini force field - Martini 3 - a new iteration of its protein model was also made available. The Martini 3 protein force field is an evolution of its Martini 2 counterpart, aimed at improving many of the shortcomings that had been previously identified. In this mini-review, we first provide a general overview of the model and then focus on the successful advances made in the short time since its release, many of which would not have been possible before. Furthermore, we discuss reported limitations, potential directions for model improvement and comment on what the likely future development and application avenues are.
Collapse
Affiliation(s)
- Luís Borges-Araújo
- Laboratoire de Biologie et Modélisation de la Cellule, CNRS, UMR 5239, Inserm, U1293, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France; Centre Blaise Pascal de Simulation et de Modélisation Numérique, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France
| | - Gilberto P Pereira
- Laboratoire de Biologie et Modélisation de la Cellule, CNRS, UMR 5239, Inserm, U1293, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France; Centre Blaise Pascal de Simulation et de Modélisation Numérique, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France
| | - Mariana Valério
- Laboratoire de Biologie et Modélisation de la Cellule, CNRS, UMR 5239, Inserm, U1293, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France; Centre Blaise Pascal de Simulation et de Modélisation Numérique, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France
| | - Paulo C T Souza
- Laboratoire de Biologie et Modélisation de la Cellule, CNRS, UMR 5239, Inserm, U1293, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France; Centre Blaise Pascal de Simulation et de Modélisation Numérique, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France.
| |
Collapse
|
2
|
Weiand E, Koenig PH, Rodriguez-Ropero F, Roiter Y, Angioletti-Uberti S, Dini D, Ewen JP. Boundary Lubrication Performance of Polyelectrolyte-Surfactant Complexes on Biomimetic Surfaces. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:7933-7946. [PMID: 38573738 PMCID: PMC11025133 DOI: 10.1021/acs.langmuir.3c03737] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/03/2023] [Revised: 03/20/2024] [Accepted: 03/20/2024] [Indexed: 04/05/2024]
Abstract
Aqueous mixtures of oppositely charged polyelectrolytes and surfactants are useful in many industrial applications, such as shampoos and hair conditioners. In this work, we investigate the friction between biomimetic hair surfaces in the presence of adsorbed complexes formed from cationic polyelectrolytes and anionic surfactants in an aqueous solution. We apply nonequilibrium molecular dynamics (NEMD) simulations using the coarse-grained MARTINI model. We first developed new MARTINI parameters for cationic guar gum (CGG), a functionalized, plant-derived polysaccharide. The complexation of CGG and the anionic surfactant sodium dodecyl sulfate (SDS) on virgin and chemically damaged biomimetic hair surfaces was studied using a sequential adsorption approach. We then carried out squeeze-out and sliding NEMD simulations to assess the boundary lubrication performance of the CGG-SDS complex compressed between two hair surfaces. At low pressure, we observe a synergistic friction behavior for the CGG-SDS complex, which gives lower shear stress than either pure CGG or SDS. Here, friction is dominated by viscous dissipation in an interfacial layer comprising SDS and water. At higher pressures, which are probably beyond those usually experienced during hair manipulation, SDS and water are squeezed out, and friction increases due to interdigitation. The outcomes of this work are expected to be beneficial to fine-tune and screen sustainable hair care formulations to provide low friction and therefore a smooth feel and reduced entanglement.
Collapse
Affiliation(s)
- Erik Weiand
- Department
of Mechanical Engineering, Imperial College
London, South Kensington Campus, London SW7 2AZ, U.K.
- Institute
of Molecular Science and Engineering, Imperial
College London, South
Kensington Campus, London SW7 2AZ, U.K.
- Thomas
Young Centre for the Theory and Simulation of Materials, Imperial College London, South Kensington Campus, London SW7 2AZ, U.K.
| | - Peter H. Koenig
- Corporate
Functions Analytical and Data & Modeling Sciences, Mason Business
Center, The Procter and Gamble Company, Mason, Ohio 45040, United States
| | - Francisco Rodriguez-Ropero
- Corporate
Functions Analytical and Data & Modeling Sciences, Mason Business
Center, The Procter and Gamble Company, Mason, Ohio 45040, United States
| | - Yuri Roiter
- Corporate
Functions Analytical and Data & Modeling Sciences, Mason Business
Center, The Procter and Gamble Company, Mason, Ohio 45040, United States
| | - Stefano Angioletti-Uberti
- Institute
of Molecular Science and Engineering, Imperial
College London, South
Kensington Campus, London SW7 2AZ, U.K.
- Thomas
Young Centre for the Theory and Simulation of Materials, Imperial College London, South Kensington Campus, London SW7 2AZ, U.K.
- Department
of Materials, Imperial College London, South Kensington Campus, London SW7 2AZ, U.K.
| | - Daniele Dini
- Department
of Mechanical Engineering, Imperial College
London, South Kensington Campus, London SW7 2AZ, U.K.
- Institute
of Molecular Science and Engineering, Imperial
College London, South
Kensington Campus, London SW7 2AZ, U.K.
- Thomas
Young Centre for the Theory and Simulation of Materials, Imperial College London, South Kensington Campus, London SW7 2AZ, U.K.
| | - James P. Ewen
- Department
of Mechanical Engineering, Imperial College
London, South Kensington Campus, London SW7 2AZ, U.K.
- Institute
of Molecular Science and Engineering, Imperial
College London, South
Kensington Campus, London SW7 2AZ, U.K.
- Thomas
Young Centre for the Theory and Simulation of Materials, Imperial College London, South Kensington Campus, London SW7 2AZ, U.K.
| |
Collapse
|
3
|
Latham AP, Zhu L, Sharon DA, Ye S, Willard AP, Zhang X, Zhang B. Microphase Separation Produces Interfacial Environment within Diblock Biomolecular Condensates. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.03.30.534967. [PMID: 37034777 PMCID: PMC10081284 DOI: 10.1101/2023.03.30.534967] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
The phase separation of intrinsically disordered proteins is emerging as an important mechanism for cellular organization. However, efforts to connect protein sequences to the physical properties of condensates, i.e., the molecular grammar, are hampered by a lack of effective approaches for probing high-resolution structural details. Using a combination of multiscale simulations and fluorescence lifetime imaging microscopy experiments, we systematically explored a series of systems consisting of diblock elastin-like polypeptides (ELP). The simulations succeeded in reproducing the variation of condensate stability upon amino acid substitution and revealed different microenvironments within a single condensate, which we verified with environmentally sensitive fluorophores. The interspersion of hydrophilic and hydrophobic residues and a lack of secondary structure formation result in an interfacial environment, which explains both the strong correlation between ELP condensate stability and interfacial hydrophobicity scales, as well as the prevalence of protein-water hydrogen bonds. Our study uncovers new mechanisms for condensate stability and organization that may be broadly applicable.
Collapse
Affiliation(s)
- Andrew P Latham
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Longchen Zhu
- Department of Chemistry, School of Science and Research Center for Industries of the Future, Westlake University, 600 Dunyu Road, Hangzhou 310030, Zhejiang Province, China
- Westlake Laboratory of Life Sciences and Biomedicine, 18 Shilongshan Road, Hangzhou 310024, Zhejiang Province, China
- These authors contributed equally to this work
| | - Dina A Sharon
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
- These authors contributed equally to this work
| | - Songtao Ye
- Department of Chemistry, School of Science and Research Center for Industries of the Future, Westlake University, 600 Dunyu Road, Hangzhou 310030, Zhejiang Province, China
- Westlake Laboratory of Life Sciences and Biomedicine, 18 Shilongshan Road, Hangzhou 310024, Zhejiang Province, China
| | - Adam P Willard
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Xin Zhang
- Department of Chemistry, School of Science and Research Center for Industries of the Future, Westlake University, 600 Dunyu Road, Hangzhou 310030, Zhejiang Province, China
- Westlake Laboratory of Life Sciences and Biomedicine, 18 Shilongshan Road, Hangzhou 310024, Zhejiang Province, China
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
| |
Collapse
|
4
|
Rajpersaud T, Tabandeh S, Leon L, Loverde SM. Molecular Dynamics Simulations of Polyelectrolyte Complexes. Biomacromolecules 2024; 25:1468-1480. [PMID: 38366971 DOI: 10.1021/acs.biomac.3c01032] [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: 02/19/2024]
Abstract
Polyelectrolyte complexes (PECs) are currently of great interest due to their applications toward developing new adaptive materials and their relevance in membraneless organelles. These complexes emerge during phase separation when oppositely charged polymers are mixed in aqueous media. Peptide-based PECs are particularly useful toward developing new drug delivery methods due to their inherent biocompatibility. The underlying peptide sequence can be tuned to optimize specific material properties of the complex, such as interfacial tension and viscosity. Given their applicability, it would be advantageous to understand the underlying sequence-dependent phase behavior of oppositely charged peptides. Here, we report microsecond molecular dynamic simulations to characterize the effect of hydrophobicity on the sequence-dependent peptide conformation for model polypeptide sequences that were previously reported by Tabandeh et al. These sequences are designed with alternating chirality of the peptide backbone. We present microsecond simulations of six oppositely charged peptide pairs, characterizing the sequence-dependent effect on peptide size, degree of hydrogen bonding, secondary structure, and conformation. This analysis recapitulates sensible trends in peptide conformation and degree of hydrogen bonding, consistent with experimentally reported results. Ramachandran plots reveal that backbone conformation at the single amino acid level is highly influenced by the neighboring sequence in the chain. These results give insight into how subtle changes in hydrophobic side chain size and chirality influence the strength of hydrogen bonding between the chains and, ultimately, the secondary structure. Furthermore, principal component analysis reveals that the minimum energy structures may be subtly modulated by the underlying sequence.
Collapse
Affiliation(s)
- Tania Rajpersaud
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
| | - Sara Tabandeh
- Department of Materials Science and Engineering, University of Central Florida, 4000 Central Florida Boulevard, Orlando, FL 32816, United States
| | - Lorraine Leon
- Department of Materials Science and Engineering, University of Central Florida, 4000 Central Florida Boulevard, Orlando, FL 32816, United States
| | - Sharon M Loverde
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
- Department of Chemistry, College of Staten Island, The City University of New York, 2800 Victory Boulevard, Staten Island, NY 10314, United States
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY 10016, United States
| |
Collapse
|
5
|
Mondal S, Cui Q. Sequence Sensitivity in Membrane Remodeling by Polyampholyte Condensates. J Phys Chem B 2024; 128:2087-2099. [PMID: 38407041 DOI: 10.1021/acs.jpcb.3c08149] [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: 02/27/2024]
Abstract
Intrinsically disordered peptides (IDPs) have been found to undergo liquid-liquid phase separation (LLPS) and produce complex coacervates that play numerous regulatory roles in the cell. Recent experimental studies have discovered that LLPS at or near the membrane surface helps in the biomolecular organization during signaling events and can significantly alter the membrane morphology. However, the molecular mechanism and microscopic details of such processes still remain unclear. Here we study the effect of polyampholyte and polyelectrolyte condensation on two different anionic membranes, as they represent a majority of naturally occurring IDPs. The polyampholytes are fifty-residue polymers, made of glutamate(E) and lysine(K) with different charge patterns. The polyelectrolytes are separate chains of E25 and K25. We first calibrate the MARTINI v3.0 force field and then perform long-time-scale coarse-grained molecular dynamics simulations. We find that condensates formed by all the polyampholytes get adsorbed on the membrane. However, the strong polyampholytes (i.e., blocky sequences) can remodel the membranes more prominently than the weaker ones (i.e., scrambled sequences). Condensates formed by the blocky sequences induce a significant negative curvature (∼0.1 nm-1) and local demixing of lipids, whereas those by the scrambled sequences tend to wet the membrane to a greater extent without generating significant curvature or demixing. We perform several microscopic analyses to characterize the nature of the interaction between membranes and these condensates. Our analyses of interaction energetics reveal that membrane remodeling and/or wetting are favored by enhanced interactions between polyampholytes with lipids and the counterions.
Collapse
Affiliation(s)
- Sayantan Mondal
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, Massachusetts 02215, United States
| |
Collapse
|
6
|
Zerze GH. Optimizing the Martini 3 Force Field Reveals the Effects of the Intricate Balance between Protein-Water Interaction Strength and Salt Concentration on Biomolecular Condensate Formation. J Chem Theory Comput 2024; 20:1646-1655. [PMID: 37043540 DOI: 10.1021/acs.jctc.2c01273] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/13/2023]
Abstract
Condensation/dissolution has become a widely acknowledged biological macromolecular assembly phenomenon in subcellular compartmentalization. The MARTINI force field offers a coarse-grained protein model with a resolution that preserves molecular details with an explicit (CG) solvent. Despite its relatively higher resolution, it can still achieve condensate formation in a reasonable computing time with explicit solvent and ionic species. Therefore, it is highly desirable to tune this force field to be able to reproduce the experimentally observed properties of the condensate formation. In this work, we studied the condensate formation of the low-sequence complexity domain of fused in sarcoma protein using a MARTINI 3 force field by systematically modifying (increasing) the protein-water interaction strength and varying the salt concentration. We found that the condensate formation is sensitive both to the protein-water interaction strength and the presence of salt. While the unmodified MARTINI force field yields a complete collapse of proteins into one dense phase (i.e., no dilute phase), we reported a range of modified protein-water interaction strength that is capable of capturing the experimentally found transfer free energy between dense and dilute phases. We also found that the condensates lose their spherical shape upon the addition of salt, especially when the protein-water interactions are weak. Interchain amino acid contact map analysis showed one explanation for this observation: the protein-protein contact fraction reduces as salt is added to systems (when the protein-water interactions are weak), consistent with electrostatic screening effects. This reduction might be responsible for the condensates becoming nonspherical upon the addition of salt by reducing the need for minimizing the interfacial area. However, as the protein-water interactions become stronger to the extent that makes the transfer free energy agree well with experimentally observed transfer free energy, we found an increase in the protein-protein contact fraction upon the addition of salt, consistent with the salting-out effects. Therefore, we concluded that there is an intricate balance between screening effects and salting-out effects upon the addition of salt and this balance is highly sensitive to the strength of protein-water interactions.
Collapse
Affiliation(s)
- Gül H Zerze
- William A. Brookshire Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas 77204, United States
| |
Collapse
|
7
|
Staňo R, van Lente J, Lindhoud S, Košovan P. Sequestration of Small Ions and Weak Acids and Bases by a Polyelectrolyte Complex Studied by Simulation and Experiment. Macromolecules 2024; 57:1383-1398. [PMID: 38370910 PMCID: PMC10867894 DOI: 10.1021/acs.macromol.3c01209] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Revised: 11/27/2023] [Accepted: 11/29/2023] [Indexed: 02/20/2024]
Abstract
Mixing of oppositely charged polyelectrolytes can result in phase separation into a polymer-poor supernatant and a polymer-rich polyelectrolyte complex (PEC). We present a new coarse-grained model for the Grand-reaction method that enables us to determine the composition of the coexisting phases in a broad range of pH and salt concentrations. We validate the model by comparing it to recent simulations and experimental studies, as well as our own experiments on poly(acrylic acid)/poly(allylamine hydrochloride) complexes. The simulations using our model predict that monovalent ions partition approximately equally between both phases, whereas divalent ones accumulate in the PEC phase. On a semiquantitative level, these results agree with our own experiments, as well as with other experiments and simulations in the literature. In the sequel, we use the model to study the partitioning of a weak diprotic acid at various pH values of the supernatant. Our results show that the ionization of the acid is enhanced in the PEC phase, resulting in its preferential accumulation in this phase, which monotonically increases with the pH. Currently, this effect is still waiting to be confirmed experimentally. We explore how the model parameters (particle size, charge density, permittivity, and solvent quality) affect the measured partition coefficients, showing that fine-tuning of these parameters can make the agreement with the experiments almost quantitative. Nevertheless, our results show that charge regulation in multivalent solutes can potentially be exploited in engineering the partitioning of charged molecules in PEC-based systems at various pH values.
Collapse
Affiliation(s)
- Roman Staňo
- Faculty
of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
- Vienna
Doctoral School in Physics, University of
Vienna, Boltzmanngasse
5, 1090 Vienna, Austria
| | - Jéré
J. van Lente
- Department
of Molecules & Materials, University
of Twente, Drienerlolaan 5, 7522 NB Enschede, The Netherlands
| | - Saskia Lindhoud
- Department
of Molecules & Materials, University
of Twente, Drienerlolaan 5, 7522 NB Enschede, The Netherlands
| | - Peter Košovan
- Department
of Physical and Macromolecular Chemistry, Faculty of Science, Charles University, Hlavova 8, 128 40 Prague 2, Czech Republic
| |
Collapse
|
8
|
Hemmasi E, Tohidian M, Makki H. Morphology and Transport Study of Acid-Base Blend Proton Exchange Membranes by Molecular Simulations: Case of Chitosan/Nafion. J Phys Chem B 2023; 127:10624-10635. [PMID: 38037344 PMCID: PMC10726362 DOI: 10.1021/acs.jpcb.3c05332] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2023] [Revised: 10/28/2023] [Accepted: 11/20/2023] [Indexed: 12/02/2023]
Abstract
Blending a basic polymer (e.g., chitosan) with Nafion can modify some membrane properties in direct methanol fuel cell applications, e.g., controlling methanol crossover, by regulating the morphology of hydrophilic channels. Unraveling the mechanisms by which the channel morphology is modified is essential to formulate design strategies for acid-base blend membrane development. Thus, we use molecular simulations to analyze the morphological features of a blend membrane (at 75/25 chitosan/Nafion wt %), i.e., (i) water/polymer phase organizations, (ii) number and size of water clusters, and (iii) quantitative morphological measures of hydrophilic channels, and compare them to the pure Nafion in a wide range of water contents. It is found that the affinity of water to different hydrophilic groups in the blend membrane can result in more distorted and dispersed hydrophilic phase and fewer bulk water-like features compared to pure Nafion. Also, the width of the hydrophilic network bottleneck, i.e., pore limiting diameter (PLD), is found to be almost five times smaller for the blend membrane compared to Nafion at their maximum water contents. Moreover, by changing the chitosan/Nafion weight ratio from 75/25 to 0/100, we show that as Nafion content increases, all channel morphological characteristics alter monotonically except PLD. This is mainly due to the strong acid-base interactions between Nafion and chitosan, which hinder the monotonic growth of PLD. Interestingly, water and methanol diffusion coefficients are strongly correlated with PLD, suggesting that PLD can be used as a single parameter for tailoring the blending ratio for achieving the desired diffusion properties of acid-base membranes.
Collapse
Affiliation(s)
- Ehsan Hemmasi
- Department
of Polymer and Color Engineering, Amirkabir
University of Technology, 424 Hafez Avenue, Tehran 59163-4311, Iran
| | - Mahdi Tohidian
- Department
of Polymer and Color Engineering, Amirkabir
University of Technology, 424 Hafez Avenue, Tehran 59163-4311, Iran
| | - Hesam Makki
- Department
of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L69 7ZD, U.K.
| |
Collapse
|
9
|
Lin Z, Beneyton T, Baret JC, Martin N. Coacervate Droplets for Synthetic Cells. SMALL METHODS 2023; 7:e2300496. [PMID: 37462244 DOI: 10.1002/smtd.202300496] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Revised: 06/15/2023] [Indexed: 12/24/2023]
Abstract
The design and construction of synthetic cells - human-made microcompartments that mimic features of living cells - have experienced a real boom in the past decade. While many efforts have been geared toward assembling membrane-bounded compartments, coacervate droplets produced by liquid-liquid phase separation have emerged as an alternative membrane-free compartmentalization paradigm. Here, the dual role of coacervate droplets in synthetic cell research is discussed: encapsulated within membrane-enclosed compartments, coacervates act as surrogates of membraneless organelles ubiquitously found in living cells; alternatively, they can be viewed as crowded cytosol-like chassis for constructing integrated synthetic cells. After introducing key concepts of coacervation and illustrating the chemical diversity of coacervate systems, their physicochemical properties and resulting bioinspired functions are emphasized. Moving from suspensions of free floating coacervates, the two nascent roles of these droplets in synthetic cell research are highlighted: organelle-like modules and cytosol-like templates. Building the discussion on recent studies from the literature, the potential of coacervate droplets to assemble integrated synthetic cells capable of multiple life-inspired functions is showcased. Future challenges that are still to be tackled in the field are finally discussed.
Collapse
Affiliation(s)
- Zi Lin
- Université de Bordeaux, CNRS, Centre de Recherche Paul Pascal, UMR5031, 115 avenue du Dr. Schweitzer, 33600, Pessac, France
| | - Thomas Beneyton
- Université de Bordeaux, CNRS, Centre de Recherche Paul Pascal, UMR5031, 115 avenue du Dr. Schweitzer, 33600, Pessac, France
| | - Jean-Christophe Baret
- Université de Bordeaux, CNRS, Centre de Recherche Paul Pascal, UMR5031, 115 avenue du Dr. Schweitzer, 33600, Pessac, France
| | - Nicolas Martin
- Université de Bordeaux, CNRS, Centre de Recherche Paul Pascal, UMR5031, 115 avenue du Dr. Schweitzer, 33600, Pessac, France
| |
Collapse
|
10
|
Ingólfsson HI, Rizuan A, Liu X, Mohanty P, Souza PCT, Marrink SJ, Bowers MT, Mittal J, Berry J. Multiscale simulations reveal TDP-43 molecular-level interactions driving condensation. Biophys J 2023; 122:4370-4381. [PMID: 37853696 PMCID: PMC10720261 DOI: 10.1016/j.bpj.2023.10.016] [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: 12/08/2022] [Revised: 05/27/2023] [Accepted: 10/16/2023] [Indexed: 10/20/2023] Open
Abstract
The RNA-binding protein TDP-43 is associated with mRNA processing and transport from the nucleus to the cytoplasm. TDP-43 localizes in the nucleus as well as accumulating in cytoplasmic condensates such as stress granules. Aggregation and formation of amyloid-like fibrils of cytoplasmic TDP-43 are hallmarks of numerous neurodegenerative diseases, most strikingly present in >90% of amyotrophic lateral sclerosis (ALS) patients. If excessive accumulation of cytoplasmic TDP-43 causes, or is caused by, neurodegeneration is presently not known. In this work, we use molecular dynamics simulations at multiple resolutions to explore TDP-43 self- and cross-interaction dynamics. A full-length molecular model of TDP-43, all 414 amino acids, was constructed from select structures of the protein functional domains (N-terminal domain, and two RNA recognition motifs, RRM1 and RRM2) and modeling of disordered connecting loops and the low complexity glycine-rich C-terminus domain. All-atom CHARMM36m simulations of single TDP-43 proteins served as guides to construct a coarse-grained Martini 3 model of TDP-43. The Martini model and a coarser implicit solvent C⍺ model, optimized for disordered proteins, were subsequently used to probe TDP-43 interactions; self-interactions from single-chain full-length TDP-43 simulations, cross-interactions from simulations with two proteins and simulations with assemblies of dozens to hundreds of proteins. Our findings illustrate the utility of different modeling scales for accessing TDP-43 molecular-level interactions and suggest that TDP-43 has numerous interaction preferences or patterns, exhibiting an overall strong, but dynamic, association and driving the formation of biomolecular condensates.
Collapse
Affiliation(s)
- Helgi I Ingólfsson
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California.
| | - Azamat Rizuan
- Artie McFerrin Department of Chemical Engineering, Texas A&M College of Engineering, College Station, Texas
| | - Xikun Liu
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California; Department of Chemistry & Biochemistry, University of California Santa Barbara, Santa Barbara, California
| | - Priyesh Mohanty
- Artie McFerrin Department of Chemical Engineering, Texas A&M College of Engineering, College Station, Texas
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS & University of Lyon, Lyon, France; Laboratory of Biology and Modeling of the Cell, École Normale Supérieure de Lyon, Université Claude Bernard Lyon 1, CNRS UMR 5239 and Inserm U1293, 46 Allée d'Italie, Lyon, France
| | - Siewert J Marrink
- Groningen Biomolecular Science and Biotechnology Institute, University of Groningen, Groningen, the Netherlands
| | - Michael T Bowers
- Department of Chemistry & Biochemistry, University of California Santa Barbara, Santa Barbara, California
| | - Jeetain Mittal
- Artie McFerrin Department of Chemical Engineering, Texas A&M College of Engineering, College Station, Texas; Department of Chemistry, Texas A&M University, College Station, Texas; Interdisciplinary Graduate Program in Genetics and Genomics, Texas A&M University, College Station, Texas
| | - Joel Berry
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California
| |
Collapse
|
11
|
Borges-Araújo L, Patmanidis I, Singh AP, Santos LHS, Sieradzan AK, Vanni S, Czaplewski C, Pantano S, Shinoda W, Monticelli L, Liwo A, Marrink SJ, Souza PCT. Pragmatic Coarse-Graining of Proteins: Models and Applications. J Chem Theory Comput 2023; 19:7112-7135. [PMID: 37788237 DOI: 10.1021/acs.jctc.3c00733] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.
Collapse
Affiliation(s)
- Luís Borges-Araújo
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Ilias Patmanidis
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Akhil P Singh
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
| | - Lucianna H S Santos
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
- Institut de Pharmacologie Moléculaire et Cellulaire, Université Côte d'Azur, Inserm, CNRS, 06560 Valbonne, France
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Sergio Pantano
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita, Okayama 700-8530, Japan
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Adam Liwo
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| |
Collapse
|
12
|
Sinha NJ, Cunha KC, Murphy R, Hawker CJ, Shea JE, Helgeson ME. Competition between β-Sheet and Coacervate Domains Yields Diverse Morphologies in Mixtures of Oppositely Charged Homochiral Polypeptides. Biomacromolecules 2023; 24:3580-3588. [PMID: 37486022 DOI: 10.1021/acs.biomac.3c00361] [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: 07/25/2023]
Abstract
Biomolecular assembly processes involving competition between specific intermolecular interactions and thermodynamic phase instability have been implicated in a number of pathological states and technological applications of biomaterials. As a model for such processes, aqueous mixtures of oppositely charged homochiral polypeptides such as poly-l-lysine and poly-l-glutamic acid have been reported to form either β-sheet-rich solid-like precipitates or liquid-like coacervate droplets depending on competing hydrogen bonding interactions. Herein, we report studies of polypeptide mixtures that reveal unexpectedly diverse morphologies ranging from partially coalescing and aggregated droplets to bulk precipitates, as well as a previously unreported re-entrant liquid-liquid phase separation at high polypeptide concentration and ionic strength. Combining our experimental results with all-atom molecular dynamics simulations of folded polypeptide complexes reveals a concentration dependence of β-sheet-rich secondary structure, whose relative composition correlates with the observed macroscale morphologies of the mixtures. These results elucidate a crucial balance of interactions that are important for controlling morphology during coacervation in these and potentially similar biologically relevant systems.
Collapse
Affiliation(s)
- Nairiti J Sinha
- Department of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106, United States
- Materials Department and Materials Research Laboratory, University of California Santa Barbara, Santa Barbara, California 93106, United States
- Department of Chemistry and Biochemistry, University of California Santa Barbara, Santa Barbara, California 93106, United States
| | - Keila Cristina Cunha
- Department of Chemistry and Biochemistry, University of California Santa Barbara, Santa Barbara, California 93106, United States
| | - Robert Murphy
- Department of Chemistry and Biochemistry, University of California Santa Barbara, Santa Barbara, California 93106, United States
| | - Craig J Hawker
- Materials Department and Materials Research Laboratory, University of California Santa Barbara, Santa Barbara, California 93106, United States
- Department of Chemistry and Biochemistry, University of California Santa Barbara, Santa Barbara, California 93106, United States
| | - Joan-Emma Shea
- Department of Chemistry and Biochemistry, University of California Santa Barbara, Santa Barbara, California 93106, United States
| | - Matthew E Helgeson
- Department of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106, United States
| |
Collapse
|
13
|
Szała-Mendyk B, Phan TM, Mohanty P, Mittal J. Challenges in studying the liquid-to-solid phase transitions of proteins using computer simulations. Curr Opin Chem Biol 2023; 75:102333. [PMID: 37267850 PMCID: PMC10527940 DOI: 10.1016/j.cbpa.2023.102333] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2023] [Revised: 04/24/2023] [Accepted: 04/30/2023] [Indexed: 06/04/2023]
Abstract
"Membraneless organelles," also referred to as biomolecular condensates, perform a variety of cellular functions and their dysregulation is implicated in cancer and neurodegeneration. In the last two decades, liquid-liquid phase separation (LLPS) of intrinsically disordered and multidomain proteins has emerged as a plausible mechanism underlying the formation of various biomolecular condensates. Further, the occurrence of liquid-to-solid transitions within liquid-like condensates may give rise to amyloid structures, implying a biophysical link between phase separation and protein aggregation. Despite significant advances, uncovering the microscopic details of liquid-to-solid phase transitions using experiments remains a considerable challenge and presents an exciting opportunity for the development of computational models which provide valuable, complementary insights into the underlying phenomenon. In this review, we first highlight recent biophysical studies which provide new insights into the molecular mechanisms underlying liquid-to-solid (fibril) phase transitions of folded, disordered and multi-domain proteins. Next, we summarize the range of computational models used to study protein aggregation and phase separation. Finally, we discuss recent computational approaches which attempt to capture the underlying physics of liquid-to-solid transitions along with their merits and shortcomings.
Collapse
Affiliation(s)
- Beata Szała-Mendyk
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, TAMU 3127, College Station, 77843, Texas, United States.
| | - Tien Minh Phan
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, TAMU 3127, College Station, 77843, Texas, United States.
| | - Priyesh Mohanty
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, TAMU 3127, College Station, 77843, Texas, United States.
| | - Jeetain Mittal
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, TAMU 3127, College Station, 77843, Texas, United States; Department of Chemistry, Texas A&M University, TAMU 3255, College Station, 77843, Texas, United States; Interdisciplinary Graduate Program in Genetics and Genomics, Texas A&M University, TAMU 3255, College Station, 77843, Texas, United States.
| |
Collapse
|
14
|
Fagerberg E, Skepö M. Comparative Performance of Computer Simulation Models of Intrinsically Disordered Proteins at Different Levels of Coarse-Graining. J Chem Inf Model 2023. [PMID: 37339604 DOI: 10.1021/acs.jcim.3c00113] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/22/2023]
Abstract
Coarse-graining is commonly used to decrease the computational cost of simulations. However, coarse-grained models are also considered to have lower transferability, with lower accuracy for systems outside the original scope of parametrization. Here, we benchmark a bead-necklace model and a modified Martini 2 model, both coarse-grained models, for a set of intrinsically disordered proteins, with the different models having different degrees of coarse-graining. The SOP-IDP model has earlier been used for this set of proteins; thus, those results are included in this study to compare how models with different levels of coarse-graining compare. The sometimes naive expectation of the least coarse-grained model performing best does not hold true for the experimental pool of proteins used here. Instead, it showed the least good agreement, indicating that one should not necessarily trust the otherwise intuitive notion of a more advanced model inherently being better in model choice.
Collapse
Affiliation(s)
- Eric Fagerberg
- Theoretical Chemistry, Lund University, POB 124, SE-221 00 Lund, Sweden
| | - Marie Skepö
- Theoretical Chemistry, Lund University, POB 124, SE-221 00 Lund, Sweden
- LINXS - Institute of Advanced Neutron and X-ray Science, Scheelevägen 19, SE-223 70 Lund, Sweden
| |
Collapse
|
15
|
Hilpert C, Beranger L, Souza PCT, Vainikka PA, Nieto V, Marrink SJ, Monticelli L, Launay G. Facilitating CG Simulations with MAD: The MArtini Database Server. J Chem Inf Model 2023; 63:702-710. [PMID: 36656159 DOI: 10.1021/acs.jcim.2c01375] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The MArtini Database (MAD - https://mad.ibcp.fr) is a web server designed for the sharing of structures and topologies of molecules parametrized with the Martini coarse-grained (CG) force field. MAD can also convert atomistic structures into CG structures and prepare complex systems (including proteins, lipids, etc.) for molecular dynamics (MD) simulations at the CG level. It is dedicated to the generation of input files for Martini 3, the most recent version of this popular CG force field. Specifically, the MAD server currently includes tools to submit or retrieve CG models of a wide range of molecules (lipids, carbohydrates, nanoparticles, etc.), transform atomistic protein structures into CG structures and topologies, with fine control on the process and assemble biomolecules into large systems, and deliver all files necessary to start simulations in the GROMACS MD engine.
Collapse
Affiliation(s)
- Cécile Hilpert
- Microbiologie Moléculaire et Biochimie Structurale (MMSB), UMR 5086 CNRS & University of Lyon. 7 passage du Vercors, 69367 Lyon, France
| | - Louis Beranger
- Microbiologie Moléculaire et Biochimie Structurale (MMSB), UMR 5086 CNRS & University of Lyon. 7 passage du Vercors, 69367 Lyon, France
| | - Paulo C T Souza
- Microbiologie Moléculaire et Biochimie Structurale (MMSB), UMR 5086 CNRS & University of Lyon. 7 passage du Vercors, 69367 Lyon, France
| | - Petteri A Vainikka
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Vincent Nieto
- Microbiologie Moléculaire et Biochimie Structurale (MMSB), UMR 5086 CNRS & University of Lyon. 7 passage du Vercors, 69367 Lyon, France
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Luca Monticelli
- Microbiologie Moléculaire et Biochimie Structurale (MMSB), UMR 5086 CNRS & University of Lyon. 7 passage du Vercors, 69367 Lyon, France
| | - Guillaume Launay
- Microbiologie Moléculaire et Biochimie Structurale (MMSB), UMR 5086 CNRS & University of Lyon. 7 passage du Vercors, 69367 Lyon, France
| |
Collapse
|
16
|
Paloncýová M, Šrejber M, Čechová P, Kührová P, Zaoral F, Otyepka M. Atomistic Insights into Organization of RNA-Loaded Lipid Nanoparticles. J Phys Chem B 2023; 127:1158-1166. [PMID: 36602349 DOI: 10.1021/acs.jpcb.2c07671] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
RNA-based therapies have shown promise in a wide range of applications, from cancer therapy, treatment of inherited diseases to vaccination. Encapsulation of RNA into ionizable lipid (IL) containing lipid nanoparticles (LNPs) has enabled its safe and targeted delivery. We present here the simulations of the self-assembly process of pH-sensitive RNA-carrying LNPs and their internal morphology. At low pH, the simulations confirm a lipid core encapsulating RNA in the hexagonal phase. Our all-atom and coarse-grained simulations show that an RNA molecule inside an LNP is protected from interactions with ions by being enveloped in the charged ILs. At neutral pH, representing the environment after LNP administration into human tissues, LNPs expelled most of the encapsulated RNA and water and formed separate bulk IL-rich and ordered the helper-lipid-rich phase. Helper lipids arranged themselves to be in contact with RNA or water. The presented models provide atomistic understanding of the LNP structure and open a way to investigate them in silico, varying the LNP composition or interacting with other biostructures aiming at increasing the efficiency of RNA-based medicine.
Collapse
Affiliation(s)
- Markéta Paloncýová
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc 779 00, Czech Republic
| | - Martin Šrejber
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc 779 00, Czech Republic
| | - Petra Čechová
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc 779 00, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc 779 00, Czech Republic
| | - Filip Zaoral
- IT4Innovations, VŠB - Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba 708 00, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc 779 00, Czech Republic.,IT4Innovations, VŠB - Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba 708 00, Czech Republic
| |
Collapse
|
17
|
Liu Y, Wang X, Wan Z, Ngai T, Tse YLS. Capturing coacervate formation and protein partition by molecular dynamics simulation. Chem Sci 2023; 14:1168-1175. [PMID: 36756326 PMCID: PMC9891350 DOI: 10.1039/d2sc01164f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2022] [Accepted: 12/22/2022] [Indexed: 12/25/2022] Open
Abstract
Biomolecules localize and function in microenvironments where their local concentration, spatial organization, and biochemical reactivity are regulated. To compartmentalize and control the local properties of the native microenvironment, cellular mimics and artificial bioreactors have been developed to study the properties of membraneless organelles or mimic the bio-environment for life origin. Here, we carried out molecular dynamics simulation with the Martini 3.0 model to reproduce the experimental salt concentration and pH dependency of different complex coacervates. We showed that coacervates inside vesicles are able to change their shape. In addition, we used these coacervate systems to explore the partitioning of the ubiquitous cytoskeletal protein actin and found that actin spontaneously partitions to all the coacervate peripheries. Therefore, we believe that our study can provide a better understanding of the versatile coacervate platform, where biomolecules partition and gather to fulfill their biological duties.
Collapse
Affiliation(s)
- Yang Liu
- College of Polymer Science and Engineering, State Key Laboratory of Polymer Materials Engineering, Sichuan University Chengdu 610065 China.,Department of Chemistry, The Chinese University of Hong Kong Sha Tin Hong Kong China
| | - Xinyan Wang
- Department of Chemistry, The Chinese University of Hong Kong Sha Tin Hong Kong China
| | - Zhili Wan
- School of Food Science and Engineering, South China University of TechnologyGuangzhouChina
| | - To Ngai
- Department of Chemistry, The Chinese University of Hong Kong Sha Tin Hong Kong China
| | - Ying-Lung Steve Tse
- Department of Chemistry, The Chinese University of Hong Kong Sha Tin Hong Kong China
| |
Collapse
|
18
|
Mabuchi T, Kijima J, Yamashita Y, Miura E, Muraoka T. Coacervate Formation of Elastin-like Polypeptides in Explicit Aqueous Solution Using Coarse-Grained Molecular Dynamics Simulations. Macromolecules 2023. [DOI: 10.1021/acs.macromol.2c02195] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Takuya Mabuchi
- Frontier Research Institute for Interdisciplinary Sciences, Tohoku University, 2-1-1 Katahira,
Aoba-ku, Sendai, Miyagi980-8577, Japan
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira,
Aoba-ku, Sendai, Miyagi980-8577, Japan
| | - Junko Kijima
- Frontier Research Institute for Interdisciplinary Sciences, Tohoku University, 2-1-1 Katahira,
Aoba-ku, Sendai, Miyagi980-8577, Japan
| | - Yukino Yamashita
- Department of Applied Chemistry, Graduate School of Engineering, Tokyo University of Agriculture and Technology, 2-24-16, Naka-cho,
Koganei, Tokyo184-8588, Japan
| | - Erika Miura
- Department of Applied Chemistry, Graduate School of Engineering, Tokyo University of Agriculture and Technology, 2-24-16, Naka-cho,
Koganei, Tokyo184-8588, Japan
| | - Takahiro Muraoka
- Department of Applied Chemistry, Graduate School of Engineering, Tokyo University of Agriculture and Technology, 2-24-16, Naka-cho,
Koganei, Tokyo184-8588, Japan
| |
Collapse
|
19
|
Valdes-Garcia G, Heo L, Lapidus LJ, Feig M. Modeling Concentration-dependent Phase Separation Processes Involving Peptides and RNA via Residue-Based Coarse-Graining. J Chem Theory Comput 2023; 19:10.1021/acs.jctc.2c00856. [PMID: 36607820 PMCID: PMC10323037 DOI: 10.1021/acs.jctc.2c00856] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Biomolecular condensation, especially liquid-liquid phase separation, is an important physical process with relevance for a number of different aspects of biological functions. Key questions of what drives such condensation, especially in terms of molecular composition, can be addressed via computer simulations, but the development of computationally efficient yet physically realistic models has been challenging. Here, the coarse-grained model COCOMO is introduced that balances the polymer behavior of peptides and RNA chains with their propensity to phase separate as a function of composition and concentration. COCOMO is a residue-based model that combines bonded terms with short- and long-range terms, including a Debye-Hückel solvation term. The model is highly predictive of experimental data on phase-separating model systems. It is also computationally efficient and can reach the spatial and temporal scales on which biomolecular condensation is observed with moderate computational resources.
Collapse
Affiliation(s)
- Gilberto Valdes-Garcia
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824, USA
| | - Lim Heo
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824, USA
| | - Lisa J. Lapidus
- Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA
| | - Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824, USA
| |
Collapse
|
20
|
Grünewald F, Punt MH, Jefferys EE, Vainikka PA, König M, Virtanen V, Meyer TA, Pezeshkian W, Gormley AJ, Karonen M, Sansom MSP, Souza PCT, Marrink SJ. Martini 3 Coarse-Grained Force Field for Carbohydrates. J Chem Theory Comput 2022; 18:7555-7569. [PMID: 36342474 DOI: 10.1021/acs.jctc.2c00757] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The Martini 3 force field is a full reparametrization of the Martini coarse-grained model for biomolecular simulations. Due to the improved interaction balance, it allows for a more accurate description of condensed phase systems. In the present work, we develop a consistent strategy to parametrize carbohydrate molecules accurately within the framework of Martini 3. In particular, we develop a canonical mapping scheme which decomposes arbitrarily large carbohydrates into a limited number of fragments. Bead types for these fragments have been assigned by matching physicochemical properties of mono- and disaccharides. In addition, guidelines for assigning bonds, angles, and dihedrals were developed. These guidelines enable a more accurate description of carbohydrate conformations than in the Martini 2 force field. We show that models obtained with this approach are able to accurately reproduce osmotic pressures of carbohydrate water solutions. Furthermore, we provide evidence that the model differentiates correctly the solubility of the polyglucoses dextran (water-soluble) and cellulose (water insoluble but soluble in ionic liquids). Finally, we demonstrate that the new building blocks can be applied to glycolipids. We show they are able to reproduce membrane properties and induce binding of peripheral membrane proteins. These test cases demonstrate the validity and transferability of our approach.
Collapse
Affiliation(s)
- Fabian Grünewald
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Groningen 9747 AG, The Netherlands
| | - Mats H Punt
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Groningen 9747 AG, The Netherlands
| | - Elizabeth E Jefferys
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Petteri A Vainikka
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Groningen 9747 AG, The Netherlands
| | - Melanie König
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Groningen 9747 AG, The Netherlands
| | - Valtteri Virtanen
- Natural Chemistry Research Group, Department of Chemistry, University of Turku, FI-20014 Turku, Finland
| | - Travis A Meyer
- Department of Biomedical Engineering, Rutgers, The State University of New Jersey, Piscataway, New Jersey 08854, United States
| | - Weria Pezeshkian
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Groningen 9747 AG, The Netherlands.,The Niels Bohr International Academy, Niels Bohr Institute, University of Copenhagen, Copenhagen 2100, Denmark
| | - Adam J Gormley
- Department of Biomedical Engineering, Rutgers, The State University of New Jersey, Piscataway, New Jersey 08854, United States
| | - Maarit Karonen
- Natural Chemistry Research Group, Department of Chemistry, University of Turku, FI-20014 Turku, Finland
| | - Mark S P Sansom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry, UMR 5086 CNRS and University of Lyon, Lyon 69367, France
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Groningen 9747 AG, The Netherlands
| |
Collapse
|
21
|
Paloncýová M, Pykal M, Kührová P, Banáš P, Šponer J, Otyepka M. Computer Aided Development of Nucleic Acid Applications in Nanotechnologies. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2022; 18:e2204408. [PMID: 36216589 DOI: 10.1002/smll.202204408] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/23/2022] [Revised: 09/12/2022] [Indexed: 06/16/2023]
Abstract
Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.
Collapse
Affiliation(s)
- Markéta Paloncýová
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Martin Pykal
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Petra Kührová
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Pavel Banáš
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Jiří Šponer
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, v. v. i., Královopolská 135, Brno, 612 65, Czech Republic
| | - Michal Otyepka
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
- IT4Innovations, VŠB - Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba, 708 00, Czech Republic
| |
Collapse
|
22
|
Shillcock JC, Lagisquet C, Alexandre J, Vuillon L, Ipsen JH. Model biomolecular condensates have heterogeneous structure quantitatively dependent on the interaction profile of their constituent macromolecules. SOFT MATTER 2022; 18:6674-6693. [PMID: 36004748 DOI: 10.1039/d2sm00387b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Biomolecular condensates play numerous roles in cells by selectively concentrating client proteins while excluding others. These functions are likely to be sensitive to the spatial organization of the scaffold proteins forming the condensate. We use coarse-grained molecular simulations to show that model intrinsically-disordered proteins phase separate into a heterogeneous, structured fluid characterized by a well-defined length scale. The proteins are modelled as semi-flexible polymers with punctate, multifunctional binding sites in good solvent conditions. Their dense phase is highly solvated with a spatial structure that is more sensitive to the separation of the binding sites than their affinity. We introduce graph theoretic measures to quantify their heterogeneity, and find that it increases with increasing binding site number, and exhibits multi-timescale dynamics. The model proteins also swell on passing from the dilute solution to the dense phase. The simulations predict that the structure of the dense phase is modulated by the location and affinity of binding sites distant from the termini of the proteins, while sites near the termini more strongly affect its phase behaviour. The relations uncovered between the arrangement of weak interaction sites on disordered proteins and the material properties of their dense phase can be experimentally tested to give insight into the biophysical properties, pathological effects, and rational design of biomolecular condensates.
Collapse
Affiliation(s)
- Julian C Shillcock
- Blue Brain Project and Laboratory of Molecular and Chemical Biology of Neurodegeneration, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland.
| | - Clément Lagisquet
- LAMA, Univ. Savoie Mont Blanc, CNRS, LAMA, 73376 Le Bourget du Lac, France.
| | - Jérémy Alexandre
- Brain Mind Institute, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Laurent Vuillon
- LAMA, Univ. Savoie Mont Blanc, CNRS, LAMA, 73376 Le Bourget du Lac, France.
| | - John H Ipsen
- Dept. of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| |
Collapse
|
23
|
Mondal S, Cui Q. Coacervation of poly-electrolytes in the presence of lipid bilayers: mutual alteration of structure and morphology. Chem Sci 2022; 13:7933-7946. [PMID: 35865903 PMCID: PMC9258347 DOI: 10.1039/d2sc02013k] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2022] [Accepted: 06/13/2022] [Indexed: 11/21/2022] Open
Abstract
Many intrinsically disordered peptides have been shown to undergo liquid-liquid phase separation and form complex coacervates, which play various regulatory roles in the cell. Recent experimental studies found that such phase separation processes may also occur at the lipid membrane surface and help organize biomolecules during signaling events; in some cases, phase separation of proteins at the membrane surface was also observed to lead to significant remodeling of the membrane morphology. The molecular mechanisms that govern the interactions between complex coacervates and lipid membranes and the impacts of such interactions on their structure and morphology, however, remain unclear. Here we study the coacervation of poly-glutamate (E30) and poly-lysine (K30) in the presence of lipid bilayers of different compositions. We carry out explicit-solvent coarse-grained molecular dynamics simulations by using the MARTINI (v3.0) force-field. We find that more than 20% anionic lipids are required for the coacervate to form stable contact with the bilayer. Upon wetting, the coacervate induces negative curvature to the bilayer and facilitates local lipid demixing, without any peptide insertion. The magnitude of negative curvature, extent of lipid demixing, and asphericity of the coacervate increase with the concentration of anionic lipids. Overall, we observe a decrease in the number of contacts among the polyelectrolytes as the droplet spreads over the bilayer. Therefore, unlike previous suggestions, interactions among polyelectrolytes do not constitute a driving force for the membrane bending upon wetting by the coacervate. Rather, analysis of interaction energy components suggests that bending of the membrane is favored by enhanced interactions between polyelectrolytes with lipids as well as with counterions. Kinetic studies reveal that, at the studied polyelectrolyte concentrations, the coacervate formation precedes bilayer wetting.
Collapse
Affiliation(s)
- Sayantan Mondal
- Department of Chemistry, Boston University 590 Commonwealth Avenue Boston MA 02215 USA (+1)-617-353-6189
| | - Qiang Cui
- Department of Chemistry, Boston University 590 Commonwealth Avenue Boston MA 02215 USA (+1)-617-353-6189
- Department of Physics, Boston University 590 Commonwealth Avenue Boston MA 02215 USA
- Department of Biomedical Engineering, Boston University 44 Cummington Mall Boston MA 02215 USA
| |
Collapse
|
24
|
Marrink SJ, Monticelli L, Melo MN, Alessandri R, Tieleman DP, Souza PCT. Two decades of Martini: Better beads, broader scope. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1620] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Siewert J. Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials University of Groningen Groningen The Netherlands
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry (MMSB ‐ UMR 5086) CNRS & University of Lyon Lyon France
| | - Manuel N. Melo
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa Oeiras Portugal
| | - Riccardo Alessandri
- Pritzker School of Molecular Engineering University of Chicago Chicago Illinois USA
| | - D. Peter Tieleman
- Centre for Molecular Simulation and Department of Biological Sciences University of Calgary Alberta Canada
| | - Paulo C. T. Souza
- Molecular Microbiology and Structural Biochemistry (MMSB ‐ UMR 5086) CNRS & University of Lyon Lyon France
| |
Collapse
|
25
|
Dai Y, Xie Z, Liang L. Pore Formation Mechanism of A-Beta Peptide on the Fluid Membrane: A Combined Coarse-Grained and All-Atomic Model. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27123924. [PMID: 35745043 PMCID: PMC9231318 DOI: 10.3390/molecules27123924] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/06/2022] [Revised: 06/05/2022] [Accepted: 06/13/2022] [Indexed: 11/27/2022]
Abstract
In Alzheimer’s disease, ion permeability through the ionic channel formed by Aβ peptides on cellular membranes appears to underlie neuronal cell death. An understanding of the formation mechanism of the toxic ionic channel by Aβ peptides is very important, but remains unclear. Our simulation results demonstrated the dynamics and mechanism of channel formation by Aβ1-28 peptides on the DPPC and POPC membrane by the coarse-grained method. The ionic channel formation is driven by the gyration of the radius and solvent accessible molecular surface area of Aβ1-28 peptides. The ionic channel formation mechanism was explored by the free energy profile based on the distribution of the gyration of the radius and solvent accessible molecular surface area of Aβ1-28 peptides on the fluid membrane. The stability and water permeability of the ionic channel formed by Aβ peptides was investigated by all-atomic model simulation. Our simulation showed that the ionic channel formed by Aβ1-28 peptides is very stable and has a good water permeability. This could help us to understand the pore formation mechanism by Aβ1-28 peptides on the fluidic membrane. It also provides us with a guideline by which to understand the toxicity of Aβ1-28 peptides’ pores to the cell.
Collapse
Affiliation(s)
- Yuxi Dai
- College of Automation, Hangzhou Dianzi University, Hangzhou 310018, China;
| | - Zhexing Xie
- College of Accounting, Hangzhou Dianzi University, Hangzhou 310018, China;
| | - Lijun Liang
- College of Automation, Hangzhou Dianzi University, Hangzhou 310018, China;
- Correspondence:
| |
Collapse
|
26
|
Thomasen FE, Pesce F, Roesgaard MA, Tesei G, Lindorff-Larsen K. Improving Martini 3 for Disordered and Multidomain Proteins. J Chem Theory Comput 2022; 18:2033-2041. [PMID: 35377637 DOI: 10.1021/acs.jctc.1c01042] [Citation(s) in RCA: 35] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Coarse-grained molecular dynamics simulations are a useful tool to determine conformational ensembles of proteins. Here, we show that the coarse-grained force field Martini 3 underestimates the global dimensions of intrinsically disordered proteins (IDPs) and multidomain proteins when compared with small-angle X-ray scattering (SAXS) data and that increasing the strength of protein-water interactions favors more expanded conformations. We find that increasing the strength of interactions between protein and water by ca. 10% results in improved agreement with the SAXS data for IDPs and multidomain proteins. We also show that this correction results in a more accurate description of self-association of IDPs and folded proteins and better agreement with paramagnetic relaxation enhancement data for most IDPs. While simulations with this revised force field still show deviations to experiments for some systems, our results suggest that it is overall a substantial improvement for coarse-grained simulations of soluble proteins.
Collapse
Affiliation(s)
- F Emil Thomasen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Francesco Pesce
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Mette Ahrensback Roesgaard
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Giulio Tesei
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Kresten Lindorff-Larsen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| |
Collapse
|
27
|
Lebold KM, Best RB. Tuning Formation of Protein-DNA Coacervates by Sequence and Environment. J Phys Chem B 2022; 126:2407-2419. [PMID: 35317553 DOI: 10.1021/acs.jpcb.2c00424] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The high concentration of nucleic acids and positively charged proteins in the cell nucleus provides many possibilities for complex coacervation. We consider a prototypical mixture of nucleic acids together with the polycationic C-terminus of histone H1 (CH1). Using a minimal coarse-grained model that captures the shape, flexibility, and charge distributions of the molecules, we find that coacervates are readily formed at physiological ionic strengths, in agreement with experiment, with a progressive increase in local ordering at low ionic strength. Variation of the positions of charged residues in the protein tunes phase separation: for well-mixed or only moderately blocky distributions of charge, there is a modest increase of local ordering with increasing blockiness that is also associated with an increased propensity to phase separate. This ordering is also associated with a slowdown of rotational and translational diffusion in the dense phase. However, for more extreme blockiness (and consequently higher local charge density), we see a qualitative change in the condensed phase to become a segregated structure with a dramatically increased ordering of the DNA. Naturally occurring proteins with these sequence properties, such as protamines in sperm cells, are found to be associated with very dense packing of nucleic acids.
Collapse
Affiliation(s)
- Kathryn M Lebold
- Laboratory of Chemical Physics, National Institutes of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Robert B Best
- Laboratory of Chemical Physics, National Institutes of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892, United States
| |
Collapse
|
28
|
Gallo P, Bachler J, Bove LE, Böhmer R, Camisasca G, Coronas LE, Corti HR, de Almeida Ribeiro I, de Koning M, Franzese G, Fuentes-Landete V, Gainaru C, Loerting T, de Oca JMM, Poole PH, Rovere M, Sciortino F, Tonauer CM, Appignanesi GA. Advances in the study of supercooled water. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2021; 44:143. [PMID: 34825973 DOI: 10.1140/epje/s10189-021-00139-1] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Accepted: 10/17/2021] [Indexed: 06/13/2023]
Abstract
In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.
Collapse
Affiliation(s)
- Paola Gallo
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy.
| | - Johannes Bachler
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | - Livia E Bove
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale A. Moro 5, 00185, Roma, Italy
- Sorbonne Université, CNRS UMR 7590, IMPMC, 75005, Paris, France
| | - Roland Böhmer
- Fakultät Physik, Technische Universität Dortmund, 44221, Dortmund, Germany
| | - Gaia Camisasca
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy
| | - Luis E Coronas
- Secció de Física Estadística i Interdisciplinària-Departament de Física de la Matèria Condensada, Universitat de Barcelona, & Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, C. Martí i Franquès 1, 08028, Barcelona, Spain
| | - Horacio R Corti
- Departamento de Física de la Materia Condensada, Centro Atómico Constituyentes, Comisión Nacional de Energía Atómica, B1650LWP, Buenos Aires, Argentina
| | - Ingrid de Almeida Ribeiro
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859, Campinas, São Paulo, Brazil
| | - Maurice de Koning
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859, Campinas, São Paulo, Brazil
- Center for Computing in Engineering & Sciences, Universidade Estadual de Campinas, UNICAMP, 13083-861, Campinas, São Paulo, Brazil
| | - Giancarlo Franzese
- Secció de Física Estadística i Interdisciplinària-Departament de Física de la Matèria Condensada, Universitat de Barcelona, & Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, C. Martí i Franquès 1, 08028, Barcelona, Spain
| | - Violeta Fuentes-Landete
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | - Catalin Gainaru
- Fakultät Physik, Technische Universität Dortmund, 44221, Dortmund, Germany
| | - Thomas Loerting
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | | | - Peter H Poole
- Department of Physics, St. Francis Xavier University, Antigonish, NS, B2G 2W5, Canada
| | - Mauro Rovere
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy
| | - Francesco Sciortino
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale A. Moro 5, 00185, Roma, Italy
| | - Christina M Tonauer
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | - Gustavo A Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000, Bahía Blanca, Argentina
| |
Collapse
|
29
|
Dhamankar S, Webb MA. Chemically specific coarse‐graining of polymers: Methods and prospects. JOURNAL OF POLYMER SCIENCE 2021. [DOI: 10.1002/pol.20210555] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Satyen Dhamankar
- Department of Chemical and Biological Engineering Princeton University Princeton New Jersey USA
| | - Michael A. Webb
- Department of Chemical and Biological Engineering Princeton University Princeton New Jersey USA
| |
Collapse
|