Machine learning and protein allostery

β-sheets mediate the conformational change and allosteric signal transmission between the AsLOV2 termini

Avena sativa phototropin 1 light-oxygen-voltage 2 domain (AsLOV2) is a model protein of Per-Arnt-Sim (PAS) superfamily, characterized by conformational changes in response to external environmental stimuli. This conformational change begins with the unfolding of the N-terminal A'α helix in the dark state followed by the unfolding of the C-terminal Jα helix. The light state is characterized by the unfolded termini and the subsequent modifications in hydrogen bond patterns. In this photoreceptor, β-sheets are identified as crucial components for mediating allosteric signal transmission between the two termini. Through combined experimental and computational investigations, the Hβ and Iβ strands are recognized as the most critical and influential β-sheets in AsLOV2's allosteric mechanism. To elucidate the role of these β-sheets, we introduced 13 distinct mutations (F490L, N492A, L493A, F494L, H495L, L496F, Q497A, R500A, F509L, Q513A, L514A, D515V, and T517V) and conducted comprehensive molecular dynamics simulations. In-depth hydrogen bond analyses emphasized the role of two hydrogen bonds, Asn482-Leu453 and Gln479-Val520, in the observed distinct behaviors of L493A, L496F, Q497A, and D515V mutants. This illustrates the role of β-sheets in the transmission of the allosteric signal upon the photoactivation of the light state.

Microsecond Molecular Dynamics Simulations and Markov State Models of Mutation-Induced Allosteric Mechanisms for the Light-Oxygen-Voltage 2 Protein : Revealing Structural Basis of Signal Transmission Induced by Photoactivation of the Light Protein State

Avena Sativa phototropin 1 Light-oxygen-voltage 2 domain (AsLOV2) is the model protein of Per-Arnt-Sim (PAS) superfamily, characterized by conformational changes in response to external environmental stimuli. This conformational change is initiated by the unfolding of the N-terminal helix in the dark state followed by the unfolding of the C-terminal helix. The light state is defined by the unfolded termini and the subsequent modifications in hydrogen bond patterns. In this photoreceptor, β-sheets have been identified as crucial components for mediating allosteric signal transmission between the two termini. In this study, we combined microsecond all-atm molecular dynamics simulations and Markov state modeling of conformational states to quantify molecular basis of mutation-induced allostery in the AsLOV2 protein. Through a combination of computational investigations, we determine that the Hβ and Iβ strands are the most critical structural elements involved in the allosteric mechanism. To elucidate the role of these β-sheets, we introduced 13 distinct mutations (F490L, N492A, L493A, F494L, H495L, L496F, Q497A, R500A, F509L, Q513A, L514A, D515V, and T517V) and conducted comprehensive simulation analysis. The results highlighted the role of two hydrogen bond Asn482-Leu453 and Gln479-Val520 in the observed distinct behaviors of L493A, L496F, Q497A, and D515V mutants. The comprehensive atomistic-level analysis of the conformational landscapes revealed the critical functional role of β-sheet segments in the transmission of the allosteric signal upon the photoactivation of the light state.

Advances of Predicting Allosteric Mechanisms Through Protein Contact in New Technologies and Their Application

