Two pathways mediate interdomain allosteric regulation in pin1

Elucidating the cleavage patterns of allergenic epitopes of Parvalbumin by L. helveticus fermentation through cell envelope proteinase hydrolysis and computer simulation

Parvalbumin (PV) is recognized as a primary allergen in fish. Previous studies showed that L.helveticus Lh191404 reduces the immunogenicity of Atlantic cod through fermentation, primarily by producing cell envelope proteinase (CEP) to hydrolyze the substrate. However, the cleavage sites of allergenic epitopes during fermentation remained unclear. Peptidomic analysis revealed that the enzymatic hydrolysis process can serve as a proxy for the fermentation process to some extent. Both the fermentation and enzymatic hydrolysis groups contained cleavage sites located within the epitopes. Computer simulations demonstrated that the flexible region of PV is the optimal binding site for CEP. Ala4, Asp52, Gln69, Arg76, and Ala77 located within this flexible region were the cleavage sites in a simulated fermentation environment. The majority of the predicted cleavage sites aligned closely with those observed experimentally. MD simulation offers a framework for the targeted screening of LABs that can compromise the structural integrity of target proteins.

Synergistic insights into positive allosteric modulator and agonist using Gaussian accelerated and tau random acceleration simulations in the metabotropic glutamate receptor 2

Schizophrenia is a severe brain disorder that usually produces a lifetime of disability. Related research shows activating metabotropic glutamate receptors holds therapeutic potential. Agonist-positive allosteric modulations (ago-PAMs) not only activate metabotropic glutamate receptors but also enhance glutamate-induced responses, offering a promising treatment strategy. However, the molecular mechanisms by which ago-PAM enhances glutamate-induced responses remain unclear, as does the potential influence of glutamate on ago-PAM. In this study, Gaussian accelerated molecular dynamics and tau random acceleration molecular dynamics simulations were employed to investigate the molecular mechanism between ago-PAM and glutamate in full-length mGlu2. Results suggest that the ago-PAM JNJ-46281222 enhances the binding affinity and residence time of glutamates in the Venus flytrap (VFT) domains by initiating a variant reverse communication from the heptahelical transmembrane (7TM) domains to VFTs via the cysteine-rich domains. Meanwhile, glutamate facilitates the interaction between Trp676 and Glu701 to further induce the relaxation of TM5, promoting the opening of the PAM-binding pocket. Glutamate can also promote the upward rotation of the cyclopropylmethyl group of the JNJ-46281222 to bring the TM6-TM6 distance closer. Nevertheless, it remains uncertain how the binding between mGlu2 and G protein differs when induced by small molecules binding in allosteric sites, orthosteric sites, or both. In conclusion, this study shed new light on the positive coordination relationship between ago-PAM and glutamate in the full-length mGlu2 receptor, which could help develop novel and more effective ago-PAM to treat schizophrenia.

AlloViz: A tool for the calculation and visualisation of protein allosteric communication networks

Allostery, the presence of functional interactions between distant parts of proteins, is a critical concept in the field of biochemistry and molecular biology, particularly in the context of protein function and regulation. Understanding the principles of allosteric regulation is essential for advancing our knowledge of biology and developing new therapeutic strategies. This paper presents AlloViz, an open-source Python package designed to quantitatively determine, analyse, and visually represent allosteric communication networks on the basis of molecular dynamics (MD) simulation data. The software integrates well-known techniques for understanding allosteric properties simplifying the process of accessing, rationalising, and representing protein allostery and communication routes. It overcomes the inefficiency of having multiple methods with heterogeneous implementations and showcases the advantages of using MD simulations and multiple replicas to obtain statistically sound information on protein dynamics; it also enables the calculation of "consensus-like" scores aggregating methods that consider multiple structural aspects of allosteric networks. We demonstrate the features of AlloViz on two proteins: β-arrestin 1, a key player for regulating G protein-coupled receptor (GPCR) signalling, and the protein tyrosine phosphatase 1B, an important pharmaceutical target for allosteric inhibitors. The software includes comprehensive documentation and examples, tutorials, and a user-friendly graphical interface.

A novel bivalent interaction mode underlies a non-catalytic mechanism for Pin1-mediated protein kinase C regulation

Regulated hydrolysis of the phosphoinositide phosphatidylinositol(4,5)-bis-phosphate to diacylglycerol and inositol-1,4,5-P3 defines a major eukaryotic pathway for translation of extracellular cues to intracellular signaling circuits. Members of the lipid-activated protein kinase C isoenzyme family (PKCs) play central roles in this signaling circuit. One of the regulatory mechanisms employed to downregulate stimulated PKC activity is via a proteasome-dependent degradation pathway that is potentiated by peptidyl-prolyl isomerase Pin1. Here, we show that contrary to prevailing models, Pin1 does not regulate conventional PKC isoforms α and βII via a canonical cis-trans isomerization of the peptidyl-prolyl bond. Rather, Pin1 acts as a PKC binding partner that controls PKC activity via sequestration of the C-terminal tail of the kinase. The high-resolution structure of full-length Pin1 complexed to the C-terminal tail of PKCβII reveals that a novel bivalent interaction mode underlies the non-catalytic mode of Pin1 action. Specifically, Pin1 adopts a conformation in which it uses the WW and PPIase domains to engage two conserved phosphorylated PKC motifs, the turn motif and hydrophobic motif, respectively. Hydrophobic motif is a non-canonical Pin1-interacting element. The structural information combined with the results of extensive binding studies and experiments in cultured cells suggest that non-catalytic mechanisms represent unappreciated modes of Pin1-mediated regulation of AGC kinases and other key enzymes/substrates.

An allosteric mechanism for potent inhibition of SARS-CoV-2 main proteinase

The main proteinase (Mpro) of SARS-CoV-2 plays a critical role in cleaving viral polyproteins into functional proteins required for viral replication and assembly, making it a prime drug target for COVID-19. It is well known that noncompetitive inhibition offers potential therapeutic options for treating COVID-19, which can effectively reduce the likelihood of cross-reactivity with other proteins and increase the selectivity of the drug. Therefore, the discovery of allosteric sites of Mpro has both scientific and practical significance. In this study, we explored the binding characteristics and inhibiting process of Mpro activity by two recently reported allosteric inhibitors, pelitinib and AT7519 which were obtained by the X-ray screening experiments, to probe the allosteric mechanism via molecular dynamic (MD) simulations. We found that pelitinib and AT7519 can stably bind to Mpro far from the active site. The binding affinity is estimated to be -24.37 ± 4.14 and - 26.96 ± 4.05 kcal/mol for pelitinib and AT7519, respectively, which is considerably stable compared with orthosteric drugs. Furthermore, the strong binding caused clear changes in the catalytic site of Mpro, thus decreasing the substrate accessibility. The community network analysis also validated that pelitinib and AT7519 strengthened intra- and inter-domain communication of Mpro dimer, resulting in a rigid Mpro, which could negatively impact substrate binding. In summary, our findings provide the detailed working mechanism for the two experimentally observed allosteric sites of Mpro. These allosteric sites greatly enhance the 'druggability' of Mpro and represent attractive targets for the development of new Mpro inhibitors.

