Biomedical subjects

Christofer S Tautermann

Publications and source records attributed to Christofer S Tautermann.

5 recordsLinked to original sources

PubMed2003

Influence of backbone conformations of human carbonic anhydrase II on carbon dioxide hydration: hydration pathways and binding of bicarbonate.

In this study, the hydration of carbon dioxide and the formation of bicarbonate in human carbonic anhydrase II have been examined. From semiempirical QM/MM molecular dynamics studies, dominant conformations of the protein backbone, possibly contributing to the catalytic activity, have been isolated and further examined by means of density functional QM/MM methods. In agreement with experimental observations, a binding site for cyanate, which acts as an inhibitor, has been located, whereas for carbon dioxide, depending on the conformation of the protein environment, either a different binding site or no binding site has been found. In the latter case, carbon dioxide diffuses barrierless to the zinc-bound oxygen, and then a weakly bound bicarbonate complex is formed. The formed complex is characterized by a long C-O bond to the zinc-bound hydroxide. The nature of the calculated stationary points was verified by determination of vibrational frequencies. Finally, the dissociation of the formed bicarbonate from zinc has been considered. Therefore, a water molecule was included in the QM zone of the QM/MM hybrid potential, and minimization yielded a pentacoordinated intermediate. From a potential energy scan, an activation energy of 6.2 kcal/mol for dissociation of bicarbonate from Zn has been found.

Bicarbonates↗

PubMed2003

Modeling anhydrous and aqua copper(II) amino acid complexes: a new molecular mechanics force field parametrization based on quantum chemical studies and experimental crystal data.

This paper presents the vacuum structures of aquacopper(II) bis(amino acid) complexes with glycine, sarcosine, N,N-dimethylglycine, and N-tert-butyl-N-methylglycine estimated using the B3LYP method. The differences between the B3LYP vacuum structures and experimental crystal structures suggested considerable influence of crystal lattice packing effects on the changes in the complexes' geometries. A previously developed molecular mechanics force field for modeling anhydrous copper(II) amino acidates was reoptimized to simulate these changes and predict the properties of both trans and cis anhydrous and aqua copper(II) amino acid complexes. The modeling included experimental molecular and crystal structures of 13 anhydrous and 10 aqua copper(II) amino acidates with the same atom types (Cu(II), C, H, N, and O) but various copper(II) coordination polyhedron geometries, crystal symmetries, and intermolecular interactions. The empirical parameters of the selected potential energy functions were optimized on the B3LYP vacuum copper(II) coordination geometries of three anhydrous copper(II) amino acidates and on experimental crystalline internal coordinates and unit cell dimensions of six anhydrous and six aqua copper(II) amino acid complexes. The respective equilibrium structures were calculated in vacuo and in simulated crystalline environment. The efficacy of the final force field, FFW, was examined. The total root-mean-square deviations between the experimental and theoretical crystal values were 0.018 A in the bond lengths, 2.2 degrees in the valence angles, 5.5 degrees in the torsion angles, and 0.395 A in the unit cell lengths. FFW reproduced the unit cell volumes in the range from -8.1 to 9.6%. The means of Cu to axial water oxygen distances were 2.4 +/- 0.1 A (experiment) and 2.6 +/- 0.1 A (FFW). This paper describes the ability of the molecular mechanics model and FFW force field to simulate the flexibility of the metal coordination polyhedron. The new force field proved effective in predicting the most stable molecular conformation of copper(II) amino acidato systems in vacuo.

Journal Article↗

PubMed2003

Extended method for adiabatic mode reordering.

