diff --git a/bindings/python/algorithm/expose-kinematics-derivatives.cpp b/bindings/python/algorithm/expose-kinematics-derivatives.cpp
index af254fa889..d7b004e85f 100644
--- a/bindings/python/algorithm/expose-kinematics-derivatives.cpp
+++ b/bindings/python/algorithm/expose-kinematics-derivatives.cpp
@@ -113,9 +113,21 @@ namespace pinocchio
               "\tmodel: model of the kinematic tree\n"
               "\tdata: data related to the model\n");
 
+      bp::def("computeJointKinematicHessians",
+          &computeJointKinematicHessians<double, 0, JointCollectionDefaultTpl,
+                                         VectorXd>,
+          bp::args("model", "data", "q"),
+          "Computes the full model Jacobian, i.e. the stack of all the "
+          "motion subspaces expressed in the coordinate world frame.\n"
+          "The result is accessible through data.J. This function computes "
+          "also the forward kinematics of the model.\n\n"
+          "Parameters:\n"
+          "\tmodel: model of the kinematic tree\n"
+          "\tdata: data related to the model\n"
+          "\tq: the joint configuration vector (size model.nq)\n",
+          bp::return_value_policy<bp::return_by_value>());
+
     }
-    
-    
-    
+
   } // namespace python
 } // namespace pinocchio
diff --git a/include/pinocchio/bindings/python/multibody/data.hpp b/include/pinocchio/bindings/python/multibody/data.hpp
index 841c4022af..05deb4a5ab 100644
--- a/include/pinocchio/bindings/python/multibody/data.hpp
+++ b/include/pinocchio/bindings/python/multibody/data.hpp
@@ -173,6 +173,8 @@ namespace pinocchio
         .ADD_DATA_PROPERTY(dq_after,"Generalized velocity after the impact.")
         .ADD_DATA_PROPERTY(staticRegressor,"Static regressor.")
         .ADD_DATA_PROPERTY(jointTorqueRegressor,"Joint torque regressor.")
+
+        .ADD_DATA_PROPERTY(kinematic_hessians,"Joint Kinematic Hessians.")
         
         .def(bp::self == bp::self)
         .def(bp::self != bp::self)