Rational Coupled Dynamics Network Manipulation Rescues Disease-Relevant Mutant Cystic Fibrosis Transmembrane Conductance Regulator

Cracking the Code: Enhancing Molecular Tools for Progress in Nanobiotechnology

Nature continually refines its processes for optimal efficiency, especially within biological systems. This article explores the collaborative efforts of researchers worldwide, aiming to mimic nature's efficiency by developing smarter and more effective nanoscale technologies and biomaterials. Recent advancements highlight progress and prospects in leveraging engineered nucleic acids and proteins for specific tasks, drawing inspiration from natural functions. The focus is developing improved methods for characterizing, understanding, and reprogramming these materials to perform user-defined functions, including personalized therapeutics, targeted drug delivery approaches, engineered scaffolds, and reconfigurable nanodevices. Contributions from academia, government agencies, biotech, and medical settings offer diverse perspectives, promising a comprehensive approach to broad nanobiotechnology objectives. Encompassing topics from mRNA vaccine design to programmable protein-based nanocomputing agents, this work provides insightful perspectives on the trajectory of nanobiotechnology toward a future of enhanced biomimicry and technological innovation.

Allosteric regulation of kinase activity in living cells

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

Allosteric regulation of kinase activity in living cells

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

Comparison of allosteric signaling in DnaK and BiP using mutual information between simulated residue conformations

The heat shock protein 70 kDa (Hsp70) chaperone system serves as a critical component of protein quality control across a wide range of prokaryotic and eukaryotic organisms. Divergent evolution and specialization to particular organelles have produced numerous Hsp70 variants which share similarities in structure and general function, but differ substantially in regulatory aspects, including conformational dynamics and activity modulation by cochaperones. The human Hsp70 variant BiP (also known as GRP78 or HSPA5) is of therapeutic interest in the context of cancer, neurodegenerative diseases, and viral infection, including for treatment of the pandemic virus SARS-CoV-2. Due to the complex conformational rearrangements and high sequential variance within the Hsp70 protein family, it is in many cases poorly understood which amino acid mutations are responsible for biochemical differences between protein variants. In this study, we predicted residues associated with conformational regulation of human BiP and Escherichia coli DnaK. Based on protein structure networks obtained from molecular dynamics simulations, we analyzed the shared information between interaction timelines to highlight residue positions with strong conformational coupling to their environment. Our predictions, which focus on the binding processes of the chaperone's substrate and cochaperones, indicate residues filling potential signaling roles specific to either DnaK or BiP. By combining predictions of individual residues into conformationally coupled chains connecting ligand binding sites, we predict a BiP specific secondary signaling pathway associated with substrate binding. Our study sheds light on mechanistic differences in signaling and regulation between Hsp70 variants, which provide insights relevant to therapeutic applications of these proteins.

Spatial covariance analysis reveals the residue-by-residue thermodynamic contribution of variation to the CFTR fold

Although the impact of genome variation on the thermodynamic properties of function on the protein fold has been studied in vitro, it remains a challenge to assign these relationships across the entire polypeptide sequence in vivo. Using the Gaussian process regression based principle of Spatial CoVariance, we globally assign on a residue-by-residue basis the biological thermodynamic properties that contribute to the functional fold of CFTR in the cell. We demonstrate the existence of a thermodynamically sensitive region of the CFTR fold involving the interface between NBD1 and ICL4 that contributes to its export from endoplasmic reticulum. At the cell surface a new set of residues contribute uniquely to the management of channel function. These results support a general 'quality assurance' view of global protein fold management as an SCV principle describing the differential pre- and post-ER residue interactions contributing to compartmentalization of the energetics of the protein fold for function. Our results set the stage for future analyses of the quality systems managing protein sequence-to-function-to-structure broadly encompassing genome design leading to protein function in complex cellular relationships responsible for diversity and fitness in biology in response to the environment.

Design and engineering of allosteric communications in proteins

Allostery in proteins plays an important role in regulating protein activities and influencing many biological processes such as gene expression, enzyme catalysis, and cell signaling. The process of allostery takes place when a signal detected at a site on a protein is transmitted via a mechanical pathway to a functional site and, thus, influences its activity. The pathway of allosteric communication consists of amino acids that form a network with covalent and non-covalent bonds. By mutating residues in this allosteric network, protein engineers have successfully established novel allosteric pathways to achieve desired properties in the target protein. In this review, we highlight the most recent and state-of-the-art techniques for allosteric communication engineering. We also discuss the challenges that need to be overcome and future directions for engineering protein allostery.

Mapping co-regulatory interactions among ligand-binding sites in ryanodine receptor 1

