PubMed Health⌕ Search

Biomedical subjects

Bernd Ensing

Publications and source records attributed to Bernd Ensing.

13 recordsLinked to original sources
PubMed2007

Proton shuttles and phosphatase activity in soluble epoxide hydrolase.

Recently, a novel metal Mg2+-dependent phosphatase activity has been discovered in the N-terminal domain of the soluble epoxide hydrolase (sEH), opening a new branch of fatty acid metabolism and providing an additional site for drug targeting. Importantly, the sEH N-terminal fold belongs to the haloacid dehalogenase (HAD) superfamily, which comprises a vast majority of phosphotransferases. Herein, we present the results of a computational study of the sEH phosphatase activity, which includes classical molecular dynamics (MD) simulations and mixed quantum mechanical/molecular mechanics (QM/MM) calculations. On the basis of experimental results, a two-step mechanism has been proposed and herein investigated: (1) phosphoenzyme intermediate formation and (2) phosphoenzyme intermediate hydrolysis. Building on our earlier work, we now provide a detailed description of the reaction mechanism for the whole catalytic cycle along with its free energy profile. The present computations suggest metaphosphate-like transition states for these phosphoryl transfers. They also reveal that the enzyme promotes water deprotonation and facilitates shuttling of protons via a metal-ligand connecting water bridge (WB). These WB-mediated proton shuttles are crucial for the activation of the solvent nucleophile and for the stabilization of the leaving group. Moreover, due to the conservation of structural features in the N-terminal catalytic site of sEH and other members of the HAD superfamily, we suggest a generalization of our findings to these other metal-dependent phosphatases.

Catalysis↗
PubMed2006

Self-assembling cyclic peptides: molecular dynamics studies of dimers in polar and nonpolar solvents.

The self-assembly of cyclic D,L-alpha-peptides into hollow nanotubes is a crucial mechanistic step in their application as antibacterial and drug-delivery agents. To understand this process, molecular dynamics (MD) simulations were performed on dimers of cyclic peptides formed from cyclo [(-L-Trp-D-N-MeLeu-)4-]2 and cyclo [(-L-Trp-D-Leu-)4-]2 subunits in nonpolar (nonane) and polar (water) solvent. The dimers were observed to be stable only in nonpolar solvent over the full 10 ns length of the MD trajectory. The behavior of the dimers in different solvents is rationalized in terms of the intersubunit hydrogen bonding, hydrogen bonding with the solvent, and planarity of the rings. It is shown that the phi and psi dihedral angles of a single uncapped ring in nonane lie in the beta-sheet region of the Ramachandran plot, and the ring stays in a flat conformation. Steered MD (SMD) simulations based on Jarzynski's equality were performed to obtain the potential of mean force as a function of the distance between the two rings of the capped dimer in nonane. It is also shown that a single peptide subunit prefers to reside close to the nonane/water interface rather than in bulk solvent because of the amphiphilic character of the peptide ring. The present MD results build the foundation for using MD simulations to study the mechanism of the formation of cyclic peptide nanotubes in lipid bilayers.

Dimerization↗
PubMed2006

Metadynamics as a tool for exploring free energy landscapes of chemical reactions.

The metadynamics or hills method is a relatively new molecular dynamics technique aimed to enhance the sampling of separated regions in phase space and map out the underlying free energy landscape as a function of a small number of order parameters or collective variables. The high efficiency allows for the application of metadynamics in combination with first principles dynamics methods, in particular with Car-Parrinello molecular dynamics, to study processes in which changes in the electronic structure play a dominant role, such as chemical reactions. The option to choose several independent collective variables is important to tackle complex and concerted transformations that lack an obvious a priori choice for a single reaction coordinate. In this Account, we discuss the role of metadynamics in the search of transition states, local minima, reaction paths, free energy profiles, and reaction coordinates among a growing list of alternative methods.

Journal Article↗
PubMed2005

Computational study of phosphatase activity in soluble epoxide hydrolase: high efficiency through a water bridge mediated proton shuttle.

