Oleic Acid Could Act as a Channel Blocker in the Inhibition of nAChR: Insights from Molecular Dynamics Simulations.

The activation of the muscular nicotinic acetylcholine receptor (nAChR) produces the opening of the channel, with the consequent increase in the permeability of cations, triggering an excitatory signal. Free fatty acids (FFA) are known to modulate the activity of the receptor as noncompetitive antagonists, acting at the membrane-AChR interface. We present molecular dynamics simulations of a model of nAChR in a desensitized closed state embedded in a lipid bilayer in which distinct membrane phospholipids were replaced by two different monounsaturated FFA that differ in the position of a double bond. This allowed us to detect and describe that the cis-18:1ω-9 FFA were located at the interface between the transmembrane segments of α2 and γ subunits diffused into the channel lumen with the consequent potential ability to block the channel to the passage of ions.

IDP Force Fields Applied to Model PPII-Rich 33-mer Gliadin Peptides.

The 33-mer gliadin peptide and its deamidated metabolite, 33-mer DGP, are the immunodominant peptides responsible for the adaptive immune response in celiac disease (CD). CD is a complex autoimmune chronic disorder triggered by gluten ingestion that affects the small intestine and affects ∼1% of the global population. The 33-mers are polyproline II-rich (PPII) and intrinsically disordered peptides (IDPs), whose structures remain elusive. We sampled the conformational ensembles of both 33-mer peptides via molecular dynamics simulations employing two force fields (FFs) (Amber ff03ws and Amber ff99SB-disp) specifically validated for other IDPs. Our results show that both FFs allow the extensive exploration of the conformational landscape, which was not possible with the standard FF GROMOS53A6 reported before. Clustering analysis of the trajectories showed that the five largest clusters (78-88% of the total structures) present elongated, semielongated, and curved conformations in both FFs. Large average radius of gyration and solvent-exposed surfaces characterized these structures. While the structures sampled are similar, the Amber ff99SB-disp trajectories explored folded conformations with a higher probability. In addition, PPII secondary structure was preserved throughout the trajectories (58-73%) together with a non-negligible content of ß structures (11-23%), in agreement with previous experimental results. This work represents the initial step in studying further the interaction of these peptides with other biologically relevant molecules, which could lead to finally disclose the molecular events that lead to CD.

Molecular mechanisms of 33-mer gliadin peptide oligomerisation.

The proteolytic resistant 33-mer gliadin peptide is an immunodominant fragment in gluten and responsible for the celiac disease and other gluten-related disorders. Meanwhile, the primary structure of the 33-mer is associated with the adaptive immune response in celiac patients, and the structural transformation of the 33-mer into protofilaments activates a primordial innate immune response in human macrophages. This means that accumulation, oligomerisation and structural transformation of the 33-mer could be the unknown first event that triggers the disease. Herein, we reveal the early stepwise mechanism of 33-mer oligomerisation by combining multiple computational simulations, tyrosine cross-linking, fluorescence spectroscopy and circular dichroism experiments. Our theoretical findings demonstrated that the partial charge distribution along the 33-mer molecule and the presence of glutamine that favours H-bonds between the oligomers are the driving forces that trigger oligomerisation. The high content of proline is critical for the formation of the flexible PPII secondary structure that led to a ß structure transition upon oligomerisation. Experimentally, we stabilised the 33-mer small oligomers by dityrosine cross-linking, detecting from dimers to higher molecular weight oligomers, which confirmed our simulations. The relevance of 33-mer oligomers as a trigger of the disease as well as its inhibition may be a novel therapeutic strategy for the treatment of gluten-related disorders.

Multiscale modelization in a small virus: Mechanism of proton channeling and its role in triggering capsid disassembly.

In this work, we assess a previously advanced hypothesis that predicts the existence of ion channels in the capsid of small and non-enveloped icosahedral viruses. With this purpose we examine Triatoma Virus (TrV) as a case study. This virus has a stable capsid under highly acidic conditions but disassembles and releases the genome in alkaline environments. Our calculations range from a subtle sub-atomic proton interchange to the dismantling of a large-scale system representing several million of atoms. Our results provide structure-based explanations for the three roles played by the capsid to enable genome release. First, we observe, for the first time, the formation of a hydrophobic gate in the cavity along the five-fold axis of the wild-type virus capsid, which can be disrupted by an ion located in the pore. Second, the channel enables protons to permeate the capsid through a unidirectional Grotthuss-like mechanism, which is the most likely process through which the capsid senses pH. Finally, assuming that the proton leak promotes a charge imbalance in the interior of the capsid, we model an internal pressure that forces shell cracking using coarse-grained simulations. Although qualitatively, this last step could represent the mechanism of capsid opening that allows RNA release. All of our calculations are in agreement with current experimental data obtained using TrV and describe a cascade of events that could explain the destabilization and disassembly of similar icosahedral viruses.

