How does heparin prevent the pH inactivation of cathepsin B? Allosteric mechanism elucidated by docking and molecular dynamics

In silico approaches for better understanding cysteine cathepsin-glycosaminoglycan interactions

Cysteine cathepsins constitute the largest cathepsin family, with 11 proteases in human that are present primarily within acidic endosomal and lysosomal compartments. They are involved in the turnover of intracellular and extracellular proteins. They are synthesized as inactive procathepsins that are converted to mature active forms. Cathepsins play important roles in physiological and pathological processes and, therefore, receive increasing attention as potential therapeutic targets. Their maturation and activity can be regulated by glycosaminoglycans (GAGs), long linear negatively charged polysaccharides composed of recurring dimeric units. In this review, we summarize recent computational progress in the field of (pro)cathepsin-GAG complexes analyses.

Characterization of a novel dipeptidyl peptidase 1 of Trichinella spiralis and its participation in larval invasion

The research aimed to describe a new Trichinella spiralis dipeptidyl peptidase 1 (TsDPP1) and investigate its functions in the larval invasion of intestinal epithelial cells (IECs). The gene TsDPP1 was successfully replicated and produced in Escherichia coli BL21 (DE3), showing a strong immune response. TsDPP1 was detected in diverse stages of T. spiralis and showed significant expression in the intestine infective larvae (IIL) and adult worms at 6 days post infection, as confirmed by qPCR and Western blot analysis. The primary localization of TsDPP1 in this parasite was observed in cuticles, stichosomes, and embryos by using the indirect immunofluorescence assay (IIFA). rTsDPP1 exhibited the enzymatic function of natural dipeptidyl peptidase and showed specific binding to IECs, and the binding site was found to be localized on cell membrane. Following transfection with dsRNA-TsDPP1, the expression of TsDPP1 mRNA and protein in muscle larvae (ML) were decreased by approximately 63.52 % and 58.68 %, correspondingly. The activity of TsDPP1 in the ML and IIL treated with dsRNA-TsDPP1 was reduced by 42.98 % and 45.07 %, respectively. The acceleration of larval invasion of IECs was observed with rTsDPP1, while the invasion was suppressed by anti-rTsDPP1 serum. The ability of the larvae treated with dsRNA-TsDPP1 to invade IECs was hindered by 31.23 %. In mice infected with dsRNA-treated ML, the intestinal IIL, and adults experienced a significant decrease in worm burdens and a noticeable reduction in adult female length and fecundity compared to the PBS group. These findings indicated that TsDPP1 significantly impedes the invasion, growth, and reproductive capacity of T. spiralis in intestines, suggesting its potential as a target for anti-Trichinella vaccines.

Microglial P2X4 receptors promote ApoE degradation and contribute to memory deficits in Alzheimer's disease

Numerous evidences support that microglia contributes to the progression of Alzheimer's disease. P2X4 receptors are ATP-gated channels with high calcium permeability, which are de novo expressed in a subset of reactive microglia associated with various pathological contexts, contributing to microglial functions. P2X4 receptors are mainly localized in lysosomes and trafficking to the plasma membrane is tightly regulated. Here, we investigated the role of P2X4 in the context of Alzheimer's disease (AD). Using proteomics, we identified Apolipoprotein E (ApoE) as a specific P2X4 interacting protein. We found that P2X4 regulates lysosomal cathepsin B (CatB) activity promoting ApoE degradation; P2rX4 deletion results in higher amounts of intracellular and secreted ApoE in both bone-marrow-derived macrophage (BMDM) and microglia from APPswe/PSEN1dE9 brain. In both human AD brain and APP/PS1 mice, P2X4 and ApoE are almost exclusively expressed in plaque-associated microglia. In 12-month-old APP/PS1 mice, genetic deletion of P2rX4 reverses topographical and spatial memory impairment and reduces amount of soluble small aggregates of Aß1-42 peptide, while no obvious alteration of plaque-associated microglia characteristics is observed. Our results support that microglial P2X4 promotes lysosomal ApoE degradation, indirectly altering Aß peptide clearance, which in turn might promotes synaptic dysfunctions and cognitive deficits. Our findings uncover a specific interplay between purinergic signaling, microglial ApoE, soluble Aß (sAß) species and cognitive deficits associated with AD.

