Pliability in the m6A-Binding Region Extends Druggability of YTH Domains.

Epitranscriptomic mRNA modifications affect gene expression, with their altered balance detected in various cancers. YTHDF proteins contain the YTH reader domain recognizing the m6A mark on mRNA and represent valuable drug targets. Crystallographic structures have been determined for all three family members; however, discrepancies are present in the organization of the m6A-binding pocket. Here, we present new crystallographic structures of the YTH domain of YTHDF1, accompanied by computational studies, showing that this domain can exist in different stable conformations separated by a significant energetic barrier. During the transition, additional conformations are explored, with peculiar druggable pockets appearing and offering new opportunities for the design of YTH-interfering small molecules.

Probing the conformational dynamics of an Ago-RNA complex in water/methanol solution.

Argonaute (Ago) proteins mediate target recognition guiding miRNA to bind complementary mRNA primarily in the seed region. However, additional pairing can occur beyond the seed, forming a supplementary duplex that can contribute to the guide-target affinity. In order to shed light on the connection, between protein-RNA interactions and miRNA-mRNA seed and supplementary duplex mobility, we carried out molecular dynamics simulations at the microsecond time-scale using a different approach compared to the ones normally used. Until now, theoretical investigations with classical MD on Ago-RNA complexes have been focused primarily on pure water solvent, which mimics the natural environment of biological molecules. Here, we explored the conformational space of a human Ago2 (hAgo2) bound to the seed + supplementary miRNA-mRNA duplex, using the solvent environment as a molecular probe. MD simulations have been performed in a mixture of water/MeOH at a molar ratio of 70 : 30 as well as in pure water for comparison. Our findings revealed that the mixed solvent promotes protein RNA association, principally enhancing salt-linkages between basic amino acid side-chains and acidic phosphates of the sugar-phosphate backbone. The primary effect registered was the restriction of supplementary duplex flexibility and the stabilization of the miRNA 3' terminus. Interestingly, we observed that the influence of the solvent appears to have almost no impact on the conformation of the seed duplex.

Role of the membrane anchor in the regulation of Lck activity.

Theoretical work suggests that collective spatiotemporal behavior of integral membrane proteins should be modulated by boundary lipids sheathing their membrane anchors. Here, we show evidence for this prediction while investigating the mechanism for maintaining a steady amount of the active form of integral membrane protein Lck kinase (LckA) by Lck trans-autophosphorylation regulated by the phosphatase CD45. We used super-resolution microscopy, flow cytometry, and pharmacological and genetic perturbation to gain insight into the spatiotemporal context of this process. We found that LckA is generated exclusively at the plasma membrane, where CD45 maintains it in a ceaseless dynamic equilibrium with its unphosphorylated precursor. Steady LckA shows linear dependence, after an initial threshold, over a considerable range of Lck expression levels. This behavior fits a phenomenological model of trans-autophosphorylation that becomes more efficient with increasing LckA. We then challenged steady LckA formation by genetically swapping the Lck membrane anchor with structurally divergent ones, such as that of Src or the transmembrane domains of LAT, CD4, palmitoylation-defective CD4 and CD45 that were expected to drastically modify Lck boundary lipids. We observed small but significant changes in LckA generation, except for the CD45 transmembrane domain that drastically reduced LckA due to its excessive lateral proximity to CD45. Comprehensively, LckA formation and maintenance can be best explained by lipid bilayer critical density fluctuations rather than liquid-ordered phase-separated nanodomains, as previously thought, with "like/unlike" boundary lipids driving dynamical proximity and remoteness of Lck with itself and with CD45.

Conformational response to ligand binding of TMPRSS2, a protease involved in SARS-CoV-2 infection: Insights through computational modeling.

