Search NASA⌕ Search

SEARCH · Search NASA

Results for “Drug Design”

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 199 records · Page 11

Chapter Twelve - Proton transfer and drug binding details revealed in neutron diffraction studies of wild-type and drug resistant HIV-1 protease

HIV-1 protease is an essential therapeutic target for the design and development of antiviral inhibitors to treat AIDS. We used room temperature neutron crystallography to accurately determine hydrogen atom positions in several protease complexes with clinical drugs, amprenavir and darunavir. Hydrogen bonding interactions were carefully mapped to provide an unprecedented picture of drug binding to the protease target. We demonstrate that hydrogen atom positions within the enzyme catalytic site can be altered by introducing drug resistant mutations and by protonating surface residues that trigger proton transfer reactions between the catalytic Asp residues and the hydroxyl group of darunavir. When protein perdeuteration is not feasible, we validate the use of initial H/D exchange with unfolded protein and partial deuteration in pure D2O with hydrogenous glycerol to maximize deuterium incorporation into the protein, with no detrimental effects on the growth of quality crystals suitable for neutron diffraction experiments.

60 APPLIED LIFE SCIENCES↗

Studying Host-Pathogen Interactions In 3-D: Organotypic Models For Infectious Disease And Drug Development

Representative, reproducible and high-throughput models of human cells and tissues are critical for a meaningful evaluation of host-pathogen interactions and are an essential component of the research developmental pipeline. The most informative infection models - animals, organ explants and human trials - are not suited for extensive evaluation of pathogenesis mechanisms and screening of candidate drugs. At the other extreme, more cost effective and accessible infection models such as conventional cell culture and static co-culture may not capture physiological and three-dimensional aspects of tissue biology that are important in assessing pathogenesis, and effectiveness and cytotoxicity of therapeutics. Our lab has used innovative bioengineering technology to establish biologically meaningful 3-D models of human tissues that recapitulate many aspects of the differentiated structure and function of the parental tissue in vivo, and we have applied these models to study infectious disease. We have established a variety of different 3-D models that are currently being used in infection studies - including small intestine, colon, lung, placenta, bladder, periodontal ligament, and neuronal models. Published work from our lab has shown that our 3-D models respond to infection with bacterial and viral pathogens in ways that reflect the infection process in vivo. By virtue of their physiological relevance, 3-D cell cultures may also hold significant potential as models to provide insight into the neuropathogenesis of HIV infection. Furthermore, the experimental flexibility, reproducibility, cost-efficiency, and high throughput platform afforded by these 3-D models may have important implications for the design and development of drugs with which to effectively treat neurological complications of HIV infection.

Nickerson, Cheryl A.↗

Drugsniffer: An Open Source Workflow for Virtually Screening Billions of Molecules for Binding Affinity to Protein Targets

The SARS-CoV2 pandemic has highlighted the importance of efficient and effective methods for identification of therapeutic drugs, and in particular has laid bare the need for methods that allow exploration of the full diversity of synthesizable small molecules. While classical high-throughput screening methods may consider up to millions of molecules, virtual screening methods hold the promise of enabling appraisal of billions of candidate molecules, thus expanding the search space while concurrently reducing costs and speeding discovery. Here, we describe a new screening pipeline, called drugsniffer, that is capable of rapidly exploring drug candidates from a library of billions of molecules, and is designed to support distributed computation on cluster and cloud resources. As an example of performance, our pipeline required ~40,000 total compute hours to screen for potential drugs targeting three SARS-CoV2 proteins among a library of ~3.7 billion candidate molecules.

59 BASIC BIOLOGICAL SCIENCES↗

Microfluidic synthesis of multilayered lipid–polymer hybrid nanoparticles for the formulation of low solubility drugs

