SEARCH · Engineering Papers

Results for “interface dynamics”

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

SourceResearch connectionFrom yearThrough year

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

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Hydration structure of flat and stepped MgO surfaces

We investigate the solvation structure of flat and stepped MgO(001) in neutral liquid water using ab initio molecular dynamics based on a hybrid density functional with dispersion corrections. Our simulations show that the MgO surface is covered by a densely packed layer of mixed intact and dissociated adsorbed water molecules in a planar arrangement with strong intermolecular H-bonds. The water dissociation fractions in this layer are >20% and >30% on the flat and stepped surfaces, respectively. Slightly above the first water layer, we observe metastable OH groups perpendicular to the interface, similar to those reported in low temperature studies of water monolayers on MgO. Here, these species receive hydrogen bonds from four nearby water molecules in the first layer and have their hydrophobic H end directed toward bulk water, while their associated protons are bound to surface oxygens. The formation of these OH species is attributed to the strong basicity of the MgO surface and can be relevant for understanding various phenomena from morphology evolution and growth of (nano)crystalline MgO particles to heterogeneous catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Detailed studies of the processes in low energy H irradiation of Li and Li-compound surfaces

We have used a combination of pico-to-nano temporal/spatial scale computational physics and chemistry modeling of plasma–material interfaces in the tokamak fusion plasma edges to unravel the evolving characteristics, not readily accessible by empirical means, of lithium-, oxygen-, and hydrogen-containing materials of plasma-facing components under irradiation by hydrogen and its isotopes. In the present calculation, amorphous lithium compound surfaces containing oxygen, Li 2 O, and LiOH were irradiated by 1-100 eV particles at incident angles on the surface ranging from perpendicular to almost grazing angles. Consequential surface processes, reflection, retention, and sputtering were studied at “the same footing” and compared to earlier results from amorphous Li and LiH surfaces. In conclusion, the critical role of charging dynamics of lithium, oxygen, and hydrogen atoms in the surface chemistry during hydrogen-fuel irradiation was found to drive the kinetics and dynamics of these surfaces in unexpected ways that ultimately could have profound effects on fusion plasma confinement behavior and surface erosion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Nano-Engineered Interfaces in Dual-Layer Electrodes for Protonic Ceramic Cells with Enhanced Stability and Kinetics

Enhancing interfacial stability and charge transfer in protonic ceramic cells (PCCs) remains a critical challenge, as structural degradation and interfacial resistance often compromise durability and efficiency. Here, we report a nanoengineered dual-layer oxygen electrode architecture designed to address these limitations by introducing a fine-grained nanoparticle interfacial contact layer beneath a porous catalytic backbone. The nanoscale powders, through enhanced sintering activity, densify into a robust interfacial layer that promotes strong chemical bonding, uniform adhesion, and continuous ionic/electronic pathways with the BCZYYb electrolyte. This hierarchical architecture mitigates delamination, redistributes mechanical stress, and establishes efficient charge and mass transport channels without relying on corrosive surface treatments. Electrochemical evaluation demonstrates that the dual-layer design markedly reduces interfacial polarization resistance and accelerates electrode kinetics. Compared to the single-layer counterpart, the architecture achieves a peel strength of 44.53 N/cm 2 , a 40% improvement in peak power density (0.96 W cm –2 at 600 °C), and a 130% enhancement in electrolysis current density (4.78 A cm –2 at 1.57 V). Faradaic efficiency remains as high as 88% under high steam concentrations, underscoring minimal charge loss during practical operation. Notably, the electrode retains stability across 450–600 °C and under transient voltage cycling, with impedance spectra confirming suppressed interfacial resistance growth over prolonged use. These results highlight nanoscale interface engineering as a powerful route to enhance both mechanical robustness and electrochemical kinetics in PCCs. The demonstrated scalability and durability of this architecture provide a versatile platform for advancing solid-state electrochemical systems, including reversible fuel cells and high-efficiency hydrogen production technologies.

Faradaic efficiency↗

DOE OSTI2024

SPARC-X: Quantum simulations at extreme scale - reactive dynamics from first principles

