Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Protein interaction network”

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.

179 records · Page 10

Driving Electron Transfer in Photosystem I Using Far-Red Light: Overall Perspectives

Photosystem I (PSI) is a photosynthetic protein–pigment complex that, upon photoexcitation, transfers electrons to ferredoxin, facilitating the production of NADPH. Isolated PSI reaction centers (RCs) have also been used in hybrid systems to reduce protons and produce ‘biohydrogen’. This review article examines how various cyanobacteria with similar photosynthetic machinery utilize different wavelengths of light to execute photosynthetic electron transport through PSI. Key factors, such as, the structure of the electron transfer cofactors, the protein environment surrounding the primary donor pigments and hydrogen-bonding interactions with the surrounding protein matrix are analyzed to understand their roles in maintaining efficient electron transfer when it is driven using photons of different energies. We compare PSI complexes with known atomic structures from four species of cyanobacteria, Thermosynechococcus elongatus, Acaryochloris marina, Halomicronema hongdechloris, and Fischerella thermalis. T. elongatus is typical of most oxygenic photosynthetic organisms in that it requires visible light and uses only chlorophyll a (Chl a ) in PSI. In contrast, H. hongdechloris and F. thermalis are photoacclimating species capable of producing Chl f and Chl d that use red light when little visible light is available. A. marina , on the other hand, is adapted to red light conditions and consistently utilizes Chl d as its primary photosynthetic pigment, maintaining a stable pigment composition. Here, we explore the structural and functional differences between the PSI RCs of these organisms and the impact of these differences on electron transport. The structural differences in the cofactors influence both the absorption wavelengths of the cofactors and the energy levels of the intermediate states of electron transfer. An analysis of the surrounding protein shows how it has been adapted and underscores the interplay between the pigment structure, protein environment, and hydrogen bonding networks in tuning the efficiency and adaptability of photosynthetic mechanisms across different species of cyanobacteria.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Heterogenous catalysis for oxygen tolerant photoredox atom transfer radical polymerization and small-molecule dehalogenation

Heterogeneous photocatalysts (PCs) have garnered attention for their sustainability and cost-effectiveness. Despite the existence of various types of these PCs, their synthesis often involves complex, multi-step procedures and laborious purification. Herein, we propose a simple method for attaching small-molecule photocatalytic species onto crosslinked 3-D polymer networks as insoluble scaffolds to create robust heterogeneous PCs. The highly swellable poly(ethylene glycol)-based ChemMatrix (CM) resin, known for its amphiphilic properties and high functional group loading, facilitated the covalent immobilization of the photoredox dye Eosin Y (EY), but also streamlined functionalization with Ir( III ) complexes. The resulting heterogeneous CM-EY demonstrated efficient photocatalytic performance in open-to-air dual photoredox catalysis of atom transfer radical polymerization (photo-ATRP) under green light. This was confirmed by the well-controlled synthesis of polymers with molecular masses ranging from 20 kDa to 300 kDa and low dispersities. Furthermore, CM-EY exhibited excellent photostability and recyclability over multiple cycles of ATRP. The heterogeneous catalysis of photo-ATRP provided high temporal control and enabled benign conditions for synthesizing protein-polymer hybrids (PPH). When combined with the initiator-modified CM (CM-BIB), CM-EY facilitated the solid-phase synthesis of homopolymers and block copolymers with recyclable performance. However, the coordinatively bound Ir@CM showed decreased catalytic activity and efficiency toward photoredox dehalogenation due to the leaching of active species during recycling. This study highlights the advantages of the covalent linking of catalysts to solid supports over non-covalent interactions, underscoring the potential of functionalized polymer resin as a promising scaffold. Such an approach offers customization and tunability, presenting opportunities for innovation in green chemistry.

Kapil, Kriti↗

A global lipid map defines a network essential for Zika virus replication

Zika virus (ZIKV), an emerging mosquito-borne pathogen of global concern, remodels intra-cellular lipid membranes to form replication factories. How ZIKV regulates lipid networks to allow infection, and consequences for disease is poorly understood. Here, we performed comprehensive lipidomics to create a lipid network map during ZIKV infection. We found that ZIKV significantly alters host lipid composition with striking changes seen within sub-classes of sphingolipids. We also demonstrated that ectopic expression of NS4B, one of the non-structural proteins of ZIKV with known roles in ZIKV replication resulted in similar changes on cellular lipidome. Ceramide, a key bioactive molecule with functions in signaling and apoptosis, redistributes to ZIKV replication factories. Increase of ceramide through knockout of sphingomyelin synthase or via sphingomyelinase-catalyzed degradation of sphingomyelin to ceramide sensitizes cells to ZIKV infection. Remarkably, we observed an enrichment of pro-apoptotic ceramide sub-species in infected cells, providing a link between ZIKV biogenesis and apoptosis. Thus, we identify a sphingolipid metabolic network with a critical role in ZIKV replication and show that ceramide flux is a key mediator of ZIKV infection.

