Biomedical subjects

David A Case

Publications and source records attributed to David A Case.

At least 19 recordsLinked to original sources

PubMed2005

Biomolecular simulations at constant pH.

Like temperature and pressure, the solution pH is an important thermodynamic variable that is commonly varied in experiments and is used by cells to influence biochemical function. It is now becoming feasible to carry out practical molecular dynamics simulations that mimic the thermodynamics of such experiments, by allowing proton transfer between the system of interest and a hypothetical bath of protons at a given pH. These are demanding calculations, because the energetics of charge changes upon protonation or deprotonation must be accurately modeled, and because such simulations must sample both molecular configurations and the large number of protonation states that are possible for a molecule with many titrating sites.

Algorithms↗

PubMed2004

Predicting interactions of winged-helix transcription factors with DNA.

Determining protein-DNA interactions is important for understanding gene regulation, DNA repair and chromatin structure. Unfortunately, the structures of DNA-bound complexes are often difficult to obtain experimentally, so the development of computational methods that provide good models of these complexes would be valuable. Here, we present a rigid-body docking approach using the computer program DOT. DOT performs a complete, six-dimensional search of all orientations for two rigid molecules and calculates the interaction energy as the sum of electrostatic and van der Waals terms. DOT was applied to three winged-helix transcription factors that share similar DNA-binding structural motifs but bind DNA in different ways. Docking with linear B-form DNA models accomplished several objectives; it (1) distinguished the different ways the transcription factors bind DNA, (2) identified each protein's DNA-binding site and the DNA orientation at the site and (3) gave at least one solution among the three best-ranked that shows the protein side chain-DNA base interactions responsible for recognition. Furthermore, the ensemble of top-ranked, docked linear B-DNA fragments indicated the DNA bending induced upon protein binding. Docking linear B-DNA to structures of the transcription factor FadR suggests that the allosteric, conformational change induced upon effector binding results in loss of the ability to bend DNA as well as loss of sequence-specific interactions with DNA. The electrostatic energy term calculated by DOT is comparable to the electrostatic binding energy calculated by Poisson-Boltzmann methods. Our results show rigid-body docking that includes a rigorous treatment of the electrostatic interaction energy can be effective in predicting protein-DNA interactions.

Bacterial Proteins↗

PubMed2004

Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. I. Research design and results on d(CpG) steps.

We describe herein a computationally intensive project aimed at carrying out molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide base sequences. This initiative was undertaken by an international collaborative effort involving nine research groups, the "Ascona B-DNA Consortium" (ABC). Calculations were carried out on the 136 cases imbedded in 39 DNA oligomers with repeating tetranucleotide sequences, capped on both ends by GC pairs and each having a total length of 15 nucleotide pairs. All MD simulations were carried out using a well-defined protocol, the AMBER suite of programs, and the parm94 force field. Phase I of the ABC project involves a total of approximately 0.6 mus of simulation for systems containing approximately 24,000 atoms. The resulting trajectories involve 600,000 coordinate sets and represent approximately 400 gigabytes of data. In this article, the research design, details of the simulation protocol, informatics issues, and the organization of the results into a web-accessible database are described. Preliminary results from 15-ns MD trajectories are presented for the d(CpG) step in its 10 unique sequence contexts, and issues of stability and convergence, the extent of quasiergodic problems, and the possibility of long-lived conformational substates are discussed.

Base Sequence↗

PubMed2004

Change in protein flexibility upon complex formation: analysis of Ras-Raf using molecular dynamics and a molecular framework approach.

