Protein-protein Complexes Research Articles

Engineering the Mechanical Stability of a Therapeutic Complex between Affibody and Programmed Death-Ligand 1 by Anchor Point Selection.

Protein-protein complexes can vary in mechanical stability depending on the direction from which force is applied. Here, we investigated the mechanical stability of a complex between a binding scaffold called Affibody and an immune checkpoint protein Programmed Death-Ligand 1 (PD-L1). We used AFM single-molecule force spectroscopy with bioorthogonal clickable peptide handles, shear stress bead adhesion assays, molecular modeling, and steered molecular dynamics (SMD) to understand the pulling point dependency of the mechanostability of the Affibody:(PD-L1) complex. We observed a wide range of rupture forces depending on the anchor point. Pulling from residue #22 on Affibody generated an intermediate state attributed to partially unfolded PD-L1, while pulling from Affibody's N-terminus generated a force-activated catch bond. Pulling from residue #22 or #47 on Affibody generated high rupture forces, with the complex breaking at up to ∼190 pN under loading rates of ∼104-105 pN/s, representing a ∼4-fold increase as compared with low-force N-terminal pulling. SMD simulations showed relative tendencies in rupture forces that were consistent with experiments and, through visualization of force propagation networks, provided mechanistic insights. These results demonstrate how the mechanical properties of protein-protein interfaces can be controlled by informed choice of site-specific bioconjugation points within molecules, with implications for optimal bioconjugation strategies in drug delivery vehicles.

Integration of molecular coarse-grained model into geometric representation learning framework for protein-protein complex property prediction

Structure-based machine learning algorithms have been utilized to predict the properties of protein-protein interaction (PPI) complexes, such as binding affinity, which is critical for understanding biological mechanisms and disease treatments. While most existing algorithms represent PPI complex graph structures at the atom-scale or residue-scale, these representations can be computationally expensive or may not sufficiently integrate finer chemical-plausible interaction details for improving predictions. Here, we introduce MCGLPPI, a geometric representation learning framework that combines graph neural networks (GNNs) with MARTINI molecular coarse-grained (CG) models to predict PPI overall properties accurately and efficiently. Extensive experiments on three types of downstream PPI property prediction tasks demonstrate that at the CG-scale, MCGLPPI achieves competitive performance compared with the counterparts at the atom- and residue-scale, but with only a third of computational resource consumption. Furthermore, CG-scale pre-training on protein domain-domain interaction structures enhances its predictive capabilities for PPI tasks. MCGLPPI offers an effective and efficient solution for PPI overall property predictions, serving as a promising tool for the large-scale analysis of biomolecular interactions.

Chemical Tools for Probing the Ub/Ubl Conjugation Cascades.

Conjugation of ubiquitin (Ub) and structurally related ubiquitin-like proteins (Ubls), essential for many cellular processes, employs multi-step reactions orchestrated by specific E1, E2 and E3 enzymes. The E1 enzyme activates the Ub/Ubl C-terminus in an ATP-dependent process that results in the formation of a thioester linkage with the E1 active site cysteine. The thioester-activated Ub/Ubl is transferred to the active site of an E2 enzyme which then interacts with an E3 enzyme to promote conjugation to the target substrate. The E1-E2-E3 enzymatic cascades utilize labile intermediates, extensive conformational changes, and vast combinatorial diversity of short-lived protein-protein complexes to conjugate Ub/Ubl to various substrates in a regulated manner. In this review, we discuss various chemical tools and methods used to study the consecutive steps of Ub/Ubl activation and conjugation, which are often too elusive for direct studies. We focus on methods developed to probe enzymatic activities and capture and characterize stable mimics of the transient intermediates and transition states, thereby providing insights into fundamental mechanisms in the Ub/Ubl conjugation pathways.

Mechanisms of Copper Selectivity and Release by the Metallochaperone CusF: Insights from CO-Binding, Rapid-Freeze-Quench EXAFS, and Unnatural Amino Acid Substitution.

