Neural relational inference to learn long-range allosteric interactions in proteins from molecular dynamics simulations

Exploration of novel human neutrophil elastase inhibitors from natural compounds: Virtual screening, in vitro, molecular dynamics simulation and in vivo study

BACKGROUND

Human neutrophil elastase (HNE) is an important contributor to lung diseases such as acute lung injury (ALI) or acute respiratory distress syndrome. Therefore, this study aimed to identify natural HNE inhibitors with anti-inflammatory activity through machine learning algorithms, in vitro assays, molecular dynamic simulation, and an in vivo ALI assay.

METHODS

Based on the optimized Discovery Studio two-dimensional molecular descriptors, combined with different molecular fingerprints, six machine learning models were established using the Naïve Bayesian (NB) method to identify HNE inhibitors. Subsequently, the optimal model was utilized to screen 6925 drug-like compounds obtained from the Traditional Chinese Medicine Systems Pharmacy Database and Analysis Platform (TCMSP), followed by ADMET analysis. Finally, 10 compounds with reported anti-inflammatory activity were selected to determine their inhibitory activities against HNE in vitro, and the compounds with the best activity were selected for a 100 ns molecular dynamics simulation and its anti-inflammatory effect was evaluated using Poly (I:C)-induced ALI model.

RESULTS

The evaluation of the in vitro HNE inhibition efficiency of the 10 selected compounds showed that the flavonoid tricetin had the strongest inhibitory effect on HNE. The molecular dynamics simulation indicated that the binding of tricetin to HNE was relatively stable throughout the simulation. Importantly, in vivo experiments indicated that tricetin treatment substantially improved the Poly (I:C)-induced ALI.

CONCLUSION

The proposed NB model was proved valuable for exploring novel HNE inhibitors, and natural tricetin was screened out as a novel HNE inhibitor, which was confirmed by in vitro and in vivo assays for its inhibitory activities.

Integrating Molecular Dynamics and Machine Learning Algorithms to Predict the Functional Profile of Kinase Ligands

The modulation of protein function via designed small molecules is providing new opportunities in chemical biology and medicinal chemistry. While drugs have traditionally been developed to block enzymatic activities through active site occupation, a growing number of strategies now aim to control protein functions in an allosteric fashion, allowing for the tuning of a target's activation or deactivation via the modulation of the populations of conformational ensembles that underlie its function. In the context of the discovery of new active leads, it would be very useful to generate hypotheses for the functional impact of new ligands. Since the discovery and design of allosteric modulators (inhibitors/activators) is still a challenging and often serendipitous target, the development of a rapid and robust approach to predict the functional profile of a new ligand would significantly speed up candidate selection. Herein, we present different machine learning (ML) classifiers to distinguish between potential orthosteric and allosteric binders. Our approach integrates information on the chemical fingerprints of the ligands with descriptors that recapitulate ligand effects on protein functional motions. The latter are derived from molecular dynamics (MD) simulations of the target protein in complex with orthosteric or allosteric ligands. In this framework, we train and test different ML architectures, which are initially probed on the classification of orthosteric versus allosteric ligands for cyclin-dependent kinases (CDKs). The results demonstrate that different ML methods can successfully partition allosteric versus orthosteric effectors (although to different degrees). Next, we further test the models with FDA-approved CDK drugs, not included in the original dataset, as well as ligands that target other kinases, to test the range of applicability of these models outside of the domain on which they were developed. Overall, the results show that enriching the training dataset with chemical physics-based information on the protein-ligand dynamic cross-talk can significantly expand the reach and applicability of approaches for the prediction and classification of the mode of action of small molecules.

NRIMD, a Web Server for Analyzing Protein Allosteric Interactions Based on Molecular Dynamics Simulation

Long-range allosteric communication between distant sites and active sites in proteins is central to biological regulation but still poorly characterized, limiting the development of protein engineering and drug design. Addressing this gap, NRIMD is an open-access web server for analyzing long-range interactions in proteins from molecular dynamics (MD) simulations, such as the effect of mutations at distal sites or allosteric ligand binding at allosteric sites on the active center. Based on our recent works on neural relational inference using graph neural networks, this cloud-based web server accepts MD simulation data on any length of residues in the alpha-carbon skeleton format from mainstream MD software. The input trajectory data are validated at the frontend deployed on the cloud and then processed on the backend deployed on a high-performance computer system with a collection of complementary tools. The web server provides a one-stop-shop MD analysis platform to predict long-range interactions and their paths between distant sites and active sites. It provides a user-friendly interface for detailed analysis and visualization. To the best of our knowledge, NRIMD is the first-of-its-kind online service to provide comprehensive long-range interaction analysis on MD simulations, which significantly lowers the barrier of predictions on protein long-range interactions using deep learning. The NRIMD web server is publicly available at https://nrimd.luddy.indianapolis.iu.edu/.

The structural biology and dynamics of malate dehydrogenases

Malate dehydrogenase (MDH) enzymes catalyze the reversible oxidoreduction of malate to oxaloacetate using NAD(P) as a cofactor. This reaction is vital for metabolism and the exchange of reducing equivalents between cellular compartments. There are more than 100 structures of MDH in the Protein Data Bank, representing species from archaea, bacteria, and eukaryotes. This conserved family of enzymes shares a common nucleotide-binding domain, substrate-binding domain, and subunits associate to form a dimeric or a tetrameric enzyme. Despite the variety of crystallization conditions and ligands in the experimental structures, the conformation and configuration of MDH are similar. The quaternary structure and active site dynamics account for most conformational differences in the experimental MDH structures. Oligomerization appears essential for activity despite each subunit having a structurally independent active site. There are two dynamic regions within the active site that influence substrate binding and possibly catalysis, with one of these regions adjoining the subunit interface. In this review, we introduce the reader to the general structural framework of MDH highlighting the conservation of certain features and pointing out unique differences that regulate MDH enzyme activity.

Integrative residue-intuitive machine learning and MD Approach to Unveil Allosteric Site and Mechanism for β2AR

