Quantum-mechanical studies on the conformation of nucleic acids and their constituents

Structural insight into antisense gapmer-RNA oligomer duplexes through molecular dynamics simulations

There is an extensive research carrying out on antisense technology and the molecules entering into clinical trials are increasing rapidly. Phosphorothioate (PS) is a chemical modification in which nonbridged oxygen is replaced with a sulfur, consequently providing resistance against nuclease activity. The 2'-4' conformationally restricted nucleoside has the structural features of both 2'-O-methoxy ethyl RNA (MOE), which shows good toxicity profile, and locked nucleic acid (LNA), which shows good binding affinity towards the target RNA. These modifications have been studied and suggested that they can be a potential therapeutic agents in antisense therapy. Mipomersen (ISIS 301012), which contains the novel nucleoside modification has been used to target to apolipoprotein (Apo B), which reduces LDL cholesterol by 6-41%. In this study, classical molecular dynamics (MD) simulations were performed on six different antisense gapmer/target-RNA oligomer duplexes (LNA-PS-LNA/RNA, RcMOE-PS-RcMOE/RNA, ScMOE-PS-ScMOE/RNA, MOE-PS-MOE/RNA, PS-DNA/RNA and DNA/RNA) to investigate the structural dynamics, stability and solvation properties. The LNA, MOE nucleotides present in respective duplexes are showing the structure of A-form and the PS-DNA nucleotides resemble the structure of B-form helix with respect to some of the helical parameters. Free energy calculations suggest that the oligomer, which contains LNA binds to the RNA strongly than other modifications as shown in experimental results. The MOE modified nucleotide, which although had a lower binding affinity but higher solvent accessible surface area (SASA) compared to the other modifications, may be influencing the toxicity and hence may be used it in Mipomersen, the second antisense molecule which is approved by FDA. Communicated by Ramaswamy H. Sarma.

A mechanochemical switch to control radical intermediates

B₁₂-dependent enzymes employ radical species with exceptional prowess to catalyze some of the most chemically challenging, thermodynamically unfavorable reactions. However, dealing with highly reactive intermediates is an extremely demanding task, requiring sophisticated control strategies to prevent unwanted side reactions. Using hybrid quantum mechanical/molecular mechanical simulations, we follow the full catalytic cycle of an AdoB₁₂-dependent enzyme and present the details of a mechanism that utilizes a highly effective mechanochemical switch. When the switch is "off", the 5'-deoxyadenosyl radical moiety is stabilized by releasing the internal strain of an enzyme-imposed conformation. Turning the switch "on," the enzyme environment becomes the driving force to impose a distinct conformation of the 5'-deoxyadenosyl radical to avoid deleterious radical transfer. This mechanochemical switch illustrates the elaborate way in which enzymes attain selectivity of extremely chemically challenging reactions.

Structural significance of hypermodified nucleic acid base hydroxywybutine (OHyW) which occur at 37th position in the anticodon loop of yeast tRNA(Phe)

Conformational preferences of hypermodified nucleic acid base hydroxywybutine (OHyW) have been studied using quantum chemical single point semi-empirical PM3 method. Automated geometry optimization using semi-empirical RM1, molecular mechanics force field (MMFF) along with ab-initio HF-SCF (6-31G** basis set) and DFT (B3LYP/6-31G** basis set) calculations have also been made to compare the salient features. Molecular electrostatic potentials (MEPs) depict the polarities of hydroxywybutine (OHyW) side chain. Another conformational study showed that hydroxywybutosine side chain interacts with adjacent bases within the anticodon loop of tRNA(Phe). The solvent accessible surface area (SASA) calculations revealed the structural role of hydroxywybutine in anticodon loop. Explicit molecular dynamics (MD) simulation has been done over the PM3 most stable structure of OHyW. The hydroxywybutine side chain prefers 'distal' conformation i.e. spreads away from the cyclic five membered imidazole moiety of modified tricyclic guanine base. The predicted preferred conformation of hydroxywybutine may prevent extended Watson-Crick base pairing during protein biosynthesis process. This conformation of OHyW stabilized by intramolecular interactions between O(6)⋯HO(16), O(6)⋯HC(15) and O(20)⋯HC(17). Further stabilization is also expected from interactions between O(22)⋯HC(16) and O(23)⋯HC(15). Explicit molecular dynamics (MD) simulation over the PM3 most stable structure of OHyW support the preferred geometry by preserving the 'distal' orientation of hydroxywybutine side chain and intramolecular hydrogen bonding interactions. MD simulation study revealed the role of hydroxyl group of OHyW to avoid fluctuations and prevent multiple iso-energetic conformations of hydroxywybutine side chain as compared to wybutine (yW).