60 APPLIED LIFE SCIENCES↗

Charged Particle Tracking via Edge-Classifying Interaction Networks

Recent work has demonstrated that geometric deep learning methods such as graph neural networks (GNNs) are well suited to address a variety of reconstruction problems in high-energy particle physics. In particular, particle tracking data are naturally represented as a graph by identifying silicon tracker hits as nodes and particle trajectories as edges, given a set of hypothesized edges, edge-classifying GNNs identify those corresponding to real particle trajectories. In this work, we adapt the physics-motivated interaction network (IN) GNN toward the problem of particle tracking in pileup conditions similar to those expected at the high-luminosity Large Hadron Collider. Assuming idealized hit filtering at various particle momenta thresholds, we demonstrate the IN’s excellent edge-classification accuracy and tracking efficiency through a suite of measurements at each stage of GNN-based tracking: graph construction, edge classification, and track building. The proposed IN architecture is substantially smaller than previously studied GNN tracking architectures; this is particularly promising as a reduction in size is critical for enabling GNN-based tracking in constrained computing environments. Furthermore, the IN may be represented as either a set of explicit matrix operations or a message passing GNN. Efforts are underway to accelerate each representation via heterogeneous computing resources towards both high-level and low-latency triggering applications.

accelerator physics↗

Bond order predictions using deep neural networks

Machine learning is an extremely powerful tool for the modern theoretical chemist since it provides a method for bypassing costly algorithms for solving the Schrödinger equation. Already, it has proven able to infer molecular and atomic properties such as charges, enthalpies, dipoles, excited state energies, and others. Most of these machine learning algorithms proceed by inferring properties of individual atoms, even breaking down total molecular energy into individual atomic contributions. In this paper, we introduce a modified version of the Hierarchically Interacting Particle Neural Network (HIP-NN) capable of making predictions on the bonds between atoms rather than on the atoms themselves. We train the modified HIP-NN to infer bond orders for a large number of small organic molecules as computed via the Natural Bond Orbital package. We demonstrate that the trained model is extensible to molecules much larger than those in the training set by studying its performance on the COMP6 dataset. This method has applications in cheminformatics and force field parameterization and opens a promising future for machine learning models to predict other quantities that are defined between atoms such as density matrix elements, Hamiltonian parameters, and molecular reactivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Identification of genetic interactions with priB links the PriA/PriB DNA replication restart pathway to double-strand DNA break repair in Escherichia coli

Collisions between DNA replication complexes (replisomes) and impediments such as damaged DNA or proteins tightly bound to the chromosome lead to premature dissociation of replisomes at least once per cell cycle in Escherichia coli. Left unrepaired, these events produce incompletely replicated chromosomes that cannot be properly partitioned into daughter cells. DNA replication restart, the process that reloads replisomes at prematurely terminated sites, is therefore essential in E. coli and other bacteria. Three replication restart pathways have been identified in E. coli: PriA/PriB, PriA/PriC, and PriC/Rep. A limited number of genetic interactions between replication restart and other genome maintenance pathways have been defined, but a systematic study placing replication restart reactions in a broader cellular context has not been performed. We have utilized transposon-insertion sequencing to identify new genetic interactions between DNA replication restart pathways and other cellular systems. Known genetic interactors with the priB replication restart gene (uniquely involved in the PriA/PriB pathway) were confirmed and several novel priB interactions were discovered. Targeted genetic and imaging-based experiments with priB and its genetic partners revealed significant double-strand DNA break accumulation in strains with mutations in dam, rep, rdgC, lexA, or polA. Modulating the activity of the RecA recombinase partially suppressed the detrimental effects of rdgC or lexA mutations in ΔpriB cells. Taken together, our results highlight roles for several genes in double-strand DNA break homeostasis and define a genetic network that facilitates DNA repair/processing upstream of PriA/PriB-mediated DNA replication restart in E. coli.

