Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes

Functional Loop Dynamics and Characterization of the Inactive State of the NS2B-NS3 Dengue Protease due to Allosteric Inhibitor Binding

Dengue virus, a flavivirus that causes dengue shock syndrome and dengue hemorrhagic fever, is currently prevalent worldwide. A two-component protease (NS2B-NS3) is essential for maturation, representing an important target for designing anti-flavivirus drugs. Previously, consideration has been centered on developing active-site inhibitors of NS2B-NS3pro. However, the flat and charged nature of its active site renders difficulties in developing inhibitors, suggesting an alternative strategy for identifying allosteric inhibitors. The allosterically sensitive site of the dengue protease is located near Ala125, between the 120s loop and 150s loop. Using atomistic molecular dynamics simulations, we have explored the protease's conformational dynamics upon binding of an allosteric inhibitor. Furthermore, characterization of the inherent flexible loops (71-75s loop, 120s loop, and 150s loop) is carried out for allosteric-inhibitor-bound wild-type and mutant A125C variants and a comparison is performed with its unbound state to extract the structural changes describing the inactive state of the protease. Our study reveals that compared to the unliganded system, the inhibitor-bound system shows large structural changes in the 120s loop and 150s loop in contrast to the rigid 71-75s loop. The unliganded system shows a closed-state pocket in contrast to the open state for the wild-type complex that locks the protease into the open and inactive-state conformations. However, the mutant complex fluctuates between open and closed states. Also, we tried to see how mutation and binding of an allosteric inhibitor perturb the connectivity in a protein structure network (PSN) at contact levels. Altogether, our study reveals the mechanism of conformational rearrangements of loops at the molecular level, locking the protein in an inactive conformation, which may be useful for developing allosteric inhibitors.

In Silico Screening of Quorum Sensing Inhibitor Candidates Obtained by Chemical Similarity Search

Quorum sensing (QS) is a bacterial communication using signal molecules, by which they sense population density of their own species, leading to group behavior such as biofilm formation and virulence. Autoinducer-2 (AI2) is a QS signal molecule universally used by both gram-positive and gram-negative bacteria. Inhibition of QS mediated by AI2 is important for various practical applications, including prevention of gum-disease caused by biofilm formation of oral bacteria. In this research, molecular docking and molecular dynamics (MD) simulations were performed for molecules that are chemically similar to known AI2 inhibitors that might have a potential to be quorum sensing inhibitors. The molecules that form stable complexes with the AI2 receptor protein were found, suggesting that they could be developed as a novel AI2 inhibitors after further in vitro validation. The result suggests that combination of ligand-based drug design and computational methods such as MD simulation, and experimental verification, may lead to development of novel AI inhibitor, with a broad range of practical applications.

A strategy to optimize the peptide-based inhibitors against different mutants of the spike protein of SARS-CoV-2

The SARS-CoV-2 virus has caused high-priority health concerns at a global level. Vaccines have stalled the proliferation of viruses to some extent. Yet, the emergence of newer, potentially more infectious, and dangerous mutants such as Delta and Omicron are among the major challenges in finding a more permanent solution for this pandemic. The effectiveness of antivirals Molnupiravir and Paxlovid, authorized for emergency use by the FDA, are yet to be assessed on a larger population. Patients with a high risk of disease progression have received treatment with antibody-cocktail. Most of the mutations leading to the new lineage of SARS-CoV-2 are found in the spike protein of this virus that plays a key role in facilitating host entry. The current study has investigated how to modify a promising peptide-based inhibitor of spike protein, LCB3, against common mutations, N501Y and K417N in the target protein so that it retains its efficacy against the spike protein. LCB3 being a prototype for protein-based inhibitors is an ideal testing system to learn about protein-based inhibitors. This study proposes the substitutions of amino acid residues of LCB3 inhibitor using a structure-based approach that considers free energy decomposition of residues, the distance between atoms, and the interaction among amino acids. The binding free energy calculations suggest a possible improvement in the binding affinity of existing inhibitor LCB3 to the mutant forms of the S-protein using simple substitutions at specific positions of the inhibitor. This approach, being general, can be used in different inhibitors and other mutations and help in fighting against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Molecular modelling study to discover novel JAK2 signaling pathway inhibitor

The JAK2/STAT signaling cascades facilitates receptor signals which is responsible for cell growth, survival and homeostasis. Ligand binding to JAKs causes phosphorylation other proteins known as STATs, which translocate to the nucleus and regulate transcription of several important proteins. Growth hormone, prolactin and γ-interferon known agonists of JAK STAT receptors, signal to the nucleus by a more direct manner than the receptor tyrosine kinases. Mutations in JAKs may be responsible for immunodeficiency and myeloproliferative disorders because of its important role in cytokine signaling and making the pathway a therapeutic target for various disease. The present study screened Zinc database to find novel JAK2 inhibitors using virtual high throughput screening techniques. Selection of compound for further study was on the basis of docking score, free energy and binding pattern of the compound. Molecular simulation and MM/GBSA free energy was evaluated for the binding interactions and the stability of docked conformations. Several parameters which determine protein ligand interaction like RMSD, RMSF, Rg and binding pattern were observed. Hydrogen bonds (Glu 930, 932 and Asp 994) after 150 ns simulation were observed between identified compound INC000096136346 and it was similar to known inhibitor ruxolitinib. MM/GBSA free energy was comparable to known inhibitor ruxolitinib. ZINC000096136346 qualify Lipinski's rule of five, rule of three, WDI like rule and there is one violation in lead like rule.Communicated by Ramaswamy H. Sarma.

Molecular investigation of the dual inhibition mechanism for targeted P53 regulator MDM2/MDMX inhibitors

Inhibitors that competitively bind MDM2/MDMX can block the inhibition of P53 by MDM2/MDMX and restart its tumor-suppressive effect. Molecular studies targeting MDM2/MDMX inhibitors have always been a hot topic in anticancer drug design. Although numerous inhibitors have been designed previously against MDM2/MDMX, their dual inhibition efficacy has not been demonstrated, and few studies assessed the general causes affecting the dual inhibition of MDM2/MDMX by these inhibitors. Here, molecular dynamics simulations and alanine scanning combined with the interaction entropy method were employed to precisely investigate whether 16 inhibitors could dually inhibit MDM2/MDMX and the similarities and differences in the interaction modes. Thereby addressing the key residue sites affecting dual inhibition. Residues L54/M53, I61/60, M62/61, Y67/66, and V93/92 of MDM2/MDMX, which are in corresponding positions in both protein structures, provide significant conditions for these inhibitors to bind to MDM2/MDMX tightly. In addition, most of these inhibitors prefer to bind MDM2 than MDMX, and residues H96 and I99 in MDM2 are attractive targets for inhibitors, resulting in inhibitors binding to MDM2/MDMX with different affinity. These key residues should be considered in the development of dual inhibitors. For these 16 inhibitors, most have dual inhibitory potential for MDM2/MDMX based on the binding affinity of the complexes. Still, it is questionable whether they can exert excellent dual inhibition considering the assessment of the hot-spots. At least their binding affinity for MDMX is not superior to that for MDM2 due to the difference in energy of the van der Waals interactions at the key sites. Furthermore, based on the analysis of three representative inhibitors (TUZ/HRH and HRQ with different binding preferences for MDM2/MDMX), 3-chloropyridine in TUZ leads to the differential binding affinity between the inhibitor and MDM2/MDMX. It readily forms hydrophobic interactions with the surrounding residues H96 and I99. But this phenomenon does not occur in the TUZ-MDMX system, implying the critical role of residues H96/P95 and I99/L98. And the completely different binding mechanism of HRQ binding to MDM2/MDMX explains its inability to inhibit MDM2 well. Thus, we are cautious about its dual inhibitory ability. Besides, HRH is more prone to strong van der Waals interactions with MDM2 than MDMX whereas its 2-chlorofluorobenzene is detrimental to this. We hope that these findings will provide reliable molecular insights for the screening and optimization of targeting MDM2/MDMX dual inhibitors.

