Phosphopeptide binding to the N-SH2 domain of tyrosine phosphatase SHP2 correlates with the unzipping of its central ß-sheet.

SHP2 is a tyrosine phosphatase that plays a regulatory role in multiple intracellular signaling cascades and is known to be oncogenic in certain contexts. In the absence of effectors, SHP2 adopts an autoinhibited conformation with its N-SH2 domain blocking the active site. Given the key role of N-SH2 in regulating SHP2, this domain has been extensively studied, often by X-ray crystallography. Using a combination of structural analyses and molecular dynamics (MD) simulations we show that the crystallographic environment can significantly influence the structure of the isolated N-SH2 domain, resulting in misleading interpretations. As an orthogonal method to X-ray crystallography, we use a combination of NMR spectroscopy and MD simulations to accurately determine the conformation of apo N-SH2 in solution. In contrast to earlier reports based on crystallographic data, our results indicate that apo N-SH2 in solution primarily adopts a conformation with a fully zipped central ß-sheet, and that partial unzipping of this ß-sheet is promoted by binding of either phosphopeptides or even phosphate/sulfate ions.

Atomistic ensemble of active SHP2 phosphatase.

SHP2 phosphatase plays an important role in regulating several intracellular signaling pathways. Pathogenic mutations of SHP2 cause developmental disorders and are linked to hematological malignancies and cancer. SHP2 comprises two tandemly-arranged SH2 domains, a catalytic PTP domain, and a disordered C-terminal tail. Under physiological, non-stimulating conditions, the catalytic site of PTP is occluded by the N-SH2 domain, so that the basal activity of SHP2 is low. Whereas the autoinhibited structure of SHP2 has been known for two decades, its active, open structure still represents a conundrum. Since the oncogenic mutant SHP2E76K almost completely populates the active, open state, this mutant has been extensively studied as a model for activated SHP2. By molecular dynamics simulations and accurate explicit-solvent SAXS curve predictions, we present the heterogeneous atomistic ensemble of constitutively active SHP2E76K in solution, encompassing a set of conformational arrangements and radii of gyration in agreement with experimental SAXS data.

Revealing Allostery in PTPN11 SH2 Domains from MD Simulations.

Src-homology 2 (SH2) domains are protein interaction domains that bind to specific peptide motifs containing phosphotyrosine. SHP2, a tyrosine phosphatase encoded by PTPN11 gene, which has been emerged as positive or negative modulator in multiple signaling pathways, contains two SH2 domains, respectively, called N-SH2 and C-SH2. These domains play a relevant role in regulating SHP2 activity, either by recognizing its binding partners or by blocking its catalytic site. Considering the multiple functions that these domains carry out in SHP2, N-SH2 and C-SH2 represent an interesting case of study. In this chapter, we present a methodology that permits, by means of the principal component analysis (PCA), to study and to rationalize the structures adopted by the SH2 domains, in terms of the conformations of their binding sites. The structures can be distinguished, grouped, classified, and reported in a diagram. This approach permits to identify the accessible conformations of the SH2 domains in different binding conditions and to eventually reveal allosteric interactions. The method further reveals that the conformation dynamics of N-SH2 and C-SH2 strongly differ, which likely reflects their distinct functional roles.

Experiment-guided molecular simulations define a heterogeneous structural ensemble for the PTPN11 tandem SH2 domains.

SHP2 plays an important role in regulating cellular processes, and its pathogenic mutations cause developmental disorders and are linked to cancer. SHP2 is a multidomain protein, comprising two SH2 domains arranged in tandem, a catalytic PTP domain, and a disordered C-terminal tail. SHP2 is activated upon binding two linked phosphopeptides to its SH2 domains, and the peptide orientation and spacing between binding sites are critical for enzymatic activation. For decades, the tandem SH2 has been extensively studied to identify the relative orientation of the two SH2 domains that most effectively binds effectors. So far, neither crystallography nor experiments in solution have provided conclusive results. Using experiment-guided molecular simulations, we determine the heterogeneous structural ensemble of the tandem SH2 in solution in agreement with experimental data from small-angle X-ray scattering and NMR residual dipolar couplings. In the solution ensemble, N-SH2 adopts different orientations and positions relative to C-SH2. We suggest that the intrinsic structural plasticity of the tandem SH2 allows SHP2 to respond to external stimuli and is essential for its functional activity.

