Q61 mutant-mediated dynamics changes of the GTP-KRAS complex probed by Gaussian accelerated molecular dynamics and free energy landscapes

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.

Expanding the Genetic Code of Xenopus laevis Embryos

The incorporation of unnatural amino acids into proteins through genetic code expansion has been successfully adapted to African claw-toed frog embryos. Six unique unnatural amino acids are incorporated site-specifically into proteins and demonstrate robust and reliable protein expression. Of these amino acids, several are caged analogues that can be used to establish conditional control over enzymatic activity. Using light or small molecule triggers, we exhibit activation and tunability of protein functions in live embryos. This approach was then applied to optical control over the activity of a RASopathy mutant of NRAS, taking advantage of generating explant cultures from Xenopus. Taken together, genetic code expansion is a robust approach in the Xenopus model to incorporate novel chemical functionalities into proteins of interest to study their function and role in a complex biological setting.

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.

Conformational States of the GDP- and GTP-Bound HRAS Affected by A59E and K117R: An Exploration from Gaussian Accelerated Molecular Dynamics

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.

GTP-Bound N-Ras Conformational States and Substates Are Modulated by Membrane and Point Mutation

Oncogenic Ras proteins are known to present multiple conformational states, as reported by the great variety of crystallographic structures. The GTP-bound states are grouped into two main states: the "inactive" state 1 and the "active" state 2. Recent reports on H-Ras have shown that state 2 exhibits two substates, directly related to the orientation of Tyr32: toward the GTP-bound pocket and outwards. In this paper, we show that N-Ras exhibits another substate of state 2, related to a third orientation of Tyr32, toward Ala18 and parallel to the GTP-bound pocket. We also show that this substate is highly sampled in the G12V mutation of N-Ras and barely present in its wild-type form, and that the G12V mutation prohibits the sampling of the GTPase-activating protein (GAP) binding substate, rendering this mutation oncogenic. Furthermore, using molecular dynamics simulations, we explore the importance of the membrane on N-Ras' conformational state dynamics and its strong influence on Ras protein stability. Moreover, the membrane has a significant influence on the conformational (sub)states sampling of Ras. This, in turn, is of crucial importance in the activation/deactivation cycle of Ras, due to the binding of guanine nucleotide exchange factor proteins (GEFs)/GTPase-activating proteins (GAPs).

Binding selectivity analysis of AURKs inhibitors through molecular dynamics simulation studies

Aurora kinases (AURKs) have been identified as promising biological targets for the treatment of cancer. In this study, molecular dynamics simulations were employed to investigate the binding selectivity of three inhibitors (HPM, MPY, and VX6) towards AURKA and AURKB by predicting their binding free energies. The results show that the inhibitors HPM, MPY, and VX6 have more favorable interactions with AURKB as compared to AURKA. The binding energy decomposition analysis revealed that four common residue pairs (L139, L83), (V147, V91), (L210, L154), and (L263, L207) showed significant binding energies with HPM, MPY, and VX6, hence responsible for the binding selectivity of AURKA and AURKB to the inhibitors. The MD trajectory analysis also revealed that the inhibitors affect the dynamic flexibility of protein structure, which is also responsible for the partial selectivity of HPM, MPY, and VX6 towards AURKA and AURKB. As expected, this study provides useful insights for the design of potential inhibitors with high selectivity for AURKA and AURKB.

Mechanism-Based Redesign of GAP to Activate Oncogenic Ras

Ras GTPases play a crucial role in cell signaling pathways. Mutations of the Ras gene occur in about one third of cancerous cell lines and are often associated with detrimental clinical prognosis. Hot spot residues Gly12, Gly13, and Gln61 cover 97% of oncogenic mutations, which impair the enzymatic activity in Ras. Using QM/MM free energy calculations, we present a two-step mechanism for the GTP hydrolysis catalyzed by the wild-type Ras.GAP complex. We found that the deprotonation of the catalytic water takes place via the Gln61 as a transient Brønsted base. We also determined the reaction profiles for key oncogenic Ras mutants G12D and G12C using QM/MM minimizations, matching the experimentally observed loss of catalytic activity, thereby validating our reaction mechanism. Using the optimized reaction paths, we devised a fast and accurate procedure to design GAP mutants that activate G12D Ras. We replaced GAP residues near the active site and determined the activation barrier for 190 single mutants. We furthermore built a machine learning for ultrafast screening, by fast prediction of the barrier heights, tested both on the single and double mutations. This work demonstrates that fast and accurate screening can be accomplished via QM/MM reaction path optimizations to design protein sequences with increased catalytic activity. Several GAP mutations are predicted to re-enable catalysis in oncogenic G12D, offering a promising avenue to overcome aberrant Ras-driven signal transduction by activating enzymatic activity instead of inhibition. The outlined computational screening protocol is readily applicable for designing ligands and cofactors analogously.

Genetically Encoded Aminocoumarin Lysine for Optical Control of Protein-Nucleotide Interactions in Zebrafish Embryos