Fisetin 8-C-glucoside as entry inhibitor in SARS CoV-2 infection: molecular modelling study

Coronaviruses are RNA viruses that infect varied species including humans. TMPRSS2 is gateway for SARS CoV-2 entry into the host cell. It causes proteolytic activation of spike protein and discharge of the peptide into host cell. The TMPRSS2 inhibition could be one of the approaches to stop the viral entry, therefore, interaction pattern and binding energies for Fisetin and TMPRSS2 have been explored in the present study. TMPRSS2 peptide was used for homology modelling and then for further study. Molecular docking score and MMGBSA Binding energy of Fisetin was better than Nafamostat, a known inhibitor of TMPRSS2. Post docking MM-GBSA free energy for Fisetin and Nafamostat was -42.78 and -21.11 kcal/mol, respectively. Fisetin forms H bond with Val 25, His 41, Lys 42, Lys 45, Glu 44, Ser186. Nafamostat formed H bonds with Lys 85, Asp 90, Asp 203. RMSD plots of TMPRSS2, TMPRSS2-Fisetin and TMPRSS2-Nafamostat complex showed stable profile with very small fluctuation during entire simulation of 150 ns. Significant decrease in TMPRSS2-Fisetin and TMPRSS2-Nafamostat complex fluctuation occurred around His 41, Glu 44, Gly 136, Ser 186 in RMSF study. During simulation Fisetin interaction was observed with residues Val 25, His 41, Glu 44, Lys 45, Lys 87, Gly 136, Gln 183, Ser 186 likewise interaction of Nafamostat with Lys 85, Asp 90, Asn 163, Asp 203 and Ser 205. Post simulation MM-GBSA free energy was found to be -51.87 ± 4.3 and -48.23 ± 4.39 kcal/mol for TMPRSS2 with Fisetin and Nafamostat, respectively.Communicated by Ramaswamy H. Sarma.

Rational Design Yields Molecular Insights on Leaf-Binding of Anchor Peptides

In times of a constantly growing world population and increasing demand for food, sustainable agriculture is crucial. The rainfastness of plant protection agents is of pivotal importance to reduce the amount of applied nutrients, herbicides, and fungicides. As a result of protective agent wash-off, plant protection is lost, and soils and groundwater are severely polluted. To date, rainfastness of plant protection products has been achieved by adding polymeric adjuvants to the agrochemicals. However, polymeric adjuvants will be regarded as microplastics in the future, and environmentally friendly alternatives are needed. Anchor peptides (APs) are promising biobased and biodegradable adhesion promoters. Although the adhesion of anchor peptides to artificial surfaces, such as polymers, has already been investigated in theory and experimentally, exploiting the adhesion to biological surfaces remains challenging. The complex nature and composition of biological surfaces such as plant leaves and fruit surfaces complicate the generation of accurate models. Here, we present the first detailed three-layered atomistic model of the surface of apple leaves and use it to compute free energy profiles of the adhesion and desorption of APs to and from that surface. Our model is validated by a novel fluorescence-based microtiter plate (MTP) assay that mimics these complex processes and allows for quantifying them. For the AP Macaque Histatin, we demonstrate that aromatic and positively charged amino acids are essential for binding to the waxy apple leaf surface. The established protocols should generally be applicable for tailoring the binding properties of APs to biological interfaces.

A Physics-Guided Neural Network for Predicting Protein–Ligand Binding Free Energy: From Host–Guest Systems to the PDBbind Database

Calculation of protein–ligand binding affinity is a cornerstone of drug discovery. Classic implicit solvent models, which have been widely used to accomplish this task, lack accuracy compared to experimental references. Emerging data-driven models, on the other hand, are often accurate yet not fully interpretable and also likely to be overfitted. In this research, we explore the application of Theory-Guided Data Science in studying protein–ligand binding. A hybrid model is introduced by integrating Graph Convolutional Network (data-driven model) with the GBNSR6 implicit solvent (physics-based model). The proposed physics-data model is tested on a dataset of 368 complexes from the PDBbind refined set and 72 host–guest systems. Results demonstrate that the proposed Physics-Guided Neural Network can successfully improve the “accuracy” of the pure data-driven model. In addition, the “interpretability” and “transferability” of our model have boosted compared to the purely data-driven model. Further analyses include evaluating model robustness and understanding relationships between the physical features.

A Computational Approach Applied to the Study of Potential Allosteric Inhibitors Protease NS2B/NS3 from Dengue Virus

Dengue virus (DENV) is a danger to more than 400 million people in the world, and there is no specific treatment. Thus, there is an urgent need to develop an effective method to combat this pathology. NS2B/NS3 protease is an important biological target due it being necessary for viral replication and the fact that it promotes the spread of the infection. Thus, this study aimed to design DENV NS2B/NS3pro allosteric inhibitors from a matrix compound. The search was conducted using the Swiss Similarity tool. The compounds were subjected to molecular docking calculations, molecular dynamics simulations (MD) and free energy calculations. The molecular docking results showed that two compounds, ZINC000001680989 and ZINC000001679427, were promising and performed important hydrogen interactions with the Asn152, Leu149 and Ala164 residues, showing the same interactions obtained in the literature. In the MD, the results indicated that five residues, Lys74, Leu76, Asn152, Leu149 and Ala166, contribute to the stability of the ligand at the allosteric site for all of the simulated systems. Hydrophobic, electrostatic and van der Waals interactions had significant effects on binding affinity. Physicochemical properties, lipophilicity, water solubility, pharmacokinetics, druglikeness and medicinal chemistry were evaluated for four compounds that were more promising, showed negative indices for the potential penetration of the Blood Brain Barrier and expressed high human intestinal absorption, indicating a low risk of central nervous system depression or drowsiness as the the side effects. The compound ZINC000006694490 exhibited an alert with a plausible level of toxicity for the purine base chemical moiety, indicating hepatotoxicity and chromosome damage in vivo in mouse, rat and human organisms. All of the compounds selected in this study showed a synthetic accessibility (SA) score lower than 4, suggesting the ease of new syntheses. The results corroborate with other studies in the literature, and the computational approach used here can contribute to the discovery of new and potent anti-dengue agents.

Molecular Dynamics Simulation Investigation of the Binding and Interaction of the EphA6-Odin Protein Complex

