Computational tools for exploring peptide-membrane interactions in gram-positive bacteria

The vital cellular functions in Gram-positive bacteria are controlled by signaling molecules known as quorum sensing peptides (QSPs), considered promising therapeutic interventions for bacterial infections. In the bacterial system QSPs bind to membrane-coupled receptors, which then auto-phosphorylate and activate intracellular response regulators. These response regulators induce target gene expression in bacteria. One of the most reliable trends in drug discovery research for virulence-associated molecular targets is the use of peptide drugs or new functionalities. In this perspective, computational methods act as auxiliary aids for biologists, where methodologies based on machine learning and in silico analysis are developed as suitable tools for target peptide identification. Therefore, the development of quick and reliable computational resources to identify or predict these QSPs along with their receptors and inhibitors is receiving considerable attention. The databases such as Quorumpeps and Quorum Sensing of Human Gut Microbes (QSHGM) provide a detailed overview of the structures and functions of QSPs. The tools and algorithms such as QSPpred, QSPred-FL, iQSP, EnsembleQS and PEPred-Suite have been used for the generic prediction of QSPs and feature representation. The availability of compiled key resources for utilizing peptide features based on amino acid composition, positional preferences, and motifs as well as structural and physicochemical properties, including biofilm inhibitory peptides, can aid in elucidating the QSP and membrane receptor interactions in infectious Gram-positive pathogens. Herein, we present a comprehensive survey of diverse computational approaches that are suitable for detecting QSPs and QS interference molecules. This review highlights the utility of these methods for developing potential biomarkers against infectious Gram-positive pathogens.

Repurposing of drugs for combined treatment of COVID-19 cytokine storm using machine learning

Severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) induced cytokine storm is the major cause of COVID-19 related deaths. Patients have been treated with drugs that work by inhibiting a specific protein partly responsible for the cytokines production. This approach provided very limited success, since there are multiple proteins involved in the complex cell signaling disease mechanisms. We targeted five proteins: Angiotensin II receptor type 1 (AT1R), A disintegrin and metalloprotease 17 (ADAM17), Nuclear Factor‑Kappa B (NF‑κB), Janus kinase 1 (JAK1) and Signal Transducer and Activator of Transcription 3 (STAT3), which are involved in the SARS‑CoV‑2 induced cytokine storm pathway. We developed machine-learning (ML) models for these five proteins, using known active inhibitors. After developing the model for each of these proteins, FDA-approved drugs were screened to find novel therapeutics for COVID‑19. We identified twenty drugs that are active for four proteins with predicted scores greater than 0.8 and eight drugs active for all five proteins with predicted scores over 0.85. Mitomycin C is the most active drug across all five proteins with an average prediction score of 0.886. For further validation of these results, we used the PyRx software to conduct protein-ligand docking experiments and calculated the binding affinity. The docking results support findings by the ML model. This research study predicted that several drugs can target multiple proteins simultaneously in cytokine storm-related pathway. These may be useful drugs to treat patients because these therapies can fight cytokine storm caused by the virus at multiple points of inhibition, leading to synergistically effective treatments.

Bioinformatics evaluation of anticancer properties of GP63 protein-derived peptides on MMP2 protein of melanoma cancer

Background

GP63, also known as Leishmanolysin, is a multifunctional virulence factor abundant on the surface of Leishmania spp. small peptides with anticancer capabilities that are selective and toxic to cancer cells are known as anticancer peptides. We aimed to demonstrate the activity of GP63 and its anticancer properties on melanoma using a range of in silico tools and screening methods to identify predicted and designed anticancer peptides.

Methods

Various in silico modeling methodologies are used to establish the three-dimensional (3D) structure of GP63. Refinement and re-evaluation of the modeled structures and the built models' quality evaluated using the different docking used to find the interacting amino acids between MMP2 and GP63 and its anticancer peptides. AntiCP2.0 is used for screening anticancer peptides. 2D interaction plots of protein-ligand complexes evaluated by Protein-Ligand Interaction Profiler server. It is for the first time that used anticancer peptides of GP63 and the predicted and designed peptides.

Results

We used 3 peptides of GP63 based on the AntiCP 2.0 server with scores of 0.63, 0.53, and 0.49, and common peptides of GP63/MMP2 (continues peptide: mean the completely selected peptide after docking with non-anticancer effect, predicted with 0.58 score and designed peptides with 0.47 and 0.45 scores by AntiCP 2.0 server).

Conclusions

The antileishmanial and anticancer peptide research topics exemplify the multidisciplinary nature of peptide research. The advancement of therapeutics targeting cancer and/or Leishmania requires an interconnected research strategy shown in this work.

Recent developments on UDP-N-acetylmuramoyl-L-alanine-D-gutamate ligase (Mur D) enzyme for antimicrobial drug development: An emphasis on in-silico approaches

Introduction

The rapid emergence of antibiotic resistance among various bacterial pathogens has been one of the major concerns of health organizations across the world. In this context, for the development of novel inhibitors against antibiotic-resistant bacterial pathogens, UDP-N-Acetylmuramoyl-L-Alanine-D-Glutamate Ligase (MurD) enzyme represents one of the most apposite targets.

Body

The present review focuses on updated advancements on MurD-targeted inhibitors in recent years along with genetic regulation, structural and functional characteristics of the MurD enzyme from various bacterial pathogens. A concise account of various crystal structures of MurD enzyme, submitted into Protein Data Bank is also discussed.

Discussion

MurD, an ATP dependent cytoplasmic enzyme is an important target for drug discovery. The genetic organization of MurD enzyme is well elucidated and many crystal structures of MurD enzyme are submitted into Protein Data bank. Various inhibitors against MurD enzyme have been developed so far with an increase in the use of in-silico methods in the recent past. But cell permeability barriers and conformational changes of MurD enzyme during catalytic reaction need to be addressed for effective drug development. So, a combination of in-silico methods along with experimental work is proposed to counter the catalytic machinery of MurD enzyme.

Molecular docking studies of phytochemicals from Terminalia chebula for identification of potential multi-target inhibitors of SARS-CoV-2 proteins

Background

The COVID-19 caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has emerged as a global pandemic claiming more than 6 million lives worldwide as of 16 March 2022. Till date, no medicine has been developed which is proved to have 100% efficiency in combating against this deadly disease. We focussed on ayurvedic medicines to identify drug-like candidates for treatment and management of COVID-19. Among all ayurvedic medicines, we were interested in Terminalia chebula (T. chebula), as it is known to have antibacterial, antifungal, antiviral, antioxidant and anti-inflammatory properties.

Objectives

