Molecular dynamics simulations help to understand the complex behavior of molecules. The output of such a simulation describes the trajectories of individual atoms as snapshots of atom positions in time. Many compression schemes were developed to reduce the amount of data needed for storing long trajectories. This is achieved by limiting the precision of coordinates, encoding differences instead of absolute values, dimensionality reduction by principal component analysis, or by using polynomials approximating vertex trajectories. However, compression schemes using actual bonds between atoms have not been utilized to their full potential. Therefore, we developed a lossy compression method that captures the local, mostly rotational movement of atoms with respect to their bonded neighbors and predicts their positions in each frame. This allows full control over the data distortion. In our experiments, the method achieves data rates which are substantially better than the rates achieved by competing methods at the same error level.

• Klíčová slova
• 2010 MSC: 92C40, 68P30, Compression, Encoding, Graph traversal, Molecular dynamics, Molecular simulations, Trajectory,
• MeSH
• algoritmy * MeSH
• analýza hlavních komponent MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH

The delineation of protein-lipid interfaces is essential for understanding the mechanisms of various membrane-associated processes crucial to plant development and growth, including signalling, trafficking, and membrane transport. Due to their highly dynamic nature, the precise characterization of lipid-protein interactions by experimental techniques is challenging. Molecular dynamics simulations provide a powerful computational alternative with a spatial-temporal resolution allowing the atomistic-level description. In this review, we aim to introduce plant scientists to molecular dynamics simulations. We describe different steps of performing molecular dynamics simulations and provide a broad survey of molecular dynamics studies investigating plant protein-lipid interfaces. Our aim is also to illustrate that combining molecular dynamics simulations with artificial intelligence-based protein structure determination opens up unprecedented possibilities for future investigations of dynamic plant protein-lipid interfaces.

• Klíčová slova
• Integral membrane protein, membrane, molecular dynamics simulations, peripheral membrane protein, protein–lipid interactions, structural modelling,
• MeSH
• rostlinné proteiny * metabolismus   chemie   MeSH
• rostliny metabolismus   MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• přehledy MeSH
• Názvy látek
• rostlinné proteiny * MeSH

We report the first complete description of the molecular mechanisms behind the transition of the N-methyl-d-aspartate (NMDA) receptor from the state where the transmembrane domain (TMD) and the ion channel are in the open configuration to the relaxed unliganded state where the channel is closed. Using an aggregate of nearly 1 µs of unbiased all-atom implicit membrane and solvent molecular dynamics (MD) simulations we identified distinct structural states of the NMDA receptor and revealed functionally important residues (GluN1/Glu522, GluN1/Arg695, and GluN2B/Asp786). The role of the "clamshell" motion of the ligand binding domain (LBD) lobes in the structural transition is supplemented by the observed structural similarity at the level of protein domains during the structural transition, combined with the overall large rearrangement necessary for the opening and closing of the receptor. The activated and open states of the receptor are structurally similar to the liganded crystal structure, while in the unliganded receptor the extracellular domains perform rearrangements leading to a clockwise rotation of up to 45 degrees around the longitudinal axis of the receptor, which closes the ion channel. The ligand-induced rotation of extracellular domains transferred by LBD-TMD linkers to the membrane-anchored ion channel is responsible for the opening and closing of the transmembrane ion channel, revealing the properties of NMDA receptor as a finely tuned molecular machine.

• Klíčová slova
• NMDA receptor transition, glutamate receptor gating, molecular dynamics simulations, molecular modeling, open and closed state,
• MeSH
• krysa rodu Rattus MeSH
• receptory N-methyl-D-aspartátu chemie   metabolismus   MeSH
• simulace molekulární dynamiky * MeSH
• zvířata MeSH
• Check Tag
• krysa rodu Rattus MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• receptory N-methyl-D-aspartátu MeSH

A computational method of modeling random coils of hyaluronan was developed based on the molecular-dynamics simulations. An oligosaccharide of 48 monosaccharide units was equilibrated within a 70-100ns simulation and randomly chosen pieces of this molecule from different simulation frames were combined to constitute a long polysaccharide chain, both for hyaluronan and its non-ionic analog containing glucose instead of glucuronic acid. The dihedral angles of the glycoside connections of the pieces obeyed the statistics deduced from the simulation. The simulations were performed at various concentrations of NaCl and MgCl2. The calculated radii of gyration show a striking agreement with experimental data from the literature and indicate a key importance of the polymer-ion interactions for the random-coil conformation, but a low influence of the excluded volume of the chain and the carboxylate-groups repulsion. The method has thus the potential to become a versatile tool of modeling macromolecules of various semirigid polymers.

• Klíčová slova
• Dihedral angle, Electrolyte, Hyaluronan, Molecular dynamics, Radius of gyration,
• MeSH
• elektrolyty chemie   MeSH
• kyselina hyaluronová chemie   MeSH
• molekulární konformace MeSH
• simulace molekulární dynamiky * MeSH
• software MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• elektrolyty MeSH
• kyselina hyaluronová MeSH

