Search

Bioinformatics is an important research area for statisticians. This discussion provides some additional topics to the paper, namely on statistical contributions to detect differential expressed genes, for protein structure prediction, and for the analysis of highly correlated features in Glycomics datasets.

Polymerase δ‐interacting protein 2 (POLDIP2, PDIP38) is a multifaceted, “moonlighting” protein, involved in binding protein partners from many different cellular processes, including mitochondrial metabolism and DNA replication and repair. How POLDIP2 interacts with many different proteins is unknown. Towards this goal, we present the crystal structure of POLDIP2 to 2.8 Å, which exhibited a compact two‐domain β‐strand‐rich globular structure, confirmed by circular dichroism and small angle X‐ray scattering approaches. POLDIP2 comprised canonical DUF525 and YccV domains, but with a conserved domain linker packed tightly, resulting in an “extended” YccV module. A central channel was observed, which we hypothesize could influence structural changes potentially mediated by redox conditions, following observation of a modified cysteine residue in the channel. Unstructured regions were rebuilt by ab initio modelling to generate a model of full‐length POLDIP2. Molecular dynamics simulations revealed a highly dynamic N‐terminal region tethered to the YccV‐domain by an extended linker, potentially facilitating interactions with distal binding partners. Models of POLDIP2 complexed with two of its partners, PrimPol and PCNA, indicated that dynamic flexibility of the POLDIP2 N‐terminus and loop regions likely mediate protein interactions., PDB Code(s): 6Z9C

Improvements in cryo-trapping technology enable single-crystal time-resolved crystallography of an enzyme with 40 ms time resolution., Time-resolved crystallography of biomolecules in action has advanced rapidly as methods for serial crystallography have improved, but the large number of crystals and the complex experimental infrastructure that are required remain serious obstacles to its widespread application. Here, millisecond mix-and-quench crystallography (MMQX) has been developed, which yields millisecond time-resolved data using far fewer crystals and routine remote synchrotron data collection. To demonstrate the capabilities of MMQX, the conversion of oxaloacetic acid to phosphoenolpyruvate by phosphoenolpyruvate carboxy­kinase (PEPCK) is observed with a time resolution of 40 ms. By lowering the entry barrier to time-resolved crystallography, MMQX should enable a broad expansion in structural studies of protein dynamics.

Amyloid‐β peptide (Aβ) oligomers are pathogenic species of amyloid aggregates in Alzheimer's disease. Like certain protein toxins, Aβ oligomers permeabilize cellular membranes, presumably through a pore formation mechanism. Owing to their structural and stoichiometric heterogeneity, the structure of these pores remains to be characterized. We studied a functional Aβ42‐pore equivalent, created by fusing Aβ42 to the oligomerizing, soluble domain of the α‐hemolysin (αHL) toxin. Our data reveal Aβ42‐αHL oligomers to share major structural, functional, and biological properties with wild‐type Aβ42‐pores. Single‐particle cryo‐EM analysis of Aβ42‐αHL oligomers (with an overall 3.3 Å resolution) reveals the Aβ42‐pore region to be intrinsically flexible. The Aβ42‐αHL oligomers will allow many of the features of the wild‐type amyloid oligomers to be studied that cannot be otherwise, and may be a highly specific antigen for the development of immuno‐base diagnostics and therapies., A functional and structural Aβ42‐oligomer equivalent, created by fusing Aβ42 to the oligomerizing, soluble domain of the α‐hemolysin (αHL) toxin, shares major structural, functional, and biological properties with Alzheimer's wild‐type Aβ42 oligomers. Cryo‐EM studies reveal that the αHL‐displayed oligomers structurally remain in a well‐defined stoichiometry, providing a stable mimetic antigen for the development of conformation‐specific antibodies.

Summary Coordinating stress defense and plant growth is a survival strategy for adaptation to different environments that contains a series of processes, such as, cell growth, division and differentiation. However, little is known about the coordination mechanism for protein conformation change.A cyclophilin OsCYP20‐2 with a variant interacts with SLENDER RICE1 (SLR1) and OsFSD2 in the nucleus and chloroplasts, respectively, to integrate chilling tolerance and cell elongation in rice (Oryza sativa) (FSD2, Fe-superoxide dismutase 2).Mass spectrum assay showed that OsNuCYP20‐2 localized at the nucleus (nuclear located OsCYP20‐2) was a new variant of OsCYP20‐2 that truncated 71 amino‐acid residues in N‐terminal. The loss‐of function OsCYP20‐2 mutant showed sensitivity to chilling stress with accumulation of extra reactive oxygen species (ROS). In chloroplasts, the full‐length OsCYP20‐2 promotes OsFSD2 forming homodimers which enhance its activity, eliminating the accumulation of ROS under chilling stress. However, the mutant had shorter epidermal cells in comparison with wild‐type Hwayoung (HY). In the nucleus, OsCYP20‐2 caused conformation change of SLR1 to promote its degradation for cell elongation.Our data reveal a cyclophilin with a variant with dual‐localization in chloroplasts and the nucleus, which mediate chilling tolerance and cell elongation.

