Electrostatic Complementarity Research Articles

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

Discovery and Computational Studies of Potent Covalent Kinase Inhibitors with α-Substituent Electrophiles Targeting Cysteine.

Both reversible noncovalent inhibitors and irreversible covalent inhibitors targeting tyrosine kinases have their disadvantages. The reversible covalent inhibitors with electrophilic group cyanoacrylamide as warheads reacting with cysteine residues could solve the dilemmas. However, there are still several unresolved issues regarding the electrophilic groups. In this manuscript, a series of EGFR inhibitors with double electron-withdrawing substituents introduced into the Cα position on the olefin bond were designed and synthesized. The binding structures and characteristics of inhibitors with the kinase in both the first noncovalent binding phase and the second covalent binding step were explored and combined with molecular docking and molecular dynamics simulations. Then, the reverse β-elimination reactions of the thiol-Michael adducts were investigated by applying density functional theory calculations. In addition, the effects of different electrophilic substituents of Cα on the binding between the inhibitors and kinase were elucidated. The results suggested that the electrophilicity and size of the electron-withdrawing groups play an important role in the specific interactions during the reaction. The compounds with the electron-withdrawing groups that had medium electrostatic and steric complementarity to the kinase active site could cooperatively stabilize the complexes and showed relatively good potent activities in the kinase assay experiment. The mechanical and structural information in this study could enhance our understanding of the functioning of the electron-withdrawing groups in the covalent inhibitors. The results might help to design efficient cysteine targeting inhibitors in the future.

Identification of Potential Druggable Targets and Structure-Based Virtual Screening for Drug-like Molecules against the Shrimp Pathogen Enterocytozoon hepatopenaei.

Enterocytozoon hepatopenaei (EHP) causes slow growth syndrome in shrimp, resulting in huge economic losses for the global shrimp industry. Despite worldwide reports, there are no effective therapeutics for controlling EHP infections. In this study, five potential druggable targets of EHP, namely, aquaporin (AQP), cytidine triphosphate (CTP) synthase, thymidine kinase (TK), methionine aminopeptidase2 (MetAP2), and dihydrofolate reductase (DHFR), were identified via functional classification of the whole EHP proteome. The three-dimensional structures of the proteins were constructed using the artificial-intelligence-based program AlphaFold 2. Following the prediction of druggable sites, the ZINC15 and ChEMBL databases were screened against targets using docking-based virtual screening. Molecules with affinity scores ≥ 7.5 and numbers of interactions ≥ 9 were initially selected and subsequently enriched based on their ADMET properties and electrostatic complementarities. Five compounds were finally selected against each target based on their complex stabilities and binding energies. The compounds CHEMBL3703838, CHEMBL2132563, and CHEMBL133039 were selected against AQP; CHEMBL1091856, CHEMBL1162979, and CHEMBL525202 against CTP synthase; CHEMBL4078273, CHEMBL1683320, and CHEMBL3674540 against TK; CHEMBL340488, CHEMBL1966988, and ZINC000828645375 against DHFR; and CHEMBL3913373, ZINC000016682972, and CHEMBL3142997 against MetAP2.The compounds exhibited high stabilities and low binding free energies, indicating their abilities to suppress EHP infections; however, further validation is necessary for determining their efficacy.

Energy decomposition and waterswapping analysis to investigate the SNP associated DPD mediated 5-FU resistance

ABSTRACT 5-fluorouracil is an essential component of systemic chemotherapy for colon, breast, head, and neck cancer patients. However, tumoral overexpression of the dihydropyrimidine dehydrogenase has rendered 5-FU clinically ineffective by inactivating it to 5'-6'-dihydro fluorouracil. The responses to 5-FU in terms of efficacy and toxicity greatly differ depending upon the population group, because of variability in the DPD activity levels. In the current study, key active site amino acids involved in the 5-FU inactivation were investigated by modelling the 3D structure of human DPD in a complex with 5-FU. The identified amino acids were analyzed for their possible missense mutations available in dbSNP database. Out of 12 missense SNPs, four were validated either by sequencing in the 1000 Genomes project or frequency/genotype data. The recorded validated missense SNPs were further considered to analyze the effect of their respective alterations on 5-FU binding. Overall findings suggested that population bearing the Glu611Val DPD mutation (rs762523739) is highly vulnerable to 5-FU resistance. From the docking, electrostatic complementarity, dynamics, and energy decomposition analyses it was found that the above mutation showed superior scores than the wild DPD -5FU complex. Therefore, prescribing prodrug NUC-3373 or DPD inhibitors (Gimeracil/3-Cyano-2,6-Dihydroxypyridines) as adjuvant therapy may overcome the 5-FU resistance.

3D-QSAR assisted identification of selective CYP1B1 inhibitors: an effective bioisosteric replacement/molecular docking/electrostatic complementarity analysis.

Cytochrome P450-1B1 is a majorly overexpressed drug-metabolizing enzyme in tumors and is responsible for inactivation and subsequent resistance to a variety of anti-cancer drugs, i.e., docetaxel, tamoxifen, and cisplatin. In the present study, a 3D quantitative structure-activity relationship (3D-QSAR) model has been constructed for the identification, design, and optimization of novel CYP1B1 inhibitors. The model has been built using a set of 148 selective CYP1B1 inhibitors. The developed model was evaluated based on certain statistical parameters including q2 and r2 which showed the acceptable predictive and descriptive capability of the generated model. The developed 3D-QSAR model assisted in understanding the key molecular fields which were firmly related to the selective CYP1B1 inhibition. A theoretic approach for the generation of new lead compounds with optimized CYP1B1 receptor affinity has been performed utilizing bioisosteric replacement analysis. These generated molecules were subjected to a developed 3D-QSAR model to predict the inhibitory activity potentials. Furthermore, these compounds were scrutinized through the activity atlas model, molecular docking, electrostatic complementarity, molecular dynamics, and waterswap analysis. The final hits might act as selective CYP1B1 inhibitors which could address the issue of resistance. This 3D-QSAR includes several chemically diverse selective CYP1B1 receptor ligands and well accounts for the individual ligand's inhibition affinities. These features of the developed 3D-QSAR model will ensure future prospective applications of the model to speed up the identification of new potent and selective CYP1B1 receptor ligands.

A Vaccine Strategy Based on the Identification of an Annular Ganglioside Binding Motif in Monkeypox Virus Protein E8L.

The recent outbreak of Monkeypox virus requires the development of a vaccine specifically directed against this virus as quickly as possible. We propose here a new strategy based on a two-step analysis combining (i) the search for binding domains of viral proteins to gangliosides present in lipid rafts of host cells, and (ii) B epitope predictions. Based on previous studies of HIV and SARS-CoV-2 proteins, we show that the Monkeypox virus cell surface-binding protein E8L possesses a ganglioside-binding motif consisting of several subsites forming a ring structure. The binding of the E8L protein to a cluster of gangliosides GM1 mimicking a lipid raft domain is driven by both shape and electrostatic surface potential complementarities. An induced-fit mechanism unmasks selected amino acid side chains of the motif without significantly affecting the secondary structure of the protein. The ganglioside-binding motif overlaps three potential linear B epitopes that are well exposed on the unbound E8L surface that faces the host cell membrane. This situation is ideal for generating neutralizing antibodies. We thus suggest using these three sequences derived from the E8L protein as immunogens in a vaccine formulation (recombinant protein, synthetic peptides or genetically based) specific for Monkeypox virus. This lipid raft/ganglioside-based strategy could be used for developing therapeutic and vaccine responses to future virus outbreaks, in parallel to existing solutions.