Revealing Atomic-Level Mechanisms of Protein Allostery with Molecular Dynamics Simulations

Exploring Binding Sites in Chagas Disease Protein TcP21 Using Integrated Mixed Solvent Molecular Dynamics Approaches

Chagas disease, caused by the protozoan Trypanosoma cruzi, remains a significant global health burden, particularly in Latin America, where millions are at risk. This disease predominantly affects socioeconomically vulnerable populations, aggravating economic inequality, marginalization, and low political visibility. Despite extensive research, effective treatments are still lacking, partly due to the complex biology of the parasite and its infection mechanisms. This study focuses on TcP21, a novel 21 kDa protein secreted by extracellular amastigotes, which has been implicated in T. cruzi infection via an alternative infective pathway. Although the potential of TcP21 for understanding Chagas disease is promising, further exploration is necessary, particularly in identifying potential binding sites on its surface. Computational tools offer a versatile and effective strategy for preliminary binding site assessment, facilitating a more cost-efficient allocation of experimental resources. In this study, we employed three independent computational approaches─mixed solvent molecular dynamics simulations (MSMD), fragment-based molecular docking, and pharmacophore model docking coupled with molecular dynamics simulations─to identify potential binding sites and provide comprehensive insights into TcP21. The three methodologies converged on a common site located on the external surface of the protein, characterized by key residues such as GLU55, ASP52, VAL70, ILE62, and TRP77. The protonated amino, acetamido, and phenyl groups of the pharmacophore probe were consistently observed to interact with the site via a network of salt bridges, hydrogen bonds, charge-charge interactions, and alkyl-π interactions, suggesting these groups play a significant role in ligand binding. This study does not aim to propose specific therapeutic hits but to highlight a still unknown and unexplored protein involved in T. cruzi cell invasion. In this regard, given the strong correlation between the three distinct approaches used for mapping, we consider this study offers valuable insights for further research into P21 and its role in Chagas disease.

FDA-approved drugs as PIM-1 kinase inhibitors: A drug repurposed approach for cancer therapy

PIM-1 kinase, a member of the Serine/Threonine kinase family, has emerged as a promising therapeutic target in various cancers due to its role in promoting tumor growth and resistance to conventional therapies. In this study, we employed a structure-based approach to screen 3800 FDA-approved drugs to discover potential inhibitors of PIM-1. After an initial selection of 50 candidates based on high docking scores, four drugs, stanozolol, alfaxalone, rifaximin, and telmisartan, were identified as strong PIM-1 binders, interacting with key residues in the ATP-binding pocket of the kinase. To assess the stability of these interactions, we conducted all-atom molecular dynamic simulations, confirming favorable dynamics. Experimental validation via a kinase inhibition assay on recombinant PIM-1 showed that rifaximin significantly inhibited PIM-1 activity, with an IC50 of ∼26 μM. Fluorescence binding assays further demonstrated a strong binding affinity for rifaximin, with a binding constant, corroborated by isothermal titration calorimetry studies. Our findings suggest that rifaximin may serve as a potential repurposed drug for targeting PIM-1 in cancer treatment. However, further validations are required in a clinical setting before the final therapeutic implications.

Multiple Functional Protein-Protein Interaction Interfaces Allosterically Regulate ATP-Binding in Cyclin-Dependent Kinase-1

The ATP-dependent phosphorylation activity of cyclin-dependent kinase 1 (CDK1), an essential enzyme for cell cycle progression, is regulated by interactions with Cyclin-B, substrate, and Cks proteins. We have recently shown that active site acetylation in CDK1 abrogated binding to Cyclin-B which posits an intriguing long-range communication between the catalytic site and the protein-protein interaction (PPI) interface. Now, we demonstrate a general allosteric link between the CDK1 active site and all three of its PPI interfaces through atomistic molecular dynamics (MD) simulations. Specifically, we examined ATP binding free energies to CDK1 in native nonacetylated (K33wt) and acetylated (K33Ac) forms as well as the acetyl-mimic K33Q and the acetyl-null K33R mutant forms, which are accessible in vitro. In agreement with experiments, ATP binding is stronger in K33wt relative to the other three perturbed states. Free energy decomposition reveals, in addition to expected local changes, significant and selective nonlocal entropic responses to ATP binding/perturbation of K33 from theαC -helix, activation loop (A-loop), andαG -α H segments in CDK1 which interface with Cyclin-B, substrate, and Cks proteins, respectively. Statistical analysis reveals that while entropic responses of protein segments to active site perturbations are on average correlated with their dynamical changes, such correlations are lost in about 9%-48% of the dataset depending on the segment. Besides proving the bi-directional communication between the active site and the CDK1:Cyclin-B interface, our study uncovers a hitherto unknown mode of ATP binding regulation by multiple PPI interfaces in CDK1.

Special Issue: "Molecular Dynamics Simulations and Structural Analysis of Protein Domains"

The 3D protein structure is the basis for all their biological functions [...].

Computational insights into the mechanisms underlying structural destabilization and recovery in trafficking-deficient hERG mutants

Cardiovascular diseases are a major global health concern, responsible for a significant number of deaths each year, often linked to cardiac arrhythmias resulting from dysfunction in ion channels. Hereditary Long QT Syndrome (LQTS) is a condition characterized by a prolonged QT interval on ECG, increasing the risk of sudden cardiac death. The most common type of LQTS, LQT2, is caused by mutations in the hERG gene, affecting a potassium ion channel. The majority of these mutations disrupt the channel's trafficking to the cell membrane, leading to intracellular retention. Specific high-affinity hERG blockers (e.g., E-4031) can rescue this mutant phenotype, but the exact mechanism is unknown. This study used accelerated molecular dynamics simulations to investigate how these mutations affect the hERG channel's structure, folding, endoplasmic reticulum (ER) retention, and trafficking. We reveal that these mutations induce structural changes in the channel, narrowing its central pore and altering the conformation of the intracellular domains. These changes expose internalization signals that contribute to ER retention and degradation of the mutant hERG channels. Moreover, the study found that the trafficking rescue drug E-4031 can inhibit these structural changes, potentially rescuing the mutant channels. This research offers valuable insights into the structural issues responsible for the degradation of rescuable transmembrane trafficking mutants. Understanding the defective trafficking structure of the hERG channel could help identify binding sites for small molecules capable of restoring proper folding and facilitating channel trafficking. This knowledge has the potential to lead to mechanism-based therapies that address the condition at the cellular level, which may prove more effective than treating clinical symptoms, ultimately offering hope for individuals with hereditary Long QT Syndrome.