Changes in flexibility upon protein-protein complex formation of H-Ras and the Ras-binding domain of C-Raf1 have been investigated using the molecular framework approach FIRST (Floppy Inclusion and Rigid Substructure Topology) and molecular dynamics simulations (MD) of in total approximately 35 ns length. In a computational time of about one second, FIRST identifies flexible and rigid regions in a single, static three-dimensional molecular framework, whose vertices represent protein atoms and whose edges represent covalent and non-covalent (hydrogen bond and hydrophobic) constraints and fixed bond angles within the protein. The two methods show a very good agreement with respect to the identification of changes in flexibility in both binding partners on a local scale. This implies that flexibility can be successfully predicted by identifying which bonds limit motion within a molecule and how they are coupled. In particular, as identified by MD, the beta-sheet in Raf shows considerably more pronounced orientational correlations in the bound state compared to the unbound state. Similarly, FIRST assigns the beta-sheet to the largest rigid cluster of the complex. Interestingly, FIRST allows us to identify that interactions across the interface (but not conformational changes upon complex formation) result in the observed rigidification. Since regions of the beta-sheet of Raf that do not interact directly with Ras become rigidified, this also demonstrates the long-range aspect to rigidity percolation. Possible implications of the change of flexibility of the Ras-binding domain of Raf on the activation of Raf upon complex formation are discussed. Finally, the sensitivity of FIRST results with respect to the representation of non-covalent interactions used as constraints is probed.

Algorithms↗

PubMed2004

Development and testing of a general amber force field.

We describe here a general Amber force field (GAFF) for organic molecules. GAFF is designed to be compatible with existing Amber force fields for proteins and nucleic acids, and has parameters for most organic and pharmaceutical molecules that are composed of H, C, N, O, S, P, and halogens. It uses a simple functional form and a limited number of atom types, but incorporates both empirical and heuristic models to estimate force constants and partial atomic charges. The performance of GAFF in test cases is encouraging. In test I, 74 crystallographic structures were compared to GAFF minimized structures, with a root-mean-square displacement of 0.26 A, which is comparable to that of the Tripos 5.2 force field (0.25 A) and better than those of MMFF 94 and CHARMm (0.47 and 0.44 A, respectively). In test II, gas phase minimizations were performed on 22 nucleic acid base pairs, and the minimized structures and intermolecular energies were compared to MP2/6-31G* results. The RMS of displacements and relative energies were 0.25 A and 1.2 kcal/mol, respectively. These data are comparable to results from Parm99/RESP (0.16 A and 1.18 kcal/mol, respectively), which were parameterized to these base pairs. Test III looked at the relative energies of 71 conformational pairs that were used in development of the Parm99 force field. The RMS error in relative energies (compared to experiment) is about 0.5 kcal/mol. GAFF can be applied to wide range of molecules in an automatic fashion, making it suitable for rational drug design and database searching.

Nucleic Acids↗

PubMed2004

NMR structure of a cyclic polyamide-DNA complex.

The solution structure of a cyclic polyamide ligand complexed to a DNA oligomer, derived from NMR restrained molecular mechanics, is presented. The polyamide, cyclo-gamma-ImPyPy-gamma-PyPyPy-, binds to target DNA with a nanomolar dissociation constant as characterized by quantitative footprinting previously reported. 2D (1)H NMR data were used to generate distance restraints defining the structure of this cyclic polyamide with the DNA duplex d(5'-GCCTGTTAGCG-3'):d(5'-CGCTAACAGGC-3'). Data interpretation used complete relaxation matrix analysis of the NOESY cross-peak intensities with the program MARDIGRAS. The NMR-based distance restraints (276 total) were applied in restrained molecular dynamics calculations using a solvent model, yielding structures with an rmsd for the ligand and binding site of approximately 1 A. The resulting structures indicate some distortion of the DNA in the binding site. The constraints from cyclization lead to altered stacking of the rings in the halves of the cyclic ligand relative to unlinked complexes. Despite this, the interactions with DNA are very similar to what has been found in unlinked complexes. Measurements of ligand amide and DNA imino proton exchange rates indicate very slow dissociation of the ligand and show that the DNA can undergo opening fluctuations while the ligand is bound although the presence of the ligand decreases their frequency relative to the free DNA.

Amides↗

PubMed2004

Exploring protein native states and large-scale conformational changes with a modified generalized born model.

