Computational study of protein specificity: the molecular basis of HIV-1 protease drug resistance

Binding Mechanism of Inhibitors to BRD4 and BRD9 Decoded by Multiple Independent Molecular Dynamics Simulations and Deep Learning

Bromodomain 4 and 9 (BRD4 and BRD9) have been regarded as important targets of drug designs in regard to the treatment of multiple diseases. In our current study, molecular dynamics (MD) simulations, deep learning (DL) and binding free energy calculations are integrated to probe the binding modes of three inhibitors (H1B, JQ1 and TVU) to BRD4 and BRD9. The MD trajectory-based DL successfully identify significant functional function domains, such as BC-loop and ZA-loop. The information from the post-processing analysis of MD simulations indicates that inhibitor binding highly influences the structural flexibility and dynamic behavior of BRD4 and BRD9. The results of the MM-GBSA calculations not only suggest that the binding ability of H1B, JQ1 and TVU to BRD9 are stronger than to BRD4, but they also verify that van der Walls interactions are the primary forces responsible for inhibitor binding. The hot spots of BRD4 and BRD9 revealed by residue-based free energy estimation provide target sites of drug design in regard to BRD4 and BRD9. This work is anticipated to provide useful theoretical aids for the development of selective inhibitors over BRD family members.

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.

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.

Identification of potential inhibitors targeting yellow fever virus helicase through ligand and structure-based computational studies

Yellow fever is a flavivirus having plus-sensed RNA which encodes a single polyprotein. Host proteases cut this polyprotein into seven nonstructural proteins including a vital NS3 protein. The present study aims to identify the most effective inhibitor against the helicase (NS3) using different advanced ligand and structure-based computational studies. A set of 300 ligands was selected against helicase by chemical structural similarity model, which are similar to S-adenosyl-l-cysteine using infiniSee. This tool screens billions of compounds through a similarity search from in-built chemical spaces (CHEMriya, Galaxi, KnowledgeSpace and REALSpace). The pharmacophore was designed from ligands in the library that showed same features. According to the sequence of ligands, six compounds (29, 87, 99, 116, 148, and 208) were taken for pharmacophore designing against helicase protein. Subsequently, compounds from the library which showed the best pharmacophore shared-features were docked using FlexX functionality of SeeSAR and their optibrium properties were analyzed. Afterward, their ADME was improved by replacing the unfavorable fragments, which resulted in the generation of new compounds. The selected best compounds (301, 302, 303 and 304) were docked using SeeSAR and their pharmacokinetics and toxicological properties were evaluated using SwissADME. The optimal inhibitor for yellow fever helicase was 2-amino-N-(4-(dimethylamino)thiazol-2-yl)-4-methyloxazole-5-carboxamide (302), which exhibits promising potential for drug development.Communicated by Ramaswamy H. Sarma.

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.

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.

Searching for potential inhibitors of SARS-COV-2 main protease using supervised learning and perturbation calculations

Inhibiting the biological activity of SARS-CoV-2 Mpro can prevent viral replication. In this context, a hybrid approach using knowledge- and physics-based methods was proposed to characterize potential inhibitors for SARS-CoV-2 Mpro. Initially, supervised machine learning (ML) models were trained to predict a ligand-binding affinity of ca. 2 million compounds with the correlation on a test set of R=0.748±0.044. Atomistic simulations were then used to refine the outcome of the ML model. Using LIE/FEP calculations, nine compounds from the top 100 ML inhibitors were suggested to bind well to the protease with the domination of van der Waals interactions. Furthermore, the binding affinity of these compounds is also higher than that of nirmatrelvir, which was recently approved by the US FDA to treat COVID-19. In addition, the ligands altered the catalytic triad Cys145 - His41 - Asp187, possibly disturbing the biological activity of SARS-CoV-2.

Finding potential inhibitors against RNA-dependent RNA polymerase (RdRp) of bovine ephemeral fever virus (BEFV): an in-silico study

