Binding of Inhibitors to BACE1 Affected by pH-Dependent Protonation: An Exploration from Multiple Replica Gaussian Accelerated Molecular Dynamics and MM-GBSA Calculations

Rational design, docking, simulation, synthesis, and in vitro studies of small benzothiazole molecules as selective BACE1 inhibitors

BACE-1 is an encouraging target for the development of AD therapeutics. However, many BACE-1 inhibitors failed clinical trials due to their non-selectivity towards BACE-2 or adverse effects. Herein, a set of 96 benzothiazoles were designed based on the structural features of Atabecestat and Riluzole to find a promising selective BACE-1 inhibitor. Out of the 96 designed compounds, compound 72 showed comparable binding affinity with BACE-1 as compared to Atabecestat, and more selective towards BACE-1 as compared to BACE-2. The BACE-1 docking score of Atabecestat and compound 72 were found to be -7.76 and -7.49, respectively, while their corresponding MM-GBSA ΔGbind energy were -70.39 and -68.97 kcal/mol. In contrast, the BACE-2 docking score of Atabecestat and compound 72 were found to be -6.24 and -5.32, respectively, while their corresponding MM-GBSA ΔGbind energy were -56.02 and -43.46 kcal/mol. The strong binding affinity of compound 72 was further validated by 100 ns dynamics study. The physicochemical and pharmacokinetic (ADME) profile of compound 72 predicted it as an excellent orally bioavailable brain-penetrant molecule. To confirm these results, compound 72 was synthesized and spectrally characterized. The selectivity and inhibitory potential (IC50) of compound 72 was estimated by in vitro BACE-1 and BACE-2 FRET assay. Compound 72 was found to inhibit BACE-1 with IC50 121.65 nM, while it was found to be less potent on BACE-2 (IC50 480.92 nM), as compared to Atabecestat (BACE-1, IC50 13.25 nM and BACE-2, IC50 7.15 nM).

Molecular mechanism underlying effect of D93 and D289 protonation states on inhibitor-BACE1 binding: exploration from multiple independent Gaussian accelerated molecular dynamics and deep learning

BACE1 has been regarded as an essential drug design target for treating Alzheimer's disease (AD). Multiple independent Gaussian accelerated molecular dynamics simulations (GaMD), deep learning (DL), and molecular mechanics general Born surface area (MM-GBSA) method are integrated to elucidate the molecular mechanism underlying the effect of D93 and D289 protonation on binding of inhibitors OV6 and 4B2 to BACE1. The GaMD trajectory-based DL successfully identifies significant function domains. Dynamic analysis shows that the protonation of D93 and D289 strongly affects the structural flexibility and dynamic behaviour of BACE1. Free energy landscapes indicate that inhibitor-bound BACE1s have more conformational states in the protonated states than the wild-type (WT) BACE1, and show more binding poses of inhibitors. Binding affinities calculated using the MM-GBSA method indicate that the protonation of D93 and D289 highly disturbs the binding ability of inhibitors to BACE1. In addition, the protonation of two residues significantly affects the hydrogen bonding interactions (HBIs) of OV6 and 4B2 with BACE1, altering their binding activity to BACE1. The binding hot spots of BACE1 recognized by residue-based free energy estimations provide rational targeting sites for drug design towards BACE1. This study is anticipated to provide theoretical aids for drug development towards treatment of AD.

Structure and dynamics of whole-sequence homology model of ORF3a protein of SARS-CoV-2: An insight from microsecond molecular dynamics simulations

