With the use of temporal derivative of flux function, a two-stage temporal discretization has been recently proposed in the design of fourth-order schemes based on the generalized Riemann problem (GRP) [21] and gas-kinetic scheme (GKS) [28]. In this paper, the fourth-order gas-kinetic scheme will be extended to solve the compressible multicomponent flow equations, where the two-stage temporal discretization and fifth-order WENO reconstruction will be used in the construction of the scheme. Based on the simplified two-species BGK model [41], the coupled Euler equations for individual species will be solved. Each component has its individual gas distribution function and the equilibrium states for each component are coupled by the physical requirements of total momentum and energy conservation in particle collisions. The second-order flux function is used to achieve the fourth-order temporal accuracy, and the robustness is as good as the second-order schemes. At the same time, the source terms, such as the gravitational force and the chemical reaction, will be explicitly included in the two-stage temporal discretization through their temporal derivatives. Many numerical tests from the shock-bubble interaction to ZND detonative waves are presented to validate the current approach.