Search NASA⌕ Search

SEARCH · Search NASA

Results for “protein 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Structure-based development of a subtype-selective orexin 1 receptor antagonist

Orexin receptors belong to the superfamily of G-protein–coupled receptors (GPCRs) which represent the largest class of drug targets in humans. Despite the recent progress in structural biology, the development of subtype-selective orexin receptor and GPCR ligands in general remains challenging, due to the high sequence similarity among individual receptor subtypes. However, subtype-selective molecules are key to discerning the individual contributions of receptor subtypes to (patho-)physiology. Starting from the clinically used, non–subtype-selective orexin receptor antagonist suvorexant, we demonstrate how docking, crystallography, medicinal chemistry, and in vitro pharmacology can be combined to exploit single amino acid sequence differences for the discovery of subtype-selective probes. Such compounds will help to better understand orexin receptor pharmacology and develop promising drug candidates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Putting AlphaFold models to work with phenix.process_predicted_model and ISOLDE

AlphaFold has recently become an important tool in providing models for experimental structure determination by X-ray crystallography and cryo-EM. Large parts of the predicted models typically approach the accuracy of experimentally determined structures, although there are frequently local errors and errors in the relative orientations of domains. Importantly, residues in the model of a protein predicted by AlphaFold are tagged with a predicted local distance difference test score, informing users about which regions of the structure are predicted with less confidence. AlphaFold also produces a predicted aligned error matrix indicating its confidence in the relative positions of each pair of residues in the predicted model. The phenix.process_predicted_model tool downweights or removes low-confidence residues and can break a model into confidently predicted domains in preparation for molecular replacement or cryo-EM docking. These confidence metrics are further used in ISOLDE to weight torsion and atom–atom distance restraints, allowing the complete AlphaFold model to be interactively rearranged to match the docked fragments and reducing the need for the rebuilding of connecting regions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structure of coagulation factor VIII bound to a patient-derived anti-C1 domain antibody inhibitor

The development of pathogenic antibody inhibitors against coagulation factor VIII (FVIII) occurs in ~30% of congenital hemophilia A patients receiving FVIII replacement therapy as well as in all cases of acquired hemophilia A. KM33 is an anti-C1 domain antibody inhibitor previously isolated from a severe hemophilia A patient. In addition to potently blocking FVIII binding to von Willebrand factor and phospholipid surfaces, KM33 disrupts FVIII binding to lipoprotein receptor-related protein 1 (LRP1), which drives FVIII hepatic clearance and antigen presentation in dendritic cells. Here, we report on the structure of FVIII bound to NB33, a recombinant derivative of KM33, by single-particle cryo-electron microscopy. Furthermore, structural analysis reveals the NB33 epitope localizes to FVIII residues R2090-S2094 and I2158-R2159 which constitute membrane-binding loops in the C1 domain. Further analysis reveals multiple FVIII lysine and arginine residues, previously shown to mediate binding to LRP1, dock onto an acidic cleft at the NB33 variable domain interface, thus blocking a putative LRP1 binding site. Together, these results demonstrate a novel mechanism of FVIII inhibition by a patient-derived antibody inhibitor and provide structural evidence toward engineering FVIII with reduced LRP1-mediated clearance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Viral Nuclease Inhibitors: Small molecule disruptors of the UL12 alkaline nuclease display broad anti-herpes virus activity

Herpes simplex virus 1 (HSV-1) UL12 encodes a highly conserved 5′ → 3′ alkaline exonuclease that is essential for the production of infectious virus. Together with the viral single-stranded DNA-binding/annealing protein ICP8, UL12 functions as a two-component recombinase that mediates recombination-dependent viral DNA replication. Here, we present the crystal structure of the catalytic domain of the HSV alkaline nuclease (UL12), which provides the first view of an α-herpesvirus alkaline nuclease. Using this structure, we optimized a series of small-molecule viral nuclease inhibitors (VNIs) that target the UL12 active site and potently inhibit UL12 exonuclease activity in vitro. We have thus established a robust platform for structure-based docking, SAR analysis and rational inhibitor design. Because UL12 orthologs are conserved across all human herpesviruses, we examined the activity of these compounds against the β- and γ-herpesvirus alkaline nucleases UL98 and SOX and found that they inhibit all three enzymes. The VNIs also exhibit antiviral activity against HSV-1 and HCMV in cell culture. EC 50 and IC 50 values were in the nanomolar to low micromolar range. Together, these findings establish herpesvirus alkaline nucleases as conserved, druggable antiviral targets and provide a foundation for the development of broad-spectrum anti-herpesvirus therapeutics, either as standalone agents or in combination with existing nucleoside analogs.