Characterization of Heparin's Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin's solution conformational dynamics and structure-function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule's underlying conformational ensemble, which then provides important information necessary for understanding structure-function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin's conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. The anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

Leveraging 3D Model Systems to Understand Viral Interactions with the Respiratory Mucosa

Respiratory viruses remain a significant cause of morbidity and mortality in the human population, underscoring the importance of ongoing basic research into virus-host interactions. However, many critical aspects of infection are difficult, if not impossible, to probe using standard cell lines, 2D culture formats, or even animal models. In vitro systems such as airway epithelial cultures at air-liquid interface, organoids, or 'on-chip' technologies allow interrogation in human cells and recapitulate emergent properties of the airway epithelium-the primary target for respiratory virus infection. While some of these models have been used for over thirty years, ongoing advancements in both culture techniques and analytical tools continue to provide new opportunities to investigate airway epithelial biology and viral infection phenotypes in both normal and diseased host backgrounds. Here we review these models and their application to studying respiratory viruses. Furthermore, given the ability of these systems to recapitulate the extracellular microenvironment, we evaluate their potential to serve as a platform for studies specifically addressing viral interactions at the mucosal surface and detail techniques that can be employed to expand our understanding.

Cysteine-type cathepsins promote the effector phase of acute cutaneous delayed-type hypersensitivity reactions

Cysteine-type cathepsins such as cathepsin B are involved in various steps of inflammatory processes such as antigen processing and angiogenesis. Here, we uncovered the role of cysteine-type cathepsins in the effector phase of T cell-driven cutaneous delayed-type hypersensitivity reactions (DTHR) and the implication of this role on therapeutic cathepsin B-specific inhibition.

Methods: Wild-type, cathepsin B-deficient (Ctsb^-/-) and cathepsin Z-deficient (Ctsz^-/-) mice were sensitized with 2,4,6-trinitrochlorobenzene (TNCB) on the abdomen and challenged with TNCB on the right ear to induce acute and chronic cutaneous DTHR. The severity of cutaneous DTHR was assessed by evaluating ear swelling responses and histopathology. We performed fluorescence microscopy on tissue from inflamed ears and lymph nodes of wild-type mice, as well as on biopsies from psoriasis patients, focusing on cathepsin B expression by T cells, B cells, macrophages, dendritic cells and NK cells. Cathepsin activity was determined noninvasively by optical imaging employing protease-activated substrate-like probes. Cathepsin expression and activity were validated ex vivo by covalent active site labeling of proteases and Western blotting.

Results: Noninvasive in vivo optical imaging revealed strong cysteine-type cathepsin activity in inflamed ears and draining lymph nodes in acute and chronic cutaneous DTHR. In inflamed ears and draining lymph nodes, cathepsin B was expressed by neutrophils, dendritic cells, macrophages, B, T and natural killer (NK) cells. Similar expression patterns were found in psoriatic plaques of patients. The biochemical methods confirmed active cathepsin B in tissues of mice with cutaneous DTHR. Topically applied cathepsin B inhibitors significantly reduced ear swelling in acute but not chronic DTHR. Compared with wild-type mice, Ctsb^-/- mice exhibited an enhanced ear swelling response during acute DTHR despite a lack of cathepsin B expression. Cathepsin Z, a protease closely related to cathepsin B, revealed compensatory expression in inflamed ears of Ctsb^-/- mice, while cathepsin B expression was reciprocally elevated in Ctsz^-/- mice.

Conclusion: Cathepsin B is actively involved in the effector phase of acute cutaneous DTHR. Thus, topically applied cathepsin B inhibitors might effectively limit DTHR such as contact dermatitis or psoriasis. However, the cathepsin B and Z knockout mouse experiments suggested a complementary role for these two cysteine-type proteases.

Dissecting a novel allosteric mechanism of cruzain: A computer-aided approach