Implicit solvation models provide, for many applications, a reasonably accurate and computationally effective way to describe the electrostatics of aqueous solvation. Here, a popular analytical Generalized Born (GB) solvation model is modified to improve its accuracy in calculating the solvent polarization part of free energy changes in large-scale conformational transitions, such as protein folding. In contrast to an earlier GB model (implemented in the AMBER-6 program), the improved version does not overstabilize the native structures relative to the finite-difference Poisson-Boltzmann continuum treatment. In addition to improving the energy balance between folded and unfolded conformers, the algorithm (available in the AMBER-7 and NAB molecular modeling packages) is shown to perform well in more than 50 ns of native-state molecular dynamics (MD) simulations of thioredoxin, protein-A, and ubiquitin, as well as in a simulation of Barnase/Barstar complex formation. For thioredoxin, various combinations of input parameters have been explored, such as the underlying gas-phase force fields and the atomic radii. The best performance is achieved with a previously proposed modification to the torsional potential in the Amber ff99 force field, which yields stable native trajectories for all of the tested proteins, with backbone root-mean-square deviations from the native structures being approximately 1.5 A after 6 ns of simulation time. The structure of Barnase/Barstar complex is regenerated, starting from an unbound state, to within 1.9 A relative to the crystal structure of the complex.

Bacterial Proteins↗

PubMed2004

Proton binding to proteins: pK(a) calculations with explicit and implicit solvent models.

Ionizable residues play important roles in protein structure and activity, and proton binding is a valuable reporter of electrostatic interactions in these systems. We use molecular dynamics free energy simulations (MDFE) to compute proton pKa shifts, relative to a model compound in solution, for three aspartate side chains in two proteins. Simulations with explicit solvent and with an implicit, dielectric continuum solvent are reported. The implicit solvent simulations use the generalized Born (GB) model, which provides an approximate, analytical solution to Poisson's equation. With explicit solvent, the direction of the pKa shifts is correct in all three cases with one force field (AMBER) and in two out of three cases with another (CHARMM). For two aspartates, the dielectric response to ionization is found to be linear, even though the separate protein and solvent responses can be nonlinear. For thioredoxin Asp26, nonlinearity arises from the presence of two substates that correspond to the two possible orientations of the protonated carboxylate. For this side chain, which is partly buried and has a large pKa upshift, very long simulations are needed to correctly sample several slow degrees of freedom that reorganize in response to the ionization. Thus, nearby Lys57 rotates to form a salt bridge and becomes buried, while three waters intercalate along the opposite edge of Asp26. Such strong and anisotropic reorganization is very difficult to predict with Poisson-Boltzmann methods that only consider electrostatic interactions and employ a single protein structure. In contrast, MDFE with a GB dielectric continuum solvent, used for the first time for pKa calculations, can describe protein reorganization accurately and gives encouraging agreement with experiment and with the explicit solvent simulations.

Animals↗

PubMed2004

High-resolution field-cycling NMR studies of a DNA octamer as a probe of phosphodiester dynamics and comparison with computer simulation.