Allosteric drugs offer a new avenue for modern drug design. However, the identification of cryptic allosteric sites presents a formidable challenge. Following the allostery nature of residue-driven conformation transition, we propose a state-of-the-art computational pipeline by developing a residue-intuitive hybrid machine learning (RHML) model coupled with molecular dynamics (MD) simulation, through which we can efficiently identify the allosteric site and allosteric modulator as well as reveal their regulation mechanism. For the clinical target β2-adrenoceptor (β2AR), we discover an additional allosteric site located around residues D792.50, F2826.44, N3187.45 and S3197.46 and one putative allosteric modulator ZINC5042. Using Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and protein structure network (PSN), the allosteric potency and regulation mechanism are probed to further improve identification accuracy. Benefiting from sufficient computational evidence, the experimental assays then validate our predicted allosteric site, negative allosteric potency and regulation pathway, showcasing the effectiveness of the identification pipeline in practice. We expect that it will be applicable to other target proteins.

Engineering salt-tolerant Cas12f1 variants for gene-editing applications

CRISPR has revolutionized the field of genome editing in life sciences by serving as a versatile and state-of-the-art tool. Cas12f1 is a small nuclease of the bacterial immunity CRISPR system with an ideal size for cellular delivery, in contrast to CRISPR-associated (Cas) proteins like Cas9 or Cas12. However, Cas12f1 works best at low salt concentrations. In this study, we find that the plasticity of certain Cas12f1 regions (K196-Y202 and I452-L515) is negatively affected by increased salt concentrations. On this line, key protein domains (REC1, WED, Nuc, lid) that are involved in the DNA-target recognition and the activation of the catalytic RuvC domain are in turn also affected. We suggest that salt concentration should be taken in to consideration for activity assessments of Cas engineered variants, especially if the mutations are on the protospacer adjacent motif interacting domain. The results can be exploited for the engineering of Cas variants and the assessment of their activity at varying salt concentrations. We propose that the K198Q mutation can restore at great degree the compromised plasticity and could potentially lead to salt-tolerant Cas12f1 variants. The methodology can be also employed for the study of biomolecules in terms of their salinity tolerance.Communicated by Ramaswamy H. Sarma.

Molecular Dynamics Simulation Combined with Neural Relationship Inference and Markov Model to Reveal the Relationship between Conformational Regulation and Bioluminescence Properties of Gaussia Luciferase

Gaussia luciferase (Gluc) is currently known as the smallest naturally secreted luciferase. Due to its small molecular size, high sensitivity, short half-life, and high secretion efficiency, it has become an ideal reporter gene and is widely used in monitoring promoter activity, studying protein-protein interactions, protein localization, high-throughput drug screening, and real-time monitoring of tumor occurrence and development. Although studies have shown that different Gluc mutations exhibit different bioluminescent properties, their mechanisms have not been further investigated. The purpose of this study is to reveal the relationship between the conformational changes of Gluc mutants and their bioluminescent properties through molecular dynamics simulation combined with neural relationship inference (NRI) and Markov models. Our results indicate that, after binding to the luciferin coelenterazine (CTZ), the α-helices of the 109-119 residues of the Gluc Mutant2 (GlucM2, the flash-type mutant) are partially unraveled, while the α-helices of the same part of the Gluc Mutant1 (GlucM1, the glow-type mutant) are clearly formed. The results of Markov flux analysis indicate that the conformational differences between glow-type and flash-type mutants when combined with luciferin substrate CTZ mainly involve the helicity change of α7. The most representative conformation and active pocket distance analysis indicate that compared to the flash-type mutant GlucM2, the glow-type mutant GlucM1 has a higher degree of active site closure and tighter binding. In summary, we provide a theoretical basis for exploring the relationship between the conformational changes of Gluc mutants and their bioluminescent properties, which can serve as a reference for the modification and evolution of luciferases.

The early communication stages between serine proteases and enterovirus capsids in the race for viral disintegration

Serine proteases are important environmental contributors of enterovirus biocontrol. However, the structural features of molecular interaction accounting for the susceptibility of enteroviruses to proteases remains unexplained. Here, we describe the molecular mechanisms involved in the recruitment of serine proteases to viral capsids. Among the virus types used, coxsackievirus A9 (CVA9), but not CVB5 and echovirus 11 (E11), was inactivated by Subtilisin A in a host-independent manner, while Bovine Pancreatic Trypsin (BPT) only reduced CVA9 infectivity in a host-dependent manner. Predictive interaction models of each protease with capsid protomers indicate the main targets as internal disordered protein (IDP) segments exposed either on the 5-fold vertex (DE loop VP1) or at the 5/2-fold intersection (C-terminal end VP1) of viral capsids. We further show that a functional binding protease/capsid depends on both the strength and the evolution over time of protease-VP1 complexes, and lastly on the local adaptation of proteases on surrounding viral regions. Finally, we predicted three residues on CVA9 capsid that trigger cleavage by Subtilisin A, one of which may act as a sensor residue contributing to enzyme recognition on the DE loop. Overall, this study describes an important biological mechanism involved in enteroviruses biocontrol.

Active Discovery of the Allosteric Inhibitor Targeting Botrytis cinerea Chitinase Based on Neural Relational Inference for Food Preservation

Currently, allosteric inhibitors have emerged as an effective strategy in the development of preservatives against the drug-resistant Botrytis cinerea (B. cinerea). However, their passively driven development efficiency has proven challenging to meet the practical demands. Here, leveraging the deep learning Neural Relational Inference (NRI) framework, we actively identified an allosteric inhibitor targeting B. cinerea Chitinase, namely, 2-acetonaphthone. 2-Acetonaphthone binds to the crucial domain of Chitinase, forming the strong interaction with the allosteric sites. Throughout the interaction process, 2-acetonaphthone diminished the overall connectivity of the protein, inducing conformational changes. These findings align with the results obtained from Chitinase activity experiments, revealing an IC50 value of 67.6 μg/mL. Moreover, 2-acetonaphthone exhibited outstanding anti-B. cinerea activity by inhibiting Chitinase. In the gray mold infection model, 2-acetonaphthone significantly extended the preservation time of cherry tomatoes, positioning it as a promising preservative for fruit storage.

