Affinity enhancement of nanobody binding to EGFR: in silico site-directed mutagenesis and molecular dynamics simulation approaches

CDR grafting and site-directed mutagenesis approach for the generation and affinity maturation of Anti-CD20 nanobody

BACKGROUND

Recently, new and advanced techniques have been adopted to design and produce nanobodies, which are used in diagnostic and immunotherapy treatments. Traditionally, nanobodies are prepared from camelid immune libraries that require animal treatments. However, such approaches require large library sizes and complicated selection procedures. The current study has employed CDR grafting and site-directed mutagenesis techniques to create genetically engineered nanobodies against the tumor marker CD20 (anti-CD20 nanobodies) used in leukemia treatment.

METHODS AND RESULTS

In this study, we utilized the swapping method to graft CDRs from the VH Rituximab antibody to VHH CDRs. We aimed to enhance the binding affinity of the nanobodies by substituting the amino acids (Y101R-Y102R-Y107R) in the VHH-CDR3. To assess the binding capacity of the mutated nanobodies, we conducted an ELISA test. Moreover, through flow cytometry analysis, we compared the fluorescence intensity of the grafted CD20 and mutant nanobodies with that of the commercially available human anti-CD20 in Raji cells. The results showed a significant difference in the fluorescence intensity of the grafted nanobodies and mutant nanobodies when compared to the commercially available human anti-CD20.

CONCLUSION

The approach we followed in this study makes it possible to create multiple anti-CD20 nanobodies with varying affinities without the need for extensive selection efforts. Additionally, our research has demonstrated that computational tools are highly reliable in designing functional nanobodies.

New insights into the inhibitory effect of phenol carboxylic acid antioxidants on mushroom tyrosinase by molecular dynamic studies and experimental assessment

The inhibitory effects of ferulic and chlorogenic acids on tyrosinase activity were investigated through multi-spectroscopic and molecular docking techniques. Ferulic and chlorogenic acids, flavonoid compounds, demonstrated inhibitory monophenolase activities of tyrosinase. The inhibitor effects against monophenolase activity were in a reversible and competitive manner with ki value equal to 6.8 and 7.5 µM respectively. The affinity between tyrosinase and L-DOPA decreased when fatty acids were added to the solution. The multi-spectroscopic techniques like UV-vis, fluorescence, and isothermal calorimetry are employed to investigate changes. Intrinsic fluorescence quenching and conformational changes of tyrosinase by hydrophobic interaction were confirmed. Tyrosinase had two and three binding sites for ferulic and chlorogenic acids with a binding constant in the order of magnitude of -6.8 and -7.2 kcal/mol. In addition, the secondary structural changes with Circular dichroism (CD) analysis, secondary structure (DSSP), radius of gyration (Rg) and analysis of hydrogen bonds (H-bonds) confirmed. Ferulic acid effect can be observed obviously and also content of α-helix decreased. Thermodynamic parameters indicated that the interaction between enzyme and ferulic and chlorogenic acids followed a spontaneous reaction dynamic manner with ΔG = -14.78 kJ/mol and ΔG = -14.61 kJ/mol (298k). The findings highlighted the potential applications of ferulic acid and chlorogenic acids in food and drug industries as potent inhibitors of tyrosinase.Communicated by Ramaswamy H. Sarma.

Affinity enhancement of CR3022 binding to RBD; in silico site directed mutagenesis using molecular dynamics simulation approaches

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a disease which caused by a novel beta coronavirus. Structural and non-structural proteins are expressed by the virus gene fragments. The RBD of the S1 protein of the virus has the ability to interact with potent antibodies including CR3022, which was characterized to target the S protein of the virus which can efficiently neutralize the SARS-CoV in vitro and in vivo. In current study, we aimed to design CR3022 based antibody with high affinity compared with wild-type CR3022 using MD simulation method. Two variants were designed based on the amino acid binding conformation and the free binding energy of the critical amino acids which involved in CR3022-RBD interactions were evaluated. In this study three complexes were evaluated; CR3022-RBD, V1-RBD and V2-RBD using molecular dynamics simulations carried out for 100 ns in each case. Then, all the complexes were simulated for 100 ns. In the next step, to calculate the free binding affinity of the wild CR3022 and mutant antibody (V1 and V2) with RBD, the PMF method was performed. The RMSD profile demonstrated that all three complexes were equilibrated after 85 ns. Furthermore, the free binding energy results indicated that the V2-RBD complex has the higher binding affinity than V1-RBD and CR3022-RBD complexes. It should be noted that in above variants, the electrostatic energy and the number of H-bonds between the antibody and RBD increased. Thus, it is suggested that both designed antibodies could be considered as appropriate candidates for covid-19 disease treatment.Communicated by Ramaswamy H. Sarma.

