Search NASA⌕ Search

SEARCH · Search NASA

Results for “drug binding”

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 91 records · Page 5

Characterization of Lactobacilli Phage Endolysins and Their Functional Domains–Potential Live Biotherapeutic Testing Reagents

Phage endolysin-specific binding characteristics and killing activity support their potential use in biotechnological applications, including potency and purity testing of live biotherapeutic products (LBPs). LBPs contain live organisms, such as lactic acid bacteria (LAB), and are intended for use as drugs. Our approach uses the endolysin cell wall binding domains (CBD) for LBP potency assays and the endolysin killing activity for purity assays. CBDs of the following five lactobacilli phage lysins were characterized: CL1, Jlb1, Lj965, LL-H, and ΦJB. They exhibited different bindings to 27 LAB strains and were found to bind peptidoglycan or surface polymers. Flow cytometry based on CBD binding was used to enumerate viable counts of two strains in the mixture. CL1-lys, jlb1-lys, and ΦJB-lys and their enzymatic domains (EADs) exhibited cell wall digestive activity and lytic activity against LAB. Jlb1-EAD and ΦJB-EAD were more sensitive than their respective hololysins to buffer pH and NaCl changes. The ΦJB-EAD exhibited stronger lytic activity than ΦJB-lys, possibly due to ΦJB-CBD-mediated sequestration of ΦJB-lys by cell debris. CBD multiplex assays indicate that these proteins may be useful LBP potency reagents, and the lytic activity suggests that CL1-lys, jlb1-lys, and ΦJB-lys and their EADs are good candidates for LBP purity reagent development.

60 APPLIED LIFE SCIENCES↗

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↗

Revealing the mechanism of action of a first-in-class covalent inhibitor of $\mathrm{KRASG12C}$ $\mathrm{(ON)}$ and other functional properties of oncogenic $\text{KRAS}$ by 31 P $\mathrm{NMR}$