A method for determining the sample flow rate and concomitant light condition in time-resolved serial femtosecond crystallography is developed to analyze the intermediate-state structures of photosystem II., Photosystem II (PSII) catalyzes light-induced water oxidation through an Si-state cycle, leading to the generation of di-oxygen, protons and electrons. Pump–probe time-resolved serial femtosecond crystallography (TR-SFX) has been used to capture structural dynamics of light-sensitive proteins. In this approach, it is crucial to avoid light contamination in the samples when analyzing a particular reaction intermediate. Here, a method for determining a condition that avoids light contamination of the PSII microcrystals while minimizing sample consumption in TR-SFX is described. By swapping the pump and probe pulses with a very short delay between them, the structural changes that occur during the S1-to-S2 transition were examined and a boundary of the excitation region was accurately determined. With the sample flow rate and concomitant illumination conditions determined, the S2-state structure of PSII could be analyzed at room temperature, revealing the structural changes that occur during the S1-to-S2 transition at ambient temperature. Though the structure of the manganese cluster was similar to previous studies, the behaviors of the water molecules in the two channels (O1 and O4 channels) were found to be different. By comparing with the previous studies performed at low temperature or with a different delay time, the possible channels for water inlet and structural changes important for the water-splitting reaction were revealed.

The pandemic outbreak of coronavirus disease 2019 (COVID‐19) across the world has led to millions of infection cases and caused a global public health crisis. Current research suggests that SARS‐CoV‐2 is a highly contagious coronavirus that spreads rapidly through communities. To understand the mechanisms of viral replication, it is imperative to investigate coronavirus viral replicase, a huge protein complex comprising up to 16 viral nonstructural and associated host proteins, which is the most promising antiviral target for inhibiting viral genome replication and transcription. Recently, several components of the viral replicase complex in SARS‐CoV‐2 have been solved to provide a basis for the design of new antiviral therapeutics. Here, we report the crystal structure of the SARS‐CoV‐2 nsp7+8 tetramer, which comprises two copies of each protein representing nsp7's full‐length and the C‐terminus of nsp8 owing to N‐terminus proteolysis during the process of crystallization. We also identified a long helical extension and highly flexible N‐terminal domain of nsp8, which is preferred for interacting with single‐stranded nucleic acids., PDB Code(s): 7DCD

Irradiation of proteins with intense X-ray radiation generates specific structural and chemical alterations, such as metal-ion reduction, reduction of disulfide bonds and decarboxylation reactions. An analysis of a set of disulfide bonds in known structures suggests that the most susceptible disulfide bonds feature a carbonyl O atom positioned along the extension of the S—S bond vector, forming a stabilizing S⋯O interaction, which upon irradiation leads to polarization of the disulfide bond and its eventual reduction., Radiation-induced damage to protein crystals during X-ray diffraction data collection is a major impediment to obtaining accurate structural information on macromolecules. Some of the specific impairments that are inflicted upon highly brilliant X-ray irradiation are metal-ion reduction, disulfide-bond cleavage and a loss of the integrity of the carboxyl groups of acidic residues. With respect to disulfide-bond reduction, previous results have indicated that not all disulfide bridges are equally susceptible to damage. A careful analysis of the chemical environment of disulfide bonds in the structures of elastase, lysozyme, acetylcholinesterase and other proteins suggests that S—S bonds which engage in a close contact with a carbonyl O atom along the extension of the S—S bond vector are more susceptible to reduction than the others. Such an arrangement predisposes electron transfer to occur from the O atom to the disulfide bond, leading to its reduction. The interaction between a nucleophile and an electrophile, akin to hydrogen bonding, stabilizes protein structures, but it also provides a pathway of electron transfer to the S—S bond, leading to its reduction during exposure of the protein crystal to an intense X-ray beam. An otherwise stabilizing interaction can thus be the cause of destabilization under the condition of radiation exposure.

