Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chain conformation”

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 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↗
NASA NTRS1992

Structure of a human monoclonal antibody Fab fragment against gp41 of human immunodeficiency virus type

The three-dimensional structure of a human monoclonal antibody (Fab), which binds specifically to a major epitope of the transmembrane protein gp41 of the human immunodeficiency virus type 1, has been determined by crystallographic methods to a resolution of 2.7 A. It has been previously determined that this antibody recognizes the epitope SGKLICTTAVPWNAS, belongs to the subclass IgG1 (kappa), and exhibits antibody-dependent cellular cytotoxicity. The quaternary structure of the Fab is in an extended conformation with an elbow bend angle between the constant and variable domains of 175 degrees. Structurally, four of the hypervariable loops can be classified according to previously recognized canonical structures. The third hypervariable loops of the heavy (H3) and light chain (L3) are structurally distinct. Hypervariable loop H3, residues 102H-109H, is unusually extended from the surface. The complementarity-determining region forms a hydrophobic binding pocket that is created primarily from hypervariable loops L3, H3, and H2.

He, X. M.↗
NASA NTRS1992

Structure of a human monoclonal antibody Fab fragment against gp41 of human immunodeficiency virus type 1

The three-dimensional structure of a human monoclonal antibody (Fab), which binds specifically to a major epitope of the transmembrane protein gp41 of the human immunodeficiency virus type 1, has been determined by crystallographic methods to a resolution of 2.7 A. It has been previously determined that this antibody recognizes the epitope SGKLICTTAVPWNAS, belongs to the subclass IgG1 (kappa), and exhibits antibody-dependent cellular cytotoxicity. The quaternary structure of the Fab is in an extended conformation with an elbow bend angle between the constant and variable domains of 175 deg. Structurally, four of the hypervariable loops can be classified according to previously recognized canonical structures. The third hypervariable loops of the heavy (H3) and light chain (L3) are structurally distinct. Hypervariable loop H3, residues 102H-109H, is unusually extended from the surface. The complementarity-determining region forms a hydrophobic binding pocket that is created primarily from hypervariable loops L3, H3, and H2.

He, Xiao M.↗
DOE OSTI2023

Molecular Cage Reports on Its Contents: Spectroscopic Signatures of Cryo-Cooled K + - and Ba 2+ -Benzocryptand Complexes

UV photofragment spectroscopy and IR–UV double resonance methods are used to determine the structure and spectroscopic responses of a three-dimensional [2.2.2]-benzocryptand cage to the incorporation of a single K + or Ba 2+ imbedded inside it (labeled as K + -BzCrypt, Ba 2+ -BzCrypt). We studied the isolated ion-cryptand complex under cryo-cooled conditions, brought into the gas phase by nano-electrospray ionization. Incorporation of a phenyl ring in place of the central ethyl group in one of the three N-CH 2 -CH 2 -O-CH 2 -CH 2 -O-CH 2 -CH 2 -N chains provides a UV chromophore whose S 0 –S 1 transition we probe. K + -BzCrypt and Ba 2+ -BzCrypt have their S 0 –S 1 origin transitions at 35,925 and 36,446 cm –1 , respectively, blue-shifted by 174 and 695 cm –1 from that of 1,2-dimethoxybenzene. These origins are used to excite a single conformation of each complex selectively and record their IR spectra using IR–UV dip spectroscopy. The alkyl CH stretch region (2800–3000 cm –1 ) is surprisingly sensitive to the presence and nature of the encapsulated ion. We carried out an exhaustive conformational search of cage conformations for K + -BzCrypt and Ba 2+ -BzCrypt, identifying two conformations (A and B) that lie below all others in energy. Here, we extend our local mode anharmonic model of the CH stretch region to these strongly bound ion-cage complexes to predict conformation-specific alkyl CH stretch spectra, obtaining quantitative agreement with experiment for conformer A, the gas-phase global minimum. The large electrostatic effect of the charge on the O- and N-lone pairs affects the local mode frequencies of the CH 2 groups adjacent to these atoms. The localized CH 2 scissors modes are pushed up in frequency by the adjacent O/N-atoms so that their overtones have little effect on the alkyl CH stretch region. However, the localized CH 2 wags are nearly degenerate and strongly coupled to one another, producing an array of delocalized wag normal modes, whose highest frequency members reach up above 1400 cm –1 . As such, their overtones mix significantly with the CH stretch modes, most notably involving the CH 2 symmetric stretch fundamentals of the central ethyl groups in the all-alkyl chains and the CH stretches adjacent to the N-atoms and antiperiplanar to the nitrogen lone pair.

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

