Search NASA⌕ Search

SEARCH · Search NASA

Results for “Protein structure determination”

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

Legionella metaeffector MavL reverses ubiquitin ADP-ribosylation via a conserved arginine-specific macrodomain

ADP-ribosylation is a reversible post-translational modification involved in various cellular activities. Removal of ADP-ribosylation requires (ADP-ribosyl)hydrolases, with macrodomain enzymes being a major family in this category. The pathogen Legionella pneumophila mediates atypical ubiquitination of host targets using the SidE effector family in a process that involves ubiquitin ADP-ribosylation on arginine 42 as an obligatory step. Here, we show that the Legionella macrodomain effector MavL regulates this pathway by reversing the arginine ADP-ribosylation, likely to minimize potential detrimental effects caused by the modified ubiquitin. We determine the crystal structure of ADP-ribose-bound MavL, providing structural insights into recognition of the ADP-ribosyl group and catalytic mechanism of its removal. Further analyses reveal DUF4804 as a class of MavL-like macrodomain enzymes whose representative members show unique selectivity for mono-ADP-ribosylated arginine residue in synthetic substrates. We find such enzymes are also present in eukaryotes, as exemplified by two previously uncharacterized (ADP-ribosyl)hydrolases in Drosophila melanogaster . Crystal structures of several proteins in this class provide insights into arginine specificity and a shared mode of ADP-ribose interaction distinct from previously characterized macrodomains. Collectively, our study reveals a new regulatory layer of SidE-catalyzed ubiquitination and expands the current understanding of macrodomain enzymes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Structures of methane and ammonia monooxygenases in native membranes

Methane- and ammonia-oxidizing bacteria play key roles in the global carbon and nitrogen cycles, respectively. These bacteria use homologous copper membrane monooxygenases to accomplish the defining chemical transformations of their metabolisms: the oxidations of methane to methanol by particulate methane monooxygenase (pMMO) and ammonia to hydroxylamine by ammonia monooxygenase (AMO), enzymes of prime interest for applications in mitigating climate change. However, investigations of these enzymes have been hindered by the need for disruptive detergent solubilization prior to structure determination, confounding studies of pMMO and precluding studies of AMO. Here, we overcome these challenges by using cryoEM to visualize pMMO and AMO directly in their native membrane arrays at 2.4 to 2.8 Å resolution. These structures reveal details of the copper centers, numerous bound lipids, and previously unobserved components, including identifiable and distinct supernumerary helices interacting with pMMO and AMO, suggesting a widespread role for these helices in copper membrane monooxygenases. Comparisons between these structures, their metallocofactors, and their unexpected protein–protein interactions highlight features that may govern activity or the formation of higher-order arrays in native membranes. The ability to obtain molecular insights within the native membrane will enable further understanding of these environmentally important enzymes.

Science & Technology - Other Topics↗
DOE OSTI2020

A deubiquitylase with an unusually high-affinity ubiquitin-binding domain from the scrub typhus pathogen Orientia tsutsugamushi

Ubiquitin mediated signaling contributes critically to host cell defenses during pathogen infection. Many pathogens manipulate the ubiquitin system to evade these defenses. Here we characterize a likely effector protein bearing a deubiquitylase (DUB) domain from the obligate intracellular bacterium Orientia tsutsugamushi, the causative agent of scrub typhus. The Ulp1-like DUB prefers ubiquitin substrates over ubiquitin-like proteins and efficiently cleaves polyubiquitin chains of three or more ubiquitins. The co-crystal structure of the DUB (OtDUB) domain with ubiquitin revealed three bound ubiquitins: one engages the S1 site, the second binds an S2 site contributing to chain specificity and the third binds a unique ubiquitin-binding domain (UBD). The UBD modulates OtDUB activity, undergoes a pronounced structural transition upon binding ubiquitin, and binds monoubiquitin with an unprecedented ~5 nM dissociation constant. The characterization and high-resolution structure determination of this enzyme should aid in its development as a drug target to counter Orientia infections.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Site specific N- and O-glycosylation mapping of the spike proteins of SARS-CoV-2 variants of concern

