Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conformation dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Protein folding from heterogeneous unfolded state revealed by time-resolved X-ray solution scattering

One of the most challenging tasks in biological science is to understand how a protein folds. In theoretical studies, the hypothesis adopting a funnel-like free-energy landscape has been recognized as a prominent scheme for explaining protein folding in views of both internal energy and conformational heterogeneity of a protein. Despite numerous experimental efforts, however, comprehensively studying protein folding with respect to its global conformational changes in conjunction with the heterogeneity has been elusive. Here we investigate the redox-coupled folding dynamics of equine heart cytochrome c (cyt-c) induced by external electron injection by using time-resolved X-ray solution scattering. A systematic kinetic analysis unveils a kinetic model for its folding with a stretched exponential behavior during the transition toward the folded state. With the aid of the ensemble optimization method combined with molecular dynamics simulations, we found that during the folding the heterogeneously populated ensemble of the unfolded state is converted to a narrowly populated ensemble of folded conformations. These observations obtained from the kinetic and the structural analyses of X-ray scattering data reveal that the folding dynamics of cyt-c accompanies many parallel pathways associated with the heterogeneously populated ensemble of unfolded conformations, resulting in the stretched exponential kinetics at room temperature. This finding provides direct evidence with a view to microscopic protein conformations that the cyt-c folding initiates from a highly heterogeneous unfolded state, passes through still diverse intermediate structures, and reaches structural homogeneity by arriving at the folded state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Investigation of the Effect of Temperature on the Structure of SARS-CoV-2 Spike Protein by Molecular Dynamics Simulations

Statistical and epidemiological data imply temperature sensitivity of the SARS-CoV2 coronavirus. However, the molecular level understanding of the virus structure at different temperature is still not clear. Spike protein is the outermost structural protein of the SARS-CoV-2 virus which interacts with the Angiotensin Converting Enzyme 2 (ACE2), a human receptor, and enters the respiratory system. In this study, we performed an all atom molecular dynamics simulation to study the effect of temperature on the structure of the Spike protein. After 200 ns of simulation at different temperatures, we came across some interesting phenomena exhibited by the protein. We found that the solvent exposed domain of Spike protein, namely S1, is more mobile than the transmembrane domain, S2. Structural studies implied the presence of several charged residues on the surface of N-terminal Domain of S1 which are optimally oriented at 10– 30°C. Bioinformatics analyses indicated that it is capable of binding to other human receptors and should not be disregarded. Additionally, we found that receptor binding motif (RBM), present on the receptor binding domain (RBD) of S1, begins to close around temperature of 40°C and attains a completely closed conformation at 50°C. We also found that the presence of glycan moieties did not influence the observed protein dynamics. Nevertheless, the closed conformation disables its ability to bind to ACE2, due to the burying of its receptor binding residues. Our results clearly show that there are active and inactive states of the protein at different temperatures. This would not only prove beneficial for understanding the fundamental nature of the virus, but would be also useful in the development of vaccines and therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Microscopic insights into the solvation of polyethylene glycol chains in water: A machine learning potential approach

Polyethylene glycol (PEG) is a structurally simple, nontoxic, and water-soluble polymer widely utilized in medical and pharmaceutical applications. Notably, when a PEG chain is immersed in water, the surrounding water molecules play a key role in driving conformational changes of this macromolecule. In this study, we explore the solvation behavior of PEG under mechanical strain using molecular dynamics simulations, with an interatomic potential obtained from machine learning. Our focus is on the transition from the favored coil-like conformation to an extended one under external force. Through analyses of radial distribution functions, hydrogen bonding, and solvation dynamics, we uncover how mechanical stretching influences the local hydration environment. Furthermore, we disentangle the enthalpic and entropic contributions to the conformational stability of PEG in water. Surprisingly, our neural network potential model identifies dewetting of PEG C-atoms, and not water H-bonding with PEG O-atoms, as the main enthalpic driving force for the coiling of PEG in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Coordination Sphere of Lanthanide Aqua Ions Resolved with Ab Initio Molecular Dynamics and X-ray Absorption Spectroscopy