Insights into the Assembly and Conformation of Nanoparticle Organic Hybrid Materials (NOHMs) in Solution with Varying Grafting Type

Recently, nanoparticle organic hybrid materials (NOHMs) have been considered promising electrolytes for energy storage applications due to their unique combination of properties, which include high thermal stability, negligible vapor pressure, and easy synthesis. However, the structural features of NOHMs in electrochemical systems are not well understood. Hence, we investigate the impact of grafting type (ionic vs covalent) on the structure and performance of NOHMs that consist of a silica core and grafted HPE polymer (Jeffamine M2070) in aqueous solution with and without the presence of the supporting electrolyte (0.1 M KHCO 3 ) using small-angle neutron scattering. Careful analyses of the results indicate that even at low grafting density (0.8 chains nm –2 ) and concentration (1–3 wt %), ionic NOHMs solutions contain a free polymer that may interact with the grafted polymer, altering the assembly of NOHMs in solution. Our investigation also elucidates that the covalent NOHM solutions consist of both aggregated and dispersed NOHMs in the solution. Moreover, the addition of the supporting electrolyte affects the assembly and structure of the polymer in the ionic NOHMs solution significantly, whereas only a slight change in the conformation of the grafted polymer is observed for covalent NOHMs. These conformational changes alter the performance of the NOHMs solution, resulting in a dramatic change in viscosity. Careful analysis shows that this significant alteration of solution viscosity in the ionic NOHMs can be traced back to the presence of the free polymer that interacts with the grafted layer and the attenuation of this interaction between the free and grafted polymers with the addition of salt. Furthermore, these studies elucidate the impact of grafting type on the overall structure, assembly, and performance of functionalized NOHMs in aqueous solution and explicitly correlate the structure and assembly of NOHMs to their transport properties and must be considered when designing functional systems that incorporate NOHMs.

25 ENERGY STORAGE↗
DOE OSTI2024

CryoEM structures of the human CLC-2 voltage-gated chloride channel reveal a ball-and-chain gating mechanism

CLC-2 is a voltage-gated chloride channel that contributes to electrical excitability and ion homeostasis in many different tissues. Among the nine mammalian CLC homologs, CLC-2 is uniquely activated by hyperpolarization, rather than depolarization, of the plasma membrane. The molecular basis for the divergence in polarity of voltage gating among closely related homologs has been a long-standing mystery, in part because few CLC channel structures are available. Here, we report cryoEM structures of human CLC-2 at 2.46 – 2.76 Å, in the presence and absence of the selective inhibitor AK-42. AK-42 binds within the extracellular entryway of the Cl – -permeation pathway, occupying a pocket previously proposed through computational docking studies. In the apo structure, we observed two distinct conformations involving rotation of one of the cytoplasmic C-terminal domains (CTDs). In the absence of CTD rotation, an intracellular N-terminal 15-residue hairpin peptide nestles against the TM domain to physically occlude the Cl – -permeation pathway. This peptide is highly conserved among species variants of CLC-2 but is not present in other CLC homologs. Previous studies suggested that the N-terminal domain of CLC-2 influences channel properties via a “ball-and-chain” gating mechanism, but conflicting data cast doubt on such a mechanism, and thus the structure of the N-terminal domain and its interaction with the channel has been uncertain. Through electrophysiological studies of an N-terminal deletion mutant lacking the 15-residue hairpin peptide, we support a model in which the N-terminal hairpin of CLC-2 stabilizes a closed state of the channel by blocking the cytoplasmic Cl – -permeation pathway.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Effective description of sub-maximal chaos: stringy effects for SYK scrambling

