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

A novel Imidazo[1,2-a]pyridine derivative modulates active KRASG12D through off-like conformational shifts in switch-I and switch-II regions, mimicking inactive KRASG12D

KRASG12D are the most prevalent oncogenic mutations and a promising target for solid tumor therapies. However, its inhibition exhibits tremendous challenge due to the necessity of high binding affinity to obviate the need for covalent binders. Here we report the evidence of a novel class of Imidazo[1,2-a]pyridine derivative as potentially significant novel inhibitors of KRASG12D, discovered through extensive ligand-based screening against 2-[(2R)-piperidin-2-yl]-1H-indole, an important scaffold for KRASG12D inhibition via switch-I/II (S-I/II) pocket. The proposed compounds exhibited similar binding affinities and overlapped pose configurations to 2-[(2R)-piperidin-2-yl]-1H-indole, serving as a reliable starting point for drug discovery. Comparative free energy profiles demonstrated that C4 [2-methyl-3-((5-phenyl-1H-1,2,4-triazol-3-yl)methyl)imidazo[1,2-a]pyridine] effectively shifted the protein to a stable low-energy conformation via a prominent transition state. The conformational changes across the transition revealed the conformational shift of switch-I and II to a previously known off-like conformation of inactive KRASG12D with rmsd of 0.91 Å. These conformations were even more prominent than the privileged scaffold 2-[(2R)-piperidin-2-yl]-1H-indole. The representative structure overlay of C4 and another X-ray crystallography solved BI-2852 bound inactive KRASG12D revealed that Switch-I and II exhibited off-like conformations. The cumulative variance across the first eigenvalue that accounted for 57 % of the collective variance validated this on-to-off transition. In addition, the relative interaction of C4 binding showed consistent patterns with BI-2852. Taken together, our results support the inhibitory activity of [2-methyl-3-((5-phenyl-1H-1,2,4-triazol-3-yl)methyl)imidazo[1,2-a]pyridine] by shifting active KRASG12D to an inactive conformation.

Binding selectivity analysis of new delhi metallo-beta-lactamase-1 inhibitors using molecular dynamics simulations: Exploring possibilities for decoding antimicrobial drug resistance

BACKGROUND

New Delhi metallo-beta-lactamase-1 (NDM1) confers resistance to several bacterial species against a broad range of beta-lactam antibiotics and turning them into superbugs that pose a significant threat to healthcare systems worldwide. As such, it is a potentially relevant biological target for counteracting bacterial infections. Given the lack of effective treatment options against NDM1 producing bacteria, finding a reliable inhibitor for the NDM1 enzyme is crucial.

METHODS

Using molecular dynamics simulations, the binding selectivities and affinities of three ligands, viz. PNK, 3S0, and N1G were investigated against NDM1.

RESULTS

The results indicate that N1G binds with more affinity to NDM1 than PNK and 3S0. The binding energy decomposition analysis revealed that residues I35, W93, H189, K211, and N220 showed significant binding energies with PNK, 3S0, and N1G, and hence are crucially involved in the binding of the ligands to NDM1. Molecular dynamics trajectory analysis further elicited that the ligands influence dynamic flexibility of NDM1 morphology, which contributes to the partial selectivities of PNK, 3S0, and N1G.

CONCLUSIONS

This in silico study offers a vital information for developing potential NDM1 inhibitors with high selectivity. Nevertheless, in vitro and in vivo experimental validation is mandated to extend the possible applications of these ligands as NDM1 inhibitors that succor in combating antimicrobial resistance.

The switch states of the GDP-bound HRAS affected by point mutations: a study from Gaussian accelerated molecular dynamics simulations and free energy landscapes

Point mutations play a vital role in the conformational transformation of HRAS. In this work, Gaussian accelerated molecular dynamics (GaMD) simulations followed by constructions of free energy landscapes (FELs) were adopted to explore the effect of mutations D33K, A59T and L120A on conformation states of the GDP-bound HRAS. The results from the post-processing analyses on GaMD trajectories suggest that mutations alter the flexibility and motion modes of the switch domains from HRAS. The analyses from FELs show that mutations induce more disordered states of the switch domains and affect interactions of GDP with HRAS, implying that mutations yield a vital effect on the binding of HRAS to effectors. The GDP-residue interaction network revealed by our current work indicates that salt bridges and hydrogen bonding interactions (HBIs) play key roles in the binding of GDP to HRAS. Furthermore, instability in the interactions of magnesium ions and GDP with the switch SI leads to the extreme disorder of the switch domains. This study is expected to provide the energetic basis and molecular mechanism for further understanding the function of HRAS.Communicated by Ramaswamy H. Sarma.

