Celia Fonseca Guerra

Find an error

Name: Célia Fonseca Guerra
Organization: VU University , Belgium
Department: Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling
Title: (PhD)

TOPICS

Co-reporter:Biswarup Jash;Philipp Scharf;Nikolas Sandmann;Célia Fonseca Guerra;Dominik A. Megger;Jens Müller
Chemical Science (2010-Present) 2017 vol. 8(Issue 2) pp:1337-1343
Publication Date(Web):2017/01/30
DOI:10.1039/C6SC03482A
A nucleoside analogue comprising the ligand 1H-imidazo[4,5-f][1,10]phenanthroline (P) was applied to develop a molecular beacon capable of discriminating the canonical nucleobases cytosine and thymine. The beacon is based on the formation of a stable Ag+-mediated base pair between P and cytosine, whereas the presence of Ag+ strongly destabilizes nucleic acids comprising an artificial base pair between P and thymine. Metal-mediated base pair formation was investigated by temperature-dependent UV spectroscopy and CD spectroscopy and complemented by extensive DFT calculations. The molecular beacon significantly extends the application spectrum of nucleic acids with metal-mediated base pairs. It is of potential use in the detection of single-nucleotide polymorphisms.
Co-reporter:Abril C. Castro;Marcel Swart;Célia Fonseca Guerra
Physical Chemistry Chemical Physics 2017 vol. 19(Issue 21) pp:13496-13502
Publication Date(Web):2017/05/31
DOI:10.1039/C7CP00397H
In the present study, we have theoretically analyzed supramolecular complexes based on the Watson–Crick A–T and A–U base pairs using dispersion-corrected density functional theory (DFT). Hydrogen atoms H8 and/or H6 in the natural adenine and thymine/uracil bases were replaced, respectively, by substituents X8, Y6 = NH−, NH2, NH3+ (N series), O−, OH, OH2+ (O series), F, Cl or Br (halogen series). We examined the effect of the substituents on the hydrogen-bond lengths, strength and bonding mechanism, and the NMR shielding constants of the C2-adenine and C2-thymine/uracil atoms in the base pairs. The general belief in the literature that there is a direct connection between changes in the hydrogen-bond strength and the C2-adenine shielding constant is conclusively rejected by our computations.
Co-reporter:Gábor Paragi;Zoltán Kupihár;Gábor Endre;Célia Fonseca Guerra;Lajos Kovács
Organic & Biomolecular Chemistry 2017 vol. 15(Issue 10) pp:2174-2184
Publication Date(Web):2017/03/08
DOI:10.1039/C6OB02574A
5-Substituted uracils (NH2 or OH groups in position 5) have been examined theoretically and experimentally as potential building blocks in quadruplex structures. Our high level Density Functional Theory (DFT) calculations showed that the tetramer formation and stacking energies for 5-substituted uracils are similar to the energies of purine-based xanthine (X) or guanine (G) structures. As tetrads of 5-substituted uracils cover almost exactly the same area as purine tetrads, mixed tetrads or quadruplex structures based on X or G and 5-substituted uracil motifs are possible. According to the calculations, 5-hydroxyuracil-based structures are the best candidates for experimental implementation which was corroborated by the existence of higher complexes in the mass spectra of 1-benzyl-5-hydroxyuracil. These pyrimidine-based molecules can be used as efficient building blocks in different applications including aptamers, bio-sensors or – taking into account the larger cavity in the central region of 5-hydroxyuracil structures – as an artificial ion channel.
Co-reporter:F. Zaccaria, G. Paragi and C. Fonseca Guerra  
Physical Chemistry Chemical Physics 2016 vol. 18(Issue 31) pp:20895-20904
Publication Date(Web):03 May 2016
DOI:10.1039/C6CP01030J
The alkali metal ion affinity of guanine quadruplexes has been studied using dispersion-corrected density functional theory (DFT-D). We have done computational investigations in aqueous solution that mimics artificial supramolecular conditions where guanine bases assemble into stacked quartets as well as biological environments in which telomeric quadruplexes are formed. In both cases, an alkali metal cation is needed to assist self-assembly. Our quantum chemical computations on these supramolecular systems are able to reproduce the experimental order of affinity of the guanine quadruplexes for the cations Li+, Na+, K+, Rb+, and Cs+. The strongest binding is computed between the potassium cation and the quadruplex as it occurs in nature. The desolvation and the size of alkali metal cations are thought to be responsible for the order of affinity. Until now, the relative importance of these two factors has remained unclear and debated. By assessing the quantum chemical ‘size’ of the cation, determining the amount of deformation of the quadruplex needed to accommodate the cation and through the energy decomposition analysis (EDA) of the interaction energy between the cation and the guanines, we reveal that the desolvation and size of the alkali metal cation are both almost equally responsible for the order of affinity.
Co-reporter:O. A. Stasyuk, H. Szatylowicz, T. M. Krygowski and C. Fonseca Guerra  
Physical Chemistry Chemical Physics 2016 vol. 18(Issue 17) pp:11624-11633
Publication Date(Web):15 Jan 2016
DOI:10.1039/C5CP07483E
The substituent effect of the amino and nitro groups on the electronic system of benzene has been investigated quantum chemically using quantitative Kohn–Sham molecular orbital theory and a corresponding energy decomposition analysis (EDA). The directionality of electrophilic substitution in aniline can accurately be explained with the amount of contribution of the 2pz orbitals on the unsubstituted carbon atoms to the highest occupied π orbital. For nitrobenzene, the molecular π orbitals cannot explain the regioselectivity of electrophilic substitution as there are two almost degenerate π orbitals with nearly the same 2pz contributions on the unsubstituted carbon atoms. The Voronoi deformation density analysis has been applied to aniline and nitrobenzene to obtain an insight into the charge rearrangements due to the substituent. This analysis method identified the orbitals involved in the C–N bond formation of the π system as the cause for the π charge accumulation at the ortho and para positions in the case of the NH2 group and the largest charge depletion at these same positions for the NO2 substituent. Furthermore, we showed that it is the repulsive interaction between the πHOMO of the phenyl radical and the πHOMO of the NH2 radical that is responsible for pushing up the πHOMO of aniline and therefore activating this π orbital of the phenyl ring towards electrophilic substitution.
Co-reporter:Matú&x161; Durec;Francesco Zaccaria;Dr. Célia FonsecaGuerra;Dr. Radek Marek
Chemistry - A European Journal 2016 Volume 22( Issue 31) pp:10912-10922
Publication Date(Web):
DOI:10.1002/chem.201601608

