Functional Characterization of Two PLP-Dependent Enzymes Involved in Capsular Polysaccharide Biosynthesis from Campylobacter jejuni.

Campylobacter jejuni is a Gram-negative, pathogenic bacterium that causes campylobacteriosis, a form of gastroenteritis. C. jejuni is the most frequent cause of food-borne illness in the world, surpassing Salmonella and E. coli. Coating the surface of C. jejuni is a layer of sugar molecules known as the capsular polysaccharide that, in C. jejuni NCTC 11168, is composed of a repeating unit of d-glycero-l-gluco-heptose, d-glucuronic acid, d-N-acetyl-galactosamine, and d-ribose. The d-glucuronic acid moiety is further amidated with either serinol or ethanolamine. It is unknown how these modifications are synthesized and attached to the polysaccharide. Here, we report the catalytic activities of two previously uncharacterized, pyridoxal phosphate (PLP)-dependent enzymes, Cj1436 and Cj1437, from C. jejuni NCTC 11168. Using a combination of mass spectrometry and nuclear magnetic resonance, we determined that Cj1436 catalyzes the decarboxylation of l-serine phosphate to ethanolamine phosphate. Cj1437 was shown to catalyze the transamination of dihydroxyacetone phosphate to (S)-serinol phosphate in the presence of l-glutamate. The probable routes to the ultimate formation of the glucuronamide substructures in the capsular polysaccharides of C. jejuni are discussed.

Building better polymerases: Engineering the replication of expanded genetic alphabets.

DNA polymerases are today used throughout scientific research, biotechnology, and medicine, in part for their ability to interact with unnatural forms of DNA created by synthetic biologists. Here especially, natural DNA polymerases often do not have the "performance specifications" needed for transformative technologies. This creates a need for science-guided rational (or semi-rational) engineering to identify variants that replicate unnatural base pairs (UBPs), unnatural backbones, tags, or other evolutionarily novel features of unnatural DNA. In this review, we provide a brief overview of the chemistry and properties of replicative DNA polymerases and their evolved variants, focusing on the Klenow fragment of Taq DNA polymerase (Klentaq). We describe comparative structural, enzymatic, and molecular dynamics studies of WT and Klentaq variants, complexed with natural or noncanonical substrates. Combining these methods provides insight into how specific amino acid substitutions distant from the active site in a Klentaq DNA polymerase variant (ZP Klentaq) contribute to its ability to replicate UBPs with improved efficiency compared with Klentaq. This approach can therefore serve to guide any future rational engineering of replicative DNA polymerases.

Improving the Treatment of Acute Lymphoblastic Leukemia.

l-Asparaginase (EC 3.5.1.1) was first used as a component of combination drug therapies to treat acute lymphoblastic leukemia (ALL), a cancer of the blood and bone marrow, almost 50 years ago. Administering this enzyme to reduce asparagine levels in the blood is a cornerstone of modern clinical protocols for ALL; indeed, this remains the only successful example of a therapy targeted against a specific metabolic weakness in any form of cancer. Three problems, however, constrain the clinical use of l-asparaginase. First, a type II bacterial variant of l-asparaginase is administered to patients, the majority of whom are children, which produces an immune response thereby limiting the time over which the enzyme can be tolerated. Second, l-asparaginase is subject to proteolytic degradation in the blood. Third, toxic side effects are observed, which may be correlated with the l-glutaminase activity of the enzyme. This Perspective will outline how asparagine depletion negatively impacts the growth of leukemic blasts, discuss the structure and mechanism of l-asparaginase, and briefly describe the clinical use of chemically modified forms of clinically useful l-asparaginases, such as Asparlas, which was recently given FDA approval for use in children (babies to young adults) as part of multidrug treatments for ALL. Finally, we review ongoing efforts to engineer l-asparaginase variants with improved therapeutic properties and briefly detail emerging, alternate strategies for the treatment of forms of ALL that are resistant to asparagine depletion.

Consecutive non-natural PZ nucleobase pairs in DNA impact helical structure as seen in 50 µs molecular dynamics simulations.

