Energy landscape quantifications of histone H3.3 recognition by chaperone DAXX reveal an uncoupled binding specificity and affinity

Histone variant H3.3 differs from the canonical histone H3.1 by only five amino acids, yet its chaperone death domain-associated protein (DAXX) can specifically recognize H3.3 over H3.1, despite having a large DAXX-interacting surface on the H3.3-H4 heterodimer common to that on the H3.1-H4 complex. This observation gives rise to the question of, from the binding energy point view, how high binding specificity may be achieved with small differences of the overall binding energy for protein-protein interactions in general. Here we investigate the mechanism of coupling of binding specificity and affinity in protein-protein interactions using the DAXX-H3.3-H4 complex as a model. Using a multi-scale method, we found that the hydrophobic interactions between DAXX and the H3.3-specific region contributed to their initial binding process. And the structural flexibility of the interacting partners contributed to the binding affinity after their encounter. By quantifying the free energy landscape, we revealed that the interaction between the specific residues of H3.3 and DAXX decreased the encounter barrier height while the folding of H3.3-H4 and DAXX increased the depth of the free energy basin of the final binding state. The encounter barrier height, which is not coupled to the thermodynamic stability of the final binding state, had a marked effect on the initial binding rate of flexible histones and chaperones. Based on the energy landscape theory, we found that the intrinsic binding energy funnel of this uncoupled recognition process was affected by the structural flexibility and the flexibility modulated the degree of coupling between binding specificity and affinity. Our work offers a biophysical explanation of the specific recognition between the histones and their chaperones, and also extends the use of energy landscape theory for understanding molecular recognitions in general.

Comparison of ATP-binding pockets and discovery of homologous recombination inhibitors

The ATP binding sites of many enzymes are structurally related, which complicates their development as therapeutic targets. In this work, we explore a diverse set of ATPases and compare their ATP binding pockets using different strategies, including direct and indirect structural methods, in search of pockets attractive for drug discovery. We pursue different direct and indirect structural strategies, as well as ligandability assessments to help guide target selection. The analyses indicate human RAD51, an enzyme crucial in homologous recombination, as a promising, tractable target. Inhibition of RAD51 has shown promise in the treatment of certain cancers but more potent inhibitors are needed. Thus, we design compounds computationally against the ATP binding pocket of RAD51 with consideration of multiple criteria, including predicted specificity, drug-likeness, and toxicity. The molecules designed are evaluated experimentally using molecular and cell-based assays. Our results provide two novel hit compounds against RAD51 and illustrate a computational pipeline to design new inhibitors against ATPases.

Perspectives on the landscape and flux theory for describing emergent behaviors of the biological systems

We give a review on the landscape theory of the equilibrium biological systems and landscape-flux theory of the nonequilibrium biological systems as the global driving force. The emergences of the behaviors, the associated thermodynamics in terms of the entropy and free energy and dynamics in terms of the rate and paths have been quantitatively demonstrated. The hierarchical organization structures have been discussed. The biological applications ranging from protein folding, biomolecular recognition, specificity, biomolecular evolution and design for equilibrium systems as well as cell cycle, differentiation and development, cancer, neural networks and brain function, and evolution for nonequilibrium systems, cross-scale studies of genome structural dynamics and experimental quantifications/verifications of the landscape and flux are illustrated. Together, this gives an overall global physical and quantitative picture in terms of the landscape and flux for the behaviors, dynamics and functions of biological systems.

Physics of biomolecular recognition and conformational dynamics

Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.

Rapid Assessment of T-Cell Receptor Specificity of the Immune Repertoire

Accurate assessment of TCR-antigen specificity at the whole immune repertoire level lies at the heart of improved cancer immunotherapy, but predictive models capable of high-throughput assessment of TCR-peptide pairs are lacking. Recent advances in deep sequencing and crystallography have enriched the data available for studying TCR-p-MHC systems. Here, we introduce a pairwise energy model, RACER, for rapid assessment of TCR-peptide affinity at the immune repertoire level. RACER applies supervised machine learning to efficiently and accurately resolve strong TCR-peptide binding pairs from weak ones. The trained parameters further enable a physical interpretation of interacting patterns encoded in each specific TCR-p-MHC system. When applied to simulate thymic selection of an MHC-restricted T-cell repertoire, RACER accurately estimates recognition rates for tumor-associated neoantigens and foreign peptides, thus demonstrating its utility in helping address the large computational challenge of reliably identifying the properties of tumor antigen-specific T-cells at the level of an individual patient's immune repertoire.

