Surface energetics and protein-protein interactions: analysis and mechanistic implications

Molecular recognition of SARS-CoV-2 spike protein with three essential partners: exploring possible immune escape mechanisms of viral mutants

OBJECTIVE

The COVID-19 epidemic is raging around the world, with the emergence of viral mutant strains such as Delta and Omicron, posing severe challenges to people's health and quality of life. A full understanding life cycle of the virus in host cells helps to reveal inactivation mechanism of antibody and provide inspiration for the development of a new-generation vaccines.

METHODS

In this work, molecular recognitions and conformational changes of SARS-CoV-2 spike protein mutants (i.e., Delta, Mu, and Omicron) and three essential partners (i.e., membrane receptor hACE2, protease TMPRSS2, and antibody C121) both were compared and analyzed using molecular simulations.

RESULTS

Water basin and binding free energy calculations both show that the three mutants possess higher affinity for hACE2 than WT, exhibiting stronger virus transmission. The descending order of cleavage ability by TMPRSS2 is Mu, Delta, Omicron, and WT, which is related to the new S1/S2 cutting site induced by transposition effect. The inefficient utilization of TMPRSS2 by Omicron is consistent with its primary entry into cells via the endosomal pathway. In addition, RBD-directed antibody C121 showed obvious resistance to Omicron, which may have originated from high fluctuation of approaching angles, high flexibility of I472-F490 loop, and reduced binding ability.

CONCLUSIONS

According to the overall characteristics of the three mutants, high infectivity, high immune escape, and low virulence may be the future evolutionary selection of SARS-CoV-2. In a word, this work not only proposes the possible resistance mechanism of SARS-CoV-2 mutants, but also provides theoretical guidance for the subsequent drug design against COVID-19 based on S protein structure.

Spatial organization of hydrophobic and charged residues affects protein thermal stability and binding affinity

What are the molecular determinants of protein–protein binding affinity and whether they are similar to those regulating fold stability are two major questions of molecular biology, whose answers bring important implications both from a theoretical and applicative point of view. Here, we analyze chemical and physical features on a large dataset of protein–protein complexes with reliable experimental binding affinity data and compare them with a set of monomeric proteins for which melting temperature data was available. In particular, we probed the spatial organization of protein (1) intramolecular and intermolecular interaction energies among residues, (2) amino acidic composition, and (3) their hydropathy features. Analyzing the interaction energies, we found that strong Coulombic interactions are preferentially associated with a high protein thermal stability, while strong intermolecular van der Waals energies correlate with stronger protein–protein binding affinity. Statistical analysis of amino acids abundances, exposed to the molecular surface and/or in interaction with the molecular partner, confirmed that hydrophobic residues present on the protein surfaces are preferentially located in the binding regions, while charged residues behave oppositely. Leveraging on the important role of van der Waals interface interactions in binding affinity, we focused on the molecular surfaces in the binding regions and evaluated their shape complementarity, decomposing the molecular patches in the 2D Zernike basis. For the first time, we quantified the correlation between local shape complementarity and binding affinity via the Zernike formalism. In addition, considering the solvent interactions via the residue hydropathy, we found that the hydrophobicity of the binding regions dictates their shape complementary as much as the correlation between van der Waals energy and binding affinity. In turn, these relationships pave the way to the fast and accurate prediction and design of optimal binding regions as the 2D Zernike formalism allows a rapid and superposition-free comparison between possible binding surfaces.

Interactions of proteins with metal-based nanoparticles from a point of view of analytical chemistry - Challenges and opportunities

Interactions of proteins with nanomaterials draw attention of many research groups interested in fundamental phenomena. However, alongside with valuable information regarding physicochemical aspects of such processes and their mechanisms, they more and more often prove to be useful from a point of view of bioanalytics. Deliberate use of processes based on adsorption of proteins on nanoparticles (or vice versa) allows for a development of new analytical methods and improvement of the existing ones. It also leads to obtaining of nanoparticles of desired properties and functionalities, which can be used as elements of analytical tools for various applications. Due to interactions with nanoparticles, proteins can also gain new functionalities or lose their interfering potential, which from perspective of bioanalytics seems to be very inviting and attractive. In the framework of this article we will discuss the bioanalytical potential of interactions of proteins with a chosen group of nanoparticles, and implementation of so driven processes for biosensing. Moreover, we will show both positive and negative (opportunities and challenges) aspects resulting from the presence of proteins in media/samples containing metal-based nanoparticles or their precursors.

Modeling the Dynamics of Protein-Protein Interfaces, How and Why?

Protein-protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein-protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein-protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein-protein interface.

Protein Surface Printer for Exploring Protein Domains