Abstract

Repetitive guanine-rich nucleic acid sequences play a crucial role in maintaining genome stability and the cell life cycle and represent potential targets for regulatory drugs. Recently, it has been demonstrated that guanine-based ligands with a porphyrin core can be used as markers of G-quadruplex assemblies in cell tissues. Herein, model systems of guanine-based ligands are explored by DFT methods. The energies of formation of modified guanine tetrads and those of modified tetrads stacked on the top of natural guanine tetrads have been calculated. The interaction energy has been decomposed into contributions from hydrogen bonding, stacking, and ion coordination and a twist–rise potential energy scan has been performed to find the individual local minima. Energy decomposition analysis reveals the impact of various substituents (F, Cl, Br, I, Me, NMe2) on individual energy terms. In addition, cooperative reinforcement in forming the modified and stacked tetrads, as well as the frontier orbitals participating in the hydrogen-bonding framework involving the HOMO–LUMO gap between the occupied σHOMO on the proton-accepting C=O and =N− groups and unoccupied σLUMO on the N−H groups, has been studied. The investigated systems are demonstrated to have a potential in ligand development, mainly due to stacking enhancement compared with natural guanine, which is used as a reference.

Co-reporter:Raquel Gavara, Elisabet Aguiló, Célia Fonseca Guerra, Laura Rodríguez, and João Carlos Lima
Inorganic Chemistry 2015 Volume 54(Issue 11) pp:5195-5203
Publication Date(Web):May 20, 2015
DOI:10.1021/acs.inorgchem.5b00025
The complexes [Au(4-pyridylethynyl)(phosph)] (phosph = PTA (1), DAPTA (2)) are known to produce supramolecular aggregates and gels in water. We studied the impact of these aggregation processes in the absorption spectra, 1H NMR (at different temperatures and concentrations), and DLS and estimated the equilibrium constant for a single step aggregation of the molecule (K = 26760 and 2590 M–1 for 1 and 2, respectively, at 25 °C). We present spectroscopic evidence for the presence of Au···Au contacts in the aggregates: the recorded changes on 1H NMR and the appearance of new absorption bands assigned to (σ*Au···Au–π*) have been attributed to the short (Au···Au) average distances in the aggregates. Relativistic density functional theory computations support the existence of short Au···Au distances and reveal charge-transfer in the aurophilic interactions. The free energy for a single step aggregation was calculated from the experimental data, and the value obtained (ΔG ∼ −20 kJ/mol) is in good agreement with the expected values in the order of the energies found for hydrogen bonds. The DFT computations confirm the experimental findings that aggregation of monomer 1 is stronger than the aggregation of monomer 2 and the existence of aurophilic interactions.
Co-reporter:Lando P. Wolters, Nicole W. G. Smits and Célia Fonseca Guerra  
Physical Chemistry Chemical Physics 2015 vol. 17(Issue 3) pp:1585-1592
Publication Date(Web):12 Nov 2014
DOI:10.1039/C4CP03740E
Halogen bonds are shown to possess the same characteristics as hydrogen bonds: charge transfer, resonance assistance and cooperativity. This follows from the computational analyses of the structure and bonding in N-halo-base pairs and quartets. The objective was to achieve an understanding of the nature of resonance-assisted halogen bonds (RAXB): how they resemble or differ from the better understood resonance-assisted hydrogen bonds (RAHB) in DNA. We present an accurate physical model of the RAXB based on the molecular orbital theory, which is derived from the corresponding energy decomposition analyses and study of the charge distribution. We show that the RAXB arise from classical electrostatic interaction and also receive strengthening from donor–acceptor interactions within the σ-electron system. Similar to RAHB, there is also a small stabilization by π-electron delocalization. This resemblance leads to prove cooperativity in N-halo-guanine quartets, which originates from the charge separation that occurs with donor–acceptor orbital interactions in the σ-electron system.
Co-reporter:Laia Guillaumes;Dr. Sílvia Simon;Dr. Célia FonsecaGuerra
ChemistryOpen 2015 Volume 4( Issue 3) pp:
Publication Date(Web):
DOI:10.1002/open.201500022