Recently, a new branch of fatty acid metabolism has been opened by the novel phosphatase activity found in the N-terminal domain of the, hence bifunctional, soluble epoxide hydrolase (sEH). Importantly, this finding has also provided a new site for drug targeting in sEH's activity regulation. Classical MD and hybrid Car-Parrinello QM/MM calculations have been performed to investigate the reaction mechanism of the phosphoenzyme intermediate formation in the first step of the catalysis. The results support a concerted multi-event reaction mechanism: (1) a dissociative in-line nucleophilic substitution for the phosphoryl transfer reaction; (2) a double proton transfer involved in the formation of a good leaving group in the transition state. The presence of a water bridge in the substrate/enzyme complex allowed an efficient proton shuttle, showing its key role in speeding up the catalysis. The calculated free energy of the favored catalytic pathway is approximately 19 kcal/mol, in excellent agreement with experimental data.

Animals↗
PubMed2005

Perspective on the reactions between F- and CH3CH2F: the free energy landscape of the E2 and SN2 reaction channels.

Recently, we computed the 3D free energy surface of the base-induced elimination reaction between F(-) and CH(3)CH(2)F by using a powerful technique within Car-Parrinello molecular dynamics simulation. Here, the set of three order parameters is expanded to six, which allows the study of the competing elimination and substitution reactions simultaneously. The power of the method is exemplified by the exploration of the six-dimensional free energy landscape, sampling, and mapping out the eight stable states as well as the connecting bottlenecks. The free energy profile and barrier along the E2 and S(N)2 reaction channels are refined by using umbrella sampling. The two mechanisms do not share a common "E2C-like" transition state. Comparison with the zero temperature profiles shows a particularly significant entropy contribution to the S(N)2 channel.

Algorithms↗
PubMed2005

A recipe for the computation of the free energy barrier and the lowest free energy path of concerted reactions.

The recently introduced hills method (Proc. Natl. Acad. Sci. U.S.A. 2002, 99, 12562) is a powerful tool to compute the multidimensional free energy surface of intrinsically concerted reactions. We have extended this method by focusing our attention on localizing the lowest free energy path that connects the stable reactant and product states. This path represents the most probable reaction mechanism, similar to the zero temperature intrinsic reaction coordinate, but also includes finite temperature effects. The transformation of the multidimensional problem to a one-dimensional reaction coordinate allows for accurate convergence of the free energy profile along the lowest free energy path using standard free energy methods. Here we apply the hills method, our lowest free energy path search algorithm, and umbrella sampling to the prototype S(N)2 reaction. The hills method replaces the in many cases difficult problem of finding a good reaction coordinate with choosing relatively simple collective variables, such as the bond lengths of the broken and formed chemical bonds. The second part of the paper presents a guide to using the hills method, in which we test and fine-tune the method for optimal accuracy and efficiency using the umbrella sampling results as a reference.

Journal Article↗
PubMed2005

Lipid bilayer perturbations around a transmembrane nanotube: a coarse grain molecular dynamics study.

The perturbations induced in a lipid bilayer by the presence of a transmembrane nanotube are investigated using coarse grained molecular dynamics. Meniscus formation by the lipids and tilting of the nanotube occur in response to hydrophobic mismatch, although these two effects do not compensate completely for the total mismatch. The lipid head-to-tail vector field is examined and shows strong ordering in the membrane plane regardless of the nanotube length. Molecular layering at the lipid-nanotube interface is reported. This study extends previous theoretical approaches to a more realistic setting.

Biophysical Phenomena↗
PubMed2005

Structure and dynamics of model pore insertion into a membrane.

A cylindrical transmembrane molecule is constructed by linking hydrophobic sites selected from a coarse grain model. The resulting hollow tube assembly serves as a representation of a transmembrane channel, pore, or a carbon nanotube. The interactions of a coarse grain di-myristoyl-phosphatidyl-choline hydrated bilayer with both a purely hydrophobic tube and a tube with hydrophilic caps are studied. The hydrophobic tube rotates in the membrane and becomes blocked by lipid tails after a few tens of nanoseconds. The hydrophilic sites of the capped tube stabilize it by anchoring the tube in the lipid headgroup/water interfacial region of each membrane leaflet. The capped tube remains free of lipid tails. The capped tube spontaneously conducts coarse grain water sites; the free-energy profile of this process is calculated using three different methods and is compared to the barrier for water permeation through the lipid bilayer. Spontaneous tube insertion into an undisturbed lipid bilayer is also studied, which we reported briefly in a previous publication. The hydrophobic tube submerges into the membrane core in a carpetlike manner. The capped tube laterally fuses with the closest leaflet, and then, after plunging into the membrane interior, rotates to assume a transbilayer orientation. Two lipids become trapped at the end of the tube as it penetrates the membrane. The hydrophilic headgroups of these lipids associate with the lower tube cap and assist the tube in crossing the interior of the membrane. When the rotation is complete these lipids detach from the tube caps and fuse with the lower leaflet lipids.

