Design, synthesis and biological evaluation of PD-1 derived peptides as inhibitors of PD-1/PD-L1 complex formation for cancer therapy.

Over the past few years, many molecules such as monoclonal antibodies, affibodies, nanobodies, and small compounds have been designed and tested as inhibitors of PD-1/PD-L1 complex formation. Some of them have been successfully implemented into clinical oncology practice. However, the majority of these compounds have disadvantages and limitations, such as high production price, potential for immunogenicity and/or prolonged clearance. Thus, new inhibitors of the PD-1/PD-L1 immune checkpoints are needed. Recently, peptides emerged as potential novel approach for blocking receptor/ligand interaction. In the presented studies we have designed, synthesised and tested peptides, which are potential inhibitors of the PD-1/PD-L1 axis. The amino acid sequences of the designed peptides were based on the binding sites of PD-1 to PD-L1, as determined by the crystal structure of the protein complex and also based on MM/GBSA analysis. Interactions of the peptides with PD-L1 protein were confirmed using SPR, while their inhibitory properties were studied using cell-based PD-1/PD-L1 immune checkpoint blockade assays. The characterization of the peptides has shown that the peptides PD-1(119-142)T120C-E141C, PD-1(119-142)C123-S137C and PD-1(122-138)C123-S137C strongly bind to PD-L1 protein and disrupt the interaction of the proteins. PD-1(122-138)C123-S137C peptide was shown to have the best inhibitory potential from the panel of peptides. Its 3D NMR structure was determined and the binding site to PD-L1 was established using molecular modelling methods. Our results indicate that the PD-1 derived peptides are able to mimic the PD-1 protein and inhibit PD-1/PD-L1 complex formation.

Impact of calcium ions on the structural and dynamic properties of heparin oligosaccharides by computational analysis.

Heparin (HP) belongs to glycosaminoglycans (GAGs), anionic linear polysaccharides composed of repetitive disaccharide units. They are key players in many biological processes occurring in the extracellular matrix and at the cell surface. GAGs are challenging molecules for computational research due to their high chemical heterogeneity, flexibility, periodicity, pseudosymmetry, predominantly electrostatics-driven nature of interactions with their protein partners and potential multipose binding. The molecular mechanisms underlying GAG interactions mediated by divalent ions, which are important for GAG binding to several proteins, are not well understood. The goal of this study was to characterize the binding of Ca2+ to two HP oligosaccharides of different lengths (dp10 and dp18, dp: degree of polymerization) and their impact on HP conformational space and their dynamic behavior with the use of molecular dynamics (MD)-based approaches with two Ca2+ parameter sets. MD data suggested that the flexibility of the monosaccharides, the glycosidic linkages and ring puckering were not affected by the presence of Ca2+, in contrast to H-bond propensities and the calculated Rg for a fraction of the oligosaccharide populations in both dp10 and dp18. Moreover, the essential differences in the data obtained by using two Ca2+ parameter sets were reported.

Molecular dynamics-based descriptors of 3-O-Sulfated Heparan sulfate as contributors of protein binding specificity.

Glycosaminoglycans are linear periodic and anionic polysaccharides found in the extracellular matrix, involved in a range of key biochemical processes as a result of their interactions with a variety of protein partners. Due to the template-less synthesis, high flexibility and charge of GAGs, as well as the multipose binding of GAG ligands to receptors, the specificity of GAG-protein interactions can be difficult to elucidate. In this study we propose a set of MD-based descriptors of unbound Heparan Sulfate hexasaccharides that can be used to characterize GAGs and explain their binding affinity to a set of protein receptors. With the help of experimental data on GAG-protein binding affinity, we were able to further characterize the nature of this interaction in addition to providing a basis for predictor functions of GAG-protein binding specificity.

Modeling glycosaminoglycan-protein complexes.

Glycosaminoglycans are long linear and complex polysaccharides that are fundamental components of the mammalian extracellular matrix. Therefore, it is crucial to appropriately characterize molecular structure, dynamics, and interactions of protein-glycosaminoglycans complexes for improving understanding of molecular mechanisms underlying GAG biological function. Nevertheless, this proved challenging experimentally, and theoretical techniques are beneficial to construct new hypotheses and aid the interpretation of experimental data. The scope of this mini-review is to summarize four specific aspects of the current theoretical approaches for investigating noncovalent protein-glycosaminoglycan complexes such as molecular docking, free binding energy calculations, modeling ion impact, and addressing the phenomena of multipose binding of glycosaminoglycans to proteins.

Effect of Tetraphenylborate on Physicochemical Properties of Bovine Serum Albumin.

