Structure, stability, and dynamics of canonical and noncanonical base pairs: quantum chemical studies

The importance of non-Watson-Crick base pairs in the three-dimensional structure of RNA is now well established. The structure and stability of these noncanonical base pairs are, however, poorly understood. We have attempted to understand structural features of 33 frequently occurring base pairs using density functional theory. These are of three types, namely (i) those stabilized by two or more polar hydrogen bonds between the bases, (ii) those having one polar and another C-H...O/N type interactions, and (iii) those having one H-bond between the bases and another involving one of the sugars linked to the bases. We found that the base pairs having two polar H-bonds are very stable as compared to those having one C-H...O/N interaction. Our quantitatively analysis of structures of these optimized base pairs indicates that they possess a different amount of nonplanarity with large propeller or buckle values as also observed in the crystal structures. We further found that geometry optimization does not modify the hydrogen-bonding pattern, as values of shear and open angle of the base pairs remain conserved. The structures of initial crystal geometry and final optimized geometry of some base pairs having only one polar H-bond and a C-H...O/N interaction, however, are significantly different, indicating the weak nature of the nonpolar interaction. The base pair flexibility, as measured from normal-mode analysis, in terms of the intrinsic standard deviations of the base pair structural parameters are in conformity with those calculated from RNA crystal structures. We also noticed that deformation of a base pair along the stretch direction is impossible for all of the base pairs, and movements of the base pairs along shear and open are also quite restricted. The base pair opening mode through alteration of propeller or buckle is considerably less restricted for most of the base pairs.     

Structures for controlling synthesis design and reaction modelling via knowledge acquisition

The performance of programs for computer-assisted organic synthesis (CAOS) is directly related to the knowledge base and to the procedures for deciding how the knowledge base should be used. The architecture of a knowledge base (text base, data base, rule base, model base and commonsense base) is discussed. The problem of knowledge acquisition is described with regard to pruning the productivity of CAOS systems. The commercial program 1stCLASS (for symbolic processing) is used with the SAAO (devised here for numeric processing) to form a commonsense base. For this purpose, knowledge was acquired by induction from the solutions provided by human experts to given problems. The general procedure is applied to the recognition of reaction sites which are crucial in the simulation of organic reactions.     

Sequence dependence in base flipping: experimental and computational studies

Base flipping is the movement of a DNA base from an intrahelical, base-stacked position to an extrahelical, solvent-exposed position. As there are favorable interactions for an intrahelical base, both hydrogen bonding and base stacking, base flipping is expected to be energetically prohibitive for an undamaged DNA duplex. For damaged DNA bases, however, the energetic cost of base flipping may be considerably lower. Using a selective, non-covalent assay for base flipping, the sequence dependence of base flipping in DNA sequences containing an abasic site has been studied. The dissociation constants of the zinc-cyclen complex to small molecules and single strands of DNA as well as the equilibrium constants for base flipping have been determined for these sequences. Molecular dynamics simulations of the zinc-cyclen complex bound to both single- and double-stranded DNA have been performed in an attempt to rationalize the differences in the dissociation constants obtained for the two systems. The results are compared to previous studies of base flipping in DNA containing an abasic site.     

Effect of loop distortion on the stability and structural dynamics of DNA hairpin and dumbbell conjugates

The thermal stability and conformational dynamics of DNA hairpin and dumbbell conjugates having short A-tract base pair domains connected by tri- or hexa(ethylene glycol) linkers is reported. The formation of stable base-paired A-tract hairpins having oligo(ethylene glycol) linkers requires a minimum of four or five A-T base pairs. The formation of base-paired dumbbells having oligo(ethylene glycol) linkers by means of chemical ligation of nicked dumbbells requires a minimum of two A-T base pairs on either side of the nick. Molecular modeling indicates that the hexa(ethylene glycol) linker is sufficiently long to permit formation of strain-free loop regions and B-DNA base pair domains. In contrast, the tri(ethylene glycol) is too short to permit Watson-Crick base pairing between the bases attached to the linker. The shorter linker distorts the duplex, resulting in fluxional behavior in which the base pairs adjacent to the linker and at the open end of the hairpin dissociate on the nanosecond time scale. The loss of interstrand binding energy caused by these fluctuations leads to a difference of approximately 5 degrees C in melting temperature between EG3 and EG6 hairpins. An analysis of the fluxional behavior of the EG3 adjacent base-pair has been used to study the pathways for base flipping and base stacking, including the identification of rotated base (partially flipped) intermediates that have not been described previously for A-T base pairs.     

Vibrational dynamics of DNA. I. Vibrational basis modes and couplings

