Structural and energetic insights into the selective inhibition of PKMYT1 against WEE1

Protein kinase, membrane-associated tyrosine/threonine 1 (PKMYT1), a member of the WEE family and responsible for the regulation of CDK1 phosphorylation, has been considered a promising therapeutic target for cancer therapy. However, the highly structural conservation of the ATP-binding sites of the WEE family poses a challenge to the design of selective inhibitors for PKMYT1. Here, molecular docking, multiple microsecond-length molecular dynamics (MD) simulations and end-point free energy calculations were performed to uncover the molecular mechanism of the binding selectivity of RP-6306 toward PKMYT1 over its highly homologous kinase WEE1. The binding specificity of RP-6306 reported in previous experimental bioassays was clarified by MD simulations and binding free energy calculations. Further, the binding free energy prediction indicated that the binding selectivity of RP-6306 largely derived from the difference in the protein-ligand electrostatic interactions. The per-residue free energy decomposition suggested that the non-conserved gatekeeper residue in the hinge domain of PKMYT1/WEE1, Thr187/Asn376, is the critical factor responsible for the binding selectivity of RP-6306 toward PKMYT1. In addition, a water-mediated hydrogen bond was formed between RP-6306 and Gly191 at the hinge domain in the PKMYT1/RP-6306 complex, which was absent in the WEE1/RP-6306 complex. This study is expected to offer useful information for the design of more potent and selective PKMYT1 inhibitors.Communicated by Ramaswamy H. Sarma.

MD investigation on the differences in the dynamic interactions between the specific ligand azamulin and two CYP3A isoforms, 3A4 and 3A5

The unmarked potential drug molecule azamulin has been found to be a specific inhibitor of CYP3A4 and CYP3A5 in recent years, but this molecule also shows different binding ability and affinity to the two CYP3A isoforms. In order to explore the microscopic mechanism, conventional molecular dynamics (MD) simulation methods were performed to study the dynamic interactions between two isoforms and azamulin. The simulation results show that the binding of the ligand leads to different structural properties of two CYP3A proteins. First of all, compared with apo-CYP3A4, the binding of the ligand azamulin can lead to changes in the structural flexibility of CYP3A4, i.e., holo-CYP3A4 is more flexible than apo-CYP3A4. The structural changes of CYP3A5 are just the opposite. The ligand binding increases the rigidity of CYP3A5. Furthermore, the representative structures of the production phase in the MD simulation were in details analyzed to obtain the microscopic interactions between the ligand azamulin and two CYP3A isoforms at the atomic level. It is speculated that the difference of composition and interaction of the active sites is the fundamental cause of the change of structural properties of the two proteins.Communicated by Ramaswamy H. Sarma.

A λ-dynamics investigation of insulin Wakayama and other A3 variant binding affinities to the insulin receptor

Insulin Wakayama is a clinical insulin variant where a conserved valine at the third residue on insulin's A chain (ValA3) is replaced with a leucine (LeuA3), impairing insulin receptor (IR) binding by 140-500 fold. This severe impact on binding from such a subtle modification has posed an intriguing problem for decades. Although experimental investigations of natural and unnatural A3 mutations have highlighted the sensitivity of insulin-IR binding to minor changes at this site, an atomistic explanation of these binding trends has remained elusive. We investigate this problem computationally using λ-dynamics free energy calculations to model structural changes in response to perturbations of the ValA3 side chain and to calculate associated relative changes in binding free energy (ΔΔGbind). The Wakayama LeuA3 mutation and seven other A3 substitutions were studied in this work. The calculated ΔΔGbind results showed high agreement compared to experimental binding potencies with a Pearson correlation of 0.88 and a mean unsigned error of 0.68 kcal/mol. Extensive structural analyses of λ-dynamics trajectories revealed that critical interactions were disrupted between insulin and the insulin receptor as a result of the A3 mutations. This investigation also quantifies the effect that adding an A3 Cδ atom or losing an A3 Cγ atom has on insulin's binding affinity to the IR. Thus, λ-dynamics was able to successfully model the effects of subtle modifications to insulin's A3 side chain on its protein-protein interactions with the IR and shed new light on a decades-old mystery: the exquisite sensitivity of hormone-receptor binding to a subtle modification of an invariant insulin residue.

Identification of Potential Inhibitors of the SARS-CoV-2 NSP13 Helicase via Structure-Based Ligand Design, Molecular Docking and Nonequilibrium Alchemical Simulations

We have assembled a computational pipeline based on virtual screening, docking techniques, and nonequilibrium molecular dynamics simulations, with the goal of identifying possible inhibitors of the SARS-CoV-2 NSP13 helicase, catalyzing by ATP hydrolysis the unwinding of double or single-stranded RNA in the viral replication process inside the host cell. The druggable sites for broad-spectrum inhibitors are represented by the RNA binding sites at the 5' entrance and 3' exit of the central channel, a structural motif that is highly conserved across coronaviruses. Potential binders were first generated using structure-based ligand techniques. Their potency was estimated by using four popular docking scoring functions. Common docking hits for NSP13 were finally tested using advanced nonequilibrium alchemical techniques for binding free energy calculations on a high-performing parallel cluster. Four potential NSP13 inhibitors with potency from submicrimolar to nanomolar were finally identified.

Molecular modeling of novel 2-aminopyridine derivatives as potential JAK2 inhibitors: a rational strategy for promising anticancer agents

Janus kinase 2(JAK2) is a potential target for anticancer drugs in the treatment of numerous myeloproliferative diseases due to its central role in the JAK/STAT signaling cascade. In this study, the binding behavior of 2 amino-pyridine derivatives as JAK2 inhibitors was investigated by using multifaceted strategies including 3D-QSAR, molecular docking, Fingerprint analysis, MD simulations, and MM-PBSA calculations. A credible COMFA (q2 = 0.606 and r2 = 0.919) and COMSIA (q2 = 0.641 and r2 = 0.992) model was developed, where the internal and external validation revealed that the obtained 3D-QSAR models could be capable of predicting bioactivities of JAK2 inhibitors. The structural criteria provided by the contour maps of model were used to computationally develop more potent 100 new JAK2 inhibitors. Docking studies were conducted on the model data set and newly developed compounds (in-house library) to demonstrate their binding mechanism and highlight the key interacting residues within JAK2 active site. The selected docked complexes underwent MD simulation (100 ns), which contributed in the further study of the binding interactions. Binding free energy analyses (MMGB/PBSA) revealed that key residues such as Glu930, Leu932 (hinge region), Asp939 (solvent accessible region), Arg980, Asn981and Asp994 (catalytic site) have a significantly facilitate ligand-protein interactions through H-bonding and van der Waals interactions. The preliminary in-silico ADMET evaluation revealed encouraging results for all the modeled and in-house library compounds. The findings of this research have the potential to offer valuable recommendations for the advancement of novel, potent, and efficacious JAK2 inhibitors. Overall, this work has successfully employed a wide range of computer-based methodologies to understand the interaction dynamics between 2-amino-pyridine derivatives and the JAK2 enzyme, which is a crucial target in myeloproliferative disorders.Communicated by Ramaswamy H. Sarma.

Identification of indole-based natural compounds as inhibitors of PARP-1 against triple-negative breast cancer: a computational study

Triple-negative breast cancer (TNBC) is the most aggressive kind of breast cancer known to mankind. It is a heterogeneous disease that is formed due to the missing estrogen, progesterone and human epidermal growth factor 2 receptors. Poly(ADP-ribose) polymerase-1 (PARP-1) protein helps in the development of TNBC by repairing the cancer cells, which proliferate and spread metastatically. To determine the potential PARP-1 inhibitors (PARPi), 0.2 million natural products from Universal Natural Product Database were screened using molecular docking and six hit compounds were selected based on their binding affinity towards PARP-1. The bio-availability and drug-like properties of these natural products were evaluated using ADMET analysis. Molecular dynamics simulations were conducted for these complexes for 200 ns to examine their structural stability and dynamic behaviour and further compared with the complex of talazoparib (TALA), an FDA-approved PARPi. Using MM/PBSA calculations, we conclude that the complexes HIT-3 and HIT-5 (-25.64 and -23.14 kcal/mol, respectively) show stronger binding energies with PARP-1 than TALA with PARP-1 (-10.74 kcal/mol). Strong interactions were observed between the compounds and hotspot residues, Asp770, Ala880, Tyr889, Tyr896, Ala898, Asp899 and Tyr907, of PARP-1 due to the existence of various types of non-covalent interactions between the compounds and PARP-1. This research offers critical information about PARPi, which could potentially be incorporated into the treatment of TNBC. Moreover, these findings were validated by comparing them with an FDA-approved PARPi.

Mechanistic and structural insight into R2R3-MYB transcription factor in plants: molecular dynamics based binding free energy analysis

