diff --git a/Project.toml b/Project.toml
index 842cef8..1a60996 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "SBML"
 uuid = "e5567a89-2604-4b09-9718-f5f78e97c3bb"
-authors = ["The developers of SBML.jl"]
 version = "1.6"
+authors = ["The developers of SBML.jl"]
 
 [deps]
 DocStringExtensions = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"
@@ -16,7 +16,7 @@ ConstructionBase = "1.3"
 DocStringExtensions = "0.8, 0.9"
 IfElse = "0.1"
 SBML_jll = "5.19.5"
-Symbolics = "5,6"
+Symbolics = "7"
 Unitful = "1"
 julia = "1.6"
 
diff --git a/test/loadmodels.jl b/test/loadmodels.jl
index a843fdf..7029ba8 100644
--- a/test/loadmodels.jl
+++ b/test/loadmodels.jl
@@ -345,7 +345,14 @@ end
     m = readSBML(joinpath(@__DIR__, "data", "01565-sbml-l3v1.xml"))
 
     if TEST_SYMBOLICS
-        @test interpret_as_num(m.reactions["J23"].kinetic_math) == 0.0
+        # Older versions of symbolics somehow auto-simplify this, newer ones do not.
+        @variables S23 S23b
+        @test (
+            isequal(interpret_as_num(m.reactions["J23"].kinetic_math), 0.0) || isequal(
+                interpret_as_num(m.reactions["J23"].kinetic_math),
+                S23 * S23b * log(Num(1)) / log(Num(10)),
+            )
+        )
 
         @variables S29 S29b
         @test isequal(interpret_as_num(m.reactions["J29"].kinetic_math), 2.0 * S29 * S29b)