How to compute a bounding box of a TopoDS_Shape in a specific coordinate system (not global)

Hello everyone,

I am trying to compute the bounding box of a TopoDS_Shape, but I need the box to be expressed in a specific local coordinate system, not in the global one.
For example, I have a transformation (represented by a gp_Ax2 or a gp_Trsf) that defines the local coordinate system in which I want the bounding box to be calculated.
My question is:
What is the recommended way in OpenCascade to compute the bounding box of a shape with respect to a custom coordinate system?
I would like to get the bounding box dimensions and position relative to that local system, not the global one.
Any advice or example code would be greatly appreciated.
Thank you!