Molecular basis of the activation and dissociation of dimeric PYL2 receptor in abscisic acid signaling

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.

Rationalizing diverse binding mechanisms to the same protein fold: insights for ligand recognition and biosensor design

The engineering of novel protein-ligand binding interactions, particularly for complex drug-like molecules, is an unsolved problem which could enable many practical applications of protein biosensors. In this work, we analyzed two engineer ed biosensors, derived from the plant hormone sensor PYR1, to recognize either the agrochemical mandipropamid or the synthetic cannabinoid WIN55,212-2. Using a combination of quantitative deep mutational scanning experiments and molecular dynamics simulations, we demonstrated that mutations at common positions can promote protein-ligand shape complementarity and revealed prominent differences in the electrostatic networks needed to complement diverse ligands. MD simulations indicate that both PYR1 protein-ligand complexes bind a single conformer of their target ligand that is close to the lowest free energy conformer. Computational design using a fixed conformer and rigid body orientation led to new WIN55,212-2 sensors with nanomolar limits of detection. This work reveals mechanisms by which the versatile PYR1 biosensor scaffold can bind diverse ligands. This work also provides computational methods to sample realistic ligand conformers and rigid body alignments that simplify the computational design of biosensors for novel ligands of interest.

Contribution of Signaling Partner Association to Strigolactone Receptor Selectivity

The parasitic plant witchweed, Striga hermonthica, results in agricultural losses of billions of dollars per year. It perceives its host via plant hormones called strigolactones, which act as germination stimulants for witchweed. Strigolactone signaling involves substrate binding to the strigolactone receptor, followed by substrate hydrolysis and a conformational change from an inactive, or open state, to an active, or closed state. In the active state, the receptor associates with a signaling partner, MAX2. Recently, it was shown that this MAX2 association process acts as a strong contributor to the uniquely high signaling activity observed in ShHTL7; however, it is unknown why ShHTL7 has enhanced MAX2 association affinity. Using an umbrella sampling molecular dynamics approach, we characterized the association processes of AtD14, ShHTL7, a mutant of ShHTL7, and ShHTL6 with MAX2 homologue OsD3. From these results, we show that ShHTL7 has an enhanced standard binding free energy of OsD3 compared to those of the other receptors. Additionally, our results suggest that the overall topology of the T2/T3 helix region is likely an important modulator of MAX2 binding. Thus, differences in MAX2 association, modulated by differences in the T2/T3 helix region, are a contributor to differences in signaling activity between different strigolactone receptors.

Resolving binding pathways and solvation thermodynamics of plant hormone receptors

Plant hormones are small molecules that regulate plant growth, development, and responses to biotic and abiotic stresses. They are specifically recognized by the binding site of their receptors. In this work, we resolved the binding pathways for eight classes of phytohormones (auxin, jasmonate, gibberellin, strigolactone, brassinosteroid, cytokinin, salicylic acid, and abscisic acid) to their canonical receptors using extensive molecular dynamics simulations. Furthermore, we investigated the role of water displacement and reorganization at the binding site of the plant receptors through inhomogeneous solvation theory. Our findings predict that displacement of water molecules by phytohormones contributes to free energy of binding via entropy gain and is associated with significant free energy barriers for most systems analyzed. Also, our results indicate that displacement of unfavorable water molecules in the binding site can be exploited in rational agrochemical design. Overall, this study uncovers the mechanism of ligand binding and the role of water molecules in plant hormone perception, which creates new avenues for agrochemical design to target plant growth and development.

A clathrin-related protein FaRRP1/SCD2 integrates ABA trafficking and signaling to regulate strawberry fruit ripening

