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

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

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

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

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

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

• MeSH
• DNA chemie   MeSH
• katalýza MeSH
• konformace nukleové kyseliny * MeSH
• počítačová simulace MeSH
• RNA chemie   MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• přehledy MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• DNA MeSH
• RNA 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

Abstract Aromatic stacking of nucleic acid bases is one of the key players in determining the structure and dynamics of nucleic acids. The arrangement of nucleic acid bases with extensive overlap of their aromatic rings gave rise to numerous often contradictory suggestions about the physical origins of stacking and the possible role of delocalized electrons in stacked aromatic π systems, leading to some confusion about the issue. The recent advance of computer hardware and software finally allowed the application of state of the art quantum-mechanical approaches with inclusion of electron correlation effects to study aromatic base stacking, now providing an ultimitate qualitative description of the phenomenon. Base stacking is determined by an interplay of the three most commonly encountered molecular interactions: dispersion attraction, electrostatic interaction, and short-range repulsion. Unusual (aromatic- stacking specific) energy contributions were in fact not evidenced and are not necessary to describe stacking. The currently used simple empirical potential form, relying on atom-centered constant point charges and Lennard-Jones van der Waals terms, is entirely able to reproduce the essential features of base stacking. Thus, we can conclude that base stacking is in principle one of the best described interactions in current molecular modeling and it allows to study base stacking in DNA using large-scale classical molecular dynamics simulations. Neglect of cooperativity of stacking appears to be the most serious approximation of the currently used force field form. This review summarizes recent developments in the field. It is written for an audience that is not necessarily expert in computational quantum chemistry and follows up on our previous contribution (Sponer et. al., J. Biomol. Struct. Dyn. 14, 117, (1997)). First, the applied methodology, its accuracy, and the physical nature of base stacking is briefly overviewed, including a comment on the accuracy of other molecular orbital methods and force fields. Then, base stacking is contrasted with hydrogen bonding, the other dominant force in nucleic acid structure. The sequence dependence and cooperativity of base stacking is commented on, and finally a brief introduction into recent progress in large-scale molecular dynamics simulations of nucleic acids is provided. Using four stranded DNA assemblies as an example, we demonstrate the efficacy of current molecular dynamics techniques that utilize refined and verified force fields in the study of stacking in nucleic acid molecules.

• MeSH
• DNA * chemie   MeSH
• konformace nukleové kyseliny MeSH
• molekulární modely MeSH
• nukleové kyseliny chemie   MeSH
• simulace molekulární dynamiky * MeSH
• statická elektřina MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• přehledy MeSH
• Názvy látek
• DNA * MeSH
• nukleové kyseliny MeSH

A recent study described an allosteric effect in which the binding of a protein to DNA is influenced by another protein bound nearby. The effect shows a periodicity of ∼10 basepairs and decays with increasing protein-protein distance. As a mechanistic explanation, the authors reported a similar periodic, decaying pattern of the correlation coefficient between major groove widths inferred from a shorter molecular dynamics simulation. Here we show that in a state-of-the-art, microsecond-long simulation of the same DNA sequence, the periodicity of the correlation coefficient is not observed. To study the problem further, we extend an earlier mechanical model of DNA allostery based on constrained minimization of effective quadratic deformation energy of the DNA. We demonstrate that, if the constraints mimicking the bound proteins are properly applied, the periodicity in the binding energy is indeed recovered.

• MeSH
• alosterická regulace MeSH
• DNA chemie   MeSH
• konformace nukleové kyseliny * MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA MeSH

Although there were experimental indications that phospholipid bilayers hydrated with D(2)O express different biophysical properties compared with hydration by ordinary H(2)O, a molecular concept for this behavior difference was only recently proposed by a molecular dynamics simulations study [T. Róg et al., J. Phys. Chem. B, 2009, 113, 2378-2387]. Here we attempt to verify those theoretical predictions by fluorescence measurements on 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membranes. Specifically, we determine the water isotope effect on headgroup hydration and mobility, lateral lipid diffusion and lipid backbone packing. Time-dependent fluorescence shift experiments show significantly slower dynamics and lower hydration of the headgroup region for a bilayer hydrated with D(2)O, an observation in good agreement with the calculated predicted differences in duration of lipid-lipid and lipid-water bridges and extent of water penetration into the bilayer, respectively. The water isotope effect on the lipid order parameter of the bilayer core (measured by fluorescence anisotropy) and lateral diffusion of lipid molecules (determined by two-focus fluorescence correlation spectroscopy) is close to the experimental errors of the experiments, however also refers to slightly more rigid organization of phospholipid bilayers in heavy water. This study confirms the view that the water isotope effect can be particularly found in time-resolved physicochemical properties of the membrane. Together with the simulations our experiments provide a comprehensive, molecular view on the effect of D(2)O on phospholipid bilayers.

• MeSH
• časové faktory MeSH
• fluorescence MeSH
• fosfolipidy chemie   MeSH
• lipidové dvojvrstvy chemie   MeSH
• molekulární konformace MeSH
• oxid deuteria chemie   MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• fosfolipidy MeSH
• lipidové dvojvrstvy MeSH
• oxid deuteria MeSH