A new molecular mechanics force field for the oxidized form of blue copper proteins

HNO to NO Conversion Mechanism with Copper Zinc Superoxide Dismutase, Comparison with Heme Protein Mediated Conversions, and the Origin of Questionable Reversibility

The interconversion of NO and HNO, via copper zinc superoxide dismutase (CuZnSOD), is important in biomedicine and for HNO detection. Many mechanistic questions, including the decades-long debate on reversibility, were resolved in this work. Calculations of various active-site and full-protein models show that the basic mechanism is proton-coupled electron transfer with a computed barrier of 10.98 kcal mol^-1 , which is in excellent agreement with experimental results (10.62 kcal mol^-1 ), and this nonheme protein-mediated reaction has many significant mechanistic differences compared with the conversions mediated by heme proteins due to geometric and electronic factors. The reasons for the irreversible nature of this conversion and models with the first thermodynamically favorable and kinetically feasible mechanism for the experimental reverse reaction were discovered. Such results are the first for nonheme enzyme mediated HNO to NO conversions, which shall facilitate other related studies and HNO probe development.

Systematic Parameterization and Simulation of Boronic Acid-β-Lactamase Aqueous Solution in Developing the ABEEMσπ Polarizable Force Field

Boronic acid, an inhibitor of β-lactamase, has begun to be applied to the treatment of biological infections and tumors. Scientists are working to develop new and more effective boronic acid. Molecular dynamics (MD) simulation provides a powerful auxiliary tool for drug design. However, the current force fields have no boron-related parameters. In this work, an atom-bond electronegativity equalization method at the σπ level (ABEEMσπ) polarizable force field (ABEEMσπ PFF) of boronic acid and β-lactamase has been developed to determine the potential functions and parameters. The interaction between boron and serine in β-lactamase is regarded as a bonded mode. The interaction between them is simulated by the Morse potential energy function, which is close to the experimental change of the stretching potential energy in a large range. The potential energy surfaces of the bond length, bond angle, and dihedral angle of boronic acid-β-lactamase have the same stability point and change trend as M06-2X/6-311G**. For 47 boronic acid-β-lactamase training molecules, the linear correlation coefficient (R) of the charge distribution between the ABEEMσπ PFF and HF/STO-3G is greater than 0.96. Attributed to the fact that the charge distribution of the ABEEMσπ PFF can fluctuate with the change of geometry and environment, the polarization effect and charge-transfer effect are well reflected. The binding ability of different boronic acids with the same β-lactamase is different. A total of 10 boronic acid-β-lactamase model molecules and 10 boronic acid-β-lactamase and water complexes are simulated. The order of binding energy of five large model molecules calculated by the ABEEMσπ PFF is consistent with that of the MP2 method. The binding energies of boronic acid-β-lactamase and water complexes are close to those of the MP2 method. The results of MD simulation of five aqueous boronic acid-β-lactamase complexes in the NVT ensemble verify the rationality of boron-related parameters of the ABEEMσπ PFF, which have a good application prospect. This study lays a solid theoretical foundation for further study of the inhibition of boronic acid on β-lactamase.

Evolved Fusarium oxysporum laccase expressed in Saccharomyces cerevisiae

Fusarium oxysporum laccase was functionally expressed in Saccharomyces cerevisiae and engineered towards higher expression levels and higher reactivity towards 2,6-dimethoxyphenol, that could be used as a mediator for lignin modification. A combination of classical culture optimization and protein engineering led to around 30 times higher activity in the culture supernatant. The winner mutant exhibited three times lower Km, four times higher kcat and ten times higher catalytic efficiency than the parental enzyme. The strategy for laccase engineering was composed of a combination of random methods with a rational approach based on QM/MM MD studies of the enzyme complex with 2,6-dimethoxyphenol. Laccase mediator system with 2,6-dimethoxyphenol caused fulvic acids release from biosolubilized coal.

AMBER Force Field Parameters for Cobalt-Containing Biological Systems: A Systematic Derivation Study

In the present work, the parameterization of a set of cobalt-containing systems has been performed to create a comprehensive library for bonding parameters of biological Co-containing systems. A standard process for the extraction and validation of parameters was employed, which could be used to create force field parameters for the other metal-containing systems. All protein data banks were searched to extract common chemical groups in bonding with cobalt, and finally, 16 structures were designed to represent the binding model of the chemical moieties with cobalt. The Hessian matrix of each structure was computed at the B3LYP/6-311++G(2d,2p) level of theory and the Seminario method was employed to compute cobalt bond stretching and angle bending parameters. Validation of the derived parameters was performed using structural minimization and molecular dynamics (MD) simulations of four models. Further validation was performed using an extensive MD simulation on carbonic anhydrase II as a common cobalt-containing metalloprotein. The results demonstrated that among models, the bonded model in combination with the RESP charges can produce the most reliable and accurate structural conformations for the metal site of cobalt-containing systems.

