ABSTRACT We decomposed the molecular gas in the spiral galaxy NGC 628 (M74) into multiscale hub-filament structures using the CO (2$-$1) line by the dendrogram algorithm. All leaf structures as potential hubs were classified into three categories, i.e. leaf-HFs-A, leaf-HFs-B and leaf-HFs-C. Leaf-HFs-A exhibit the best hub-filament morphology, which also have the highest density contrast, the largest mass and the lowest virial ratio. We employed the filfinder algorithm to identify and characterize filaments within 185 leaf-HFs-A structures, and fitted the velocity gradients around the intensity peaks. Measurements of velocity gradients provide evidence for gas inflow within these structures, which can serve as a kinematic evidence that these structures are hub-filament structures. The numbers of the associated 21 μm and H α structures and the peak intensities of 7.7 μm, 21 μm, and H α emissions decrease from leaf-HFs-A to leaf-HFs-C. The spatial separations between the intensity peaks of CO and 21 μm structures of leaf-HFs-A are larger than those of leaf-HFs-C. These evidence indicate that leaf-HFs-A are more evolved than leaf-HFs-C. There may be an evolutionary sequence from leaf-HFs-C to leaf-HFs-A. Currently, leaf-HFs-C lack a distinct gravitational collapse process that would result in a significant density contrast. The density contrast can effectively measure the extent of the gravitational collapse and the depth of the gravitational potential of the structure which, in turn, shapes the hub-filament morphology. Combined with the kinematic analysis presented in previous studies, a picture emerges that molecular gas in spiral galaxies is organized into network structures through the gravitational coupling of multiscale hub-filament structures. Molecular clouds, acting as knots within these networks, serve as hubs, which are local gravitational centres and the main sites of star formation.