Leveraging Bidirectional Nature of Allostery To Inhibit Protein-Protein Interactions (PPIs): A Case Study of PCSK9-LDLR Interaction

The protein PCSK9 (proprotein convertase subtilisin/Kexin type 9) negatively regulates the recycling of LDLR (low-density lipoprotein receptor), leading to an elevated plasma level of LDL. Inhibition of PCSK9-LDLR interaction has emerged as a promising therapeutic strategy to manage hypercholesterolemia. However, the large interaction surface area between PCSK9 and LDLR makes it challenging to identify a small molecule competitive inhibitor. An alternative strategy would be to identify distal cryptic sites as targets for allosteric inhibitors that can remotely modulate PCSK9-LDLR interaction. Using several microseconds long molecular dynamics (MD) simulations, we demonstrate that on binding with LDLR, there is a significant conformational change (population shift) in a distal loop (residues 211-222) region of PCSK9. Consistent with the bidirectional nature of allostery, we establish a clear correlation between the loop conformation and the binding affinity with LDLR. Using a thermodynamic argument, we establish that the loop conformations predominantly present in the apo state of PCSK9 would have lower LDLR binding affinity, and they would be potential targets for designing allosteric inhibitors. We elucidate the molecular origin of the allosteric coupling between this loop and the LDLR binding interface in terms of the population shift in a set of salt bridges and hydrogen bonds. Overall, our work provides a general strategy toward identifying allosteric hotspots: compare the conformational ensemble of the receptor between the apo and bound states of the protein and identify distal conformational changes, if any. The inhibitors should be designed to bind and stabilize the apo-specific conformations.

G protein-coupled receptors (GPCRs): advances in structures, mechanisms, and drug discovery

G protein-coupled receptors (GPCRs), the largest family of human membrane proteins and an important class of drug targets, play a role in maintaining numerous physiological processes. Agonist or antagonist, orthosteric effects or allosteric effects, and biased signaling or balanced signaling, characterize the complexity of GPCR dynamic features. In this study, we first review the structural advancements, activation mechanisms, and functional diversity of GPCRs. We then focus on GPCR drug discovery by revealing the detailed drug-target interactions and the underlying mechanisms of orthosteric drugs approved by the US Food and Drug Administration in the past five years. Particularly, an up-to-date analysis is performed on available GPCR structures complexed with synthetic small-molecule allosteric modulators to elucidate key receptor-ligand interactions and allosteric mechanisms. Finally, we highlight how the widespread GPCR-druggable allosteric sites can guide structure- or mechanism-based drug design and propose prospects of designing bitopic ligands for the future therapeutic potential of targeting this receptor family.

Mechanism of drug resistance in HIV-1 protease subtype C in the presence of Atazanavir

AIDS is one of the deadliest diseases in the history of humankind caused by HIV. Despite the technological development, curtailing the viral infection inside human host still remains a challenge. Therapies such as HAART uses a combination of drugs to inhibit the viral activity. One of the important targets includes HIV protease and inhibiting its activity will minimize the production of mature structural proteins. However, the genetic diversity and the occurrence of drug resistant mutations adds complexity to effective drug design. In this study, we aimed at understanding the drug binding mechanism of one such subtype, namely subtype C and its insertion variant L38HL. We performed multiple molecular dynamics simulations along with binding free energy analysis of wild-type and L38HL bound to Atazanavir (ATV). From the analysis, we revealed that the insertion alters the hydrogen bond and hydrophobic interaction networks. The alterations in the interaction networks increase flexibility at the hinge-fulcrum interface. Further, the effects of these changes affect flap tip curling. Moreover, the changes in the hinge-fulcrum-cantilever interface alters the concerted motion of the functional regions leading to change in the direction of flap movement thus causing a subtle change in the active site volume. Additionally, formation of intramolecular hydrogen bonds in the ATV docked to L38HL restricted the movement of R1 and R2 groups thereby altering the interactions. Overall, the changes in the flexibility of flap together with the changes in the active site volume and compactness of the ligand provide insights for increased binding affinity of ATV with L38HL.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

Allosteric regulation of kinase activity in living cells

The dysregulation of protein kinases is associated with multiple diseases due to the kinases' involvement in a variety of cell signaling pathways. Manipulating protein kinase function, by controlling the active site, is a promising therapeutic and investigative strategy to mitigate and study diseases. Kinase active sites share structural similarities, making it difficult to specifically target one kinase, and allosteric control allows specific regulation and study of kinase function without directly targeting the active site. Allosteric sites are distal to the active site but coupled via a dynamic network of inter-atomic interactions between residues in the protein. Establishing an allosteric control over a kinase requires understanding the allosteric wiring of the protein. Computational techniques offer effective and inexpensive mapping of the allosteric sites on a protein. Here, we discuss the methods to map and regulate allosteric communications in proteins, and strategies to establish control over kinase functions in live cells and organisms. Protein molecules, or 'sensors,' are engineered to function as tools to control allosteric activity of the protein as these sensors have high spatiotemporal resolution and help in understanding cell phenotypes after immediate activation or inactivation of a kinase. Traditional methods used to study protein functions, such as knockout, knockdown, or mutation, cannot offer a sufficiently high spatiotemporal resolution. We discuss the modern repertoire of tools to regulate protein kinases as we enter a new era in deciphering cellular signaling and developing novel approaches to treat diseases associated with signal dysregulation.

Structural basis of dimerization of chemokine receptors CCR5 and CXCR4