The surface of proteins is vital in determining protein functions. Herein, a program, Protein Surface Printer (PSP), is built that performs multiple functions in quantifying protein surface domains. Two proteins, PETase and cytochrome P450, are used to validate that the program supports atomistic simulations with different combinations of programs and force fields. A case study is conducted on the structural analysis of the spike proteins of SARS-CoV-2 and SARS-CoV and the human cell receptor ACE2. Although the surface domains of both spike proteins are highly similar, their receptor-binding domains (RBDs) and the O-linked glycan domains are structurally different. The O-linked glycan domain of SARS-CoV-2 is highly positively charged, which may promote binding to negatively charged human cells.

Development of peptide aptamers as alternatives for antibody in the detection of amyloid-beta 42 aggregates

Alzheimer's disease (AD) causes cognitive impairment and serious social isolation. However, there are no effective treatments and even no established confirmatory diagnostic tools for the disease. Amyloid beta (Aβ) aggregation in the brain is the best-known pathognomonic mechanism of AD, so various methods for Aβ detection have been developed for the diagnosis of this disease. We synthesized two novel, ultra-sensitive peptide probes specialized in detecting Aβ aggregates, and examined their potential for future diagnostic application. The peptides are produced through phage high-throughput screening (HTS) and amplified through a serial process called biopanning, which is a repeating method of elution and amplification of probes. We picked phages specific for amyloid from two kinds of phage display. The synthesized peptides were confirmed to have excellent binding affinity to Aβ aggregates, by immunohistochemical staining and western blotting using the brains of 3X transgenic (Tg) AD mice at different stages (5-7, 12-17 months old) of AD severity. In the present study, it was confirmed that newly developed amyloid-binding peptides could be used as novel probes for the detection of Aβ aggregates, which can be used for clinical diagnosis of AD in the future.

Surface charge of Merkel cell polyomavirus small T antigen determines cell transformation through allosteric FBW7 WD40 domain targeting

Merkel cell polyomavirus (MCV) small T (sT) is the main oncoprotein in Merkel cell carcinoma (MCC) development. A unique domain of sT, LT stabilization domain (LSD), has been reported to bind and inactivate multiple SCF (Skp1-Cullin-F-box) E3 ligases. These interactions impede the turnover of MCV large T (LT) antigen and cellular oncoproteins such as c-Myc and cyclin E, thereby promoting viral replication and cell transformation. However, it is currently unclear how this LSD region contributes to multiple transforming activities of sT. Structural docking simulation of sT and F-box and WD repeat domain-containing 7 (FBW7) revealed a novel allosteric interaction between sT and FBW7 WD40 domain. This model is supported by experimental evidence confirming that charge engineering in the LSD alters sT-WD40 binding. Specifically, loss of net positive charge in the LSD prevents sT-FBW7 binding by abrogating the electrostatic interaction, thus impeding inhibition of FBW7 by sT. Furthermore, positively charged mutations in the LSD significantly restored the sT function and its ability to transform rodent fibroblast cells. We infer that the surface charge of sT is a major determinant for targeting E3 ligases, which leads to sT-induced cell transformation, an observation that could be used to develop targeted therapeutics for MCC.

A spatial similarity of stereochemical environments formed by amino acid residues defines a common epitope of two non-homologous proteins

It is critical for development of high-quality antibodies in research and diagnostics to predict accurately their cross-reactivities with "off-target" molecules, which potentially induce false results. Herein, we report a good example of such a cross-reactivity for an off-target due to a stereochemical environment of epitopes, which does not simply depend on amino acid sequences. We found that significant subpopulation of a polyclonal peptide antibody against Bcnt (Bucentaur) (anti-BCNT-C antibody) cross-reacted with a completely different protein, glutamine synthetase (GS), and identified four amino acids, GYFE, in its C-terminal region as the core amino acids for the cross-reaction. Consistent with this finding, the anti-BCNT-C antibody strongly recognized endogenously and exogenously expressed GS in tissues and cultured cells by Western blotting and immunohistochemistry. Furthermore, we elucidated that the cross-reaction is caused by a spatial similarity between the stereochemical environments formed by amino acid residues, including the GYFE of GS and the GYIE of Bcnt, rather than by their primary sequences. These results suggest it is critical to comprehensively analyze antibody interactions with target molecules including off-targets with special attention to the physicochemical environments of epitope-paratope interfaces to decrease the risk of false interpretations of results using antibodies in science and clinical applications.

Water follows polar and nonpolar protein surface domains

