Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural variant”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A self-consistent Hamiltonian model of the ponderomotive force and its structure preserving discretization

In the presence of an inhomogeneous oscillatory electric field, charged particles experience a net force, averaged over the oscillatory timescale, known as the ponderomotive force. We derive a one-dimensional Hamiltonian model which self-consistently couples the electromagnetic field to a plasma which experiences the ponderomotive force. We derive a family of structure preserving discretizations of the model of varying order in space and time using conforming and broken finite element exterior calculus spectral element methods. In all variants of our discretization framework, the method is found to conserve the Casimir invariants of the continuous model to machine precision and the energy to the order of the splitting method used.

Physics↗
DOE OSTI2024

Anomalous Retrograde Drifts in Obstructed Magnetron Microdischarges: a Consequence of a Field Reversal in the Anode Sheath?

The goal of this project was to carry out experiments at the Plasma Collaborative Research Facility (Princeton Plasma Physics Laboratory) aimed at confirming or refuting the hypothesis that an electric field reversal occurs in the near-anode region of a strongly obstructed planar magnetron microdischarges. Magnetron discharges are magnetized plasmas with E x B field topologies, and variants of such discharges are widely used in materials processing. The topologies are also similar to that of Hall thrusters used in space propulsion. Understanding the structure of these discharges will aid in developing reliable simulations that are intended to facilitate the further development of such technologies, as the field reversals, presumed to be due to strong axial-driven electron transport, provides further insight into transport properties. The possible existence of such a field reversal is hinted at because of the retrograde-propagating (i.e., negative E x B) spokes or gradient-driven drift instabilities seen in the high speed videos and segmented anode probing of the plasma. The presence of this field reversal should be reflected in the ion dynamics, particularly, in an ion drift towards the anode. The proposed experiments to be conducted at the PCRF are to characterize this anode-directed ion drift using state-of-the-art laser-induced fluorescence velocimetry. In particular, we seek to see if the near-anode ion velocity reverses direction when the external circuit is varied to alter the electron current demand imposed at the anode. Such an ion velocity reversal would confirm the mechanism, which is that the highly obstructed (small gapped) discharge produces strong axial gradients in plasma density, driving more electrons towards the anode then that demanded by the external circuit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Long-Range Allosteric Communication Modulated by Active Site Mn(II) Coordination Drives Catalysis in Xanthobacter autotrophicus Acetone Carboxylase

Acetone carboxylase (AC) from Xanthobacter autotrophicus is a 360 KDa α2β2γ2 heterohexamer that catalyzes the ATP-dependent formation of phosphorylated acetone and bicarbonate intermediates that react at Mn(II) metal active sites to form acetoacetate. Structural models of X. autotrophicus AC (XaAC) with and without nucleotides reveal that the binding and phosphorylation of the two substrates occurs ~40 Å from the Mn(II) active sites where acetoacetate is formed. Based on the crystal structures, a significant conformational change was proposed to open and close a tunnel that facilitates the passage of reaction intermediates between the sites for nucleotide binding and phosphorylation of substrates and Mn(II) sites of acetoacetate formation. We have employed electron paramagnetic resonance (EPR), kinetic assays, and hydrogen/deuterium exchange mass spectrometry (HDX-MS) of poised ligand-bound states and site-specific amino acid variants to complete an in-depth analysis of Mn(II) coordination and allosteric communication throughout the catalytic cycle. In contrast with the established paradigms for carboxylation, our analyses of XaAC suggested a carboxylate shift that couples both local and long-range structural transitions. Shifts in the coordination mode of a single carboxylic acid residue (αE89) mediate both catalysis proximal to a Mn(II) center and communication with an ATP active site in a separate subunit of a 180 kDa α2β2γ2 complex at a distance of 40 Å. This work demonstrates the power of combining structural models from X-ray crystallography with solution-phase spectroscopy and biophysical techniques to elucidate functional aspects of a multi-subunit enzyme.

Biochemistry & Molecular Biology↗
DOE OSTI2022

Structural definition of a pan-sarbecovirus neutralizing epitope on the spike S2 subunit