We have developed the massively parallel electronic structure code SPARC-X: a computational framework for performing Kohn-Sham Density Functional Theory (DFT) calculations that can scale linearly with the number of atoms in the system, while being able to leverage petascale and emerging exascale parallel computers to study chemical phenomena at unprecedented length and time scales. SPARC-X exploits a recent breakthrough in electronic structure methodologies: systematically improvable, strictly local, orthonormal, discontinuous real-space bases that efficiently and systematically capture the local chemistry of the system. With further adaptation using new machine-learning techniques and the use of the massively parallel Spectral Quadrature (SQ) electronic structure method, the algorithmic complexity and prefactor associated with DFT calculations involving semilocal as well as hybrid functionals are dramatically reduced. Using petascale computational resources, SPARC-X enables quantum mechanical simulations at length and time scales previously accessible only by empirical approaches, e.g., 1,000,000 atoms for a few picoseconds using semilocal functionals or 1,000 atoms for a few picoseconds using hybrid functionals. Using exascale resources, the sizes and times targeted are two orders of magnitude larger. Such a capability has applications in a wide variety of chemical sciences, including reactive interfaces where large length- and/or long time-scales are needed and traditional force fields fail. This is particularly important in dynamic catalysis, where bond breaking and formation must be understood in detail. We developed, tested, and employed the SPARC-X framework to understand the photocatalytic properties of TiO 2 nanoparticles, revealing finite size effects that cannot be captured with standard model systems or functionals. This integrated development and application strategy ensures that SPARC-X remains a robust, efficient, and scalable software package for quantum simulations on current petascale and emerging exascale computing resources.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Methanol at Water–TiO 2 Interfaces: Free Energies of Water and Methanol Dissociation

Methanol adsorption on TiO 2 surfaces has long been studied due to its role in enhancing photocatalytic hydrogen evolution, yet how it modulates surface chemistry under aqueous conditions remains little understood. Using molecular dynamics with an ab initio-based deep neural network potential, we find that methanol adsorption induces markedly different effects on the aqueous surfaces of anatase and rutile, the two common phases of TiO 2 . In anatase, methanol adsorption significantly enhances water dissociation, which is otherwise rare at the neat water interface. This enhancement arises from an alternative dissociation pathway mediated by surface-bound methoxyl groups. In contrast, methanol adsorption tends to suppress water dissociation on rutile, replacing it with thermodynamically favored methanol dissociation. Overall, methanol adsorption in an aqueous environment alters not only the availability of key reactive intermediates involved in hydrogen evolution but also the hydrogen source, which turns out to be primarily methanol on rutile, whereas both water and methanol are consumed on anatase. These results provide mechanistic insights into the coupled roles of organic adsorbates and water at photocatalytic interfaces, with implications on how methanol enhances the activity of H 2 evolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Dynamic and weak electric double layers in ultrathin nanopores

The unique properties of aqueous electrolytes in ultrathin nanopores have drawn a great deal of attention in a variety of applications, such as power generation, water desalination, and disease diagnosis. Inside the nanopore, at the interface, properties of ions differ from those predicted by the classical ionic layering models (e.g., Gouy–Chapman electric double layer) when the thickness of the nanopore approaches the size of a single atom (e.g., nanopores in a single-layer graphene membrane). Here, using extensive molecular dynamics simulations, the structure and dynamics of aqueous ions inside nanopores are studied for different thicknesses, diameters, and surface charge densities of carbon-based nanopores [ultrathin graphene and finite-thickness carbon nanotubes (CNTs)]. The ion concentration and diffusion coefficient in ultrathin nanopores show no indication of the formation of a Stern layer (an immobile counter-ionic layer) as the counter-ions and nanopore atoms are weakly correlated in time compared to the strong correlation observed in thick nanopores. The weak correlation observed in ultrathin nanopores is indicative of a weak adsorption of counter-ions onto the surface compared to that of thick pores. The vanishing counter-ion adsorption (ion–wall correlation) in ultrathin nanopores leads to several orders of magnitude shorter ionic residence times (picoseconds) compared to the residence times in thick CNTs (seconds). The results of this study will help better understand the structure and dynamics of aqueous ions in ultrathin nanopores.

2D materials↗

DOE OSTI2021

Prediction of Bi 2 Te 3 -Sb 2 Te 3 Interfacial Conductance and Superlattice Thermal Conductivity Using Molecular Dynamics Simulations

