Molecular Dynamics and Free Energy Perturbation Analysis of Guava (Psidium guajava) Leaf-Derived Inhibitors Against Coagulation Factor IXa

Joseph De Luna,Shanahi Chelledie Gonzales,Jimuel Jan Nuqui,Evalyn Capinding,Corazon Sacdalan
DOI: https://doi.org/10.26434/chemrxiv-2024-zrhp9
2024-12-03
Abstract:There is an ambiguity whether guava (Psidium guajava) leaves support or inhibits coagulation as it was commonly used as wound dressing. Our previous study extracted, identified, and screened using docking, bioactive compounds from guava leaves that has the potential to be coagulation factor IXa inhibitor. Docking has static limitations that require further evaluation [1]. To investigate the potential of bioactive compounds as inhibitors of coagulation factor IXa (CFIXa) using molecular dynamics (MD) simulations and binding free energy perturbation analyses.Top three docking-scored CFIXa-ligand complexes from our previous study were evaluated using MD simulations with GROMACS 2024 and CHARMM force fields. Simulations included minimization, NVT, NPT with and without restraints, followed by a 10 ns production run. Structural analysis was performed using RMSD, RMSF, SASA, and secondary structure analysis in UCSF ChimeraX. Binding free energies (ΔGbind) were calculated using gmx_mmpbsa with PB and GB models, and decomposition analysis identified key residues contributing to binding. Compound 7 showed the strongest binding with CFIXa, with a binding free energy ranging from -10 to -20 kcal/mol and stable interactions, particularly at H:THR245, while compound 2 demonstrated moderate binding with a binding free energy of -5 to -15 kcal/mol. Compound 3 exhibited instability, with RMSD fluctuations from 0.06 to 7.71 nm, and binding free energy approaching 0 kcal/mol, confirming its unbinding during the simulation. Compound 7 demonstrated the strongest binding interactions and stability, making it the most promising candidate for further development, while compound 2 showed moderate binding and potential for optimization through structural refinements. In contrast, compound 3 exhibited instability and weak binding, suggesting that modifications could improve its binding efficiency and specificity.
Chemistry
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to evaluate whether the bioactive compounds extracted from guava (Psidium guajava) leaves can be effective inhibitors of coagulation factor IXa (CFIXa). Specifically, the research aims to further verify the dynamic interactions and stability between these compounds and CFIXa through molecular dynamics (MD) simulations and binding free - energy perturbation analysis, thereby providing a theoretical basis for the development of new anticoagulant drugs. ### Research Background 1. **Coagulation Disorders**: Thrombotic diseases are an important global health challenge, characterized by the abnormal formation of blood clots within blood vessels, which may lead to serious complications such as heart attacks and strokes. Coagulation factor IXa plays a key role in the endogenous coagulation pathway, and its inhibition helps to reduce thrombosis while reducing the risk of severe bleeding caused by direct thrombin inhibitors. 2. **Potential of Natural Compounds**: Natural compounds of plant origin have received extensive attention for their potential therapeutic applications. Guava leaves have multiple pharmacological properties, such as anti - inflammatory, antioxidant, and antibacterial effects, but their role in the coagulation process is not clear. Some studies suggest that they may have pro - coagulant activity, while others point out their potential anticoagulant effects, so more in - depth research is needed to understand their specific impact on hemostasis. ### Research Methods - **Preliminary Research**: The authors previously screened out several possible CFIXa inhibitors through docking studies and discovered some preliminary binding site information. However, docking studies are limited to static structural analysis and cannot capture the dynamic behavior of protein - ligand complexes. - **Molecular Dynamics Simulation**: To overcome this limitation, this study used GROMACS to conduct MD simulations to explore the dynamic properties, flexibility, and stability of these compounds under near - physiological conditions. - **Binding Free - Energy Calculation**: The MM - PB/GBSA method was used to calculate the binding free - energy to quantitatively evaluate the thermodynamic stability of the complex, and the contribution of key residues to binding was determined through decomposition analysis. ### Main Conclusions - **Compound 7** exhibits the strongest binding ability and stability, especially with a stable interaction at H:THR245, and the binding free - energy is between - 10 and - 20 kcal/mol. This makes it the most promising candidate drug. - **Compound 2** shows a moderate binding ability, with a binding free - energy between - 5 and - 15 kcal/mol, indicating that its specificity and binding strength can be improved through structural optimization. - **Compound 3** shows significant instability, with large RMSD fluctuations and a binding free - energy close to 0 kcal/mol, and finally dissociates during the simulation process. Nevertheless, its initial interaction with catalytic residues suggests that its binding efficiency and specificity can be enhanced through improvement. In conclusion, this study not only provides a more comprehensive understanding of the CFIXa inhibition mechanism but also lays the foundation for the development of new anticoagulant drugs based on guava leaf extracts.