Hydrogen bond equilibrium constants of some unusual nucleotide base pairs

Approximate hydrogen bond association constants were determined for base pairs formed by an adenine derivative and a number of unusual pyrimidine bases. A series is found in which the H-bond strength in the base-pairs varies. In certain cases the H-bond equilibrium constant is larger than in the adenine-thymine pair. Inosine derivatives seem to have a non-negligible chance of replacing guanosine in the guanosine-cytosine pair. Infrared, near-infrared (overtone) and NMR spectra were used to determine the equilibrium constants.     

Observation of H-bond mediated 3hJH2H3 coupling constants across Watson-Crick AU base pairs in RNA

^3hJ_H2H3trans-hydrogen bond scalar coupling constants have been observed for the first time in Watson-Crick AU base pairs in uniformly ¹⁵N-labeled RNA oligonucleotides using a new ^2hJ_NN-HNN-E. COSY experiment. The experiment utilizes adenosine H2 (AH2) for original polarization and detection, while employing ^2hJ_NNcouplings for coherence transfer across the hydrogen bonds (H-bonds). The H3 protons of uracil bases are unperturbed throughout the experiment so that these protons appear as passive spins in E. COSY patterns. ^3hJ_H2H3coupling constants can therefore be accurately measured in the acquisition dimension from the displacement of the E. COSY multiplet components, which are separated by the relatively large ¹J_H3N3coupling constants in the indirect dimension of the two-dimensional experiment. The ^3hJ_H2H3scalar coupling constants determined for AU base pairs in the two RNA hairpins examined here have been found to be positive and range in magnitude up to 1.8 Hz. Using a molecular fragment representation of an AU base pair, density functional theory/finite field perturbation theory (DFT/FPT) methods have been applied to attempt to predict the relative contributions of H-bond length and angular geometry to the magnitude of ^3hJ_H2H3coupling constants. Although the DFT/FPT calculations did not reproduce the full range of magnitude observed experimentally for the ^3hJ_H2H3coupling constants, the calculations do predict the correct sign and general trends in variation in size of these coupling constants. The calculations suggest that the magnitude of the coupling constants depends largely on H-bond length, but can also vary with differences in base pair geometry. The dependency of the ^3hJ_H2H3coupling constant on H-bond strength and geometry makes it a new probe for defining base pairs in NMR studies of nucleic acids.     

H···N hydrogen bond lengths in double stranded DNA from internucleotide dipolar couplings

The ratio of the internucleotide dipolar coupling and the corresponding one-bond imino ¹⁵N-¹H dipolar coupling provides a measure for the N···H/H-N distance ratio. Measurements were carried out for a dodecamer, d(CGCGAATTCGCG)₂, in which a C-G and an A-T basepair were uniformly enriched in ¹⁵N. When assuming H-bonds to be perfectly linear, dipolar data indicate time-averaged hydrogen bond lengths of 1.80±0.03 Å for A-T and 1.86±0.02 Å for C-G. When using H-bond orientations from high resolution X-ray data, H-bond lengths are about 0.1 Å shorter.     

Synthesis,biological evaluation,and molecular modelling studies of potent human neutrophil elastase (HNE) inhibitors

We report the synthesis and biological evaluation of a new series of 3- or 4-(substituted)phenylisoxazolones as HNE inhibitors. Due to tautomerism of the isoxazolone nucleus, two isomers were obtained as final compounds (2-NCO and 5-OCO) and the 2-NCO derivatives were the most potent with IC₅₀ values in the nanomolar range (20–70?nM). Kinetic experiments indicated that 2-NCO 7d and 5-OCO 8d are both competitive HNE inhibitors. Molecular modelling on 7d and 8d suggests for the latter a more crowded region about the site of the nucleophilic attack, which could explain its lowered activity. In addition molecular dynamics (MD) simulations showed that the isomer 8d appears more prone to form H-bond interactions which, however, keep the reactive sites quite distant for the attack by Ser195. By contrast the amide 7d appears more mobile within the active pocket, since it makes single H-bond interactions affording a favourable orientation for the nucleophilic attack.     

Hydrogen bonding interactions in noradrenaline-DMSO complexes: DFT and QTAIM studies of structure,properties and topology