In silico approach to probe the binding affinity between OMVs harboring the ZEGFR affibody and the EGF receptor

There is a growing interest in designing a nanocarrier containing an EGFR targeting affibody to direct toward cancer cells. Here, cytolysin A was cloned at the N-terminus of Z_EGFR:1907 affibody to guarantee its surface presentation on the OMVs while targeting the epidermal growth factor receptors (EGFRs). A separate construct including a fusogenic peptide (GALA) was also designed for the endosomal escape of the nanocarrier. Binding of the two constructs ClyA-affi_EGFR and ClyA-affi_EGFR-GALA to domain III of EGFR was investigated using molecular docking and molecular dynamic simulations. The higher stability of the ClyA-affi_EGFR-GALA/EGFR as compared to the ClyA-affi_EGFR/EGFR complex was evident. The ClyA-affi_EGFR-GALA structure showed a higher RMSD during the first half of the simulation time implying a much less stable behavior. Plateau state of the radius of gyration plot of ClyA-affi_EGFR-GALA confirmed a well-folded structure in the presence of the GALA sequence. Solvent accessible surface area for both proteins was in the same range. The data obtained from hydrogen bond analysis revealed a more equilibrated and stable form of the ClyA-affi_EGFR-GALA structure upon interaction with EGFR. The data provided here was a requisite for our biological evaluation of the synthesized constructs as a component of a novel drug delivery system.

Enfuvirtide, an HIV-1 fusion inhibitor peptide, can act as a potent SARS-CoV-2 fusion inhibitor: an in silico drug repurposing study

Regarding the urgency of therapeutic measures for coronavirus disease 2019 (COVID-19) pandemic, the use of available drugs with FDA approval is preferred because of the less time and cost required for their development. In silico drug repurposing is an accurate way to speed up the screening of the existing FDA-approved drugs to find a therapeutic option for COVID-19. The similarity in SARS-CoV-2 and HIV-1 fusion mechanism to host cells can be a key point for Inhibit SARS-CoV-2 entry into host cells by HIV fusion inhibitors. Accordingly, in this study, an HIV-1 fusion inhibitor called Enfuvirtide (Enf) was selected. The affinity and essential residues involving in the Enf binding to the S2 protein of SARS-CoV-2, HIV-1 gp41 protein and angiotensin-converting enzyme 2 (ACE-2) as a negative control, was evaluated using molecular docking. Eventually, Enf-S2 and Enf-gp41 protein complexes were simulated by molecular dynamics (MD) in terms of binding affinity and stability. Based on the most important criteria such as docking score, cluster size, energy and dissociation constant, the strongest interaction was observed between Enf with the S2 protein. In addition, MD results confirmed that Enf-S2 protein interaction was remarkably stable and caused the S2 protein residues to undergo the fewest fluctuations. In conclusion, it can be stated that Enf can act as a strong SARS-CoV-2 fusion inhibitor and demonstrates the potential to enter the clinical trial phase of COVID-19.

Communicated by Ramaswamy H. Sarma

Immunoinformatics and molecular dynamics studies to predict T-cell-specific epitopes of four Klebsiella pneumoniae fimbriae antigens

Klebsiella pneumoniae (K. pneumoniae) is a causative agent of severe infections in humans. There is no publically available vaccine for K. pneumoniae infections yet. Here, using comprehensive immunoinformatics methods, T-cell-specific epitopes of four type 1 fimbriae antigens of K. pneumoniae were predicted and evaluated as potential vaccine candidates. Both CD8+ (class I) and CD4+ (class II) T-cell-specific epitopes were predicted and the epitopes similar to human proteome were excluded. Subsequently, the windows of class-II epitopes containing class-I epitopes were determined. The immunogenicity, IFN-γ production and population coverage were also estimated. Using the 3D structure of HLA and epitopes, molecular docking was carried out. Two best epitopes were selected for molecular dynamics studies. Our prediction and analyses resulted in the several dominant epitopes for each antigen. The docking results showed that all selected epitopes can bind to their restricted HLA molecules with high affinity. The molecular dynamics results indicated the stability of system with minimum possible deviation, suggesting the selected epitopes can be promising candidates for stably binding to HLA molecules. Altogether, our results suggest that the selected T-cell-specific epitopes of K. pneumoniae fimbriae antigens, particularly the two epitopes confirmed by molecular dynamics, can be applied for vaccine development. However, the in vitro and in vivo studies are required to authenticate the results of the present study.Communicated by Ramaswamy H. Sarma.

