Free energy profile and mechanism of self-assembly of peptide amphiphiles based on a collective assembly coordinate

Research output: Contribution to journalArticle

19 Citations (Scopus)

Abstract

By combining targeted molecular dynamics (TMD) simulations, umbrella sampling, and the weighted histogram analysis method (WHAM), we have calculated the potential of mean force (PMF) for the transition between the bound and free states of 90 peptide amphiphiles (PAs) in aqueous solution, with the bound state corresponding to a cylindrical micelle fiber. We specifically consider a collective reaction coordinate, the radius of gyration of the PAs, to describe assembly in this work. It is found that the free energy, enthalpy, and entropy differences between the free and bound states are -126 kcal/mol, -185 kcal/mol, and -190 cal/(mol K), respectively, for the self-assembly process. This indicates that the driving force to form the micelle structure is enthalpic. The enthalpic driving forces originate from several factors, including the conformational energy of PAs and the electrostatic and van der Waals interaction energy between solvent molecules and between solvent and PAs. Among these interactions, the solvent electrostatic interaction is the dominating one, contributing 54% of the total driving force. The PMF profile can be recognized as involving two stages of assembly: (1) PAs initially approach each other in mostly random configurations and loosely aggregate, resulting in significant desolvation and initiation of head-tail conformational reorganization; (2) PAs undergo a conformational disorder-to-order transition, including forming secondary structures that include more β-sheets and fewer random coils, along with tail-head core-shell alignment and condensation that leads to total exclusion of water from the core. The PMF decreases slowly in the first stage, but rapidly in the second. This study demonstrates a hierarchy of assembly steps in which PA structural changes, solvation, and redistribution of solvent molecules play significant roles in the PA self-assembly process.

Original languageEnglish
Pages (from-to)9004-9013
Number of pages10
JournalJournal of Physical Chemistry B
Volume117
Issue number30
DOIs
Publication statusPublished - Aug 1 2013

Fingerprint

Amphiphiles
Self assembly
Peptides
Free energy
peptides
self assembly
assembly
free energy
profiles
Micelles
micelles
electrostatics
Molecules
gyration
Solvation
interactions
Coulomb interactions
exclusion
histograms
hierarchies

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry
  • Materials Chemistry
  • Surfaces, Coatings and Films

Cite this

Free energy profile and mechanism of self-assembly of peptide amphiphiles based on a collective assembly coordinate. / Yu, Tao; Schatz, George C.

In: Journal of Physical Chemistry B, Vol. 117, No. 30, 01.08.2013, p. 9004-9013.

Research output: Contribution to journalArticle

@article{740f022f96df4ccea295cacef9b8497f,
title = "Free energy profile and mechanism of self-assembly of peptide amphiphiles based on a collective assembly coordinate",
abstract = "By combining targeted molecular dynamics (TMD) simulations, umbrella sampling, and the weighted histogram analysis method (WHAM), we have calculated the potential of mean force (PMF) for the transition between the bound and free states of 90 peptide amphiphiles (PAs) in aqueous solution, with the bound state corresponding to a cylindrical micelle fiber. We specifically consider a collective reaction coordinate, the radius of gyration of the PAs, to describe assembly in this work. It is found that the free energy, enthalpy, and entropy differences between the free and bound states are -126 kcal/mol, -185 kcal/mol, and -190 cal/(mol K), respectively, for the self-assembly process. This indicates that the driving force to form the micelle structure is enthalpic. The enthalpic driving forces originate from several factors, including the conformational energy of PAs and the electrostatic and van der Waals interaction energy between solvent molecules and between solvent and PAs. Among these interactions, the solvent electrostatic interaction is the dominating one, contributing 54{\%} of the total driving force. The PMF profile can be recognized as involving two stages of assembly: (1) PAs initially approach each other in mostly random configurations and loosely aggregate, resulting in significant desolvation and initiation of head-tail conformational reorganization; (2) PAs undergo a conformational disorder-to-order transition, including forming secondary structures that include more β-sheets and fewer random coils, along with tail-head core-shell alignment and condensation that leads to total exclusion of water from the core. The PMF decreases slowly in the first stage, but rapidly in the second. This study demonstrates a hierarchy of assembly steps in which PA structural changes, solvation, and redistribution of solvent molecules play significant roles in the PA self-assembly process.",
author = "Tao Yu and Schatz, {George C}",
year = "2013",
month = "8",
day = "1",
doi = "10.1021/jp404835q",
language = "English",
volume = "117",
pages = "9004--9013",
journal = "Journal of Physical Chemistry B",
issn = "1520-6106",
publisher = "American Chemical Society",
number = "30",

}

