Accumulating clinical studies have consistently demonstrated that the microbes in the human body closely interact with the human host, actively participating in the regulation of drug effectiveness. Identifying the associations between microbes and drugs can facilitate the development of drug discovery, and microbes have become a new target in antimicrobial drug development. However, the discovery of microbe-drug associations relies on clinical or biological experiments, which are not only time-consuming but also financially burdensome. Thus, the utilization of computational methods to predict microbe-drug associations holds promise for reducing costs and enhancing the efficiency of biological experiments. Here, we introduce a new computational method, called HKFGCN (Heterogeneous information Kernel Fusion Graph Convolution Network), to predict the microbe-drug associations. Instead of extracting feature from a single network in previous studies, HKFGCN separately extracts topological information features from different networks, and further refines them by generating Gaussian kernel features. HKFGCN consists of three main steps. Firstly, we constructed two similarity networks and a microbe-drug association network based on numerous biological data. Second, we employed two types of encoders to extract features from these networks. Next, Gaussian kernel features were obtained from the drug and microbe features at each layer. Finally, we reconstructed the bipartite microbe-drug graph based on the learned representations. Experimental results demonstrate the excellent performance of the HKFGCN model across different datasets using the cross-validation scheme. Additionally, we conduced case studies on human immunodeficiency virus, and the results were corroborated by existing literatures. The prediction model’s code is available at https://github.com/roll-of-bubble/HKFGCN.