Energetic basis of molecular recognition in a DNA aptamer

The thermal stability and ligand binding properties of the L-argininamide-binding DNA aptamer (5'-GATCGAAACGTAGCGCCTTCGATC-3') were studied by spectroscopic and calorimetric methods. Differential calorimetric studies showed that the uncomplexed aptamer melted in a two-state reaction with a melting temperature T(m)=50.2+/-0.2 degrees C and a folding enthalpy DeltaH(0)(fold)=-49.0+/-2.1 kcal mol(-1). These values agree with values of T(m)=49.6 degrees C and DeltaH(0)(fold)=-51.2 kcal mol(-1) predicted for a simple hairpin structure. Melting of the uncomplexed aptamer was dependent upon salt concentration, but independent of strand concentration. The T(m) of aptamer melting was found to increase as L-argininamide concentrations increased. Analysis of circular dichroism titration data using a single-site binding model resulted in the determination of a binding free energy DeltaG(0)(bind)=-5.1 kcal mol(-1). Isothermal titration calorimetry studies revealed an exothermic binding reaction with DeltaH(0)(bind)=-8.7 kcal mol(-1). Combination of enthalpy and free energy produce an unfavorable entropy of -TDeltaS(0)=+3.6 kcal mol(-1). A molar heat capacity change of -116 cal mol(-1) K(-1) was determined from calorimetric measurements at four temperatures over the range of 15-40 degrees C. Molecular dynamics simulations were used to explore the structures of the unligated and ligated aptamer structures. From the calculated changes in solvent accessible surface areas of these structures a molar heat capacity change of -125 cal mol(-1) K(-1) was calculated, a value in excellent agreement with the experimental value. The thermodynamic signature, along with the coupled CD spectral changes, suggest that the binding of L-argininamide to its DNA aptamer is an induced-fit process in which the binding of the ligand is thermodynamically coupled to a conformational ordering of the nucleic acid.     

Ab initio modeling of glycosyl torsions and anomeric effects in a model carbohydrate: 2-ethoxy tetrahydropyran

A range of ab initio calculations were carried out on the axial and equatorial anomers of the model carbohydrate 2-ethoxy tetrahydropyran to evaluate the level of theory required to accurately evaluate the glycosyl dihedral angle and the anomeric ratio. Vacuum CCSD(T)/CBS extrapolations at the global minimum yield DeltaE = E(equatorial) - E(axial) = 1.42 kcal/mol. When corrected for solvent (by the IEFPCM model), zero-point vibrations and entropy, DeltaG(298) = 0.49 kcal/mol, in excellent agreement with the experimental value of 0.47 +/- 0.3 kcal/mol. A new additivity scheme, the layered composite method (LCM), yields DeltaE to within 0.1 kcal/mol of the CCSD(T)/CBS result at a fraction of the computer requirements. Anomeric ratios and one-dimensional torsional surfaces generated by LCM and the even more efficient MP2/cc-pVTZ level of theory are in excellent agreement, indicating that the latter is suitable for force-field parameterization of carbohydrates. Hartree-Fock and density functional theory differ from CCSD(T)/CBS for DeltaE by approximately 1 kcal/mol; they show similar deviations in torsional surfaces evaluated from LCM. A comparison of vacuum and solvent-corrected one- and two-dimensional torsional surfaces indicates the equatorial form of 2-ethoxy tetrahydropyran is more sensitive to solvent than the axial.     

Linear free-energy analysis of mercury(II) and cadmium(II) binding to three-stranded coiled coils