The MYB transcription factor (TF) family is essential for various plant growth and development processes, including responses to biotic and abiotic stresses. This study investigated the R2R3-MYB protein structure from five plants, including cereal crops. The R2R3-MYB protein structure was docked with the DNA structure, and the best complexes were selected for two runs of molecular dynamics (MD) simulations to investigate the key interacting residues and conformational changes in the R2R3-MYB proteins caused by DNA binding. The MM/PBSA method calculated the binding free energy for each R2R3-MYB protein-DNA complex, showing strong interaction. Hydrophobic and hydrogen bonds significantly stabilized the R2R3-MYB protein-DNA complexes. The principal component analysis showed high restrictions on the movement of protein atoms in the phase space. A similar MD simulation analysis was performed using the crystal structure of the R2R3-MYB protein-DNA complex from Arabidopsis thaliana, and the generated complexes resembled the X-ray crystal structure. This is the first detailed study on the R2R3-MYB protein-DNA complex in cereal crops, providing a cost-effective solution to identify the key interacting residues and analyze the conformational changes in the MYB domain before and after DNA binding.Communicated by Ramaswamy H. Sarma.

3D-QSAR, docking and molecular dynamics simulations of novel Pyrazolo-pyridazinone derivatives as covalent inhibitors of FGFR1: a scientific approach for possible anticancer agents

Developing highly potent covalent inhibitors of Fibroblast growth factor receptors 1 (FGFR1) has always been a challenging task. In the current study, various computational techniques, such as 3D-QSAR, covalent docking, fingerprinting analysis, MD simulation followed by MMGB/PBSA, and per-residue energy decomposition analysis were used to explore the binding mechanism of pyrazolo[3,4-d]pyridazinone derivatives to FGFR1. The high q2 and r2 values for the CoMFA and CoMSIA models, suggest that the constructed 3D-QSAR models could reliably predict the bioactivities of FGFR1 inhibitors. The structural requirements revealed by the model's contour maps were strategically used to computationally create an in-house library of more than 100 new FGFR1 inhibitors using the R-group exploration technique implemented in the SparkTM software. The compounds from the in-house library were also mapped in the 3D-QSAR model that predicts comparable pIC50 values with the experimental values. A comparison between 3D-QSAR generated contours and molecular docking conformation of ligands was performed to reveal the fundamentals to design potent FGFR1 covalent inhibitors. The estimated binding free energies (MMGB/PBSA) for the selected compounds were in agreement with the experimental value ranking of their binding affinities towards FGFR1. Furthermore, per-residue energy decomposition analysis has identified Arg627 and Glu531 to contribute significantly in improved binding affinity of compound W16. During ADME analysis, the majority of in-house library compounds exhibited pharmacokinetic properties superior to those of experimentally produced compounds. These new compounds may help researchers better understand FGFR1 inhibition and lead to the creation of novel, potent FGFR1 inhibitors.Communicated by Ramaswamy H. Sarma.

Efficiency of membrane fusion inhibitors on different hemagglutinin subtypes: insight from a molecular dynamics simulation perspective

The challenge in vaccine development, along with drug resistance issues, has encouraged the search for new anti-influenza drugs targeting different viral proteins. Hemagglutinin (HA) glycoprotein, crucial in the viral replication cycle, has emerged as a promising therapeutic target. CBS1117 and JNJ4796 were reported to exhibit similar potencies against infectious group 1 influenza, which included H1 and H5 HAs; however, their potencies were significantly reduced against group 2 HA. This study aims to explore the molecular binding mechanisms and group specificity of these fusion inhibitors against both group 1 (H5) and group 2 (H3) HA influenza viruses using molecular dynamics simulations. CBS1117 and JNJ4796 exhibit stronger interactions with key residues within the H5 HA binding pocket compared to H3-ligand complexes. Hydrogen bonding and hydrophobic interactions involving residues, such as H381, Q401, T3251 (H5-CBS1117), T3181 (H5-JNJ4796), W212, I452, V482, and V522 predominantly contribute to stabilizing H5-ligand systems. In contrast, these interactions are notably weakened in H3-inhibitor complexes. Predicted protein-ligand binding free energies align with experimental data, indicating CBS1117 and JNJ4796's preference for heterosubtypic group 1 HA binding. Understanding the detailed atomistic mechanisms behind the varying potencies of these inhibitors against the two HA groups can significantly contribute to the development and optimization of effective HA fusion inhibitors. To accomplish this, the knowledge of the transition of HA from its pre- to post-fusion states, the molecular size of ligands, and their potential binding regions, could be carefully considered.Communicated by Ramaswamy H. Sarma.

Revealing the binding mechanism of BACE1 inhibitors through molecular dynamics simulations

Alzheimer's disease is a debilitating neurodegenerative disorder, and the Beta-Site Amyloid Precursor Protein Cleaving Enzyme 1 (BACE1) is a key therapeutic target in its treatment. This study employs molecular dynamics simulations and binding energy analysis to investigate the binding interactions between BACE1 and four selected small molecules: CNP520, D9W, NB641, and NB360. The binding model analysis indicates that the binding of BACE1 with four molecules are stable, except the loop regions show significant fluctuation. The binding free energy analyses reveal that NB360 exhibits the highest binding affinity with BACE1, surpassing other molecules (CNP520, D9W, and NB641). Detailed energy component assessments highlight the critical roles of electrostatic interactions and van der Waals forces in the binding process. Furthermore, residue contribution analysis identifies key amino acids influencing the binding process across all systems. Hydrogen bond analysis reveals a limited number of bonds between BACE1 and each small molecule, highlighting the importance of structural modifications to enable more stable hydrogen bonds. This research provides valuable insights into the molecular mechanisms of potential Alzheimer's disease therapeutics, guiding the way for improved drug design and the development of effective treatments targeting BACE1.Communicated by Ramaswamy H. Sarma.

An Efficient Approach to the Accurate Prediction of Mutational Effects in Antigen Binding to the MHC1

The major histocompatibility complex (MHC) can recognize and bind to external peptides to generate effective immune responses by presenting the peptides to T cells. Therefore, understanding the binding modes of peptide-MHC complexes (pMHC) and predicting the binding affinity of pMHCs play a crucial role in the rational design of peptide vaccines. In this study, we employed molecular dynamics (MD) simulations and free energy calculations with an Alanine Scanning with Generalized Born and Interaction Entropy (ASGBIE) method to investigate the protein-peptide interaction between HLA-A*02:01 and the G9209 peptide derived from the melanoma antigen gp100. The energy contribution of individual residue was calculated using alanine scanning, and hotspots on both the MHC and the peptides were identified. Our study shows that the pMHC binding is dominated by the van der Waals interactions. Furthermore, we optimized the ASGBIE method, achieving a Pearson correlation coefficient of 0.91 between predicted and experimental binding affinity for mutated antigens. This represents a significant improvement over the conventional MM/GBSA method, which yields a Pearson correlation coefficient of 0.22. The computational protocol developed in this study can be applied to the computational screening of antigens for the MHC1 as well as other protein-peptide binding systems.

Molecular Dynamics Simulation of an Iron(III) Binding Site on the Fc Domain of IgG1 Relevant for Visible Light-Induced Protein Fragmentation

Molecular dynamics simulations were employed to investigate the interaction between Fe(III) and an iron-binding site composed of THR259, ASP252, and GLU261 on the Fc domain of an IgG1. The goal was to provide microscopic mechanistic information for the photochemical, iron-dependent site-specific oxidative fragmentation of IgG1 at THR259 reported in our previous paper. The distance between Fe(III) and residues of interest as well as the binding pocket size was examined for both protonated and deprotonated THR259. The Fe(III) binding free energy (ΔG) was estimated by using an umbrella sampling approach. The pKa shift of the THR259 hydroxyl group caused by the presence of nearby Fe(III) was estimated based on a thermodynamic cycle. The simulation results show that Fe(III) resides inside the proposed binding pocket and profoundly changes the pocket configuration. The ΔG values indicate that the pocket possesses a strong binding affinity for Fe(III). Furthermore, Fe(III) profoundly lowers the pKa value of the THR259 hydroxyl group by 5.4 pKa units.

Structural basis for the mechanism of interaction of SARS-CoV-2 B.1.640.2 variant RBD with the host receptors hACE2 and GRP78

The inflicted chaos instigated by the SARS-CoV-2 (Severe acute respiratory syndrome coronavirus 2) globally continues with the emergence of novel variants. The current global outbreak is aggravated by the manifestation of novel variants, which affect the effectiveness of the vaccine, attachment with hACE2 (human Angiotensin-converting enzyme 2) and immune evasion. Recently, a new variant named University Hospital Institute (IHU) (B.1.640.2) was reported in France in November 2021 and is spreading globally affecting public healthcare. The B.1.640.2 SARS-CoV-2 strain revealed 14 mutations and 9 deletions in spike protein. Thus, it is important to understand how these variations in the spike protein impact the communication with the host. A protein coupling approach along with molecular simulation protocols was used to interpret the variation in the binding of the wild type (WT) and B.1.640.2 variant with hACE2 and Glucose-regulating protein 78 (GRP78) receptors. The initial docking scores revealed a stronger binding of the B.1.640.2-RBD with both the hACE2 and GRP78. To further understand the crucial dynamic changes, we looked at the structural and dynamic characteristics and also explored the variations in the bonding networks between the WT and B.1.640.2-RBD (receptor-binding domain) in association with hACE2 and GRP78, respectively. Our findings revealed that the variant complex demonstrated distinct dynamic properties in contrast to the wild type due to the acquired mutations. Finally, to provide conclusive evidence on the higher binding by the B.1.640.2 variant the TBE was computed for each complex. For the WT with hACE2 the TBE was quantified to be-61.38 ± 0.96 kcal/mol and for B.1.640.2 variant the TBE was estimated to be -70.47 ± 1.00 kcal/mol. For the WT-RBD-GRP78 the TBE -was computed to be 32.32 ± 0.56 kcal/mol and for the B.1.640.2-RBD a TBE of -50.39 ± 0.88 kcal/mol was reported. This show that these mutations are the basis for higher binding and infectivity produced by B.1.640.2 variant and can be targeted for drug designing against it.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulation study on the binding mechanism between carbon nanotubes and RNA-dependent RNA polymerase

