A novel approach of dynamic cross correlation analysis on molecular dynamics simulations and its application to Ets1 dimer-DNA complex

Assessment of the in vitro anti-diabetic activity with molecular dynamic simulations of limonoids isolated from Adalia lemon peels

Limonoids are important constituents of citrus that have a significant impact on promoting human health. Therefore, the primary focus of this research was to assess the overall limonoid content and isolate limonoids from Adalia lemon (Citrus limon L.) peels for their potential use as antioxidants and anti-diabetic agents. The levels of limonoid aglycones in the C. limon peel extract were quantified through a colorimetric assay, revealing a concentration of 16.53 ± 0.93 mg/L limonin equivalent. Furthermore, the total concentration of limonoid glucosides was determined to be 54.38 ± 1.02 mg/L. The study successfully identified five isolated limonoids, namely limonin, deacetylnomilin, nomilin, obacunone 17-O-β-D-glucopyranoside, and limonin 17-O-β-D-glucopyranoside, along with their respective yields. The efficacy of the limonoids-rich extract and the five isolated compounds was evaluated at three different concentrations (50, 100, and 200 µg/mL). It was found that both obacunone 17-O-β-D-glucopyranoside and limonin 17-O-β-D-glucopyranoside possessed the highest antioxidant, free radical scavenging, and anti-diabetic activities, followed by deacetylnomilin, and then the limonoids-rich extract. The molecular dynamic simulations were conducted to predict the behavior of the isolated compounds upon binding to the protein's active site, as well as their interaction and stability. The results revealed that limonin 17-O-β-D-glucopyranoside bound to the protein complex system exhibited a relatively more stable conformation than the Apo system. The analysis of Solvent Accessible Surface Area (SASA), in conjunction with the data obtained from Root-Mean-Square Deviation (RMSD), Root-Mean-Square Fluctuation (RMSF), and Radius of Gyration (ROG) computations, provided further evidence that the limonin 17-O-β-D-glucopyranoside complex system remained stable within the catalytic domain binding site of the human pancreatic alpha-amylase (HPA)-receptor. The research findings suggest that the limonoids found in Adalia lemon peels have the potential to be used as effective natural substances in creating innovative therapeutic treatments for conditions related to oxidative stress and disorders in carbohydrate metabolism.

In silico analysis of potential inhibitors for breast cancer targeting 17beta-hydroxysteroid dehydrogenase type 1 (17beta-HSD1) catalyses

Breast cancer (BC) is still one of the major issues in world health, especially for women, which necessitates innovative therapeutic strategies. In this study, we investigated the efficacy of retinoic acid derivatives as inhibitors of 17beta-hydroxysteroid dehydrogenase type 1 (17beta-HSD1), which plays a crucial role in the biosynthesis and metabolism of oestrogen and thereby influences the progression of BC and, the main objective of this investigation is to identify the possible drug candidate against BC through computational drug design approach including PASS prediction, molecular docking, ADMET profiling, molecular dynamics simulations (MD) and density functional theory (DFT) calculations. The result has reported that total eight derivatives with high binding affinity and promising pharmacokinetic properties among 115 derivatives. In particular, ligands 04 and 07 exhibited a higher binding affinity with values of -9.9 kcal/mol and -9.1 kcal/mol, respectively, than the standard drug epirubicin hydrochloride, which had a binding affinity of -8.2 kcal/mol. The stability of the ligand-protein complexes was further confirmed by MD simulations over a 100-ns trajectory, which included assessments of hydrogen bonds, root mean square deviation (RMSD), root mean square Fluctuation (RMSF), dynamic cross-correlation matric (DCCM) and principal component analysis. The study emphasizes the need for experimental validation to confirm the therapeutic utility of these compounds. This study enhances the computational search for new BC drugs and establishes a solid foundation for subsequent experimental and clinical research.

Facile and Novel Synthetic Approach, Molecular Docking, Molecular Dynamics, and Drug-Likeness Evaluation of 9-Substituted Acridine Derivatives as Dual Anticancer and Antimicrobial Agents

In the present study, numerous acridine derivatives A1-A20 were synthesized via aromatic nucleophilic substitution (SNAr) reaction of 9-chloroacridine with carbonyl hydrazides, amines, or phenolic derivatives depending upon facile, novel, and eco-friendly approaches (Microwave and ultrasonication assisted synthesis). The structures of the new compounds were elucidated using spectroscopic methods. The title products were assessed for their antimicrobial, antioxidant, and antiproliferative activities using numerous assays. Promisingly, the investigated compounds mainstream revealed promising antibacterial and anticancer activities. Thereafter, the investigated compounds' expected mode of action was debated by using an array of in silico studies. Compounds A2 and A3 were the most promising antimicrobial agents, while compounds A2, A5, and A7 revealed the most cytotoxic activities. Accordingly, RMSD, RMSF, Rg, and SASA analyses of compounds A2 and A3 were performed, and MMPBSA was calculated. Lastly, the ADMET (absorption, distribution, metabolism, excretion, and toxicity) analyses of the novel acridine derivatives were investigated. The tested compounds' existing screening results afford an inspiring basis leading to developing new compelling antimicrobial and anticancer agents based on the acridine scaffold.

A computational strategy for systematic virtual screening of plasmodium falciparum heme detoxification protein inhibitors from the Drugbank database

Antimalarial drug resistance poses one of the greatest threats to malaria treatment, resulting in increased morbidity and mortality. Heme Detoxification Protein (HDP) is among the essential hemoglobinases of P. falciparum (Pf), a vital molecular target for the treatment of malaria. In this study, we utilized the virtual screening workflow tool of the Schrodinger suite to find the best hits for the PfHDP from the DrugBank library. A total of 14,942 compounds were identified against the PfHDP. The top compounds with the highest docking scores and least energy scores were subjected to molecular simulations for 500 nanosecond to check the stability of the protein-drug complexes. The top three DrugBank compounds were found to be stable over 500 ns, namely DB09298 (silibinin), DB07426 (1-Hydroxy-2-(1,1':3',1''-Terphenyl-3-Yloxy) Ethane-1,1-Diyl] Bis (Phosphonic Acid), and DB07410 [(2-(3-Dibenzofuran-4-yl-Phenyl)-1-Hydroxy-1-Phosphono-Ethyl]-Phosphonic Acid). Overall analysis suggests that the top three compounds, DB09298, DB07426, and DB07410, have good stability for 500 ns. Their scaffolds can be used to design and develop new analogs of the target HDP protein. Silibinin, the anti-cancer drug, was found to be highly stable for the entire simulation period as compared to the other compounds. DB07426 shows its therapeutic effect on bones, especially in the treatment of osteoporosis, and DB07410 has anti-tumor, antibacterial, anti-oxidative, and anti-viral activities. All three compounds can be considered for repurposing as antimalarial drugs to evaluate the binding capacity or inhibition potential of these compounds. Further in-vivo and in-vitro analysis against the PfHDP protein should be conducted.Communicated by Ramaswamy H. Sarma.

Exploring the mechanism of the PTP1B inhibitors by molecular dynamics and experimental study

Protein tyrosine phosphatase 1B (PTP1B) has proven to be an attractive target for the treatment of cancer, diabetes and other diseases. Although many PTP1B inhibitors with various scaffolds have been developed, there is still a lack of PTP1B inhibitor with high specificity and acceptable pharmacological properties. Therefore, it is urgent to develop more methods to explore complex action mode of PTP1B and ligands for designing ideal PTP1B modulators. In this work, we developed a potential molecular dynamics (MD) analytic mode to analyze the mechanism of active compounds 6a and 6e against PTP1B from different perspectives, including the stable ability, interactions and binding site of ligand and protein, the binding energy, relative movement between residues and changes in protein internal interactions. The simulated results demonstrated that compound 6a bound more stably to the active pocket of PTP1B than 6e due to its smaller molecular volume (326 Å3), matched electronegativity, and enhanced the positive correlation motion of residues, especially for WPD loop and P loop. Lastly, compound 6a as a competitive inhibitor for PTP1B was verified by enzyme kinetic assay. This work successfully studied the mechanism of compound 6a against PTP1B from various aspects, enriched the analysis of interaction mode between PTP1B and inhibitors. In summary, we hope that this work could provide more theoretical information for designing and developing more novel and ideal PTP1B inhibitors in the future.

Structure-based virtual screening and molecular docking approaches to identify potential inhibitors against KIF2C to combat glioma