Probing the Role of Murine Neuroglobin CDloop-D-Helix Unit in CO Ligand Binding and Structural Dynamics.

We produced a neuroglobin variant, namely, Ngb CDless, with the excised CDloop- and D-helix, directly joining the C- and E-helices. The CDless variant retained bis-His hexacoordination, and we investigated the role of the CDloop-D-helix unit in controlling the CO binding and structural dynamics by an integrative approach based on X-ray crystallography, rapid mixing, laser flash photolysis, resonance Raman spectroscopy, and molecular dynamics simulations. Rapid mixing and laser flash photolysis showed that ligand affinity was unchanged with respect to the wild-type protein, albeit with increased on and off constants for rate-limiting heme iron hexacoordination by the distal His64. Accordingly, resonance Raman spectroscopy highlighted a more open distal pocket in the CO complex that, in agreement with MD simulations, likely involves His64 swinging inward and outward of the distal heme pocket. Ngb CDless displays a more rigid overall structure with respect to the wild type, abolishing the structural dynamics of the CDloop-D-helix hypothesized to mediate its signaling role, and it retains ligand binding control by distal His64. In conclusion, this mutant may represent a tool to investigate the involvement of CDloop-D-helix in neuroprotective signaling in a cellular or animal model.

Conformational dynamics promotes disordered regions from function-dispensable to essential in evolved site-specific DNA recombinases.

Protein intrinsically disordered regions (IDRs) play pivotal roles in molecular recognition and regulatory processes through structural disorder-to-order transitions. To understand and exploit the distinctive functional implications of IDRs and to unravel the underlying molecular mechanisms, structural disorder-to-function relationships need to be deciphered. The DNA site-specific recombinase system Cre/loxP represents an attractive model to investigate functional molecular mechanisms of IDRs. Cre contains a functionally dispensable disordered N-terminal tail, which becomes indispensable in the evolved Tre/loxLTR recombinase system. The difficulty to experimentally obtain structural information about this tail has so far precluded any mechanistic study on its involvement in DNA recombination. Here, we use in vitro and in silico evolution data, conformational dynamics, AI-based folding simulations, thermodynamic stability calculations, mutagenesis and DNA recombination assays to investigate how evolution and the dynamic behavior of this IDR may determine distinct functional properties. Our studies suggest that partial conformational order in the N-terminal tail of Tre recombinase and its packing to a conserved hydrophobic surface on the protein provide thermodynamic stability. Based on our results, we propose a link between protein stability and function, offering new plausible atom-detailed mechanistic insights into disorder-function relationships. Our work highlights the potential of N-terminal tails to be exploited for regulation of the activity of Cre-like tyrosine-type SSRs, which merits future investigations and could be of relevance in future rational engineering for their use in biotechnology and genomic medicine.

The loops of the N-SH2 binding cleft do not serve as allosteric switch in SHP2 activation.

The Src-homology-2 domain-containing phosphatase SHP2 is a critical regulator of signal transduction, being implicated in cell growth and differentiation. Activating mutations cause developmental disorders and act as oncogenic drivers in hematologic cancers. SHP2 is activated by phosphopeptide binding to the N-SH2 domain, triggering the release of N-SH2 from the catalytic PTP domain. Based on early crystallographic data, it has been widely accepted that opening of the binding cleft of N-SH2 serves as the key "allosteric switch" driving SHP2 activation. To test the putative coupling between binding cleft opening and SHP2 activation as assumed by the allosteric switch model, we critically reviewed structural data of SHP2, and we used extensive molecular dynamics (MD) simulation and free energy calculations of isolated N-SH2 in solution, SHP2 in solution, and SHP2 in a crystal environment. Our results demonstrate that the binding cleft in N-SH2 is constitutively flexible and open in solution and that a closed cleft found in certain structures is a consequence of crystal contacts. The degree of opening of the binding cleft has only a negligible effect on the free energy of SHP2 activation. Instead, SHP2 activation is greatly favored by the opening of the central ß-sheet of N-SH2. We conclude that opening of the N-SH2 binding cleft is not the key allosteric switch triggering SHP2 activation.

