Evaluation of the impact of single-nucleotide polymorphisms on treatment response, survival and toxicity with cytarabine and anthracyclines in patients with acute myeloid leukaemia: a systematic review protocol

Butterfly Effect in Cytarabine: Combined NMR-NQR Experiment, Solid-State Computational Modeling, Quantitative Structure-Property Relationships and Molecular Docking Study

Cytarabine (Ara-C) is a synthetic isomer of cytidine that differs from cytidine and deoxycytidine only in the sugar. The use of arabinose instead of deoxyribose hinders the formation of phosphodiester linkages between pentoses, preventing the DNA chain from elongation and interrupting the DNA synthesis. The minor structural alteration (the inversion of hydroxyl at the 2' positions of the sugar) leads to change of the biological activity from anti-depressant and DNA/RNA block builder to powerful anti-cancer. Our study aimed to determine the molecular nature of this phenomenon. Three 1H-14N NMR-NQR experimental techniques, followed by solid-state computational modelling (Quantum Theory of Atoms in Molecules, Reduced Density Gradient and 3D Hirshfeld surfaces), Quantitative Structure-Property Relationships, Spackman's Hirshfeld surfaces and Molecular Docking were used. Multifaceted analysis-combining experiments, computational modeling and molecular docking-provides deep insight into three-dimensional packing at the atomic and molecular levels, but is challenging. A spectrum with nine lines indicating the existence of three chemically inequivalent nitrogen sites in the Ara-C molecule was recorded, and the lines were assigned to them. The influence of the structural alteration on the NQR parameters was modeled in the solid (GGA/RPBE). For the comprehensive description of the nature of these interactions several factors were considered, including relative reactivity and the involvement of heavy atoms in various non-covalent interactions. The binding modes in the solid state and complex with dCK were investigated using the novel approaches: radial plots, heatmaps and root-mean-square deviation of the binding mode. We identified the intramolecular OH···O hydrogen bond as the key factor responsible for forcing the glycone conformation and strengthening NH···O bonds with Gln97, Asp133 and Ara128, and stacking with Phe137. The titular butterfly effect is associated with both the inversion and the presence of this intramolecular hydrogen bond. Our study elucidates the differences in the binding modes of Ara-C and cytidine, which should guide the design of more potent anti-cancer and anti-viral analogues.

Role of Pharmacogenetics in the Treatment of Acute Myeloid Leukemia: Systematic Review and Future Perspectives

Acute myeloid leukemia (AML) is a heterogeneous disease characterized by remarkable toxicity and great variability in response to treatment. Plenteous pharmacogenetic studies have already been published for classical therapies, such as cytarabine or anthracyclines, but such studies remain scarce for newer drugs. There is evidence of the relevance of polymorphisms in response to treatment, although most studies have limitations in terms of cohort size or standardization of results. The different responses associated with genetic variability include both increased drug efficacy and toxicity and decreased response or resistance to treatment. A broad pharmacogenetic understanding may be useful in the design of dosing strategies and treatment guidelines. The aim of this study is to perform a review of the available publications and evidence related to the pharmacogenetics of AML, compiling those studies that may be useful in optimizing drug administration.

Polygenic Ara-C Response Score Identifies Pediatric Patients With Acute Myeloid Leukemia in Need of Chemotherapy Augmentation

PURPOSE

To establish a patient-specific polygenic score derived from cytarabine (ara-C) pathway pharmacogenomic evaluation to personalize acute myeloid leukemia (AML) treatment.

MATERIALS AND METHODS

Single nucleotide polymorphisms (SNPs) in the ara-C-pathway genes were analyzed with outcome in patients from the multicenter-AML02 trial (N = 166). Multi-SNP predictor modeling was used to develop 10-SNP Ara-C_SNP score (ACS10) using top SNPs predictive of minimal residual disease and event-free survival (EFS) from the AML02-cohort and four SNPs previously associated with ara-C triphosphate levels in the AML97 trial. ACS10 was evaluated for association with outcomes in each clinical trial arms: the standard low-dose ara-C (LDAC, n = 91) and augmented high-dose ara-C (HDAC, n = 75) arms of AML02 and the standard Ara-C, daunorubicin and etoposide (ADE) (n = 465) and the augmented ADE + gemtuzumab ozogamicin (GO; n = 466) arms of AAML0531 trial.

RESULTS

In the standard LDAC-arm of AML02 cohort, the low-ACS10 score group (≤ 0) had significantly worse EFS (ACS10 low v high hazard ratio [HR] = 2.81; 95% CI, 1.45 to 5.43; P = .002) and overall survival (OS; HR = 2.98; 95% CI, 1.32 to 6.75; P = .009) compared with the high-ACS10 group (score > 0). These results were validated in the standard-ADE arm of AAML0531, with poor outcome in the low-ASC10 group compared with the high-ACS10 group (EFS: HR = 1.35, 95% CI, 1.04 to 1.75, P = .026; OS: HR = 1.64, 95% CI, 1.2 to 2.22, P = .002). Within the augmented arms (AML02-HDAC and AAML0531-ADE + GO), EFS and OS did not differ between low- and high-ACS10 score groups. In both cohorts, patients with low-ACS10 consistently showed a 10-percentage point improvement in 5-year EFS with augmented therapy (AML02-HDAC or AAML0531-ADE + GO arms) than with standard therapy (AML02-LDAC or AAML0531-ADE arms).

CONCLUSION

Patients with low-ACS10 score experienced significantly poor outcome when treated on standard regimen. Augmentation with either high-dose ara-C or GO addition improved outcome in low-ACS10 group. A polygenic ACS10 score can identify patients with unfavorable pharmacogenetic characteristics and offers a potential for an elective augmented therapy option.

Sequence analysis for SNP detection and phylogenetic reconstruction of SARS-cov-2 isolated from Nigerian COVID-19 cases

Background

Coronaviruses are a group of viruses that belong to the Family Coronaviridae, genus Betacoronavirus. In December 2019, a new coronavirus disease (COVID-19) characterized by severe respiratory symptoms was discovered. The causative pathogen was a novel coronavirus known as 2019-nCoV and later as SARS-CoV-2. Within two months of its discovery, COVID-19 became a pandemic causing widespread morbidity and mortality.

Methodology

Whole genome sequence data of SARS-CoV-2 isolated from Nigerian COVID-19 cases were retrieved by downloading from GISAID database. A total of 18 sequences that satisfied quality assurance (length ≥29,700 nts and number of unknown bases denoted as “N” ≤ 5%) were used for the study. In addition, genome sequence of SARS-CoV-2 obtained from Nigeria's COVID-19 index case (Accession ID: EPI_ISL_413550) and the reference genome (Accession NC_ 045512.2) were obtained from GISAID and the GenBank databases, respectively. Multiple sequence alignment (MSA) was done in MAFFT (Version 7.471) while SNP calling was implemented in DnaSP (Version 6.12.03), respectively and then visualized in Jalview (Version 2.11.1.0). Phylogenetic analysis was with MEGA X software.

Results

Nigerian SARS-CoV-2 had 99.9% genomic similarity with four large conserved genomic regions. A total of 66 SNPs were identified out of which 31 were informative. Nucleotide diversity assessment gave Pi = 0.00048 and average SNP frequency of 2.22 SNPs per 1000 nts. Non-coding genomic regions particularly 5′UTR and 3′UTR had a SNP density of 3.77 and 35.4, respectively. The region with the highest SNP density was ORF10 with a frequency of 8.55 SNPs/1000 nts). This value was significantly higher (P < 0.01) than that of the spike gene, the region of greatest interest in SARS-CoV-2 genomics. Majority (72.2%) of viruses in Nigeria are of L lineage with preponderance of D614G mutation which accounted for 11 (61.1%) out of the 18 viral sequences. Nigeria SARS-CoV-2 revealed 3 major clades namely Oyo, Ekiti and Osun on a maximum likelihood phylogenetic tree.

