Topological and dynamical properties of Azurin anchored to a gold substrate as investigated by molecular dynamics simulation

Reorganization Free Energy of Copper Proteins in Solution, in Vacuum, and on Metal Surfaces

Metalloproteins, known to efficiently transfer electronic charge in biological systems, recently found their utilization in nanobiotechnological devices where the protein is placed into direct contact with metal surfaces. The feasibility of oxidation/reduction of the protein redox sites is affected by the reorganization free energies, one of the key parameters determining the transfer rates. While their values have been measured and computed for proteins in their native environments, i.e., in aqueous solution, the reorganization free energies of dry proteins or proteins adsorbed to metal surfaces remain unknown. Here, we investigate the redox properties of blue copper protein azurin, a prototypical redox-active metalloprotein previously probed by various experimental techniques both in solution and on metal/vacuum interfaces. We used a hybrid QM/MM computational technique based on DFT to explore protein dynamics, flexibility, and corresponding reorganization free energies in aqueous solution, vacuum, and on vacuum gold interfaces. Somewhat surprisingly, the reorganization free energy only slightly decreases when azurin is dried because the loss of the hydration shell leads to larger flexibility of the protein near its redox site. At the vacuum gold surfaces, the energetics of the structure relaxation depends on the adsorption geometry, however, significant reduction of the reorganization free energy was not observed. These findings have important consequences for the charge transport mechanism in vacuum devices, showing that the free energy barriers for protein oxidation remain significant even under ultra-high vacuum conditions.

Tuning Structure and Dynamics of Blue Copper Azurin Junctions via Single Amino-Acid Mutations

In the growing field of biomolecular electronics, blue-copper Azurin stands out as one of the most widely studied protein in single-molecule contacts. Interestingly, despite the paramount importance of the structure/dynamics of molecular contacts in their transport properties, these factors remain largely unexplored from the theoretical point of view in the context of single Azurin junctions. Here we address this issue using all-atom Molecular Dynamics (MD) of Pseudomonas Aeruginosa Azurin adsorbed to a Au(111) substrate. In particular, we focus on the structure and dynamics of the free/adsorbed protein and how these properties are altered upon single-point mutations. The results revealed that wild-type Azurin adsorbs on Au(111) along two well defined configurations: one tethered via cysteine groups and the other via the hydrophobic pocket surrounding the Cu 2 + . Surprisingly, our simulations revealed that single amino-acid mutations gave rise to a quenching of protein vibrations ultimately resulting in its overall stiffening. Given the role of amino-acid vibrations and reorientation in the dehydration process at the protein-water-substrate interface, we suggest that this might have an effect on the adsorption process of the mutant, giving rise to new adsorption configurations.

Experimental characterization and simulation of amino acid and peptide interactions with inorganic materials

Inspired by nature, many applications and new materials benefit from the interplay of inorganic materials and biomolecules. A fundamental understanding of complex organic-inorganic interactions would improve the controlled production of nanomaterials and biosensors to the development of biocompatible implants for the human body. Although widely exploited in applications, the interaction of amino acids and peptides with most inorganic surfaces is not fully understood. To date, precisely characterizing complex surfaces of inorganic materials and analyzing surface-biomolecule interactions remain challenging both experimentally and computationally. This article reviews several approaches to characterizing biomolecule-surface interactions and illustrates the advantages and disadvantages of the methods presented. First, we explain how the adsorption mechanism of amino acids/peptides to inorganic surfaces can be determined and how thermodynamic and kinetic process constants can be obtained. Second, we demonstrate how this data can be used to develop models for peptide-surface interactions. The understanding and simulation of such interactions constitute a basis for developing molecules with high affinity binding domains in proteins for bioprocess engineering and future biomedical technologies.

Structure, Dynamics, and Electron Transfer of Azurin Bound to a Gold Electrode