Thanks to the considerable research which has been undertaken in the last few years to improve our understanding of the biology and mechanism of action of SARS-CoV-2, we know how the virus uses its surface spike protein to infect host cells. The transmembrane prosthesis, serine 2 (TMPRSS2) protein, located on the surface of human cells, recognizes the cleavage site in the spike protein, leading to the release of the fusion peptide and entry of the virus into the host cells. Because of its role, TMPRSS2 has been proposed as a drug target to prevent infection by the virus. In this study, we aim to increase our understanding of TMPRSS2 using long scale microsecond atomistic molecular dynamics simulations, focusing on the conformational changes over time. The comparison between simulations conducted on the protein in the native (apo) and inhibited form (holo), has shown that in the holo form the inhibitor stabilizes the catalytic site and induces rearrangements in the extracellular domain of the protein. In turn, it leads to the formation of a new cavity in the vicinity of the ligand binding pocket that is stable in the microsecond time scale. Given the low specificity of known protease inhibitors, these findings suggest a new potential drug target site that can be used to improve TMPRSS2 specific recognition by newly designed inhibitors.

Computational Study of Helicase from SARS-CoV-2 in RNA-Free and Engaged Form.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the pandemic that broke out in 2020 and continues to be the cause of massive global upheaval. Coronaviruses are positive-strand RNA viruses with a genome of ~30 kb. The genome is replicated and transcribed by RNA-dependent RNA polymerase together with accessory factors. One of the latter is the protein helicase (NSP13), which is essential for viral replication. The recently solved helicase structure revealed a tertiary structure composed of five domains. Here, we investigated NSP13 from a structural point of view, comparing its RNA-free form with the RNA-engaged form by using atomistic molecular dynamics (MD) simulations at the microsecond timescale. Structural analyses revealed conformational changes that provide insights into the contribution of the different domains, identifying the residues responsible for domain-domain interactions in both observed forms. The RNA-free system appears to be more flexible than the RNA-engaged form. This result underlies the stabilizing role of the nucleic acid and the functional core role of these domains.

Altered Local Interactions and Long-Range Communications in UK Variant (B.1.1.7) Spike Glycoprotein.

The COVID-19 pandemic is caused by SARS-CoV-2. Currently, most of the research efforts towards the development of vaccines and antibodies against SARS-CoV-2 were mainly focused on the spike (S) protein, which mediates virus entry into the host cell by binding to ACE2. As the virus SARS-CoV-2 continues to spread globally, variants have emerged, characterized by multiple mutations of the S glycoprotein. Herein, we employed microsecond-long molecular dynamics simulations to study the impact of the mutations of the S glycoprotein in SARS-CoV-2 Variant of Concern 202012/01 (B.1.1.7), termed the "UK variant", in comparison with the wild type, with the aim to decipher the structural basis of the reported increased infectivity and virulence. The simulations provided insights on the different dynamics of UK and wild-type S glycoprotein, regarding in particular the Receptor Binding Domain (RBD). In addition, we investigated the role of glycans in modulating the conformational transitions of the RBD. The overall results showed that the UK mutant experiences higher flexibility in the RBD with respect to wild type; this behavior might be correlated with the increased transmission reported for this variant. Our work also adds useful structural information on antigenic "hotspots" and epitopes targeted by neutralizing antibodies.

Molecular dynamics recipes for genome research.

Molecular dynamics (MD) simulation allows one to predict the time evolution of a system of interacting particles. It is widely used in physics, chemistry and biology to address specific questions about the structural properties and dynamical mechanisms of model systems. MD earned a great success in genome research, as it proved to be beneficial in sorting pathogenic from neutral genomic mutations. Considering their computational requirements, simulations are commonly performed on HPC computing devices, which are generally expensive and hard to administer. However, variables like the software tool used for modeling and simulation or the size of the molecule under investigation might make one hardware type or configuration more advantageous than another or even make the commodity hardware definitely suitable for MD studies. This work aims to shed lights on this aspect.

Revisiting the "satisfaction of spatial restraints" approach of MODELLER for protein homology modeling.

