Biomedical subjects

Bernd Wellenzohn

Publications and source records attributed to Bernd Wellenzohn.

17 recordsLinked to original sources

PubMed2006

Identification of atherosclerosis-associated conformational heat shock protein 60 epitopes by phage display and structural alignment.

Autoimmune reactions to HSP60 are believed to play a key role during development of early atherosclerosis. Due to the high degree of phylogenetic conservation between microbial and human HSP60, bacterial infections might be responsible for inducing cross-reactivity to self HSP60, which is expressed on the surface of arterial endothelial cells stressed by classical atherosclerosis risk factors. Conformational epitopes recognized by polyclonal anti-mycobacterial HSP60 antibodies from subjects with atherosclerosis were identified using a phage displayed random library of cyclic constrained 7mer peptides. After five rounds of selection, DNA sequencing of strongly binding clones revealed that one peptide motif (CIGSPSTNC) was present in 64% of all clones, and a second motif (CSFHYQNRC) in 14%. Using a newly developed method for structural alignment of small constrained peptides onto a protein surface, we located the motif present in 14% of all clones on the surface of mycobacterial HSP60. The motif present in 64% of all clones was found on the surface of mycobacterial HSP60 as well as in the homologous region of human HSP60, which makes this epitope a promising candidate for further investigations on cross-reactive epitopes involved in early atherogenesis.

Amino Acid Sequence↗

PubMed2005

The N6-methyl group of adenine further increases the BI stability of DNA compared to C5-methyl groups.

Methylated DNA bases are natural modifications which play an important role in protein-DNA interactions. Recent experimental and theoretical results have shown an influence of the base modification on the conformational behavior of the DNA backbone. MD simulations of four different B-DNA dodecamers (d(GC)(6), d(AT)(6), d(G(5mCG)(5)C), and d(A(T6mA)(5)T)) have been performed with the aim to examine the influence of methyl groups on the B-DNA backbone behavior. An additional control simulation of d(AU)(6) has also been performed to examine the further influence of the C5-methyl group in thymine. Methyl groups in the major groove (as in C5-methylcytosine, thymine, or N6-methyladenine) decrease the BII substate population of RpY steps. Due to methylation a clearer distinction of the BI substate stability between YpR and RpY (CpG/GpC or TpA/ApT) steps arises. A positive correlation between the BII substate population and base stacking distances is seen only for poly(GC). A methyl group added into the major groove increases mean water residence times around the purine N7 atom, which may stabilize the BI substate by improving the hydration network between the DNA backbone and the major groove. The N6-methyl group also forms a water molecule bridge between the N6 and O4 atoms, and thus further stabilizes the BI substate.

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

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

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

3D-QSAR studies on natural acetylcholinesterase inhibitors of Sarcococca saligna by comparative molecular field analysis (CoMFA).

We have derived a comprehensive structure-activity relationship (SAR) picture for a new series of natural acetylcholinesterase inhibitors isolated from Sarcococca saligna. A set of 32 previously isolated and tested pregnane-type steroidal alkaloids inhibitors were investigated with respect to their IC(50) values (pIC(50)) against the AChE enzyme in order to derive CoMFA models using atom-based alignment. A highly significant CoMFA model was obtained with r(2) value of 0.974. The q(2) (cross validation r(2)) value also confirms the statistical significance of our model.

Acetylcholinesterase↗

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

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

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

The conformer substates of nonoriented B-type DNA in double helical poly(dG-dC).

