Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational docking”

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 19 records

Fragment binding to the Nsp3 macrodomain of SARS-CoV-2 identified through crystallographic screening and computational docking

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) macrodomain within the nonstructural protein 3 counteracts host-mediated antiviral adenosine diphosphate–ribosylation signaling. This enzyme is a promising antiviral target because catalytic mutations render viruses nonpathogenic. Here, we report a massive crystallographic screening and computational docking effort, identifying new chemical matter primarily targeting the active site of the macrodomain. Crystallographic screening of 2533 diverse fragments resulted in 214 unique macrodomain-binders. An additional 60 molecules were selected from docking more than 20 million fragments, of which 20 were crystallographically confirmed. X-ray data collection to ultra-high resolution and at physiological temperature enabled assessment of the conformational heterogeneity around the active site. Several fragment hits were confirmed by solution binding using three biophysical techniques (differential scanning fluorimetry, homogeneous time-resolved fluorescence, and isothermal titration calorimetry). The 234 fragment structures explore a wide range of chemotypes and provide starting points for development of potent SARS-CoV-2 macrodomain inhibitors.

59 BASIC BIOLOGICAL SCIENCES↗

Integration of Computational Docking into Anti-Cancer Drug Response Prediction Models

Cancer is a heterogeneous disease in that tumors of the same histology type can respond differently to a treatment. Anti-cancer drug response prediction is of paramount importance for both drug development and patient treatment design. Although various computational methods and data have been used to develop drug response prediction models, it remains a challenging problem due to the complexities of cancer mechanisms and cancer-drug interactions. To better characterize the interaction between cancer and drugs, we investigate the feasibility of integrating computationally derived features of molecular mechanisms of action into prediction models. Specifically, we add docking scores of drug molecules and target proteins in combination with cancer gene expressions and molecular drug descriptors for building response models. The results demonstrate a marginal improvement in drug response prediction performance when adding docking scores as additional features, through tests on large drug screening data. We discuss the limitations of the current approach and provide the research community with a baseline dataset of the large-scale computational docking for anti-cancer drugs.

60 APPLIED LIFE SCIENCES↗

Dynamic Docking Test System (DDTS) active table computer program NASA Advanced Docking System (NADS)

A computer program was developed to describe the three-dimensional motion of the Dynamic Docking Test System active table. The input consists of inertia and geometry data, actuator structural data, forcing function data, hydraulics data, servo electronics data, and integration control data. The output consists of table responses, actuator bending responses, and actuator responses.

Gates, R. M.↗

GPU-Accelerated Drug Discovery with Docking on the Summit Supercomputer: Porting, Optimization, and Application to COVID-19 Research

Protein-ligand docking is an in silico tool used to screen potential drug compounds for their ability to bind to a given protein receptor within a drug-discovery campaign. Experimental drug screening is expensive and time consuming, and it is desirable to carry out large scale docking calculations in a high-throughput manner to narrow the experimental search space. Few of the existing computational docking tools were designed with high performance computing in mind. Therefore, optimizations to maximize use of high-performance computational resources available at leadership-class computing facilities enables these facilities to be leveraged for drug discovery. Here we present the porting, optimization, and validation of the AutoDock-GPU program for the Summit supercomputer, and its application to initial compound screening efforts to target proteins of the SARS-CoV-2 virus responsible for the current COVID-19 pandemic.

LeGrand, Scott↗

Structure of 3-mercaptopropionic acid dioxygenase with a substrate analog reveals bidentate substrate binding at the iron center