Given by χ torsional angles, rotamers describe the side-chain conformations of amino acid residues in a protein based on the rotational isomers (hence the word rotamer). Constructed rotamer libraries, based on either protein crystal structures or dynamics studies, are the tools for classifying rotamers (torsional angles) in a way that reflect their frequency in nature. Rotamer libraries are routinely used in structure modeling and evaluation. In this perspective article, we would like to encourage researchers to apply rotamer analyses beyond their traditional use. Molecular dynamics (MD) of proteins highlight the in silico behavior of molecules in solution and thus can identify favorable side-chain conformations. In this article, we used simple computational tools to study rotamer dynamics (RD) in MD simulations. First, we isolated each frame in the MD trajectories in separate Protein Data Bank files via the cpptraj module in AMBER. Then, we extracted torsional angles via the Bio3D module in R language. The classification of torsional angles was also done in R according to the penultimate rotamer library. RD analysis is useful for various applications such as protein folding, study of rotamer-rotamer relationship in protein-protein interaction, real-time correlation between secondary structures and rotamers, study of flexibility of side chains in binding site for molecular docking preparations, use of RD as guide in functional analysis and study of structural changes caused by mutations, providing parameters for improving coarse-grained MD accuracy and speed, and many others. Major challenges facing RD to emerge as a new scientific field involve the validation of results via easy, inexpensive wet-lab methods. This realm is yet to be explored.

• MeSH
• isomerie MeSH
• konformace proteinů MeSH
• proteiny chemie   MeSH
• rotace * MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• přehledy MeSH
• Názvy látek
• proteiny MeSH

Riboswitches often occur in the 5'-untranslated regions of bacterial mRNA where they regulate gene expression. The preQ(1) riboswitch controls the biosynthesis of a hypermodified nucleoside queuosine in response to binding the queuosine metabolic intermediate. Structures of the ligand-bound and ligand-free states of the preQ(1) riboswitch from Thermoanaerobacter tengcongensis were determined recently by X-ray crystallography. We used multiple, microsecond-long molecular dynamics simulations (29 μs in total) to characterize the structural dynamics of preQ(1) riboswitches in both states. We observed different stabilities of the stem in the bound and free states, resulting in different accessibilities of the ribosome-binding site. These differences are related to different stacking interactions between nucleotides of the stem and the associated loop, which itself adopts different conformations in the bound and free states. We suggest that the loop not only serves to bind preQ(1) but also transmits information about ligand binding from the ligand-binding pocket to the stem, which has implications for mRNA accessibility to the ribosome. We explain functional results obscured by a high salt crystallization medium and help to refine regions of disordered electron density, which demonstrates the predictive power of our approach. Besides investigating the functional dynamics of the riboswitch, we have also utilized this unique small folded RNA system for analysis of performance of the RNA force field on the μs time scale. The latest AMBER parmbsc0χ(OL3) RNA force field is capable of providing stable trajectories of the folded molecule on the μs time scale. On the other hand, force fields that are not properly balanced lead to significant structural perturbations on the sub-μs time scale, which could easily lead to inappropriate interpretation of the simulation data.

• MeSH
• bakteriální RNA chemie   MeSH
• krystalografie rentgenová MeSH
• molekulární modely MeSH
• riboswitch * MeSH
• simulace molekulární dynamiky * MeSH
• Thermoanaerobacter chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• bakteriální RNA MeSH
• riboswitch * MeSH

Highly Ca2+ selective channels trigger a large variety of cellular signaling processes in both excitable and non-excitable cells. Among these channels, the Orai channel is unique in its activation mechanism and its structure. It mediates Ca2+ influx into the cytosol with an extremely small unitary conductance over longer time-scales, ranging from minutes up to several hours. Its activation is regulated by the Ca2+ content of the endoplasmic reticulum (ER). Depletion of luminal [Ca2+]ER is sensed by the STIM1 single transmembrane protein that directly binds and gates the Orai1 channel. Orai mediated Ca2+ influx increases cytosolic Ca2+ from 100 nM up to low micromolar range close to the pore and thereby forms Ca2+ microdomains. Hence, these features of the Orai channel can trigger long-term signaling processes without affecting the overall Ca2+ content of a single living cell. Here we focus on the architecture and dynamic conformational changes within the Orai channel. This review summarizes current achievements of molecular dynamics simulations in combination with live cell recordings to address gating and permeation of the Orai channel with molecular precision.

• Klíčová slova
• Calcium, Gating, MD-simulations, Molecular dynamics, Orai, Orai1, Pore, STIM,
• MeSH
• lidé MeSH
• nádorové proteiny chemie   metabolismus   MeSH
• protein ORAI1 chemie   metabolismus   MeSH
• protein STIM1 chemie   metabolismus   MeSH
• simulace molekulární dynamiky * MeSH
• vápník metabolismus   MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• přehledy MeSH
• Názvy látek
• nádorové proteiny MeSH
• ORAI1 protein, human MeSH Prohlížeč
• protein ORAI1 MeSH
• protein STIM1 MeSH
• STIM1 protein, human MeSH Prohlížeč
• vápník MeSH