Three betacoronaviruses have crossed the species barrier and established human-to-human transmission causing significant morbidity and mortality in the past 20 years. The most current and widespread of these is SARS-CoV-2. The identification of CoVs with zoonotic potential in animal reservoirs suggests that additional outbreaks could occur. Monoclonal antibodies targeting conserved neutralizing epitopes on diverse CoVs can form the basis for prophylaxis and therapeutic treatments and enable the design of vaccines aimed at providing pan-CoV protection. We previously identified a neutralizing monoclonal antibody, CV3-25 that binds to the SARS-CoV-2 spike, neutralizes the SARS-CoV-2 Beta variant comparably to the ancestral Wuhan Hu-1 strain, cross neutralizes SARS-CoV-1 and binds to recombinant proteins derived from the spike-ectodomains of HCoV-OC43 and HCoV-HKU1. Here, we show that the neutralizing activity of CV3-25 is maintained against the Alpha, Delta, Gamma and Omicron variants of concern as well as a SARS-CoV-like bat coronavirus with zoonotic potential by binding to a conserved linear peptide in the stem-helix region. Negative stain electron microscopy and a 1.74 Å crystal structure of a CV3-25/peptide complex demonstrates that CV3-25 binds to the base of the stem helix at the HR2 boundary to an epitope that is distinct from other stem-helix directed neutralizing mAbs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Towards Quantum-Chemical Level Calculations of SARS-CoV-2 Spike Protein Variants of Concern by First Principles Density Functional Theory

The spike protein (S-protein) is a crucial part of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with its many domains responsible for binding, fusion, and host cell entry. In this review we use the density functional theory (DFT) calculations to analyze the atomic-scale interactions and investigate the consequences of mutations in S-protein domains. We specifically describe the key amino acids and functions of each domain, which are essential for structural stability as well as recognition and fusion processes with the host cell; in addition, we speculate on how mutations affect these properties. Such unprecedented large-scale ab initio calculations, with up to 5000 atoms in the system, are based on the novel concept of amino acid–amino acid-bond pair unit (AABPU) that allows for an alternative description of proteins, providing valuable information on partial charge, interatomic bonding and hydrogen bond (HB) formation. In general, our results show that the S-protein mutations for different variants foster an increased positive partial charge, alter the interatomic interactions, and disrupt the HB networks. We conclude by outlining a roadmap for future computational research of biomolecular virus-related systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Phosphoserine Charge State Drives Ion Condensation and Spatial Polyamine Presentation in Multirepeat Silaffin

Diatom silaffins direct silica biomineralization through heavily post-translationally modified repeat domains, yet how these modifications reshape the multirepeat conformational ensemble remains unknown. We report all-atom MD simulations of a 195-residue construct spanning repeats R1–R7 of Sil1p from Cylindrotheca fusiformis, carrying the full complement of native PTMs: phosphoserine (pSer), long-chain polyamines (LCPA), dimethyllysine (MLY), and trimethylhydroxylysine phosphate (TPL). We simulated three variants (Native, P1/singly deprotonated phosphate, and P2/doubly deprotonated phosphate) at two NaCl concentrations in triplicate for 500 ns each. All systems disorder from the AlphaFold 3 starting structure. Phosphate charge state, not ionic strength, is the dominant control of ensemble compaction and ion organization. Doubly deprotonated phosphate organizes an extensive Na + condensation shell (∼100 ions, 20% of box Na + ) and a heterogeneous bridging network that integrates both pSer and TPL phosphate groups. The resulting ensemble is compact with LCPA side chains exhibiting above-median solvent accessibility in 84% of simulation frames in P2 at 300 mM. This is higher than any other condition we simulated and indicates that polyamine groups are preferentially surface-presented in the most compact, ion-organized state. A charge-neutralization control confirms that this compact state is a structured intermediate maintained by the bridging network, not a simple collapsed globule. Here, this repeat-scale spatial organization is not captured by single-repeat peptide studies. Understanding the dynamics, mechanism, and spatial organization of PTM-rich silaffin at the repeat scale is a step closer to hierarchical biomimetic materials beyond simple silica morphologies.

Amines↗
DOE OSTI2025

Supramolecular Assembly of Lanthanide-Binding Tag Peptides for Aqueous Separation of Rare Earth Elements

