Vicia faba SV channel VfTPC1 is a hyperexcitable variant of plant vacuole Two Pore Channels.

To fire action-potential-like electrical signals, the vacuole membrane requires the two-pore channel TPC1, formerly called SV channel. The TPC1/SV channel functions as a depolarization-stimulated, non-selective cation channel that is inhibited by luminal Ca2+. In our search for species-dependent functional TPC1 channel variants with different luminal Ca2+ sensitivity, we found in total three acidic residues present in Ca2+ sensor sites 2 and 3 of the Ca2+-sensitive AtTPC1 channel from Arabidopsis thaliana that were neutral in its Vicia faba ortholog and also in those of many other Fabaceae. When expressed in the Arabidopsis AtTPC1-loss-of-function background, wild-type VfTPC1 was hypersensitive to vacuole depolarization and only weakly sensitive to blocking luminal Ca2+. When AtTPC1 was mutated for these VfTPC1-homologous polymorphic residues, two neutral substitutions in Ca2+ sensor site 3 alone were already sufficient for the Arabidopsis At-VfTPC1 channel mutant to gain VfTPC1-like voltage and luminal Ca2+ sensitivity that together rendered vacuoles hyperexcitable. Thus, natural TPC1 channel variants exist in plant families which may fine-tune vacuole excitability and adapt it to environmental settings of the particular ecological niche.

Merging Signaling with Structure: Functions and Mechanisms of Plant Glutamate Receptor Ion Channels.

Plant glutamate receptor-like (GLR) genes encode ion channels with demonstrated roles in electrical and calcium (Ca2+) signaling. The expansion of the GLR family along the lineage of land plants, culminating in the appearance of a multiclade system among flowering plants, has been a topic of interest since their discovery nearly 25 years ago. GLRs are involved in many physiological processes, from wound signaling to transcriptional regulation to sexual reproduction. Emerging evidence supports the notion that their fundamental functions are conserved among different groups of plants as well. In this review, we update the physiological and genetic evidence for GLRs, establishing their role in signaling and cell-cell communication. Special emphasis is given to the recent discussion of GLRs' atomic structures. Along with functional assays, a structural view of GLRs' molecular organization presents a window for novel hypotheses regarding the molecular mechanisms underpinning signaling associated with the ionic fluxes that GLRs regulate. Newly uncovered transcriptional regulations associated with GLRs-which propose the involvement of genes from all clades ofArabidopsis thaliana in ways not previously observed-are discussed in the context of the broader impacts of GLR activity. We posit that the functions of GLRs in plant biology are probably much broader than anticipated, but describing their widespread involvement will only be possible with (a) a comprehensive understanding of the channel's properties at the molecular and structural levels, including protein-protein interactions, and (b) the design of new genetic approaches to explore stress and pathogen responses where precise transcriptional control may result in more precise testable hypotheses to overcome their apparent functional redundancies.

The Surprising Dynamics of Electrochemical Coupling at Membrane Sandwiches in Plants.

Transport processes across membranes play central roles in any biological system. They are essential for homeostasis, cell nutrition, and signaling. Fluxes across membranes are governed by fundamental thermodynamic rules and are influenced by electrical potentials and concentration gradients. Transmembrane transport processes have been largely studied on single membranes. However, several important cellular or subcellular structures consist of two closely spaced membranes that form a membrane sandwich. Such a dual membrane structure results in remarkable properties for the transport processes that are not present in isolated membranes. At the core of membrane sandwich properties, a small intermembrane volume is responsible for efficient coupling between the transport systems at the two otherwise independent membranes. Here, we present the physicochemical principles of transport coupling at two adjacent membranes and illustrate this concept with three examples. In the supplementary material, we provide animated PowerPoint presentations that visualize the relationships. They could be used for teaching purposes, as has already been completed successfully at the University of Talca.

TPC1-Type Channels in Physcomitrium patens: Interaction between EF-Hands and Ca2.

