PubMed Health⌕ Search

Biomedical subjects

Klaus R Liedl

Publications and source records attributed to Klaus R Liedl.

At least 19 recordsLinked to original sources
PubMed2005

M.TaqI facilitates the base flipping via an unusual DNA backbone conformation.

MD simulations have been carried out to understand the dynamical behavior of the DNA substrate of the Thermus aquaticus DNA methyltransferase (M.TaqI) in the methylation process at N6 of adenine. As starting structures, an x-ray structure of M.TaqI in complex with DNA and cofactor analogue (PDB code: 1G 38) and free decamer d(GTTCGATGTC)(2) were taken. The x-ray structure shows two consecutive BII substates that are not observed in the free decamer. These consecutive BII substates are also observed during our simulation. Additionally, their facing backbones adopt the same conformations. These double facing BII substates are stable during the last 9 ns of the trajectories and result in a stretched DNA structure. On the other hand, protein-DNA contacts on 5' and 3' phosphodiester groups of the partner thymine of flipped adenine have changed. The sugar and phosphate parts of thymine have moved further into the empty space left by the flipping base without the influence of protein. Furthermore, readily high populated BII substates at the GpA step of palindromic tetrad TCGA rather than CpG step are observed in the free decamer. On the contrary, the BI substate at the GpA step is observed on the flipped adenine strand. A restrained MD simulation, reproducing the BI/BII pattern in the complex, demonstrated the influence of the unusual backbone conformation on the dynamical behavior of the target base. This finding along with the increased nearby interstrand phosphate distance is supportive to the N6-methylation mechanism.

Adenine↗
PubMed2005

Z-DNA's conformer substates revealed by FT-IR difference spectroscopy of nonoriented left-handed double helical poly(dG-dC).

Nonoriented hydrated films of double helical poly(dG-dC) in the Z-form were studied by Fourier transform infrared (FT-IR) spectroscopy either as equilibrated slow-cooled samples between 290 and 220 K or, after quenching into the glassy state, as nonequilibrated film isothermally at 200, 220, and 240 K. IR spectral changes on isothermal relaxation at 200 and 220 K toward equilibrium, caused by interconversion of two conformer substates (CS) called Z1 and Z2, are revealed by IR difference spectra. Pronounced spectral changes on Z1-to-Z2 interconversion occur between approximately 750-1250 cm(-1) and these are attributed to structural changes of the phosphodiester-sugar backbone caused by changes of torsion angles, and to decreasing hydrogen-bonding of the ionic phosphate group with water molecules. These spectral changes on Z1-to-Z2 transition can be related to structural differences between ZI and ZII CS observed in single crystals. ZI/ZII CS occurs only at (dGpdC) base steps, and similar behavior is assumed for Z1/Z2. The Z1/Z2 population ratio was determined via curve resolution of marker bands for Z1 and Z2 centered at 785 and 779 cm(-1). This ratio is 0.64 at 290 K, corresponding to 39% of the phosphates of the (dGpdC) base steps in Z1 and 61% in Z2, and it increases to 1.24 on cooling to 220 K. For the Z2<=>Z1 equilibrium, an enthalpy change of -4.9 +/- 0.2 kJ mol(dGpdC)(-1) is obtained from the temperature dependence of the equilibrium constant. Z1 interconverts into Z2 at isothermal relaxation at 200 and 220 K, whereas on slow cooling from ambient temperature, Z2 interconverts into Z1. This unexpected reversal of CS interconversion is attributed to slow restructuring of hydration shells of the CS on quenching, in the same manner reported by Pichler et al. for the BI and BII CS of B-DNA (J. Phys. Chem. B 106, 3263-3274 (2002)). IR difference curves demonstrate two time scales on isothermal relaxation of Z1-->Z2 interconversion, a fast one for structural relaxation of the sugar-phosphate backbone, and a slow one for relaxation of the hydration shells. This slowing down of restructuring of CS hydration shells at approximately 220-240 K could be the cause for the suppression of biological functions at low temperatures.