Z: Little is known about the influence of multiple consecutive 'non-standard' ( , 6-amino-5-nitro-2(1H)-pyridone, and , 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one) nucleobase pairs on the structural parameters of duplex DNA. nucleobase pairs follow standard rules for Watson-Crick base pairing but have rearranged hydrogen bonding donor and acceptor groups. Using the X-ray crystal structure as a starting point, we have modeled the motions of a DNA duplex built from a self-complementary oligonucleotide (5΄-CTTATPPPZZZATAAG-3΄) in water over a period of 50 µs and calculated DNA local parameters, step parameters, helix parameters, and major/minor groove widths to examine how the presence of multiple, consecutive nucleobase pairs might impact helical structure. In these simulations, the -containing DNA duplex exhibits a significantly wider major groove and greater average values of stagger, slide, rise, twist and h-rise than observed for a 'control' oligonucleotide in which nucleobase pairs are replaced by . The molecular origins of these structural changes are likely associated with at least two differences between and . First, the electrostatic properties of differ from in terms of density distribution and dipole moment. Second, differences are seen in the base stacking of pairs in dinucleotide steps, arising from energetically favorable stacking of the nitro group in with π-electrons of the adjacent base.

Second-Shell Hydrogen Bond Impacts Transition-State Structure in Bacillus subtilis Oxalate Decarboxylase.

There is considerable interest in how "second-shell" interactions between protein side chains and metal ligands might modulate Mn(II) ion redox properties and reactivity in metalloenzymes. One such Mn-dependent enzyme is oxalate decarboxylase (OxDC), which catalyzes the disproportionation of oxalate monoanion into formate and CO2. Electron paramagnetic resonance (EPR) studies have shown that a mononuclear Mn(III) ion is formed in OxDC during catalytic turnover and that the removal of a hydrogen bond between one of the metal ligands (Glu101) and a conserved, second-shell tryptophan residue (Trp132) gives rise to altered zero-field splitting parameters for the catalytically important Mn(II) ion. We now report heavy-atom kinetic isotope effect measurements on the W132F OxDC variant, which test the hypothesis that the Glu101/Trp132 hydrogen bond modulates the stability of the Mn(III) ion during catalytic turnover. Our results suggest that removing the Glu101/Trp132 hydrogen bond increases the energy of the oxalate radical intermediate from which decarboxylation takes place. This finding is consistent with a model in which the Glu101/Trp132 hydrogen bond in WT OxDC modulates the redox properties of the Mn(II) ion.

Toward an Expanded Genome: Structural and Computational Characterization of an Artificially Expanded Genetic Information System.

Although the fundamental properties of DNA as first proposed by Watson and Crick in 1953 provided a basic understanding of how duplex DNA was organized and might be replicated, it was not until the first crystal structures of DNA (Z-DNA in 1979, B-DNA in 1980, and A-DNA in 1982) that the true complexity of the molecule began to be appreciated. Many crystal structures of oligonucleotides have since shed light on the helical forms that "Watson-Crick" DNA can adopt, their associated groove widths, and the properties of the nucleobase pairs and their interactions in all three helical forms. Additional understanding of the properties of Watson-Crick DNA has been provided by computational studies employing a variety of theoretical methods. Together with these studies devoted to understanding Watson-Crick DNA, recent efforts to expand the genetic alphabet have founded a new field in synthetic biology. One of these efforts, the artificially expanded genetic information system (AEGIS) developed by Steven Benner and co-workers, takes advantage of orthogonal hydrogen bonding to produce DNA comprised of six nucleobase pairs, of which the most extensively studied is referred to as P:Z with P being 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one) and Z being 6-amino-5-nitro-2(1H)-pyridone. P:Z forms three edge-on hydrogen bonds that differ from standard Watson-Crick pairs in the arrangement of acceptors and donor groups; P presents acceptor, acceptor, donor, and Z presents donor, donor, acceptor. Z is unique among the AEGIS nucleobases in having a nitro group present in the major groove. PZ-containing DNA has been exploited in a number of clinical applications and is being used to develop receptors and catalysts. Ultimately, the grand challenge will be to create a semisynthetic organism with an expanded genome. Furthermore, just as our understanding of the properties of natural DNA have benefited from structural and computational characterization, so too will our understanding of artificial DNA. This Account focuses on the structural and biophysical properties of AEGIS DNA containing P:Z pairs. We begin with the fundamental properties of P:Z nucleobase pairs, including their electrostatic potential and hydrogen-bonding energies, as elucidated by quantum mechanical calculations. We then examine the impact of including multiple consecutive P:Z pairs into duplex DNA providing an opportunity to investigate stacking interactions between P:Z pairs. The self-complementary 5'-CTTATPPTAZZATAAG was crystallized in B-form using the host-guest system along with analogous natural sequences including Gs or As. Use of the host-guest system to characterize B-DNA obviates a number of limitations on the structural characterization of sequences of interest; these include the ability to crystallize the desired sequences and to distinguish structural effects imparted by the lattice constraints from those inherent in the sequence itself. On the other hand, 3/6ZP, 5'-CTTATPPPZZZATAAG, was crystallized in A-form in a DNA-only lattice allowing a comparative analysis of P:Z pairs in two of the biologically relevant helical forms: A- and B-DNA. Computational studies on the 3/6ZP sequence starting in A-form provide additional evidence for a more energetically favorable stacking interaction, which we term the "slide" conformer, observed in the A-form crystal structure; this unusual stacking interaction plays a major role in altering the conformational dynamics observed for the PZ-containing duplex as compared to a GC-containing "control" duplex in long time scale molecular dynamics simulations. This combined use of structural and computational strategies paves the way for obtaining a detailed description of artificial DNA, both in how it differs from Watson-Crick DNA and in the rational discovery of proteins, such as endonucleases, transcription factors, and polymerases, which can specifically manipulate DNA containing AEGIS nucleobase pairs.

