In this article we show how the quasi-Gaussian entropy (QGE) theory can be used to treat systems of polyatomic flexible molecules, where the usual semirigid description is not always appropriate. We describe a completely general derivation of the QGE theory which does not make use of any semirigid approximation, and therefore it is very suited for large and flexible molecules. Using molecular dynamics simulations of flexible molecules in vacuo, we investigated the ability of the theory to describe intramolecular energy fluctuations and conformational equilibria of purely classical molecules in the ideal gas condition. Results show that the gamma state level of the theory and its generalization for treating conformational equilibria (multi-gamma state model) provide excellent theoretical models when applied to three polyatomic molecules of increasing conformational freedom.