Oligomerization of Monoamine Transporters.

Transporters of the monoamine transporter (MAT) family regulate the uptake of important neurotransmitters like dopamine, serotonin, and norepinephrine. The MAT family functions using the electrochemical gradient of ions across the membrane and comprises three transporters, dopamine transporter (DAT), serotonin transporter (SERT), and norepinephrine transporter (NET). MAT transporters have been observed to exist in monomeric states to higher-order oligomeric states. Structural features, allosteric modulation, and lipid environment regulate the oligomerization of MAT transporters. NET and SERT oligomerization are regulated by levels of PIP2 present in the membrane. The kink present in TM12 in the MAT family is crucial for dimer interface formation. Allosteric modulation in the dimer interface hinders dimer formation. Oligomerization also influences the transporters' function, trafficking, and regulation. This chapter will focus on recent studies on monoamine transporters and discuss the factors affecting their oligomerization and its impact on their function.

SWEET family transporters act as water conducting carrier proteins in plants.

Dedicated water channels are involved in the facilitated diffusion of water molecules across the cell membrane in plants. Transporter proteins are also known to transport water molecules along with substrates, however the molecular mechanism of water permeation is not well understood in plant transporters. Here, we show plant sugar transporters from the SWEET (Sugar Will Eventually be Exported Transporter) family act as water-conducting carrier proteins via a variety of passive and active mechanisms that allow diffusion of water molecules from one side of the membrane to the other. This study provides a molecular perspective on how plant membrane transporters act as water carrier proteins, a topic that has not been extensively explored in literature. Water permeation in membrane transporters could occur via four distinct mechanisms which form our hypothesis for water transport in SWEETs. These hypothesis are tested using molecular dynamics simulations of the outward-facing, occluded, and inward-facing state of AtSWEET1 to identify the water permeation pathways and the flux associated with them. The hydrophobic gates at the center of the transport tunnel act as a barrier that restricts water permeation. We have performed in silico single and double mutations of the hydrophobic gate residues to examine the changes in the water conductivity. Surprisingly, the double mutant allows the water permeation to the intracellular half of the membrane and forms a continuous water channel. These computational results are validated by experimentally examining the transport of hydrogen peroxide molecules by the AtSWEET family of transporters. We have also shown that the transport of hydrogen peroxide follows the similar mechanism as water transport in AtSWEET1. Finally, we conclude that similar water-conduction states are also present in other SWEET transporters due to the high sequence and structure conservation exhibited by this transporter family.

Deep mutational scanning reveals sequence to function constraints for SWEET family transporters.

Protein science is entering a transformative phase enabled by deep mutational scans that provide an unbiased view of the residue level interactions that mediate function. However, it has yet to be extensively used to characterize the mutational and evolutionary landscapes of plant proteins. Here, we apply the method to explore sequence-function relationships within the sugar transporter AtSWEET13. DMS results describe how mutational interrogation throughout different regions of the protein affects AtSWEET13 abundance and transport function. Our results identify novel transport-enhancing mutations that are validated using the FRET sensor assays. Extending DMS results to phylogenetic analyses reveal the role of transmembrane helix 4 (TM4) which makes the SWEET family transporters distinct from prokaryotic SemiSWEETs. We show that transmembrane helix 4 is intolerant to motif swapping with other clade-specific SWEET TM4 compositions, despite accommodating single point-mutations towards aromatic and charged polar amino acids. We further show that the transfer learning approaches based on physics and ML based In silico variant prediction tools have limited utility for engineering plant proteins as they were unable to reproduce our experimental results. We conclude that DMS can produce datasets which, when combined with the right predictive computational frameworks, can direct plant engineering efforts through derivative phenotype selection and evolutionary insights.

Interplay between phosphorylation and oligomerization tunes the conformational ensemble of SWEET transporters.

SWEET sugar transporters are desirable biotechnological targets for improving plant growth. One engineering strategy includes modulating how SWEET transporters are regulated. Phosphorylation and oligomerization have been shown to positively regulate SWEET function, leading to increased sugar transport activity. However, constitutive phosphorylation may not be beneficial to plant health under basal conditions. Structural and mechanistic understanding of the interplay between phosphorylation and oligomerization in functional regulation of SWEETs remains limited. Using extensive molecular dynamics simulations coupled with Markov state models, we demonstrate the thermodynamic and kinetic effects of SWEET phosphorylation and oligomerization using OsSWEET2b as a model. We report that the beneficial effects of these SWEET regulatory mechanisms bias outward-facing states and improved extracellular gating, which complement published experimental findings. Our results offer molecular insights to SWEET regulation and may guide engineering strategies throughout the SWEET transport family.