The bovine ephemeral fever virus (BEFV) is an enzootic agent that affects millions of bovines and causes major economic losses. Though the virus is seasonally reported with a very high morbidity rate (80-100%) from African, Australian, and Asiatic continents, it remains a neglected pathogen in many of its endemic areas, with no proper therapeutic drugs or vaccines presently available for treatment. The RNA-dependent RNA polymerase (RdRp) catalyzes the viral RNA synthesis and is an appropriate candidate for antiviral drug developments. We utilized integrated computational tools to build the 3D model of BEFV-RdRp and then predicted its probable active binding sites. The virtual screening and optimization against these active sites, using several small-molecule inhibitors from a different category of Life Chemical database and FDA-approved drugs from the ZINC database, was performed. We found nine molecules that have docking scores varying between -6.84 to -10.43 kcal/mol. Furthermore, these complexes were analyzed for their conformational dynamics and thermodynamic stability using molecular dynamics simulations in conjunction with the molecular mechanics generalized Born surface area (MM-GBSA) scheme. The binding free energy calculations depict that the electrostatic interactions play a dominant role in the RdRp-inhibitor binding. The hot spot residues, such as Arg565, Asp631, Glu633, Asp740, and Glu707, were found to control the RdRp-inhibitor interaction. The ADMET analysis strongly suggests favorable pharmacokinetics of these compounds that may prove useful for treating the BEFV ailment. Overall, we anticipate that these findings would help explore and develop a wide range of anti-BEFV therapy.Communicated by Ramaswamy H. Sarma.

The evolution of the HIV-1 protease folding stability

The evolution of structural proteins is generally constrained by the folding stability. However, little is known about the particular capacity of viral proteins to accommodate mutations that can potentially affect the protein stability and, in general, the evolution of the protein stability over time. As an illustrative model case, here, we investigated the evolution of the stability of the human immunodeficiency virus (HIV-1) protease (PR), which is a common HIV-1 drug target, under diverse evolutionary scenarios that include (1) intra-host virus evolution in a cohort of seventy-five patients sampled over time, (2) intra-host virus evolution sampled before and after specific PR-based treatments, and (3) inter-host evolution considering extant and ancestral (reconstructed) PR sequences from diverse HIV-1 subtypes. We also investigated the specific influence of currently known HIV-1 PR resistance mutations on the PR folding stability. We found that the HIV-1 PR stability fluctuated over time within a constant and wide range in any studied evolutionary scenario, accommodating multiple mutations that partially affected the stability while maintaining activity. We did not identify relationships between change of PR stability and diverse clinical parameters such as viral load, CD4⁺ T-cell counts, and a surrogate of time from infection. Counterintuitively, we predicted that nearly half of the studied HIV-1 PR resistance mutations do not significantly decrease stability, which, together with compensatory mutations, would allow the protein to adapt without requiring dramatic stability changes. We conclude that the HIV-1 PR presents a wide structural plasticity to acquire molecular adaptations without affecting the overall evolution of stability.

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.

Identifying Possible AChE Inhibitors from Drug-like Molecules via Machine Learning and Experimental Studies

Acetylcholinesterase (AChE) is one of the most important drug targets for Alzheimer's disease (AD) treatment. In this work, a machine learning model was trained to rapidly and accurately screen large chemical databases for the potential inhibitors of AChE. The obtained results were then validated via in vitro enzyme assay. Moreover, atomistic simulations including molecular docking and molecular dynamics simulations were then used to understand molecular insights into the binding process of ligands to AChE. In particular, two compounds including benzyl trifluoromethyl ketone and trifluoromethylstyryl ketone were indicated as highly potent inhibitors of AChE because they established IC₅₀ values of 0.51 and 0.33 μM, respectively. The obtained IC₅₀ of two compounds is significantly lower than that of galantamine (2.10 μM). The predicted log(BB) suggests that the compounds may be able to traverse the blood-brain barrier. A good agreement between computational and experimental studies was observed, indicating that the hybrid approach can enhance AD therapy.

Losartan Interactions with 2-Hydroxypropyl-β-CD

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

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.

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.

Using a comprehensive approach to investigate the interaction between Kinesin-5/Eg5 and the microtubule

