Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protein dynamics”

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

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

At least 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Permeability-Engineered Compartmentalization Enables In Vitro Reconstitution of Sustained Synthetic Biology Systems

In nature, biological compartments such as cells rely on dynamically controlled permeability for matter exchange and complex cellular activities. Likewise, the ability to engineer compartment permeability is crucial for in vitro systems to gain sustainability, robustness, and complexity. However, rendering in vitro compartments such a capability is challenging. Here, a facile strategy is presented to build permeability-configurable compartments, and marked advantages of such compartmentalization are shown in reconstituting sustained synthetic biology systems in vitro. Through microfluidics, the strategy produces micrometer-sized layered microgels whose shell layer serves as a sieving structure for biomolecules and particles. In this configuration, the transport of DNAs, proteins, and bacteriophages across the compartments can be controlled an guided by a physical model. Through permeability engineering, a compartmentalized cell-free protein synthesis system sustains multicycle protein production; ≈100 000 compartments are repeatedly used in a five-cycle synthesis, featuring a yield of 2.2 mg mL -1 . Further, the engineered bacteria-enclosing compartments possess near-perfect phage resistance and enhanced environmental fitness. In a complex river silt environment, compartmentalized whole-cell biosensors show maintained activity throughout the 32 h pollutant monitoring. It is anticipated that permeability-engineered compartmentalization should pave the way for practical synthetic biology applications such as green bioproduction, environmental sensing, and bacteria-based therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Liquid Structure and Transport Properties of the Deep Eutectic Solvent Ethaline

A range of techniques including physical property measurements, neutron scattering experiments, ab initio molecular dynamics, and classical molecular dynamics simulations are used to probe the structural, thermodynamic, and transport properties of a deep eutectic solvent comprised of a 1:2 molar ratio of choline chloride and ethylene glycol. This mixture, known as Ethaline, has many desirable properties for use in a range of applications, and therefore, understanding its liquid structure and transport properties is of interest. Simulation results are able to capture experimental densities, diffusivities, viscosities, and structure factors extremely well. Here, the solvation environment is dynamic and dominated by different hydrogen bonding interactions. Dynamic heterogeneities resulting from hydrogen bonding interactions are quantified. Rotational dynamics of molecular dipole moments of choline and ethylene glycol are computed and found to exhibit a fast and slow mode.

25 ENERGY STORAGE↗
DOE OSTI2025

Switchable client specificity in a dual functional chaperone coordinates light-harvesting complex biogenesis

The proper assembly of light-harvesting complexes (LHCs) is critical for photosynthesis and requires the biogenesis of light-harvesting chlorophylla,b-binding proteins (LHCPs) to be coordinated with chlorophyll (Chl) biosynthesis. The mechanism underlying this coordination is not well understood. Here, we show that a conserved molecular chaperone, chloroplast signal recognition particle 43-kDa protein (cpSRP43), provides a molecular thermostat that helps maintain this coordination. cpSRP43 undergoes a conformational rearrangement between a well-folded closed state and a partially disordered open state. Closed cpSRP43 is dedicated to the biogenesis of LHCPs, whereas open cpSRP43 protects multiple Chl biosynthesis enzymes from heat-induced destabilization. Rising temperature shifts cpSRP43 to the open state, enabling it to protect heat-destabilized Chl biosynthesis enzymes. Our results reveal the molecular basis of a posttranslational mechanism for the thermoadaptation of LHC biogenesis. They also demonstrate how an adenosine triphosphate–independent chaperone uses conformational dynamics to switch its activity and client selectivity, thereby adapting to different proteostatic demands under shifting environmental conditions.

Science & Technology - Other Topics↗
DOE OSTI2020

FLIMJ: An open-source ImageJ toolkit for fluorescence lifetime image data analysis

In the field of fluorescence microscopy, there is continued demand for dynamic technologies that can exploit the complete information from every pixel of an image. One imaging technique with proven ability for yielding additional information from fluorescence imaging is Fluorescence Lifetime Imaging Microscopy (FLIM). FLIM allows for the measurement of how long a fluorophore stays in an excited energy state, and this measurement is affected by changes in its chemical microenvironment, such as proximity to other fluorophores, pH, and hydrophobic regions. This ability to provide information about the microenvironment has made FLIM a powerful tool for cellular imaging studies ranging from metabolic measurement to measuring distances between proteins. The increased use of FLIM has necessitated the development of computational tools for integrating FLIM analysis with image and data processing. To address this need, we have created FLIMJ, an ImageJ plugin and toolkit that allows for easy use and development of extensible image analysis workflows with FLIM data. Built on the FLIMLib decay curve fitting library and the ImageJ Ops framework, FLIMJ offers FLIM fitting routines with seamless integration with many other ImageJ components, and the ability to be extended to create complex FLIM analysis workflows. Building on ImageJ Ops also enables FLIMJ’s routines to be used with Jupyter notebooks and integrate naturally with science-friendly programming in, e.g., Python and Groovy. We show the extensibility of FLIMJ in two analysis scenarios: lifetime-based image segmentation and image colocalization. We also validate the fitting routines by comparing them against industry FLIM analysis standards.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Microbial biofilms as living photoconductors due to ultrafast electron transfer in cytochrome OmcS nanowires

