On tautomerism of the cytosine molecule

Tautomerism of the cytosine molecule is discussed in connection with recent experimental matrix-isolation infrared spectroscopic measurements and recent ab initio calculations of relative stabilities of tautomers and of IR spectra for different tautomeric forms of the compound. Experimental IR spectra in the N-H and O-H stretching regions and in the C = O stretching region are presented for cytosine and for its several derivatives considered as model compounds. This experimental evidence, as well as the quantum-mechanical calculations (including both electron correlation and zero-point vibrational contributions), clearly indicate that two tautomers of cytosine, i.e. the amino-hydroxy and amino-oxo forms with the hydrogen atom at the N(1) position, exist in equilibrium when the cytosine molecule is isolated in an inert environment. The effect of the environment on the relative stabilities of several tautomers is also discussed briefly.     

Ab initio prediction of optical rotation: comparison of density functional theory and Hartree-Fock methods for three 2,7,8-trioxabicyclo[3.2.1]octanes

We report ab initio calculations of the frequency-dependent electric dipole-magnetic dipole polarizabilities, beta(nu), at the sodium D line frequency and, thence, of the specific rotations, [alpha](D), of 2,7,8-trioxabicyclo[3.2.1]octane, 1, and its 1-methyl derivative, 2, using the Density Functional Theory (DFT) and Hartree-Fock/Self-Consistent Field (HF/SCF) methodologies. Gauge-invariant (including) atomic orbitals (GIAOs) are used to ensure origin-independent [alpha](D) values. Using large basis sets which include diffuse functions DFT [alpha](D) values are in good agreement with experimental values (175.8 degrees and 139.2 degrees for (1S,5R)-1 and -2, respectively); errors are in the range 25-35 degrees. HF/SCF [alpha](D) values, in contrast, are much less accurate; errors are in the range 75-95 degrees. The use of small basis sets which do not include diffuse functions substantially lowers the accuracy of predicted [alpha](D) values, as does the use of the static limit approximation: beta(nu) approximately beta(o). The use of magnetic-field-independent atomic orbitals, FIAOs, instead of GIAOs, leads to origin-dependent, and therefore nonphysical, [alpha](D) values. We also report DFT calculations of [alpha](D) for the 1-phenyl derivative of 1, 3. DFT calculations find two stable conformations, differing in the orientation of the phenyl group, of very similar energy, and separated by low barriers. Values of [alpha](D) predicted using two different algorithms for averaging over phenyl group orientations are in good agreement with experiment. In principle, the absolute configuration (AC) of a chiral molecule can be assigned by comparison of the optical rotation predicted ab initio to the experimental value. Our results demonstrate the critical importance of the choice of ab initio methodology in obtaining reliable optical rotations and, hence, ACs, and show that, at the present time, DFT constitutes the method of choice.     

Active site contacts in the purine nucleoside phosphorylase--hypoxanthine complex by NMR and ab initio calculations

Hypoxanthine (Hx) with specific (15)N labels has been used to probe hydrogen-bonding interactions with purine nucleoside phosphorylase (PNP) by NMR spectroscopy. Hx binds to human PNP as the N-7H tautomer, and the N-7H (1)H and (15)N chemical shifts are located at 13.9 and 156.5 ppm, respectively, similar to the solution values. In contrast, the (1)H and (15)N chemical shifts of N-1H in the PNP.Hx complex are shifted downfield by 3.5 and 7.5 ppm to 15.9 and 178.8 ppm, respectively, upon binding. Thus, hydrogen bonding at N-1H is stronger than at N-7H in the complex. Ab initio chemical shift calculations on model systems that simulate Hx in solution and bound to PNP are used to interpret the NMR data. The experimental N-7H chemical shift changes are caused by competing effects of two active site contacts. Hydrogen bonding of Glu201 to N-1H causes upfield shifts of the N-7H group, while the local hydrogen bond (C=O to N-7H from Asn243) causes downfield shifts. The observed N-7H chemical shift can be reproduced by a hydrogen bond distance approximately 0.13 A shorter (but within experimental error) of the experimental value found in the X-ray crystal structure of the bovine PNP.Hx complex. The combined use of NMR and ab initio chemical shift computational analysis provides a novel approach to understand enzyme-ligand interactions in PNP, a target for anticancer agents. This approach has the potential to become a high-resolution tool for structural determination.     

