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.

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.

Role of internal loop dynamics in antibiotic permeability of outer membrane porins.

Gram-negative bacteria pose a serious public health concern due to resistance to many antibiotics, caused by the low permeability of their outer membrane (OM). Effective antibiotics use porins in the OM to reach the interior of the cell; thus, understanding permeation properties of OM porins is instrumental to rationally develop broad-spectrum antibiotics. A functionally important feature of OM porins is undergoing open-closed transitions that modulate their transport properties. To characterize the molecular basis of these transitions, we performed an extensive set of molecular dynamics (MD) simulations of Escherichia coli OM porin OmpF. Markov-state analysis revealed that large-scale motion of an internal loop, L3, underlies the transition between energetically stable open and closed states. The conformation of L3 is controlled by H bonds between highly conserved acidic residues on the loop and basic residues on the OmpF ß-barrel. Mutation of key residues important for the loop's conformation shifts the equilibrium between open and closed states and regulates translocation of permeants (ions and antibiotics), as observed in the simulations and validated by our whole-cell accumulation assay. Notably, one mutant system G119D, which we find to favor the closed state, has been reported in clinically resistant bacterial strains. Overall, our accumulated ∼200 µs of simulation data (the wild type and mutants) along with experimental assays suggest the involvement of internal loop dynamics in permeability of OM porins and antibiotic resistance in Gram-negative bacteria.

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.

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.

Engineered ACE2 decoy mitigates lung injury and death induced by SARS-CoV-2 variants.

Vaccine hesitancy and emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) escaping vaccine-induced immune responses highlight the urgency for new COVID-19 therapeutics. Engineered angiotensin-converting enzyme 2 (ACE2) proteins with augmented binding affinities for SARS-CoV-2 spike (S) protein may prove to be especially efficacious against multiple variants. Using molecular dynamics simulations and functional assays, we show that three amino acid substitutions in an engineered soluble ACE2 protein markedly augmented the affinity for the S protein of the SARS-CoV-2 WA-1/2020 isolate and multiple VOCs: B.1.1.7 (Alpha), B.1.351 (Beta), P.1 (Gamma) and B.1.617.2 (Delta). In humanized K18-hACE2 mice infected with the SARS-CoV-2 WA-1/2020 or P.1 variant, prophylactic and therapeutic injections of soluble ACE22.v2.4-IgG1 prevented lung vascular injury and edema formation, essential features of CoV-2-induced SARS, and above all improved survival. These studies demonstrate broad efficacy in vivo of an engineered ACE2 decoy against SARS-CoV-2 variants in mice and point to its therapeutic potential.

Effect of histidine covalent modification on strigolactone receptor activation and selectivity.

The parasitic weed Striga has led to billions of dollars' worth of agricultural productivity loss worldwide. Striga detects host plants using compounds of the strigolactone class of phytohormones. Early steps in the strigolactone signaling pathway involve substrate binding and hydrolysis followed by a conformational change to an "active" or "closed" state, after which it associates with a MAX2-family downstream signaling partner. The structures of the inactive and active states of strigolactone receptors are known through X-ray crystallography, and the transition pathway from the inactive to active state in apo receptors has previously been characterized using molecular dynamics simulations. However, it also has been suggested that a covalent butenolide modification of the receptor on the catalytic histidine through substrate hydrolysis promotes formation of the active state. Using molecular dynamics simulations, we show that the presence of the covalent butenolide enhances activation in both AtD14, a receptor found in Arabidopsis, and ShHTL7, a receptor found in Striga, but the enhancement is ∼50 times greater in ShHTL7. We also show that several conserved interactions with the covalent butenolide modification promote transition to the active state in both AtD14 (non-parasite) and ShHTL7 (parasite). Finally, we demonstrate that the enhanced activation of ShHTL7 likely results from disruption of ShHTL7-specific histidine interactions that inhibited activation in the apo case. These results provide a possible explanation for difference in strigolactone sensitivity seen between different strigolactone-sensitive proteins and can be used to aid the design of selective modulators to control Striga parasites.

Activation mechanism of the human Smoothened receptor.

