Search

Recent advances in targeted protein degradation have enabled chemical hijacking of the ubiquitin-proteasome system to treat disease. The catalytic rate of cereblon (CRBN)-dependent bifunctional degradation activating compounds (BiDAC), which recruit CRBN to a chosen target protein, resulting in its ubiquitination and proteasomal degradation, is an important parameter to consider during the drug discovery process. In this work, an in vitro system was developed to measure the kinetics of BRD4 bromodomain 1 (BD1) ubiquitination by fitting an essential activator kinetic model to these data. The affinities between BiDACs, BD1, and CRBN in the binary complex, ternary complex, and full ubiquitination complex were characterized. Together, this work provides a new tool for understanding and optimizing the catalytic and thermodynamic properties of BiDACs.

With early life likely to have existed in a hot environment, enzymes had to cope with an inherent drop in catalytic speed caused by lowered temperature. Here we characterize the molecular mechanisms underlying thermoadaptation of enzyme catalysis in adenylate kinase using ancestral sequence reconstruction spanning 3 billion years of evolution. We show that evolution solved the enzyme’s key kinetic obstacle—how to maintain catalytic speed on a cooler Earth—by exploiting transition-state heat capacity. Tracing the evolution of enzyme activity and stability from the hot-start toward modern hyperthermophilic, mesophilic, and psychrophilic organisms illustrates active pressure versus passive drift in evolution on a molecular level, refutes the debated activity/stability trade-off, and suggests that the catalytic speed of adenylate kinase is an evolutionary driver for organismal fitness.

Macromolecular function is rooted in energy landscapes, where sequence determines not a single structure but an ensemble of conformations. Hence, evolution modifies a protein’s function by altering its energy landscape. Here, we recreate the evolutionary pathway between two modern human oncogenes, Src and Abl, by reconstructing their common ancestors. Our evolutionary reconstruction combined with x-ray structures of the common ancestor and pre–steady-state kinetics reveals a detailed atomistic mechanism for selectivity of the successful cancer drug Gleevec. Gleevec affinity is gained during the evolutionary trajectory toward Abl and lost toward Src, primarily by shifting an induced-fit equilibrium that is also disrupted in the clinical T315I resistance mutation. This work reveals the mechanism of Gleevec specificity while offering insights into how energy landscapes evolve.

Kinases perform phosphoryl-transfer reactions in milliseconds; without enzymes, these reactions would take about 8,000 years under physiological conditions. Despite extensive studies, a comprehensive understanding of kinase energy landscapes, including both chemical and conformational steps, is lacking. Here we scrutinize the microscopic steps in the catalytic cycle of adenylate kinase, through a combination of NMR measurements during catalysis, pre-steady-state kinetics, molecular-dynamics simulations and crystallography of active complexes. We find that the Mg(2+) cofactor activates two distinct molecular events: phosphoryl transfer (>10(5)-fold) and lid opening (10(3)-fold). In contrast, mutation of an essential active site arginine decelerates phosphoryl transfer 10(3)-fold without substantially affecting lid opening. Our results highlight the importance of the entire energy landscape in catalysis and suggest that adenylate kinases have evolved to activate key processes simultaneously by precise placement of a single, charged and very abundant cofactor in a preorganized active site.

Protein kinases are obvious drug targets against cancer, owing to their central role in cellular regulation. Since the discovery of Gleevec, a potent and specific inhibitor of Abl kinase, as a highly successful cancer therapeutic, the ability of this drug to distinguish between Abl and other tyrosine kinases such as Src has been intensely investigated but without much success. Using NMR and fast kinetics, we establish a new model that solves this longstanding question of how the two tyrosine kinases adopt almost identical structures when bound to Gleevec but have vastly different affinities. We show that, in contrast to all other proposed models, the origin of Abl's high affinity lies predominantly in a conformational change after binding. An energy landscape providing tight affinity via an induced fit and binding plasticity via a conformational-selection mechanism is likely to be general for many inhibitors.

Sophisticated protein kinase networks, empowering complexity in higher organisms, are also drivers of devastating diseases such as cancer. Accordingly, these enzymes have become major drug targets of the twenty-first century. However, the holy grail of designing specific kinase inhibitors aimed at specific cancers has not been found. Can new approaches in cancer drug design help win the battle with this multi-faced and quickly evolving enemy? In this perspective we discuss new strategies and ideas that were born out of a recent breakthrough in understanding the molecular basis underlying the clinical success of the cancer drug Gleevec. An “old” method, stopped-flow kinetics, combined with old enzymes, the ancestors dating back up to about billion years, provides an unexpected outlook for future intelligent design of drugs.

Protein phosphorylation is an essential regulatory mechanism that affects all aspect of cellular life from division and growth to aging and death. Misregulation of the signaling cascades leads to severe detrimental effects, and in humans often associated with cancer and other diseases. Phosphorylation is performed by a class of protein called kinases. Activation and deactivation of kinases is normally under tight control and is regulated via different mechanisms that are incredibly complex. In this work we combine phylogenetic resurrection techniques with biophysical and chemical approaches to analyze the regulatory mechanisms of modern tyrosine oncokinases Src and Abl, their common ancestor and the common ancestors between several other families of tyrosin kinases. Our results show how the regulatory elements appeared and developed throughout the evolution enabling selective regulation of complex modern cascades.