The conformation of water around proteins is of paramount importance, as it determines protein interactions. Although the average water properties around the surface of proteins have been provided experimentally and computationally, protein surfaces are highly heterogeneous. Therefore, it is crucial to determine the correlations of water to the local distributions of polar and nonpolar protein surface domains to understand functions such as aggregation, mutations, and delivery. By using atomistic simulations, we investigate the orientation and dynamics of water molecules next to 4 types of protein surface domains: negatively charged, positively charged, and charge-neutral polar and nonpolar amino acids. The negatively charged amino acids orient around 98% of the neighboring water dipoles toward the protein surface, and such correlation persists up to around 16 Å from the protein surface. The positively charged amino acids orient around 94% of the nearest water dipoles against the protein surface, and the correlation persists up to around 12 Å. The charge-neutral polar and nonpolar amino acids are also orienting the water neighbors in a quantitatively weaker manner. A similar trend was observed in the residence time of the nearest water neighbors. These findings hold true for 3 technically important enzymes (PETase, cytochrome P450, and organophosphorus hydrolase). Our results demonstrate that the water-amino acid degree of correlation follows the same trend as the amino acid contribution in proteins solubility, namely, the negatively charged amino acids are the most beneficial for protein solubility, then the positively charged amino acids, and finally the charge-neutral amino acids.

A Recombinant Affinity Reagent Specific for a Phosphoepitope of Akt1

The serine/threonine-protein kinase, Akt1, plays an important part in mammalian cell growth, proliferation, migration and angiogenesis, and becomes activated through phosphorylation. To monitor phosphorylation of threonine 308 in Akt1, we developed a recombinant phosphothreonine-binding domain (pTBD) that is highly selective for the Akt1 phosphopeptide. A phage-display library of variants of the Forkhead-associated 1 (FHA1) domain of yeast Rad53p was screened by affinity selection to the phosphopeptide, 301-KDGATMKpTFCGTPEY-315, and yielded 12 binding clones. The strongest binders have equilibrium dissociation constants of 160–180 nanomolar and are phosphothreonine-specific in binding. The specificity of one Akt1-pTBD was compared to commercially available polyclonal antibodies (pAbs) generated against the same phosphopeptide. The Akt1-pTBD was either equal to or better than three pAbs in detecting the Akt1 pT308 phosphopeptide in ELISAs.

PD-L1 Nanobody Competitively Inhibits the Formation of the PD-1/PD-L1 Complex: Comparative Molecular Dynamics Simulations

The anti-PD-L1 monoclonal antibody (mAb) targeting PD-1/PD-L1 immune checkpoint has achieved outstanding results in clinical application and has become one of the most popular anti-cancer drugs. The mechanism of molecular recognition and inhibition of PD-L1 mAbs is not yet clear, which hinders the subsequent antibody design and modification. In this work, the trajectories of PD-1/PD-L1 and nanobody/PD-L1 complexes were obtained via comparative molecular dynamics simulations. Then, a series of physicochemical parameters including hydrogen bond, dihedral angle distribution, pKa value and binding free energy, and so forth, were all comparatively analyzed to investigate the recognition difference between PD-L1 and PD-1 and nanobody. Both LR113 (the amino acid residues in PD-L1 are represented by the lower left sign of L) and LR125 residues of PD-L1 undergo significant conformational change after association with mAbs, which dominates a strong electrostatic interaction. Solvation effect analysis revealed that solvent-water enhanced molecular recognition between PD-L1 and nanobody. By combining the analyses of the time-dependent root mean squared fluctuation (RMSF), free energy landscape, clustering and energy decomposition, the potential inhibition mechanism was proposed that the nanobody competitively and specifically bound to the β-sheet groups of PD-L1, reduced the PD-L1’s flexibility and finally blocked the formation of PD-1/PD-L1 complex. Based on the simulation results, site-directed mutagenesis of ND99 (the amino acid residues in Nano are displayed by the lower left sign of N) and NQ116 in the nanobody may be beneficial for improving antibody activity. This work offers some structural guidance for the design and modification of anticancer mAbs based on the structure of the PD-1/PD-L1 complex.

A structural perspective on the interactions of TRAF6 and Basigin during the onset of melanoma: A molecular dynamics simulation study

Metastatic melanoma is the most fatal type of skin cancer. The roles of matrix metalloproteinases (MMPs) have well been established in the onset of melanoma. Basigin (BSG) belongs to the immunoglobulin superfamily and is critical for induction of extracellular MMPs during the onset of various cancers including melanoma. Tumor necrosis factor receptor-associated factor 6 (TRAF6) is an E3-ligase that interacts with BSG and mediates its membrane localization, which leads to MMP expression in melanoma cells. This makes TRAF6 a potential therapeutic target in melanoma. We here conducted protein-protein interaction studies on TRAF6 and BSG to get molecular level insights of the reactions. The structure of human BSG was constructed by protein threading. Molecular-docking method was applied to develop the TRAF6-BSG complex. The refined docked complex was further optimized by molecular dynamics simulations. Results from binding free energy, surface properties, and electrostatic interaction analysis indicate that Lys340 and Glu417 of TRAF6 play as the anchor residues in the protein interaction interface. The current study will be helpful in designing specific modulators of TRAF6 to control melanoma metastasis.