The crystal structures of HigA3 antitoxin and DNA-bound HigA3 allow detailed structural insights into the HigBA toxin–antitoxin system of Mycobacterium tuberculosis. The unique dimerization of HigA3 induces DNA bending., The bacterial toxin–antitoxin (TA) system regulates cell growth under various environmental stresses. Mycobacterium tuberculosis, the causative pathogen of tuberculosis (TB), has three HigBA type II TA systems with reverse gene organization, consisting of the toxin protein HigB and labile antitoxin protein HigA. Most type II TA modules are transcriptionally autoregulated by the antitoxin itself. In this report, we first present the crystal structure of the M. tuberculosis HigA3 antitoxin (MtHigA3) and MtHigA3 bound to its operator DNA complex. We also investigated the interaction between MtHigA3 and DNA using NMR spectroscopy. The MtHigA3 antitoxin structure is a homodimer that contains a structurally well conserved DNA-binding domain at the N-terminus and a dimerization domain at the C-terminus. Upon comparing the HigA homologue structures, a distinct difference was found in the C-terminal region that possesses the β-lid, and diverse orientations of two helix–turn–helix (HTH) motifs from HigA homologue dimers were observed. The structure of MtHigA3 bound to DNA reveals that the promoter DNA is bound to two HTH motifs of the MtHigA3 dimer presenting 46.5° bending, and the distance between the two HTH motifs of each MtHigA3 monomer was increased in MtHigA3 bound to DNA. The β-lid, which is found only in the tertiary structure of MtHigA3 among the HigA homologues, causes the formation of a tight dimerization network and leads to a unique arrangement for dimer formation that is related to the curvature of the bound DNA. This work could contribute to the understanding of the HigBA system of M. tuberculosis at the atomic level and may contribute to the development of new antibiotics for TB treatment.

Motivation Halides are negatively charged ions of halogens, forming fluorides (F−), chlorides (Cl−), bromides (Br−) and iodides (I−). These anions are quite reactive and interact both specifically and non-specifically with proteins. Despite their ubiquitous presence and important roles in protein function, little is known about the preferences of halides binding to proteins. To address this problem, we performed the analysis of halide–protein interactions, based on the entries in the Protein Data Bank. Results We have compiled a pipeline for the quick analysis of halide-binding sites in proteins using the available software. Our analysis revealed that all of halides are strongly attracted by the guanidinium moiety of arginine side chains, however, there are also certain preferences among halides for other partners. Furthermore, there is a certain preference for coordination numbers in the binding sites, with a correlation between coordination numbers and amino acid composition. This pipeline can be used as a tool for the analysis of specific halide–protein interactions and assist phasing experiments relying on halides as anomalous scatters. Availability and implementation All data described in this article can be reproduced via complied pipeline published at https://github.com/rostkick/Halide_sites/blob/master/README.md. Supplementary information Supplementary data are available at Bioinformatics online.

The relationship between GTP binding and hydrolysis on the one hand and septin filament assembly and function on the other has yet to be fully elucidated. The crystal structures of SEPT3 subgroup members bound to either GDP or to a GTP analogue shed light on this question by indicating how the squeezing of one of the inter-subunit interfaces may lead to the exposure of a polybasic region important for membrane association., Human septins 3, 9 and 12 are the only members of a specific subgroup of septins that display several unusual features, including the absence of a C-terminal coiled coil. This particular subgroup (the SEPT3 septins) are present in rod-like octameric protofilaments but are lacking in similar hexameric assemblies, which only contain representatives of the three remaining subgroups. Both hexamers and octamers can self-assemble into mixed filaments by end-to-end association, implying that the SEPT3 septins may facilitate polymerization but not necessarily function. These filaments frequently associate into higher order complexes which associate with biological membranes, triggering a wide range of cellular events. In the present work, a complete compendium of crystal structures for the GTP-binding domains of all of the SEPT3 subgroup members when bound to either GDP or to a GTP analogue is provided. The structures reveal a unique degree of plasticity at one of the filamentous interfaces (dubbed NC). Specifically, structures of the GDP and GTPγS complexes of SEPT9 reveal a squeezing mechanism at the NC interface which would expel a polybasic region from its binding site and render it free to interact with negatively charged membranes. On the other hand, a polyacidic region associated with helix α5′, the orientation of which is particular to this subgroup, provides a safe haven for the polybasic region when retracted within the interface. Together, these results suggest a mechanism which couples GTP binding and hydrolysis to membrane association and implies a unique role for the SEPT3 subgroup in this process. These observations can be accounted for by constellations of specific amino-acid residues that are found only in this subgroup and by the absence of the C-terminal coiled coil. Such conclusions can only be reached owing to the completeness of the structural studies presented here.

