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.

90 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Maritime Fuel Cell Generator Project: 2018 – 2023

This report summarizes activity in the Maritime Fuel Cell (MarFC) Generator Project from 2018 – 2023. FY 2018 saw the implementation of upgrades and repairs, making the unit more reliable and operator friendly. In FY2019 the team engaged the Scripps Institution of Oceanography (SIO) to use the MarFC to provide zero-emission shore power to the research vessel R/V Robert Gordon Sproul while in port at the Nimitz Marine Facility in San Diego, CA. In FY2020, the MarFC unit was shipped to San Diego, CA. A fueling contract with IGX was established to support MarFC operations at SIO, with renewable hydrogen provided by the California State University Los Angeles (CSULA) hydrogen station. The project team (Sandia, Cummins/Hydrogenics) provided training to SIO staff on the technical details of the MarFC, the safe use of hydrogen in general and the MarFC in particular. The first fueling of the unit at the SIO pier was successfully completed by IGX. The first powering of a vessel with fuel-cell shore power was conducted with the R/V Robert Gordon Sproul. While the mechanical systems (lights, AC, ventilation, hydraulics, pumps and cranes) were powered without incident, problems arose when powering the computer systems. Inspections carried out in FY2021 revealed the MarFC needed routine maintenance. Maintenance was performed and the unit was upgraded. The MarFC was turned on after the year pause, and initial test data on power levels and stability were collected. FY2022 was a year spent repairing, upgrading and testing the MarFC unit. Spikes in power and voltage were observed above 60 kW that could potentially extend below 60 kW with time. Such spikes could cause problems with the Sproul electrical systems. These age-related problems, the extended time for the Sproul spent in dry dock for scheduled upgrades, and the associated need to reschedule the vessel’s high-priority science missions made it no longer possible to deploy the unit at SIO. After due consideration, the decision was made by DOE, MARAD and the project team to cease the deployment, remove the MarFC from SIO, and not pursue further deployment activities. On December 2, 2022, the MarFC unit was removed from the Scripps Nimitz Marine Facility and shipped to Fridley, Michigan. The Cummins/Hydrogenics plan for the unit is to assess the condition of the MarFC subcomponents, and then use it as a training/learning system for technical employees new to hydrogen fuel-cell technology. After summarizing project activity from 2018 – 2023, this report provides a review of lessons learned. This report provides next steps in contemplating a follow-on project that would further advance the use of fuel-cell-based shore power in a marine setting. A comparison is made of the project results to the original objectives. This report ends with an accounting of presentations stemming from the project, and a list of references.

08 HYDROGEN↗
DOE OSTI2021

Rational Enzyme Design for Controlled Functionalization of Acetylated Xylan for Cell-Free Polymer Biosynthesis

Xylan O-acetyltransferase 1 (XOAT1) is involved in O-acetylating the backbone of hemicellulose xylan. Recent structural analysis of XOAT1 showed two unequal lobes forming a cleft that is predicted to accommodate and position xylan acceptors into proximity with the catalytic triad. Here, we used docking and molecular dynamics simulations to investigate the optimal orientation of xylan in the binding cleft of XOAT1 and identify putative key residues (Gln445 and Arg444 on Minor lobe & Asn312, Met311 and Asp403 on Major lobe) involved in substrate interactions. Site-directed mutagenesis coupled with biochemical analyses revealed the major lobe of XOAT1 is important for xylan binding. Mutation of single key residues yielded XOAT1 variants with various enzymatic efficiencies that are applicable to one-pot synthesis of xylan polymers with different degrees of O-acetylation. Taken together, our results demonstrate the effectiveness of computational modeling in guiding enzyme engineering aimed at modulating xylan and redesigning plant cell walls.

