New paradigms in actomyosin energy transduction: Critical evaluation of non-traditional models for orthophosphate release.

Release of the ATP hydrolysis product ortophosphate (Pi) from the active site of myosin is central in chemo-mechanical energy transduction and closely associated with the main force-generating structural change, the power-stroke. Despite intense investigations, the relative timing between Pi-release and the power-stroke remains poorly understood. This hampers in depth understanding of force production by myosin in health and disease and our understanding of myosin-active drugs. Since the 1990s and up to today, models that incorporate the Pi-release either distinctly before or after the power-stroke, in unbranched kinetic schemes, have dominated the literature. However, in recent years, alternative models have emerged to explain apparently contradictory findings. Here, we first compare and critically analyze three influential alternative models proposed previously. These are either characterized by a branched kinetic scheme or by partial uncoupling of Pi-release and the power-stroke. Finally, we suggest critical tests of the models aiming for a unified picture.

Trans vs. cis: a computational study of enasidenib resistance due to IDH2 mutations.

Isocitrate dehydrogenase 2 (IDH2) is a homodimeric enzyme that plays an important role in energy production. A mutation R140Q in one monomer makes the enzyme tumourigenic. Enasidenib is an effective inhibitor of IDH2/R140Q. A secondary mutation Q316E leads to enasidenib resistance. This mutation was hitherto only found in trans, i.e. where one monomer has the R140Q mutation and the other carries the Q316E mutation. It is not clear if the mutation only leads to resistance when in trans or if it has been discovered in trans only by chance, since it was only reported in two patients. Using molecular dynamics (MD) simulations we show that the binding of enasidenib to IDH2 is indeed much weaker when the Q316E mutation takes place in trans not in cis, which provides a molecular explanation for the clinical finding. This is corroborated by non-covalent interaction (NCI) analysis and DFT calculations. Whereas the MD simulations show a loss of one hydrogen bond upon the resistance mutation, NCI and energy decomposition analysis (EDA) reveal that a multitude of interactions are weakened.

Structural and thermodynamic characterization of allosteric transitions in human serum albumin with metadynamics simulations.

Human serum albumin (HSA) is the most prominent protein in blood plasma, responsible for the maintenance of blood viscosity and transport of endogenous and exogenous molecules. Fatty acids (FA) are the most common ligands of HSA and their binding can modify the protein's structure. The protein can assume two well-defined conformations, referred to as 'Neutral' and 'Basic'. The Neutral (N) state occurs at pH close to 7.0 and in the absence of bound FA. The Basic (B) state occurs at pH higher than 8.0 or when the protein is bound to long-chain FA. HSA's allosteric behaviour is dependent on the number on FA bound to the structure. However, the mechanism of this allosteric regulation is not clear. To understand how albumin changes its conformation, we compared a series of HSA structures deposited in the protein data bank to identify the minimum amount of FA bound to albumin, which is enough to drive the allosteric transition. Thereafter, non-biased molecular dynamics (MD) simulations were used to track protein's dynamics. Surprisingly, running an ensemble of relatively short MD simulations, we observed rapid transition from the B to the N state. These simulations revealed differences in the mobilities of the protein's subdomains, with one domain unable to fully complete its transition. To track the transition dynamics in full, we used these results to choose good geometrical collective variables for running metadynamics simulations. The metadynamics calculations showed that there was a low energy barrier for the transition from the B to the N state, while a higher energy barrier was observed for the N to the B transition. These calculations also offered valuable insights into the transition process.

QM/MM-Based Energy Decomposition Analysis Method for Large Systems.

In this work, a QM/MM-based EDA method, called GKS-EDA(QM/MM), is proposed. As an extension of GKS-EDA, this scheme divides the total interaction energy into electrostatic, exchange-repulsion, polarization, and correlation/dispersion terms. GKS-EDA(QM/MM) can be applied to describe the interactions of large-scale systems combined with various QM/MM platforms. By using the examples of a hydrated hydronium ion complex in water solution, the barnase-barstar complex, and MMP-13-pyrimidinetrione in a metalloprotein, the capability of GKS-EDA(QM/MM) for various interactions in large systems is validated.

Allosteric enhancement of the BCR-Abl1 kinase inhibition activity of nilotinib by cobinding of asciminib.

