Many naturally occurring elastomers are intrinsically disordered proteins (IDPs) built up of repeating units and they can demonstrate two types of thermoresponsive phase behavior. Systems characterized by lower critical solution temperatures (LCST) undergo phase separation above the LCST whereas systems characterized by upper critical solution temperatures (UCST) undergo phase separation below the UCST. There is congruence between thermoresponsive coil-globule transitions and phase behavior. Specifically, the theta temperatures above or below which the IDPs transition from coils to globules serve as useful proxies for the LCST / UCST values. This implies that one can design sequences with desired values for the theta temperature with either increasing (UCST) or decreasing radii of gyration (LCST) above the theta temperature. Here, we show that the Monte Carlo simulations performed in the so-called intrinsic solvation (IS) limit version of the temperature-dependent ABSINTH implicit solvation model, yields a robust heuristic for discriminating between sequences with known LCST versus UCST phase behavior. Accordingly, we use this heuristic in a supervised approach, integrate it with a genetic algorithm, combine this with IS limit simulations, and show how novel sequences can be designed that have LCST phase behavior. These calculations are aided by direct estimates of temperature dependent free energies of solvation for model compounds that are derived using the polarizable AMOEBA forcefield. To demonstrate the validity of our designs, we calculate coil-globule transition profiles using the full ABSINTH model and combine these with the Gaussian Cluster Theory to show that the designed IDPs do show LCST phase behavior.