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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS2001

Structural models of the MscL gating mechanism

Three-dimensional structural models of the mechanosensitive channel of large conductance, MscL, from the bacteria Mycobacterium tuberculosis and Escherichia coli were developed for closed, intermediate, and open conformations. The modeling began with the crystal structure of M. tuberculosis MscL, a homopentamer with two transmembrane alpha-helices, M1 and M2, per subunit. The first 12 N-terminal residues, not resolved in the crystal structure, were modeled as an amphipathic alpha-helix, called S1. A bundle of five parallel S1 helices are postulated to form a cytoplasmic gate. As membrane tension induces expansion, the tilts of M1 and M2 are postulated to increase as they move away from the axis of the pore. Substantial expansion is postulated to occur before the increased stress in the S1 to M1 linkers pulls the S1 bundle apart. During the opening transition, the S1 helices and C-terminus amphipathic alpha-helices, S3, are postulated to dock parallel to the membrane surface on the perimeter of the complex. The proposed gating mechanism reveals critical spatial relationships between the expandable transmembrane barrel formed by M1 and M2, the gate formed by S1 helices, and "strings" that link S1s to M1s. These models are consistent with numerous experimental results and modeling criteria.

NASA Discipline Cell Biology↗
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↗
DOE OSTI2024

Identification and characterization of the low molecular mass ferredoxins involved in central metabolism in Heliomicrobium modesticaldum

The homodimeric Type I reaction center (RC) from Heliomicrobium modesticaldum lacks the PsaC subunit found in Photosystem I and instead uses the interpolypeptide [4Fe–4S] cluster F X as the terminal electron acceptor. Our goal was to identify which of the small mobile dicluster ferredoxins encoded by the H. modesticaldum genome are capable of accepting electrons from the heliobacterial RC (HbRC) and pyruvate:ferredoxin oxidoreductase (PFOR), a key metabolic enzyme. Analysis of the genome revealed seven candidates: HM1_1462 (PshB1), HM1_1461 (PshB2), HM1_2505 (Fdx3), HM1_0869 (FdxB), HM1_1043, HM1_0357, and HM1_2767. Heterologous expression in Escherichia coli and studies using time-resolved optical spectroscopy revealed that only PshB1, PshB2, and Fdx3 are capable of accepting electrons from the HbRC and PFOR. Modeling studies using AlphaFold show that only PshB1, PshB2, and Fdx3 should be capable of docking on PFOR at a positively charged patch that overlays a surface-proximal [4Fe–4S] cluster. Proteomic analysis of wild-type and gene deletion strains ΔpshB1, ΔpshB2, ΔpshB1pshB2, and Δfdx3 grown under nitrogen-replete conditions revealed that Fdx3 is undetectable in the wild-type, ΔpshB1, and Δfdx3 strains, but it is present in the ΔpshB2 and ΔpshB1pshB2 strains, implying that Fdx3 may substitute for PshB2. When grown under nitrogen-deplete conditions, Fdx3 is present in the wild-type and all deletion strains except for Δfdx3. None of the knockout strains demonstrated significant impairment during chemotrophic dark growth on pyruvate, photoheterotrophic light growth on pyruvate, or phototrophic growth on acetate+CO 2 , indicating a high degree of redundancy among these three electron transfer proteins. As a result, loss of both PshB1 and PshB2, but not FdxB, resulted in poor growth under N 2 -fixing conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

DNA-PAINT Imaging with Hydrogel Imprinting and Clearing

Hydrogel-embedding is a versatile technique in fluorescence microscopy, offering stabilization, optical clearing, and the physical expansion of biological specimens. DNA-PAINT is a super-resolution microscopy approach based on the diffusion and transient binding of fluorescently labeled oligos, but its feasibility in hydrogels has not yet been explored. In this study, we demonstrate that polyacrylamide hydrogels support sufficient diffusion for effective DNA-PAINT imaging. Using acrydite-anchored oligonucleotides imprinted from patterned DNA origami nanostructures and microtubule filaments in fixed cells, we find that hydrogel embedding preserves docking strand positioning at the nanoscale. Sample clearing via protease treatment had minor structural effects on the microtubule structure and enhanced diffusion and accessibility to hydrogel-imprinted docking strands. Our work demonstrates promising potential for diffusion and binding-based fluorescence imaging applications in hydrogel-embedded samples.

DNA origami↗
DOE OSTI2024

Structural and Interactional Analysis of the Flavonoid Pathway Proteins: Chalcone Synthase, Chalcone Isomerase and Chalcone Isomerase-like Protein

