The activation or inactivation of B-cell lymphoma-2 (Bcl-2) antagonist/killer (Bak) is critical for controlling mitochondrial outer membrane permeabilization-dependent apoptosis. Its pro-apoptotic activity is controlled by intermolecular interactions with the Bcl-2 homology 3 (BH3) domain, which is accommodated in the hydrophobic pocket of Bak. Bcl-2-interacting protein 5 (Bnip5) is a noncanonical BH3 domain-containing protein that interacts with Bak. Bnip5 is characterized by its controversial effects on the regulation of the pro-apoptotic activity of Bak. In the present study, we determined the crystal structure of Bak bound to Bnip5 BH3. The intermolecular association appeared to be typical at first glance, but we found that it is maintained by tight hydrophobic interactions together with hydrogen/ionic bonds, which accounts for their high binding affinity with a dissociation constant of 775 nM. Structural analysis of the complex showed that Bnip5 interacts with Bak in a manner similar to that of the Bak-activating pro-apoptotic factor peroxisomal testis-enriched protein 1, particularly in the destabilization of the intramolecular electrostatic network of Bak. Our structure is considered to reflect the initial point of drastic and consecutive conformational and stoichiometric changes in Bak induced by Bnip5 BH3, which helps in explaining the effects of Bnip5 in regulating Bak-mediated apoptosis.

This article reports and analyzes the results of protein complex model accuracy estimation by our methods (DeepUMQA3 and GraphGPSM) in the 15th Critical Assessment of techniques for protein Structure Prediction (CASP15). The new deep learning-based multimeric complex model accuracy estimation methods are proposed based on the ensemble of three-level features coupling with deep residual/graph neural networks. For the input multimeric complex model, we describe it from three levels: overall complex features, intra-monomer features, and inter-monomer features. We designed an overall ultrafast shape recognition (USR) to characterize the relationship between local residues and the overall complex topology, and an inter-monomer USR to characterize the relationship between the residues of one monomer and the topology of other monomers. DeepUMQA3 (Group name: GuijunLab-RocketX) ranked first in the interface residue accuracy estimation of CASP15. The Pearson correlation between the interface residue Local Distance Difference Test (lDDT) predicted by DeepUMQA3 and the real lDDT is 0.570, the only method that exceeds 0.5. Among the top 5 methods, DeepUMQA3 achieved the highest Pearson correlation of lDDT on 25 out of 39 targets. GraphGPSM (Group name: GuijunLab-PAthreader) has TM-score Pearson correlations greater than 0.9 on 14 targets, showing a good ability to estimate the overall fold accuracy. The DeepUMQA3 server is available at http://zhanglab-bioinf.com/DeepUMQA/ and the GraphGPSM server is available at http://zhanglab-bioinf.com/GraphGPSM/.

AlphaFold2 has revolutionized structure prediction by achieving high accuracy comparable to experimentally determined structures. However, there is still room for improvement, especially for challenging cases like multimers. A key to the success of AlphaFold is its ability to assess and rank its own predictions. Our basic idea for the Wallner group in CASP15 was to exploit this excellent scoring function in AlphaFold by massive sampling. To achieve this goal, we conducted AlphaFold runs using six different settings, using templates, without templates, and with an increased number of recycles for both multimer v1 and v2 weights. In all instances, we enabled dropout layers during inference, allowing for sampling of uncertainty and enhancing the diversity of the generated models. In total, 274 289 models were generated for the 38 targets in CASP15, with a median of 4810 models per target. Of these 38 targets, 10 were high quality, 11 were medium quality, 11 were acceptable, and only 6 were incorrect. The improvement over the baseline method, NBIS-AF2-multimer, is substantial, with the mean DockQ increasing from 0.43 to 0.56, with several targets showing a DockQ score increase of +0.6 units. Remarkable, considering Wallner and NBIS-AF2-multimer were using identical input data. The success can be attributed to the diversified sampling using dropout with different settings and, in particular, the use of multimer v1, which is much more susceptible to sampling compared with v2. The method is available here: http://wallnerlab.org/AFsample/.