Confocal Raman microspectral imaging was adopted to elucidate the cellular drug responses of osteosarcoma cells (OC) to N-[N-(3, 5-difluorophenyl acetyl)-L-alanyl]-sphenylglycine butyl ester (DAPT), a γ-secretase inhibitor, by identifying the drug induced subcellular compositional and structural changes. Methods: Spectral information were acquired from cultured osteosarcoma cells treated with 0 (Untreated Group, UT), 10 (10 μM DAPT treated, 10T), 20 μM (20 μM DAPT treated, 20T) DAPT for 24 hours. A one-way ANOVA and Tukey's honest significant difference (HSD) post hoc multiple test were sequentially applied to address spectral features among three groups. Multivariate algorithms such as K-means clustering analysis (KCA) and Principal component analysis (PCA) were used to highlight the structural and compositional differences, while, univariate imaging was applied to illustrate the distribution pattern of certain cellular components after drug treatment. Results: Major biochemical changes in DAPT-induced apoptosis came from changes in the content and structure of proteins, lipids, and nucleic acids. By adopted multivariate algorithms, the drug induced cellular changes was identified by the morphology and spectral characteristics between untreated cells and treated cells, testified that DAPT mainly acted in the nuclear region. With the increase of the drug concentration, the content of main subcellular compositions, such nucleic acid, protein, and lipid decreased. In an addition, DAPT-induced nuclear fragmentation and apoptosis was depicted by the univariate Raman image of major cellular components (nucleic acids, proteins and lipids). Conclusions: The achieved Raman spectral and imaging results illustrated detailed DAPT-induced subcellular compositional and structural variations as a function of drug dose. Such observations can not only explain drug therapeutic mechanisms of OC DAPT treatment, and also provide new insights for accessing the medicine curative efficacy and predicting prognosis.

The structure of human MICAL3 provides information about its mechanism., MICAL is an oxidoreductase that participates in cytoskeleton reorganization via actin disassembly in the presence of NADPH. Although three MICALs (MICAL1, MICAL2 and MICAL3) have been identified in mammals, only the structure of mouse MICAL1 has been reported. Here, the first crystal structure of human MICAL3, which contains the flavin-containing monooxygenase (FMO) and calponin-homology (CH) domains, is reported. MICAL3 has an FAD/NADP-binding Rossmann-fold domain for mono­oxygenase activity like MICAL1. The FMO and CH domains of both MICAL3 and MICAL1 are highly similar in structure, but superimposition of the two structures shows a different relative position of the CH domain in the asymmetric unit. Based on kinetic analyses, the catalytic efficiency of MICAL3 dramatically increased on adding F-actin only when the CH domain was available. However, this did not occur when two residues, Glu213 and Arg530, were mutated in the FMO and CH domains, respectively. Overall, MICAL3 is structurally highly similar to MICAL1, which suggests that they may adopt the same catalytic mechanism, but the difference in the relative position of the CH domain produces a difference in F-actin substrate specificity.

MiR-10b, miR-335, and miR-21 are classes of microRNAs (miRNAs) that are overexpressed in breast cancer. Thus, in our study we aimed to test the hypothesis that miRNAs may have direct interactions with proteins and the possibility to inhibit/activate the functional site of proteins and enzymes. For this purpose, we choose three miRNAs involved in breast cancer to study interactions between some proteins and genes, including BRCA1 and PTEN, by processing the docking and matching tools using the Hex8 and HADDOCK server. Mathematically, the hidden Markov models were created by using MATLAB script according to the algorithm in order to study and validate the interactions and bonds between proteins and miRNAs. The main results demonstrate the ability of miR-10b, miR-335, and miR-21 to create direct interactions with 3D protein structures. Furthermore, these results may lead to another pathway of research, i.e. the direct interaction between proteins and their sub-units, to highlight the data obtained previously and demonstrate that proteins may directly interact with ncRNA instead of mRNA. Moreover, our study suggests developing research on different pathways of association proteins-miRNAs as a part of epigenetic extra-nuclear regulation. Taken together, our study provides the first evidence of direct interactions between miRNAs and proteins.