Ryanodine receptor 1 (RyR1) is an intracellular calcium ion (Ca²⁺ ) release channel required for skeletal muscle contraction. Although cryo-electron microscopy identified binding sites of three coactivators Ca²⁺ , ATP, and caffeine (CFF), the mechanism of co-regulation and synergy of these activators is unknown. Here, we report allosteric connections among the three ligand-binding sites and pore region in (i) Ca²⁺ bound-closed, (ii) ATP/CFF bound-closed, (iii) Ca²⁺ /ATP/CFF bound-closed, and (iv) Ca²⁺ /ATP/CFF bound-open RyR1 states. We identified two dominant networks of interactions that mediate communication between the Ca²⁺ -binding site and pore region in Ca²⁺ bound-closed state, which partially overlapped with the pore communications in ATP/CFF bound-closed RyR1 state. In Ca²⁺ /ATP/CFF bound-closed and -open RyR1 states, co-regulatory interactions were analogous to communications in the Ca²⁺ bound-closed and ATP/CFF bound-closed states. Both ATP- and CFF-binding sites mediate communication between the Ca²⁺ -binding site and the pore region in Ca²⁺ /ATP/CFF bound-open RyR1 structure. We conclude that Ca²⁺ , ATP, and CFF propagate their effects to the pore region through a network of overlapping interactions that mediate allosteric control and molecular synergy in channel regulation.

Comparative Analysis of Structural and Dynamical Features of Ribosome Upon Association With mRNA Reveals Potential Role of Ribosomal Proteins

Ribosomes play a critical role in maintaining cellular proteostasis. The binding of messenger RNA (mRNA) to the ribosome regulates kinetics of protein synthesis. To generate an understanding of the structural, mechanistic, and dynamical features of mRNA recognition in the ribosome, we have analysed mRNA-protein interactions through a structural comparison of the ribosomal complex in the presence and absence of mRNA. To do so, we compared the 3-Dimensional (3D) structures of components of the two assembly structures and analysed their structural differences because of mRNA binding, using elastic network models and structural network-based analysis. We observe that the head region of 30S ribosomal subunit undergoes structural displacement and subunit rearrangement to accommodate incoming mRNA. We find that these changes are observed in proteins that lie far from the mRNA-protein interface, implying allostery. Further, through perturbation response scanning, we show that the proteins S13, S19, and S20 act as universal sensors that are sensitive to changes in the inter protein network, upon binding of 30S complex with mRNA and other initiation factors. Our study highlights the significance of mRNA binding in the ribosome complex and identifies putative allosteric sites corresponding to alterations in structure and/or dynamics, in regions away from mRNA binding sites in the complex. Overall, our work provides fresh insights into mRNA association with the ribosome, highlighting changes in the interactions and dynamics of the ribosome assembly because of the binding.

A topological switch in CFTR modulates channel activity and sensitivity to unfolding

The cystic fibrosis transmembrane conductance regulator (CFTR) anion channel is essential to maintain fluid homeostasis in key organs. Functional impairment of CFTR due to mutations in the cftr gene leads to cystic fibrosis. Here, we show that the first nucleotide-binding domain (NBD1) of CFTR can spontaneously adopt an alternate conformation that departs from the canonical NBD fold previously observed. Crystallography reveals that this conformation involves a topological reorganization of NBD1. Single-molecule fluorescence resonance energy transfer microscopy shows that the equilibrium between the conformations is regulated by adenosine triphosphate binding. However, under destabilizing conditions, such as the disease-causing mutation F508del, this conformational flexibility enables unfolding of the β-subdomain. Our data indicate that, in wild-type CFTR, this conformational transition of NBD1 regulates channel function, but, in the presence of the F508del mutation, it allows domain misfolding and subsequent protein degradation. Our work provides a framework to design conformation-specific therapeutics to prevent noxious transitions.

Nanoscale programming of cellular and physiological phenotypes: inorganic meets organic programming

The advent of protein design in recent years has brought us within reach of developing a "nanoscale programing language," in which molecules serve as operands with their conformational states functioning as logic gates. Combining these operands into a set of operations will result in a functional program, which is executed using nanoscale computing agents (NCAs). These agents would respond to any given input and return the desired output signal. The ability to utilize natural evolutionary processes would allow code to "evolve" in the course of computation, thus enabling radically new algorithmic developments. NCAs will revolutionize the studies of biological systems, enable a deeper understanding of human biology and disease, and facilitate the development of in situ precision therapeutics. Since NCAs can be extended to novel reactions and processes not seen in biological systems, the growth of this field will spark the growth of biotechnological applications with wide-ranging impacts, including fields not typically considered relevant to biology. Unlike traditional approaches in synthetic biology that are based on the rewiring of signaling pathways in cells, NCAs are autonomous vehicles based on single-chain proteins. In this perspective, I will introduce and discuss this new field of biological computing, as well as challenges and the future of the NCA. Addressing these challenges will provide a significant leap in technology for programming living cells.

Engineering an Allosteric Control of Protein Function

Allosteric regulation in proteins is fundamental to many important biological processes. Allostery has been employed to control protein functions by regulating protein activity. Engineered allosteric regulation allows controlling protein activity in subsecond time scale and has a broad range of applications, from dissecting spatiotemporal dynamics in biochemical cascades to applications in biotechnology and medicine. Here, we review the concept of allostery in proteins and various approaches to identify allosteric sites and pathways. We then provide an overview of strategies and tools used in allosteric protein regulation and their utility in biological applications. We highlight various classes of proteins, where regulation is achieved through allostery. Finally, we analyze the current problems, critical challenges, and future prospective in achieving allosteric regulation in proteins.