We report the results of the "UM-TBM" and "Zheng" groups in CASP15 for protein monomer and complex structure prediction. These prediction sets were obtained using the D-I-TASSER and DMFold-Multimer algorithms, respectively. For monomer structure prediction, D-I-TASSER introduced four new features during CASP15: (i) a multiple sequence alignment (MSA) generation protocol that combines multi-source MSA searching and a structural modeling-based MSA ranker; (ii) attention-network based spatial restraints; (iii) a multi-domain module containing domain partition and arrangement for domain-level templates and spatial restraints; (iv) an optimized I-TASSER-based folding simulation system for full-length model creation guided by a combination of deep learning restraints, threading alignments, and knowledge-based potentials. For 47 free modeling targets in CASP15, the final models predicted by D-I-TASSER showed average TM-score 19% higher than the standard AlphaFold2 program. We thus showed that traditional Monte Carlo-based folding simulations, when appropriately coupled with deep learning algorithms, can generate models with improved accuracy over end-to-end deep learning methods alone. For protein complex structure prediction, DMFold-Multimer generated models by integrating a new MSA generation algorithm (DeepMSA2) with the end-to-end modeling module from AlphaFold2-Multimer. For the 38 complex targets, DMFold-Multimer generated models with an average TM-score of 0.83 and Interface Contact Score of 0.60, both significantly higher than those of competing complex prediction tools. Our analyses on complexes highlighted the critical role played by MSA generating, ranking, and pairing in protein complex structure prediction. We also discuss future room for improvement in the areas of viral protein modeling and complex model ranking.

Cell invasion is an important process in cancer progression and recurrence. Invasion and implantation of cancer cells from their original place to other tissues, by disabling vital organs, challenges the treatment of cancer patients. Given the importance of the matter, many molecular treatments have been developed to inhibit cancer cell invasion. Because of their low production cost and ease of production, peptides are valuable therapeutic molecules for inhibiting cancer cell invasion. In recent years, advances in the field of computational biology have facilitated the design of anti-cancer peptides. In our investigation, using computational biology approaches such as evolutionary analysis, residue scanning, protein-peptide interaction analysis, molecular dynamics, and free energy analysis, our team designed a peptide library with about 100 000 candidates based on A6 (acetyl-KPSSPPEE-amino) sequence which is an anti-invasion peptide. During computational studies, two of the designed peptides that give the highest scores and showed the greatest sequence similarity to A6 were entered into the experimental analysis workflow for further analysis. In experimental analysis steps, the anti-metastatic potency and other therapeutic effects of designed peptides were evaluated using MTT assay, RT-qPCR, zymography analysis, and invasion assay. Our study disclosed that the IK1 (acetyl-RPSFPPEE-amino) peptide, like A6, has great potency to inhibit the invasion of cancer cells.

Cyanobacteriochromes (CBCRs) are linear tetrapyrrole bilin-binding photoreceptors of cyanobacteria that exhibit high spectral diversity, gaining attention in optogenetics and bioimaging applications. Several engineering studies on CBCRs were attempted, especially for designing near-infrared (NIR) fluorescent proteins with longer fluorescence wavelengths. However, despite continuous efforts, a key component regulating fluorescence emission property in CBCRs is still poorly understood. As a model system, we focused on red/green CBCR Slr1393g3, from the unicellular cyanobacterium Synechocystis sp. PCC 6803 to engineer Pr to get far-red light-emitting property. Energy profiling and pairwise structural comparison of Slr1393g3 variants effectively reveal the mutations that are critical to the fluorescence changes. H497 seems to play a key role in stabilizing the chromophore environment, especially the α3 helix, while H495, T499, and Q502 are potential key residues determining fluorescence emission peak wavelength. We also found that mutations of α2 and α4 helical regions are closely related to the chromophore binding stability and likely affect fluorescence properties. Taken together, our computational analysis suggests that the fluorescence of Slr1393g3 is mainly controlled by the stabilization of the chromophore binding pocket. The predicted key residues potentially regulating the fluorescence emission property of a red/green CBCR will be advantageous for designing improved NIR fluorescent protein when combined with in vitro molecular evolution approaches.

