Synthetically derived BiAux modulates auxin co-receptor activity to stimulate lateral root formation.

The root system plays an essential role in plant growth and adaptation to the surrounding environment. The root clock periodically specifies lateral root prebranch sites (PBS), where a group of pericycle founder cells (FC) is primed to become lateral root founder cells and eventually give rise to lateral root primordia or lateral roots (LRs). This clock-driven organ formation process is tightly controlled by modulation of auxin content and signaling. Auxin perception entails the physical interaction of TRANSPORT INHIBITOR RESPONSE 1 (TIR1) or AUXIN SIGNALING F-BOX (AFBs) proteins with AUXIN/INDOLE-3-ACETIC ACID (Aux/IAA) repressors to form a co-receptor system. Despite the apparent simplicity, the understanding of how specific auxin co-receptors are assembled remains unclear. We identified the compound bis-methyl auxin conjugated with N-glucoside, or BiAux, in Arabidopsis (Arabidopsis thaliana) that specifically induces the formation of PBS and the emergence of LR, with a slight effect on root elongation. Docking analyses indicated that BiAux binds to F-box proteins, and we showed that BiAux function depends on TIR1 and AFB2 F-box proteins and AUXIN RESPONSE FACTOR 7 activity, which is involved in FC specification and LR formation. Finally, using a yeast (Saccharomyces cerevisiae) heterologous expression system, we showed that BiAux favors the assemblage of specific co-receptors subunits involved in LR formation and enhances AUXIN/INDOLE-3-ACETIC ACID 28 protein degradation. These results indicate that BiAux acts as an allosteric modulator of specific auxin co-receptors. Therefore, BiAux exerts a fine-tune regulation of auxin signaling aimed to the specific formation of LR among the many development processes regulated by auxin.

Intrinsic disorder in the dynamic evolution of structure, stability, and flexibility of potyviral VLP assemblies: A computational study.

An all-atom Molecular Dynamics (MD) study was applied to three viral nanoparticles (VLPs) of Turnip mosaic virus (TuMV), a potyvirus: the particles genetically functionalized with two peptides, VIP (human vasoactive intestinal peptide) and VEGFR (peptide derived from the human receptor 3 of the vascular endothelial growth factor), and the non-functionalized VLP. Previous experimental results showed that VIP-VLP was the only construct of the three that was not viable. VLPs subjected to our MD study were modeled by four complete turns of the particle involving 35 subunits of the coat protein (CP). The MD simulations showed differences in structures and interaction energies associated to the crucial contribution of the disordered N-terminal arms of CP to the global stability of the particle. These differences suggested an overall stability greater in VEGFR-VLP and smaller in VIP-VLP as compared to the unfunctionalized VLP. Our novel MD study of potyviral VLPs revealed essential clues about structure and interactions of these assembled protein particles and suggests that the computational prediction of the viability of VLPs can be a valuable contribution in the field of viral nanobiotechnology.

Structural intrinsic disorder in a functionalized potyviral coat protein as a main viability determinant of its assembled nanoparticles.

The viability of viral-derived nanoparticles (virions and VLPs) aimed to nanobiotechnological functionalizations of the coat protein (CP) of turnip mosaic virus has been studied by means of advanced computational methodologies that include molecular dynamics. The study has allowed to model the structure of the complete CP and its functionalization with three different peptides and obtain essential structural features such as order/disorder, interactions, and electrostatic potentials of their constituent domains. The results provide for the first time a dynamic view of a complete potyvirus CP, since experimental available structures so far obtained lack N- and C-terminal segments. The relevance of disorder in the most distal N-terminal subdomain, and the interaction of the less distal N-terminal subdomain with the highly ordered CP core, stand out as crucial characteristic for a viable CP. Preserving them proved of outmost importance to obtain viable potyviral CPs presenting peptides at their N-terminus.

Arabidopsis immune responses triggered by cellulose- and mixed-linked glucan-derived oligosaccharides require a group of leucine-rich repeat malectin receptor kinases.

