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.

N-Terminally Lipidated Sialorphin Analogs-Synthesis, Molecular Modeling, In Vitro Effect on Enkephalins Degradation by NEP and Treatment of Intestinal Inflammation in Mice.

Pharmacotherapy for inflammatory bowel disease (IBD) is difficult, and some patients do not respond to currently available treatments. Therefore, the discovery of novel anti-IBD agents is imperative. Our aim was the synthesis of lipidated analogs of sialorphin and the in vitro characterization of their effect on the degradation of Met-enkephalin by neutral endopeptidase (NEP). We also investigated in vivo whether the most active inhibitor (peptide VIII) selected in the in vitro studies could be a potential candidate for the treatment of colitis. Peptides were synthesized by the solid-phase method. Molecular modeling technique was used to explain the effect of fatty acid chain length in sialorphin analogs on the ligand-enzyme interactions. The anti-inflammatory effect was evaluated in the dextran sulphate sodium (DSS)-induced model of colitis in mice. Peptide VIII containing stearic acid turned out to be in vitro the strongest inhibitor of NEP. We have also shown that the length of the chain of stearic acid fits the size of the grove of NEP. Peptides VII and VIII exhibited in vivo similar anti-inflammatory activity. Our results suggest that lipidation of sialorphin molecule is a promising direction in the search for NEP inhibitors that protect enkephalins.

Theoretical Investigation of the Coronavirus SARS-CoV-2 (COVID-19) Infection Mechanism and Selectivity.

The SARS-CoV-2 virus, commonly known as COVID-19, first occurred in December 2019 in Wuhan, Hubei Province, China. Since then, it has become a tremendous threat to human health. With a pandemic threat, it is in the significant interest of the scientific world to establish its method of infection. In this manuscript, we combine knowledge of the infection mechanism with theoretical methods to answer the question of the virus's selectivity. We proposed a two-stage infection mechanism. In the first step, the virus interacts with the ACE2 receptor, with the "proper strength". When the interaction is too strong, the virus will remain in an "improper position"; if the interaction is too weak, the virus will "run away" from the cell. We also indicated three residues (positions 30, 31, and 353) located on the ACE2 protein-binding interface, which seems to be crucial for successful infection. Our results indicate that these residues are necessary for the initiation of the infection process.

Low-Molecular Pyrazine-Based DNA Binders: Physicochemical and Antimicrobial Properties.

Pyrazine and its derivatives are a large group of compounds that exhibit broad biological activity, the changes of which can be easily detected by a substituent effect or a change in the functional group. The present studies combined theoretical research with the density functional theory (DFT) approach (B3LYP/6-311+G**) and experimental (potentiometric and spectrophotometric) analysis for a thorough understanding of the structure of chlorohydrazinopyrazine, its physicochemical and cytotoxic properties, and the site and nature of interaction with DNA. The obtained results indicated that 2-chloro-3-hydrazinopyrazine (2Cl3HP) displayed the highest affinity to DNA. Cytotoxicity studies revealed that the compound did not exhibit toxicity toward human dermal keratinocytes, which supported the potential application of 2Cl3HP in clinical use. The study also attempted to establish the possible equilibria occurring in the aqueous solution and, using both theoretical and experimental methods, clearly showed the hydrophilic nature of the compound. The experimental and theoretical results of the study confirmed the quality of the compound, as well as the appropriateness of the selected set of methods for similar research.

Prediction of protein assemblies, the next frontier: The CASP14-CAPRI experiment.

We present the results for CAPRI Round 50, the fourth joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of twelve targets, including six dimers, three trimers, and three higher-order oligomers. Four of these were easy targets, for which good structural templates were available either for the full assembly, or for the main interfaces (of the higher-order oligomers). Eight were difficult targets for which only distantly related templates were found for the individual subunits. Twenty-five CAPRI groups including eight automatic servers submitted ~1250 models per target. Twenty groups including six servers participated in the CAPRI scoring challenge submitted ~190 models per target. The accuracy of the predicted models was evaluated using the classical CAPRI criteria. The prediction performance was measured by a weighted scoring scheme that takes into account the number of models of acceptable quality or higher submitted by each group as part of their five top-ranking models. Compared to the previous CASP-CAPRI challenge, top performing groups submitted such models for a larger fraction (70-75%) of the targets in this Round, but fewer of these models were of high accuracy. Scorer groups achieved stronger performance with more groups submitting correct models for 70-80% of the targets or achieving high accuracy predictions. Servers performed less well in general, except for the MDOCKPP and LZERD servers, who performed on par with human groups. In addition to these results, major advances in methodology are discussed, providing an informative overview of where the prediction of protein assemblies currently stands.