Methyl parathion hydrolase (MPH) is an enzyme of the metallo-β-lactamase superfamily, which hydrolyses a wide range of organophosphates (OPs). Recently, MPH has attracted attention as a promising enzymatic bioremediator. The crystal structure of MPH enzyme shows a dimeric form, with each subunit containing a binuclear metal ion center. MPH also demonstrates metal ion-dependent selectivity patterns. The origins of these patterns remain unclear but are linked to open questions about the more general role of metal ions in functional evolution and divergence within enzyme superfamilies. We aimed to investigate and compare the binding of different OP pesticides to MPH with cobalt(II) metal ions. In this study, MPH was modeled from Ochrobactrum sp. with different OP pesticides bound, including methyl paraoxon and dichlorvos and profenofos. The docked structures for each substrate optimized by DFT calculation were selected and subjected to atomistic molecular dynamics simulations for 500 ns. It was found that alpha metal ions did not coordinate with all the pesticides. Rather, the pesticides coordinated with less buried beta metal ions. It was also observed that the coordination of beta metal ions was perturbed to accommodate the pesticides. The binding free energy calculations and structure-based pharmacophore model revealed that all the three substrates could bind well at the active site. However, profenofos exhibit a stronger binding affinity to MPH in comparison to the other two substrates. Therefore, our findings provide molecular insight on the binding of different OP pesticides which could help us design the enzyme for OP pesticides degradation.

We describe the modeling method for RNA tertiary structures employed by team AIchemy_RNA2 in the 15th Critical Assessment of Structure Prediction (CASP15). The method consists of the following steps. Firstly, secondary structure information was derived from various manually-verified sources. With this information, the full length RNA was fragmented into structural modules. The structures of each module were predicted and then assembled into the full structure. To reduce the searching conformational space, an RNA structure was organized into an optimal base folding tree. And to further improve the sampling efficiency, the energy surface was smoothed at high temperatures during the Monte Carlo sampling to make it easier to move across the energy barrier. The statistical potential energy function BRiQ was employed during Monte Carlo energy optimization.

Proteins are played key roles in different functionalities in our daily life. All functional roles of a protein are a bit enhanced in interaction compared to individuals. Identification of essential proteins of an organism is a time consume and costly task during observation in the wet lab. The results of observation in wet lab always ensure high reliability and accuracy in the biological ground. Essential protein prediction using computational approaches is an alternative choice in research. It proves its significance rapidly in day-to-day life as well as reduces the experimental cost of wet lab effectively. Existing computational methods were implemented using Protein interaction networks (PPIN), Sequence, Gene Expression Dataset (GED), Gene Ontology (GO), Orthologous groups, and Subcellular localized datasets. Machine learning has diverse categories of features that enable to model and predict essential macromolecules of understudied organisms. A novel methodology MEM-FET (membership feature) is predicted based on features, that is, edge clustering coefficient, Average clustering coefficient, subcellular localization, and Gene Ontology within a compartment of common neighbors. The accuracy (ACC) values of the predicted true positive (TP) essential proteins are 0.79, 0.74, 0.78, and 0.71 for YHQ, YMIPS, YDIP, and YMBD datasets. An enriched set of essential proteins are also predicted using the MEM-FET algorithm. Ensemble ML also validated the proposed model with an accuracy of 60%. It has been predicted that MEM-FET algorithms outperform other existing algorithms with an ACC value of 80% for the yeast dataset.

Ribonucleic acid (RNA) molecules serve as master regulators of cells by encoding their biological function in the ribonucleotide sequence, particularly their ability to interact with other molecules. To understand how RNA molecules perform their biological tasks and to design new sequences with specific functions, it is of great benefit to be able to computationally predict how RNA folds and interacts in the cellular environment. Our workflow for computational modeling of the 3D structures of RNA and its interactions with other molecules uses a set of methods developed in our laboratory, including MeSSPredRNA for predicting canonical and non-canonical base pairs, PARNASSUS for detecting remote homology based on comparisons of sequences and secondary structures, ModeRNA for comparative modeling, the SimRNA family of programs for modeling RNA 3D structure and its complexes with other molecules, and QRNAS for model refinement. In this study, we present the results of testing this workflow in predicting RNA 3D structures in the CASP15 experiment. The overall high score of the computational models predicted by our group demonstrates the robustness of our workflow and its individual components in terms of predicting RNA 3D structures of acceptable quality that are close to the target structures. However, the variance in prediction quality is still quite high, and the results are still too far from the level of protein 3D structure predictions. This exercise led us to consider several improvements, especially to better predict and enforce stacking interactions and non-canonical base pairs.

