Considerations Around Structure-Based Drug Discovery for KRAS Using DOCK.

Molecular docking is a popular computational tool in drug discovery. Leveraging structural information, docking software predicts binding poses of small molecules to cavities on the surfaces of proteins. Virtual screening for ligand discovery is a useful application of docking software. In this chapter, using the enigmatic KRAS protein as an example system, we endeavor to teach the reader about best practices for performing molecular docking with UCSF DOCK. We discuss methods for virtual screening and docking molecules on KRAS. We present the following six points to optimize our docking setup for prosecuting a virtual screen: protein structure choice, pocket selection, optimization of the scoring function, modification of sampling spheres and sampling procedures, choosing an appropriate portion of chemical space to dock, and the choice of which top scoring molecules to pick for purchase.

Structural dynamics of RAF1-HSP90-CDC37 and HSP90 complexes reveal asymmetric client interactions and key structural elements.

RAF kinases are integral to the RAS-MAPK signaling pathway, and proper RAF1 folding relies on its interaction with the chaperone HSP90 and the cochaperone CDC37. Understanding the intricate molecular interactions governing RAF1 folding is crucial for comprehending this process. Here, we present a cryo-EM structure of the closed-state RAF1-HSP90-CDC37 complex, where the C-lobe of the RAF1 kinase domain binds to one side of the HSP90 dimer, and an unfolded N-lobe segment of the RAF1 kinase domain threads through the center of the HSP90 dimer. CDC37 binds to the kinase C-lobe, mimicking the N-lobe with its HxNI motif. We also describe structures of HSP90 dimers without RAF1 and CDC37, displaying only N-terminal and middle domains, which we term the semi-open state. Employing 1 µs atomistic simulations, energetic decomposition, and comparative structural analysis, we elucidate the dynamics and interactions within these complexes. Our quantitative analysis reveals that CDC37 bridges the HSP90-RAF1 interaction, RAF1 binds HSP90 asymmetrically, and that HSP90 structural elements engage RAF1's unfolded region. Additionally, N- and C-terminal interactions stabilize HSP90 dimers, and molecular interactions in HSP90 dimers rearrange between the closed and semi-open states. Our findings provide valuable insight into the contributions of HSP90 and CDC37 in mediating client folding.

Hydrophobic interactions dominate the recognition of a KRAS G12V neoantigen.

Specificity remains a major challenge to current therapeutic strategies for cancer. Mutation associated neoantigens (MANAs) are products of genetic alterations, making them highly specific therapeutic targets. MANAs are HLA-presented (pHLA) peptides derived from intracellular mutant proteins that are otherwise inaccessible to antibody-based therapeutics. Here, we describe the cryo-EM structure of an antibody-MANA pHLA complex. Specifically, we determine a TCR mimic (TCRm) antibody bound to its MANA target, the KRASG12V peptide presented by HLA-A*03:01. Hydrophobic residues appear to account for the specificity of the mutant G12V residue. We also determine the structure of the wild-type G12 peptide bound to HLA-A*03:01, using X-ray crystallography. Based on these structures, we perform screens to validate the key residues required for peptide specificity. These experiments led us to a model for discrimination between the mutant and the wild-type peptides presented on HLA-A*03:01 based exclusively on hydrophobic interactions.

Sodium/Iodide Symporter Metastable Intermediates Provide Insights into Conformational Transition between Principal Thermodynamic States.

The Sodium/Iodide Symporter (NIS), a 13-helix transmembrane protein found in the thyroid and other tissues, transports iodide, a required constituent of thyroid hormones T3 and T4. Despite extensive experimental information and clinical data, structural details of the intermediate microstates comprising the conformational transition of NIS between its inwardly and outwardly open states remain unresolved. We present data from a combination of enhanced sampling and transition path molecular dynamics (MD) simulations that elucidate the principal intermediate states comprising the inwardly to outwardly open transition of fully bound and apo NIS under an enforced ionic gradient. Our findings suggest that in both the absence and presence of bound physiological ions, NIS principally occupies a proximally inward to inwardly open state. When fully bound, NIS is also found to occupy a rare "inwardly occluded" state. The results of this work provide novel insight into the populations of NIS intermediates and the free energy landscape comprising the conformational transition, adding to a mechanistic understanding of NIS ion transport. Moreover, the knowledge gained from this approach can serve as a basis for studies of NIS mutants to target therapeutic interventions.

Class I PI3K Biology.

