1
|
Natsrita P, Charoenkwan P, Shoombuatong W, Mahalapbutr P, Faksri K, Chareonsudjai S, Rungrotmongkol T, Pipattanaboon C. Machine-learning-assisted high-throughput identification of potent and stable neutralizing antibodies against all four dengue virus serotypes. Sci Rep 2024; 14:17165. [PMID: 39060292 PMCID: PMC11282219 DOI: 10.1038/s41598-024-67487-8] [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: 12/15/2023] [Accepted: 07/11/2024] [Indexed: 07/28/2024] Open
Abstract
Several computational methods have been developed to identify neutralizing antibodies (NAbs) covering four dengue virus serotypes (DENV-1 to DENV-4); however, limitations of the dataset and the resulting performance remain. Here, we developed a new computational framework to predict potent and stable NAbs against DENV-1 to DENV-4 using only antibody (CDR-H3) and epitope sequences as input. Specifically, our proposed computational framework employed sequence-based ML and molecular dynamic simulation (MD) methods to achieve more accurate identification. First, we built a novel dataset (n = 1108) by compiling the interactions of CDR-H3 and epitope sequences with the half maximum inhibitory concentration (IC50) values, which represent neutralizing activities. Second, we achieved an accurately predictive ML model that showed high AUC values of 0.879 and 0.885 by tenfold cross-validation and independent tests, respectively. Finally, our computational framework could be applied to filter approximately 2.5 million unseen antibodies into two final candidates that showed strong and stable binding to all four serotypes. In addition, the most potent and stable candidate (1B3B9_V21) was evaluated for its development potential as a therapeutic agent by molecular docking and MD simulations. This study provides an antibody computational approach to facilitate the high-throughput identification of NAbs and accelerate the development of therapeutic antibodies.
Collapse
Affiliation(s)
- Piyatida Natsrita
- Department of Microbiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, 40002, Thailand
| | - Phasit Charoenkwan
- Modern Management and Information Technology, College of Arts, Media and Technology, Chiang Mai University, Chiang Mai, 50200, Thailand
| | - Watshara Shoombuatong
- Center for Research Innovation and Biomedical Informatics, Faculty of Medical Technology, Mahidol University, Bangkok, 10700, Thailand
| | - Panupong Mahalapbutr
- Department of Biochemistry, Faculty of Medicine, Khon Kaen University, Khon Kaen, 40002, Thailand
| | - Kiatichai Faksri
- Department of Microbiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, 40002, Thailand
- Research and Diagnostic Center for Emerging Infectious Diseases, Khon Kaen University, Khon Kaen, 40002, Thailand
| | - Sorujsiri Chareonsudjai
- Department of Microbiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, 40002, Thailand
| | - Thanyada Rungrotmongkol
- Center of Excellent in Biocatalyst and Sustainable Biotechnology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Chonlatip Pipattanaboon
- Department of Microbiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, 40002, Thailand.
| |
Collapse
|
2
|
Hernández González JE, de Araujo AS. Alchemical Calculation of Relative Free Energies for Charge-Changing Mutations at Protein-Protein Interfaces Considering Fixed and Variable Protonation States. J Chem Inf Model 2023; 63:6807-6822. [PMID: 37851531 DOI: 10.1021/acs.jcim.3c00972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2023]
Abstract
The calculation of relative free energies (ΔΔG) for charge-changing mutations at protein-protein interfaces through alchemical methods remains challenging due to variations in the system's net charge during charging steps, the possibility of mutated and contacting ionizable residues occurring in various protonation states, and undersampling issues. In this study, we present a set of strategies, collectively termed TIRST/TIRST-H+, to address some of these challenges. Our approaches combine thermodynamic integration (TI) with the prediction of pKa shifts to calculate ΔΔG values. Moreover, special sets of restraints are employed to keep the alchemically transformed molecules separated. The accuracy of the devised approaches was assessed on a large and diverse data set comprising 164 point mutations of charged residues (Asp, Glu, Lys, and Arg) to Ala at the protein-protein interfaces of complexes with known three-dimensional structures. Mean absolute and root-mean-square errors ranging from 1.38 to 1.66 and 1.89 to 2.44 kcal/mol, respectively, and Pearson correlation coefficients of ∼0.6 were obtained when testing the approaches on the selected data set using the GPU-TI module of Amber18 suite and the ff14SB force field. Furthermore, the inclusion of variable protonation states for the mutated acid residues improved the accuracy of the predicted ΔΔG values. Therefore, our results validate the use of TIRST/TIRST-H+ in prospective studies aimed at evaluating the impact of charge-changing mutations to Ala on the stability of protein-protein complexes.
Collapse
|
3
|
Beard H, Cholleti A, Pearlman D, Sherman W, Loving KA. Applying physics-based scoring to calculate free energies of binding for single amino acid mutations in protein-protein complexes. PLoS One 2013; 8:e82849. [PMID: 24340062 PMCID: PMC3858304 DOI: 10.1371/journal.pone.0082849] [Citation(s) in RCA: 164] [Impact Index Per Article: 14.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2013] [Accepted: 10/28/2013] [Indexed: 11/18/2022] Open
Abstract
Predicting changes in protein binding affinity due to single amino acid mutations helps us better understand the driving forces underlying protein-protein interactions and design improved biotherapeutics. Here, we use the MM-GBSA approach with the OPLS2005 force field and the VSGB2.0 solvent model to calculate differences in binding free energy between wild type and mutant proteins. Crucially, we made no changes to the scoring model as part of this work on protein-protein binding affinity--the energy model has been developed for structure prediction and has previously been validated only for calculating the energetics of small molecule binding. Here, we compare predictions to experimental data for a set of 418 single residue mutations in 21 targets and find that the MM-GBSA model, on average, performs well at scoring these single protein residue mutations. Correlation between the predicted and experimental change in binding affinity is statistically significant and the model performs well at picking "hotspots," or mutations that change binding affinity by more than 1 kcal/mol. The promising performance of this physics-based method with no tuned parameters for predicting binding energies suggests that it can be transferred to other protein engineering problems.
Collapse
Affiliation(s)
- Hege Beard
- Schrödinger, New York, New York, United States of America
| | | | - David Pearlman
- Schrödinger, New York, New York, United States of America
| | - Woody Sherman
- Schrödinger, New York, New York, United States of America
| | | |
Collapse
|
4
|
Abstract
ERM (Ezrin-Radixin-Moesin) proteins are key cross-linkers of the plasma membrane and the actin cytoskeleton. They are regulated by the intramolecular association of the N-terminal FERM (band-four point one, Ezrin, Radixin, Moesin) and C-terminal CERMAD (ERM association domain) domains (N/C interaction), which masks the binding surfaces of the domains for other molecules. The N/C interface is characterized by the highly distributed binding of CERMAD through a beta-strand and four alpha-helices to a globular FERM. Though it is a target for multiple regulatory signals, little is known about the dynamics/thermodynamics governing this interface. Recent implications of Ezrin in cancer metastasis have increased the necessity to understand this regulatory switch. In this study, we report residue-specific stabilities of Ezrin CERMAD at the Ezrin N/C interface obtained using hydrogen-deuterium exchange NMR. These stabilities vary across secondary structural elements and identify F583 and L586 as key anchor residues for the most stable element, alphaD. Macroscopic N/C binding energetics, obtained using isothermal titration calorimetry (ITC) reveals a high affinity (Kd =176 nM) enthalpy-driven binding (DeltaH = -26 kcal/mol, TDeltaS = -17 kcal/mol) at 25 degrees C at pH 7 in MES and phosphate buffers. A 10-fold increase in affinity was observed for measurements in acetate buffer, suggesting that an acetate-like molecule might promote the repressed form of the complex, possibly through interaction with the F2 subdomain of FERM, which resembles the acyl-CoA binding protein. In summary, our results have illustrated the dynamic nature of this regulatory interface and provide a foundation for investigating the role of regulatory signals on the stability of this interface.
Collapse
Affiliation(s)
- Bhargavi Jayaraman
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, New York 14853, USA
| | | |
Collapse
|
5
|
Fenalti G, Hampe CS, O'connor K, Banga JP, Mackay IR, Rowley MJ, El-Kabbani O. Molecular characterization of a disease associated conformational epitope on GAD65 recognised by a human monoclonal antibody b96.11. Mol Immunol 2007; 44:1178-89. [PMID: 16930708 DOI: 10.1016/j.molimm.2006.06.025] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2006] [Accepted: 06/20/2006] [Indexed: 10/24/2022]
Abstract
Autoantibodies to the 65kDa isoform of glutamate decarboxylase (GAD65) are associated with type I diabetes and recognise highly conformational epitope(s) that remain to be defined. The human recombinant Fab from mAb b96.11 inhibits binding of most GAD65 antibody positive sera from patients and its epitope has previously been localized to the middle region of GAD65. Recent studies indicate that b96.11 antibody specificity predicts the risk of developing type 1 diabetes in prediabetic individuals. We describe the use homology modelling, protein-protein docking simulations and biopanning of random peptide phage displayed libraries with b96.11 to predict contact amino acids on the interface of GAD65/Fab b96.11 complex. Further analysis by in vitro mutagenesis of GAD65 followed by radioimmunoprecipitation refined the amino acids contributing to the b96.11 epitope. Our studies show an interface characterized by a protruding antibody-combining site centered on the long heavy chain CDR3 loop of Fab b96.11 establishing interactions with the critical residue Phe(344) in the core of the epitope on GAD65, surrounded by charged sites within (375)RK(376) and (305)DER(307). The epitope requires residues from both middle and the C-terminal domains, and is the first precise definition of an epitope on GAD65. The nature of the b96.11 epitope leads to considerations of potential structural variations for differences in antigenicity between the isoforms GAD65 and GAD67. The study shows the utility of using a combination of in silico techniques and experimental data for molecular characterization and localization of conformational epitopes for which crystal structures are lacking.
Collapse
Affiliation(s)
- Gustavo Fenalti
- Department of Biochemistry and Molecular Biology, Monash University, Clayton, Victoria 3800, Australia
| | | | | | | | | | | | | |
Collapse
|
6
|
Streltsov VA, Carmichael JA, Nuttall SD. Structure of a shark IgNAR antibody variable domain and modeling of an early-developmental isotype. Protein Sci 2005; 14:2901-9. [PMID: 16199666 PMCID: PMC2253229 DOI: 10.1110/ps.051709505] [Citation(s) in RCA: 56] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
The new antigen receptor (IgNAR) antibodies from sharks are disulphide bonded dimers of two protein chains, each containing one variable and five constant domains. Three types of IgNAR variable domains have been discovered, with Type 3 appearing early in shark development and being overtaken by the antigen-driven affinity-matured Type 1 and 2 response. Here, we have determined the first structure of a naturally occurring Type 2 IgNAR variable domain, and identified the disulphide bond that links and stabilizes the CDR1 and CDR3 loops. This disulphide bridge locks the CDR3 loop in an "upright" conformation in contrast to other shark antibody structures, where a more lateral configuration is observed. Further, we sought to model the Type 3 isotype based on the crystallographic structure reported here. This modeling indicates (1) that internal Type 3-specific residues combine to pack into a compact immunoglobulin core that supports the CDR loop regions, and (2) that despite apparent low-sequence variability, there is sufficient plasticity in the CDR3 loop to form a conformationally diverse antigen-binding surface.
Collapse
Affiliation(s)
- Victor A Streltsov
- CSIRO Molecular and Health Technologies, 343 Royal Parade, Parkville, Victoria 3052, Australia.
| | | | | |
Collapse
|
7
|
NH···π Interactions Stabilize the Most-Hindered Rotamer of the(S)-Tyrosinato Side Group in Bis[(S)-tyrosinato](diamine)cobalt(III) Complexes: An NMR Spectroscopic and DFT Study. Eur J Inorg Chem 2005. [DOI: 10.1002/ejic.200400873] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
8
|
Abstract
Tyrosine residues of neuroendocrine peptides are frequently the targets of oxidation reactions, one of which involves hydroxylation to peptidyl-3, 4-dihydroxy-phenyl-L-alanine (DOPA). The reactivity in vitro of peptidyl-DOPA in two neuroendocrine peptides, a neurotensin fragment (pELYENK) and proctolin (RYLPT), was investigated using ultraviolet-visible scanning spectrophotometry and matrix-assisted laser desorption ionization mass spectrometry following oxidation by tyrosinase and periodate. The peptides form covalently coupled dimers and trimers, and their masses are consistent with the presence of diDOPA cross-links. Lysine does not appear to participate in multimer formation because it is efficiently recovered in fragmentation ladders using subtilisin. While multimer formation in the neurotensin-derived peptide can be blocked effectively by adding N-acetyl-DOPA-ethylester to the reaction medium, the DOPA ethylester couples itself four to five times to each peptide.
Collapse
Affiliation(s)
- L A Burzio
- Surgical Sealants, Inc., Woburn, Massachusetts 01801, USA.
| | | |
Collapse
|
9
|
Burgess K, Han I, Zhang A, Zheng WH, Shanmugam K, Featherstone MS, Saragovi HU. DiSSiMiL: Diverse Small Size Mini-Libraries applied to simple and rapid epitope mapping of a monoclonal antibody. THE JOURNAL OF PEPTIDE RESEARCH : OFFICIAL JOURNAL OF THE AMERICAN PEPTIDE SOCIETY 2001; 57:68-76. [PMID: 11168890 DOI: 10.1034/j.1399-3011.2001.00786.x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Methods for screening protein-protein interactions are useful in protein science and for the generation of drug leads. We set out to develop a simplified assay to rapidly test protein-protein interactions, with a library of 400 pentapeptides comprising the 20 natural amino acids at two variable positions followed by three glycines (NH2-X1X2GGG). The library was used to identify the epitope of monoclonal antibody (mAb) 10D11 directed against the HOXD4 protein. Three pentapeptide 'hits' were selected (VYGGG, PWGGG and WKGGG) from direct binding assays screening for pentapeptide-mAb interactions; and from assays using pentapeptides in solution to competitively block HOXD4-mAb interactions. Alignment of the three 'hit' pentapeptides to the HOXD4 sequence predicts the mAb 10D11 epitope as NH2-VYPWMK. Synthesis of NH2-VYPWMK hexapeptide confirmed this prediction; and an alanine scan of HOXD4 ablated binding by mAb 10D11 when amino acids in the putative epitope were mutated. We propose that these simplified but diverse libraries can be used for rapid epitope mapping of some mAbs, and for generating lead small peptide analogs that interfere with receptor-ligand or other protein-protein interactions, or with enzymatic activity.
Collapse
Affiliation(s)
- K Burgess
- Department Chemistry, Texas A&M University, College Station, USA
| | | | | | | | | | | | | |
Collapse
|
10
|
Lang S, Xu J, Stuart F, Thomas RM, Vrijbloed JW, Robinson JA. Analysis of antibody A6 binding to the extracellular interferon gamma receptor alpha-chain by alanine-scanning mutagenesis and random mutagenesis with phage display. Biochemistry 2000; 39:15674-85. [PMID: 11123892 DOI: 10.1021/bi000838z] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The monoclonal antibody A6 binds a conformational epitope comprising mainly the CC' surface loop on the N-terminal fibronectin type-III domain of the extracellular interferon gamma receptor (IFNgammaR). The crystal structure of an A6 Fab-IFNgammaR complex revealed an interface rich in the aromatic side chains of Trp, Tyr, and His residues. These aromatic side chains appear to interact with both polar and hydrophobic groups at the interface, a property which, in general, may be advantageous for ligand binding. To analyze these interactions in more detail, the affinities of 19 A6 alanine-scanning mutants for the IFNgammaR have been measured, using engineered A6 single chain variable region fragments, and a surface plasmon resonance biosensor. Energetically important side chains (DeltaG(mutant) - DeltaG(wt) > 2.4 kcal/mol), that form distinct hot spots in the binding interface, have been identified on both proteins. These include V(L)W92 in A6, whose benzenoid ring appears well situated for a pi-cation (or pi-amine) interaction with the side chain of receptor residue K47 and simultaneously for T-stacking onto the indole ring of W82 in the receptor. At another site, energetically important residues V(H)W52 and V(H)W53, as well as V(H)D54 and V(H)D56, surround the aliphatic side chain of the hot receptor residue K52. Taken together, the results show that side chains distributed across the interface, including many aromatic ones, make key energetic contributions to binding. In addition, the receptor CC' loop has been subjected to random mutagenesis, and receptor mutants with high affinity for A6 have been selected by phage display. Residues previously identified as important for receptor binding to A6 were conserved in the clones isolated. Some mutants, however, showed a much improved affinity for A6, due to changes at Glu55, a residue that appeared to be energetically unimportant for binding the antibody by alanine-scanning mutagenesis. An E55P receptor mutant bound A6 with a 600-fold increase in affinity (K(D) approximately 20 pM), which is one of the largest improvements in affinity from a single point mutation reported so far at any protein-protein interface.
Collapse
MESH Headings
- Alanine/genetics
- Amino Acid Sequence
- Antibodies, Monoclonal/biosynthesis
- Antibodies, Monoclonal/chemistry
- Antibodies, Monoclonal/genetics
- Antibodies, Monoclonal/metabolism
- Bacteriophage M13/genetics
- Bacteriophage M13/isolation & purification
- Base Sequence
- Binding Sites, Antibody/genetics
- Biosensing Techniques
- Circular Dichroism
- Cloning, Molecular
- Enzyme-Linked Immunosorbent Assay
- Extracellular Space/genetics
- Extracellular Space/metabolism
- Genomic Library
- Immunoglobulin Fragments/biosynthesis
- Immunoglobulin Fragments/chemistry
- Immunoglobulin Fragments/genetics
- Immunoglobulin Fragments/metabolism
- Molecular Sequence Data
- Mutagenesis, Insertional
- Mutagenesis, Site-Directed
- Receptors, Interferon/biosynthesis
- Receptors, Interferon/genetics
- Receptors, Interferon/immunology
- Receptors, Interferon/metabolism
- Interferon gamma Receptor
Collapse
Affiliation(s)
- S Lang
- Institute of Organic Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | | | | | | | | | | |
Collapse
|
11
|
Abstract
The application of surface plasmon resonance biosensors in life sciences and pharmaceutical research continues to increase. This review provides a comprehensive list of the commercial 1999 SPR biosensor literature and highlights emerging applications that are of general interest to users of the technology. Given the variability in the quality of published biosensor data, we present some general guidelines to help increase confidence in the results reported from biosensor analyses.
Collapse
Affiliation(s)
- R L Rich
- Center for Biomolecular Interaction Analysis, University of Utah School of Medicine, Salt Lake City 84132, USA
| | | |
Collapse
|
12
|
Jouanguy E, Dupuis S, Pallier A, Döffinger R, Fondanèche MC, Fieschi C, Lamhamedi-Cherradi S, Altare F, Emile JF, Lutz P, Bordigoni P, Cokugras H, Akcakaya N, Landman-Parker J, Donnadieu J, Camcioglu Y, Casanova JL. In a novel form of IFN-gamma receptor 1 deficiency, cell surface receptors fail to bind IFN-gamma. J Clin Invest 2000; 105:1429-36. [PMID: 10811850 PMCID: PMC315467 DOI: 10.1172/jci9166] [Citation(s) in RCA: 123] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/1999] [Accepted: 04/06/2000] [Indexed: 11/17/2022] Open
Abstract
Complete IFN-gamma receptor ligand-binding chain (IFNgammaR1) deficiency is a life-threatening autosomal recessive immune disorder. Affected children invariably die of mycobacterial infection, unless bone marrow transplantation is undertaken. Pathogenic IFNGR1 mutations identified to date include nonsense and splice mutations and frameshift deletions and insertions. All result in a premature stop codon upstream from the segment encoding the transmembrane domain, precluding cell surface expression of the receptors. We report herein two sporadic and two familial cases of a novel form of complete IFNgammaR1 deficiency in which normal numbers of receptors are detected at the cell surface. Two in-frame deletions and two missense IFNGR1 mutations were identified in the segment encoding the extracellular ligand-binding domain of the receptor. Eight independent IFNgammaR1-specific mAb's, including seven blocking antibodies, gave recognition patterns that differed between patients, suggesting that different epitopes were altered by the mutations. No specific binding of (125)I-IFN-gamma to cells was observed in any patient, however, and the cells failed to respond to IFN-gamma. The mutations therefore cause complete IFNgammaR1 deficiency by disrupting the IFN-gamma-binding site without affecting surface expression. The detection of surface IFNgammaR1 molecules by specific antibodies, including blocking antibodies, does not exclude a diagnosis of complete IFNgammaR1 deficiency.
Collapse
Affiliation(s)
- E Jouanguy
- Laboratoire de Génétique Humaine des Maladies Infectieuses, Faculté de Médecine Necker, Paris, France
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|