Allosteric pathways of SARS and SARS-CoV-2 spike protein identified by neural relational inference

The receptor binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein must undergo a crucial conformational transition to invade human cells. It is intriguing that this transition is accompanied by a synchronized movement of the entire spike protein. Therefore, it is possible to design allosteric regulators targeting non-RBD but hindering the conformational transition of RBD. To understand the allosteric mechanism in detail, we establish a computational framework by integrating coarse-grained molecular dynamic simulations and a state-of-the-art neural network model called neural relational inference. Leveraging this framework, we have elucidated the allosteric pathway of the SARS-CoV-2 spike protein at the residue level and identified the molecular mechanisms involved in the transmission of allosteric signals. The movement of D614 is coupled with that of Q321. This interaction subsequently influences the movement of K528/K529, ultimately coupling with the movement of RBD during conformational changes. Mutations that weaken the interactions within this pathway naturally block the allosteric signal transmission, thereby modulating the conformational transitions. This observation also offers a rationale for the distinct allosteric patterns observed in the SARS-CoV spike protein. Our result provides a useful method for analyzing the dynamics of potential viral variants in the future.

Combined experiments and molecular simulations for understanding the thermo-responsive behavior and gelation of methylated glucans with different glycosidic linkages

HYPOTHESIS

Elucidation of the micro-mechanisms of sol-gel transition of gelling glucans with different glycosidic linkages is crucial for understanding their structure-property relationship and for various applications. Glucans with distinct molecular chain structures exhibit unique gelation behaviors. The disparate gelation phenomena observed in two methylated glucans, methylated (1,3)-β-d-glucan of curdlan (MECD) and methylated (1,4)-β-d-glucan of cellulose (MC), notwithstanding their equivalent degrees of substitution, are intricately linked to their unique molecular architectures and interactions between glucan and water.

EXPERIMENTS

Density functional theory and molecular dynamics simulations focused on the electronic property distinctions between MECD and MC, alongside conformational variations during thermal gelation. Inline attenuated total reflection Fourier transform infrared spectroscopy tracked secondary structure alterations in MECD and MC. To corroborate the simulation results, additional analyses including circular dichroism, rheology, and micro-differential scanning calorimetry were performed.

FINDINGS

Despite having similar thermally induced gel networks, MECD and MC display distinct physical gelation patterns and molecular-level conformational changes during gelation. The network of MC gel was formed via a "coil-to-ring" transition, followed by ring stacking. In contrast, the MECD gel comprised compact irregular helices accompanied by notable volume shrinkage. These variations in gelation behavior are ascribed to heightened hydrophobic interactions and diminished hydrogen bonding in both systems upon heating, resulting in gelation. These findings provide valuable insights into the microstructural changes during gelation and the thermo-gelation mechanisms of structurally similar polysaccharides.

Adaptive Workflows of Machine Learning Illuminate the Sequential Operation Mechanism of the TAK1's Allosteric Network

Allostery is a fundamental mechanism driving biomolecular processes that holds significant therapeutic concern. Our study rigorously investigates how two distinct machine-learning algorithms uniquely classify two already close-to-active DFG-in states of TAK1, differing just by the presence or absence of its allosteric activator TAB1, from an ensemble mixture of conformations (obtained from 2.4 μs molecular dynamics (MD) simulations). The novelty, however, lies in understanding the deeper algorithmic potentials to systematically derive a diverse set of differential residue connectivity features that reconstruct the essential mechanistic architecture for TAK1-TAB1 allostery in such a close-to-active biochemical scenario. While the recursive, random forest-based workflow displays the potential of conducting discretized, hierarchical derivation of allosteric features, a multilayer perceptron-based approach gains considerable efficacy in revealing fluid connected patterns of features when hybridized with mutual information scoring. Interestingly, both pipelines benchmark similar directions of functional conformational changes for TAK1's activation. The findings significantly advance the depth of mechanistic understanding by highlighting crucial activation signatures along a directed C-lobe → activation loop → ATP pocket channel of information flow, including (1) the αF-αE biterminal alignments and (2) the "catalytic" drift of the activation loop toward kinase active site. Besides, some novel allosteric hotspots (K253, Y206, N189, etc.) are further recognized as TAB1 sensors, transducers, and responders, including a benchmark E70 mutation site, precisely mapping the important structural segments for sequential allosteric execution. Hence, our work demonstrates how to navigate through greater structural depths and dimensions of dynamic allosteric machineries just by leveraging standard ML methods in suitable streamlined workflows adaptive to the specific system and objectives.

A comprehensive survey on protein-ligand binding site prediction

Protein-ligand binding site prediction is critical for protein function annotation and drug discovery. Biological experiments are time-consuming and require significant equipment, materials, and labor resources. Developing accurate and efficient computational methods for protein-ligand interaction prediction is essential. Here, we summarize the key challenges associated with ligand binding site (LBS) prediction and introduce recently published methods from their input features, computational algorithms, and ligand types. Furthermore, we investigate the specificity of allosteric site identification as a particular LBS type. Finally, we discuss the prospective directions for machine learning-based LBS prediction in the near future.

Characterization of binding kinetics and intracellular signaling of new psychoactive substances targeting cannabinoid receptor using transition-based reweighting method

New psychoactive substances (NPS) targeting cannabinoid receptor 1 pose a significant threat to society as recreational abusive drugs that have pronounced physiological side effects. These greater adverse effects compared to classical cannabinoids have been linked to the higher downstream β-arrestin signaling. Thus, understanding the mechanism of differential signaling will reveal important structure-activity relationship essential for identifying and potentially regulating NPS molecules. In this study, we simulate the slow (un)binding process of NPS MDMB-Fubinaca and classical cannabinoid HU-210 from CB1 using multi-ensemble simulation to decipher the effects of ligand binding dynamics on downstream signaling. The transition-based reweighing method is used for the estimation of transition rates and underlying thermodynamics of (un)binding processes of ligands with nanomolar affinities. Our analyses reveal major interaction differences with transmembrane TM7 between NPS and classical cannabinoids. A variational autoencoder-based approach, neural relational inference (NRI), is applied to assess the allosteric effects on intracellular regions attributable to variations in binding pocket interactions. NRI analysis indicate a heightened level of allosteric control of NPxxY motif for NPS-bound receptors, which contributes to the higher probability of formation of a crucial triad interaction (Y7.53-Y5.58-T3.46) necessary for stronger β-arrestin signaling. Hence, in this work, MD simulation, data-driven statistical methods, and deep learning point out the structural basis for the heightened physiological side effects associated with NPS, contributing to efforts aimed at mitigating their public health impact.