The Arabidopsis AtSWEET13 transporter discriminates sugars by selective facial and positional substrate recognition.

Transporters are targeted by endogenous metabolites and exogenous molecules to reach cellular destinations, but it is generally not understood how different substrate classes exploit the same transporter's mechanism. Any disclosure of plasticity in transporter mechanism when treated with different substrates becomes critical for developing general selectivity principles in membrane transport catalysis. Using extensive molecular dynamics simulations with an enhanced sampling approach, we select the Arabidopsis sugar transporter AtSWEET13 as a model system to identify the basis for glucose versus sucrose molecular recognition and transport. Here we find that AtSWEET13 chemical selectivity originates from a conserved substrate facial selectivity demonstrated when committing alternate access, despite mono-/di-saccharides experiencing differing degrees of conformational and positional freedom throughout other stages of transport. However, substrate interactions with structural hallmarks associated with known functional annotations can help reinforce selective preferences in molecular transport.

Energetics of substrate transport in proton-dependent oligopeptide transporters.

The PepT So transporter mediates the transport of peptides across biological membranes. Despite advancements in structural biology, including cryogenic electron microscopy structures resolving PepT So in different states, the molecular basis of peptide recognition and transport by PepT So is not fully elucidated. In this study, we employed molecular dynamics simulations, Markov State Models (MSMs), and Transition Path Theory (TPT) to investigate the transport mechanism of an alanine-alanine peptide (Ala-Ala) through the PepT So transporter. Our simulations revealed conformational changes and key intermediate states involved in peptide translocation. We observed that the presence of the Ala-Ala peptide substrate lowers the free energy barriers associated with transition to the inward-facing state. Furthermore, we elucidated the proton transport model and analyzed the pharmacophore features of intermediate states, providing insights for rational drug design. These findings highlight the significance of substrate binding in modulating the conformational dynamics of PepT So and identify critical residues that facilitate transport.

Leveraging machine learning models for peptide-protein interaction prediction.

Peptides play a pivotal role in a wide range of biological activities through participating in up to 40% protein-protein interactions in cellular processes. They also demonstrate remarkable specificity and efficacy, making them promising candidates for drug development. However, predicting peptide-protein complexes by traditional computational approaches, such as docking and molecular dynamics simulations, still remains a challenge due to high computational cost, flexible nature of peptides, and limited structural information of peptide-protein complexes. In recent years, the surge of available biological data has given rise to the development of an increasing number of machine learning models for predicting peptide-protein interactions. These models offer efficient solutions to address the challenges associated with traditional computational approaches. Furthermore, they offer enhanced accuracy, robustness, and interpretability in their predictive outcomes. This review presents a comprehensive overview of machine learning and deep learning models that have emerged in recent years for the prediction of peptide-protein interactions.

Structure-function analysis of plant G-protein regulatory mechanisms identifies key Gα-RGS protein interactions.

Heterotrimeric GTP-binding protein alpha subunit (Gα) and its cognate regulator of G-protein signaling (RGS) protein transduce signals in eukaryotes spanning protists, amoeba, animals, fungi, and plants. The core catalytic mechanisms of the GTPase activity of Gα and the interaction interface with RGS for the acceleration of GTP hydrolysis seem to be conserved across these groups; however, the RGS gene is under low selective pressure in plants, resulting in its frequent loss. Our current understanding of the structural basis of Gα:RGS regulation in plants has been shaped by Arabidopsis Gα, (AtGPA1), which has a cognate RGS protein. To gain a comprehensive understanding of this regulation beyond Arabidopsis, we obtained the x-ray crystal structures of Oryza sativa Gα, which has no RGS, and Selaginella moellendorffi (a lycophyte) Gα that has low sequence similarity with AtGPA1 but has an RGS. We show that the three-dimensional structure, protein-protein interaction with RGS, and the dynamic features of these Gα are similar to AtGPA1 and metazoan Gα. Molecular dynamic simulation of the Gα-RGS interaction identifies the contacts established by specific residues of the switch regions of GTP-bound Gα, crucial for this interaction, but finds no significant difference due to specific amino acid substitutions. Together, our data provide valuable insights into the regulatory mechanisms of plant G-proteins but do not support the hypothesis of adaptive co-evolution of Gα:RGS proteins in plants.