Protein-protein interaction plays an important role in the development of almost all cells. Elucidating the dynamic binding and affinity of a protein-protein complex is essential for understanding the biological functions of proteins. EphA6 and Odin proteins are members of the Eph (erythropoietin-producing hepatocyte) receptor family and the Anks (ankyrin repeat and sterile α motif domain-containing) family, respectively. Odin significantly functions in regulating endocytosis, degradation, and stability of EphA receptors. In this work, the key residues of the interaction interface were determined through a hydrogen bond, root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), and dynamic correlation analysis of the conventional molecular dynamics (MD) simulations. The calculated standard binding free energy, -7.92 kcal/mol, between EphA6 and Odin is quite consistent with the experimental measurement value, -8.73 kcal/mol. By the combination of several MD simulation techniques, our investigation of the binding process reveals the detailed representative characteristics of the entire binding pathway at the molecular level. Based on the obtained potential of the mean force (PMF) curve, the analysis of the simulation trajectories shows that the residue Arg1013 in the receptor EphA6 is responsible for capturing Asp739 and Asp740 in the ligand Odin during the initial stage of binding. In the later stage of binding, the hydrogen bonds and salt bridges between a series of residues Lys973, Leu1007, Gly1009, His1010, and Arg1012 in the receptor and residues Leu735, Asn736, Asp739, Asp740, and Asp753 in the ligand mainly contribute to the stability of the protein complex. In addition, the specific change process of the receptor-ligand-binding mode is also clarified during the binding process. Our present simulation will promote a deep understanding of the protein-protein interaction, and the identified key interresidue interaction will be theoretical guidance for the design of protein drugs.

Structure-based design and synthesis of a novel long-chain 4''-alkyl ether derivative of EGCG as potent EGFR inhibitor: in vitro and in silico studies

Herein, we report the discovery of a novel long-chain ether derivative of (−)-epigallocatechin-3-gallate (EGCG), a major green tea polyphenol as a potent EGFR inhibitor. A series of 4′′-alkyl EGCG derivatives have been synthesized via regio-selectively alkylating the 4′′ hydroxyl group in the D-ring of EGCG and tested for their antiproliferative activities against high (A431), moderate (HeLa), and low (MCF-7) EGFR-expressing cancer cell lines. The most potent compound, 4′′-C₁₄ EGCG showed the lowest IC₅₀ values across all the tested cell lines. 4′′-C₁₄ EGCG was also found to be significantly more stable than EGCG under physiological conditions (PBS at pH 7.4). Further western blot analysis and imaging data revealed that 4′′-C₁₄ EGCG induced cell death in A431 cells with shrunken nuclei, nuclear fragmentation, membrane blebbing, and increased population of apoptotic cells where BAX upregulation and BCL_XL downregulation were observed. In addition, autophosphorylation of EGFR and its downstream signalling proteins Akt and ERK were markedly inhibited by 4′′-C₁₄ EGCG. MD simulation and the MM/PBSA analysis disclosed the binding mode of 4′′-C₁₄ EGCG in the ATP-binding site of EGFR kinase domain. Taken together, our findings demonstrate that 4′′-C₁₄ EGCG can act as a promising potent EGFR inhibitor with enhanced stability.

Among the synthesized 4′′-alkyl EGCG derivatives, 4′′-C₁₄ EGCG inhibited EGF stimulated phosphorylation of EGFR and its downstream signaling pathways, ERK and Akt. 4′′-C₁₄ EGCG showed significantly improved stability than EGCG and induced apoptosis.

Identification of PRMT5 inhibitors with novel scaffold structures through virtual screening and biological evaluations

Protein arginine methyltransferase 5 (PRMT5), an important member in PRMT family, has been validated as a promising anticancer target. In this study, through the combination of virtual screening and biological experiments, we have identified two PRMT5 inhibitors with novel scaffold structures. Among them, compound Y2431 showed moderate activity with IC₅₀ value of 10.09 μM and displayed good selectivity against other methyltransferases. The molecular docking analysis and molecular dynamics (MD) simulations suggested that the compound occupied the substrate-arginine binding site. Furthermore, Y2431 exhibited anti-proliferative activity to leukemia cells by inducing cell cycle arrest. Overall, the hit compound could provide a novel scaffold for further optimization of small-molecule PRMT5 inhibitors.

Discovery of Histone Deacetylase Inhibitor Using Molecular Modeling and Free Energy Calculations

The histone acetylation-deacetylation at lysine regulates the functions of many cellular proteins. An increased expression of HDAC6 can cause an increased amount of deacetylated histones, which leads to an inhibition of gene expression and has been associated with cancer cell proliferation. The present study screened the ZINC database to find novel HDAC6 inhibitors using virtual high-throughput screening techniques. The docking score, free energy, and binding pattern of the complexes were used to select a best ligand for further study. Molecular dynamic simulations, binding interactions, and the stability of docked conformations were investigated. Several parameters that determine protein-ligand interactions, such as root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), radius of gyration (Rg), and binding pattern, were observed. Hydrogen bonds were observed at His 573 and Gly 582 after a 150 ns simulation with identified compound ZINC000002845205, and they were similar to known inhibitor Panobinostat. The molecular mechanics with generalised Born and surface area solvation (MM/GBSA) free energy was comparable to known inhibitor Panobinostat. ZINC000002845205 qualifies drug-likeness according to Lipinski's rule-of-five, rule-of-three, and the World Drug Index (WDI)-like rule, but there is one violation in the lead-like rule.

Withasomniferol C, a new potential SARS-CoV-2 main protease inhibitor from the Withania somnifera plant proposed by in silico approaches

Exploring potent herbal medicine candidates is a promising strategy for combating a pandemic in the present global health crisis. In Ayurveda (a traditional medicine system in India), Withania somnifera (WS) is one of the most important herbs and it has been used for millennia as Rasayana (a type of juice) for its wide-ranging health benefits. WS phytocompounds display a broad spectrum of biological activities (such as antioxidant, anticancer and antimicrobial) modulate detoxifying enzymes, and enhance immunity. Inspired by the numerous biological actions of WS phytocompounds, the present investigation explored the potential of the WS phytocompounds against the SARS-CoV-2 main protease (3CL^pro). We selected 11 specific withanolide compounds, such as withaphysalin, withasomniferol, and withafastuosin, through manual literature curation against 3CL^pro. A molecular similarity analysis showed their similarity with compounds that have an established inhibitory activity against the SARS-CoV-2. In silico molecular docking and molecular dynamics simulations elucidated withasomniferol C (WS11) as a potential candidate against SARS-CoV-2 3CL^pro. Additionally, the present work also presents a new method of validating docking poses using the AlteQ method.

Quorum Sensing and Biofilm Disrupting Potential of Imidazole Derivatives in Chromobacterium violaceum Using Antimicrobial and Drug Discovery Approaches

Population of drug-resistant bacteria have increased at an alarming rate in the past few decades. The major reason for increasing drug resistance is the lack of new antibiotics and limited drug targets. It has therefore been a vital task to develop new antibiotics with different drug targets. Two such targets are biofilm formation and quorum sensing. Quorum sensing is cell to cell communication used by bacteria that initiates many important survival processes and aids in establishing pathogenesis. Both biofilm and quorum sensing are inter-related processes and play a major role in physiological and pathogenesis processes. In this study, five novel imidazole derivatives (IMA-1-IMA-5) were synthesised and tested for their antibacterial and anti-quorum sensing activities against Chromobacterium violaceum using different in silico and in vitro techniques following the standard protocols. In silico results revealed that all compounds were able to effectively bind to and interact sufficiently with the target protein CviR. CviR is a protein to which autoinducers bind to initiate the quorum sensing process. In silico results also revealed that the compounds generated favourable structural dynamics implying that the compounds would be able to effectively bind to CviR and inhibit quorum sensing. Susceptibility results revealed that IMA-1 is the most active of all the derivatives against both planktonic cells and biofilms. Qualitative and quantitative evaluation of anti-quorum sensing activity at sub-inhibitory concentrations of these compounds also revealed high activity for IMA-1. Down-regulation of most of the quorum sensing genes when cells were treated with the test compounds affirmed the high anti-quorum sensing activities of these compounds. The results from this study are promising and urges on the use of anti-quorum sensing and biofilm disrupting molecules to combat multi-drug resistance problem.

