A stable hyponitrite-bridged iron porphyrin complex.

The coupling of two nitric oxide (NO) molecules in heme active sites is an important contributor to the conversion of NO to nitrous oxide (N(2)O) by heme-containing enzymes. Several formulations for the presumed heme-Fe{N(2)O(2)}(n-) intermediates have been proposed previously, however, no crystal structures of heme-Fe{N(2)O(2)}(n-) systems have been reported to date. We report the first isolation and characterization of a stable bimetallic hyponitrite iron porphyrin, [(OEP)Fe](2)(mu-N(2)O(2)), prepared from the reaction of [(OEP)Fe](2)(mu-O) with hyponitrous acid. Density functional theoretical calculations were performed on the model compound [(porphine)Fe](2)(mu-N(2)O(2)) to characterize its electronic structure and properties.

Revisiting magnesium chelation by teichoic acid with phosphorus solid-state NMR and theoretical calculations.

Teichoic acids are essential components of the Gram-positive bacterial cell wall. One of their many functions is metal binding, a vital process for bacterial growth. With the combination of phosphorus-31 solid-state NMR spectroscopy and theoretical calculations using density functional theory (DFT), we have determined that the binding mode between teichoic acids and magnesium involves bidentate coordination by the phosphate groups of teichoic acid. Measurement of chemical shift anisotropy tensors gave a reduced anisotropy (delta) of 49.25 ppm and an asymmetry (eta) of 0.7. DFT calculations with diglycerol phosphate and triglycerol diphosphate model compounds were completed with Mg(2+) in anhydrous as well as partially hydrated bidentate and fully hydrated monodentate, bidentate, and bridging binding modes. (31)P CSA tensors were calculated from the energy-minimized model compounds using the combined DFT and GIAO methods, resulting in dramatically different tensor values for each binding mode. The anhydrous bidentate chelation mode was found to be a good approximation of the experimental data, an observation that alters the current monodentate paradigm for metal chelation by teichoic acids.

Evidence for a catalytic dyad in the active site of homocitrate synthase from Saccharomyces cerevisiae.

Homocitrate synthase (acetyl-coenzyme A: 2-ketoglutarate C-transferase; E.C. 2.3.3.14) (HCS) catalyzes the condensation of acetyl-CoA (AcCoA) and alpha-ketoglutarate (alpha-KG) to give homocitrate and CoA. Although the structure of an HCS has not been solved, the structure of isopropylmalate synthase (IPMS), a homologue, has been solved (Koon, N., Squire, C. J., and Baker, E. N. (2004) Proc. Natl. Acad. Sci. U.S.A. 101, 8295-8300). Three active site residues in IPMS, Glu-218, His-379, and Tyr-410, were proposed as candidates for catalytic residues involved in deprotonation of the methyl group of AcCoA prior to the Claisen condensation to give homocitrylCoA. All three of the active site residues in IPMS are conserved in the HCS from Saccharomyces cerevisiae. Site-directed mutagenesis has been carried out to probe the role of the homologous residues, Glu-155, His-309, and Tyr-320, in the S. cerevisiae HCS. No detectable activity was observed for the H309A and H309N mutant enzyme, but a slight increase in activity was observed for H309A in the presence of 300 mM imidazole, which is still 1000-fold lower than that of wild type (wt). The E155Q and E155A mutant enzymes exhibited 1000-fold lower activity than wt. The activity of E155A, but not of E155Q, could be partially rescued by formate; a K act of 60 mM with a modest 4-fold maximum activation was observed. In the presence of formate, E155A gives k cat, K AcCoA, and K alpha-KG values of 0.0031 s (-1), 13 muM, and 39 microM, respectively, while a primary kinetic deuterium isotope effect of about 1.4 was obtained on V, with deuterium in the methyl of AcCoA. The pH dependence of k cat for E155A in the presence of formate gave a p K a of 7.9 for a group that must be protonated for optimum activity, similar to that observed for the wt enzyme. However, a partial change was observed on the acid side of the profile, compared to the all or none change observed for wt giving a p K a of about 6.7. The k cat for E155Q decreased at high pH, similar to the wt enzyme, but was pH independent at low pH. The Y320F mutant enzyme only lost 25-fold activity compared to that of the wt, giving k cat, K AcCoA, and K alpha-KG values of 0.039 s (-1), 33 microM, and 140 microM, respectively, and a primary kinetic deuterium isotope effect of 1.3 and 1.8 on V/ K AcCoA and V, respectively; the pH dependence of k cat was similar to that of the wt. These data, combined with a constant pH molecular dynamics simulation study, suggest that a catalytic dyad comprising Glu-155 and His-309 acts to deprotonate the methyl group of AcCoA, while Tyr320 is likely not directly involved in catalysis, but may aid in orienting the reactant and/or the catalytic dyad.

Recent advances in implicit solvent-based methods for biomolecular simulations.

