Computational optimization of angiotensin-converting enzyme 2 for SARS-CoV-2 Spike molecular recognition

Zepyros: a webserver to evaluate the shape complementarity of protein-protein interfaces

Motivation

Shape complementarity of molecular surfaces at the interfaces is a well-known characteristic of protein-protein binding regions, and it is critical in influencing the stability of the complex. Measuring such complementarity is of great importance for a number of theoretical and practical implications; however, only a limited number of tools are currently available to efficiently and rapidly assess it.

Results

Here, we introduce Zepyros (ZErnike Polynomials analYsis of pROtein Shapes), a webserver for fast measurement of the shape complementarity between two molecular interfaces of a given protein-protein complex using structural information. Zepyros is implemented as a publicly available tool with a user-friendly interface.

Availability and implementation

Our server can be found at the following link (all major browser supported): https://zepyros.bio-groups.com.

Compact Assessment of Molecular Surface Complementarities Enhances Neural Network-Aided Prediction of Key Binding Residues

Predicting interactions between proteins is fundamental for understanding the mechanisms underlying cellular processes, since protein-protein complexes are crucial in physiological conditions but also in many diseases, for example by seeding aggregates formation. Despite the many advancements made so far, the performance of docking protocols is deeply dependent on their capability to identify binding regions. From this, the importance of developing low-cost and computationally efficient methods in this field. We present an integrated novel protocol mainly based on compact modeling of protein surface patches via sets of orthogonal polynomials to identify regions of high shape/electrostatic complementarity. By incorporating both hydrophilic and hydrophobic contributions, we define new binding matrices, which serve as effective inputs for training a neural network. In this work, we propose a new Neural Network (NN)-based architecture, Core Interacting Residues Network (CIRNet), which achieves a performance in terms of Area Under the Receiver Operating Characteristic Curve (ROC AUC) of approximately 0.87 in identifying pairs of core interacting residues on a balanced data set. In a blind search for core interacting residues, CIRNet distinguishes them from random decoys with an ROC AUC of 0.72. We test this protocol to enhance docking algorithms by filtering the proposed poses, addressing one of the still open problems in computational biology. Notably, when applied to the top ten models from three widely used docking servers, CIRNet improves docking outcomes, significantly reducing the average RMSD between the selected poses and the native state. Compared to another state-of-the-art tool for rescaling docking poses, CIRNet more efficiently identified the worst poses generated by the three docking servers under consideration and achieved superior rescaling performance in two cases.

The impact of hypertension on clinical manifestations of Omicron variant BA.1 infection in adult patients

COVID-19 caused by Omicron BA.1 has resulted in a global humanitarian crisis. In this COVID-19 pandemic era, hypertension has been receiving increased attention. Omicron BA.1 infection combined with hypertension created a serious public health problem and complicated the treatment and prognosis of COVID-19. The aim of our study was to assess the implications of hypertension for the clinical manifestations of adult patients (APs) infected with Omicron BA.1. This single-center retrospective cohort study enrolled consecutive COVID-19 APs, who were admitted to Tianjin First Central Hospital from 01 August 2022 to 30 November 2022. All included APs were divided into two groups: hypertension and non-hypertension group. The APs' baseline demographic, laboratory, clinical, and radiological characteristics were collected and analyzed. Of 512 APs admitted with PCR proven COVID-19, 161 (31.45%) APs had comorbid hypertension. Hypertension APs have older age, higher body mass index, lower Ct-values of the viral target genes at admission, and longer hospital stay than non-hypertension APs. Furthermore, hypertension aggravates the clinical classification, impairs liver, kidney, and myocardium function, and abnormalizes the coagulation system in Omicron BA.1- infected APs. Moreover, hypertension elevates inflammation levels and lung lesion involvement while weakened virus-specific IgM level in APs with Omicron BA.1 infection. Hypertension APs tend to have worse clinical conditions at baseline than those non-hypertension APs. This study indicates that hypertension is a contributor to the poor clinical manifestations of Omicron BA.1-infected APs and supports that steps to control blood pressure should be a vital consideration for reducing the burden of Omicron BA.1 infection in hypertension individuals.

IMPORTANCE

This study provided inclusive insight regarding the relationship between hypertension and Omicron BA.1 infection and supported that hypertension was an adverse factor for COVID-19 APs. In conclusion, this study showed that hypertension was considered to be associated with severe conditions, and a contributor to poor clinical manifestations. Proper medical management of hypertension patients is an imperative step in mitigating the severity of Omicron BA.1 variant infection.

Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

Design of protein-binding peptides with controlled binding affinity: the case of SARS-CoV-2 receptor binding domain and angiotensin-converting enzyme 2 derived peptides

