Collective Dynamics Differentiates Functional Divergence in Protein Evolution

Protein evolution is most commonly studied by analyzing related protein sequences and generating ancestral sequences through Bayesian and Maximum Likelihood methods, and/or by resurrecting ancestral proteins in the lab and performing ligand binding studies to determine function. Structural and dynamic evolution have largely been left out of molecular evolution studies. Here we incorporate both structure and dynamics to elucidate the molecular principles behind the divergence in the evolutionary path of the steroid receptor proteins. We determine the likely structure of three evolutionarily diverged ancestral steroid receptor proteins using the Zipping and Assembly Method with FRODA (ZAMF). Our predictions are within ∼2.7 Å all-atom RMSD of the respective crystal structures of the ancestral steroid receptors. Beyond static structure prediction, a particular feature of ZAMF is that it generates protein dynamics information. We investigate the differences in conformational dynamics of diverged proteins by obtaining the most collective motion through essential dynamics. Strikingly, our analysis shows that evolutionarily diverged proteins of the same family do not share the same dynamic subspace, while those sharing the same function are simultaneously clustered together and distant from those, that have functionally diverged. Dynamic analysis also enables those mutations that most affect dynamics to be identified. It correctly predicts all mutations (functional and permissive) necessary to evolve new function and ∼60% of permissive mutations necessary to recover ancestral function

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

SH3 Domain-Peptide Binding Energy Calculations Based on Structural Ensemble and Multiple Peptide Templates

SH3 domains mediate signal transduction by recognizing short peptides. Understanding of the driving forces in peptide recognitions will help us to predict the binding specificity of the domain-peptide recognition and to understand the molecular interaction networks of cells. However, accurate calculation of the binding energy is a tough challenge. In this study, we propose three ideas for improving our ability to predict the binding energy between SH3 domains and peptides: (1) utilizing the structural ensembles sampled from a molecular dynamics simulation trajectory, (2) utilizing multiple peptide templates, and (3) optimizing the sequence-structure mapping. We tested these three ideas on ten previously studied SH3 domains for which SPOT analysis data were available. The results indicate that calculating binding energy using the structural ensemble was most effective, clearly increasing the prediction accuracy, while the second and third ideas tended to give better binding energy predictions. We applied our method to the five SH3 targets in DREAM4 Challenge and selected the best performing method

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