Hybrid phospholipid/block copolymer membranes where polymers and lipids are molecularly mixed or phase-separated into polymer-rich and lipid-rich domains are promising drug delivery materials. Harnessing the chemical diversity of polymers and the biocompatability of lipids is a compelling approach to design the next generation of drug carriers. Here, we report on the development of a microfluidics-based strategy analogous to produce lipid nanoparticles (LNPs) for the nanomanufacturing of multilayered hybrid nanoparticles (HNPs). In conclusion, using X-ray scattering, Cryo-electron, and polarized microscopy we show that phosphatidylcholine (PC) and PBD-b-PEO (poly(butadiene-block-ethylene oxide)) hybrid membranes can be nanomanufactured by microfluidics into HNPs with dense and multilayered cores which are ideal carriers of low-solubility drugs of the Biopharmaceutical Classification System (BCS) II and IV such as antimalarial DSM265 and Paclitaxel, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Heterotropic roles of divalent cations in the establishment of allostery and affinity maturation of integrin αXβ2

Allosteric activation and silencing of leukocyte β2-integrins transpire through cation-dependent structural changes, which mediate integrin biosynthesis and recycling, and are essential to designing leukocyte-specific drugs. Stepwise addition of Mg 2+ reveals two mutually coupled events for the αXβ2 ligand-binding domain—the αX I-domain—corresponding to allostery establishment and affinity maturation. Electrostatic alterations in the Mg 2+ -binding site establish long-range couplings, leading to both pH- and Mg 2+ -occupancy-dependent biphasic stability change in the αX I-domain fold. The ligand-binding sensorgrams show composite affinity events for the αX I-domain accounting for the multiplicity of the αX I-domain conformational states existing in the solution. On cell surfaces, increasing Mg 2+ concentration enhanced adhesiveness of αXβ2. This work highlights how intrinsically flexible pH- and cation-sensitive architecture endows a unique dynamic continuum to the αI-domain structure on the intact integrin, thereby revealing the importance of allostery establishment and affinity maturation in both extracellular and intracellular integrin events.

36 MATERIALS SCIENCE↗

Predicting receptor-ligand pairing preferences in plant-microbe interfaces via molecular dynamics and machine learning

Microbiome assembly, structure, and dynamics significantly influence plant health. Secreted microbial signaling molecules initiate and mediate symbiosis by binding to structurally compatible plant receptors. For example, lipo-chitooligosaccharides (LCOs), produced by nitrogen-fixing rhizobial bacteria and various fungi, are recognized by plant lysin motif receptor-like kinases (LysM-RLKs), which activate the common symbiotic pathway. Accurately predicting these molecular interactions could reveal complementary signatures underlying the initial stages of endosymbiosis. Despite the breakthrough in protein-ligand structure prediction with deep learning-based tools, such as AlphaFold3, the large size and highly flexible nature of signaling compounds like LCOs present major challenges for detailed structural characterization and binding-affinity prediction. Typical structure-/physics-based methods of ligand virtual screening are designed for small, drug-like molecules, often rely on high-resolution, experimentally determined structures of the protein receptors, and rarely achieve sufficient sampling to obtain converged thermodynamic quantities with large ligands. In this study, we developed a hybrid molecular dynamics/machine learning (MD/ML) approach capable of predicting binding affinity rankings with high accuracy in systems involving large, flexible ligands, despite limited experimental structural information. Using coarse initial structural models, the predictions using the MD/ML workflow achieved strong alignment with experimental trends, particularly in the top-affinity tier for four legume LysM-RLKs (LYR3) binding to LCOs and a chitooligosaccharide. Furthermore, the MD-based conformation selection protocol provided critical structural insights into substrate specificity and binding mechanisms. This study demonstrates a powerful method to screen for challenging cognate ligand-receptors and advance our understanding of the molecular basis of microbial colonization in plants.

Lipo-chitooligosaccharides↗

Machine Learning Models for Predicting Molecular UV–Vis Spectra with Quantum Mechanical Properties

Accurate understanding of Ultraviolet–visible (UV–Vis) spectra is critical for highthroughput design of compounds for drug discovery. Experimentally determining UV–Vis spectra can become expensive when dealing with a large quantity of novel molecules. This provides us an opportunity to drive computational advances in molecular property predictions using quantum mechanics and machine learning. In this work, we use both Quantum Mechanically (QM) predicted and measured UV–Vis spectra as input to modify four different machine learning architectures: UVvis-SchNet, UVvis- DTNN, UVvis-Transformer, and UVvis-MPNN. Here we find that the UVvis-MPNN model outperforms the other models when using optimized 3D coordinates and QM predicted spectra as input features. This model has the highest performance for predicting UVVisible spectra with a training RMSE of 0.06 and validation RMSE of 0.08.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural, Electronic, and Electrostatic Determinants for Inhibitor Binding to Subsites S1 and S2 in SARS-CoV-2 Main Protease