Kinesins are microtubule-based motor proteins that play important roles ranging from intracellular transport to cell division. Human Kinesin-5 (Eg5) is essential for mitotic spindle assembly during cell division. By combining molecular dynamics (MD) simulations with other multi-scale computational approaches, we systematically studied the interaction between Eg5 and the microtubule. We find the electrostatic feature on the motor domains of Eg5 provides attractive interactions to the microtubule. Additionally, the folding and binding energy analysis reveals that the Eg5 motor domain performs its functions best when in a weak acidic environment. Molecular dynamics analyses of hydrogen bonds and salt bridges demonstrate that, on the binding interfaces of Eg5 and the tubulin heterodimer, salt bridges play the most significant role in holding the complex. The salt bridge residues on the binding interface of Eg5 are mostly positive, while salt bridge residues on the binding interface of tubulin heterodimer are mostly negative. Such salt bridge residue distribution is consistent with electrostatic potential calculations. In contrast, the interface between α and β-tubulins is dominated by hydrogen bonds rather than salt bridges. Compared to the Eg5/α-tubulin interface, the Eg5/β-tubulin interface has a greater number of salt bridges and higher occupancy for salt bridges. This asymmetric salt bridge distribution may play a significant role in Eg5′s directionality. The residues involved in hydrogen bonds and salt bridges are identified in this work and may be helpful for anticancer drug design.

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.

Investigating potency of TMC-126 against wild-type and mutant variants of HIV-1 protease: a molecular dynamics and free energy study

A detailed computational study was performed to investigate the conformational changes of flap region and the mechanism underlying the binding of the inhibitor TMC-126 to HIV-1 protease (PR1) and its mutant variants through molecular dynamics simulations in conjunction with the molecular mechanics Poisson-Boltzmann (MM-PBSA) free energy calculation. Further, we have studied the effectiveness of the inhibitor against HIV-2 protease (PR2). The MM-PBSA calculation suggests that TMC-126 loses its potency against mutant variants and PR2 compared to wild-type PR1 mainly due to the loss in intermolecular electrostatic interactions. The potency of the inhibitor decreases in the order: wild type PR1 > M46L > MDR20 > I50V > PR2 > V32I > A28S. Our study reveals that the flap of PR1 adopts a semi-open conformation due to the mutation I50V or MDR20. The dissimilar nature of the movement of the flap tip of both monomers is evident from the dynamic cross-correlation map. The protein structural network analysis displays that mutation causes structural rearrangements and changes the communication path between residues. Overall, we believe our study may help explore and accelerate the development of novel HIV-1/HIV-2 protease inhibitors with better potency.

Identification of Potential Inhibitors against Epstein-Barr Virus Nuclear Antigen 1 (EBNA1): An Insight from Docking and Molecular Dynamic Simulations

Epstein-Barr virus (EBV), a known tumorigenic virus, is associated with various neuropathies, including multiple sclerosis (MS). However, there is no anti-EBV FDA-approved drug available in the market. Our study targeted EBV protein EBV nuclear antigen 1 (EBNA1), crucial in virus replication and expressed in all the stages of viral latencies. This dimeric protein binds to an 18 bp palindromic DNA sequence and initiates the process of viral replication. We chose phytochemicals and FDA-approved MS drugs based on literature survey followed by their evaluation efficacies as anti-EBNA1 molecules. Molecular docking revealed FDA drugs ozanimod, siponimod, teriflunomide, and phytochemicals; emodin; protoapigenone; and EGCG bound to EBNA1 with high affinities. ADMET and Lipinski's rule analysis of the phytochemicals predicted favorable druggability. We supported our assessments of pocket druggability with molecular dynamics simulations and binding affinity predictions by the molecular mechanics generalized Born surface area (MM/GBSA) method. Our results establish a stable binding for siponimod and ozanimod with EBNA1 mainly via van der Waals interactions. We identified hot spot residues like I481', K477', L582', and K586' in the binding of ligands. In particular, K477' at the amino terminal of EBNA1 is known to establish interaction with two bases at the major groove of the DNA. Siponimod bound to EBNA1 engaging K477', thus plausibly making it unavailable for DNA interaction. Computational alanine scanning further supported the significant roles of K477', I481', and K586' in the binding of ligands with EBNA1. Conclusively, the compounds showed promising results to be used against EBNA1.

Structural determinants for subnanomolar inhibition of the secreted aspartic protease Sapp1p from Candida parapsilosis

