A computational study of expanded heterocyclic nucleosides in DNA

Probing the Effects of Pyrimidine Functional Group Switches on Acyclic Fleximer Analogues for Antiviral Activity

Due to their ability to inhibit viral DNA or RNA replication, nucleoside analogues have been used for decades as potent antiviral therapeutics. However, one of the major limitations of nucleoside analogues is the development of antiviral resistance. In that regard, flexible nucleoside analogues known as “fleximers” have garnered attention over the years due to their ability to survey different amino acids in enzyme binding sites, thus overcoming the potential development of antiviral resistance. Acyclic fleximers have previously demonstrated antiviral activity against numerous viruses including Middle East Respiratory Syndrome coronavirus (MERS-CoV), Ebola virus (EBOV), and, most recently, flaviviruses such as Dengue (DENV) and Yellow Fever Virus (YFV). Due to these interesting results, a Structure Activity Relationship (SAR) study was pursued in order to analyze the effect of the pyrimidine functional group and acyl protecting group on antiviral activity, cytotoxicity, and conformation. The results of those studies are presented herein.

How do hydrophobic nucleobases differ from natural DNA nucleobases? Comparison of structural features and duplex properties from QM calculations and MD simulations

Computational (DFT and MD simulation) methods are employed to systematically characterize the structural and energetic properties of five hydrophobic nucleobases (FEMO, MMO2, NaM, 5SICS and TPT3) that constitute four unnatural base pairs (FEMO:5SICS, MMO2:5SICS, NaM:5SICS and TPT3:NaM). These hydrophobic bases have been recently shown to be replicated when present between natural bases in DNA duplexes, with the highest replication fidelity and efficiency occuring for the TPT3:NaM pair. Our QM calculations suggest that the preferred (anti) glycosidic orientations of nucleosides containing hydrophobic bases are similar to the natural DNA nucleosides despite differences in their chemical structures. However, due to the inability to form interbase hydrogen bonds, hydrophobic base pairs intrinsically prefer nonplanar, distorted geometries, many of which are stabilized through π-π stacking interactions. Furthermore, the intrinsic stacking potential between a hydrophobic and a natural base is similar to that between two natural bases, indicating that the strength of stacking interactions in DNA duplexes containing hydrophobic bases is likely comparable to natural DNA. However, in contrast to the isolated base-pair geometries, our MD simulations suggest that the hydrophobic base pairs adopt variable geometries within DNA, which range from stacked (5SICS:FEMO) to nearly planar (5SICS:NaM and SICS:MMO2) to planar (TPT3:NaM). As a result, the duplex structural features at the site of modification depend on the identity of the hydrophobic base pair, where the TPT3:NaM pair causes the least structural changes compared to natural DNA. Overall, the structural insight obtained from our calculations on DNA containing hydrophobic base pairs explains the experimentally-observed higher fidelity and efficiency during replication of TPT3:NaM compared to other hydrophobic nucleobase pairs. By providing valuable structural information that explains the intrinsic and duplex properties of this class of unnatural nucleobases, the present work may aid the future design of improved hydrophobic analogues.

Thiophene-expanded guanosine analogues of Gemcitabine

The chemotherapeutic drug Gemcitabine, 2',2'-difluoro-2'-deoxycytidine, has long been the standard of care for a number of cancers. Gemcitabine's chemotherapeutic properties stem from its 2',2'-difluoro-2'-deoxyribose sugar, which mimics the natural nucleoside, but also disrupts nucleic acid synthesis, leading to cell death. As a result, numerous analogues have been prepared to further explore the biological implications for this structural modification. In that regard, a thieno-expanded guanosine analogue was of interest due to biological activity previously observed for the tricyclic heterobase scaffold. Several analogues were prepared, including the McGuigan ProTide, however the parent nucleoside exhibited the best chemotherapeutic activity, specifically against breast cancer cell lines (89.53% growth inhibition).

Bicyclic and Tricyclic "Expanded" Nucleobase Analogues of Sofosbuvir: New Scaffolds for Hepatitis C Therapies

Given the impressive success of Gilead's Sofosbuvir, many laboratories, including ours, have explored the unique 2'-sugar modification (2'-Me, 2'-F) of nucleoside analogues in the hopes of exploiting the biological activity that this unique modification has imparted to the nucleoside scaffold. In that regard, we have combined our tricyclic "expanded" purine base motif with the 2'-Me, 2'-F sugar modification. Although the synthesis of these complex molecules proved to be nontrivial, with the best results coming from a linear approach, the overall strategy resulted in highly promising biological results for several of the target compounds, including their corresponding McGuigan ProTides. Modest activity against HCV was observed with inhibitory concentrations of as low as 20 μM.

