Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “binding affinities”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Structure–Function Relationships in Sequence-Controlled Copolymers for Rare Earth Element Chelation

The ability to tune material function through primary sequence is a defining feature of biological macromolecules, allowing precise control over structure and target interactions in complex aqueous environments. However, translating sequence–structure–function relationships to synthetic macromolecules is challenging due to their dispersity in sequence, conformation, and composition. Here, we report systematic studies of amphiphilic polymer chelators designed to probe how composition and patterning influence binding affinity and selectivity for rare earth elements (REEs), a series of technologically relevant metals with challenging separation profiles. A library of copolymers varying hydrophobic monomer composition and patterning was synthesized via reversible addition–fragmentation chain transfer (RAFT) polymerization, spanning statistical, gradient, and block architectures. REE binding was quantified using a high-throughput colorimetric assay, and reconstruction of polymer ensembles using kinetic stochastic simulations enabled quantitative comparisons of sequence heterogeneity, linking local monomer colocalization to emergent REE binding. Further, we investigated the role of different hydrophobic comonomers in tuning metal coordination, with binding trends linked to structural features that influence binding site desolvation. Complementary dynamic light scattering (DLS) and small-angle X-ray scattering (SAXS) measurements showed that both polymer and monomer architecture modulate metal-induced conformational changes, and that multichain assembly behavior emerges beyond critical hydrophobic thresholds. Sequence control also altered REE selectivity, with nonmonotonic differences observed across compositionally identical polymers with different sequence architectures. Together, these findings establish design principles that connect polymer sequence and structure to binding performance, guiding the design of macromolecular chelators with enhanced affinity and selectivity for applications in separations, sensing, and catalysis.

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

The effect of protein mutations on drug binding suggests ensuing personalised drug selection

The advent of personalised medicine promises a deeper understanding of mechanisms and therefore therapies. However, the connection between genomic sequences and clinical treatments is often unclear. We studied 50 breast cancer patients belonging to a population-cohort in the state of Qatar. From Sanger sequencing, we identified several new deleterious mutations in the estrogen receptor 1 gene (ESR1). The effect of these mutations on drug treatment in the protein target encoded by ESR1, namely the estrogen receptor, was achieved via rapid and accurate protein–ligand binding affinity interaction studies which were performed for the selected drugs and the natural ligand estrogen. Four nonsynonymous mutations in the ligand-binding domain were subjected to molecular dynamics simulation using absolute and relative binding free energy methods, leading to the ranking of the efficacy of six selected drugs for patients with the mutations. Our study shows that a personalised clinical decision system can be created by integrating an individual patient’s genomic data at the molecular level within a computational pipeline which ranks the efficacy of binding of particular drugs to variant proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Computing the Relative Affinity of Chlorophylls a and b to Light-Harvesting Complex II

In plants and algae, the primary antenna protein bound to photosystem II is light-harvesting complex II (LHCII), a pigment–protein complex that binds eight chlorophyll (Chl) a molecules and six Chl b molecules. Chl a and Chl b differ only in that Chl a has a methyl group (–CH 3 ) on one of its pyrrole rings, while Chl b has a formyl group (–CHO) at that position. This blue-shifts the Chl b absorbance relative to Chl a . It is not known how the protein selectively binds the right Chl type at each site. Knowing the selection criteria would allow the design of light-harvesting complexes that bind different Chl types, modifying an organism to utilize the light of different wavelengths. The difference in the binding affinity of Chl a and Chl b in pea and spinach LHCII was calculated using multiconformation continuum electrostatics and free energy perturbation. Both methods have identified some Chl sites where the bound Chl type ( a or b ) has a significantly higher affinity, especially when the protein provides a hydrogen bond for the Chl b formyl group. However, the Chl a sites often have little calculated preference for one Chl type, so they are predicted to bind a mixture of Chl a and b . The electron density of the spinach LHCII was reanalyzed, which, however, confirmed that there is negligible Chl b in the Chl a -binding sites. Finally, it is suggested that the protein chooses the correct Chl type during folding, segregating the preferred Chl to the correct binding site.

chemical calculations↗
DOE OSTI2025

Cross-species analysis of FcγRIIa/b (CD32a/b) polymorphisms at position 131: structural and functional insights into the mechanism of IgG- mediated phagocytosis in human and macaque

