The NMR Structure and Dynamics of the Two-Domain Tick Carboxypeptidase Inhibitor Reveal Flexibility in Its Free Form and Stiffness upon Binding to Human Carboxypeptidase B

Carboxypeptidase inhibitors from Solanaceae as a new subclass of pathogenesis related peptide aiming biotechnological targets for plant defense

Background: Plant protease inhibitors play a crucial role in inhibiting proteases produced by phytopathogens and exhibiting inhibitory effects on nematodes, fungi, and insects, making them promising candidates for crop protection. Specifically, carboxypeptidase inhibitors, a subset of proteinase inhibitors, have been extensively studied in potato and tomato of Solanaceae plant family. However, further research is needed to fully understand the functions and biotechnological potential of those inhibitors in plants. This work aimed to in silico characterize carboxypeptidase inhibitors from Solanaceae as potential antimicrobial and defense agents focused on biotechnological targets. Methods: The methodology employed involved search in UniProt, PDB, KNOTTIN, NCBI, and MEROPS databases for solanaceous carboxypeptidase inhibitors, phylogenetic relationships and conservation patterns analyzes using MEGA-X software and Clustal Omega/MView tools, physicochemical properties and antimicrobial potential prediction using ProtParam, ToxinPred, iAMPred, and APD3 tools, and structural features prediction using PSIPRED. Results and discussion: A systematic literature search was conducted to identify relevant studies on Solanaceae carboxypeptidase inhibitors and their activities against pathogens. The selected studies were reviewed and the main findings compiled. The characterization of Solanaceae carboxypeptidase inhibitors proposed for the first time the global sequence consensus motif CXXXCXXXXDCXXXXXCXXC, shedding light on carboxypeptidase inhibitors distribution, sequence variability, and conservation patterns. Phylogenetic analysis showed evolutionary relationships within the Solanaceae family, particularly in Capsicum, Nicotiana, and Solanum genera. Physicochemical characteristics of those peptides indicated their similarity to antimicrobial peptides. Predicted secondary structures exhibited variations, suggesting a broad spectrum of action, and studies had been demonstrated their activities against various pathogens. Conclusion: Carboxypeptidase inhibitors are being proposed here as a new subclass of PR-6 pathogenesis-related proteins, which will aid in a focused understanding of their functional roles in plant defense mechanisms. These findings confirm the Solanaceae carboxypeptidase inhibitors potential as defense agents and highlight opportunities for their biotechnological applications in pathogen control.

Structural Analysis of the Black-Legged Tick Saliva Protein Salp15

Salp15 is one of the proteins in the saliva of the tick Ixodes scapularis. Together with other biomolecules injected into the mammalian host at the biting site, it helps the tick to sustain its blood meal for days. Salp15 interferes with the cellular immune response of the mammalian host by inhibiting the activation of CD4⁺ T-lymphocytes. This function is co-opted by pathogens that use the tick as a vector and invade the host when the tick bites, such as Borrelia burgdorferi, the causative agent of Lyme borreliosis. Because of the immunity-suppressing role of Salp15, it has been proposed as a candidate for therapeutic applications in disorders of the immune system. The protein is produced as a 135-residue long polypeptide and secreted without its N-terminal signal 1–21 sequence. Detailed structural studies on Salp15 are lacking because of the difficulty in producing large amounts of the folded protein. We report the production of Salp15 and its structural analysis by NMR. The protein is monomeric and contains a flexible N-terminal region followed by a folded domain with mixed α + β secondary structures. Our results are consistent with a three-dimensional structural model derived from AlphaFold, which predicts the formation of three disulfide bridges and a free C-terminal cysteine.

Structure, function, and evolution of Gga-AvBD11, the archetype of the structural avian-double-β-defensin family

Out of the 14 avian β-defensins identified in the Gallus gallus genome, only 3 are present in the chicken egg, including the egg-specific avian β-defensin 11 (Gga-AvBD11). Given its specific localization and its established antibacterial activity, Gga-AvBD11 appears to play a protective role in embryonic development. Gga-AvBD11 is an atypical double-sized defensin, predicted to possess 2 motifs related to β-defensins and 6 disulfide bridges. The 3-dimensional NMR structure of the purified Gga-AvBD11 is a compact fold composed of 2 packed β-defensin domains. This fold is the archetype of a structural family, dubbed herein as avian-double-β-defensins (Av-DBD). We speculate that AvBD11 emanated from a monodomain gene ancestor and that similar events might have occurred in arthropods, leading to another structural family of less compact DBDs. We show that Gga-AvBD11 displays antimicrobial activities against gram-positive and gram-negative bacterial pathogens, the avian protozoan Eimeria tenella, and avian influenza virus. Gga-AvBD11 also shows cytotoxic and antiinvasive activities, suggesting that it may not only be involved in innate protection of the chicken embryo, but also in the (re)modeling of embryonic tissues. Finally, the contribution of either of the 2 Gga-AvBD11 domains to these biological activities was assessed, using chemically synthesized peptides. Our results point to a critical importance of the cationic N-terminal domain in mediating antibacterial, antiparasitic, and antiinvasive activities, with the C-terminal domain potentiating the 2 latter activities. Strikingly, antiviral activity in infected chicken cells, accompanied by marked cytotoxicity, requires the full-length protein.

