Identification of key residues for protein conformational transition using elastic network model

Insights into Activation Dynamics and Functional Sites of Inwardly Rectifying Potassium Channel Kir3.2 by an Elastic Network Model Combined with Perturbation Methods

The inwardly rectifying potassium channel Kir3.2, a member of the inward rectifier potassium (Kir) channel family, exerts important biological functions through transporting potassium ions outside of the cell, during which a large-scale synergistic movement occurs among its different domains. Currently, it is not fully understood how the binding of the ligand to the Kir3.2 channel leads to the structural changes and which key residues are responsible for the channel gating and allosteric dynamics. Here, we construct the Gaussian network model (GNM) of the Kir3.2 channel with the secondary structure and covalent interaction information considered (sscGNM), which shows a better performance in reproducing the channel's flexibility compared with the traditional GNM. In addition, the sscANM-based perturbation method is used to simulate the channel's conformational transition caused by the activator PIP2's binding. By applying certain forces to the PIP2 binding pocket, the coarse-grained calculations generate the similar conformational changes to the experimental observation, suggesting that the topology structure as well as PIP2 binding are crucial to the allosteric activation of the Kir3.2 channel. We also utilize the sscGNM-based thermodynamic cycle method developed by us to identify the key residues whose mutations significantly alter the channel's binding free energy with PIP2. We identify not only the residues important for the specific binding but also the ones critical for the allosteric transition coupled with PIP2 binding. This study is helpful for understanding the working mechanism of Kir3.2 channels and can provide important information for related drug design.

Dynamics of binding interactions of TDP-43 and RNA: An equally weighted multiscale elastic network model study

Transactive response DNA binding protein 43 (TDP-43), an alternative-splicing regulator, can specifically bind long UG-rich RNAs, associated with a range of neurodegenerative diseases. Upon binding RNA, TDP-43 undergoes a large conformational change with two RNA recognition motifs (RRMs) connected by a long linker rearranged, strengthening the binding affinity of TDP-43 with RNA. We extend the equally weighted multiscale elastic network model (ewmENM), including its Gaussian network model (ewmGNM) and Anisotropic network model (ewmANM), with the multiscale effect of interactions considered, to the characterization of the dynamics of binding interactions of TDP-43 and RNA. The results reveal upon RNA binding a loss of flexibility occurs to TDP-43's loop3 segments rich in positively charged residues and C-terminal of high flexibility, suggesting their anchoring RNA, induced fit and conformational adjustment roles in recognizing RNA. Additionally, based on movement coupling analyses, it is found that RNA binding strengthens the interactions among intra-RRM β-sheets and between RRMs partially through the linker's mediating role, which stabilizes RNA binding interface, facilitating RNA binding efficiency. In addition, utilizing our proposed thermodynamic cycle method combined with ewmGNM, we identify the key residues for RNA binding whose perturbations induce a large change in binding free energy. We identify not only the residues important for specific binding, but also the ones critical for the conformational rearrangement between RRMs. Furthermore, molecular dynamics simulations are also performed to validate and further interpret the ENM-based results. The study demonstrates a useful avenue to utilize ewmENM to investigate the protein-RNA interaction dynamics characteristics.

Study on functional sites in human multidrug resistance protein 1 (hMRP1)