Pathogenic Candida albicans yeasts frequently cause infections in hospitals. Antifungal drugs lose effectiveness due to other Candida species and resistance. New medications are thus required. Secreted aspartic protease of C. parapsilosis (Sapp1p) is a promising target. We have thus solved the crystal structures of Sapp1p complexed to four peptidomimetic inhibitors. Three potent inhibitors (K_i: 0.1, 0.4, 6.6 nM) resembled pepstatin A (K_i: 0.3 nM), a general aspartic protease inhibitor, in terms of their interactions with Sapp1p. However, the weaker inhibitor (K_i: 14.6 nM) formed fewer nonpolar contacts with Sapp1p, similarly to the smaller HIV protease inhibitor ritonavir (K_i: 1.9 µM), which, moreover, formed fewer H-bonds. The analyses have revealed the structural determinants of the subnanomolar inhibition of C. parapsilosis aspartic protease. Because of the high similarity between Saps from different Candida species, these results can further be used for the design of potent and specific Sap inhibitor-based antimycotic drugs.

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.

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

Molecular dynamics simulations data of six compounds F3J-BRD4/CBP, EX1-BRD4/CBP, and E2T-BRD4/CBP

The data here described are related to the research article entitled “Molecular dynamics insights into binding selectivity of inhibitors toward BRD4 and CBP” [1]. Bromodomain-containing protein 4 (BRD4) and CREB binding protein (CBP) play important roles in tumorigenesis and development. We performed 200-ns molecular dynamics (MD) simulations on three pairs of inhibitor-BRD4 and inhibitor-CBP complexes to clarify binding selectivity of inhibitors toward BRD4 and CBP. Principal component (PC) analysis was used to probe changes in internal dynamics and conformations of BRD4 and CBP due to inhibitor bindings. Analysis of residue-based free energy decomposition was employed to explore the roles of separate residues in binding selectivity of inhibitors to BRD4 versus CBP.

Characterizing an allosteric inhibitor-induced inactive state in with-no-lysine kinase 1 using Gaussian accelerated molecular dynamics simulations

The binding of an allosteric inhibitor in WNK1 leads to the inactive state.

Molecular Dynamics Protocols for the Study of Cyclodextrin Drug Delivery Systems

Hypertension treatment is a current therapeutic priority as there is a constantly increasing part of the population that suffers from this risk factor, which may lead to cardiovascular and encephalic episodes and eventually to death. A number of marketed medicines consist of active ingredients that may be relatively potent; however, there is plenty of room to enhance their pharmacological profile and therapeutic index by improving specific physicochemical properties. In this work, we focus on a class of blood pressure regulators, called sartans, and we present the computational scheme for the pharmacological improvement of irbesartan (IRB) as a representative example. IRB has been shown to exert increased pharmacological action compared with other sartans, but it appears to be highly lipophilic and violates Lipinski rule (MLogP >4.15). To circumvent this drawback, proper hydrophilic molecules, such as cyclodextrins, can be used as drug carriers. This chapter describes the combinatory use of computational methods, namely molecular docking, quantum mechanics, molecular dynamics, and free energy calculations, to study the interactions and the energetic contributions that govern the IRB:cyclodextrin association. We provide a detailed computational protocol, which aims to assist the improvement of the pharmacological properties of sartans. This protocol can also be applied to any other drug molecule with diminished hydrophilic character.

Accurate Prediction of Inhibitor Binding to HIV-1 Protease Using CANDOCK

The human immunodeficiency virus 1 (HIV-1) protease is an important target for treating HIV infection. Our goal was to benchmark a novel molecular docking protocol and determine its effectiveness as a therapeutic repurposing tool by predicting inhibitor potency to this target. To accomplish this, we predicted the relative binding scores of various inhibitors of the protease using CANDOCK, a hierarchical fragment-based docking protocol with a knowledge-based scoring function. We first used a set of 30 HIV-1 protease complexes as an initial benchmark to optimize the parameters for CANDOCK. We then compared the results from CANDOCK to two other popular molecular docking protocols Autodock Vina and Smina. Our results showed that CANDOCK is superior to both of these protocols in terms of correlating predicted binding scores to experimental binding affinities with a Pearson coefficient of 0.62 compared to 0.48 and 0.49 for Vina and Smina, respectively. We further leveraged the Database of Useful Decoys: Enhanced (DUD-E) HIV protease set to ascertain the effectiveness of each protocol in discriminating active versus decoy ligands for proteases. CANDOCK again displayed better efficacy over the other commonly used molecular docking protocols with area under the receiver operating characteristic curve (AUROC) of 0.94 compared to 0.71 and 0.74 for Vina and Smina. These findings support the utility of CANDOCK to help discover novel therapeutics that effectively inhibit HIV-1 and possibly other retroviral proteases.