The ORF3a is a large accessory protein in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which plays an important role in virulence and viral replication; especially in inflammasome activation and apoptosis. However,, the existing cryo-EM structure of SARS-CoV-2 ORF3a is incomplete, . making it challenging to understand its structural and functional features. The aim of this study is to investigate the dynamic behaviors of the full-sequence homology model of ORF3a and compare it with the cryo-EM structure using microsecond molecular dynamics simulations. The previous studies indicated that the unresolved residues of the cryo-EM structure are not only involved in the pathogenesis of the SARS-CoV-2 but also exhibit a significant antigenicity. The dynamics scenario of homology model revealed higher RMSD, Rg, and SASA values with stable pattern when compared to the cryo-EM structure. Moreover, the RMSF analysis demonstrated higher fluctuations at specific positions (1-43, 97-110, 172-180, 219-243) in the model structure, whereas the cryo-EM structure displayed lower overall drift (except 1-43) in comparison to the model structure.Secondary structural features indicated that a significant unfolding in the transmembrane domains and β-strand at positions 166 to 172, affecting the stability and compactness of the cryo-EM structure , whereas the model exhibited noticeable unfolding in transmembrane domains and small-coiled regions in the N-terminal. , The results from molecular docking and steered molecular dynamics investigations showed the model structure had a greater number of non-bonding interactions, leading to enhanced stability when compared to the cryo-EM structure. Consequently, higher forces were necessary for unbinding of the baricitinib and ruxolitinib inhibitors from the model structure.. Our findings can help better understanding of the significance of unresolved residues at the molecular level. Additionally, this information can guide researchers for experimental endeavors aimed at completing the full-sequence structure of the ORF3a.Communicated by Ramaswamy H. Sarma.

Flavonoid and Chalcone Scaffolds as Inhibitors of BACE1: Recent Updates

Flavonoids and chalcones are two major classes of chemical moieties that have a vast background of pharmacological activities. Chalcone is a subclass of flavonoids whose therapeutic potential has been implicated due to an array of bioactivities. A lot of research works have shown interest in investigating the neuroprotective effect of these molecules, and have revealed them to be much more potent molecules that can be used to treat neurodegenerative disorders. Beta-site APP cleaving enzyme (BACE1), which is majorly found in the brain, is one of the reasons behind the development of Alzheimer's disease (AD). Flavonoids and chalcones have proven clinical data that they inhibit the production of Aβ plaques that are involved in the progression of AD. In this article, we have provided a detailed chronological review of the research work on the BACE1 inhibiting potency of both flavonoids and chalcones. Almost all the flavonoids and chalcones mentioned in this article have shown very good in vitro and in vivo BACE1 inhibiting activity. The docking studies and the structural importance of some BACE1-inhibiting flavonoids, as well as chalcones, are also mentioned here.

Binding selectivity analysis of AURKs inhibitors through molecular dynamics simulation studies

Aurora kinases (AURKs) have been identified as promising biological targets for the treatment of cancer. In this study, molecular dynamics simulations were employed to investigate the binding selectivity of three inhibitors (HPM, MPY, and VX6) towards AURKA and AURKB by predicting their binding free energies. The results show that the inhibitors HPM, MPY, and VX6 have more favorable interactions with AURKB as compared to AURKA. The binding energy decomposition analysis revealed that four common residue pairs (L139, L83), (V147, V91), (L210, L154), and (L263, L207) showed significant binding energies with HPM, MPY, and VX6, hence responsible for the binding selectivity of AURKA and AURKB to the inhibitors. The MD trajectory analysis also revealed that the inhibitors affect the dynamic flexibility of protein structure, which is also responsible for the partial selectivity of HPM, MPY, and VX6 towards AURKA and AURKB. As expected, this study provides useful insights for the design of potential inhibitors with high selectivity for AURKA and AURKB.

Molecular insights and optimization strategies for the competitive binding of engineered ACE2 proteins: a multiple replica molecular dynamics study

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) continues to spread globally, and rapid viral evolution and the emergence of new variants pose challenges to pandemic control. During infection, the spike protein of SARS-CoV-2 interacts with the human ACE2 protein via its receptor binding domain (RBD), and it is known that engineered forms of ACE2 can compete with wild-type (WT) ACE2 for binding to inhibit infection. Here, we conducted multiple replica molecular dynamics (MRMD) simulations to study the mechanisms of the engineered ACE2 variants 3N39 and 3N94 and provide directions for optimization. Our findings reveal that engineered ACE2 is notably more efficacious in systems that show weaker binding to WT ACE2 (i.e., WT and BA.1 RBD), but also faces immune escape as the virus evolves. Moreover, by modifying residue types near the binding interface, engineered ACE2 alters the electrostatic potential distribution and reconfigures the hydrogen bonding network, which results in modified binding to the RBD. However, this structural rearrangement does not occur in all RBD variants. In addition, we identified potentially engineerable beneficial residues and potentially engineerable detrimental residues in both ACE2 and RBD. Functional conservation can thus enable the optimization of these residues and improve the binding competitiveness of engineered ACE2, which therefore provides additional immune escape prevention. Finally, we conclude that these findings have implications for understanding the mechanisms responsible for engineered ACE2 and can help us to develop engineered ACE2 proteins that show superior performance.