An allosteric interaction controls the activation mechanism of SHP2 tyrosine phosphatase.

SHP2 is a protein tyrosine phosphatase (PTP) involved in multiple signaling pathways. Mutations of SHP2 can result in Noonan syndrome or pediatric malignancies. Inhibition of wild-type SHP2 represents a novel strategy against several cancers. SHP2 is activated by binding of a phosphopeptide to the N-SH2 domain of SHP2, thereby favoring dissociation of the N-SH2 domain and exposing the active site on the PTP domain. The conformational transitions controlling ligand affinity and PTP dissociation remain poorly understood. Using molecular simulations, we revealed an allosteric interaction restraining the N-SH2 domain into a SHP2-activating and a stabilizing state. Only ligands selecting for the activating N-SH2 conformation, depending on ligand sequence and binding mode, are effective activators. We validate the model of SHP2 activation by rationalizing modified basal activity and responsiveness to ligand stimulation of several N-SH2 variants. This study provides mechanistic insight into SHP2 activation and may open routes for SHP2 regulation.

Structural Determinants of Phosphopeptide Binding to the N-Terminal Src Homology 2 Domain of the SHP2 Phosphatase.

SH2 domain-containing tyrosine phosphatase 2 (SHP2), encoded by PTPN11, plays a fundamental role in the modulation of several signaling pathways. Germline and somatic mutations in PTPN11 are associated with different rare diseases and hematologic malignancies, and recent studies have individuated SHP2 as a central node in oncogenesis and cancer drug resistance. The SHP2 structure includes two Src homology 2 domains (N-SH2 and C-SH2) followed by a catalytic protein tyrosine phosphatase (PTP) domain. Under basal conditions, the N-SH2 domain blocks the active site, inhibiting phosphatase activity. Association of the N-SH2 domain with binding partners containing short amino acid motifs comprising a phosphotyrosine residue (pY) leads to N-SH2/PTP dissociation and SHP2 activation. Considering the relevance of SHP2 in signaling and disease and the central role of the N-SH2 domain in its allosteric regulation mechanism, we performed microsecond-long molecular dynamics (MD) simulations of the N-SH2 domain complexed to 12 different peptides to define the structural and dynamical features determining the binding affinity and specificity of the domain. Phosphopeptide residues at position -2 to +5, with respect to pY, have significant interactions with the SH2 domain. In addition to the strong interaction of the pY residue with its conserved binding pocket, the complex is stabilized hydrophobically by insertion of residues +1, +3, and +5 in an apolar groove of the domain and interaction of residue -2 with both the pY and a protein surface residue. Additional interactions are provided by hydrogen bonds formed by the backbone of residues -1, +1, +2, and +4. Finally, negatively charged residues at positions +2 and +4 are involved in electrostatic interactions with two lysines (Lys89 and Lys91) specific for the SHP2 N-SH2 domain. Interestingly, the MD simulations illustrated a previously undescribed conformational flexibility of the domain, involving the core ß sheet and the loop that closes the pY binding pocket.

Pathogenic PTPN11 variants involving the poly-glutamine Gln255 -Gln256 -Gln257 stretch highlight the relevance of helix B in SHP2's functional regulation.