A Time-Dependent Quantum Approach to Allostery and a Comparison With Light-Harvesting in Photosynthetic Phenomenon

The allosteric effect is one of the most important processes in regulating the function of proteins, and the elucidation of this phenomenon plays a significant role in understanding emergent behaviors in biological regulation. In this process, a perturbation, generated by a ligand in a part of the macromolecule (the allosteric site), moves along this system and reaches a specific (active) site, dozens of Ångströms away, with a great efficiency. The dynamics of this perturbation in the macromolecule can model precisely the allosteric process. In this article, we will be studying the general characteristics of allostery, using a time-dependent quantum approach to obtain rules that apply to this kind of process. Considering the perturbation as a wave that moves within the molecular system, we will characterize the allosteric process with three of the properties of this wave in the active site: (1) t_a, the characteristic time for reaching that site, (2) A_a, the amplitude of the wave in this site, and (3) B_a, its corresponding spectral broadening. These three parameters, together with the process mechanism and the perturbation efficiency in the process, can describe the phenomenon. One of the main purposes of this paper is to link the parameters t_a, A_a, and B_a and the perturbation efficiency to the characteristics of the system. There is another fundamental process for life that has some characteristics similar to allostery: the light-harvesting (LH) process in photosynthesis. Here, as in allostery, two distant macromolecular sites are involved—two sites dozens of Ångströms away. In both processes, it is particularly important that the perturbation is distributed efficiently without dissipating in the infinite degrees of freedom within the macromolecule. The importance of considering quantum effects in the LH process is well documented in literature, and the quantum coherences are experimentally proven by time-dependent spectroscopic techniques. Given the existing similarities between these two processes in macromolecules, in this work, we suggest using Quantum Mechanics (QM) to study allostery.

Mapping allosteric communications within individual proteins

Allostery in proteins influences various biological processes such as regulation of gene transcription and activities of enzymes and cell signaling. Computational approaches for analysis of allosteric coupling provide inexpensive opportunities to predict mutations and to design small-molecule agents to control protein function and cellular activity. We develop a computationally efficient network-based method, Ohm, to identify and characterize allosteric communication networks within proteins. Unlike previously developed simulation-based approaches, Ohm relies solely on the structure of the protein of interest. We use Ohm to map allosteric networks in a dataset composed of 20 proteins experimentally identified to be allosterically regulated. Further, the Ohm allostery prediction for the protein CheY correlates well with NMR CHESCA studies. Our webserver, Ohm.dokhlab.org, automatically determines allosteric network architecture and identifies critical coupled residues within this network.

Residue-Level Allostery Propagates through the Effective Coarse-Grained Hessian

The long-ranged coupling between residues that gives rise to allostery in a protein is built up from short-ranged physical interactions. Computational tools used to predict this coupling and its functional relevance have relied on the application of graph theoretical metrics to residue-level correlations measured from all-atom molecular dynamics simulations. The short-ranged interactions that yield these long-ranged residue-level correlations are quantified by the effective coarse-grained Hessian. Here we compute an effective harmonic coarse-grained Hessian from simulations of a benchmark allosteric protein, IGPS, and demonstrate the improved locality of this graph Laplacian over two other connectivity matrices. Additionally, two centrality metrics are developed that indicate the direct and indirect importance of each residue at producing the covariance between the effector binding pocket and the active site. The residue importance indicated by these two metrics is corroborated by previous mutagenesis experiments and leads to unique functional insights; in contrast to previous computational analyses, our results suggest that fP76-hK181 is the most important contact for conveying direct allosteric paths across the HisF-HisH interface. The connectivity around fD98 is found to be important at affecting allostery through indirect means.

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.

Bridging Genomics to Phenomics at Atomic Resolution through Variation Spatial Profiling

To understand the impact of genome sequence variation (the genotype) responsible for biological diversity and human health (the phenotype) including cystic fibrosis and Alzheimer's disease, we developed a Gaussian-process-based machine learning (ML) approach, variation spatial profiling (VSP). VSP uses a sparse collection of known variants found in the population that perturb the protein fold to define unknown variant function based on the emergent general principle of spatial covariance (SCV). SCV quantitatively captures the role of proximity in genotype-to-phenotype spatial-temporal relationships. Phenotype landscapes generated through SCV provide a platform that can be used to describe the functional properties that drive sequence-to-function-to-structure design of the polypeptide fold at atomic resolution. We provide proof of principle that SCV can enable the use of population-based genomic platforms to define the origins and mechanism of action of genotype-to-phenotype transformations contributing to the health and disease of an individual.

Disease-relevant mutations alter amino acid co-evolution networks in the second nucleotide binding domain of CFTR

