Allostery and cooperativity revisited

Sparse group selection and analysis of function-related residue for protein-state recognition

Machine learning methods have helped to advance wide range of scientific and technological field in recent years, including computational chemistry. As the chemical systems could become complex with high dimension, feature selection could be critical but challenging to develop reliable machine learning based prediction models, especially for proteins as bio-macromolecules. In this study, we applied sparse group lasso (SGL) method as a general feature selection method to develop classification model for an allosteric protein in different functional states. This results into a much improved model with comparable accuracy (Acc) and only 28 selected features comparing to 289 selected features from a previous study. The Acc achieves 91.50% with 1936 selected feature, which is far higher than that of baseline methods. In addition, grouping protein amino acids into secondary structures provides additional interpretability of the selected features. The selected features are verified as associated with key allosteric residues through comparison with both experimental and computational works about the model protein, and demonstrate the effectiveness and necessity of applying rigorous feature selection and evaluation methods on complex chemical systems.

Prediction of Protein Allosteric Signalling Pathways and Functional Residues Through Paths of Optimised Propensity

Allostery commonly refers to the mechanism that regulates protein activity through the binding of a molecule at a different, usually distal, site from the orthosteric site. The omnipresence of allosteric regulation in nature and its potential for drug design and screening render the study of allostery invaluable. Nevertheless, challenges remain as few computational methods are available to effectively predict allosteric sites, identify signalling pathways involved in allostery, or to aid with the design of suitable molecules targeting such sites. Recently, bond-to-bond propensity analysis has been shown successful at identifying allosteric sites for a large and diverse group of proteins from knowledge of the orthosteric sites and its ligands alone by using network analysis applied to energy-weighted atomistic protein graphs. To address the identification of signalling pathways, we propose here a method to compute and score paths of optimised propensity that link the orthosteric site with the identified allosteric sites, and identifies crucial residues that contribute to those paths. We showcase the approach with three well-studied allosteric proteins: h-Ras, caspase-1, and 3-phosphoinositide-dependent kinase-1 (PDK1). Key residues in both orthosteric and allosteric sites were identified and showed agreement with experimental results, and pivotal signalling residues along the pathway were also revealed, thus providing alternative targets for drug design. By using the computed path scores, we were also able to differentiate the activity of different allosteric modulators.

Kinetic solvent viscosity effects uncover an internal isomerization of the enzyme-substrate complex in Pseudomonas aeruginosa PAO1 NADH:Quinone oxidoreductase

NAD(P)H:quinone oxidoreductases (NQOs) play an essential protective role as antioxidants in the detoxification of quinones in both Prokaryotes and Eukaryotes. NQO from Pseudomonas aeruginosa PAO1 uses FMN to catalyze the two-electron reduction of various quinones with NADH. In this study, steady-state kinetics, kinetic solvent viscosity effects, and rapid reaction kinetics were used to determine which kinetic steps control the overall turnover of the enzyme with benzoquinone or juglone. The rate constant for flavin reduction (k_red) at pH 6.0 was 12.9 ± 0.3 s^-1, and the K_d for NADH was at least an order of magnitude lower than 90 μM. With benzoquinone, the k_cat value was 11.7 ± 0.3 s^-1, consistent with flavin reduction being almost entirely rate-limiting for overall turnover. With juglone, a k_cat value of 10.0 ± 0.5 s^-1 was recorded. The normalized plot of the relative solvent viscosity effects on the k_cat values established that hydride transfer from NADH to the FMN and quinol product release, with a calculated rate constant (k_P-rel) of 52 s^-1, are partially rate-limiting for the overall turnover of NQO. Kinetic solvent viscosity effects with glucose or sucrose revealed a hyperbolic dependence on the k_cat and k_cat/K_m values with benzoquinone or juglone, respectively, consistent with the presence of a solvent-sensitive internal isomerization of the enzyme-substrate complex (ES). The data demonstrate opposing effects of benzoquinone and juglone on the equilibrium of the NQO ES isomerization with glucose or sucrose. Thus, our study demonstrates how quinol substrate properties alter the equilibrium of NQO ES isomerization.

Molecular Dynamics Simulations Establish the Molecular Basis for the Broad Allostery Hotspot Distributions in the Tetracycline Repressor

It is imperative to identify the network of residues essential to the allosteric coupling for the purpose of rationally engineering allostery in proteins. Deep mutational scanning analysis has emerged as a function-centric approach for identifying such allostery hotspots in a comprehensive and unbiased fashion, leading to observations that challenge our understanding of allostery at the molecular level. Specifically, a recent deep mutational scanning study of the tetracycline repressor (TetR) revealed an unexpectedly broad distribution of allostery hotspots throughout the protein structure. Using extensive molecular dynamics simulations (up to 50 μs) and free energy computations, we establish the molecular and energetic basis for the strong anticooperativity between the ligand and DNA binding sites. The computed free energy landscapes in different ligation states illustrate that allostery in TetR is well described by a conformational selection model, in which the apo state samples a broad set of conformations, and specific ones are selectively stabilized by either ligand or DNA binding. By examining a range of structural and dynamic properties of residues at both local and global scales, we observe that various analyses capture different subsets of experimentally identified hotspots, suggesting that these residues modulate allostery in distinct ways. These results motivate the development of a thermodynamic model that qualitatively explains the broad distribution of hotspot residues and their distinct features in molecular dynamics simulations. The multifaceted strategy that we establish here for hotspot evaluations and our insights into their mechanistic contributions are useful for modulating protein allostery in mechanistic and engineering studies.

Cooperative Protein Allosteric Transition Mediated by a Fluctuating Transmission Network

Allosteric communication between distant protein sites represents a key mechanism of biomolecular regulation and signal transduction. Compared to other processes such as protein folding, however, the dynamical evolution of allosteric transitions is still not well understood. As an example of allosteric coupling between distant protein regions, we consider the global open-closed motion of the two domains of T4 lysozyme, which is triggered by local motion in the hinge region. Combining extensive molecular dynamics simulations with a correlation analysis of interresidue contacts, we identify a network of interresidue distances that move in a concerted manner. The cooperative process originates from a cogwheel-like motion of the hydrophobic core in the hinge region, which constitutes an evolutionary conserved and flexible transmission network. Through rigid contacts and the protein backbone, the small local changes of the hydrophobic core are passed on to the distant terminal domains and lead to the emergence of a rare global conformational transition. As in an Ising-type model, the cooperativity of the allosteric transition can be explained via the interaction of local fluctuations.

Prediction of Allosteric Communication Pathways in Proteins

MOTIVATION

