Search NASA⌕ Search

DOE OSTI · 1778008

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

Also available from

Abstract

In this work, we present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Acharya, Atanu, Agarwal, Rupesh, Baker, M. B., Baudry, Jerome Y., Bhowmik, Debsindhu, Boehm, Swen, Byler, Kendall G., Chen, S. Y., Coates, Leighton, Cooper, Connor J., Demerdash, Omar, Daidone, Isabella, Eblen, John D., Ellingson, Sally R., Forli, Stefano, Glaser, Jens, Gumbart, James C., Gunnels, John, Hernandez, Oscar, Irle, Stephan, Kneller, Daniel W., Kovalevsky, Andrii Y., Larkin, Jeffrey, Lawrence, Travis J., LeGrand, Scott, Liu, S.-H. Kevin, Mitchell, Julie C., Park, Gilchan, Parks, Jerry M., Pavlova, Anna, Petridis, Loukas, Poole, Duncan, Pouchard, Line, Ramanathan, Arvind, Rogers, David M., Santos-Martins, Diogo, Scheinberg, Aaron, Sedova, Ada, Shen, Yue, Smith, Jeremy C., Smith, Micholas D., Soto, Carlos, Tsaris, Aristeidis (Aris), Thavappiragasam, Mathialakan, Tillack, Andreas F., Vermaas, Josh V., Vuong, Van Quan, Yin, Junqi, Yoo, Shinjae, Zahran, Mai, Zanetti-Polzi, Laura. 2020-12-16. Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19. https://doi.org/10.1021/acs.jcim.0c01010

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

Metabolic flux and resource balance in the oleaginous yeast Rhodotorula toruloides

The yeast Rhodotorula toruloides is a promising bioproduction organism due to its high lipid yields and ability to grow on cheap and abundant substrates. Quantitative, systems-level assessment of its metabolic activity is accordingly merited. Resource-balance analysis (RBA) models capture not only reaction stoichiometry but also enzyme requirements for catalysis, providing valuable tools for understanding metabolic trade-offs and optimizing metabolic engineering strategies. Here, in this work, we present systems-level measurements of R. toruloides metabolic flux based on isotope tracing and metabolic flux analysis. In combination with new proteomic measurements, these flux data are used to parameterize a genome-scale resource balance model rtRBA. We find that S. cerevisiae and R. toruloides grow at nearly indistinguishable rates using similar biosynthetic but dramatically different central metabolic programs. R. toruloides consumes one-fifth as much glucose, which it metabolizes primarily via the pentose phosphate pathway and TCA cycle unlike primarily glycolysis in S. cerevisiae . Overall, across these two divergent yeasts, protein abundances aligned more closely than metabolic flux. Resource balance modeling of these metabolic programs predicts superior theoretical yields but lower productivities in R. toruloides than S. cerevisiae for industrial chemicals, highlighting the value of rapid glucose uptake for productivity but respiratory metabolism for yields.

60 APPLIED LIFE SCIENCES↗

Opportunities for bioenergy crops to support transitions from irrigated agriculture and conserve the U.S. High Plains Aquifer

This study investigates potential economic and groundwater driven transitions from irrigated maize production—the dominant irrigated cropping system in the High Plains Aquifer (HPA) region—to alternative crops such as sorghum and switchgrass, two common bioenergy feedstocks. Unsustainable groundwater extraction in the U.S. High Plains presents critical challenges including reduced irrigation capacities, diminished crop yields, lower land values, and escalating energy costs. Using a spatially explicit optimization framework combined with comprehensive economic and hydrogeological data, we evaluate optimal land-use strategies and irrigation system investments over a 30-year planning period. Results indicate significant regional variability in the future economic viability of irrigated agriculture, driven by differences in aquifer recharge rates, groundwater availability, and market conditions. Nebraska and parts of northern Texas can sustain irrigated maize profitability due to relatively favorable groundwater conditions and lower land rents, respectively. In contrast, many portions of Kansas and southern Texas are more likely to transition to dryland agriculture within two decades. Colorado and New Mexico show potential for significant adoption of switchgrass production as an alternative biomass-based energy crop. Overall, over the 30-year horizon, our model implies that approximately 23% of currently irrigated maize area across the HPA region may transition to dryland farming under business-as-usual conditions. This transition is complemented by a threefold increase in non-irrigated sorghum production, from 0.20 to 0.77 Mt yr−1, indicating that groundwater-driven shifts in agricultural production may support a larger regional base for bioenergy feedstocks. The study also reveals opportunities for producers to optimize economic returns and biomass production potential.

60 APPLIED LIFE SCIENCES↗