A method is presented for efficient co-crystal structure determination for G protein-coupled receptors taking advantage of serial femtosecond crystallography., Rational structure-based drug design (SBDD) relies on the availability of a large number of co-crystal structures to map the ligand-binding pocket of the target protein and use this information for lead-compound optimization via an iterative process. While SBDD has proven successful for many drug-discovery projects, its application to G protein-coupled receptors (GPCRs) has been limited owing to extreme difficulties with their crystallization. Here, a method is presented for the rapid determination of multiple co-crystal structures for a target GPCR in complex with various ligands, taking advantage of the serial femtosecond crystallography approach, which obviates the need for large crystals and requires only submilligram quantities of purified protein. The method was applied to the human β2-adrenergic receptor, resulting in eight room-temperature co-crystal structures with six different ligands, including previously unreported structures with carvedilol and propranolol. The generality of the proposed method was tested with three other receptors. This approach has the potential to enable SBDD for GPCRs and other difficult-to-crystallize membrane proteins.

Human eyes absent (EYA) proteins possess Tyr phosphatase activity which is critical for numerous cancer and metastasis promoting activities, making it an attractive target for cancer therapy. In this work, we demonstrate that the inhibitor-bound form of EYA2 does not favour binding to Mg2+ , which is indispensable for the Tyr phosphatase activity. We further describe characterization and optimization of this class of allosteric inhibitors. A series of analogues were synthesized to improve potency of the inhibitors and to elucidate structure-activity relationships. Two co-crystal structures confirm the binding modes of this class of inhibitors. Our medicinal chemical, structural, biochemical and biophysical studies provide insight into the molecular interactions of EYA2 with these allosteric inhibitors. The compounds derived from this study are useful for exploring the function of the Tyr phosphatase activity of EYA2 in normal and cancerous cells and serve as reference compounds for screening or developing allosteric phosphatase inhibitors. Finally, the co-crystal structures reported in this study will aid in structure-based drug discovery against EYA2. This article is protected by copyright. All rights reserved.

JLP-G and ALP were supported by the ERDF/Spanish Ministry of Science, Innovation and Universities-State Research Agency (Grant RTI2018-096246-B-I00) and Consejeria de Economia, Conocimiento, Empresas y Universidad, Junta de Andalucia (Grants P11-CTS-7187 and P18-RT-2413) . NM-T was supported by Aula FUNCANIS-UGR. ES was supported by the ERDF/Spanish Ministry of Science, Innovation and Universities-State Research Agency (Grant SAF2015-69796) . Access to an EU_FT-ICR_MS network installation was funded by the EU Horizon 2020 grant 731077. EA-C and MM were supported by the Spanish Ministry of Science and Innovation-State Research Agency (Grant PID2019-103901 GB-I00) and Gobierno de Aragon-FEDER (Grant E35_20R) . Support of the BioCeV center (CZ.1.05/1.1.00/02.0109) and the CMS/CIISB facility (MEYS CZ-LM2018127) is also gratefully acknowledged. ANN was supported by grants BT/PR26099/BID/7/811/2017 from Department of Biotechnology (DBT, India) and MTR/2019/000392 from Science, Engineering and Research Board (SERB, India) ., The multifunctional nature of human flavoproteins is critically linked to their ability to populate multiple conformational states. Ligand binding, post-translational modifications and disease-associated mutations can reshape this functional landscape, although the structure-function relationships of these effects are not well understood. Herein, we characterized the structural and functional consequences of two mutations (the cancer associated P187S and the phosphomimetic S82D) on different ligation states which are relevant to flavin binding, intracellular stability and catalysis of the disease-associated NQO1 flavoprotein. We found that these mutations affected the stability locally and their effects propagated differently through the protein structure depending both on the nature of the mutation and the ligand bound, showing directional preference from the mutated site and leading to specific phenotypic manifestations in different functional traits (FAD binding, catalysis and inhibition, intracellular stability and pharmacological response to ligands). Our study thus supports that pleitropic effects of disease-causing mutations and phosphorylation events on human flavoproteins may be caused by longrange structural propagation of stability effects to different functional sites that depend on the ligation-state and site-specific perturbations. Our approach can be of general application to investigate these pleiotropic effects at the flavoproteome scale in the absence of high-resolution structural models., ERDF/Spanish Ministry of Science, Innovation and Universities-State Research Agency RTI2018-096246-B-I00- SAF2015-69796, Junta de Andalucia P11-CTS-7187- P18-RT-2413, Aula FUNCANIS-UGR, European Commission 731077, Spanish Ministry of Science and Innovation-State Research Agency PID2019-103901 GB-I00, Gobierno de Aragon-FEDER E35_20R, BioCeV center CZ.1.05/1.1.00/02.0109, CMS/CIISB facility MEYS CZ-LM2018127, Department of Biotechnology (DBT) India BT/PR26099/BID/7/811/2017, Science, Engineering and Research Board (SERB, India) MTR/2019/000392

