Intermolecular insights into allosteric inhibition of histone lysine-specific demethylase 1

BackgroundHistone lysine-specific demethylase 1 (LSD1) has become a potential anticancer target for the novel drug discovery. Recent reports have shown that SP2509 and its derivatives strongly inhibit LSD1 as allosteric inhibitors. However, the binding mechanism of these allosteric inhibitors in the allosteric site of LSD1 is not known yet.MethodsThe stability and binding mechanism of allosteric inhibitors in the binding site of LSD1 were evaluated by molecular docking, ligand-based pharmacophore, molecular dynamics (MD) simulations, molecular mechanics generalized born surface area (MM/GBSA) analysis, quantum mechanics/molecular mechanics (QM/MM) calculation and Hirshfeld surface analysis.ResultsThe conformational geometry and the intermolecular interactions of allosteric inhibitors showed high binding affinity towards allosteric site of LSD1 with the neighboring amino acids (Gly358, Cys360, Leu362, Asp375 and Glu379). Meanwhile, MD simulations and MM/GBSA analysis were performed on selected allosteric inhibitors in complex with LSD1 protein, which confirmed the high stability and binding affinity of these inhibitors in the allosteric site of LSD1.ConclusionThe simulation results revealed the crucial factors accounting for allosteric inhibitors of LSD1, including different protein–ligand interactions, the positions and conformations of key residues, and the ligands flexibilities. Meanwhile, a halogen bond interaction between chlorine atom of ligand and key residues Trp531 and His532 was recurrent in our analysis confirming its importance.General significanceOverall, our research analyzed in depth the binding modes of allosteric inhibitors with LSD1 and could provide useful information for the design of novel allosteric inhibitors.     

Homology modeling,docking, and molecular dynamics simulation of the receptor GALR2 and its interactions with galanin and a positive allosteric modulator

Galanin receptor type 2 (GALR2) is a class A G-protein-coupled receptor (GPCR), and it has been reported that orthosteric ligands and positive allosteric modulators (PAMs) of GALR2 could potentially be used to treat epilepsy. So far, the X-ray structure of this receptor has not been resolved, and knowledge of the 3D structure of GALR2 may prove informative in attempts to design novel ligands and to explore the mechanism for the allosteric modulation of this receptor. In this study, homology modeling was used to obtain several GALR2 models using known templates. ProSA-web Z-scores and Ramachandran plots as well as pre-screening against a test dataset of known compounds were all utilized to select the best model of GALR2. Molecular dockings of galanin (a peptide) and a nonpeptide ligand were carried out to choose the (GALR2 model)–galanin complex that showed the closest agreement with the corresponding experimental data. Finally, a 50-ns MD simulation was performed to study the interactions between the GALR2 model and the synthetic and endogenous ligands. The results from docking and MD simulation showed that, besides the reported residues, Tyr160^4.60, Ile105^3.32, Ala274^7.35, and Tyr163^ECL2 also appear to play important roles in the binding of galanin. The potential allosteric binding pockets in the GALR2 model were then investigated via MD simulation. The results indicated that the mechanism for the allosteric modulation caused by PAMs is the binding of the PAM at pocket III, which is formed by galanin, ECL2, TM2, TM3, and ECL1; this results in the disruption of the Na⁺-binding site and/or the Na⁺ ion pathway, leading to GALR2 agonism.     

Negative allosteric modulators of cannabinoid receptor 2: protein modeling,binding site identification and molecular dynamics simulations in the presence of an orthosteric agonist

Abstract