Carrying out density functional theory calculations of four DNA bases, base derivatives, Watson-Crick (WC) base pairs, and multiple-layer base pair stacks, we studied vibrational dynamics of delocalized modes with frequency ranging from 1400 to 1800 cm(-1). These modes have been found to be highly sensitive to structure fluctuation and base pair conformation of DNA. By identifying eight fundamental basis modes, it is shown that the normal modes of base pairs and multilayer base pair stacks can be described by linear combinations of these vibrational basis modes. By using the Hessian matrix reconstruction method, vibrational coupling constants between the basis modes are determined for WC base pairs and multilayer systems and are found to be most strongly affected by the hydrogen bonding interaction between bases. It is also found that the propeller twist and buckle motions do not strongly affect vibrational couplings and basis mode frequencies. Numerically simulated IR spectra of guanine-cytosine and adenine-thymine bases pairs as well as of multilayer base pair stacks are presented and described in terms of coupled basis modes. It turns out that, due to the small interlayer base-base vibrational interactions, the IR absorption spectrum of multilayer base pair system does not strongly depend on the number of base pairs.     

Mixed base of hydrophilic ointment and purified lanolin to improve the drug release rate and absorption of water of minocycline hydrochloride ointment for treatment of bedsores.

A desired ointment bases for better treatment of bedsores was developed to improve the release rate of minocycline hydrochloride (MH) and the water absorption capacity using various types of hydrophobic to hydrophilic ointment base. The influence of purified lanolin (PL) on the release behavior of MH from hydrophilic ointment (HO) base was primarily focused on. It was found that the release rate of drug increased with increase in the hydrophilicity of the base. A linear correlation between the apparent release rate constant of drug from the HO and PL mixed ointment base at various combination ratios and the elution of ointment base was noted. The HO ointment base containing 30% PL had the highest apparent release rate constant of MH. The mixed ointment base with the lowest viscosity showed the highest absorption of water and elution of ointment base. In conclusion, it was found that HO (70%) and PL (30%) mixed ointment base was a promising candidate for better treatment of bedsores.     

Ion trap collisional activation of the deprotonated deoxymononucleoside and deoxydinucleoside monophosphates

Deoxymononucleoside and deoxydinucleoside monophosphate anions formed by electrospray have been subjected to ion trap collisional activation. The threshold for decomposition via loss of base is significantly lower for the deoxymononucleoside 3′-monophosphates than for the corresponding 5′-monophosphates, which indicates that the presence of a charged 3′ phosphate group facilitates base loss. The behavior of the bases among each class of isomers shows slight variation in threshold and tandem mass spectrometry efficiency with tile notable exception of 2′-deoxyguanosine 5′-monophosphate. This ion is exceptionally stable toward decomposition via base loss, which reflects a strong hydrogen bonding interaction between the base and the phosphate group. All dinucleotides fragment via similar mechanisms, but the propensity for neutral base loss relative to loss of a charged base is highly dependent on the identities of both the 5′ and 3′ bases. The behavior of the dinucleotides under collisional activation conditions supports the proposal that base loss proceeds via a proton-bound dimer intermediate in which loss of the charged base directly competes with loss of the neutral base. Application of the kinetic method allows for quantitative predictions of the differences of the gas-phase acidities of the dimer components.     

Theoretical study of the scalar coupling constants across the noncovalent contacts in RNA base pairs: the cis- and trans-watson-crick/sugar edge base pair family

The structure and function of RNA molecules are substantially affected by non-Watson-Crick base pairs actively utilizing the 2'-hydroxyl group of ribose. Here we correlate scalar coupling constants across the noncovalent contacts calculated for the cis- and trans-WC/SE (Watson-Crick/sugar edge) RNA base pairs with the geometry of base to base and sugar to base hydrogen bond(s). 23 RNA base pairs from the 32 investigated were found in RNA crystal structures, and the calculated scalar couplings are therefore experimentally relevant with regard to the binding patterns occurring in this class of RNA base pairs. The intermolecular scalar couplings 1hJ(N,H), 2hJ(N,N), 2hJ(C,H), and 3hJ(C,N) were calculated for the N-H...N and N-H...O=C base to base contacts and various noncovalent links between the sugar hydroxyl and RNA base. Also, the intramolecular 1J(N,H) and 2J(C,H) couplings were calculated for the amino or imino group of RNA base and the ribose 2'-hydroxyl group involved in the noncovalent interactions. The calculated scalar couplings have implications for validation of local geometry, show specificity for the amino and imino groups of RNA base involved in the linkage, and can be used for discrimination between the cis- and trans-WC/SE base pairs. The RNA base pairs within an isosteric subclass of the WC/SE binding patterns can be further sorted according to the scalar couplings calculated across different local noncovalent contacts. The effect of explicit water inserted in the RNA base pairs on the magnitude of the scalar couplings was calculated, and the data for discrimination between the water-inserted and direct RNA base pairs are presented. The calculated NMR data are significant for structural interpretation of the scalar couplings in the noncanonical RNA base pairs.     