Metallochaperones are small proteins that shuttle essential metal ions such as Cu selectively to their cellular targets. CusF has unusual Cu(I) coordination, bound by two methionines, one histidine and a capping tryptophan residue, W44. Here we compare the CO binding reactivity of the wild type (WT) protein and its W44A, F, and M variants. Fourier-transform infrared (FTIR) indicates that W44A provides unhindered access for CO, while W44M is unreactive. WT is also largely unreactive to CO suggesting that the tryptophan cap is effective in shielding the Cu(I) center from exogenous adduct formation, while the Phe variant shows partial reactivity suggestive of an equilibrium between cap-on and cap-off conformers. Rates of metal transfer to the partner CusB are consistent with the π-cation cap providing both selectivity and redox protection. Unnatural amino acid substitutions of the W44 ligand with cyano-Phe and Br-Phe underpin the conclusion that the Phe ligand is a less effective capping residue. Finally, density functional theory (DFT) calculations validate the CO-binding strategy. Overall, the study suggests that CusF uses the tryptophan cap to protect against exogenous ligand (O2) attack while the mechanism of protein-protein complex formation allows the cap to swing out of the way, and thus have minimal effect on the rates of metal transfer.

Assessment of Binding Affinity in the Complexes of CoV-S-Protein’s RBD and the ACE2 Using Convolutional Neural Networks

The experimentally obtained structures of 48 complexes of the ACE2 receptor with the S protein’s RBD of the coronaviruses SARS-CoV and SARS-CoV-2 (including mutant forms of the latter) were assessed and the dissociation constant was calculated for them. Prediction of binding affinity was carried out using ProBAN, a neural network algorithm, previously developed by the authors, and a number of other algorithms for estimating the Gibbs free energy such as Prodigy, FoldX, DFIRE and RosettaDock. A comparison of the evaluation results shows that ProBAN has the best prediction quality (Pearson correlation − 0.56, MAE − 0.7 kcal/mol) of all the analyzed algorithms. The results obtained suggest better quality of affinity prediction for other protein-protein complexes. Information about the complexes under study and prediction results are available in the repository at the link: https://github.com/EABogdanova/ProBAN_RBD-ACE2.

Plastocyanin and Cytochrome f Complex Structures Obtained by NMR, Molecular Dynamics, and AlphaFold 3 Methods Compared to Cryo-EM Data.

Plastocyanin is a small mobile protein that facilitates electron transfer through the formation of short-lived protein-protein complexes with cytochrome bf and photosystem 1. Due to the transient nature of plastocyanin-cytochrome f complex, the lack of a long-lived tight complex makes it impossible to determine its structure by X-ray diffraction analysis. Up to today, a number of slightly different structures of such complexes have been obtained by experimental and computer methods. Now, artificial intelligence gives us the possibility to predict the structures of intermolecular complexes. In this study, we compare encounter and final complexes obtained by Brownian and molecular dynamics methods, as well as the structures predicted by AlphaFold 3, with NMR and cryo-EM data. Surprisingly, the best match for the plastocyanin electron density obtained by cryo-EM was demonstrated by an AlphaFold 3 structure. The orientation of plastocyanin in this structure almost completely coincides with its orientation obtained by molecular dynamics calculation, and, at the same time, it is different from the orientation of plastocyanin predicted on the basis of NMR data. This is even more unexpected given that only NMR structures for the plastocyanin-cytochrome f complex are available in the PDB database, which was used to train AlphaFold 3.

Proteome-wide copy-number estimation from transcriptomics

Protein copy numbers constrain systems-level properties of regulatory networks, but proportional proteomic data remain scarce compared to RNA-seq. We related mRNA to protein statistically using best-available data from quantitative proteomics and transcriptomics for 4366 genes in 369 cell lines. The approach starts with a protein’s median copy number and hierarchically appends mRNA–protein and mRNA–mRNA dependencies to define an optimal gene-specific model linking mRNAs to protein. For dozens of cell lines and primary samples, these protein inferences from mRNA outmatch stringent null models, a count-based protein-abundance repository, empirical mRNA-to-protein ratios, and a proteogenomic DREAM challenge winner. The optimal mRNA-to-protein relationships capture biological processes along with hundreds of known protein-protein complexes, suggesting mechanistic relationships. We use the method to identify a viral-receptor abundance threshold for coxsackievirus B3 susceptibility from 1489 systems-biology infection models parameterized by protein inference. When applied to 796 RNA-seq profiles of breast cancer, inferred copy-number estimates collectively re-classify 26–29% of luminal tumors. By adopting a gene-centered perspective of mRNA–protein covariation across different biological contexts, we achieve accuracies comparable to the technical reproducibility of contemporary proteomics.

MPA-MutPred: a novel strategy for accurately predicting the binding affinity change upon mutation in membrane protein complexes.