Sharma, Nidhi↗
DOE OSTI2022

AlphaFold and Structural Mass Spectrometry Enable Interrogations on the Intrinsically Disordered Regions in Cyanobacterial Light-harvesting Complex Phycobilisome

Intrinsically disordered proteins/regions (IDPRs) are a very large and functionally important class of proteins that participate in weak multivalent interactions in protein complexes. They are recalcitrant for interrogations using X-ray crystallography and cryo-EM. The IDPRs observed at the interface of the photosynthetic pigment protein complexes (PPCs) remain much less clear, e.g., the major cyanobacterial light-harvesting complex (PBS) contains an unstructured PB-loop insertion in the phycocyanobilin domain (PB domain) of ApcE (the largest polypeptide in PBS). Here, a joint platform is built to probe such structural domains. This platform is characterized by two-round progressive justifications of in silico models by using the structural mass spectrometry data. First, the AlphaFold-generated 3D structure of the PB domain (containing PB-loop) was justified in the context of PBS. Second, docking the AlphaFold-generated ApcG (a ligand) into the first-step justified structure (a receptor). The final ligand-receptor complex was then subjected to a second-round justification, again, by using unequivocal isotopically-encoded cross-links identified in LC-MS/MS. This work reveals a full-length PB-loop structure modelled in the PBS basal cylinder, free from any spatial conflicts against the other subunits in PBS. The structure of PB domain highlights the close associations of the intrinsically disordered PB-loop with its binding partners in PBS, including ApcG, another IDPR. The PB-loop region involved in the binding of photosystem II (PSII) is also discussed in the context of excitation energy transfer regulation. Finally, this work calls attention to the highly disordered, yet interrogatable interface between the light-harvesting antenna complexes and the reaction centers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Machine Learning-driven Molecular Design for Therapeutic Discovery

The ongoing novel coronavirus pandemic (COVID-19) has highlighted the need for new therapeutics to counter the threat of emerging viral pathogens. The main proteases are a promising target for developing antiviral inhibitors. In this work, we utilized a novel combination of artificial intelligence-driven iterative design of covalent inhibitor candidates, physics-based computational modeling of protein-inhibitor interactions, and “All in One” Native MS biophysical assay screening and characterization of therapeutic candidates. With our existing expertise in hit generation using a particular scaffold as a starting point, we first generated tens of thousands of compounds that preserve the key scaffold. In order to optimize the candidates, we calculated about 136 descriptors consisting of 2D and 3D features for molecules targeting the SARS-CoV-2 Main protease (Mpro). These compounds were initially filtered according to properties and further sorted by predicted binding affinity using our automated docking modeling and machine learning methods. We tested a handful of candidates and identified two as inhibitors of Mpro with micromolar affinities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

The N‐terminal intrinsically disordered region of Ncb5or docks with the cytochrome b 5 core to form a helical motif that is of ancient origin

