We employ microsecond atomic-scale molecular dynamics simulations to get insight into the structural and thermal properties of heat-resistant bulk polyimides. As electrostatic interactions are essential for the polyimides considered, we propose a two-step equilibration protocol that includes long (microsecond-scale) MD simulations of polymer melt with partial atomic charges switched off, followed by relatively short runs (100 ns) of the polymer system with full electrostatics. We demonstrate that macroscopic characteristics of a polyimide sample (such as the glass transition temperature and density) are not particularly sensitive to the degree of equilibration. However, great caution should be paid when local structural characteristics are considered: proper equilibration of the local polymer structure (monitored through the radius of gyration and the end-to-end distance of individual chains) is found to require simulations on a microsecond time scale. Finally, we found a dramatic impact of electrostatic interactions on the properties of the bulk polyimides considered: when intra- and intermolecular dipole-dipole interactions come into play we witness compaction of individual polymer coils and eventually an increase in the glass transition temperature and polymer density.