Abstract The glycosylation on the spike (S) protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, modulates the viral infection by altering conformational dynamics, receptor interaction and host immune responses. Several variants of concern (VOCs) of SARS-CoV-2 have evolved during the pandemic, and crucial mutations on the S protein of the virus have led to increased transmissibility and immune escape. In this study, we compare the site-specific glycosylation and overall glycomic profiles of the wild type Wuhan-Hu-1 strain (WT) S protein and five VOCs of SARS-CoV-2: Alpha, Beta, Gamma, Delta and Omicron. Interestingly, both N- and O-glycosylation sites on the S protein are highly conserved among the spike mutant variants, particularly at the sites on the receptor-binding domain (RBD). The conservation of glycosylation sites is noteworthy, as over 2 million SARS-CoV-2 S protein sequences have been reported with various amino acid mutations. Our detailed profiling of the glycosylation at each of the individual sites of the S protein across the variants revealed intriguing possible association of glycosylation pattern on the variants and their previously reported infectivity. While the sites are conserved, we observed changes in the N- and O-glycosylation profile across the variants. The newly emerged variants, which showed higher resistance to neutralizing antibodies and vaccines, displayed a decrease in the overall abundance of complex-type glycans with both fucosylation and sialylation and an increase in the oligomannose-type glycans across the sites. Among the variants, the glycosylation sites with significant changes in glycan profile were observed at both the N -terminal domain and RBD of S protein, with Omicron showing the highest deviation. The increase in oligomannose-type happens sequentially from Alpha through Delta. Interestingly, Omicron does not contain more oligomannose-type glycans compared to Delta but does contain more compared to the WT and other VOCs. O-glycosylation at the RBD showed lower occupancy in the VOCs in comparison to the WT. Our study on the sites and pattern of glycosylation on the SARS-CoV-2 S proteins across the VOCs may help to understand how the virus evolved to trick the host immune system. Our study also highlights how the SARS-CoV-2 virus has conserved both N - and O - glycosylation sites on the S protein of the most successful variants even after undergoing extensive mutations, suggesting a correlation between infectivity/ transmissibility and glycosylation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Structural and biophysical insights into targeting of claudin-4 by a synthetic antibody fragment

Claudins are a 27-member family of ~25 kDa membrane proteins that integrate into tight junctions to form molecular barriers at the paracellular spaces between endothelial and epithelial cells. As the backbone of tight junction structure and function, claudins are attractive targets for modulating tissue permeability to deliver drugs or treat disease. However, structures of claudins are limited due to their small sizes and physicochemical properties—these traits also make therapy development a challenge. Here we report the development of a synthetic antibody fragment (sFab) that binds human claudin-4 and the determination of a high-resolution structure of it bound to claudin-4/enterotoxin complexes using cryogenic electron microscopy. Structural and biophysical results reveal this sFabs mechanism of select binding to human claudin-4 over other homologous claudins and establish the ability of sFabs to bind hard-to-target claudins to probe tight junction structure and function. The findings provide a framework for tight junction modulation by sFabs for tissue-selective therapies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Csx3 is a cyclic oligonucleotide phosphodiesterase associated with type III CRISPR–Cas that degrades the second messenger cA 4

