Search

Toll-like receptors (TLRs) play an important role in the innate immune response. While a lot is known about the structures of their extracellular parts, many questions are still left unanswered, when the structural basis of TLR activation is analyzed for the TLR intracellular domains. Here we report the structure and dynamics of TLR1 toll-interleukin like (TIR) cytoplasmic domain in crystal and in solution. We found that the TLR1-TIR domain is capable of specific binding of Zn with nanomolar affinity. Interactions with Zn are mediated by cysteine residues 667 and 686 and C667 is essential for the Zn binding. Potential structures of the TLR1-TIR/Zn complex were predicted in silico. Using the functional assays for the heterodimeric TLR1/2 receptor, we found that both Zn addition and Zn depletion affect the activity of TLR1, and C667A mutation disrupts the receptor activity. Analysis of C667 position in the TLR1 structure and possible effects of C667A mutation, suggests that zinc-binding ability of TLR1-TIR domain is critical for the receptor activation., Lushpa et al report the structure and dynamics of the TLR1 toll-interleukin like (TIR) cytoplasmic domain in both crystal and solution. They demonstrate that the TLR1 TIR domain is capable of specific binding of Zn with nanomolar affinity, which appears to be critical for receptor activation, and provide potential structures TLR1-TIR/Zn complex based on in silico data.

Neuronal calcium sensor-1 (NCS-1) is a four-EF-hand ubiquitous signaling protein modulating neuronal function and survival, which participates in neurodegeneration and carcinogenesis. NCS-1 recognizes specific sites on cellular membranes and regulates numerous targets, including G-protein coupled receptors and their kinases (GRKs). Here, with the use of cellular models and various biophysical and computational techniques, we demonstrate that NCS-1 is a redox-sensitive protein, which responds to oxidizing conditions by the formation of disulfide dimer (dNCS-1), involving its single, highly conservative cysteine C38. The dimer content is unaffected by the elevation of intracellular calcium levels but increases to 10–30% at high free zinc concentrations (characteristic of oxidative stress), which is accompanied by accumulation of the protein in punctual clusters in the perinuclear area. The formation of dNCS-1 represents a specific Zn2+-promoted process, requiring proper folding of the protein and occurring at redox potential values approaching apoptotic levels. The dimer binds Ca2+ only in one EF-hand per monomer, thereby representing a unique state, with decreased α-helicity and thermal stability, increased surface hydrophobicity, and markedly improved inhibitory activity against GRK1 due to 20-fold higher affinity towards the enzyme. Furthermore, dNCS-1 can coordinate zinc and, according to molecular modeling, has an asymmetrical structure and increased conformational flexibility of the subunits, which may underlie their enhanced target-binding properties. In HEK293 cells, dNCS-1 can be reduced by the thioredoxin system, otherwise accumulating as protein aggregates, which are degraded by the proteasome. Interestingly, NCS-1 silencing diminishes the susceptibility of Y79 cancer cells to oxidative stress-induced apoptosis, suggesting that NCS-1 may mediate redox-regulated pathways governing cell death/survival in response to oxidative conditions.

S100A1 is a member of the S100 family of small ubiquitous Ca2+-binding proteins, which participates in the regulation of cell differentiation, motility, and survival. It exists as homo- or heterodimers. S100A1 has also been shown to bind Zn2+, but the molecular mechanisms of this binding are not yet known. In this work, using ESI-MS and ITC, we demonstrate that S100A1 can coordinate 4 zinc ions per monomer, with two high affinity (KD~4 and 770 nm) and two low affinity sites. Using competitive binding experiments between Ca2+ and Zn2+ and QM/MM molecular modeling we conclude that Zn2+ high affinity sites are located in the EF-hand motifs of S100A1. In addition, two lower affinity sites can bind Zn2+ even when the EF-hands are saturated by Ca2+, resulting in a 2Ca2+:S100A1:2Zn2+ conformer. Finally, we show that, in contrast to calcium, an excess of Zn2+ produces a destabilizing effect on S100A1 structure and leads to its aggregation. We also determined a higher affinity to Ca2+ (KD~0.16 and 24 μm) than was previously reported for S100A1, which would allow this protein to function as a Ca2+/Zn2+-sensor both inside and outside cells, participating in diverse signaling pathways under normal and pathological conditions.

