protein dynamics modeling
Dotaz
Zobrazit nápovědu
The purpose of this quick guide is to help new modelers who have little or no background in comparative modeling yet are keen to produce high-resolution protein 3D structures for their study by following systematic good modeling practices, using affordable personal computers or online computational resources. Through the available experimental 3D-structure repositories, the modeler should be able to access and use the atomic coordinates for building homology models. We also aim to provide the modeler with a rationale behind making a simple list of atomic coordinates suitable for computational analysis abiding to principles of physics (e.g., molecular mechanics). Keeping that objective in mind, these quick tips cover the process of homology modeling and some postmodeling computations such as molecular docking and molecular dynamics (MD). A brief section was left for modeling nonprotein molecules, and a short case study of homology modeling is discussed.
- MeSH
- algoritmy MeSH
- aminokyseliny chemie MeSH
- biologické modely MeSH
- databáze proteinů MeSH
- internet MeSH
- ionty MeSH
- koncentrace vodíkových iontů MeSH
- ligandy MeSH
- počítačová simulace MeSH
- posttranslační úpravy proteinů MeSH
- proteiny chemie MeSH
- rozpouštědla MeSH
- sbalování proteinů MeSH
- simulace molekulární dynamiky MeSH
- simulace molekulového dockingu MeSH
- software MeSH
- strojové učení MeSH
- strukturní homologie proteinů MeSH
- voda MeSH
- výpočetní biologie metody MeSH
- zobrazování trojrozměrné metody MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem 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
Recently, two independent (15)N NMR relaxation studies indicated that in contrast to the decreased flexibility expected for induced-fit interactions, the backbone flexibility of major urinary protein isoform I (MUP-I) slightly increased upon complex formation with its natural pheromone 2-sec-butyl-4,5-dihydrothiazol. We have investigated the subtle details of molecular interactions by molecular dynamics simulations in explicit solvent. The calculated order parameters S(2) for a free- and ligand-bound protein supply evidence that mobility in various regions of MUP-I can be directly related to small conformational changes of the free- and complexed protein resulting from modifications of the hydrogen bonding network.
The serine/threonine protein phosphatase type 5 (PP5) is a promising target for designing new antitumor drugs. This enzyme is a member of the PPP phosphatases gene family, which catalyzes a dephosphorylation reaction: a regulatory process in the signal transduction pathway that controls various biological processes. The aim of this work is to study and compare the inhibition of PP5 by ten cantharidin-like inhibitors in order to bring about contributions relevant to the better comprehension of their inhibitory activity. In this theoretical investigation, we used molecular dynamics techniques to understand the role of key interactions that occur in the protein active site; QM calculations were employed to study the interaction mode of these inhibitors in the enzyme. In addition, atoms in molecules (AIM) calculations were carried out to characterize the chemical bonds among the atoms involved and investigate the orbital interactions with their respective energy values. The obtained results suggest that the Arg275, Asn303, His304, His352, Arg400, His427, Glu428, Val429, Tyr451, and Phe446 residues favorably contribute to the interactions between inhibitors and PP5. However, the Asp271 and Asp244 amino acid residues do not favor such interactions for some inhibitors. Through the QM calculations, we can suggest that the reactional energy of the coordination mechanism of these inhibitors in the PP5 active site is quite important and is responsible for the inhibitory activity. The AIM technique employed in this work was essential to get a better comprehension of the transition states acquired from the mechanism simulation. This work offers insights of how cantharidin-like inhibitors interact with human PP5, potentially allowing the design of more specific and even less cytotoxic drugs for cancer treatments. Graphical Abstract Interactions of cantharidin-like inhibitors with human protein phosphatase-5 in a Mg2+ system.
- MeSH
- hořčík chemie MeSH
- inhibitory enzymů farmakologie MeSH
- jaderné proteiny antagonisté a inhibitory chemie MeSH
- kantharidin analogy a deriváty chemie farmakologie MeSH
- katalytická doména MeSH
- kationty dvojmocné chemie MeSH
- lidé MeSH
- proteinfosfatasy antagonisté a inhibitory chemie MeSH
- racionální návrh léčiv MeSH
- simulace molekulární dynamiky * MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
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
Molecular motions of free and pheromone-bound mouse major urinary protein I, previously investigated by NMR relaxation, were simulated in 30 ns molecular dynamics (MD) runs. The backbone flexibility was described in terms of order parameters and correlation times, commonly used in the NMR relaxation analysis. Special attention was paid to the effect of conformational changes on the nanosecond time scale. Time-dependent order parameters were determined in order to separate motions occurring on different time scales. As an alternative approach, slow conformational changes were identified from the backbone torsion angle variances, and "conformationally filtered" order parameters were calculated for well-defined conformation states. A comparison of the data obtained for the free and pheromone-bound protein showed that some residues are more rigid in the bound form, but a larger portion of the protein becomes more flexible upon the pheromone binding. This finding is in general agreement with the NMR results. The higher flexibility observed on the fast (fs-ps) time scale was typically observed for the residues exhibiting higher conformational freedom on the ns time scale. An inspection of the hydrogen bond network provided a structural explanation for the flexibility differences between the free and pheromone-bound proteins in the simulations.
The 14-3-3 protein family performs regulatory functions in eukaryotic organisms by binding to a large number of phosphorylated protein partners. Whilst the binding mode of the phosphopeptides within the primary 14-3-3 binding site is well established based on the crystal structures of their complexes, little is known about the binding process itself. We present a computational study of the process by which phosphopeptides bind to the 14-3-3ζ protein. Applying a novel scheme combining Hamiltonian replica exchange molecular dynamics and distancefield restraints allowed us to map and compare the most likely phosphopeptide-binding pathways to the 14-3-3ζ protein. The most important structural changes to the protein and peptides involved in the binding process were identified. In order to bind phosphopeptides to the primary interaction site, the 14-3-3ζ adopted a newly found wide-opened conformation. Based on our findings we additionally propose a secondary interaction site on the inner surface of the 14-3-3ζ dimer, and a direct interference on the binding process by the flexible C-terminal tail. A minimalistic model was designed to allow for the efficient calculation of absolute binding affinities. Binding affinities calculated from the potential of mean force along the binding pathway are in line with the available experimental estimates for two of the studied systems.
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.
BACKGROUND: Protein function is determined by many factors, namely by its constitution, spatial arrangement, and dynamic behavior. Studying these factors helps the biochemists and biologists to better understand the protein behavior and to design proteins with modified properties. One of the most common approaches to these studies is to compare the protein structure with other molecules and to reveal similarities and differences in their polypeptide chains. RESULTS: We support the comparison process by proposing a new visualization technique that bridges the gap between traditionally used 1D and 3D representations. By introducing the information about mutual positions of protein chains into the 1D sequential representation the users are able to observe the spatial differences between the proteins without any occlusion commonly present in 3D view. Our representation is designed to serve namely for comparison of multiple proteins or a set of time steps of molecular dynamics simulation. CONCLUSIONS: The novel representation is demonstrated on two usage scenarios. The first scenario aims to compare a set of proteins from the family of cytochromes P450 where the position of the secondary structures has a significant impact on the substrate channeling. The second scenario focuses on the protein flexibility when by comparing a set of time steps our representation helps to reveal the most dynamically changing parts of the protein chain.
The effect of non-denaturing concentrations of three different organic solvents, formamide, acetone and isopropanol, on the structure of haloalkane dehalogenases DhaA, LinB, and DbjA at the protein-solvent interface was studied using molecular dynamics simulations. Analysis of B-factors revealed that the presence of a given organic solvent mainly affects the dynamical behavior of the specificity-determining cap domain, with the exception of DbjA in acetone. Orientation of organic solvent molecules on the protein surface during the simulations was clearly dependent on their interaction with hydrophobic or hydrophilic surface patches, and the simulations suggest that the behavior of studied organic solvents in the vicinity of hyrophobic patches on the surface is similar to the air/water interface. DbjA was the only dimeric enzyme among studied haloalkane dehalogenases and provided an opportunity to explore effects of organic solvents on the quaternary structure. Penetration and trapping of organic solvents in the network of interactions between both monomers depends on the physico-chemical properties of the organic solvents. Consequently, both monomers of this enzyme oscillate differently in different organic solvents. With the exception of LinB in acetone, the structures of studied enzymes were stabilized in water-miscible organic solvents.
- MeSH
- 2-propanol chemie farmakologie MeSH
- aceton chemie farmakologie MeSH
- formamidy chemie farmakologie MeSH
- hydrofobní a hydrofilní interakce MeSH
- hydrolasy chemie MeSH
- krystalografie rentgenová MeSH
- kvarterní struktura proteinů účinky léků MeSH
- molekulární modely MeSH
- rozpouštědla chemie MeSH
- simulace molekulární dynamiky MeSH
- terciární struktura proteinů účinky léků MeSH
- voda chemie MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