Carbon nanotubes (CNTs) have potential prospects in disease treatment, so it is of great significance to study CNTs as the possible inhibitors of RNA-dependent RNA polymerase (RdRp). Through the way of using the RdRp of SARS-COV-2 as a model, five armchair single-walled carbon nanotubes (SWCNTs) (namely Dn, which stands for CNTs (n, m = n), n = 3-7) and RdRp have been selected to study the interactions by means of molecular docking and molecular dynamics simulation. After five SWCNT-RdRp complex systems have been subjected to the molecular dynamics simulations of 100 ns, and Molecular Mechanics Poisson - Boltzmann Surface Area (MMPBSA) has been used to calculate the binding free energy, it is found that the binding free energy of the D6 system (-189.541 kJ/mol) is significantly higher than that of the other four systems, and most of the amino acids with strong positive effects on binding are usually basic amino acids. What's more, in the further investigation of the specific interaction mechanism between CNT (6,6) and RdRp, it is revealed that the three amino acid residues LYS545, ARG553 and ARG555 located in the nucleoside triphosphate (NTP) entry channel all have strong effects. In addition, it is also observed that when ARG555 has been inserted into SWCNT, a stable structure will be formed, which will break the original NTP entry channel structure and inhibit virus replication. Therefore, it can be concluded that certain specific types of SWCNT, such as CNT (6,6), could be potential small molecule inhibitors in the treatment of coronavirus.Communicated by Ramaswamy H. Sarma.

Exploring binding potential of two new indole alkaloids from Nauclea officinalis against third and fourth generation EGFR: druglikeness, in silico ADMET, docking, DFT, molecular dynamics simulation, and MMGBSA study

This study investigates the anti-cancer potential of recently discovered indole alkaloids from Nauclea Officinalis against third and fourth-generation EGFR mutations using computational tools. Through ADMET profiling, druglikeness prediction, docking, and simulations, we assessed their pharmacokinetics, binding interactions, and stability. Promising druglikeness and binding affinity were observed, particularly for (±)-19-O-butylangustoline, which demonstrated stronger binding against both EGFR mutants. MD simulations confirmed stable interactions, with (±)-19-O-butylangustoline exhibiting the highest stability. These findings highlight these indole alkaloids as potential anti-cancer agents, with (±)-19-O-butylangustoline warranting further optimisation for therapeutic development. This study informs their potential through insights into molecular properties and binding energetics.

Exploring the natural products chemical space to abrogate the F3L-dsRNA interface of monkeypox virus to enhance the immune responses using molecular screening and free energy calculations

Amid the ongoing monkeypox outbreak, there is an urgent need for the rapid development of effective therapeutic interventions capable of countering the immune evasion mechanisms employed by the monkeypox virus (MPXV). The evasion strategy involves the binding of the F3L protein to dsRNA, resulting in diminished interferon (IFN) production. Consequently, our current research focuses on utilizing virtual drug screening techniques to target the RNA binding domain of the F3L protein. Out of the 954 compounds within the South African natural compound database, only four demonstrated notable docking scores: -6.55, -6.47, -6.37, and -6.35 kcal/mol. The dissociation constant (KD) analysis revealed a stronger binding affinity of the top hits 1-4 (-5.34, -5.32, -5.29, and -5.36 kcal/mol) with the F3L in the MPXV. All-atom simulations of the top-ranked hits 1 to 4 consistently exhibited stable dynamics, suggesting their potential to interact effectively with interface residues. This was further substantiated through analyses of parameters such as radius of gyration (Rg), Root Mean Square Fluctuation, and hydrogen bonding. Cumulative assessments of binding free energy confirmed the top-performing candidates among all the compounds, with values of -35.90, -52.74, -28.17, and -32.11 kcal/mol for top hits 1-4, respectively. These results indicate that compounds top hit 1-4 could hold significant promise for advancing innovative drug therapies, suggesting their suitability for both in vivo and in vitro experiments.

Binding Free Energy Calculation Based on the Fragment Molecular Orbital Method and Its Application in Designing Novel SHP-2 Allosteric Inhibitors

The accurate prediction of binding free energy is a major challenge in structure-based drug design. Quantum mechanics (QM)-based approaches show promising potential in predicting ligand-protein binding affinity by accurately describing the behavior and structure of electrons. However, traditional QM calculations face computational limitations, hindering their practical application in drug design. Nevertheless, the fragment molecular orbital (FMO) method has gained widespread application in drug design due to its ability to reduce computational costs and achieve efficient ab initio QM calculations. Although the FMO method has demonstrated its reliability in calculating the gas phase potential energy, the binding of proteins and ligands also involves other contributing energy terms, such as solvent effects, the 'deformation energy' of a ligand's bioactive conformations, and entropy. Particularly in cases involving ionized fragments, the calculation of solvation free energy becomes particularly crucial. We conducted an evaluation of some previously reported implicit solvent methods on the same data set to assess their potential for improving the performance of the FMO method. Herein, we develop a new QM-based binding free energy calculation method called FMOScore, which enhances the performance of the FMO method. The FMOScore method incorporates linear fitting of various terms, including gas-phase potential energy, deformation energy, and solvation free energy. Compared to other widely used traditional prediction methods such as FEP+, MM/PBSA, MM/GBSA, and Autodock vina, FMOScore showed good performance in prediction accuracies. By constructing a retrospective case study, it was observed that incorporating calculations for solvation free energy and deformation energy can further enhance the precision of FMO predictions for binding affinity. Furthermore, using FMOScore-guided lead optimization against Src homology-2-containing protein tyrosine phosphatase 2 (SHP-2), we discovered a novel and potent allosteric SHP-2 inhibitor (compound 8).

Fascin - F-actin interaction studied by molecular dynamics simulation and protein network analysis

Actin bundles are an important component of cellular cytoskeleton and participate in the movement of cells. The formation of actin bundles requires the participation of many actin binding proteins (ABPs). Fascin is a member of ABPs, which plays a key role in bundling filamentous actin (F-actin) to bundles. However, the detailed interactions between fascin and F-actin are unclear. In this study, we construct an atomic-level structure of fascin - F-actin complex based on a rather poor cryo-EM data with resolution of 20 nm. We first optimized the geometries of the complex by molecular dynamics (MD) simulation and analyzed the binding site and pose of fascin which bundles two F-actin chains. Next, binding free energy of fascin was calculated by MM/GBSA method. Finally, protein structure network analysis (PSNs) was performed to analyze the key residues for fascin binding. Our results show that residues of K22, E27, E29, K41, K43, R110, R149, K358, R408 and K471 on fascin are important for its bundling, which are in good agreement with the experimental data. On the other hand, the consistent results indicate that the atomic-level model of fascin - F-actin complex is reliable. In short, this model can be used to understand the detailed interactions between fascin and F-actin, and to develop novel potential drugs targeting fascin.Communicated by Ramaswamy H. Sarma.

Structure-based virtual screening and molecular dynamics simulations for detecting novel candidates for allosteric inhibition of EGFRT790M

Structure-based virtual screening (SBVS) was applied to predict lead compounds for the allosteric inhibition of epidermal growth factor receptor (EGFR) by screening the library of chemical compounds prepared from the e-molecules chemical database. The library of chemical compounds consisting of 133,083 ligands was composed by evaluating the chemical and physical properties of e-molecules chemicals. The prepared library was screened by CCDC Gold software in the allosteric binding site of EGFRT790M using the library and virtual screening default parameters to filter out, respectively. The GOLD fitness scores 75 and 80 were selected as threshold values for the library and virtual screening processes, respectively. After the docking study, molecular dynamics simulations (MDS) of the top 25 compounds were built for calculating binding free energies from their MDS trajectories. MM-GBSA binding free energies for the compounds were computed from 20 ns MDS, 50 ns MDS and 200 ns MDS trajectories to filter out the candidates. Following MM-GBSA/MM-PBSA binding free energy calculations, six compounds were detected as the most promising candidates for allosteric inhibition of EGFRT790M. The dynamic behaviors of final compounds inside EGFR T790M were searched using structure stability, binding modes and energy decomposition analysis. Besides, the estimated inhibitors were exposed to docking study and MM-GBSA/MM-PBSA binding free energy calculations inside wild-type EGFR, respectively, to be determined their selectivity towards mutant form. Five of the estimated inhibitors displayed estimated selectivity towards EGFRT790M. Besides the ADMET properties of the estimated inhibitors were predicted by PreAdmet tools.Communicated by Ramaswamy H. Sarma.

Identification and stability analysis of potential ADP-ribose modification sites on vascular endothelial growth factor (VEGF) through molecular dynamics simulation