Selective activation of the cannabinoid receptor subtype 2 (CB2) shows promise for treating pain, inflammation, multiple sclerosis, cancer, ischemic/reperfusion injury and osteoporosis. Target selectivity and off-target side effects are two major limiting factors for orthosteric ligands, and therefore, the search for allosteric modulators (AMs) is a widely used drug discovery approach. To date, only a limited number of negative CB2 AMs have been identified, possessing only micromolar activity at best, and the CB2 receptor’s allosteric site(s) are not well characterized. Herein, we used computational approaches including receptor modeling, site mapping, docking, molecular dynamics (MD) simulations and binding free energy calculations to predict, characterize and validate allosteric sites within the complex of the CB2 receptor with bound orthosteric agonist CP55,940. After docking of known negative CB2 allosteric modulators (NAMs), dihydro-gambogic acid (DHGA) and trans-β-caryophyllene (TBC) (note that TBC also shows agonist activity), at the predicted allosteric sites, the best total complex with CB2, CP55,940 and NAM was embedded into a hydrated lipid bilayer and subjected to a 200 ns MD simulation. The presence of an AM affected the CB2–CP55,940 complex, altering the relative positioning of the toggle switch residues and promoting a strong π–π interaction between Phe117^3.36 and Trp258^6.48. Binding of either TBC or DHGA to a putative allosteric pocket directly adjacent to the orthosteric ligand reduced the binding free energy of CP55,940, which is consistent with the expected effect of a negative AM. The identified allosteric sites present immense scope for the discovery of novel classes of CB2 AMs.     

Construction and optimization of a family of genetically encoded metabolite sensors by semirational protein engineering

A family of genetically-encoded metabolite sensors has been constructed using bacterial periplasmic binding proteins (PBPs) linearly fused to protein fluorophores. The ligand-induced conformational change in a PBP allosterically regulates the relative distance and orientation of a fluorescence resonance energy transfer (FRET)-compatible protein pair. Ligand binding is transduced into a macroscopic FRET observable, providing a reagent for in vitro and in vivo ligand-measurement and visualization. Sensors with a higher FRET signal change are required to expand the dynamic range and allow visualization of subtle analyte changes under high noise conditions. Various observations suggest that factors other than inter-fluorophore separation contribute to FRET transfer efficiency and the resulting ligand-dependent spectral changes. Empirical and rational protein engineering leads to enhanced allosteric linkage between ligand binding and chromophore rearrangement; modifications predicted to decrease chromophore rotational averaging enhance the signal change, emphasizing the importance of the rotational freedom parameter kappa2 to FRET efficiency. Tighter allosteric linkage of the PBP and the fluorophores by linker truncation or by insertion of chromophores into the binding protein at rationally designed sites gave rise to sensors with improved signal change. High-response sensors were obtained with fluorescent proteins attached to the same binding PBP lobe, suggesting that indirect allosteric regulation during the hinge-bending motion is sufficient to give rise to a FRET response. The optimization of sensors for glucose and glutamate, ligands of great clinical interest, provides a general framework for the manipulation of ligand-dependent allosteric signal transduction mechanisms.     

A new allosteric site in glycogen phosphorylase b as a target for drug interactions

BACKGROUND: In muscle and liver, glycogen concentrations are regulated by the coordinated activities of glycogen phosphorylase (GP) and glycogen synthase. GP exists in two forms: the dephosphorylated low-activity form GPb and the phosphorylated high-activity form GPa. In both forms, allosteric effectors can promote equilibrium between a less active T state and a more active R state. GP is a possible target for drugs that aim to prevent unwanted glycogen breakdown and to stimulate glycogen synthesis in non-insulin-dependent diabetes. As a result of a data bank search, 5-chloro-1H-indole-2-carboxylic acid (1-(4-fluorobenzyl)-2-(4-hydroxypiperidin-1-yl)-2-oxoethy l)amide, CP320626, was identified as a potent inhibitor of human liver GP. Structural studies have been carried out in order to establish the mechanism of this unusual inhibitor. RESULTS: The structure of the cocrystallised GPb-CP320626 complex has been determined to 2.3 A resolution. CP320626 binds at a site located at the subunit interface in the region of the central cavity of the dimeric structure. The site has not previously been observed to bind ligands and is some 15 A from the AMP allosteric site and 33 A from the catalytic site. The contacts between GPb and CP320626 comprise six hydrogen bonds and extensive van der Waals interactions that create a tight binding site in the T-state conformation of GPb. In the R-state conformation of GPa these interactions are significantly diminished. CONCLUSIONS: CP320626 inhibits GPb by binding at a new allosteric site. Although over 30 A from the catalytic site, the inhibitor exerts its effects by stabilising the T state at the expense of the R state and thereby shifting the allosteric equilibrium between the two states. The new allosteric binding site offers a further recognition site in the search for improved GP inhibitors.     