In this study, we evaluated potential inhibitory effects of phytochemicals from T. chebula against eight structural and functional proteins of SARS-CoV-2.

Material and methods

We performed blind molecular docking studies using fifteen phytochemicals from T. chebula against the proteins of SARS-CoV-2. The three-dimensional proteins structures were analysed and potential drug-binding sites were identified. The drug-likeness properties of the ligands were assessed as well.

Results

Analysing the docking results by comparing Atomic Contact Energy (ACE) and intermolecular interactions along with assessment of ADME/T properties identified 1,3,6-Trigalloyl glucose (−332.14 ± 55.74 kcal/mol), Beta-Sitosterol (−324.75 ± 36.98 kcal/mol) and Daucosterol (−335.67 ± 104.79 kcal/mol) as most promising candidates which exhibit significantly high inhibition efficiency against all eight protein targets.

Conclusions

We believe that our study has the potential to help the scientific communities to develop multi-target drugs from T. chebula to combat against the deadly pathogen of COVID-19, with the support of extensive wet lab analysis.

Neomenthol prevents the proliferation of skin cancer cells by restraining tubulin polymerization and hyaluronidase activity

Introduction

Neomenthol, a cyclic monoterpenoid, is a stereoisomer of menthol present in the essential oil of Mentha spp. It is used in food as a flavoring agent, in cosmetics and medicines because of its cooling effects. However, neomenthol has not been much explored for its anticancer potential. Additionally, targeting hyaluronidase, Cathepsin-D, and ODC by phytochemicals is amongst the efficient approach for cancer prevention and/or treatment.

Objectives

To investigate the molecular and cell target-based antiproliferative potential of neomenthol on human cancer (A431, PC-3, K562, A549, FaDu, MDA-MB-231, COLO-205, MCF-7, and WRL-68) and normal (HEK-293) cell lines.

Methods

The potency of neomenthol was evaluated on human cancer and normal cell line using SRB, NRU and MTT assays. The molecular target based study of neomenthol was carried out in cell-free and cell-based test systems. Further, the potency of neomenthol was confirmed by quantitative real-time PCR analysis and molecular docking studies. The in vivo anticancer potential of neomenthol was performed on mice EAC model and the toxicity examination was accomplished through in silico, ex vivo and in vivo approaches.

Results

Neomenthol exhibits a promising activity (IC₅₀ 17.3 ± 6.49 μM) against human epidermoid carcinoma (A431) cells by arresting the G2/M phase and increasing the number of sub-diploid cells. It significantly inhibits hyaluronidase activity (IC₅₀ 12.81 ± 0.01 μM) and affects the tubulin polymerization. The expression analysis and molecular docking studies support the in vitro molecular and cell target based results. Neomenthol prevents EAC tumor formation by 58.84% and inhibits hyaluronidase activity up to 10% at 75 mg/kg bw, i.p. dose. The oral dose of 1000 mg/kg bw was found safe in acute oral toxicity studies.

Conclusion

Neomenthol delayed the growth of skin carcinoma cells by inhibiting the tubulin polymerization and hyaluronidase activity, which are responsible for tumor growth, metastasis, and angiogenesis.

Dataset for dynamics and conformational changes in human NEIL2 protein analyzed by integrative structural biology approach

This work presents new data on human endonuclease VIII-like 2 protein (hNEIL2), a part of DNA glycosylases of the helix–two-turn–helix structural superfamily. While X-ray structure of oNEIL2 (opossum Monodelphis) was resolved partially [1], the structure of hNEIL2 has not yet been determined. This data article describes a powerful combination of hydrogen-deuterium exchange mass spectrometry, homology modeling, and molecular dynamics simulations for protein conformational dynamics analysis.

The data supplied in this work are related to the research article entitled “Dynamics and Conformational Changes in Human NEIL2 DNA Glycosylase Analyzed by Hydrogen/Deuterium Exchange Mass Spectrometry”.

Analysis of SARS-COV2 spike protein variants among Iraqi isolates

The ongoing pandemic of COVID-19 caused by the SARS-COV2 virus has triggered millions of deaths around the globe. Emerging several variants of the virus with increased transmissibility, the severity of disease, and the ability of the virus to escape from the immune system has a cause for concerns. Here, we compared the spike protein sequence of 91 human SARS CoV2 strains of Iraq to the first reported sequence of SARS-CoV2 isolate from Wuhan Hu-1/China. The strains were isolated between June 2020 and March 2021. Twenty-two distinct mutations were identified within the spike protein regions which were: L5F, L18F, T19R, S151T, G181A, A222V, A348S, L452 (Q or M), T478K, N501Y, A520S, A522V, A570D, S605A, D614G, Q675H, N679K, P681H, T716I, S982A, A1020S, D1118H. The most frequently mutations occurred at the D614G (87/91), followed by S982A (50/91), and A570D (48/91), respectively. In addition, a distinct shift was observed in the type of SARS-COV2 variants present in 2020 compared to 2021 isolates. In 2020, B.1.428.1 lineage was appeared to be a dominant variant (85%). However, the diversity of the variants increased in 2021, and the majority (73%) of the isolated were appeared to belong to B.1.1.7 lineage (VOC/alpha variants). To our knowledge, this is the first major genome analysis of SARS-CoV2 in Iraq. The data from this research could provide insights into SARS-CoV2 evolution, and can be potentially used to recognize the effective vaccine against the disease.

Recent developments in the medicinal chemistry of single boron atom-containing compounds

Various boron-containing drugs have been approved for clinical use over the past two decades, and more are currently in clinical trials. The increasing interest in boron-containing compounds is due to their unique binding properties to biological targets; for example, boron substitution can be used to modulate biological activity, pharmacokinetic properties, and drug resistance. In this perspective, we aim to comprehensively review the current status of boron compounds in drug discovery, focusing especially on progress from 2015 to December 2020. We classify these compounds into groups showing anticancer, antibacterial, antiviral, antiparasitic and other activities, and discuss the biological targets associated with each activity, as well as potential future developments.

MCDB: A comprehensive curated mitotic catastrophe database for retrieval, protein sequence alignment, and target prediction

