Drug targeting is a progressive area of research with folate receptor alpha (FRα) receiving significant attention as a biological marker in cancer drug delivery. The binding affinity of folic acid (FA) to the FRα active site provides a basis for recognition of FRα. In this study, FA was conjugated to beta-cyclodextrin (βCD) and subjected to in silico analysis (molecular docking and molecular dynamics (MD) simulation (100 ns)) to investigate the affinity and stability for the conjugated system compared to unconjugated and apo systems (ligand free). Docking studies revealed that the conjugated FA bound into the active site of FRα with a docking score (free binding energy < -15 kcal/mol), with a similar binding pose to that of unconjugated FA. Subsequent analyses from molecular dynamics (MD) simulations, root mean square deviation (RMSD), root mean square fluctuation (RMSF), and radius of gyration (Rg) demonstrated that FA and FA-βCDs created more dynamically stable systems with FRα than the apo-FRα system. All systems reached equilibrium with stable RMSD values ranging from 1.9-2.4 Å and the average residual fluctuation values of the FRα backbone atoms for all residues (except for terminal residues ARG8, THR9, THR214, and LEU215) were less than 2.1 Å with a consistent Rg value of around 16.8 Å throughout the MD simulation time (0-100 ns). The conjugation with βCD improved the stability and decreased the mobility of all the residues (except residues 149-151) compared to FA-FRα and apo-FRα systems. Further analysis of H-bonds, binding free energy (MM-PBSA), and per residue decomposition energy revealed that besides APS81, residues HIS20, TRP102, HIS135, TRP138, TRP140, and TRP171 were shown to have more favourable energy contributions in the holo systems than in the apo-FRα system, and these residues might have a direct role in increasing the stability of holo systems.

Folate receptor alpha (FRα) is known as a biological marker for many cancers due to its overexpression in cancerous epithelial tissue. The folic acid (FA) binding affinity to the FRα active site provides a basis for designing more specific targets for FRα. Heterocyclic rings have been shown to interact with many receptors and are important to the metabolism and biological processes within the body. Nineteen FA analogs with substitution with various heterocyclic rings were designed to have higher affinity toward FRα. Molecular docking was used to study the binding affinity of designed analogs compared to FA, methotrexate (MTX), and pemetrexed (PTX). Out of 19 FA analogs, analogs with a tetrazole ring (FOL03) and benzothiophene ring (FOL08) showed the most negative binding energy and were able to interact with ASP81 and SER174 through hydrogen bonds and hydrophobic interactions with amino acids of the active site. Hence, 100 ns molecular dynamics (MD) simulations were carried out for FOL03, FOL08 compared to FA, MTX, and PTX. The root mean square deviation (RMSD) and root mean square fluctuation (RMSF) of FOL03 and FOL08 showed an apparent convergence similar to that of FA, and both of them entered the binding pocket (active site) from the pteridine part, while the glutamic part was stuck at the FRα pocket entrance during the MD simulations. Molecular mechanics Poisson-Boltzmann surface accessible (MM-PBSA) and H-bond analysis revealed that FOL03 and FOL08 created more negative free binding and electrostatic energy compared to FA and PTX, and both formed stronger H-bond interactions with ASP81 than FA with excellent H-bond profiles that led them to become bound tightly in the pocket. In addition, pocket volume calculations showed that the volumes of active site for FOL03 and FOL08 inside the FRα pocket were smaller than the FA-FRα system, indicating strong interactions between the protein active site residues with these new FA analogs compared to FA during the MD simulations.