The application of double electron‐electron resonance (DEER) with site‐directed spin labeling (SDSL) to measure distances in proteins and protein complexes in living cells puts rigorous restraints on the spin‐label. The linkage and paramagnetic centers need to resist the reducing conditions of the cell. Rigid attachment of the probe to the protein improves precision of the measured distances. Here, three two‐armed GdIII complexes, GdIII‐CLaNP13a/b/c were synthesized. Rather than the disulfide linkage of most other CLaNP molecules, a thioether linkage was used to avoid reductive dissociation of the linker. The doubly GdIII labeled N55C/V57C/K147C/T151C variants of T4Lysozyme were measured by 95 GHz DEER. The constructs were measured in vitro, in cell lysate and in Dictyostelium discoideum cells. Measured distances were 4.5 nm, consistent with results from paramagnetic NMR. A narrow distance distribution and typical modulation depth, also in cell, indicate complete and durable labeling and probe rigidity due to the dual attachment sites., Two armed and ready! A two armed, thus rigid, GdIII‐CLaNP13 probe for double electron–electron resonance (DEER) measurements on proteins can be used to measure nanometer distances in cells with high precision.

Values of S2CH and S2NH order parameters derived from NMR relaxation measurements on proteins cannot be used straightforwardly to determine protein structure because they cannot be related to a single protein structure, but are defined in terms of an average over a conformational ensemble. Molecular dynamics simulation can generate a conformational ensemble and thus can be used to restrain S2CH and S2NH order parameters towards experimentally derived target values S2CH (exp) and S2NH (exp). Application of S2CH and S2NH order-parameter restraining MD simulation to bond vectors in 63 side chains of the protein hen egg white lysozyme using 51 S2CH (exp) target values and 28 S2NH (exp) target values shows that a conformational ensemble compatible with the experimentally derived data can be obtained by using this technique. It is observed that S2CH order-parameter restraining of C-H bonds in methyl groups is less reliable than S2NH order-parameter restraining because of the possibly less valid assumptions and approximations used to derive experimental S2CH (exp) values from NMR relaxation measurements and the necessity to adopt the assumption of uniform rotational motion of methyl C-H bonds around their symmetry axis and of the independence of these motions from each other. The restrained simulations demonstrate that side chains on the protein surface are highly dynamic. Any hydrogen bonds they form and that appear in any of four different crystal structures, are fluctuating with short lifetimes in solution., ChemBioChem, 22 (6), ISSN:1439-4227, ISSN:1439-7633

Predicting protein structure from the amino acid sequence has been a challenge with theoretical and practical significance in biophysics. Despite the recent progresses elicited by improved inter‐residue contact prediction, contact‐based structure prediction has gradually reached the performance ceiling. New methods have been proposed to predict the inter‐residue distance, but unanimously by simplifying the real‐valued distance prediction into a multiclass classification problem. Here, a lightweight regression‐based distance prediction method is shown, which adopts the generative adversarial network to capture the delicate geometric relationship between residue pairs and thus could predict the continuous, real‐valued inter‐residue distance rapidly and satisfactorily. The predicted residue distance map allows quick structure modeling by the CNS suite, and the constructed models approach the same level of quality as the other state‐of‐the‐art protein structure prediction methods when tested on CASP13 targets. Moreover, this method can be used directly for the structure prediction of membrane proteins without transfer learning., A lightweight GAN‐based protein inter‐residue distance prediction method is reported. This method is capable of capturing the delicate geometric relationship between residue pairs and thus could predict the continuous, real‐valued distance rapidly and satisfactorily. Pipelined with structure modeling software like the CNS suite, this method could construct protein structure models with the state‐of‐the‐art quality.

The main protease (3CLp) of the SARS‐CoV‐2, the causative agent for the COVID‐19 pandemic, is one of the main targets for drug development. To be active, 3CLp relies on a complex interplay between dimerization, active site flexibility, and allosteric regulation. The deciphering of these mechanisms is a crucial step to enable the search for inhibitors. In this context, using NMR spectroscopy, we studied the conformation of dimeric 3CLp from the SARS‐CoV‐2 and monitored ligand binding, based on NMR signal assignments. We performed a fragment‐based screening that led to the identification of 38 fragment hits. Their binding sites showed three hotspots on 3CLp, two in the substrate binding pocket and one at the dimer interface. F01 is a non‐covalent inhibitor of the 3CLp and has antiviral activity in SARS‐CoV‐2 infected cells. This study sheds light on the complex structure‐function relationships of 3CLp and constitutes a strong basis to assist in developing potent 3CLp inhibitors., We report the liquid‐sate NMR spectroscopy analysis of the dimeric SARS‐CoV‐2 main protease (3CLp), including its backbone assignments, to study its complex conformational regulation. Using fragment‐based NMR screening, we highlighted three hotspots on the protein, two in the substrate binding pocket and one at the dimer interface, and we identified a non‐covalent reversible inhibitor of 3CLp that has antiviral activity in infected cells.