The plant immune system perceives a diversity of carbohydrate ligands from plant and microbial cell walls through the extracellular ectodomains (ECDs) of pattern recognition receptors (PRRs), which activate pattern-triggered immunity (PTI). Among these ligands are oligosaccharides derived from mixed-linked ß-1,3/ß-1,4-glucans (MLGs; e.g. ß-1,4-D-(Glc)2 -ß-1,3-D-Glc, MLG43) and cellulose (e.g. ß-1,4-D-(Glc)3 , CEL3). The mechanisms behind carbohydrate perception in plants are poorly characterized except for fungal chitin oligosaccharides (e.g. ß-1,4-d-(GlcNAc)6 , CHI6), which involve several receptor kinase proteins (RKs) with LysM-ECDs. Here, we describe the isolation and characterization of Arabidopsis thaliana mutants impaired in glycan perception (igp) that are defective in PTI activation mediated by MLG43 and CEL3, but not by CHI6. igp1-igp4 are altered in three RKs - AT1G56145 (IGP1), AT1G56130 (IGP2/IGP3) and AT1G56140 (IGP4) - with leucine-rich-repeat (LRR) and malectin (MAL) domains in their ECDs. igp1 harbors point mutation E906K and igp2 and igp3 harbor point mutation G773E in their kinase domains, whereas igp4 is a T-DNA insertional loss-of-function mutant. Notably, isothermal titration calorimetry (ITC) assays with purified ECD-RKs of IGP1 and IGP3 showed that IGP1 binds with high affinity to CEL3 (with dissociation constant KD  = 1.19 ± 0.03 µm) and cellopentaose (KD  = 1.40 ± 0.01 µM), but not to MLG43, supporting its function as a plant PRR for cellulose-derived oligosaccharides. Our data suggest that these LRR-MAL RKs are components of a recognition mechanism for both cellulose- and MLG-derived oligosaccharide perception and downstream PTI activation in Arabidopsis.

Differences in the Elastomeric Behavior of Polyglycine-Rich Regions of Spidroin 1 and 2 Proteins.

Two different polyglycine-rich fragments were selected as representatives of major ampullate gland spidroins (MaSp) 1 and 2 types, and their behavior in a water-saturated environment was simulated within the framework of molecular dynamics (MD). The selected fragments are found in the sequences of the proteins MaSp1a and MaSp2.2a of Argiope aurantia with respective lengths of 36 amino acids (MaSp1a) and 50 amino acids (MaSp2.2s). The simulation took the fully extended ß-pleated conformation as reference, and MD was used to determine the equilibrium configuration in the absence of external forces. Subsequently, MD were employed to calculate the variation in the distance between the ends of the fragments when subjected to an increasing force. Both fragments show an elastomeric behavior that can be modeled as a freely jointed chain with links of comparable length, and a larger number of links in the spidroin 2 fragment. It is found, however, that the maximum recovery force recorded from the spidroin 2 peptide (Fmax ≈ 400 pN) is found to be significantly larger than that of the spidroin 1 (Fmax ≈ 250 pN). The increase in the recovery force of the spidroin 2 polyglycine-rich fragment may be correlated with the larger values observed in the strain at breaking of major ampullate silk fibers spun by Araneoidea species, which contain spidroin 2 proteins, compared to the material produced by spider species that lack these spidroins (RTA-clade).

Secretory protein beta-lactoglobulin in cattle stable dust may contribute to the allergy-protective farm effect.