The most frequently used approach for protein structure prediction is currently homology modeling. The 3D model building phase of this methodology is critical for obtaining an accurate and biologically useful prediction. The most widely employed tool to perform this task is MODELLER. This program implements the "modeling by satisfaction of spatial restraints" strategy and its core algorithm has not been altered significantly since the early 1990s. In this work, we have explored the idea of modifying MODELLER with two effective, yet computationally light strategies to improve its 3D modeling performance. Firstly, we have investigated how the level of accuracy in the estimation of structural variability between a target protein and its templates in the form of σ values profoundly influences 3D modeling. We show that the σ values produced by MODELLER are on average weakly correlated to the true level of structural divergence between target-template pairs and that increasing this correlation greatly improves the program's predictions, especially in multiple-template modeling. Secondly, we have inquired into how the incorporation of statistical potential terms (such as the DOPE potential) in the MODELLER's objective function impacts positively 3D modeling quality by providing a small but consistent improvement in metrics such as GDT-HA and lDDT and a large increase in stereochemical quality. Python modules to harness this second strategy are freely available at https://github.com/pymodproject/altmod. In summary, we show that there is a large room for improving MODELLER in terms of 3D modeling quality and we propose strategies that could be pursued in order to further increase its performance.

Computational Studies of SARS-CoV-2 3CLpro: Insights from MD Simulations.

Given the enormous social and health impact of the pandemic triggered by severe acute respiratory syndrome 2 (SARS-CoV-2), the scientific community made a huge effort to provide an immediate response to the challenges posed by Coronavirus disease 2019 (COVID-19). One of the most important proteins of the virus is an enzyme, called 3CLpro or main protease, already identified as an important pharmacological target also in SARS and Middle East respiratory syndrome virus (MERS) viruses. This protein triggers the production of a whole series of enzymes necessary for the virus to carry out its replicating and infectious activities. Therefore, it is crucial to gain a deeper understanding of 3CLpro structure and function in order to effectively target this enzyme. All-atoms molecular dynamics (MD) simulations were performed to examine the different conformational behaviors of the monomeric and dimeric form of SARS-CoV-2 3CLpro apo structure, as revealed by microsecond time scale MD simulations. Our results also shed light on the conformational dynamics of the loop regions at the entry of the catalytic site. Studying, at atomic level, the characteristics of the active site and obtaining information on how the protein can interact with its substrates will allow the design of molecules able to block the enzymatic function crucial for the virus.

Cysteine 180 Is a Redox Sensor Modulating the Activity of Human Pyridoxal 5'-Phosphate Histidine Decarboxylase.

Histidine decarboxylase is a pyridoxal 5'-phosphate enzyme catalyzing the conversion of histidine to histamine, a bioactive molecule exerting its role in many modulatory processes. The human enzyme is involved in many physiological functions, such as neurotransmission, gastrointestinal track function, cell growth, and differentiation. Here, we studied the functional properties of the human enzyme and, in particular, the effects exerted at the protein level by two cysteine residues: Cys-180 and Cys-418. Surprisingly, the enzyme exists in an equilibrium between a reduced and an oxidized form whose extent depends on the redox state of Cys-180. Moreover, we determined that (i) the two enzymatic redox species exhibit modest structural changes in the coenzyme microenvironment and (ii) the oxidized form is slightly more active and stable than the reduced one. These data are consistent with the model proposed by bioinformatics analyses and molecular dynamics simulations in which the Cys-180 redox state could be responsible for a structural transition affecting the C-terminal domain reorientation leading to active site alterations. Furthermore, the biochemical properties of the purified C180S and C418S variants reveal that C180S behaves like the reduced form of the wild-type enzyme, while C418S is sensitive to reductants like the wild-type enzyme, thus allowing the identification of Cys-180 as the redox sensitive switch. On the other hand, Cys-418 appears to be a residue involved in aggregation propensity. A possible role for Cys-180 as a regulatory switch in response to different cellular redox conditions could be suggested.