Cystic Fibrosis (CF) is an inherited disease caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) ion channel. Mutations in CFTR cause impaired chloride ion transport in the epithelial tissues of patients leading to cardiopulmonary decline and pancreatic insufficiency in the most severely affected patients. CFTR is composed of twelve membrane-spanning domains, two nucleotide-binding domains (NBDs), and a regulatory domain. The most common mutation in CFTR is a deletion of phenylalanine at position 508 (ΔF508) in NBD1. Previous research has primarily concentrated on the structure and dynamics of the NBD1 domain; However numerous pathological mutations have also been found in the lesser-studied NBD2 domain. We have investigated the amino acid co-evolved network of interactions in NBD2, and the changes that occur in that network upon the introduction of CF and CF-related mutations (S1251N(T), S1235R, D1270N, N1303K(T)). Extensive coupling between the α- and β-subdomains were identified with residues in, or near Walker A, Walker B, H-loop and C-loop motifs. Alterations in the predicted residue network varied from moderate for the S1251T perturbation to more severe for N1303T. The S1235R and D1270N networks varied greatly compared to the wildtype, but these CF mutations only affect ion transport preference and do not severely disrupt CFTR function, suggesting dynamic flexibility in the network of interactions in NBD2. Our results also suggest that inappropriate interactions between the β-subdomain and Q-loop could be detrimental. We also identified mutations predicted to stabilize the NBD2 residue network upon introduction of the CF and CF-related mutations, and these predicted mutations are scored as benign by the MUTPRED2 algorithm. Our results suggest the level of disruption of the co-evolution predictions of the amino acid networks in NBD2 does not have a straightforward correlation with the severity of the CF phenotypes observed.

Intrinsically disordered regions regulate the activities of ATP binding cassette transporters

ATP binding cassette (ABC) proteins are a large family of membrane proteins present in all kingdoms of life. These multi-domain proteins are comprised, at minimum, of two membrane-spanning domains (MSD1, MSD2) and two cytosolic nucleotide binding domains (NBD1, NBD2). ATP binding and hydrolysis at the NBDs enables ABC proteins to actively transport solutes across membranes, regulate activities of other proteins, or function as channels. Like most eukaryotic membrane proteins, ABC proteins contain intrinsically disordered regions (IDRs). These conformationally dynamic regions in ABC proteins possess residual structure, are sites of phosphorylation, and mediate protein-protein interactions. Here, we review the role of IDRs in regulating ABC protein activity.

Inter-Active Site Communication Mediated by the Dimer Interface β-Sheet in the Half-the-Sites Enzyme, Thymidylate Synthase

Thymidylate synthase (TS) is a dimeric enzyme conserved in all life forms that exhibits the allosteric feature of half-the-sites activity. Neither the reason for nor the mechanism of this phenomenon is understood. We used a combined nuclear magnetic resonance (NMR) and molecular dynamics approach to study a stable intermediate preceding hydride transfer, which is the rate-limiting and half-the-sites step. In NMR titrations with ligands leading to this intermediate, we measured chemical shifts of the apoenzyme (lig₀), the saturated holoenzyme (lig₂), and the typically elusive singly bound (lig₁) states. Approximately 40 amides showed quartet patterns providing direct NMR evidence of coupling between the active site and probes >30 Å away in the distal subunit. Quartet peak patterns have symmetrical character, indicating reciprocity in communicating the first and second binding events to the distal protomer. Quartets include key catalytic residues and map to the dimer interface β-sheet, which also represents the shortest path between the two active sites. Simulations corroborate the coupling observed in solution in that there is excellent overlap between quartet residues and main-chain atoms having intersubunit cross-correlated motions. Simulations identify five hot spot residues, three of which lie at the kink in the unique β-bulge abutting the active sites on either end of the sheet. Interstrand cross-correlated motions become more organized and pronounced as the enzyme progresses from lig₀ to lig₁ and ultimately lig₂. Coupling in the apparently symmetrical complex has implications for half-the-sites reactivity and potentially resolves the paradox of inequivalent TS active sites despite the vast majority of X-ray structures appearing to be symmetrical.

β-Methylamino-L-alanine substitution of serine in SOD1 suggests a direct role in ALS etiology

Exposure to the environmental toxin β-methylamino-L-alanine (BMAA) is linked to amyotrophic lateral sclerosis (ALS), but its disease-promoting mechanism remains unknown. We propose that incorporation of BMAA into the ALS-linked protein Cu,Zn superoxide dismutase (SOD1) upon translation promotes protein misfolding and aggregation, which has been linked to ALS onset and progression. Using molecular simulation and predictive energetic computation, we demonstrate that substituting any serine with BMAA in SOD1 results in structural destabilization and aberrant dynamics, promoting neurotoxic SOD1 aggregation. We propose that translational incorporation of BMAA into SOD1 is directly responsible for its toxicity in neurodegeneration, and BMAA modification of SOD1 may serve as a biomarker of ALS.

Stabilization of μ-opioid receptor facilitates its cellular translocation and signaling

The G protein-coupled μ-opioid receptor (μ-OR) mediates the majority of analgesia effects for morphine and other pain relievers. Despite extensive studies of its structure and activation mechanisms, the inherently low maturation efficiency of μ-OR represents a major hurdle to understanding its function. Here we computationally designed μ-OR mutants with altered stability to probe the relationship between cell-surface targeting, signal transduction, and agonist efficacy. The stabilizing mutation T315Y enhanced μ-OR trafficking to the plasma membrane and significantly promoted the morphine-mediated inhibition of downstream signaling. In contrast, the destabilizing mutation R165Y led to intracellular retention of μ-OR and reduced the response to morphine stimulation. These findings suggest that μ-OR stability is an important factor in regulating receptor signaling and provide a viable avenue to improve the efficacy of analgesics.