Transition Metal Ion Interactions with Disordered Amyloid-β Peptides in the Pathogenesis of Alzheimer's Disease: Insights from Computational Chemistry Studies

Monomers and oligomers of the amyloid-β peptide aggregate to form the fibrils found in the brains of Alzheimer's disease patients. These monomers and oligomers are largely disordered and can interact with transition metal ions, affecting the mechanism and kinetics of amyloid-β aggregation. Due to the disordered nature of amyloid-β, its rapid aggregation, as well as solvent and paramagnetic effects, experimental studies face challenges in the characterization of transition metal ions bound to amyloid-β monomers and oligomers. The details of the coordination chemistry between transition metals and amyloid-β obtained from experiments remain debated. Furthermore, the impact of transition metal ion binding on the monomeric or oligomeric amyloid-β structures and dynamics are still poorly understood. Computational chemistry studies can serve as an important complement to experimental studies and can provide additional knowledge on the binding between amyloid-β and transition metal ions. Many research groups conducted first-principles calculations, ab initio molecular dynamics simulations, quantum mechanics/classical mechanics simulations, and classical molecular dynamics simulations for studying the interplay between transition metal ions and amyloid-β monomers and oligomers. This review summarizes the current understanding of transition metal interactions with amyloid-β obtained from computational chemistry studies. We also emphasize the current view of the coordination chemistry between transition metal ions and amyloid-β. This information represents an important foundation for future metal ion chelator and drug design studies aiming to combat Alzheimer's disease.

Key Role of the Local Hydrophobicity in the East Patch of Plastocyanins on Their Thermal Stability and Redox Properties

Understanding the molecular basis of the thermal stability and functionality of redox proteins has important practical applications. Here, we show a distinct thermal dependence of the spectroscopic and electrochemical properties of two plastocyanins from the thermophilic cyanobacterium Phormidium laminosum and their mesophilic counterpart from Synechocystis sp. PCC 6803, despite the similarity of their molecular structures. To explore the origin of these differences, we have mimicked the local hydrophobicity in the east patch of the thermophilic protein by replacing a valine of the mesophilic plastocyanin by isoleucine. Interestingly, the resulting mutant approaches the thermal stability, redox thermodynamics, and dynamic coupling of the flexible site motions of the thermophilic protein, indicating the existence of a close connection between the hydrophobic packing of the east patch region of plastocyanin and the functional control and stability of the oxidized and reduced forms of the protein.

Copper Binding Induces Polymorphism in Amyloid-β Peptide: Results of Computational Models

Amyloid-β (Aβ) peptides are intrinsically disordered peptides, and their aggregation is the hallmark of Alzheimer's disease development. The propensity of the Aβ peptide to intermolecular interactions, the latter favoring different types of oligomers and aggregated forms, has been the object of a huge number of studies. Several facts are now established: the presence of large amount of d-block (M) ions (Zn, Cu, and Fe) in the aggregated forms; the 1:1 M/Aβ ratio favors the formation of amorphous aggregates, with an aggregation rate lower than that in the absence of such ions. In particular, statistical models describing the interactions between copper and amyloid peptides are mandatory to explain the relationship between neurodegeneration, copper dyshomeostasis, and overproduction of reactive oxygen species, the latter event occurring with aging. In this work, we show, by replica-exchange molecular dynamics simulations, that a copper ion (Cu²⁺) bound as in the experimentally observed prevailing coordination enhances the probability of closed structures that hinder the formation of extended intermolecular hydrogen bonds that stabilize fibrillar ordered aggregated forms. On the other hand, this effect enhances the catalytic role of the complex during the lifetime of soluble forms.

Dinuclear ZnII and mixed CuII-ZnII complexes of artificial patellamides as phosphatase models

The patellamides (cyclic pseudo-octapeptides) are produced by Prochloron, a symbiont of the ascidians, marine invertebrate filter feeders. These pseudo-octapeptides are present in the cytoplasm and a possible natural function of putative metal complexes of these compounds is hydrolase activity, however the true biological role is still unknown. The dinuclear Cu^II complexes of synthetic patellamide derivatives have been shown in in vitro experiments to be efficient hydrolase model catalysts. Many hydrolase enzymes, specifically phosphatases and carboanhydrases, are Zn^II-based enzymes and therefore, we have studied the Zn^II and mixed Zn^II/Cu^II solution chemistry of a series of synthetic patellamide derivatives, including solution structural and computational work, with the special focus on model phosphatase chemistry with bis-(2,4-dinitrophenyl)phosphate (BDNPP) as the substrate. The Zn^II complexes of a series of ligands are shown to form complexes of similar structure and stability compared to the well-studied Cu^II analogues and the phosphatase reactivities are also similar. Since the complex stabilities and phosphatase activities are generally a little lower compared to those of Cu^II and since the concentration of Zn^II in Prochloron cells is slightly smaller, we conclude that the Cu^II complexes of the patellamides are more likely to be of biological importance.