Lethal digenic mutations in the K+ channels Kir4.1 (KCNJ10) and SLACK (KCNT1) associated with severe-disabling seizures and neurodevelopmental delay.

A 2-yr-old boy presented profound developmental delay, failure to thrive, ataxia, hypotonia, and tonic-clonic seizures that caused the death of the patient. Targeted and whole exome sequencing revealed two heterozygous missense variants: a novel mutation in the KCNJ10 gene that encodes for the inward-rectifying K+ channel Kir4.1 and another previously characterized mutation in KCNT1 that encodes for the Na+-activated K+ channel known as Slo2.2 or SLACK. The objectives of this study were to perform the clinical and genetic characterization of the proband and his family and to examine the functional consequence of the Kir4.1 mutation. The mutant and wild-type KCNJ10 constructs were generated and heterologously expressed in Xenopus laevis oocytes, and whole cell K+ currents were measured using the two-electrode voltage-clamp technique. The KCNJ10 mutation c.652C>T resulted in a p.L218F substitution at a highly conserved residue site. Wild-type KCNJ10 expression yielded robust Kir current, whereas currents from oocytes expressing the mutation were reduced, remarkably. Western Blot analysis revealed reduced protein expression by the mutation. Kir5.1 subunits display selective heteromultimerization with Kir4.1 constituting channels with unique kinetics. The effect of the mutation on Kir4.1/5.1 channel activity was twofold: a reduction in current amplitudes and an increase in the pH-dependent inhibition. We thus report a novel loss-of-function mutation in Kir4.1 found in a patient with a coexisting mutation in SLACK channels that results in a fatal disease.NEW & NOTEWORTHY We present and characterize a novel mutation in KCNJ10 Unlike previously reported EAST/SeSAME patients, our patient was heterozygous, and contrary to previous studies, mimicking the heterozygous state by coexpression resulted in loss of channel function. We report in the same patient co-occurrence of a KCNT1 mutation resulting in a more severe phenotype. This study provides new insights into the phenotypic spectrum and to the genotype-phenotype correlations associated with EAST/SeSAME and MMFSI.

Thermal compaction of the intrinsically disordered protein tau: entropic, structural, and hydrophobic factors.

Globular denatured proteins have structural properties similar to those of random coils. Experiments on denatured proteins have shown that when the temperature is increased thermal compaction may take place, resulting in a reduction of their radius of gyration Rg to range between 5% and 35% of its initial value. This phenomenon has been attributed to various causes, namely entropic, hydrophobic, and structural factors. The intrinsically disordered protein tau, which helps in nucleating and stabilizing microtubules in the axons of the neurons, also undergoes a relevant compaction process: when its temperature is increased from 293 K to 333 K its gyration radius decreases by 18%. We have performed an atomistic simulation of this molecule, at the lowest and highest temperatures of the mentioned interval, using both standard molecular dynamics and metadynamics, in parallel with small-angle X-ray scattering experiments. Using the fit of the experimental data and a genetic algorithm to select the most probable configurations among those produced in both atomistic simulations (standard MD and metadynamics), we were able to compute relevant changes, related to the temperature increase, in the average angles between residues, in the transient secondary structures, in the solvent accessible surface area, and in the number of intramolecular H-bonds. The analysis of the data showed how to decompose the compaction phenomenon into three contributions. An estimate of the entropic contribution to the compaction was obtained using the changes in the mean values of the angles between contiguous residues. The computation of the solvent accessible surface at the two temperatures allowed an estimation of the second factor contributing to the compaction, namely the increase in the hydrophobic interaction. We also measured the change in the average number of residues temporarily being in α-helices, 3-helices, PP II helices, ß-sheets and ß-turns. Those changes in the secondary structure population produce a reduction in the contour length of the protein, yielding a structural contribution to the reduction of Rg. This analysis shows that in tau the entropic factor accounts for about 60% of the compaction, the hydrophobic factor for about 25%, and the change in the secondary structure for about 15%.

