Search NASA⌕ Search

SEARCH · Search NASA

Results for “Drug discovery”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Alchemical Free Energy Estimators and Molecular Dynamics Engines: Accuracy, Precision, and Reproducibility

The binding free energy between a ligand and its target protein is an essential quantity to know at all stages of the drug discovery pipeline. Assessing this value computationally can offer insight into where efforts should be focused in the pursuit of effective therapeutics to treat a myriad of diseases. In this work, we examine the computation of alchemical relative binding free energies with an eye for assessing reproducibility across popular molecular dynamics packages and free energy estimators. The focus of this work is on 54 ligand transformations from a diverse set of protein targets: MCL1, PTP1B, TYK2, CDK2, and thrombin. These targets are studied with three popular molecular dynamics packages: OpenMM, NAMD2, and NAMD3 alpha. Trajectories collected with these packages are used to compare relative binding free energies calculated with thermodynamic integration and free energy perturbation methods. The resulting binding free energies show good agreement between molecular dynamics packages with an average mean unsigned error between them of 0.50 kcal/mol. The correlation between packages is very good, with the lowest Spearman’s, Pearson’s and Kendall’s tau correlation coefficients being 0.92, 0.91, and 0.76, respectively. Agreement between thermodynamic integration and free energy perturbation is shown to be very good when using ensemble averaging.

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

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the lambda-exchange of free energy perturbation (FEP). Further, two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/lambda-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

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

Mechanistic Insights toward Accelerated Selenium-Catalyzed Allylic and Propargylic C–H Amination

Catalytic allylic and propargylic C–H aminations present valuable opportunities for the late-stage modification of pharmacophores in drug discovery. However, modern methodology is limited by reliance on expensive transition metal catalysts or slow reactivity. While selenium-catalyzed methods avoid some of these issues, in their current state, they require high catalyst loadings (15 mol %), superstoichiometric quantities of the amine and oxidant source, and long reaction times. Furthermore, our understanding of the mechanism remains incomplete: e.g., what is the catalytically active species, how is the precatalyst converted to it, and what is the role of the ligand? Here, in this paper, we report an N-heterocyclic carbene selenide (NHC-Se) that allows for substantially reduced loadings of selenium without sacrificing reaction time, improves conversions and yields for challenging substrates, and even exhibits the capacity to catalyze 1,4-allylic diamination of alkenes. To understand the origin of the enhanced activity compared to the state-of-the-art NHC-Se precatalyst, we conducted a mechanistic study that supports ligand-free selenium diimide as the active enophile in these systems, while the NHC-derived byproducts, like the corresponding urea, facilitate turnover. We also isolate and structurally characterize NHC-Se mono- and diimido species and explore their mechanistic role. Thus, this work advances C–H amination methodology, our understanding of its mechanistic underpinnings, and selenium chemistry at large.

allylic amination↗
DOE OSTI2021

Potent Noncovalent Inhibitors of the Main Protease of SARS-CoV-2 from Molecular Sculpting of the Drug Perampanel Guided by Free Energy Perturbation Calculations

Starting from our previous finding of 14 known drugs as inhibitors of the main protease (M pro ) of SARS-CoV-2, the virus responsible for COVID-19, we have redesigned the weak hit perampanel to yield multiple noncovalent, nonpeptidic inhibitors with ca. 20 nM IC 50 values in a kinetic assay. Free-energy perturbation (FEP) calculations for M pro -ligand complexes provided valuable guidance on beneficial modifications that rapidly delivered the potent analogues. The design efforts were confirmed and augmented by determination of high-resolution X-ray crystal structures for five analogues bound to M pro . Results of cell-based antiviral assays further demonstrated the potential of the compounds for treatment of COVID-19. In addition to the possible therapeutic significance, the work clearly demonstrates the power of computational chemistry for drug discovery, especially FEP-guided lead optimization.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

S -Adenosylhomocysteine Analogs Selectively Suppress Pan-Coronavirus Replication by Inhibition of nsp14 Methyltransferase

To address the ongoing threat of SARS-CoV-2 and potential emergence of novel coronaviruses, we employed a comprehensive strategy to identify and synthesize inhibitors of coronavirus methyltransferases with chemical analogs of S-adenosylhomocysteine (SAH). Two analogs, designated 4h and 4p, inhibit both mouse hepatitis virus and SARS-CoV-2 replication. Compound 4p was the most potent with half-maximal inhibition of biochemical activity at 0.2 μM and antiviral activity at ∼20 μM. This compound also has low cytotoxicity and preferentially inhibits nsp14 over nsp16 and human methyltransferases. Furthermore, molecular docking based on a newly determined crystal structure of the apo nsp16−nsp10 complex predicts that 4p occupies both the Sadenosylmethione and Gppp binding pockets of nsp14 and nsp16. Selectivity of 4p for nsp14 is likely due to the enhanced structural stability of the nsp14 binding pocket relative to nsp16. These findings highlight SAH analogs as scaffolds for pan-coronavirus therapeutics and underscore the value of structure-guided design in antiviral drug discovery.

