Distance-dependent hydrophobic-hydrophobic contacts in protein folding simulations

A primer for junior trainees: Recognition of RNA modifications by RNA-binding proteins

The complexity of RNA cannot be fully expressed with the canonical A, C, G, and U alphabet. To date, over 170 distinct chemical modifications to RNA have been discovered in living systems. RNA modifications can profoundly impact the cellular outcomes of messenger RNAs (mRNAs), transfer and ribosomal RNAs, and noncoding RNAs. Additionally, aberrant RNA modifications are associated with human disease. RNA modifications are a rising topic within the fields of biochemistry and molecular biology. The role of RNA modifications in gene regulation, disease pathogenesis, and therapeutic applications increasingly captures the attention of the scientific community. This review aims to provide undergraduates, junior trainees, and educators with an appreciation for the significance of RNA modifications in eukaryotic organisms, alongside the skills required to identify and analyze fundamental RNA-protein interactions. The pumilio RNA-binding protein and YT521-B homology (YTH) family of modified RNA-binding proteins serve as examples to highlight the fundamental biochemical interactions that underlie the specific recognition of both unmodified and modified ribonucleotides, respectively. By instilling these foundational, textbook concepts through practical examples, this review contributes an analytical toolkit that facilitates engagement with RNA modifications research at large.

Cimicifugoside H-2 as an Inhibitor of IKK1/Alpha: A Molecular Docking and Dynamic Simulation Study

One of the most challenging issues scientists face is finding a suitable non-invasive treatment for cancer, as it is widespread around the world. The efficacy of phytochemicals that target oncogenic pathways appears to be quite promising and has gained attention over the past few years. We investigated the effect of docking phytochemicals isolated from the rhizomes of the Cimicifuga foetida plant on different domains of the IκB kinase alpha (IKK1/alpha) protein. The Cimicifugoside H-2 phytochemical registered a high docking score on the activation loop of IKK1/alpha amongst the other phytochemicals compared to the positive control. The interaction of the protein with Cimicifugoside H-2 was mostly stabilized by hydrogen bonds and hydrophobic interactions. A dynamic simulation was then performed with the Cimicifugoside H-2 phytochemical on the activation loop of IKK1/alpha, revealing that Cimicifugoside H-2 is a possible inhibitor of this protein. The pharmacokinetic properties of the drug were also examined to assess the safety of administering the drug. Therefore, in this in silico study, we discovered that the Cimicifugoside H-2 phytochemical inhibits the actively mutated conformation of IKK1/alpha, potentially suppressing the nuclear factor kappa light chain enhancer of activated B cells (NF-κB) pathway.

Electrochemical and structural characterization of recombinant respiratory proteins of the acidophilic iron oxidizer Ferrovum sp. PN-J47-F6 suggests adaptations to the acidic pH at protein level

The tendency of the periplasmic redox proteins in acidophiles to have more positive redox potentials (Em) than their homologous counterparts in neutrophiles suggests an adaptation to acidic pH at protein level, since thermodynamics of electron transfer processes are also affected by acidic pH. Since this conclusion is mainly based on the electrochemical characterization of redox proteins from extreme acidophiles of the genus Acidithiobacillus, we aimed to characterize three recombinant redox proteins of the more moderate acidophile Ferrovum sp. PN-J47-F6. We applied protein film voltammetry and linear sweep voltammetry coupled to UV/Vis spectroscopy to characterize the redox behavior of HiPIP-41, CytC-18, and CytC-78, respectively. The Em-values of HiPIP-41 (571 ± 16 mV), CytC-18 (276 ± 8 mV, 416 ± 2 mV), and CytC-78 (308 ± 7 mV, 399 ± 7 mV) were indeed more positive than those of homologous redox proteins in neutrophiles. Moreover, our findings suggest that the adaptation of redox proteins with respect to their Em occurs more gradually in response to the pH, since there are also differences between moderate and more extreme acidophiles. In order to address structure function correlations in these redox proteins with respect to structural features affecting the Em, we conducted a comparative structural analysis of the Ferrovum-derived redox proteins and homologs of Acidithiobacillus spp. and neutrophilic proteobacteria. Hydrophobic contacts in the redox cofactor binding pockets resulting in a low solvent accessibility appear to be the major factor contributing to the more positive Em-values in acidophile-derived redox proteins. While additional cysteines in HiPIPs of acidophiles might increase the effective shielding of the [4Fe-4S]-cofactor, the tight shielding of the heme centers in acidophile-derived cytochromes is achieved by a drastic increase in hydrophobic contacts (A.f. Cyc41), and by a larger fraction of aromatic residues in the binding pockets (CytC-18, CytC-78).

Identification of the flavivirus conserved residues in the envelope protein hinge region for the rational design of a candidate West Nile live-attenuated vaccine

The flavivirus envelope protein is a class II fusion protein that drives flavivirus-cell membrane fusion. The membrane fusion process is triggered by the conformational change of the E protein from dimer in the virion to trimer, which involves the rearrangement of three domains, EDI, EDII, and EDIII. The movement between EDI and EDII initiates the formation of the E protein trimer. The EDI-EDII hinge region utilizes four motifs to exert the hinge effect at the interdomain region and is crucial for the membrane fusion activity of the E protein. Using West Nile virus (WNV) NY99 strain derived from an infectious clone, we investigated the role of eight flavivirus-conserved hydrophobic residues in the EDI-EDII hinge region in the conformational change of E protein from dimer to trimer and viral entry. Single mutations of the E-A54, E-I130, E-I135, E-I196, and E-Y201 residues affected infectivity. Importantly, the E-A54I and E-Y201P mutations fully attenuated the mouse neuroinvasive phenotype of WNV. The results suggest that multiple flavivirus-conserved hydrophobic residues in the EDI-EDII hinge region play a critical role in the structure-function of the E protein and some contribute to the virulence phenotype of flaviviruses as demonstrated by the attenuation of the mouse neuroinvasive phenotype of WNV. Thus, as a proof of concept, residues in the EDI-EDII hinge region are proposed targets to engineer attenuating mutations for inclusion in the rational design of candidate live-attenuated flavivirus vaccines.