Trypanosoma cruzi is the causative agent of Chagas disease, a neglected infection affecting millions of people in tropical regions. There are several chemotherapeutic agents for the treatment of this disease, but most of them are highly toxic and generate resistance. Currently, the development of allosteric inhibitors constitutes a promising research field, since it can improve the accessibility to more selective and less toxic medicines. To date, the allosteric drugs prediction is a state-of-the-art topic in rational structure-based computational design. In this work, a simulation strategy was developed for computational discovery of allosteric inhibitors, and it was applied to cruzain, a promising target and the major cysteine protease of T. cruzi. Molecular dynamics simulations, binding free energy calculations and network-based modelling of residue interactions were combined to characterize and compare molecular distinctive features of the apo form and the cruzain-allosteric inhibitor complexes. By using geometry-based criteria on trajectory snapshots, we predicted two main allosteric sites suitable for drug targeting. The results suggest dissimilar mechanisms exerted by the same allosteric site when binding different potential allosteric inhibitors. Finally, we identified the residues involved in suboptimal paths linking the identified site and the orthosteric site. The present study constitutes the first approximation to the design of cruzain allosteric inhibitors and may serve for future pharmacological intervention. Here, no major effects on active site structure were observed due to compound binding (modification of distance and angles between catalytic residues), which indicates that allosteric regulation in cruzain might be mediated via alterations of its dynamical properties similarly to allosteric inhibition of human cathepsin K (HCatK). The current findings are particularly relevant for the design of allosteric modulators of papain-like cysteine proteases.

Cathepsins: Proteases that are vital for survival but can also be fatal

The state of enzymes in the human body determines the normal physiology or pathology, so all the six classes of enzymes are crucial. Proteases, the hydrolases, can be of several types based on the nucleophilic amino acid or the metal cofactor needed for their activity. Cathepsins are proteases with serine, cysteine, or aspartic acid residues as the nucleophiles, which are vital for digestion, coagulation, immune response, adipogenesis, hormone liberation, peptide synthesis, among a litany of other functions. But inflammatory state radically affects their normal roles. Released from the lysosomes, they degrade extracellular matrix proteins such as collagen and elastin, mediating parasite infection, autoimmune diseases, tumor metastasis, cardiovascular issues, and neural degeneration, among other health hazards. Over the years, the different types and isoforms of cathepsin, their optimal pH and functions have been studied, yet much information is still elusive. By taming and harnessing cathepsins, by inhibitors and judicious lifestyle, a gamut of malignancies can be resolved. This review discusses these aspects, which can be of clinical relevance.

A novel allosteric mechanism in the cysteine peptidase cathepsin K discovered by computational methods

Allosteric modifiers have the potential to fine-tune enzyme activity. Therefore, targeting allosteric sites is gaining increasing recognition as a strategy in drug design. Here we report the use of computational methods for the discovery of the first small-molecule allosteric inhibitor of the collagenolytic cysteine peptidase cathepsin K, a major target for the treatment of osteoporosis. The molecule NSC13345 is identified by high-throughput docking of compound libraries to surface sites on the peptidase that are connected to the active site by an evolutionarily conserved network of residues (protein sector). The crystal structure of the complex shows that NSC13345 binds to a novel allosteric site on cathepsin K. The compound acts as a hyperbolic mixed modifier in the presence of a synthetic substrate, it completely inhibits collagen degradation and has good selectivity for cathepsin K over related enzymes. Altogether, these properties qualify our methodology and NSC13345 as promising candidates for allosteric drug design.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

Towards the computational design of protein post-translational regulation

Protein post-translational modifications (PTMs) are a fast and versatility mechanism used by the cell to regulate the function of proteins in response to changing conditions. PTMs can alter the activity of proteins by allosteric regulation or by controlling protein interactions, localization and abundance. Recent advances in proteomics have revealed the extent of regulation by PTMs and the different mechanisms used in nature to exert control over protein function via PTMs. These developments can serve as the foundation for the rational design of protein regulation. Here we review the advances in methods to determine the function of PTMs, protein allosteric control and examples of rational design of PTM regulation. These advances create an opportunity to move synthetic biology forward by making use of a level of regulation that is of yet unexplored.