Effects of unsaturated fatty acids (Arachidonic/Oleic Acids) on stability and structural properties of Calprotectin using molecular docking and molecular dynamics simulation approach

Calprotectin is a heterodimeric protein complex with two subunits called S100A8/A9. The protein has an essential role in inflammation process and various human diseases. It has the ability to bind to unsaturated fatty acids including Arachidonic acid, Oleic acid and etc., which could be considered as a major carrier for fatty acids. In this study we aimed to appraise the thermodynamics and structural changes of Calprotectin in presence of Arachidonic acid/Oleic acid) using docking and molecular dynami simulation method. To create the best conformation of Calprotectin-Oleic acid/Arachidonic acid complexes, the docking process was performed. The complexes with the best binding energy were selected as the models for molecular dynamics simulation process. Furthermore, the structural and thermodynamics properties of the complexes were evaluated too. The Root Mean Square Deviation and Root Mean Square Fluctuation results showed that the binding of Arachidonic acid/Oleic acid to Calprotectin can cause the protein structural changes which was confirmed by Define Secondary Structure of Proteins results. Accordingly, the binding free energy results verified that binding of Oleic acid to Calprotectin leads to instability of S100A8/A9 subunits in the protein. Moreover, the electrostatic energy contribution of the complexes (Calprotectin-Oleic acid/Arachidonic acid) was remarkably higher than van der Waals energy. Thus, the outcome of this study confirm that Oleic acid has a stronger interaction with Calprotectin in comparison with Arachidonic acid. Our findings indicated that binding of unsaturated fatty acids to Calprotectin leads to structural changes of the S100A8/A9 subunits which could be beneficial to play a biological role in inflammation process.

Study of HSA interactions with arachidonic acid using spectroscopic methods revealing molecular dynamics of HSA-AA interactions

The interaction between human serum albumin (HSA) and arachidonic acid (AA) as an unsaturated fatty acid were investigated in the present study using methods including UV-VIS spectrophotometry, fluorescence and circular dichroism (CD) spectroscopy, lifetime measurements, fluorescence anisotropy measurements and visual molecular dynamics (MD). The thermodynamic parameters were assessed from HSA thermal and chemical denaturation in the presence and absence of AA. From the thermal denaturation, the Tm and ΔG˚(298K) magnitudes obtained were 327.7 K and 88 kJ/mol, respectively, for HSA alone, and 323.4 K and 85 kJ/mol, respectively, following treatment with a 10 µM AA concentration. The same manner of reduction in Gibbs free energy as a criterion of protein stability was achieved during chemical denaturation by urea in the presence of AA. The present study investigates HSA binding nature through MD approaches, and the results indicated that the binding affinity of AA to the subdomain IIA of HSA is greater compared with that of subdomain IIIA. Although the HSA regular secondary structure evaluation by CD exhibited a minor change following incubation with AA, its tertiary structure revealed an observable fluctuation. Thus, it appears that the interaction between AA and HSA requires minor instability and partial structural changes.

In silico assessment of human Calprotectin subunits (S100A8/A9) in presence of sodium and calcium ions using Molecular Dynamics simulation approach

Calprotectin is a heterodimeric protein complex which consists of two subunits including S100A8 and S100A9. This protein has a major role in different inflammatory disease and various types of cancers. In current study we aimed to evaluate the structural and thermodynamic changes of the subunits and the complex in presence of sodium and calcium ions using molecular dynamics (MD) simulation. Therefore, the residue interaction network (RIN) was visualized in Cytoscape program. In next step, to measure the binding free energy, the potential of mean force (PMF) method was performed. Finally, the molecular mechanics Poisson-Boltzmann surface area (MMPBSA) method was applied as an effective tool to calculate the molecular model affinities. The MD simulation results of the subunits represented their structural changes in presence of Ca2+. Moreover, the RIN and Hydrogen bond analysis demonstrated that cluster interactions between Calprotectin subunits in presence of Ca2+ were greater in comparison with Na+. Our findings indicated that the binding free energy of the subunits in presence of Ca2+ was significantly greater than Na+. The results revealed that Ca2+ has the ability to induce structural changes in subunits in comparison with Na+ which lead to create stronger interactions between. Hence, studying the physical characteristics of the human proteins could be considered as a powerful tool in theranostics and drug design purposes.

Rational MD simulations for improvement the affinity of nanobody against PlGF (placenta growth factor): mutagenesis based on electrostatic interactions

AbbreviationsCOMcenter of mass distanceMDmolecular dynamicsMM-PBSAMolecular Mechanics Poisson-Boltzmann Surface AreaNbnanobodyPlGFplacenta growth factorRgradius of gyrationRMSDroot mean-square deviationSASAsolvent-accessible surface areaVEGFvascular endothelial growth factor.

