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.
- MeSH
- Dinucleoside Phosphates * chemistry MeSH
- DNA MeSH
- Humans MeSH
- Ligands MeSH
- Oligonucleotides MeSH
- Molecular Dynamics Simulation * MeSH
- Check Tag
- Humans MeSH
- Publication type
- Journal Article 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 family of stromal interaction molecules (STIM) includes two widely expressed single-pass endoplasmic reticulum (ER) transmembrane proteins and additional splice variants that act as precise ER-luminal Ca2+ sensors. STIM proteins mainly function as one of the two essential components of the so-called Ca2+ release-activated Ca2+ (CRAC) channel. The second CRAC channel component is constituted by pore-forming Orai proteins in the plasma membrane. STIM and Orai physically interact with each other to enable CRAC channel opening, which is a critical prerequisite for various downstream signalling pathways such as gene transcription or proliferation. Their activation commonly requires the emptying of the intracellular ER Ca2+ store. Using their Ca2+ sensing capabilities, STIM proteins confer this Ca2+ content-dependent signal to Orai, thereby linking Ca2+ store depletion to CRAC channel opening. Here we review the conformational dynamics occurring along the entire STIM protein upon store depletion, involving the transition from the quiescent, compactly folded structure into an active, extended state, modulation by a variety of accessory components in the cell as well as the impairment of individual steps of the STIM activation cascade associated with disease.
- MeSH
- Calcium Release Activated Calcium Channels * MeSH
- Membrane Proteins metabolism MeSH
- ORAI1 Protein MeSH
- Stromal Interaction Molecule 1 metabolism MeSH
- Stromal Interaction Molecules * metabolism MeSH
- Calcium metabolism MeSH
- Calcium Signaling physiology MeSH
- Publication type
- Journal Article MeSH
- Review MeSH
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.
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.
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.
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
- Nucleic Acid Conformation MeSH
- Humans MeSH
- RNA-Binding Proteins chemistry metabolism MeSH
- RNA chemistry metabolism MeSH
- Molecular Dynamics Simulation * MeSH
- Computational Biology methods MeSH
- Animals MeSH
- Check Tag
- Humans MeSH
- Animals MeSH
- Publication type
- Journal Article MeSH
- Research Support, Non-U.S. Gov't MeSH
- Review 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.
- MeSH
- Cytochrome P-450 CYP1A2 chemistry MeSH
- Cytochromes b5 chemistry MeSH
- Protein Conformation MeSH
- Humans MeSH
- Molecular Dynamics Simulation MeSH
- Molecular Docking Simulation MeSH
- Static Electricity MeSH
- Check Tag
- Humans MeSH
- Publication type
- Journal Article MeSH
- Research Support, Non-U.S. Gov't 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
- Adenosine Triphosphate chemistry metabolism MeSH
- Ion Transport MeSH
- Protein Conformation * MeSH
- Fatty Acids chemistry metabolism MeSH
- Membrane Proteins chemistry MeSH
- Mitochondrial Proteins chemistry metabolism MeSH
- Mice MeSH
- Amino Acid Sequence MeSH
- Molecular Dynamics Simulation * MeSH
- Protein Stability MeSH
- Uncoupling Protein 2 chemistry metabolism MeSH
- Protein Binding MeSH
- Structure-Activity Relationship MeSH
- Animals MeSH
- Check Tag
- Mice MeSH
- Animals MeSH
- Publication type
- Journal Article MeSH
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 chemistry pharmacology MeSH
- Acetone chemistry pharmacology MeSH
- Formamides chemistry pharmacology MeSH
- Hydrophobic and Hydrophilic Interactions MeSH
- Hydrolases chemistry MeSH
- Crystallography, X-Ray MeSH
- Protein Structure, Quaternary drug effects MeSH
- Models, Molecular MeSH
- Solvents chemistry MeSH
- Molecular Dynamics Simulation MeSH
- Protein Structure, Tertiary drug effects MeSH
- Water chemistry MeSH
- Publication type
- Journal Article MeSH
- Research Support, Non-U.S. Gov't MeSH