Conformational preferences of modified nucleoside N(2)-methylguanosine (m(2)G) and its derivative N(2), N(2)-dimethylguanosine (m(2)(2)G) occur at 26th position (hinge region) in tRNA

Conformational preferences of the modified nucleosides N(2)-methylguanosine (m(2)G) and N(2), N(2)-dimethylguanosine (m(2)(2)G) have been studied theoretically by using quantum chemical perturbative configuration interaction with localized orbitals (PCILO) method. Automated complete geometry optimization using semiempirical quantum chemical RM1, along with ab initio molecular orbital Hartree-Fock (HF-SCF), and density functional theory (DFT) calculations has also been made to compare the salient features. Single-point energy calculation studies have been made on various models of m(2)G26:C/A/U44 and m(2)(2)G26:C/A/U44. The glycosyl torsion angle prefers "syn" (χ = 286°) conformation for m(2)G and m(2)(2)G molecules. These conformations are stabilized by N(3)-HC2' and N(3)-HC3' by replacing weak interaction between O5'-HC(8). The N(2)-methyl substituent of (m(2)G26) prefers "proximal" or s-trans conformation. It may also prefer "distal" or s-cis conformation that allows base pairing with A/U44 instead of C at the hinge region. Thus, N(2)-methyl group of m(2)G may have energetically two stable s-trans m(2)G:C/A/U or s-cis m(2)G:A/U rotamers. This could be because of free rotations around C-N bond. Similarly, N(2), N(2)-dimethyl substituent of (m(2)(2)G) prefers "distal" conformation that may allow base pairing with A/U instead of C at 44th position. Such orientations of m(2)G and m(2)(2)G could play an important role in base-stacking interactions at the hinge region of tRNA during protein biosynthesis process.

Computational insights into the mechanism of radical generation in B12-dependent methylmalonyl-CoA mutase

ONIOM calculations have provided novel insights into the mechanism of homolytic Co-C5' bond cleavage in the 5'-deoxyadenosylcobalamin cofactor catalyzed by methylmalonyl-CoA mutase. We have shown that it is a stepwise process in which conformational changes in the 5'-deoxyadenosine moiety precede the actual homolysis step. In the transition state structure for homolysis, the Co-C5' bond elongates by approximately 0.5 Angstroms from the value found in the substrate-bound reactant complex. The overall barrier to homolysis is approximately 10 kcal/mol, and the radical products are approximately 2.5 kcal/mol less stable than the initial ternary complex of enzyme, substrate, and cofactor. The movement of the deoxyadenosine moiety during the homolysis step positions the resulting 5'-deoxyadenosyl radical for the subsequent hydrogen atom transfer from the substrate, methylmalonyl-CoA.

Conformational analysis of canonical 2-deoxyribonucleotides. 2. Purine nucleotides

The molecular structure of different conformers of isolated canonical purine 2'-deoxyribonucleotides 2-deoxyadenosine-5'-phosphate (pdA) and 2'-deoxyguanosine-5'-phosphate (pdG) was optimized using the B3LYP/6-31G(d) method. The results of the calculations reveal that the geometrical parameters and relative stability of the conformers significantly depend on the nature of the nucleobase, its orientation, the conformation of the furanose ring, the charge of the phosphate group and the character of the intramolecular hydrogen bonds. Analysis of the electron density distribution in purine nucleotides reveals the existence of a number of intramolecular hydrogen bonds. In general, the south conformer has a lower energy at the anti orientation of the base, and both conformers occur as the most stable for the syn orientation of the nucleobases. The results of the calculations reveal that the geometry and relative energy of the conformers of purine DNTs may be easily tuned by the charge of the phosphate group.