Two-pore channels (TPCs) are members of the superfamily of ligand-gated and voltage-sensitive ion channels in the membranes of intracellular organelles of eukaryotic cells. The evolution of ordinary plant TPC1 essentially followed a very conservative pattern, with no changes in the characteristic structural footprints of these channels, such as the cytosolic and luminal regions involved in Ca2+ sensing. In contrast, the genomes of mosses and liverworts encode also TPC1-like channels with larger variations at these sites (TPC1b channels). In the genome of the model plant Physcomitrium patens we identified nine non-redundant sequences belonging to the TPC1 channel family, two ordinary TPC1-type, and seven TPC1b-type channels. The latter show variations in critical amino acids in their EF-hands essential for Ca2+ sensing. To investigate the impact of these differences between TPC1 and TPC1b channels, we generated structural models of the EF-hands of PpTPC1 and PpTPC1b channels. These models were used in molecular dynamics simulations to determine the frequency with which calcium ions were present in a coordination site and also to estimate the average distance of the ions from the center of this site. Our analyses indicate that the EF-hand domains of PpTPC1b-type channels have a lower capacity to coordinate calcium ions compared with those of common TPC1-like channels.

Antifungal Activity and In Silico Studies on 2-Acylated Benzo- and Naphthohydroquinones.

The high rates of morbidity and mortality due to fungal infections are associated with a limited antifungal arsenal and the high toxicity of drugs. Therefore, the identification of novel drug targets is challenging due to the several resemblances between fungal and human cells. Here, we report the in vitro antifungal evaluation of two acylphenols series, namely 2-acyl-1,4-benzo- and 2-acyl-1,4-naphthohydroquinones. The antifungal properties were assessed on diverse Candida and filamentous fungi strains through the halo of inhibition (HOI) and minimal inhibitory concentration (MIC). The antifungal activities of 2-acyl-1,4-benzohydroquinone derivatives were higher than those of the 2-acyl-1,4-naphthohydroquinone analogues. The evaluation indicates that 2-octanoylbenzohydroquinone 4 is the most active member of the 2-acylbenzohydroquinone series, with MIC values ranging from 2 to 16 µg/mL. In some fungal strains (i.e., Candida krusei and Rhizopus oryzae), such MIC values of compound 4 (2 and 4 µg/mL) were comparable to that obtained by amphotericin B (1 µg/mL). The compound 4 was evaluated for its antioxidant activity by means of FRAP, ABTS and DPPH assays, showing moderate activity as compared to standard antioxidants. Molecular docking studies of compound 4 and ADMET predictions make this compound a potential candidate for topical pharmacological use. The results obtained using the most active acylbenzohydroquinones are promising because some evaluated Candida strains are known to have decreased sensitivity to standard antifungal treatments.

Computational Analyses of the AtTPC1 (Arabidopsis Two-Pore Channel 1) Permeation Pathway.

Two Pore Channels (TPCs) are cation-selective voltage- and ligand-gated ion channels in membranes of intracellular organelles of eukaryotic cells. In plants, the TPC1 subtype forms the slowly activating vacuolar (SV) channel, the most dominant ion channel in the vacuolar membrane. Controversial reports about the permeability properties of plant SV channels fueled speculations about the physiological roles of this channel type. TPC1 is thought to have high Ca2+ permeability, a conclusion derived from relative permeability analyses using the Goldman-Hodgkin-Katz (GHK) equation. Here, we investigated in computational analyses the properties of the permeation pathway of TPC1 from Arabidopsis thaliana. Using the crystal structure of AtTPC1, protein modeling, molecular dynamics (MD) simulations, and free energy calculations, we identified a free energy minimum for Ca2+, but not for K+, at the luminal side next to the selectivity filter. Residues D269 and E637 coordinate in particular Ca2+ as demonstrated in in silico mutagenesis experiments. Such a Ca2+-specific coordination site in the pore explains contradicting data for the relative Ca2+/K+ permeability and strongly suggests that the Ca2+ permeability of SV channels is largely overestimated from relative permeability analyses. This conclusion was further supported by in silico electrophysiological studies showing a remarkable permeation of K+ but not Ca2+ through the open channel.

Acidosis-induced activation of anion channel SLAH3 in the flooding-related stress response of Arabidopsis.

