Probing light chain mutation effects on thrombin via molecular dynamics simulations and machine learning

The Light Chain Allosterically Enhances the Protease Activity of Murine Urokinase-Type Plasminogen Activator

The active form of the murine urokinase-type plasminogen activator (muPA) is formed by a 27-residue disordered light chain connecting the amino-terminal fragment (ATF) with the serine protease domain. The two chains are tethered by a disulfide bond between C1CT in the disordered light chain and C122CT in the protease domain. Previous work showed that the presence of the disordered light chain affected the inhibition of the protease domain by antibodies. Here we show that the disordered light chain induced a 3.7-fold increase in kcat of the protease domain of muPA. In addition, hydrogen-deuterium exchange mass spectrometry (HDX-MS) and accelerated molecular dynamics (AMD) were performed to identify the interactions between the disordered light chain and the protease domain. HDX-MS revealed that the light chain is contacting the 110s, the turn between the β10- and β11-strand, and the β7-strand. A reduction in deuterium uptake was also observed in the activation loop, the 140s loop and the 220s loop, which forms the S1-specificty pocket where the substrate binds. These loops are further away from where the light chain seems to be interacting with the protease domain. Our results suggest that the light chain most likely increases the activity of muPA by allosterically favoring conformations in which the specificity pocket is formed. We propose a model by which the allostery would be transmitted through the β-strands of the β-barrels to the loops on the other side of the protease domain.

Allosteric Modulation of Thrombin by Thrombomodulin: Insights from Logistic Regression and Statistical Analysis of Molecular Dynamics Simulations

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.

Unraveling the Role of Hydrogen Bonds in Thrombin via Two Machine Learning Methods

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.

Machine learning and protein allostery

The fundamental biological importance and complexity of allosterically regulated proteins stem from their central role in signal transduction and cellular processes. Recently, machine-learning approaches have been developed and actively deployed to facilitate theoretical and experimental studies of protein dynamics and allosteric mechanisms. In this review, we survey recent developments in applications of machine-learning methods for studies of allosteric mechanisms, prediction of allosteric effects and allostery-related physicochemical properties, and allosteric protein engineering. We also review the applications of machine-learning strategies for characterization of allosteric mechanisms and drug design targeting SARS-CoV-2. Continuous development and task-specific adaptation of machine-learning methods for protein allosteric mechanisms will have an increasingly important role in bridging a wide spectrum of data-intensive experimental and theoretical technologies.

Discovery of novel VEGFR2-TK inhibitors by phthalimide pharmacophore based virtual screening, molecular docking, MD simulation and DFT

Currently, numerous potent chemotherapeutic agents are available in the market but most of them show poor pharmacokinetics, lethal effects and drug resistance during their enduring use. The increased cancer cases, deaths and need of better treatment stimulates us to give newer lifesaving anticancer drugs. The phthalimide derivatives are structurally diverse and exert potential anticancer activity. In this regard, the 3D QSAR Pharmacophore model was developed and validated using fifty-eight phthalimide derivatives. The validation parameters corroborated the reliability and statistical robustness of CEASER Hypo 1. Three databases-NCI Open, Drug Bank, and Asinex were submitted to ADMET and drug-like filtering; 117893 drug-like compounds were mapped on CEASER Hypo 1; and 362 hits with IC50 <1 µM were discovered. These hits were docked on VEGFR2-TK, and in the form of results fifteen hits exhibited greater affinity than sorafenib. The top lead ASN 03206926 was subjected for MD simulation (100 ns) and RMSD, Rg, RMSF, number of hydrogen bonds, and SASA verified that the complex was stable, rigid and highly compact. Results demonstrated GLU885, PHE918, CYS919, LYS920, HIS1026, CYS1045, ASP1046 are the essential residues for favourable interactions. The binding free energy calculations support the affinity and stability revealed by docking and MD simulation. The DFT calculations, negative binding energy and lower HOMO-LUMO band gap revealed that the process is spontaneous and ASN 03206926 is very reactive. Following extensive analysis we suggest that the ASN 03206926 might be employed as a new VEGFR2-TK inhibitor for the treatment of breast and VEGFR2-TK associated cancers.Communicated by Ramaswamy H. Sarma.

Simulations suggest double sodium binding induces unexpected conformational changes in thrombin

Thrombin is a Na[Formula: see text]-activated serine protease existing in two forms targeted to procoagulant and anticoagulant activities, respectively. There is one Na[Formula: see text]-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[Formula: see text]-binding sites in thrombin and that Na[Formula: see text] 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[Formula: see text] binding modes on thrombin. From the root-mean-square fluctuations (RMSF) for the [Formula: see text]-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[Formula: see text] binding to thrombin might be a multiple-step mechanism rather than a 2-step mechanism.