Motivation The investigation of the structure of biological systems at the molecular level gives insights about their functions and dynamics. Shape and surface of biomolecules are fundamental to molecular recognition events. Characterizing their geometry can lead to more adequate predictions of their interactions. In the present work, we assess the performance of reference shape retrieval methods from the computer vision community on protein shapes. Results Shape retrieval methods are efficient in identifying orthologous proteins and tracking large conformational changes. This work illustrates the interest for the protein surface shape as a higher-level representation of the protein structure that (i) abstracts the underlying protein sequence, structure or fold, (ii) allows the use of shape retrieval methods to screen large databases of protein structures to identify surficial homologs and possible interacting partners and (iii) opens an extension of the protein structure–function paradigm toward a protein structure-surface(s)-function paradigm. Availabilityand implementation All data are available online at http://datasetmachat.drugdesign.fr. Supplementary information Supplementary data are available at Bioinformatics online.

Dicer is a member of the ribonuclease III enzyme family and processes double-stranded RNA into small functional RNAs. The variation in the domain architecture of Dicer among different species whilst preserving its biological dicing function is intriguing. Here, we describe the structure and function of a novel catalytically active RNase III protein, a non-canonical Dicer (PsDCR1), found in budding yeast Pichia stipitis. The structure of the catalytically active region (the catalytic RNase III domain and double-stranded RNA-binding domain 1 [dsRBD1]) of DCR1 showed that RNaseIII domain is structurally similar to yeast RNase III (Rnt1p) but uniquely presents dsRBD1 in a diagonal orientation, forming a catalytic core made of homodimer and large RNA-binding surface. The second dsRNA binding domain at C-terminus, which is absent in Rnt1, enhances the RNA cleavage activity. Although the cleavage pattern of PsDCR1 anchors an apical loop similar to Rnt1, the cleavage activity depended on the sequence motif at the lower stem, not the apical loop, of hairpin RNA. Through RNA sequencing and RNA mutations, we showed that RNA cleavage by PsDCR1 is determined by the stem-loop structure of the RNA substrate, suggesting the possibility that stem-loop RNA-guided gene silencing pathway exists in budding yeast.

The COVID-19 main protease (Mpro), one of the conserved proteins of the novel coronavirus is crucial for its replication and so is a very lucrative drug target. Till now, there is no drug molecule that has been convincingly identified as the inhibitor of the function of this protein. The current pandemic situation demands a shortcut to quickly reach to a lead compound or a drug, which may not be the best but might serve as an interim solution at least. Following this notion, the present investigation uses virtual screening to find a molecule which is alraedy approved as a drug for some other disease but could be repurposed to inhibit Mpro. The potential of the present method of work to identify such a molecule, which otherwise would have been missed out, lies in the fact that instead of just using the crystallographically identified conformation of the receptor’s ligand binding pocket, molecular dynamics generated ensemble of conformations has been used. It implicitly included the possibilities of “induced-fit” and/or “population shift” mechanisms of ligand fitting. As a result, the investigation has not only identified antiviral drugs like ribavirin, ritonavir, etc., but it has also captured a wide variety of drugs for various other diseases like amrubicin, cangrelor, desmopressin, diosmin, etc. as the potent possibilities. Some of these ligands are versatile to form stable interactions with various different conformations of the receptor and therefore have been statistically surfaced in the investigation. Overall the investigation offers a wide range of compounds for further testing to confirm their scopes of applications to combat the COVID-19 pandemic. Supplementary Information The online version contains supplementary material available at 10.1007/s00894-021-04732-1.