Abstract

Invited for this issues cover are Dr. Célia Fonseca Guerra from the VU University of Amsterdam and her collaborators at the University of Girona. The cover picture shows H-bonds in the adenine–thymine Watson–Crick base pair. An essential part of these H-bonds is their covalent component arising from donor–acceptor interactions between N or O lone pairs and the N−H antibonding σ* acceptor orbital. This charge-transfer interaction is represented by green figures walking on the pedestrian crossing, connecting the bases. This covalent component is the reason why H-bonds between DNA and/or unsaturated model bases are significantly stronger than those between analogous saturated bases. This contrasts sharply with the classical picture of predominantly electrostatic H-bonds which is not only incomplete in terms of a proper bonding mechanism, but also fails to explain the trend in stability. For more details, see the Full Paper on p. 318 ff.

Co-reporter:L. Guillaumes;Dr. S. Simon;Dr. C. FonsecaGuerra
ChemistryOpen 2015 Volume 4( Issue 3) pp:
Publication Date(Web):
DOI:10.1002/open.201580301
Co-reporter:Indranil Sinha;Dr. Célia FonsecaGuerra;Dr. Jens Müller
Angewandte Chemie 2015 Volume 127( Issue 12) pp:3674-3677
Publication Date(Web):
DOI:10.1002/ange.201411931

Abstract