The hydrogen bonding interactions between noradrenaline (NA) and DMSO were studied with density functional theory (DFT) regarding their geometries, energies, vibrational frequencies, and topological features of the electron density. The quantum theory of atoms in molecules (QTAIM) and the natural bond orbital (NBO) analyses were employed to elucidate the hydrogen bonding interaction characteristics in noradrenaline-DMSO complexes. The H-bonds involving the hydroxyls hydrogen in NA and the O atom in DMSO are dominant intermolecular H-bonds and are stronger than other H-bonds involving the methyl hydrogen of DMSO as a H-donor. The weak H-bonds also include a π H-bond which involves the benzene ring as a H-donor or H-acceptor. QTAIM identified the weak H-bonds formed between the methyl hydrogen of DMSO and the N atom in NA in some complexes (AB5, AB6 and AB7), which cannot be further confirmed by NBO and other methods, so there are probably no interactions between hydrogen and nitrogen atoms among these complexes. A good linear relationship between logarithmic electron density (lnρ _b ) at the bond critical point (BCP) and structural parameter (δR _H···Y) was found. The formations of new H-bonds in some complexes are helpful to strengthen the original intramolecular H-bond, this is attributed to the cooperativity of H-bonds in complexes and can be learned from the structure results and the NBO and QTAIM analyses. Analysis of various physically meaningful contributions arising from the energy decomposition procedures show that the orbital interactions of H-bond is predominant during the formation of the complex, moreover, both the hydrogen bonding interaction and the structural deformation are responsible for the stability of the complexes.     

Quantum-Chemical Analysis of C-H…O and C-H…N Interactions in RNA Base Pairs—H-Bond Versus Anti-H-Bond Pattern

Abstract

Geometries and interaction energies of unusual UU and AA base pairs with one standard hydrogen bond (H-bond) and additional C-H…O or C-H…N contacts have been determined by quantum-chemical methods taking into account electron correlation. Whereas the C-H bond length in the UU C-H…O contact increases upon complex formation (H-bond pattern), the C-H bond of the AA C-H….N interaction is shortened (anti-H-bond pattern). The same properties are found for model complexes between U or A and formaldehyde that have intermolecular C-H…acceptor contacts but no standard H-bonds. Both the C-H…acceptor H-bond and anti-H-bond interactions are attractive. A possible influence of the donor CH group charge distribution on the interaction pattern is discussed.     

Exploring the nature of the H-bonds between the human class II MHC protein,HLA-DR1 (DRB*0101) and the influenza virus hemagglutinin peptide,HA306-318, using the quantum theory of atoms in molecules

The nature of the H-bonds between the human protein HLA-DR1 (DRB*0101) and the hemagglutinin peptide HA306-318 has been studied using the Quantum Theory of Atoms in Molecules for the first time. We have found four H-bond groups: one conventional CO··HN bond group and three nonconventional CO··HC, π··HC involving aromatic rings and HN··HC_aliphatic groups. The calculated electron density at the determined H-bond critical points suggests the follow protein pocket binding trend: P1 (2,311) >> P9 (1.109) > P4 (0.950) > P6 (0.553) > P7 (0.213) which agrees and reveal the nature of experimental findings, showing that P1 produces by a long way the strongest binding of the HLA-DR1 human protein molecule with the peptide backbone as consequence of the vast number of H-bonds in the P1 area and at the same time the largest specific binding of the peptide Tyr308 residue with aromatic residues located at the binding groove floor. The present results suggest the topological analysis of the electronic density as a valuable tool that allows a non-arbitrary partition of the pockets binding energy via the calculated electron density at the determined critical points.     

Design,synthesis and docking study of novel coumarin ligands as potential selective acetylcholinesterase inhibitors

New coumaryl-thiazole derivatives with the acetamide moiety as a linker between the alkyl chains and/or the heterocycle nucleus were synthesized and in vitro tested as acetylcholinesterase (AChE) inhibitors. 2-(diethylamino)-N-(4-(2-oxo-2H-chromen-3-yl)thiazol-2-yl)acetamide (6c, IC₅₀ value of 43?nM) was the best AChE inhibitor with a selectivity index of 4151.16 over BuChE. Kinetic study of AChE inhibition revealed that 6c was a mixed-type inhibitor. Moreover, the result of H4IIE hepatoma cell toxicity assay for 6c showed negligible cell death. Molecular docking studies were also carried out to clarify the inhibition mode of the more active compounds. Best pose of compound 6c is positioned into the active site with the coumarin ring wedged between the residues of the CAS and catalytic triad of AChE. In addition, the coumarin ring is anchored into the gorge of the enzyme by H-bond with Tyr130.     

Identification of novel urease inhibitors: pharmacophore modeling,virtual screening and molecular docking studies

Abstract