59 BASIC BIOLOGICAL SCIENCES↗

The SEEL motif and members of the MYB-related REVEILLE transcription factor family are important for the expression of LORELEI in the synergid cells of the Arabidopsis female gametophyte

Synergid cells in the micropylar end of the female gametophyte are required for critical cell–cell signaling interactions between the pollen tube and the ovule that precede double fertilization and seed formation in flowering plants. LORELEI (LRE) encodes a putative GPI-anchored protein that is expressed primarily in the synergid cells, and together with FERONIA, a receptor-like kinase, it controls pollen tube reception by the receptive synergid cell. Still, how LRE expression is controlled in synergid cells remains poorly characterized. We identified candidate cis-regulatory elements enriched in LRE and other synergid cell-expressed genes. One of the candidate motifs (‘TAATATCT’) in the LRE promoter was an uncharacterized variant of the Evening Element motif that we named as the Short Evening Element-like (SEEL) motif. Deletion or point mutations in the SEEL motif of the LRE promoter resulted in decreased reporter expression in synergid cells, demonstrating that the SEEL motif is important for expression of LRE in synergid cells. Additionally, we found that LRE expression is decreased in the loss of function mutants of REVEILLE (RVE) transcription factors, which are clock genes known to bind the SEEL and other closely related motifs. Here, we propose that RVE transcription factors regulate LRE expression in synergid cells by binding to the SEEL motif in the LRE promoter. Identification of cis-regulatory elements and transcription factors involved in the expression of LRE will serve as a foundation to characterize the gene regulatory networks in synergid cells.

59 BASIC BIOLOGICAL SCIENCES↗

Structure and Sulfur: Tuning the Viscoelastic and Surface Properties of Natural Keratin Fibers

Natural keratin fibers, such as wool, possess a complex hierarchical structure that governs their mechanical properties and surface energy. However, the extent to which these characteristics are influenced by combined contributions of structural variations (e.g., fiber diameter, intermediate filament (IF) packing) and chemical composition (e.g., disulfide bond density) remains poorly understood. In this study, we investigate wool fibers from five sheep breeds (Merino, Polwarth, Cheviot, Eider, and Devon) to elucidate how these factors influence viscoelasticity and surface interactions. Using a multimodal approach integrating interfacial and bulk characterization methods, including inverse gas chromatography (IGC), atomic force microscopy-infrared spectroscopy (AFM-IR), X-ray photoelectron spectroscopy (XPS), uniaxial tensile testing, and synchrotron small-angle X-ray scattering (SAXS), we show that the nanometer-thick 18-methyleicosanoic acid (18-MEA) layer is consistently present across all wool types and plays a key role in governing hydrophobicity and surface heterogeneity. A controlled isothermal treatment at 200 °C, designed to cleave disulfide bonds, results in a nearly 40% reduction in specific surface area across all fiber types, accompanied by a significant decrease in tensile strength and 80% reduction in elongation at break for Merino and Devon wool, but limited influence on the mechanical properties of Eider fibers. Furthermore, rate-dependent tensile testing within the elastic regime reveals distinct viscoelastic responses among the fiber types, suggesting that the sulfur-rich protein matrix surrounding IFs and its structure contribute actively to stress partitioning. Altogether, when combined with conclusions from SAXS measurements of IF spacing, our work offers compelling insights into the role of the keratin-associated protein (KAP) matrix in shaping wool fiber mechanics. Differences in mechanical behavior among wool types, despite similar IF spacing or sulfur content, highlight the importance of matrix composition and cross-linking density, suggesting that the molecular architecture of the KAP network may be a dominant factor in determining fiber performance.

X-ray scattering↗

Multifunctional Silk Fibroin Hydrogels with Strong Adhesion for Tissue Sealing and Wearable Electronic Sensors

Multifunctional hydrogels with excellent adhesion, biodegradability, and conductivity are essential for overcoming the obstacles of postoperative secondary injury, flexible sensing instability, and so on. Herein, we develop a multifunctional silk fibroin (SF) hydrogel modified with poly(acrylic acid). Owing to the stable chemical cross-linking network and the abundant carboxylic acid groups of the SF network, the SF hydrogel exhibits a high tensile strength of 74.34 kPa due to sufficient cohesion and interfacial interactions. Additionally, the tensile strain reaches a maximum of 414.6%, the compressive strength is 0.9 MPa, and the shear adhesive strength for pig skin tissues is as high as 64 kPa. Compared with most hydrogels, our multifunctional SF hydrogel with a low swelling ratio provides excellent adhesion, biodegradation, and conductivity, which shows advantages in terms of invasive tissue sealing. Furthermore, the use of self-adhesive SF hydrogels as conductive hydrogels in flexible sensors also benefits the collection of physiological electricity and human motion signals in the field of wearable and implantable electronic devices.