Recent advances of Pin1 inhibitors as potential anticancer agents

Pin1 (proline isomerase peptidyl-prolyl isomerase NIMA-interacting-1), as a member of PPIase family, catalyzes cis-trans isomerization of pThr/Ser-Pro amide bonds of its substrate proteins, further regulating cell proliferation, division, apoptosis, and transformation. Pin1 is overexpressed in various cancers and is positively correlated with tumor initiation and progression. Pin1 inhibition can effectively reduce tumor growth and cancer stem cell expansion, block metastatic spread, and restore chemosensitivity, suggesting that targeting Pin1 may be an effective strategy for cancer treatment. Considering the promising therapeutic effects of Pin1 inhibitors on cancers, we herein are intended to comprehensively summarize the reported Pin1 inhibitors, mainly highlighting their structures, biological functions and binding modes, in hope of providing a reference for the future drug discovery.

Dissecting the Allosteric Fine-Tuning of Enzyme Catalysis

Fully understanding the mechanism of allosteric regulation in biomolecules requires separating and examining all of the involved factors. In enzyme catalysis, allosteric effector binding shifts the structure and dynamics of the active site, leading to modified energetic (e.g., energy barrier) and dynamical (e.g., diffusion coefficient) factors underlying the catalyzed reaction rate. Such modifications can be subtle and dependent on the type of allosteric effector, representing a fine-tuning of protein function. The microscopic description of allosteric regulation at the level of function-dictating factors has prospective applications in fundamental and pharmaceutical sciences, which is, however, largely missing so far. Here, we characterize the allosteric fine-tuning of enzyme catalysis, using human Pin1 as an example, by performing more than half-millisecond all-atom molecular dynamics simulations. Changes of reaction kinetics and the dictating factors, including the free energy surface along the reaction coordinate and the diffusion coefficient of the reaction dynamics, under various enzyme and allosteric effector binding conditions are examined. Our results suggest equal importance of the energetic and dynamical factors, both of which can be modulated allosterically, and the combined effect determines the final allosteric output. We also reveal the potential dynamic basis for allosteric modulation using an advanced statistical technique to detect function-related conformational dynamics. Methods developed in this work can be applied to other allosteric systems.

Allosteric inhibition of myosin by phenamacril: a synergistic mechanism revealed by computational and experimental approaches

BACKGROUND

Myosin plays a crucial role in cellular processes, while its dysfunction can lead to organismal malfunction. Phenamacril (PHA), a highly species-specific and non-competitive inhibitor of myosin I (FgMyoI) from Fusarium graminearum, has been identified as an effective fungicide for controlling plant diseases caused by partial Fusarium pathogens, such as wheat scab and rice bakanae. However, the molecular basis of its action is still unclear.

RESULTS

This study used multiple computational approaches first to elucidate the allosteric inhibition mechanism of FgMyoI by PHA at the atomistic level. The results indicated the increase of adenosine triphosphate (ATP) binding affinity upon PHA binding, which might impede the release of hydrolysis products. Furthermore, simulations revealed a broadened outer cleft and a significantly more flexible interface for actin binding, accompanied by a decrease in signaling transduction from the catalytic center to the actin-binding interface. These various effects might work together to disrupt the actomyosin cycle and hinder the ability of motor to generate force. Our experimental results further confirmed that PHA reduces the enzymatic activity of myosin and its binding with actin.

CONCLUSION

Therefore, our findings demonstrated that PHA might suppress the function of myosin through a synergistic mechanism, providing new insights into myosin allostery and offering new avenues for drug/fungicide discovery targeting myosin. © 2023 Society of Chemical Industry.

Exploring the molecular mechanism of cold-adaption of an alkaline protease mutant by molecular dynamics simulations and residue interaction network

Psychrophilic proteases have attracted enormous attention in past decades, due to their high catalytic activity at low temperatures in a wide range of industrial processes, especially in the detergent and leather industries. Among them, H5 is an alkaline protease mutant, which featuring psychrophilic-like behavior, but the reasons that H5 with higher activity at low temperatures are still poorly understood. Herein, the molecular dynamics (MD) simulations combined with residue interaction network (RIN) were utilized to investigate the mechanisms of the cold-adaption of mutant H5. The results demonstrated that two loops involved in the substrate binding G100-S104 and S125-S129 in H5 had higher mobility, and the distance enlargement between the two loops modulated the substrate's accessibility compared with wild type counterpart. Besides, H5 enhanced conformational flexibility by weakening salt bridges and increasing interaction with the solvent. In particular, the absence of Lys251-Asp197-Arg247 salt bridge network may contribute to the structural mobility. Based on the free energy landscape and molecular mechanics Poisson-Boltzmann surface area of the wild type and H5, it was elucidated that H5 possessed a large population of interconvertible conformations, resulting in the weaker substrate binding free energy. The calculated RIN topology parameters such as the average degree, average cluster coefficient, and average path length further verified that the mutant H5 attenuated residue-to-residue interactions. The investigation of the mechanisms by which how the residue mutation affects the stability and activity of enzymes provides a theoretical basis for the development of cold-adapted protease.

A novel bivalent interaction mode underlies a non-catalytic mechanism for Pin1-mediated Protein Kinase C regulation

Regulated hydrolysis of the phosphoinositide phosphatidylinositol(4,5)-bis-phosphate to diacylglycerol and inositol-1,4,5-P3 defines a major eukaryotic pathway for translation of extracellular cues to intracellular signaling circuits. Members of the lipid-activated protein kinase C isoenzyme family (PKCs) play central roles in this signaling circuit. One of the regulatory mechanisms employed to downregulate stimulated PKC activity is via a proteasome-dependent degradation pathway that is potentiated by peptidyl-prolyl isomerase Pin1. Here, we show that contrary to prevailing models, Pin1 does not regulate conventional PKC isoforms α and βII via a canonical cis-trans isomerization of the peptidyl-prolyl bond. Rather, Pin1 acts as a PKC binding partner that controls PKC activity via sequestration of the C-terminal tail of the kinase. The high-resolution structure of Pin1 complexed to the C-terminal tail of PKCβII reveals that a novel bivalent interaction mode underlies the non-catalytic mode of Pin1 action. Specifically, Pin1 adopts a compact conformation in which it engages two conserved phosphorylated PKC motifs, the turn motif and hydrophobic motif, the latter being a non-canonical Pin1-interacting element. The structural information, combined with the results of extensive binding studies and in vivo experiments suggest that non-catalytic mechanisms represent unappreciated modes of Pin1-mediated regulation of AGC kinases and other key enzymes/substrates.

Dynamical characterization and multiple unbinding paths of two PreQ1 ligands in one pocket