Importance of oligoelectrolyte end effects for the thermodynamics of conformational transitions of nucleic acid oligomers: a grand canonical Monte Carlo analysis.

Effects of salt concentration on the stabilities of oligonucleotide helices are analyzed directly in terms of delta gamma N----yN identical to gamma denyN - gamma natN, the difference in the salt-nucleotide phosphate preferential interaction coefficients for the denatured state, having yN phosphate charges, and for the native state, having N phosphate charges (y = 1 for hairpin denaturation and y = 0.5 for dimer denaturation). Previous experimental studies of the denaturation of hairpin oligonucleotides (having 18 less than N less than 44) indicate significant differences between delta gamma N----N and delta gamma infinity, the value determined for the denaturation of the corresponding polynucleotide. These differences are thermodynamic manifestations of the oligoelectrolyte end effect. In contrast, the available data on the denaturation of oligonucleotide dimer helices (N less than or equal to 22) imply that differences between delta gamma infinity and delta gamma N----0.5N, and hence oligoelectrolyte end effects, are small or negligible. To determine the origin of these apparently conflicting implications concerning the importance of oligoelectrolyte end effects, we have calculated the N dependence of gamma N from grand canonical Monte Carlo simulations for an idealized model of the structure and charge distribution of each oligomer conformation. Our calculations are in quantitative agreement with the experimental finding for d(TA) hairpin oligomers that -delta gamma N----N decreases linearly as N-1 increases, and with the extant experimental determinations of delta gamma N----0.5N. These results provide an illustration of how the large electrostatic end effects exhibited by the hairpin denaturation data are masked when delta gamma infinity is compared with values of delta gamma N----0.5N for short dimer helices (N less than or equal to 22). For 0.5N greater than 24, -delta gamma N----0.5N is predicted to be a linear function of N-1 whose slope has the opposite sign from, and is more salt-concentration dependent than, the corresponding slope of -delta gamma N----N as a function of N-1. Our calculations also yield predictions about the N dependences of the individual values of gamma N that can be tested by determining Donnan coefficients from membrane dialysis equilibrium experiments. For long enough hairpin and dimer oligonucleotides (yN greater than or equal to 24), in either native or denatured forms, we predict that the (positive) difference gamma infinity - gamma N increases linearly with increasing N-1. For smaller values of N the difference gamma infinity - gamma N continues to increase with increasing N-1.     

Low-barrier hydrogen bond hypothesis in the catalytic triad residue of serine proteases: correlation between structural rearrangement and chemical shifts in the acylation process

To elucidate the catalytic advantage of the low-barrier hydrogen bond (LBHB), we analyze the hydrogen bonding network of the catalytic triad (His57-Asp102-Ser195) of serine protease trypsin, one of the best examples of the LBHB reaction mechanism. Especially, we focus on the correlation between the change of the chemical shifts and the structural rearrangement of the active site in the acylation process. To clarify LBHB, we evaluate the two complementary properties. First, we calculate the NMR chemical shifts of the imidazole ring of His57 by the gauge-including atomic orbital (GIAO) approach within the ab initio QM/MM framework. Second, the free energy profile of the proton transfer from His57 to Asp102 in the tetrahedral intermediate is obtained by ab initio QM/MM calculations combined with molecular dynamics free energy perturbation (MD-FEP) simulations. The present analyses reveal that the calculated shifts reasonably reproduce the observed values for (1)H chemical shift of H(epsilon)(1) and H(delta)(1) in His57. The (15)N and (13)C chemical shifts are also consistent with the experiments. It is also shown that the proton between His57 and Asp102 is localized at the His57 side. This largely downfield chemical shift is originated from the strong electrostatic interaction, not a covalent-like bonding character between His57 and Asp102. Also, it is proved that a slight downfield character of H(epsilon)(1) is originated from a electrostatic interaction between His57 and the backbone carbonyl group of Val213 and Ser214. These downfield chemical shifts are observed only when the tetrahedral intermediate is formed in the acylation process.     

An Investigation of Structural Stability and Internal Rotation in 3-Cyclopropenecarboxaldehyde and 3-Cyclopropenecarboxylic Acid Fluoride by ab initio Calculations