Impact of Cu(II) Binding on Structures and Dynamics of Aβ42 Monomer and Dimer: Molecular Dynamics Study

The classical force field, which is compatible with the Amber force field 99SB, has been obtained for the interaction of Cu(II) with monomer and dimers of amyloid-β peptides using the coordination where Cu(II) is bound to His6, His13 (or His14), and Asp1 with distorted planar geometry. The newly developed force field and molecular dynamics simulation were employed to study the impact of Cu(II) binding on structures and dynamics of Aβ₄₂ monomer and dimers. It was shown that in the presence of Cu(II) the β content of monomer is reduced substantially compared with the wild-type Aβ₄₂ suggesting that, in accord with experiments, metal ions facilitate formation of amorphous aggregates rather than amyloid fibrils with cross-β structures. In addition, one possible mechanism for amorphous assembly is that the Asp23-Lys28 salt bridge, which plays a crucial role in β sheet formation, becomes more flexible upon copper ion binding to the Aβ N-terminus. The simulation of dimers was conducted with the Cu(II)/Aβ stoichiometric ratios of 1:1 and 1:2. For the 1:1 ratio Cu(II) delays the Aβ dimerization process as observed in a number of experiments. The mechanism underlying this phenomenon is associated with slow formation of interchain salt bridges in dimer as well as with decreased hydrophobicity of monomer upon Cu-binding.

Water exit pathways and proton pumping mechanism in B-type cytochrome c oxidase from molecular dynamics simulations

Cytochrome c oxidase (CcO) is a vital enzyme that catalyzes the reduction of molecular oxygen to water and pumps protons across mitochondrial and bacterial membranes. While proton uptake channels as well as water exit channels have been identified for A-type CcOs, the means by which water and protons exit B-type CcOs remain unclear. In this work, we investigate potential mechanisms for proton transport above the dinuclear center (DNC) in ba3-type CcO of Thermus thermophilus. Using long-time scale, all-atom molecular dynamics (MD) simulations for several relevant protonation states, we identify a potential mechanism for proton transport that involves propionate A of the active site heme a3 and residues Asp372, His376 and Glu126(II), with residue His376 acting as the proton-loading site. The proposed proton transport process involves a rotation of residue His376 and is in line with experimental findings. We also demonstrate how the strength of the salt bridge between residues Arg225 and Asp287 depends on the protonation state and that this salt bridge is unlikely to act as a simple electrostatic gate that prevents proton backflow. We identify two water exit pathways that connect the water pool above the DNC to the outer P-side of the membrane, which can potentially also act as proton exit transport pathways. Importantly, these water exit pathways can be blocked by narrowing the entrance channel between residues Gln151(II) and Arg449/Arg450 or by obstructing the entrance through a conformational change of residue Tyr136, respectively, both of which seem to be affected by protonation of residue His376.

Data for molecular dynamics simulations of B-type cytochrome c oxidase with the Amber force field

Cytochrome c oxidase (CcO) is a vital enzyme that catalyzes the reduction of molecular oxygen to water and pumps protons across mitochondrial and bacterial membranes. This article presents parameters for the cofactors of ba₃-type CcO that are compatible with the all-atom Amber ff12SB and ff14SB force fields. Specifically, parameters were developed for the Cu_A pair, heme b, and the dinuclear center that consists of heme a₃ and Cu_B bridged by a hydroperoxo group. The data includes geometries in XYZ coordinate format for cluster models that were employed to compute proton transfer energies and derive bond parameters and point charges for the force field using density functional theory. Also included are the final parameter files that can be employed with the Amber leap program to generate input files for molecular dynamics simulations with the Amber software package. Based on the high resolution (1.8 Å) X-ray crystal structure of the ba₃-type CcO from Thermus thermophilus (Protein Data Bank ID number PDB: 3S8F), we built a model that is embedded in a POPC lipid bilayer membrane and solvated with TIP3P water molecules and counterions. We provide PDB data files of the initial model and the equilibrated model that can be used for further studies.

Extending the essential dynamics analysis to investigate molecular properties: application to the redox potential of proteins

Here, a methodology is proposed to investigate the collective fluctuation modes of an arbitrary set of observables, maximally contributing to the fluctuation of another functionally relevant observable. The methodology, based on the analysis of fully classical molecular dynamics (MD) simulations, exploits the essential dynamics (ED) method, originally developed to analyse the collective motions in proteins. We apply this methodology to identify the residues that are more relevant for determining the reduction potential (E(0)) of a redox-active protein. To this aim, the fluctuation modes of the single-residue electrostatic potentials mostly contributing to the fluctuations of the total electrostatic potential (the main determinant of E(0)) are investigated for wild-type azurin and two of its mutants with a higher E(0). By comparing the results here obtained with a previous study on the same systems [Zanetti-Polzi et al., Org. Biomol. Chem., 2015, 13, 11003] we show that the proposed methodology is able to identify the key sites that determine E(0). This information can be used for a general deeper understanding of the molecular mechanisms on the basis of the redox properties of the proteins under investigation, as well as for the rational design of mutants with a higher or lower E(0). From the results of the present analysis we propose a new azurin mutant that, according to our calculations, shows a further increase of E(0).