Mutations influence the conformational dynamics of the GDP/KRAS complex

Mutations near allosteric sites can have a significant impact on the function of KRAS. Three specific mutations, K104Q, G12D/K104Q, and G12D/G75A, which are located near allosteric positions, were selected to investigate the molecular mechanisms behind mutation-induced influences on the activity of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations followed by the principal component analysis (PCA) were performed to improve the sampling of conformational states. The results revealed that these mutations significantly alter the structural flexibility, correlated motions, and dynamic behavior of the switch regions that are essential for KRAS binding to effectors or regulators. Furthermore, the mutations have a significant impact on the hydrogen bonding interactions between GDP and the switch regions, as well as on the electrostatic interactions of magnesium ions (Mg2+) with these regions. Our results verified that these mutations strongly influence the binding of KRAS to its effectors or regulators and allosterically regulate the activity. We believe that this work can provide valuable theoretical insights into a deeper understanding of KRAS function.Communicated by Ramaswamy H. Sarma.

Unveiling the State Transition Mechanisms of Ras Proteins through Enhanced Sampling and QM/MM Simulations

In cells, wild-type RasGTP complexes exist in two distinct states: active State 2 and inactive State 1. These complexes regulate their functions by transitioning between the two states. However, the mechanisms underlying this state transition have not been clearly elucidated. To address this, we conducted a detailed simulation study to characterize the energetics of the stable states involved in the state transitions of the HRasGTP complex, specifically from State 2 to State 1. This was achieved by employing multiscale quantum mechanics/molecular mechanics and enhanced sampling molecular dynamics methods. Based on the simulation results, we constructed the two-dimensional free energy landscapes that provide crucial information about the conformational changes of the HRasGTP complex from State 2 to State 1. Furthermore, we also explored the conformational changes from the intermediate state to the product state during guanosine triphosphate hydrolysis. This study on the conformational changes involved in the HRas state transitions serves as a valuable reference for understanding the corresponding events of both KRas and NRas as well.

Identification of novel scaffolds targeting SIRT3 through molecular modeling techniques for the treatment of Hepatocellular carcinoma

Hepatocellular carcinoma is one of the top causes of cancer-related death globally. SIRT3 belongs to the Sirtuin family of proteins, a collection of NAD+-dependent enzymes that play a role in controlling several cellular functions, including metabolism, aging, and stress response. SIRT3 expression has been discovered to be often downregulated in HCC tissues relative to normal liver tissues. Hence, SIRT3 may function as a tumor suppressor in HCC. In the present study, pharmacophore-based virtual screening of a small molecule database was performed initially, and then the screened hits were docked to the active site of SIRT3 to choose the best binding modes. One co-crystal ligand (PDB name: 1NQ) was utilized as a template to generate pharmacophore model query. A total of 0.2 million compounds from the VITAS-M Lab database were downloaded and prepared for virtual screening. Following database preparation, ligand-based virtual screening was performed using the pharmacophore query model generated in the previous phase. The compounds with the same pharmacophoric characteristics as the query at the same distance were screened. There were a total of 74 hits that matched the query model. These compounds were then docked to the SIRT3 using the standard precision protocol of the glide tool. To select hits with high binding affinities, a threshold of -8 kcal/mol was used. Based on the glide gscore, two hits were chosen. These two hits were selected to investigate the stability of the protein-ligand complex by molecular dynamics simulation. All of these findings indicate that the selected hit compounds C1 and C2 can serve as lead compounds in inhibiting the biological activity of SIRT3 requiring further detailed investigations.Communicated by Ramaswamy H. Sarma.

Impacts of Mutations in the P-Loop on Conformational Alterations of KRAS Investigated with Gaussian Accelerated Molecular Dynamics Simulations

G12 mutations heavily affect conformational transformation and activity of KRAS. In this study, Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the GDP-bound wild-type (WT), G12A, G12D, and G12R KRAS to probe mutation-mediated impacts on conformational alterations of KRAS. The results indicate that three G12 mutations obviously affect the structural flexibility and internal dynamics of the switch domains. The analyses of the free energy landscapes (FELs) suggest that three G12 mutations induce more conformational states of KRAS and lead to more disordered switch domains. The principal component analysis shows that three G12 mutations change concerted motions and dynamics behavior of the switch domains. The switch domains mostly overlap with the binding region of KRAS to its effectors. Thus, the high disorder states and concerted motion changes of the switch domains induced by G12 mutations affect the activity of KRAS. The analysis of interaction network of GDP with KRAS signifies that the instability in the interactions of GDP and magnesium ion with the switch domain SW1 drives the high disordered state of the switch domains. This work is expected to provide theoretical aids for understanding the function of KRAS.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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