93 SC - Biological and Environmental Research (BER↗
DOE OSTI2023

Design of Diverse Asymmetric Pockets in De Novo Homo-oligomeric Proteins

A challenge for design of protein–small-molecule recognition is that incorporation of cavities with size, shape, and composition suitable for specific recognition can considerably destabilize protein monomers. This challenge can be overcome through binding pockets formed at homo-oligomeric interfaces between folded monomers. Interfaces surrounding the central homo-oligomer symmetry axes necessarily have the same symmetry and so may not be well suited to binding asymmetric molecules. To enable general recognition of arbitrary asymmetric substrates and small molecules, we developed an approach to designing asymmetric interfaces at off-axis sites on homo-oligomers, analogous to those found in native homo-oligomeric proteins such as glutamine synthetase. We symmetrically dock curved helical repeat proteins such that they form pockets at the asymmetric interface of the oligomer with sizes ranging from several angstroms, appropriate for binding a single ion, to up to more than 20 Å across. Of the 133 proteins tested, 84 had soluble expression in E. coli, 47 had correct oligomeric states in solution, 35 had small-angle X-ray scattering (SAXS) data largely consistent with design models, and 8 had negative-stain electron microscopy (nsEM) 2D class averages showing the structures coming together as designed. Both an X-ray crystal structure and a cryogenic electron microscopy (cryoEM) structure are close to the computational design models. The nature of these proteins as homo-oligomers allows them to be readily built into higher-order structures such as nanocages, and the asymmetric pockets of these structures open rich possibilities for small-molecule binder design free from the constraints associated with monomer destabilization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Computational studies reveal Fluorine based quinolines to be potent inhibitors for proteins involved in SARS-CoV-2 assembly

World is witnessing one of the worst pandemics of this century caused by SARS-CoV-2 virus which has affected millions of individuals. Despite rapid efforts to develop vaccines and drugs for COVID-19, the disease is still not under control. Chloroquine (CQ) and Hydroxychloroquine (HCQ) are two very promising inhibitors which have shown positive effect in combating the disease in preliminary experimental studies, but their use was reduced due to severe side-effects. Here, we performed a theoretical investigation of the same by studying the binding of the molecules with SARS-COV-2 Spike protein, the complex formed by Spike and ACE2 human receptor and a human serine protease TMPRSS2 which aids in cleavage of the Spike protein to initiate the viral activation in the body. Both the molecules had shown very good docking energies in the range of -6kcal/mol. Subsequently, we did a high throughput screening for other potential quinoline candidates which could be used as inhibitors. From the large pool of ligand candidates, we shortlisted the top three ligands (binding energy -8kcal/mol). We tested the stability of the docked complexes by running Molecular Dynamics (MD) simulations where we observed the stability of the quinoline analogues with the Spike-ACE2 and TMPRSS2 nevertheless the quinolines were not stable with the Spike protein alone. Thus, although the inhibitors bond very well with the protein molecules their intrinsic binding affinity depends on the protein dynamics. Moreover, the quinolines were stable when bound to electronegative pockets of Spike-ACE2 or TMPRSS2 but not with Viral Spike protein. We also observed that a Fluoride based compound: 3-[3-(Trifluoromethyl)phenyl]quinoline helps the inhibitor to bind with both Spike-ACE2 and TMPRSS2 with equal probability. The molecular details presented in this study would be very useful for developing quinoline based drugs for COVID-19 treatment.

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

Structural Characterization of TRAF6 N-Terminal for Therapeutic Uses and Computational Studies on New Derivatives

Tumor necrosis factor receptor-associated factors (TRAFs) are a protein family with a wide variety of roles and binding partners. Among them, TRAF6, a ubiquitin ligase, possesses unique receptor binding specificity and shows diverse functions in immune system regulation, cellular signaling, central nervous system, and tumor formation. TRAF6 consists of an N-terminal Really Interesting New Gene (RING) domain, multiple zinc fingers, and a C-terminal TRAF domain. TRAF6 is an important therapeutic target for various disorders and structural studies of this protein are crucial for the development of next-generation therapeutics. Here, we presented a TRAF6 N-terminal structure determined at the Turkish light source “Turkish DeLight” to be 3.2 Å resolution at cryogenic temperature (PDB ID: 8HZ2). This structure offers insight into the domain organization and zinc-binding, which are critical for protein function. Since the RING domain and the zinc fingers are key targets for TRAF6 therapeutics, structural insights are crucial for future research. Separately, we rationally designed numerous new compounds and performed molecular docking studies using this template (PDB ID:8HZ2). According to the results, 10 new compounds formed key interactions with essential residues and zinc ion in the N-terminal region of TRAF6. Molecular dynamic (MD) simulations were performed for 300 ns to evaluate the stability of three docked complexes (compounds 256, 322, and 489). Compounds 256 and 489 was found to possess favorable bindings with TRAF6. These new compounds also showed moderate to good pharmacokinetic profiles, making them potential future drug candidates as TRAF6 inhibitors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

On the Rapid Calculation of Binding Affinities for Antigen and Antibody Design and Affinity Maturation Simulations

The accurate and efficient calculation of protein-protein binding affinities is an essential component in antibody and antigen design and optimization, and in computer modeling of antibody affinity maturation. Such calculations remain challenging despite advances in computer hardware and algorithms, primarily because proteins are flexible molecules, and thus, require explicit or implicit incorporation of multiple conformational states into the computational procedure. The astronomical size of the amino acid sequence space further compounds the challenge by requiring predictions to be computed within a short time so that many sequence variants can be tested. In this study, we compare three classes of methods for antibody/antigen (Ab/Ag) binding affinity calculations: (i) a method that relies on the physical separation of the Ab/Ag complex in equilibrium molecular dynamics (MD) simulations, (ii) a collection of 18 scoring functions that act on an ensemble of structures created using homology modeling software, and (iii) methods based on the molecular mechanics-generalized Born surface area (MM-GBSA) energy decomposition, in which the individual contributions of the energy terms are scaled to optimize agreement with the experiment. When applied to a set of 49 antibody mutations in two Ab/HIV gp120 complexes, all of the methods are found to have modest accuracy, with the highest Pearson correlations reaching about 0.6. In particular, the most computationally intensive method, i.e., MD simulation, did not outperform several scoring functions. The optimized energy decomposition methods provided marginally higher accuracy, but at the expense of requiring experimental data for parametrization. Within each method class, we examined the effect of the number of independent computational replicates, i.e., modeled structures or reinitialized MD simulations, on the prediction accuracy. We suggest using about ten modeled structures for scoring methods, and about five simulation replicates for MD simulations as a rule of thumb for obtaining reasonable convergence. We anticipate that our study will be a useful resource for practitioners working to incorporate binding affinity calculations within their protein design and optimization process.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

PTM‐Psi : A python package to facilitate the computational investigation of p ost‐ t ranslational m odification on p rotein s tructures and their i mpacts on dynamics and functions

Abstract Post‐translational modification (PTM) of a protein occurs after it has been synthesized from its genetic template, and involves chemical modifications of the protein's specific amino acid residues. Despite of the central role played by PTM in regulating molecular interactions, particularly those driven by reversible redox reactions, it remains challenging to interpret PTMs in terms of protein dynamics and function because there are numerous combinatorially enormous means for modifying amino acids in response to changes in the protein environment. In this study, we provide a workflow that allows users to interpret how perturbations caused by PTMs affect a protein's properties, dynamics, and interactions with its binding partners based on inferred or experimentally determined protein structure. This Python‐based workflow, called PTM‐Psi , integrates several established open‐source software packages, thereby enabling the user to infer protein structure from sequence, develop force fields for non‐standard amino acids using quantum mechanics, calculate free energy perturbations through molecular dynamics simulations, and score the bound complexes via docking algorithms. Using the S ‐nitrosylation of several cysteines on the GAP2 protein as an example, we demonstrated the utility of PTM‐Psi for interpreting sequence–structure–function relationships derived from thiol redox proteomics data. We demonstrate that the S ‐nitrosylated cysteine that is exposed to the solvent indirectly affects the catalytic reaction of another buried cysteine over a distance in GAP2 protein through the movement of the two ligands. Our workflow tracks the PTMs on residues that are responsive to changes in the redox environment and lays the foundation for the automation of molecular and systems biology modeling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

HD-Bind: Encoding of Molecular Structure with Low Precision, Hyperdimensional Binary Representations

Publicly available collections of drug-like molecules have grown to comprise tens of billions of compounds due to advances in combinatorial chemistry. Traditional methods for identifying "hit" molecules from a large collection of potential drug-like candidates have relied on biophysical theory to compute approximations to the Gibbs free energy of the binding interaction between the drug and its protein target. These approaches have the major drawback that they require exceptional computing capabilities for even relatively small collections of molecules. Hyperdimensional Computing (HDC) is a recently-proposed learning paradigm that represents data with high-dimension binary vectors; this allows the use of low-precision binary vector arithmetic to create models of the data that can be learned without the need for the gradient-based optimization required in many conventional machine learning and deep learning methods. This algorithmic simplicity allows for acceleration in hardware that has been previously demonstrated in a range of application areas. We consider existing HDC approaches for molecular property classification and introduce two novel encodings of a commonly-used molecular representation, the extended connectivity fingerprint (ECFP). We show that HDC-based inference methods are as much as 91 times more efficient than traditional machine learning methods, and achieve an acceleration of nearly nine orders of magnitude compared to molecular docking. Our results show that HDC accelerated methods retain competitive accuracy on a number of well-studied tasks such as molecular property predictions using the MoleculeNet dataset, and bind/no-bind activity classification using the DUD-E and LIT-PCBA datasets. Our work thus motivates further investigation into molecular representation learning to develop ultraefficient pre-screening tools.

Jones, William↗
DOE OSTI2021

Structure and assembly of the diiron cofactor in the heme-oxygenase–like domain of the N-nitrosourea–producing enzyme SznF

In biosynthesis of the pancreatic cancer drug streptozotocin, the tridomain nonheme-iron oxygenase SznF hydroxylates N δ and N ω ' of N ω -methyl- L -arginine before oxidatively rearranging the triply modified guanidine to the N-methyl-N-nitrosourea pharmacophore. A previously published structure visualized the monoiron cofactor in the enzyme’s C-terminal cupin domain, which promotes the final rearrangement, but exhibited disorder and minimal metal occupancy in the site of the proposed diiron cofactor in the N-hydroxylating heme-oxygenase–like (HO-like) central domain. We leveraged our recent observation that the N-oxygenating µ-peroxodiiron(III/III) intermediate can form in the HO-like domain after the apo protein self-assembles its diiron(II/II) cofactor to solve structures of SznF with both of its iron cofactors bound. These structures of a biochemically validated member of the emerging heme-oxygenase–like diiron oxidase and oxygenase (HDO) superfamily with intact diiron cofactor reveal both the large-scale conformational change required to assemble the O 2 -reactive Fe 2 (II/II) complex and the structural basis for cofactor instability—a trait shared by the other validated HDOs. During cofactor (dis)assembly, a ligand-harboring core helix dynamically (un)folds. The diiron cofactor also coordinates an unanticipated Glu ligand contributed by an auxiliary helix implicated in substrate binding by docking and molecular dynamics simulations. The additional carboxylate ligand is conserved in another N-oxygenating HDO but not in two HDOs that cleave carbon–hydrogen and carbon–carbon bonds to install olefins. Among ~9,600 sequences identified bioinformatically as members of the emerging HDO superfamily, ~25% conserve this additional carboxylate residue and are thus tentatively assigned as N-oxygenases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Targeting Aspergillus allergen oryzin with a chemical probe at atomic precision

We report the molecular basis of Aspergillus fumigatus oryzin, allergen Asp f 13, or alkaline proteinase ALP1, containing the sequence motif His–Asp–Ser of the subtilisin family, structure, and function at atomic detail. Given the resolution of the data (1.06 Å), we use fragment molecular replacement with ideal polyalanine α-helices to determine the first crystal structure of oryzin. We probe the catalytic serine through formation of an irreversible bond to a small molecule compound, specifically labeling it, describing the amino acid residues performing the catalytic function. Defined by a self-processed pro-peptide, the active site architecture shapes up pocket-like subsites that bind to and unveil the S1′–S4′ substrate binding preferences. We use molecular modeling to dock a model of the pro-peptide in the S1–S4 region and to dock collagen along the active site cleft. Opposite to the face harboring the catalytic serine, the enzyme binds to a calcium ion in a binding site created by backbone flipping. We use thermal unfolding to show that this metal ion provides structural stability. With no known host inhibitor identified thus far, this structure may hasten the progress of developing new therapeutic agents for diseases caused by pathogenic fungi.

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

High-throughput virtual laboratory for drug discovery using massive datasets

Time-to-solution for structure-based screening of massive chemical databases for COVID-19 drug discovery has been decreased by an order of magnitude, and a virtual laboratory has been deployed at scale on up to 27,612 GPUs on the Summit supercomputer, allowing an average molecular docking of 19,028 compounds per second. Over one billion compounds were docked to two SARS-CoV-2 protein structures with full optimization of ligand position and 20 poses per docking, each in under 24 hours. GPU acceleration and high-throughput optimizations of the docking program produced 350× mean speedup over the CPU version (50× speedup per node). GPU acceleration of both feature calculation for machine-learning based scoring and distributed database queries reduced processing of the 2.4 TB output by orders of magnitude. The resulting 50× speedup for the full pipeline reduces an initial 43 day runtime to 21 hours per protein for providing high-scoring compounds to experimental collaborators for validation assays.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Integrated data-driven and experimental approaches to accelerate lead optimization targeting SARS-CoV- 2 main protease

Identification of potential therapeutic candidates can be expedited by integrating computational modeling with domain aware machine learning (ML) approaches followed by experimental validation. Generative deep learning models have been recently developed that can generate thousands of new candidates, but their physiochemical properties are typically not optimized. Using our deep learning models and a scaffold as a starting point, we generated tens of thousands of compounds for SARS-CoV-2 M pro that preserve the core scaffold. Here we utilized and implemented several computational tools such as structural alert and toxicity analysis, high throughput virtual screening, ML-based 3D quantitative structure–activity relationships, multi-parameter optimization, and graph neural networks on libraries of generated candidates to predict biological activity and binding affinity a priori. From these collective computational results, eight promising candidates were identified and tested experimentally using Native Mass Spectrometry (MS) and FRET-based functional assays. Two compounds, with quinazoline-2-thiol and acetylpiperidine core moiety showed IC 50 values in the low micromolar range: 2.95±0.0017 µM and 3.41±0.0015 µM, respectively. The molecular dynamics simulations further highlight that binding of these compounds results in allosteric modulations in the chain B and the interface domains of the M pro . The key fragments from these top hits can be used as input for closed loop lead optimization in the integrated pipeline.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Full Scale 3D Computational Model of the Industrial -Scale Coal Fired Boiler Performance for Temperature Sensor Installation Guidance

Abstract Nearly 30% of the electricity is generated by using coal as the primary fuel in the US. One of the major concerns in coal-fired power plants is the failure of boiler tubes that leads to unscheduled maintenance and has a huge economical and societal impact. High temperature flue gas along with ash pass over the boiler tubes, which over time leads to tube failure. Therefore, developing temperature sensors for harsh environments and install them for temperature sensing and boiler tube lifetime prediction is an urgent need. On the side of sensor development, the location of the sensor installation is important for stable sensing performance and easy calibration. In this study, computational fluid dynamics and heat transfer modeling are adopted to establish a full-scale 3-dimensional model of a coal-fired boiler to investigate the flue gas temperature distribution within the boiler and identify the proper locations for sensor installation. We proposed three criteria to select the temperature sensor installation location: (1) select the boiler tube panel away from the sidewalls, (2) select the boiler tube section closer to the top wall of the boiler; and (3) select the boiler tube on the back of the boiler panel (not directly facing the flue gas flow). In these regions, the flue gas temperature is stable, providing an ideal environment for stable temperature sensing and calibration.

Gupta, Tanuj↗
DOE OSTI2020

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

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.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

tinyIFD: A High-Throughput Binding Pose Refinement Workflow Through Induced-Fit Ligand Docking

A critical step in structure-based drug discovery is predicting whether and how a candidate molecule binds to a model of a therapeutic target. However, substantial protein side chain movements prevent current screening methods, such as docking, from accurately predicting the ligand conformations and require expensive refinements to produce viable candidates. Here, we present the development of a high-throughput and flexible ligand pose refinement workflow, called “tinyIFD”. The main features of the workflow include the use of specialized high-throughput, small-system MD simulation code mdgx.cuda and an actively learning model zoo approach. We show the application of this workflow on a large test set of diverse protein targets, achieving 66% and 76% success rates for finding a crystal-like pose within the top-2 and top-5 poses, respectively. We also applied this workflow to the SARS-CoV-2 main protease (M pro ) inhibitors, where we demonstrate the benefit of the active learning aspect in this workflow.

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

Comparative Analysis of TCR and TCR-pMHC Complex Structure Prediction Tools

The rapid development of computational approaches for predicting the structures of T cell receptors (TCRs) and TCR-peptide-major histocompatibility (TCR-pMHC) complexes, accelerated by AI breakthroughs such as AlphaFold, has made it feasible to calculate these structures with increasing accuracy. Although these tools show great potential, their relative accuracy and limitations remain unclear due to the lack of standardized benchmarks. Here, we systematically evaluate seven tools for predicting isolated TCR structures together with six tools for predicting TCR-pMHC complex structures. The methods include homology-based approaches, general prediction tools using AlphaFold, TCR-specific tools derived from AlphaFold2, and the newly developed tFold-TCR model. The evaluation uses a post-training data set comprising 40 αβ TCRs and 27 TCR-pMHC complexes (21 Class I and 6 Class II). Model accuracy is assessed at global, local, and interface levels using a variety of metrics. We find that each tool offers distinct advantages in various aspects of its predictions. AlphaFold2, AlphaFold3, and tFold-TCR excel in overall accuracy of TCR structure prediction, and TCRmodel2 and AlphaFold2 perform well in overall accuracy of TCR-pMHC structure prediction. However, TCR-specific tools derived from AlphaFold2 show lower accuracy in the framework region than both homology-based methods and general-purpose tools such as AlphaFold, and challenges remain for all in modeling CDR3 loops, docking orientations, TCR-peptide interfaces, and Class II MHC-peptide interfaces. Furthermore, these findings will guide researchers in selecting appropriate tools, emphasize the importance of using multiple evaluation metrics to assess model performance, and offer suggestions for improving TCR and TCR-pMHC structure prediction tools.

Chemical structure↗
DOE OSTI2024

Predicting permeation of compounds across the outer membrane of P. aeruginosa using molecular descriptors

The ability Gram-negative pathogens have at adapting and protecting themselves against antibiotics has increasingly become a public health threat. Data-driven models identifying molecular properties that correlate with outer membrane (OM) permeation and growth inhibition while avoiding efflux could guide the discovery of novel classes of antibiotics. Here we evaluate 174 molecular descriptors in 1260 antimicrobial compounds and study their correlations with antibacterial activity in Gram-negative Pseudomonas aeruginosa. The descriptors are derived from traditional approaches quantifying the compounds’ intrinsic physicochemical properties, together with, bacterium-specific from ensemble docking of compounds targeting specific MexB binding pockets, and all-atom molecular dynamics simulations in different subregions of the OM model. Using these descriptors and the measured inhibitory concentrations, we design a statistical protocol to identify predictors of OM permeation/inhibition. We find consistent rules across most of our data highlighting the role of the interaction between the compounds and the OM. An implementation of the rules uncovered in our study is shown, and it demonstrates the accuracy of our approach in a set of previously unseen compounds. Our analysis sheds new light on the key properties drug candidates need to effectively permeate/inhibit P. aeruginosa, and opens the gate to similar data-driven studies in other Gram-negative pathogens.

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

AI-Accelerated Design of Targeted Covalent Inhibitors for SARS-CoV-2

Direct-acting antivirals for the treatment of the COVID-19 pandemic caused by the SARS-CoV-2 virus are needed to complement vaccination efforts. Given the ongoing emergence of new variants, automated experimentation, and active learning based fast workflows for antiviral lead discovery remain critical to our ability to address the pandemic’s evolution in a timely manner. While several such pipelines have been introduced to discover candidates with noncovalent interactions with the main protease (M pro ), here we developed a closed-loop artificial intelligence pipeline to design electrophilic warhead-based covalent candidates. Here, this work introduces a deep learning-assisted automated computational workflow to introduce linkers and an electrophilic “warhead” to design covalent candidates and incorporates cutting-edge experimental techniques for validation. Using this process, promising candidates in the library were screened, and several potential hits were identified and tested experimentally using native mass spectrometry and fluorescence resonance energy transfer (FRET)-based screening assays. We identified four chloroacetamide-based covalent inhibitors of M pro with micromolar affinities (K I of 5.27 μM) using our pipeline. Experimentally resolved binding modes for each compound were determined using room-temperature X-ray crystallography, which is consistent with the predicted poses. The induced conformational changes based on molecular dynamics simulations further suggest that the dynamics may be an important factor to further improve selectivity, thereby effectively lowering KI and reducing toxicity. These results demonstrate the utility of our modular and data-driven approach for potent and selective covalent inhibitor discovery and provide a platform to apply it to other emerging targets.

60 APPLIED LIFE SCIENCES↗