It has been proposed that the exponential decay and subsequent power law saturation of out-of-time-order correlation functions can be universally described by collective ‘scramblon’ modes. We develop this idea from a path integral perspective in several examples, thereby establishing a general formalism. After reformulating previous work on the Schwarzian theory and identity conformal blocks in two-dimensional CFTs relevant for systems in the infinite coupling limit with maximal quantum Lyapunov exponent, we focus on theories with sub-maximal chaos: we study the large-q limit of the SYK quantum dot and chain, both of which are amenable to analytical treatment at finite coupling. In both cases we identify the relevant scramblon modes, derive their effective action, and find bilocal vertex functions, thus constructing an effective description of chaos. The final results can be matched in detail to stringy corrections to the gravitational eikonal S-matrix in holographic CFTs, including a stringy Regge trajectory, bulk to boundary propagators, and multi-string effects that are unexplored holographically.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Dynamic Implications of Noncovalent Interactions in Amphiphilic Single-Chain Polymer Nanoparticles

Single-chain polymer nanoparticles (SCNPs) combine the chemical diversity of synthetic polymers with the intricate structure of biopolymers, generating versatile biomimetic materials. The mobility of polymer chain segments at length scales similar to secondary structural elements in proteins is critical to SCNP structure and thus function. However, the influence of noncovalent interactions used to form SCNPs (e.g., hydrogen-bonding and biomimetic secondary-like structure) on these conformational dynamics is challenging to quantitatively assess. To isolate the effects of noncovalent interactions on SCNP structure and conformational dynamics, we synthesized a series of amphiphilic copolymers containing dimethylacrylamide and monomers capable of forming these different interactions: (1) di(phenylalanine) acrylamide that forms intramolecular β-sheet-like cross-links, (2) phenylalanine acrylamide that forms hydrogen-bonds but lacks a defined local structure, and (3) benzyl acrylamide that has the lowest propensity for hydrogen-bonding. Each SCNP formed folded structures comparable to those of intrinsically disordered proteins, as observed by size exclusion chromatography and small angle neutron scattering. The dynamics of these polymers, as characterized by a combination of dynamic light scattering and neutron spin echo spectroscopy, was well described using the Zimm with internal friction (ZIF) model, highlighting the role of each noncovalent interaction to additively restrict the internal relaxations of SCNPs. These results demonstrate the utility of local scale interactions to control SCNP polymer dynamics, guiding the design of functional biomimetic materials with refined binding sites and tunable kinetics.

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

Variational quantum simulation of the critical Ising model with symmetry averaging

Here we investigate the use of deep multiscale entanglement renormalization ansatz (DMERA) circuits as a variational ansatz. We use the exactly solvable one-dimensional critical transverse-field Ising model as a test bed. Numerically exact simulation of the quantum circuit ansatz can in this case be carried out to hundreds of qubits by exploiting efficient classical algorithms for simulating matchgate circuits. We find that, for this system, the DMERA strongly outperforms a standard quantum approximate optimization algorithm (QAOA)–style ansatz, and that a major source of systematic error in correlation functions approximated using the DMERA is the breaking of the translational and Kramers-Wannier symmetries of the transverse-field Ising model. We are able to reduce this error by up to four orders of magnitude by symmetry averaging, without incurring additional cost in qubits or circuit depth. Here, we propose that this technique for mitigating systematic error could be applied to noisy intermediate-scale quantum (NISQ) simulations of physical systems with other symmetries.

1-dimensional spin chains↗
DOE OSTI2022

Eigenstate entanglement scaling for critical interacting spin chains

With increasing subsystem size and energy, bipartite entanglement entropies of energy eigenstates cross over from the groundstate scaling to a volume law. In previous work, we pointed out that, when strong or weak eigenstate thermalization (ETH) applies, the entanglement entropies of all or, respectively, almost all eigenstates follow a single crossover function. The crossover functions are determined by the subsystem entropy of thermal states and assume universal scaling forms in quantum-critical regimes. This was demonstrated by field-theoretical arguments and the analysis of large systems of non-interacting fermions and bosons. Here, we substantiate such scaling properties for integrable and non-integrable interacting spin-1/2 chains at criticality using exact diagonalization. In particular, we analyze XXZ and transverse-field Ising models with and without next-nearest-neighbor interactions. Indeed, the crossover of thermal subsystem entropies can be described by a universal scaling function following from conformal field theory. Furthermore, we analyze the validity of ETH for entanglement in these models. Even for the relatively small system sizes that can be simulated, the distributions of eigenstate entanglement entropies are sharply peaked around the subsystem entropies of the corresponding thermal ensembles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The impact of the electric field of metal ions on the vibrations and internal hydrogen bond strength in alkali metal ion di- and triglycine complexes

