skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A perspective on coarse-graining methodologies for biomolecules: resolving self-assembly over extended spatiotemporal scales
The process of self-assembly of biomolecules underlies the formation of macromolecular assemblies, biomolecular materials and protein folding, and thereby is critical in many disciplines and related applications. This process typically spans numerous spatiotemporal scales and hence, is well suited for scientific interrogation via coarse-grained (CG) models used in conjunction with a suitable computational approach. This perspective provides a discussion on different coarse-graining approaches which have been used to develop CG models that resolve the process of self-assembly of biomolecules.  more » « less
Award ID(s):
1654325 2118860
PAR ID:
10503350
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Frontiers
Date Published:
Journal Name:
Frontiers in Soft Matter
Volume:
4
ISSN:
2813-0499
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Protein mimics such as peptoids form self-assembled nanostructures whose shape and function are governed by the side chain chemistry and secondary structure. Experiments have shown that a peptoid sequence with a helical secondary structure assembles into microspheres that are stable under various conditions. The conformation and organization of the peptoids within the assemblies remains unknown and is elucidated in this study via a hybrid, bottom-up coarse-graining approach. The resultant coarse-grained (CG) model preserves the chemical and structural details that are critical for capturing the secondary structure of the peptoid. The CG model accurately captures the overall conformation and solvation of the peptoids in an aqueous solution. Furthermore, the model resolves the assembly of multiple peptoids into a hemispherical aggregate that is in qualitative agreement with the corresponding results from experiments. The mildly hydrophilic peptoid residues are placed along the curved interface of the aggregate. The composition of the residues on the exterior of the aggregate is determined by two conformations adopted by the peptoid chains. Hence, the CG model simultaneously captures sequence-specific features and the assembly of a large number of peptoids. This multiscale, multiresolution coarse-graining approach could help in predicting the organization and packing of other tunable oligomeric sequences of relevance to biomedicine and electronics. 
    more » « less
  2. Coarse-grained (CG) models facilitate an efficient exploration of complex systems by reducing the unnecessary degrees of freedom of the fine-grained (FG) system while recapitulating major structural correlations. Unlike structural properties, assessing dynamic properties in CG modeling is often unfeasible due to the accelerated dynamics of the CG models, which allows for more efficient structural sampling. Therefore, the ultimate goal of the present series of articles is to establish a better correspondence between the FG and CG dynamics. To assess and compare dynamical properties in the FG and the corresponding CG models, we utilize the excess entropy scaling relationship. For Paper I of this series, we provide evidence that the FG and the corresponding CG counterpart follow the same universal scaling relationship. By carefully reviewing and examining the literature, we develop a new theory to calculate excess entropies for the FG and CG systems while accounting for entropy representability. We demonstrate that the excess entropy scaling idea can be readily applied to liquid water and methanol systems at both the FG and CG resolutions. For both liquids, we reveal that the scaling exponents remain unchanged from the coarse-graining process, indicating that the scaling behavior is universal for the same underlying molecular systems. Combining this finding with the concept of mapping entropy in CG models, we show that the missing entropy plays an important role in accelerating the CG dynamics. 
    more » « less
  3. Bottom-up coarse-grained (CG) molecular dynamics models are parameterized using complex effective Hamiltonians. These models are typically optimized to approximate high dimensional data from atomistic simulations. However, human validation of these models is often limited to low dimensional statistics that do not necessarily differentiate between the CG model and said atomistic simulations. We propose that classification can be used to variationally estimate high dimensional error and that explainable machine learning can help convey this information to scientists. This approach is demonstrated using Shapley additive explanations and two CG protein models. This framework may also be valuable for ascertaining whether allosteric effects at the atomistic level are accurately propagated to a CG model. 
    more » « less
  4. Reconstituted photosynthetic proteins which are activated upon exposure to solar energy hold enormous potential for powering future solid state devices and solar cells. The functionality and integration of these proteins into such devices has been successfully enabled by lipid-like peptides. Yet, a fundamental understanding of the organization of these peptides with respect to the photosynthetic proteins and themselves remains unknown and is critical for guiding the design of such light-activated devices. This study investigates the relative organization of one such peptide sequence V 6 K 2 (V: valine and K: lysine) within assemblies. Given the expansive spatiotemporal scales associated with this study, a hybrid coarse-grained (CG) model which captures the structure, conformation and aggregation of the peptide is adopted. The CG model uses a combination of iterative Boltzmann inversion and force matching to provide insight into the relative organization of V 6 K 2 in assemblies. The CG model reproduces the structure of a V 6 K 2 peptide sequence along with its all atom (AA) solvation structure. The relative organization of multiple peptides in an assembly, as captured by CG simulations, is in agreement with corresponding results from AA simulations. Also, a backmapping procedure reintroduces the AA details of the peptides within the aggregates captured by the CG model to demonstrate the relative organization of the peptides. Furthermore, a large number of peptides self-assemble into an elongated micelle in the CG simulation, which is consistent with experimental findings. The coarse-graining procedure is tested for transferability to longer peptide sequences, and hence can be extended to other amphiphilic peptide sequences. 
    more » « less
  5. Biomolecular adsorption has great significance in medical, environmental, and technological processes. Understanding adsorption equilibrium and binding kinetics is essential for advanced process implementation. This requires identifying intrinsic determinants that predict optimal adsorption properties at bio–hybrid interfaces. Solid-binding peptides (SBPs) have targetable intrinsic properties involving peptide–peptide and peptide–solid interactions, which result in high-affinity material-selective binding. Atomic force microscopy investigations confirmed this complex interplay of multi-step peptide assemblies in a cooperative modus. Yet, most studies report adsorption properties of SBPs using non-cooperative or single-step adsorption models. Using non-cooperative kinetic models for predicting cooperative self-assembly behavior creates an oversimplified view of peptide adsorption, restricting implementing SBPs beyond their current use. To address these limitations and provide insight into surface-level events during self-assembly, a novel method, the Frequency Response Cooperativity model, was developed. This model iteratively fits adsorption data through spectral analysis of several time-dependent kinetic parameters. The model, applied to a widely used gold-binding peptide data obtained using a quartz crystal microbalance with dissipation, verified multi-step assembly. Peak deconvolution of spectral plots revealed distinct differences in the size and distribution of the kinetic rates present during adsorption across the concentrations. This approach provides new fundamental insights into the intricate dynamics of self-assembly of biomolecules on surfaces. 
    more » « less