A novel virus, severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2), causing coronavirus disease 2019 (COVID-19) worldwide appeared in 2019. Currently, we do not have a medicament that treats the disease. One of the rea-sons for the absence of treatment is related to the scarcity of detailed scientific knowledge of the members of the Coro-naviridae family, including the Middle East Respiratory Syndrome Coronavirus (MERS-CoV). Structural studies of the MERS-CoV proteins in the current literature are extremely limited. We present here detailed characterization of the struc-tural properties of MERS-CoV macro domain in aqueous solution at the atomic level with dynamics. For this study, we conducted extensive replica exchange molecular dynamics simulations linked to a generative neural networks and we use the resulting trajectories for structural analysis. We perform structural clustering based on the radius of gyration and end-to-end distance of MERS-CoV macro domain in aqueous solution with dynamics at the atomic level. We also report and analyze the residue-level intrinsic disorder features, flexibility and secondary structure. Furthermore, we study the pro-pensities of this macro domain for protein-protein interactions and for the RNA and DNA binding. Results are in agree-ment with available nuclear magnetic resonance spectroscopy findings and present more detailed insights into the struc-tural properties of MERS CoV macro domain. Overall, this work further shows that neural networks can be used as an exploratory tool for the studies of CoV family molecular conformational space at the nano level.

When two or more amino acid mutations occur in protein systems, they can interact in a non-additive fashion termed epistasis. One way to quantify epistasis between mutation pairs in protein systems is by using free energy differences: ϵ = 𝚫𝚫G1,2 - (𝚫𝚫G1 + 𝚫𝚫G2) where 𝚫𝚫G refers to the change in the Gibbs free energy, subscripts 1 and 2 refer to single mutations in arbitrary order and 1,2 refers to the double mutant. In this study, we explore possible biophysical mechanisms that drive pairwise epistasis in both protein-protein binding affinity and protein folding stability. Using the largest available datasets containing experimental protein structures and free energy data, we derived statistical models for both binding and folding epistasis (ϵ) with similar explanatory power (R2) of 0.299 and 0.258, respectively. These models contain terms and interactions that are consistent with intuition. For example, increasing the Cartesian separation between mutation sites leads to a decrease in observed epistasis for both folding and binding. Our results provide insight into factors that contribute to pairwise epistasis in protein systems and their importance in explaining epistasis. However, the low explanatory power indicates that more study is needed to fully understand this phenomenon.

CASP (Critical Assessment of Structure prediction) conducts community experiments to determine the state of the art in computing protein structure from amino acid sequence. The process relies on the experimental community providing information about not yet public or about to be solved structures, for use as targets. For some targets, the experimental structure is not solved in time for use in CASP. Calculated structure accuracy improved dramatically in this round, implying that models should now be much more useful for resolving many sorts of experimental difficulty. To test this, selected models for seven unsolved targets were provided to the experimental groups. These models were from the AlphaFold2 group, who overall submitted the most accurate predictions in CASP14. Four targets were solved with the aid of the models, and, additionally, the structure of an already solved target was improved. An a-posteriori analysis showed that in some cases models from other groups would also be effective. This paper provides accounts of the successful application of models to structure determination, including molecular replacement for X-ray crystallography, backbone tracing and sequence positioning in a Cryo-EM structure, and correction of local features. The results suggest that in future there will be greatly increased synergy between computational and experimental approaches to structure determination.

The Protein Data Bank (PDB) file format remains a popular format used and supported by many software to represent coordinates of macromolecular structures. It however suffers from drawbacks such as error-prone manual editing. Because of that, various software toolkits have been developed to facilitate its editing and manipulation, but, to date, there is no online tool available for this purpose. Here we present PDB-Tools Web, a flexible online service for manipulating PDB files. It offers a rich and user-friendly graphical user interface that allows users to mix-and-match more than 40 individual tools from the pdb-tools suite. Those can be combined in a few clicks to perform complex pipelines, which can be saved and uploaded. The resulting processed PDB files can be visualized online and downloaded. The web server is freely available at https://wenmr.science.uu.nl/pdbtools.