DNA, Z-Form↗
PubMed2005

Towards an understanding of DNA recognition by the methyl-CpG binding domain 1.

CpG methylation determines a variety of biological functions of DNA. The methylation signal is interpreted by proteins containing a methyl-CpG binding domain (MBDs). Based on the NMR structure of MBD1 complexed with methylated DNA we analysed the recognition mode by means of molecular dynamics simulations. As the protein is monomeric and recognizes a symmetrically methylated CpG step, the recognition mode is an asymmetric one. We find that the two methyl groups do not contribute equally to the binding energy. One methyl group is associated with the major part of the binding energy and the other one nearly does not contribute at all. The contribution of the two cytosine methyl groups to binding energy is calculated to be -3.6 kcal/mol. This implies a contribution of greater than two orders of magnitude to the binding constant. The conserved amino acid Asp32 is known to be essential for DNA binding by MBD1, but so far no direct contact with DNA has been observed. We detected a direct DNA base contact to Asp32. This could be the main reason for the importance of this amino acid. MBD contacts DNA exclusively in the major groove, the minor groove is reserved for histone contacts. We found a deformation of the minor groove shape due to complexation by MBD1, which indicates an information transfer between the major and the minor groove.

5-Methylcytosine↗
PubMed2004

Cooperative effects on the formation of intercalation sites.

Daunomycin is one of the most important agents used in anticancer chemotherapy. It interacts with DNA through intercalation of its planar chromophore between successive base pairs. The effect of intercalation on structure, dynamics and energetics is the topic of a wealth of scientific studies. In the present study, we report a computational examination of the energetics of the intercalation process. In detail, we concentrate on the energetic penalty that intercalation of daunomycin introduces into DNA by disturbing it from its unbound conformation. For these means, we are analyzing already published molecular dynamics simulations of daunomycin-DNA complexes and present novel simulations of a bisdaunomycin-DNA and a 9-dehydroxydaunomycin-DNA intercalated complex using the MM-GBSA module implemented in the AMBER suite of programs. Using this molecular dynamics based, continuum solvent method we were able to calculate the energy required to form an intercalation site. Consequently, we compare the free energy of the duplex d(CGCGCGATCGCGCG)(2) in the B-form conformation with the respective conformations when intercalated with daunomycin and a bisintercalating analog. Our results show that the introduction of one single intercalation site costs approximately 32 kcal/mol. For double intercalation, or intercalation of the bisintercalator, the respective value for one intercalation site decreases to 27 and 24 kcal/mol, respectively, at a theoretical salt concentration of 0.15 M. This proposes that at least in these cases, a synergistic effect takes place. Although it is well known that intercalation leads to substantial disturbance of the DNA conformation, already performed investigations suggest a lower energetic penalty. Nevertheless to the best of our knowledge the calculations presented here are the most complete ones and consider hydration effects for the first time. The interaction energy between the ligand and the DNA certainly over-compensates this penalty for introducing the intercalation site and thus favors complexation. Such analyses are helpful for the description of allosteric effects in protein ligand binding.

Antibiotics, Antineoplastic↗
PubMed2004

Double hydrogen tunneling revisited: the breakdown of experimental tunneling criteria.

Formic acid dimer was chosen as a model system to investigate synchronous double proton transfer by means of variational transition state theory (VTST) for various isotopically modified hydrogen species. The electronic barrier for the double proton transfer was evaluated to be 7.9 kcal/mol, thus being significantly lower than it was determined in previous studies. The tunneling probabilities were evaluated at temperatures from 100 up to 400 K and typical Arrhenius behavior with enhancement by tunneling is observed. When comparing the transmission factors kappa in dependence of the mass of the tunneling hydrogen, it was found that there are two maxima, one at very low masses (e.g., 0.114 amu, corresponding to the muonium entity) and one maximum at around 2 amu (corresponding to deuterium). With the knowledge of the VTST-hydrogen transfer rates and the corresponding tunneling corrections, various tunneling criteria were tested (e.g., Swain-Schaad exponents) and were shown to fail in this reaction in predicting the extent of tunneling. This finding adds another aspect in the ongoing "Tunneling-Enhancement by Enzymes" discussion, as the used tunneling criteria based on experimental reaction rates may fail to predict tunneling behavior correctly.

