molecular simulations
Dotaz
Zobrazit nápovědu
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.
Fluorescence solvent relaxation experiments are based on the characterization of time-dependent shifts in the fluorescence emission of a chromophore, yielding polarity and viscosity information about the chromophore's immediate environment. A chromophore applied to a phospholipid bilayer at a well-defined location (with respect to the z-axis of the bilayer) allows monitoring of the hydration and mobility of the probed segment of the lipid molecules. Specifically, time-resolved fluorescence experiments, fluorescence quenching data and molecular dynamic (MD) simulations show that 6-lauroyl-2-dimethylaminonaphthalene (Laurdan) probes the hydration and mobility of the sn-1 acyl groups in a phosphatidylcholine bilayer. The time-dependent fluorescence shift (TDFS) of Laurdan provides information on headgroup compression and expansion induced by the addition of different amounts of cationic lipids to phosphatidylcholine bilayers. Those changes were predicted by previous MD simulations. Addition of truncated oxidized phospholipids leads to increased mobility and hydration at the sn-1 acyl level. This experimental finding can be explained by MD simulations, which indicate that the truncated chains of the oxidized lipid molecules are looping back into aqueous phase, hence creating voids below the glycerol level. Fluorescence solvent relaxation experiments are also useful in understanding salt effects on the structure and dynamics of lipid bilayers. For example, such experiments demonstrate that large anions increase hydration and mobility at the sn-1 acyl level of phosphatidylcholine bilayers, an observation which could not be explained by standard MD simulations. If polarizability is introduced into the applied force field, however, MD simulations show that big soft polarizable anions are able to interact with the hydrophilic/hydrophobic interface of the lipid bilayer, penetrating to the level probed by Laurdan, and that they expand and destabilize the bilayer making it more hydrated and mobile.
We provide a critical assessment of explicit-solvent atomistic molecular dynamics (MD) simulations of RNA and protein/RNA complexes, written primarily for non-specialists with an emphasis to explain the limitations of MD. MD simulations can be likened to hypothetical single-molecule experiments starting from single atomistic conformations and investigating genuine thermal sampling of the biomolecules. The main advantage of MD is the unlimited temporal and spatial resolution of positions of all atoms in the simulated systems. Fundamental limitations are the short physical time-scale of simulations, which can be partially alleviated by enhanced-sampling techniques, and the highly approximate atomistic force fields describing the simulated molecules. The applicability and present limitations of MD are demonstrated on studies of tetranucleotides, tetraloops, ribozymes, riboswitches and protein/RNA complexes. Wisely applied simulations respecting the approximations of the model can successfully complement structural and biochemical experiments. WIREs RNA 2017, 8:e1405. doi: 10.1002/wrna.1405 For further resources related to this article, please visit the WIREs website.
- MeSH
- konformace nukleové kyseliny MeSH
- lidé MeSH
- proteiny vázající RNA chemie metabolismus MeSH
- RNA chemie metabolismus MeSH
- simulace molekulární dynamiky * MeSH
- výpočetní biologie metody MeSH
- zvířata MeSH
- Check Tag
- lidé MeSH
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- přehledy MeSH
Explicit solvent molecular dynamics simulations have been used to complement preceding experimental and computational studies of folding of guanine quadruplexes (G-DNA). We initiate early stages of unfolding of several G-DNAs by simulating them under no-salt conditions and then try to fold them back using standard excess salt simulations. There is a significant difference between G-DNAs with all-anti parallel stranded stems and those with stems containing mixtures of syn and anti guanosines. The most natural rearrangement for all-anti stems is a vertical mutual slippage of the strands. This leads to stems with reduced numbers of tetrads during unfolding and a reduction of strand slippage during refolding. The presence of syn nucleotides prevents mutual strand slippage; therefore, the antiparallel and hybrid quadruplexes initiate unfolding via separation of the individual strands. The simulations confirm the capability of G-DNA molecules to adopt numerous stable locally and globally misfolded structures. The key point for a proper individual folding attempt appears to be correct prior distribution of syn and anti nucleotides in all four G-strands. The results suggest that at the level of individual molecules, G-DNA folding is an extremely multi-pathway process that is slowed by numerous misfolding arrangements stabilized on highly variable timescales.
- MeSH
- DNA chemie MeSH
- G-kvadruplexy * MeSH
- jednovláknová DNA chemie MeSH
- lidé MeSH
- simulace molekulární dynamiky * MeSH
- telomery chemie MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Biomolecular simulations are routinely used in biochemistry and molecular biology research; however, they often fail to match expectations of their impact on pharmaceutical and biotech industry. This is caused by the fact that a vast amount of computer time is required to simulate short episodes from the life of biomolecules. Several approaches have been developed to overcome this obstacle, including application of massively parallel and special purpose computers or non-conventional hardware. Methodological approaches are represented by coarse-grained models and enhanced sampling techniques. These techniques can show how the studied system behaves in long time-scales on the basis of relatively short simulations. This review presents an overview of new simulation approaches, the theory behind enhanced sampling methods and success stories of their applications with a direct impact on biotechnology or drug design.
Molecular modeling in combination with powder X-ray diffraction (XRD) provided new information on the organization of the interlayer space of Mg-Al layered double hydroxide (LDH) containing intercalated porphyrin anions [5,10,15,20-tetrakis(4-sulfonatophenyl)porphyrin (TPPS)]. Anion-exchange and rehydration procedures were used for the preparation of TPPS-containing LDH with an Mg/Al molar ratio of 2. Molecular modeling was carried out in the Cerius(2) and Materials Studio modeling environment. Three types of models were created in order to simulate the experimental XRD patterns of LDH intercalates with a TPPS loading of 70-80% with respect to the theoretical anion exchange capacity (AEC). The models represent single-phase systems with a 100% TPPS loading in the interlayer space (Type 1) and models represent the coexistence of two phases corresponding to the total exchange from 75 to 92% (Type 2). To cover other possible arrangements, models with the coexistence of both TPPS and NO(3)(-) anions in the same interlayer space were calculated (Type 3). The models are described and compared with experimental data. In all cases, guest TPPS anions are tilted with respect to the hydroxide layers, and are horizontally shifted to each other by up to one-half of the TPPS diameter. According to the energy characteristics and simulated XRD, the most probable arrangement is of Type 2, where some layers are saturated with TPPS anions and others are filled with original NO(3)(-) anions.
We have carried out extended set of μs-scale explicit solvent MD simulations of all possible G-triplexes which can participate in folding pathways of the human telomeric quadruplex. Our study accumulates almost 60 μs of simulation data, which is by about three orders of magnitude larger sampling compared to the earlier simulations of human telomeric G-DNA triplexes. Starting structures were obtained from experimental quadruplex structures by deleting either the first or the last strand. The life-times of antiparallel triplexes with lateral and diagonal loops are at least on μs-scale, which should be sufficient to contribute to the folding pathways. However, the triplex states may involve structures with various local deviations from the ideal triplexes, such as strand tilting and various alternative and incomplete triads. The simulations reveal easy rearrangements between lateral and diagonal loop triplex topologies. Propeller loops of antiparallel triplexes may to certain extent interfere with the G-triplexes but these structures are still viable candidates to participate in the folding. In contrast, all-parallel all-anti triplexes are very unstable and are unlikely to contribute to the folding. Although our simulations demonstrate that antiparallel G-triplexes, if folded, would have life-times sufficient to participate in the quadruplex folding, the results do not rule out the possibility that the G-triplexes are out-competed by other structures not included in our study. Among them, numerous possible misfolded structures containing guanine quartets can act as off-path intermediates with longer life-times than the triplexes. Besides analyzing the structural dynamics of a diverse set of G-DNA triplexes, we also provide a brief discussion of the limitations of the simulation methodology, which is necessary for proper understanding of the simulation data.
- MeSH
- DNA chemie MeSH
- G-kvadruplexy * MeSH
- konformace nukleové kyseliny * MeSH
- lidé MeSH
- nukleové kyseliny chemie genetika MeSH
- sekvence nukleotidů MeSH
- simulace molekulární dynamiky MeSH
- telomery chemie genetika MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
BACKGROUND: Innovations in computer hardware and software capabilities have paved the way for advances in molecular modelling techniques and methods, leading to an unprecedented expansion of their potential applications. In contrast to the docking technique, which usually identifies the most stable selector-selectand (SO-SA) complex for each enantiomer, the molecular dynamics (MD) technique enables the consideration of a distribution of the SO-SA complexes based on their energy profile. This approach provides a more truthful representation of the processes occurring within the column. However, benchmark procedures and focused guidelines for computational treatment of enantioselectivity at the molecular level are still missing. RESULTS: Twenty-eight molecular dynamics simulations were performed to study the enantiorecognition mechanisms of seven N-3,5-dinitrobenzoylated α- and β-amino acids (DNB-AAs), occurring with the two quinine- and quinidine-based (QN-AX and QD-AX) chiral stationary phases (CSPs), under polar-ionic conditions. The MD protocol was optimized in terms of box size, simulation run time, and frame recording frequency. Subsequently, all the trajectories were analyzed by calculating both the type and amount of the interactions engaged by the selectands (SAs) with the two chiral selectors (SOs), as well as the conformational and interaction energy profiles of the formed SA-SO associates. All the MDs were in strict agreement with the experimental enantiomeric elution order and allowed to establish (i) that salt-bridge and H-bond interactions play a pivotal role in the enantiorecognition mechanisms, and (ii) that the π-cation and π-π interactions are the discriminant chemical features between the two SOs in ruling the chiral recognition mechanism. SIGNIFICANCE: The results of this work clearly demonstrate the high contribution given by MD simulations in the comprehension of the enantiorecognition mechanism with Cinchona alkaloid-based CSPs. However, from this research endeavor it clearly emerged that the MD protocol optimization is crucial for the quality of the produced results.
The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.
- MeSH
- DNA chemie MeSH
- G-kvadruplexy * MeSH
- guanin chemie MeSH
- konformace nukleové kyseliny MeSH
- ligandy MeSH
- RNA chemie MeSH
- rozpouštědla chemie MeSH
- simulace molekulární dynamiky * MeSH
- telomery chemie MeSH
- vodíková vazba MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- přehledy MeSH
- Research Support, N.I.H., Extramural MeSH
- Research Support, U.S. Gov't, Non-P.H.S. MeSH
The flexibility, active site volume, solvation, and access path dynamics of six metabolically active mammalian cytochromes P450 (human 2A6, 2C9, 2D6, 2E1, 3A4 and rabbit 2B4) are extensively studied using molecular dynamics (MD) simulations. On average, the enzymes' overall structures equilibrate on a 50+ ns timescale. The very open CYP2B4 structure closes slowly over the course of the simulation. The volumes of the active sites fluctuate by more than 50% during the MD runs; these fluctuations are mainly due to movements of the main chains, with only a handful of amino acid residues in CYP2B4, CYP2D6, CYP2A6 and CYP2C9 showing significant independent side chain movement. The volume of the active site of CYP2E1 fluctuates heavily, ranging from 220 to 1310 A(3), due to the opening and closing of gates to two adjacent cavities. CYP2E1 has the least hydrated active site of the studied CYPs; this is consistent with its preference for non-polar substrates. The CYP2A6 and CYP2E1 active sites are deeply buried, with access paths that are narrower than the radius of a water molecule. However, waters are still able to access these active sites due to local adaptations of the channel to accommodate their passage. This finding may imply that the access paths of the CYPs never fully open prior to contact with the substrate; instead, the substrate may induce adaptive conformational changes during its passage to the active site. This may also explain why some substrate recognition sites are localized along individual enzymes' access paths.
- MeSH
- katalytická doména MeSH
- králíci MeSH
- lidé MeSH
- simulace molekulární dynamiky MeSH
- systém (enzymů) cytochromů P-450 chemie MeSH
- zvířata MeSH
- Check Tag
- králíci MeSH
- lidé MeSH
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
