Molecular Dynamics Simulations of HLA-CW4-Β2M-KIR2DL1 Protein and Homology Modeling of a Complex Associated with Psoriasis Disease (HLA-CW6-Β 2M-KIR2DS1)

2021 ◽  
Author(s):  
Mansour H Almatarneh ◽  
Ahmad M Alqaisi ◽  
Enas K Ibrahim ◽  
Ghada G Kayed ◽  
Joshua W Hollett
Keyword(s):  
Molecular Dynamics ◽  
Conformational Changes ◽  
Complex Structure ◽  
3D Structure ◽  
Md Simulations ◽  
Salt Bridges ◽  
Small Peptide ◽  
Β2 Microglobulin ◽  
Dynamics Simulations ◽  
The Right

Molecular dynamics (MD) simulation was used to study the interactions of two immune proteins of HLA-Cw4-β2m-KIR2DL1 complex with small peptide QYDDAVYKL (nine amino acids) in an aqueous solution. This study aims to gain a detailed information about the conformational changes and the dynamics of the complex. The right parameters and force field for performing the MD simulations that was needed to calibrate the complex structure were determined. The non-bonded interactions (Electrostatic and van der Waals contributions), H-bond formation, and salt bridges between the ligand HLA-Cw4 and the receptor KIR2DL1 were estimated using the obtained MD trajectories. The buried surface area due to binding was calculated to get insight into the causes of specificity of receptor to ligand and explains mutations experiment. The study concluded that β2-microglobulin, one part of the complex, is not directly interacting with the peptide at the groove; therefore, it could be neglected from simulation. Our results showed that β2-microglobulin does not have any significant effect on the dynamics of the 3D-structure of the complex. This project will help in understanding to optimize candidate drug design, a small peptide that disrupts the interaction, for the optimal biological effect.

scholarly journals The β-cyclodextrin/benzene complex and its hydrogen bonds – a theoretical study using molecular dynamics, quantum mechanics and COSMO-RS

2013 ◽  
Vol 9 ◽  
pp. 118-134 ◽  
Author(s):  
Jutta Erika Helga Köhler ◽  
Nicole Grczelschak-Mick
Keyword(s):  
Molecular Dynamics ◽  
Quantum Mechanics ◽  
Hydrogen Bonds ◽  
Md Simulations ◽  
Benzene Molecule ◽  
Glucose Unit ◽  
Hand Orientation ◽  
Different Temperatures ◽  
Dynamics Simulations ◽  
The Right

Four highly ordered hydrogen-bonded models of β-cyclodextrin (β-CD) and its inclusion complex with benzene were investigated by three different theoretical methods: classical quantum mechanics (QM) on AM1 and on the BP/TZVP-DISP3 level of approximation, and thirdly by classical molecular dynamics simulations (MD) at different temperatures (120 K and 273 to 300 K). The hydrogen bonds at the larger O2/O3 rim of empty β-CDs prefer the right-hand orientation, e.g., O3-H…O2-H in the same glucose unit and bifurcated towards …O4 and O3 of the next glucose unit on the right side. On AM1 level the complex energy was −2.75 kcal mol−1 when the benzene molecule was located parallel inside the β-CD cavity and −2.46 kcal mol−1 when it was positioned vertically. The AM1 HOMO/LUMO gap of the empty β-CD with about 12 eV is lowered to about 10 eV in the complex, in agreement with data from the literature. AM1 IR spectra displayed a splitting of the O–H frequencies of cyclodextrin upon complex formation. At the BP/TZVP-DISP3 level the parallel and vertical positions from the starting structures converged to a structure where benzene assumes a more oblique position (−20.16 kcal mol−1 and −20.22 kcal mol−1, resp.) as was reported in the literature. The character of the COSMO-RS σ-surface of β-CD was much more hydrophobic on its O6 rim than on its O2/O3 side when all hydrogen bonds were arranged in a concerted mode. This static QM picture of the β-CD/benzene complex at 0 K was extended by MD simulations. At 120 K benzene was mobile but always stayed inside the cavity of β-CD. The trajectories at 273, 280, 290 and 300 K certainly no longer displayed the highly ordered hydrogen bonds of β-CD and benzene occupied many different positions inside the cavity, before it left the β-CD finally at its O2/O3 side.

