Protein function and allostery: a dynamic relationship

Probing conformational dynamics of DNA binding by CO-sensing transcription factor, CooA

The transcription factor CooA is a CRP/FNR (cAMP receptor protein/ fumarate and nitrate reductase) superfamily protein that uses heme to sense carbon monoxide (CO). Allosteric activation of CooA in response to CO binding is currently described as a series of discrete structural changes, without much consideration for the potential role of protein dynamics in the process of DNA binding. This work uses site-directed spin-label electron paramagnetic resonance spectroscopy (SDSL-EPR) to probe slow timescale (μs-ms) conformational dynamics of CooA with a redox-stable nitroxide spin label, and IR spectroscopy to probe the environment at the CO-bound heme. A series of cysteine substitution variants were created to selectively label CooA in key functional regions, the heme-binding domain, the 4/5-loop, the hinge region, and the DNA binding domain. The EPR spectra of labeled CooA variants are compared across three functional states: Fe(III) "locked off", Fe(II)-CO "on", and Fe(II)-CO bound to DNA. We observe changes in the multicomponent EPR spectra at each location; most notably in the hinge region and DNA binding domain, broadening the description of the CooA allosteric mechanism to include the role of protein dynamics in DNA binding. DNA-dependent changes in IR vibrational frequency and band broadening further suggest that there is conformational heterogeneity in the active WT protein and that DNA binding alters the environment of the heme-bound CO.

Dimeric allostery mechanism of the plant circadian clock photoreceptor ZEITLUPE

In Arabidopsis thaliana, the Light-Oxygen-Voltage (LOV) domain containing protein ZEITLUPE (ZTL) integrates light quality, intensity, and duration into regulation of the circadian clock. Recent structural and biochemical studies of ZTL indicate that the protein diverges from other members of the LOV superfamily in its allosteric mechanism, and that the divergent allosteric mechanism hinges upon conservation of two signaling residues G46 and V48 that alter dynamic motions of a Gln residue implicated in signal transduction in all LOV proteins. Here, we delineate the allosteric mechanism of ZTL via an integrated computational approach that employs atomistic simulations of wild type and allosteric variants of ZTL in the functional dark and light states, together with Markov state and supervised machine learning classification models. This approach has unveiled key factors of the ZTL allosteric mechanisms, and identified specific interactions and residues implicated in functional allosteric changes. The final results reveal atomic level insights into allosteric mechanisms of ZTL function that operate via a non-trivial combination of population-shift and dynamics-driven allosteric pathways.

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.

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.

Binding and structural asymmetry governs ligand sensitivity in a cyclic nucleotide-gated ion channel

HCN channel opening is facilitated by cyclic nucleotides, but what determines the sensitivity of these channels to cAMP or cGMP is unclear. Ng et al. propose that ligand sensitivity depends on negative cooperativity and the asymmetric effects of ligand binding on channel structure and pore opening.

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels open more easily when cAMP or cGMP bind to a domain in the intracellular C-terminus in each of four identical subunits. How sensitivity of the channels to these ligands is determined is not well understood. Here, we apply a mathematical model, which incorporates negative cooperativity, to gating and mutagenesis data available in the literature and combine the results with binding data collected using isothermal titration calorimetry. This model recapitulates the concentration–response data for the effects of cAMP and cGMP on wild-type HCN2 channel opening and, remarkably, predicts the concentration–response data for a subset of mutants with single-point amino acid substitutions in the binding site. Our results suggest that ligand sensitivity is determined by negative cooperativity and asymmetric effects on structure and channel opening, which are tuned by ligand-specific interactions and residues within the binding site.

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.

Evolutionary conservation of the allosteric activation of factor VIIa by tissue factor in lamprey