Light-induced microbial electron transfer has potential for efficient production of value-added chemicals, biofuels and biodegradable materials owing to diversified metabolic pathways. However, most microbes lack photoactive proteins and require synthetic photosensitizers that suffer from photocorrosion, photodegradation, cytotoxicity, and generation of photoexcited radicals that are harmful to cells, thus severely limiting the catalytic performance. Therefore, there is a pressing need for biocompatible photoconductive materials for efficient electronic interface between microbes and electrodes. Here we show that living biofilms of Geobacter sulfurreducens use nanowires of cytochrome OmcS as intrinsic photoconductors. Photoconductive atomic force microscopy shows up to 100-fold increase in photocurrent in purified individual nanowires. Photocurrents respond rapidly (<100 ms) to the excitation and persist reversibly for hours. Femtosecond transient absorption spectroscopy and quantum dynamics simulations reveal ultrafast (~200 fs) electron transfer between nanowire hemes upon photoexcitation, enhancing carrier density and mobility. Our work reveals a new class of natural photoconductors for whole-cell catalysis.

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

Structure and dynamics of hydrodynamically interacting finite-size Brownian particles in a spherical cavity: Spheres and cylinders

The structure and dynamics of confined suspensions of particles of arbitrary shape are of interest in multiple disciplines from biology to engineering. Theoretical studies are often limited by the complexity of long-range particle-particle and particle-wall forces, including many-body fluctuating hydrodynamic interactions. Here, we report a computational study on the diffusion of spherical and cylindrical particles confined in a spherical cavity. We rely on an immersed-boundary general geometry Ewald-like method to capture lubrication and long-range hydrodynamics and include appropriate non-slip conditions at the confining walls. A Chebyshev polynomial approximation is used to satisfy the fluctuation-dissipation theorem for the Brownian suspension. We explore how lubrication, long-range hydrodynamics, particle volume fraction, and shape affect the equilibrium structure and the diffusion of the particles. It is found that once the particle volume fraction is greater than 10%, the particles start to form layered aggregates that greatly influence particle dynamics. Hydrodynamic interactions strongly influence the particle diffusion by inducing spatially dependent short-time diffusion coefficients, stronger wall effects on the particle diffusion toward the walls, and a sub-diffusive regime-caused by crowding-in the long-time particle mobility. The level of asymmetry of the cylindrical particles considered here is enough to induce an orientational order in the layered structure, decreasing the diffusion rate and facilitating a transition to the crowded mobility regime at low particle concentrations. Our results offer fundamental insights into the diffusion and distribution of globular and fibrillar proteins inside cells.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Consequences of Altered Root Nutrient Uptake for Soil Carbon Stabilization (Final Report)

Objectives: The primary objectives of this project are to improve our understanding of tropical forest belowground processes, and to increase representation of this biome in ecosystem-scale and global C cycle models. Belowground dynamics present a major source of uncertainty inhibiting our ability to predict C cycle responses to climate change. Representation of tropical forests in these models is particularly lacking, even though these ecosystems hold >25% of terrestrial C stocks. While the majority of humid tropical forests exist on relatively infertile soils (i.e. poor in rock-derived nutrients), enormous gradients in soil nutrient availability exist at landscape-scales due to shifts in geology and soil order. Soil fertility is likely to greatly influence how tropical soil C storage will respond to the declines in precipitation predicted for the tropics. In particular, root dynamics vary across soil fertility gradients, with roots representing the major input of C to subsurface soils. Root characteristics related to soil fertility include biomass, turnover, C exudates, tissue chemistry, and nutrient uptake rates, with each of these likely sensitive to changes in moisture. The proposed project will fully integrate field research with an existing ecosystem model to assess potential effects of drying on belowground tropical C dynamics across a range of soil fertilities. Project Description:. This project will identify key linkages among soil fertility, root nutrient uptake, root dynamics, and soil C storage. The hypothesized framework for understanding these linkages is: Tropical rainforest roots in relatively fertile soils have lesser root biomass, faster turnover, fewer exudates, and improved root tissue quality relative to infertile sites, all resulting from increased root nutrient uptake. Ultimately, these root characteristics in fertile soils promote proportionally greater long-term soil C storage in organo-mineral associations. Specifically, lower root exudates and improved root tissue quality in fertile soils reduce microbial respiration of extant soil C (i.e. “priming) and increase microbial C use efficiency, leading to sorption of protein-rich microbial and plant organic matter on soil mineral surfaces. However, the larger pools of mineral-associated C in fertile tropical soils are more vulnerable to decreased rainfall than smaller C pools in infertile soils, in part because of the lesser root biomass, faster turnover, and softer, more labile root tissues. Also, the shallower lateral rooting structures typical in fertile soils experience greater desiccation and death than deeper roots in infertile sites, with little short-term adaptation to drying.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Comparing Zwitterionic and PEG Exteriors of Polyelectrolyte Complex Micelles