Using infrared predissociation spectroscopy of cryogenic ions, we revisit the vibrational spectra of alkali metal ion (Li+, Na+, K+) di- and triglycine complexes. We assign their most stable conformation, which involves metal ion coordination to all C=O groups and an internal NH⋯NH2 hydrogen bond in the peptide backbone. An analysis of the spectral shifts of the OH and C=O stretching vibrations across the different metal ions and peptide chain lengths shows that these are largely caused by the electric field of the metal ion, which varies in strength as a function of the square of the distance. The metal ion–peptide interaction also remotely modulates the strength of internal hydrogen bonding in the peptide backbone via the weakening of the amide C=O bond, resulting in a decrease in internal hydrogen bond strength from Li+ > Na+ > K+.

Chemistry↗
DOE OSTI2020

NMR Studies of Block Copolymer-Based Supramolecules in Solution

Hierarchical assemblies from block copolymer (BCP)-based supramolecules have shown immense potential as programmable materials owing to their versatility for incorporating functional molecules and provide access to arrays of hierarchical structures. However, there remains a knowledge gap on the formation of the supramolecule in solution. Here, we applied NMR techniques to investigate the solution-phase behavior of the most studied supramolecular systems, polystyrene-block-poly(4-vinylpyridine)(3-pentadecylphenol) (PS-b-P4VP(PDP) r ). The results show that the supramolecule likely adopts a coil-comb conformation, despite the small molecule's (PDP) rapid exchange between the bonded and free states. The exchange rate (>10 4 s -1 ) exceeds the NMR time scale at the frequency of interest. The supramolecules form under dilute conditions (~2 vol %) and are attributed to the enthalpic gain of the hydrogen bonding between the PDP and 4VP. As the solute concentration increases (>10 vol %), the supramolecule forms micelle-like aggregates with PDP accumulated within the comb-block's pervaded volume based on analysis of the apparent molecular weight, viscosity, and chain dynamics. Overall, this work sheds light on the long-standing question regarding the evolution of the constituents in the BCP-based supramolecule in solution and provides valuable guidance toward their solution-based processing and morphological control.

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

Temporal and spatial resolution of distal protein motions that activate hydrogen tunneling in soybean lipoxygenase

The enzyme soybean lipoxygenase (SLO) provides a prototype for deep tunneling mechanisms in hydrogen transfer catalysis. This work combines room temperature X-ray studies with extended hydrogen–deuterium exchange experiments to define a catalytically-linked, radiating cone of aliphatic side chains that connects an active site iron center of SLO to the protein–solvent interface. Employing eight variants of SLO that have been appended with a fluorescent probe at the identified surface loop, nanosecond fluorescence Stokes shifts have been measured. We report a remarkable identity of the energies of activation ( E a ) for the Stokes shifts decay rates and the millisecond C–H bond cleavage step that is restricted to side chain mutants within an identified thermal network. These findings implicate a direct coupling of distal protein motions surrounding the exposed fluorescent probe to active site motions controlling catalysis. While the role of dynamics in enzyme function has been predominantly attributed to a distributed protein conformational landscape, the presented data implicate a thermally initiated, cooperative protein reorganization that occurs on a timescale faster than nanosecond and represents the enthalpic barrier to the reaction of SLO.

Science & Technology - Other Topics↗
DOE OSTI2024

Secondary structure determines electron transport in peptides

