Two-Stage Recognition Mechanism of the SARS-CoV-2 Receptor-Binding Domain to Angiotensin-Converting Enzyme-2 (ACE2)

The SARS-CoV-2 virus, commonly known as COVID-19, occurred in 2019. It is a highly contagious illness with effects ranging from mild symptoms to severe illness. It is also one of the best-known pathogens since more than 200,000 scientific papers occurred in the last few years. With the publication of the SARS-CoV-2 (SARS-CoV-2-CTD) spike (S) protein in a complex with human ACE2 (hACE2) (PDB (6LZG)), the molecular analysis of one of the most crucial steps on the infection pathway was possible. The aim of this manuscript is to simulate the most widely spread mutants of SARS-CoV-2, namely Alpha, Beta, Gamma, Delta, Omicron, and the first recognized variant (natural wild type). With the wide search of the hypersurface of the potential energy performed using the UNRES force field, the intermediate state of the ACE2-RBD complex was found. R403, K/N/T417, L455, F486, Y489, F495, Y501, and Y505 played a crucial role in the protein recognition mechanism. The intermediate state cannot be very stable since it will prevent the infection cascade.

Impact of AlphaFold on structure prediction of protein complexes: The CASP15-CAPRI experiment

We present the results for CAPRI Round 54, the 5th joint CASP-CAPRI protein assembly prediction challenge. The Round offered 37 targets, including 14 homodimers, 3 homo-trimers, 13 heterodimers including 3 antibody-antigen complexes, and 7 large assemblies. On average ~70 CASP and CAPRI predictor groups, including more than 20 automatics servers, submitted models for each target. A total of 21 941 models submitted by these groups and by 15 CAPRI scorer groups were evaluated using the CAPRI model quality measures and the DockQ score consolidating these measures. The prediction performance was quantified by a weighted score based on the number of models of acceptable quality or higher submitted by each group among their five best models. Results show substantial progress achieved across a significant fraction of the 60+ participating groups. High-quality models were produced for about 40% of the targets compared to 8% two years earlier. This remarkable improvement is due to the wide use of the AlphaFold2 and AlphaFold2-Multimer software and the confidence metrics they provide. Notably, expanded sampling of candidate solutions by manipulating these deep learning inference engines, enriching multiple sequence alignments, or integration of advanced modeling tools, enabled top performing groups to exceed the performance of a standard AlphaFold2-Multimer version used as a yard stick. This notwithstanding, performance remained poor for complexes with antibodies and nanobodies, where evolutionary relationships between the binding partners are lacking, and for complexes featuring conformational flexibility, clearly indicating that the prediction of protein complexes remains a challenging problem.

A Novel Cryptic Clostridial Peptide That Kills Bacteria by a Cell Membrane Permeabilization Mechanism

This work reports detailed characteristics of the antimicrobial peptide Intestinalin (P30), which is derived from the LysC enzyme of Clostridium intestinale strain URNW. The peptide shows a broader antibacterial spectrum than the parental enzyme, showing potent antimicrobial activity against clinical strains of Gram-positive staphylococci and Gram-negative pathogens and causing between 3.04 ± 0.12 log kill for Pseudomonas aeruginosa PAO1 and 7.10 ± 0.05 log kill for multidrug-resistant Acinetobacter baumannii KPD 581 at a 5 μM concentration. Moreover, Intestinalin (P30) prevents biofilm formation and destroys 24-h and 72-h biofilms formed by Acinetobacter baumannii CRAB KPD 205 (reduction levels of 4.28 and 2.62 log CFU/mL, respectively). The activity of Intestinalin is combined with both no cytotoxicity and little hemolytic effect against mammalian cells. The nuclear magnetic resonance and molecular dynamics (MD) data show a high tendency of Intestinalin to interact with the bacterial phospholipid cell membrane. Although positively charged, Intestinalin resides in the membrane and aggregates into small oligomers. Negatively charged phospholipids stabilize peptide oligomers to form water- and ion-permeable pores, disrupting the integrity of bacterial cell membranes. Experimental data showed that Intestinalin interacts with negatively charged lipoteichoic acid (logK based on isothermal titration calorimetry, 7.45 ± 0.44), causes membrane depolarization, and affects membrane integrity by forming large pores, all of which result in loss of bacterial viability. IMPORTANCE Antibiotic resistance is rising rapidly among pathogenic bacteria, becoming a global public health problem that threatens the effectiveness of therapies for many infectious diseases. In this respect, antimicrobial peptides appear to be an interesting alternative to combat bacterial pathogens. Here, we report the characteristics of an antimicrobial peptide (of 30 amino acids) derived from the clostridial LysC enzyme. The peptide showed killing activity against clinical strains of Gram-positive and Gram-negative pathogens. Experimental data and computational modeling showed that this peptide forms transmembrane pores, directly engaging the negatively charged phospholipids of the bacterial cell membrane. Consequently, dissipation of the electrochemical gradient across cell membranes affects many vital processes, such as ATP synthesis, motility, and transport of nutrients. This kind of dysfunction leads to the loss of bacterial viability. Our firm conviction is that the presented study will be a helpful resource in searching for novel antimicrobial peptides that could have the potential to replace conventional antibiotics.

