Structural and dynamical properties of manganese catalase and the synthetic protein DF1 and their implication for reactivity from classical molecular dynamics calculations

Comparison of Experimental and Broken Symmetry Density Functional Theory Calculated Electron Paramagnetic Resonance Parameters for the Manganese Catalase Active Site in the Superoxidized MnIII/MnIV State

Broken symmetry density functional theory has been used to calculate g-tensor, ⁵⁵Mn, ¹⁴N, and ¹⁷O hyperfine couplings for active site models of superoxidized Mn^III/Mn^IV manganese catalase both in its native and azide-inhibited form. While a good agreement is found between the calculated and experimental g-tensor and ⁵⁵Mn hyperfine couplings for all models, the active site geometry and Mn ion oxidation state can only be readily distinguished based on a comparison of the calculated and experimental ¹⁴N azide and ¹⁷O HFCs. This comparison shows that only models containing a Jahn-Teller distorted 5-coordinate (Mn^III)₂ site and a 6-coordinate (Mn^IV)₁ site can satisfactorily reproduce the experimental ¹⁴N and ¹⁷O hyperfine couplings.

Biochemical and spectroscopic characterization of dinuclear Mn-sites in artificial four-helix bundle proteins

To better understand metalloproteins with Mn-clusters, we have designed artificial four-helix bundles to have one, two, or three dinuclear metal centers able to bind Mn(II). Circular dichroism measurements showed that the Mn-proteins have substantial α-helix content, and analysis of electron paramagnetic resonance spectra is consistent with the designed number of bound Mn-clusters. The Mn-proteins were shown to catalyze the conversion of hydrogen peroxide into molecular oxygen. The loss of hydrogen peroxide was dependent upon the concentration of protein with bound Mn, with the proteins containing multiple Mn-clusters showing greater activity. Using an oxygen sensor, the oxygen concentration was found to increase with a rate up to 0.4μM/min, which was dependent upon the concentrations of hydrogen peroxide and the Mn-protein. In addition, the Mn-proteins were shown to serve as electron donors to bacterial reaction centers using optical spectroscopy. Similar binding of the Mn-proteins to reaction centers was observed with an average dissociation constant of 2.3μM. The Mn-proteins with three metal centers were more effective at this electron transfer reaction than the Mn-proteins with one or two metal centers. Thus, multiple Mn-clusters can be incorporated into four-helix bundles with the capability of performing catalysis and electron transfer to a natural protein.

Conserved methionine dictates substrate preference in Nramp-family divalent metal transporters

Natural resistance-associated macrophage protein (Nramp) family transporters catalyze uptake of essential divalent transition metals like iron and manganese. To discriminate against abundant competitors, the Nramp metal-binding site should favor softer transition metals, which interact either covalently or ionically with coordinating molecules, over hard calcium and magnesium, which interact mainly ionically. The metal-binding site contains an unusual, but conserved, methionine, and its sulfur coordinates transition metal substrates, suggesting a vital role in their transport. Using a bacterial Nramp model system, we show that, surprisingly, this conserved methionine is dispensable for transport of the physiological manganese substrate and similar divalents iron and cobalt, with several small amino acid replacements still enabling robust uptake. Moreover, the methionine sulfur's presence makes the toxic metal cadmium a preferred substrate. However, a methionine-to-alanine substitution enables transport of calcium and magnesium. Thus, the putative evolutionary pressure to maintain the Nramp metal-binding methionine likely exists because it-more effectively than any other amino acid-increases selectivity for low-abundance transition metal transport in the presence of high-abundance divalents like calcium and magnesium.

Manganese Oxidation State Assignment for Manganese Catalase

The oxidation state assignment of the manganese ions present in the superoxidized manganese (III/IV) catalase active site is determined by comparing experimental and broken symmetry density functional theory calculated (14)N, (17)O, and (1)H hyperfine couplings. Experimental results have been interpreted to indicate that the substrate water is coordinated to the Mn(III) ion. However, by calculating hyperfine couplings for both scenarios we show that water is coordinated to the Mn(IV) ion and that the assigned oxidation states of the two manganese ions present in the site are the opposite of that previously proposed based on experimental measurements alone.