Computational Determination of Potential Inhibitors of SARS-CoV-2 Main Protease

The novel coronavirus (SARS-CoV-2) has infected several million people and caused thousands of deaths worldwide since December 2019. As the disease is spreading rapidly all over the world, it is urgent to find effective drugs to treat the virus. The main protease (Mpro) of SARS-CoV-2 is one of the potential drug targets. Therefore, in this context, we used rigorous computational methods, including molecular docking, fast pulling of ligand (FPL), and free energy perturbation (FEP), to investigate potential inhibitors of SARS-CoV-2 Mpro. We first tested our approach with three reported inhibitors of SARS-CoV-2 Mpro, and our computational results are in good agreement with the respective experimental data. Subsequently, we applied our approach on a database of ∼4600 natural compounds, as well as 8 available HIV-1 protease (PR) inhibitors and an aza-peptide epoxide. Molecular docking resulted in a short list of 35 natural compounds, which was subsequently refined using the FPL scheme. FPL simulations resulted in five potential inhibitors, including three natural compounds and two available HIV-1 PR inhibitors. Finally, FEP, the most accurate and precise method, was used to determine the absolute binding free energy of these five compounds. FEP results indicate that two natural compounds, cannabisin A and isoacteoside, and an HIV-1 PR inhibitor, darunavir, exhibit a large binding free energy to SARS-CoV-2 Mpro, which is larger than that of 13b, the most reliable SARS-CoV-2 Mpro inhibitor recently reported. The binding free energy largely arises from van der Waals interaction. We also found that Glu166 forms H-bonds to all of the inhibitors. Replacing Glu166 by an alanine residue leads to ∼2.0 kcal/mol decreases in the affinity of darunavir to SARS-CoV-2 Mpro. Our results could contribute to the development of potential drugs inhibiting SARS-CoV-2.

Prediction of Drug Potencies of BACE1 Inhibitors: A Molecular Dynamics Simulation and MM_GB(PB)SA Scoring

Alzheimer’s disease (AD) is a progressive neurodegenerative brain disorder. One of the important therapeutic approaches of AD is the inhibition of β-site APP cleaving enzyme-1 (BACE1). This enzyme plays a central role in the synthesis of the pathogenic β-amyloid peptides (Aβ) in Alzheimer’s disease. A group of potent BACE1 inhibitors with known X-ray structures (PDB ID 5i3X, 5i3Y, 5iE1, 5i3V, 5i3W, 4LC7, 3TPP) were studied by molecular dynamics simulation and binding energy calculation employing MM_GB(PB)SA. The calculated binding energies gave Kd values of 0.139 µM, 1.39 nM, 4.39 mM, 24.3 nM, 1.39 mM, 29.13 mM, and 193.07 nM, respectively. These inhibitors showed potent inhibitory activities in enzymatic and cell assays. The Kd values are compared with experimental values and the structures are discussed in view of the energy contributions to binding. Drug likeness of these inhibitors is also discussed. Accommodation of ligands in the catalytic site of BACE1 is discussed depending on the type of fragment involved in each structure. Molecular dynamics (MD) simulations and energy studies were used to explore the recognition of the selected BACE1 inhibitors by Asp32, Asp228, and the hydrophobic flap. The results show that selective BACE1 inhibition may be due to the formation of strong electrostatic interactions with Asp32 and Asp228 and a large number of hydrogen bonds, in addition to π–π and van der Waals interactions with the amino acid residues located inside the catalytic cavity. Interactions with the ligands show a similar binding mode with BACE1. These results help to rationalize the design of selective BACE1 inhibitors.

Discovery of a stable tripeptide targeting the N-domain of CRF1 receptor