Essentials Tissue factor (TF) enhances factor VIIa (FVIIa) activity through structural and dynamic changes. We analyzed conservation of TF-activated FVIIa allosteric networks in extant vertebrate lamprey. Lamprey Tf/FVIIa molecular dynamics show conserved Tf-induced structural/dynamic FVIIa changes. Lamprey Tf activation of FVIIa allosteric networks follows molecular pathways similar to human.

SUMMARY

Background Previous studies have provided insight into the molecular basis of human tissue factor (TF) activation of activated factor VII (FVIIa). TF-induced allosteric networks of FVIIa activation have been rationalized through analysis of the dynamic changes and residue connectivities in the human soluble TF (sTF)/FVIIa complex structure during molecular dynamics (MD) simulation. Evolutionary conservation of the molecular mechanisms for TF-induced allosteric FVIIa activation between humans and extant vertebrate jawless fish (lampreys), where blood coagulation emerged more than 500 million years ago, is unknown and of considerable interest. Objective To model the sTf/FVIIa complex from cloned Petromyzon marinus lamprey sequences, and with comparisons to human sTF/FVlla investigate conservation of allosteric mechanisms of FVIIa activity enhancement by soluble TF using MD simulations. Methods Full-length cDNAs of lamprey tf and f7 were cloned and characterized. Comparative models of lamprey sTf/FVIIa complex and free FVIIa were determined based on constructed human sTF/FVIIa complex and free FVIIa models, used in full-atomic MD simulations, and characterized using dynamic network analysis approaches. Results Allosteric paths of correlated motion from Tf contact points in lamprey sTf/FVIIa to the FVIIa active site were determined and quantified, and were found to encompass residue-residue interactions along significantly similar paths compared with human. Conclusions Despite low conservation of residues between lamprey and human proteins, 30% TF and 39% FVII, the structural and protein dynamic effects of TF activation of FVIIa appear conserved and, moreover, present in extant vertebrate proteins from 500 million years ago when TF/FVIIa-initiated extrinsic pathway blood coagulation emerged.

Interfacial water effect on cooperativity and signal communication in Scapharca dimeric hemoglobin

Cooperativity is important in controlling the biological functions of allosteric proteins.

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.

Entrapment and characterization of functional allosteric conformers of hemocyanin in sol–gel matrices

Entrapment of hemocyanin in sol–gel stabilizes conformations scarcely populated in solution, allowing for their structural and functional analysis.

Dancing through Life: Molecular Dynamics Simulations and Network-Centric Modeling of Allosteric Mechanisms in Hsp70 and Hsp110 Chaperone Proteins

Hsp70 and Hsp110 chaperones play an important role in regulating cellular processes that involve protein folding and stabilization, which are essential for the integrity of signaling networks. Although many aspects of allosteric regulatory mechanisms in Hsp70 and Hsp110 chaperones have been extensively studied and significantly advanced in recent experimental studies, the atomistic picture of signal propagation and energetics of dynamics-based communication still remain unresolved. In this work, we have combined molecular dynamics simulations and protein stability analysis of the chaperone structures with the network modeling of residue interaction networks to characterize molecular determinants of allosteric mechanisms. We have shown that allosteric mechanisms of Hsp70 and Hsp110 chaperones may be primarily determined by nucleotide-induced redistribution of local conformational ensembles in the inter-domain regions and the substrate binding domain. Conformational dynamics and energetics of the peptide substrate binding with the Hsp70 structures has been analyzed using free energy calculations, revealing allosteric hotspots that control negative cooperativity between regulatory sites. The results have indicated that cooperative interactions may promote a population-shift mechanism in Hsp70, in which functional residues are organized in a broad and robust allosteric network that can link the nucleotide-binding site and the substrate-binding regions. A smaller allosteric network in Hsp110 structures may elicit an entropy-driven allostery that occurs in the absence of global structural changes. We have found that global mediating residues with high network centrality may be organized in stable local communities that are indispensable for structural stability and efficient allosteric communications. The network-centric analysis of allosteric interactions has also established that centrality of functional residues could correlate with their sensitivity to mutations across diverse chaperone functions. This study reconciles a wide spectrum of structural and functional experiments by demonstrating how integration of molecular simulations and network-centric modeling may explain thermodynamic and mechanistic aspects of allosteric regulation in chaperones.