We present the results for CAPRI Round 54, the 5th joint CASP-CAPRI protein assembly prediction challenge. The Round offered 37 targets, including 14 homo-dimers, 3 homo-trimers, 13 hetero-dimers including 3 antibody-antigen complexes, and 7 large assemblies. On average ~70 CASP and CAPRI predictor groups, including more than 20 automatics servers, submitted models for each target. A total of 21941 models submitted by these groups and by 15 CAPRI scorer groups were evaluated using the CAPRI model quality measures and the DockQ score consolidating these measures. The prediction performance was quantified by a weighted score based on the number of models of acceptable quality or higher submitted by each group among their 5 best models. Results show substantial progress achieved across a significant fraction of the 60+ participating groups. High-quality models were produced for about 40% for the targets compared to 8% two years earlier, a remarkable improvement resulting from the wide use of the AlphaFold2 and AlphaFold-Multimer software. Creative use was made of the deep learning inference engines affording the sampling of a much larger number of models and enriching the multiple sequence alignments with sequences from various sources. Wide use was also made of the AlphaFold confidence metrics to rank models, permitting top performing groups to exceed the results of the public AlphaFold-Multimer version used as a yard stick. This notwithstanding, performance remained poor for complexes with antibodies and nanobodies, where evolutionary relationships between the binding partners are lacking, and for complexes featuring conformational flexibility, clearly indicating that the prediction of protein complexes remains a challenging problem.

Functional regulation via conformational dynamics is well known in structured proteins, but less well characterized in intrinsically disordered proteins and their complexes. Using NMR spectroscopy we have identified a dynamic regulatory mechanism in the human insulin-like growth factor (IGF) system involving the central, intrinsically disordered linker domain of human IGF-binding protein-2 ( hIGFBP2). The bioavailability of IGFs is regulated by the proteolysis of IGF-binding proteins. In the case of hIGFBP2, the linker domain (L- hIGFBP2) retains its intrinsic disorder upon binding IGF-1 but its dynamics are significantly altered, both in the IGF binding region and distantly located protease cleavage sites. The increase in flexibility of the linker domain upon IGF-1 binding may explain the IGF-dependent modulation of proteolysis of IGFBP2 in this domain. As IGF homeostasis is important for cell growth and function, and its dysregulation is a key contributor to several cancers, our findings open up new avenues for the design of IGFBP analogs inhibiting IGF-dependent tumors.

Several plastic degrading enzymes have been described in the literature, most notably PETases that are capable of hydrolyzing polyethylene terephthalate (PET) plastic. One of them, the PETase from Ideonella sakaiensis, a bacterium isolated from environmental samples within a PET bottle recycling site, was the subject of extensive studies. To test how widespread PETase functionality is in other bacterial communities, we used a cascade of BLAST searches in the JGI metagenomic datasets and showed that PETases can also be found in other metagenomic environmental samples from both human affected and relatively pristine sites. To confirm their classification as PETases, we verified that the newly identified proteins have the PETase sequence signatures common to all PETases and that phylogenetic analyses group them with the experimentally characterized PETases. Additionally, docking analysis was performed in order to further confirm the functional assignment of the putative environmental PETases.