Exosite Binding in Thrombin: A Global Structural/Dynamic Overview of Complexes with Aptamers and Other Ligands

Thrombin is the key enzyme of the entire hemostatic process since it is able to exert both procoagulant and anticoagulant functions; therefore, it represents an attractive target for the developments of biomolecules with therapeutic potential. Thrombin can perform its many functional activities because of its ability to recognize a wide variety of substrates, inhibitors, and cofactors. These molecules frequently are bound to positively charged regions on the surface of protein called exosites. In this review, we carried out extensive analyses of the structural determinants of thrombin partnerships by surveying literature data as well as the structural content of the Protein Data Bank (PDB). In particular, we used the information collected on functional, natural, and synthetic molecular ligands to define the anatomy of the exosites and to quantify the interface area between thrombin and exosite ligands. In this framework, we reviewed in detail the specificity of thrombin binding to aptamers, a class of compounds with intriguing pharmaceutical properties. Although these compounds anchor to protein using conservative patterns on its surface, the present analysis highlights some interesting peculiarities. Moreover, the impact of thrombin binding aptamers in the elucidation of the cross-talk between the two distant exosites is illustrated. Collectively, the data and the work here reviewed may provide insights into the design of novel thrombin inhibitors.

Light Chain Mutation Effects on the Dynamics of Thrombin

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.

Identify thrombin inhibitor with novel skeleton based on virtual screening study

Virtual screening refers to the screening of active compounds based on a small-molecule database. This procedure can rapidly select active compounds with pharmaceutical properties from millions of molecules, thus considerably reducing the number of experimental screening compounds and cost of drug development and shortening the research cycle. In this paper, a pharmacophore screening method was used for virtual screening to determine new scaffold compounds with potential anticoagulant activities. The pharmacophore model (Model_01-20) was constructed in SYBYL-X 2.0 based on dabigatran derivatives (D1-D9) with micromolar to nanomolar activities and tested by decoy test method. Model_01 was selected to screen more than 1600 million compounds in the Zinc 12.0 database. Furtherly, molecular docking analysis and ADME prediction were conducted on more than 100,000 screened compounds. Finally, two compounds (Z-19 and Z-29) were selected for anticoagulant activity test in vitro, Compound Z-29 with tryptophan aurone structure was found possess anticoagulant effect and its IC₅₀ = 22.9 ± 6.88 μM. ADME prediction results show that compound Z-29 features a high intestinal absorption rate, which is valuable for further in-depth research. The research results of this paper can be used for further structural modification and optimisation to guide the design and provide new ideas and methods for the discovery of new thrombin inhibitors.Communicated by Ramaswamy H. Sarma.

An investigation into the allosteric mechanism of GPCR A2A adenosine receptor with trajectory-based information theory and complex network model

G protein-coupled receptors (GPCRs), a large superfamily of transmembrane (TM) proteins, allosterically transduce the signal of ligand binding in the extracellular (EC) domain to couple to effector proteins in the intracellular (IC) domain, therefore forming the largest class of drug targets. The A_2A adenosine receptor (A_2AAR), a class-A GPCR, has been extensively studied as it offers numerous possibilities for therapeutic applications. However, the mechanism of allosteric communication between EC and IC domains is not completely clear. In this work, we utilize torsional mutual information to quantify the correlated motions of residue pairs from its molecular dynamics (MD) simulation trajectories, and further use the complex network model to obtain allosteric pipelines and hubs. The identified allosteric communication pipelines mainly transmit the signal from EC domain to the cytoplasmic ends of TM helix 5 (TM5), TM6 and TM7. The allosteric hubs, mostly located at TM5, TM6 and TM7, play an important role in mediating allosteric signal transmission to keep the receptor rigid and prevent G protein from binding to IC domain, which can explain the reason why their mutations distant from ligand-binding site do not affect the ligand binding affinity but affect the ligand efficacy. Additionally, we identify the key residues located in antagonist ZM241385 binding pocket which mediate multiple allosteric pathways and have been experimentally proven to play a critical role in affecting the ligand potency. This study is helpful for understanding the allosteric communication mechanism of A_2AAR, and can provide valuable information for the structure-based drug design of GPCRs.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulations of human α-thrombin in different structural contexts: evidence for an aptamer-guided cooperation between the two exosites