To resolve the fleeting structures of lanthanide Ln 3+ aqua ions in solution, we (i) performed the first ab initio molecular dynamics (AIMD) simulations of the entire series of Ln 3+ aqua ions in explicit water solvent using pseudopotentials and basis sets recently optimized for lanthanides and (ii) measured the symmetry of the hydrating waters about Ln 3+ ions (Nd 3+ , Dy 3+ , Er 3+ , Lu 3+ ) for the first time with extended X-ray absorption fine structure (EXAFS). EXAFS spectra were measured experimentally and generated from AIMD trajectories to directly compare simulation, which concurrently considers the electronic structure and the atomic dynamics in solution, with experiment. In this study, we performed a comprehensive evaluation of EXAFS multiple-scattering analysis (up to 6.5 Å) to measure Ln–O distances and angular correlations (i.e., symmetry) and elucidate the molecular geometry of the first hydration shell. This evaluation, in combination with symmetry-dependent L 3 - and L 1 -edge spectral analysis, shows that the AIMD simulations remarkably reproduces the experimental EXAFS data. The error in the predicted Ln–O distances is less than 0.07 Å for the later lanthanides, while we observed excellent agreement with predicted distances within experimental uncertainty for the early lanthanides. Our analysis revealed a dynamic, symmetrically disordered first coordination shell, which does not conform to a single molecular geometry for most lanthanides. This work sheds critical light on the highly elusive coordination geometry of the Ln 3+ aqua ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Uncovering Backbone Conformation for Rigid DPP-Based Donor–Acceptor Conjugated Polymer Using Deuterium Labeling and Neutron Scattering

The conjugated polymer’s backbone conformation dictates the delocalization of electrons, ultimately affecting its optoelectronic properties. Most conjugated polymers can be viewed as semirigid rods with their backbone embedded among long alkyl side chains. Thus, it is challenging to experimentally quantify the conformation of a conjugated backbone. Here, we performed contrast variation neutron scattering on rigid conjugated donor–acceptor (D–A) diketopyrrolopyrrole (DPP) polymers with selectively deuterated side chains to measure the conjugated backbone conformation. We first synthesized DPPbased polymers with deuterated side chains, confirmed by NMR and FTIR. Using contrast variation neutron scattering, we found that the DPP-based conjugated polymers are much more rigid than poly(3-alkylthiophenes), with persistence length (L p ) at 16–18 nm versus 2–3 nm. More importantly, in contrast to the relatively flexible poly(3-alkylthiophenes) whose backbone is more flexible than the whole polymer, we found that the backbone of DPP-based polymers has the same L p value compared to the whole polymer chain. This indicates that side chain interference on backbone conformation is not present for the semirigid polymer, which is further confirmed by coarse-grained molecular dynamics (CG-MD) simulations. Our work provides a novel protocol to probe polymer’s backbone conformation and paradigm-shifting understanding of the backbone conformation of semirigid conjugated polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Short-wavelength-sensitive 2 (Sws2) visual photopigment models combined with atomistic molecular simulations to predict spectral peaks of absorbance

For many species, vision is one of the most important sensory modalities for mediating essential tasks that include navigation, predation and foraging, predator avoidance, and numerous social behaviors. The vertebrate visual process begins when photons of the light interact with rod and cone photoreceptors that are present in the neural retina. Vertebrate visual photopigments are housed within these photoreceptor cells and are sensitive to a wide range of wavelengths that peak within the light spectrum, the latter of which is a function of the type of chromophore used and how it interacts with specific amino acid residues found within the opsin protein sequence. Minor differences in the amino acid sequences of the opsins are known to lead to large differences in the spectral peak of absorbance (i.e. the λ max value). In our prior studies, we developed a new approach that combined homology modeling and molecular dynamics simulations to gather structural information associated with chromophore conformation, then used it to generate statistical models for the accurate prediction of λ max values for photopigments derived from Rh1 and Rh2 amino acid sequences. In the present study, we test our novel approach to predict the λ max of phylogenetically distant Sws2 cone opsins. To build a model that can predict the λ max using our approach presented in our prior studies, we selected a spectrally-diverse set of 11 teleost Sws2 photopigments for which both amino acid sequence information and experimentally measured λ max values are known. The final first-order regression model, consisting of three terms associated with chromophore conformation, was sufficient to predict the λ max of Sws2 photopigments with high accuracy. This study further highlights the breadth of our approach in reliably predicting λ max values of Sws2 cone photopigments, evolutionary-more distant from template bovine RH1, and provided mechanistic insights into the role of known spectral tuning sites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Solvent organization in the ultrahigh-resolution crystal structure of crambin at room temperature