The task of vibrational mode reordering is very important for reaction valley studies and for the determination of small curvature tunneling effects. An extended algorithm for adiabatic mode reordering is presented. It is based on the method introduced by Konkoli et al. [J Comput Chem 1997, 18, 1282], which is shown to suffer from numerical problems in the region of frequency-crossings and avoided crossings. One improvement is the use of cubic splines for the interpolation of the projected matrix of force constants, which allows larger step sizes between the discrete points along the reaction path, where vibrational analysis is performed. The main improvement is the use of perturbation theory to resolve crossings and avoided crossings. Within this theoretical framework it becomes clear why the method of the maximal overlap between the normal modes cannot work properly, as eigenvectors associated with nearby eigenvalues tend to become "wobbly". Thus a perturbative procedure is designed that is used for all cases where two harmonic frequencies approach each other and the overlap of the associated normal modes is of no practical use. Advantages of the new procedure are the use of larger step sizes along the minimum energy path and the much more reliable resolution of mode-crossings and avoided crossings independent of the systems symmetry. In addition to that it is shown that one should be very cautious in all computational situations when working with eigenvectors associated with nearby eigenvalues.

Journal Article↗

PubMed2002

About the stability of sulfurous acid (H2SO3) and its dimer.

The characterization and isolation of sulfurous acid (H2SO3) have never been accomplished and thus still remain one of the greatest open challenges of inorganic chemistry. It is known that H2SO3 is thermodynamically unstable. In this study, however, we show that a Ci-symmetric dimer of sulfurous acid (H2SO3)2 is 3.5 kcal mol-1 more stable than its dissociation products SO2 and H2O at 77 K. Additionally, we have investigated the kinetic stability of the sulfurous acid monomer with respect to dissociation into SO2 and H2O and the kinetic isotope effect (KIE) on this reaction by transition-state theory. At 77 K, the half-life of H2SO3 is 15 x 10(9) years, but for the deuterated molecule (D2SO3) it increases to 7.9 x 10(26) years. At room temperature, the half-life of sulfurous acid is only 24 hours; however, a KIE of 3.2 x 10(4) increases it to a remarkable 90 years. Water is an efficient catalyst for the dissociation reaction since it reduces the reaction barrier tremendously. With the aid of two water molecules, one can observe a change in the reaction mechanism for sulfurous acid decomposition with increasing temperature. The most likely mechanism below 170 K is via an eight-membered transition-state ring; yet, above 170 K, a mechanism with a six-membered transition state ring becomes the predominant one. For deuterated sulfurous acid, this change in reaction mechanism can be observed at 120 K. Consequently, between 120 and 170 K, different predominant reaction mechanisms occur for the decomposition of normal and deuterated sulfurous acid when assisted by two water molecules. However, the much longer half-life of deuterated sulfurous acid and the stability of the sulfurous acid dimer at 77 K are encouraging for future synthesis and characterization under laboratory conditions.

Journal Article↗

PubMed2002

Towards the experimental decomposition rate of carbonic acid (H2CO3) in aqueous solution.

Dry carbonic acid has recently been shown to be kinetically stable even at room temperature. Addition of water molecules reduces this stability significantly, and the decomposition (H2CO3 + nH2O --> (n+1)H2O + CO2) is extremely accelerated for n = 1, 2, 3. By including two water molecules, a reaction rate that is a factor of 3000 below the experimental one (10 s(-1)) at room temperature was found. In order to further remove the gap between experiment and theory, we increased the number of water molecules involved to 3 and took into consideration different mechanisms for thorough elucidation of the reaction. A mechanism whereby the reaction proceedes via a six-membered transition state turns out to be the most efficient one over the whole examined temperature range. The determined reaction rates approach experimental values in aqueous solution reasonably well; most especially, a significant increase in the rates in comparison to the decomposition reaction with fewer water molecules is found. Further agreement with experiment is found in the kinetic isotope effects (KIE) for the deuterated species. For water-free carbonic acid, the KIE (i.e., kH2CO3/kD2CO3) for the decomposition reaction is predicted to be 220 at 300 K, whereas it amounts to 2.2-3.0 for the investigated mechanisms including three water molecules. This result is therefore reasonably close to the experimental value of 2 (at 300 K). These KIEs are in much better accordance with the experiment than the KIE for decomposition with fewer water entities.

Carbonic Acid↗