Selective and eco-friendly separation and purification methods for rare earth elements (REEs) are necessary to meet the increasing demand for these valuable metals, which are extensively used in modern electronics and clean energy technologies. Mining feedstocks consist of REE mixtures as stable trivalent cations (Ln 3+ ) that are difficult to separate due to their identical charge and similar size. Lanthanide-binding tags (LBTs), peptide chelates that coordinate Ln 3+ in binding pockets, show promise as selective, high-affinity extractants. We demonstrate that the LBT variant LBTLLA 5– , designed for high selectivity for Tb 3+ , is an effective extractant, forming complexes with REEs in solution that subsequently organize into self-assembling structures rich in Ln 3+ . These structures condense into aggregates that can be separated, enabling an efficient, all-aqueous, eco-friendly separation process. The self-assembled structures are studied using dynamic light scattering, ζ-potential measurements, transmission electron microscopy, anomalous small-angle X-ray scattering, inductively coupled plasma optical emission spectroscopy, and ultraviolet–visible absorption spectroscopy, which confirm LBTLLA 5– peptide-REE ion binding and the further assembly of micron-scale structures rich in REEs. Molecular dynamics simulations reveal the interactions promoting aggregation as well as the integrity of the binding pocket upon self-assembly. We find that LBTLLA 5– :Ln 3+ complexes recruit excess cations within the macrostructures, and we demonstrate that aggregation and selective separation can be controlled by manipulating the metal-peptide ratio in solution. Furthermore, we demonstrate separation from equimolar mixtures of REE pairs Tb 3+ -Lu 3+ and Tb 3+ -La 3+ , supporting the application of LBT peptides as a platform for the selective separation of REEs.

LBT peptides↗
DOE OSTI2023

Field-theoretic analysis of hadronization using soft drop jet mass

One of the greatest challenges in quantum chromodynamics is understanding the hadronization mechanism, which is also crucial for carrying out precision physics with jet substructure. In this paper, we bring together recent advancements in our understanding of nonperturbative structure of the soft drop jet mass based on field theory, with precise perturbative calculations at next-to-next-to-leading logarithmic accuracy of its multidifferential variants. This allows for a model-independent analysis of power corrections associated with hadronization in a systematic manner. We test and calibrate hadronization models and their interplay with parton showers by comparing our universality predictions with various event generators for quark and gluon initiated jets in both lepton-lepton and hadron-hadron collisions. Our findings reveal that hadronization models perform better for quark jets relative to gluon jets. Our results provide a valuable toolbox for precision studies with the soft drop jet mass and pave the way for future analyses using real-world collider data. The stringent constraints derived in our framework are useful for improving the modeling of hadronization and its interplay with parton showers in next-generation event generators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Altered Cardiac Energetics and Mitochondrial Dysfunction in Hypertrophic Cardiomyopathy

Background: Hypertrophic cardiomyopathy (HCM) is a complex disease partly explained by the effects of individual gene variants on sarcomeric protein biomechanics. At the cellular level, HCM mutations most commonly enhance force production, leading to higher energy demands. Despite significant advances in elucidating sarcomeric structure–function relationships, there is still much to be learned about the mechanisms that link altered cardiac energetics to HCM phenotypes. Here, we test the hypothesis that changes in cardiac energetics represent a common pathophysiologic pathway in HCM. Methods: We performed a comprehensive multiomics profile of the molecular (transcripts, metabolites, and complex lipids), ultrastructural, and functional components of HCM energetics using myocardial samples from 27 HCM patients and 13 normal controls (donor hearts). Results: Integrated omics analysis revealed alterations in a wide array of biochemical pathways with major dysregulation in fatty acid metabolism, reduction of acylcarnitines, and accumulation of free fatty acids. HCM hearts showed evidence of global energetic decompensation manifested by a decrease in high energy phosphate metabolites (ATP, ADP, and phosphocreatine) and a reduction in mitochondrial genes involved in creatine kinase and ATP synthesis. Accompanying these metabolic derangements, electron microscopy showed an increased fraction of severely damaged mitochondria with reduced cristae density, coinciding with reduced citrate synthase activity and mitochondrial oxidative respiration. These mitochondrial abnormalities were associated with elevated reactive oxygen species and reduced antioxidant defenses. However, despite significant mitochondrial injury, HCM hearts failed to upregulate mitophagic clearance. Conclusions: Overall, our findings suggest that perturbed metabolic signaling and mitochondrial dysfunction are common pathogenic mechanisms in patients with HCM. These results highlight potential new drug targets for attenuation of the clinical disease through improving metabolic function and reducing mitochondrial injury.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Chickpea NCR13 disulfide cross-linking variants exhibit profound differences in antifungal activity and modes of action