A few key residues determine the high redox potential shift in azurin mutants

The wide range of variability of the reduction potential (E(0)) of blue-copper proteins has been the subject of a large number of studies in the past several years. In particular, a series of azurin mutants have been recently rationally designed tuning E(0) over a very broad range (700 mV) without significantly altering the redox-active site [Marshall et al., Nature, 2009, 462, 113]. This clearly suggests that interactions outside the primary coordination sphere are relevant to determine E(0) in cupredoxins. However, the molecular determinants of the redox potential variability are still undisclosed. Here, by means of atomistic molecular dynamics simulations and hybrid quantum/classical calculations, the mechanisms that determine the E(0) shift of two azurin mutants with high potential shifts are unravelled. The reduction potentials of native azurin and of the mutants are calculated obtaining results in good agreement with the experiments. The analysis of the simulations reveals that only a small number of residues (including non-mutated ones) are relevant in determining the experimentally observed E(0) variation via site-specific, but diverse, mechanisms. These findings open the path to the rational design of new azurin mutants with different E(0).

Modeling the met form of human tyrosinase: a refined and hydrated pocket for antagonist design

Tyrosinases are type 3 copper proteins involved in melanin biosynthesis, responsible for skin and hair color in mammals. To steer tyrosinase inhibitor discovery for therapeutic and cosmetic purposes, structural information about human tyrosinase is necessary. As this protein has never been crystallized so far, we derived a robust homology model built using structural information from Streptomyces castaneoglobisporus and Ipomea batata catecholoxidase enzymes. The active site containing two copper atoms in co-ordination with six histidine residues was refined through an optimization protocol based on molecular mechanics parameters for copper co-ordination and charges calculated by quantum mechanics methods. Five structural water molecules and a hydroxyl ion were found to be essential for optimization. The superimposition of the human homology model on crystallographic structures of tyrosinases from other species revealed similar overall backbone topologies, active site conformations, and conserved water molecules. Phenylthiourea (PTU), the tyrosinase inhibitor of reference, was then docked into the solvated human active pocket. A binding mode consistent with crystallographic information was obtained. Taken together, these findings demonstrated that the human tyrosinase model, deposited in the Protein Model Database, is a reliable structure for future rational inhibitor design projects.

Importance of polarization effect in the study of metalloproteins: application of polarized protein specific charge scheme in predicting the reduction potential of azurin

Molecular dynamics (MD) simulation is commonly used in the study of protein dynamics, and in recent years, the extension of MD simulation to the study of metalloproteins is gaining much interest. Choice of force field is crucial in MD studies, and the inclusion of metal centers complicates the process of accurately describing the electrostatic environment that surrounds the redox centre. Herein, we would like to explore the importance of including electrostatic contribution from both protein and solvent in the study of metalloproteins. MD simulations with the implementation of thermodynamic integration will be conducted to model the reduction process of azurin from Pseudomonas aeruginosa. Three charge schemes will be used to derive the partial charges of azurin. These charge schemes differ in terms of the amount of immediate environment, respective to copper, considered during charge fitting, which ranges from the inclusion of copper and residues in the first coordination sphere during density functional theory charge fitting to the comprehensive inclusion of protein and solvent effect surrounding the metal centre using polarized protein-specific charge scheme. From the simulations conducted, the relative reduction potential of the mutated azurins respective to that of wild-type azurin (ΔEcal) were calculated and compared with experimental values. The ΔEcal approached experimental value with increasing consideration of environmental effect hence substantiating the importance of polarization effect in the study of metalloproteins. This study also attests the practicality of polarized protein-specific charge as a computational tool capable of incorporating both protein environment and solvent effect into MD simulations.

An Angular Overlap Model for Cu(II) Ion in the AMOEBA Polarizable Force Field

An extensible polarizable force field for transition metal ion was developed based on AMOEBA and the angular overlap model (AOM) with consistent treatment of electrostatics for all atoms. Parameters were obtained by fitting molecular mechanics (MM) energies to various ab initio gas-phase calculations. The results of parameterization were presented for copper (II) ion ligated to water and model fragments of amino acid residues involved in the copper binding sites of type 1 copper proteins. Molecular dynamics (MD) simulations were performed on aqueous copper (II) ion at various temperatures, as well as plastocyanin (1AG6) and azurin (1DYZ). Results demonstrated that the AMOEBA-AOM significantly improves the accuracy of classical MM in a number of test cases when compared to ab initio calculations. The Jahn-Teller distortion for hexa-aqua copper (II) complex was handled automatically without specifically designating axial and in-plane ligands. Analyses of MD trajectories resulted in a 6-coordination first solvation shell for aqueous copper (II) ion and a 1.8ns average residence time of water molecules. The ensemble average geometries of 1AG6 and 1DYZ copper binding sites were in general agreement with X-ray and previous computational studies.

