Capture and quality control mechanisms for adenosine-5'-triphosphate binding

Discovery of the Xenon-Protein Interactome Using Large-Scale Measurements of Protein Folding and Stability

The intermolecular interactions of noble gases in biological systems are associated with numerous biochemical responses, including apoptosis, inflammation, anesthesia, analgesia, and neuroprotection. The molecular modes of action underlying these responses are largely unknown. This is in large part due to the limited experimental techniques to study protein-gas interactions. The few techniques that are amenable to such studies are relatively low-throughput and require large amounts of purified proteins. Thus, they do not enable the large-scale analyses that are useful for protein target discovery. Here, we report the application of stability of proteins from rates of oxidation (SPROX) and limited proteolysis (LiP) methodologies to detect protein-xenon interactions on the proteomic scale using protein folding stability measurements. Over 5000 methionine-containing peptides and over 5000 semi-tryptic peptides, mapping to ∼1500 and ∼950 proteins, respectively, in the yeast proteome, were assayed for Xe-interacting activity using the SPROX and LiP techniques. The SPROX and LiP analyses identified 31 and 60 Xe-interacting proteins, respectively, none of which were previously known to bind Xe. A bioinformatics analysis of the proteomic results revealed that these Xe-interacting proteins were enriched in those involved in ATP-driven processes. A fraction of the protein targets that were identified are tied to previously established modes of action related to xenon's anesthetic and organoprotective properties. These results enrich our knowledge and understanding of biologically relevant xenon interactions. The sample preparation protocols and analytical methodologies developed here for xenon are also generally applicable to the discovery of a wide range of other protein-gas interactions in complex biological mixtures, such as cell lysates.

Assessing the performance of MM/PBSA and MM/GBSA methods. 7. Entropy effects on the performance of end-point binding free energy calculation approaches

Entropy effects play an important role in drug-target interactions, but the entropic contribution to ligand-binding affinity is often neglected by end-point binding free energy calculation methods, such as MM/GBSA and MM/PBSA, due to the expensive computational cost of normal mode analysis (NMA). Here, we systematically investigated entropy effects on the prediction power of MM/GBSA and MM/PBSA using >1500 protein-ligand systems and six representative AMBER force fields. Two computationally efficient methods, including NMA based on truncated structures and the interaction entropy approach, were used to estimate the entropic contributions to ligand-target binding free energies. In terms of the overall accuracy, we found that, for the minimized structures, in most cases the inclusion of the conformational entropies predicted by truncated NMA (enthalpynmode_min_9Å) compromises the overall accuracy of MM/GBSA and MM/PBSA compared with the enthalpies calculated based on the minimized structures (enthalpymin). However, for the MD trajectories, the binding free energies can be improved by the inclusion of the conformation entropies predicted by either truncated-NMA for a relatively high dielectric constant (εin = 4) or the interaction entropy method for εin = 1-4. In terms of reproducing the absolute binding free energies, the binding free energies estimated by including the truncated-NMA entropies based on the MD trajectories (ΔGnmode_md_9Å) give the lowest average absolute deviations against the experimental data among all the tested strategies for both MM/GBSA and MM/PBSA. Although the inclusion of the truncated NMA based on the MD trajectories (ΔGnmode_md_9Å) for a relatively high dielectric constant gave the overall best result and the lowest average absolute deviations against the experimental data (for the ff03 force field), it needs too much computational time. Alternatively, considering that the interaction entropy method does not incur any additional computational cost and can give comparable (at high dielectric constant, εin = 4) or even better (at low dielectric constant, εin = 1-2) results than the truncated-NMA entropy (ΔGnmode_md_9Å), the interaction entropy approach is recommended to estimate the entropic component for MM/GBSA and MM/PBSA based on MD trajectories, especially for a diverse dataset. Furthermore, we compared the predictions of MM/GBSA with six different AMBER force fields. The results show that the ff03 force field (ff03 for proteins and gaff with AM1-BCC charges for ligands) performs the best, but the predictions given by the tested force fields are comparable, implying that the MM/GBSA predictions are not very sensitive to force fields.