Binding Sites↗
PubMed2004

A minimum free energy reaction path for the E2 reaction between fluoro ethane and a fluoride ion.

The prototype binuclear elimination (E2) reaction illustrates the mechanism of a large number of biochemical and industrial applied processes but has received surprisingly little attention in theoretical studies compared to, for example, the substitution (SN2) reaction. This is due to its concerted mechanism, which requires an independent description of the three bonds that are being formed or broken. In this work, we have taken the advantage of a new and promising methodology to efficiently sample intrinsically multidimensional free-energy surfaces. We locate the lowest free-energy reaction path in the 3D configurational space and use this finite-temperature intrinsic reaction coordinate in an umbrella sampling scheme to access the temperature contributions to high accuracy. The small increase of the barrier and the decrease of the overall endothermicity for the E2 reaction due to entropic contributions is non-trivial. Moreover, our strategy to efficiently handle multiple reaction coordinates could be a great benefit to many chemistry-related fields, such as enzyme catalysis, reactions in solution, and nucleation processes.

Journal Article↗
PubMed2004

Methane-to-methanol oxidation by the hydrated iron(IV) oxo species in aqueous solution: a combined DFT and car-parrinello molecular dynamics study.

Previously, we have shown that the ferryl ion ([FeIVO]2+) is easily produced from Fenton's reagent (i.e., a mixture of Fe2+ ions and H2O2 in aqueous solution), using DFT and Car-Parrinello MD calculations. To verify that the ferryl ion can indeed act as the active species in oxidation reactions with Fenton's reagent, we study in the present paper the reactivity of the ferryl ion toward an organic substrate, in particular the oxidation of methane to methanol. In the first part of this paper, we perform static DFT calculations on the reaction of CH4 with the [(H2O)5FeIVO]2+ complex in vacuo that show a strong prevalence of the oxygen-rebound mechanism over the methane coordination mechanism. This is in agreement with the static DFT results for methane oxidation by biocatalysts MMO and P450, but not with those for methane oxidation by bare metal-oxo ions, where the methane coordination mechanism prevails. The highest energy barrier in the oxygen-rebound mechanism is only 3 kcal/mol in vacuo, whereas in the methane coordination mechanism the highest barrier is 23 kcal/mol. Overall the oxidation reaction energy is downhill by 47 kcal/mol. We conclude that the ferryl ion can indeed act as the oxidative intermediate in the Fenton oxidation of organic species. In the second part of this paper, we perform a preliminary assessment of solvent effects on the oxidation by the ferryl ion in aqueous solution using the method of constrained (first principles) molecular dynamics. The free energy barrier of the H-abstraction reaction from methane by the ferryl ion (i.e., the first step in the rebound mechanism) in aqueous solution is, with 22 kcal/mol in solution, significantly higher than in vacuo. Given the fact that methane has a relatively strong C-H bond (ca. 10 kcal/mol stronger than the C-H bonds in the more typical Fenton's reagent substrates), we infer that for many organic substrates oxidation with the ferryl ion as an active intermediate may be a perfectly viable route.

Journal Article↗
PubMed2003

O2 evolution in the Fenton reaction.

We investigated the mechanism involved in the oxygen production in the Fenton chemistry by means of density functional theory calculations. This study extends previous work in which we proposed that the Fe(IV)O2+ complex is the key active intermediate in the Fenton reaction. Here we provide a consistent picture of the entire reaction cycle by analyzing how the active species, Fe(IV)O2+, can react with hydrogen peroxide to produce O2 and regenerate the Fe2+ catalyst. These results are also relevant in view of the analogies with important enzyme-catalyzed oxidation reactions.

Journal Article↗