The binding interactions of bovine serum albumin (BSA) with tetraphenylborate ions ([B(Ph)4]-) have been investigated by a set of experimental methods (isothermal titration calorimetry, steady-state fluorescence spectroscopy, differential scanning calorimetry and circular dichroism spectroscopy) and molecular dynamics-based computational approaches. Two sets of structurally distinctive binding sites in BSA were found under the experimental conditions (10 mM cacodylate buffer, pH 7, 298.15 K). The obtained results, supported by the competitive interactions experiments of SDS with [B(Ph)4]- for BSA, enabled us to find the potential binding sites in BSA. The first site is located in the subdomain I A of the protein and binds two [B(Ph)4]- ions (logK(ITC)1 = 7.09 ± 0.10; ΔG(ITC)1 = -9.67 ± 0.14 kcal mol-1; ΔH(ITC)1 = -3.14 ± 0.12 kcal mol-1; TΔS(ITC)1 = -6.53 kcal mol-1), whereas the second site is localized in the subdomain III A and binds five ions (logK(ITC)2 = 5.39 ± 0.06; ΔG(ITC)2 = -7.35 ± 0.09 kcal mol-1; ΔH(ITC)2 = 4.00 ± 0.14 kcal mol-1; TΔS(ITC)2 = 11.3 kcal mol-1). The formation of the {[B(Ph)4]-}-BSA complex results in an increase in the thermal stability of the alfa-helical content, correlating with the saturation of the particular BSA binding sites, thus hindering its thermal unfolding.

Modeling protein structures with the coarse-grained UNRES force field in the CASP14 experiment.

The UNited RESidue (UNRES) force field was tested in the 14th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP14), in which larger oligomeric and multimeric targets were present compared to previous editions. Three prediction modes were tested (i) ab initio (the UNRES group), (ii) contact-assisted (the UNRES-contact group), and (iii) template-assisted (the UNRES-template group). For most of the targets, the contact restraints were derived from the server models top-ranked by the DeepQA method, while the DNCON2 method was used for 11 targets. Our consensus-fragment procedure was used to run template-assisted predictions. Each group also processed the Nuclear Magnetic Resonance (NMR)- and Small Angle X-Ray Scattering (SAXS)-data assisted targets. The average Global Distance Test Total Score (GDT_TS) of the 'Model 1' predictions were 29.17, 39.32, and 56.37 for the UNRES, UNRES-contact, and UNRES-template predictions, respectively, increasing by 0.53, 2.24, and 3.76, respectively, compared to CASP13. It was also found that the GDT_TS of the UNRES models obtained in ab initio mode and in the contact-assisted mode decreases with the square root of chain length, while the exponent in this relationship is 0.20 for the UNRES-template group models and 0.11 for the best performing AlphaFold2 models, which suggests that incorporation of database information, which stems from protein evolution, brings in long-range correlations, thus enabling the correction of force-field inaccuracies.

Evaluation of replica exchange with repulsive scaling approach for docking glycosaminoglycans.

Glycosaminoglycans (GAGs), long linear periodic anionic polysaccharides, are key molecules in the extracellular matrix (ECM). Therefore, deciphering their role in the biologically relevant context is important for fundamental understanding of the processes ongoing in ECM and for establishing new strategies in the regenerative medicine. Although GAGs represent a number of computational challenges, molecular docking is a powerful tool for analysis of their interactions. Despite the recent development of GAG-specific docking approaches, there is plenty of room for improvement. Here, replica exchange molecular dynamics with repulsive scaling (REMD-RS) recently proved to be a successful approach for protein-protein complexes, was applied to dock GAGs. In this method, effective pairwise radii are increased in different Hamiltonian replicas. REMD-RS is shown to be an attractive alternative to classical docking approaches for GAGs. This work contributes to setting up of GAG-specific computational protocols and provides new insights into the nature of these biological systems.

Computational insights into the role of calcium ions in protein-glycosaminoglycan systems.

Glycosaminoglycans (GAGs) are anionic, periodic, linear polysaccharides which are composed of periodic disaccharide units. They play a vital role in many biological processes ongoing in the extracellular matrix. In terms of computational approaches, GAGs are very challenging molecules due to their high flexibility, periodicity, predominantly electrostatic-driven nature of interactions with their protein counterparts and potential multipose binding. Furthermore, the molecular mechanisms underlying GAG-mediated interactions are not fully known yet, and experimental techniques alone are not always sufficient to gain insights into them. The aim of this study was to characterize protein-ion-GAG complexes for the systems where ions are directly involved in GAG binding. Molecular docking, molecular dynamics and free energy calculation approaches were applied to model and rigorously analyse the interactions between annexins (II and V), calcium ions (Ca2+) and heparin (HP). The computational data were examined and discussed in the context of the structural data previously reported by the crystallographic studies. The computational results confirm that the presence of Ca2+ has a tremendous impact on the annexin-HP binding site. This study provides a general computational pipeline to discover the complexity of protein-GAG interactions and helps to understand the role of ions involved at the atomic level. The limitations of the applied protocols are described and discussed pointing at the challenges persisting in the state-of-the-art in silico tools to study protein-ion-GAG systems.