Riboswitches naturally regulate gene expression in bacteria by binding to specific small molecules. Class 1 preQ1 riboswitch aptamer is an important model not only for RNA folding but also as a target for designing small molecule antibiotics due to its well-known minimal aptamer domain. Here, we ran a total of 62.4 μs conventional and enhanced-sampling molecular dynamics (MD) simulations to characterize the determinants underlying the binding of the preQ1-II riboswitch aptamer to two preQ1 ligands in one binding pocket. Decomposition of binding free energy suggested that preQ1 ligands at α and β sites interact with four nucleotides (G5, C17, C18, and A30) and two nucleotides (A12 and C31), respectively. Mg2+ ions play a crucial role in both stabilizing the binding pocket and facilitating ligand binding. The flexible preQ1 ligand at the β site leads to the top of the binding pocket loosening and thus pre-organizes the riboswitch for ligand entry. Enhanced sampling simulations further revealed that the preQ1 ligand at the α site unbinds through two orthogonal pathways, which are dependent on whether or not a β site preQ1 ligand is present. One of the two preQ1 ligands has been identified in the binding pocket, which will aid to identify the second preQ1 Ligand. Our work provides new information for designing robust ligands.

Regulation of eukaryotic protein kinases by Pin1, a peptidyl-prolyl isomerase

The peptidyl-prolyl isomerase Pin1 cooperates with proline-directed kinases and phosphatases to regulate multiple oncogenic pathways. Pin1 specifically recognizes phosphorylated Ser/Thr-Pro motifs in proteins and catalyzes their cis-trans isomerization. The Pin1-catalyzed conformational changes determine the stability, activity, and subcellular localization of numerous protein substrates. We conducted a survey of eukaryotic protein kinases that are regulated by Pin1 and whose Pin1 binding sites have been identified. Our analyses reveal that Pin1 target sites in kinases do not fall exclusively within the intrinsically disordered regions of these enzymes. Rather, they fall into three groups based on their location: (i) within the catalytic kinase domain, (ii) in the C-terminal kinase region, and (iii) in regulatory domains. Some of the kinases downregulated by Pin1 activity are tumor-suppressing, and all kinases upregulated by Pin1 activity are functionally pro-oncogenic. These findings further reinforce the rationale for developing Pin1-specific inhibitors as attractive pharmaceuticals for cancer therapy.

Dynamics and Allosteric Information Pathways of Unphosphorylated c-Cbl

Human c-Cbl is a RING-type ligase and plays a central role in the protein degradation cascade. To elucidate its conformational changes related to substrate binding, we performed molecular dynamics simulations of different variants/states of c-Cbl for a cumulative time of 68 μs. Our simulations demonstrate that before the substrate binds, the RING domain samples a broad set of conformational states at a biologically relevant salt concentration, including the closed, partially open, and fully open states, whereas substrate binding leads to a restricted conformational sampling. Phe378 and the C-terminal region play an essential role in stabilizing the partially open state. To visualize the allosteric signal transmission pathways from the substrate-binding site to the 40 Å apart RING domain and identify the critical residues for allostery, we have created a subgraph from the optimal and suboptimal paths. Redundant paths are seen in the SH2 domain where the substrate binds, while the major bottlenecks are found at the junction between the SH2 domain and the linker helix region as well as that between the SH2 domain and the 4H bundle. These bottlenecks separate the paths into two overall routes. The nodes/residues at the bottlenecks on the subgraph are considered allosteric hot spots. This subgraph approach provides a general tool for network visualization and determination of critical residues for allostery. The structurally and allosterically critical residues identified in our work are testable and would provide valuable insights into the emerging strategies for drug discovery, such as targeted protein degradation.

Ligand-specific conformational change drives interdomain allostery in Pin1

Pin1 is a two-domain cell regulator that isomerizes peptidyl-prolines. The catalytic domain (PPIase) and the other ligand-binding domain (WW) sample extended and compact conformations. Ligand binding changes the equilibrium of the interdomain conformations, but the conformational changes that lead to the altered domain sampling were unknown. Prior evidence has supported an interdomain allosteric mechanism. We recently introduced a magnetic resonance-based protocol that allowed us to determine the coupling of intra- and interdomain structural sampling in apo Pin1. Here, we describe ligand-specific conformational changes that occur upon binding of pCDC25c and FFpSPR. pCDC25c binding doubles the population of the extended states compared to the virtually identical populations of the apo and FFpSPR-bound forms. pCDC25c binding to the WW domain triggers conformational changes to propagate via the interdomain interface to the catalytic site, while FFpSPR binding displaces a helix in the PPIase that leads to repositioning of the PPIase catalytic loop.

Molecular Insights into the Intrinsic Dynamics and Their Roles During Catalysis in Pin1 Peptidyl-prolyl Isomerase

Proteins are intrinsically dynamic and change conformations over a wide range of time scales. While the conformational dynamics have been realized to be important for protein functions, e.g., in activity-stability trade-offs, how they play a role during enzyme catalysis has been of debate over decades. By studying Pin1 peptidyl-prolyl isomerase using extensive molecular dynamics simulations, here we discuss how the slow intrinsic dynamics of Pin1 observed in the NMR relaxation dispersion experiment occur and couple to isomerization reactions in molecular detail. In particular, we analyze the angular correlation functions of the backbone N-H bonds and find that slow conformational transitions occur at about the 3₁₀ helix in the apo state. These events at the helical region further affect the residues at about the ligand binding site. Unfolding of this helix leads to a tight hydrogen bond between the helical region and the ligand binding loop, thus forming a stable coiled structure. The helical and coiled structures are found to be characteristic of the Pin1-ligand complex with the ligand in the trans and cis states, respectively. These results indicate that the changes in the slow dynamics of Pin1 by the isomerization reaction occur via the shift in populations of the helical and coiled states, where the balance is dependent on the ligand isomerization states.

Critical Extracellular Ca2+ Dependence of the Binding between PTH1R and a G-Protein Peptide Revealed by MD Simulations

The parathyroid hormone type 1 receptor (PTH1R), a canonical class B GPCR, is regulated by a positive allosteric modulator, extracellular Ca²⁺. Calcium ions prolong the residence time of PTH on the PTH1R, leading to increased receptor activation and duration of cAMP signaling. But the essential mechanism of the allosteric behavior of PTH1R is not fully understood. Here, extensive molecular dynamics (MD) simulations are performed for the PTH1R-G-protein combinations with and without Ca²⁺ to describe how calcium ions allosterically engage receptor-G-protein coupling. We find that the binding of Ca²⁺ stabilizes the conformation of the PTH1R-PTH-spep (the α5 helix of Gs protein) complex, especially the extracellular loop 1 (ECL1). Moreover, the MM-GBSA result indicates that Ca²⁺ allosterically promotes the interaction between PTH1R and spep, consistent with the observation of steered molecular dynamics (SMD) simulations. We further illuminate the possible allosteric signaling pathway from the stable Ca²⁺-coupling site to the intracellular G-protein binding site. These results unveil structural determinants for Ca²⁺ allosterism in the PTH1R-PTH-spep complex and give insights into pluridimensional GPCR signaling regulated by calcium ions.

A Specific pSer/Thr-Pro Motif Generates Interdomain Communication Bifurcations of Two Modes of Pin1 in Solution Nuclear Magnetic Resonance