Identification of bio-active food compounds as potential SARS-CoV-2 PLpro inhibitors-modulators via negative image-based screening and computational simulations

Despite significant studies on the COVID-19 pandemic, scientists around the world are still battling to find a definitive therapy against the ongoing severe global health crisis. In this study, advanced computational approaches have been employed to identify bioactive food constituents as potential SARS-CoV-2 PLpro inhibitors-modulators. As a validated antiviral drug target, PLpro has gained tremendous attention for therapeutics developments. Therefore, targeting the multifunctional SARS-CoV-2 PLpro protein, ∼1039 bioactive dietary compounds have been screened extensively through novel techniques like negative image-based (NIB) screening and molecular docking approaches. In particular, the three different models of NIB screening have been generated and used to re-score the dietary compounds based on the negative image which is created by reversing the shape and electrostatics features of PLpro protein's ligand-binding cavity. Further, 100 ns molecular dynamics simulation has been performed and MM-GBSA based binding free energies have been estimated for the final proposed four dietary compounds (PC000550, PC000361, PC000558, and PC000573) as potential inhibitors/modulators of SARS-CoV-2 PLpro protein. Employed computational study outcome also has been compared with respect to the earlier experimentally investigated compound GRL0617 against SARS-CoV-2 PLpro protein, which suggests much greater interaction potential in terms of binding affinity and other energetic contributions for the proposed dietary compounds. Hence, the present study suggests that proposed dietary compounds can be suitable chemical entities for modulating the activity of PLpro protein or can be further utilized for optimizing or screening of novel chemical surrogates, however also needs experimental evaluation for entry in clinical studies for better assessment.

Synthetic Heparan Sulfate Mimetic Pixatimod (PG545) Potently Inhibits SARS-CoV-2 by Disrupting the Spike-ACE2 Interaction

Heparan sulfate (HS) is a cell surface polysaccharide recently identified as a coreceptor with the ACE2 protein for the S1 spike protein on SARS-CoV-2 virus, providing a tractable new therapeutic target. Clinically used heparins demonstrate an inhibitory activity but have an anticoagulant activity and are supply-limited, necessitating alternative solutions. Here, we show that synthetic HS mimetic pixatimod (PG545), a cancer drug candidate, binds and destabilizes the SARS-CoV-2 spike protein receptor binding domain and directly inhibits its binding to ACE2, consistent with molecular modeling identification of multiple molecular contacts and overlapping pixatimod and ACE2 binding sites. Assays with multiple clinical isolates of SARS-CoV-2 virus show that pixatimod potently inhibits the infection of monkey Vero E6 cells and physiologically relevant human bronchial epithelial cells at safe therapeutic concentrations. Pixatimod also retained broad potency against variants of concern (VOC) including B.1.1.7 (Alpha), B.1.351 (Beta), B.1.617.2 (Delta), and B.1.1.529 (Omicron). Furthermore, in a K18-hACE2 mouse model, pixatimod significantly reduced SARS-CoV-2 viral titers in the upper respiratory tract and virus-induced weight loss. This demonstration of potent anti-SARS-CoV-2 activity tolerant to emerging mutations establishes proof-of-concept for targeting the HS-Spike protein-ACE2 axis with synthetic HS mimetics and provides a strong rationale for clinical investigation of pixatimod as a potential multimodal therapeutic for COVID-19.

Identification of the probable structure of the sAPPα-GABABR1a complex and theoretical solutions for such cases

Amyloid precursor protein (APP) is the core of the pathogenesis of Alzheimer's disease (AD). Existing studies have shown that the soluble secreted APP (sAPPα) fragment obtained from the hydrolysis of APP by α-secretase has a synaptic function. Thereinto, a nine-residue fragment (APP9mer) of the extension domain region of sAPPα can bind directly and selectively to the N-terminal sushi1 domain (SD1) of the γ-aminobutyric acid type B receptor subunit 1a (GABA_BR1a) protein, which can influence synaptic transmission and plasticity by changing the GABA_BR1a conformation. APP9mer is a highly flexible, disordered region, and as such it is difficult to experimentally determine the optimal APPmer-SD1 binding complex. In this study we constructed two types of APP9mer-SD1 complexes through molecular docking and molecular dynamics simulation, aiming to explore the recognition function and mechanism of the specific binding of APP9mer with SD1, from which the most probable APPmer-SD1 model conformation is predicted. All the data from the analyses of RMSD, RMSF, PCA, DCCM and MM/PBSA binding energy as well as comparison with the experimental dissociation constant K_d suggest that 2NC is the most likely conformation to restore the crystal structure of the experimental APP9mer-SD1 complex. Of note, the key recognition residues of APP9mer are D24, D25, D27, W29 and W30, which mainly act on the 9-45 residue domain of SD1 (consisting of two loops and three short β-chains at the N-terminus of SD1). The mini-model with key residues identified establishes the molecular basis with deep insight into the interaction between APP and GABA_BR1a and provides a target for the development of therapeutic strategies for modulating GABA_BR1a-specific signaling in neurological and psychiatric disorders. More importantly, the study offers a theoretical solution for how to determine a biomolecular structure with a highly flexible, disordered fragment embedded within. The flexible fragment involved in a protein structure has to be deserted usually during the structural determination with experimental methods (e.g. X-ray crystallography, etc.).

Calcium Ion Binding to the Mutants of Calmodulin: A Structure-Based Computational Predictive Model of Binding Affinity Using a Charge Scaling Approach in Molecular Dynamics Simulation

The binding of calcium ions (Ca²⁺) to the calcium-binding proteins (CBPs) controls a plethora of regulatory processes. Among the roles played by CBPs in several diseases, the onset and progress of some cardiovascular diseases are caused by mutations in calmodulin (CaM), an important member of CBPs. Rationalization and prediction of the binding affinity of Ca²⁺ ions to the CaM can play important roles in understanding the origin of cardiovascular diseases. However, there is no robust structure-based computational method for predicting the binding affinity of Ca²⁺ ions to the different forms of CBPs in general and CaM in particular. In the current work, we have devised a fast yet accurate computational technique to accurately calculate the binding affinity of Ca²⁺ to the different forms of CaM. This method combines the well-known molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method and a charge scaling approach developed by previous authors that takes care of the polarization of CaM and Ca²⁺ ions. Our detailed analysis of the different components of binding free energy shows that subtle changes in electrostatics and van der Waals contribute to the difference in the binding affinity of mutants from that of the wild type (WT), and the charge scaling approach is superior in calculating these subtle changes in electrostatics as compared to the nonpolarizable force field used in this work. A statistically significant regression model made from our binding free energy calculations gives a correlation coefficient close to 0.8 to the experimental results. This structure-based predictive model can open up a new strategy to understand and predict the binding of Ca²⁺ to the mutants of CBPs, in general.

DNA Sequence Specificity Reveals a Role of the HLTF HIRAN Domain in the Recognition of Trinucleotide Repeats

