The natural mutation by deletion of Lys9 in the thrombin A-chain affects the pKa value of catalytic residues, the overall enzyme's stability and conformational transitions linked to Na+ binding

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.

A heterozygous pathogenic variant in the ATP6V1A gene triggering epilepsy in a large Chinese pedigree

Epilepsy is one of the most common disorders in children, with an incidence rate of approximately 5%. Although an increasing number of genes have been demonstrated to be pathogenic factors in epilepsy, evidence for a potential pathogenic role of ATP6V1A remains limited. Herein, the clinical and genetic data of a 5-year-old boy who experienced seizures at 9 months of age are collected. Genetic variants are screened using whole-exome sequencing (WES), and the effects of the candidate variants are further validated at both the RNA and protein levels. WES reveals a heterozygous variant [NM_001690.4: c .1132 C>T, p.Leu378Phe] of the ATP6V1A gene. This variant is not reported in the public database, but is predicted to be deleterious by multiple software packages, and classified as a variant of unknown significance following the American College of Medical Genetics and Genomics guidelines. Quantitative PCR and western blotting further confirm its down-regulatory role in both the RNA and protein expression of ATP6V1A. This case report confirms the pathogenicity of ATP6V1A in epilepsy with solid experimental evidence, thereby expanding the phenotype spectrum of ATP6V1A variants. More importantly, we show that seizures triggered by ATP6V1A variants could be controlled by Levetiracetam, crucially rescuing the development of the patient.

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.

Functional impact of cancer patient-associated Bcl-xL mutations

Bcl-xL, an antiapoptotic protein, is frequently overexpressed in cancer to promote survival of tumor cells. However, we have previously shown that Bcl-xL promotes migration, invasion, and metastasis independent of its antiapoptotic function in mitochondria. The pro-metastatic function of Bcl-xL may require its translocation into the nucleus. Besides overexpression, patient-associated mutations of Bcl-xL have been identified in large-scale cancer genomics projects. Understanding the functions of these mutations will guide the development of precision medicine. Here, we selected four patient-associated Bcl-xL mutations, R132W, N136K, R165W, and A201T, to investigate their impacts on antiapoptosis, migration, and nuclear translocation. We found that all four mutation proteins could be detected in both the nucleus and cytosol. Although all four mutations disrupted the antiapoptosis function, one of these mutants, N136K, significantly improved the ability to promote cell migration. These data suggest the importance of developing novel Bcl-xL inhibitors to ablate both antiapoptotic and pro-metastatic functions of Bcl-xL in cancer.

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

Thrombin is a key component for chemotherapeutic and antithrombotic therapy development. As the physiologic and pathologic roles of the light chain still remain vague, here, we continue previous efforts to understand the impacts of the disease-associated single deletion of LYS9 in the light chain. By combining supervised and unsupervised machine learning methodologies and more traditional structural analyses on data from 10 μs molecular dynamics simulations, we show that the conformational ensemble of the ΔK9 mutant is significantly perturbed. Our analyses consistently indicate that LYS9 deletion destabilizes both the catalytic cleft and regulatory functional regions and result in some conformational changes that occur in tens to hundreds of nanosecond scaled motions. We also reveal that the two forms of thrombin each prefer a distinct binding mode of a Na+ ion. We expand our understanding of previous experimental observations and shed light on the mechanisms of the LYS9 deletion associated bleeding disorder by providing consistent but more quantitative and detailed structural analyses than early studies in literature. With a novel application of supervised learning, i.e. the decision tree learning on the hydrogen bonding features in the wild-type and ΔK9 mutant forms of thrombin, we predict that seven pairs of critical hydrogen bonding interactions are significant for establishing distinct behaviors of wild-type thrombin and its ΔK9 mutant form. Our calculations indicate the LYS9 in the light chain has both localized and long-range allosteric effects on thrombin, supporting the opinion that light chain has an important role as an allosteric effector.

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.

Homology modeling and in silico prediction of Ulcerative colitis associated polymorphisms of NOD1