G protein-coupled receptors (GPCRs) are prominent drug targets responsible for extracellular-to-intracellular signal transduction. GPCRs can form functional dimers that have been poorly characterized so far. Here, we show the dimerization mechanism of the chemokine receptors CCR5 and CXCR4 by means of an advanced free-energy technique named coarse-grained metadynamics. Our results reproduce binding events between the GPCRs occurring in the minute timescale, revealing a symmetric and an asymmetric dimeric structure for each of the three investigated systems, CCR5/CCR5, CXCR4/CXCR4, and CCR5/CXCR4. The transmembrane helices TM4-TM5 and TM6-TM7 are the preferred binding interfaces for CCR5 and CXCR4, respectively. The identified dimeric states differ in the access to the binding sites of the ligand and G protein, indicating that dimerization may represent a fine allosteric mechanism to regulate receptor activity. Our study offers structural basis for the design of ligands able to modulate the formation of CCR5 and CXCR4 dimers and in turn their activity, with therapeutic potential against HIV, cancer, and immune-inflammatory diseases.

Allosteric regulation of kinase activity in living cells

The dysregulation of protein kinases is associated with multiple diseases due to the kinases' involvement in a variety of cell signaling pathways. Manipulating protein kinase function, by controlling the active site, is a promising therapeutic and investigative strategy to mitigate and study diseases. Kinase active sites share structural similarities making it difficult to specifically target one kinase, allosteric control allows specific regulation and study of kinase function without directly targeting the active site. Allosteric sites are distal to the active site but coupled via a dynamic network of inter-atomic interactions between residues in the protein. Establishing an allosteric control over a kinase requires understanding the allosteric wiring of the protein. Computational techniques offer effective and inexpensive mapping of the allosteric sites on a protein. Here, we discuss methods to map and regulate allosteric communications in proteins, and strategies to establish control over kinase functions in live cells and organisms. Protein molecules, or "sensors" are engineered to function as tools to control allosteric activity of the protein as these sensors have high spatiotemporal resolution and help in understanding cell phenotypes after immediate activation or inactivation of a kinase. Traditional methods used to study protein functions, such as knockout, knockdown, or mutation, cannot offer a sufficiently high spatiotemporal resolution. We discuss the modern repertoire of tools to regulate protein kinases as we enter a new era in deciphering cellular signaling and developing novel approaches to treat diseases associated with signal dysregulation.

Structure-driven protein engineering for production of valuable natural products

Proteins are the most frequently used biocatalysts, and their structures determine their functions. Modifying the functions of proteins on the basis of their structures lies at the heart of protein engineering, opening a new horizon for metabolic engineering by efficiently generating stable enzymes. Many attempts at classical metabolic engineering have focused on improving specific metabolic fluxes and producing more valuable natural products by increasing gene expression levels and enzyme concentrations. However, most naturally occurring enzymes show limitations, and such limitations have hindered practical applications. Here we review recent advances in protein engineering in synthetic biology, chemoenzymatic synthesis, and plant metabolic engineering and describe opportunities for designing and constructing novel enzymes or proteins with desirable properties to obtain more active natural products.

Allosteric Signaling in PDZ Energetic Networks: Embedding Error Analysis

Allosteric signaling in proteins has been known for some half a century, yet how the signal traverses the protein remains an active area of research. Recently, the importance of electrostatics to achieve long-range signaling has become increasingly appreciated. Our laboratory has been working on developing network approaches to capture such interactions. In this study, we turn our attention to the well-studied allosteric model protein, PDZ. We study the allosteric dynamics on a per-residue basis in key constructs involving the PDZ domain, its allosteric effector, and its peptide ligand. We utilize molecular dynamics trajectories to create the networks for the constructs to explore the allosteric effect by plotting the heat kernel results onto axes defined by principal components. We introduce a new metric to quantitate the volume sampled by a residue in the latent space. We relate our findings to PDZ and the greater field of allostery.

Variable Regions of p53 Isoforms Allosterically Hard Code DNA Interaction

Allosteric regulation of protein activity pervades biology as the "second secret of life." We have been examining the allosteric regulation and mutant reactivation of the tumor suppressor protein p53. We have found that generalizing the definition of allosteric effector to include entire proteins and expanding the meaning of binding site to include the interface of a transcription factor with its DNA to be useful in understanding the modulation of protein activity. Here, we cast the variable regions of p53 isoforms as allosteric regulators of p53 interactions with its consensus DNA. We implemented molecular dynamics simulations and our lab's new techniques of molecular dynamics (MD) sectors and MD-Markov state models to investigate the effects of nine naturally occurring splice variant isoforms of p53. We find that all of the isoforms differ from wild type in their dynamic properties and how they interact with the DNA. We consider the implications of these findings on allostery and cancer treatment.

Emerging Methods and Applications to Decrypt Allostery in Proteins and Nucleic Acids

Many large protein-nucleic acid complexes exhibit allosteric regulation. In these systems, the propagation of the allosteric signaling is strongly coupled to conformational dynamics and catalytic function, challenging state-of-the-art analytical methods. Here, we review established and innovative approaches used to elucidate allosteric mechanisms in these complexes. Specifically, we report network models derived from graph theory and centrality analyses in combination with molecular dynamics (MD) simulations, introducing novel schemes that implement the synergistic use of graph theory with enhanced simulations methods and ab-initio MD. Accelerated MD simulations are used to construct "enhanced network models", describing the allosteric response over long timescales and capturing the relation between allostery and conformational changes. "Ab-initio network models" combine graph theory with ab-initio MD and quantum mechanics/molecular mechanics (QM/MM) simulations to describe the allosteric regulation of catalysis by following the step-by-step dynamics of biochemical reactions. This approach characterizes how the allosteric regulation changes from reactants to products and how it affects the transition state, revealing a tense-to-relaxed allosteric regulation along the chemical step. Allosteric models and applications are showcased for three paradigmatic examples of allostery in protein-nucleic acid complexes: (i) the nucleosome core particle, (ii) the CRISPR-Cas9 genome editing system and (iii) the spliceosome. These methods and applications create innovative protocols to determine allosteric mechanisms in protein-nucleic acid complexes that show tremendous promise for medicine and bioengineering.

