We introduce PySCES, an LSC-TeraChem Python interface code to enable GPU-accelerated NAMD simulations using classical dynamics trajectories within an extended phase space of nuclear and electronic variables. Our interface makes use of the TeraChem protocol buffer client-server model and can be parallelized over multiple GPU nodes, opening the door for semiclassical NAMD simulations of excited state dynamics in the condensed phase.