adhesive↗

Preclinical Risk Evaluation of Normal Tissue Injury With Novel Radiosensitizers

Genotoxic damage induced by radiation triggers a highly coordinated DNA damage response, and molecular inhibitors of key nodes within this complex response network can profoundly enhance the antitumor efficacy of radiation. This is especially true for drugs targeting the catalytic subunit of DNA-dependent protein kinase, which is a core component of the nonhomologous end-joining DNA repair pathway, and ataxia telangiectasia mutated, which coordinates cell cycle arrest, apoptosis, and DNA repair functionalities after radiation exposure. Unlike the more modest in vitro radiosensitizing effects seen with classic sensitizing agents such as cisplatin, 5-fluorouracil, or taxanes, DNA-dependent protein kinase or ataxia telangiectasia mutated inhibitors provide much more robust sensitizing effects in vitro, as might be anticipated from targeting these key DNA repair modulators. However, patients with homozygous inactivating mutations of ataxia telangiectasia mutated or mice with homozygous defects in DNA-dependent protein kinase (severe combined immunodeficiency) have profoundly enhanced acute normal tissue radiation reactions. Therefore, there is significant potential that the combination of small molecule inhibitors of these kinases with radiation could cause similar dose-limiting acute normal tissue toxicities. Similarly, although less understood, inhibition of these DNA repair response pathways could markedly increase the risk of late radiation toxicities. Because these potent radiosensitizers could be highly useful to improve local control of otherwise radiation-resistant tumors, understanding the potential for elevated risks of radiation injury is essential for optimizing therapeutic ratio and developing safe and informative clinical trials. In this review, we will discuss 2 straightforward models to assess the potential for enhanced mucosal toxicity in the oral cavity and small intestine established in our laboratories. We also will discuss similar strategies for evaluating potential drug–radiation interactions with regard to increased risks of debilitating late effects.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Integrated lipidomic and proteomic profiling reveals metabolic network disruption by SARS-CoV-2 variants

The rapid evolution of SARS-CoV-2 has produced myriad viral strains with increasing transmissibility and capacity for immune evasion. While effective vaccination campaigns have reduced the fatalities associated with SARS-CoV-2, infections continue, and a detailed understanding of how this virus manipulates host biochemical pathways remains elusive. We asked both whether the patterns of host lipid rewiring remained consistent across variants and whether the changes in the abundance of lipid classes are related to changes in the expression of the enzymes involved in their biosynthesis. We compared global nontargeted lipidomics on A549-ACE2 cells infected with the delta variant (B.1.617.2), or the omicron (B.1.1.529) variant to our previous results of global nontargeted lipidomics on A549-ACE2 cells infected with the original WA1 strain and further performed quantitative proteomics to assess changes in the host proteome. We found that metabolic rewiring, both on the lipid and the enzymatic level, is remarkably consistent across all three variants. We further mapped changes in the expression of host metabolic enzymes, linking enzyme expression to alterations in the abundance of specific lipids during infection. This analysis identified key proteins related to virus-mediated changes in lipid abundance, including fatty acid synthase (FASN), lysosomal acid lipase (LIPA), and ORMDL, a regulator of sphingolipid biosynthesis. These integrated lipidomic and proteomic experiments shed light on the importance of the complex network of host metabolism networks that support SARS-CoV-2 infection and suggest that lipid metabolism may be a promising avenue for uncovering conserved therapeutic targets.

SARS-CoV-2↗

HDBind: encoding of molecular structure with hyperdimensional binary representations