Mutations in the interface of membrane protein (MP) complexes are key contributors to a broad spectrum of human diseases, primarily due to changes in their binding affinities. While various methods exist for predicting the mutation-induced changes in binding affinity (ΔΔG) in protein-protein complexes, none are specific to MP complexes. This study proposes a novel strategy for ΔΔG prediction in MP complexes, which combines linear and nonlinear models, to obtain a more robust model with improved prediction accuracy. We used multiple linear regression to extract informative features that influence the binding affinity in MP complexes, which included changes in the stability of the complex, conservation score, electrostatic interaction, relatively accessible surface area, and interface contacts. Further, using gradient boosting regressor on the selected features, we developed MPA-MutPred, a novel method specific for predicting the ΔΔG of membrane protein-protein complexes, and it is freely accessible at https://web.iitm.ac.in/bioinfo2/MPA-MutPred/. Our method achieved a correlation of 0.75 and a mean absolute error (MAE) of 0.73kcal/mol in the jack-knife test conducted on a dataset of 770 mutants. We further validated the method using a blind test set of 86 mutations, obtaining a correlation of 0.85 and an MAE of 0.77kcal/mol. We anticipate that this method can be used for large-scale studies to understand the influence of binding affinity change on disease-causing mutations in MP complexes, thereby aiding in the understanding of disease mechanisms and the identification of potential therapeutic targets.

Computational screening of the effects of mutations on protein-protein off-rates and dissociation mechanisms by τRAMD

The dissociation rate, or its reciprocal, the residence time (τ), is a crucial parameter for understanding the duration and biological impact of biomolecular interactions. Accurate prediction of τ is essential for understanding protein-protein interactions (PPIs) and identifying potential drug targets or modulators for tackling diseases. Conventional molecular dynamics simulation techniques are inherently constrained by their limited timescales, making it challenging to estimate residence times, which typically range from minutes to hours. Building upon its successful application in protein-small molecule systems, τ-Random Acceleration Molecular Dynamics (τRAMD) is here investigated for estimating dissociation rates of protein-protein complexes. τRAMD enables the observation of unbinding events on the nanosecond timescale, facilitating rapid and efficient computation of relative residence times. We tested this methodology for three protein-protein complexes and their extensive mutant datasets, achieving good agreement between computed and experimental data. By combining τRAMD with MD-IFP (Interaction Fingerprint) analysis, dissociation mechanisms were characterized and their sensitivity to mutations investigated, enabling the identification of molecular hotspots for selective modulation of dissociation kinetics. In conclusion, our findings underscore the versatility of τRAMD as a simple and computationally efficient approach for computing relative protein-protein dissociation rates and investigating dissociation mechanisms, thereby aiding the design of PPI modulators.

Effect of Direct and Water-Mediated Interactions on the Identification of Hotspots in Biomolecular Complexes with Multiple Subsystems.

Identification of important residues in biochemical complexes is often a crucial step for many problems in molecular biology and biochemistry. A method is proposed to identify hotspots in biomolecular complexes based on a new metric, derived from networks representing molecular subunits (residues, bridging solvent molecules, ligands etc.) connected by interactions. A singular value decomposition of the weighted adjacency matrix is used to construct a scalar rank for each subunit that reflects its importance in the residue interaction network. This metric is called the singular value centrality. In addition, a new formalism is proposed to account for water-mediated interactions in the ranking of residues. Interactions for a residue network can be provided by various computational methods. In this work interactions are obtained from full quantum-mechanical calculations of protein-protein complexes using the fragment molecular orbital method. The ranking results are shown to be in good agreement with earlier computational and experimental studies. The developed method can be used to gain a deeper insight into the role of subunits in complex biomolecular systems.

Machine learning, network pharmacology, and molecular dynamics reveal potent cyclopeptide inhibitors against dengue virus proteins.

