Structural basis of single molecular heparin-FX06 interaction revealed by SPM measurements and molecular simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Guo C
PROVIDER: S-EPMC3508726 | biostudies-literature | 2012 Dec
REPOSITORIES: biostudies-literature
Publications
Structural basis of single molecular heparin-FX06 interaction revealed by SPM measurements and molecular simulations.
Guo Cunlan C Wang Bin B Wang Lianchun L Xu Bingqian B
Chemical communications (Cambridge, England) 20121201 100
Heparin, a functionalized polysaccharide, is observed under a scanning tunneling microscope, which shows atomic scale conformational details. The peptide FX06 is found to bind to five consecutive sugar units of heparin and this interaction is directly revealed by atomic force microscopy and dynamic force spectroscopy measurements. The determined free energy change agrees well with the dynamic calculation result. ...[more]
Similar Datasets
Project description:The WRKY transcription factors are a class of DNA-binding proteins involved in diverse plant processes play critical roles in response to abiotic and biotic stresses. Genome-wide divergence analysis of WRKY gene family in Hordeum vulgare provided a framework for molecular evolution and functional roles. So far, the crystal structure of WRKY from barley has not been resolved; moreover, knowledge of the three-dimensional structure of WRKY domain is pre-requisites for exploring the protein-DNA recognition mechanisms. Homology modelling based approach was used to generate structures for WRKY DNA binding domain (DBD) and its variants using AtWRKY1 as a template. Finally, the stability and conformational changes of the generated model in unbound and bound form was examined through atomistic molecular dynamics (MD) simulations for 100 ns time period.In this study, we investigated the comparative binding pattern of WRKY domain and its variants with W-box cis-regulatory element using molecular docking and dynamics (MD) simulations assays. The atomic insight into WRKY domain exhibited significant variation in the intermolecular hydrogen bonding pattern, leading to the structural anomalies in the variant type and differences in the DNA-binding specificities. Based on the MD analysis, residual contribution and interaction contour, wild-type WRKY (HvWRKY46) were found to interact with DNA through highly conserved heptapeptide in the pre- and post-MD simulated complexes, whereas heptapeptide interaction with DNA was missing in variants (I and II) in post-MD complexes. Consequently, through principal component analysis, wild-type WRKY was also found to be more stable by obscuring a reduced conformational space than the variant I (HvWRKY34). Lastly, high binding free energy for wild-type and variant II allowed us to conclude that wild-type WRKY-DNA complex was more stable relative to variants I.The results of our study revealed complete dynamic and structural information about WRKY domain-DNA interactions. However, no structure base information reported to date for WRKY variants and their mechanism of interaction with DNA. Our findings highlighted the importance of selecting a sequence to generate newer transgenic plants that would be increasingly tolerance to stress conditions.
| S-EPMC5810047 | biostudies-literature
Project description:Bcl-X(L), an antiapoptotic Bcl-2 family protein, plays a central role in the regulation of the apoptotic pathway. Heterodimerization of the antiapoptotic Bcl-2 family proteins with the proapoptotic family members such as Bad, Bak, Bim and Bid is a crucial step in the apoptotic regulation. In addition to these conventional binding partners, recent evidences reveal that the Bcl-2 family proteins also interact with noncanonical binding partners such as p53. Our previous NMR studies showed that Bcl-X(L): BH3 peptide and Bcl-X(L): SN15 peptide (a peptide derived from residues S15-N29 of p53) complex structures share similar modes of bindings. To further elucidate the molecular basis of the interactions, here we have employed molecular dynamics simulations coupled with MM/PBSA approach. Bcl-X(L) and other Bcl-2 family proteins have 4 hydrophobic pockets (p1-p4), which are occupied by four systematically spaced hydrophobic residues (h1-h4) of the proapoptotic Bad and Bak BH3 peptides. We observed that three conserved hydrophobic residues (F19, W23 and L26) of p53 (SN15) peptide anchor into three hydrophobic pockets (p2-p4) of Bcl-X(L) in a similar manner as BH3 peptide. Our results provide insights into the novel molecular recognition by Bcl-X(L) with p53.
| S-EPMC3198449 | biostudies-literature
Project description:The first KRAS(G12C) targeting inhibitor in clinical development, AMG 510, has shown promising antitumor activity in clinical trials. On the molecular level, however, the interaction dynamics of this covalently bound drug-protein complex has been undetermined. Here, we disclose the interaction dynamics of the KRAS(G12C)-AMG 510 complex by long timescale all-atom molecular dynamics (MD) simulations (total of 75 μs). Moreover, we investigated the influence of the recently reported post-translational modification (PTM) of KRAS' N-terminus, removal of initiator methionine (iMet1) with acetylation of Thr2, to this complex. Our results demonstrate that AMG 510 does not entrap KRAS into a single conformation, as one would expect based on the crystal structure, but rather into an ensemble of conformations. AMG 510 binding is extremely stable regardless of highly dynamic interface of KRAS' switches. Overall, KRAS(G12C)-AMG 510 complex partially mimic the native dynamics of GDP bound KRAS; however, AMG 510 stabilizes the α3-helix region. N-terminally modified KRAS displays similar interaction dynamics with AMG 510 as when Met1 is present, but this PTM appears to stabilize β2-β3-loop. These results provide novel conformational insights on the molecular level to KRAS(G12C)-AMG 510 interactions and dynamics, providing new perspectives to RAS related drug discovery.
| S-EPMC7371895 | biostudies-literature
Project description:Relating ion channel (iCh) structural dynamics to physiological function remains a challenge. Current experimental and computational techniques have limited ability to explore this relationship in atomistic detail over physiological timescales. A framework associating iCh structure to function is necessary for elucidating normal and disease mechanisms. We formulated a modeling schema that overcomes the limitations of current methods through applications of artificial intelligence machine learning. Using this approach, we studied molecular processes that underlie human IKs voltage-mediated gating. IKs malfunction underlies many debilitating and life-threatening diseases. Molecular components of IKs that underlie its electrophysiological function include KCNQ1 (a pore-forming tetramer) and KCNE1 (an auxiliary subunit). Simulations, using the IKs structure-function model, reproduced experimentally recorded saturation of gating-charge displacement at positive membrane voltages, two-step voltage sensor (VS) movement shown by fluorescence, iCh gating statistics, and current-voltage relationship. Mechanistic insights include the following: 1) pore energy profile determines iCh subconductance; 2) the entire protein structure, not limited to the pore, contributes to pore energy and channel subconductance; 3) interactions with KCNE1 result in two distinct VS movements, causing gating-charge saturation at positive membrane voltages and current activation delay; and 4) flexible coupling between VS and pore permits pore opening at lower VS positions, resulting in sequential gating. The new modeling approach is applicable to atomistic scale studies of other proteins on timescales of physiological function.
| S-EPMC6129186 | biostudies-literature
Project description:Human hemoglobin (Hb) is a benchmark protein of structural biology that shaped our view of allosterism over 60 years ago, with the introduction of the MWC model based on Perutz structures of the oxy(R) and deoxy(T) states and the more recent Tertiary Two-State model that proposed the existence of individual subunit states -"r" and "t"-, whose structure is yet unknown. Cooperative oxygen binding is essential for Hb function, and despite decades of research there are still open questions related to how tertiary and quaternary changes regulate oxygen affinity. In the present work, we have determined the free energy profiles of oxygen migration and for HisE7 gate opening, with QM/MM calculations of the oxygen binding energy in order to address the influence of tertiary differences in the control of oxygen affinity. Our results show that in the α subunit the low to high affinity transition is achieved by a proximal effect that mostly affects oxygen dissociation and is the driving force of the allosteric transition, while in the β subunit the affinity change results from a complex interplay of proximal and distal effects, including an increase in the HE7 gate opening, that as shown by free energy profiles promotes oxygen uptake.
| S-EPMC5589765 | biostudies-literature
Project description:Polyplexes composed of polyethyleneimine (PEI) and DNA or siRNA have attracted great attention for their use in gene therapy. Although many physicochemical characteristics of these polyplexes remain unknown, PEI/DNA complexes have been repeatedly shown to be more stable than their PEI/siRNA counterparts. Here, we examine potential causes for this difference using atomistic molecular dynamics simulations of complexation between linear PEI and DNA or siRNA duplexes containing the same number of bases. The two types of polyplexes are stabilized by similar interactions, as PEI amines primarily interact with nucleic acid phosphate groups but also occasionally interact with groove atoms of both nucleic acids. However, the number of interactions in PEI/DNA complexes is greater than in comparable PEI/siRNA complexes, with interactions between protonated PEI amines and DNA being particularly enhanced. These results indicate that structural differences between DNA and siRNA may play a role in the increased stability of PEI/DNA complexes. In addition, we investigate the binding of PEI chains to polyplexes that have a net positive charge. The binding of PEI to these overcharged complexes involves interactions between PEI and areas on the nucleic acid surface that have maintained a negative electrostatic potential and is facilitated by the release of water from the nucleic acid.
| S-EPMC5677264 | biostudies-literature
Project description:A mechanical view provides an attractive alternative for predicting the behavior of complex systems since it circumvents the resource-intensive requirements of atomistic models; however, it remains extremely challenging to characterize the mechanical responses of a system at the molecular level. Here, the structural distribution is proposed to be an effective means to extracting the molecular mechanical properties. End-to-end distance distributions for a series of short poly-L-proline peptides with the sequence P(n)CG(3)K-biotin (n = 8, 12, 15 and 24) were used to experimentally illustrate this new approach. High-resolution single-molecule Förster-type resonance energy transfer (FRET) experiments were carried out and the conformation-resolving power was characterized and discussed in the context of the conventional constant-time binning procedure for FRET data analysis. It was shown that the commonly adopted theoretical polymer models-including the worm-like chain, the freely jointed chain, and the self-avoiding chain-could not be distinguished by the averaged end-to-end distances, but could be ruled out using the molecular details gained by conformational distribution analysis because similar polymers of different sizes could respond to external forces differently. Specifically, by fitting the molecular conformational distribution to a semi-flexible polymer model, the effective persistence lengths for the series of short poly-L-proline peptides were found to be size-dependent with values of ~190 Å, ~67 Å, ~51 Å, and ~76 Å for n = 8, 12, 15, and 24, respectively. A comprehensive computational modeling was carried out to gain further insights for this surprising discovery. It was found that P(8) exists as the extended all-trans isomaer whereas P(12) and P(15) predominantly contained one proline residue in the cis conformation. P(24) exists as a mixture of one-cis (75%) and two-cis (25%) isomers where each isomer contributes to an experimentally resolvable conformational mode. This work demonstrates the resolving power of the distribution-based approach, and the capacity of integrating high-resolution single-molecule FRET experiments with molecular modeling to reveal detailed structural information about the conformation of molecules on the length scales relevant to the study of biological molecules.
| S-EPMC3361908 | biostudies-literature
Project description:The hemostatic function of von Willebrand factor is downregulated by the metalloprotease ADAMTS13, which cleaves at a unique site normally buried in the A2 domain. Exposure of the proteolytic site is induced in the wild-type by shear stress as von Willebrand factor circulates in blood. Mutations in the A2 domain, which increase its susceptibility to cleavage, cause type 2A von Willebrand disease. In this study, molecular dynamics simulations suggest that the A2 domain unfolds under tensile force progressively through a series of steps. The simulation results also indicated that three type 2A mutations in the C-terminal half of the A2 domain, L1657I, I1628T and E1638K, destabilize the native state fold of the protein. Furthermore, all three type 2A mutations lowered in silico the tensile force necessary to undock the C-terminal helix α6 from the rest of the A2 domain, the first event in the unfolding pathway. The mutations F1520A, I1651A and A1661G were also predicted by simulations to destabilize the A2 domain and facilitate exposure of the cleavage site. Recombinant A2 domain proteins were expressed and cleavage assays were performed with the wild-type and single-point mutants. All three type 2A and two of the three predicted mutations exhibited increased rate of cleavage by ADAMTS13. These results confirm that destabilization of the helix α6 in the A2 domain facilitates exposure of the cleavage site and increases the rate of cleavage by ADAMTS13.
| S-EPMC3479114 | biostudies-literature
Project description:Given that immunoproteasome inhibitors are currently being developed for a variety of potent therapeutic purposes, the unique specificity of an ?',?'-epoxyketone peptide (UK101) toward the LMP2 subunit of the immunoproteasome (analogous to ?5 subunit of the constitutive proteasome) has been investigated in this study for the first time by employing homology modeling, molecular docking, molecular dynamics simulation, and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) binding free energy calculations. On the basis of the simulated binding structures, the calculated binding free energies are in qualitative agreement with the corresponding experimental data, and the selectivity of UK101 is explained reasonably. The observed selectivity of UK101 for the LMP2 subunit is rationalized by the requirement for both a linear hydrocarbon chain at the N terminus and a bulky group at the C terminus of the inhibitor, because the LMP2 subunit has a much more favorable hydrophobic pocket interacting with the linear hydrocarbon chain, and the bulky group at the C terminus has a steric clash with the Tyr 169 in ?5 subunit. Finally, our results help to clarify why UK101 is specific to the LMP2 subunit of immunoproteasome, and this investigation should be valuable for rational design of more potent LMP2-specific inhibitors.
| S-EPMC2945215 | biostudies-literature
Project description:The point mutations M205S and M205R have been demonstrated to severely disturb the folding and maturation process of the cellular prion protein (PrP(C)). These disturbances have been interpreted as consequences of mutation-induced structural changes in PrP, which are suggested to involve helix 1 and its attachment to helix 3, because the mutated residue M205 of helix 3 is located at the interface of these two helices. Furthermore, current models of the prion protein scrapie (PrP(Sc)), which is the pathogenic isoform of PrP(C) in prion diseases, imply that helix 1 disappears during refolding of PrP(C) into PrP(Sc). Based on molecular-dynamics simulations of wild-type and mutant PrP(C) in aqueous solution, we show here that the native PrP(C) structure becomes strongly distorted within a few nanoseconds, once the point mutations M205S and M205R have been applied. In the case of M205R, this distortion is characterized by a motion of helix 1 away from the hydrophobic core into the aqueous environment and a subsequent structural decay. Together with experimental evidence on model peptides, this decay suggests that the hydrophobic attachment of helix 1 to helix 3 at M205 is required for its correct folding into its stable native structure.
| S-EPMC1459491 | biostudies-literature