Post-translational modifications (PTMs) are crucial covalent processes that alter protein properties, achieved through proteolytic cleavage or addition of modifying groups like acetyl, phosphoryl, glycosyl, or methyl to amino acids. ADP-ribosylation is a reversible post-translational modification, where ADP-ribose units are covalently attached to target protein side chains. Vascular endothelial growth factor (VEGF) is a potent angiogenic factor that plays a key role in physiological and pathological conditions. Studies have reported that ADP-ribosylation affects VEGF's ability to bind to VEGF receptors, impacting angiogenesis signalling. However, the specific amino acid undergoing ADP-ribosylation on VEGF remained unknown. To understand the mechanism of ADP-ribose addition to VEGF, an in silico study was designed. The study initially checked for the presence of any conserved motif where ADP-ribosylation could potentially occur and identified the presence of the EIE motif in VEGF, a probable site for ADP-ribosylation for many proteins. Subsequently, the amino acids near this motif were selected and their structural properties were analyzed. Surface-exposed amino acids were chosen, and ADP-ribose was then added to their side chains. The results revealed that the amino acids ASP (67) and GLU (70) underwent glycosidic linkage with ADP-ribose, indicating that they are the most probable modification sites. Subsequently, Molecular dynamic simulation analysis such as RMSD, RMSF, Rg, PCA, and FEL, along with MM-PBSA binding free energy calculations were performed to understand the stability of the VEGF-ADP-ribose complexes. The analysis revealed that amino acid at position 67 (ASP67) is the most probable site for ADP-ribosylation in VEGF.Communicated by Ramaswamy H. Sarma.

A computational investigation of potential plant-based bioactive compounds against drug-resistant Staphylococcus aureus of multiple target proteins

Drug-resistant Staphylococcus aureus (DRSA) poses a significant global health threat, like bacteremia, endocarditis, skin, soft tissue, bone, and joint infections. Nowadays, the resistance against conventional drugs has been a prompt and focused medical concern. The present study aimed to explore the inhibitory potential of plant-based bioactive compounds (PBBCs) against effective target proteins using a computational approach. We retrieved and verified 22 target proteins associated with DRSA and conducted a screening process that involved testing 87 PBBCs. Molecular docking was performed between screened PBBCs and reference drugs with selected target proteins via AutoDock. Subsequently, we filtered the target proteins and top PBBCs based on their binding affinity scores. Furthermore, molecular dynamic simulation was carried out through GROMACS for a duration of 100 ns, and the binding free energy was calculated using the gmx_MMPBSA. The result showed consistent hydrogen bonding interactions among the amino acid residues Ser 149, Arg 151, Thr 165, Thr 216, Glu 239, Ser 240, Ile 14, as well as Asn 18, Gln 19, Lys 45, Thr 46, Tyr 109, with their respective target proteins of the penicillin-binding protein and dihydrofolate reductase complex. Additionally, we assessed the pharmacokinetic properties of screened PBBCs via SwissADME and AdmetSAR. The findings suggest that β-amyrin, oleanolic acid, kaempferol, quercetin, and friedelin have the potential to inhibit the selected target proteins. In future research, both in vitro and in vivo, experiments will be needed to establish these PBBCs as potent antimicrobial drugs for DRSA.Communicated by Ramaswamy H. Sarma.

Expanding the therapeutic arsenal against cancer: a computational investigation of hybrid xanthone derivatives as selective Topoisomerase 2α ATPase inhibitors

The DNA topoisomerase II (topo II) enzyme plays an important role in the replication, recombination, and repair of DNA. Despite their widespread applications in cancer therapy, new, selective, and potent topo II inhibitors with better pharmaceutical profiles are needed to handle drug resistance and severe adverse effects. In this respect, an array of 36 new anticancer compounds was designed based on a Xanthone core tethered to multifunctional Pyridine-amines and Imidazole scaffold via alkyl chain linkers. An integrated in silico approach was used to understand the structural basis and mechanism of inhibition of the hybrid xanthone derivatives. In this study, we established an initial virtual screening workflow based on pharmacophore mapping, docking, and cancer target association to validate the target selection process. Next, a simulation-based docking was conducted along with pharmacokinetic analysis to filter out the five best compounds (7, 10, 25, 27, and 30) having binding energies within the range of -60.45 to -40.97 kcal/mol. The screened compounds were further subjected to molecular dynamics simulation for 200 ns followed by MM-GBSA and ligand properties analysis to assess the stability and binding affinity to hTOP2α. The top-ranking hits 3,7-bis(3-(2-aminopyridin-3-ylhydroxy)propoxy)-1-hydroxy-9H-xanthen-9-one (ligand 7) and 3,8-bis(3-(2-aminopyridin-3-ylhydroxy)propoxy)-1-hydroxy-9H-xanthen-9-one (ligand 25) were found to have no toxicity, optimum pharmacokinetic and, DFT properties and stable intermolecular interactions with the active site of hTopo IIα protein. In conclusion, further in vitro and in vivo experimental validation of the identified lead molecules is warranted for the discovery of new human Topoisomerase 2 alpha inhibitors.Communicated by Ramaswamy H. Sarma.

Improved catalytic performance and molecular insight for lipoxygenase from Enterovibrio norvegicus via directed evolution

Lipoxygenase (LOX) holds significant promise for food and pharmaceutical industries. However, albeit its application has been hampered by low catalytic activity and suboptimal thermostability. To address the drawbacks, a directed evolution strategy was explored to enhance the catalytic activity and thermostability of LOX from Enterovibrio norvegicus (EnLOX) for the first time. After two rounds of error-prone polymerase chain reaction (error-prone PCR) and one generations of sequential DNA shuffling, all of four different mutants showed a significant increase in the specific activity of EnLOX, ranging from 132.07 ± 9.34 to 330.17 ± 18.54 U/mg. Among these mutants, D95E/T99A/A121H/S142N/N444W/S613G (EAHNWG) exhibited the highest specific activity, which was 8.25-fold higher than the wild-type enzyme (WT). Meanwhile, the catalytic efficiency (K cat /K m) of EAHNWG was also improved, which was 13.61 ± 1.67 s-1 μM-1, in comparison to that of WT (4.83 ± 0.38 s-1 μM-1). In addition, mutant EAHNWG had a satisfied thermostability with the t 1/2,50 °C value of 6.44 ± 0.24 h, which was 0.4 h longer than that of the WT. Furthermore, the molecular dynamics simulation and structural analysis demonstrated that the reduction of hydrogen bonds number, the enhancement of hydrophobic interactions in the catalytic pocket, and the improvement of flexibility of the lid domain facilitated structural stability and the strength of substrate binding capacity for improved thermal stability and catalytic efficiency of mutant LOX after directed evolution. Overall, these results could provide the guidance for further enzymatic modification of LOX with high catalytic performance for industrial application.

The stability and dynamics of the Aβ40/Aβ42 interlaced mixed fibrils

The accumulation of fibrillar amyloid-β (Aβ) aggregates in the brain, predominantly comprising 40- and 42-residue amyloid-β (Aβ40 and Aβ42), is a major pathological hallmark of Alzheimer's disease (AD). Aβ40 and Aβ42 naturally coexist in the brain under normal physiological conditions, and their interplay is generally considered to be a critical factor in the progression of AD. In addition to forming homogeneous oligomers and fibrils, Aβ40 and Aβ42 are also reported to co-assemble into hetero-oligomers and interlaced mixed fibrils, as evidenced by solid-state nuclear magnetic resonance spectroscopy (NMR), high molecular weight mass spectrometry and cross-seeding experiments. However, the exact molecular mechanisms underlying these processes remain unclear. In this study, we have used a recently resolved structurally uniform 1:1 mixture of Aβ40/Aβ42 interlaced mixed fibril as a prototype to gain insights into the molecular-level interactions between Aβ40 and Aβ42. We employed fully atomistic molecular dynamics simulation and compared the results with a homogeneous U-shaped Aβ40 fibrillar model. Our simulations using two different force fields provide conclusive evidence that the Aβ40/Aβ42 interlaced mixed fibril is energetically more favorable than the homogeneous Aβ40 fibrillar model. Furthermore, we also show that the increase in stability observed in the mixed model stems primarily from the packing interfaces and the stacking interfaces between C-termini. Our simulation results provide valuable mechanistic insights that are not readily accessible in experiment and could have significant implications for both the pathogenesis of AD and the development of current therapeutic strategies.Communicated by Ramaswamy H. Sarma.

Identification of potential molecular targets and repurposed drugs for tuberculosis using network-based screening approach, molecular docking, and simulation

The spread of drug-resistant strains of tuberculosis has hampered efforts to control the disease worldwide. The Mycobacterium tuberculosis cell wall envelope is dynamic, with complex features that protect it from the host immunological response. As a result, the bacterial cell wall components represent a potential target for drug discovery. Protein-protein interaction networks (PPIN) are critical for understanding disease conditions and identifying precise therapeutic targets. We used a rational theoretical approach by constructing a PPIN with the proteins involved in cell wall biosynthesis. The PPIN was constructed through the STRING database and embB was identified as a key protein by using four topological measures, betweenness, closeness, degree, and eigenvector, in the CytoNCA tool in Cytoscape. The 'Drug repurposing' approach was employed to find suitable inhibitors against embB. We used the Schrödinger suites for molecular docking, molecular dynamics simulation, and binding free energy calculations to validate the binding of protein with the ligand. FDA-approved drugs from the ZINC database and DrugBank were screened against embB (PDB ID: 7BVF) using high-throughput virtual screening, standard precision, and extra precision docking. The drugs were screened based on the XP docking score of the standard drug ethambutol. Accordingly, from the top five hits, azilsartan and dihydroergotamine were selected based on the binding free energy values and were further subjected to Molecular Dynamics Simulation studies for 100 ns. Our study confirms that Azilsartan and Dihydroergotamine form stable complexes with embB and can be used as potential lead molecules based on further in vitro and in vivo experimental validation.Communicated by Ramaswamy H. Sarma.

