The development of methods able to modulate the binding affinity between proteins and peptides is of paramount biotechnological interest in view of a vast range of applications that imply designed polypeptides capable to impair or favour Protein-Protein Interactions. Here, we applied a peptide design algorithm based on shape complementarity optimization and electrostatic compatibility and provided the first experimental in vitro proof of the efficacy of the design algorithm. Focusing on the interaction between the SARS-CoV-2 Spike Receptor-Binding Domain (RBD) and the human angiotensin-converting enzyme 2 (ACE2) receptor, we extracted a 23-residues long peptide that structurally mimics the major interacting portion of the ACE2 receptor and designed in silico five mutants of such a peptide with a modulated affinity. Remarkably, experimental KD measurements, conducted using biolayer interferometry, matched the in silico predictions. Moreover, we investigated the molecular determinants that govern the variation in binding affinity through molecular dynamics simulation, by identifying the mechanisms driving the different values of binding affinity at a single residue level. Finally, the peptide sequence with the highest affinity, in comparison with the wild type peptide, was expressed as a fusion protein with human H ferritin (HFt) 24-mer. Solution measurements performed on the latter constructs confirmed that peptides still exhibited the expected trend, thereby enhancing their efficacy in RBD binding. Altogether, these results indicate the high potentiality of this general method in developing potent high-affinity vectors for hindering/enhancing protein-protein associations.

Numerous aptamers against various targets have been identified through the technology of systematic evolution of ligands by exponential enrichment (SELEX), but the affinity of these aptamers are often insufficient due to the limitations of SELEX. Therefore, a more rational in silico screening strategy (ISS) was developed for efficient screening of high affinity aptamers, which took shape complementarity and thermodynamic stability into consideration. Neuron specific enolase (NSE), a tumor marker, was selected as the target molecule. In the screening process, three aptamer candidates with good shape complementarity, lower ΔG values, and higher ZDOCK scores were produced. The dissociation constant (Kd) of these candidates to NSE was determined to be 10.13 nM, 14.82 nM, and 2.76 nM, respectively. Each of them exhibited higher affinity to NSE than the parent aptamer (Kd = 23.83 nM). Finally, an antibody-free fluorescence aptasensor assay, based on the aptamer with the highest affinity, P-5C8G, was conducted, resulting in a limit of detection (LOD) value of 1.8 nM, which was much lower than the parental aptamer (P, LOD = 12.6 nM). The proposed ISS approach provided an efficient and universal strategy to improve the aptamer to have a high affinity and good analytical utility.

Many factors influence biomolecule binding, and its assessment constitutes an elusive challenge in computational structural biology. In this aspect, the evaluation of shape complementarity at molecular interfaces is one of the main factors to be considered. We focus on the particular case of antibody-antigen complexes to quantify the complementarities occurring at molecular interfaces. We relied on a method we recently developed, which employs the 2D Zernike descriptors, to characterize the investigated regions with an ordered set of numbers summarizing the local shape properties. Collecting a structural dataset of antibody-antigen complexes, we applied this method and we statistically distinguished, in terms of shape complementarity, pairs of the interacting regions from the non-interacting ones. Thus, we set up a novel computational strategy based on in silico mutagenesis of antibody-binding site residues. We developed a Monte Carlo procedure to increase the shape complementarity between the antibody paratope and a given epitope on a target protein surface. We applied our protocol against several molecular targets in SARS-CoV-2 spike protein, known to be indispensable for viral cell invasion. We, therefore, optimized the shape of template antibodies for the interaction with such regions. As the last step of our procedure, we performed an independent molecular docking validation of the results of our Monte Carlo simulations.

The structure of a protein plays a pivotal role in determining its function. Often, the protein surface\'s shape and curvature dictate its nature of interaction with other proteins and biomolecules. However, marked by corrugations and roughness, a protein\'s surface representation poses significant challenges for its curvature-based characterization. In the present study, we employ unsupervised machine learning to segment the protein surface into patches. To measure the surface curvature of a patch, we present an algebraic sphere fitting method that is fast, accurate, and robust. Moreover, we use local curvatures to show the existence of \"shape complementarity\" in protein-protein, antigen-antibody, and protein-ligand interfaces. We believe that the current approach could help understand the relationship between protein structure and its biological function and can be used to find binding partners of a given protein.

We report a design strategy for integrative assembly of heteromeric [2]catenanes. The design focuses on the shape and functional group match of two different metalla-rectangles. A series of dipyridyl ligands with different lengths, widths and functional groups were designed and used for assembly experiments. Six heteromeric [2]catenanes were obtained both by direct mixture of two pre-assembled metalla-rectangles and one-pot three-component self-assembly. Multiple analytic methods were employed to characterize the catenanes, including single crystal X-ray diffraction analysis, NMR spectroscopy, mass spectroscopy and elemental analysis.