Peptides mediate the interdomain communication of Pin1 (peptidyl-prolyl cis-trans isomerase) and can regulate its conformation and biochemical functions, providing an idea for drug design using Pin1. Two template peptide sequences have been widely used in the extended or compact state of Pin1 (Cdc25C, E-Q-P-L-pT-P-V-T-D-L; Pintide, W-F-Y-pS-P-R). The way in which specific pSer/Thr-Pro peptides regulate interdomain communication to achieve the opposite state is not clear. In this study, we subdivided the sequence composition of eight types of modified peptides and investigated the interaction with Pin1 by solution nuclear magnetic resonance and molecular dynamics. Demonstrating sequence dependence on the pSer-Pro or pThr-Pro motif and different residues in anchoring the WW domain, the Pin peptide (Pintide, PintideT, Pin25C, and Pin25CT) transmits this concentration accumulation to the PPIase domain, thus exhibiting two anchoring tendencies. However, the Cdc peptide (Cdc25C, Cdc25CS, Cdctide, and CdctideS) has a low binding energy that makes it difficult for the conformation to reach a steady state. In addition, Pin1 is influenced by both compact and extended states, regulated precisely by the sequence as well as by threonine or serine. These results provide new insight into the interdomain communication of Pin1 via pSer/Thr-Pro peptide binding.

Quantitative Analysis of Dynamic Allostery

Dynamic allostery refers to one important class of allosteric regulation that does not involve noticeable conformational changes upon effector binding. In recent years, many "quasi"-dynamic allosteric proteins have been found to only experience subtle conformational changes during allosteric regulation. However, as enthalpic and entropic contributions are coupled to each other and even tiny conformational changes could bring in noticeable free energy changes, a quantitative description is essential to understand the contribution of pure dynamic allostery. Here, by developing a unified anisotropic elastic network model (uANM) considering both side-chain information and ligand heavy atoms, we quantitatively estimated the contribution of pure dynamic allostery in a dataset of known allosteric proteins by excluding the conformational changes upon ligand binding. We found that the contribution of pure dynamic allostery is generally small (much weaker than previously expected) and robustly exhibits an allosteric activation effect, which exponentially decays with the distance between the substrate and the allosteric ligand. We further constructed toy models to study the determinant factors of dynamic allostery in monomeric and oligomeric proteins using the uANM. Analysis of the toy models revealed that a short distance, a small angle between the two ligands, strong protein-ligand interactions, and weak protein internal interactions lead to strong dynamic allostery. Our study provides a quantitative estimation of pure dynamic allostery and facilitates the understanding of dynamic-allostery-controlled biological processes and the design of allosteric drugs and proteins.

Neural relational inference to learn long-range allosteric interactions in proteins from molecular dynamics simulations

Protein allostery is a biological process facilitated by spatially long-range intra-protein communication, whereby ligand binding or amino acid change at a distant site affects the active site remotely. Molecular dynamics (MD) simulation provides a powerful computational approach to probe the allosteric effect. However, current MD simulations cannot reach the time scales of whole allosteric processes. The advent of deep learning made it possible to evaluate both spatially short and long-range communications for understanding allostery. For this purpose, we applied a neural relational inference model based on a graph neural network, which adopts an encoder-decoder architecture to simultaneously infer latent interactions for probing protein allosteric processes as dynamic networks of interacting residues. From the MD trajectories, this model successfully learned the long-range interactions and pathways that can mediate the allosteric communications between distant sites in the Pin1, SOD1, and MEK1 systems. Furthermore, the model can discover allostery-related interactions earlier in the MD simulation trajectories and predict relative free energy changes upon mutations more accurately than other methods.

Mechanistic insight into the destabilization of p53TD tetramer by cancer-related R337H mutation: a molecular dynamics study

The p53 protein is a tumor suppressor crucial for cell cycle and genome integrity. In a very large proportion of human cancers, p53 is frequently inactivated by mutations located in its DNA-binding domain (DBD). Some experimental studies reported that the inherited R337H mutation located in the p53 tetramerization domain (p53TD) can also result in destabilization of the p53 protein, and consequently lead to an organism prone to cancer setup. However, the underlying R337H mutation-induced structural destabilization mechanism is not well understood. Herein, we investigate the structural stability and dynamic property of the wild type p53TD tetramer and its cancer-related R337H mutant by performing multiple microsecond molecular dynamics simulations. It is found that R337H mutation destroys the R337-D352 hydrogen bonds, weakens the F341-F341 π-π stacking interaction and the hydrophobic interaction between aliphatic hydrocarbons of R337 and M340, leading to more solvent exposure of all the hydrophobic cores, and thus disrupting the structural integrity of the tetramer. Importantly, our simulations show for the first time that R337H mutation results in unfolding of the α-helix starting from the N-terminal region (residues ³³⁵RER(H)FEM³⁴⁰). Consistently, community network analyses reveal that R337H mutation reduces dynamical correlation and global connectivity of p53TD tetramer, which destabilizes the structure of the p53TD tetramer. This study provides the atomistic mechanism of R337H mutation-induced destabilization of p53TD tetramer, which might be helpful for in-depth understanding of the p53 loss-of-function mechanism.

Reconstruction of Coupled Intra- and Interdomain Protein Motion from Nuclear and Electron Magnetic Resonance

Proteins composed of multiple domains allow for structural heterogeneity and interdomain dynamics that may be vital for function. Intradomain structures and dynamics can influence interdomain conformations and vice versa. However, no established structure determination method is currently available that can probe the coupling of these motions. The protein Pin1 contains separate regulatory and catalytic domains that sample "extended" and "compact" states, and ligand binding changes this equilibrium. Ligand binding and interdomain distance have been shown to impact the activity of Pin1, suggesting interdomain allostery. In order to characterize the conformational equilibrium of Pin1, we describe a novel method to model the coupling between intra- and interdomain dynamics at atomic resolution using multistate ensembles. The method uses time-averaged nuclear magnetic resonance (NMR) restraints and double electron-electron resonance (DEER) data that resolve distance distributions. While the intradomain calculation is primarily driven by exact nuclear Overhauser enhancements (eNOEs), J couplings, and residual dipolar couplings (RDCs), the relative domain distribution is driven by paramagnetic relaxation enhancement (PREs), RDCs, interdomain NOEs, and DEER. Our data support a 70:30 population of the compact and extended states in apo Pin1. A multistate ensemble describes these conformations simultaneously, with distinct conformational differences located in the interdomain interface stabilizing the compact or extended states. We also describe correlated conformations between the catalytic site and interdomain interface that may explain allostery driven by interdomain contact.

Substitutions at Nonconserved Rheostat Positions Modulate Function by Rewiring Long-Range, Dynamic Interactions