Nonequilibrium Thermodynamics in Cell Biology: Extending Equilibrium Formalism to Cover Living Systems

We discuss new developments in the nonequilibrium dynamics and thermodynamics of living systems, giving a few examples to demonstrate the importance of nonequilibrium thermodynamics for understanding biological dynamics and functions. We study single-molecule enzyme dynamics, in which the nonequilibrium thermodynamic and dynamic driving forces of chemical potential and flux are crucial for the emergence of non-Michaelis-Menten kinetics. We explore single-gene expression dynamics, in which nonequilibrium dissipation can suppress fluctuations. We investigate the cell cycle and identify the nutrition supply as the energy input that sustains the stability, speed, and coherence of cell cycle oscillation, from which the different vital phases of the cell cycle emerge. We examine neural decision-making processes and find the trade-offs among speed, accuracy, and thermodynamic costs that are important for neural function. Lastly, we consider the thermodynamic cost for specificity in cellular signaling and adaptation.

Position-, disorder-, and salt-dependent diffusion in binding-coupled-folding of intrinsically disordered proteins

Successful extensions of protein-folding energy landscape theory to intrinsically disordered proteins' (IDPs') binding-coupled-folding transition can enormously simplify this biomolecular process into diffusion along a limited number of reaction coordinates, and the dynamics subsequently is described by Kramers' rate theory. As the critical pre-factor, the diffusion coefficient D has direct implications on the binding kinetics. Here, we employ a structure-based model (SBM) to calculate D in the binding-folding of an IDP prototype. We identify a strong position-dependent D during binding by applying a reaction coordinate that directly measures the fluctuations in a Cartesian configuration space. Using the malleability of the SBM, we modulate the degree of conformational disorder in an isolated IDP and determine complex effects of intrinsic disorder on D varying for different binding stages. Here, D tends to increase with disorder during initial binding but shows a non-monotonic relationship with disorder in terms of a decrease-followed-by-increase in D during the late binding stage. The salt concentration, which correlates with electrostatic interactions via Debye-Hückel theory in our SBM, also modulates D in a stepwise way. The speeding up of diffusion by electrostatic interactions is observed during the formation of the encounter complex at the beginning of binding, while the last diffusive binding dynamics is hindered by non-native salt bridges. Because D describes the diffusive speed locally, which implicitly reflects the roughness of the energy landscape, we are eventually able to portray the binding energy landscape, including that from IDPs' binding, then to binding with partial folding, and finally to rigid docking, as well as that under different environmental salt concentrations. Our theoretical results provide key mechanistic insights into IDPs' binding-folding, which is internally conformation- and externally salt-controlled with respect to diffusion.

Quantifying the Kinetic Residence Time as a Potential Complement to Affinity for the Aptamer Selection

Aptamers have been widely developed for biotechnological and therapeutic applications in recent years. Increasing evidence shows that the quality of the aptamer is determined not only by the thermodynamic stability but also by the kinetic residence time when binding with its target. However, both experimental and computational selection methods of aptamers concentrate solely on the binding affinity optimization. Here, we propose a computational method for the quantification of the residence time by describing the kinetics on the underlying funneled binding energy landscape of aptamer-target complex as a diffusion process. The quantified residence time is examined to have the capacity to discriminate native aptamer-target complexes against non-native ones. It is found that the residence time is correlated with the binding affinity but with significant dispersion, suggesting that the residence time can be a potential complement in selecting the aptamer. On the basis of the results, a two-dimensional selection method with both the thermodynamic binding affinity and the residence time as kinetic specificity is suggested for diverse applications of aptamers. This alters conventional aptamer selection methods by considering the kinetic residence time in addition to the affinity.

Quantifying the Intrinsic Conformation Energy Landscape Topography of Proteins with Large-Scale Open-Closed Transition

