I read it as they are using a molecular dynamics package to define a function (potential energy of the protein structure) and are using gradient descent to (slowly) reach stable extrema (static configurations of the protein).
That makes more sense, but the energy function isn't really a feature of molecular dynamics, anything that does molecular optimization or modeling will have this information built in.