Cytosolic pattern recognition receptors play key roles in innate immune response. Nucleotide binding and oligomerisation domain containing protein 1 (NOD1) belonging to the Nod-like receptor C (NLRC) sub-family of Nod-like receptors (NLRs) is important for detection and clearance of intra-cellular Gram negative bacteria. NOD1 is involved in activation of pro-inflammatory pathways. Limited structural data is available for NOD1. Using different templates for each domain of NOD1, we determined the full-length homology model of NOD1. ADP binding amino acids within the nucleotide binding domain (NBD) of NOD1 were also predicted. Key residues in inter-domain interaction were identified by sequence comparison with Oryctolagus cuniculus NOD2, a related protein. Interactions between NBD and winged helix domain (WHD) were found to be conserved in NOD1. Functional and structural effect of single nucleotide polymorphisms within the NOD1 NBD domain associated with susceptibility risk to Ulcerative colitis (UC), an inflammatory disorder of the colon was evaluated by in silico studies. Mutations W219R and L349P were predicted to be damaging and disease associated by prediction programs SIFT, PolyPhen2, PANTHER, SNP&GO, PhD SNP and SNAP2. We further validated the effect of W219R and L349P mutation on NOD1 function in vitro. Elevated mRNA expression of pro-inflammatory cytokines IL8 and IL-1β was seen as compared to the wild type NOD1 in intestinal epithelial cell line HT29 when stimulated with NOD1 ligand. Thus, these mutations may indeed have a bearing on pathogenesis of inflammation during UC.

Identification of Deleterious Mutations in Myostatin Gene of Rohu Carp (Labeo rohita) Using Modeling and Molecular Dynamic Simulation Approaches

The myostatin (MSTN) is a known negative growth regulator of skeletal muscle. The mutated myostatin showed a double-muscular phenotype having a positive significance for the farmed animals. Consequently, adequate information is not available in the teleosts, including farmed rohu carp, Labeo rohita. In the absence of experimental evidence, computational algorithms were utilized in predicting the impact of point mutation of rohu myostatin, especially its structural and functional relationships. The four mutations were generated at different positions (p.D76A, p.Q204P, p.C312Y, and p.D313A) of MSTN protein of rohu. The impacts of each mutant were analyzed using SIFT, I-Mutant 2.0, PANTHER, and PROVEAN, wherein two substitutions (p.D76A and p.Q204P) were predicted as deleterious. The comparative structural analysis of each mutant protein with the native was explored using 3D modeling as well as molecular-dynamic simulation techniques. The simulation showed altered dynamic behaviors concerning RMSD and RMSF, for either p.D76A or p.Q204P substitution, when compared with the native counterpart. Interestingly, incorporated two mutations imposed a significant negative impact on protein structure and stability. The present study provided the first-hand information in identifying possible amino acids, where mutations could be incorporated into MSTN gene of rohu carp including other carps for undertaking further in vivo studies.

Analyzing effects of naturally occurring missense mutations

Single-point mutation in genome, for example, single-nucleotide polymorphism (SNP) or rare genetic mutation, is the change of a single nucleotide for another in the genome sequence. Some of them will produce an amino acid substitution in the corresponding protein sequence (missense mutations); others will not. This paper focuses on genetic mutations resulting in a change in the amino acid sequence of the corresponding protein and how to assess their effects on protein wild-type characteristics. The existing methods and approaches for predicting the effects of mutation on protein stability, structure, and dynamics are outlined and discussed with respect to their underlying principles. Available resources, either as stand-alone applications or webservers, are pointed out as well. It is emphasized that understanding the molecular mechanisms behind these effects due to these missense mutations is of critical importance for detecting disease-causing mutations. The paper provides several examples of the application of 3D structure-based methods to model the effects of protein stability and protein-protein interactions caused by missense mutations as well.

Predicting the functional impact of protein mutations: application to cancer genomics