Allostery in proteins is an essential phenomenon in biological processes. In this paper, we present a computational model to predict paths of maximum information transfer between active and allosteric sites. In this information theoretic study, we use mutual information as the measure of information transfer, where transition probability of information from one residue to its contacting neighbors is proportional to the magnitude of mutual information between the two residues. Starting from a given residue and using a Hidden Markov Model, we successively determine the neighboring residues that eventually lead to a path of optimum information transfer. The Gaussian approximation of mutual information between residue pairs is adopted. The limits of validity of this approximation are discussed in terms of a nonlinear theory of mutual information and its reduction to the Gaussian form.

RESULTS

Predictions of the model are tested on six widely studied cases, CheY Bacterial Chemotaxis, B-cell Lymphoma extra-large Bcl-xL, Human proline isomerase cyclophilin A (CypA), Dihydrofolate reductase DHFR, HRas GTPase, and Caspase-1. The communication transmission rendering the propagation of local fluctuations from the active sites throughout the structure in multiple paths correlate well with the known experimental data. Distinct paths originating from the active site may likely represent a multi functionality such as involving more than one allosteric site and/or preexistence of some other functional states. Our model is computationally fast and simple, and can give allosteric communication pathways, which are crucial for the understanding and control of protein functionality.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mechanical Allosteric Couplings of Redox-Induced Conformational Changes in Respiratory Complex I

We apply linear response theory to calculate mechanical allosteric couplings in respiratory complex I between the iron sulfur cluster N2, located in the catalytic cavity, and the membrane part of the enzyme, separated from it by more than 50 Å. According to our hypothesis, the redox reaction of ubiquinone in the catalytic cavity of the enzyme generates an unbalanced charge that via repulsion of the charged redox center N2 produces local mechanical stress that transmits into the membrane part of the enzyme where it induces proton pumping. Using coarse-grained simulations of the enzyme, we calculated mechanistic allosteric couplings that reveal the pathways of the mechanical transmission of the stress along the enzyme. The results shed light on the recent experimental studies where a stabilization of the enzyme with an introduced disulfide bridge resulted in the abolishing of proton pumping. Simulation of the disulfide bond action indicates a dramatic change of the mechanistic coupling pathways in line with experimental findings.

Validation of an allosteric binding site of Src kinase identified by unbiased ligand binding simulations

Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often "cryptic" without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a v-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein-ligand systems.

Allostery, and how to define and measure signal transduction

Here we ask: What is productive signaling? How to define it, how to measure it, and most of all, what are the parameters that determine it? Further, what determines the strength of signaling from an upstream to a downstream node in a specific cell? These questions have either not been considered or not entirely resolved. The requirements for the signal to propagate downstream to activate (repress) transcription have not been considered either. Yet, the questions are pivotal to clarify, especially in diseases such as cancer where determination of signal propagation can point to cell proliferation and to emerging drug resistance, and to neurodevelopmental disorders, such as RASopathy, autism, attention-deficit/hyperactivity disorder (ADHD), and cerebral palsy. Here we propose a framework for signal transduction from an upstream to a downstream node addressing these questions. Defining cellular processes, experimentally measuring them, and devising powerful computational AI-powered algorithms that exploit the measurements, are essential for quantitative science.

From Distinct to Differential Conformational Dynamics to Map Allosteric Communication Pathways in Proteins

Initiation of biological processes involving protein-ligand binding, transient protein-protein interactions, or amino acid modifications alters the conformational dynamics of proteins. Accompanying these biological processes are ensuing coupled atomic level conformational changes within the proteins. These conformational changes collectively connect multiple amino acid residues at distal allosteric, binding, and/or active sites. Local changes due to, for example, binding of a regulatory ligand at an allosteric site initiate the allosteric regulation. The allosteric signal propagates throughout the protein structure, causing changes at distal sites, activating, deactivating, or modifying the function of the protein. Hence, dynamical responses within protein structures to stimuli contain critical information on protein function. In this Perspective, we examine the description of allosteric regulation from protein dynamical responses and associated alternative and emerging computational approaches to map allosteric communication pathways between distal sites in proteins at the atomic level.

Bacterial Transcriptional Regulators: A Road Map for Functional, Structural, and Biophysical Characterization

The different niches through which bacteria move during their life cycle require a fast response to the many environmental queues they encounter. The sensing of these stimuli and their correct response is driven primarily by transcriptional regulators. This kind of protein is involved in sensing a wide array of chemical species, a process that ultimately leads to the regulation of gene transcription. The allosteric-coupling mechanism of sensing and regulation is a central aspect of biological systems and has become an important field of research during the last decades. In this review, we summarize the state-of-the-art techniques applied to unravel these complex mechanisms. We introduce a roadmap that may serve for experimental design, depending on the answers we seek and the initial information we have about the system of study. We also provide information on databases containing available structural information on each family of transcriptional regulators. Finally, we discuss the recent results of research about the allosteric mechanisms of sensing and regulation involving many transcriptional regulators of interest, highlighting multipronged strategies and novel experimental techniques. The aim of the experiments discussed here was to provide a better understanding at a molecular level of how bacteria adapt to the different environmental threats they face.

Altered cyclic nucleotide-binding and pore opening in a diseased human HCN4 channel

A growing number of nonsynonymous mutations in the human HCN4 channel gene, the major component of the funny channel of the sinoatrial node, are associated with disease but how they impact channel structure and function, and, thus, how they result in disease, is not clear for any of them. Here, we study the S672R mutation, in the cyclic nucleotide-binding domain of the channel, which has been associated with an inherited bradycardia in an Italian family. This may be the best studied of all known mutations, yet the underlying molecular and atomistic mechanisms remain unclear and controversial. We combine measurements of binding by isothermal titration calorimetry to a naturally occurring tetramer of the HCN4 C-terminal region with a mathematical model to show that weaker binding of cAMP to the mutant channel contributes to a lower level of facilitation of channel opening at submicromolar ligand concentrations but that, in general, facilitation occurs over a range that is similar between the mutant and wild-type because of enhanced opening of the mutant channel when liganded. We also show that the binding affinity for cGMP, which produces the same maximum facilitation of HCN4 opening as cAMP, is weaker in the mutant HCN4 channel but that, for both wild-type and mutant, high-affinity binding of cGMP occurs in a range of concentrations below 1 μM. Thus, binding of cGMP to the HCN4 channel may be relevant normally in vivo and reduced binding of cGMP, as well as cAMP, to the mutant channel may contribute to the reduced resting heart rate observed in the affected family.

The Inherent Coupling of Intrinsically Disordered Regions in the Multidomain Receptor Tyrosine Kinase KIT

