Peptide biosensors for anticancer drugs: Design in silico to work in denaturizing environment

Replica-exchange optimization of antibody fragments

In the framework of the rational design of macromolecules capable of binding to a specific target for biosensing applications, we here further develop an evolutionary protocol designed to optimize the binding affinity of protein binders. In particular we focus on the optimization of the binding portion of small antibody fragments known as nanobodies (or VHH) and choose the hen egg white lysozyme (HEWL) as our target. By implementing a replica exchange scheme for this optimization, we show that an initial hit is not needed and similar solutions can be found by either optimizing an already known anti-HEWL VHH or a randomly selected binder (here a VHH selective towards another macromolecule). While we believe that exhaustive searches of the mutation space are most appropriate when only few key residues have to be optimized, in case a lead binder is not available the proposed evolutionary algorithm should be instead the method of choice.

Antibody Affinity Maturation Using Computational Methods: From an Initial Hit to Small-Scale Expression of Optimized Binders

Nanobodies (VHHs) are engineered fragments of the camelid single-chain immunoglobulins. The VHH domain contains the highly variable segments responsible for antigen recognition. VHHs can be easily produced as recombinant proteins. Their small size is a good advantage for in silico approaches. Computer methods represent a valuable strategy for the optimization and improvement of their binding affinity. They also allow for epitope selection offering the possibility to design new VHHs for regions of a target protein that are not naturally immunogenic. Here we present an in silico mutagenic protocol developed to improve the binding affinity of nanobodies together with the first step of their in vitro production. The method, already proven successful in improving the low Kd of a nanobody hit obtained by panning, can be employed for the ex novo design of antibody fragments against selected protein target epitopes.

Protocol for iterative optimization of modified peptides bound to protein targets

Peptides are commonly used as therapeutic agents. However, they suffer from easy degradation and instability. Replacing natural by non-natural amino acids can avoid these problems, and potentially improve the affinity towards the target protein. Here, we present a computational pipeline to optimize peptides based on adding non-natural amino acids while improving their binding affinity. The workflow is an iterative computational evolution algorithm, inspired by the PARCE protocol, that performs single-point mutations on the peptide sequence using modules from the Rosetta framework. The modifications can be guided based on the structural properties or previous knowledge of the biological system. At each mutation step, the affinity to the protein is estimated by sampling the complex conformations and applying a consensus metric using various open protein-ligand scoring functions. The mutations are accepted based on the score differences, allowing for an iterative optimization of the initial peptide. The sampling/scoring scheme was benchmarked with a set of protein-peptide complexes where experimental affinity values have been reported. In addition, a basic application using a known protein-peptide complex is also provided. The structure- and dynamic-based approach allows users to optimize bound peptides, with the option to personalize the code for further applications. The protocol, called mPARCE, is available at: https://github.com/rochoa85/mPARCE/.

Bicyclic peptide-based assay for uPA cancer biomarker

The use of synthetic bioreceptors to develop biosensing platforms has been recently gaining momentum. This case study compares the performance of a biosensing platform for the human biomarker urokinase-type plasminogen activator (h-uPA) when using two bicyclic peptides (P₁ and P₂) with different affinities for the target protein. The bioreceptors P₁ and P₂ were immobilized on magnetic microbeads and tested within a sandwich-type affinity electrochemical assay. Apart from enabling h-uPA quantification at nanomolar levels (105.8 ng/mL for P₁ and 32.5 ng/mL for P₂), this case study showed the potential of synthetic bicyclic peptides applicability and how bioreceptor affinity can influence the performance of the final sensing platform.

Computational Evolution Protocol for Peptide Design