Uncovered Dynamic Coupling Resolves the Ambiguous Mechanism of Phenylalanine Hydroxylase Oxygen Binding

Phenylalanine hydroxylase (PAH) is an iron enzyme catalyzing the oxidation of l-Phe to l-Tyr during phenylalanine catabolism. Dysfunction of PAH leads to the debilitating condition phenylketonuria (PKU), which prompted research into the structure and function of PAH over the last 50 years. Despite intensive study, there is no consensus on the atomistic details of the mechanism of O₂ binding and splitting by wild-type (WT) PAH and how it varies with PKU-inducing mutations, Arg158Gln and Glu280Lys. We studied structures involved in a proposed mechanism for the WT and mutants using extensive mixed quantum-classical molecular dynamics simulations. Simulations reveal a previously unobserved dynamic coupling between the active site and the mutation sites, suggesting how they can affect the catalytic performance of PAH. Furthermore, the effect of the coupling on the PAH structure agrees with and expands our understanding of the experimentally observed differences in activity between the WT and mutants.

Allostery in Its Many Disguises: From Theory to Applications

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

Chokepoints in Mechanical Coupling Associated with Allosteric Proteins: The Pyruvate Kinase Example

Although the critical role of allostery in controlling enzymatic processes is well appreciated, there is a current dearth in our understanding of its underlying mechanisms, including communication between binding sites. One potential key aspect of intersite communication is the mechanical coupling between residues in a protein. Here, we introduce a graph-based computational approach to investigate the mechanical coupling between distant parts of a protein, highlighting effective pathways via which protein motion can transfer energy between sites. In this method, each residue is treated as a node on a weighted, undirected graph, in which the edges are defined by locally correlated motions of those residues and weighted by the strength of the correlation. The method was validated against experimental data on allosteric regulation in the human liver pyruvate kinase as obtained from full-protein alanine-scanning mutagenesis (systematic mutation) studies, as well as computational data on two G-protein-coupled receptors. The method provides semiquantitative information on the regulatory importance of specific structural elements. It is shown that these elements are key for the mechanical coupling between distant parts of the protein by providing effective pathways for energy transfer. It is also shown that, although there are a multitude of energy transfer pathways between distant parts of a protein, these pathways share a few common nodes that represent effective "chokepoints" for the communication.

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

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

Toward Comprehensive Allosteric Control over Protein Activity

Universality of allosteric signaling in proteins, molecular machines, and receptors complemented by the great advantages of prospected allosteric drugs in the highly specific, non-competitive, and modulatory nature of their actions calls for deeper theoretical understanding of allosteric communication. We present a computational model that makes it possible to tackle the problem of modulating the energetics of protein allosteric communication. In the context of the energy landscape paradigm, allosteric signaling is always a result of perturbations, such as ligand binding, mutations, and intermolecular interactions. The calculation of local partition functions in the protein harmonic model with perturbations allows us to evaluate the energetics of allosteric communication at the single-residue level. In this framework, Allosteric Signaling Maps are proposed as a tool to exhaustively describe allosteric communication in the protein, to tune already existing signaling, and to design new elements of regulation for taking the protein activity under allosteric control.

Site-Mutation of Hydrophobic Core Residues Synchronically Poise Super Interleukin 2 for Signaling: Identifying Distant Structural Effects through Affordable Computations

A superkine variant of interleukin-2 with six site mutations away from the binding interface developed from the yeast display technique has been previously characterized as undergoing a distal structure alteration which is responsible for its super-potency and provides an elegant case study with which to get insight about how to utilize allosteric effect to achieve desirable protein functions. By examining the dynamic network and the allosteric pathways related to those mutated residues using various computational approaches, we found that nanosecond time scale all-atom molecular dynamics simulations can identify the dynamic network as efficient as an ensemble algorithm. The differentiated pathways for the six core residues form a dynamic network that outlines the area of structure alteration. The results offer potentials of using affordable computing power to predict allosteric structure of mutants in knowledge-based mutagenesis.

Mapping allosteric linkage to channel gating by extracellular domains in the human epithelial sodium channel

The epithelial sodium channel (ENaC) mediates sodium absorption in lung, kidney, and colon epithelia. Channels in the ENaC/degenerin family possess an extracellular region that senses physicochemical changes in the extracellular milieu and allosterically regulates the channel opening. Proteolytic cleavage activates the ENaC opening, by the removal of specific segments in the finger domains of the α- and γ ENaC-subunits. Cleavage causes perturbations in the extracellular region that propagate to the channel gate. However, it is not known how the channel structure mediates the propagation of activation signals through the extracellular sensing domains. Here, to identify the structure-function determinants that mediate allosteric ENaC activation, we performed MD simulations, thiol modification of residues substituted by cysteine, and voltage-clamp electrophysiology recordings. Our simulations of an ENaC heterotetramer, α₁βα₂γ, in the proteolytically cleaved and uncleaved states revealed structural pathways in the α-subunit that are responsible for ENaC proteolytic activation. To validate these findings, we performed site-directed mutagenesis to introduce cysteine substitutions in the extracellular domains of the α-, β-, and γ ENaC-subunits. Insertion of a cysteine at the α-subunit Glu⁵⁵⁷ site, predicted to stabilize a closed state of ENaC, inhibited ENaC basal activity and retarded the kinetics of proteolytic activation by 2-fold. Our results suggest that the lower palm domain of αENaC is essential for ENaC activation. In conclusion, our integrated computational and experimental approach suggests key structure-function determinants for ENaC proteolytic activation and points toward a mechanistic model for the allosteric communication in the extracellular domains of the ENaC/degenerin family channels.