Creating small-molecule antivirals specific for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins is crucial to battle coronavirus disease 2019 (COVID-19). SARS-CoV-2 main protease (Mpro) is an established drug target for the design of protease inhibitors. We performed a structure–activity relationship (SAR) study of noncovalent compounds that bind in the enzyme’s substrate-binding subsites S1 and S2, revealing structural, electronic, and electrostatic determinants of these sites. The study was guided by the X-ray/neutron structure of Mpro complexed with Mcule-5948770040 (compound 1), in which protonation states were directly visualized. Virtual reality-assisted structure analysis and small-molecule building were employed to generate analogues of 1. In vitro enzyme inhibition assays and room-temperature X-ray structures demonstrated the effect of chemical modifications on Mpro inhibition, showing that (1) maintaining correct geometry of an inhibitor’s P1 group is essential to preserve the hydrogen bond with the protonated His163; (2) a positively charged linker is preferred; and (3) subsite S2 prefers nonbulky modestly electronegative groups.

60 APPLIED LIFE SCIENCES↗

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↗

Assessment of enzyme active site positioning and tests of catalytic mechanisms through X-ray–derived conformational ensembles

How enzymes achieve their enormous rate enhancements remains a central question in biology, and our understanding to date has impacted drug development, influenced enzyme design, and deepened our appreciation of evolutionary processes. While enzymes position catalytic and reactant groups in active sites, physics requires that atoms undergo constant motion. Numerous proposals have invoked positioning or motions as central for enzyme function, but a scarcity of experimental data has limited our understanding of positioning and motion, their relative importance, and their changes through the enzyme’s reaction cycle. To examine positioning and motions and test catalytic proposals, we collected “room temperature” X-ray crystallography data for Pseudomonas putida ketosteroid isomerase (KSI), and we obtained conformational ensembles for this and a homologous KSI from multiple PDB crystal structures. Here, ensemble analyses indicated limited change through KSI’s reaction cycle. Active site positioning was on the 1- to 1.5-Å scale, and was not exceptional compared to noncatalytic groups. The KSI ensembles provided evidence against catalytic proposals invoking oxyanion hole geometric discrimination between the ground state and transition state or highly precise general base positioning. Instead, increasing or decreasing positioning of KSI’s general base reduced catalysis, suggesting optimized Ångstrom-scale conformational heterogeneity that allows KSI to efficiently catalyze multiple reaction steps. Ensemble analyses of surrounding groups for WT and mutant KSIs provided insights into the forces and interactions that allow and limit active-site motions. Most generally, this ensemble perspective extends traditional structure–function relationships, providing the basis for a new era of “ensemble–function” interrogation of enzymes.

59 BASIC BIOLOGICAL SCIENCES↗

Biophysical and X-ray structural studies of the (GGGTT)3GGG G-quadruplex in complex with N-methyl mesoporphyrin IX

The G-quadruplex (GQ) is a well-studied non-canonical DNA structure formed by G-rich sequences found at telomeres and gene promoters. Biological studies suggest that GQs may play roles in regulating gene expression, DNA replication, and DNA repair. Small molecule ligands were shown to alter GQ structure and stability and thereby serve as novel therapies, particularly against cancer. In this work, we investigate the interaction of a G-rich sequence, 5’-GGGTTGGGTTGGGTTGGG-3’ (T1), with a water-soluble porphyrin, N- methyl mesoporphyrin IX (NMM) via biophysical and X-ray crystallographic studies. UV-vis and fluorescence titrations, as well as a Job plot, revealed a 1:1 binding stoichiometry with an impressively tight binding constant of 30–50 μM -1 and ΔG 298 of -10.3 kcal/mol. Eight extended variants of T1 (named T2 –T9) were fully characterized and T7 was identified as a suitable candidate for crystallographic studies. We solved the crystal structures of the T1- and T7-NMM complexes at 2.39 and 2.34 Å resolution, respectively. Both complexes form a 5’-5’ dimer of parallel GQs capped by NMM at the 3’ G-quartet, supporting the 1:1 binding stoichiometry. Our work provides invaluable details about GQ-ligand binding interactions and informs the design of novel anticancer drugs that selectively recognize specific GQs and modulate their stability for therapeutic purposes.