As large-scale re-sequencing of genomes reveals many protein mutations, especially in human cancer tissues, prediction of their likely functional impact becomes important practical goal. Here, we introduce a new functional impact score (FIS) for amino acid residue changes using evolutionary conservation patterns. The information in these patterns is derived from aligned families and sub-families of sequence homologs within and between species using combinatorial entropy formalism. The score performs well on a large set of human protein mutations in separating disease-associated variants (∼19 200), assumed to be strongly functional, from common polymorphisms (∼35 600), assumed to be weakly functional (area under the receiver operating characteristic curve of ∼0.86). In cancer, using recurrence, multiplicity and annotation for ∼10 000 mutations in the COSMIC database, the method does well in assigning higher scores to more likely functional mutations ('drivers'). To guide experimental prioritization, we report a list of about 1000 top human cancer genes frequently mutated in one or more cancer types ranked by likely functional impact; and, an additional 1000 candidate cancer genes with rare but likely functional mutations. In addition, we estimate that at least 5% of cancer-relevant mutations involve switch of function, rather than simply loss or gain of function.

Thrombin a-chain: activation remnant or allosteric effector?

Although prothrombin is one of the most widely studied enzymes in biology, the role of the thrombin A-chain has been neglected in comparison to the other domains. This paper summarizes the current data on the prothrombin catalytic domain A-chain region and the subsequent thrombin A-chain. Attention is given to biochemical characterization of naturally occurring prothrombin A-chain mutations and alanine scanning mutants in this region. While originally considered to be simply an activation remnant with little physiologic function, the thrombin A-chain is now thought to play a role as an allosteric effector in enzymatic reactions and may also be a structural scaffold to stabilize the protease domain.

Modeling effects of human single nucleotide polymorphisms on protein-protein interactions

A large set of three-dimensional structures of 264 protein-protein complexes with known nonsynonymous single nucleotide polymorphisms (nsSNPs) at the interface was built using homology-based methods. The nsSNPs were mapped on the proteins' structures and their effect on the binding energy was investigated with CHARMM force field and continuum electrostatic calculations. Two sets of nsSNPs were studied: disease annotated Online Mendelian Inheritance in Man (OMIM) and nonannotated (non-OMIM). It was demonstrated that OMIM nsSNPs tend to destabilize the electrostatic component of the binding energy, in contrast with the effect of non-OMIM nsSNPs. In addition, it was shown that the change of the binding energy upon amino acid substitutions is not related to the conservation of the net charge, hydrophobicity, or hydrogen bond network at the interface. The results indicate that, generally, the effect of nsSNPs on protein-protein interactions cannot be predicted from amino acids' physico-chemical properties alone, since in many cases a substitution of a particular residue with another amino acid having completely different polarity or hydrophobicity had little effect on the binding energy. Analysis of sequence conservation showed that nsSNP at highly conserved positions resulted in a large variance of the binding energy changes. In contrast, amino acid substitutions corresponding to nsSNPs at nonconserved positions, on average, were not found to have a large effect on binding affinity. pKa calculations were performed and showed that amino acid substitutions could change the wild-type proton uptake/release and thus resulting in different pH-dependence of the binding energy.

Predicting potentially functional SNPs in drug-response genes

SNPs are known to contribute to variations in drug response and there are more than 14 million polymorphisms spanning the human genome. However, not all of these SNPs are functional. It would be impractical and costly to evaluate every individual SNP for functionality experimentally. Consequently, one of the major challenges for researchers has been to seek out functional SNPs from all the SNPs in the human genome. In silico or bioinformatic methods are economical, less labor intensive, yet powerful approaches to filter out potentially functional SNPs in drug-response genes for further study. This allows researchers to prioritize which SNPs to subsequently evaluate experimentally for drug-response studies, as well as potentially providing insights into possible mechanisms underlying how SNPs may affect drug-response genes.

Role of the A chain in thrombin function