DNA damage tolerance (DDT) pathways enable cells to cope with a variety of replication blocks that threaten their ability to complete DNA replication. Helicase-like transcription factor (HLTF) plays a central role in the error-free DDT pathway, template switching (TS), by serving as a ubiquitin ligase to polyubiquitinate the DNA sliding clamp PCNA, which promotes TS initiation. HLTF also serves as an ATP-dependent DNA translocase facilitating replication fork remodeling. The HIP116, Rad5p N-terminal (HIRAN) domain of HLTF specifically recognizes the unmodified 3'-end of single-stranded DNA (ssDNA) at stalled replication forks to promote fork regression. Several crystal structures of the HIRAN domain in complex with ssDNA have been reported; however, optimal ssDNA sequences for high-affinity binding with the domain have not been described. Here we elucidated DNA sequence preferences of HLTF HIRAN through systematic studies of its binding to ssDNA substrates using fluorescence polarization assays and a computational analysis of the ssDNA:HIRAN interaction. These studies reveal that the HLTF HIRAN domain preferentially recognizes a (T/C)TG sequence at the 3'-hydroxyl ssDNA end, which occurs in the CTG trinucleotide repeat (TNR) regions that are susceptible to expansion and deletion mutations identified in neuromuscular and neurodegenerative disorders. These findings support a role for HLTF in maintaining the stability of difficult to replicate TNR microsatellite regions.

Plant-derived active compounds as a potential nucleocapsid protein inhibitor of SARS-CoV-2: an in-silico study

The coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2. This virus has a high mismatch repair proofreading ability due to its unique exonuclease activity, making it knotty to treat. The nucleocapsid protein can serve as a potential antiviral drug target, as this protein is responsible for multiple captious functions during the viral life cycle. Herein, we have investigated the potential to repurpose active antiviral compounds of plant origins for treating the SARS-CoV-2 infection. In the present study, we followed the molecular docking methodology to screen druggable natural plants' active compounds against the nucleocapsid protein of SARS-CoV-2. The virtual screening of all 68 compounds revealed that the top seven active compounds, such as withanolide D, hypericin, silymarin, oxyacanthine, withaferin A, Acetyl aleuritolic acid, and rhein, exhibit good binding affinity with druggable ADME properties, toxicity, and Pass prediction. The stability of the docked complexes was studied by conducting molecular simulations of 100 ns. MM-GBSA calculated the binding free energy uncovered that withanolide D, hypericin, and silymarin result in highly stable binding conformations in three different sites of the nucleocapsid protein. However, further investigation is needed in order to validate the candidacy of these inhibitors for clinical trials. HighlightsNatural plants' active compounds may aid in the inhibition of SARS-CoV-2 replication and COVID-19 therapeutics.Hypericin, silymarin, withanolide D, oxyacanthine, withaferin A, Acetyl aleuritolic acid, and rhein are effective against SARS-CoV-2 N protein.Studied natural plants' active compounds could be useful against COVID-19 and its associated organs comorbidities.ADMET properties of selected compounds favor these compounds as druggable candidates.Communicated by Ramaswamy H. Sarma.

Investigation of molecular mechanism leading to gefitinib and osimertinib resistance against EGFR tyrosine kinase: molecular dynamics and binding free energy calculation

Tyrosine kinase (TK) is an important protein responsible for phosphorylation of variety of proteins that helps in signal transduction process in transferring signal to regulate various physiological and biochemical processes. Drugs inhibiting signal transduction pathways can be a very rational approach to inhibit cellular physiological and biochemical process. Tyrosine kinase inhibitors are a wide family of drugs that have been used successfully in cancer chemotherapy. Certain mutations around the catalytic cleft may cause conformational changes at binding site and leads to decrease in inhibitor sensitivity to TK mutants. EGFR^T790M mutation is the first recognized acquired resistance after tyrosine kinase inhibitor therapy that leads to resistant to first generation TKI in about 50% of non-small cell lung carcinoma patients. Third generation EGFR-TKIs bind irreversibly to the C797, which is present in the ATP-binding pocket. The present work provides a molecular mechanism for understanding the Gefitinib and Osimertinib sensitivities with the EGFR^WILD, EGFR^L858R, EGFR^T790M, EGFR^T790M+C797S mutants using molecular modelling techniques. Changes in response against Gefitinib and Osimertinib were observed with the change of amino acids at the tyrosine kinase domain of EGFR^WILD and its mutants (EGFR^L858R, EGFR^T790M, EGFR^T790M+C797S). RMSD, RMSF and binding energies calculation well correlates with the change in clinical observation.Communicated by Ramaswamy H. Sarma.

Structural insights into the discrepant synergistic activity of Codlemone and (Z)-8-dodecenol towards Grapholita molesta pheromones

BACKGROUND

Insect pheromone synergists have been widely used to produce potent pheromone products for environment-friendly pest control. Codlemone (Cod) and (Z)-8-dodecenol (Dod) are two major Grapholita molesta pheromone synergists, with Cod having greater synergism and affinity for G. molesta pheromone binding protein 2 (GmolPBP2). Uncovering structural information key to the different binding affinity of Cod and Dod to GmolPBP2 would gain insights into what causes their synergy activity discrepancy.

RESULTS

Binding modes of the two synergists in the binding pocket of GmolPBP2 were analyzed and compared by molecular dynamics-based approaches. Although Cod and Dod were stabilized in a similar hydrophobic pocket, their interaction details with GmolPBP2 were divergent due to the extra double bond (C10═C11) in Cod. The C10═C11 improved the hydrophobic interactions of Cod with around residues. Such hydrophobic interaction improvement was also reflected in the raised importance of Phe11 in the GmolPBP2-Cod interaction. Not only that, the increased hydrophobic forces introduced by the C10═C11 changed the CH2-OH orientation in the GmolPBP2-Cod complex, which improved the H-bond interaction. Electrostatic complementarity analysis further indicated the positive role of C10═C11 in optimizing GmolPBP2-Cod interaction.

CONCLUSION

The C10═C11 is thought to contribute greatly to Cod's stronger synergy as a group key to the higher GmolPBP2-affinity, based on which the improvement directions for Cod and Dod were addressed as well. Our findings will aid in the development and optimization of more effective pheromone synergists, resulting in more effective pheromone-based pest management.

Computational Insight into Biotransformation Profiles of Organophosphorus Flame Retardants to Their Diester Metabolites by Cytochrome P450

Biotransformation of organophosphorus flame retardants (OPFRs) mediated by cytochrome P450 enzymes (CYPs) has a potential correlation with their toxicological effects on humans. In this work, we employed five typical OPFRs including tris(1,3-dichloro-2-propyl) phosphate (TDCIPP), tris(1-chloro-2-propyl) phosphate (TCIPP), tri(2-chloroethyl) phosphate (TCEP), triethyl phosphate (TEP), and 2-ethylhexyl diphenyl phosphate (EHDPHP), and performed density functional theory (DFT) calculations to clarify the CYP-catalyzed biotransformation of five OPFRs to their diester metabolites. The DFT results show that the reaction mechanism consists of Cα-hydroxylation and O-dealkylation steps, and the biotransformation activities of five OPFRs may follow the order of TCEP ≈ TEP ≈ EHDPHP > TCIPP > TDCIPP. We further performed molecular dynamics (MD) simulations to unravel the binding interactions of five OPFRs in the CYP3A4 isoform. Binding mode analyses demonstrate that CYP3A4-mediated metabolism of TDCIPP, TCIPP, TCEP, and TEP can produce the diester metabolites, while EHDPHP metabolism may generate para-hydroxyEHDPHP as the primary metabolite. Moreover, the EHDPHP and TDCIPP have higher binding potential to CYP3A4 than TCIPP, TCEP, and TEP. This work reports the biotransformation profiles and binding features of five OPFRs in CYP, which can provide meaningful clues for the further studies of the metabolic fates of OPFRs and toxicological effects associated with the relevant metabolites.