Introduction Antibodies play a critical role in immunity in part by mediating clearance of pathogens and infected cells by antibody-dependent cellular phagocytosis (ADCP) through engagement of Fc gamma receptors (FcγRs) on innate immune cells. Among these, FcγRIIa (CD32a) is a key activating receptor expressed on macrophages, dendritic cells, and other antigen-presenting cells. Its affinity for IgG and ability to mediate ADCP is influenced by allelic polymorphisms. In humans, a single amino acid polymorphism at position 131, where histidine (H) is substituted with arginine (R), leads to decreased IgG1 and IgG2 subclass binding affinity and, consequently, lower efficiency of phagocytic responses. Rhesus macaques ( Macaca mulatta ), which are widely used as nonhuman primate models, exhibit a similar polymorphism at position 131 of FcγRIIa, but with arginine replaced by proline (P). Here, we investigated structure-function relationships associated with the FcγRIIa polymorphism at position 131 in both species, specifically with respect to IgG1 and IgG2. Methods We determined the structures of complexes formed by each variant with IgG1 Fc and those formed by the higher affinity variant with IgG2 Fc for both species by x-ray crystallography and linked these structures to affinity and activity using SPR and an ADCP assay. We also determined the structure of human inhibitory FcγRIIb (CD32b) in complex with IgG1 Fc by x-ray crystallography. Results Through analysis of these structures, our studies reveal that FcγRIIa engagement is minimally influenced by Fc glycan composition, distinguishing it from FcγRIIIa whose affinity is strongly influenced by glycan-composition. Comparative structures of human and macaque FcγRIIa variants demonstrate species- and allele-specific differences in Fc binding, but our functional assays showed only minimal allele-specific effects in humans. In contrast, allele-specific effects in macaques were highly significant; the macaque P 131 variant showing uniformly reduced IgG affinity. Conclusion These insights highlight fundamental interspecies and allelic distinctions that are critical for interpreting FcγRIIa-mediated effector functions in macaque models and for optimizing translational antibody and vaccine design.

Tolbert, William D.↗
DOE OSTI2022

An Aluminum-Based Metal–Organic Cage for Cesium Capture

Metal–organic cages are a class of supramolecular structures that often require the careful selection of organic linkers and metal nodes. Of this class, few examples of metal–organic cages exist where the nodes are composed of main group metals. Herein, we have prepared an aluminum-based metal–organic cage, H 8 [Al 8 (pdc) 8 (OAc) 8 O 4 ] (Al-pdc-AA), using inexpensive and commercially available materials. The cage formation was achieved via solvothermal self-assembly of solvated aluminum and pyridine-dicarboxylic linkers in the presence of a capping agent, acetic acid. The obtained supramolecular structure was characterized by single-crystal X-ray diffraction (SCXRD), thermogravimetric analysis, and NMR spectroscopy. Based on crystal structure and computational analyses, the cage has a 3.7 Å diameter electron-rich cavity suitable for the binding of cations such as cesium (ionic radius of 1.69 Å). Here, the host–guest interactions were probed with 1 H and 133 Cs NMR spectroscopy in DMSO, where at low concentrations, Cs + binds to Al-pdc-AA in a 1:1 ratio. The binding site was identified from the crystal structure of CsH 7 [Al 8 (pdc) 8 (OAc) 8 O 4 ] (Cs + Al-pdc-AA), and a binding affinity of ~10 6 –10 7 M –1 was determined from NMR titration experiments. The Al-pdc-AA showed improved selectivity for cesium binding over alkali metal cations (Cs + > Rb + > K + >> Na + ~ Li + ). Collectively, the study reports a novel aluminum cage that can serve as a promising host for efficient and selective cesium removal.

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

Effect of Ionic Strength and Anion Identity on Rare‐Earth Element Chelation by Poly(acrylic acid)