Insights into Phosphorylation-Induced Protein Allostery and Conformational Dynamics of Glycogen Phosphorylase via Integrative Structural Mass Spectrometry and In Silico Modeling

Allosteric regulation plays a fundamental role in innumerable biological processes. Understanding its dynamic mechanism and impact at the molecular level is of great importance in disease diagnosis and drug discovery. Glycogen phosphorylase (GP) is a phosphoprotein responding to allosteric regulation and has significant biological importance to glycogen metabolism. Although the atomic structures of GP have been previously solved, the conformational dynamics of GP related to allostery regulation remain largely elusive due to its macromolecular size (∼196 kDa). Here, we integrated native top-down mass spectrometry (nTD-MS), hydrogen-deuterium exchange MS (HDX-MS), protection factor (PF) analysis, molecular dynamics (MD) simulations, and allostery signaling analysis to examine the structural basis and dynamics for the allosteric regulation of GP by phosphorylation. nTD-MS reveals differences in structural stability as well as oligomeric state between the unphosphorylated (GPb) and phosphorylated (GPa) forms. HDX-MS, PF analysis, and MD simulations further pinpoint the structural differences between GPb and GPa involving the binding interfaces (the N-terminal and tower-tower helices), catalytic site, and PLP-binding region. More importantly, it also allowed us to complete the missing link of the long-range communication process from the N-terminal tail to the catalytic site caused by phosphorylation. This integrative MS and in silico-based platform is highly complementary to biophysical approaches and yields valuable insights into protein structures and dynamic regulation.

Application of molecular dynamics simulation in biomedicine

Molecular dynamics (MD) simulation has been widely used in the field of biomedicine to study the conformational transition of proteins caused by mutation or ligand binding/unbinding. It provides some perspectives those are difficult to find in traditional biochemical or pathological experiments, for example, detailed effects of mutations on protein structure and protein-protein/ligand interaction at the atomic level. In this review, a broad overview on conformation changes and drug discovery by MD simulation is given. We first discuss the preparation of protein structure for MD simulation, which is a key step that determines the accuracy of the simulation. Then, we summarize the applications of commonly used force fields and MD simulations in scientific research. Finally, enhanced sampling methods and common applications of these methods are introduced. In brief, MD simulation is a powerful tool and it can be used to guide experimental study. The combination of MD simulation and experimental techniques is an a priori means to solve the biomedical problems and give a deep understanding on the relationship between protein structure and function.

Integration of machine learning with computational structural biology of plants

Computational structural biology of proteins has developed rapidly in recent decades with the development of new computational tools and the advancement of computing hardware. However, while these techniques have widely been used to make advancements in human medicine, these methods have seen less utilization in the plant sciences. In the last several years, machine learning methods have gained popularity in computational structural biology. These methods have enabled the development of new tools which are able to address the major challenges that have hampered the wide adoption of the computational structural biology of plants. This perspective examines the remaining challenges in computational structural biology and how the development of machine learning techniques enables more in-depth computational structural biology of plants.

SenseNet, a tool for analysis of protein structure networks obtained from molecular dynamics simulations

Computational methods play a key role for investigating allosteric mechanisms in proteins, with the potential of generating valuable insights for innovative drug design. Here we present the SenseNet (“Structure ENSEmble NETworks”) framework for analysis of protein structure networks, which differs from established network models by focusing on interaction timelines obtained by molecular dynamics simulations. This approach is evaluated by predicting allosteric residues reported by NMR experiments in the PDZ2 domain of hPTP1e, a reference system for which previous computational predictions have shown considerable variance. We applied two models based on the mutual information between interaction timelines to estimate the conformational influence of each residue on its local environment. In terms of accuracy our prediction model is comparable to the top performing model published for this system, but by contrast benefits from its independence from NMR structures. Our results are complementary to experimental data and the consensus of previous predictions, demonstrating the potential of our new analysis tool SenseNet. Biochemical interpretation of our model suggests that allosteric residues in the PDZ2 domain form two distinct clusters of contiguous sidechain surfaces. SenseNet is provided as a plugin for the network analysis software Cytoscape, allowing for ease of future application and contributing to a system of compatible tools bridging the fields of system and structural biology.

Molecular Mechanism for the Allosteric Inhibition of the Human Serotonin Transporter by Antidepressant Escitalopram

Human serotine transporter (hSERT) is one of the most influential drug targets, and its allosteric modulators (e.g., escitalopram) have emerged to be the next-generation medication for psychiatric disorders. However, the molecular mechanism underlying the allosteric modulation of hSERT is still elusive. Here, the simulation strategies of conventional (cMD) and steered (SMD) molecular dynamics were applied to investigate this molecular mechanism from distinct perspectives. First, cMD simulations revealed that escitalopram's binding to hSERT's allosteric site simultaneously enhanced its binding to the orthosteric site. Then, SMD simulation identified that the occupation of hSERT's allosteric site by escitalopram could also block its dissociation from the orthosteric site. Finally, by comparing the simulated structures of two hSERT-escitalopram complexes with and without allosteric modulation, a new conformational coupling between an extracellular (Arg104-Glu494) and an intracellular (Lys490-Glu494) salt bridge was identified. In summary, this study explored the mechanism underlying the allosteric modulation of hSERT by collectively applying two MD simulation strategies, which could facilitate our understanding of the allosteric modulations of not only hSERT but also other clinically important therapeutic targets.

Correlated Motions in Structural Biology

Correlated motions in proteins arising from the collective movements of residues have long been proposed to be fundamentally important to key properties of proteins, from allostery and catalysis to evolvability. Recent breakthroughs in structural biology have made it possible to capture proteins undergoing complex conformational changes, yet intrinsic correlated motions within a conformation remain one of the least understood facets of protein structure. For many decades, the analysis of total X-ray scattering held the promise of animating crystal structures with correlated motions. With recent advances in both X-ray detectors and data interpretation methods, this long-held promise can now be met. In this Perspective, we will introduce how correlated motions are captured in total scattering and provide guidelines for the collection, interpretation, and validation of data. As structural biology continues to push the boundaries, we see an opportunity to gain atomistic insight into correlated motions using total scattering as a bridge between theory and experiment.