Exploring biased activation characteristics by molecular dynamics simulation and machine learning for the μ-opioid receptor

Biased ligands selectively activating specific downstream signaling pathways (termed as biased activation) exhibit significant therapeutic potential. However, the conformational characteristics revealed are very limited for the biased activation, which is not conducive to biased drug development. Motivated by the issue, we combine extensive accelerated molecular dynamics simulations and an interpretable deep learning model to probe the biased activation features for two complex systems constructed by the inactive μOR and two different biased agonists (G-protein-biased agonist TRV130 and β-arrestin-biased agonist endomorphin2). The results indicate that TRV130 binds deeper into the receptor core compared to endomorphin2, located between W2936.48 and D1142.50, and forms hydrogen bonding with D1142.50, while endomorphin2 binds above W2936.48. The G protein-biased agonist induces greater outward movements of the TM6 intracellular end, forming a typical active conformation, while the β-arrestin-biased agonist leads to a smaller extent of outward movements of TM6. Compared with TRV130, endomorphin2 causes more pronounced inward movements of the TM7 intracellular end and more complex conformational changes of H8 and ICL1. In addition, important residues determining the two different biased activation states were further identified by using an interpretable deep learning classification model, including some common biased activation residues across Class A GPCRs like some key residues on the TM2 extracellular end, ECL2, TM5 intracellular end, TM6 intracellular end, and TM7 intracellular end, and some specific important residues of ICL3 for μOR. The observations will provide valuable information for understanding the biased activation mechanism for GPCRs.

Machine learning approaches in predicting allosteric sites

Allosteric regulation is a fundamental biological mechanism that can control critical cellular processes via allosteric modulator binding to protein distal functional sites. The advantages of allosteric modulators over orthosteric ones have sparked the development of numerous computational approaches, such as the identification of allosteric binding sites, to facilitate allosteric drug discovery. Building on the success of machine learning (ML) models for solving complex problems in biology and chemistry, several ML models for predicting allosteric sites have been developed. In this review, we provide an overview of these models and discuss future perspectives powered by the field of artificial intelligence such as protein language models.

Opportunities and Challenges for Machine Learning-Assisted Enzyme Engineering

Enzymes can be engineered at the level of their amino acid sequences to optimize key properties such as expression, stability, substrate range, and catalytic efficiency-or even to unlock new catalytic activities not found in nature. Because the search space of possible proteins is vast, enzyme engineering usually involves discovering an enzyme starting point that has some level of the desired activity followed by directed evolution to improve its "fitness" for a desired application. Recently, machine learning (ML) has emerged as a powerful tool to complement this empirical process. ML models can contribute to (1) starting point discovery by functional annotation of known protein sequences or generating novel protein sequences with desired functions and (2) navigating protein fitness landscapes for fitness optimization by learning mappings between protein sequences and their associated fitness values. In this Outlook, we explain how ML complements enzyme engineering and discuss its future potential to unlock improved engineering outcomes.

Protein dynamics underlying allosteric regulation

Allostery is the mechanism by which information and control are propagated in biomolecules. It regulates ligand binding, chemical reactions, and conformational changes. An increasing level of experimental resolution and control over allosteric mechanisms promises a deeper understanding of the molecular basis for life and powerful new therapeutics. In this review, we survey the literature for an up-to-date biological and theoretical understanding of protein allostery. By delineating five ways in which the energy landscape or the kinetics of a system may change to give rise to allostery, we aim to help the reader grasp its physical origins. To illustrate this framework, we examine three systems that display these forms of allostery: allosteric inhibitors of beta-lactamases, thermosensation of TRP channels, and the role of kinetic allostery in the function of kinases. Finally, we summarize the growing power of computational tools available to investigate the different forms of allostery presented in this review.

Atomic surface achieved through a novel cross-scale model from macroscale to nanoscale

Chemical mechanical polishing (CMP) is widely used to achieve an atomic surface globally, yet its cross-scale polishing mechanisms are elusive. Moreover, traditional CMP normally employs toxic and corrosive slurries, resulting in potential pollution to the environment. To overcome these challenges, a novel cross-scale model from the millimeter to nanometer scale is proposed, which was confirmed by a newly developed green CMP process. The developed CMP slurry consisted of hydrogen peroxide, sodium carbonate, sodium hydroxycellulose, and silica. Prior to CMP, fused silica was polished by a ceria slurry. After CMP, the surface roughness (Sa) was 0.126 nm, the material-removal rate was 88.3 nm min-1, and the thickness of the damaged layer was 8.8 nm. The proposed model was built by fibers, through integrating Eulerian and Lagrangian models and reactive force field-molecular dynamics. The results predicted by the model were in good agreement with those of CMP experimentally. A model for large-sized fibers revealed that a direct contact area of 11.12% was obtained for a non-woven polishing pad during the CMP experiments. Another model constructed via combining Eulerian and Lagrangian functions showed that the stress at the intersections of the fibers varied mainly from 0.1 to 0.01 MPa and was higher than the stress at other parts. An increase in viscosity led to a decrease in the areas with low stress, demonstrating that viscosity enhanced the stress and facilitated the removal of material. At the microscale and nanoscale, the stress of the abrasive surface exposed to the workpiece changed from 2.21 to 6.43 GPa. Stress at the interface contributed to the formation of bridging bonds, further promoting the removal of material. With increasing the compressive stress, the material-removal form was transformed from a single atom to molecular chains. The proposed model and developed green CMP offer new insights to understand the cross-scale polishing mechanism, as well as for designing and manufacturing novel polishing slurries, pads, and setups.