Improved Consensus-Fragment Selection in Template-Assisted Prediction of Protein Structures with the UNRES Force Field in CASP13.

The method for protein-structure prediction, which combines the physics-based coarse-grained UNRES force field with knowledge-based modeling, has been developed further and tested in the 13th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP13). The method implements restraints from the consensus fragments common to server models. In this work, the server models to derive fragments have been chosen on the basis of quality assessment; a fully automatic fragment-selection procedure has been introduced, and Dynamic Fragment Assembly pseudopotentials have been fully implemented. The Global Distance Test Score (GDT_TS), averaged over our "Model 1" predictions, increased by over 10 units with respect to CASP12 for the free-modeling category to reach 40.82. Our "Model 1" predictions ranked 20 and 14 for all and free-modeling targets, respectively (upper 20.2% and 14.3% of all models submitted to CASP13 in these categories, respectively), compared to 27 (upper 21.1%) and 24 (upper 18.9%) in CASP12, respectively. For oligomeric targets, the Interface Patch Similarity (IPS) and Interface Contact Similarity (ICS) averaged over our best oligomer models increased from 0.28 to 0.36 and from 12.4 to 17.8, respectively, from CASP12 to CASP13, and top-ranking models of 2 targets (H0968 and T0997o) were obtained (none in CASP12). The improvement of our method in CASP13 over CASP12 was ascribed to the combined effect of the overall enhancement of server-model quality, our success in selecting server models and fragments to derive restraints, and improvements of the restraint and potential-energy functions.

A Peptidomimetic Fluorescent Probe to Detect the Trypsin ß2 Subunit of the Human 20S Proteasome.

This work describes the chemical synthesis, combinatorial selection, and enzymatic evaluation of peptidomimetic fluorescent substrates specific for the trypsin-like (ß2) subunit of the 20S human proteasome. After deconvolution of a library comprising nearly 6000 compounds composed of peg substituted diaminopropionic acid DAPEG building blocks, the sequence ABZ-Dap(O2(Cbz))-Dap(GO1)-Dap(O2(Cbz))-Arg-ANB-NH2, where ABZ is 2-aminobenzoic acid, and ANB- 5 amino 2- nitro benzoic acid was selected. Its cleavage followed sigmoidal kinetics, characteristic for allosteric enzymes, with Km = 3.22 ± 0.02 µM, kcat = 245 s-1, and kcat/Km = 7.61 × 107 M-1 s-1. This process was practically halted when a selective inhibitor of the ß2 subunit of the 20S human proteasome was supplemented to the reaction system. Titration of the substrate resulting in decreased amounts of proteasome 20S produced a linear signal up to 10-11 M. Using this substrate, we detected human proteasome 20S in human urine samples taken from the bladders of cancer patients. This observation could be useful for the noninvasive diagnosis of this severe disease.

1-Substituted sialorphin analogues-synthesis, molecular modelling and in vitro effect on enkephalins degradation by NEP.

Rat sialorphin (Gln-His-Asn-Pro-Arg) is a natural blocker of neprilysin (NEP) that belongs to the family of endogenous opioid peptide-degrading enzymes. Studies have confirmed the efficiency of sialorphin in blocking the activity of NEP, both in vitro and in vivo. It has been demonstrated that this inhibitor has a strong analgesic, anti-inflammatory, immunological and metabolic effect either directly or indirectly by affecting the level of Met/Leu-enkephalins. In this work, sialorphin and their 12 analogues were synthesised using the solid-phase method. The effect of the peptides on the degradation of Met-enkephalin by NEP and metabolic degradation in human plasma was investigated in vitro. We show that the change in the N-terminal amino acid configuration from L to D in almost all peptides, except D-Arg-His-Asn-Pro-Arg (peptide XI), led to the abolition of their inhibitory activity. With molecular modelling technique we explained the structural properties of the L and D-arginine located on the N-terminal part of the peptide. The detailed analysis of the protein binding pocket allowed us to explain why D-arginine is so unique among all D residues. Peptide XI showed the highest stability among the tested peptides in human plasma. For instance sialorphin after a 2-hour incubation in human plasma was almost completely decomposed, while the level of peptide XI dropped to 45% after 48 h under these conditions.

