Small-molecule RNA binders have emerged as an important pharmacological modality. A profound understanding of the ligand selectivity, binding mode, and influential factors governing ligand engagement with RNA targets is the foundation for rational ligand design. Here, we report a novel class of coumarin derivatives exhibiting selective binding affinity towards single G RNA bulges. Harnessing the computational power of all-atom Gaussian accelerated Molecular Dynamics (GaMD) simulations, we unveiled a rare minor groove binding mode of the ligand with a key interaction between the coumarin moiety and the G bulge. This predicted binding mode is consistent with results obtained from structure-activity-relationship (SAR) studies and transverse relaxation measurements by NMR spectroscopy. We further generated 444 molecular descriptors from 69 coumarin derivatives and identified key contributors to the binding events, such as charge state and planarity, by lasso (least absolute shrinkage and selection operator) regression. Strikingly, small structure perturbations on these key contributors, such as the addition of a methyl group that disrupts the planarity of the ligand resulted in > 100-fold reduction in the binding affinity. Our work deepened the understanding of RNA-small molecule interactions and integrated a new generalizable platform for the rational design of selective small-molecule RNA binders.