Search NASA⌕ Search

SEARCH · Search NASA

Results for “MD simulation”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Performance Advantaged ..beta..-Ketone Containing Polymers - Benefits in Manufacturing, Performance, and End-of-Life

Biomass derived chemicals can possess extended functionality, often in the form of heteroatoms like oxygen and nitrogen, that can enable performance advantaged properties in emergent materials. These performance advantages can manifest in multiple ways such as reduced manufacturing intensity, better performance, or enhanced end-of-life options. One illustrative example of this concept is ..beta..-ketoadipic acid (..beta..KA), a C6 dicarboxylic acid with a ketone in its backbone. ..beta..KA can be obtained via biological cultivation from aromatic compounds, sugars, and even waste plastics. When ..beta..KA is implemented into polymers, namely nylons, in place of adipic acid the overall material properties are increased. Specifically, the glass transition temperature (which is the measurement at which temperature a plastic softens) is increased by 69 degrees C and the water permeability is reduced by 20%. Molecular dynamic (MD) simulations reveal that the ..beta..-ketone in ..beta..KA leads to restricted movement across multiple backbone carbons when compared to adipic acid. Furthermore, MD simulations reveal that the ..beta..-ketone, when compared to no-ketone or an a-ketone, can lead to an increase in hydrogen bonding between neighboring chains. Both phenomena may help explain the molecular origin of the increases in polymer performance. When implemented into polyesters, namely PET, ..beta..KA can maintain the requisite performance while enabling enhanced degradation. Aside from performance advantages in thermomechanical performance, the production of ..beta..KA also is advantaged in its manufacturing. Although estimates for the minimum selling price (MSP) of ..beta..KA exceed adipic acid the manufacture of ..beta..KAwould reduce supply chain energy and GHG emissions by >50% and >30% respectively. Further TEA analysis reveals that the MSP of ..beta..KA is primarily driven by the feedstock cost while the capital expenditures are driven largely by fermentation costs. Overall, the work presented within will be illustrative of how biomass derived molecules, namely ..beta..KA, can manifest performance advantages across multiple benchmarks and may enable a path to market for biomass-derived materials.

ADVANCED PROPULSION SYSTEMS,BIOMASS FUELS↗
DOE OSTI2022

..beta..-Ketoadipic Acid Production in P. putida for Performance-Advantaged Nylons

Biomass derived chemicals can possess extended functionality, often in the form of heteroatoms like oxygen and nitrogen, that can enable performance advantaged properties in emergent materials. One illustrative example of this concept is ..beta..-ketoadipic acid (..beta..KA), a C6 dicarboxylic acid with a ketone in its backbone. ..beta..KA can be obtained via biological cultivation from aromatic compounds, sugars, and even waste plastics. When ..beta..KA is implemented into polymers, namely nylons, in place of adipic acid the overall material properties are increased. Specifically, the glass transition temperature (which is the measurement at which temperature a plastic softens) is increased by 69 degrees C and the water permeability is reduced by 20%. Molecular dynamic (MD) simulations reveal that the ..beta..-ketone in ..beta..KA leads to restricted movement across multiple backbone carbons when compared to adipic acid. Furthermore, MD simulations reveal that the ..beta..-ketone, when compared to no-ketone or an a-ketone, can lead to an increase in hydrogen bonding between neighboring chains. Both phenomena may help explain the molecular origin of the increases in polymer performance. Aside from performance advantages in thermomechanical performance, the production of ..beta..KA also is advantaged in its manufacturing. Although estimates for the minimum selling price (MSP) of ..beta..KA exceed adipic acid the manufacture of ..beta..KA would reduce supply chain energy and GHG emissions by >50% and >30% respectively. Further TEA analysis reveals that the MSP of ..beta..KA is primarily driven by the feedstock cost while the capital expenditures are driven largely by fermentation costs. Overall, the work presented within will be illustrative of how biomass derived molecules, namely ..beta..KA, can manifest performance advantages across multiple benchmarks and may enable a path to market for biomass-derived materials.

beta-ketoadipic acid↗
DOE OSTI2024

Hydrogen Bonding Inside Anionic Polymeric Brush Layer: Machine Learning-Driven Exploration of the Relative Roles of the Polymer Steric Effect, Charging, and Type of Screening Counterions