Substrate Prediction for RiPP Biosynthetic Enzymes via Masked Language Modeling and Transfer Learning.

Ribosomally synthesized and post-translationally modified peptide (RiPP) biosynthetic enzymes often exhibit promiscuous substrate preferences that cannot be reduced to simple rules. Large language models are promising tools for predicting such peptide fitness landscapes. However, state-of-the-art protein language models are trained on relatively few peptide sequences. A previous study comprehensively profiled the peptide substrate preferences of LazBF (a two-component serine dehydratase) and LazDEF (a three-component azole synthetase) from the lactazole biosynthetic pathway. We demonstrated that masked language modeling of LazBF substrate preferences produced language model embeddings that improved downstream classification models of both LazBF and LazDEF substrates. Similarly, masked language modelling of LazDEF substrate preferences produced embeddings that improved the performance of classification models of both LazBF and LazDEF substrates. Our results suggest that the models learned functional forms that are transferable between distinct enzymatic transformations that act within the same biosynthetic pathway. Our transfer learning method improved performance and data efficiency in data-scarce scenarios. We then fine-tuned models on each data set and showed that the fine-tuned models provided interpretable insight that we anticipate will facilitate the design of substrate libraries that are compatible with desired RiPP biosynthetic pathways.

Functional regulation of aquaporin dynamics by lipid bilayer composition.

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics are examined. We demonstrate that SoPIP2;1's structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

Binding Position Dependent Modulation of Smoothened Activity by Cyclopamine.

Cyclopamine is a natural alkaloid that is known to act as an agonist when it binds to the Cysteine Rich Domain (CRD) of the Smoothened receptor and as an antagonist when it binds to the Transmembrane Domain (TMD). To study the effect of cyclopamine binding to each binding site experimentally, mutations in the other site are required. Hence, simulations are critical for understanding the WT activity due to binding at different sites. Additionally, there is a possibility that cyclopamine could bind to both sites simultaneously especially at high concentration, the implications of which remain unknown. We performed three independent sets of simulations to observe the receptor activation with cyclopamine bound to each site independently (CRD, TMD) and bound to both sites simultaneously. Using multi-milliseconds long aggregate MD simulations combined with Markov state models and machine learning, we explored the dynamic behavior of cyclopamine's interactions with different domains of WT SMO. A higher population of the active state at equilibrium, a lower activation free energy barrier of ~ 2 kcal/mol, and expansion of the hydrophobic tunnel to facilitate cholesterol transport agrees with the cyclopamine's agonistic behavior when bound to the CRD of SMO. A higher population of the inactive state at equilibrium, a higher free energy barrier of ~ 4 kcal/mol and restricted the hydrophobic tunnel to impede cholesterol transport showed cyclopamine's antagonistic behavior when bound to TMD. With cyclopamine bound to both sites, there was a slightly larger inactive population at equilibrium and an increased free energy barrier (~ 3.5 kcal/mol). The tunnel was slightly larger than when solely bound to TMD, and showed a balance between agonism and antagonism with respect to residue movements exhibiting an overall weak antagonistic effect.

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.

The Experimentalist's Guide to Machine Learning for Small Molecule Design.

Initially part of the field of artificial intelligence, machine learning (ML) has become a booming research area since branching out into its own field in the 1990s. After three decades of refinement, ML algorithms have accelerated scientific developments across a variety of research topics. The field of small molecule design is no exception, and an increasing number of researchers are applying ML techniques in their pursuit of discovering, generating, and optimizing small molecule compounds. The goal of this review is to provide simple, yet descriptive, explanations of some of the most commonly utilized ML algorithms in the field of small molecule design along with those that are highly applicable to an experimentally focused audience. The algorithms discussed here span across three ML paradigms: supervised learning, unsupervised learning, and ensemble methods. Examples from the published literature will be provided for each algorithm. Some common pitfalls of applying ML to biological and chemical data sets will also be explained, alongside a brief summary of a few more advanced paradigms, including reinforcement learning and semi-supervised learning.

Leveraging Machine Learning Models for Peptide-Protein Interaction Prediction.

