RORγT agonists as immune modulators in anticancer therapy

RORγT is a transcription factor that directs the development of Th17 lymphocytes and other IL-17-expressing cells (e.g., Tc17 and ILC3 cells). These cells are involved in the body's defense against pathogenic bacteria and fungi, but they also participate in maintaining the proinflammatory environment in some autoimmune diseases and play a role in the immune system's response to cancer. Similar to other members of the nuclear receptor superfamily, the activity of RORγT is regulated by low-molecular-weight ligands. Therefore, extensive efforts have been dedicated to identifying inverse agonists that diminish the activity of this receptor and subsequently inhibit the development of autoimmune diseases. Unfortunately, in the pursuit of an ideal inverse agonist, the development of agonists has been overlooked. It is important to remember that these types of compounds, by stimulating lymphocytes expressing RORγT (Th17 and Tc17), can enhance the immune system's response to tumors. In this review, we present recent advancements in the biology of RORγT agonists and their potential application in anticancer therapy.

The application of machine learning methods to the prediction of novel ligands for RORγ/RORγT receptors

In this work, we developed and applied a computational procedure for creating and validating predictive models capable of estimating the biological activity of ligands. The combination of modern machine learning methods, experimental data, and the appropriate setup of molecular descriptors led to a set of well-performing models. We thoroughly inspected both the methodological space and various possibilities for creating a chemical feature space. The resulting models were applied to the virtual screening of the ZINC20 database to identify new, biologically active ligands of RORγ receptors, which are a subfamily of nuclear receptors. Based on the known ligands of RORγ, we selected candidates and calculate their predicted activities with the best-performing models. We chose two candidates that were experimentally verified. One of these candidates was confirmed to induce the biological activity of the RORγ receptors, which we consider proof of the efficacy of the proposed methodology.

Anti-hepatocellular carcinoma activity of the cyclin-dependent kinase inhibitor AT7519

Hepatocellular carcinoma (HCC) is one of the most common cancerous tumors and one of the leading causes of death among cancer-related disorders. Chemotherapy is ineffective in HCC patients, and the number of drugs that are in use is limited. Thus, new molecules are needed that could increase the effectiveness of anti-HCC regimens. Here, we show that AT7519, a CDK inhibitor, exerts positive effects on HCC cells: it inhibits proliferation, migration and clonogenicity. Detailed analysis of the transcriptomes of cells treated with this compound indicated that AT7519 affects a substantial portion of genes that are associated with HCC development and progression. Moreover, we showed that the concomitant use of AT7519 with gefitinib or cabozantinib sensitized HCC cells to these drugs. Thus, our research indicates that AT7519 is worth considering in monotherapy for hepatocellular carcinoma patients or in combination with other drugs, e.g., gefitinib or cabozantinib.

Identification of Corosolic and Oleanolic Acids as Molecules Antagonizing the Human RORγT Nuclear Receptor Using the Calculated Fingerprints of the Molecular Similarity

RORγT is a protein product of the RORC gene belonging to the nuclear receptor subfamily of retinoic-acid-receptor-related orphan receptors (RORs). RORγT is preferentially expressed in Th17 lymphocytes and drives their differentiation from naive CD4+ cells and is involved in the regulation of the expression of numerous Th17-specific cytokines, such as IL-17. Because Th17 cells are implicated in the pathology of autoimmune diseases (e.g., psoriasis, inflammatory bowel disease, multiple sclerosis), RORγT, whose activity is regulated by ligands, has been recognized as a drug target in potential therapies against these diseases. The identification of such ligands is time-consuming and usually requires the screening of chemical libraries. Herein, using a Tanimoto similarity search, we found corosolic acid and other pentacyclic tritepenes in the library we previously screened as compounds highly similar to the RORγT inverse agonist ursolic acid. Furthermore, using gene reporter assays and Th17 lymphocytes, we distinguished compounds that exert stronger biological effects (ursolic, corosolic, and oleanolic acid) from those that are ineffective (asiatic and maslinic acids), providing evidence that such combinatorial methodology (in silico and experimental) might help wet screenings to achieve more accurate results, eliminating false negatives.

Digoxin, an Overlooked Agonist of RORγ/RORγT

Digoxin was one of the first identified RORγT receptor inverse agonists inhibiting the differentiation of Th17 cells. However, this compound exhibits inhibitory activity at relatively high concentrations that mediate cytotoxic effects. We previously identified several cardenolides that are structurally similar to digoxin that were able to induce RORγ/RORγT-dependent transcription. These observations encouraged us to reanalyze the effects of digoxin on RORγ/RORγT-dependent transcription at low, noncytotoxic concentrations. Digoxin induced RORγ/RORγT-dependent transcription in HepG2 and Th17 cells. Furthermore, analysis of the transcriptomes of Th17 cells cultured in the presence of digoxin revealed the induction of the expression of numerous Th17-specific genes, including IL17A/F, IL21, IL22, IL23R, CCR4, and CCR6. Thus, our study, which includes data obtained from intact cells, indicates that digoxin, similar to other cardenolides, is a potent RORγ/RORγT receptor activator and that its structure may serve as a starting point for the design of dedicated molecules that can be used in the development of adoptive cell therapy (ACT).