Identification of HLA-I restricted epitopes in six vaccine candidates of Leishmania tropica using immunoinformatics and molecular dynamics simulation approaches

In spite of numerous studies on vaccination for various species of Leishmania, research on the development of an effective vaccine for L. tropica is very scarce. In silico epitope prediction is a new way to survey the best vaccine candidates. Here, we predicted the best epitopes of six L. tropica antigens with vaccine capability against this pathogen, using highly frequent HLA-I alleles. Based on the frequent HLA alleles, the protein sequences were screened individually using four different MHC prediction applications, namely SYFPEITHI, ProPredI, BIMAS, and IEDB. Several in silico assays including clustering, human similarity exclusion, epitope conservancy prediction, investigating in experimental records, immunogenicity prediction, and prediction of population coverage were performed to narrow the results and to find the best epitopes. The selected epitopes and their restricted HLA-I alleles were docked and the final epitopes with the lowest binding energy (the highest binding affinity) were chosen. Finally, the stability and the binding properties of the best epitope-HLA-I combinations were analyzed using molecular dynamics simulation studies. We found ten potential peptides with strong binding affinity to highly frequent HLA-I alleles that can be further evaluated as vaccine targets against L. tropica.

Mutational pressure and natural selection in epidermal growth factor receptor gene during germline and somatic mutagenesis in cancer cells

In this study we investigated nucleotide usage biases along the length of a gene encoding human epidermal growth factor receptor (EGFR) and found out that there had been mutational GC-pressure with stronger asymmetric C-pressure in that gene before the preferable direction of nucleotide mutations changed. Current preferable direction of germline mutations in EGFR gene has been estimated with the help of Ensembl data base of gene variations. Preferable direction of somatic mutations in EGFR gene from cancer cells has been estimated with the help of COSMIC data base. Both germline and somatic mutations in cancer cells have the same GC to AT preferable direction in EGFR gene. These data have been used with the aim to find fragments of EGFR gene that have lower probability of missense C to T and G to A transitions to occur. So, the less mutable parts of extracellular EGFR domain are: C-terminal part of the first beta barrel and the central part of the second beta barrel. The less mutable parts of tyrosine kinase EGFR domain are: ATP-binding site (partially), regulatory alpha helix, and fragments that change their secondary structure during the activation process. These parts of EGFR should be considered as the best targets for new types of therapy development. Such criterion as low mutability is especially important for the selection of targets for anti-tumor therapy, since we have detected positive selection of amino acid replacements during somatic mutagenesis of EGFR gene in cancer cells.

Recent advances in the selection and identification of antigen-specific nanobodies

Nanobodies represent the next-generation antibody-derived biologics with significant advances over conventional antibodies. Several rapid and robust techniques for isolating highly specific nanobodies have been developed. Antigen specific nanobodies are selected from constructed nanobody libraries, which can be classified into 3 main types: immune library, naïve library, and semisynthetic/synthetic library. The immune library is the most widely used strategy for nanobody screening. Target specific nanobodies are highly enriched in immune libraries than in non-immune libraries; however, it is largely limited by the natural antigenicity of antigens. The naïve library is thus developed. Despite the lack of somatic maturation, protein engineering can be employed to significantly increase the affinities of selected binders. However, a substantial amount of blood samples collected from a large number of individual animals is a prerequisite to ensure the diversity of the naïve library. With this issue considered, the semisynthetic/synthetic library may be a promising path toward obtaining a limitless source of nanobodies against a variety of antigens without the need of animals. In this review, we summarize the state-of-the-art screening technologies with different libraries. The approaches presented here can further boost the diverse applications of nanobodies in biomedicine and biotechnology.

Nanobodies As Novel Agents for Targeting Angiogenesis in Solid Cancers

Solid cancers are dependent on angiogenesis for sustenance. The FDA approval of Bevacizumab in 2004 inspired many scientists to develop more inhibitors of angiogenesis. Although several monoclonal antibodies (mAbs) are being administered to successfully combat various pathologies, the complexity and large size of mAbs seem to narrow the therapeutic applications. To improve the performance of cancer therapeutics, including those blocking tumor angiogenesis, attractive strategies such as miniaturization of the antibodies have been introduced. Nanobodies (Nbs), small single-domain antigen-binding antibody fragments, are becoming promising therapeutic and diagnostic proteins in oncology due to their favorable unique structural and functional properties. This review focuses on the potential and state of the art of Nbs to inhibit the angiogenic process for therapy and the use of labeled Nbs for non-invasive in vivo imaging of the tumors.