Inhibitors that bind competitively to the ATP binding pocket in the kinase domain of the oncogenic fusion protein BCR-Abl1 are used successfully in targeted therapy of chronic myeloid leukemia (CML). Such inhibitors provided the first proof of concept that kinase inhibition can succeed in a clinical setting. However, emergence of drug resistance and dose-dependent toxicities limit the effectiveness of these drugs. Therefore, treatment with a combination of drugs without overlapping resistance mechanisms appears to be an appropriate strategy. In the present work, we explore the effectiveness of combination therapies of the recently developed allosteric inhibitor asciminib with the ATP-competitive inhibitors nilotinib and dasatinib in inhibiting the BCR-Abl1 kinase activity in CML cell lines. Through these experiments, we demonstrate that asciminib significantly enhances the inhibition activity of nilotinib, but not of dasatinib. Exploring molecular mechanisms for such allosteric enhancement via systematic computational investigation incorporating molecular dynamics, metadynamics simulations, and density functional theory calculations, we found two distinct contributions. First, binding of asciminib triggers conformational changes in the inactive state of the protein, thereby making the activation process less favorable by ∼4 kcal/mol. Second, the binding of asciminib decreases the binding free energies of nilotinib by ∼3 and ∼7 kcal/mol for the wildtype and T315I-mutated protein, respectively, suggesting the possibility of reducing nilotinib dosage and lowering risk of developing resistance in the treatment of CML.

Combination strategies to overcome drug resistance in FLT+ acute myeloid leukaemia.

BACKGROUND: Acute myeloid leukaemia (AML) remains difficult to treat despite the development of novel formulations and targeted therapies. Activating mutations in the FLT3 gene are common among patients and make the tumour susceptible to FLT3 inhibitors, but resistance to such inhibitors develops quickly. METHODS: We examined combination therapies aimed at FLT3+-AML, and studied the development of resistance using a newly developed protocol. Combinations of FLT3, CDK4/6 and PI3K inhibitors were tested for synergism. RESULTS: We show that AML cells express CDK4 and that the CDK4/6 inhibitors palbociclib and abemaciclib inhibit cellular growth. PI3K inhibitors were also effective in inhibiting the growth of AML cell lines that express FLT3-ITD. Whereas resistance to quizartinib develops quickly, the combinations overcome such resistance. CONCLUSIONS: This study suggests that a multi-targeted intervention involving a CDK4/6 inhibitor with a FLT3 inhibitor or a pan-PI3K inhibitor might be a valuable therapeutic strategy for AML to overcome drug resistance. Moreover, many patients cannot tolerate high doses of the drugs that were studied (quizartinib, palbociclib and PI3K inhibitors) for longer periods, and it is therefore of high significance that the drugs act synergistically and lower doses can be used.

Estimating the Gibbs Hydration Energies of Actinium and Trans-Plutonium Actinides.

The use of actinides for medical, scientific and technological purposes has gained momentum in the recent years. This creates a need to understand their interactions with biomolecules, both at the interface and as they become complexed. Calculation of the Gibbs binding energies of the ions to biomolecules, i. e., the Gibbs energy change associated with a transfer of an ion from the water phase to its binding site, could help to understand the actinides' toxicities and to design agents that bind them with high affinities. To this end, there is a need to obtain accurate reference values for actinide hydration, that for most actinides are not available from experiment. In this study, a set of ionic radii is developed that enables future calculations of binding energies for Pu3+ and five actinides with renewed scientific and technological interest: Ac3+ , Am3+ , Cm3+ , Bk3+ and Cf3+ . Reference hydration energies were calculated using quantum chemistry and ion solvation theory and agree well for all ions except Ac3+ , where ion solvation theory seems to underestimate the magnitude of the Gibbs hydration energy. The set of radii and reference energies that are presented here provide means to calculate binding energies for actinides and biomolecules.

Multiscale QM/MM Simulations Identify the Roles of Asp239 and 1-OH···Nucleophile in Transition State Stabilization in Arabidopsis thaliana Cell-Wall Invertase 1.

