Use of chemotherapy and loco-regional therapy in stage IA triple-negative breast cancer and their association with oncologic outcomes: A cancer registry study.

We aimed to evaluate the role of adjuvant chemotherapy and loco-regional therapy for stage IA (pT1, pN0) triple-negative breast cancer (TNBC) in a real-world setting. We identified patients with pT1, pN0 TNBC diagnosed between 2009 and 2021 within the Baden-Württemberg cancer registry (BWCR), Germany. Overall survival (OS) was assessed using Kaplan-Meier statistics and multivariate Cox regression models (adjusted for age, use of chemotherapy, local therapy (breast conserving therapy [breast conserving surgery + radiotherapy] vs. mastectomy), and tumor histologic subtype). A total of 1231 patients with a median follow-up of 45.9 months were identified: 1.0% (12 of 1231) with pT1mi stage, 9.5% (117 of 1231) with pT1a, 23.7% (292 of 1231) with pT1b, and 65.8% (810 of 1231) with pT1c. Multivariate Cox regression analysis revealed no significant influence for the use of chemotherapy on OS in pT1b patients (HR 0.90, 95% CI 0.43-1.90). For pT1c patients with Grade 1-2 tumors, the use of chemotherapy was not significantly associated OS (HR 1.01, 95% CI 0.48-2.11) but breast conserving therapy was associated with improved OS (HR 0.41, 95% CI 0.18-0.93). For pT1c patients with Grade 3 tumors, the use of chemotherapy (HR 0.51, 95% CI 0.33-0.78) as well as breast conserving therapy (HR 0.42, 95% CI 0.23-0.76) was associated with OS. This data suggests that OS in stage IA TNBC is strongly influenced by local therapy rather than the use of chemotherapy, except for pT1c patients with Grade 3 tumors. Larger studies with longer-term follow-up are welcomed to fully inform this discussion.

Oncologic Outcomes for Different Axillary Staging Techniques in Patients with Nodal-Positive Breast Cancer Undergoing Neoadjuvant Systematic Treatment: A Cancer Registry Study.

BACKGROUND: Targeted approaches such as targeted axillary dissection (TAD) or sentinel lymph node biopsy (SLNB) showed false-negative rates of < 10% compared with axillary lymph node dissection (ALND) in patients with nodal-positive breast cancer undergoing neoadjuvant systemic treatment (NAST). We aimed to evaluate real-world oncologic outcomes for different axillary staging techniques. METHODS: We identified nodal-positive breast cancer patients undergoing NAST from 2016 to 2021 from the state cancer registry of Baden-Wuerttemberg, Germany. Invasive disease-free survival (iDFS) was assessed using Kaplan-Meier statistics and multivariate Cox regression models (adjusted for age, ypN stage, ypT stage, and tumor biologic subtype). RESULTS: A total of 2698 patients with a median follow-up of 24.7 months were identified: 2204 underwent ALND, 460 underwent SLNB (255 with ≥ 3 sentinel lymph nodes [SLNs] removed, 205 with 1-2 SLNs removed), and 34 underwent TAD. iDFS 3 years after surgery was 69.7% (ALND), 76.6% (SLNB with ≥ 3 SLNs removed), 76.7% (SLNB with < 3 SLNs removed), and 78.7% (TAD). Multivariate Cox regression analysis showed no significant influence of different axillary staging techniques on iDFS (hazard ratio [HR] for SLNB with < 3 SLNs removed 0.96, 95% confidence interval [CI] 0.62-1.50; HR for SLNB with ≥ 3 SLNs removed 0.86, 95% CI 0.56-1.3; HR for TAD 0.23, 95% CI 0.03-1.64; ALND reference), and for ypN+ (HR 1.92, 95% CI 1.49-2.49), triple-negative breast cancer (HR 2.35, 95% CI 1.80-3.06), and ypT3-4 (HR 2.93, 95% CI 2.02-4.24). CONCLUSION: These real-world data provide evidence that patient selection for de-escalated axillary surgery for patients with nodal-positive breast cancer undergoing NAST was successfully adopted and no early alarm signals of iDFS detriment were detected.