Mitotic catastrophe (MC) is a form of programmed cell death induced by mitotic process disorders, which is very important in tumor prevention, development, and drug resistance. Because rapidly increased data for MC is vigorously promoting the tumor-related biomedical and clinical study, it is urgent for us to develop a professional and comprehensive database to curate MC-related data. Mitotic Catastrophe Database (MCDB) consists of 1214 genes/proteins and 5014 compounds collected and organized from more than 8000 research articles. Also, MCDB defines the confidence level, classification criteria, and uniform naming rules for MC-related data, which greatly improves data reliability and retrieval convenience. Moreover, MCDB develops protein sequence alignment and target prediction functions. The former can be used to predict new potential MC-related genes and proteins, and the latter can facilitate the identification of potential target proteins of unknown MC-related compounds. In short, MCDB is such a proprietary, standard, and comprehensive database for MC-relate data that will facilitate the exploration of MC from chemists to biologists in the fields of medicinal chemistry, molecular biology, bioinformatics, oncology and so on. The MCDB is distributed on http://www.combio-lezhang.online/MCDB/index_html/.

Co-crystallization and structure determination: An effective direction for anti-SARS-CoV-2 drug discovery

Safer and more-effective drugs are urgently needed to counter infections with the highly pathogenic SARS-CoV-2, cause of the COVID-19 pandemic. Identification of efficient inhibitors to treat and prevent SARS-CoV-2 infection is a predominant focus. Encouragingly, using X-ray crystal structures of therapeutically relevant drug targets (PLpro, Mpro, RdRp, and S glycoprotein) offers a valuable direction for anti-SARS-CoV-2 drug discovery and lead optimization through direct visualization of interactions. Computational analyses based primarily on MMPBSA calculations have also been proposed for assessing the binding stability of biomolecular structures involving the ligand and receptor. In this study, we focused on state-of-the-art X-ray co-crystal structures of the abovementioned targets complexed with newly identified small-molecule inhibitors (natural products, FDA-approved drugs, candidate drugs, and their analogues) with the assistance of computational analyses to support the precision design and screening of anti-SARS-CoV-2 drugs.

Toxicity evaluation of Eleutherine plicata Herb. extracts and possible cell death mechanism

•

Ethanol extract of Eleutherine plicata showed low in vitro and in vivo cytotoxic potential.

•

The dichloromethane fraction was cytotoxic to HepG2 and caused DNA. However, no toxicity was observed in vivo.

•

Isoeleutherin caused DNA damage by the comet method and activated caspase-8 in the in silico study.

Eleutherine plicata has been shown to be a promising medicinal plant, and its activity has been associated with naphthoquinones. The present study aimed at evaluating the cytotoxicity, genotoxicity, and oral toxicity of the ethanol extract (EEEp), dichloromethane fraction (FDMEp) of E. plicata, and isoeleutherin. For the cytotoxicity evaluation, the viability test (MTT) was used. Genotoxicity was accessed through the Comet assay (alkaline version), acute and subacute oral toxicities were also evaluated. The antioxidant capacity of the samples in the wells where the cells were treated with E. plicata was evaluated. Furthermore, the participation of caspase-8 in the possible mechanism of action of isoeleutherin, eleutherin, and eleutherol was also investigated through a docking study. FDMEp and isoeleutherin were cytotoxic, with higher rates of DNA fragmentation observed for FDMEp and isoeleutherin, and all samples displayed higher antioxidant potential than the control. In the acute oral toxicity test, EEEp, FDMEp, and isoeleutherin did not cause significant clinical changes. In the subacute toxicity assay, EEEp and FDMEp also did not cause clinical, hematological, or biochemical changes. The three compounds bound similarly to caspase-8. Despite the results of cytotoxicity, in vitro studies demonstrated that the use of EEEp appears to be safe and cell death may involve its binding to caspase-8.

Structural characterization of two solute-binding proteins for N,N'-diacetylchitobiose/N,N',N''-triacetylchitotoriose of the gram-positive bacterium, Paenibacillus sp. str. FPU-7

The chitinolytic bacterium Paenibacillus sp. str. FPU-7 efficiently degrades chitin into oligosaccharides such as N-acetyl-D-glucosamine (GlcNAc) and disaccharides (GlcNAc)₂ through multiple secretory chitinases. Transport of these oligosaccharides by P. str. FPU-7 has not yet been clarified. In this study, we identified nagB1, predicted to encode a sugar solute-binding protein (SBP), which is a component of the ABC transport system. However, the genes next to nagB1 were predicted to encode two-component regulatory system proteins rather than transmembrane domains (TMDs). We also identified nagB2, which is highly homologous to nagB1. Adjacent to nagB2, two genes were predicted to encode TMDs. Binding experiments of the recombinant NagB1 and NagB2 to several oligosaccharides using differential scanning fluorimetry and surface plasmon resonance confirmed that both proteins are SBPs of (GlcNAc)₂ and (GlcNAc)₃. We determined their crystal structures complexed with and without chitin oligosaccharides at a resolution of 1.2 to 2.0 Å. The structures shared typical SBP structural folds and were classified as subcluster D-I. Large domain motions were observed in the structures, suggesting that they were induced by ligand binding via the "Venus flytrap" mechanism. These structures also revealed chitin oligosaccharide recognition mechanisms. In conclusion, our study provides insight into the recognition and transport of chitin oligosaccharides in bacteria.

Zeta-carbonic anhydrases show CS2 hydrolase activity: A new metabolic carbon acquisition pathway in diatoms?

•

CDCA1 is a carbonic anhydrase that can utilize Zn(II) or Cd(II) as catalytic metal.

•

CDCA1 has been the first enzyme to show an efficient utilization of Cd(II) ions in Nature.

•

By using a multidisciplinary approach, we discovered that CS₂ is a substrate for this enzyme.

•

CDCA1 is the unique enzyme, known so far, able to use both CS₂ and CO₂ as substrates.

CDCA1 is a very peculiar member of the Carbonic Anhydrase (CA) family. It has been the first enzyme to show an efficient utilization of Cd(II) ions in Nature and a unique adaptation capability to live on the surface ocean. Indeed, in this environment, which is extremely depleted in essential metal ions, CDCA1 can utilize Zn(II) or Cd(II) as catalytic metal to support the metabolic needs of fast growing diatoms. In this paper we demonstrate a further catalytic versatility of this enzyme by using a combination of X-ray crystallography, molecular dynamics simulations and enzymatic experiments. First we identified the CO₂ binding site and the way in which this substrate travels from the environment to the enzyme active site. Then, starting from the observation of a structural similarity with the substrate entry route of CS₂ hydrolase from Acidanius A1-3, we hypothesized and demonstrated that also CS₂ is a substrate for CDCA1. This finding is new and unexpected since until now only few CS₂ hydrolases have been characterized, and none of them is reported to have any CO₂ hydratase action. The physiological implications of this supplementary catalytic activity still remain to be unveiled. We suggest here that it could represent another ability of diatoms expressing CDCA1 to adapt to the external environment. Indeed, the ability of this enzyme to convert CS₂ could represent an alternative source of carbon acquisition for diatoms, in addition to CO₂.