Conclusion and Recommendation

There was a preponderance of L lineage (to include the new lineage scheme) and D614G mutants. Nigerian SARS-CoV-2 genome revealed ORF1ab as the region containing the highest SNP density as compared to the spike gene. The implication of this distribution of SNPs for the empirical lower infectivity of SARS-CoV-2 in Nigeria is discussed. This also underscores the need for more aggressive testing and treatment of COVID-19 in Nigeria. Additionally, attempt to produce testing kits for COVID-19 in Nigeria should consider the conserved regions identified in this study. Strict adherence to COVID-19 preventive measure is recommended in view of Nigerian SARS-CoV-2 phylogenetic clustering pattern, which suggests intensive community transmission possibly rooted in communal culture characteristic of many ethnicities in Nigeria.

Overwhelming mutations or SNPs of SARS-CoV-2: A point of caution

The morbidity of SARS-CoV-2 (COVID-19) is reaching 3 Million landmark causing and a serious public health concern globally and it is enigmatic how several antiviral and antibody treatments were not effective in the different period across the globe. With the drastic increasing number of positive cases around the world WHO raised the importance in the assessment of the risk of spread and understanding genetic modifications that could have occurred in the SARS-CoV-2. Using all available deep sequencing data of complete genome from all over the world (NCBI repository), we identified several hundreds of point mutations or SNPs in SARS-CoV-2 all across the genome. This could be the cause for the constant change and differed virulence with an increase in mortality and morbidity. Among the 12 different countries (one sequence from each country) with complete genome sequencing data, we noted the 47 key point mutations or SNPs located along the entire genome that might have impact in the virulence and response to different antivirals against SARS-CoV-2. In this regard, key viral proteins of spike glycoprotein, Nsp1, RdRp and the ORF8 region got heavily mutated within these 3 months via person-to-person passage. We also discuss what could be the possible cause of this rapid mutation in the SARS-CoV-2.