Pharmacophore modeling and atom-based three-dimensional quantitative structure–activity relationship (3D-QSAR) have been developed on N-acylglycino- and hippurohydroxamic acid derivatives, which are known potential inhibitors of urease. This is followed by virtual screening and ADMET (absorption, distribution, metabolism, excretion and toxicity) studies on a large library of known drugs in order to get lead molecules as Helicobacter pylori urease inhibitors. A suitable three-featured pharmacophore model comprising one H-bond acceptor and two H-bond donor features (ADD.10) has been found to be the best QSAR model. An external library of compounds (～3000 molecules), pre-filtered using Lipinski’s rule of five, has been further screened using the pharmacophore model ADD.10. By analyzing the fitness of the hits with respect to the pharmacophore model and their binding interaction inside the urease active site, four molecules have been predicted to be extremely good urease inhibitors. Two of these have significant potential and should be taken up for further drug-designing process.     

Study of the phase transition in lysozyme crystals by Raman spectroscopy

BackgroundRecently, it has been revealed that tetragonal lysozyme crystals show a phase transition at 307 K upon heating. The underlying mechanisms of the phase transition are still not fully understood. Here we focus on the study of high-frequency vibrational modes arising from the protein and their temperature evolution in the vicinity of T_ph as well as on the detailed study of crystalline water dynamics near T_ph.MethodsRaman experiments have been performed at temperatures 295–323 K including T_ph. The low-frequency modes and the modes of fingerprint region, CH- and OH-stretching regions have been analyzed.Results and conclusionsIn spite of the absence of noticeable rearrangements in protein structure, the high-frequency vibrational modes of lysozyme located in the fingerprint region have been found to exhibit the features of critical dynamics near T_ph. Pronounced changes in the dynamics of α-helixes and Tyr residues exposed on the protein surface point to the important role of H-bond rearrangements at the phase transition. Additionally the study of temperature evolution of OH-stretching modes has shown an increase in distortions of tertahedral H-bond network of crystalline water above T_ph. These changes in water dynamics could play a crucial role in the mechanisms of the phase transition.General significanceThe present results shed light on the mechanisms of the phase transition in lysozyme crystals.     

Pyrazolo[3,4-d]pyrimidine Nucleic Acids: Adjustment of the dA-dT to the dG-dC Base Pair Stability

Abstract

The pyrazolo[3,4-d]pyrimidine-4,6-diamine nucleosides 2b-d stabilize the dA-dT base pair significantly when the dA-residue is replaced. Oligonucleotide duplexes incorporating 2b-d show a 4–6°C T _m increase per modification. The 7-bromo compound 2b harmonizes the stability of the dA-dT vs. the dG-dC pair. According to this the stability of such duplexes depends no longer on the base pair composition of a DNA molecule.     

Synthesis and Triplex Binding Properties of Oligonucleotides Containing a Novel Nucleobase

Abstract

The thiazolo-indole compound 1 bearing the complementary donor-acceptor-donor sites (dad) was designed for specific recognition of an AT inverted base pair in pyrimidine triple helix motif. It was successfully incorporated into 14-mer oligonucleotide using a serinol unit as sugar derivative. The triple helix hybridization studies were examined by means of thermal denaturation experiments with a 26-mer DNA duplex containing the AT inverted base pair.     

Structural Requirements at the 15.1–16.1 Position of the Hammerhead Ribozyme

Abstract

Apart from the A^15.1?U^16.1 -> I^15.1?C^16.1 change which reverses the polarity of an important H-bond in the hammerhead structure no other functional group changes at the 15.1 purine residue seem to be compatible with the requirements of efficient catalysis. The I^15.1 and A^15.1 ribozymes are equally suitable for practical applications because there are only minor differences in the acceptance of stabilising residues.     

Description of Base Motion and Role of Excitations in the Melting of Poly(dG) · Poly(dC)

Abstract

We study the contribution of various vibrational modes to the melting of poly(dG) · poly(dC). We find that the principal contribution comes from the H-bond breathing modes that have been observed in Raman scattering and that we have associated with helix melting. We show the softening of these modes on approach to melting in agreement with the observed behavior. We also describe the contribution to melting from base rotation modes that others have suggested are important in melting.     

Perturbation of Hydrogen Bonds in the Adenine…Thymine Base Pair by Na+, Mg2+, Ca2+ and NH4 + Cations

Abstract