This paper employs a combination of all-atom molecular dynamics (MD) simulations and unsupervised machine learning (ML) for studying the water-water hydrogen bonds (HBs) inside the anionic poly-acrylic acid (PAA) brushes modeled using all-atom MD simulations. PAA brush layer with different charge fraction (f), namely f=0, f=0.25, and f=1, is considered. Water-water interactions, both inside and outside the brush layer, are represented through distinct clusters of tupules of variables representing distances associated with the interacting water molecules. While clusters representing the HBs are present for water inside and outside the brushes, several clusters representing the long-range water-water interactions are missing for the water molecules inside the highly charged (f=1) PAA brushes. More importantly, inside highly charged brushes, the edge of the clusters representing the water-water HBs is progressively shortened, as compared to that in the bulk. Both these results stem from the presence of the PAA brushes imparting the steric effect and the charge effect, or the effect associated with enhanced interactions of water molecules with PE charges and counterions, thereby disrupting the water connectivity. This water-charged-species interaction also increases the water-water HB angle, i.e., makes the water-water HBs less stable inside the highly charged PAA brush layer. The narrowing of the clusters representing the HBs and the alteration of the angle characterizing the HBs confirm that the conditions defining the water-water HBs change inside the PAA brush layer as a function of the charges on the PAA brush layer. Furthermore, we show that the use of the generic definition of HBs, as compared to using our simulation-motivated modified definition of water-water HBs, overpredict the number of water-water HBs inside the PAA brush layer. Finally, we employ this all-atom-MD-ML framework to quantify the effect of other types of screening counterions (Li + , Ca 2+ , and Y 3+ ions) in determining the water-water interactions and water-water HB properties inside the PAA brush layer. Furthermore, the findings of the present study, confirming the weakening of water-water HBs inside the PAA brush layer, points to the possibility that the water molecules will be more available for hydrating the brush layer and counterions, thereby leading to a more pronounced wetting of the PAA brush layer.

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

The importance of accounting for the Tolman correction to surface tension for nucleation and growth modeling of Fe clusters

Gibbs free energies of clusters are required for predictive modeling of cluster growth during condensation of a cooling vapor. Here, we present a straightforward method of calculating free energies of cluster formation using the data from molecular dynamics (MD) simulations. We apply this method to iron clusters having from 2 to 100 atoms. The energies obtained are verified by comparing to an MD-simulated equilibrium cluster size distribution in a sub-saturated vapor. We show that these free energies differ significantly from those obtained with a commonly used spherical cluster approximation, which relies on a surface tension coefficient of a flat surface, as it is used in the classical nucleation theory (CNT). We show that the spherical cluster approximation in CNT can be improved by using a cluster-size-dependent Tolman correction for the surface tension. The Tolman length and effective surface tension values were derived for iron clusters, and they significantly differ from the commonly used experimentally measured values. This improved approximation does not account for geometric magic number effects responsible for spikes and troughs in densities of neighbor cluster sizes. Nonetheless, it allows to more accurately model cluster formation from a cooling vapor. It better reproduces the condensation timeline, overall shape of the cluster size distribution, average cluster size, and the distribution width. In contrast, using a constant surface tension coefficient (as done in CNT) resulted in incorrect condensation dynamics and cluster size distributions. The analytical expression for cluster nucleation rate from CNT was updated to account for the size-dependence of cluster surface tension.

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

Probing High-Pressure Structural Evolution in Polyurea with In Situ Energy-Dispersive X-ray Diffraction and Molecular Dynamics Simulations

Polyurea, an elastomer with a phase-segregated microstructure, has been proven as an effective coating in defense applications. To gain a more complete understanding of the high-pressure atomic-level morphology of these phases and to validate molecular dynamics (MD) simulations, multi-angle energy-dispersive X-ray diffraction experiments were performed in situ up to pressures of ~6 GPa at room temperature. Structure factors were obtained and compared to MD simulations with an average error of less than 5% between major peak positions. The first sharp diffraction peak shifted from 4.56 Å to lower d-spacing with pressure, indicating compression between hard segments. This was further supported by the behavior of a peak at ~3.86 Å from the pair distribution function (PDF), suspected to represent π-stacking and separation between soft segments. Compression within the hard segments themselves is minimal as low-r peaks in the PDF are not greatly affected by pressure.

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

Evaluation of polyanionic cyclodextrins as high affinity binding scaffolds for fentanyl