Bioinformatics insights into CENP-T and CENP-W protein-protein interaction disruptive amino acid substitution in the CENP-T-W complex

Kinetochores are multi-protein assemblies present at the centromere of the human chromosome and play a crucial role in cellular mitosis. The CENP-T and CENP-W chains form a heterodimer, which is an integral part of the inner kinetochore, interacting with the linker DNA on one side and the outer kinetochore on the other. Additionally, the CENP-T-W dimer interacts with other regulatory proteins involved in forming inner kinetochores. The specific roles of different amino acids in the CENP-W at the protein-protein interaction (PPI) interface during the CENP-T-W dimer formation remain incompletely understood. Since cell division goes awry in diseases like cancer, this CENP-T-W partnership is a potential target for new drugs that could restore healthy cell division. We employed molecular docking, binding free energy calculations, and molecular dynamics (MD) simulations to investigate the disruptive effects of amino acids substitutions in the CENP-W chain on CENP-T-W dimer formation. By conducting a molecular docking study and analysing hydrogen bonding interactions, we identified key residues in CENP-W (ASN-46, ARG-53, LEU-83, SER-86, ARG-87, and GLY-88) for further investigation. Through site-directed mutagenesis and subsequent binding free energy calculations, we refined the selection of mutant. We chose four mutants (N46K, R53K, L83K, and R87E) of CENP-W to assess their comparative potential in forming CENP-T-W dimer. Our analysis from 250 ns long revealed that the substitution of LEU83 and ARG53 residues in CENP-W with the LYS significantly disrupts the formation of CENP-T-W dimer. In conclusion, LEU83 and ARG53 play a critical role in CENP-T and CENP-W dimerization which is ultimately required for cellular mitosis. Our findings not only deepen our understanding of cell division but also hint at exciting drug-target possibilities.

In silico discovery of potential sodium-glucose cotransporter-2 inhibitors from natural products for treatment of heart failure via molecular docking and molecular dynamics simulation approach

Heart failure (HF) is the end stage of cardiovascular disease. Because of its complex condition and poor prognosis, HF has become an important public health problem in the world. Sodium-glucose cotransporter-2 (SGLT2) is a member of the glucose transporter family. Recently, SGLT2 inhibitors have been applied to treat HF. In this study, the main aim was to discover natural SGLT2 inhibitor from Chinese herbs through docking-based virtual screening. Totally 113 natural compounds of potential SGLT2 inhibitor were identified, which displayed docking affinity higher than six approved inhibitors (dapagliflozin (IC50 = 4.9 nM), canagliflozin (IC50 = 4.4 nM 6.7), ipragliflozin (IC50 = 7.4 nM), empagliflozin (IC50 = 3.1 nM), tofogliflozin (IC50 = 4 nM) and luseogliflozin (IC50 = 2.3 nM)) through docking-based virtual screening. Then, the top three hits (ZINC70455591, ZINC85594065 and ZINC14588133) and six known inhibitors were selected for molecular dynamics simulation and the binding free energy calculation using molecular mechanics Poisson-Boltzmann surface area to demonstrate the stability and affinity of docked complexes. These results pointed out that the three docked complexes were stabilized and the chosen compounds were tightly adhering to the binding site of SGLT2. Besides, pharmacokinetic properties of the selected compounds showed those natural compounds may be potential drug candidates. This study may be contributed to further in vitro and in vivo validation and the development of novel SGLT2 inhibitor for treating HF.Communicated by Ramaswamy H. Sarma.

Differences of Atomic-Level Interactions between Midazolam and Two CYP Isoforms 3A4 and 3A5

CYP 3A4 and CYP 3A5 are two important members of the human cytochrome P450 family. Although their overall structures are similar, the local structures of the active site are different, which directly leads to obvious individual differences in drug metabolic efficacy and toxicity. In this work, midazolam (MDZ) was selected as the probe substrate, and its interaction with two proteins, CYP 3A4 and CYP 3A5, was studied by molecular dynamics simulation (MD) along with the calculation of the binding free energy. The results show that two protein-substrate complexes have some similarities in enzyme-substrate binding; that is, in both complexes, Ser119 forms a high occupancy hydrogen bond with MDZ, which plays a key role in the stability of the interaction between MDZ and the enzymes. However, the complex formed by CYP 3A4 and MDZ is more stable, which may be attributed to the sandwich structure formed by the fluorophenyl group of the substrate with Leu216 and Leu482. Our study interprets the binding differences between two isoform-substrate complexes and reveals a structure-function relationship from the atomic perspective, which is expected to provide a theoretical basis for accurately measuring the effectiveness and toxicity of drugs for individuals in the era of precision medicine.

A multi-tier computational screening framework to effectively search the mutational space of SARS-CoV-2 receptor binding motif to identify mutants with enhanced ACE2 binding abilities

SARS-CoV-2 gained crucial mutations at the receptor binding domain (RBD) that often changed the course of the pandemic leading to new waves with increased case fatality. Variants are observed with enhanced transmission and immune invasion abilities. Thus, predicting future variants with enhanced transmission ability is a problem of utmost research interest. Here, we have developed a multi-tier exhaustive SARS-CoV-2 mutation screening platform combining MM/GBSA, extensive molecular dynamics simulations, and steered molecular dynamics to identify RBD mutants with enhanced ACE2 binding capability. We have identified four RBM mutations (F490K, S494K, G504F, and the P499L) with significantly higher ACE2 binding abilities than wild-type RBD. Compared to wild-type RBD, they all form stable complexes with more hydrogen bonds and salt-bridge interactions with ACE2. Our simulation data suggest that these mutations allosterically alter the packing of the RBM interface of the RBD-ACE2 complex. As a result, the rupture force required to break the RBD-ACE2 contacts is significantly higher for these mutants.

The binding mechanism of failed, in processing and succeed inhibitors target SARS-CoV-2 main protease

Since the outbreak of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), several variants have caused a persistent pandemic. Consequently, it is crucial to develop new potential anti-SARS-CoV-2 drugs with specificity. To minimize potential failures and preserve valuable clinical resources for the development of other useful drugs, researchers must enhance their understanding of the interactions between drugs and SARS-CoV-2. While numerous crystal structures of the SARS-CoV-2 main protease (SCM) and its inhibitors have been reported, they provide only static snapshots and fail to capture the dynamic nature of SCM/inhibitor interactions. Herein, we conducted molecular dynamics simulations for five SCM complexes: ritonavir (SCM/RTV), lopinavir (SCM/LPV), the identified inhibitor N3 (SCM/N3), the approved inhibitor ensitrelvir (SCM/ESV), and the approved drug nirmatrelvir (SCM/NMV). Additionally, we explored the potential for covalent bond formation in the N3 and NMV inhibitors through QM/MM calculations using Umbrella sampling. The results show that the binding site is highly flexible to fit those five different inhibitors and each compound has its unique binding mode at the same binding site. Moreover, the binding affinities of positive and negative inhibitors to SCM exhibit significant differences. By gaining insights into the dynamics, we can potentially elucidate why lopinavir/ritonavir, initially considered promising, failed to effectively treat COVID-19. Furthermore, understanding the mechanistic aspects of N3 and NMV inhibition on SCM not only contributes to rational drug discovery against COVID-19 but also aids future studies on the catalytic mechanisms of main proteases in other novel coronaviruses.Communicated by Ramaswamy H. Sarma.

Mutational analyses, pharmacophore-based inhibitor design and in silico validation for Zika virus NS3-helicase

Zika virus is responsible for causing Zika infections and was declared as a public health emergency of international concern in February 2016. The Zika virus NS3-helicase is a viable drug target for the design of inhibitors due to its essential role in the replication of viral genome. The viral RNA is unwound by the NS3-helicase in order to enable the reproduction of viral genome by the NS5 protein. Zika virus infections in humans are being reported for the last 15 years. We have therefore carried out amino acid mutational analyses of NS3-helicase. NS3-helicase has two crucial binding sites: the ATP and RNA binding sites. The cofactor-ATP based pharmacophore was generated for virtual screening of ZINC database using Pharmit server, that is followed by molecular docking and molecular dynamics simulations of potential hits as probable Zika virus NS3-helicase inhibitors at the cofactor binding site. The drug-like properties of the molecules were analysed and, DFT calculations were performed on the five best molecules to reveal their stability in solvent phase compared to gas phase, the HOMO and LUMO and electrostatic potential maps to analyze the electronic and geometric characteristics. These are significant findings towards the discovery of new inhibitors of Zika virus NS3-helicase, a promising drug target to treat the Zika virus infection.Communicated by Ramaswamy H. Sarma.

Docking and Molecular Dynamics Simulations Clarify Binding Sites for Interactions of Novel Marine Sulfated Glycans with SARS-CoV-2 Spike Glycoprotein