Thiol dioxygenases are a subset of nonheme iron oxygenases that catalyze the formation of sulfinic acids from sulfhydryl-containing substrates and dioxygen. Among this class, cysteine dioxygenases (CDOs) and 3-mercaptopropionic acid dioxygenases (3MDOs) are the best characterized, and the mode of substrate binding for CDOs is well understood. However, the manner in which 3-mercaptopropionic acid (3MPA) coordinates to the nonheme iron site in 3MDO remains a matter of debate. A model for bidentate 3MPA coordination at the 3MDO Fe-site has been proposed on the basis of computational docking, whereas steady-state kinetics and EPR spectroscopic measurements suggest a thiolate-only coordination of the substrate. To address this gap in knowledge, we determined the structure of Azobacter vinelandii 3MDO (Av3MDO) in complex with the substrate analog and competitive inhibitor, 3-hydroxypropionic acid (3HPA). The structure together with DFT computational modeling demonstrates that 3HPA and 3MPA associate with iron as chelate complexes with the substrate-carboxylate group forming an additional interaction with Arg168 and the thiol bound at the same position as in CDO. A chloride ligand was bound to iron in the coordination site assigned as the O 2 -binding site. Supporting HYSCORE spectroscopic experiments were performed on the (3MPA/NO)-bound Av3MDO iron nitrosyl (S = 3/2) site. In combination with spectroscopic simulations and optimized DFT models, this work provides an experimentally verified model of the Av3MDO enzyme–substrate complex, effectively resolving a debate in the literature regarding the preferred substrate-binding denticity. These results elegantly explain the observed 3MDO substrate specificity, but leave unanswered questions regarding the mechanism of substrate-gated reactivity with dioxygen.

36 MATERIALS SCIENCE↗

Novel small molecule FGF 23 inhibitors increase serum phosphate and improve skeletal abnormalities in Hyp mice

We report that excess fibroblast growth factor 23 (FGF23) causes hereditary hypophosphatemic rickets, such as X-linked hypophosphatemia (XLH) and tumor induced osteomalacia (TIO). A small molecule that specifically binds to FGF23 to prevent activation of the FGFR/a-Klotho complex has potential advantages over the currently approved systemically administered FGF23 blocking antibody. Using structure-based drug design we previously identified ZINC13407541 (N-[[2-(2-phenylethenyl)cyclopenten-1-yl]methylidene]hydroxylamine) as a small molecule antagonist for FGF23. Additional structure-activity studies developed a series of ZINC13407541 analogues with enhanced drug-like properties. In this study, we tested in a pre-clinical Hyp mouse homologue of XLH a direct connect analogue (8n) [(E)-2-(4-(tert-butyl)phenyl)cyclopent-1-ene-1-carbaldehyde oxime] that exhibited the greatest stability in microsomal assays, and 13a [(E)-2-((E)-4-methylstyryl)benzaldehyde oxime] that exhibited increased in vitro potency. Using cryo-electron microscopy (Cryo-EM) structure and computational docking, we identified a key binding residue (Q156) of the FGF23 antagonists, ZINC13407541 and its analogues (8n and 13a) in the N-terminal domain of FGF23 protein. Site-directed mutagenesis and bimolecular fluorescence complementation (BiFC)-fluorescence resonance energy transfer (FRET) assay confirmed the binding site of these three antagonists. We found that pharmacological inhibition of FGF23 with either of these compounds blocked FGF23 signaling and increased serum phosphate and 1,25(OH)2D concentrations in Hyp mice. Long-term parenteral treatment with 8n or 13a also enhanced linear bone growth, increased mineralization of bone, and narrowed the growth plate in Hyp mice. The more potent 13a compound had greater therapeutic effects in Hyp mice. Further optimization of these FGF23 inhibitors may lead to versatile drugs to treat excess FGF23-mediated disorders.

60 APPLIED LIFE SCIENCES↗

Carbonic Anhydrase Variants Catalyze the Reduction of Dialkyl Ketones with High Enantioselectivity

Abstract Human carbonic anhydrase II (hCAII) naturally catalyzes the reaction between two achiral molecules—water and carbon dioxide—to yield the achiral product carbonic acid through a zinc hydroxide intermediate. We have previously shown that a zinc hydride, instead of a hydroxide, can be generated in this enzyme to create a catalyst for the reduction of aryl ketones. Dialkyl ketones are more challenging to reduce, and the enantioselective reduction of dialkyl ketones with two alkyl groups that are similar in size and electronic properties, is a particularly challenging transformation to achieve with high activity and selectivity. Here, we show that hCAII, as well as a double mutant of it, catalyzes the enantioselective reduction of dialkyl ketones with high yields and enantioselectivities, even when the two alkyl groups are similar in size. We also show that variants of hCAII catalyze the site‐selective reduction of one ketone over the other in an unsymmetrical aliphatic diketone. Computational docking of a dialkyl ketone to variants of hCAII containing the zinc hydride provides insights into the origins of the reactivity of various substrates and the high enantioselectivity of the transformations and show how a confined environment can control the enantioselectivity of an abiological intermediate.