Amino acid substitutions at nonconserved protein positions can have noncanonical and "long-distance" outcomes on protein function. Such outcomes might arise from changes in the internal protein communication network, which is often accompanied by changes in structural flexibility. To test this, we calculated flexibilities and dynamic coupling for positions in the linker region of the lactose repressor protein. This region contains nonconserved positions for which substitutions alter DNA-binding affinity. We first chose to study 11 substitutions at position 52. In computations, substitutions showed long-range effects on flexibilities of DNA-binding positions, and the degree of flexibility change correlated with experimentally measured changes in DNA binding. Substitutions also altered dynamic coupling to DNA-binding positions in a manner that captured other experimentally determined functional changes. Next, we broadened calculations to consider the dynamic coupling between 17 linker positions and the DNA-binding domain. Experimentally, these linker positions exhibited a wide range of substitution outcomes: Four conserved positions tolerated hardly any substitutions ("toggle"), ten nonconserved positions showed progressive changes from a range of substitutions ("rheostat"), and three nonconserved positions tolerated almost all substitutions ("neutral"). In computations with wild-type lactose repressor protein, the dynamic couplings between the DNA-binding domain and these linker positions showed varied degrees of asymmetry that correlated with the observed toggle/rheostat/neutral substitution outcomes. Thus, we propose that long-range and noncanonical substitutions outcomes at nonconserved positions arise from rewiring long-range communication among functionally important positions. Such calculations might enable predictions for substitution outcomes at a range of nonconserved positions.

Conformational dynamics is critical for the allosteric inhibition of cGAS upon acetyl-mimic mutations

Detection of cytosolic dsDNA by cyclic GMP-AMP synthase (cGAS) is critical for the immune system to sense and fight against infection, but chronic activation of cGAS by self-DNA leads to autoimmune diseases without effective treatment yet. It was found that acetylation on either Lys384, Lys394, or Lys414 could inhibit the catalytic production of cGAMP by cGAS, and further suppressed self-DNA-induced autoimmunity. However, the implied mechanism remains unclear. Here, extensive molecular dynamics simulations combined with multiple analytical approaches were employed to uncover the allosteric inhibition mechanisms by using the K-to-Q mutations to mimic acetylation. Results suggested that the exterior loops contributed most to the conformational dynamics of cGAS, and two concerted intrinsic motions were observed: the inward/outward or twisting movement for the outer appendage of lobe 1 and the open/closed swing of the active-site loops. Mutations slightly affected the binding of dsDNA and cGAMP. The shift of the conformational sampling of the active-site loops or residues around cGAMP upon mutation might potentially explain the inhibition of cGAS activity. Moreover, the intra- and inter-molecular coupling was weakened upon mutations more or less but via distinct pathways. Hence, conformational dynamics play a vital role in the allosteric inhibition of cGAS upon the studied acetyl-mimic mutations. As the studied acetyl-mimic mutations are located at either the inter-lobe or inter-molecular interfaces, hence except for acetylation, our findings might help the development of new therapeutics against autoimmune diseases due to abnormal cGAS activation by designing inter-lobe or intermolecular allosteric inhibitors.

Deciphering the T790M/L858R-Selective Inhibition Mechanism of an Allosteric Inhibitor of EGFR: Insights from Molecular Simulations

Allosteric inhibitors have lately received great attention because of their unique advantages, representing a more suitable choice for combinatory therapeutics targeting resistance-relevant signaling cascades. Among the various inhibitors, an allosteric small-molecule inhibitor, JBJ-04-125-02, has been proven to be effective against EGFR^T790M/L858R mutant in vivo and in vitro. Herein, an in silico approach was adopted to shed light on the deep understanding of the higher selectivity of JBJ-04-125-02 against EGFR^T790M/L858R mutant than wild-type EGFR. Our results indicate that JBJ-04-125-02 prefers to bind with the EGFR^T790M/L858R mutant, stabilizes the inactive conformation, and further allosterically affects the conformations and dynamics of the interlobe cleft, including both the allosteric site and the ATP-binding site. Furthermore, docking results confirm that the binding of JBJ-04-125-02 at the allosteric site decreases the binding affinity of ANP (an ATP analogue) at the orthosteric site, especially for the Mut-holo one, which might further inhibit the function of EGFR. The present work provides a clear picture of the mutant-selective inhibition mechanism of an allosteric inhibitor of EGFR. The findings might pave the way for designing allosteric drugs targeting EGFR mutant lung cancer patients, which also takes a step forward in terms of drug resistance caused by protein mutations.

Protein Allostery at Atomic Resolution

Protein allostery is a phenomenon involving the long range coupling between two distal sites in a protein. In order to elucidate allostery at atomic resoluion on the ligand-binding WW domain of the enzyme Pin1, multistate structures were calculated from exact nuclear Overhauser effect (eNOE). In its free form, the protein undergoes a microsecond exchange between two states, one of which is predisposed to interact with its parent catalytic domain. In presence of the positive allosteric ligand, the equilibrium between the two states is shifted towards domain-domain interaction, suggesting a population shift model. In contrast, the allostery-suppressing ligand decouples the side-chain arrangement at the inter-domain interface thereby reducing the inter-domain interaction. As such, this mechanism is an example of dynamic allostery. The presented distinct modes of action highlight the power of the interplay between dynamics and function in the biological activity of proteins.

Identification of potential allosteric binding sites in cathepsin K based on intramolecular communication

Network theory methods and molecular dynamics (MD) simulations are accepted tools to study allosteric regulation. Indeed, dynamic networks built upon correlation analysis of MD trajectories provide detailed information about communication paths between distant sites. In this context, we aimed to understand whether the efficiency of intramolecular communication could be used to predict the allosteric potential of a given site. To this end, we performed MD simulations and network theory analyses in cathepsin K (catK), whose allosteric sites are well defined. To obtain a quantitative measure of the efficiency of communication, we designed a new protocol that enables the comparison between properties related to ensembles of communication paths obtained from different sites. Further, we applied our strategy to evaluate the allosteric potential of different catK cavities not yet considered for drug design. Our predictions of the allosteric potential based on intramolecular communication correlate well with previous catK experimental and theoretical data. We also discuss the possibility of applying our approach to other proteins from the same family.

Allostery and Epistasis: Emergent Properties of Anisotropic Networks

Understanding the underlying mechanisms behind protein allostery and non-additivity of substitution outcomes (i.e., epistasis) is critical when attempting to predict the functional impact of mutations, particularly at non-conserved sites. In an effort to model these two biological properties, we extend the framework of our metric to calculate dynamic coupling between residues, the Dynamic Coupling Index (DCI) to two new metrics: (i) EpiScore, which quantifies the difference between the residue fluctuation response of a functional site when two other positions are perturbed with random Brownian kicks simultaneously versus individually to capture the degree of cooperativity of these two other positions in modulating the dynamics of the functional site and (ii) DCI_asym, which measures the degree of asymmetry between the residue fluctuation response of two sites when one or the other is perturbed with a random force. Applied to four independent systems, we successfully show that EpiScore and DCI_asym can capture important biophysical properties in dual mutant substitution outcomes. We propose that allosteric regulation and the mechanisms underlying non-additive amino acid substitution outcomes (i.e., epistasis) can be understood as emergent properties of an anisotropic network of interactions where the inclusion of the full network of interactions is critical for accurate modeling. Consequently, mutations which drive towards a new function may require a fine balance between functional site asymmetry and strength of dynamic coupling with the functional sites. These two tools will provide mechanistic insight into both understanding and predicting the outcome of dual mutations.

PIN1 Provides Dynamic Control of MYC in Response to Extrinsic Signals

