Refine
Has Fulltext
- yes (366) (remove)
Is part of the Bibliography
- yes (366)
Year of publication
Document Type
- Journal article (237)
- Doctoral Thesis (113)
- Preprint (14)
- Book article / Book chapter (1)
- Report (1)
Language
- English (366) (remove)
Keywords
- Organische Chemie (68)
- Supramolekulare Chemie (21)
- Selbstorganisation (18)
- self-assembly (16)
- Farbstoff (13)
- perylene bisimide (13)
- fluorescence (12)
- water oxidation (12)
- Merocyanine (11)
- RNA (11)
Institute
- Institut für Organische Chemie (366) (remove)
Schriftenreihe
Sonstige beteiligte Institutionen
- International Max Planck Research School Molecular Biology, University of Göttingen, Germany (2)
- Agricultural Center, BASF SE, 67117 Limburgerhof, Germany (1)
- Center for Nanoscale Microscopy and Molecular Physiology of the Brain (CNMPB), Göttingen, Germany (1)
- Center for Nanosystems Chemistry (1)
- Center for Nanosystems Chemistry (CNC), University of Würzburg (1)
- Center for Nanosystems Chemistry (CNC), Universität Würzburg, Am Hubland, 97074 Würzburg, Germany (1)
- Charles University, Faculty of Mathematics and Physics, Ke Karlovu 5, 121 16 Prague, Czech Republic (1)
- Cluster of Excellence "Multiscale Bioimaging: from Molecular Machines to Networks of Excitable Cells, Göttingen (1)
- Department of Cellular Biochemistry, University Medical Center Göttingen (1)
- Department of Cellular Biochemistry, University Medical Centre Göttingen (1)
In this work the influence of “active” bridge units on the electron transfer (ET) mechanism within organic donor-bridge-electrode arrays in self-assembled monolayers (SAMs) was studied by spectroscopic and electrochemical methods. In the first part of this work ferrocenealkanethiols 1 – 3 and the ferrocenearylthiols 4, 5 were investigated to get experience in the monolayer preparation for measuring ET rates. Cyclic voltammetry of the monolayers indicates that homogeneously mixed monolayers containing redox active molecules and dummy molecules were formed. For the known ferrocenealkanethiols 1 – 3 the ET rates could be confirmed compared to the ones measured by Creager et al. [206]. As expected the ET rate decreases by increasing chain length of the alkane spacer from 2 to 3. Changing the bonding between the redox centre and the alkane spacer with the same bridge lenght, e. g. by using a carboxy-group in case of 1, does not influence the ET behaviour very strong. The aromatic ferrocenethiols 4 and 5 show very high ET rates due to the strong conjugated system although the distance between the redox centre and the electrode is comparable to the C8-alkyl compound 2. The electronic coupling factors all indicate a nonadiabatic ET between the redox centre and the electrode. As expected the electronic coupling factors increase with decreasing spacer length or with an enlarged conjugated system. To sum up, experience in monolayer preparation could be obtained, the measured ET rates for well known ferrocenealkane-compounds 1 - 3 could be verified and the information could be transferred to the conjugated systems 4 and 5. In the second part the triarylamine- 29, 32 and the phenothiazinealkanethiol 35 have been examined relative to their ET behaviour in mixed monolayers. The cyclic voltammograms of the diluted monolayers indicate that homogeneously formed monolayers are present. The ET rates of triarylamine- 29, 32 and phenothiazinealkanethiols 35 are 10 to 100 times higher than compared to ferrocenealkanethiols with equal chain length[183, 206], whereas in a [Ru(bpy)2(pp)]+-containing monolayer the same value was observed [177]. Almost two parameters influence the ET rate constant: the electronic coupling matrix element and the reorganisation energy  [209]. The ET rate in donor substituted alkanethiols is mainly influenced by the reorganisation energy  [177] and even small changes have a dramatic effect on the observed processes, therefore an increasing ET rate from the ferrocene (high reorganisation energy) over the phenothiazine 35 and the [Ru(bpy)2(pp)]+ to the triarylamine chromophores 29 and 32 (low reorganisation energy) is observed. Furthermore the bonding between the redox centres and the alkane spacer plays an important role on the ET rate in case of the triarylamines 29 and 32 opposite to the assumption made by Creager et al. that the connection does not play any role. For the electron rich ether connected compound 29 the ET is not only dominated by the reorganisation energy but also by mesomeric effects where the positive charge of the electron rich derivative 29 is more located at the ether function so that the chain is formally shortend by one atom resulting in higher ET rates than compared to 32. In the third part of the thesis a series of “molecular wires” consisting of methoxy- or chloro-substituted triarylamines and phenothiazines with different bridge units and bridge length between the redox centre and the anchor thiol function have been prepared in order to investigate their ET-behaviour. Cyclic voltammetry and UV/vis-spectroscopy show that the oxidation potential and the energetic states could be controlled very well by introducing different redox centres and bridge units resulting in a decreasing oxidation potential of the redox centres and a bathochromic shift of the absorption bands in the UV/vis-spectra. Also the densitiy of the chromophores in mixed monolayers could be controlled very well for only three compounds (49, 52 and 87) with nitrile-substituted bridges reliable ET rates could be obtained. In these chromophores the ET rate decreases by increasing the density of the redox active molecules in the mixed monolayers indicating that the adsorption geometry changes with coverage with the chromophores tilting to a more upright orientation as the surface becomes more crowded. For all other compounds the measurements were limited by the fast ET rates. Conformational, as well as a very weak distance dependence of the ET resulting in very high ET rates [172] or unfavourable HOMO-LUMO energies of the donor, bridge and the electrode are reasons for this behaviour. The fact that compound 49 shows almost the same rate constant independent of the length (n = 2 or n = 3) may indicate that a hopping process is operating for which a much weaker length dependence is expected than in the case of a superexchange.
Complexation properties of 2,2':6',2''-terpyridine (tpy) have been studied with a series of first row transition metal ions by UV-vis, 1H NMR and isothermal titration calorimetry and ƒ´H values for the tpy complexation processes have been determined. These studies reveal that Zn2+ is the best suited metal ion for the reversible coordination of the terpyridine ligand. Thus, supramolecular coordination polymerization of perylene bisimide fluorophores containing terpyridine functionalities have been investigated by using Zn2+ as metal ion. The formation of the dimeric complexes in the case of monotopic model comounds and coordination polymerization of ditopic functional building blocks have been confirmed by 1H NMR studies. The optical properties of dimeric and polymeric complexes have been investigated by UV-vis and fluorescence spectroscopy. The Zn2+ coordination to the terpyridine unit does not effect the advantageous fluorescence properties of perylene bisimide moieties. The reversibility of the formation of coordination polymers has been established by 1H NMR and additionally by DOSY NMR and fluorescence anisotropy measurements. Coordination polymer strands can be visualized by atomic force microscopy (AFM), which also reveals the formation of an ordered monolayer film at higher concentration. The average polymer length has been determined by AFM to 15 repeat units, which correlates well with the value estimated by 1H NMR to >10 repeat units.
The present thesis encompasses two parts. The first supramolecular part focuses on the development of new flexible self-assembling zwitterions as building blocks for supramolecular polymers. In the second part, the aim was to develop bioorganic receptors for amino acids and dipeptides in aqueous media. Both research projects are based on the guanidiniocarbonyl pyrrole 1 as a new efficient binding motif for the complexation of carboxylates in polar solution.A necessary requirement for the realization of these research projects was to develop an efficient and mild synthetic approach for the cationic guanidiniocarbonyl pyrroles in general. The harsh reaction conditions of the previously used method and the problematic purification of the cationic guanidinocarbonyl pyrroles so far prevented a more extensive exploration in bioorganic and supramolecular research. In the course of this work I successfully developed a new synthesis starting with mono tBoc-protected guanidine that was coupled with a benzyl protected pyrrole carboxylic acid. After deprotection of the benzyl group, a key intermediate in the newly developed synthesis, the tBoc-protected guanidinocarbonyl pyrrole acid, was obtained. This new, mild and extremely efficient synthetic approach for the introduction of acyl guanidines is now the standard procedure in our group for the preparation of both solution and solid-phase guanidiniocarbonyl pyrroles. With this facile method at hand, a new class of flexible zwitterions, in which a carboxylate is linked via an alkyl chain to a guanidiniocarbonyl pyrrole cation was synthesized. The self-aggregation and the influence of the length and therefore flexibility of the alkyl spacer on the structure and stability of the formed aggregates were studied in solution and gas phase. In solution the aggregation was studied by NMR-dilution experiments in DMSO which suggest that flexible zwitterions with n = 1, 3 and 5 form oligomers. For n = 1, highly stable helical aggregates with nanometer size are formed. In the gas phase studies the stability and the fragmentation kinetics of a series of sodiated dimeric zwitterions with n = 2, 3 and 5 were investigated. This was done by infrared multiphoton dissociation Fourier transform ion cyclotron resonance mass spectrometry (IRMPD-FT-ICR-MS). These kinds of studies can be used in the future for a more directed design of supramolecular building blocks The bioorganic research part comprises three different projects. In a first project I synthesized four new arginine analogues which can be implemented in peptides as a substitute for arginine. Therefore, I developed the new multi-step synthesis shown below for these arginine analogues. As a test for their application in normal solid phase synthesis, I successfully prepared a tripeptide sequence Ala-AA1-Val (AA: arginine analogue. In a second project I studied the influence of additional ionic interactions within our binding motif. I synthesized a di-cationic and a tris-cationic receptor and evaluated the binding properties via NMR titration experiments against a variety of amino acids. Especially, the tris-cationic receptor was capable to strongly complex amino acids. The association constants were about a factor of 100 higher than those for the guanidiniocarbonyl pyrroles known so far. Even in 90 %water/10 % DMSO the association constants determined by NMR titration were extremely high with values around Kass = 2000 M-1. In the third project I developed a de-novo designed receptor for C-terminal dipeptides in a beta-sheet conformation based on molecular calculations. This receptor was studied in NMR and also UV titration experiments. In 40 % water/ 60 % DMSO the association constants were too strong to be measured by NMR titration experiments. Therefore, the complexation properties of 12 were studied by UV titration in water (with 10 % DMSO added for solubility reasons) with various dipeptides and amino acids as substrates. The data show that 12 binds dipeptides very efficiently even in water with association constants Kass > 10000 M-1, making 12 one of the most effective dipeptide receptors known so far. In contrast to that, simple amino acids are bound up to ten times less efficiently (Kass > 1000 M-1) than dipeptides. In the series of dipeptides studied the complex stability increases depending on the side chains present in the order Gly < Ala < Val which is a result of the decreasing flexibility of the peptide and the increasing hydrophobicity of the side chains. The binding properties of this receptor are superior to any other dipeptide receptor reported so far. Within my thesis I have not only developed an essential, mild and efficient synthetic approach for guanidiniocarbonyl pyrroles in general, but also a new binding motif for the complexation of amino acids 15, 11 and in addition a dipeptide receptor 12 that is superior to all dipeptides receptors known so far.
Summary The nature of the chemical bond is a topic under constant debate. What is known about individual molecular properties and functional groups is often taught and rationalized by explaining Lewis structures, which, in turn, make extensive use of the valence concept. The valence concept distinguishes between electrons, which do not participate in chemical interactions (core electrons) and those, which do (single, double, triple bonds, lone-pair electrons, etc.). Additionally, individual electrons are assigned to atomic centers. The valence concept is of paramount success: It allows the successful planning of chemical syntheses and analyses, it explains the behavior of individual functional groups, and, moreover, it provides the “language” to think of and talk about molecular structure and chemical interactions. The resounding success of the valence concept may be misleading to forget its approximative character. On the other hand, quantum mechanics provide in principle a quantitative description of all chemical phenomena, but there is no discrimination between electrons in quantum mechanics. From the quantum mechanical point of view there are only indistinguishable electrons in the field of the nuclei, i.e., it is impossible to assign a given electron to a particular center or to ascribe a particular purpose to individual electrons. The concept of indistinguishability of micro particles is founded on the Heisenberg uncertainty relation, which states, that wavepackets diverge in the 6N dimensional phase space, such that individual trajectories can not be identified. Hence it is a deep-rooted and approved physical concept. As an introduction to the present work density partitioning schemes were discussed, which divide the total molecular density into chemically meaningful areas. These partitioning schemes are intimately related to either the concepts of bound atoms in a molecule (as in the Atoms In Molecules theory (AIM) according to Bader or as in the Hirshfeld partitioning scheme) or to the concept of chemical structure in the sense of Lewis structures, which divide the total molecular density into core and valence density, where the valence density is split up again into bonding and non-bonding electron densities. Examples are early and recent loge theories, the topological analysis by means of the Electron Localization Function (ELF), and the Natural Bond Orbital (NBO) approach. Of these partitioning schemes, the theories according to Bader (AIM), to Becke and Edgecomb (ELF) and according to Weinhold (NBO and Natural Resonance Theory, NRT), respectively, were reviewed in detail critically. Points of criticism were explicated for each of the mentioned theories. Since theoretically derived electron densities are to be compared to experimentally derived densities, a brief introduction into the theory of X-ray di®raction experiments was given and the multipole formalism was introduced. The procedure of density refinement was briefly discussed. Various suggestions for improvements were developed: One strategy would be the employment of model parameters, which are to a maximum degree mutually orthogonal, with the object of minimizing correlations among the model parameters, e.g., to introduce nodal planes into the radial functions of the multipole model. A further suggestion involves the guidance of the iterative refinement procedure by an extremum principle, which states, that when di®erent solutions to the least squares minimization problem are available with about the same statistical measures of quality and with about the same residual density, then the solution is to prefer, which yields a minimum density at the bond critical point (BCP) and a maximum polarity in terms of the ratio of distances between the BCP and the nuclei. This suggestion is based on the well known fact, that the bond polarity (in terms of the ratio of distances between the BCP and the respective nuclei) is underestimated in the experiment. Another suggestion for including physical constraints is the explicit consideration of the virial theorem, e.g., by evaluating the integration of the Laplacian over the entire atomic basins and comparing this value to zero and to the value obtained from the integration of the electron gradient field over the atomic surface. The next suggestion was to explicitly use the electrostatic theorem of Feynman (often also denoted as Hellmann-Feynman theorem), which states, that the forces onto the nuclei can be calculated from the purely classical electrostatic forces of the electron distribution and the nuclei distribution. For a stationary system, these forces must add to zero. This also provides an internal quality criterion of the density model. This can be performed in an iterative way during the refinement procedure or as a test of the final result. The use of the electrostatic theorem is expected to reduce significantly correlations among static density parameters and parameters describing vibrations, since it is a valuable tool to discriminate between physically reasonable and artificial static electron densities. All of these mentioned suggestions can be applied as internal quality criteria. The last suggestion is based on the idea to initiate the experimental refinement with a set of model parameters, which is, as much as possible close to the final solution. This can be achieved by performing periodic boundary conditions calculations, from which theoretically created files are obtained, which contain the Miller indices (h, k, l) and the respective intensity I. This file is used for a model parameter estimation (refinement), which excludes vibrations. The resulting parameters can be used for the experimental refinement, where, in a first step, the density parameters are fixed to determine the parameters describing vibrations. For a fine tuning, again the electrostatic theorem and the other above mentioned suggestions could be applied. Theoretical predictions should not be biased by the method of computation. Therefore the dependence of the density analyzing tools on the level of calculation (method of calculation/basis set) and on the substituents in complex chemical bonding situations were evaluated in the second part of the present work. A number of compounds containing formal single and double sulfur nitrogen bonds was investigated. For these compounds, experimental data were also available. The calculated data were compared internally and with the experimental results. The internal comparison was drawn with regard to questions of convergency as well as with regard to questions of consistency: The resulting molecular properties from NBO/NRT analyses were found to be very stable, when the geometries were optimized at the respective level of theory. This stability is valid for variations in the methods of calculation as well as for variations in the basis set. Only the individual resonance weights of the contributing Natural Lewis Structures differed considerably depending on the level of calculation and depending on the substituents. However, the deviations were in both cases to a large extent within a limit which preserves the descending order of the leading resonance structure weights. The resulting bond orders, i.e., the total, covalent and ionic bond order from NRT calculations, were not affected by the shift in the resonance weights. The analysis of the bond topological parameters resulted in a discrimination between insensitive parameters and sensitive parameters. The stable parameters do neither depend strongly on the method of calculation nor on the basis set. Only minor variation occurs in the numerical values of these parameters, when the level of calculation is changed or even when other functional groups (H, Me, or tBu) are employed, as long as the methods of calculation do not drop considerably below a standard level. The bond descriptors of the sulfur nitrogen bonds were found to be also stable with respect to the functional groups R = H, R = Me, and R = tBu. Stable parameters are the bond distance, the density at the bond critical point (BCP) and the ratio of distances between the BCP and the nuclei A and B, which varies clearly when considering the formal bond type. For very small basis sets like the 3-21G basis set, this characteristic stability collapses. The sensitive parameters are based on the second derivatives of the density with respect to the coordinates. This is in accordance with the well known fact, that the total second derivative of the density with respect to the coordinates is a strongly oscillating function with positive as well as negative values. A profound deviation has to be anticipated as a consequence of strong oscillations. lambda3, which describes the local charge depletion in the direction of the interaction line, is the most varying parameter. A detailed analysis revealed that the position of the BCP in the rampant edge of the Laplacian distribution is responsible for the sensitivity of the numerical value of lambda3 in formal double bonds. Since the slope of the Laplacian assumes very high values in its rampant edge, a tiny displacement of the BCP leads already to a considerable change in lambda3. This instability is not a failure of the underlying theory, but it yields de facto to a considerable dependence of sensitive bond topological properties on the method of calculation and on the applied basis sets. Since the total second derivative is important to judge on the nature of the bond in the AIM theory (closed shell interactions versus shared interactions), the changes in lambda3 can lead to differing chemical interpretations. The comparison of theoretically derived bond topological properties of various sulfur nitrogen bonds provides the possibility to measure the self consistency of this data set. All data sets clearly exhibit a linear correlation between the bond distances and the density at the BCP on one hand and between the bond distances and the Laplacian values at the BCP on the other hand. These correlations were almost independent of the basis set size. In this context, the linear regression has to be regarded exclusively as a descriptive statistics tool. There is no correlation anticipated a priori. The formal bond type was found to be readily deducible from the theoretically obtained bond topological descriptors of the model systems. In this sense, the bond topological properties are self consistent despite of the numerical sensitivity of the derivatives, as exemplified above. Often, calculations are performed with the experimentally derived equilibrium geometries and not with optimized ones. Applying this approach, the computationally costly geometry optimizations are saved. Following this approach the bond topological properties were calculated using very flexible basis sets and employing the fixed experimental geometry (which, of course, includes the application of tBu groups). Regression coe±cients similar to those from optimized geometries were obtained for correlations between bond distances and the densities at the BCP as well as for the correlation between bond distances and the Laplacian at the BCP, i.e. the approach is valid. However, the data points scattered less and the coe±cient of correlation was clearly increased when geometry optimizations were performed beforehand. The comparison between data obtained from theory and experiment revealed fundamental discrepancies: In the data set of bond topological parameters from the experiment, the behavior of only 2 out of 3 insensitive parameters was comparable to the behavior of the theoretically obtained values, i.e. theoretical and experimental bond distances as well as theoretical and experimental densities at the BCP correlate. From the theoretically obtained data it was easy to deduce the formal bond type from the position of the BCP, since it changed in a systematic manner. The respective experimentally obtained values were almost constant and did not change systematically. For the SN bonds containing compounds, the total second derivative assumes exclusively negative values in the experiment. Due to the different internal behavior, experimentally and theoretically sensitive bond topological values could not be compared directly. The qualitative agreement in the Laplacian distribution, however, was excellent. In the third and last part of this work, the application to chemical systems follows. Formal hypervalent molecules, i.e. molecules where some atoms are considered to hold more than 8 electrons in their valence shell, were investigated. These were compounds containing sulfur nitrogen bonds (H(NtBu)2SMe, H2C{S(NtBu)2(NHtBu)}2, S(NtBu)2 and S(NtBu)3) and a highly coordinated silicon compound. The set of sulfur nitrogen compounds also contained a textbook example for valence expansion, the sulfur triimide. For these molecules, experimental reference values were available from high resolution X-ray experiments. The experimental results were in the case of the sulfur triimide not unique. Furthermore, from the experimental bond topological data no definite conclusion about the formal bonding type could be drawn. The situation of sulfur nitrogen bonds in the above mentioned set of molecules was analyzed in terms of a geometry discussion and by means of a topological analysis. The methyl-substituted isolated molecules served as model compounds. For the interpretation of the bonding situation additional NBO/NRT calculations were preformed for the sulfur nitrogen compounds and an ELF calculation and analysis was performed for the silicon compound. The ELF analysis included not only the presentation and discussion of the ELF-isosurfaces (eta = 0.85), but also the investigation of populations of disynaptic valence basins and the percentage contributions to these populations of the individual atoms when the disynaptic valence basins are split into atomic contributions according to Bader’s partitioning scheme. The question of chemical interest was whether hypervalency is present in the set of molecules or not. In the first case the octet rule would be violated, in the second case Pauling’s verdict would be violated. While the concept of hypervalency is well established in chemistry, the violation of Pauling’s verdict is not. The quantitative numbers of the sensitive bond topological values from theory and experiment were not comparable, since no systematic relationship between the experimentally and theoretically determined sensitive bond descriptors was found. However, the insensitive parameters are in good agreement and the qualitative Laplacian distribution is, with few exceptions, in excellent agreement. The formal bonding type was deduced from experimental and theoretical topological data by considering the number and shape of valence shell charge concentrations in proximity to the sulfur and nitrogen centers. The results from NBO/NRT calculations confirmed the findings. All employed density analyzing tools AIM, ELF and NBO/NRT coincided in describing the bonding situation in the formally hypervalent molecules as highly polar. A comparison and analysis of experimentally and theoretically derived electron densities led consistently to the result, that regarding this set of molecules, hypervalency has to be excluded unequivocally.
The present work consists of two parts. The first one deals with theoretical questions and tests the performance of orbitals obtained from a self-interaction free KS method, the LHFapproach, in multireference ab initio methods. The purpose of this part is to enable a more efficient computation of excitation energies, which is important for the spectroscopic characterization of many organic and bioorganic molecules. The second part focuses on bioorganic questions and studies the base pairing properties of the purine base xanthine in order to explain, e.g., the unusually high stability of selfpairing xanthine alanyl-PNA double strands and the mutagenicity of xanthine formed in DNA. Part1: In contrast to HF- and standard DFT-methods, the LHF-approach leads to a fully bound virtual orbital spectrum, because Coulomb self interactions are exactly canceled in the LHFansatz. Furthermore, the energies of the occupied orbitals are not upshifted, like it is the case for standard DFT-methods, so that Koopmans' theorem remains valid. In line with this, also the occupied LHF-orbitals are somewhat more compact than standard DFT-orbitals. The present work shows that both properties are of great benefit for MR methods. The virtual LHF-orbitals are well optimized and allow an efficient description of excited states and static correlation in both MRCI- and MRPT2-approaches. Furthermore, the higher compactness of the occupied LHF- compared to standard DFT-orbitals leads to a better description of the center ion of Rydberg states. However, for each of the two advantages mentioned at least one example molecule has been found, for which LHF-orbitals actually perform worse than HF-and/or standard DFT-orbitals. This shows, that even though LHF virtual orbitals allow an excellent MRCI- and MRPT2-description for the electronically excited states of a large number of molecules, this cannot be generalized and their performance needs to be tested for each individual case. In the second part of the present work, the base pairing properties of xanthine and xanthine derivatives were studied. The purpose of this part was to find an explanation for the unexpectedly high stability of the xanthine alanyl PNA double strand. Furthermore, it was analyzed, why xanthine, that is formed from guanine in DNA under chemical stress, is able to form mismatched base pairs with the pyrimidine base thymine. Stability of xanthine alanyl PNA: In the first step, the regioisomer present in the considered alanyl PNA was identified to be the N7-regioisomer of xanthine by a theoretical analysis of the 13C-NMR spectrum. To analyze the stability of the xanthine self-pairing, a simplified model was set up, in which the stability of the PNA double strand was explained solely by the energy contributions from H-bonding and base stacking. For that purpose, the dimerization and stacking energies for the xanthine-xanthine, guaninecytosine, adenine-thymine and xanthine-2,6-diaminopurine base pairs were computed using DFT and MP2 methods. Solvent effects were taken into account by the conductor like screening model. The influence of the peptide backbone on the stacking geometry was considered by force field optimizations. While the individual contributions from hydrogen bonding and stacking do not correlate with the melting temperature Tm, the sum of both correlates linearly with Tm. This correlation is somewhat surprising, because this means that the effects of the entropy and the molecular water environment either cancel or are similar for all systems compared. In this model, the stability of the xanthine selfpairing mainly stems from an enlarged stacking interaction, while the H-bonds give only minor contributions to the stability of the xanthine selfpaired double strand of alanyl-PNA. Base pairing properties of N9-Xanthine: The computation of the base pairing properties of N9-xanthine revealed a strong variation in the individual H-bond strengths for the selfpairing of xanthine, that range from -4 to -11 kcal/mol in the gas phase and -2.5 to -5 kcal/mol in polar solvent. By comparison with model systems it was shown that the strong variance of the H-bond strength is mainly due to attractive or repulsive secondary electrostatic interactions. For the homodimer of hypoxanthine it was shown that the increase of aromaticity in the pyrimidine ring upon dimer formation leads to a strengthening of the hydrogen bonds. Mutagenicity of hypoxanthine and xanthine: Several neutral and anionic Watson-Crick base pairs of xanthine were computed with MP2- and DFT-methods in order to explain the mutagenicity of hypoxanthine and xanthine. Also basepairs involving tautomeric forms of xanthine and hypoxanthine were considered. To evaluate the dimerization energies found, the dimers were classified into pairings that have the exact geometry of the canonical base pairs and those that realize a distorted Watson-Crick pairing mode. The computations show that a stable pairing which realizes the exact geometry of a canonical Watson Crick base pairing is only possible for the pairing of xanthine to cytosine, however, the base pairs are only weakly bound. The dimerization energies of both the neutral and the anionic pairing is around 0 kcal/mol, so that the xanthine-cytosine base pairs are incorporated into DNA solely because the base pairs fulfill the geometric demands of DNA polymerase, but it does not profit from any additional stabilization due to hydrogen bonding. The bonding that in the Watson-Crick pairing mode xanthine has almost no affinity to cytosine is in correspondence with the experimental result that the cytosine-xanthine base pair is incorporated into DNA at a much lower rate than the cytosine-guanine base pair, which has a very strong hydrogen bonding. While the affinity of xanthine to cytosine is very low, the computations predict that xanthine is able to form a stable Watson-Crick pairing with thymine. However, the pairing has a somewhat distorted Watson-Crick geometry, so that its high stability is outbalanced by the worsened fit to the binding pocket of DNA-polymerase. As a consequence, the xanthinethymine pairing is incorporated into DNA not at a faster, but only at a rate comparable to that of the xanthine-cytosine pairing.
The first part of this work focuses on the characterization of systems which complex electronic structures require the application of multi-reference methods. The anti-tumor efficacy of the natural product Neocarzinostatin is based on the formation of diradicals and causes DNA cleavage and finally cytolysis. Computations on model systems performed in the present work show the influence of structural features on the mode of action and the efficacy of this antitumor-antibiotic. The cyclization of systems related to the enyne-cumulene framework like the enyne-allenes was investigated earlier and relations to the more unusual class of enyne-ketenes are analyzed. The class of enyne-ketenes (and also the enyne-allenes) show a broad spectrum of possible intermediates (diradicals, zwitterions, allenes). The electronic structures of these intermediates are also possible for the (heteroatom substituted) 1,2,4-cyclohexatriene and a model for their energetic sequence based on high-level multi-reference computations is proposed. In all three projects the application of multi-reference approaches is necessary to obtain a comprehensive picture of the reactivity and electronic structure but also shows up the limits inherently existing in the currently available programs with respect to the size of the molecules. In the second part, algorithms for a multi-reference Moller-Plesset perturbation theory (MR-MP2) program, designed to perform large-scale computations, were developed and implemented. The MR-MP2 approach represents the most cost-effective multireference ansatz and requires an efficient evaluation of the Hamilton matrix for which an algorithm is designed to instantly recognize only non-vanishing matrix elements and to employ the recurring interaction patterns of the Hamilton matrix. The direct construction of the Hamilton matrix is additionally parallelized to work on cluster environments.