Coronavirus↗
DOE OSTI2020

Chemoinformatic-Guided Engineering of Polyketide Synthases

Polyketide synthase (PKS) engineering is an attractive method to generate new molecules such as commodity, fine and specialty chemicals. A significant challenge is re-engineering a partially reductive PKS module to produce a saturated β-carbon through a reductive loop (RL) exchange. In this work, we sought to establish that chemoinformatics, a field traditionally used in drug discovery, offers a viable strategy for RL exchanges. We first introduced a set of donor RLs of diverse genetic origin and chemical substrates into the first extension module of the lipomycin PKS (LipPKS1). Product titers of these engineered unimodular PKSs correlated with chemical structure similarity between the substrate of the donor RLs and recipient LipPKS1, reaching a titer of 165 mg/L of short-chain fatty acids produced by the host Streptomyces albus J1074. Finally, expanding this method to larger intermediates that require bimodular communication, we introduced RLs of divergent chemosimilarity into LipPKS2 and determined triketide lactone production. Collectively, we observed a statistically significant correlation between atom pair chemosimilarity and production, establishing a new chemoinformatic method that may aid in the engineering of PKSs to produce desired, unnatural products.

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

A Unified Approach to Decarboxylative Halogenation of (Hetero)aryl Carboxylic Acids

Aryl halides are a fundamental motif in synthetic chemistry, playing a critical role in metal-mediated cross-coupling reactions and serving as important scaffolds in drug discovery. Although thermal decarboxylative functionalization of aryl carboxylic acids has been extensively explored, the scope of existing halodecarboxylation methods remains limited, and there currently exists no unified strategy that provides access to any type of aryl halide from an aryl carboxylic acid precursor. Herein, we report a general catalytic method for direct decarboxylative halogenation of (hetero)aryl carboxylic acids via ligand-to-metal charge transfer. This strategy accommodates an exceptionally broad scope of substrates. Here, we leverage an aryl radical intermediate toward divergent functionalization pathways: (1) atom transfer to access bromo- or iodo(hetero)arenes or (2) radical capture by copper and subsequent reductive elimination to generate chloro- or fluoro(hetero)arenes. The proposed ligand-to-metal charge transfer mechanism is supported through an array of spectroscopic studies.

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

Predicting the structural basis of targeted protein degradation by integrating molecular dynamics simulations with structural mass spectrometry

Targeted protein degradation (TPD) is a promising approach in drug discovery for degrading proteins implicated in diseases. A key step in this process is the formation of a ternary complex where a heterobifunctional molecule induces proximity of an E3 ligase to a protein of interest (POI), thus facilitating ubiquitin transfer to the POI. In this work, we characterize 3 steps in the TPD process. (1) We simulate the ternary complex formation of SMARCA2 bromodomain and VHL E3 ligase by combining hydrogen-deuterium exchange mass spectrometry with weighted ensemble molecular dynamics (MD). (2) We characterize the conformational heterogeneity of the ternary complex using Hamiltonian replica exchange simulations and small-angle X-ray scattering. (3) We assess the ubiquitination of the POI in the context of the full Cullin-RING Ligase, confirming experimental ubiquitinomics results. Differences in degradation efficiency can be explained by the proximity of lysine residues on the POI relative to ubiquitin.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

High-throughput platform for yeast morphological profiling predicts the targets of bioactive compounds

Abstract Morphological profiling is an omics-based approach for predicting intracellular targets of chemical compounds in which the dose-dependent morphological changes induced by the compound are systematically compared to the morphological changes in gene-deleted cells. In this study, we developed a reliable high-throughput (HT) platform for yeast morphological profiling using drug-hypersensitive strains to minimize compound use, HT microscopy to speed up data generation and analysis, and a generalized linear model to predict targets with high reliability. We first conducted a proof-of-concept study using six compounds with known targets: bortezomib, hydroxyurea, methyl methanesulfonate, benomyl, tunicamycin, and echinocandin B. Then we applied our platform to predict the mechanism of action of a novel diferulate-derived compound, poacidiene. Morphological profiling of poacidiene implied that it affects the DNA damage response, which genetic analysis confirmed. Furthermore, we found that poacidiene inhibits the growth of phytopathogenic fungi, implying applications as an effective antifungal agent. Thus, our platform is a new whole-cell target prediction tool for drug discovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Artificial intelligence in cancer research, diagnosis and therapy