Blue copper redox protein azurin (AZ) constitutes an ideal active element for building bionano-optoelectronic devices based on the intriguing interplay among its electron transfer (ET), vibrational, and optical properties. A full comprehension of its dynamical and functional behavior is required for efficient applications. Here, AZ bound to gold electrode via its disulfide bridge was investigated by a molecular dynamics simulation approach taking into account for gold electron polarization which provides a more realistic description of the protein-gold interaction. Upon binding to gold, AZ undergoes slight changes in its secondary structure with the preservation of the copper-containing active site structure. Binding of AZ to gold promotes new collective motions, with respect to free AZ, as evidenced by essential dynamics. Analysis of the ET from the AZ copper ion to the gold substrate, performed by the Pathways model, put into evidence the main residues and structural motifs of AZ involved in the ET paths. During the dynamical evolution of the bionanosystem, transient contacts between some lateral protein atoms and the gold substrate occurred; concomitantly, the opening of additional ET channels with much higher rates was registered. These results provide new and detailed insights on the dynamics and ET properties of the AZ-gold system, by also helping to rationalize some imaging and conductive experimental evidences and also to design new bionanodevices with tailored features.

Modeling and simulation of protein-surface interactions: achievements and challenges

Understanding protein-inorganic surface interactions is central to the rational design of new tools in biomaterial sciences, nanobiotechnology and nanomedicine. Although a significant amount of experimental research on protein adsorption onto solid substrates has been reported, many aspects of the recognition and interaction mechanisms of biomolecules and inorganic surfaces are still unclear. Theoretical modeling and simulations provide complementary approaches for experimental studies, and they have been applied for exploring protein-surface binding mechanisms, the determinants of binding specificity towards different surfaces, as well as the thermodynamics and kinetics of adsorption. Although the general computational approaches employed to study the dynamics of proteins and materials are similar, the models and force-fields (FFs) used for describing the physical properties and interactions of material surfaces and biological molecules differ. In particular, FF and water models designed for use in biomolecular simulations are often not directly transferable to surface simulations and vice versa. The adsorption events span a wide range of time- and length-scales that vary from nanoseconds to days, and from nanometers to micrometers, respectively, rendering the use of multi-scale approaches unavoidable. Further, changes in the atomic structure of material surfaces that can lead to surface reconstruction, and in the structure of proteins that can result in complete denaturation of the adsorbed molecules, can create many intermediate structural and energetic states that complicate sampling. In this review, we address the challenges posed to theoretical and computational methods in achieving accurate descriptions of the physical, chemical and mechanical properties of protein-surface systems. In this context, we discuss the applicability of different modeling and simulation techniques ranging from quantum mechanics through all-atom molecular mechanics to coarse-grained approaches. We examine uses of different sampling methods, as well as free energy calculations. Furthermore, we review computational studies of protein-surface interactions and discuss the successes and limitations of current approaches.

A combined atomic force microscopy imaging and docking study to investigate the complex between p53 DNA binding domain and Azurin

The tumor suppressor p53 interacts with the redox copper protein Azurin (AZ) forming a complex which is of some relevance in biomedicine and cancer therapy. To obtain information on the spatial organization of this complex when it is immobilized on a substrate, we have used tapping mode-atomic force microscopy (TM-AFM) imaging combined with computational docking. The vertical dimension and the bearing volume of the DNA binding domain (DBD) of p53, anchored to functionalized gold substrate through exposed lysine residues, alone and after deposing AZ, have been measured by TM-AFM. By a computational docking approach, a three-dimensional model for the DBD of p53, before and after addition of AZ, have been predicted. Then we have calculated the possible arrangements of these biomolecular systems on gold substrate by finding a good agreement with the related experimental distribution of the height. The potentiality of the approach combining TM-AFM imaging and computational docking for the study of biomolecular complexes immobilized on substrates is briefly discussed.

Modeling the interaction between the N-terminal domain of the tumor suppressor p53 and azurin

