Molecular Mechanism of Binding Selectivity of Inhibitors toward BACE1 and BACE2 Revealed by Multiple Short Molecular Dynamics Simulations and Free-Energy Predictions

Molecular recognition between volatile molecules and odorant binding proteins 7 by homology modeling, molecular docking and molecular dynamics simulation

BACKGROUND

Odorant-binding proteins (OBPs) in insects are key to detection and recognition of external chemical signals associated with survival. OBP7 in Spodoptera frugiperda's larval stage (SfruOBP7) may search for host plants by sensing plant volatiles, which are important sources of pest attractants and repellents. However, the atomic-level basis of binding modes remains elusive.

RESULTS

SfruOBP7 structure was constructed through homology modeling, and complex models of six plant volatiles ((E)-2-hexenol, α-pinene, (Z)-3-hexenyl acetate, lauric acid, O-cymene and 1-octanol) and SfruOBP7 were obtained through molecular docking. To study the detailed interactions between the six plant volatile molecules and SfruOBP7, we conducted three 300 ns molecular dynamics simulations for each study object. The correlation coefficients between binding free energy obtained by molecular mechanics/generalized Born surface area together with solvated interaction energy methods and experimental values are 0.90 and 0.88, respectively, showing a good correlation. By comparing binding free energy along with interaction patterns between SfruOBP7 and the six volatile molecules, hotspot residues of SfruOBP7 when binding with different volatile molecules were determined. Hydrophobic interactions stemming from van der Waals interactions play a significant role in SfruOBP7 and these plant volatile systems.

CONCLUSION

The optimized three-dimensional structure of SfruOBP7 and its binding modes with six plant volatiles revealed their interactions, thus providing a means for estimating the binding energies of other plant volatiles. Our study will help to guide the rational design of effective and selective insect attractants. © 2024 Society of Chemical Industry.

Insight into the molecular recognition of human and polar bear pregnane X receptor by three organic pollutants using molecular docking and molecular dynamics simulations

Pregnane X receptor (PXR) is a heterologous biosensor that is involved in the metabolic pathway of environmental pollutants, regulating the transcription of genes involved in biotransformation. There are significant differences in the selectivity and specificity of organic pollutants (OPs) toward polar bear PXR (pbPXR) and human PXR (hPXR), but the detailed dynamical characteristics of their interactions are unclear. Homology Modeling, molecular docking, molecular dynamics simulation, and free energy calculation were used to analyze the recognition of pbPXR and hPXR by three OPs: BPA, chlordane and toxaphene. Comparing interaction patterns along with binding free energy of pbPXR and hPXR with these three OPs revealed that although pbPXR and hPXR interact similar with these three OPs, these OPs have different effects on the internal dynamics of pbPXR and hPXR. This results in significant alterations in the interaction of key residues near Leu209, Met243, Phe288, Met323, and His407 with OPs, thereby influencing their binding energy. Non-polar interactions, especially van der Waals interactions, were found to be the dominating factors in interacting of these OPs with PXRs. The region surrounding these key residues facilitates hydrophobic contacts with PXR, which are crucial for the selective activation of PXRs in different species by these three OPs. These findings are of significant guidance in understanding the impacts of environmental endocrine disruptors on different organisms.

Practical Three-Component Regioselective Synthesis of Drug-Like 3-Aryl(or heteroaryl)-5,6-dihydrobenzo[h]cinnolines as Potential Non-Covalent Multi-Targeting Inhibitors To Combat Neurodegenerative Diseases

Neurodegenerative diseases (NDs) are one of the prominent health challenges facing contemporary society, and many efforts have been made to overcome and (or) control it. In this research paper, we described a practical one-pot two-step three-component reaction between 3,4-dihydronaphthalen-1(2H)-one (1), aryl(or heteroaryl)glyoxal monohydrates (2a-h), and hydrazine monohydrate (NH2NH2•H2O) for the regioselective preparation of some 3-aryl(or heteroaryl)-5,6-dihydrobenzo[h]cinnoline derivatives (3a-h). After synthesis and characterization of the mentioned cinnolines (3a-h), the in silico multi-targeting inhibitory properties of these heterocyclic scaffolds have been investigated upon various Homo sapiens-type enzymes, including hMAO-A, hMAO-B, hAChE, hBChE, hBACE-1, hBACE-2, hNQO-1, hNQO-2, hnNOS, hiNOS, hPARP-1, hPARP-2, hLRRK-2(G2019S), hGSK-3β, hp38α MAPK, hJNK-3, hOGA, hNMDA receptor, hnSMase-2, hIDO-1, hCOMT, hLIMK-1, hLIMK-2, hRIPK-1, hUCH-L1, hPARK-7, and hDHODH, which have confirmed their functions and roles in the neurodegenerative diseases (NDs), based on molecular docking studies, and the obtained results were compared with a wide range of approved drugs and well-known (with IC50, EC50, etc.) compounds. In addition, in silico ADMET prediction analysis was performed to examine the prospective drug properties of the synthesized heterocyclic compounds (3a-h). The obtained results from the molecular docking studies and ADMET-related data demonstrated that these series of 3-aryl(or heteroaryl)-5,6-dihydrobenzo[h]cinnolines (3a-h), especially hit ones, can really be turned into the potent core of new drugs for the treatment of neurodegenerative diseases (NDs), and/or due to the having some reactionable locations, they are able to have further organic reactions (such as cross-coupling reactions), and expansion of these compounds (for example, with using other types of aryl(or heteroaryl)glyoxal monohydrates) makes a new avenue for designing novel and efficient drugs for this purpose.

Molecular insights into inhibiting effects of lignin on cellulase investigated by molecular dynamics simulation

The hydrolysis of lignocellulose into fermentable monosaccharides using cellulases represents a critical stage in lignocellulosic bioconversion. However, the inactivation of cellulase in the presence of lignin is attributed to the high cost of biofinery. To address this challenge, a comprehensive investigation into the structure-function relationship underlying lignin-driven cellulase inactivation is essential. In this study, molecular docking and molecular dynamics (MD) simulations were employed to explore the impacts of lignin fragments on the catalytic efficiency of cellulase at the atomic level. The findings revealed that soluble lignin fragments and cellulose could spontaneously form stable complexes with cellulase, indicating a competitive binding scenario. The enzyme's structure remained unchanged upon binding to lignin. Furthermore, specific amino acid residues have been identified as involved in interactions with lignin and cellulose. Hydrophobic interactions were found to dominate the binding of lignin to cellulase. Based on the mechanisms underlying the interactions between lignin fragments and cellulase, decreased hydrophobicity and change in the charge of lignin may mitigate the inhibition of cellulase. Furthermore, site mutations and chemical modification are also feasible to improve the efficiency of cellulase. This study may contribute valuable insights into the design of more lignin-resistant enzymes and the optimization of lignocellulosic pretreatment technologies.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.