Characterizing Drug-Target Residence Time with Metadynamics: How To Achieve Dissociation Rate Efficiently without Losing Accuracy against Time-Consuming Approaches

Drug-target residence time plays a vital role in drug efficacy. However, there is still no effective strategy to predict drug residence time. Here, we propose to use the optimized (or minimized) structures derived from holo-state proteins to calculate drug residence time, which could give a comparable or even better prediction accuracy compared with those calculated utilizing a large number of molecular dynamics (MD) structures based on the Poisson process. Besides, in addition to the Poisson process, one may use fewer samples for predicting residence time due to the reason that, in a large extent, the calculated drug residence time is stable and independent of the number of samples used for the prediction. With remarkably reduced computational load, the proposed strategy may be promising for large-scale drug residence time prediction, such as post-processing in virtual screening (VS) and lead compound optimization.

The structural basis of ATP as an allosteric modulator

Adenosine-5'-triphosphate (ATP) is generally regarded as a substrate for energy currency and protein modification. Recent findings uncovered the allosteric function of ATP in cellular signal transduction but little is understood about this critical behavior of ATP. Through extensive analysis of ATP in solution and proteins, we found that the free ATP can exist in the compact and extended conformations in solution, and the two different conformational characteristics may be responsible for ATP to exert distinct biological functions: ATP molecules adopt both compact and extended conformations in the allosteric binding sites but conserve extended conformations in the substrate binding sites. Nudged elastic band simulations unveiled the distinct dynamic processes of ATP binding to the corresponding allosteric and substrate binding sites of uridine monophosphate kinase, and suggested that in solution ATP preferentially binds to the substrate binding sites of proteins. When the ATP molecules occupy the allosteric binding sites, the allosteric trigger from ATP to fuel allosteric communication between allosteric and functional sites is stemmed mainly from the triphosphate part of ATP, with a small number from the adenine part of ATP. Taken together, our results provide overall understanding of ATP allosteric functions responsible for regulation in biological systems.

P-loop conformation governed crizotinib resistance in G2032R-mutated ROS1 tyrosine kinase: clues from free energy landscape

Tyrosine kinases are regarded as excellent targets for chemical drug therapy of carcinomas. However, under strong purifying selection, drug resistance usually occurs in the cancer cells within a short term. Many cases of drug resistance have been found to be associated with secondary mutations in drug target, which lead to the attenuated drug-target interactions. For example, recently, an acquired secondary mutation, G2032R, has been detected in the drug target, ROS1 tyrosine kinase, from a crizotinib-resistant patient, who responded poorly to crizotinib within a very short therapeutic term. It was supposed that the mutation was located at the solvent front and might hinder the drug binding. However, a different fact could be uncovered by the simulations reported in this study. Here, free energy surfaces were characterized by the drug-target distance and the phosphate-binding loop (P-loop) conformational change of the crizotinib-ROS1 complex through advanced molecular dynamics techniques, and it was revealed that the more rigid P-loop region in the G2032R-mutated ROS1 was primarily responsible for the crizotinib resistance, which on one hand, impaired the binding of crizotinib directly, and on the other hand, shortened the residence time induced by the flattened free energy surface. Therefore, both of the binding affinity and the drug residence time should be emphasized in rational drug design to overcome the kinase resistance.

Cancers can eventually confer drug resistance to the continued medication. In most cases, mutations occurred in a drug target can attenuate the binding affinity of the drugs. Here, we studied the drug resistance mechanisms of the mutations G2032R in the ROS1 tyrosine kinase in fusion-type NSCLC. It is well known that the phosphate-binding loop (P-loop) plays a vital role in the binding of competitive inhibitors in tyrosine kinases, and numerous mutations have been found occurred around the P-loop, which may affect the binding/unbinding process of a drug. Free energy surfaces were constructed to characterize the impact of the mutation to the binding/unbinding process of a well-known NSCLC drug, crizotinib. Two advanced free energy calculation methods, namely funnel based well-tempered metadynamics and umbrella sampling based absolute binding free energy calculation achieved consistent results with the experimental data, suggesting that the rigid P-loop of the mutated target was mainly responsible for the crizotinib resistance to ROS1 tyrosine kinase.