Molecular Dynamics Simulations of the Human Ecto-5'-Nucleotidase (h-ecto-5'-NT, CD73): Insights into Protein Flexibility and Binding Site Dynamics.

Human ecto-5'-nucleotidase (h-ecto-5'-NT, CD73) is a homodimeric Zn2+-binding metallophosphoesterase that hydrolyzes adenosine 5'-monophosphate (5'-AMP) to adenosine and phosphate. h-Ecto-5'-NT is a key enzyme in purinergic signaling pathways and has been recognized as a promising biological target for several diseases, including cancer and inflammatory, infectious, and autoimmune diseases. Despite its importance as a biological target, little is known about h-ecto-5'-NT dynamics, which poses a considerable challenge to the design of inhibitors of this target enzyme. Here, to explore h-ecto-5'-NT flexibility, all-atom unbiased molecular dynamics (MD) simulations were performed. Remarkable differences in the dynamics of the open (catalytically inactive) and closed (catalytically active) conformations of the apo-h-ecto-5'-NT were observed during the simulations, and the nucleotide analogue inhibitor AMPCP was shown to stabilize the protein structure in the closed conformation. Our results suggest that the large and complex domain motion that enables the h-ecto-5'-NT open/closed conformational switch is slow, and therefore, it could not be completely captured within the time scale of our simulations. Nonetheless, we were able to explore the faster dynamics of the h-ecto-5'-NT substrate binding site, which is mainly located at the C-terminal domain and well conserved among the protein's open and closed conformations. Using the TRAPP ("Transient Pockets in Proteins") approach, we identified transient subpockets close to the substrate binding site. Finally, conformational states of the substrate binding site with higher druggability scores than the crystal structure were identified. In summary, our study provides valuable insights into h-ecto-5'-NT structural flexibility, which can guide the structure-based design of novel h-ecto-5'-NT inhibitors.

Druggability Assessment in TRAPP Using Machine Learning Approaches.

Accurate protein druggability predictions are important for the selection of drug targets in the early stages of drug discovery. Because of the flexible nature of proteins, the druggability of a binding pocket may vary due to conformational changes. We have therefore developed two statistical models, a logistic regression model (TRAPP-LR) and a convolutional neural network model (TRAPP-CNN), for predicting druggability and how it varies with changes in the spatial and physicochemical properties of a binding pocket. These models are integrated into TRAnsient Pockets in Proteins (TRAPP), a tool for the analysis of binding pocket variations along a protein motion trajectory. The models, which were trained on publicly available and self-augmented datasets, show equivalent or superior performance to existing methods on test sets of protein crystal structures and have sufficient sensitivity to identify potentially druggable protein conformations in trajectories from molecular dynamics simulations. Visualization of the evidence for the decisions of the models in TRAPP facilitates identification of the factors affecting the druggability of protein binding pockets.

A workflow for exploring ligand dissociation from a macromolecule: Efficient random acceleration molecular dynamics simulation and interaction fingerprint analysis of ligand trajectories.

The dissociation of ligands from proteins and other biomacromolecules occurs over a wide range of timescales. For most pharmaceutically relevant inhibitors, these timescales are far beyond those that are accessible by conventional molecular dynamics (MD) simulation. Consequently, to explore ligand egress mechanisms and compute dissociation rates, it is necessary to enhance the sampling of ligand unbinding. Random Acceleration MD (RAMD) is a simple method to enhance ligand egress from a macromolecular binding site, which enables the exploration of ligand egress routes without prior knowledge of the reaction coordinates. Furthermore, the τRAMD procedure can be used to compute the relative residence times of ligands. When combined with a machine-learning analysis of protein-ligand interaction fingerprints (IFPs), molecular features that affect ligand unbinding kinetics can be identified. Here, we describe the implementation of RAMD in GROMACS 2020, which provides significantly improved computational performance, with scaling to large molecular systems. For the automated analysis of RAMD results, we developed MD-IFP, a set of tools for the generation of IFPs along unbinding trajectories and for their use in the exploration of ligand dynamics. We demonstrate that the analysis of ligand dissociation trajectories by mapping them onto the IFP space enables the characterization of ligand dissociation routes and metastable states. The combined implementation of RAMD and MD-IFP provides a computationally efficient and freely available workflow that can be applied to hundreds of compounds in a reasonable computational time and will facilitate the use of τRAMD in drug design.