Molecular Basis Explanation of Imatinib Resistance of Bcr-Abl Due to T315I and P-Loop Mutations from Molecular Dynamics Simulations.

Blood ◽  
2007 ◽  
Vol 110 (11) ◽  
pp. 2917-2917
Author(s):  
Tai-Sung Lee ◽  
Steven Potts ◽  
Hagop Kantarjian ◽  
Jorge Cortes ◽  
Francis Giles ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Conformational Changes ◽  
Electrostatic Interactions ◽  
Large Scale ◽  
Resistance Mechanism ◽  
Md Simulations ◽  
Resistance Mechanisms ◽  
Static Structure ◽  
Imatinib Resistance ◽  
Dynamics Simulations

Abstract Molecular dynamics (MD) simulations on the complex of imatinib with the wild-type, T315I, and other 10 P-loop mutants of the tyrosine kinase Bcr-Abl have been performed to study the imatinib resistance mechanism at the atomic level. MD simulations show that large scale computational simulations could offer insight information that a static structure or simple homology modeling methods cannot provide for studying the Bcr-Abl imatinib resistance problem, especially in the case of conformational changes due to remote mutations. By utilizing the Molecular Mechanics/Poisson-Boltzmann surface area (MM-PBSA) techniques and analyzing the interactions between imatinib and individual residues, imatinib resistance mechanisms not previously thought have been revealed. Non-directly contacted P-loop mutations either unfavorably change the direct electrostatic interactions with imatinib, or cause the conformational changes influencing the contact energies between imatinib and other non-P-loop residues. We demonstrate that imatinib resistance of T315I mainly comes from the breakdown of the interactions between imatinib and E286 and M290, contradictory to previously suggested that the missing hydrogen bonding is the main contribution. We also demonstrate that except for the mutations of the direct contact residues, such as L248 and Y253, the unfavorable electrostatic interaction between P-loop and imatinib is the main reason for resistance for the P-loop mutations. Furthermore, in Y255H, protonation of the histidin is essential for rendering this mutation resistant to Gleevec. Our results demonstrate that MD is a powerful way to verify and predict clinical response or resistance to imatinib and other potential drugs.

scholarly journals Molecular Dynamics Simulations of a Catalytic Multivalent Peptide–Nanoparticle Complex

2021 ◽  
Vol 22 (7) ◽  
pp. 3624
Author(s):  
Sutapa Dutta ◽  
Stefano Corni ◽  
Giorgia Brancolini
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Conformational Changes ◽  
Self Assembly ◽  
Brownian Dynamics ◽  
Rational Design ◽  
Complex Structure ◽  
Gold Nanocluster ◽  
Field Development ◽  
Dynamics Simulations

Molecular modeling of a supramolecular catalytic system is conducted resulting from the assembling between a small peptide and the surface of cationic self-assembled monolayers on gold nanoparticles, through a multiscale iterative approach including atomistic force field development, flexible docking with Brownian Dynamics and µs-long Molecular Dynamics simulations. Self-assembly is a prerequisite for the catalysis, since the catalytic peptides do not display any activity in the absence of the gold nanocluster. Atomistic simulations reveal details of the association dynamics as regulated by defined conformational changes of the peptide due to peptide length and sequence. Our results show the importance of a rational design of the peptide to enhance the catalytic activity of peptide–nanoparticle conjugates and present a viable computational approach toward the design of enzyme mimics having a complex structure–function relationship, for technological and nanomedical applications.

scholarly journals Effects of number of parallel runs and frequency of bias-strength replacement in generalized ensemble molecular dynamics simulations