Evolution shapes interaction patterns for epistasis and specific protein binding in a two-component signaling system

The elegant design of protein sequence/structure/function relationships arises from the interaction patterns between amino acid positions. A central question is how evolutionary forces shape the interaction patterns that encode long-range epistasis and binding specificity. Here, we combined family-wide evolutionary analysis of natural homologous sequences and structure-oriented evolution simulation for two-component signaling (TCS) system. The magnitude-frequency relationship of coupling conservation between positions manifests a power-law-like distribution and the positions with highly coupling conservation are sparse but distributed intensely on the binding surfaces and hydrophobic core. The structure-specific interaction pattern involves further optimization of local frustrations at or near the binding surface to adapt the binding partner. The construction of family-wide conserved interaction patterns and structure-specific ones demonstrates that binding specificity is modulated by both direct intermolecular interactions and long-range epistasis across the binding complex. Evolution sculpts the interaction patterns via sequence variations at both family-wide and structure-specific levels for TCS system.

Using deep learning and molecular dynamics simulations to unravel the regulation mechanism of peptides as noncompetitive inhibitor of xanthine oxidase

Xanthine oxidase (XO) is a crucial enzyme in the development of hyperuricemia and gout. This study focuses on LWM and ALPM, two food-derived inhibitors of XO. We used molecular docking to obtain three systems and then conducted 200 ns molecular dynamics simulations for the Apo, LWM, and ALPM systems. The results reveal a stronger binding affinity of the LWM peptide to XO, potentially due to increased hydrogen bond formation. Notable changes were observed in the XO tunnel upon inhibitor binding, particularly with LWM, which showed a thinner, longer, and more twisted configuration compared to ALPM. The study highlights the importance of residue F914 in the allosteric pathway. Methodologically, we utilized the perturbed response scan (PRS) based on Python, enhancing tools for MD analysis. These findings deepen our understanding of food-derived anti-XO inhibitors and could inform the development of food-based therapeutics for reducing uric acid levels with minimal side effects.

Modifying the amino acids in conformational motion pathway of the α-amylase of Geobacillus stearothermophilus improved its activity and stability

Amino acids along the conformational motion pathway of the enzyme molecule correlated to its flexibility and rigidity. To enhance the enzyme activity and thermal stability, the motion pathway of Geobacillus stearothermophilus α-amylase has been identified and molecularly modified by using the neural relational inference model and deep learning tool. The significant differences in substrate specificity, enzymatic kinetics, optimal temperature, and thermal stability were observed among the mutants with modified amino acids along the pathway. Mutants especially the P44E demonstrated enhanced hydrolytic activity and catalytic efficiency (kcat/KM) than the wild-type enzyme to 95.0% and 93.8% respectively, with the optimum temperature increased to 90°C. This mutation from proline to glutamic acid has increased the number and the radius of the bottleneck of the channels, which might facilitate transporting large starch substrates into the enzyme. The mutation could also optimize the hydrogen bonding network of the catalytic center, and diminish the spatial hindering to the substrate entry and exit from the catalytic center.

Progressive computational approaches to facilitate decryption of allosteric mechanism and drug discovery

Allostery is a ubiquitous biological phenomenon where perturbation at topologically distal areas of a protein serves as a trigger to fine-tune the orthosteric site and thus regulate protein function. The investigation of allosteric regulation greatly enhances our understanding of human diseases and broadens avenue for drug discovery. For decades, owing to the difficulty in allostery characterization through serendipitous experimental screening, researchers have developed several innovative computational approaches, which proves to accelerate the elucidation of allostery. Herein, we review the state-of-the-art advance of computational methodologies for allostery study, with particular emphasis on promising trends emerging over the past two years. We expect this review will outline the latest landscape of allostery study and inspire researchers to further facilitate this field.

Revealing Allosteric Mechanism of Amino Acid Binding Proteins from Open to Closed State

Amino acid binding proteins (AABPs) undergo significant conformational closure in the periplasmic space of Gram-negative bacteria, tightly binding specific amino acid substrates and then initiating transmembrane transport of nutrients. Nevertheless, the possible closure mechanisms after substrate binding, especially long-range signaling, remain unknown. Taking three typical AABPs-glutamine binding protein (GlnBP), histidine binding protein (HisJ) and lysine/arginine/ornithine binding protein (LAOBP) in Escherichia coli (E. coli)-as research subjects, a series of theoretical studies including sequence alignment, Gaussian network model (GNM), anisotropic network model (ANM), conventional molecular dynamics (cMD) and neural relational inference molecular dynamics (NRI-MD) simulations were carried out. Sequence alignment showed that GlnBP, HisJ and LAOBP have high structural similarity. According to the results of the GNM and ANM, AABPs' Index Finger and Thumb domains exhibit closed motion tendencies that contribute to substrate capture and stable binding. Based on cMD trajectories, the Index Finger domain, especially the I-Loop region, exhibits high molecular flexibility, with residues 11 and 117 both being potentially key residues for receptor-ligand recognition and initiation of receptor allostery. Finally, the signaling pathway of AABPs' conformational closure was revealed by NRI-MD training and trajectory reconstruction. This work not only provides a complete picture of AABPs' recognition mechanism and possible conformational closure, but also aids subsequent structure-based design of small-molecule oncology drugs.

Protein Structure Prediction: Challenges, Advances, and the Shift of Research Paradigms

Protein structure prediction is an interdisciplinary research topic that has attracted researchers from multiple fields, including biochemistry, medicine, physics, mathematics, and computer science. These researchers adopt various research paradigms to attack the same structure prediction problem: biochemists and physicists attempt to reveal the principles governing protein folding; mathematicians, especially statisticians, usually start from assuming a probability distribution of protein structures given a target sequence and then find the most likely structure, while computer scientists formulate protein structure prediction as an optimization problem - finding the structural conformation with the lowest energy or minimizing the difference between predicted structure and native structure. These research paradigms fall into the two statistical modeling cultures proposed by Leo Breiman, namely, data modeling and algorithmic modeling. Recently, we have also witnessed the great success of deep learning in protein structure prediction. In this review, we present a survey of the efforts for protein structure prediction. We compare the research paradigms adopted by researchers from different fields, with an emphasis on the shift of research paradigms in the era of deep learning. In short, the algorithmic modeling techniques, especially deep neural networks, have considerably improved the accuracy of protein structure prediction; however, theories interpreting the neural networks and knowledge on protein folding are still highly desired.