Ependymin-related proteins fold into a dimer using a subunit formed by stacking antiparallel β-sheets and contain a hydrophobic pocket that is found in the similar folds of bacterial VioE and LolA/LolB. As in VioE and LolA/LolB, this hydrophobic pocket may be used to bind hydrophobic molecules., Ependymin was first discovered as a predominant protein in brain extracellular fluid in fish and was suggested to be involved in functions mostly related to learning and memory. Orthologous proteins to ependymin called ependymin-related proteins (EPDRs) have been found to exist in various tissues from sea urchins to humans, yet their functional role remains to be revealed. In this study, the structures of EPDR1 from frog, mouse and human were determined and analyzed. All of the EPDR1s fold into a dimer using a monomeric subunit that is mostly made up of two stacking antiparallel β-sheets with a curvature on one side, resulting in the formation of a deep hydrophobic pocket. All six of the cysteine residues in the monomeric subunit participate in the formation of three intramolecular disulfide bonds. Other interesting features of EPDR1 include two asparagine residues with glycosylation and a Ca2+-binding site. The EPDR1 fold is very similar to the folds of bacterial VioE and LolA/LolB, which also use a similar hydrophobic pocket for their respective functions as a hydrophobic substrate-binding enzyme and a lipoprotein carrier, respectively. A further fatty-acid binding assay using EPDR1 suggests that it indeed binds to fatty acids, presumably via this pocket. Additional interactome analysis of EPDR1 showed that EPDR1 interacts with insulin-like growth factor 2 receptor and flotillin proteins, which are known to be involved in protein and vesicle translocation.

The tryptophan synthases from three human pathogens show remarkable structural conservation, but at the same time display local differences in both their catalytic and allosteric sites that may be responsible for the observed differences in catalysis and inhibitor binding. This functional dissimilarity may be exploited in the design of species-specific enzyme inhibitors., Tryptophan biosynthesis is one of the most characterized processes in bacteria, in which the enzymes from Salmonella typhimurium and Escherichia coli serve as model systems. Tryptophan synthase (TrpAB) catalyzes the final two steps of tryptophan biosynthesis in plants, fungi and bacteria. This pyridoxal 5′-phosphate (PLP)-dependent enzyme consists of two protein chains, α (TrpA) and β (TrpB), functioning as a linear αββα heterotetrameric complex containing two TrpAB units. The reaction has a complicated, multistep mechanism resulting in the β-replacement of the hydroxyl group of l-serine with an indole moiety. Recent studies have shown that functional TrpAB is required for the survival of pathogenic bacteria in macrophages and for evading host defense. Therefore, TrpAB is a promising target for drug discovery, as its orthologs include enzymes from the important human pathogens Streptococcus pneumoniae, Legionella pneumophila and Francisella tularensis, the causative agents of pneumonia, legionnaires’ disease and tularemia, respectively. However, specific biochemical and structural properties of the TrpABs from these organisms have not been investigated. To fill the important phylogenetic gaps in the understanding of TrpABs and to uncover unique features of TrpAB orthologs to spearhead future drug-discovery efforts, the TrpABs from L. pneumophila, F. tularensis and S. pneumoniae have been characterized. In addition to kinetic properties and inhibitor-sensitivity data, structural information gathered using X-ray crystallo­graphy is presented. The enzymes show remarkable structural conservation, but at the same time display local differences in both their catalytic and allosteric sites that may be responsible for the observed differences in catalysis and inhibitor binding. This functional dissimilarity may be exploited in the design of species-specific enzyme inhibitors.

Background: Hydrogen Sulfide (H2S), a third member of gasotransmitter family along with nitric oxide (NO) and carbon monoxide (CO), exerts a wide range of cellular and molecular actions in our body. There is a large body of evidence suggesting that H2S plays an important role in cancer metastasis; however, the molecular mechanisms of H2S-mediated acceleration of cancer metastasis remain unknown. Methods: We examined the promote effects of H2S on phenotype of gastric cancer (GC) cells (including those of express wild type CD36 and mutant CD36) in vitro and in vivo. GC patients' samples were used for clinical translational significance evaluation. Findings: H2S triggered lipid metabolism reprogramming by significantly up-regulating the expression of the fatty-acid receptor CD36 (CD36) and directly activating CD36 in GC cells. Mechanistically, a disulfide bond located between cysteine (Cys)333 and Cys272 within the CD36 protein structure that was labile to H2S-mediated modification. The long chain-fatty acid (LC-FA) binding pocket was capped by a turn in the CD36 protein, located between helical and sheet structures that were stabilized by the Cys333-Cys272. This limited the secondary binding between LC-FAs and lysine (Lys)334. Breaking the Cys333-Cys272 disulfide bond restored the second LC-FA binding conformation of CD36. Targeting CD36 in vivo blocked H2S-promoted metastasis and improved animal survival. Interpretation: These findings identify that the Cys333-Cys272 disulfide bond disrupted the integrity of the second LC-FA binding conformation of CD36. Therefore, CD36 can directly activate LC-FA access to the cytoplasm by acting as a direct target molecule for H2S. Keywords: H2S, CD36, Disulfide bond, Long chain-fatty acid, Nrf2, Metastasis