Human α-thrombin (thrombin) is a multifunctional enzyme that plays a pivotal role in the coagulation pathway. Thrombin activity can be effectively modulated by G-quadruplex-based oligonucleotide aptamers that specifically interact with the two positively charged regions (exosites I and II) on the protein surface. Although insightful atomic-level snapshots of the recognition between thrombin and aptamers have been recently achieved through crystallographic analyses, some dynamic aspects of this interaction have not been fully characterized. We here report molecular dynamics simulations of thrombin in different association states: ligand-free and binary/ternary complexes with the aptamers TBA (at exosite I) and HD22_27mer (at exosite II). The simulations carried out on the binary and ternary complexes formed by thrombin with these aptamers provide a dynamic view of the interactions that stabilize them in a crystal-free environment. Interestingly, the analysis of the dynamics of the exosites in different thrombin binding states clearly indicates that the HD22_27mer binding at the exosite II favours conformations of exosite I that are prone to the TBA binding. Similar effects are observed upon the binding of TBA to the exosite I. These observations provide an atomic-level picture of the exosite inter-communication in thrombin and explain the experimentally detected cooperativity of the TBA/HD22_27mer binding.

Mechanism and inhibition kinetics of peptide P13 as thrombin inhibitor

Excessive coagulation can easily lead to arterial and venous thrombosis, which is the main reason for the evolution of myocardial infarction and cerebrovascular accidents. As a key coagulation factor for the coagulation pathway, thrombin has become a remarkable target for the control of thrombosis. The synthesized peptide P13 with amino acid sequence of N-RGDAGFAGDDAPR was expected to be an inhibitor with higher antithrombotic activity. The results showed that the IC₅₀ (50% inhibition of thrombin activity) of the peptide P13 was determined by colorimetric method to be 115 µM. And enzyme kinetic experiments showed that P13 was a competitive inhibitor of thrombin with K_i = 106 µM. Fluorescence spectra and three-dimensional fluorescence showed that P13 could alter the secondary structure of thrombin and the microenvironment of certain chromogenic amino acids. P13 can spontaneously bind with thrombin exosite 1 in the form of 1:1 mainly through hydrogen bonding and van der Waals force. And the optimal docking mode of P13 and thrombin was revealed by molecular docking with "-CDOCKER_Energy" of 178.679 kcal mol^-1. This study revealed P13 may become a potential anticoagulant drug widely used after further studies in preclinical and clinical trials.

Interaction of the synthetic antithrombotic peptide P10 with thrombin: a spectroscopy study

Thrombin is a critical serine protease in the coagulation system and is widely used as a target protein for antithrombotics. Spectroscopic analysis is a simple and effective method that is used to study the interaction between small molecules and proteins. In this study, the interactions of a potential antithrombotic peptide AGFAGDDAPR (P10) with thrombin were investigated by fluorescence spectroscopy, UV-vis spectroscopy, circular dichroism, Fourier-transform infrared spectroscopy and Raman spectroscopy, respectively. The results showed that the peptide P10 bonded to thrombin via hydrogen bonding and van der Waals forces, resulting in fluorescence quenching. And, the secondary structure of thrombin changed, the β-sheet decreased, and the random coil increased. The peptide P10 bonded to proline and lysine, and changed the space structure of thrombin, resulting in inhibition of thrombin activity. The results contributed to exploration of the mechanism of this potential antithrombotic drug interaction with thrombin in order to provide a preliminary understanding of the pharmacodynamic properties of P10.

Molecular interaction between bi-antennary phenylboronic acid and sialic acid using density functional theory and multi-time scale trajectories

Communicated by Ramaswamy H. Sarma.

Na+-binding modes involved in thrombin's allosteric response as revealed by molecular dynamics simulations, correlation networks and Markov modeling

The monovalent sodium ion (Na+) is a critical modulator of thrombin. However, the mechanism of thrombin's activation by Na+ has been widely debated for more than twenty years. Details of the linkage between thrombin and Na+ remain vague due to limited temporal and spatial resolution in experiments. In this work, we combine microsecond scale atomic-detailed molecular dynamics simulations with correlation network analyses and hidden Markov modeling to probe the detailed thermodynamic and kinetic picture of Na+-binding events and their resulting allosteric responses in thrombin. We reveal that ASP189 and ALA190 comprise a stable Na+-binding site (referred as "inner" Na+-binding site) along with the previously known one (referred as "outer" Na+-binding site). The corresponding newly identified Na+-binding mode introduces significant allosteric responses in thrombin's regulatory regions by stabilizing selected torsion angles of residues responsive to Na+-binding. Our Markov model indicates that the bound Na+ prefers to transfer between the two Na+-binding sites when an unbinding event takes place. These results suggest a testable hypothesis of a substrate-driven Na+ migration (ΔG ∼ 1.7 kcal mol-1) from the "inner" Na+-binding site to the "outer" one during thrombin's catalytic activities. The binding of a Na+ ion at the "inner" Na+-binding site should be inferred as a prerequisite for thrombin's efficient recognition to the substrate, which opens a new angle for our understanding of Na+-binding's allosteric activation on thrombin and sheds light on detailed processes in thrombin's activation.