Abstract
L-amino acids play a vital role in various biological processes, including hormone production, adipose tissue growth and repair, and protein synthesis. Consequently, the design of novel L-amino acid-based drugs with fewer side effects has recently garnered significant attention. Ketorolac, while effective for pain reduction, presents several adverse side effects, motivating the search for alternatives. An ideal alternative should structurally resemble ketorolac but offer improved chemical and thermal stability, as well as enhanced hydrogen bonding capability. In this project, we designed an L-amino acid-based structure as a potential alternative to ketorolac. Using computer-aided drug discovery methods, we generated a detailed, atomic-level description of the drug ligand and evaluated its selectivity properties, ADMET prediction, and Pre-Tox analysis. Molecular docking calculations were then performed to identify the most favorable active site for the newly designed ligand, mimicking ketorolac's binding pose within the protein. The resulting AVH structure exhibited a close resemblance to ketorolac. After these initial assessments, we performed molecular dynamics simulations. The RMSD and RMSF plots indicate the stability of the dynamic structure of the substituted ligand. Analysis of the rotation and hydrogen bond plots revealed strong interactions between the protein and the ligand. Furthermore, analysis of the protein secondary structure and residue-wise secondary structure using the DSSP algorithm confirmed the suitability of the proposed structure and, in some cases, showed improved stability compared to its interaction with ketorolac. In the DCCM analysis, a stronger correlation pattern was observed in AVH, indicating greater stability in this pattern. Principal component analysis based on a free energy perspective revealed closely spaced data points representing low-energy and relatively stable states that are energetically favorable for the system. The protein binding cavity in the electrostatic potential map shows negative and neutral regions, likely indicating greater stabilization of the ligand in the active site. Finally, the MM-PBSA results suggest that although the binding energy of AVH is slightly higher than that of KTR, it remains a negative value, which is thermodynamically favorable, and the binding affinity of AVH to the protein remains strong. Based on these results, we selected and synthesized the best structure compared to ketorolac. We predict that this new structure, designed using a natural product and showing similar chemical effects, will be a potential alternative to ketorolac.