The role of strong electrostatic interactions at the dimer interface of human glutathione synthetase

The obligate homodimer human glutathione synthetase (hGS) provides an ideal system for exploring the role of protein-protein interactions in the structural stability, activity and allostery of enzymes. The two active sites of hGS, which are 40 Å apart, display allosteric modulation by the substrate γ-glutamylcysteine (γ-GC) during the synthesis of glutathione, a key cellular antioxidant. The two subunits interact at a relatively small dimer interface dominated by electrostatic interactions between S42, R221, and D24. Alanine scans of these sites result in enzymes with decreased activity, altered γ-GC affinity, and decreased thermal stability. Molecular dynamics simulations indicate these mutations disrupt interchain bonding and impact the tertiary structure of hGS. While the ionic hydrogen bonds and salt bridges between S42, R221, and D24 do not mediate allosteric communication in hGS, these interactions have a dramatic impact on the activity and structural stability of the enzyme.

The role of slow and fast protein motions in allosteric interactions

Allostery is fundamentally thermodynamic in nature. Long-range communication in proteins may be mediated not only by changes in the mean conformation with enthalpic contribution but also by changes in dynamic fluctuations with entropic contribution. The important role of protein motions in mediating allosteric interactions has been established by NMR spectroscopy. By using CAP as a model system, we have shown how changes in protein structure and internal dynamics can allosterically regulate protein function and activity. The results indicate that changes in conformational entropy can give rise to binding enhancement, binding inhibition, or have no effect in the expected affinity, depending on the magnitude and sign of enthalpy-entropy compensation. Moreover, allosteric interactions can be regulated by the modulation a low-populated conformation states that serve as on-pathway intermediates for ligand binding. Taken together, the interplay between fast internal motions, which are intimately related to conformational entropy, and slow internal motions, which are related to poorly populated conformational states, can regulate protein activity in a way that cannot be predicted on the basis of the protein's ground-state structure.

Disordered allostery: lessons from glucocorticoid receptor

Allostery is a biological regulation mechanism of significant importance in cell signaling, metabolism, and disease. Although the ensemble basis of allostery has been known for years, only recently has emphasis shifted from interpreting allosteric mechanism in terms of discrete structural pathways to ones that focus on the statistical nature of the signal propagation process, providing a vehicle to unify allostery in structured, dynamic, and disordered systems. In particular, intrinsically disordered (ID) proteins (IDPs), which lack a unique, stable structure, have been directly demonstrated to exhibit allostery in numerous systems, a reality that challenges traditional structure-based models that focus on allosteric pathways. In this chapter, we will discuss the historical context of allostery and focus on studies from human glucocorticoid receptor (GR), a member of the steroid hormone receptor (SHR) family. The numerous translational isoforms of the disordered N-terminal domain of GR consist of coupled thermodynamic domains that contribute to the delicate balance of states in the ensemble and hence in vivo activity. The data are quantitatively interpreted using the ensemble allosteric model (EAM) that considers only the intrinsic and measurable energetics of allosteric systems. It is demonstrated that the EAM provides mechanistic insight into the distribution of states in solution and provides an interpretation for how certain translational isoforms of GR display enhanced and repressed transcriptional activities. The ensemble nature of allostery illuminated from these studies lends credence to the EAM and provides ground rules for allostery in all systems.

Conformational dynamics of nonsynonymous variants at protein interfaces reveals disease association