Traditional methods for identifying “hit” molecules from a large collection of potential drug-like candidates rely on biophysical theory to compute approximations to the Gibbs free energy of the binding interaction between the drug and its protein target. These approaches have a significant limitation in that they require exceptional computing capabilities for even relatively small collections of molecules. Increasingly large and complex state-of-the-art deep learning approaches have gained popularity with the promise to improve the productivity of drug design, notorious for its numerous failures. However, as deep learning models increase in their size and complexity, their acceleration at the hardware level becomes more challenging. Hyperdimensional Computing (HDC) has recently gained attention in the computer hardware community due to its algorithmic simplicity relative to deep learning approaches. The HDC learning paradigm, which represents data with high-dimension binary vectors, allows the use of low-precision binary vector arithmetic to create models of the data that can be learned without the need for the gradient-based optimization required in many conventional machine learning and deep learning methods. This algorithmic simplicity allows for acceleration in hardware that has been previously demonstrated in a range of application areas (computer vision, bioinformatics, mass spectrometery, remote sensing, edge devices, etc.). To the best of our knowledge, our work is the first to consider HDC for the task of fast and efficient screening of modern drug-like compound libraries. We also propose the first HDC graph-based encoding methods for molecular data, demonstrating consistent and substantial improvement over previous work. We compare our approaches to alternative approaches on the well-studied MoleculeNet dataset and the recently proposed LIT-PCBA dataset derived from high quality PubChem assays. We demonstrate our methods on multiple target hardware platforms, including Graphics Processing Units (GPUs) and Field Programmable Gate Arrays (FPGAs), showing at least an order of magnitude improvement in energy efficiency versus even our smallest neural network baseline model with a single hidden layer. Our work thus motivates further investigation into molecular representation learning to develop ultra-efficient pre-screening tools. We make our code publicly available at https://github.com/LLNL/hdbind.

59 BASIC BIOLOGICAL SCIENCES↗

Cryo-EM Visualization of Intermolecular π-Electron Interactions within π-Conjugated Peptidic Supramolecular Polymers

The self-assembly of “π-peptides” – molecules with π-electron cores substituted with two or more oligopeptide chains – brings organic electronic function into biologically relevant nanomaterials. π-Peptides assemble into fibrillar nanomaterials as driven by enthalpic peptide-based hydrogen bonding networks and pi-core-based quadrupolar interactions. A large body of spectroscopic, morphological and computational studies informs on the nature of the self-assembly process and the resulting nanostructures, but detailed structural information has remained elusive. Here, inspired by the recent use of cryogenic electron microscopy (cryo-EM) to provide high-resolution structures for synthetic peptide nanomaterials, we present here the use of cryo-EM to offer ca. 3 Å resolution of π-peptide nanomaterial assemblies, visualizing for the first time the nature of the intermolecular π-core electronic interactions responsible for energy transport through these supramolecular materials.

Group theory↗

Single-cell proteomics of Arabidopsis leaf mesophyll reveals dynamic protein responses to water-deficit stress

Background The application of single-cell omics tools to biological systems can provide unique insights into diverse cellular populations and their heterogeneous responses to internal and external perturbations. Thus far, most single-cell studies in plant systems have been limited to RNA-sequencing approaches, which only provide indirect readouts of cellular functions. Results Here, we present a single-cell proteomics workflow for plant cells that integrates tape-sandwich protoplasting, piezoelectric cell sorting, nanoPOTS sample preparation, and ion mobility-based MS data acquisition method for label-free single-cell proteomics analysis of Arabidopsis leaf mesophyll cells. From a single leaf protoplast, over 3,000 proteins were quantified with high precision. The workflow is demonstrated to identify stress associated changes in protein abundance by analyzing 117 protoplasts from well-watered and water-deficit stressed plants. Additionally, we describe a new approach for constructing covarying protein networks at the single-cell level and demonstrate how single-cell protein covariation analysis can reveal previously unrecognized protein functions while also capturing stress-induced changes in protein–protein dynamics. Conclusions The label-free scProteomic approach presented here represents a significant advance through the demonstration of a facile protoplast isolation method combined with deep and precise proteomic coverage of Arabidopsis leaf mesophyll cell types. We believe this study will serve as an informative reference to future plant scProteomic investigations.

Arabidopsis↗

Hydrogen Bonding and Its Effect on the Orientational Dynamics of Water Molecules inside Polyelectrolyte Brush-Induced Soft and Active Nanoconfinement