This chapter is an introduction to phosphoinositide 3-kinases (PI3K), with class I PI3Ks as the central focus. First, the various PI3K isoforms in class I are presented with emphasis on their overall structure, subunits, subunit constitutive domains, domain-domain interactions, and functional relevance. This structural analysis is followed by a comprehensive history of seminal investigations into PI3K activity. Next, we highlight the divergent roles of the isoforms: PI3Kα, PI3Kß, PI3Kδ, and PI3Kγ. This section details signaling pathways in which these PI3K isoforms are involved, including the key upstream regulators of PI3K activity and some downstream cellular effects. Nodes of the PI3K pathway are also presented. Inhibitors of some isoforms are discussed to give an overview of the basis of some immunotherapies that are being used to target cell signaling. Finally, the chapter ends with a discussion of the dysregulation of PI3Ks in diseases including APDS, asthma, arthritis, and oncogenic mutations.

The Iodide Transport Defect-Causing Y348D Mutation in the Na+/I- Symporter Renders the Protein Intrinsically Inactive and Impairs Its Targeting to the Plasma Membrane.

Background: The sodium/iodide (Na+/I-) symporter (NIS) mediates active transport of I- into the thyroid gland. Mutations in the SLC5A5 gene, which encodes NIS, cause I- transport defects (ITDs)-which, if left untreated, lead to congenital hypothyroidism and consequent cognitive and developmental deficiencies. The ITD-causing NIS mutation Y348D, located in transmembrane segment (TMS) 9, was reported in three Sudanese patients. Methods: We generated cDNAs coding for Y348D NIS and mutants with other hydrophilic and hydrophobic amino acid substitutions at position 348 and transfected them into cells. The activity of the resulting mutants was quantitated by radioiodide transport assays. NIS glycosylation was investigated by Western blotting after endoglycosidase H (Endo H) and PNGase-F glycosidase treatment. Subcellular localization of the mutant proteins was ascertained by flow cytometry analysis, cell surface biotinylation, and immunofluorescence. The intrinsic activity of Y348D was studied by measuring radioiodide transport in membrane vesicles prepared from Y348D-NIS-expressing cells. Our NIS homology models and molecular dynamics simulations were used to identify residues that interact with Y348 and investigate possible interactions between Y348 and the membrane. The sequences of several Slc5 family transporters were aligned, and a phylogenetic tree was generated in ClustalX. Results: Cells expressing Y348D NIS transport no I-. Furthermore, Y348D NIS is only partially glycosylated, is retained intracellularly, and is intrinsically inactive. Hydrophilic residues other than Asp at position 348 also yield NIS proteins that fail to be targeted to the plasma membrane (PM), whereas hydrophobic residues at this position, which we show do not interact with the membrane, rescue PM targeting and function. Conclusions: Y348D NIS does not reach the PM and is intrinsically inactive. Hydrophobic amino acid substitutions at position 348, however, preserve NIS activity. Our findings are consistent with our homology model's prediction that Y348 should face the side opposite the TMS9 residues that coordinate Na+ and participate in Na+ transport, and with the notion that Y348 interacts only with hydrophobic residues. Hydrophilic or charged residues at position 348 have deleterious effects on NIS PM targeting and activity, whereas a hydrophobic residue at this position rescues NIS activity.

Allosteric Activation of PI3Kα Results in Dynamic Access to Catalytically Competent Conformations.

Class I phosphoinositide-3-kinases (PI3Ks) phosphorylate PIP2 at its 3' inositol position to generate PIP3, a second messenger that influences signaling cascades regulating cellular growth, survival, and proliferation. Previous studies have suggested that PI3Kα activation involves dislodging the p85α nSH2 domain from the p110α catalytic subunit by binding activated receptor tyrosine kinases. We carried out molecular dynamics simulations to determine, mechanistically and structurally, how PI3Kα conformations are influenced by physiological effectors and the nSH2 domain. We demonstrate that changes in protein dynamics mediated by allosteric regulation significantly increase the population of catalytically competent states without changing the enzyme ground-state structure. Furthermore, we demonstrate that modulation of active-site residue interactions with enzyme substrates can reciprocally influence nSH2 domain dynamics. Together, these results suggest that dynamic allostery plays a role in populating the catalytically competent conformation of PI3Kα, and provide a key platform for the design of novel chemotherapeutic PI3Kα inhibitors.

Comparison of the Conformations of KRAS Isoforms, K-Ras4A and K-Ras4B, Points to Similarities and Significant Differences.

