{"doi":"10.1371/journal.pone.0275759","title":"Rosetta FlexPepDock to predict peptide-MHC binding: An approach for non-canonical amino acids","abstract":"Computation methods that predict the binding of peptides to MHC-I are important tools for screening and identifying immunogenic antigens and have the potential to accelerate vaccine and drug development. However, most available tools are sequence-based and optimized only for peptides containing the twenty canonical amino acids. This omits a large number of peptides containing non-canonical amino acids (NCAA), or residues that undergo varied post-translational modifications such as glycosylation or phosphorylation. These modifications fundamentally alter peptide immunogenicity. Similarly, existing structure-based methods are biased towards canonical peptide backbone structures, which may or may not be preserved when NCAAs are present. Rosetta FlexPepDock ab-initio is a structure-based computational protocol able to evaluate peptide-receptor interaction where no prior information of the peptide backbone is known. We benchmarked FlexPepDock ab-initio for docking canonical peptides to MHC-I, and illustrate for the first time the method's ability to accurately model MHC-I bound epitopes containing NCAAs. FlexPepDock ab-initio protocol was able to recapitulate near-native structures (≤1.5Å) in the top lowest-energy models for 20 out of 25 cases in our initial benchmark. Using known experimental binding affinities of twenty peptides derived from an influenza-derived peptide, we showed that FlexPepDock protocol is able to predict relative binding affinity as Rosetta energies correlate well with experimental values (r = 0.59, p = 0.006). ROC analysis revealed 80% true positive and a 40% false positive rate, with a prediction power of 93%. Finally, we demonstrate the protocol's ability to accurately recapitulate HLA-A*02:01 bound phosphopeptide backbone structures and relative binding affinity changes, the theoretical structure of the lymphocytic choriomeningitis derived glycosylated peptide GP392 bound to MHC-I H-2Db, and isolevuglandin-adducted peptides. The ability to use non-canonical amino acids in the Rosetta FlexPepDock protocol may provide useful insight into critical amino acid positions where the post-translational modification modulates immunologic responses.","journal":"PLoS ONE","year":2022,"id":254947,"datarank":0.4636563680037475,"base_score":3.091042453358316,"endowment":3.091042453358316,"self_citation_contribution":0.4636563680037475,"citation_network_contribution":0.0,"self_endowment_contribution":0.4636563680037475,"citer_contribution":0.0,"corpus_percentile":null,"corpus_rank":null,"citation_count":21,"citer_count":0,"citers_with_citation_signal":0,"citers_with_endowment":0,"datacite_reuse_total":0,"is_dataset":false,"is_dataset_confidence":0.9594,"is_data_producer":false,"deposit_databanks":null,"is_oa":true,"file_count":0,"downloads":0,"has_version_chain":false,"published_date":"2022-01-01","fair_score":null,"fair_percentile":null,"algorithm_id":"datarank_citation_only_1hop_v6","ranking_scope":"data_only","authors":[{"id":298178,"name":"Natália R. Barbaro","orcid":"0000-0003-1596-1085","position":1,"is_corresponding":false},{"id":36760,"name":"Rocco Moretti","orcid":"0000-0003-2162-1116","position":2,"is_corresponding":false},{"id":225956,"name":"David G. Harrison","orcid":"0000-0001-9363-6451","position":3,"is_corresponding":false},{"id":36757,"name":"Jens Meiler","orcid":"0000-0001-8945-193X","position":4,"is_corresponding":false},{"id":902605,"name":"Nathaniel Bloodworth","orcid":"0000-0002-5758-003X","position":0,"is_corresponding":true}],"reference_count":32,"raw_metadata":{"citation_network_status":"fetched"},"created_at":"2026-07-19T00:25:12.987829Z","pmid":"36512534","pmcid":null,"fwci":null,"citation_percentile":null,"influential_citations":0,"oa_status":null,"license":null,"views":0,"total_file_size_bytes":0,"version_count":0,"fair_f":null,"fair_a":null,"fair_i":null,"fair_r":null,"fair_zscore":null,"fair_rationale":null,"fair_model":null,"fair_agent_version":null,"fair_fulltext_source":null,"fair_has_llm":null,"fair_computed_at":null,"clinical_trials":[],"software_tools":[],"db_accessions":[],"linked_datasets":[],"topics":[]}