RTK KIT regulates a variety of crucial cellular processes via its cytoplasmic domain (CD), which is composed of the tyrosine kinase domain, crowned by the highly flexible domains—the juxtamembrane region, kinase insertion domain, and C-tail, which are key recruitment regions for downstream signalling proteins. To prepare a structural basis for the characterization of the interactions of KIT with its signalling proteins (KIT INTERACTOME), we generated the 3D model of the full-length CD attached to the transmembrane helix. This generic model of KIT in inactive state was studied by molecular dynamics simulation under conditions mimicking the natural environment of KIT. With the accurate atomistic description of the multidomain KIT dynamics, we explained its intrinsic (intra-domain) and extrinsic (inter-domain) disorder and represented the conformational assemble of KIT through free energy landscapes. Strongly coupled movements within each domain and between distant domains of KIT prove the functional interdependence of these regions, described as allosteric regulation, a phenomenon widely observed in many proteins. We suggested that KIT, in its inactive state, encodes all properties of the active protein and its post-transduction events.

Allosterism in the PDZ Family

Dynamic allosterism allows the propagation of signal throughout a protein. The PDZ (PSD-95/Dlg1/ZO-1) family has been named as a classic example of dynamic allostery in small modular domains. While the PDZ family consists of more than 200 domains, previous efforts have primarily focused on a few well-studied PDZ domains, including PTP-BL PDZ2, PSD-95 PDZ3, and Par6 PDZ. Taken together, experimental and computational studies have identified regions of these domains that are dynamically coupled to ligand binding. These regions include the αA helix, the αB lower-loop, and the αC helix. In this review, we summarize the specific residues on the αA helix, the αB lower-loop, and the αC helix of PTP-BL PDZ2, PSD-95 PDZ3, and Par6 PDZ that have been identified as participants in dynamic allostery by either experimental or computational approaches. This review can serve as an index for researchers to look back on the previously identified allostery in the PDZ family. Interestingly, our summary of previous work reveals clear consistencies between the domains. While the PDZ family has a low sequence identity, we show that some of the most consistently identified allosteric residues within PTP-BL PDZ2 and PSD-95 PDZ3 domains are evolutionarily conserved. These residues include A46/A347, V61/V362, and L66/L367 on PTP-BL PDZ2 and PSD-95 PDZ3, respectively. Finally, we expose a need for future work to explore dynamic allostery within (1) PDZ domains with multiple binding partners and (2) multidomain constructs containing a PDZ domain.

Microscopic mechanisms of cooperative communications within single nanocatalysts

Catalysis is an experimental approach to accelerate chemical reactions. It plays a critical role in modern industries. Recent experimental studies uncovered striking observations of cooperative communications for reactions on nanocatalysts. In these experiments, it was shown that the chemical reactions observed at specific active sites might effectively stimulate the same reactions at the neighboring sites. We developed a theoretical model to investigate the microscopic mechanisms of these phenomena. Our idea is that the catalytic communication is the result of the complex dynamics of charged holes. Explicit calculations are able to quantitatively explain all experimental observations, clarifying the molecular origin of cooperative communications. The presented theoretical framework might be utilized for developing efficient catalytic systems with better control over chemical reactions.

Catalysis is a method of accelerating chemical reactions that is critically important for fundamental research as well as for industrial applications. It has been recently discovered that catalytic reactions on metal nanoparticles exhibit cooperative effects. The mechanism of these observations, however, remains not well understood. In this work, we present a theoretical investigation on possible microscopic origin of cooperative communications in nanocatalysts. In our approach, the main role is played by positively charged holes on metal surfaces. A corresponding discrete-state stochastic model for the dynamics of holes is developed and explicitly solved. It is shown that the observed spatial correlation lengths are given by the average distances migrated by the holes before they disappear, while the temporal memory is determined by their lifetimes. Our theoretical approach is able to explain the universality of cooperative communications as well as the effect of external electric fields. Theoretical predictions are in agreement with experimental observations. The proposed theoretical framework quantitatively clarifies some important aspects of the microscopic mechanisms of heterogeneous catalysis.

Structural and (Pseudo-)Enzymatic Properties of Neuroglobin: Its Possible Role in Neuroprotection

Neuroglobin (Ngb), the third member of the globin family, was discovered in human and murine brains in 2000. This monomeric globin is structurally similar to myoglobin (Mb) and hemoglobin (Hb) α and β subunits, but it hosts a bis-histidyl six-coordinated heme-Fe atom. Therefore, the heme-based reactivity of Ngb is modulated by the dissociation of the distal HisE7-heme-Fe bond, which reflects in turn the redox state of the cell. The high Ngb levels (~100–200 μM) present in the retinal ganglion cell layer and in the optic nerve facilitate the O₂ buffer and delivery. In contrast, the very low levels of Ngb (~1 μM) in most tissues and organs support (pseudo-)enzymatic properties including NO/O₂ metabolism, peroxynitrite and free radical scavenging, nitrite, hydroxylamine, hydrogen sulfide reduction, and the nitration of aromatic compounds. Here, structural and (pseudo-)enzymatic properties of Ngb, which are at the root of tissue and organ protection, are reviewed, envisaging a possible role in the protection from neuronal degeneration of the retina and the optic nerve.

Protein Assembly by Design

Proteins are nature's primary building blocks for the construction of sophisticated molecular machines and dynamic materials, ranging from protein complexes such as photosystem II and nitrogenase that drive biogeochemical cycles to cytoskeletal assemblies and muscle fibers for motion. Such natural systems have inspired extensive efforts in the rational design of artificial protein assemblies in the last two decades. As molecular building blocks, proteins are highly complex, in terms of both their three-dimensional structures and chemical compositions. To enable control over the self-assembly of such complex molecules, scientists have devised many creative strategies by combining tools and principles of experimental and computational biophysics, supramolecular chemistry, inorganic chemistry, materials science, and polymer chemistry, among others. Owing to these innovative strategies, what started as a purely structure-building exercise two decades ago has, in short order, led to artificial protein assemblies with unprecedented structures and functions and protein-based materials with unusual properties. Our goal in this review is to give an overview of this exciting and highly interdisciplinary area of research, first outlining the design strategies and tools that have been devised for controlling protein self-assembly, then describing the diverse structures of artificial protein assemblies, and finally highlighting the emergent properties and functions of these assemblies.

The multi-specific human 17 beta-hydroxysteroid dehydrogenase type 7: Non-competitive inhibitors can target different catalyses to facilitate breast cancer treatment

