A subscription to JoVE is required to view this content. Sign in or start your free trial.

Summary

This antibody homology modeling prediction protocol is followed by antibody-receptor Pyrx docking and molecular dynamic simulation. These three primary methods are used to visualize the accurate antibody-receptor binding areas and the binding stability of the final structure.

Abstract

Single-chain fragment variable (scFv) antibodies were previously constructed of variable light and heavy chains joined by a (Gly4-Ser) 3 linker. The linker was created using molecular modeling software as a loop structure. Here, we introduce a protocol forin silico analysis of a complete scFv antibody that interacts with the epidermal growth factor receptor (EGFR). The homology modeling, with Pyrx of protein-protein docking and molecular dynamic simulation of the interacting scFv antibody and EGFR First, the authors used a protein structure modeling program and Python for homology modeling, and the antibody scFv structure was modeled for homology. The investigators downloaded Pyrx software as a platform in the docking study. The Molecular dynamic simulation was run using modeling software. Results show that when the MD simulation was subjected to energy minimization, the protein model had the lowest binding energy (-5.4 kcal/M). In addition, the MD simulation in this study showed that the docked EGFR-scFv antibody was stable for 20-75 ns when the movement of the structure increased sharply to 7.2 Å. In conclusion, in silicoanalysiswas performed, and the molecular docking and molecular dynamics simulations of the scFv antibody proved the effectiveness of the designed immune-therapeutic drug scFv as a specific drug therapy for EGFR.

Introduction

Conformational changes in the protein (ligand and receptor) always occur based on structure-based functions. The study of the possible binding grooves of the protein and prediction of the stable binding interaction is an advanced method to prepare drugs for better use in the human body. Homology modeling followed by docking and molecular dynamic simulation is a straightforward method for accurate prediction of stable interactions of binding between the residues of receptors and constructed antibodies that are used as specific personalized medicine1,2. The predicted model structure can show conformational changes and rearrangements in ligand-receptor binding sites, particularly at the antibody-receptor interface. There are many reasons for these changes, such as the rotation of side chains, global structural transformation, or more complex modifications. The main reason for homology modeling is to distinguish a protein's tertiary structure from its primary structure2,3.

A tyrosine kinase receptor called epidermal growth factor receptor (EGFR) plays many biological roles in cancer cells, including apoptosis4,5, differentiation6,7, cell cycle progression8,9, development9,10, and transcription11. EGFR is one of the well-known therapeutic targets for breast cancer12. The overexpression of regular kinase activity such as EGFR usually leads to cancer cell progression, which can be repressed by many kinds of cancer inhibitors13. The epidermal growth factor receptor (EGFR) was used as a receptor for the single chain fragment variable specifically constructed to work against this receptor. Its predicted structure was used to test the antibody binding activity.

In this paper, the scFv antibody structure was modeled using modeling software with Python script and the homology modeling method14,15. A homology model can be built from the protein and amino acid sequences of receptors and ligands16,17. Additionally, advanced bioinformatics technologies such as molecular docking were employed to predict how small molecule ligands will bind to the correct target binding site. The docking would balance the development of novel drugs directed toward multiple diseases. The binding behavior is taken into consideration5,18.

Furthermore, molecular docking is a critical technique to facilitate and speed up ligand-receptor binding development. Molecular docking enables scientists to virtually screen a library of ligands against a target protein and predict the binding conformations and affinities of the ligands to the target receptor protein. Molecular dynamic simulation (MNS) demonstrates how the residues move in space, simulates the antibody motions toward their receptors, and finally informs antibody design efforts. This study is a novel prediction of grid box dimensions that decided how the scFv antibody binds to EGFR and the detection of the energy and time of that binding in MDsimulation.

Protocol

1. Secondary structure predictions of a single chain fragment variable (scFv) protein

  1. Build the single-chain fragment variable (scFv) protein's 3D structure with BLAST protein data bank (PDB), KABAT numbering, and the modeling software. The scFv consists of a linker (Gly4-Ser) that connects a variable heavy chain (VH) and a variable light chain (VL).
  2. Use the molecular modeling software to build the linker as a loop structure, and perform all these methods as described in previous studies2,19,20.