Chemoenzymatic Synthesis, Nanotization, and Anti-Aspergillus Activity of Optically Enriched Fluconazole Analogues.

Despite recent advances in diagnostic and therapeutic methods in antifungal research, aspergillosis still remains a leading cause of morbidity and mortality. One strategy to address this problem is to enhance the activity spectrum of known antifungals, and we now report the first successful application of Candida antarctica lipase (CAL) for the preparation of optically enriched fluconazole analogues. Anti-Aspergillus activity was observed for an optically enriched derivative, (-)-S-2-(2',4'-difluorophenyl)-1-hexyl-amino-3-(1‴,2‴,4‴)triazol-1‴-yl-propan-2-ol, which exhibits MIC values of 15.6 µg/ml and 7.8 µg/disc in broth microdilution and disc diffusion assays, respectively. This compound is tolerated by mammalian erythrocytes and cell lines (A549 and U87) at concentrations of up to 1,000 µg/ml. When incorporated into dextran nanoparticles, the novel, optically enriched fluconazole analogue exhibited improved antifungal activity against Aspergillus fumigatus (MIC, 1.63 µg/ml). These results not only demonstrate the ability of biocatalytic approaches to yield novel, optically enriched fluconazole derivatives but also suggest that enantiomerically pure fluconazole derivatives, and their nanotized counterparts, exhibiting anti-Aspergillus activity may have reduced toxicity.

Synthesis of macromolecular systems via lipase catalyzed biocatalytic reactions: applications and future perspectives.

Enzymes, being remarkable catalysts, are capable of accepting a wide range of complex molecules as substrates and catalyze a variety of reactions with a high degree of chemo-, stereo- and regioselectivity in most of the reactions. Biocatalysis can be used in both simple and complex chemical transformations without the need for tedious protection and deprotection chemistry that is very common in traditional organic synthesis. This current review highlights the applicability of one class of biocatalysts viz."lipases" in synthetic transformations, the resolution of pharmaceutically important small molecules including polyphenols, amides, nucleosides and their precursors, the development of macromolecular systems (and their applications as drug/gene carriers), flame retardants, polymeric antioxidants and nanocrystalline solar cells, etc.

Assessing the Influence of Mutation on GTPase Transition States by Using X-ray Crystallography, 19 F†NMR, and DFT Approaches.

We report X-ray crystallographic and 19 F NMR studies of the G-protein RhoA complexed with MgF3- , GDP, and RhoGAP, which has the mutation Arg85'Ala. When combined with DFT calculations, these data permit the identification of changes in transition state (TS) properties. The X-ray data show how Tyr34 maintains solvent exclusion and the core H-bond network in the active site by relocating to replace the missing Arg85' sidechain. The 19 F NMR data show deshielding effects that indicate the main function of Arg85' is electronic polarization of the transferring phosphoryl group, primarily mediated by H-bonding to O3G and thence to PG . DFT calculations identify electron-density redistribution and pinpoint why the TS for guanosine 5'-triphosphate (GTP) hydrolysis is higher in energy when RhoA is complexed with RhoGAPArg85'Ala relative to wild-type (WT) RhoGAP. This study demonstrates that 19 F NMR measurements, in combination with X-ray crystallography and DFT calculations, can reliably dissect the response of small GTPases to site-specific modifications.