PIN1 is a phosphorylation-directed member of the peptidyl-prolyl cis/trans isomerase (PPIase) family that facilitates conformational changes in phosphorylated targets such as c-MYC (MYC). Following signaling events that mediate phosphorylation of MYC at Serine 62, PIN1 establishes structurally distinct pools of MYC through its trans-cis and cis-trans isomerization activity at Proline 63. Through these isomerization steps, PIN1 functionally regulates MYC's stability, the molecular timing of its DNA binding and transcriptional activity, and its subnuclear localization. Recently, our group showed that Serine 62 phosphorylated MYC can associate with the inner basket of the nuclear pore (NP) in a PIN1-dependent manner. The poised euchromatin at the NP basket enables rapid cellular response to environmental signals and cell stress, and PIN1-mediated trafficking of MYC calibrates this response. In this perspective, we describe the molecular aspects of PIN1 target recognition and PIN1's function in the context of its temporal and spatial regulation of MYC.

Integrated Computational Approaches and Tools forAllosteric Drug Discovery

Understanding molecular mechanisms underlying the complexity of allosteric regulationin proteins has attracted considerable attention in drug discovery due to the benefits and versatilityof allosteric modulators in providing desirable selectivity against protein targets while minimizingtoxicity and other side effects. The proliferation of novel computational approaches for predictingligand-protein interactions and binding using dynamic and network-centric perspectives has ledto new insights into allosteric mechanisms and facilitated computer-based discovery of allostericdrugs. Although no absolute method of experimental and in silico allosteric drug/site discoveryexists, current methods are still being improved. As such, the critical analysis and integration ofestablished approaches into robust, reproducible, and customizable computational pipelines withexperimental feedback could make allosteric drug discovery more efficient and reliable. In this article,we review computational approaches for allosteric drug discovery and discuss how these tools can beutilized to develop consensus workflows for in silico identification of allosteric sites and modulatorswith some applications to pathogen resistance and precision medicine. The emerging realization thatallosteric modulators can exploit distinct regulatory mechanisms and can provide access to targetedmodulation of protein activities could open opportunities for probing biological processes and insilico design of drug combinations with improved therapeutic indices and a broad range of activities.

Detecting Functional Dynamics in Proteins with Comparative Perturbed-Ensembles Analysis

Recent advances have made all-atom molecular dynamics (MD) a powerful tool to sample the conformational energy landscape. There are still however three major challenges in the application of MD to biological systems: accuracy of force field, time scale, and the analysis of simulation trajectories. Significant progress in addressing the first two challenges has been made and extensively reviewed previously. This Account focuses on strategies of analyzing simulation data of biomolecules that also covers ways to properly design simulations and validate simulation results. In particular, we examine an approach named comparative perturbed-ensembles analysis, which we developed to efficiently detect dynamics in protein MD simulations that can be linked to biological functions. In our recent studies, we implemented this approach to understand allosteric regulations in several disease-associated human proteins. The central task of a comparative perturbed-ensembles analysis is to compare two or more conformational ensembles of a system generated by MD simulations under distinct perturbation conditions. Perturbations can be different sequence variations, ligand-binding conditions, and other physical/chemical modifications of the system. Each simulation is long enough (e.g., microsecond-long) to ensure sufficient sampling of the local substate. Then, sophisticated bioinformatic and statistical tools are applied to extract function-related information from the simulation data, including principal component analysis, residue-residue contact analysis, difference contact network analysis (dCNA) based on the graph theory, and statistical analysis of side-chain conformations. Computational findings are further validated with experimental data. By comparing distinct conformational ensembles, functional micro- to millisecond dynamics can be inferred. In contrast, such a time scale is difficult to reach in a single simulation; even when reached for a single condition of a system, it is elusive as to what dynamical motions are related to functions without, for example, comparing free and substrate-bound proteins at the minimum. We illustrate our approach with three examples. First, we discuss using the approach to identify allosteric pathways in cyclophilin A (CypA), a member of a ubiquitous class of peptidyl-prolyl cis-trans isomerase enzymes. By comparing side-chain torsion-angle distributions of CypA in wild-type and mutant forms, we identified three pathways: two are consistent with recent nuclear magnetic resonance experiments, whereas the third is a novel pathway. Second, we show how the approach enables a dynamical-evolution analysis of the human cyclophilin family. In the analysis, both conserved and divergent conformational dynamics across three cyclophilin isoforms (CypA, CypD, and CypE) were summarized. The conserved dynamics led to the discovery of allosteric networks resembling those found in CypA. A residue wise determinant underlying the unique dynamics in CypD was also detected and validated with additional mutational MD simulations. In the third example, we applied the approach to elucidate a peptide sequence-dependent allosteric mechanism in human Pin 1, a phosphorylation-dependent peptidyl-prolyl isomerase. We finally present our outlook of future directions. Especially, we envisage how the approach could help open a new avenue in drug discovery.

Mutational effects of Cys113 on structural dynamics of Pin1

Pin1 is a peptidyl-prolyl isomerase (PPIase) which catalyzes cis/trans isomerization of pS/pT-P bond. Its activity is related to various cellular functions including suppression of Alzheimer's disease. A cysteine residue C113 is known to be important for its PPIase activity; a mutation C113A reduced the activity by 130-fold. According to various nuclear magnetic resonance experiments for mutants of C113 and molecular dynamics (MD) simulation of wild-type Pin1, the protonation sate of Sγ of C113 regulates the hydrogen-bonding network of the dual-histidine motif (H59, H157) whose dynamics may affect substrate binding ability. However, it was still unclear why such local dynamic changes altered the PPIase activity of Pin1. In this study, we performed 500 ns of MD simulations of full-length wild-type Pin1 and C113A mutant in order to elucidate why the mutation C113A drastically reduced the PPIase activity of Pin1. The principal component analysis for both MD trajectories clearly elucidated that the mutation C113A suppressed the dynamics of Pin1 because it stabilized a hydrogen-bond between Nɛ of H59 and Oγ of S115. In the dynamics of wild-type protein, the phosphate binding loop (K63-S71) as well as the interdomain hinge showed the closed-open dynamics which correlated with the change of the hydrogen-bonding network of the dual-histidine motif. In contrast, in the dynamics of C113A mutant, the phosphate binding loop took only the closed conformation together with the interdomain hinge. Such closed-open dynamics must be essential for the PPIase activity of Pin1.

Probing Protein Allostery as a Residue-Specific Concept via Residue Response Maps

Allosteric regulation is a well-established phenomenon defined as a distal conformational or dynamical change of the protein upon allosteric effector binding. Here, we developed a novel approach to delineate allosteric effects in proteins. In this approach, we applied robust machine learning methods, including deep neural network and random forest, on extensive molecular dynamics (MD) simulations to distinguish otherwise similar allosteric states of proteins. Using the PDZ3 domain of PDS-95 as a model protein, we demonstrated that the allosteric effects could be represented as residue-specific properties through two-dimensional property-residue maps, which we refer to as "residue response maps". These maps were constructed through two machine learning methods and could accurately describe how different properties of various residues are affected upon allosteric perturbation on protein. Based on the "residue response maps", we propose allostery as a residue-specific concept, suggesting that all residues could be considered as allosteric residues because each residue "senses" the allosteric events through changing its single or multiple attributes in a quantitatively unique way. The "residue response maps" could be used to fingerprint a protein based on the unique patterns of residue responses upon binding events, providing a novel way to systematically describe the protein allosteric effects of each residue upon perturbation.