The novel coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) still has serious negative effects on health, social life, and economics. Recently, vaccines from various companies have been urgently approved to control SARS-CoV-2 infections. However, any specific antiviral drug has not been confirmed so far for regular treatment. An important target is the main protease (Mpro), which plays a major role in replication of the virus. In this study, Gaussian and residue network models are employed to reveal two distinct potential allosteric sites on Mpro that can be evaluated as drug targets besides the active site. Then, FDA-approved drugs are docked to three distinct sites with flexible docking using AutoDock Vina to identify potential drug candidates. 14 best molecule hits for the active site of Mpro are determined. 6 of these also exhibit high docking scores for the potential allosteric regions. Full-atom molecular dynamics simulations with MM-GBSA method indicate that compounds docked to active and potential allosteric sites form stable interactions with high binding free energy (∆Gbind) values. ∆Gbind values reach -52.06 kcal/mol for the active site, -51.08 kcal/mol for the potential allosteric site 1, and -42.93 kcal/mol for the potential allosteric site 2. Energy decomposition calculations per residue elucidate key binding residues stabilizing the ligands that can further serve to design pharmacophores. This systematic and efficient computational analysis successfully determines ivermectine, diosmin and selinexor currently subjected to clinical trials, and further proposes bromocriptine, elbasvir as Mpro inhibitor candidates to be evaluated against SARS-CoV-2 infection

Degradation of solid polyethylene terephthalate (PET) by leaf branch compost cutinase (LCC) produces various PET-derived degradation intermediates (DIs), in addition to terephthalic acid (TPA), which is the recyclable terminal product of all PET degradation. Although DIs can also be converted into TPA, in solution, by LCC, the TPA that is obtained through enzymatic degradation of PET, in practice, is always contaminated by DIs. Here, we demonstrate that the primary reason for non-degradation of DIs into TPA in solution is the efficient binding of LCC onto the surface of solid PET. Although such binding enhances the degradation of solid PET, it depletes the surrounding solution of enzyme that could otherwise have converted DIs into TPA. To retain a sub-population of enzyme in solution that would mainly degrade DIs, we introduced mutations to reduce the hydrophobicity of areas surrounding LCC’s active site, with the express intention of reducing LCC’s binding to solid PET. Despite the consequent reduction in invasion and degradation of solid PET, overall levels of production of TPA were ~3.6-fold higher, due to the partitioning of enzyme between solid PET and the surrounding solution, and the consequent heightened production of TPA from DIs. Further, synergy between such mutated LCC (F125L/F243I LCC) and wild-type LCC resulted in even higher yields, and TPA of nearly ~100% purity.

The structure of a protein plays a pivotal role in determining its function. Often, the protein surface’s shape and curvature dictate its nature of interaction with other proteins and biomolecules. However, marked by corrugations and roughness, a protein’s surface representation poses significant challenges for its curvature-based characterization. In the present study, we employ unsupervised machine learning to segment the protein surface into patches. To measure the surface curvature of a patch, we present an algebraic sphere fitting method that is fast, accurate, and robust. Moreover, we use local curvatures to show the existence of “shape complementarity” in protein-protein, antigen-antibody, and protein-ligand interfaces. We believe that the current approach could help understand the relationship between protein structure and its biological function and can be used to find binding partners of a given protein.

Order and disorder govern protein functions, but there is a great diversity in disorder, from regions that are – and stay – fully disordered to conditional order. This diversity is still difficult to decipher even though it is encoded in the amino acid sequences. Here, we developed an analytic Python package, named pyHCA, to estimate the foldability of a protein segment from the only information of its amino acid sequence and based on a measure of its density in regular secondary structures associated with hydrophobic clusters, as defined by the Hydrophobic Cluster Analysis (HCA) approach. The tool was designed by optimizing the separation between foldable segments from databases of disorder (DisProt) and order (SCOPe (soluble domains) and OPM (transmembrane domains)). It allows to specify the ratio between order, embodied by regular secondary structures (either participating in the hydrophobic core of well-folded 3D structures or conditionally formed in intrinsically disordered regions) and disorder. We illustrated the relevance of pyHCA with several examples and applied it to the sequences of the proteomes of 21 species ranging from prokaryotes and archaea to unicellular and multicellular eukaryotes, for which structure models are provided in the AlphaFold2 databases. Cases of low-confidence scores related to disorder were distinguished from those of sequences that we identified as foldable but are still excluded from accurate modeling by AlphaFold2 due to a lack of sequence homologs or to compositional biases. Overall, our approach is complementary to AlphaFold2, providing guides to map structural innovations through evolutionary processes, at proteome and gene scales.

