Search PubMed⌕ Search

Biomedical subjects

Walter Thiel

Publications and source records attributed to Walter Thiel.

At least 37 records · Page 2Linked to original sources
PubMed2005

Two states and two more in the mechanisms of hydroxylation and epoxidation by cytochrome P450.

Past studies have shown that oxidation reactions by P450 Compound I (Cpd I) can be described by two competing quartet and doublet spin states, which possess three unpaired electrons, hence tri-radicals. One electron excitation from the delta orbital to sigma* xy generates two states that possess five unpaired electrons, so-called penta-radicals, in sextet and quartet situations, and which were shown by theory to lie only approximately 12-14 kcal/mol higher in energy than the tri-radical ground states (ref 7). The present study focuses on the C-H hydroxylation and C=C epoxidation of propene by these penta-radical states. It is shown that the initial energy differences, between the penta-radical and tri-radical states, diminish along the reaction pathway, due to the favorable and cumulative exchange stabilization of the more open-shell species. Furthermore, theory suggests that hydrogen bonding to the thiolate ligand, and general polarity of the environment, reduce these gaps further, thereby making the penta-radical states accessible to ground-state reactivity. The interconversion between the tri-radical and penta-radical states along the reaction coordinate will depend on the dynamics of spin-flips and energy barriers between the states. Especially interesting should be the region of the reaction intermediates; for both epoxidation and hydroxylation, this region is typified by a dense manifold of spin states and electromeric states (that differ by the oxidation state of iron), such that the total reactivity would be expected to reflect the interplay of these states, giving rise to multistate reactivity.

Alkenes↗
PubMed2005

Preparation of [5,6]- and [6,6]-oxahomofullerene derivatives and their interconversion by Lewis acid assisted reactions of fullerene mixed peroxides.

[60]Fullerene mixed peroxides C60(O)(OOtBu)4 exhibit chemo- and regioselective reactions under mild conditions. The epoxy moiety is opened by ferric chloride to form vicinal hydroxy chloride C60Cl(OH)(OOtBu)4. BF3 is also effective in opening the epoxy moiety. The O-O bond of the fullerene mixed peroxide is cleaved by aluminum chloride to form both [5,6]- and [6,6]-fullerene hemiketals (oxohomo[60]fullerenes). A Hock-type rearrangement is proposed for the formation of the hemiketals, in which a fullerene C-C bond is cleaved. Lewis acids and/or visible light can initiate isomerization of the hemiketal isomers. Single-crystal X-ray analysis and theoretical calculations confirmed the results.

Journal Article↗
PubMed2005

Unusual temperature effects in propene polymerization using stereorigid zirconocene catalysts.

Free Car-Parrinello molecular dynamics (CPMD) simulations of four diastereomers of the zirconium-propene complexes [{iPr(3-iPr-CpFlu)}ZriBu(C3H6)]+ (Cp=cyclopentadienyl; Flu=fluorenyl) provide valuable insight into the mechanism and stereocontrol of propene polymerization with stereorigid metallocenes. Spontaneous insertion of propene into the zirconium-isobutyl bond is not observed, and propene is found to be weakly bound and to rotate relatively freely around the C--C bond to be formed. Large-amplitude rotation of the isopropyl substituent around the Cp--iPr bond may play a role in triggering dissociation of propene. Three of the four diastereomers eliminate propene during the course of the simulations, which makes dissociation the dominating event on a 20-ps timescale. The CPMD simulations thus support the validity of the assumption, fundamental to statistical propagation models, that each insertion is independent of the preceding insertions. Using insertion barriers from static density functional calculations, the statistical model predicts the polypropene microstructure in good agreement with experiment at low polymerization temperatures for the catalysts {iPr(3-R-CpFlu)}ZrCl2 (R=H, iPr, tBu). The predictions become less accurate at higher temperatures, probably due to the onset of the competing back-skip reaction, which is not included in the model.

Alkenes↗
PubMed2005

The palladium-catalyzed cross-coupling reaction of carboxylic anhydrides with arylboronic acids: a DFT study.

