diff --git a/python/sdist/amici/setup.template.py b/python/sdist/amici/setup.template.py
index 1f2d099110..80fd3fde2a 100644
--- a/python/sdist/amici/setup.template.py
+++ b/python/sdist/amici/setup.template.py
@@ -72,7 +72,6 @@ def get_extension() -> CMakeExtension:
     ext_modules=[MODEL_EXT],
     packages=find_namespace_packages(),
     install_requires=["amici==TPL_AMICI_VERSION"],
-    extras_require={"wurlitzer": ["wurlitzer"]},
     python_requires=">=3.10",
     package_data={},
     zip_safe=False,
diff --git a/python/sdist/amici/swig_wrappers.py b/python/sdist/amici/swig_wrappers.py
index 5762f50bec..2e90891df3 100644
--- a/python/sdist/amici/swig_wrappers.py
+++ b/python/sdist/amici/swig_wrappers.py
@@ -1,9 +1,7 @@
 """Convenience wrappers for the swig interface"""
 
 import logging
-import sys
 import warnings
-from contextlib import contextmanager, suppress
 from typing import Any
 
 import amici
@@ -30,22 +28,6 @@
     "get_model_settings",
 ]
 
-try:
-    from wurlitzer import sys_pipes
-except ModuleNotFoundError:
-    sys_pipes = suppress
-
-
-@contextmanager
-def _capture_cstdout():
-    """Redirect C/C++ stdout to python stdout if python stdout is redirected,
-    e.g. in ipython notebook"""
-    if sys.stdout == sys.__stdout__:
-        yield
-    else:
-        with sys_pipes():
-            yield
-
 
 def runAmiciSimulation(
     model: AmiciModel,
@@ -82,10 +64,9 @@ def runAmiciSimulation(
             stacklevel=1,
         )
 
-    with _capture_cstdout():
-        rdata = amici_swig.runAmiciSimulation(
-            _get_ptr(solver), _get_ptr(edata), _get_ptr(model)
-        )
+    rdata = amici_swig.runAmiciSimulation(
+        _get_ptr(solver), _get_ptr(edata), _get_ptr(model)
+    )
     _log_simulation(rdata)
     if solver.getReturnDataReportingMode() == amici.RDataReporting.full:
         _ids_and_names_to_rdata(rdata, model)
@@ -124,15 +105,14 @@ def runAmiciSimulations(
             stacklevel=1,
         )
 
-    with _capture_cstdout():
-        edata_ptr_vector = amici_swig.ExpDataPtrVector(edata_list)
-        rdata_ptr_list = amici_swig.runAmiciSimulations(
-            _get_ptr(solver),
-            edata_ptr_vector,
-            _get_ptr(model),
-            failfast,
-            num_threads,
-        )
+    edata_ptr_vector = amici_swig.ExpDataPtrVector(edata_list)
+    rdata_ptr_list = amici_swig.runAmiciSimulations(
+        _get_ptr(solver),
+        edata_ptr_vector,
+        _get_ptr(model),
+        failfast,
+        num_threads,
+    )
     for rdata in rdata_ptr_list:
         _log_simulation(rdata)
         if solver.getReturnDataReportingMode() == amici.RDataReporting.full:
diff --git a/python/sdist/pyproject.toml b/python/sdist/pyproject.toml
index 7c903273b5..10a3b05374 100644
--- a/python/sdist/pyproject.toml
+++ b/python/sdist/pyproject.toml
@@ -31,7 +31,6 @@ dependencies = [
     "python-libsbml",
     "pandas>=2.0.2",
     "pyarrow",
-    "wurlitzer",
     "toposort",
     "setuptools>=48",
     "mpmath",