Abstract Rare‐earth elements (REEs) are critical components of many modern technologies, but their efficient separation and purification continue to pose significant challenges. Metal‐binding polymers may enable the development of promising alternative approaches to conventional separation techniques (e.g., solvent extraction), and insights into the interactions between the metal ions and polymers in solution are crucial to such development. In this work, the effect of ionic strength on the chelation of REEs by poly(acrylic acid) is investigated using isothermal titration calorimetry (ITC) to measure the metal‐binding thermodynamics. ITC experiments revealed that when ionic strength increases from increasing the sodium acetate buffer strength, polymer‐REE binding becomes increasingly disfavored. However, when ionic strength is increased by the addition of sodium chloride, binding again becomes increasingly disfavored, but to a lesser extent than when ionic strength increases from increasing the buffer strength. Additionally, the reduction in binding affinity is enthalpically driven when ionic strength is raised with sodium chloride, while it is entropically driven when raised with sodium acetate buffer, suggesting competive binding by the acetate ions. These results highlight the importance of ionic strength in polymer‐metal binding thermodynamics and will guide the development of new processes for the purification of critical metals.

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

Integrating functional scoring and regulatory data to predict the effect of non-coding SNPs in a complex neurological disease

Abstract Most SNPs associated with complex diseases seem to lie in non-coding regions of the genome; however, their contribution to gene expression and disease phenotype remains poorly understood. Here, we established a workflow to provide assistance in prioritising the functional relevance of non-coding SNPs of candidate genes as susceptibility loci in polygenic neurological disorders. To illustrate the applicability of our workflow, we considered the multifactorial disorder migraine as a model to follow our step-by-step approach. We annotated the overlap of selected SNPs with regulatory elements and assessed their potential impact on gene expression based on publicly available prediction algorithms and functional genomics information. Some migraine risk loci have been hypothesised to reside in non-coding regions and to be implicated in the neurotransmission pathway. In this study, we used a set of 22 non-coding SNPs from neurotransmission and synaptic machinery-related genes previously suggested to be involved in migraine susceptibility based on our candidate gene association studies. After prioritising these SNPs, we focused on non-reported ones that demonstrated high regulatory potential: (1) VAMP2_rs1150 (3′ UTR) was predicted as a target of hsa-mir-5010-3p miRNA, possibly disrupting its own gene expression; (2) STX1A_rs6951030 (proximal enhancer) may affect the binding affinity of zinc-finger transcription factors (namely ZNF423) and disturb TBL2 gene expression; and (3) SNAP25_rs2327264 (distal enhancer) expected to be in a binding site of ONECUT2 transcription factor. This study demonstrated the applicability of our practical workflow to facilitate the prioritisation of potentially relevant non-coding SNPs and predict their functional impact in multifactorial neurological diseases.

Felício, Daniela↗
DOE OSTI2025

Mechanistic studies of small molecule ligands selective to RNA single G bulges

Abstract Small-molecule RNA binders have emerged as an important pharmacological modality. A profound understanding of the ligand selectivity, binding mode, and influential factors governing ligand engagement with RNA targets is the foundation for rational ligand design. Here, we report a novel class of coumarin derivatives exhibiting selective binding affinity towards single G RNA bulges. Harnessing the computational power of all-atom Gaussian accelerated molecular dynamics simulations, we unveiled a rare minor groove binding mode of the ligand with a key interaction between the coumarin moiety and the G bulge. This predicted binding mode is consistent with results obtained from structure-activity relationship studies and transverse relaxation measurements by nuclear magnetic resonance spectroscopy. We further generated 444 molecular descriptors from 69 coumarin derivatives and identified key contributors to the binding events, such as charge state and planarity, by lasso (least absolute shrinkage and selection operator) regression. Our work deepened the understanding of RNA-small molecule interactions and integrated a new framework for the rational design of selective small-molecule RNA binders.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Combining computational modeling and experimental library screening to affinity-mature VEEV-neutralizing antibody F5