Journal Article↗
PubMed2004

Water-mediated contacts in the trp-repressor operator complex recognition process.

Water-mediated contacts are known as an important recognition tool in trp-repressor operator systems. One of these contacts involves two conserved base pairs (G(6).C(-6) and A(5). T(-5)) and three amino acids (Lys 72, Ile 79, and Ala 80). To investigate the nature of these contacts, we analyzed the X-ray structure (PDB code: 1TRO) of the trp-repressor operator complex by means of molecular dynamics simulations. This X-ray structure contains two dimers that exhibit structural differences. From these two different starting structures, two 10 ns molecular dynamics simulations have been performed. Both of our simulations show an increase of water molecules in the major groove at one side of the dimer, while the other side remains unchanged compared to the X-ray structure. Though the maximum residence time of the concerned water molecules decreases with an increase of solvent at the interface, these water molecules continue to play an important role in mediating DNA-protein contacts. This is shown by new stable amino acids-DNA distances and a long water residence time compared to free DNA simulation. To maintain stability of the new contacts, the preferential water binding site on O6(G6) is extended. This extension agrees with mutation experiment data on A5 and G6, which shows different relative affinity due to mutation on these bases [A. Joachimiak, T. E. Haran, P. B. Sigler, EMBO Journal 1994, Vol. 13, No. (2) pp. 367-372]. Due to the rearrangements in the system, the phosphate of the base G6 is able to interconvert to the B(II) substate, which is not observed on the other half side of the complex. The decrease of the number of hydrogen bonds between protein and DNA backbone could be the initial step of the dissociation process of the complex, or in other words an intermediate complex conformation of the association process. Thus, we surmise that these features show the importance of water-mediated contacts in the trp-repressor operator recognition process.

Base Pairing↗
PubMed2004

The ground-state tunneling splitting of various carboxylic acid dimers.

Carboxylic acid dimers in gas phase reveal ground-state tunneling splittings due to a double proton transfer between the two subunits. In this study we apply a recently developed accurate semiclassical method to determine the ground-state tunneling splittings of eight different carboxylic acid derivative dimers (formic acid, benzoic acid, carbamic acid, fluoro formic acid, carbonic acid, glyoxylic acid, acrylic acid, and N,N-dimethyl carbamic acid) and their fully deuterated analogs. The calculated splittings range from 5.3e-4 to 0.13 cm(-1) (for the deuterated species from 2.8e-7 to 3.3e-4 cm(-1)), thus indicating a strong substituent dependence of the splitting, which varies by more than two orders of magnitude. One reason for differences in the splittings could be addressed to different barriers heights, which vary from 6.3 to 8.8 kcal/mol, due to different mesomeric stabilization of the various transition states. The calculated splittings were compared to available experimental data and good agreement was found. A correlation could be found between the tunneling splitting and the energy barrier of the double proton transfer, as the splitting increases with increased strength of the hydrogen bonds. From this correlation an empirical formula was derived, which allows the prediction of the ground-state tunneling splitting of carboxylic acid dimers at a very low cost and the tunneling splittings for parahalogen substituted benzoic acid dimers is predicted.

Journal Article↗
PubMed2004

Daunomycin intercalation stabilizes distinct backbone conformations of DNA.

