Electrostatic Complementarity Research Articles

Bluues_cplx: Electrostatics at Protein-Protein and Protein-Ligand Interfaces.

(1) Background: Electrostatics plays a capital role in protein-protein and protein-ligand interactions. Implicit solvent models are widely used to describe electrostatics and complementarity at interfaces. Electrostatic complementarity at the interface is not trivial, involving surface potentials rather than the charges of surfacial contacting atoms. (2) Results: The program bluues_cplx, here used in conjunction with the software NanoShaper to compute molecular surfaces, has been used to compute the electrostatic properties of 756 protein-protein and 189 protein-ligand complexes along with the corresponding isolated molecules. (3) Methods: The software we make available here uses Generalized Born (GB) radii, computed by a molecular surface integral, to output several descriptors of electrostatics at protein (and in general, molecular) interfaces. We illustrate the usage of the software by analyzing a dataset of protein-protein and protein-ligand complexes, thus extending and refining previous analyses of electrostatic complementarity at protein interfaces. (4) Conclusions: The complete analysis of a molecular complex is performed in tens of seconds on a PC, and the results include the list of surfacial contacting atoms, their charges and Pearson correlation coefficient, the list of contacting surface points with the electrostatic potential (computed for the isolated molecules) and Pearson correlation coefficient, the electrostatic and hydrophobic free energy with different contributions for the isolated molecules, their complex and the difference for all terms. The software is readily usable for any molecular complex in solution.

The Determination of Molecular Electrostatic Potential and Anticancer Properties of Eugenol: A Theoretical Study

Cancer is one of the deadliest diseases worldwide, and for this reason, it is a prominent field of study in drug development. It has been reported in various studies that some of the plants and essential oils obtained from plants have high anticancer activities. This situation is related to the compound groups found in plants and essential oils. Studies on using essential oils in combination with synthetic drugs or aromatherapy are ongoing. Essential oils show cytotoxic properties and may play a role in the death of cancer cells. Eugenol is an important compound found in clove, laurel, and cinnamon essential oils that has anticancer activity in various types of cancer. Eugenol has the ability to reduce cyclooxygenase-2 (COX-2) activity and to inhibit cell proliferation through NF-κB suppression in various types of cancer. In this study, the binding profiles of eugenol with COX-2 and Human inhibitor of nuclear transcription factor κB (IkB) kinase beta, which plays a crucial role in the NF-κB signaling pathway, were examined by molecular docking study, which is one of the methods used in computer-aided drug design. A supporting study was performed to understand the electrostatic complementarity between ligand and receptor by molecular electrostatic potential (MEP) analysis. As a result of the study, it was comparatively presented that eugenol has similar interaction profiles with reference compounds.

Determining the Electrostatic Contributions of GTPase-GEF Complexes on Interfacial Drug Binding Specificity: A Case Study of a Protein-Drug-Protein Complex.

Understanding the factors that contribute to specificity of protein-protein interactions allows for design of orthosteric small molecules. Within this environment, a small molecule requires both structural and electrostatic complementarity. While the structural contribution to protein-drug-protein specificity is well characterized, electrostatic contributions require more study. To this end, we used a series of protein complexes involving Arf1 bound to guanine nucleotide exchange factors (GEFs) that are sensitive or resistant to the small molecule brefeldin A (BFA). By comparing BFA-sensitive Arf1-Gea1p and Arf1-ARNO with different combinations of four BFA sensitizing ARNO mutations (ARNOwt, ARNO1M, ARNO3M, and ARNO4M), we describe how electrostatic environments at each interface guide BFA binding specificity. We labeled Arf1 with cyanocysteine at several interfacial sites and measured by nitrile adsorption frequencies to map changes in electric field at each interface using the linear Stark equation. Temperature dependence of nitrile vibrational spectra was used to investigate differences in hydrogen bonding environments. These comparisons showed that interfacial electric field at the surface of Arf1 varied substantially depending on the GEF. The greatest differences were seen between Arf1-ARNOwt and Arf1-ARNO4M, suggesting a greater change in electric field is required for BFA binding to Arf1-ARNO. Additionally, rigidity of the interface of the Arf1-ARNO complex correlated strongly with BFA sensitivity, indicating that flexible interfaces are sensitive to disruption upon orthosteric small molecule binding. These findings demonstrate a qualitatively consistent electrostatic environment for Arf1 binding and more subtle differences preventing BFA specificity. We discuss how these results will guide improved design of other small molecules that can target protein-protein interfaces.

