ABSTRACT
Translocator protein (TSPO) is an 18 kDa transmembrane protein, localized primarily on the outer mitochondrial membrane. It has been found to be involved in various physiological processes and pathophysiological conditions. Though studies on its structure have been performed only recently, there is little information on the nature of dynamics and doubts about some structures referenced in the literature, especially the NMR structure of mouse TSPO. In the present work, we thoroughly study the dynamics of mouse TSPO protein by means of atomistic molecular dynamics simulations, in presence as well as in absence of the diagnostic ligand PKA. We considered two starting structures: the NMR structure and a homology model (HM) generated on the basis of X-ray structures from bacterial TSPO. We examine the conformational landscape in both the modes for both starting points, in presence and absence of the ligand, in order to measure its impact for both structures. The analysis highlights high flexibility of the protein globally, but NMR simulations show a surprisingly flexibility even in the presence of the ligand. Interestingly, this is not the case for HM calculations, to the point that the ligand seems not so stable as in the NMR system and an unbinding event process is partially sampled. All those results tend to show that the NMR structure of mTSPO seems not deficient but is just in another portion of the global conformation space of TSPO.
ABSTRACT
We report here an instructive case referred at 16 months-old for exploration of hemolysis without anemia (compensated anemia with reticulocytosis). The biology tests confirmed the hemolysis with increased total and indirect bilirubin. The usual hemolysis diagnosis tests were normal (DAT, G6PD, PK, Hb electrophoresis) except cytology and ektacytometry suggesting an association of multiple red blood cell (RBC) membrane disorders. This led us to propose a molecular screening analysis using targeted-Next Generation Sequencing (t-NGS) with a capture technique on 93 genes involved in RBC and erythropoiesis defects. We identified 4 missense heterozygous allelic variations, all of them were described without any significance (VUS) in the SLC4A1, RhAG, PIEZO1 and SPTB genes. The study of the familial cosegregation and research functional tests allowed to decipher the role of at least two by two genes in the phenotype and the hemolytic disease of this young patient. Specialized t-NGS panel (or virtual exome/genome sequencing) in a disease-referent laboratory and the motivated collaboration of clinicians, biologists and scientists should be the gold standard for improving the diagnosis of the patients affected with RBC diseases or rare inherited anemias.
Subject(s)
Hematologic Diseases , Spherocytosis, Hereditary , Humans , Spherocytosis, Hereditary/diagnosis , Spherocytosis, Hereditary/genetics , Spectrin/genetics , High-Throughput Nucleotide Sequencing , Hemolysis , Mutation , Erythrocytes , Phenotype , Anion Exchange Protein 1, Erythrocyte/genetics , Ion Channels/geneticsABSTRACT
Conformational flexibility plays an essential role in antibodies' functional and structural stability. They facilitate and determine the strength of antigen-antibody interactions. Camelidae express an interesting subtype of single-chain antibody, named Heavy Chain only Antibody. They have only one N-terminal Variable domain (VHH) per chain, composed of Frameworks (FRs) and Complementarity Determining regions (CDRs) like their VH and VL counterparts in IgG. Even when expressed independently, VHH domains display excellent solubility and (thermo)stability, which helps them to retain their impressive interaction capabilities. Sequence and structural features of VHH domains contributing to these abilities have already been studied compared to classical antibodies. To have the broadest view and understand the changes in dynamics of these macromolecules, large-scale molecular dynamics simulations for a large number of non-redundant VHH structures have been performed for the first time. This analysis reveals the most prevalent movements in these domains. It reveals the four main classes of VHHs dynamics. Diverse local changes were observed in CDRs with various intensities. Similarly, different types of constraints were observed in CDRs, while FRs close to CDRs were sometimes primarily impacted. This study sheds light on the changes in flexibility in different regions of VHH that may impact their in silico design.
Subject(s)
Camelidae , Immunoglobulin Variable Region , Animals , Immunoglobulin Variable Region/chemistry , Complementarity Determining Regions/chemistry , Immunoglobulin Heavy Chains/chemistry , Molecular Dynamics SimulationABSTRACT
Zika virus (ZIKV) from Uganda (UG) expresses a phenotype related to fetal loss, whereas the variant from Brazil (BR) induces microcephaly in neonates. The differential virulence has a direct relation to biomolecular mechanisms that make one strain more aggressive than the other. The nonstructural protein 1 (NS1) is a key viral toxin to comprehend these viral discrepancies because of its versatility in many processes of the virus life cycle. Here, we aim to examine through coarse-grained models and molecular dynamics simulations the protein-membrane interactions for both NS1ZIKV-UG and NS1ZIKV-BR dimers. A first evaluation allowed us to establish that the NS1 proteins, in the membrane presence, explore new conformational spaces when compared to systems simulated without a lipid bilayer. These events derive from both differential coupling patterns and discrepant binding affinities to the membrane. The N-terminal domain, intertwined loop, and greasy finger proposed previously as binding membrane regions were also computationally confirmed by us. The anchoring sites have aromatic and ionizable residues that manage the assembly of NS1 toward the membrane, especially for the Ugandan variant. Furthermore, in the presence of the membrane, the difference in the dynamic cross-correlation of residues between the two strains is particularly pronounced in the putative epitope regions. This suggests that the protein-membrane interaction induces changes in the distal part related to putative epitopes. Taken together, these results open up new strategies for the treatment of flaviviruses that would specifically target these dynamic differences.
Subject(s)
Zika Virus Infection , Zika Virus , Humans , Virulence/genetics , Viral Nonstructural Proteins/chemistry , Antibodies, Viral , EpitopesABSTRACT
BACKGROUND: An antibody directed against a high-prevalence red blood cell (RBC) antigen was detected in a 67-year-old female patient of North African ancestry with a history of a single pregnancy and blood transfusion. So far, the specificity of the proband's alloantibody remained unknown in our immunohematology reference laboratory. STUDY DESIGN AND METHODS: Whole-exome sequencing (WES) was performed on the proband's DNA. The reactivity to the SLC29A1-encoded ENT1 adenosine transporter was investigated by flow cytometry analyses of ENT1-expressing HEK293 cells, and RBCs from Augustine-typed individuals. Erythrocyte protein expression level, nucleoside-binding capacity, and molecular structure of the proband's ENT1 variant were further explored by western blot, flow cytometry, and molecular dynamics calculations, respectively. RESULTS: A missense variant was identified in the SLC29A1 gene, which encodes the Augustine blood group system. It arises from homozygosity for a rare c.242A > G missense mutation that results in a nonsynonymous p.Asn81Ser substitution within the large extracellular loop of ENT1. Flow cytometry analyses demonstrated that the proband's antibody was reactive against HEK-293 cells transfected with control but not proband's SLC29A1 cDNA. Consistent with this finding, proband's antibody was found to be reactive with At(a-) (AUG:-2), but not AUG:-1 (null phenotype) RBCs. Data from structural analysis further supported that the proband's p.Asn81Ser variation does not alter ENT1 binding of its specific inhibitor NBMPR. CONCLUSION: Our study provides evidence for a novel high-prevalence antigen, AUG4 (also called ATAM after the proband's name) in the Augustine blood group system, encoded by the rare SLC29A1 variant allele AUG*04 (c.242A > G, p.Asn81Ser).
Subject(s)
Blood Group Antigens , Pregnancy , Female , Humans , HEK293 Cells , Prevalence , Blood Group Antigens/genetics , Isoantibodies , Molecular StructureABSTRACT
Heavy Chain Only Antibodies are specific to Camelid species. Despite the lack of the light chain variable domain, their heavy chain variable domain (VH) domain, named VHH or nanobody, has promising potential applications in research and therapeutic fields. The structural study of VHH is therefore of great interest. Unfortunately, considering the huge amount of sequences that might be produced, only about one thousand of VHH experimental structures are publicly available in the Protein Data Bank, implying that structural model prediction of VHH is a necessary alternative to obtaining 3D information besides its sequence. The present study aims to assess and compare the quality of predictions from different modelling methodologies. Established comparative & homology modelling approaches to recent Deep Learning-based modelling strategies were applied, i.e. Modeller using single or multiple structural templates, ModWeb, SwissModel (with two evaluation schema), RoseTTAfold, AlphaFold 2 and NanoNet. The prediction accuracy was evaluated using RMSD, TM-score, GDT-TS, GDT-HA and Protein Blocks distance metrics. Besides the global structure assessment, we performed specific analyses of Frameworks and CDRs structures. We observed that AlphaFold 2 and especially NanoNet performed better than the other evaluated softwares. Importantly, we performed molecular dynamics simulations of an experimental structure and a NanoNet predicted model of a VHH in order to compare the global structural flexibility and local conformations using Protein Blocks. Despite rather similar structures, substantial differences in dynamical properties were observed, which underlies the complexity of the task of model evaluation.Communicated by Ramaswamy H. Sarma.
Subject(s)
Immunoglobulin Heavy Chains , Immunoglobulin Variable Region , Immunoglobulin Variable Region/chemistry , Immunoglobulin Heavy Chains/chemistryABSTRACT
Flaviviruses comprise a large group of arboviral species that are distributed in several countries of the tropics, neotropics, and some temperate zones. Since they can produce neurological pathologies or vascular damage, there has been intense research seeking better diagnosis and treatments for their infections in the last decades. The flavivirus NS1 protein is a relevant clinical target because it is involved in viral replication, immune evasion, and virulence. Being a key factor in endothelial and tissue-specific modulation, NS1 has been largely studied to understand the molecular mechanisms exploited by the virus to reprogram host cells. A central part of the viral maturation processes is the NS1 oligomerization because many stages rely on these protein-protein assemblies. In the present study, the self-associations of NS1 proteins from Zika, Dengue, and West Nile viruses are examined through constant-pH coarse-grained biophysical simulations. Free energies of interactions were estimated for different oligomeric states and pH conditions. Our results show that these proteins can form both dimers and tetramers under conditions near physiological pH even without the presence of lipids. Moreover, pH plays an important role mainly controlling the regimes where van der Waals interactions govern their association. Finally, despite the similarity at the sequence level, we found that each flavivirus has a well-characteristic protein-protein interaction profile. These specific features can provide new hints for the development of binders both for better diagnostic tools and the formulation of new therapeutic drugs.
Subject(s)
Flavivirus , West Nile virus , Zika Virus Infection , Zika Virus , Humans , Viral Nonstructural Proteins/metabolism , Virus Replication , Zika Virus/metabolismABSTRACT
VHH, i.e., VH domains of camelid single-chain antibodies, are very promising therapeutic agents due to their significant physicochemical advantages compared to classical mammalian antibodies. The number of experimentally solved VHH structures has significantly improved recently, which is of great help, because it offers the ability to directly work on 3D structures to humanise or improve them. Unfortunately, most VHHs do not have 3D structures. Thus, it is essential to find alternative ways to get structural information. The methods of structure prediction from the primary amino acid sequence appear essential to bypass this limitation. This review presents the most extensive overview of structure prediction methods applied for the 3D modelling of a given VHH sequence (a total of 21). Besides the historical overview, it aims at showing how model software programs have been shaping the structural predictions of VHHs. A brief explanation of each methodology is supplied, and pertinent examples of their usage are provided. Finally, we present a structure prediction case study of a recently solved VHH structure. According to some recent studies and the present analysis, AlphaFold 2 and NanoNet appear to be the best tools to predict a structural model of VHH from its sequence.
Subject(s)
Camelids, New World , Immunoglobulin Heavy Chains , Amino Acid Sequence , Animals , Antibodies , Immunoglobulin Heavy Chains/chemistry , Models, StructuralABSTRACT
The transmembrane glycoprotein CD36, which is responsible of the metabolic disorders, and the elevated intake of fat induces lipid buildup, is a multifunctional scavenger receptor signaling those functions in high-affinity tissue uptake of long-chain fatty acids. In this study, we used series of molecular dynamics simulations of the wild type and mutants types K164A CD36 protein interacting with one palmitic acid (PLM) besides simulations of the wild type interacting with the three PLM to find out the mechanism of the functioning of the complex CD36/Fatty acids and the unraveling of the role of the mutation. Additionally we determined whether Lys164, mostly exposed to protein surface, played important roles in fatty acid uptake. These simulations revealed, the conformational changes induced by Lys164 residue and the altered interactions induced by the mutagenesis of surface lysine that was badly influencing the folding, utility, solubility, and stability form of the variant. Furthermore, Lys164 residue provided the structural basis of forming an opening at the region of principal portal for the dissociation of palmitic acid. The results of our simulations revealed hole two fatty acids found in CD36 cavity structure and it was the most preferred to CD36 structure stabilization.
Subject(s)
CD36 Antigens/metabolism , Fatty Acids/metabolism , Biological Transport , CD36 Antigens/chemistry , CD36 Antigens/genetics , Humans , Molecular Dynamics Simulation , Palmitic Acid/metabolism , Point Mutation , Protein ConformationABSTRACT
Viruses can impact and affect human populations in a severe way. The appropriate differentiation among several species or strains of viruses is one of the biggest challenges for virology and infectiology studies. The detection of measurables-quantified discrepancies allows for more accurate clinical diagnoses and treatments for viral diseases. In the present study, we have used a computational approach to explore the dynamical properties of the nonstructural protein 1 from two strains of Zika virus. Our results show that despite a high sequence similarity, the two viral proteins from different origins can exhibit significant dissimilar structural dynamics, which complement their reported differential virulence. The present study opens up new ways in the understanding of the infectivity for these biological entities.
Subject(s)
Zika Virus Infection , Zika Virus , Humans , Viral Nonstructural Proteins , VirulenceABSTRACT
Translocator protein (TSPO), a mitochondrial membrane protein, has been extensively studied, and its role is still debated and continues to be enigmatic. From a structural perspective, despite availability of atomic structures from different species, the possible oligomeric state and its 3D structure remains elusive. In the present study, we attempted to study dynamics of TSPO from the perspective of oligomerization. In this aim, we examined if and how TSPO monomers could assemble to form a dimer. Accordingly, we performed several coarse-grained molecular dynamics simulations considering two different initial configurations, one with a pair of TSPO monomers distantly placed in a model of a bilayer composed of DMPC/cholesterol mixture and the other with preformed dimer models with different starting interactions. We identify stable TSPO dimers with diverse interfaces, some of which were consistent with earlier experimental observations on putative TSPO oligomer interfaces. For most of the stable ones, interactions between aromatic residues were significantly overrepresented in diverse oligomeric organizations. Interestingly, we identified different communication pathways that involve dimer interfaces. Additionally, we observed that cholesterol molecules in close interaction with the TSPO dimer were able to translocate through the bilayer. This phenomenon might be related to the putative mechanism of cholesterol transport and could be increased and favored by the dimer formation. Overall, our observations shed new light on TSPO oligomerization and bring new perspectives on its dynamics, as well its interactions with protein and ligand partners.
Subject(s)
Molecular Dynamics Simulation , Receptors, GABA , Carrier Proteins , Cholesterol , Dimerization , Receptors, GABA/metabolismABSTRACT
Viruses are enthusiastically studied due to the great impact that these organisms can have on human health. Computational approaches can contribute offering tools that can shed light on important molecular mechanisms that help to design new diagnostic procedures. Several cellular processes between the immune-host system and the pathogenic organism are dependent on specific intermolecular interactions. In this study, we evaluated theoretical approaches to understand some properties of the antigen-antibody interactions considering the titratable properties of all ionizable residues of the nonstructural viral protein 1 (NS1) of the West Nile virus (WNV) and the Zika virus (ZIKV). Constant-pH Monte Carlo simulations were performed to estimate electrostatic properties such as the pKa shifts (ΔpKa). We proposed an alternative criterion for the discrimination of antigenic residues based on ΔpKas. Our outcomes were analyzed by an evaluation of the sensitivity and specificity through a receiver operating characteristic (ROC). As a starting point, we used the known crystallographic structure for the complex of NS1WNV(176-352) and the specific antibody 22NS1 (PDB ID 4OII ) to differentiate the residues belonging to that interface. With an optimal threshold for the absolute value of the pKa shifts, we found that is possible to predict antigenic epitopes reproducing the interfaces as defined by the X-ray structure. After this validation, we evaluated theoretical predictions based on protein-protein (PP) complexation simulations. From them, we observe amino acids with an antigenic potential and defined the optimum threshold that was applied for two strains of ZIKV (i.e., Uganda and Brazil). Several ionizable residues with antigenic capacity were identified. This is favorably related to some studies that show the high immunogenicity of secreted NS1. This approach opens up an important discussion about what are termed here "electrostatic epitopes" and how they work as an important reference in the paratope-epitope interaction for viral systems.
Subject(s)
Epitopes/chemistry , Epitopes/immunology , Flavivirus/immunology , Models, Molecular , Static Electricity , Antibodies, Viral/chemistry , Antibodies, Viral/immunology , Monte Carlo Method , Protein ConformationABSTRACT
Protein structures are highly dynamic macromolecules. This dynamics is often analysed through experimental and/or computational methods only for an isolated or a limited number of proteins. Here, we explore large-scale protein dynamics simulation to observe dynamics of local protein conformations using different perspectives. We analysed molecular dynamics to investigate protein flexibility locally, using classical approaches such as RMSf, solvent accessibility, but also innovative approaches such as local entropy. First, we focussed on classical secondary structures and analysed specifically how ß-strand, ß-turns, and bends evolve during molecular simulations. We underlined interesting specific bias between ß-turns and bends, which are considered as the same category, while their dynamics show differences. Second, we used a structural alphabet that is able to approximate every part of the protein structures conformations, namely protein blocks (PBs) to analyse (i) how each initial local protein conformations evolve during dynamics and (ii) if some exchange can exist among these PBs. Interestingly, the results are largely complex than simple regular/rigid and coil/flexible exchange. AbbreviationsNeqnumber of equivalentPBProtein BlocksPDBProtein DataBankRMSfroot mean square fluctuationsCommunicated by Ramaswamy H. Sarma.
Subject(s)
Molecular Dynamics Simulation , Proteins , Entropy , Protein Conformation , Protein Structure, Secondary , Proteins/geneticsABSTRACT
Flexibility is an intrinsic essential feature of protein structures, directly linked to their functions. To this day, most of the prediction methods use the crystallographic data (namely B-factors) as the only indicator of protein's inner flexibility and predicts them as rigid or flexible. PredyFlexy stands differently from other approaches as it relies on the definition of protein flexibility (i) not only taken from crystallographic data, but also (ii) from Root Mean Square Fluctuation (RMSFs) observed in Molecular Dynamics simulations. It also uses a specific representation of protein structures, named Long Structural Prototypes (LSPs). From Position-Specific Scoring Matrix, the 120 LSPs are predicted with a good accuracy and directly used to predict (i) the protein flexibility in three categories (flexible, intermediate and rigid), (ii) the normalized B-factors, (iii) the normalized RMSFs, and (iv) a confidence index. Prediction accuracy among these three classes is equivalent to the best two class prediction methods, while the normalized B-factors and normalized RMSFs have a good correlation with experimental and in silico values. Thus, PredyFlexy is a unique approach, which is of major utility for the scientific community. It support parallelization features and can be run on a local cluster using multiple cores.
Subject(s)
Molecular Dynamics Simulation , Protein Conformation , Proteins/chemistry , Databases, Protein , Datasets as Topic , SoftwareABSTRACT
Glucose plays a crucial role in the mammalian cell metabolism. In the erythrocytes and endothelial cells of the blood-brain barrier, glucose uptake is mediated by the glucose transporter type 1 (GluT1). GluT1 deficiency or mutations cause severe physiological disorders. GluT1 is also an important target in cancer therapy as it is overexpressed in tumor cells. Previous studies have suggested that GluT1 mediates solute transfer through a cycle of conformational changes. However, the corresponding 3D structures adopted by the transporter during the transfer process remain elusive. In the present work, we first elucidate the whole conformational landscape of GluT1 in the absence of glucose, using long molecular dynamics simulations and show that the transitions can be accomplished through thermal fluctuations. Importantly, we highlight a strong coupling between intracellular and extracellular domains of the protein that contributes to the transmembrane helices reorientation during the transition. The conformations adopted during the simulations differ from the known 3D bacterial homologs structures resolved in similar states. In holo state simulations, we find that glucose transits along the pathway through significant rotational motions, while maintaining hydrogen bonds with the protein. These persistent motions affect side chains orientation, which impacts protein mechanics and allows glucose progression.
Subject(s)
Glucose Transporter Type 1/metabolism , Glucose/metabolism , Biological Transport , Glucose Transporter Type 1/chemistry , Molecular Dynamics Simulation , Principal Component Analysis , Protein Domains , Protein Structure, Secondary , TemperatureABSTRACT
The interactions between viruses and actin cytoskeleton have been widely studied. We showed that rotaviruses remodel microfilaments in intestinal cells and demonstrated that this was due to the VP4 spike protein. Microfilaments mainly occur in the apical domain of infected polarized enterocytes and favor the polarized apical exit of viral progeny. The present work aims at the identification of molecular determinants of actin-VP4 interactions. We used various deletion mutants of VP4 that were transfected into Cos-7 cells and analyzed interactions by immunofluorescence confocal microscopy. It has been established that the C-terminal part of VP4 is embedded within viral particles when rotavirus assembles. The use of specific monoclonal antibodies demonstrated that VP4 is expressed in different forms in infected cells: classically as spike on the outer layer of virus particles, but also as free soluble protein in the cytosol. The C terminus of free VP4 was identified as interacting with actin microfilaments. The VP4 actin binding domain is unable to promote microfilament remodeling by itself; the coiled-coil domain is also required in this process. This actin-binding domain was shown to dominate a previously identified peroxisomal targeting signal, located in the three last amino acids of VP4. The newly identified actin-binding domain is highly conserved in rotavirus strains from species A, B, and C, suggesting that actin binding and remodeling is a general strategy for rotavirus exit. This provides a novel mechanism of protein-protein interactions, not involving cell signaling pathways, to facilitate rotavirus exit.IMPORTANCE Rotaviruses are causal agents of acute infantile viral diarrhea. In intestinal cells, in vitro as well as in vivo, virus assembly and exit do not imply cell lysis but rely on an active process in which the cytoskeleton plays a major role. We describe here a novel molecular mechanism by which the rotavirus spike protein VP4 drives actin remodeling. This relies on the fact that VP4 occurs in different forms. Besides its structural function within the virion, a large proportion of VP4 is expressed as free protein. Here, we show that free VP4 possesses a functional actin-binding domain. This domain, in coordination with a coiled-coil domain, promotes actin cytoskeleton remodeling, thereby providing the capacity to destabilize the cell membrane and allow efficient rotavirus exit.
Subject(s)
Actins/chemistry , Actins/metabolism , Capsid Proteins/chemistry , Capsid Proteins/metabolism , Rotavirus/metabolism , Animals , Binding Sites , COS Cells , Capsid Proteins/genetics , Cell Line , Chlorocebus aethiops , Microscopy, Confocal , Microscopy, Fluorescence , Mutation , Protein Binding , Protein Domains , Rotavirus/chemistry , Rotavirus/geneticsABSTRACT
The flavivirus genus has several organisms responsible for generating various diseases in humans. Recently, especially in tropical regions, Zika virus (ZIKV) has raised great health concerns due to the high number of cases affecting the area during the last years that has been accompanied by a rise in the cases of the Guillain-Barré syndrome and fetal and neonatal microcephaly. Diagnosis is still difficult since the clinical symptoms between ZIKV and other flaviviruses (e.g., dengue and yellow fever) are highly similar. The understanding of their common physicochemical properties that are pH-dependent and biomolecular interaction features and their differences sheds light on the relation strain-virulence and might suggest alternative strategies toward differential serological diagnostics and therapeutic intervention. Due to their immunogenicity, the primary focus of this study was on the ZIKV nonstructural proteins 1 (NS1). By means of computational studies and semiquantitative theoretical analyses, we calculated the main physicochemical properties of this protein from different strains that are directly responsible for the biomolecular interactions and, therefore, can be related to the differential infectivity of the strains. We also mapped the electrostatic differences at both the sequence and structural levels for the strains from Uganda to Brazil, which could suggest possible molecular mechanisms for the increase of the virulence of ZIKV in Brazil. Exploring the interfaces used by NS1 to self-associate in some different oligomeric states and interact with membranes and the antibody, we could map the strategy used by the ZIKV during its evolutionary process. This indicates possible molecular mechanisms that can be correlated with the different immunological responses. By comparing with the known antibody structure available for the West Nile virus, we demonstrated that this antibody would have difficulties to neutralize the NS1 from the Brazilian strain. The present study also opens up perspectives to computationally design high-specificity antibodies.
ABSTRACT
This chapter describes a protocol to establish a three-dimensional (3D) model of a protein and to explore its conformational landscape. It combines predictions from up-to-date bioinformatics methods with low-resolution experimental data. It also proposes to examine rapidly the dynamics of the protein using molecular dynamics simulations with a coarse-grained force field. Tools for analyzing these trajectories are suggested as well as those for constructing all-atoms models. Thus, starting from a protein sequence and using free software, the user can get important conformational information, which might improve the knowledge about the protein function.
Subject(s)
Cryoelectron Microscopy/methods , Receptors, GABA/chemistry , Amino Acid Sequence , Computational Biology/methods , Humans , Membrane Transport Proteins/chemistry , Models, Molecular , Molecular Dynamics Simulation , Protein Conformation , Receptors, GABA/genetics , SoftwareABSTRACT
Serotonin (5-hydroxytryptamine, 5-HT) is a well-known neurotransmitter that is involved in a growing number of functions in peripheral tissues. Recent studies have shown nonpharmacological functions of 5-HT linked to its chemical properties. Indeed, it was reported that 5-HT may, on the one hand, bind lipid membranes and, on the other hand, protect red blood cells through a mechanism independent of its specific receptors. To better understand these underevaluated properties of 5-HT, we combined biochemical, biophysical, and molecular dynamics simulations approaches to characterize, at the molecular level, the antioxidant capacity of 5-HT and its interaction with lipid membranes. To do so, 5-HT was added to red blood cells and lipid membranes bearing different degrees of unsaturation. Our results demonstrate that 5-HT acts as a potent antioxidant and binds with a superior affinity to lipids with unsaturation on both alkyl chains. We show that 5-HT locates at the hydrophobic-hydrophilic interface, below the glycerol group. This interfacial location is stabilized by hydrogen bonds between the 5-HT hydroxyl group and lipid headgroups and allows 5-HT to intercept reactive oxygen species, preventing membrane oxidation. Experimental and molecular dynamics simulations using membrane enriched with oxidized lipids converge to further reveal that 5-HT contributes to the termination of lipid peroxidation by direct interaction with active groups of these lipids and could also contribute to limit the production of new radicals. Taken together, our results identify 5-HT as a potent inhibitor of lipid peroxidation and offer a different perspective on the role of this pleiotropic molecule.
Subject(s)
Antioxidants/metabolism , Cell Membrane/metabolism , Membrane Lipids/metabolism , Serotonin/metabolism , Antioxidants/administration & dosage , Antioxidants/chemistry , Cell Membrane/chemistry , Erythrocytes/chemistry , Erythrocytes/metabolism , Flow Cytometry , Humans , Hydrogen Bonding , Hydrophobic and Hydrophilic Interactions , Lipid Bilayers/chemistry , Lipid Bilayers/metabolism , Lipid Peroxidation , Liposomes/chemistry , Liposomes/metabolism , Microscopy, Confocal , Molecular Dynamics Simulation , Oxidation-Reduction , Serotonin/administration & dosage , Serotonin/chemistryABSTRACT
Knowing the position of protein structures within the membrane is crucial for fundamental and applied research in the field of molecular biology. Only few web resources propose coordinate files of oriented transmembrane proteins, and these exclude predicted structures, although they represent the largest part of the available models. In this article, we present TMPL (http://www.dsimb.inserm.fr/TMPL/), a database of transmembrane protein structures (α-helical and ß-sheet) positioned in the lipid bilayer. It is the first database to include theoretical models of transmembrane protein structures, making it a large repository with more than 11 000 entries. The TMPL database also contains experimentally solved protein structures, which are available as either atomistic or coarse-grained models. A unique feature of TMPL is the possibility for users to update the database by uploading, through an intuitive web interface, the membrane assignments they can obtain with our recent OREMPRO web server.