Peptides play a pivotal role in a wide range of biological activities through participating in up to 40% protein-protein interactions in cellular processes. They also demonstrate remarkable specificity and efficacy, making them promising candidates for drug development. However, predicting peptide-protein complexes by traditional computational approaches, such as Docking and Molecular Dynamics simulations, still remains a challenge due to high computational cost, flexible nature of peptides, and limited structural information of peptide-protein complexes. In recent years, the surge of available biological data has given rise to the development of an increasing number of machine learning models for predicting peptide-protein interactions. These models offer efficient solutions to address the challenges associated with traditional computational approaches. Furthermore, they offer enhanced accuracy, robustness, and interpretability in their predictive outcomes. This review presents a comprehensive overview of machine learning and deep learning models that have emerged in recent years for the prediction of peptide-protein interactions.

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.

Adaptive Sampling Methods for Molecular Dynamics in the Era of Machine Learning.

Molecular dynamics (MD) simulations are fundamental computational tools for the study of proteins and their free energy landscapes. However, sampling protein conformational changes through MD simulations is challenging due to the relatively long time scales of these processes. Many enhanced sampling approaches have emerged to tackle this problem, including biased sampling and path-sampling methods. In this Perspective, we focus on adaptive sampling algorithms. These techniques differ from other approaches because the thermodynamic ensemble is preserved and the sampling is enhanced solely by restarting MD trajectories at particularly chosen seeds rather than introducing biasing forces. We begin our treatment with an overview of theoretically transparent methods, where we discuss principles and guidelines for adaptive sampling. Then, we present a brief summary of select methods that have been applied to realistic systems in the past. Finally, we discuss recent advances in adaptive sampling methodology powered by deep learning techniques, as well as their shortcomings.

Sequence basis for selectivity of ephrin-B2 ligand for Eph receptors and pathogenic henipavirus G glycoproteins.

IMPORTANCE: Ephrin-B2 (EFNB2) is a ligand for six Eph receptors in humans and regulates multiple cell developmental and signaling processes. It also functions as the cell entry receptor for Nipah virus and Hendra virus, zoonotic viruses that can cause respiratory and/or neurological symptoms in humans with high mortality. Here, we investigate the sequence basis of EFNB2 specificity for binding the Nipah virus attachment G glycoprotein over Eph receptors. We then use this information to engineer EFNB2 as a soluble decoy receptor that specifically binds the attachment glycoproteins of the Nipah virus and other related henipaviruses to neutralize infection. These findings further mechanistic understanding of protein selectivity and may facilitate the development of diagnostics or therapeutics against henipavirus infection.

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.

Functional Regulation of Aquaporin Dynamics by Lipid Bilayer Composition.

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics were examined. We demonstrate that SoPIP2;1s structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

Sequence controlled secondary structure is important for the site-selectivity of lanthipeptide cyclization.

Lanthipeptides are ribosomally synthesized and post-translationally modified peptides that are generated from precursor peptides through a dehydration and cyclization process. ProcM, a class II lanthipeptide synthetase, demonstrates high substrate tolerance. It is enigmatic that a single enzyme can catalyze the cyclization process of many substrates with high fidelity. Previous studies suggested that the site-selectivity of lanthionine formation is determined by substrate sequence rather than by the enzyme. However, exactly how substrate sequence contributes to site-selective lanthipeptide biosynthesis is not clear. In this study, we performed molecular dynamic simulations for ProcA3.3 variants to explore how the predicted solution structure of the substrate without enzyme correlates to the final product formation. Our simulation results support a model in which the secondary structure of the core peptide is important for the final product's ring pattern for the substrates investigated. We also demonstrate that the dehydration step in the biosynthesis pathway does not influence the site-selectivity of ring formation. In addition, we performed simulation for ProcA1.1 and 2.8, which are well-suited candidates to investigate the connection between order of ring formation and solution structure. Simulation results indicate that in both cases, C-terminal ring formation is more likely which was supported by experimental results. Our findings indicate that the substrate sequence and its solution structure can be used to predict the site-selectivity and order of ring formation, and that secondary structure is a crucial factor influencing the site-selectivity. Taken together, these findings will facilitate our understanding of the lanthipeptide biosynthetic mechanism and accelerate bioengineering efforts for lanthipeptide-derived products.