Formation of Hexacoordinate Mn(III) in Bacillus subtilis Oxalate Decarboxylase Requires Catalytic Turnover.

Oxalate decarboxylase (OxDC) catalyzes the disproportionation of oxalic acid monoanion into CO2 and formate. The enzyme has long been hypothesized to utilize dioxygen to form mononuclear Mn(III) or Mn(IV) in the catalytic site during turnover. Recombinant OxDC, however, contains only tightly bound Mn(II), and direct spectroscopic detection of the metal in higher oxidation states under optimal catalytic conditions (pH 4.2) has not yet been reported. Using parallel mode electron paramagnetic resonance spectroscopy, we now show that substantial amounts of Mn(III) are indeed formed in OxDC, but only in the presence of oxalate and dioxygen under acidic conditions. These observations provide the first direct support for proposals in which Mn(III) removes an electron from the substrate to yield a radical intermediate in which the barrier to C-C bond cleavage is significantly decreased. Thus, OxDC joins a small list of enzymes capable of stabilizing and controlling the reactivity of the powerful oxidizing species Mn(III).

Substrate Binding Mode and Molecular Basis of a Specificity Switch in Oxalate Decarboxylase.

Oxalate decarboxylase (OxDC) catalyzes the conversion of oxalate into formate and carbon dioxide in a remarkable reaction that requires manganese and dioxygen. Previous studies have shown that replacing an active-site loop segment Ser(161)-Glu(162)-Asn(163)-Ser(164) in the N-terminal domain of OxDC with the cognate residues Asp(161)-Ala(162)-Ser-(163)-Asn(164) of an evolutionarily related, Mn-dependent oxalate oxidase gives a chimeric variant (DASN) that exhibits significantly increased oxidase activity. The mechanistic basis for this change in activity has now been investigated using membrane inlet mass spectrometry (MIMS) and isotope effect (IE) measurements. Quantitative analysis of the reaction stoichiometry as a function of oxalate concentration, as determined by MIMS, suggests that the increased oxidase activity of the DASN OxDC variant is associated with only a small fraction of the enzyme molecules in solution. In addition, IE measurements show that C-C bond cleavage in the DASN OxDC variant proceeds via the same mechanism as in the wild-type enzyme, even though the Glu(162) side chain is absent. Thus, replacement of the loop residues does not modulate the chemistry of the enzyme-bound Mn(II) ion. Taken together, these results raise the possibility that the observed oxidase activity of the DASN OxDC variant arises from an increased level of access of the solvent to the active site during catalysis, implying that the functional role of Glu(162) is to control loop conformation. A 2.6 Å resolution X-ray crystal structure of a complex between oxalate and the Co(II)-substituted ΔE162 OxDC variant, in which Glu(162) has been deleted from the active site loop, reveals the likely mode by which the substrate coordinates the catalytically active Mn ion prior to C-C bond cleavage. The "end-on" conformation of oxalate observed in the structure is consistent with the previously published V/K IE data and provides an empty coordination site for the dioxygen ligand that is thought to mediate the formation of Mn(III) for catalysis upon substrate binding.

Substrate Distortion and the Catalytic Reaction Mechanism of 5-Carboxyvanillate Decarboxylase.

5-Carboxyvanillate decarboxylase (LigW) catalyzes the conversion of 5-carboxyvanillate to vanillate in the biochemical pathway for the degradation of lignin. This enzyme was shown to require Mn(2+) for catalytic activity and the kinetic constants for the decarboxylation of 5-carboxyvanillate by the enzymes from Sphingomonas paucimobilis SYK-6 (kcat = 2.2 s(-1) and kcat/Km = 4.0 × 10(4) M(-1) s(-1)) and Novosphingobium aromaticivorans (kcat = 27 s(-1) and kcat/Km = 1.1 × 10(5) M(-1) s(-1)) were determined. The three-dimensional structures of both enzymes were determined in the presence and absence of ligands bound in the active site. The structure of LigW from N. aromaticivorans, bound with the substrate analogue, 5-nitrovanillate (Kd = 5.0 nM), was determined to a resolution of 1.07 Å. The structure of this complex shows a remarkable enzyme-induced distortion of the nitro-substituent out of the plane of the phenyl ring by approximately 23°. A chemical reaction mechanism for the decarboxylation of 5-carboxyvanillate by LigW was proposed on the basis of the high resolution X-ray structures determined in the presence ligands bound in the active site, mutation of active site residues, and the magnitude of the product isotope effect determined in a mixture of H2O and D2O. In the proposed reaction mechanism the enzyme facilitates the transfer of a proton to C5 of the substrate prior to the decarboxylation step.