Engineered monoclonal antibodies have proven to be highly effective therapeutics in recent viral outbreaks. However, despite technical advancements, an ability to rapidly adapt or increase antibody affinity and by extension, therapeutic efficacy, has yet to be fully realized. We endeavored to stand-up such a pipeline using molecular modeling combined with experimental library screening to increase the affinity of F5, a monoclonal antibody with potent neutralizing activity against Venezuelan Equine Encephalitis Virus (VEEV), to recombinant VEEV (IAB) E1E2 antigen. We modeled the F5/E1E2 binding interface and generated predictions for mutations to improve binding using a Rosetta-based approach and dTERMen, an informatics approach. The modeling was complicated by the fact that a high-resolution structure of F5 is not available and the H3 loop of F5 exceeds the length for which current modeling approaches can determine a unique structure. A subset of the predicted mutations from both methods were incorporated into a phage display library of scFvs. This library and a library generated by error-prone PCR were screened for binding affinity to the recombinant antigen. Results from the screens identified favorable mutations which were incorporated into 12 human-IgG1 variants. The best variant, containing eight mutations, improved KD from 0.63 nM (parental) to 0.01 nM. While this did not improve neutralization or therapeutic potency of F5 against IAB, it did increase cross-reactivity to other closely related VEEV epizootic and enzootic strains, demonstrating the potential of this method to rapidly adapt existing therapeutics to emerging viral strains.

affinity-maturation↗
DOE OSTI2023

Designer molecules of the synaptic organizer MDGA1 reveal 3D conformational control of biological function

MDGAs (MAM domain-containing glycosylphosphatidylinositol anchors) are synaptic cell surface molecules that regulate the formation of trans-synaptic bridges between neurexins (NRXNs) and neuroligins (NLGNs), which promote synaptic development. Mutations in MDGAs are implicated in various neuropsychiatric diseases. MDGAs bind NLGNs in cis on the postsynaptic membrane and physically block NLGNs from binding to NRXNs. In crystal structures, the six immunoglobulin (Ig) and single fibronectin III domains of MDGA1 reveal a striking compact, triangular shape, both alone and in complex with NLGNs. Whether this unusual domain arrangement is required for biological function or other arrangements occur with different functional outcomes is unknown. Here, we show that WT MDGA1 can adopt both compact and extended 3D conformations that bind NLGN2. Designer mutants targeting strategic molecular elbows in MDGA1 alter the distribution of 3D conformations while leaving the binding affinity between soluble ectodomains of MDGA1 and NLGN2 intact. In contrast, in a cellular context, these mutants result in unique combinations of functional consequences, including altered binding to NLGN2, decreased capacity to conceal NLGN2 from NRXN1β, and/or suppressed NLGN2-mediated inhibitory presynaptic differentiation, despite the mutations being located far from the MDGA1-NLGN2 interaction site. Thus, the 3D conformation of the entire MDGA1 ectodomain appears critical for its function, and its NLGN-binding site on Ig1-Ig2 is not independent of the rest of the molecule. As a result, global 3D conformational changes to the MDGA1 ectodomain via strategic elbows may form a molecular mechanism to regulate MDGA1 action within the synaptic cleft.

MDGAs↗
DOE OSTI2020

Enhanced Binding of SARS-CoV-2 Spike Protein to Receptor by Distal Polybasic Cleavage Sites

The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein plays a crucial role in binding the human cell receptor ACE2 that is required for viral entry. Many studies have been conducted to target the structures of RBD–ACE2 binding and to design RBD-targeting vaccines and drugs. Nevertheless, mutations distal from the SARS-CoV-2 RBD also impact its transmissibility and antibody can target non-RBD regions, suggesting the incomplete role of the RBD region in the spike protein–ACE2 binding. Here, in order to elucidate distant binding mechanisms, we analyze complexes of ACE2 with the wild-type spike protein and with key mutants via large-scale all-atom explicit solvent molecular dynamics simulations. We find that though distributed approximately 10 nm away from the RBD, the SARS-CoV-2 polybasic cleavage sites enhance, via electrostatic interactions and hydration, the RBD–ACE2 binding affinity. A negatively charged tetrapeptide (GluGluLeuGlu) is then designed to neutralize the positively charged arginine on the polybasic cleavage sites. We find that the tetrapeptide GluGluLeuGlu binds to one of the three polybasic cleavage sites of the SARS-CoV-2 spike protein lessening by 34% the RBD–ACE2 binding strength. This significant binding energy reduction demonstrates the feasibility to neutralize RBD–ACE2 binding by targeting this specific polybasic cleavage site. Our work enhances understanding of the binding mechanism of SARS-CoV-2 to ACE2, which may aid the design of therapeutics for COVID-19 infection.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Directed Synthesis of Gold Nanoparticle Superstructures Using Self-Assembling Peptoids Containing Metal-Bonding N-Heterocyclic Carbenes