Isoflavonoid is one of the groups of flavonoids that play pivotal roles in the survival of land plants. Chalcone synthase (CHS), the first enzyme of the isoflavonoid biosynthetic pathway, catalyzes the formation of a common isoflavonoid precursor. We have previously reported that an isozyme of soybean CHS (termed GmCHS1) is a key component of the isoflavonoid metabolon, a protein complex to enhance efficiency of isoflavonoid production. Here, we determined the crystal structure of GmCHS1 as a first step of understanding the metabolon structure, as well as to better understand the catalytic mechanism of GmCHS1.

Petroleum based plastics are durable and accumulate in all ecological niches. Knowledge on enzymatic degradation is sparse. Today, less than 50 verified plastics-active enzymes are known. First examples of enzymes acting on the polymers polyethylene terephthalate (PET) and polyurethane (PUR) have been reported together with a detailed biochemical and structural description. Further, very few polyamide (PA) oligomer active enzymes are known. In this paper, the current known enzymes acting on the synthetic polymers PET and PUR are briefly summarized, their published activity data were collected and integrated into a comprehensive open access database. The Plastics-Active Enzymes Database (PAZy) represents an inventory of known and experimentally verified plastics-active enzymes. Almost 3000 homologues of PET-active enzymes were identified by profile hidden Markov models. Over 2000 homologues of PUR-active enzymes were identified by BLAST. Based on multiple sequence alignments, conservation analysis identified the most conserved amino acids, and sequence motifs for PET- and PUR-active enzymes were derived.

The recognition of Cannabis as a source of new compounds suitable for medical use has attracted strong interest from the scientific community in its research, and substantial progress has accumulated regarding cannabinoids’ activity; however, a thorough description of their molecular mechanisms of action remains a task to complete. Highlighting their complex pharmacology, the list of cannabinoids’ interactors has vastly expanded beyond the canonical cannabinoid receptors. Among those, we have focused our study on the glycine receptor (GlyR), an ion channel involved in the modulation of nervous system responses, including, to our interest, sensitivity to peripheral pain. Here, we report the use of computational methods to investigate possible binding modes between the GlyR and Δ 9-tetrahydrocannabinol (THC). After obtaining a first pose for the THC binding from a biased molecular docking simulation and subsequently evaluating it by molecular dynamic simulations, we found a dynamic system with an identifiable representative binding mode characterized by the specific interaction with two transmembrane residues (Phe293 and Ser296). Complementarily, we assessed the role of membrane cholesterol in this interaction and positively established its relevance for THC binding to GlyR. Lastly, the use of restrained molecular dynamics simulations allowed us to refine the description of the binding mode and of the cholesterol effect. Altogether, our findings contribute to the current knowledge about the GlyR-THC mode of binding and propose a new starting point for future research on how cannabinoids in general, and THC in particular, modulate pain perception in view of its possible clinical applications.

SARS-CoV-2 is neutralized by proteins that block receptor-binding sites on spikes that project from the viral envelope. In particular, substantial research investment has advanced monoclonal antibody therapies to the clinic where there are signs of partial efficacy in reducing viral burden and hospitalization. An alternative is to use the host entry receptor, ACE2, as a soluble decoy that broadly blocks SARS-associated coronaviruses with limited potential for viral escape. Here, we summarize efforts to engineer higher affinity variants of soluble ACE2 that rival the potency of affinity-matured antibodies. Strategies have also been used to increase the valency of ACE2 decoys for avid spike interactions and to improve pharmacokinetics via IgG fusions. Finally, the intrinsic catalytic activity of ACE2 for the turnover of the vasoconstrictor angiotensin II may directly address COVID-19 symptoms and protect against lung and cardiovascular injury, conferring dual mechanisms of action unachievable by monoclonal antibodies. Soluble ACE2 derivatives therefore have the potential to be next generation therapeutics for addressing the immediate needs of the current pandemic and possible future outbreaks.