Genetically induced dysfunctions of Kir2.1 channels: implications for short QT3 syndrome and autism-epilepsy phenotype.

Short QT3 syndrome (SQT3S) is a cardiac disorder characterized by a high risk of mortality and associated with mutations in Kir2.1 (KCNJ2) channels. The molecular mechanisms leading to channel dysfunction, cardiac rhythm disturbances and neurodevelopmental disorders, potentially associated with SQT3S, remain incompletely understood. Here, we report on monozygotic twins displaying a short QT interval on electrocardiogram recordings and autism-epilepsy phenotype. Genetic screening identified a novel KCNJ2 variant in Kir2.1 that (i) enhanced the channel's surface expression and stability at the plasma membrane, (ii) reduced protein ubiquitylation and degradation, (iii) altered protein compartmentalization in lipid rafts by targeting more channels to cholesterol-poor domains and (iv) reduced interactions with caveolin 2. Importantly, our study reveals novel physiological mechanisms concerning wild-type Kir2.1 channel processing by the cell, such as binding to both caveolin 1 and 2, protein degradation through the ubiquitin-proteasome pathway; in addition, it uncovers a potential multifunctional site that controls Kir2.1 surface expression, protein half-life and partitioning to lipid rafts. The reported mechanisms emerge as crucial also for proper astrocyte function, suggesting the need for a neuropsychiatric evaluation in patients with SQT3S and offering new opportunities for disease management.

Effect of DNA on the conformational dynamics of the endonucleases I-DmoI as provided by molecular dynamics simulations.

The conformational behavior of the wild-type endonucleases I-DmoI and two of its mutants has been studied in the presence and in the absence of DNA target sequences by means of extended molecular dynamics simulations. Our results show that in the absence of DNA, the three protein forms explore a similar essential conformational space, whereas when bound to the same DNA target sequence of 25 base pairs, they diversify and restrain the subspace explored. In addition, the differences in the essential subspaces explored by the residues near the catalytic site for both the bound and unbound forms are discussed in background of the experimental protein activity.

Isoamphipathic antibacterial molecules regulating activity and toxicity through positional isomerism.

Peptidomimetic antimicrobials exhibit a selective interaction with bacterial cells over mammalian cells once they have achieved an optimum amphiphilic balance (hydrophobicity/hydrophilicity) in the molecular architecture. To date, hydrophobicity and cationic charge have been considered the crucial parameters to attain such amphiphilic balance. However, optimization of these properties is not enough to circumvent unwanted toxicity towards mammalian cells. Hence, herein, we report new isoamphipathic antibacterial molecules (IAMs: 1-3) where positional isomerism was introduced as one of the guiding factors for molecular design. This class of molecules displayed good (MIC = 1-8 µg mL-1 or µM) to moderate [MIC = 32-64 µg mL-1 (32.2-64.4 µM)] antibacterial activity against multiple Gram-positive and Gram-negative bacteria. Positional isomerism showed a strong influence on regulating antibacterial activity and toxicity for ortho [IAM-1: MIC = 1-32 µg mL-1 (1-32.2 µM), HC50 = 650 µg mL-1 (654.6 µM)], meta [IAM-2: MIC = 1-16 µg mL-1 (1-16.1 µM), HC50 = 98 µg mL-1 (98.7 µM)] and para [IAM-3: MIC = 1-16 µg mL-1 (1-16.1 µM), HC50 = 160 µg mL-1 (161.1 µM)] isomers. Co-culture studies and investigation of membrane dynamics indicated that ortho isomer, IAM-1 exerted more selective activity towards bacterial over mammalian membranes, compared to meta and para isomers. Furthermore, the mechanism of action of the lead molecule (IAM-1) has been characterized through detailed molecular dynamics simulations. In addition, the lead molecule displayed substantial efficacy against dormant bacteria and mature biofilms, unlike conventional antibiotics. Importantly, IAM-1 exhibited moderate in vivo activity against MRSA wound infection in a murine model with no detectable dermal toxicity. Altogether, the report explored the design and development of isoamphipathic antibacterial molecules to establish the role of positional isomerism in achieving selective and potential antibacterial agents.