Gas phase RDX decomposition pathways using coupled cluster theory.

Electronic and free energy barriers for a series of gas-phase RDX decomposition mechanisms have been obtain using coupled cluster singles, doubles, and perturbative triples with complete basis set (CCSD(T)/CBS) electronic energies for MBPT(2)/cc-pVTZ structures. Importantly, we have located a well-defined transition state for NN homolysis, in the initial RDX decomposition step, thereby obtaining a true barrier for this reaction. These calculations support the view that HONO elimination is preferred at STP over other proposed mechanisms, including NN homolysis, "triple whammy" and NONO isomerization. Indeed, our calculated values of Arrhenius parameters are in agreement with experimental findings for gas phase RDX decomposition. We also investigate a number of new pathways leading to breakdown of the intermediate formed by the initial HONO elimination, and find that NN homolysis in this intermediate has an activation energy barrier comparable with that computed for HONO elimination.

Synthesis and anti-inflammatory activity evaluation of novel triazolyl-isatin hybrids.

New isatin-triazole based hybrids have been synthesized and evaluated for their inhibitory activity of TNF-α induced expression of Intercellular Adhesion Molecule-1 (ICAM-1) on the surface of human endothelial cells. Structure-activity relationship (SAR) studies revealed that the presence of the electron-attracting bromo substituent at position-5 of the isatin moiety played an important role in enhancing the anti-inflammatory potential of the synthesized compounds. Z-1-[3-(1H-1,2,4-Triazol-1-yl)propyl]-5-bromo-3-[2-(4-methoxyphenyl)hydrazono]indolin-2-one (19) with an IC50 = 20 µM and 89% ICAM-1 inhibition with MTD at 200 µM was found to be the most potent of all the synthesized derivatives. Introduction of 1,2,4-triazole ring and electron-donating methoxy group on the phenylhydrazone moiety resulted in four-fold increase of the anti-inflammatory activity.

(19)F†NMR and DFT Analysis Reveal Structural and Electronic Transition State Features for RhoA-Catalyzed GTP Hydrolysis.

Molecular details for RhoA/GAP catalysis of the hydrolysis of GTP to GDP are poorly understood. We use (19)F NMR chemical shifts in the MgF3(-) transition state analogue (TSA) complex as a spectroscopic reporter to indicate electron distribution for the γ-PO3(-) oxygens in the corresponding TS, implying that oxygen coordinated to Mg has the greatest electron density. This was validated by QM calculations giving a picture of the electronic properties of the transition state (TS) for nucleophilic attack of water on the γ-PO3(-) group based on the structure of a RhoA/GAP-GDP-MgF3(-) TSA complex. The TS model displays a network of 20 hydrogen bonds, including the GAP Arg85' side chain, but neither phosphate torsional strain nor general base catalysis is evident. The nucleophilic water occupies a reactive location different from that in multiple ground state complexes, arising from reorientation of the Gln-63 carboxamide by Arg85' to preclude direct hydrogen bonding from water to the target γ-PO3(-) group.

Structural basis for a six nucleotide genetic alphabet.