Ultrahigh-resolution structures provide unprecedented details about protein dynamics, hydrogen bonding and solvent networks. The reported 0.70 Å, room-temperature crystal structure of crambin is the highest-resolution ambient-temperature structure of a protein achieved to date. Sufficient data were collected to enable unrestrained refinement of the protein and associated solvent networks using SHELXL . Dynamic solvent networks resulting from alternative side-chain conformations and shifts in water positions are revealed, demonstrating that polypeptide flexibility and formation of clathrate-type structures at hydrophobic surfaces are the key features endowing crambin crystals with extraordinary diffraction power.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Membrane Proteins Have Distinct Fast Internal Motion and Residual Conformational Entropy

Abstract The internal motions of integral membrane proteins have largely eluded comprehensive experimental characterization. Here the fast side‐chain dynamics of the α‐helical sensory rhodopsin II and the β‐barrel outer membrane protein W have been investigated in lipid bilayers and detergent micelles by solution NMR relaxation techniques. Despite their differing topologies, both proteins have a similar distribution of methyl‐bearing side‐chain motion that is largely independent of membrane mimetic. The methyl‐bearing side chains of both proteins are, on average, more dynamic in the ps–ns timescale than any soluble protein characterized to date. Accordingly, both proteins retain an extraordinary residual conformational entropy in the folded state, which provides a counterbalance to the absence of the hydrophobic effect. Furthermore, the high conformational entropy could greatly influence the thermodynamics underlying membrane‐protein functions, including ligand binding, allostery, and signaling.

O'Brien, Evan S.↗
DOE OSTI2024

Design of intrinsically disordered protein variants with diverse structural properties

Intrinsically disordered proteins (IDPs) perform a broad range of functions in biology, suggesting that the ability to design IDPs could help expand the repertoire of proteins with novel functions. Computational design of IDPs with specific conformational properties has, however, been difficult because of their substantial dynamics and structural complexity. We describe a general algorithm for designing IDPs with specific structural properties. We demonstrate the power of the algorithm by generating variants of naturally occurring IDPs that differ in compaction, long-range contacts, and propensity to phase separate. We experimentally tested and validated our designs and analyzed the sequence features that determine conformations. We show how our results are captured by a machine learning model, enabling us to speed up the algorithm. Our work expands the toolbox for computational protein design and will facilitate the design of proteins whose functions exploit the many properties afforded by protein disorder.

Science & Technology - Other Topics↗
DOE OSTI2025

Machine Learning‐Augmented Molecular Dynamics Simulations (MD) Reveal Insights Into the Disconnect Between Affinity and Activation of ZTP Riboswitch Ligands

Abstract The challenge of targeting RNA with small molecules necessitates a better understanding of RNA–ligand interaction mechanisms. However, the dynamic nature of nucleic acids, their ligand‐induced stabilization, and how conformational changes influence gene expression pose significant difficulties for experimental investigation. This work employs a combination of computational and experimental methods to address these challenges. By integrating structure‐informed design, crystallography, and machine learning‐augmented all‐atom molecular dynamics simulations (MD), we synthesized, biophysically and biochemically characterized, and studied the dissociation of a library of small molecule activators of the 5‐aminoimidazole–4–carboxamide ribonucleotide triphosphate (ZTP) riboswitch, a ligand‐binding RNA motif that regulates bacterial gene expression. We uncovered key interaction mechanisms, revealing valuable insights into the role of ligand binding kinetics on riboswitch activation. Further, we established that ligand on‐rates determine activation potency as opposed to binding affinity and elucidated RNA structural differences, which provide mechanistic insights into the interplay of RNA structure on riboswitch activation.

Chemistry↗
DOE OSTI2023

Tuning Organic Semiconductor Packing and Morphology through Non-equilibrium Solution Processing (Final Report)