TRAPP webserver: predicting protein binding site flexibility and detecting transient binding pockets.

The TRAnsient Pockets in Proteins (TRAPP) webserver provides an automated workflow that allows users to explore the dynamics of a protein binding site and to detect pockets or sub-pockets that may transiently open due to protein internal motion. These transient or cryptic sub-pockets may be of interest in the design and optimization of small molecular inhibitors for a protein target of interest. The TRAPP workflow consists of the following three modules: (i) TRAPP structure- generation of an ensemble of structures using one or more of four possible molecular simulation methods; (ii) TRAPP analysis-superposition and clustering of the binding site conformations either in an ensemble of structures generated in step (i) or in PDB structures or trajectories uploaded by the user; and (iii) TRAPP pocket-detection, analysis, and visualization of the binding pocket dynamics and characteristics, such as volume, solvent-exposed area or properties of surrounding residues. A standard sequence conservation score per residue or a differential score per residue, for comparing on- and off-targets, can be calculated and displayed on the binding pocket for an uploaded multiple sequence alignment file, and known protein sequence annotations can be displayed simultaneously. The TRAPP webserver is freely available at http://trapp.h-its.org.

Modeling and simulation of protein-surface interactions: achievements and challenges.

Understanding protein-inorganic surface interactions is central to the rational design of new tools in biomaterial sciences, nanobiotechnology and nanomedicine. Although a significant amount of experimental research on protein adsorption onto solid substrates has been reported, many aspects of the recognition and interaction mechanisms of biomolecules and inorganic surfaces are still unclear. Theoretical modeling and simulations provide complementary approaches for experimental studies, and they have been applied for exploring protein-surface binding mechanisms, the determinants of binding specificity towards different surfaces, as well as the thermodynamics and kinetics of adsorption. Although the general computational approaches employed to study the dynamics of proteins and materials are similar, the models and force-fields (FFs) used for describing the physical properties and interactions of material surfaces and biological molecules differ. In particular, FF and water models designed for use in biomolecular simulations are often not directly transferable to surface simulations and vice versa. The adsorption events span a wide range of time- and length-scales that vary from nanoseconds to days, and from nanometers to micrometers, respectively, rendering the use of multi-scale approaches unavoidable. Further, changes in the atomic structure of material surfaces that can lead to surface reconstruction, and in the structure of proteins that can result in complete denaturation of the adsorbed molecules, can create many intermediate structural and energetic states that complicate sampling. In this review, we address the challenges posed to theoretical and computational methods in achieving accurate descriptions of the physical, chemical and mechanical properties of protein-surface systems. In this context, we discuss the applicability of different modeling and simulation techniques ranging from quantum mechanics through all-atom molecular mechanics to coarse-grained approaches. We examine uses of different sampling methods, as well as free energy calculations. Furthermore, we review computational studies of protein-surface interactions and discuss the successes and limitations of current approaches.

Protein Binding Pocket Dynamics.

