diff --git a/Cargo.toml b/Cargo.toml
index 93a499a..9108b3e 100644
--- a/Cargo.toml
+++ b/Cargo.toml
@@ -30,6 +30,7 @@ num-bigint = "0.4.6"
 num-traits = "0.2.19"
 either = "1.13.0"
 itertools = "0.13.0"
+clap = "4.5.21"
 
 # Include Z3 dependencies as strictly as possible, we don't want
 # this to change because it might break our release builds.
diff --git a/conda/meta.yaml b/conda/meta.yaml
index 6cb6a09..81162e9 100644
--- a/conda/meta.yaml
+++ b/conda/meta.yaml
@@ -1,5 +1,5 @@
 {% set name = "biodivine_aeon" %}
-{% set version = "1.2.0" %}
+{% set version = "1.2.1" %}
 
 package:
   name: "{{ name|lower }}"
diff --git a/src/bindings/pbn_control/control.rs b/src/bindings/pbn_control/control.rs
index 03d8348..a7434c3 100644
--- a/src/bindings/pbn_control/control.rs
+++ b/src/bindings/pbn_control/control.rs
@@ -138,6 +138,10 @@ impl Control {
     /// all network colors). When this option is active, other results that have been computed
     /// so far are still returned.
     ///
+    /// Optionally, you can also provide a set of `initial_states` that will be considered as the
+    /// only potential initial states of the network (e.g. in case of some initial states are not
+    /// considered valid).
+    ///
     #[staticmethod]
     #[pyo3(signature = (graph, phenotype, oscillation_type = None, size_limit = None, stop_when_found = false, initial_states = None))]
     pub fn phenotype_permanent(