Study on the allosteric activation mechanism of SHP2 via elastic network models and neural relational inference molecular dynamics simulation

As a ubiquitous protein tyrosine phosphatase, SHP2 is involved in PD-1/PD-L1 mediated tumor immune escape and undergoes substantial conformational changes. Therefore, it is considered an ideal target for tumor intervention. However, the allosteric mechanisms of SHP2 binding PD-1 intracellular ITIM/ITSM phosphopeptides remain unclear, which greatly hinders the development of novel structure-based anticancer allosteric inhibitors. In this work, the open and closed structural models of SHP2 are first constructed based on this knowledge; next their motion modes are investigated via elastic network models such as the Gaussian network model (GNM), anisotropic network model (ANM) and adaptive anisotropic network model (aANM); and finally, a possible allosteric signaling pathway is proposed using a neural relational inference molecular dynamics (NRI-MD) simulation embedded with an artificial intelligence (AI) strategy. In GNM and ANM, the N-SH2, C-SH2 and PTP domains all exhibit distinct dynamics partitions, and the N-SH2/C-SH2 regions show a rigid rotation relative to PTP. According to a series of intermediate snapshots given by aANM, N-SH2 is first identified with pY223 specifically, inducing a D'E-loop to change from β-sheets to random coils, and then, C-SH2 serves as a fulcrum to drive N-SH2 to rotate 110° completely away from the original active sites of PTP. Finally, a possible allosteric signaling-transfer path for SHP2, namely R220-R138-T108-R32, is proposed based on NRI-MD sampling. This work provides a possible allosteric mechanism of SHP2, which is helpful for the following design of novel allosteric inhibitors and is expected to be used in clinical synergies with PD-1 monoclonal antibody.

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

The Motion Paradigm of Pre-Dock Zearalenone Hydrolase Predictions with Molecular Dynamics and the Docking Phase with Umbrella Sampling

Zearalenone (ZEN) is one of the most prevalent estrogenic mycotoxins, is produced mainly by the Fusarium family of fungi, and poses a risk to the health of animals. Zearalenone hydrolase (ZHD) is an important enzyme capable of degrading ZEN into a non-toxic compound. Although previous research has investigated the catalytic mechanism of ZHD, information on its dynamic interaction with ZEN remains unknown. This study aimed to develop a pipeline for identifying the allosteric pathway of ZHD. Using an identity analysis, we identified hub genes whose sequences can generalize a set of sequences in a protein family. We then utilized a neural relational inference (NRI) model to identify the allosteric pathway of the protein throughout the entire molecular dynamics simulation. The production run lasted 1 microsecond, and we analyzed residues 139-222 for the allosteric pathway using the NRI model. We found that the cap domain of the protein opened up during catalysis, resembling a hemostatic tape. We used umbrella sampling to simulate the dynamic docking phase of the ligand-protein complex and found that the protein took on a square sandwich shape. Our energy analysis, using both molecular mechanics/Poisson-Boltzmann (Generalized-Born) surface area (MMPBSA) and Potential Mean Force (PMF) analysis, showed discrepancies, with scores of -8.45 kcal/mol and -1.95 kcal/mol, respectively. MMPBSA, however, obtained a similar score to that of a previous report.

Structural Changes, Biological Consequences, and Repurposing of Colchicine Site Ligands

Microtubule-targeting agents (MTAs) bind to one of several distinct sites in the tubulin dimer, the subunit of microtubules. The binding affinities of MTAs may vary by several orders of magnitude, even for MTAs that specifically bind to a particular site. The first drug binding site discovered in tubulin was the colchicine binding site (CBS), which has been known since the discovery of the tubulin protein. Although highly conserved throughout eukaryotic evolution, tubulins show diversity in their sequences between tubulin orthologs (inter-species sequence differences) and paralogs (intraspecies differences, such as tubulin isotypes). The CBS is promiscuous and binds to a broad range of structurally distinct molecules that can vary in size, shape, and affinity. This site remains a popular target for the development of new drugs to treat human diseases (including cancer) and parasitic infections in plants and animals. Despite the rich knowledge about the diversity of tubulin sequences and the structurally distinct molecules that bind to the CBS, a pattern has yet to be found to predict the affinity of new molecules that bind to the CBS. In this commentary, we briefly discuss the literature evidencing the coexistence of the varying binding affinities for drugs that bind to the CBS of tubulins from different species and within species. We also comment on the structural data that aim to explain the experimental differences observed in colchicine binding to the CBS of β-tubulin class VI (TUBB1) compared to other isotypes.

Inverse problem beyond two-body interaction: The cubic mean-field Ising model

In this paper, we solve the inverse problem for the cubic mean-field Ising model. Starting from configuration data generated according to the distribution of the model, we reconstruct the free parameters of the system. We test the robustness of this inversion procedure both in the region of uniqueness of the solutions and in the region where multiple thermodynamics phases are present.

Accelerating therapeutic protein design with computational approaches toward the clinical stage

Therapeutic protein, represented by antibodies, is of increasing interest in human medicine. However, clinical translation of therapeutic protein is still largely hindered by different aspects of developability, including affinity and selectivity, stability and aggregation prevention, solubility and viscosity reduction, and deimmunization. Conventional optimization of the developability with widely used methods, like display technologies and library screening approaches, is a time and cost-intensive endeavor, and the efficiency in finding suitable solutions is still not enough to meet clinical needs. In recent years, the accelerated advancement of computational methodologies has ushered in a transformative era in the field of therapeutic protein design. Owing to their remarkable capabilities in feature extraction and modeling, the integration of cutting-edge computational strategies with conventional techniques presents a promising avenue to accelerate the progression of therapeutic protein design and optimization toward clinical implementation. Here, we compared the differences between therapeutic protein and small molecules in developability and provided an overview of the computational approaches applicable to the design or optimization of therapeutic protein in several developability issues.