Cas10 is the signature gene for type III CRISPR–Cas surveillance complexes. Unlike type I and type II systems, type III systems do not require a protospacer adjacent motif and target nascent RNA associated with transcriptionally active DNA. Further, target RNA recognition activates the cyclase domain of Cas10, resulting in the synthesis of cyclic oligoadenylate second messengers. These second messengers are recognized by ancillary Cas proteins harboring CRISPR-associated Rossmann fold (CARF) domains and regulate the activities of these proteins in response to invading nucleic acid. Csx3 is a distant member of the CARF domain superfamily previously characterized as a Mn 2+ -dependent deadenylation exoribonuclease. However, its specific role in CRISPR–Cas defense remains to be determined. Here we show that Csx3 is strongly associated with type III systems and that Csx3 binds cyclic tetra-adenylate (cA 4 ) second messenger with high affinity. Further, Csx3 harbors cyclic oligonucleotide phosphodiesterase activity that quickly degrades this cA 4 signal. In addition, structural analysis identifies core elements that define the CARF domain fold, and the mechanistic basis for ring nuclease activity is discussed. Overall, the work suggests that Csx3 functions within CRISPR–Cas as a counterbalance to Cas10 to regulate the duration and amplitude of the cA 4 signal, providing an off ramp from the programmed cell death pathway in cells that successfully cure viral infection.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Cryo-EM structures of engineered active bc 1- cbb 3 type CIII 2 CIV super-complexes and electronic communication between the complexes

Respiratory electron transport complexes are organized as individual entities or combined as large supercomplexes (SC). Gram-negative bacteria deploy a mitochondrial-like cytochrome (cyt) bc 1 (Complex III, CIII 2 ), and may have specific cbb 3 -type cyt c oxidases (Complex IV, CIV) instead of the canonical aa 3 -type CIV. Electron transfer between these complexes is mediated by soluble ( c 2 ) and membrane-anchored ( c y ) cyts. Here, we report the structure of an engineered bc 1 - cbb 3 type SC (CIII 2 CIV, 5.2Å resolution) and three conformers of native CIII 2 (3.3 Å resolution). The SC is active in vivo and in vitro, contains all catalytic subunits and cofactors, and two extra transmembrane helices attributed to cyt c y and the assembly factor CcoH. The cyt c y is integral to SC, its cyt domain is mobile and it conveys electrons to CIV differently than cyt c 2 . The successful production of a native-like functional SC and determination of its structure illustrate the characteristics of membrane-confined and membrane-external respiratory electron transport pathways in Gram-negative bacteria.

42 ENGINEERING↗
DOE OSTI2023

Unusual modifications of protein biomarkers expressed by plasmid, prophage, and bacterial host of pathogenic Escherichia coli identified using top‐down proteomic analysis

Rationale Pathogenic bacteria often carry prophage (bacterial viruses) and plasmids (small circular pieces of DNA) that may harbor toxin, antibacterial, and antibiotic resistance genes. Proteomic characterization of pathogenic bacteria should include the identification of host proteins and proteins produced by prophage and plasmid genomes. Methods Protein biomarkers of two strains of Shiga toxin–producingEscherichia coli(STEC) were identified using antibiotic induction, matrix‐assisted laser desorption/ionization tandem time‐of‐flight (MALDI‐TOF‐TOF) tandem mass spectrometry (MS/MS) with post‐source decay (PSD), top‐down proteomic (TDP) analysis, and plasmid sequencing. Alphafold2 was also used to compare predicted in silico structures of the identified proteins to prominent fragment ions generated using MS/MS‐PSD. Strain samples were also analyzed with and without chemical reduction treatment to detect the attachment of pendant groups bound by thioester or disulfide bonds. Results Shiga toxin was detected and/or identified in both STEC strains. For the first time, we also identified the osmotically inducible protein (OsmY) whose sequence unexpectedly had two forms: a full and a truncated sequence. The truncated OsmY terminates in the middle of an α‐helix as determined by Alphafold2. A plasmid‐encoded colicin immunity protein was also identified with and without attachment of an unidentified cysteine‐bound pendant group (~307 Da). Plasmid sequencing confirmed top‐down analysis and the identification of a promoter upstream of the immunity gene that is activated by antibiotic induction, that is, SOS box. Conclusions TDP analysis, coupled with other techniques (e.g., antibiotic induction, chemical reduction, plasmid sequencing, and in silico protein modeling), is a powerful tool to identify proteins (and their modifications), including prophage‐ and plasmid‐encoded proteins, produced by pathogenic microorganisms.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Biomolecular condensates form spatially inhomogeneous network fluids