Plants, as sessile organisms, gained the ability to sense and respond to biotic and abiotic stressors to survive severe changes in their environments. The change in our climate comes with extreme dry periods but also episodes of flooding. The latter stress condition causes anaerobiosis-triggered cytosolic acidosis and impairs plant function. The molecular mechanism that enables plant cells to sense acidity and convey this signal via membrane depolarization was previously unknown. Here, we show that acidosis-induced anion efflux from Arabidopsis (Arabidopsis thaliana) roots is dependent on the S-type anion channel AtSLAH3. Heterologous expression of SLAH3 in Xenopus oocytes revealed that the anion channel is directly activated by a small, physiological drop in cytosolic pH. Acidosis-triggered activation of SLAH3 is mediated by protonation of histidine 330 and 454. Super-resolution microscopy analysis showed that the increase in cellular proton concentration switches SLAH3 from an electrically silent channel dimer into its active monomeric form. Our results show that, upon acidification, protons directly switch SLAH3 to its open configuration, bypassing kinase-dependent activation. Moreover, under flooding conditions, the stress response of Arabidopsis wild-type (WT) plants was significantly higher compared to SLAH3 loss-of-function mutants. Our genetic evidence of SLAH3 pH sensor function may guide the development of crop varieties with improved stress tolerance.

Transforming Non-Selective Angiotensin-Converting Enzyme Inhibitors in C- and N-domain Selective Inhibitors by Using Computational Tools.

The two-domain dipeptidylcarboxypeptidase Angiotensin-I-converting enzyme (EC 3.4.15.1; ACE) plays an important physiological role in blood pressure regulation via the reninangiotensin and kallikrein-kinin systems by converting angiotensin I to the potent vasoconstrictor angiotensin II, and by cleaving a number of other substrates including the vasodilator bradykinin and the anti-inflammatory peptide N-acetyl-SDKP. Therefore, the design of ACE inhibitors is within the priorities of modern medical sciences for treating hypertension, heart failures, myocardial infarction, and other related diseases. Despite the success of ACE inhibitors for the treatment of hypertension and congestive heart failure, they have some adverse effects, which could be attenuated by selective domain inhibition. Crystal structures of both ACE domains (nACE and cACE) reported over the last decades could facilitate the rational drug design of selective inhibitors. In this review, we refer to the history of the discovery of ACE inhibitors, which has been strongly related to the development of molecular modeling methods. We stated that the design of novel selective ACE inhibitors is a challenge for current researchers which requires a thorough understanding of the structure of both ACE domains and the help of molecular modeling methodologies. Finally, we performed a theoretical design of potential selective derivatives of trandolaprilat, a drug approved to treat critical conditions of hypertension, to illustrate how to use molecular modeling methods such as de novo design, docking, Molecular Dynamics (MD) simulations, and free energy calculations for creating novel potential drugs with specific interactions inside nACE and cACE binding sites.

Docking, Interaction Fingerprint, and Three-Dimensional Quantitative Structure-Activity Relationship (3D-QSAR) of Sigma1 Receptor Ligands, Analogs of the Neuroprotective Agent RC-33.

The human Sigma1 receptor (S1R), which has been identified as a target with an important role in neuropsychological disorders, was first crystallized 3 years ago. Since S1R structure has no relation with another previous crystallized structures, the presence of the new crystal is an important hallmark for the design of agonists and antagonists against this important target. Some years ago, our group identified RC-33, a potent and selective S1R agonist, endowed with neuroprotective properties. In this work, drawing on new structural information, we studied the interactions of RC-33 and its analogs with the S1R binding site by using computational methods such as docking, interaction fingerprints, and receptor-guided alignment three dimensional quantitative structure-activity relationship (3D-QSAR). We found that RC-33 and its analogs adopted similar orientations within S1R binding site, with high similitude with orientations of the crystallized ligands; such information was used for identifying the residues involved in chemical interactions with ligands. Furthermore, the structure-activity relationship of the studied ligands was adequately described considering classical QSAR tests. All relevant aspects of the interactions between the studied compounds and S1R were covered here, through descriptions of orientations, binding interactions, and features that influence differential affinities. In this sense, the present results could be useful in the future design of novel S1R modulators.

Nutrient exchange in arbuscular mycorrhizal symbiosis from a thermodynamic point of view.