Internal degrees of freedom, structural motifs, and conformational energetics of the 5'-deoxyadenosyl radical: implications for function in adenosylcobalamin-dependent enzymes. A computational study

The potential energy surface of the free 5'-deoxyadenosyl radical in the gas phase is explored using density functional and second-order Møller-Plesset perturbation theories with 6-31G(d) and 6-31++G(d,p) basis sets and interpreted in terms of attractive and repulsive interactions. The 5',8-cyclization is found to be exothermic by approximately 20 kcal/mol but kinetically unfavorable; the lowest cyclization transition state (TS) lies about 7 kcal/mol higher than the highest TS for conversion between most of the open isomers. In open isomers, the two energetically most important attractive interactions are the hydrogen bonds (a) between the 2'-OH group and the N3 adenine center and (b) between the 2'-OH and 3'-OH groups. The relative ribose-adenine rotation about the C1'-N9 glycosyl bond in a certain range changes the energy by as much as 10-15 kcal/mol, the origin being (i) the repulsive 2'-H.H-C8 and O1'.N3 and (ii) the attractive 2'-OH.N3 ribose-adenine interactions. The hypothetical synergy between the glycosyl rotation and the Co-C bond scission may contribute to the experimentally established labilization of the Co-C bond in enzyme-bound adenosylcobalamin. The computational results are not inconsistent with the rotation about the C1'-N9 glycosyl bond being the principal coordinate for long-range radical migration in coenzyme B(12)-dependent enzymes. The effect of the protein environment on the model system results reported here remains an open question.

A pentapeptide model for an early folding step in the refolding of staphylococcal nuclease: the role of its turn propensity

Recently the folding of a staphylococcal nuclease (P117G) variant was examined with the hydrogen-deuterium (H-D) exchange technique. Many of the residues that showed significant protection are located in beta-sheet regions. About half the residues protected belong to an antiparallel beta-hairpin structure (residues 21-35) in the native structure. The beta-hairpin structure is formed by strands 2 and 3 of sheet 2 connected by the sequence 27Y KGQP31 in a type 1' reverse turn conformation with a 4-->1 hydrogen bonding between Q30 NH and Y27 C = O. We have targeted the conformational characterization of the peptide model Ac-YKGQP-NH2 with 1H two-dimensional nmr techniques in aqueous solution with a view to assessing its propensity to sample turn conformational forms and thus initiate the formation of beta-hairpin structure. Based upon the observed d alpha n (i, i + 1), d alpha n (i, i + 3), and dnn (i, i + 1) nuclear Overhauser effect connectivities, temperature coefficients for amide protons and conformational analysis with quantum mechanical perturbative configuration interaction over localized orbitals method, we conclude that the model peptide samples turn conformational forms with reduced conformational entropy. We suggest that the turn can nucleate the formation of the beta-hairpin structure in the refolding of nuclease. Observation of turn propensity for this sequence is consistent with the folding mechanism of the Greek key motif (present in Staphylococcal nuclease) proposed in the literature.

A study into the effects of protein binding on nucleotide conformation

In this study, we examine the effects of binding to protein upon nucleotide conformation, by the comparison of X-ray crystal structures of free and protein-bound nucleotides. A dataset of structurally non-homologous protein-nucleotide complexes was derived from the Brookhaven Protein Data Bank by a novel protocol of dual sequential and structural alignments, and a dataset of native nucleotide structures was obtained from the Cambridge Structural Database. The nucleotide torsion angles and sugar puckers, which describe nucleotide conformation, were analysed in both datasets and compared. Differences between them are described and discussed. Overall, the nucleotides were found to bind in low energy conformations, not significantly different from their 'free' conformations except that they adopted an extended conformation in preference to the 'closed' structure predominantly observed by free nucleotide. The archetypal conformation of a protein-bound nucleotide is derived from these observations.

UV photoelectron spectroscopy and ab initio characterization of valence orbital structures and conformations of neutral phosphate esters