Human 17β-hydroxysteroid dehydrogenase type 7 (17β-HSD7), a special multifunctional enzyme, activates the estrogen estrone while inactivating the potent androgen dihydrotestosterone. Thus, this enzyme has become an ideal target for hormone-dependent breast cancer treatment, as its inhibition leads to estradiol reduction and dihydrotestosterone restoration. However, a particular concern has arisen related to an additional role in cholesterol biosynthesis, as inhibition of the enzyme may lead to undesirable side effects. Our findings demonstrate that the available enzyme inhibitors are non-competitive. Among these, many such as INH81, are specific toward sex-hormone conversion, whereas others represented by 4-bromo-ethynylestradiol, are more specific for zymosterone reduction occurring during cholesterol biosynthesis. The binding of non-competitive inhibitors does not affect the substrate binding on the enzyme. This is the first demonstration of non-competitive inhibitors acting selectively on different catalyses, thereby facilitating inhibitor uses for breast cancer treatment. We aim to quickly communicate the novel results.

The functional role of the hemoglobin-water interface

The interface between hemoglobin (Hb) and its environment, in particular water, is of great physiological relevance. Here, results from in vitro, in vivo, and computational experiments (molecular dynamics simulations) are summarized and put into perspective. One of the main findings from the computations is that the stability of the deoxy, ligand-free T-state (T₀) can be stabilized relative to the deoxy R-state (R₀) only in sufficiently large simulation boxes for the hydrophobic effect to manifest itself. This effect directly influences protein stability and is operative also under physiological conditions. Furthermore, molecular simulations provide a dynamical interpretation of the Perutz model for Hb function. Results from experiments using higher protein concentrations and realistic cellular environments are also discussed. One of the next great challenges for computational studies, which as we show is likely to be taken up in the near future, is to provide a molecular-level understanding of the dynamics of proteins in such crowded environments.

Nicotinic receptors: From protein allostery to computational neuropharmacology

We propose an extension and further development of the Monod-Wyman-Changeux model for allosteric transitions of regulatory proteins to brain communications and specifically to neurotransmitters receptors, with the nicotinic acetylcholine receptor (nAChR) as a model of ligand-gated ion channels. The present development offers an expression of the change of the gating isomerization constant caused by pharmacological ligand binding in terms of its value in the absence of ligands and several "modulation factors", which vary with orthosteric ligand binding (agonists/antagonists), allosteric ligand binding (positive allosteric modulators/negative allosteric modulators) and receptor desensitization. The new - explicit - formulation of such "modulation factors", provides expressions for the pharmacological attributes of potency, efficacy, and selectivity for the modulatory ligands (including endogenous neurotransmitters) in terms of their binding affinity for the active, resting, and desensitized states of the receptor. The current formulation provides ways to design neuroactive compounds with a controlled pharmacological profile, opening the field of computational neuro-pharmacology.

Conformational equilibria in allosteric control of Hsp70 chaperones

Heat-shock proteins of 70 kDa (Hsp70s) are vital for all life and are notably important in protein folding. Hsp70s use ATP binding and hydrolysis at a nucleotide-binding domain (NBD) to control the binding and release of client polypeptides at a substrate-binding domain (SBD); however, the mechanistic basis for this allostery has been elusive. Here, we first characterize biochemical properties of selected domain-interface mutants in bacterial Hsp70 DnaK. We then develop a theoretical model for allosteric equilibria among Hsp70 conformational states to explain the observations: a restraining state, Hsp70_R-ATP, restricts ATP hydrolysis and binds peptides poorly, whereas a stimulating state, Hsp70_S-ATP, hydrolyzes ATP rapidly and has high intrinsic substrate affinity but rapid binding kinetics. We support this model for allosteric regulation with DnaK structures obtained in the postulated stimulating state S with biochemical tests of the S-state interface and with improved peptide-binding-site definition in an R-state structure.

Rigidifying a De Novo Enzyme Increases Activity and Induces a Negative Activation Heat Capacity

Conformational sampling profoundly impacts the overall activity and temperature dependence of enzymes. Peroxidases have emerged as versatile platforms for high-value biocatalysis owing to their broad palette of potential biotransformations. Here, we explore the role of conformational sampling in mediating activity in the de novo peroxidase C45. We demonstrate that 2,2,2-triflouoroethanol (TFE) affects the equilibrium of enzyme conformational states, tending toward a more globally rigid structure. This is correlated with increases in both stability and activity. Notably, these effects are concomitant with the emergence of curvature in the temperature-activity profile, trading off activity gains at ambient temperature with losses at high temperatures. We apply macromolecular rate theory (MMRT) to understand enzyme temperature dependence data. These data point to an increase in protein rigidity associated with a difference in the distribution of protein dynamics between the ground and transition states. We compare the thermodynamics of the de novo enzyme activity to those of a natural peroxidase, horseradish peroxidase. We find that the native enzyme resembles the rigidified de novo enzyme in terms of the thermodynamics of enzyme catalysis and the putative distribution of protein dynamics between the ground and transition states. The addition of TFE apparently causes C45 to behave more like the natural enzyme. Our data suggest robust, generic strategies for improving biocatalytic activity by manipulating protein rigidity; for functional de novo protein catalysts in particular, this can provide more enzyme-like catalysts without further rational engineering, computational redesign, or directed evolution.

Dynamic Cellular Cartography: Mapping the Local Determinants of Oligodendrocyte Transcription Factor 2 (OLIG2) Function in Live Cells Using Massively Parallel Fluorescence Correlation Spectroscopy Integrated with Fluorescence Lifetime Imaging Microscopy (mpFCS/FLIM)

