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.

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.

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.

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.

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.

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins.

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.

Orthology-based analysis helps map evolutionary diversification and predict substrate class use of BAHD acyltransferases.

Large enzyme families catalyze metabolic diversification by virtue of their ability to use diverse chemical scaffolds. How enzyme families attain such functional diversity is not clear. Furthermore, duplication and promiscuity in such enzyme families limits their functional prediction, which has produced a burgeoning set of incompletely annotated genes in plant genomes. Here, we address these challenges using BAHD acyltransferases as a model. This fast-evolving family expanded drastically in land plants, increasing from one to five copies in algae to approximately 100 copies in diploid angiosperm genomes. Compilation of >160 published activities helped visualize the chemical space occupied by this family and define eight different classes based on structural similarities between acceptor substrates. Using orthologous groups (OGs) across 52 sequenced plant genomes, we developed a method to predict BAHD acceptor substrate class utilization as well as origins of individual BAHD OGs in plant evolution. This method was validated using six novel and 28 previously characterized enzymes and helped improve putative substrate class predictions for BAHDs in the tomato genome. Our results also revealed that while cuticular wax and lignin biosynthetic activities were more ancient, anthocyanin acylation activity was fixed in BAHDs later near the origin of angiosperms. The OG-based analysis enabled identification of signature motifs in anthocyanin-acylating BAHDs, whose importance was validated via molecular dynamic simulations, site-directed mutagenesis and kinetic assays. Our results not only describe how BAHDs contributed to evolution of multiple chemical phenotypes in the plant world but also propose a biocuration-enabled approach for improved functional annotation of plant enzyme families.

Impact of Increased Membrane Realism on Conformational Sampling of Proteins.

The realism and accuracy of lipid bilayer simulations through molecular dynamics (MD) are heavily dependent on the lipid composition. While the field is pushing toward implementing more heterogeneous and realistic membrane compositions, a lack of high-resolution lipidomic data prevents some membrane protein systems from being modeled with the highest level of realism. Given the additional diversity of real-world cellular membranes and protein-lipid interactions, it is still not fully understood how altering membrane complexity affects modeled membrane protein functions or if it matters over long-timescale simulations. This is especially true for organisms whose membrane environments have little to no computational study, such as the plant plasma membrane. Tackling these issues in tandem, a generalized, realistic, and asymmetric plant plasma membrane with more than 10 different lipid species is constructed herein. Classical MD simulations of pure membrane constructs were performed to evaluate how altering the compositional complexity of the membrane impacted the plant membrane properties. The apo form of a plant sugar transporter, OsSWEET2b, was inserted into membrane models where lipid diversity was calculated in either a size-dependent or size-independent manner. An adaptive sampling simulation regime validated by Markov-state models was performed to capture the gating dynamics of OsSWEET2b in each of these membrane constructs. In comparison to previous OsSWEET2b simulations performed in a pure POPC bilayer, we confirm that simulations performed within a native-like membrane composition alter the stabilization of apo OsSWEET2b conformational states by ∼1 kcal/mol. The free-energy barriers of intermediate conformational states decrease when realistic membrane complexity is simplified, albeit roughly within sampling error, suggesting that protein-specific responses to membranes differ due to altered packing caused by compositional fluctuations. This work serves as a case study where a more realistic bilayer composition makes unbiased conformational sampling easier to achieve than with simplified bilayers.

CYP2J2 Molecular Recognition: A New Axis for Therapeutic Design.