Investigators have studied how proteins enforce nonstandard geometries on metal centers to assess the question of how protein structures can define the coordination geometry and binding affinity of an active-site metal cofactor. We have shown that cysteine-substituted versions of the TRI peptide series [AcG-(LKALEEK)(4)G-NH(2)] bind Hg(II) and Cd(II) in geometries that are different from what is normally found with thiol ligands in aqueous solution. A fundamental question has been whether this structural perturbation is due to protein influence or a change in the metal geometry preference. To address this question, we have completed linear free-energy analyses that correlate the association of three-stranded coiled coils in the absence of a metal with the binding affinity of the peptides to the heavy metals, Hg(II) and Cd(II). In this paper, six new members of this family have been synthesized, replacing core leucine residues with smaller and less hydrophobic residues, consequently leading to varying degrees of self-association affinities. At the same time, studies with some smaller and longer sequenced peptides have also been examined. All of these peptides are seen to sequester Hg(II) and Cd(II) in an uncommon trigonal environment. For both metals, the binding is strong with micromolar dissociation constants. For binding of Hg(II) to the peptides, the dissociation constants range from 2.4 x 10(-)(5) M for Baby L12C to 2.5 x 10(-)(9) M for Grand L9C for binding of the third thiolate to a linear Hg(II)(pep)(2) species. The binding of Hg(II) to the peptide Grand L9C is similar in energetics for metal binding in the metalloregulatory protein, mercury responsive (merR), displaying approximately 50% trigonal Hg(II) formation at nanomolar metal concentrations. Approximately, 11 kcal/mol of the Hg(II)(Grand L9C)(3)(-) stability is due to peptide interactions, whereas only 1-4 kcal/mol stabilization results from Hg(II)(RS)(2) binding the third thiolate ligand. This further validates the hypothesis that the favorable tertiary interactions in protein systems such as merR go a long way in stabilizing nonnatural coordination environments in biological systems. Similarly, for the binding of Cd(II) to the TRI family, the dissociation constants range from 1.3 x 10(-)(6) M for Baby L9C to 8.3 x 10(-)(9) M for TRI L9C, showing a similar nature of stable aggregate formation.     

Reaction pathway and free energy profiles for butyrylcholinesterase-catalyzed hydrolysis of acetylthiocholine

The catalytic mechanism for butyrylcholineserase (BChE)-catalyzed hydrolysis of acetylthiocholine (ATCh) has been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy (QM/MM-FE) calculations on both acylation and deacylation of BChE. Additional quantum mechanical (QM) calculations have been carried out, along with the QM/MM-FE calculations, to understand the known substrate activation effect on the enzymatic hydrolysis of ATCh. It has been shown that the acylation of BChE with ATCh consists of two reaction steps including the nucleophilic attack on the carbonyl carbon of ATCh and the dissociation of thiocholine ester. The deacylation stage includes nucleophilic attack of a water molecule on the carboxyl carbon of substrate and dissociation between the carboxyl carbon of substrate and hydroxyl oxygen of Ser198 side chain. QM/MM-FE calculation results reveal that the acylation of BChE is rate-determining. It has also been demonstrated that an additional substrate molecule binding to the peripheral anionic site (PAS) of BChE is responsible for the substrate activation effect. In the presence of this additional substrate molecule at PAS, the calculated free energy barrier for the acylation stage (rate-determining step) is decreased by ~1.7 kcal/mol. All of our computational predictions are consistent with available experimental kinetic data. The overall free energy barriers calculated for BChE-catalyzed hydrolysis of ATCh at regular hydrolysis phase and substrate activation phase are ~13.6 and ~11.9 kcal/mol, respectively, which are in reasonable agreement with the corresponding experimentally derived activation free energies of 14.0 kcal/mol (for regular hydrolysis phase) and 13.5 kcal/mol (for substrate activation phase).     

A density functional investigation of the extradiol cleavage mechanism in non-heme iron catechol dioxygenases

The mechanism for extradiol cleavage in non-heme iron catechol dioxygenase was modelled theoretically via density functional theory. Based on the Fe(II)-His,His,Glu motif observed in enzymes, an active site model complex, [Fe(acetate)(imidazole)(2)(catecholate)(O(2))](-), was optimized for states with six, four and two unpaired electrons (U6, U4 and U2, respectively). The transfer of the terminal atom of the coordinated dioxygen leading to "ferryl" Fe=O intermediates spontaneously generates an extradiol epoxide. The computed barriers range from 19 kcal mol(-1) on the U6 surface to approximately 25 kcal mol(-1) on the U4 surface, with overall reaction energies of +11.6, 6.3 and 7.1 kcal mol(-1) for U6, U4 and U2, respectively. The calculations for a protonated process reveal the terminal oxygen of O(2) to be the thermodynamically favoured site but subsequent oxygen transfer to the catechol has a barrier of approximately 30-40 kcal mol(-1), depending on the spin state. Instead, protonating the acetate group gives a slightly higher energy species but a subsequent barrier on the U4 surface of only 7 kcal mol(-1) relative to the hydroperoxide complex. The overall exoergicity increases to 13 kcal mol(-1). The favoured proton-assisted pathway does not involve significant radical character and has features reminiscent of a Criegee rearrangement which involves the participation of the aromatic ring pi-orbitals in the formation of the new carbon-oxygen bond. The subsequent collapse of the epoxide, attack by the coordinated hydroxide and final product formation proceeds with an overall exoergicity of approximately 75 kcal mol(-1) on the U4 surface.     