Glioma, a kind of malignant brain tumor, is extremely lethal. Kinesin family member 2C (KIF2C) was found to have an aberrant expression in several cancer types, including lung cancer and glioma. KIF2C may therefore be a useful therapeutic target for the treatment of glioma. In the current study, new drug candidates that may function as KIF2C enzyme inhibitors were discovered. MTi OpenScreen was used to carry out the structure-based virtual screening of an inbuilt drug library containing 150,000 compounds. These compounds belong to different classes, such as natural product-based compounds (NP-lib), purchasable approved drugs (Drugs-lib), and food constituents compound collection (FOOD-lib). Based on their binding affinities, a total of 84 compounds were further pushed to calculate ADMET properties. The compounds (16) meeting the ADMET cutoff ranges were then further docked to the receptor to find their plausible binding modes using the Glide tool's standard precision (SP) technique. The docking results were examined using the Glide gscore, and the best binding compounds (Rimacalib and Sarizotan) were chosen to test their stability with KIF2C protein through molecular dynamics (MD) simulation. Similarly, Principal Component Analysis and cross-correlation matrix were also examined. The MM/GBSA binding free energies showed a considerable energy contribution in the binding of hits with the KIF2C. Collectively, these findings strongly suggest the potential of the lead compounds to inhibit the biological function of KIF2C, emphasizing the need for further investigation in this area.Communicated by Ramaswamy H. Sarma.

Free energy landscape and thermodynamics properties of novel mutations in PncA of pyrazinamide resistance isolates of Mycobacterium tuberculosis

Pyrazinamide (PZA) is one of the first-line antituberculosis therapy, active against non-replicating Mycobacterium tuberculosis (Mtb). The conversion of PZA into pyrazinoic acid (POA), the active form, required the activity of pncA gene product pyrazinamidase (PZase) activity. Mutations occurred in pncA are the primary cause behind the PZA resistance. However, the resistance mechanism is important to explore using high throughput computational approaches. Here we aimed to explore the mechanism of PZA resistance behind novel P62T, L120R, and V130M mutations in PZase using 200 ns molecular dynamics (MD) simulations. MD simulations were performed to observe the structural changes for these three mutants (MTs) compared to the wild types (WT). Root means square fluctuation, the radius of gyration, free energy landscape, root means square deviation, dynamic cross-correlation motion, and pocket volume were found in variation between WT and MTs, revealing the effects of P62T, L120R, and V130M. The free energy conformational landscape of MTs differs significantly from the WT system, lowering the binding of PZA. The geometric shape complementarity of the drug (PZA) and target protein (PZase) further confirmed that P62T, L120R, and V130M affect the protein structure. These effects on PZase may cause vulnerability to convert PZA into POA.Communicated by Ramaswamy H. Sarma.

How do the mutations in PfK13 protein promote anti-malarial drug resistance?

Plasmodium falciparum develops resistance to artemisinin upon exposure to the anti-malarial drug. Various mutations in the Plasmodium falciparum Kelch13 (PfK13) protein such as Y493H, R539T, I543T and C580Y have been associated with anti-malarial drug resistance. These mutations impede the regular ubiquitination process that eventually invokes drug resistance. However, the relationship between the mutation and the mechanism of drug resistance has not yet been fully elucidated. The comparative protein dynamics are studied by performing the classical molecular dynamics (MD) simulations and subsequent analysis of the trajectories adopting root-mean-square fluctuations, the secondary-structure predictions and the dynamical cross-correlation matrix analysis tools. Here, we observed that the mutations in the Kelch-domain do not have any structural impact on the mutated site; however, it significantly alters the overall dynamics of the protein. The loop-region of the BTB-domain especially for Y493H and C580Y mutants is found to have the enhanced dynamical fluctuations. The enhanced fluctuations in the BTB-domain could affect the protein-protein (PfK13-Cullin) binding interactions in the ubiquitination process and eventually lead to anti-malarial drug resistance.Communicated by Ramaswamy H. Sarma.

Identification of a new Kir6 potassium channel and comparison of properties of Kir6 subtypes by structural modelling and molecular dynamics

ATP-sensitive potassium ion channels (KATP) are transmembrane proteins that modulate insulin release and muscle contraction. KATP channels are composed of two types of subunit, Kir6 and SUR, which exist in two and three isoforms respectively with different tissue distribution. In this work, we identify a previously undescribed ancestral vertebrate gene encoding a Kir6-related protein that we have named Kir6.3, which may not have a SUR binding partner, unlike the other two Kir6 proteins. Whereas Kir6.3 was lost in amniotes including mammals, it is still present in several early-diverging vertebrate lineages such as frogs, coelacanth, and rayfinned fishes. Molecular dynamics (MD) simulations using homology models of Kir6.1, Kir6.2, and Kir6.3 from the coelacanth Latimeria chalumnae showed that the three proteins exhibit subtle differences in their dynamics. Steered MD simulations of Kir6-SUR pairs suggest that Kir6.3 has a lower binding affinity for the SUR proteins than either Kir6.1 or Kir6.2. As we found no additional SUR gene in the genomes of the species that have Kir6.3, it most likely forms a lone tetramer. These findings invite studies of the tissue distribution of Kir6.3 in relation to the other Kir6 as well as SUR proteins to determine the functional roles of Kir6.3.

Molecular insights into the inhibition mechanism of harringtonine against essential proteins associated with SARS-CoV-2 entry

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has recently posed a serious threat to global public health. Harringtonine (HT), as a small-molecule antagonist, has antiviral activity against a variety of viruses. There is evidence that HT can inhibit the SARS-CoV-2 entry into host cells by blocking the Spike protein and transmembrane protease serine 2 (TMPRSS2). However, the molecular mechanism underlying the inhibition effect of HT is largely elusive. Here, docking and all-atom molecular dynamics simulations were used to investigate the mechanism of HT against the receptor binding domain (RBD) of Spike, TMPRSS2, as well as the complex of RBD and angiotensin-converting enzyme 2 complex (RBD-ACE2). The results reveal that HT binds to all proteins primarily through hydrogen bond and hydrophobic interactions. Binding with HT influences the structural stability and dynamic motility processes of each protein. The interactions of HT with residues N33, H34 and K353 of ACE2, and residue K417 and Y453 of RBD contribute to disrupting the binding affinity between RBD and ACE2, which may hinder the virus entry into host cells. Our research provides molecular insights into the inhibition mechanism of HT against SARS-CoV-2 associated proteins, which will help for the novel antiviral drugs development.

MDiGest: A Python package for describing allostery from molecular dynamics simulations

Many biological processes are regulated by allosteric mechanisms that communicate with distant sites in the protein responsible for functionality. The binding of a small molecule at an allosteric site typically induces conformational changes that propagate through the protein along allosteric pathways regulating enzymatic activity. Elucidating those communication pathways from allosteric sites to orthosteric sites is, therefore, essential to gain insights into biochemical processes. Targeting the allosteric pathways by mutagenesis can allow the engineering of proteins with desired functions. Furthermore, binding small molecule modulators along the allosteric pathways is a viable approach to target reactions using allosteric inhibitors/activators with temporal and spatial selectivity. Methods based on network theory can elucidate protein communication networks through the analysis of pairwise correlations observed in molecular dynamics (MD) simulations using molecular descriptors that serve as proxies for allosteric information. Typically, single atomic descriptors such as α-carbon displacements are used as proxies for allosteric information. Therefore, allosteric networks are based on correlations revealed by that descriptor. Here, we introduce a Python software package that provides a comprehensive toolkit for studying allostery from MD simulations of biochemical systems. MDiGest offers the ability to describe protein dynamics by combining different approaches, such as correlations of atomic displacements or dihedral angles, as well as a novel approach based on the correlation of Kabsch-Sander electrostatic couplings. MDiGest allows for comparisons of networks and community structures that capture physical information relevant to allostery. Multiple complementary tools for studying essential dynamics include principal component analysis, root mean square fluctuation, as well as secondary structure-based analyses.

Identifying non-nucleoside inhibitors of RNA-dependent RNA-polymerase of SARS-CoV-2 through per-residue energy decomposition-based pharmacophore modeling, molecular docking, and molecular dynamics simulation

BACKGROUND AND OBJECTIVE

The current coronavirus disease-2019 (COVID-19) pandemic has triggered a worldwide health and economic crisis. The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes the disease and completes its life cycle using the RNA-dependent RNA-polymerase (RdRp) enzyme, a prominent target for antivirals. In this study, we have computationally screened ∼690 million compounds from the ZINC20 database and 11,698 small molecule inhibitors from DrugBank to find existing and novel non-nucleoside inhibitors for SARS-CoV-2 RdRp.

METHODS

Herein, a combination of the structure-based pharmacophore modeling and hybrid virtual screening methods, including per-residue energy decomposition-based pharmacophore screening, molecular docking, pharmacokinetics, and toxicity evaluation were employed to retrieve novel as well as existing RdRp non-nucleoside inhibitors from large chemical databases. Besides, molecular dynamics simulation and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) method were used to investigate the binding stability and calculate the binding free energy of RdRp-inhibitor complexes.

RESULTS

Based on docking scores and significant binding interactions with crucial residues (Lys553, Arg557, Lys623, Cys815, and Ser816) in the RNA binding site of RdRp, three existing drugs, ZINC285540154, ZINC98208626, ZINC28467879, and five compounds from ZINC20 (ZINC739681614, ZINC1166211307, ZINC611516532, ZINC1602963057, and ZINC1398350200) were selected, and the conformational stability of RdRp due to their binding was confirmed through molecular dynamics simulation. The free energy calculations revealed these compounds possess strong binding affinities for RdRp. In addition, these novel inhibitors exhibited drug-like features, good absorption, distribution, metabolism, and excretion profile and were found to be non-toxic.