Quantum Mechanics/Molecular Mechanics Studies on the Catalytic Mechanism of a Novel Esterase (FmtA) of Staphylococcus aureus

FmtA is a novel esterase that shares the penicillin-binding protein (PBP) core structural folding but found to hydrolyze the removal of d-Ala from teichoic acids. Molecular docking, dynamics, and MM-GBSA of FmtA and its variants S127A, K130A, Y211A, D213A, and K130AY211A, in the presence or absence of wall teichoic acid (WTA), suggest that active site residues S127, K130, Y211, D213, N343, and G344 play a role in substrate binding. Quantum mechanics (QM)/molecular mechanics (MM) calculations reveal that during WTA catalysis, K130 deprotonates S127, and the nucleophilic S127 attacks the carbonyl carbon of d-Ala bound to WTA. The tetrahedral intermediate (TI) complex is stabilized by hydrogen bonding to the oxyanion holes. The TI complex displays a high energy gap and collapses to an energetically favorable acyl-enzyme complex.

Application of amino acid ionic liquids for increasing the stability of DNA in long term storage

The structural stability of DNA is important because of its biological activity. DNAs due to their inherent chemical properties are not stable in an aqueous solution, therefore, a long period of storage of DNA at the ambient condition in bioscience is of importance. Ionic liquids (ILs) as interesting alternatives compared to organic solvents and water due to their considerable properties can be used as new agents to increase the stability of DNA for a long period of storage. In this article, molecular dynamics (MD) simulations and quantum chemistry calculations were applied to investigate the effects of amino acid ionic liquids ([BMIM][Ala], [BMIM][Gly], [BMIM][Val], [BMIM][Pro] and [BMIM][Leu]) on the dynamical behavior and the structural stability of calf thymus DNA. Based on the obtained MD results ILs enter into the solvation shell of the DNA and push away the water molecules from the DNA surface. Structural analysis shows that [BMIM]⁺ cations can occupy the DNA minor groove without disturbing the double-helical structure of DNA. ILs due to strong electrostatic and van der Waals (vdW) interactions with the DNA structure contribute to the stability of the double-helical structure. Quantum chemistry calculations indicate that the interactions between the [BMIM]⁺ cation and DNA structure has an electrostatic character. Moreover, this cation forms a more stable complex with the CGCG region of the DNA in comparison with AATT base pairs. Overall, the results of this study can provide new insight into the application of ILs for maintaining DNA stability during long-term storage.Communicated by Ramaswamy H. Sarma.

Origin of the tight binding mode to ACE2 triggered by multi-point mutations in the omicron variant: a dynamic insight

The continuous spread of the newly emerged SARS-CoV-2 Omicron variant (B.1.1.529) has become an important reason for the surge in COVID-19 infections. Its numerous mutated residues containing key sites on the receptor-binding domain (RBD) undoubtedly pose new challenges for epidemic control. Although the preventive measures are becoming more sophisticated, the effects of mutations on the binding of the virus to the receptor protein remain to be elucidated. Here, we used molecular dynamics (MD) simulations to investigate the differences in the binding mode between the Omicron variant and the angiotensin-converting enzyme 2 (ACE2) compared to the wild-type strain (WT). Multi-point mutations in the Omicron variant RBD could cause the conformation shift in the large Loop (where T478K and E484A are located), which makes it easier to wrap the N-terminal helix of ACE2 and form tighter contacts. The stronger electrostatic interaction was the main reason for its enhanced binding affinity as compared to WT. This was due to the large number of positively charged patches (N440K, T478K, Q493R, Q498R, and Y505H) formed by the substitution of neutral amino acids at multiple sites. The appearance of these highly polar hydrophilic amino acids may cause local perturbations and affect the electrostatic complementarity of RBD with the ACE2, and further mediate conformational changes. Thus, a more extensive interaction network was found in the mutation system and the complex interaction cluster was formed near E37@ACE2, which was essential for the stable binding of the two. In addition, we speculated that these mutations may affect the electrostatic complementarity with the four potential antibodies to reduce the sensitivity of the virus to antibodies. This study reveals the key details of the Omicron variant binding to ACE2 and provides important theoretical views for the enhanced infectivity of this variant. We hope that these observations can provide timely molecular insights for responding to the Omicron variant pandemic.

Losartan Interactions with 2-Hydroxypropyl-β-CD

Losartan potassium salt (LSR) is a well-known antihypertensive drug with proven beneficial effects on human health. Its formulation with the non-toxic 2-hydroxypropyl-β-cyclodextrin (2-HP-β-CD) could improve its pharmacological profile. Thus, its molecular interactions are studied using a combination of Differential Scanning Calorimetry (DSC), Nuclear Magnetic Resonance (NMR) and Molecular Dynamics (MD). First, its complexation is shown through Differential Scanning Calorimetry as lyophilization provided distinct thermal properties in comparison to the mixture. The complexation is further proved by utilizing the chemical shift changes in the complexation and T1 values. Furthermore, the reversible favorable complexation was shown by MD calculations. Such physical chemical properties provide evidence that this formulation must be further explored through biological experiments.

Identification of Trovafloxacin, Ozanimod, and Ozenoxacin as potent c-Myc G-quadruplex stabilizers to suppress c-Myc transcription and myeloma growth

c-Myc is a major oncogene that is estimated to result in almost all human cancers and the c-Myc downregulation has become an attractive strategy for cancer treatment. For it is hard to design compounds that can directly interact with the c-Myc protein, the DNA G-quadruplex (G4) was discovered in its promoter region which was referred to as a potential drug target for controlling c-Myc expression. In this study, a combined strategy of molecular docking-based virtual screening, molecular dynamics (MD) simulation, and molecular mechanics/generalized Born surface area (MM/GBSA) free energy calculation was conducted on the existing FDA-Approved Drugs Library, eight compounds were selected for further experimental assay. Among them, five compounds exhibited dose-dependently anticancer activities against RPMI-8226 cells with IC50 values less than 18.4 μM. Further experiments showed that Trovafloxacin, Ozanimod, and Ozenoxacin decreased c-Myc mRNA level obviously and downregulated c-Myc expression significantly. In summary, compounds Trovafloxacin, Ozanimod, and Ozenoxacin might be regarded as new c-Myc G4 stabilizers for the treatment of c-Myc related cancers in the future.

sVmKTx, a transcriptome analysis-based synthetic peptide analogue of Vm24, inhibits Kv1.3 channels of human T cells with improved selectivity

Kv1.3 K⁺ channels play a central role in the regulation of T cell activation and Ca²⁺ signaling under physiological and pathophysiological conditions. Peptide toxins targeting Kv1.3 have a significant therapeutic potential in the treatment of autoimmune diseases; thus, the discovery of new toxins is highly motivated. Based on the transcriptome analysis of the venom gland of V. mexicanus smithi a novel synthetic peptide, sVmKTx was generated, containing 36 amino acid residues. sVmKTx shows high sequence similarity to Vm24, a previously characterized peptide from the same species, but contains a Glu at position 32 as opposed to Lys32 in Vm24. Vm24 inhibits Kv1.3 with high affinity (K_d = 2.9 pM). However, it has limited selectivity (~1,500-fold) for Kv1.3 over hKv1.2, hKCa3.1, and mKv1.1. sVmKTx displays reduced Kv1.3 affinity (K_d = 770 pM) but increased selectivity for Kv1.3 over hKv1.2 (~9,000-fold) as compared to Vm24, other channels tested in the panel (hKCa3.1, hKv1.1, hKv1.4, hKv1.5, rKv2.1, hKv11.1, hKCa1.1, hNav1.5) were practically insensitive to the toxin at 2.5 μM. Molecular dynamics simulations showed that introduction of a Glu instead of Lys at position 32 led to a decreased structural fluctuation of the N-terminal segment of sVmKTx, which may explain its increased selectivity for Kv1.3. sVmKTx at 100 nM concentration decreased the expression level of the Ca²⁺ -dependent T cell activation marker, CD40 ligand. The high affinity block of Kv1.3 and increased selectivity over the natural peptide makes sVmKTx a potential candidate for Kv1.3 blockade-mediated treatment of autoimmune diseases.