The dynamics of protein binding pockets are crucial for their interaction specificity. Structural flexibility allows proteins to adapt to their individual molecular binding partners and facilitates the binding process. This implies the necessity to consider protein internal motion in determining and predicting binding properties and in designing new binders. Although accounting for protein dynamics presents a challenge for computational approaches, it expands the structural and physicochemical space for compound design and thus offers the prospect of improved binding specificity and selectivity. A cavity on the surface or in the interior of a protein that possesses suitable properties for binding a ligand is usually referred to as a binding pocket. The set of amino acid residues around a binding pocket determines its physicochemical characteristics and, together with its shape and location in a protein, defines its functionality. Residues outside the binding site can also have a long-range effect on the properties of the binding pocket. Cavities with similar functionalities are often conserved across protein families. For example, enzyme active sites are usually concave surfaces that present amino acid residues in a suitable configuration for binding low molecular weight compounds. Macromolecular binding pockets, on the other hand, are located on the protein surface and are often shallower. The mobility of proteins allows the opening, closing, and adaptation of binding pockets to regulate binding processes and specific protein functionalities. For example, channels and tunnels can exist permanently or transiently to transport compounds to and from a binding site. The influence of protein flexibility on binding pockets can vary from small changes to an already existent pocket to the formation of a completely new pocket. Here, we review recent developments in computational methods to detect and define binding pockets and to study pocket dynamics. We introduce five different classes of protein pocket dynamics: (1) appearance/disappearance of a subpocket in an existing pocket; (2) appearance/disappearance of an adjacent pocket on the protein surface in the direct vicinity of an already existing pocket; (3) pocket breathing, which may be caused by side-chain fluctuations or backbone or interdomain vibrational motion; (4) opening/closing of a channel or tunnel, connecting a pocket inside the protein with solvent, including lid motion; and (5) the appearance/disappearance of an allosteric pocket at a site on a protein distinct from an already existing pocket with binding of a ligand to the allosteric binding site affecting the original pocket. We suggest that the class of pocket dynamics, as well as the type and extent of protein motion affecting the binding pocket, should be factors considered in choosing the most appropriate computational approach to study a given binding pocket. Furthermore, we examine the relationship between pocket dynamics classes and induced fit, conformational selection, and gating models of ligand binding on binding kinetics and thermodynamics. We discuss the implications of protein binding pocket dynamics for drug design and conclude with potential future directions for computational analysis of protein binding pocket dynamics.

When the Label Matters: Adsorption of Labeled and Unlabeled Proteins on Charged Surfaces.

Fluorescent labels are often attached to proteins to monitor binding and adsorption processes. Docking simulations for native hen egg white lysozyme (HEWL) and HEWL labeled with fluorescein isothiocyanate show that these adsorb differently on charged surfaces. Attachment of even a small label can significantly change the interaction properties of a protein. Thus, the results of experiments with fluorescently labeled proteins should be interpreted by modeling the structures and computing the interaction properties of both labeled and unlabeled species.

SDA 7: A modular and parallel implementation of the simulation of diffusional association software.

The simulation of diffusional association (SDA) Brownian dynamics software package has been widely used in the study of biomacromolecular association. Initially developed to calculate bimolecular protein-protein association rate constants, it has since been extended to study electron transfer rates, to predict the structures of biomacromolecular complexes, to investigate the adsorption of proteins to inorganic surfaces, and to simulate the dynamics of large systems containing many biomacromolecular solutes, allowing the study of concentration-dependent effects. These extensions have led to a number of divergent versions of the software. In this article, we report the development of the latest version of the software (SDA 7). This release was developed to consolidate the existing codes into a single framework, while improving the parallelization of the code to better exploit modern multicore shared memory computer architectures. It is built using a modular object-oriented programming scheme, to allow for easy maintenance and extension of the software, and includes new features, such as adding flexible solute representations. We discuss a number of application examples, which describe some of the methods available in the release, and provide benchmarking data to demonstrate the parallel performance.

Computational screening of the effects of mutations on protein-protein off-rates and dissociation mechanisms by τRAMD.