The functions of biomolecular condensates are thought to be influenced by their material properties, and these will be determined by the internal organization of molecules within condensates. However, structural characterizations of condensates are challenging, and rarely reported. Here, we deploy a combination of small angle neutron scattering, fluorescence recovery after photobleaching, and coarse-grained molecular dynamics simulations to provide structural descriptions of model condensates that are formed by macromolecules from nucleolar granular components (GCs). We show that these minimal facsimiles of GCs form condensates that are network fluids featuring spatial inhomogeneities across different length scales that reflect the contributions of distinct protein and peptide domains. The network-like inhomogeneous organization is characterized by a coexistence of liquid- and gas-like macromolecular densities that engenders bimodality of internal molecular dynamics. These insights suggest that condensates formed by multivalent proteins share features with network fluids formed by systems such as patchy or hairy colloids.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Mechanistic insights into a heterobifunctional degrader-induced PTPN2/N1 complex

PTPN2 (protein tyrosine phosphatase non-receptor type 2, or TC-PTP) and PTPN1 are attractive immuno-oncology targets, with the deletion of Ptpn1 and Ptpn2 improving response to immunotherapy in disease models. Targeted protein degradation has emerged as a promising approach to drug challenging targets including phosphatases. We developed potent PTPN2/N1 dual heterobifunctional degraders (Cmpd-1 and Cmpd-2) which facilitate efficient complex assembly with E3 ubiquitin ligase CRL4 CRBN , and mediate potent PTPN2/N1 degradation in cells and mice. To provide mechanistic insights into the cooperative complex formation introduced by degraders, we employed a combination of structural approaches. Our crystal structure reveals how PTPN2 is recognized by the tri-substituted thiophene moiety of the degrader. We further determined a high-resolution structure of DDB1-CRBN/Cmpd-1/PTPN2 using single-particle cryo-electron microscopy (cryo-EM). This structure reveals that the degrader induces proximity between CRBN and PTPN2, albeit the large conformational heterogeneity of this ternary complex. The molecular dynamic (MD)-simulations constructed based on the cryo-EM structure exhibited a large rigid body movement of PTPN2 and illustrated the dynamic interactions between PTPN2 and CRBN. Together, our study demonstrates the development of PTPN2/N1 heterobifunctional degraders with potential applications in cancer immunotherapy. Furthermore, the developed structural workflow could help to understand the dynamic nature of degrader-induced cooperative ternary complexes.

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

In situ potassium and hydrogen ion exchange into a cubic zirconium silicate microporous material

The selective separation of ions from aqueous systems, and even in the human body, is a crucial to overall environmental management and health. Nanoporous materials are widely known for their selective removal of cations from aqueous media, and therefore have been targeted for use as a pharmaceutical to treat hyperkalemia. This study investigated the detailed crystallographic molecular mechanisms that control the potassium ion selectivity in the nanoporous cubic zirconium silicate (CZS) related materials. Using time-resolved in situ Raman spectroscopy and time-resolved in situ X-ray diffraction, the selectivity mechanisms were determined to involve a synchronous cation-cation repulsion process that serves to open a favorable coordination bonding environment for potassium, not unlike the ion selectivity filter process found in potassium ion channels in proteins. Enhancement of ion exchange was observed when the CZS material was in a partial protonated state (≈3:1 Na:H), causing an expansion of the unit-cell volume, enlargement of the 7 member-ring window, and distortion of framework polyhedra, which allowed increased accessibility to the cage structures and resulted in rapid irreversible potassium ion exchange.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

RNA polymerase II trapped on a molecular treadmill: Structural basis of persistent transcriptional arrest by a minor groove DNA binder

