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.

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.

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.

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.

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.

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.

Non-branched ß-1,3-glucan oligosaccharides trigger immune responses in Arabidopsis.

Fungal cell walls, which are essential for environmental adaptation and host colonization by the fungus, have been evolutionarily selected by plants and animals as a source of microbe-associated molecular patterns (MAMPs) that, upon recognition by host pattern recognition receptors (PRRs), trigger immune responses conferring disease resistance. Chito-oligosaccharides [ß-1,4-N-acetylglucosamine oligomers, (GlcNAc)n ] are the only glycosidic structures from fungal walls that have been well-demonstrated to function as MAMPs in plants. Perception of (GlcNAc)4-8 by Arabidopsis involves CERK1, LYK4 and LYK5, three of the eight members of the LysM PRR family. We found that a glucan-enriched wall fraction from the pathogenic fungus Plectosphaerella cucumerina which was devoid of GlcNAc activated immune responses in Arabidopsis wild-type plants but not in the cerk1 mutant. Using this differential response, we identified the non-branched 1,3-ß-d-(Glc) hexasaccharide as a major fungal MAMP. Recognition of 1,3-ß-d-(Glc)6 was impaired in cerk1 but not in mutants defective in either each of the LysM PRR family members or in the PRR-co-receptor BAK1. Transcriptomic analyses of Arabidopsis plants treated with 1,3-ß-d-(Glc)6 further demonstrated that this fungal MAMP triggers the expression of immunity-associated genes. In silico docking analyses with molecular mechanics and solvation energy calculations corroborated that CERK1 can bind 1,3-ß-d-(Glc)6 at effective concentrations similar to those of (GlcNAc)4 . These data support that plants, like animals, have selected as MAMPs the linear 1,3-ß-d-glucans present in the walls of fungi and oomycetes. Our data also suggest that CERK1 functions as an immune co-receptor for linear 1,3-ß-d-glucans in a similar way to its proposed function in the recognition of fungal chito-oligosaccharides and bacterial peptidoglycan MAMPs.

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.

Pleiotropic Effects of Resistance-Breaking Mutations on Particle Stability Provide Insight into Life History Evolution of a Plant RNA Virus.

In gene-for-gene host-virus interactions, virus evolution to infect and multiply in previously resistant host genotypes, i.e., resistance breaking, is a case of host range expansion, which is predicted to be associated with fitness penalties. Negative effects of resistance-breaking mutations on within-host virus multiplication have been documented for several plant viruses. However, understanding virus evolution requires analyses of potential trade-offs between different fitness components. Here we analyzed whether coat protein (CP) mutations in Pepper mild mottle virus that break L-gene resistance in pepper affect particle stability and, thus, survival in the environment. For this purpose, CP mutations determining the overcoming of L 3 and L 4 resistance alleles were introduced in biologically active cDNA clones. The kinetics of the in vitro disassembly of parental and mutant particles were compared under different conditions. Resistance-breaking mutations variously affected particle stability. Structural analyses identified the number and type of axial and side interactions of adjacent CP subunits in virions, which explained differences in particle stability and contribute to understanding of tobamovirus disassembly. Resistance-breaking mutations also affected virus multiplication and virulence in the susceptible host, as well as infectivity. The sense and magnitude of the effects of resistance-breaking mutations on particle stability, multiplication, virulence, or infectivity depended on the specific mutation rather than on the ability to overcome the different resistance alleles, and effects on different traits were not correlated. Thus, the results do not provide evidence of links or trade-offs between particle stability, i.e., survival, and other components of virus fitness or virulence.IMPORTANCE The effect of survival on virus evolution remains underexplored, despite the fact that life history trade-offs may constrain virus evolution. We approached this topic by analyzing whether breaking of L-gene resistance in pepper by Pepper mild mottle virus, determined by coat protein (CP) mutations, is associated with reduced particle stability and survival. Resistance-breaking mutations affected particle stability by altering the interactions between CP subunits. However, the sense and magnitude of these effects were unrelated to the capacity to overcome different resistance alleles. Thus, resistance breaking was not traded with survival. Resistance-breaking mutations also affected virus fitness within the infected host, virulence, and infectivity in a mutation-specific manner. Comparison of the effects of CP mutations on these various traits indicates that there are neither trade-offs nor positive links between survival and other life history traits. These results demonstrate that trade-offs between life history traits may not be a general constraint in virus evolution.

Expression and Interaction Analysis among Saffron ALDHs and Crocetin Dialdehyde.