Protein intrinsic disorder (ID) is a complex and context-dependent phenomenon that covers a continuum between fully disordered states and folded states with long dynamic regions. The lack of a ground truth that fits all ID flavors and the potential for order-to-disorder transitions depending on specific conditions makes ID prediction challenging. The CAID2 challenge aimed to evaluate the performance of different prediction methods across different benchmarks, leveraging the annotation provided by the DisProt database, which stores the coordinates of ID regions when there is experimental evidence in the literature. The CAID2 challenge demonstrated varying performance of different prediction methods across different benchmarks, highlighting the need for continued development of more versatile and efficient prediction software. Depending on the application, researchers may need to balance performance with execution time when selecting a predictor. Methods based on AlphaFold2 seem to be good ID predictors but they are better at detecting absence of order rather than ID regions as defined in DisProt. The CAID2 predictors can be freely used through the CAID Prediction Portal, and CAID has been integrated into OpenEBench, which will become the official platform for running future CAID challenges.

Fatty acids are used in fundamental cellular processes, such as membrane biogenesis, energy generation, post-translational modification of proteins, and so forth. These processes require the activation of fatty acids by adenosine triphosphate (ATP), followed by condensation with coenzyme-A (CoA), catalyzed by the omnipresent enzyme called Fatty acyl-CoA ligases (FACLs). However, Fatty acyl-AMP ligases (FAALs), the structural homologs of FACLs, operate in an unprecedented CoA-independent manner. FAALs transfer fatty acids to the acyl carrier protein (ACP) domain of polyketide synthases (PKS) and non-ribosomal peptide synthetases (NRPS) for the biosynthesis of various antibiotics, lipopeptides, virulent complex lipids, and so forth in bacteria. Recent structural and biochemical insights from our group provide a detailed understanding of the mode of CoA rejection and ACP acceptance by FAALs. In this review, we have discussed advances in the mechanistic, evolutionary, and functional understanding of FAALs and FAAL-like domains across life forms. Here, we are proposing a "Five-tier" mechanistic model to explain the specificity of FAALs. We further demonstrate how FAAL-like domains have been repurposed into a new family of proteins in eukaryotes with a novel function in lipid metabolism.

As CASP15 participants, in the new category of 3D RNA structure prediction, we applied expert modeling with the support of our proprietary system RNAComposer. Although RNAComposer is primarily known as an automated web server, its features allow it to be used interactively, for example, for homology-based modeling or assembling models from user-provided structural elements. In the paper, we present various scenarios of applying the system to predict the 3D RNA structures that we employed. Their combination with expert input, comparative analysis of models, and routines to select representative resultant structures form a ready-for-reuse workflow. With selected examples, we demonstrate its application for the in silico modeling of natural and synthetic RNA molecules targeted in CASP15.

In CASP15, we used an integrated hierarchical and hybrid approach to predict RNA structures. The approach involves three steps. First, with the use of physics-based methods, Vfold2D-MC and VfoldMCPX, we predict the 2D structures from the sequence. Second, we employ template-based methods, Vfold3D and VfoldLA, to build 3D scaffolds for the predicted 2D structures. Third, using the 3D scaffolds as initial structures and the predicted 2D structures as constraints, we predict the 3D structure from coarse-grained molecular dynamics simulations, IsRNA and RNAJP. Our approach was evaluated on 12 RNA targets in CASP15 and ranked second among all the 34 participating teams. The result demonstrated the reliability of our method in predicting RNA 2D structures with high accuracy and RNA 3D structures with moderate accuracy. Further improvements in RNA structure prediction for the next round of CASP may come from the incorporation of the physics-based method with machine learning techniques.