Elongating RNA polymerase II (Pol II) can be paused or arrested by a variety of obstacles. These obstacles include DNA lesions, DNA-binding proteins, and small molecules. Hairpin pyrrole-imidazole (Py-Im) polyamides bind to the minor groove of DNA in a sequence-specific manner and induce strong transcriptional arrest. Remarkably, this Py-Im–induced Pol II transcriptional arrest is persistent and cannot be rescued by transcription factor TFIIS. In contrast, TFIIS can effectively rescue the transcriptional arrest induced by a nucleosome barrier. The structural basis of Py-Im–induced transcriptional arrest and why TFIIS cannot rescue this arrest remain elusive. Here we determined the X-ray crystal structures of four distinct Pol II elongation complexes (Pol II ECs) in complex with hairpin Py-Im polyamides as well as of the hairpin Py-Im polyamides–dsDNA complex. We observed that the Py-Im oligomer directly interacts with RNA Pol II residues, introduces compression of the downstream DNA duplex, prevents Pol II forward translocation, and induces Pol II backtracking. These results, together with biochemical studies, provide structural insight into the molecular mechanism by which Py-Im blocks transcription. Our structural study reveals why TFIIS fails to promote Pol II bypass of Py-Im–induced transcriptional arrest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Thermally stable manganese(Ⅲ) peroxido complexes with hindered N3 tripodal ligands: Structures and their physicochemical properties

Mononuclear manganese(III) peroxido complexes are candidates for the reaction intermediates in manganese containing proteins, such as manganese superoxide dismutase (Mn-SOD) etc. In this study, manganese(III) peroxido complexes [Mn(O 2 )(L3)] and [Mn(O 2 )(L10)] ligated by anionic N3 type ligands with sterically hindered substituents, hydrotris(3-tertiary butyl-5-isopropyl-1-pyrazolyl)borate (L3 - ) and hydrotris(3-adamantyl-5-isopropyl-1-pyrazolyl)borate (L10 - ), respectively, were structurally characterized. These complexes are the first examples of structurally characterized five-coordinate manganese(III) peroxido complexes. Their characteristic ν(O—O) and ν(Mn—O) stretchings were determined by using H 2 18 O 2 for the first time. Theoretical calculations were performed to obtain further insight into their structural parameters. The decomposed products were obtained as [{Mn III (μ-O)(L3)} 2 Mn IV ] and [Mn III (OH){L10(O)}] from [Mn(O 2 )(L3)] and [Mn(O 2 )(L10)], respectively.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Organic matter distribution in the icy environments of Taylor Valley, Antarctica

Glaciers can accumulate and release organic matter affecting the structure and function of associated terrestrial and aquatic ecosystems. Here we analyzed 18 ice cores collected from six locations in Taylor Valley (McMurdo Dry Valleys), Antarctica to determine the spatial abundance and quality of organic matter, and the spatial distribution of bacterial density and community structure from the terminus of the Taylor Glacier to the coast (McMurdo Sound). Our results showed that dissolved and particulate organic carbon (DOC and POC) concentrations in the ice core samples increased from the Taylor Glacier to McMurdo Sound, a pattern also shown by bacterial cell density. Fluorescence Excitation Emission Matrices Spectroscopy (EEMs) and multivariate parallel factor (PARAFAC) modeling identified one humic-like (C1) and one protein-like (C2) component in ice cores whose fluorescent intensities all increased from the Polar Plateau to the coast. The fluorescence index showed that the bioavailability of dissolved organic matter (DOM) also decreased from the Polar Plateau to the coast. Partial least squares path modeling analysis revealed that bacterial abundance was the main positive biotic factor influencing both the quantity and quality of organic matter. Marine aerosol influenced the spatial distribution of DOC more than katabatic winds in the ice cores. Certain bacterial taxa showed significant correlations with DOC and POC concentrations. Collectively, our results show the tight connectivity among organic matter spatial distribution, bacterial abundance and meteorology in the McMurdo Dry Valley ecosystem.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Plant-specific features of respiratory supercomplex I + III2 from Vigna radiata