Proteins populate a landscape of conformations that changes upon binding native ligands or drugs, altering protein function. Different conformations reveal distinct potential drug target sites, including allosteric sites distal to the active site of the protein. It has been historically difficult to identify allosteric modulators early in the screening process since the structural methods required to reveal binding sites are relatively low throughput.

Design of specific kinase inhibitors is an appealing approach for developing new anticancer treatments. However, only a few success stories have been reported to date. Here we demonstrate how the combination of old-fashioned and new biophysical tools together with recent advances in genomics and molecular evolution can aid in overcoming existing limitations.

We have engineered a mutant of Dictyostelium discoideum ( Dicty ) myosin II that contains the same fast-reacting “SH1” thiol as in muscle myosin, spin-labeled it, and performed electron paramagnetic resonance (EPR) to compare the structure of the force-generating region of the two myosins. Dicty myosin serves as a model system for muscle myosin because of greater ease of mutagenesis, expression, and crystallization. The catalytic domains of these myosins have nearly identical crystal structures in the apo state, but there are significant differences in ATPase kinetics, and there are no crystal structures of skeletal muscle myosin with bound nucleotides, so another structural technique is needed. Previous EPR studies, with a spin label attached to SH1 in muscle myosin, have resolved the key structural states of this region. Therefore, we have performed identical experiments on both myosins spin-labeled at equivalent sites. Spectra were identical for the two myosins in the apo and ADP-bound states. With bound ADP and phosphate analogs, ( i ) both proteins exhibit two resolved structural states (prepowerstroke, postpowerstroke) in a single biochemical state (defined by the bound nucleotide), and ( ii ) these structural states are essentially identical in the two myosins but ( iii ) are occupied to different extents as a function of the biochemical state. We conclude that ( i ) myosin structural and biochemical states do not have a one-to-one correspondence, and ( ii ) Dicty myosin can serve as a good analog for structural studies of muscle myosin only if differences in the coupling between biochemical and structural states are taken into account.

Nitroxide radicals are widely used as molecular probes in different fields of chemistry and biology. In this work, we describe pH-sensitive imidazoline- and imidazolidine-based nitroxides with pK values in the range 4.7-7.6 (2,2,3,4,5,5-hexamethylperhydroimidazol-1-oxyl, 4-amino-2,2,5,5-tetramethyl-2,5-dihydro-1H-imidazol-1-oxyl, 4-dimethylamino-2,2-diethyl-5,5-dimethyl-2,5-dihydro-1H-imidazol-1-oxyl, and 2,2-diethyl-5,5-dimethyl-4-pyrrolidyline-1-yl-2,5-dihydro-1H-imidazol-1-oxyl), which allow the pH-monitoring inside chloroplasts. We have demonstrated that EPR spectra of these spin-probes localized in the thylakoid lumen markedly change with the light-induced acidification of the thylakoid lumen in chloroplasts. Comparing EPR spectrum parameters of intrathylakoid spin-probes with relevant calibrating curves, we could estimate steady-state values of lumen pHin established during illumination of chloroplasts with continuous light. For isolated bean (Vicia faba) chloroplasts suspended in a medium with pHout=7.8, we found that pHin approximately 5.4-5.7 in the state of photosynthetic control, and pHin approximately 5.7-6.0 under photophosphorylation conditions. Thus, ATP synthesis occurs at a moderate acidification of the thylakoid lumen, corresponding to transthylakoid pH difference DeltapH approximately 1.8-2.1. These values of DeltapH are consistent with a point of view that under steady-state conditions the proton gradient DeltapH is the main contributor to the proton motive force driving the operation of ATP synthesis, provided that stoichiometric ratio H+/ATP is nor =4-4.7.

We have used transient time-resolved FRET (TR2FRET) to monitor the conformation of the relay helix in a myosin II functional mutant during the recovery stroke in real time. Myosin was perturbed with the F506A mutation (Dictyostelium discoideum sequence), located within the relay loop in the force-generating region. F506 is a highly conserved residue in myosin II and is a hypertrophic cardiomyopathy mutation site. Previous studies [Tsiavaliaris, EMBO Rep, 2002, 3(11), 1099] showed a significant effect of the F506A mutation on myosin function. Actin affinity in the presence of ATP was increased, and the mutant did not move actin filaments in in vitro motility assays. A small decrease in intrinsic fluorescence was observed upon addition of excess ATP, but ATP binding and hydrolysis were not affected by the mutation. It was proposed that the F506A disrupts the communication between the active site and the lever arm. We engineered a double-Cys myosin mutant (A639C:K498C) in the Cys-less background with the F506A functional mutation, and labeled the mutant with optical probes. We used TR-FRET to determine the interprobe distance, and TR2FRET measurements after rapid mixing with ATP revealed changes in the relay helix conformation during the recovery stroke in real time. The mutation induced significant disorder of the relay helix in the force-generating region, but myosin still produces a recovery stroke, changing the relay helix conformation from straight to bent. We conclude that (a) the relay helix is disordered in myosin functional mutant F506A, which demonstrates the importance of the relay loop - relay helix interaction in the relay helix stabilization, and (b) the relay helix is the major structural element in the force-generating region of myosin, responsible for communication from the active site to the converter domain and the lever arm.

