An official website of the United States government
Molecular dynamics (MD) simulations provide a powerful means of exploring the dynamic behavior of biomolecular systems at the atomic level. However, analyzing the vast data sets generated by MD simulations poses significant challenges. This article discusses the energy landscape visualization method (ELViM), a multidimensional reduction technique inspired by the energy landscape theory. ELViM transcends one-dimensional representations, offering a comprehensive analysis of the effective conformational phase space without the need for predefined reaction coordinates. We apply the ELViM to study the folding landscape of the antimicrobial peptide Polybia-MP1, showcasing its versatility in capturing complex biomolecular dynamics. Using dissimilarity matrices and a force-scheme approach, the ELViM provides intuitive visualizations, revealing structural correlations and local conformational signatures. The method is demonstrated to be adaptable, robust, and applicable to various biomolecular systems.
more »
« less
Latent Models of Molecular Dynamics Data: Automatic Order Parameter Generation for Peptide Fibrillization
Variational autoencoders are artificial neural networks with the capability to reduce highly dimensional sets of data to smaller dimensional, latent representations. In this work, these models are applied to molecular dynamics simulations of the self-assembly of coarse-grained peptides to obtain a singled-valued order parameter for amyloid aggregation. This automatically learned order parameter is constructed by time-averaging the latent parameterizations of internal coordinate representations and compared to the nematic order parameter which is commonly used to study ordering of similar systems in literature. It is found that the latent space value provides more tailored insight into the aggregation mechanism’s details, correctly identifying fibril formation in instances where the nematic order parameter fails to do so. A means is provided by which the latent space value can be analyzed so that the major contributing internal coordinates are identified, allowing for a direct interpretation of the latent space order parameter in terms of the behavior of the system. The latent model is found to be an effective and convenient way of representing the data from the dynamic ensemble and provides a means of reducing the dimensionality of a system whose scale exceeds molecular systems so-far considered with similar tools. This bypasses a need for researcher speculation on what elements of a system best contribute to summarizing major transitions, and suggests latent models are effective and insightful when applied to large systems with a diversity of complex behavior.
more »
« less
- Award ID(s):
- 1716956
- PAR ID:
- 10184397
- Date Published:
- Journal Name:
- The Journal of Physical Chemistry B
- ISSN:
- 1520-6106
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In high-resolution adiabatic scanning calorimetry (ASC) experiments, data for the temperature dependence of the specific enthalpy, h(T), and of the specific heat capacity, c(p)(T), are simultaneously obtained, from which the order of the phase transition and critical behaviour can be evaluated. ASC was applied to study the nematic to ferroelectric nematic phase transition (N-N-F) in the liquid crystal molecule 4-[(4-nitrophenoxy)carbonyl]phenyl 2,4-dimethoxybenzoate (RM734). The N-N-F was found to be very weakly first order with a latent heat Delta h = 0.115 +/- 0.005 J/g. The pretransitional specific heat capacity behaviour is substantially larger in the high-temperature N phase than in the low-temperature N-F phase. In both phases the power-law analysis of c(p)(T) resulted in a critical exponent alpha = 0.50 +/- 0.05 and amplitude ratio A(NF)/A(N) = 0.42 +/- 0.03. The very small latent heat and the value of alpha indicate that the N-N-F transition is close to a tricritical point. This is confirmed by a value of the order parameter exponent beta approximate to 0.25, recently obtained from electric polarisation measurements. Invoking two-scale-factor universality, it follows from the low value of A(NF)/A(N) ratio that the size of the critical fluctuations is much larger in the N-F phase than in the N phase.more » « less
-
Abstract Modern scanning probe techniques, such as scanning tunneling microscopy, provide access to a large amount of data encoding the underlying physics of quantum matter. In this work, we show how convolutional neural networks can be used to learn effective theoretical models from scanning tunneling microscopy data on correlated moiré superlattices. Moiré systems are particularly well suited for this task as their increased lattice constant provides access to intra-unit-cell physics, while their tunability allows for the collection of high-dimensional data sets from a single sample. Using electronic nematic order in twisted double-bilayer graphene as an example, we show that incorporating correlations between the local density of states at different energies allows convolutional neural networks not only to learn the microscopic nematic order parameter, but also to distinguish it from heterostrain. These results demonstrate that neural networks are a powerful method for investigating the microscopic details of correlated phenomena in moiré systems and beyond.more » « less
-
null (Ed.)Recent experiments in various cell types have shown that two-dimensional tissues often display local nematic order, with evidence of extensile stresses manifest in the dynamics of topological defects. Using a mesoscopic model where tissue flow is generated by fluctuating traction forces coupled to the nematic order parameter, we show that the resulting tissue dynamics can spontaneously produce local nematic order and an extensile internal stress. A key element of the model is the assumption that in the presence of local nematic alignment, cells preferentially crawl along the nematic axis, resulting in anisotropy of fluctuations. Our work shows that activity can drive either extensile or contractile stresses in tissue, depending on the relative strength of the contractility of the cortical cytoskeleton and tractions by cells on the extracellular matrix.more » « less
-
Abstract In a material prone to a nematic instability, anisotropic strain in principle provides a preferred symmetry-breaking direction for the electronic nematic state to follow. This is consistent with experimental observations, where electronic nematicity and structural anisotropy typically appear hand-in-hand. In this work, we discover that electronic nematicity can be locally decoupled from the underlying structural anisotropy in strain-engineered iron-selenide (FeSe) thin films. We use heteroepitaxial molecular beam epitaxy to grow FeSe with a nanoscale network of modulations that give rise to spatially varying strain. We map local anisotropic strain by analyzing scanning tunneling microscopy topographs, and visualize electronic nematic domains from concomitant spectroscopic maps. While the domains form so that the energy of nemato-elastic coupling is minimized, we observe distinct regions where electronic nematic ordering fails to flip direction, even though the underlying structural anisotropy is locally reversed. The findings point towards a nanometer-scale stiffness of the nematic order parameter.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1021/acs.jpcb.0c05763