Pre- and post-docking sampling of conformational changes using ClustENM and HADDOCK for protein-protein and protein-DNA systems

Incorporating the dynamic nature of biomolecules in the modeling of their complexes is a challenge, especially when the extent and direction of the conformational changes taking place upon binding is unknown. Estimating whether the binding of a biomolecule to its partner(s) occurs in a conformational state accessible to its unbound form (“conformational selection”) and/or the binding process induces conformational changes (“induced‐fit”) is another challenge. We propose here a method combining conformational sampling using ClustENM—an elastic network‐based modeling procedure—with docking using HADDOCK, in a framework that incorporates conformational selection and induced‐fit effects upon binding. The extent of the applied deformation is estimated from its energetical costs, inspired from mechanical tensile testing on materials. We applied our pre‐ and post‐docking sampling of conformational changes to the flexible multidomain protein‐protein docking benchmark and a subset of the protein‐DNA docking benchmark. Our ClustENM‐HADDOCK approach produced acceptable to medium quality models in 7/11 and 5/6 cases for the protein‐protein and protein‐DNA complexes, respectively. The conformational selection (sampling prior to docking) has the highest impact on the quality of the docked models for the protein‐protein complexes. The induced‐fit stage of the pipeline (post‐sampling), however, improved the quality of the final models for the protein‐DNA complexes. Compared to previously described strategies to handle conformational changes, ClustENM‐HADDOCK performs better than two‐body docking in protein‐protein cases but worse than a flexible multidomain docking approach. However, it does show a better or similar performance compared to previous protein‐DNA docking approaches, which makes it a suitable alternative.

Peptidase inhibitors in tick physiology

Peptidase inhibitors regulate a wide range of physiological processes involved in the interaction between hematophagous parasites and their hosts, including tissue remodeling, the immune response and blood coagulation. In tick physiology, peptidase inhibitors have a crucial role in adaptation to improve parasitism mechanisms, facilitating blood feeding by interfering with defense-related host peptidases. Recently, a larger number of studies on this topic led to the description of several new tick inhibitors displaying interesting novel features, for example a role in pathogen transmission to the host. A comprehensive review discussing these emerging concepts can therefore shed light on peptidase inhibitor functions, their relevance to tick physiology and their potential applications. Here, we summarize and examine the general characteristics, functional diversity and action of tick peptidase inhibitors with known physiological roles in the tick-host-pathogen interaction.

Design, synthesis, and functional evaluation of leukocyte function associated antigen-1 antagonists in early and late stages of cancer development

The integrin leukocyte function associated antigen 1 (LFA-1) binds the intercellular adhesion molecule 1 (ICAM-1) by its α(L)-chain inserted domain (I-domain). This interaction plays a key role in cancer and other diseases. We report the structure-based design, small-scale synthesis, and biological activity evaluation of a novel family of LFA-1 antagonists. The design led to the synthesis of a family of highly substituted homochiral pyrrolidines with antiproliferative and antimetastatic activity in a murine model of colon carcinoma, as well as potent antiadhesive properties in several cancer cell lines in the low micromolar range. NMR analysis of their binding to the isolated I-domain shows that they bind to the I-domain allosteric site (IDAS), the binding site of other allosteric LFA-1 inhibitors. These results provide evidence of the potential therapeutic value of a new set of LFA-1 inhibitors, whose further development is facilitated by a synthetic strategy that is versatile and fully stereocontrolled.

Proliferating cell nuclear antigen (PCNA) interactions in solution studied by NMR

PCNA is an essential factor for DNA replication and repair. It forms a ring shaped structure of 86 kDa by the symmetric association of three identical protomers. The ring encircles the DNA and acts as a docking platform for other proteins, most of them containing the PCNA Interaction Protein sequence (PIP-box). We have used NMR to characterize the interactions of PCNA with several other proteins and fragments in solution. The binding of the PIP-box peptide of the cell cycle inhibitor p21 to PCNA is consistent with the crystal structure of the complex. A shorter p21 peptide binds with reduced affinity but retains most of the molecular recognition determinants. However the binding of the corresponding peptide of the tumor suppressor ING1 is extremely weak, indicating that slight deviations from the consensus PIP-box sequence dramatically reduce the affinity for PCNA, in contrast with a proposed less stringent PIP-box sequence requirement. We could not detect any binding between PCNA and the MCL-1 or the CDK2 protein, reported to interact with PCNA in biochemical assays. This suggests that they do not bind directly to PCNA, or they do but very weakly, with additional unidentified factors stabilizing the interactions in the cell. Backbone dynamics measurements show three PCNA regions with high relative flexibility, including the interdomain connector loop (IDCL) and the C-terminus, both of them involved in the interaction with the PIP-box. Our work provides the basis for high resolution studies of direct ligand binding to PCNA in solution.