Length and Charge of the N-terminus Regulate the Lifetime of the Signaling State of Photoactive Yellow Protein

Photoactive yellow protein (PYP) is one of the most extensively studied photoreceptors. Nevertheless, the role of the N-terminus in the photocycle and structural transitions is still elusive. Here, we attached additional amino acids to the N-terminus of PYP and investigated the effect of the length and charge of additional N-terminal residues using circular dichroism, two-dimensional nuclear magnetic resonance (2D-NMR), transient absorption (TA), and transient grating (TG) spectroscopic techniques. TA experiments showed that, except for negatively charged residues (5D-PYP), additional N-terminal residues of PYP generally enable faster dark recovery from the putative signaling state (pB2) to the ground state (pG). TG data showed that although the degree of structural changes can be controlled by adjusting specific amino acid residues in the extended N-terminus of N-terminal extended PYPs (NE-PYPs), the dark recovery times of wt-PYP and NE-PYPs, except for 5D-PYP, are independent of the structural differences between pG and pB2 states. These results demonstrate that the recovery time and the degree of structural change can be regulated by controlling the length and sequence of N-terminal residues of PYP. The findings in this study emphasize the need for careful attention to the remaining amino acid residues when designing recombinant proteins for genetic engineering purposes.

Rinmaker: a fast, versatile and reliable tool to determine residue interaction networks in proteins

BACKGROUND

Residue Interaction Networks (RINs) map the crystallographic description of a protein into a graph, where amino acids are represented as nodes and non-covalent bonds as edges. Determination and visualization of a protein as a RIN provides insights on the topological properties (and hence their related biological functions) of large proteins without dealing with the full complexity of the three-dimensional description, and hence it represents an invaluable tool of modern bioinformatics.

RESULTS

We present RINmaker, a fast, flexible, and powerful tool for determining and visualizing RINs that include all standard non-covalent interactions. RINmaker is offered as a cross-platform and open source software that can be used either as a command-line tool or through a web application or a web API service. We benchmark its efficiency against the main alternatives and provide explicit tests to show its performance and its correctness.

CONCLUSIONS

RINmaker is designed to be fully customizable, from a simple and handy support for experimental research to a sophisticated computational tool that can be embedded into a large computational pipeline. Hence, it paves the way to bridge the gap between data-driven/machine learning approaches and numerical simulations of simple, physically motivated, models.

Three-helix bundle and SH3-type barrels: autonomously stable structural motifs in small and large proteins

In this study, we investigated two variants of a three-helix bundle and SH3-type barrel, compact in space, present in small and large proteins of various living organisms. Using a neural graph network, proteins with three-helix bundle (n = 1377) and SH3-type barrels (n = 1914) spatial folds were selected. Molecular experiments were performed for small proteins with these folds, and motifs were studied autonomously outside the protein environment at 300, 340, and 370 K. A comparative analysis of the main parameters of the structures in the course of the experiment was performed, including gyration radius, area accessible to the solvent, number of hydrophobic and hydrogen bonds, and root-mean-square deviation of atomic positions (RMSD). We exhibited an autonomous stability of the studied folds outside the protein environment in an aquatic medium. We aimed to demonstrate the possibility of analyzing three-helix bundle and SH3-type barrels autonomously outside the protein globule, thereby reducing the computational time and increasing performance without significant loss of information.Communicated by Ramaswamy H. Sarma.

Divergent evolution drives high diversity of toll-like receptors (TLRs) in passerine birds: Buntings and finches

Toll-like receptors (TLRs) form a key component of animal innate immunity, being responsible for recognition of conserved microbial structures. As such, TLRs may be subject to diversifying and balancing selection, which maintains allelic variation both within and between populations. However, most research on TLRs in non-model avian species is focused on bottlenecked populations with depleted genetic variation. Here, we assessed variation at the extracellular domains of three TLR genes (TLR1LA, TLR3, TLR4) across eleven species from two passerine families of buntings (Emberizidae) and finches (Fringillidae), all having large breeding population sizes (millions of individuals). We found extraordinary TLR polymorphism in our study taxa, with >100 alleles detected at TLR1LA and TLR4 across species and high haplotype diversity (>0.75) in several species. Despite recent species divergence, no nucleotide allelic variants were shared between species, suggesting rapid TLR evolution. Higher variation at TLR1LA and TLR4 than TLR3 was associated with a stronger signal of diversifying selection, as measured with nucleotide substitutions rates and the number of positively selected sites (PSS). Structural protein modelling of TLRs showed that some PSS detected within TLR1LA and TLR4 were previously recognized as functionally important sites or were located in their proximity, possibly affecting ligand recognition. Furthermore, we identified PSS responsible for major surface electrostatic charge clustering, which may indicate their adaptive importance. Our study provides compelling evidence for the divergent evolution of TLR genes in buntings and finches and indicates that high TLR variation may be adaptively maintained via diversifying selection acting on functional ligand binding sites.

Conformational Dynamics of the Activated GLP-1 Receptor-Gs Complex Revealed by Cross-Linking Mass Spectrometry and Integrative Structure Modeling

Despite advances in characterizing the structures and functions of G protein-coupled receptors (GPCRs), our understanding of GPCR activation and signaling is still limited by the lack of information on conformational dynamics. It is particularly challenging to study the dynamics of GPCR complexes with their signaling partners because of their transient nature and low stability. Here, by combining cross-linking mass spectrometry (CLMS) with integrative structure modeling, we map the conformational ensemble of an activated GPCR-G protein complex at near-atomic resolution. The integrative structures describe heterogeneous conformations for a high number of potential alternative active states of the GLP-1 receptor-G_s complex. These structures show marked differences from the previously determined cryo-EM structure, especially at the receptor-G_s interface and in the interior of the G_s heterotrimer. Alanine-scanning mutagenesis coupled with pharmacological assays validates the functional significance of 24 interface residue contacts only observed in the integrative structures, yet absent in the cryo-EM structure. Through the integration of spatial connectivity data from CLMS with structure modeling, our study provides a new approach that is generalizable to characterizing the conformational dynamics of GPCR signaling complexes.

Simultaneous functioning of different light-harvesting complexes-a strategy of adaptation of purple bacterium Rhodopseudomonas palustris to low illumination conditions