Unsymmetrical binding modes of the HOPNO inhibitor of tyrosinase: from model complexes to the enzyme

The deciphering of the binding mode of tyrosinase (Ty) inhibitors is essential to understand how to regulate the tyrosinase activity. In this paper, by combining experimental and theoretical methods, we studied an unsymmetrical tyrosinase functional model and its interaction with 2-hydroxypyridine-N-oxide (HOPNO), a new and efficient competitive inhibitor for bacterial Ty. The tyrosinase model was a dinuclear copper complex bridged by a chelated ring with two different complexing arms (namely (bis(2-ethylpyridyl)amino)methyl and (bis(2-methylpyridyl)amino)methyl). The geometrical asymmetry of the complex induces an unsymmetrical binding of HOPNO. Comparisons have been made with the binding modes obtained on similar symmetrical complexes. Finally, by using quantum mechanics/molecular mechanics (QM/MM) calculations, we studied the binding mode in tyrosinase from a bacterial source. A new unsymmetrical binding mode was obtained, which was linked to the second coordination sphere of the enzyme.

A molecular simulation based assessment of binding of metal ions on micelles

An assessment of the ability of a micellar surface to bind different metal ions using molecular simulation is presented in this study. Sodium dodecyl sulfate (SDS) is considered as the anionic surfactant. Various relevant characteristics of SDS-metal ion systems are estimated to quantify preferential binding of metal ions. These are electrostatic energy, total potential energy of the system, radial distribution function, and entropy and free energy change of the system. By examining these parameters, the relative extents of binding of different metal ions to the micellar surface are assessed.

Structural Survey of Zinc Containing Proteins and the Development of the Zinc AMBER Force Field (ZAFF)

Currently the Protein Data Bank (PDB) contains over 18,000 structures that contain a metal ion including Na, Mg, K, Ca, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Pd, Ag, Cd, Ir, Pt, Au, and Hg. In general, carrying out classical molecular dynamics (MD) simulations of metalloproteins is a convoluted and time consuming process. Herein, we describe MCPB (Metal Center Parameter Builder), which allows one, to conveniently and rapidly incorporate metal ions using the bonded plus electrostatics model (Hoops et al., J. Am. Chem. Soc. 1991, 113, 8262-8270) into the AMBER Force Field (FF). MCPB was used to develop a Zinc FF, ZAFF, which is compatible with the existing AMBER FFs. The PDB was mined for all Zn containing structures with most being tetrahedrally bound. The most abundant primary shell ligand combinations were extracted and FFs were created. These include Zn bound to CCCC, CCCH, CCHH, CHHH, HHHH, HHHO, HHOO, HOOO, HHHD, and HHDD (O = water and the remaining are 1 letter amino acid codes). Bond and angle force constants and RESP charges were obtained from B3LYP/6-31G* calculations of model structures from the various primary shell combinations. MCPB and ZAFF can be used to create FFs for MD simulations of metalloproteins to study enzyme catalysis, drug design and metalloprotein crystal refinement.

Surface residues dynamically organize water bridges to enhance electron transfer between proteins

Cellular energy production depends on electron transfer (ET) between proteins. In this theoretical study, we investigate the impact of structural and conformational variations on the electronic coupling between the redox proteins methylamine dehydrogenase and amicyanin from Paracoccus denitrificans. We used molecular dynamics simulations to generate configurations over a duration of 40 ns (sampled at 100-fs intervals) in conjunction with an ET pathway analysis to estimate the ET coupling strength of each configuration. In the wild-type complex, we find that the most frequently occurring molecular configurations afford superior electronic coupling due to the consistent presence of a water molecule hydrogen-bonded between the donor and acceptor sites. We attribute the persistence of this water bridge to a "molecular breakwater" composed of several hydrophobic residues surrounding the acceptor site. The breakwater supports the function of nearby solvent-organizing residues by limiting the exchange of water molecules between the sterically constrained ET region and the more turbulent surrounding bulk. When the breakwater is affected by a mutation, bulk solvent molecules disrupt the water bridge, resulting in reduced electronic coupling that is consistent with recent experimental findings. Our analysis suggests that, in addition to enabling the association and docking of the proteins, surface residues stabilize and control interprotein solvent dynamics in a concerted way.

Changes in non-core regions stabilise plastocyanin from the thermophilic cyanobacterium Phormidium laminosum

We report a theoretical investigation on the different stabilities of two plastocyanins. The first one belongs to the thermophilic cyanobacterium Phormidium laminosum and the second one belongs to its mesophilic relative Synechocystis sp. These proteins share the same topology and secondary-structure elements; however, the melting temperatures of their oxidised species differ by approximately 15 K. Long-time-scale molecular dynamics simulations, performed at different temperatures, show that the thermophilic protein optimises a set of intramolecular interactions (interstrand hydrogen bonding, salt bridging and hydrophobic clustering) within the region that comprises the strands beta 5 and beta 6, loop L5 and the helix. This region exhibits most of the differences in the primary sequence between the two proteins and, in addition, it is involved in the interaction with known physiological partners. Further work is in progress to unveil the specific structural features responsible for the different thermal stability of the two proteins.