Expanded genetic systems are most likely to work with natural enzymes if the added nucleotides pair with geometries that are similar to those displayed by standard duplex DNA. Here, we present crystal structures of 16-mer duplexes showing this to be the case with two nonstandard nucleobases (Z, 6-amino-5-nitro-2(1H)-pyridone and P, 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)one) that were designed to form a Z:P pair with a standard "edge on" Watson-Crick geometry, but joined by rearranged hydrogen bond donor and acceptor groups. One duplex, with four Z:P pairs, was crystallized with a reverse transcriptase host and adopts primarily a B-form. Another contained six consecutive Z:P pairs; it crystallized without a host in an A-form. In both structures, Z:P pairs fit canonical nucleobase hydrogen-bonding parameters and known DNA helical forms. Unique features include stacking of the nitro group on Z with the adjacent nucleobase ring in the A-form duplex. In both B- and A-duplexes, major groove widths for the Z:P pairs are approximately 1 Å wider than those of comparable G:C pairs, perhaps to accommodate the large nitro group on Z. Otherwise, ZP-rich DNA had many of the same properties as CG-rich DNA, a conclusion supported by circular dichroism studies in solution. The ability of standard duplexes to accommodate multiple and consecutive Z:P pairs is consistent with the ability of natural polymerases to biosynthesize those pairs. This, in turn, implies that the GACTZP synthetic genetic system can explore the entire expanded sequence space that additional nucleotides create, a major step forward in this area of synthetic biology.

Facile C(sp(2))-C(sp(2)) bond cleavage in oxalic acid-derived radicals.

Oxalate decarboxylase (OxDC) catalyzes the Mn-dependent conversion of the oxalate monoanion into CO2 and formate. Many questions remain about the catalytic mechanism of OxDC although it has been proposed that the reaction proceeds via substrate-based radical intermediates. Using coupled cluster theory combined with implicit solvation models we have examined the effects of radical formation on the structure and reactivity of oxalic acid-derived radicals in aqueous solution. Our results show that the calculated solution-phase free-energy barrier for C-C bond cleavage to form CO2 is decreased from 34.2 kcal/mol for oxalic acid to only 9.3 kcal/mol and a maximum of 3.5 kcal/mol for the cationic and neutral oxalic acid-derived radicals, respectively. These studies also show that the C-C σ bonding orbital of the radical cation contains only a single electron, giving rise to an elongated C-C bond distance of 1.7 Å; a similar lengthening of the C-C bond is not observed for the neutral radical. This study provides new chemical insights into the structure and stability of plausible intermediates in the catalytic mechanism of OxDC, and suggests that removal of an electron to form a radical (with or without the concomitant loss of a proton) may be a general strategy for cleaving the unreactive C-C bonds between adjacent sp(2)-hybridized carbon atoms.

A mechanochemical switch to control radical intermediates.

B12-dependent enzymes employ radical species with exceptional prowess to catalyze some of the most chemically challenging, thermodynamically unfavorable reactions. However, dealing with highly reactive intermediates is an extremely demanding task, requiring sophisticated control strategies to prevent unwanted side reactions. Using hybrid quantum mechanical/molecular mechanical simulations, we follow the full catalytic cycle of an AdoB12-dependent enzyme and present the details of a mechanism that utilizes a highly effective mechanochemical switch. When the switch is "off", the 5'-deoxyadenosyl radical moiety is stabilized by releasing the internal strain of an enzyme-imposed conformation. Turning the switch "on," the enzyme environment becomes the driving force to impose a distinct conformation of the 5'-deoxyadenosyl radical to avoid deleterious radical transfer. This mechanochemical switch illustrates the elaborate way in which enzymes attain selectivity of extremely chemically challenging reactions.

Assigning the EPR fine structure parameters of the Mn(II) centers in Bacillus subtilis oxalate decarboxylase by site-directed mutagenesis and DFT/MM calculations.

Oxalate decarboxylase (OxDC) catalyzes the Mn-dependent conversion of the oxalate monoanion into CO2 and formate. EPR-based strategies for investigating the catalytic mechanism of decarboxylation are complicated by the difficulty of assigning the signals associated with the two Mn(II) centers located in the N- and C-terminal cupin domains of the enzyme. We now report a mutational strategy that has established the assignment of EPR fine structure parameters to each of these Mn(II) centers at pH 8.5. These experimental findings are also used to assess the performance of a multistep strategy for calculating the zero-field splitting parameters of protein-bound Mn(II) ions. Despite the known sensitivity of calculated D and E values to the computational approach, we demonstrate that good estimates of these parameters can be obtained using cluster models taken from carefully optimized DFT/MM structures. Overall, our results provide new insights into the strengths and limitations of theoretical methods for understanding electronic properties of protein-bound Mn(II) ions, thereby setting the stage for future EPR studies on the electronic properties of the Mn(II) centers in OxDC and site-specific variants.