During evolution, viruses had to adapt to an increasingly complex environment of eukaryotic cells. Viral proteins that need to enter the cell nucleus or associate with nucleoli possess nuclear localization signals (NLSs) and nucleolar localization signals (NoLSs) for nuclear and nucleolar accumulation, respectively. As viral proteins are relatively small, acquisition of novel sequences seems to be a more complicated task for viruses than for eukaryotes. Here, we carried out a comprehensive analysis of the basic domain (BD) of HIV-1 Tat to show how viral proteins might evolve with NLSs and NoLSs without an increase in protein size. The HIV-1 Tat BD is involved in several functions, the most important being the transactivation of viral transcription. The BD also functions as an NLS, although it is substantially longer than a typical NLS. It seems that different regions in the BD could function as NLSs due to its enrichment with positively charged amino acids. Additionally, the high positive net charge inevitably causes the BD to function as an NoLS through a charge-specific mechanism. The integration of NLSs and NoLSs into functional domains enriched with positively charged amino acids might be a mechanism that allows the condensation of different functional sequences in small protein regions and, as a result, reduces protein size, influencing the origin and evolution of NLSs and NoLSs in viruses. IMPORTANCE Here, we investigated the molecular mechanism of nuclear localization signal (NLS) and nucleolar localization signal (NoLS) integration into the basic domain of HIV-1 Tat ((49)RKKRRQRRR(57)) and found that these two supplementary functions (i.e., function of NLS and function of NoLS) are embedded in the basic domain amino acid sequence. The integration of NLSs and NoLSs into functional domains of viral proteins enriched with positively charged amino acids is a mechanism that allows the concentration of different functions within small protein regions. Integration of NLS and NoLS into functional protein domains might have influenced the viral evolution, as this could prevent an increase in the protein size.

In this paper, we, for the first time, describe the interaction between the butyrylcholinesterase enzyme and echothiophate, a popular model compound and an analogue of the chemical warfare agents VX and VR, at the atomistic level. Competition between the two echothiophate conformations in the active site was found using molecular modeling techniques. The first one is close to the mode of binding of the substrates of choline series (butyrylcholine and butyrylthiocholine) and is inhibitory, since it is unable to react with the enzyme. The second one is characterized by a significantly worse estimated binding affinity and is reactive. Thus, echothiophate combines the features of two types of inhibitors: competitive and suicidal. This observation will help clarify the kinetic reaction scheme in order to accurately assess the kinetic constants, which is especially important when designing new butyrylcholinesterase variants capable of full-cycle hydrolysis of organophosphorus compounds.

Most proteins fold into 3D structures that determine how they function and orchestrate the biological processes of the cell. Recent developments in computational methods have led to protein structure predictions that have reached the accuracy of experimentally determined models. While this has been independently verified, the implementation of these methods across structural biology applications remains to be tested. Here, we evaluate the use of AlphaFold 2 (AF2) predictions in the study of characteristic structural elements; the impact of missense variants; function and ligand binding site predictions; modelling of interactions; and modelling of experimental structural data. For 11 proteomes, an average of 25% additional residues can be confidently modelled when compared to homology modelling, identifying structural features rarely seen in the PDB. AF2-based predictions of protein disorder and protein complexes surpass state-of-the-art tools and AF2 models can be used across diverse applications equally well compared to experimentally determined structures, when the confidence metrics are critically considered. In summary, we find that these advances are likely to have a transformative impact in structural biology and broader life science research.