59 BASIC BIOLOGICAL SCIENCES↗

Free Energy Landscapes for Elucidating the Structural Consequences of Exon-20 mutations on the ErbB Family of Protein Kinases

The ErbB family of protein kinases plays an important role in major cellular functions and consequently mutations in the functional regions of these proteins are implicated in several types of cancer growths. To envision rational design of small molecule drugs that target the diseased proteins it is important to quantify the structural effects of the mutations, as some of these mutants render the protein resistant to tyrosine kinase inhibitors (TKIs). Herein we use advanced sampling techniques and long-timescale molecular dynamics simulations to predict the effect of major exon 20 mutations on the ErbB family, specifically EGFR and HER2 proteins. Exon 20 mutations have been clinically known to induce TKI resistance, though the mechanisms of such an effect is poorly understood. By mapping out the free energy landscape of the mutants and comparing them against the wild-type, we elucidate the structural differences in the binding pocket region that alter the nature of drug-protein interactions. We believe that these insights will play a pivotal role in developing small molecule drugs that overcome the TKI resistance.

Ashwin Ravichandran↗

A multi-scale pipeline linking drug transcriptomics with pharmacokinetics predicts in vivo interactions of tuberculosis drugs

Tuberculosis (TB) is the deadliest infectious disease worldwide. The design of new treatments for TB is hindered by the large number of candidate drugs, drug combinations, dosing choices, and complex pharmaco-kinetics/dynamics (PK/PD). Here we study the interplay of these factors in designing combination therapies by linking a machine-learning model, INDIGO-MTB, which predicts in vitro drug interactions using drug transcriptomics, with a multi-scale model of drug PK/PD and pathogen-immune interactions called GranSim. We calculate an in vivo drug interaction score (iDIS) from dynamics of drug diffusion, spatial distribution, and activity within lesions against various pathogen sub-populations. The iDIS of drug regimens evaluated against non-replicating bacteria significantly correlates with efficacy metrics from clinical trials. Our approach identifies mechanisms that can amplify synergistic or mitigate antagonistic drug interactions in vivo by modulating the relative distribution of drugs. Our mechanistic framework enables efficient evaluation of in vivo drug interactions and optimization of combination therapies.

59 BASIC BIOLOGICAL SCIENCES↗

Developing a SARS-CoV-2 main protease binding prediction random forest model for drug repurposing for COVID-19 treatment

The coronavirus disease 2019 (COVID-19) global pandemic resulted in millions of people becoming infected with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus and close to seven million deaths worldwide. It is essential to further explore and design effective COVID-19 treatment drugs that target the main protease of SARS-CoV-2, a major target for COVID-19 drugs. In this study, machine learning was applied for predicting the SARS-CoV-2 main protease binding of Food and Drug Administration (FDA)-approved drugs to assist in the identification of potential repurposing candidates for COVID-19 treatment. Ligands bound to the SARS-CoV-2 main protease in the Protein Data Bank and compounds experimentally tested in SARS-CoV-2 main protease binding assays in the literature were curated. These chemicals were divided into training (516 chemicals) and testing (360 chemicals) data sets. To identify SARS-CoV-2 main protease binders as potential candidates for repurposing to treat COVID-19, 1188 FDA-approved drugs from the Liver Toxicity Knowledge Base were obtained. A random forest algorithm was used for constructing predictive models based on molecular descriptors calculated using Mold2 software. Model performance was evaluated using 100 iterations of fivefold cross-validations which resulted in 78.8% balanced accuracy. The random forest model that was constructed from the whole training dataset was used to predict SARS-CoV-2 main protease binding on the testing set and the FDA-approved drugs. Model applicability domain and prediction confidence on drugs predicted as the main protease binders discovered 10 FDA-approved drugs as potential candidates for repurposing to treat COVID-19. Our results demonstrate that machine learning is an efficient method for drug repurposing and, thus, may accelerate drug development targeting SARS-CoV-2.