Chen, Reichi↗

Specificity determinants revealed by the structure of glycosyltransferase Campylobacter concisus PglA

Abstract In selected Campylobacter species, the biosynthesis of N‐linked glycoconjugates via the pgl pathway is essential for pathogenicity and survival. However, most of the membrane‐associated GT‐B fold glycosyltransferases responsible for diversifying glycans in this pathway have not been structurally characterized which hinders the understanding of the structural factors that govern substrate specificity and prediction of resulting glycan composition. Herein, we report the 1.8 Å resolution structure of Campylobacter concisus PglA, the glycosyltransferase responsible for the transfer of N ‐acetylgalatosamine (GalNAc) from uridine 5′‐diphospho‐ N ‐acetylgalactosamine (UDP‐GalNAc) to undecaprenyl‐diphospho‐ N , N ′‐diacetylbacillosamine (UndPP‐diNAcBac) in complex with the sugar donor GalNAc. This study identifies distinguishing characteristics that set PglA apart within the GT4 enzyme family. Computational docking of the structure in the membrane in comparison to homologs points to differences in interactions with the membrane‐embedded acceptor and the structural analysis of the complex together with bioinformatics and site‐directed mutagenesis identifies donor sugar binding motifs. Notably, E113, conserved solely among PglA enzymes, forms a hydrogen bond with the GalNAc C6″‐OH. Mutagenesis of E113 reveals activity consistent with this role in substrate binding, rather than stabilization of the oxocarbenium ion transition state, a function sometimes ascribed to the corresponding residue in GT4 homologs. The bioinformatic analyses reveal a substrate‐specificity motif, showing that Pro281 in a substrate binding loop of PglA directs configurational preference for GalNAc over GlcNAc. This proline is replaced by a conformationally flexible glycine, even in distant homologs, which favor substrates with the same stereochemistry at C4, such as glucose. The signature loop is conserved across all Campylobacter PglA enzymes, emphasizing its importance in substrate specificity.

Vuksanovic, Nemanja↗

Discovery of multidrug efflux pump inhibitors with a novel chemical scaffold

Multidrug efflux is a major contributor to antibiotic resistance in Gram-negative bacterial pathogens. Inhibition of multidrug efflux pumps is a promising approach for reviving the efficacy of existing antibiotics. Previously, inhibitors targeting both the efflux transporter AcrB and the membrane fusion protein AcrA in the Escherichia coli AcrAB-TolC efflux pump were identified. In this work we use existing physicochemical property guidelines to generate a filtered library of compounds for computational docking. We then experimentally test the top candidate coumpounds using in vitro binding assays and in vivo potentiation assays in bacterial strains with controllable permeability barriers. We thus identify a new class of inhibitors of E. coli AcrAB-TolC. Six molecules with a shared scaffold were found to potentiate the antimicrobial activity of erythromycin and novobiocin in hyperporinated E. coli cells. Importantly, these six molecules were also active in wild-type strains of both Acinetobacter baumannii and Klebsiella pneumoniae, potentiating the activity of erythromycin and novobiocin up to 8-fold.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design, Synthesis, and Biological Activity of Novel Ornithine Decarboxylase (ODC) Inhibitors

We here describe the design, synthesis, and biological activity of novel ornithine decarboxylase (ODC) inhibitors that show significantly higher potency in vitro than α-difluoromethylornithine (DFMO), a U.S. Food and Drug Administration (FDA) approved drug. We report two X-ray structures of ODC complexed with new ODC inhibitors, computational docking, molecular dynamics, and binding free energy calculations to validate the experimental models. The X-ray structures reveal that covalent adducts with pyridoxal phosphate (PLP) are formed in the active site of the human ODC enzyme, as verified by their preparation and enzymatic testing. Finally, we verified that the cellular activity of endogenous ODC was inhibited, and polyamine levels were reduced. Given that ODC is a clinically validated target, combined with the fact that DFMO is currently the only ODC inhibitor in clinical use for several indications, the further development of more potent ODC inhibitors with superior activity and physical properties is warranted.

60 APPLIED LIFE SCIENCES↗