The dengue virus is a major global health hazard responsible for an estimated 390 million diseases yearly. This study focused on identifying cyclopeptide inhibitors for envelope structural proteins E, NS1, NS3, and NS5. Additionally, 5579 cyclopeptides were individually screened against the four target proteins using a machine learning-based quantitative structure-activity relationship model. Subsequently, the best 10 cyclopeptides from each protein were selected for molecular docking with their corresponding proteins. Moreover, the protein-peptide complexes with the highest affinity were subjected to a 100-ns molecular dynamics simulation. The protein-protein complexes exhibited superior structural stability and binding interactions. Based on the results of the MD simulation analyses, which included checking values for Root Mean Square Deviation, Root Mean Square Fluctuation, Principal Component Analysis (PCA), free energy landscapes, and energetic components, it was found that NS5-CP03714 complex is more stable and has stronger binding interactions than NS3-CP02054. PCA and free energy landscape plots have confirmed the higher conformational stability of NS5-CP03714. Analysis of the energetic components revealed that NS5-CP03714 (total binding energy = -47.19kcal/mol) exhibits more favorable interaction energies and overall binding energy compared to NS3-CP02054 (total binding energy = -27.36kcal/mol), suggesting a stronger and more stable formation of the complex. In addition, the drug-target network of two specific peptides (CP02950 and CP05582) and their associated target proteins were analyzed. This analysis revealed valuable information about their ability to target several proteins and their potential for broad-spectrum activity. Additional experimental investigations are necessary to validate these computational results and assess the efficacy of identified peptide inhibitors in biological systems.

SAAMBE-MEM: a sequence-based method for predicting binding free energy change upon mutation in membrane protein-protein complexes.

Mutations in protein-protein interactions can affect the corresponding complexes, impacting function and potentially leading to disease. Given the abundance of membrane proteins, it is crucial to assess the impact of mutations on the binding affinity of these proteins. Although several methods exist to predict the binding free energy change due to mutations in protein-protein complexes, most require structural information of the protein complex and are primarily trained on the SKEMPI database, which is composed mainly of soluble proteins. A novel sequence-based method (SAAMBE-MEM) for predicting binding free energy changes (ΔΔG) in membrane protein-protein complexes due to mutations has been developed. This method utilized the MPAD database, which contains binding affinities for wild-type and mutant membrane protein complexes. A machine learning model was developed to predict ΔΔG by leveraging features such as amino acid indices and position-specific scoring matrices (PSSM). Through extensive dataset curation and feature extraction, SAAMBE-MEM was trained and validated using the XGBoost regression algorithm. The optimal feature set, including PSSM-related features, achieved a Pearson correlation coefficient of 0.64, outperforming existing methods trained on the SKEMPI database. Furthermore, it was demonstrated that SAAMBE-MEM performs much better when utilizing evolution-based features in contrast to physicochemical features. The method is accessible via a web server and standalone code at http://compbio.clemson.edu/SAAMBE-MEM/. The cleaned MPAD database is available at the website.

Evaluation of AlphaFold 3's Protein-Protein Complexes for Predicting Binding Free Energy Changes upon Mutation.

AlphaFold 3 (AF3), the latest version of protein structure prediction software, goes beyond its predecessors by predicting protein-protein complexes. It could revolutionize drug discovery and protein engineering, marking a major step toward comprehensive, automated protein structure prediction. However, independent validation of AF3's predictions is necessary. In this work, we evaluate AF3 complex structures using the SKEMPI 2.0 database which involves 317 protein-protein complexes and 8338 mutations. AF3 complex structures when applied to the most advanced TDL model, MT-TopLap (MultiTask-Topological Laplacian), give rise to a very good Pearson correlation coefficient of 0.86 for predicting protein-protein binding free energy changes upon mutation, which is slightly less than the 0.88 achieved earlier with the Protein Data Bank (PDB) structures. Nonetheless, AF3 complex structures led to a 8.6% increase in the prediction RMSE compared to original PDB complex structures. Additionally, some of AF3's complex structures have large errors, which were not captured in its ipTM performance metric. Finally, it is found that AF3's complex structures are not reliable for intrinsically flexible regions or domains.

Quantitative Characterization of the Impact of Protein-Protein Interactions on Ligand-Protein Binding: A Multi-Chain Dynamics Perturbation Analysis Method.