Novel peripheral light-harvesting (LH) complex designated as LL LH2 was isolated along with LH4 complex from Rhodopseudomonas palustris cells grown under low light intensity (LL). FPLC-MS/MS allowed to reveal PucABd and PucBabc apoproteins in LL LH2 complex, which is different from previously described LH4 complex containing PucABd, PucABa and PucBb. The main carotenoids in LL LH2 complex were rhodopin and 3,4-didehydrorhodopin. Three-dimensional modeling demonstrated which amino acid residues of all the β-subunits could interact with carotenoids (Car) and bacteriochlorophyll a (BChl a). Analysis of amino acid sequences of α-subunits of both LL complexes showed presence of different C-terminal motifs, IESSVNVG in αa subunit and IESSIKAV in αd subunit, in the same positions of C-termini, which could reflect different retention force of LL LH2 and LH4 on hydroxyl apatite, facilitating successful isolation of these complexes. Differences of these LL complexes in protein and carotenoid composition, in efficiency of energy transfer from Car to BChl a, which is two times lower in LL LH2 than in LH4, allow to assign it to a novel type of light-harvesting complex in Rhodopseudomonas palustris.

Dynamic Evolution of Avian RNA Virus Sensors: Repeated Loss of RIG-I and RIPLET

Retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated protein 5 (MDA5) are key RNA virus sensors belonging to the RIG-I-like receptor (RLR) family. The activation of the RLR inflammasome leads to the establishment of antiviral state, mainly through interferon-mediated signaling. The evolutionary dynamics of RLRs has been studied mainly in mammals, where rare cases of RLR gene losses were described. By in silico screening of avian genomes, we previously described two independent disruptions of MDA5 in two bird orders. Here, we extend this analysis to approximately 150 avian genomes and report 16 independent evolutionary events of RIG-I inactivation. Interestingly, in almost all cases, these inactivations are coupled with genetic disruptions of RIPLET/RNF135, an ubiquitin ligase RIG-I regulator. Complete absence of any detectable RIG-I sequences is unique to several galliform species, including the domestic chicken (Gallus gallus). We further aimed to determine compensatory evolution of MDA5 in RIG-I-deficient species. While we were unable to show any specific global pattern of adaptive evolution in RIG-I-deficient species, in galliforms, the analyses of positive selection and surface charge distribution support the hypothesis of some compensatory evolution in MDA5 after RIG-I loss. This work highlights the dynamic nature of evolution in bird RNA virus sensors.

Identification of the Shigella flexneri Wzy Domain Modulating WzzpHS-2 Interaction and Detection of the Wzy/Wzz/Oag Complex

Shigella flexneri implements the Wzy-dependent pathway to biosynthesize the O antigen (Oag) component of its surface lipopolysaccharide. The inner membrane polymerase Wzy_SF catalyzes the repeat addition of undecaprenol-diphosphate-linked Oag (Und-PP-RUs) to produce a polysaccharide, the length of which is tightly regulated by two competing copolymerase proteins, Wzz_SF (short-type Oag; 10 to 17 RUs) and Wzz_pHS-2 (very-long-type Oag; >90 RUs). The nature of the interaction between Wzy_SF and Wzz_SF/Wzz_pHS-2 in Oag polymerization remains poorly characterized, with the majority of the literature characterizing the individual protein constituents of the Wzy-dependent pathway. Here, we report instead a major investigation into the specific binding interactions of Wzy_SF with its copolymerase counterparts. For the first time, a region of Wzy_SF that forms a unique binding site for Wzz_pHS-2 has been identified. Specifically, this work has elucidated key Wzy_SF moieties at the N- and C-terminal domains (NTD and CTD) that form an intramolecular pocket modulating the Wzz_pHS-2 interaction. Novel copurification data highlight that disruption of residues within this NTD-CTD pocket impairs the interaction with Wzz_pHS-2 without affecting Wzz_SF binding, thereby specifically disrupting polymerization of longer polysaccharide chains. This study provides a novel understanding of the molecular interaction of Wzy_SF with Wzz_SF/Wzz_pHS-2 in the Wzy-dependent pathway and, furthermore, detects the Wzy/Wzz/Und-PP-Oag complex for the first time. Beyond S. flexneri, this work may be extended to provide insight into the interactions between protein homologues expressed by related species, especially members of Enterobacteriaceae, that produce dual Oag chain length determinants. IMPORTANCE Shigella flexneri is a pathogen causing significant morbidity and mortality, predominantly devastating the pediatric age group in developing countries. A major virulence factor contributing to S. flexneri pathogenesis is its surface lipopolysaccharide, which is comprised of three domains: lipid A, core oligosaccharide, and O antigen (Oag). The Wzy-dependent pathway is the most common biosynthetic mechanism implemented for Oag biosynthesis by Gram-negative bacteria, including S. flexneri. The nature of the interaction between the polymerase, Wzy_SF, and the polysaccharide copolymerases, Wzz_SF and Wzz_pHS-2, in Oag polymerization is poorly characterized. This study investigates the molecular interplay between Wzy_SF and its copolymerases, deciphering key interactions in the Wzy-dependent pathway that may be extended beyond S. flexneri, providing insight into Oag biosynthesis in Gram-negative bacteria.

Proline/Glycine residues of the PG-levels guide conformational changes along the transport cycle in the mitochondrial carnitine/acylcarnitine carrier (SLC25A20)