BACKGROUND: Growing up on a cattle farm and consuming raw cow's milk protects against asthma and allergies. We expect a cattle-specific protein as active component in this farm effect. METHODS: Dust was collected from cattle and poultry stables and from mattresses of households. Urine was obtained from cattle, and ambient aerosols were sampled. Samples were analysed for BLG by SDS PAGE/immunoblot and mass spectrometry, and for association with metals by SEC-ICP-MS. PBMC of healthy donors were incubated with BLG +/- zinc, and proliferation and cytokines determined. BALB/c mice were pre-treated intranasally with stable dust extract containing BLG or depleted of BLG, and subsequent allergy response after sensitization was evaluated on antibody and symptom level. RESULTS: A major protein in dust from cattle farms and ambient air was identified as BLG. Urine from female and male cattle is a major source of BLG. In dust samples, BLG was associated with zinc. In vitro, zinc-BLG provoked significantly lower proliferation of CD4+ and CD8+ cells while inducing significantly higher levels of IFN-γ and IL-6 than the apo-BLG devoid of zinc. In vivo, pre-treatment of mice with dust extract containing BLG resulted in lower allergy symptom scores to BLG and unrelated Bet v 1 than pre-treatment with extract depleted of BLG. These in vitro and in vivo effects were independent of endotoxin. CONCLUSION: The lipocalin BLG is found in large amounts in cattle urine, accumulates in bovine dust samples and is aerosolized around farms. Its association with zinc favorably shapes the human cellular immune response towards Th1-cytokines in vitro. BLG together with zinc in stable dust protects mice from allergic sensitization. BLG with its associated ligands may in an innate manner contribute to the allergy-protective farm effect.

Micronutritional supplementation with a holoBLG-based FSMP (food for special medical purposes)-lozenge alleviates allergic symptoms in BALB/c mice: Imitating the protective farm effect.

BACKGROUND: Previously, the protective farm effect was imitated using the whey protein beta-lactoglobulin (BLG) that is spiked with iron-flavonoid complexes. Here, we formulated for clinical translation a lozenge as food for special medical purposes (FSMP) using catechin-iron complexes as ligands for BLG. The lozenge was tested in vitro and in a therapeutical BALB/c mice model. METHODS: Binding of iron-catechin into BLG was confirmed by spectroscopy and docking calculations. Serum IgE binding of children allergic or tolerating milk was assessed to loaded (holo-) versus empty (apo-) BLG and for human mast cell degranulation. BLG and Bet v 1 double-sensitized mice were orally treated with the holoBLG or placebo lozenge, and immunologically analysed after systemic allergen challenge. Human PBMCs of pollen allergic subjects were flow cytometrically assessed after stimulation with apoBLG or holoBLG using catechin-iron complexes as ligands. RESULTS: One major IgE and T cell epitope were masked by catechin-iron complexes, which impaired IgE binding of milk-allergic children and degranulation of mast cells. In mice, only supplementation with the holoBLG lozenge reduced clinical reactivity to BLG and Bet v 1, promoted Tregs, and suppressed antigen presentation. In allergic subjects, stimulation of PBMCs with holoBLG led to a significant increase of intracellular iron in circulating CD14+ cells with significantly lower expression of HLADR and CD86 compared to their stimulation with apoBLG. CONCLUSION: The FSMP lozenge targeted antigen presenting cells and dampened immune activation in human immune cells and allergic mice in an antigen-non-specific manner, thereby conferring immune resilience against allergic symptoms.

Environment and coordination of FeMo-co in the nitrogenase metallochaperone NafY.

In nitrogenase biosynthesis, the iron-molybdenum cofactor (FeMo-co) is externally assembled at scaffold proteins and delivered to the NifDK nitrogenase component by the NafY metallochaperone. Here we have used nuclear magnetic resonance, molecular dynamics, and functional analysis to elucidate the environment and coordination of FeMo-co in NafY. H121 stands as the key FeMo-co ligand. Regions near FeMo-co diverge from H121 and include the η1, α1, α2 helical lobe and a narrow path between H121 and C196.

Machine Learning of Analytical Electron Density in Large Molecules Through Message-Passing.

Machine learning milestones in computational chemistry are overshadowed by their unaccountability and the overwhelming zoo of tools for each specific task. A promising path to tackle these problems is using machine learning to reproduce physical magnitudes as a basis to derive many other properties. By using a model of the electron density consisting of an analytical expansion on a linear set of isotropic and anisotropic functions, we implemented in this work a message-passing neural network able to reproduce electron density in molecules with just a 2.5% absolute error in complex cases. We also adapted our methodology to describe electron density in large biomolecules (proteins) and to obtain atomic charges, interaction energies, and DFT energies. We show that electron density learning is a new promising avenue with a variety of forthcoming applications.