Kallistatin as an inhibitory protein against colorectal cancer cells through binding to LRP6

Kallistatin (KL) is a member of the serine proteinase inhibitor (serpin) family regulating oxidative stress, vascular relaxation, inflammation, angiogenesis, cell proliferation, and invasion. The heparin-binding site of Kallistatin has an important role in the interaction with LRP6 leading to the blockade of the Wnt signaling pathway. In this study, we aimed to explore the structural basis of the Kallistatin-LRP6E1E4 complex using in silico approaches and evaluating the anti-proliferative, apoptotic, and cell cycle arrest activities of Kallistatin in colon cancer lines. The molecular docking showed Kallistatin could bind to the LRP6E3E4 much stronger than LRP6E1E2. The Kallistatin-LRP6E1E2 and Kallistatin-LRP6E3E4 complexes were stable during Molecular Dynamics (MD) simulation. The Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) showed that the Kallistatin-LRP6E3E4 has a higher binding affinity compared to Kallistatin-LRP6E1E2. Kallistatin induced higher cytotoxicity and apoptosis in HCT116 compared to the SW480 cell line. This protein-induced cell-cycle arrest in both cell lines at the G1 phase. The B-catenin, cyclin D1, and c-Myc expression levels were decreased in response to treatment with Kallistatin in both cell lines while the LRP6 expression level was decreased in the HCT116 cell line. Kallistatin has a greater effect on the HCT116 cell line compared to the SW480 cell line. Kallistatin can be used as a cytotoxic and apoptotic-inducing agent in colorectal cancer cell lines.

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.

Roles of Accelerated Molecular Dynamics Simulations in Predictions of Binding Kinetic Parameters

Rational predictions on binding kinetics parameters of drugs to targets play significant roles in future drug designs. Full conformational samplings of targets are requisite for accurate predictions of binding kinetic parameters. In this review, we mainly focus on the applications of enhanced sampling technologies in calculations of binding kinetics parameters and residence time of drugs. The methods involved in molecular dynamics simulations are applied to not only probe conformational changes of targets but also reveal calculations of residence time that is significant for drug efficiency. For this review, special attention are paid to accelerated molecular dynamics (aMD) and Gaussian aMD (GaMD) simulations that have been adopted to predict the association or disassociation rate constant. We also expect that this review can provide useful information for future drug design.

Insights into the selectivity of a brain-penetrant CDK4/6 vs CDK1/2 inhibitor for glioblastoma used in multiple replica molecular dynamics simulations

Cyclin dependent kinases (CDKs) play an important role in cell cycle regulation and their dysfunction is associated with many cancers. That is why CDKs have been attractive targets for the treatment of cancer. Glioblastoma is a cancer caused by the aberrant expression of CDK4/6, so exploring the mechanism of the selection of CDK4/6 toward inhibitors relative to the other family members CDK1/2 is essential. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), molecular mechanics Poisson-Boltzmann/Generalized Born surface area (MM-PB/GBSA) and other methods were integrated to decipher the selectively binding mechanism of the inhibitor N1J to CDK4/6 and CDK1/2. Molecular electrostatic potential (MESP) analysis provides an explanation for the N1J selectivity. Residue-based free energy decomposition reveals that most of the hot residues are located at the same location of CDKs proteins, but the different types of residues in different proteins cause changes in binding energy, which is considered as a potential developmental direction to improve the selectivity of inhibitors to CDK4/6. These results provide insights into the source of inhibitor and CDK4/6 selectivity for the future development of more selective inhibitors.Communicated by Ramaswamy H. Sarma.

Virtual screening of novel mTOR inhibitors for the potential treatment of human colorectal cancer

The abnormal activation of the mTOR pathway is closely related to the occurrence and progression of cancer, especially colorectal cancer. In this study, a rational virtual screening strategy has been established and MT-5, a novel mTOR inhibitor with a quinoline scaffold, was obtained from the ChemDiv database. MT-5 showed potent kinase inhibitory activity (IC50: 8.90 μM) and antiproliferative effects against various cancer cell lines, especially HCT-116 cells (IC50: 4.61 μM), and this was 2.2-fold more potent than that of the cisplatin control (IC50: 9.99 μM). Western blot, cell migration, cycle arrest, and apoptosis assays were performed with HCT-116 cells to investigate the potential anticancer mechanism of MT-5. Metabolic stability results in vitro indicated that MT-5 exhibited good stability profiles in artificial gastrointestinal fluids, rat plasma, and liver microsomes. In addition, the key contribution of the residues around the binding pocket of MT-5 in binding to the mTOR protein was also investigated from a computational perspective.

Structure-Based Drug Design of Novel Triaminotriazine Derivatives as Orally Bioavailable IDH2R140Q Inhibitors with High Selectivity and Reduced hERG Inhibitory Activity for the Treatment of Acute Myeloid Leukemia

Neomorphic IDH2R140Q mutation is commonly found in acute myeloid leukemia (AML), and inhibiting its activity has been validated as an effective treatment for AML. Herein, we report a series of highly potent and selective IDH2R140Q inhibitors. Among them, compound 36 was identified as the most promising inhibitor, with an IC50 value of 29 nM and more than 490-fold selectivity over wild-type IDH2. The compound significantly suppressed D2HG production (IC50 = 10 nM) and induced differentiation in TF-1/IDH2R140Q cells. Furthermore, it showed reasonable pharmacokinetic properties with high bioavailability (F = 90.3%) and an appropriate half-life (T1/2 = 6.4 h). In vivo, oral administration of compound 36 at a dose of 25 mg/kg effectively reduced D2HG levels in the tumor of TF-1/IDH2R140Q xenograft mouse model. Besides, compound 36 displayed little effect on the hERG current. These results suggest that compound 36 has the potential to be an efficacious treatment for AML.

Binding Mechanism of Inhibitors to Heat Shock Protein 90 Investigated by Multiple Independent Molecular Dynamics Simulations and Prediction of Binding Free Energy

