The present article describes techniques for classical simulations of proteins and protein-nucleic acid complexes, revealing their dynamics and protein-substrate binding energies. The approach is based on classical atomistic molecular dynamics (MD) simulations of the experimentally determined structures of the complexes. MD simulations can provide dynamics of complexes in realistic solvents on microsecond timescales, and the free energy methods are able to provide Gibbs free energies of binding of substrates, such as nucleic acids, to proteins. The chapter describes methodologies for the preparation of computer models of biomolecular complexes and free energy perturbation methodology for evaluating Gibbs free energies of binding. The applications are illustrated with examples of snapshots of proteins and their complexes with nucleic acids, as well as the precise Gibbs free energies of binding.