Cathepsin C inhibition as a potential treatment strategy in cancer

Epidemiological studies established an association between chronic inflammation and higher risk of cancer. Inhibition of proteolytic enzymes represents a potential treatment strategy for cancer and prevention of cancer metastasis. Cathepsin C (CatC) is a highly conserved lysosomal cysteine dipeptidyl aminopeptidase required for the activation of pro-inflammatory neutrophil serine proteases (NSPs, elastase, proteinase 3, cathepsin G and NSP-4). NSPs are locally released by activated neutrophils in response to pathogens and non-infectious danger signals. Activated neutrophils also release neutrophil extracellular traps (NETs) that are decorated with several neutrophil proteins, including NSPs. NSPs are not only NETs constituents but also play a role in NET formation and release. Although immune cells harbor large amounts of CatC, additional cell sources for this protease exists. Upregulation of CatC expression was observed in different tissues during carcinogenesis and correlated with metastasis and poor patient survival. Recent mechanistic studies indicated an important interaction of tumor-associated CatC, NSPs, and NETs in cancer development and metastasis and suggested CatC as a therapeutic target in a several cancer types. Cancer cell-derived CatC promotes neutrophil recruitment in the inflammatory tumor microenvironment. Because the clinical consequences of genetic CatC deficiency in humans resulting in the elimination of NSPs are mild, small molecule inhibitors of CatC are assumed as safe drugs to reduce the NSP burden. Brensocatib, a nitrile CatC inhibitor is currently tested in a phase 3 clinical trial as a novel anti-inflammatory therapy for patients with bronchiectasis. However, recently developed CatC inhibitors possibly have protective effects beyond inflammation. In this review, we describe the pathophysiological function of CatC and discuss molecular mechanisms substantiating pharmacological CatC inhibition as a potential strategy for cancer treatment.

Structure-based design and in vivo anti-arthritic activity evaluation of a potent dipeptidyl cyclopropyl nitrile inhibitor of cathepsin C