Recent studies have shown that the protein interface sites between individual monomeric units in biological assemblies are enriched in disease-associated non-synonymous single nucleotide variants (nsSNVs). To elucidate the mechanistic underpinning of this observation, we investigated the conformational dynamic properties of protein interface sites through a site-specific structural dynamic flexibility metric (dfi) for 333 multimeric protein assemblies. dfi measures the dynamic resilience of a single residue to perturbations that occurred in the rest of the protein structure and identifies sites contributing the most to functionally critical dynamics. Analysis of dfi profiles of over a thousand positions harboring variation revealed that amino acid residues at interfaces have lower average dfi (31%) than those present at non-interfaces (50%), which means that protein interfaces have less dynamic flexibility. Interestingly, interface sites with disease-associated nsSNVs have significantly lower average dfi (23%) as compared to those of neutral nsSNVs (42%), which directly relates structural dynamics to functional importance. We found that less conserved interface positions show much lower dfi for disease nsSNVs as compared to neutral nsSNVs. In this case, dfi is better as compared to the accessible surface area metric, which is based on the static protein structure. Overall, our proteome-wide conformational dynamic analysis indicates that certain interface sites play a critical role in functionally related dynamics (i.e., those with low dfi values), therefore mutations at those sites are more likely to be associated with disease.

Insight into the allosteric mechanism of Scapharca dimeric hemoglobin

Allosteric regulation is an essential function of many proteins that control a variety of different processes such as catalysis, signal transduction, and gene regulation. Structural rearrangements have historically been considered the main means of communication between different parts of a protein. Recent studies have highlighted the importance, however, of changes in protein flexibility as an effective way to mediate allosteric communication across a protein. Scapharca dimeric hemoglobin (HbI) is the simplest possible allosteric system, with cooperative ligand binding between two identical subunits. Thermodynamic equilibrium studies of the binding of oxygen to HbI have shown that cooperativity is an entropically driven effect. The change in entropy of the system observed upon ligand binding may arise from changes in the protein, the ligand, or the water of the system. The goal of this study is to determine the contribution of the change in entropy of the protein backbone to HbI cooperative binding. Molecular dynamics simulations and nuclear magnetic resonance relaxation techniques have revealed that the fast internal motions of HbI contribute to the cooperative binding to carbon monoxide in two ways: (1) by contributing favorably to the free energy of the system and (2) by participating in the cooperative mechanism at the HbI subunit interface. The internal dynamics of the weakly cooperative HbI mutant, F97Y, were also investigated with the same methods. The changes in backbone NH dynamics observed for F97Y HbI upon ligand binding are not as large as for the wild type, in agreement with the reduced cooperativity observed for this mutant. The results of this study indicate that interface flexibility and backbone conformational entropy of HbI participate in and are important for the cooperative mechanism of carbon monoxide binding.

Energy landscapes of functional proteins are inherently risky

Evolutionary pressure for protein function leads to unavoidable sampling of conformational states that are at risk of misfolding and aggregation. The resulting tension between functional requirements and the risk of misfolding and/or aggregation in the evolution of proteins is becoming more and more apparent. One outcome of this tension is sensitivity to mutation, in which only subtle changes in sequence that may be functionally advantageous can tip the delicate balance toward protein aggregation. Similarly, increasing the concentration of aggregation-prone species by reducing the ability to control protein levels or compromising protein folding capacity engenders increased risk of aggregation and disease. In this Perspective, we describe examples that epitomize the tension between protein functional energy landscapes and aggregation risk. Each case illustrates how the energy landscapes for the at-risk proteins are sculpted to enable them to perform their functions and how the risks of aggregation are minimized under cellular conditions using a variety of compensatory mechanisms.

Solution NMR structures of Pyrenophora tritici-repentis ToxB and its inactive homolog reveal potential determinants of toxin activity