Sexual reproduction is almost ubiquitous among extant eukaryotes. As most asexual lineages are short-lived, abandoning sex is commonly regarded as an evolutionary dead end. Still, putative anciently asexual lineages challenge this view. One of the most striking examples are bdelloid rotifers, microscopic freshwater invertebrates believed to have completely abandoned sexual reproduction tens of Myr ago. Here, we compare whole genomes of 11 wild-caught individuals of the bdelloid rotifer Adineta vaga and present evidence that some patterns in its genetic variation are incompatible with strict clonality and lack of genetic exchange. These patterns include genotype proportions close to Hardy-Weinberg expectations within loci, lack of linkage disequilibrium between distant loci, incongruent haplotype phylogenies across the genome, and evidence for hybridization between divergent lineages. Analysis of triallelic sites independently corroborates these findings. Our results provide evidence for interindividual genetic exchange and recombination in A. vaga, a species previously thought to be anciently asexual., Ancient, asexual lineages are rare as a lack of recombination is usually an evolutionary dead end. Here, authors compare complete genomes of 11 individual bdelloid rotifers that suggest evidence of regular genetic exchange between individuals in a species that was previously thought to be asexual.

Transactive response DNA and RNA binding protein 43 kDa (TDP-43) is a highly conserved heterogeneous nuclear ribonucleoprotein (hnRNP), which is involved in several steps of protein production including transcription and splicing. Its aggregates are frequently observed in motor neurons from amyotrophic lateral sclerosis patients and in the most common variant of frontotemporal lobar degeneration. Recently it was shown that TDP-43 is able to bind Zn2+ by its RRM domain. In this work, we have investigated Zn2+ binding to a short peptide 256&ndash, 264 from C-terminus of RRM2 domain using isothermal titration calorimetry, electrospray ionization mass spectrometry, QM/MM simulations, and NMR spectroscopy. We have found that this peptide is able to bind zinc ions with a Ka equal to 1.6 &times, 105 M&minus, 1. Our findings suggest the existence of a zinc binding site in the C-terminal region of RRM2 domain. Together with the existing structure of the RRM2 domain of TDP-43 we propose a model of its complex with Zn2+ which illustrates how zinc might regulate DNA/RNA binding.

The search for compounds that can inhibit the interaction of certain viral proteins with their cellular partners is a promising trend in the development of antiviral drugs. We have previously shown that binding of HIV-1 integrase with human Ku70 protein is essential for viral replication. Here, we present a novel, cheap, and fast assay to search for inhibitors of these proteins&rsquo, binding based on the usage of genetically encoded fluorescent tags linked to both integrase and Ku70. Using this approach, we have elucidated structure-activity relationships for a set of oligonucleotide conjugates with eosin and shown that their inhibitory activity is primarily achieved through interactions between the conjugate nucleic bases and integrase. Molecular modeling of HIV-1 integrase in complex with the conjugates suggests that they can shield E212/L213 residues in integrase, which are crucial for its efficient binding to Ku70, in a length-dependent manner. Using the developed system, we have found the 11-mer phosphorothioate bearing 3&rsquo, end eosin-Y to be the most efficient inhibitor among the tested conjugates.