Investigation of the Inhibition Mechanism of Xanthine Oxidoreductase by Oxipurinol: A Computational Study

Xanthine oxidoreductase (XOR) is an enzyme found in various organisms. It converts hypoxanthine to xanthine and urate, which are crucial steps in purine elimination in humans. Elevated uric acid levels can lead to conditions like gout and hyperuricemia. Therefore, there is significant interest in developing drugs that target XOR for treating these conditions and other diseases. Oxipurinol, an analogue of xanthine, is a well-known inhibitor of XOR. Crystallographic studies have revealed that oxipurinol directly binds to the molybdenum cofactor (MoCo) in XOR. However, the precise details of the inhibition mechanism are still unclear, which would be valuable for designing more effective drugs with similar inhibitory functions. In this study, molecular dynamics and quantum mechanics/molecular mechanics calculations are employed to investigate the inhibition mechanism of XOR by oxipurinol. The study examines the structural and dynamic effects of oxipurinol on the pre-catalytic structure of the metabolite-bound system. Our results provide insights on the reaction mechanism catalyzed by the MoCo center in the active site, which aligns well with experimental findings. Furthermore, the results provide insights into the residues surrounding the active site and propose an alternative mechanism for developing alternative covalent inhibitors.

Identification Mechanism of BACE1 on Inhibitors Probed by Using Multiple Separate Molecular Dynamics Simulations and Comparative Calculations of Binding Free Energies

β-amyloid cleaving enzyme 1 (BACE1) is regarded as an important target of drug design toward the treatment of Alzheimer's disease (AD). In this study, three separate molecular dynamics (MD) simulations and calculations of binding free energies were carried out to comparatively determine the identification mechanism of BACE1 for three inhibitors, 60W, 954 and 60X. The analyses of MD trajectories indicated that the presence of three inhibitors influences the structural stability, flexibility and internal dynamics of BACE1. Binding free energies calculated by using solvated interaction energy (SIE) and molecular mechanics generalized Born surface area (MM-GBSA) methods reveal that the hydrophobic interactions provide decisive forces for inhibitor-BACE1 binding. The calculations of residue-based free energy decomposition suggest that the sidechains of residues L91, D93, S96, V130, Q134, W137, F169 and I179 play key roles in inhibitor-BACE1 binding, which provides a direction for future drug design toward the treatment of AD.

Deciphering Selectivity Mechanism of BRD9 and TAF1(2) toward Inhibitors Based on Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations

BRD9 and TAF1(2) have been regarded as significant targets of drug design for clinically treating acute myeloid leukemia, malignancies, and inflammatory diseases. In this study, multiple short molecular dynamics simulations combined with the molecular mechanics generalized Born surface area method were employed to investigate the binding selectivity of three ligands, 67B, 67C, and 69G, to BRD9/TAF1(2) with IC50 values of 230/59 nM, 1400/46 nM, and 160/410 nM, respectively. The computed binding free energies from the MM-GBSA method displayed good correlations with that provided by the experimental data. The results indicate that the enthalpic contributions played a critical factor in the selectivity recognition of inhibitors toward BRD9 and TAF1(2), indicating that 67B and 67C could more favorably bind to TAF1(2) than BRD9, while 69G had better selectivity toward BRD9 over TAF1(2). In addition, the residue-based free energy decomposition approach was adopted to calculate the inhibitor–residue interaction spectrum, and the results determined the gatekeeper (Y106 in BRD9 and Y1589 in TAF1(2)) and lipophilic shelf (G43, F44, and F45 in BRD9 and W1526, P1527, and F1528 in TAF1(2)), which could be identified as hotspots for designing efficient selective inhibitors toward BRD9 and TAF1(2). This work is also expected to provide significant theoretical guidance and insightful molecular mechanisms for the rational designs of efficient selective inhibitors targeting BRD9 and TAF1(2).

Binding modes of GDP, GTP and GNP to NRAS deciphered by using Gaussian accelerated molecular dynamics simulations

