Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 OSTI2020

CHARMM-GUI Free Energy Calculator for Absolute and Relative Ligand Solvation and Binding Free Energy Simulations

Alchemical free energy simulations have long been utilized to predict free energy changes for binding affinity and solubility of small molecules. However, while the theoretical foundation of these methods is well established, seamlessly handling many of the practical aspects regarding the preparation of the different thermodynamic end states of complex molecular systems and the numerous processing scripts often remains a burden for successful applications. Here, we present CHARMM-GUI Free Energy Calculator (http://www.charmm-gui.org/input/fec) that provides various alchemical free energy perturbation molecular dynamics (FEP/MD) systems with input and post-processing scripts for NAMD and GENESIS. Four submodules are available: Absolute Ligand Binder (for absolute ligand binding FEP/MD), Relative Ligand Binder (for relative ligand binding FEP/MD), Absolute Ligand Solvator (for absolute ligand solvation FEP/MD), and Relative Ligand Solvator (for relative ligand solvation FEP/MD). Each module is designed to build multiple systems of a set of selected ligands at once for high-throughput FEP/MD simulations. The capability of Free Energy Calculator is illustrated by absolute and relative solvation FEP/MD of a set of ligands and absolute and relative binding FEP/MD of a set of ligands for T4-lysozyme in solution and the adenosine A 2A receptor in a membrane. The calculated free energy values are overall consistent with the experimental and published free energy results (within ~1 kcal/mol). We hope that Free Energy Calculator is useful to carry out high-throughput FEP/MD simulations in the field of biomolecular sciences and drug discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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

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↗