We propose a computational investigation on the interaction mechanisms between SARS-CoV-2 spike protein and possible human cell receptors. In particular, we make use of our newly developed numerical method able to determine efficiently and effectively the relationship of complementarity between portions of protein surfaces. This innovative and general procedure, based on the representation of the molecular isoelectronic density surface in terms of 2D Zernike polynomials, allows the rapid and quantitative assessment of the geometrical shape complementarity between interacting proteins, which was unfeasible with previous methods. Our results indicate that SARS-CoV-2 uses a dual strategy: in addition to the known interaction with angiotensin-converting enzyme 2, the viral spike protein can also interact with sialic-acid receptors of the cells in the upper airways.

Despite the huge effort to contain the infection, the novel SARS-CoV-2 coronavirus has rapidly become pandemic, mainly due to its extremely high human-to-human transmission capability, and a surprisingly high viral charge of symptom-less people. While the seek for a vaccine is still ongoing, promising results have been obtained with antiviral compounds. In particular, lactoferrin is regarded to have beneficial effects both in preventing and soothing the infection. Here, we explore the possible molecular mechanisms with which lactoferrin interferes with SARS-CoV-2 cell invasion, preventing attachment and/or entry of the virus. To this aim, we search for possible interactions lactoferrin may have with virus structural proteins and host receptors. Representing the molecular iso-electron surface of proteins in terms of 2D-Zernike descriptors, we 1) identified putative regions on the lactoferrin surface able to bind sialic acid present on the host cell membrane, sheltering the cell from the virus attachment; 2) showed that no significant shape complementarity is present between lactoferrin and the ACE2 receptor, while 3) two high complementarity regions are found on the N- and C-terminal domains of the SARS-CoV-2 spike protein, hinting at a possible competition between lactoferrin and ACE2 for the binding to the spike protein.

We present a method for efficiently and effectively assessing whether and where two proteins can interact with each other to form a complex. This is still largely an open problem, even for those relatively few cases where the 3D structure of both proteins is known. In fact, even if much of the information about the interaction is encoded in the chemical and geometric features of the structures, the set of possible contact patches and of their relative orientations are too large to be computationally affordable in a reasonable time, thus preventing the compilation of reliable interactome. Our method is able to rapidly and quantitatively measure the geometrical shape complementarity between interacting proteins, comparing their molecular iso-electron density surfaces expanding the surface patches in term of 2D Zernike polynomials. We first test the method against the real binding region of a large dataset of known protein complexes, reaching a success rate of 0.72. We then apply the method for the blind recognition of binding sites, identifying the real region of interaction in about 60 % of the analyzed cases. Finally, we investigate how the efficiency in finding the right binding region depends on the surface roughness as a function of the expansion order.

Expansion of the amino-acid repertoire with synthetic derivatives introduces novel structures and functionalities into proteins. In this study, we improved the antigen binding of antibodies by incorporating halogenated tyrosines at multiple selective sites. Tyrosines in the Fab fragment of an anti-EGF-receptor antibody 059-152 were systematically replaced with 3-bromo- and 3-chlorotyrosines, and simultaneous replacements at four specific sites were found to cause a tenfold increase in the affinity toward the antigen. Structure modeling suggested that this effect was due to enhanced shape complementarity between the antigen and antibody molecules. On the other hand, we showed that chlorination in the constant domain, far from the binding interface, of Rituximab Fab also increased the affinity significantly (up to 17-fold). Our results showed that antigen binding is tunable with the halogenation in and out of the binding motifs.

BACKGROUND: Protein-protein docking is a valuable computational approach for investigating protein-protein interactions. Shape complementarity is the most basic component of a scoring function and plays an important role in protein-protein docking. Despite significant progresses, shape representation remains an open question in the development of protein-protein docking algorithms, especially for grid-based docking approaches.
RESULTS: We have proposed a new pairwise shape-based scoring function (LSC) for protein-protein docking which adopts an exponential form to take into account long-range interactions between protein atoms. The LSC scoring function was incorporated into our FFT-based docking program and evaluated for both bound and unbound docking on the protein docking benchmark 4.0. It was shown that our LSC achieved a significantly better performance than four other similar docking methods, ZDOCK 2.1, MolFit/G, GRAMM, and FTDock/G, in both success rate and number of hits. When considering the top 10 predictions, LSC obtained a success rate of 51.71% and 6.82% for bound and unbound docking, respectively, compared to 42.61% and 4.55% for the second-best program ZDOCK 2.1. LSC also yielded an average of 8.38 and 3.94 hits per complex in the top 1000 predictions for bound and unbound docking, respectively, followed by 6.38 and 2.96 hits for the second-best ZDOCK 2.1.
CONCLUSIONS: The present LSC method will not only provide an initial-stage docking approach for post-docking processes but also have a general implementation for accurate representation of other energy terms on grids in protein-protein docking. The software has been implemented in our HDOCK web server at http://hdock.phys.hust.edu.cn/.