Mitochondrial carnitine/acylcarnitine carrier (CAC) is a member of the mitochondrial carrier (MC) family and imports acylcarnitine into the mitochondrial matrix in exchange for carnitine, playing a pivotal role in carnitine shuttle, crucial for fatty acid oxidation. The crystallized structure of CAC has not been solved yet, however, the availability of several in vitro/in silico studies, also based on the crystallized structures of the ADP/ATP carrier in the cytosolic-conformation and in the matrix-conformation, has made possible to confirm the hypothesis of the single-binding centered-gated pore mechanism for all the members of the MC family. In addition, our recent bioinformatics analyses allowed quantifying in silico the importance of protein residues of MC substrate binding region, of those involved in the formation of the matrix and cytosolic gates, and of those belonging to the Pro/Gly (PG) levels, proposed to be crucial for the tilting/kinking/bending of the six MC transmembrane helices, funneling the substrate translocation pathway. Here we present a combined in silico/in vitro analysis employed for investigating the role played by a group of 6 proline residues and 6 glycine residues, highly conserved in CAC, belonging to MC PG-levels. Residues of the PG-levels surround the similarly located MC common substrate binding region, and were proposed to lead conformational changes and substrate translocation, following substrate binding. For our analysis, we employed 3D molecular modeling approaches, alanine scanning site-directed mutagenesis and in vitro transport assays. Our analysis reveals that P130 (H3), G268 (H6) and G220 (H5), mutated in alanine, affect severely CAC transport activity (mutant catalytic efficiency lower than 5 % compared to the wild type CAC), most likely due to their major role in triggering CAC conformational changes, following carnitine binding. Notably, P30A (H1) and G121A (H3) CAC mutants, increase the carnitine uptake up to 217 % and 112 %, respectively, compared to the wild type CAC.

Profiling the Site of Protein CoAlation and Coenzyme A Stabilization Interactions

Coenzyme A (CoA) is a key cellular metabolite known for its diverse functions in metabolism and regulation of gene expression. CoA was recently shown to play an important antioxidant role under various cellular stress conditions by forming a disulfide bond with proteins, termed CoAlation. Using anti-CoA antibodies and liquid chromatography tandem mass spectrometry (LC-MS/MS) methodologies, CoAlated proteins were identified from various organisms/tissues/cell-lines under stress conditions. In this study, we integrated currently known CoAlated proteins into mammalian and bacterial datasets (CoAlomes), resulting in a total of 2093 CoAlated proteins (2862 CoAlation sites). Functional classification of these proteins showed that CoAlation is widespread among proteins involved in cellular metabolism, stress response and protein synthesis. Using 35 published CoAlated protein structures, we studied the stabilization interactions of each CoA segment (adenosine diphosphate (ADP) moiety and pantetheine tail) within the microenvironment of the modified cysteines. Alternating polar-non-polar residues, positively charged residues and hydrophobic interactions mainly stabilize the pantetheine tail, phosphate groups and the ADP moiety, respectively. A flexible nature of CoA is observed in examined structures, allowing it to adapt its conformation through interactions with residues surrounding the CoAlation site. Based on these findings, we propose three modes of CoA binding to proteins. Overall, this study summarizes currently available knowledge on CoAlated proteins, their functional distribution and CoA-protein stabilization interactions.

PIF - A Java library for finding atomic interactions and extracting geometric features supporting the analysis of protein structures

Proteins play an essential role in the functioning of living organisms. The enormity of the atomic interactions in proteins is essential in controlling their spatial structures and dynamics. It can also provide scientists with valuable information that help to determine the native structures of proteins. This paper presents the PIF (Protein Interaction Finder) library for the Java language, enabling the identification of selected atomic interactions (hydrogen and disulfide bonds, ionic, hydrophobic, aromatic-aromatic, sulfur-aromatic, and amino-aromatic interactions) based on the three-dimensional structure of proteins. The interaction calculation rules applied in PIF rely on documented theoretical foundations gathered from experimental studies of interactions in native protein structures. The library has a universal purpose, supporting drug discovery and development processes and protein structure modeling. Finding the atomic interactions can also deliver numerical features for various Artificial Intelligence (AI) models built for protein analysis. The conducted research comparing the results obtained with the use of the PIF library and competing tools has shown that our solution can effectively determine the interactions occurring in protein structures for entire collections of proteins. Moreover, as a solution that provides a programming interface, the PIF library can be used in any Java project, making it a universal tool.

Personalized Medicine in Mitochondrial Health and Disease: Molecular Basis of Therapeutic Approaches Based on Nutritional Supplements and Their Analogs

Mitochondrial diseases (MDs) may result from mutations affecting nuclear or mitochondrial genes, encoding mitochondrial proteins, or non-protein-coding mitochondrial RNA. Despite the great variability of affected genes, in the most severe cases, a neuromuscular and neurodegenerative phenotype is observed, and no specific therapy exists for a complete recovery from the disease. The most used treatments are symptomatic and based on the administration of antioxidant cocktails combined with antiepileptic/antipsychotic drugs and supportive therapy for multiorgan involvement. Nevertheless, the real utility of antioxidant cocktail treatments for patients affected by MDs still needs to be scientifically demonstrated. Unfortunately, clinical trials for antioxidant therapies using α-tocopherol, ascorbate, glutathione, riboflavin, niacin, acetyl-carnitine and coenzyme Q have met a limited success. Indeed, it would be expected that the employed antioxidants can only be effective if they are able to target the specific mechanism, i.e., involving the central and peripheral nervous system, responsible for the clinical manifestations of the disease. Noteworthily, very often the phenotypes characterizing MD patients are associated with mutations in proteins whose function does not depend on specific cofactors. Conversely, the administration of the antioxidant cocktails might determine the suppression of endogenous oxidants resulting in deleterious effects on cell viability and/or toxicity for patients. In order to avoid toxicity effects and before administering the antioxidant therapy, it might be useful to ascertain the blood serum levels of antioxidants and cofactors to be administered in MD patients. It would be also worthwhile to check the localization of mutations affecting proteins whose function should depend (less or more directly) on the cofactors to be administered, for estimating the real need and predicting the success of the proposed cofactor/antioxidant-based therapy.

Modeling SARS-CoV-2 spike/ACE2 protein-protein interactions for predicting the binding affinity of new spike variants for ACE2, and novel ACE2 structurally related human protein targets, for COVID-19 handling in the 3PM context

Aims

The rapid spread of new SARS-CoV-2 variants has highlighted the crucial role played in the infection by mutations occurring at the SARS-CoV-2 spike receptor binding domain (RBD) in the interactions with the human ACE2 receptor. In this context, it urgently needs to develop new rapid tools for quickly predicting the affinity of ACE2 for the SARS-CoV-2 spike RBD protein variants to be used with the ongoing SARS-CoV-2 genomic sequencing activities in the clinics, aiming to gain clues about the transmissibility and virulence of new variants, to prevent new outbreaks and to quickly estimate the severity of the disease in the context of the 3PM.

