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

Selectivity Mechanism of Pyrrolopyridone Analogues Targeting Bromodomain 2 of Bromodomain-Containing Protein 4 from Molecular Dynamics Simulations

Bromodomain and extra-terminal domain (BET) proteins play an important role in epigenetic regulation and are linked to several diseases; therefore, they are interesting targets. BET has two bromodomains: bromodomain 1 (BD1) and BD2. Selective targeting of BD1 or BD2 may produce different activities and greater effects than pan-BD inhibitors. However, the selective mechanism of the specific core must be studied at the atomic level. This study determined the effectiveness of pyrrolopyridone analogues to selectively inhibit BD2 using a pan-BD inhibitor (ABBV-075) and a selective-BD2 inhibitor (ABBV-744). Molecular dynamics simulations and calculations of binding free energies were used to systematically study the selectivity of BD2 inhibition by the pyrrolopyridone analogues. Overall, the pyrrolopyridone analogue inhibitors targeting BD2 interacted mainly with the following amino acid pairs between bromodomain-containing protein 4 (BRD4)-BD1 and BRD4-BD2 complexes: I146/V439, N140/N433, D144/H437, P82/P375, V87/V380, D88/D381, and Y139/Y432. The pyrrolopyridone analogues targeting BRD4-BD2 were divided into five regions based on selectivity mechanism. These results suggest that the R3 and R5 regions of pyrrolopyridone analogues can be modified to improve the selectivity between BRD4-BD1 and BRD4-BD2. The selectivity of BD2 inhibition by pyrrolopyridone analogues can be used to design novel BD2 inhibitors based on a pyrrolopyridone core.

Unveiling the Inhibitory Potentials of Peptidomimetic Azanitriles and Pyridyl Esters towards SARS-CoV-2 Main Protease: A Molecular Modelling Investigation

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is responsible for COVID-19, which was declared a global pandemic in March 2020 by the World Health Organization (WHO). Since SARS-CoV-2 main protease plays an essential role in the virus’s life cycle, the design of small drug molecules with lower molecular weight has been a promising development targeting its inhibition. Herein, we evaluated the novel peptidomimetic azatripeptide and azatetrapeptide nitriles against SARS-CoV-2 main protease. We employed molecular dynamics (MD) simulations to elucidate the selected compounds’ binding free energy profiles against SARS-CoV-2 and further unveil the residues responsible for the drug-binding properties. Compound 8 exhibited the highest binding free energy of −49.37 ± 0.15 kcal/mol, followed by compound 7 (−39.83 ± 0.19 kcal/mol), while compound 17 showed the lowest binding free energy (−23.54 ± 0.19 kcal/mol). In addition, the absorption, distribution, metabolism, and excretion (ADME) assessment was performed and revealed that only compound 17 met the drug-likeness parameters and exhibited high pharmacokinetics to inhibit CYP1A2, CYP2C19, and CYP2C9 with better absorption potential and blood-brain barrier permeability (BBB) index. The additional intermolecular evaluations suggested compound 8 as a promising drug candidate for inhibiting SARS-CoV-2 Mpro. The substitution of isopropane in compound 7 with an aromatic benzene ring in compound 8 significantly enhanced the drug’s ability to bind better at the active site of the SARS-CoV-2 Mpro.

Structural Investigations and Binding Mechanisms of Oseltamivir Drug Resistance Conferred by the E119V Mutation in Influenza H7N9 Virus