Processing and Maturation of Cathepsin C Zymogen: A Biochemical and Molecular Modeling Analysis.

Cysteine cathepsin C (CatC) is a ubiquitously expressed, lysosomal aminopeptidase involved in the activation of zymogens of immune-cell-associated serine proteinases (elastase, cathepsin G, proteinase 3, neutrophil serine proteinase 4, lymphocyte granzymes, and mast cell chymases). CatC is first synthetized as an inactive zymogen containing an intramolecular chain propeptide, the dimeric form of which is processed into the mature tetrameric form by proteolytic cleavages. A molecular modeling analysis of proCatC indicated that its propeptide displayed a similar fold to those of other lysosomal cysteine cathepsins, and could be involved in dimer formation. Our in vitro experiments revealed that human proCatC was processed and activated by CatF, CatK, and CatV in two consecutive steps of maturation, as reported for CatL and CatS previously. The unique positioning of the propeptide domains in the proCatC dimer complex allows this order of cleavages to be understood. The missense mutation Leu172Pro within the propeptide region associated with the Papillon-Lefèvre and Haim-Munk syndrome altered the proform stability as well as the maturation of the recombinant Leu172Pro proform.

Prediction of protein structure with the coarse-grained UNRES force field assisted by small X-ray scattering data and knowledge-based information.

A new approach to assisted protein-structure prediction has been proposed, which is based on running multiplexed replica exchange molecular dynamics simulations with the coarse-grained UNRES force field with restraints derived from knowledge-based models and distance distribution from small angle X-ray scattering (SAXS) measurements. The latter restraints are incorporated into the target function as a maximum-likelihood term that guides the shape of the simulated structures towards that defined by SAXS. The approach was first verified with the 1KOY protein, for which the distance distribution was calculated from the experimental structure, and subsequently used to predict the structures of 11 data-assisted targets in the CASP12 experiment. Major improvement of the GDT_TS was obtained for 2 targets, minor improvement for other 2 while, for 6 target GDT_TS deteriorated compared with that calculated for predictions without the SAXS data, partly because of assuming a wrong multimeric state (for Ts866) or because the crystal conformation was more compact than the solution conformation (for Ts942). Particularly good results were obtained for Ts909, in which use of SAXS data resulted in the selection of a correctly packed trimer and, subsequently, increased the GDT_TS of monomer prediction. It was found that running simulations with correct oligomeric state is essential for the success in SAXS-data-assisted prediction.

A new protein nucleic-acid coarse-grained force field based on the UNRES and NARES-2P force fields.

Based on the coarse-grained UNRES and NARES-2P models of proteins and nucleic acids, respectively, developed in our laboratory, in this work we have developed a coarse-grained model of systems containing proteins and nucleic acids. The UNRES and NARES-2P effective energy functions have been applied to the protein and nucleic-acid components of a system, respectively, while protein-nucleic-acid interactions have been described by the respective coarse-grained potentials developed in our recent work (Yin et al., J. Chem Theory Comput. 2015, 11, 1792). The Debye-Hückel screening has been applied to the electrostatic-interaction energy between the phosphate groups and charged amino-acid side chains. The model has been integrated into the UNRES package for coarse-grained molecular dynamics simulations of proteins and the implementation has been tested for energy conservation in microcanonical molecular dynamics runs and for temperature conservation in canonical molecular dynamics runs. Two case studies were performed: (i) the dynamics of the Ku protein heterodimer bound to DNA, for which it was found that the Ku70/Ku80 protein complex plays an active role in DNA repairing and (ii) conformational changes of the multiple antibiotic resistance (MarA) protein occurring during DNA binding, for which the functionally important motions occurring during this process were identified. © 2018 Wiley Periodicals, Inc.

Alanine scan of sialorphin and its hybrids with opiorphin: synthesis, molecular modelling and effect on enkephalins degradation.