From Deep Mutational Mapping of Allosteric Protein Landscapes to Deep Learning of Allostery and Hidden Allosteric Sites: Zooming in on "Allosteric Intersection" of Biochemical and Big Data Approaches

The recent advances in artificial intelligence (AI) and machine learning have driven the design of new expert systems and automated workflows that are able to model complex chemical and biological phenomena. In recent years, machine learning approaches have been developed and actively deployed to facilitate computational and experimental studies of protein dynamics and allosteric mechanisms. In this review, we discuss in detail new developments along two major directions of allosteric research through the lens of data-intensive biochemical approaches and AI-based computational methods. Despite considerable progress in applications of AI methods for protein structure and dynamics studies, the intersection between allosteric regulation, the emerging structural biology technologies and AI approaches remains largely unexplored, calling for the development of AI-augmented integrative structural biology. In this review, we focus on the latest remarkable progress in deep high-throughput mining and comprehensive mapping of allosteric protein landscapes and allosteric regulatory mechanisms as well as on the new developments in AI methods for prediction and characterization of allosteric binding sites on the proteome level. We also discuss new AI-augmented structural biology approaches that expand our knowledge of the universe of protein dynamics and allostery. We conclude with an outlook and highlight the importance of developing an open science infrastructure for machine learning studies of allosteric regulation and validation of computational approaches using integrative studies of allosteric mechanisms. The development of community-accessible tools that uniquely leverage the existing experimental and simulation knowledgebase to enable interrogation of the allosteric functions can provide a much-needed boost to further innovation and integration of experimental and computational technologies empowered by booming AI field.

Machine learning and protein allostery

The fundamental biological importance and complexity of allosterically regulated proteins stem from their central role in signal transduction and cellular processes. Recently, machine-learning approaches have been developed and actively deployed to facilitate theoretical and experimental studies of protein dynamics and allosteric mechanisms. In this review, we survey recent developments in applications of machine-learning methods for studies of allosteric mechanisms, prediction of allosteric effects and allostery-related physicochemical properties, and allosteric protein engineering. We also review the applications of machine-learning strategies for characterization of allosteric mechanisms and drug design targeting SARS-CoV-2. Continuous development and task-specific adaptation of machine-learning methods for protein allosteric mechanisms will have an increasingly important role in bridging a wide spectrum of data-intensive experimental and theoretical technologies.

The quaternary question: Determining allostery in spastin through dynamics classification learning and bioinformatics

The nanomachine from the ATPases associated with various cellular activities superfamily, called spastin, severs microtubules during cellular processes. To characterize the functionally important allostery in spastin, we employed methods from evolutionary information, to graph-based networks, to machine learning applied to atomistic molecular dynamics simulations of spastin in its monomeric and the functional hexameric forms, in the presence or absence of ligands. Feature selection, using machine learning approaches, for transitions between spastin states recognizes all the regions that have been proposed as allosteric or functional in the literature. The analysis of the composition of the Markov State Model macrostates in the spastin monomer, and the analysis of the direction of change in the top machine learning features for the transitions, indicate that the monomer favors the binding of ATP, which primes the regions involved in the formation of the inter-protomer interfaces for binding to other protomer(s). Allosteric path analysis of graph networks, built based on the cross-correlations between residues in simulations, shows that perturbations to a hub specific for the pre-hydrolysis hexamer propagate throughout the structure by passing through two obligatory regions: the ATP binding pocket, and pore loop 3, which connects the substrate binding site to the ATP binding site. Our findings support a model where the changes in the terminal protomers due to the binding of ligands play an active role in the force generation in spastin. The secondary structures in spastin, which are found to be highly degenerative within the network paths, are also critical for feature transitions of the classification models, which can guide the design of allosteric effectors to enhance or block allosteric signaling.

Exploring and Learning the Universe of Protein Allostery Using Artificial Intelligence Augmented Biophysical and Computational Approaches

Allosteric mechanisms are commonly employed regulatory tools used by proteins to orchestrate complex biochemical processes and control communications in cells. The quantitative understanding and characterization of allosteric molecular events are among major challenges in modern biology and require integration of innovative computational experimental approaches to obtain atomistic-level knowledge of the allosteric states, interactions, and dynamic conformational landscapes. The growing body of computational and experimental studies empowered by emerging artificial intelligence (AI) technologies has opened up new paradigms for exploring and learning the universe of protein allostery from first principles. In this review we analyze recent developments in high-throughput deep mutational scanning of allosteric protein functions; applications and latest adaptations of Alpha-fold structural prediction methods for studies of protein dynamics and allostery; new frontiers in integrating machine learning and enhanced sampling techniques for characterization of allostery; and recent advances in structural biology approaches for studies of allosteric systems. We also highlight recent computational and experimental studies of the SARS-CoV-2 spike (S) proteins revealing an important and often hidden role of allosteric regulation driving functional conformational changes, binding interactions with the host receptor, and mutational escape mechanisms of S proteins which are critical for viral infection. We conclude with a summary and outlook of future directions suggesting that AI-augmented biophysical and computer simulation approaches are beginning to transform studies of protein allostery toward systematic characterization of allosteric landscapes, hidden allosteric states, and mechanisms which may bring about a new revolution in molecular biology and drug discovery.

Protein structure and folding pathway prediction based on remote homologs recognition using PAthreader

Recognition of remote homologous structures is a necessary module in AlphaFold2 and is also essential for the exploration of protein folding pathways. Here, we propose a method, PAthreader, to recognize remote templates and explore folding pathways. Firstly, we design a three-track alignment between predicted distance profiles and structure profiles extracted from PDB and AlphaFold DB, to improve the recognition accuracy of remote templates. Secondly, we improve the performance of AlphaFold2 using the templates identified by PAthreader. Thirdly, we explore protein folding pathways based on our conjecture that dynamic folding information of protein is implicitly contained in its remote homologs. The results show that the average accuracy of PAthreader templates is 11.6% higher than that of HHsearch. In terms of structure modelling, PAthreader outperform AlphaFold2 and ranks first on the CAMEO blind test for the latest three months. Furthermore, we predict protein folding pathways for 37 proteins, in which the results of 7 proteins are almost consistent with those of biological experiments, and the other 30 human proteins have yet to be verified by biological experiments, revealing that folding information can be exploited from remote homologous structures.

