A detailed picture of a protein-carbohydrate hydrogen-bonding network revealed by NMR and MD simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Nestor G
PROVIDER: S-EPMC8091458 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The gap junction protein connexin43 (Cx43) binds to the second PDZ domain of Zonula occludens-1 (ZO-1) through its C-terminal tail, mediating the regulation of gap junction plaque size and dynamics. Biochemical study demonstrated that the very C-terminal 12 residues of Cx43 are necessary and sufficient for ZO-1 PDZ2 binding and phosphorylation at residues Ser (-9) and Ser (-10) of the peptide can disrupt the association. However, only a crystal structure of ZO-1 PDZ2 in complex with a shorter 9 aa peptide of connexin43 was solved experimentally. Here, the interactions between ZO-1 PDZ2 and the short, long and phosphorylated Cx43 peptides were studied using molecular dynamics (MD) simulations and free energy calculation. The short peptide bound to PDZ2 exhibits large structural variations, while the extension of three upstream residues stabilizes the peptide conformation and enhanced the interaction. Phosphorylation at Ser(-9) significantly weakens the binding and results in conformational flexibility of the peptide. Glu210 of ZO-1 PDZ2 was found to be a key regulatory point in Cx43 binding and phosphorylation induced dissociation.
| S-EPMC3121883 | biostudies-literature
Project description:BACKGROUND:The role of dynamics in protein functions including signal transduction is just starting to be deciphered. Eph receptors with 16 members divided into A- and B- subclasses are respectively activated by 9 A- and B-ephrin ligands. EphA4 is the only receptor capable of binding to all 9 ephrins and small molecules with overlapped interfaces. RESULTS:We first determined the structures of the EphA4 ligand binding domain (LBD) in two crystals of P1 space group. Noticeably, 8 EphA4 molecules were found in one asymmetric unit and consequently from two crystals we obtained 16 structures, which show significant conformational variations over the functionally critical A-C, D-E, G-H and J-K loops. The 16 new structures, together with previous 9 ones, can be categorized into two groups: closed and open forms which resemble the uncomplexed and complexed structures of the EphA4 LBD respectively. To assess whether the conformational diversity over the loops primarily results from the intrinsic dynamics, we initiated 30-ns molecular dynamics (MD) simulations for both closed and open forms. The results indicate that the loops do have much higher intrinsic dynamics, which is further unravelled by NMR H/D exchange experiments. During simulations, the open form has the RMS deviations slightly larger than those of the closed one, suggesting the open form may be less stable in the absence of external contacts. Furthermore, no obvious exchange between two forms is observed within 30 ns, implying that they are dynamically separated. CONCLUSIONS:Our study provides the first experimental and computational result revealing that the intrinsic dynamics are most likely underlying the conformational diversity observed for the EphA4 LBD loops mediating the binding affinity and specificity. Interestingly, the open conformation of the EphA4 LBD is slightly unstable in the absence of it natural ligand ephrins, implying that the conformational transition from the closed to open has to be driven by the high-affinity interaction with ephrins because the weak interaction with small molecule was found to be insufficient to trigger the transition. Our results therefore highlight the key role of protein dynamics in Eph-ephrin signalling and would benefit future design of agonists/antagonists targeting Eph receptors.
| S-EPMC3274464 | biostudies-literature
Project description:Hydrogen exchange (HX) studies have provided critical insight into our understanding of protein folding, structure, and dynamics. More recently, hydrogen exchange mass spectrometry (HX-MS) has become a widely applicable tool for HX studies. The interpretation of the wealth of data generated by HX-MS experiments as well as other HX methods would greatly benefit from the availability of exchange predictions derived from structures or models for comparison with experiment. Most reported computational HX modeling studies have employed solvent-accessible-surface-area based metrics in attempts to interpret HX data on the basis of structures or models. In this study, a computational HX-MS prediction method based on classification of the amide hydrogen bonding modes mimicking the local unfolding model is demonstrated. Analysis of the NH bonding configurations from molecular dynamics (MD) simulation snapshots is used to determine partitioning over bonded and nonbonded NH states and is directly mapped into a protection factor (PF) using a logistics growth function. Predicted PFs are then used for calculating deuteration values of peptides and compared with experimental data. Hydrogen exchange MS data for fatty acid synthase thioesterase (FAS-TE) collected for a range of pHs and temperatures was used for detailed evaluation of the approach. High correlation between prediction and experiment for observable fragment peptides is observed in the FAS-TE and additional benchmarking systems that included various apo/holo proteins for which literature data were available. In addition, it is shown that HX modeling can improve experimental resolution through decomposition of in-exchange curves into rate classes, which correlate with prediction from MD. Successful rate class decompositions provide further evidence that the presented approach captures the underlying physical processes correctly at the single residue level. This assessment is further strengthened in a comparison of residue resolved protection factor predictions for staphylococcal nuclease with NMR data, which was also used to compare prediction performance with other algorithms described in the literature. The demonstrated transferable and scalable MD based HX prediction approach adds significantly to the available tools for HX-MS data interpretation based on available structures and models.
| S-EPMC4587094 | biostudies-literature
Project description:Independent experimental and computational approaches show agreement concerning arginine/membrane interactions when a single arginine is introduced at selected positions within the membrane-spanning region of acetyl-GGALW(5)LALALAL(12)AL(14)ALALW(19)LAGA-ethanolamide, designated GWALP23. Peptide sequence isomers having Arg in position 12 or position 14 display markedly different behaviors, as deduced by both solid-state NMR experiments and coarse-grained molecular dynamics (CG-MD) simulations. With respect to the membrane normal of DOPC or DPPC lipid bilayer membranes, GWALP23-R14 shows one major state whose apparent average tilt is approximately 10 degrees greater than that of GWALP23. The presence of R14 furthermore induces bilayer thinning and peptide displacement to "lift" the charged guanidinium toward the bilayer surface. By contrast, GWALP23-R12 exhibits multiple states that are in slow exchange on the NMR time scale, with CG-MD simulations indicating two distinct positions with different screw rotation angles in the membrane, along with an increased tendency to exit the lipid bilayer.
| S-EPMC2888137 | biostudies-literature
Project description:Ceramides play a key modulatory role in many cellular processes, which results from their effect on the structure and dynamics of biological membranes. In this study, we investigate the influence of C16-ceramide (C16) on the biophysical properties of DMPC lipid bilayers using solid-state NMR and atomistic molecular dynamics (MD) simulations. MD simulations and NMR measurements were carried out for a pure DMPC bilayer and for a 20% DMPC-C16 mixture. Calculated key structural properties, namely area per lipid, chain order parameters, and mass density profiles, indicate that C16 has an ordering effect on the DMPC bilayer. Furthermore, the simulations predict that specific hydrogen-bonds form between DMPC and C16 molecules. Multi-nuclear solid-state NMR was used to verify these theoretical predictions. Chain order parameters extracted from (13)C(1)H dipole couplings were measured for both lipid and ceramide and follow the trend suggested by the MD simulations. Furthermore, (1)H-MAS NMR experiments showed a direct contact between ceramide and lipids.
| S-EPMC4137758 | biostudies-literature
Project description:Variable domains of camelid antibodies (so-called nanobodies or VHH) are the smallest antibody fragments that retain complete functionality and therapeutic potential. Understanding of the nanobody-binding interface has become a pre-requisite for rational antibody design and engineering. The nanobody-binding interface consists of up to three hypervariable loops, known as the CDR loops. Here, we structurally and dynamically characterize the conformational diversity of an anti-GFP-binding nanobody by using molecular dynamics simulations in combination with experimentally derived data from nuclear magnetic resonance (NMR) spectroscopy. The NMR data contain both structural and dynamic information resolved at various timescales, which allows an assessment of the quality of protein MD simulations. Thus, in this study, we compared the ensembles for the anti-GFP-binding nanobody obtained from MD simulations with results from NMR. We find excellent agreement of the NOE-derived distance maps obtained from NMR and MD simulations and observe similar conformational spaces for the simulations with and without NOE time-averaged restraints. We also compare the measured and calculated order parameters and find generally good agreement for the motions observed in the ps-ns timescale, in particular for the CDR3 loop. Understanding of the CDR3 loop dynamics is especially critical for nanobodies, as this loop is typically critical for antigen recognition.
| S-EPMC9141556 | biostudies-literature
Project description:Molecular dynamics (MD) simulations and nuclear magnetic resonance spin-relaxation measurements provide detailed insights into ps-ns structural dynamics of proteins. An analysis of discrepancies between the two methods is presented for the B3 immunoglobulin-binding domain of streptococcal protein G. MD simulations using three MD force fields (OPLS-AA, AMBER ff99SB, and AMBER ff03) overestimate the flexibility of backbone N--H vectors at the borders of secondary structure and in loops when compared with experimentally determined backbone amide generalized order parameters (Hall and Fushman, J Am Chem Soc 2006; 12:7855-7870). Comparison with a previous study of residual dipolar coupling constants (Bouvignies et al., Proc Natl Acad Sci USA 2005;102:13885-13890) indicates that slower timescale motions do not account for the discrepancies. Structural analysis reveals that relative imbalance between the description of hydrogen bonding and other terms of modern force fields may be responsible for disagreement.
| S-EPMC2750878 | biostudies-literature
Project description:BackgroundHow proteins approach surrounding molecules is fundamental to our understanding of the specific interactions that occur at the surface of proteins. The enhanced surface accessibility of small molecules such as organic solvents and paramagnetic probes to protein binding sites has been observed; however, the molecular basis of this finding has not been fully established. Recently, it has been suggested that hydration dynamics play a predominant role in controlling the distribution of hot spots on surface of proteins.ResultsIn the present study, the hydration of the archaeal multifunctional protein Sso7d from Solfolobus solfataricus was investigated using a combination of computational and experimental data derived from molecular dynamics simulations and ePHOGSY NMR spectroscopy.ConclusionsWe obtained a convergent protein hydration landscape that indicated how the shape and stability of the Sso7d hydration shell could modulate the function of the protein. The DNA binding domain overlaps with the protein region involved in chaperon activity and this domain is hydrated only in a very small central region. This localized hydration seems to favor intermolecular approaches from a large variety of ligands. Conversely, high water density was found in surface regions of the protein where the ATP binding site is located, suggesting that surface water molecules play a role in protecting the protein from unspecific interactions.
| S-EPMC3207888 | biostudies-literature
Project description:Peptides play a critical role in the life of organisms, performing completely different functions. The biological activity of some peptides, such as cyclosporins, can be determined by the degree of membrane permeability. Thus, it becomes important to study how the molecule interacts with lipid bilayers. Cyclosporins C, E, H and L were characterised molecular dynamics simulation; NMR spectroscopy studies were also carried out for cyclosporins C and E. The comparison of one- and two-dimensional spectra revealed certain similarities between spatial structures of the studied cyclosporin variants. Upon dissolving in water containing DPC micelles, which serve as model membranes, subtle changes in the NMR spectra appear, but in a different way for different cyclosporins. In order to understand whether observed changes are related to any structural modifications, simulation of the interaction of the peptide with the phospholipid micelle was performed. The onset of the interaction was observed, when the peptide is trapped to the surface of the micelle. Simulations of this kind are also of interest in the light of the well-known membrane permeability of cyclosporin, which is important for its biological action.
| S-EPMC9965255 | biostudies-literature
Project description:Interactions of charged molecules with biomembranes regulate many of their biological activities, but their binding affinities to lipid bilayers are difficult to measure experimentally and model theoretically. Classical molecular dynamics (MD) simulations have the potential to capture the complex interactions determining how charged biomolecules interact with membranes, but systematic overbinding of sodium and calcium cations in standard MD simulations raises the question of how accurately force fields capture the interactions between lipid membranes and charged biomolecules. Here, we evaluate the binding of positively charged small molecules, etidocaine, and tetraphenylphosphonium to a phosphatidylcholine (POPC) lipid bilayer using the changes in lipid head-group order parameters. We observed that these molecules behave oppositely to calcium and sodium ions when binding to membranes: (i) their binding affinities are not overestimated by standard force field parameters, (ii) implicit inclusion of electronic polarizability increases their binding affinity, and (iii) they penetrate into the hydrophobic membrane core. Our results can be explained by distinct binding mechanisms of charged small molecules with hydrophobic moieties and monoatomic ions. The binding of the former is driven by hydrophobic effects, while the latter has direct electrostatic interactions with lipids. In addition to elucidating how different kinds of charged biomolecules bind to membranes, we deliver tools for further development of MD simulation parameters and methodology.
| S-EPMC9483918 | biostudies-literature