CONCLUSION

The compounds identified in the study by multifold computational strategy can be validated in vitro as potential non-nucleoside inhibitors of SARS-CoV-2 RdRp and holds promise for the discovery of novel drugs against COVID-19 in future.

Understanding Drug Resistance of Wild-Type and L38HL Insertion Mutant of HIV-1 C Protease to Saquinavir

Acquired immunodeficiency syndrome (AIDS) is one of the most challenging infectious diseases to treat on a global scale. Understanding the mechanisms underlying the development of drug resistance is necessary for novel therapeutics. HIV subtype C is known to harbor mutations at critical positions of HIV aspartic protease compared to HIV subtype B, which affects the binding affinity. Recently, a novel double-insertion mutation at codon 38 (L38HL) was characterized in HIV subtype C protease, whose effects on the interaction with protease inhibitors are hitherto unknown. In this study, the potential of L38HL double-insertion in HIV subtype C protease to induce a drug resistance phenotype towards the protease inhibitor, Saquinavir (SQV), was probed using various computational techniques, such as molecular dynamics simulations, binding free energy calculations, local conformational changes and principal component analysis. The results indicate that the L38HL mutation exhibits an increase in flexibility at the hinge and flap regions with a decrease in the binding affinity of SQV in comparison with wild-type HIV protease C. Further, we observed a wide opening at the binding site in the L38HL variant due to an alteration in flap dynamics, leading to a decrease in interactions with the binding site of the mutant protease. It is supported by an altered direction of motion of flap residues in the L38HL variant compared with the wild-type. These results provide deep insights into understanding the potential drug resistance phenotype in infected individuals.

Structural Insights into ATP-Sensitive Potassium Channel Mechanics: A Role of Intrinsically Disordered Regions

Commonly used techniques, such as CryoEM or X-ray, are not able to capture the structural reorganizations of disordered regions of proteins (IDR); therefore, it is difficult to assess their functions in proteins based exclusively on experiments. To fill this gap, we used computational molecular dynamics (MD) simulation methods to capture IDR dynamics and trace biological function-related interactions in the Kir6.2/SUR1 potassium channel. This ATP-sensitive octameric complex, one of the critical elements in the insulin secretion process in human pancreatic β-cells, has four to five large, disordered fragments. Using unique MD simulations of the full Kir6.2/SUR1 channel complex, we present an in-depth analysis of the dynamics of the disordered regions and discuss the possible functions they could have in this system. Our MD results confirmed the crucial role of the N-terminus of the Kir6.2 fragment and the L0-loop of the SUR1 protein in the transfer of mechanical signals between domains that trigger insulin release. Moreover, we show that the presence of IDRs affects natural ligand binding. Our research takes us one step further toward understanding the action of this vital complex.

In silico scanning of structural and functional deleterious nsSNPs in Arabidopsis thaliana's SOG1 protein, using molecular dynamic simulation approaches

Suppressor of gamma response 1 (SOG1) is a member of the NAC domain family transcription factors of the DNA damage response (DDR) signaling in the plant's genome. SOG1 is directly involved in transcriptional response to DNA damage, cell cycle checkpoints and ATR or ATM-mediated activation of the DNA damage responses and repair functioning in programmed cell death and regulation of end reduplication. Different mutations in the SOG1 protein lead to severe diseases and, ultimately, cell death. Single nucleotide polymorphisms (SNPs) are an important type of genetic alteration that cause different diseases or programmed cell death. The current study applied different computational approaches to Arabidopsis thaliana L. SOG1 protein to identify the potential deleterious nsSNPs and monitor their impact on the structure, function and protein stability. Various bioinformatics tools were applied to analyze the retrieved 34 nsSNPs and interestingly extracted four deleterious nsSNPs, that is, ensvath13968004 (Q166L), tmp18998388 (P159L), ensvath01103049 (K199N) and tmp18998295 (Y190F). For example, homology modeling, conservation and conformational analysis of the mutant's models were considered to scrutinize the deviations of these variants from the native SOG1 structure. All atoms molecular dynamic simulation confirmed the significance of these mutations on the protein stability, residual and structural conformation, compactness, surface conformation, dominant motion, Gibbs free energy distribution and dynamic effects. Similarly, protein-protein interaction revealed that SOG1 operates as a hub-linking cluster of various proteins, and any changes in the SOG1 might result in the disassociation of several signal transduction cascades.Communicated by Ramaswamy H. Sarma.

Structural basis of constitutive c-Src kinase activity due to R175L and W118A mutations

Cellular Src (c-Src) belongs to a non-receptor membrane-associated tyrosine kinase family that plays essential roles in cellular processes. Growing evidence suggests that R175L and W118A mutations in SH2/SH3 domains of c-Src functionally inactivate these domains leading to constitutive activation of kinase domain (KD). Here we modeled c-Src^R175L, c-Src^W118A and c-Src^W118A+R175L structures by inducing phosphorylation at Y416 or Y527, respectively to characterize the comparative dynamics in the active versus inactive states through molecular dynamics simulation assay. We observed more conformational readjustments in c-Src^open than its close variants. In particular, C-terminal tail residues of c-Src^W118A-open and c-Src^{W118A+R175L-open} demonstrate significantly higher transitions. The cross-correlation analysis revealed an anticorrelation behavior in the motion of KD with respect to SH2, SH3 and the linker region of Src^{W118A+R175L-open}, while in c-Src^WT-open, SH2 and SH3 domains were anticorrelated, while KD and C-terminal tail motions were correlated. Due to these conformational differences, c-Src open forms exhibited lower interaction between pY527 and SH2 domain. Through detailed structural analysis, we observed a uniform myristate binding cavity in c-Src^WT-open, while the myristoyl pockets of mutant forms were deformed. We propose that constitutive activation of mutant Src forms may presumably be achieved by the prolonged membrane binding due to unusual conformations of C-terminal and myristoyl switch residues that may result in a higher dephosphorylation rate at pY527 in the myristoylated c-Src. Thus, our study establishes novel clues to decipher the constitutive activation status of c-Src in response to known mutations that may help in devising novel therapeutic strategies for cancer metastasis treatment.Communicated by Ramaswamy H. Sarma.

SHP-1 tyrosine phosphatase binding to c-Src kinase phosphor-dependent conformations: A comparative structural framework

SHP-1 is a cytosolic tyrosine phosphatase that is primarily expressed in hematopoietic cells. It acts as a negative regulator of numerous signaling pathways and controls multiple cellular functions involved in cancer pathogenesis. This study describes the binding preferences of SHP-1 (pY536) to c-Srcopen (pY416) and c-Srcclose (pY527) through in silico approaches. Molecular dynamics simulation analysis revealed more conformational changes in c-Srcclose upon binding to SHP-1, as compared to its active/open conformation that is stabilized by the cooperative binding of the C-SH2 domain and C-terminal tail of SHP-1 to c-Src SH2 and KD. In contrast, c-Srcclose and SHP-1 interaction is mediated by PTP domain-specific WPD-loop (WPDXGXP) and Q-loop (QTXXQYXF) binding to c-Srcclose C-terminal tail residues. The dynamic correlation analysis demonstrated a positive correlation for SHP-1 PTP with KD, SH3, and the C-terminal tail of c-Srcclose. In the case of the c-Srcopen-SHP-1 complex, SH3 and SH2 domains of c-Srcopen were correlated to C-SH2 and the C-terminal tail of SHP-1. Our findings reveal that SHP1-dependent c-Src activation through dephosphorylation relies on the conformational shift in the inhibitory C-terminal tail that may ease the recruitment of the N-SH2 domain to phosphotyrosine residue, resulting in the relieving of the PTP domain. Collectively, this study delineates the intermolecular interaction paradigm and underlying conformational readjustments in SHP-1 due to binding with the c-Src active and inactive state. This study will largely help in devising novel therapeutic strategies for targeting cancer development.

Non-Ewald methods for evaluating the electrostatic interactions of charge systems: similarity and difference

In molecular simulations, it is essential to properly calculate the electrostatic interactions of particles in the physical system of interest. Here we consider a method called the non-Ewald method, which does not rely on the standard Ewald method with periodic boundary conditions, but instead relies on the cutoff-based techniques. We focus on the physicochemical and mathematical conceptual aspects of the method in order to gain a deeper understanding of the simulation methodology. In particular, we take into account the reaction field (RF) method, the isotropic periodic sum (IPS) method, and the zero-multipole summation method (ZMM). These cutoff-based methods are based on different physical ideas and are completely distinguishable in their underlying concepts. The RF and IPS methods are "additive" methods that incorporate information outside the cutoff region, via dielectric medium and isotropic boundary condition, respectively. In contrast, the ZMM is a "subtraction" method that tries to remove the artificial effects, generated near the boundary, from the cutoff sphere. Nonetheless, we find physical and/or mathematical similarities between these methods. In particular, the modified RF method can be derived by the principle of neutralization utilized in the ZMM, and we also found a direct relationship between IPS and ZMM.