The dissociation rate, or its reciprocal, the residence time (τ), is a crucial parameter for understanding the duration and biological impact of biomolecular interactions. Accurate prediction of τ is essential for understanding protein-protein interactions (PPIs) and identifying potential drug targets or modulators for tackling diseases. Conventional molecular dynamics simulation techniques are inherently constrained by their limited timescales, making it challenging to estimate residence times, which typically range from minutes to hours. Building upon its successful application in protein-small molecule systems, τ-Random Acceleration Molecular Dynamics (τRAMD) is here investigated for estimating dissociation rates of protein-protein complexes. τRAMD enables the observation of unbinding events on the nanosecond timescale, facilitating rapid and efficient computation of relative residence times. We tested this methodology for three protein-protein complexes and their extensive mutant datasets, achieving good agreement between computed and experimental data. By combining τRAMD with MD-IFP (Interaction Fingerprint) analysis, dissociation mechanisms were characterized and their sensitivity to mutations investigated, enabling the identification of molecular hotspots for selective modulation of dissociation kinetics. In conclusion, our findings underscore the versatility of τRAMD as a simple and computationally efficient approach for computing relative protein-protein dissociation rates and investigating dissociation mechanisms, thereby aiding the design of PPI modulators.

TRAPP: a tool for analysis of transient binding pockets in proteins.

We present TRAPP (TRAnsient Pockets in Proteins), a new automated software platform for tracking, analysis, and visualization of binding pocket variations along a protein motion trajectory or within an ensemble of protein structures that may encompass conformational changes ranging from local side chain fluctuations to global backbone motions. TRAPP performs accurate grid-based calculations of the shape and physicochemical characteristics of a binding pocket for each structure and detects the conserved and transient regions of the pocket in an ensemble of protein conformations. It also provides tools for tracing the opening of a particular subpocket and residues that contribute to the binding site. TRAPP thus enables an assessment of the druggability of a disease-related target protein taking its flexibility into account.

Clinical Characterization of HCC/CCA Mixed Cancers in a Population-based Cohort.

BACKGROUND AND AIMS: Primary liver cancer (PLC) ranks among of the most common cancers worldwide. Within this group, a minority of cases displays characteristics of both hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA), known as combined hepatocellular cholangiocarcinoma (cHCC- CCA). Currently, there is no specific standardized therapy for these mixed tumors. Therefore, the aim of our study was to analyze the clinical course, treatment and outcome of cHCC-CCA patients in a European population-based registry. METHODS: We investigated 9,144 patients with PLC (6,622 HCC, 2,356 iCCA, and 166 cHCC-CCA) diagnosed between 2009 and 2020. All data were obtained from Clinical Cancer Registry of Baden-Württemberg (BW), Germany. RESULTS: In all three groups patients were predominantly male (82%, 57%, and 68% for HCC, iCCA and cHCC-CCA groups, respectively). 48% of cHCC-CCA patients were diagnosed as stage IV cancers, which was more than for HCC (31%) but less compared to CCA (64%). Overall median survival of cHCC-CCA patients was worse compared to HCC (9-13 months vs. 15.5 months, p<0.001) and rather comparable to CCA (11.8 months). CONCLUSIONS: Our data demonstrated that cHCC-CCA tumors appear to have a distinct clinical course with worse overall survival compared to HCC. Thus, identification of these cancers by histopathology is essential in order to further characterize this tumor entity and to provide accurate treatment to these patients.

A quantitative, real-time assessment of binding of peptides and proteins to gold surfaces.