Small cysteine-rich antifungal peptides with multi-site modes of action (MoA) have potential for development as biofungicides. In particular, legumes of the inverted repeat-lacking clade express a large family of nodule-specific cysteine-rich (NCR) peptides that orchestrate differentiation of nitrogen-fixing bacteria into bacteroids. These NCRs can form two or three intramolecular disulfide bonds and a subset of these peptides with high cationicity exhibits antifungal activity. However, the importance of intramolecular disulfide pairing and MoA against fungal pathogens for most of these plant peptides remains to be elucidated. Our study focused on a highly cationic chickpea NCR13, which has a net charge of +8 and contains six cysteines capable of forming three disulfide bonds. NCR13 expression in Pichia pastoris resulted in formation of two peptide folding variants, NCR13_PFV1 and NCR13_PFV2, that differed in the pairing of two out of three disulfide bonds despite having an identical amino acid sequence. The NMR structure of each PFV revealed a unique three-dimensional fold with the PFV1 structure being more compact but less dynamic. Surprisingly, PFV1 and PFV2 differed profoundly in the potency of antifungal activity against several fungal plant pathogens and their multi-faceted MoA. PFV1 showed significantly faster fungal cell-permeabilizing and cell entry capabilities as well as greater stability once inside the fungal cells. Additionally, PFV1 was more effective in binding fungal ribosomal RNA and inhibiting protein translation in vitro. Furthermore, when sprayed on pepper and tomato plants, PFV1 was more effective in reducing disease symptoms caused by Botrytis cinerea, causal agent of gray mold disease in fruits, vegetables, and flowers. In conclusion, our work highlights the significant impact of disulfide pairing on the antifungal activity and MoA of NCR13 and provides a structural framework for design of novel, potent antifungal peptides for agricultural use.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Resolving Lipid Structures for Threat-Agnostic Signatures

Development of threat-agnostic biosignatures is critical for the detection of and response to biological threats that go beyond the historical list-based approach. Lipids are a class of structurally diverse biomolecules that hold great promise as relevant threat-agnostic biosignatures, yet the fundamental challenge of resolving and decoding structural differences among lipids remains an obstacle in the pursuit of signature development. Structural analysis of complex lipids remains stymied due to spectral complexity and the inability to resolve low abundance lipids. To address these challenges, we integrated ozonolysis with ion mobility mass spectrometry (IMS-MS) analysis to determine the structural information for complex lipid species and (2) reduced lipid sample complexity prior to introduction into the IMS-MS using normal phase high pressure liquid chromatography (HPLC) to enable the ozonolysis of low abundance lipids. We applied these new capabilities to define the lipid structures in antimicrobial resistant (AMR) and antimicrobial susceptible (AMS) biothreat organism Yersinia pestis and cells infected with Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2) variants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Determining Compression Characteristics of Honeycomb Material - 19674