TY - JOUR

T1 - Free energy profile and mechanism of self-assembly of peptide amphiphiles based on a collective assembly coordinate

AU - Yu, Tao

AU - Schatz, George C

PY - 2013/8/1

Y1 - 2013/8/1

N2 - By combining targeted molecular dynamics (TMD) simulations, umbrella sampling, and the weighted histogram analysis method (WHAM), we have calculated the potential of mean force (PMF) for the transition between the bound and free states of 90 peptide amphiphiles (PAs) in aqueous solution, with the bound state corresponding to a cylindrical micelle fiber. We specifically consider a collective reaction coordinate, the radius of gyration of the PAs, to describe assembly in this work. It is found that the free energy, enthalpy, and entropy differences between the free and bound states are -126 kcal/mol, -185 kcal/mol, and -190 cal/(mol K), respectively, for the self-assembly process. This indicates that the driving force to form the micelle structure is enthalpic. The enthalpic driving forces originate from several factors, including the conformational energy of PAs and the electrostatic and van der Waals interaction energy between solvent molecules and between solvent and PAs. Among these interactions, the solvent electrostatic interaction is the dominating one, contributing 54% of the total driving force. The PMF profile can be recognized as involving two stages of assembly: (1) PAs initially approach each other in mostly random configurations and loosely aggregate, resulting in significant desolvation and initiation of head-tail conformational reorganization; (2) PAs undergo a conformational disorder-to-order transition, including forming secondary structures that include more β-sheets and fewer random coils, along with tail-head core-shell alignment and condensation that leads to total exclusion of water from the core. The PMF decreases slowly in the first stage, but rapidly in the second. This study demonstrates a hierarchy of assembly steps in which PA structural changes, solvation, and redistribution of solvent molecules play significant roles in the PA self-assembly process.

AB - By combining targeted molecular dynamics (TMD) simulations, umbrella sampling, and the weighted histogram analysis method (WHAM), we have calculated the potential of mean force (PMF) for the transition between the bound and free states of 90 peptide amphiphiles (PAs) in aqueous solution, with the bound state corresponding to a cylindrical micelle fiber. We specifically consider a collective reaction coordinate, the radius of gyration of the PAs, to describe assembly in this work. It is found that the free energy, enthalpy, and entropy differences between the free and bound states are -126 kcal/mol, -185 kcal/mol, and -190 cal/(mol K), respectively, for the self-assembly process. This indicates that the driving force to form the micelle structure is enthalpic. The enthalpic driving forces originate from several factors, including the conformational energy of PAs and the electrostatic and van der Waals interaction energy between solvent molecules and between solvent and PAs. Among these interactions, the solvent electrostatic interaction is the dominating one, contributing 54% of the total driving force. The PMF profile can be recognized as involving two stages of assembly: (1) PAs initially approach each other in mostly random configurations and loosely aggregate, resulting in significant desolvation and initiation of head-tail conformational reorganization; (2) PAs undergo a conformational disorder-to-order transition, including forming secondary structures that include more β-sheets and fewer random coils, along with tail-head core-shell alignment and condensation that leads to total exclusion of water from the core. The PMF decreases slowly in the first stage, but rapidly in the second. This study demonstrates a hierarchy of assembly steps in which PA structural changes, solvation, and redistribution of solvent molecules play significant roles in the PA self-assembly process.

UR - http://www.scopus.com/inward/record.url?scp=84881256326&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84881256326&partnerID=8YFLogxK

U2 - 10.1021/jp404835q

DO - 10.1021/jp404835q

M3 - Article

VL - 117

SP - 9004

EP - 9013

JO - Journal of Physical Chemistry B

JF - Journal of Physical Chemistry B

SN - 1520-6106

IS - 30

ER -