Pyrenophora tritici-repentis Ptr ToxB (ToxB) is a proteinaceous host-selective toxin produced by Pyrenophora tritici-repentis (P. tritici-repentis), a plant pathogenic fungus that causes the disease tan spot of wheat. One feature that distinguishes ToxB from other host-selective toxins is that it has naturally occurring homologs in non-pathogenic P. tritici-repentis isolates that lack toxic activity. There are no high-resolution structures for any of the ToxB homologs, or for any protein with >30% sequence identity, and therefore what underlies activity remains an open question. Here, we present the NMR structures of ToxB and its inactive homolog Ptr toxb. Both proteins adopt a β-sandwich fold comprising three strands in each half that are bridged together by two disulfide bonds. The inactive toxb, however, shows higher flexibility localized to the sequence-divergent β-sandwich half. The absence of toxic activity is attributed to a more open structure in the vicinity of one disulfide bond, higher flexibility, and residue differences in an exposed loop that likely impacts interaction with putative targets. We propose that activity is regulated by perturbations in a putative active site loop and changes in dynamics distant from the site of activity. Interestingly, the new structures identify AvrPiz-t, a secreted avirulence protein produced by the rice blast fungus, as a structural homolog to ToxB. This homology suggests that fungal proteins involved in either disease susceptibility such as ToxB or resistance such as AvrPiz-t may have a common evolutionary origin.

Harnessing allostery: a novel approach to drug discovery

Allostery is the most direct and efficient way for regulation of biological macromolecule function, ranging from the control of metabolic mechanisms to signal transduction pathways. Allosteric modulators target to allosteric sites, offering distinct advantages compared to orthosteric ligands that target to active sites, such as greater specificity, reduced side effects, and lower toxicity. Allosteric modulators have therefore drawn increasing attention as potential therapeutic drugs in the design and development of new drugs. In recent years, advancements in our understanding of the fundamental principles underlying allostery, coupled with the exploitation of powerful techniques and methods in the field of allostery, provide unprecedented opportunities to discover allosteric proteins, detect and characterize allosteric sites, design and develop novel efficient allosteric drugs, and recapitulate the universal features of allosteric proteins and allosteric modulators. In the present review, we summarize the recent advances in the repertoire of allostery, with a particular focus on the aforementioned allosteric compounds.

Epitope fluctuations in the human papillomavirus are under dynamic allosteric control: a computational evaluation of a new vaccine design strategy

The dynamic properties of the capsid of the human papillomavirus (HPV) type 16 were examined using classical molecular dynamics simulations. By systematically comparing the structural fluctuations of the capsid protein, a strong dynamic allosteric connection between the epitope containing loops and the h4 helix located more than 50 Å away is identified, which was not recognized thus far. Computer simulations show that restricting the structural fluctuations of the h4 helix is key to rigidifying the epitopes, which is thought to be required for eliciting a proper immune response. The allostery identified in the components of the HPV is nonclassical because the mean structure of the epitope carrying loops remains unchanged, but as a result of allosteric effect the structural fluctuations are altered significantly, which in turn changes the biochemical reactivity profile of the epitopes. Exploiting this novel insight, a new vaccine design strategy is proposed wherein a relatively small virus capsid fragment is deposited on a silica nanoparticle in such a way that the fluctuations of the h4 helix are suppressed. The structural and dynamic properties of the epitope carrying loops on this hybrid nanoparticle match the characteristics of epitopes found on the full virus-like particle precisely, suggesting that these nanoparticles may serve as potent, cost-effective, and safe alternatives to traditionally developed vaccines. The structural and dynamic properties of the hybrid nanoparticle are examined in detail to establish the general concepts of the proposed new design.

Allosteric communication in the KIX domain proceeds through dynamic repacking of the hydrophobic core

The KIX domain of the transcriptional coactivator CREB binding protein (CBP) co-operatively mediates interactions between transcription factors. Binding of the transcription factor mixed-lineage leukemia (MLL) induces the formation of a low-populated conformer of KIX that resembles the conformation of the KIX domain in the presence of a second transcription factor molecule. NMR spin relaxation studies have previously shown that allosteric coupling proceeds through a network of hydrophobic core residues that bridge the two binding sites. Here we describe high-resolution NMR solution structures of the binary complex of KIX with MLL and the ternary complex of KIX formed with MLL and phosphorylated kinase inducible domain of CREB (pKID) as a second ligand. We show that binding of pKID to the binary complex of KIX with MLL is accompanied by a defined repacking of the allosteric network in the hydrophobic core of the protein. Rotamer populations derived from methyl group ¹³C chemical shifts reveal a dynamic contribution to the repacking process that is not captured by the structural coordinates and exemplify the dynamic nature of allosteric communication in the KIX domain.