Allostery is an intriguing phenomenon wherein the binding activity of a biological macromolecule is modulated via non-canonical binding site, resulting in synchronized functional changes. The mechanics underlying allostery are relatively complex and this review is focused on common methodologies used to study allostery, such as X-ray crystallography, NMR spectroscopy, and HDXMS. Different methodological approaches are used to generate data in different scenarios. For example, X-ray crystallography provides high-resolution structural information, NMR spectroscopy offers dynamic insights into allosteric interactions in solution, and HDXMS provides information on protein dynamics. The residue transition state (RTS) approach has emerged as a critical tool in understanding the energetics and conformational changes associated with allosteric regulation. Allostery has significant implications in drug discovery, gene transcription, disease diagnosis, and enzyme catalysis. Enzymes' catalytic activity can be modulated by allosteric regulation, offering opportunities to develop novel therapeutic alternatives. Understanding allosteric mechanisms associated with infectious organisms like SARS-CoV and bacterial pathogens can aid in the development of new antiviral drugs and antibiotics. Allosteric mechanisms are crucial in the regulation of a variety of signal transduction and cell metabolism pathways, which in turn govern various cellular processes. Despite progress, challenges remain in identifying allosteric sites and characterizing their contribution to a variety of biological processes. Increased understanding of these mechanisms can help develop allosteric systems specifically designed to modulate key biological mechanisms, providing novel opportunities for the development of targeted therapeutics. Therefore, the current review aims to summarize common methodologies that are used to further our understanding of allosteric mechanisms. In conclusion, this review provides insights into the methodologies used for the study of allostery, its applications in in silico modeling, the mechanisms underlying antibody allostery, and the ongoing challenges and prospects in advancing our comprehension of this intriguing phenomenon.

PASSerRank: Prediction of allosteric sites with learning to rank

Allostery plays a crucial role in regulating protein activity, making it a highly sought-after target in drug development. One of the major challenges in allosteric drug research is the identification of allosteric sites. In recent years, many computational models have been developed for accurate allosteric site prediction. Most of these models focus on designing a general rule that can be applied to pockets of proteins from various families. In this study, we present a new approach using the concept of Learning to Rank (LTR). The LTR model ranks pockets based on their relevance to allosteric sites, that is, how well a pocket meets the characteristics of known allosteric sites. After the training and validation on two datasets, the Allosteric Database (ASD) and CASBench, the LTR model was able to rank an allosteric pocket in the top three positions for 83.6% and 80.5% of test proteins, respectively. The model outperforms other common machine learning models with higher F1 scores (0.662 in ASD and 0.608 in CASBench) and Matthews correlation coefficients (0.645 in ASD and 0.589 in CASBench). The trained model is available on the PASSer platform (https://passer.smu.edu) to aid in drug discovery research.

Markov State Models and Perturbation-Based Approaches Reveal Distinct Dynamic Signatures and Hidden Allosteric Pockets in the Emerging SARS-Cov-2 Spike Omicron Variant Complexes with the Host Receptor: The Interplay of Dynamics and Convergent Evolution Modulates Allostery and Functional Mechanisms

The new generation of SARS-CoV-2 Omicron variants displayed a significant growth advantage and increased viral fitness by acquiring convergent mutations, suggesting that the immune pressure can promote convergent evolution leading to the sudden acceleration of SARS-CoV-2 evolution. In the current study, we combined structural modeling, microsecond molecular dynamics simulations, and Markov state models to characterize conformational landscapes and identify specific dynamic signatures of the SARS-CoV-2 spike complexes with the host receptor ACE2 for the recently emerged highly transmissible XBB.1, XBB.1.5, BQ.1, and BQ.1.1 Omicron variants. Microsecond simulations and Markovian modeling provided a detailed characterization of the functional conformational states and revealed the increased thermodynamic stabilization of the XBB.1.5 subvariant, which can be contrasted to more dynamic BQ.1 and BQ.1.1 subvariants. Despite considerable structural similarities, Omicron mutations can induce unique dynamic signatures and specific distributions of the conformational states. The results suggested that variant-specific changes of the conformational mobility in the functional interfacial loops of the receptor-binding domain in the SARS-CoV-2 spike protein can be fine-tuned through crosstalk between convergent mutations which could provide an evolutionary path for modulation of immune escape. By combining atomistic simulations and Markovian modeling analysis with perturbation-based approaches, we determined important complementary roles of convergent mutation sites as effectors and receivers of allosteric signaling involved in modulation of conformational plasticity and regulation of allosteric communications. This study also revealed hidden allosteric pockets and suggested that convergent mutation sites could control evolution and distribution of allosteric pockets through modulation of conformational plasticity in the flexible adaptable regions.