The development of methods able to modulate the binding affinity between proteins and peptides is of paramount biotechnological interest in view of a vast range of applications that imply designed polypeptides capable to impair or favour Protein-Protein Interactions. Here, we applied a peptide design algorithm based on shape complementarity optimization and electrostatic compatibility and provided the first experimental in vitro proof of the efficacy of the design algorithm. Focusing on the interaction between the SARS-CoV-2 Spike Receptor-Binding Domain (RBD) and the human angiotensin-converting enzyme 2 (ACE2) receptor, we extracted a 23-residues long peptide that structurally mimics the major interacting portion of the ACE2 receptor and designed in silico five mutants of such a peptide with a modulated affinity. Remarkably, experimental KD measurements, conducted using biolayer interferometry, matched the in silico predictions. Moreover, we investigated the molecular determinants that govern the variation in binding affinity through molecular dynamics simulation, by identifying the mechanisms driving the different values of binding affinity at a single residue level. Finally, the peptide sequence with the highest affinity, in comparison with the wild type peptide, was expressed as a fusion protein with human H ferritin (HFt) 24-mer. Solution measurements performed on the latter constructs confirmed that peptides still exhibited the expected trend, thereby enhancing their efficacy in RBD binding. Altogether, these results indicate the high potentiality of this general method in developing potent high-affinity vectors for hindering/enhancing protein-protein associations.

Computational structural-based GPCR optimization for user-defined ligand: Implications for the development of biosensors

Organisms have developed effective mechanisms to sense the external environment. Human-designed biosensors exploit this natural optimization, where different biological machinery have been adapted to detect the presence of user-defined molecules. Specifically, the pheromone pathway in the model organism Saccharomyces cerevisiae represents a suitable candidate as a synthetic signaling system. Indeed, it expresses just one G-Protein Coupled Receptor (GPCR), Ste2, able to recognize pheromone and initiate the expression of pheromone-dependent genes. To date, the standard procedure to engineer this system relies on the substitution of the yeast GPCR with another one and on the modification of the yeast G-protein to bind the inserted receptor. Here, we propose an innovative computational procedure, based on geometrical and chemical optimization of protein binding pockets, to select the amino acid substitutions required to make the native yeast GPCR able to recognize a user-defined ligand. This procedure would allow the yeast to recognize a wide range of ligands, without a-priori knowledge about a GPCR recognizing them or the corresponding G protein. We used Monte Carlo simulations to design on Ste2 a binding pocket able to recognize epinephrine, selected as a test ligand. We validated Ste2 mutants via molecular docking and molecular dynamics. We verified that the amino acid substitutions we identified make Ste2 able to accommodate and remain firmly bound to epinephrine. Our results indicate that we sampled efficiently the huge space of possible mutants, proposing such a strategy as a promising starting point for the development of a new kind of S.cerevisiae-based biosensors.

Investigating the competition between ACE2 natural molecular interactors and SARS-CoV-2 candidate inhibitors

The SARS-CoV-2 pandemic still poses a threat to the global health as the virus continues spreading in most countries. Therefore, the identification of molecules capable of inhibiting the binding between the ACE2 receptor and the SARS-CoV-2 spike protein is of paramount importance. Recently, two DNA aptamers were designed with the aim to inhibit the interaction between the ACE2 receptor and the spike protein of SARS-CoV-2. Indeed, the two molecules interact with the ACE2 receptor in the region around the K353 residue, preventing its binding of the spike protein. If on the one hand this inhibition process hinders the entry of the virus into the host cell, it could lead to a series of side effects, both in physiological and pathological conditions, preventing the correct functioning of the ACE2 receptor. Here, we discuss through a computational study the possible effect of these two very promising DNA aptamers, investigating all possible interactions between ACE2 and its experimentally known molecular partners. Our in silico predictions show that some of the 10 known molecular partners of ACE2 could interact, physiologically or pathologically, in a region adjacent to the K353 residue. Thus, the curative action of the proposed DNA aptamers could recruit ACE2 from its biological functions.

Reversal of the unique Q493R mutation increases the affinity of Omicron S1-RBD for ACE2