Abstract NADH cytochromeb 5 oxidoreductase (Ncb5or) is a cytosolic ferric reductase implicated in diabetes and neurological conditions. Ncb5or comprises cytochromeb 5 (b 5 ) and cytochromeb 5 reductase (b 5 R) domains separated by a CHORD‐Sgt1 (CS) linker domain. Ncb5or redox activity depends on proper inter‐domain interactions to mediate electron transfer from NADH or NADPH via FAD to heme. While full‐length human Ncb5or has proven resistant to crystallization, we have succeeded in obtaining high‐resolution atomic structures of theb 5 domain and a construct containing the CS andb 5 R domains (CS/b 5 R). Ncb5or also contains an N‐terminal intrinsically disordered region of 50 residues that has no homologs in other protein families in animals but features a distinctive, conserved L 34 MDWIRL 40 motif also present in reduced lateral root formation (RLF) protein in rice and increased recombination center 21 in baker's yeast, all attaching to ab 5 domain. After unsuccessful attempts at crystallizing a human Ncb5or construct comprising the N‐terminal region naturally fused to theb 5 domain, we were able to obtain a high‐resolution atomic structure of a recombinant rice RLF construct corresponding to residues 25–129 of human Ncb5or (52% sequence identity; 74% similarity). The structure reveals Trp 120 (corresponding to invariant Trp 37 in Ncb5or) to be part of an 11‐residue α‐helix (S 116 QMDWLKLTRT 126 ) packing against two of the four helices in theb 5 domain that surround heme (α2 and α5). The Trp 120 side chain forms a network of interactions with the side chains of four highly conserved residues corresponding to Tyr 85 and Tyr 88 (α2), Cys 124 (α5), and Leu 47 in Ncb5or. Circular dichroism measurements of human Ncb5or fragments further support a key role of Trp 37 in nucleating the formation of the N‐terminal helix, whose location in the N/b 5 module suggests a role in regulating the function of this multi‐domain redox enzyme. This study revealed for the first time an ancient origin of a helical motif in the N/b 5 module as reflected by its existence in a class of cytochromeb 5 proteins from three kingdoms among eukaryotes.

Biochemistry & Molecular Biology↗
DOE OSTI2020

Molecular basis of P[II] major human rotavirus VP8* domain recognition of histo-blood group antigens

Initial cell attachment of rotavirus (RV) to specific cell surface glycan receptors, which is the essential first step in RV infection, is mediated by the VP8* domain of the spike protein VP4. Recently, human histo-blood group antigens (HBGAs) have been identified as receptors or attachment factors for human RV strains. RV strains in the P[4] and P[8] genotypes of the P[II] genogroup share common recognition of the Lewis b (Le b ) and H type 1 antigens, however, the molecular basis of receptor recognition by the major human P[8] RVs remains unknown due to lack of experimental structural information. Here, we used nuclear magnetic resonance (NMR) spectroscopy-based titration experiments and NMR-derived high ambiguity driven docking (HADDOCK) methods to elucidate the molecular basis for P[8] VP8* recognition of the Leb (LNDFH I) and type 1 HBGAs. We also used X-ray crystallography to determine the molecular details underlying P[6] recognition of H type 1 HBGAs. Unlike P[6]/P[19] VP8*s that recognize H type 1 HBGAs in a binding surface composed of an α-helix and a β-sheet, referred as the “βα binding site”, the P[8] and P[4] VP8*s bind Le b HBGAs in a previously undescribed pocket formed by the edges of two β-sheets, referred to as the “ββ binding site”. Importantly, the P[8] and P[4] VP8*s retain binding capability to non-Le b type 1 HBGAs using the βα binding site. The presence of two distinct binding sites for Le b and non-Le b HBGA glycans in the P[8] and P[4] VP8* domains suggests host-pathogen co-evolution under structural and functional adaptation of RV pathogens to host glycan polymorphisms. Assessment and understanding of the precise impact of this co-evolutionary process in determining RV host ranges and cross-species RV transmission should facilitate improved RV vaccine development and prediction of future RV strain emergence and epidemics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Structural basis for the specificity of PPM1H phosphatase for Rab GTPases

LRRK2 serine/threonine kinase is associated with inherited Parkinson’s disease. LRRK2 phosphorylates a subset of Rab GTPases within their switch 2 motif to control their interactions with effectors. Recent work has shown that the metal-dependent protein phosphatase PPM1H counteracts LRRK2 by dephosphorylating Rabs. PPM1H is highly selective for LRRK2 phosphorylated Rabs, and closely related PPM1J exhibits no activity towards substrates such as Rab8a phosphorylated at Thr72 (pThr72). Here, we have identified the molecular determinant of PPM1H specificity for Rabs. The crystal structure of PPM1H reveals a structurally conserved phosphatase fold that strikingly has evolved a 110-residue flap domain adjacent to the active site. The flap domain distantly resembles tudor domains that interact with histones in the context of epigenetics. Cellular assays, crosslinking and 3-D modelling suggest that the flap domain encodes the docking motif for phosphorylated Rabs. Consistent with this hypothesis, a PPM1J chimaera with the PPM1H flap domain dephosphorylates pThr72 of Rab8a both in vitro and in cellular assays. Therefore, PPM1H has acquired a Rab-specific interaction domain within a conserved phosphatase fold.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Mechanistic Duality of Bacterial Efflux Substrates and Inhibitors: Example of Simple Substituted Cinnamoyl and Naphthyl Amides