Conformational Stability and Denaturation Processes of Proteins Investigated by Electrophoresis under Extreme Conditions

The functional structure of proteins results from marginally stable folded conformations. Reversible unfolding, irreversible denaturation, and deterioration can be caused by chemical and physical agents due to changes in the physicochemical conditions of pH, ionic strength, temperature, pressure, and electric field or due to the presence of a cosolvent that perturbs the delicate balance between stabilizing and destabilizing interactions and eventually induces chemical modifications. For most proteins, denaturation is a complex process involving transient intermediates in several reversible and eventually irreversible steps. Knowledge of protein stability and denaturation processes is mandatory for the development of enzymes as industrial catalysts, biopharmaceuticals, analytical and medical bioreagents, and safe industrial food. Electrophoresis techniques operating under extreme conditions are convenient tools for analyzing unfolding transitions, trapping transient intermediates, and gaining insight into the mechanisms of denaturation processes. Moreover, quantitative analysis of electrophoretic mobility transition curves allows the estimation of the conformational stability of proteins. These approaches include polyacrylamide gel electrophoresis and capillary zone electrophoresis under cold, heat, and hydrostatic pressure and in the presence of non-ionic denaturing agents or stabilizers such as polyols and heavy water. Lastly, after exposure to extremes of physical conditions, electrophoresis under standard conditions provides information on irreversible processes, slow conformational drifts, and slow renaturation processes. The impressive developments of enzyme technology with multiple applications in fine chemistry, biopharmaceutics, and nanomedicine prompted us to revisit the potentialities of these electrophoretic approaches. This feature review is illustrated with published and unpublished results obtained by the authors on cholinesterases and paraoxonase, two physiologically and toxicologically important enzymes.

Probing Alterations in MDM2 Catalytic Core Structure Effect of Garcinia Mangostana Derivatives: Insight from Molecular Dynamics Simulations

The MDM2-p53 protein-protein interaction is a promising model for researchers to design, study, and discover new anticancer drugs. The design of therapeutically active compounds that can maintain or restore the binding of MDM2 to p53 has been found to limit the oncogenic activities of both. This led to the current development of a group of xanthone-core and cis-imidazoline analogs compounds, among which γ-Mangostin (GM), α-Mangostin (AM), and Nutlin exhibited their MDM2-p53 interaction inhibitory effects. Therefore, in this study, we seek to determine the mechanisms by which these compounds elicit MDM2-p53 interaction targeting. Unique to the binding of GM, AM, and Nutlin, from our findings, they share the same three active site residues Val76, Tyr50, and Gly41, which represent the top active side residues that contribute to high electrostatic energy. Consequently, the free binding energy contributed enormously to the binding of these compounds, which culminated in the high binding affinities of GM, AM, and Nutlin with high values. Furthermore, GM, AM, and Nutlin commonly interrupted the stable and compact conformation of MDM2 coupled with its active site, where Cα deviations were relatively high. We believe that our findings would assist in the design of more potent active anticancer drugs.

CB1 as a novel target for Ginkgo biloba's terpene trilactone for controlling chemotherapy-induced peripheral neuropathy (CIPN)

The application of antineoplastic chemotherapeutic agents causes a common side effect known as chemotherapy-induced peripheral neuropathy (CIPN) that leads to reducing the quality of patient's life. This research involves the performance of molecular docking and molecular dynamic (MD) simulation studies to explore the impact of terpenoids of Ginkgo biloba on the targets (CB-1, TLR4, FAAH-1, COX-1, COX-2) that can significantly affect the controlling of CIPN's symptoms. According to the in-vitro and in-vivo investigations, terpenoids, particularly ginkgolides B, A, and bilobalide, can cause significant effects on neuropathic pain. The molecular docking results disclosed the tendency of our ligands to interact with mainly CB1 and FAAH-1, as well as partly with TLR4, throughout their interactions with targets. Terpene trilactone can exhibit a lower rate of binding energy than CB1's inhibitor (7dy), while being precisely located in the CB1's active site and capable of inducing stable interactions by forming hydrogen bonds. The analyses of MD simulation proved that ginkgolide B was a more suitable activator and inhibitor for CB1 and TLR4, respectively, when compared to bilobalide and ginkgolide A. Moreover, bilobalide is capable of inhibiting FAAH-1 more effectively than the two other ligands. According to the analyses of ADME, every three ligands followed the Lipinski's rule of five. Considering these facts, the exertion of three ligands is recommended for their anti-inflammatory, neuroprotective, and anti-nociception influences caused by primarily activating CB1 and inhibiting FAAH-1 and TLR4; in this regard, these compounds can stand as potential candidates for the control and treatment of CIPN's symptoms.

Recognition between CD147 and cyclophilin A deciphered by accelerated molecular dynamics simulations

CD147 functions as the receptor of extracellular cyclophilin A (CypA) in various diseases, and CD147-CypA binding ulteriorly underlies the pathological process of various viral infections including HIV-1, SARS, and SARS-CoV-2. Although CyPA has been identified as a key intermediate pro-inflammatory factor, the mechanism by which CD147 cooperates with CypA in the development of the cytokine storm remains largely unknown, and the binding profile of CD147 with CypA remains to be elucidated as well. Here, we prepared three binding models of the CD147-CypA complex, including the active site of CypA severally binding to the groove bound by the Ig1 and Ig2 domains (model-0), P180-G181 (model-1), and P211 (model-2) of CD147, as well as introducing mutations P180A-G181A and P211A individually in each model. All systems were studied using accelerated molecular dynamics simulations and the molecular mechanics generalized Born surface area (MM/GBSA) method. For model-0, CypA bound to the ectodomain of CD147 with the highest binding affinity. Moreover, mutations P180A-G181A of CD147 in model-0 decreased the binding affinity and weakened the dynamic correlation between CD147 and CypA, which resulted in CypA shifting from the initial binding location. Other residue mutations of CD147 did not significantly affect the CD147-CypA binding, as reflected by the energy and structural analyses. Compared with surface plasmon resonance results and nuclear magnetic resonance shift signals, CypA should tend to reciprocally bind to the groove of CD147, and the binding process might be modulated by P180-G181 rather than P211. Besides, residue R201 of CD147 is critical for CD147-CypA binding and needs further experimental verification. These findings further our understanding of the recruitment between CD147 and CypA and its potential role in the development of inflammation and viral infection.

Novel insights into the structural changes induced by disease-associated mutations in TDP-43: a computational approach

Over the last two decades, the pathogenic aggregation of TAR DNA-binding protein 43 (TDP-43) is found to be strongly associated with several fatal neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTD), etc. While the mutations and truncation in TDP-43 protein have been suggested to be responsible for TDP-43 pathogenesis by accelerating the aggregation process, the effects of these mutations on the bio-mechanism of pathological TDP-43 protein remained poorly understood. Investigating this at the molecular level, we formulized an integrated workflow of molecular dynamic simulation and machine learning models (MD-ML). By performing an extensive structural analysis of three disease-related mutations (i.e., I168A, D169G, and I168A-D169G) in the conserved RNA recognition motifs (RRM1) of TDP-43, we observed that the I168A-D169G double mutant delineates the highest packing of the protein inner core as compared to the other mutations, which may indicate more stability and higher chances of pathogenesis. Moreover, through our MD-ML workflow, we identified the biological descriptors of TDP-43 which includes the interacting residue pairs and individual protein residues that influence the stability of the protein and could be experimentally evaluated to develop potential therapeutic strategies.Communicated by Ramaswamy H. Sarma.

In Silico and In Vitro Evaluations of Fluorophoric Thiazolo-[2,3-b]quinazolinones as Anti-cancer Agents Targeting EGFR-TKD

Epidermal growth factor receptor tyrosine kinase domain (EGFR-TKD) plays a pivotal role in cellular signaling, growth, and metabolism. The EGFR-TKD is highly expressed in cancer cells and was endorsed as a therapeutic target for cancer management to overcome metastasis, cell proliferation, and angiogenesis. The novel thiazolo-[2,3-b]quinazolinones series were strategically developed by microwave-assisted organic synthesis and multi dominos reactions aimed to identify the potent thiazolo-[2,3-b]quinazolinone inhibitor against EGFR-TKD. This study explores the binding stability and binding strength of newly developed series via molecular docking, molecular dynamics simulation, and MM/PBSA and MM/GBSA calculations. The binding interaction was observed to be through the functional groups on aryl substituents at positions 3 and 5 of the thiazolo-[2, 3-b]quinazolinone scaffold. The methyl substituents at position 8 of the ligands had prominent hydrophobic interactions corroborating their bindings similar to the reference FDA-approved drug erlotinib in the active site. ADMET predictions reveal that derivatives 5ab, 5aq, and 5bq are drug-like and may be effective in in vitro study. Molecular dynamics simulation for 100 ns of docked complexes revealed their stability at the atomistic level. The ΔGbinding of thiazolo-[2,3-b]quinazolinone was found to be 5ab - 22.45, 5aq - 22.23, and 5bq - 20.76 similar to standard drug, and erlotinib - 24.11 kcal/mol was determined by MM/GBSA method. Furthermore, the anti-proliferative activity of leads of thiazolo-[2,3-b]quinazolinones (n = 3) was studied against breast cancer cell line (MCF-7) and non-small lung carcinoma cell line (H-1299). The highest inhibitions in cell proliferation were shown by 5bq derivatives, and the IC₅₀ was found to be 6.5 ± 0.67 µM against MCF-7 and 14.8 µM against H-1299. The noscapine was also taken as a positive control and showed IC₅₀ at higher concentrations 37 ± 1 against MCF-7 and 46.5 ± 1.2 against H-1299.