Research & Experimental Medicine↗

Binding mechanism and binding free energy of amino acids and citrate to hydroxyapatite surfaces as a function of crystallographic facet, pH, and electrolytes

Hydroxyapatite (HAP) is the major mineral phase in bone and teeth. The interaction of individual amino acids and citrate ions with different crystallographic HAP surfaces has remained uncertain for decades, creating a knowledge gap to rationally design interactions with peptides, proteins, and drugs. In this contribution, we quantify the binding mechanisms and binding free energies of the 20 end-capped natural amino acids and citrate ions on the basal (001) and prismatic (010)/(020) planes of hydroxyapatite at pH values of 7 and 5 for the first time at the molecular scale. We utilized over 1500 steered molecular dynamics simulations with highly accurate potentials that reproduce surface and hydration energies of (hkl) hydroxyapatite surfaces at different pH values. Charged residues demonstrate a much higher affinity to HAP than charge-neutral species due to the formation of superficial ion pairs and ease of penetration into layers of water molecules on the mineral surface. Binding free energies range from 0 to -60 kJ/mol and were determined with ~10% uncertainty. The highest affinity was found for citrate, followed by Asp(-) and Glu(-), and followed after a gap by Arg(+), Lys(+), as well as by His(+) at pH 5. The (hkl)-specific area density of calcium ions, the protonation state of phosphate ions, and subsurface directional order of the ions in HAP lead to surface-specific binding patterns. Amino acids without ionic side groups exhibit weak binding, between -3 and 0 kJ/mol, due to difficulties to penetrate the first layer of water molecules on the apatite surfaces. We explain recognition processes that remained elusive in experiments, in prior simulations, discuss agreement with available data, and reconcile conflicting interpretations. In conclusion, the findings can serve as useful input for the design of peptides, proteins, and drug molecules for the modification of bone and teeth-related materials, as well as control of apatite mineralization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding interactions between biomolecules and two-dimensional nanomaterials using in silico microscopes

We report two-dimensional (2D) nanomaterials such as graphene are increasingly used in research and industry for various biomedical applications. Extensive experimental and theoretical studies have revealed that 2D nanomaterials are promising drug delivery vehicles, yet certain materials exhibit toxicity under biological conditions. So far, it is known that 2D nanomaterials possess strong adsorption propensities for biomolecules. To mitigate potential toxicity and retain favorable physical and chemical properties of 2D nanomaterials, it is necessary to explore the underlying mechanisms of interactions between biomolecules and nanomaterials for the subsequent design of biocompatible 2D nanomaterials for nanomedicine. The purpose of this review is to integrate experimental findings with theoretical observations and facilitate the study of 2D nanomaterial interaction with biomolecules at the molecular level. We discuss the current understanding and progress of 2D nanomaterial interaction with proteins, lipid membranes, and DNA based on molecular dynamics (MD) simulation. In this review, we focus on the 2D graphene nanosheet and briefly discuss other 2D nanomaterials. With the ever-growing computing power, we can image nanoscale processes using MD simulation that are otherwise not observable in experiment. We expect that molecular characterization of the complex behavior between 2D nanomaterials and biomolecules will help fulfill the goal of designing effective 2D nanomaterials as drug delivery platforms.

2D nanomaterials↗

Flow matching meets biology and life science: a survey

Over the past decade, advances in generative modeling, such as generative adversarial networks, masked autoencoders, and diffusion models, have significantly transformed biological research and discovery, enabling breakthroughs in molecule design, protein generation, catalysis discovery, drug discovery, and beyond. At the same time, biological applications have served as valuable testbeds for evaluating the capabilities of generative models. Recently, flow matching has emerged as a powerful and efficient alternative to diffusion-based generative modeling, with growing interest in its application to problems in biology and life sciences. This paper presents the first comprehensive survey of recent developments in flow matching and its applications in biological domains. We begin by systematically reviewing the foundations and variants of flow matching, and then categorize its applications into three major areas: biological sequence modeling, molecule generation and design, and peptide and protein generation. For each, we provide an in-depth review of recent progress. We also summarize commonly used datasets and software tools, and conclude with a discussion of potential future directions.

59 BASIC BIOLOGICAL SCIENCES↗