Smoothened (SMO) is a membrane protein of the class F subfamily of G protein-coupled receptors (GPCRs) and maintains homeostasis of cellular differentiation. SMO undergoes conformational change during activation, transmitting the signal across the membrane, making it amenable to bind to its intracellular signaling partner. Receptor activation has been studied at length for class A receptors, but the mechanism of class F receptor activation remains unknown. Agonists and antagonists bound to SMO at sites in the transmembrane domain (TMD) and the cysteine-rich domain have been characterized, giving a static view of the various conformations SMO adopts. Although the structures of the inactive and active SMO outline the residue-level transitions, a kinetic view of the overall activation process remains unexplored for class F receptors. We describe SMO's activation process in atomistic detail by performing 300 µs of molecular dynamics simulations and combining it with Markov state model theory. A molecular switch, conserved across class F and analogous to the activation-mediating D-R-Y motif in class A receptors, is observed to break during activation. We also show that this transition occurs in a stage-wise movement of the transmembrane helices: TM6 first, followed by TM5. To see how modulators affect SMO activity, we simulated agonist and antagonist-bound SMO. We observed that agonist-bound SMO has an expanded hydrophobic tunnel in SMO's core TMD, whereas antagonist-bound SMO shrinks this tunnel, further supporting the hypothesis that cholesterol travels through a tunnel inside Smoothened to activate it. In summary, this study elucidates the distinct activation mechanism of class F GPCRs and shows that SMO's activation process rearranges the core TMD to open a hydrophobic conduit for cholesterol transport.

Mechanistic origin of partial agonism of tetrahydrocannabinol for cannabinoid receptors.

Cannabinoid receptor 1 (CB1) is a therapeutically relevant drug target for controlling pain, obesity, and other central nervous system disorders. However, full agonists and antagonists of CB1 have been reported to cause serious side effects in patients. Therefore, partial agonists have emerged as a viable alternative as they can mitigate overstimulation and side effects. One of the key bottlenecks in the design of partial agonists, however, is the lack of understanding of the molecular mechanism of partial agonism itself. In this study, we examine two mechanistic hypotheses for the origin of partial agonism in cannabinoid receptors and predict the mechanistic basis of partial agonism exhibited by Δ9-Tetrahydrocannabinol (THC) against CB1. In particular, we inspect whether partial agonism emerges from the ability of THC to bind in both agonist and antagonist-binding poses or from its ability to only partially activate the receptor. We used extensive molecular dynamics simulations and Markov state modeling to capture the THC binding in both antagonist and agonist-binding poses in the CB1 receptor. Furthermore, we predict that binding of THC in the agonist-binding pose leads to rotation of toggle switch residues and causes partial outward movement of intracellular transmembrane helix 6 (TM6). Our simulations also suggest that the alkyl side chain of THC plays a crucial role in determining partial agonism by stabilizing the ligand in the agonist and antagonist-like poses within the pocket. Taken together, this study provides important insights into the mechanistic origin of the partial agonism of THC.

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.

Imparting Selective Fluorophilic Interactions in Redox Copolymers for the Electrochemically Mediated Capture of Short-Chain Perfluoroalkyl Substances.

With increasing regulations on per- and polyfluoroalkyl substances (PFAS) across the world, understanding the molecular level interactions that drive their binding by functional adsorbent materials is key to effective PFAS removal from water streams. With the phaseout of legacy long-chain PFAS, the emergence of short-chain PFAS has posed a significant challenge for material design due to their higher mobility and hydrophilicity and inefficient removal by conventional treatment methods. Here, we demonstrate how cooperative molecular interactions are essential to target short-chain PFAS (from C4 to C7) by tailoring structural units to enhance affinity while modulating the electrochemical control of capture and release of PFAS. We report a new class of fluorinated redox-active amine-functionalized copolymers to leverage both fluorophilic and electrostatic interactions for short-chain PFAS binding. We combine molecular dynamics (MD) simulations and electrosorption to elucidate the role of the designer functional groups in enabling affinity toward short-chain PFAS. Preferential interaction coefficients from MD simulations correlated closely with experimental trends: fluorination enhanced the overall PFAS uptake and promoted the capture of less hydrophobic short-chain PFAS (C ≤ 5), while electrostatic interactions provided by secondary amine groups were sufficient to capture PFAS with higher hydrophobicity (C ≥ 6). The addition of an induced electric field showed favorable kinetic enhancement for the shortest PFAS and increased the reversibility of release from the electrode. Integration of these copolymers with electrochemical separations showed potential for removing these contaminants at environmentally relevant conditions while eliminating the need for chemical regeneration.