Host receptor-targeted therapeutic approach to counter pathogenic New World mammarenavirus infections

Five New World mammarenaviruses (NWMs) cause life-threatening hemorrhagic fever (HF). Cellular entry by these viruses is mediated by human transferrin receptor 1 (hTfR1). Here, we demonstrate that an antibody (ch128.1/IgG1) which binds the apical domain of hTfR1, potently inhibits infection of attenuated and pathogenic NWMs in vitro. Computational docking of the antibody Fab crystal structure onto the known structure of hTfR1 shows an overlapping receptor-binding region shared by the Fab and the viral envelope glycoprotein GP1 subunit that binds hTfR1, and we demonstrate competitive inhibition of NWM GP1 binding by ch128.1/IgG1 as the principal mechanism of action. Importantly, ch128.1/IgG1 protects hTfR1-expressing transgenic mice against lethal NWM challenge. Additionally, the antibody is well-tolerated and only partially reduces ferritin uptake. Our findings provide the basis for the development of a novel, host receptor-targeted antibody therapeutic broadly applicable to the treatment of HF of NWM etiology.

59 BASIC BIOLOGICAL SCIENCES↗

An artificial intelligence accelerated virtual screening platform for drug discovery

Abstract Structure-based virtual screening is a key tool in early drug discovery, with growing interest in the screening of multi-billion chemical compound libraries. However, the success of virtual screening crucially depends on the accuracy of the binding pose and binding affinity predicted by computational docking. Here we develop a highly accurate structure-based virtual screen method, RosettaVS, for predicting docking poses and binding affinities. Our approach outperforms other state-of-the-art methods on a wide range of benchmarks, partially due to our ability to model receptor flexibility. We incorporate this into a new open-source artificial intelligence accelerated virtual screening platform for drug discovery. Using this platform, we screen multi-billion compound libraries against two unrelated targets, a ubiquitin ligase target KLHDC2 and the human voltage-gated sodium channel Na V 1.7. For both targets, we discover hit compounds, including seven hits (14% hit rate) to KLHDC2 and four hits (44% hit rate) to Na V 1.7, all with single digit micromolar binding affinities. Screening in both cases is completed in less than seven days. Finally, a high resolution X-ray crystallographic structure validates the predicted docking pose for the KLHDC2 ligand complex, demonstrating the effectiveness of our method in lead discovery.

Science & Technology - Other Topics↗

Structure of saguaro cactus virus 3′ translational enhancer mimics 5′ cap for eIF4E binding

The genomes of several plant viruses contain RNA structures at their 3' ends called cap-independent translation enhancers (CITEs) that bind the host protein factors such as mRNA 5' cap-binding protein eIF4E for promoting cap-independent genome translation. However, the structural basis of such 5' cap-binding protein recognition by the uncapped RNA remains largely unknown. Here, we have determined the crystal structure of a 3' CITE, panicum mosaic virus-like translation enhancer (PTE) from the saguaro cactus virus (SCV), using a Fab crystallization chaperone. The PTE RNA folds into a three-way junction architecture with a pseudoknot between the purine-rich R domain and pyrimidine-rich Y domain, which organizes the overall structure to protrude out a specific guanine nucleotide, G18, from the R domain that comprises a major interaction site for the eIF4E binding. The superimposable crystal structures of the wild-type, G18A, G18C, and G18U mutants suggest that the PTE scaffold is preorganized with the flipped-out G18 ready to dock into the eIF4E 5' cap-binding pocket. The binding studies with wheat and human eIF4Es using gel electrophoresis and isothermal titration calorimetry, and molecular docking computation for the PTE–eIF4E complex demonstrated that the PTE structure essentially mimics the mRNA 5' cap for eIF4E binding. Such 5' cap mimicry by the uncapped and structured viral RNA highlights how viruses can exploit RNA structures to mimic the host protein-binding partners and bypass the canonical mechanisms for their genome translation, providing opportunities for a better understanding of virus-host interactions and non-canonical translation mechanisms found in many pathogenic RNA viruses.

3' cap-independent translation enhancers↗

CryoEM structures of the human CLC-2 voltage-gated chloride channel reveal a ball-and-chain gating mechanism