Mössbauer characterization of an unusual high-spin side-on peroxo-Fe3+ species in the active site of superoxide reductase from Desulfoarculus Baarsii. Density functional calculations on related models

Superoxide reductase (SOR) is an Fe protein that catalyzes the reduction of superoxide to give H(2)O(2). Recently, the mutation of the Glu47 residue into alanine (E47A) in the active site of SOR from Desulfoarculus baarsii has allowed the stabilization of an iron-peroxo species when quickly reacted with H(2)O(2) [Mathé et al. (2002) J. Am. Chem. Soc. 124, 4966-4967]. To further investigate this non-heme peroxo-iron species, we have carried out a M?ssbauer study of the (57)Fe-enriched E47A SOR from D. baarsii reacted quickly with H(2)O(2). Considering the M?ssbauer data, we conclude, in conjunction with the other spectroscopic data available and with the results of density functional calculations on related models, that this species corresponds to a high-spin side-on peroxo-Fe(3+) complex. This is one of the first examples of such a species in a biological system for which M?ssbauer parameters are now available: delta(/Fe) = 0.54 (1) mm/s, DeltaE(Q) = -0.80 (5) mm/s, and the asymmetry parameter eta = 0.60 (5) mm/s. The M?ssbauer and spin Hamiltonian parameters have been evaluated on a model from the side-on peroxo complex (model 2) issued from the oxidized iron center in SOR from Pyrococcus furiosus, for which structural data are available in the literature [Yeh et al. (2000) Biochemistry 39, 2499-2508]. For comparison, similar calculations have been carried out on a model derived from 2 (model 3), where the [CH(3)-S](1)(-) group has been replaced by the neutral [NH(3)](0) group [Neese and Solomon (1998) J. Am. Chem. Soc. 120, 12829-12848]. Both models 2 and 3 contain a formally high-spin Fe(3+) ion (i.e., with empty minority spin orbitals). We found, however, a significant fraction ( approximately 0.6 for 2, approximately 0.8 for 3) of spin (equivalently charge) spread over two occupied (minority spin) orbitals. The quadrupole splitting value for 2 is found to be negative and matches quite well the experimental value. The computed quadrupole tensors are rhombic in the case of 2 and axial in the case of 3. This difference originates directly from the presence of the thiolate ligand in 2. A correlation between experimental isomer shifts for Fe(3+) mononuclear complexes with computed electron densities at the iron nucleus has been built and used to evaluate the isomer shift values for 2 and 3 (0.56 and 0.63 mm/s, respectively). A significant increase of isomer shift value is found upon going from a methylthiolate to a nitrogen ligand for the Fe(3+) ion, consistent with covalency effects due to the presence of the axial thiolate ligand. Considering that the isomer shift value for 3 is likely to be in the 0.61-0.65 mm/s range [Horner et al. (2002) Eur. J. Inorg. Chem., 3278-3283], the isomer shift value for a high-spin eta(2)-O(2) Fe(3+) complex with an axial thiolate group can be estimated to be in the 0.54-0.58 mm/s range. The occurrence of a side-on peroxo intermediate in SOR is discussed in relation to the recent data published for a side-on peroxo-Fe(3+) species in another biological system [Karlsson et al. (2003) Science 299, 1039-1042].     

The enthalpy change for the reduction of nicotinamide–adenine dinucleotide

The heat of the reaction NAD(+)+propan-2-ol=NADH+acetone+H(+) was determined to be 42.5+/-0.6kJ/mol (10.17+/-0.15kcal/mol) from equilibrium measurements at 9-42 degrees C catalysed by yeast alcohol dehydrogenase. With the aid of thermochemical data for acetone and propan-2-ol the values of DeltaH=-29.2kJ/mol (-6.99kcal/mol) and DeltaG(0)=22.1kJ/mol (5.28kcal/mol) are derived for the reduction of NAD (NAD(+)+H(2)=NADH+H(+)). These values are consistent with analogous but less accurate data for the ethanol-acetaldehyde reaction. Thermodynamic data for the reduction of NAD and NADP are summarized.     

Volume and enthalpy profiles of CO binding to Fe(II) tetrakis-(4-sulfonatophenyl)porphyrin.