The HeI UV photoelectron spectrum of trimethyl phosphate (TMP) has been measured and interpreted with the aid of SCF molecular orbital calculations carried out with STO-3G, STO-3G* and 4-31G basis functions. The photoelectron spectrum of TMP is more accurately reproduced by results from 4-31G calculations than by results from STO-3G or STO-3G* calculations. However, all three basis sets yield results which predict the same assignment of the photoelectron spectrum. Results at the 4-31G level indicate that whether calculations are based on crystallographic bond angles and bond lengths or on STO-3G optimized geometries has little effect on the energetic ordering of the upper occupied orbitals. The energetic ordering of orbitals is also found to be only weakly dependent upon the torsional angle phi, describing rotation of ester groups about P-O bonds and upon the torsional angle psi, describing rotation of methyl groups about C-O bonds. For trimethyl phosphate, with C3 symmetry, the vertical ionization potentials of the upper occupied orbitals are 10.81 eV (8e), 11.4 eV (9a), 11.93 eV (7e), 12.6-12.9 eV (8a and 6e), 14.4 eV (7a) and 15.0-16.0 eV (5e and 6a). Calculations at the 4-31G level indicate that many of the highest occupied orbitals in neutral dimethyl phosphate and methyl phosphate have energies and electron distributions similar to orbitals in TMP. For TMP, a search for optimized values of phi and psi has been carried out at the STO-3G*level. In agreement with previous NMR studies and with classical potential calculations, the STO-3G* results indicate that both the gauche (phi = 53.1 degrees) and anticlinal (phi = 141.9 degrees) conformations are thermally accessible. Also in agreement with the classical potential calculations, the STO-3G* results predict that in the all gauche conformation energy is minimized when the methyl groups assume a staggered geometry (psi = 60 degrees to 80 degrees) and that an energy maximum occurs for an eclipsed geometry (phi = 0 degrees to 20 degrees). A study of the dependence of optimized values of O-P-O ester bond angles on the torsional angles, phi, was carried out at the STO-3G, STO-3G* and 4-31G levels. The results demonstrate that for C3 symmetry, the coupling of O-P-O angles to phi is influence by repulsive steric interactions.

Conformational studies of nucleic acids: IV. The conformational energetics of oligonucleotides: d(ApApApA) and ApApApA

