Difluoro-dioxolo-benzoimidazol-benzamides as potent inhibitors of CK1δ and ε with nanomolar inhibitory activity on cancer cell proliferation.

Deregulation of CK1 (casein kinase 1) activity can be involved in the development of several pathological disorders and diseases such as cancer. Therefore, research interest in identifying potent CK1-specific inhibitors is still increasing. A previously published potent and selective benzimidazole-derived CK1δ/ε-specific inhibitor compound with significant effects on several tumor cell lines was further modified to difluoro-dioxolo-benzoimidazole derivatives displaying remarkable inhibitory effects and increased intracellular availability. In the present study, we identified two heterocyclic molecules as new CK1-specific inhibitor compounds with favorable physicochemical properties and notable selectivity in a kinome-wide screen. Being compared to other CK1 isoforms, these compounds exhibited advanced isoform selectivity toward CK1δ. Moreover, newly designed compounds showed increased growth inhibitory activity in a panel of different tumor cell lines as determined by analyses of cell viability and cell cycle distribution. In summary, presented lead optimization resulted in new highly selective CK1δ-specific small molecule inhibitors with increased biological activity.

2-Benzamido-N-(1H-benzo[d]imidazol-2-yl)thiazole-4-carboxamide derivatives as potent inhibitors of CK1δ/ε.

In this study we identified two heterocyclic compounds (5 and 6) as potent and specific inhibitors of CK1δ (IC(50) = 0.040 and 0.042 µM, respectively). Whereas compound 5 exhibited fivefold higher affinity towards CK1δ than to CK1ε (IC(50) CK1ε = 0.199 µM), compound 6 also inhibited CK1ε (IC(50) = 0.0326 µM) in the same range as CK1δ. Selected compound 5 was screened over 442 kinases identifying 5 as a highly potent and selective inhibitor of CK1δ. X-ray analysis of 5 bound to CK1δ demonstrated its binding mode. In addition, characterization of 5 and 6 in a cell biological approach revealed the ability of both compounds to inhibit proliferation of tumor cell lines in a dose and cell line specific manner. In summary, our optimizations lead to the development of new highly selective CK1δ and ε specific inhibitors with biological activity.

Application of information theory to feature selection in protein docking.

In the era of structural genomics, the prediction of protein interactions using docking algorithms is an important goal. The success of this method critically relies on the identification of good docking solutions among a vast excess of false solutions. We have adapted the concept of mutual information (MI) from information theory to achieve a fast and quantitative screening of different structural features with respect to their ability to discriminate between physiological and nonphysiological protein interfaces. The strategy includes the discretization of each structural feature into distinct value ranges to optimize its mutual information. We have selected 11 structural features and two datasets to demonstrate that the MI is dimensionless and can be directly compared for diverse structural features and between datasets of different sizes. Conversion of the MI values into a simple scoring function revealed that those features with a higher MI are actually more powerful for the identification of good docking solutions. Thus, an MI-based approach allows the rapid screening of structural features with respect to their information content and should therefore be helpful for the design of improved scoring functions in future. In addition, the concept presented here may also be adapted to related areas that require feature selection for biomolecules or organic ligands.

Monitoring biological membrane-potential changes: a CI QM/MM study.

In recent decades, new less-invasive, nonlinear optical methods have been proposed and optimized for monitoring fast physiological processes in biological cells. One of these methods allows the action potential (AP) in cardiomyocytes or neurons to be monitored by means of second-harmonic generation (SHG). We now present the first, to our knowledge, simulations of the dependency of the intensity of the second harmonic (I(SHG)) on variations of the transmembrane potential (TMP) in a cardiomyocyte during an action potential (AP). For this, an amphiphilic potential-sensitive styryl dye molecule with nonlinear optical properties was embedded in a dipalmitoylphosphatidylcholine (DPPC) bilayer, replacing one of the phospholipid molecules. External electrical fields with different strengths were applied across the membrane to simulate the AP of a heart-muscle cell. We used a combined classical/quantum mechanical approach to model the structure and the spectroscopic properties of the embedded chromophore. Two 10 ns molecular dynamics (MD) simulations provided input geometries for semiempirical molecular orbital (QM/MM) single-point configuration interaction (CI) calculations, which were used to calculate the wavelengths and oscillator strengths of electronic transitions in the di-8-ANEPPS dye molecule. The results were then used in a sum-over-states treatment to calculate the second-order hyperpolarizability. The square of the hyperpolarizability scales with the intensity of the second harmonic, which is used to monitor the action potentials of cardiomyocytes experimentally. Thus, we computed changes in the intensity of the second harmonic (DeltaI(SHG)) as function of TMP changes. Our results agree well with experimental measurements.

Molecular dynamics characterization of the structures and induction mechanisms of a reverse phenotype of the tetracycline receptor.

