Search NASA⌕ Search

SEARCH · Search NASA

Results for “docking simulations”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

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↗
DOE OSTI2022

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↗
DOE OSTI2024

Novel sterol binding domains in bacteria

Sterol lipids are widely present in eukaryotes and play essential roles in signaling and modulating membrane fluidity. Although rare, some bacteria also produce sterols, but their function in bacteria is not known. Moreover, many more species, including pathogens and commensal microbes, acquire or modify sterols from eukaryotic hosts through poorly understood molecular mechanisms. The aerobic methanotroph Methylococcus capsulatus was the first bacterium shown to synthesize sterols, producing a mixture of C-4 methylated sterols that are distinct from those observed in eukaryotes. C-4 methylated sterols are synthesized in the cytosol and localized to the outer membrane, suggesting that a bacterial sterol transport machinery exists. Until now, the identity of such machinery remained a mystery. In this study, we identified three novel proteins that may be the first examples of transporters for bacterial sterol lipids. The proteins, which all belong to well-studied families of bacterial metabolite transporters, are predicted to reside in the inner membrane, periplasm, and outer membrane of M. capsulatus, and may work as a conduit to move modified sterols to the outer membrane. Quantitative analysis of ligand binding revealed their remarkable specificity for 4-methylsterols, and crystallographic structures coupled with docking and molecular dynamics simulations revealed the structural bases for substrate binding by two of the putative transporters. Their striking structural divergence from eukaryotic sterol transporters signals that they form a distinct sterol transport system within the bacterial domain. Finally, bioinformatics revealed the widespread presence of similar transporters in bacterial genomes, including in some pathogens that use host sterol lipids to construct their cell envelopes. The unique folds of these bacterial sterol binding proteins should now guide the discovery of other proteins that handle this essential metabolite.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

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↗
DOE OSTI2022

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↗
DOE OSTI2025

Unraveling the Hsp70-ROS-autophagy axis in pentachlorophenol-challenged lung and liver epithelial cells

Pentachlorophenol (PCP) was extensively utilized as an organochlorine pesticide and wood preservative in the United States from the 1930s until the Environmental Protection Agency (EPA) imposed restrictions due to concerns about its toxicity and potential carcinogenic properties. Although it is no longer widely used, PCP remains a concern due to its environmental persistence and potential for long-term health effects. Significant occupational and environmental exposures have likely occurred, with the health and economic costs of PCP exposure potentially being substantial given its known toxicity. Notably, PCP exhibits rapid absorption through both the skin and respiratory system and has been shown to cause hepatotoxicity, developmental toxicity, immunotoxicity, irritation, and carcinogenicity in laboratory animal studies. PCP exposure induces oxidative stress, a key mechanism underlying its inflammatory and toxic effects, which can activate cellular stress responses including upregulation of heat shock protein 70 (Hsp70). Previous studies in lung and liver epithelial cells have shown that Hsp70 and oxidative stress play pivotal roles in triggering autophagy. This study establishes the critical role of the Hsp70-reactive oxygen species (ROS)-autophagy axis in regulating cellular responses to PCP exposure in human alveolar (A549) and liver carcinoma (HepG2) epithelial cells. Our research elucidated the molecular mechanisms underlying PCP's cellular effects, demonstrating that its exposure resulted in increased expression of autophagy-related proteins (Beclin-1, LC3B, ATG12, and ATG16), subunits of NADPH oxidase (NCF-1, NCF-2, NOX2, and Rac), and antioxidant proteins (SOD and GPx) in both lung and liver cell types. Notably, PCP augmented the interaction between Hsp70 and the autophagy regulator Beclin-1. Pretreatment with the ROS inhibitor N-acetylcysteine or Hsp70 knockdown markedly reversed PCP-induced responses. Our in-silico protein–protein docking analysis and molecular dynamics simulation studies revealed enhanced interactions and/or stable confirmations maintained throughout the simulations for TLR4-Hsp70 and Hsp70-Beclin-1 complexes in the presence of PCP. These findings provide a strong foundation for future studies, employing in vivo experimental models and human populations to identify promising targets for PCP-induced toxicity and cellular injury. As a result, these findings may have far-reaching implications for public health and environmental policy, ultimately leading to the identification of biomarkers and the development of more effective interventions for environmentally induced toxicity and diseases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Screening of Therapeutic Agents for COVID-19 Using Machine Learning and Ensemble Docking Studies

The current pandemic demands a search for therapeutic agents against the novel coronavirus SARS-CoV-2. In this work, we present an efficient computational strategy that combines machine learning (ML)-based models and high-fidelity ensemble docking studies to enable rapid screening of possible therapeutic ligands. Targeting the binding affinity of molecules for either the isolated SARS-CoV-2 S-protein at its host receptor region or the S-protein:human ACE2 interface complex, we screen ligands from drug and biomolecule data sets that can potentially limit and/or disrupt the host-virus interactions. Top scoring one hundred eighty-seven ligands (with 75 approved by the Food and Drug Administration) are further validated by all atom docking studies. Important molecular descriptors ( 2 χ n , topological surface area, and ring count) and promising chemical fragments (oxolane, hydroxy, and imidazole) are identified to guide future experiments. Overall, this work expands our knowledge of small-molecule treatment against COVID-19 and provides a general screening pathway (combining quick ML models with expensive high-fidelity simulations) for targeting several chemical/biochemical problems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Identification of structural transitions in bacterial fatty acid binding proteins that permit ligand entry and exit at membranes

