Calculating e (base of the natural log) to high precision in Python?

Is it possible to calculate the value of the mathematical constant, e with high precision (2000+ decimal places) using Python?

I am particularly interested in a solution either in or that integrates with NumPy or SciPy.