Germline PTPN11 mutations cause Noonan syndrome (NS), the most common disorder among RASopathies. PTPN11 encodes SHP2, a protein tyrosine-phosphatase controlling signaling through the RAS-MAPK and PI3K-AKT pathways. Generally, NS-causing PTPN11 mutations are missense changes destabilizing the inactive conformation of the protein or enhancing its binding to signaling partners. Here, we report on two PTPN11 variants resulting in the deletion or duplication of one of three adjacent glutamine residues (Gln255 -to-Gln257 ). While p.(Gln257dup) caused a typical NS phenotype in carriers of a first family, p.(Gln257del) had incomplete penetrance in a second family. Missense mutations involving Gln256 had previously been reported in NS. This poly-glutamine stretch is located on helix B of the PTP domain, a region involved in stabilizing SHP2 in its autoinhibited state. Molecular dynamics simulations predicted that changes affecting this motif perturb the SHP2's catalytically inactive conformation and/or substrate recognition. Biochemical data showed that duplication and deletion of Gln257 variably enhance SHP2's catalytic activity, while missense changes involving Gln256 affect substrate specificity. Expression of mutants in HEK293T cells documented their activating role on MAPK signaling, uncoupling catalytic activity and modulation of intracellular signaling. These findings further document the relevance of helix B in the regulation of SHP2's function.

Proximal and distal control for ligand binding in neuroglobin: role of the CD loop and evidence for His64 gating.

Neuroglobin (Ngb) is predominantly expressed in neurons of the central and peripheral nervous systems and it clearly seems to be involved in neuroprotection. Engineering Ngb to observe structural and dynamic alterations associated with perturbation in ligand binding might reveal important structural determinants, and could shed light on key features related to its mechanism of action. Our results highlight the relevance of the CD loop and of Phe106 as distal and proximal controls involved in ligand binding in murine neuroglobin. We observed the effects of individual and combined mutations of the CD loop and Phe106 that conferred to Ngb higher CO binding velocities, which we correlate with the following structural observations: the mutant F106A shows, upon CO binding, a reduced heme sliding hindrance, with the heme present in a peculiar double conformation, whereas in the CD loop mutant "Gly-loop", the original network of interactions between the loop and the heme was abolished, enhancing binding via facilitated gating out of the distal His64. Finally, the double mutant, combining both mutations, showed a synergistic effect on CO binding rates. Resonance Raman spectroscopy and MD simulations support our findings on structural dynamics and heme interactions in wild type and mutated Ngbs.

Dynamic multiple-scattering treatment of X-ray absorption: Parameterization of a new molecular dynamics force field for myoglobin.

We present a detailed analysis of the X-ray absorption near-edge structure (XANES) data on the Fe K-edge of CO Myoglobin based on a combined procedure of Molecular Dynamics (MD) calculations and MXAN (Minuit XANes) data analysis that we call D-MXAN. The ability of performing quantitative XANES data analysis allows us to refine classical force field MD parameters, thus obtaining a reliable tool for the atomic investigation of this important model system for biological macromolecules. The iterative procedure here applied corrects the greatest part of the structural discrepancy between classical MD sampling and experimental determinations. Our procedure, moreover, is able to discriminate between different heme conformational basins visited during the MD simulation, thus demonstrating the necessity of a sampling on the order of tens of nanoseconds, even for an application such X-ray absorption spectroscopy data analysis.

Clinical and functional characterization of two novel ZBTB20 mutations causing Primrose syndrome.

