We aim to investigate the molecular environment of the supernova remnant (SNR) G150.3+4.5, and explore its association with ambient molecular clouds (MCs). We find that MCs within the velocity range of [-14, -2]km/s are spatially distributed along the radio shell of the SNR. Line broadening and asymmetries are observed in the CO spectra of the clouds. Moreover, we find that the molecular clouds around the shell have systematic velocity gradients in the position-velocity (PV) diagram. Both the morphology alignment and gas kinematics suggest that the SNR is associated with the ambient MCs at approximately 740 pc. Based on the CO gas distance, the dimension and the age of the SNR are estimated to be 40pcx33pc and 38000 years, respectively. The very high energy emission of 1LHAASO J0428+5531 toward the SNR may originate from the interaction between the SNR and the surrounding MCs.