Pyruvate Kinase Regulates the Pentose-Phosphate Pathway in Response to Hypoxia in Mycobacterium tuberculosis

In response to the stress of infection, Mycobacterium tuberculosis (Mtb) reprograms its metabolism to accommodate nutrient and energetic demands in a changing environment. Pyruvate kinase (PYK) is an essential glycolytic enzyme in the phosphoenolpyruvate-pyruvate-oxaloacetate node that is a central switch point for carbon flux distribution. Here we show that the competitive binding of pentose monophosphate inhibitors or the activator glucose 6-phosphate (G6P) to MtbPYK tightly regulates the metabolic flux. Intriguingly, pentose monophosphates were found to share the same binding site with G6P. The determination of a crystal structure of MtbPYK with bound ribose 5-phosphate (R5P), combined with biochemical analyses and molecular dynamic simulations, revealed that the allosteric inhibitor pentose monophosphate increases PYK structural dynamics, weakens the structural network communication, and impairs substrate binding. G6P, on the other hand, primes and activates the tetramer by decreasing protein flexibility and strengthening allosteric coupling. Therefore, we propose that MtbPYK uses these differences in conformational dynamics to up- and down-regulate enzymic activity. Importantly, metabolome profiling in mycobacteria reveals a significant increase in the levels of pentose monophosphate during hypoxia, which provides insights into how PYK uses dynamics of the tetramer as a competitive allosteric mechanism to retard glycolysis and facilitate metabolic reprogramming toward the pentose-phosphate pathway for achieving redox balance and an anticipatory metabolic response in Mtb.

Tuning the "violin" of protein kinases: The role of dynamics-based allostery

The intricacies of allosteric regulation of protein kinases continue to engage the research community. Allostery, or control from a distance, is seen as a fundamental biomolecular mechanism for proteins. From the traditional methods of conformational selection and induced fit, the field has grown to include the role of protein motions in defining a dynamics-based allosteric approach. Harnessing of these continuous motions in the protein to exert allosteric effects can be defined by a "violin" model that focuses on distributions of protein vibrations as opposed to concerted pathways. According to this model, binding of an allosteric modifier causes global redistribution of dynamics in the protein kinase domain that leads to changes in its catalytic properties. This model is consistent with the "entropy-driven allostery" mechanism proposed by Cooper and Dryden in 1984 and does not require, but does not exclude, any major structural changes. We provide an overview of practical implementation of the violin model and how it stands amidst the other known models of protein allostery. Protein kinases have been described as the biomolecules of interest. © 2019 IUBMB Life, 71(6):685-696, 2019.

Disordered linkers in multidomain allosteric proteins: Entropic effect to favor the open state or enhanced local concentration to favor the closed state?

There are many multidomain allosteric proteins where an allosteric signal at the allosteric domain modifies the activity of the functional domain. Intrinsically disordered regions (linkers) are widely involved in this kind of regulation process, but the essential role they play therein is not well understood. Here, we investigated the effect of linkers in stabilizing the open or the closed states of multidomain proteins using combined thermodynamic deduction and coarse-grained molecular dynamics simulations. We revealed that the influence of linker can be fully characterized by an effective local concentration [B]₀ . When K_d is smaller than [B]₀ , the closed state would be favored; while the open state would be preferred when K_d is larger than [B]₀ . We used four protein systems with markedly different domain-domain binding affinity and structural order/disorder as model systems to understand the relationship between [B]₀ and the linker length as well as its flexibility. The linker length is the main practical determinant of [B]₀ . [B]₀ of a flexible linker with 40-60 residues was determined to be in a narrow range of 0.2-0.6 mM, while a too short or too long length would dramatically decrease [B]₀ . With the revealed [B]₀ range, the introduction of a flexible linker makes the regulation of weakly interacting partners possible.

Allostery in Its Many Disguises: From Theory to Applications

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

Integration of network models and evolutionary analysis into high-throughput modeling of protein dynamics and allosteric regulation: theory, tools and applications

Proteins are dynamical entities that undergo a plethora of conformational changes, accomplishing their biological functions. Molecular dynamics simulation and normal mode analysis methods have become the gold standard for studying protein dynamics, analyzing molecular mechanism and allosteric regulation of biological systems. The enormous amount of the ensemble-based experimental and computational data on protein structure and dynamics has presented a major challenge for the high-throughput modeling of protein regulation and molecular mechanisms. In parallel, bioinformatics and systems biology approaches including genomic analysis, coevolution and network-based modeling have provided an array of powerful tools that complemented and enriched biophysical insights by enabling high-throughput analysis of biological data and dissection of global molecular signatures underlying mechanisms of protein function and interactions in the cellular environment. These developments have provided a powerful interdisciplinary framework for quantifying the relationships between protein dynamics and allosteric regulation, allowing for high-throughput modeling and engineering of molecular mechanisms. Here, we review fundamental advances in protein dynamics, network theory and coevolutionary analysis that have provided foundation for rapidly growing computational tools for modeling of allosteric regulation. We discuss recent developments in these interdisciplinary areas bridging computational biophysics and network biology, focusing on promising applications in allosteric regulations, including the investigation of allosteric communication pathways, protein–DNA/RNA interactions and disease mutations in genomic medicine. We conclude by formulating and discussing future directions and potential challenges facing quantitative computational investigations of allosteric regulatory mechanisms in protein systems.

Uncorrelated Effect of Interdomain Contact on Pin1 Isomerase Activity Reveals Positive Catalytic Cooperativity

Pin1 is a two-domain peptidyl-prolyl isomerase (PPIase) associated with neurodegeneration and tumorigenesis. The two domains, a WW and a PPIase domain, are connected by a flexible linker, making Pin1 adopt various conformations ranging from compact to extended, wherein Pin1 exhibits different extents of interdomain contact. Previous studies have shown that weakening interdomain contact increases the isomerase activity of Pin1. Here, we propose an NMR chemical shift correlation-analysis-based method that will be general for two-domain proteins to gauge two-state populations of Pin1, and we report a linker-modified mutant of Pin1 with enhanced interdomain contact and increased isomerase activity, with the latter suggesting an uncorrelated effect of interdomain contact on isomerase activity. Thus, although bindings of different substrates in the WW domain impose opposite effects on interdomain contact, in both cases, it may promote isomerization, implying cooperativity between substrate binding in the WW domain and isomerization in the PPIase domain.

Molecular Mechanism of the Pin1-Histone H1 Interaction