Dynamic fluctuations provide the basis of a conformational switch mechanism in apo cyclic AMP receptor protein

Escherichia coli cyclic AMP Receptor Protein (CRP) undergoes conformational changes with cAMP binding and allosterically promotes CRP to bind specifically to the DNA. In that, the structural and dynamic properties of apo CRP prior to cAMP binding are of interest for the comprehension of the activation mechanism. Here, the dynamics of apo CRP monomer/dimer and holo CRP dimer were studied by Molecular Dynamics (MD) simulations and Gaussian Network Model (GNM). The interplay of the inter-domain hinge with the cAMP and DNA binding domains are pre-disposed in the apo state as a conformational switch in the CRP's allosteric communication mechanism. The hinge at L134-D138 displaying intra- and inter-subunit coupled fluctuations with the cAMP and DNA binding domains leads to the emergence of stronger coupled fluctuations between the two domains and describes an on state. The flexible regions at K52-E58, P154/D155 and I175 maintain the dynamic coupling of the two domains. With a shift in the inter-domain hinge position towards the N terminus, nevertheless, the latter correlations between the domains loosen and become disordered; L134-D138 dynamically interacts only with the cAMP and DNA binding domains of its own subunit, and an off state is assumed. We present a mechanistic view on how the structural dynamic units are hierarchically built for the allosteric functional mechanism; from apo CRP monomer to apo-to-holo CRP dimers.

Phosphorylation of the retinoic acid receptor alpha induces a mechanical allosteric regulation and changes in internal dynamics

Nuclear receptor proteins constitute a superfamily of proteins that function as ligand dependent transcription factors. They are implicated in the transcriptional cascades underlying many physiological phenomena, such as embryogenesis, cell growth and differentiation, and apoptosis, making them one of the major signal transduction paradigms in metazoans. Regulation of these receptors occurs through the binding of hormones, and in the case of the retinoic acid receptor (RAR), through the binding of retinoic acid (RA). In addition to this canonical scenario of RAR activity, recent discoveries have shown that RAR regulation also occurs as a result of phosphorylation. In fact, RA induces non-genomic effects, such as the activation of kinase signaling pathways, resulting in the phosphorylation of several targets including RARs themselves. In the case of RARα, phosphorylation of Ser369 located in loop L9–10 of the ligand-binding domain leads to an increase in the affinity for the protein cyclin H, which is part of the Cdk-activating kinase complex of the general transcription factor TFIIH. The cyclin H binding site in RARα is situated more than 40 Å from the phosphorylated serine. Using molecular dynamics simulations of the unphosphorylated and phosphorylated forms of the receptor RARα, we analyzed the structural implications of receptor phosphorylation, which led to the identification of a structural mechanism for the allosteric coupling between the two remote sites of interest. The results show that phosphorylation leads to a reorganization of a local salt bridge network, which induces changes in helix extension and orientation that affects the cyclin H binding site. This results in changes in conformation and flexibility of the latter. The high conservation of the residues implicated in this signal transduction suggests a mechanism that could be applied to other nuclear receptor proteins.