The mechanism of the cross-coupling of phenylboronic acid with acetic anhydride, a viable model of the widely used Suzuki reaction, has been studied by DFT calculations at the BP86/6-31G level of theory. Two alternative catalytic cycles have been investigated, one starting from a neutral Pd(0)L(2) complex, the other from an anionic "Jutand-type" [Pd(0)L(2)X](-) species. The reaction profiles are in good agreement with the experimental findings, as both pathways require only moderate activation energies. Both pathways are dominated by cis-configured square-planar palladium(II)diphosphine intermediates. Despite careful investigations, we did not find in this model reaction any evidence for five-coordinate palladium(II) intermediates, which are commonly believed to cause the profound effects of counterions in palladium-catalyzed transformations. Instead, our calculations suggest that the higher catalytic activity of anionic complexes, such as [Pd(PMe(3))(2)OAc](-), may arise from their stronger ability to coordinate to carbon electrophiles. The transmetalation sequence is the same for both catalytic cycles, involving the dissociation of one phosphine ligand from the palladium. In the decisive transition state, in which the phenyl group is transferred from boron to palladium, the acetate base is found to be in a bridging coordination between these two atoms.

Journal Article↗
PubMed2005

Why are BINOL-based monophosphites such efficient ligands in Rh-catalyzed asymmetric olefin hydrogenation?

Whereas recent synthetic studies concerning Rh-catalyzed olefin hydrogenation based on BINOL-derived monodentate phosphites have resulted in an efficient and economically attractive preparative method, very little is known concerning the source of the unexpectedly high levels of enantioselectivity (ee often 90-99%). The present mechanistic study, which includes the NMR characterization of the precatalysts, kinetic measurements with focus on nonlinear effects, and DFT calculations, constitutes a first step in understanding this hydrogenation system. The two most important features which have emerged from these efforts are the following: (1) two monodentate P-ligands are attached to rhodium, and (2) the lock-and-key mechanism holds, in which the thermodynamics of Rh/olefin complexation with formation of the major and minor diastereomeric intermediates dictates the stereochemical outcome. The major diastereomer leads to the favored enantiomeric product, which is opposite to the state of affairs in classical Rh-catalyzed olefin hydrogenation based on chiral chelating diphosphines (anti lock-and-key mechanism as proposed by Halpern).

Journal Article↗
PubMed2005

Ring-closing olefin metathesis on ruthenium carbene complexes: model DFT study of stereochemistry.

Ring-closing metathesis (RCM) is the key step in a recently reported synthesis of salicylihalamide and related model compounds. Experimentally, the stereochemistry of the resulting cycloolefin (cis/trans) depends strongly on the substituents that are present in the diene substrate. To gain insight into the factors that govern the observed stereochemistry, density functional theory (DFT) calculations have been carried out for a simplified dichloro(2-propylidene)(imidazole-2-ylidene)ruthenium catalyst I, as well as for the real catalyst II with two mesityl substituents on the imidazole ring. Four model substrates are considered, which are closely related to the systems studied experimentally, and in each case, two pathways A and B are possible since the RCM reaction can be initiated by coordination of either of the two diene double bonds to the metal center. The first metathesis yields a carbene intermediate, which can then undergo a second metathesis by ring closure, metallacycle formation, and metallacycle cleavage to give the final cycloolefin complex. According to the DFT calculations, the stereochemistry is always determined in the second metathesis reaction, but the rate-determining step may be different for different catalysts, substrates, and pathways. The ancillary N-heterocyclic carbene ligand lies in the Ru-Cl-Cl plane in the simplified catalyst I, but is perpendicular to it in the real catalyst II, and this affects the relative energies of the relevant intermediates and transition states. Likewise, the introduction of methyl substituents in the diene substrates influences these relative energies appreciably. Good agreement with the experimentally observed stereochemistry is only found when using the real catalyst II and the largest model substrates in the DFT calculations.

Journal Article↗
PubMed2005

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.

Algorithms↗
PubMed2005

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.

Algorithms↗
PubMed2005

Toward identification of the compound I reactive intermediate in cytochrome P450 chemistry: a QM/MM study of its EPR and Mössbauer parameters.

