Discovery of CDK8/CycC Ligands with a New Virtual Screening Tool

The General Atomic and Molecular Electronic Structure System (GAMESS): Novel Methods on Novel Architectures

The primary focus of GAMESS over the last 5 years has been the development of new high-performance codes that are able to take effective and efficient advantage of the most advanced computer architectures, both CPU and accelerators. These efforts include employing density fitting and fragmentation methods to reduce the high scaling of well-correlated (e.g., coupled-cluster) methods as well as developing novel codes that can take optimal advantage of graphical processing units and other modern accelerators. Because accurate wave functions can be very complex, an important new functionality in GAMESS is the quasi-atomic orbital analysis, an unbiased approach to the understanding of covalent bonds embedded in the wave function. Best practices for the maintenance and distribution of GAMESS are also discussed.

Uncovering water effects in protein-ligand recognition: importance in the second hydration shell and binding kinetics

Developing a ligand with high affinity for a specific protein target is essential for drug design, and water molecules are well known to play a key role in protein-drug recognition. However, predicting the role of particularly ordered water molecules in drug binding remains challenging. Furthermore, hydration free energy contributed from the water network, including the second shell of water molecules, is far from being well studied. In this research we focused on these aspects to accurately and efficiently evaluate water effects in protein-ligand binding affinity. We developed a new strategy using a free-energy calculation method, VM2. We successfully predicted the stable ordered water molecules in a number of protein systems: PDE 10a, HSP90, tryptophan synthase (TRPS), CDK2 and Factor Xa. In some of these, the second shell of water molecules appeared to be critical in protein-ligand binding. We also applied the strategy to largely improve binding free energy calculation using the MM/PBSA method. When applying MM/PBSA alone for two systems, CDK2 and Factor Xa, the computed binding free energy resulted in poor to moderate R² values with experimental data. However, including water free energy correction greatly improved the free energy calculation. Furthermore, our work helped to explain how xk263 is a 1000 times faster binder to HIVp than ritonavir, a potentially useful tool for investigating binding kinetics. Our studies reveal the importance of fully considering water effects in therapeutic developments in pharmaceutical and biotechnology industries and for fundamental research in protein-ligand recognition.

Capture of activated dioxygen intermediates at the copper-active site of a lytic polysaccharide monooxygenase

Metalloproteins perform a diverse array of redox-related reactions facilitated by the increased chemical functionality afforded by their metallocofactors. Lytic polysaccharide monooxygenases (LPMOs) are a class of copper-dependent enzymes that are responsible for the breakdown of recalcitrant polysaccharides via oxidative cleavage at the glycosidic bond. The activated copper-oxygen intermediates and their mechanism of formation remains to be established. Neutron protein crystallography which permits direct visualization of protonation states was used to investigate the initial steps of oxygen activation directly following active site copper reduction in Neurospora crassa LPMO9D. Herein, we cryo-trap an activated dioxygen intermediate in a mixture of superoxo and hydroperoxo states, and we identify the conserved second coordination shell residue His157 as the proton donor. Density functional theory calculations indicate that both superoxo and hydroperoxo active site states are stable. The hydroperoxo formed is potentially an early LPMO catalytic reaction intermediate or the first step in the mechanism of hydrogen peroxide formation in the absence of substrate. We observe that the N-terminal amino group of the copper coordinating His1 remains doubly protonated directly following molecular oxygen reduction by copper. Aided by molecular dynamics and mining minima free energy calculations we establish that the conserved second-shell His161 in MtPMO3* displays conformational flexibility in solution and that this flexibility is also observed, though to a lesser extent, in His157 of NcLPMO9D. The imidazolate form of His157 observed in our structure following oxygen intermediate protonation can be attributed to abolished His157 flexibility due steric hindrance in the crystal as well as the solvent-occluded active site environment due to crystal packing. A neutron crystal structure of NcLPMO9D at low pH further supports occlusion of the active site since His157 remains singly protonated even at acidic conditions.

Computation of host-guest binding free energies with a new quantum mechanics based mining minima algorithm

A new method called QM-VM2 is presented that efficiently combines statistical mechanics with quantum mechanical (QM) energy potentials in order to calculate noncovalent binding free energies of host-guest systems. QM-VM2 efficiently couples the use of semi-empirical QM (SEQM) energies and geometry optimizations with an underlying molecular mechanics (MM) based conformational search, to find low SEQM energy minima, and allows for processing of these minima at higher levels of ab initio QM theory. A progressive geometry optimization scheme is introduced as a means to increase conformational sampling efficiency. The newly implemented QM-VM2 is used to compute the binding free energies of the host molecule cucurbit[7]uril and a set of 15 guest molecules. The results are presented along with comparisons to experimentally determined binding affinities. For the full set of 15 host-guest complexes, which have a range of formal charges from +1 to +3, SEQM-VM2 based binding free energies show poor correlation with experiment, whereas for the ten +1 complexes only, a significant correlation (R² = 0.8) is achieved. SEQM-VM2 generation of conformers followed by single-point ab initio QM calculations at the dispersion corrected restricted Hartree-Fock-D3(BJ) and TPSS-D3(BJ) levels of theory, as post-processing corrections, yields a reasonable correlation with experiment for the full set of host-guest complexes (R² = 0.6 and R² = 0.7, respectively) and an excellent correlation for the +1 formal charge set (R² = 1.0 and R² = 0.9, respectively), as long as a sufficiently large basis set (triple-zeta quality) is employed. The importance of the inclusion of configurational entropy, even at the MM level, for the achievement of good correlation with experiment was demonstrated by comparing the calculated ΔE values with experiment and finding a considerably poorer correlation with experiment than for the calculated free energy ΔE - TΔS. For the complete set of host-guest systems with the range of formal charges, it was observed that the deviation of the predicted binding free energy from experiment correlates somewhat with the net charge of the systems. This observation leads to a simple empirical interpolation scheme to improve the linear regression of the full set.