Organic semiconductors (OSCs) are a promising candidate to produce low-cost, large area, and flexible electronics. There has previously been successful development of organic field effect transistors (OFETs), photovoltaics (OPV), and bioelectronics using OSCs. Solution processing of these materials allows for the fabrication of large-area devices in the kinetic crystallization regime. The charge transport capabilities have been shown to depend on the morphology and molecular packing of the OSC within thin films. Our hypothesis is that solution processing conditions may significantly impact OSC morphology and as a result the charge transport. Therefore, this work has focused on obtaining a fundamental understanding of the morphology and molecular packing of OSCs for optimal device performance. Through our work, we have gained a better understanding of how these structural conditions were influenced by the solution processing conditions. Our studies have led to a more systematic understanding of the various parameters that impact OSC morphology. Solution processing leads to non-equilibrium films, thus allowing for the formation of diverse morphologies that are inaccessible by other fabrication methods. Previously, our group has focused on tuning the morphology of small-molecular OSCs [53, 55-57]. However, little work had been devoted to polymer OSCs. There was a lack of detailed studies characterizing the solution-state of polymer OSCs in terms of their conformation, degree of entanglement, polymer aggregation, and chain relaxation dynamics. The solution state properties are also likely to be influenced by the rigidity and molecular weight of the polymer OSCs investigated. Thus, we have focused our attention on gaining a better insight of polymer OSC films and solution processing methods. Through this proposal, our approach is to investigate the correlation between solution-state properties and the morphology of the resulting polymer OSC films. We worked three specific aims: investigate the effects of 1) polymer OSC solution-state properties on final film morphology; 2) molecular additives on the solution-state and final film properties; 3) controlled pre-aggregation in the solution-state in the final film morphology. More rigid and planar polymer backbones should promote interchain charge transport and more efficient interchain hopping between polymer chains. Through tailoring the polymer backbone, polymer sidechains, and molecular weight, we expected the altered solution-state properties to affect the final film morphology. In addition, reducing the entanglements and promoting chain alignment will likely prevent charge carrier trapping through conformation disorders. We thus studied different mechanisms, such as molecular additives, to reduce entanglements in solution. Devices fabricated from these solutions were expected to have improved charge transport abilities. In addition to tailoring the solution-state characteristics of the polymer OSCs, we investigated the effect of solution processing methods on the molecular packing and morphology of polymer OSC films. Two main solution processing techniques, e.g., spin-coating and solution shearing, were employed to fabricate OFETs. Spin coating was employed to prepare OSC films, which produce isotropic films. This technique creates several parameters to tune such as spin coating speed, acceleration, and time. Additionally, our research group developed the solution shearing method, which consists of the solution initially sandwiched between two plates. By sliding the top plate, the solution front is exposed and drying begins. This technique allows for the creation of aligned large crystalline domains. Solution shearing consists of different processing parameters which can affect the final morphology, such as shearing speed, substrate temperature and temperature gradient, distance between both plates, and the tilt angle of the top plate. Due to the complexity of the polymer systems investigated, numerous techniques were employed to characterize the polymer OSC solution state and films. All the materials were characterized using the DOE supported synchrotron X-ray scattering facilities at the Stanford Synchrotron Radiation Lightsource (SSRL). Using grazing incidence X-ray scattering (GIXS) and Near Edge X-ray Absorption Fine Structure (NEXAFS) techniques, the crystalline structure and molecular orientations of the thin films were measured. Optical absorption (UV-Vis) spectroscopy was employed to determine the aggregation state in solution and films of the polymer systems. Polarized UV-Vis also allowed for the determination of the relative degree of polymer chain alignment for solution sheared films. Additionally, various other techniques were used to investigate other properties within the film, such as atomic force microscopy (AFM) and solution rheology. Finally, the device performance is quantified through the fabrication and characterization of OFETs, which will highlight the effects of morphology on charge transport.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cryo-EM structure of transmembrane AAA+ protease FtsH in the ADP state