Erstmals wurde eine parallele DNA-Doppelhelix mit Hoogsteen-Basenpaarung erzeugt, die bereitwillig ein Ag+-Ion in ein internes Fehlpaar unter Bildung eines metallvermittelten Basenpaares aufnimmt. Hierzu wurde das hochgradig stabilisierende 6 FP-Ag+-6 FP-Basenpaar entwickelt, das die künstliche Nucleobase 6-Furylpurin (6 FP) enthält. Eine Kombination aus temperaturabhängiger UV-Spektroskopie, CD-Spektroskopie und DFT-Rechnungen wurde angewendet, um die Bildung dieses Basenpaares zu bestätigen. Die Nucleobase 6 FP kann metallvermittelte Basenpaare sowohl über die Watson-Crick-Seite (d. h. in regulärer antiparalleler DNA) als auch über die Hoogsteen-Seite (d. h. in paralleler DNA) bilden, abhängig von der Oligonucleotidsequenz und den experimentellen Bedingungen. Bei dem 6 FP-Ag+-6 FP-Basenpaar in paralleler DNA handelt es sich um das am stärksten stabilisierende Ag+-vermittelte Basenpaar in einer Nucleinsäure, über das bisher berichtet wurde, mit einem Anstieg der Schmelztemperatur um fast 15 °C beim Binden eines Ag+-Ions.

Co-reporter:Indranil Sinha;Dr. Célia FonsecaGuerra;Dr. Jens Müller
Angewandte Chemie International Edition 2015 Volume 54( Issue 12) pp:3603-3606
Publication Date(Web):
DOI:10.1002/anie.201411931

Abstract

The first parallel-stranded DNA duplex with Hoogsteen base pairing that readily incorporates an Ag+ ion into an internal mispair to form a metal-mediated base pair has been created. Towards this end, the highly stabilizing 6 FP-Ag+-6 FP base pair comprising the artificial nucleobase 6-furylpurine (6 FP) was devised. A combination of temperature-dependent UV spectroscopy, CD spectroscopy, and DFT calculations was used to confirm the formation of this base pair. The nucleobase 6 FP is capable of forming metal-mediated base pairs both by the Watson–Crick edge (i.e. in regular antiparallel-stranded DNA) and by the Hoogsteen edge (i.e. in parallel-stranded DNA), depending on the oligonucleotide sequence and the experimental conditions. The 6 FP-Ag+-6 FP base pair within parallel-stranded DNA is the most strongly stabilizing Ag+-mediated base pair reported to date for any type of nucleic acid, with an increase in melting temperature of almost 15 °C upon the binding of one Ag+ ion.

Co-reporter:L. Guillaumes;Dr. S. Simon;Dr. C. FonsecaGuerra
ChemistryOpen 2015 Volume 4( Issue 3) pp:318-327
Publication Date(Web):
DOI:10.1002/open.201402132

Abstract

Hydrogen bonds play a crucial role in many biochemical processes and in supramolecular chemistry. In this study, we show quantum chemically that neither aromaticity nor other forms of π assistance are responsible for the enhanced stability of the hydrogen bonds in adenine–thymine (AT) DNA base pairs. This follows from extensive bonding analyses of AT and smaller analogs thereof, based on dispersion-corrected density functional theory (DFT). Removing the aromatic rings of either A or T has no effect on the Watson–Crick bond strength. Only when the smaller mimics become saturated, that is, when the hydrogen-bond acceptor and donor groups go from sp2 to sp3, does the stability of the resulting model complexes suddenly drop. Bonding analyses based on quantitative Kohn–Sham molecular orbital theory and corresponding energy decomposition analyses (EDA) show that the stronger hydrogen bonds in the unsaturated model complexes and in AT stem from stronger electrostatic interactions as well as enhanced donor–acceptor interactions in the σ-electron system, with the covalency being responsible for shortening the hydrogen bonds in these dimers.