Phosphorus-spin longitudinal relaxation rates of the DNA duplex octamer [d(GGAATTCC)](2) have been measured from 0.1 to 17.6 T by means of conventional and new field-cycling NMR methods. The high-resolution field-cycling method is identical to a conventional relaxation experiment, except that after preparation the sample is moved pneumatically from its usual position at the center of the high-resolution magnet upward to a lower field above its normal position and then returned to the center for readout after it has relaxed for the programmed relaxation delay at the low field. This is the first measurement of all longitudinal relaxation rates R(1) of a nuclear species in a macromolecule over virtually the entire accessible magnetic field range. For detailed analysis, three magnetic field regions can be delineated: (i) dipolar relaxation dominates at fields below 2 T, (ii) chemical shift anisotropy (CSA) relaxation is roughly constant from 2 to 6 T, and (iii) a square-law increasing dependence is seen at fields higher than approximately 6 T due to internal motion CSA relaxation. The analysis provides a rotational correlation time (tau(r) = 4.1 +/- 0.3 ns) for the duplex at both 1.5 and 0.25 mM concentrations (of duplex) at 22 degrees C. For comparison, extraction of tau(r) in the conventional way from the ratio of T(1)/T(2) at 14 T yields 3.2 ns. The tau(r) discrepancy disappears when we exclude the contribution of internal motion from the R(1) in the ratio. The low-field dipolar relaxation provides a weighted inverse sixth power sum of the distances from the phosphorus to the protons responsible for relaxation. This average is similar for all phosphates in the octamer and similar to that in previous B-DNA structures (its inverse sixth root is about 2.40 A for two different concentrations of octamer). The CSA relaxation at intermediate field provides an estimate of the order parameter squared, S(c)(2), for each phosphorus. S(c)(2) is about 0.7-1, clearly different for different phosphate linkages in the octamer duplex. The increasing R(1) at high fields reflects CSA relaxation due to internal motions, for which a correlation time, tau(hf), can be approximately extracted with the aid of additional measurements at 14.0 and 17.6 T. We conclude that tau(hf) values are relatively large, in the range of about 150 ps. Insight into the motions leading to this correlation time was gained by a 28 ns molecular dynamics simulation of the molecule. S(2) and tau(s) (corresponding to tau(hf)) predicted by this simulation were in good agreement with the experimental values from the field-cycling data. Both the effect of Mg(2+) on the dynamic parameters extracted from (31)P relaxation rates and the field dependence of relaxation rates for several protons of the octamer were measured. High-resolution field cycling opens up the possibility of monitoring residue-specific dipolar interactions and dynamics for the phosphorus nuclei of diverse oligonucleotides.

Cations, Divalent↗

PubMed2004

Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures.

This study compares generalized Born (GB) and Poisson (PB) methods for calculating electrostatic solvation energies of proteins. A large set of GB and PB implementations from our own laboratories as well as others is applied to a series of protein structure test sets for evaluating the performance of these methods. The test sets cover a significant range of native protein structures of varying size, fold topology, and amino acid composition as well as nonnative extended and misfolded structures that may be found during structure prediction and folding/unfolding studies. We find that the methods tested here span a wide range from highly accurate and computationally demanding PB-based methods to somewhat less accurate but more affordable GB-based approaches and a few fast, approximate PB solvers. Compared with PB solvation energies, the latest, most accurate GB implementations were found to achieve errors of 1% for relative solvation energies between different proteins and 0.4% between different conformations of the same protein. This compares to accurate PB solvers that produce results with deviations of less than 0.25% between each other for both native and nonnative structures. The performance of the best GB methods is discussed in more detail for the application for force field-based minimizations or molecular dynamics simulations.

Hydrogen Bonding↗

PubMed2004

Converging free energy estimates: MM-PB(GB)SA studies on the protein-protein complex Ras-Raf.

Estimating protein-protein interaction energies is a very challenging task for current simulation protocols. Here, absolute binding free energies are reported for the complex H-Ras/C-Raf1 using the MM-PB(GB)SA approach, testing the internal consistency and model dependence of the results. Averaging gas-phase energies (MM), solvation free energies as determined by Generalized Born models (GB/SA), and entropic contributions calculated by normal mode analysis for snapshots obtained from 10 ns explicit-solvent molecular dynamics in general results in an overestimation of the binding affinity when a solvent-accessible surface area-dependent model is used to estimate the nonpolar solvation contribution. Applying the sum of a cavity solvation free energy and explicitly modeled solute-solvent van der Waals interaction energies instead provides less negative estimates for the nonpolar solvation contribution. When the polar contribution to the solvation free energy is determined by solving the Poisson-Boltzmann equation (PB) instead, the calculated binding affinity strongly depends on the atomic radii set chosen. For three GB models investigated, different absolute deviations from PB energies were found for the unbound proteins and the complex. As an alternative to normal-mode calculations, quasiharmonic analyses have been performed to estimate entropic contributions due to changes of solute flexibility upon binding. However, such entropy estimates do not converge after 10 ns of simulation time, indicating that sampling issues may limit the applicability of this approach. Finally, binding free energies estimated from snapshots of the unbound proteins extracted from the complex trajectory result in an underestimate of binding affinity. This points to the need to exercise caution in applying the computationally cheaper "one-trajectory-alternative" to systems where there may be significant changes in flexibility and structure due to binding. The best estimate for the binding free energy of Ras-Raf obtained in this study of -8.3 kcal mol(-1) is in good agreement with the experimental result of -9.6 kcal mol(-1), however, further probing the transferability of the applied protocol that led to this result is necessary.