Methods

In our study, we used a computational pipeline for calculating the interaction energies at the SARS-CoV-2 spike RBD/ACE2 protein-protein interface for a selected group of characterized infectious variants of concern/interest (VoC/VoI). By using our pipeline, we built 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for the VoC B.1.1.7-United Kingdom (carrying the mutations of concern/interest N501Y, S494P, E484K at the RBD), P.1-Japan/Brazil (RBD mutations: K417T, E484K, N501Y), B.1.351-South Africa (RBD mutations: K417N, E484K, N501Y), B.1.427/B.1.429-California (RBD mutations: L452R), the B.1.141 (RBD mutations: N439K), and the recent B.1.617.1-India (RBD mutations: L452R; E484Q) and the B.1.620 (RBD mutations: S477N; E484K). Then, we used the obtained 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for predicting the interaction energies at the protein-protein interface.

Results

Along SARS-CoV-2 mutation database screening and mutation localization analysis, it was ascertained that the most dangerous mutations at VoC/VoI spike proteins are located mainly at three regions of the SARS-CoV-2 spike "boat-shaped" receptor binding motif, on the RBD domain. Notably, the P.1 Japan/Brazil variant present three mutations, K417T, E484K, N501Y, located along the entire receptor binding motif, which apparently determines the highest interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein-protein interface, among those calculated. Conversely, it was also observed that the replacement of a single acidic/hydrophilic residue with a basic residue (E484K or N439K) at the "stern" or "bow" regions, of the boat-shaped receptor binding motif on the RBD, appears to determine an interaction energy with ACE2 receptor higher than that observed with single mutations occurring at the "hull" region or with other multiple mutants. In addition, our pipeline allowed searching for ACE2 structurally related proteins, i.e., THOP1 and NLN, which deserve to be investigated for their possible involvement in interactions with the SARS-CoV-2 spike protein, in those tissues showing a low expression of ACE2, or as a novel receptor for future spike variants. A freely available web-tool for the in silico calculation of the interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein-protein interface, starting from the sequences of the investigated spike and/or ACE2 variants, was made available for the scientific community at: https://www.mitoairm.it/covid19affinities.

Conclusion

In the context of the PPPM/3PM, the employment of the described pipeline through the provided webservice, together with the ongoing SARS-CoV-2 genomic sequencing, would help to predict the transmissibility of new variants sequenced from future patients, depending on SARS-CoV-2 genomic sequencing activities and on the specific amino acid replacement and/or on its location on the SARS-CoV-2 spike RBD, to put in play all the possible counteractions for preventing the most deleterious scenarios of new outbreaks, taking into consideration that a greater transmissibility has not to be necessarily related to a more severe manifestation of the disease.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-021-00267-w.

Enhanced Molecular Dynamics Method to Efficiently Increase the Discrimination Capability of Computational Protein-Protein Docking

Protein–protein docking typically consists of the generation of putative binding conformations, which are subsequently ranked by fast heuristic scoring functions. The simplicity of these functions allows for computational efficiency but has severe repercussions on their discrimination capabilities. In this work, we show the effectiveness of suitable descriptors calculated along short scaled molecular dynamics runs in recognizing the nearest-native bound conformation among a set of putative structures generated by the HADDOCK tool for eight protein–protein systems.

Repeated MDA5 Gene Loss in Birds: An Evolutionary Perspective

Two key cytosolic receptors belonging to the retinoic acid-inducible gene I (RIG-I)-like receptor (RLR) family sense the viral RNA-derived danger signals: RIG-I and melanoma differentiation-associated protein 5 (MDA5). Their activation establishes an antiviral state by downstream signaling that ultimately activates interferon-stimulated genes (ISGs). While in rare cases RIG-I gene loss has been detected in mammalian and avian species, most notably in the chicken, MDA5 pseudogenization has only been detected once in mammals. We have screened over a hundred publicly available avian genome sequences and describe an independent disruption of MDA5 in two unrelated avian lineages, the storks (Ciconiiformes) and the rallids (Gruiformes). The results of our RELAX analysis confirmed the absence of negative selection in the MDA5 pseudogene. In contrast to our prediction, we have shown, using multiple dN/dS-based approaches, that the MDA5 loss does not appear to have resulted in any compensatory evolution in the RIG-I gene, which may partially share its ligand-binding specificity. Together, our results indicate that the MDA5 pseudogenization may have important functional effects on immune responsiveness in these two avian clades.

Quantum Speedup for Protein Structure Prediction

Protein structure prediction (PSP) predicts the native conformation for a given protein sequence. Classically, the problem has been shown to belong to the NP-complete complexity class. Its applications range from physics, through bioinformatics to medicine and quantum biology. It is possible however to speed it up with quantum computational methods, as we show in this paper. Here we develop a fast quantum algorithm for PSP in three-dimensional hydrophobic-hydrophilic model on body-centered cubic lattice with quadratic speedup over its classical counterparts. Given a protein sequence of n amino acids, our algorithm reduces the temporal and spatial complexities to, respectively, [Formula: see text] and O(n² logn) . With respect to oracle-related quantum algorithms for the NP-complete problems, we identify our algorithm as optimal. To justify the feasibility of the proposed algorithm we successfully solve the problem on IBM quantum simulator involving 21 and 25 qubits. We confirm the experimentally obtained high probability of success in finding the desired conformation by calculating the theoretical probability estimations.

The stability of R-spine defines RAF inhibitor resistance: A comprehensive analysis of oncogenic BRAF mutants with in-frame insertion of αC-β4 loop