Shared functional attributes between the mecA gene product of Staphylococcus sciuri and penicillin-binding protein 2a of methicillin-resistant Staphylococcus aureus

The genome of Staphylococcus aureus is constantly in a state of flux, acquiring genes that enable the bacterium to maintain resistance in the face of antibiotic pressure. The acquisition of the mecA gene from an unknown origin imparted S. aureus with broad resistance to beta-lactam antibiotics, with the resultant strain designated as methicillin-resistant S. aureus (MRSA). Epidemiological and genetic evidence suggests that the gene encoding PBP 2a of MRSA might have originated from Staphylococcus sciuri, an animal pathogen, where it exists as a silent gene of unknown function. We synthesized, cloned, and expressed the mecA gene of S. sciuri in Escherichia coli, and the protein product was purified to homogeneity. Biochemical characterization and comparison of the protein to PBP 2a of S. aureus revealed them to be highly similar. These characteristics start with sequence similarity but extend to biochemical behavior in inhibition by beta-lactam antibiotics, to the existence of an allosteric site for binding of bacterial peptidoglycan, to the issues of the sheltered active site, and to the need for conformational change in making the active site accessible to the substrate and the inhibitors. Altogether, the evidence strongly argues that the kinship between the two proteins is deep-rooted on the basis of many biochemical attributes quantified in this study.     

Time-resolved fluorescence studies of heterotropic ligand binding to cytochrome P450 3A4

Cytochrome P450 3A4 (CYP3A4) is a major enzymatic determinant of drug and xenobiotic metabolism that demonstrates remarkable substrate diversity and complex kinetic properties. The complex kinetics may result, in some cases, from multiple binding of ligands within the large active site or from an effector molecule acting at a distal allosteric site. Here, the fluorescent probe TNS (2-p-toluidinylnaphthalene-6-sulfonic acid) was characterized as an active site fluorescent ligand. UV-vis difference spectroscopy revealed a TNS-induced low-spin heme absorbance spectrum with an apparent K(d) of 25.4 +/- 2 microM. Catalytic turnover using 7-benzyloxyquinoline (7-BQ) as a substrate demonstrated TNS-dependent inhibition with an IC(50) of 9.9 +/- 0.1 microM. These results suggest that TNS binds in the CYP3A4 active site. The steady-state fluorescence of TNS increased upon binding to CYP3A4, and fluorescence titrations yielded a K(d) of 22.8 +/- 1 microM. Time-resolved frequency-domain measurement of TNS fluorescence lifetimes indicates a testosterone (TST)-dependent decrease in the excited-state lifetime of TNS, concomitant with a decrease in the steady-state fluorescence intensity. In contrast, the substrate erythromycin (ERY) had no effect on TNS lifetime, while it decreased the steady-state fluorescence intensity. Together, the results suggest that TNS binds in the active site of CYP3A4, while the first equivalent of TST binds at a distant allosteric effector site. Furthermore, the results are the first to indicate that TST bound to the effector site can modulate the environment of the heterotropic ligand.     

How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics

The tunnel region at triosephosphate isomerase (TIM)’s dimer interface, distant from its catalytic site, is a target site for certain benzothiazole derivatives that inhibit TIM’s catalytic activity in Trypanosoma cruzi, the parasite that causes Chagas disease. We performed multiple 100-ns molecular-dynamics (MD) simulations and elastic network modeling (ENM) on both apo and complex structures to shed light on the still unclear inhibitory mechanism of one such inhibitor, named bt10. Within the time frame of our MD simulations, we observed stabilization of aromatic clusters at the dimer interface and enhancement of intersubunit hydrogen bonds in the presence of bt10, which point to an allosteric effect rather than destabilization of the dimeric structure. The collective dynamics dictated by the topology of TIM is known to facilitate the closure of its catalytic loop over the active site that is critical for substrate entrance and product release. We incorporated the ligand’s effect on vibrational dynamics by applying mixed coarse-grained ENM to each one of 54,000 MD snapshots. Using this computationally efficient technique, we observed altered collective modes and positive shifts in eigenvalues due to the constraining effect of bt10 binding. Accordingly, we observed allosteric changes in the catalytic loop’s dynamics, flexibility, and correlations, as well as the solvent exposure of catalytic residues. A newly (to our knowledge) introduced technique that performs residue-based ENM scanning of TIM revealed the tunnel region as a key binding site that can alter global dynamics of the enzyme.     

Ligand modifications to reduce the relative resistance of multi-drug resistant HIV-1 protease

Proper proteolytic processing of the HIV-1 Gag/Pol polyprotein is required for HIV infection and viral replication. This feature has made HIV-1 protease an attractive target for antiretroviral drug design for the treatment of HIV-1 infected patients. To examine the role of the P1 and P1′positions of the substrate in inhibitory efficacy of multi-drug resistant HIV-1 protease 769 (MDR 769), we performed a series of structure–function studies. Using the original CA/p2 cleavage site sequence, we generated heptapeptides containing one reduced peptide bond with an L to F and A to F double mutation at P1 and P1′ (F-r-F), and an A to F at P1′ (L-r-F) resulting in P1/P1′ modified ligands. Here, we present an analysis of co-crystal structures of CA/p2 F-r-F, and CA/p2 L-r-F in complex with MDR 769. To examine conformational changes in the complex structure, molecular dynamic (MD) simulations were performed with MDR769–ligand complexes. MD trajectories show the isobutyl group of both the lopinavir analog and the CA/p2 L-r-F substrate cause a conformational change of in the active site of MDR 769. IC₅₀ measurements suggest the non identical P1/P1′ ligands (CA/p2 L-r-F and lopinavir analog) are more effective against MDR proteases as opposed to identical P1/P1′ligands. Our results suggest that a non identical P1/P1′composition may be more favorable for the inhibition of MDR 769 as they induce conformational changes in the active site of the enzyme resulting in disruption of the two-fold symmetry of the protease, thus, stabilizing the inhibitor in the active site.     

Overexpression and enzymatic characterization of Neisseria gonorrhoeae penicillin-binding protein 4.

The penicillin-binding proteins (PBPs) are ubiquitous bacterial enzymes involved in cell wall biosynthesis, and are the targets of the beta-lactam antibiotics. The low molecular mass Neisseria gonorrhoeae PBP 4 (NG PBP 4) is the fourth PBP revealed in the gonococcal genome. NG PBP 4 was cloned, overexpressed, purified, and characterized for beta-lactam binding, DD-carboxypeptidase activity, acyl-donor substrate specificity, transpeptidase activity, inhibition by a number of active site directed reagents, and pH profile. NG PBP 4 was efficiently acylated by penicillin (30,000 m-1.s-1). Against a set of five alpha- and epsilon-substituted l-Lys-D-Ala-D-Ala substrates, NG PBP 4 exhibited wide variation in specificity with a preference for N epsilon-acylated substrates, suggesting a possible preference for crosslinked pentapeptide substrates in the cell wall. Substrates with an N epsilon-Cbz group demonstrated pronounced substrate inhibition. NG PBP 4 showed 30-fold higher activity against the depsipeptide Lac-ester substrate than against the analogous peptide substrate, an indication that k2 (acylation) is rate determining for carboxypeptidase activity. No transpeptidase activity was apparent in a model transpeptidase reaction. Among a number of active site-directed agents, N-chlorosuccinimide, elastinal, iodoacetamide, iodoacetic acid, and phenylglyoxal gave substantial inhibition, and methyl boronic acid gave modest inhibition. The pH profile for activity against Ac2-l-Lys-D-Ala-d-Ala (kcat/Km) was bell-shaped, with pKa values at 6.9 and 10.1. Comparison of the enzymatic properties of NG PBP 4 with other DD-carboxypeptidases highlights both similarities and differences within these enzymes, and suggests the possibility of common mechanistic roles for the two highly conserved active site lysines in Class A and C low molecular mass PBPs.     