The Effects of N-Linked Glycosylation on SLC6 Transporters.

Membrane transporters of the solute carrier 6 (SLC6) family mediate various physiological processes by facilitating the translocation of amino acids, neurotransmitters, and other metabolites. In the body, the activity of these transporters is tightly controlled through various post-translational modifications with implications on protein expression, stability, membrane trafficking, and dynamics. While N-linked glycosylation is a universal regulatory mechanism among eukaryotes, a consistent mechanism of how glycosylation affects the SLC6 transporter family remains elusive. It is generally believed that glycans influence transporter stability and membrane trafficking; however, the role of glycosylation on transporter dynamics remains disputable, with differing conclusions among individual transporters across the SLC6 family. In this study, we collected over 1 ms of aggregated all-atom molecular dynamics (MD) simulation data to systematically identify the impact of N-glycans on SLC6 transporter dynamics. We modeled four human SLC6 transporters, the serotonin, dopamine, glycine, and B0AT1 transporters, by first simulating all possible combinations of a glycan attached to each glycosylation site followed by investigating the effect of larger, oligo-N-linked glycans to each transporter. The simulations reveal that glycosylation does not significantly affect the transporter structure but alters the dynamics of the glycosylated extracellular loop and surrounding regions. The structural consequences of glycosylation on the loop dynamics are further emphasized with larger glycan molecules attached. However, no apparent differences in ligand stability or movement of the gating helices were observed, and as such, the simulations suggest that glycosylation does not have a profound effect on conformational dynamics associated with substrate transport.

Reconciling membrane protein simulations with experimental DEER spectroscopy data.

Spectroscopy experiments are crucial to study membrane proteins for which traditional structure determination methods still prove challenging. Double electron-electron resonance (DEER) spectroscopy experiments provide protein residue-pair distance distributions that are indicative of their conformational heterogeneity. Atomistic molecular dynamics (MD) simulations are another tool that have been proven to be vital to study the structural dynamics of membrane proteins such as to identify inward-open, occluded, and outward-open conformations of transporter membrane proteins, among other partially open or closed states of the protein. Yet, studies have reported that there is no direct consensus between the distributional data from DEER experiments and MD simulations, which has challenged validation of structures obtained from long-timescale simulations and using simulations to design experiments. Current coping strategies for comparisons rely on heuristics, such as mapping the nearest matching peaks between two ensembles or biased simulations. Here we examine the differences in residue-pair distance distributions arising due to the choice of membranes around the protein and covalent modification of a pair of residues to nitroxide spin labels in DEER experiments. Through comparing MD simulations of two proteins, PepTSo and LeuT-both of which have been characterized using DEER experiments previously-we show that the proteins' dynamics are similar despite the choice of the detergent micelle as a membrane mimetic in DEER experiments. On the other hand, covalently modified residues show slight local differences in their dynamics and a huge divergence when the oxygen atom pair distances between spin labeled residues are measured rather than protein backbone distances. Given the computational expense associated with pairwise MTSSL labeled MD simulations, we examine the use of biased simulations to explore the conformational dynamics of the spin labels only to reveal that such simulations alter the underlying protein dynamics. Our study identifies the main cause for the mismatch between DEER experiments and MD simulations and will accelerate the development of potential mitigation strategies to improve the match.

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.

The substrate import mechanism of the human serotonin transporter.

The serotonin transporter (SERT) initiates the reuptake of extracellular serotonin in the synapse to terminate neurotransmission. The cryogenic electron microscopy structures of SERT bound to ibogaine and the physiological substrate serotonin resolved in different states have provided a glimpse of the functional conformations at atomistic resolution. However, the conformational dynamics and structural transitions to intermediate states are not fully understood. Furthermore, the molecular basis of how serotonin is recognized and transported remains unclear. In this study, we performed unbiased microsecond-long simulations of the human SERT to investigate the structural dynamics to various intermediate states and elucidated the complete substrate import pathway. Using Markov state models, we characterized a sequential order of conformational-driven ion-coupled substrate binding and transport events and calculated the free energy barriers of conformation transitions associated with the import mechanism. We find that the transition from the occluded to inward-facing state is the rate-limiting step for substrate import and that the substrate decreases the free energy barriers to achieve the inward-facing state. Our study provides insights on the molecular basis of dynamics-driven ion-substrate recognition and transport of SERT that can serve as a model for other closely related neurotransmitter transporters.

Role of substrate recognition in modulating strigolactone receptor selectivity in witchweed.

Witchweed, or Striga hermonthica, is a parasitic weed that destroys billions of dollars' worth of crops globally every year. Its germination is stimulated by strigolactones exuded by its host plants. Despite high sequence, structure, and ligand-binding site conservation across different plant species, one strigolactone receptor in witchweed, ShHTL7, uniquely exhibits a picomolar EC50 for downstream signaling. Previous biochemical and structural analyses have hypothesized that this unique ligand sensitivity can be attributed to a large binding pocket volume in ShHTL7 resulting in enhanced ability to bind substrates, but additional structural details of the substrate-binding process would help explain its role in modulating the ligand selectivity. Using long-timescale molecular dynamics simulations, we demonstrate that mutations at the entrance of the binding pocket facilitate a more direct ligand-binding pathway to ShHTL7, whereas hydrophobicity at the binding pocket entrance results in a stable "anchored" state. We also demonstrate that several residues on the D-loop of AtD14 stabilize catalytically inactive conformations. Finally, we show that strigolactone selectivity is not modulated by binding pocket volume. Our results indicate that while ligand binding is not the sole modulator of strigolactone receptor selectivity, it is a significant contributing factor. These results can be used to inform the design of selective antagonists for strigolactone receptors in witchweed.

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.

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

Phytohormone abscisic acid (ABA) is essential for plant responses to biotic and abiotic stresses. Dimeric receptors are a class of PYR1/PYL/RCAR (pyrabactin resistance 1/PYR1-like/regulatory component of ABA receptors) ABA receptors that are important for various ABA responses. While extensive experimental and computational studies have investigated these receptors, it remains not fully understood how ABA leads to their activation and dissociation for interaction with downstream protein phosphatase 2C (PP2C). Here, we study the activation and the homodimeric association processes of the PYL2 receptor as well as its heterodimeric association with protein phosphatase 2C 16 (HAB1) using molecular dynamics simulations. Free energy landscapes from ∼223 µs simulations show that dimerization substantially constrains PYL2 conformational plasticity and stabilizes the inactive state, resulting in lower ABA affinity. Also, we establish the thermodynamic model for competitive binding between homodimeric PYL2 association and heterodimeric PYL2-HAB1 association in the absence and presence of ABA. Our results suggest that the binding of ABA destabilizes the PYL2 complex and further stabilizes PYL2-HAB1 association, thereby promoting PYL2 dissociation. Overall, this study explains several key aspects on the activation of dimeric ABA receptors, which provide new avenues for selective regulation of these receptors.

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.

Divalent cations promote TALE DNA-binding specificity.

Recent advances in gene editing have been enabled by programmable nucleases such as transcription activator-like effector nucleases (TALENs) and CRISPR-Cas9. However, several open questions remain regarding the molecular machinery in these systems, including fundamental search and binding behavior as well as role of off-target binding and specificity. In order to achieve efficient and specific cleavage at target sites, a high degree of target site discrimination must be demonstrated for gene editing applications. In this work, we studied the binding affinity and specificity for a series of TALE proteins under a variety of solution conditions using in vitro fluorescence methods and molecular dynamics (MD) simulations. Remarkably, we identified that TALEs demonstrate high sequence specificity only upon addition of small amounts of certain divalent cations (Mg2+, Ca2+). However, under purely monovalent salt conditions (K+, Na+), TALEs bind to specific and non-specific DNA with nearly equal affinity. Divalent cations preferentially bind to DNA over monovalent cations, which attenuates non-specific interactions between TALEs and DNA and further stabilizes specific interactions. Overall, these results uncover new mechanistic insights into the binding action of TALEs and further provide potential avenues for engineering and application of TALE- or TALEN-based systems for genome editing and regulation.