The structural stability and internal rotation in 3-cyclopropenecarboxaldehyde and 3-cyclo-propenecarboxylic acid fluoride were investigated by ab initio calculations with a 6-31G* atomic basis in the latter and a 6-311G* atomic basis set in the former case. For the sake of comparison also results obtained with a 3-21G basis are given in the paper. As expected, it turned out that this basis set is not large enough for three-membered rings. The calculations were carried out both at the Restricted Hartree-Fock (HF) and the second order Moller-Plesset (MP2) levels. The trans-form is predicted to be the lower energy conformer for both molecules. However, in case of the fluoride the two conformers have nearly the same energy. Full optimization was performed at the transition states and the asymmetric potential function for the CXO internal rotations was predicted for both molecules.     

A modified version of the Cornell et al. force field with improved sugar pucker phases and helical repeat

We have examined some subtle parameter modifications to the Cornell et al. force field, which has proven quite successful in reproducing nucleic acid properties, but whose C2'-endo sugar pucker phase and helical repeat for B DNA appear to be somewhat underestimated. Encouragingly, the addition of a single V2 term involving the atoms C(sp3)-O-(sp3)-C(sp3)-N(sp2), which can be nicely rationalized because of the anomeric effect (lone pairs on oxygen are preferentially oriented relative to the electron withdrawing N), brings the sugar pucker phase of C2'-endo sugars to near perfect agreement with ab initio calculations (W near 162 degrees). Secondly, the use of high level ab initio calculations on entire nucleosides (in contrast to smaller model systems necessitated in 1994-95 by computer limitations) lets one improve the chi torsional potential for nucleic acids. Finally, the O(sp3)-C(sp3)- C(sp3)-O(sp3) V2 torsional potential has been empirically adjusted to reproduce the ab initio calculated relative energy of C2'-endo and C3'-endo nucleosides. These modifications are tested in molecular dynamics simulations of mononucleosides (to assess sugar pucker percentages) and double helices of DNA and RNA (to assess helical and sequence specific structural properties). In both areas, the modified force field leads to improved agreement with experimental data.     

High-resolution structure determination of the CylR2 homodimer using paramagnetic relaxation enhancement and structure-based prediction of molecular alignment

Structure determination of homooligomeric proteins by NMR spectroscopy is difficult due to the lack of chemical shift perturbation data, which is very effective in restricting the binding interface in heterooligomeric systems, and the difficulty of obtaining a sufficient number of intermonomer distance restraints. Here we solved the high-resolution solution structure of the 15.4 kDa homodimer CylR2, the regulator of cytolysin production from Enterococcus faecalis, which deviates by 1.1 angstroms from the previously determined X-ray structure. We studied the influence of different experimental information such as long-range distances derived from paramagnetic relaxation enhancement, residual dipolar couplings, symmetry restraints and intermonomer Nuclear Overhauser Effect restraints on the accuracy of the derived structure. In addition, we show that it is useful to combine experimental information with methods of ab initio docking when the available experimental data are not sufficient to obtain convergence to the correct homodimeric structure. In particular, intermonomer distances may not be required when residual dipolar couplings are compared to values predicted on the basis of the charge distribution and the shape of ab initio docking solutions.     

Coordination properties of 2-aminocyclopentene-1-dithiocarboxylic acid to transition metal ions as studied by ab initio calculations

The question of the(N, S) vs. (S, S) coordination mode on M x (ACDA)(2) complexes (ACDA=2-aminocyclopentene-1-dithiocarboxylic acid, M=Ni(2+), Pd(2+), Pt(2+)) was assessed through an extensive ab initio study, using the hybrid B3LYP density functional approach. The (S,S)coordination was found to be the most stable one, with an energy difference of ca. 50 kJ mol(-1) relative to the(N, S) coordination mode. Detailed analysis of the ab initio results indicates that this preference is a result of the combined effect of geometry constraints and electron distribution within the complex.     

Discriminative learning for protein conformation sampling