One for All, All for One: The Peculiar Dynamics of TNF-Receptor-Associated Factor (TRAF2) Subunits

TNF Receptor-Associated Factor 2 (TRAF2) is a homo-trimer belonging to the TNF-receptor-associated factor family (TRAFs). The TRAF2 oligomeric state is crucial for receptor binding, the interaction with other proteins (involved in the TNFR signaling), and the interaction with biological membranes. In this study, we present a computational analysis of the Molecular Dynamics of TRAF2-C (a truncated and soluble TRAF2 form) to identify patterns in the interactions between the three chains. We have performed a canonical analysis of the motion applied to molecular dynamics starting from the available crystal structure to identify correlated motions in TRAF2 dynamics. We have computed the displacement matrix, providing a frame-by-frame displacement for each residue in the dynamic. We provide the results in terms of the correlation matrix, which represents a detailed map of the correlated motions of residues. Eventually, we computed the so-called dynamical clusters, based on the Principal Component Analysis (PCA) of the motion (displacement) and the k means application on the first two principal components space. The results clearly indicate that, most of the time, two chains move in a strongly correlated motion, while the third chain follows a freer motion. A detailed analysis of the correlation matrix also shows that a few specific interface residues characterize the interaction of the more independent subunit with the other two. These findings suggest that the equilibrium between the trimer and the dissociated species (dimers and monomers) might be finely tuned by controlling a few critical residues in the protein quaternary structure, probably facilitating the regulation of oligomerization and dissociation in vivo.

Enhanced Conformational Sampling with an Adaptive Coarse-Grained Elastic Network Model Using Short-Time All-Atom Molecular Dynamics

Compared to all-atom molecular dynamics (AA-MD) simulations, coarse-grained (CG) MD simulations can significantly reduce calculation costs. However, existing CG-MD methods are unsuitable for sampling structures that depart significantly from the initial structure without any biased force. In this study, we developed a new adaptive CG elastic network model (ENM), in which the dynamic cross-correlation coefficient based on short-time AA-MD of at most ns order is considered. By applying Bayesian optimization to search for a suitable parameter among the vast parameter space of adaptive CG-ENM, we succeeded in reducing the searching cost to approximately 10% of those for random sampling and exhaustive sampling. To evaluate the performance of adaptive CG-ENM, we applied the new methodology to adenylate kinase (ADK) and glutamine binding protein (GBP) in the apo state. The results showed that the structural ensembles explored by adaptive CG-ENM could be considerably more diverse than those by conventional ENMs with enhanced sampling such as temperature replica exchange MD and long-time AA-MD of 1 μs. In particular, some of the structures sampled by adaptive ENM are relatively close to the holo-type structures of ADK and GBP. Furthermore, as a challenging task, to demonstrate the advantages of the CG model with lower calculation cost, we applied our new methodology to a larger biomolecule, integrin (αV) in the inactive state. Then, we sampled various structural ensembles, including extended structures that are apparently different from inactive ones.

Identification of potential target endoribonuclease NSP15 inhibitors of SARS-COV-2 from natural products through high-throughput virtual screening and molecular dynamics simulation

SARS‐CoV‐2 wreaks havoc around the world, triggering the COVID‐19 pandemic. It has been confirmed that the endoribonuclease NSP15 is crucial to the viral replication, and thus identified as a potential drug target against COVID‐19. The NSP15 protein was used as the target to conduct high‐throughput virtual screening on 30,926 natural products from the NPASS database to identify potential NSP15 inhibitors. And 100 ns molecular dynamics simulations were performed on the NSP15 and NSP15‐NPC198199 system. In all, 10 natural products with high docking scores with NSP15 protein were obtained, among which compound NPC198199 scored the highest. The analysis of the binding mode between NPC198199 and NSP15 found that NPC198199 would form H‐bond interactions with multiple key residues at the catalytic site. Subsequently, a series of post‐dynamics simulation analyses (including RMSD, RMSF, PCA, DCCM, RIN, binding free energy, and H‐bond occupancy) were performed to further explore inhibitory mechanism of compound NPC198199 on NSP15 protein at the molecular level. The research strongly indicates that the 10 natural compounds screened can be used as potential inhibitors of NSP15, and provides valuable information for the subsequent drug discovery of anti‐SARS‐CoV‐2.

Mechanistic insight into the impact of a bivalent ligand on the structure and dynamics of a GPCR oligomer

Development of effective bivalent ligands has become the focus of intensive research toward modulation of G protein-coupled receptor (GPCR) oligomers, particularly in the field of GPCR pharmacology. Experimental studies have shown that they increased binding affinity and signaling potency compared to their monovalent counterparts, yet underlying molecular mechanism remains elusive. To address this, we performed accelerated molecular dynamics simulations on bivalent-ligand bound Adenosine 2A receptor (A2AR) dimer in the context of a modeled tetramer, which consists of A2AR and dopamine 2 receptor (D2R) homodimers and their cognate G proteins. Our results demonstrate that bivalent ligand impacted interactions between pharmacophore groups and ligand binding residues, thus modulating allosteric communication network and water channel formed within the receptor. Moreover, it also strengthens contacts between receptor and G protein, by modulating the volume of ligand binding pocket and intracellular domain of the receptor. Importantly, we showed that impact evoked by the bivalent ligand on A2AR dimer was also transmitted to apo D2R, which is part of the neighboring D2R dimer. To the best of our knowledge, this is the first study that provides a mechanistic insight into the impact of a bivalent ligand on dynamics of a GPCR oligomer. Consequently, this will pave the way for development of effective ligands for modulation of GPCR oligomers and hence treatment of crucial diseases such as Parkinson's disease and cancer.

In Silico Drug Repurposing Approach: Investigation of Mycobacterium tuberculosis FadD32 Targeted by FDA-Approved Drugs

Background: Despite the enormous efforts made towards combating tuberculosis (TB), the disease remains a major global threat. Hence, new drugs with novel mechanisms against TB are urgently needed. Fatty acid degradation protein D32 (FadD32) has been identified as a promising drug target against TB, the protein is required for the biosynthesis of mycolic acids, hence, essential for the growth and multiplication of the mycobacterium. However, the FadD32 mechanism upon the binding of FDA-approved drugs is not well established. Herein, we applied virtual screening (VS), molecular docking, and molecular dynamic (MD) simulation to identify potential FDA-approved drugs against FadD32. Methodology/Results: VS technique was found promising to identify four FDA-approved drugs (accolate, sorafenib, mefloquine, and loperamide) with higher molecular docking scores, ranging from -8.0 to -10.0 kcal/mol. Post-MD analysis showed that the accolate hit displayed the highest total binding energy of -45.13 kcal/mol. Results also showed that the accolate hit formed more interactions with FadD32 active site residues and all active site residues displayed an increase in total binding contribution. RMSD, RMSF, Rg, and DCCM analysis further supported that the presence of accolate exhibited more structural stability, lower bimolecular flexibility, and more compactness into the FadD32 protein. Conclusions: Our study revealed accolate as the best potential drug against FadD32, hence a prospective anti-TB drug in TB therapy. In addition, we believe that the approach presented in the current study will serve as a cornerstone to identifying new potential inhibitors against a wide range of biological targets.

Photo-Switchable Sulfonylureas Binding to ATP-Sensitive Potassium Channel Reveal the Mechanism of Light-Controlled Insulin Release