Molecular dynamics simulation study on the structural instability of the most common cystic fibrosis-associated mutant ΔF508-CFTR

Cystic fibrosis transmembrane conductance regulator (CFTR) is an anion channel that belongs to the ATP binding cassette protein superfamily. Deletion of phenylalanine at position 508 (ΔF508) is the most common CF-associated mutation and is present in nearly 90% of CF patients. Currently, atomistic level studies are insufficient for understanding the mechanism by which the deletion of a single amino acid causes greatly reduced folding as well as trafficking and gating defects. To clarify this mechanism, we first constructed an atomic model of the inward-facing ΔF508-CFTR and performed allatom molecular dynamics (MD) simulations of the protein in a membrane environment. All of the computational methodologies used are based on those developed in our previous study for wild-type CFTR. Two important findings were obtained. First, consistent with several previous computational results, the deletion of F508 causes a disruption of a hydrophobic cluster located at the interface between the nucleotide binding domain 1 (NBD1) and intracellular loop 4 (ICL4). This exerts unfavorable influences on the correlated motion between ICLs and transmembrane domains (TMDs), likely resulting in gating defects. Second, the F508 deletion affected the NBD1-NBD2 interface via allosteric communication originating from the correlated motion between NBDs and ICLs. As a result, several unusual inter-residue interactions are caused at the NBD1-NBD2 interface. In other words, their correct dimerization is impaired. This study provided insight into the atomic-level details of structural and dynamics changes caused by the ΔF508 mutation and thus provides good insight for drug design.

The implications of CFTR structural studies for cystic fibrosis drug development

Development of Cystic Fibrosis Transmembrane conductance Regulator (CFTR) modulators, targeting the root cause of cystic fibrosis (CF), represents a challenge in the era of personalized medicine, as CFTR mutations lead to a variety of phenotypes, which likely require different, specific treatments. CF drug development is also complicated by the need to preserve the right balance between stability and flexibility, required for optimal function of the CFTR protein. In this review, we highlight how structural data can be exploited in this context to understand the molecular mechanisms of disease-associated mutations, to characterize the mechanisms of action of known modulators and to rationalize the search for novel, specific compounds.

Identification of potential allosteric communication pathways between functional sites of the bacterial ribosome by graph and elastic network models

BACKGROUND

Accumulated evidence indicates that bacterial ribosome employs allostery throughout its structure for protein synthesis. The nature of the allosteric communication between remote functional sites remains unclear, but the contact topology and dynamics of residues may play role in transmission of a perturbation to distant sites.

METHODS/RESULTS

We employ two computationally efficient approaches - graph and elastic network modeling to gain insights about the allosteric communication in ribosome. Using graph representation of the structure, we perform k-shortest pathways analysis between peptidyl transferase center-ribosomal tunnel, decoding center-peptidyl transferase center - previously reported functional sites having allosteric communication. Detailed analysis on intact structures points to common and alternative shortest pathways preferred by different states of translation. All shortest pathways capture drug target sites and allosterically important regions. Elastic network model further reveals that residues along all pathways have the ability of quickly establishing pair-wise communication and to help the propagation of a perturbation in long-ranges during functional motions of the complex.

CONCLUSIONS

Contact topology and inherent dynamics of ribosome configure potential communication pathways between functional sites in different translation states. Inter-subunit bridges B2a, B3 and P-tRNA come forward for their high potential in assisting allostery during translation. Especially B3 emerges as a potential druggable site.

GENERAL SIGNIFICANCE

This study indicates that the ribosome topology forms a basis for allosteric communication, which can be disrupted by novel drugs to kill drug-resistant bacteria. Our computationally efficient approach not only overlaps with experimental evidence on allosteric regulation in ribosome but also proposes new druggable sites.

Insights into channel dysfunction from modelling and molecular dynamics simulations

Developments in structural biology mean that the number of different ion channel structures has increased significantly in recent years. Structures of ion channels enable us to rationalize how mutations may lead to channelopathies. However, determining the structures of ion channels is still not trivial, especially as they necessarily exist in many distinct functional states. Therefore, the use of computational modelling can provide complementary information that can refine working hypotheses of both wild type and mutant ion channels. The simplest but still powerful tool is homology modelling. Many structures are available now that can provide suitable templates for many different types of ion channels, allowing a full three-dimensional interpretation of mutational effects. These structural models, and indeed the structures themselves obtained by X-ray crystallography, and more recently cryo-electron microscopy, can be subjected to molecular dynamics simulations, either as a tool to help explore the conformational dynamics in detail or simply as a means to refine the models further. Here we review how these approaches have been used to improve our understanding of how diseases might be linked to specific mutations in ion channel proteins. This article is part of the Special Issue entitled 'Channelopathies.'