Although targeting BRAF mutants with RAF inhibitors has achieved promising outcomes in cancer therapy, drug resistance remains a remarkable challenge, and underlying molecular mechanisms are not fully understood. Here, we characterized a previously unknown group of oncogenic BRAF mutants with in-frame insertions (LLR^ins506 or VLR^ins506) of αC-β4 loop. Using structure modeling and molecular dynamics simulation, we found that these insertions formed a large hydrophobic network that stabilizes R-spine and thus triggers the catalytic activity of BRAF. Furthermore, these insertions disrupted BRAF dimer interface and impaired dimerization. Unlike BRAF(V600E), these BRAF mutants with low dimer affinity were strongly resistant to all RAF inhibitors in clinic or clinical trials, which arises from their stabilized R-spines. As predicted by molecular docking, the stabilized R-spines in other BRAF mutants also conferred drug resistance. Together, our data indicated that the stability of R-spine but not dimer affinity determines the RAF inhibitor resistance of oncogenic BRAF mutants.

5-Methylcytosine Substantially Enhances the Thermal Stability of DNA Minidumbbells

Minidumbbell (MDB) is a recently identified non-B DNA structure that has been proposed to associate with genetic instabilities. It also serves as a functional structural motif in DNA nanotechnology. DNA molecular switches constructed using MDBs show instant and complete structural conversions with easy manipulations. The availability of stable MDBs can broaden their applications. In this work, we found that substitutions of cytosine with 5-methylcytosine could lead to a significant enhancement in the thermal stabilities of MDBs. Consecutive methylations of cytosine in MDBs brought about cumulative stabilization with a drastic increase in the melting temperature by 23 °C. NMR solution structures of two MDBs containing 5-methylcytosine residues have been successfully determined and revealed that the enhanced stabilities resulted primarily from favorable hydrophobic contacts, more stable base pairs and enhanced base-base stackings involving the methyl group of 5-methylcytosine.

Effects of Adenine Methylation on the Structure and Thermodynamic Stability of a DNA Minidumbbell

DNA methylation is a prevalent regulatory modification in prokaryotes and eukaryotes. N1-methyladenine (m1A) and N6-methyladenine (m6A) have been found to be capable of altering DNA structures via disturbing Watson-Crick base pairing. However, little has been known about their influences on non-B DNA structures, which are associated with genetic instabilities. In this work, we investigated the effects of m1A and m6A on both the structure and thermodynamic stability of a newly reported DNA minidumbbell formed by two TTTA tetranucleotide repeats. As revealed by the results of nuclear magnetic resonance spectroscopic studies, both m1A and m6A favored the formation of a T·m1A and T·m6A Hoogsteen base pair, respectively. More intriguingly, the m1A and m6A modifications brought about stabilization and destabilization effects on the DNA minidumbbell, respectively. This work provides new biophysical insights into the effects of adenine methylation on the structure and thermodynamic stability of DNA.

Polyphenols Weaken Pea Protein Gel by Formation of Large Aggregates with Diminished Noncovalent Interactions

Polyphenols are well-known native cross-linkers and gel strengthening agents for many animal proteins. However, their role in modifying plant protein gels remains unclear. In this study, multiple techniques were applied to unravel the influence of green tea polyphenols (GTP) on pea protein gels and the underlying mechanisms. We found that the elasticity and viscosity of pea protein gels decreased with increased GTP. The protein backbone became less rigid when GTP was present based on shortened T_1ρ^H in relaxation solid-state NMR measurements. Electron microscopy and small-angle X-ray scattering showed that gels weakened by GTP possessed disrupted networks with the presence of large protein aggregates. Solvent extraction and molecular dynamic simulation revealed a reduction in hydrophobic interactions and hydrogen bonds among proteins in gels containing GTP. The current findings may be applicable to other plant proteins for greater control of gel structures in the presence of polyphenols, expanding their utilization in food and biomedical applications.

Targeting mitochondrial metabolite transporters in Penicillium expansum for reducing patulin production

There is an increasing need of alternative treatments to control fungal infection and consequent mycotoxin accumulation in harvested fruits and vegetables. Indeed, only few biological targets of antifungal agents have been characterized and can be used for limiting fungal spread from decayed fruits/vegetables to surrounding healthy ones during storage. On this concern, a promising target of new antifungal treatments may be represented by mitochondrial proteins due to some species-specific functions played by mitochondria in fungal morphogenesis, drug resistance and virulence. One of the most studied mycotoxins is patulin produced by several species of Penicillium and Aspergillus genera. Patulin is toxic to many biological systems including bacteria, higher plants and animalia. Although precise biochemical mechanisms of patulin toxicity in humans are not completely clarified, its high presence in fresh and processed apple fruits and other apple-based products makes necessary developing a strategy for limiting its presence/accumulation. Patulin biosynthetic pathway consists of an enzymatic cascade, whose first step is represented by the synthesis of 6-methylsalicylic acid, obtained from the condensation of one acetyl-CoA molecule with three malonyl-CoA molecules. The most abundant acetyl-CoA precursor is represented by citrate produced by mitochondria. In the present investigation we report about the possibility to control patulin production through the inhibition of mitochondrial/peroxisome transporters involved in the export of acetyl-CoA precursors from mitochondria and/or peroxisomes, with specific reference to the predicted P. expansum mitochondrial Ctp1p, DTC, Sfc1p, Oac1p and peroxisomal PXN carriers.

Targeting Penicillium expansum GMC Oxidoreductase with High Affinity Small Molecules for Reducing Patulin Production

Simple Summary

With the urgent necessity of potential treatments for limiting mycotoxin production and postharvest fungal rots, we propose a combined in silico/in vitro/in vivo strategy for the rapid and effective identification of bioactive small molecules, chosen among a chemical library hosting approved drugs and phytochemicals, to be used after harvest. The molecular target of our analysis was the GMC oxidoreductase from Penicillium expansum involved in the biosynthesis of patulin, a mycotoxin that can contaminate many foods, especially fruits and fruit-based products. The employed in silico/in vitro/in vivo assays described in our study proved the effectiveness of our strategy and in particular of two small molecules, 6-hydroxycoumarin (structurally related to umbelliferon, an already characterized patulin synthase inhibitor) and meticrane (an already approved drug) in reducing patulin accumulation. Our findings highly recommend the mentioned ligands to be subjected to further analysis for being used in the next future in place of other more toxic compounds, in postharvest treatments based on dipping or drenching methods.

Abstract