Quantum mechanical/molecular mechanical (QM/MM) methods have been used in conjunction with density functional theory (DFT) and correlated ab initio methods to predict the electron paramagnetic resonance (EPR) and Mossbauer (MB) properties of Compound I in P450(cam). For calibration purposes, a small Fe(IV)-oxo complex [Fe(O)(NH(3))(4)(H(2)O)](2+) was studied. The (3)A(2) and (5)A(1) states (in C(4)(v)() symmetry) are found to be within 0.1-0.2 eV. The large zero-field splitting (ZFS) of the (FeO)(2+) unit in the (3)A(2) state arises from spin-orbit coupling with the low-lying quintet and singlet states. The intrinsic g-anisotropy is very small. The spectroscopic properties of the model complex [Fe(O)(TMC)(CH(3)CN)](2+) (TMC = 1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane) are well reproduced by theory. In the model complexes [Fe(O)(TMP)(X)](+) (TMP = tetramesitylporphyrin, X = nothing or H(2)O) the computations again account for the observed spectroscopic properties and predict that the coupling of the (5)A(1) state of the (FeO)(2+) unit to the porphyrin radical leads to a low-lying sextet/quartet manifold approximately 12 kcal/mol above the quartet ground state. The calculations on cytochrome P450(cam), with and without the simulation of the protein environment by point charges, predict a small antiferromagnetic coupling (J approximately -13 to -16 cm(-)(1); H(HDvV) = - 2JS(A)S(B)) and a large ZFS > 15 cm(-)(1) (with E/D approximately 1/3) which will compete with the exchange coupling. This leads to three Kramers doublets of mixed multiplicity which are all populated at room temperature and may therefore contribute to the observed reactivity. The MB and ligand hyperfine couplings ((14)N, (1)H) are fairly sensitive to the protein environment which controls the spin density distribution between the porphyrin ring and the axial cysteinate ligand.

Cytochrome P-450 Enzyme System↗
PubMed2005

Dipole moment and rovibrational intensities in the electronic ground state of NH3: bridging the gap between ab initio theory and spectroscopic experiment.

We report theoretical values for the transition moments of an extensive set of vibrational bands in the electronic ground state of (14)NH(3). For selected bands, we have further made detailed simulations of the rotational structure. The calculations are carried out by means of recently developed computational procedures for describing the nuclear motion and are based on a high-level ab initio potential energy surface, and high-level dipole moment surfaces, for the electronic ground state of NH(3). The reported theoretical intensity values are compared to, and found to agree very well with, corresponding experimental results. It is believed that the computational method, in conjunction with high-quality ab initio potential energy and dipole moment surfaces, can simulate rotation-vibration spectra of XY(3) pyramidal molecules prior to observation with sufficient accuracy to facilitate the observation of these spectra. By degrading the accuracy of selected elements of the calculations, we have also investigated the influence of customary approximations on the computed intensity values.

Journal Article↗
PubMed2005

Combined quantum mechanical/molecular mechanical study on the pentacoordinated ferric and ferrous cytochrome P450cam complexes.

The pentacoordinated ferric and ferrous cytochrome P450(cam) complexes have been investigated by combined quantum mechanical/molecular mechanical (QM/MM) calculations in the presence of a protein/solvent environment and by QM calculations on the isolated QM regions with use of density functional theory. The B3LYP functional has been found more reliable than the BLYP and BHLYP functionals for estimating the relative state energies. The B3LYP/CHARMM calculations with an all-electron basis set for iron give high-spin ground states for the title complexes, in agreement with experiment. The comparison of the B3LYP/CHARMM results of the entire protein system with the B3LYP calculations on the naked QM regions shows that the amount of stabilization by the protein environment is largest for the intermediate-spin states, followed by the high-spin states of the complexes. The calculation of Mössbauer parameters in the presence of the enzyme environment confirms the double occupation of the d(xz) orbital in the quintet spin state of the ferrous complex, consistent with the computed QM/MM energies in the enzyme environment, while the d(x)2(-)(y)2 orbital is doubly occupied in the gas-phase quintet state.

Camphor 5-Monooxygenase↗
PubMed2004

Properties of fullerene[50] and D5h decachlorofullerene[50]: a computational study.