Perturbation of the hydrogen bonds in the adenine…thymine base pair by Na⁺, Mg²⁺, Ca²⁺ and NH₄ ⁺ cations has been investigated by means of ab initio SCF calculations with the STO-3G basis set. The geometry of adenine…thymine, as well as those of the perturbed pairs were optimized. Approach of any cation to thymine at 06 leads to destabilization of the adenine…thy mine pair; divalent cations (Mg²⁺, Ca²⁺) have a profound effect on the structure of the base pair. The approach of a cation to other available sites (thymine: O2, adenine N1 and N3) leads, on the other hand, to stabilization of the base pair. If a water molecule is placed between the cation and the base pair, the structure and stability of the base pair are changed only negligibly.     

Observation of the closing of individual hydrogen bonds during TFE-induced helix formation in a peptide

Helix formation of an S-peptide analog, comprising the first 20 residues of Ribonuclease A and two additional N-terminal residues, was studied by measuring hydrogen bond (H-bond) (h3)J(NC') scalar couplings as a function of 2,2,2-trifluoroethanol (TFE) concentration. The (h3)J(NC') couplings give direct evidence for the closing of individual backbone N-H***O = C H-bonds during the TFE-induced formation of secondary structure. Whereas no (h3)J(NC') correlations could be detected without TFE, alpha-helical (i,i +4) H-bond correlations were observed for the amides of residues A5 to M15 in the presence of TFE. The analysis of individual coupling constants indicates that alpha-helix formation starts at the center of the S-peptide around residue E11 and proceeds gradually from there to both peptide ends as the TFE concentration is increased. At 60% to 90% TFE, well-formed alpha-helical H-bonds were observed for the amides hydrogens of residues K9 to Q13, whereas H-bonds of residues T5 to A8, H14, and M15 are affected by fraying. No intramolecular backbone H-bonds are present at and beyond the putative helix stop signal D16. As the (h3)J(NC') constants represent ensemble averages and the dependence of (h3)J(NC') on H-bond lengths is very steep, the size of the individual (h3)J(NC') coupling constants can be used as a measure for the population of a closed H-bond. These individual populations are in agreement with results derived from the Lifson-Roig theory for coil-to-helix transitions. The present work shows that the closing of individual H-bonds during TFE-induced helix formation can be monitored by changes in the size of H-bond scalar couplings.     

DNA Models for A,B, C and D Conformations Related to Fiber X-ray,Infrared and NMR Measurements

Abstract

A conformational analysis of the A, B, C and D DNA forms was made in order to establish molecular models presenting a good agreement with experimental data obtained from fiber X-ray, infrared linear dichroism and ³¹P NMR. The proposed models have been refined and do present good stereochemistry and optimized H-bond distances between bases associated with the Watson-Crick pairing. The DNA conformations proposed are a left handed double helix for the C form and right handed helices for A, B and D. Relations to conformational transitions between these forms are discussed.     

Comprehensive determination of <Superscript>3</Superscript>J<Subscript>HNHα</Subscript> for unfolded proteins using <Superscript>13</Superscript>C′-resolved spin-echo difference spectroscopy

An experiment is presented to determine ³J_HNHα coupling constants, with significant advantages for applications to unfolded proteins. The determination of coupling constants for the peptide chain using 1D ¹H, or 2D and 3D ¹H-¹⁵N correlation spectroscopy is often hampered by extensive resonance overlap when dealing with flexible, disordered proteins. In the experiment detailed here, the overlap problem is largely circumvented by recording ¹H-¹³C′ correlation spectra, which demonstrate superior resolution for unfolded proteins. J-coupling constants are extracted from the peak intensities in a pair of 2D spin-echo difference experiments, affording rapid acquisition of the coupling data. In an application to the cytoplasmic domain of human neuroligin-3 (hNlg3cyt) data were obtained for 78 residues, compared to 54 coupling constants obtained from a 3D HNHA experiment. The coupling constants suggest that hNlg3cyt is intrinsically disordered, with little propensity for structure.     

Tyrosine deprotonation and associated hydrogen bond rearrangements in a photosynthetic reaction center

Photosynthetic reaction centers from Blastochloris viridis possess Tyr-L162 located mid-way between the special pair chlorophyll (P) and the heme (heme3). While mutation of the tyrosine does not affect the kinetics of electron transfer from heme3 to P, recent time-resolved Laue diffraction studies reported displacement of Tyr-L162 in response to the formation of the photo-oxidized P^+•, implying a possible tyrosine deprotonation event. pK _a values for Tyr-L162 were calculated using the corresponding crystal structures. Movement of deprotonated Tyr-L162 toward Thr-M185 was observed in P^+• formation. It was associated with rearrangement of the H-bond network that proceeds to P via Thr-M185 and His-L168.