ATP-sensitive potassium (KATP) channels are present in numerous organs, including the heart, brain, and pancreas. Physiological opening and closing of KATPs present in pancreatic β-cells, in response to changes in the ATP/ADP concentration ratio, are correlated with insulin release into the bloodstream. Sulfonylurea drugs, commonly used in type 2 diabetes mellitus treatment, bind to the octamer KATP channels composed of four pore-forming Kir6.2 and four SUR1 subunits and increase the probability of insulin release. Azobenzene-based derivatives of sulfonylureas, such as JB253 inspired by well-established antidiabetic drug glimepiride, allow for control of this process by light. The mechanism of that phenomenon was not known until now. In this paper, we use molecular docking, molecular dynamics, and metadynamics to reveal structural determinants explaining light-controlled insulin release. We show that both trans- and cis-JB253 bind to the same SUR1 cavity as antidiabetic sulfonylurea glibenclamide (GBM). Simulations indicate that, in contrast to trans-JB253, the cis-JB253 structure generated by blue light absorption promotes open structures of SUR1, in close similarity to the GBM effect. We postulate that in the open SUR1 structures, the N-terminal tail from Kir6.2 protruding into the SUR1 pocket is stabilized by flexible enough sulfonylureas. Therefore, the adjacent Kir6.2 pore is more often closed, which in turn facilitates insulin release. Thus, KATP conductance is regulated by peptide linkers between its Kir6.2 and SUR1 subunits, a phenomenon present in other biological signaling pathways. Our data explain the observed light-modulated activity of photoactive sulfonylureas and widen a way to develop new antidiabetic drugs having reduced adverse effects.

Inter-Subunit Dynamics Controls Tunnel Formation During the Oxygenation Process in Hemocyanin Hexamers

Hemocyanin from horseshoe crab in its active form is a homo-hexameric protein. It exists in open and closed conformations when transitioning between deoxygenated and oxygenated states. Here, we present a detailed dynamic atomistic investigation of the oxygenated and deoxygenated states of the hexameric hemocyanin using explicit solvent molecular dynamics simulations. We focus on the variation in solvent cavities and the formation of tunnels in the two conformational states. By employing principal component analysis and CVAE-based deep learning, we are able to differentiate between the dynamics of the deoxy- and oxygenated states of hemocyanin. Finally, our results identify the deoxygenated open conformation, which adopts a stable, closed conformation after the oxygenation process.

Structure-based identification of natural compound inhibitor against M. tuberculosis thioredoxin reductase: insight from molecular docking and dynamics simulation

Antioxidant systems of M. tuberculosis (Mtb) play an important role in providing resistance in the hostile environment of mononuclear phagocytes. Thioredoxin system is a known antioxidant system that consists of three copies of thioredoxins (Trxs) and a single copy of thioredoxin reductase (TrxR). TrxR has been validated as an essential gene known to be involved in the reduction of peroxides, dinitrobenzenes and hydroperoxides, and is crucial in maintaining the survival of Mtb in macrophages. Recently, it has been demonstrated to be a druggable target. In this study, molecular docking was applied to screen more than 20,000 natural compounds from the Traditional Chinese Medicine database. Theoretical calculation of ΔGbinding by the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) methods indicated two top-hit compounds that bind with a high affinity to the allosteric site, consisting of a hinge region, of TrxR. Further, stability and binding analysis of both compounds were carried out with molecular dynamics simulation. An analysis of conformational variation by principal component analysis (PCA) and protein contact network (PCN) uncovered the conformational changes in the compound-bound forms of protein. The NADPH domain formed many new interactions with the FAD domain in the compound-bound form, signifying that the binding may render an effect on the protein structure and function. Our results suggest that these two compounds could potentially be used for structure-based lead inhibitors against TrxR. The inhibitor selected as lead compound will be used further as a scaffold to optimize as novel anti-tuberculosis therapeutic.Communicated by Ramaswamy H. Sarma.

Camptothecin shows better promise than Curcumin in the inhibition of the Human Telomerase: A computational study

OBJECTIVES

The Human Telomerase enzyme has become a drug target in the treatment of cancers and age-related disorders. This study aims to identify potential natural inhibitors of the Human Telomerase from compounds derived from edible African plants.

MATERIALS AND METHODS

A library of 1,126 natural compounds was molecularly docked against the Telomerase Reverse Transcriptase (PDB ID: 5ugw), the catalytic subunit of the target protein. Curcumin, a known Telomerase inhibitor was used as the standard. The front-runner compounds were screened for bioavailability, pharmacokinetic properties, and bioactivity using the SWISSADME, PKCSM, and Molinspiration webservers respectively. The molecular dynamic simulation and analyses of the apo and holo proteins were performed by the Galaxy supercomputing webserver.

RESULTS

The results of the molecular docking and virtual screening reveal Augustamine and Camptothecin as lead compounds. Augustamine has better drug-likeness and pharmacokinetic properties while Camptothecin showed better bioactivity and stronger binding affinity (-8.2 kcal/mol) with the target. The holo structure formed by Camptothecin showed greater inhibitory activity against the target with a total RMSF of 169.853, B-Factor of 20.164, and 108 anti-correlating residues.

CONCLUSION

Though they both act at the same binding site, Camptothecin induces greater Telomerase inhibition and better molecular stability than the standard, Curcumin. Further tests are required to investigate the inhibitory activities of the lead compounds.

Computational Insights into the Deleterious Impacts of Missense Variants on N-Acetyl-d-glucosamine Kinase Structure and Function

An enzyme of the mammalian amino-sugar metabolism pathway, N-acetylglucosamine kinase (NAGK), that synthesizes N-acetylglucosamine (GlcNAc)-6-phosphate, is reported to promote dynein functions during mitosis, axonal and dendritic growth, cell migration, and selective autophagy, which all are unrelated to its enzyme activity. As non-enzymatic structural functions can be altered by genetic variation, we made an effort in this study aimed at deciphering the pathological effect of nonsynonymous single-nucleotide polymorphisms (nsSNPs) in NAGK gene. An integrated computational approach, including molecular dynamics (MD) simulation and protein-protein docking simulation, was used to identify the damaging nsSNPs and their detailed structural and functional consequences. The analysis revealed the four most damaging variants (G11R, G32R, G120E, and A156D), which are highly conserved and functional, positioned in both small (G11R and G32R) and large (G120E and A156D) domains of NAGK. G11R is located in the ATP binding region, while variants present in the large domain (G120E and A156D) were found to induce substantial alterations in the structural organizations of both domains, including the ATP and substrate binding sites. Furthermore, all variants were found to reduce binding energy between NAGK and dynein subunit DYNLRB1, as revealed by protein-protein docking and MM-GBSA binding energy calculation supporting their deleteriousness on non-canonical function. We hope these findings will direct future studies to gain more insight into the role of these variants in the loss of NAGK function and their role in neurodevelopmental disorders.

Insight into the Therapeutic Potential of a Bicyclic Hydroxypyridone Compound 2-[(2,4-Dichlorophenyl)methyl]-7-hydroxy-1,2,3,4-tetrahydro-8H-pyrido[1,2-a]pyrazin-8-one as COMT Inhibitor in the Treatment of Parkinson's Disease: A Molecular Dynamic Simulation Approach

Parkinson's disease (PD) is one of the most targeted neurodegenerative diseases in clinical research. Awareness of research is due to its increasing number of affected people worldwide. The pathology of PD has been linked to several key proteins upregulation such as the catechol O-Methyltransferase (COMT). Hence, the synthesis of compounds possessing inhibitory capacity has been the frontline of research in recent years. Several compounds have been synthesized among which is the nitrocatechol. However, major limitations associated with the nitrocatechol scaffold include the inability to possess adequate CNS penetration properties and hepatic toxicity associated with the compounds. However, a series of bicyclic hydroxypyridones compounds were synthesized to evaluate their inhibitory potentials on COMT protein with compound 38 (c38) 2-[(2,4-dichlorophenyl)methyl]-7-hydroxy-1,2,3,4-tetrahydro-8H-pyrido[1,2-a]pyrazin-8-one shown to have a 40 fold increase level coverage in its IC₅₀ over brain exposure when compared to the other synthesized compound. The molecular dynamics method was employed to understand the nature of interaction exhibited by c38. Molecular mechanics of c38 revealed a disruptive effect on the secondary structure of COMT protein. Per residue decomposition analysis revealed similar crucial residues involved in the favorable binding of c38 and tolcapone implicated its increased inhibitory capacity on COMT in preventing PD. Free binding energy (ΔG_bind ) of c38 further revealed the inhibitory capacity towards COMT protein in comparison to the FDA approved tolcapone. Ligand mobility analysis of both compounds showed a timewise different mobility pattern across the simulation time frame at the active site pocket of the protein connoting the different inhibitory potency exhibited by c38 and tolcapone. Findings from this study revealed optimization of c38 could facilitate the discovery of new compounds with enhanced inhibitory properties towards COMT in treating PD.

A Comparative Study to Decipher the Structural and Dynamics Determinants Underlying the Activity and Thermal Stability of GH-11 Xylanases