Dynamic allostery controls the peptide sensitivity of the Ly49C natural killer receptor

Using a variety of activating and inhibitory receptors, natural killer (NK) cells protect against disease by eliminating cells that have downregulated class I major histocompatibility complex (MHC) proteins, such as in response to cell transformation or viral infection. The inhibitory murine NK receptor Ly49C specifically recognizes the class I MHC protein H-2K^b. Unusual among NK receptors, Ly49C exhibits a peptide-dependent sensitivity to H-2K^b recognition, which has not been explained despite detailed structural studies. To gain further insight into Ly49C peptide sensitivity, we examined Ly49C recognition biochemically and through the lens of dynamic allostery. We found that the peptide sensitivity of Ly49C arises through small differences in H-2K^b-binding affinity. Although molecular dynamics simulations supported a role for peptide-dependent protein dynamics in producing these differences in binding affinity, calorimetric measurements indicated an enthalpically as opposed to entropically driven process. A quantitative linkage analysis showed that this emerges from peptide-dependent dynamic tuning of electrostatic interactions across the Ly49C–H-2K^b interface. We propose a model whereby different peptides alter the flexibility of H-2K^b, which in turn changes the strength of electrostatic interactions across the protein–protein interface. Our results provide a quantitative assessment of how peptides alter Ly49C-binding affinity, suggest the underlying mechanism, and demonstrate peptide-driven allostery at work in class I MHC proteins. Lastly, our model provides a solution for how dynamic allostery could impact binding of some, but not all, class I MHC partners depending on the structural and chemical composition of the interfaces.

A novel RNA pol II CTD interaction site on the mRNA capping enzyme is essential for its allosteric activation

Recruitment of the mRNA capping enzyme (CE/RNGTT) to the site of transcription is essential for the formation of the 5' mRNA cap, which in turn ensures efficient transcription, splicing, polyadenylation, nuclear export and translation of mRNA in eukaryotic cells. The CE GTase is recruited and activated by the Serine-5 phosphorylated carboxyl-terminal domain (CTD) of RNA polymerase II. Through the use of molecular dynamics simulations and enhanced sampling techniques, we provide a systematic and detailed characterization of the human CE-CTD interface, describing the effect of the CTD phosphorylation state, length and orientation on this interaction. Our computational analyses identify novel CTD interaction sites on the human CE GTase surface and quantify their relative contributions to CTD binding. We also identify, for the first time, allosteric connections between the CE GTase active site and the CTD binding sites, allowing us to propose a mechanism for allosteric activation. Through binding and activity assays we validate the novel CTD binding sites and show that the CDS2 site is essential for CE GTase activity stimulation. Comparison of the novel sites with cocrystal structures of the CE-CTD complex in different eukaryotic taxa reveals that this interface is considerably more conserved than previous structures have indicated.

Exploring the putative mechanism of allosteric modulations by mixed-action kappa/mu opioid receptor bitopic modulators

The modulation and selectivity mechanisms of seven mixed-action kappa opioid receptor (KOR)/mu opioid receptor (MOR) bitopic modulators were explored. Molecular modeling results indicated that the 'message' moiety of seven bitopic modulators shared the same binding mode with the orthosteric site of the KOR and MOR, whereas the 'address' moiety bound with different subdomains of the allosteric site of the KOR and MOR. The 'address' moiety of seven bitopic modulators bound to different subdomains of the allosteric site of the KOR and MOR may exhibit distinguishable allosteric modulations to the binding affinity and/or efficacy of the 'message' moiety. Moreover, the 3-hydroxy group on the phenolic moiety of the seven bitopic modulators induced selectivity to the KOR over the MOR.

Optogenetic Analysis of Allosteric Control in Protein Tyrosine Phosphatases

Allosteric regulation enables dynamic adjustments to protein function that permit tight control over cellular biochemistry. Discrepancies in the allosteric systems of related proteins can thus reveal important differences in their susceptibilities to influential stimuli (e.g., allosteric ligands, mutations, or post-translational modifications). This study uses an optogenetic actuator as a tool to compare the allosteric systems of two structurally related regulatory proteins: protein tyrosine phosphatase 1B (PTP1B) and T-cell protein tyrosine phosphatase (TCPTP). It begins with an interesting observation: The fusion of a protein light switch to the allosterically influential α7 helix of PTP1B permits optical modulation of its catalytic activity, but a similar fusion to TCPTP does not. A subsequent analysis of different PTP chimeras shows that replacing regions of TCPTP with homologous regions from PTP1B can enhance photocontrol; as TCPTP becomes more "PTP1B-like", its photosensitivity increases. Interestingly, the structural changes required for photocontrol also enhance the sensitivity of TCPTP to other allosteric inputs, notably, an allosteric inhibitor and a newly reported activating mutation. Our findings indicate that the allosteric functionality of the α7 helix of PTP1B is not conserved across the PTP family and highlight residues necessary to transfer this functionality to other PTPs. More broadly, our results suggest that simple gene fusion events can strengthen allosteric communication within individual protein domains and describe an intriguing application for optogenetic actuators as structural probes-a sort of physically disruptive "ratchet"-for studying protein allostery.

Synergistic Allostery in Multiligand-Protein Interactions

Amide hydrogen-deuterium exchange mass spectrometry is powerful for describing combinatorial coupling effects of a cooperative ligand pair binding at noncontiguous sites: adenosine at the ATP-pocket and a docking peptide (PIFtide) at the PIF-pocket, on a model protein kinase PDK1. Binding of two ligands to PDK1 reveal multiple hotspots of synergistic allostery with cumulative effects greater than the sum of individual effects mediated by each ligand. We quantified this synergism and ranked these hotspots using a difference in deuteration-based approach, which showed that the strongest synergistic effects were observed at three of the critical catalytic loci of kinases: the αB-αC helices, and HRD-motif loop, and DFG-motif. Additionally, we observed weaker synergistic effects at a distal GHI-subdomain locus. Synergistic changes in deuterium exchange observed at a distal site but not at the intermediate sites of the large lobe of the kinase reveals allosteric propagation in proteins to operate through two modes. Direct electrostatic interactions between polar and charged amino acids that mediate targeted relay of allosteric signals, and diffused relay of allosteric signals through soft matter-like hydrophobic core amino acids. Furthermore, we provide evidence that the conserved β-3 strand lysine of protein kinases (Lys111 of PDK1) functions as an integrator node to coordinate allosteric coupling of the two ligand-binding sites. It maintains indirect interactions with the ATP-pocket and mediates a critical salt bridge with a glutamate (Glu130) of αC helix, which is conserved across all kinases. In summary, allosteric propagation in cooperative, dual-liganded enzyme targets is bidirectional and synergistic and offers a strategy for combinatorial drug development.