Implicit solvent-based methods play an increasingly important role in molecular modeling of biomolecular structure and dynamics. Recent methodological developments have mainly focused on the extension of the generalized Born (GB) formalism for variable dielectric environments and accurate treatment of nonpolar solvation. Extensive efforts in parameterization of GB models and implicit solvent force fields have enabled ab initio simulation of protein folding to native or near-native structures. Another exciting area that has benefited from the advances in implicit solvent models is the development of constant pH molecular dynamics methods, which have recently been applied to the calculations of protein pK(a) values and the studies of pH-dependent peptide and protein folding.

Linking folding with aggregation in Alzheimer's beta-amyloid peptides.

Growing evidence suggests that the beta-amyloid (Abeta) peptides of Alzheimer's disease are generated in early endosomes and that small oligomers are the principal toxic species. We sought to understand whether and how the solution pH, which is more acidic in endosomes than the extracellular environment, affects the conformational processes of Abeta. Using constant pH molecular dynamics simulations of two model peptides, Abeta(1-28) and Abeta(10-42), we found that the folding landscape of Abeta is strongly modulated by pH and is most favorable for hydrophobically driven aggregation at pH 6. Thus, our theoretical findings substantiate the possibility that Abeta oligomers develop intracellularly before secretion into the extracellular milieu, where they may disrupt synaptic activity or act as seeds for plaque formation.

Exploring atomistic details of pH-dependent peptide folding.

Modeling pH-coupled conformational dynamics allows one to probe many important pH-dependent biological processes, ranging from ATP synthesis, enzyme catalysis, and membrane fusion to protein folding/misfolding and amyloid formation. This work illustrates the strengths and capabilities of continuous constant pH molecular dynamics in exploring pH-dependent conformational transitions in proteins by revisiting an experimentally well studied model protein fragment, the C peptide from ribonuclease A. The simulation data reveal a bell-shaped pH profile for the total helix content, in agreement with experiment, and several pairs of electrostatic interactions that control the relative populations of unfolded and partially folded states of various helical lengths. The latter information greatly complements and extends that attainable by current experimental techniques. The present work paves the way for new and exciting applications, such as the study of pH-dependent molecular mechanism in the formation of amyloid comprising peptides from Alzheimer's and Parkinson's diseases.

Toward the accurate first-principles prediction of ionization equilibria in proteins.

The calculation of pK(a) values for ionizable sites in proteins has been traditionally based on numerical solutions of the Poisson-Boltzmann equation carried out using a high-resolution protein structure. In this paper, we present a method based on continuous constant pH molecular dynamics (CPHMD) simulations, which allows the first-principles description of protein ionization equilibria. Our method utilizes an improved generalized Born implicit solvent model with an approximate Debye-Hückel screening function to account for salt effects and the replica-exchange (REX) protocol for enhanced conformational and protonation state sampling. The accuracy and robustness of the present method are demonstrated by 1 ns REX-CPHMD titration simulations of 10 proteins, which exhibit anomalously large pK(a) shifts for the carboxylate and histidine side chains. The experimental pK(a) values of these proteins are reliably reproduced with a root-mean-square error ranging from 0.6 unit for proteins containing few buried ionizable side chains to 1.0 unit or slightly higher for proteins containing ionizable side chains deeply buried in the core and experiencing strong charge-charge interactions. This unprecedented level of agreement with experimental benchmarks for the de novo calculation of pK(a) values suggests that the CPHMD method is maturing into a practical tool for the quantitative prediction of protein ionization equilibria, and this, in turn, opens a door to atomistic simulations of a wide variety of pH-coupled conformational phenomena in biological macromolecules such as protein folding or misfolding, aggregation, ligand binding, membrane interaction, and catalysis.

Constant pH molecular dynamics with proton tautomerism.

The current article describes a new two-dimensional lambda-dynamics method to include proton tautomerism in continuous constant pH molecular dynamics (CPHMD) simulations. The two-dimensional lambda-dynamics framework is used to devise a tautomeric state titration model for the CPHMD simulations involving carboxyl and histidine residues. Combined with the GBSW implicit solvent model, the new method is tested on titration simulations of blocked histidine and aspartic acid as well as two benchmark proteins, turkey ovomucoid third domain (OMTKY3) and ribonuclease A (RNase A). A detailed analysis of the errors inherent to the CPHMD methodology as well as those due to the underlying solvation model is given. The average absolute error for the computed pKa values in OMTKY3 is 1.0 pK unit. In RNase A the average absolute errors for the carboxyl and histidine residues are 1.6 and 0.6 pK units, respectively. In contrast to the previous work, the new model predicts the correct sign for all the pKa shifts, but one, in the benchmark proteins. The predictions of the tautomeric states of His12 and His48 and the conformational states of His48 and His119 are in agreement with experiment. Based on the simulations of OMTKY3 and RNase A, the current work has demonstrated the capability of the CPHMD technique in revealing pH-coupled conformational dynamics of protein side chains.

Smooth solvation method for d-orbital semiempirical calculations of biological reactions. 1. Implementation.