Many protein-protein interactions (PPIs) affect the ways in which small molecules bind to their constituent proteins, which can impact drug efficacy and regulatory mechanisms. While recent advances have improved our ability to independently predict both PPIs and ligand-protein interactions (LPIs), a comprehensive understanding of how PPIs affect LPIs is still lacking. Here, we examined 63 pairs of ligand-protein complexes in a benchmark dataset for protein-protein docking studies and quantified six typical effects of PPIs on LPIs. A multi-chain dynamics perturbation analysis method, called mcDPA, was developed to model these effects and used to predict small-molecule binding regions in protein-protein complexes. Our results illustrated that the mcDPA can capture the impact of PPI on LPI to varying degrees, with six similar changes in its predicted ligand-binding region. The calculations showed that 52% of the examined complexes had prediction accuracy at or above 50%, and 55% of the predictions had a recall of not less than 50%. When applied to 33 FDA-approved protein-protein-complex-targeting drugs, these numbers improved to 60% and 57% for the same accuracy and recall rates, respectively. The method developed in this study may help to design drug-target interactions in complex environments, such as in the case of protein-protein interactions.

Do All Paths Lead to Rome? How Reliable is Umbrella Sampling Along a Single Path?

Molecular dynamics (MD) simulations are widely applied to estimate absolute binding free energies of protein-ligand and protein-protein complexes. A routinely used method for binding free energy calculations with MD is umbrella sampling (US), which calculates the potential of mean force (PMF) along a single reaction coordinate. Surprisingly, in spite of its widespread use, few validation studies have focused on the convergence of the free energy computed along a single path for specific cases, not addressing the reproducibility of such calculations in general. In this work, we therefore investigate the reproducibility and convergence of US along a standard distance-based reaction coordinate for various protein-protein and protein-ligand complexes, following commonly used guidelines for the setup. We show that repeating the complete US workflow can lead to differences of 2-20 kcal/mol in computed binding free energies. We attribute those discrepancies to small differences in the binding pathways. While these differences are unavoidable in the established US protocol, the popularity of the latter could hint at a lack of awareness of such reproducibility problems. To test if the convergence of PMF profiles can be improved if multiple pathways are sampled simultaneously, we performed additional simulations with an adaptive-biasing method, here the accelerated weight histogram (AWH) approach. Indeed, the PMFs obtained from AHW simulations are consistent and reproducible for the systems tested. To the best of our knowledge, our work is the first to attempt a systematic assessment of the pitfalls in one the most widely used protocols for computing binding affinities. We anticipate therefore that our results will provide an incentive for a critical reassessment of the validity of PMFs computed with US, and make a strong case to further benchmark the performance of adaptive-biasing methods for computing binding affinities.

Machine-learning-based structural analysis of interactions between antibodies and antigens

Computational analysis of paratope-epitope interactions between antibodies and their corresponding antigens can facilitate our understanding of the molecular mechanism underlying humoral immunity and boost the design of new therapeutics for many diseases. The recent breakthrough in artificial intelligence has made it possible to predict protein-protein interactions and model their structures. Unfortunately, detecting antigen-binding sites associated with a specific antibody is still a challenging problem. To tackle this challenge, we implemented a deep learning model to characterize interaction patterns between antibodies and their corresponding antigens. With high accuracy, our model can distinguish between antibody-antigen complexes and other types of protein-protein complexes. More intriguingly, we can identify antigens from other common protein binding regions with an accuracy of higher than 70% even if we only have the epitope information. This indicates that antigens have distinct features on their surface that antibodies can recognize. Additionally, our model was unable to predict the partnerships between antibodies and their particular antigens. This result suggests that one antigen may be targeted by more than one antibody and that antibodies may bind to previously unidentified proteins. Taken together, our results support the precision of antibody-antigen interactions while also suggesting positive future progress in the prediction of specific pairing.

ProBASS - a language model with sequence and structural features for predicting the effect of mutations on binding affinity.

Protein-protein interactions (PPIs) govern virtually all cellular processes. Even a single mutation within PPI can significantly influence overall protein functionality and potentially lead to various types of diseases. To date, numerous approaches have emerged for predicting the change in free energy of binding (ΔΔGbind) resulting from mutations, yet the majority of these methods lack precision. In recent years, protein language models (PLMs) have been developed and shown powerful predictive capabilities by leveraging both sequence and structural data from protein-protein complexes. Yet, PLMs have not been optimized specifically for predicting ΔΔGbind. We developed an approach to predict effects of mutations on PPI binding affinity based on two most advanced protein language models ESM2 and ESM-IF1 that incorporate PPI sequence and structural features, respectively. We used the two models to generate embeddings for each PPI mutant and subsequently fine-tuned our model by training on a large dataset of experimental ΔΔGbind values. Our model, ProBASS (Protein Binding Affinity from Structure and Sequence) achieved a correlation with experimental ΔΔGbind values of 0.83 ± 0.05 for single mutations and 0.69 ± 0.04 for double mutations when model training and testing was done on the same PDB. Moreover, ProBASS exhibited very high correlation (0.81 ± 0.02) between prediction and experiment when training and testing was performed on a dataset containing 2325 single mutations in 132 PPIs. ProBASS surpasses the state-of-the-art methods in correlation with experimental data and could be further trained as more experimental data becomes available. Our results demonstrate that the integration of extensive datasets containing ΔΔGbind values across multiple PPIs to refine the pre-trained PLMs represents a successful approach for achieving a precise and broadly applicable model for ΔΔGbind prediction, greatly facilitating future protein engineering and design studies.