In saffron, the cleavage of zeaxanthin by means of CCD2 generates crocetin dialdehyde, which is then converted by an unknown aldehyde dehydrogenase to crocetin. A proteome from saffron stigma was released recently and, based on the expression pattern and correlation analyses, five aldehyde dehydrogenases (ALDHs) were suggested as possible candidates to generate crocetin from crocetin dialdehydes. We selected four of the suggested ALDHs and analyzed their expression in different tissues, determined their activity over crocetin dialdehyde, and performed structure modeling and docking calculation to find their specificity. All the ALDHs were able to convert crocetin dialdehyde to crocetin, but two of them were stigma tissue-specific. Structure modeling and docking analyses revealed that, in all cases, there was a high coverage of residues in the models. All of them showed a very close conformation, indicated by the low root-mean-square deviation (RMSD) values of backbone atoms, which indicate a high similarity among them. However, low affinity between the enzymes and the crocetin dialdehyde were observed. Phylogenetic analysis and binding affinities calculations, including some ALDHs from Gardenia jasmonoides, Crocus sieberi, and Buddleja species that accumulate crocetin and Bixa orellana synthetizing the apocarotenoid bixin selected on their expression pattern matching with the accumulation of either crocins or bixin, pointed out that family 2 C4 members might be involved in the conversion of crocetin dialdehyde to crocetin with high specificity.

A Comparative Study of Human Saposins.

Saposins are small proteins implicated in trafficking and loading of lipids onto Cluster of Differentiation 1 (CD1) receptor proteins that in turn present lipid antigens to T cells and a variety of T-cell receptors, thus playing a crucial role in innate and adaptive immune responses in humans. Despite their low sequence identity, the four types of human saposins share a similar folding pattern consisting of four helices linked by three conserved disulfide bridges. However, their lipid-binding abilities as well as their activities in extracting, transporting and loading onto CD1 molecules a variety of sphingo- and phospholipids in biological membranes display two striking characteristics: a strong pH-dependence and a structural change between a compact, closed conformation and an open conformation. In this work, we present a comparative computational study of structural, electrostatic, and dynamic features of human saposins based upon their available experimental structures. By means of structural alignments, surface analyses, calculation of pH-dependent protonation states, Poisson-Boltzmann electrostatic potentials, and molecular dynamics simulations at three pH values representative of biological media where saposins fulfill their function, our results shed light into their intrinsic features. The similarities and differences in this class of proteins depend on tiny variations of local structural details that allow saposins to be key players in triggering responses in the human immune system.

Identification of the ligand of Pru p 3, a peach LTP.

KEY MESSAGE: Pru p 3, a peach LTP, is located in pollinated flower styles and secreting downy hairs, transporting a derivative of camptothecin bound to phytosphingosine. Pru p 3 may inhibit a second pollination and may keep away herbivores until seed maturation. The allergen Pru p 3, a peach lipid transfer protein, has been well studied. However, its physiological function remains to be elucidated. Our results showed that Pru p 3 usually carries a lipid ligand that play an essential role in its function in plants. Using ESI-qToF, we observed that the ligand was a derivative of camptothecin binding to phytosphingosine, wich that is inserted into the hydrophobic tunnel of the protein. In addition, the described ligand displayed topoisomerase I activity inhibition and self-fluorescence, both recognized as camptothecin properties. During flower development, the highest expression of Pru p 3 was detected in the styles of pollinated flowers, in contrast to its non-expression in unpollinated pistils, where expression decreased after anthesis. During ripening, the expression of Pru p 3 were observed mainly in peel but not in pulp. In this sense, Pru p 3 protein was also localized in trichomes covering the fruit epidermis.

Computational study of pH-dependent oligomerization and ligand binding in Alt a 1, a highly allergenic protein with a unique fold.

Alt a 1 is a highly allergenic protein from Alternaria fungi responsible for several respiratory diseases. Its crystal structure revealed a unique ß-barrel fold that defines a new family exclusive to fungi and forms a symmetrical dimer in a butterfly-like shape as well as tetramers. Its biological function is as yet unknown but its localization in cell wall of Alternaria spores and its interactions in the onset of allergy reactions point to a function to transport ligands. However, at odds with binding features in ß-barrel proteins, monomeric Alt a 1 seems unable to harbor ligands because the barrel is too narrow. Tetrameric Alt a 1 is able to bind the flavonoid quercetin, yet the stability of the aggregate and the own ligand binding are pH-dependent. At pH 6.5, which Alt a 1 would meet when secreted by spores in bronchial epithelium, tetramer-quercetin complex is stable. At pH 5.5, which Alt a 1 would meet in apoplast when infecting plants, the complex breaks down. By means of a combined computational study that includes docking calculations, empirical pKa estimates, Poisson-Boltzmann electrostatic potentials, and Molecular Dynamics simulations, we identified a putative binding site at the dimeric interface between subunits in tetramer. We propose an explanation on the pH-dependence of both oligomerization states and protein-ligand affinity of Alt a 1 in terms of electrostatic variations associated to distinct protonation states at different pHs. The uniqueness of this singular protein can thus be tracked in the combination of all these features.

A theoretical study on the reaction of ozone with aqueous iodide.

