Open access peer-reviewed chapter

Design of Bioelectrochemical Interfaces Assisted by Molecular Dynamics Simulations

Written By

Abraham Vidal-Limon, Guillermo Antonio Huerta-Miranda, Wendy I. García-García and Margarita Miranda-Hernández

Submitted: July 7th, 2020 Reviewed: September 4th, 2020 Published: September 25th, 2020

DOI: 10.5772/intechopen.93884

From the Edited Volume

Homology Molecular Modeling

Edited by Rafael Trindade Maia, Rômulo Maciel de Moraes Filho and Magnólia Campos

Chapter metrics overview

614 Chapter Downloads

View Full Metrics


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.


  • bioelectrochemical interfaces
  • homology modeling
  • covalent immobilization
  • gold nanoclusters
  • molecular dynamics

1. Introduction

For some years now, the design and construction of bioelectrochemical interfaces (BEI), ranging from electrochemical biosensors (EC) for analytical applications, biofuel cells (BFC) to the development of biocomputing systems for information processing, have been topics where the scientific community has considerable participation. The implementation and integration of biomolecules with electronic elements or conductive surfaces to produce bioelectronic devices is a current topic and gathers great importance for developing biosensors that have an essential role in clinical applications, food quality control analysis and forensic medicine [1, 2, 3]. These applications always involve high sensitivity, low-sample volume, and low-cost production. In most BEI, charge transfer processes are carried out in environments with high ionic concentration, necessary for their function. These interfaces can be composed of macromolecules such as proteins, peptides, or more complex systems such as bacteria. Since the activity and lifetime depend on the correct interaction between the conductive surfaces (current collector) and the biomaterials, to achieve good bioelectrochemical responses, the enzymes need to improve the orientation towards the surface of the current collector [4, 5].

1.1 Covalent immobilization through alkanethiol linkers

The immobilization of various biomaterials is an important issue for BEI design, the covalent immobilization of enzymes stands out among the different enzymatic coupling strategies since the bonds formed through the linkers anchored to the current collector promotes direct charge transfer responses [6, 7]. Alcakanothiols are organic molecules widely used to establish self-assembled monolayers (SAMs) on the surface of gold electrodes [8, 9, 10, 11]. The thiol groups chemisorb on gold electrodes forming strong thiolate-gold bonds (Au-S). The resulting monolayer may be used to design molecular scaffolds to couple enzymes [8, 12]. The main drawback of covalent immobilization systems is the risk of having a low surface concentration of enzymes that are either active or correctly oriented for direct charge transfer; which may lead to a low-efficiency bioelectrochemical response [13]. Therefore, having active enzymes and a good bioelectrochemical response becomes an important task. Aromatic molecules like 4-aminothiophenol, structurally mimics enzymatic substrates and can promote stable and direct contacts, necessary for efficient bioelectrochemical reactions.

Researches have proposed nanomaterials as support matrices for the design of BEI to increase the surface concentration of active enzymes on the current collector. During the last decade, nanomaterials, coupled with enzymes, have had significant relevance in the design of biosensors [5, 14, 15, 16]. Current advances in synthesis methodologies of these materials allow having a wide variety of nanomaterials with different sizes, shapes, surface charges, and physicochemical characteristics [4, 14, 17].

Some of these nanomaterials can be modified in different ways to improve biocompatibility. Biocompatibility refers to the biological nature events that do not interfere with those of electronic signal transduction and vice versa. Metallic nanoparticles are promising materials that increase the electroactive area and improve the sensitivity and stability of the attached enzymes on the electrodes, bringing the enzymatic active site close to the electrode (the redox cofactor should not exceed a 20 Å distance) to achieve direct charge transfer reactions [14, 18, 19, 20, 21, 22, 23].

In our group, we have established the conditions to improve the electrochemical response of horseradish peroxidase (HRP) coupled to gold nanoparticles modified with the aromatic alkanothiol 4-aminotiophenol, using glassy carbon (GC) as current collector [24]. We proposed that a critical factor in HRP bioelectrochemical response is to promote conformational changes and proper enzyme orientation when coupled to the electrode. The protein conformation around the prosthetic group will determine the redox potential of the Fe3+/Fe2+ couple. The activity towards H2O2 also depends on the protein structural arrangements [25, 26]. We found that a proper environment for the enzyme activity was achieved by increasing the distance between gold nanoparticles. The voltammetric studies of the prosthetic heme group showed significant differences between the enzyme immobilized on randomly deposited gold nanoparticles and the enzyme immobilized on a well-dispersed gold nanoparticles deposit (Figure 1). The dispersed electrode improved the electrochemical response of the enzyme; this fact showed that the distance between enzymes is important, probably because longer distances decreased the steric impediments between the enzymes, and it appears to be a better immobilization strategy. With these results, we proposed that an essential part of BEI design is related to the structure of the biological systems. The electrochemical sensing is dependent on the recognition or transformation of the substrate, as previously observed [27].

