SEARCH · Engineering Papers

Results for “protein function prediction”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Predicting Mössbauer Parameters of Nonheme Diiron Complexes with Density Functional Theory

Mössbauer spectroscopy provides significant insights into the electronic structure and environment of the metal centers. In this work, we investigate the electronic structures of a set of nonheme diiron complexes by evaluating two key parameters pertaining to Mössbauer spectroscopy, namely, the isomer shift (δ) and quadrupole splitting (|ΔE Q |), using different levels of density functional theory (DFT). The diiron systems investigated here span diverse oxidation states, bridging motifs, and spin coupling patterns, which present a challenging case for theoretical predictions. We demonstrate that the combination of B97-D3/def2-TZVP is an efficient approach in modeling both the δ and |ΔEQ| values with high accuracy for the representative nonheme diiron complexes. We also show that δ is accurately predicted irrespective of the choice of approximate density functional while the |ΔE Q | is sensitive to the level of theory employed. Further investigation shows that the present methodology assessed using synthetic nonheme diiron complexes could be extended to nonheme diiron enzyme active sites, featuring both ferromagnetic and antiferromagnetic coupling between the iron centers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Structure and function of the divalent anion/Na + symporter from Vibrio cholerae and a humanized variant

Integral membrane proteins of the divalent anion/Na + symporter (DASS) family translocate dicarboxylate, tricarboxylate or sulphate across cell membranes, typically by utilizing the preexisting Na + gradient. The molecular determinants for substrate recognition by DASS remain obscure, largely owing to the absence of any substrate-bound DASS structure. Here we present 2.8-Å resolution X-ray structures of VcINDY, a DASS from Vibrio cholerae that catalyses the co-transport of Na + and succinate. These structures portray the Na + -bound VcINDY in complexes with succinate and citrate, elucidating the binding sites for substrate and two Na + ions. Furthermore, we report the structures of a humanized variant of VcINDY in complexes with succinate and citrate, which predict how a human citrate-transporting DASS may interact with its bound substrate. Our findings provide insights into metabolite transport by DASS, establishing a molecular basis for future studies on the regulation of this transport process.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2019

X-ray Scattering and Coarse-Grained Simulations for Clustering and Interactions of Monoclonal Antibodies at High Concentrations

Attractive protein-protein interactions (PPI) in concentrated monoclonal antibody (mAb) solutions may lead to reversible oligomers (clusters) that impact colloidal stability and viscosity. Herein, the PPI are tuned for two mAbs via the addition of arginine (Arg), NaCl, or ZnSO 4 as characterized by the structure factor ( S eff (q)) with small angle x-ray scattering (SAXS). The SAXS data are fit with molecular dynamics simulations by placing a physically relevant short-ranged attractive interaction on selected beads in coarse-grained 12 bead models of the mAb shape. The optimized 12-bead models are then used to differentiate key microstructural properties, including center of mass radial distribution functions ( g COM (r)), coordination numbers and cluster size distributions (CSD). The addition of co-solutes results in more attractive Seff(q) relative to the no co-solute control for all systems tested, with the most attractive systems showing an upturn at low q . Only the All1 model with an attractive site in each Fab and Fc region (possessing Fab-Fab, Fab-Fc and Fc-Fc interactions) can reproduce this upturn, and the corresponding CSDs show the presence of larger clusters compared to the control. In general, for models with similar net attractions, i.e., second osmotic virial coefficients, the size of the clusters increases as the attraction is concentrated on a smaller number of evenly distributed beads. Here, the cluster size distributions from simulations are used to improve the understanding and prediction of experimental viscosities. The ability to discriminate between models with bead interactions at particular Fab and Fc bead sites from SAXS simulations, and to provide real space properties (CSD and g COM (r)), will be of interest in engineering protein sequence and formulating protein solutions for weak PPI to minimize aggregation and viscosities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Uncovering Sequence and Structural Characteristics of Fungal Expansin‐Related Proteins With Potential to Drive Substrate Targeting