Calorimetry↗

PubMed2004

Constant pH molecular dynamics in generalized Born implicit solvent.

A new method is proposed for constant pH molecular dynamics (MD), employing generalized Born (GB) electrostatics. Protonation states are modeled with different charge sets, and titrating residues sample a Boltzmann distribution of protonation states as the simulation progresses, using Monte Carlo sampling based on GB-derived energies. The method is applied to four different crystal structures of hen egg-white lysozyme (HEWL). pK(a) predictions derived from the simulations have root-mean-square (RMS) error of 0.82 relative to experimental values. Similarity of results between the four crystal structures shows the method to be independent of starting crystal structure; this is in contrast to most electrostatics-only models. A strong correlation between conformation and protonation state is noted and quantitatively analyzed, emphasizing the importance of sampling protonation states in conjunction with dynamics.

Algorithms↗

PubMed2003

Reintroducing electrostatics into protein X-ray structure refinement: bulk solvent treated as a dielectric continuum.

Structural refinement of proteins involves the minimization of a target function that combines X-ray data with a set of restraints enforcing stereochemistry and packing. Electrostatic interactions are not ordinarily included in the target function, partly because they cannot be calculated reliably without a description of dielectric screening by solvent in the crystal. With the recent development of accurate implicit solvent models to describe this screening, the question arises as to whether a more detailed target function including electrostatic and solvation terms can yield more accurate structures or somewhat different structures of equivalent accuracy. The Generalized Born (GB) model is one such model that describes the solvent as a dielectric continuum, taking into account its heterogeneous distribution within the crystal. It is used here for X-ray refinements of three protein structures with experimental diffraction data to 2.4, 2.9 and 3.2 A, respectively. In each case, a higher resolution structure is available for comparison. The new target function includes stereochemical restraints, van der Waals, Coulomb and solvation interactions, along with the usual X-ray pseudo-energy term, which employs the likelihood estimator of Pannu and Read. Multiple simulated-annealing refinements were performed in torsion-angle space with a conventional target function and the new GB target function, yielding ensembles of refined structures. The new target function yields structures of similar accuracy, as measured by the free R factor, map/model correlations and deviations from the high-resolution structures. About 10% of side-chain conformations differ between the two sets of refinements, in the sense that the two ensembles of conformations do not completely overlap. Over 75% of the differences correspond to surface side chains. For one of the proteins, the GB set has a greater dispersion, indicating that for this case the conventional target function overestimates the true precision. As GB parameterization continues to improve, we expect that this approach will become increasingly useful.

Aspartate-tRNA Ligase↗

PubMed2003

Design and synthesis of highly constrained factor Xa inhibitors: amidine-substituted bis(benzoyl)--diazepan-2-ones and bis(benzylidene)-bis(gem-dimethyl)cycloketones.

Two conformationally constrained templates have been designed to provide selective inhibitors of the coagulation cascade serine protease, Factor Xa (FXa). The most active inhibitor, 2,7-bis[(Z)-p-amidinobenzylidene)]-3,3,6,6-tetramethylcycloheptanone, exhibits a K(i) of 42 nM against FXa, with strong selectivity against thrombin (1000-fold), trypsin (300-fold) and plasmin (900-fold). With only two freely rotatable bonds, molecular modeling suggests that one amidine group is positioned into the S1 pocket, forming hydrogen bonds with the side chain of Asp189, similar to other amidine-based inhibitors, with the second benzamidine positioned into the S4 pocket in a position to form strong cation-pi bonding with the S4 aryl cage. We suggest that this interaction plays an important role in the specificity of these inhibitors against other serine proteases.

