NMR-Chemical-Shift-Driven Protocol Reveals the Cofactor-Bound, Complete Structure of Dynamic Intermediates of the Catalytic Cycle of Oncogenic KRAS G12C Protein and the Significance of the Mg2+ Ion.

In this work, catalytically significant states of the oncogenic G12C variant of KRAS, those of Mg2+-free and Mg2+-bound GDP-loaded forms, have been determined using CS-Rosetta software and NMR-data-driven molecular dynamics simulations. There are several Mg2+-bound G12C KRAS/GDP structures deposited in the Protein Data Bank (PDB), so this system was used as a reference, while the structure of the Mg2+-free but GDP-bound state of the RAS cycle has not been determined previously. Due to the high flexibility of the Switch-I and Switch-II regions, which also happen to be the catalytically most significant segments, only chemical shift information could be collected for the most important regions of both systems. CS-Rosetta was used to derive an "NMR ensemble" based on the measured chemical shifts, which, however, did not contain the nonprotein components of the complex. We developed a torsional restraint set for backbone torsions based on the CS-Rosetta ensembles for MD simulations, overriding the force-field-based parametrization in the presence of the reinserted cofactors. This protocol (csdMD) resulted in complete models for both systems that also retained the structural features and heterogeneity defined by the measured chemical shifts and allowed a detailed comparison of the Mg2+-bound and Mg2+-free states of G12C KRAS/GDP.

Covalent fragment mapping of KRasG12C revealed novel chemotypes with in vivo potency.

G12C mutant KRas is considered druggable by allele-specific covalent inhibitors due to the nucleophilic character of the oncogenic mutant cysteine at position 12. Discovery of these inhibitors requires the optimization of both covalent and noncovalent interactions. Here, we report covalent fragment screening of our electrophilic fragment library of diverse non-covalent scaffolds equipped with 40 different electrophilic functionalities to identify fragments as suitable starting points targeting Cys12. Screening the library against KRasG12C using Ellman's free thiol assay, followed by protein NMR and cell viability assays, resulted in two potential inhibitor chemotypes. Characterization of these scaffolds in in vitro cellular- and in vivo xenograft models revealed them as promising starting points for covalent drug discovery programs.

The Importance of Mg2+ -Free State in Nucleotide Exchange of Oncogenic K-Ras Mutants.

For efficient targeting of oncogenic K-Ras interaction sites, a mechanistic picture of the Ras-cycle is necessary. Herein, we used NMR relaxation techniques and molecular dynamics simulations to decipher the role of slow dynamics in wild-type and three oncogenic P-loop mutants of K-Ras. Our measurements reveal a dominant two-state conformational exchange on the ms timescale in both GDP- and GTP-bound K-Ras. The identified low-populated higher energy state in GDP-loaded K-Ras has a conformation reminiscent of a nucleotide-bound/Mg2+ -free state characterized by shortened ß2/ß3-strands and a partially released switch-I region preparing K-Ras for the interaction with the incoming nucleotide exchange factor and subsequent reactivation. By providing insight into mutation-specific differences in K-Ras structural dynamics, our systematic analysis improves our understanding of prolonged K-Ras signaling and may aid the development of allosteric inhibitors targeting nucleotide exchange in K-Ras.

Investigating thermal stability based on the structural changes of lactase enzyme by several orthogonal methods.

Thermal stability of lactase (ß-galactosidase) enzyme has been studied by a variety of physico-chemical methods. ß-galactosidase is the main active ingredient of medications for lactose intolerance. It is typically produced industrially by the Aspergillus oryzae filamentous fungus. Lactase was used as a model to help understand thermal stability of enzyme-type biopharmaceuticals. Enzyme activity (hydrolyzation of lactose) of ß-galactosidase was determined after storing the solid enzyme substance at various temperatures. For a better understanding of the relationship between structure and activity changes we determined the mass and size of the molecules with gel electrophoresis and dynamic light scattering and detected aggregation processes. A bottom-up proteomic procedure was used to determine the primary amino acid sequence and to investigate changes in the N-glycosylation pattern of the protein. NMR and CD spectroscopic methods were used to observe changes in higher order structures and to reveal relationships between structural and functional changes.

Bacterial fermentation and isotope labelling optimized for amyloidogenic proteins.