Conformational switch of a flexible loop in human laminin receptor determines laminin-1 interaction.

The 37/67-kDa human laminin receptor(LamR) is a cell surface protein that interacts with molecules located in the extra-cellular matrix. In particular, interactions between LamR and laminins play a major role in mediating changes in the cellular environment that affect cell adhesion, neurite outgrowth, tumor growth and metastasis. The exact interaction mode of laminin-1 and LamR is not fully understood. Laminin-1 is thought to bind to LamR through interaction with the so-called peptide G (residues 161­180) and the C-terminal helix (residues 205­229). Here we performed 100-ns atomistic force field based molecular dynamics simulations to explore the structure and dynamics of LamR related to laminin-1 interactions. Our main finding is that loop 188­197 in the C-terminal region is highly flexible. It undergoes a major change resulting in a conformational switch that partially solvent exposes the R180 residue in the final part of the G peptide. So, R180 could contribute to laminin-1 binding. Projection of the simulations along the first two principal components also confirms the importance of this conformational switch in the LamR. This may be a basic prerequisite to clarify the key structural determinants of the interaction of LamR with laminin-1.

Episodic ataxia type 1 mutations affect fast inactivation of K+ channels by a reduction in either subunit surface expression or affinity for inactivation domain.

Episodic ataxia type 1 (EA1) is an autosomal dominant disorder characterized by continuous myokymia and episodic attacks of ataxia. Mutations in the gene KCNA1 that encodes the voltage-gated potassium channel Kv1.1 are responsible for EA1. In several brain areas, Kv1.1 coassembles with Kv1.4, which confers N-type inactivating properties to heteromeric channels. It is therefore likely that the rate of inactivation will be determined by the number of Kv1.4 inactivation particles, as set by the precise subunit stoichiometry. We propose that EA1 mutations affect the rate of N-type inactivation either by reduced subunit surface expression, giving rise to a reduced number of Kv1.1 subunits in heterotetramer Kv1.1-Kv1.4 channels, or by reduced affinity for the Kv1.4 inactivation domain. To test this hypothesis, quantified amounts of mRNA for Kv1.4 or Kv1.1 containing selected EA1 mutations either in the inner vestibule of Kv1.1 on S6 or in the transmembrane regions were injected into Xenopus laevis oocytes and the relative rates of inactivation and stoichiometry were determined. The S6 mutations, V404I and V408A, which had normal surface expression, reduced the rate of inactivation by a decreased affinity for the inactivation domain while the mutations I177N in S1 and E325D in S5, which had reduced subunit surface expression, increased the rate of N-type inactivation due to a stoichiometric increase in the number of Kv1.4 subunits.

Autism with seizures and intellectual disability: possible causative role of gain-of-function of the inwardly-rectifying K+ channel Kir4.1.

The inwardly-rectifying potassium channel Kir4.1 is a major player in the astrocyte-mediated regulation of [K(+)](o) in the brain, which is essential for normal neuronal activity and synaptic functioning. KCNJ10, encoding Kir4.1, has been recently linked to seizure susceptibility in humans and mice, and is a possible candidate gene for Autism Spectrum Disorders (ASD). In this study, we performed a mutational screening of KCNJ10 in 52 patients with epilepsy of "unknown cause" associated with impairment of either cognitive or communicative abilities, or both. Among them, 14 patients fitted the diagnostic criteria for ASD. We identified two heterozygous KCNJ10 mutations (p.R18Q and p.V84M) in three children (two unrelated families) with seizures, ASD, and intellectual disability. The mutations replaced amino acid residues that are highly conserved throughout evolution and were undetected in about 500 healthy chromosomes. The effects of mutations on channel activity were functionally assayed using a heterologous expression system. These studies indicated that the molecular mechanism contributing to the disorder relates to an increase in either surface-expression or conductance of the Kir4.1 channel. Unlike previous syndromic associations of genetic variants in KCNJ10, the pure neuropsychiatric phenotype in our patients suggests that the new mutations affect K(+) homeostasis mainly in the brain, by acting through gain-of-function defects. Dysfunction in astrocytic-dependent K(+) buffering may contribute to autism/epilepsy phenotype, by altering neuronal excitability and synaptic function, and may represent a new target for novel therapeutic approaches.