Molecular Dynamics-Derived Descriptor Informs the Impact of Mutation on the Catalytic Turnover Number in Lactonase Across Substrates

Molecular dynamics simulations have been extensively employed to reveal the roles of protein dynamics in mediating enzyme catalysis. However, simulation-derived predictive descriptors that inform the impacts of mutations on catalytic turnover numbers remain largely unexplored. In this work, we report the identification of molecular modeling-derived descriptors to predict mutation effect on the turnover number of lactonase SsoPox with both native and non-native substrates. The study consists of 10 enzyme-substrate complexes resulting from a combination of five enzyme variants with two substrates. For each complex, we derived 15 descriptors from molecular dynamics simulations and applied principal component analysis to rank the predictive capability of the descriptors. A top-ranked descriptor was identified, which is the solvent-accessible surface area (SASA) ratio of the substrate to the active site pocket. A uniform volcano-shaped plot was observed in the distribution of experimental activation free energy against the SASA ratio. To achieve efficient lactonase hydrolysis, a non-native substrate-bound enzyme variant needs to involve a similar range of the SASA ratio to the native substrate-bound wild-type enzyme. The descriptor reflects how well the enzyme active site pocket accommodates a substrate for reaction, which has the potential of guiding optimization of enzyme reaction turnover for non-native chemical transformations.

Folate Conjugated Polyethylene Glycol Probe for Tumor-Targeted Drug Delivery of 5-Fluorouracil

A targeted delivery system is primarily intended to carry a potent anticancer drug to specific tumor sites within the bodily tissues. In the present study, a carrier system has been designed using folic acid (FA), bis-amine polyethylene glycol (PEG), and an anticancer drug, 5-fluorouracil (5-FU). FA and PEG were joined via an amide bond, and the resulting FA-PEG-NH₂ was coupled to 5-FU producing folate-polyethylene glycol conjugated 5-fluorouracil (FA-PEG-5-FU). Spectroscopic techniques (UV-Vis, ¹HNMR, FTIR, and HPLC) were used for the characterization of products. Prodrug (FA-PEG-5-FU) was analyzed for drug release profile (in vitro) up to 10 days and compared to a standard anticancer drug (5-FU). Folate conjugate was also analyzed to study its folate receptors (FR) mediated transport and in vitro cytotoxicity assays using HeLa cancer cells/Vero cells, respectively, and antitumor activity in tumor-bearing mice models. Folate conjugate showed steady drug release patterns and improved uptake in the HeLa cancer cells than Vero cells. Folate conjugate treated mice group showed smaller tumor volumes; specifically after the 15th day post-treatment, tumor sizes were decreased significantly compared to the standard drug group (5-FU). Molecular docking findings demonstrated importance of Trp138, Trp140, and Lys136 in the stabilization of flexible loop flanking the active site. The folic acid conjugated probe has shown the potential of targeted drug delivery and sustained release of anticancer drug to tumor lesions with intact antitumor efficacy.

Reorienting Mechanism of Harderoheme in Coproheme Decarboxylase-A Computational Study

Coproheme decarboxylase (ChdC) is an important enzyme in the coproporphyrin-dependent pathway (CPD) of Gram-positive bacteria that decarboxylates coproheme on two propionates at position 2 and position 4 sequentially to generate heme b by using H₂O₂ as an oxidant. This work focused on the ChdC from Geobacillus stearothermophilus (GsChdC) to elucidate the mechanism of its sequential two-step decarboxylation of coproheme. The models of GsChdC in a complex with substrate and reaction intermediate were built to investigate the reorienting mechanism of harderoheme. Targeted molecular dynamics simulations on these models validated that harderoheme is able to rotate in the active site of GsChdC with a 19.06-kcal·mol^-1 energy barrier after the first step of decarboxylation to bring the propionate at position 4 in proximity of Tyr145 to continue the second decarboxylation step. The harderoheme rotation mechanism is confirmed to be much easier than the release-rebinding mechanism. In the active site of GsChdC, Trp157 and Trp198 comprise a "gate" construction to regulate the clockwise rotation of the harderoheme. Lys149 plays a critical role in the rotation mechanism, which not only keeps the Trp157-Trp198 "gate" from being closed but also guides the propionate at position 4 through the gap between Trp157 and Trp198 through a salt bridge interaction.

Thermodynamic and structural insights into the repurposing of drugs that bind to SARS-CoV-2 main protease

Although researchers have been working tirelessly since the COVID-19 outbreak, so far only three drugs - remdesivir, ronapreve and molnupiravir - have been approved for use in some countries which directly target the SARS-CoV-2 virus. Given the slow pace and substantial costs of new drug discovery and development, together with the urgency of the matter, repurposing of existing drugs for the ongoing disease is an attractive proposition. In a recent study, a high-throughput X-ray crystallographic screen was performed for a selection of drugs which have been approved or are in clinical trials. Thirty-seven compounds have been identified from drug libraries all of which bind to the SARS-CoV-2 main protease (3CL^pro). In the current study, we use molecular dynamics simulation and an ensemble-based free energy approach, namely, enhanced sampling of molecular dynamics with approximation of continuum solvent (ESMACS), to investigate a subset of the aforementioned compounds. The drugs studied here are highly diverse, interacting with different binding sites and/or subsites of 3CL^pro. The predicted free energies are compared with experimental results wherever they are available and they are found to be in excellent agreement. Our study also provides detailed energetic insights into the nature of the associated drug-protein binding, in turn shedding light on the design and discovery of potential drugs.

Discovery of potential competitive inhibitors against With-No-Lysine kinase 1 for treating hypertension by virtual screening, inverse pharmacophore-based lead optimization, and molecular dynamics simulations

The With-No-Lysine (WNK) has received attention because of its involvement in hypertension. Genetic mutation in the genes of WNK, leading to its overexpression, has been reported in Familial Hyperkalaemic Hypertension, and thus WNK is considered a potential drug target. Herein, we have performed a high-throughput virtual screening of ~11,000 compounds, mainly the natural phytochemical compounds and kinase inhibitory libraries, to find potential competitive inhibitors against WNK1. Initially, candidates with a docking score of ~ -10.0 kcal/mol or less were selected to further screen their good pharmacological properties by applying absorption, distribution, metabolism, excretion, and toxicity (ADMET). Finally, six docked compounds bearing appreciable binding affinities and WNK1 selectivity were complimented with 500 ns long all-atom molecular dynamic simulations. Subsequently, the MMPBSA scheme (Molecular Mechanics Poisson Boltzmann Surface Area) suggested three phytochemical compounds, C00000947, C00020451, and C00005031, with favourable binding affinity against WNK1. Among them, C00000947 acts as the most potent competitive inhibitor of WNK1. Further, inverse pharmacophore-based lead optimization of the C00000947 leads to one potential compound, meciadanol, which shows better binding affinity and specificity than C00000947 towards WNK1, which may be further exploited to develop effective therapeutics against WNK1-associated hypertension after in vitro and in vivo validation.