Computational peptide design is useful for therapeutics, diagnostics, and vaccine development. To select the most promising peptide candidates, the key is describing accurately the peptide-target interactions at the molecular level. We here review a computational peptide design protocol whose key feature is the use of all-atom explicit solvent molecular dynamics for describing the different peptide-target complexes explored during the optimization. We describe the milestones behind the development of this protocol, which is now implemented in an open-source code called PARCE. We provide a basic tutorial to run the code for an antibody fragment design example. Finally, we describe three additional applications of the method to design peptides for different targets, illustrating the broad scope of the proposed approach.

Magnetic-Field-Driven Extraction of Bioreceptors into Polymeric Membranes for Label-Free Potentiometric Biosensing

We report here the concept of a magnetically controlled extraction of hydrophilic bioreceptors into polymeric membranes for bioassays. The potentiometric assay relies on the intrinsic charges of an antimicrobial peptide and its unique recognition abilities, which can eliminate the probe labeling and indicator addition. The target binding event could effectively prevent the extraction of the peptide into the polymeric membrane doped with an ion exchanger, thus resulting in a potential change. The potentiometric response properties of the peptide assembled on magnetic beads can be dynamically controlled and modulated by applying a magnetic field. Staphylococcus aureus, as a model of food-borne pathogens, was measured at levels down to 10 CFU mL^-1 . Based on this sensing strategy, a potentiometric array was developed for the pattern recognition of bacteria. The proposed general platform can be used for potentiometric biosensing using other hydrophilic bioreceptors.

Computational Evolution of Beta-2-Microglobulin Binding Peptides for Nanopatterned Surface Sensors

The bottom-up design of smart nanodevices largely depends on the accuracy by which each of the inherent nanometric components can be functionally designed with predictive methods. Here, we present a rationally designed, self-assembled nanochip capable of capturing a target protein by means of pre-selected binding sites. The sensing elements comprise computationally evolved peptides, designed to target an arbitrarily selected binding site on the surface of beta-2-Microglobulin (β2m), a globular protein that lacks well-defined pockets. The nanopatterned surface was generated by an atomic force microscopy (AFM)-based, tip force-driven nanolithography technique termed nanografting to construct laterally confined self-assembled nanopatches of single stranded (ss)DNA. These were subsequently associated with an ssDNA-peptide conjugate by means of DNA-directed immobilization, therefore allowing control of the peptide's spatial orientation. We characterized the sensitivity of such peptide-containing systems against β2m in solution by means of AFM-based differential topographic imaging and surface plasmon resonance (SPR) spectroscopy. Our results show that the confined peptides are capable of specifically capturing β2m from the surface-liquid interface with micromolar affinity, hence providing a viable proof-of-concept for our approach to peptide design.

Recent advances in computational methods for biosensor design

Biosensors are analytical tools with a great application in healthcare, food quality control, and environmental monitoring. They are of considerable interest to be designed by using cost-effective and efficient approaches. Designing biosensors with improved functionality or application in new target detection has been converted to a fast-growing field of biomedicine and biotechnology branches. Experimental efforts have led to valuable successes in the field of biosensor design; however, some deficiencies restrict their utilization for this purpose. Computational design of biosensors is introduced as a promising key to eliminate the gap. A set of reliable structure prediction of the biosensor segments, their stability, and accurate descriptors of molecular interactions are required to computationally design biosensors. In this review, we provide a comprehensive insight into the progress of computational methods to guide the design and development of biosensors, including molecular dynamics simulation, quantum mechanics calculations, molecular docking, virtual screening, and a combination of them as the hybrid methodologies. By relying on the recent advances in the computational methods, an opportunity emerged for them to be complementary or an alternative to the experimental methods in the field of biosensor design.

In silico design of peptides as potential ligands to resistin

