Electronic and geometric structures of the blue copper site of azurin investigated by QM/MM hybrid calculations

Converged Structural and Spectroscopic Properties for Refined QM/MM Models of Azurin

Blue copper proteins continue to challenge experiment and theory with their electronic structure and spectroscopic properties that respond sensitively to the coordination environment of the copper ion. In this work, we report state-of-the art electronic structure studies for geometric and spectroscopic properties of the archetypal "Type I" copper protein azurin in its Cu(II) state. A hybrid quantum mechanics/molecular mechanics (QM/MM) approach is used, employing both density functional theory (DFT) and coupled cluster with singles, doubles, and perturbative triples (CCSD(T)) methods for the QM region, the latter method making use of the domain-based local pair natural orbital (DLPNO) approach. Models of increasing QM size are employed to investigate the convergence of critical geometric parameters. It is shown that convergence is slow and that a large QM region is critical for reproducing the short experimental Cu-SCys112 distance. The study of structural convergence is followed by investigation of spectroscopic parameters using both DFT and DLPNO-CC methods and comparing these to the experimental spectrum using simulations. The results allow us to examine for the first time the distribution of spin densities and hyperfine coupling constants at the coupled cluster level, leading us to revisit the experimental assignment of the 33S hyperfine splitting. The wavefunction-based approach to obtain spin-dependent properties of open-shell systems demonstrated here for the case of azurin is transferable and applicable to a large array of bioinorganic systems.

COMPARATIVE INVESTIGATION OF THE EFFECT OF TYPE OF DENSITY FUNCTIONAL IN THE DETERMINATION OF GEOMETRICAL PARAMETERS IN A Cu COMPLEX

The effects of short-range electron correlation, long-range electron exchange, local and nonlocal parts of density, higher order gradients of density, and adding some percentage of Hartree–Fock exchange to the functional on the prediction of geometrical parameters were investigated. A copper complex namely 1,2-bis(1,4,7-triaza-1-cyclononyl) ethane copper (II) with Jahn–Teller distortion in octahedral geometry was used to evaluate the performance of 50 commonly available density functionals. The standard 3-21G basis set was used for all light elements, while pseudo potential LANL2DZ was used for the copper atom. The best bond lengths and bond angles were obtained using M05-2x and OP functionals respectively. Also in order to more accurate survey the performance of B3LYP, we used this functional with two all-electron basis sets (6-31G and 3-21G) and three basis sets involving effective core potentials (LANL2DZ/3-21G, LANL2DZ, and LACVP).

Spectroscopic and DFT studies of second-sphere variants of the type 1 copper site in azurin: covalent and nonlocal electrostatic contributions to reduction potentials

The reduction potentials (E(0)) of type 1 (T1) or blue copper (BC) sites in proteins and enzymes with identical first coordination spheres around the redox active copper ion can vary by ~400 mV. Here, we use a combination of low-temperature electronic absorption and magnetic circular dichroism, electron paramagnetic resonance, resonance Raman, and S K-edge X-ray absorption spectroscopies to investigate a series of second-sphere variants--F114P, N47S, and F114N in Pseudomonas aeruginosa azurin--which modulate hydrogen bonding to and protein-derived dipoles nearby the Cu-S(Cys) bond. Density functional theory calculations correlated to the experimental data allow for the fractionation of the contributions to tuning E(0) into covalent and nonlocal electrostatic components. These are found to be significant, comparable in magnitude, and additive for active H-bonds, while passive H-bonds are mostly nonlocal electrostatic in nature. For dipoles, these terms can be additive to or oppose one another. This study provides a methodology for uncoupling covalency from nonlocal electrostatics, which, when coupled to X-ray crystallographic data, distinguishes specific local interactions from more long-range protein/active interactions, while affording further insight into the second-sphere mechanisms available to the protein to tune the E(0) of electron-transfer sites in biology.

Biological applications of hybrid quantum mechanics/molecular mechanics calculation

Since in most cases biological macromolecular systems including solvent water molecules are remarkably large, the computational costs of performing ab initio calculations for the entire structures are prohibitive. Accordingly, QM calculations that are jointed with MM calculations are crucial to evaluate the long-range electrostatic interactions, which significantly affect the electronic structures of biological macromolecules. A UNIX-shell-based interface program connecting the quantum mechanics (QMs) and molecular mechanics (MMs) calculation engines, GAMESS and AMBER, was developed in our lab. The system was applied to a metalloenzyme, azurin, and PU.1-DNA complex; thereby, the significance of the environmental effects on the electronic structures of the site of interest was elucidated. Subsequently, hybrid QM/MM molecular dynamics (MD) simulation using the calculation system was employed for investigation of mechanisms of hydrolysis (editing reaction) in leucyl-tRNA synthetase complexed with the misaminoacylated tRNA(Leu), and a novel mechanism of the enzymatic reaction was revealed. Thus, our interface program can play a critical role as a powerful tool for state-of-the-art sophisticated hybrid ab initio QM/MM MD simulations of large systems, such as biological macromolecules.

Modulation of electronic structures of bases through DNA recognition of protein

The effects of environmental structures on the electronic states of functional regions in a fully solvated DNA·protein complex were investigated using combined ab initio quantum mechanics/molecular mechanics calculations. A complex of a transcriptional factor, PU.1, and the target DNA was used for the calculations. The effects of solvent on the energies of molecular orbitals (MOs) of some DNA bases strongly correlate with the magnitude of masking of the DNA bases from the solvent by the protein. In the complex, PU.1 causes a variation in the magnitude among DNA bases by means of directly recognizing the DNA bases through hydrogen bonds and inducing structural changes of the DNA structure from the canonical one. Thus, the strong correlation found in this study is the first evidence showing the close quantitative relationship between recognition modes of DNA bases and the energy levels of the corresponding MOs. Thus, it has been revealed that the electronic state of each base is highly regulated and organized by the DNA recognition of the protein. Other biological macromolecular systems can be expected to also possess similar modulation mechanisms, suggesting that this finding provides a novel basis for the understanding for the regulation functions of biological macromolecular systems.