Co-reporter:Erica C.M. Ording-Wenker, Martijn van der Plas, Maxime A. Siegler, Sylvestre Bonnet, F. Matthias Bickelhaupt, Célia Fonseca Guerra, and Elisabeth Bouwman
Inorganic Chemistry 2014 Volume 53(Issue 16) pp:8494-8504
Publication Date(Web):August 4, 2014
DOI:10.1021/ic501060w
The redox equilibrium between dinuclear CuII μ-thiolate and CuI disulfide structures has been analyzed experimentally and via DFT calculations. Two new ligands, L2SSL2 and L4SSL4, and their CuII μ-thiolate and CuI disulfide complexes were synthesized. For L2SSL2, these two redox-isomeric copper species are shown to be in equilibrium, which depends on both temperature and solvent. For L4SSL4 the μ-thiolate species forms as the kinetic product and further evolves into the disulfide complex under thermodynamic control, which creates the unprecedented possibility to compare both species under the same reaction conditions. The energies of the μ-thiolate and disulfide complexes for two series of related ligands have been calculated with DFT; the results rationalize the experimentally observed structures, and emphasize the important role that steric requirements play in the formation of the CuII thiolate structure.
Co-reporter:Jordi Poater, Marcel Swart, F. Matthias Bickelhaupt and Célia Fonseca Guerra  
Organic & Biomolecular Chemistry 2014 vol. 12(Issue 26) pp:4691-4700
Publication Date(Web):20 May 2014
DOI:10.1039/C4OB00427B
We have computationally investigated the structure and stability of B-DNA. To this end, we have analyzed the bonding in a series of 47 stacks consisting of two base pairs, in which the base pairs cover the full range of natural Watson–Crick pairs, mismatched pairs, and artificial DNA base pairs. Our analyses provide detailed insight into the role and relative importance of the various types of interactions, such as, hydrogen bonding, π–π stacking interactions, and solvation/desolvation. Furthermore, we have analyzed the functionality of the twist-angle on the stability of the structure. Interestingly, we can show that all stacked base pairs benefit from a stabilization by 6 to 12 kcal mol−1 if stacked base pairs are twisted from 0° to 36°, that is, if they are mutually rotated from a congruent superposition to the mutually twisted stacking configuration that occurs in B-DNA. This holds especially for stacked AT pairs but also for other stacked base pairs, including GC. The electronic mechanism behind this preference for a twisted arrangement depends on the base pairs involved. We also show that so-called “diagonal interactions” (or cross terms) in the stacked base pairs are crucial for understanding the stability of B-DNA, in particular, in GC-rich sequences.
Co-reporter:Erica C. M. Ording-Wenker;Martijn van der Plas;Dr. Maxime A. Siegler;Dr. Célia FonsecaGuerra;Dr. Elisabeth Bouwman
Chemistry - A European Journal 2014 Volume 20( Issue 51) pp:16913-16921
Publication Date(Web):
DOI:10.1002/chem.201403918

Abstract

The proton-induced electron-transfer reaction of a CuII μ-thiolate complex to a CuI-containing species has been investigated, both experimentally and computationally. The CuII μ-thiolate complex [CuII2(LMeS)2]2+ is isolated with the new pyridyl-containing ligand LMeSSLMe, which can form both CuII thiolate and CuI disulfide complexes, depending on the solvent. Both the CuII and the CuI complexes show reactivity upon addition of protons. The multivalent tetranuclear complex [CuI2CuII2(LS)2(CH3CN)6]4+ crystallizes after addition of two equivalents of strong acid to a solution containing the μ-thiolate complex [CuII2(LS)2]2+ and is further analyzed in solution. This study shows that, upon addition of protons to the CuII thiolate compound, the ligand dissociates from the copper centers, in contrast to an earlier report describing redox isomerization to a CuI disulfide species that is protonated at the pyridyl moieties. Computational studies of the protonated CuII μ-thiolate and CuI disulfide species with LSSL show that already upon addition of two equivalents of protons, ligand dissociation forming [CuI(CH3CN)4]+ and protonated ligand is energetically favored over conversion to a protonated CuI disulfide complex.