The use of vaccinations and antiviral medications have gained popularity in the therapeutic management of avian influenza H7N9 virus lately. Antiviral medicines are more popular due to being readily available. The presence of the neuraminidase protein in the avian influenza H7N9 virus and its critical role in the cleavage of sialic acid have made it a target drug in the development of influenza virus drugs. Generally, the neuraminidase proteins have common conserved amino acid residues and any mutation that occurs around or within these conserved residues affects the susceptibility and replicability of the influenza H7N9 virus. Herein, we investigated the interatomic and intermolecular dynamic impacts of the experimentally reported E119V mutation on the oseltamivir resistance of the influenza H7N9 virus. We extensively employed molecular dynamic (MD) simulations and subsequent post-MD analyses to investigate the binding mechanisms of oseltamivir-neuraminidase wildtype and E119V mutant complexes. The results revealed that the oseltamivir-wildtype complex was more thermodynamically stable than the oseltamivir-E119V mutant complex. Oseltamivir exhibited a greater binding affinity for wildtype (−15.46 ± 0.23 kcal/mol) relative to the E119V mutant (−11.72 ± 0.21 kcal/mol). The decrease in binding affinity (−3.74 kcal/mol) was consistent with RMSD, RMSF, SASA, PCA, and hydrogen bonding profiles, confirming that the E119V mutation conferred lower conformational stability and weaker protein–ligand interactions. The findings of this oseltamivir-E119V mutation may further assist in the design of compounds to overcome E119V mutation in the treatment of influenza H7N9 virus patients.

Computational studies with flavonoids and terpenoids as BRPF1 inhibitors: in silico biological activity prediction, molecular docking, molecular dynamics simulations, MM/PBSA calculations

The BRPF1 protein is encoded by the BRPF1 gene. In addition, the BRPF1 gene is known to be upregulated in leukaemia. Recent studies have shown that it is also overexpressed in hepatocellular carcinoma (HCC) as well. Therefore, BRPF1 is a significant target for anti-cancer drug development studies, especially on HCC. 40 terpenoids and flavonoids were chosen because of their anticancer properties given in the literature. In this study, the biological activity of molecules was also investigated with in silico structure-activity relationship analysis. In addition, interactions between a series of terpenoids and flavonoids and the BRPF1 protein were investigated by molecular docking and molecular dynamics simulations. The energy change caused by the interactions of BRPF1 with different compounds was also evaluated by MM/PBSA calculations. It has been revealed that compound 5 (-9.2 kcal/mol), a kind of secoclerodane type diterpenoid, has a higher affinity both compared to other flavonoids and terpenoids, and 9F9 (-7.9 kcal/mol), a selective BRPF1 inhibitor. The study presented in this article demonstrates that compound 5, as a natural product, could form a chemical scaffold for the development of selective BRPF1 bromodomain inhibitors.

Impact of the R292K Mutation on Influenza A (H7N9) Virus Resistance towards Peramivir: A Molecular Dynamics Perspective

In March 2013, a novel avian influenza A (H7N9) virus emerged in China. By March 2021, it had infected more than 1500 people, raising concerns regarding its epidemic potential. Similar to the highly pathogenic H5N1 virus, the H7N9 virus causes severe pneumonia and acute respiratory distress syndrome in most patients. Moreover, genetic analysis showed that this avian H7N9 virus carries human adaptation markers in the hemagglutinin and polymerase basic 2 (PB2) genes associated with cross-species transmissibility. Clinical studies showed that a single mutation, neuraminidase (NA) R292K (N2 numbering), induces resistance to peramivir in the highly pathogenic H7N9 influenza A viruses. Therefore, to evaluate the risk for human public health and understand the possible source of drug resistance, we assessed the impact of the NA-R292K mutation on avian H7N9 virus resistance towards peramivir using various molecular dynamics approaches. We observed that the single point mutation led to a distorted peramivir orientation in the enzyme active site which, in turn, perturbed the inhibitor's binding. The R292K mutation induced a decrease in the interaction among neighboring amino acid residues when compared to its wild-type counterpart, as shown by the high degree of fluctuations in the radius of gyration. MM/GBSA calculations revealed that the mutation caused a decrease in the drug binding affinity by 17.28 kcal/mol when compared to the that for the wild-type enzyme. The mutation caused a distortion of hydrogen bond-mediated interactions with peramivir and increased the accessibility of water molecules around the K292 mutated residue.

Molecular Dynamic Simulations of Bromodomain and Extra-Terminal Protein 4 Bonded to Potent Inhibitors