yDNA versus yyDNA pyrimidines: computational analysis of the effects of unidirectional ring expansion on the preferred sugar-base orientation, hydrogen-bonding interactions and stacking abilities

The properties of natural, y- and yy-pyrimidines are compared using computational (B3LYP, MP2) methods. Ring expansion upon incorporation of benzene or naphthalene into the natural pyrimidines affects the preferred orientation of the base about the glycosidic bond in the corresponding nucleoside to a similar extent. Specifically, although the natural pyrimidines preferentially adopt the anti orientation with respect to the 2'-deoxyribose moiety, the expanded analogues will likely display (anti/syn) conformational heterogeneity, which may lead to alternate hydrogen-bonding modes in double-stranded duplexes. Nevertheless, the A:T Watson-Crick hydrogen-bond strengths do not significantly change upon base expansion, while the G:C interaction energy is slightly strengthened upon incorporation of either expanded pyrimidine. The largest effect of base expansion occurs in the stacking energies. Specifically, the maximum (most negative) stacking energies in isolated dimers formed by aligning the nucleobase centers of mass can be increased up to 45% by inclusion of a single y-pyrimidine and up to 55% by consideration of a yy-pyrimidine. Similar increases in the stacking interactions are found when a simplified duplex model composed of two stacked (hydrogen-bonded) base pairs is considered, where both the intrastrand and interstrand stacking interactions can be increased and the effects are more pronounced for the yy-pyrimidines. Moreover, the total stability (sum of all hydrogen-bonding and stacking interactions) is greater for duplexes containing expanded yy-pyrimidines compared to y-pyrimidines, which is mainly due to enhanced stacking interactions. Thus, our calculations suggest that multiple unidirectional increases in the size of the nucleobase spacer can continuously enhance the stability of expanded duplexes.

On the structure and dynamics of duplex GNA

Glycol nucleic acid (GNA), with a nucleotide backbone comprising of just three carbons and the stereocenter derived from propylene glycol (1,2-propanediol), is a structural analog of nucleic acids with intriguing biophysical properties, such as formation of highly stable antiparallel duplexes with high Watson-Crick base pairing fidelity. Previous crystallographic studies of double stranded GNA (dsGNA) indicated two forms of backbone conformations, an elongated M-type (containing metallo-base pairs) and the condensed N-type (containing brominated base pairs). A herein presented new crystal structure of a GNA duplex at 1.8 Å resolution from self-complementary 3'-CTC(Br)UAGAG-2' GNA oligonucleotides reveals an N-type conformation with alternating gauche-anti torsions along its (O3'-C3'-C2'-O2') backbone. To elucidate the conformational state of dsGNA in solution, molecular dynamic simulations over a period of 20 ns were performed with the now available repertoire of structural information. Interestingly, dsGNA adopts conformational states in solution intermediate between experimentally observed backbone conformations: simulated dsGNA shows the all-gauche conformation characteristic of M-type GNA with the higher helical twist common to N-type GNA structures. The so far counterintuitive, smaller loss of entropy upon duplex formation as compared to DNA can be traced back to the conformational flexibility inherent to dsGNA but missing in dsDNA. Besides extensive interstrand base stacking and conformational preorganization of single strands, this flexibility contributes to the extraordinary thermal stability of GNA.

Mechanistic studies in the synthesis of a series of thieno-expanded xanthosine and guanosine nucleosides

During the synthetic pursuit of guanosine (triG) and xanthosine (triX) tricyclic nucleosides analogues, an interesting side product was discovered. In an effort to uncover the mechanistic factors leading to this result, a series of reaction conditions were investigated. It was found that by varying the conditions, the appearance of the side product could be controlled. In addition, the yield of the desired products could be manipulated to afford either a 50:50 mix of both triG and triX, or a majority of one or the other. To demonstrate the broad utility of the method, it was also adapted to the synthesis of guanosine and xanthosine from 5-amino-1-beta-D-ribofuranosyl-4-imidazolecarboxyamide (AICAR). The mechanistic details surrounding the synthetic efforts are reported herein.