Problems with accessing ODE parameters after optimisation when modelling with ModelingToolkit

Thanks! Yes, I tried this option but I unfortunately got the following error messages:

println(results.ps[protein_model.x1_act])
ERROR: LoadError: MethodError: no method matching parameter_values(::SciMLBase.NullParameters)
An arithmetic operation was performed on a NullParameters object. This means no parameters were passed
into the AbstractSciMLProblem (e.x.: ODEProblem) but the parameters object `p` was used in an arithmetic
expression. Two common reasons for this issue are:

1. Forgetting to pass parameters into the problem constructor. For example, `ODEProblem(f,u0,tspan)` should
be `ODEProblem(f,u0,tspan,p)` in order to use parameters.

2. Using the wrong function signature. For example, with `ODEProblem`s the function signature is always
`f(du,u,p,t)` for the in-place form or `f(u,p,t)` for the out-of-place form. Note that the `p` argument
will always be in the function signature regardless of if the problem is defined with parameters!



Closest candidates are:
  parameter_values(::Any, ::Any)
   @ SymbolicIndexingInterface ~/.julia/packages/SymbolicIndexingInterface/Pz8Cf/src/parameter_indexing.jl:16
  parameter_values(::SciMLBase.AbstractJumpProblem)
   @ SciMLBase ~/.julia/packages/SciMLBase/Dwomw/src/problems/problem_interface.jl:12
  parameter_values(::ModelingToolkit.MTKParameters, ::ModelingToolkit.ParameterIndex)
   @ ModelingToolkit ~/.julia/packages/ModelingToolkit/Mxj1Q/src/systems/parameter_buffer.jl:216
  ...

Stacktrace:
 [1] parameter_values(prob::SciMLBase.NullParameters, i::Nothing)
   @ SymbolicIndexingInterface ~/.julia/packages/SymbolicIndexingInterface/Pz8Cf/src/parameter_indexing.jl:16
 [2] parameter_values(A::SciMLBase.OptimizationSolution{…}, i::Nothing)
   @ SciMLBase ~/.julia/packages/SciMLBase/Dwomw/src/solutions/solution_interface.jl:44
 [3] _getter
   @ ~/.julia/packages/SymbolicIndexingInterface/Pz8Cf/src/parameter_indexing.jl:117 [inlined]
 [4] getter
   @ ~/.julia/packages/SymbolicIndexingInterface/Pz8Cf/src/parameter_indexing.jl:139 [inlined]
 [5] getindex(::SymbolicIndexingInterface.ParameterIndexingProxy{SciMLBase.OptimizationSolution{…}}, ::Num)
   @ SymbolicIndexingInterface ~/.julia/packages/SymbolicIndexingInterface/Pz8Cf/src/parameter_indexing_proxy.jl:14
 [6] top-level scope
   @ ~/Documents/Master/MasterThesis/RhoAModelling/Notebooks/parameters.jl:84
 [7] include(fname::String)
   @ Base.MainInclude ./client.jl:489
 [8] top-level scope
   @ REPL[2]:1

Or on the ODE solution:

sol = solve(remake(simul_data_prob, p=results.u), saveat=0.5, dtmax = 0.5)
println(sol.ps[protein_model.x1_act])
ERROR: LoadError: ArgumentError: invalid index: ModelingToolkit.ParameterIndex{SciMLStructures.Tunable, Tuple{Int64, Int64}}(SciMLStructures.Tunable(), (1, 2)) of type ModelingToolkit.ParameterIndex{SciMLStructures.Tunable, Tuple{Int64, Int64}}
Stacktrace:
  [1] to_index(i::ModelingToolkit.ParameterIndex{SciMLStructures.Tunable, Tuple{Int64, Int64}})
    @ Base ./indices.jl:300
  [2] to_index(A::Vector{Float64}, i::ModelingToolkit.ParameterIndex{SciMLStructures.Tunable, Tuple{Int64, Int64}})
    @ Base ./indices.jl:277
  [3] _to_indices1(A::Vector{Float64}, inds::Tuple{Base.OneTo{…}}, I1::ModelingToolkit.ParameterIndex{SciMLStructures.Tunable, Tuple{…}})
    @ Base ./indices.jl:359
  [4] to_indices
    @ ./indices.jl:354 [inlined]
  [5] to_indices
    @ ./indices.jl:345 [inlined]
  [6] getindex
    @ ./abstractarray.jl:1291 [inlined]
  [7] parameter_values(arr::Vector{Float64}, i::ModelingToolkit.ParameterIndex{SciMLStructures.Tunable, Tuple{Int64, Int64}})
    @ SymbolicIndexingInterface ~/.julia/packages/SymbolicIndexingInterface/Pz8Cf/src/parameter_indexing.jl:15
  [8] parameter_values(A::ODESolution{…}, i::ModelingToolkit.ParameterIndex{…})
    @ SciMLBase ~/.julia/packages/SciMLBase/Dwomw/src/solutions/solution_interface.jl:44
  [9] (::SymbolicIndexingInterface.var"#_getter#9"{…})(::SymbolicIndexingInterface.Timeseries, prob::ODESolution{…})
    @ SymbolicIndexingInterface ~/.julia/packages/SymbolicIndexingInterface/Pz8Cf/src/parameter_indexing.jl:120
 [10] (::SymbolicIndexingInterface.var"#getter#10"{…})(::ODESolution{…})
    @ SymbolicIndexingInterface ~/.julia/packages/SymbolicIndexingInterface/Pz8Cf/src/parameter_indexing.jl:139
 [11] getindex(::SymbolicIndexingInterface.ParameterIndexingProxy{ODESolution{…}}, ::Num)
    @ SymbolicIndexingInterface ~/.julia/packages/SymbolicIndexingInterface/Pz8Cf/src/parameter_indexing_proxy.jl:14
 [12] top-level scope
    @ ~/Documents/Master/MasterThesis/RhoAModelling/Notebooks/parameters.jl:88
 [13] include(fname::String)
    @ Base.MainInclude ./client.jl:489
 [14] top-level scope
    @ REPL[4]:1