Abstract AAA+ proteases regulate numerous physiological and cellular processes through tightly regulated proteolytic cleavage of protein substrates driven by ATP hydrolysis. FtsH is the only known family of membrane-anchored AAA+ proteases essential for membrane protein quality control. Although a spiral staircase rotation mechanism for substrate translocation across the FtsH pore has been proposed, the detailed conformational changes among various states have not been clear due to absence of FtsH structures in these states. We report here the cryo-EM structure for Thermotoga maritima FtsH ( Tm FtsH) in a fully ADP-bound symmetric state. Comparisons of the ADP-state structure with its apo-state and a substrate-engaged yeast YME1 structure show conformational changes in the ATPase domains, rather than the protease domains. A reconstruction of the full-length Tm FtsH provides structural insights for the dynamic transmembrane and the periplasmic domains. Our structural analyses expand the understanding of conformational switches between different nucleotide states in ATP hydrolysis by FtsH.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Applying Crystallography and 19 F NMR to investigate dynamics and partner protein interactions in a long chain Flavodoxin

Flavodoxin (Fld) is a small FMN containing protein that is involved in single electron transfer. The long‐chain flavodoxin in Rhodopseudomonas palustris bacteria replaces ferredoxin as a low‐potential electron carrier when iron is scarce. Thus it is proposed to interact with the bifurcating electron transfer flavoprotein (ETF) that yields low‐potential electrons. A surface loop on Fld interacts with another of Fld's partner proteins, so we hypothesize that it also mediates Fld's interaction with ETF. To monitor interactions with ETF directly and investigate dynamics in this loop, we are using 19 F NMR in solution. 19 F is hyperresponsive to changes in its chemical environment with a chemical shift range of >300 ppm. To provide a static reference point and assess structural heterogeneity, we are also exploiting X‐ray crystallography. In this study we selectively fluorinated the five tyrosine residues in Fld. We obtained resonance assignments from 19 F spectra of Fld variants in which individual tyrosine residues have been replaced. We obtain well resolved signals for each residue, but the resonances' linewidths indicate dynamics that affects some resonances more than others. Y90 residue displays two resonances demonstrating two different conformations that interconvert slowly on an NMR time scale. Meanwhile the crystal structure solved at 2.1Å resolution reveals two molecules per asymmetric unit providing two perspectives on the details of the structure. The crystal symmetry is monoclinic in contrast to most of the other Flds, which are orthorhombic. Interestingly, the long loop bearing Y121 and Y123 is not well resolved in chain B of the crystal structure, and the NMR line of Y123 is exceptionally broad, both indicating that the loop is dynamic and capable of altering its conformation to accomodate binding to a partner protein. Future directions include monitoring the changes in the 19 F NMR of the Fld when titrated with the partner protein, temperature dependence of the NMR spectrum and relaxation studies to evaluate time scales of motions.

Khan, Sharique↗
DOE OSTI2022

Paired Simulations and Experimental Investigations into the Calcium-Dependent Conformation of Albumin

Serum albumin is the most abundant protein in blood plasma, and it is involved in multiple biological processes. Serum albumin has recently been adapted for improving biomaterial integration with bone tissue, and studies have shown the importance of this protein in bone repair and regeneration. However, the mechanism of action is not yet clear. In stark contrast, other studies have demonstrated that albumin blocks cell adhesion to surfaces, which is seen as a limitation to its bone healing role. These apparent contradictions suggest that the conformation of albumin facilitates its bioactivity, leading to enhanced bone repair. Serum albumin is known to play a major role in maintaining the calcium ion concentration in blood plasma. Due to the prevalence of calcium at bone repair and regeneration sites, it has been hypothesized that calcium binding to serum albumin triggers a conformational change, leading to bioactivity. In the current study, molecular modeling approaches including molecular docking, atomic molecular dynamics (MD) simulation, and coarse-grained MD simulation were used to test this hypothesis by investigating the conformational changes induced in bovine serum albumin by interaction with calcium ions. The computational results were qualitatively validated with experimental Fourier-transform infrared spectroscopy analysis. We find that free calcium ions in solution transiently bind with the three major loops in albumin, triggering a conformational change where N-terminal and C-terminal domains separate from each other in a partial unfolding process. The separation distance between these domains was found to correlate with the calcium ion concentration. The experimental data support the simulation results showing that albumin has enhanced conformational heterogeneity upon exposure to intermediate levels of calcium, without any significant secondary structure changes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Free Energy and Flexibility Analysis of Autoinhibited Human BRAF