Flavine adenine dinucleotide (FAD) dependent glucose methanol choline oxidoreductase (GMC oxidoreductase) is the terminal key enzyme of the patulin biosynthetic pathway. GMC oxidoreductase catalyzes the oxidative ring closure of (E)-ascladiol to patulin. Currently, no protein involved in the patulin biosynthesis in Penicillium expansum has been experimentally characterized or solved by X-ray diffraction. Consequently, nothing is known about P. expansum GMC oxidoreductase substrate-binding site and mode of action. In the present investigation, a 3D comparative model for P. expansum GMC oxidoreductase has been described. Furthermore, a multistep computational approach was used to identify P. expansum GMC oxidoreductase residues involved in the FAD binding and in substrate recognition. Notably, the obtained 3D comparative model of P. expansum GMC oxidoreductase was used for performing a virtual screening of a chemical/drug library, which allowed to predict new GMC oxidoreductase high affinity ligands to be tested in in vitro/in vivo assays. In vitro assays performed in presence of 6-hydroxycoumarin and meticrane, among the highly affinity predicted binders, confirmed a dose-dependent inhibition (17–81%) of patulin production by 6-hydroxycoumarin (10 µM–1 mM concentration range), whereas the approved drug meticrane inhibited patulin production by 43% already at 10 µM. Furthermore, 6-hydroxycoumarin and meticrane caused a 60 and 41% reduction of patulin production, respectively, in vivo on apples at 100 µg/wound.

Minidumbbell structures formed by ATTCT pentanucleotide repeats in spinocerebellar ataxia type 10

Spinocerebellar ataxia type 10 (SCA10) is a progressive genetic disorder caused by ATTCT pentanucleotide repeat expansions in intron 9 of the ATXN10 gene. ATTCT repeats have been reported to form unwound secondary structures which are likely linked to large-scale repeat expansions. In this study, we performed high-resolution nuclear magnetic resonance spectroscopic investigations on DNA sequences containing two to five ATTCT repeats. Strikingly, we found the first two repeats of all these sequences well folded into highly compact minidumbbell (MDB) structures. The 3D solution structure of the sequence containing two ATTCT repeats was successfully determined, revealing the MDB comprises a regular TTCTA and a quasi TTCT/A pentaloops with extensive stabilizing loop-loop interactions. We further carried out in vitro primer extension assays to examine if the MDB formed in the primer could escape from the proofreading function of DNA polymerase. Results showed that when the MDB was formed at 5-bp or farther away from the priming site, it was able to escape from the proofreading by Klenow fragment of DNA polymerase I and thus retained in the primer. The intriguing structural findings bring about new insights into the origin of genetic instability in SCA10.

Positive selection and convergent evolution shape molecular phenotypic traits of innate immunity receptors in tits (Paridae)

Despite widespread variability and redundancy abounding animal immunity, little is currently known about the rate of evolutionary convergence (functionally analogous traits not inherited from a common ancestor) in host molecular adaptations to parasite selective pressures. Toll-like receptors (TLRs) provide the molecular interface allowing hosts to recognize pathogenic structures and trigger early danger signals initiating an immune response. Using a novel combination of bioinformatic approaches, here we explore genetic variation in ligand-binding regions of bacteria-sensing TLR4 and TLR5 in 29 species belonging to the tit family of passerine birds (Aves: Paridae). Three out of the four consensual positively selected sites in TLR4 and six out of 14 positively selected positions in TLR5 were located on the receptor surface near the functionally important sites, and based on the phylogenetic pattern evolved in a convergent (parallel) manner. This type of evolution was also seen at one N-glycosylation site and two positively selected phosphorylation sites, providing the first evidence of convergence in post-translational modifications in evolutionary immunology. Finally, the overall mismatch between phylogeny and the clustering of surface charge distribution demonstrates that convergence is common in overall TLR4 and TLR5 molecular phenotypes involved in ligand binding. Our analysis did not reveal any broad ecological traits explaining the convergence observed in electrostatic potentials, suggesting that information on microbial symbionts may be needed to explain TLR evolution. Adopting state-of-the-art predictive structural bionformatics, we have outlined a new broadly applicable methodological approach to estimate the functional significance of positively selected variation and test for the adaptive molecular convergence in protein-coding polymorphisms.

High-Resolution Structures of DNA Minidumbbells Comprising Type II Tetraloops with a Purine Minor Groove Residue

Minidumbbell (MDB) is a newly discovered DNA structure formed by native sequences, which serves as a possible structural intermediate causing repeat expansion mutations in the genome and also a functional structural motif in constructing DNA-based molecular switches. Until now, all the reported MDBs containing two adjacent type II tetraloops were formed by pyrimidine-rich sequences 5'-YYYR YYYR-3' (Y and R represent pyrimidine and purine, respectively), wherein the second and sixth residues folded into the minor groove and interacted with each other. In this study, we have conducted a high-resolution nuclear magnetic resonance (NMR) spectroscopic investigation on alternative MDB-forming sequences and discovered that an MDB could also be formed stably with a purine in the minor groove, which has never been observed in any previously reported DNA type II tetraloops. Our refined NMR solution structures of the two MDBs formed by 5'-CTTG CATG-3' and 5'-CTTG CGTG-3' reveal that the sixth purine residue was driven into the minor groove via base-base stacking with the second thymine residue and adenine stacked better than guanine. The results of our present research work expand the sequence criteria for the formation of MDBs and shed light to explore the significance of MDBs.

PathMolD-AB: Spatiotemporal pathways of protein folding using parallel molecular dynamics with a coarse-grained model

Solving the protein folding problem (PFP) is one of the grand challenges still open in computational biophysics. Globular proteins are believed to evolve from initial configurations through folding pathways connecting several thermodynamically accessible states in a free energy landscape until reaching its minimum, inhabited by the stable native structures. Despite its huge computational burden, molecular dynamics (MD) is the leading approach in the PFP studies by preserving the Newtonian temporal evolution in the canonical ensemble. Non-trivial improvements are provided by highly parallel implementations of MD in cost-effective GPUs, concomitant to multiscale descriptions of proteins by coarse-grained minimalist models. In this vein, we present the PathMolD-AB framework, a comprehensive software package for massively parallel MD simulations using the canonical ensemble, structural analysis, and visualization of the folding pathways using the minimalist AB-model. It has, also, a tool to compare the results with proteins re-scaled from the PDB. We simulate and analyze, as case studies, the folding of four proteins: 13FIBO, 2GB1, 1PLC and 5ANZ, with 13, 55, 99 and 223 amino acids, respectively. The datasets generated from simulations correspond to the MD evolution of 3500 folding pathways, encompassing 35×106 states, which contains the spatial amino acid positions, the protein free energies and radii of gyration at each time step. Results indicate that the speedup of our approach grows logarithmically with the protein length and, therefore, it is suited for most of the proteins in the PDB. The predicted structures simulated by PathMolD-AB were similar to the re-scaled biological structures, indicating that it is promising for the study of the PFP study.