Pin1 is an essential peptidyl-prolyl isomerase (PPIase) that catalyzes cis-trans prolyl isomerization in proteins containing pSer/Thr-Pro motifs. It has an N-terminal WW domain that targets these motifs and a C-terminal PPIase domain that catalyzes isomerization. Recently, Pin1 was shown to modify the conformation of phosphorylated histone H1 and stabilize the chromatin-H1 interaction by increasing its residence time. This Pin1-histone H1 interaction plays a key role in pathogen response, in infection, and in cell cycle control; therefore, anti-Pin1 therapeutics are an important focus for treating infections as well as cancer. Each of the H1 histones (H1.0-H1.5) contains several potential Pin1 recognition pSer/pThr-Pro motifs. To understand the Pin1-histone H1 interaction fully, we investigated how both the isolated WW domain and full-length Pin1 interact with three H1 histone substrate peptide sequences that were previously identified as important binding partners (H1.1, H1.4, and H1.5). NMR spectroscopy was used to measure the binding affinities and the interdomain dynamics upon binding to these sequences. We observed different K_D values depending on the histone binding site, suggesting that energetics play a role in guiding the Pin1-histone interaction. While interdomain interactions vary between the peptides, we find no evidence for allosteric activation for the histone H1 substrates.

Gears-In-Motion: The Interplay of WW and PPIase Domains in Pin1

Pin1 belongs to the family of the peptidyl-prolyl cis-trans isomerase (PPIase), which is a class of enzymes that catalyze the cis/trans isomerization of the Proline residue. Pin1 is unique and only catalyzes the phosphorylated Serine/Threonine-Proline (S/T-P) motifs of a subset of proteins. Since the discovery of Pin1 as a key protein in cell cycle regulation, it has been implicated in numerous diseases, ranging from cancer to neurodegenerative diseases. The main features of Pin1 lies in its two main domains: the WW (two conserved tryptophan) domain and the PPIase domain. Despite extensive studies trying to understand the mechanisms of Pin1 functions, how these two domains contribute to the biological roles of Pin1 in cellular signaling requires more investigations. The WW domain of Pin1 is known to have a higher affinity to its substrate than that of the PPIase domain. Yet, the WW domain seems to prefer the trans configuration of phosphorylated S/T-P motif, while the PPIase catalyzes the cis to trans isomerasion. Such contradicting information has generated much confusion as to the actual mechanism of Pin1 function. In addition, dynamic allostery has been suggested to be important for Pin1 function. Henceforth, in this review, we will be looking at the progress made in understanding the function of Pin1, and how these understandings can aid us in overcoming the diseases implicated by Pin1 such as cancer during drug development.

Molecular mechanism of 15-lipoxygenase allosteric activation and inhibition

Human reticulocyte 15-lipoxygenase (15-LOX) plays an important role in inflammation resolution and is also involved in many cancer-related processes. Both an activator and an inhibitor will serve as research tools for understanding the biological functions of 15-LOX and provide opportunities for drug discovery. In a previous study, both allosteric activators and inhibitors of 15-LOX were discovered through a virtual screening based computational approach. However, why molecules binding to the same site causes different effects remains to be disclosed. In the present study, we used previously reported activator and inhibitor molecules as probes to elucidate the mechanism of allosteric regulation of 15-LOX. We measured the influences of the allosteric activator and inhibitor on the enzymatic reaction rate and found that the activator increases 15-LOX activity by preventing substrate inhibition instead of increasing the turnover number. The inhibitor can also prevent substrate inhibition but decreases the turnover number at the same time, resulting in inhibition. Molecular dynamics simulations were conducted to help explain the underlying mechanism of allostery. Both the activator and inhibitor were demonstrated to be able to prevent 15-LOX from transforming into potentially inactive conformations. Compared to the activator, the inhibitor molecule restrains the motions of residues around the substrate binding site and reduces the flexibility of 15-LOX. These results explained the different effects between the activator and the inhibitor and shed light on how to effectively design novel activator molecules.

Same but not alike: Structure, flexibility and energetics of domains in multi-domain proteins are influenced by the presence of other domains

The majority of the proteins encoded in the genomes of eukaryotes contain more than one domain. Reasons for high prevalence of multi-domain proteins in various organisms have been attributed to higher stability and functional and folding advantages over single-domain proteins. Despite these advantages, many proteins are composed of only one domain while their homologous domains are part of multi-domain proteins. In the study presented here, differences in the properties of protein domains in single-domain and multi-domain systems and their influence on functions are discussed. We studied 20 pairs of identical protein domains, which were crystallized in two forms (a) tethered to other proteins domains and (b) tethered to fewer protein domains than (a) or not tethered to any protein domain. Results suggest that tethering of domains in multi-domain proteins influences the structural, dynamic and energetic properties of the constituent protein domains. 50% of the protein domain pairs show significant structural deviations while 90% of the protein domain pairs show differences in dynamics and 12% of the residues show differences in the energetics. To gain further insights on the influence of tethering on the function of the domains, 4 pairs of homologous protein domains, where one of them is a full-length single-domain protein and the other protein domain is a part of a multi-domain protein, were studied. Analyses showed that identical and structurally equivalent functional residues show differential dynamics in homologous protein domains; though comparable dynamics between in-silico generated chimera protein and multi-domain proteins were observed. From these observations, the differences observed in the functions of homologous proteins could be attributed to the presence of tethered domain. Overall, we conclude that tethered domains in multi-domain proteins not only provide stability or folding advantages but also influence pathways resulting in differences in function or regulatory properties.

High prevalence of multi-domain proteins in proteomes has been attributed to higher stability and functional and folding advantages of the multi-domain proteins. Influence of tethering of domains on the overall properties of proteins has been well studied but its influence on the properties of the constituent domains is largely unaddressed. Here, we investigate the influence of tethering of domains in multi-domain proteins on the structural, dynamics and energetics properties of the constituent domains and its implications on the functions of proteins. To this end, comparative analyses were carried out for identical protein domains crystallized in tethered and untethered forms. Also, comparative analyses of single-domain proteins and their homologous multi-domain proteins were performed. The analyses suggest that tethering influences the structural, dynamic and energetic properties of constituent protein domains. Our observations hint at regulation of protein domains by tethered domains in multi-domain systems, which may manifest at the differential function observed between single-domain and homologous multi-domain proteins.

Hinge-Shift Mechanism Modulates Allosteric Regulations in Human Pin1

Allostery, which is regulation from distant sites, plays a major role in biology. While traditional allostery is described in terms of conformational change upon ligand binding as an underlying principle, it is possible to have allosteric regulations without significant conformational change through modulating the conformational dynamics by altering the local effective elastic modulus of the protein upon ligand binding. Pin1 utilizes this dynamic allostery to regulate its function. It is a modular protein containing a WW domain and a larger peptidyl prolyl isomerase domain (PPIase) that isomerizes phosphoserine/threonine-proline (pS/TP) motifs. The WW domain serves as a docking module, whereas catalysis solely takes place within the PPIase domain. Here, we analyze the change in the dynamic flexibility profile of the PPIase domain upon ligand binding to the WW domain. Substrate binding to the WW domain induces the formation of a new rigid hinge site around the interface of the two domains and loosens the flexibility of a rigid site existing in the Apo form around the catalytic site. This hinge-shift mechanism enhances the dynamic coupling of the catalytic positions with the PPIase domain, where the rest of the domain can cooperatively respond to the local conformational changes around the catalytic site, leading to an increase in catalytic efficiency.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.