2. Template selection and scFv and EGFR 3D structure prediction and homology modeling

  1. Choose template 1ivo for EGFR structures (based on its high resolution). Download the 1ivo.pdb file from the pdb website, as shown in Figure 1B.
  2. Prepare the input 1ivo.pdb file as described below.
    1. In the 1ivo.pdb file, remove all external ligands by opening the pdb.org website and selecting the 1ivo. Structure, and looking for the name of the ligands under the small molecule title on the 1ivo structure page of the pdb website.
    2. Find the ligand name NAG. Open the 1ivo.pdb file downloaded from the pdb website and find the termination residue (TER.).
    3. Delete the residues of the external ligands in the 1ivo structure, starting from the residue after TER. and before the residue ends. Save the 1ivo.pdb file on the system.
  3. Prepare the saved 1ivo.pdb file as described below.
    1. Download the Autodock docking software (autodock.scripps.edu) from the window selection area. Click on the Open 1ivo.pdb File.
    2. Use the Edit command to choose Add Hydrogen > Add, then select Polar Only, and then press Ok.
    3. Use the Edit command to add Kollman charges (Supplementary Figure 1). Use the Edit command to delete water. Save the 1ivo.pdb file on the pc.
  4. Minimize the energy of the 1ivo.pdb structure as described below.
    1. Download SPDBV. software from http://spdbv.vital-it.ch/disclaim.html. Open the 1ivo.pdb file.
    2. Select all. Select the command Perf > Energy minimization > Ok (Supplementary Figure 2). Save the 1ivo file on the pc.
  5. Prepare the full model scFv using homology modeling as described below.
    1. Download the modeling software17 and the Python script 3.7.9 shell from the Window- 64. Keep the downloaded software files in the D drive.
  6. Prepare the input files as described below.
    1. Load the scFv Pdb file in fasta format from the NCBI website and rename the file TARGET.ali. as described in Supplementary Coding File 1. Choose the template using the Blast section in NCBI, paste the sequenced file, select in pdb format 7det.pdb as described in Supplementary Coding File 2, and then submit. Then, use the pdb.org website to obtain the template file.
    2. Prepare the third input file as align2d.py (Python) as described in Supplementary Coding File 3, which opens as shown in Supplementary Figure 3A. Press the Show More Option, then go to Edit with IDLE > Edit with EDLE (64-bit). Run using the run module 5 command in the align2d.py to obtain two output files: Tar- 7det.ali and Tar- 7det.pap.
  7. Complete the previous three steps to use the command in the last input file.
  8. Add a new input file model-single.py (command python file) as shown in Supplementary Coding File 4 and as described below.
    1. Press the Show More Option, then go to Edit with IDLE > Edit with EDLE3.7 (64-bit). Run using the (run module 5) command as shown in Supplementary Figure 3B.
      ​NOTE: The resulting output files are the six files of the homology models shown in Supplementary Figure 3C.

3. Receptor secondary structure prediction and evaluation

  1. Detect the homology models' correction and accuracy as described below.
    1. Create the Ramachandran plot for the scFv models and EGFR models by downloading the visualization tool from  https://discover.3ds.com/discovery-studio-visualizer-download.
    2. Open the file, then right-click with the mouse and select the display sequence (Supplementary Figure 4). Copy the sequence and paste it into the Pictorial database of 3D structures (pdbsum) www.ebi.ac.uk/thornton-srv/databases/pdbsum/.
    3. Select search by sequence, paste the copy of the sequence, and then submit it. Create the plot as shown in Figure 1B,D.

