Protein dynamics Dotaz Zobrazit nápovědu
Cold Spring Harbor symposia on quantitative biology ; Vol. 60
[1st ed.] xxiv, 843 s. : il.
- Konspekt
- Biochemie. Molekulární biologie. Biofyzika
- NLK Obory
- biochemie
- biologie
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
Protein dynamics are often invoked in explanations of enzyme catalysis, but their design has proven elusive. Here we track the role of dynamics in evolution, starting from the evolvable and thermostable ancestral protein AncHLD-RLuc which catalyses both dehalogenase and luciferase reactions. Insertion-deletion (InDel) backbone mutagenesis of AncHLD-RLuc challenged the scaffold dynamics. Screening for both activities reveals InDel mutations localized in three distinct regions that lead to altered protein dynamics (based on crystallographic B-factors, hydrogen exchange, and molecular dynamics simulations). An anisotropic network model highlights the importance of the conformational flexibility of a loop-helix fragment of Renilla luciferases for ligand binding. Transplantation of this dynamic fragment leads to lower product inhibition and highly stable glow-type bioluminescence. The success of our approach suggests that a strategy comprising (i) constructing a stable and evolvable template, (ii) mapping functional regions by backbone mutagenesis, and (iii) transplantation of dynamic features, can lead to functionally innovative proteins.
- MeSH
- buňky NIH 3T3 MeSH
- katalýza MeSH
- kinetika MeSH
- konformace proteinů MeSH
- luciferasy renil chemie genetika metabolismus MeSH
- luciferasy chemie genetika metabolismus MeSH
- mutace MeSH
- mutageneze MeSH
- myši MeSH
- proteinové inženýrství * MeSH
- savci MeSH
- simulace molekulární dynamiky * MeSH
- stabilita enzymů MeSH
- teplota MeSH
- vazebná místa MeSH
- zvířata MeSH
- Check Tag
- myši MeSH
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem 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.
Motivation: Studying the transport paths of ligands, solvents, or ions in transmembrane proteins and proteins with buried binding sites is fundamental to the understanding of their biological function. A detailed analysis of the structural features influencing the transport paths is also important for engineering proteins for biomedical and biotechnological applications. Results: CAVER Analyst 2.0 is a software tool for quantitative analysis and real-time visualization of tunnels and channels in static and dynamic structures. This version provides the users with many new functions, including advanced techniques for intuitive visual inspection of the spatiotemporal behavior of tunnels and channels. Novel integrated algorithms allow an efficient analysis and data reduction in large protein structures and molecular dynamic simulations. Availability and implementation: CAVER Analyst 2.0 is a multi-platform standalone Java-based application. Binaries and documentation are freely available at www.caver.cz. Supplementary information: Supplementary data are available at Bioinformatics online.
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
Formation of transient complexes of cytochrome P450 (P450) with another protein of the endoplasmic reticulum membrane, cytochrome b5 (cyt b5), dictates the catalytic activities of several P450s. Therefore, we examined formation and binding modes of the complex of human P450 1A2 with cyt b5. Docking of soluble domains of these proteins was performed using an information-driven flexible docking approach implemented in HADDOCK. Stabilities of the five unique binding modes of the P450 1A2-cyt b5 complex yielded by HADDOCK were evaluated using explicit 10 ns molecular dynamics (MD) simulations in aqueous solution. Further, steered MD was used to compare the stability of the individual P450 1A2-cyt b5 binding modes. The best binding mode was characterized by a T-shaped mutual orientation of the porphyrin rings and a 10.7 Å distance between the two redox centers, thus satisfying the condition for a fast electron transfer. Mutagenesis studies and chemical cross-linking, which, in the absence of crystal structures, were previously used to deduce specific P450-cyt b5 interactions, indicated that the negatively charged convex surface of cyt b5 binds to the positively charged concave surface of P450. Our simulations further elaborate structural details of this interface, including nine ion pairs between R95, R100, R138, R362, K442, K455, and K465 side chains of P450 1A2 and E42, E43, E49, D65, D71, and heme propionates of cyt b5. The universal heme-centric system of internal coordinates was proposed to facilitate consistent classification of the orientation of the two porphyrins in any protein complex.
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.