The MP2-F12 method in the Turbomole program package

A detailed description of the explicitly correlated second-order Møller-Plesset perturbation theory (MP2-F12) method, as implemented in the TURBOMOLE program package, is presented. The TURBOMOLE implementation makes use of density fitting, which greatly reduces the prefactor for integral evaluation. Methods are available for the treatment of ground states of open- and closed-shell species, using unrestricted as well as restricted (open-shell) Hartree-Fock reference determinants. Various methodological choices and approximations are discussed. The performance of the TURBOMOLE implementation is illustrated by example calculations of the molecules leflunomide, prednisone, methotrexate, ethylenedioxytetrafulvalene, and a cluster model for the adsorption of methanol on the zeolite H-ZSM-5. Various basis sets are used, including the correlation-consistent basis sets specially optimized for explicitly correlated calculations (cc-pVXZ-F12).

Stabilization of very rare tautomers of uracil by an excess electron

We characterized valence-type and dipole-bound anionic states of uracil using various electronic structure methods. We found that the most stable anion is related to neither the canonical 2,4-dioxo nor a rare imino-hydroxy tautomer. Instead, it is related to an imino-oxo tautomer, in which the N1H proton is transferred to the C5 atom. This valence anion is characterized by an electron vertical detachment energy (VDE) of 1267 meV and it is adiabatically stable with respect to the canonical neutral by 3.93 kcal mol(-1). It is also more stable by 2.32 and 5.10 kcal mol(-1) than the dipole-bound and valence anion, respectively, of the canonical tautomer. The VDE values for the former and the latter are 73 and 506 meV, respectively. Another, anionic, low-lying imino-oxo tautomer with a VDE of 2499 meV has a proton transferred from N3H to C5. It is less stable than the neutral canonical tautomer by 1.38 kcal mol(-1). The mechanism of formation of anionic tautomers with the carbon C5 protonated may involve intermolecular proton transfer or dissociative electron attachment to the canonical neutral tautomer followed by a barrier-free attachment of a hydrogen atom to C5. The six-member ring structure of anionic tautomers with carbon atoms protonated might be unstable upon an excess electron detachment. Indeed, the neutral systems resulting from electron detachment from anionic tautomers with carbon atoms protonated evolve along barrier-free decomposition pathways to a linear or a bicyclo structure, which might be viewed as lesions to RNA. Within the PCM hydration model, the low-lying valence anions become adiabatically bound with respect to the canonical neutral and the two most stable tautomers have carbon atoms protonated.

Rovibrational energy levels of H3(+) with energies above the barrier to linearity

The H(3)(+) potential energy surface is sampled at 5900 geometries with the emphasis on the nonequilibrium and asymptotic points. Apart from the Born-Oppenheimer energy converged to the accuracy better than 0.02 cm(-1), the adiabatic and the leading relativistic corrections are computed at each geometry. To represent analytically the potential energy surface, the parameters of a power series are determined by fitting to the computed energy points. Possible choice of nuclear masses simulating the nonadiabatic effects in solving the nuclear Schrodinger equation is analyzed. A set of theoretically predicted rovibrational transitions is confronted with the experimental data in the 10,700-13,700 cm(-1) window of the spectra.

Electrostatic potential maps of damaged DNA studied by image analysis tools. 8-Oxoguanine and abasic site lesions

Changes of electrostatic potential around the DNA molecule resulting from chemical modifications of nucleotides may play a role in enzymatic recognition of damaged sites. The electrostatic potential around the DNA fragments containing either the intact guanine-cytosine pair or 8-oxoguanine-cytosine or the guanine-abasic site was projected on a cylindrical surface around the double helix. The 2D maps of EP of intact and damaged DNA fragments were compared using image analysis methods. Occurrence of abasic site and 8-oxoguanine lesions were found to be reflected in the EP maps. In the case of the 8-oxoguanine lesion, the two phosphate groups and countercations of the damaged strand are moved away from the lesion in opposite directions, whereas they are moved in the same direction in the case of the abasic site lesion. The characteristic features of 8-oxoguanine lesion might be identified in the major groove, whereas the features of abasic site lesion the minor groove.

Photoelectron spectrum of valence anions of uracil and first-principles calculations of excess electron binding energies