Factors determining the deriving force of DNA formation: geometrical differences of base pairs, dehydration of bases, and the arginine assisting

The mechanism of the fidelity synthesis of DNA associated with the process of dGTP combination to the DNA template was explored. The exclusion of water molecules from the hydrated DNA bases can amplify the energy difference between the correct and incorrect base pairs, but the effect of the water molecules on the Gibbs free energy of formation is dependent on the binding sites for the water molecules. The water detachment from the incoming dNTP is not the only factor but the first step for the successful replication of DNA. The second step is the selection of the DNA polymerase on the DNA base pair through the comparison between the correct DNA base and the incorrect DNA base. The bonding of the Arg668 with the incoming dNTP can enlarge the Gibbs free energies of formation of the base pairs, especially the correct base pairs, thus increasing the driving force of DNA formation. When the DNA base of the primer terminus is correct, the extension of the guanine and the adenine is quicker than that of the cytosine and the thymine because of the hydrogen bonding fork formation of Arg668 with the minor groove of the primer terminus and the ring oxygen of the deoxyribose moiety of the incoming dNTP. Because of the geometry differences of the incorrect base pairs with the correct base pairs, the effect from the DNA polymerase is smaller on the incorrect base pair than on the correct base pair, and the extension of a mispair is slower than that of a correct base pair. This decreases the extension rate of the base pair and thus allows proofreading exonuclease activity to excise the incorrect base pair. Arg668 cannot prevent the extension of the GT mispair, as well as the GC correct base pair, and GA and GG mispairs. This may be attributed to the small geometry difference between the GT base pair and the correct AT base pair.     

固体碱催化剂

综述了最近30年来固体碱催化剂的研究现状，包括固体碱催化剂的种类、优缺点以及各类催化剂中影响其催化性能的因素。重点探讨了碱性前驱体和载体对负载型无机固化碱碱强度的影响，展望了未来固体碱剂的发展方向。     

Conformational specificity of non-canonical base pairs and higher order structures in nucleic acids: crystal structure database analysis

Non-canonical base pairs contribute immensely to the structural and functional variability of RNA, which calls for a detailed characterization of their spatial conformation. Intra-base pair parameters, namely propeller, buckle, open-angle, stagger, shear and stretch describe structure of base pairs indicating planarity and proximity of association between the two bases. In order to study the conformational specificities of non-canonical base pairs occurring in RNA crystal structures, we have upgraded NUPARM software to calculate these intra-base pair parameters using a new base pairing edge specific axis system. Analysis of base pairs and base triples with the new edge specific axis system indicate the presence of specific structural signatures for different classes of non-canonical pairs and triples. Differentiating features could be identified for pairs in cis or trans orientation, as well as those involving sugar edges or C-H-mediated hydrogen bonds. It was seen that propeller for all types of base pairs in cis orientation are generally negative, while those for trans base pairs do not have any preference. Formation of a base triple is seen to reduce propeller of the associated base pair along with reduction of overall flexibility of the pairs. We noticed that base pairs involving sugar edge are generally more non-planar, with large propeller or buckle values, presumably to avoid steric clash between the bulky sugar moieties. These specific conformational signatures often provide an insight into their role in the structural and functional context of RNA.     

Stable, specific, and reversible base pairing via Schiff base

We here present a novel covalently linked base pair via Schiff base formation between 5-formyluracil (fU) and 5-aminocytosine (AmC). Formation of the Schiff base linkage proceeds reversibly and does not require any additives. The cross-linked DNA is very stable under denaturing conditions, whereas it completely dissociates upon heating at 90 degrees C. The pairing ability of AmC and fU is very specific and is applicable to the detection of fU, which is the major oxidative lesion of T in DNA. We propose the Schiff base linkage as a new artificial base pairing scheme to create functional DNAs.     

Synthesis of Fluorobenzene and Benzimidazole Nucleic‐Acid Analogues and Their Influence on Stability of RNA Duplexes