Bismuth telluride (Bi 2 Te 3 ) and its alloys with antimony telluride (Sb 2 Te 3 ) have long been considered to be the best room-temperature bulk thermoelectric (TE) materials. In recent decades, proof-of-concept demonstrations on Bi 2 Te 3 -Sb 2 Te 3 nanostructures have shown high TE performance due to reduction in lattice thermal conductivities. Particularly, ultra-low thermal conductivities have been observed in Bi 2 Te 3 -Sb 2 Te 3 1D superlattices, leading to thermoelectric figures of merit (ZT) as high as 2.4. In contrast, very few computational studies have been performed to provide insight into the phonon transport across these nanostructures. In this work, we use non-equilibrium molecular dynamics simulations with previously developed force fields to simulate thermal transport across Bi 2 Te 3 -Sb 2 Te 3 interfaces and superlattices. We first calculate the thermal conductance associated with a Bi 2 Te 3 -Sb 2 Te 3 interface across a temperature range of 200–400 K. Furthermore, the values are also compared with thermal conductances calculated by a modified Landauer transport formalism using phonon transmission coefficients obtained from the diffuse mismatch model. Our results show that inelastic scattering processes contribute to an increase in interfacial thermal conductance at higher temperatures. Finally, we calculate the thermal conductivities of Bi 2 Te 3 -Sb 2 Te 3 superlattices with varying period lengths from 2 to 18 nm. A minimum thermal conductivity of 0.27 W/mK is observed at a period length of 4 nm, which is attributed to the competition between incoherent and coherent phonon transport regimes. In comparison with previous experimental measurements in the literature, our results show good agreement with respect to the range of thermal conductivity values and the period length corresponding to the minimum superlattice thermal conductivity.

42 ENGINEERING↗

DOE OSTI2022

Foundational Science to Accelerate Nuclear Energy Innovation [Brochure]

The foundational science gaps inhibiting the advancement of nuclear energy technologies are identified and tackled in five priority research opportunities. These opportunities pave the way to accelerate the development and ultimately the adoption of new nuclear energy systems. They include the fundamental aspects of ion-electron interactions; novel properties of next-generation coolants and solvents; interfacial dynamics, not only in solids, but in other aspects of nuclear reactors; novel operando and in-situ monitoring and sensing; and artificial intelligence to accelerate condensed phases discovery. Building on the foundation established by previous BES workshops, these opportunities encompass recent advances in fundamental knowledge and focus on the experimental and computational methods needed to resolve major technical challenges for nuclear energy technologies. Through developing fundamental scientific insight as well as pushing the frontiers of modeling complex systems and probing the operation of materials and chemical systems in extreme environments, research motivated by the priorities identified here will further develop the promise, potential, and utilization of nuclear energy for a clean energy future. The PROs are as follows: (1) Master complex electronic structures to tailor thermochemical reactivity, transport, and microstructural evolution; (2) Interrogate and direct the physics and chemistry underpinning next-generation coolants and solvents; (3) Elucidate and control the underlying physics and chemistry of interfaces in complex nuclear environments; (4) Bridge multi-fidelity multi-resolution experiments, computational modeling, and data science to control dynamic behavior; and (5) Harness artificial intelligence to design inherently resilient condensed phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

A Fortran–Python interface for integrating machine learning parameterization into earth system models

Abstract. Parameterizations in earth system models (ESMs) are subject to biases and uncertainties arising from subjective empirical assumptions and incomplete understanding of the underlying physical processes. Recently, the growing representational capability of machine learning (ML) in solving complex problems has spawned immense interests in climate science applications. Specifically, ML-based parameterizations have been developed to represent convection, radiation, and microphysics processes in ESMs by learning from observations or high-resolution simulations, which have the potential to improve the accuracies and alleviate the uncertainties. Previous works have developed some surrogate models for these processes using ML. These surrogate models need to be coupled with the dynamical core of ESMs to investigate the effectiveness and their performance in a coupled system. In this study, we present a novel Fortran–Python interface designed to seamlessly integrate ML parameterizations into ESMs. This interface showcases high versatility by supporting popular ML frameworks like PyTorch, TensorFlow, and scikit-learn. We demonstrate the interface's modularity and reusability through two cases: an ML trigger function for convection parameterization and an ML wildfire model. We conduct a comprehensive evaluation of memory usage and computational overhead resulting from the integration of Python codes into the Fortran ESMs. By leveraging this flexible interface, ML parameterizations can be effectively developed, tested, and integrated into ESMs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Stability and Fidelity of Power-Hardware-in-the-Loop Using NREL's CGI