Chalcone synthase (CHS) and chalcone isomerase (CHI) catalyze the first two committed steps of the flavonoid pathway that plays a pivotal role in the growth and reproduction of land plants, including UV protection, pigmentation, symbiotic nitrogen fixation, and pathogen resistance. Based on the obtained X-ray crystal structures of CHS, CHI, and chalcone isomerase-like protein (CHIL) from the same monocotyledon, Panicum virgatum, along with the results of the steady-state kinetics, spectroscopic/thermodynamic analyses, intermolecular interactions, and their effect on each catalytic step are proposed. In addition, PvCHI’s unique activity for both naringenin chalcone and isoliquiritigenin was analyzed, and the observed hierarchical activity for those type-I and -II substrates was explained with the intrinsic characteristics of the enzyme and two substrates. The structure of PvCHS complexed with naringenin supports uncompetitive inhibition. PvCHS displays intrinsic catalytic promiscuity, evident from the formation of p-coumaroyltriacetic acid lactone (CTAL) in addition to naringenin chalcone. In the presence of PvCHIL, conversion of p-coumaroyl-CoA to naringenin through PvCHS and PvCHI displayed ~400-fold increased Vmax with reduced formation of CTAL by 70%. Supporting this model, molecular docking, ITC (Isothermal Titration Calorimetry), and FRET (Fluorescence Resonance Energy Transfer) indicated that both PvCHI and PvCHIL interact with PvCHS in a non-competitive manner, indicating the plausible allosteric effect of naringenin on CHS. Significantly, the presence of naringenin increased the affinity between PvCHS and PvCHIL, whereas naringenin chalcone decreased the affinity, indicating a plausible feedback mechanism to minimize spontaneous incorrect stereoisomers. These are the first findings from a three-body system from the same species, indicating the importance of the macromolecular assembly of CHS-CHI-CHIL in determining the amount and type of flavonoids produced in plant cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Structural characterization of a GNAT family acetyltransferase from Elizabethkingia anophelis bound to acetyl-CoA reveals a new dimeric interface

General control non-repressible 5 (GCN5)-related N-acetyltransferases (GNATs) catalyse the acetylation of a diverse range of substrates, thereby orchestrating a variety of biological processes within prokaryotes and eukaryotes. GNAT enzymes can catalyze the transfer of an acetyl group from acetyl coenzyme A to substrates such as aminoglycoside antibiotics, amino acids, polyamines, peptides, vitamins, catecholamines, and large macromolecules including proteins. Although GNATs generally exhibit low to moderate sequence identity, they share a conserved catalytic fold and conserved structural motifs. In this current study we characterize the high-resolution X-ray crystallographic structure of a GNAT enzyme bound with acetyl-CoA from Elizabethkingia anophelis , an important multi-drug resistant bacterium. The tertiary structure is comprised of six α-helices and nine β-strands, and is similar with other GNATs. We identify a new and uncharacterized GNAT dimer interface, which is conserved in at least two other unpublished GNAT structures. This suggests that GNAT enzymes can form at least five different types of dimers, in addition to a range of other oligomers including trimer, tetramer, hexamer, and dodecamer assemblies. The high-resolution structure presented in this study is suitable for future in-silico docking and structure–activity relationship studies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Correlated conformational dynamics of the human GluN1-GluN2A type N-methyl-D-aspartate (NMDA) receptor

N-Methyl-D-aspartate receptors (NMDARs) are glutamate-gated ion channels found in the nerve cell membranes. As a result of overexcitation of NMDARs, neuronal death occurs and may lead to diseases such as epilepsy, stroke, Alzheimer’s disease, and Parkinson’s disease. In this study, human GluN1- GluN2A type NMDAR structure is modeled based on the X-ray structure of the Xenopus laevis template and missing loops are added by ab-initio loop modeling. The final structure is chosen according to two different model assessment scores. To be able to observe the structural changes upon ligand binding, glycine and glutamate molecules are docked into the corresponding binding sites of the receptor. Subsequently, molecular dynamics simulations of 1.3 μs are performed for both apo and ligand-bound structures. Structural parameters, which have been considered to show functionally important changes in previous NMDAR studies, are monitored as conformational rulers to understand the dynamics of the conformational changes. Moreover, principal component analysis (PCA) is performed for the equilibrated part of the simulations. Additionally, from these analyses, the differences in between apo and ligand-bound simulations can be summarized as the following: The girdle right at the beginning of the pore loop, which connects M2 and M3 helices of the ion channel, partially opens. Ligands act like an adhesive for the ligand-binding domain (LBD) by keeping the bi-lobed structure together and consequently this is reflected to the overall dynamics of the protein as an increased correlation of the LBD with especially the amino-terminal domain (ATD) of the protein.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Data Science-Driven Analysis of Substrate-Permissive Diketopiperazine Reverse Prenyltransferase NotF: Applications in Protein Engineering and Cascade Biocatalytic Synthesis of (-)-Eurotiumin A