Phosphorylation, Mg-ADP, and Inhibitors Differentially Shape the Conformational Dynamics of the A-Loop of Aurora-A.

The conformational state of the activation loop (A-loop) is pivotal for the activity of most protein kinases. Hence, the characterization of the conformational dynamics of the A-loop is important to increase our understanding of the molecular processes related to diseases and to support the discovery of small molecule kinase inhibitors. Here, we carry out a combination of molecular dynamics (MD) and essential dynamics (ED) analyses to fully map the effects of phosphorylation, ADP, and conformation disrupting (CD) inhibitors (i.e., CD532 and MLN8054) on the dynamics of the A-loop of Aurora-A. MD revealed that the stability of the A-loop in an open conformation is enhanced by single phospho-Thr-288, while paradoxically, the presence of a second phosphorylation at Thr-287 decreases such stability and renders the A-loop more fluctuant in time and space. Moreover, we found that this post-translational modification has a significant effect on the direction of the A-loop motions. ED analysis suggests that the presence of the phosphate moiety induces the dynamics of Aurora-A to sample two distinct energy minima, instead of a single large minimum, as in unphosphorylated Aurora-A states. This observation indicates that the conformational distributions of Aurora-A with both single and double phospho-threonine modifications are remarkably different from the unphosphorylated state. In the closed states, binding of CD532 and MLN8054 inhibitors has the effect of increasing the distance of the N- and C-lobes of the kinase domain of Aurora-A, and the angle analysis between those two lobes during MD simulations showed that the N- and C-lobes are kept more open in presence of CD532, compared to MLN8054. As the A-loop is a common feature of Aurora protein kinases, our studies provide a general description of the conformational dynamics of this structure upon phosphorylation and different ligands binding.

Diverse driving forces underlie the invariant occurrence of the T42A, E139D, I282V and T468M SHP2 amino acid substitutions causing Noonan and LEOPARD syndromes.

Missense PTPN11 mutations cause Noonan and LEOPARD syndromes (NS and LS), two developmental disorders with pleiomorphic phenotypes. PTPN11 encodes SHP2, an SH2 domain-containing protein tyrosine phosphatase functioning as a signal transducer. Generally, different substitutions of a particular amino acid residue are observed in these diseases, indicating that the crucial factor is the residue being replaced. For a few codons, only one substitution is observed, suggesting the possibility of specific roles for the residue introduced. We analyzed the biochemical behavior and ligand-binding properties of all possible substitutions arising from single-base changes affecting codons 42, 139, 279, 282 and 468 to investigate the mechanisms underlying the invariant occurrence of the T42A, E139D and I282V substitutions in NS and the Y279C and T468M changes in LS. Our data demonstrate that the isoleucine-to-valine change at codon 282 is the only substitution at that position perturbing the stability of SHP2's closed conformation without impairing catalysis, while the threonine-to-alanine change at codon 42, but not other substitutions of that residue, promotes increased phosphopeptide-binding affinity. The recognition specificity of the C-SH2 domain bearing the E139D substitution differed substantially from its wild-type counterpart acquiring binding properties similar to those observed for the N-SH2 domain, revealing a novel mechanism of SHP2's functional dysregulation. Finally, while functional selection does not seem to occur for the substitutions at codons 279 and 468, we point to deamination of the methylated cytosine at nucleotide 1403 as the driving factor leading to the high prevalence of the T468M change in LS.