The corticotropin-releasing factor (CRF) and its CRF1 receptor (CRF₁R) play a central role in the maintenance of homeostasis. Malfunctioning of the CRF/CRF₁R unit is associated with several disorders, such as anxiety and depression. Non-peptide CRF₁R-selective antagonists have been shown to exert anxiolytic and antidepressant effects on experimental animals. However, none of them is in clinical use today because of several side effects, thus demonstrating the need for the development of other more suitable CRF₁R antagonists. In an effort to develop novel CRF₁R antagonists we designed, synthesized and chemically characterized two tripeptide analogues of CRF, namely (R)-LMI and (S)-LMI, having their Leu either in R (or D) or in S (or L) configuration, respectively. Their design was based on the crystal structure of the N-extracellular domain (N-domain) of CRF₁R/CRF complex, using a relevant array of computational methods. Experimental evaluation of the stability of synthetic peptides in human plasma has revealed that (R)-LMI is proteolytically more stable than (S)-LMI. Based on this finding, (R)-LMI was selected for pharmacological characterization. We have found that (R)-LMI is a CRF antagonist, inhibiting (1) the CRF-stimulated accumulation of cAMP in HEK 293 cells expressing the CRF₁R, (2) the production of interleukins by adipocytes and (3) the proliferation rate of RAW 264.7 cells. (R)-LMI likely blocked agonist actions by interacting with the N-domain of CRF₁R as suggested by data using a constitutively active chimera of CRF₁R. We propose that (R)-LMI can be used as an optimal lead compound in the rational design of novel CRF antagonists.

Antihypertensive activity and molecular interactions of irbesartan in complex with 2-hydroxypropyl-β-cyclodextrin

Irbesartan (IRB) exerts beneficial effects either alone or in combination with other drugs on numerous diseases, such as cancer, diabetes, and hypertension. However, due to its high lipophilicity, IRB does not possess the optimum pharmacological efficiency. To circumvent this problem, a drug delivery system with 2-hydroxypropyl-β-cyclodextrin (2-HP-β-CD) was explored. The 1:1 complex between IRB and 2-HP-β-CD was identified through ESI QTF HRMS. Dissolution studies showed a higher dissolution rate of the lyophilized IRB-2-HP-β-CD complex than the tablet containing IRB at pH = 1.2. DSC results revealed the differences of the thermal properties between the complex and various mixtures consisting of the two components, namely IRB and 2-HP-β-CD. Interestingly, depending on the way the mixture preparation was conducted, different association between the two components was observed. Molecular dynamics (MD) simulations predicted the favorable formation of the above complex and identified the dominant interactions between IRB and 2-HP-β-CD. In vitro pharmacological results verified that the inclusion complex not only preserves the binding affinity of IRB for AT1R receptor, but also it slightly increases it. As the complex formulation lacks the problems of the tablet, our approach is a promising new way to improve the efficiency of IRB.

Probing molecular mechanism of inhibitor bindings to bromodomain-containing protein 4 based on molecular dynamics simulations and principal component analysis

It is well known that bromodomain-containing protein 4 (BRD4) has been thought as a promising target utilized for treating various human diseases, such as inflammatory disorders, malignant tumours, acute myelogenous leukaemia (AML), bone diseases, etc. For this study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were integrated together to uncover binding modes of inhibitors 8P9, 8PU, and 8PX to BRD4(1). The results obtained from binding free energy calculations show that van der Waals interactions act as the main regulator in bindings of inhibitors to BRD4(1). The information stemming from PCA reveals that inhibitor associations extremely affect conformational changes, internal dynamics, and movement patterns of BRD4(1). Residue-based free energy decomposition method was wielded to unveil contributions of independent residues to inhibitor bindings and the data signify that hydrogen bonding interactions and hydrophobic interactions are decisive factors affecting bindings of inhibitors to BRD4(1). Meanwhile, eight residues Trp81, Pro82, Val87, Leu92, Leu94, Cys136, Asn140, and Ile146 are recognized as the common hot interaction spots of three inhibitors with BRD4(1). The results from this work are expected to provide a meaningfully theoretical guidance for design and development of effective inhibitors inhibiting of the activity of BRD4.

Structure-based virtual screening of dipeptidyl peptidase 4 inhibitors and their in vitro analysis