Fatty acid (FA) transfer proteins extract FA from membranes and sequester them to facilitate their movement through the cytosol. Detailed structural information is available for these soluble protein–FA complexes, but the structure of the protein conformation responsible for FA exchange at the membrane is unknown. Staphylococcus aureus FakB1 is a prototypical bacterial FA transfer protein that binds palmitate within a narrow, buried tunnel. Here, we define the conformational change from a “closed” FakB1 state to an “open” state that associates with the membrane and provides a path for entry and egress of the FA. Using NMR spectroscopy, we identified a conformationally flexible dynamic region in FakB1, and X-ray crystallography of FakB1 mutants captured the conformation of the open state. In addition, molecular dynamics simulations show that the new amphipathic α-helix formed in the open state inserts below the phosphate plane of the bilayer to create a diffusion channel for the hydrophobic FA tail to access the hydrocarbon core and place the carboxyl group at the phosphate layer. The membrane binding and catalytic properties of site-directed mutants were consistent with the proposed membrane docked structure predicted by our molecular dynamics simulations. Finally, the structure of the bilayer-associated conformation of FakB1 has local similarities with mammalian FA binding proteins and provides a conceptual framework for how these proteins interact with the membrane to create a diffusion channel from the FA location in the bilayer to the protein interior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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

Identification of inhibitors against SARS-CoV-2 variants of concern using virtual screening and metadynamics-based enhanced sampling

Among the variants of SARS-CoV-2, some are more infectious than the Wild-type. Interestingly, these mutations enable the virus to evade the therapeutic efforts. Hence, there is a need for candidate drug molecules that can potently bind with all the variants. Here we have adopted a strategy combining virtual screening, molecular docking followed by rigorous sampling by metadynamics simulations to find candidate molecules. From our results we found four highly potent drug candidates that can bind to the Spike-RBD of all the variants of the virus. Additionally, we also found that certain signature residues on the RBM region commonly bind to each of these inhibitors. Thus, our study not only gives information on the chemical compounds, but also residues on the proteins which could be targeted for future drug and vaccine development studies.

60 APPLIED LIFE 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 OSTI2021

Functional Characterization and Structural Insights Into Stereoselectivity of Pulegone Reductase in Menthol Biosynthesis

Monoterpenoids are the main components of plant essential oils and the active components of some traditional Chinese medicinal herbs like Mentha haplocalyx Briq ., Nepeta tenuifolia Briq ., Perilla frutescens (L.) Britt and Pogostemin cablin (Blanco) Benth. Pulegone reductase is the key enzyme in the biosynthesis of menthol and is required for the stereoselective reduction of the Δ 2,8 double bond of pulegone to produce the major intermediate menthone, thus determining the stereochemistry of menthol. However, the structural basis and mechanism underlying the stereoselectivity of pulegone reductase remain poorly understood. In this study, we characterized a novel (-)-pulegone reductase from Nepeta tenuifolia ( Nt PR), which can catalyze (-)-pulegone to (+)-menthone and (-)-isomenthone through our RNA-seq, bioinformatic analysis in combination with in vitro enzyme activity assay, and determined the structure of (+)-pulegone reductase from ( Mp PR) by using X-ray crystallography, molecular modeling and docking, site-directed mutagenesis, molecular dynamics simulations, and biochemical analysis. We identified and validated the critical residues in the crystal structure of Mp PR involved in the binding of the substrate pulegone. We also further identified that residues Leu56, Val282, and Val284 determine the stereoselectivity of the substrate pulegone, and mainly contributes to the product stereoselectivity. This work not only provides a starting point for the understanding of stereoselectivity of pulegone reductases, but also offers a basis for the engineering of menthone/menthol biosynthetic enzymes to achieve high-titer, industrial-scale production of enantiomerically pure products.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Molecular Dynamics Simulations of a Catalytic Multivalent Peptide–Nanoparticle Complex

Molecular modeling of a supramolecular catalytic system is conducted resulting from the assembling between a small peptide and the surface of cationic self-assembled monolayers on gold nanoparticles, through a multiscale iterative approach including atomistic force field development, flexible docking with Brownian Dynamics and µs-long Molecular Dynamics simulations. Self-assembly is a prerequisite for the catalysis, since the catalytic peptides do not display any activity in the absence of the gold nanocluster. Atomistic simulations reveal details of the association dynamics as regulated by defined conformational changes of the peptide due to peptide length and sequence. Our results show the importance of a rational design of the peptide to enhance the catalytic activity of peptide–nanoparticle conjugates and present a viable computational approach toward the design of enzyme mimics having a complex structure–function relationship, for technological and nanomedical applications.

59 BASIC BIOLOGICAL 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

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 OSTI2025

Structure-Based Design of Small-Molecule Inhibitors of Human Interleukin-6

Human Interleukin-6 (hIL-6) is a pro inflammatory cytokine that binds to its receptor, IL-6Rα followed by binding to gp130 and subsequent dimerization to form a hexamer signaling complex. As a critical inflammation mediator, hIL-6 is associated with a diverse range of diseases and monoclonal antibodies in clinical use that either target IL-6Rα or hIL-6 to inhibit signaling. Here, we perform high-throughput structure-based computational screening using ensemble docking for small-molecule antagonists for which the target conformations were taken from 600 ns long molecular dynamics simulations of the apo protein. Prior knowledge of the contact sites from binary complex studies and experimental work was incorporated into the docking studies. The top 20 scoring ligands from the in silico studies after post analysis were subjected to in vitro functional assays. Among these compounds, the ligand with the second-highest calculated binding affinity experimentally showed an ~84% inhibitory effect on IL6-induced STAT3 reporter activity at 10 μM concentration. This finding may pave the way for designing small-molecule inhibitors of hIL-6 of therapeutic significance.

Human Interleukin-6↗
DOE OSTI2022

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↗
DOE OSTI2022

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↗