Conformational disorder in quercetin dihydrate revealed from ultrahigh-resolution synchrotron diffraction.

Quercetin, a bioflavonoid abundant in plants, boasts antioxidant properties and plays a crucial role in various biological systems. The diffraction data of a quercetin dihydrate crystal have been measured at 20 (2) K to ultrahigh resolution (0.30 Å) using a synchrotron X-ray source. After meticulous multipolar refinement of the charge density, Fourier residual electron density peaks were identified, particularly at the position of hydrogen atom H15 of the catechol ring. This observation revealed a subtle disorder in the molecule, prompting the modelling of the catechol ring in two positions with occupancy percentages of 98.4% and 1.6% in the anti and syn conformations, respectively. Intermolecular interactions are analysed using Hirshfeld fingerprint plots and enrichment ratios. With the presence of numerous O-H...O hydrogen bonds, the packing shows good electrostatic complementarity between the quercetin molecule and its surroundings. The parallel displaced stacking interaction between two anti-quercetin molecules related by a translation along the a axis is, however, not attractive for its electrostatic contribution. The syn conformation shows more attractive quercetin dimers than the anti one. On the other hand, electrostatic interactions between quercetin and the two water molecules are stronger in the anti conformation. The electrostatic interactions of quercetin with human inositol polyphosphate multikinase were analysed in the structure of the complex found in the Protein Data Bank and compared with those the take place in the quercetin crystal packing.

Integrated computational approaches for identification of potent pyrazole-based glycogen synthase kinase-3β (GSK-3β) inhibitors: 3D-QSAR, virtual screening, docking, MM/GBSA, EC, MD simulation studies.

Glycogen synthase kinase-3β (GSK-3β) has emerged as a crucial target due to its substantial contribution in various cellular processes. Dysfunctional GSK-3β activity can lead to ion channel disturbances, sustain abnormal excitability, and contribute to the pathogenesis of epilepsy and other GSK-3β-related disorders. A set of 82 pyrazole analogs was utilized to study its structural features using a three-dimensional quantitative structure-activity relationship (3D-QSAR), virtual screening, molecular docking, and molecular dynamics. The QSAR model, validated using internal and external methods, demonstrated robustness with a high correlation coefficient r2training = 0.99, cross-validation coefficient q2 = 0.79, r2test = 0.69, and r2external = 0.74. The "Average of Actives" in the Activity Atlas model identified 17 molecules as active. Subsequent pharmacophore-based virtual screening of 17 actives yielded 70 compounds, which were selected as the prediction set to determine the potential GSK-3β inhibitors. Docking studies pinpointed compound P66 as the promising lead compound, with a docking score of - 10.555kcal/mol. These findings were further supported by electrostatic potential (ESP), electrostatic complementarity (EC), and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) analyses. Furthermore, a 500ns molecular dynamics (MD) simulation confirmed the structural and conformational stability of the lead complex throughout the simulation period. As a result, this study suggests that compound P66 holds the potential to be a potent lead candidate for the inhibition of GSK-3β, offering a novel therapeutic approach for GSK-3β related disorders, including epilepsy.

Combining Complementarity and Binding Energetics in the Assessment of Protein Interactions: EnCPdock-A Practical Manual.