Compartmentalization and integration of molecular processes through diffusion are basic mechanisms through which cells perform biological functions. To characterize these mechanisms in live cells, quantitative and ultrasensitive analytical methods with high spatial and temporal resolution are needed. Here, we present quantitative scanning-free confocal microscopy with single-molecule sensitivity, high temporal resolution (∼10 μs/frame), and fluorescence lifetime imaging capacity, developed by integrating massively parallel fluorescence correlation spectroscopy with fluorescence lifetime imaging microscopy (mpFCS/FLIM); we validate the method, use it to map in live cell location-specific variations in the concentration, diffusion, homodimerization, DNA binding, and local environment of the oligodendrocyte transcription factor 2 fused with the enhanced Green Fluorescent Protein (OLIG2-eGFP), and characterize the effects of an allosteric inhibitor of OLIG2 dimerization on these determinants of OLIG2 function. In particular, we show that cytoplasmic OLIG2-eGFP is largely monomeric and freely diffusing, with the fraction of freely diffusing OLIG2-eGFP molecules being f_D,free^cyt = (0.75 ± 0.10) and the diffusion time τ_D,free^cyt = (0.5 ± 0.3) ms. In contrast, OLIG2-eGFP homodimers are abundant in the cell nucleus, constituting ∼25% of the nuclear pool, some f_D,bound^nuc = (0.65 ± 0.10) of nuclear OLIG2-eGFP is bound to chromatin DNA, whereas freely moving OLIG2-eGFP molecules diffuse at the same rate as those in the cytoplasm, as evident from the lateral diffusion times τ_D,free^nuc = τ_D,free^cyt = (0.5 ± 0.3) ms. OLIG2-eGFP interactions with chromatin DNA, revealed through their influence on the apparent diffusion behavior of OLIG2-eGFP, τ_D,bound^nuc (850 ± 500) ms, are characterized by an apparent dissociation constant K_d,app^OLIG2-DNA = (45 ± 30) nM. The apparent dissociation constant of OLIG2-eGFP homodimers was estimated to be K_d,app^{(OLIG2-eGFP)2} ≈ 560 nM. The allosteric inhibitor of OLIG2 dimerization, compound NSC 50467, neither affects OLIG2-eGFP properties in the cytoplasm nor does it alter the overall cytoplasmic environment. In contrast, it significantly impedes OLIG2-eGFP homodimerization in the cell nucleus, increasing five-fold the apparent dissociation constant, K_{d,app,NSC50467}^{(OLIG2-eGFP)2} ≈ 3 μM, thus reducing homodimer levels to below 7% and effectively abolishing OLIG2-eGFP specific binding to chromatin DNA. The mpFCS/FLIM methodology has a myriad of applications in biomedical research and pharmaceutical industry. For example, it is indispensable for understanding how biological functions emerge through the dynamic integration of location-specific molecular processes and invaluable for drug development, as it allows us to quantitatively characterize the interactions of drugs with drug targets in live cells.

Molecular dynamics simulations of allosteric motions and competitive inhibition of the Zika virus helicase

The 2015 Zika outbreak sparked major global concern and emphasized the reality and dangers still posed by mosquito borne pathogens. While efforts have been made to develop a vaccine and other therapeutics, there is still a great demand for antiviral drugs targeting Zika and other flaviviruses. The non-structural protein 3 (NS3) helicase is a vital component of the viral replication complex, tasked with unwinding the viral dsRNA molecule into single strands. Given this critical function, the Zika virus helicase is a potential therapeutic target and the focus of many ongoing research efforts. Using a combination of drug docking and molecular dynamics simulations, we have identified a list of competitive helicase inhibitors targeting the ATP hydrolysis site and have discovered a potential allosteric site capable of distorting both of the protein's active sites.

Impacts of mutations on dynamic allostery of adenylate kinase

Escherichia coli adenylate kinase (AK) is composed of CORE domain and two branch domains: LID and AMP-binding domain (AMPbd). AK exhibits considerable allostery in a reversible phosphoryl transfer reaction, which is largely attributed to the relative motion of LID and AMPbd with respect to CORE. Such an allosteric conformational change is also evident in the absence of ligands. Recent studies showed that the mutations in branch domains can adjust dynamic allostery and alter the substrate affinity and enzyme activity. In this work, we use all-atom molecular dynamics simulation to study the impacts of mutations in branch domains on AK's dynamic allostery by comparing two double mutants, i.e., LID mutant (Val135Gly, Val142Gly) and AMPbd mutant (Ala37Gly, Ala55Gly), with wild-type. Two mutants undergo considerable conformational fluctuation and exhibit deviation from the initially extended apo state to more compact structures. The LID domain in the LID mutant adjusts its relative position and firmly adheres to CORE. More strikingly, AMPbd mutations affect the relative positions of both the AMPbd domain and remote LID domain. Both domains undergo considerable movement, especially the inherent hinge swing motion of the flexible LID domain. In both mutants, the mutations can enhance the inter-domain interaction. The results about the conformation change of AK in both mutants are in line with the experiment of AK's affinity and activity. As revealed by our findings, the flexibility of branch domains and their inherent motions, especially LID domain, is highly relevant to dynamic allostery in the AK system.

From Selection to Instruction and Back: Competing Conformational Selection and Induced Fit Pathways in Abiotic Hosts

Two limiting cases of molecular recognition, induced fit (IF) and conformational selection (CS), play a central role in allosteric regulation of natural systems. The IF paradigm states that a substrate "instructs" the host to change its shape after complexation, while CS asserts that a guest "selects" the optimal fit from an ensemble of preexisting host conformations. With no studies that quantitatively address the interplay of two limiting pathways in abiotic systems, we herein and for the first time describe the way by which twisted capsule M-1, encompassing two conformers M-1(+) and M-1(-), trap CX₄ (X=Cl, Br) to give CX₄ ⊂M-1(+) and CX₄ ⊂M-1(-), with all four states being in thermal equilibrium. With the assistance of 2D EXSY, we found that CBr₄ would, at its lower concentrations, bind M-1 via a M-1(+)→M-1(-)→CBr₄ ⊂M-1(-) pathway corresponding to conformational selection. For M-1 complexing CCl₄ though, data from 2D EXSY measurements and 1D NMR line-shape analysis suggested that lower CCl₄ concentrations would favor CS while the IF pathway prevailed at higher proportions of the guest. Since CS and IF are not mutually exclusive, we reason that our work sets the stage for characterizing the dynamics of a wide range of already existing hosts to broaden our fundamental understanding of their action. The objective is to master the way in which encapsulation takes place for designing novel and allosteric sequestering agents, catalysts and chemosensors akin to those found in nature.

First-principles prediction of the information processing capacity of a simple genetic circuit

Given the stochastic nature of gene expression, genetically identical cells exposed to the same environmental inputs will produce different outputs. This heterogeneity has been hypothesized to have consequences for how cells are able to survive in changing environments. Recent work has explored the use of information theory as a framework to understand the accuracy with which cells can ascertain the state of their surroundings. Yet the predictive power of these approaches is limited and has not been rigorously tested using precision measurements. To that end, we generate a minimal model for a simple genetic circuit in which all parameter values for the model come from independently published data sets. We then predict the information processing capacity of the genetic circuit for a suite of biophysical parameters such as protein copy number and protein-DNA affinity. We compare these parameter-free predictions with an experimental determination of protein expression distributions and the resulting information processing capacity of E. coli cells. We find that our minimal model captures the scaling of the cell-to-cell variability in the data and the inferred information processing capacity of our simple genetic circuit up to a systematic deviation.

Cooperativity between the orthosteric and allosteric ligand binding sites of RORγt

RORγt is a nuclear receptor associated with several diseases. Various synthetic ligands have been developed that target the canonical orthosteric or a second, allosteric pocket of RORγt. We show that orthosteric and allosteric ligands can simultaneously bind to RORγt and that their potency is positively influenced by the other ligand, a phenomenon called cooperative dual ligand binding. The mechanism behind cooperative binding in proteins is poorly understood, primarily due to the lack of structural data. We solved 12 crystal structures of RORγt, simultaneously bound to various orthosteric and allosteric ligands. In combination with molecular dynamics, we reveal a mechanism responsible for the cooperative binding behavior. Our comprehensive structural studies provide unique insights into how cooperative binding occurs in proteins.