It is known that the half life of the tumor suppressor p53 can be increased by the interaction with the bacterial protein azurin, resulting in an enhanced anti-tumoral activity. The understanding of the molecular mechanisms on the basis of this phenomenon can open the way to new anti-cancer strategies. Some experimental works have given evidence of an interaction between p53 and azurin (AZ); however the binding regions of the proteins are still unknown. Recently, fluorescence studies have shown that p53 partakes in the binding with the bacterial protein by its N-terminal (NT) domain. Here we have used a computational method to get insight into this interacting mode. The model that we propose for the best complex between AZ and p53 has been obtained from a rigid-body docking, coupled with a molecular dynamics (MD) simulation, a free energy calculation, and validated by mutagenesis analysis. We have found a high degree of geometric fit between the two proteins that are kept together by several hydrophobic interactions and numerous hydrogen bonds. Interestingly, it has emerged that AZ binds essentially to the helices H(I) and H(III) of the p53 NT domain, which are also interacting regions for the foremost inhibitor of p53, MDM2.

Molecular simulation of bovine beta-lactoglobulin adsorbed onto a positively charged solid surface

To obtain detailed insight into the mechanism of beta-lactoglobulin (beta-Lg) adsorption to a stainless steel surface at acidic pH, the adsorption of positively charged beta-Lg to a positively charged surface (Au (100) surface with virtual positive charge) was simulated using classical molecular dynamics. The initial orientation and position of beta-Lg on the surface were determined using Monte Carlo simulation using the implicit water system. Molecular dynamics simulation with the explicit water system was conducted for a 5 ns simulation time to monitor beta-Lg adsorption. To investigate surface charge density effects on adsorption of beta-Lg, the positive charge number per Au atom on the (100) surface, C, was varied from 0 to +0.0250|e|. Stable adsorption occurred in MD simulations when C was equal to or less than +0.0200|e|. Among these surface Au charge conditions, no large difference was observed in the vertical separation distance between the surface and the protein's center of mass, and the orientation angle. This fact indicates that the main interactions contributing to the adsorption were van der Waals interactions. The protein domain contacting the surface was near Thr125, agreeing with previous experimental studies. Considering simulation results and those previous experimental studies suggests a detailed adsorption mechanism of beta-Lg at acidic pH: beta-Lg molecule is adsorbed initially with the specific part of 125-135th residues close to the surface by van der Waals interactions. Simultaneously or subsequently, side carboxylic groups of acidic amino acid residues near the surface in 125-135th residues dissociate, leading to firmer adsorption by attractive electrostatic residue-surface interaction.

Docking study and free energy simulation of the complex between p53 DNA-binding domain and azurin

Molecular interaction between p53 tumor suppressor and the copper protein azurin (AZ) has been demonstrated to enhance p53 stability and hence antitumoral function, opening new perspectives in cancer treatment. While some experimental work has provided evidence for AZ binding to p53, no crystal structure for the p53-AZ complex was solved thus far. In this work the association between AZ and the p53 DNA-binding domain (DBD) was investigated by computational methods. Using a combination of rigid-body protein docking, experimental mutagenesis information, and cluster analysis 10 main p53 DBD-AZ binding modes were generated. The resulting structures were further characterized by molecular dynamics (MD) simulations and free energy calculations. We found that the highest scored docking conformation for the p53 DBD-AZ complex also yielded the most favorable free energy value. This best three-dimensional model for the complex was validated by using a computational mutagenesis strategy. In this structure AZ binds to the flexible L(1) and s(7)-s(8) loops of the p53 DBD and stabilizes them through protein-protein tight packing interactions, resulting in high degree of both surface matching and electrostatic complementarity.

Electronic coupling between azurin and gold at different protein/substrate orientations

By means of constrained classical molecular dynamics simulations, we have computed the structure of azurin deposited on a Au(111) surface at different possible orientations and the azimuthal forces acting on the protein at each sampled conformation. We have then evaluated the effect of the angular variation on the speed of electron tunneling between the protein redox site and the metal surface. We find that the azurin/gold electronic coupling has a strong dependence on the molecular orientation and is greatly enhanced by inclining the protein to lie as flat as possible on the surface. We discuss the implications of our results for scanning probe microscopy experiments in which tunneling currents are measured while the protein is subjected to mechanical forces exerted by the tip of the instrument.