Expansins loosen plant cell wall networks through disrupting non-covalent bonds between cellulose microfibrils and matrix polysaccharides. Whereas expansins were first discovered in plants, expansin-related proteins have since been identified in bacteria and fungi. The biological function of microbial expansins remains unclear; however, several studies have shown distinct binding preferences toward different structural polysaccharides. Earlier studies of bacterial expansin-related proteins uncovered sequence and structural features that correlate to substrate binding. Herein, 20 fungal expansin-related sequences were recombinantly produced in Komagataella phaffii, and the purified proteins were compared in terms of substrate binding to cellulosic and chitinous substrates. The impact of pH on the zeta potential of prioritized substrates was also measured, and Principal Component Analysis was performed to uncover correlations between protein characteristics (e.g., pI, hydrophobicity, surface charge distribution) and measured substrate binding preferences. Whereas acidic proteins with a predicted pI less than 5.0 preferentially bound to chitin, basic proteins with pI greater than 8.0 preferentially bound to xylan and xylan-containing fiber. Similar to many cellulases, binding to cellulose was correlated to relatively high aromatic amino acid content in the protein sequence and presence of a carbohydrate binding module (CBM), which in the case of expansins is a C-terminal CBM63. Whereas overall sequence characteristics could be correlated to substrate binding preference, the identity of amino acids occupying conserved positions that impact protein activity was better correlated with loosenin versus expansin classifications.

chitin↗

DOE OSTI2021

tSFM 1.0: tRNA Structure–Function Mapper