Bromodomain and extra-terminal domain (BET) subfamily is the most studied subfamily of bromodomain-containing proteins (BCPs) family which can modulate acetylation signal transduction and produce diverse physiological functions. Thus, the BET family can be treated as an alternative strategy for targeting androgen-receptor (AR)-driven cancers. In order to explore the effect of inhibitors binding to BRD4 (the most studied member of BET family), four 150 ns molecular dynamic simulations were performed (free BRD4, Cpd4-BRD4, Cpd9-BRD4 and Cpd19-BRD4). Docking studies showed that Cpd9 and Cpd19 were located at the active pocket, as well as Cpd4. Molecular dynamics (MD) simulations indicated that only Cpd19 binding to BRD4 can induce residue Trp81-Ala89 partly become α-helix during MD simulations. MM-GBSA calculations suggested that Cpd19 had the best binding effect with BRD4 followed by Cpd4 and Cpd9. Computational alanine scanning results indicated that mutations in Phe83 made the greatest effects in Cpd9-BRD4 and Cpd19-BRD4 complexes, showing that Phe83 may play crucial roles in Cpd9 and Cpd19 binding to BRD4. Our results can provide some useful clues for further BCPs family search.

Theoretically exploring selective-binding mechanisms of BRD4 through integrative computational approaches

The origin of cancer is related to the dysregulation of multiple signal pathways and of physiological processes. Bromodomain-containing protein 4 (BRD4) has become an attractive target for the development of anticancer and anti-inflammatory agents since it can epigenetically regulate the transcription of growth-promoting genes. The synthesized BRD4 inhibitors with new chemical structures can reduce the drug resistance, but their binding modes and the inhibitory mechanism remain unclear. Here, we initially constructed robust QSAR models based on 68 reported tetrahydropteridin analogues using topomer CoMFA and HQSAR. On the basis of QSAR results, we designed 16 novel tetrahydropteridin analogues with modified structures and carried out docking studies. Instead of significant hydrogen bondings with amino acid residue Asn140 as reported in previous research, the molecular docking modelling suggested a novel docking pose that involves the amino acid residues (Trp81, Pro82, Val87, Leu92, Leu94, Cys136, Asp144, and Ile146) at the active site of BRD4. The MD simulations, free energy calculations, and residual energy contributions all indicate that hydrophobic interactions are decisive factors affecting bindings between inhibitors and BRD4. The current study provides new insights that can aid the discovery of BRD4 inhibitors with enhanced anti-cancer ability.

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 dynamics study of CDC25BR492L mutant causing the activity decrease of CDC25B

Cell division cycle 25B (CDC25B) was responsible for regulating the various stages of cell division in the cell cycle. R492L was one of the common types of CDC25B mutants. Researches showed that compared to CDC25B^WT, CDC25B^R492L mutant had a ∼100-fold reduction in the rate constant for forming phosphatase intermediate (k₂). However, the molecular basis of how the CDC25B^R492L mutant influenced the process of binding between CDC25B and CDK2/CyclinA was not yet known. Therefore, the optimizations of three-dimensional structure of the CDC25B^WT-CDK2/CyclinA system and the CDC25B^R492L-CDK2/CyclinA system were constructed by ZDOCK and RDOCK, and five methods were employed to verify the reasonability of the docking structure. Then the molecular dynamics simulations on the two systems were performed to explore the reason why CDC25B^R492L mutant caused the weak interactions between CDC25B^R492L and CDK2/CyclinA, respectively. The remote docking site (Arg488-Tyr497) and the second active site (Lys538-Arg544) of CDC25B were observed to have high fluctuations in the CDC25B^R492L-CDK2/CyclinA system with post-analysis, where the high fluctuation of these two regions resulted in weak interactions between CD25B and CDK2. In addition, Asp38-Glu42 and Asp206-Asp210 of CDK2 showed the slightly descending fluctuation, and CDK2 revealed an enhanced the self-interaction, which made CDK2 keep a relatively stable state in the CDC25B^R492L-CDK2/CyclinA system. Finally, Leu492 of CDC25B was speculated to be the key residue, which had great effects on the binding between CDC25B^R492L and CDK2 in the CDC25B^R492L-CDK2/CyclinA system. Consequently, overall analyses appeared in this study ultimately offered a helpful understanding of the weak interactions between CDC25B^R492L and CDK2.

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.

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.

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.

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.

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.