Abstract Cyclodextrins (CDs) have been previously shown to display modest equilibrium binding affinities ( K a ~ 100–200 M -1 ) for the synthetic opioid analgesic fentanyl. In this work, we describe the synthesis of new CDs possessing extended thioalkylcarboxyl or thioalkylhydroxyl moieties and assess their binding affinity towards fentanyl hydrochloride. The optimal CD studied displays a remarkable affinity for the opioid of K a = 66,500 M −1 , the largest value reported for such an inclusion complex to date. One dimensional 1 H Nuclear Magnetic Resonance (NMR) as well as Rotational Frame Overhauser Spectroscopy (2D-ROESY) experiments supported by molecular dynamics (MD) simulations suggest an unexpected binding behavior, with fentanyl able to bind the CD interior in one of two distinct orientations. Binding energies derived from the MD simulations work correlate strongly with NMR-derived affinities highlighting its utility as a predictive tool for CD candidate optimization. The performance of these host molecules portends their utility as platforms for medical countermeasures for opioid exposure, as biosensors, and in other forensic science applications.

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

Machine Learning Reveals the Critical Interactions for SARS-CoV-2 Spike Protein Binding to ACE2

SARS-CoV and SARS-CoV-2 bind to the human ACE2 receptor in practically identical conformations, although several residues of the receptor-binding domain (RBD) differ between them. Herein, we have used molecular dynamics (MD) simulations, machine learning (ML), and free-energy perturbation (FEP) calculations to elucidate the differences in binding by the two viruses. Although only subtle differences were observed from the initial MD simulations of the two RBD–ACE2 complexes, ML identified the individual residues with the most distinctive ACE2 interactions, many of which have been highlighted in previous experimental studies. FEP calculations quantified the corresponding differences in binding free energies to ACE2, and examination of MD trajectories provided structural explanations for these differences. Lastly, the energetics of emerging SARS-CoV-2 mutations were studied, showing that the affinity of the RBD for ACE2 is increased by N501Y and E484K mutations but is slightly decreased by K417N.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Role of α → ε → α phase transformation on the spall behavior of iron at atomic scales

Shock compression of iron microstructures above a threshold stress results in a α (BCC) → ε (HCP) transformation and the propagation of the release wave brings the metal back to the α phase following the ε→α transformation. Predicting failure behavior under shock loading conditions (spallation) relies on understanding the evolution of defects in the microstructure as it undergoes the α→ε→α phase transformation. This study uses molecular dynamics (MD) simulations to investigate the role of defect evolution during α→ε→α phase transformation on the spall strength values of single-crystal (sc) Fe microstructures. The MD simulations aim to characterize the ε phase fraction formed during shock compression and the defects during shock release for variations in loading orientations and shock stresses. The simulations are carried out for loading along the [100], [110], [111], and [112] orientations and for impact velocities ranging from 600 m/s to 1 km/s. The ε phase fractions during compression and defects (dislocations, twins) characterized during spall failure show an orientation dependence that affects the spall strength values. Here, the lowest value for spall strength is observed for the <110> loading orientation that shows a high density of twinning at the spall plane, whereas the highest value is observed for the <100> orientation and is associated with a α (BCC) → γ (FCC) transformation at the spall plane. The correlations of the spall strength values with the strain rates and with the ε phase fractions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Water Structure and Dynamics Near the Surfaces of Silicalite-1

Zeolites are crystalline microporous aluminosilicates that are commonly used as industrial sorbents and membranes. Solvent structuring and dynamics near zeolite crystal surfaces are thought to influence interfacial transport and molecular recognition processes, as well as fundamental aspects of their crystallization from solution. Here, in this work, we use molecular dynamics (MD) simulations to investigate the behavior of interfacial water near the exposed (010), (100), and (101) crystal faces of silicalite-1, one of the most widely studied zeolites. The MD simulations reveal that water’s translational and orientational order is strongly influenced by the distinct corrugations, pore apertures, and functional group distributions presented on each surface. Specifically, we observe two distinct hydration layers near each surface. Water molecules in the contact layer are hydrogen bonded to two or three exposed surface silanol groups. The relative populations of water molecules that are doubly and triply hydrogen bonded to the surface strongly depend on the distribution of the exposed silanols. The interactions with exposed surface silanols also influence water’s local orientational order and distribution across each surface. We also show that the structuring of the solvent near the different faces of silicalite-1 strongly impacts the interfacial dynamics. The translational and orientational relaxation dynamics of water are slowest in the contact layers and correlate with structural ordering near each face.

42 ENGINEERING↗
DOE OSTI2026

Disentangling cation effects on ion mobility and structure in ionic liquid electrolytes