In silico evaluation of flavonoids as effective antiviral agents on the spike glycoprotein of SARS-CoV-2

The novel coronavirus pandemic has spread over in 213 countries as of July 2020. Approximately 12 million people have been infected so far according to the reports from World Health Organization (WHO). Preventive measures are being taken globally to avoid the rapid spread of virus. In the current study, an in silico approach is carried out as a means of inhibiting the spike protein of the novel coronavirus by flavonoids from natural sources that possess both antiviral and anti-inflammatory properties. The methodology is focused on molecular docking of 10 flavonoid compounds that are docked with the spike protein of SARS-CoV-2, to determine the highest binding affinity at the binding site. Molecular dynamics simulation was carried out with the flavonoid-protein complex showing the highest binding affinity and highest interactions. The flavonoid naringin showed the least binding energy of -9.8 Kcal/mol with the spike protein which was compared with the standard drug, dexamethasone which is being repurposed to treat critically ill patients. MD simulation was carried out on naringin-spike protein complex for their conformational stability in the active site of the novel coronavirus spike protein. The RMSD of the complex appeared to be more stable when compared to that of the protein from 0.2 nm to 0.4 nm. With the aid of this in silico approach further in vitro studies can be carried out on these flavonoids against the novel coronavirus as a means of viral protein inhibitors.

Nicotinic cholinergic system and COVID-19: In silico evaluation of nicotinic acetylcholine receptor agonists as potential therapeutic interventions

SARS-CoV-2 infection was announced as a pandemic in March 2020. Since then, several scientists have focused on the low prevalence of smokers among hospitalized COVID-19 patients. These findings led to our hypothesis that the Nicotinic Cholinergic System (NCS) plays a crucial role in the manifestation of COVID-19 and its severe symptoms. Molecular modeling revealed that the SARS-CoV-2 Spike glycoprotein might bind to nicotinic acetylcholine receptors (nAChRs) through a cryptic epitope homologous to snake toxins, substrates well documented and known for their affinity to the nAChRs. This binding model could provide logical explanations for the acute inflammatory disorder in patients with COVID-19, which may be linked to severe dysregulation of NCS. In this study, we present a series of complexes with cholinergic agonists that can potentially prevent SARS-CoV-2 Spike glycoprotein from binding to nAChRs, avoiding dysregulation of the NCS and moderating the symptoms and clinical manifestations of COVID-19. If our hypothesis is verified by in vitro and in vivo studies, repurposing agents currently approved for smoking cessation and neurological conditions could provide the scientific community with a therapeutic option in severe COVID-19.

Crystal structure of the atypically adhesive SpaB basal pilus subunit: Mechanistic insights about its incorporation in lactobacillar SpaCBA pili

To successfully colonize a host or environment, certain genera and species of Gram-positive bacteria have evolved to utilize the so-called sortase-dependent pilus, a long multi-subunit and non-flagellar surface adhesin. One example of this is Lactobacillus rhamnosus GG, a gut-adapted probiotic strain that produces SpaCBA pili. These structures are covalent hetero-oligomers built from three types of pilin subunit, each with a specific location and function (i.e., backbone SpaA for length, tip SpaC for adhesion, and basal SpaB for anchoring). Functionally, the SpaCBA pilus exhibits a promiscuous affinity for components on intestinal surfaces (e.g., mucus, collagen, and epithelial cells), which is largely attributed to the SpaC subunit. Then again, the basal SpaB pilin, in addition to acting as the terminal subunit during pilus assembly, displays an out of character mucoadhesive function. To address the structural basis of this unusual dual functionality, we reveal the 2.39 ​Å resolution crystal structure of SpaB. SpaB consists of one immunoglobulin-like CnaB domain and contains a putative intermolecular isopeptide bond-linking lysine and internal isopeptide bond-asparagine in an FPKN pilin motif within the C-terminal end. Remarkably, we found that a C-terminal stretch of positively charged lysine and arginine residues likely accounts for the atypical mucoadhesiveness of SpaB. Although harboring an autocatalytic triad of residues for a potential internal isopeptide interaction, the SpaB crystal structure lacked the visible electron density for intact bond formation, yet its presence was subsequently confirmed by mass spectral analysis. Finally, we propose a structural model that captures the exclusive basal positioning of SpaB in the SpaCBA pilus.

Advances in integrative structural biology: Towards understanding protein complexes in their cellular context

Microorganisms rely on protein interactions to transmit signals, react to stimuli, and grow. One of the best ways to understand these protein interactions is through structural characterization. However, in the past, structural knowledge was limited to stable, high-affinity complexes that could be crystallized. Recent developments in structural biology have revolutionized how protein interactions are characterized. The combination of multiple techniques, known as integrative structural biology, has provided insight into how large protein complexes interact in their native environment. In this mini-review, we describe the past, present, and potential future of integrative structural biology as a tool for characterizing protein interactions in their cellular context.

Unlocking COVID therapeutic targets: A structure-based rationale against SARS-CoV-2, SARS-CoV and MERS-CoV Spike

There are no approved target therapeutics against SARS-CoV-2 or other beta-CoVs. The beta-CoV Spike protein is a promising target considering the critical role in viral infection and pathogenesis and its surface exposed features. We performed a structure-based strategy targeting highly conserved druggable regions resulting from a comprehensive large-scale sequence analysis and structural characterization of Spike domains across SARSr- and MERSr-CoVs. We have disclosed 28 main consensus druggable pockets within the Spike. The RBD and SD1 (S1 subunit); and the CR, HR1 and CH (S2 subunit) represent the most promising conserved druggable regions. Additionally, we have identified 181 new potential hot spot residues for the hSARSr-CoVs and 72 new hot spot residues for the SARSr- and MERSr-CoVs, which have not been described before in the literature. These sites/residues exhibit advantageous structural features for targeted molecular and pharmacological modulation. This study establishes the Spike as a promising anti-CoV target using an approach with a potential higher resilience to resistance development and directed to a broad spectrum of Beta-CoVs, including the new SARS-CoV-2 responsible for COVID-19. This research also provides a structure-based rationale for the design and discovery of chemical inhibitors, antibodies or other therapeutic modalities successfully targeting the Beta-CoV Spike protein.

Spectroscopic and molecular docking studies for characterizing binding mechanism and conformational changes of human serum albumin upon interaction with Telmisartan