Co-reporter:Erica C. M. Ording-Wenker;Martijn van der Plas;Dr. Maxime A. Siegler;Dr. Célia FonsecaGuerra;Dr. Elisabeth Bouwman
Chemistry - A European Journal 2014 Volume 20( Issue 51) pp:
Publication Date(Web):
DOI:10.1002/chem.201490211
Co-reporter:Dr. Célia FonsecaGuerra;Dr. Pablo J. SanzMiguel;Andrea Cebollada;Dr. F. Matthias Bickelhaupt;Dr. Bernhard Lippert
Chemistry - A European Journal 2014 Volume 20( Issue 31) pp:9494-9499
Publication Date(Web):
DOI:10.1002/chem.201403066

Abstract

The exocyclic amino groups of cytosine and adenine nucleobases are normally almost flat, with the N atoms essentially sp2 hybridized and the lone pair largely delocalized into the heterocyclic rings. However, a change to marked pyramidality of the amino group (N then sp3 hybridized, lone pair essentially localized at N) occurs during i) involvement of an amino proton in strong hydrogen bonding donor conditions or ii) with monofunctional metal coordination following removal of one of the two protons.

Co-reporter:Dr. Célia FonsecaGuerra;Dr. Zsolt Szekeres;Dr. F. Matthias Bickelhaupt
Chemistry - A European Journal 2011 Volume 17( Issue 32) pp:8816-8818
Publication Date(Web):
DOI:10.1002/chem.201101335
Co-reporter:Dr. Célia FonsecaGuerra;Hester Zijlstra;Dr. Gábor Paragi;Dr. F. Matthias Bickelhaupt
Chemistry - A European Journal 2011 Volume 17( Issue 45) pp:12612-12622
Publication Date(Web):
DOI:10.1002/chem.201102234

Abstract

We show that the cooperative reinforcement between hydrogen bonds in guanine quartets is not caused by resonance-assisted hydrogen bonding (RAHB). This follows from extensive computational analyses of guanine quartets (G4) and xanthine quartets (X4) based on dispersion-corrected density functional theory (DFT-D). Our investigations cover the situation of quartets in the gas phase, in aqueous solution as well as in telomere-like stacks. A new mechanism for cooperativity between hydrogen bonds in guanine quartets emerges from our quantitative Kohn–Sham molecular orbital (MO) and corresponding energy decomposition analyses (EDA). Our analyses reveal that the intriguing cooperativity originates from the charge separation that goes with donor–acceptor orbital interactions in the σ-electron system, and not from the strengthening caused by resonance in the π-electron system. The cooperativity mechanism proposed here is argued to apply, beyond the present model systems, also to other hydrogen bonds that show cooperativity effects.