N-Heterocyclic carbene (NHC) ligands, with strong metal-binding affinity, offer a robust platform for constructing organic–inorganic nanohybrids with high stability and tunable properties. However, achieving precise structural control in a simple manner remains challenging. Here, we report a one-pot synthesis of nanohybrids using self-assembling peptoids functionalized with histidine-2-ylidene, which simultaneously enable peptoid assembly and NHC–metal binding. Histidine-2-ylidene-functionalized peptoids were designed to self-assemble and provide NHC binding sites, while AuNPs served as the inorganic component due to NHCs’ strong affinity for gold surfaces. The resulting peptoid-NHC@AuNPs form well-defined vesicles that are characterized by UV–vis spectroscopy, X-ray photoelectron spectroscopy, and electron microscopy. Importantly, the vesicle size and morphology can be tuned via the peptoid sequence or environmental conditions. Further experiments highlight the crucial role of the NHC sites in the formation and stabilization of these nanohybrids. This modular strategy offers a versatile route to fabricating functional NHC-based nanohybrids for potential applications in sensing.

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

Structure of VanS from vancomycin-resistant enterococci: A sensor kinase with weak ATP binding

The VanRS two-component system regulates the resistance phenotype of vancomycin-resistant enterococci. VanS is a sensor histidine kinase that responds to the presence of vancomycin by autophosphorylating and subsequently transferring the phosphoryl group to the response regulator, VanR. The phosphotransfer activates VanR as a transcription factor, which initiates the expression of resistance genes. Structural information about VanS proteins has remained elusive, hindering the molecular-level understanding of their function. Here, we present X-ray crystal structures for the catalytic and ATP-binding (CA) domains of two VanS proteins, derived from vancomycin-resistant enterococci types A and C. Both proteins adopt the canonical Bergerat fold that has been observed for CA domains of other prokaryotic histidine kinases. We attempted to determine structures for the nucleotide-bound forms of both proteins; however, despite repeated efforts, these forms could not be crystallized, prompting us to measure the proteins' binding affinities for ATP. Unexpectedly, both CA domains displayed low affinities for the nucleotide, with K D values in the low millimolar range. Since these K D values are comparable to intracellular ATP concentrations, this weak substrate binding could reflect a way of regulating expression of the resistance phenotype.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Ensemble Simulations and Experimental Free Energy Distributions: Evaluation and Characterization of Isoxazole Amides as SMYD3 Inhibitors

Optimization of binding affinities for ligands to their target protein is a primary objective in rational drug discovery. Herein, we report on a collaborative study that evaluates various compounds designed to bind to the SET and MYND domain-containing protein 3 (SMYD3). SMYD3 is a histone methyltransferase and plays an important role in transcriptional regulation in cell proliferation, cell cycle, and human carcinogenesis. Experimental measurements using the scintillation proximity assay show that the distributions of binding free energies from a large number of independent measurements exhibit non-normal properties. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to predict the binding free energies and to provide a detailed chemical insight into the nature of ligand-protein binding. Our results show that the 1-trajectory ESMACS protocol works well for the set of ligands studied here. Although one unexplained outlier exists, we obtain excellent statistical ranking across the set of compounds from the ESMACS protocol and good agreement between calculations and experiments for the relative binding free energies from the TIES protocol. ESMACS and TIES are again found to be powerful protocols for the accurate comparison of the binding free energies.

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

Extending quantum-mechanical benchmark accuracy to biological ligand-pocket interactions

Predicting the binding affinity of ligands to protein pockets is key in the drug design pipeline. The flexibility of ligand-pocket motifs arises from a range of attractive and repulsive electronic interactions during binding. Accurately accounting for all interactions requires robust quantum-mechanical (QM) benchmarks, which are scarce for ligand-pocket systems. Additionally, disagreement between “gold standard” Coupled Cluster (CC) and Quantum Monte Carlo (QMC) methods casts doubt on many benchmarks for larger non-covalent systems. We introduce the “QUantum Interacting Dimer” (QUID) benchmark framework containing 170 non-covalent (non-)equilibrium systems modeling chemically and structurally diverse ligand-pocket motifs. Symmetry-adapted perturbation theory shows that QUID broadly covers non-covalent binding motifs and energetic contributions. Robust binding energies are obtained using complementary CC and QMC methods, achieving agreement of 0.5 kcal/mol. The benchmark data analysis reveals that several dispersion-inclusive density functional approximations provide accurate energy predictions, though their atomic van der Waals forces differ in magnitude and orientation. Contrarily, semiempirical methods and empirical force fields require improvements in capturing non-covalent interactions (NCIs) for out-of-equilibrium geometries. The wide span of NCIs, highly accurate interaction energies, and analysis of molecular properties take QUID beyond the “gold standard” for QM benchmarks of ligand-protein systems.