4. Protein-protein docking

  1. Download the virtual screening tool software.
  2. Go to File> Read Molecules > Load 1ivo.pdb. Right-click on the protein in the autodock panel to make a macromolecule. Right-click again to make a ligand (Supplementary Figure 5).
  3. Click in the autodock panel to select protein and then select ligand.
  4. Open the protein list. Then, from the list, select scFv Protein.
  5. Go to Toggle Selection Spheres. Adjust the grid box to the center of the receptor. Click Forward when the round pink button appears.
  6. To prepare the pdbqt files for both scFv-antibody and EGFR (1ivo) structures, use the following steps.
    1. Go to C Drive > Program Files (86) > Users, then choose the pyrx file that contains macromolecules and protein output files that were saved as a pdbqt file.
    2. Then, save the single-chain fragment variable (scFv) antibody pdbqt file.
  7. Download PyMOL software at PyMOL | pymol.org. Use PyMOL software to show the scFv antibody-receptor EGFR configurations.
    1. Go to file and open C:\Users\ilham\.mgltools\PyRx\Macromolecules\protein. Prepare the docking configurations of the scFvantibodyinteracting with the receptor in Figure 2A as described below.
      1. Use the display option to show the (1ivo)-receptor file as sequence residues with a white background shown in Supplementary Figure 6.
      2. Display the docking configuration file with higher resolution to see the ligand color in green and red residue colors. Display the (1ivo)-receptor rigid surface in yellow.
  8. Prepare the docking configurations of the scFvantibodyinteracting with the receptor in Figure 2B, as described below.
    1. Download the docking software from the window selection area. Use Autodock to show the scFv antibody-receptor EGFR configurations and conformations.
      1. In Autodock, choose the Analyze Option, then open the Autodock Vina result. Go to File and open C:\Users\ilham\.mgltools\PyRx\Macromolecules\protein.
      2. Select the protein receptor pdb file, then select the area of the ligand configuration (scFv antibody structure). Connect the rigid surface of the receptor with the configuration ofthe docked structure and hide the rest of the receptor. Hide the far residues of the receptor from the connected residues to the ligand, as shown in Figure 2B.
    2. The protein-protein complex was then considered ready to perform MD simulation.

5. Molecular dynamic simulation (MD simulation) of the EGFR-scFv antibody docking complex

  1. Download the MD simulation software and use it as follows.
    1. Prepare the EGFR (1ivo) pdb file using the reparation wizard as in Supplementary Figure 7A. Operate the preprocess section to refine the file. Send the refined file to be set in the system builder.
    2. Load the molecular dynamics simulation software from the working directory. Add the ions and upgrade the refined file to reach 20 Å to submit the job (Table 1), also shown in Supplementary Figure 7B.
    3. Load the EGFR (1ivo) pdb from the imported file, then choose 100 ns timesteps to run it (Supplementary Figure 7C).
  2. Start the analysis of the simulation after the completion of the MDsimulation as described below.
    1. Create a job folder and save it in the cms file . Load the cms file to perform this step in the MD simulation.
    2. Create a working directory for project folders and report the energy values. Use S.I.D. pdf to report the simulation, as shown in Figure 3A, and the interaction diagram and H bond, as shown in Figure 3B.
    3. Load the pdf file of cms by browsing from the folder and use TIP3P as the model for file volume minimization.
    4. Create the solvation file to perform this step shown in Supplementary Figure 7D. Save pdf file through the software, and analyze the data, resulting in Figure 4, Figure 5, and Figure 6.
  3. Generate an MDsimulation finalization setup by creating the resolve file. Find the results in the boundary box, as shown in Figure 7.

Results

Using phage display technology, the scFv gene anti-EGFR was created from the mouse B-cell hybridoma line C3A820,21. The single chain fragment variable (scFv) structure models of the VH and VL structures were built separately, according to Chua et al.22. After that, the models were visible as ribbons produced using RasMol. Then, using molecular modeling software, a synthetic peptide [Gly4Ser)3 was used to join the separately modeled VH and ...

Discussion

EGFR is the primary target receptor of breast cancer. EGFR overexpression increases breast cancer cases around the world. Meanwhile, specific antibodies such as single chain fragment variables are antibodies that move easily via blood circulation and have a fast clearance rate in the body. Antibodies are a wise solution and an effective immunotherapy drug37. Therefore, structure-based drug design must identify inhibitory medicines, such as scFv antibodies, that work specifically against a target r...