Cooperative ligand binding is an important phenomenon in biological systems where ligand binding influences the binding of another ligand at an alternative site of the protein via an intramolecular network of interactions. The underlying mechanisms behind cooperative binding remain poorly understood, primarily due to the lack of structural data of these ternary complexes. Using time-resolved fluorescence resonance energy transfer (TR-FRET) studies, we show that cooperative ligand binding occurs for RORγt, a nuclear receptor associated with the pathogenesis of autoimmune diseases. To provide the crucial structural insights, we solved 12 crystal structures of RORγt simultaneously bound to various orthosteric and allosteric ligands. The presence of the orthosteric ligand induces a clamping motion of the allosteric pocket via helices 4 to 5. Additional molecular dynamics simulations revealed the unusual mechanism behind this clamping motion, with Ala355 shifting between helix 4 and 5. The orthosteric RORγt agonists regulate the conformation of Ala355, thereby stabilizing the conformation of the allosteric pocket and cooperatively enhancing the affinity of the allosteric inverse agonists.

Roadmap on biology in time varying environments

Biological organisms experience constantly changing environments, from sudden changes in physiology brought about by feeding, to the regular rising and setting of the Sun, to ecological changes over evolutionary timescales. Living organisms have evolved to thrive in this changing world but the general principles by which organisms shape and are shaped by time varying environments remain elusive. Our understanding is particularly poor in the intermediate regime with no separation of timescales, where the environment changes on the same timescale as the physiological or evolutionary response. Experiments to systematically characterize the response to dynamic environments are challenging since such environments are inherently high dimensional. This roadmap deals with the unique role played by time varying environments in biological phenomena across scales, from physiology to evolution, seeking to emphasize the commonalities and the challenges faced in this emerging area of research.

Molecular basis for the allosteric activation mechanism of the heterodimeric imidazole glycerol phosphate synthase complex

Imidazole glycerol phosphate synthase (HisFH) is a heterodimeric bienzyme complex operating at a central branch point of metabolism. HisFH is responsible for the HisH-catalyzed hydrolysis of glutamine to glutamate and ammonia, which is then used for a cyclase reaction by HisF. The HisFH complex is allosterically regulated but the underlying mechanism is not well understood. Here, we elucidate the molecular basis of the long range, allosteric activation of HisFH. We establish that the catalytically active HisFH conformation is only formed when the substrates of both HisH and HisF are bound. We show that in this conformation an oxyanion hole in the HisH active site is established, which rationalizes the observed 4500-fold allosteric activation compared to the inactive conformation. In solution, the inactive and active conformations are in a dynamic equilibrium and the HisFH turnover rates correlate with the population of the active conformation, which is in accordance with the ensemble model of allostery.

The allosteric regulation of the bienzyme complex imidazole glycerol phosphate synthase (HisFH) remains to be elucidated. Here, the authors provide structural insights into the dynamic allosteric mechanism by which ligand binding to the cyclase and glutaminase active sites of HisFH regulate enzyme activation.

Establishing the allosteric mechanism in CRISPR-Cas9

Allostery is a fundamental property of proteins, which regulates biochemical information transfer between spatially distant sites. Here, we report on the critical role of molecular dynamics (MD) simulations in discovering the mechanism of allosteric communication within CRISPR-Cas9, a leading genome editing machinery with enormous promises for medicine and biotechnology. MD revealed how allostery intervenes during at least three steps of the CRISPR-Cas9 function: affecting DNA recognition, mediating the cleavage and interfering with the off-target activity. An allosteric communication that activates concerted DNA cleavages was found to led through the L1/L2 loops, which connect the HNH and RuvC catalytic domains. The identification of these "allosteric transducers" inspired the development of novel variants of the Cas9 protein with improved specificity, opening a new avenue for controlling the CRISPR-Cas9 activity. Discussed studies also highlight the critical role of the recognition lobe in the conformational activation of the catalytic HNH domain. Specifically, the REC3 region was found to modulate the dynamics of HNH by sensing the formation of the RNA:DNA hybrid. The role of REC3 was revealed to be particularly relevant in the presence of DNA mismatches. Indeed, interference of REC3 with the RNA:DNA hybrid containing mismatched pairs at specific positions resulted in locking HNH in an inactive "conformational checkpoint" conformation, thereby hampering off-target cleavages. Overall, MD simulations established the fundamental mechanisms underlying the allosterism of CRISPR-Cas9, aiding engineering strategies to develop new CRISPR-Cas9 variants for improved genome editing.

Dynamic energy correlation analysis of E. coli aspartokinase III and alteration of allosteric regulation by manipulating energy transduction pathways

Conformational change associated with allosteric regulation in a protein is ultimately driven by energy transformation. However, little is known about the latter process. In this work, we combined steered molecular dynamics simulations and sequence conservation analysis to investigate the conformational changes and energy transformation in the allosteric enzyme aspartokinase III (AK III) from Escherichia coli. Correlation analysis of energy change at residue level indicated significant transformation between electrostatic energy and dihedral angle energy during the allosteric regulation. Key amino acid residues located in the corresponding energy transduction pathways were identified by dynamic energy correlation analysis. To verify their functions, residues with a high energy correlation in the pathways were altered and their effects on allosteric regulation of AKIII were determined. This study sheds new insights into energy transformation during allosteric regulation of AK III and proposes a strategy to identify key residues that are involved in intramolecular energy transduction and thus in driving the allosteric process.

Allosteric regulation of histone lysine methyltransferases: from context-specific regulation to selective drugs

Histone lysine methyltransferases (HKMTs) are key regulators of many cellular processes. By definition, HKMTs catalyse the methylation of lysine residues in histone proteins. The enzymatic activities of HKMTs are under precise control, with their allosteric regulation emerging as a prevalent paradigm. We review the molecular mechanisms of allosteric regulation of HKMTs using well-studied histone H3 (K4, K9, K27 and K36) methyltransferases as examples. We discuss the current advances and future potential in targeting allosteric sites of HKMTs for drug development.

On the Emergence of Orientational Order in Folded Proteins with Implications for Allostery