Purpose Primary open-angle glaucoma (POAG) is a common ocular disease, associated with abnormalities in aqueous humor circulation and an increase in intraocular pressure (IOP), leading to progressive optical neuropathy and loss of vision. POAG pathogenesis includes alterations of the structural properties of the sclera, especially in the optic nerve head area, contributing to the degeneration of the retinal ganglion cells. Abnormal sclera biomechanics hinder adequate compensation of IOP fluctuations, thus aggravating POAG progression. The proteomic basis of biomechanical disorders in glaucomatous sclera remains poorly understood. This study is aimed at revealing alterations in major scleral proteins, associated with POAG, at different stages of the disease and with different IOP conditions. Methods Samples of sclera were collected from 67 patients with POAG during non-penetrating deep sclerectomy and from nine individuals without POAG. Scleral proteins were extracted with a strong lysis buffer, containing a combination of an ionic detergent, a chaotropic agent, and a disulfide reducing agent, and were separated using sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS–PAGE). The major scleral proteins were selected, subjected to in-gel digestion, and identified using matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF)/TOF mass spectrometry (MS), coupled with tandem mass spectrometry (MS/MS). The specific POAG-associated alterations of the selected proteins were analyzed with SDS–PAGE and confirmed with western blotting of the scleral extracts, using the respective antibodies. The group of POAG-associated proteins was analyzed using Gene Ontology and genome-wide association study enrichment and protein–protein interaction network prediction. Results A total of 11 proteins were identified, among which six proteins, namely, vimentin, angiopoietin-related protein 7, annexin A2, serum amyloid P component, serum albumin, and thrombospondin-4, were found to be upregulated in the sclera of patients with advanced and terminal POAG. In the early stages of the disease, thrombospondin-4 level was, on the contrary, reduced when compared with the control, whereas the concentration of vimentin varied, depending on the IOP level. Moreover, angiopoietin-related protein 7 manifested as two forms, exhibiting opposite behavior: The common 45 kDa form grew with the progression of POAG, whereas the 35 kDa (apparently non-glycosylated) form was absent in the control samples, appeared in patients with early POAG, and decreased in concentration over the course of the disease. Functional bioinformatics analysis linked the POAG-associated proteins with IOP alterations and predicted their secretion into extracellular space and their association with extracellular vesicles and a collagen-containing extracellular matrix. Conclusions POAG is accompanied by alterations of the scleral proteome, which represent a novel hallmark of the disease and can reflect pathological changes in scleral biochemistry and biomechanics. The potential mechanisms underlying these changes relate mainly to the structure of the extracellular matrix, protein glycosylation, and calcium binding, and may involve fibroblast cytoskeleton regulation, as well as oxidative and inflammatory responses.

Neuronal calcium sensors are a family of N-terminally myristoylated membrane-binding proteins possessing a different intracellular localization and thereby targeting unique signaling partner(s). Apart from the myristoyl group, the membrane attachment of these proteins may be modulated by their N-terminal positively charged residues responsible for specific recognition of the membrane components. Here, we examined the interaction of neuronal calcium sensor-1 (NCS-1) with natural membranes of different lipid composition as well as individual phospholipids in form of multilamellar liposomes or immobilized monolayers and characterized the role of myristoyl group and N-terminal lysine residues in membrane binding and phospholipid preference of the protein. NCS-1 binds to photoreceptor and hippocampal membranes in a Ca2+-independent manner and the binding is attenuated in the absence of myristoyl group. Meanwhile, the interaction with photoreceptor membranes is less dependent on myristoylation and more sensitive to replacement of K3, K7, and/or K9 of NCS-1 by glutamic acid, reflecting affinity of the protein to negatively charged phospholipids. Consistently, among the major phospholipids, NCS-1 preferentially interacts with phosphatidylserine and phosphatidylinositol with micromolar affinity and the interaction with the former is inhibited upon mutating of N-terminal lysines of the protein. Remarkably, NCS-1 demonstrates pronounced specific binding to phosphoinositides with high preference for phosphatidylinositol-3-phosphate. The binding does not depend on myristoylation and, unexpectedly, is not sensitive to the charge inversion mutations. Instead, phosphatidylinositol-3-phosphate can be recognized by a specific site located in the N-terminal region of the protein. These data provide important novel insights into the general mechanism of membrane binding of NCS-1 and its targeting to specific phospholipids ensuring involvement of the protein in phosphoinositide-regulated signaling pathways.