The A chain of thrombin is covalently linked to the catalytic B chain but is separate from any known epitope for substrate recognition. In this study we present the results of the Ala replacement of 12 charged residues controlling the stability of the A chain and its interaction with the B chain. Residues Arg4 and Glu8 play a significant role in substrate recognition, even though they are located > 20 A away from residues of the catalytic triad, the primary specificity pocket and the Na+ site. The R4A mutation causes significant perturbation of Na+ binding, fibrinogen clotting and PAR1 cleavage, but modest reduction of protein C activation in the presence of thrombomodulin. These findings challenge our current paradigm of thrombin structure-function relations focused exclusively on the properties of the catalytic B chain, and explain why certain naturally occurring mutations of the A chain cause serious bleeding.

Role of A-chain in functioning of the active site of human alpha-thrombin

This review summarizes current data suggesting that A-chain of the human alpha-thrombin molecule plays a role of allosteric effector in catalytic reactions with various substrates. Special attention is paid to the relationship between A-chain structure and catalytic activity of thrombin. The existence of this relationship is based on studies of natural mutation of A-chain of the alpha-thrombin molecule. Use of molecular and essential dynamics confirmed the role of A-chain in changes of conformation and catalytic properties of this enzyme; these changes involve residues located in the specificity sites and some inserting loops. Current knowledge on structure and properties of thrombin can be used for the development of new antithrombin agents.

Thrombin

Thrombin is a Na+-activated, allosteric serine protease that plays opposing functional roles in blood coagulation. Binding of Na+ is the major driving force behind the procoagulant, prothrombotic and signaling functions of the enzyme, but is dispensable for cleavage of the anticoagulant protein C. The anticoagulant function of thrombin is under the allosteric control of the cofactor thrombomodulin. Much has been learned on the mechanism of Na+ binding and recognition of natural substrates by thrombin. Recent structural advances have shed light on the remarkable molecular plasticity of this enzyme and the molecular underpinnings of thrombin allostery mediated by binding to exosite I and the Na+ site. This review summarizes our current understanding of the molecular basis of thrombin function and allosteric regulation. The basic information emerging from recent structural, mutagenesis and kinetic investigation of this important enzyme is that thrombin exists in three forms, E*, E and E:Na+, that interconvert under the influence of ligand binding to distinct domains. The transition between the Na+ -free slow from E and the Na+ -bound fast form E:Na+ involves the structure of the enzyme as a whole, and so does the interconversion between the two Na+ -free forms E* and E. E* is most likely an inactive form of thrombin, unable to interact with Na + and substrate. The complexity of thrombin function and regulation has gained this enzyme pre-eminence as the prototypic allosteric serine protease. Thrombin is now looked upon as a model system for the quantitative analysis of biologically important enzymes.

Conformational Analysis of γ‘ Peptide (410−427) Interactions with Thrombin Anion Binding Exosite II

Thrombin utilizes two anion binding exosites to supplement binding of fibrinogen to this serine protease. Approximately 7-15% of the fibrinogen gamma chain exists as the highly anionic gamma' variant (408VRPEHPAETEY(S)DSLY(S)PEDDL427). This segment has been demonstrated to target thrombin ABE-II and can accommodate sites of phosphorylation in place of sulfonation without sacrificing binding affinity. The present work employed 1D and 2D solution NMR to characterize the structural features of the bound gamma' peptide (410-427) and to evaluate the requirement of sulfonation for effective thrombin interaction. The results indicate the gamma' residues 414-427 make significant contact with the enzyme, a beta-turn exists between residues 422-425 in the presence of thrombin, and there is a large cluster of through-space interactions involving residues 418-422. Effective contact with ABE-II requires the presence of at least one phosphotyrosine residue with Y(P)422 being the more important player. Hydrogen-deuterium exchange (HDX) coupled with MALDI-TOF MS was implemented to examine the location of the gamma' peptide-thrombin interface and to screen for changes in solvent exposure at distant sites. The HDX results demonstrate that the gamma' peptide interacts with or is in close proximity to thrombin residues R93, R97, R173, and R175. The binding of the gamma' peptide also protects other regions of thrombin from deuterium exchange. Affected regions include segments of ABE-I, the autolysis loop, the edge of the active site region, and the A-chain. Finally, thrombin forms a ternary complex with the gamma' peptide and PPACK, generating an enzyme whose solvent-exposed regions are even further stabilized from HDX.