Direct Binding of the Corrector VX-809 to Human CFTR NBD1: Evidence of an Allosteric Coupling between the Binding Site and the NBD1:CL4 Interface

Understanding the mechanism of action of modulator compounds for the cystic fibrosis transmembrane conductance regulator (CFTR) is key for the optimization of therapeutics as well as obtaining insights into the molecular mechanisms of CFTR function. We demonstrate the direct binding of VX-809 to the first nucleotide-binding domain (NBD1) of human CFTR. Disruption of the interaction between C-terminal helices and the NBD1 core upon VX-809 binding is observed from chemical shift changes in the NMR spectra of residues in the helices and on the surface of β-strands S3, S9, and S10. Binding to VX-809 leads to a significant negative shift in NBD1 thermal melting temperature (T_m), pointing to direct VX-809 interaction shifting the NBD1 conformational equilibrium. An inter-residue correlation analysis of the chemical shift changes provides evidence of allosteric coupling between the direct binding site and the NBD1:CL4 interface, thus enabling effects on the interface in the absence of direct binding in that location. These NMR binding data and the negative T_m shifts are very similar to those previously reported by us for binding of the dual corrector-potentiator CFFT-001 to NBD1 (Hudson et al., 2012), suggesting that the two compounds may share some aspects of their mechanisms of action. Although previous studies have shown an important role for VX-809 in modulating the conformation of the first membrane spanning domain (Aleksandrov et al., 2012; Ren et al., 2013), this additional mode of VX-809 binding provides insight into conformational dynamics and allostery within CFTR.

Molecular modelling and molecular dynamics of CFTR

The cystic fibrosis transmembrane conductance regulator (CFTR) protein is a member of the ATP-binding cassette (ABC) transporter superfamily that functions as an ATP-gated channel. Considerable progress has been made over the last years in the understanding of the molecular basis of the CFTR functions, as well as dysfunctions causing the common genetic disease cystic fibrosis (CF). This review provides a global overview of the theoretical studies that have been performed so far, especially molecular modelling and molecular dynamics (MD) simulations. A special emphasis is placed on the CFTR-specific evolution of an ABC transporter framework towards a channel function, as well as on the understanding of the effects of disease-causing mutations and their specific modulation. This in silico work should help structure-based drug discovery and design, with a view to develop CFTR-specific pharmacotherapeutic approaches for the treatment of CF in the context of precision medicine.

Computational approaches to investigating allostery

Allosteric regulation plays a key role in many biological processes, such as signal transduction, transcriptional regulation, and many more. It is rooted in fundamental thermodynamic and dynamic properties of macromolecular systems that are still poorly understood and are moreover modulated by the cellular context. Here we review the computational approaches used in the investigation of allosteric processes in protein systems. We outline how the models of allostery have evolved from their initial formulation in the sixties to the current views, which more fully account for the roles of the thermodynamic and dynamic properties of the system. We then describe the major classes of computational approaches employed to elucidate the mechanisms of allostery, the insights they have provided, as well as their limitations. We complement this analysis by highlighting the role of computational approaches in promising practical applications, such as the engineering of regulatory modules and identifying allosteric binding sites.

Pharmacological Chaperones: Design and Development of New Therapeutic Strategies for the Treatment of Conformational Diseases

Errors in protein folding may result in premature clearance of structurally aberrant proteins, or in the accumulation of toxic misfolded species or protein aggregates. These pathological events lead to a large range of conditions known as conformational diseases. Several research groups have presented possible therapeutic solutions for their treatment by developing novel compounds, known as pharmacological chaperones. These cell-permeable molecules selectively provide a molecular scaffold around which misfolded proteins can recover their native folding and, thus, their biological activities. Here, we review therapeutic strategies, clinical potentials, and cost-benefit impacts of several classes of pharmacological chaperones for the treatment of a series of conformational diseases.

Controlling Allosteric Networks in Proteins

Allosteric transition, defined as conformational changes induced by ligand binding, is one of the fundamental properties of proteins. Allostery has been observed and characterized in many proteins, and has been recently utilized to control protein function via regulation of protein activity. Here, we review the physical and evolutionary origin of protein allostery, as well as its importance to protein regulation, drug discovery, and biological processes in living systems. We describe recently developed approaches to identify allosteric pathways, connected sets of pairwise interactions that are responsible for propagation of conformational change from the ligand-binding site to a distal functional site. We then present experimental and computational protein engineering approaches for control of protein function by modulation of allosteric sites. As an example of application of these approaches, we describe a synergistic computational and experimental approach to rescue the cystic-fibrosis-associated protein cystic fibrosis transmembrane conductance regulator, which upon deletion of a single residue misfolds and causes disease. This example demonstrates the power of allosteric manipulation in proteins to both elucidate mechanisms of molecular function and to develop therapeutic strategies that rescue those functions. Allosteric control of proteins provides a tool to shine a light on the complex cascades of cellular processes and facilitate unprecedented interrogation of biological systems.