Isolation of an individual allosteric interaction in tetrameric phosphofructokinase from Bacillus stearothermophilus

Phosphofructokinase from Bacillus stearothermophilus (BsPFK) is a model allosteric enzyme system in which the interactions between substrates and allosteric effectors have been extensively studied. However, the oligomeric nature of BsPFK has made it difficult to determine the molecular basis of the allosteric regulation because of the multitude of different types of heterotropic and homotropic interactions that are possible between the four active sites and four allosteric sites in the native tetramer. In an attempt to alleviate the complexity of the system and thereby allow the quantitation of a single interaction between one active site and one allosteric site, site-directed mutagenesis has been coupled with a hybrid-forming scheme to create and isolate a tetramer of BsPFK in which only a single active site and a single allosteric site are capable of binding their respective ligands with high (i.e., near wild type) affinity. Characterization of this single allosteric interaction indicates that the free energy involved in the inhibition by the allosteric effector phosphoenolpyruvate (PEP) is 1.48 +/- 0.15 kcal/mol compared to the 3.58 +/- 0.02 kcal/mol measured for the enzyme.     

Molecular Mechanisms of Bitopic Ligand Engagement with the M1 Muscarinic Acetylcholine Receptor

TBPB and 77-LH-28-1 are selective agonists of the M₁ muscarinic acetylcholine receptor (mAChR) that may gain their selectivity through a bitopic mechanism, interacting concomitantly with the orthosteric site and part of an allosteric site. The current study combined site-directed mutagenesis, analytical pharmacology,and molecular modeling to gain further insights into the structural basis underlying binding and signaling by these agonists. Mutations within the orthosteric binding site caused similar reductions in affinity and signaling efficacy for both selective and prototypical orthosteric ligands. In contrast, the mutation of residues within transmembrane helix (TM) 2 and the second extracellular loop (ECL2) discriminated between the different classes of ligand. In particular, ECL2 appears to be involved in the selective binding of bitopic ligands and in coordinating biased agonism between intracellular calcium mobilization and ERK1/2 phosphorylation. Molecular modeling of the interaction between TBPB and the M₁ mAChR revealed a binding pose predicted to extend from the orthosteric site up toward a putative allosteric site bordered by TM2, TM3, and TM7, thus consistent with a bitopic mode of binding. Overall, these findings provide valuable structural and mechanistic insights into bitopic ligand actions and receptor activation and support a role for ECL2 in dictating the active states that can be adopted by a G protein-coupled receptor. This may enable greater selective ligand design and development for mAChRs and facilitate improved identification of bitopic ligands.     

Distances between structural metal ion, substrates, and allosteric modifier of fructose bisphosphatase

The binding of two paramagnetic probes within a subunit of fructose bisphosphatase, viz., Mn2+ at a structural site and a nitroxide spin-label at a sulfhydryl site, has permitted the measurement of NMR and electron paramagnetic resonance (EPR) relaxation rates to map the active and allosteric site topography. Distances from these loci to the phosphoryl of fructose 6-phosphate (Fru-6-P) and inorganic phosphate (Pi) and four nuclei of adenosine 5'-phosphate (AMP) (the phosphorus nucleus, H-8, H-2, and H-1') were obtained. These measurements located the Mn2+ approximately equidistant from the two phosphoryl moieties of the product ligands Fru-6-P and Pi and in close proximity to the AMP. The adenosine moiety of the latter is oriented anti. Analysis of EPR data revealed that the nitroxide group is approximately 16 A from the structural Mn2+ site. Calculation of the residence times for the hydrolysis reaction products suggests that their dissociation should not be rate limiting in the overall reaction cycle.     

