Synthesis of 3′,4′-difluoro-3′-deoxyribonucleosides and its evaluation of the biological activities: Discovery of a novel type of anti-HCV agent 3′,4′-difluorocordycepin

Upon reacting 3′,4′-unsaturated cytosine (8 and 9) and adenine nucleosides (13 and 14) with XeF₂/BF₃·OEt₂, the respective novel 3′,4′-difluoro-3′-deoxyribofuranosyl nucleosides (10–12 and 15–18) could be obtained. Formation of anti-adducts (11, 16 and 18) revealed that the fluorination involved oxonium ions as incipient intermediates. TBDMS-protected 3′,4′-unsaturated adenosine provided the β-face adducts as sole stereoisomers whereas α-face-selectivity was observed with the TBDPS-protected adenosine 14. The evaluation of the novel 3′-deoxy-3′,4′-difluororibofuranosylcytosine-(19–21) and adenine nucleosides (22–25) against antitumor and antiviral activities revealed that 3′,4′-difluorocordycepin (24) was found to possess anti-HCV activity. The SI of 24 was comparable to that of the anti-HCV drug ribavirin. However, sofosbuvir, FDA-approved novel anti-HCV drug, showed better SI value. Our finding revealed that the introduction of the fluoro-substituent into the 4′-position of cordycepin derivatives decreased the cytotoxicity to the host cell with retention of the antiviral activity.     

A simulation approach for power calculation in large cohort studies based on multistate models

Realistic power calculations for large cohort studies and nested case control studies are essential for successfully answering important and complex research questions in epidemiology and clinical medicine. For this, we provide a methodical framework for general realistic power calculations via simulations that we put into practice by means of an R‐based template. We consider staggered recruitment and individual hazard rates, competing risks, interaction effects, and the misclassification of covariates. The study cohort is assembled with respect to given age‐, gender‐, and community distributions. Nested case‐control analyses with a varying number of controls enable comparisons of power with a full cohort analysis. Time‐to‐event generation under competing risks, including delayed study‐entry times, is realized on the basis of a six‐state Markov model. Incidence rates, prevalence of risk factors and prefixed hazard ratios allow for the assignment of age‐dependent transition rates given in the form of Cox models. These provide the basis for a central simulation‐algorithm, which is used for the generation of sample paths of the underlying time‐inhomogeneous Markov processes. With the inclusion of frailty terms into the Cox models the Markov property is specifically biased. An “individual Markov process given frailty” creates some unobserved heterogeneity between individuals. Different left‐truncation‐ and right‐censoring patterns call for the use of Cox models for data analysis. p‐values are recorded over repeated simulation runs to allow for the desired power calculations. For illustration, we consider scenarios with a “testing” character as well as realistic scenarios. This enables the validation of a correct implementation of theoretical concepts and concrete sample size recommendations against an actual epidemiological background, here given with possible substudy designs within the German National Cohort.     

Potential pitfalls of the PTV concept in dose-to-medium planning optimization

In typical treatment planning of 3D IMRT, the incident energy fluence is optimized to achieve a homogeneous dose distribution to the PTV. The PTV includes the tumour but also healthy tissues that may have a different dose response for the same incident energy fluence, like bony structures included in the PTV (mandibles in head and neck tumours or femoral bones in sarcomas). Dose to medium optimization compensates for this heterogeneous response, leading to a non-homogeneous energy fluence in the PTV and a non-homogeneous dose in the CTV in the presence of geometric errors. We illustrate qualitatively this statement in a cylindrical geometry where the PTV includes a CTV (7 cm diameter) made of water surrounded by ICRU compact bone (1.2 cm thickness); such configuration was chosen to exaggerate the aforementioned effect. Optimization was performed assuming dose equals photon energy fluence times mass energy absorption coefficient. Bone has a 4% lower dose response in a 6 MV flattening filter free spectrum. After optimization either in medium or assuming everything as water composition, the geometry was shifted by 1.2 cm and dose recomputed. As expected, compensating for the under-response of the bone material during optimization in medium leads to an overdosage of the CTV when patient geometric errors are taken into account. Optimization in dose assuming everything as water composition leads to a uniform coverage. Robust optimization or forcing a uniform atomic composition in the PTV margin may resolve this incompatibility between the PTV concept and dose to medium optimization.     