The beautiful structures of single- and multi-domain proteins are clearly ordered in some fashion but cannot be readily classified using group theory methods that are successfully used to describe periodic crystals. For this reason, protein structures are considered to be aperiodic, and may have evolved this way for functional purposes, especially in instances that require a combination of softness and rigidity within the same molecule. By analyzing the solved protein structures, we show that orientational symmetry is broken in the aperiodic arrangement of the secondary structure elements (SSEs), which we deduce by calculating the nematic order parameter, P2. We find that the folded structures are nematic droplets with a broad distribution of P2. We argue that a non-zero value of P2, leads to an arrangement of the SSEs that can resist external forces, which is a requirement for allosteric proteins. Such proteins, which resist mechanical forces in some regions while being flexible in others, transmit signals from one region of the protein to another (action at a distance) in response to binding of ligands (oxygen, ATP, or other small molecules).

Stability and Folding of the Unusually Stable Hemoglobin from Synechocystis is Subtly Optimized and Dependent on the Key Heme Pocket Residues

AIMS

The aim of our study is to understand the biophysical traits that govern the stability and folding of Synechocystis hemoglobin, a unique cyanobacterial globin that displays unusual traits not observed in any of the other globins discovered so far.

BACKGROUND

For the past few decades, classical hemoglobins such as vertebrate hemoglobin and myoglobin have been extensively studied to unravel the stability and folding mechanisms of hemoglobins. However, the expanding wealth of hemoglobins identified in all life forms with novel properties, like heme coordination chemistry and globin fold, have added complexity and challenges to the understanding of hemoglobin stability, which has not been adequately addressed. Here, we explored the unique truncated and hexacoordinate hemoglobin from the freshwater cyanobacterium Synechocystis sp. PCC 6803 known as "Synechocystis hemoglobin (SynHb)". The "three histidines" linkages to heme are novel to this cyanobacterial hemoglobin.

OBJECTIVE

Mutational studies were employed to decipher the residues within the heme pocket that dictate the stability and folding of SynHb.

METHODS

Site-directed mutants of SynHb were generated and analyzed using a repertoire of spectroscopic and calorimetric tools.

RESULTS

The results revealed that the heme was stably associated to the protein under all denaturing conditions with His117 playing the anchoring role. The studies also highlighted the possibility of existence of a "molten globule" like intermediate at acidic pH in this exceptionally thermostable globin. His117 and other key residues in the heme pocket play an indispensable role in imparting significant polypeptide stability.

CONCLUSION

Synechocystis hemoglobin presents an important model system for investigations of protein folding and stability in general. The heme pocket residues influenced the folding and stability of SynHb in a very subtle and specific manner and may have been optimized to make this Hb the most stable known as of date. Other: The knowledge gained hereby about the influence of heme pocket amino acid side chains on stability and expression is currently being utilized to improve the stability of recombinant human Hbs for efficient use as oxygen delivery vehicles.

Exploring the Allosteric Territory of Protein Function

While the pervasiveness of allostery in proteins is commonly accepted, we further show the generic nature of allosteric mechanisms by analyzing here transmembrane ion-channel viroporin 3a and RNA-dependent RNA polymerase (RdRp) from SARS-CoV-2 along with metabolic enzymes isocitrate dehydrogenase 1 (IDH1) and fumarate hydratase (FH) implicated in cancers. Using the previously developed structure-based statistical mechanical model of allostery (SBSMMA), we share our experience in analyzing the allosteric signaling, predicting latent allosteric sites, inducing and tuning targeted allosteric response, and exploring the allosteric effects of mutations. This, yet incomplete list of phenomenology, forms a complex and unique allosteric territory of protein function, which should be thoroughly explored. We propose a generic computational framework, which not only allows one to obtain a comprehensive allosteric control over proteins but also provides an opportunity to approach the fragment-based design of allosteric effectors and drug candidates. The advantages of allosteric drugs over traditional orthosteric compounds, complemented by the emerging role of the allosteric effects of mutations in the expansion of the cancer mutational landscape and in the increased mutability of viral proteins, leave no choice besides further extensive studies of allosteric mechanisms and their biomedical implications.

Biomolecular QM/MM Simulations: What Are Some of the "Burning Issues"?

QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Chemistry, are there still "burning challenges" in QM/MM methods, especially for biomolecular systems? In this short Perspective, we discuss several issues that we believe greatly impact the robustness and quantitative applicability of QM/MM simulations to many, if not all, biomolecules. We highlight these issues with observations and relevant advances from recent studies in our group and others in the field. Despite such limited scope, we hope the discussions are of general interest and will stimulate additional developments that help push the field forward in meaningful directions.

Control of Oxygen Affinity in Mammalian Hemoglobins: Implications for a System Biology Description of the Respiratory Properties of the Red Blood Cell

Hemoglobin and myoglobin have been considered for a long time the paradigmatic model systems for protein function, to the point of being defined the "hydrogen atom[s] of biology". Given this privileged position and the huge amount of quantitative information available on these proteins, the red blood cell might appear as the model system and"hydrogen atom" of system biology. Indeed, since the red cell's main function is O2 transport by hemoglobin, the gap between the protein and the cell may appear quite small. Yet, a surprisingly large amount of detailed biochemical information is required for the modelization of the respiratory properties of the erythrocyte. This problem is compounded if modelization aims at uncovering or explaining evolutionarily selected functional properties of hemoglobin. The foremost difficulty lies in the fact that hemoglobins having different intrinsic properties and relatively ancient evolutionary divergence may behave similarly in the complex milieu of blood, whereas very similar hemoglobins sharing a substantial sequence similarity may present important functional differences because of the mutation of a few key residues. Thus, the functional properties of hemoglobin and blood may reflect more closely the recent environmental challenges than the remote evolutionary history of the animal. We summarize in this review the case of hemoglobins from mammals, in an attempt to provide a reasoned summary of their complexity that, we hope, may be of help to scientists interested in the quantitative exploration of the evolutionary physiology of respiration. Indeed the basis of a meaningful modelization of the red cell requires a large amount of information collected in painstaking and often forgotten studies of the biochemical properties of hemoglobin carried out over more than a century.

Dynamical nonequilibrium molecular dynamics reveals the structural basis for allostery and signal propagation in biomolecular systems

ABSTRACT

A dynamical approach to nonequilibrium molecular dynamics (D-NEMD), proposed in the 1970s by Ciccotti et al., is undergoing a renaissance and is having increasing impact in the study of biological macromolecules. This D-NEMD approach, combining MD simulations in stationary (in particular, equilibrium) and nonequilibrium conditions, allows for the determination of the time-dependent structural response of a system using the Kubo-Onsager relation. Besides providing a detailed picture of the system's dynamic structural response to an external perturbation, this approach also has the advantage that the statistical significance of the response can be assessed. The D-NEMD approach has been used recently to identify a general mechanism of inter-domain signal propagation in nicotinic acetylcholine receptors, and allosteric effects in β -lactamase enzymes, for example. It complements equilibrium MD and is a very promising approach to identifying and analysing allosteric effects. Here, we review the D-NEMD approach and its application to biomolecular systems, including transporters, receptors, and enzymes.