Mutually beneficial confluence of structure-based modeling of protein dynamics and machine learning methods

Proteins sample an ensemble of conformers under physiological conditions, having access to a spectrum of modes of motions, also called intrinsic dynamics. These motions ensure the adaptation to various interactions in the cell, and largely assist in, if not determine, viable mechanisms of biological function. In recent years, machine learning frameworks have proven uniquely useful in structural biology, and recent studies further provide evidence to the utility and/or necessity of considering intrinsic dynamics for increasing their predictive ability. Efficient quantification of dynamics-based attributes by recently developed physics-based theories and models such as elastic network models provides a unique opportunity to generate data on dynamics for training ML models towards inferring mechanisms of protein function, assessing pathogenicity, or estimating binding affinities.

DNNGP, a deep neural network-based method for genomic prediction using multi-omics data in plants

Genomic prediction is an effective way to accelerate the rate of agronomic trait improvement in plants. Traditional methods typically use linear regression models with clear assumptions; such methods are unable to capture the complex relationships between genotypes and phenotypes. Non-linear models (e.g., deep neural networks) have been proposed as a superior alternative to linear models because they can capture complex non-additive effects. Here we introduce a deep learning (DL) method, deep neural network genomic prediction (DNNGP), for integration of multi-omics data in plants. We trained DNNGP on four datasets and compared its performance with methods built with five classic models: genomic best linear unbiased prediction (GBLUP); two methods based on a machine learning (ML) framework, light gradient boosting machine (LightGBM) and support vector regression (SVR); and two methods based on a DL framework, deep learning genomic selection (DeepGS) and deep learning genome-wide association study (DLGWAS). DNNGP is novel in five ways. First, it can be applied to a variety of omics data to predict phenotypes. Second, the multilayered hierarchical structure of DNNGP dynamically learns features from raw data, avoiding overfitting and improving the convergence rate using a batch normalization layer and early stopping and rectified linear activation (rectified linear unit) functions. Third, when small datasets were used, DNNGP produced results that are competitive with results from the other five methods, showing greater prediction accuracy than the other methods when large-scale breeding data were used. Fourth, the computation time required by DNNGP was comparable with that of commonly used methods, up to 10 times faster than DeepGS. Fifth, hyperparameters can easily be batch tuned on a local machine. Compared with GBLUP, LightGBM, SVR, DeepGS and DLGWAS, DNNGP is superior to these existing widely used genomic selection (GS) methods. Moreover, DNNGP can generate robust assessments from diverse datasets, including omics data, and quickly incorporate complex and large datasets into usable models, making it a promising and practical approach for straightforward integration into existing GS platforms.

Origin of Correlations between Local Conformational States of Consecutive Amino Acid Residues and Their Role in Shaping Protein Structures and in Allostery

By analyzing the Kubo-cluster-cumulant expansion of the potential of mean force of polypeptide chains corresponding to backbone-local interactions averaged over the rotation of the peptide groups about the Cα···Cα virtual bonds, we identified two important kinds of "along-chain" correlations that pertain to extended chain segments bordered by turns (usually the β-strands) and to the folded spring-like segments (usually α-helices), respectively, and are expressed as multitorsional potentials. These terms affect the positioning of structural elements with respect to each other and, consequently, contribute to determining their packing. Additionally, for extended chain segments, the correlation terms contribute to propagating the conformational change at one end to the other end, which is characteristic of allosteric interactions. We confirmed both findings by statistical analysis of the virtual-bond geometry of 77 950 proteins. Augmenting coarse-grained and, possibly, all-atom force fields with these correlation terms could improve their capacity to model protein structure and dynamics.

Deciphering the QR Code of the CRISPR-Cas9 System: Synergy between Gln768 (Q) and Arg976 (R)

Markov state models (MSMs) and machine learning (ML) algorithms can extrapolate the long-time-scale behavior of large biomolecules from molecular dynamics (MD) trajectories. In this study, an MD-MSM-ML scheme has been applied to probe the large endonuclease (Cas9) in the bacterial adaptive immunity CRISPR-Cas9 system. CRISPR has become a programmable and state-of-the-art powerful genome editing tool that has already revolutionized life sciences. CRISPR-Cas9 is programmed to process specific DNA sequences in the genome. However, human/biomedical applications are compromised by off-target DNA damage. Characterization of Cas9 at the structural and biophysical levels is a prerequisite for the development of efficient and high-fidelity Cas9 variants. The Cas9 wild type and two variants (R63A-R66A-R70A, R69A-R71A-R74A-R78A) are studied herein. The configurational space of Cas9 is provided with a focus on the conformations of the side chains of two residues (Gln768 and Arg976). A model for the synergy between those two residues is proposed. The results are discussed within the context of experimental literature. The results and methodology can be exploited for the study of large biomolecules in general and for the engineering of more efficient and safer Cas9 variants for applications.

Protein Function Analysis through Machine Learning

Machine learning (ML) has been an important arsenal in computational biology used to elucidate protein function for decades. With the recent burgeoning of novel ML methods and applications, new ML approaches have been incorporated into many areas of computational biology dealing with protein function. We examine how ML has been integrated into a wide range of computational models to improve prediction accuracy and gain a better understanding of protein function. The applications discussed are protein structure prediction, protein engineering using sequence modifications to achieve stability and druggability characteristics, molecular docking in terms of protein-ligand binding, including allosteric effects, protein-protein interactions and protein-centric drug discovery. To quantify the mechanisms underlying protein function, a holistic approach that takes structure, flexibility, stability, and dynamics into account is required, as these aspects become inseparable through their interdependence. Another key component of protein function is conformational dynamics, which often manifest as protein kinetics. Computational methods that use ML to generate representative conformational ensembles and quantify differences in conformational ensembles important for function are included in this review. Future opportunities are highlighted for each of these topics.