RESUMO
Oxidative stress caused by pollution and lifestyle changes causes an excess of free radicals that react chemically with cell constituents leading to irreversible damage. There are molecules known as antioxidants that reduce the levels of free radicals. Some pigments of fruits and vegetables known as anthocyanins have antioxidant properties. Their interaction with the cell membrane becomes a crucial step in studying these substances. In this research, molecular dynamics simulations, particularly, coarse-grained molecular dynamics (CGMD) were used. This technique aims to replace functional groups with corresponding beads that represent their level of polarity and affinities to other chemical groups. Also, umbrella sampling was carried out to obtain free energy profiles that describe well the orientation and location of antioxidants in a membrane considering Trolox, Cyanidin, Gallic Acid, and Resveratrol molecules to study the structural effects they cause on it. It was concluded in this study that an antioxidant when crossing the membrane does not cause either damage to the structural properties or the loss of packing and stratification of phospholipids. it was also observed that the most reactive part of the molecules could easily approach area A prone to lipid oxidation, which can describe the antioxidant capacity of these molecules.Communicated by Ramaswamy H. Sarma.
RESUMO
Protein-protein association events are involved in many physiological and pathological processes. Cataract disease is a pathology that manifests protein aggregation of crystallins. ß-Crystallins are present in a high proportion in the eye lens. Therefore, the structural study of the dimerization properties of crystallins can shed light on the first stages of protein aggregation. In the present work, we examine the protein-protein association profiles of the human ßB2-crystallin by employing extensive coarse-grained molecular dynamics (CG-MD) and the Markov state analysis. Interestingly, our results clearly show important changes in the protein dimerization kinetics between wt-HßB2C and the deamidated systems. The two systems show dimeric conformations. However, the association and dissociation rates are very different. Our results show that the deamidated system can associate faster and dissociate slower than the wt- HßB2C. The deamidated system is in a slightly opened conformation with the Greek-key motifs well folded, suggesting that a complete unfolding of the protein is not required for aggregation. Our results describe the first stages of crystallin aggregation due to post-translational modifications.
RESUMO
The small soluble aggregates of Aß1-42 are broadly documented as potential targets for the development of new compounds with the capacity to inhibit the early stages of Alzheimer´s disease. Nevertheless, Aß1-42 peptides show an intrinsically disordered character with a high propensity for aggregation, which complicates the identification of conserved structural patterns. Because of this, experimental techniques find substantial difficulties in the characterization of such soluble oligomers. Theoretical techniques, such as molecular dynamics (MD) simulations, provide a possible workaround for this problem. However, the computational cost associated with comprehensively sampling the vast conformational space accessible to these peptides might become prohibitive. In this sense, coarse-grained (CG) simulations can effectively overcome that hurdle at a fraction of the computational cost. In this dataset, we furnish an extensive collection of Aß1-42 peptides in dimeric conformation generated with the SIRAH force field for CG MD simulation. It comprises 25 independent trajectories in .xtc (gromacs) format of Aß1-42 couples of peptides that evolve towards dimeric states along eleven µs-long unbiased simulations. Thanks to the backmapping capabilities of our force field, pseudo atomistic coordinates can be straightforwardly recovered from MD trajectories reported here and analyzed with popular molecular editing programs. This set of simulations performed at room conditions and physiological salt concentrations may furnish a complete collection of inter-peptide interfaces that can be used in high-throughput docking or as new starting states for peptide oligomerization seeding of Aß1-42 dimerization.
RESUMO
Oxidative stress plays an essential role in the regulation of vital processes in living organisms. Reactive oxygen species can react chemically with the constituents of the cells leading to irreversible damage. The first structure of the cell in contact with the environment that surrounds it is the membrane, which protects it and allows the exchange of substances. Some signals manifest when the components of a bilayer are undergoing oxidation, like an increase in the lipid area, decrease in the thickness of the bilayer, and exchange of the oxidized groups toward the bilayer surface. In this investigation, a molecular dynamics simulation was done on a set of Dioleoylphosphatidylcholine membranes with different percentage of oxidized lipids, in order to observe the effect of the oxidation degree on the membrane structure. It was found that, as higher the concentration of oxidized lipids is, the larger the damage of the membrane. This is reflected in the increase in the lipid area and the decrease in the thickness and membrane packing. Also, it was observed that hydrophobicity inside the membrane decreases as the oxidation percentage increases.Communicated by Ramaswamy H. Sarma.
Assuntos
Bicamadas Lipídicas , Simulação de Dinâmica Molecular , Membrana Celular/química , Interações Hidrofóbicas e Hidrofílicas , Bicamadas Lipídicas/química , OxirreduçãoRESUMO
This dataset contains a collection of molecular dynamics (MD) simulations of polyglutamine (polyQ) and glutamine-rich (Q-rich) peptides in the multi-microsecond timescale. Primary data from coarse-grained simulations performed using the SIRAH force field has been processed to provide fully atomistic coordinates. The dataset encloses MD trajectories of polyQs of 4 (Q4), 11 (Q11), and 36 (Q36) amino acids long. In the case of Q11, simulations in presence of Q5 and QEQQQ peptides, which modulate aggregation, are also included. The dataset also comprises MD trajectories of the gliadin related p31-43 peptide, and Insulin's C-peptide at pH=7 and pH=3.2, which constitute examples of Q-rich and Q-poor aggregating peptides. The dataset grants molecular insights on the role of glutamines in spontaneous and unbiased ab-initio aggregation of a series of peptides using a homogeneous set of simulations [1]. The trajectory files are provided in Protein Data Bank (PDB) format containing the Cartesian coordinates of all heavy atoms in the aggregating peptides. Further analyses of the trajectories can be performed directly using any molecular visualization/analysis software suites.
RESUMO
Poly glutamine and glutamine-rich peptides play a central role in a plethora of pathological aggregation events. However, biophysical characterization of soluble oligomers -the most toxic species involved in these processes- remains elusive due to their structural heterogeneity and dynamical nature. Here, we exploit the high spatio-temporal resolution of coarse-grained simulations as a computational microscope to characterize the aggregation propensity and morphology of a series of polyglutamine and glutamine-rich peptides. Comparative analysis of ab-initio aggregation pinpointed a double role for glutamines. In the first phase, glutamines mediate seeding by pairing monomeric peptides, which serve as primers for higher-order nucleation. According to the glutamine content, these low molecular-weight oligomers may then proceed to create larger aggregates. Once within the aggregates, buried glutamines continue to play a role in their maturation by optimizing solvent-protected hydrogen bonds networks.
RESUMO
CaMKII is a protein kinase whose function is regulated by the binding of the Calcium/Calmodulin complex (Ca2+/CaM). It is a major player in the Long Term Potentiation process where it acts as a molecular switch, oscillating between inhibited and active conformations. The mechanism for the switching is thought to be initiated by Ca2+/CaM binding, which allows the trans-phosphorylation of a subunit of CaMKII by a neighboring kinase, leading to the active state of the system. A combination of all-atom and coarse-grained MD simulations with free energy calculations, led us to reveal an interplay of electrostatic forces exerted by Ca2+/CaM on CaMKII, which initiate the activation process. The highly electrically charged Ca2+/CaM neutralizes basic regions in the linker domain of CaMKII, facilitating its opening and consequent activation. The emerging picture of CaMKII's behavior highlights the preponderance of electrostatic interactions, which are modulated by the presence of Ca2+/CaM and the phosphorylation of key sites.
Assuntos
Proteína Quinase Tipo 2 Dependente de Cálcio-Calmodulina , Cálcio , Cálcio/metabolismo , Proteína Quinase Tipo 2 Dependente de Cálcio-Calmodulina/metabolismo , Calmodulina , Fosforilação , Ligação ProteicaRESUMO
The detection of small molecules in living cells using genetically encoded FRET sensors has revolutionized our understanding of signaling pathways at the sub-cellular level. However, engineering fluorescent proteins and specific binding domains to create new sensors remains challenging because of the difficulties associated with the large size of the polypeptides involved, and their intrinsically huge conformational variability. Indeed, FRET sensors' design still relies on vague structural notions, and trial and error combinations of linkers and protein modules. We recently designed a FRET sensor for the second messenger cAMP named CUTie (Cyclic nucleotide Universal Tag for imaging experiments), which granted sub-micrometer resolution in living cells. Here we apply a combination of sequence/structure analysis to produce a new-generation FRET sensor for the second messenger cGMP based on Protein kinase G I (PKGI), which we named CUTie2. Coarse-grained molecular dynamics simulations achieved an exhaustive sampling of the relevant spatio-temporal coordinates providing a quasi-quantitative prediction of the FRET efficiency, as confirmed by in vitro experiments. Moreover, biochemical characterization showed that the cGMP binding module maintains virtually the same affinity and selectivity for its ligand thant the full-length protein. The computational approach proposed here is easily generalizable to other allosteric protein modules, providing a cost effective-strategy for the custom design of FRET sensors.
RESUMO
The Tumor Necrosis Factor (TNF) and the TNF receptor (TNFR) superfamilies are composed of 19 ligands and 30 receptors, respectively. The oligomeric properties of ligands, both membrane bound and soluble, has been studied most. However, less is known about the oligomeric properties of TNFRs. Earlier reports identified the extracellular, membrane-distal, cysteine-rich domain as a pre-ligand assembly domain which stabilizes receptor dimers and/or trimers in the absence of ligand. Nevertheless, recent reports based on structural nuclear magnetic resonance (NMR) highlight the intrinsic role of the transmembrane domains to form dimers (p75NTR), trimers (Fas), or dimers of trimers (DR5). Thus, understanding the structural basis of transmembrane oligomerization may shed light on the mechanism for signal transduction and the impact of disease-associated mutations in this region. To this end, here we used an in silico coarse grained molecular dynamics approach with Martini force field to study TNFR transmembrane homotypic interactions. We have first validated this approach studying the three TNFR described by NMR (p75NTR, Fas, and DR5). We have simulated membrane patches composed of 36 helices of the same receptor equidistantly distributed in order to get unbiassed information on spontaneous proteins assemblies. Good agreement was found in the specific residues involved in homotypic interactions and we were able to observe dimers, trimers, and higher-order oligomers corresponding to those reported in NMR experiments. We have, applied this approach to study the assembly of disease-related mutations being able to assess their impact on oligomerization stability. In conclusion, our results showed the usefulness of coarse grained simulations with Martini force field to study in an unbiased manner higher order transmembrane oligomerization.
RESUMO
We present a new coarse-grained (CG) model of cholesterol (CHOL) for the electrostatic-based ELBA force field. A distinguishing feature of our CHOL model is that the electrostatics is modeled by an explicit point dipole which interacts through an ideal vacuum permittivity. The CHOL model parameters were optimized in a systematic fashion, reproducing the electrostatic and nonpolar partitioning free energies of CHOL in lipid/water mixtures predicted by full-detailed atomistic molecular dynamics simulations. The CHOL model has been validated by comparison to structural, dynamic and thermodynamic properties with experimental and atomistic simulation reference data. The simulation of binary DPPC/cholesterol mixtures covering the relevant biological content of CHOL in mammalian membranes is shown to correctly predict the main lipid behavior as observed experimentally.
Assuntos
Colesterol/química , Simulação de Dinâmica Molecular , Membrana Celular/química , Difusão , Bicamadas Lipídicas/química , Estrutura Molecular , Eletricidade Estática , TermodinâmicaRESUMO
In the present work we investigated the differential interactions of the antimicrobial peptides (AMPs) aurein 1.2 and maculatin 1.1 with a bilayer composed of a mixture of the lipids 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (POPG) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE). We carried out molecular dynamics (MD) simulations using a coarse-grained approach within the MARTINI force field. The POPE/POPG mixture was used as a simple model of a bacterial (prokaryotic cell) membrane. The results were compared with our previous findings for structures of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), a representative lipid of mammalian cells. We started the simulations of the peptide-lipid system from two different initial conditions: peptides in water and peptides inside the hydrophobic core of the membrane, employing a pre-assembled lipid bilayer in both cases. Our results show similarities and differences regarding the molecular behavior of the peptides in POPE/POPG in comparison to their behavior in a POPC membrane. For instance, aurein 1.2 molecules can adopt similar pore-like structures on both POPG/POPE and POPC membranes, but the peptides are found deeper in the hydrophobic core in the former. Maculatin 1.1 molecules, in turn, achieve very similar structures in both kinds of bilayers: they have a strong tendency to form clusters and induce curvature. Therefore, the results of this study provide insight into the mechanisms of action of these two peptides in membrane leakage, which allows organisms to protect themselves against potentially harmful bacteria. Graphical Abstract Aurein pore structure (green) in a lipid bilayer composed by POPE (blue) and POPG (red) mixture. It is possible to see water beads (light blue) inside the pore.
Assuntos
Proteínas de Anfíbios/química , Antibacterianos/química , Peptídeos Catiônicos Antimicrobianos/química , Bicamadas Lipídicas/química , Fosfatidiletanolaminas/química , Fosfatidilgliceróis/química , Sequência de Aminoácidos , Proteínas de Anfíbios/isolamento & purificação , Proteínas de Anfíbios/metabolismo , Animais , Antibacterianos/isolamento & purificação , Antibacterianos/metabolismo , Peptídeos Catiônicos Antimicrobianos/isolamento & purificação , Peptídeos Catiônicos Antimicrobianos/metabolismo , Anuros/metabolismo , Sítios de Ligação , Interações Hidrofóbicas e Hidrofílicas , Bicamadas Lipídicas/metabolismo , Simulação de Dinâmica Molecular , Fosfatidilcolinas/química , Fosfatidilcolinas/metabolismo , Fosfatidiletanolaminas/metabolismo , Fosfatidilgliceróis/metabolismo , Ligação Proteica , Conformação Proteica em alfa-Hélice , Domínios e Motivos de Interação entre ProteínasRESUMO
Amphiphilic Janus dendrimers (JDs) are repetitively branched molecules with hydrophilic and hydrophobic components that self-assemble in water to form a variety of morphologies, including vesicles analogous to liposomes with potential pharmaceutical and medical application. To date, the self-assembly of JDs has not been fully investigated thus it is important to gain insight into its mechanism and dependence on JDs' molecular structure. In this study, the aggregation behavior in water of a second-generation bis-MPA JD was evaluated using experimental and computational methods. Dispersions of JDs in water were carried out using the thin-film hydration and ethanol injection methods. Resulting assemblies were characterized by dynamic light scattering, confocal microscopy, and atomic force microscopy. Furthermore, a coarse-grained molecular dynamics (CG-MD) simulation was performed to study the mechanism of JDs aggregation. The obtaining of assemblies in water with no interdigitated bilayers was confirmed by the experimental characterization and CG-MD simulation. Assemblies with dendrimersome characteristics were obtained using the ethanol injection method. The results of this study establish a relationship between the molecular structure of the JD and the properties of its aggregates in water. Thus, our findings could be relevant for the design of novel JDs with tailored assemblies suitable for drug delivery systems.
Assuntos
Dendrímeros/química , Simulação de Dinâmica Molecular , Água/química , Interações Hidrofóbicas e Hidrofílicas , Estrutura MolecularRESUMO
In this work, we developed a coarse-grained model of sumatriptan suitable for extensive molecular dynamics simulations. First, we confirmed the interfacial distribution of this drug in bilayers through cryogenic transmission electron microscopy and small-angle X-ray scattering techniques, as was predicted by our previous atomistic simulations. Based on these simulations, we developed a coarse-grained model for sumatriptan able to reproduce its overall molecular behavior, captured by atomistic simulations and experiments. We then tested the sumatriptan model in a micellar environment along with experimental characterization of sumatriptan-loaded micelles. The simulation results showed good agreement with photon correlation spectroscopy and electrophoretic mobility experiments performed in this work. The particle size of the obtained micelles was comparable with the simulated ones; meanwhile, zeta-potential results suggest adsorption of the drug on the micellar surface. This model is a step forward in the search for a suitable drug-delivery system for sumatriptan.
Assuntos
Simulação de Dinâmica Molecular , Sumatriptana/química , Bicamadas Lipídicas/química , Lipossomos/química , Micelas , Microscopia Eletrônica , Conformação Molecular , Poloxâmero/química , Espalhamento a Baixo Ângulo , Difração de Raios XRESUMO
Coarse-grained simulation schemes are increasingly gaining popularity in the scientific community because of the significant speed up granted, allowing a considerable expansion of the accessible time and size scales accessible to molecular simulations. However, the number of compatible force fields capable of representing ensembles containing different molecular species (i.e., Protein, DNA, etc) is still limited. Here, we present a set of parameters and simplified representation for lipids compatible with the SIRAH force field for coarse-grained simulations ( http://www.sirahff.com ). We show that the present model not only achieves a correct reproduction of structural parameters as area per lipid and thickness, but also dynamic descriptors such as diffusion coefficient, order parameters, and proper temperature driven variations. Adding phospholipid membranes to the existing aqueous solution, protein and DNA representations of the SIRAH force field permit considering the most common problems tackled by the biomolecular simulation community.
Assuntos
Dimiristoilfosfatidilcolina , Bicamadas Lipídicas/química , Simulação de Dinâmica Molecular , DNA/química , Proteínas/químicaRESUMO
Coarse-grained (CG) models allow enlarging the size and time scales that are reachable by atomistic molecular dynamics simulations. A CG force field (FF) for lipids and amino acids that possesses a polarizable water model has been developed following the MARTINI parametrization strategy, the BMW-MARTINI [1]. We tested the BMW-MARTINI FF capability to describe some structural and thermodynamical properties of lipid monolayers and bilayers. We found that, since the surface tension values of oil/water interfaces calculated with the model are not correct, compression isotherms of lipid monolayers present artifacts. Also, this FF predicts DPPC and DAPC bilayers to remain in the Lα phase at temperatures as low as 283K, contrary to the expected from their experimental Tm values. Finally, simulations at constant temperature of bilayers of saturated lipids belonging to PC homologous, showed an increase in the mean molecular area (Mma) upon increasing the chain length, inversely to the experimental observation. We refined BMW-MARTINI FF by modifying as few parameters as possible in order to bring simulated and experimental measurements closer. We have also modified structural parameters of the lipid geometry that do not have direct influence in global properties of the bilayer membranes or monolayers, but serve to approach the obtained CG geometry to atomistic reference values. The refined FF is able to better reproduce phase transition temperatures and Mma for saturated PC bilayers than BMW-MARTINI and MARTINI FF. Finally, the simulated surface pressure-Mma isotherms of PC monolayers resemble the experimental ones and eliminate serious artifacts of previous models.
Assuntos
1,2-Dipalmitoilfosfatidilcolina/química , Bicamadas Lipídicas/química , Simulação de Dinâmica Molecular , Fosfatidilcolinas/química , Transição de Fase , Tensão Superficial , Temperatura , Termodinâmica , Água/químicaRESUMO
Biological membranes and model lipid systems containing high amounts of unsaturated lipids and sterols are subject to chemical and/or photo-induced lipid oxidation, which leads to the creation of exotic oxidized lipid products (OxPLs). OxPLs are known to have significant physiological impact in cellular systems and also affect physical properties of both biological and model lipid bilayers. In this paper we (i) provide a perspective on the existing literature on simulations of lipid bilayer systems containing oxidized lipid species as well as the main related experimental results, (ii) describe our new data of all-atom and coarse-grained simulations of hydroperoxidized lipid monolayer and bilayer systems and (iii) provide a comparison of the MARTINI and ELBA coarse grained force fields for lipid bilayer systems. We show that the better electrostatic treatment of interactions in ELBA is able to resolve previous conflicts between experiments and simulations. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.