Human multidrug resistance protein 1 (hMRP1) is an important member of the ATP-binding cassette (ABC) transporter superfamily. It can extrude a variety of anticancer drugs and physiological organic anions across the plasma membrane, which is activated by substrate binding, and is accompanied by large-scale cooperative movements between different domains. Currently, it remains unclear completely about how the specific interactions between hMRP1 and its substrate are and which critical residues are responsible for allosteric signal transduction. To the end, we first construct an inward-facing state of hMRP1 using homology modeling method, and then dock substrate proinflammatory agent leukotriene C4 (LTC4) to hMRP1 pocket. The result manifests LTC4 interacts with two parts of hMRP1 pocket, namely the positively charged pocket (P pocket) and hydrophobic pocket (H pocket), similar to its binding mode with bMRP1 (bovine MRP1). Additionally, we use the Gaussian network model (GNM)-based thermodynamic method proposed by us to identify the key residues whose perturbations markedly alter their binding free energy. Here the conventional GNM is improved with covalent/non-covalent interactions and secondary structure information considered (denoted as sscGNM). In the result, sscGNM improves the flexibility prediction, especially for the nucleotide binding domains with rich kinds of secondary structures. The 46 key residue clusters located in different subdomains are identified which are highly consistent with experimental observations. Furtherly, we explore the long-range cooperation within the transporter. This study is helpful for strengthening the understanding of the work mechanism in ABC exporters and can provide important information to scientists in drug design studies.

Identification of key sites controlling protein functional motions by using elastic network model combined with internal coordinates

The elastic network model (ENM) is an effective method to extract the intrinsic dynamical properties encoded in protein tertiary structures. We have proposed a new ENM-based analysis method to reveal the motion modes directly responsible for a specific protein function, in which an internal coordinate related to the specific function was introduced to construct the internal/Cartesian hybrid coordinate space. In the present work, the function-related internal coordinates combined with a linear perturbation method were applied to identify the key sites controlling specific protein functional motions. The change in the fluctuations of the internal coordinate in response to residue perturbation was calculated in the hybrid coordinate space by using the linear response theory. The residues with the large fluctuation changes were identified to be the key sites that allosterically control the specific protein function. Two proteins, i.e., human DNA polymerase β and the chaperonin from Methanococcus maripaludis, were investigated as case studies, in which several collective and local internal coordinates were applied to identify the functionally key residues of these two studied proteins. The calculation results are consistent with the experimental observations. It is found that different collective internal coordinates lead to similar results, where the predicted functionally key sites are located at similar positions in the protein structure. While for the local internal coordinates, the predicted key sites tend to be situated at the region near to the coordinate-involving residues. Our studies provide a starting point for further exploring other function-related internal coordinates for other interesting proteins.

Interpreting the Dynamics of Binding Interactions of snRNA and U1A Using a Coarse-Grained Model

The binding interactions of small nuclear RNAs (snRNA) and the associated protein factors are critical to the function of spliceosomes in alternatively splicing primary RNA transcripts. Although molecular dynamics simulations are a powerful tool to interpret the mechanism of biological processes, the atomic-level simulations are, however, too expensive and with limited accuracy for the large-size systems, such as snRNA-protein complexes. We extend the coarse-grained Gaussian network model, which models the RNA-protein complexes as a harmonic chain of Cα, P, and O4' atoms, to investigating the impact of the snRNA-binding interaction on the dynamic stability of the human U1A protein, which is a major component of the spliceosomal U1 small nuclear ribonucleoprotein particle. The results reveal that the first and third loops and the C-terminal helix regions of the U1A domain undergo a significant loss of flexibility upon the RNA binding due to the forming of mostly electrostatic and hydrogen bond interactions with RNA 5' stem and loop. By examining the residues whose mutations significantly change the binding free energy between U1A and snRNA, the Gaussian network model-based calculations show that not only the residues at the binding sites that are traditionally considered to play a major role in U1A-RNA association but also those residues that are far away from the RNA-binding interface can participate in the long-range allosteric signal transmission; these calculations are quantitatively consistent with the data observed in the recent snRNA binding experiments. The study demonstrates a useful avenue to utilize the simplified elastic network model to investigate the dynamics characteristics of the biologically important macromolecular interactions.

Energy transport pathway in proteins: Insights from non-equilibrium molecular dynamics with elastic network model

