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.

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.

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.

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.

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.

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.

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.

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.

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.

Modeling protein structures with the coarse-grained UNRES force field in the CASP14 experiment.

The UNited RESidue (UNRES) force field was tested in the 14th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP14), in which larger oligomeric and multimeric targets were present compared to previous editions. Three prediction modes were tested (i) ab initio (the UNRES group), (ii) contact-assisted (the UNRES-contact group), and (iii) template-assisted (the UNRES-template group). For most of the targets, the contact restraints were derived from the server models top-ranked by the DeepQA method, while the DNCON2 method was used for 11 targets. Our consensus-fragment procedure was used to run template-assisted predictions. Each group also processed the Nuclear Magnetic Resonance (NMR)- and Small Angle X-Ray Scattering (SAXS)-data assisted targets. The average Global Distance Test Total Score (GDT_TS) of the 'Model 1' predictions were 29.17, 39.32, and 56.37 for the UNRES, UNRES-contact, and UNRES-template predictions, respectively, increasing by 0.53, 2.24, and 3.76, respectively, compared to CASP13. It was also found that the GDT_TS of the UNRES models obtained in ab initio mode and in the contact-assisted mode decreases with the square root of chain length, while the exponent in this relationship is 0.20 for the UNRES-template group models and 0.11 for the best performing AlphaFold2 models, which suggests that incorporation of database information, which stems from protein evolution, brings in long-range correlations, thus enabling the correction of force-field inaccuracies.

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.

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.

Scale-consistent approach to the derivation of coarse-grained force fields for simulating structure, dynamics, and thermodynamics of biopolymers.

In this chapter the scale-consistent approach to the derivation of coarse-grained force fields developed in our laboratory is presented, in which the effective energy function originates from the potential of mean force of the system under consideration and embeds atomistically detailed interactions in the resulting energy terms through use of Kubo's cluster-cumulant expansion, appropriate selection of the major degrees of freedom to be averaged out in the derivation of analytical approximations to the energy terms, and appropriate expression of the interaction energies at the all-atom level in these degrees of freedom. Our approach enables the developers to find correct functional forms of the effective coarse-grained energy terms, without having to import them from all-atom force fields or deriving them on a heuristic basis. In particular, the energy terms derived in such a way exhibit correct dependence on coarse-grained geometry, in particular on site orientation. Moreover, analytical formulas for the multibody (correlation) terms, which appear to be crucial for coarse-grained modeling of many of the regular structures such as, e.g., protein α-helices and ß-sheets, can be derived in a systematic way. Implementation of the developed theory to the UNIfied COarse-gRaiNed (UNICORN) model of biological macromolecules, which consists of the UNRES (for proteins), NARES-2P (for nucleic acids), and SUGRES-1P (for polysaccharides) components, and is being developed in our laboratory is described. Successful applications of UNICORN to the prediction of protein structure, simulating the folding and stability of proteins and nucleic acids, and solving biological problems are discussed.

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.

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.

Evaluation of the scale-consistent UNRES force field in template-free prediction of protein structures in the CASP13 experiment.

The recent NEWCT-9P version of the coarse-grained UNRES force field for proteins, with scale-consistent formulas for the local and correlation terms, has been tested in the CASP13 experiment of the blind-prediction of protein structure, in the ab initio, contact-assisted, and data-assisted modes. Significant improvement of the performance has been observed with respect to the CASP11 and CASP12 experiments (by over 10 GDT_TS units for the ab initio mode predictions and by over 15 GDT_TS units for the contact-assisted prediction, respectively), which is a result of introducing scale-consistent terms and improved handling of contact-distance restraints. As in previous CASP exercises, UNRES ranked higher in the free modeling category than in the general category that included template based modeling targets. Use of distance restraints from the predicted contacts, albeit many of them were wrong, resulted in the increase of GDT_TS by over 8 units on average and introducing sparse restraints from small-angle X-ray/neutron scattering and chemical cross-link-mass-spectrometry experiments, and ambiguous restraints from nuclear magnetic resonance experiments has also improved the predictions by 8.6, 9.7, and 10.7 GDT_TS units on average, respectively.

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.

Introduction of Phosphorylated Residues into the UNRES Coarse-Grained Model: Toward Modeling of Signaling Processes.

Phosphorylated proteins take part in many signaling pathways and play a key role in homeostasis regulation. The all-atom force fields enable us to study the systems containing phosphorylated proteins, but they are limited to short time scales. In this paper, we report the extension of the physics-based coarse-grained UNRES force field to treat systems with phosphorylated amino-acid residues. To derive the respective potentials, appropriate physics-based analytical expressions were fitted to the potentials of mean force of systems modeling phosphorylated amino-acid residues computed in our previous work and implemented in UNRES. The extended UNRES performed well in ab initio simulations of two miniproteins containing phosphorylated residues, strongly suggesting that realistic large-scale simulations of processes involving phosphorylated proteins, especially signaling processes, are now possible.

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, IcatCXPZ-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 IcatCXPZ-01 accumulated in the bone marrow reaching levels suitable for CatC inhibition. Subcutaneous administration of IcatCXPZ-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.

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.