I am trying to use stochastic modelling in python to model the changes in concentrations of RNA monomers A, B and their dimers AB, and I need to find the rate constants of A+B<->AB. Are there any functions in the ViennaRNA library that can output the intermediate state of the process of two RNA molecules cofolding, thus letting us find the activation energy and find the rate constants of the two RNA molecules folding into a duplex and the reverse reaction? I tried to look for such functions but was not successful in doing so. If such function(s) do not exist in ViennaRNA, are there any other ways to find the rate of RNA cofolding? Also I would like to ask if there is a full python library for ViennaRNA as some of the C functions do not seem to have python counterparts in the library provided by the ViennaRNA package, or are there ways to call C functions using python?