Investigations Into Chemically Stabilized Four-Letter DNA for DNA-Encoded Chemistry

DNA-encoded libraries are a prime technology for target-based small molecule screening. Native DNA used as genetic compound barcode is chemically vulnerable under many reaction conditions. DNA barcodes that are composed of pyrimidine nucleobases, 7-deazaadenine, and 7-deaza-8-azaguanine have been investigated for their suitability for encoded chemistry both experimentally and computationally. These four-letter barcodes were readily ligated by T4 ligation, amplifiable by Taq polymerase, and the resultant amplicons were correctly sequenced. Chemical stability profiling showed a superior chemical stability compared to native DNA, though higher susceptibility to depurination than a three-letter code based on pyrimidine DNA and 7-deazaadenine.

Chemically Stabilized DNA Barcodes for DNA-Encoded Chemistry

DNA‐encoded compound libraries are a widely used small molecule screening technology. One important aim in library design is the coverage of chemical space through structurally diverse molecules. Yet, the chemical reactivity of native DNA barcodes limits the toolbox of reactions for library design. Substituting the chemically vulnerable purines by 7‐deazaadenine, which exhibits tautomerization stability similar to natural adenine with respect to the formation of stable Watson–Crick pairs, yielded ligation‐competent, amplifiable, and readable DNA barcodes for encoded chemistry with enhanced stability against protic acid‐ and metal ion‐promoted depurination. The barcode stability allowed for straightforward translation of 16 exemplary reactions that included isocyanide multicomponent reactions, acid‐promoted Pictet–Spengler and Biginelli reactions, and metal‐promoted pyrazole syntheses on controlled pore glass‐coupled barcodes for diverse DEL design. The Boc protective group of reaction products offered a convenient handle for encoded compound purification.

Quantum-mechanical property prediction of solvated drug molecules: what have we learned from a decade of SAMPL blind prediction challenges?

