Molecular Dynamics
Dotaz
Zobrazit nápovědu
... Numerical Methods for Stochastic Molecular Dynamics -- 8. ...
Interdisciplinary Applied Mathematics, ISSN 0939-6047 39
1st edition XXII, 443 stran : ilustrace ; 24 cm
- MeSH
- matematika MeSH
- simulace molekulární dynamiky MeSH
- Publikační typ
- monografie MeSH
- Konspekt
- Matematika
- NLK Obory
- přírodní vědy
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.
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.
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.
In this study, quantitative structure-activity relationship studies which make use of molecular dynamics trajectories were performed on a set of 54 glucokinase protein activators. The conformations obtained by molecular dynamics simulation were superimposed according to the twelve alignments tested in a virtual three-dimensional box comprised of 2 Å cells. The models were generated by the technique that combines genetic algorithms and partial least squares. The best alignment models generated with a determination coefficient (r(2)) between 0.674 and 0.743 and cross-validation (q(2)) between 0.509 and 0.610, indicating good predictive capacity. The 4D-QSAR models developed in this study suggest novel molecular regions to be explored in the search for better glucokinase activators.
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.
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.
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.
Ricin is a toxin found in the castor seeds and listed as a chemical weapon by the Chemical Weapons Convention (CWC) due to its high toxicity combined with the easiness of obtention and lack of available antidotes. The relatively frequent episodes of usage or attempting to use ricin in terrorist attacks reinforce the urge to develop an antidote for this toxin. In this sense, we selected in this work the current RTA (ricin catalytic subunit) inhibitor with the best experimental performance, as a reference molecule for virtual screening in the PubChem database. The selected molecules were then evaluated through docking studies, followed by drug-likeness investigation, molecular dynamics simulations and Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) calculations. In every step, the selection of molecules was mainly based on their ability to occupy both the active and secondary sites of RTA, which are located right next to each other, but are not simultaneously occupied by the current RTA inhibitors. Results show that the three PubChem compounds 18309602, 18498053, and 136023163 presented better overall results than the reference molecule itself, showing up as new hits for the RTA inhibition, and encouraging further experimental evaluation.
- MeSH
- algoritmy MeSH
- chemické bojové látky chemie MeSH
- ligandy MeSH
- molekulární konformace MeSH
- molekulární struktura MeSH
- objevování léků MeSH
- ricin antagonisté a inhibitory chemie MeSH
- simulace molekulární dynamiky MeSH
- simulace molekulového dockingu MeSH
- vazebná místa MeSH
- vodíková vazba MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Molecular dynamics (MD) simulations of uncoupling proteins (UCP), a class of transmembrane proteins relevant for proton transport across inner mitochondrial membranes, represent a complicated task due to the lack of available structural data. In this work, we use a combination of homology modelling and subsequent microsecond molecular dynamics simulations of UCP2 in the DOPC phospholipid bilayer, starting from the structure of the mitochondrial ATP/ADP carrier (ANT) as a template. We show that this protocol leads to a structure that is impermeable to water, in contrast to MD simulations of UCP2 structures based on the experimental NMR structure. We also show that ATP binding in the UCP2 cavity is tight in the homology modelled structure of UCP2 in agreement with experimental observations. Finally, we corroborate our results with conductance measurements in model membranes, which further suggest that the UCP2 structure modeled from ANT protein possesses additional key functional elements, such as a fatty acid-binding site at the R60 region of the protein, directly related to the proton transport mechanism across inner mitochondrial membranes.
- MeSH
- adenosintrifosfát chemie metabolismus MeSH
- iontový transport MeSH
- konformace proteinů * MeSH
- mastné kyseliny chemie metabolismus MeSH
- membránové proteiny chemie MeSH
- mitochondriální proteiny chemie metabolismus MeSH
- myši MeSH
- sekvence aminokyselin MeSH
- simulace molekulární dynamiky * MeSH
- stabilita proteinů MeSH
- uncoupling protein 2 chemie metabolismus MeSH
- vazba proteinů MeSH
- vztahy mezi strukturou a aktivitou MeSH
- zvířata MeSH
- Check Tag
- myši MeSH
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
