Agnieszka Bronowska

A reliable docking/scoring scheme based on the semiempirical quantum mechanical PM6-DH2 method accurately covering dispersion and H-bonding: HIV-1 protease with 22 ligands.

In this study, we introduce a fast and reliable rescoring scheme for docked complexes based on a semiempirical quantum mechanical PM6-DH2 method. The method utilizes a PM6-based Hamiltonian with corrections for dispersion energy and hydrogen bonds. The total score is constructed as the sum of the PM6-DH2 interaction enthalpy, the empirical force field (AMBER) interaction entropy, and the sum of the deformation (PM6-DH2, SMD) and the desolvation (SMD) energies of the ligand. The main advantage of the procedure is the fact that we do not add any empirical parameter for either an individual component of the total score or an individual protein-ligand complex. This rescoring method is applied to a very challenging system, namely, the HIV-1 protease with a set of ligands. As opposed to the conventional DOCK procedure, the PM6-DH2 rescoring based on all of the terms distinguishes between binders and nonbinders and provides a reliable correlation of the theoretical and experimental binding free energies. Such a dramatic improvement, resulting from the PM6-DH2 rescoring of all the complexes, provides a valuable yet inexpensive tool for rational drug discovery and de novo ligand design.

A proton relay enhances H2O2 sensitivity of GAPDH to facilitate metabolic adaptation.

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is sensitive to reversible oxidative inactivation by hydrogen peroxide (H2O2). Here we show that H2O2 reactivity of the active site thiolate (C152) is catalyzed by a previously unrecognized mechanism based on a dedicated proton relay promoting leaving group departure. Disruption of the peroxidatic reaction mechanism does not affect the glycolytic activity of GAPDH. Therefore, specific and separate mechanisms mediate the reactivity of the same thiolate nucleophile toward H2O2 and glyceraldehyde 3-phosphate, respectively. The generation of mutants in which the glycolytic and peroxidatic activities of GAPDH are comprehensively uncoupled allowed for a direct assessment of the physiological relevance of GAPDH H2O2 sensitivity. Using yeast strains in which wild-type GAPDH was replaced with H2O2-insensitive mutants retaining full glycolytic activity, we demonstrate that H2O2 sensitivity of GAPDH is a key component of the cellular adaptive response to increased H2O2 levels.

Drug profiling: knowing where it hits.

Off-target hits of drugs can lead to serious adverse effects or, conversely, to unforeseen alternative medical utility. Selectivity profiling against large panels of potential targets is essential for the drug discovery process to minimize attrition and maximize therapeutic utility. Lately, it has become apparent that drug promiscuity (polypharmacology) goes well beyond target families; therefore, lowering the profiling costs and expanding the coverage of targets is an industry-wide challenge to improve predictions. Here, we review current and promising drug profiling alternatives and commercial solutions in these exciting emerging fields.

Phosphatidylinositol 4,5-bisphosphate triggers activation of focal adhesion kinase by inducing clustering and conformational changes

Significance Nonreceptor tyrosine kinases are major players in cell signaling. Among them, focal adhesion kinase (FAK) is the key integrator of signals from growth factors and cell adhesion. In cancer, FAK is frequently overexpressed, and by promoting adhesion to the tumor stroma and ECM, FAK provides important signals for tumor invasion and metastasis. Although autoinhibitory mechanisms have previously been described and the players involved in FAK regulation are largely known, on a mechanistic level, FAK activation is currently not understood. Here, we present a multidisciplinary approach demonstrating a multistep mechanism resulting in FAK activation. This mechanistic insight enables the design of alternative strategies for the discovery of potential anticancer drugs that inhibit both catalytic and scaffolding functions of FAK with high specificity. Focal adhesion kinase (FAK) is a nonreceptor tyrosine kinase (NRTK) with key roles in integrating growth and cell matrix adhesion signals, and FAK is a major driver of invasion and metastasis in cancer. Cell adhesion via integrin receptors is well known to trigger FAK signaling, and many of the players involved are known; however, mechanistically, FAK activation is not understood. Here, using a multidisciplinary approach, including biochemical, biophysical, structural, computational, and cell biology approaches, we provide a detailed view of a multistep activation mechanism of FAK initiated by phosphatidylinositol-4,5-bisphosphate [PI(4,5)P2]. Interestingly, the mechanism differs from canonical NRTK activation and is tailored to the dual catalytic and scaffolding function of FAK. We find PI(4,5)P2 induces clustering of FAK on the lipid bilayer by binding a basic region in the regulatory 4.1, ezrin, radixin, moesin homology (FERM) domain. In these clusters, PI(4,5)P2 induces a partially open FAK conformation where the autophosphorylation site is exposed, facilitating efficient autophosphorylation and subsequent Src recruitment. However, PI(4,5)P2 does not release autoinhibitory interactions; rather, Src phosphorylation of the activation loop in FAK results in release of the FERM/kinase tether and full catalytic activation. We propose that PI(4,5)P2 and its generation in focal adhesions by the enzyme phosphatidylinositol 4-phosphate 5-kinase type Iγ are important in linking integrin signaling to FAK activation.