The objective is to determine compression test characteristics of stainless steel honeycomb material to be able to represent honeycomb structures accurately in analytical models used to simulate hypothetical accident scenarios of shipping packages. Honeycomb is a material used primary in the aerospace industry due to its high strength-to-weight ratio. Because it doesn't have a shelf life and can withstand high heats, it is an excellent candidate for a structural material in package designs. Honeycomb has orthotropic material properties. The honeycomb currently being evaluated is made of metal ribbons spot welded together to form hexagon pattern between metal plates. The hexagon ribbons are brazed to the metal plates. The complexity of the honeycomb significantly lengthens the simulation time, which is further complicated when the brazing, welding, and imperfections of the material is considered. This means to affectively represent the honeycomb in simulation programs, such as Abaqus, the material needs to be approximated as a uniform orthotropic material. This requires material properties in each of the three directions, T, W, and L. The T direction is defined as the 'strong' direction, perpendicular to the honeycomb sheet. The L direction is parallel the ribbon and the W direction is perpendicular to the ribbon. Honeycomb has 3 different compression stages. Stage 1 is the initial compression where the honeycomb maintains its structural integrity and does not permanently deform from loads from normal operations. Stage 2 is when the initial buckling causes irreversible damage to the honeycomb. Stage 2 is the one we are most interested in because it absorbs the most energy from a hypothetical accident scenario. During Stage 2, the honeycomb fails layer-by-layer, indicating that the more layers, the longer material crushing is sustained, and thus the more energy that is absorbed. Stage 3 is final compression, similar to compressing solid metal. During stage 3 the stress increases with a diminishing rate of elongation and the honeycomb is completely failed where the plates between the honeycomb core sandwich the crushed honeycomb ribbon. The stress-strain graph bellow shows compression test of 3 different samples in the T direction. The red lines divide the three stages. The far left is Stage 1, the middle stage 2, and the right stage 3. The data collected is displayed in the table below. The expressions in the table represent the regression formula that represents each stage of the compression on a stress-strain graph. Stage 1 is assumed to intersect with the origin; data for stage 3 for the crush in W and L directions where unable to be gathered due to the nature of failure for those directions. Stage 1 and 2 are linear regressions while stage 3 is a degree 2, to best match the curve. Linear regression is used on stage 2 because ultimately that represents the energy absorbed. With more data, a degree n x 2 regression would be a more appropriate regression for stage 2, where n is the number of layers. The table below is the approximation of each scenario and stage. Each stage starts and ends at the intersection of the next stage. There where two main objectives for these tests: firstly to understand how honeycomb performs under extreme compression, and secondly to be able to numerically represent the honeycomb structure. Both of these where completed to varying degrees. It is important understood how the honeycomb fails. If it is crushed in the T direction it retains its integrity even after being crushed, however when crushed in the L and W crush direction, if it fails, it disintegrates, and loses all integrity. Also, the more layers the more time the material spends in stage 2. Failure is started by buckling, thus if there is any imperfection, the stress will not spike but transition straight into stage 2. The data collected and aggregated can be used for initial simulation of honeycomb used in packages. The initial testing has set the ground work for more data to be collected in order to verify results and to allow more confidence in the simulation results. This is only the initial data collected. The next steps is to continue to collect more data to verify results and to test more variants of honeycomb, with different brazing, layers, and shape.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Rapid identification of neutralizing antibodies against SARS-CoV-2 variants by mRNA display

The increasing prevalence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants with the ability to escape existing humoral protection conferred by previous infection and/or immunization necessitates the discovery of broadly reactive neutralizing antibodies (nAbs). Utilizing mRNA display, we identify a set of antibodies against SARS-CoV-2 spike (S) proteins and characterize the structures of nAbs that recognize epitopes in the S1 subunit of the S glycoprotein. These structural studies reveal distinct binding modes for several antibodies, including the targeting of rare cryptic epitopes in the receptor-binding domain (RBD) of S that interact with angiotensin-converting enzyme 2 (ACE2) to initiate infection, as well as the S1 subdomain 1. Further, we engineer a potent ACE2-blocking nAb to sustain binding to S RBD with the E484K and L452R substitutions found in multiple SARS-CoV-2 variants. We demonstrate that mRNA display is an approach for the rapid identification of nAbs that can be used in combination to combat emerging SARS-CoV-2 variants.

SARS-CoV-2↗
DOE OSTI2022

Latent Mechanisms of Polarization Switching from In Situ Electron Microscopy Observations

In situ scanning transmission electron microscopy enables observation of the domain dynamics in ferroelectric materials as a function of externally applied bias and temperature. The resultant data sets contain a wealth of information on polarization switching and phase transition mechanisms. However, identification of these mechanisms from observational data sets has remained a problem due to a large variety of possible configurations, many of which are degenerate. Here, an approach based on a combination of deep learning-based semantic segmentation, rotationally invariant variational autoencoder (VAE), and non-negative matrix factorization to enable learning of a latent space representation of the data with multiple real-space rotationally equivalent variants mapped to the same latent space descriptors is introduced. By varying the size of training sub-images in the VAE, the degree of complexity in the structural descriptors is tuned from simple domain wall detection to the identification of switching pathways. Importantly, this yields a powerful tool for the exploration of the dynamic data in mesoscopic electron, scanning probe, optical, and chemical imaging. Moreover, this work adds to the growing body of knowledge of incorporating physical constraints into the machine and deep-learning methods to improve learned descriptors of physical phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Lithium Iron Aluminum Nickelate, LiNi x Fe y Al z O 2 —New Sustainable Cathodes for Next–Generation Cobalt–Free Li–Ion Batteries

