Python Libraries: Numerical Propagators and Access Calculations

I have written some of my own code in Matlab that performs two body unperturbed as well as J2 perturbed orbit propagation. I am shifting to a new environment that requires the use of Python, and can't seem to see if there are any libraries to that do numerical propagation and access calculations with Kepler orbital elements.

I am aware of libraries like astropy, but they all seem limited to just utilizing TLE's and no way of defining a custom orbit. I do not need super high fidelity.