Conserved charged amino acids are key determinants for fatty acid binding proteins (FABPs)-membrane interactions. A multi-methodological computational approach.

Based on the analysis of the mechanism of ligand transfer to membranes employing in vitro methods, Fatty Acid Binding Protein (FABP) family has been divided in two subgroups: collisional and diffusional FABPs. Although the collisional mechanism has been well characterized employing in vitro methods, the structural features responsible for the difference between collisional and diffusional mechanisms remain uncertain. In this work, we have identified the amino acids putatively responsible for the interaction with membranes of both, collisional and diffusional, subgroups of FABPs. Moreover, we show how specific changes in FABPs' structure could change the mechanism of interaction with membranes. We have computed protein-membrane interaction energies for members of each subgroup of the family, and performed Molecular Dynamics simulations that have shown different configurations for the initial interaction between FABPs and membranes. In order to generalize our hypothesis, we extended the electrostatic and bioinformatics analysis over FABPs of different mammalian genus. Also, our methodological approach could be used for other systems involving protein-membrane interactions.

Echinococcus granulosus antigen B: a Hydrophobic Ligand Binding Protein at the host-parasite interface.

Lipids are mainly solubilized by various families of lipid binding proteins which participate in their transport between tissues as well as cell compartments. Among these families, Hydrophobic Ligand Binding Proteins (HLBPs) deserve special consideration since they comprise intracellular and extracellular members, are able to bind a variety of fatty acids, retinoids and some sterols, and are present exclusively in cestodes. Since these parasites have lost catabolic and biosynthetic pathways for fatty acids and cholesterol, HLBPs are likely relevant for lipid uptake and transportation between parasite and host cells. Echinococcus granulosus antigen B (EgAgB) is a lipoprotein belonging to the HLBP family, which is very abundant in the larval stage of this parasite. Herein, we review the literature on EgAgB composition, structural organization and biological properties, and propose an integrated scenario in which this parasite HLBP contributes to adaptation to mammalian hosts by meeting both metabolic and immunomodulatory parasite demands.

Structure of the Triatoma virus capsid.

The members of the Dicistroviridae family are non-enveloped positive-sense single-stranded RNA (+ssRNA) viruses pathogenic to beneficial arthropods as well as insect pests of medical importance. Triatoma virus (TrV), a member of this family, infects several species of triatomine insects (popularly named kissing bugs), which are vectors for human trypanosomiasis, more commonly known as Chagas disease. The potential use of dicistroviruses as biological control agents has drawn considerable attention in the past decade, and several viruses of this family have been identified, with their targets covering honey bees, aphids and field crickets, among others. Here, the crystal structure of the TrV capsid at 2.5 Å resolution is reported, showing that as expected it is very similar to that of Cricket paralysis virus (CrPV). Nevertheless, a number of distinguishing structural features support the introduction of a new genus (Triatovirus; type species TrV) under the Dicistroviridae family. The most striking differences are the absence of icosahedrally ordered VP4 within the infectious particle and the presence of prominent projections that surround the fivefold axis. Furthermore, the structure identifies a second putative autoproteolytic DDF motif in protein VP3, in addition to the conserved one in VP1 which is believed to be responsible for VP0 cleavage during capsid maturation. The potential meaning of these new findings is discussed.

Similar structures but different mechanisms: Prediction of FABPs-membrane interaction by electrostatic calculation.

The role of fatty acid binding proteins as intracellular fatty acid transporters may require their direct interaction with membranes. In this way different mechanisms have been previously characterized through experimental studies suggesting different models for FABPs-membrane association, although the process in which the molecule adsorbs to the membrane remains to be elucidated. To estimate the importance of the electrostatic energy in the FABP-membrane interaction, we computationally modeled the interaction of different FABPs with both anionic and neutral membranes. Free Electrostatic Energy of Binding (dE), was computed using Finite Difference Poisson Boltzmann Equation (FDPB) method as implemented in APBS (Adaptive Poisson Boltzmann Solver). Based on the computational analysis, it is found that recruitment to membranes is facilitated by non-specific electrostatic interactions. Also energetic analysis can quantitatively differentiate among the mechanisms of membrane association proposed and determinate the most energetically favorable configuration for the membrane-associated states of different FABPs. This type of calculations could provide a starting point for further computational or experimental analysis.