Cysteine cathepsin activity regulation by glycosaminoglycans

Cysteine cathepsins are a group of enzymes normally found in the endolysosomes where they are primarily involved in intracellular protein turnover but also have a critical role in MHC II-mediated antigen processing and presentation. However, in a number of pathologies cysteine cathepsins were found to be heavily upregulated and secreted into extracellular milieu, where they were found to degrade a number of extracellular proteins. A major role in modulating cathepsin activities play glycosaminoglycans, which were found not only to facilitate their autocatalytic activation including at neutral pH, but also to critically modulate their activities such as in the case of the collagenolytic activity of cathepsin K. The interaction between cathepsins and glycosaminoglycans will be discussed in more detail.

Impact of M36I polymorphism on the interaction of HIV-1 protease with its substrates: insights from molecular dynamics

Background

Over the last decades, a vast structural knowledge has been gathered on the HIV-1 protease (PR). Noticeably, most of the studies focused the B-subtype, which has the highest prevalence in developed countries. Accordingly, currently available anti-HIV drugs target this subtype, with considerable benefits for the corresponding patients.

However, in developing countries, there is a wide variety of HIV-1 subtypes carrying PR polymorphisms related to reduced drug susceptibility. The non-active site mutation, M36I, is the most frequent polymorphism, and is considered as a non-B subtype marker.

Yet, the structural impact of this substitution on the PR structure and on the interaction with natural substrates remains poorly documented.

Results

Herein, we used molecular dynamics simulations to investigate the role of this polymorphism on the interaction of PR with six of its natural cleavage-sites substrates.

Free energy analyses by MMPB/SA calculations showed an affinity decrease of M36I-PR for the majority of its substrates. The only exceptions were the RT-RH, with equivalent affinity, and the RH-IN, for which an increased affinity was found. Furthermore, molecular simulations suggest that, unlike other peptides, RH-IN induced larger structural fluctuations in the wild-type enzyme than in the M36I variant.

Conclusions

With multiple approaches and analyses we identified structural and dynamical determinants associated with the changes found in the binding affinity of the M36I variant. This mutation influences the flexibility of both PR and its complexed substrate. The observed impact of M36I, suggest that combination with other non-B subtype polymorphisms, could lead to major effects on the interaction with the 12 known cleavage sites, which should impact the virion maturation.

Gene-centric meta-analyses for central adiposity traits in up to 57 412 individuals of European descent confirm known loci and reveal several novel associations

Waist circumference (WC) and waist-to-hip ratio (WHR) are surrogate measures of central adiposity that are associated with adverse cardiovascular events, type 2 diabetes and cancer independent of body mass index (BMI). WC and WHR are highly heritable with multiple susceptibility loci identified to date. We assessed the association between SNPs and BMI-adjusted WC and WHR and unadjusted WC in up to 57 412 individuals of European descent from 22 cohorts collaborating with the NHLBI's Candidate Gene Association Resource (CARe) project. The study population consisted of women and men aged 20-80 years. Study participants were genotyped using the ITMAT/Broad/CARE array, which includes ∼50 000 cosmopolitan tagged SNPs across ∼2100 cardiovascular-related genes. Each trait was modeled as a function of age, study site and principal components to control for population stratification, and we conducted a fixed-effects meta-analysis. No new loci for WC were observed. For WHR analyses, three novel loci were significantly associated (P < 2.4 × 10(-6)). Previously unreported rs2811337-G near TMCC1 was associated with increased WHR (β ± SE, 0.048 ± 0.008, P = 7.7 × 10(-9)) as was rs7302703-G in HOXC10 (β = 0.044 ± 0.008, P = 2.9 × 10(-7)) and rs936108-C in PEMT (β = 0.035 ± 0.007, P = 1.9 × 10(-6)). Sex-stratified analyses revealed two additional novel signals among females only, rs12076073-A in SHC1 (β = 0.10 ± 0.02, P = 1.9 × 10(-6)) and rs1037575-A in ATBDB4 (β = 0.046 ± 0.01, P = 2.2 × 10(-6)), supporting an already established sexual dimorphism of central adiposity-related genetic variants. Functional analysis using ENCODE and eQTL databases revealed that several of these loci are in regulatory regions or regions with differential expression in adipose tissue.