Comparison of Entropic Contributions to Binding in a “Hydrophilic” versus “Hydrophobic” Ligand−Protein Interaction

In the present study we characterize the thermodynamics of binding of histamine to recombinant histamine-binding protein (rRaHBP2), a member of the lipocalin family isolated from the brown-ear tick Rhipicephalus appendiculatus. The binding pocket of this protein contains a number of charged residues, consistent with histamine binding, and is thus a typical example of a “hydrophilic” binder. In contrast, a second member of the lipocalin family, the recombinant major urinary protein (rMUP), binds small hydrophobic ligands, with a similar overall entropy of binding in comparison with rRaHBP2. Having extensively studied ligand binding thermodynamics for rMUP previously, the data we obtained in the present study for HBP enables a comparison of the driving forces for binding between these classically distinct binding processes in terms of entropic contributions from ligand, protein, and solvent. In the case of rRaHBP2, we find favorable entropic contributions to binding from desolvation of the ligand; however, the overall entropy of binding is unfavorable due to a dominant unfavorable contribution arising from the loss of ligand degrees of freedom, together with the sequestration of solvent water molecules into the binding pocket in the complex. This contrasts with binding in rMUP where desolvation of the protein binding pocket makes a minor contribution to the overall entropy of binding given that the pocket is substantially desolvated prior to binding.

Mechanism of Focal Adhesion Kinase Mechanosensing

Mechanosensing at focal adhesions regulates vital cellular processes. Here, we present results from molecular dynamics (MD) and mechano-biochemical network simulations that suggest a direct role of Focal Adhesion Kinase (FAK) as a mechano-sensor. Tensile forces, propagating from the membrane through the PIP2 binding site of the FERM domain and from the cytoskeleton-anchored FAT domain, activate FAK by unlocking its central phosphorylation site (Tyr576/577) from the autoinhibitory FERM domain. Varying loading rates, pulling directions, and membrane PIP2 concentrations corroborate the specific opening of the FERM-kinase domain interface, due to its remarkably lower mechanical stability compared to the individual alpha-helical domains and the PIP2-FERM link. Analyzing downstream signaling networks provides further evidence for an intrinsic mechano-signaling role of FAK in broadcasting force signals through Ras to the nucleus. This distinguishes FAK from hitherto identified focal adhesion mechano-responsive molecules, allowing a new interpretation of cell stretching experiments.

Quantum Mechanical Scoring: Structural and Energetic Insights into Cyclin-Dependent Kinase 2 Inhibition by Pyrazolo[1,5-a]pyrimidines

A quantum mechanics (QM)-based scoring function has been applied to complexes of cyclin-dependent kinase 2 (CDK2) and thirty-one pyrazolo[1,5-a]pyrimidine-based inhibitors and their bioisosteres. A hybrid three-layer QM/MM setup (DFT-D/PM6-D3H4X/AMBER in generalized Born solvent) was used here for the first time as an extension of our previous full QM and SQM/MM (SQM means semiempirical QM) approaches. Two approaches to obtain the structures of the CDK2/inhibitor complexes were examined: i) building the modifications from one X-ray structure available coupled with a conformational search and ii) docking the compounds into CDK2. The QM-based scoring entailed a QM/SQM/MM optimization followed by calculations of the binding scores which were subsequently correlated with the experimental binding free energies. The correlation for the building protocol was good (r(2) = 0.64, predictive index = 0.81), whereas the docking approach failed. A decomposition of the interaction energies to ligand fragments enabled us to rationalize the differences in the binding affinities. In conclusion, we have developed and refined a QM-based scoring protocol and successfully applied it to reproduce the binding affinities in congeneric series of CDK2 inhibitors and to rationalize their potency. We thus propose that such a tool can be used in computer-aided rational drug design.