The presentation focuses on dynamic analysis of power hardware-in-the-loop (PHIL) experiments using National Renewable Energy Laboratory Flatirons Campus' grid simulator - 7MVA Controllable Grid Interface (CGI). A control diagram of typical PHIL experiment is described, and transfer functions of critical components are explained. An impedance-based definition of PHIL experiment bandwidth is defined using open loop transfer functions of the system. Further an example resistive divider circuit is introduced that is used to explain implication of PHIL on accuracy of simulation by comparing measured results against ideal expected signals. At the end the real-life PHIL simulation using Maui power system model is described showing various oscillation modes observed and discussing the meaning of these in context of dynamic bandwidth of simulation.

grid simulator↗

DOE OSTI2023

Hemp fiber reinforced dual dynamic network vitrimer biocomposites with direct incorporation of amino silane

Natural fiber composites are inexpensive and renewable alternatives to traditional fiber reinforced polymers (FRPs). However, existing natural fiber composites primarily rely on petrochemical based thermosetting matrices, which are difficult to recycle due to their stable crosslinked network structures. To address this challenge, it is desirable to develop natural fiber composites with inherently recyclable biobased matrices. In this study, we developed a dual dynamic network vitrimer matrix from hempseed oil and limonene derivatives and demonstrated its application in hemp fiber reinforced composites. To improve the weak interface between hemp fiber and the polymer matrix, a common challenge in the realm of natural fiber composites, we directly incorporated amino silane into the vitrimer matrix. The amino silane participated in the polymer network, increasing the crosslink density and toughening the matrix, as evidenced by a significantly improved impact strength from 3.5 kJ/m 2 to 10.3 kJ/m 2 . Moreover, the incorporation of amino silane resulted in a lower water absorption by 11% in 7-day soaking for the composites, demonstrating improved fiber/matrix interfacial interaction. Furthermore, our biobased vitrimer matrix exhibits both imine and hydroxy-ester dynamic bonds, which enable the recycling of the biocomposite through a mild and cost-effective aminolysis process (100°C, 3h, ambient pressure). The decomposed polymer matrix was successfully reused as a polyol for polyurethane adhesives, and the surface morphology of recovered hemp fibers was analyzed and compared with that of the original fibers. Furthermore, these findings will help broaden the use of vitrimers for practical natural fiber composite applications, raise awareness of the problem of recycling natural fiber composite waste, and shed light on the interfacial challenge of natural fiber reinforced vitrimer composites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Self-Assembled Oligomers Facilitate Amino Acid-Driven CO 2 Capture at the Air–Aqueous Interface

Direct air capture of CO 2 using amino acid absorbents, such as glycine or sarcosine, is constrained by the relatively slow mass transfer of CO 2 through the air–aqueous interface. Our recent study showed a marked improvement in CO 2 capture by introducing CO 2 -permeable oligo-dimethylsiloxane (ODMS-MIM + ) oligomers with cationic (imidazolium, MIM + ) headgroups. Here, in this work, we have employed all-atom molecular dynamics simulations in combination with subensemble analysis using network theory to provide a detailed molecular picture of the behavior of CO 2 and the glycinate anions (Gly – ) at the ODMS-MIM + decorated air–aqueous interfaces. We show that the cationic head groups of the surfactants enhance the concentration and lifetime of Gly – in the interfacial region, while ODMS tails promote the physisorption of CO 2 in the interfacial region. Together, these two factors increase the effective region of contact and the probability of interactions between CO 2 and Gly – compared to that of the pure air–aqueous interface. The fundamental insights gained in this work establish essential foundations for developing hybrid systems with oligomer-decorated interfaces to maximize the overall CO 2 capture rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

PigmentHunter: A point-and-click application for automated chlorophyll-protein simulations