Quercetin protected against isoniazide‐induced HepG2 cell apoptosis by activating the SIRT1/ERK pathway

Isoniazid (INH) is one of the most commonly used antituberculosis drugs, but its clinical applications have been limited by severe hepatic toxicity. Quercetin (Que), a natural flavonoid, has been proved to have many medicinal properties. This study aimed to clarify the possible protective effects of Que against INH‐induced hepatotoxicity using HepG2 cells. Our results indicated that Que significantly increased cell viability, superoxide dismutase, and GSH levels, while decreased alanine aminotransferase/aspartate aminotransferase levels. Besides, Que significantly abrogated INH‐induced cell apoptosis by upregulating the expression levels of Bcl‐2 and decreasing the levels of Bax, cleaved caspase‐3, and cleaved caspase‐9. Furthermore, Que obviously reversed the inhibition of INH on Sirtuin 1 (SIRT1) expression and extracellular signal‐regulated kinase (ERK) phosphorylation. Next, the SIRT1 inhibitor EX527 blocked the enhancement of Que upon ERK phosphorylation. Notably, EX527 partially abolished the beneficial effects of Que. In brief, our results provided the first evidence that Que protected against INH‐induced HepG2 cells by regulating the SIRT1/ERK pathway.     

On the role of the conformational flexibility of the active-site lid on the allosteric kinetics of glucosamine-6-phosphate deaminase

The active site of glucosamine-6-phosphate deaminase from Escherichia coli (GlcN6P deaminase, EC 3.5.99.6) has a complex lid formed by two antiparallel beta-strands connected by a helix-loop segment (158-187). This motif contains Arg172, which is a residue involved in binding the substrate in the active-site, and three residues that are part of the allosteric site, Arg158, Lys160 and Thr161. This dual binding role of the motif forming the lid suggests that it plays a key role in the functional coupling between active and allosteric sites. Previous crystallographic work showed that the temperature coefficients of the active-site lid are very large when the enzyme is in its T allosteric state. These coefficients decrease in the R state, thus suggesting that this motif changes its conformational flexibility as a consequence of the allosteric transition. In order to explore the possible connection between the conformational flexibility of the lid and the function of the deaminase, we constructed the site-directed mutant Phe174-Ala. Phe174 is located at the C-end of the lid helix and its side-chain establishes hydrophobic interactions with the remainder of the enzyme. The crystallographic structure of the T state of Phe174-Ala deaminase, determined at 2.02 A resolution, shows no density for the segment 162-181, which is part of the active-site lid (PDB 1JT9). This mutant form of the enzyme is essentially inactive in the absence of the allosteric activator, N-acetylglucosamine-6-P although it recovers its activity up to the wild-type level in the presence of this ligand. Spectrometric and binding studies show that inactivity is due to the inability of the active-site to bind ligands when the allosteric site is empty. These data indicate that the conformational flexibility of the active-site lid critically alters the binding properties of the active site, and that the occupation of the allosteric site restores the lid conformational flexibility to a functional state.     

Molecular dynamics simulation of the complex PBP-2x with drug cefuroxime to explore the drug resistance mechanism of Streptococcus suis R61

