Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular 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

Elucidation of Agonist and Antagonist Dynamic Binding Patterns in ER-α by Integration of Molecular Docking, Molecular Dynamics Simulations and Quantum Mechanical Calculations

Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.

59 BASIC BIOLOGICAL SCIENCES↗

Portability for GPU-accelerated molecular docking applications for cloud and HPC: can portable compiler directives provide performance across all platforms?

High-throughput structure-based screening of drug-like molecules has become a common tool in biomedical research. Recently, acceleration with graphics processing units (GPUs) has provided a large performance boost for molecular docking programs. Both cloud and high-performance computing (HPC) resources have been used for large screens with molecular docking programs; while NVIDIA GPUs have dominated cloud and HPC resources, new vendors such as AMD and Intel are now entering the field, creating the problem of software portability across different GPUs. Ideally, software productivity could be maximized with portable programming models that are able to maintain high performance across architectures. While in many cases compiler directives have been used as an easy way to offload parallel regions of a CPU-based program to a GPU accelerator, they may also be an attractive programming model for providing portability across different GPU vendors, in which case the porting process may proceed in the reverse direction: from low-level, architecture-specific code to higher-level directive-based abstractions. MiniMDock is a new mini-application (miniapp) designed to capture the essential computational kernels found in molecular docking calculations, such as are used in phar-maceutical drug discovery efforts, in order to test different solutions for porting across GPU architectures. Here we extend MiniMDock to GPU offloading with OpenMP directives, and compare to performance of kernels using CUDA and HIP on NVIDIA and AMD GPUs, respectively, as well as across different compilers, exploring performance bottlenecks. We document this reverse-porting process, from highly optimized device code to a higher-level version using directives, compare code structure, and describe barriers that were overcome in this effort.

Thavappiragasam, Mathialakan↗

Resveratrol Ameliorates Lipopolysaccharide-Induced Sudden Sensorineural Hearing Loss in In Vitro Model through Multitarget Antiapoptotic Mechanism Based on Network Pharmacology and Molecular Docking

Objective. To explore the effects of resveratrol (RSV) on hair cell apoptosis caused by sudden sensorineural hearing loss (SSNHL) and its effect on lipopolysaccharide-induced apoptosis of HEI-OC1 cells. Methods. We used the network pharmacology method to screen molecules related to RSV for the treatment of SSNHL and analyzed these molecules and their enriched biological processes and signaling pathways through Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analysis. We selected hub genes related to apoptosis using protein-protein interaction (PPI) analysis for in vitro and molecular docking verification. Results. Eighty overlapping genes were identified as potential targets for RSV treatment of SSNHL. Further GO analysis showed that the biological processes were mainly related to toxicity, cell proliferation, and lipopolysaccharide reactions. KEGG analysis showed that the AGE-RAGE signaling pathway in diabetic complications, Kaposi’s sarcoma-associated herpesvirus infection, FoxO signaling pathway, PI3K-Akt signaling pathway, and other inflammatory signaling pathways were concentrated. AKT1, STAT3, JUN, TNF, TP53, MAPK3, CASP3, and VEGFA were screened as HUB genes using PPI analysis. The apoptosis-related proteins TNF, CASP3, AKT1, and TP53 were selected for in vitro experiments, which showed that mRNA was significantly different before and after RSV intervention, confirming that the corresponding protein receptors could bind well with RSV. Conclusion. RSV mainly affects the prognosis of SSNHL through anti-inflammatory effects and may improve hair cell apoptosis caused by inflammatory factors through multitargeted interventions involving TNF, CASP3, AKT1, and TP53.

Ye, Shiming↗

Identification of Small-Molecule Inhibitors of Fibroblast Growth Factor 23 Signaling via In Silico Hot Spot Prediction and Molecular Docking to α-Klotho