We developed a cost sensitive isotope labelling procedure using a fed-batch fermentation method and tested its efficiency producing the 15 N-, 13 C- and 15 N/13 C-labelled variants of an amyloidogenic miniprotein (E5: EEEAVRLYIQWLKEGGPSSGRPPPS). E5 is a surface active protein, which forms amyloids in solution. Here, we confirm, using both PM-IRRAS and AFM measurements, that the air-water interface triggers structural rearrangement and promotes the amyloid formation of E5, and thus it is a suitable test protein to work out efficient isotope labelling schemes even for such difficult sequences. E. coli cells expressing the recombinant, ubiquitin-fused miniprotein were grown in minimal media containing either unlabelled nutrients, or 15 N-NH4 Cl and/or 13 C-D-Glc. The consumption rates of NH4 Cl and D-Glc were quantitatively monitored during fermentation and their ratio was established to be 1:5 (for NH4 Cl: D-Glc). One- and two-step feeding schemes were custom-optimized to enhance isotope incorporation expressing five different E5 miniprotein variants. With the currently optimized protocols we could achieve a 1.5- to 5-fold increase of yields of several miniproteins coupled to a similar magnitude of cost reduction as compared to flask labelling protocols.

Assessment of Tractable Cysteines for Covalent Targeting by Screening Covalent Fragments.

Targeted covalent inhibition and the use of irreversible chemical probes are important strategies in chemical biology and drug discovery. To date, the availability and reactivity of cysteine residues amenable for covalent targeting have been evaluated by proteomic and computational tools. Herein, we present a toolbox of fragments containing a 3,5-bis(trifluoromethyl)phenyl core that was equipped with chemically diverse electrophilic warheads showing a range of reactivities. We characterized the library members for their reactivity, aqueous stability and specificity for nucleophilic amino acids. By screening this library against a set of enzymes amenable for covalent inhibition, we showed that this approach experimentally characterized the accessibility and reactivity of targeted cysteines. Interesting covalent fragment hits were obtained for all investigated cysteine-containing enzymes.

Dynamically encoded reactivity of Ras enzymes: opening new frontiers for drug discovery.

Decoding molecular flexibility in order to understand and predict biological processes-applying the principles of dynamic-structure-activity relationships (DSAR)-becomes a necessity when attempting to design selective and specific inhibitors of a protein that has overlapping interaction surfaces with its upstream and downstream partners along its signaling cascade. Ras proteins are molecular switches that meet this definition perfectly. The close-lying P-loop and the highly flexible switch I and switch II regions are the site of nucleotide-, assisting-, and effector-protein binding. Oncogenic mutations that also appear in this region do not cause easily characterized overall structural changes, due partly to the inherent conformational heterogeneity and pliability of these segments. In this review, we present an overview of the results obtained using approaches targeting Ras dynamics, such as nuclear magnetic resonance (NMR) measurements and experiment-based modeling calculations (mostly molecular dynamics (MD) simulations). These methodologies were successfully used to decipher the mutant- and isoform-specific nature of certain transient states, far-lying allosteric sites, and the internal interaction networks, as well as the interconnectivity of the catalytic and membrane-binding regions. This opens new therapeutic potential: the discovered interaction hotspots present hitherto not targeted, selective sites for drug design efforts in diverse locations of the protein matrix.

Tumor-Suppressor p53TAD1-60 Forms a Fuzzy Complex with Metastasis-Associated S100A4: Structural Insights and Dynamics by an NMR/MD Approach.

Conformationally flexible protein complexes represent a major challenge for structural and dynamical studies. We present herein a method based on a hybrid NMR/MD approach to characterize the complex formed between the disordered p53TAD1-60 and the metastasis-associated S100A4. Disorder-to-order transitions of both TAD1 and TAD2 subdomains upon interaction is detected. Still, p53TAD1-60 remains highly flexible in the bound form, with residues L26, M40, and W53 being anchored to identical hydrophobic pockets of the S100A4 monomer chains. In the resulting "fuzzy" complex, the clamp-like binding of p53TAD1-60 relies on specific hydrophobic anchors and on the existence of extended flexible segments. Our results demonstrate that structural and dynamical NMR parameters (cumulative Δδ, SSP, temperature coefficients, relaxation time, hetNOE) combined with MD simulations can be used to build a structural model even if, due to high flexibility, the classical solution structure calculation is not possible.

Structural impact of GTP binding on downstream KRAS signaling.