The heat shock protein (HSP90) has been an import target of drug design in the treatment of human disease. An exploration of the conformational changes in HSP90 can provide useful information for the development of efficient inhibitors targeting HSP90. In this work, multiple independent all-atom molecular dynamics (AAMD) simulations followed by calculations of the molecular mechanics generalized Born surface area (MM-GBSA) were performed to explore the binding mechanism of three inhibitors (W8Y, W8V, and W8S) to HSP90. The dynamics analyses verified that the presence of inhibitors impacts the structural flexibility, correlated movements, and dynamics behavior of HSP90. The results of the MM-GBSA calculations suggest that the selection of GB models and empirical parameters has important influences on the predicted results and verify that van der Waals interactions are the main forces that determine inhibitor-HSP90 binding. The contributions of separate residues to the inhibitor-HSP90 binding process indicate that hydrogen-bonding interactions (HBIs) and hydrophobic interactions play important roles in HSP90-inhibitor identifications. Moreover, residues L34, N37, D40, A41, D79, I82, G83, M84, F124, and T171 are recognized as hot spots of inhibitor-HSP90 binding and provide significant target sites of for the design of drugs related to HSP90. This study aims to contribute to the development of efficient inhibitors that target HSP90 by providing an energy-based and theoretical foundation.

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).

Structural and energetic basis of interaction between human estrogen-related receptor γ and environmental endocrine disruptors from multiple molecular dynamics simulations and free energy predictions

Environmental endocrine disruptors (EEDs), a class of molecules that are widespread in our environment, may adversely affect the endocrine system. Exploring the interactions between these compounds and their potential targets is important for assessing their role in the organism. Focused on the human estrogen-related receptor γ (hERRγ) with BPA, BPB, HPTE, BPE, BP(2,2)(Et), and BP(2,2)(MeO) complexes, respectively, we groped for the mechanisms of conformational changes and interactions of hERRγ when binding to these six EEDs by combining multiple molecular dynamics (MD) simulations with energy prediction (MM-PBSA and solvated interaction energy (SIE)). Dynamics analysis results revealed these six EEDs have different effects on the internal dynamics of hERRγ, resulting in significant changes in the interaction of key residues around Leu268, Val313, Leu345, and Phe435 with EEDs, and thus affected its binding energy with these EEDs. The energy calculations further demonstrated that van der Waals interactions are critical for these EEDs binding to hERRγ. These results present detailed molecular insight into the interaction features between EEDs and hERRγ and help guide the search for safer alternatives to BPA.

Molecular basis of ssDNA recognition by RBM45 protein of neurodegenerative disease from multiple molecular dynamics simulations and energy predictions

Neurodegenerative diseases (NDD) are a group of cognitive and behavioral disorders characterized by progressive loss of neuronal structure and function. As the population ages, the incidence is getting higher and higher, but there is currently no effective treatment. The details of RNA/DNA recognition by the RNA-binding protein RBM45 closely related to neurodegenerative diseases through its two tandem RNA-recognition domains at its N-terminus have important implications for structure-based drug discovery against degenerative diseases. To explore the key characteristics of interaction between ssDNA and RBM45, we performed multiple molecular dynamics (MD) simulations along with MM-PBSA energy prediction on the complexes of wild type (WT) and three mutant RBM45s (K100A, F124A/Y165A, and F29A/F70A/F124A/Y165A) with ssDNA, respectively. The findings suggest that these mutated residues of RBM45 modify the interaction of their surrounding residues with ssDNA, thereby affecting RBM45 protein binding to ssDNA. In contrast with WT RBM45 protein, variations in van der Waals and electrostatic interactions with ssDNA caused by these three RBM45 mutants are critical to affect binding between them. In addition, energy analysis showed that RBM45 is a specific ssDNA-binding protein. The results of our work provide valuable theoretical guidance for the design effective drugs of NDD.

Insights into pralsetinib resistance to the non-gatekeeper RET kinase G810C mutation through molecular dynamics simulations

OBJECTIVE

RET (rearranged during transfection) kinase, as a transmembrane receptor tyrosine kinase, is a therapeutic target for several human cancer such as non-small cell lung cancer (NSCLC) and thyroid cancer. Pralsetinib is a recently approved drug for the treatment of RET-driven NSCLC and thyroid cancers. A single point mutation G810C at the C-lobe of the RET kinase causes pralsetinib resistance to this non-gatekeeper variant. However, the detailed mechanism remains poorly understood.

METHODS

Here, multiple microsecond molecular dynamics (MD) simulations, molecular mechanics/generalized born surface area (MM/GBSA) binding free energy calculations, and community network analysis were performed to reveal the mechanism of pralsetinib resistance to the RET G810C mutant.

RESULTS

The simulations showed that the G810C mutation had a minor effect on the overall conformational dynamics of the RET kinase domain. Energetic analysis suggested that the G810C mutation reduced the binding affinity of pralsetinib to the mutant. Per-residue energy contribution and structural analyses revealed that the hydrogen bonding interactions between pralsetinib and the hinge residues Glu805 and Ala807 were disrupted in the G810C mutant, which were responsible for the decreased binding affinity of pralsetinib to the mutant.

CONCLUSIONS

The obtained results may provide understanding of the mechanism of pralsetinib resistance to the non-gatekeeper RET G810C mutant.

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.

Multiple Molecular Dynamics Simulations and Energy Analysis Unravel the Dynamic Properties and Binding Mechanism of Mutants HIV-1 Protease with DRV and CA-p2

PRS17, a variant of human immunodeficiency virus type I protease (HIV-1 PR), has 17 mutated residues showing high levels of multidrug resistance. To describe the effects of these mutated residues on the dynamic properties and the binding mechanism of PR with substrate and inhibitor, focused on six systems (two complexes of WT PR and PRS17 with inhibitor Darunavir (DRV), two complexes of WT PR and PRS17 with substrate analogue CA-p2, two unligand WT PR and PRS17), we performed multiple molecular dynamics (MD) simulations combined with MM-PBSA and solvated interaction energy (SIE) methods. For both the unligand PRs and ligand-PR complexes, the results from simulations revealed 17 mutated residues alter the flap-flap distance, the distance from flap regions to catalytic sites, and the curling degree of the flap tips. These mutated residues changed the flexibility of the flap region in PR, and thus affected its binding energy with DRV and CA-p2, resulting in differences in sensitivity. Hydrophobic cavity makes an important contribution to the binding of PR and ligands. And most noticeable of all, the binding of the guanidine group in CA-p2 and Arg8' of PRS17 is useful for increasing their binding ability. These results have important guidance for the further design of drugs against multidrug resistant PR. IMPORTANCE Developing effective anti-HIV inhibitors is the current requirement to cope with the emergence of the resistance of mutants. Compared with the experiments, MD simulations along with energy calculations help reduce the time and cost of designing new inhibitors. Based on our simulation results, we propose two factors that may help design effective inhibitors against HIV-1 PR: (i) importance of hydrophobic cavity, and (ii) introduction of polar groups similar to the guanidine group.