Molecular-dynamics simulations have been used to investigate the mechanism of induction of a mutant (revTetR) of the tetracycline repressor protein (TetR) that shows the reverse phenotype (i.e., it is induced in the absence of tetracyclines and not in their presence). Low-frequency, normal-mode analyses demonstrate that the reverse phenotype is reproduced by the simulations on the basis of criteria established for wild-type TetR. The reverse phenotype is caused by the fact that the DNA-binding heads in revTetR are closer than the ideal distance needed for DNA-binding when no inducer is present. This distance increases on binding an inducer. Whereas this distance increase makes the interhead distance too large in wild-type TetR, it increases to the ideal value in revTetR. Thus, the mechanism of induction is the same for the two proteins, but the consequences are reversed because of the smaller interhead distance in revTetR when no inducer is present.

SCRF-DFT and NMR comparison of tetracycline and 5a,6-anhydrotetracycline in solution.

A combination of structures, energies, and spectral data calculated using density functional theory (DFT) with experimental NMR data has been used to assign conformational equilibria for tetracycline and 5a,6-anhydrotetracycline in water at pH 1, 7, and 10 and in chloroform (5a,6-anhydrotetracycline) and methanol (tetracycline). The results suggest that tetracycline always prefers the extended conformation but that 5a,6-anhydrotetracycline exists in water as a mixture of the two conformers and in chloroform exclusively in the twisted conformation. The conformational equilibria are also shown to be pH dependent.

Structural changes and binding characteristics of the tetracycline-repressor binding site on induction.

The binding motif (pharmacophore) for induction and the changes in the structure of the binding site that accompany induction have been determined from molecular-dynamics simulations on the tetracycline-repressor signal-transduction protein. The changes and the induction mechanism are discussed and compared with conclusions drawn from earlier X-ray structures. The differences in inducer strength of tetracycline and 5a,6-anhydrotetracycline are discussed with respect to their interaction in the MD simulations.

The structure of 5a,6-anhydrotetracycline and its Mg2+ complexes in aqueous solution.

Semiempirical molecular orbital theory has been used for a systematic scan of the binding positions for a Mg2+ ion with 5a,6-anhydrotetracycline taking both conformational flexibility and possible different tautomeric forms into account. The magnesium ion has been calculated alone and with four or five complexed water molecules in order to simulate the experimental situation more closely. The results are analyzed by comparing the behavior of the title compound with that of tetracycline itself and possible causes for the stronger induction of the Tetracycline Receptor (TetR) by 5a,6-anhydrotetracycline than by tetracycline are considered. Energetically favored 3D-structure of the zwitteranionic 5a,6-anhydrotetracycline magnesium complex in solution.

Molecular dynamics simulations of the tetracycline-repressor protein: the mechanism of induction.

Molecular dynamics simulations on the tetracycline-repressor (TetR) protein, both in the absence of an inducer and complexed with the inducers tetracycline and 5a,6-anhydrotetracycline, show significant differences in the structures and dynamics of the induced and non-induced forms of the protein. Calpha-density-difference plots, low-frequency normal vibrations and inter-residue interaction energies all point to a common mechanism of induction. The inducer displaces Asp156 from the magnesium ion in the binding pocket, leading to a short cascade of rearrangements of salt bridges that results in the allosteric change. The increased flexibility of the induced form of the protein is suggested to contribute to the decrease in binding affinity to DNA on induction.

Simulating FRET from tryptophan: is the rotamer model correct?

We present a computational model study designed to simulate the results of time-resolved fluorescence spectra of tryptophan in proteins. In such measurements, the occurrence of more than one fluorescence lifetime is generally attributed to the existence of several tryptophan rotamers and/or structural conformations of the protein structure. The protein system we chose for this initial study is the tetracycline repressor (TetR), an interesting model system for the investigation of the mechanisms of transcriptional regulation. Fluorescence resonance energy transfer (FRET) from tryptophan to tetracycline is frequently observed in complexes of the TetR with the antibiotic tetracycline. We use a combined classical/quantum mechanical approach to model the structure and the spectroscopic properties of the TetR-tetracycline complex. A classical molecular dynamics simulation provides input geometries for semiempirical quantum mechanical/molecular mechanical (QM/MM) single-point configuration interaction (CI) calculations, which are used to calculate tryptophan vertical absorption and fluorescence energies and intensities as well as relative FRET rate constants. These rate constants together with the Einstein coefficients for spontaneous emission and an assumed rate for nonradiative deactivation allow us to simulate fluorescence decay curves with and without FRET and for the entire ensemble as well as for individual rotamers. Our results indicate that the classical "rotamer model", used to explain the multiexponential fluorescence-decay curves of time-resolved tryptophan emission spectra, can be extended to systems with FRET acceptors present in the protein matrix but that the interpretation of the fitted lifetimes is different to that usually used.

Systematic surface scan of the most favorable interaction sites of magnesium ions with tetracycline.

AM1 semiempirical molecular orbital calculations have been used to probe the complexation sites for naked and hydrated magnesium ions to the different conformations and protonation states of tetracycline. The calculations reveal a wealth of possible magnesium complexation sites within a small energy range, but also indicate that magnesium complexation does not change the conformational behavior of tetracycline significantly. A hitherto unknown solvated conformation is suggested for deprotonated tetracycline.