Oncogenic RAS proteins, involved in ∼30% of human tumors, are molecular switches of various signal transduction pathways. Here we apply a new protocol for the NMR study of KRAS in its (inactive) GDP- and (activated) GTP-bound form, allowing a comprehensive analysis of the backbone dynamics of its WT-, G12C- and G12D variants. We found that Tyr32 shows opposite mobility with respect to the backbone of its surroundings: it is more flexible in the GDP-bound form while more rigid in GTP-complexes (especially in WT- and G12D-GTP). Using the G12C/Y32F double mutant, we showed that the presence of the hydroxyl group of Tyr32 has a marked effect on the G12C-KRAS-GTP system as well. Molecular dynamics simulations indicate that Tyr32 is linked to the γ-phosphate of GTP in the activated states - an arrangement shown, using QM/MM calculations, to support catalysis. Anchoring Tyr32 to the γ-phosphate contributes to the capture of the catalytic waters participating in the intrinsic hydrolysis of GTP and supports a simultaneous triple proton transfer step (catalytic water → assisting water → Tyr32 → O1G of the γ-phosphate) leading to straightforward product formation. The coupled flip of negatively charged residues of switch I toward the inside of the effector binding pocket potentiates ligand recognition, while positioning of Thr35 to enter the coordination sphere of the Mg2+ widens the pocket. Position 12 mutations do not disturb the capture of Tyr32 by the γ-phosphate, but (partially) displace Gln61, which opens up the catalytic pocket and destabilizes catalytic water molecules thus impairing intrinsic hydrolysis.

1H, 15N backbone assignment and comparative analysis of the wild type and G12C, G12D, G12V mutants of K-Ras bound to GDP at physiological pH.

K-Ras protein is a membrane-bound small GTPase acting as a molecular switch. It plays a key role in many signal transduction pathways regulating cell proliferation, differentiation, survival, etc. It alternates between its GTP-bound active and the GDP-bound inactive conformers regulated by guanine nucleotide exchange factors and GTPase activating proteins. Its most frequent oncogenic mutants are G12C, G12D, and G12V that have impaired GTPase activity, thus induce malignant tumors. Here we report the resonance assignment of the backbone 1H and 15N nuclei of K-Ras wildtype, G12C, G12D and G12V proteins' catalytic G domain (1-169 residues) in GDP-bound state, and 13C of backbone and side chains of G12C mutant at physiological pH 7.4. Triple resonance data were used to get secondary structure information and backbone dynamics of G12C, the best-known drug target among K-Ras mutants. Simultaneous investigation of G12C, G12D and G12V mutants, along with the wild type form at the very same conditions allowed us to perform a comprehensive analysis based on the combined chemical shifts to reveal the effect of mutation at G12 position on structure. Intriguingly, the G12C and G12V mutants found to be structurally very similar at the three most important regions of K-Ras (P-loop, Switch-I, Switch-II), while the G12D mutant significantly differs at P-loop and Switch-II from the wildtype as well as G12C and G12V mutants. However, in Switch-I it hardly deviates from the wildtype protein.

[Allele-specific inhibitors of mutant KRAS are in the focus of RASopathy consortium]. / Allélspecifikus inhibitorok nyomában: a RASopátia konzorcium célpontjában a KRAS fehérje onkogén mutációi.

The RASopathy consortium was built from research groups of the Budapest University of Technology and Economics, Eötvös Loránd University, Semmelweis University and two startups: KINETO Lab Ltd. and Fototronic Ltd. The goal was to design and test novel covalent and allele-specific KRAS small molecular inhibitors. KRAS is the most frequently mutated human oncogene which was unsuccessfully targeted until recently. The consortium established G12C-expressing bacterial and human cancer cell models (homo- and heterozygous variants) of lung, colorectal and pancreatic tumors. Using covalent fragment and acrylamide warhead libraries we were able to select novel candidates of small molecular G12C-specific inhibitors which were compared to published best-in-class drug candidates.

DUckCov: a Dynamic Undocking-Based Virtual Screening Protocol for Covalent Binders.

Thanks to recent guidelines, the design of safe and effective covalent drugs has gained significant interest. Other than targeting non-conserved nucleophilic residues, optimizing the noncovalent binding framework is important to improve potency and selectivity of covalent binders toward the desired target. Significant efforts have been made in extending the computational toolkits to include a covalent mechanism of protein targeting, like in the development of covalent docking methods for binding mode prediction. To highlight the value of the noncovalent complex in the covalent binding process, here we describe a new protocol using tethered and constrained docking in combination with Dynamic Undocking (DUck) as a tool to privilege strong protein binders for the identification of novel covalent inhibitors. At the end of the protocol, dedicated covalent docking methods were used to rank and select the virtual hits based on the predicted binding mode. By validating the method on JAK3 and KRas, we demonstrate how this fast iterative protocol can be applied to explore a wide chemical space and identify potent targeted covalent inhibitors.