Six different ribonucleoside phosphoramidites with fluorobenzenes or fluorobenzimidazoles as base analogues, one abasic site, and inosine were synthesized and incorporated into oligoribonucleotides. The oligomers were investigated by means of UV and CD spectroscopy to assess the contribution of H‐bonding, base stacking, and solvation to the stability of the RNA duplex. CD Spectra show that the incorporation of modified nucleosides does not lead to changes in the structure of RNA. The T_m differences determined are based on changes in base stacking and solvation. Individual contributions of base stacking and solvation of the modified nucleosides could be determined. In fluorobenzene⋅fluorobenzimidazole‐modified base pairs, a duplex‐stabilizing force was found that points to a weak F⋅⋅⋅H H‐bond.     

Thermodynamic coupling of the loop and stem in unusually stable DNA hairpins closed by CG base pairs

For certain DNA hairpin loops, a CG closing base pair has enhanced stability over other closing base pairs, which cannot be explained by the current nearest-neighbor model. We report the use of three-carbon (C3) spacers to investigate the expandability of DNA hairpin loops and the coupling between the loop and closing base pair. Inserting the C3-spacers at most positions in these model loops produced only a modest stabilization or destabilization except for insertion between the 5' end of the loop and the CG closing base pair, which gave a large destabilization. Further investigation on tetraloops and triloops with other closing base pairs established that this destabilization is specific to the unusually stable CG closing base pair. Studies with the nucleotide analogues 2-aminopurine and 2,6-diaminopurine indicated that this stabilization may be due to coupling between functional groups on the first base of the loop and the CG closing base pair. The C3-spacers provide a simple way to interrupt potential interactions and thereby probe loop/stem coupling.     

Hydration and stability of nucleic acid bases and base pairs

Empirical, quantum chemical calculations and molecular dynamics simulations of the role of a solvent on tautomerism of nucleic acid bases and structure and properties of nucleic acid base pairs are summarized. Attention was paid to microhydrated (by one and two water molecules) complexes, for which structures found by scanning of empirical potential surfaces were recalculated at a correlated ab initio level. Additionally, isolated as well as mono- and dihydrated H-bonded, T-shaped and stacked structures of all possible nucleic acid base pairs were studied at the same theoretical levels. We demonstrate the strong influence of a solvent on the tautomeric equilibrium between the tautomers of bases and on the spatial arrangement of the bases in a base pair. The results provide clear evidence that the prevalence of either the stacked or hydrogen-bonded structures of the base pairs in the solvent is not determined only by its bulk properties, but rather by specific hydrophilic interactions of the base pair with a small number of solvent molecules.     

Location and base selectivity on fragmentation of brominated oligodeoxynucleotides

Bromine-modified oligodeoxynucleotides (ODNs) were fragmented in the electrospray source to study the influence of brominated bases on fragmentation. Several 13-mer ODNs containing a brominated pyrimidine base, BrdU (5-bromodeoxyuridine) or BrdC (5-bromodeoxycytidine), were analyzed. Low cone voltage fragmentation yields a loss of the brominated base with a preferential loss for the brominated base closer to the 5'-end (2-position > 4-position > 12-position) as well as a preferential loss of BrdU over BrdC. Higher cone voltage produces backbone fragmentation with complementary a(n)-base and w(m) ions close to the brominated base. On the basis of these observations, we located the brominated base in the sequence for all of the ODNs studied.     

Analysis of stacking overlap in nucleic acid structures: algorithm and application

RNA contains different secondary structural motifs like pseudo-helices, hairpin loops, internal loops, etc. in addition to anti-parallel double helices and random coils. The secondary structures are mainly stabilized by base-pairing and stacking interactions between the planar aromatic bases. The hydrogen bonding strength and geometries of base pairs are characterized by six intra-base pair parameters. Similarly, stacking can be represented by six local doublet parameters. These dinucleotide step parameters can describe the quality of stacking between Watson–Crick base pairs very effectively. However, it is quite difficult to understand the stacking pattern for dinucleotides consisting of non canonical base pairs from these parameters. Stacking interaction is a manifestation of the interaction between two aromatic bases or base pairs and thus can be estimated best by the overlap area between the planar aromatic moieties. We have calculated base pair overlap between two consecutive base pairs as the buried van der Waals surface between them. In general, overlap values show normal distribution for the Watson–Crick base pairs in most double helices within a range from 45 to 50 Å² irrespective of base sequence. The dinucleotide steps with non-canonical base pairs also are seen to have high overlap value, although their twist and few other parameters are rather unusual. We have analyzed hairpin loops of different length, bulges within double helical structures and pseudo-continuous helices using our algorithm. The overlap area analyses indicate good stacking between few looped out bases especially in GNRA tetraloop, which was difficult to quantitatively characterise from analysis of the base pair or dinucleotide step parameters. This parameter is also seen to be capable to distinguish pseudo-continuous helices from kinked helix junctions.     

