1
|
Dong XY, Liu R, Seroski DT, Hudalla GA, Hall CK. Programming co-assembled peptide nanofiber morphology via anionic amino acid type: Insights from molecular dynamics simulations. PLoS Comput Biol 2023; 19:e1011685. [PMID: 38048311 PMCID: PMC10729967 DOI: 10.1371/journal.pcbi.1011685] [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: 06/13/2023] [Revised: 12/19/2023] [Accepted: 11/13/2023] [Indexed: 12/06/2023] Open
Abstract
Co-assembling peptides can be crafted into supramolecular biomaterials for use in biotechnological applications, such as cell culture scaffolds, drug delivery, biosensors, and tissue engineering. Peptide co-assembly refers to the spontaneous organization of two different peptides into a supramolecular architecture. Here we use molecular dynamics simulations to quantify the effect of anionic amino acid type on co-assembly dynamics and nanofiber structure in binary CATCH(+/-) peptide systems. CATCH peptide sequences follow a general pattern: CQCFCFCFCQC, where all C's are either a positively charged or a negatively charged amino acid. Specifically, we investigate the effect of substituting aspartic acid residues for the glutamic acid residues in the established CATCH(6E-) molecule, while keeping CATCH(6K+) unchanged. Our results show that structures consisting of CATCH(6K+) and CATCH(6D-) form flatter β-sheets, have stronger interactions between charged residues on opposing β-sheet faces, and have slower co-assembly kinetics than structures consisting of CATCH(6K+) and CATCH(6E-). Knowledge of the effect of sidechain type on assembly dynamics and fibrillar structure can help guide the development of advanced biomaterials and grant insight into sequence-to-structure relationships.
Collapse
Affiliation(s)
- Xin Y. Dong
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina, United States of America
| | - Renjie Liu
- Department of Biomedical Engineering, University of Florida, Gainesville, Florida, United States of America
| | - Dillon T. Seroski
- Department of Biomedical Engineering, University of Florida, Gainesville, Florida, United States of America
| | - Gregory A. Hudalla
- Department of Biomedical Engineering, University of Florida, Gainesville, Florida, United States of America
| | - Carol K. Hall
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina, United States of America
| |
Collapse
|
2
|
Franke L, Peter C. Visualizing the Residue Interaction Landscape of Proteins by Temporal Network Embedding. J Chem Theory Comput 2023; 19:2985-2995. [PMID: 37122117 DOI: 10.1021/acs.jctc.2c01228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
Characterizing the structural dynamics of proteins with heterogeneous conformational landscapes is crucial to understanding complex biomolecular processes. To this end, dimensionality reduction algorithms are used to produce low-dimensional embeddings of the high-dimensional conformational phase space. However, identifying a compact and informative set of input features for the embedding remains an ongoing challenge. Here, we propose to harness the power of Residue Interaction Networks (RINs) and their centrality measures, established tools to provide a graph theoretical view on molecular structure. Specifically, we combine the closeness centrality, which captures global features of the protein conformation at residue-wise resolution, with EncoderMap, a hybrid neural-network autoencoder/multidimensional-scaling like dimensionality reduction algorithm. We find that the resulting low-dimensional embedding is a meaningful visualization of the residue interaction landscape that resolves structural details of the protein behavior while retaining global interpretability. This feature-based graph embedding of temporal protein graphs makes it possible to apply the general descriptive power of RIN formalisms to the analysis of protein simulations of complex processes such as protein folding and multidomain interactions requiring no protein-specific input. We demonstrate this on simulations of the fast folding protein Trp-Cage and the multidomain signaling protein FAT10. Due to its generality and modularity, the presented approach can easily be transferred to other protein systems.
Collapse
Affiliation(s)
- Leon Franke
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, Konstanz 78457, Germany
- Konstanz Research School Chemical Biology, University of Konstanz, Universitätsstraße 10, Konstanz 78457, Germany
| | - Christine Peter
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, Konstanz 78457, Germany
| |
Collapse
|
3
|
Hunkler S, Buhl T, Kukharenko O, Peter C. Generating a conformational landscape of ubiquitin chains at atomistic resolution by back-mapping based sampling. Front Chem 2023; 10:1087963. [PMID: 36704619 PMCID: PMC9871295 DOI: 10.3389/fchem.2022.1087963] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Accepted: 12/23/2022] [Indexed: 01/12/2023] Open
Abstract
Ubiquitin chains are flexible multidomain proteins that have important biological functions in cellular signalling. Computational studies with all-atom molecular dynamics simulations of the conformational spaces of polyubiquitins can be challenging due to the system size and a multitude of long-lived meta-stable states. Coarse graining is an efficient approach to overcome this problem-at the cost of losing high-resolution details. Recently, we proposed the back-mapping based sampling (BMBS) approach that reintroduces atomistic information into a given coarse grained (CG) sampling based on a two-dimensional (2D) projection of the conformational landscape, produces an atomistic ensemble and allows to systematically compare the ensembles at the two levels of resolution. Here, we apply BMBS to K48-linked tri-ubiquitin, showing its applicability to larger systems than those it was originally introduced on and demonstrating that the algorithm scales very well with system size. In an extension of the original BMBS we test three different seeding strategies, i.e. different approaches from where in the CG landscape atomistic trajectories are initiated. Furthermore, we apply a recently introduced conformational clustering algorithm to the back-mapped atomistic ensemble. Thus, we obtain insight into the structural composition of the 2D landscape and illustrate that the dimensionality reduction algorithm separates different conformational characteristics very well into different regions of the map. This cluster analysis allows us to show how atomistic trajectories sample conformational states, move through the projection space and in sum converge to an atomistic conformational landscape that slightly differs from the original CG map, indicating a correction of flaws in the CG template.
Collapse
Affiliation(s)
- Simon Hunkler
- Department of Chemistry, University of Konstanz, Konstanz, Germany
| | - Teresa Buhl
- Department of Chemistry, University of Konstanz, Konstanz, Germany
| | | | - Christine Peter
- Department of Chemistry, University of Konstanz, Konstanz, Germany,*Correspondence: Christine Peter,
| |
Collapse
|
4
|
Blasius J, Zaby P, Dölz J, Kirchner B. Uncertainty quantification of phase transition quantities from cluster weighting calculations. J Chem Phys 2022; 157:014505. [DOI: 10.1063/5.0093057] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we investigate how uncertainties in experimental input data influence the results of quantum cluster equilibrium calculations. In particular, we focus on the calculation of vaporization enthalpies and entropies of seven organic liquids, compare two computational approaches for their calculation and investigate how these properties are affected by changes in the experimental input data. It is observed that the vaporization enthalpies and entropies show a smooth dependence on changes in the reference density and boiling point. The reference density is found to have only a small influence on the vaporization thermodynamics, whereas the boiling point has a large influence on the vaporization enthalpy but only a small influence on the vaporization entropy. Furthermore we employed the Gauss--Hermite estimator in order to quantify the uncertainty in the thermodynamic functions that stems from inaccuracies in the experimental reference data at the example of the vaporization enthalpy of (\textit{R})-butan-2-ol. We quantify the uncertainty as 30.95~$\cdot$10$^{-3}$~kJ~mol$^{-1}$. Additionally we compare the convergence behaviour and computational effort of the Gauss--Hermite estimator with the Monte Carlo approach and show the superiority of the former. By this, we present how uncertainty quantification can be applied to examples from theoretical chemistry.
Collapse
Affiliation(s)
- Jan Blasius
- University of Bonn Institute of Physical and Theoretical Chemistry, Germany
| | - Paul Zaby
- University of Bonn Institute of Physical and Theoretical Chemistry, Germany
| | - Jürgen Dölz
- Institute for Numerical Simulation, University of Bonn, Friedrich-Hirzebruch-Allee 7 53115 Bonn, Germany, Germany
| | - Barbara Kirchner
- Theoretical and Physical Chemistry, University of Bonn Institute of Physical and Theoretical Chemistry, Germany
| |
Collapse
|
5
|
Lemke T, Edte M, Gebauer D, Peter C. Three Reasons Why Aspartic Acid and Glutamic Acid Sequences Have a Surprisingly Different Influence on Mineralization. J Phys Chem B 2021; 125:10335-10343. [PMID: 34473925 DOI: 10.1021/acs.jpcb.1c04467] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding the role of polymers rich in aspartic acid (Asp) and glutamic acid (Glu) is the key to gaining precise control over mineralization processes. Despite their chemical similarity, experiments revealed a surprisingly different influence of Asp and Glu sequences. We conducted molecular dynamics simulations of Asp and Glu peptides in the presence of calcium and chloride ions to elucidate the underlying phenomena. In line with experimental differences, in our simulations, we indeed find strong differences in the way the peptides interact with ions in solution. The investigated Asp pentapeptide tends to pull a lot of ions into its vicinity, and many structures with clusters of calcium and chloride ions on the surface of the peptide can be observed. Under the same conditions, comparatively fewer ions can be found in proximity of the investigated Glu pentapeptide, and the structures are characterized by single calcium ions bound to multiple carboxylate groups. Based on our simulation data, we identified three reasons contributing to these differences, leading to a new level of understanding additive-ion interactions.
Collapse
Affiliation(s)
- Tobias Lemke
- Theoretical Chemistry, University of Konstanz, 78547 Konstanz, Germany
| | - Moritz Edte
- Theoretical Chemistry, University of Konstanz, 78547 Konstanz, Germany
| | - Denis Gebauer
- Institute of Inorganic Chemistry, Leibniz University Hannover, 30167 Hannover, Germany
| | - Christine Peter
- Theoretical Chemistry, University of Konstanz, 78547 Konstanz, Germany
| |
Collapse
|
6
|
Kubincová A, Riniker S, Hünenberger PH. Solvent-scaling as an alternative to coarse-graining in adaptive-resolution simulations: The adaptive solvent-scaling (AdSoS) scheme. J Chem Phys 2021; 155:094107. [PMID: 34496576 DOI: 10.1063/5.0057384] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new approach termed Adaptive Solvent-Scaling (AdSoS) is introduced for performing simulations of a solute embedded in a fine-grained (FG) solvent region itself surrounded by a coarse-grained (CG) solvent region, with a continuous FG ↔ CG switching of the solvent resolution across a buffer layer. Instead of relying on a distinct CG solvent model, the AdSoS scheme is based on CG models defined by a dimensional scaling of the FG solvent by a factor s, accompanied by an s-dependent modulation of the atomic masses and interaction parameters. The latter changes are designed to achieve an isomorphism between the dynamics of the FG and CG models, and to preserve the dispersive and dielectric solvation properties of the solvent with respect to a solute at FG resolution. This scaling approach offers a number of advantages compared to traditional coarse-graining: (i) the CG parameters are immediately related to those of the FG model (no need to parameterize a distinct CG model); (ii) nearly ideal mixing is expected for CG variants with similar s-values (ideal mixing holding in the limit of identical s-values); (iii) the solvent relaxation timescales should be preserved (no dynamical acceleration typical for coarse-graining); (iv) the graining level NG (number of FG molecules represented by one CG molecule) can be chosen arbitrarily (in particular, NG = s3 is not necessarily an integer); and (v) in an adaptive-resolution scheme, this level can be varied continuously as a function of the position (without requiring a bundling mechanism), and this variation occurs at a constant number of particles per molecule (no occurrence of fractional degrees of freedom in the buffer layer). By construction, the AdSoS scheme minimizes the thermodynamic mismatch between the different regions of the adaptive-resolution system, leading to a nearly homogeneous scaled solvent density s3ρ. Residual density artifacts in and at the surface of the boundary layer can easily be corrected by means of a grid-based biasing potential constructed in a preliminary pure-solvent simulation. This article introduces the AdSoS scheme and provides an initial application to pure atomic liquids (no solute) with Lennard-Jones plus Coulomb interactions in a slab geometry.
Collapse
Affiliation(s)
- Alžbeta Kubincová
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zurich, Vladimir Prelog-Weg 2, CH-8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zurich, Vladimir Prelog-Weg 2, CH-8093 Zürich, Switzerland
| | - Philippe H Hünenberger
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zurich, Vladimir Prelog-Weg 2, CH-8093 Zürich, Switzerland
| |
Collapse
|
7
|
Chen M. Collective variable-based enhanced sampling and machine learning. THE EUROPEAN PHYSICAL JOURNAL. B 2021; 94:211. [PMID: 34697536 PMCID: PMC8527828 DOI: 10.1140/epjb/s10051-021-00220-w] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/02/2021] [Accepted: 10/03/2021] [Indexed: 05/14/2023]
Abstract
ABSTRACT Collective variable-based enhanced sampling methods have been widely used to study thermodynamic properties of complex systems. Efficiency and accuracy of these enhanced sampling methods are affected by two factors: constructing appropriate collective variables for enhanced sampling and generating accurate free energy surfaces. Recently, many machine learning techniques have been developed to improve the quality of collective variables and the accuracy of free energy surfaces. Although machine learning has achieved great successes in improving enhanced sampling methods, there are still many challenges and open questions. In this perspective, we shall review recent developments on integrating machine learning techniques and collective variable-based enhanced sampling approaches. We also discuss challenges and future research directions including generating kinetic information, exploring high-dimensional free energy surfaces, and efficiently sampling all-atom configurations.
Collapse
Affiliation(s)
- Ming Chen
- Department of Chemistry, Purdue University, West Lafayette, IN 47907 USA
| |
Collapse
|
8
|
Shen K, Sherck N, Nguyen M, Yoo B, Köhler S, Speros J, Delaney KT, Fredrickson GH, Shell MS. Learning composition-transferable coarse-grained models: Designing external potential ensembles to maximize thermodynamic information. J Chem Phys 2020; 153:154116. [DOI: 10.1063/5.0022808] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Affiliation(s)
- Kevin Shen
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, USA
| | - Nicholas Sherck
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| | - My Nguyen
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| | - Brian Yoo
- BASF Corporation, Tarrytown, New York 10591, USA
| | | | - Joshua Speros
- California Research Alliance (CARA) by BASF, Berkeley, California 94720, USA
| | - Kris T. Delaney
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, USA
| | - Glenn H. Fredrickson
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, USA
- Department of Materials Engineering, University of California, Santa Barbara, California 93106, USA
| | - M. Scott Shell
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| |
Collapse
|
9
|
Wu K, Xu S, Wan B, Xiu P, Zhou X. A novel multiscale scheme to accelerate atomistic simulations of bio-macromolecules by adaptively driving coarse-grained coordinates. J Chem Phys 2020; 152:114115. [PMID: 32199430 DOI: 10.1063/1.5135309] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
All-atom molecular dynamics (MD) simulations of bio-macromolecules can yield relatively accurate results while suffering from the limitation of insufficient conformational sampling. On the other hand, the coarse-grained (CG) MD simulations efficiently accelerate conformational changes in biomolecules but lose atomistic details and accuracy. Here, we propose a novel multiscale simulation method called the adaptively driving multiscale simulation (ADMS)-it efficiently accelerates biomolecular dynamics by adaptively driving virtual CG atoms on the fly while maintaining the atomistic details and focusing on important conformations of the original system with irrelevant conformations rarely sampled. Herein, the "adaptive driving" is based on the short-time-averaging response of the system (i.e., an approximate free energy surface of the original system), without requiring the construction of the CG force field. We apply the ADMS to two peptides (deca-alanine and Ace-GGPGGG-Nme) and one small protein (HP35) as illustrations. The simulations show that the ADMS not only efficiently captures important conformational states of biomolecules and drives fast interstate transitions but also yields, although it might be in part, reliable protein folding pathways. Remarkably, a ∼100-ns explicit-solvent ADMS trajectory of HP35 with three CG atoms realizes folding and unfolding repeatedly and captures the important states comparable to those from a 398-µs standard all-atom MD simulation.
Collapse
Affiliation(s)
- Kai Wu
- Department of Engineering Mechanics, Zhejiang University, Hangzhou 310027, China
| | - Shun Xu
- Computer Network Information Center, Chinese Academy of Sciences, Beijing 100190, China
| | - Biao Wan
- Beijing Computational Science Research Center, Beijing 1100193, China
| | - Peng Xiu
- Department of Engineering Mechanics, Zhejiang University, Hangzhou 310027, China
| | - Xin Zhou
- School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
10
|
Gkeka P, Stoltz G, Barati Farimani A, Belkacemi Z, Ceriotti M, Chodera JD, Dinner AR, Ferguson AL, Maillet JB, Minoux H, Peter C, Pietrucci F, Silveira A, Tkatchenko A, Trstanova Z, Wiewiora R, Lelièvre T. Machine Learning Force Fields and Coarse-Grained Variables in Molecular Dynamics: Application to Materials and Biological Systems. J Chem Theory Comput 2020; 16:4757-4775. [PMID: 32559068 PMCID: PMC8312194 DOI: 10.1021/acs.jctc.0c00355] [Citation(s) in RCA: 82] [Impact Index Per Article: 20.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Machine learning encompasses tools and algorithms that are now becoming popular in almost all scientific and technological fields. This is true for molecular dynamics as well, where machine learning offers promises of extracting valuable information from the enormous amounts of data generated by simulation of complex systems. We provide here a review of our current understanding of goals, benefits, and limitations of machine learning techniques for computational studies on atomistic systems, focusing on the construction of empirical force fields from ab initio databases and the determination of reaction coordinates for free energy computation and enhanced sampling.
Collapse
Affiliation(s)
- Paraskevi Gkeka
- Integrated Drug Discovery, Sanofi R&D, 91385 Chilly-Mazarin, France
| | - Gabriel Stoltz
- CERMICS, Ecole des Ponts, Marne-la-Vallée, France
- Matherials Project-Team, Inria Paris, 75012 Paris, France
| | | | - Zineb Belkacemi
- Integrated Drug Discovery, Sanofi R&D, 91385 Chilly-Mazarin, France
- CERMICS, Ecole des Ponts, Marne-la-Vallée, France
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - John D Chodera
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, United States
| | - Aaron R Dinner
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, United States
| | | | - Hervé Minoux
- Integrated Drug Discovery, Sanofi R&D, 94403 Vitry-sur-Seine, France
| | | | - Fabio Pietrucci
- UMR CNRS 7590, MNHN, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Sorbonne Université, 75005 Paris, France
| | - Ana Silveira
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, United States
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Zofia Trstanova
- School of Mathematics, The University of Edinburgh, Edinburgh EH9 3FD, U.K
| | - Rafal Wiewiora
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, United States
| | - Tony Lelièvre
- CERMICS, Ecole des Ponts, Marne-la-Vallée, France
- Matherials Project-Team, Inria Paris, 75012 Paris, France
| |
Collapse
|
11
|
Li W, Burkhart C, Polińska P, Harmandaris V, Doxastakis M. Backmapping coarse-grained macromolecules: An efficient and versatile machine learning approach. J Chem Phys 2020; 153:041101. [DOI: 10.1063/5.0012320] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Wei Li
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Craig Burkhart
- The Goodyear Tire and Rubber Company, Akron, Ohio 44305, USA
| | - Patrycja Polińska
- Goodyear Innovation Center Luxembourg, Avenue Gordon Smith, L-7750 Colmar-Berg, Luxembourg
| | - Vagelis Harmandaris
- Department of Applied Mathematics, University of Crete, and IACM FORTH, GR-71110 Heraklion, Greece
- Computation-Based Science and Technology Research Center, The Cyprus Institute, Nicosia 2121, Cyprus
| | - Manolis Doxastakis
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, Tennessee 37996, USA
| |
Collapse
|