Joint academic-industrial projects supporting drug discovery are frequently pursued to deploy and benchmark cutting-edge methodical developments from academia in a real-world industrial environment at different scales. The dimensionality of tasks ranges from small molecule physicochemical property assessment over protein-ligand interaction up to statistical analyses of biological data. This way, method development and usability both benefit from insights gained at both ends, when predictiveness and readiness of novel approaches are confirmed, but the pharmaceutical drug makers get early access to novel tools for the quality of drug products and benefit of patients. Quantum-mechanical and simulation methods particularly fall into this group of methods, as they require skills and expense in their development but also significant resources in their application, thus are comparatively slowly dripping into the realm of industrial use. Nevertheless, these physics-based methods are becoming more and more useful. Starting with a general overview of these and in particular quantum-mechanical methods for drug discovery we review a decade-long and ongoing collaboration between Sanofi and the Kast group focused on the application of the embedded cluster reference interaction site model (EC-RISM), a solvation model for quantum chemistry, to study small molecule chemistry in the context of joint participation in several SAMPL (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenges. Starting with early application to tautomer equilibria in water (SAMPL2) the methodology was further developed to allow for challenge contributions related to predictions of distribution coefficients (SAMPL5) and acidity constants (SAMPL6) over the years. Particular emphasis is put on a frequently overlooked aspect of measuring the quality of models, namely the retrospective analysis of earlier datasets and predictions in light of more recent and advanced developments. We therefore demonstrate the performance of the current methodical state of the art as developed and optimized for the SAMPL6 pKa and octanol-water log P challenges when re-applied to the earlier SAMPL5 cyclohexane-water log D and SAMPL2 tautomer equilibria datasets. Systematic improvement is not consistently found throughout despite the similarity of the problem class, i.e. protonation reactions and phase distribution. Hence, it is possible to learn about hidden bias in model assessment, as results derived from more elaborate methods do not necessarily improve quantitative agreement. This indicates the role of chance or coincidence for model development on the one hand which allows for the identification of systematic error and opportunities toward improvement and reveals possible sources of experimental uncertainty on the other. These insights are particularly useful for further academia-industry collaborations, as both partners are then enabled to optimize both the computational and experimental settings for data generation.

Tautomeric Equilibria of Nucleobases in the Hachimoji Expanded Genetic Alphabet

Evolution has yielded biopolymers that are constructed from exactly four building blocks and are able to support Darwinian evolution. Synthetic biology aims to extend this alphabet, and we recently showed that 8-letter (hachimoji) DNA can support rule-based information encoding. One source of replicative error in non-natural DNA-like systems, however, is the occurrence of alternative tautomeric forms, which pair differently. Unfortunately, little is known about how structural modifications impact free-energy differences between tautomers of the non-natural nucleobases used in the hachimoji expanded genetic alphabet. Determining experimental tautomer ratios is technically difficult, and so, strategies for improving hachimoji DNA replication efficiency will benefit from accurate computational predictions of equilibrium tautomeric ratios. We now report that high-level quantum-chemical calculations in aqueous solution by the embedded cluster reference interaction site model, benchmarked against free-energy molecular simulations for solvation thermodynamics, provide useful quantitative information on the tautomer ratios of both Watson-Crick and hachimoji nucleobases. In agreement with previous computational studies, all four Watson-Crick nucleobases adopt essentially only one tautomer in water. This is not the case, however, for non-natural nucleobases and their analogues. For example, although the enols of isoguanine and a series of related purines are not populated in water, these heterocycles possess N₁-H and N₃-H keto tautomers that are similar in energy, thereby adversely impacting accurate nucleobase pairing. These robust computational strategies offer a firm basis for improving experimental measurements of tautomeric ratios, which are currently limited to studying molecules that exist only as two tautomers in solution.

The SAMPL6 challenge on predicting octanol-water partition coefficients from EC-RISM theory

Results are reported for octanol–water partition coefficients (log P) of the neutral states of drug-like molecules provided during the SAMPL6 (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenge from applying the “embedded cluster reference interaction site model” (EC-RISM) as a solvation model for quantum-chemical calculations. Following the strategy outlined during earlier SAMPL challenges we first train 1- and 2-parameter water-free (“dry”) and water-saturated (“wet”) models for n-octanol solvation Gibbs energies with respect to experimental values from the “Minnesota Solvation Database” (MNSOL), yielding a root mean square error (RMSE) of 1.5 kcal mol⁻¹ for the best-performing 2-parameter wet model, while the optimal water model developed for the pK_a part of the SAMPL6 challenge is kept unchanged (RMSE 1.6 kcal mol⁻¹ for neutral compounds from a model trained on both neutral and ionic species). Applying these models to the blind prediction set yields a log P RMSE of less than 0.5 for our best model (2-parameters, wet). Further analysis of our results reveals that a single compound is responsible for most of the error, SM15, without which the RMSE drops to 0.2. Since this is the only compound in the challenge dataset with a hydroxyl group we investigate other alcohols for which Gibbs energy of solvation data for both water and n-octanol are available in the MNSOL database to demonstrate a systematic cause of error and to discuss strategies for improvement.

High pressure response of 1H NMR chemical shifts of purine nucleotides

The study of the pressure response by NMR spectroscopy provides information on the thermodynamics of conformational equilibria in proteins and nucleic acids. For obtaining a database for expected pressure effects on free nucleotides and nucleotides bound in macromolecular complexes, the pressure response of ¹H chemical shifts and J-coupling constants of the purine 5'-ribonucleotides AMP, ADP, ATP, GMP, GDP, and GTP were studied in the absence and presence of Mg²⁺-ions. Experiments are supported by quantum-chemical calculations of populations and chemical shift differences in order to corroborate structural interpretations and to estimate missing data for AMP. The preference of the ribose S puckering obtained from the analysis of the experimental J-couplings is also confirmed by the calculations. In addition, the pressure response of the non-hydrolysable GTP analogues GppNHp, GppCH₂p, and GTPγS was examined within a pressure range up to 200 MPa. As observed earlier for ³¹P NMR chemical shifts of these nucleotides the pressure dependence of chemical shifts is clearly non-linear in most cases. In di- and tri-phospho nucleosides, the resonances of the two protons bound to the ribose 5' carbon are non-equivalent and can be observed separately. The gg-rotamer at C4'- C5' bond is strongly preferred and the downfield shifted resonance can be assigned to the H5″ proton in the nucleotides. In contrast, in adenosine itself the frequencies of the two resonances are interchanged.

pKa calculations for tautomerizable and conformationally flexible molecules: partition function vs. state transition approach

Calculations of acidities of molecules with multiple tautomeric and/or conformational states require adequate treatment of the relative energetics of accessible states accompanied by a statistical-mechanical formulation of their contribution to the macroscopic pK_a value. Here, we demonstrate rigorously the formal equivalence of two such approaches: a partition function treatment and statistics over transitions between molecular tautomeric and conformational states in the limit of a theory that does not require adjustment by empirical parameters correcting energetic values. However, for a frequently employed correction scheme, linear scaling of (free) energies and regression with respect to reference data taking an additive constant into account, this equivalence breaks down if more than one acid or base state is involved. The consequences of the resulting inconsistency are discussed on our datasets developed for aqueous pK_a predictions during the recent SAMPL6 challenge, where molecular state energetics were computed based on the "embedded cluster reference interaction site model" (EC-RISM). This method couples integral equation theory as a solvation model to quantum-chemical calculations and yielded a test set root mean square error of 1.1 pK units from a partition function ansatz. For all practical purposes, the present results indicate that a state transition approach yields comparable accuracy despite the formal theoretical inconsistency, and that an additive regression intercept, which is strictly constant in the limit of large compound mass only, is a valid approximation. Graphical abstract Embedded cluster reference interaction site model-derived vs. experimental pK_a for the test set calculated with either the partition function (blue) or the state transition approach (red), using m as a free parameter.

[Estimation of free cortisol in urine (author's transl)]

Radio-immunological estimation of free cortisol in urine in a group of normal persons showed a median daily excretion of 42,5 +/- 22,8 mug (n = 31) using radio-loss-labelling, extraction, and thin layer chromatography. After simple extraction the median value was 54.1 +/- 15.3 mug/d (n = 10). Radioimmunological estimation in simple urine extract seems to be sufficient for routine clinical purposes. The estimation of free cortisol in urine permits a better differentiation between normal adrenal function and adrenal hyperactivity than the determination of 17-hydroxycorticosteroids (17-OHCS). This is shown by comparison of 1) the frequency distribution of the basal values of 17-OHCS and free cortisol in normal persons and in Cushing's disease, 2) the suppression effect of dexamethasone on 17-OHCS and free cortisol in Cushing patients, and 3) the stimulatory effect of ACTH on 17-OHCS and free cortisol in normal persons.

[Diurnal variation in the urinary excretion of vanilmandelic acid and aldosterone in normal and hypertensive persons (author's transl)]

The excretion of creatinine, VMA, and aldosterone was measured in the urine of 14 normal and 25 persons with stable essential hypertension during a day time (7.00-23.00 hrs) and a night time (23.00-7.00 hrs) period. Excretion of VMA and aldosterone per hour in the two groups was lower during the night period than during the day. No significant difference existed between the excretion values of both groups. Night/day rations for the excretion of VMA and of aldosterone respectively, per mg of creatinine excreted were calculated for each subject. In the normal group the night/day ratios of aldosterone excretion showed a significant correlation to the corresponding ratios of VMA excretion; such correlation was not found in the hypertensive group. The fact that a correlation was observed in the normal group is compatible with the hypothesis that catecholamines are involved in the regulation of the renin-angiotensin-aldosterone system. Lack of correlation in the group of hypertensive does not as yet allow any conclusions. In two patients with primary hyperaldosteronism, the aldosterone excretion during the night was relatively more elevated than in the day.