Large-scale conformational changes of proteins, including the open-closed transitions, are crucial for a variety of protein functions. These open-closed transitions are often associated with ligand binding. However, the understandings of the underlying mechanisms of the conformational changes within proteins during the open-closed transitions are still challenging at present. In this study, we quantified the intrinsic underlying conformational energy landscapes of five different proteins with large-scale open-closed transitions. This is realized by exploring the underlying density of states and the intrinsic conformational energy landscape topography measure Λ. Λ is a dimensionless ratio of conformational energy gap δE versus conformational energy roughness δE and configurational entropy S or size of the intrinsic conformational energy landscape. By quantifying the Λ of intrinsic open-closed conformational (Λ^oc) and intrinsic global folding (Λ^global) energy landscapes, we show that both intrinsic open-closed conformation energy and entropy landscapes are funneled toward the closed state. Furthermore, our results indicate the strong correlations between Λ and thermodynamics (conformational state transition temperature against trapping temperature) as well as between Λ and kinetics (open-closed kinetic time) of these proteins. This shows that the intrinsic conformational landscape topography determines both the conformational thermodynamic stability and kinetic speed of the conformational dynamics. Our investigations provide important insights for understanding the fundamental mechanisms of the protein conformational dynamics in a physical and global way.

SPA-LN: a scoring function of ligand-nucleic acid interactions via optimizing both specificity and affinity

Nucleic acids have been widely recognized as potential targets in drug discovery and aptamer selection. Quantifying the interactions between small molecules and nucleic acids is critical to discover lead compounds and design novel aptamers. Scoring function is normally employed to quantify the interactions in structure-based virtual screening. However, the predictive power of nucleic acid–ligand scoring functions is still a challenge compared to other types of biomolecular recognition. With the rapid growth of experimentally determined nucleic acid–ligand complex structures, in this work, we develop a knowledge-based scoring function of nucleic acid–ligand interactions, namely SPA-LN. SPA-LN is optimized by maximizing both the affinity and specificity of native complex structures. The development strategy is different from those of previous nucleic acid–ligand scoring functions which focus on the affinity only in the optimization. The native conformation is stabilized while non-native conformations are destabilized by our optimization, making the funnel-like binding energy landscape more biased toward the native state. The performance of SPA-LN validates the development strategy and provides a relatively more accurate way to score the nucleic acid–ligand interactions.

Scoring Functions of Protein-Ligand Interactions

Scoring function of protein-ligand interactions is used to recognize the “native” binding pose of a ligand on the protein and to predict the binding affinity, so that the active small molecules can be discriminated from the non-active ones. Scoring function is widely used in computationally molecular docking and structure-based drug discovery. The development and improvement of scoring functions have broad implications in pharmaceutical industry and academic research. During the past three decades, much progress have been made in methodology and accuracy for scoring functions, and many successful cases have be witnessed in virtual database screening. In this chapter, the authors introduced the basic types of scoring functions and their derivations, the commonly-used evaluation methods and benchmarks, as well as the underlying challenges and current solutions. Finally, the authors discussed the promising directions to improve and develop scoring functions for future molecular docking-based drug discovery.

Energy Landscape Topography Reveals the Underlying Link Between Binding Specificity and Activity of Enzymes

Enzyme activity (often quantified by k_cat/K_m) is the main function of enzyme when it is active against the specific substrate. Higher or lower activities are highly desired for the design of novel enzyme and drug resistance. However, it is difficult to measure the activities of all possible variants and find the “hot-spot” within the limit of experimental time. In this study, we explore the underlying energy landscape of enzyme-substrate interactions and introduce the intrinsic specificity ratio (ISR), which reflects the landscape topography. By studying two concrete systems, we uncover the statistical correlation between the intrinsic specificity and the enzyme activity k_cat/K_m. This physics-based concept and method show that the energy landscape topography is valuable for understanding the relationship between enzyme specificity and activity. In addition, it can reveal the underlying mechanism of enzyme-substrate actions and has potential applications on enzyme design.

Universal statistical fluctuations in thermodynamics and kinetics of single molecular recognition

We investigated the main universal statistical distributions of single molecular recognition. The distributions of the single molecule binding free energy spectrum or density of states were characterized in the ligand-receptor binding energy landscape. The analytical results are consistent with the microscopic molecular simulations. The free energy distribution of different binding modes or states for a single molecule ligand receptor pair is approximately Gaussian near the mean and exponential at the tail. The equilibrium constant of single molecule binding is log-normal distributed near the mean and power law distributed near the tail. Additionally, we found that the kinetics distribution of single molecule ligand binding can be characterized by log-normal around the mean and power law distribution near the tail. This distribution is caused by exploration of the underlying inhomogeneous free energy landscape. Different ligand-receptor binding complexes have the same universal form of distribution but differ in parameters.