The key to the allergenicity of lipid transfer protein (LTP) ligands: A structural characterization.

Plant lipid transfer proteins are a large family that can be found in all land plants. They have a hydrophobic cavity that allows them to harbor lipids and facilitates their traffic between membranes. However, in humans, this plant protein family is responsible for the main food allergies in the Mediterranean area. Nevertheless, not only the protein itself but also its ligand is relevant for allergic sensitization. The main aim of the present work is to analyse the natural ligands carried by four allergenic LTPs (Tri a 14, Art v 3, Par j 2, and Ole e 7), compared with the previously identified ligand of Pru p 3 (CPT-PHS ligand), and clarify their role within the immunological reactions. Results showed that the ligands of the LTPs studied shared a chemical identity, in which the presence of a polar head was essential to the protein-ligand binding. This ligand was transported through a skin cellular model, and phosphorylated phytosphingosine could be detected as result of cell metabolism. Since sphingosine kinase 1 was overexpressed in keratinocytes incubated with the LTP-ligand complex, this enzyme might be responsible for the phosphorylation of the phytosphingosine fraction of the CPT-PHS ligand. This way, phytosphingosine-1-phosphate could be mimicking the role of the human inflammatory mediator sphingosine-1-phosphate, explaining why LTPs are associated with more severe allergic responses. In conclusion, this work contributes to the understanding of the chemical nature and behavior of lipid ligands carried by allergens, which would help to gain insight into their role during allergic sensitization.

Cow's milk protein ß-lactoglobulin confers resilience against allergy by targeting complexed iron into immune cells.

BACKGROUND: Beta-lactoglobulin (BLG) is a bovine lipocalin in milk with an innate defense function. The circumstances under which BLG is associated with tolerance of or allergy to milk are not understood. OBJECTIVE: Our aims were to assess the capacity of ligand-free apoBLG versus loaded BLG (holoBLG) to protect mice against allergy by using an iron-quercetin complex as an exemplary ligand and to study the molecular mechanisms of this protection. METHODS: Binding of iron-quercetin to BLG was modeled and confirmed by spectroscopy and docking calculations. Serum IgE binding to apoBLG and holoBLG in children allergic to milk and children tolerant of milk was assessed. Mice were intranasally treated with apoBLG versus holoBLG and analyzed immunologically after systemic challenge. Aryl hydrocarbon receptor (AhR) activation was evaluated with reporter cells and Cyp1A1 expression. Treated human PBMCs and human mast cells were assessed by fluorescence-activated cell sorting and degranulation, respectively. RESULTS: Modeling predicted masking of major IgE and T-cell epitopes of BLG by ligand binding. In line with this modeling, IgE binding in children allergic to milk was reduced toward holoBLG, which also impaired degranulation of mast cells. In mice, only treatments with holoBLG prevented allergic sensitization and anaphylaxis, while sustaining regulatory T cells. BLG facilitated quercetin-dependent AhR activation and, downstream of AhR, lung Cyp1A1 expression. HoloBLG shuttled iron into monocytic cells and impaired their antigen presentation. CONCLUSION: The cargo of holoBLG is decisive in preventing allergy in vivo. BLG without cargo acted as an allergen in vivo and further primed human mast cells for degranulation in an antigen-independent fashion. Our data provide a mechanistic explanation why the same proteins can act either as tolerogens or as allergens.

Analytical Model of Electron Density and Its Machine Learning Inference.

We present an analytical model representation of the electron density ρ(r) in molecules in the form of expansions of a few functions (exponentials and a Gaussian) per atom. Based on a former analytical model of ρ(r) in atoms, we devised its molecular implementation by introducing the anisotropy inherent in the electron distribution of atoms in molecules by means of proper anisotropic functions. The resulting model named A2MD (anisotropic analytical model of density) takes an analytical form highly suitable for obtaining the electron density in large biomolecules as its computational cost scales linearly with the number of atoms. To obtain the parameters of the model, we first devised a fitting procedure to reference electron densities obtained in ab initio correlated quantum calculations. Second, in order to skip costly ab initio calculations, we also developed a machine learning (ML)-based predictor that used neural networks trained on broad molecular datasets to determine the parameters of the model. The resulting ML methodology that we named A2MDnet (A2MD network-trained) was able to provide reliable electron densities as a basis to predict molecular features without requiring quantum calculations. The results presented together with the low computational scaling associated to the A2MD representation of ρ(r) suggest potential applications to obtain reliable electron densities and ρ(r)-based molecular properties in biomacromolecules.