To obtain insights into the dynamics of nutrient exchange in arbuscular mycorrhizal (AM) symbiosis, we modelled mathematically the two-membrane system at the plant-fungus interface and simulated its dynamics. In computational cell biology experiments, the full range of nutrient transport pathways was tested for their ability to exchange phosphorus (P)/carbon (C)/nitrogen (N) sources. As a result, we obtained a thermodynamically justified, independent and comprehensive model of the dynamics of the nutrient exchange at the plant-fungus contact zone. The predicted optimal transporter network coincides with the transporter set independently confirmed in wet-laboratory experiments previously, indicating that all essential transporter types have been discovered. The thermodynamic analyses suggest that phosphate is released from the fungus via proton-coupled phosphate transporters rather than anion channels. Optimal transport pathways, such as cation channels or proton-coupled symporters, shuttle nutrients together with a positive charge across the membranes. Only in exceptional cases does electroneutral transport via diffusion facilitators appear to be plausible. The thermodynamic models presented here can be generalized and adapted to other forms of mycorrhiza and open the door for future studies combining wet-laboratory experiments with computational simulations to obtain a deeper understanding of the investigated phenomena.

Insights into the Structural Requirements of 2(S)-Amino-6-Boronohexanoic Acid Derivatives as Arginase I Inhibitors: 3D-QSAR, Docking, and Interaction Fingerprint Studies.

Human arginase I (hARGI) is an important enzyme involved in the urea cycle; its overexpression has been associated to cardiovascular and cerebrovascular diseases. In the last years, several congeneric sets of hARGI inhibitors have been reported with possible beneficial roles for the cardiovascular system. At the same time, crystallographic data have been reported including hARGI⁻inhibitor complexes, which can be considered for the design of novel inhibitors. In this work, the structure⁻activity relationship (SAR) of Cα substituted 2(S)-amino-6-boronohexanoic acid (ABH) derivatives as hARGI inhibitors was studied by using a three-dimensional quantitative structure⁻activity relationships (3D-QSAR) method. The predictivity of the obtained 3D-QSAR model was demonstrated by using internal and external validation experiments. The best model revealed that the differential hARGI inhibitory activities of the ABH derivatives can be described by using steric and electrostatic fields; the local effects of these fields in the activity are presented. In addition, binding modes of the above-mentioned compounds inside the hARGI binding site were obtained by using molecular docking. It was found that ABH derivatives adopted the same orientation reported for ABH within the hARGI active site, with the substituents at Cα exposed to the solvent with interactions with residues at the entrance of the binding site. The hARGI residues involved in chemical interactions with inhibitors were identified by using an interaction fingerprints (IFPs) analysis.

Mycobacterium tuberculosis serine/threonine protein kinases: structural information for the design of their specific ATP-competitive inhibitors.

In the last decades, human protein kinases (PKs) have been relevant as targets in the development of novel therapies against many diseases, but the study of Mycobacterium tuberculosis PKs (MTPKs) involved in tuberculosis pathogenesis began much later and has not yet reached an advanced stage of development. To increase knowledge of these enzymes, in this work we studied the structural features of MTPKs, with focus on their ATP-binding sites and their interactions with inhibitors. PknA, PknB, and PknG are the most studied MTPKs, which were previously crystallized; ATP-competitive inhibitors have been designed against them in the last decade. In the current work, reported PknA, PknB, and PknG inhibitors were extracted from literature and their orientations inside the ATP-binding site were proposed by using docking method. With this information, interaction fingerprints were elaborated, which reveal the more relevant residues for establishing chemical interactions with inhibitors. The non-crystallized MTPKs PknD, PknF, PknH, PknJ, PknK, and PknL were also studied; their three-dimensional structural models were developed by using homology modeling. The main characteristics of MTPK ATP-binding sites (the non-crystallized and crystallized MTPKs, including PknE and PknI) were accounted; schemes of the main polar and nonpolar groups inside their ATP-binding sites were constructed, which are suitable for a major understanding of these proteins as antituberculotic targets. These schemes could be used for establishing comparisons between MTPKs and human PKs in order to increase selectivity of MTPK inhibitors. As a key tool for guiding medicinal chemists interested in the design of novel MTPK inhibitors, our work provides a map of the structural elements relevant for the design of more selective ATP-competitive MTPK inhibitors.

Folding and Lipid Composition Determine Membrane Interaction of the Disordered Protein COR15A.