Intra-molecular energy transport between distant functional sites plays important roles in allosterically regulating the biochemical activity of proteins. How to identify the specific intra-molecular signaling pathway from protein tertiary structure remains a challenging problem. In the present work, a non-equilibrium dynamics method based on the elastic network model (ENM) was proposed to simulate the energy propagation process and identify the specific signaling pathways within proteins. In this method, a given residue was perturbed and the propagation of energy was simulated by non-equilibrium dynamics in the normal modes space of ENM. After that, the simulation results were transformed from the normal modes space to the Cartesian coordinate space to identify the intra-protein energy transduction pathways. The proposed method was applied to myosin and the third PDZ domain (PDZ3) of PSD-95 as case studies. For myosin, two signaling pathways were identified, which mediate the energy transductions form the nucleotide binding site to the 50 kDa cleft and the converter subdomain, respectively. For PDZ3, one specific signaling pathway was identified, through which the intra-protein energy was transduced from ligand binding site to the distant opposite side of the protein. It is also found that comparing with the commonly used cross-correlation analysis method, the proposed method can identify the anisotropic energy transduction pathways more effectively.

An efficient perturbation method to predict the functionally key sites of glutamine binding protein

Glutamine-Binding Protein (GlnBP) of Escherichia coli, an important member of the periplasmic binding protein family, is responsible for the first step in the active transport of glutamine across the cytoplasmic membrane. In this work, the functionally key regulation sites of GlnBP were identified by utilizing a perturbation method proposed by our group, in which the residues whose perturbations markedly change the binding free energy between GlnBP and glutamine are considered to be functionally key residues. The results show that besides the substrate binding sites, some other residues distant from the binding pocket, including the ones in the hinge regions between the two domains, the front- and back- door channels and the exposed region, are important for the function of glutamine binding and transport. The predicted results are well consistent with the theoretical and experimental data, which indicates that our method is an effective approach to identify the key residues important for both ligand binding and long-range allosteric signal transmission. This work can provide some insights into the function performance of GlnBP and the physical mechanism of its allosteric regulation.

Evolutionary Conserved Positions Define Protein Conformational Diversity

Conformational diversity of the native state plays a central role in modulating protein function. The selection paradigm sustains that different ligands shift the conformational equilibrium through their binding to highest-affinity conformers. Intramolecular vibrational dynamics associated to each conformation should guarantee conformational transitions, which due to its importance, could possibly be associated with evolutionary conserved traits. Normal mode analysis, based on a coarse-grained model of the protein, can provide the required information to explore these features. Herein, we present a novel procedure to identify key positions sustaining the conformational diversity associated to ligand binding. The method is applied to an adequate refined dataset of 188 paired protein structures in their bound and unbound forms. Firstly, normal modes most involved in the conformational change are selected according to their corresponding overlap with structural distortions introduced by ligand binding. The subspace defined by these modes is used to analyze the effect of simulated point mutations on preserving the conformational diversity of the protein. We find a negative correlation between the effects of mutations on these normal mode subspaces associated to ligand-binding and position-specific evolutionary conservations obtained from multiple sequence-structure alignments. Positions whose mutations are found to alter the most these subspaces are defined as key positions, that is, dynamically important residues that mediate the ligand-binding conformational change. These positions are shown to be evolutionary conserved, mostly buried aliphatic residues localized in regular structural regions of the protein like β-sheets and α-helix.

Domain Motions and Functionally-Key Residues of L-Alanine Dehydrogenase Revealed by an Elastic Network Model

Mycobacterium tuberculosisl-alanine dehydrogenase (l-MtAlaDH) plays an important role in catalyzing l-alanine to ammonia and pyruvate, which has been considered to be a potential target for tuberculosis treatment. In the present work, the functional domain motions encoded in the structure of l-MtAlaDH were investigated by using the Gaussian network model (GNM) and the anisotropy network model (ANM). The slowest modes for the open-apo and closed-holo structures of the enzyme show that the domain motions have a common hinge axis centered in residues Met133 and Met301. Accompanying the conformational transition, both the 1,4-dihydronicotinamide adenine dinucleotide (NAD)-binding domain (NBD) and the substrate-binding domain (SBD) move in a highly coupled way. The first three slowest modes of ANM exhibit the open-closed, rotation and twist motions of l-MtAlaDH, respectively. The calculation of the fast modes reveals the residues responsible for the stability of the protein, and some of them are involved in the interaction with the ligand. Then, the functionally-important residues relevant to the binding of the ligand were identified by using a thermodynamic method. Our computational results are consistent with the experimental data, which will help us to understand the physical mechanism for the function of l-MtAlaDH.