Multiple S100 protein isoforms and C-terminal phosphorylation contribute to the paralog-selective regulation of nonmuscle myosin 2 filaments.

Nonmuscle myosin 2 (NM2) has three paralogs in mammals, NM2A, NM2B, and NM2C, which have both unique and overlapping functions in cell migration, formation of cell-cell adhesions, and cell polarity. Their assembly into homo- and heterotypic bipolar filaments in living cells is primarily regulated by phosphorylation of the N-terminally bound regulatory light chain. Here, we present evidence that the equilibrium between these filaments and single NM2A and NM2B molecules can be controlled via S100 calcium-binding protein interactions and phosphorylation at the C-terminal end of the heavy chains. Furthermore, we show that in addition to S100A4, other members of the S100 family can also mediate disassembly of homotypic NM2A filaments. Importantly, these proteins can selectively remove NM2A molecules from heterotypic filaments. We also found that tail phosphorylation (at Ser-1956 and Ser-1975) of NM2B by casein kinase 2, as well as phosphomimetic substitutions at sites targeted by protein kinase C (PKC) and transient receptor potential cation channel subfamily M member 7 (TRPM7), down-regulates filament assembly in an additive fashion. Tail phosphorylation of NM2A had a comparatively minor effect on filament stability. S100 binding and tail phosphorylation therefore preferentially disassemble NM2A and NM2B, respectively. These two distinct mechanisms are likely to contribute to the temporal and spatial sorting of the two NM2 paralogs within heterotypic filaments. The existence of multiple NM2A-depolymerizing S100 paralogs offers the potential for diverse regulatory inputs modulating NM2A filament disassembly in cells and provides functional redundancy under both physiological and pathological conditions.

Ezrin interacts with S100A4 via both its N- and C-terminal domains.

Ezrin belongs to the ERM (ezrin, radixin, moesin) protein family that has a role in cell morphology changes, adhesion and migration as an organizer of the cortical cytoskeleton by linking actin filaments to the apical membrane of epithelial cells. It is highly expressed in a variety of human cancers and promotes metastasis. Members of the Ca2+-binding EF-hand containing S100 proteins have similar pathological properties; they are overexpressed in cancer cells and involved in metastatic processes. In this study, using tryptophan fluorescence and stopped-flow kinetics, we show that S100A4 binds to the N-terminal ERM domain (N-ERMAD) of ezrin with a micromolar affinity. The binding involves the F2 lobe of the N-ERMAD and follows an induced fit kinetic mechanism. Interestingly, S100A4 binds also to the unstructured C-terminal actin binding domain (C-ERMAD) with similar affinity. Using NMR spectroscopy, we characterized the complex of S100A4 with the C-ERMAD and demonstrate that no ternary complex is simultaneously formed with the two ezrin domains. Furthermore, we show that S100A4 co-localizes with ezrin in HEK-293T cells. However, S100A4 very weakly binds to full-length ezrin in vitro indicating that the interaction of S100A4 with ezrin requires other regulatory events such as protein phosphorylation and/or membrane binding, shifting the conformational equilibrium of ezrin towards the open state. As both proteins play an important role in promoting metastasis, the characterization of their interaction could shed more light on the molecular events contributing to this pathological process.

Multilevel Changes in Protein Dynamics upon Complex Formation of the Calcium-Loaded S100A4 with a Nonmuscle Myosin IIA Tail Fragment.

Dysregulation of Ca2+ -binding S100 proteins plays important role in various diseases. The asymmetric complex of Ca2+ -bound S100A4 with nonmuscle myosin IIA has high stability and highly increased Ca2+ affinity. Here we investigated the possible causes of this allosteric effect by NMR spectroscopy. Chemical shift-based secondary-structure analysis did not show substantial changes for the complex. Backbone dynamics revealed slow-timescale local motions in the H1 helices of homodimeric S100A4; these were less pronounced in the complex form and might be accompanied by an increase in dimer stability. Different mobilities in the Ca2+ -coordinating EF-hand sites indicate that they communicate by an allosteric mechanism operating through changes in protein dynamics; this must be responsible for the elevated Ca2+ affinity. These multilevel changes in protein dynamics as conformational adaptation allow S100A4 fine-tuning of its protein-protein interactions inside the cell during Ca2+ signaling.