2019 ◽  
Vol 1 ◽  
pp. e4 ◽  
Author(s):  
Takuya Shimato ◽  
Kota Kasahara ◽  
Junichi Higo ◽  
Takuya Takahashi
Keyword(s):  
Molecular Dynamics ◽  
Potential Energy ◽  
Molecular Dynamics Simulations ◽  
Conformational Changes ◽  
High Frequency ◽  
Md Simulations ◽  
Energy Space ◽  
Bias Potential ◽  
Ensemble Approach ◽  
Dynamics Simulations

Background The generalized ensemble approach with the molecular dynamics (MD) method has been widely utilized. This approach usually has two features. (i) A bias potential, whose strength is replaced during a simulation, is applied. (ii) Sampling can be performed by many parallel runs of simulations. Although the frequency of the bias-strength replacement and the number of parallel runs can be adjusted, the effects of these settings on the resultant ensemble remain unclear. Method In this study, we performed multicanonical MD simulations for a foldable mini-protein (Trp-cage) and two unstructured peptides (8- and 20-residue poly-glutamic acids) with various settings. Results As a result, running many short simulations yielded robust results for the Trp-cage model. Regarding the frequency of the bias-potential replacement, although using a high frequency enhanced the traversals in the potential energy space, it did not promote conformational changes in all the systems.

scholarly journals Molecular dynamics simulations disclose early stages of the photo-activation of cryptochrome 4

2018 ◽  
Author(s):  
D. R. Kattnig ◽  
C. Nielsen ◽  
I. A. Solov’yov
Keyword(s):  
Molecular Dynamics ◽  
Conformational Changes ◽  
Large Scale ◽  
Conformational Dynamics ◽  
Md Simulations ◽  
Quantum Processes ◽  
Helical Segment ◽  
High Betweenness ◽  
Dynamics Simulations ◽  
Fad Binding

AbstractBirds appear to be equipped with a light-dependent, radical-pair-based magnetic compass that relies on truly quantum processes. While the identity of the sensory protein has remained speculative, cryptochrome 4 has recently been identified as the most auspicious candidate. Here, we report on allatom molecular dynamics (MD) simulations addressing the structural reorganisations that accompany the photoreduction of the flavin cofactor in a model of the European robin cryptochrome 4 (ErCry4). Extensive MD simulations reveal that the photo-activation of ErCry4 induces large-scale conformational changes on short (hundreds of nanoseconds) timescales. Specifically, the photo-reduction is accompanied with the release of the C-terminal tail, structural rearrangements in the vicinity of the FAD-binding site, and the noteworthy formation of an α-helical segment at the N-terminal part. Some of these rearrangements appear to expose potential phosphorylation sites. We describe the conformational dynamics of the protein using a graph-based approach that is informed by the adjacency of residues and the correlation of their local motions. This approach reveals densely coupled reorganisation communities, which facilitate an efficient signal transduction due to a high density of hubs. These communities are interconnected by a small number of highly important residues characterized by high betweenness centrality. The network approach clearly identifies the sites restructuring upon photoactivation, which appear as protrusions or delicate bridges in the reorganisation network. We also find that, unlike in the homologous cryptochrome from D. melanogaster, the release of the C-terminal domain does not appear to be correlated with the transposition of a histidine residue close to the FAD cofactor.

scholarly journals An equilibrium constant pH molecular dynamics method for accurate prediction of pH-dependence in protein systems: Theory and application

2020 ◽  
Author(s):  
Ekaterina Kots ◽  
Derek M. Shore ◽  
Harel Weinstein
Keyword(s):  
Molecular Dynamics ◽  
Equilibrium Constant ◽  
Conformational Changes ◽  
Ph Dependence ◽  
Md Simulations ◽  
Water Soluble ◽  
Conformational Sampling ◽  
Ph Change ◽  
Constant Ph ◽  
Dynamics Simulations