Drug resistance of Streptococcus suis strains is a worldwide problem for both humans and pigs. Previous studies have noted that penicillin-binding protein (PBPs) mutation is one important cause of β-lactam antibiotic resistance. In this study, we used the molecular dynamics (MD) method to study the interaction differences between cefuroxime (CES) and PBP2x within two newly sequenced Streptococcus suis: drug-sensitive strain A7, and drug-resistant strain R61. The MM-PBSA results proved that the drug bound much more tightly to PBP2x in A7 (PBP2x-A7) than to PBP2x in R61 (PBP2x-R61). This is consistent with the evidently different resistances of the two strains to cefuroxime. Hydrogen bond analysis indicated that PBP2x-A7 preferred to bind to cefuroxime rather than to PBP2x-R61. Three stable hydrogen bonds were formed by the drug and PBP2x-A7, while only one unstable bond existed between the drug and PBP2x-R61. Further, we found that the Gln569, Tyr594, and Gly596 residues were the key mutant residues contributing directly to the different binding by pair wise energy decomposition comparison. By investigating the binding mode of the drug, we found that mutant residues Ala320, Gln553, and Thr595 indirectly affected the final phenomenon by topological conformation alteration. Above all, our results revealed some details about the specific interaction between the two PBP2x proteins and the drug cefuroxime. To some degree, this explained the drug resistance mechanism of Streptococcus suis and as a result could be helpful for further drug design or improvement.     

Dynamic properties of the N-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two 3-dimensional domain-swapped dimers with different quaternary structures. One dimer is characterized by the swapping of the C-terminal region (C-Dimer) and presents a rather loose structure. The other dimer (N-Dimer) exhibits a very compact structure with exchange of the N-terminal helix. Here we report the results of a molecular dynamics/essential dynamics (MD/ED) study carried out on the N-Dimer. This investigation, which represents the first MD/ED analysis on a three-dimensional domain-swapped enzyme, provides information on the dynamic properties of the active site residues as well as on the global motions of the dimer subunits. In particular, the analysis of the flexibility of the active site residues agrees well with recent crystallographic and site-directed mutagenesis studies on monomeric RNase A, thus indicating that domain swapping does not affect the dynamics of the active sites. A slight but significant rearrangement of N-Dimer quaternary structure, favored by the formation of additional hydrogen bonds at subunit interface, has been observed during the MD simulation. The analysis of collective movements reveals that each subunit of the dimer retains the functional breathing motion observed for RNase A. Interestingly, the breathing motion of the two subunits is dynamically coupled, as they open and close in phase. These correlated motions indicate the presence of active site intercommunications in this dimer. On these bases, we propose a speculative mechanism that may explain negative cooperativity in systems preserving structural symmetry during the allosteric transitions.     

Community analysis of large-scale molecular dynamics simulations elucidated dynamics-driven allostery in tyrosine kinase 2

TYK2 is a nonreceptor tyrosine kinase, member of the Janus kinases (JAK), with a central role in several diseases, including cancer. The JAKs' catalytic domains (KD) are highly conserved, yet the isolated TYK2-KD exhibits unique specificities. In a previous work, using molecular dynamics (MD) simulations of a catalytically impaired TYK2-KD variant (P1104A) we found that this amino acid change of its JAK-characteristic insert (αFG), acts at the dynamics level. Given that structural dynamics is key to the allosteric activation of protein kinases, in this study we applied a long-scale MD simulation and investigated an active TYK2-KD form in the presence of adenosine 5′-triphosphate and one magnesium ion that represents a dynamic and crucial step of the catalytic cycle, in other protein kinases. Community analysis of the MD trajectory shed light, for the first time, on the dynamic profile and dynamics-driven allosteric communications within the TYK2-KD during activation and revealed that αFG and amino acids P1104, P1105, and I1112 in particular, hold a pivotal role and act synergistically with a dynamically coupled communication network of amino acids serving intra-KD signaling for allosteric regulation of TYK2 activity. Corroborating our findings, most of the identified amino acids are associated with cancer-related missense/splice-site mutations of the Tyk2 gene. We propose that the conformational dynamics at this step of the catalytic cycle, coordinated by αFG, underlie TYK2-unique substrate recognition and account for its distinct specificity. In total, this work adds to knowledge towards an in-depth understanding of TYK2 activation and may be valuable towards a rational design of allosteric TYK2-specific inhibitors.