To elucidate the consequences of the saturated-unsaturated nature of lipid surface films, monolayers formed by an equimolar mixture of 1-palmitoyl-2-oleyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipids are investigated in a wide range of surface pressures. As such mixtures share some features with naturally-occurring surfactants, for example the lung surfactant, the systems are studied at the temperature relevant for human body. All-atom molecular dynamics simulations and Langmuir trough experiments are employed. The binary lipid mixture is compared with the corresponding one-component systems. Atomistic-level alterations of monolayer molecular properties upon lateral compression are scrutinized. These involve elevation of lateral ordering of lipid chains, modulation of chain and headgroup orientation, and reduction of lipid hydration. The presence of the unsaturated POPC in the DPPC/POPC mixture reduces the liquid expanded-liquid condensed coexistence region and moderates the phase transition. Simulations predict that nanoscale lipid de-mixing occurs with small transient DPPC clusters emerging due to local fluctuations of the lateral lipid arrangement. A vertical sorting of lipids induced by lateral compression is also observed, with DPPC transferred toward the water phase. Both the conformational lipid alterations due to monolayer compression as well as the existence of lateral dynamic inhomogeneities of the lipid film are potentially pertain to dynamic and non-homogeneous lipid interfacial systems.

• Klíčová slova
• Langmuir trough, Lung surfactant, Molecular dynamics, Phospholipid monolayers,
• MeSH
• 1,2-dipalmitoylfosfatidylcholin analogy a deriváty   chemie   MeSH
• fosfatidylcholiny chemie   MeSH
• lipidy chemie   MeSH
• molekulární konformace MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• 1,2-dipalmitoyl-sn-glycero-3-ethylphosphocholine MeSH Prohlížeč
• 1,2-dipalmitoylfosfatidylcholin MeSH
• 1,2-oleoylphosphatidylcholine MeSH Prohlížeč
• fosfatidylcholiny MeSH
• lipidy MeSH

Human stimulator of interferon genes (hSTING) is a signaling adaptor protein that triggers innate immune system by response to cytosolic DNA and second messenger cyclic dinucleotides (CDNs). Natural CDNs contain purine nucleobase with different phosphodiester linkage types (3'-3', 2'-2' or mixed 2'-3'-linkages) and exhibit different binding affinity towards hSTING, ranging from micromolar to nanomolar. High-affinity CDNs are considered as suitable candidates for treatment of chronic hepatitis B and cancer. We have used molecular dynamics simulations to investigate dynamical aspects of binding of natural CDNs (specifically, 2'-2'-cGAMP, 2'-3'-cGAMP, 3'-3'-cGAMP, 3'-3'-c-di-AMP, and 3'-3'-c-di-GMP) with hSTINGwt protein. Our results revealed that CDN/hSTINGwt interactions are controlled by the balance between fluctuations (conformational changes) in the CDN ligand and the protein dynamics. Binding of different CDNs induces different degrees of conformational/dynamics changes in hSTINGwt ligand binding cavity, especially in α1-helices, the so-called lid region and α2-tails. The ligand residence time in hSTINGwt protein pocket depends on different contribution of R232 and R238 residues interacting with oxygen atoms of phosphodiester groups in ligand, water distribution around interacting charged centers (in protein residues and ligand) and structural stability of closed conformation state of hSTINGwt protein. These findings may perhaps guide design of new compounds modulating hSTING activity.Communicated by Ramaswamy H. Sarma.

• Klíčová slova
• Human stimulator of interferon genes STING, Molecular Dynamics, cyclic dinucleotides CDNs,
• MeSH
• dinukleosidfosfáty * chemie   MeSH
• DNA MeSH
• lidé MeSH
• ligandy MeSH
• oligonukleotidy MeSH
• simulace molekulární dynamiky * MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• dinukleosidfosfáty * MeSH
• DNA MeSH
• ligandy MeSH
• oligonukleotidy MeSH

Glutathione (GSH) is a common antioxidant and its biological activity depends on the conformation and protonation state. We used molecular dynamics, Raman and Raman optical activity (ROA) spectroscopies to investigate GSH structural changes in a broad pH range. Factor analysis of the spectra provided protonation constants (2.05, 3.45, 8.62, 9.41) in good agreement with previously published values. Following the analysis, spectra of differently protonated forms were obtained by extrapolation. The complete deprotonation of the thiol group above pH 11 was clearly visible in the spectra; however, many spectral features did not change much with pH. Experimental spectra at various pH values were decomposed into the simulated ones, which allowed us to study the conformer populations and quality of molecular dynamics (MD). According to this combined ROA/MD analysis conformation of the GSH backbone is affected by the pH changes only in a limited way. The combination of ROA with the computations thus has the potential to improve the MD force field and obtain more accurate populations of the conformer species. The methodology can be used for any molecule, but for a more detailed insight better computational techniques are needed in the future.

• Klíčová slova
• Raman optical activity, density functional theory, glutathione, molecular dynamics, spectra modeling,
• MeSH
• glutathion MeSH
• molekulární konformace MeSH
• optická otáčivost MeSH
• Ramanova spektroskopie * metody   MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• glutathion MeSH