Figure 1.

Cyclic voltammogram of HRP response at different scan rates, immobilized on GC electrode modified with gold nanoparticles functionalized with the alkanethiol molecule 4-aminothiophenol. Phosphate buffer pH 6.8.

The advent of high-resolution and robust techniques as X-ray crystallography [28, 29] and cryo-EM [30, 31, 32] has enormously contributed to the baggage of structural information available for bioelectrochemical applications. These techniques describe the atomic positions of enzymes, as well as the conformational dynamics resulted from recognition and binding processes, allowing the prediction of the possible electrochemical behavior of enzymes immobilized on electrodes.

1.2 Molecular modeling methodologies

The use of computational tools such as molecular mechanics and quantum mechanics for the study of chemical and biological reactions, involves a mathematical treatment of a large number of particles (hundreds of thousands of atoms that build up molecules like proteins, nucleic acids, lipids, and sugars) [33, 34]. In general terms, these methods are ideal for obtaining chemical and physical properties from three-dimensional molecular models [35, 36, 37]. It is necessary to correlate the structural models with their corresponding chemical, catalytic or biological properties; towards the rationalization and design of molecules for specific applications. Mainly, molecular Dynamics (MD) methods are based on the solution of Newton’s second law for all the atoms in the system, and help to predict the behavior of a biomolecule during a specific time. The integration of motion equations allows to analyze the trajectories corresponding to position, velocity, and acceleration of each particle of the system in any fraction of time [38]. Unfortunately, the simplification of considering atoms as spheres, and bonds as springs or tensors (with different parameters depending on the type of bond), makes difficult to study bond-breaking and formation, or reaction mechanisms. Quantum mechanics (QM) mathematically describes the fundamental behavior of matter on an electronic scale by solving the Schrödinger equation for each atom in the system. QM also describes the behavior of atoms and molecules, in terms of their chemical reactivity, geometry, and their optical, electrical, magnetic, and mechanical properties [39, 40]. The coupled implementation of these techniques could describe the biomolecules involved in charge transfer processes, and design strategies for efficient BEI.

1.3 Molecular dynamics as a tool to predict the electrochemical activity in proteins

A challenging task in BEI design is the selection of robust biomolecules capable of tolerate non-physiological conditions, without loosing efficiency and conformation [2]. The protein folding is governed by a series of molecular interactions between the amino acids and the surrounding chemical environment [34]. Therefore, predicting their behavior with molecular models contributes to the optimization of resources before the development of BEI.

In our previous work the VP6 rotavirus capsid protein was encapsulated with the ionic polymer Nafion on GC. We demonstrated that this electrode could transfer electric charge applying an external potential, when using the redox probe potassium ferricyanide (K3[Fe(CN)6]). In parallel, we modeled by MD the electrostatic and conformational states of VP6 (Figure 2). This model suggested a route for ionic conductivity, where the electrostatic protein surface displayed a negative charge which interacted with the ferricyanide redox probe, promoting the charge transfer reaction. In order to show if this electrochemical activity was particular to proteins in general, under the same conditions as the experiments with VP6, bovine serum albumin (BSA), whose primary biological function is to transport different molecules through the bloodstream, was tested. However, we did not observe the charge transfer reaction of potassium ferricyanide when this protein was encapsulated (Figure 2). These results demonstrated that VP6 protein could be used as conductive scaffold for the development of different BEI applications [41].

Figure 2.

Comparison of the cyclic Voltammogram of VP6 (pink) and BSA (blue) coupled to GC, at 20 mV/s in 0.01 M K3[Fe(CN)6], 1 M KCl.

Calculations of polarized electrostatic surface potentials from homology models of viral proteins, surprisingly showed that several capsids could display field effects as the recorded on VP6. In contrast, this effect cannot be displayed on non-charged proteins as BSA (Figure 3). Therefore, our previous data on viral capsids can serve as a workflow to identify candidate proteins or enzymes for construction of BEI devices. Altogether, theoretical and experimental reports on BEI [41, 42, 43] are examples of how structural information could help to elucidate the behavior of biological systems during electrochemical reactions. Nonetheless, the aforementioned studies depend on structural information from biomolecules, and certain biomolecules are very complex and highly labile, which difficult elucidation of atomic positions. However, since new protein sequences are continuously available, predictive methodologies as homology modeling becomes a valuable tool to asset structural information.

Figure 3.

Polarized electrostatic Isosurface potentials of different viral capsids. (A) Bovine serum albumin protein, (B) rotavirus VP6 trimer capsid protein, (C) influenza virus capsid protein, (D) hepatitis B virus capsid protein, (E) HIV-1 virus capsid protein, (F) HPV virus capsid protein. All polarized electrostatic Isosurfaces were calculated with PME approximation on capsid arrangements and displayed as positive (blue) or negative (red) mesh isosurfaces with VMD v.1.9.3 [44].