How RNA-Binding Proteins Interact with RNA: Molecules and Mechanisms

RNA-binding proteins (RBPs) comprise a large class of over 2,000 proteins that interact with transcripts in all manner of RNA-driven processes. The structures and mechanisms that RBPs use to bind and regulate RNA are incredibly diverse. In this review, we take a look at the components of protein-RNA interaction, from the molecular level to multi-component interaction. We first summarize what is known about protein-RNA molecular interactions based on analyses of solved structures. We additionally describe software currently available for predicting protein-RNA interaction and other resources useful for the study of RBPs. We then review the structure and function of seventeen known RNA-binding domains and analyze the hydrogen bonds adopted by protein-RNA structures on a domain-by-domain basis. We conclude with a summary of the higher-level mechanisms that regulate protein-RNA interactions.

A benchmark of optimally folded protein structures using integer programming and the 3D-HP-SC model

The Protein Structure Prediction (PSP) problem comprises, among other issues, forecasting the three-dimensional native structure of proteins using only their primary structure information. Most computational studies in this area use synthetic data instead of real biological data. However, the closer to the real-world, the more the impact of results and their applicability. This work presents 17 real protein sequences extracted from the Protein Data Bank for a benchmark to the PSP problem using the tri-dimensional Hydrophobic-Polar with Side-Chains model (3D-HP-SC). The native structure of these proteins was found by maximizing the number of hydrophobic contacts between the side-chains of amino acids. The problem was treated as an optimization problem and solved by means of an Integer Programming approach. Although the method optimally solves the problem, the processing time has an exponential trend. Therefore, due to computational limitations, the method is a proof-of-concept and it is not applicable to large sequences. For unknown sequences, an upper bound of the number of hydrophobic contacts (using this model) can be found, due to a linear relationship with the number of hydrophobic residues. The comparison between the predicted and the biological structures showed that the highest similarity between them was found with distance thresholds around 5.2-8.2 Å. Both the dataset and the programs developed will be freely available to foster further research in the area.

Ground state instabilities of protein shells are eliminated by buckling

We propose a hybrid discrete-continuum model to study the ground state of protein shells. The model allows for shape transformation of the shell and buckling transitions as well as the competition between states with different symmetries that characterize discrete particle models with radial pair potentials. Our main results are as follows. For large Föppl-von Kármán (FvK) numbers the shells have stable isometric ground states. As the FvK number is reduced, shells undergo a buckling transition resembling that of thin-shell elasticity theory. When the width of the pair potential is reduced below a critical value, then buckling coincides with the onset of structural instability triggered by over-stretched pair potentials. Chiral shells are found to be more prone to structural instability than achiral shells. It is argued that the well-width appropriate for protein shells lies below the structural instability threshold. This means that the self-assembly of protein shells with a well-defined, stable structure is possible only if the bending energy of the shell is sufficiently low so that the FvK number of the assembled shell is above the buckling threshold.

Scrutiny of electrostatic-driven conformational ordering of polypeptide chains in DMSO: a study with a model oligopeptide

The molecular mechanism of DMSO-induced stabilisation of β-sheets is attributed to the combination of polar electrostatic interactions among side chains, and backbone desolvation through bulky side chains which promotes backbone hydrogen bonding.

Probing the role of electrostatics of polypeptide main-chain in protein folding by perturbing N-terminal residue stereochemistry: DFT study with oligoalanine models

Energetics of folding (ΔH_E→F, in kcal mol⁻¹) from the extended (E) structure to the folded (F) structure for Ia and Ib critically depend on the geometrical relationship between the backbone peptide units of the polypeptide structure.

Contribution to the prediction of the fold code: application to immunoglobulin and flavodoxin cases

Background

Folding nucleus of globular proteins formation starts by the mutual interaction of a group of hydrophobic amino acids whose close contacts allow subsequent formation and stability of the 3D structure. These early steps can be predicted by simulation of the folding process through a Monte Carlo (MC) coarse grain model in a discrete space. We previously defined MIRs (Most Interacting Residues), as the set of residues presenting a large number of non-covalent neighbour interactions during such simulation. MIRs are good candidates to define the minimal number of residues giving rise to a given fold instead of another one, although their proportion is rather high, typically [15-20]% of the sequences. Having in mind experiments with two sequences of very high levels of sequence identity (up to 90%) but different folds, we combined the MIR method, which takes sequence as single input, with the “fuzzy oil drop” (FOD) model that requires a 3D structure, in order to estimate the residues coding for the fold. FOD assumes that a globular protein follows an idealised 3D Gaussian distribution of hydrophobicity density, with the maximum in the centre and minima at the surface of the “drop”. If the actual local density of hydrophobicity around a given amino acid is as high as the ideal one, then this amino acid is assigned to the core of the globular protein, and it is assumed to follow the FOD model. Therefore one obtains a distribution of the amino acids of a protein according to their agreement or rejection with the FOD model.

Results

We compared and combined MIR and FOD methods to define the minimal nucleus, or keystone, of two populated folds: immunoglobulin-like (Ig) and flavodoxins (Flav). The combination of these two approaches defines some positions both predicted as a MIR and assigned as accordant with the FOD model. It is shown here that for these two folds, the intersection of the predicted sets of residues significantly differs from random selection. It reduces the number of selected residues by each individual method and allows a reasonable agreement with experimentally determined key residues coding for the particular fold. In addition, the intersection of the two methods significantly increases the specificity of the prediction, providing a robust set of residues that constitute the folding nucleus.