Computational study reveals substituted benzimidazole derivatives' binding selectivity to PI3Kδ and PI3Kγ

Phosphatidylinositol 3-kinase (PI3K) is a key regulatory kinase in the PI3K/AKT/mTOR signaling pathway, which is involved in the regulation of cell proliferation, differentiation, apoptosis, and angiogenesis. Class IA PI3K isoforms γ and δ share a highly homologous ATP binding site and are distinguished by only a few residues around the binding site. Subtype-selective inhibitors have been proven to have great advantages in tumor treatment. Preliminary studies have obtained PI3K inhibitors bearing a benzimidazole structural motif with a certain selectivity for PI3Kδ and PI3Kγ subtypes. On this basis, we investigated the selective inhibitory mechanism of PI3Kδ and PI3Kγ using four developed inhibitors via molecular docking, molecular dynamics, binding free energy calculations, and residue energy decomposition. This study could provide references for the further development of PI3K-isoform-selective inhibitors.

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.

Multiple Molecular Dynamics Simulations and Free-Energy Predictions Uncover the Susceptibility of Variants of HIV-1 Protease against Inhibitors Darunavir and KNI-1657

HIV-1 protease (PR) is considered to be the main targets of anti-AIDS drug design because of its role in the proteolytic processing of viral polyproteins. However, the emergence of drug-resistant HIV has become a major problem in the therapy of HIV-1-infected patients. Focused on the complexes of wild type (WT) PR and two mutant PRs (V32I/L33F/I54M/V82I and V32I/L33F/I54M/I84 V) with inhibitors Darunavir (DRV) and KNI-1657 (KNI), respectively, we have conducted research on the conformational dynamics and the resistance mechanism caused by residue mutations through multiple molecular dynamics (MD) simulations combined with an energy (MM-PBSA and solvated interaction energy (SIE)) prediction. The results indicate that mutated residues of PR alter the distance between flap regions and catalytic sites, the volume of the inner catalytic site, and the curling degree of the flap tips, thereby affecting DRV and KNI inhibitor binding to PR. These mutated residues reduced the binding affinity of the two mutant PRs to DRV, resulting in drug resistance, whereas the two mutant PRs increase the binding affinity with KNI, indicating they enhance the sensitivity to KNI. Compared with the WT PR, the changes in van der Waals interaction and electrostatic interaction in the two variant PRs play a vital part in the binding of PR with DRV and KNI. These results may supply valuable guidance for the design of anti-AIDS drugs targeting PR.

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.

Molecular Mechanism Exploration of Potent Fluorinated PI3K Inhibitors with a Triazine Scaffold: Unveiling the Unusual Synergistic Effect of Pyridine-to-Pyrimidine Ring Interconversion and CF3 Defluorination

The phosphatidylinostitol-3-kinase (PI3K)/AKT/mammalian target of rapamycin signaling pathway is a vital regulator of cell proliferation, growth, and survival, which is frequently overactivated in many human cancers. To this effect, PI3K, which is an important mediator of this pathway, has been pinpointed as a crucial target in cancer therapy and hence the importance of PI3K inhibitors. It was recently reported that defluorination and pyridine-to-pyrimidine ring interconversion increase the potency of specific small-molecule inhibitors of PI3K. Compound 4, an inhibitor with the difluorinated pyrimidine motif, was found to be eight times more potent against PI3K than compound 1, an inhibitor with the trifluorinated pyridine motif. This observation presents the need to rationally resolve the differential inhibitory mechanisms exhibited by both compounds. In this present work, we employed multiple computational approaches to investigate and distinguish the binding modes of 1 and 4 in addition to the effects they mediate on the secondary structure of PI3K. Likewise, we evaluated two other derivatives, compounds 2 with the difluorinated pyridine motif and 3 with the trifluorinated pyrimidine motif, to investigate the cooperativity effect between the defluorination of CF₃ and pyridine-to-pyrimidine ring interconversion. Findings revealed that PI3K, upon interaction with 4, exhibited a series of structural changes that favored the binding of the inhibitor at the active-site region. Furthermore, a positive (synergistic) cooperativity effect was observed between CF₃ defluorination and pyridine-to-pyrimidine ring interconversion. Moreover, there was a good correlation between the binding free energy estimated and the biological activity reported experimentally. Energy decomposition analysis revealed that the major contributing force to binding affinity variations between 1 and 4 is the electrostatic energy. Per-residue energy-based hierarchical clustering analysis further identified four hot-spot residues ASP841, TYR867, ASP964, and LYS833 and four warm-spot residues ASP836, SER806, ASP837, and LYS808, which essentially mediate the optimal and higher-affinity binding of compound 4 to PI3K relative to 1. This study therefore provides rational insights into the mechanisms by which 4 exhibited superior PI3K-inhibitory activities over 1, which is vital for future structure-based drug discovery efforts in PI3K targeting.

Molecular mechanism related to the binding of fluorophores to Mango-II revealed by multiple-replica molecular dynamics simulations

Recently, RNA aptamers activating small-molecule fluorophores have been successfully applied to tag and track RNAs in vivo. It is of significance to investigate the molecular mechanism of the fluorophore-RNA aptamer bindings at the atomic level to seek a possible pathway to enhance the fluorescence efficiency of fluorophores. In this work, multiple replica molecular dynamics (MRMD) simulations, essential dynamics (ED) analysis, and hierarchical clustering analysis were coupled to probe the effect of A22U mutation on the binding of two fluorophores, TO1-Biotin (TO1) and TO3-Biotin (TO3), to the Mango-II RNA aptamer (Mango-II). ED analysis reveals that A22U induces alterations in the binding pocket and sites of TO1 and TO3 to the Mango-II, which in turn tunes the fluorophore-RNA interface and changes the interactions of TO1 and TO3 with separate nucleotides of Mango-II. Dynamics analyses also uncover that A22U exerts the opposite impact on the molecular surface areas of the Mango-II and sugar puckers of nucleotides 22 and 23 in Mango-II complexed with TO1 and TO3. Moreover, the calculations of binding free energies suggest that A22U strengthens the binding ability of TO1 to the mutated Mango-II but weakens TO3 to the mutated Mango-II when compared with WT. These findings imply that point mutation in nucleotides possibly tune the fluorescence of fluorophores binding to RNA aptamers, providing a possible scheme to enhance the fluorescence of fluorophores.

Insights into effect of the Asp25/Asp25' protonation states on binding of inhibitors Amprenavir and MKP97 to HIV-1 protease using molecular dynamics simulations and MM-GBSA calculations

The protonation states of two aspartic acids in the catalytic strands of HIV-1 protease (PR) remarkably affect bindings of inhibitors to PR. It is requisite for the design of potent inhibitors towards PR to investigate the influences of Asp25/Asp25' protonated states on dynamics behaviour of PR and binding mechanism of inhibitors to PR. In this work, molecular dynamics (MD) simulations, MM-GBSA method and principal component (PC) analysis were coupled to explore the effect of Asp25/Asp25' protonation states on conformational changes of PR and bindings of Amprenavir and MKP97 to PR. The results show that the Asp25/Asp25' protonation states exert different impacts on structural fluctuations, flexibility and motion modes of PR. Dynamics analysis verifies that Asp25/Asp25' protonated states highly affect conformational dynamics of two flaps in PR. The binding free energy calculations results suggest that the Asp25/Asp25' protonated states obviously strengthen bindings of inhibitors to PR compared to the non-protonation state. Calculations of residue-based free energy decomposition indicate that the Asp25/Asp25' protonation not only disturbs the interaction network of inhibitors with PR but also stabilizes bindings of inhibitors to PR by cancelling the electrostatic repulsive interaction. Therefore, special attentions should be paid to the Asp25/Asp25' protonation in the design of potent inhibitors towards PR.

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

To date, inhibiting the activity of β-amyloid cleaving enzyme 1 (BACE1) has been considered an efficient approach for treating Alzheimer's disease (AD). In the current work, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations and the molecular mechanics general Born surface area (MM-GBSA) method were combined to investigate the effect of pH-dependent protonation on the binding of the inhibitors CS9, C6U, and 6WE to BACE1. Dynamic analyses based on the MR-GaMD trajectory show that pH-dependent protonation strongly affects the structural flexibility, correlated motions, and dynamic behavior of inhibitor-bound BACE1. According to the constructed free energy profiles, in the protonated state at low pH, inhibitor-bound BACE1 tends to populate at more conformations than in high pH. The binding free energies calculated by MM-GBSA suggest that inhibitors possess stronger binding abilities under the protonation conditions at high pH than under the protonation conditions at low pH. Moreover, pH-dependent protonation exerts a significant effect on the hydrogen bonding interactions of CS9, C6U, and 6WE to BACE1, which correspondingly alters the binding abilities of the three inhibitors to BACE1. Furthermore, in different protonated environments, three inhibitors share common interaction clusters and similar binding sites in BACE1, which are reliably used as efficient targets for the design of potent inhibitors of BACE1.

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.

Molecular mechanism concerning conformational changes of CDK2 mediated by binding of inhibitors using molecular dynamics simulations and principal component analysis

Cyclin-dependent kinase 2 (CDK2) has been regarded as a promising drug target for anti-tumour agents. In this study, molecular dynamics (MD) simulations and principal component (PC) analysis were used to explore binding mechanism of three inhibitors 1PU, CDK, 50Z to CDK2 and influences of their bindings on conformational changes of CDK2. The results show that bindings of inhibitors yield obvious impacts on internal dynamics, movement patterns and conformational changes of CDK2. In addition, molecular mechanics generalized Born surface area (MM-GBSA) was applied to calculate binding free energies between three inhibitors and CDK2 and evaluate their binding ability to CDK2. The results show that CDK has the strongest binding to CDK2 among the current three inhibitors. Residue-based free energy decomposition method was further utilized to decode the contributions of a single residue to binding of inhibitors, and it was found that three inhibitors not only produce hydrogen bonding interactions and hydrophobic interactions with key residues of CDK2, which promotes binding of three inhibitors to CDK2, but also share similar binding modes. This work is expected to be helpful for design of efficient drugs targeting CDK2.

An Overview of β-Amyloid Cleaving Enzyme 1 (Bace1) in Alzheimer's Disease Therapy Elucidating its Exosite-Binding Antibody and Allosteric Inhibitor

Over decades of its identification, numerous past and ongoing researches have focused on the therapeutic roles of β-amyloid cleaving enzyme 1 (BACE1) as a target in treating Alzheimer's disease (AD). Although the initial BACE1 inhibitors at phase-3 clinical trials tremendously reduced β-amyloid-associated plaques in patients with AD, the researchers eventually discontinued the tests due to the lack of potency. This discontinuation has resulted in limited drug development and discovery targeted at BACE1, despite the high demand for dementia and AD therapies. It is, therefore, imperative to describe the detailed underlying biological basis of the BACE1 therapeutic option in neurological diseases. Herein, we highlight BACE1 bioactivity, genetic properties, and role in neurodegenerative therapy. We review research contributions to BACE1 exosite-binding antibody and allosteric inhibitor development as AD therapies. The review also covers BACE1 biological function, the disease-associated mechanisms, and the enzyme conditions for amyloid precursor protein sites splitting. Based on the present review, we suggest further studies on anti-BACE1 exosite antibodies and BACE1 allosteric inhibitors. Non-active site inhibition might be the way forward to BACE1 therapy in Alzheimer's neurological disorder.

Free energy calculations of ALS-causing SOD1 mutants reveal common perturbations to stability and dynamics along the maturation pathway

With over 150 heritable mutations identified as disease-causative, superoxide dismutase 1 (SOD1) has been a main target of amyotrophic lateral sclerosis (ALS) research and therapeutic efforts. However, recent evidence has suggested that neither loss of function nor protein aggregation is responsible for promoting neurotoxicity. Furthermore, there is no clear pattern to the nature or the location of these mutations that could suggest a molecular mechanism behind SOD1-linked ALS. Here, we utilize reliable and accurate computational techniques to predict the perturbations of 10 such mutations to the free energy changes of SOD1 as it matures from apo monomer to metallated dimer. We find that the free energy perturbations caused by these mutations strongly depend on maturational progress, indicating the need for state-specific therapeutic targeting. We also find that many mutations exhibit similar patterns of perturbation to native and non-native maturation, indicating strong thermodynamic coupling between the dynamics at various sites of maturation within SOD1. These results suggest the presence of an allosteric network in SOD1 which is vulnerable to disruption by these mutations. Analysis of these perturbations may contribute to uncovering a unifying molecular mechanism which explains SOD1-linked ALS and help to guide future therapeutic efforts.

Understanding the Polypharmacological Profiles of Triple Reuptake Inhibitors by Molecular Simulation

The triple reuptake inhibitors (TRIs) class is a class of effective inhibitors of human monoamine transporters (hMATs), which includes dopamine, norepinephrine, and serotonin transporters (hDATs, hNETs, and hSERTs). Due to the high degree of structural homology of the binding sites of those transporters, it is a great challenge to design potent TRIs with fine-tuned binding profiles. The molecular determinants responsible for the binding selectivity of TRIs to hDATs, hNETs, and hSERTs remain elusive. In this study, the solved X-ray crystallographic structure of hSERT in complex with escitalopram was used as a basis for modeling nine complexes of three representative TRIs (SEP225289, NS2359, and EB1020) bound to their corresponding targets. Molecular dynamics (MD) and effective post-trajectory analysis were performed to estimate the drug binding free energies and characterize the selective profiles of each TRI to hMATs. The common binding mode of studied TRIs to hMATs was revealed by hierarchical clustering analysis of the per-residue energy. Furthermore, the combined protein-ligand interaction fingerprint and residue energy contribution analysis indicated that several conserved and nonconserved "Warm Spots" such as S149, V328, and M427 in hDAT, F317, F323, and V325 in hNET and F335, F341, and V343 in hSERT were responsible for the TRI-binding selectivity. These findings provided important information for rational design of a single drug with better polypharmacological profiles through modulating multiple targets.

Binding Selectivity of Inhibitors toward Bromodomains BAZ2A and BAZ2B Uncovered by Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations

Two Bromodomain-Containing proteins BAZ2A and BAZ2B are responsible for remodeling chromatin and regulating noncoding RNAs. As for our current studies, integration of multiple short molecular dynamics simulations (MSMDSs) with molecular mechanics generalized Born surface area (MM-GBSA) method is adopted for insights into binding selectivity of three small molecules D8Q, D9T and UO1 to BAZ2A against BAZ2B. The calculations of MM-GBSA unveil that selectivity of inhibitors toward BAZ2A and BAZ2B highly depends on the enthalpy changes and the details uncover that D8Q has better selectivity toward BAZ2A than BAZ2B, D9T more favorably bind to BAZ2B than BAZ2A, and UO1 does not show obvious selectivity toward these two proteins. The analysis of interaction network between residues and inhibitors indicates that seven residues are mainly responsible for the selectivity of D8Q, six residues for D9T and four residues provide significant contributions to associations of UO1 with two proteins. Moreover the analysis of interaction network not only reveals warm spots of inhibitor bindings to BAZ2A and BAZ2B but also unveils that common residue pairs, including (W1816, W1887), (P1817, P1888), (F1818, F1889), (V1822, V1893), (N1823, N1894),(L1826, L1897), (V1827, V1898), (F1872, F1943), (N1873, N1944) and (V1879, I1950) belonging to (BAZ2A, BAZ2B), induce mainly binding differences of inhibitors to BAZ2A and BAZ2B. Hence, insights from our current studies offer useful dynamics information relating with conformational alterations and structure-affinity relationship at atomistic levels for novel therapeutic strategies toward BAZ2A and BAZ2B.

Binding mechanism of inhibitors to p38α MAP kinase deciphered by using multiple replica Gaussian accelerated molecular dynamics and calculations of binding free energies

The p38α MAP Kinase has been an important target of drug design for treatment of inflammatory diseases and cancers. This work applies multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations and the molecular mechanics generalized Born surface area (MM-GBSA) method to probe the binding mechanism of inhibitors L51, R24 and 1AU to p38α. Dynamics analyses show that inhibitor bindings exert significant effect on conformational changes of the active helix α2 and the conserved DFG loop. The rank of binding free energies calculated with MM-GBSA not only agrees well with that determined by the experimental IC50 values but also suggests that mutual compensation between the enthalpy and entropy changes can improve binding of inhibitors to p38α. The analyses of free energy landscapes indicate that the L51, R24 and 1AU bound p38α display a DFG-out conformation. The residue-based free energy decomposition method is used to evaluate contributions of separate residues to the inhibitor-p38α binding and the results imply that residues V30, V38, L74, L75, I84, T106, H107, L108, M109, L167, F169 and D168 can be utilized as efficient targets of potent inhibitors toward p38α.

Screening Platform Based on Inductively Coupled Plasma Mass Spectrometry for β-Site Amyloid Protein Cleaving Enzyme 1 (BACE1) Inhibitors

β-Site amyloid protein cleaving enzyme 1 (BACE1) is a promising therapeutic target for developing inhibitors to alleviate Alzheimer's disease (AD). Herein, we established an inductively coupled plasma mass spectrometry (ICPMS)-based inhibitor screening platform. A biotin-labeled lanthanide-coded peptide probe (LCPP; biotin-PEG₂-EVNLDAEC-DOTA-Ln) was designed to determine the activity of BACE1 and evaluate the degree of inhibition of inhibitors. The platform was first validated with two commercially available inhibitors (BSI I and BSI IV) in terms of IC₅₀ values and then applied to two newly designed inhibitors (inhibitors II and III) based on the crystal structure of BACE1 interacting with inhibitor I, and each of them contained an acylguanidine core structure. We found that their inhibition effects were improved as evaluated by the sensitive and accurate LCPP-ICPMS platform, demonstrating its ability for new drug screening.

Molecular mechanisms of inhibitor bindings to A-FABP deciphered by using molecular dynamics simulations and calculations of MM-GBSA

Adipocyte fatty-acid binding protein (A-FABP) plays a central role in many aspects of metabolic diseases. It is an important target in drug design for treatment of FABP-related diseases. In this study, molecular dynamics (MD) simulations followed by calculations of molecular mechanics generalized Born surface area (MM-GBSA) and principal components analysis (PCA) were implemented to decipher molecular mechanism correlating with binding of inhibitors 57Q, 57P and L96 to A-FABP. The results show that van der Waals interactions are the leading factors to control associations of 57Q, 57P, and L96 with A-FABP, which reveals an energetic basis for designing of clinically available inhibitors towards A-FABP. The information from PCA and cross-correlation analysis rationally unveils that inhibitor bindings affect conformational changes of A-FABP and change relative movements between residues. Decomposition of binding affinity into contributions of individual residues not only detects hot spots of inhibitor/A-FABP binding but also shows that polar interactions of the positively charged residue Arg126 with three inhibitors provide a significant contribution for stabilization of the inhibitor/A-FABP bindings. Furthermore, the binding strength of L96 to residues Ser55, Phe57 and Lys58 are stronger than that of inhibitors 57Q and 57P to these residues.

Mutation-Induced Impacts on the Switch Transformations of the GDP- and GTP-Bound K-Ras: Insights from Multiple Replica Gaussian Accelerated Molecular Dynamics and Free Energy Analysis

Mutations yield significant effect on the structural flexibility of two switch domains, SW1 and SW2, in K-Ras, which is considered as an important target of anticancer drug design. To unveil a molecular mechanism with regard to mutation-mediated tuning on the activity of K-Ras, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations followed by analysis of free energy landscapes (FELs) are performed on the GDP- and GTP-bound wild-type (WT), G12V, and D33E K-Ras. The results suggest that G12V and D33E not only evidently change the flexibility of SW1 and SW2 but also greatly affect correlated motions of SW1 and SW2 separately relative to the P-loop and SW1, which exerts a certain tuning on the activity of K-Ras. The information stemming from the analyses of FELs reveals that the conformations of SW1 and SW2 are in high disorders in the GDP- and GTP-associated WT and mutated K-Ras, possibly producing significant effect on binding of guanine nucleotide exchange factors or effectors to K-Ras. The interaction networks of GDP and GTP with K-Ras are identified and the results uncover that the instability in hydrogen-bonding interactions of SW1 with GDP and GTP is mostly responsible for conformational disorder of SW1 and SW2 as well as tunes the activity of oncogenic K-Ras.

Sensitivity of Binding Free Energy Calculations to Initial Protein Crystal Structure

Binding free energy calculations using alchemical free energy (AFE) methods are widely considered to be the most rigorous tool in the computational drug discovery arsenal. Despite this, the calculations suffer from accuracy, precision, and reproducibility issues. In this publication, we perform a high-throughput study of more than a thousand AFE calculations, utilizing over 220 μs of total sampling time, on three different protein systems to investigate the impact of the initial crystal structure on the resulting binding free energy values. We also consider the influence of equilibration time and discover that the initial crystal structure can have a significant effect on free energy values obtained at short timescales that can manifest itself as a free energy difference of more than 1 kcal/mol. At longer timescales, these differences are largely overtaken by important rare events, such as torsional ligand motions, typically resulting in a much higher uncertainty in the obtained values. This work emphasizes the importance of rare event sampling and long-timescale dynamics in free energy calculations even for routinely performed alchemical perturbations. We conclude that an optimal protocol should not only concentrate computational resources on achieving convergence in the alchemical coupling parameter (λ) space but also on longer simulations and multiple repeats.

Insights into interaction mechanism of inhibitors E3T, E3H and E3B with CREB binding protein by using molecular dynamics simulations and MM-GBSA calculations

CREB binding protein (CBP) and its paralog E1A binding protein (p300) are related to the development of inflammatory diseases, cancers and other diseases, and have been potential targets for the treatment of human diseases. In this work, interaction mechanism of three small molecules E3T, E3H, and E3B with CBP was investigated by employing molecular dynamics (MD) simulations, principal component analysis (PCA), and molecular mechanics/generalized born surface area (MM-GBSA) method. The results indicate that inhibitor bindings cause the changes of movement modes and structural flexibility of CBP, and van der Waals interactions mostly drive associations of inhibitors with CBP. In the meantime, the results based on inhibitor-residue interactions not only show that eight residues of CBP can strongly interact with E3T, E3H and E3B but also verify that the CH-CH, CH-π, and π-π interactions are responsible for vital contributions in associations of E3T, E3H and E3B with CBP. In addition, the H-O radial distribution functions (RDFs) were computed to assess the stability of hydrogen bonding interactions between inhibitors and CBP, and the obtained information identifies several key hydrogen bonds playing key roles in bindings of E3T, E3H and E3B to CBP. Potential new inhibitors have been proposed.

A Structure-Based Discovery Platform for BACE2 and the Development of Selective BACE Inhibitors

The ability to perform routine structure-guided drug design for selective BACE inhibitors has been limited because of the lack of robust platform for BACE2 expression, purification, and crystallization. To overcome this limitation, we developed a platform that produces 2-3 mg of pure BACE2 protein per liter of E. coli culture, and we used this protein to design macrocyclic compounds that potently and selectively inhibit BACE1 over BACE2. Compound 2 was found to potently inhibit BACE 1 (K_i = 5 nM) with a selectivity of 214-fold over BACE2. The X-ray crystal structures of unbound BACE2 (2.2 Å) and BACE2 bound to compound 3 (3.0 Å and K_i = 7 nM) were determined and compared to the X-ray structures of BACE1 revealing the S1-S3 subsite as a selectivity determinant. This platform should enable a more rapid development of new and selective BACE inhibitors for the treatment of Alzheimer's disease or type II diabetes.

Identification of new BACE1 inhibitors for treating Alzheimer's disease

Alzheimer's disease (AD) is a type of brain disorder, wherein a person experiences gradual memory loss, state of confusion, hallucination, agitation, and personality change. AD is marked by the presence of extracellular amyloid plaques and intracellular neurofibrillary tangles (NFTs) and synaptic losses. Increased cases of AD in recent times created a dire need to discover or identify chemical compounds that can cease the development of AD. This study focuses on finding potential drug molecule(s) active against β-secretase, also known as β-site amyloid precursor protein cleaving enzyme 1 (BACE1). Clustering analysis followed by phylogenetic studies on microarray datasets retrieved from GEO browser showed that BACE1 gene has genetic relatedness with the RCAN1 gene. A ligand library comprising 60 natural compounds retrieved from literature and 25 synthetic compounds collected from DrugBank were screened. Further, 350 analogues of potential parent compounds were added to the library for the docking purposes. Molecular docking studies identified 11-oxotigogenin as the best ligand molecule. The compound showed the binding affinity of - 11.1 Kcal/mole and forms three hydrogen bonds with Trp124, Ile174, and Arg176. The protein-ligand complex was subjected to 25 ns molecular dynamics simulation and the potential energy of the complex was found to be - 1.24579e+06 Kcal/mole. In this study, 11-oxotigogenin has shown promising results against BACE1, which is a leading cause of AD, hence warrants for in vitro and in vivo validation of the same. In addition, in silico identification of 11-oxotigogenin as a potential anti-AD compound paves the way for designing of chemical scaffolds to discover more potent BACE1 inhibitors.Graphical abstract.