Conformational Motions and Functionally Key Residues for Vitamin B12 Transporter BtuCD-BtuF Revealed by Elastic Network Model with a Function-Related Internal Coordinate

BtuCD-BtuF from Escherichia coli is a binding protein-dependent adenosine triphosphate (ATP)-binding cassette (ABC) transporter system that uses the energy of ATP hydrolysis to transmit vitamin B12 across cellular membranes. Experimental studies have showed that during the transport cycle, the transporter undergoes conformational transitions between the "inward-facing" and "outward-facing" states, which results in the open-closed motions of the cytoplasmic gate of the transport channel. The opening-closing of the channel gate play critical roles for the function of the transporter, which enables the substrate vitamin B12 to be translocated into the cell. In the present work, the extent of opening of the cytoplasmic gate was chosen as a function-related internal coordinate. Then the mean-square fluctuation of the internal coordinate, as well as the cross-correlation between the displacement of the internal coordinate and the movement of each residue in the protein, were calculated based on the normal mode analysis of the elastic network model to analyze the function-related motions encoded in the structure of the system. In addition, the key residues important for the functional motions of the transporter were predicted by using a perturbation method. In order to facilitate the calculations, the internal coordinate was introduced as one of the axes of the coordinate space and the conventional Cartesian coordinate space was transformed into the internal/Cartesian space with linear approximation. All the calculations were carried out in this internal/Cartesian space. Our method can successfully identify the functional motions and key residues for the transporter BtuCD-BtuF, which are well consistent with the experimental observations.

Prediction of protein structural features from sequence data based on Shannon entropy and Kolmogorov complexity

While the genome for a given organism stores the information necessary for the organism to function and flourish it is the proteins that are encoded by the genome that perhaps more than anything else characterize the phenotype for that organism. It is therefore not surprising that one of the many approaches to understanding and predicting protein folding and properties has come from genomics and more specifically from multiple sequence alignments. In this work I explore ways in which data derived from sequence alignment data can be used to investigate in a predictive way three different aspects of protein structure: secondary structures, inter-residue contacts and the dynamics of switching between different states of the protein. In particular the use of Kolmogorov complexity has identified a novel pathway towards achieving these goals.

Study, by use of coarse-grained models, of the functionally crucial residues and allosteric pathway of anesthetic regulation of the Gloeobacter violaceus ligand-gated ion channel

Although pentameric ligand-gated ion channels (pLGICs) have been found to be the targets of general anesthetics, the mechanism of the effects of anesthetics on pLGICs remains elusive. pLGICs from Gloeobacter violaceus (GLIC) can be inhibited by the anesthetic ketamine. X-ray crystallography has shown that the ketamine binding site is distant from the channel gate of the GLIC. It is still not clear how ketamine controls the function of the GLIC by long-range allosteric regulation. In this work, the functionally crucial residues and allosteric pathway of anesthetic regulation of the GLIC were identified by use of a coarse-grained thermodynamic method developed by our group. In our method, the functionally crucial sites were identified as the residues thermodynamically coupled with binding of ketamine. The results from calculation were highly consistent with experimental data. Our study aids understanding of the mechanism of the anesthetic action of ketamine on the GLIC by long-range allosteric modulation.

Prediction of allosteric sites on protein surfaces with an elastic-network-model-based thermodynamic method