Arabidopsis thaliana cell-wall invertase 1 (AtCWIN1), a key enzyme in sucrose metabolism in plants, catalyzes the hydrolysis of sucrose into fructose and glucose. AtCWIN1 belongs to the glycoside hydrolase GH-J clan, where two carboxylate residues (Asp23 and Glu203 in AtCWIN1) are well documented as a nucleophile and an acid/base catalyst. However, details at the atomic level about the role of neighboring residues and enzyme-substrate interactions during catalysis are not fully understood. Here, quantum mechanical/molecular mechanical (QM/MM) free-energy simulations were carried out to clarify the origin of the observed decreased rates in Asp239Ala, Asp239Asn, and Asp239Phe in AtCWIN1 compared to the wild type and delineate the role of Asp239 in catalysis. The glycosylation and deglycosylation steps were considered in both wild type and mutants. Deglycosylation is predicted to be the rate-determining step in the reaction, with a calculated overall free-energy barrier of 15.9 kcal/mol, consistent with the experimental barrier (15.3 kcal/mol). During the reaction, the -1 furanosyl ring underwent a conformational change corresponding to 3E ↔ [E2]⧧ ↔ 1E according to the nomenclature of saccharide structures along the full catalytic reaction. Asp239 was found to stabilize not only the transition state but also the fructosyl-enzyme intermediate, which explains findings from previous structural and mutagenesis experiments. The 1-OH···nucleophile interaction has been found to provide an important contribution to the transition state stabilization, with a contribution of ∼7 kcal/mol, and affected glycosylation more significantly than deglycosylation. This study provides molecular insights that improve the current understanding of sucrose binding and hydrolysis in members of clan GH-J, which may benefit protein engineering research. Finally, a rationale on the sucrose inhibitor configuration in chicory 1-FEH IIa, proposed a long time ago in the literature, is also provided based on the QM/MM calculations.

Resistance to a tyrosine kinase inhibitor mediated by changes to the conformation space of the kinase.

Gilteritinib is a highly selective and effective inhibitor of the FLT3/ITD mutated protein, and is used successfully in treating acute myeloid leukaemia (AML). Unfortunately, tumour cells gradually develop resistance to gilteritinib due to mutations in the molecular drug target. The atomistic details behind this observed resistance are not clear, since the protein structure of the complex is only available in the inactive state, while the drug binds better to the active state. To overcome this limitation, we used a computer-aided approach where we docked gilteritinib to the active site of FLT3/ITD and calculated the Gibbs free energy difference between the binding energies of the parental and mutant enzymes. These calculations agreed with experimental estimations for one mutation (F691L) but not the other (D698N). To further understand how these mutations operate, we used metadynamics simulations to study the conformational landscape of the activation process. Both mutants show a lower activation energy barrier which suggests that they are more likely to adopt an active state until inhibited, making the mutant enzymes more active. This suggests that a higher efficiency of tyrosine kinases contributes to resistance not only against type 2 but also against type 1 kinase inhibitors.

A general tight-binding based energy decomposition analysis scheme for intermolecular interactions in large molecules.

In this work, a general tight-binding based energy decomposition analysis (EDA) scheme for intermolecular interactions is proposed. Different from the earlier version [Xu et al., J. Chem. Phys. 154, 194106 (2021)], the current tight-binding based density functional theory (DFTB)-EDA is capable of performing interaction analysis with all the self-consistent charge (SCC) type DFTB methods, including SCC-DFTB2/3 and GFN1/2-xTB, despite their different formulas and parameterization schemes. In DFTB-EDA, the total interaction energy is divided into frozen, polarization, and dispersion terms. The performance of DFTB-EDA with SCC-DFTB2/3 and GFN1/2-xTB for various interaction systems is discussed and assessed.

Combating drug resistance in acute myeloid leukaemia by drug rotations: the effects of quizartinib and pexidartinib.