Modeling the zing finger protein SmZF1 from Schistosoma mansoni: Insights into DNA binding and gene regulation

Zinc finger proteins are widely found in eukaryotes, representing an important class of DNA-binding proteins frequently involved in transcriptional regulation. Zinc finger motifs are composed by two antiparallel β-strands and one α-helix, stabilized by a zinc ion coordinated by conserved histidine and cysteine residues. In Schistosoma mansoni, these regulatory proteins are known to modulate morphological and physiological changes, having crucial roles in parasite development. A previously described C(2)H(2) zinc finger protein, SmZF1, was shown to be present in cell nuclei of different life stages of S. mansoni and to activate gene transcription in a heterologous system. A high-quality SmZF1 tridimensional structure was generated using comparative modeling. Molecular dynamics simulations of the obtained structure revealed stability of the zinc fingers motifs and high flexibility on the terminals, comparable to the profile observed on the template X-ray structure based on thermal b-factors. Based on the protein tridimensional features and amino acid composition, we were able to characterize four C(2)H(2) zinc finger motifs, the first involved in protein-protein interactions while the three others involved in DNA binding. We defined a consensus DNA binding sequence using three distinct algorithms and further carried out docking calculations, which revealed the interaction of fingers 2-4 with the predicted DNA. A search for S. mansoni genes presenting putative SmZF1 binding sites revealed 415 genes hypothetically under SmZF1 control. Using an automatic annotation and GO assignment approach, we found that the majority of those genes code for proteins involved in developmental processes. Taken together, these results present a consistent base to the structural and functional characterization of SmZF1.

Cysteine cathepsins: from structure, function and regulation to new frontiers

It is more than 50 years since the lysosome was discovered. Since then its hydrolytic machinery, including proteases and other hydrolases, has been fairly well identified and characterized. Among these are the cysteine cathepsins, members of the family of papain-like cysteine proteases. They have unique reactive-site properties and an uneven tissue-specific expression pattern. In living organisms their activity is a delicate balance of expression, targeting, zymogen activation, inhibition by protein inhibitors and degradation. The specificity of their substrate binding sites, small-molecule inhibitor repertoire and crystal structures are providing new tools for research and development. Their unique reactive-site properties have made it possible to confine the targets simply by the use of appropriate reactive groups. The epoxysuccinyls still dominate the field, but now nitriles seem to be the most appropriate “warhead”. The view of cysteine cathepsins as lysosomal proteases is changing as there is now clear evidence of their localization in other cellular compartments. Besides being involved in protein turnover, they build an important part of the endosomal antigen presentation. Together with the growing number of non-endosomal roles of cysteine cathepsins is growing also the knowledge of their involvement in diseases such as cancer and rheumatoid arthritis, among others. Finally, cysteine cathepsins are important regulators and signaling molecules of an unimaginable number of biological processes. The current challenge is to identify their endogenous substrates, in order to gain an insight into the mechanisms of substrate degradation and processing. In this review, some of the remarkable advances that have taken place in the past decade are presented. This article is part of a Special Issue entitled: Proteolysis 50 years after the discovery of lysosome.

Free Energy Profiles along Consensus Normal Modes Provide Insight into HIV-1 Protease Flap Opening

Describing biological macromolecular energetics from computer simulations can pose major challenges, and often necessitates enhanced conformational sampling. We describe the calculation of conformational free-energy profiles along carefully chosen collective coordinates: "consensus" normal modes, developed recently as robust alternatives to conventional normal modes. In an application to the HIV-1 protease, we obtain efficient sampling of significant flap opening movements governing inhibitor binding from relatively short simulations, in close correspondence with experimental results.