With the growing need for renewable sources of energy, the interest for enzymes capable of biomass degradation has been increasing. In this paper, we consider two different xylanases from the GH-11 family: the particularly active GH-11 xylanase from Neocallimastix patriciarum, NpXyn11A, and the hyper-thermostable mutant of the environmentally isolated GH-11 xylanase, EvXyn11^TS. Our aim is to identify the molecular determinants underlying the enhanced capacities of these two enzymes to ultimately graft the abilities of one on the other. Molecular dynamics simulations of the respective free-enzymes and enzyme–xylohexaose complexes were carried out at temperatures of 300, 340, and 500 K. An in-depth analysis of these MD simulations showed how differences in dynamics influence the activity and stability of these two enzymes and allowed us to study and understand in greater depth the molecular and structural basis of these two systems. In light of the results presented in this paper, the thumb region and the larger substrate binding cleft of NpXyn11A seem to play a major role on the activity of this enzyme. Its lower thermal stability may instead be caused by the higher flexibility of certain regions located further from the active site. Regions such as the N-ter, the loops located in the fingers region, the palm loop, and the helix loop seem to be less stable than in the hyper-thermostable EvXyn11^TS. By identifying molecular regions that are critical for the stability of these enzymes, this study allowed us to identify promising targets for engineering GH-11 xylanases. Eventually, we identify NpXyn11A as the ideal host for grafting the thermostabilizing traits of EvXyn11^TS.

Propeptide in Rhizopus chinensis Lipase: New Insights into Its Mechanism of Activity and Substrate Selectivity by Computational Design

Most fungal lipases contain a propeptide, which is very important for their function and substrate selectivity. In the present study, Rhizopus chinensis lipase (RCL) was used as a research model to explore the mechanism of the propeptide of the lipase. Conventional molecular dynamics (MD) and metadynamics simulations were used to explore the mechanism by which the propeptide affects the activity of the lipase, which was subsequently verified by mutation experiments. MD simulations indicated that the propeptide had an inhibitory effect on the lid movement of RCL and found a key region (Val5-Thr10) on the propeptide. Subsequently, site-directed mutations were created in this region. The mutations enhanced the lipase catalytic efficiency to 700% and showed the potential for the propeptide to shift the substrate specificity of RCL. The specificity and activity of RCL mutants also had similar trends to wild-type RCL toward triglycerides with varying chain lengths. The mutual corroboration of simulation and site-directed mutagenesis results revealed the vital role of the key propeptide region in the catalytic activity and substrate specificity of the lipase.

Exploring the conformational dynamics and flexibility of intrinsically disordered HIV-1 Nef protein using molecular dynamic network approaches

Intrinsically disordered proteins represent a class of proteins that lack fixed and well-defined three-dimensional structures in solution. HIV-1 Nef is an intrinsically disordered peripheral membrane protein involved in the replication and pathogenesis of HIV-1 infection. Nef controls expression levels of cell surface CD4 molecules that are essential for adaptive immunity. Despite the lack of fixed and stable structures, Nef physically interacts with the host cellular proteins (AP-1/MHC-I) and modulates intracellular trafficking pathways. Therefore, it is essential to understand how this dynamic conformational flexibility affects Nef structures and function. In this study, we combined all-atom molecular dynamics (MD) simulations and dynamic network approaches to better understand the structure and dynamics of Nef in two different forms, the free unbound and the bound state. Using the MD simulation approach, we show that the intrinsically disordered Nef exhibit a large dynamic field with more atomic fluctuations and lesser thermodynamic stability in the unbound conditions. The conformations of Nef change over time, and this protein remains more compact, folded, and stable in the bound form. The dynamic network analysis revealed regions of the protein capable of modulating the conformational behavior of the disordered Nef. The average betweenness centrality (BC) unveiled residues that are critical for mediating protein-protein interactions. The average shortest path length (L) and the perturbation response scanning exposed residues that are likely to be important in steering protein conformational changes. Overall, the study demonstrates how all-atom MD simulations combined with the dynamic network approach can be used to gain further insights into the structure and dynamics-function relationship of intrinsically disordered HIV-1 Nef.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02698-8.

Distinguishing the optimal binding mechanism through reversible and irreversible inhibition analysis of HSP72 protein in cancer therapy

Over the past two decades, covalent inhibitors have gained much interest and are living up to their reputation as a powerful tool in drug discovery. Covalent inhibitors possess several significant advantages, including increased biochemical efficiency, prolonged duration and the ability to target shallow, solvent-exposed substrate-binding domains. One of the enzymes that have been both covalently and non-covalently targeted is the heat shock protein 72 (HSP72). This elevated enzyme expression in cancer cells may be responsible for tumorigenesis and tumor progression by providing chemotherapy resistance. A critical gap remains in the molecular understanding of the structural mechanism's covalent and non-covalent binding to HSP72. In this study, we explore the most optimal binding mechanism in the inhibition of the HSP72. Based on the molecular dynamic analyses, it was evident that the non-covalent complex showed more stability than the covalent complex. The covalent ligand, however, was more able to induce and stabilize the sealed conformation of the HSP72-NBD ATP binding domain throughout the. Also, the non-covalent ligand does not induce any significant conformational change as it remained close to the shape of the unbound complex; and the affinity is only dependent on the multiple hydrogen bonds in contrast to the covalent ligand. This is supported by the secondary structure elements and principal component analysis that was more dominant in the covalently inhibited complex. Covalent bond induced the α-helices sealed conformation of the HSP72-NBD; based on our findings, this will prevent other small molecules from interacting at the ATP binding site domain. Moreover, inhibition of the ATP binding domain can directly affect the ATPs protein folding mechanism of the HSP72 enzyme. The essential dynamic analysis presented in this report compliments the binding mechanism of HSP72, establishing covalent inhibition as the preferred method of inhibiting the HSP72 protein. The findings from this study may assist in the design of more target-specific HSP72 covalent inhibitors exploring the surface-exposed lysine residues.

Structural Determination, Biological Function, and Molecular Modelling Studies of Sulfoaildenafil Adulterated in Herbal Dietary Supplement

Unapproved ingredients included in herbal medicines and dietary supplements have been detected as adulterated synthetic drugs used for erectile dysfunction. Extraction from a dietary supplement was performed to isolate the compounds by HPLC analysis. The structural characterization was confirmed using mass spectrometry (ESI-TOF/MS and LC-MS/MS), 1H NMR, and 13C NMR spectroscopy techniques. Results identified the thus-obtained compound to be sulfoaildenafil, a thioketone analogue of sildenafil. The biological activities of this active compound have been focused for the first time by the experimental point of view performance in vitro. The results revealed that sulfoaildenafil can affect the therapeutic level of nitric oxide through the upregulation of nitric oxide synthase and phosphodiesterase type 5 (PDE5) gene expressions. This bulk material, which displays structural similarity to sildenafil, was analyzed for the presence of a PDE5 inhibitor using a theoretical calculation. These unique features of the potential activity of PDE5 protein and its inhibitors, sildenafil and sulfoaildenafil, may play a key consideration for understanding the mode of actions and predicting the biological activities of PDE5 inhibitors.

Targeting ebola virus VP40 protein through novel inhibitors: exploring the structural and dynamic perspectives on molecular landscapes

Ebola filovirus (EBOV) is one of the deadliest known infectious agents, and a cause of Western African epidemics from 2013 to 2016. The virus has infected nearly 3000 humans and almost 1900 have died. In the past few years, various small molecules have been discovered to display efficiency against EBOV and some of them have progressed towards clinical trials. Even though continuous attempts have been made to find antiEBOV therapeutics, no potential drugs are yet approved against this viral infection. The development of small antiviral inhibitors has gained tremendous attention in the attempt to overcome EVD. With this background, we seek to offer molecular insights into EBOV VP40 protein inhibition, using all atom molecular mechanics methodology and binding free energy calculations. We have selected five novel reported inhibitors against VP40 protein, namely Comp1, Comp2, Comp3, Comp4, and Comp5, and explored their binding against the same target. It was evident from the analysis that all the inhibitors displayed stability in complex with VP40 protein; however, Comp1 exhibited enhanced stability and compactness. Comp1 unveiled favorable binding, which accounted for positive correlation motions in the active site residues. Likewise, Comp1 revealed the most promising binding (ΔG_bind - 40.3504 kcal/mol) as compared to the other four inhibitors, which disclosed relatively less favorable ΔG_bind. The highest binding energy of Comp1 to VP40 protein can be primarily endorsed to the upsurge in van der Waals energy by ΔE_vdW - 37.1609 kcal/mol and Coulomb energy by ΔE_ele - 52.7332 kcal/mol. Also, the hydrogen bond network is robust in Comp1-VP40 complex, with four hydrogen bonds, whilst it is less in other inhibitors. The outcomes from this report may assist in the advancement of novel VP40 inhibitors with high selectivity and potency for EVD therapeutics.

A phenylalanine dynamic switch controls the interfacial activation of Rhizopus chinensis lipase

