diff --git a/src/pyscipopt/scip.pyi b/src/pyscipopt/scip.pyi
index 61ecf3073..e1926d7f8 100644
--- a/src/pyscipopt/scip.pyi
+++ b/src/pyscipopt/scip.pyi
@@ -236,12 +236,10 @@ class Expr:
     def __next__(self): ...
     def __pow__(self, other): ...
     def __radd__(self, other): ...
-    def __reduce_cython__(self, *args, **kwargs): ...
     def __rmul__(self, other): ...
     def __rpow__(self, other): ...
     def __rsub__(self, other): ...
     def __rtruediv__(self, other): ...
-    def __setstate_cython__(self, *args, **kwargs): ...
     def __sub__(self, other): ...
     def __truediv__(self, other): ...
 
@@ -274,12 +272,10 @@ class GenExpr:
     def __neg__(self): ...
     def __pow__(self, other): ...
     def __radd__(self, other): ...
-    def __reduce_cython__(self, *args, **kwargs): ...
     def __rmul__(self, other): ...
     def __rpow__(self, other): ...
     def __rsub__(self, other): ...
     def __rtruediv__(self, other): ...
-    def __setstate_cython__(self, *args, **kwargs): ...
     def __sub__(self, other): ...
     def __truediv__(self, other): ...
 
@@ -357,6 +353,7 @@ class MatrixExpr(numpy.ndarray):
     def __ge__(self, other: object) -> bool: ...
     def __iadd__(self, other): ...
     def __le__(self, other: object) -> bool: ...
+    def __matmul__(self, *args, **kwargs): ...
     def __mul__(self, other): ...
     def __pow__(self, other): ...
     def __radd__(self, other): ...
@@ -536,6 +533,7 @@ class Model:
     def getCurrentNode(self, *args, **kwargs): ...
     def getCutEfficacy(self, *args, **kwargs): ...
     def getCutLPSolCutoffDistance(self, *args, **kwargs): ...
+    def getCutoffbound(self, *args, **kwargs): ...
     def getDepth(self, *args, **kwargs): ...
     def getDualMultiplier(self, *args, **kwargs): ...
     def getDualSolVal(self, *args, **kwargs): ...
@@ -559,7 +557,9 @@ class Model:
     def getLhs(self, *args, **kwargs): ...
     def getLinearConsIndicator(self, *args, **kwargs): ...
     def getLocalEstimate(self, *args, **kwargs): ...
+    def getLowerbound(self, *args, **kwargs): ...
     def getMajorVersion(self, *args, **kwargs): ...
+    def getMaxDepth(self, *args, **kwargs): ...
     def getMinorVersion(self, *args, **kwargs): ...
     def getNBestSolsFound(self, *args, **kwargs): ...
     def getNBinVars(self, *args, **kwargs): ...
@@ -581,12 +581,14 @@ class Model:
     def getNLeaves(self, *args, **kwargs): ...
     def getNLimSolsFound(self, *args, **kwargs): ...
     def getNNlRows(self, *args, **kwargs): ...
+    def getNNodeLPIterations(self, *args, **kwargs): ...
     def getNNodes(self, *args, **kwargs): ...
     def getNReaders(self, *args, **kwargs): ...
     def getNSepaRounds(self, *args, **kwargs): ...
     def getNSiblings(self, *args, **kwargs): ...
     def getNSols(self, *args, **kwargs): ...
     def getNSolsFound(self, *args, **kwargs): ...
+    def getNStrongbranchLPIterations(self, *args, **kwargs): ...
     def getNTotalNodes(self, *args, **kwargs): ...
     def getNVars(self, *args, **kwargs): ...
     def getNVarsAnd(self, *args, **kwargs): ...
@@ -602,6 +604,7 @@ class Model:
     def getOpenNodes(self, *args, **kwargs): ...
     def getParam(self, *args, **kwargs): ...
     def getParams(self, *args, **kwargs): ...
+    def getPlungeDepth(self, *args, **kwargs): ...
     def getPresolvingTime(self, *args, **kwargs): ...
     def getPrimalRay(self, *args, **kwargs): ...
     def getPrimalRayVal(self, *args, **kwargs): ...