Probing binding modes of GDP, GTP and GNP to NRAS are of significance for understanding the regulation mechanism on the activity of RAS proteins. Four separate Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the apo, GDP-, GTP- and GNP-bound NRAS. Dynamics analyses suggest that binding of three ligands highly affects conformational states of the switch domains from NRAS, which disturbs binding of NRAS to its effectors. The analyses of free energy landscapes (FELs) indicate that binding of GDP, GTP and GNP induces more energetic states of NRAS compared to the apo NRAS but the presence of GNP makes the switch domains more ordered than binding of GDP and GNP. The information of interaction networks of ligands with NRAS reveals that the π-π interaction of residue F28 and the salt bridge interactions of K16 and D119 with ligands stabilize binding of GDP, GTP and GNP to NRAS. Meanwhile magnesium ion plays a bridge role in interactions of ligands with NRAS, which is favourable for associations of GDP, GTP and GNP with NRAS. This work is expected to provide useful information for deeply understanding the function and activity of NRAS.

Decoding the Identification Mechanism of an SAM-III Riboswitch on Ligands through Multiple Independent Gaussian-Accelerated Molecular Dynamics Simulations

S-Adenosyl-l-methionine (SAM)-responsive riboswitches play a central role in the regulation of bacterial gene expression at the level of transcription attenuation or translation inhibition. In this study, multiple independent Gaussian-accelerated molecular dynamics simulations were performed to decipher the identification mechanisms of SAM-III (S_MK) on ligands SAM, SAH, and EEM. The results reveal that ligand binding highly affects the structural flexibility, internal dynamics, and conformational changes of SAM-III. The dynamic analysis shows that helices P3 and P4 as well as two junctions J23 and J24 of SAM-III are highly susceptible to ligand binding. Analyses of free energy landscapes suggest that ligand binding induces different free energy profiles of SAM-III, which leads to the difference in identification sites of SAM-III on ligands. The information on ligand-nucleotide interactions not only uncovers that the π-π, cation-π, and hydrogen bonding interactions drive identification of SAM-III on the three ligands but also reveals that different electrostatic properties of SAM, SAH, and EEM alter the active sites of SAM-III. Meanwhile, the results also verify that the adenine group of SAM, SAH, and EEM is well recognized by conserved nucleotides G7, A29, U37, A38, and G48. We expect that this study can provide useful information for understanding the applications of SAM-III in chemical, synthetic RNA biology, and biomedical fields.

Investigating the Dynamic Binding Behavior of PMX53 Cooperating with Allosteric Antagonist NDT9513727 to C5a Anaphylatoxin Chemotactic Receptor 1 through Gaussian Accelerated Molecular Dynamics and Free-Energy Perturbation Simulations

C5a anaphylatoxin chemotactic receptor 1 (C5aR1) is an important target in anti-inflammatory therapeutics. The cyclic peptide antagonist PMX53 binds to the orthosteric site located in the extracellular vestibule of C5aR1, and the non-peptide antagonist NDT9513727 binds to the allosteric site formed by the middle region of TM3 (trans-membrane helix), TM4, and TM5. We catch a sight of the variational binding mode of PMX53 during the Gaussian accelerated molecular dynamic (GaMD) simulations. In the binary complex of C5aR1 and PMX53, the PMX53 takes a dynamic binding mechanism during the simulation. Namely, the side chain of Arg⁶ of PMX53 extends to TM6-TM7 (pose 1) or swings to TM5 (pose 2), forming a salt bridge with Glu199. Meanwhile, in the ternary complex of C5aR1 with PMX53 and NDT9513727, the side chain of Arg⁶ of PMX53 swings to TM5 (pose 2) from extending to TM6-TM7 (pose 1) at the beginning of the GaMD simulation. In subsequent simulation, PMX53 stabilizes in the pose 2 binding mode by forming a stable salt bridge with Glu199. The free-energy perturbation (FEP) calculations demonstrate that pose 1 (ΔG_binding = -10.94 kcal/mol) is more stable in the binary complex and pose 2 (ΔG_binding = -7.91 kcal/mol) is unstable because of highly dynamic TM5. NDT9513727 interacts directly with TM4 and TM5 and stabilizes the hydrophobic stack between the extracellular sides of the two helices. Therefore, pose 2 (ΔG_binding = -16.27 kcal/mol) is notably stable than pose 1 (ΔG_binding = -9.78 kcal/mol) in the ternary complex. The identification of a novel binding mode of PMX53 and the detailed structural information of PMX53 interacting with a receptor obtained by GaMD simulations will be helpful in designing potent antagonists of C5aR1.

Multitargeted Molecular Docking and Dynamic Simulation Studies of Bioactive Compounds from Rosmarinus officinalis against Alzheimer's Disease

Alzheimer's disease (AD) has been associated with the hallmark features of cholinergic dysfunction, amyloid beta (Aβ) aggregation and impaired synaptic transmission, which makes the associated proteins, such as β-site amyloid precursor protein cleaving enzyme 1 (BACE I), acetylcholine esterase (AChE) and synapsin I, II and III, major targets for therapeutic intervention. The present study investigated the therapeutic potential of three major phytochemicals of Rosmarinus officinalis, ursolic acid (UA), rosmarinic acid (RA) and carnosic acid (CA), based on their binding affinity with AD-associated proteins. Detailed docking studies were conducted using AutoDock vina followed by molecular dynamic (MD) simulations using Amber 20. The docking analysis of the selected molecules showed the binding energies of their interaction with the target proteins, while MD simulations comprising root mean square deviation (RMSD), root mean square fluctuation (RMSF) and molecular mechanics/generalized born surface area (MM/GBSA) binding free energy calculations were carried out to check the stability of bound complexes. The drug likeness and the pharmacokinetic properties of the selected molecules were also checked through the Lipinski filter and ADMETSAR analysis. All these bioactive compounds demonstrated strong binding affinity with AChE, BACE1 and synapsin I, II and III. The results showed UA and RA to be potential inhibitors of AChE and BACE1, exhibiting binding energies comparable to those of donepezil, used as a positive control. The drug likeness and pharmacokinetic properties of these compounds also demonstrated drug-like characteristics, indicating the need for further in vitro and in vivo investigations to ascertain their therapeutic potential for AD.

Decoding drug resistant mechanism of V32I, I50V and I84V mutations of HIV-1 protease on amprenavir binding by using molecular dynamics simulations and MM-GBSA calculations

Mutations V32I, I50V and I84V in the HIV-1 protease (PR) induce drug resistance towards drug amprenavir (APV). Multiple short molecular dynamics (MSMD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) method were utilized to investigate drug-resistant mechanism of V32I, I50V and I84V towards APV. Dynamic information arising from MSMD simulations suggest that V32I, I50V and I84V highly affect structural flexibility, motion modes and conformational behaviours of two flaps in the PR. Binding free energies calculated by MM-GBSA method suggest that the decrease in binding enthalpy and the increase in binding entropy induced by mutations V32I, I50V and I84V are responsible for drug resistance of the mutated PRs on APV. The energetic contributions of separate residues on binding of APV to the PR show that V32I, I50V and I84V highly disturb the interactions of two flaps with APV and mostly drive the decrease in binding ability of APV to the PR. Thus, the conformational changes of two flaps in the PR caused by V32I, I50V and I84V play key roles in drug resistance of three mutated PR towards APV. This study can provide useful dynamics information for the design of potent inhibitors relieving drug resistance.

Deciphering Conformational Changes of the GDP-Bound NRAS Induced by Mutations G13D, Q61R, and C118S through Gaussian Accelerated Molecular Dynamic Simulations

The conformational changes in switch domains significantly affect the activity of NRAS. Gaussian-accelerated molecular dynamics (GaMD) simulations of three separate replicas were performed to decipher the effects of G13D, Q16R, and C118S on the conformational transformation of the GDP-bound NRAS. The analyses of root-mean-square fluctuations and dynamics cross-correlation maps indicated that the structural flexibility and motion modes of the switch domains involved in the binding of NRAS to effectors are highly altered by the G13D, Q61R, and C118Smutations. The free energy landscapes (FELs) suggested that mutations induce more energetic states in NRAS than the GDP-bound WT NRAS and lead to high disorder in the switch domains. The FELs also indicated that the different numbers of sodium ions entering the GDP binding regions compensate for the changes in electrostatic environments caused by mutations, especially for G13D. The GDP-residue interactions revealed that the disorder in the switch domains was attributable to the unstable hydrogen bonds between GDP and two residues, V29 and D30. This work is expected to provide information on the energetic basis and dynamics of conformational changes in switch domains that can aid in deeply understanding the target roles of NRAS in anticancer treatment.