Allosteric regulation of proteins is critically important in many biological processes. Here we focused on the allosteric pathway of communication within a ligand-regulated transcription factor, the Retinoic Acid Receptor (RAR). Recent experimental studies performed with the RARα subtype have shown that phosphorylation of a residue located at one extremity of an α-helix in RAR, leads to a changes in binding affinity at the other extremity of the same helix for cyclin H, a binding partner that is necessary for gene transcription activation. The purpose of our study was to understand the conformational changes occurring within the receptor upon phosphorylation. Molecular dynamics simulations are well suited for this sort of study. Through this approach, we were able to show that although the overall structure of the phosphorylated RAR shows no distinct difference from the unphosphorylated form, evidence is provided for an allosteric regulation pathway that implicates more subtle changes, such as changes in side chain orientations, which affect the internal dynamics of the receptor.

Ligand-induced dynamic changes in extended PDZ domains from NHERF1

The multi-domain scaffolding protein NHERF1 modulates the assembly and intracellular trafficking of various transmembrane receptors and ion-transport proteins. The two PDZ (postsynaptic density 95/disk large/zonula occluden 1) domains of NHERF1 possess very different ligand-binding capabilities: PDZ1 recognizes a variety of membrane proteins with high affinity, while PDZ2 only binds limited number of target proteins. Here using NMR, we have determined the structural and dynamic mechanisms that differentiate the binding affinities of the two PDZ domains, for the type 1 PDZ-binding motif (QDTRL) in the carboxyl terminus of cystic fibrosis transmembrane regulator. Similar to PDZ2, we have identified a helix-loop-helix subdomain coupled to the canonical PDZ1 domain. The extended PDZ1 domain is highly flexible with correlated backbone motions on fast and slow timescales, while the extended PDZ2 domain is relatively rigid. The malleability of the extended PDZ1 structure facilitates the transmission of conformational changes at the ligand-binding site to the remote helix-loop-helix extension. By contrast, ligand binding has only modest effects on the conformation and dynamics of the extended PDZ2 domain. The study shows that ligand-induced structural and dynamic changes coupled with sequence variation at the putative PDZ binding site dictate ligand selectivity and binding affinity of the two PDZ domains of NHERF1.

Structural dynamics flexibility informs function and evolution at a proteome scale

Protein structures are dynamic entities with a myriad of atomic fluctuations, side-chain rotations, and collective domain movements. Although the importance of these dynamics to proper functioning of proteins is emerging in the studies of many protein families, there is a lack of broad evidence for the critical role of protein dynamics in shaping the biological functions of a substantial fraction of residues for a large number of proteins in the human proteome. Here, we propose a novel dynamic flexibility index (dfi) to quantify the dynamic properties of individual residues in any protein and use it to assess the importance of protein dynamics in 100 human proteins. Our analyses involving functionally critical positions, disease-associated and putatively neutral population variations, and the rate of interspecific substitutions per residue produce concordant patterns at a proteome scale. They establish that the preservation of dynamic properties of residues in a protein structure is critical for maintaining the protein/biological function. Therefore, structural dynamics needs to become a major component of the analysis of protein function and evolution. Such analyses will be facilitated by the dfi, which will also enable the integrative use of structural dynamics with evolutionary conservation in genomic medicine as well as functional genomics investigations.

GTPases IF2 and EF-G bind GDP and the SRL RNA in a mutually exclusive manner

Translational GTPases (trGTPases) are involved in all four stages of protein biosynthesis: initiation, elongation, termination and ribosome recycling. The trGTPases Initiation Factor 2 (IF2) and Elongation Factor G (EF-G) respectively orchestrate initiation complex formation and translocation of the peptidyl-tRNA:mRNA complex through the bacterial ribosome. The ribosome regulates the GTPase cycle and efficiently discriminates between the GDP- and GTP-bound forms of these proteins. Using Isothermal Titration Calorimetry, we have investigated interactions of IF2 and EF-G with the sarcin-ricin loop of the 23S rRNA, a crucial element of the GTPase-associated center of the ribosome. We show that binding of IF2 and EF-G to a 27 nucleotide RNA fragment mimicking the sarcin-ricin loop is mutually exclusive with that of GDP, but not of GTP, providing a mechanism for destabilization of the ribosome-bound GDP forms of translational GTPases.