Structural investigation of zymogenic and activated forms of human blood coagulation factor VIII: a computational molecular dynamics study

BACKGROUND

Human blood coagulation factor VIII (fVIII) is a large plasma glycoprotein with sequential domain arrangement in the order A1-a1-A2-a2-B-a3-A3-C1-C2. The A1, A2 and A3 domains are interconnected by long linker peptides (a1, a2 and a3) that possess the activation sites. Proteolysis of fVIII zymogen by thrombin or factor Xa results in the generation of the activated form (fVIIIa) which serves as a critical co-factor for factor IXa (fIXa) enzyme in the intrinsic coagulation pathway.

RESULTS

In our efforts to elucidate the structural differences between fVIII and fVIIIa, we developed the solution structural models of both forms, starting from an incomplete 3.7 A X-ray crystal structure of fVIII zymogen, using explicit solvent MD simulations. The full assembly of B-domainless single-chain fVIII was built between the A1-A2 (Ala1-Arg740) and A3-C1-C2 (Ser1669-Tyr2332) domains. The structural dynamics of fVIII and fVIIIa, simulated for over 70 ns of time scale, enabled us to evaluate the integral motions of the multi-domain assembly of the co-factor and the possible coordination pattern of the functionally important calcium and copper ion binding in the protein.

CONCLUSIONS

MD simulations predicted that the acidic linker peptide (a1) between the A1 and A2 domains is largely flexible and appears to mask the exposure of putative fIXa enzyme binding loop (Tyr555-Asp569) region in the A2 domain. The simulation of fVIIIa, generated from the zymogen structure, predicted that the linker peptide (a1) undergoes significant conformational reorganization upon activation by relocating completely to the A1-domain. The conformational transition led to the exposure of the Tyr555-Asp569 loop and the surrounding region in the A2 domain. While the proposed linker peptide conformation is predictive in nature and warrants further experimental validation, the observed conformational differences between the zymogen and activated forms may explain and support the large body of experimental data that implicated the critical importance of the cleavage of the peptide bond between the Arg372 and Ser373 residues for the full co-factor activity of fVIII.

Structural and mechanistic insights into the oxy form of tyrosinase from molecular dynamics simulations

The first, long time scale (16-ns) ligand field molecular dynamics (LFMD) simulations of the oxy form of tyrosinase are reported. The calculations use our existing type 3 copper force field for the peroxido-bridged [Cu(2)O(2)](2+) unit which is here translated from MMFF into the AMBER format together with a new charge scheme. The protein secondary and tertiary structures are not significantly altered by removing the 'caddie' protein, ORF378, which must be bound to tyrosinase before crystals will grow. A comprehensive principal component analysis of the Cartesian coordinates from the final 8 ns shows that the protein backbone is relatively rigid. However, the significant butterfly fold of the [Cu(2)O(2)](2+) moiety observed in the X-ray structure, presumably due to the caddie protein tyrosine at the active site, is absent in the simulations. LFMD gives a clear and persistent distinction between equatorial and axial Cu-N distances, with the latter about 0.2 A longer and remaining syn to each other. However, the two coordination spheres display important differences. LFMD simulations of the symmetric model complex [mu-eta(2):mu(2)-O(2){Cu(Meim)(3)}(2)](2+) (Meim is 5-methyl-1H-imidazole) provide a mechanism for syn-anti interchange of axial ligands which suggests, in combination with the old experimental X-ray data, the new LFMD simulations and traditional coordination chemistry arguments, that His(54) on Cu(A) is 'insipiently axial' and that a combination of a butterfly distortion of the [Cu(2)O(2)](2+) group and a rotation of the Cu(A)(His)(3) moiety converts the vacant, initially axial, binding site on Cu(A) into a much more favourable equatorial site.

Electrostatic strain and concerted motions in the transient complex between plastocyanin and cytochrome f from the cyanobacterium Phormidium laminosum

Many fleeting macromolecular interactions, like those being involved in electron transport, are essential in biology. However, little is known about the behaviour of the partners and their dynamics within their short-lived complex. To tackle such issue, we have performed molecular dynamics simulations on an electron transfer complex formed by plastocyanin and cytochrome f from the cyanobacterium Phormidium laminosum. Besides simulations of the isolated partners, two independent trajectories of the complex were calculated, starting from the two different conformations in the NMR ensemble. The first one leads to a more stable ensemble with a shorter distance between the metal sites of the two partners. The second experiences a significant drift of the complex conformation. Analyses of the distinct calculations show that the conformation of cytochrome f is strained upon binding of its partner, and relaxes upon its release. Interestingly, the principal component analysis of the trajectories indicates that plastocyanin displays a concerted motion with the small domain of cytochrome f that can be attributed to electrostatic interactions between the two proteins.