A series of model polyelectrolyte complex micelles (PCMs) was prepared to investigate the consequences of neutral and zwitterionic chemistries and distinct charged cores on the size and stability of nanocarriers. Using aqueous reversible addition-fragmentation chain transfer (RAFT) polymerization, we synthesized a well-defined diblock polyelectrolyte system, poly(2-methacryloyloxyethyl phosphorylcholine methacrylate)-block-poly((vinylbenzyl) trimethylammonium) (PMPC-PVBTMA), at various neutral and charged block lengths to compare directly against PCM structure–property relationships centered on poly(ethylene glycol)-block-poly((vinylbenzyl) trimethylammonium) (PEG-PVBTMA) and poly(ethylene glycol)-block-poly(l-lysine) (PEG-PLK). After complexation with a common polyanion, poly(sodium acrylate), the resulting PCMs were characterized by dynamic light scattering (DLS) and small angle X-ray scattering (SAXS). We observed uniform assemblies of spherical micelles with a diameter ~1.5–2× larger when PMPC-PVBTMA was used compared to PEG-PLK and PEG-PVBTMA via SAXS and DLS. In addition, PEG-PLK PCMs proved most resistant to dissolution by both monovalent and divalent salt, followed by PEG-PVBTMA then PMPC-PVBTMA. All micelle systems were serum stable in 100% fetal bovine serum over the course of 8 h by time-resolved DLS, demonstrating minimal interactions with serum proteins and potential as in vivo drug delivery vehicles. This thorough study of the synthesis, assembly, and characterization of zwitterionic polymers in PCMs advances the design space for charge-driven micelle assemblies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Molecular basis for lipid recognition by the prostaglandin D 2 receptor CRTH2

Significance Many amphipathic lipids, including cannabinoids, lysophospholipids, prostanoids, and leukotrienes, bind to G protein–coupled receptors (GPCRs) through the lipid bilayer. How these lipid ligands get anchored to their receptors and enter ligand-binding pockets is largely elusive. Our results revealed a distinctive lipid recognition mechanism for the prostaglandin D 2 (PGD 2 ) receptor CRTH2 that is not shared by other prostanoid receptors. Based on the structural analysis, we propose two representative binding modes of lipids in GPCRs. We also identified critical structural motifs and conformational changes of CRTH2 for capturing lipids from the lipidic environment and lipid transition to the binding pocket, thereby offering unprecedented molecular insight into the dynamic lipid-binding processes for GPCRs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Data-driven prediction of α IIb β 3 integrin activation paths using manifold learning and deep generative modeling

The integrin heterodimer is a transmembrane protein critical for driving cellular process and is a therapeutic target in the treatment of multiple diseases linked to its malfunction. Activation of integrin involves conformational transitions between bent and extended states. Some of the conformations that are intermediate between bent and extended states of the heterodimer have been experimentally characterized, but the full activation pathways remain unresolved both experimentally due to their transient nature and computationally due to the challenges in simulating rare barrier crossing events in these large molecular systems. An understanding of the activation pathways can provide new fundamental understanding of the biophysical processes associated with the dynamic interconversions between bent and extended states and unveil new putative therapeutic targets. In this work, we apply nonlinear manifold learning to coarse-grained molecular dynamics simulations of bent, extended, and two intermediate states of αI I b β3 integrin to learn a low-dimensional embedding of the configurational phase space. We then train deep generative models to learn an inverse mapping between the low-dimensional embedding and high-dimensional molecular space and use these models to interpolate the molecular configurations constituting the activation pathways between the experimentally characterized states. Furthermore, this work furnishes plausible predictions of integrin activation pathways and reports a generic and transferable multi-scale technique to predict transition pathways for biomolecular systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Molecular dynamics analysis of the binding of human interleukin-6 with interleukin-6 α-receptor