We introduce a deep learning-based ligand pose scoring model called zPoseScore for predicting protein-ligand complexes in the 15th Critical Assessment of Protein Structure Prediction (CASP15). Our contributions are threefold: first, we generate six training and evaluation data sets by employing advanced data augmentation and sampling methods. Second, we redesign the "zFormer" module, inspired by AlphaFold2's Evoformer, to efficiently describe protein-ligand interactions. This module enables the extraction of protein-ligand paired features that lead to accurate predictions. Finally, we develop the zPoseScore framework with zFormer for scoring and ranking ligand poses, allowing for atomic-level protein-ligand feature encoding and fusion to output refined ligand poses and ligand per-atom deviations. Our results demonstrate excellent performance on various testing data sets, achieving Pearson's correlation R  = 0.783 and 0.659 for ranking docking decoys generated based on experimental and predicted protein structures of CASF-2016 protein-ligand complexes. Additionally, we obtain an averaged local distance difference test (lDDT pli = 0.558) of AIchemy LIG2 in CASP15 for de novo protein-ligand complex structure predictions. Detailed analysis shows that accurate ligand binding site prediction and side-chain orientation are crucial for achieving better prediction performance. Our proposed model is one of the most accurate protein-ligand pose prediction models and could serve as a valuable tool in small molecule drug discovery.

Recently, we presented a framework for understanding protein structure based on the idea that simple constructs of holding hands or touching of objects can be used to rationalize the common characteristics of globular proteins. We developed a consistent approach for understanding the formation of the two key common building blocks of helices and sheets as well as the compatible assembly of secondary structures into the tertiary structure through the notion of poking pairwise interactions. Here we benchmark our predictions with a detailed analysis of structural data of over 4000 proteins from the Protein Data Bank. We also present the results of detailed computer simulations of a simplified model demonstrating a pre-sculpted free energy landscape, determined by geometry and symmetry, comprising numerous minima corresponding to putative native state structures. We explore the consequences of our model. Our results suggest that symmetry and geometry are a powerful guide to capture the simplicity underlying protein complexity.

Before the controversial approval of humanized monoclonal antibody lecanemab, which binds to the soluble amyloid-β protofibrils, all the treatments available earlier, for Alzheimer's disease (AD) were symptomatic. The researchers are still struggling to find a breakthrough in AD therapeutic medicine, which is partially attributable to lack in understanding of the structural information associated with the intrinsically disordered proteins and amyloids. One of the major challenges in this area of research is to understand the structural diversity of intrinsically disordered proteins under in vitro conditions. Therefore, in this review, we have summarized the in vitro applications of biophysical methods, which are aimed to shed some light on the heterogeneity, pathogenicity, structures and mechanisms of the intrinsically disordered protein aggregates associated with proteinopathies including AD. This review will also rationalize some of the strategies in modulating disease-relevant pathogenic protein entities by small molecules using structural biology approaches and biophysical characterization. We have also highlighted tools and techniques to simulate the in vivo conditions for native and cytotoxic tau/amyloids assemblies, urge new chemical approaches to replicate tau/amyloids assemblies similar to those in vivo conditions, in addition to designing novel potential drugs.

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), 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.

Diversity in the biochemical workhorses of the cell-that is, proteins-is achieved by the innumerable permutations offered primarily by the 20 canonical L-amino acids prevalent in all biological systems. Yet, proteins are known to additionally undergo unusual modifications for specialized functions. Of the various post-translational modifications known to occur in proteins, the recently identified non-disulfide cross-links are unique, residue-specific covalent modifications that confer additional structural stability and unique functional characteristics to these biomolecules. We review an exclusive class of amino acid cross-links encompassing aromatic and sulfur-containing side chains, which not only confer superior biochemical characteristics to the protein but also possess additional spectroscopic features that can be exploited as novel chromophores. Studies of their in vivo reaction mechanism have facilitated their specialized in vitro applications in hydrogels and protein anchoring in monolayer chips. Furthering the discovery of unique canonical cross-links through new chemical, structural, and bioinformatics tools will catalyze the development of protein-specific hyperstable nanostructures, superfoods, and biotherapeutics.