Proteins play a key role in biological electron transport, but the structure–function relationships governing the electronic properties of peptides are not fully understood. Despite recent progress, understanding the link between peptide conformational flexibility, hierarchical structures, and electron transport pathways has been challenging. Here, we use single-molecule experiments, molecular dynamics (MD) simulations, nonequilibrium Green’s function-density functional theory (NEGF-DFT), and unsupervised machine learning to understand the role of secondary structure on electron transport in peptides. Our results reveal a two-state molecular conductance behavior for peptides across several different amino acid sequences. MD simulations and Gaussian mixture modeling are used to show that this two-state molecular conductance behavior arises due to the conformational flexibility of peptide backbones, with a high-conductance state arising due to a more defined secondary structure (beta turn or 3 10 helices) and a low-conductance state occurring for extended peptide structures. These results highlight the importance of helical conformations on electron transport in peptides. Conformer selection for the peptide structures is rationalized using principal component analysis of intramolecular hydrogen bonding distances along peptide backbones. Molecular conformations from MD simulations are used to model charge transport in NEGF-DFT calculations, and the results are in reasonable qualitative agreement with experiments. Projected density of states calculations and molecular orbital visualizations are further used to understand the role of amino acid side chains on transport. Overall, our results show that secondary structure plays a key role in electron transport in peptides, which provides broad avenues for understanding the electronic properties of proteins.

Science & Technology - Other Topics↗
DOE OSTI2025

Allosteric inhibition of the IZUMO1–JUNO fertilization complex by the naturally occurring antisperm antibody OBF13

Sperm IZUMO1 binds to egg JUNO, and this interaction is essential for mammalian fertilization. Isolated from a female mouse immunized with syngeneic sperm, the antisperm antibody OBF13 recognizes IZUMO1 and inhibits murine fertilization. How OBF13 interferes with sperm–egg interactions was unknown. Here, we present the X-ray crystal structure of IZUMO1 in complex with OBF13. OBF13 binds to the apex of the four-helix domain of IZUMO1, distant from the JUNO-binding site. Our crystal structure of OBF13-bound IZUMO1 resembles apo-IZUMO1 and differs from the structure of IZUMO1 in complex with JUNO. We identify that OBF13 carries a low level of somatic hypermutation, and through deep mutational scanning, we engineer an affinity-enhanced OBF13 variant. This OBF13 variant single-chain fragment variable decreases the apparent affinity of IZUMO1 for membrane-bound murine JUNO and blocks the binding of acrosome-reacted sperm to eggs, thereby preventing fertilization. We propose allostery between the OBF13 epitope and the JUNO-binding site. OBF13 inhibits a conformational change in IZUMO1, preventing fusion-competent sperm from adhering to murine eggs during fertilization. Surprisingly, murine IZUMO1 binds to hamster JUNO with an affinity ~20-fold higher than to murine JUNO. The decreased affinity caused by OBF13 of murine IZUMO1 for hamster JUNO is sufficient for murine sperm to bind to and fuse with hamster eggs. Our studies provide a structural and mechanistic framework for species-specific, allosteric inhibition of IZUMO1 by a naturally occurring antisperm antibody and offer insights into the development of immunocontraceptives.

Science & Technology - Other Topics↗
DOE OSTI2025

Fluorinated polythiophenes with ester side chains for boosting the V OC and efficiency in non-fullerene polymer solar cells

As the power conversion efficiency (PCE) of organic photovoltaics (OPVs) approximates the 19 % threshold, wide bandgap (WBG) polythiophenes (PTs) have gained increasing attention due to their simple chemical structures and ease of synthesis, making them promising candidates for large-scale production. However, the benchmark polymer poly(3-hexylthiophene) (P3HT) is limited by its high-lying HOMO energy level, which restricts the open-circuit voltage (V OC ) in solar cells. Here, in this study, we introduce a novel series of PT derivatives (PDC8-T, PDC8-3T, PDC16-3T, and PDC16-3T-2F) featuring ester side chains designed to fine-tune electronic properties through a streamlined three-step synthesis. Additionally, we incorporated a π-spacer to reduce steric hindrance and elongated alkyl side chains to improve solubility and processability. Compared to P3HT, these PT derivatives demonstrate a significant reduction in HOMO energy levels, lowering by approximately 0.3–0.4 eV. Among them, PDC16-3T-2F—with fluorine atom substitution—achieves the lowest HOMO energy level, induces a coplanar molecular conformation, and enhances polymer aggregation behavior. We evaluated these PT derivatives in inverted non-fullerene bulk-heterojunction (NFA BHJ) OPVs. The PDC8-3T device showed a relatively low PCE of 0.69 %, with a V OC of 0.76 V, a short-circuit current density (J SC ) of 3.32 mA/cm 2 , and a fill factor (FF) of 27.3 %. In contrast, the PDC16-3T-2F device achieved an impressive PCE of 7.21 %, with a V OC of 0.85 V, a J SC of 14.60 mA/cm 2 , and an FF of 58.4 %. This remarkable improvement is attributed to the fluorine substitution, which not only enhances molecular orientation but also downshifts the HOMO energy level and further boosts the V OC . Hence, these molecular design strategies have led to a fibrillar bicontinuous interpenetrating network with optimal nanoscale phase separation within the active layer.