Human HRAS, KRAS, and NRAS genes encode four isoforms of Ras, a p21 GTPase. Mutations in KRAS account for the majority of RAS-driven cancers. The KRAS has two splice variants, K-Ras4A and K-Ras4B. Due to their reversible palmitoylation, K-Ras4A and N-Ras have bimodal signaling states. K-Ras4A and K-Ras4B differ in four catalytic domain residues (G151R/D153E/K165Q/H166Y) and in their disordered C-terminal hypervariable region (HVR). In K-Ras4A, the HVR is not as strongly positively charged as in K-Ras4B (+6e vs +9e). Here, we performed all-atom molecular dynamics simulations to elucidate isoform-specific differences between the two splice variants. We observe that the catalytic domain of GDP-bound K-Ras4A has a more exposed nucleotide binding pocket than K-Ras4B, and the dynamic fluctuations in switch I and II regions also differ; both factors may influence guanine-nucleotide exchange. We further observe that like K-Kas4B, full-length K-Ras4A exhibits nucleotide-dependent HVR fluctuations; however, these fluctuations differ between the GDP-bound forms of K-Ras4A and K-Ras4B. Unlike K-Ras4B where the HVR tends to cover the effector binding region, in K-Ras4A, autoinhibited states are unstable. With lesser charge, the K-Ras4A HVR collapses on itself, making it less available for binding the catalytic domain. Since the HVRs of N- and H-Ras are weakly charged (+1e and +2e, respectively), autoinhibition may be a unique feature of K-Ras4B.

A New View of Ras Isoforms in Cancers.

Does small GTPase K-Ras4A have a single state or two states, one resembling K-Ras4B and the other N-Ras? A recent study of K-Ras4A made the remarkable observation that even in the absence of the palmitoyl, K-Ras4A can be active at the plasma membrane. Importantly, this suggests that K-Ras4A may exist in two distinct signaling states. In state 1, K-Ras4A is only farnesylated, like K-Ras4B; in state 2, farnesylated and palmitoylated, like N-Ras. The K-Ras4A hypervariable region sequence is positively charged, in between K-Ras4B and N-Ras. Taken together, this raises the possibility that the farnesylated but nonpalmitoylated state 1, like K-Ras4B, binds calmodulin and is associated with colorectal and other adenocarcinomas like lung cancer and pancreatic ductal adenocarcinoma. On the other hand, state 2 may be associated with melanoma and other cancers where N-Ras is a major contributor, such as acute myeloid leukemia. Importantly, H-Ras has two, singly and doubly, palmitoylated states that may also serve distinct functional roles. The multiple signaling states of palmitoylated Ras isoforms question the completeness of small GTPase Ras isoform statistics in different cancer types and call for reevaluation of concepts and protocols. They may also call for reconsideration of oncogenic Ras therapeutics.

Prediction of passive blood-brain partitioning: straightforward and effective classification models based on in silico derived physicochemical descriptors.

The distribution of compounds between blood and brain is a very important consideration for new candidate drug molecules. In this paper, we describe the derivation of two linear discriminant analysis (LDA) models for the prediction of passive blood-brain partitioning, expressed in terms of logBB values. The models are based on computationally derived physicochemical descriptors, namely the octanol/water partition coefficient (logP), the topological polar surface area (TPSA) and the total number of acidic and basic atoms, and were obtained using a homogeneous training set of 307 compounds, for all of which the published experimental logBB data had been determined in vivo. In particular, since molecules with logBB>0.3 cross the blood-brain barrier (BBB) readily while molecules with logBB<-1 are poorly distributed to the brain, on the basis of these thresholds we derived two distinct models, both of which show a percentage of good classification of about 80%. Notably, the predictive power of our models was confirmed by the analysis of a large external dataset of compounds with reported activity on the central nervous system (CNS) or lack thereof. The calculation of straightforward physicochemical descriptors is the only requirement for the prediction of the logBB of novel compounds through our models, which can be conveniently applied in conjunction with drug design and virtual screenings.

Cyclin-dependent kinase 5 modulates the transcriptional activity of the mineralocorticoid receptor and regulates expression of brain-derived neurotrophic factor.

Glucocorticoids, major end effectors of the stress response, play an essential role in the homeostasis of the central nervous system (CNS) and contribute to memory consolidation and emotional control through their intracellular receptors, the glucocorticoid and mineralocorticoid receptors. Cyclin-dependent kinase 5 (CDK5), on the other hand, plays important roles in the morphogenesis and functions of the central nervous system, and its aberrant activation has been associated with development of neurodegenerative disorders. We previously reported that CDK5 phosphorylated the glucocorticoid receptor and modulated its transcriptional activity. Here we found that CDK5 also regulated mineralocorticoid receptor-induced transcriptional activity by phosphorylating multiple serine and threonine residues located in its N-terminal domain through physical interaction. Aldosterone and dexamethasone, respectively, increased and suppressed mRNA/protein expression of brain-derived neurotrophic factor (BDNF) in rat cortical neuronal cells, whereas the endogenous glucocorticoid corticosterone showed a biphasic effect. CDK5 enhanced the effect of aldosterone and dexamethasone on BDNF expression. Because this neurotrophic factor plays critical roles in neuronal viability, synaptic plasticity, consolidation of memory, and emotional changes, we suggest that aberrant activation of CDK5 might influence these functions through corticosteroid receptors/BDNF.