The aim of this multidisciplinary research is to perform molecular simulations on horseradish peroxidase as a model system, to generate data that can be used as selection criteria for design and further experimental validation in BEI development.


2. General techniques: Horseradish modeling and system building

The X-ray crystal structure from horseradish peroxidase (HRP) was collected from Protein Data Bank (, PDB code 1H5A) [45], and reconstructed using i-tasser metaserver [46]. The best model was selected based on lowest RMSD and higher C-value, to achieve the best initial coordinates. To improve the model system, all missing atoms and sidechain positions were modeled against available plant peroxidases on the PDB. A penta-coordinated resting state heme group was crystallized as Fe3+, and axial coordination with a histidine groups was retained.

The Protein Preparation Wizard routine (Schrödinger Maestro v2019-4, New York, 2019), was applied as a preparation step to correctly assign missing hydrogen atoms and protonation states of ionizable residues. Finally, PROPKA2.0 sub-routine was applied to whole model system [47] and pH 7 was fixed accordingly to our previous studies.

Enzyme structural parameters were described using the ff14SB force field [48]. The full systems comprised ~60,000 atoms in a cubic box of 15 Å length including TIP3P explicit water and Na + ions to ensure overall charge neutrality. Non-bonded interactions were calculated within a 12 Å cutoff, and long-range electrostatics were treated using the Particle-Mesh Ewald method [49]. SHAKE algorithm was enabled to constrain all bonds involving hydrogen during simulations. Conventional MD protocol for each system replica comprised: 1) 5000 steepest descent minimization steps followed by 10,000 conjugate gradient minimization steps; 2) 500 ps of progressive NVT heating from 0 to 300 K; 3) 5 ns of NVT equilibration, and 4) 10 ns of NPT dynamics at 300 K and 1 bar.

2.1 Gaussian accelerated molecular dynamics (GaMD)

Conventional MD simulations were carried out to initially relax full system coordinates, using AMBER18 software package [50]. Then, the GaMD module implemented in AMBER v. 18 was applied to perform extra 50 ns of short classic MD simulation to collect the statistics for calculating GaMD acceleration parameters. Extra 50 ns of short equilibration was applied after adding the boost potential, and finally three independent 500 ns GaMD production simulations with randomized initial atomic velocities. This method applies Gaussian functions (boost potentials) on the total potential energy term, which disturb the potential energy surface allowing the exploration of different energy states. In addition, functions that directly disturb the dihedral angles of the amide bonds were applied, promoting conformational changes. All GaMD simulations were performed with a dual-boost level by setting the reference energy to the lower bound. One boost potential is applied to the dihedral energetic term and the other to the potential energy term. The average and SD of the system potential energies were calculated every 500,000 steps (1 ns) for all simulation systems. The upper limit of the boost potential SD, σ0, was set to 6.0 kcal·mol−1 for both the dihedral and the total potential energetic terms [51, 52, 53]. The system temperature was ~298 K and 1 atm pressure, with integration steps of 2 fs. Three different system coordinates were extracted from GaMD trajectories, each one representing a different energy and conformational state.

2.2 Building of nanomaterial structures

The selected nanomaterials were Au nanoclusters (AuNCs) modified with two types of linkers: 4-aminothiophenol (ATP) and mercaptobenzoic acid (MBA); each cluster had 2.1 nm diameter and a nucleus composed of 96 Au atoms (Au314 (SR) 96). The nanostructures were modeled as reported elsewhere [54]. Six gold nanoclusters (6 AuNCs) were placed as a thin layer, emulating our previous experimental system (~ 20 nm gold nanoparticles) [24].

2.3 Molecular coupling

Adaptive Steered Molecular Dynamics (ASMD) simulations [55, 56] started from the earlier selected conformational states, derived from the GaMD analysis. Due to the chemical complexity of the enzyme surface, the simulations were driven by the electrostatic surface (PME) [49] for each conformational state and the 6 AuNCs layer. Hence, the scans were performed at speeds of 10 Å / ns, which comprised a 10-step profile of 1 Å between the Fe(III) of the active site and the ATP amino group or MBA group carboxylic of the 6 AuNCs. Three independent replicas were performed through gradual scans, involving measurements of the free energy of coupling between HRP and 6 AuNCs. The free energy was measured by shortening the distance between the central Fe(III) atom and the amino groups of the linkers, i.e., ATP amino groups and central Fe(III) atom of HRP heme (NH3+➔ Fe3+).


3. Results and discussion

3.1 HRP model system

The direct applications for protein structures produced by molecular modeling techniques such as homology modeling, include identification of structural and functional regions within a protein, and can be exploted as a lead for further experimental studies such as mutation analysis, catalytic and electrochemical characterization [27]. Furthermore, if homology modeling is combined with other computational methods such as molecular dynamics or quantum mechanics, the produced model can be used to screen different applications.