Human serum albumin (HSA), one of the most copious plasma proteins is responsible for binding and transportation of many exogenous and endogenous ligands including drugs. In this study, we intended to explore the extent and types of binding interaction present between HSA and the antihypertensive drug, telmisartan (TLM). The conformational changes in HSA due to this binding were also studied using different spectroscopic and molecular docking techniques. The spectral shifting and intensity variations upon interaction with TLM were studied using FT-IR spectroscopy. Binding constant and the change in absorption of HSA at its λ_max was analyzed using absorption spectroscopy. Eventually, the types and extent of binding interactions were confirmed using molecular docking technique. Results have shown that TLM significantly interacts with the binding site-1 of HSA utilizing strong hydrogen bonding with Glu292, and Lys195 residues. The UV-absorption intensities were found to be decreased serially as the drug concentration increased with a binding constant of 1.01 × 10³ M⁻¹. The secondary structure analysis using FT-IR spectroscopy also revealed a marked reduction in the α-helix (56%) component of HSA on interaction. This study gives critical insights into the interaction of TLM with HSA protein which eventually affects the concentration of TLM reaching the site of action and ultimately its therapeutic profile.

Structural insights into human Arginase-1 pH dependence and its inhibition by the small molecule inhibitor CB-1158

Arginase-1 is a manganese-dependent metalloenzyme that catalyzes the hydrolysis of L-arginine into L-ornithine and urea. Arginase-1 is abundantly expressed by tumor-infiltrating myeloid cells that promote tumor immunosuppression, which is relieved by inhibition of Arginase-1. We have characterized the potencies of the Arginase-1 reference inhibitors (2S)-2-amino-6-boronohexanoic acid (ABH) and N ^ω-hydroxy-nor-L-arginine (nor-NOHA), and studied their pH-dependence and binding kinetics. To gain a better understanding of the structural changes underlying the high pH optimum of Arginase-1 and its pH-dependent inhibition, we determined the crystal structure of the human Arginase-1/ABH complex at pH 7.0 and 9.0. These structures revealed that at increased pH, the manganese cluster assumes a more symmetrical coordination structure, which presumably contributes to its increase in catalytic activity. Furthermore, we show that binding of ABH involves the presence of a sodium ion close to the manganese cluster. We also studied the investigational new drug CB-1158 (INCB001158). This inhibitor has a low-nanomolar potency at pH 7.4 and increases the thermal stability of Arginase-1 more than ABH and nor-NOHA. Moreover, CB-1158 displays slow association and dissociation kinetics at both pH 9.5 and 7.4, as indicated by surface plasmon resonance. The potent character of CB-1158 is presumably due to its increased rigidity compared to ABH as well as the formation of an additional hydrogen-bond network as observed by resolution of the Arginase-1/CB-1158 crystal structure.

MCC1019, a selective inhibitor of the Polo-box domain of Polo-like kinase 1 as novel, potent anticancer candidate

Polo-like kinase (PLK1) has been identified as a potential target for cancer treatment. Although a number of small molecules have been investigated as PLK1 inhibitors, many of which showed limited selectivity. PLK1 harbors a regulatory domain, the Polo box domain (PBD), which has a key regulatory function for kinase activity and substrate recognition. We report on 3-bromomethyl-benzofuran-2-carboxylic acid ethyl ester (designated: MCC1019) as selective PLK1 inhibitor targeting PLK1 PBD. Cytotoxicity and fluorescence polarization-based screening were applied to a library of 1162 drug-like compounds to identify potential inhibitors of PLK1 PBD. The activity of compound MC1019 against the PLK1 PBD was confirmed using fluorescence polarization and microscale thermophoresis. This compound exerted specificity towards PLK1 over PLK2 and PLK3. MCC1019 showed cytotoxic activity in a panel of different cancer cell lines. Mechanistic investigations in A549 lung adenocarcinoma cells revealed that MCC1019 induced cell growth inhibition through inactivation of AKT signaling pathway, it also induced prolonged mitotic arrest—a phenomenon known as mitotic catastrophe, which is followed by immediate cell death via apoptosis and necroptosis. MCC1019 significantly inhibited tumor growth in vivo in a murine lung cancer model without affecting body weight or vital organ size, and reduced the growth of metastatic lesions in the lung. We propose MCC1019 as promising anti-cancer drug candidate.

Human Albumin Impairs Amyloid β-peptide Fibrillation Through its C-terminus: From docking Modeling to Protection Against Neurotoxicity in Alzheimer's disease

Alzheimer's disease (AD) is a neurodegenerative process characterized by the accumulation of extracellular deposits of amyloid β-peptide (Aβ), which induces neuronal death. Monomeric Aβ is not toxic but tends to aggregate into β-sheets that are neurotoxic. Therefore to prevent or delay AD onset and progression one of the main therapeutic approaches would be to impair Aβ assembly into oligomers and fibrils and to promote disaggregation of the preformed aggregate. Albumin is the most abundant protein in the cerebrospinal fluid and it was reported to bind Aβ impeding its aggregation. In a previous work we identified a 35-residue sequence of clusterin, a well-known protein that binds Aβ, that is highly similar to the C-terminus (CTerm) of albumin. In this work, the docking experiments show that the average binding free energy of the CTerm-Aβ1-42 simulations was significantly lower than that of the clusterin-Aβ1-42 binding, highlighting the possibility that the CTerm retains albumin's binding properties. To validate this observation, we performed in vitro structural analysis of soluble and aggregated 1 μM Aβ1-42 incubated with 5 μM CTerm, equimolar to the albumin concentration in the CSF. Reversed-phase chromatography and electron microscopy analysis demonstrated a reduction of Aβ1-42 aggregates when the CTerm was present. Furthermore, we treated a human neuroblastoma cell line with soluble and aggregated Aβ1-42 incubated with CTerm obtaining a significant protection against Aβ-induced neurotoxicity. These in silico and in vitro data suggest that the albumin CTerm is able to impair Aβ aggregation and to promote disassemble of Aβ aggregates protecting neurons.

Structural insights into chemokine CCL17 recognition by antibody M116

The homeostatic chemokine CCL17, also known as thymus and activation regulated chemokine (TARC), has been associated with various diseases such as asthma, idiopathic pulmonary fibrosis, atopic dermatitis and ulcerative colitis. Neutralization of CCL17 by antibody treatment ameliorates the impact of disease by blocking influx of T cells. Monoclonal antibody M116 derived from a combinatorial library shows potency in neutralizing CCL17-induced signaling. To gain insight into the structural determinants of antigen recognition, the crystal structure of M116 Fab was determined in complex with CCL17 and in the unbound form. Comparison of the structures revealed an unusual induced-fit mechanism of antigen recognition that involves cis-trans isomerization in two CDRs. The structure of the CCL17-M116 complex revealed the antibody binding epitope, which does not overlap with the putative receptor epitope, suggesting that the current model of chemokine-receptor interactions, as observed in the CXCR4-vMIP-II system, may not be universal.