Type 2 diabetes mellitus (T2DM) is one of the most widely prevalent metabolic disorders with no cure to date thus remains the most challenging task in the current drug discovery. Therefore, the only strategy to control diabetes prevalence is to develop novel efficacious therapeutics. Dipeptidyl Peptidase 4 (DPP-4) inhibitors are currently used as anti-diabetic drugs for the inhibition of incretins. This study aims to construct the chemical feature based on pharmacophore models for dipeptidyl peptidase IV. The structure-based pharmacophore modeling has been employed to evaluate new inhibitors of DPP-4. A four-featured pharmacophore model was developed from crystal structure of DPP-4 enzyme with 4-(2-aminoethyl) benzenesulfonyl fluoride in its active site via pharmacophore constructing tool of Molecular Operating Environment (MOE) consisting F1 Hyd (hydrophobic region), F2 Hyd|Cat|Don (hydrophobic cationic and donor region), F3 Acc (acceptor region) and F4 Hyd (hydrophobic region). The generated pharmacophore model was used for virtual screening of in-house compound library (the available compounds which were used for initial screening to get the few compounds for the current studies). The resultant selected compounds, after virtual screening were further validated using in vitro assay. Furthermore, structure-activity relationship was carried out for the compounds possessing significant inhibition potential after docking studies. The binding free energy of analogs was evaluated via molecular mechanics generalized Born surface area (MM-GBSA) and Poisson-Boltzmann surface area (MM-PBSA) methods using AMBER 16 as a molecular dynamics (MD) simulation package. Based on potential findings, we report that selected candidates are more likely to be used as DPP-4 inhibitors or as starting leads for the development of novel and potent DPP-4 inhibitors.

Effects of Disulfide Bonds on Binding of Inhibitors to β-Amyloid Cleaving Enzyme 1 Decoded by Multiple Replica Accelerated Molecular Dynamics Simulations

The β-amyloid cleaving enzyme 1 (BACE1) has been thought to be an efficient target for treatment of Alzheimer's disease (AD). Deep insight into inhibitor-BACE1 binding mechanism is of significance for design of potent drugs toward BACE1. In this work, multiple replica accelerated molecular dynamics (MR-aMD) simulations, principal component (PC) analysis, and free energy landscapes were integrated to decode the effect of disulfide bonds (SSBs) in BACE1 on bindings of three inhibitors 3KO, 3KT, and 779 to BACE1. The results from cross-correlation analysis suggest that the breaking of SSBs exerts significant influence on structural flexibility and internal dynamics of inhibitor-bound BACE1. PC analysis and free energy landscapes reveal that the breaking of SSBs not only evidently induces the conformational rearrangement of BACE1 but also highly changes binding poses of three inhibitors in BACE1 and leads to more disordered binding of three inhibitors to BACE1, which is further supported by the increase in binding entropy of inhibitors to BACE1 due to the breaking of SSBs. Residue-based free energy decomposition method was utilized to evaluate contributions of separate residues to inhibitor-BACE1 binding. The results suggest that although the breaking of SSBs in BACE1 does not destroy the interaction network of inhibitors with BACE1 it changes interaction strength of some residues with inhibitors. Meanwhile, the information from residue-based free energy decomposition indicates that residues L91, S96, V130, Y132, Q134, W137, F169, I171, and I179 can be used as efficient targets of drug design toward BACE1.

Gibbs Free Energy Calculation of Mutation in PncA and RpsA Associated With Pyrazinamide Resistance

A central approach for better understanding the forces involved in maintaining protein structures is to investigate the protein folding and thermodynamic properties. The effect of the folding process is often disturbed in mutated states. To explore the dynamic properties behind mutations, molecular dynamic (MD) simulations have been widely performed, especially in unveiling the mechanism of drug failure behind mutation. When comparing wild type (WT) and mutants (MTs), the structural changes along with solvation free energy (SFE), and Gibbs free energy (GFE) are calculated after the MD simulation, to measure the effect of mutations on protein structure. Pyrazinamide (PZA) is one of the first-line drugs, effective against latent Mycobacterium tuberculosis isolates, affecting the global TB control program 2030. Resistance to this drug emerges due to mutations in pncA and rpsA genes, encoding pyrazinamidase (PZase) and ribosomal protein S1 (RpsA) respectively. The question of how the GFE may be a measure of PZase and RpsA stabilities, has been addressed in the current review. The GFE and SFE of MTs have been compared with WT, which were already found to be PZA-resistant. WT structures attained a more stable state in comparison with MTs. The physiological effect of a mutation in PZase and RpsA may be due to the difference in energies. This difference between WT and MTs, depicted through GFE plots, might be useful in predicting the stability and PZA-resistance behind mutation. This study provides useful information for better management of drug resistance, to control the global TB problem.

Characterizing Protein-Ligand Binding Using Atomistic Simulation and Machine Learning: Application to Drug Resistance in HIV-1 Protease

