Change in Allosteric Network Affects Binding Affinities of PDZ Domains: Analysis through Perturbation Response Scanning

The allosteric mechanism plays a key role in cellular functions of several PDZ domain proteins (PDZs) and is directly linked to pharmaceutical applications; however, it is a challenge to elaborate the nature and extent of these allosteric interactions. One solution to this problem is to explore the dynamics of PDZs, which may provide insights about how intramolecular communication occurs within a single domain. Here, we develop an advancement of perturbation response scanning (PRS) that couples elastic network models with linear response theory (LRT) to predict key residues in allosteric transitions of the two most studied PDZs (PSD-95 PDZ3 domain and hPTP1E PDZ2 domain). With PRS, we first identify the residues that give the highest mean square fluctuation response upon perturbing the binding sites. Strikingly, we observe that the residues with the highest mean square fluctuation response agree with experimentally determined residues involved in allosteric transitions. Second, we construct the allosteric pathways by linking the residues giving the same directional response upon perturbation of the binding sites. The predicted intramolecular communication pathways reveal that PSD-95 and hPTP1E have different pathways through the dynamic coupling of different residue pairs. Moreover, our analysis provides a molecular understanding of experimentally observed hidden allostery of PSD-95. We show that removing the distal third alpha helix from the binding site alters the allosteric pathway and decreases the binding affinity. Overall, these results indicate that (i) dynamics plays a key role in allosteric regulations of PDZs, (ii) the local changes in the residue interactions can lead to significant changes in the dynamics of allosteric regulations, and (iii) this might be the mechanism that each PDZ uses to tailor their binding specificities regulation

Beyond the Binding Site: The Role of the β2 – β3 Loop and Extra-Domain Structures in PDZ Domains

A general paradigm to understand protein function is to look at properties of isolated well conserved domains, such as SH3 or PDZ domains. While common features of domain families are well understood, the role of subtle differences among members of these families is less clear. Here, molecular dynamics simulations indicate that the binding mechanism in PSD95-PDZ3 is critically regulated via interactions outside the canonical binding site, involving both the poorly conserved loop and an extra-domain helix. Using the CRIPT peptide as a prototypical ligand, our simulations suggest that a network of salt-bridges between the ligand and this loop is necessary for binding. These contacts interconvert between each other on a time scale of a few tens of nanoseconds, making them elusive to X-ray crystallography. The loop is stabilized by an extra-domain helix. The latter influences the global dynamics of the domain, considerably increasing binding affinity. We found that two key contacts between the helix and the domain, one involving the loop, provide an atomistic interpretation of the increased affinity. Our analysis indicates that both extra-domain segments and loosely conserved regions play critical roles in PDZ binding affinity and specificity

Detailed Regulatory Mechanism of the Interaction between ZO-1 PDZ2 and Connexin43 Revealed by MD Simulations

The gap junction protein connexin43 (Cx43) binds to the second PDZ domain of Zonula occludens-1 (ZO-1) through its C-terminal tail, mediating the regulation of gap junction plaque size and dynamics. Biochemical study demonstrated that the very C-terminal 12 residues of Cx43 are necessary and sufficient for ZO-1 PDZ2 binding and phosphorylation at residues Ser (-9) and Ser (-10) of the peptide can disrupt the association. However, only a crystal structure of ZO-1 PDZ2 in complex with a shorter 9 aa peptide of connexin43 was solved experimentally. Here, the interactions between ZO-1 PDZ2 and the short, long and phosphorylated Cx43 peptides were studied using molecular dynamics (MD) simulations and free energy calculation. The short peptide bound to PDZ2 exhibits large structural variations, while the extension of three upstream residues stabilizes the peptide conformation and enhanced the interaction. Phosphorylation at Ser(-9) significantly weakens the binding and results in conformational flexibility of the peptide. Glu210 of ZO-1 PDZ2 was found to be a key regulatory point in Cx43 binding and phosphorylation induced dissociation

Binding Free Energy Landscape of Domain-Peptide Interactions

Peptide recognition domains (PRDs) are ubiquitous protein domains which mediate large numbers of protein interactions in the cell. How these PRDs are able to recognize peptide sequences in a rapid and specific manner is incompletely understood. We explore the peptide binding process of PDZ domains, a large PRD family, from an equilibrium perspective using an all-atom Monte Carlo (MC) approach. Our focus is two different PDZ domains representing two major PDZ classes, I and II. For both domains, a binding free energy surface with a strong bias toward the native bound state is found. Moreover, both domains exhibit a binding process in which the peptides are mostly either bound at the PDZ binding pocket or else interact little with the domain surface. Consistent with this, various binding observables show a temperature dependence well described by a simple two-state model. We also find important differences in the details between the two domains. While both domains exhibit well-defined binding free energy barriers, the class I barrier is significantly weaker than the one for class II. To probe this issue further, we apply our method to a PDZ domain with dual specificity for class I and II peptides, and find an analogous difference in their binding free energy barriers. Lastly, we perform a large number of fixed-temperature MC kinetics trajectories under binding conditions. These trajectories reveal significantly slower binding dynamics for the class II domain relative to class I. Our combined results are consistent with a binding mechanism in which the peptide C terminal residue binds in an initial, rate-limiting step

Predictive Power of Molecular Dynamics Receptor Structures in Virtual Screening

