In silico study of the binding of daunomycin and phenylalanine transfer RNA: probe molecular recognition for structure-based drug design

Rational designs of pharmaceutical compounds targeting specific RNAs require a comprehensive understanding of molecular recognition mechanisms. Knowledge of binding affinity and specificity can be gained via computational modeling and simulation techniques. In this work, an integrated computational strategy combining QM calculation, molecular docking, conventional and adaptive steered MD simulations, and the var-MM/GBSA approach was proposed to probe the binding behaviors of daunomycin (DAU) and phenylalanine transfer RNA (tRNA(Phe)) at a micro-scale level. Gathering experimental information enables us to eliminate improper predictions for the binding of DAU and tRNA(Phe), and the calculations of PMF and Delta G(binding) lead to the identification of the binding structure of the complex. Further, structural and energetic analysis of the DAU:tRNA(Phe) complex revealed that daunomycinone of DAU contributes the intermolecular VDW energies to nucleotides G(15), C-48 and U-59 on tRNA(Phe), responsible for the binding specificity; meanwhile daunosamine contributes the intermolecular ELE + EGB energies to U-50, responsible for the binding affinity.