Exploring the deactivation mechanism of human β2 adrenergic receptor by accelerated molecular dynamic simulations

The β2 adrenergic receptor (β2AR), one of important members of the G protein coupled receptors (GPCRs), has been suggested as an important target for cardiac and asthma drugs. Two replicas of Gaussian accelerated molecular dynamics (GaMD) simulations are performed to explore the deactivation mechanism of the active β2AR bound by three different substrates, including the agonist (P0G), antagonist (JTZ) and inverse agonist (JRZ). The simulation results indicate that the Gs protein is needed to stabilize the active state of the β2AR. Without the Gs protein, the receptor could transit from the active state toward the inactive state. During the transition process, helix TM6 moves toward TM3 and TM5 in geometric space and TM5 shrinks upwards. The intermediate state is captured during the transition process of the active β2AR toward the inactive one, moreover the changes in hydrophobic interaction networks between helixes TM3, TM5, and TM6 and the formation of a salt bridge between residues Arg3.50 and Glu6.30 drive the transition process. We expect that this finding can provide energetic basis and molecular mechanism for further understanding the function and target roles of the β2AR.

Free Energy Profiles Relating With Conformational Transition of the Switch Domains Induced by G12 Mutations in GTP-Bound KRAS

Mutations of G12 in KRAS have been involved in different cancers. Multiple replica-Gaussian accelerated molecular dynamics (MR-GaMD) simulations are applied to investigate conformational changes of the switch domains caused by G12C, G12D and G12R. Free energy landscapes suggest that G12C, G12D and G12R induce more energetic states compared to the GTP-bound WT KRAS and make the conformations of the switch domains more disordered, which disturbs bindings of KRAS to effectors. Dynamics analyses based on MR-GaMD trajectory show that G12C, G12D and G12R not only change structural flexibility of the switch domains but also affect their motion behavior, indicating that these three mutations can be used to tune the activity of KRAS. The analyses of interaction networks verify that the instability in interactions of the GTP with the switch SⅠ plays an important role in the high disorder states of the switch domain. This work is expected to provide useful information for deeply understanding the function of KRAS.

Binding mechanism of inhibitors to SARS-CoV-2 main protease deciphered by multiple replica molecular dynamics simulations

The outbreak caused by SARS-CoV-2 has received extensive worldwide attention. As the main protease (M^pro) in SARS-CoV-2 has no human homologues, it is feasible to reduce the possibility of targeting the host protein by accidental drugs. Thus, M^pro has been an attractive target of efficient drug design for anti-SARS-CoV-2 treatment. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), and the molecular mechanics-generalized Born surface area (MM-GBSA) method were integrated together to decipher the binding mechanism of four inhibitors masitinib, O6K, FJC and GQU to M^pro. The results indicate that the binding of four inhibitors clearly affects the structural flexibility and internal dynamics of M^pro along with dihedral angle changes of key residues. The analysis of FELs unveils that the stability in the relative orientation and geometric position of inhibitors to M^pro is favorable for inhibitor binding. Residue-based free energy decomposition reveals that the inhibitor-M^pro interaction networks involving hydrogen bonding interactions and hydrophobic interactions provide significant information for the design of potent inhibitors against M^pro. The hot spot residues including H41, M49, F140, N142, G143, C145, H163, H164, M165, E166 and Q189 identified by computational alanine scanning are considered as reliable targets of clinically available inhibitors inhibiting the activities of M^pro.

Q61 mutant-mediated dynamics changes of the GTP-KRAS complex probed by Gaussian accelerated molecular dynamics and free energy landscapes