Structural distortions due to missense mutations in human formylglycine-generating enzyme leading to multiple sulfatase deficiency

The major candidate for multiple sulfatase deficiency is a defective formylglycine-generating enzyme (FGE). Though adequately produced, mutations in FGE stall the activation of sulfatases and prevent their activity. Missense mutations, viz. E130D, S155P, A177P, W179S, C218Y, R224W, N259I, P266L, A279V, C336R, R345C, A348P, R349Q and R349W associated with multiple sulfatase deficiency are yet to be computationally studied. Aforementioned mutants were initially screened through ws-SNPs&GO^3D program. Mutant R345C acquired the highest score, and hence was studied in detail. Discrete molecular dynamics explored structural distortions due to amino acid substitution. Therein, comparative analyses of wild type and mutant were carried out. Changes in structural contours were observed between wild type and mutant. Mutant had low conformational fluctuation, high atomic mobility and more compactness than wild type. Moreover, free energy landscape showed mutant to vary in terms of its conformational space as compared to wild type. Subsequently, wild type and mutant were subjected to single-model analyses. Mutant had lesser intra molecular interactions than wild type suggesting variations pertaining to its secondary structure. Furthermore, simulated thermal denaturation showed dissimilar pattern of hydrogen bond dilution. Effects of these variations were observed as changes in elements of secondary structure. Docking studies of mutant revealed less favourable binding energy towards its substrate as compared to wild type. Therefore, theoretical explanations for structural distortions of mutant R345C leading to multiple sulfatase deficiency were revealed. The protocol of the study could be useful to examine the effectiveness of pharmacological chaperones prior to experimental studies.

Camelid Ig V genes reveal significant human homology not seen in therapeutic target genes, providing for a powerful therapeutic antibody platform

Camelid immunoglobulin variable (IGV) regions were found homologous to their human counterparts; however, the germline V repertoires of camelid heavy and light chains are still incomplete and their therapeutic potential is only beginning to be appreciated. We therefore leveraged the publicly available HTG and WGS databases of Lama pacos and Camelus ferus to retrieve the germline repertoire of V genes using human IGV genes as reference. In addition, we amplified IGKV and IGLV genes to uncover the V germline repertoire of Lama glama and sequenced BAC clones covering part of the Lama pacos IGK and IGL loci. Our in silico analysis showed that camelid counterparts of all human IGKV and IGLV families and most IGHV families could be identified, based on canonical structure and sequence homology. Interestingly, this sequence homology seemed largely restricted to the Ig V genes and was far less apparent in other genes: 6 therapeutically relevant target genes differed significantly from their human orthologs. This contributed to efficient immunization of llamas with the human proteins CD70, MET, interleukin (IL)-1β and IL-6, resulting in large panels of functional antibodies. The in silico predicted human-homologous canonical folds of camelid-derived antibodies were confirmed by X-ray crystallography solving the structure of 2 selected camelid anti-CD70 and anti-MET antibodies. These antibodies showed identical fold combinations as found in the corresponding human germline V families, yielding binding site structures closely similar to those occurring in human antibodies. In conclusion, our results indicate that active immunization of camelids can be a powerful therapeutic antibody platform.

Resolving self-association of a therapeutic antibody by formulation optimization and molecular approaches

A common challenge encountered during development of high concentration monoclonal antibody formulations is preventing self-association. Depending on the antibody and its formulation, self-association can be seen as aggregation, precipitation, opalescence or phase separation. Here we report on an unusual manifestation of self-association, formation of a semi-solid gel or "gelation." Therapeutic monoclonal antibody C4 was isolated from human B cells based on its strong potency in neutralizing bacterial toxin in animal models. The purified antibody possessed the unusual property of forming a firm, opaque white gel when it was formulated at concentrations >30 mg/mL and the temperature was <6°C. Gel formation was reversible with temperature. Gelation was affected by salt concentration or pH, suggesting an electrostatic interaction between IgG monomers. A comparison of the C4 amino acid sequences to consensus germline sequences revealed differences in framework regions. A C4 variant in which the framework sequence was restored to the consensus germline sequence did not gel at 100 mg/mL at temperatures as low as 1°C. Additional genetic analysis was used to predict the key residue(s) involved in the gelation. Strikingly, a single substitution in the native antibody, replacing heavy chain glutamate 23 with lysine (E23K), was sufficient to prevent gelation. These results indicate that the framework region is involved in intermolecular interactions. The temperature dependence of gelation may be related to conformational changes near glutamate 23 or the regions it interacts with. Molecular engineering of the framework can be an effective approach to resolve the solubility issues of therapeutic antibodies.

Structural insight into activation of homoserine dehydrogenase from the archaeon Sulfolobus tokodaii via reduction

Homoserine dehydrogenase (HSD; 305 amino acid residues) catalyzes an NAD(P)-dependent reversible reaction between l-homoserine and aspartate 4-semialdehyde and is involved in the aspartate pathway. HSD from the hyperthermophilic archaeon Sulfolobus tokodaii was markedly activated (2.5-fold) by the addition of 0.8 mM dithiothreitol. The crystal structure of the homodimer indicated that the activation was caused by cleavage of the disulfide bond formed between two cysteine residues (C303) in the C-terminal regions of the two subunits.

Multiple binding modes of a small molecule to human Keap1 revealed by X-ray crystallography and molecular dynamics simulation

•

Keap1 is useful target for the design of drugs that regulate the response to oxidative stresses.

•

We determined two complex crystal structures of Keap1 with a small molecule ligand.

•

The ligand binds to Keap1 so as to mimic the physiological substrate Nrf2.

•

From molecular dynamics simulation results, the binding modes observed may be atypical in solution.

•

Key residues for ligand binding are common between crystal and MD structures.