Investigating Conformational Dynamics and Allostery in the p53 DNA-Binding Domain Using Molecular Simulations

The p53 tumor suppressor is a multifaceted context-dependent protein, which is involved in multiple cellular pathways, with the ability to either keep the cells alive or to kill them through mechanisms such as apoptosis. To complicate this picture, cancer cells that express mutant p53 becomes addicted to the mutant activity, so that the mutant variant features a myriad of gain-of-function activities, opening different venues for therapy. This makes essential to think outside the box and apply new approaches to the study of p53 structure-(mis)function relationship to find new critical components of its pathway or to understand how known parts are interconnected, compete, or cooperate. In this context, I will here illustrate how to integrate different computational methods to the identification of possible allosteric effects transmitted from the DNA binding interface of p53 to regions for cofactor recruitment. The protocol can be extended to any other cases of study. Indeed, it does not necessarily apply only to the study of DNA-induced effects, but more broadly to the investigation of long-range effects induced by a biological partner that binds to a biomolecule of interest.

Surveying the Side-Chain Network Approach to Protein Structure and Dynamics: The SARS-CoV-2 Spike Protein as an Illustrative Case

Network theory-based approaches provide valuable insights into the variations in global structural connectivity between different dynamical states of proteins. Our objective is to review network-based analyses to elucidate such variations, especially in the context of subtle conformational changes. We present technical details of the construction and analyses of protein structure networks, encompassing both the non-covalent connectivity and dynamics. We examine the selection of optimal criteria for connectivity based on the physical concept of percolation. We highlight the advantages of using side-chain-based network metrics in contrast to backbone measurements. As an illustrative example, we apply the described network approach to investigate the global conformational changes between the closed and partially open states of the SARS-CoV-2 spike protein. These conformational changes in the spike protein is crucial for coronavirus entry and fusion into human cells. Our analysis reveals global structural reorientations between the two states of the spike protein despite small changes between the two states at the backbone level. We also observe some differences at strategic locations in the structures, correlating with their functions, asserting the advantages of the side-chain network analysis. Finally, we present a view of allostery as a subtle synergistic-global change between the ligand and the receptor, the incorporation of which would enhance drug design strategies.

A hybrid approach reveals the allosteric regulation of GTP cyclohydrolase I

Guanosine triphosphate (GTP) cyclohydrolase I (GCH1) catalyzes the conversion of GTP to dihydroneopterin triphosphate (H2NTP), the initiating step in the biosynthesis of tetrahydrobiopterin (BH4). Besides other roles, BH4 functions as cofactor in neurotransmitter biosynthesis. The BH4 biosynthetic pathway and GCH1 have been identified as promising targets to treat pain disorders in patients. The function of mammalian GCH1s is regulated by a metabolic sensing mechanism involving a regulator protein, GCH1 feedback regulatory protein (GFRP). GFRP binds to GCH1 to form inhibited or activated complexes dependent on availability of cofactor ligands, BH4 and phenylalanine, respectively. We determined high-resolution structures of human GCH1-GFRP complexes by cryoelectron microscopy (cryo-EM). Cryo-EM revealed structural flexibility of specific and relevant surface lining loops, which previously was not detected by X-ray crystallography due to crystal packing effects. Further, we studied allosteric regulation of isolated GCH1 by X-ray crystallography. Using the combined structural information, we are able to obtain a comprehensive picture of the mechanism of allosteric regulation. Local rearrangements in the allosteric pocket upon BH4 binding result in drastic changes in the quaternary structure of the enzyme, leading to a more compact, tense form of the inhibited protein, and translocate to the active site, leading to an open, more flexible structure of its surroundings. Inhibition of the enzymatic activity is not a result of hindrance of substrate binding, but rather a consequence of accelerated substrate binding kinetics as shown by saturation transfer difference NMR (STD-NMR) and site-directed mutagenesis. We propose a dissociation rate controlled mechanism of allosteric, noncompetitive inhibition.

Transient association between proteins elicits alteration of dynamics at sites far away from interfaces

Proteins are known to undergo structural changes upon binding to partner proteins. However, the prevalence, extent, location, and function of change in protein dynamics due to transient protein-protein interactions is not well documented. Here, we have analyzed a dataset of 58 protein-protein complexes of known three-dimensional structure and structures of their corresponding unbound forms to evaluate dynamics changes induced by binding. Fifty-five percent of cases showed significant dynamics change away from the interfaces. This change is not always accompanied by an observed structural change. Binding of protein partner is found to alter inter-residue communication within the tertiary structure in about 90% of cases. Also, residue motions accessible to proteins in unbound form were not always maintained in the bound form. Further analyses revealed functional roles for the distant site where dynamics change was observed. Overall, the results presented here strongly suggest that alteration of protein dynamics due to binding of a partner protein commonly occurs.

Protonation-Induced Dynamic Allostery in PDZ Domain: Evidence of Perturbation-Independent Universal Response Network

Dynamic allostery is a relatively new paradigm where certain external perturbations may lead to modulation of conformational dynamics at a distant part of a protein without significant changes in the overall structure. While most well-characterized examples of dynamic allostery involve binding with other entities like small molecules, peptides, or nucleic acids, in this work we demonstrate that chemical modifications like protonation may lead to significant dynamical allosteric response in a PDZ domain protein. Tuning the protonation states of two histidine residues (H317 and H372), we identify the allosteric pathways responsible for the dynamic response. Interestingly, the same set of residues that constitute the allosteric response network upon ligand binding seem to be responsible for protonation-induced dynamic allostery. Thus, we propose the existence of an inherent universal response network in signaling proteins, where the same set of residues can respond to varying types of external perturbations in terms of rearrangement of hydrogen-bonded network and redistribution of electrostatic interaction energies.

Insight to Functional Conformation and Noncovalent Interactions of Protein-Protein Assembly Using MALDI Mass Spectrometry

Noncovalent interactions are the keys to the structural organization of biomolecule e.g., proteins, glycans, lipids in the process of molecular recognition processes e.g., enzyme-substrate, antigen-antibody. Protein interactions lead to conformational changes, which dictate the functionality of that protein-protein complex. Besides biophysics techniques, noncovalent interaction and conformational dynamics, can be studied via mass spectrometry (MS), which represents a powerful tool, due to its low sample consumption, high sensitivity, and label-free sample. In this review, the focus will be placed on Matrix-Assisted Laser Desorption Ionization Mass Spectrometry (MALDI-MS) and its role in the analysis of protein-protein noncovalent assemblies exploring the relationship within noncovalent interaction, conformation, and biological function.