Individual oncogenic KRAS mutants confer distinct differences in biochemical properties and signaling for reasons that are not well understood. KRAS activity is closely coupled to protein dynamics and is regulated through two interconverting conformations: state 1 (inactive, effector binding deficient) and state 2 (active, effector binding enabled). Here, we use 31 P NMR to delineate the differences in state 1 and state 2 populations present in WT and common KRAS oncogenic mutants (G12C, G12D, G12V, G13D, and Q61L) bound to its natural substrate GTP or a commonly used nonhydrolyzable analog GppNHp (guanosine-5'-[(β,γ)-imido] triphosphate). Our results show that GppNHp-bound proteins exhibit significant state 1 population, whereas GTP-bound KRAS is primarily (90% or more) in state 2 conformation. This observation suggests that the predominance of state 1 shown here and in other studies is related to GppNHp and is most likely nonexistent in cells. We characterize the impact of this differential conformational equilibrium of oncogenic KRAS on RAF1 kinase effector RAS-binding domain and intrinsic hydrolysis. Through a KRAS G12C drug discovery, we have identified a novel small-molecule inhibitor, BBO-8956, which is effective against both GDP- and GTP-bound KRAS G12C. We show that binding of this inhibitor significantly perturbs state 1–state 2 equilibrium and induces an inactive state 1 conformation in GTP-bound KRAS G12C. In the presence of BBO-8956, RAF1–RAS-binding domain is unable to induce a signaling competent state 2 conformation within the ternary complex, demonstrating the mechanism of action for this novel and active-conformation inhibitor.

31P NMR↗

Defining the Antitumor Mechanism of Action of a Clinical-stage Compound as a Selective Degrader of the Nuclear Pore Complex

Cancer cells are acutely dependent on nuclear transport due to elevated transcriptional activity, suggesting an unrealized opportunity for selective therapeutic inhibition of the nuclear pore complex (NPC). Through large-scale phenotypic profiling of cancer cell lines, genome-scale functional genomic modifier screens, and mass spectrometry–based proteomics, we discovered that the clinical drug PRLX-93936 is a molecular glue that binds and reprograms the TRIM21 ubiquitin ligase to degrade the NPC. Upon compound-induced TRIM21 recruitment, the nuclear pore is ubiquitylated and degraded, resulting in the loss of short-lived cytoplasmic mRNA transcripts and the induction of cancer cell apoptosis. Direct compound binding to TRIM21 was confirmed via surface plasmon resonance and X-ray crystallography, whereas compound-induced TRIM21–nucleoporin complex formation was demonstrated through multiple orthogonal approaches in cells and in vitro. Phenotype-guided optimization yielded compounds with 10-fold greater potency and drug-like properties, along with robust pharmacokinetics and efficacy against pancreatic cancer xenografts and patient-derived organoids.

Yuan, Linjie [Stanford School of Medicine, CA (Uni↗

The performance of ensemble-based free energy protocols in computing binding affinities to ROS1 kinase

Optimization of binding affinities for compounds to their target protein is a primary objective in drug discovery. Herein we report on a collaborative study that evaluates a set of compounds binding to ROS1 kinase. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to rank the binding free energies. The predicted binding free energies from ESMACS simulations show good correlations with experimental data for subsets of the compounds. Consistent binding free energy differences are generated for TIES and ESMACS. Although an unexplained overestimation exists, we obtain excellent statistical rankings across the set of compounds from the TIES protocol, with a Pearson correlation coefficient of 0.90 between calculated and experimental activities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of SARS-CoV-2 Main Protease Mutations at Positions L50, E166, and L167 Rendering Resistance to Covalent and Noncovalent Inhibitors

SARS-CoV-2 propagation under nirmatrelvir and ensitrelvir pressure selects for main protease (MPro) drug-resistant mutations E166V (DRM2), L50F/E166V (DRM3), E166A/L167F (DRM4), and L50F/E166A/L167F (DRM5). DRM2-DRM5 undergoes N-terminal autoprocessing to produce mature MPro with dimer dissociation constants (K dimer ) 2–3 times larger than that of the wildtype. Co-selection of L50F restores catalytic activity of DRM2 and DRM4 from ~10 to 30%, relative to that of the wild-type enzyme, without altering K dimer . Binding affinities and thermodynamic profiles that parallel the drug selection pressure, exhibiting significant decreases in affinity through entropy/enthalpy compensation, were compared with GC373. Reorganization of the active sites due to mutations observed in the inhibitor-free DRM3 and DRM4 structures as compared to MPro WT may account for the reduced binding affinities, although DRM2 and DRM3 complexes with ensitrelvir are almost identical to MPro WT -ensitrelvir. In conclusion, chemical reactivity changes of the mutant active sites due to differences in electrostatic and protein dynamics effects likely contribute to losses in binding affinities.

60 APPLIED LIFE SCIENCES↗

Challenges for Implementing FAIR Digital Objects with High Performance Workflows

New types of workflows are being used in science that couple traditional distributed and high-performance computing (HPC) with data-intensive approaches, and orchestrate ensembles of numerical simulations and artificial intelligence (AI) models. Such workflows may use AI models to supplement computation where numerical simulations may be too computationally expensive, to automate trivial yet time consuming operations, to perform preliminary selections among intractable numbers of combinations in domains as diverse as protein binding, fine-grid climate simulations, and drug discovery.

97 MATHEMATICS AND COMPUTING↗

Transferrin receptor targeting by de novo sheet extension

The de novo design of polar protein–protein interactions is challenging because of the thermodynamic cost of stripping water away from the polar groups. Here, we describe a general approach for designing proteins which complement exposed polar backbone groups at the edge of beta sheets with geometrically matched beta strands. We used this approach to computationally design small proteins that bind to an exposed beta sheet on the human transferrin receptor (hTfR), which shuttles interacting proteins across the blood–brain barrier (BBB), opening up avenues for drug delivery into the brain. We describe a design which binds hTfR with a 20 nM K d , is hyperstable, and crosses an in vitro microfluidic organ-on-a-chip model of the human BBB. Our design approach provides a general strategy for creating binders to protein targets with exposed surface beta edge strands.

59 BASIC BIOLOGICAL SCIENCES↗

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↗

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↗

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↗

A non-traditional crystal-based compound screening method targeting the ATP binding site of Plasmodium falciparum GRP78 for identification of novel nucleoside analogues

Drug resistance to front-line malarial treatments represents an ongoing threat to control malaria, a vector borne infectious disease. The malarial parasite, Plasmodium falciparum has developed genetic variants, conferring resistance to the current standard therapeutic artemisinin and its derivatives commonly referred to as artemisinin-combination therapies (ACTs). Emergence of multi-drug resistance parasite genotypes is a warning of potential treatment failure, reaffirming the urgent and critical need to find and validate alternate drug targets to prevent the spread of disease. An attractive and novel drug target includes glucose-regulated protein 78 kDa (GRP78, or BiP), an essential molecular chaperone protein involved in the unfolded protein response that is upregulated in ACT treated P. falciparum parasites. We have shown that both sequence and structure are closely related to human GRP78 (hGRP78), a chaperone belonging to the HSP70 class of ATPase proteins, which is often upregulated in cellular stress responses and cancer. By screening a library of nucleoside analogues, we identified eight ‘hit’ compounds binding at the active site of the ATP binding domain of P. falciparum GRP78 using a high-throughput ligand soaking screen using x-ray crystallography. These compounds were further evaluated using protein thermal shift assays to assess target binding activity. The nucleoside analogues identified from our screen provide a starting point for the development of more potent and selective antimalarial inhibitors. In addition, we have established a well-defined, high-throughput crystal-based screening approach that can be applied to many crystallizable P. falciparum proteins for generating anti-Plasmodium specific compounds.

59 BASIC BIOLOGICAL SCIENCES↗

Structure- and Interaction-Based Design of Anti-SARS-CoV-2 Aptamers

Aptamer selection against novel infections is a complicated and time-consuming approach. Synergy can be achieved by using computational methods together with experimental procedures. In this study, we aim to develop a reliable methodology for a rational aptamer in silico et vitro design. The new approach combines multiple steps: (1) Molecular design, based on screening in a DNA aptamer library and directed mutagenesis to fit the protein tertiary structure; (2) 3D molecular modeling of the target; (3) Molecular docking of an aptamer with the protein; (4) Molecular dynamics (MD) simulations of the complexes; (5) Quantum-mechanical (QM) evaluation of the interactions between aptamer and target with further analysis; (6) Experimental verification at each cycle for structure and binding affinity by using small-angle X-ray scattering, cytometry, and fluorescence polarization. By using a new iterative design procedure, structure- and interaction-based drug design (SIBDD), a highly specific aptamer to the receptor-binding domain of the SARS-CoV-2 spike protein, was developed and validated. The SIBDD approach enhances speed of the high-affinity aptamers development from scratch, using a target protein structure. The method could be used to improve existing aptamers for stronger binding. This approach brings to an advanced level the development of novel affinity probes, functional nucleic acids. It offers a blueprint for the straightforward design of targeting molecules for new pathogen agents and emerging variants.

60 APPLIED LIFE SCIENCES↗

Subcutaneous administration of insulin-like growth factor (IGF)-II/IGF binding protein-2 complex stimulates bone formation and prevents loss of bone mineral density in a rat model of disuse osteoporosis

Elevated serum levels of insulin-like growth factor binding protein-2 (IGFBP-2) and a precursor form of IGF-II are associated with marked increases in bone formation and skeletal mass in patients with hepatitis C-associated osteosclerosis. In vitro studies indicate that IGF-II in complex with IGFBP-2 has high affinity for bone matrix and is able to stimulate osteoblast proliferation. The purpose of this study was to determine the ability of the IGF-II/IGFBP-2 complex to increase bone mass in vivo. Osteopenia of the femur was induced by unilateral sciatic neurectomy in rats. At the time of surgery, 14-day osmotic minipumps containing vehicle or 2 microg IGF-II+9 microg IGFBP-2/100g body weight/day were implanted subcutaneously in the neck. Bone mineral density (BMD) measurements were taken the day of surgery and 14 days later using a PIXImus small animal densitometer. Neurectomy of the right hindlimb resulted in a 9% decrease in right femur BMD (P<0.05 vs. baseline). This loss in BMD was completely prevented by treatment with IGF-II/IGFBP-2. On the control limb, there was no loss of BMD over the 14 days and IGF-II/IGFBP-2 treatment resulted in a 9% increase in left femur BMD (P<0.05). Bone histomorphometry indicated increases in endocortical and cancellous bone formation rates and in trabecular thickness. These results demonstrate that short-term administration of the IGF-II/IGFBP-2 complex can prevent loss of BMD associated with disuse osteoporosis and stimulate bone formation in adult rats. Furthermore, they provide proof of concept for a novel anabolic approach to increasing bone mass in humans with osteoporosis.

NASA Discipline Regulatory Physiology↗

Pathways and Mechanism of Caffeine Binding to Human Adenosine A 2A Receptor

Caffeine (CFF) is a common antagonist to the four subtypes of adenosine G-protein-coupled receptors (GPCRs), which are critical drug targets for treating heart failure, cancer, and neurological diseases. However, the pathways and mechanism of CFF binding to the target receptors remain unclear. In this study, we have performed all-atom-enhanced sampling simulations using a robust Gaussian-accelerated molecular dynamics (GaMD) method to elucidate the binding mechanism of CFF to human adenosine A 2A receptor (A 2A AR). Multiple 500–1,000 ns GaMD simulations captured both binding and dissociation of CFF in the A 2A AR. The GaMD-predicted binding poses of CFF were highly consistent with the x-ray crystal conformations with a characteristic hydrogen bond formed between CFF and residue N6.55 in the receptor. In addition, a low-energy intermediate binding conformation was revealed for CFF at the receptor extracellular mouth between ECL2 and TM1. While the ligand-binding pathways of the A 2A AR were found similar to those of other class A GPCRs identified from previous studies, the ECL2 with high sequence divergence serves as an attractive target site for designing allosteric modulators as selective drugs of the A 2A AR.

59 BASIC BIOLOGICAL SCIENCES↗

Adaptive language model training for molecular design

Abstract The vast size of chemical space necessitates computational approaches to automate and accelerate the design of molecular sequences to guide experimental efforts for drug discovery. Genetic algorithms provide a useful framework to incrementally generate molecules by applying mutations to known chemical structures. Recently, masked language models have been applied to automate the mutation process by leveraging large compound libraries to learn commonly occurring chemical sequences (i.e., using tokenization) and predict rearrangements (i.e., using mask prediction). Here, we consider how language models can be adapted to improve molecule generation for different optimization tasks. We use two different generation strategies for comparison, fixed and adaptive. The fixed strategy uses a pre-trained model to generate mutations; the adaptive strategy trains the language model on each new generation of molecules selected for target properties during optimization. Our results show that the adaptive strategy allows the language model to more closely fit the distribution of molecules in the population. Therefore, for enhanced fitness optimization, we suggest the use of the fixed strategy during an initial phase followed by the use of the adaptive strategy. We demonstrate the impact of adaptive training by searching for molecules that optimize both heuristic metrics, drug-likeness and synthesizability, as well as predicted protein binding affinity from a surrogate model. Our results show that the adaptive strategy provides a significant improvement in fitness optimization compared to the fixed pre-trained model, empowering the application of language models to molecular design tasks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