Find p-value of model in python

Below is my code and I'm looking to find the p-value for the model that was created, kindly help me with this!

combined_data = combined_data.fillna(0)
predictors_1 = combined_data[['total_cases','new_cases','new_deaths','Transit Mobility 7-Day Avg']]
response_1 = combined_data['Reproduction 7-Day Avg']

regression_1 = linear_model.LinearRegression()
model_1 = regression_1.fit(predictors_1,response_1)

print('Coefficients: ',model_1.coef_,'Intercept: ',model_1.intercept_,'\nR-Squared: ',model_1.score(predictors_1,response_1))