The entry of SARS-CoV-2 into the host cell is mediated by its S-glycoprotein (SGP). Sulfated glycans bind to the SGP receptor-binding domain (RBD), which forms a ternary complex with its receptor angiotensin converting enzyme 2. Here, we have conducted a thorough and systematic computational study of the binding of four oligosaccharide building blocks from novel marine sulfated glycans (isolated from Pentacta pygmaea and Isostichopus badionotus) to the non-glycosylated and glycosylated RBD. Blind docking studies using three docking programs identified five potential cryptic binding sites. Extensive site-targeted docking and molecular dynamics simulations using two force fields confirmed only two binding sites (Sites 1 and 5) for these novel, highly charged sulfated glycans, which were also confirmed by previously published reports. This work showed the structural features and key interactions driving ligand binding. A previous study predicted Site 2 to be a potential binding site, which was not observed here. The use of several molecular modeling approaches gave a comprehensive assessment. The detailed comparative study utilizing multiple modeling approaches is the first of its kind for novel glycan-SGP interaction characterization. This study provided insights into the key structural features of these novel glycans as they are considered for development as potential therapeutics.

Designing new drug candidates as inhibitors against wild and mutant type neuraminidases: molecular docking, molecular dynamics and binding free energy calculations

Influenza virus is the cause of the death of millions of people with about 3-4 pandemics every hundred years in history. It also turns into a seasonal disease, bringing about approximately 5-15% of the population to be infected and 290,000-650,000 people to die every year. These numbers reveal that it is necessary to be on the alert to work towards influenza in order to protect public health. There are FDA-approved antiviral drugs such as oseltamivir and zanamivir recommended by the World Center for Disease Prevention. However, after the recent outbreaks such as bird flu and swine flu, increasing studies have shown that the flu virus has gained resistance to these drugs. So, there is an urgent need to find new drugs effective against this virus. This study aims to investigate new drug candidates targeting neuraminidase (NA) for the treatment of influenza by using computer aided drug design approaches. They involve virtual scanning, de novo design, rational design, docking, MD, MMGB/PBSA. The investigation includes H1N1, H5N1, H2N2 and H3N2 neuraminidase proteins and their mutant variants possessing resistance to FDA-approved drugs. Virtual screening consists of approximately 30 thousand molecules while de novo and rational designs produced over a hundred molecules. These approaches produced three lead molecules with binding energies for both non-mutant (-34.84, -59.99 and -60.66 kcal/mol) and mutant (-40.40, -58.93, -76.19 kcal/mol) H2N2 NA calculated by MM-PBSA compared with those of oseltamivir -25.64 and -18.40 respectively. The results offer new drug candidates against influenza infection.Communicated by Ramaswamy H. Sarma.

Binding of Mammea A/AA (MA) to calf thymus DNA revealed by the ratiometric absorbance of MA in the UV-visible range molecular dynamic simulations and TD-DFT calculations

The in vitro anti-proliferative activity of MA (5,7-dihydroxy-8-(3-methylbut-2-enyl)-6-(3-methyl-1-oxobutyl)-4-phenyl[1]2H-[1]benzopyran-2-one)on a variety of cancer cells was previously demonstrated. This work strives to understand the mechanisms by which MA exerts this biological activity. Thereafter, the binding of MA to calf thymus DNA was studied by monitoring the change in the UV-visible absorbance of MA. It was found that, the response of MA to binding with calf thymus DNA is characterised by an increase in the AS/AL ratio of the absorbance of the longest wavelength absorption band to the shortest one, and the appearance of a new band at about 377 nm assigned to S0→S1 transition, which is red shifted as compared to free MA. From the bands ratio, the binding constant is found to be 4.3x105 M-1, indicating strong binding. The deduced binding free energy, enthalpy and entropy are -7.7 kcal/mol, -10.89 ± 0.28 kcal/mol and -54.46 ± 4 J/K, respectively, indicating that MA binds to DNA by a non-bonding Van der Waals type interactions and hydrogen bonds. Further study with classical molecular dynamics shows that MA binds to DNA by intercalation, where it is positioned between two AT base pairs. Unlike isolated MA, TDDFT calculations on ten images extracted from the MD trajectory show that, the frontier molecular orbitals of the complex are distributed over the DNA and MA. This indicates a strong stacking interaction and then explains the hypochromism and the red shift of the S0→S1 transition. The present work demonstrates the potency of MA as antitumor compound and as absorbance-based molecular probe.Communicated by Ramaswamy H. Sarma.

Structure based exploration of potential lead molecules against the extracellular cysteine protease (EcpA) of Staphylococcus epidermidis: a therapeutic halt

Nosocomial infection caused by Staphylococcus epidermidis is one of the most widely spread diseases affecting the world's population. No strategies have been developed to overcome this infection and inhibit its spread in immunocompromised patients or patients with indwelling medical devices. EcpA is an extracellular cysteine protease protein involved in biofilm formation on medical devices. Thus, blocking this mechanism may be viable for developing a drug against S. epidermidis. The current research aimed to find new, potent inhibitors that could stop the S. epidermidis EcpA protein from functioning. This study attempted to identify the most promising drug candidates using structure-based virtual screening (SBVS) from libraries of natural ligands. The top-scored molecules were shortlisted based on their IC50 values and pharmacophore properties and further validated through density functional theory (DFT) studies. We found five inhibitors using virtual screening, and the results indicate that these drugs had the highest energy binding potential towards the EcpA targets when compared to the reference molecule E-64, a known cysteine protease inhibitor. In order to evaluate the binding conformational stability of protein-ligand complexes, molecular dynamics (MD) simulations were performed in triplicate for 100 ns, revealing the significant stability of anticipated molecules at the docked site. Furthermore, principal component analysis and binding free energy calculations were performed to understand the dynamics and stability of the complexes. The current study indicated that these compounds looked to be suitable novel inhibitors of the EcpA protein and pave the path for further discovery of novel inhibitors of EcpA.Communicated by Ramaswamy H. Sarma.

Exploring the natural products chemical space through a molecular search to discover potential inhibitors that target the hypoxia-inducible factor (HIF) prolyl hydroxylase domain (PHD)

Introduction: Hypoxia-inducible factor (HIF) prolyl hydroxylase domain (PHD) enzymes are major therapeutic targets of anemia and ischemic/hypoxia diseases. To overcome safety issues, liver failure, and problems associated with on-/off-targets, natural products due to their novel and unique structures offer promising alternatives as drug targets. Methods: In the current study, the Marine Natural Products, North African, South African, East African, and North-East African chemical space was explored for HIF-PHD inhibitors discovery through molecular search, and the final hits were validated using molecular simulation and free energy calculation approaches. Results: Our results revealed that CMNPD13808 with a docking score of -8.690 kcal/mol, CID15081178 with a docking score of -8.027 kcal/mol, CID71496944 with a docking score of -8.48 kcal/mol and CID11821407 with a docking score of -7.78 kcal/mol possess stronger activity than the control N-[(4-hydroxy-8-iodoisoquinolin-3-yl)carbonyl]glycine, 4HG (-6.87 kcal/mol). Interaction analysis revealed that the target compounds interact with Gln239, Tyr310, Tyr329, Arg383 and Trp389 residues, and chelate the active site iron in a bidentate manner in PHD2. Molecular simulation revealed that these target hits robustly block the PHD2 active site by demonstrating stable dynamics. Furthermore, the half-life of the Arg383 hydrogen bond with the target ligands, which is an important factor for PHD2 inhibition, remained almost constant in all the complexes during the simulation. Finally, the total binding free energy of each complex was calculated as CMNPD13808-PHD2 -72.91 kcal/mol, CID15081178-PHD2 -65.55 kcal/mol, CID71496944-PHD2 -68.47 kcal/mol, and CID11821407-PHD2 -62.06 kcal/mol, respectively. Conclusion: The results show the compounds possess good activity in contrast to the control drug (4HG) and need further in vitro and in vivo validation for possible usage as potential drugs against HIF-PHD2-associated diseases.

Computational investigation on the binding modes of PET polymer to PETase

Poly(ethylene terephthalate) (PET) has been widely utilized in daily life, but its non-degradability has induced severe environmental and health problems. Recently, PETase, which has been isolated from bacterium Ideonella sakaiensisis, was reported to have the highest PET degradation activity and specificity under room temperature, but no crystal structure for PET in complex with PETase has been reported. To provide deep insight into the binding mode of PET polymer on PETase and the binding interactions, we employed molecular docking and molecular dynamics simulations to study the substrate binding at the atomic level. Different PET oligomers have been studied with chain lengths varying from 2 to 8. In addition, the binding energies and hot-spot residues were analyzed to gain better insights into the binding mechanism by MM/GBSA approach. The PET oligomers adopt stable and reactive conformations in a shallow cleft on a flat surface of PETase. The binding cleft can only accommodate four moieties, and others beyond the region will be stabilized by the π-stacking interactions with Trp156 at the terephthalic acid terminal. Our studies provide a clear picture of how the binding mode of PET polymer and its interactions with PETase change with the chain length. Those studies would provide useful information for the rational design of catalytically more efficient PETase variants toward plastic degradation.Communicated by Ramaswamy H. Sarma.

Binding Affinity of Trastuzumab and Pertuzumab Monoclonal Antibodies to Extracellular HER2 Domain