Primrose syndrome (PS) is a rare disorder characterized by macrocephaly, tall stature, intellectual disability, autistic traits, and disturbances of glucose metabolism with insulin-resistant diabetes and distal muscle wasting occurring in adulthood. The disorder is caused by functional dysregulation of ZBTB20, a transcriptional repressor controlling energetic metabolism and developmental programs. ZBTB20 maps in a genomic region that is deleted in the 3q13.31 microdeletion syndrome, which explains the clinical overlap between the two disorders. A narrow spectrum of amino acid substitutions in a restricted region of ZBTB20 encompassing the first and second zinc-finger motifs have been reported thus far. Here, we characterize clinically and functionally the first truncating mutation [(c.1024delC; p.(Gln342Serfs*42)] and a missense change affecting the third zinc-finger motif of the protein [(c.1931C > T; p.(Thr644Ile)]. Our data document that both mutations have dominant negative impact on wild-type ZBTB20, providing further evidence of the specific behavior of PS-causing mutations on ZBTB20 function.

Structural, Functional, and Clinical Characterization of a Novel PTPN11 Mutation Cluster Underlying Noonan Syndrome.

Germline mutations in PTPN11, the gene encoding the Src-homology 2 (SH2) domain-containing protein tyrosine phosphatase (SHP2), cause Noonan syndrome (NS), a relatively common, clinically variable, multisystem disorder. Here, we report on the identification of five different PTPN11 missense changes affecting residues Leu261 , Leu262 , and Arg265 in 16 unrelated individuals with clinical diagnosis of NS or with features suggestive for this disorder, specifying a novel disease-causing mutation cluster. Expression of the mutant proteins in HEK293T cells documented their activating role on MAPK signaling. Structural data predicted a gain-of-function role of substitutions at residues Leu262 and Arg265 exerted by disruption of the N-SH2/PTP autoinhibitory interaction. Molecular dynamics simulations suggested a more complex behavior for changes affecting Leu261 , with possible impact on SHP2's catalytic activity/selectivity and proper interaction of the PTP domain with the regulatory SH2 domains. Consistent with that, biochemical data indicated that substitutions at codons 262 and 265 increased the catalytic activity of the phosphatase, while those affecting codon 261 were only moderately activating but impacted substrate specificity. Remarkably, these mutations underlie a relatively mild form of NS characterized by low prevalence of cardiac defects, short stature, and cognitive and behavioral issues, as well as less evident typical facial features.

Exploring Multiple Binding Modes Using Confined Replica Exchange Molecular Dynamics.

Molecular docking is extensively applied to determine the position of a ligand on its receptor despite the rather poor correspondence between docking scores and experimental binding affinities found in several studies, especially for systems structurally unrelated with those used in the scoring functions' training sets. Here, we present a method for the prediction of binding modes and binding free energies, which uses replica exchange molecular dynamics in combination with a receptor-shaped piecewise potential, confining the ligand in the proximity of the receptor surface and limiting the accessible conformational space of interest. We assess our methodology with a set of protein receptor-ligand test cases. In every case studied, the method is able to locate the ligand on the experimentally known receptor binding site, and it gives as output the binding free energy. The added value of our approach with respect to other available methods is that it quickly performs a conformational space search, providing a set of bound (or unbound) configurations, which can be used to determine phenomenological structural and energetic properties of an experimental binding state as a result of contributions provided by diversified multiple binding poses.

Mutations Impairing GSK3-Mediated MAF Phosphorylation Cause Cataract, Deafness, Intellectual Disability, Seizures, and a Down Syndrome-like Facies.

Transcription factors operate in developmental processes to mediate inductive events and cell competence, and perturbation of their function or regulation can dramatically affect morphogenesis, organogenesis, and growth. We report that a narrow spectrum of amino-acid substitutions within the transactivation domain of the v-maf avian musculoaponeurotic fibrosarcoma oncogene homolog (MAF), a leucine zipper-containing transcription factor of the AP1 superfamily, profoundly affect development. Seven different de novo missense mutations involving conserved residues of the four GSK3 phosphorylation motifs were identified in eight unrelated individuals. The distinctive clinical phenotype, for which we propose the eponym Aymé-Gripp syndrome, is not limited to lens and eye defects as previously reported for MAF/Maf loss of function but includes sensorineural deafness, intellectual disability, seizures, brachycephaly, distinctive flat facial appearance, skeletal anomalies, mammary gland hypoplasia, and reduced growth. Disease-causing mutations were demonstrated to impair proper MAF phosphorylation, ubiquitination and proteasomal degradation, perturbed gene expression in primary skin fibroblasts, and induced neurodevelopmental defects in an in vivo model. Our findings nosologically and clinically delineate a previously poorly understood recognizable multisystem disorder, provide evidence for MAF governing a wider range of developmental programs than previously appreciated, and describe a novel instance of protein dosage effect severely perturbing development.

Mutations in ZBTB20 cause Primrose syndrome.

Primrose syndrome and 3q13.31 microdeletion syndrome are clinically related disorders characterized by tall stature, macrocephaly, intellectual disability, disturbed behavior and unusual facial features, with diabetes, deafness, progressive muscle wasting and ectopic calcifications specifically occurring in the former. We report that missense mutations in ZBTB20, residing within the 3q13.31 microdeletion syndrome critical region, underlie Primrose syndrome. This finding establishes a genetic link between these disorders and delineates the impact of ZBTB20 dysregulation on development, growth and metabolism.

Sulfated glycosaminoglycans exploit the conformational plasticity of bone morphogenetic protein-2 (BMP-2) and alter the interaction profile with its receptor.

Sulfated glycosaminoglycans (GAGs) can direct cellular processes by interacting with proteins of the extracellular matrix (ECM). In this study we characterize the interaction profiles of chemically sulfated hyaluronan (HA) and chondroitin sulfate (CS) derivatives with bone morphogenetic protein-2 (BMP-2) and investigate their relevance for complex formation with the receptor BMPR-IA. These goals were addressed by surface plasmon resonance (SPR) and ELISA in combination with molecular modeling and dynamics simulation. We found not only the interaction of BMP-2 with GAGs to be dependent on the type and sulfation of GAGs but also BMP-2/GAG/BMPR-IA complex formation. The conformational plasticity of the BMP-2 N-termini plays a key role in the structural and thermodynamic characteristics of the BMP-2/GAG/BMPR-IA system. Hence we propose a model that provides direct insights into the importance of the structural and dynamical properties of the BMP-2/BMPR-IA system for its regulation by sulfated GAGs, in which structural asymmetry plays a key role.

Analysis of infrared spectra of ß-hairpin peptides as derived from molecular dynamics simulations.

Infrared temperature-dependent spectroscopy is a well-known tool to characterize folding/unfolding transitions in peptides and proteins, assuming that the higher the temperature, the higher the unfolded population. The infrared spectra at different temperatures of two ß-hairpin peptides (gramicidin S analogues GS6 and GS10) are here reconstructed by means of molecular dynamics (MD) simulations and a theoretical-computational method based on the perturbed matrix method. The calculated temperature-dependent spectra result in good agreement with the experimental available spectra. The same methodology has been then used to reconstruct the spectra corresponding to the pure unfolded and folded states, as defined from the MD simulations, in order to better understand the temperature-dependent spectra and to help the interpretation of the experimental spectra. For example, our results show that in the case of the GS6 peptide the analysis of the temperature-dependent spectra cannot be used to investigate the folding/unfolding kinetics within the usual assumption that the higher the temperature, the higher the probability of the unfolded state.

Structure of the lipodepsipeptide syringomycin E in phospholipids and sodium dodecylsulphate micelle studied by circular dichroism, NMR spectroscopy and molecular dynamics.

Syringomycin E (SRE) is a member of a family of lipodepsipeptides that characterize the secondary metabolism of the plant-associated bacteria Pseudomonas syringae pv. syringae. It displays phytotoxic, antifungal and haemolytic activities, due to the membrane interaction and ion channel formation. To gain an insight into the conformation of SRE in the membrane environment, we studied the conformation of SRE bound to SDS micelle, a suitable model for the membrane-bound SRE. In fact, highly similar circular dichroism (CD) spectra were obtained for SRE bound to sodium dodecylsulphate (SDS) and to a phospholipid bilayer, indicating the conformational equivalence of SRE in these two media, at difference with the CD spectrum of SRE in water solution. The structure of SDS-bound SRE was determined by NMR spectroscopy combined with molecular dynamics calculations in octane environment. The results of this study highlight the influence of the interaction with lipids in determining the three-dimensional structure of SRE and provide the basis for further investigations on structural determinants of syringomycin E-membrane interaction.