Thermal unfolding simulations of NBD1 domain variants reveal structural motifs associated with the impaired folding of F508del-CFTR

The deletion of phenylalanine 508 reshapes the conformational space of the NBD1 domain that populates unique intermediate states that provide insights into the molecular events that underlie the impaired folding of F508del-NBD1.

Gain-of-Function Mutation W493R in the Epithelial Sodium Channel Allosterically Reconfigures Intersubunit Coupling

Sodium absorption in epithelial cells is rate-limited by the epithelial sodium channel (ENaC) activity in lung, kidney, and the distal colon. Pathophysiological conditions, such as cystic fibrosis and Liddle syndrome, result from water-electrolyte imbalance partly due to malfunction of ENaC regulation. Because the quaternary structure of ENaC is yet undetermined, the bases of pathologically linked mutations in ENaC subunits α, β, and γ are largely unknown. Here, we present a structural model of heterotetrameric ENaC α1βα2γ that is consistent with previous cross-linking results and site-directed mutagenesis experiments. By using this model, we show that the disease-causing mutation αW493R rewires structural dynamics of the intersubunit interfaces α1β and α2γ. Changes in dynamics can allosterically propagate to the channel gate. We demonstrate that cleavage of the γ-subunit, which is critical for full channel activation, does not mediate activation of ENaC by αW493R. Our molecular dynamics simulations led us to identify a channel-activating electrostatic interaction between α2Arg-493 and γGlu-348 at the α2γ interface. By neutralizing a sodium-binding acidic patch at the α1β interface, we reduced ENaC activation of αW493R by more than 2-fold. By combining homology modeling, molecular dynamics, cysteine cross-linking, and voltage clamp experiments, we propose a dynamics-driven model for the gain-of-function in ENaC by αW493R. Our integrated computational and experimental approach advances our understanding of structure, dynamics, and function of ENaC in its disease-causing state.

Applications of Discrete Molecular Dynamics in biology and medicine

Discrete Molecular Dynamics (DMD) is a physics-based simulation method using discrete energetic potentials rather than traditional continuous potentials, allowing microsecond time scale simulations of biomolecular systems to be performed on personal computers rather than supercomputers or specialized hardware. With the ongoing explosion in processing power even in personal computers, applications of DMD have similarly multiplied. In the past two years, researchers have used DMD to model structures of disease-implicated protein folding intermediates, study assembly of protein complexes, predict protein-protein binding conformations, engineer rescue mutations in disease-causative protein mutants, design a protein conformational switch to control cell signaling, and describe the behavior of polymeric dispersants for environmental cleanup of oil spills, among other innovative applications.

Hallmarks of therapeutic management of the cystic fibrosis functional landscape

The cystic fibrosis (CF) transmembrane conductance regulator (CFTR) protein does not operate in isolation, rather in a dynamic network of interacting components that impact its synthesis, folding, stability, intracellular location and function, referred to herein as the 'CFTR Functional Landscape (CFFL)'. For the prominent F508del mutation, many of these interactors are deeply connected to a protein fold management system, the proteostasis network (PN). However, CF encompasses an additional 2000 CFTR variants distributed along its entire coding sequence (referred to as CFTR2), and each variant contributes a differential liability to PN management of CFTR and to a protein 'social network' (SN) that directs the probability of the (patho)physiologic events that impact ion transport in each cell, tissue and patient in health and disease. Recognition of the importance of the PN and SN in driving the unique patient CFFL leading to disease highlights the importance of precision medicine in therapeutic management of disease progression. We take the view herein that it is not CFTR, rather the PN/SN, and their impact on the CFFL, that are the key physiologic forces driving onset and clinical progression of CF. We posit that a deep understanding of each patients PN/SN gained by merging genomic, proteomic (mass spectrometry (MS)), and high-content microscopy (HCM) technologies in the context of novel network learning algorithms will lead to a paradigm shift in CF clinical management. This should allow for generation of new classes of patient specific PN/SN directed therapeutics for personalized management of the CFFL in the clinic.

REMD Simulations Reveal the Dynamic Profile and Mechanism of Action of Deleterious, Rescuing, and Stabilizing Perturbations to NBD1 from CFTR

Cystic Fibrosis (CF) is a lethal, genetic disease caused by mutations to the CFTR chloride channel. The most common CF causing mutation is the deletion of F508 from the first Nucleotide Binding Domain (F508del-NBD1). This mutation leads to a thermally unstable domain and a misfolded, nonfunctioning CFTR. Replica Exchange MD simulations were used to simulate seven NBD1 constructs including wt and F508del-NBD1 both alone and in the presence of known rescuing mutations as well as F508del-NBD1 in complex with a known small (ligand) stabilizer. Analyzing the resulting trajectories suggests that differences in the biochemical properties of the constructs result from local and coupled differences in their dynamic profiles. A comparative analysis of these profiles as well as of the resulting trajectories reveals how the different perturbations exert their deleterious, rescuing, and stabilizing effects on NBD1. These simulations may therefore be useful for the design and mechanism-of-action analysis of new NBD1 stabilizers.