Molecular dynamics (MD) simulation is a well-established method for understanding protein dynamics. Conformations from unrestrained MD simulations have yet to be assessed for blind virtual screening (VS) by docking. This study presents a critical analysis of the predictive power of MD snapshots to this regard, evaluating two well-characterized systems of varying flexibility in ligand-bound and unbound configurations. Results from such VS predictions are discussed with respect to experimentally determined structures. In all cases, MD simulations provide snapshots that improve VS predictive power over known crystal structures, possibly due to sampling more relevant receptor conformations. Additionally, MD can move conformations previously not amenable to docking into the predictive range

Effect of the COVID-19 pandemic on surgery for indeterminate thyroid nodules (THYCOVID): a retrospective, international, multicentre, cross-sectional study

Background Since its outbreak in early 2020, the COVID-19 pandemic has diverted resources from non-urgent and elective procedures, leading to diagnosis and treatment delays, with an increased number of neoplasms at advanced stages worldwide. The aims of this study were to quantify the reduction in surgical activity for indeterminate thyroid nodules during the COVID-19 pandemic; and to evaluate whether delays in surgery led to an increased occurrence of aggressive tumours.Methods In this retrospective, international, cross-sectional study, centres were invited to participate in June 22, 2022; each centre joining the study was asked to provide data from medical records on all surgical thyroidectomies consecutively performed from Jan 1, 2019, to Dec 31, 2021. Patients with indeterminate thyroid nodules were divided into three groups according to when they underwent surgery: from Jan 1, 2019, to Feb 29, 2020 (global prepandemic phase), from March 1, 2020, to May 31, 2021 (pandemic escalation phase), and from June 1 to Dec 31, 2021 (pandemic decrease phase). The main outcomes were, for each phase, the number of surgeries for indeterminate thyroid nodules, and in patients with a postoperative diagnosis of thyroid cancers, the occurrence of tumours larger than 10 mm, extrathyroidal extension, lymph node metastases, vascular invasion, distant metastases, and tumours at high risk of structural disease recurrence. Univariate analysis was used to compare the probability of aggressive thyroid features between the first and third study phases. The study was registered on ClinicalTrials.gov, NCT05178186.Findings Data from 157 centres (n=49 countries) on 87 467 patients who underwent surgery for benign and malignant thyroid disease were collected, of whom 22 974 patients (18 052 [78 center dot 6%] female patients and 4922 [21 center dot 4%] male patients) received surgery for indeterminate thyroid nodules. We observed a significant reduction in surgery for indeterminate thyroid nodules during the pandemic escalation phase (median monthly surgeries per centre, 1 center dot 4 [IQR 0 center dot 6-3 center dot 4]) compared with the prepandemic phase (2 center dot 0 [0 center dot 9-3 center dot 7]; p<0 center dot 0001) and pandemic decrease phase (2 center dot 3 [1 center dot 0-5 center dot 0]; p<0 center dot 0001). Compared with the prepandemic phase, in the pandemic decrease phase we observed an increased occurrence of thyroid tumours larger than 10 mm (2554 [69 center dot 0%] of 3704 vs 1515 [71 center dot 5%] of 2119; OR 1 center dot 1 [95% CI 1 center dot 0-1 center dot 3]; p=0 center dot 042), lymph node metastases (343 [9 center dot 3%] vs 264 [12 center dot 5%]; OR 1 center dot 4 [1 center dot 2-1 center dot 7]; p=0 center dot 0001), and tumours at high risk of structural disease recurrence (203 [5 center dot 7%] of 3584 vs 155 [7 center dot 7%] of 2006; OR 1 center dot 4 [1 center dot 1-1 center dot 7]; p=0 center dot 0039).Interpretation Our study suggests that the reduction in surgical activity for indeterminate thyroid nodules during the COVID-19 pandemic period could have led to an increased occurrence of aggressive thyroid tumours. However, other compelling hypotheses, including increased selection of patients with aggressive malignancies during this period, should be considered. We suggest that surgery for indeterminate thyroid nodules should no longer be postponed even in future instances of pandemic escalation.Funding None.Copyright (c) 2023 Published by Elsevier Ltd. All rights reserved

Manipulation of conformational change in proteins by single-residue perturbations

Using the perturbation-response scanning (PRS) technique, we study a set of 25 proteins that display a variety of conformational motions upon ligand binding (e.g., shear, hinge, allosteric). In most cases, PRS determines single residues that may be manipulated to achieve the resulting conformational change. PRS reveals that for some proteins, binding-induced conformational change may be achieved through the perturbation of residues scattered throughout the protein, whereas in others, perturbation of specific residues confined to a highly specific region is necessary. Overlaps between the experimental and PRS-calculated atomic displacement vectors are usually more descriptive of the conformational change than those obtained from a modal analysis of elastic network models. Furthermore, the largest overlaps obtained by the latter approach do not always appear in the most collective modes; there are cases where more than one mode yields comparable overlap sizes. We show that success of the modal analysis depends on an absence of redundant paths in the protein. PRS thus demonstrates that several relevant modes can be induced simultaneously by perturbing a single select residue on the protein. We also illustrate the biological relevance of applying PRS to the GroEL, adenylate kinase, myosin, and kinesin structures in detail by showing that the residues whose perturbation leads to precise conformational changes usually correspond to those experimentally determined to be functionally important

Plant-expressed cocaine hydrolase variants of butyrylcholinesterase exhibit altered allosteric effects of cholinesterase activity and increased inhibitor sensitivity

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper