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.

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.

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) is an ancient protein critical for CO2-fixation and global biogeochemistry. Form-I RuBisCO complexes uniquely harbor small subunits that form a hexadecameric complex together with their large subunits. The small subunit protein is thought to have significantly contributed to RuBisCO's response to the atmospheric rise of O2 ∼2.5 billion years ago, marking a pivotal point in the enzyme's evolutionary history. Here, we performed a comprehensive evolutionary analysis of extant and ancestral RuBisCO sequences and structures to explore the impact of the small subunit's earliest integration on the molecular dynamics of the overall complex. Our simulations suggest that the small subunit restricted the conformational flexibility of the large subunit early in its history, impacting the evolutionary trajectory of the Form-I RuBisCO complex. Molecular dynamics investigations of CO2 and O2 gas distribution around predicted ancient RuBisCO complexes suggest that a proposed “CO2-reservoir” role for the small subunit is not conserved throughout the enzyme's evolutionary history. The evolutionary and biophysical response of RuBisCO to changing atmospheric conditions on ancient Earth showcase multi-level and trackable responses of enzymes to environmental shifts over long timescales.

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.

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.

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.