In our case, the available information comes from the genetic sequence of HRP enzyme and incomplete structures deposited on PDB, hence, it is necessary to reconstruct and evaluate our final model system based on homology models.

Our initial model was reconstructed with i-tasser metaserver, using as template the resting state horseradish peroxidase (1H5A), and building a waterbox for further molecular dynamics analysis (Figure 4).

Figure 4.

Horseradish peroxidase homology model reconstruction and MD system. (A) Modeling pipeline applied to reconstruction of HRP, (B) final solvated HRP model for MD simulations.

3.2 Conformational sampling

The homology model system of HRP was subjected to an initial conformational search with GaMD method for 500 ns of simulation which is comparable to conventional MD simulations on the order of microseconds. This data allow us to reconstruct an energy hypersurface from three structural and energetic variables: total energy of protein, RMSD of backbone and superficial hydrogen bonds count. The total energy of the proteins indicates how close is our model system to a minimum energy state; RMSD of the alpha carbons is related to changes in relative positions of the atoms, regarding to the X-ray crystal model (Figure 5); and the hydrogen bonds count between aspartic acid (D), glutamic acid (E), lysine (K) and arginine (R) residues, indicate how prone the enzyme is to bond to AuNCs linkers functional groups (NH3+ or COOH). With this approach, energetically more favorable conformations were selected under two criteria: (1) Low RMSD values, i.e., preserved HRP structure necessary for catalysis; (2) High solvent exposure of groups necessary for the esterification reaction, i.e., selecting those conformations where the probability of having the previous residues exposed to the solvent is higher. At the end of this sampling, three HRP conformations with the necessary intermolecular interactions were selected; with this approach we try to theoretically predict the formation of the amide bond in the solvent exposed residues, while the over-all enzyme structure remains correctly folded (Figure 5).

Figure 5.

Conformational sampling of HRP. The hypersurface was reconstructed using the second order McLauren cumulative expansion method. The free energy values were extracted from each conformation, depending on its structural variation and the solvent-accessible hydrogen bonds of amino acids D, E, K, and R.

Using the GaMD approach, we were able to explore the conformational diversity along different reaction coordinates of the HRP enzyme, and to predict more efficient couplings to AuNCs. This approximation over the boost potentials allows a reconstruction with less noise on the potential energy, which results in a more robust method to deal with artifacts during variations in the potential energy [57, 58]. Hence, our data suggest that during simulation, there are different enzyme conformations which preferentially bind to AuNCs and increase the interactions with the linker functional groups.

3.3 Interaction between modified 6 AuNCs and HRP

The main objective of this work is to give insights using molecular modeling methodologies, for the design of more efficient BEI, based on HRP enzyme and AuNCs. The structural description of HRP with molecular dynamics methodologies can help to design experimental strategies to increase the success of the immobilization and simultaneously preserve structure and reactivity of the enzyme.

In order to estimate the immobilization efficiency between GaMD selected structures of HRP and AuCNs functionalized with linkers, free energy calculations were performed with adaptive steered molecular dynamics. The prediction model system is intended to recreate the interface previously described in Section 1.2 (Figure 1). GC was used as current collector and HRP was coupled to functionalized-gold nano particles modified with 4-aminothiophenol (Figure 6A) [24].

Figure 6.

Gold nanoparticles experimental and theoretical model systems. (A) SEM image of gold nanoparticle electrodeposit (average particle diameter ≈ 18 nm) in glassy carbon electrode (GC), (B) AuNCs modified with 4-aminothiophenol linker used for free energy calculations.

For all the calculations, the average electrostatic surface potential was directly extracted from GaMD, since electrostatic term is computed for each integration step. HRP showed two well characterized electrostatic isosurfaces, with a prominent negatively charged surface derived from the presence of 12 exposed acidic residues. Each electrostatic potential was plotted as isosurface over each HRP conformation.

The free energy calculations between negatively charged isosurface of HRP and AuNCs modified with ATP linker (positively charged NH3+) showed coupling values of 60 Kcal · mol−1, with a minimum distance of 22 Å from the Fe(III) of the heme group (~ 2 Å to the AuNCs NH3+ functional group) (Figure 7A). On the other hand, the positively charged isosurface of HRP showed higher coupling energies of ~120 Kcal · mol−1 at a minimum distance of ~21.5 Å to the same AuNCs functional group. For the AuNCs modified with MBA, the coupling energy of the negatively charged isosurface of HRP was 130 Kcal · mol−1 at 33 Å distance, (~3 Å from surface functional groups), while the positive HRP surface showed free coupling energies of ~75 Kcal · mol−1 and minimum distances of 32.5 Å from Fe(III) heme group (3.5 Å from surface functional group) (Figure 7B).

Figure 7.