Enkephalins are involved in a number of physiological processes. However, these peptides are quickly degraded by peptidases, e.g. the neutral endopeptidase (NEP). Inhibition of the enzymatic degradation of enkephalins is one of the possible approaches to prolong their activity. Selective inhibitor of NEP, sialorphin, is the attractive lead compound for enkephalins degradation studies. In this work, an alanine scan of sialorphin and a series of its hybrids with opiorphin, synthesised by the solid phase method, were performed. The effect of the peptides on degradation of Met-enkephalin by NEP in vitro was investigated. Molecular modelling technique was used to identify residues responsible for protein-ligand interactions. We showed that substitution of amino acids Gln1, Pro4 and Arg5 of sialorphin for Ala significantly reduced the half-life of Met-enkephalin in the presence of NEP. [Ala3]sialorphin displayed a higher inhibitory potency against NEP than sialorphin. Substitution of His2 for Ala led to a compound which was as active as lead compound. Sialorphin has a structure which hardly tolerates substitution in its sequence at positions 1, 4 and 5. The conversion of His2 for alanine in sialorphin is tolerated very well. The higher inhibitory potency of [Ala3]sialorphin than sialorphin against NEP is caused by removal of the hydrophilic residue (Asn) and a better fit of the peptide to the enzyme-binding pocket. The role of side chains of sialorphin in degradation of enkephalin by NEP has been explored. This study also provides an important SAR information essential for further drug design.

Intra- and intersubunit changes accompanying thermal activation of the HtrA2(Omi) protease homotrimer.

HtrA2(Omi) protease is involved in the maintenance of mitochondrial homeostasis and stimulation of apoptosis as well as in development of cancer and neurodegenerative disorders. The protein is a homotrimer whose subunits comprise serine protease domain (PD) and PDZ regulatory domain. In the basal, inactive state, a tight interdomain interface limits access both to the PDZ peptide (carboxylate) binding site and to the PD catalytic center. The molecular mechanism of activation is not well understood. To further the knowledge of HtrA2 thermal activation we monitored the dynamics of the PDZ-PD interactions during temperature increase using tryptophan-induced quenching (TrIQ) method. The TrIQ results suggested that during activation the PDZ domain changed its position versus PD inside a subunit, including a prominent change affecting the L3 regulatory loop of PD, and also changed its interactions with the PD of the adjacent subunit (PD*), specifically with its L1* regulatory loop containing the active site serine. The α5 helix of PDZ was involved in both, the intra- and intersubunit changes of interactions and thus seems to play an important role in HtrA2 activation. The amino acid substitutions designed to decrease the PDZ interactions with the PD or PD* promoted protease activity at a wide range of temperatures, which supports the conclusions based on the TrIQ analysis. The model presented in this work describes PDZ movement in relation to PD and PD*, resulting in an increased access to the peptide binding and active sites, and conformational changes of the L3 and L1* loops.

Performance of protein-structure predictions with the physics-based UNRES force field in CASP11.

Participating as the Cornell-Gdansk group, we have used our physics-based coarse-grained UNited RESidue (UNRES) force field to predict protein structure in the 11th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP11). Our methodology involved extensive multiplexed replica exchange simulations of the target proteins with a recently improved UNRES force field to provide better reproductions of the local structures of polypeptide chains. All simulations were started from fully extended polypeptide chains, and no external information was included in the simulation process except for weak restraints on secondary structure to enable us to finish each prediction within the allowed 3-week time window. Because of simplified UNRES representation of polypeptide chains, use of enhanced sampling methods, code optimization and parallelization and sufficient computational resources, we were able to treat, for the first time, all 55 human prediction targets with sizes from 44 to 595 amino acid residues, the average size being 251 residues. Complete structures of six single-domain proteins were predicted accurately, with the highest accuracy being attained for the T0769, for which the CαRMSD was 3.8 Å for 97 residues of the experimental structure. Correct structures were also predicted for 13 domains of multi-domain proteins with accuracy comparable to that of the best template-based modeling methods. With further improvements of the UNRES force field that are now underway, our physics-based coarse-grained approach to protein-structure prediction will eventually reach global prediction capacity and, consequently, reliability in simulating protein structure and dynamics that are important in biochemical processes. AVAILABILITY AND IMPLEMENTATION: Freely available on the web at http://www.unres.pl/ CONTACT: has5@cornell.edu.

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.

Molecular dynamics of protein A and a WW domain with a united-residue model including hydrodynamic interaction.