Engineering and application of a biosensor with focused ligand specificity

Cell factories converting bio-based precursors to chemicals present an attractive avenue to a sustainable economy, yet screening of genetically diverse strain libraries to identify the best-performing whole-cell biocatalysts is a low-throughput endeavor. For this reason, transcriptional biosensors attract attention as they allow the screening of vast libraries when used in combination with fluorescence-activated cell sorting (FACS). However, broad ligand specificity of transcriptional regulators (TRs) often prohibits the development of such ultra-high-throughput screens. Here, we solve the structure of the TR LysG of Corynebacterium glutamicum, which detects all three basic amino acids. Based on this information, we follow a semi-rational engineering approach using a FACS-based screening/counterscreening strategy to generate an l-lysine insensitive LysG-based biosensor. This biosensor can be used to isolate l-histidine-producing strains by FACS, showing that TR engineering towards a more focused ligand spectrum can expand the scope of application of such metabolite sensors.

Transcriptional biosensors represent powerful tools for the screening of vast strain libraries, but the broad ligand specificity of some transcriptional regulators (TRs) can prohibit such applications. Here authors present the engineering of a LysG-based biosensor with a focused ligand specificity to isolate L-histidine-producing strains.

How Do Molecular Dynamics Data Complement Static Structural Data of GPCRs

G protein-coupled receptors (GPCRs) are implicated in nearly every physiological process in the human body and therefore represent an important drug targeting class. Advances in X-ray crystallography and cryo-electron microscopy (cryo-EM) have provided multiple static structures of GPCRs in complex with various signaling partners. However, GPCR functionality is largely determined by their flexibility and ability to transition between distinct structural conformations. Due to this dynamic nature, a static snapshot does not fully explain the complexity of GPCR signal transduction. Molecular dynamics (MD) simulations offer the opportunity to simulate the structural motions of biological processes at atomic resolution. Thus, this technique can incorporate the missing information on protein flexibility into experimentally solved structures. Here, we review the contribution of MD simulations to complement static structural data and to improve our understanding of GPCR physiology and pharmacology, as well as the challenges that still need to be overcome to reach the full potential of this technique.

Allosteric Regulation at the Crossroads of New Technologies: Multiscale Modeling, Networks, and Machine Learning

Allosteric regulation is a common mechanism employed by complex biomolecular systems for regulation of activity and adaptability in the cellular environment, serving as an effective molecular tool for cellular communication. As an intrinsic but elusive property, allostery is a ubiquitous phenomenon where binding or disturbing of a distal site in a protein can functionally control its activity and is considered as the "second secret of life." The fundamental biological importance and complexity of these processes require a multi-faceted platform of synergistically integrated approaches for prediction and characterization of allosteric functional states, atomistic reconstruction of allosteric regulatory mechanisms and discovery of allosteric modulators. The unifying theme and overarching goal of allosteric regulation studies in recent years have been integration between emerging experiment and computational approaches and technologies to advance quantitative characterization of allosteric mechanisms in proteins. Despite significant advances, the quantitative characterization and reliable prediction of functional allosteric states, interactions, and mechanisms continue to present highly challenging problems in the field. In this review, we discuss simulation-based multiscale approaches, experiment-informed Markovian models, and network modeling of allostery and information-theoretical approaches that can describe the thermodynamics and hierarchy allosteric states and the molecular basis of allosteric mechanisms. The wealth of structural and functional information along with diversity and complexity of allosteric mechanisms in therapeutically important protein families have provided a well-suited platform for development of data-driven research strategies. Data-centric integration of chemistry, biology and computer science using artificial intelligence technologies has gained a significant momentum and at the forefront of many cross-disciplinary efforts. We discuss new developments in the machine learning field and the emergence of deep learning and deep reinforcement learning applications in modeling of molecular mechanisms and allosteric proteins. The experiment-guided integrated approaches empowered by recent advances in multiscale modeling, network science, and machine learning can lead to more reliable prediction of allosteric regulatory mechanisms and discovery of allosteric modulators for therapeutically important protein targets.

Revisiting allostery in CREB-binding protein (CBP) using residue-based interaction energy

CREB-binding protein (CBP) is a multi-subunit scaffold protein complex in transcription regulation process, binding and interacting with ligands such as mixed-lineage leukemia (MLL) and c-Myb allosterically. Here in this study, we have revisited the concept of allostery in CBP via residue-based interaction energy calculation based on molecular dynamics (MD) simulations. To this end, we conducted MD simulations of KIX:MLL:c-Myb ternary complex, its binary components and kinase-inducible domain (KID) interacting domain (KIX) backbone. Interaction energy profiles and cross correlation analysis were performed and the results indicated that KIX:MLL and KIX:c-Myb:MLL complexes demonstrate significant similarities according to both analysis methods. Two regions in the KIX backbone were apparent from the interaction energy and cross correlation maps that hold a key to allostery phenomena observed in CBP. While one of these regions are related to the ligand binding residues, the other comprises of L₁₂-G₂ loop and α₃ helix regions that have been found to have a significant role in allosteric signal propagation. All in all, residue-based interaction energy calculation method is demonstrated to be a valuable calculation technique for the detection of allosteric signal propagation and ligand interaction regions.

Spectral analysis of molecular dynamics simulations on PDZ: MD sectors