Co-reporter:Gábor Paragi, Zoltán Kupihár, Gábor Endre, Célia Fonseca Guerra and Lajos Kovács
Organic & Biomolecular Chemistry 2017 - vol. 15(Issue 10) pp:NaN2184-2184
Publication Date(Web):2016/12/20
DOI:10.1039/C6OB02574A
5-Substituted uracils (NH2 or OH groups in position 5) have been examined theoretically and experimentally as potential building blocks in quadruplex structures. Our high level Density Functional Theory (DFT) calculations showed that the tetramer formation and stacking energies for 5-substituted uracils are similar to the energies of purine-based xanthine (X) or guanine (G) structures. As tetrads of 5-substituted uracils cover almost exactly the same area as purine tetrads, mixed tetrads or quadruplex structures based on X or G and 5-substituted uracil motifs are possible. According to the calculations, 5-hydroxyuracil-based structures are the best candidates for experimental implementation which was corroborated by the existence of higher complexes in the mass spectra of 1-benzyl-5-hydroxyuracil. These pyrimidine-based molecules can be used as efficient building blocks in different applications including aptamers, bio-sensors or – taking into account the larger cavity in the central region of 5-hydroxyuracil structures – as an artificial ion channel.
Co-reporter:Biswarup Jash, Philipp Scharf, Nikolas Sandmann, Célia Fonseca Guerra, Dominik A. Megger and Jens Müller
Chemical Science (2010-Present) 2017 - vol. 8(Issue 2) pp:NaN1343-1343
Publication Date(Web):2016/10/12
DOI:10.1039/C6SC03482A
A nucleoside analogue comprising the ligand 1H-imidazo[4,5-f][1,10]phenanthroline (P) was applied to develop a molecular beacon capable of discriminating the canonical nucleobases cytosine and thymine. The beacon is based on the formation of a stable Ag+-mediated base pair between P and cytosine, whereas the presence of Ag+ strongly destabilizes nucleic acids comprising an artificial base pair between P and thymine. Metal-mediated base pair formation was investigated by temperature-dependent UV spectroscopy and CD spectroscopy and complemented by extensive DFT calculations. The molecular beacon significantly extends the application spectrum of nucleic acids with metal-mediated base pairs. It is of potential use in the detection of single-nucleotide polymorphisms.
Co-reporter:Lando P. Wolters, Nicole W. G. Smits and Célia Fonseca Guerra
Physical Chemistry Chemical Physics 2015 - vol. 17(Issue 3) pp:NaN1592-1592
Publication Date(Web):2014/11/12
DOI:10.1039/C4CP03740E
Halogen bonds are shown to possess the same characteristics as hydrogen bonds: charge transfer, resonance assistance and cooperativity. This follows from the computational analyses of the structure and bonding in N-halo-base pairs and quartets. The objective was to achieve an understanding of the nature of resonance-assisted halogen bonds (RAXB): how they resemble or differ from the better understood resonance-assisted hydrogen bonds (RAHB) in DNA. We present an accurate physical model of the RAXB based on the molecular orbital theory, which is derived from the corresponding energy decomposition analyses and study of the charge distribution. We show that the RAXB arise from classical electrostatic interaction and also receive strengthening from donor–acceptor interactions within the σ-electron system. Similar to RAHB, there is also a small stabilization by π-electron delocalization. This resemblance leads to prove cooperativity in N-halo-guanine quartets, which originates from the charge separation that occurs with donor–acceptor orbital interactions in the σ-electron system.
Co-reporter:F. Zaccaria, G. Paragi and C. Fonseca Guerra
Physical Chemistry Chemical Physics 2016 - vol. 18(Issue 31) pp:NaN20904-20904
Publication Date(Web):2016/05/03
DOI:10.1039/C6CP01030J
The alkali metal ion affinity of guanine quadruplexes has been studied using dispersion-corrected density functional theory (DFT-D). We have done computational investigations in aqueous solution that mimics artificial supramolecular conditions where guanine bases assemble into stacked quartets as well as biological environments in which telomeric quadruplexes are formed. In both cases, an alkali metal cation is needed to assist self-assembly. Our quantum chemical computations on these supramolecular systems are able to reproduce the experimental order of affinity of the guanine quadruplexes for the cations Li+, Na+, K+, Rb+, and Cs+. The strongest binding is computed between the potassium cation and the quadruplex as it occurs in nature. The desolvation and the size of alkali metal cations are thought to be responsible for the order of affinity. Until now, the relative importance of these two factors has remained unclear and debated. By assessing the quantum chemical ‘size’ of the cation, determining the amount of deformation of the quadruplex needed to accommodate the cation and through the energy decomposition analysis (EDA) of the interaction energy between the cation and the guanines, we reveal that the desolvation and size of the alkali metal cation are both almost equally responsible for the order of affinity.
Co-reporter:O. A. Stasyuk, H. Szatylowicz, T. M. Krygowski and C. Fonseca Guerra
Physical Chemistry Chemical Physics 2016 - vol. 18(Issue 17) pp:NaN11633-11633
Publication Date(Web):2016/01/15
DOI:10.1039/C5CP07483E
The substituent effect of the amino and nitro groups on the electronic system of benzene has been investigated quantum chemically using quantitative Kohn–Sham molecular orbital theory and a corresponding energy decomposition analysis (EDA). The directionality of electrophilic substitution in aniline can accurately be explained with the amount of contribution of the 2pz orbitals on the unsubstituted carbon atoms to the highest occupied π orbital. For nitrobenzene, the molecular π orbitals cannot explain the regioselectivity of electrophilic substitution as there are two almost degenerate π orbitals with nearly the same 2pz contributions on the unsubstituted carbon atoms. The Voronoi deformation density analysis has been applied to aniline and nitrobenzene to obtain an insight into the charge rearrangements due to the substituent. This analysis method identified the orbitals involved in the C–N bond formation of the π system as the cause for the π charge accumulation at the ortho and para positions in the case of the NH2 group and the largest charge depletion at these same positions for the NO2 substituent. Furthermore, we showed that it is the repulsive interaction between the πHOMO of the phenyl radical and the πHOMO of the NH2 radical that is responsible for pushing up the πHOMO of aniline and therefore activating this π orbital of the phenyl ring towards electrophilic substitution.
Co-reporter:Abril C. Castro, Marcel Swart and Célia Fonseca Guerra
Physical Chemistry Chemical Physics 2017 - vol. 19(Issue 21) pp:NaN13502-13502
Publication Date(Web):2017/04/25
DOI:10.1039/C7CP00397H
In the present study, we have theoretically analyzed supramolecular complexes based on the Watson–Crick A–T and A–U base pairs using dispersion-corrected density functional theory (DFT). Hydrogen atoms H8 and/or H6 in the natural adenine and thymine/uracil bases were replaced, respectively, by substituents X8, Y6 = NH−, NH2, NH3+ (N series), O−, OH, OH2+ (O series), F, Cl or Br (halogen series). We examined the effect of the substituents on the hydrogen-bond lengths, strength and bonding mechanism, and the NMR shielding constants of the C2-adenine and C2-thymine/uracil atoms in the base pairs. The general belief in the literature that there is a direct connection between changes in the hydrogen-bond strength and the C2-adenine shielding constant is conclusively rejected by our computations.
Co-reporter:Jordi Poater, Marcel Swart, F. Matthias Bickelhaupt and Célia Fonseca Guerra
Organic & Biomolecular Chemistry 2014 - vol. 12(Issue 26) pp:NaN4700-4700
Publication Date(Web):2014/05/20
DOI:10.1039/C4OB00427B
We have computationally investigated the structure and stability of B-DNA. To this end, we have analyzed the bonding in a series of 47 stacks consisting of two base pairs, in which the base pairs cover the full range of natural Watson–Crick pairs, mismatched pairs, and artificial DNA base pairs. Our analyses provide detailed insight into the role and relative importance of the various types of interactions, such as, hydrogen bonding, π–π stacking interactions, and solvation/desolvation. Furthermore, we have analyzed the functionality of the twist-angle on the stability of the structure. Interestingly, we can show that all stacked base pairs benefit from a stabilization by 6 to 12 kcal mol−1 if stacked base pairs are twisted from 0° to 36°, that is, if they are mutually rotated from a congruent superposition to the mutually twisted stacking configuration that occurs in B-DNA. This holds especially for stacked AT pairs but also for other stacked base pairs, including GC. The electronic mechanism behind this preference for a twisted arrangement depends on the base pairs involved. We also show that so-called “diagonal interactions” (or cross terms) in the stacked base pairs are crucial for understanding the stability of B-DNA, in particular, in GC-rich sequences.
Benzo[lmn][3,8]phenanthroline-1,3,6,8(2H,7H)-tetrone, 4,9-dichloro-
Methanetricarbonitrile, fluoro-
Benzo[lmn][3,8]phenanthroline-1,3,6,8(2H,7H)-tetrone, 4,5,9,10-tetrachloro-
GERMANE, (1,1-DIMETHYLETHYL)METHYL-
Germane, (1,1-dimethylethyl)dimethyl-
Bismuthine, dimethyl-
Arsine, (1,1-dimethylethyl)dimethyl-
CHLORONIUM, (1,1-DIMETHYLETHYL)METHYL-
Argon(1+), methyl-
PROPANE, 2-METHYL-2-(METHYLTELLURO)-