A multidomain flexible docking approach to deal with large conformational changes in the modeling of biomolecular complexes

Binding-induced backbone and large-scale conformational changes represent one of the major challenges in the modeling of biomolecular complexes by docking. To address this challenge, we have developed a flexible multidomain docking protocol that follows a "divide-and-conquer" approach to model both large-scale domain motions and small- to medium-scale interfacial rearrangements: the flexible binding partner is treated as an assembly of subparts/domains that are docked simultaneously making use of HADDOCK's multidomain docking ability. For this, the flexible molecules are cut at hinge regions predicted using an elastic network model. The performance of this approach is demonstrated on a benchmark covering an unprecedented range of conformational changes of 1.5 to 19.5 Å. We show from a statistical survey of known complexes that the cumulative sum of eigenvalues obtained from the elastic network has some predictive power to indicate the extent of the conformational change to be expected.

Reduced stability and increased dynamics in the human proliferating cell nuclear antigen (PCNA) relative to the yeast homolog

Proliferating Cell Nuclear Antigen (PCNA) is an essential factor for DNA replication and repair. PCNA forms a toroidal, ring shaped structure of 90 kDa by the symmetric association of three identical monomers. The ring encircles the DNA and acts as a platform where polymerases and other proteins dock to carry out different DNA metabolic processes. The amino acid sequence of human PCNA is 35% identical to the yeast homolog, and the two proteins have the same 3D crystal structure. In this report, we give evidence that the budding yeast (sc) and human (h) PCNAs have highly similar structures in solution but differ substantially in their stability and dynamics. hPCNA is less resistant to chemical and thermal denaturation and displays lower cooperativity of unfolding as compared to scPCNA. Solvent exchange rates measurements show that the slowest exchanging backbone amides are at the β-sheet, in the structure core, and not at the helices, which line the central channel. However, all the backbone amides of hPCNA exchange fast, becoming undetectable within hours, while the signals from the core amides of scPCNA persist for longer times. The high dynamics of the α-helices, which face the DNA in the PCNA-loaded form, is likely to have functional implications for the sliding of the PCNA ring on the DNA since a large hole with a flexible wall facilitates the establishment of protein-DNA interactions that are transient and easily broken. The increased dynamics of hPCNA relative to scPCNA may allow it to acquire multiple induced conformations upon binding to its substrates enlarging its binding diversity.

Protease inhibitors as models for the study of oxidative folding

The correct balance between proteases and their natural protein inhibitors is of great importance in living systems. Protease inhibitors usually comprise small folds that are crosslinked by a high number of disulfide bonds, making them perfect models for the study of oxidative folding. To date, the oxidative folding of numerous protease inhibitors has been analyzed, revealing a great diversity of folding pathways that differ mainly in the heterogeneity and native disulfide-bond content of their intermediates. The two extremes of this diversity are represented by bovine pancreatic trypsin inhibitor and hirudin, which fold, respectively, via few native intermediates and heterogeneous scrambled isomers. Other proteins, such as leech carboxypeptidase inhibitor, share characteristics of both models displaying mixed folding pathways. The study of the oxidative folding of two-domain inhibitors, such as secretory leukocyte protease inhibitor, tick carboxypeptidase inhibitor, and Ascaris carboxypeptidase inhibitor, has provided some clues about how two-domain protease inhibitors may fold, that is, either by folding each domain autonomously or with one domain assisting in the folding of the other. Finally, the recent determination of the structures of the major intermediates of protease inhibitors has shed light on the molecular mechanisms guiding the oxidative folding of small disulfide-rich proteins.

Catalytic residues in hydrolases: analysis of methods designed for ligand-binding site prediction

The comparison of eight tools applicable to ligand-binding site prediction is presented. The methods examined cover three types of approaches: the geometrical (CASTp, PASS, Pocket-Finder), the physicochemical (Q-SiteFinder, FOD) and the knowledge-based (ConSurf, SuMo, WebFEATURE). The accuracy of predictions was measured in reference to the catalytic residues documented in the Catalytic Site Atlas. The test was performed on a set comprising selected chains of hydrolases. The results were analysed with regard to size, polarity, secondary structure, accessible solvent area of predicted sites as well as parameters commonly used in machine learning (F-measure, MCC). The relative accuracies of predictions are presented in the ROC space, allowing determination of the optimal methods by means of the ROC convex hull. Additionally the minimum expected cost analysis was performed. Both advantages and disadvantages of the eight methods are presented. Characterization of protein chains in respect to the level of difficulty in the active site prediction is introduced. The main reasons for failures are discussed. Overall, the best performance offers SuMo followed by FOD, while Pocket-Finder is the best method among the geometrical approaches.