We have used dipolar electron-electron resonance (DEER) and time-resolved fluorescence resonance energy transfer (TR-FRET) to investigate the role of the myosin relay helix in coupling between the active site and the force-generating region of myosin II. Two double-Cys Dictyostelium myosin constructs have been engineered, and the structure of the relay helix was monitored by measuring interprobe distances in MSL/MSL or IAEDANS/Dabcyl-labeled myosin. Experiments were performed on WT myosin and on F506A, a functional mutant that has close to normal enzymatic activity but completely lacks motor functions (e.g. unable to move actin filaments or support cell development). We found that the WT myosin relay helix adopts two distinct states (straight and bent), with the bent conformation populated when ATP and ADP.Pi analogs are bound at the active site. In contrast, binding of nucleotide analogs had very little effect on relay helix conformation in the F506A mutant. In addition, the width of the distance distribution was significantly larger in the F506A compared to WT myosin, indicating loss of structural organization. Our results demonstrate that the relay helix plays a key role in coupling of myosin ATPase and motor activities. Loss of functionality observed in F506A myosin can be explained by the disruption of the relay helix-relay loop interactions that normally stabilize well-defined conformations of the myosin force-generating region allowing it to switch between distinct structural states.

We have used DEER (double electron-electron resonance) and TR-FRET (time-resolved fluorescence resonance energy transfer) to study conformational changes within the myosin II relay helix. Major structural changes during the myosin II ATPase cycle take place in the force-generating domain. Crystal structures show that the converter domain, the relay helix, and SH1 helix have different conformations in the proposed pre- and post- powerstroke structural states of myosin. This dramatic structural change is detected in solution studies by intrinsic fluorescence [1, 2] and EPR of a spin label attached to SH1 [3]. In the present study, we focus on the relay helix as a crucial structural element involved in coupling between the force-generating domain and the nucleotide-binding pocket. Cysteine mutations were introduced into a Cys-lite construct of Dictyostelium discoideum (Dicty) myosin in the lower 50k domain (either D515C or A639C) and the C-terminal end of the relay helix (K498C). These constructs were selectively modified with either MSL/MSL or IAEDANS/DABCYL pairs, and the distance between probes was measured in different myosin conformations trapped with nucleotides or nucleotide analogs Two conformations of the relay helix (with distinct probe-to-probe distances, presumably corresponding to the “straight” and “bent“ states of the relay helix) were resolved. Observed distances were in good agreement with existing crystal structures, but at least two distinct structural states were present in certain biochemical states (e.g, with bound ADP.BeF3, ADP.Vi, ADP.AlF4). The mole fraction of the “bent” conformation was higher with post-hydrolysis analogs (ADP.Vi, ADP.AlF4) bound at the active site. Our results reveal structural rearrangements within a single subdomain of myosin and provide insights into the coupling between ATP binding and changes in the force-generating region.

SummaryMolecular recognition plays a central role in biology, and protein dynamics has been acknowledged to be important in this process. However, it is highly debated whether conformational changes happen before ligand binding to produce a binding-competent state (conformational selection) or are caused in response to ligand binding (induced fit). Proposals for both mechanisms in protein/protein recognition have been primarily based on structural arguments. However, the distinction between them is a question of the probabilities of going via these two opposing pathways. Here, we present a direct demonstration of exclusive conformational selection in protein/protein recognition by measuring the flux for rhodopsin kinase binding to its regulator recoverin, an important molecular recognition in the vision system. Using nuclear magnetic resonance (NMR) spectroscopy, stopped-flow kinetics, and isothermal titration calorimetry, we show that recoverin populates a minor conformation in solution that exposes a hydrophobic binding pocket responsible for binding rhodopsin kinase. Protein dynamics in free recoverin limits the overall rate of binding.

Cancer is a serious illness and the second leading cause of death in the US. It is caused by misregulation of the signaling cascades and uncontrolled cell proliferation. Protein kinases are the key players in the regulatory machinery, and in the last 20 years they have become an increasingly important drug targets. The main requirement for such drugs is selectivity: ideally, only one kinase should be affected by an individual drug. The prominent success story is an FDA-approved chronic myeloid leukemia therapeutic Gleevec. Gleevec is a highly specific inhibitor of the Abl kinase and has very little potency towards a closely related Src subfamily of kinases. This specificity is particularly intriguing considering that the drug-binding pockets look nearly identical in these kinases. We combined NMR and pre-steady-state kinetic experiments to monitor Gleevec binding in real time and with residue-specific precision. Our data establish a novel model that fully and quantitatively explains the observed 3000-fold difference in Gleevec's affinities to Abl and Src, solving a long-standing paradox in the field. These results reveal the general principles of kinase-drug interactions and highlight the molecular mechanisms behind nanomolar affinities found in clinically relevant therapeutics.