Daunomycin is a widely used antibiotic of the anthracycline family. In the present study we reveal the structural properties and important intercalator-DNA interactions by means of molecular dynamics. As most of the X-ray structures of DNA-daunomycin intercalated complexes are short hexamers or octamers of DNA with two drug molecules per doublehelix we calculated a self complementary 14-mer oligodeoxyribonucleotide duplex d(CGCGCGATCGCGCG)2 in the B-form with two putative intercalation sites at the 5'-CGA-3' step on both strands. Consequently we are able to look at the structure of a 1:1 complex and exclude crystal packing effects normally encountered in most of the X-ray crystallographic studies conducted so far. We performed different 10 to 20 ns long molecular dynamics simulations of the uncomplexed DNA structure, the DNA-daunomycin complex and a 1:2 complex of DNA-daunomycin where the two intercalator molecules are stacked into the two opposing 5'-CGA-3' steps. Thereby--in contrast to X-ray structures--a comparison of a complex of only one with a complex of two intercalators per doublehelix is possible. The chromophore of daunomycin is intercalated between the 5'-CG-3' bases while the daunosamine sugar moiety is placed in the minor groove. We observe a flexibility of the dihedral angle at the glycosidic bond, leading to three different positions of the ammonium group responsible for important contacts in the minor groove. Furthermore a distinct pattern of BI and BII around the intercalation site is induced and stabilized. This indicates a transfer of changes in the DNA geometry caused by intercalation to the DNA backbone.

Binding Sites↗
PubMed2003

C5-methylation of cytosine in B-DNA thermodynamically and kinetically stabilizes BI.

We have examined the backbone dynamics of two alternating purine-pyrimidine dodecamers. One sequence consists of "pure" GC bases; the other one contains 5-methylcytosines. The effect of the methyl groups on the backbone substates BI/BII was investigated by means of molecular dynamics. The methylation influences, on one hand, the transition barrier between BI and BII and, on the other hand, the state of equilibrium. The kinetic consequences are an increase of the DeltaG of Gp5mC steps by 1.5 kcal/mol and a decrease of the DeltaG of 5mCpG steps by 0.8 kcal/mol (compared with the nonmethylated DNA). Thus, the additive group differentiates between the two occurring dinucleotide steps and renders the phosphate of the 5-methylcytosine more rigid, as proposed by experimental studies. The thermodynamic consequences are an increase of the DeltaG of Gp5mC steps by 1.1 kcal/mol and a decrease of the DeltaG of 5mCpG steps by 0.8 kcal/mol. The reason for this shift in equilibrium is still not completely clear on a molecular basis. But we can conclude that the indirect readout of DNA is influenced by methylation.

Computer Simulation↗
PubMed2003

Molecular docking studies of natural cholinesterase-inhibiting steroidal alkaloids from Sarcococca saligna.

Alkaloids isolated from Sarcococca saligna significantly inhibit acetyl- and butyrylcholinesterase enzyme, suggesting discovery of inhibitors for nervous-system disorders. Studying interactions with the active site of the AChE enzyme from Torpedo californica, we have identified hydrophobic interactions inside the aromatic gorge area as the major stabilizing factor in enzyme-inhibitor complexes of these alkaloids. Molecular Dynamics simulation of a predicted complex indicates that ligand binding does not extensively alter enzyme structure, but reduces flexibility at the gorge.

Acetylcholinesterase↗
PubMed2003

A QM-MM interface between CHARMM and TURBOMOLE: implementation and application to systems in bulk phase and biologically active systems.

The implementation of a hybrid QM-MM approach combining ab initio and density functional methods of TURBOMOLE with the molecular mechanics program package CHARMM is described. An interface has been created to allow data exchange between the two applications. With this method the efficient multiprocessor capabilities of TURBOMOLE can be utilized with CHARMM running as a single processor application. Therefore, features of nonparallel running code in CHARMM like the TRAVEL module for locating saddle points or VIBRAN for the calculation of second derivatives can be exploited by running the CPU intensive QM calculations in parallel. To test the methodology, several small systems are studied with both Hartree-Fock and density functional methods and varying QM-MM boundaries. Also, the computationally efficient RI-J method has been examined for use in QM-MM applications. A B(12) cofactor containing cobalt has been studied, to examine systems with a large QM region and transition metals. All tested methods perform satisfactory in comparison with pure quantum calculations. Additionally, algorithms for the characterization of saddle points have been tested for their potential use in QM-MM problems. The TRAVEL module of CHARMM has been applied to the Menshutkin reaction in the condensed phase, and a saddle point was located. This saddle point was verified by calculation of a steepest descent path connecting educt, transition state, and product, and by calculation of vibrational modes.