Understanding the molecular mechanism of the GTP-KRAS binding is significant for improving the target roles of KRAS in cancer treatment. In this work, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations were applied to decode the effect of Q61A, Q61H and Q61L on the activity of KRAS. Dynamics analyses based on MR-GaMD trajectory reveal that motion modes and dynamics behavior of the switch domain in KRAS are heavily affected by the three Q61 mutants. Information of free energy landscapes (FELs) shows that Q61A, Q61H and Q61L induce structural disorder of the switch domain and disturb the activity of KRAS. Analysis of the interaction network uncovers that the decrease in the stability of hydrogen bonding interactions (HBIs) of GTP with residues V29 and D30 induced by Q61A, Q61H and Q61L is responsible for the structural disorder of the switch-I and that in the occupancy of the hydrogen bond between GTP and residue G60 leads to the structural disorder of the switch-II. Thus, the high disorder of the switch domain caused by three current Q61 mutants produces a significant effect on binding of KRAS to its effectors. This work is expected to provide useful information for further understanding function and target roles of KRAS in anti-cancer drug development.

Theoretical exploration of the binding selectivity of inhibitors to BRD7 and BRD9 with multiple short molecular dynamics simulations

Bromodomain-containing proteins 7 and 9 (BRD7 and BRD9) have been considered as potential targets of clinical drug design toward treatment of human cancers and other diseases. Multiple short molecular dynamics simulations and binding free energy predictions were carried out to decipher the binding selectivity of three inhibitors 4L2, 5U6, and 6KT toward BRD7 and BRD9. The results show that 4L2 has more favorable binding ability to BRD7 over BRD9 compared to 5U6 and 6KT, while 5U6 and 6KT possess more favorable associations with BRD9 than BRD7. Furthermore, estimations of residue-based free energy decompositions further identify that four common residue pairs, including (F155, F44), (V160, V49), (Y168, Y57) and (Y217, Y106) in (BRD7, BRD9) generate obvious binding differences with 4L2, 5U6, and 6KT, which mostly drives the binding selectivity of 4L2, 5U6, and 6KT to BRD7 and BRD9. Dynamic information arising from trajectory analysis also suggests that inhibitor bindings affect structural flexibility and motion modes, which is responsible for the partial selectivity of 4L2, 5U6, and 6KT toward BRD7 and BRD9. As per our expectation, this study theoretically provides useful hints for design of dual inhibitors with high selectivity on BRD7 and BRD9.

Bromodomains (BRDs) are structurally conserved epigenetic reader modules observed in numerous chromatin- and transcription-associated proteins that have a capability to identify acetylated lysine residues.

Deciphering binding mechanism of inhibitors to SARS-COV-2 main protease through multiple replica accelerated molecular dynamics simulations and free energy landscapes

The pneumonia outbreak caused by the SARS-CoV-2 virus poses a serious threat to human health and the world economy. Development of safe and highly effective antiviral drugs is of great...

Binding Mechanism between Acetylcholinesterase and Drugs Pazopanib and Lapatinib: Biochemical and Biophysical Studies

Tyrosine kinase inhibitors (TKIs) are antitumor compounds that prevent the phosphorylation of proteins in a biological environment. However, the multitarget performance of TKIs promotes them as possible candidates for drug repositioning. In this work, interaction and inhibition studies through spectroscopic and computational techniques to evaluate the binding effectiveness of lapatinib and pazopanib TKIs to acetylcholinesterase (AChE) are reported. The results indicated potent inhibition at the μM level. The types of inhibition were identified, with pazopanib acting through non-competitive inhibition and lapatinib through acompetitive inhibition. The fluorescence suppression studies indicate a static mechanism for lapatinib-AChE and pazopanib-AChE systems, with a binding constant in the order of 10⁵ M^-1. The obtained thermodynamic parameters reveal interactions driven by van der Waals forces and hydrogen bonds in the lapatinib-AChE system (ΔH° and ΔS° < 0). In contrast, the pazopanib-AChE system shows positive ΔH° and ΔS°, characteristic of hydrophobic interactions. The Foster resonance energy transfer study supports the fluorescence studies performed. The 3D fluorescence studies suggest changes in the microenvironment of the tryptophan and tyrosine residues of the protein in contact with lapatinib and pazopanib. The results suggest effective inhibition and moderate interaction of the drugs with AChE, making them interesting for conducting more in-depth repositioning studies as AChE inhibitors.

Decoding molecular mechanism underlying binding of drugs to HIV-1 protease with molecular dynamics simulations and MM-GBSA calculations