Proceedings of the National Academy of Sciences of the United States of America 117(37), 22841 - 22848 (2020). doi:10.1073/pnas.2010317117, Quantum mechanics/molecular mechanics (QM/MM) maturation of an immunoglobulin (Ig) powered by supercomputation delivers novel functionality to this catalytic template and facilitates artificial evolution of biocatalysts. We here employ density functional theory-based (DFT-b) tight binding and funnel metadynamics to advance our earlier QM/MM maturation of A17 Ig-paraoxonase (WTIgP) as a reactibody for organophosphorus toxins. It enables regulation of biocatalytic activity for tyrosine nucleophilic attack on phosphorus. The single amino acid substitution l-Leu47Lys results in 340-fold enhanced reactivity for paraoxon. The computed ground-state complex shows substrate-induced ionization of the nucleophilic l-Tyr37, now H-bonded to l-Lys47, resulting from repositioning of l-Lys47. Multiple antibody structural homologs, selected by phenylphosphonate covalent capture, show contrasting enantioselectivities for a P-chiral phenylphosphonate toxin. That is defined by crystallographic analysis of phenylphosphonylated reaction products for antibodies A5 and WTIgP. DFT-b analysis using QM regions based on these structures identifies transition states for the favored and disfavored reactions with surprising results. This stereoselection analysis is extended by funnel metadynamics to a range of WTIgP variants whose predicted stereoselectivity is endorsed by experimental analysis. The algorithms used here offer prospects for tailored design of highly evolved, genetically encoded organophosphorus scavengers and for broader functionalities of members of the Ig superfamily, including cell surface-exposed receptors., Published by National Acad. of Sciences, Washington, DC

BackgroundHIV-1 integration results in genomic DNA gaps that are repaired by cellular DNA repair pathways. This step of the lentiviral life cycle remains poorly understood despite its crucial importance for successful replication. We and others reported that Ku70 protein of the non-homologous end joining pathway (NHEJ) directly binds HIV-1 integrase (IN). Here, we studied the importance of this interaction for post-integrational gap repair and the recruitment of NHEJ factors in this process.ResultsWe engineered HIV-based pseudovirus with mutant IN defective in Ku70 binding and generated heterozygous Ku70, Ku80 and DNA-PKcs human knockout (KO) cells using CRISPR/Cas9. KO of either of these proteins or inhibition of DNA-PKcs catalytic activity substantially decreased the infectivity of HIV-1 with native IN but not with the mutant one. We used a recently developed qPCR assay for the measurement of gap repair efficiency to show that HIV-1 with mutant IN was defective in DNA post-integrational repair, whereas the wild type virus displayed such a defect only when NHEJ system was disrupted in any way. This effect was present in CRISPR/Cas9 modified 293T cells, in Jurkat and CEM lymphoid lines and in primary human PBMCs.ConclusionsOur data provide evidence that IN recruits DNA-PK to the site of HIV-1 post-integrational repair due to Ku70 binding—a novel finding that explains the involvement of DNA-PK despite the absence of free double stranded DNA breaks. In addition, our data clearly indicate the importance of interactions between HIV-1 IN and Ku70 in HIV-1 replication at the post-integrational repair step.

Science advances 6(26), eaaz9861 - (2020). doi:10.1126/sciadv.aaz9861, Microbial communities are self-controlled by repertoires of lethal agents, the antibiotics. In their turn, these antibiotics are regulated by bioscavengers that are selected in the course of evolution. Kinase-mediated phosphorylation represents one of the general strategies for the emergence of antibiotic resistance. A new subfamily of AmiN-like kinases, isolated from the Siberian bear microbiome, inactivates antibiotic amicoumacin by phosphorylation. The nanomolar substrate affinity defines AmiN as a phosphotransferase with a unique catalytic efficiency proximal to the diffusion limit. Crystallographic analysis and multiscale simulations revealed a catalytically perfect mechanism providing phosphorylation exclusively in the case of a closed active site that counteracts substrate promiscuity. AmiN kinase is a member of the previously unknown subfamily representing the first evidence of a specialized phosphotransferase bioscavenger., Published by Assoc., Washington, DC [u.a.]