Ionic liquids (ILs) are low-temperature molten salts, where ion transport is primarily governed by ion–ion interactions. Yet, the impact of organic IL cations on critical electrolyte properties such as ion dissociation and overall transport behavior in lithium-salt-doped ILs remains poorly understood. Moreover, despite their critical role in designing IL-based electrolytes for energy storage applications, ion–ion interactions and ion-specific transport under an applied electrical potential are seldom quantified, largely due to the unique experimental and computational challenges involved. Herein, we compare transport properties obtained using 1 H, 7 Li, and 19 F pulsed-field gradient nuclear magnetic resonance (NMR) and electrophoretic NMR (eNMR) with those measured by electrochemical impedance spectroscopy. Non-equilibrium molecular dynamics (MD) simulations and eNMR confirm the presence of negatively charged [Li(TFSI) n ] (1−n) aggregates that migrate towards the positive electrode, resulting in negative lithium transference numbers. Equilibrium MD simulations reveal a vehicular Li ion transport mechanism facilitated by long-lived aggregates with Li + cations strongly bound to multiple TFSI − anions. Finally, we observe an inverse relationship between the apparent charge of the TFSI − anion in the neat IL, which is dictated by the IL cation, and Li + transport in the salt-doped systems. This highlights the opportunity to tune electrolyte performance by tailoring cation chemistry.

Li-ion batteries↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
DOE OSTI2021

Unfolding bovine α-lactalbumin with T-jump: Characterizing disordered intermediates via time-resolved x-ray solution scattering and molecular dynamics simulations

The protein folding process often proceeds through partially folded transient states. Therefore, a structural understanding of these disordered states is crucial for developing mechanistic models of the folding process. Characterization of unfolded states remains challenging due to their disordered nature, and incorporating multiple methods is necessary. Combining the time-resolved x-ray solution scattering (TRXSS) signal with molecular dynamics (MD), we are able to characterize transient partially folded states of bovine α-lactalbumin, a model system widely used for investigation of molten globule states, during its unfolding triggered by a temperature jump. We track the unfolding process between 20 μs and 70 ms and demonstrate that it passes through three distinct kinetic states. The scattering signals associated with these transient species are then analyzed with TRXSS constrained MD simulations to produce protein structures that are compatible with the input signals. Without utilizing any experimentally extracted kinetic information, the constrained MD simulation successfully drove the protein to an intermediate molten globule state; signals for two later disordered states are refined to terminal unfolded states. From our examination of the structural characteristics of these disordered states, we discuss the implications disordered states have on the folding process, especially on the folding pathway. Finally, we discuss the potential applications and limitations of this method.

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

Characterization of Heparin’s Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin’s solution conformational dynamics and structure–function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule’s underlying conformational ensemble, which then provides important information necessary for understanding structure–function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin’s conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. In this work, the anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

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

Solvent-Mediated Control of Nanocellulose Dispersion: An Integrated Computational and Experimental Investigation

Fibrillated cellulose derived from forestry feedstocks represents a renewable and high-strength materials platform for circular bioeconomies. However, its practical implementation is hindered by the irreversible aggregation of nanocellulose architectures, including cellulose nanofibers (CNFs). Solvent-based dispersion offers a simple and practical route to prevent CNF aggregation. Here, in this work, we integrate classical and enhanced sampling molecular dynamics (MD) simulations with experimental suspension rheology and atomic force microscopy (AFM) to elucidate how solvent environments tune CNF–CNF interactions and dispersion stability. CNF–CNF contact free energies computed from MD simulations reveal reduced aggregation in acetone/water, γ-valerolactone (GVL)/water, and tetrahydrofuran (THF)/water and pure acetone compared with pure water, reflecting stronger CNF-solvent relative to inter-CNF interactions. Correspondingly, CNF-solvent suspensions in these solvent systems exhibit stronger inter-fibril network structures and enhanced recovery compared to water, indicating improved CNF-solvent affinity. Liquid cell AFM imaging in acetone–water mixtures and in pure acetone further confirm the presence of well-dispersed CNFs. By combining multiscale computation with targeted experiments, this study establishes a rational framework for solvent design to achieve stable nanocellulose dispersions for high-strength biobased materials and efficient bioenergy conversion.

cellulose↗
DOE OSTI2022

Understanding interactions between biomolecules and two-dimensional nanomaterials using in silico microscopes