Optimizing the affinity and specificity of ligand binding with the inclusion of solvation effect

Solvation effect is an important factor for protein-ligand binding in aqueous water. Previous scoring function of protein-ligand interactions rarely incorporates the solvation model into the quantification of protein-ligand interactions, mainly due to the immense computational cost, especially in the structure-based virtual screening, and nontransferable application of independently optimized atomic solvation parameters. In order to overcome these barriers, we effectively combine knowledge-based atom-pair potentials and the atomic solvation energy of charge-independent implicit solvent model in the optimization of binding affinity and specificity. The resulting scoring functions with optimized atomic solvation parameters is named as specificity and affinity with solvation effect (SPA-SE). The performance of SPA-SE is evaluated and compared to 20 other scoring functions, as well as SPA. The comparative results show that SPA-SE outperforms all other scoring functions in binding affinity prediction and "native" pose identification. Our optimization validates that solvation effect is an important regulator to the stability and specificity of protein-ligand binding. The development strategy of SPA-SE sets an example for other scoring function to account for the solvation effect in biomolecular recognitions.

The universal statistical distributions of the affinity, equilibrium constants, kinetics and specificity in biomolecular recognition

We uncovered the universal statistical laws for the biomolecular recognition/binding process. We quantified the statistical energy landscapes for binding, from which we can characterize the distributions of the binding free energy (affinity), the equilibrium constants, the kinetics and the specificity by exploring the different ligands binding with a particular receptor. The results of the analytical studies are confirmed by the microscopic flexible docking simulations. The distribution of binding affinity is Gaussian around the mean and becomes exponential near the tail. The equilibrium constants of the binding follow a log-normal distribution around the mean and a power law distribution in the tail. The intrinsic specificity for biomolecular recognition measures the degree of discrimination of native versus non-native binding and the optimization of which becomes the maximization of the ratio of the free energy gap between the native state and the average of non-native states versus the roughness measured by the variance of the free energy landscape around its mean. The intrinsic specificity obeys a Gaussian distribution near the mean and an exponential distribution near the tail. Furthermore, the kinetics of binding follows a log-normal distribution near the mean and a power law distribution at the tail. Our study provides new insights into the statistical nature of thermodynamics, kinetics and function from different ligands binding with a specific receptor or equivalently specific ligand binding with different receptors. The elucidation of distributions of the kinetics and free energy has guiding roles in studying biomolecular recognition and function through small-molecule evolution and chemical genetics.

Uncovering the principles and underlying mechanisms of biomolecular recognition and molecular binding process is crucial for understanding the function and evolution, yet challenging. We meet the challenge by quantifying the statistical natures of the relevant physical variables of biomolecular recognition using the analytical model combined with microscopic flexible docking simulation methods. We uncovered the universal statistical laws obeyed by the affinity, equilibrium constant, intrinsic specificity and kinetics for biomolecular recognition. The general statistical laws based on energy landscape theory can serve as a conceptual framework for molecular recognition in biological repertoires. They can be applied to molecular selection, in vitro evolution process, high throughput screening and virtual screening for drug discovery. The statistical laws in combinations with experiments provide quantitative signatures of a specific ligand binding to a specific receptor, these resultant laws as a guideline will contribute to drug design against a specific target. Our developed statistical methodology is general and applicable for all other biomolecular recognitions.

Computational discovery and experimental verification of tyrosine kinase inhibitor pazopanib for the reversal of memory and cognitive deficits in rat model neurodegeneration

Pazopanib, a tyrosine kinase inhibitor marketed for cancer treatment, abrogates the course of neurodegeneration.

Cognition and memory impairment are hallmarks of the pathological cascade of various neurodegenerative disorders. Herein, we developed a novel computational strategy with two-dimensional virtual screening for not only affinity but also specificity. We integrated the two-dimensional virtual screening with ligand screening for 3D shape, electrostatic similarity and local binding site similarity to find existing drugs that may reduce the signs of cognitive deficits. For the first time, we found that pazopanib, a tyrosine kinase inhibitor marketed for cancer treatment, inhibits acetylcholinesterase (AchE) activities at sub-micromolar concentration. We evaluated and compared the effects of intragastrically-administered pazopanib with donepezil, a marketed AchE inhibitor, in cognitive and behavioral assays including the novel object recognition test, Y maze and Morris water maze test. Surprisingly, we found that pazopanib can restore memory loss and cognitive dysfunction to a similar extent as donepezil in a dosage of 15 mg kg^–1, only one fifth of the equivalent clinical dosage for cancer treatment. Furthermore, we demonstrated that pazopanib dramatically enhances the hippocampal Ach levels and increases the expression of the synaptic marker SYP. These findings suggest that pazopanib may become a viable treatment option for memory and cognitive deficits with a good safety profile in humans.