In recent years, cobalt has become a critical constraint on the supply chain of the Li-ion battery industry. Here, with the ever-increasing projections for electric vehicles, the dependency of current Li-ion batteries on the ever-fluctuating cobalt prices poses serious environmental and sustainability issues. To address these challenges, a new class of cobalt-free materials with general formula of LiNi x Fe y Al z O 2 ( x + y + z = 1), termed as the lithium iron aluminum nickelate (NFA) class of cathodes, is introduced. These cobalt-free materials are synthesized using the sol–gel process to explore their compositional landscape by varying aluminum and iron. These NFA variants are characterized using electron microscopy, neutron and X-ray diffraction, and Mössbauer and X-ray photoelectron spectroscopy to investigate their morphological, physical, and crystal-structure properties. Operando experiments by X-ray diffraction, Mössbauer spectroscopy, and galvanostatic intermittent titration have been also used to study the crystallographic transitions, electrochemical activity, and Li-ion diffusivity upon lithium removal and uptake in the NFA cathodes. NFA compositions yield specific capacities of ≈200 mAh g -1 , demonstrating reasonable rate capability and cycling stability with ≈80% capacity retention after 100 charge/discharge cycles. While this is an early stage of research, the potential that these cathodes could have as viable candidates in next-generation cobalt-free lithium-ion batteries is highlighted here.

25 ENERGY STORAGE↗
DOE OSTI2023

Potent and selective covalent inhibition of the papain-like protease from SARS-CoV-2

Abstract Direct-acting antivirals are needed to combat coronavirus disease 2019 (COVID-19), which is caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The papain-like protease (PLpro) domain of Nsp3 from SARS-CoV-2 is essential for viral replication. In addition, PLpro dysregulates the host immune response by cleaving ubiquitin and interferon-stimulated gene 15 protein from host proteins. As a result, PLpro is a promising target for inhibition by small-molecule therapeutics. Here we design a series of covalent inhibitors by introducing a peptidomimetic linker and reactive electrophile onto analogs of the noncovalent PLpro inhibitor GRL0617. The most potent compound inhibits PLpro with k inact /K I = 9,600 M −1 s −1 , achieves sub-μM EC 50 values against three SARS-CoV-2 variants in mammalian cell lines, and does not inhibit a panel of human deubiquitinases (DUBs) at >30 μM concentrations of inhibitor. An X-ray co-crystal structure of the compound bound to PLpro validates our design strategy and establishes the molecular basis for covalent inhibition and selectivity against structurally similar human DUBs. These findings present an opportunity for further development of covalent PLpro inhibitors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Large-scale application of free energy perturbation calculations for antibody design

Abstract Alchemical free energy perturbation (FEP) is a rigorous and powerful technique to calculate the free energy difference between distinct chemical systems. Here we report our implementation of automated large-scale FEP calculations, using the Amber software package, to facilitate antibody design and evaluation. In combination with Hamiltonian replica exchange, our FEP simulations aim to predict the effect of mutations on both the binding affinity and the structural stability. Importantly, we incorporate multiple strategies to faithfully estimate the statistical uncertainties in the FEP results. As a case study, we apply our protocols to systematically evaluate variants of the m396 antibody for their conformational stability and their binding affinity to the spike proteins of SARS-CoV-1 and SARS-CoV-2. By properly adjusting relevant parameters, the particle collapse problems in the FEP simulations are avoided. Furthermore, large statistical errors in a small fraction of the FEP calculations are effectively reduced by extending the sampling, such that acceptable statistical uncertainties are achieved for the vast majority of the cases with a modest total computational cost. Finally, our predicted conformational stability for the m396 variants is qualitatively consistent with the experimentally measured melting temperatures. Our work thus demonstrates the applicability of FEP in computational antibody design.

59 BASIC BIOLOGICAL SCIENCES↗