Keap1 protein acts as a cellular sensor for oxidative stresses and regulates the transcription level of antioxidant genes through the ubiquitination of a corresponding transcription factor, Nrf2. A small molecule capable of binding to the Nrf2 interaction site of Keap1 could be a useful medicine. Here, we report two crystal structures, referred to as the soaking and the cocrystallization forms, of the Kelch domain of Keap1 with a small molecule, Ligand1. In these two forms, the Ligand1 molecule occupied the binding site of Keap1 so as to mimic the ETGE motif of Nrf2, although the mode of binding differed in the two forms. Because the Ligand1 molecule mediated the crystal packing in both the forms, the influence of crystal packing on the ligand binding was examined using a molecular dynamics (MD) simulation in aqueous conditions. In the MD structures from the soaking form, the ligand remained bound to Keap1 for over 20 ns, whereas the ligand tended to dissociate in the cocrystallization form. The MD structures could be classified into a few clusters that were related to but distinct from the crystal structures, indicating that the binding modes observed in crystals might be atypical of those in solution. However, the dominant ligand recognition residues in the crystal structures were commonly used in the MD structures to anchor the ligand. Therefore, the present structural information together with the MD simulation will be a useful basis for pharmaceutical drug development.

A word of caution about biological inference - Revisiting cysteine covalent state predictions

•

High prediction accuracy is often believed to validate implicit biological assumptions.

•

Cys redox state predictions assume a local sequence environmental effect.

•

Removing local sequence signals did not reduce prediction accuracy.

•

Cys redox state predictions apparently correlate with subcellular localization.

•

Subcellular localization depends on global sequence composition.

The success of methods for predicting the redox state of cysteine residues from the sequence environment seemed to validate the basic assumption that this state is mainly determined locally. However, the accuracy of predictions on randomized sequences or of non-cysteine residues remained high, suggesting that these predictions rather capture global features of proteins such as subcellular localization, which depends on composition. This illustrates that even high prediction accuracy is insufficient to validate implicit assumptions about a biological phenomenon. Correctly identifying the relevant underlying biochemical reasons for the success of a method is essential to gain proper biological insights and develop more accurate and novel bioinformatics tools.

Molecular mechanisms in the selective basal activation of pyrabactin receptor 1: Comparative analysis of mutants

Pyrabactin receptors (PYR) play a central role in abscisic acid (ABA) signal transduction; they are ABA receptors that inhibit type 2C protein phosphatases (PP2C). Molecular aspects contributing to increased basal activity of PYR against PP2C are studied by molecular dynamics (MD) simulations. An extensive series of MD simulations of the apo-form of mutagenized PYR1 as a homodimer and in complex with homology to ABA-insensitive 1 (HAB1) phosphatase are reported. In order to investigate the detailed molecular mechanisms mediating PYR1 activity, the MD data was analyzed by essential collective dynamics (ECD), a novel approach that allows the identification, with atomic resolution, of persistent dynamic correlations based on relatively short MD trajectories. Employing the ECD method, the effects of select mutations on the structure and dynamics of the PYR1 complexes were investigated and considered in the context of experimentally determined constitutive activities against HAB1. Approaches to rationally design constitutively active PYR1 constructs to increase PP2C inhibition are discussed.

Functional characterization of diverse ring-hydroxylating oxygenases and induction of complex aromatic catabolic gene clusters in Sphingobium sp. PNB

Sphingobium sp. PNB, like other sphingomonads, has multiple ring-hydroxylating oxygenase (RHO) genes. Three different fosmid clones have been sequenced to identify the putative genes responsible for the degradation of various aromatics in this bacterial strain. Comparison of the map of the catabolic genes with that of different sphingomonads revealed a similar arrangement of gene clusters that harbors seven sets of RHO terminal components and a sole set of electron transport (ET) proteins. The presence of distinctly conserved amino acid residues in ferredoxin and in silico molecular docking analyses of ferredoxin with the well characterized terminal oxygenase components indicated the structural uniqueness of the ET component in sphingomonads. The predicted substrate specificities, derived from the phylogenetic relationship of each of the RHOs, were examined based on transformation of putative substrates and their structural homologs by the recombinant strains expressing each of the oxygenases and the sole set of available ET proteins. The RHO AhdA1bA2b was functionally characterized for the first time and was found to be capable of transforming ethylbenzene, propylbenzene, cumene, p-cymene and biphenyl, in addition to a number of polycyclic aromatic hydrocarbons. Overexpression of aromatic catabolic genes in strain PNB, revealed by real-time PCR analyses, is a way forward to understand the complex regulation of degradative genes in sphingomonads.

Crystal structures of Phanerochaete chrysosporium pyranose 2-oxidase suggest that the N-terminus acts as a propeptide that assists in homotetramer assembly

The flavin-dependent homotetrameric enzyme pyranose 2-oxidase (P2O) is found mostly, but not exclusively, in lignocellulose-degrading fungi where it catalyzes the oxidation of β-d-glucose to the corresponding 2-keto sugar concomitantly with hydrogen peroxide formation during lignin solubilization. Here, we present crystal structures of P2O from the efficient lignocellulolytic basidiomycete Phanerochaete chrysosporium. Structures were determined of wild-type PcP2O from the natural fungal source, and two variants of recombinant full-length PcP2O, both in complex with the slow substrate 3-deoxy-3-fluoro-β-d-glucose. The active sites in PcP2O and P2O from Trametes multicolor (TmP2O) are highly conserved with identical substrate binding. Our structural analysis suggests that the 17 °C higher melting temperature of PcP2O compared to TmP2O is due to an increased number of intersubunit salt bridges. The structure of recombinant PcP2O expressed with its natural N-terminal sequence, including a proposed propeptide segment, reveals that the first five residues of the propeptide intercalate at the interface between A and B subunits to form stabilizing, mainly hydrophobic, interactions. In the structure of mature PcP2O purified from the natural source, the propeptide segment in subunit A has been replaced by a nearby loop in the B subunit. We propose that the propeptide in subunit A stabilizes the A/B interface of essential dimers in the homotetramer and that, upon maturation, it is replaced by the loop in the B subunit to form the mature subunit interface. This would imply that the propeptide segment of PcP2O acts as an intramolecular chaperone for oligomerization at the A/B interface of the essential dimer.

•

Structures of pyranose 2-oxidase from Phanerochaete chrysosporium were determined.

•

The N-terminus may act as a propeptide with a role in homotetramer assembly.

•

A large number of salt bridges between subunits provides thermostability.

•

The substrate is bound in the productive binding mode for oxidation at C2.

Enzymatic production of mono-ubiquitinated proteins for structural studies: The example of the Josephin domain of ataxin-3