Specificity and affinity quantification of flexible recognition from underlying energy landscape topography

Flexibility in biomolecular recognition is essential and critical for many cellular activities. Flexible recognition often leads to moderate affinity but high specificity, in contradiction with the conventional wisdom that high affinity and high specificity are coupled. Furthermore, quantitative understanding of the role of flexibility in biomolecular recognition is still challenging. Here, we meet the challenge by quantifying the intrinsic biomolecular recognition energy landscapes with and without flexibility through the underlying density of states. We quantified the thermodynamic intrinsic specificity by the topography of the intrinsic binding energy landscape and the kinetic specificity by association rate. We found that the thermodynamic and kinetic specificity are strongly correlated. Furthermore, we found that flexibility decreases binding affinity on one hand, but increases binding specificity on the other hand, and the decreasing or increasing proportion of affinity and specificity are strongly correlated with the degree of flexibility. This shows more (less) flexibility leads to weaker (stronger) coupling between affinity and specificity. Our work provides a theoretical foundation and quantitative explanation of the previous qualitative studies on the relationship among flexibility, affinity and specificity. In addition, we found that the folding energy landscapes are more funneled with binding, indicating that binding helps folding during the recognition. Finally, we demonstrated that the whole binding-folding energy landscapes can be integrated by the rigid binding and isolated folding energy landscapes under weak flexibility. Our results provide a novel way to quantify the affinity and specificity in flexible biomolecular recognition.

Flexibility in biomolecular recognition is crucial for the function. Flexibility often leads to moderate binding affinity but high binding specificity, challenging the conventional wisdom that high specificity is guaranteed by high affinity. Currently, understanding of the relationship between affinity and specificity in flexible biomolecular recognition is still obscure, even in a qualitative way. By exploring the intrinsic biomolecular recognition energy landscapes, we provided a novel way to quantify the thermodynamic intrinsic specificity by energy landscape topography and kinetic specificity by association rate. We show quantitatively that flexibility decreases binding affinity while increases binding specificity, and the relative changes in affinity and specificity are strongly correlated with the degree of flexibility. Our results show that more (less) flexibility leads to weaker (stronger) coupling between affinity and specificity. Importantly, we demonstrated that flexibility modulates affinity and specificity through the underlying energy landscape. Our study establishes the quantitative relationship among flexibility, affinity and specificity, bridging the gap between theory and experiments.

Optimizing scoring function of protein-nucleic acid interactions with both affinity and specificity

Protein-nucleic acid (protein-DNA and protein-RNA) recognition is fundamental to the regulation of gene expression. Determination of the structures of the protein-nucleic acid recognition and insight into their interactions at molecular level are vital to understanding the regulation function. Recently, quantitative computational approach has been becoming an alternative of experimental technique for predicting the structures and interactions of biomolecular recognition. However, the progress of protein-nucleic acid structure prediction, especially protein-RNA, is far behind that of the protein-ligand and protein-protein structure predictions due to the lack of reliable and accurate scoring function for quantifying the protein-nucleic acid interactions. In this work, we developed an accurate scoring function (named as SPA-PN, SPecificity and Affinity of the Protein-Nucleic acid interactions) for protein-nucleic acid interactions by incorporating both the specificity and affinity into the optimization strategy. Specificity and affinity are two requirements of highly efficient and specific biomolecular recognition. Previous quantitative descriptions of the biomolecular interactions considered the affinity, but often ignored the specificity owing to the challenge of specificity quantification. We applied our concept of intrinsic specificity to connect the conventional specificity, which circumvents the challenge of specificity quantification. In addition to the affinity optimization, we incorporated the quantified intrinsic specificity into the optimization strategy of SPA-PN. The testing results and comparisons with other scoring functions validated that SPA-PN performs well on both the prediction of binding affinity and identification of native conformation. In terms of its performance, SPA-PN can be widely used to predict the protein-nucleic acid structures and quantify their interactions.