The idea of protein "sectors" posits that sparse subsets of amino acid residues form cooperative networks that are key elements of protein stability, ligand binding, and allosterism. To date, protein sectors have been calculated by the statistical coupling analysis (SCA) method of Ranganathan and co-workers via the spectral analysis of conservation-weighted evolutionary covariance matrices obtained from a multiple sequence alignments of homologous families of proteins. SCA sectors, a knowledge-based protocol, have been indentified with functional properties and allosterism for a number of systems. In this study, we investigate the utility of the sector idea for the analysis of physics-based molecular dynamics (MD) trajectories of proteins. Our test case for this procedure is PSD95- PDZ3, one of the smallest proteins for which allosterism has been observed. It has served previously as a model system for a number of prediction algorithms, and is well characterized by X-ray crystallography, NMR spectroscopy and site specific mutagenisis. All-atom MD simulations were performed for a total of 500 nanoseconds using AMBER, and MD-calculated covariance matrices for the fluctuations of residue displacements and non-bonded interaction energies were subjected to spectral analysis in a manner analogous to that of SCA. The composition of MD sectors was compared with results from SCA, site specific mutagenesis, and allosterism. The concordance indicates that MD sectors are a viable protocol for analyzing MD trajectories and provide insight into the physical origin of the phenomenon.Communicated by Ramaswamy H. Sarma.

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.

Mechanism of allosteric activation of human mRNA cap methyltransferase (RNMT) by RAM: insights from accelerated molecular dynamics simulations

The RNA guanine-N7 methyltransferase (RNMT) in complex with RNMT-activating miniprotein (RAM) catalyses the formation of a N7-methylated guanosine cap structure on the 5' end of nascent RNA polymerase II transcripts. The mRNA cap protects the primary transcript from exonucleases and recruits cap-binding complexes that mediate RNA processing, export and translation. By using microsecond standard and accelerated molecular dynamics simulations, we provide for the first time a detailed molecular mechanism of allosteric regulation of RNMT by RAM. We show that RAM selects the RNMT active site conformations that are optimal for binding of substrates (AdoMet and the cap), thus enhancing their affinity. Furthermore, our results strongly suggest the likely scenario in which the cap binding promotes the subsequent AdoMet binding, consistent with the previously suggested cooperative binding model. By employing the network community analyses, we revealed the underlying long-range allosteric networks and paths that are crucial for allosteric regulation by RAM. Our findings complement and explain previous experimental data on RNMT activity. Moreover, this study provides the most complete description of the cap and AdoMet binding poses and interactions within the enzyme's active site. This information is critical for the drug discovery efforts that consider RNMT as a promising anti-cancer target.

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.

Data-driven computational analysis of allosteric proteins by exploring protein dynamics, residue coevolution and residue interaction networks

BACKGROUND

Computational studies of allosteric interactions have witnessed a recent renaissance fueled by the growing interest in modeling of the complex molecular assemblies and biological networks. Allosteric interactions in protein structures allow for molecular communication in signal transduction networks.

METHODS

In this work, we performed a large scale comprehensive and multi-faceted analysis of >300 diverse allosteric proteins and complexes with allosteric modulators. By modeling and exploring coarse-grained dynamics, residue coevolution, and residue interaction networks for allosteric proteins, we have determined unifying molecular signatures shared by allosteric systems.

RESULTS

The results of this study have suggested that allosteric inhibitors and allosteric activators may differentially affect global dynamics and network organization of protein systems, leading to diverse allosteric mechanisms. By using structural and functional data on protein kinases, we present a detailed case study that that included atomic-level analysis of coevolutionary networks in kinases bound with allosteric inhibitors and activators.

CONCLUSIONS

We have found that coevolutionary networks can form direct communication pathways connecting functional regions and can recapitulate key regulatory sites and interactions responsible for allosteric signaling in the studied protein systems. The results of this computational investigation are compared with the experimental studies and reveal molecular signatures of known regulatory hotspots in protein kinases.

GENERAL SIGNIFICANCE

This study has shown that allosteric inhibitors and allosteric activators can have a different effect on residue interaction networks and can exploit distinct regulatory mechanisms, which could open up opportunities for probing allostery and new drug combinations with broad range of activities.

Allostery in G protein-coupled receptors investigated by molecular dynamics simulations

G-protein-coupled receptors (GPCRs) are allosteric signaling machines that trigger distinct functional responses depending on the particular conformational state they adopt upon binding. This so-called GPCR functional selectivity is prompted by ligands of different efficacy binding at orthosteric or allosteric sites on the receptor, as well as by interactions with intracellular protein partners or other receptor types. Molecular dynamics (MD) simulations can provide important mechanistic, thermodynamic, and kinetic insights into these interactions at a level of molecular detail that is necessary to rightly inform modern drug discovery. Here, we review the most recent MD contributions to understanding GPCR allostery, with an emphasis on their strengths and limitations.

The Influence of Chemical Change on Protein Dynamics: A Case Study with Pyruvate Formate-Lyase

Pyruvate formate-lyase (PFL) catalyzes the reversible conversion of pyruvate and coenzyme A (CoA) into formate and acetyl-CoA in two half-reactions. For the second half-reaction to take place, the S-H group of CoA must enter the active site of the enzyme to retrieve a protein-bound acetyl group. However, CoA is bound at the protein surface, whereas the active site is buried in the protein interior, some 20-30 Å away. The PFL system was therefore subjected to a series of extensive molecular dynamics simulations (in the μs range) and a host of advanced analysis procedures. Models representing PFL before and after the first half-reaction were used to examine the possible effect of enzyme acetylation. All simulated structures were found to be relatively stable compared to the initial crystal structure. Although the adenine portion of CoA remained predominantly bound at the protein surface, the binding of the S-H group was significantly more labile. A potential entry channel for CoA, which would allow the S-H group to reach the active site, was identified and characterized. The channel was found to be associated with accentuated fluctuations and a higher probability of being in an open state in acetylated systems. This result suggests that the acetylation of the enzyme assumes a prominent functional role, whereby the formation of the acyl intermediate serves to initiate a subtle signaling cascade that influences the protein dynamics and facilitates the entry of the second substrate.

Atomistic Modeling of the ABL Kinase Regulation by Allosteric Modulators Using Structural Perturbation Analysis and Community-Based Network Reconstruction of Allosteric Communications