The SARS-CoV-2 Omicron variant containing 15 mutations, including the unique Q493R, in the spike protein receptor binding domain (S1-RBD) is highly infectious. While comparison with previously reported mutations provide some insights, the mechanism underlying the increased infections and the impact of the reversal of the unique Q493R mutation seen in BA.4, BA.5, BA.2.75, BQ.1 and XBB lineages is not yet completely understood. Here, using structural modelling and molecular dynamics (MD) simulations, we show that the Omicron mutations increases the affinity of S1-RBD for ACE2, and a reversal of the unique Q493R mutation further increases the ACE2-S1-RBD affinity. Specifically, we performed all atom, explicit solvent MD simulations using a modelled structure of the Omicron S1-RBD-ACE2 and compared the trajectories with the WT complex revealing a substantial reduction in the Cα-atom fluctuation in the Omicron S1-RBD and increased hydrogen bond and other interactions. Residue level analysis revealed an alteration in the interaction between several residues including a switch in the interaction of ACE2 D38 from S1-RBD Y449 in the WT complex to the mutated R residue (Q493R) in Omicron complex. Importantly, simulations with Revertant (Omicron without the Q493R mutation) complex revealed further enhancement of the interaction between S1-RBD and ACE2. Thus, results presented here not only provide insights into the increased infectious potential of the Omicron variant but also a mechanistic basis for the reversal of the Q493R mutation seen in some Omicron lineages and will aid in understanding the impact of mutations in SARS-CoV-2 evolution.

Optimizing variant-specific therapeutic SARS-CoV-2 decoys using deep-learning-guided molecular dynamics simulations

Treatment of COVID-19 with a soluble version of ACE2 that binds to SARS-CoV-2 virions before they enter host cells is a promising approach, however it needs to be optimized and adapted to emerging viral variants. The computational workflow presented here consists of molecular dynamics simulations for spike RBD-hACE2 binding affinity assessments of multiple spike RBD/hACE2 variants and a novel convolutional neural network architecture working on pairs of voxelized force-fields for efficient search-space reduction. We identified hACE2-Fc K31W and multi-mutation variants as high-affinity candidates, which we validated in vitro with virus neutralization assays. We evaluated binding affinities of these ACE2 variants with the RBDs of Omicron BA.3, Omicron BA.4/BA.5, and Omicron BA.2.75 in silico. In addition, candidates produced in Nicotiana benthamiana, an expression organism for potential large-scale production, showed a 4.6-fold reduction in half-maximal inhibitory concentration (IC50) compared with the same variant produced in CHO cells and an almost six-fold IC50 reduction compared with wild-type hACE2-Fc.

Multi-task learning for predicting SARS-CoV-2 antibody escape

The coronavirus pandemic has revolutionized our world, with vaccination proving to be a key tool in fighting the disease. However, a major threat to this line of attack are variants that can evade the vaccine. Thus, a fundamental problem of growing importance is the identification of mutations of concern with high escape probability. In this paper we develop a computational framework that harnesses systematic mutation screens in the receptor binding domain of the viral Spike protein for escape prediction. The framework analyzes data on escape from multiple antibodies simultaneously, creating a latent representation of mutations that is shown to be effective in predicting escape and binding properties of the virus. We use this representation to validate the escape potential of current SARS-CoV-2 variants.

Shape Complementarity Optimization of Antibody–Antigen Interfaces: The Application to SARS-CoV-2 Spike Protein

Many factors influence biomolecule binding, and its assessment constitutes an elusive challenge in computational structural biology. In this aspect, the evaluation of shape complementarity at molecular interfaces is one of the main factors to be considered. We focus on the particular case of antibody–antigen complexes to quantify the complementarities occurring at molecular interfaces. We relied on a method we recently developed, which employs the 2D Zernike descriptors, to characterize the investigated regions with an ordered set of numbers summarizing the local shape properties. Collecting a structural dataset of antibody–antigen complexes, we applied this method and we statistically distinguished, in terms of shape complementarity, pairs of the interacting regions from the non-interacting ones. Thus, we set up a novel computational strategy based on in silico mutagenesis of antibody-binding site residues. We developed a Monte Carlo procedure to increase the shape complementarity between the antibody paratope and a given epitope on a target protein surface. We applied our protocol against several molecular targets in SARS-CoV-2 spike protein, known to be indispensable for viral cell invasion. We, therefore, optimized the shape of template antibodies for the interaction with such regions. As the last step of our procedure, we performed an independent molecular docking validation of the results of our Monte Carlo simulations.

In-Silico Evidence for a Two Receptor Based Strategy of SARS-CoV-2

We propose a computational investigation on the interaction mechanisms between SARS-CoV-2 spike protein and possible human cell receptors. In particular, we make use of our newly developed numerical method able to determine efficiently and effectively the relationship of complementarity between portions of protein surfaces. This innovative and general procedure, based on the representation of the molecular isoelectronic density surface in terms of 2D Zernike polynomials, allows the rapid and quantitative assessment of the geometrical shape complementarity between interacting proteins, which was unfeasible with previous methods. Our results indicate that SARS-CoV-2 uses a dual strategy: in addition to the known interaction with angiotensin-converting enzyme 2, the viral spike protein can also interact with sialic-acid receptors of the cells in the upper airways.