A nonoriented hydrated film of poly(dG-dC) with ?20 water molecules per nucleotide (called B* by Loprete and Hartman (Biochem. 32, 4077-4082 (1993)) was studied by Fourier transform infrared (FT-IR) spectroscopy either as equilibrated sample between 290 and 270 K or, after quenching into the glassy state, as nonequilibrated film isothermally at 200 and 220 K. IR spectral changes on isothermal relaxation at 200 and 220 K, caused by interconversion of two conformer substates, are revealed by difference spectra. Comparison with difference curves obtained in the same manner from two classical B-DNA forms, namely the d(CGCGAATTCGCG)(2) dodecamer and polymeric NaDNA from salmon testes, revealed that the spectral changes on B(I)-to-B(II) interconversion in the classical B-DNA forms are very similar to those in the B*-form, and that the spectroscopic differences between the B(I) and B(II) features from classical B-DNA and those from the modified B*-form are minor. Nonexponential kinetics of the B(I)-->B(II) transition in the B*-form of poly(dG-dC) at 200 K showed that the structural relaxation time is about three times of that in the classical B-DNA forms (approximately equal to 30 versus approximately equal to 10 min at 200 K). The unexpected reversal of conformer substates interconversion (that is B(II)-->B(I) transition on cooling from 290 K and B(I)-->B(II) transition on isothermal relaxation at 200 K) observed for classical B-DNA occurs also in the modified B*-form. We therefore conclude that restructuring of hydration shells rules the low-temperature dynamics of the B*-form via its two conformer substates in the same manner reported for classical B-DNA by Pichler et al. (J. Phys. Chem. B 106, 3263-3274 (2002)).

Animals↗

PubMed2002

Characterization of the vitamin E-binding properties of human plasma afamin.

Human plasma afamin, the fourth member of the albumin gene family, is shown to be a specific binding protein for vitamin E. A radio ligand-binding assay followed by Scatchard and Hill analysis are used to show that afamin has a binding affinity for both alpha-tocopherol and gamma-tocopherol, two of the most important forms of vitamin E, in vitro. The binding-dissociation constant was determined to be 18 microM, indicating that afamin plays a role as vitamin E carrier in body fluids such as human plasma and follicular fluid under physiological conditions. Additionally, we demonstrate that afamin has multiple binding sites for both alpha- and gamma-tocopherol. Due to the large binding capacity of afamin for vitamin E, it might take over the role of vitamin E transport in body fluids under conditions where the lipoprotein system is not sufficient for vitamin E transport. To confirm the experimental results, we performed homology modeling and docking calculations on the predicted tertiary structure, which showed coincidence between calculated and in vitro results.

Binding Sites↗

PubMed2002

PvuII-endonuclease induces structural alterations at the scissile phosphate group of its cognate DNA.

We investigated the PvuII endonuclease with its cognate DNA by means of molecular dynamics simulations. Comparing the complexed DNA with a reference simulation of free DNA, we saw structural changes at the scissile phosphodiester bond. At this GpC step, the enzyme induces the highest twist and axial rise, inclination is increased and the minor groove widened. The distance between the scissile phosphate group and the phosphate group of the following thymine base is shortened significantly, indicating a substrate-assisted catalysis. A feasible reason for this vicinity is the catalytically important amino acid residue lysine 70, which bridges the free oxygen atoms of the successive phosphate groups. Due to this geometry, a compact reaction pocket is formed where a water molecule can be held, thus bringing the reaction partners for hydrolysis into contact. The O1-P-O2 angle of the scissile nucleotide is decreased, probably due to a complexation of the negative oxygen atoms through protein and solvent contacts.

Binding Sites↗

PubMed2002

Indirect readout of the trp-repressor-operator complex by B-DNA's backbone conformation transitions.

Although the trp-repressor-operator complex is one of the best studied transcriptional controlling systems, some questions regarding the specific recognition of the operator by the repressor remain. We performed a 2.35 ns long molecular dynamics simulation to clarify the influence of the two B-DNA backbone conformational substates B(I) and B(II) on complexation. The trp-repressor-operator is an ideal biological system for this study because experimental results have already figured out that the interaction between the internucleotide phosphates and the protein is essential for the formation of the high affinity complex. Our simulation supports these results, but more important it shows a strong correlation between the B(I)/B(II) phosphate substate and the number of interactions with this phosphate. In particular the B(I) <==> B(II) transitions occur synchronous to hydrogen bond breaking or formation. To the best of our knowledge, this was observed for the first time. Thus, we conclude that the sequence specific B(I)/B(II) behavior contributes via indirect readout to sequence specific recognition. These results have implication for the design of transcription-controlling drugs in view of the recently published influence of minor groove binders on the B(I)/B(II) pattern. The simulation also agrees with crystallographically observed hydration sites. This is consistent with experimental results and indicates the correctness of the model used.

Bacterial Proteins↗

PubMed2002

Comparative sequencing of nucleic acids by liquid chromatography-tandem mass spectrometry.

An algorithm was developed for the computer-aided interpretation of fragment ion spectra from collision-induced dissociation of multiply charged oligodeoxynucleotide ions generated by electrospray ionization. The method compares the experimental spectrum to the m/z values predicted by employing established fragmentation pathways from a known reference sequence. The closeness of matching between the measured spectrum and the predicted set of fragment ions is characterized by the fitness, which takes into account the difference between measured and predicted m/z values, the intensity of the fragment ions, the number of fragments assigned, and the number of nucleotide positions not covered by fragment ions in the experimental spectrum. Smaller values for the fitness indicate a closer match between measured spectrum and predicted m/z values. To substantiate the identity of investigated sequence and reference sequence, or to identify point mutations or insertions/deletions, the reference sequence is systematically varied by incorporating all four possible nucleotides A, T, G, and C at each position in the sequence followed by identification of the correct sequence by the lowest fitness value. Collision energy was shown to have a major impact on the interpretability of the tandem mass spectra by the comparative sequencing algorithm, and the optimal collision energy depended on the length of the fragmented oligodeoxynucleotide. The analytical system was successfully applied to verify DNA sequences as well as to detect and localize point mutations or insertions/deletions in 5-51-mer oligodeoxynucleotides.

Algorithms↗