Chlorophyll proteins (CPs) are the workhorses of biological photosynthesis, working together to absorb solar energy, transfer it to chemically active reaction centers, and control the charge-separation process that drives its storage as chemical energy. Yet predicting CP optical and electronic properties remains a serious challenge, driven by the computational difficulty of treating large, electronically coupled molecular pigments embedded in a dynamically structured protein environment. To address this challenge, we introduce here an analysis tool called PigmentHunter, which automates the process of preparing CP structures for molecular dynamics (MD), running short MD simulations on the nanoHUB.org science gateway, and then using electrostatic and steric analysis routines to predict optical absorption, fluorescence, and circular dichroism spectra within a Frenkel exciton model. Inter-pigment couplings are evaluated using point-dipole or transition-charge coupling models, while site energies can be estimated using both electrostatic and ring-deformation approaches. The package is built in a Jupyter Notebook environment, with a point-and-click interface that can be used either to manually prepare individual structures or to batch-process many structures at once. Here, we illustrate PigmentHunter’s capabilities with example simulations on spectral line shapes in the light harvesting 2 complex, site energies in the Fenna–Matthews–Olson protein, and ring deformation in photosystems I and II.

14 SOLAR ENERGY↗

DOE OSTI2026

Theoretical Insights into Reaction-Induced Transformation and Tuning of Catalytic Behavior in Heterogenous Catalysis

Reaction-induced transformations in heterogenous catalysis represent diverse phenomena that challenge traditional views of static catalyst surfaces. From surface adsorbate dynamics, atomic rearrangements, to composition and phase transitions, these processes reveal the profound differences between idealized model systems under ultrahigh vacuum and the complex, evolving interfaces that govern real catalytic behaviors under reaction conditions. Here, this perspective reviews recent theoretical efforts to provide atomic-level mechanistic insights into significant reaction-induced transformations and their impact on catalytic activity and selectivity. It underscores the need for an integrated framework that combines predictive simulations with operando characterization to uncover active sites and mechanisms under realistic operating conditions. Achieving this requires accelerating existing simulations to fully capture diverse reaction-induced surface dynamics, enabling scalable and accurate modeling of catalysts as condition-dependent, dynamically evolving systems. Such approaches are critical to bridge the gap between theory and practice, offering a pathway to more impactful and predictive catalyst design.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

Local conformations and heterogeneities in structures and dynamics of isotactic polypropylene adsorbed onto carbon fiber

Carbon fiber (CF) reinforced polymers (CFRPs) have experienced widespread use in various industries. One of the most important parameters that controls the macroscopic property of CFRPs is the interface between a polymer matrix and CF. There is growing evidence to suggest the formation of a bound polymer layer (BPL), i.e., polymer chains that physically adsorb on a filler surface. However, this interface is always in contact with the thicker part of a polymer matrix, rendering its understanding a difficult task. Therein, we use CF-reinforced isotactic polypropylene (iPP) as a rational CFRP. To characterize the BPL on the CF surface, we extracted it from the CFRP using solvent-rinsing with p-xylene. The physical and thermal properties of the BPL were characterized by differential scanning calorimetry and thermogravimetric analysis, while its microscopic structures and dynamics were probed by small-angle neutron scattering and quasi-elastic neutron scattering (QENS) techniques. Subsequently, we employed atomistic molecular dynamics (MD) simulations to complement the QENS results above the bulk melting temperature and reveal details that were experimentally inaccessible. In this work, we observed that the degree of crystallinity of the BPL was quite lower than the bulk, while the melting temperature of the BPL remained the same as the bulks. Within the given length and time scales probed by QENS and MD, we also observed that most of the bound chains were mobile, with the formation of a high-density region (less than 1 nm in thickness) near the CF surface. The segmental dynamics of the bound chains probed by both QENS and MD were also much faster than those of the free chains, possibly due to the presence of a free surface region at the topmost surface of the BPL. Furthermore, the MD results demonstrated that the backbone chains and side groups lie nearly flat on the CF surface, which is the driving force for the flattening process of the iPP BPL to overcome the conformational entropy loss in the total free energy.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Mechanism of ion pairing–induced co-ion penetration into the stern layer