14 SOLAR ENERGY↗
DOE OSTI2024

Correction to “Frank–Kasper Phases of Diblock Copolymer Melts Studied with the DPD Model: SCF Results”

In a recent study, (1) we performed the self-consistent field (SCF) calculations of a dissipative particle dynamics (DPD) model of conformationally asymmetric diblock copolymer A-B melts to study the stability of various Frank–Kasper (FK) phases formed by such A-B melts. While all the equations in ref (1) are correct, we did not calculate the dimensionless (mean-field) Helmholtz free energies per chain βf c = βu c,χ + βu c,κ – s c /k B and the dimensionless stresses dβf c /dθ of ordered phases accurately. For example, the dimensionless (nonbonded) internal energy per chain due to the system compressibility βu c,κ = (N/2κ)∑ q βû 0 (|q|)(ϕ̂ A (q) + ϕ̂ B (q) – δ q,0 )(ϕ̂ A (−q) + ϕ̂ B (−q) – δ– q,0 ) should be calculated in the reciprocal space as given (note that βû 0 (|q|) is obtained analytically, and that since ϕ A (r) is real 𝜙̂ A (−𝐪)= $\overline{𝜙̂ A}$(𝐪), i.e., the complex conjugate of ϕ̂ A (q)); this corresponds to numerically evaluating βu c,κ = (N/2κ)∫dr(ϕ A (r) + ϕ B (r) – 1)ψ(r) with ψ(r) ≡ (∫dr′βu 0 (|r–r′|)(ϕ A (r′) + ϕ B (r′) – 1)/V in the real space via the (composite) trapezoidal rule. While the trapezoidal rule is less accurate than the Romberg integration (2) in general, the opposite occurs when the integration domain is periodic; this was pointed out by Matsen in ref (3) and explained in detail in refs (4and5). With our ignorance of this exceptional convergence property of the trapezoidal rule for periodic functions, we calculated all convolutions having the form of ψ(r) in the reciprocal space via the fast Fourier transforms (6) but those having the form of βu c,κ in the real space via Romberg integration; the latter leads to inaccurate results and unfortunately incorrect phase diagrams reported in ref (1), which are corrected here.

He, Juntong↗
DOE OSTI2023

Atomistic Simulations of Polydisperse Lignin Melts Using Simple Polydisperse Residue Input Generator

Understanding the physics of lignin will help rationalize its function in plant cell walls as well as aiding practical applications such as deriving biofuels and bioproducts. Here, in this work, we present SPRIG (Simple Polydisperse Residue Input Generator), a program for generating atomic-detail models of random polydisperse lignin copolymer melts i.e., the state most commonly found in nature. Using these models, we use all-atom molecular dynamics (MD) simulations to investigate the conformational and dynamic properties of polydisperse melts representative of switchgrass (Panicum virgatum L.) lignin. Polydispersity, branching and monolignol sequence are found to not affect the calculated glass transition temperature, T g . The Flory–Huggins scaling parameter for the segmental radius of gyration is 0.42 ± 0.02, indicating that the chains exhibit statistics that lie between a globular chain and an ideal Gaussian chain. Below T g the atomic mean squared displacements are independent of molecular weight. In contrast, above T g , they decrease with increasing molecular weight. Therefore, a monodisperse lignin melt is a good approximation to this polydisperse lignin when only static properties are probed, whereas the molecular weight distribution needs to be considered while analyzing lignin dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