In this work, we have examined the molecular mechanisms of allosteric regulation of the ABL tyrosine kinase at the atomic level. Atomistic modeling of the ABL complexes with a panel of allosteric modulators has been performed using a combination of molecular dynamics simulations, structural residue perturbation scanning, and a novel community analysis of the residue interaction networks. Our results have indicated that allosteric inhibitors and activators may exert a differential control on allosteric signaling between the kinase binding sites and functional regions. While the inhibitor binding can strengthen the closed ABL state and induce allosteric communications directed from the allosteric pocket to the ATP binding site, the DPH activator may induce a more dynamic open form and activate allosteric couplings between the ATP and substrate binding sites. By leveraging a network-centric theoretical framework, we have introduced a novel community analysis method and global topological parameters that have unveiled the hierarchical modularity and the intercommunity bridging sites in the residue interaction network. We have found that allosteric functional hotspots responsible for the kinase regulation may serve the intermodular bridges in the global interaction network. The central conclusion from this analysis is that the regulatory switch centers play a fundamental role in the modular network organization of ABL as the unique intercommunity bridges that connect the SH2 and SH3 domains with the catalytic core into a functional kinase assembly. The hierarchy of network organization in the ABL regulatory complexes may allow for the synergistic action of dense intercommunity links required for the robust signal transfer in the catalytic core and sparse network bridges acting as the regulatory control points that orchestrate allosteric transitions between the inhibited and active kinase forms.

Dissecting a novel allosteric mechanism of cruzain: A computer-aided approach

Trypanosoma cruzi is the causative agent of Chagas disease, a neglected infection affecting millions of people in tropical regions. There are several chemotherapeutic agents for the treatment of this disease, but most of them are highly toxic and generate resistance. Currently, the development of allosteric inhibitors constitutes a promising research field, since it can improve the accessibility to more selective and less toxic medicines. To date, the allosteric drugs prediction is a state-of-the-art topic in rational structure-based computational design. In this work, a simulation strategy was developed for computational discovery of allosteric inhibitors, and it was applied to cruzain, a promising target and the major cysteine protease of T. cruzi. Molecular dynamics simulations, binding free energy calculations and network-based modelling of residue interactions were combined to characterize and compare molecular distinctive features of the apo form and the cruzain-allosteric inhibitor complexes. By using geometry-based criteria on trajectory snapshots, we predicted two main allosteric sites suitable for drug targeting. The results suggest dissimilar mechanisms exerted by the same allosteric site when binding different potential allosteric inhibitors. Finally, we identified the residues involved in suboptimal paths linking the identified site and the orthosteric site. The present study constitutes the first approximation to the design of cruzain allosteric inhibitors and may serve for future pharmacological intervention. Here, no major effects on active site structure were observed due to compound binding (modification of distance and angles between catalytic residues), which indicates that allosteric regulation in cruzain might be mediated via alterations of its dynamical properties similarly to allosteric inhibition of human cathepsin K (HCatK). The current findings are particularly relevant for the design of allosteric modulators of papain-like cysteine proteases.

Fluctuation correlations as major determinants of structure- and dynamics-driven allosteric effects

Both structure- and dynamics-driven allosteric effects are determined by the correlation of distance fluctuations in proteins.

Polyphosphoinositide-Binding Domains: Insights from Peripheral Membrane and Lipid-Transfer Proteins

Within eukaryotic cells, biochemical reactions need to be organized on the surface of membrane compartments that use distinct lipid constituents to dynamically modulate the functions of integral proteins or influence the selective recruitment of peripheral membrane effectors. As a result of these complex interactions, a variety of human pathologies can be traced back to improper communication between proteins and membrane surfaces; either due to mutations that directly alter protein structure or as a result of changes in membrane lipid composition. Among the known structural lipids found in cellular membranes, phosphatidylinositol (PtdIns) is unique in that it also serves as the membrane-anchored precursor of low-abundance regulatory lipids, the polyphosphoinositides (PPIn), which have restricted distributions within specific subcellular compartments. The ability of PPIn lipids to function as signaling platforms relies on both non-specific electrostatic interactions and the selective stereospecific recognition of PPIn headgroups by specialized protein folds. In this chapter, we will attempt to summarize the structural diversity of modular PPIn-interacting domains that facilitate the reversible recruitment and conformational regulation of peripheral membrane proteins. Outside of protein folds capable of capturing PPIn headgroups at the membrane interface, recent studies detailing the selective binding and bilayer extraction of PPIn species by unique functional domains within specific families of lipid-transfer proteins will also be highlighted. Overall, this overview will help to outline the fundamental physiochemical mechanisms that facilitate localized interactions between PPIn lipids and the wide-variety of PPIn-binding proteins that are essential for the coordinate regulation of cellular metabolism and membrane dynamics.

Molecular dynamics simulations of nucleotide release from the circadian clock protein KaiC reveal atomic-resolution functional insights

The cyanobacterial clock proteins KaiA, KaiB, and KaiC form a powerful system to study the biophysical basis of circadian rhythms, because an in vitro mixture of the three proteins is sufficient to generate a robust ∼24-h rhythm in the phosphorylation of KaiC. The nucleotide-bound states of KaiC critically affect both KaiB binding to the N-terminal domain (CI) and the phosphotransfer reactions that (de)phosphorylate the KaiC C-terminal domain (CII). However, the nucleotide exchange pathways associated with transitions among these states are poorly understood. In this study, we integrate recent advances in molecular dynamics methods to elucidate the structure and energetics of the pathway for Mg·ADP release from the CII domain. We find that nucleotide release is coupled to large-scale conformational changes in the KaiC hexamer. Solvating the nucleotide requires widening the subunit interface leading to the active site, which is linked to extension of the A-loop, a structure implicated in KaiA binding. These results provide a molecular hypothesis for how KaiA acts as a nucleotide exchange factor. In turn, structural parallels between the CI and CII domains suggest a mechanism for allosteric coupling between the domains. We relate our results to structures observed for other hexameric ATPases, which perform diverse functions.