Allostery is a rapid and efficient way in many biological processes to regulate protein functions, where binding of an effector at the allosteric site alters the activity and function at a distant active site. Allosteric regulation of protein biological functions provides a promising strategy for novel drug design. However, how to effectively identify the allosteric sites remains one of the major challenges for allosteric drug design. In the present work, a thermodynamic method based on the elastic network model was proposed to predict the allosteric sites on the protein surface. In our method, the thermodynamic coupling between the allosteric and active sites was considered, and then the allosteric sites were identified as those where the binding of an effector molecule induces a large change in the binding free energy of the protein with its ligand. Using the proposed method, two proteins, i.e., the 70 kD heat shock protein (Hsp70) and GluA2 alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA) receptor, were studied and the allosteric sites on the protein surface were successfully identified. The predicted results are consistent with the available experimental data, which indicates that our method is a simple yet effective approach for the identification of allosteric sites on proteins.

Identification of functionally key residues in AMPA receptor with a thermodynamic method

AMPA receptor mediates the fast excitatory synaptic transmission in the central nervous system, and it is activated by the binding of glutamate that results in the opening of the transmembrane ion channel. In the present work, the thermodynamic method developed by our group was improved and then applied to identify the functionally key residues that regulate the glutamate-binding affinity of AMPA receptor. In our method, the key residues are identified as those whose perturbation largely changes the ligand binding free energy of the protein. It is found that besides the ligand binding sites, other residues distant from the binding cleft can also influence the glutamate binding affinity through a long-range allosteric regulation. These allosteric sites include the hinge region of the ligand binding cleft, the dimer interface of the ligand binding domain, the linkers between the ligand binding domain and the transmembrane domain, and the interface between the N-terminal domain and the ligand binding domain. Our calculation results are consistent with the available experimental data. The results are helpful for our understanding of the mechanism of long-range allosteric communication in the AMPA receptor and the mechanism of channel opening triggered by glutamate binding.

Structure and dynamics of molecular networks: a novel paradigm of drug discovery: a comprehensive review

Despite considerable progress in genome- and proteome-based high-throughput screening methods and in rational drug design, the increase in approved drugs in the past decade did not match the increase of drug development costs. Network description and analysis not only give a systems-level understanding of drug action and disease complexity, but can also help to improve the efficiency of drug design. We give a comprehensive assessment of the analytical tools of network topology and dynamics. The state-of-the-art use of chemical similarity, protein structure, protein-protein interaction, signaling, genetic interaction and metabolic networks in the discovery of drug targets is summarized. We propose that network targeting follows two basic strategies. The "central hit strategy" selectively targets central nodes/edges of the flexible networks of infectious agents or cancer cells to kill them. The "network influence strategy" works against other diseases, where an efficient reconfiguration of rigid networks needs to be achieved by targeting the neighbors of central nodes/edges. It is shown how network techniques can help in the identification of single-target, edgetic, multi-target and allo-network drug target candidates. We review the recent boom in network methods helping hit identification, lead selection optimizing drug efficacy, as well as minimizing side-effects and drug toxicity. Successful network-based drug development strategies are shown through the examples of infections, cancer, metabolic diseases, neurodegenerative diseases and aging. Summarizing >1200 references we suggest an optimized protocol of network-aided drug development, and provide a list of systems-level hallmarks of drug quality. Finally, we highlight network-related drug development trends helping to achieve these hallmarks by a cohesive, global approach.

The emergence of protein complexes: quaternary structure, dynamics and allostery. Colworth Medal Lecture

All proteins require physical interactions with other proteins in order to perform their functions. Most of them oligomerize into homomers, and a vast majority of these homomers interact with other proteins, at least part of the time, forming transient or obligate heteromers. In the present paper, we review the structural, biophysical and evolutionary aspects of these protein interactions. We discuss how protein function and stability benefit from oligomerization, as well as evolutionary pathways by which oligomers emerge, mostly from the perspective of homomers. Finally, we emphasize the specificities of heteromeric complexes and their structure and evolution. We also discuss two analytical approaches increasingly being used to study protein structures as well as their interactions. First, we review the use of the biological networks and graph theory for analysis of protein interactions and structure. Secondly, we discuss recent advances in techniques for detecting correlated mutations, with the emphasis on their role in identifying pathways of allosteric communication.