Proteins often function as part of permanent or transient multimeric complexes, and understanding function of these assemblies requires knowledge of their three-dimensional structures. While the ability of AlphaFold to predict structures of individual proteins with unprecedented accuracy has revolutionized structural biology, modeling structures of protein assemblies remains challenging. To address this challenge, we developed a protocol for predicting structures of protein complexes involving model sampling followed by scoring focused on the subunit-subunit interaction interface. In this protocol, we diversified AlphaFold models by varying construction and pairing of multiple sequence alignments as well as increasing the number of recycles. In cases when AlphaFold failed to assemble a full protein complex or produced unreliable results, additional diverse models were constructed by docking of monomers or subcomplexes. All the models were then scored using a newly developed method, VoroIF-jury, which relies only on structural information. Notably, VoroIF-jury is independent of AlphaFold self-assessment scores and therefore can be used to rank models originating from different structure prediction methods. We tested our protocol in CASP15 and obtained top results, significantly outperforming the standard AlphaFold-Multimer pipeline. Analysis of our results showed that the accuracy of our assembly models was capped mainly by structure sampling rather than model scoring. This observation suggests that better sampling, especially for the antibody-antigen complexes, may lead to further improvement. Our protocol is expected to be useful for modeling and/or scoring protein assemblies.

We present a model, based on symmetry and geometry, for proteins. Using elementary ideas from mathematics and physics, we derive the geometries of discrete helices and sheets. We postulate a compatible solvent-mediated emergent pairwise attraction that assembles these building blocks, while respecting their individual symmetries. Instead of seeking to mimic the complexity of proteins, we look for a simple abstraction of reality that yet captures the essence of proteins. We employ analytic calculations and detailed Monte Carlo simulations to explore some consequences of our theory. The predictions of our approach are in accord with experimental data. Our framework provides a rationalization for understanding the common characteristics of proteins. Our results show that the free energy landscape of a globular protein is pre-sculpted at the backbone level, sequences and functionalities evolve in the fixed backdrop of the folds determined by geometry and symmetry, and that protein structures are unique in being simultaneously characterized by stability, diversity, and sensitivity.

We present the monomer and multimer structure prediction results of our methods in CASP15. We first designed an elaborate pipeline that leverages complementary sequence databases and advanced database searching algorithms to generate high-quality multiple sequence alignments (MSAs). Top MSAs were then selected for the subsequent step of structure prediction. We utilized trRosettaX2 and AlphaFold2 for monomer structure prediction (group name Yang-Server), and AlphaFold-Multimer for multimer structure prediction (group name Yang-Multimer). Yang-Server and Yang-Multimer are ranked at the top and the fourth, respectively, for monomer and multimer structure prediction. For 94 monomers, the average TM-score of the predicted structure models by Yang-Server is 0.876, compared to 0.798 by the default AlphaFold2 (i.e., the group NBIS-AF2-standard). For 42 multimers, the average DockQ score of the predicted structure models by Yang-Multimer is 0.464, compared to 0.389 by the default AlphaFold-Multimer (i.e., the group NBIS-AF2-multimer). Detailed analysis of the results shows that several factors contribute to the improvement, including improved MSAs, iterated modeling for large targets, interplay between monomer and multimer structure prediction for intertwined structures, etc. However, the structure predictions for orphan proteins and multimers remain challenging, and breakthroughs in this area are anticipated in the future.

Kainate receptors are a subtype of ionotropic glutamate receptors that form transmembrane channels upon binding glutamate. Here, we have investigated the mechanism of partial agonism in heteromeric GluK2/K5 receptors, where the GluK2 and GluK5 subunits have distinct agonist binding profiles. Using single-molecule Förster resonance energy transfer, we found that at the bi-lobed agonist-binding domain, the partial agonist AMPA-bound receptor occupied intermediate cleft closure conformational states at the GluK2 cleft, compared to the more open cleft conformations in apo form and more closed cleft conformations in the full agonist glutamate-bound form. In contrast, there is no significant difference in cleft closure states at the GluK5 agonist-binding domain between the partial agonist AMPA- and full agonist glutamate-bound states. Additionally, unlike the glutamate-bound state, the dimer interface at the agonist-binding domain is not decoupled in the AMPA-bound state. Our findings suggest that partial agonism observed with AMPA binding is mediated primarily due to differences in the GluK2 subunit, highlighting the distinct contributions of the subunits towards activation.