Fibroblast growth factor 23 (FGF23) is a therapeutic target for treating hereditary and acquired hypophosphatemic disorders, such as X-linked hypophosphatemic (XLH) rickets and tumor-induced osteomalacia (TIO), respectively. FGF23-induced hypophosphatemia is mediated by signaling through a ternary complex formed by FGF23, the FGF receptor (FGFR), and α-Klotho. Currently, disorders of excess FGF23 are treated with an FGF23-blocking antibody, burosumab. Small-molecule drugs that disrupt protein/protein interactions necessary for the ternary complex formation offer an alternative to disrupting FGF23 signaling. In this study, the FGF23:α-Klotho interface was targeted to identify small-molecule protein/protein interaction inhibitors since it was computationally predicted to have a large fraction of hot spots and two druggable residues on α-Klotho. We further identified Tyr433 on the KL1 domain of α-Klotho as a promising hot spot and α-Klotho as an appropriate drug-binding target at this interface. Subsequently, we performed in silico docking of ~5.5 million compounds from the ZINC database to the interface region of α-Klotho from the ternary crystal structure. Following docking, 24 and 20 compounds were in the final list based on the lowest binding free energies to α-Klotho and the largest number of contacts with Tyr433, respectively. Herein, five compounds were assessed experimentally by their FGF23-mediated extracellular signal-regulated kinase (ERK) activities in vitro, and two of these reduced activities significantly. Both these compounds were predicted to have favorable binding affinities to α-Klotho but not have a large number of contacts with the hot spot Tyr433. ZINC12409120 was found experimentally to disrupt FGF23:α-Klotho interaction to reduce FGF23-mediated ERK activities by 70% and have a half maximal inhibitory concentration (IC50) of 5.0 ± 0.23 μM. Molecular dynamics (MD) simulations of the ZINC12409120:α-Klotho complex starting from in silico docking poses reveal that the ligand exhibits contacts with residues on the KL1 domain, the KL1–KL2 linker, and the KL2 domain of α-Klotho simultaneously, thereby possibly disrupting the regular function of α-Klotho and impeding FGF23:α-Klotho interaction. ZINC12409120 is a candidate for lead optimization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Discovery of Small-Molecule Inhibitors of SARS-CoV-2 Proteins Using a Computational and Experimental Pipeline

A rapid response is necessary to contain emergent biological outbreaks before they can become pandemics. The novel coronavirus (SARS-CoV-2) that causes COVID-19 was first reported in December of 2019 in Wuhan, China and reached most corners of the globe in less than two months. In just over a year since the initial infections, COVID-19 infected almost 100 million people worldwide. Although similar to SARS-CoV and MERS-CoV, SARS-CoV-2 has resisted treatments that are effective against other coronaviruses. Crystal structures of two SARS-CoV-2 proteins, spike protein and main protease, have been reported and can serve as targets for studies in neutralizing this threat. We have employed molecular docking, molecular dynamics simulations, and machine learning to identify from a library of 26 million molecules possible candidate compounds that may attenuate or neutralize the effects of this virus. The viability of selected candidate compounds against SARS-CoV-2 was determined experimentally by biolayer interferometry and FRET-based activity protein assays along with virus-based assays. In the pseudovirus assay, imatinib and lapatinib had IC 50 values below 10 μM, while candesartan cilexetil had an IC 50 value of approximately 67 µM against M pro in a FRET-based activity assay. Comparatively, candesartan cilexetil had the highest selectivity index of all compounds tested as its half-maximal cytotoxicity concentration 50 (CC 50 ) value was the only one greater than the limit of the assay (>100 μM).

59 BASIC BIOLOGICAL SCIENCES↗

Crystal structure of dihydrofolate reductase from the filarial nematode W. bancrofti in complex with NADPH and folate

Lymphatic filariasis is a debilitating illness with an estimated 50 million cases as of 2018. The majority of cases are caused by the parasitic worm W. bancrofti and additional cases by the worms B. malayi and B. timori. Dihydrofolate reductase (DHFR) is an established target in the treatment of cancer, bacterial, and protozoal infections and may be a potential target for drugs targeting parasitic worm infections, including filariasis. Recent studies have shown that known antifolate compounds, including methotrexate, inhibit the activity of W. bancrofti DHFR (WbDHFR). However, the absence of structural information for filarial DHFRs has limited the study of more in-depth structure-function relationships. We report the structure of WbDHFR complexed with NADPH and folate using X-ray diffraction data measured to 2.47 Å resolution. The structure of WbDHFR reveals the usual DHFR fold and is currently only the second nematode DHFR structure in the Protein Data Bank. The equilibrium dissociation constants for NADPH (90 ± 29 nM) and folate (23 ± 4 nM) were determined by equilibrium titrations. The interactions of known antifolates with WbDHFR were analyzed using molecular docking programs and molecular dynamics simulations. Antifolates with a hydrophobic core and extended linker formed favorable interactions with WbDHFR. These combined data should now facilitate the rational design of filarial DHFR inhibitors, which in turn can be used to determine whether DHFR is a viable drug target for filariasis and whether existing antifolates may be repurposed for its treatment.

60 APPLIED LIFE SCIENCES↗

Paired Simulations and Experimental Investigations into the Calcium-Dependent Conformation of Albumin