The catalytic mechanism of most lipases involves a step called "interfacial activation" which significantly increases lipases activity beyond the critical micellar concentration (CMC) of substrate. In the present study, Rhizopus chinensis lipase (RCL) was used as a research model to explore the mechanism of lipase interfacial activation beyond the CMC. Molecular dynamic (MD) simulations indicated the open- and closed-lid transitions and revealed that Phe113 was the critical site for RCL activation by its dynamic flipping. Such putative switch affecting interfacial activation has not been reported in lipase so far. The function of Phe113 was subsequently verified by mutation experiments. The F113W mutant increases the lipase catalytic efficiency (1.9 s^-1·μM^-1) to 280% at the optimum temperature (40 °C) and pH 8.5 with the addition of 0.12 μg protein in the 200 μL reaction system. MD simulations indicated that the fast flipping rate from the closed to the open state, the high open state proportion, and the exposure of the catalytic triad are the main reasons for the lipase activation. The mutual corroboration of simulations and site-directed mutagenesis results revealed the vital role of Phe113 in the lipase activation.

Molecular modelling investigation for drugs and nutraceuticals against protease of SARS-CoV-2

The widespread problem of a 2019-novel coronavirus (SARS-CoV-2) strain outbreak in Wuhan, China has prompted a search for new drugs to protect against and treat this disease. It is necessary to immediately investigate this due to the mutation of the viral genome and there being no current protective vaccines or therapeutic drugs. Molecular modelling and molecular docking based on in silico screening strategies were employed to determine the potential activities of seven HIV protease (HIV-PR) inhibitors, two flu drugs, and eight natural compounds. The computational approach was carried out to discover the structural modes with a high binding affinity for these drugs on the homology structure of the Wuhan coronavirus protease (SARS-CoV-2 PR). From the theoretical calculations, all the drugs and natural compounds demonstrated various favorable binding affinities. An interesting finding was that the natural compounds tested had a higher potential binding activity with the pocket sites of SARS-CoV-2 PR compared to the groups of HIV-PR inhibitors. The binding modes of each complex illustrated between the drugs and compounds interacted with the functional group of amino acids in the binding pocket via hydrophilic, hydrophobic, and hydrogen bond interactions using the molecular dynamics simulation technique. This result supports the idea that existing protease inhibitors and natural compounds could be used to treat the new coronavirus. This report sought to provide fundamental knowledge as preliminary experimental data to propose an existing nutraceutical material against viral infection. Collectively, it is suggested that molecular modelling and molecular docking are suitable tools to search and screen for new drugs and natural compounds that can be used as future treatments for viral diseases.

Therapeutic Path to Double Knockout: Investigating the Selective Dual-Inhibitory Mechanisms of Adenosine Receptors A1 and A2 by a Novel Methoxy-Substituted Benzofuran Derivative in the Treatment of Parkinson's Disease

The dual inhibition of adenosine receptors A1 (A₁ AR) and A2 (A_2A AR) has been considered as an efficient strategy in the treatment of Parkinson's disease (PD). This led to the recent development of a series of methoxy-substituted benzofuran derivatives among which compound 3j exhibited dual-inhibitory potencies in the micromolar range. Therefore, in this study, we seek to resolve the mechanisms by which this novel compound elicits its selective dual targeting against A₁ AR and A_2A AR. Unique to the binding of 3j in both proteins, from our findings, is the ring-ring interaction elicited by _A1Phe275 (→ _A2Phe170) with the benzofuran ring of the compound. As observed, this π-stacking interaction contributes notably to the stability of 3j at the active sites of A₁ and A_2A AR. Besides, conserved active site residues in the proteins such as _A1Ala170 (→ _A2Ala65), _A1Ile173 (→ _A2Ile68), _A1Val191 (→ _A2Val86), _A1Leu192 (→ _A2Leu87), _A1Ala195 (→ _A2Ala90), _A1Met284 (→ _A2Met179), _A1Tyr375 (→ _A2Tyr369), _A1Ile378 (→ _A2Ile372), and _A1His382 (→ _A2His376) were commonly involved with other ring substituents which further complement the dual binding and stability of 3j. This reflects a similar interaction mechanism that involved aromatic (π) interactions. Consequentially, vdW energies contributed immensely to the dual binding of the compound, which culminated in high ΔG_binds that were homogenous in both proteins. Furthermore, 3j commonly disrupted the stable and compact conformation of A₁ and A_2A AR, coupled with their active sites where Cα deviations were relatively high. Ligand mobility analysis also revealed that both compounds exhibited a similar motion pattern at the active site of the proteins relative to their optimal dual binding. We believe that findings from this study with significantly aid the structure-based design of highly selective dual-inhibitors of A₁ and A_2A AR.

A quantum biochemistry investigation of the protein-protein interactions for the description of allosteric modulation on biomass-degrading chimera

The worldwide dependence of population on fossil fuels continues to have several harmful implications for the environment. Bioethanol is an excellent option for renewable fuel to replace the current greenhouse gas emitters. In addition, its production by enzymatic route has gained space among the industrial processes because it replaces the traditional acid treatment. Due to its high versatility, the xylanase family is used in this process as an accessory enzyme for degrading the lignocellulosic substrate of biomass. A chimera built by a xylanolytic domain (Xyl) and a xylose-binding protein (XBP) showed an experimentally improved catalytic efficiency and interdomain allosteric modulation after xylose binding. In this context, we performed a quantum biochemistry characterization of the interactions between these domains and dynamic cross-correlation (DCC) analysis after performing molecular dynamics (DM) simulations of the systems in the presence and absence of xylose in the XBP active site. We used the density functional theory (DFT) within the molecular fractionation with the conjugated caps (MFCC) approach to describe the pair energies, and the corresponding energy difference between the chimera domains responsible for the allosteric effect and amino acid DCC to evaluate the interdomain coupling differences between the energy states. The detailed energetic investigation together with the related structural and dynamics counterparts revealed the molecular mechanisms of chimeric improvement of the xylanase activity observed experimentally. This mechanism was correlated with greater stability and high connectivity at the interdomain interface in the xylose bound relative to the free chimera. We identify the contributions of hydrogen bonds, hydrophobic interactions and water-mediated interactions in the interdomain region responsible for stability together with the structural and dynamical elements related to the allosteric effect. Taken together, these observations led to a comprehensive understanding of the chimera's modulatory action that occurs through the formation of a highly connected interface that makes the essential movements related to xylanolytic activity in xylanase correlated to those of the xylose-binding protein.

Structural Determinants of Insulin Release: Disordered N-Terminal Tail of Kir6.2 Affects Potassium Channel Dynamics through Interactions with Sulfonylurea Binding Region in a SUR1 Partner

Inward rectifying potassium ion channels (KATP), sensitive to the ATP/ADP concentration ratio, play an important, control role in pancreatic β cells. The channels close upon the increase of this ratio, which, in turn, triggers insulin release to blood. Numerous mutations in KATP lead to severe and widespread medical conditions such as diabetes. The KATP system consists of a pore made of four Kir6.2 subunits and four accompanying large SUR1 proteins belonging to the ABCC transporters group. How SUR1 affects KATP function is not yet known; therefore, we created simplified models of the Kir6.2 tetramer based on recently determined cryo-EM KATP structures. Using all-atom molecular dynamics (MD) with the CHARMM36 force field, targeted MD, and molecular docking, we revealed functionally important rearrangements in the Kir6.2 pore, induced by the presence of the SUR1 protein. The cytoplasmic domain of Kir6.2 (CTD) is brought closer to the membrane due to interactions with SUR1. Each Kir6.2 subunit has a conserved, functionally important, disordered N-terminal tail. Using molecular docking, we found that the Kir6.2 tail easily docks to the sulfonylurea drug binding region located in the adjacent SUR1 protein. We reveal, for the first time, dynamical behavior of the Kir6.2/SUR1 system, confirming a physiological role of the Kir6.2 disordered tail, and we indicate structural determinants of KATP-dependent insulin release from pancreatic β cells.

Long-range Regulation of Partially Folded Amyloidogenic Peptides

Neurodegeneration involves abnormal aggregation of intrinsically disordered amyloidogenic peptides (IDPs), usually mediated by hydrophobic protein-protein interactions. There is mounting evidence that formation of α-helical intermediates is an early event during self-assembly of amyloid-β42 (Aβ42) and α-synuclein (αS) IDPs in Alzheimer’s and Parkinson’s disease pathogenesis, respectively. However, the driving force behind on-pathway molecular assembly of partially folded helical monomers into helical oligomers assembly remains unknown. Here, we employ extensive molecular dynamics simulations to sample the helical conformational sub-spaces of monomeric peptides of both Aβ42 and αS. Our computed free energies, population shifts, and dynamic cross-correlation network analyses reveal a common feature of long-range intra-peptide modulation of partial helical folds of the amyloidogenic central hydrophobic domains via concerted coupling with their charged terminal tails (N-terminus of Aβ42 and C-terminus of αS). The absence of such inter-domain fluctuations in both fully helical and completely unfolded (disordered) states suggests that long-range coupling regulates the dynamicity of partially folded helices, in both Aβ42 and αS peptides. The inter-domain coupling suggests a form of intra-molecular allosteric regulation of the aggregation trigger in partially folded helical monomers. This approach could be applied to study the broad range of amyloidogenic peptides, which could provide a new path to curbing pathogenic aggregation of partially folded conformers into oligomers, by inhibition of sites far from the hydrophobic core.