The photoelectron spectrum (PES) of the uracil anion is reported and discussed from the perspective of quantum chemical calculations of the vertical detachment energies (VDEs) of the anions of various tautomers of uracil. The PES peak maximum is found at an electron binding energy of 2.4 eV, and the width of the main feature suggests that the parent anions are in a valence rather than a dipole-bound state. The canonical tautomer as well as four tautomers that result from proton transfer from an NH group to a C atom were investigated computationally. At the Hartree-Fock and second-order Moller-Plesset perturbation theory levels, the adiabatic electron affinity (AEA) and the VDE have been converged to the limit of a complete basis set to within +/-1 meV. Post-MP2 electron-correlation effects have been determined at the coupled-cluster level of theory including single, double, and noniterative triple excitations. The quantum chemical calculations suggest that the most stable valence anion of uracil is the anion of a tautomer that results from a proton transfer from N1H to C5. It is characterized by an AEA of 135 meV and a VDE of 1.38 eV. The peak maximum is as much as 1 eV larger, however, and the photoelectron intensity is only very weak at 1.38 eV. The PES does not lend support either to the valence anion of the canonical tautomer, which is the second most stable anion, and whose VDE is computed at about 0.60 eV. Agreement between the peak maximum and the computed VDE is only found for the third most stable tautomer, which shows an AEA of approximately -0.1 eV and a VDE of 2.58 eV. This tautomer results from a proton transfer from N3H to C5. The results illustrate that the characteristics of biomolecular anions are highly dependent on their tautomeric form. If indeed the third most stable anion is observed in the experiment, then it remains an open question why and how this species is formed under the given conditions.

Coupled-cluster and explicitly correlated perturbation-theory calculations of the uracil anion

A valence-type anion of the canonical tautomer of uracil has been characterized using explicitly correlated second-order Moller-Plesset perturbation theory (RI-MP2-R12) in conjunction with conventional coupled-cluster theory with single, double, and perturbative triple excitations. At this level of electron-correlation treatment and after inclusion of a zero-point vibrational energy correction, determined in the harmonic approximation at the RI-MP2 level of theory, the valence anion is adiabatically stable with respect to the neutral molecule by 40 meV. The anion is characterized by a vertical detachment energy of 0.60 eV. To obtain accurate estimates of the vertical and adiabatic electron binding energies, a scheme was applied in which electronic energy contributions from various levels of theory were added, each of them extrapolated to the corresponding basis-set limit. The MP2 basis-set limits were also evaluated using an explicitly correlated approach, and the results of these calculations are in agreement with the extrapolated values. A remarkable feature of the valence anionic state is that the adiabatic electron binding energy is positive but smaller than the adiabatic electron binding energy of the dipole-bound state.

On the Unusual Stability of Valence Anions of Thymine Based on Very Rare Tautomers: A Computational Study

We characterized anionic states of thymine using various electronic structure methods, with the most accurate results obtained at the CCSD(T)/aug-cc-pVDZ level of theory followed by extrapolations to complete basis set limits. We found that the most stable anion in the gas phase is related to an imino-oxo tautomer, in which the N1H proton is transferred to the C5 atom. This valence anion, aT(c5)(nl), is characterized by an electron vertical detachment energy (VDE) of 1251 meV and it is adiabatically stable with respect to the canonical neutral nT(can) by 2.4 kcal/mol. It is also more stable than the dipole-bound (aT(dbs)(can)), and valence anion aT(val)(can) of the canonical tautomer. The VDE values for aT(dbs)(can)and T(val)(can) are 55 and 457 meV, respectively. Another, anionic, low-lying imino-oxo tautomer with a VDE of 2458 meV has a proton transferred from N3H to C5 aT(c5)(n3). It is less stable than aT(val)(can) by 3.3 kcal/mol. The mechanism of formation of anionic tautomers with the carbons C5 or C6 protonated may involve intermolecular proton transfer or dissociative electron attachment to the canonical neutral tautomer followed by a barrier-free attachment of a hydrogen atom to C5. The six-member ring structure of the anionic tautomers with carbon atoms protonated is unstable upon an excess electron detachment. Within the PCM hydration model, the low-lying valence anions become adiabatically bound with respect to the canonical neutral; becomes the most stable, being followed by aT(c5)(nl), aT(c5)(n3), aT(can), and aT(c5)(nl).

Anion of the formic acid dimer as a model for intermolecular proton transfer induced by a π* excess electron

The neutral and anionic formic acid dimers have been studied at the second-order Moller-Plesset and coupled-cluster level of theory with single, double, and perturbative triple excitations with augmented, correlation-consistent basis sets of double- and triple-zeta quality. Scans of the potential-energy surface for the anion were performed at the density-functional level of theory with a hybrid B3LYP functional and a high-quality basis set. Our main finding is that the formic acid dimer is susceptible to intermolecular proton transfer upon an excess electron attachment. The unpaired electron occupies a pi(*) orbital, the molecular moiety that accommodates an excess electron "buckles," and a proton is transferred to the unit where the excess electron is localized. As a consequence of these geometrical transformations, the electron vertical detachment energy becomes substantial, 2.35 eV. The anion is barely adiabatically unstable with respect to the neutral at 0 K. However, at standard conditions and in terms of Gibbs free energy, the anion is more stable than the neutral by +37 meV. The neutral and anionic dimers display different IR characteristics. In summary, the formic acid dimer can exist in two quasidegenerate states (neutral and anionic), which can be viewed as "zero" and "one" in the binary system. These two states are switchable and distinguishable.