Exploring the Conformational Ensembles of Protein-Protein Complex with Transformer-Based Generative Model.

Protein-protein interactions are the basis of many protein functions, and understanding the contact and conformational changes of protein-protein interactions is crucial for linking the protein structure to biological function. Although difficult to detect experimentally, molecular dynamics (MD) simulations are widely used to study the conformational ensembles and dynamics of protein-protein complexes, but there are significant limitations in sampling efficiency and computational costs. In this study, a generative neural network was trained on protein-protein complex conformations obtained from molecular simulations to directly generate novel conformations with physical realism. We demonstrated the use of a deep learning model based on the transformer architecture to explore the conformational ensembles of protein-protein complexes through MD simulations. The results showed that the learned latent space can be used to generate unsampled conformations of protein-protein complexes for obtaining new conformations complementing pre-existing ones, which can be used as an exploratory tool for the analysis and enhancement of molecular simulations of protein-protein complexes.

Engineering the Mechanical Stability of a Therapeutic Affibody/PD-L1 Complex by Anchor Point Selection.

Protein-protein complexes can vary in mechanical stability depending on the direction from which force is applied. Here we investigated the anisotropic mechanical stability of a molecular complex between a therapeutic non-immunoglobulin scaffold called Affibody and the extracellular domain of the immune checkpoint protein PD-L1. We used a combination of single-molecule AFM force spectroscopy (AFM-SMFS) with bioorthogonal clickable peptide handles, shear stress bead adhesion assays, molecular modeling, and steered molecular dynamics (SMD) simulations to understand the pulling point dependency of mechanostability of the Affibody:(PD-L1) complex. We observed diverse mechanical responses depending on the anchor point. For example, pulling from residue #22 on Affibody generated an intermediate unfolding event attributed to partial unfolding of PD-L1, while pulling from Affibody's N-terminus generated force-activated catch bond behavior. We found that pulling from residue #22 or #47 on Affibody generated the highest rupture forces, with the complex breaking at up to ~ 190 pN under loading rates of ~104-105 pN/sec, representing a ~4-fold increase in mechanostability as compared with low force N-terminal pulling. SMD simulations provided consistent tendencies in rupture forces, and through visualization of force propagation networks provided mechanistic insights. These results demonstrate how mechanostability of therapeutic protein-protein interfaces can be controlled by informed selection of anchor points within molecules, with implications for optimal bioconjugation strategies in drug delivery vehicles.

Binding Mechanism between Platelet Glycoprotein and Cyclic Peptide Elucidated by McMD-Based Dynamic Docking.

The cyclic peptide OS1 (amino acid sequence: CTERMALHNLC), which has a disulfide bond between both termini cysteine residues, inhibits complex formation between the platelet glycoprotein Ibα (GPIbα) and the von Willebrand factor (vWF) by forming a complex with GPIbα. To study the binding mechanism between GPIbα and OS1 and, therefore, the inhibition mechanism of the protein-protein GPIbα-vWF complex, we have applied our multicanonical molecular dynamics (McMD)-based dynamic docking protocol starting from the unbound state of the peptide. Our simulations have reproduced the experimental complex structure, although the top-ranking structure was an intermediary one, where the peptide was bound in the same location as in the experimental structure; however, the β-switch of GPIbα attained a different conformation. Our analysis showed that subsequent refolding of the β-switch results in a more stable binding configuration, although the transition to the native configuration appears to take some time, during which OS1 could dissociate. Our results show that conformational changes in the β-switch are crucial for successful binding of OS1. Furthermore, we identified several allosteric binding sites of GPIbα that might also interfere with vWF binding, and optimization of the peptide to target these allosteric sites might lead to a more effective inhibitor, as these are not dependent on the β-switch conformation.