Artificial intelligence and machine learning techniques are breaking into biomedical research and health care, which importantly includes cancer research and oncology, where the potential applications are vast. These include detection and diagnosis of cancer, subtype classification, optimization of cancer treatment and identification of new therapeutic targets in drug discovery. While big data used to train machine learning models may already exist, leveraging this opportunity to realize the full promise of artificial intelligence in both the cancer research space and the clinical space will first require significant obstacles to be surmounted. In this Viewpoint article, we asked four experts for their opinions on how we can begin to implement artificial intelligence while ensuring standards are maintained so as transform cancer diagnosis and the prognosis and treatment of patients with cancer and to drive biological discovery.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Chemical screening by time-resolved X-ray scattering to discover allosteric probes

Abstract Drug discovery relies on efficient identification of small-molecule leads and their interactions with macromolecular targets. However, understanding how chemotypes impact mechanistically important conformational states often remains secondary among high-throughput discovery methods. Here, we present a conformational discovery pipeline integrating time-resolved, high-throughput small-angle X-ray scattering (TR-HT-SAXS) and classic fragment screening applied to allosteric states of the mitochondrial import oxidoreductase apoptosis-inducing factor (AIF). By monitoring oxidized and X-ray-reduced AIF states, TR-HT-SAXS leverages structure and kinetics to generate a multidimensional screening dataset that identifies fragment chemotypes allosterically stimulating AIF dimerization. Fragment-induced dimerization rates, quantified with time-resolved SAXS similarity analysis ( k VR ), capture structure–activity relationships (SAR) across the top-ranked 4-aminoquinoline chemotype. Crystallized AIF–aminoquinoline complexes validate TR-SAXS-guided SAR, supporting this conformational chemotype for optimization. AIF–aminoquinoline structures and mutational analysis reveal active site F482 as an underappreciated allosteric stabilizer of AIF dimerization. This conformational discovery pipeline illustrates TR-HT-SAXS as an effective technology for targeting chemical leads to important macromolecular states.

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

Cryo2StructData: A Large Labeled Cryo-EM Density Map Dataset for AI-based Modeling of Protein Structures

The advent of single-particle cryo-electron microscopy (cryo-EM) has brought forth a new era of structural biology, enabling the routine determination of large biological molecules and their complexes at atomic resolution. The high-resolution structures of biological macromolecules and their complexes significantly expedite biomedical research and drug discovery. However, automatically and accurately building atomic models from high-resolution cryo-EM density maps is still time-consuming and challenging when template-based models are unavailable. Artificial intelligence (AI) methods such as deep learning trained on limited amount of labeled cryo-EM density maps generate inaccurate atomic models. To address this issue, we created a dataset called Cryo2StructData consisting of 7,600 preprocessed cryo-EM density maps whose voxels are labelled according to their corresponding known atomic structures for training and testing AI methods to build atomic models from cryo-EM density maps. Cryo2StructData is larger than existing, publicly available datasets for training AI methods to build atomic protein structures from cryo-EM density maps. We trained and tested deep learning models on Cryo2StructData to validate its quality showing that it is ready for being used to train and test AI methods for building atomic models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Correlation between hemolytic activity, cytotoxicity and systemic in vivo toxicity of synthetic antimicrobial peptides

The use of non-standard toxicity models is a hurdle in the early development of antimicrobial peptides towards clinical applications. Herein we report an extensive in vitro and in vivo toxicity study of a library of 24 peptide-based antimicrobials with narrow spectrum activity towards veterinary pathogens. The haemolytic activity of the compounds was evaluated against four different species and the relative sensitivity against the compounds was highest for canine erythrocytes, intermediate for rat and human cells and lowest for bovine cells. Selected peptides were additionally evaluated against HeLa, HaCaT and HepG2 cells which showed increased stability towards the peptides. Therapeutic indexes of 50–500 suggest significant cellular selectivity in comparison to bacterial cells. Three peptides were administered to rats in intravenous acute dose toxicity studies up to 2– 8 × MIC. None of the injected compounds induced any systemic toxic effects in vivo at the concentrations employed illustrating that the correlation between the different assays is not obvious. This work sheds light on the in vitro and in vivo toxicity of this class of promising compounds and provides insights into the relationship between the different toxicity models often employed in different manners to evaluate the toxicity of novel bioactive compounds in general.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Adaptive hyperparameter updating for training restricted Boltzmann machines on quantum annealers