Abstract The last steps of cellular respiration—an essential metabolic process in plants—are carried out by mitochondrial oxidative phosphorylation. This process involves a chain of multi-subunit membrane protein complexes (complexes I–V) that form higher-order assemblies called supercomplexes. Although supercomplexes are the most physiologically relevant form of the oxidative phosphorylation complexes, their functions and structures remain mostly unknown. Here we present the cryogenic electron microscopy structure of the supercomplex I + III 2 from Vigna radiata (mung bean). The structure contains the full subunit complement of complex I, including a newly assigned, plant-specific subunit. It also shows differences in the mitochondrial processing peptidase domain of complex III 2 relative to a previously determined supercomplex with complex IV. The supercomplex interface, while reminiscent of that in other organisms, is plant specific, with a major interface involving complex III 2 ’s mitochondrial processing peptidase domain and no participation of complex I’s bridge domain. The complex I structure suggests that the bridge domain sets the angle between the enzyme’s two arms, limiting large-scale conformational changes. Moreover, complex I’s catalytic loops and its response in active-to-deactive assays suggest that, in V. radiata , the resting complex adopts a non-canonical state and can sample deactive- or open-like conformations even in the presence of substrate. This study widens our understanding of the possible conformations and behaviour of complex I and supercomplex I + III 2 . Further studies of complex I and its supercomplexes in diverse organisms are needed to determine the universal and clade-specific mechanisms of respiration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Structure Function Studies of Photosystem II Using X-Ray Free Electron Lasers

The structure and mechanism of the water-oxidation chemistry that occurs in photosystem II have been subjects of great interest. The advent of X-ray free electron lasers allowed the determination of structures of the stable intermediate states and of steps in the transitions between these intermediate states, bringing a new perspective to this field. The room-temperature structures collected as the photosynthetic water oxidation reaction proceeds in real time have provided important novel insights into the structural changes and the mechanism of the water oxidation reaction. The time-resolved measurements have also given us a view of how this reaction—which involves multielectron, multiproton processes—is facilitated by the interaction of the ligands and the protein residues in the oxygen-evolving complex. These structures have also provided a picture of the dynamics occurring in the channels within photosystem II that are involved in the transport of the substrate water to the catalytic center and protons to the bulk.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

The catalytic domains of Streptococcus mutans glucosyltransferases: a structural analysis

Streptococcus mutans, found in the human oral cavity, is a significant contributor to the pathogenesis of dental caries. This bacterium expresses three genetically distinct types of glucosyltransferases named GtfB (GTF-I), GtfC (GTF-SI) and GtfD (GTF-S) that play critical roles in the development of dental plaque. The catalytic domains of GtfB, GtfC and GtfD contain conserved active-site residues for the overall enzymatic activity that relate to hydrolytic glycosidic cleavage of sucrose to glucose and fructose, release of fructose and generation of a glycosyl-enzyme intermediate in the reducing end. In a subsequent transglycosylation step, the glucosyl moiety is transferred to the nonreducing end of an acceptor to form a growing glucan polymer chain made up of glucose molecules. It has been proposed that both sucrose breakdown and glucan synthesis occur in the same active site of the catalytic domain, although the active site does not appear to be large enough to accommodate both functions. These three enzymes belong to glycoside hydrolase family 70 (GH70), which shows homology to glycoside hydrolase family 13 (GH13). GtfC synthesizes both soluble and insoluble glucans (α-1,3 and α-1,6 glycosidic linkages), while GtfB and GtfD synthesize only insoluble or soluble glucans, respectively. Here, crystal structures of the catalytic domains of GtfB and GtfD are reported. These structures are compared with previously determined structures of the catalytic domain of GtfC. With this work, apo structures and inhibitor-complex structures with acarbose are now available for the catalytic domains of GtfC and GtfB. The structure of GtfC with maltose allows further identification and comparison of active-site residues. A model of sucrose binding to GtfB is also included. The new structure of the catalytic domain of GtfD affords a structural comparison of the three S. mutans glycosyltransferases. Unfortunately, the catalytic domain of GtfD is not complete since crystallization resulted in the structure of a truncated protein lacking approximately 200 N-terminal residues of domain IV.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Assessment of software methods for estimating protein-protein relative binding affinities

A growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < -0.5 kcal/mol) with high (83–98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.

59 BASIC BIOLOGICAL SCIENCES↗