The conformational preference of tetrahomodioxacalix[4]arenes with three different para substituents on the phenolic ring has been investigated by using ab initio molecular orbital theory (RHF/6-31$G^{ast}$) and density functional theory (B3LYP/6-31$G^{ast}$). The stability order is predicted to be cone > C-1,2-alternate > partial cone > 1,3-alternate > COC-1,2-alternate. The distorted cone conformation is found to be most stable in a gas phase and the calculated results are in agreement with the reported $^1$H NMR and X-ray experimental observations. The substitution of methylene with dimethyleneoxa bridges increases the size of the annulus of the molecule, its conformational mobility, and the number of hydrogen bonding patterns. The thermodynamic stability and the conformational characteristics of tetrahomodioxacalix[4]arenes are discussed in regards of the number of phenolic hydrogen bonding patterns and the polarity of a molecule. The substituent effects on the para position of the phenolic ring are also introduced.