Human interleukin-6 (hIL-6) is a multifunctional cytokine that regulates immune and inflammatory responses in addition to metabolic and regenerative processes and cancer. hIL-6 binding to the IL-6 receptor (IL-6Rα) induces homodimerization and recruitment of the glycoprotein (gp130) to form a hexameric signaling complex. Anti-IL-6 and IL-6R antibodies are clinically approved inhibitors of IL-6 signaling pathway for treating rheumatoid arthritis and Castleman's disease, respectively. There is a potential to develop novel small molecule IL-6 antagonists derived from understanding the structural basis for IL-6/IL-6Rα interactions. In this work, we combine homology modeling with extensive molecular dynamics (MD) simulations to examine the association of hIL-6 with IL-6Rα. A comparison with MD of apo hIL-6 reveals that the binding of hIL-6 to IL-6Rα induces structural and dynamic rearrangements in the AB loop region of hIL-6, disrupting intraprotein contacts and increasing the flexibility of residues 48 to 58 of the AB loop. In contrast, due to the involvement of residues 59 to 78 in forming contacts with the receptor, these residues of the AB loop are observed to rigidify in the presence of the receptor. The binary complex is primarily stabilized by two pairs of salt bridges, Arg181 (hIL-6)- Glu182 (IL-6Rα) and Arg184 (hIL-6)- Glu183 (IL-6Rα) as well as hydrophobic and aromatic stacking interactions mediated essentially by Phe residues in both proteins. An interplay of electrostatic, hydrophobic, hydrogen bonding, and aromatic stacking interactions facilitates the formation of the hIL-6/IL-6Rα complex.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Dynamic nitrogen fixation in an aerobic endophyte of Populus

Biological nitrogen fixation by microbial diazotrophs can contribute significantly to nitrogen availability in non-nodulating plant species. In this study of molecular mechanisms and gene expression relating to biological nitrogen fixation, the aerobic nitrogen-fixing endophyte Burkholderia vietnamiensis, strain WPB, isolated from Populus trichocarpa served as a model for endophyte–poplar interactions. Nitrogen-fixing activity was observed to be dynamic on nitrogen-free medium with a subset of colonies growing to form robust, raised globular like structures. Secondary ion mass spectrometry (NanoSIMS) confirmed that N-fixation was uneven within the population. A fluorescent transcriptional reporter (GFP) revealed that the nitrogenase subunit nifH is not uniformly expressed across genetically identical colonies of WPB and that only ~11% of the population was actively expressing the nifH gene. Higher nifH gene expression was observed in clustered cells through monitoring individual bacterial cells using single-molecule fluorescence in situ hybridization. Through 15 N 2 enrichment, we identified key nitrogenous metabolites and proteins synthesized by WPB and employed targeted metabolomics in active and inactive populations. We cocultivated WPB Pnif-GFP with poplar within a RhizoChip, a synthetic soil habitat, which enabled direct imaging of microbial nifH expression within root epidermal cells. We observed that nifH expression is localized to the root elongation zone where the strain forms a unique physical interaction with the root cells. This work employed comprehensive experimentation to identify novel mechanisms regulating both biological nitrogen fixation and beneficial plant–endophyte interactions.