ABSTRACTComputational modeling and simulation of biomolecular systems at their functional pH ranges requires an accurate approach to exploring the pH dependence of conformations and interactions. Here we present a new approach – the Equilibrium Constant pH (ECpH) method – to perform conformational sampling of protein systems in the framework of molecular dynamics simulations in an N, P, T-thermodynamic ensemble. The performance of ECpH is illustrated for two proteins with experimentally determined conformational responses to pH change: the small globular water-soluble bovine b-lactoglobulin (BBL), and the dimer transmembrane antiporter CLC-ec1 Cl−/H+. We show that with computational speeds comparable to equivalent canonical MD simulations we performed, the ECpH trajectories reproduce accurately the pH-dependent conformational changes observed experimentally in these two protein systems, some of which were not seen in the corresponding canonical MD simulations.Abstract FigureTable of Contents artwork

scholarly journals Salt bridges gate α-catenin activation at intercellular junctions

2018 ◽  
Vol 29 (2) ◽  
pp. 111-122 ◽  
Author(s):  
Samantha Barrick ◽  
Jing Li ◽  
Xinyu Kong ◽  
Alokananda Ray ◽  
Emad Tajkhorshid ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Mechanical Activation ◽  
Conformational Changes ◽  
Intercellular Junctions ◽  
Salt Bridge ◽  
Salt Bridges ◽  
Equilibrium Binding ◽  
Binding Data ◽  
Dynamics Simulations ◽  
Intercellular Adhesions

Molecular dynamics simulations, equilibrium binding measurements, and fluorescence imaging reveal the influence of a key salt bridge in the mechanical activation of α-catenin at intercellular adhesions. Simulations reveal possible α-catenin conformational changes underlying experimental fluorescence and equilibrium binding data.

scholarly journals Structural Changes of Sarco/Endoplasmic Reticulum Ca2+-ATPase Induced by Rutin Arachidonate: A Molecular Dynamics Study

Biomolecules ◽  
2020 ◽  
Vol 10 (2) ◽  
pp. 214
Author(s):  
Yoel Rodríguez ◽  
Magdaléna Májeková
Keyword(s):  
Molecular Dynamics ◽  
Endoplasmic Reticulum ◽  
Conformational Changes ◽  
Calcium Ions ◽  
Structural Changes ◽  
Acid Group ◽  
Bilayer Membrane ◽  
Salt Bridges ◽  
Drug Candidates ◽  
Dynamics Simulations

Sarco/endoplasmic reticulum Ca2+-ATPase (SERCA) maintains the level of calcium concentration in cells by pumping calcium ions from the cytoplasm to the lumen while undergoing substantial conformational changes, which can be stabilized or prevented by various compounds. Here we attempted to clarify the molecular mechanism of action of new inhibitor rutin arachidonate, one of the series of the acylated rutin derivatives. We performed molecular dynamics simulations of SERCA1a protein bound to rutin arachidonate positioned in a pure dipalmitoylphosphatidylcholine bilayer membrane. Our study predicted the molecular basis for the binding of rutin arachidonate towards SERCA1a in the vicinity of the binding site of calcium ions and near the location of the well-known inhibitor thapsigargin. The stable hydrogen bond between Glu771 and rutin arachidonate plays a key role in the binding. SERCA1a is kept in the E2 conformation preventing the formation of important salt bridges between the side chains of several residues, primarily Glu90 and Lys297. All in all, the structural changes induced by the binding of rutin arachidonate to SERCA1a may shift proton balance near the titrable residues Glu771 and Glu309 into neutral species, hence preventing the binding of calcium ions to the transmembrane binding sites and thus affecting calcium homeostasis. Our results could lead towards the design of new types of inhibitors, potential drug candidates for cancer treatment, which could be anchored to the transmembrane region of SERCA1a by a lipophilic fatty acid group.

scholarly journals Exploring Large Domain Motions in Proteins Using Atomistic Molecular Dynamics with Enhanced Conformational Sampling