Antithrombin III↗

PubMed2003

Insights into the mobility of methyl-bearing side chains in proteins from (3)J(CC) and (3)J(CN) couplings.

Side-chain dynamics in proteins can be characterized by the NMR measurement of (13)C and (2)H relaxation rates. Evaluation of the corresponding spectral densities limits the slowest motions that can be studied quantitatively to the time scale on which the overall molecular tumbling takes place. A different measure for the degree of side-chain order about the C(alpha)-C(beta) bond (chi(1) angle) can be derived from (3)J(C)(')(-)(C)(gamma) and (3)J(N)(-)(C)(gamma) couplings. These couplings can be measured at high accuracy, in particular for Thr, Ile, and Val residues. In conjunction with the known backbone structures of ubiquitin and the third IgG-binding domain of protein G, and an extensive set of (13)C-(1)H side-chain dipolar coupling measurements in oriented media, these (3)J couplings were used to parametrize empirical Karplus relationships for (3)J(C)(')(-)(C)(gamma) and (3)J(N)(-)(C)(gamma). These Karplus curves agree well with results from DFT calculations, including an unusual phase shift, which causes the maximum (3)J(CC) and (3)J(CN) couplings to occur for dihedral angles slightly smaller than 180 degrees, particularly noticeable in Thr residues. The new Karplus curves permit determination of rotamer populations for the chi(1) torsion angles. Similar rotamer populations can be derived from side-chain dipolar couplings. Conversion of these rotamer populations into generalized order parameters, S(J)(2) and S(D)(2), provides a view of side-chain dynamics that is complementary to that obtained from (13)C and (2)H relaxation. On average, results agree well with literature values for (2)H-relaxation-derived S(rel)(2) values in ubiquitin and HIV protease, but also identify a fraction of residues for which S(J,D)(2) < S(rel)(2). This indicates that some of the rotameric averaging occurs on a time scale too slow to be observable in traditional relaxation measurements.

HIV Protease↗

PubMed2003

Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes.

Absolute binding free energy calculations and free energy decompositions are presented for the protein-protein complexes H-Ras/C-Raf1 and H-Ras/RalGDS. Ras is a central switch in the regulation of cell proliferation and differentiation. In our study, we investigate the capability of the molecular mechanics (MM)-generalized Born surface area (GBSA) approach to estimate absolute binding free energies for the protein-protein complexes. Averaging gas-phase energies, solvation free energies, and entropic contributions over snapshots extracted from trajectories of the unbound proteins and the complexes, calculated binding free energies (Ras-Raf: -15.0(+/-6.3)kcal mol(-1); Ras-RalGDS: -19.5(+/-5.9)kcal mol(-1)) are in fair agreement with experimentally determined values (-9.6 kcal mol(-1); -8.4 kcal mol(-1)), if appropriate ionic strength is taken into account. Structural determinants of the binding affinity of Ras-Raf and Ras-RalGDS are identified by means of free energy decomposition. For the first time, computationally inexpensive generalized Born (GB) calculations are applied in this context to partition solvation free energies along with gas-phase energies between residues of both binding partners. For selected residues, in addition, entropic contributions are estimated by classical statistical mechanics. Comparison of the decomposition results with experimentally determined binding free energy differences for alanine mutants of interface residues yielded correlations with r(2)=0.55 and 0.46 for Ras-Raf and Ras-RalGDS, respectively. Extension of the decomposition reveals residues as far apart as 25A from the binding epitope that can contribute significantly to binding free energy. These "hotspots" are found to show large atomic fluctuations in the unbound proteins, indicating that they reside in structurally less stable regions. Furthermore, hotspot residues experience a significantly larger-than-average decrease in local fluctuations upon complex formation. Finally, by calculating a pair-wise decomposition of interactions, interaction pathways originating in the binding epitope of Raf are found that protrude through the protein structure towards the loop L1. This explains the finding of a conformational change in this region upon complex formation with Ras, and it may trigger a larger structural change in Raf, which is considered to be necessary for activation of the effector by Ras.