BACKGROUND: Acute myeloid leukaemia (AML) is an aggressive blood cancer. In approximately 30% of the cases, driver mutations in the FLT3 gene are identified. FLT3 inhibitors are used in treatment of such patients together with cytotoxic drugs or (in refractory AML) as single agents. Unfortunately, resistance to FLT3 inhibitors limits their efficacy. Resistance is often due to secondary mutations in the gene encoding the molecular target. The gatekeeper mutation F691L confers resistance to specific FLT3 inhibitors such as quizartinib, but pexidartinib is much less resistance to this mutation. Pexidartinib alone is however sensitive to many other resistance mutations. In chronic myeloid leukaemia (CML), it has been suggested that rotation between drugs with a different landscape of resistance mutations might postpone the emergence of resistance. METHODS: We studied the effect of quizartinib and pexidartinib in AML cell lines that express FLT3 (MOLM-14 and MV4-11). Using a rotation protocol, we further examined whether the emergence of resistance could be postponed. Computational modelling was used to analyse the onset of resistance and suggest which mutations are most likely to occur in a quantitative fashion. RESULTS: The cells were sensitive to both inhibitors but quickly developed resistance that could be inherited, suggesting a genetic origin. Rotation protocols were not useful to postpone the emergence of resistance, which implies that such protocols, or changing from pexidartinib to quizartinib (or vice-versa) should not be used in patients. The computational modelling led to similar conclusions and suggested that F691L is the most common mutation to occur with quizartinib, and also when both drugs are used in rotation. CONCLUSIONS: AML patients are not likely to benefit from a quizartinib/pexidartinib rotation protocol. A combination of tyrosine kinase inhibitors (with different molecular targets) might be more useful in the future. Development of specific FLT3 inhibitors that are less sensitive to resistance mutations might also lead to a better outcome.

Understanding intermolecular interactions of large systems in ground state and excited state by using density functional based tight binding methods.

A novel energy decomposition analysis scheme, named DFTB-EDA, is proposed based on the density functional based tight-binding method (DFTB/TD-DFTB), which is a semi-empirical quantum mechanical method based on Kohn-Sham-DFT for large-scale calculations. In DFTB-EDA, the total interaction energy is divided into three terms: frozen density, polarization, and dispersion. Owing to the small cost of DFTB/TD-DFTB, DFTB-EDA is capable of analyzing intermolecular interactions in large molecular systems containing several thousand atoms with high computational efficiency. It can be used not only for ground states but also for excited states. Test calculations, involving the S66 and L7 databases, several large molecules, and non-covalent bonding complexes in their lowest excited states, demonstrate the efficiency, usefulness, and capabilities of DFTB-EDA. Finally, the limits of DFTB-EDA are pointed out.

Inferring time-dependent population growth rates in cell cultures undergoing adaptation.

BACKGROUND: The population growth rate is an important characteristic of any cell culture. During sustained experiments, the growth rate may vary due to competition or adaptation. For instance, in presence of a toxin or a drug, an increasing growth rate indicates that the cells adapt and become resistant. Consequently, time-dependent growth rates are fundamental to follow on the adaptation of cells to a changing evolutionary landscape. However, as there are no tools to calculate the time-dependent growth rate directly by cell counting, it is common to use only end point measurements of growth rather than tracking the growth rate continuously. RESULTS: We present a computer program for inferring the growth rate over time in suspension cells using nothing but cell counts, which can be measured non-destructively. The program was tested on simulated and experimental data. Changes were observed in the initial and absolute growth rates, betraying resistance and adaptation. CONCLUSIONS: For experiments where adaptation is expected to occur over a longer time, our method provides a means of tracking growth rates using data that is normally collected anyhow for monitoring purposes. The program and its documentation are freely available at https://github.com/Sandalmoth/ratrack under the permissive zlib license.

Reverse engineering directed gene regulatory networks from transcriptomics and proteomics data of biomining bacterial communities with approximate Bayesian computation and steady-state signalling simulations.

BACKGROUND: Network inference is an important aim of systems biology. It enables the transformation of OMICs datasets into biological knowledge. It consists of reverse engineering gene regulatory networks from OMICs data, such as RNAseq or mass spectrometry-based proteomics data, through computational methods. This approach allows to identify signalling pathways involved in specific biological functions. The ability to infer causality in gene regulatory networks, in addition to correlation, is crucial for several modelling approaches and allows targeted control in biotechnology applications. METHODS: We performed simulations according to the approximate Bayesian computation method, where the core model consisted of a steady-state simulation algorithm used to study gene regulatory networks in systems for which a limited level of details is available. The simulations outcome was compared to experimentally measured transcriptomics and proteomics data through approximate Bayesian computation. RESULTS: The structure of small gene regulatory networks responsible for the regulation of biological functions involved in biomining were inferred from multi OMICs data of mixed bacterial cultures. Several causal inter- and intraspecies interactions were inferred between genes coding for proteins involved in the biomining process, such as heavy metal transport, DNA damage, replication and repair, and membrane biogenesis. The method also provided indications for the role of several uncharacterized proteins by the inferred connection in their network context. CONCLUSIONS: The combination of fast algorithms with high-performance computing allowed the simulation of a multitude of gene regulatory networks and their comparison to experimentally measured OMICs data through approximate Bayesian computation, enabling the probabilistic inference of causality in gene regulatory networks of a multispecies bacterial system involved in biomining without need of single-cell or multiple perturbation experiments. This information can be used to influence biological functions and control specific processes in biotechnology applications.