Restricted Boltzmann Machines (RBMs) have been proposed for developing neural networks for a variety of unsupervised machine learning applications such as image recognition, drug discovery, and materials design. The Boltzmann probability distribution is used as a model to identify network parameters by optimizing the likelihood of predicting an output given hidden states trained on available data. Training such networks often requires sampling over a large probability space that must be approximated during gradient based optimization. Quantum annealing has been proposed as a means to search this space more efficiently which has been experimentally investigated on D-Wave hardware. D-Wave implementation requires selection of an effective inverse temperature or hyperparameter (β) within the Boltzmann distribution which can strongly influence optimization. Here, we show how this parameter can be estimated as a hyperparameter applied to D-Wave hardware during neural network training by maximizing the likelihood or minimizing the Shannon entropy. We find both methods improve training RBMs based upon D-Wave hardware experimental validation on an image recognition problem. Neural network image reconstruction errors are evaluated using Bayesian uncertainty analysis which illustrate more than an order magnitude lower image reconstruction error using the maximum likelihood over manually optimizing the hyperparameter. The maximum likelihood method is also shown to out-perform minimizing the Shannon entropy for image reconstruction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

AI-accelerated protein-ligand docking for SARS-CoV-2 is 100-fold faster with no significant change in detection

Protein-ligand docking is a computational method for identifying drug leads. The method is capable of narrowing a vast library of compounds down to a tractable size for downstream simulation or experimental testing and is widely used in drug discovery. While there has been progress in accelerating scoring of compounds with artificial intelligence, few works have bridged these successes back to the virtual screening community in terms of utility and forward-looking development. We demonstrate the power of high-speed ML models by scoring 1 billion molecules in under a day (50 k predictions per GPU seconds). We showcase a workflow for docking utilizing surrogate AI-based models as a pre-filter to a standard docking workflow. Our workflow is ten times faster at screening a library of compounds than the standard technique, with an error rate less than 0.01% of detecting the underlying best scoring 0.1% of compounds. Our analysis of the speedup explains that another order of magnitude speedup must come from model accuracy rather than computing speed. In order to drive another order of magnitude of acceleration, we share a benchmark dataset consisting of 200 million 3D complex structures and 2D structure scores across a consistent set of 13 million “in-stock” molecules over 15 receptors, or binding sites, across the SARS-CoV-2 proteome. We believe this is strong evidence for the community to begin focusing on improving the accuracy of surrogate models to improve the ability to screen massive compound libraries 100 × or even 1000 × faster than current techniques and reduce missing top hits. The technique outlined aims to be a fast drop-in replacement for docking for screening billion-scale molecular libraries.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Global characterization of biosynthetic gene clusters in non-model eukaryotes using domain architectures

Abstract The majority of pharmaceuticals are derived from natural products, bioactive compounds naturally synthesized by organisms to provide evolutionary advantages. Although the rich evolutionary history of eukaryotic algal species implicates a high potential for natural product-based drug discovery, it remains largely untouched. This study investigates 2762 putative biosynthetic gene clusters (BGCs) from 212 eukaryotic algal genomes. To analyze a vast set of structurally diverse BGCs, we employed comparative analysis based on the vectorization of biosynthetic domains, referred to as biosynthetic domain architecture (BDA). By characterizing core biosynthetic machineries through BDA, we identified key BDAs of modular BGCs in diverse eukaryotes and introduced 16 candidate modular BGCs with similar BDAs to previously validated BGCs. This study provides a global characterization of eukaryotic algal BGCs, offering an alternative to laborious manual curation for BGC prioritization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Huntingtin structure is orchestrated by HAP40 and shows a polyglutamine expansion-specific interaction with exon 1

Huntington’s disease results from expansion of a glutamine-coding CAG tract in the huntingtin (HTT) gene, producing an aberrantly functioning form of HTT. Both wildtype and disease-state HTT form a hetero-dimer with HAP40 of unknown functional relevance. Here, we demonstrate in vivo and in cell models that HTT and HAP40 cellular abundance are coupled. Integrating data from a 2.6 Å cryo-electron microscopy structure, cross-linking mass spectrometry, small-angle X-ray scattering, and modeling, we provide a near-atomic-level view of HTT, its molecular interaction surfaces and compacted domain architecture, orchestrated by HAP40. Native mass spectrometry reveals a remarkably stable hetero-dimer, potentially explaining the cellular inter-dependence of HTT and HAP40. The exon 1 region of HTT is dynamic but shows greater conformational variety in the polyglutamine expanded mutant than wildtype exon 1. Our data provide a foundation for future functional and drug discovery studies targeting Huntington’s disease and illuminate the structural consequences of HTT polyglutamine expansion.

59 BASIC BIOLOGICAL SCIENCES↗