Residual ligand entropy in the binding of p-substituted benzenesulfonamide ligands to bovine carbonic anhydrase II.

In studies on the thermodynamics of ligand-protein interactions, it is often assumed that the configurational and conformational entropy of the ligand is zero in the bound state (i.e., the ligand is rigidly fixed in the binding pocket). However, there is little direct experimental evidence for this assumption, and in the case of binding of p-substituted benzenesulfonamide inhibitors to bovine carbonic anhydrase II (BCA II), the observed thermodynamic binding signature derived from isothermal titration calorimetry experiments leads indirectly to the conclusion that a considerable degree of residual entropy remains in the bound ligand. Specifically, the entropy of binding increases with glycine chain length n, and strong evidence exists that this thermodynamic signature is not driven by solvent reorganization. By use of heteronuclear (15)N NMR relaxation measurements in a series (n = 1-6) of (15)N-glycine-enriched ligands, we find that the observed thermodynamic binding signature cannot be explained by residual ligand dynamics in the bound state, but rather results from the indirect influence of ligand chain length on protein dynamics.

Force-Sensitive Autoinhibition of the von Willebrand Factor Is Mediated by Interdomain Interactions

Von Willebrand factor (VWF) plays a central role in hemostasis. Triggered by shear-stress, it adheres to platelets at sites of vascular injury. Inactivation of VWF has been associated to the shielding of its adhesion sites and proteolytic cleavage. However, the molecular nature of this shielding and its coupling to cleavage under shear-forces in flowing blood remain unknown. In this study, we describe, to our knowledge, a new force-sensory mechanism for VWF-platelet binding, which addresses these questions, based on a combination of molecular dynamics (MD) simulations, atomic force microscopy (AFM), and microfluidic experiments. Our MD simulations demonstrate that the VWF A2 domain targets a specific region at the VWF A1 domain, corresponding to the binding site of the platelet glycoprotein Ibα (GPIbα) receptor, thereby causing its blockage. This implies autoinhibition of the VWF for the binding of platelets mediated by the A1-A2 protein-protein interaction. During force-probe MD simulations, a stretching force dissociated the A1A2 complex, thereby unblocking the GPIbα binding site. Dissociation was found to be coupled to the unfolding of the A2 domain, with dissociation predominantly occurring before exposure of the cleavage site in A2, an observation that is supported by our AFM experiments. This suggests that the A2 domain prevents platelet binding in a force-dependent manner, ensuring that VWF initiates hemostasis before inactivation by proteolytic cleavage. Microfluidic experiments with an A2-deletion VWF mutant resulted in increased platelet binding, corroborating the key autoinhibitory role of the A2 domain within VWF multimers. Overall, autoinhibition of VWF mediated by force-dependent interdomain interactions offers the molecular basis for the shear-sensitive growth of VWF-platelet aggregates, and might be similarly involved in shear-induced VWF self-aggregation and other force-sensing functions in hemostasis.

Allosteric Regulation of Focal Adhesion Kinase by PIP2 and ATP

Focal adhesion kinase (FAK) is a nonreceptor tyrosine kinase that regulates cell signaling, proliferation, migration, and development. A major mechanism of regulation of FAK activity is an intramolecular autoinhibitory interaction between two of its domains--the catalytic and FERM domains. Upon cell adhesion to the extracellular matrix, FAK is being translocated toward focal adhesion sites and activated. Interactions of FAK with phosphoinositide phosphatidylinsositol-4,5-bis-phosphate (PIP₂) are required to activate FAK. However, the molecular mechanism of the activation remains poorly understood. Recent fluorescence resonance energy transfer experiments revealed a closure of the FERM-kinase interface upon ATP binding, which is reversed upon additional binding of PIP₂. Here, we addressed the allosteric regulation of FAK by performing all-atom molecular-dynamics simulations of a FAK fragment containing the catalytic and FERM domains, and comparing the dynamics in the absence or presence of ATP and PIP₂. As a major conformational change, we observe a closing and opening motion upon ATP and additional PIP₂ binding, respectively, in good agreement with the fluorescence resonance energy transfer experiments. To reveal how the binding of the regulatory PIP₂ to the FERM F2 lobe is transduced to the very distant F1/N-lobe interface, we employed force distribution analysis. We identified a network of mainly charged residue-residue interactions spanning from the PIP₂ binding site to the distant interface between the kinase and FERM domains, comprising candidate residues for mutagenesis to validate the predicted mechanism of FAK activation.