Non-heme manganese catalase--the 'other' catalase

Non-heme manganese catalases are widely distributed over microbial life and represent an environmentally important alternative to heme-containing catalases in antioxidant defense. Manganese catalases contain a binuclear manganese complex as their catalytic active site rather than a heme, and cycle between Mn(2)(II,II) and Mn(2)(III,III) states during turnover. X-ray crystallography has revealed the key structural elements of the binuclear manganese active site complex that can serve as the starting point for computational studies on the protein. Four manganese catalase enzymes have been isolated and characterized, and the enzyme appears to have a broad phylogenetic distribution including both bacteria and archae. More than 100 manganese catalase genes have been annotated in genomic databases, although the assignment of many of these putative manganese catalases needs to be experimentally verified. Iron limitation, exposure to low levels of peroxide stress, thermostability and cyanide resistance may provide the biological and environmental context for the occurrence of manganese catalases.

Structure and dynamics of a dizinc metalloprotein: effect of charge transfer and polarization

Structures and dynamics of a recently designed dizinc metalloprotein (DFsc) (J. Mol. Biol. 2003, 334, 1101) are studied by molecular dynamics simulation using a dynamically adapted polarized force field derived from fragment quantum calculation for protein in solvent. To properly describe the effect of charge transfer and polarization in the present approach, quantum chemistry calculation of the zinc-binding group is periodically performed (on-the-fly) to update the atomic charges of the zinc-binding group during the MD simulation. Comparison of the present result with those obtained from simulations under standard AMBER force field reveals that charge transfer and polarization are critical to maintaining the correct asymmetric metal coordination in the DFsc. Detailed analysis of the result also shows that dynamic fluctuation of the zinc-binding group facilitates solvent interaction with the zinc ions. In particular, the dynamic fluctuation of the zinc-zinc distance is shown to be an important feature of the catalytic function of the di-ion zinc-binding group. Our study demonstrates that the dynamically adapted polarization approach is computationally practical and can be used to study other metalloprotein systems.

Dissection of the difference between the group I metal ions in inhibiting GSK3β: a computational study

Glycogen synthase kinase 3β (GSK3β) is a serine/threonine kinase that requires two cofactor Mg(2+) ions for catalysis in regulating many important cellular signals. Experimentally, Li(+) is a competitive inhibitor of GSK3β relative to Mg(2+), while this mechanism is not experienced with other group I metal ions. Herein, we use native Mg(2)(2+)-Mg(1)(2+) GSK3β and its Mg(2)(2+)-M(1)(+) (M = Li, Na, K, and Rb) derivatives to investigate the effect of metal ion substitution on the mechanism of inhibition through two-layer ONIOM-based quantum mechanics/molecular mechanics (QM/MM) calculations and molecular dynamics (MD) simulations. The results of ONIOM calculations elucidate that the interaction of Na(+), K(+), and Rb(+) with ATP is weaker compared to that of Mg(2+) and Li(+) with ATP, and the critical triphosphate moiety of ATP undergoes a large conformational change in the Na(+), K(+), and Rb(+) substituted systems. As a result, the three metal ions (Na(+), K(+), and Rb(+)) are not stable and depart from the active site, while Mg(2+) and Li(+) can stabilize in the active site, evident in MD simulations. Comparisons of Mg(2)(2+)-Mg(1)(2+) and Mg(2)(2+)-Li(1)(+) systems reveal that the inline phosphor-transfer of ATP and the two conserved hydrogen bonds between Lys85 and ATP, together with the electrostatic potential at the Li(1)(+) site, are disrupted in the Mg(2)(2+)-Li(1)(+) system. These computational results highlight the possible mechanism why Li(+) inhibits GSK3β.

Computational design of protein-ligand interfaces: potential in therapeutic development