Structure-conditioned information statistics have proven useful to predict and visualize tRNA Class-Informative Features (CIFs) and their evolutionary divergences. Although permutation P-values can quantify the significance of CIF divergences between two taxa, their naive Monte Carlo approximation is slow and inaccurate. The Peaks-over-Threshold approach of Knijnenburg et al. (2009) promises improvements to both speed and accuracy of permutation P-values, but has no publicly available API. Here, we present tRNA Structure–Function Mapper (tSFM) v1.0, an open-source, multi-threaded application that efficiently computes, visualizes and assesses significance of single- and paired-site CIFs and their evolutionary divergences for any RNA, protein, gene or genomic element sequence family. Multiple estimators of permutation P-values for CIF evolutionary divergences are provided along with confidence intervals. tSFM is implemented in Python 3 with compiled C extensions and is freely available through GitHub (https://github.com/tlawrence3/tSFM) and PyPI.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2016

Hepatic cytochrome P450 activity, abundance, and expression throughout human development

Cytochrome P450s are Phase I metabolic enzymes that play critical roles in the biotransformation of endogenous compounds and xenobiotics. The expression and activity of P450 enzymes can vary considerably throughout human development, especially when comparing fetal development to neonates, children, and adults. In an effort to develop a more comprehensive understanding of the ontogeny of P450 expression and activity we employed a multi-omic characterization of P450 transcript expression, protein abundance, and functional activity. To quantify the functional activity of individual P450s we employ activity-based protein profiling, which uses modified mechanism-based inhibitors of P450s as chemical probes, in tandem with proteomic analyses to quantify activity. Our results reveal life-stage-dependent variability in P450 expression, abundance, and activity throughout human development and frequent discordant relationships between expression and activity. The results were used to distribute P450s into three general classes based upon developmental stage of expression and activity. We have significantly expanded the knowledge of P450 ontogeny, particularly at the level of individual P450 activity. We anticipate that our ontogeny results will be useful for enabling predictive therapeutic dosing, and for avoiding potentially adverse and harmful reactions during maturation from both therapeutic drugs and environmental xenobiotics.

Sadler, Natalie C.↗

DOE OSTI2019

Snails In Silico: A Review of Computational Studies on the Conopeptides

Marine cone snails are carnivorous gastropods that use peptide toxins called conopeptides both as a defense mechanism and as a means to immobilize and kill their prey. These peptide toxins exhibit a large chemical diversity that enables exquisite specificity and potency for target receptor proteins. This diversity arises in terms of variations both in amino acid sequence and length, and in posttranslational modifications, particularly the formation of multiple disulfide linkages. Most of the functionally characterized conopeptides target ion channels of animal nervous systems, which has led to research on their therapeutic applications. Many facets of the underlying molecular mechanisms responsible for the specificity and virulence of conopeptides, however, remain poorly understood. In this review, we will explore the chemical diversity of conopeptides from a computational perspective. First, we discuss current approaches used for classifying conopeptides. Next, we review different computational strategies that have been applied to understanding and predicting their structure and function, from machine learning techniques for predictive classification to docking studies and molecular dynamics simulations for molecular-level understanding. We then review recent novel computational approaches for rapid high-throughput screening and chemical design of conopeptides for particular applications. We close with an assessment of the state of the field, emphasizing important questions for future lines of inquiry.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Analysis of the Active Site Cysteine Residue of the Sacrificial Sulfur Insertase LarE from Lactobacillus plantarum

LarE from Lactobacillus plantarum is an ATP-dependent sulfur transferase that sacrifices its Cys176 sulfur atom to form a dehydroalanine (Dha) side chain during biosynthesis of the covalently linked nickel-pincer nucleotide (NPN) cofactor (pyridinium 3-thioamide-5-thiocarboxylic acid mononucleotide) of lactate racemase. Coenzyme A (CoA) stabilizes LarE and forms a CoA-Cys176 mixed disulfide with the protein. This study presents the crystal structure of the LarE/CoA complex, revealing protein interactions with CoA that mimic those for binding ATP. CoA weakly inhibits LarE activity, and the persulfide of CoA is capable of partially regenerating functional LarE from the Dha176 form of the protein. The physiological relevance of this cycling reaction is unclear. A new form of LarE was discovered, an NPN-LarE covalent adduct, explaining prior results in which activation of the lactate racemase apoprotein required only the isolated LarE. The crystal structure of the inactive C176A variant revealed a fold essentially identical to that of wild-type LarE. Additional active site variants of LarE were created and their activities characterized, with all LarE variants analyzed in terms of the structure. Finally, the L. plantarum LarE structure was compared to a homology model of Thermoanaerobacterium thermosaccharolyticum LarE, predicted to contain three cysteine residues at the active site, and to other proteins with a similar fold and multiple active site cysteine residues. Lastly, these findings suggest that some LarE orthologs may not be sacrificial but instead may catalyze sulfur transfer by using a persulfide mechanism or from a labile site on a [4Fe-4S] cluster at this position.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Boost of the Bio-memristor Performance for Artificial Electronic Synapses by Surface Reconstruction

Biomaterial-based memristors (bio-memristors) are often adopted to emulate biological synapse functions and applied to construct neural computing networks in brain-inspired chip systems. However, the randomness of conductive filament formation in bio-memristors inhibits their switching performance by causing the dispersion of the device-switching parameters. In this case, a facile porous silk fibroin (p-SF) memristor was obtained through a protein surface reconstruction strategy, in which the size of the hole can be adjusted by the density of hybrid nanoseeds. Here, the porous SF memristors exhibit greatly enhanced electrical characteristics, including uniform I-V cycles, centralized distribution of the switching voltages, and both high and low resistances, compared to devices without pores. The results of three-dimensional (3D) simulations based on classical density functional theory (cDFT) suggest that the reconstructed pores in the SF layers guide the formation and fracture of Ag filaments under an electric field and enhance the overall conductivity by separating Ag + ion and electron diffusion pathways. Ag + ions are predicted to preferentially diffuse through pores, whereas electrons diffuse through the SF network. Interestingly, the device conductance can be bidirectionally modulated gradually by positive and negative voltages, can faithfully simulate short-term and long-term plasticity, and can even realize the triplet-spike-timing-dependent plasticity (triplet-STDP) rule, which can be used for pattern recognition in biological systems. The simulation results reveal that a memristor network of this type has an accuracy of ~95.78% in memory learning and the capability of pattern learning. This work provides a facile technology route to improve the performance of bionic-material memristors.

36 MATERIALS SCIENCE↗

NASA NTRS2016

Integrating Large Scale Data Sets to Develop Predictive Hypotheses of Low-Dose Radiation-Induced Health Effects

Over one hundred years of radiation biology research has revealed much about the DNA damages induced by the deposition of energy from exposure to ionizing radiation and the subsequent cellular responses. However, there are still significant gaps in our understanding of how these might lead to detrimental health effects, particularly at low doses (100 mGy (milligray)). Recent advances in high throughput omics technologies enable interrogation of induced radiation effects at the genomic, proteomic and metabolomic levels. These include changes in gene expression, protein modifications, e.g., phosphorylation, acetylation, and methylation, and metabolic changes. We will discuss the integration of data obtained from multiple omics platforms to understand radiation dose, and dose rate effects in a complex human tissue model as a function of time. We will use as an example our results on the low dose responses in a 3D human skin model.

ionizing radiation↗

DOE OSTI2014

The crystal structure of S. cerevisiae Sad1, a catalytically inactive deubiquitinase that is broadly required for pre-mRNA splicing

Sad1 is an essential splicing factor initially identified in a genetic screen in Saccharomyces cerevisiae for snRNP assembly defects. Based on sequence homology, Sad1, or USP39 in humans, is predicted to comprise two domains: a zinc finger ubiquitin binding domain (ZnF-UBP) and an inactive ubiquitin-specific protease (iUSP) domain, both of which are well conserved. The role of these domains in splicing and their interaction with ubiquitin are unknown. We first used splicing microarrays to analyze Sad1 function in vivo and found that Sad1 is critical for the splicing of nearly all yeast intron-containing genes. By using in vitro assays, we then showed that it is required for the assembly of the active spliceosome. To gain structural insights into Sad1 function, we determined the crystal structure of the full-length protein at 1.8 Å resolution. In the structure, the iUSP domain forms the characteristic ubiquitin binding pocket, though with an amino acid substitution in the active site that results in complete inactivation of the enzymatic activity of the domain. The ZnF-UBP domain of Sad1 shares high structural similarly to other ZnF-UBPs; however, Sad1’s ZnF-UBP does not possess the canonical ubiquitin binding motif. Given the precedents for ZnF-UBP domains to function as activators for their neighboring USP domains, we propose that Sad1’s ZnF-UBP acts in a ubiquitin-independent capacity to recruit and/or activate Sad1’s iUSP domain to interact with the spliceosome.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Integrative Modeling and Analysis of Fungal Central Carbon Metabolism

Over a thousand fungal genomes have been sequenced, yet manually curated genome-scale metabolic models (GEMs) are available for only a limited number of species. Moreover, these models have often been developed independently, leading to inconsistencies in namespaces, compartment definitions, and pathway representations that hinder comparative analysis, the systematic reuse of prior curation efforts, and the integration of consolidated metabolic knowledge. Here, we present the Consolidated Fungal Core Metabolism Model (CFCMM), constructed by integrating thirteen published fungal models spanning Ascomycota, Mucoromycota, and both Crabtree-positive and Crabtree-negative yeasts. We harmonized metabolites and reactions into a non-redundant shared ModelSEED ontological space, standardized compartmentalization, and refined gene–protein–reaction (GPR) rules. Using pathway-level visualization and systematic gap detection, we further improved the integrated network through literature-guided curation to correct stoichiometry, stereospecificity, and pathway architecture. Orthologous protein family reconstruction and functional annotation workflows were used to validate and inform GPR associations, with particular emphasis on ambiguous enzyme superfamilies and membrane-associated components. Using the resulting CFCMM, we built high-quality central carbon core models for each fungus and performed flux balance analysis to quantify ATP-yield variation under aerobic and anaerobic conditions, explicitly evaluating scenarios driven by differences in electron transport chain (ETC) composition. Simulations reproduced the expected fermentative yield of approximately 2 mmol ATP per mmol glucose under anaerobic conditions and separated the thirteen fungi into two bioenergetic groups under aerobic respiration based on Complex I status, with predicted yields of approximately 30 versus 22 mmol ATP per mmol glucose. Forcing flux through the alternative oxidase bypass further reduced ATP yields to approximately 12 and 4 mmol ATP per mmol glucose in Complex I-containing and Complex I-lacking fungi, respectively. Collectively, this work provides a manually curated, ModelSEED-consistent, and extensible fungal core metabolic template, deployed in DOE KBase as a resource for automated reconstruction of central carbon core models from any sequenced fungal genome. In addition, the CFCMM provides modular components for developing GEMs with more accurate energy predictions and enables robust comparative analyses of fungal bioenergetics and core metabolic diversity

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Structural Heterogeneity and Hydrodynamics of an Intrinsically Disordered Protein Condensate

Biology demonstrates precise control over the free-energy landscape through the selective partitioning of biomacromolecules into membraneless organelles, enabling essential functions such as biochemical transformations, signaling cascades, and mechanical reinforcement. Although the function of these condensates depends on their underlying structure and hydrodynamics, molecular-scale information on these systems remains sparse. Here, in this study, neutron scattering is used to probe the organization and dynamics of the intrinsically disordered N-terminal domain of Galectin-3, an extracellular lectin responsible for facilitating liquid–liquid phase separation on the cellular surface, in both dilute and condensed phases. Dilute solutions contain isolated protein chains in equilibrium with mesoscopic clusters, whereas the condensed phase adopts a bicontinuous, microemulsion-like morphology. The dilute phase behavior is quantitatively described by coarse-grained polymer models from soft-matter physics, demonstrating their predictive power for complex biological proteins. At elevated concentrations, the proteins self-assemble akin to block copolymers, microphase separating through the aggregation of hydrophobic domains along the protein contour. The resulting condensate remains fluid-like despite a 25-fold increase in concentration; its internal hydrodynamics slow by only a factor of 3 relative to dilute protein chains. These results provide a molecular-level framework for how disordered proteins achieve both the structural complexity and dynamic fluidity of biomolecular condensates.

Carrick, Brian R. [Massachusetts Inst. of Technolo↗

DOE OSTI2024

Characterization of a widespread sugar phosphate-processing bacterial microcompartment

Many prokaryotes form Bacterial Microcompartments (BMCs) that encapsulate segments of specialized metabolic pathways to enhance catalysis. The various functions of metabolosomes, catabolic BMCs, are dictated by the signature enzyme that processes initial substrates of the confined pathway. The components and native functions of several metabolosomes have been experimentally characterized; however one of the most prevalent across all bacteria has yet to be studied. Sugar Phosphate Utilizing (SPU) BMC loci encode enzymes predicted to be involved in sugar phosphate metabolism. The SPU genetic loci are found in organisms occupying habitats ranging from soils to hot springs, highlighting the ubiquity of the SPU BMC. We bioinformatically characterized seven SPU subtypes, all which contain an enzyme unique to SPU BMCs, a deoxyribose 5-phosphate aldolase (DERA). Here, we define the fundamental characteristics of SPU BMCs and have expressed, purified, and characterized a set of SPU core enzymes. These include a protein-protein complex formed between a SPU BMC DERA and a predicted ribose 5-phosphate isomerase. Further, we show that the SPU BMC DERA is catalytically active and propose that it acts as the universal signature enzyme for the SPU BMC, with implications for fundamental understanding and biotechnological applications of SPU BMCs.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

ZMPY3D: accelerating protein structure volume analysis through vectorized 3D Zernike moments and Python-based GPU integration

Abstract Motivation Volumetric 3D object analyses are being applied in research fields such as structural bioinformatics, biophysics, and structural biology, with potential integration of artificial intelligence/machine learning (AI/ML) techniques. One such method, 3D Zernike moments, has proven valuable in analyzing protein structures (e.g., protein fold classification, protein–protein interaction analysis, and molecular dynamics simulations). Their compactness and efficiency make them amenable to large-scale analyses. Established methods for deriving 3D Zernike moments, however, can be inefficient, particularly when higher order terms are required, hindering broader applications. As the volume of experimental and computationally-predicted protein structure information continues to increase, structural biology has become a “big data” science requiring more efficient analysis tools. Results This application note presents a Python-based software package, ZMPY3D, to accelerate computation of 3D Zernike moments by vectorizing the mathematical formulae and using graphical processing units (GPUs). The package offers popular GPU-supported libraries such as CuPy and TensorFlow together with NumPy implementations, aiming to improve computational efficiency, adaptability, and flexibility in future algorithm development. The ZMPY3D package can be installed via PyPI, and the source code is available from GitHub. Volumetric-based protein 3D structural similarity scores and transform matrix of superposition functionalities have both been implemented, creating a powerful computational tool that will allow the research community to amalgamate 3D Zernike moments with existing AI/ML tools, to advance research and education in protein structure bioinformatics. Availability and implementation ZMPY3D, implemented in Python, is available on GitHub (https://github.com/tawssie/ZMPY3D) and PyPI, released under the GPL License.

Lai, Jhih-Siang (ORCID:0000000156775890)↗

DOE OSTI2012

Improvements of biomass deconstruction enzymes

Sandia National Laboratories and DSM Innovation, Inc. collaborated on the investigation of the structure and function of cellulases from thermophilic fungi. Sandia's role was to use its expertise in protein structure determination and X-ray crystallography to solve the structure of these enzymes in their native state and in their substrate and product bound states. Sandia was also tasked to work with DSM to use the newly solved structure to, using computational approaches, analyze enzyme interactions with both bound substrate and bound product; the goal being to develop approaches for rationally designing improved cellulases for biomass deconstruction. We solved the structures of five cellulases from thermophilic fungi. Several of these were also solved with bound substrate/product, which allowed us to predict mutations that might enhance activity and stability.

09 BIOMASS FUELS↗

DOE OSTI2023

A High-Quality Genome-Scale Model for Rhodococcus opacus Metabolism

Rhodococcus opacus is a bacterium that has a high tolerance to aromatic compounds and can produce significant amounts of triacylglycerol (TAG). Here, we present iGR1773, the first genome-scale model (GSM) of R. opacus PD630 metabolism based on its genomic sequence and associated data. The model includes 1773 genes, 3025 reactions, and 1956 metabolites, was developed in a reproducible manner using CarveMe, and was evaluated through Metabolic Model tests (MEMOTE). We combine the model with two Constraint-Based Reconstruction and Analysis (COBRA) methods that use transcriptomics data to predict growth rates and fluxes: E-Flux2 and SPOT (Simplified Pearson Correlation with Transcriptomic data). Growth rates are best predicted by E-Flux2. Flux profiles are more accurately predicted by E-Flux2 than flux balance analysis (FBA) and parsimonious FBA (pFBA), when compared to 44 central carbon fluxes measured by 13C-Metabolic Flux Analysis (13C-MFA). Under glucose-fed conditions, E-Flux2 presents an R2 value of 0.54, while predictions based on pFBA had an inferior R2 of 0.28. We attribute this improved performance to the extra activity information provided by the transcriptomics data. For phenol-fed metabolism, in which the substrate first enters the TCA cycle, E-Flux2’s flux predictions display a high R2 of 0.96 while pFBA showed an R2 of 0.93. We also show that glucose metabolism and phenol metabolism function with similar relative ATP maintenance costs. These findings demonstrate that iGR1773 can help the metabolic engineering community predict aromatic substrate utilization patterns and perform computational strain design.

Roell, Garrett W.↗

DOE OSTI2023

HLA-Clus: HLA class I clustering based on 3D structure

In a previous paper, we classified populated HLA class I alleles into supertypes and subtypes based on the similarity of 3D landscape of peptide binding grooves, using newly defined structure distance metric and hierarchical clustering approach. Compared to other approaches, our method achieves higher correlation with peptide binding specificity, intra-cluster similarity (cohesion), and robustness. Here we introduce HLA-Clus, a Python package for clustering HLA Class I alleles using the method we developed recently and describe additional features including a new nearest neighbor clustering method that facilitates clustering based on user-defined criteria. The HLA-Clus pipeline includes three stages: First, HLA Class I structural models are coarse grained and transformed into clouds of labeled points. Second, similarities between alleles are determined using a newly defined structure distance metric that accounts for spatial and physicochemical similarities. Finally, alleles are clustered via hierarchical or nearest-neighbor approaches. We also interfaced HLA-Clus with the peptide:HLA affinity predictor MHCnuggets. By using the nearest neighbor clustering method to select optimal allele-specific deep learning models in MHCnuggets, the average accuracy of peptide binding prediction of rare alleles was improved. The HLA-Clus package offers a solution for characterizing the peptide binding specificities of a large number of HLA alleles. This method can be applied in HLA functional studies, such as the development of peptide affinity predictors, disease association studies, and HLA matching for grafting. HLA-Clus is freely available at our GitHub repository (https://github.com/yshen25/HLA-Clus).

59 BASIC BIOLOGICAL SCIENCES↗