Super-Enhancers in the Regulation of Gene Transcription: General Aspects and Antitumor Targets

Super-enhancers (genome elements that activate gene transcription) are DNA regions with an elevated concentration of transcriptional complexes. These multiprotein structures contain, among other components, the cyclin-dependent kinases 8 and 19. These and other transcriptional protein kinases are regarded as novel targets for pharmacological inhibition by antitumor drug candidates.

Pyrazolyl-Ureas as Interesting Scaffold in Medicinal Chemistry

The pyrazole nucleus has long been known as a privileged scaffold in the synthesis of biologically active compounds. Within the numerous pyrazole derivatives developed as potential drugs, this review is focused on molecules characterized by a urea function directly linked to the pyrazole nucleus in a different position. In the last 20 years, the interest of numerous researchers has been especially attracted by pyrazolyl-ureas showing a wide spectrum of biological activities, ranging from the antipathogenic activities (bacteria, plasmodium, toxoplasma, and others) to the anticarcinogenic activities. In particular, in the anticancer field, pyrazolyl-ureas have been shown to interact at the intracellular level on many pathways, in particular on different kinases such as Src, p38-MAPK, TrKa, and others. In addition, some of them evidenced an antiangiogenic potential that deserves to be explored. This review therefore summarizes all these biological data (from 2000 to date), including patented compounds.

Synthesis and biological evaluation of small molecule modulators of CDK8/Cyclin C complex with phenylaminoquinoline scaffold

Background

CDK8/CycC complex has kinase activity towards the carboxyterminal domain of RNA polymerase II, and contributes to the regulation of transcription via association with the mediator complex. Different human malignancies, mainly colorectal and gastric cancers, were produced as a result of overexpression of CDK8/CycC in the mediator complex. Therefore, CDK8/CycC complex represents as a cancer oncogene and it has become a potential target for developing CDK8/CycC modulators.

Methods

A series of nine 4-phenylaminoquinoline scaffold-based compounds 5a-i was synthesized, and biologically evaluated as potential CDK8/CycC complex inhibitors.

Results

The scaffold substituent effects on the intrinsic inhibitory activity toward CDK8/CycC complex are addressed trying to present a novel outlook of CDK8/CycC Complex inhibitors with 4-phenylaminoquinoline scaffold in cancer therapy. The secondary benzenesulfonamide analogues proved to be the most potent compounds in suppressing CDK8/CycC enzyme, whereas, their primary benzenesulfonamide analogues showed inferior activity. Moreover, the benzene reversed sulfonamide analogues were totally inactive.

Discussion

The titled scaffold showed promising inhibitory activity data and there is a crucial role of un/substituted sulfonamido group for CDK8/CycC complex inhibitory activity. Compound 5d showed submicromolar potency against CDK8/CycC (IC₅₀ = 0.639 µM) and it can be used for further investigations and to design another larger library of phenylaminoquinoline scaffold-based analogues in order to establish detailed SARs.

Transient States and Barriers from Molecular Simulations and the Milestoning Theory: Kinetics in Ligand-Protein Recognition and Compound Design

This study presents a novel computational approach to study molecular recognition and binding kinetics for drug-like compounds dissociating from a flexible protein system. The intermediates and their free energy profile during ligand association and dissociation processes control ligand-protein binding kinetics and bring a more complete picture of ligand-protein binding. The method applied the milestoning theory to extract kinetics and thermodynamics information from running short classical molecular dynamics (MD) simulations for frames from a given dissociation path. High-dimensional ligand-protein motions (3N-6 degrees of freedom) during ligand dissociation were reduced by use of principal component modes for assigning more than 100 milestones, and classical MD runs were allowed to travel multiple milestones to efficiently obtain ensemble distribution of initial structures for MD simulations and estimate the transition time and rate during ligand traveling between milestones. We used five pyrazolourea ligands and cyclin-dependent kinase 8 with cyclin C (CDK8/CycC) as our model system as well as metadynamics and a pathway search method to sample dissociation pathways. With our strategy, we constructed the free energy profile for highly mobile biomolecular systems. The computed binding free energy and residence time correctly ranked the pyrazolourea ligand series, in agreement with experimental data. Guided by a barrier of a ligand passing an αC helix and activation loop, we introduced one hydroxyl group to parent compounds to design our ligands with increased residence time and validated our prediction by experiments. This work provides a novel and robust approach to investigate dissociation kinetics of large and flexible systems for understanding unbinding mechanisms and designing new small-molecule drugs with desired binding kinetics.