Background: Biomolecular interactions that modulate biological processes occur mainly in cavities throughout the surface of biomolecular structures. In the data science era, structural biology has benefited from the increasing availability of biostructural data due to advances in structural determination and computational methods. In this scenario, data-intensive cavity analysis demands efficient scripting routines built on easily manipulated data structures. To fulfill this need, we developed pyKVFinder, a Python package to detect and characterize cavities in biomolecular structures for data science and automated pipelines. Results: pyKVFinder efficiently detects cavities in biomolecular structures and computes their volume, area, depth and hydropathy, storing these cavity properties in NumPy arrays. Benefited from Python ecosystem interoperability and data structures, pyKVFinder can be integrated with third-party scientific packages and libraries for mathematical calculations, machine learning and 3D visualization in automated workflows. As proof of pyKVFinder’s capabilities, we successfully identified and compared ADRP substrate-binding site of SARS-CoV-2 and a set of homologous proteins with pyKVFinder, showing its integrability with data science packages such as matplotlib, NGL Viewer, SciPy and Jupyter notebook. Conclusions: We introduce an efficient, highly versatile and easily integrable software for detecting and characterizing biomolecular cavities in data science applications and automated protocols. pyKVFinder facilitates biostructural data analysis with scripting routines in the Python ecosystem and can be building blocks for data science and drug design applications.

The nucleocapsid (N) protein of the SARS-CoV-2 virus, the causal agent of COVID-19, is a multifunction phosphoprotein that plays critical roles in the virus life cycle, including transcription and packaging of the viral RNA. To play such diverse roles, the N protein has two globular RNA-binding modules, the N- (NTD) and C-terminal (CTD) domains, which are connected by an intrinsically disordered region. Despite the wealth of structural data available for the isolated NTD and CTD, how these domains are arranged in the full-length protein and how the oligomerization of N influences its RNA-binding activity remains largely unclear. Herein, using experimental data from electron microscopy and biochemical/biophysical techniques combined with molecular modeling and molecular dynamics simulations, we show that, in the absence of RNA, the N protein formed structurally dynamic dimers, with the NTD and CTD arranged in extended conformations. However, in the presence of RNA, the N protein assumed a more compact conformation where the NTD and CTD are packed together. We also provided an octameric model for the full-length N bound to RNA that is consistent with electron microscopy images of the N protein in the presence of RNA. Together, our results shed new light on the dynamics and higher-order oligomeric structure of this versatile protein.

Deep learning methods, trained on the increasing set of available protein 3D structures and sequences, have substantially impacted the protein modeling and design field. These advancements have facilitated the creation of novel proteins, or the optimization of existing ones designed for specific functions, such as binding a target protein. Despite the demonstrated potential of such approaches in designing general protein binders, their application in designing immunotherapeutics remains relatively underexplored. A relevant application is the design of T cell receptors (TCRs). Given the crucial role of T cells in mediating immune responses, redirecting these cells to tumor or infected target cells through the engineering of TCRs has shown promising results in treating diseases, especially cancer. However, the computational design of TCR interactions presents challenges for current physics-based methods, particularly due to the unique natural characteristics of these interfaces, such as low affinity and cross-reactivity. For this reason, in this study, we explored the potential of two structure-based deep learning protein design methods, ProteinMPNN and ESM-IF1, in designing fixed-backbone TCRs for binding target antigenic peptides presented by the MHC through different design scenarios. To evaluate TCR designs, we employed a comprehensive set of sequence- and structure-based metrics, highlighting the benefits of these methods in comparison to classical physics-based design methods and identifying deficiencies for improvement.

Coronavirus disease 2019 is a multi-systemic syndrome that caused a pandemic. Proteomic studies have shown changes in protein expression and interaction involved in signaling pathways related to SARS-CoV-2 infections. Protein phosphatases play a crucial role in regulating cell signaling. In this study, we assessed the potential involvement of protein phosphatases and their associated signaling pathways during SARS-CoV-2 infection by conducting a meta-analysis of proteome databases from COVID-19 patients. We identified both direct and indirect interactions between human protein phosphatases and viral proteins, as well as the expression levels and phosphorylation status of intermediate proteins. Our analyses revealed that PPP2CA and PTEN are key phosphatases involved in cell cycle and apoptosis regulation during SARS-CoV-2 infection. We also highlighted the direct involvement of PPP2CA in the cell division throughout its interaction with CDC20 protein (cell division cycle protein 20 homolog). This evidence strongly suggests that both proteins play critical roles during SARS-CoV-2 infection and represent potential targets for COVID-19 treatment.

Transition-state (TS) theory has provided the theoretical framework to explain the enormous rate accelerations of chemical reactions by enzymes. Given that proteins display large ensembles of conformations, unique TSs would pose a huge entropic bottleneck for enzyme catalysis. To shed light on this question, we studied the nature of the enzymatic TS for the phosphoryl-transfer step in adenylate kinase by quantum-mechanics/molecular-mechanics calculations. We find a structurally wide set of energetically equivalent configurations that lie along the reaction coordinate and hence a broad transition-state ensemble (TSE). A conformationally delocalized ensemble, including asymmetric TSs, is rooted in the macroscopic nature of the enzyme. The computational results are buttressed by enzyme kinetics experiments that confirm the decrease of the entropy of activation predicted from such wide TSE. TSEs as a key for efficient enzyme catalysis further boosts a unifying concept for protein folding and conformational transitions underlying protein function.