A new approach for the theoretical exploration of the potential energy surface and thermodynamic properties evaluation of structurally-nonrigid molecular complexes explicitly accounting the rovibrational dynamics of monomers within the complex has been developed and applied to the estimation of thermodynamic properties of the hydrogen bonded complexes (H2O)2, H2OO3 and H2OSO2. The new approach consists of the classical consideration of intermolecular motions of monomers, direct ab initio calculation of intermolecular potential, and the Monte Carlo integration during the partition function evaluation. The full six-dimensional intermolecular potential energy surfaces (PES) of these complexes required for the thermodynamic calculations and consisted of 14016 (H2OO3) and 9672 (H2OSO2) unique points were calculated by the ab initio (MP2/6-311++G(2d,2p)) method. The local minima found at the calculated PESs were used as starting points for the farther high-level geometry optimization (up to QCISD/aug-ccpVTZ and MP2/aug-ccpVQZ). It was found that the structure of the global minimum of H2OSO2 complex is in the agreement with the available experimental data whereas the global minimum of H2OO3 corresponds to the twisted asymmetric hydrogen-bonded structure (C1) different from that one proposed on the basis of microwave spectra and previous quantum chemical calculations. By using the new approach, the thermodynamic properties and the equilibrium constants of H2OO3, and H2OSO2 were obtained for the first time with explicit accounting their rovibrational dynamics. The best estimates for the K0(298) of the complexes H2OO3, and H2OSO2 are 1.0510-2 and 3.1510-2, respectively (empirically corrected values are 8.710-3 and 2.510-2). The corresponding uncertainty of these values is estimated as 25-35% (based on the known values of (H2O)2 and D2O)2 equilibrium constants).