Serum albumin is the most abundant protein in blood plasma, and it is involved in multiple biological processes. Serum albumin has recently been adapted for improving biomaterial integration with bone tissue, and studies have shown the importance of this protein in bone repair and regeneration. However, the mechanism of action is not yet clear. In stark contrast, other studies have demonstrated that albumin blocks cell adhesion to surfaces, which is seen as a limitation to its bone healing role. These apparent contradictions suggest that the conformation of albumin facilitates its bioactivity, leading to enhanced bone repair. Serum albumin is known to play a major role in maintaining the calcium ion concentration in blood plasma. Due to the prevalence of calcium at bone repair and regeneration sites, it has been hypothesized that calcium binding to serum albumin triggers a conformational change, leading to bioactivity. In the current study, molecular modeling approaches including molecular docking, atomic molecular dynamics (MD) simulation, and coarse-grained MD simulation were used to test this hypothesis by investigating the conformational changes induced in bovine serum albumin by interaction with calcium ions. The computational results were qualitatively validated with experimental Fourier-transform infrared spectroscopy analysis. We find that free calcium ions in solution transiently bind with the three major loops in albumin, triggering a conformational change where N-terminal and C-terminal domains separate from each other in a partial unfolding process. The separation distance between these domains was found to correlate with the calcium ion concentration. The experimental data support the simulation results showing that albumin has enhanced conformational heterogeneity upon exposure to intermediate levels of calcium, without any significant secondary structure changes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evolution of carbapenemase activity in the class C β-lactamase ADC-1

Antibiotic resistance in bacteria poses a significant threat to public health. Among dozens of available antimicrobial agents, carbapenems are used as drugs of choice for the treatment of serious infections caused by pathogens resistant to other antibiotics. However, their usefulness has been severely compromised due to the emergence and wide spread of carbapenem-resistant clinical isolates worldwide. High-level resistance to carbapenems in bacteria is mediated by the production of β-lactamases from three molecular classes, A, B, and D, but not by class C enzymes. In this study, we selected a triple mutant of the intrinsic class C Acinetobacter-derived cephalosporinase ADC-1 (ADC-1 TM ) that confers high-level resistance to the carbapenems meropenem, ertapenem, and doripenem. Kinetic experiments demonstrated that the apparent binding affinity, along with the acylation and deacylation rates, were all improved for the mutant enzyme. X-ray crystallography, molecular docking, and molecular dynamics simulations revealed that the amino acid substitutions in ADC-1 TM produce significant changes in the enzyme active site architecture and binding mode of the carbapenem ertapenem. These changes allow for better positioning of a deacylating water for nucleophilic attack, thus explaining the significantly improved rate of ertapenem deacylation by ADC-1 TM . In this study, we showed for the first time that a class C β-lactamase can produce high-level resistance to carbapenem antibiotics, which underlines the potential for enzymes of this class to evolve such resistance and could further exacerbate the problem of antibiotic resistance in bacteria.

59 BASIC BIOLOGICAL SCIENCES↗

C6 Hydroxymethyl-Substituted Carbapenem MA-1-206 Inhibits the Major Acinetobacter baumannii Carbapenemase OXA-23 by Impeding Deacylation

Acinetobacter baumannii has become a major nosocomial pathogen, as it is often multidrug-resistant, which results in infections characterized by high mortality rates. The bacterium achieves high levels of resistance to β-lactam antibiotics by producing β-lactamases, enzymes which destroy these valuable agents. Historically, the carbapenem family of b-lactam antibiotics have been the drugs of choice for treating A. baumannii infections. However, their effectiveness has been significantly diminished due to the pathogen’s production of carbapenem-hydrolyzing class D β-lactamases (CHDLs); thus, new antibiotics and inhibitors of these enzymes are urgently needed. Here, we describe a new carbapenem antibiotic, MA-1-206, in which the canonical C6 hydroxyethyl group has been replaced with hydroxymethyl. The antimicrobial susceptibility studies presented here demonstrated that this compound is more potent than meropenem and imipenem against A. baumannii producing OXA-23, the most prevalent CHDL of this pathogen, and also against strains producing the CHDL OXA-24/40 and the class B metallo-β-lactamase VIM-2. Our kinetic and mass spectrometry studies revealed that this drug is a reversible inhibitor of OXA-23, where inhibition takes place through a branched pathway. X-ray crystallographic studies, molecular docking, and molecular dynamics simulations of the OXA-23-MA-1-206 complex show that the C6 hydroxymethyl group forms a hydrogen bond with the carboxylated catalytic lysine of OXA-23, effectively preventing deacylation. These results provide a promising strategy for designing a new generation of CHDL-resistant carbapenems to restore their efficacy against deadly A. baumannii infections. Carbapenem antibiotics are the drugs of choice for treatment of deadly infections caused by Gram-negative bacteria. However, their efficacy is severely compromised by the wide spread of carbapenem-hydrolyzing class D β-lactamases (CHDLs). The importance of this research is the discovery that substitution of the canonical hydroxyethyl group of carbapenems by a hydroxymethyl significantly enhances stability against inactivation by the major CHDL of Acinetobacter baumannii, OXA-23. These results provide a novel strategy for designing next-generation, carbapenemase-stable carbapenems to fight multidrug-resistant infections caused by Gram-negative pathogens

