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.

Introducing intrinsic disorder reduces electrostatic steering in protein-protein interactions

Protein-protein interactions underlie many critical biology functions, such as cellular signaling and gene expression, in which electrostatic interactions can play a critical role in mediating the specificity and stability of protein complexes. A substantial portion of proteins are intrinsically disordered, and the influences of structural disorder on binding kinetics and thermodynamics have been widely investigated. However, whether the effect of electrostatic steering depends on structural disorder remains unexplored. In this work, we addressed the consequence of introducing intrinsic disorder in the electrostatic steering of the E3/Im3 complex using molecular dynamics simulation. Our results recapitulated the experimental observations that the responses of stability and kinetics to salt concentration for the ordered E3/Im3 complex were larger than those for the disordered E3/Im3 complex. Mechanistic analysis revealed that the native contact interactions involved in the encounter state and the transition state were essentially identical for both ordered and disordered E3. Therefore, the observed difference in electrostatic steering between ordered E3 and disordered E3 may result from their difference in conformation rather than their difference in binding mechanism. Because charged residues are frequently involved in protein-protein interactions, our results suggest that increasing structural disorder is expected to generally modulate the effect of electrostatic steering.

Exploration of nucleoprotein α-MoRE and XD interactions of Nipah and Hendra viruses

Henipavirus, including Hendra virus (HeV) and Nipah virus (NiV), is a newly discovered human pathogen genus. The nucleoprotein of Henipavirus contains an α-helical molecular recognition element (α-MoRE) that folds upon binding to the X domain (XD) of the phosphoprotein (P). In order to explore the conformational dynamics of free α-MoREs and the underlying binding-folding mechanism with XD, atomic force field-based and hybrid structure-based MD simulations were carried out. In our empirical force field-based simulations, characteristic structures and helicities of α-MoREs reveal the co-existence of partially structured and disordered conformations, as in the case of the well characterized cognate measles virus (MeV) α-MoRE. In spite of their overall similarity, the two α-MoREs display subtle helicity differences in their C-terminal region, but much different from that of MeV. For the α-MoRE/XD complexes, the results of our hybrid structure-based simulations provide the coupled binding-folding landscapes, and unveil a wide conformational selection mechanism at early binding stages, followed by a final induce-fit mechanism selection process. However, the HeV and NiV complexes have a lower binding barrier compared to that of MeV. Moreover, the HeV α-MoRE/XD complex shows much less coupling effects between binding and folding compared to that from both NiV and MeV. Our analysis revealed that contrary to NiV and MeV, the N- and C-terminal regions of the HeV α-MoRE maintains a low helicity also in the bound form.