Novel Monte Carlo dose calculation algorithm for robotic radiosurgery with multi leaf collimator: Dosimetric evaluation

PurposeAt introduction in 2014, dose calculation for the first MLC on a robotic SRS/SBRT platform was limited to a correction-based Finite-Size Pencil Beam (FSPB) algorithm. We report on the dosimetric accuracy of a novel Monte Carlo (MC) dose calculation algorithm for this MLC, included in the Precision™ treatment planning system.MethodsA phantom was built of one slab (5.0 cm) of lung-equivalent material (0.09…0.29 g/cc) enclosed by 3.5 cm (above) and 5 cm (below) slabs of solid water (1.045 g/cc). This was irradiated using rectangular (15.4 × 15.4 mm² to 53.8 × 53.7 mm²) and two irregular MLC-fields. Radiochromic film (EBT3) was positioned perpendicular to the slabs and parallel to the beam. Calculated dose distributions were compared to film measurements using line scans and 2D gamma analysis.ResultsMeasured and MC calculated percent depth dose curves showed a characteristic dose drop within the low-density region, which was not correctly reproduced by FSPB. Superior average gamma pass rates (2%/1 mm) were found for MC (91.2 ± 1.5%) compared to FSPB (55.4 ± 2.7%). However, MC calculations exhibited localized anomalies at mass density transitions around 0.15 g/cc, which were traced to a simplification in electron transport. Absence of these anomalies was confirmed in a modified build of the MC engine, which increased gamma pass rates to 96.6 ± 1.2%.ConclusionsThe novel MC algorithm greatly improves dosimetric accuracy in heterogeneous tissue, potentially expanding the clinical use of robotic radiosurgery with MLC. In-depth, independent validation is paramount to identify and reduce the residual uncertainties in any software solution.     

UV resonance Raman studies on the activation mechanism of human hematopoietic prostaglandin D2 synthase by a divalent cation, Mg

The Mg²⁺ ion-assisted activation mechanism of the active site Tyr8 of a human hematopoietic prostaglandin D₂ synthase (H-PGDS) was studied by ultraviolet resonance Raman (UVRR) spectroscopy. Addition of Mg²⁺ to the native H-PGDS at pH 8.0 resulted in the Y8a Raman band of Tyr8 shifting from 1615 cm⁻¹ to 1600 cm⁻¹. This large shift to lower energy of the tyrosine Y8a vibrational mode is caused by the deprotonation of the tyrosine phenol group promoted by binding of Mg²⁺. Upon subsequent addition of glutathione (GSH), the Mg²⁺/H-PGDS solution showed the Tyr8 Raman band shifted to 1611 cm⁻¹, which is 11 cm⁻¹ higher than the frequency of the Mg²⁺ complex of H-PGDS, but 4 cm⁻¹ lower than the Mg²⁺ free enzyme. These UVRR observations suggest that the deprotonated Tyr8 in the presence of Mg²⁺ is re-protonated by the abstraction of H⁺ from the thiol group of GSH, and that the re-protonated Tyr8 species forms a hydrogen bond with the thiolate anion of GSH. Density functional theory calculations on several model complexes of p-cresol were also performed, which suggested that the pK_a and vibrational frequencies of the Tyr8 phenol group are affected by the degree and structure of hydration of the Tyr8 residue.     

Possible existence of ion pairs at the mouths of ion channels