We report two-dimensional (2D) nanomaterials such as graphene are increasingly used in research and industry for various biomedical applications. Extensive experimental and theoretical studies have revealed that 2D nanomaterials are promising drug delivery vehicles, yet certain materials exhibit toxicity under biological conditions. So far, it is known that 2D nanomaterials possess strong adsorption propensities for biomolecules. To mitigate potential toxicity and retain favorable physical and chemical properties of 2D nanomaterials, it is necessary to explore the underlying mechanisms of interactions between biomolecules and nanomaterials for the subsequent design of biocompatible 2D nanomaterials for nanomedicine. The purpose of this review is to integrate experimental findings with theoretical observations and facilitate the study of 2D nanomaterial interaction with biomolecules at the molecular level. We discuss the current understanding and progress of 2D nanomaterial interaction with proteins, lipid membranes, and DNA based on molecular dynamics (MD) simulation. In this review, we focus on the 2D graphene nanosheet and briefly discuss other 2D nanomaterials. With the ever-growing computing power, we can image nanoscale processes using MD simulation that are otherwise not observable in experiment. We expect that molecular characterization of the complex behavior between 2D nanomaterials and biomolecules will help fulfill the goal of designing effective 2D nanomaterials as drug delivery platforms.

2D nanomaterials↗
DOE OSTI2025

Unraveling the role of temperature on the onset of ejecta formation at atomic scales

Shock wave interactions with perturbations on a free surface can lead to the inversion and growth of the perturbation and eventual ejection of a jet of material, referred to as an ejecta microjet, from shocked surfaces. Here, this study carries out large-scale molecular dynamics (MD) simulations using single-crystal Cu and Sn systems with a pre-existing groove to characterize the localized gradients in temperatures and pressures generated that render the microjet formation. MD simulations are carried out for three loading orientations (along the [001], [011], and [111] directions) and shock pressures ranging from 16 GPa to 100 GPa to understand the role of temperatures generated at the groove vertex on the formation of ejecta. The simulations suggest that the interaction of the shock wave with the groove results in a localized increase in temperatures, leading to localized softening at the groove vertex and the generation of an ejecta microjet. For Cu systems, the simulations suggest that the orientation effects on shock wave structures, velocities, and localized softening affect the ejecta formation at low pressures. Jetting is only observed when temperatures at the groove vertex are high enough to induce localized softening (close to melting temperature), and the jet velocity increases with shock pressure. In contrast, the loading orientation is rendered inconsequential for Sn systems due to the melting of the material during shock compression at the pressures chosen. The jet velocities are similar, regardless of crystal orientation in Sn systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermal boundary conductance between high thermal conductivity boron arsenide and silicon

Thermal boundary conductance (TBC) is important for heat dissipation in light-emitting diodes (LEDs). In this study, we predicted the TBC between the high thermal conductivity boron arsenide (BAs) and silicon (Si) by nonequilibrium molecular dynamics (MD) simulations. From the thermal conductivity accumulation function with respect to phonon frequency, the dominant phonon frequencies for heat conduction in BAs are extremely different from those in Si. However, our nonequilibrium MD simulations indicated that the TBC of the BAs/Si interface was still high compared to most other interfaces, even though there was a major frequency mismatch in the thermal conductivity accumulation function between BAs and Si. The primary reason for the high TBC is the overlap of phonon density of states between BAs and Si in the frequency range of 5–8 THz. The range of predicted TBC of the BAs/Si interface was between 200 and 300 MW/m2 K in the temperature range of 300–700 K, and the values of the TBC were not sensitive to the temperature. We also found that the TBCs in Si/BAs and Si/Ge interfaces were close to each other considering the simulation uncertainty. This work indicates BAs as an excellent material for heat dissipation across the interfaces.

Wei, Zhiyong (ORCID:000000018686451X)↗
DOE OSTI2025

Binding Free Energy Analysis of Colicin D, E3 and E8 to Their Respective Cognate Immunity Proteins Using Computational Simulations

Colicins are antimicrobial proteins produced by bacteria for the purpose of destroying neighboring bacteria. Colicin activity is neutralized by a specific cognate immunity protein in order to protect the host. This study investigates the structural and binding mechanisms underlying the interaction of colicin-D, -E3 and -E8 to their respective immunity proteins (ImD, Im3 and Im8) using structure prediction, molecular dynamics (MD) simulations and MM-PBSA approach of free energy calculations. High-confidence colicin-immunity (Col-Im) complex structures predicted using AlphaFold2 were subjected to MD simulations of 150 ns with GROMACS and were analyzed for the binding free energy calculation using gmx_MMPBSA. Results showed that the complex of Col_E3-Im3 exhibited the most favorable binding free energy, driven by strong van der Waals and electrostatic interactions. Col_D-ImD and Col_E8-Im8 also showed the favorable binding. Electrostatics and hydrogen bonding emerged as a key factor driving binding and stability, while polar solvation acted as a destabilizing factor across all systems. These outcomes provide an understanding of the molecular mechanisms of Col-Im systems, with potential applications for developing natural antimicrobials for food safety.

Biochemistry & Molecular Biology↗