The focus of this study is to examine volume and enthalpy profiles of ligand binding associated with CO-Fe(II) tetrakis-(4-sulfonato phenyl)-porphyrin (COFe(II)4SP) in aqueous solution. Temperature dependent photothermal beam deflection was employed to probe the overall enthalpy and volume changes associated with CO-photolysis and recombination. The analysis demonstrates that ligand recombination occurs with a pseudo first order rate constant of (2.5+/-0.2)x10(4) s(-1) (at 25 degrees C) with a corresponding volume decrease of 6+/-1 ml/mol. The activation enthalpy (DeltaH(double dagger)) and volume (DeltaV(double dagger)) change for CO recombination (determined from temperature/pressure dependent transient absorption spectroscopy) are found to be 3.9 kcal/mol and 8.2 ml/mol, respectively. These data are consistent with a mechanism in which photolysis yields a five-coordinate high spin (H(2)O)Fe(II)4SP complex that recombines in a single step to form the low spin (CO)(H(2)O)Fe(II)4SP complex. Base elimination, often associated with CO photolysis from hemes, is not observed in this system. The overall volume changes suggest a transition state with significant high spin character. Furthermore, these results demonstrate the utility of coupling photothermal techniques with variable pressure/temperature transient absorption spectroscopy to probe heme reaction dynamics.     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

Temperature dependence of collagen fluorescence.

Dermal collagens have several fluorescent moieties in the UV and visible spectral regions that may serve as molecular probes of collagen. We studied the temperature-dependence of a commercial calf skin collagen and acid-extracted Skh-1 hairless mouse collagen at temperatures from 9 degrees C to 60 degrees C for excitation/emission wavelengths 270/305 nm (tyrosine), 270/360 nm (excimer-like aggregated species), 325/400 nm (dityrosine) and 370/450 nm (glycation adduct). L-tyrosine (1 x 10(-5) M in 0.5 M HOAc) acted as a "reference compound" devoid of any collagen structural effects. In general, the fluorescence efficiency of these fluorophores decreases with increasing temperature. Assuming that rate constant for fluorescence deactivation has the form k(d)(T) = k(d) degrees exp (-DeltaE/RT), an Arrhenius plot of log[(1/Phi) - 1] vs. 1/T affords a straight line whose (negative) slope is proportional to the activation energy, DeltaE, of the radiationless process(es) that compete with fluorescence. Because it is difficult to accurately measure Phi(f) for collagen-bound fluorophores, we derived an approximate formula for an activation parameter, DeltaE*, evaluated from an Arrhenius-like plot of log 1/I(N)vs. 1/T, (1/I(N)vs. is the reciprocal normalized fluorescence intensity). Tyrosine in dilute solution affords a linear Arrhenius plot in both of the above cases. Using the known value of Phi(f) = 0.21 for free tyrosine at room temperature, we determined that DeltaE* is accurate to approximately 25% in the present instance. Collagen curves are non-linear, but they are quasi-linear below approximately 20 degrees C, where the helical form predominates. Values of DeltaE* determined from the data at T < 20 degrees C ranged from 6.2-8.4 kJ mol(-1) (1.5-2.0 kcal mol(-1)) for mouse collagen and 10.3-11.4 kJ mol(-1) (2.5-2.7 kcal mol(-1)) for calf skin collagen, consistent with collisional deactivation of the fluorescent state via thermally enhanced molecular vibrations and rotations. Above 20 degrees C, log 1/I(N)vs. 1/T plots from Skh-1 hairless mouse collagen are concave-downward, suggesting that fluorescence deactivation from the denatured coil has a significant temperature-independent component. For calf skin collagen, these plots are concave-upward, suggesting an increase in activation energy above Tm. These results suggest that collagen backbone and supramolecular structure can influence the temperature dependence of the bound fluorophores, indicating the future possibility of using activation data as a probe of supramolecular structure and conformation.     

Electron transfer kinetics and mechanistic study of the thionicotinamide coordinated to the pentacyanoferrate(III)/(II) complexes: a model system for the in vitro activation of thioamides anti-tuberculosis drugs