CLC-2 is a voltage-gated chloride channel that contributes to electrical excitability and ion homeostasis in many different tissues. Among the nine mammalian CLC homologs, CLC-2 is uniquely activated by hyperpolarization, rather than depolarization, of the plasma membrane. The molecular basis for the divergence in polarity of voltage gating among closely related homologs has been a long-standing mystery, in part because few CLC channel structures are available. Here, we report cryoEM structures of human CLC-2 at 2.46 – 2.76 Å, in the presence and absence of the selective inhibitor AK-42. AK-42 binds within the extracellular entryway of the Cl – -permeation pathway, occupying a pocket previously proposed through computational docking studies. In the apo structure, we observed two distinct conformations involving rotation of one of the cytoplasmic C-terminal domains (CTDs). In the absence of CTD rotation, an intracellular N-terminal 15-residue hairpin peptide nestles against the TM domain to physically occlude the Cl – -permeation pathway. This peptide is highly conserved among species variants of CLC-2 but is not present in other CLC homologs. Previous studies suggested that the N-terminal domain of CLC-2 influences channel properties via a “ball-and-chain” gating mechanism, but conflicting data cast doubt on such a mechanism, and thus the structure of the N-terminal domain and its interaction with the channel has been uncertain. Through electrophysiological studies of an N-terminal deletion mutant lacking the 15-residue hairpin peptide, we support a model in which the N-terminal hairpin of CLC-2 stabilizes a closed state of the channel by blocking the cytoplasmic Cl – -permeation pathway.

59 BASIC BIOLOGICAL SCIENCES↗

Performance Portability of Molecular Docking Miniapp On Leadership Computing Platforms

Rapidly changing computer architectures, such as those found at high-performance computing (HPC) facilities, present the need for mini-applications (miniapps) that capture essential algorithms used in large applications to test program performance and portability, aiding transitions to new systems. The COVID-19 pandemic has fueled a flurry of activity in computational drug discovery, including the use of supercomputers and GPU acceleration for massive virtual screens for therapeutics. Recent work targeting COVID-19 at the Oak Ridge Leadership Computing Facility (OLCF) used the GPU-accelerated program AutoDock-GPU to screen billions of compounds on the Summit supercomputer. In this paper we present the development of a new miniapp, miniAutoDock-GPU, that can be used to evaluate the performance and portability of GPU-accelerated protein-ligand docking programs on different computer architectures. These tests are especially relevant as facilities transition from petascale systems and prepare for upcoming exascale systems that will use a variety of GPU vendors. The key calculations, namely, the Lamarckian genetic algorithm combined with a local search using a Solis-Wets based random optimization algorithm, are implemented. We developed versions of the miniapp using several different programming models for GPU acceleration, including a version using the CUDA runtime API for NVIDIA GPUs, and the Kokkos middle-ware API which is facilitated by C++ template libraries. A third version, currently in progress, uses the HIP programming model. These efforts will help facilitate the transition to exascale systems for this important emerging HPC application, as well as its use on a wide range of heterogeneous platforms.

Thavappiragasam, Mathialakan↗

Spacecraft flight simulation: A human factors investigation into the man-machine interface between an astronaut and a spacecraft performing docking maneuvers and other proximity operations

The anticipated increase in rendezvous and docking activities in the various space programs in the Space Station era necessitates a renewed interest in manual docking procedures. Ten test subjects participated in computer simulated docking missions in which the influence of initial velocity was examined. All missions started from a resting position of 304.8 meters (1000 feet) along the space station's +V-bar axis. Test subjects controlled their vehicle with a translational hand controller and digital auto pilot which are both virtually identical to their space shuttle counterparts. While the 0.1 percent rule (range rate is equal to 0.1 percent of the range) used by space shuttle pilots is comfortably safe, it is revealed to be extremely inefficient in terms of time and not justifiable in terms of marginal safety. Time is worth money, not only because of training and launch costs, but because the sooner a pilot and spacecraft return from a mission, the sooner they can begin the next one. Inexperienced test subjects reduced the costs of simulated docking by close to a factor of 2 and achieved safe dockings in less than 4 percent of the time the baseline approach would entail. This reduction in time can be used to save lives in the event of an accident on orbit, and can tremendously reduce docking costs if fuel is produced from waste water on orbit.

Brody, Adam R.↗