Long Distance Electron-Transfer Mechanism in Peptidylglycine α-Hydroxylating Monooxygenase: A Perfect Fitting for a Water Bridge

The active sites of copper enzymes have been the subject of many theoretical and experimental investigations from a number of years. Such studies have embraced topics devoted to the modeling of the first coordination sphere at the metallic cations up to the development of biomimetic, or bioinspired, catalytic systems. At least from the theoretical viewpoint, fewer efforts have been dedicated to elucidate how the two copper cations act concertedly in noncoupled dicopper enzymes such as peptidylglycine alpha-hydroxylating monooxygenase (PHM) and dopamine beta-monooxygenase (DbetaM). In these metalloenzymes, an electronic transfer is assumed between the two distant copper cations (11 A). Recent experimental results suggest that this transfer occurs through water molecules, a phenomenon which has been theoretically evidenced to be of high efficiency in the case of cytochrome b5 (Science, 2005, 310, 1311). In the present contribution dedicated to PHM, we overpass the common theoretical approaches dedicated to the electronic and geometrical structures of sites CuM or CuH restricted to their first coordination spheres and aim at directly comparing theoretical results to the experimentally measured activity of the PHM enzyme. To achieve this goal, molecular dynamics simulations were performed on wild-type and various mutants of PHM. More precisely, we provide an estimate of the electron-transfer efficiency between the CuM and CuH sites by means of such molecular dynamics simulations coupled to Marcus theory joined to the Beratan model to approximate the required coupling matrix elements. The theoretical results are compared to the kinetics measurements performed on wild and mutated PHM. The present work, the dynamic aspects of which are essential, accounts for the experimental results issued from mutagenesis. It supports the conclusion that an electronic transfer can occur between two copper(I) sites along a bridge involving a set of hydrogen and chemical bonds. Residue Gln170 is evidenced to be the keystone of this water-mediated pathway.

Molecular dynamics simulations of the enzyme Cu, Zn superoxide dismutase

The enzyme Cu, Zn superoxide dismutase (Cu,Zn-SOD) is a ubiquitous oxireductase, which is responsible for the cellular defense against oxidative stress caused by the high toxicity of the superoxide radical, and has been also linked to some cases of familiar amyotrophic lateral sclerosis. In the present study a set of molecular mechanics parameters for the active site of Cu,Zn-SOD has been derived. Afterward, an extensive molecular dynamics simulation has been carried out in an aqueous environment. The obtained results shed a further light on the structural flexibility of the backbone, where the active site is nested, and the solvation shell occupancy. The relatively small backbone deviation, shown by a root-mean-square deviation below 1.0 A, confirms the accuracy of the parameters. The solvent shell analysis has shown that the first solvation shell is located at about 5 A from the copper ion, generating an empty cavity with enough space to accommodate the superoxide radical. The low residence time means that a high permutation rate of water molecules in both solvation shells is consistent with the efficiency of this catalytic mechanism. Hybrid studies using ONIOM methodologies can now be done to evaluate the mechanistic implications of the explicit inclusion of the whole system.

Insights into Cu(I) exchange in HAH1 using quantum mechanical and molecular simulations

The human antioxidant protein, HAH1, is an important participant in a Cu(I) transport chain, delivering one Cu(I) ion to the Wilson's (WND) or Menkes disease protein (MNK). Full geometry optimizations and second-derivative calculations were performed on several binding site models using the B3LYP functional to derive parameters for the construction of a novel molecular mechanical (MM) force field for Cu(I) and its ligating residues in HAH1. MM minimization and molecular dynamics (MD) calculations were then performed using the AMBER suite to validate the newly generated force field. The X-ray crystal structure of the protein and the geometry of the Cu(I) binding site within the protein were reproduced by the MD simulations on the protein based on rmsd and visual inspection, validating the new force field parameters. The results from the quantum mechanical (QM) and MD simulations suggest that either a two- or three-coordinate exchange reaction is preferred and that it is unlikely that a four-coordinate Cu(I) species plays a role in copper exchange.

Comprehensive Molecular Mechanics Model for Oxidized Type I Copper Proteins: Active Site Structures, Strain Energies, and Entatic Bulging