Hydrogen-bonding interactions in peptide nucleic acid and deoxyribonucleic acid: a comparative study

Peptide nucleic acid (PNA) is a synthetic analogue of deoxyribonucleic acid (DNA) capable of tightly binding to itself and DNA with high specificity. Using hybrid density functional methods, hydrogen-bond (H-bond) strengths have been evaluated for isolated Watson-Crick base pairs, PNA base pairs, and charged as well as neutral DNA base pairs. Heterogeneous base pairs of PNA with charged and neutral DNA have also been investigated. The competing effects of short-range H-bonding and long-range Coulombic repulsions in charged DNA base pairs have been analyzed. Polarizable continuum models have been employed to evaluate solvation effects on the binding energies.     

Base pairing in RNA structures: A computational analysis of structural aspects and interaction energies

The base pairing patterns in RNA structures are more versatile and completely different as compared to DNA. We present here results of ab-initio studies of structures and interaction energies of eight selected RNA base pairs reported in literature. Interaction energies, including BSSE correction, of hydrogen added crystal geometries of base pairs have been calculated at the HF/6-31G** level. The structures and interaction energies of the base pairs in the crystal geometry are compared with those obtained after optimization of the base pairs. We find that the base pairs become more planar on full optimization. No change in the hydrogen bonding pattern is seen. It is expected that the inclusion of appropriate considerations of many of these aspects of RNA base pairing would significantly improve the accuracy of RNA secondary structure prediction.     

Polymer Supported Schiff Base Complexes of Iron(III), Cobalt(II) and Nickel(II) Ions and their Catalytic Activity in Oxidation of Phenol and Cyclohexene

The polymer supported transition metal complexes of N,N′‐bis (o‐hydroxy acetophenone) hydrazine (HPHZ) Schiff base were prepared by immobilization of N,N′‐bis(4‐amino‐o‐hydroxyacetophenone)hydrazine (AHPHZ) Schiff base on chloromethylated polystyrene beads of a constant degree of crosslinking and then loading iron(III), cobalt(II) and nickel(II) ions in methanol. The complexation of polymer anchored HPHZ Schiff base with iron(III), cobalt(II) and nickel(II) ions was 83.30%, 84.20% and 87.80%, respectively, whereas with unsupported HPHZ Schiff base, the complexation of these metal ions was 80.3%, 79.90% and 85.63%. The unsupported and polymer supported metal complexes were characterized for their structures using I.R, UV and elemental analysis. The iron(III) complexes of HPHZ Schiff base were octahedral in geometry, whereas cobalt(II) and nickel(II) complexes showed square planar structures as supported by UV and magnetic measurements. The thermogravimetric analysis (TGA) of HPHZ Schiff base and its metal complexes was used to analyze the variation in thermal stability of HPHZ Schiff base on complexation with metal ions. The HPHZ Schiff base showed a weight loss of 58% at 500°C, but its iron(III), cobalt(II) and nickel(II) ions complexes have shown a weight loss of 30%, 52% and 45% at same temperature. The catalytic activity of metal complexes was tested by studying the oxidation of phenol and epoxidation of cyclohexene in presence of hydrogen peroxide as an oxidant. The supported HPHZ Schiff base complexes of iron(III) ions showed 64.0% conversion for phenol and 81.3% conversion for cyclohexene at a molar ratio of 1∶1∶1 of substrate to catalyst and hydrogen peroxide, but unsupported complexes of iron(III) ions showed 55.5% conversion for phenol and 66.4% conversion for cyclohexene at 1∶1∶1 molar ratio of substrate to catalyst and hydrogen peroxide. The product selectivity for catechol (CTL) and epoxy cyclohexane (ECH) was 90.5% and 96.5% with supported HPHZ Schiff base complexes of iron(III) ions, but was found to be low with cobalt(II) and nickel(II) ions complexes of Schiff base. The selectivity for catechol (CTL) and epoxy cyclohexane (ECH) was different with studied metal ions and varied with molar ratio of metal ions in the reaction mixture. The selectivity was constant on varying the molar ratio of hydrogen peroxide and substrate. The energy of activation for epoxidation of cyclohexene and phenol conversion in presence of polymer supported HPHZ Schiff base complexes of iron(III) ions was 8.9 kJ mol^?1 and 22.8 kJ mol^?1, respectively, but was high with Schiff base complexes of cobalt(II) and nickel(II) ions and with unsupported Schiff base complexes.