Knot polynomials have been used to detect and classify knots in biomolecules. Computation of knot polynomials in DNA and protein molecules have revealed the existence of knotted structures, and provided important insight into their topological structures. However, conventional knot polynomials are not well suited to study RNA molecules, as RNA structures are determined by stem regions which are not taken into account in conventional knot polynomials. In this study, we develop a new class of knot polynomials specifically designed to study RNA molecules, which considers stem regions. We demonstrate that our knot polynomials have direct structural relation with RNA molecules, and can be used to classify the topology of RNA secondary structures. Furthermore, we point out that these knot polynomials can be used to model the topological effects of disulfide bonds in protein molecules.