Cytochrome P450 (CYP) epoxygenases are a special subset of heme-containing CYP enzymes capable of performing the epoxidation of polyunsaturated fatty acids (PUFA) and the metabolism of xenobiotics. This dual functionality positions epoxygenases along a metabolic crossroad. Therefore, structure-function studies are critical for understanding their role in bioactive oxy-lipid synthesis, drug-PUFA interactions, and for designing therapeutics that directly target the epoxygenases. To better exploit CYP epoxygenases as therapeutic targets, there is a need for improved understanding of epoxygenase structure-function. Of the characterized epoxygenases, human CYP2J2 stands out as a potential target because of its role in cardiovascular physiology. In this review, the early research on the discovery and activity of epoxygenases is contextualized to more recent advances in CYP epoxygenase enzymology with respect to PUFA and drug metabolism. Additionally, this review employs CYP2J2 epoxygenase as a model system to highlight both the seminal works and recent advances in epoxygenase enzymology. Herein we cover CYP2J2's interactions with PUFAs and xenobiotics, its tissue-specific physiological roles in diseased states, and its structural features that enable epoxygenase function. Additionally, the enumeration of research on CYP2J2 identifies the future needs for the molecular characterization of CYP2J2 to enable a new axis of therapeutic design.

Extensive CRISPR RNA modification reveals chemical compatibility and structure-activity relationships for Cas9 biochemical activity.

CRISPR (clustered regularly interspaced short palindromic repeat) endonucleases are at the forefront of biotechnology, synthetic biology and gene editing. Methods for controlling enzyme properties promise to improve existing applications and enable new technologies. CRISPR enzymes rely on RNA cofactors to guide catalysis. Therefore, chemical modification of the guide RNA can be used to characterize structure-activity relationships within CRISPR ribonucleoprotein (RNP) enzymes and identify compatible chemistries for controlling activity. Here, we introduce chemical modifications to the sugar-phosphate backbone of Streptococcus pyogenes Cas9 CRISPR RNA (crRNA) to probe chemical and structural requirements. Ribose sugars that promoted or accommodated A-form helical architecture in and around the crRNA 'seed' region were tolerated best. A wider range of modifications were acceptable outside of the seed, especially D-2'-deoxyribose, and we exploited this property to facilitate exploration of greater chemical diversity within the seed. 2'-fluoro was the most compatible modification whereas bulkier O-methyl sugar modifications were less tolerated. Activity trends could be rationalized for selected crRNAs using RNP stability and DNA target binding experiments. Cas9 activity in vitro tolerated most chemical modifications at predicted 2'-hydroxyl contact positions, whereas editing activity in cells was much less tolerant. The biochemical principles of chemical modification identified here will guide CRISPR-Cas9 engineering and enable new or improved applications.

Cross-talk of cannabinoid and endocannabinoid metabolism is mediated via human cardiac CYP2J2.

Phytocannabinoids have well-known cardiovascular implications. For instance, Δ9-tetrahydrocannabinol (Δ9-THC), the principal component of cannabis, induces tachycardia in humans. In order to understand the impact of phytocannabinoids on human cardiovascular health, there is a need to study the metabolism of phytocannabinoids by cardiac cytochromes p450 (CYPs). CYP2J2, the primary CYP of cardiomyocytes, is responsible for the metabolism of the endocannabinoid, anandamide (AEA), into cardioprotective epoxides (EET-EAs). Herein, we have investigated the kinetics of the direct metabolism of six phytocannabinoids (Δ9-THC, Δ8-tetrahydrocannabinol, cannabinol, cannabidiol, cannabigerol, and cannabichromene) by CYP2J2. CYP2J2 mainly produces 1'/1″-OH metabolites of these phytocannabinoids. These phytocannabinoids are metabolized with greater catalytic efficiency compared to the metabolism of AEA by CYP2J2. We have also determined that the phytocannabinoids are potent inhibitors of CYP2J2-mediated AEA metabolism, with Δ9-THC being the strongest inhibitor. Most of the inhibition of CYP2J2 by the phytocannabinoids follow a noncompetitive inhibition model, and therefore dramatically reduce the formation of EET-EAs by CYP2J2. Taken together, these data demonstrate that phytocannabinoids are directly metabolized by CYP2J2 and inhibit human cardiac CYP2J2, leading to a reduction in the formation of cardioprotective EET-EAs.