Utilizing a new method for modeling furanose pseudorotation (D. A. Pearlman and S.-H. Kim, J. Biomol. Struct. Dyn. 3, 85 (1985)) and the empirical multiple correlations between nucleic acid torsion angles we derived in the previous report (D. A. Pearlman and S.-H. Kim, previous paper in this issue), we have made an energetic examination of the entire conformational spaces available to two nucleic acid oligonucleotides: d(ApApApA) and ApApApA. The energies are calculated using a semi-empirical potential function. From the resulting body of data, energy contour map pairs (one for the DNA molecule, one for the RNA structure) have been created for each of the 21 possible torsion angle pairs in a nucleotide repeating unit. Of the 21 pairs, 15 have not been reported previously. The contour plots are different from those made earlier in that for each point in a particular angle-angle plot, the remaining five variable torsion angles are rotated to the values which give a minimum energy at this point. The contour maps are overall quite consistent with the experimental distribution of oligonucleotide data. A number of these maps are of particular interest: delta (C5'-C4'-C3'-O3')-chi (O4'-C1'-N9-C4), where the energetic basis for an approximately linear delta-chi correlation can be seen: zeta (C3'-O3'-P-O5')-delta, in which the experimentally observed linear correlation between zeta and delta in DNA(220 degrees less than zeta less than 280 degrees) is clearly predicted; zeta-epsilon (C4'-C3'-O3'-P), which shows that epsilon increases with decreasing zeta less than 260 degrees; alpha (O3'-P-O5'-C5')-gamma (O5'-C5'-C4'-C3') where a clear linear correlation between these angles is also apparent, consistent with experiment; and several others. For the DNA molecule studied here, the sugar torsion delta is predicted to be the most flexible, while for the RNA molecule, the greatest amount of flexibility is expected to reside in alpha and gamma. Both the DNA and RNA molecules are predicted to be highly polymorphic. Complete energy minimization has been performed on each of the minima found in the energy searches and the results further support this prediction. Possible pathways for B-form to A-form DNA interconversion suggested by the results of this study are discussed. The results of these calculations support use of the new sugar modeling technique and torsion angle correlations in future conformational studies of nucleic acids.

Base-pairing properties of O-methylated bases of nucleic acids. Energetic and steric considerations

Base-pairing properties of O-methylated nucleic-acid bases have been systematically investigated using both semi-empirical quantum-mechanical methods and a second-order perturbation formalism. The energetic, steric and electronic properties of (a) the individual methylated bases, (b) possible base-pairs formed between O-methylated and normal bases, and (c) mini-helices incorporating O-methylated bases were calculated. Two types of base-paired complexes were obtained: Those involving classical linear hydrogen bonds, and those involving bifurcated hydrogen-donor-hydrogen-acceptor interactions. In most complexes the presence of mispairs in the helical structure of nucleic acids is expected to create a local perturbation in the structure of the helix. Even though the most stable planar configurations of the mispairs may deviate markedly from those in the regular double helix, the induced deformations in the structure of the backbone are relatively small. Internal energies and geometries of mispairs are strongly affected by the conformation of the exocyclic group of the methylated bases. Another important contribution to the stability of various base-pairing schemes comes from stacking interactions.

A thermodynamic theory of codon bias in viral genes

The relationship between degeneracy in the genetic code and the occurrence of a strong codon bias is examined, with particular reference to a group of viral genomes. The present paper shows how codon bias may have been imposed by thermodynamic considerations at the time the primitive DNA first formed in the primordial soup. Using a four-state Ising-like model with stacking interactions between successive base pairs, we show how primeval periodic DNA polymers could have arisen the remnants of which are still observed in codon biases today.

Molecular electrostatic potential of the nucleic acids

It is generally acknowledged that geometrical and conformational properties of biopolymers have an important effect on their biochemical behaviour. It is less easily recognized that these properties depend also on their macromolecular electronic characteristics.

The aim of this review is to demonstrate the significance of such macromolecular electronic effects. Particularly useful for this sake is the recently much developed concept of ‘molecular electrostatic potential’ (MEP) (Scrocco & Tomasi, 1973, 1978) by which is defined the electrostatic (Coulomb) potential created in the neighbouring space by the nuclear charges and the eletronic distribution of a molecule.

Quantum chemical calculations on the two-step mechanism of proflavin binding to DNA

Quantum chemical calculations on the binding of proflavin to DNA lead to a model in which the outside binding to a phosphate group leads to an induced fit in the intercalation receptor site. The calculations suggest hydrogen bonding of the amine groups of the outside bound proflavin to the anionic oxygen of the backbone phosphate. The resulting partial neutralization facilitates the conformational transitions required for intercalation. The results are consistent with the observed preference of proflavin for dCpdG over dGpdC sequences and with the observed kinetics of the binding reaction.

Comparison of theoretical and experimental approaches to determination of conformation of nucleosides about the glycosidic bond

A study has been made by means of 1H-NMR spectroscopy of the syn in equilibrium anti dynamic equilibrium about the glycosidic bond for 5'-deoxyadenosine and some 8-substituted analogues, in different solvents. The results are compared with those previously obtained for the parent adenosine and its 8-substituted analogues. Quantum chemical calculations, with the aid of the Classical Potential and PCILO procedures, were applied to obtain the energies for different conformations of the base in adenosine and 5'-deoxyadenosine, and their 8-methyl and 8-halogeno derivatives. Good agreement was found between experimentally determined conformations in solution and those corresponding theoretically to the energy minima, particularly those calculated by the PCILO method. Comparison of the quantitative experimental data with the theoretical results was used to evaluate the validity of the latter and their applicability to studies of nucleoside conformation. The experimental and theoretical findings pointed to the existence of a marked flexibility about the glycosidic bond of the parent nucleosides and their 8-substituted analogues, when the 8-substituents were not too bulky, such as methyl or bromine. Considerations is given to possible correlations between conformational parameters in nucleosides and their 5'-deoxy analogues. It is shown that the proposed stabilization of the conformation syn by intramolecular hydrogen bonding, 5'-OH...N(3), is not in accord with the results of the present study.

The electrostatic molecular potential of yeast tRNAPhe. (I). The potential due to the phosphate backbone

We present a preliminary theoretical study of the electrostatic potential surrounding yeast tRNAPhe by computing the component of this potential due to all the 76 phosphate groups of the molecule. The general features of the results obtained are discussed and deductions concerning the binding of Mg2+ ions to the molecule in relation to the potential are presented.

A conformational rationale for the origin of the mechanism of nucleicacid-directed protein synthesis of 'living' organisms

The physical basis for the natural evolution of a primitive decoding system is presented using the concepts of molecular interactions. Oligoribonucleotides of five residues having U at the 5'-end, a purine at the 3'-end and any combination of three bases in the middle is taken as a primitive tRNA (PIT). From conformational considerations PIT is expected to have U-turn conformation wherein, N3-H3 of base U hydrogen-bonds with phosphate, three residues ahead leaving triplet bases called primitive anticodons (PAC) into a helical conformation, and this creates a cleft between U and PAC. An amino acid can be comfortably nestled into the cleft with the amide hydrogens and carboxyl oxygen hydrogen-bonded to the last purine and the first uridine, while the side-chain can interact with the cleft side of PAC. The other side of PAC is free to base-pair with triplet codons on a longer RNA. Also two PACs can 'recognize' consecutive triplet codons, and this leads to a dynamic interaction in which the amino and carboxyl ends are brought into proximity, making the formation of peptide bond feasible. The cleft formed by different anticodon triplets, broadly speaking, shows preferences for the corresponding amino acids of the presently known codon assignment. Thus the nucleicacid-directed protein synthesis, which is a unique feature of all 'living' organisms is shown to be a natural consequence of a particular way of favourable interaction between nucleic acids and amino acids, and our model provides the missing link between the chemical evolution of small organic molecules and biological evolution through the process of mutations in nucleicacids and nucleicacid-directed protein synthesis.

Origins of the specificity in the intercalation of ethidium into nucleic acids. A theoretical analysis

The sequence preferences observed in the intercalative binding of ethidium to dinucleoside phosphates have been theoretically examined. This specificity is for pyrimidine (3'--5') purine sequences as compared to their purine (3'--5') pyrimidine sequence isomers. It is shown that the stacking energies between the ethidium cation and the base pairs are fairly constant for all combinations of bases at the intercalation site. In contrast, the energy of unwinding the double helix to assume the geometry of the intercalation complex shows substantial sequence differences. Thus, the specificity observed is more readily explained in terms of these conformation energy changes than by preferential stacking interactions. These results imply that there may be a large class of intercalating molecules which exhibit similar pyrimidine (3'--5') purine sequence specificity.

Molecular orbital studies on nucleoside analogs. II. Conformation of 6-azapyrimidine nucleosides

PCILO (Perturbative Configuration Interaction using Localised Orbitals) computations have been carried out for three 6-azapyrimidine nucleosides, 6-azauridine, 6-azacytidine and 6-azathymidine, for both C(2')-endo and C(3')-endo pucker of the sugar ring. The results indicate a syn (chiCN=180 degrees) conformation followed by chiCN=90 degrees and gg conformation for C(3')-endo 6-aza analogs as compareed to the anti (chiCN=0 degrees) and gg conformation preferred by the corresponding pyrimidine nucleosides. For C(2')-endo sugar geometry, 6-azauridine and 6-azacytidine prefer, respectively, chiCN=0 degrees (anti) and phi C(4')-C(5')=60 degrees C (gg) and chiCN-240 degrees (syn) and phi C(4')-C(5')=120 degrees. The corresponding nucleosides, uridine and cytidine, show a preference for syn (chiCN=240 degrees) and gg and anti(chiCN=0 degrees) and gg , respectively. The X-ray crystallographic conformations of 6-azauridine and 6-azacytidine have been attributed to intermolecular hydrogen bonding and crystal packing forces. The results of PMR, CD and ORD studies on 6-azauridine and 6-azacytidine in aqueous solutions are in agreement with the PCILO predictions.