1
|
Alvarez JAE, Dean SN. TEMPRO: nanobody melting temperature estimation model using protein embeddings. Sci Rep 2024; 14:19074. [PMID: 39154093 PMCID: PMC11330463 DOI: 10.1038/s41598-024-70101-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2024] [Accepted: 08/13/2024] [Indexed: 08/19/2024] Open
Abstract
Single-domain antibodies (sdAbs) or nanobodies have received widespread attention due to their small size (~ 15 kDa) and diverse applications in bio-derived therapeutics. As many modern biotechnology breakthroughs are applied to antibody engineering and design, nanobody thermostability or melting temperature (Tm) is crucial for their successful utilization. In this study, we present TEMPRO which is a predictive modeling approach for estimating the Tm of nanobodies using computational methods. Our methodology integrates various nanobody biophysical features to include Evolutionary Scale Modeling (ESM) embeddings, NetSurfP3 structural predictions, pLDDT scores per sdAb region from AlphaFold2, and each sequence's physicochemical characteristics. This approach is validated with our combined dataset containing 567 unique sequences with corresponding experimental Tm values from a manually curated internal data and a recently published nanobody database, NbThermo. Our results indicate the efficacy of protein embeddings in reliably predicting the Tm of sdAbs with mean absolute error (MAE) of 4.03 °C and root mean squared error (RMSE) of 5.66 °C, thus offering a valuable tool for the optimization of nanobodies for various biomedical and therapeutic applications. Moreover, we have validated the models' performance using experimentally determined Tms from nanobodies not found in NbThermo. This predictive model not only enhances nanobody thermostability prediction, but also provides a useful perspective of using embeddings as a tool for facilitating a broader applicability of downstream protein analyses.
Collapse
Affiliation(s)
- Jerome Anthony E Alvarez
- Naval Research Laboratory, Center for Bio/Molecular Science and Engineering, Washington, DC, USA
| | - Scott N Dean
- Naval Research Laboratory, Center for Bio/Molecular Science and Engineering, Washington, DC, USA.
| |
Collapse
|
2
|
Sencanski M, Glisic S, Kubale V, Cotman M, Mavri J, Vrecl M. Computational Modeling and Characterization of Peptides Derived from Nanobody Complementary-Determining Region 2 (CDR2) Targeting Active-State Conformation of the β 2-Adrenergic Receptor (β 2AR). Biomolecules 2024; 14:423. [PMID: 38672440 PMCID: PMC11048008 DOI: 10.3390/biom14040423] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2024] [Revised: 03/20/2024] [Accepted: 03/28/2024] [Indexed: 04/28/2024] Open
Abstract
This study assessed the suitability of the complementarity-determining region 2 (CDR2) of the nanobody (Nb) as a template for the derivation of nanobody-derived peptides (NDPs) targeting active-state β2-adrenergic receptor (β2AR) conformation. Sequences of conformationally selective Nbs favoring the agonist-occupied β2AR were initially analyzed by the informational spectrum method (ISM). The derived NDPs in complex with β2AR were subjected to protein-peptide docking, molecular dynamics (MD) simulations, and metadynamics-based free-energy binding calculations. Computational analyses identified a 25-amino-acid-long CDR2-NDP of Nb71, designated P4, which exhibited the following binding free-energy for the formation of the β2AR:P4 complex (ΔG = -6.8 ± 0.8 kcal/mol or a Ki = 16.5 μM at 310 K) and mapped the β2AR:P4 amino acid interaction network. In vitro characterization showed that P4 (i) can cross the plasma membrane, (ii) reduces the maximum isoproterenol-induced cAMP level by approximately 40% and the isoproterenol potency by up to 20-fold at micromolar concentration, (iii) has a very low affinity to interact with unstimulated β2AR in the cAMP assay, and (iv) cannot reduce the efficacy and potency of the isoproterenol-mediated β2AR/β-arrestin-2 interaction in the BRET2-based recruitment assay. In summary, the CDR2-NDP, P4, binds preferentially to agonist-activated β2AR and disrupts Gαs-mediated signaling.
Collapse
Affiliation(s)
- Milan Sencanski
- Laboratory for Plant Molecular Biology, Institute of Molecular Genetics and Genetic Engineering, University of Belgrade, 11000 Belgrade, Serbia
- Laboratory for Bioinformatics and Computational Chemistry, Institute of Nuclear Sciences VINCA, National Institute of Serbia, University of Belgrade, 11000 Belgrade, Serbia;
| | - Sanja Glisic
- Laboratory for Bioinformatics and Computational Chemistry, Institute of Nuclear Sciences VINCA, National Institute of Serbia, University of Belgrade, 11000 Belgrade, Serbia;
| | - Valentina Kubale
- Institute of Preclinical Sciences, Veterinary Faculty, University of Ljubljana, 1000 Ljubljana, Slovenia; (V.K.); (M.C.)
| | - Marko Cotman
- Institute of Preclinical Sciences, Veterinary Faculty, University of Ljubljana, 1000 Ljubljana, Slovenia; (V.K.); (M.C.)
| | - Janez Mavri
- Department of Computational Biochemistry and Drug Design, National Institute of Chemistry, 1000 Ljubljana, Slovenia;
| | - Milka Vrecl
- Institute of Preclinical Sciences, Veterinary Faculty, University of Ljubljana, 1000 Ljubljana, Slovenia; (V.K.); (M.C.)
| |
Collapse
|
3
|
Cheng MH, Porritt RA, Rivas MN, Krieger JM, Ozdemir AB, Garcia G, Arumugaswami V, Fries BC, Arditi M, Bahar I. A monoclonal antibody against staphylococcal enterotoxin B superantigen inhibits SARS-CoV-2 entry in vitro. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2020. [PMID: 33269352 PMCID: PMC7709177 DOI: 10.1101/2020.11.24.395079] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
We recently discovered a superantigen-like motif, similar to Staphylococcal enterotoxin B (SEB), near the S1/S2 cleavage site of SARS-CoV-2 Spike protein, which might explain the multisystem-inflammatory syndrome (MIS-C) observed in children and cytokine storm in severe COVID-19 patients. We show here that an anti-SEB monoclonal antibody (mAb), 6D3, can bind this viral motif, and in particular its PRRA insert, to inhibit infection by blocking the access of host cell proteases, TMPRSS2 or furin, to the cleavage site. The high affinity of 6D3 for the furin-cleavage site originates from a poly-acidic segment at its heavy chain CDR2, a feature shared with SARS-CoV-2-neutralizing mAb 4A8. The affinity of 6D3 and 4A8 for this site points to their potential utility as therapeutics for treating COVID-19, MIS-C, or common cold caused by human coronaviruses (HCoVs) that possess a furin-like cleavage site.
Collapse
Affiliation(s)
- Mary Hongying Cheng
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| | - Rebecca A Porritt
- Department of Pediatrics, Division of Pediatric Infectious Diseases and Immunology, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA.,Biomedical Sciences, Infectious and Immunologic Diseases Research Center, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA
| | - Magali Noval Rivas
- Department of Pediatrics, Division of Pediatric Infectious Diseases and Immunology, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA.,Biomedical Sciences, Infectious and Immunologic Diseases Research Center, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA
| | - James M Krieger
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| | - Asli Beyza Ozdemir
- Department of Pediatrics, Division of Pediatric Infectious Diseases and Immunology, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA.,Biomedical Sciences, Infectious and Immunologic Diseases Research Center, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA
| | - Gustavo Garcia
- Department of Molecular and Medical Pharmacology, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA 90095, USA.,Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, University of California, Los Angeles, Los Angeles, CA90095, USA
| | - Vaithilingaraja Arumugaswami
- Department of Molecular and Medical Pharmacology, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA 90095, USA.,Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, University of California, Los Angeles, Los Angeles, CA90095, USA
| | - Bettina C Fries
- Department of Medicine, Stony Brook University Hospital, Stony Brook, New York, 11794, USA
| | - Moshe Arditi
- Department of Pediatrics, Division of Pediatric Infectious Diseases and Immunology, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA.,Biomedical Sciences, Infectious and Immunologic Diseases Research Center, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA
| | - Ivet Bahar
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| |
Collapse
|
4
|
Abstract
In this report, we used hydrogen exchange-mass spectrometry (HX-MS) to identify the epitopes recognized by 21 single-domain camelid antibodies (VHHs) directed against the ribosome-inactivating subunit (RTA) of ricin toxin, a biothreat agent of concern to military and public health authorities. The VHHs, which derive from 11 different B-cell lineages, were binned together based on competition ELISAs with IB2, a monoclonal antibody that defines a toxin-neutralizing hotspot ("cluster 3") located in close proximity to RTA's active site. HX-MS analysis revealed that the 21 VHHs recognized four distinct epitope subclusters (3.1-3.4). Sixteen of the 21 VHHs grouped within subcluster 3.1 and engage RTA α-helices C and G. Three VHHs grouped within subcluster 3.2, encompassing a-helices C and G, plus α-helix B. The single VHH in subcluster 3.3 engaged RTA α-helices B and G, while the epitope of the sole VHH defining subcluster 3.4 encompassed α-helices C and E, and β-strand h. Modeling these epitopes on the surface of RTA predicts that the 20 VHHs within subclusters 3.1-3.3 physically occlude RTA's active site cleft, while the single antibody in subcluster 3.4 associates on the active site's upper rim.
Collapse
|