Despite being the most ubiquitous compound on Earth, the fundamental properties of water are not fully understood, especially in nanoconfinement. Densely grafted polyelectrolyte (PE) molecules attain the configuration of a “brush”: these PE brushes, due to their ability to form hydrogen bonds (HBs) with water via the PE functional groups, act as a source of soft and active nanoconfinement for the brush-trapped water molecules. In this paper, we study the effects of PE brush-induced confinement on the structure, dynamics and energetics of the water-water and water-PE HBs. Our results indicate a significant weakening of the HBs from bulk to sparsely grafted to densely grafted brushes. i.e., by increasing the degree of brush-induced nanoconfinement. We explain that this weakening of water-water HBs is caused by the disruption of the extended network of water molecules within the brush-induced nanoconfinement. This is confirmed by performing a ring structure analysis of the water molecules, which yields a reduction in the average ring size at higher degrees of brush-induced nanoconfinement (i.e., at higher brush grafting densities). Furthermore, we investigate the role of HB on the orientational dynamics of the water molecules. Here, we observe that the rotational motion of the water molecules becomes sluggish inside the PE brushes. Recent findings have indicated that the water and counterions trapped in brush-induced nanoconfinement demonstrate structures (in combination with the PE functional group) analogous to that in “water-in-salt” electrolytes that have seen extensive recent uses for Li-ion battery applications. However, the rotational dynamics of water molecules inside the brush-induced nanoconfinement is found to be distinctly different from that of conventional “water-in-salt” electrolytes in absence of any confinement; therefore, the present study will provide the necessary platform towards conceptualizing polymer-based nanoconfinement for battery applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thiol post-translational modifications modulate allosteric regulation of the OpcA–G6PDH complex through conformational gate control

In cyanobacteria, the redox-sensitive protein OpcA acts as a metabolic switch for G6PDH, enabling rapid adjustment of reducing power generation from glycogen catabolism and thereby precisely regulating carbon flux between anabolic and catabolic pathways. Although redox-sensitive cysteines in OpcA are known to regulate G6PDH, the mechanisms by which redox post-translational modifications (PTMs) on OpcA control G6PDH structure and activity remain unclear. Here, we combine computational modeling with experimental redox proteomics in Synechococcus elongatus PCC 7942 to dissect this mechanism. Experimentally, redox proteome analysis revealed differential redox PTM patterns, particularly on cysteines within the G6PDH-binding site of OpcA. These environmentally sensitive PTM changes at the interface suggest that thiol modifications in this region form a key regulatory node. More broadly, redox proteomics identified site-specific cysteine modifications under light/dark transitions and circadian cycling, linking distinct redox regimes to discrete PTM states. We employed PTM-Psi simulations to show that thiol PTMs near the OpcA–G6PDH interface are critical for allosteric regulation of G6PDH. The thiol PTMs on OpcA affect a putative gate region in G6PDH for substrate ingress and product egress as well as key hydrogen-bond networks within the active site. We infer that PTMs on OpcA tune the conformational landscapes of individual G6PDH subunits toward functionally relevant configurations according to environmental gradients, biasing the enzyme toward catalytically favorable states. Together, our results reveal a molecular mechanism in which thiol PTMs on OpcA modulate G6PDH structure and function through PTM-induced reorganization of conformational dynamics and allosteric communication. These findings demonstrate that PTM-level regulation provides a critical control layer from genotypes to phenotypes that enables cyanobacteria to rapidly adapt to environmental fluctuations through precise metabolic fine-tuning.

Allosteric regulation↗

Modeling of Peptides with Classical and Novel Machine Learning Force Fields: A Comparison

The replacement of classical force fields (FFs) with novel neural-network-based frameworks is an emergent topic in molecular dynamics (MD) simulations. In contrast to classical FFs, which have proven their capability to provide insights into complex soft matter systems at an atomistic resolution, the machine learning (ML) potentials have yet to demonstrate their applicability for soft materials. However, the underlying philosophy, which is learning the energy of an atom in its surrounding chemical environment, makes this approach a promising tool. In particular for the exploration of novel chemical compounds, which have not been considered in the original parametrization of classical FFs. In this article, we study the performance of the ANI-2x ML model and compare the results with those of two classical FFs, namely, CHARMM27 and the GROMOS96 43a1 FF. We explore the performance of these FFs for bulk water and two model peptides, trialanine and a 9-mer of the α-aminoisobutyric acid, in vacuum and water. The results for water describe a highly ordered water structure, with a structure similar to those using ab initio molecular dynamics simulations. The energy landscape of the peptides described by Ramachandran maps show secondary structure basins similar to those of the classical FFs but differ in the position and relative stability of the basins. Details of the sampled structures show a divergent performance of the different models, which can be related either to the short-ranged nature of the ML potentials or to shortcomings of the underlying data set used for training. These findings highlight the current state of the applicability of ANI-2x ML potential for MD simulations of soft matter systems. Simultaneously, they provide insights for future improvements of current ML potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