Atmospheric iodine chemistry plays a key role in tropospheric ozone catalytic destruction, new particle formation, and as one of the possible sinks of gaseous polar elemental mercury. Moreover, it has been recently proposed that reaction of ozone with iodide on the sea surface could be the major contributor to the chemical loss of atmospheric ozone. However, the mechanism of the reaction between aqueous iodide and ozone is not well known. The aim of this paper is to improve the understanding of such a mechanism. In this paper, an ab initio study of the reaction of aqueous iodide and ozone is presented, evaluating thermodynamic data of the different reactions proposed in previous experimental studies. In addition, the structures, energetics and possible evolution of the key IOOO(-) intermediate are discussed for the first time.

Maturation of Rhizobium leguminosarum hydrogenase in the presence of oxygen requires the interaction of the chaperone HypC and the scaffolding protein HupK.

[NiFe] hydrogenases are key enzymes for the energy and redox metabolisms of different microorganisms. Synthesis of these metalloenzymes involves a complex series of biochemical reactions catalyzed by a plethora of accessory proteins, many of them required to synthesize and insert the unique NiFe(CN)2CO cofactor. HypC is an accessory protein conserved in all [NiFe] hydrogenase systems and involved in the synthesis and transfer of the Fe(CN)2CO cofactor precursor. Hydrogenase accessory proteins from bacteria-synthesizing hydrogenase in the presence of oxygen include HupK, a scaffolding protein with a moderate sequence similarity to the hydrogenase large subunit and proposed to participate as an intermediate chaperone in the synthesis of the NiFe cofactor. The endosymbiotic bacterium Rhizobium leguminosarum contains a single hydrogenase system that can be expressed under two different physiological conditions: free-living microaerobic cells (∼ 12 µm O2) and bacteroids from legume nodules (∼ 10-100 nm O2). We have used bioinformatic tools to model HupK structure and interaction of this protein with HypC. Site-directed mutagenesis at positions predicted as critical by the structural analysis have allowed the identification of HupK and HypC residues relevant for the maturation of hydrogenase. Mutant proteins altered in some of these residues show a different phenotype depending on the physiological condition tested. Modeling of HypC also predicts the existence of a stable HypC dimer whose presence was also demonstrated by immunoblot analysis. This study widens our understanding on the mechanisms for metalloenzyme biosynthesis in the presence of oxygen.

Bet v 1 from birch pollen is a lipocalin-like protein acting as allergen only when devoid of iron by promoting Th2 lymphocytes.

It is hypothesized that allergens are at the borderline of self and non-self and, through as yet elusive circumstances, mount a Th2 response for allergic sensitization. The major birch pollen allergen Bet v 1 is considered the prototype for the PR-10 protein family causing respiratory allergy. Here, we give structural evidence that Bet v 1 is a lipocalin-like protein with a striking structural resemblance to human lipocalin 2. Lipocalin 2 is highly expressed in the lung where it exerts immunoregulatory functions dependent on being loaded with siderophore-bound iron (holo-form) or not (apo-form). We demonstrate that similar to lipocalin 2, Bet v 1 is capable of binding iron via catechol-based siderophores. Thereby, calculated Kd values of 66 nm surpassed affinities to known ligands nearly by a power of 10. Moreover, we give functional evidence of the immunomodulatory capacity of Bet v 1 being dependent on its iron-loaded state. When incubated to human immune cells, only the apo-form of Bet v 1, but not the holo-form, was able to promote Th2 cells secreting IL13. These results provide for the first time a functional understanding on the allergenicity of Bet v 1 and a basis for future allergen immunotherapies counteracting Th2 immune responses on a molecular basis.

Host resistance selects for traits unrelated to resistance-breaking that affect fitness in a plant virus.

The acquisition by parasites of the capacity to infect resistant host genotypes, that is, resistance-breaking, is predicted to be hindered by across-host fitness trade-offs. All analyses of costs of resistance-breaking in plant viruses have focused on within-host multiplication without considering other fitness components, which may limit understanding of virus evolution. We have reported that host range expansion of tobamoviruses on L-gene resistant pepper genotypes was associated with severe within-host multiplication penalties. Here, we analyze whether resistance-breaking costs might affect virus survival in the environment by comparing tobamovirus pathotypes differing in infectivity on L-gene resistance alleles. We predicted particle stability from structural models, analyzed particle stability in vitro, and quantified virus accumulation in different plant organs and virus survival in the soil. Survival in the soil differed among tobamovirus pathotypes and depended on differential stability of virus particles. Structure model analyses showed that amino acid changes in the virus coat protein (CP) responsible for resistance-breaking affected the strength of the axial interactions among CP subunits in the rod-shaped particle, thus determining its stability and survival. Pathotypes ranked differently for particle stability/survival and for within-host accumulation. Resistance-breaking costs in survival add to, or subtract from, costs in multiplication according to pathotype. Hence, differential pathotype survival should be considered along with differential multiplication to understand the evolution of the virus populations. Results also show that plant resistance, in addition to selecting for resistance-breaking and for decreased multiplication, also selects for changes in survival, a trait unrelated to the host-pathogen interaction that may condition host range expansion.