Plants from temperate climates, such as the model plant Arabidopsis thaliana, are challenged with seasonal low temperatures that lead to increased freezing tolerance in fall in a process termed cold acclimation. Among other adaptations, this involves the accumulation of cold-regulated (COR) proteins, such as the intrinsically disordered chloroplast-localized protein COR15A. Together with its close homolog COR15B, it stabilizes chloroplast membranes during freezing. COR15A folds into amphipathic α-helices in the presence of high concentrations of low-molecular-mass crowders or upon dehydration. Under these conditions, the (partially) folded protein binds peripherally to membranes. In our study, we have used coarse-grained molecular dynamics simulations to elucidate the details of COR15A-membrane binding and its effects on membrane structure and dynamics. Simulation results indicate that at least partial folding of COR15A and the presence of highly unsaturated galactolipids in the membranes are necessary for efficient membrane binding. The bound protein is stabilized on the membrane by interactions of charged and polar amino acids with galactolipid headgroups and by interactions of hydrophobic amino acids with the upper part of the fatty acyl chains. Experimentally, the presence of liposomes made from a mixture of lipids mimicking chloroplast membranes induces additional folding in COR15A under conditions of partial dehydration, in agreement with the simulation results.

Energetic differences between non-domain-swapped and domain-swapped chain connectivities in the K2P potassium channel TRAAK.

Two-pore domain (K2P) channels are twofold symmetric K+ channels which control cell excitability by enabling the leak of potassium ions from cells in response to physicochemical stimuli. Crystallization of K2P channels revealed the presence of several structural features, which include an external cap. In the available crystallographic structures, the cap is present as non-domain-swapped (NDS) and domain-swapped (DS) chain conformations, where DS chain conformation exchanges two opposing outer helices 180° around the channel. In this work, energy differences between the residues located at the highest point of the cap in NDS and DS conformations were evaluated for TRAAK, a K2P channel that was crystallized in both conformations. Results indicated a preference for DS conformation, but this result is not extensible to TASK K2P channels.

Direct and Auger Electron-Induced, Single- and Double-Strand Breaks on Plasmid DNA Caused by 99mTc-Labeled Pyrene Derivatives and the Effect of Bonding Distance.

It is evident that 99mTc causes radical-mediated DNA damage due to Auger electrons, which were emitted simultaneously with the known γ-emission of 99mTc. We have synthesized a series of new 99mTc-labeled pyrene derivatives with varied distances between the pyrene moiety and the radionuclide. The pyrene motif is a common DNA intercalator and allowed us to test the influence of the radionuclide distance on damages of the DNA helix. In general, pUC 19 plasmid DNA enables the investigation of the unprotected interactions between the radiotracers and DNA that results in single-strand breaks (SSB) or double-strand breaks (DSB). The resulting DNA fragments were separated by gel electrophoresis and quantified by fluorescent staining. Direct DNA damage and radical-induced indirect DNA damage by radiolysis products of water were evaluated in the presence or absence of the radical scavenger DMSO. We demonstrated that Auger electrons directly induced both SSB and DSB in high efficiency when 99mTc was tightly bound to the plasmid DNA and this damage could not be completely prevented by DMSO, a free radical scavenger. For the first time, we were able to minimize this effect by increasing the carbon chain lengths between the pyrene moiety and the 99mTc nuclide. However, a critical distance between the 99mTc atom and the DNA helix could not be determined due to the significantly lowered DSB generation resulting from the interaction which is dependent on the type of the 99mTc binding motif. The effect of variable DNA damage caused by the different chain length between the pyrene residue and the Tc-core as well as the possible conformations of the applied Tc-complexes was supplemented with molecular dynamics (MD) calculations. The effectiveness of the DNA-binding 99mTc-labeled pyrene derivatives was demonstrated by comparison to non-DNA-binding 99mTcO4-, since nearly all DNA damage caused by 99mTcO4- was prevented by incubating with DMSO.

Flavonoids as CDK1 Inhibitors: Insights in Their Binding Orientations and Structure-Activity Relationship.

In the last years, the interactions of flavonoids with protein kinases (PKs) have been described by using crystallographic experiments. Interestingly, different orientations have been found for one flavonoid inside different PKs and different chemical substitutions lead to different orientations of the flavonoid scaffold inside one PK. Accordingly, orientation predictions of novel analogues could help to the design of flavonoids with high PK inhibitory activities. With this in mind, we studied the binding modes of 37 flavonoids (flavones and chalcones) inside the cyclin-dependent PK CDK1 using docking experiments. We found that the compounds under study adopted two different orientations into the active site of CDK1 (orientations I and II in the manuscript). In addition, quantitative structure-activity relationship (QSAR) models using CoMFA and CoMSIA methodologies were constructed to explain the trend of the CDK1 inhibitory activities for the studied flavonoids. Template-based and docking-based alignments were used. Models developed starting from docking-based alignment were applied for describing the whole dataset and compounds with orientation I. Adequate R2 and Q2 values were obtained by each method; interestingly, only hydrophobic and hydrogen bond donor fields describe the differential potency of the flavonoids as CDK1 inhibitors for both defined alignments and subsets. Our current application of docking and QSAR together reveals important elements to be drawn for the design of novel flavonoids with increased PK inhibitory activities.

Molecular dynamics simulations and CD spectroscopy reveal hydration-induced unfolding of the intrinsically disordered LEA proteins COR15A and COR15B from Arabidopsis thaliana.

The LEA (late embryogenesis abundant) proteins COR15A and COR15B from Arabidopsis thaliana are intrinsically disordered under fully hydrated conditions, but obtain α-helical structure during dehydration, which is reversible upon rehydration. To understand this unusual structural transition, both proteins were investigated by circular dichroism (CD) and molecular dynamics (MD) approaches. MD simulations showed unfolding of the proteins in water, in agreement with CD data obtained with both HIS-tagged and untagged recombinant proteins. Mainly intramolecular hydrogen bonds (H-bonds) formed by the protein backbone were replaced by H-bonds with water molecules. As COR15 proteins function in vivo as protectants in leaves partially dehydrated by freezing, unfolding was further assessed under crowded conditions. Glycerol reduced (40%) or prevented (100%) unfolding during MD simulations, in agreement with CD spectroscopy results. H-bonding analysis indicated that preferential exclusion of glycerol from the protein backbone increased stability of the folded state.

Structural and Affinity Determinants in the Interaction between Alcohol Acyltransferase from F. x ananassa and Several Alcohol Substrates: A Computational Study.

Aroma and flavor are important factors of fruit quality and consumer preference. The specific pattern of aroma is generated during ripening by the accumulation of volatiles compounds, which are mainly esters. Alcohol acyltransferase (AAT) (EC 2.3.1.84) catalyzes the esterification reaction of aliphatic and aromatic alcohols and acyl-CoA into esters in fruits and flowers. In Fragaria x ananassa, there are different volatiles compounds that are obtained from different alcohol precursors, where octanol and hexanol are the most abundant during fruit ripening. At present, there is not structural evidence about the mechanism used by the AAT to synthesize esters. Experimental data attribute the kinetic role of this enzyme to 2 amino acidic residues in a highly conserved motif (HXXXD) that is located in the middle of the protein. With the aim to understand the molecular and energetic aspects of volatiles compound production from F. x ananassa, we first studied the binding modes of a series of alcohols, and also different acyl-CoA substrates, in a molecular model of alcohol acyltransferase from Fragaria x ananassa (SAAT) using molecular docking. Afterwards, the dynamical behavior of both substrates, docked within the SAAT binding site, was studied using routine molecular dynamics (MD) simulations. In addition, in order to correlate the experimental and theoretical data obtained in our laboratories, binding free energy calculations were performed; which previous results suggested that octanol, followed by hexanol, presented the best affinity for SAAT. Finally, and concerning the SAAT molecular reaction mechanism, it is suggested from molecular dynamics simulations that the reaction mechanism may proceed through the formation of a ternary complex, in where the Histidine residue at the HXXXD motif deprotonates the alcohol substrates. Then, a nucleophilic attack occurs from alcohol charged oxygen atom to the carbon atom at carbonyl group of the acyl CoA. This mechanism is in agreement with previous results, obtained in our group, in alcohol acyltransferase from Vasconcellea pubescens (VpAAT1).