The mechanism of activation thioamide-pyridine anti-tuberculosis prodrugs is poorly described in the literature. It has recently been shown that ethionamide, an important component of second-line therapy for the treatment of multi-drug-resistant tuberculosis, is activated through an enzymatic electron transfer (ET) reaction. In an attempt to shed light on the activation of thioamide drugs, we have mimicked a redox process involving the thionicotinamide (thio) ligand, investigating its reactivity through coordination to the redox reversible [Fe(III/II)(CN)(5)(H(2)O)](2-/3-) metal center. The reaction of the Fe(III) complex with thionicotinamide leads to the ligand conversion to the 3-cyanopyridine species coordinated to a Fe(II) metal center. The rate constant, k(et)=10 s(-1), was determined for this intra-molecular ET reaction. A kinetic study for the cross-reaction of thionicotinamide and [Fe(CN)(6)](3-) was also carried out. The oxidation of thionicotinamide by [Fe(CN)(6)](3-) leads to formation of mainly 3-cyanopyridine and [Fe(CN)(6)](4-) with a k(et)=(5.38+/-0.03) M(-1)s(-1) at 25 degrees C, pH 12.0. The rate of this reaction is strongly dependent on pH due to an acid-base equilibrium related to the deprotonation of the R-SH functional group of the imidothiol form of thionicotinamide. The kinetic results reinforced the assignment of an intra-molecular mechanism for the ET reaction of [Fe(III)(CN)(5)(H(2)O)](2-) and the thioamide ligand. These results can be valuable for the design of new thiocarbonyl-containing drugs against resistant strains of Mycobacterium tuberculosis by a self-activating mechanism.     

The contribution of heme propionate groups to the conformational dynamics associated with CO photodissociation from horse heart myoglobin

Photoacoustic calorimetry and transient absorption spectroscopy were used to study conformational dynamics associated with CO photodissociation from horse heart myoglobin (Mb) reconstituted with either Fe protoporphyrin IX dimethylester (FePPDME), Fe octaethylporphyrin (FeOEP), or with native Fe protoporphyrin IX (FePPIX). The volume and enthalpy changes associated with the Fe-CO bond dissociation and formation of a transient deoxyMb intermediate for the reconstituted Mbs were found to be similar to those determined for native Mb (DeltaV1 = -2.5+/-0.6 ml mol(-1) and DeltaH1 = 8.1+/-3.0 kcal mol(-1)). The replacement of FePPIX by FeOEP significantly alters the conformational dynamics associated with CO release from protein. Ligand escape from FeOEP reconstituted Mb was determined to be roughly a factor of two faster (tau=330 ns) relative to native protein (tau=700 ns) and accompanying reaction volume and enthalpy changes were also found to be smaller (DeltaV2 = 5.4+/-2.5 ml mol(-1) and DeltaH2 = 0.7+/-2.2 kcal mol(-1)) than those for native Mb (DeltaV2 = 14.3+/-0.8 ml mol(-1) and DeltaH2 = 7.8+/-3.5 kcal mol(-1)). On the other hand, volume and enthalpy changes for CO release from FePPIX or FePPDME reconstituted Mb were nearly identical to those of the native protein. These results suggest that the hydrogen bonding network between heme propionate groups and nearby amino acid residues likely play an important role in regulating ligand diffusion through protein matrix. Disruption of this network leads to a partially open conformation of protein with less restricted ligand access to the heme binding pocket.     

Redox thermodynamics of mutant forms of the rubredoxin from Clostridium pasteurianum: identification of a stable FeIII(S-Cys)3(OH) centre in the C6S mutant

Redox thermodynamic data provide a detailed insight into control of the reduction potential E degrees' of the [Fe(S-Cys)4] site in rubredoxin. Mutant forms were studied in which specific structural changes were made in both the primary and secondary coordination spheres. Those changes have been probed by resonance Raman spectroscopy. The decrease of approximately 200 mV in E degrees' observed for the [Fe(S-Cys)3(O-Ser)]-/2- couples in the surface ligand mutants C9S and C42S is essentially enthalpic in origin and associated with the substitution of ligand thiolate by ligand olate. However, the pH dependence of the potentials below characteristic pKa(red) approximately equals 7 is an entropic contribution, plausibly associated with increased conformational flexibility induced by a longer Fe(II)-O(H)-Ser bond in the reduced form. The presence of a second surface Ser ligand in the new double mutant protein C9S/C42S affects the enthalpic term primarily for pH>pKa(red) > or = 9.3, but for pHpKa approximately 9: [Fe(III)(S-Cys)3(OH)]- + e- --> [Fe(II)(S-Cys)3(OH)]2-. pH [Fe(II)(S-Cys)3(OH2)]-.     