Interactions of peptides and proteins with inorganic surfaces are important to both natural and artificial systems; however, a detailed understanding of such interactions is lacking. In this study, we applied new approaches to quantitatively measure the binding of amino acids and proteins to gold surfaces. Real-time surface plasmon resonance (SPR) measurements showed that TEM1-ß-lactamase inhibitor protein (BLIP) interacts only weakly with Au nanoparticles (NPs). However, fusion of three histidine residues to BLIP (3H-BLIP) resulted in a significant increase in the binding to the Au NPs, which further increased when the histidine tail was extended to six histidines (6H-BLIP). Further increasing the number of His residues had no effect on the binding. A parallel study using continuous (111)-textured Au surfaces and single-crystalline, (111)-oriented, Au islands by ellipsometry, FTIR, and localized surface plasmon resonance (LSPR) spectroscopy further confirmed the results, validating the broad applicability of Au NPs as model surfaces. Evaluating the binding of all other natural amino acid homotripeptides fused to BLIP (except Cys and Pro) showed that aromatic and positively-charged residues bind preferentially to Au with respect to small aliphatic and negatively charged residues, and that the rate of association is related to the potency of binding. The binding of all fusions was irreversible. These findings were substantiated by SPR measurements of synthesized, free, soluble tripeptides using Au-NP-modified SPR chips. Here, however, the binding was reversible allowing for determination of binding affinities that correlate with the binding potencies of the related BLIP fusions. Competition assays performed between 3H-BLIP and the histidine tripeptide (3 His) suggest that Au binding residues promote the adsorption of proteins on the surface, and by this facilitate the irreversible interaction of the polypeptide chain with Au. The binding of amino acids to Au was simulated by using a continuum solvent model, showing agreement with the experimental values. These results, together with the observed binding potencies and kinetics of the BLIP fusions and free peptides, suggest a binding mechanism that is markedly different from biological protein-protein interactions.

G Protein-Coupled Receptor-Ligand Dissociation Rates and Mechanisms from τRAMD Simulations.

There is a growing appreciation of the importance of drug-target binding kinetics for lead optimization. For G protein-coupled receptors (GPCRs), which mediate signaling over a wide range of time scales, the drug dissociation rate is often a better predictor of in vivo efficacy than binding affinity, although it is more challenging to compute. Here, we assess the ability of the τ-Random Acceleration Molecular Dynamics (τRAMD) approach to reproduce relative residence times and reveal dissociation mechanisms and the effects of allosteric modulation for two important membrane-embedded drug targets: the ß2-adrenergic receptor and the muscarinic acetylcholine receptor M2. The dissociation mechanisms observed in the relatively short RAMD simulations (in which molecular dynamics (MD) simulations are performed using an additional force with an adaptively assigned random orientation applied to the ligand) are in general agreement with much more computationally intensive conventional MD and metadynamics simulations. Remarkably, although decreasing the magnitude of the random force generally reduces the number of egress routes observed, the ranking of ligands by dissociation rate is hardly affected and agrees well with experiment. The simulations also reproduce changes in residence time due to allosteric modulation and reveal associated changes in ligand dissociation pathways.

Brownian Dynamics Simulations of Proteins in the Presence of Surfaces: Long-Range Electrostatics and Mean-Field Hydrodynamics.

Simulations of macromolecular diffusion and adsorption in confined environments can offer valuable mechanistic insights into numerous biophysical processes. In order to model solutes at atomic detail on relevant time scales, Brownian dynamics simulations can be carried out with the approximation of rigid body solutes moving through a continuum solvent. This allows the precomputation of interaction potential grids for the solutes, thereby allowing the computationally efficient calculation of forces. However, hydrodynamic and long-range electrostatic interactions cannot be fully treated with grid-based approaches alone. Here, we develop a treatment of both hydrodynamic and electrostatic interactions to include the presence of surfaces by modeling grid-based and long-range interactions. We describe its application to simulate the self-association and many-molecule adsorption of the well-characterized protein hen egg-white lysozyme to mica-like and silica-like surfaces. We find that the computational model can recover a number of experimental observables of the adsorption process and provide insights into their determinants. The computational model is implemented in the Simulation of Diffusional Association (SDA) software package.

Ligand unbinding mechanisms and kinetics for T4 lysozyme mutants from τRAMD simulations.

The protein-ligand residence time, τ, influences molecular function in biological networks and has been recognized as an important determinant of drug efficacy. To predict τ, computational methods must overcome the problem that τ often exceeds the timescales accessible to conventional molecular dynamics (MD) simulation. Here, we apply the τ-Random Acceleration Molecular Dynamics (τRAMD) method to a set of kinetically characterized complexes of T4 lysozyme mutants with small, engineered binding cavities. τRAMD yields relative ligand dissociation rates in good accordance with experiments across this diverse set of complexes that differ with regard to measurement temperature, ligand identity, protein mutation and binding cavity. τRAMD thereby allows a comprehensive characterization of the ligand egress routes and determinants of τ. Although ligand dissociation by multiple egress routes is observed, we find that egress via the predominant route determines the value of τ. We also find that the presence of a greater number of metastable states along egress pathways leads to slower protein-ligand dissociation. These physical insights could be exploited in the rational optimization of the kinetic properties of drug candidates.