Modeling tools provide a valuable support for DNA origami design. However, current solutions have limited application for conformational analysis of the designs. In this work we present a tool for a thorough study of DNA origami structure and dynamics. The tool is based on a novel coarse-grained model dedicated to geometry optimization and conformational analysis of DNA origami. We explored the ability of the model to predict dynamic behavior, global shapes, and fine details of two single-layer systems designed in hexagonal and square lattices using atomic force microscopy, Förster resonance energy transfer spectroscopy, and all-atom molecular dynamic simulations for validation of the results. We also examined the performance of the model for multilayer systems by simulation of DNA origami with published cryo-electron microscopy and atomic force microscopy structures. A good agreement between the simulated and experimental data makes the model suitable for conformational analysis of DNA origami objects. The tool is available at http://vsb.fbb.msu.ru/cosm as a web-service and as a standalone version.

Nucleic acid aptamers are prospective molecular recognizing elements. Similar to antibodies, aptamers are capable of providing specific recognition due to their spatial structure. However, the apparent simplicity of oligonucleotide folding is often elusive, as there is a balance between several conformations and, in some cases, oligomeric structures. This research is focused on establishing a thermodynamic background and the conformational heterogeneity of aptamers taking a series of thrombin DNA aptamers having G-quadruplex and duplex modules as an example. A series of aptamers with similar modular structures was characterized with spectroscopic and chromatographic techniques, providing examples of the conformational homogeneity of aptamers with high inhibitory activity, as well as a mixture of monomeric and oligomeric species for aptamers with low inhibitory activity. Thermodynamic parameters for aptamer unfolding were calculated, and their correlation with aptamer functional activity was found. Detailed analysis of thrombin complexes with G-quadruplex aptamers bound to exosite I revealed the similarity of the interfaces of aptamers with drastically different affinities to thrombin. It could be suggested that there are some events during complex formation that have a larger impact on the affinity than the states of initial and final macromolecules. Possible mechanisms of the complex formation and a role of the duplex module in the association process are discussed.

Hybrid QM/MM simulations augmented with enhanced sampling techniques proved to be advantageous in different usage scenarios, from studies of biological systems to drug and enzyme design. However, there are several factors that limit the applicability of the approach. First, typical biologically relevant systems are too large and hence computationally expensive for many QM methods. Second, a majority of fast non ab initio QM methods contain parameters for a very limited set of elements, which restrains their usage for applications involving radionuclides and other unusual compounds. Therefore, there is an incessant need for new tools which will expand both type and size of simulated objects. Here we present a novel combination of widely accepted molecular modelling packages GROMACS and MOPAC2012 and demonstrate its applicability for design of a catalytic antibody capable of organophosphorus compound hydrolysis.

Sexual reproduction which involves alternation of meiosis and syngamy is the ancestral condition of extant eukaryotes. Transitions to asexual reproduction were numerous, but most of the resulting eukaryotic lineages are rather short-lived. Still, there are several exceptions to this rule including darwinulid ostracods1,2and timema stick insects3. The most striking of them is bdelloid rotifers4–6, microscopic freshwater invertebrates which underwent an extensive adaptive radiation after apparently losing meiosis over 10 Mya. Indeed, both the lack of males in numerous bdelloid species and the lack of proper homology between chromosomes6rule out ordinary sex. However, this does not exclude the possibility of some other mode of interindividual genetic exchange and recombination in their populations7. Recent analyses based on a few loci suggested genetic exchanges in this group8,9, although this has been controversial10. Here, we compare complete genomes of 11 individuals from the wild population of the bdelloid rotiferAdineta vaga,and show that its genetic structure, which involves Hardy-Weinberg proportions of genotypes within loci and lack of linkage disequilibrium between distant loci, is incompatible with strictly clonal reproduction. Instead, it can emerge only under ongoing recombination between different individuals within this species, possibly through transformation. Such a genetic structure makes the population immune to negative long-term consequences of the loss of conventional meiosis11, although this does not necessarily imply that interindividual genetic exchanges inA. vagaare directly maintained by natural selection.