Cathepsin C (CatC) is a dipeptidyl-exopeptidase which activates neutrophil serine protease precursors (elastase, proteinase 3, cathepsin G and NSP4) by removing their N-terminal propeptide in bone marrow cells at the promyelocytic stage of neutrophil differentiation. The resulting active proteases are implicated in chronic inflammatory and autoimmune diseases. Hence, inhibition of CatC represents a therapeutic strategy to suppress excessive protease activities in various neutrophil mediated diseases. We designed and synthesized a series of dipeptidyl cyclopropyl nitrile compounds as putative CatC inhibitors. One compound, IcatC_XPZ-01 ((S)-2-amino-N-((1R,2R)-1-cyano-2-(4'-(4-methylpiperazin-1-ylsulfonyl)biphenyl-4-yl)cyclopropyl)butanamide)) was identified as a potent inhibitor of both human and rodent CatC. In mice, pharmacokinetic studies revealed that IcatC_XPZ-01 accumulated in the bone marrow reaching levels suitable for CatC inhibition. Subcutaneous administration of IcatC_XPZ-01 in a monoclonal anti-collagen antibody induced mouse model of rheumatoid arthritis resulted in statistically significant anti-arthritic activity with persistent decrease in arthritis scores and paw thickness.

Biochemical properties of the HtrA homolog from bacterium Stenotrophomonas maltophilia

The HtrA proteins due to their proteolytic, and in many cases chaperone activity, efficiently counteract consequences of stressful conditions. In the environmental bacterium and nosocomial pathogen Stenotrophomonas maltophilia HtrA (HtrA_Sm) is induced as a part of adaptive response to host temperature (37°C). We examined the biochemical properties of HtrA_Sm and compared them with those of model HtrA_Ec from Escherichia coli. We found that HtrA_Sm is a protease and chaperone that operates over a wide range of pH and is highly active at temperatures between 35 and 37°C. The temperature-sensitive activity corresponded well with the lower thermal stability of the protein and weaker stability of the oligomer. Interestingly, the enzyme shows slightly different substrate cleavage specificity when compared to other bacterial HtrAs. A computational model of the three-dimensional structure of HtrA_Sm indicates differences in the S1 substrate specificity pocket and suggests weaker inter-trimer interactions when compared to HtrA_Ec. The observed features of HtrA_Sm suggest that this protein may play a protective role under stressful conditions acting both as a protease and a chaperone. The optimal temperatures for the protein activity may reflect the evolutionary adaptation of S. maltophilia to life in soil or aqueous environments, where the temperatures are usually much below 37°C.

Structural insights into the activation mechanisms of human HtrA serine proteases

Human HtrA1-4 proteins belong to the HtrA family of evolutionarily conserved serine proteases and function as important modulators of many physiological processes, including maintenance of mitochondrial homeostasis, cell signaling and apoptosis. Disturbances in their action are linked to severe diseases, including oncogenesis and neurodegeneration. The HtrA1-4 proteins share structural and functional features of other members of the HtrA protein family, however there are several significant differences in structural architecture and mechanisms of action which makes each of them unique. Our goal is to present recent studies regarding human HtrAs. We focus on their physiological functions, structure and regulation, and describe current models of activation mechanisms. Knowledge of molecular basis of the human HtrAs' action is a subject of great interest; it is crucial for understanding their relevance in cellular physiology and pathogenesis as well as for using them as targets in future therapies of diseases such as neurodegenerative disorders and cancer.

The LD loop as an important structural element required for transmission of the allosteric signal in the HtrA (DegP) protease from Escherichia coli

High-temperature requirement A (HtrA; DegP) from Escherichia coli, an important element of the extracytoplasmic protein quality-control system, is a member of the evolutionarily conserved family of serine proteases. The characteristic feature of this protein is its allosteric mode of activation. The regulatory loops, L3, L2, L1 and LD, play a crucial role in the transmission of the allosteric signal. Yet, the role of LD has not been fully elucidated. Therefore, we undertook a study to explain the role of the individual LD residues in inducing and maintaining the proteolytic activity of HtrA. We investigated the influence of amino acid substitutions located within the LD loop on the kinetics of a model substrate cleavage as well as on the dynamics of the oligomeric structure of HtrA. We found that the mutations that were expected to disturb the loop's structure and/or interactions with the remaining regulatory loops severely diminished the proteolytic activity of HtrA. The opposite effect, that is, increased activity, was observed for G174S substitution, which was predicted to strengthen the interactions mediated by LD. HtrAG174S protein had an equilibrium shifted toward the active enzyme and formed preferentially high-order oligomeric forms.

RASMOL AB - new functionalities in the program for structure analysis

For many years RasMol was one of the most used programs for molecular visualization. It was an excellent tool due to its simplicity and its low demand of computer power. Today it is replaced by OpenGL programs, which have excellent graphics that new computers can additionally handle. Molecular graphics is one of the best tools for the analysis of biomolecular data. With high efficiency and a low demand of computer power, RasMol can still be used as a quick and handy tool used for the analysis of biomolecular structures with good results. In this paper, we describe modifications to the RasMol program, as implemented on the base of RasMol AB 2. We introduced several new functions, namely: the identification of histidine isomers, and advanced structural selection and macro capabilities (as implemented in the point-click menu), which result in an increase in the speed and accuracy of structural analyses. The program can be downloaded from the project page: http://etoh.chem.univ.gda.pl/rasmol/.

The LA loop as an important regulatory element of the HtrA (DegP) protease from Escherichia coli: structural and functional studies

Bacterial HtrAs are serine proteases engaged in extracytoplasmic protein quality control and are required for the virulence of several pathogenic species. The proteolytic activity of HtrA (DegP) from Escherichia coli, a model prokaryotic HtrA, is stimulated by stressful conditions; the regulation of this process is mediated by the LA, LD, L1, L2, and L3 loops. The precise mechanism of action of the LA loop is not known due to a lack of data concerning its three-dimensional structure as well as its mode of interaction with other regulatory elements. To address these issues we generated a theoretical model of the three-dimensional structure of the LA loop as per the resting state of HtrA and subsequently verified its correctness experimentally. We identified intra- and intersubunit contacts that formed with the LA loops; these played an important role in maintaining HtrA in its inactive conformation. The most significant proved to be the hydrophobic interactions connecting the LA loops of the hexamer and polar contacts between the LA' (the LA loop on an opposite subunit) and L1 loops on opposite subunits. Disturbance of these interactions caused the stimulation of HtrA proteolytic activity. We also demonstrated that LA loops contribute to the preservation of the integrity of the HtrA oligomer and to the stability of the monomer. The model presented in this work explains the regulatory role of the LA loop well; it should also be applicable to numerous Enterobacteriaceae pathogenic species as the amino acid sequences of the members of this bacterial family are highly conserved.

Temperature-induced changes of HtrA2(Omi) protease activity and structure

HtrA2(Omi), belonging to the high-temperature requirement A (HtrA) family of stress proteins, is involved in the maintenance of mitochondrial homeostasis and in the stimulation of apoptosis, as well as in cancer and neurodegenerative disorders. The protein comprises a serine protease domain and a postsynaptic density of 95 kDa, disk large, and zonula occludens 1 (PDZ) regulatory domain and functions both as a protease and a chaperone. Based on the crystal structure of the HtrA2 inactive trimer, it has been proposed that PDZ domains restrict substrate access to the protease domain and that during protease activation there is a significant conformational change at the PDZ-protease interface, which removes the inhibitory effect of PDZ from the active site. The crystal structure of the HtrA2 active form is not available yet. HtrA2 activity markedly increases with temperature. To understand the molecular basis of this increase in activity, we monitored the temperature-induced structural changes using a set of single-Trp HtrA2 mutants with Trps located at the PDZ-protease interface. The accessibility of each Trp to aqueous medium was assessed by fluorescence quenching, and these results, in combination with mean fluorescence lifetimes and wavelength emission maxima, indicate that upon an increase in temperature the HtrA2 structure relaxes, the PDZ-protease interface becomes more exposed to the solvent, and significant conformational changes involving both domains occur at and above 30 °C. This conclusion correlates well with temperature-dependent changes of HtrA2 proteolytic activity and the effect of amino acid substitutions (V226K and R432L) located at the domain interface, on HtrA2 activity. Our results experimentally support the model of HtrA2 activation and provide an insight into the mechanism of temperature-induced changes in HtrA2 structure.

Molecular dynamics simulations of the growth of poly(chloro-para-xylylene) films

Parylene C, poly(chloro-para-xylylene) is the most widely used member of the parylene family due to its excellent chemical and physical properties. In this work we analyzed the formation of the parylene C film using molecular mechanics and molecular dynamics methods. A five unit chain is necessary to create a stable hydrophobic cluster and to adhere to a covered surface. Two scenarios were deemed to take place. The obtained results are consistent with a polymer film scaling growth mechanism and contribute to the description of the dynamic growth of the parylene C polymer.

Theoretical study of the human bradykinin-bradykinin B2 receptor complex

The interaction of bradykinin (BK) with the bradykinin B2 receptor (B2R) was analyzed by using molecular modeling (MM) and molecular dynamics (MD) simulations. A homology model for B2R has been generated and the recently determined receptor-bound solid-state NMR spectroscopic structure of BK (Lopez et al., Angew. Chem. 2008, 120, 1692-1695; Angew. Chem. Int. Ed. 2008, 47, 1668-1671) has been modeled into the binding pocket of the receptor to probe the putative ligand-receptor interface. The experimental hormone structure fitted well into the binding pocket of the receptor model and remained stable during the MD simulation. We propose a parallel orientation of the side chains for Arg1 and Arg9 in BK that is bound to B2R. The MD simulation study also allows the conformational changes that lead to the activated form of B2R to be analyzed. The hydrogen bond between N140 (3.35) and W283 (6.48) is the key interaction that keeps the receptor in its inactive form. This hydrogen bond is broken during the MD simulation due to rotation of transmembrane helix 3 (TM3) and is replaced by a new hydrogen bond between W283 (6.48) and N324 (7.45). We propose that this interaction is specific for the activated form of the bradykinin B2 receptor. Additionally, we compared and discussed our putative model in the context of the structural model of the partially activated rhodopsin (Rh*) and with the known biochemical and structural data.

Theoretical study on binding of S100B protein

S100B protein is one of the factors involved in the down-regulation of tumor suppressor protein p53, a transcription activator that signals for cycle arrest and apoptosis. As the inactivation of normal p53 functions is found in over half of human cancers, restoration of normal p53 functions through the destruction or prevention of S100B--p53 complexes represents a possible approach for the development of anti-cancer drugs. The aim of this work was to propose the S100B binding interface through an examination of the literature and use of molecular modeling (MM) techniques with AutoDock program and the AMBER force field. We propose two residues in the S100B binding pocket (Val56, Phe76) and two residues on the protein surface (Val52, Ala83) are essential for ligand binding. The data presented here indicate that interactions with these four residues are necessary for a reduction in the incidence of the S100B--p53 complex. Additionally, we have tried to explain a mechanism for the action of pentamidine, the best-known S100B ligand, and have proposed two S100B--pentamidine structures. The results presented here may be useful for the efficient design of new S100B ligands.