The RAF serine/threonine protein kinases function as direct effectors of RAS in the intracellular transmission of extracellular growth signals, and they are key targets for drug discovery, given the high incidence of oncogenic mutations in RAF and other components of this signaling pathway. In its inactive state, RAF is held in an autoinhibited conformation in the cytosol through a combination of intramolecular interactions and binding to a regulatory 14−3−3 protein dimer. Activation of RAF is initiated by its interaction with membrane-localized GTP-bound RAS, which induces conformational changes that release RAF from its autoinhibited state. However, the molecular mechanisms governing RAF activation remain incomplete, largely due to the challenges in experimentally capturing the intermediate conformational states in this process. To address this gap, we developed a comprehensive all-atom model of BRAF based on existing cryo-EM structures. Using this model, we performed extensive molecular dynamics simulations to evaluate the stability and free energy landscape of autoinhibited BRAF in solution. Our analysis reveals conformational flexibility within the autoinhibited complex, suggesting that this dynamic behavior may play a role in facilitating BRAF activation upon engagement with the membrane-bound RAS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Probing DNA structural heterogeneity by identifying conformational subensembles of a bicovalently bound cyanine dye

DNA is a re-configurable, biological information-storage unit, and much remains to be learned about its heterogeneous structural dynamics. For example, while it is known that molecular dyes templated onto DNA exhibit increased photostability, the mechanism by which the structural dynamics of DNA affect the dye photophysics remains unknown. Here, we use femtosecond, two-dimensional electronic spectroscopy measurements of a cyanine dye, Cy5, to probe local conformations in samples of single-stranded DNA (ssDNA–Cy5), double-stranded DNA (dsDNA–Cy5), and Holliday junction DNA (HJ–DNA–Cy5). A line shape analysis of the 2D spectra reveals a strong excitation–emission correlation present in only the dsDNA–Cy5 complex, which is a signature of inhomogeneous broadening. Molecular dynamics simulations support the conclusion that this inhomogeneous broadening arises from a nearly degenerate conformer found only in the dsDNA–Cy5 complex. These insights will support future studies on DNA’s structural heterogeneity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Direct Modulators of K-Ras–Membrane Interactions

Protein–membrane interactions (PMIs) are ubiquitous in cellular signaling. Initial steps of signal transduction cascades often rely on transient and dynamic interactions with the inner plasma membrane leaflet to populate and regulate signaling hotspots. Methods to target and modulate these interactions could yield attractive tool compounds and drug candidates. Here, we demonstrate that the conjugation of a medium-chain lipid tail to the covalent K-Ras(G12C) binder MRTX849 at a solvent-exposed site enables such direct modulation of PMIs. The conjugated lipid tail interacts with the tethered membrane and changes the relative membrane orientation and conformation of K-Ras(G12C), as shown by molecular dynamics (MD) simulation-supported NMR studies. In cells, this PMI modulation restricts the lateral mobility of K-Ras(G12C) and disrupts nanoclusters. The described strategy could be broadly applicable to selectively modulate transient PMIs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Learning Latent Representations to Bridge Coarse-Grained and Atomistic Resolutions in Polymer Simulations

We present a machine-learning-based framework for learning reduced-order representations of polymer chain conformations across coarse-grained (CG) and united-atom (UA) fidelities. By employing linear singular value decomposition and nonlinear autoencoders, we compress high-dimensional polymer configurations into latent spaces with minimal loss of structural accuracy. Crucially, we demonstrate a near-perfect linear mapping between CG and UA latent spaces, enabling an efficient super-resolution back-mapping procedure that reconstructs high-fidelity UA configurations from CG simulations. While minor structural inaccuracies occur, they are effectively corrected through a brief molecular dynamics relaxation, forming a practical hybrid machine learning−physics scheme. This approach establishes the key structural prerequisites for accelerated polymer dynamics simulations: a compact and accurate latent encoding of polymer chain conformations and a validated multi-fidelity mapping that permits reconstruction of UA structures from CG configurations. The extension of this framework to explicit time evolution within the latent space, enabling dynamics to be propagated at CG fidelity and decoded to UA resolution only when required, represents a natural and well-motivated direction for future work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