Antibiotic resistance poses an immediate and growing threat to human health. Multidrug efflux pumps are promising targets for overcoming antibiotic resistance with small-molecule therapeutics. Previously, we identified a diaminoquinoline acrylamide, NSC-33353, as a potent inhibitor of the AcrAB–TolC efflux pump in Escherichia coli. This inhibitor potentiates the antibacterial activities of novobiocin and erythromycin upon binding to the membrane fusion protein AcrA. It is also a substrate for efflux and lacks appreciable intrinsic antibacterial activity of its own in wild-type cells. Here, we have modified the substituents of the cinnamoyl group of NSC-33353, giving rise to analogs that retain the ability to inhibit efflux, lost the features of the efflux substrates, and gained antibacterial activity in wild-type cells. The replacement of the cinnamoyl group with naphthyl isosteres generated compounds that lack antibacterial activity but are both excellent efflux pump inhibitors and substrates. Surprisingly, these inhibitors potentiate the antibacterial activity of novobiocin but not erythromycin. Surface plasmon resonance experiments and molecular docking suggest that the replacement of the cinnamoyl group with naphthyl shifts the affinity of the compounds away from AcrA to the AcrB transporter, making them better efflux substrates and changing their mechanism of inhibition. These results provide new insights into the duality of efflux substrate/inhibitor features in chemical scaffolds that will facilitate the development of new efflux pump inhibitors.

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

Antihypertensive effect of quinoa protein under simulated gastrointestinal digestion and peptide characterization

Abstract BACKGROUND Quinoa protein is a potential source of bioactive peptides. Although some studies have demonstrated its angiotensin converting enzyme (ACE) inhibitory properties, research into its in vivo effect on blood‐pressure regulation and peptide characterization remains limited. RESULTS Quinoa protein hydrolyzate (QPH) was prepared by simulated gastrointestinal digestion. QPH lowered the systolic blood pressure (SBP) and diastolic blood pressure (DBP) in spontaneously hypertensive model rats (SHRs) from 2 h to10 h after oral administration, effectively controlling blood pressure in these SHRs. An in vitro study showed that QPH is capable of inhibiting ACE activity. This was attributed to the activity of a number of low‐molecular‐weight peptides. With relatively high scores predicted by PeptideRanker, three promising bioactive peptides, FHPFPR, NWFPLPR, and NIFRPF, were further studied and their ACE‐inhibition effects were confirmed with IC 50 values of 34.92, 16.77, and 32.40 μM, respectively. A molecular docking study provided insights into the binding of ACE with peptides, and revealed that the presence of specific amino acids in the peptide sequence (Pro, Phe, and Arg at the C‐terminal, and Asn at the N‐terminal) could contribute to the interaction between ACE and peptides. CONCLUSION These results demonstrated the potential of QPH for the management of hypertension, which indicates that it could be a good candidate for inclusion in functional foods to control high blood pressure. © 2020 Society of Chemical Industry

Guo, Huimin↗
DOE OSTI2020

Ligand-Binding-Site Refinement to Generate Reliable Holo Protein Structure Conformations from Apo Structures

The first important step in a structure-based virtual screening is the judicious selection of a receptor protein. In cases where the holo protein receptor structure is unavailable, significant reduction in virtual screening performance has been reported. In this work, we present a robust method to generate reliable holo protein structure conformations from apo structures using molecular dynamics (MD) simulation with restraints derived from holo structure binding-site templates. We perform benchmark tests on two different datasets: 40 structures from a directory of useful decoy-enhanced (DUD-E) and 84 structures from the Gunasekaran dataset. Our results show successful refinement of apo binding-site structures toward holo conformations in 82% of the test cases. In addition, virtual screening performance of 40 DUD-E structures is significantly improved using our MD-refined structures as receptors with an average enrichment factor (EF), an EF1% value of 6.2 compared to apo structures with 3.5. Docking of native ligands to the refined structures shows an average ligand root mean square deviation (RMSD) of 1.97 Å (DUD-E dataset and Gunasekaran dataset) relative to ligands in the holo crystal structures, which is comparable to the self-docking (i.e., docking of the native ligand back to its crystal structure receptor) average, 1.34 Å (DUD-E dataset) and 1.36 Å (Gunasekaran dataset). On the other hand, docking to the apo structures yields an average ligand RMSD of 3.65 Å (DUD-E) and 2.90 Å (Gunasekaran). These results indicate that our method is robust and can be useful to improve virtual screening performance of apo structures.

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