59 BASIC BIOLOGICAL SCIENCES↗

Enhancing Docking Accuracy with PECAN2, a 3D Atomic Neural Network Trained without Co-Complex Crystal Structures

Decades of drug development research have explored a vast chemical space for highly active compounds. The exponential growth of virtual libraries enables easy access to billions of synthesizable molecules. Computational modeling, particularly molecular docking, utilizes physics-based calculations to prioritize molecules for synthesis and testing. Nevertheless, the molecular docking process often yields docking poses with favorable scores that prove to be inaccurate with experimental testing. To address these issues, several approaches using machine learning (ML) have been proposed to filter incorrect poses based on the crystal structures. However, most of the methods are limited by the availability of structure data. Here, we propose a new pose classification approach, PECAN2 (Pose Classification with 3D Atomic Network 2), without the need for crystal structures, based on a 3D atomic neural network with Point Cloud Network (PCN). The new approach uses the correlation between docking scores and experimental data to assign labels, instead of relying on the crystal structures. We validate the proposed classifier on multiple datasets including human mu, delta, and kappa opioid receptors and SARS-CoV-2 Mpro. Our results demonstrate that leveraging the correlation between docking scores and experimental data alone enhances molecular docking performance by filtering out false positives and false negatives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Protein-ligand binding affinity prediction using multi-instance learning with docking structures

Recent advances in 3D structure-based deep learning approaches demonstrate improved accuracy in predicting protein-ligand binding affinity in drug discovery. These methods complement physics-based computational modeling such as molecular docking for virtual high-throughput screening. Despite recent advances and improved predictive performance, most methods in this category primarily rely on utilizing co-crystal complex structures and experimentally measured binding affinities as both input and output data for model training. Nevertheless, co-crystal complex structures are not readily available and the inaccurate predicted structures from molecular docking can degrade the accuracy of the machine learning methods. We introduce a novel structure-based inference method utilizing multiple molecular docking poses for each complex entity. Our proposed method employs multi-instance learning with an attention network to predict binding affinity from a collection of docking poses. We validate our method using multiple datasets, including PDBbind and compounds targeting the main protease of SARS-CoV-2. The results demonstrate that our method leveraging docking poses is competitive with other state-of-the-art inference models that depend on co-crystal structures. This method offers binding affinity prediction without requiring co-crystal structures, thereby increasing its applicability to protein targets lacking such data.

97 MATHEMATICS AND COMPUTING↗

PECAN2

PECAN2, Pose Classification with 3D Atomic Network 2, presents an innovative approach to pose classification in computational modeling, particularly in the context of molecular docking. Traditional methods like docking, which rely on physics-based calculations, are prone to inaccuracies in predicting binding poses during experimental testing. While machine learning (ML) approaches, such as ML-driven pose classification, have been introduced to address these issues, they often depend on the availability of crystal structure data. In contrast, PECAN2 introduces a novel pose classification method that operates independently of crystal structures. Instead, it utilizes a 3D atomic neural network with Point Cloud Network (PCN) to establish correlations between docking scores and experimental data. This approach marks a departure from previous studies that heavily relied on crystal structures for labeling. The key innovation lies in the ability of PECAN2 to enhance the performance of molecular docking by filtering out false positives and false negatives through the correlation between docking scores and experimental data.

Shim, Heesung↗

Structure- and Interaction-Based Design of Anti-SARS-CoV-2 Aptamers

Aptamer selection against novel infections is a complicated and time-consuming approach. Synergy can be achieved by using computational methods together with experimental procedures. In this study, we aim to develop a reliable methodology for a rational aptamer in silico et vitro design. The new approach combines multiple steps: (1) Molecular design, based on screening in a DNA aptamer library and directed mutagenesis to fit the protein tertiary structure; (2) 3D molecular modeling of the target; (3) Molecular docking of an aptamer with the protein; (4) Molecular dynamics (MD) simulations of the complexes; (5) Quantum-mechanical (QM) evaluation of the interactions between aptamer and target with further analysis; (6) Experimental verification at each cycle for structure and binding affinity by using small-angle X-ray scattering, cytometry, and fluorescence polarization. By using a new iterative design procedure, structure- and interaction-based drug design (SIBDD), a highly specific aptamer to the receptor-binding domain of the SARS-CoV-2 spike protein, was developed and validated. The SIBDD approach enhances speed of the high-affinity aptamers development from scratch, using a target protein structure. The method could be used to improve existing aptamers for stronger binding. This approach brings to an advanced level the development of novel affinity probes, functional nucleic acids. It offers a blueprint for the straightforward design of targeting molecules for new pathogen agents and emerging variants.

60 APPLIED LIFE SCIENCES↗

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↗

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↗