Reactivity pathways for nitric oxide and nitrosonium with iron complexes in biologically relevant sulfur coordination spheres

The interaction of nitric oxide (NO) with iron-sulfur cluster proteins results in the formation of dinitrosyl iron complexes (DNICs) coordinated by cysteine residues from the peptide backbone or with low molecular weight sulfur-containing molecules like glutathione. Such DNICs are among the modes available in biology to store, transport, and deliver NO to its relevant targets. In order to elucidate the fundamental chemistry underlying the formation of DNICs and to characterize possible intermediates in the process, we have investigated the interaction of NO (g) and NO(+) with iron-sulfur complexes having the formula [Fe(SR)(4)](2-), where R=(t)Bu, Ph, or benzyl, chosen to mimic sulfur-rich iron sites in biology. The reaction of NO (g) with [Fe(S(t)Bu)(4)](2-) or [Fe(SBz)(4)](2-) cleanly affords the mononitrosyl complexes (MNICs), [Fe(S(t)Bu)(3)(NO)](-) (1) and [Fe(SBz)(3)(NO)](-) (3), respectively, by ligand displacement. Mononitrosyl species of this kind were previously unknown. These complexes further react with NO (g) to generate the corresponding DNICs, [Fe(SPh)(2)(NO)(2)](-) (4) and [Fe(SBz)(2)(NO)(2)](-) (5), with concomitant reductive elimination of the coordinated thiolate donors. Reaction of [Fe(SR)(4)](2-) complexes with NO(+) proceeds by a different pathway to yield the corresponding dinitrosyl S-bridged Roussin red ester complexes, [Fe(2)(mu-S(t)Bu)(2)(NO)(4)] (2), [Fe(2)(mu-SPh)(2)(NO)(4)] (7) and [Fe(2)(mu-SBz)(2)(NO)(4)] (8). The NO/NO(+) reactivity of an Fe(II) complex with a mixed nitrogen/sulfur coordination sphere was also investigated. The DNIC and red ester species, [Fe(S-o-NH(2)C(6)H(4))(2)(NO)(2)](-) (6) and [Fe(2)(mu-S-o-NH(2)C(6)H(4))(2)(NO)(4)] (9), were generated. The structures of 8 and 9 were verified by X-ray crystallography. The MNIC complex 1 can efficiently deliver NO to iron-porphyrin complexes like [Fe(TPP)Cl], a reaction that is aided by light. Removal of the coordinated NO ligand of 1 by photolysis and addition of elemental sulfur generates higher nuclearity Fe/S clusters.     

Carbinolamine formation and dehydration in a DNA repair enzyme active site