The binding affinity of trastuzumab and pertuzumab to HER2 has been studied using both experimental and in silico methods. The experiments were conducted using the antibodies in their complete IgG form, as used in clinical therapy, and the extracellular domain of the HER2 protein in solution. This approach provides a precise, reproducible, and reliable view of the interaction between them in physicochemical conditions similar to those found in the tumoral environment. Dynamic light scattering and size exclusion chromatography coupled with tetra detection were utilized to characterize the protein complexes, measure their concentrations, and calculate the equilibrium-free binding energy, ΔGbind. In addition, PRODIGY, a QSAR-like model with excellent predictive ability, was employed to obtain in silico ΔGbind estimations. The results obtained indicate that pertuzumab exhibits a slightly higher binding affinity to HER2 than trastuzumab. The difference in binding affinity was explained based on the contribution of the different interfacial contact (IC) descriptors to the ΔGbind value estimated by the PRODIGY model. Furthermore, experiments revealed that the pertuzumab IgG antibody binds preferentially to two HER2 proteins, one per Fab fragment, while trastuzumab mainly forms a monovalent complex. This finding was interpreted based on a geometrical model that identified steric crowding in the trastuzumab-HER2 complex as compared with the pertuzumab-HER2 complex.

Theoretical and Experimental Studies on Plant Light-Dependent Protochlorophyllide Oxidoreductase as a Novel Target for Searching Potential Herbicides

Herbicide resistance is a prevalent problem that has posed a foremost challenge to crop production worldwide. Light-dependent enzyme NADPH: protochlorophyllide oxidoreductase (LPOR) in plants is a metabolic target that could satisfy this unmet demand. Herein, for the first time, we embarked on proposing a new mode of action of herbicides by performing structure-based virtual screening targeting multiple LPOR binding sites, with the determination of further bioactivity on the lead series. The feasibility of exploiting high selectivity and safety herbicides targeting LPOR was discussed from the perspective of the origin and phylogeny. Besides, we revealed the structural rearrangement and the selection key for NADPH cofactor binding to LPOR. Based on these, multitarget virtual screening was performed and the result identified compounds 2 affording micromolar inhibition, in which the IC₅₀ reached 4.74 μM. Transcriptome analysis revealed that compound 2 induced more genes related to chlorophyll synthesis in Arabidopsis thaliana, especially the LPOR genes. Additionally, we clarified that these compounds binding to the site enhanced the overall stability and local rigidity of the complex systems from molecular dynamics simulation. This study delivers a guideline on how to assess activity-determining features of inhibitors to LPOR and how to translate this knowledge into the design of novel and effective inhibitors against malignant weed that act by targeting LPOR.

Application of Ligand- and Structure-Based Prediction Models for the Design of Alkylhydrazide-Based HDAC3 Inhibitors as Novel Anti-Cancer Compounds

Histone deacetylases (HDAC) represent promising epigenetic targets for several diseases including different cancer types. The HDAC inhibitors approved to date are pan-HDAC inhibitors and most show a poor selectivity profile, side effects, and in particular hydroxamic-acid-based inhibitors lack good pharmacokinetic profiles. Therefore, the development of isoform-selective non-hydroxamic acid HDAC inhibitors is a highly regarded field in medicinal chemistry. In this study, we analyzed different ligand-based and structure-based drug design techniques to predict the binding mode and inhibitory activity of recently developed alkylhydrazide HDAC inhibitors. Alkylhydrazides have recently attracted more attention as they have shown promising effects in various cancer cell lines. In this work, pharmacophore models and atom-based quantitative structure-activity relationship (QSAR) models were generated and evaluated. The binding mode of the studied compounds was determined using molecular docking as well as molecular dynamics simulations and compared with known crystal structures. Calculated free energies of binding were also considered to generate QSAR models. The created models show a good explanation of in vitro data and were used to develop novel HDAC3 inhibitors.

Exploring the Traditional Chinese Medicine (TCM) database chemical space to target I7L protease from monkeypox virus using molecular screening and simulation approaches

In the current study, we used molecular screening and simulation approaches to target I7L protease from monkeypox virus (mpox) from the Traditional Chinese Medicines (TCM) database. Using molecular screening, only four hits TCM27763, TCM33057, TCM34450 and TCM31564 demonstrated better pharmacological potential than TTP6171 (control). Binding of these molecules targeted Trp168, Asn171, Arg196, Cys237, Ser240, Trp242, Glu325, Ser326, and Cys328 residues and may affect the function of I7L protease in in vitro assay. Moreover, molecular simulation revealed stable dynamics, tighter structural packing and less flexible behaviour for all the complexes. We further reported that the average hydrogen bonds in TCM27763, TCM33057, TCM34450 and TCM31564I7L complexes remained higher than the control drug. Finally, the BF energy results revealed -62.60 ± 0.65 for the controlI7L complex, for the TCM27763I7L complex -71.92 ± 0.70 kcal/mol, for the TCM33057I7L complex the BF energy was -70.94 ± 0.70 kcal/mol, for the TCM34450I7L the BF energy was -69.94 ± 0.85 kcal/mol while for the TCM31564I7L complex the BF energy was calculated to be -69.16 ± 0.80 kcal/mol. Although, we used stateoftheart computational methods, these are theoretical insights that need further experimental validation.

Building polarization into protein-inhibitor binding dynamics in rational drug design for rheumatoid arthritis

Standard force field-based simulations to accomplish structure-based evaluations of lead molecules is a powerful tool. Combining protein fragmentation into tractable sub-systems with continuum solvation method is envisaged to enable quantum mechanics-based electronic structure calculations of macromolecules in their realistic environment. This along with incorporation of many-body polarization effect in molecular dynamics simulations may augment an accurate description of electrostatics of protein-inhibitor systems for effective drug design. Rheumatoid arthritis (RA) is a complex autoimmune disorder plagued by the ceiling effect of current targeted therapies, encouraging identification of new druggable targets and corresponding drug design to tackle the refractory form of disease. In this study, polarization-inclusive force field approach has been used to model protein solvation and ligand binding for 'Mitogen-activated protein kinase' (MAP3K8), a regulatory node of notable pharmacological relevance in RA synovial biology. For MAP3K8 inhibitors belonging to different scaffold series, the calculations illustrated differential electrostatic contribution to their relative binding affinities and successfully explained examples from available structure-activity relationship studies. Results from this study exemplified i) the advantage of this approach in reliably ranking inhibitors having close nanomolar range activities for the same target; and ii) its prospective application in lead molecule identification aiding drug discovery efforts in RA.Communicated by Ramaswamy H. Sarma.

Binding Affinity and Mechanisms of Potential Antidepressants Targeting Human NMDA Receptors

Depression, a mental disorder that plagues the world, is a burden on many families. There is a great need for new, fast-acting antidepressants to be developed. N-methyl-D-aspartic acid (NMDA) is an ionotropic glutamate receptor that plays an important role in learning and memory processes and its TMD region is considered as a potential target to treat depression. However, due to the unclear binding sites and pathways, the mechanism of drug binding lacks basic explanation, which brings great complexity to the development of new drugs. In this study, we investigated the binding affinity and mechanisms of an FDA-approved antidepressant (S-ketamine) and seven potential antidepressants (R-ketamine, memantine, lanicemine, dextromethorphan, Ro 25-6981, ifenprodil, and traxoprodil) targeting the NMDA receptor by ligand-protein docking and molecular dynamics simulations. The results indicated that Ro 25-6981 has the strongest binding affinity to the TMD region of the NMDA receptor among the eight selected drugs, suggesting its potential effective inhibitory effect. We also calculated the critical binding-site residues at the active site and found that residues Leu124 and Met63 contributed the most to the binding energy by decomposing the free energy contributions on a per-residue basis. We further compared S-ketamine and its chiral molecule, R-ketamine, and found that R-ketamine had a stronger binding capacity to the NMDA receptor. This study provides a computational reference for the treatment of depression targeting NMDA receptors, and the proposed results will provide potential strategies for further antidepressant development and is a useful resource for the future discovery of fast-acting antidepressant candidates.

Potential Binding Sites of Pharmacological Chaperone NCGC00241607 on Mutant β-Glucocerebrosidase and Its Efficacy on Patient-Derived Cell Cultures in Gaucher and Parkinson's Disease

Mutations in the GBA1 gene, encoding the lysosomal enzyme glucocerebrosidase (GCase), cause Gaucher disease (GD) and are the most common genetic risk factor for Parkinson's disease (PD). Pharmacological chaperones (PCs) are being developed as an alternative treatment approach for GD and PD. To date, NCGC00241607 (NCGC607) is one of the most promising PCs. Using molecular docking and molecular dynamics simulation we identified and characterized six allosteric binding sites on the GCase surface suitable for PCs. Two sites were energetically more preferable for NCGC607 and located nearby to the active site of the enzyme. We evaluated the effects of NCGC607 treatment on GCase activity and protein levels, glycolipids concentration in cultured macrophages from GD (n = 9) and GBA-PD (n = 5) patients as well as in induced human pluripotent stem cells (iPSC)-derived dopaminergic (DA) neurons from GBA-PD patient. The results showed that NCGC607 treatment increased GCase activity (by 1.3-fold) and protein levels (by 1.5-fold), decreased glycolipids concentration (by 4.0-fold) in cultured macrophages derived from GD patients and also enhanced GCase activity (by 1.5-fold) in cultured macrophages derived from GBA-PD patients with N370S mutation (p < 0.05). In iPSC-derived DA neurons from GBA-PD patients with N370S mutation NCGC607 treatment increased GCase activity and protein levels by 1.1-fold and 1.7-fold (p < 0.05). Thus, our results showed that NCGC607 could bind to allosteric sites on the GCase surface and confirmed its efficacy on cultured macrophages from GD and GBA-PD patients as well as on iPSC-derived DA neurons from GBA-PD patients.