Resistin is a hormone of biological interest due to its connection with several diseases of worldwide concern. This work aims to design a series of cyclic peptides as "lead compounds" to identify potential ligands to resistin. To this end, we propose an approach based on a peptide design algorithm plus a two-stage selection which accounts for selectivity, one of the most forgotten steps in the design of ligands. Following this approach, we have been able to identify several peptides as strong candidates for the design of elements of bio-recognition. Those peptides present low scoring binding energy to albumin, good water solubility, stability in water at 300 K, and high scoring binding energy to resistin. Among those peptides, two were chosen, to perform a more rigorous calculation of binding free energy based on the Alchemical Absolute Binding Free Energy method. We were able to establish a methodological route for the development of strong candidates for the design of ligands to resistin. Graphical Abstract Combined MD + MC + AABFE approach to design and screening of high-affinity binders to resistin.

A consensus protocol for the in silico optimisation of antibody fragments

We present an in silico mutagenetic protocol for improving the binding affinity of single domain antibodies (or nanobodies, VHHs). The method iteratively attempts random mutations in the interacting region of the protein and evaluates the resulting binding affinity towards the target by scoring, with a collection of scoring functions, short explicit solvent molecular dynamics trajectories of the binder-target complexes. The acceptance/rejection of each attempted mutation is carried out by a consensus decision-making algorithm, which considers all individual assessments derived from each scoring function. The method was benchmarked by evolving a single complementary determining region (CDR) of an anti-HER2 VHH hit obtained by direct panning of a phage display library. The optimised VHH mutant showed significantly enhanced experimental affinity with respect to the original VHH it matured from. The protocol can be employed as it is for the optimization of peptides, antibody fragments, and (given enough computational power) larger antibodies.

Predicting the Affinity of Peptides to Major Histocompatibility Complex Class II by Scoring Molecular Dynamics Simulations

Predicting the binding affinity of peptides able to interact with major histocompatibility complex (MHC) molecules is a priority for researchers working in the identification of novel vaccines candidates. Most available approaches are based on the analysis of the sequence of peptides of known experimental affinity. However, for MHC class II receptors, these approaches are not very accurate, due to the intrinsic flexibility of the complex. To overcome these limitations, we propose to estimate the binding affinity of peptides bound to an MHC class II by averaging the score of the configurations from finite-temperature molecular dynamics simulations. The score is estimated for 18 different scoring functions, and we explored the optimal manner for combining them. To test the predictions, we considered eight peptides of known binding affinity. We found that six scoring functions correlate with the experimental ranking of the peptides significantly better than the others. We then assessed a set of techniques for combining the scoring functions by linear regression and logistic regression. We obtained a maximum accuracy of 82% for the predicted sign of the binding affinity using a logistic regression with optimized weights. These results are potentially useful to improve the reliability of in silico protocols to design high-affinity binding peptides for MHC class II receptors.

FeptideDB: A web application for new bioactive peptides from food protein

BACKGROUND

Bioactive peptides derived from food are important sources for alternative medicine and possess therapeutic activity. Several biochemical methods have been achieved to isolate bioactive peptides from food, which are tedious and time consuming. In silico methods are an alternative process to reduce cost and time with respect to bioactive peptide production. In this paper, FeptideDB was used to collect bioactive peptide (BP) data from both published research articles and available bioactive peptide databases. FeptideDB was developed to assist in forecasting bioactive peptides from food by combining peptide cleavage tools and database matching. Furthermore, this application was able to predict the potential of cleaved peptides from 'enzyme digestion module' to identify new ACE (angiotensin converting enzyme) inhibitors using an automatic molecular docking approach.

RESULTS

The FeptideDB web application contains tools for generating all possible peptides cleaved from input protein by various available enzymes. This database was also used for analysis and visualization to assist in bioactive peptide discovery. One module of FeptideDB has the ability to create 3-dimensional peptide structures to further predict inhibitors for the target protein, ACE (angiotensin converting enzyme).

CONCLUSIONS

FeptideDB is freely available to researchers who are interested in exploring bioactive peptides. The FeptideDB interface is easy to use, allowing users to rapidly retrieve data based on desired search criteria. FeptideDB is freely available at http://www4g.biotec.or.th/FeptideDB/. Ultimately, FeptideDB is a computational aid for assessing peptide bioactivities.