Puleva, Mirela [University of Luxembourg, Luxembou↗
DOE OSTI2020

Structural basis for potent neutralization of SARS-CoV-2 and role of antibody affinity maturation

SARS-CoV-2 is a betacoronavirus virus responsible for the COVID-19 pandemic. Here, we determine the X-ray crystal structure of a potent neutralizing monoclonal antibody, CV30, isolated from a patient infected with SARS-CoV-2, in complex with the receptor binding domain. The structure reveals that CV30 binds to an epitope that overlaps with the human ACE2 receptor binding motif providing a structural basis for its neutralization. CV30 also induces shedding of the S1 subunit, indicating an additional mechanism of neutralization. A germline reversion of CV30 results in a substantial reduction in both binding affinity and neutralization potential indicating the minimal somatic mutation is needed for potently neutralizing antibodies against SARS-CoV-2.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Large Scale Study of Ligand–Protein Relative Binding Free Energy Calculations: Actionable Predictions from Statistically Robust Protocols

The accurate and reliable prediction of protein–ligand binding affinities can play a central role in the drug discovery process as well as in personalized medicine. Of considerable importance during lead optimization are the alchemical free energy methods that furnish an estimation of relative binding free energies (RBFE) of similar molecules. Recent advances in these methods have increased their speed, accuracy, and precision. This is evident from the increasing number of retrospective as well as prospective studies employing them. However, such methods still have limited applicability in real-world scenarios due to a number of important yet unresolved issues. Here, we report the findings from a large data set comprising over 500 ligand transformations spanning over 300 ligands binding to a diverse set of 14 different protein targets which furnish statistically robust results on the accuracy, precision, and reproducibility of RBFE calculations. We use ensemble-based methods which are the only way to provide reliable uncertainty quantification given that the underlying molecular dynamics is chaotic. These are implemented using TIES (Thermodynamic Integration with Enhanced Sampling). Results achieve chemical accuracy in all cases. Ensemble simulations also furnish information on the statistical distributions of the free energy calculations which exhibit non-normal behavior. We find that the “enhanced sampling” method known as replica exchange with solute tempering degrades RBFE predictions. We also report definitively on numerous associated alchemical factors including the choice of ligand charge method, flexibility in ligand structure, and the size of the alchemical region including the number of atoms involved in transforming one ligand into another. Our findings provide a key set of recommendations that should be adopted for the reliable application of RBFE methods.

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

High-throughput native mass spectrometry as experimental validation for in silico drug design

In this project, we developed automated workflows for both experimental validation and computational prediction of protein-ligand interactions. The ultimate goal is to establish an integrated pipeline for high throughput design of inhibitors to enzymes relevant to all areas of biological research. Our experimental approach is based on native mass spectrometry (native MS), which measures accurate masses and quantify the relative abundance of protein-ligand complexes to define binding affinity. We set up an in-house built autosampler with highly flexible configurations to minimize the manual steps for high throughput native MS. In parallel, we also performed manual native MS to characterize the binding of substrates and inhibitors of SARS-Cov-2 nonstructural protein nsp10/16 in order to optimize the experimental parameters for future automation. On the computational side, we streamlined the pipeline to achieve minimal manual intervention for predicting enzyme inhibitors via simulation, using the same nsp10/16 system as an example. Using the native MS method we examined 8 top-ranked designed compounds, 2 of which showed weak binding of ~50 µM. The information from native MS experiment provided critical insights and the foundation for a fully integrated workflow for enzyme inhibitor design.

59 BASIC BIOLOGICAL SCIENCES↗