The design of bioelectrochemical interfaces (BEI) is an interesting topic that recently demands attention. The synergy between biomolecules and chemical components is necessary to achieve high molecular selectivity and sensitivity for the development of biosensors, synthesis of different compounds, or catalytic processes. For most BEI, the charge transfer process occurs in environments with particular chemical conditions; modeling these environments is a challenging task and requires multidisciplinary efforts. These interfaces can be composed of biomolecules, such as proteins, DNA, or more complex systems like microorganisms. Oxidoreductases enzymes are good candidates, among others, due to their catalytic activities and structural characteristics. In BEI, enzymes are immobilized on conductive surfaces to improve charge transfer processes. Covalent immobilization is the most common method to prolong lifetime or modulate the detection process. However, it is necessary to implement new methodologies that allow the selection of the best candidates for a more efficient design. Homology modeling of oxidoreductases combined with Molecular Dynamics (MD) simulation methods are alternative and already routinely used tools to investigate the structure, dynamics, and thermodynamics of biological molecules. Our motivation is to show different techniques of molecular modeling (Homology Modeling, Gaussian accelerated molecular dynamics, directed adaptive molecular dynamics and electrostatic surface calculations), and using horseradish peroxidase as a model to understand the interactions between biomolecules and gold nanoclusters (as current collector). Additionally, we present our previous studies considering molecular simulations and we discuss recent advances in biomolecular simulations aimed at biosensor design.