Over the past several decades, atomistic simulations of biomolecules, whether carried out using molecular dynamics or Monte Carlo techniques, have provided detailed insights into their function. Comparing the results of such simulations for a few closely related systems has guided our understanding of the mechanisms by which changes such as ligand binding or mutation can alter the function. The general problem of detecting and interpreting such mechanisms from simulations of many related systems, however, remains a challenge. This problem is addressed here by applying supervised and unsupervised machine learning techniques to a variety of thermodynamic observables extracted from molecular dynamics simulations of different systems. As an important test case, these methods are applied to understand the evasion by human immunodeficiency virus type-1 (HIV-1) protease of darunavir, a potent inhibitor to which resistance can develop via the simultaneous mutation of multiple amino acids. Complex mutational patterns have been observed among resistant strains, presenting a challenge to developing a mechanistic picture of resistance in the protease. In order to dissect these patterns and gain mechanistic insight into the role of specific mutations, molecular dynamics simulations were carried out on a collection of HIV-1 protease variants, chosen to include highly resistant strains and susceptible controls, in complex with darunavir. Using a machine learning approach that takes advantage of the hierarchical nature in the relationships among the sequence, structure, and function, an integrative analysis of these trajectories reveals key details of the resistance mechanism, including changes in the protein structure, hydrogen bonding, and protein-ligand contacts.

Effect of mutations on binding of ligands to guanine riboswitch probed by free energy perturbation and molecular dynamics simulations

Riboswitches can regulate gene expression by direct and specific interactions with ligands and have recently attracted interest as potential drug targets for antibacterial. In this work, molecular dynamics (MD) simulations, free energy perturbation (FEP) and molecular mechanics generalized Born surface area (MM-GBSA) methods were integrated to probe the effect of mutations on the binding of ligands to guanine riboswitch (GR). The results not only show that binding free energies predicted by FEP and MM-GBSA obtain an excellent correlation, but also indicate that mutations involved in the current study can strengthen the binding affinity of ligands GR. Residue-based free energy decomposition was applied to compute ligand-nucleotide interactions and the results suggest that mutations highly affect interactions of ligands with key nucleotides U22, U51 and C74. Dynamics analyses based on MD trajectories indicate that mutations not only regulate the structural flexibility but also change the internal motion modes of GR, especially for the structures J12, J23 and J31, which implies that the aptamer domain activity of GR is extremely plastic and thus readily tunable by nucleotide mutations. This study is expected to provide useful molecular basis and dynamics information for the understanding of the function of GR and possibility as potential drug targets for antibacterial.

Molecular mechanism with regard to the binding selectivity of inhibitors toward FABP5 and FABP7 explored by multiple short molecular dynamics simulations and free energy analyses

Recently, fatty acid binding proteins 5 and 7 (FABP5 and FABP7) have been regarded as the prospective targets for clinically treating multiple diseases related to FABPs.

Effective estimation of the inhibitor affinity of HIV-1 protease via a modified LIE approach

The inhibition of the Human Immunodeficiency Virus Type 1 Protease (HIV-1 PR) can prevent the synthesis of new viruses. Computer-aided drug design (CADD) would enhance the discovery of new therapies, through which the estimation of ligand-binding affinity is critical to predict the most efficient inhibitor. A time-consuming binding free energy method would reduce the usefulness of CADD. The modified linear interaction energy (LIE) approach emerges as an appropriate protocol that performs this task. In particular, the polar interaction free energy, which is obtained via numerically resolving the linear Poisson–Boltzmann equation, plays as an important role in driving the binding mechanism of the HIV-1 PR + inhibitor complex. The electrostatic interaction energy contributes to the attraction between two molecules, but the vdW interaction acts as a repulsive factor between the ligand and the HIV-1 PR. Moreover, the ligands were found to adopt a very strong hydrophobic interaction with the HIV-1 PR. Furthermore, the results obtained corroborate the high accuracy and precision of computational studies with a large correlation coefficient value R = 0.83 and a small RMSE δ_RMSE = 1.25 kcal mol⁻¹. This method is less time-consuming than the other end-point methods, such as the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) and free energy perturbation (FEP) approaches. Overall, the modified LIE approach would provide ligand-binding affinity with HIV-1 PR accurately, precisely, and rapidly, resulting in a more efficient design of new inhibitors.

The inhibition of the Human Immunodeficiency Virus Type 1 Protease (HIV-1 PR) can prevent the synthesis of new viruses.