The core metabolic reactions of life drive electrons through a class of redox protein enzymes, the oxidoreductases. The energetics of electron flow is determined by the redox potentials of organic and inorganic cofactors as tuned by the protein environment. Understanding how protein structure affects oxidation-reduction energetics is crucial for studying metabolism, creating bioelectronic systems, and tracing the history of biological energy utilization on Earth. We constructed ProtReDox ([https://protein-redox-potential.web.app](https://protein-redox-potential.web.app)), a manually curated database of experimentally determined redox potentials. With over 500 measurements, we can begin to identify how proteins modulate oxidation-reduction energetics across the tree of life. By mapping redox potentials onto networks of oxidoreductase fold evolution, we can infer the evolution of electron transfer energetics over deep-time. ProtReDox is designed to include user-contributed submissions with the intention of making it a valuable resource for researchers in this field.

Dissimilatory sulfite reductase is an ancient enzyme that has linked the global sulfur and carbon biogeochemical cycles since at least 3.47 Gya. While much has been learned about the phylogenetic distribution and diversity of DsrAB across environmental gradients, far less is known about the structural changes that occurred to maintain DsrAB function as the enzyme accompanied diversification of sulfate/sulfite reducing organisms (SRO) into new environments. Analyses of available crystal structures of DsrAB from Archaeoglobus fulgidus and Desulfovibrio vulgaris, representing early and late evolving lineages, respectively, show that certain features of DsrAB are structurally conserved, including active siro-heme binding motifs. Whether such structural features are conserved among DsrAB recovered from varied environments, including hot spring environments that host representatives of the earliest evolving SRO lineage (e.g., MV2-Eury), is not known. To begin to overcome these gaps in our understanding of the evolution of DsrAB, structural models from MV2.Eury were generated and evolutionary sequence co-variance analyses were conducted on a curated DsrAB database. Phylogenetically diverse DsrAB harbor many conserved functional residues including those that ligate active siro-heme(s). However, evolutionary co-variance analysis of monomeric DsrAB subunits revealed several False Positive Evolutionary Couplings (FPEC) that correspond to residues that have co-evolved despite being too spatially distant in the monomeric structure to allow for direct contact. One set of FPECs corresponds to residues that form a structural path between the two active siro-heme moieties across the interface between heterodimers, suggesting the potential for allostery or electron transfer within the enzyme complex. Other FPECs correspond to structural loops and gaps that may have been selected to stabilize enzyme function in different environments. These structural bioinformatics results suggest that DsrAB has maintained allosteric communication pathways between subunits as SRO diversified into new environments. The observations outlined here provide a framework for future biochemical and structural analyses of DsrAB to examine potential allosteric control of this enzyme.

The goal of CASP experiments is to monitor the progress in the protein structure prediction field. During the 14th CASP edition we aimed to test our capabilities of predicting structures of protein complexes. Our protocol for modeling protein assemblies included both template-based modeling and free docking. Structural templates were identified using sensitive sequence-based searches. If sequence-based searches failed, we performed structure-based template searches using selected CASP server models. In the absence of reliable templates we applied free docking starting from monomers generated by CASP servers. We evaluated and ranked models of protein complexes using an improved version of protein structure quality assessment method, VoroMQA, taking into account both interaction interface and global structure scores. If reliable templates could be identified, generally accurate models of protein assemblies were generated with the exception of an antibody-antigen interaction. The success of free docking mainly depended on the accuracy of initial subunit models and on the scoring of docking solutions. To put our overall results in perspective, we analyzed our performance in the context of other CASP groups. Although the subunits in our assembly models often were not of the top quality, these models had, overall, the best predicted interfaces according to several protein-protein interface accuracy measures. Since we did not use co-evolution-based prediction of inter-chain contacts, we attribute our relative success in predicting interfaces primarily to the emphasis on the interaction interface when modeling and scoring.