The present paper describes the extension of a recently developed smooth conductor-like screening model for solvation to a d-orbital semiempirical framework (MNDO/d-SCOSMO) with analytic gradients that can be used for geometry optimizations, transition state searches, and molecular dynamics simulations. The methodology is tested on the potential energy surfaces for separating ions and the dissociative phosphoryl transfer mechanism of methyl phosphate. The convergence behavior of the smooth COSMO method with respect to discretization level is examined and the numerical stability of the energy and gradient are compared to that from conventional COSMO calculations. The present method is further tested in applications to energy minimum and transition state geometry optimizations of neutral and charged metaphosphates, phosphates, and phosphoranes that are models for stationary points in transphosphorylation reaction pathways of enzymes and ribozymes. The results indicate that the smooth COSMO method greatly enhances the stability of quantum mechanical geometry optimization and transition state search calculations that would routinely fail with conventional solvation methods. The present MNDO/d-SCOSMO method has considerable computational advantages over hybrid quantum mechanical/molecular mechanical methods with explicit solvation, and represents a potentially useful tool in the arsenal of multi-scale quantum models used to study biochemical reactions.

Smooth solvation method for d-orbital semiempirical calculations of biological reactions. 2. Application to transphosphorylation thio effects in solution.

Density-functional and semiempirical quantum methods and continuum dielectric and explicit solvation models are applied to study the role of solvation on the stabilization of native and thio-substituted transphosphorylation reactions. Extensive comparison is made between results obtained from the different methods. For the semiempirical methods, explicit solvation was treated using a hybrid quantum mechanical/molecular mechanical (QM/MM) approach and the implicit solvation was treated using a recently developed smooth solvation model implemented into a d-orbital semiempirical framework (MNDO/d-SCOSMO) within CHARMM. The different quantum and solvation methods were applied to the transesterification of a 3'-ribose,5'-methyl phosphodiester that serves as a nonenzymatic model for the self-cleavage reaction catalyzed by the hammerhead and hairpin ribozymes. Thio effects were studied for a double sulfur substitution at the nonbridging phosphoryl oxygen positions. The reaction profiles of both the native and double sulfur-substituted reactions from the MNDO/d-SCOSMO calculations were similar to the QM/MM results and consistent with the experimentally observed trends. These results underscore the need for a d-orbital semiempirical representation for phosphorus and sulfur for the study of experimentally observed thio effects in enzymatic and nonenzymatic phosphoryl transfer reactions. One of the major advantages of the present approach is that it can be applied to model chemical reactions at a significantly lower computational cost than either the density-functional calculations with implicit solvation or the semiempirical QM/MM simulations with explicit solvent.

Quantum descriptors for biological macromolecules from linear-scaling electronic structure methods.

The characterization of electrostatic and chemical properties at the surface of biological macromolecules is of interest in elucidating the fundamental biological structure-function relationships as well as in problems of rational drug design. This paper presents a set of macromolecular quantum descriptors for the characterization of biological macromolecules in solution that can be obtained with modest computational cost from linear-scaling semi-empirical quantum/solvation methods. The descriptors discussed include: solvent-polarized electrostatic surface potential maps, equilibrated atomic charges, Fukui reactivity indices, approximate local hardness maps, and relative proton potentials. These properties are applied to study the conformational dependence of the electrostatic surface potential of the solvated phosphate-binding protein mutant (T141D), the regioselectivity of the zinc finger domains of HIV-1 nucleocapsid (NC) protein, and the order of pKa values of acidic residues in turkey ovomucoid third domain (OMTKY3) and of the zinc-binding residues in the carboxyl terminal zinc finger of NC. In all cases, insight beyond that obtainable from purely classical models is gained and can be used to rationalize the experimental observations. The macromolecular quantum descriptors presented here greatly extend the arsenal of tools for macromolecular characterization and offer promise in applications to modern structure-based drug design.

Insights into the regioselectivity and RNA-binding affinity of HIV-1 nucleocapsid protein from linear-scaling quantum methods.

Human immunodeficiency virus type 1 (HIV-1) nucleocapsid protein (NC) plays several important roles in the viral life-cycle and presents an attractive target for rational drug design. Here, the macromolecular reactivity of NC and its binding to RNA is characterized through determination of electrostatic and chemical descriptors derived from linear-scaling quantum calculations in solution. The computational results offer a rationale for the experimentally observed susceptibility of the Cys49 thiolate toward small-molecule electrophilic agents, and support the recently proposed stepwise protonation mechanism of the C-terminal Zn-coordination complex. The distinctive binding mode of NC to SL2 and SL3 stem-loops of the HIV-1 genomic RNA packaging signal is studied on the basis of protein side-chain contributions to the electrostatic binding energies. These results indicate the importance of several basic residues in the 3(10) helical region and the N-terminal zinc finger, and rationalize the presence of several evolutionarily conserved residues in NC. The combined reactivity and RNA-binding study provides new insights that may contribute toward the structure-based design of anti-HIV therapies.