Protein structure prediction without using templates (i.e., ab initio folding) is one of the most challenging problems in structural biology. In particular, conformation sampling poses as a major bottleneck of ab initio folding. This article presents CRFSampler, an extensible protein conformation sampler, built on a probabilistic graphical model Conditional Random Fields (CRFs). Using a discriminative learning method, CRFSampler can automatically learn more than ten thousand parameters quantifying the relationship among primary sequence, secondary structure, and (pseudo) backbone angles. Using only compactness and self-avoiding constraints, CRFSampler can efficiently generate protein-like conformations from primary sequence and predicted secondary structure. CRFSampler is also very flexible in that a variety of model topologies and feature sets can be defined to model the sequence-structure relationship without worrying about parameter estimation. Our experimental results demonstrate that using a simple set of features, CRFSampler can generate decoys with much higher quality than the most recent HMM model.     

Ab initio computational study of beta-cellobiose conformers using B3LYP/6-311++G**

The molecular structure of 27 conformers of beta-cellobiose were studied in vacuo through gradient geometry optimization using B3LYP density functionals and the 6-311++G** basis set. The conformationally dependent geometry changes and energies were explored as well as the hydrogen-bonding network. The lowest electronic energy structures found were not those suggested from available crystallographic and NMR solution data, where the glycosidic dihedral angles fall in the region (phi, psi) approximately (40 degrees, -20 degrees ). Rather, 'flipped' conformations in which the dihedral angles are in the range (phi, psi) approximately (180 degrees, 0 degrees ) are energetically more stable by approximately 2.5 kcal/mol over the 'experimentally accepted' structure. Further, when the vibrational free energy, deltaG, obtained from the calculated frequencies, is compared throughout the series, structures with (phi, psi) in the experimentally observed range still have higher free energy ( approximately 2.0 kcal/mol) than 'flipped' forms. The range of bridging dihedral angles of the 'normal' conformers, resulting from the variance in the phi dihedral is larger than that found in the 'flipped' forms. Due to this large flat energy surface for the normal conformations, we surmise that the summation of populations of these conformations will favor the 'normal' conformations, although evidence suggests that polar solvent effects may play the dominant role in providing stability for the 'normal' forms. Even though some empirical studies previously found the 'flipped' conformations to be lowest in energy, these studies have been generally discredited because they were in disagreement with experimental results. Most of the DFT/ab initio conformations reported here have not been reported previously in the ab initio literature, in part because the use of less rigorous theoretical methods, i.e. smaller basis sets, have given results in general agreement with experimental data, that is, they energetically favored the 'normal' forms. These are the first DFT/ab initio calculations at this level of theory, apparently because of the length and difficulty of carrying out optimizations at these high levels.     

A novel approach for assessing macromolecular complexes combining soft-docking calculations with NMR data

We present a novel and efficient approach for assessing protein-protein complex formation, which combines ab initio docking calculations performed with the protein docking algorithm BiGGER and chemical shift perturbation data collected with heteronuclear single quantum coherence (HSQC) or TROSY nuclear magnetic resonance (NMR) spectroscopy. This method, termed "restrained soft-docking," is validated for several known protein complexes. These data demonstrate that restrained soft-docking extends the size limitations of NMR spectroscopy and provides an alternative method for investigating macromolecular protein complexes that requires less experimental time, effort, and resources. The potential utility of this novel NMR and simulated docking approach in current structural genomic initiatives is discussed.     

Monensin A acid complexes as a model of electrogenic transport of sodium cation

New Monensin A acid complexes with water molecule, sodium chloride and sodium perchlorate were obtained and studied by X-ray and (1)H, (13)C NMR and FT-IR methods as well as ab initio calculations. The crystal structure of the complexes indicates the complexation of the water molecule and Na(+) cation in the pseudo-cycle conformation of the Monensin acid molecule stabilised by intramolecular hydrogen bonds. Important for stabilisation of this structure is also the intermolecular hydrogen bonds with water molecule or the coordination bonds with Na(+) cation. It is demonstrated that the counterions forming intermolecular hydrogen bonds with OH groups influence the strength of the intramolecular hydrogen bonds, but they have no influence on the formation of pseudo-cyclic structure. Spectroscopic studies of the complexes in dichloromethane solution have shown that the pseudo-cyclic structure of the compounds is conserved. As follows from the ab initio calculations, the interactions between the Na(+) cation and the electronegative oxygen atoms of Monensin acid totally change the molecular electrostatic potential around the supramolecular Monensin acid-Na(+) cationic complex relative to that of the neutral Monensin acid molecule.     