A computational structural study on the DNA-protecting role of the tardigrade-unique Dsup protein.

The remarkable ability of tardigrades to withstand a wide range of physical and chemical extremes has attracted a considerable interest in these small invertebrates, with a particular focus on the protective roles of proteins expressed during such conditions. The discovery that a tardigrade-unique protein named Dsup (damage suppressor) protects DNA from damage produced by radiation and radicals, has raised expectations concerning its potential applications in biotechnology and medicine. We present in this paper what might be dubbed a "computational experiment" on the Dsup-DNA system. By means of molecular modelling, calculations of electrostatic potentials and electric fields, and all-atom molecular dynamics simulations, we obtained a dynamic picture of the Dsup-DNA interaction. Our results suggest that the protein is intrinsically disordered, which enables Dsup to adjust its structure to fit DNA shape. Strong electrostatic attractions and high protein flexibility drive the formation of a molecular aggregate in which Dsup shields DNA. While the precise mechanism of DNA protection conferred by Dsup remains to be elucidated, our study provides some molecular clues of their association that could be of interest for further investigation in this line.

Association between flower stalk elongation, an Arabidopsis developmental trait, and the subcellular location and movement dynamics of the nonstructural protein P3 of Turnip mosaic virus.

Virus infections affect plant developmental traits but this aspect of the interaction has not been extensively studied so far. Two strains of Turnip mosaic virus differentially affect Arabidopsis development, especially flower stalk elongation, which allowed phenotypical, cellular, and molecular characterization of the viral determinant, the P3 protein. Transiently expressed wild-type green fluorescent protein-tagged P3 proteins of both strains and selected mutants of them revealed important differences in their behaviour as endoplasmic reticulum (ER)-associated peripheral proteins flowing along the reticulum, forming punctate accumulations. Three-dimensional (3D) model structures of all expressed P3 proteins were computationally constructed through I-TASSER protein structure predictions, which were used to compute protein surfaces and map electrostatic potentials to characterize the effect of amino acid changes on features related to protein interactions and to phenotypical and subcellular results. The amino acid at position 279 was the main determinant affecting stalk development. It also determined the speed of ER-flow of the expressed proteins and their final location. A marked change in the protein surface electrostatic potential correlated with changes in subcellular location. One single amino acid in the P3 viral protein determines all the analysed differential characteristics between strains differentially affecting flower stalk development. A model proposing a role of the protein in the intracellular movement of the viral replication complex, in association with the viral 6K2 protein, is proposed. The type of association between both viral proteins could differ between the strains.

Structural Dynamics of the Lipid Antigen-Binding Site of CD1d Protein.

CD1 molecules present lipid antigens to T-cells in early stages of immune responses. Whereas CD1‒lipid‒T-cell receptors interactions are reasonably understood, molecular details on initial trafficking and loading of lipids onto CD1 proteins are less complete. We present a molecular dynamics (MD) study of human CD1d, the isotype that activates iNKT cells. MD simulations and calculations of properties and Poisson-Boltzmann electrostatic potentials were used to explore the dynamics of the antigen-binding domain of the apo-form, CD1d complexes with three lipid-antigens that activate iNKT cells and CD1d complex with GM2AP, a protein that assists lipid loading onto CD1 molecules in endosomes/lysosomes. The study was done at pH 7 and 4.5, values representative of strongly acidic environments in endosomal compartments. Our findings revealed dynamic features of the entrance to the hydrophobic channels of CD1d modulated by two α helices with sensitivity to the type of lipid. We also found lipid- and pH-dependent dynamic changes in three exposed tryptophans unique to CD1d among the five human CD1 isotypes. On the basis of modelled structures, our data also revealed external effects produced by the helper protein GM2AP only when it interacts in its open form, thus suggesting that the own assistant protein also adapts conformation to association with CD1d.