Computational design of protein-ligand interfaces finds optimal amino acid sequences within a small-molecule binding site of a protein for tight binding of a specific small molecule. It requires a search algorithm that can rapidly sample the vast sequence and conformational space, and a scoring function that can identify low energy designs. This review focuses on recent advances in computational design methods and their application to protein-small molecule binding sites. Strategies for increasing affinity, altering specificity, creating broad-spectrum binding, and building novel enzymes from scratch are described. Future prospects for applications in drug development are discussed, including limitations that will need to be overcome to achieve computational design of protein therapeutics with novel modes of action.

Controlling and fine tuning the physical properties of two identical metal coordination sites in de novo designed three stranded coiled coil peptides

Herein we report how de novo designed peptides can be used to investigate whether the position of a metal site along a linear sequence that folds into a three-stranded α-helical coiled coil defines the physical properties of Cd(II) ions in either CdS(3) or CdS(3)O (O-being an exogenous water molecule) coordination environments. Peptides are presented that bind Cd(II) into two identical coordination sites that are located at different topological positions at the interior of these constructs. The peptide GRANDL16PenL19IL23PenL26I binds two Cd(II) as trigonal planar 3-coordinate CdS(3) structures whereas GRANDL12AL16CL26AL30C sequesters two Cd(II) as pseudotetrahedral 4-coordinate CdS(3)O structures. We demonstrate how for the first peptide, having a more rigid structure, the location of the identical binding sites along the linear sequence does not affect the physical properties of the two bound Cd(II). However, the sites are not completely independent as Cd(II) bound to one of the sites ((113)Cd NMR chemical shift of 681 ppm) is perturbed by the metalation state (apo or [Cd(pep)(Hpep)(2)](+) or [Cd(pep)(3)](-)) of the second center ((113)Cd NMR chemical shift of 686 ppm). GRANDL12AL16CL26AL30C shows a completely different behavior. The physical properties of the two bound Cd(II) ions indeed depend on the position of the metal center, having pK(a2) values for the equilibrium [Cd(pep)(Hpep)(2)](+) → [Cd(pep)(3)](-) + 2H(+) (corresponding to deprotonation and coordination of cysteine thiols) that range from 9.9 to 13.9. In addition, the L26AL30C site shows dynamic behavior, which is not observed for the L12AL16C site. These results indicate that for these systems one cannot simply assign a "4-coordinate structure" and assume certain physical properties for that site since important factors such as packing of the adjacent Leu, size of the intended cavity (endo vs exo) and location of the metal site play crucial roles in determining the final properties of the bound Cd(II).

Challenges in the computational design of proteins

Protein design has many applications not only in biotechnology but also in basic science. It uses our current knowledge in structural biology to predict, by computer simulations, an amino acid sequence that would produce a protein with targeted properties. As in other examples of synthetic biology, this approach allows the testing of many hypotheses in biology. The recent development of automated computational methods to design proteins has enabled proteins to be designed that are very different from any known ones. Moreover, some of those methods mostly rely on a physical description of atomic interactions, which allows the designed sequences not to be biased towards known proteins. In this paper, we will describe the use of energy functions in computational protein design, the use of atomic models to evaluate the free energy in the unfolded and folded states, the exploration and optimization of amino acid sequences, the problem of negative design and the design of biomolecular function. We will also consider its use together with the experimental techniques such as directed evolution. We will end by discussing the challenges ahead in computational protein design and some of their future applications.

Evolution of catalases from bacteria to humans

Excessive hydrogen peroxide is harmful for almost all cell components, so its rapid and efficient removal is of essential importance for aerobically living organisms. Conversely, hydrogen peroxide acts as a second messenger in signal-transduction pathways. H(2)O(2) is degraded by peroxidases and catalases, the latter being able both to reduce H(2)O(2) to water and to oxidize it to molecular oxygen. Nature has evolved three protein families that are able to catalyze this dismutation at reasonable rates. Two of the protein families are heme enzymes: typical catalases and catalase-peroxidases. Typical catalases comprise the most abundant group found in Eubacteria, Archaeabacteria, Protista, Fungi, Plantae, and Animalia, whereas catalase-peroxidases are not found in plants and animals and exhibit both catalatic and peroxidatic activities. The third group is a minor bacterial protein family with a dimanganese active site called manganese catalases. Although catalyzing the same reaction (2 H(2)O(2)--> 2 H(2)O+ O(2)), the three groups differ significantly in their overall and active-site architecture and the mechanism of reaction. Here, we present an overview of the distribution, phylogeny, structure, and function of these enzymes. Additionally, we report about their physiologic role, response to oxidative stress, and about diseases related to catalase deficiency in humans.