Stimulated by the recent preparation and characterization of the first [50]fullerene derivative, decachlorofullerene[50] (Science 2004, 304, 699), we have performed a systematic density functional study on the electronic and spectroscopic properties of C(50), its anions and derivatives such as C(50)Cl(10) and C(50)Cl(12). The ground state of C(50) has D(3) symmetry with a spheroid shape, and is highly aromatic; the best D(5h)C(50) singlet is nonaromatic. Both D(3)() and D(5h)() isomers of C(50) have high electron affinities and can be reduced easily. Due to the unstable fused pentagon structural features, C(50) is chemically labile and subject to addition reactions such as chlorination, dimerization and polymerization. The equatorial pentagon-pentagon fusions of D(5h)C(50) are active sites for chemical reactions; hence, D(5h)C(50) may behave as a multivalent group. The computed IR, Raman, (13)C NMR and UV-vis spectra of the D(5h)C(50)Cl(10) molecule agree well with the experimental data. Finally, D(5h)C(50)Cl(10) is predicted to have a high electron affinity and, hence, might serve as an electron-acceptor in photonic/photovoltaic applications. The geometry and (13)C NMR chemical shifts of C(50)Cl(12) were computed to assist further isolation experiments.

Journal Article↗
PubMed2004

Density functional study of valence trapping in a mixed-valent dimanganese complex.

Density functional calculations are reported for the linked system [(Fulv)[Mn(CO)(2)](2)(mu-dppm)] in the neutral and singly and doubly charged states. The neutral system and dication are characterized by well-localized valence states involving two Mn(II) and Mn(I) centers, respectively. At the BP86/AE1 level, the electronic structure of the monocation is that of a delocalized mixed-valence species. At the B3LYP/AE1 level, the spin density of the unpaired electron is largely localized on one Mn atom, in accord with the formulation of a trapped-valence ground state (Atwood, C. G.; Geiger, W. E. J. Am. Chem. Soc. 2000, 122, 5477). The qualitative IR pattern of the CO stretching vibrations is well reproduced at the B3LYP level, supporting the empirical charge distribution parameter employed for the interpretation of the IR spectra.

Computer Simulation↗
PubMed2004

A systematic ab initio study of the equilibrium geometry and vibrational wave numbers of bismuthine.

The equilibrium structure and the harmonic and anharmonic force fields of BiH(3) are determined by high-level ab initio calculations using a variety of correlation treatments, basis sets, and pseudopotentials, partly in combination with core polarization potentials. Spin-orbit effects are included by a configuration interaction treatment. This systematic study serves to establish a reliable computational protocol for such calculations and, in particular, to minimize basis set superposition errors through an improved new basis set and/or counterpoise corrections. Using the recommended procedures, the best ab initio results for the equilibrium geometry and the fundamental vibrational wave numbers are in good agreement with the available experimental data, which further supports the recent spectroscopic identification of BiH(3). The ground-state total atomization energy of BiH(3) is predicted to be 153.1 kcal/mol.

Journal Article↗
PubMed2004

Quantum mechanical/molecular mechanical investigation of the mechanism of C-H hydroxylation of camphor by cytochrome P450cam: theory supports a two-state rebound mechanism.

The stereospecific cytochrome P450-catalyzed hydroxylation of the C(5)-H((5-exo)) bond in camphor has been studied theoretically by a combined quantum mechanical/molecular mechanical (QM/MM) approach. Density functional theory is employed to treat the electronic structure of the active site (40-100 atoms), while the protein and solvent environment (ca. 24,000 atoms) is described by the CHARMM force field. The calculated energy profile of the hydrogen-abstraction oxygen-rebound mechanism indicates that the reaction takes place in two spin states (doublet and quartet), as has been suggested earlier on the basis of calculations on simpler models ("two-state reactivity"). While the reaction on the doublet potential energy surface is nonsynchronous, yet effectively concerted, the quartet pathway is truly stepwise, including formation of a distinct intermediate substrate radical and a hydroxo-iron complex. Comparative calculations in the gas phase demonstrate the effect of the protein environment on the geometry and relative stability of intermediates (in terms of spin states and redox electromers) through steric constraints and electronic polarization.

Camphor↗