Inactive-to-Active Transition of Human Thymidine Kinase 1 Revealed by Molecular Dynamics Simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Makurat S
PROVIDER: S-EPMC8757434 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Hemoglobin is one of the most widely studied proteins genetically, biochemically, and structurally. It is an oxygen carrying tetrameric protein that imparts the characteristic red color to blood. Each chain of hemoglobin harbors a heme group embedded in a hydrophobic pocket. Several studies have investigated structural variations present in mammalian hemoglobin and their functional implications. However, camel hemoglobin has not been thoroughly explored, especially from a structural perspective. Importantly, very little is known about how the heme group interacts with hemoglobin under varying conditions of osmolarity and temperature. Several experimental studies have indicated that the tense (T) state is more stable than the relaxed (R) state of hemoglobin under normal physiological conditions. Despite the fact that R state is less stable than the T state, no extensive structural dynamics studies have been performed to investigate global quaternary transitions of R state hemoglobin under normal physiological conditions. To evaluate this, several 500 ns all-atom molecular dynamics simulations were performed to get a deeper understanding of how camel hemoglobin behaves under stress, which it is normally exposed to, when compared to human hemoglobin. Notably, camel hemoglobin was more stable under physiological stress when compared to human hemoglobin. Additionally, when compared to camel hemoglobin, cofactor-binding regions of hemoglobin also exhibited more fluctuations in human hemoglobin under the conditions studied. Several differences were observed between the residues of camel and human hemoglobin that interacted with heme. Importantly, distal residues His58 of α hemoglobin and His63 of β hemoglobin formed more sustained interactions, especially at higher temperatures, in camel hemoglobin. These residues are important for oxygen binding to hemoglobin. Thus, this work provides insights into how camel and human hemoglobin differ in their interactions under stress.
| S-EPMC8741986 | biostudies-literature
Project description:To investigate the dynamics of the orexin 2 receptor, which is a class A G protein-coupled receptor, we recently performed several microsecond-scale molecular dynamics simulations of the wild-type protein, of a mutant that stabilizes the inactive state, and of constitutively active mutants of the class A G protein-coupled receptors. Herein, we review the results of these molecular dynamics simulations of the orexin 2 receptor. In these simulations, characteristic conformational changes were observed in the V3096.40Y mutant. The conformational changes were related to the outward movement of the transmembrane helix 6 and the inward movement of the transmembrane helix 7, which are common structural changes in the activation of G protein-coupled receptors. The index for the quantitative evaluation of the active and inactive states of class A G protein-coupled receptors and the mechanism of the inward movement of the transmembrane helix 7 were examined. In this review, we also discuss the activation mechanism by comparing the structures obtained from the molecular dynamics simulations with the structure of the active state of the orexin 2 receptor clarified by cryo-electron microscopy in the recent years.
| S-EPMC8921386 | biostudies-literature
Project description:Small globular proteins and peptides commonly exhibit two-state folding kinetics in which the rate limiting step of folding is the surmounting of a single free energy barrier at the transition state (TS) separating the folded and the unfolded states. An intriguing question is whether the polypeptide chain reaches, and leaves, the TS by completely random fluctuations, or whether there is a directed, stepwise process. Here, the folding TS of a 15-residue β-hairpin peptide, Peptide 1, is characterized using independent 2.5 μs-long unbiased atomistic molecular dynamics (MD) simulations (a total of 15 μs). The trajectories were started from fully unfolded structures. Multiple (spontaneous) folding events to the NMR-derived conformation are observed, allowing both structural and dynamical characterization of the folding TS. A common loop-like topology is observed in all the TS structures with native end-to-end and turn contacts, while the central segments of the strands are not in contact. Non-native sidechain contacts are present in the TS between the only tryptophan (W11) and the turn region (P7-G9). Prior to the TS the turn is found to be already locked by the W11 sidechain, while the ends are apart. Once the ends have also come into contact, the TS is reached. Finally, along the reactive folding paths the cooperative loss of the W11 non-native contacts and the formation of the central inter-strand native contacts lead to the peptide rapidly proceeding from the TS to the native state. The present results indicate a directed stepwise process to folding the peptide.
| S-EPMC3169518 | biostudies-literature
Project description:To compare ordered water positions from experiment with those from molecular dynamics (MD) simulations, a number of MD models of water structure in crystalline endoglucanase were calculated. The starting MD model was derived from a joint X-ray and neutron diffraction crystal structure, enabling the use of experimentally assigned protonation states. Simulations were performed in the crystalline state, using a periodic 2 × 2 × 2 supercell with explicit solvent. Water X-ray and neutron scattering density maps were computed from MD trajectories using standard macromolecular crystallography methods. In one set of simulations, harmonic restraints were applied to bias the protein structure toward the crystal structure. For these simulations, the recall of crystallographic waters using strong peaks in the MD water electron density was very good, and there also was substantial visual agreement between the boomerang-like wings of the neutron scattering density and the crystalline water hydrogen positions. An unrestrained simulation also was performed. For this simulation, the recall of crystallographic waters was much lower. For both restrained and unrestrained simulations, the strongest water density peaks were associated with crystallographic waters. The results demonstrate that it is now possible to recover crystallographic water structure using restrained MD simulations but that it is not yet reasonable to expect unrestrained MD simulations to do the same. Further development and generalization of MD water models for force-field development, macromolecular crystallography, and medicinal chemistry applications is now warranted. In particular, the combination of room-temperature crystallography, neutron diffraction, and crystalline MD simulations promises to substantially advance modeling of biomolecular solvation.
| S-EPMC6438627 | biostudies-literature
Project description:Aquaporins are protein channels located across the cell membrane with the role of conducting water or other small sugar alcohol molecules (aquaglyceroporins). The high-resolution X-ray structure of the human aquaporin 5 (HsAQP5) shows that HsAQP5, as all the other known aquaporins, exhibits tetrameric structure. By means of molecular dynamics simulations we analyzed the role of spontaneous fluctuations on the structural behavior of the human AQP5. We found that different conformations within the tetramer lead to a distribution of monomeric channel structures, which can be characterized as open or closed. The switch between the two states of a channel is a tap-like mechanism at the cytoplasmic end which regulates the water passage through the pore. The channel is closed by a translation of the His67 residue inside the pore. Moreover, water permeation rate calculations revealed that the selectivity filter, located at the other end of the channel, regulates the flow rate of water molecules when the channel is open, by locally modifying the orientation of His173. Furthermore, the calculated permeation rates of a fully open channel are in good agreement with the reported experimental value.
| S-EPMC3614956 | biostudies-literature
Project description:Many protein kinases are characterized by at least two structural forms corresponding to the highest level of activity (active) and low or no activity, (inactive). Further, protein dynamics is an important consideration in understanding the molecular and mechanistic basis of enzyme function. In this work, we use protein kinase A (PKA) as the model system and perform microsecond range molecular dynamics (MD) simulations on six variants which differ from one another in terms of active and inactive form, with or without bound ligands, C-terminal tail and phosphorylation at the activation loop. We find that the root mean square fluctuations in the MD simulations are generally higher for the inactive forms than the active forms. This difference is statistically significant. The higher dynamics of inactive states has significant contributions from ATP binding loop, catalytic loop, and αG helix. Simulations with and without C-terminal tail show this differential dynamics as well, with lower dynamics both in the active and inactive forms if C-terminal tail is present. Similarly, the dynamics associated with the inactive form is higher irrespective of the phosphorylation status of Thr 197. A relatively stable stature of active kinases may be better suited for binding of substrates and detachment of the product. Also, phosphoryl group transfer from ATP to the phosphosite on the substrate requires precise transient coordination of chemical entities from three different molecules, which may be facilitated by the higher stability of the active state.
| S-EPMC6371217 | biostudies-literature
Project description:Kinases control many aspects of cellular signaling and are therefore therapeutic targets for numerous disease states. Monitoring the conformational changes that drive activation and inactivation of the catalytic kinase core is a challenging experimental problem due to the dynamic nature of these enzymes. We apply [(13)C] reductive methylation to chemically introduce NMR-active nuclei into unlabeled protein kinases. The results demonstrate that solution NMR spectroscopy can be used to monitor specific changes in the chemical environment of structurally important lysines in a [(13)C]-methylated kinase as it shifts from the inactive to active state. This approach provides a solution based method to complement X-ray crystallographic data and can be applied to nearly any kinase, regardless of size or method of production.
| S-EPMC4301085 | biostudies-literature
Project description:The competitive inhibitor cocaine and the non-competitive inhibitor ibogaine induce different conformational states of the human serotonin transporter. It has been shown from accessibility experiments that cocaine mainly induces an outward-facing conformation, while the non-competitive inhibitor ibogaine, and its active metabolite noribogaine, have been proposed to induce an inward-facing conformation of the human serotonin transporter similar to what has been observed for the endogenous substrate, serotonin. The ligand induced conformational changes within the human serotonin transporter caused by these three different types of ligands, substrate, non-competitive and competitive inhibitors, are studied from multiple atomistic molecular dynamics simulations initiated from a homology model of the human serotonin transporter. The results reveal that diverse conformations of the human serotonin transporter are captured from the molecular dynamics simulations depending on the type of the ligand bound. The inward-facing conformation of the human serotonin transporter is reached with noribogaine bound, and this state resembles a previously identified inward-facing conformation of the human serotonin transporter obtained from molecular dynamics simulation with bound substrate, but also a recently published inward-facing conformation of a bacterial homolog, the leucine transporter from Aquifex Aoelicus. The differences observed in ligand induced behavior are found to originate from different interaction patterns between the ligands and the protein. Such atomic-level understanding of how an inhibitor can dictate the conformational response of a transporter by ligand binding may be of great importance for future drug design.
| S-EPMC3680404 | biostudies-literature
Project description:Investigating ligand-regulated allosteric coupling between protein domains is fundamental to understand cell-life regulation. The Hsp70 family of chaperones represents an example of proteins in which ATP binding and hydrolysis at the Nucleotide Binding Domain (NBD) modulate substrate recognition at the Substrate Binding Domain (SBD). Herein, a comparative analysis of an allosteric (Hsp70-DnaK) and a non-allosteric structural homolog (Hsp110-Sse1) of the Hsp70 family is carried out through molecular dynamics simulations, starting from different conformations and ligand-states. Analysis of ligand-dependent modulation of internal fluctuations and local deformation patterns highlights the structural and dynamical changes occurring at residue level upon ATP-ADP exchange, which are connected to the conformational transition between closed and open structures. By identifying the dynamically responsive protein regions and specific cross-domain hydrogen-bonding patterns that differentiate Hsp70 from Hsp110 as a function of the nucleotide, we propose a molecular mechanism for the allosteric signal propagation of the ATP-encoded conformational signal.
| S-EPMC3531320 | biostudies-literature
Project description:Acetylcholinesterase, with a deep, narrow active-site gorge, attracts enormous interest due to its particularly high catalytic efficiency and its inhibitors used for treatment of Alzheimer's disease. To facilitate the massive pass-through of the substrate and inhibitors, "breathing" motions to modulate the size of the gorge are an important prerequisite. However, the molecular mechanism that governs such motions is not well explored. Here, to systematically investigate intrinsic motions of the enzyme, we performed microsecond molecular dynamics simulations on the monomer and dimer of Torpedo californica acetylcholinesterase (TcAChE) as well as the complex of TcAChE bound with the drug E2020. It has been revealed that protein-ligand interactions and dimerization both keep the gorge in bulk, and opening events of the gorge increase dramatically compared to the monomer. Dynamics of three subdomains, S3, S4 and the Ω-loop, are tightly associated with variations of the gorge size while the dynamics can be changed by ligand binding or protein dimerization. Moreover, high correlations among these subdomains provide a basis for remote residues allosterically modulating the gorge motions. These observations are propitious to expand our understanding of protein structure and function as well as providing clues for performing structure-based drug design.
| S-EPMC5468367 | biostudies-literature