Disclosures

The authors have nothing to disclose.

Acknowledgements

None.

Materials

NameCompanyCatalog NumberComments
Autodock softwareCenter for Computational structural Biology AutoDock (scripps.edu)
Desmond Maestro 19.4 software Schrodingerwww.schrodinger.com 
Download Discovery Studio 2021  Dassault Systems https://discover.3ds.com/discovery-studio-visualizer-download.
Modeler Version 9.24[17] University of Californiahttps://salilab.org/modeller/9.24/release.html
Pictorial database of 3D structures (pdbsum)EMBL-EBI www.ebi.ac.uk/thornton-srv/databases/pdbsum/
PyMOL software SchrodingerPyMOL | pymol.org
Pyrx software Sourceforge Download PyRx - Virtual Screening Tool (sourceforge.net)
Python script 3.7.9 shell from the window (64)PythonPython Release Python 3.7.9 | Python.org
SPDBV software Expasyhttp://spdbv.vital-it.ch/disclaim.html

References

  1. Clark, J. J., Orban, Z. J., Carlson, H. A. Predicting binding sites from unbound versus bound protein structures. Sci Rep. 10 (1), 15856 (2020).
  2. Huang, Y., et al. A stepwise docking molecular dynamics approach for simulating antibody recognition with substantial conformational changes. Comput Struct Biotechnol J. 20, 710-720 (2022).
  3. Mahgoub, I. O., Ali, A. M., Hamid, M., Alitheen, N. M. Single chain fragment variable (scFv) secondary structure prediction and evaluation. FASEB J. 25, (2011).
  4. Kim, H., et al. Titanium dioxide nanoparticles induce apoptosis by interfering with EGFR signaling in human breast cancer cells. Environ Res. 175, 117-123 (2019).
  5. Mahgoub, E. O., Abdella, G. M. Improved exosome isolation methods from non-small lung cancer cells (NC1975) and their characterization using morphological and surface protein biomarker methods. J Cancer Res Clin Oncol. 149 (10), 7505-7514 (2023).
  6. Shaurova, T., Zhang, L., Goodrich, D. W., Hershberger, P. A. Understanding lineage plasticity as a path to targeted therapy failure in EGFR-mutant non-small cell lung cancer. Front Genet. 11, 281 (2020).
  7. Mahgoub, E. O., Bolad, A. K. Construction, expression and characterisation of a single chain variable fragment in the Escherichia coli periplasmic that recognise MCF-7 breast cancer cell line. J Cancer Res Ther. 10 (2), 265-273 (2014).
  8. Mahgoub, I. O. Design, expression and characterization of a single chain fragment variable anti-mcf-7 antibody; A humanized antibody derived from monoclonal antibody. Ann Res Conf Proc. 2014 (1), (2014).
  9. Sun, X. L., et al. Dimeric (-)-epigallocatechin-3-gallate inhibits the proliferation of lung cancer cells by inhibiting the EGFR signaling pathway. Chem Biol Interact. 365, 110084 (2022).
  10. Mahgoub, E. O., Haik, Y., Qadri, S. Comparison study of exosomes molecules driven from (NCI1975) NSCLC cell culture supernatant isolation and characterization techniques. FASEB J. 33, 647 (2019).
  11. Madeddu, C., et al. EGFR-mutated non-small cell lung cancer and resistance to immunotherapy: Role of the tumor microenvironment. Int J Mol Sci. 23 (12), 6489 (2022).
  12. Mahgoub, I., Bolad, A. K., Mergani, M. and immune-characterization of single chain fragment variable (scFv) antibody recognize breast cancer cells line (MCF-7). J. Immunother Cancer. 2, 6 (2014).
  13. Ladner, R. C., Sato, A. K., Gorzelany, J., de Souza, M. Phage display-derived peptides as therapeutic alternatives to antibodies. Drug Discov Today. 9 (12), 525-529 (2004).
  14. Khan, M. H., Manoj, K., Pramod, S. Reproductive disorders in dairy cattle under semi-intensive system of rearing in North-Eastern India. Vet World. 9 (5), 512-518 (2016).
  15. Al-Refaei, M. A., Makki, R. M., Ali, H. M. Structure prediction of transferrin receptor protein 1 (TfR1) by homology modelling, docking, and molecular dynamics simulation studies. Heliyon. 6 (1), 03221 (2020).
  16. Anbuselvam, M., et al. Structure-based virtual screening, pharmacokinetic prediction, molecular dynamics studies for the identification of novel EGFR inhibitors in breast cancer. J Biomol Struct Dyn. 39 (12), 4462-4471 (2021).
  17. Oduselu, G. O., Ajani, O. O., Ajamma, Y. U., Brors, B., Adebiyi, E. Homology modelling and molecular docking studies of selected substituted benzo d imidazol-1-yl) methyl) benzimidamide scaffolds on Plasmodium falciparum adenylosuccinate lyase receptor. Bioinform Biol Insights. 13, 1177932219865533 (2019).
  18. Mahgoub, E. O., et al. et al. of exosome isolation techniques in lung cancer. Mol Biol Rep. 47 (9), 7229-7251 (2020).
  19. Shree, P. Targeting COVID-19 (SARS-CoV-2) main protease through active phytochemicals of ayurvedic medicinal plants -Withania somnifera (Ashwagandha), Tinospora cordifolia (Giloy) and Ocimum sanctum (Tulsi) - a molecular docking study. J Biomol Struct Dyn. 40 (1), 190-203 (2022).
  20. Mahgoub, I. O. Expression and characterization of a functional single-chain variable fragment (scFv) protein recognizing MCF7 breast cancer cells in E. coli cytoplasm. Biochem Genet. 50 (7-8), 625-641 (2012).
  21. Mahgoub, E. O. Single chain fragment variables antibody binding to EGF receptor in the surface of MCF7 breast cancer cell line: Application and production review. Open J Genet. 7 (2), 84-103 (2017).
  22. Heng, C. K., Othman, R. Y. Bioinformatics in molecular immunology laboratories demonstrated: Modeling an anti-CMV scFv antibody. Bioinformation. 1 (4), 118-120 (2006).
  23. Mahgoub, E. O., Ahmed, B. Correctness and accuracy of template-based modeled single chain fragment variable (scFv) protein anti- breast cancer cell line (MCF-7). Open J Genet. 3 (3), 183-194 (2013).
  24. Khare, N., et al. Homology modelling, molecular docking and molecular dynamics simulation studies of CALMH1 against secondary metabolites of Bauhinia variegata to treat Alzheimer's disease. Brain Sci. 12 (6), 770 (2022).
  25. Hu, J., Rao, L., Fan, X., Zhang, G. Identification of ligand-binding residues using protein sequence profile alignment and query-specific support vector machine model. Anal Biochem. 604, 113799 (2020).
  26. Santos, L. H. S., Ferreira, R. S., Caffarena, E. R. Integrating molecular docking and molecular dynamics simulations. Methods Mol Biol. 2053, 13-34 (2019).
  27. Xu, D., Tsai, C. J., Nussinov, R. Hydrogen bonds and salt bridges across protein-protein interfaces. Protein Eng. 10 (9), 999-1012 (1997).
  28. Das, S., Chakrabarti, S. Classification and prediction of protein-protein interaction interface using machine learning algorithm. Sci Rep. 11 (1), 1761 (2021).
  29. Mahmoud, S. S. A., Elkaeed, E. B., Alsfouk, A. A., Abdelhafez, E. M. N. Molecular docking and dynamic simulation revealed the potential inhibitory activity of opioid compounds targeting the main protease of SARS-CoV-2. Biomed Res Int. 2022, 1672031 (2022).
  30. Villada, C., Ding, W., Bonk, A., Bauer, T. Simulation-assisted determination of the minimum melting temperature composition of MgCl2-KCl-NaCl salt mixture for next-generation molten salt thermal energy storage. Front. Energy Res. 10, 809663 (2022).
  31. Hog, S. E., Rjiba, A., Jelassi, J., Dorbez-Sridi, R. NaCl salt effect on water structure: a Monte Carlo simulation study. Phys Chem Liq. 60 (5), 767-777 (2022).
  32. Maruyama, Y., Igarashi, R., Ushiku, Y., Mitsutake, A. Analysis of protein folding simulation with moving root mean square deviation. J Chem Inf Model. 63 (5), 1529-1541 (2023).
  33. Winarski, K. L., et al. Vaccine-elicited antibody that neutralizes H5N1 influenza and variants binds the receptor site and polymorphic sites. Proc Natl Acad Sci U S A. 112 (30), 9346-9351 (2015).
  34. Awal, M. A., et al. Structural-guided identification of small molecule inhibitor of UHRF1 methyltransferase activity. Front Genet. 13, 928884 (2022).
  35. Sharma, P., Hu-Lieskovan, S., Wargo, J. A., Ribas, A. adaptive, and acquired resistance to cancer immunotherapy. Cell. 168 (4), 707-723 (2017).
  36. Geerds, C., Wohlmann, J., Haas, A., Niemann, H. H. Structure of Rhodococcus equi virulence-associated protein B (VapB) reveals an eight-stranded antiparallel β-barrel consisting of two Greek-key motifs. Acta Crystallogr F Struct Biol Commun. 70 (7), 866-871 (2014).
  37. Mahgoub, E., et al. The therapeutic effects of tumor treating fields on cancer and noncancerous cells. Arabian J Chem. 14 (10), 103386 (2021).
  38. Hospital, A., Goni, J. R., Orozco, M., Gelpi, J. L. Molecular dynamics simulations: advances and applications. Adv Appl Bioinform Chem. 8, 37-47 (2015).
  39. Kaur, J., Tiwari, R., Kumar, A., Singh, N. Bioinformatic analysis of Leishmania donovani long-chain fatty acid-CoA ligase as a novel drug target. Mol Biol Int. 2011, 278051 (2011).
  40. Páll, S., Hess, B. A flexible algorithm for calculating pair interactions on SIMD architectures. Comp Phy Comm. 184 (12), 2641-2650 (2013).
  41. Tobi, D., Bahar, I. Structural changes involved in protein binding correlate with intrinsic motions of proteins in the unbound state. Proc Natl Acad Sci U S A. 102 (52), 18908-18913 (2005).
  42. Sinha, S., Wang, S. M. of VUS and unclassified variants in BRCA1 BRCT repeats by molecular dynamics simulation. Comput Struct Biotechnol J. 18, 723-736 (2020).
  43. Patel, D., Athar, M., Jha, P. C. Exploring Ruthenium-based organometallic inhibitors against Plasmodium falciparum calcium dependent kinase 2 (PfCDPK2): A combined ensemble docking, QM/MM and molecular dynamics study. Chem Select. 6 (32), 8189-8199 (2021).
  44. Khan, A., et al. Higher infectivity of the SARS-CoV-2 new variants is associated with K417N/T, E484K, and N501Y mutants: an insight from structural data. J Cell Physiol. 236 (10), 7047-7057 (2021).
  45. Zhang, Q., Shao, d., Xu, P., Jiang, Z. Effects of an electric field on the conformational transition of the protein: Pulsed and oscillating electric fields with different frequencies. Polymers. 14 (1), 123 (2021).
  46. Mahtarin, R., et al. et al. and dynamics of membrane protein in SARS-CoV-2. Biomol Struct Dyn. 40 (10), 4725-4738 (2022).

Reprints and Permissions

Request permission to reuse the text or figures of this JoVE article

Request Permission

This article has been published

Video Coming Soon

We use cookies to enhance your experience on our website.

By continuing to use our website or clicking “Continue”, you are agreeing to accept our cookies.

Learn More