Butyrylcholinesterase (BChE) is considered as an efficient stoichiometric antidote against organophosphorus (OP) poisons. Recently we utilized combination of calculations and ultrahigh-throughput screening (uHTS) to select BChE variants capable of catalytic destruction of OP pesticide paraoxon. The purpose of this study was to elucidate the molecular mechanism underlying enzymatic hydrolysis of paraoxon by BChE variants using hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. Detailed analysis of accomplished QM/MM runs revealed that histidine residues introduced into the acyl-binding loop are always located in close proximity with aspartate residue at position 70. Histidine residue acts as general base thus leading to attacking water molecule activation and subsequent SN2 inline hydrolysis resulting in BChE reactivation. This combination resembles canonical catalytic triad found in active centers of various proteases. Carboxyl group activates histidine residue by altering its pKa, which in turn promotes the activation of water molecule in terms of its nucleophilicity. Observed re-protonation of catalytic serine residue at position 198 from histidine residue at position 438 recovers initial configuration of the enzyme’s active center, facilitating next catalytic cycle. We therefore suggest that utilization of uHTS platform in combination with deciphering of molecular mechanisms by QM/MM calculations may significantly improve our knowledge of enzyme function, propose new strategies for enzyme design and open new horizons in generation of catalytic bioscavengers against OP poisons.

In this study, we investigated the response of DNA-modified carbon nanotubes in the presence of thrombin and albumin. Interactions between DNA aptamers and carbon nanotubes were analyzed by Raman spectroscopy and by measuring the resistance of the construct after each step of biosensor assembly. It was shown that aptamer binding to carbon nanotubes decreased the G-band intensity and lowered the resistance of the sensor structure. The responses of the sensor exposed to thrombin and albumin differed substantially from each other, which provides a basis for the development of a highly selective biosensor. These results can be used in the development of a new generation of personal health monitoring systems.

Cyanea is a genus of large bloom-forming scyphozoans, including some of the most conspicuous rep- resentatives of megaplankton. Its taxonomy has been revised repeatedly throughout the last century due to the fact that most of the morphological characteristics ofCyanea species, such as color, structure of gastrovascular system and number of tentacles, may overlap greatly in different populations. Here, we report a new species ofCyanea,Cyaneatzetlinii sp. nov., from the White Sea, which is distinguishable from all previously describedCyanea species by an eye-spot-bearing bulb formed at the base of each rhopalium. This well-rec- ognizable morphological characteristic is supported at the molecular level by a substantial genetic distance in mito- chondrial (CO1: 9.6-10.6 %, 16S RNA: 3.1-3.5 %) as well as nuclear (ITS: 5.0 %, 18S RNA: 0.1 %) loci, making it the sister species to Cyanea capillata. Taking into account the young geological age of the White Sea and a substantial genetic divergence between C. tzetlinii sp. nov. and the nearest sister species, we suppose thatC.tzetlinii sp. nov. has been advected to the White Sea from elsewhere and may also inhabit other Arctic seas. Past ecological studies in the White Sea and possibly in other Arctic Seas could have conflatedC. tzetlinii sp. nov. with other species, which likely affected the analyses.

Disulfide bridges of bovine α-thrombin are studied using the analysis of the Raman spectral features in the spectral interval 500–550 cm−1. The changes under study are caused by the transitions from lyophilized state to the native solution in PBS with a protein concentration of 2.8 mM (100 mg/ml) and from native solution to solution with partially reduced disulfide bridges. The reduction takes place when dithiothreitol (DTT) is added to the sample.

Peptides are promising drug candidates due to high specificity and standout safety. Identification of bioactive peptides de novo using molecular docking is a widely used approach. However, current scoring functions are poorly optimized for peptide ligands. In this work, we present a novel algorithm PeptoGrid that rescores poses predicted by AutoDock Vina according to frequency information of ligand atoms with particular properties appearing at different positions in the target protein’s ligand binding site. We explored the relevance of PeptoGrid ranking with a virtual screening of peptide libraries using angiotensin-converting enzyme and GABAB receptor as targets. A reasonable agreement between the computational and experimental data suggests that PeptoGrid is suitable for discovering functional leads.