1
|
Wang J, Yang W, Zhao L, Wei B, Chen J. Binding Mechanism of Inhibitors to BRD4 and BRD9 Decoded by Multiple Independent Molecular Dynamics Simulations and Deep Learning. Molecules 2024; 29:1857. [PMID: 38675678 PMCID: PMC11054187 DOI: 10.3390/molecules29081857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2024] [Revised: 04/16/2024] [Accepted: 04/16/2024] [Indexed: 04/28/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Jian Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Wanchun Yang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Lu Zhao
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Benzheng Wei
- Center for Medical Artificial Intelligence, Shandong University of Traditional Chinese Medicine, Qingdao 266112, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| |
Collapse
|
2
|
Bao H, Wang W, Sun H, Chen J. The switch states of the GDP-bound HRAS affected by point mutations: a study from Gaussian accelerated molecular dynamics simulations and free energy landscapes. J Biomol Struct Dyn 2024; 42:3363-3381. [PMID: 37216340 DOI: 10.1080/07391102.2023.2213355] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Accepted: 05/04/2023] [Indexed: 05/24/2023]
Abstract
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.
Collapse
Affiliation(s)
- Huayin Bao
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
3
|
Shen C, Yin J, Wang M, Yu Z, Xu X, Zhou Z, Hu Y, Xia C, Hu G. Mutations influence the conformational dynamics of the GDP/KRAS complex. J Biomol Struct Dyn 2024:1-14. [PMID: 38529923 DOI: 10.1080/07391102.2024.2331627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2024] [Accepted: 02/20/2024] [Indexed: 03/27/2024]
Abstract
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.
Collapse
Affiliation(s)
- Congcong Shen
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou, China
| | - Jie Yin
- Qingyun People's Hospital, Dezhou, China
| | - Min Wang
- Qingyun People's Hospital, Dezhou, China
| | - Zhiping Yu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou, China
| | - Xin Xu
- School of Science, Xi'an Polytechnic University, Xi'an, China
| | - Zhongshun Zhou
- School of Science, Xi'an Polytechnic University, Xi'an, China
| | - Yingshi Hu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou, China
| | - Caijuan Xia
- School of Science, Xi'an Polytechnic University, Xi'an, China
| | - Guodong Hu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou, China
| |
Collapse
|
4
|
Basu D, Dastidar SG. Molecular Dynamics and Machine Learning reveal distinguishing mechanisms of Competitive Ligands to perturb α,β-Tubulin. Comput Biol Chem 2024; 108:108004. [PMID: 38157659 DOI: 10.1016/j.compbiolchem.2023.108004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2023] [Revised: 11/25/2023] [Accepted: 12/08/2023] [Indexed: 01/03/2024]
Abstract
The mechanisms of action of ligands competing for the Colchicine Binding Site (CBS) of the α,β-Tubulin are non-standard compared to the commonly witnessed ligand-induced inhibition of proteins. This is because their potencies are not solely judged by the binding affinity itself, but also by their capacity to bias the conformational states of the dimer. Regarding the latter requirement, it is observed that ligands competing for the same pocket that binds colchicine exhibit divergence in potential clinical outcomes. Molecular dynamics-based ∼5.2 µs sampling of α,β-Tubulin complexed with four different ligands has revealed that each ligand has its customized way of influencing the complex. Primarily, it is the proportion of twisting and/or bending characteristic of modes of the intrinsic dynamics which is revealed to be 'fundamental' to tune this variation in the mechanism. The milder influence of 'bending' makes a ligand (TUB092), better classifiable under the group of vascular disrupting agents (VDAs), which are phenotypically effective on cytoskeletons; whereas a stronger impact of 'bending' makes the classical ligand Colchicine (COL) a better Anti-Mitotic agent (AMA). Two other ligands BAL27862 (2RR) and Nocodazole (NZO) fall in the intermediate zone as they fail to explicitly induce bending modes. Random Forest Classification method and K-means Clustering is applied to reveal the efficiency of Machine Learning methods in classifying the Tubulin conformations according to their ligand-specific perturbations and to highlight the significance of specific amino acid residues, mostly positioned in the α-β and β-β interfaces involved in the mechanism. These key residues responsible to yield discriminative actions of the ligands are likely to be highly useful in future endeavours to design more precise inhibitors.
Collapse
Affiliation(s)
- Debadrita Basu
- Biological Sciences, Bose Institute, EN 80, Sector V, Bidhan Nagar, Kolkata 700091, India
| | - Shubhra Ghosh Dastidar
- Biological Sciences, Bose Institute, EN 80, Sector V, Bidhan Nagar, Kolkata 700091, India.
| |
Collapse
|
5
|
Yu Z, Wang Z, Cui X, Cao Z, Zhang W, Sun K, Hu G. Conformational States of the GDP- and GTP-Bound HRAS Affected by A59E and K117R: An Exploration from Gaussian Accelerated Molecular Dynamics. Molecules 2024; 29:645. [PMID: 38338389 PMCID: PMC10856033 DOI: 10.3390/molecules29030645] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2023] [Revised: 01/01/2024] [Accepted: 01/25/2024] [Indexed: 02/12/2024] Open
Abstract
The HRAS protein is considered a critical target for drug development in cancers. It is vital for effective drug development to understand the effects of mutations on the binding of GTP and GDP to HRAS. We conducted Gaussian accelerated molecular dynamics (GaMD) simulations and free energy landscape (FEL) calculations to investigate the impacts of two mutations (A59E and K117R) on GTP and GDP binding and the conformational states of the switch domain. Our findings demonstrate that these mutations not only modify the flexibility of the switch domains, but also affect the correlated motions of these domains. Furthermore, the mutations significantly disrupt the dynamic behavior of the switch domains, leading to a conformational change in HRAS. Additionally, these mutations significantly impact the switch domain's interactions, including their hydrogen bonding with ligands and electrostatic interactions with magnesium ions. Since the switch domains are crucial for the binding of HRAS to effectors, any alterations in their interactions or conformational states will undoubtedly disrupt the activity of HRAS. This research provides valuable information for the design of drugs targeting HRAS.
Collapse
Affiliation(s)
- Zhiping Yu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China; (Z.Y.); (Z.C.)
| | - Zhen Wang
- Pingyin People’s Hospital, Jinan 250400, China; (Z.W.); (X.C.)
| | - Xiuzhen Cui
- Pingyin People’s Hospital, Jinan 250400, China; (Z.W.); (X.C.)
| | - Zanxia Cao
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China; (Z.Y.); (Z.C.)
| | - Wanyunfei Zhang
- School of Science, Xi’an Polytechnic University, Xi’an 710048, China; (W.Z.); (K.S.)
| | - Kunxiao Sun
- School of Science, Xi’an Polytechnic University, Xi’an 710048, China; (W.Z.); (K.S.)
| | - Guodong Hu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China; (Z.Y.); (Z.C.)
| |
Collapse
|
6
|
Sun J, Liang S, Liu X, Zhang S, Li M, Zhang Q, Chen J. Insights into the selectivity of a brain-penetrant CDK4/6 vs CDK1/2 inhibitor for glioblastoma used in multiple replica molecular dynamics simulations. J Biomol Struct Dyn 2023:1-20. [PMID: 38112295 DOI: 10.1080/07391102.2023.2294175] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Accepted: 11/23/2023] [Indexed: 12/21/2023]
Abstract
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.
Collapse
Affiliation(s)
- Jiahao Sun
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Shanshan Liang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Meng Li
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
7
|
Bhadhadhara K, Jani V, Koulgi S, Sonavane U, Joshi R. Studying early structural changes in SOS1 mediated KRAS activation mechanism. Curr Res Struct Biol 2023; 7:100115. [PMID: 38188543 PMCID: PMC10765296 DOI: 10.1016/j.crstbi.2023.100115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Revised: 11/28/2023] [Accepted: 11/29/2023] [Indexed: 01/09/2024] Open
Abstract
KRAS activation is known to be modulated by a guanine nucleotide exchange factor (GEF), namely, Son of Sevenless1 (SOS1). SOS1 facilitates the exchange of GDP to GTP thereby leading to activation of KRAS. The binding of GDP/GTP to KRAS at the REM/allosteric site of SOS1 regulates the activation of KRAS at CDC25/catalytic site by facilitating its exchange. Different aspects of the allosteric activation of KRAS through SOS1 are still being explored. To understand the SOS1 mediated activation of KRAS, molecular dynamics simulations for a total of nine SOS1 complexes (KRAS-SOS1-KRAS) were performed. These nine systems comprised different combinations of KRAS-bound nucleotides (GTP/GDP) at REM and CDC25 sites of SOS1. Various conformational and thermodynamic parameters were analyzed for these simulation systems. MMPBSA free energy analysis revealed that binding at CDC25 site of SOS1 was significantly low for GDP-bound KRAS as compared to that of GTP-bound KRAS. It was observed that presence of either GDP/GTP bound KRAS at the REM site of SOS1 affected the activation related changes in the KRAS present at CDC25 site. The conformational changes at the catalytic site of SOS1 resulting from GDP/GTP-bound KRAS at the allosteric changes may hint at KRAS activation through different pathways (slow/fast/rare). The allosteric effect on activation of KRAS at CDC25 site may be due to conformations adopted by switch-I, switch-II, beta2 regions of KRAS at REM site. The effect of structural rearrangements occurring at allosteric KRAS may have led to increased interactions between SOS1 and KRAS at both the sites. The SOS1 residues involved in these important interactions with KRAS at the REM site were R694, S732 and K735. Whereas the ones interacting with KRAS at CDC25 site were S807, W809 and K814. This may suggest the crucial role of these residues in guiding the allosteric activation of KRAS at CDC25 site. The conformational shifts observed in the switch-I, switch-II and alpha3 regions of KRAS at CDC25 site may be attributed to be a part of allosteric activation. The binding affinities, interacting residues and conformational dynamics may provide an insight into development of inhibitors targeting the SOS1 mediated KRAS activation.
Collapse
Affiliation(s)
- Kirti Bhadhadhara
- High Performance Computing-Medical & Bioinformatics Applications Group, Centre for Development of Advanced Computing (C-DAC), Innovation Park, Panchawati, Pashan, Pune, 411008, India
| | - Vinod Jani
- High Performance Computing-Medical & Bioinformatics Applications Group, Centre for Development of Advanced Computing (C-DAC), Innovation Park, Panchawati, Pashan, Pune, 411008, India
| | - Shruti Koulgi
- High Performance Computing-Medical & Bioinformatics Applications Group, Centre for Development of Advanced Computing (C-DAC), Innovation Park, Panchawati, Pashan, Pune, 411008, India
| | - Uddhavesh Sonavane
- High Performance Computing-Medical & Bioinformatics Applications Group, Centre for Development of Advanced Computing (C-DAC), Innovation Park, Panchawati, Pashan, Pune, 411008, India
| | - Rajendra Joshi
- High Performance Computing-Medical & Bioinformatics Applications Group, Centre for Development of Advanced Computing (C-DAC), Innovation Park, Panchawati, Pashan, Pune, 411008, India
| |
Collapse
|
8
|
Lu H, Hu Z, Faraudo J, Martí J. In silico design of a lipid-like compound targeting KRAS4B-G12D through non-covalent bonds. NANOSCALE 2023; 15:19359-19368. [PMID: 38014474 DOI: 10.1039/d3nr04513g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
One of the most common drivers in human cancer is the peripheral membrane protein KRAS4B, able to promote oncogenic signalling. To signal, oncogenic KRAS4B not only requires a sufficient nucleotide exchange, but also needs to recruit effectors by exposing its effector-binding sites while anchoring to the phospholipid bilayer where KRAS4B-mediated signalling events occur. The enzyme phosphodiesterase-δ plays an important role in sequestering KRAS4B from the cytoplasm and targeting it to cellular membranes of different cell species. In this work, we present an in silico design of a lipid-like compound that has the remarkable feature of being able to target both an oncogenic KRAS4B-G12D mutant and the phosphodiesterase-δ enzyme. This double action is accomplished by adding a lipid tail (analogous to the farnesyl group of the KRAS4B protein) to an previously known active compound (2H-1,2,4-benzothiadiazine, 3,4-dihydro-,1,1-dioxide). The proposed lipid-like molecule was found to lock KRAS4B-G12D in its GDP-bound state by adjusting the effector-binding domain to be blocked by the interface of the lipid bilayer. Meanwhile, it can tune GTP-bound KRAS4B-G12D to shift from the active orientation state to the inactive state. The proposed compound is also observed to stably accommodate itself in the prenyl-binding pocket of phosphodiesterase-δ, which impairs KRAS4B enrichment at the lipid bilayer, potentially reducing the proliferation of KRAS4B inside the cytoplasm and its anchoring at the bilayer. In conclusion, we report a potential inhibitor of KRAS4B-G12D with a lipid tail attached to a specific warhead, a compound which has not yet been considered for drugs targeting RAS mutants. Our work provides new ways to target KRAS4B-G12D and can also foster drug discovery efforts for the targeting of oncogenes of the RAS family and beyond.
Collapse
Affiliation(s)
- Huixia Lu
- Institut de Ciencia de Materials de Barcelona (ICMAB-CSIC), Campus de la UAB, Bellaterra, Barcelona E-08193, Spain.
- Department of Physics, Technical University of Catalonia-Barcelona Tech, B5-209 Northern Campus, Jordi Girona 1-3, 08034 Barcelona, Catalonia, Spain.
| | - Zheyao Hu
- Department of Physics, Technical University of Catalonia-Barcelona Tech, B5-209 Northern Campus, Jordi Girona 1-3, 08034 Barcelona, Catalonia, Spain.
| | - Jordi Faraudo
- Institut de Ciencia de Materials de Barcelona (ICMAB-CSIC), Campus de la UAB, Bellaterra, Barcelona E-08193, Spain.
| | - Jordi Martí
- Department of Physics, Technical University of Catalonia-Barcelona Tech, B5-209 Northern Campus, Jordi Girona 1-3, 08034 Barcelona, Catalonia, Spain.
| |
Collapse
|
9
|
Sun J, Liu X, Zhang S, Li M, Zhang Q, Chen J. Molecular insights and optimization strategies for the competitive binding of engineered ACE2 proteins: a multiple replica molecular dynamics study. Phys Chem Chem Phys 2023; 25:28479-28496. [PMID: 37846774 DOI: 10.1039/d3cp03392a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2023]
Abstract
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.
Collapse
Affiliation(s)
- Jiahao Sun
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Meng Li
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250357, China.
| |
Collapse
|
10
|
Ratnasinghe BD, Haque N, Wagenknecht JB, Jensen DR, Valdivia Esparza GK, Leverence EN, Milech De Assuncao T, Mathison AJ, Lomberk G, Smith BC, Volkman BF, Urrutia R, Zimmermann MT. Beyond structural bioinformatics for genomics with dynamics characterization of an expanded KRAS mutational landscape. Comput Struct Biotechnol J 2023; 21:4790-4803. [PMID: 37841325 PMCID: PMC10570560 DOI: 10.1016/j.csbj.2023.10.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2023] [Revised: 10/04/2023] [Accepted: 10/04/2023] [Indexed: 10/17/2023] Open
Abstract
Current capabilities in genomic sequencing outpace functional interpretations. Our previous work showed that 3D protein structure calculations enhance mechanistic understanding of genetic variation in sequenced tumors and patients with rare diseases. The KRAS GTPase is among the critical genetic factors driving cancer and germline conditions. Because KRAS-altered tumors frequently harbor one of three classic hotspot mutations, nearly all studies have focused on these mutations, leaving significant functional ambiguity across the broader KRAS genomic landscape observed in cancer and non-cancer diseases. Herein, we extend structural bioinformatics with molecular simulations to study an expanded landscape of 86 KRAS mutations. We identify multiple coordinated changes strongly associated with experimentally established KRAS biophysical and biochemical properties. The patterns we observe span hotspot and non-hotspot alterations, which can all dysregulate Switch regions, producing mutation-restricted conformations with different effector binding propensities. We experimentally measured mutation thermostability and identified shared and distinct patterns with simulations. Our results indicate mutation-specific conformations, which show potential for future research into how these alterations reverberate into different molecular and cellular functions. The data we present is not predictable using current genomic tools, demonstrating the added functional information derived from molecular simulations for interpreting human genetic variation.
Collapse
Affiliation(s)
- Brian D. Ratnasinghe
- Bioinformatics Research and Development Laboratory, Linda T. and John A. Mellowes Center for Genomic Sciences and Precision Medicine, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | - Neshatul Haque
- Bioinformatics Research and Development Laboratory, Linda T. and John A. Mellowes Center for Genomic Sciences and Precision Medicine, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | - Jessica B. Wagenknecht
- Bioinformatics Research and Development Laboratory, Linda T. and John A. Mellowes Center for Genomic Sciences and Precision Medicine, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | - Davin R. Jensen
- Department of Biochemistry, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | - Guadalupe K. Valdivia Esparza
- Structural Genomics Unit, Linda T. and John A. Mellowes Center for Genomic Sciences and Precision Medicine, Medical College of Wisconsin, Milwaukee, WI 53226, USA
- Department of Biochemistry, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | - Elise N. Leverence
- Structural Genomics Unit, Linda T. and John A. Mellowes Center for Genomic Sciences and Precision Medicine, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | - Thiago Milech De Assuncao
- Structural Genomics Unit, Linda T. and John A. Mellowes Center for Genomic Sciences and Precision Medicine, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | - Angela J. Mathison
- Department of Surgery, Medical College of Wisconsin, Milwaukee, WI 53226, USA
- Structural Genomics Unit, Linda T. and John A. Mellowes Center for Genomic Sciences and Precision Medicine, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | - Gwen Lomberk
- Department of Surgery, Medical College of Wisconsin, Milwaukee, WI 53226, USA
- Structural Genomics Unit, Linda T. and John A. Mellowes Center for Genomic Sciences and Precision Medicine, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | - Brian C. Smith
- Structural Genomics Unit, Linda T. and John A. Mellowes Center for Genomic Sciences and Precision Medicine, Medical College of Wisconsin, Milwaukee, WI 53226, USA
- Department of Biochemistry, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | - Brian F. Volkman
- Structural Genomics Unit, Linda T. and John A. Mellowes Center for Genomic Sciences and Precision Medicine, Medical College of Wisconsin, Milwaukee, WI 53226, USA
- Department of Biochemistry, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | - Raul Urrutia
- Department of Surgery, Medical College of Wisconsin, Milwaukee, WI 53226, USA
- Structural Genomics Unit, Linda T. and John A. Mellowes Center for Genomic Sciences and Precision Medicine, Medical College of Wisconsin, Milwaukee, WI 53226, USA
- Department of Biochemistry, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | - Michael T. Zimmermann
- Bioinformatics Research and Development Laboratory, Linda T. and John A. Mellowes Center for Genomic Sciences and Precision Medicine, Medical College of Wisconsin, Milwaukee, WI 53226, USA
- Department of Biochemistry, Medical College of Wisconsin, Milwaukee, WI 53226, USA
- Clinical and Translational Sciences Institute, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| |
Collapse
|
11
|
Liu C, Ye D, Yang H, Chen X, Su Z, Li X, Ding M, Liu Y. RAS-targeted cancer therapy: Advances in drugging specific mutations. MedComm (Beijing) 2023; 4:e285. [PMID: 37250144 PMCID: PMC10225044 DOI: 10.1002/mco2.285] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Revised: 04/06/2023] [Accepted: 04/18/2023] [Indexed: 05/31/2023] Open
Abstract
Rat sarcoma (RAS), as a frequently mutated oncogene, has been studied as an attractive target for treating RAS-driven cancers for over four decades. However, it is until the recent success of kirsten-RAS (KRAS)G12C inhibitor that RAS gets rid of the title "undruggable". It is worth noting that the therapeutic effect of KRASG12C inhibitors on different RAS allelic mutations or even different cancers with KRASG12C varies significantly. Thus, deep understanding of the characteristics of each allelic RAS mutation will be a prerequisite for developing new RAS inhibitors. In this review, the structural and biochemical features of different RAS mutations are summarized and compared. Besides, the pathological characteristics and treatment responses of different cancers carrying RAS mutations are listed based on clinical reports. In addition, the development of RAS inhibitors, either direct or indirect, that target the downstream components in RAS pathway is summarized as well. Hopefully, this review will broaden our knowledge on RAS-targeting strategies and trigger more intensive studies on exploiting new RAS allele-specific inhibitors.
Collapse
Affiliation(s)
- Cen Liu
- Beijing University of Chinese MedicineBeijingChina
| | - Danyang Ye
- Beijing University of Chinese MedicineBeijingChina
| | - Hongliu Yang
- Beijing University of Chinese MedicineBeijingChina
| | - Xu Chen
- Beijing University of Chinese MedicineBeijingChina
| | - Zhijun Su
- Beijing University of Chinese MedicineBeijingChina
| | - Xia Li
- Institute of Genetics and Developmental BiologyChinese Academy of SciencesBeijingChina
| | - Mei Ding
- Institute of Genetics and Developmental BiologyChinese Academy of SciencesBeijingChina
| | - Yonggang Liu
- Beijing University of Chinese MedicineBeijingChina
| |
Collapse
|
12
|
Samad A, Khurshid B, Mahmood A, Rehman AU, Khalid A, Abdalla AN, Algarni AS, Wadood A. Identification of novel peptide inhibitors for oncogenic KRAS G12D as therapeutic options using mutagenesis-based remodeling and MD simulations. J Biomol Struct Dyn 2023; 41:13425-13437. [PMID: 37010994 DOI: 10.1080/07391102.2023.2192298] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Accepted: 01/22/2023] [Indexed: 04/04/2023]
Abstract
The Kirsten rat sarcoma 2 viral oncogene homolog (KRAS) serves as a molecular switch, cycling between guanosine triphosphate (GTP)-bound and inactive guanosine diphosphate (GDP)-bound states. KRAS modulates numerous signal transduction pathways including the conventional RAF-MEK-ERK pathway. Mutations in the RAS genes have been linked to the formation of malignant tumors. Human malignancies typically show mutations in the Ras gene including HRAS, KRAS, and NRAS. Among all the mutations in exon 12 and exon 13 of the KRAS gene, the G12D mutation is more prevalent in pancreatic and lung cancer and accounts for around 41% of all G12 mutations, making them potential anticancer therapeutic targets. The present study is aimed at repurposing the peptide inhibitor KD2 of the KRAS G12D mutant. We employed an in-silico mutagenesis approach to design novel peptide inhibitors from the experimentally reported peptide inhibitor, and it was found that substitutions (N8W, N8I, and N8Y) might enhance the peptide's binding affinity toward the KRAS. Molecular dynamics simulations and binding energy calculations confirmed that the newly designed peptide inhibitors are stable and that their binding affinities are stronger as compared to the wild-type peptide. The detailed analysis revealed that newly designed peptides have the potential to inhibit KRAS/Raf interaction and the oncogenic signal of the KRAS G12D mutant. Our findings strongly suggest that these peptides should be tested and clinically validated to combat the oncogenic activity of KRAS.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Abdus Samad
- Department of Biochemistry, Abdul Wali Khan University Mardan, Mardan, Pakistan
| | - Beenish Khurshid
- Department of Biochemistry, Abdul Wali Khan University Mardan, Mardan, Pakistan
| | - Arif Mahmood
- Center for Medical Genetics and Hunan Key Laboratory of Medical Genetics, School of Life Sciences, Central South University, Changsha, Hunan, China
| | - Ashfaq Ur Rehman
- Department of Molecular Biology and Biochemistry, School of Biological Sciences, University of California, Irvine, California, USA
| | - Asaad Khalid
- Substance Abuse and Toxicology Research Center, Jazan University, Jazan, Saudi Arabia
- Medicinal and Aromatic Plants and Traditional Medicine Research Institute, National Center for Research, Khartoum, Sudan
| | - Ashraf N Abdalla
- Department of Pharmacology and Toxicology, College of Pharmacy, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Alanood S Algarni
- Department of Pharmacology and Toxicology, College of Pharmacy, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Abdul Wadood
- Department of Biochemistry, Abdul Wali Khan University Mardan, Mardan, Pakistan
| |
Collapse
|
13
|
Shi S, Zheng L, Ren Y, Wang Z. Impacts of Mutations in the P-Loop on Conformational Alterations of KRAS Investigated with Gaussian Accelerated Molecular Dynamics Simulations. Molecules 2023; 28:molecules28072886. [PMID: 37049650 PMCID: PMC10095679 DOI: 10.3390/molecules28072886] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2023] [Revised: 03/19/2023] [Accepted: 03/20/2023] [Indexed: 04/14/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Shuhua Shi
- School of Science, Shandong Jianzhu University, Jinan 250101, China
| | - Linqi Zheng
- School of Science, Shandong Jianzhu University, Jinan 250101, China
| | - Yonglian Ren
- School of Science, Shandong Jianzhu University, Jinan 250101, China
| | - Ziyu Wang
- School of Science, Shandong Jianzhu University, Jinan 250101, China
- Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
| |
Collapse
|
14
|
Wang L, Wang Y, Yu Y, Liu D, Zhao J, Zhang L. Deciphering Selectivity Mechanism of BRD9 and TAF1(2) toward Inhibitors Based on Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations. Molecules 2023; 28:molecules28062583. [PMID: 36985555 PMCID: PMC10052767 DOI: 10.3390/molecules28062583] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Revised: 03/06/2023] [Accepted: 03/09/2023] [Indexed: 03/14/2023] Open
Abstract
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).
Collapse
|
15
|
Wang L, Lu D, Wang Y, Xu X, Zhong P, Yang Z. Binding selectivity-dependent molecular mechanism of inhibitors towards CDK2 and CDK6 investigated by multiple short molecular dynamics and free energy landscapes. J Enzyme Inhib Med Chem 2023; 38:84-99. [DOI: 10.1080/14756366.2022.2135511] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Affiliation(s)
- Lifei Wang
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Dan Lu
- Department of Physics, Jiangxi Agricultural University, Nanchang, PR China
| | - Yan Wang
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Xiaoyan Xu
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Peihua Zhong
- College of Computer Information and Engineering, Jiangxi Agriculture University, Nanchang, PR China
| | - Zhiyong Yang
- Department of Physics, Jiangxi Agricultural University, Nanchang, PR China
| |
Collapse
|
16
|
Bao HY, Wang W, Sun HB, Chen JZ. Binding modes of GDP, GTP and GNP to NRAS deciphered by using Gaussian accelerated molecular dynamics simulations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2023; 34:65-89. [PMID: 36762439 DOI: 10.1080/1062936x.2023.2165542] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/23/2022] [Accepted: 12/31/2022] [Indexed: 06/18/2023]
Abstract
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.
Collapse
Affiliation(s)
- H Y Bao
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - J Z Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
17
|
Chen J, Zeng Q, Wang W, Sun H, Hu G. Decoding the Identification Mechanism of an SAM-III Riboswitch on Ligands through Multiple Independent Gaussian-Accelerated Molecular Dynamics Simulations. J Chem Inf Model 2022; 62:6118-6132. [PMID: 36440874 DOI: 10.1021/acs.jcim.2c00961] [Citation(s) in RCA: 49] [Impact Index Per Article: 24.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
S-Adenosyl-l-methionine (SAM)-responsive riboswitches play a central role in the regulation of bacterial gene expression at the level of transcription attenuation or translation inhibition. In this study, multiple independent Gaussian-accelerated molecular dynamics simulations were performed to decipher the identification mechanisms of SAM-III (SMK) on ligands SAM, SAH, and EEM. The results reveal that ligand binding highly affects the structural flexibility, internal dynamics, and conformational changes of SAM-III. The dynamic analysis shows that helices P3 and P4 as well as two junctions J23 and J24 of SAM-III are highly susceptible to ligand binding. Analyses of free energy landscapes suggest that ligand binding induces different free energy profiles of SAM-III, which leads to the difference in identification sites of SAM-III on ligands. The information on ligand-nucleotide interactions not only uncovers that the π-π, cation-π, and hydrogen bonding interactions drive identification of SAM-III on the three ligands but also reveals that different electrostatic properties of SAM, SAH, and EEM alter the active sites of SAM-III. Meanwhile, the results also verify that the adenine group of SAM, SAH, and EEM is well recognized by conserved nucleotides G7, A29, U37, A38, and G48. We expect that this study can provide useful information for understanding the applications of SAM-III in chemical, synthetic RNA biology, and biomedical fields.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan250357, China
| | - Qingkai Zeng
- School of Science, Shandong Jiaotong University, Jinan250357, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan250357, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan250357, China
| | - Guodong Hu
- Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou253023, China
| |
Collapse
|
18
|
Discovering and Targeting Dynamic Drugging Pockets of Oncogenic Proteins: The Role of Magnesium in Conformational Changes of the G12D Mutated Kirsten Rat Sarcoma-Guanosine Diphosphate Complex. Int J Mol Sci 2022; 23:ijms232213865. [PMID: 36430338 PMCID: PMC9692486 DOI: 10.3390/ijms232213865] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2022] [Revised: 11/02/2022] [Accepted: 11/07/2022] [Indexed: 11/12/2022] Open
Abstract
KRAS-G12D mutations are the one of most frequent oncogenic drivers in human cancers. Unfortunately, no therapeutic agent directly targeting KRAS-G12D has been clinically approved yet, with such mutated species remaining undrugged. Notably, cofactor Mg2+ is closely related to the function of small GTPases, but no investigation has been conducted yet on Mg2+ when associated with KRAS. Herein, through microsecond scale molecular dynamics simulations, we found that Mg2+ plays a crucial role in the conformational changes of the KRAS-GDP complex. We located two brand new druggable dynamic pockets exclusive to KRAS-G12D. Using the structural characteristics of these two dynamic pockets, we designed in silico the inhibitor DBD15-21-22, which can specifically and tightly target the KRAS-G12D-GDP-Mg2+ ternary complex. Overall, we provide two brand new druggable pockets located on KRAS-G12D and suitable strategies for its inhibition.
Collapse
|
19
|
IL-2 Flexible Loops Might Play a Role in IL-2 Interaction with the High-Affinity IL-2 Receptor: A Molecular Dynamics (MD) Study. J CHEM-NY 2022. [DOI: 10.1155/2022/3646375] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
The clinical use of high-dose IL-2 in cancer immunotherapy faces several drawbacks such as toxicity and unfavorable pharmacokinetic profile. These drawbacks can be avoided by inhibiting IL-2 interaction with the CD25 subunit, which is a component of the high-affinity IL-2 receptor (IL-2Rαβγ). Several studies showed mutations of potential IL-2 residues such as R38, F42, Y45, and Y72 would produce IL-2 that is CD25-independent. In essence, structural comparison between wild-type (WT) IL-2 and CD25-independent IL-2 can be very insightful to assess the role of IL-2 flexibility and conformation in the IL-2 receptor interactions. Here, we investigated the flexibility loops and conformation of IL-2m (F24A, Y45A, and L72G), which is known to be CD25-independent, and IL-2m2 (F42Y and L72R) mutants along with WT IL-2 using MD simulations. Despite residue mutations, both IL-2m and IL-2m2 showed comparable conformational compactness and better stability than WT IL-2. Interestingly, IL-2m and IL-2m2 mutants showed rigid BC and CD loops in comparison to WT IL-2 . Also, the AB loop conformation of IL-2m was a bent structure compared to the WT IL-2 and IL-2m2. Principal component analysis (PCA) and free-energy landscape results suggested IL-2m and IL-2m2 have stable conformations compared to the WT IL-2. Therefore, these mutation sites of IL-2 produced stable and rigid loops that might prevent IL-2 from binding to the CD25 subunit. Our results can help to assess IL-2 flexibility loops to design new CD25-independent IL-2 mutants without compromising the IL-2 structure.
Collapse
|
20
|
Yu YX, Wang W, Sun HB, Zhang LL, Wang LF, Yin YY. 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. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2022; 33:805-831. [PMID: 36322686 DOI: 10.1080/1062936x.2022.2140708] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/08/2022] [Accepted: 10/20/2022] [Indexed: 06/16/2023]
Abstract
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.
Collapse
Affiliation(s)
- Y X Yu
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L L Zhang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L F Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Y Y Yin
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
21
|
Yu Z, Su H, Chen J, Hu G. Deciphering Conformational Changes of the GDP-Bound NRAS Induced by Mutations G13D, Q61R, and C118S through Gaussian Accelerated Molecular Dynamic Simulations. Molecules 2022; 27:molecules27175596. [PMID: 36080363 PMCID: PMC9457619 DOI: 10.3390/molecules27175596] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Revised: 08/24/2022] [Accepted: 08/25/2022] [Indexed: 11/29/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Zhiping Yu
- Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
| | - Hongyi Su
- Laoling People's Hospital, Dezhou 253600, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Guodong Hu
- Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
- Laoling People's Hospital, Dezhou 253600, China
| |
Collapse
|
22
|
Chen J, Wang J, Zeng Q, Wang W, Sun H, Wei B. Exploring the deactivation mechanism of human β2 adrenergic receptor by accelerated molecular dynamic simulations. Front Mol Biosci 2022; 9:972463. [PMID: 36111136 PMCID: PMC9468641 DOI: 10.3389/fmolb.2022.972463] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2022] [Accepted: 08/05/2022] [Indexed: 11/13/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
- *Correspondence: Jianzhong Chen, ; Benzheng Wei,
| | - Jian Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Qingkai Zeng
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Benzheng Wei
- Center for Medical Artificial Intelligence, Shandong University of Traditional Chinese Medicine, Qingdao, China
- *Correspondence: Jianzhong Chen, ; Benzheng Wei,
| |
Collapse
|
23
|
Mhashal AR, Yoluk O, Orellana L. Exploring the Conformational Impact of Glycine Receptor TM1-2 Mutations Through Coarse-Grained Analysis and Atomistic Simulations. Front Mol Biosci 2022; 9:890851. [PMID: 35836931 PMCID: PMC9275627 DOI: 10.3389/fmolb.2022.890851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2022] [Accepted: 05/05/2022] [Indexed: 11/13/2022] Open
Abstract
Pentameric ligand-gated ion channels (PLGICs) are a family of proteins that convert chemical signals into ion fluxes through cellular membranes. Their structures are highly conserved across all kingdoms from bacteria to eukaryotes. Beyond their classical roles in neurotransmission and neurological disorders, PLGICs have been recently related to cell proliferation and cancer. Here, we focus on the best characterized eukaryotic channel, the glycine receptor (GlyR), to investigate its mutational patterns in genomic-wide tumor screens and compare them with mutations linked to hyperekplexia (HPX), a Mendelian neuromotor disease that disrupts glycinergic currents. Our analysis highlights that cancer mutations significantly accumulate across TM1 and TM2, partially overlapping with HPX changes. Based on 3D-clustering, conservation, and phenotypic data, we select three mutations near the pore, expected to impact GlyR conformation, for further study by molecular dynamics (MD). Using principal components from experimental GlyR ensembles as framework, we explore the motions involved in transitions from the human closed and desensitized structures and how they are perturbed by mutations. Our MD simulations show that WT GlyR spontaneously explores opening and re-sensitization transitions that are significantly impaired by mutations, resulting in receptors with altered permeability and desensitization properties in agreement with HPX functional data.
Collapse
|
24
|
Chen J, Zhang S, Zeng Q, Wang W, Zhang Q, Liu X. Free Energy Profiles Relating With Conformational Transition of the Switch Domains Induced by G12 Mutations in GTP-Bound KRAS. Front Mol Biosci 2022; 9:912518. [PMID: 35586192 PMCID: PMC9108337 DOI: 10.3389/fmolb.2022.912518] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2022] [Accepted: 04/15/2022] [Indexed: 01/21/2023] Open
Abstract
Mutations of G12 in KRAS have been involved in different cancers. Multiple replica-Gaussian accelerated molecular dynamics (MR-GaMD) simulations are applied to investigate conformational changes of the switch domains caused by G12C, G12D and G12R. Free energy landscapes suggest that G12C, G12D and G12R induce more energetic states compared to the GTP-bound WT KRAS and make the conformations of the switch domains more disordered, which disturbs bindings of KRAS to effectors. Dynamics analyses based on MR-GaMD trajectory show that G12C, G12D and G12R not only change structural flexibility of the switch domains but also affect their motion behavior, indicating that these three mutations can be used to tune the activity of KRAS. The analyses of interaction networks verify that the instability in interactions of the GTP with the switch SⅠ plays an important role in the high disorder states of the switch domain. This work is expected to provide useful information for deeply understanding the function of KRAS.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
- *Correspondence: Jianzhong Chen, ; Xinguo Liu,
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Qingkai Zeng
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, China
- *Correspondence: Jianzhong Chen, ; Xinguo Liu,
| |
Collapse
|
25
|
Liang S, Liu X, Zhang S, Li M, Zhang Q, Chen J. Binding mechanism of inhibitors to SARS-CoV-2 main protease deciphered by multiple replica molecular dynamics simulations. Phys Chem Chem Phys 2022; 24:1743-1759. [PMID: 34985081 DOI: 10.1039/d1cp04361g] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The outbreak caused by SARS-CoV-2 has received extensive worldwide attention. As the main protease (Mpro) in SARS-CoV-2 has no human homologues, it is feasible to reduce the possibility of targeting the host protein by accidental drugs. Thus, Mpro 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 Mpro. The results indicate that the binding of four inhibitors clearly affects the structural flexibility and internal dynamics of Mpro 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 Mpro is favorable for inhibitor binding. Residue-based free energy decomposition reveals that the inhibitor-Mpro interaction networks involving hydrogen bonding interactions and hydrophobic interactions provide significant information for the design of potent inhibitors against Mpro. 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 Mpro.
Collapse
Affiliation(s)
- Shanshan Liang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Meng Li
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250357, China.
| |
Collapse
|
26
|
Chen J, Zeng Q, Wang W, Hu Q, Bao H. Q61 mutant-mediated dynamics changes of the GTP-KRAS complex probed by Gaussian accelerated molecular dynamics and free energy landscapes. RSC Adv 2022; 12:1742-1757. [PMID: 35425180 PMCID: PMC8978876 DOI: 10.1039/d1ra07936k] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2021] [Accepted: 01/04/2022] [Indexed: 12/19/2022] Open
Abstract
Understanding the molecular mechanism of the GTP-KRAS binding is significant for improving the target roles of KRAS in cancer treatment. In this work, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations were applied to decode the effect of Q61A, Q61H and Q61L on the activity of KRAS. Dynamics analyses based on MR-GaMD trajectory reveal that motion modes and dynamics behavior of the switch domain in KRAS are heavily affected by the three Q61 mutants. Information of free energy landscapes (FELs) shows that Q61A, Q61H and Q61L induce structural disorder of the switch domain and disturb the activity of KRAS. Analysis of the interaction network uncovers that the decrease in the stability of hydrogen bonding interactions (HBIs) of GTP with residues V29 and D30 induced by Q61A, Q61H and Q61L is responsible for the structural disorder of the switch-I and that in the occupancy of the hydrogen bond between GTP and residue G60 leads to the structural disorder of the switch-II. Thus, the high disorder of the switch domain caused by three current Q61 mutants produces a significant effect on binding of KRAS to its effectors. This work is expected to provide useful information for further understanding function and target roles of KRAS in anti-cancer drug development.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Qingkai Zeng
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Wei Wang
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Qingquan Hu
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Huayin Bao
- School of Pharmacy, Shandong University of Traditional Chinese Medicine Jinan 250355 China
| |
Collapse
|
27
|
Wang L, Wang Y, Zhao J, Yu Y, Kang N, Yang Z. Theoretical exploration of the binding selectivity of inhibitors to BRD7 and BRD9 with multiple short molecular dynamics simulations. RSC Adv 2022; 12:16663-16676. [PMID: 35754900 PMCID: PMC9169554 DOI: 10.1039/d2ra02637f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2022] [Accepted: 05/29/2022] [Indexed: 12/18/2022] Open
Abstract
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.![]()
Collapse
Affiliation(s)
- Lifei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Yan Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Juan Zhao
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Yingxia Yu
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Nianqian Kang
- Department of Physics, Jiangxi Agricultural University, Nanchang 330045, China
| | - Zhiyong Yang
- Department of Physics, Jiangxi Agricultural University, Nanchang 330045, China
| |
Collapse
|
28
|
Li M, Liu X, Zhang S, Liang S, Zhang Q, Chen J. Deciphering binding mechanism of inhibitors to SARS-COV-2 main protease through multiple replica accelerated molecular dynamics simulations and free energy landscapes. Phys Chem Chem Phys 2022; 24:22129-22143. [DOI: 10.1039/d2cp03446h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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...
Collapse
|
29
|
Yu YX, Liu WT, Li HY, Wang W, Sun HB, Zhang LL, Wu SL. Decoding molecular mechanism underlying binding of drugs to HIV-1 protease with molecular dynamics simulations and MM-GBSA calculations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2021; 32:889-915. [PMID: 34551634 DOI: 10.1080/1062936x.2021.1979647] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/21/2021] [Accepted: 09/08/2021] [Indexed: 06/13/2023]
Abstract
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.
Collapse
Affiliation(s)
- Y X Yu
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W T Liu
- Shuifa Qilu Cultural Tourism Development Co., Ltd, Shuifa Ecological Industry Group, Jinan, China
| | - H Y Li
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L L Zhang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - S L Wu
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|