Here, we have described the development of a Drude polarizable force field for PET polymers and its derived molecules. The parameter fitting protocol employed here led to accurate modeling of molecular polarizabilities, dipole moments, and water interactions by targeting QM data. Interactions between BHET and protein functional groups were evaluated, showing good agreement between gas-phase QM interaction energies and our Drude model. Lastly, new dihedral parameters were fit to QM configurational energy profiles, allowing an accurate reproduction of energy minima and barriers separating important configurations.