Coupling free energy profiles of HRP isosurfaces and AuNCs. (A) Free energy profile of HRP coupled to 6 AuNCs-ATP, (B) molecular interaction model of negative isosurface of HRP and 6 AuNCs-ATP, (C) free energy profile of HRP coupled to 6 AuNCs-MBA, (D) positive isosurface of HRP and 6 AuNCs-MBA.

The difference between both HRP isosurfaces with the positively charged surface imposed by ATP was ~70 Kcal · mol−1 which means that the enzyme electrostatic potential imposed a clear effect for the coupling to AuNCs. However, an overall difference of 15 Kcal · mol−1 between both linkers showed that the coupling assays were energetically more stable for the negative isosurfaces of the HRP enzyme than the positively charged.

The free energy calculations suggested that HRP showed an energy minimum around the bond distance between the carboxyl groups and the amino groups [59, 60]. Our data suggest that the intermolecular interactions guided by negative electrostatic surfaces of HRP, are of lower free energy for the positive charged AuNCs-ATP, than the negative charged AuNCs-MBA (Figure 8). The high electrostatic attraction between HRP and AuNCs-ATP would promote an efficient electrochemical response. On the other hand, the aforementioned interactions between HRP and negative charged linkers, like MBA, are less stable and result in higher free energy coupling profiles and lower electrostatic attraction for this interfaces.

Figure 8.

Different interactions displayed between HRP and AuNCs-linkers scaffolds.


4. Conclusions

The findings of this study on HRP coupled to gold nanoclusters, indicate that the polarized electrostatic isosurface potentials are key factors to select the most efficient linker for coupling. The resulted interaction energy and distance between HRP and AuNCs-ATP are adequate to promote the formation of covalent bonds between acidic residues and amino functional groups. The evidence from this study points towards the idea that molecular simulation methods, such as homology modeling and molecular dynamics are valuable tools to take into account for design of BEI. Our previous multidisciplinary work of VP6 capsids has led us to conclude that molecular dynamics simulations elucidate structural determinants to understand the behavior of biomolecules on BEI. An implication for using homology models coupled to molecular dynamics, is the possibility of widely sample the conformational space of biomolecules probes before electrochemical experimentation. Further theoretical and experimental studies are necessary to describe the interaction with other functional group linkers and validate by electrochemical techniques the real effect of the charge difference between AuNCs-ATP and AuNCs-MBA on the redox response of this enzyme, respectively.