An electrostatic calculation suggests that when an ion is bound near the mouth of a channel penetrating a low-dielectric membrane, a counter ion may form an ion pair with this ion. The tendency towards ion-pair formation is remarkably enhanced at channel mouths by forces (image forces) arising from the charges induced on the boundaries between different dielectrics. The binding constant for the formation of ion-pairs of monovalent ions is estimated under the assumption that local interactions between the counter ion and the channel wall are negligibly small. It is of the order of 1–10 molal^?1 or more for the binding of a Cl^? (F^?) counter ion to an $\text{Na}{}^{\mathrm{+}}$ ($\text{Li}{}^{\mathrm{+}}$) ion if appropriate conditions are fulfilled. The binding constant depends on the position of the binding site, the dimensions and geometries of the channel and channel mouth, and the state of ion loading of the channel, as well as the ionic species. The present results also indicate that when cation (anion) channels have anionic (cationic) groups as integrant parts of their channel walls, interactions between these charged groups and permeant ions are markedly enhanced by the image forces.     

The role of the medium in long-range electron transfer

 The role of the polypeptide matrix in electron transfer processes in proteins has been studied in two distinct systems: first in a protein where the induced ET is artificial, and second as part of the catalytic cycle of an enzyme. Azurins are structurally well-characterized blue single-copper proteins consisting of a rigid β-sheet polypeptide matrix. We have determined rate constants and activation parameters for intramolecular long-range electron transfer between the disulfide radical anions (generated by pulse radiolysis) and the copper(II) centre as a function of driving force and nature of the intervening medium in a large number of wild-type and single-site-mutated proteins. In ascorbate oxidase, for which the three-dimensional structure is equally well characterized, the internal ET from the type-I Cu(I) to the trinuclear Cu(II) centre has been studied. We find that the results correlate well with distance through well-defined pathways using a through-bond electron tunnelling mechanism. Received: 2 January 1997 / Accepted: 6 February 1997     

Solution structures of purine base analogues 9-deazaguanine and 9-deazahypoxanthine

Deaza analogues of nucleobases are potential drugs against infectious diseases caused by parasites. A caveat is that apart from binding their target parasite enzymes, they also bind and inhibit enzymes of the host. In order to design derivatives of deaza analogues which specifically bind target enzymes, knowledge of their molecular structure, protonation state, and predominant tautomers at physiological conditions is essential. We have employed resonance Raman spectroscopy at an excitation wavelength of 260 nm, to decipher solution structure of 9-deazaguanine (9DAG) and 9-deazahypoxanthine (9DAH). These are analogues of guanine and hypoxanthine, respectively, and have been exploited to study static complexes of nucleobase binding enzymes. Such enzymes are known to perturb pK_a of their ligands, and thus, we also determined solution structures of these analogues at two, acidic and alkaline, pH. Structure of each possible protonation state and tautomer was computed using density functional theoretical calculations. Species at various pHs were identified based on isotopic shifts in experimental wavenumbers and by comparing these shifts with corresponding computed isotopic shifts. Our results show that at physiological pH, N1 of pyrimidine ring in 9DAG and 9DAH bears a proton. At lower pH, N3 is place of protonation, and at higher pH, deprotonation occurs at N1 position. The proton at N7 of purine ring remains intact even at pH 12.5. We have further compared these results with naturally occurring nucleotides. Our results identify key vibrational modes which can report on hydrogen bonding interactions, protonation and deprotonation in purine rings upon binding to the active site of enzymes.     

Computer modeling of substrate binding to lipases from Rhizomucor miehei, Humicola lanuginosa, and Candida rugosa.

The substrate-binding sites of the triacyl glyceride lipases from Rhizomucor miehei, Humicola lanuginosa, and Candida rugosa were studied by means of computer modeling methods. The space around the active site was mapped by different probes. These calculations suggested 2 separate regions within the binding site. One region showed high affinity for aliphatic groups, whereas the other region was hydrophilic. The aliphatic site should be a binding cavity for fatty acid chains. Water molecules are required for the hydrolysis of the acyl enzyme, but are probably not readily accessible in the hydrophobic interface, in which lipases are acting. Therefore, the hydrophilic site should be important for the hydrolytic activity of the enzyme. Lipases from R. miehei and H. lanuginosa are excellent catalysts for enantioselective resolutions of many secondary alcohols. We used molecular mechanics and dynamics calculations of enzyme-substrate transition-state complexes, which provided information about molecular interactions important for the enantioselectivities of these reactions.