Unveiling conformational dynamics changes of H-Ras induced by mutations based on accelerated molecular dynamics

Uncovering molecular basis with regard to the conformational change of two switches I and II in the GppNHp (GNP)-bound H-Ras is highly significant for the understanding of Ras signaling. For this purpose, accelerated molecular dynamics (aMD) simulations and principal component (PC) analysis are integrated to probe the effect of mutations G12V, T35S and Q61K on conformational transformation between two switches of the GNP-bound H-Ras. The RMSF and cross-correlation analyses suggest that three mutations exert a vital effect on the flexibility and internal dynamics of two switches in the GNP-bound H-Ras. The results stemming from PC analysis indicate that two switches in the GNP-bound WT H-Ras tend to form a closed state in most conformations, while those in the GNP-bound mutated H-Ras display transformation between different states. This conclusion is further supported by free energy landscapes constructed by using the distances of residues 12 away from 35 and 35 away from 61 as reaction coordinates and different experimental studies. Interaction scanning is performed on aMD trajectories and the information shows that conformational transformations of two switches I and II induced by mutations extremely affect the GNP-residue interactions. Meanwhile, the scanning results also signify that residues G15, A18, F28, K117, A146 and K147 form stable contacts with GNP, while residues D30, E31, Y32, D33, P34 and E62 in two switches I and II produce unstable contacts with GNP. This study not only reveals dynamic behavior changes of two switches in H-Ras induced by mutations, but also unveils general principles and mechanisms with regard to functional conformational changes of H-Ras.

Binding selectivity of inhibitors toward the first over the second bromodomain of BRD4: theoretical insights from free energy calculations and multiple short molecular dynamics simulations

Bromodomain-containing protein 4 (BRD4) plays an important role in mediating gene transcription involved in cancers and non-cancer diseases such as acute heart failure and inflammatory diseases. In this work, multiple short molecular dynamics (MSMD) simulations are integrated with a molecular mechanics generalized Born surface area (MM-GBSA) approach to decipher binding selectivity of three inhibitors 8NS, 82Y, and 837 toward two domains BD1 and BD2 of BRD4. The results demonstrate that the enthalpy effects play critical roles in selectivity identification of inhibitors toward BD1 and BD2, determining that 8NS has better selectivity toward BD2 than BD1, while 82Y and 837 more favorably bind to BD1 than BD2. A residue-based free-energy decomposition method was used to calculate an inhibitor-residue interaction spectrum and unveil contributions of separate residues to binding selectivity. The results identify six common residues, containing (P82, P375), (V87, V380), (L92, L385), (L94, L387), (N140, N433), and (I146, V439) individually belonging to (BD1, BD2) of BRD4, and yield a considerable binding difference of inhibitors to BD1 and BD2, suggesting that these residues play key roles in binding selectivity of inhibitors toward BD1 and BD2 of BRD4. Therefore, these results provide useful dynamics information and a structure affinity relationship for the development of highly selective inhibitors targeting BD1 and BD2 of BRD4.

Mechanistic Understanding From Molecular Dynamics Simulation in Pharmaceutical Research 1: Drug Delivery

In this review, we outline the growing role that molecular dynamics simulation is able to play as a design tool in drug delivery. We cover both the pharmaceutical and computational backgrounds, in a pedagogical fashion, as this review is designed to be equally accessible to pharmaceutical researchers interested in what this new computational tool is capable of and experts in molecular modeling who wish to pursue pharmaceutical applications as a context for their research. The field has become too broad for us to concisely describe all work that has been carried out; many comprehensive reviews on subtopics of this area are cited. We discuss the insight molecular dynamics modeling has provided in dissolution and solubility, however, the majority of the discussion is focused on nanomedicine: the development of nanoscale drug delivery vehicles. Here we focus on three areas where molecular dynamics modeling has had a particularly strong impact: (1) behavior in the bloodstream and protective polymer corona, (2) Drug loading and controlled release, and (3) Nanoparticle interaction with both model and biological membranes. We conclude with some thoughts on the role that molecular dynamics simulation can grow to play in the development of new drug delivery systems.

Multiple Molecular Dynamics Simulations of the Inhibitor GRL-02031 Complex with Wild Type and Mutant HIV-1 Protease Reveal the Binding and Drug-Resistance Mechanism

Human immunodeficiency virus type 1 (HIV-1) protease is regarded as a fascinating target for drug development against HIV infection. However, mutations causing drug resistance severely limit the efficiency of the recently marketed drugs in the treatment of HIV replication. To elucidate the binding mechanism of HIV-1 protease with promising inhibitor GRL-02031 and further to probe the resistance mechanism associated with mutations (I47V, L76V, V82A, and N88D) to the inhibitor, we applied multiple molecular dynamics (MMD) simulations along with energy analysis by the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and solvated interaction energy (SIE) methodology on specific HIV-1 protease with GRL-0231 complexes. On the basis of detail analysis of the simulations, we revealed key characteristics that constitute the drug resistance of four mutation HIV-1 proteases toward GRL-02031: substitution of the side chain in these four mutation residues leads to a change in the distances between the flaps and catalytic sites, thereby reducing the affinity for GRL-02031 with these four mutation proteases, even though the L76V and N88D residues cannot directly contact GRL-02031. The results of energy analysis according to the MM-PBSA and SIE methods further indicated that hydrophobic interaction was considered to be the prime driving force for inhibitor GRL-02031 binding to protease and the decrease in van der Waals interactions between inhibitor GRL-02031 and mutant proteases as the primary cause of the drug resistance. Analyses of the hydrogen bonds and atomic interactions further provided detailed explanations for the resistance of these four mutation proteases toward inhibitor GRL-02031. The present study provides potential guidance on the structure-based inhibitors' design targeting HIV-1 protease.