The folding of the N-terminal part of the B-domain of staphylococcal protein A (PDB ID: 1BDD, a 46-residue three-α-helix bundle) and the formin-binding protein 28 WW domain (PDB ID: 1E0L, a 37-residue three-stranded anti-parallel ß protein) was studied by means of Langevin dynamics with the coarse-grained UNRES force field to assess the influence of hydrodynamic interactions on protein-folding pathways and kinetics. The unfolded, intermediate, and native-like structures were identified by cluster analysis, and multi-exponential functions were fitted to the time dependence of the fractions of native and intermediate structures, respectively, to determine bulk kinetics. It was found that introducing hydrodynamic interactions slows down both the formation of an intermediate state and the transition from the collapsed structures to the final native-like structures by creating multiple kinetic traps. Therefore, introducing hydrodynamic interactions considerably slows the folding, as opposed to the results obtained from earlier studies with the use of Go-like models.

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.

Exploring the inhibitory potential of in silico-designed small peptides on Helicobacter pylori Hp0231 (DsbK), a periplasmic oxidoreductase involved in disulfide bond formation.

Introduction: Helicobacter pylori is a bacterium that colonizes the gastric epithelium, which affects millions of people worldwide. H. pylori infection can lead to various gastrointestinal diseases, including gastric adenocarcinoma and mucosa-associated lymphoid tissue lymphoma. Conventional antibiotic therapies face challenges due to increasing antibiotic resistance and patient non-compliance, necessitating the exploration of alternative treatment approaches. In this study, we focused on Hp0231 (DsbK), an essential component of the H. pylori Dsb (disulfide bond) oxidative pathway, and investigated peptide-based inhibition as a potential therapeutic strategy. Methods: Three inhibitory peptides designed by computational modeling were evaluated for their effectiveness using a time-resolved fluorescence assay. We also examined the binding affinity between Hp0231 and the peptides using microscale thermophoresis. Results and discussion: Our findings demonstrate that in silico-designed synthetic peptides can effectively inhibit Hp0231-mediated peptide oxidation. Targeting Hp0231 oxidase activity could attenuate H. pylori virulence without compromising bacterial viability. Therefore, peptide-based inhibitors of Hp0231 could be candidates for the development of new targeted strategy, which does not influence the composition of the natural human microbiome, but deprive the bacterium of its pathogenic properties.

Influence of Temperature and Salt Concentration on the Hydrophobic Interactions of Adamantane and Hexane.

One of the definitions of hydrophobic interactions is the aggregation of nonpolar particles in a polar solvent, such as water. While this phenomenon appears to be very simple, it is crucial for many complex processes, such as protein folding, to take place. In this work, the hydrophobic association of adamantane and hexane at various temperatures and ionic strengths was studied using molecular dynamics simulations with the AMBER 16.0 program and the GAFF force field. The potentials of mean force of hydrophobic dimer formation, as well as the excess free energy, excess energy, excess entropy, and excess heat capacity corresponding to the formation of the contact minimum, were determined and analyzed. For both systems, the depth of the contact minimum in the potential of mean force was found to increase with both temperature and ionic strength. The excess heat capacity of the association at the contact minimum and T = 298 K was found to be negative and to decrease, while the excess entropy and energy were found to be positive and to increase for both systems, the changes being more pronounced for the hexane dimer. The excess heat capacity is also greater in absolute value for the hexane dimer.

Peptidomimetics Based on C-Terminus of Blm10 Stimulate Human 20S Proteasome Activity and Promote Degradation of Proteins.

Degradation of misfolded, redundant and oxidatively damaged proteins constitutes one of the cellular processes which are influenced by the 20S proteasome. However, its activity is generally thought to decrease with age which leads to the gradual accumulation of abnormal proteins in cells and their subsequent aggregation. Therefore, increasing proteasomal degradation constitutes a promising strategy to delay the onset of various age-related diseases, including neurodegenerative disorders. In this study we designed and obtained a series of peptidomimetic stimulators of 20S comprising in their sequences the C-terminal fragment of Blm10 activator. Some of the compounds were capable of enhancing the degradation of natively unfolded and oxidatively damaged proteins, such as α-synuclein and enolase, whose applicability as proteasome substrates was evaluated by microscale thermophoresis (MST). Furthermore, they increased the ChT-L activity of the proteasome in HEK293T cell extracts. Our studies indicate that the 20S proteasome-mediated protein substrates hydrolysis may be selectively increased by peptide-based stimulators acting in an allosteric manner. These compounds, after further optimization, may have the potential to counteract proteasome impairment in patients suffering from age-related diseases.