Biophysical Phenomena↗

PubMed2003

Structural, spectroscopic, and redox consequences of a central ligand in the FeMoco of nitrogenase: a density functional theoretical study.

Broken symmetry density functional and electrostatics calculations have been used to shed light on which of three proposed atoms, C, N, or O, is most likely to be present in the center of the FeMoco, the active site of nitrogenase. At the Mo(4+)4Fe(2+)3Fe(3+) oxidation level, a central N(3-) anion results in (1) calculated Fe-N bond distances that are in very good agreement with the recent high-resolution X-ray data of Einsle et al.; (2) a calculated redox potential of 0.19 eV versus the standard hydrogen electrode (SHE) for FeMoco(oxidized) + e(-) --> FeMoco(resting), in good agreement with the measured value of -0.042 V in Azotobacter vinelandii; and (3) average Mössbauer isomer shift values (IS(av) = 0.48 mm s(-1)) compatible with experiment (IS(av) = 0.40 mm s(-1)). At the more reduced Mo(4+)6Fe(2+)1Fe(3+) level, the calculated geometry around a central N(3-) anion still correlates well with the X-ray data, but the average Mössbauer isomer shift value (IS(av) = 0.54 mm s(-1)) and the redox potential of -2.21 eV show a much poorer agreement with experiment. These calculated structural, spectroscopic, and redox data indicate the most likely iron oxidation state for the resting FeMoco of nitrogenase to be 4Fe(2+)3Fe(3+). At this favored oxidation state, oxygen or carbon coordination leads to (1) Fe-O distances in poor agreement and Fe-C distances in good agreement with experiment and (2) calculated redox potentials of +0.97 eV for O(2-) and -1.31 eV for C(4-). The calculated structural parameters and/or redox data suggest either O(2-) or C(4-) is unlikely as a central anion.

Molybdoferredoxin↗

PubMed2003

Symmetry and bonding in metalloporphyrins. A modern implementation for the bonding analyses of five- and six-coordinate high-spin iron(III)-porphyrin complexes through density functional calculation and NMR spectroscopy.

Bonding interactions between the iron and the porphyrin macrocycle of five- and six-coordinate high-spin iron(III)-porphyrin complexes are analyzed within the framework of approximate density functional theory with the use of the quantitative energy decomposition scheme in combination with removal of the vacant pi orbitals of the porphyrin from the valence space. Although the relative extent of the iron-porphyrin interactions can be evaluated qualitatively through the spin population and orbital contribution analyses, the bond strengths corresponding to different symmetry representations can be only approximated quantitatively by the orbital interaction energies. In contrast to previous suggestions, there are only limited Fe --> P pi back-bonding interactions in high-spin iron(III)-porphyrin complexes. It is the symmetry-allowed bonding interaction between d(z)2 and a(2u) orbitals that is responsible for the positive pi spin densities at the meso-carbons of five-coordinate iron(III)-porphyrin complexes. Both five- and six-coordinate complexes show significant P --> Fe pi donation, which is further enhanced by the movement of the metal toward the in-plane position for six-coordinate complexes. These bonding characteristics correlate very well with the NMR data reported experimentally. The extraordinary bonding interaction between d(z)2 and a(2u) orbitals in five-coordinate iron(III)-porphyrin complexes offers a novel symmetry-controlled mechanism for spin transfer between the axial ligand sigma system and the porphyrin pi system and may be critical to the electron transfer pathways mediated by hemoproteins.

Ferric Compounds↗