In classical electric double layer (EDL) theory, co-ions—ions carrying the same charge as the electrode—are assumed to be excluded from the Stern layer because of electrostatic repulsion. However, ion pairing with adsorbed counterions may enable co-ions to access the EDL, a phenomenon that remains underexplored. Key questions therefore remain unresolved, including how deeply co-ions penetrate, how ion-specific properties govern them, and how ion pairing within the EDL differs from bulk solution. Here, in this study, we employ molecular dynamics simulations combined with well-tempered metadynamics to quantify co-ion penetration and ion-pairing thermodynamics at positively charged Mg-Al layered double hydroxide interfaces (σ = +20.6 μC/cm 2 ) in aqueous chloride electrolytes across 0.005–3 M range. We reveal a concentration-driven crossover: at dilute conditions, counterion–co-ion contact ion pairing (CIP) is thermodynamically unfavorable or only weakly favorable (ΔG CIP = +0.12 kcal/mol for Li + to −0.42 kcal/mol for K + ), leading to co-ion exclusion. As concentration increases (≥0.5 M), enhanced screening reduces surface repulsion and stabilizes CIP thermodynamics (ΔG CIP = −0.89 kcal/mol for Li + to −1.05 kcal/mol for K + at 3 M), enabling monovalent co-ions to penetrate the Stern layer following the hydration hierarchy K + > Na + > Li + , while divalent co-ions (Mg 2+ , Ca 2+ ) remain in the diffuse layer, forming only solvent-separated pairs. Comparison with bulk solution shows interfacial ion pairing is suppressed at low concentrations but converges to bulk-like thermodynamics at high concentrations, indicating surface effects vanish under strong screening. These results establish a quantitative framework linking ion-pairing thermodynamics to co-ion penetration and providing strategies to design electrochemical interfaces.

Co-ion penetration↗

DOE OSTI2026

Ionic Interdiffusion at Cathode|Solid-Electrolyte Interface: A Machine Learning–Assisted Multiscale Investigation and Mitigation Strategies

Future lithium batteries are expected to use solid electrolytes to achieve higher energy density and fast charge capabilities. However, most solid electrolytes are thermodynamically unstable against layered oxide cathodes. In this study, the stability of LiCoO2 (LCO) cathode with Li10GeP2S12 (LGPS) solid electrolyte is investigated using ab initio molecular dynamics (AIMD) and machine learning molecular dynamics (MLMD). The propensity of ionic interdiffusion, formation of a passivating interphase layer, and corresponding decay in cell performance is addressed using a continuum model. Large-scale MLMD simulations confirm that the LCO|LGPS interface permits interdiffusion of cobalt (Co) and other ionic species, leading to the formation and growth of a resistive interphase and to dramatic capacity fade even in the first cycle. We examine the literature evidence that incorporating a thin layer of LiNb0.5Ta0.5O3 (LNTO) between LCO and LGPS prevents the interdiffusion of ions. Atomistic simulations suggest that substituting lithium (Li) in LNTO with Co is thermodynamically unfavorable, thereby inhibiting ionic interdiffusion. The stable Nb5+/Ta5+ states form a rigid metal-oxide framework, which consequently also prevents the substitution of niobium (Nb) or tantalum (Ta). However, continuum-level analysis suggests that the higher mechanical stiffness of LNTO can lead to interfacial delamination between the LCO and LNTO. This phenomenon reduces the effectiveness of the protective layer. This paper, therefore, highlights the need to develop novel interlayers that balance low ionic interdiffusion with low mechanical stiffness.

Ncube, Musawenkosi K.↗

DOE OSTI2023

Integrated data-driven and experimental approaches to accelerate lead optimization targeting SARS-CoV- 2 main protease

Identification of potential therapeutic candidates can be expedited by integrating computational modeling with domain aware machine learning (ML) approaches followed by experimental validation. Generative deep learning models have been recently developed that can generate thousands of new candidates, but their physiochemical properties are typically not optimized. Using our deep learning models and a scaffold as a starting point, we generated tens of thousands of compounds for SARS-CoV-2 M pro that preserve the core scaffold. Here we utilized and implemented several computational tools such as structural alert and toxicity analysis, high throughput virtual screening, ML-based 3D quantitative structure–activity relationships, multi-parameter optimization, and graph neural networks on libraries of generated candidates to predict biological activity and binding affinity a priori. From these collective computational results, eight promising candidates were identified and tested experimentally using Native Mass Spectrometry (MS) and FRET-based functional assays. Two compounds, with quinazoline-2-thiol and acetylpiperidine core moiety showed IC 50 values in the low micromolar range: 2.95±0.0017 µM and 3.41±0.0015 µM, respectively. The molecular dynamics simulations further highlight that binding of these compounds results in allosteric modulations in the chain B and the interface domains of the M pro . The key fragments from these top hits can be used as input for closed loop lead optimization in the integrated pipeline.

60 APPLIED LIFE SCIENCES↗