Molecular dynamics simulations have become increasingly useful in studying biological systems of biomedical interest, and not just in the study of model or toy systems. In this article, the methods ...and principles of all-atom molecular dynamics will be elucidated with several examples provided of their utility to investigators interested on drug discovery.
Hydrogen bonds play a critical role in the folding and stability of proteins, such as proteins and nucleic acids, by providing strong and directional interactions. They help to maintain the secondary ...and 3D structure of proteins, and structural changes in these molecules often result from the formation or breaking of hydrogen bonds. To gain insights into these hydrogen bonding networks, we applied two machine learning models - a logistic regression model and a decision tree model - to study four variants of thrombin: wild-type, ΔK9, E8K, and R4A. Our results showed that both models have their unique advantages. The logistic regression model highlighted potential key residues (GLU295) in thrombin’s allosteric pathways, while the decision tree model identified important hydrogen bonding motifs. This information can aid in understanding the mechanisms of folding in proteins and has potential applications in drug design and other therapies. The use of these two models highlights their usefulness in studying hydrogen bonding networks in proteins.
Thrombin plays an important role in the process of hemostasis and blood coagulation. Studies in thrombin can help us find ways to treat cancer because thrombin is able to reduce the characteristic ...hypercoagulability of cancer. Thrombin is composed of two chains, the light chain and the heavy chain. The function of the heavy chain has been largely explored, while the function of the light chain was obscured until several disease-associated mutations in the light chain come to light. In this study, we want to explore the dynamic and conformation effects of mutations on the light chain further to determine possible associations between mutation, conformational changes, and disease. The study, which is a follow-up for our studies on apo thrombin and the mutant, ΔK9, mainly focuses on the mutants E8K and R4A. E8K is a disease-associated mutation, and R4A is used to study the role of Arg4, which is suggested experimentally to play a critical role for thrombin’s catalytic activities. We performed five all-atom one microsecond-scale molecular dynamics (MD) simulations for both E8K and R4A, and quantified the changes in the conformational ensemble of the mutants. From the root-mean-square fluctuations (RMSF) for the α-carbons, we find that the atomic fluctuations change in the mutants in the 60s loop and γ loop. The correlation coefficients for the α-carbons indicate that the correlation relation for atom-pairs in the protein is also impacted. The clustering analysis and the principal component analysis (PCA) consistently tell us that the catalytic pocket and the regulatory loops are destabilized by the mutations. We also find that there are two binding modes for Na+ by clustering the vector difference between the Na+ ions and the 220s loop. After further analysis, we find that there is a relation between the Na+ binding and the rigidification of the γ loop, which may shed light on the mysterious role of the γ loop in thrombin.
Mutations in the TREX1 3' → 5' exonuclease are associated with a spectrum of autoimmune disease phenotypes in humans and mice. Failure to degrade DNA activates the cGAS-STING DNA-sensing pathway ...signaling a type-I interferon (IFN) response that ultimately drives immune system activation. TREX1 and the cGAS-STING DNA-sensing pathway have also been implicated in the tumor microenvironment, where TREX1 is proposed to degrade tumor-derived DNA that would otherwise activate cGAS-STING. If tumor-derived DNA were not degraded, the cGAS-STING pathway would be activated to promote IFN-dependent antitumor immunity. Thus, we hypothesize TREX1 exonuclease inhibition as a novel immunotherapeutic strategy. We present data demonstrating antitumor immunity in the TREX1 D18N mouse model and discuss theory surrounding the best strategy for TREX1 inhibition. Potential complications of TREX1 inhibition as a therapeutic strategy are also discussed.
Thrombin is a Na
+
-activated serine protease existing in two forms targeted to procoagulant and anticoagulant activities, respectively. There is one Na
+
-binding site that has been the focus of the ...study of the thrombin. However, molecular dynamics (MD) simulations suggest that there might be actually two Na
+
-binding sites in thrombin and that Na
+
ions can even bind to two sites simultaneously. In this study, we performed 12 independent 2-µs all-atom MD simulations for the wild-type (WT) thrombin and we studied the effects of the different Na
+
binding modes on thrombin. From the root-mean-square fluctuations (RMSF) for the
α
-carbons, we see that the atomic fluctuations mainly change in the 60s, 170s, and 220s loops, and the connection (residue 167 to 170). The correlation matrices for different binding modes suggest regions that may play an important role in thrombin’s allosteric response and provide us a possible allosteric pathway for the sodium binding. Amorim-Hennig (AH) clustering tells us how the structure of the regions of interest changes on sodium binding. Principal component analysis (PCA) shows us how the different regions of thrombin change conformation together with sodium binding. Solvent-accessible surface area (SASA) exposes the conformational change in exosite I and catalytic triad. Finally, we argue that the double binding mode might be an inactive mode and that the kinetic scheme for the Na
+
binding to thrombin might be a multiple-step mechanism rather than a 2-step mechanism.
Thrombomodulin (TM), a transmembrane receptor integral to the anticoagulant pathway, governs thrombin’s substrate specificity via interaction with thrombin’s anion-binding exosite I. Despite its ...established role, the precise mechanisms underlying this regulatory function are yet to be fully unraveled. In this study, we deepen the understanding of these mechanisms through eight independent 1 μs all-atom simulations, analyzing thrombin both in its free form and when bound to TM fragments TM456 and TM56. Our investigations revealed distinct and significant conformational changes in thrombin mediated by the binding of TM56 and TM456. While TM56 predominantly influences motions within exosite I, TM456 orchestrates coordinated alterations across various loop regions, thereby unveiling a multifaceted modulatory role that extends beyond that of TM56. A highlight of our study is the identification of critical hydrogen bonds that undergo transformations during TM56 and TM456 binding, shedding light on the pivotal allosteric influence exerted by TM4 on thrombin’s structural dynamics. This work offers a nuanced appreciation of TM’s regulatory role in blood coagulation, paving the way for innovative approaches in the development of anticoagulant therapies and expanding the horizons in oncology therapeutics through a deeper understanding of molecular interactions in the coagulation pathway.
Thrombin is a multifunctional enzyme that plays an important role in blood coagulation, cell growth, and metastasis. Depending upon the binding of sodium ions, thrombin presents significantly ...different enzymatic activities. In the environment with sodium ions, thrombin is highly active in cleaving the coagulated substrates and this is referred to as the "fast" form; in the environment without sodium ions, thrombin turns catalytically less active and is in the "slow" form. Although many experimental studies over the last two decades have attempted to reveal the structural and kinetic differences between these two forms, it remains vague and disputed how the functional switch between the "fast" and "slow" forms is mediated by Na
cations. In this work, we employ microsecond-scale all-atom molecular dynamics simulations to investigate the differences in the structural ensembles in sodium-bound/unbound and potassium-bound/unbound thrombin. Our calculations indicate that the regulatory regions, including the 60s, γ loops, and exosite I and II, are primarily affected by both the bound and unbound cations. Conformational free energy surfaces, estimated from principal component analysis, further reveal the existence of multiple conformational states. The binding of a cation introduces changes in the distribution of these states. Through comparisons with potassium-binding, the binding of sodium ions appears to shift the population toward conformational states that might be catalytically favorable. Our study of thrombin in the presence of sodium/potassium ions suggests Na
-mediated generalized allostery is the mechanism of thrombin's functional switch between the "fast" and "slow" forms.