Protein ubiquitination occurs through formation of an isopeptide bond between the C-terminal glycine of ubiquitin (Ub) and the ɛ-amino group of a substrate lysine residue. This post-translational modification, which occurs through the attachment of single and/or multiple copies of mono-ubiquitin and poly-ubiquitin chains, is involved in crucial cellular events such as protein degradation, cell-cycle regulation and DNA repair. The abnormal functioning of ubiquitin pathways is also implicated in the pathogenesis of several human diseases ranging from cancer to neurodegeneration. However, despite the undoubted biological importance, understanding the molecular basis of how ubiquitination regulates different pathways has up to now been strongly limited by the difficulty of producing the amounts of highly homogeneous samples that are needed for a structural characterization by X-ray crystallography and/or NMR. Here, we report on the production of milligrams of highly pure Josephin mono-ubiquitinated on lysine 117 through large scale in vitro enzymatic ubiquitination. Josephin is the catalytic domain of ataxin-3, a protein responsible for spinocerebellar ataxia type 3. Ataxin-3 is the first deubiquitinating enzyme (DUB) reported to be activated by mono-ubiquitination. We demonstrate that the samples produced with the described method are correctly folded and suitable for structural studies. The protocol allows facile selective labelling of the components. Our results provide an important proof-of-concept that may pave the way to new approaches to the in vitro study of ubiquitinated proteins.

•

We set up a protocol for large-scale in vitro enzymatic ubiqitination.

•

This produced milligrams of highly pure mono-ubiquitinated Josephin domain of ataxin-3.

•

We applied an alternative labelling scheme for the structural characterization of the sample by NMR.

•

Ubiquitin covalently linked on lysine 117 directly interacts with Josephin but does not alter the overall fold of the protein.

Contribution of histone N-terminal tails to the structure and stability of nucleosomes

Histones are the protein components of the nucleosome, which forms the basic architecture of eukaryotic chromatin. Histones H2A, H2B, H3, and H4 are composed of two common regions, the "histone fold" and the "histone tail". Many efforts have been focused on the mechanisms by which the post-translational modifications of histone tails regulate the higher-order chromatin architecture. On the other hand, previous biochemical studies have suggested that histone tails also affect the structure and stability of the nucleosome core particle itself. However, the precise contributions of each histone tail are unclear. In the present study, we determined the crystal structures of four mutant nucleosomes, in which one of the four histones, H2A, H2B, H3, or H4, lacked the N-terminal tail. We found that the deletion of the H2B or H3 N-terminal tail affected histone-DNA interactions and substantially decreased nucleosome stability. These findings provide important information for understanding the complex roles of histone tails in regulating chromatin structure.

Proposal for a new therapy for drug-resistant malaria using Plasmodium synthetic lethality inference

Many antimalarial drugs kill malaria parasites, but antimalarial drug resistance (ADR) and toxicity to normal cells limit their usefulness. To solve this problem, we suggest a new therapy for drug-resistant malaria. The approach consists of data integration and inference through homology analysis of yeast-human-Plasmodium. If one gene of a Plasmodium synthetic lethal (SL) gene pair has a mutation that causes ADR, a drug targeting the other gene of the SL pair might be used as an effective treatment for drug-resistant strains of malaria. A simple computational tool to analyze the inferred SL genes of Plasmodium species (malaria parasites Plasmodium falciparum and Plasmodium vivax for human malarial therapy, and rodent parasite Plasmodium berghei for in vivo studies of human malarias) was established to identify SL genes that can be used as drug targets. Information on SL gene pairs with ADR genes and their first neighbors was inferred from yeast SL genes to search for pertinent antimalarial drug targets. We not only suggest drug target gene candidates for further experimental validation, but also provide information on new usage for already-described drugs. The proposed specific antimalarial drug candidates can be inferred by searching drugs that cause a fitness defect in yeast SL genes.

MH2c: Characterization of major histocompatibility α-helices - an information criterion approach

Major histocompatibility proteins share a common overall structure or peptide binding groove. Two binding groove domains, on the same chain for major histocompatibility class I or on two different chains for major histocompatibility class II, contribute to that structure that consists of two α-helices ("wall") and a sheet of eight anti-parallel beta strands ("floor"). Apart from the peptide presented in the groove, the major histocompatibility α-helices play a central role for the interaction with the T cell receptor. This study presents a generalized mathematical approach for the characterization of these helices. We employed polynomials of degree 1 to 7 and splines with 1 to 2 nodes based on polynomials of degree 1 to 7 on the α-helices projected on their principal components. We evaluated all models with a corrected Akaike Information Criterion to determine which model represents the α-helices in the best way without overfitting the data. This method is applicable for both the stationary and the dynamic characterization of α-helices. By deriving differential geometric parameters from these models one obtains a reliable method to characterize and compare α-helices for a broad range of applications.

PROGRAM SUMMARY

Program title: MH²c (MH helix curves) Catalogue identifier: AELX_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AELX_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 327 565 No. of bytes in distributed program, including test data, etc.: 17 433 656 Distribution format: tar.gz Programming language: Matlab Computer: Personal computer architectures Operating system: Windows, Linux, Mac (all systems on which Matlab can be installed) RAM: Depends on the trajectory size, min. 1 GB (Matlab) Classification: 2.1, 4.9, 4.14 External routines: Curve Fitting Toolbox and Statistic Toolbox of Matlab Nature of problem: Major histocompatibility (MH) proteins share a similar overall structure. However, identical MH alleles which present different peptides differ by subtle conformational alterations. One hypothesis is that such conformational differences could be another level of T cell regulation. By this software package we present a reliable and systematic way to compare different MH structures to each other. Solution method: We tested several fitting approaches on all available experimental crystal structures of MH to obtain an overall picture of how to describe MH helices. For this purpose we transformed all complexes into the same space and applied splines and polynomials of several degrees to them. To draw a general conclusion which method fits them best we employed the "corrected Akaike Information Criterion". The software is applicable for all kinds of helices of biomolecules. Running time: Depends on the data, for a single stationary structure the runtime should not exceed a few seconds.

Structural bioinformatics analysis of free cysteines in protein environments

Cysteine has been considered as a “hydrophilic” amino acid because of its pK_a and its ability to form (weak) hydrogen bonds. However, cysteines are found mostly in hydrophobic environments, either in S–S (disulphide) form or in free cysteine form. When free cysteines are found on the surface of proteins, they are often involved in catalytic residues, as in cysteine proteases, P-loop phosphatases, etc. Additionally, a unique property of cysteines is that their side-chain volume is different from all other amino acids. This study is focused on the discrimination between structural versus active free cysteines based on a local environment analysis which does not appear to have been attempted previously. We have demonstrated the corresponding structural positions associated with free cysteines in their three-dimensional localization environment. We examined protein samples including nine, sequenced, coronavirus proteases and cysteine-rich non-membrane proteins. Our present study shows that the sequential environments of free cysteines of coronavirus proteases are rather hydrophobic and that the free cysteines of non-membrane proteases have a higher amount of contacts to hydrophobic residues and lower amount of contacts to polar or charged residues.