HIV-1 protease (PR) is thought to be efficient targets of anti-AIDS drug design. Molecular dynamics (MD) simulations and multiple post-processing analysis technologies were applied to decipher molecular mechanism underlying binding of three drugs Lopinavir (LPV), Nelfinavir (NFV) and Atazanavir (ATV) to the PR. Binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) suggest that compensation between binding enthalpy and entropy plays a vital role in binding of drugs to PR. Dynamics analyses show that binding of LPV, NFV and ATV highly affects structural flexibility, motion modes and dynamics behaviour of the PR, especially for two flaps. Computational alanine scanning and interaction network analysis verify that although three drugs have structural difference, they share similar binding modes to the PR and common interaction clusters with the PR. The current findings also confirm that residues located interaction clusters, such as Asp25/Asp25', Gly27/Gly27', Ala28/Ala28', Asp29, Ile47/Ile47', Gly49/Gly49', Ile50/Ile50', Val82/Val82' and Ile84/Ile84, can be used as efficient targets of clinically available inhibitors towards the PR.

Cyclopentanone Derivative Attenuates Memory Impairment by Inhibiting Amyloid Plaques Formation in the 5xFAD Mice

Alzheimer's disease (AD) is a chronic neurodegenerative disorder. This study was designed to investigate the effects of cyclopentanone derivative i.e., 2-(hydroxyl-(3-nitrophenyl)methyl)cyclopentanone (3NCP) on behavior, amyloid β (Aβ) plaque deposition, and βAPP cleaving enzyme-1 (BACE-1) expression in the 5xFAD mouse brain. In this study, computational studies were conducted to predict the binding mode of the 3NCP with target sites of the β-secretase. In vivo studies were performed on the 5xFAD mice model of AD using different behavioral test models like light/dark box, elevated plus maze (EPM), and the Barnes maze tests for the assessment of anxiety, spatial learning and memory. The thioflavin-S staining, immunohistochemistry (IHC), and RT-PCR studies were carried out to find the effect of the 3NCP on the β-amyloid plaques formation and BACE-1 expression. The results of the computational studies showed that the 3NCP has excellent binding affinities for beta-secretase. The light/dark box study depicted that the 3NCP does not cause anxiety. The 3NCP treatment effects in the EPM and Barnes maze tests showed a significant effect on learning and memory. Furthermore, the results of the thioflavin staining and IHC revealed that the 3NCP significantly reduced the formation of the beta-amyloid plaques in brain tissues. Moreover, the RT-PCR study showed that 3NCP significantly reduced the BACE-1 expression in the brain. Conclusively, the results of the current study demonstrate that the 3NCP may be a potential candidate for AD treatment in the future.

Conformational transformation of switch domains in GDP/K-Ras induced by G13 mutants: An investigation through Gaussian accelerated molecular dynamics simulations and principal component analysis

Mutations in K-Ras are involved in a large number of all human cancers, thus, K-Ras is regarded as a promising target for anticancer drug design. Understanding the target roles of K-Ras is important for providing insights on the molecular mechanism underlying the conformational transformation of the switch domains in K-Ras due to mutations. In this study, multiple replica Gaussian accelerated molecular (MR-GaMD) simulations and principal component analysis (PCA) were applied to probe the effect of G13A, G13D and G13I mutations on conformational transformations of the switch domains in GDP-associated K-Ras. The results suggest that G13A, G13D and G13I enhance the structural flexibility of the switch domains, change the correlated motion modes of the switch domains and strengthen the total motion strength of K-Ras compared with the wild-type (WT) K-Ras. Free energy landscape analyses not only show that the switch domains of the GDP-bound inactive K-Ras mainly exist as a closed state but also indicate that mutations evidently alter the free energy profile of K-Ras and affect the conformational transformation of the switch domains between the closed and open states. Analyses of hydrophobic interaction contacts and hydrogen bonding interactions show that the mutations scarcely change the interaction network of GDP with K-Ras and only disturb the interaction of GDP with the switch (SW1). In summary, two newly introduced mutations, G13A and G13I, play similar adjustment roles in the conformational transformations of two switch domains to G13D and are possibly utilized to tune the activity of K-Ras and the binding of guanine nucleotide exchange factors.