In order to suggest detailed mechanistic hypotheses for the formation and dehydration of a key carbinolamine intermediate in the T4 pyrimidine dimer glycosylase (T4PDG) reaction, we have investigated these reactions using steered molecular dynamics with a coupled quantum mechanics-molecular mechanics potential (QM/MM). We carried out simulations of DNA abasic site carbinolamine formation with and without a water molecule restrained to remain within the active site quantum region. We recovered potentials of mean force (PMF) from thirty replicate reaction trajectories using Jarzynski averaging. We demonstrated feasible pathways involving water, as well as those independent of water participation. The water-independent enzyme-catalyzed reaction had a bias-corrected Jarzynski-average barrier height of approximately (6.5 kcal mol(-1) (27.2 kJ mol(-1)) for the carbinolamine formation reaction and 44.5 kcal mol(-1) (186 kJ mol(-1)) for the reverse reaction at this level of representation. When the proton transfer was facilitated with an intrinsic quantum water, the barrier height was approximately 15 kcal mol(-1) (62.8 kJ mol(-1)) in the forward (formation) reaction and 19 kcal mol(-1) (79.5 kJ mol(-1)) for the reverse. In addition, two modes of unsteered (free dynamics) carbinolamine dehydration were observed: in one, the quantum water participated as an intermediate proton transfer species, and in the other, the active site protonated glutamate hydrogen was directly transferred to the carbinolamine oxygen. Water-independent unforced proton transfer from the protonated active site glutamate carboxyl to the unprotonated N-terminal amine was also observed. In summary, complex proton transfer events, some involving water intermediates, were studied in QM/MM simulations of T4PDG bound to a DNA abasic site. Imine carbinolamine formation was characterized using steered QM/MM molecular dynamics. Dehydration of the carbinolamine intermediate to form the final imine product was observed in free, unsteered, QM/MM dynamics simulations, as was unforced acid-base transfer between the active site carboxylate and the N-terminal amine.     

Transition of hemoglobin between two tertiary conformations: determination of equilibrium and thermodynamic parameters from the reaction of 5,5'-dithiobis(2-nitrobenzoate) with the CysF9[93]beta sulfhydryl group

The equilibrium constant of the reaction of 5,5'-dithiobis(2-nitrobenzoate) with the CysF9[93]beta sulfhydryl group of hemoglobin decreases by 2 to 3 orders of magnitude between pH 5.6 and 9. The reaction is coupled to the ionizations of two groups on the protein. At 25 degrees C one group has a pK(a) of 5.31+/-0.2 when hemoglobin is in its (tertiary) r conformation, typified by the thiolate anion form of CysF9[93]beta; this changes to 7.73+/-0.4 in the (tertiary) t conformation, typified by the mixed disulfide form of the sulfhydryl. The second group ionizes with a pK(a) of 7.11+/-0.4 in the r conformation; this changes to 8.38+/-0.2 in the t conformation. K(rt), the equilibrium constant for the r<-->t isomerization process, is 0.22+/-0.06. The standard enthalpy and entropy changes for the isomerization are DeltaH(o)(rt)=24.2 kJ mol(-1) and DeltaS(o)(rt)=68.8 JK(-1)mol(-1), respectively.     

Increasing the conformational stability by replacement of heme axial ligand in c-type cytochrome

To investigate the role of the heme axial ligand in the conformational stability of c-type cytochrome, we constructed M58C and M58H mutants of the red alga Porphyra yezoensis cytochrome c(6) in which the sixth heme iron ligand (Met58) was replaced with Cys and His residues, respectively. The Gibbs free energy change for unfolding of the M58H mutant in water (DeltaG degrees (unf)=1.48 kcal/mol) was lower than that of the wild-type (2.43 kcal/mol), possibly due to the steric effects of the mutation on the apoprotein structure. On the other hand, the M58C mutant exhibited a DeltaG degrees (unf) of 5.45 kcal/mol, a significant increase by 3.02 kcal/mol compared with that of wild-type. This increase was possibly responsible for the sixth heme axial bond of M58C mutant being more stable than that of wild-type according to the heme-bound denaturation curve. Based on these observations, we propose that the sixth heme axial ligand is an important key to determine the conformational stability of c-type cytochromes, and the sixth Cys heme ligand will give stabilizing effects.     

Catalytic mechanism and product specificity of rubisco large subunit methyltransferase: QM/MM and MD investigations

Molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations have been carried out in an investigation of Rubisco large subunit methyltransferase (LSMT). It was found that the appearance of a water channel is required for the stepwise methylation by S-adenosylmethionine (AdoMet). The water channel appears in the presence of AdoMet (LSMT.Lys-NH3+.AdoMet), but is not present immediately after methyl transfer (LSMT.Lys-N(Me)H2+.AdoHcy). The water channel allows proton dissociation from both LSMT.AdoMet.Lys-NH3+ and LSMT.AdoMet.Lys-N(Me)H2+. The water channel does not appear for proton dissociation from LSMT.AdoMet.Lys-N(Me)2H+, and a third methyl transfer does not occur. By QM/MM, the calculated free energy barrier of the first methyl transfer reaction catalyzed by LSMT (Lys-NH2 + AdoMet --> Lys-N(Me)H2+ + AdoHcy) is DeltaG++ = 22.8 +/- 3.3 kcal/mol. This DeltaG++ is in remarkable agreement with the value 23.0 kcal/mol calculated from the experimental rate constant (6.2 x 10-5 s-1). The calculated DeltaG++ of the second methyl transfer reaction (AdoMet + Lys-N(Me)H --> AdoHcy + Lys-N(Me)2H+) at the QM/MM level is 20.5 +/- 3.6 kcal/mol, which is in agreement with the value 22.0 kcal/mol calculated from the experimental rate constant (2.5 x 10-4 s-1). The third methyl transfer (Lys-N(Me)2 + AdoMet --> Lys-N(Me)3+ + AdoHcy) is associated with an allowed DeltaG++ of 25.9 +/- 3.2 kcal/mol. However, this reaction does not occur because a water channel does not form to allow the proton dissociation of Lys-N(Me)2H+. Future studies will determine whether the product specificity of lysine (mono, di, and tri) methyltransferases is determined by the formation of water channels.