2020 ◽  
Vol 22 (1) ◽  
pp. 270
Author(s):  
Hisham M. Dokainish ◽  
Yuji Sugita
Keyword(s):  
Molecular Dynamics ◽  
Conformational Changes ◽  
Md Simulations ◽  
Conformational Sampling ◽  
Salt Bridges ◽  
Multidomain Proteins ◽  
Charged Residues ◽  
Ligand Interactions ◽  
Moving Domains ◽  
Domain Motions

Conformational transitions in multidomain proteins are essential for biological functions. The Apo conformations are typically open and flexible, while the Holo states form more compact conformations stabilized by protein-ligand interactions. Unfortunately, the atomically detailed mechanisms for such open-closed conformational changes are difficult to be accessed experimentally as well as computationally. To simulate the transitions using atomistic molecular dynamics (MD) simulations, efficient conformational sampling algorithms are required. In this work, we propose a new approach based on generalized replica-exchange with solute tempering (gREST) for exploring the open-closed conformational changes in multidomain proteins. Wherein, selected surface charged residues in a target protein are defined as the solute region in gREST simulation and the solute temperatures are different in replicas and exchanged between them to enhance the domain motions. This approach is called gREST selected surface charged residues (gREST_SSCR) and is applied to the Apo and Holo states of ribose binding protein (RBP) in solution. The conformational spaces sampled with gREST_SSCR are much wider than those with the conventional MD, sampling open-closed conformational changes while maintaining RBP domains’ stability. The free-energy landscapes of RBP in the Apo and Holo states are drawn along with twist and hinge angles of the two moving domains. The inter-domain salt-bridges that are not observed in the experimental structures are also important in the intermediate states during the conformational changes.

Crystal structures of Magnaporthe oryzae trehalose-6-phosphate synthase (MoTps1) suggest a model for catalytic process of Tps1

2019 ◽  
Vol 476 (21) ◽  
pp. 3227-3240 ◽  
Author(s):  
Shanshan Wang ◽  
Yanxiang Zhao ◽  
Long Yi ◽  
Minghe Shen ◽  
Chao Wang ◽  
...  
Keyword(s):  
Crystal Structures ◽  
Conformational Changes ◽  
Magnaporthe Oryzae ◽  
Structural Information ◽  
Complex Structure ◽  
Critical Role ◽  
Catalytic Process ◽  
Structural Basis ◽  
Salt Bridges ◽  
Terminal Domain

Trehalose-6-phosphate (T6P) synthase (Tps1) catalyzes the formation of T6P from UDP-glucose (UDPG) (or GDPG, etc.) and glucose-6-phosphate (G6P), and structural basis of this process has not been well studied. MoTps1 (Magnaporthe oryzae Tps1) plays a critical role in carbon and nitrogen metabolism, but its structural information is unknown. Here we present the crystal structures of MoTps1 apo, binary (with UDPG) and ternary (with UDPG/G6P or UDP/T6P) complexes. MoTps1 consists of two modified Rossmann-fold domains and a catalytic center in-between. Unlike Escherichia coli OtsA (EcOtsA, the Tps1 of E. coli), MoTps1 exists as a mixture of monomer, dimer, and oligomer in solution. Inter-chain salt bridges, which are not fully conserved in EcOtsA, play primary roles in MoTps1 oligomerization. Binding of UDPG by MoTps1 C-terminal domain modifies the substrate pocket of MoTps1. In the MoTps1 ternary complex structure, UDP and T6P, the products of UDPG and G6P, are detected, and substantial conformational rearrangements of N-terminal domain, including structural reshuffling (β3–β4 loop to α0 helix) and movement of a ‘shift region' towards the catalytic centre, are observed. These conformational changes render MoTps1 to a ‘closed' state compared with its ‘open' state in apo or UDPG complex structures. By solving the EcOtsA apo structure, we confirmed that similar ligand binding induced conformational changes also exist in EcOtsA, although no structural reshuffling involved. Based on our research and previous studies, we present a model for the catalytic process of Tps1. Our research provides novel information on MoTps1, Tps1 family, and structure-based antifungal drug design.

Sign in / Sign up

Export Citation Format

Share Document