Dynamic plasticity of the lipid antigen-binding site of CD1d is crucially favoured by acidic pH and helper proteins.

CD1 molecules present lipid antigens for recognition by T-cell receptors (TCRs). Although a reasonably detailed picture of the CD1-lipid-TCR interaction exists, the initial steps regarding lipid loading onto and exchange between CD1 proteins remain elusive. The hydrophobic nature of lipids and the fact that CD1 molecules are unable to extract lipids from membranes raise the need for the assistance of helper proteins in lipid trafficking. However, the experimental study of this traffic in the endosomal compartments at which it occurs is so challenging that computational studies can help provide mechanistic insight into the associated processes. Here we present a multifaceted computational approach to obtain dynamic structural data on the human CD1d isotype. Conformational dynamics analysis shows an intrinsic flexibility associated with the protein architecture. Electrostatic properties together with molecular dynamics results for CD1d complexes with several lipids and helper proteins unravel the high dynamic plasticity of the antigen-binding site that is crucially favoured by acidic pH and the presence of helper proteins.

Interaction of Alt a 1 with SLC22A17 in the airway mucosa.

BACKGROUND: Despite all the efforts made up to now, the reasons that facilitate a protein becoming an allergen have not been elucidated yet. Alt a 1 protein is the major fungal allergen responsible for chronic asthma, but little is known about its immunological activity. Our main purpose was to investigate the ligand-dependent interactions of Alt a 1 in the human airway epithelium. METHODS: Alt a 1 with and without its ligand (holo- and apo- forms) was incubated with the pulmonary epithelial monolayer model, Calu-3 cells. Allergen transport and cytokine production were measured. Pull-down and immunofluorescence assays were employed to identify the receptor of Alt a 1 using the epithelial cell model and mouse tissues. Receptor-allergen-ligand interactions were analyzed by computational modeling. RESULTS: The holo-form could activate human monocytes, PBMCs, and polarized airway epithelial (Calu-3) cell lines. The allergen was also transported through the monolayer, without any alteration of the epithelial integrity (TEER). Alt a 1 also induced the production of proinflammatory IL8 and specific epithelial cytokines (IL33 and IL25) by Calu-3 cells. The interaction between epithelial cells and holo-Alt a 1 was found to be mediated by the SLC22A17 receptor, and its recognition of Alt a 1 was explained in structural terms. CONCLUSIONS: Our findings identified the Alt a 1 ligand as a central player in the interaction of the allergen with airway mucosa, shedding light into its potential role in the immunological response, while unveiling its potential as a new target for therapy intervention.

Energy Landscapes of Ligand Motion Inside the Tunnel-Like Cavity of Lipid Transfer Proteins: The Case of the Pru p 3 Allergen.

Allergies are a widespread problem in western countries, affecting a large part of the population, with levels of prevalence increasingly rising due to reasons still not understood. Evidence accumulated in recent years points to an essential role played by ligands of allergen proteins in the sensitization phase of allergies. In this regard, we recently identified the natural ligand of Pru p 3, a lipid transfer protein, a major allergen from peach fruit and a model of food allergy. The ligand of Pru p 3 has been shown to play a key role in the sensitization to peach and to other plant food sources that provoke cross-reactivity in a large proportion of patients allergic to peach. However, the question of which is the binding pose of this ligand in its carrier protein, and how it can be transferred to receptors of the immune system where it develops its function as a coadjuvant was not elucidated. In this work, different molecular dynamics simulations have been considered as starting points to study the properties of the ligand⁻protein system in solution. Besides, an energy landscape based on collective variables that describe the process of ligand motion within the cavity of Pru p 3 was obtained by using well-tempered metadynamics. The simulations revealed the differences between distinct binding modes, and also revealed important aspects of the motion of the ligand throughout its carrier protein, relevant to its binding⁻unbinding process. Our findings are potentially interesting for studying protein⁻ligand systems beyond the specific case of the allergen protein dealt with here.