NUB1 traps unfolded FAT10 for ubiquitin-independent degradation by the 26S proteasome

The ubiquitin-like modifier FAT10 targets hundreds of proteins in the mammalian immune system to the 26S proteasome for degradation. This degradation pathway requires the cofactor NUB1, yet the underlying mechanisms remain unknown. Here, we reconstituted a minimal in vitro system with human components and revealed that NUB1 uses the intrinsic instability of FAT10 to trap its N-terminal ubiquitin-like domain in an unfolded state and deliver it to the 26S proteasome for engagement, allowing the degradation of FAT10-ylated substrates in a ubiquitin-independent and p97-independent manner. Using hydrogen–deuterium exchange, structural modeling and site-directed mutagenesis, we identified the formation of an intricate complex with FAT10 that activates NUB1 for docking to the 26S proteasome, and our cryo-EM studies visualized the highly dynamic NUB1 complex bound to the proteasomal Rpn1 subunit during FAT10 delivery and the early stages of ATP-dependent degradation. These findings identified a previously unknown mode of cofactor-mediated, ubiquitin-independent substrate delivery to the 26S proteasome that relies on trapping partially unfolded states for engagement by the proteasomal ATPase motor.

Arkinson, Connor↗
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 OSTI2023

Crystal structure of the CoV-Y domain of SARS-CoV-2 nonstructural protein 3

Abstract Replication of the coronavirus genome starts with the formation of viral RNA-containing double-membrane vesicles (DMV) following viral entry into the host cell. The multi-domain nonstructural protein 3 (nsp3) is the largest protein encoded by the known coronavirus genome and serves as a central component of the viral replication and transcription machinery. Previous studies demonstrated that the highly-conserved C-terminal region of nsp3 is essential for subcellular membrane rearrangement, yet the underlying mechanisms remain elusive. Here we report the crystal structure of the CoV-Y domain, the most C-terminal domain of the SARS-CoV-2 nsp3, at 2.4 Å-resolution. CoV-Y adopts a previously uncharacterized V-shaped fold featuring three distinct subdomains. Sequence alignment and structure prediction suggest that this fold is likely shared by the CoV-Y domains from closely related nsp3 homologs. NMR-based fragment screening combined with molecular docking identifies surface cavities in CoV-Y for interaction with potential ligands and other nsps. These studies provide the first structural view on a complete nsp3 CoV-Y domain, and the molecular framework for understanding the architecture, assembly and function of the nsp3 C-terminal domains in coronavirus replication. Our work illuminates nsp3 as a potential target for therapeutic interventions to aid in the on-going battle against the COVID-19 pandemic and diseases caused by other coronaviruses.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Insight into molecular basis and dynamics of full-length CRaf kinase in cellular signaling mechanisms

Raf kinases play key roles in signal transduction in cells for regulating proliferation, differentiation, and survival. Despite decades of research into functions and dynamics of Raf kinases with respect to other cytosolic proteins, understanding Raf kinases is limited by the lack of their full-length structures at the atomic resolution. Here, we present the first model of the full-length CRaf kinase obtained from artificial intelligence/machine learning algorithms with a converging ensemble of structures simulated by large-scale temperature replica exchange simulations. Our model is validated by comparing simulated structures with the latest cryo-EM structure detailing close contacts among three key domains and regions of the CRaf. Our simulations identify potentially new epitopes of intramolecule interactions within the CRaf and reveal a dynamical nature of CRaf kinases, in which the three domains can move back and forth relative to each other for regulatory dynamics. The dynamic conformations are then used in a docking algorithm to shed insight into the paradoxical effect caused by vemurafenib in comparison with a paradox breaker PLX7904. In this study, we propose a model of Raf-heterodimer/KRas-dimer as a signalosome based on the dynamics of the full-length CRaf.

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