Abscisic acid (ABA) is a critical regulator for nonclimacteric fruit ripening such as in the model plant of strawberry (Fragaria × ananassa). Although FaRRP1 is proposed to participate in clathrin-mediated endocytosis of ABA, its action molecular mechanisms in ABA signaling are not fully understood. Here, using our isolated FaRRP1 (ripening-regulation protein) and candidate ABA receptor FaPYL2 and FaABAR from strawberry fruit, a series of silico and molecular interaction analyses demonstrate that they all bind to ABA, and FaRRP1 binds both FaPYL2 and FaABAR; by contrast, the binding affinity of FaRRP1 to FaPYL2 is relatively higher. Interestingly, the binding of FaRRP1 to FaPYL2 and FaABAR affects the perception affinity to ABA. Furthermore, exogenous ABA application and FaRRP1 transgenic analyses confirm that FaRRP1 participates in clathrin-mediated endocytosis and vesicle transport. Importantly, FaRRP1, FaPYL2, and FaABAR all trigger the initiation of strawberry fruit ripening at physiological and molecular levels. In conclusion, FaRRP1 not only binds to ABA but also affects the binding affinity of FaPYL2 and FaABAR to ABA, thus promoting strawberry fruit ripening. Our findings provide novel insights into the role of FaRRP1 in ABA trafficking and signaling, at least in strawberry, a model plant for nonclimacteric fruit ripening.

Understanding Supramolecular Assembly of Supercharged Proteins

Ordered supramolecular assemblies have recently been created using electrostatic interactions between oppositely charged proteins. Despite recent progress, the fundamental mechanisms governing the assembly of oppositely supercharged proteins are not fully understood. Here, we use a combination of experiments and computational modeling to systematically study the supramolecular assembly process for a series of oppositely supercharged green fluorescent protein variants. We show that net charge is a sufficient molecular descriptor to predict the interaction fate of oppositely charged proteins under a given set of solution conditions (e.g., ionic strength), but the assembled supramolecular structures critically depend on surface charge distributions. Interestingly, our results show that a large excess of charge is necessary to nucleate assembly and that charged residues not directly involved in interprotein interactions contribute to a substantial fraction (∼30%) of the interaction energy between oppositely charged proteins via long-range electrostatic interactions. Dynamic subunit exchange experiments further show that relatively small, 16-subunit assemblies of oppositely charged proteins have kinetic lifetimes on the order of ∼10-40 min, which is governed by protein composition and solution conditions. Broadly, our results inform how protein supercharging can be used to create different ordered supramolecular assemblies from a single parent protein building block.

Integration of machine learning with computational structural biology of plants

Computational structural biology of proteins has developed rapidly in recent decades with the development of new computational tools and the advancement of computing hardware. However, while these techniques have widely been used to make advancements in human medicine, these methods have seen less utilization in the plant sciences. In the last several years, machine learning methods have gained popularity in computational structural biology. These methods have enabled the development of new tools which are able to address the major challenges that have hampered the wide adoption of the computational structural biology of plants. This perspective examines the remaining challenges in computational structural biology and how the development of machine learning techniques enables more in-depth computational structural biology of plants.

Activation Mechanism of Strigolactone Receptors and Its Impact on Ligand Selectivity between Host and Parasitic Plants

Parasitic weeds such as Striga have led to significant losses in agricultural productivity worldwide. These weeds use the plant hormone strigolactone as a germination stimulant. Strigolactone signaling involves substrate hydrolysis followed by a conformational change of the receptor to a "closed" or "active" state that associates with a signaling partner, MAX2/D3. Crystal structures of active and inactive AtD14 receptors have helped elucidate the structural changes involved in activation. However, the mechanism by which the receptor activates remains unknown. The ligand dependence of AtD14 activation has been disputed by mutagenesis studies showing that enzymatically inactive receptors are able to associate with MAX2 proteins. Furthermore, activation differences between strigolactone receptor in Striga, ShHTL7, and AtD14 could contribute to the high sensitivity to strigolactones exhibited by parasitic plants. Using molecular dynamics simulations, we demonstrate that both AtD14 and ShHTL7 could adopt an active conformation in the absence of ligand. However, ShHTL7 exhibits a higher population in the inactive apo state as compared to the AtD14 receptor. We demonstrate that this difference in inactive state population is caused by sequence differences between their D-loops and interactions with the catalytic histidine that prevent full binding pocket closure in ShHTL7. These results indicate that ligand hydrolysis would enhance the active state population by destabilizing the inactive state in ShHTL7 as compared to AtD14. We also show that the mechanism of activation is more concerted in AtD14 than in ShHTL7 and that the main barrier to activation in ShHTL7 is closing of the binding pocket.