The combined effect of shape and electrostatic complementarities (Sc, EC) at the interface of the interacting protein partners (PPI) serves as the physical basis for such associations and is a strong determinant of their binding energetics. EnCPdock (https://www.scinetmol.in/EnCPdock/) presents a comprehensive web platform for the direct conjoint comparative analyses of complementarity and binding energetics in PPIs. It elegantly interlinks the dual nature of local (Sc) and nonlocal complementarity (EC) in PPIs using the complementarity plot. It further derives an AI-based ΔGbinding with a prediction accuracy comparable to the state of the art. This book chapter presents a practical manual to conceptualize and implement EnCPdock with its various features and functionalities, collectively having the potential to serve as a valuable protein engineering tool in the design of novel protein interfaces.

The Molecular Theory of Sweet Taste: Revisit, Update, and Beyond.

The molecular origin of the sweet taste is still elusive. Herein, the canonical AH-B-X theory of sweet taste is extended by resurveying various sweeteners, which provides deeper insights into an analogous intramolecular connectivity pattern of both glucophores in sweeteners and their interaction counterparts in sweet taste receptor TAS1R2/TAS1R3: electrostatic complementarity and topochemical compatibility. Furthermore, their complementary interaction is elaborately illustrated, accounting for the common molecular feature of eliciting sweetness. Moreover, it highlights that multiple glucophores in a topological system synergistically mediate the elicitation and performance of sweetness. This perspective presents a meaningful framework for the structure-activity relationship-based molecular design and modification of sweeteners and sheds light on the mechanism of molecular evolution of TAS1R2s/TAS1R3s. The link between palatability of sweeteners and harmony relationships between their structural components via stereochemistry and network has significant implications to illuminate the underlying mechanisms by which nature designs chemical reactions to elicit the most important taste.

Novel Omicron Variants Enhance Anchored Recognition of TMEM106B: A New Pathway for SARS-CoV-2 Cellular Invasion.

The recent discovery that TMEM106B serves as a receptor mediating ACE2-independent SARS-CoV-2 entry into cells deserves attention, especially in the background of the frequent emergence of mutant strains. Here, the structure-dynamic features of this novel pathway are dissected deeply. Our investigation revealed that the large loop (RBD@471-491) could anchor TMEM106B, which was then firmly locked by another loop (RBD@444-451). The novel and widely disseminated Omicron variants (BA.2.86/EG.5.1) affect the anchoring recognition of proteins, with BA.2.86 being more likely to impact cells with limited or undetectable ACE2 expression. The large loop of the EG.5.1 variant captures TMEM106B poorly due to impaired electrostatic complementarity. Furthermore, we emphasize that antibody design against these two loops could enhance the protection of ACE2 low-expressing cells according to the alanine scanning mutagenesis of multiple antibodies. We hope this study will provide a novel perspective for the prevention and treatment against this new viral invasion pathway.

Structure of the murine CD94-NKG2A receptor in complex with Qa-1b presenting an MHC-I leader peptide.

The heterodimeric natural killer cells antigen CD94 (CD94)-NKG2-A/NKG2-B type II integral membrane protein (NKG2A) receptor family expressed on human and mouse natural killer (NK) cells monitors global major histocompatibility complex (MHC) class I cell surface expression levels through binding to MHC class Ia-derived leader sequence peptides presented by HLA class I histocompatibility antigen, alpha chain E (HLA-E; in humans) or H-2 class I histocompatibility antigen, D-37 (Qa-1b; in mice). Although the molecular basis underpinning human CD94-NKG2A recognition of HLA-E is known, the equivalent interaction in the murine setting is not. By determining the high-resolution crystal structure of murine CD94-NKG2A in complex with Qa-1b presenting the Qa-1 determinant modifier peptide (QDM), we resolved the mode of binding. Compared to the human homologue, the murine CD94-NKG2A-Qa-1b-QDM displayed alterations in the distribution of interactions across CD94 and NKG2A subunits that coincide with differences in electrostatic complementarity of the ternary complex and the lack of cross-species reactivity. Nevertheless, we show that Qa-1b could be modified through W65R + N73I mutations to mimic HLA-E, facilitating binding with both human and murine CD94-NKG2A. These data underscore human and murine CD94-NKG2A cross-species heterogeneity and provide a foundation for humanising Qa-1b in immune system models.

Unprecedented selectivity for homologous lectin targets: differential targeting of the viral receptors L-SIGN and DC-SIGN.

DC-SIGN (CD209) and L-SIGN (CD209L) are two C-type lectin receptors (CLRs) that facilitate SARS-CoV-2 infections as viral co-receptors. SARS-CoV-2 manipulates both DC-SIGN and L-SIGN for enhanced infection, leading to interest in developing receptor antagonists. Despite their structural similarity (82% sequence identity), they function differently. DC-SIGN, found in dendritic cells, shapes the immune response by recognizing pathogen-associated carbohydrate patterns. In contrast, L-SIGN, expressed in airway epithelial endothelial cells, is not directly involved in immunity. COVID-19's primary threat is the hyperactivation of the immune system, potentially reinforced if DC-SIGN engages with exogenous ligands. Therefore, L-SIGN, co-localized with ACE2-expressing cells in the respiratory tract, is a more suitable target for anti-adhesion therapy. However, designing a selective ligand for L-SIGN is challenging due to the high sequence identity of the Carbohydrate Recognition Domains (CRDs) of the two lectins. We here present Man84, a mannose ring modified with a methylene guanidine triazole at position 2. It binds L-SIGN with a K D of 12.7μM ± 1 μM (ITC) and is the first known L-SIGN selective ligand, showing 50-fold selectivity over DC-SIGN (SPR). The X-ray structure of the L-SIGN CRD/Man84 complex reveals the guanidinium group's role in achieving steric and electrostatic complementarity with L-SIGN. This allows us to trace the source of selectivity to a single amino acid difference between the two CRDs. NMR analysis confirms the binding mode in solution, highlighting Man84's conformational selection upon complex formation. Dimeric versions of Man84 achieve additional selectivity and avidity in the low nanomolar range. These compounds selectively inhibit L-SIGN dependent trans-infection by SARS-CoV-2 and Ebola virus. Man84 and its dimeric constructs display the best affinity and avidity reported to date for low-valency glycomimetics targeting CLRs. They are promising tools for competing with SARS-CoV-2 anchoring in the respiratory tract and have potential for other medical applications.

Electrostatic Impact of Brefeldin A on Thiocyanate Probes Surrounding the Interface of Arf1-BFA-ARNO4M, a Protein-Drug-Protein Complex.

Protein-protein interactions regulate many cellular processes, making them ideal drug candidates. Design of such drugs, however, is hindered by a lack of understanding of the factors that contribute to the interaction specificity. Specific protein-protein complexes possess both structural and electrostatic complementarity, and while structural complementarity of protein complexes has been extensively investigated, fundamental understanding of the complicated networks of electrostatic interactions at these interfaces is lacking, thus hindering the rational design of orthosterically binding small molecules. To better understand the electrostatic interactions at protein interfaces and how a small molecule could contribute to and fit within that environment, we used a model protein-drug-protein system, Arf1-BFA-ARNO4M, to investigate how small molecule brefeldin A (BFA) perturbs the Arf1-ARNO4M interface. By using nitrile probe labeled Arf1 sites and measuring vibrational Stark effects as well as temperature dependent infrared shifts, we measured changes in the electric field and hydrogen bonding at this interface upon BFA binding. At all five probe locations of Arf1, we found that the vibrational shifts resulting from BFA binding corroborate trends found in Poisson-Boltzmann calculations of surface potentials of Arf1-ARNO4M and Arf1-BFA-ARNO4M, where BFA contributes negative electrostatic potential to the protein interface. The data also corroborate previous hypotheses about the mechanism of interfacial binding and confirm that alternating patches of hydrophobic and polar interactions lead to BFA binding specificity. These findings demonstrate the impact of BFA on this protein-protein interface and have implications for the design of other interfacial drug candidates.

Chemical Complementarity of Tumor Resident, Adaptive Immune Receptor CDR3s and Previously Defined Hepatitis C Virus Epitopes Correlates with Improved Outcomes in Hepatocellular Carcinoma.

To better understand how adaptive immune receptors (IRs) in hepatocellular carcinoma (HCC) microenvironments are related to disease outcomes, we employed a chemical complementarity scoring algorithm to quantify electrostatic complementarity between HCC tumor TRB or IGH complementarity-determining region 3 (CDR3) amino acid (AA) sequences and previously characterized hepatitis C virus (HCV) epitopes. High electrostatic complementarity between HCC-resident CDR3s and 12 HCV epitopes was associated with greater survival probabilities, as indicated by two distinct HCC IR CDR3 datasets. Two of the HCV epitopes, HCV*71871 (TRB) and HCV*13458 (IGH), were also determined to represent significantly larger electrostatic CDR3-HCV epitope complementarity in HCV-positive HCC cases, compared with HCV-negative HCC cases, with the CDR3s representing yet a third, independent HCC dataset. Overall, the results indicated the utility of CDR3 AA sequences as biomarkers for HCC patient stratification and as potential guides for the development of therapeutic reagents.

Breakdown of Arabidopsis thaliana thioredoxins and glutaredoxins based on electrostatic similarity-Leads to common and unique interaction partners and functions.

The reversible reduction and oxidation of protein thiols was first described as mechanism to control light/dark-dependent metabolic regulation in photosynthetic organisms. Today, it is recognized as an essential mechanism of regulation and signal transduction in all kingdoms of life. Proteins of the thioredoxin (Trx) family, Trxs and glutaredoxins (Grxs) in particular, catalyze thiol-disulfide exchange reactions and are vital players in the operation of thiol switches. Various Trx and Grx isoforms are present in all compartments of the cell. These proteins have a rather broad but at the same time distinct substrate specificity. Understanding the molecular basis of their target specificity is central to the understanding of physiological and pathological redox signaling. Electrostatic complementarity of the redoxins with their target proteins has been proposed as a major reason. Here, we analyzed the electrostatic similarity of all Arabidopsis thaliana Trxs, Grxs, and proteins containing such domains. Clustering of the redoxins based on this comparison suggests overlapping and also distant target specificities and thus functions of the different sub-classes including all Trx isoforms as well as the three classes of Grxs, i.e. CxxC-, CGFS-, and CC-type Grxs. Our analysis also provides a rationale for the tuned substrate specificities of both the ferredoxin- and NADPH-dependent Trx reductases.

Molecular reshaping of phage-displayed Interleukin-2 at beta chain receptor interface to obtain potent super-agonists with improved developability profiles

Interleukin-2 (IL-2) engineered versions, with biased immunological functions, have emerged from yeast display and rational design. Here we reshaped the human IL-2 interface with the IL-2 receptor beta chain through the screening of phage-displayed libraries. Multiple beta super-binders were obtained, having increased receptor binding ability and improved developability profiles. Selected variants exhibit an accumulation of negatively charged residues at the interface, which provides a better electrostatic complementarity to the beta chain, and faster association kinetics. These findings point to mechanistic differences with the already reported superkines, characterized by a conformational switch due to the rearrangement of the hydrophobic core. The molecular bases of the favourable developability profile were tracked to a single residue: L92. Recombinant Fc-fusion proteins including our variants are superior to those based on H9 superkine in terms of expression levels in mammalian cells, aggregation resistance, stability, in vivo enhancement of immune effector responses, and anti-tumour effect.

Thermochromic Behavior of Polydiacetylene Nanomaterials Driven by Charged Peptide Amphiphiles.

The tunability of chromatic phases adapted by chromogenic polymers such as polydiacetylene (PDA) is key to their utility for robust sensing applications. Here, we investigated the influence of charged peptide interactions on the structure-dependent thermochromicity of amphiphilic PDAs. Solid-state NMR and circular dichroism analyses show that our oppositely charged peptide-PDA samples have distinct degrees of structural order, with the coassembled sample being in between the β-sheet-like positive peptide-PDA and the relatively disordered negative peptide-PDA. All solutions exhibit thermochromicity between 20 and 80 °C, whereby the hysteresis of the blue, planar phase is much larger than that of the red, twisted phase. Resonance Raman spectroscopy of films demonstrates that only coassemblies with electrostatic complementarity stabilize coexisting blue and red PDA phases. This work reveals the nature of the structural changes responsible for the thermally responsive chromatic transitions of biomolecule-functionalized polymeric materials and how this process can be directed by sequence-dictated electrostatic interactions.

Differences in the organization of interface residues tunes the stability of the SARS-CoV-2 spike-ACE2 complex.

The continuous emergence of novel variants represents one of the major problems in dealing with the SARS-CoV-2 virus. Indeed, also due to its prolonged circulation, more than ten variants of concern emerged, each time rapidly overgrowing the current viral version due to improved spreading features. As, up to now, all variants carry at least one mutation on the spike Receptor Binding Domain, the stability of the binding between the SARS-CoV-2 spike protein and the human ACE2 receptor seems one of the molecular determinants behind the viral spreading potential. In this framework, a better understanding of the interplay between spike mutations and complex stability can help to assess the impact of novel variants. Here, we characterize the peculiarities of the most representative variants of concern in terms of the molecular interactions taking place between the residues of the spike RBD and those of the ACE2 receptor. To do so, we performed molecular dynamics simulations of the RBD-ACE2 complexes of the seven variants of concern in comparison with a large set of complexes with different single mutations taking place on the RBD solvent-exposed residues and for which the experimental binding affinity was available. Analyzing the strength and spatial organization of the intermolecular interactions of the binding region residues, we found that (i) mutations producing an increase of the complex stability mainly rely on instaurating more favorable van der Waals optimization at the cost of Coulombic ones. In particular, (ii) an anti-correlation is observed between the shape and electrostatic complementarities of the binding regions. Finally, (iii) we showed that combining a set of dynamical descriptors is possible to estimate the outcome of point mutations on the complex binding region with a performance of 0.7. Overall, our results introduce a set of dynamical observables that can be rapidly evaluated to probe the effects of novel isolated variants or different molecular systems.

Electrostatic complementarity at the interface drives transient protein-protein interactions

Understanding the mechanisms driving bio-molecules binding and determining the resulting complexes’ stability is fundamental for the prediction of binding regions, which is the starting point for drug-ability and design. Characteristics like the preferentially hydrophobic composition of the binding interfaces, the role of van der Waals interactions, and the consequent shape complementarity between the interacting molecular surfaces are well established. However, no consensus has yet been reached on the role of electrostatic. Here, we perform extensive analyses on a large dataset of protein complexes for which both experimental binding affinity and pH data were available. Probing the amino acid composition, the disposition of the charges, and the electrostatic potential they generated on the protein molecular surfaces, we found that (i) although different classes of dimers do not present marked differences in the amino acid composition and charges disposition in the binding region, (ii) homodimers with identical binding region show higher electrostatic compatibility with respect to both homodimers with non-identical binding region and heterodimers. Interestingly, (iii) shape and electrostatic complementarity, for patches defined on short-range interactions, behave oppositely when one stratifies the complexes by their binding affinity: complexes with higher binding affinity present high values of shape complementarity (the role of the Lennard-Jones potential predominates) while electrostatic tends to be randomly distributed. Conversely, complexes with low values of binding affinity exploit Coulombic complementarity to acquire specificity, suggesting that electrostatic complementarity may play a greater role in transient (or less stable) complexes. In light of these results, (iv) we provide a novel, fast, and efficient method, based on the 2D Zernike polynomial formalism, to measure electrostatic complementarity without the need of knowing the complex structure. Expanding the electrostatic potential on a basis of 2D orthogonal polynomials, we can discriminate between transient and permanent protein complexes with an AUC of the ROC of sim 0.8. Ultimately, our work helps shedding light on the non-trivial relationship between the hydrophobic and electrostatic contributions in the binding interfaces, thus favoring the development of new predictive methods for binding affinity characterization.

Ligand and Gold(I) Fluorescein-AIEgens as Photosensitizers in Solution and Doped Polymers.

The synthesis of fluorescein propargyl diether (L) and two different dinuclear gold(I) derivatives containing a water-soluble phosphane [1,3,5-triaza-7-phosphatricyclo[3.3.1.13.7]decane (PTA) for complex 1 and 3,7-diacetyl-1,3,7-triaza-5-phosphabicyclo[3.3.1]nonane (DAPTA) for complex 2] has been successfully performed. All compounds display intrinsic emission from fluorescein, being less intense for gold(I) complexes due to the heavy-atom effect. All compounds aggregate in acetonitrile/water mixtures with the formation of larger aggregates for those samples containing more water content, as evidenced by dynamic light scattering and small-angle X-ray scattering experiments, in agreement with the absorption and emission data. The emission of the samples increases when they are used to obtain luminescent materials with four different organic matrices [poly(methyl methacrylate, polystyrene (PS), cellulose, and Zeonex]. The compounds display very high values of singlet oxygen (1O2) production in dichloromethane. Singlet oxygen production was also evaluated in the doped matrices, being the highest in PS and with an exciting increase on PS microspheres. Density functional theory (BP86-D3) and GFN2-xTB calculations were used to model the assembly of L and complexes 1 and 2 with the different organic matrices and rationalize the experimental findings based on the geometries, molecular electrostatic potential surfaces, and complementarity and HOMO-LUMO gaps.

Electrostatic Complementarities of Glioblastoma-Resident T-Cell Receptors and Cancer Testis Antigens Linked to Poor Outcomes and High Levels of Sphingosine Kinase-2 Expression.

Glioblastoma (GBM) is the most aggressive primary brain tumor in adults. Despite a growing understanding of glioblastoma pathology, the prognosis remains poor. In this study, we used a previously extensively benchmarked algorithm to retrieve immune receptor (IR) recombination reads from GBM exome files available from the cancer genome atlas. The T-cell receptor complementarity determining region-3 (CDR3) amino acid sequences that represent the IR recombination reads were assessed and used for the generation of chemical complementarity scores (CSs) that represent potential binding interactions with cancer testis antigens (CTAs), which is an approach particularly suited to a big data setting. The electrostatic CSs representing the TRA and TRB CDR3s and the CTAs, SPAG9, GAGE12E, and GAGE12F, indicated that an increased electrostatic CS was associated with worse disease-free survival (DFS). We also assessed the RNA expression of immune marker genes, which indicated that a high-level expression of SPHK2 and CIITA genes also correlated with high CSs and worse DFS. Furthermore, apoptosis-related gene expression was revealed to be lower when the TCR CDR3-CTA electrostatic CSs were high. Adaptive IR recombination reads from exome files have the potential to aid in GBM prognoses and may provide opportunities to detect unproductive immune responses.

Exploiting the co-crystal ligands shape, features and structure-based approaches for identification of SARS-CoV-2 Mpro inhibitors

SARS-CoV-2 enters the host cell through the ACE2 receptor and replicates its genome using an RNA-Dependent RNA Polymerase (RDRP). The functional RDRP is released from pro-protein pp1ab by the proteolytic activity of Main protease (Mpro) which is encoded within the viral genome. Due to its vital role in proteolysis of viral polyprotein chains, it has become an attractive potential drug target. We employed a hierarchical virtual screening approach to identify small synthetic protease inhibitors. Statistically optimized molecular shape and color-based features (various functional groups) from co-crystal ligands were used to screen different databases through various scoring schemes. Then, the electrostatic complementarity of screened compounds was matched with the most active molecule to further reduce the hit molecules’ size. Finally, five hundred eighty-seven molecules were docked in Mpro catalytic binding site, out of which 29 common best hits were selected based on Glide and FRED scores. Five best-fitting compounds in complex with Mpro were subjected to MD simulations to analyze their structural stability and binding affinities with Mpro using MM/GB(PB)SA models. Modeling results suggest that identified hits can act as the lead compounds for designing better active Mpro inhibitors to enhance the chemical space to combat COVID-19.Communicated by Ramaswamy H. Sarma