Cobamides↗
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

Hydration of hydroxypyrrole influences binding of ImHpPyPy-beta-Dp polyamide to DNA.

Ligands which are able to recognize DNA sequence specifically are of fundamental interest as transcription controlling drugs. Recently a polyamide ligand was developed (ImHpPyPy-beta-Dp) which differentiates in a dimeric arrangement between all four possible base pair steps in the minor groove. This is a landmark for the design of DNA binding drugs because it was believed that such a recognition could only be possible in the major groove of DNA. Although the OH groups of the hydroxypyrrole (Hp) moieties of the ligands are responsible for this sequence discrimination, experiments showed that this OH group also reduces the absolute binding constant. We performed a free energy calculation by means of thermodynamic integration in order to find out the influence of this single hydroxyl on DNA binding. In our simulation, we found that the hydroxyl group reduces binding by about 1.3 kcal/mol, which is in excellent agreement with the experimentally determined value of 1.2 kcal/mol. In further MD simulations, the structural reasons for this reduction was estimated. The results of these simulations qualitatively agree with the X-ray structures, but in contrast, in the simulations both (ImHpPyPy-beta-Dp and ImPyPyPy-beta-Dp) ligand-DNA (d(CCAGTACTGG)(2)) complexes exhibit only slight structural differences. This is consistent with a recently published second pair of similar polyamide DNA crystal structures. Thus, we believe that the explanations resulting from the X-ray structures must be modified. We attribute the large structural differences between the two polyamide DNA complexes to a buffer molecule which binds only in the case of the ImHpPyPy-beta-Dp-DNA complex at the region of interest. We propose that the differential hydration of both ligands in the unbound state is responsible for the reduction of the binding constant. Additionally, we suggest an indirect readout of DNA, because of a lengthening of the Watson-Crick base pairs, which possibly contributes to the differentiation between T.A, A.T from G.C, C.G base pairs.

Computer Simulation↗
PubMed2003

Energetic and stereochemical effects of the protein environment on substrate: a theoretical study of methylmalonyl-CoA mutase.

QM/MM methods were used to study the isomerization step from (2R)-methylmalonyl-CoA to succinyl-CoA. A pathway via a "fragmentation-recombination" mechanism is ruled out on energetic grounds. For the other radicalic pathway, involving an addition recombination step, geometries and vibrational contributions have been determined, and a barrier height of 11.70 kcal/mol was found. The effect of adjacent hydrogen-donating groups was found to reduce the energy barrier by 1-2 kcal/mol each and thus to provide a significant catalytic effect for this reaction. By means of molecular dynamics studies, the stereochemistry of the methylmalonyl-CoA mutase catalyzed reaction was examined. It is shown that TYR89 is essential for maintaining stereoselectivity of the abstraction of a hydrogen in the backreaction. The subsequent selective formation of one isomer of methylmalonyl-CoA is probably due to the presence of a bulky side chain.

Acyl Coenzyme A↗
PubMed2003

Stepwise induced fit in the pico- to nanosecond time scale governs the complexation of the even-skipped transcriptional repressor homeodomain to DNA.

Induced fit effects in the complex of a DNA decamer with two even-skipped transcriptional repressor homeodomain molecules were investigated by means of molecular dynamics simulations. Dynamics of these effects are found to be in the time scale from pico- to nanoseconds. First steps are made by the fast-moving DNA backbone phosphates, which upon binding change their B(I)/B(II) substate distribution. Further rearrangements in the DNA double helix induced upon complexation, like bending of the helix axis, changes of the minor groove width, and of different helical parameters, are slower and occur within a few nanoseconds. The flexibility of the DNA, especially of its backbone, seems thereby to play an important role for specific DNA ligand recognition.

Animals↗
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↗