Contact Map Fingerprints of Protein-Ligand Unbinding Trajectories Reveal Mechanisms Determining Residence Times Computed from Scaled Molecular Dynamics.

The binding kinetic properties of potential drugs may significantly influence their subsequent clinical efficacy. Predictions of these properties based on computer simulations provide a useful alternative to their expensive and time-consuming experimental counterparts, even at an early drug discovery stage. Herein, we perform scaled molecular dynamics (ScaledMD) simulations on a set of 27 ligands of HSP90 belonging to more than seven chemical series to estimate their relative residence times. We introduce two new techniques for the analysis and the classification of the simulated unbinding trajectories. The first technique, which helps in estimating the limits of the free energy well around the bound state, and the second one, based on a new contact map fingerprint, allow the description and the comparison of the paths that lead to unbinding. Using these analyses, we find that ScaledMD's relative residence time generally enables the identification of the slowest unbinders. We propose an explanation for the underestimation of the residence times of a subset of compounds, and we investigate how the biasing in ScaledMD can affect the mechanistic insights that can be gained from the simulations.

A Blueprint for High Affinity SARS-CoV-2 Mpro Inhibitors from Activity-Based Compound Library Screening Guided by Analysis of Protein Dynamics.

The SARS-CoV-2 coronavirus outbreak continues to spread at a rapid rate worldwide. The main protease (Mpro) is an attractive target for anti-COVID-19 agents. Unexpected difficulties have been encountered in the design of specific inhibitors. Here, by analyzing an ensemble of ∼30 000 SARS-CoV-2 Mpro conformations from crystallographic studies and molecular simulations, we show that small structural variations in the binding site dramatically impact ligand binding properties. Hence, traditional druggability indices fail to adequately discriminate between highly and poorly druggable conformations of the binding site. By performing ∼200 virtual screenings of compound libraries on selected protein structures, we redefine the protein's druggability as the consensus chemical space arising from the multiple conformations of the binding site formed upon ligand binding. This procedure revealed a unique SARS-CoV-2 Mpro blueprint that led to a definition of a specific structure-based pharmacophore. The latter explains the poor transferability of potent SARS-CoV Mpro inhibitors to SARS-CoV-2 Mpro, despite the identical sequences of the active sites. Importantly, application of the pharmacophore predicted novel high affinity inhibitors of SARS-CoV-2 Mpro, that were validated by in vitro assays performed here and by a newly solved X-ray crystal structure. These results provide a strong basis for effective rational drug design campaigns against SARS-CoV-2 Mpro and a new computational approach to screen protein targets with malleable binding sites.

Structure-kinetic relationship reveals the mechanism of selectivity of FAK inhibitors over PYK2.

There is increasing evidence of a significant correlation between prolonged drug-target residence time and increased drug efficacy. Here, we report a structural rationale for kinetic selectivity between two closely related kinases: focal adhesion kinase (FAK) and proline-rich tyrosine kinase 2 (PYK2). We found that slowly dissociating FAK inhibitors induce helical structure at the DFG motif of FAK but not PYK2. Binding kinetic data, high-resolution structures and mutagenesis data support the role of hydrophobic interactions of inhibitors with the DFG-helical region, providing a structural rationale for slow dissociation rates from FAK and kinetic selectivity over PYK2. Our experimental data correlate well with computed relative residence times from molecular simulations, supporting a feasible strategy for rationally optimizing ligand residence times. We suggest that the interplay between the protein structural mobility and ligand-induced effects is a key regulator of the kinetic selectivity of inhibitors of FAK versus PYK2.