The ligand field molecular mechanics (LFMM) model has been applied to the oxidized Type 1 copper center. In conjunction with the AMBER94 force field implemented in DommiMOE, the ligand field extension of the molecular operating environment (MOE), LFMM parameters for Cu-N(imidazole), Cu-S(thiolate), Cu-S(thioether), and Cu-O(carbonyl) interactions were developed on the basis of experimental and theoretical data for homoleptic model systems. Subsequent LFMM optimizations of the active site model complex [Cu(imidazole)2(SMe)(SMe2]+ agree with high level quantum results both structurally and energetically. Stable trigonal and tetragonal structures are located with the latter about 1.5 kcal mol-1 lower in energy. Fully optimized unconstrained structures were computed for 24 complete proteins containing T1 centers spanning four-coordinate, plastocyanin-like CuN2SS' and stellacyanin-like CuN2SO sites, plus the five-coordinate CuN2SS'O sites of the azurins. The initial structures were based on PDB coordinates augmented by a 10 A layer of water molecules. Agreement between theory and experiment is well within the experimental uncertainties. Moreover, the LFMM results for plastocyanin (Pc), cucumber basic protein (CBP) and azurin (Az) are at least as good as previously reported QM/MM structures and are achieved several orders of magnitude faster. The LFMM calculations suggest the protein provides an entatic strain of about 10 kcal mol-1. However, when combined with the intrinsic 'plasticity' of d9 Cu(II), different starting protein/solvent configurations can have a significant effect on the final optimized structure. This 'entatic bulging' results in relatively large fluctuations in the calculated metal-ligand bond lengths. For example, simply on the basis of 25 different starting configurations of the solvent molecules, the optimized Cu-S(thiolate) bond lengths in Pc vary by 0.04 A while the Cu-S(thioether) distance spans over 0.3 A. These variations are the same order of magnitude as the differences often quoted to correlate the spectroscopic properties from a set of proteins. Isolated optimizations starting from PDB coordinates (or indeed, the PDB structures themselves) may only accidentally correlate with spectroscopic measurements. The present calculations support the work of Warshel who contends that adequate configurational averaging is necessary to make proper contact with experimental properties measured in solution. The LFMM is both sufficiently accurate and fast to make this practical.

Mechanism of Nitrite Reduction at T2Cu Centers: Electronic Structure Calculations of Catalysis by Copper Nitrite Reductase and by Synthetic Model Compounds

The mechanism of nitrite reduction at the Cu(II) center of both copper nitrite reductase and a number of corresponding synthetic models has been investigated by using both QM/MM and cluster calculations employing density functional theory methods. The mechanism in both cases is found to be very similar. Initially nitrite is bound in a bidentate fashion to the Cu(II) center via the two oxygen atoms. Upon reduction of the copper center, the two possible coordination modes of the protonated nitrite, by either nitrogen or a single oxygen atom, are close in energy, with nitrogen coordination probably preferred. Further protonation of this species leads to N-O bond cleavage, and an electron transfer from the Cu(I) center to the N-O+ ligand, resulting in loss of NO and regeneration of the resting state of the enzyme having a bound water molecule.

Computer Simulation of the Interaction of Cu(I) with Cys Residues at the Binding Site of the Yeast Metallochaperone Cu(I)−Atx1

The copper binding site and electronic structure of the metallochaperone protein Atx1 were investigated using the combination of quantum mechanics methods and molecular mechanics methods in the ONIOM(QM:MM) scheme at the density functional theory (DFT) B3LYP/ 6-31G(d):AMBER level. The residues in the binding site, -Met13-Thr14-Cys15-Cu(I)-Cys18-Gly17-Ser16-, were modeled with QM and the rest of the residues with MM. Our results indicate that the structure for Cu(I)-Atx1 has the copper atom coordinated to two sulfur atoms from Cys15 (2.110 A) and Cys18 (2.141 A) with an angle S-Cu(I) -S of 166 degrees . The potential energy surface of the copper atom is used to estimate its binding energy and the force field for the copper ligands. The potential surface is shallow for the bending mode S-Cu-S, which explains the origin of the disorder observed in crystallographic and nuclear magnetic resonance studies. Using molecular dynamics for Cu(I)-Atx1 in a box of water molecules and in vacuum, with the force field derived in this work, we observed a correlated motion between the side chains of Thr14 and of Lys65 which enhances distortions in the S-Cu-S geometry. The results are compared with recent experiments and the previous models. The vibrational spectra for the copper ligands and for the residues in the binding site were computed. The localized modes for the copper ligands and the amide bands were assigned. The presence of the copper atom affects the amide bands' frequencies of the residues Cys15 and Cys18, giving resolved bands that can be used to sense changes in the binding site upon translocation of copper atom or interaction with target proteins. Furthermore, the EXAFS (extended X-ray absorption fine structure) spectrum of the proposed structure for Cu(I)-Atx1 was calculated and reproduced the experiments fairly well.

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

A classical molecular dynamics study of the electron transfer protein azurin, covalently bound to a gold substrate through its native disulphide group, is carried out at full hydration. With the aim of investigating the effects on the protein structure and dynamics as induced by the presence of an electric field, simulations are performed on neutral, positively and negatively charged substrates. A number of parameters, such as the average structure, the root mean square deviations and fluctuations, the intraprotein hydrogen bonds and solvent accessible surface of the protein, are monitored during 10 ns of run. The orientation, the height and the lateral size of the protein, with respect to the substrate are evaluated and compared with the experimental data obtained by scanning probe nanoscopies. The electron transfer properties between the copper redox center and the disulphide bridge bound to the substrate are investigated and briefly discussed.