Prediction of the most favorable configuration in the ACBP-membrane interaction based on electrostatic calculations.

Acyl-CoA binding proteins (ACBPs) are highly conserved 10 kDa cytosolic proteins that bind medium- and long-chain acyl-CoA esters. They act as intracellular carriers of acyl-CoA and play a role in acyl-CoA metabolism, gene regulation, acyl-CoA-mediated cell signaling, transport-mediated lipid synthesis, membrane trafficking and also, ACBPs were indicated as a possible inhibitor of diazepam binding to the GABA-A receptor. To estimate the importance of the non-specific electrostatic energy in the ACBP-membrane interaction, we computationally modeled the interaction of HgACBP with both anionic and neutral membranes. To compute the Free Electrostatic Energy of Binding (dE), we used the Finite Difference Poisson Boltzmann Equation (FDPB) method as implemented in APBS. In the most energetically favorable orientation, ACBP brings charged residues Lys18 and Lys50 and hydrophobic residues Met46 and Leu47 into membrane surface proximity. This conformation suggests that these four ACBP amino acids are most likely to play a leading role in the ACBP-membrane interaction and ligand intake. Thus, we propose that long range electrostatic forces are the first step in the interaction mechanism between ACBP and membranes.

A hydrophobic loop in acyl-CoA binding protein is functionally important for binding to palmitoyl-coenzyme A: a molecular dynamics study.

Acyl-CoA binding protein (ACBP) plays a key role in lipid metabolism, interacting via a partly unknown mechanism with high affinity with long chain fatty acyl-CoAs (LCFA-CoAs). At present there is no study of the microscopic way ligand binding is accomplished. We analyzed this process by molecular dynamics (MDs) simulations. We proposed a computational model of ligand, able to reproduce some evidence from nuclear magnetic resonance (NMR) data, quantitative time resolved fluorometry and X-ray crystallography. We found that a hydrophobic loop, not in the active site, is important for function. Besides, multiple sequence alignment shows hydrophobicity (and not the residues itselves) conservation.

Structure of armadillo ACBP: a new member of the acyl-CoA-binding protein family.

The X-ray structure of the tetragonal form of apo acyl-CoA-binding protein (ACBP) from the Harderian gland of the South American armadillo Chaetophractus villosus has been solved. ACBP is a carrier for activated long-chain fatty acids and has been associated with many aspects of lipid metabolism. Its secondary structure is highly similar to that of the corresponding form of bovine ACBP and exhibits the unique flattened alpha-helical bundle (up-down-down-up) motif reported for animal, yeast and insect ACBPs. Conformational differences are located in loops and turns, although these structural differences do not suffice to account for features that could be related to the unusual biochemistry and lipid metabolism of the Harderian gland.

Purification, crystallization and preliminary X-ray analysis of Triatoma virus (TrV) from Triatoma infestans.

Triatoma virus (TrV) is a viral pathogen of the blood-sucking reduviid bug Triatoma infestans, the most important vector of American human trypanosomiasis (Chagas' disease). TrV has been putatively classified as a member of the Cripavirus genus (type cricket paralysis virus) in the Dicistroviridae family. This work describes the purification of TrV particles from infected T. infestans and their crystallization and preliminary crystallographic analyses. Two different crystal forms, rhombohedral and orthorhombic, were obtained at room temperature by the hanging-drop vapour-diffusion technique using polyethylene glycol and polyethylene glycol monomethylether as precipitants. The rhombohedral crystals have unit-cell parameters a = b = 306.6, c = 788.4 A (hexagonal setting), diffract to 3.2 A resolution and contain one-third of the viral particle per asymmetric unit. The orthorhombic crystals have cell parameters a = 336, b = 351, c = 332 A, diffract to about 2.5 A resolution, and contain one-half of a virus particle in the asymmetric unit. A complete diffraction data set has been collected to 3.2 A resolution, using synchrotron radiation, from a single rhombohedral crystal under cryogenic conditions.