The strategic placement of unnatural amino acids into the active site of kinases and phosphatases has allowed for the generation of photocaged signaling proteins that offer spatiotemporal control over activation of these pathways through precise light exposure. However, deploying this technology to study cell signaling in the context of embryo development has been limited. The promise of optical control is especially useful in the early stages of an embryo where development is driven by tightly orchestrated signaling events. Here, we demonstrate light-induced activation of Protein Kinase A and a RASopathy mutant of NRAS in the zebrafish embryo using a new light-activated amino acid. We applied this approach to gain insight into the roles of these proteins in gastrulation and heart development and forge a path for further investigation of RASopathy mutant proteins in animals.

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

Engineering Small Molecule Switches of Protein Function in Zebrafish Embryos

Precise temporally regulated protein function directs the highly complex processes that make up embryo development. The zebrafish embryo is an excellent model organism to study development, and conditional control over enzymatic activity is desirable to target chemical intervention to specific developmental events and to investigate biological mechanisms. Surprisingly few, generally applicable small molecule switches of protein function exist in zebrafish. Genetic code expansion allows for site-specific incorporation of unnatural amino acids into proteins that contain caging groups that are removed through addition of small molecule triggers such as phosphines or tetrazines. This broadly applicable control of protein function was applied to activate several enzymes, including a GTPase and a protease, with temporal precision in zebrafish embryos. Simple addition of the small molecule to the media produces robust and tunable protein activation, which was used to gain insight into the development of a congenital heart defect from a RASopathy mutant of NRAS and to control DNA and protein cleavage events catalyzed by a viral recombinase and a viral protease, respectively.

Comprehensive evaluation of end-point free energy techniques in carboxylated-pillar[6]arene host-guest binding: II. regression and dielectric constant

End-point free energy calculations as a powerful tool have been widely applied in protein-ligand and protein-protein interactions. It is often recognized that these end-point techniques serve as an option of intermediate accuracy and computational cost compared with more rigorous statistical mechanic models (e.g., alchemical transformation) and coarser molecular docking. However, it is observed that this intermediate level of accuracy does not hold in relatively simple and prototypical host-guest systems. Specifically, in our previous work investigating a set of carboxylated-pillar[6]arene host-guest complexes, end-point methods provide free energy estimates deviating significantly from the experimental reference, and the rank of binding affinities is also incorrectly computed. These observations suggest the unsuitability and inapplicability of standard end-point free energy techniques in host-guest systems, and alteration and development are required to make them practically usable. In this work, we consider two ways to improve the performance of end-point techniques. The first one is the PBSA_E regression that varies the weights of different free energy terms in the end-point calculation procedure, while the second one is considering the interior dielectric constant as an additional variable in the end-point equation. By detailed investigation of the calculation procedure and the simulation outcome, we prove that these two treatments (i.e., regression and dielectric constant) are manipulating the end-point equation in a somehow similar way, i.e., weakening the electrostatic contribution and strengthening the non-polar terms, although there are still many detailed differences between these two methods. With the trained end-point scheme, the RMSE of the computed affinities is improved from the standard ~ 12 kcal/mol to ~ 2.4 kcal/mol, which is comparable to another altered end-point method (ELIE) trained with system-specific data. By tuning PBSA_E weighting factors with the host-specific data, it is possible to further decrease the prediction error to ~ 2.1 kcal/mol. These observations along with the extremely efficient optimized-structure computation procedure suggest the regression (i.e., PBSA_E as well as its GBSA_E extension) as a practically applicable solution that brings end-point methods back into the library of usable tools for host-guest binding. However, the dielectric-constant-variable scheme cannot effectively minimize the experiment-calculation discrepancy for absolute binding affinities, but is able to improve the calculation of affinity ranks. This phenomenon is somehow different from the protein-ligand case and suggests the difference between host-guest and biomacromolecular (protein-ligand and protein-protein) systems. Therefore, the spectrum of tools usable for protein-ligand complexes could be unsuitable for host-guest binding, and numerical validations are necessary to screen out really workable solutions in these 'prototypical' situations.

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.

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

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.

Conformations and binding pockets of HRas and its guanine nucleotide exchange factors complexes in the guanosine triphosphate exchange process

The human Son of Sevenless (SOS) activates the signal-transduction protein Ras by forming the complex SOS·Ras and accelerating the guanosine triphosphate (GTP) exchange in Ras. Inhibition of SOS·Ras could regulate the function of Ras in cells and has emerged as an effective strategy for battling Ras related cancers. A key factor to the success of this approach is to understand the conformational change of Ras during the GTP exchange process. In this study, we perform an extensive molecular dynamics simulation to characterize the specific conformations of Ras without and with guanine nucleotide exchange factors (GEFs) of SOS, especially for the substates of State 1 of HRasGTP∙Mg²⁺ . The potent binding pockets on the surfaces of the RasGDP∙Mg²⁺ , the S1.1 and S1.2 substates in State 1 of RasGTP∙Mg²⁺ and the ternary complexes with SOS are predicted, including the binding sites of other domains of SOS. These findings help to obtain a more thorough understanding of Ras functions in the GTP cycling process and provide a structural foundation for future drug design.

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.