The effects of combination treatments on drug resistance in chronic myeloid leukaemia: an evaluation of the tyrosine kinase inhibitors axitinib and asciminib.

BACKGROUND: Chronic myeloid leukaemia is in principle a treatable malignancy but drug resistance is lowering survival. Recent drug discoveries have opened up new options for drug combinations, which is a concept used in other areas for preventing drug resistance. Two of these are (I) Axitinib, which inhibits the T315I mutation of BCR-ABL1, a main source of drug resistance, and (II) Asciminib, which has been developed as an allosteric BCR-ABL1 inhibitor, targeting an entirely different binding site, and as such does not compete for binding with other drugs. These drugs offer new treatment options. METHODS: We measured the proliferation of KCL-22 cells exposed to imatinib-dasatinib, imatinib-asciminib and dasatinib-asciminib combinations and calculated combination index graphs for each case. Moreover, using the median-effect equation we calculated how much axitinib can reduce the growth advantage of T315I mutant clones in combination with available drugs. In addition, we calculated how much the total drug burden could be reduced by combinations using asciminib and other drugs, and evaluated which mutations such combinations might be sensitive to. RESULTS: Asciminib had synergistic interactions with imatinib or dasatinib in KCL-22 cells at high degrees of inhibition. Interestingly, some antagonism between asciminib and the other drugs was present at lower degrees on inhibition. Simulations revealed that asciminib may allow for dose reductions, and its complementary resistance profile could reduce the risk of mutation based resistance. Axitinib, however, had only a minor effect on T315I growth advantage. CONCLUSIONS: Given how asciminib combinations were synergistic in vitro, our modelling suggests that drug combinations involving asciminib should allow for lower total drug doses, and may result in a reduced spectrum of observed resistance mutations. On the other hand, a combination involving axitinib was not shown to be useful in countering drug resistance.

Pattern and Dynamics of FLT3 Duplications.

FMS-like tyrosine kinase 3 (FLT3) is mutated in ∼30% of patients that suffer from acute myeloid leukemia (AML). In about 25% of all AML patients, in-frame insertions are observed in the sequence. Most of those insertions are internal tandem duplications (ITDs) of a sequence from the protein. The characteristics of such mutations in terms of length, sequence, and location were hitherto studied in different populations, but not in a comprehensive mutation database. Here, in-frame insertions into the FLT3 gene were extracted from the Catalogue of Somatic Mutations in Cancer (COSMIC) database. These were analyzed with respect to the length, location, and sequence of the mutations. Furthermore, characteristic strings (sequences) of different lengths were identified. Mutations were shown to occur most often in the juxtamembrane zipper (JM-Z) domain of FLT3, followed by the hinge domain and first tyrosine kinase domain (TKD1), upstream of the phosphate-binding loop (P-loop). Interestingly, there are specific hot spot residues where insertions are more likely to occur. The insertions vary in length between one and 67 amino acids, with the largest insertions spanning the phosphate binding loop. Insertions that occur downstream of the P-loop are shorter. Our analysis further shows that acidic and aromatic residues are enriched in the insertions. Finally, molecular dynamics simulations were run for FLT3 with ITD insertions in the hinge and tyrosine kinase domains. On the basis of the findings, a mechanism is proposed for activation by ITDs, according to which there is no direct coupling between the length of the insertion and the activity of the mutated protein. The effect of insertions on the sensitivity of FLT3 to kinase inhibitors is discussed based on our findings.

Stochastic modelling of tyrosine kinase inhibitor rotation therapy in chronic myeloid leukaemia.

