Drug resistance impacts the effectiveness of many new therapeutics. Mutations in the therapeutic target confer resistance; however, deciphering which mutations, often remote from the enzyme active site, drive resistance is challenging. In a series of Pneumocystis jirovecii dihydrofolate reductase variants, we elucidate which interactions are key bellwethers to confer resistance to trimethoprim using homology modeling, molecular dynamics, and machine learning. Six molecular features involving mainly residues that did not vary were the best indicators of resistance.