The authors acknowledge the financial support of CONACYT A1-S-15877, PAPIIT-DGAPA-UNAM-IN204218 and Centro Nacional de Supercomputo-IPICyT, A.C. (National Supercomputing Center-CNS) grants for facilities at Thubat-Kaal 2 machine; very special thanks to Engineer Jesus Alaniz (CNS-IPICyT A.C) and Engineer Rogelio Elvira Morán (IER-UNAM) for their technical assistance (supercomputing and SEM analysis, respectively). W.I. García-García and G.A. Huerta-Miranda thank to CONCAYT for the PhD scholarship. A. Vidal-Limon thanks to CONACYT A1-S-15877 for the financial support of his academic fellowship.


  1. 1. Calabrese Barton S, Gallaway J, Atanassov P. Enzymatic Biofuel Cells for Implantable and Microscale Devices. Chem Rev. 2004;104(10):4867-86. DOI:10.1021/cr020719k
  2. 2. Cracknell JA, Vincent KA, Armstrong FA. Enzymes as Working or Inspirational Electrocatalysts for Fuel Cells and Electrolysis. Chem Rev. 2008;108(7):2439-61. DOI:10.1021/cr0680639
  3. 3. Szaciłowski K. Digital Information Processing in Molecular Systems. Chem Rev 2008;108(9):3481-548. DOI:10.1021/cr068403q
  4. 4. Jia X, Dong S, Wang E. Engineering the bioelectrochemical interface using functional nanomaterials and microchip technique toward sensitive and portable electrochemical biosensors. Biosens & Bioelectron. 2016;76:80-90. DOI: 10.1016/j.bios.2015.05.037
  5. 5. Zhou M, Dong S. Bioelectrochemical Interface Engineering: Toward the Fabrication of Electrochemical Biosensors, Biofuel Cells, and Self-Powered Logic Biosensors. Acc Chem Res . 2011;44(11):1232-43. DOI:10.1021/ar200096g
  6. 6. Nunes G., Marty JL. (2006) Immobilization of Enzymes on Electrodes. In: Guisan J.M., editors. Immobilization of Enzymes and Cells. Methods in Biotechnology™, vol 22. Humana Press. p. 239-50. DOI:10.1007/978-1-59745-053-9_21
  7. 7. Hwang ET, Lee S. Multienzymatic Cascade Reactions via Enzyme Complex by Immobilization. ACS Catalysis. 2019; (9):4402-25. DOI: 10.1021/acscatal.8b04921
  8. 8. Ganesh V, Pandey RR, Malhotra BD, Lakshminarayanan V. Electrochemical characterization of self-assembled monolayers (SAMs) of thiophenol and aminothiophenols on polycrystalline Au: Effects of potential cycling and mixed SAM formation. J Electroanal Chem. 2008;619-620(1-2):87-97. DOI: 10.1016/j.jelechem.2008.03.015
  9. 9. Łuczak T, Osińska M. New self-assembled layers composed with gold nanoparticles, cysteamine and dihydrolipoic acid deposited on bare gold template for highly sensitive and selective simultaneous sensing of dopamine in the presence of interfering ascorbic and uric acids. J Solid State Electrochem. 2017;21(3):747-58. DOI: 10.1007/s10008-016-3416-z
  10. 10. Abad JM, Vélez M, Santamaría C, Guisán JM, Matheus PR, Vázquez L, et al. Immobilization of Peroxidase Glycoprotein on Gold Electrodes Modified with Mixed Epoxy-Boronic Acid Monolayers. J Am Chem Soc. 2002;124(43):12845-53. DOI: 10.1021/ja026658p
  11. 11. Grabowska I, Maes W, Huynh Ngo T, Rohand T, Dehaen W, Radecki J, et al. Multiple redox-active sites in copper dipyrromethene-corrole self-assembled monolayers deposited onto gold electrodes. Int J Electrochem Sci. 2014;9(3):1232-49
  12. 12. Karimi Shervedani R, Samiei Foroushani M, Bagheri Dehaghi S. Functionalization of gold mercaptopropionic acid self-assembled monolayer with 5-amino-1,10-phenanthroline: Interaction with iron(II) and application for selective recognition of guanine. Electrochim Acta. 2015;164:344-52. DOI: 10.1016/j.electacta.2015.02.170
  13. 13. Mohamad NR, Marzuki NHC, Buang NA, Huyop F, Wahab RA. An overview of technologies for immobilization of enzymes and surface analysis techniques for immobilized enzymes. Biotechnol Biotechnol Equip. 2015;29(2):205-20. DOI: 10.1080/13102818.2015.1008192
  14. 14. Pinyou P, Blay V, Muresan LM, Noguer T. Enzyme-modified electrodes for biosensors and biofuel cells. Mater Horizons. 2019;1336-58. DOI: 10.1039/c9mh00013e
  15. 15. Carrara S, Baj-Rossi C, Boero C, De Micheli G. Do carbon nanotubes contribute to electrochemical biosensing? Electrochim Acta. 2014;128:102-12. DOI: 10.1016/j.electacta.2013.12.123
  16. 16. Adam C, Scodeller P, Grattieri M, Villalba M, Calvo EJ. Revisiting direct electron transfer in nanostructured carbon laccase oxygen cathodes. Bioelectrochemistry. 2016;109:101-7. DOI: 10.1016/j.bioelechem.2016.01.007
  17. 17. Chen A, Chatterjee S. Nanomaterials based electrochemical sensors for biomedical applications. Chem Soc Rev. 2013;42(12):5425-38. DOI: 10.1039/c3cs35518g
  18. 18. Konrad MP, Doherty AP, Bell SEJ. Stable and Uniform SERS Signals from Self-Assembled Two-Dimensional Interfacial Arrays of Optically Coupled Ag Nanoparticles. Anal Chem. 2013;85(14):6783-9. DOI: 10.1021/ac4008607
  19. 19. Alver E, Metin AÜ. Chitosan based metal-chelated copolymer nanoparticles: Laccase immobilization and phenol degradation studies. Int Biodeterior Biodegradation. 2017;125:235-42. DOI: 10.1016/j.ibiod.2017.07.012
  20. 20. Gao C, Guo Z, Liu J-HH, Huang X-JJ. The new age of carbon nanotubes: An updated review of functionalized carbon nanotubes in electrochemical sensors. Nanoscale. 2012;4(6):1948. DOI: 10.1039/c2nr11757f
  21. 21. Rosca V, Catalin Popescu I. Kinetic analysis of horseradish peroxidase “wiring” in redox polyelectrolyte-peroxidase multilayer assemblies. Electrochem commun. 2002;4(11):904-11. DOI: 10.1016/S1388-2481(02)00486-1
  22. 22. Raghu P, Madhusudana Reddy T, Reddaiah K, Jaidev LR, Narasimha G. A novel electrochemical biosensor based on horseradish peroxidase immobilized on Ag-nanoparticles/poly(l-arginine) modified carbon paste electrode toward the determination of pyrogallol/hydroquinone. Enzyme Microb Technol. 2013;52(6-7):377-85. DOI: 10.1016/j.enzmictec.2013.02.010
  23. 23. Rawal R, Chawla S, Pundir CS. Polyphenol biosensor based on laccase immobilized onto silver nanoparticles/multiwalled carbon nanotube/polyaniline gold electrode. Anal Biochem. 2011;419(2):196-204. DOI: 10.1016/j.ab.2011.07.028
  24. 24. Huerta-Miranda GA, Arrocha-Arcos AA, Miranda-Hernández M. Gold nanoparticles/4-aminothiophenol interfaces for direct electron transfer of horseradish peroxidase: Enzymatic orientation and modulation of sensitivity towards hydrogen peroxide detection. Bioelectrochemistry . 2018;122:77-83. DOI: 10.1016/j.bioelechem.2018.03.004
  25. 25. Smith AT, Veitch NC. Substrate binding and catalysis in heme peroxidases. Curr Opin Chem Biol. 1998;2(2):269-78. DOI: 10.1016/S1367-5931(98)80069-0
  26. 26. Battistuzzi G, Bellei M, Bortolotti CA, Sola M. Redox properties of heme peroxidases. Arch Biochem Biophys. 2010;500(1):21-36. DOI: 10.1016/
  27. 27. Fang Y-S, Huang X-J, Wang L-S, Wang J-F. An enhanced sensitive electrochemical immunosensor based on efficient encapsulation of enzyme in silica matrix for the detection of human immunodeficiency virus p24. Biosens Bioelectron. 2015;64:324-32. DOI: 10.1016/j.bios.2014.09.022
  28. 28. Su X-D, Zhang H, Terwilliger TC, Liljas A, Xiao J, Dong Y. Protein Crystallography from the Perspective of Technology Developments. Crystallogr Rev . 2015;21(1-2):122-53. DOI: 10.1080/0889311X.2014.973868
  29. 29. Helliwell JR. New developments in crystallography: exploring its technology, methods and scope in the molecular biosciences. Biosci Rep . 2017;37(4):BSR20170204. DOI: 10.1042/BSR20170204
  30. 30. Shigematsu H. Electron cryo-microscopy for elucidating the dynamic nature of live-protein complexes. Biochim Biophys acta Gen Subj. 2020;1864(2):129436. DOI: 10.1016/j.bbagen.2019.129436
  31. 31. Bai X, McMullan G, Scheres SHW. How cryo-EM is revolutionizing structural biology. Trends Biochem Sci. 2015;40(1):49-57. DOI: 10.1016/j.tibs.2014.10.005
  32. 32. Amunts A, Brown A, Toots J, Scheres SHW, Ramakrishnan V. Ribosome. The structure of the human mitochondrial ribosome. Science. 2015;348(6230):95-8. DOI: 10.1126/science.aaa1193
  33. 33. Tuszynski JA, Winter P, White D, Tseng CY, Sahu KK, Gentile F, et al. Mathematical and computational modeling in biology at multiple scales Biophysics. Theor Biol Med Model. 2014; 11(1):52. DOI: 10.1186/1742-4682-11-52
  34. 34. Dror RO, Dirks RM, Grossman JP, Xu H, Shaw DE. Biomolecular Simulation: A Computational Microscope for Molecular Biology. Annu Rev Biophys. 2012;41(1):429-52. DOI: 10.1146/annurev-biophys-042910-155245
  35. 35. Lu D, Aksimentiev A, Shih AY, Cruz-Chu E, Freddolino PL, Arkhipov A, et al. The role of molecular modeling in bionanotechnology. Physical Biology. 2006;3(1):S40-S53. DOI: 10.1088/1478-3975/3/1/S05
  36. 36. Ramachandran KI, Deepa G, Namboori K. Computational chemistry and molecular modeling: Principles and applications. 1st ed. Springer; 2008. 397 p. DOI: 10.1007/978-3-540-77304-7
  37. 37. Whitehead TA, Baker D, Fleishman SJ. Computational design of novel protein binders and experimental affinity maturation. In: Methods in Enzymology. Elsevier. 2013;523:1-19.DOI: 10.1016/B978-0-12-394292-0.00001-1
  38. 38. Heinecke A, Eckhardt W, Horsch M, Bungartz HJ. Molecular dynamics simulation. In: SpringerBriefs in Computer Science. Springer. 2015;(523): 11-29
  39. 39. Van Der Kamp MW, Mulholland AJ. Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology. Biochemistry. 2013; 52(16):2708-2728. DOI: 10.1021/bi400215w
  40. 40. van der Kamp MW, Mulholland AJ. Computational enzymology: insight into biological catalysts from modelling. Nat Prod Rep . 2008;25(6):1001-14. DOI: 10.1039/B600517A
  41. 41. Garcia-Garcia WI, Vidal-Limon A, Arrocha-Arcos AA, Palomares LA, Ramirez OT, Miranda-Hernández M. Rotavirus VP6 protein as a bio-electrochemical scaffold: Molecular dynamics and experimental electrochemistry. Bioelectrochemistry. 2019;127:180-186. DOI: 10.1016/j.bioelechem.2019.02.012
  42. 42. Plascencia-villa G, Mena JA, Castro-acosta RM, Fabián JC, Ramírez OT, Palomares LA. Strategies for the purification and characterization of protein scaffolds for the production of hybrid nanobiomaterials. 2011;879:1105-11. DOI: 10.1016/j.jchromb.2011.03.027
  43. 43. Carreño-Fuentes L, Plascencia-Villa G, Palomares LA, Moya SE, Ramírez OT. Modulating the physicochemical and structural properties of gold-functionalized protein nanotubes through thiol surface modification. Langmuir. 2014;30(49):14991-8. DOI: 10.1021/la503704a
  44. 44. Humphrey W, Dalke A, Schulten K. VMD: Visual molecular dynamics. J Mol Graph . 1996;14(1):33-8. DOI: 10.1016/0263-7855(96)00018-5
  45. 45. Berglund GI, Carlsson GH, Smith AT, Szöke H, Henriksen A, Hajdu J. The catalytic pathway of horseradish peroxidase at high resolution. Nature. 2002;417(6887):463-8. DOI: 10.1038/417463a
  46. 46. Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER Suite: protein structure and function prediction. Nat Methods . 2015;12(1):7-8. DOI: 10.1038/nmeth.3213
  47. 47. Olsson MHM, Søndergaard CR, Rostkowski M, Jensen JH. PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. J Chem Theory Comput. 2011;7(2):525-37. DOI: 10.1021/ct100578z
  48. 48. Maier JA, Martinez C, Kasavajhala K, Wickstrom L, Hauser KE, Simmerling C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J Chem Theory Comput . 2015;11(8):3696-713. DOI:10.1021/acs.jctc.5b00255
  49. 49. Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG. A smooth particle mesh Ewald method. J Chem Phys. 1995;103(19):8577-93. DOI: 10.1063/1.470117
  50. 50. Case DA, Walker RC, Cheatham TE, Simmerling C, Roitberg A, Merz KM, et al. Amber 18. Univ California, San Fr. 2018
  51. 51. Miao Y, McCammon JA. Gaussian Accelerated Molecular Dynamics: Theory, Implementation, and Applications. Annu Rep Comput Chem. 2017;13:231-78. DOI: 10.1016/bs.arcc.2017.06.005
  52. 52. Miao Y. Acceleration of biomolecular kinetics in Gaussian accelerated molecular dynamics. J Chem Phys. 2018; 149(7):072308. DOI: 10.1063/1.5024217
  53. 53. Miao Y, McCammon JA. Unconstrained enhanced sampling for free energy calculations of biomolecules: a review. Mol Simul. 2016; 42(13):1046-55. DOI: 10.1080/08927022.2015.1121541
  54. 54. Franco-Ulloa S, Riccardi L, Rimembrana F, Pini M, De Vivo M. NanoModeler: A Webserver for Molecular Simulations and Engineering of Nanoparticles. J Chem Theory Comput. 2019; 15(3):2022-32. DOI: 10.1021/acs.jctc.8b01304
  55. 55. Ozer G, Valeev EF, Quirk S, Hernandez R. Adaptive Steered Molecular Dynamics of the Long-Distance Unfolding of Neuropeptide Y. J Chem Theory Comput. 2010;6(10):3026-38. DOI:10.1021/ct100320g
  56. 56. Bureau HR, Merz DR, Hershkovits E, Quirk S, Hernandez R. Constrained unfolding of a helical peptide: Implicit versus explicit solvents. PLoS One. 2015;10(5):e0127034. DOI: 10.1371/journal.pone.0127034
  57. 57. Miao Y, Feher VA, McCammon JA. Gaussian Accelerated Molecular Dynamics: Unconstrained Enhanced Sampling and Free Energy Calculation. J Chem Theory Comput. 2015;11(8):3584-95. DOI: 10.1021/acs.jctc.5b00436
  58. 58. Miao Y, Feixas F, Eun C, McCammon JA. Accelerated molecular dynamics simulations of protein folding. J Comput Chem. 2015;36(20):1536-1549. DOI: 10.1002/jcc.23964
  59. 59. Valeur E, Bradley M. Amide bond formation: beyond the myth of coupling reagents. Chem Soc Rev. 2009;38(2):606-31. DOI: 10.1039/B701677H
  60. 60. Montalbetti CAGN, Falque V. Amide bond formation and peptide coupling. Tetrahedron . 2005;61(46):10827-52. DOI: 10.1016/j.tet.2005.08.031

Written By

Abraham Vidal-Limon, Guillermo Antonio Huerta-Miranda, Wendy I. García-García and Margarita Miranda-Hernández

Submitted: July 7th, 2020 Reviewed: September 4th, 2020 Published: September 25th, 2020