Search

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). [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Abstract Biological production and oxidation of hydrogen is mediated by hydrogenases, key enzymes for these energy-relevant reactions. Synthesis of [NiFe] hydrogenases involves a complex series of biochemical reactions to assemble protein subunits and metallic cofactors required for enzyme function. A final step in this biosynthetic pathway is the processing of a C-terminal tail (CTT) from its large subunit, thus allowing proper insertion of nickel in the unique NiFe(CN) 2 CO cofactor present in these enzymes. In silico modelling and Molecular Dynamics (MD) analyses of processed vs. unprocessed forms of Rhizobium leguminosarum bv. viciae (Rlv) hydrogenase large subunit HupL showed that its CTT (residues 582-596) is an intrinsically disordered region (IDR) that likely provides the required flexibility to the protein for the final steps of proteolytic maturation. Prediction of pKa values of ionizable side chains in both forms of the enzyme's large subunit also revealed that the presence of the CTT strongly modify the protonation state of some key residues around the active site. Furthermore, MD simulations and mutant analyses revealed that two glutamate residues (E27 in the N-terminal region and E589 inside the CTT) likely contribute to the process of nickel incorporation into the enzyme. Computational analysis also revealed structural details on the interaction of Rlv hydrogenase LSU with the endoprotease HupD responsible for the removal of CTT. Highlights • C-terminal tails of NiFe hydrogenases are Intrinsically Disordered Region (IDRs). • IDRs might provide flexibility for interaction with protease and other accessory proteins. • Modelling of [NiFe] hydrogenase unprocessed form reveals changes in the environment of key residues following processing. • Molecular dynamics and mutant analyses indicate that residues E27 and E589 are relevant for nickel incorporation. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

In the search for factors that make a protein allergenic (an issue that remains so far elusive) some common features of allergens such as small size, high stability and lipid binding are recognized in spite of their structural diversity. Other relevant but still poorly understood feature is their capability to form homodimers. We investigated by means of Molecular Dynamics (MD) calculations the stability in solution of several dimers of three major allergens from Alternaria mold, birch pollen, and peach fruit known to play essential roles in allergic disease. By running 20 ns MD simulations we found essential properties on solution that provide information of interest on their dimerization, stability of their epitopes and dynamical features of ligand binding cavities. Our results show that three essential allergen proteins display a distinct behavior on their trends to form homodimers in solution. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]