BACKGROUND: Resistance towards targeted cancer treatments caused by single nucleotide variations is a major issue in many malignancies. Currently, there are a number of available drugs for chronic myeloid leukaemia (CML), which are overcome by different sets of mutations. The main aim of this study was to explore if it can be possible to exploit this and create a treatment protocol that outperforms each drug on its own. METHODS: We present a computer program to test different treatment protocols against CML, based on available resistance mutation growth data. The evolution of a relatively stable pool of cancer stem cells is modelled as a stochastic process, with the growth of cells expressing a tumourigenic protein (here, Abl1) and any emerging mutants determined principally by the drugs used in the therapy. RESULTS: There can be some benefit to Bosutinib-Ponatinib rotation therapy even if the mutation status is unknown, whereas Imatinib-Nilotinib rotation is unlikely to improve the outcomes. Furthermore, an interplay between growth inhibition and selection effects generates a non-linear relationship between drug doses and the risk of developing resistance. CONCLUSIONS: Drug rotation therapy might be able to delay the onset of resistance in CML patients without costly ongoing observation of mutation status. Moreover, the simulations give credence to the suggestion that lower drug concentrations may achieve better results following major molecular response in CML.

Specific Ion and Concentration Effects in Acetate Solutions with Na+ , K+ and Cs.

How salt ions affect solutes and the water beyond the solvation shell is not well understood. Molecular dynamics simulations of alkali-acetate solutions were analysed here in order to examine if, and how, different cations and solute concentrations affect the water structure and the interactions between water and acetates. The results revealed that water structure is perturbed to more than 1 nm away from the acetates and that this effect is more pronounced in physiological than in molar electrolyte concentrations. Analysis of simulations of a soluble protein revealed that the water orientation is perturbed to at least 1.5 nm from the protein structure. Furthermore, modifications to the orientation of water around carboxylate side chains were shown to depend on the local environment on the protein surface, and could extend to well over 1 nm, which may have an effect on protein dynamics during MD simulations in small water boxes.

Conformational modifications induced by internal tandem duplications on the FLT3 kinase and juxtamembrane domains.

The aberrant expression of FLT3 tyrosine kinase is associated primarily with acute myeloid leukaemia. This blood malignancy is often related to the onset of internal tandem duplications (ITDs) in the native sequence of the protein. The ITDs occur mainly in the juxtamembrane domain of the protein and alter the normal activity of the enzyme. In this work, we have studied the native form of FLT3 and six mutants by molecular dynamics simulations. The catalytic activity of FLT3 is exerted by the tyrosine kinase domain (KD) and regulated by the juxtamembrane (JM) domain. Analysis of the dynamics of these two domains have shown that the introduction of ITDs in the JM domain alters both structural and dynamic parameters. The presence of ITDs allowed the protein to span a larger portion of the conformational space, particularly in the JM domain and the activation loop. The FLT3 mutants were found to adopt more stable configurations than the native enzyme. This was due to the different arrangements assumed by the JM domain. Larger fluctuations of the activation loop were found in four of the six mutants. In the native FLT3, the key residue Tyr572 is involved in a strong and stable interaction with an ion pair. This interaction, which is thought to keep the JM in place hence regulating the activity of the enzyme, was found to break in all FLT3 mutants.

Computer simulations of the signalling network in FLT3 +-acute myeloid leukaemia - indications for an optimal dosage of inhibitors against FLT3 and CDK6.

BACKGROUND: Mutations in the FMS-like tyrosine kinase 3 (FLT3) are associated with uncontrolled cellular functions that contribute to the development of acute myeloid leukaemia (AML). We performed computer simulations of the FLT3-dependent signalling network in order to study the pathways that are involved in AML development and resistance to targeted therapies. RESULTS: Analysis of the simulations revealed the presence of alternative pathways through phosphoinositide 3 kinase (PI3K) and SH2-containing sequence proteins (SHC), that could overcome inhibition of FLT3. Inhibition of cyclin dependent kinase 6 (CDK6), a related molecular target, was also tested in the simulation but was not found to yield sufficient benefits alone. CONCLUSIONS: The PI3K pathway provided a basis for resistance to treatments. Alternative signalling pathways could not, however, restore cancer growth signals (proliferation and loss of apoptosis) to the same levels as prior to treatment, which may explain why FLT3 resistance mutations are the most common resistance mechanism. Finally, sensitivity analysis suggested the existence of optimal doses of FLT3 and CDK6 inhibitors in terms of efficacy and toxicity.