Assigning stereodiversity of the 27-Me group of furostane-type steroidal saponins via NMR chemical shifts

Applicability of (13)C and (1)H NMR chemical shifts for the assignment of the 25R/25S configuration of the 27-methyl group in the case of furostane-type steroidal saponins has been investigated. A comparative study of (13)C NMR data suggest that chemical shift values for C-20, C-21, C-22, C-23, C-24, C-25, C-26 and C-27 resonances were not much influenced by R/S configuration of the 27-Me group, thus reflecting limited application of (13)C NMR chemical shifts for such stereochemical determinations. In contrast, (1)H NMR chemical shifts (delta(a), delta(b)) for geminal protons of glycosyloxy methylene (H(2)-26) exhibit pronounced dependence and the difference (Delta(ab)=delta(a)-delta(b)) among their chemical shifts [Delta(ab)= or <0.48 for 25R; Delta(ab)= or >0.57 for 25S] seems to be of general applicability for ascertaining 25R/25S orientation of the 27-methyl group of furostane-type steroidal saponins.     

Density-functional theory and ab initio Hartree-Fork studies on the structural parameters and chemical activity of the free radicals generated by benzoquinone and hydroquinone

Density-functional theory (DFT) calculations were performed for calculation of the theoretical spectra and the chemical activities of free radicals generated by benzoquinone and hydroquinone as well as the transition states, and the calculated spectra were used for the assignment of the frequencies observed in the experimental IR spectra. The calculated geometrical parameters, the predicted IR spectra, and the chemical activities of free radicals and transition states were also compared with those of benzoquinone and hydroquinone. The reactive mechanisms of free radicals generated by benzoquinone and hydroquinone are also discussed using ab initio Hartree-Fork (HF) methods.     

Conformational Equilibria of Bridled Chiroporphyrins in Solution Investigated by Vibrational Circular Dichroism

A series of bridled chiroporphyrins (BCP) and their metal complexes were prepared, in which two n‐methylene straps connect adjacent meso substituents by ester linkages. These compounds can exist as four atropisomers (αααα, αβαβ, αααβ, or ααββ) depending on the position of the meso groups relative to the macrocycle (α when above and β when below). We characterized the conformation of these chiral porphyrins and their metal (Zn, Ni, Mn) complexes by vibrational circular dichroism (VCD) associated with ab initio calculations. VCD spectra of the three metalloporphyrins were recorded in CDCl₃ and benzene solutions and ab initio calculations of their four atropoisomers were performed at the Density Functional Theory (DFT) level. The bridled chiroporphyrin with the longer straps (9 CH₂) and its nickel(II) complex can be isolated as the αβαβ atropisomer in the solid state and were found with the same conformation in CDCl₃ and benzene solutions. The bridled chiroporphyrin with the shortest straps (8 CH₂) and its zinc(II) complex can be isolated as the αααα atropisomer in the solid state, but in solution they are subject to atropisomeric equilibria, resulting in atropisomer distributions that are strongly solvent‐dependent. Comparison of the experimental VCD spectra with the predicted spectra of the four atropisomers allowed the quantification of these distributions. Finally, the manganese(III) complex also exhibits an atropisomeric equilibria in solution which is slightly solvent‐dependent. Chirality 25:480–486, 2013. © 2013 Wiley Periodicals, Inc.     

Ab initio conformational maps for disaccharides in gas phase and aqueous solution

Ab initio conformational maps for beta-lactose in both the gas phase and in aqueous solution have been constructed at the HF/6-31G(d,p) level of calculation. The results of the gas-phase ab initio calculations allow us to conclude that a rigid conformational map is able to predict the regions of the minima in the potential energy surface of beta-lactose, in full agreement with those found in the relaxed conformational map. The solvation effects do not give rise to any new local minimum in the potential energy surface of beta-lactose, but just change the relative Boltzmann populations of the conformers found in the gas-phase calculations. The values obtained for heteronuclear spin coupling constant (3J(H,C)), using the seven most stable conformers in solution are in good agreement with the available experimental values. This is a good indication that ab initio rigid conformational maps can be reliably used to sort the most stable conformers of beta-lactose.