15 N-tracking metabolomics↗
DOE OSTI2022

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. In this work, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (~98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ~1.5 million atoms is discussed in the context of MuMMI.

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

Black Liquor Valorization by Using Marine Protist Thraustochytrium striatum and the Preliminary Metabolic Mechanism Study

Black liquor that contains various phenolic compounds from lignin solubilization has been the main byproduct of alkaline pretreatment of lignocellulosic biomass in the biorefinery. In this study, black liquor from alkaline pretreatment of corn stover was used as a sole carbon source for the cultivation of a marine protist, Thraustochytrium striatum. It was found that this strain can grow on black liquor and accumulate valuable products (e.g., fatty acids and carotenoids) simultaneously. Under optimal conditions (pH = 7 and NH 4 Cl = 2 g/L), the cell mass concentration reached 5.2 g/L with total aromatics decreased from 8.18 to 3.09 g/L within 7 day incubation. Fed-batch cultivation was adopted to increase the contents of total fatty acids and carotenoids to 13% and 0.24 mg/g dry cell mass, respectively. Although various compounds such as sugars and organic acids were detected in black liquor and consumed during microbial fermentation as carbon sources, lignin-derived compounds were identified as the major substrates for T. striatum fermentation. Of total aromatics consumed, monomers including p-coumaric acid, ferulic acid, vanillin, and syringaldehyde were observed to be consumed and converted while polymeric fragments were also depolymerized and degraded. Under optimal conditions, around 50% of the total aromatics was consumed. The dynamic changes of compounds in black liquor indicated that diverse metabolic processes were involved in black liquor degradation and utilization by T. striatum. A strong adaptation of T. striatum to a wide range of pH (3–9) was also observed during black liquor fermentation.

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

Ligand-induced transmembrane conformational coupling in monomeric $\mathrm{EGFR}$

Single pass cell surface receptors regulate cellular processes by transmitting ligand-encoded signals across the plasma membrane via changes to their extracellular and intracellular conformations. This transmembrane signaling is generally initiated by ligand binding to the receptors in their monomeric form. While subsequent receptor-receptor interactions are established as key aspects of transmembrane signaling, the contribution of monomeric receptors has been challenging to isolate due to the complexity and ligand-dependence of these interactions. By combining membrane nanodiscs produced with cell-free expression, single-molecule Förster Resonance Energy Transfer measurements, and molecular dynamics simulations, we report that ligand binding induces intracellular conformational changes within monomeric, full-length epidermal growth factor receptor (EGFR). Our observations establish the existence of extracellular/intracellular conformational coupling within a single receptor molecule. We implicate a series of electrostatic interactions in the conformational coupling and find the coupling is inhibited by targeted therapeutics and mutations that also inhibit phosphorylation in cells. Collectively, these results introduce a facile mechanism to link the extracellular and intracellular regions through the single transmembrane helix of monomeric EGFR, and raise the possibility that intramolecular transmembrane conformational changes upon ligand binding are common to single-pass membrane proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Enzymatic Nitrogen Insertion into Unactivated C–H Bonds

Selective functionalization of aliphatic C–H bonds, ubiquitous in molecular structures, could allow ready access to diverse chemical products. While enzymatic oxygenation of C–H bonds is well established, the analogous enzymatic nitrogen functionalization is still unknown; nature is reliant on pre-oxidized compounds for nitrogen incorporation. Likewise, synthetic methods for selective nitrogen derivatization of unbiased C–H bonds remain elusive. In this work, new-to-nature heme-containing nitrene transferases were used as starting points for the directed evolution of enzymes to selectively aminate and amidate unactivated C(sp3)–H sites. The desymmetrization of methyl- and ethylcyclohexane with divergent site selectivity is offered as demonstration. The evolved enzymes in these lineages are highly promiscuous and show activity towards a wide array of substrates, providing a foundation for further evolution of nitrene transferase function. Computational studies and kinetic isotope effects (KIEs) are consistent with a stepwise radical pathway involving an irreversible, enantiodetermining hydrogen atom transfer (HAT), followed by a lower-barrier diastereoselectivity determining radical rebound step. In-enzyme molecular dynamics (MD) simulations reveal a predominantly hydrophobic pocket with favorable dispersion interactions with the substrate. Here, by offering a direct path from saturated precursors, these enzymes present a new biochemical logic for accessing nitrogen-containing compounds.

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

Computational prediction of the effect of amino acid changes on the binding affinity between SARS-CoV-2 spike RBD and human ACE2

The association of the receptor binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein with human angiotensin-converting enzyme 2 (hACE2) represents the first required step for cellular entry. SARS-CoV-2 has continued to evolve with the emergence of several novel variants, and amino acid changes in the RBD have been implicated with increased fitness and potential for immune evasion. Reliably predicting the effect of amino acid changes on the ability of the RBD to interact more strongly with the hACE2 can help assess the implications for public health and the potential for spillover and adaptation into other animals. Here, we introduce a two-step framework that first relies on 48 independent 4-ns molecular dynamics (MD) trajectories of RBD-hACE2 variants to collect binding energy terms decomposed into Coulombic, covalent, van der Waals, lipophilic, generalized Born solvation, hydrogen bonding, π-π packing, and self-contact correction terms. The second step implements a neural network to classify and quantitatively predict binding affinity changes using the decomposed energy terms as descriptors. The computational base achieves a validation accuracy of 82.8% for classifying single–amino acid substitution variants of the RBD as worsening or improving binding affinity for hACE2 and a correlation coefficient of 0.73 between predicted and experimentally calculated changes in binding affinities. Both metrics are calculated using a fivefold cross-validation test. Our method thus sets up a framework for screening binding affinity changes caused by unknown single– and multiple–amino acid changes offering a valuable tool to predict host adaptation of SARS-CoV-2 variants toward tighter hACE2 binding.

60 APPLIED LIFE SCIENCES↗