Prenyltransfer is an early-stage carbon-hydrogen bond (C-H) functionalization prevalent in the biosynthesis of a diverse array of biologically active bacterial, fungal, plant, and metazoan diketopiperazine (DKP) alkaloids. Toward the development of a unified strategy for biocatalytic construction of prenylated DKP indole alkaloids, we sought to identify and characterize a substrate-permissive C2 reverse prenyltransferase (PT). As the first tailoring event within the biosynthesis of cytotoxic notoamide metabolites, PT NotF catalyzes C2 reverse prenyltransfer of brevianamide F. Solving a crystal structure of NotF (in complex with native substrate and prenyl donor mimic dimethylallyl S-thiolodiphosphate (DMSPP)) revealed a large, solvent-exposed active site, intimating NotF may possess a significantly broad substrate scope. To assess the substrate selectivity of NotF, we synthesized a panel of 30 sterically and electronically differentiated tryptophanyl DKPs, the majority of which were selectively prenylated by NotF in synthetically useful conversions (2 to > 99%). Quantitative representation of this substrate library and development of a descriptive statistical model provided insight into the molecular origins of NotF's substrate promiscuity. This approach enabled the identification of key substrate descriptors (electrophilicity, size, and flexibility) that govern the rate of NotF-catalyzed prenyltransfer, and the development of an "induced fit docking (IFD)-guided" engineering strategy for improved turnover of our largest substrates. We further demonstrated the utility of NotF in tandem with oxidative cyclization using flavin monooxygenase, BvnB. This one-pot, in vitro biocatalytic cascade enabled the first chemoenzymatic synthesis of the marine fungal natural product, (-)-eurotiumin A, in three steps and 60% overall yield.

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

Structural and biochemical basis for regiospecificity of the flavonoid glycosyltransferase UGT95A1

Glycosylation is a predominant strategy plants use to fine-tune the properties of small molecule metabolites to affect their bioactivity, transport, and storage. It is also important in biotechnology and medicine as many glycosides are utilized in human health. Small molecule glycosylation is largely carried out by family 1 glycosyltransferases. Here, we report a structural and biochemical investigation of UGT95A1, a family 1 GT enzyme from Pilosella officinarum that exhibits a strong, unusual regiospecificity for the 3'-O position of flavonoid acceptor substrate luteolin. We obtained an apo crystal structure to help drive the analyses of a series of binding site mutants, revealing that while most residues are tolerant to mutations, key residues M145 and D464 are important for overall glycosylation activity. Interestingly, E347 is crucial for maintaining the strong preference for 3'-O glycosylation, while R462 can be mutated to increase regioselectivity. The structural determinants of regioselectivity were further confirmed in homologous enzymes. Our study also suggests that the enzyme contains large, highly dynamic, disordered regions. We showed that while most disordered regions of the protein have little to no implication in catalysis, the disordered regions conserved among investigated homologs are important to both the overall efficiency and regiospecificity of the enzyme. This report represents a comprehensive in-depth analysis of a family 1 GT enzyme with a unique substrate regiospecificity and may provide a basis for enzyme functional prediction and engineering.

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

Structural insights on ligand recognition at the human leukotriene B4 receptor 1

The leukotriene B4 receptor 1 (BLT1) regulates the recruitment and chemotaxis of different cell types and plays a role in the pathophysiology of infectious, allergic, metabolic, and tumorigenic human diseases. Here we present a crystal structure of human BLT1 (hBLT1) in complex with a selective antagonist MK-D-046, developed for the treatment of type 2 diabetes and other inflammatory conditions. Comprehensive analysis of the structure and structure-activity relationship data, reinforced by site-directed mutagenesis and docking studies, reveals molecular determinants of ligand binding and selectivity toward different BLT receptor subtypes and across species. The structure helps to identify a putative membrane-buried ligand access channel as well as potential receptor binding modes of endogenous agonists. These structural insights of hBLT1 enrich our understanding of its ligand recognition and open up future avenues in structure-based drug design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

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

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

Jones, William↗
DOE OSTI2025

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

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

Chemical structure↗