Identification of abemaciclib derivatives targeting cyclin-dependent kinase 4 and 6 using molecular dynamics, binding free energy calculation, synthesis, and pharmacological evaluation

CDK4/6 plays a crucial role in various cancers and is an effective anticancer drug target. However, the gap between clinical requirements and approved CDK4/6 drugs is unresolved. Thus, there is an urgent need to develop selective and oral CDK4/6 inhibitors, particularly for monotherapy. Here, we studied the interaction between abemaciclib and human CDK6 using molecular dynamics simulations, binding free energy calculations, and energy decomposition. V101 and H100 formed stable hydrogen bonds with the amine-pyrimidine group, and K43 interacted with the imidazole ring via an unstable hydrogen bond. Meanwhile, I19, V27, A41, and L152 interacted with abemaciclib through π-alkyl interactions. Based on the binding model, abemaciclib was divided into four regions. With one region modification, 43 compounds were designed and evaluated using molecular docking. From each region, three favorable groups were selected and combined with each other to obtain 81 compounds. Among them, C2231-A, which was obtained by removing the methylene group from C2231, showed better inhibition than C2231. Kinase profiling revealed that C2231-A showed inhibitory activity similar to that of abemaciclib; additionally, C2231-A inhibited the growth of MDA-MB-231 cells to a greater extent than did abemaciclib. Based on molecular dynamics simulation, C2231-A was identified as a promising candidate compound with considerable inhibitory effects on human breast cancer cell lines.

In Silico Optimization of SARS-CoV-2 Spike Specific Nanobodies

BACKGROUND

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread worldwide, caused a global pandemic, and killed millions of people. The spike protein embedded in the viral membrane is essential for recognizing human receptors and invading host cells. Many nanobodies have been designed to block the interaction between spike and other proteins. However, the constantly emerging viral variants limit the effectiveness of these therapeutic nanobodies. Therefore, it is necessary to find a prospective antibody designing and optimization approach to deal with existing or future viral variants.

METHODS

We attempted to optimize nanobody sequences based on the understanding of molecular details by using computational approaches. First, we employed a coarse-grained (CG) model to learn the energetic mechanism of the spike protein activation. Next, we analyzed the binding modes of several representative nanobodies with the spike protein and identified the key residues on their interfaces. Then, we performed saturated mutagenesis of these key residue sites and employed the CG model to calculate the binding energies.

RESULTS

Based on analysis of the folding energy of the angiotensin-converting enzyme 2 (ACE2) -spike complex, we constructed a detailed free energy profile of the activation process of the spike protein which provided a clear mechanistic explanation. In addition, by analyzing the results of binding free energy changes following mutations, we determined how the mutations can improve the complementarity with the nanobodies on spike protein. Then we chose 7KSG nanobody as a template for further optimization and designed four potent nanobodies. Finally, based on the results of the single-site saturated mutagenesis in complementarity determining regions (CDRs), combinations of mutations were performed. We designed four novel, potent nanobodies, all exhibiting higher binding affinity to the spike protein than the original ones.

CONCLUSIONS

These results provide a molecular basis for the interactions between spike protein and antibodies and promote the development of new specific neutralizing nanobodies.

A Computational Biology Study on the Structure and Dynamics Determinants of Thermal Stability of the Chitosanase from Aspergillus fumigatus

Environmentally friendly and efficient biodegradation with chitosanase for degrading chitosan to oligosaccharide has been gaining more importance. Here, we studied a chitosanase from Aspergillus fumigatus with potential for production, but does not have the ideal thermal stability. The structure predicted by the Alphafold2 model, especially the binding site and two catalytic residues, has been found to have a high similarity with the experimental structure of the chitosanase V-CSN from the same family. The effects of temperature on structure and function were studied by dynamic simulation and the results showed that the binding site had high flexibility. After heating up from 300 K to 350 K, the RMSD and RMSF of the binding site increased significantly, in particular, the downward shift of loop6 closed the binding site, resulting in the spatial hindrance of binding. The time proportions of important hydrogen bonds at the binding site decreased sharply, indicating that serious disruption of hydrogen bonds should be the main interaction factor for conformational changes. The residues contributing energetically to binding were also revealed to be in the highly flexible region, which inevitably leads to the decrease in the activity stability at high temperature. These findings provide directions for the modification of thermal stability and perspectives on the research of proteins without experimental structures.

Study of novel androgen receptor V770 variant in androgen insensitivity syndrome patients reveals the transitional state of the androgen receptor ligand binding domain homodimer

We report the discovery of the androgen receptor missense mutation V770D, that was found in two sisters suffering from complete androgen insensitivity. Experimental validation of AR V770 variants demonstrated that AR V770D was transcriptionally inactive due to the inability to dimerize and a reduced ligand binding affinity. The more conservative AR V770A variant showed a dimerization defect at low levels of DHT with a partial recovery of the transcriptional activity and of the receptor's ability to dimerize when increasing the DHT levels. With V770 located outside of the proposed LBD dimerization interface of the AR LBD homodimer crystal structure, the effects of the V770A mutation on AR dimerization were unexpected. We therefore explored whether the AR LBD dimerization interface would be better described by an alternative dimerization mode based on available human homodimeric LBD crystal structures of other nuclear receptors. Superposition of the monomeric AR LBD in the homodimeric crystal structures of GR, PR, ER, CAR, TRβ, and HNF-4α showed that the GR-like LBD dimer model was energetically the most stable. Moreover, V770 was a key energy residue in the GR-like LBD dimer while it was not involved in the stabilization of the AR LBD homodimer according to the crystal structure. Additionally, the observation that 4 AIS mutations impacted the stability of the AR LBD dimer while 16 mutations affected the GR-like LBD dimer, suggested that the AR LBD dimer crystal is a snapshot of a breathing AR LBD homodimer that can transition into the GR-like LBD dimer model.

Calculations of the binding free energies of the Comprehensive in vitro Proarrhythmia Assay (CiPA) reference drugs to cardiac ion channels

The evaluation of the inhibitory activities of drugs on multiple cardiac ion channels is required for the accurate assessment of proarrhythmic risks. Moreover, the in silico prediction of such inhibitory activities of drugs on cardiac channels can improve the efficiency of the drug-development process. Here, we performed molecular docking simulations to predict the complex structures of 25 reference drugs that were proposed by the Comprehensive in vitro Proarrhythmia Assay consortium using two cardiac ion channels, the human ether-a-go-go-related gene (hERG) potassium channel and human NaV1.5 (hNaV1.5) sodium channel, with experimentally available structures. The absolute binding free energy (ΔGbind) values of the predicted structures were calculated by a molecular dynamics-based method and compared with the experimental half-maximal inhibitory concentration (IC50) data. Furthermore, the regression analysis between the calculated values and negative of the common logarithm of the experimental IC50 values (pIC50) revealed that the calculated values of four and ten drugs deviated significantly from the regression lines of the hERG and hNaV1.5 channels, respectively. We reconsidered the docking poses and protonation states of the drugs based on the experimental data and recalculated their ΔGbind values. Finally, the calculated ΔGbind values of 24 and 19 drugs correlated with their experimental pIC50 values (coefficients of determination=0.791 and 0.613 for the hERG and hNaV1.5 channels, respectively). Thus, the regression analysis between the calculated ΔGbind and experimental IC50 data ensured the realization of an increased number of reliable complex structures.

Multidimensional virtual screening approaches combined with drug repurposing to identify potential covalent inhibitors of SARS-CoV-2 3CL protease

The outbreak of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused an unprecedented global pandemic, and new cases are still on the rise due to the absence of effective medicines. However, developing new drugs within a short time is extremely difficult. Repurposing the existing drugs provides a fast and effective strategy to identify promising inhibitors. Here we focus on the screening of drugs database for discovering potential covalent inhibitors that target 3-chymotrypsin-like protease (3CLpro), an essential enzyme mediating viral replication and transcription. Firstly, we constructed a receptor-ligand pharmacophore model and verified it through decoy set. The importance of pharmacophore features was evaluated by combining molecular dynamics simulation with interaction analyses. Then, covalent docking was used to perform further screening. According to docking score and Prime/Molecular Mechanics Generalized Born Surface Area (MM-GBSA) score, total ten compounds obtained good scores and successfully established covalent bonds with the catalytic Cys145 residue. They also formed favorable interactions with key residues in active sites and closely integrated with 3CLpro with binding modes similar to known 3CLpro inhibitor. Finally, the top four hits DB08732, DB04653, DB01871 and DB07299 were further subjected to 100 ns molecular dynamics (MD) simulation and MM-GBSA binding free energy calculations. The results suggest that the four candidates show good binding affinities for 3CLpro, which warrants further evaluation for their in-vitro/in-vivo activities. Overall, our research methods provide a valuable reference for discovering promising inhibitors against SARS-CoV-2 and help to fight against the epidemic.Communicated by Ramaswamy H. Sarma.