Non-synonymous substitution of evolutionarily conserved residue in Tau class glutathione transferases alters structural and catalytic features

Plant-specific tau glutathione transferases (GSTs) are basically involved in catalysing γ-glutathione (GSH)-dependent conjugation reactions with pesticides and herbicides, which play an important role in the detoxification of pollutants. Given the lack of systematic biochemical and structural information on tau GSTs, the study of their mediated defence mechanisms against toxic compounds has been greatly hindered. Here, we reveal the importance of the Ile residue closely interacting with GSH for the structural stability and catalytic function of GST. Evolutionary conservation analysis indicated that the crucial G-site Ile55 in the SbGSTU6 was converted to Thr53 of SbGSTU7. The comparative biochemical data on SbGSTU6, SbGSTU7 and their mutants showed that the substitution of Ile by Thr caused significant decrease in the affinity and catalytic efficiency of the GSTs. The unfavourable structural flexibility and pKa distribution of the active cavity residues were also demonstrated. Crystallography studies and molecular dynamics simulations showed that the conversion resulted in the hydrogen bond recombination with GSH and conformational rearrangement of GST active cavity, in which the Ile residue was more conducive to the formation of enzyme substrate complexes. The extensive biochemical and structural data not only reveal the critical role of the conserved G-site Ile residue in catalysing GSH-conjugate reactions but also provide valuable resources for the development of GST engineering in analytical and agricultural biotechnology.

Molecular Docking and Dynamics Studies to Explore Effective Inhibitory Peptides Against the Spike Receptor Binding Domain of SARS-CoV-2

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has become a pandemic due to the high transmission and mortality rate of this virus. The world health and economic sectors have been severely affected by this deadly virus, exacerbated by the lack of sufficient efficient vaccines. The design of effective drug candidates and their rapid development is necessary to combat this virus. In this study, we selected 23 antimicrobial peptides from the literature and predicted their structure using PEP-FOLD 3.5. In addition, we docked them to the SARS-CoV-2 spike protein receptor-binding domain (RBD) to study their capability to inhibit the RBD, which plays a significant role in virus binding, fusion and entry into the host cell. We used several docking programs including HDOCK, HPEPDOCK, ClusPro, and HawkDock to calculate the binding energy of the protein-peptide complexes. We identified four peptides with high binding free energy and docking scores. The docking results were further verified by molecular dynamics (MD) simulations to characterize the protein-peptide complexes in terms of their root-mean-square fluctuation (RMSF), root-mean-square deviation (RMSD), radius of gyration (Rg), solvent-accessible surface area (SASA), and hydrogen bond formation. Allergenicity and toxicity predictions suggested that the peptides we identified were non-allergenic and non-toxic. This study suggests that these four antimicrobial peptides could inhibit the RBD of SARS-CoV-2. Future in vitro and in vivo studies are necessary to confirm this.

A Computer-Aided Approach for the Discovery of D-Peptides as Inhibitors of SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease, also known as 3-chymotrypsin-like protease (3CLpro), is a cysteine protease responsible for the cleavage of viral polyproteins pp1a and pp1ab, at least, at eleven conserved sites, which leads to the formation of mature nonstructural proteins essential for the replication of the virus. Due to its essential role, numerous studies have been conducted so far, which have confirmed 3CLpro as an attractive drug target to combat Covid-19 and have reported a vast number of inhibitors and their co-crystal structures. Despite all the ongoing efforts, D-peptides, which possess key advantages over L-peptides as therapeutic agents, have not been explored as potential drug candidates against 3CLpro. The current work fills this gap by reporting an in silico approach for the discovery of D-peptides capable of inhibiting 3CLpro that involves structure-based virtual screening (SBVS) of an in-house library of D-tripeptides and D-tetrapeptides into the protease active site and subsequent rescoring steps, including Molecular Mechanics Generalized-Born Surface Area (MM-GBSA) free energy calculations and molecular dynamics (MD) simulations. In vitro enzymatic assays conducted for the four top-scoring D-tetrapeptides at 20 μM showed that all of them caused 55–85% inhibition of 3CLpro activity, thus highlighting the suitability of the devised approach. Overall, our results present a promising computational strategy to identify D-peptides capable of inhibiting 3CLpro, with broader application in problems involving protein inhibition.

In Silico Identification and Experimental Validation of Peptide-Based Inhibitors Targeting Clostridium difficile Toxin A

Clostridium difficile infection is mediated by two major exotoxins: toxins A (TcdA) and B (TcdB). Inhibiting the biocatalytic activities of these toxins with targeted peptide-based drugs can reduce the risk of C. difficile infection. In this work, we used a computational strategy that integrates a peptide binding design (PepBD) algorithm and explicit-solvent atomistic molecular dynamics simulation to determine promising toxin A-targeting peptides that can recognize and bind to the catalytic site of the TcdA glucosyltransferase domain (GTD). Our simulation results revealed that two out of three in silico discovered peptides, viz. the neutralizing peptides A (NPA) and B (NPB), exhibit lower binding free energies when bound to the TcdA GTD than the phage-display discovered peptide, viz. the reference peptide (RP). These peptides may serve as potential inhibitors against C. difficile infection. The efficacy of the peptides RP, NPA, and NPB to neutralize the cytopathic effects of TcdA was tested in vitro in human jejunum cells. Both phage-display peptide RP and in silico peptide NPA were found to exhibit strong toxin-neutralizing properties, thereby preventing the TcdA toxicity. However, the in silico peptide NPB demonstrates a relatively low efficacy against TcdA.

Decoding the Host-Parasite Protein Interactions Involved in Cerebral Malaria Through Glares of Molecular Dynamics Simulations

Malaria causes millions of deaths every year. The malaria parasite spends a substantial part of its life cycle inside human erythrocytes. Inside erythrocytes, it synthesizes and displays various proteins onto the erythrocyte surface, such as Plasmodium falciparum erythrocytic membrane protein-1 (PfEMP1). This protein contains cysteine-rich interdomain region (CIDR) domains which have many subtypes based on sequence diversity and can cross-talk with host molecules. The CIDRα1.4 subtype can attach host endothelial protein C receptor (EPCR). This interaction facilitates infected erythrocyte adherence to brain endothelium and subsequent development of cerebral malaria. Through molecular dynamics simulations in conjunction with the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method, we explored the mechanism of interaction in the CIDRα1-EPCR complex. We examined the structural behavior of two CIDRα1 molecules (encoded by HB3-isolate var03-gene and IT4-isolate var07-gene) with EPCR unbound and bound (complex) forms. ^HB3var03CIDRα1 in apo and complexed with EPCR was comparatively more stable than ^IT4var07CIDRα1. Both of the complexes adopted two distinct conformational energy states. The hydrophobic residues played a crucial role in the binding of both complexes. For ^HB3var03CIDRα1-EPCR, the dominant energetic components were total polar interactions, while in ^IT4var07CIDRα1-EPCR, the primary interaction was van der Waals and nonpolar solvation energy. The study also revealed details such as correlated conformational motions and secondary structure evolution. Further, it elucidated various hotspot residues involved in protein-protein recognition. Overall, our study provides additional information on the structural behavior of CIDR molecules in unbound and receptor-bound states, which will help to design potent inhibitors.