Progress in computational protein design

Current progress in computational structure-based protein design is reviewed in the areas of methodology and applications. Foundational advances include new potential functions, more efficient ways of computing energetics, flexible treatments of solvent, and useful energy function approximations, as well as ensemble-based approaches to scoring designs for inclusion of entropic effects, improvements to guaranteed and to stochastic search techniques, and methods to design combinatorial libraries for screening and selection. Applications include new approaches and successes in the design of specificity for protein folding, binding, and catalysis, in the redesign of proteins for enhanced binding affinity, and in the application of design technology to study and alter enzyme catalysis. Computational protein design continues to mature and advance.

Single Crystal 55Mn ENDOR of Concanavalin A: Detection of Two Mn2+ Sites with Different 55Mn Quadrupole Tensors

Concanavalin A is a member of the plant hemeagglutinin (or plant lectin) family that contains two metal binding sites; one, called S1, is occupied by Mn2+ and the other, S2, by Ca2+. 55Mn electron-nuclear double resonance (ENDOR) measurements were performed on a single crystal of concanavalin A at W-band (95 GHz, ~3.5 T) to determine the 55Mn nuclear quadrupole interaction in a protein binding site and its relation to structural parameters. Such measurements are easier at a high field because of the high sensitivity for size-limited samples and the reduction of second-order effects on the spectrum which simplifies spectral analysis. The analysis of the 55Mn ENDOR rotation patterns showed that two chemically inequivalent Mn2+ types are present at low temperatures, although the high-resolution X-ray structure reported only one site. Their quadrupole coupling constants, e2Qq/h, are significantly different; 10.7 +/- 0.6 MHz for Mand only -2.7 +/-0.6 MHz for M. The ENDOR data also refined the hyperfine coupling determined earlier by single-crystal EPR measurements, yielding a small but significant difference between the two: -262.5 MHz for M and -263.5 MHz for M. The principal z-axis for M is not aligned with any of the Mn-ligand directions, but is 25 off the Mn-asp10 direction, and its orientation is different than that of the zero-field splitting (ZFS) interaction. Because of the small quadrupole interaction of M the orientation dependence was very mild, leading to larger uncertainties in the asymmetry parameter. Nonetheless, there too z is not along the Mn-ligand bonds and is rotated 90 with respect to MnA. These results show, that similar to the ZFS, the quadrupolar interaction is highly sensitive to small differences in the coordination sphere of the Mn2+, and the resolution of the two types is in agreement with the earlier observation of a two-site conformational dynamic detected through the ZFS interaction, which is frozen out at low temperatures and averaged at room temperature. To account for the structural origin of the different e2Qq/h values, the electric field gradient tensor was calculated using the point-charge model. The calculations showed that a relatively small displacement of the oxygen ligand of asp10 can lead to differences on the order observed experimentally.

Modeling the charge distribution at metal sites in proteins for molecular dynamics simulations

Almost half of the proteome of living organisms is constituted of metalloproteins. Unfortunately, the ability of the current generation of molecular dynamics pairwise-additive forcefields to properly describe metal pockets is severely lacking due to the intrinsic difficulty of handling polarization and charge transfer contributions. In order to improve the description of metalloproteins, a simple reparameterization strategy is proposed herein that does not involve artificial constraints. Specifically, a non-bonded quantum mechanical-based model is used to capture the mean polarization and charge transfer contributions to the interatomic forces within the metal site. The present approach is demonstrated to provide enough accuracy to maintain the integrity of the metal pocket for a variety of metalloproteins during extended (multi-nanosecond) molecular dynamics simulations. The method enables the sampling of small conformational changes and the relaxation of local frustrations in NMR structures.