Accuracy and precision of NMR relaxation experiments and MD simulations for characterizing protein dynamics

Model-free parameters obtained from nuclear magnetic resonance (NMR) relaxation experiments and molecular dynamics (MD) simulations commonly are used to describe the intramolecular dynamical properties of proteins. To assess the relative accuracy and precision of experimental and simulated model-free parameters, three independent data sets derived from backbone ¹⁵N NMR relaxation experiments and two independent data sets derived from MD simulations of Escherichia coli ribonuclease HI are compared. The widths of the distributions of the differences between the order parameters for pairs of NMR data sets are congruent with the uncertainties derived from statistical analyses of individual data sets; thus, current protocols for analyzing NMR data encapsulate random uncertainties appropriately. Large differences in order parameters for certain residues are attributed to systematic differences between samples for intralaboratory comparisons and unknown, possibly magnetic field-dependent, experimental effects for interlaboratory comparisons. The widths of distributions of the differences between the order parameters for two NMR sets are similar to widths of distributions for an NMR and an MD set or for two MD sets. The linear correlations between the order parameters for an MD set and an NMR set are within the range of correlations observed between pairs of NMR sets. These comparisons suggest that the NMR and MD generalized order parameters for the backbone amide N—H bond vectors are of comparable accuracy for residues exhibiting motions on a fast time scale (<100 ps). Large discrepancies between NMR and MD order parameters for certain residues are attributed to the occurrence of “rare” motional events over the simulation trajectories, the disruption of an element of secondary structure in one of the simulations, and lack of consensus among the experimental data sets. Consequently, (easily detectable) severe distortions of local protein structure and infrequent motional events in MD simulations appear to be the most serious artifacts affecting the accuracy and precision, respectively, of MD order parameters relative to NMR values. In addition, MD order parameters for motions on a fast (<100 ps) timescale are more precisely determined than their NMR counterparts, thereby permitting more detailed dynamic characterization of biologically important residues by MD simulation than is sometimes possible by experimental methods. Proteins 28:481–493, 1997. © 1997 Wiley-Liss, Inc.     

An NMR spectroscopy study of bendaline-albumin interactions

The complete assignment of the 1H and 13C NMR spectra of bendaline (BNDL) was performed by mono-dimensional and homo- and hetero-correlated two-dimensional NMR experiments. The interaction between bendaline and albumin was also studied by the analysis of the motional parameters spin-lattice relaxation times, allowing the motional state of the BNDL free and bound with albumin to be defined. In absence of albumin the indazolacetic and benzylic moieties are characterized by roughly the same mobility and by positive sigma (cross-relaxation rates) values. In the presence of the macromolecule, the indazolacetic and benzylic moieties and the lysine change their motional behaviour to different extents, as indicated by correlation times. Data obtained in absence and in presence of the protein show that the molecular moiety of the bendaline most involved in the binding with albumin is the fragment H-4 H-5. The binding constant was evaluated at 2.4x10(3)M(-1).     

Conformational study on glycosylated asparagine-oligopeptides by NMR spectroscopy and molecular dynamics calculations.

The conformational properties of the homo oligomers of increasing chain length Boc-(Asn)(n)-NHMe (n = 2, 4, 5), (GlcNAc-beta-Asn)(n)-NHMe (n = 2, 4, 5, 8) and Boc-[GlcNAc(Ac)(3)-beta-Asn](n)-NHMe (n = 2, 4, 5) were studied by using NOE experiments and molecular dynamic calculations (MD). Sequential NOEs and medium range NOEs, including (i,i+2) interactions, were detected by ROESY experiments and quantified. The calculated inter-proton distances are longer than those characteristic of beta-turn secondary structures. Owing to the large conformational motions expected for linear peptides, MD simulations were performed without NMR constraints, with explicit water and by applying different treatments of the electrostatic interactions. In agreement with the NOE results, the simulations showed, for all peptides, the presence of both folded and unfolded structures. The existence of significant populations of beta-turn structures can be excluded for all the examined compounds, but two families of structures were more often recognized. The first one with sinusoidal or S-shaped forms, and another family of large turns together with some more extended conformations. Only the glycosylated pentapeptide shows in vacuo a large amount of structures with helical shaped form. The results achieved in water and in DMSO are compared and discussed, together with the effect of the glycosylation.     

Hydration study of homopolypeptides by (2)H NMR

Deuteron T(1) and T(2) was studied as a function of hydration in homopolyglycine (PG) and homopolyproline (PP). Water deuteron relaxation rates in PG conform to a hydration model involving two types of primary hydration sites where water is directly bonded to the polymer. Once these sites are filled, additional water only bonds to water molecules at the primary sites and in so doing affect their dynamics. PP exhibits an anomalous T(1) and T(2) hydration dependence which has been interpreted in terms of a cooperative water molecule-PP molecule helical conformational rearrangement which occurs once a certain hydration level is reached. The proposal of a water-PP structure is tested using molecular dynamics simulations.     

Specific DNA recognition by the Antp homeodomain: MD simulations of specific and nonspecific complexes

Four molecular dynamics simulation trajectories of complexes between the wild-type or a mutant Antennapedia homeodomain and 2 DNA sequences were generated in order to probe the mechanisms governing the specificity of DNA recognition. The starting point was published affinity measurements showing that a single protein mutation combined with a replacement of 2 base pairs yields a new high-affinity complex, whereas the other combinations, with changes on only 1 macromolecule, exhibited lower affinity. The simulations of the 4 complexes yielded fluctuating networks of interaction. On average, these networks differ significantly, explaining the switch of affinity caused by the alterations in the macromolecules. The network of mostly hydrogen-bonding interactions involving several water molecules, which was suggested both by X-ray and NMR structures of the wild-type homeodomain and its DNA operator sequence, could be reproduced in the trajectory. More interestingly, the high-affinity complex with alterations in both the protein and the DNA yielded again a dynamic but very tight network of intermolecular interactions, however, attributing a significantly stronger role to direct hydrophobic interactions at the expense of water bridges. The other 2 homeodomain-DNA complexes, with only 1 molecule altered, show on average over the trajectories a clearly reduced number of protein-DNA interactions. The observations from these simulations suggest specific experiments and thus close the circle formed by biochemical, structural, and computational studies. The shift from a water-dominated to a more "dry" interface may prove important in the design of proteins binding DNA in a specific manner.     

Conformational studies of neurohypophyseal hormones analogues with glycoconjugates by NMR spectroscopy

Two glycosylated peptides have been studied using NMR spectroscopy supported by molecular modeling. Peptide I is an oxytocin (OT) analogue in which glutamine 4 was replaced by serine with attached α‐d ‐mannose through the oxygen β atom, whereas peptide II is a lysine‐vasopressin (LVP) analogue with lysine 8 side chain modified by the attachment of glucuronic acid through an amide bond. Both peptides exhibit very weak uterotonic effect and are less susceptible to proteolytic degradation than the mother hormones. Additionally, peptide II reveals very weak pressor and antidiuretic activities. Our results have shown that the conformational preferences of glycosylated analogues are highly similar to those of their respective mother hormones. OT glycosylated analogue (I) exhibits a 3,4 β‐turn characteristic of OT‐like peptides, and vasopressin‐glycosylated analogue (II) exhibits β‐turns typical of vasopressin‐like peptides. Therefore, the lack of binding of the glycosylated analogues to the receptors can be attributed to a steric interference between the carbohydrate moieties and the receptors. We also consider this to be the reason of the very low activity of the analyzed glycopeptides. We expect that results from these studies will be helpful in designing new OT‐like and vasopressin‐like drugs. Copyright © 2014 European Peptide Society and John Wiley & Sons, Ltd.     

NMR reveals molecular interactions and dynamics of fatty acid binding to albumin

Background

The molecular details of fatty acid (FA) interactions with albumin are fundamental to understanding transport in the plasma and cellular utilization of these key nutrients and building blocks of membranes.

Scope of review

This review focuses on the development and application of NMR methods to study FA binding to albumin [bovine (BSA) and human (HSA)]. The key strategy was to use ¹³C enrichment of a specific carbon in the FA as a non-perturbing probe to permit visualization of the small ligand complexed to the very large protein. NMR contributions to illuminating molecular interactions and FA dynamics are summarized from three decades of studies.

Major conclusions

Our early studies detected multiple binding sites that we hypothesized were distinguished because of the unique tertiary structure of the protein in close proximity to the FA labeled carbon in each site. Later crystallographic structures revealed the presence of polar and charged amino acid side chains near the carboxyl carbon of the FA and unique tertiary structures lining all of the FA binding pockets. In collaboration with the crystallography group, several FA sites in the crystalline state were matched with NMR resonances in the solution state. With the newest application of NMR, 2D NMR spectroscopy detected nine binding sites, and three were located in the crystal structure through displacement of drugs with identified sites.

General significance

NMR spectroscopy utilizing the FA as a probe allows characterization of site-specific interactions, molecular motions within binding sites, the order of filling and removal of FA from sites. This article is part of a Special Issue entitled Serum Albumin.     

NMR studies on the conformation of oligomannosides and their interaction with banana lectin

The conformational and dynamic behaviour of three mannose containing oligosaccharides, a tetrasaccharide with α1→2, and α1→3, and a penta and a heptasaccharide with α1→2, α1→3, and α1→6 linkages has been evaluated by molecular mechanics and dynamics simulations and NMR spectroscopical methods. It is found that they display a fair amount of conformational freedom, with one major and one minor conformation per glycosidic linkage. The evaluation of their recognition by banana lectin has also been performed by STD NMR methods and a preliminary view of their putative interaction mode has been carried out by means of docking procedures.     

NMR structure of a concatemer of the first and second ligand-binding modules of the human low-density lipoprotein receptor

The ligand-binding domain of the human low-density lipoprotein receptor consists of seven modules, each of 40-45 residues. In the presence of calcium, these modules adopt a common polypeptide fold with three conserved disulfide bonds. A concatemer of the first and second modules (LB(1-2)) folds efficiently in the presence of calcium ions, forming the same disulfide connectivities as in the isolated modules. The three-dimensional structure of LB(1-2) has now been solved using two-dimensional 1H NMR spectroscopy and restrained molecular dynamics calculations. No intermodule nuclear Overhauser effects were observed, indicating the absence of persistent interaction between them. The near random-coil NH and H alpha chemical shifts and the low phi and psi angle order parameters of the four-residue linker suggest that it has considerable flexibility. The family of LB(1-2) structures superimposed well over LB1 or LB2, but not over both modules simultaneously. LB1 and LB2 have a similar pattern of calcium ligands, but the orientations of the indole rings of the tryptophan residues W23 and W66 differ, with the latter limiting solvent access to the calcium ion. From these studies, it appears that although most of the modules in the ligand-binding region of the receptor are joined by short segments, these linkers may impart considerable flexibility on this region.     

Lysozyme binding with sulfa group of antibiotics: Comparative binding thermodynamics and computational study

Protein–drug binding study addresses a broad domain of biological problems associating molecular functions to physiological processes composing and modifying safe and coherent drug therapeutics. Comparison of the binding and thermodynamic aspect of sulfa drugs, sulfamethazine (SMZ) and sulfadiazine (SDZ) with the protein, lysozyme (Lyz) was carried out using spectroscopic, molecular docking, and dynamic simulation studies. The fluorescence quenching and apparent binding constant for the binding reaction were calculated to be in the order of 10⁴ M⁻¹, slightly higher for SMZ as compared to that of SDZ and the binding stoichiometry values show 1:1 drug binding with each protein molecule. The binding was an enthalpy-driven spontaneous exothermic reaction favored by a negative enthalpy and a positive entropy contribution for both the complexes. The binding from the fluorescence quenching data suggests a static quenching mechanism dominated by non-polyelectrolytic components. Synchronous fluorescence denoted a conformational change in the tryptophan moiety of the protein. Molecular docking and dynamic simulation study provided a clearer view of the interaction pattern, where the drug resides on the binding pocket of the protein structure. Overall the protein, Lyz binding of SMZ was slightly more favored over SDZ.     

Molecular dynamics characterization of the active cavity of carboxypeptidase A and some of its inhibitor adducts.

Molecular dynamics (MD) calculations have been performed on carboxypeptidase A and on its adducts with inhibitors, such as d-phenylalanine (dPhe) and acetate. The catalytically essential zinc ion present in the protein was explicitly included in all the simulations. The simulation was carried out over a sphere of 15 A centered on the zinc ion. The crystallographic water molecules were explicitly taken into account; then the protein was solvated with a 18 A sphere of water molecules. MD calculations were carried out for 45-60 ps. There is no large deviation from the available X-ray structures of native and the dPhe adduct for the MD structures. Average MD structures were calculated starting from the X-ray structure of the dPhe adduct, and, from a structure obtained by docking the inhibitor in the native structure. Comparison between these two structures and with that of the native protein shows that some of the key variations produced by inhibitor binding are reproduced by MD calculations. Addition of acetate induces structural changes relevant for the understanding of the interaction network in the active cavity. The structural variations induced by different inhibitors are examined. The effects of these interactions on the catalytic mechanism and on the binding of substrate are discussed.     

NMR structure of Citrobacter freundii AmpD,comparison with bacteriophage T7 lysozyme and homology with PGRP domains

AmpD is a bacterial amidase involved in the recycling of cell-wall fragments in Gram-negative bacteria. Inactivation of AmpD leads to derepression of beta-lactamase expression, presenting a major pathway for the acquisition of constitutive antibiotic resistance. Here, we report the NMR structure of AmpD from Citrobacter freundii (PDB accession code 1J3G). A deep substrate-binding pocket explains the observed specificity for low molecular mass substrates. The fold is related to that of bacteriophage T7 lysozyme. Both proteins bind zinc at a conserved site and require zinc for amidase activity, although the enzymatic mechanism seems to differ in detail. The structure-based sequence alignment identifies conserved features that are also conserved in the eukaryotic peptidoglycan recognition protein (PGRP) domains, including the zinc-coordination site in several of them. PGRP domains thus belong to the same fold family and, where zinc-binding residues are conserved, may have amidase activity. This hypothesis is supported by the observation that human serum N-acetylmuramyl-L-alanine amidase seems to be identical with a soluble form of human PGRP-L.     

Molecular conformations of a disaccharide investigated using NMR spectroscopy

The molecular structure of -l-Rhap-(1→ 2)--l-Rhap-OMe has been investigated using conformation sensitive NMR parameters: cross-relaxation rates, scalar ³ J _CH couplings and residual dipolar couplings obtained in a dilute liquid crystalline phase. The order matrices of the two sugar residues are different, which indicates that the molecule cannot exist in a single conformation. The conformational distribution function, , related to the two glycosidic linkage torsion angles and was constructed using the APME method, valid in the low orientational order limit. The APME approach is based on the additive potential (AP) and maximum entropy (ME) models. The analyses of the trajectories generated in molecular dynamics and Langevin dynamics (LD) computer simulations gave support to the distribution functions constructed from the experimental NMR parameters. It is shown that at least two conformational regions are populated on the Ramachandran map and that these regions exhibit very different molecular order.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

NMR and MD studies on the interaction between ligand peptides and alpha-bungarotoxin

The interaction between alpha-bungarotoxin and linear synthetic peptides, mimotope of the nicotinic acetylcholine receptor binding site, has been characterised extensively by several methods and a wealth of functional, kinetic and structural data are available. Hence, this system represents a suitable model to explore in detail the dynamics of a peptide-protein interaction. Here, the solution structure of a new complex of the protein toxin with a tridecapeptide ligand exhibiting high affinity has been determined by NMR. As observed for three other previously reported mimotope-alpha-bungarotoxin complexes, also in this case correlations between biological activity and kinetic data are not fully consistent with a static discussion of structural data. Molecular dynamics simulations of the four mimotope-toxin complexes indicate that a relevant contribution to the complex stability is given by the extent of the residual flexibility that the protein maintains upon peptide binding. This feature, limiting the entropy loss caused by protein folding and binding, ought to be generally considered in a rational design of specific protein ligands.     

Phosphopeptide binding to the N-terminal SH2 domain of the p85 alpha subunit of PI 3'-kinase: a heteronuclear NMR study.

The N-terminal src-homology 2 domain of the p85 alpha subunit of phosphatidylinositol 3' kinase (SH2-N) binds specifically to phosphotyrosine-containing sequences. Notably, it recognizes phosphorylated Tyr 751 within the kinase insert of the cytoplasmic domain of the activated beta PDGF receptor. A titration of a synthetic 12-residue phosphopeptide (ESVDY*VPMLDMK) into a solution of the SH2-N domain was monitored using heteronuclear 2D and 3D NMR spectroscopy. 2D-(15N-1H) heteronuclear single-quantum correlation (HSQC) experiments were performed at each point of the titration to follow changes in both 15N and 1H chemical shifts in NH groups. When mapped onto the solution structure of the SH2-N domain, these changes indicate a peptide-binding surface on the protein. Line shape analysis of 1D profiles of individual (15N-1H)-HSQC peaks at each point of the titration suggests a kinetic exchange model involving at least 2 steps. To characterize changes in the internal dynamics of the domain, the magnitude of the (15N-1H) heteronuclear NOE for the backbone amide of each residue was determined for the SH2-N domain with and without bound peptide. These data indicate that, on a nanosecond timescale, there is no significant change in the mobility of either loops or regions of secondary structure. A mode of peptide binding that involves little conformational change except in the residues directly involved in the 2 binding pockets of the p85 alpha SH2-N domain is suggested by this study.     

Domain motions in bacteriophage T4 lysozyme: A comparison between molecular dynamics and crystallographic data

A comparison of a series of extended molecular dynamics (MD) simulations of bacteriophage T4 lysozyme in solvent with X-ray data is presented. Essential dynamics analyses were used to derive collective fluctuations from both the simulated trajectories and a distribution of crystallographic conformations. In both cases the main collective fluctuations describe domain motions. The protein consists of an N- and C-terminal domain connected by a long helix. The analysis of the distribution of crystallographic conformations reveals that the N-terminal helix rotates together with either of these two domains. The main domain fluctuation describes a closure mode of the two domains in which the N-terminal helix rotates concertedly with the C-terminal domain, while the domain fluctuation with second largest amplitude corresponds to a twisting mode of the two domains, with the N-terminal helix rotating concertedly with the N-terminal domain. For the closure mode, the difference in hinge-bending angle between the most open and most closed X-ray structure along this mode is 49 degrees. In the MD simulation that shows the largest fluctuation along this mode, a rotation of 45 degrees was observed. Although the twisting mode has much less freedom than the closure mode in the distribution of crystallographic conformations, experimental results suggest that it might be functionally important. Interestingly, the twisting mode is sampled more extensively in all MD simulations than it is in the distribution of X-ray conformations. Proteins 31:116–127, 1998. © 1998 Wiley-Liss, Inc.     

Multiple binding modes for palmitate to barley lipid transfer protein facilitated by the presence of proline 12

Molecular dynamics simulations have been used to characterise the binding of the fatty acid ligand palmitate in the barley lipid transfer protein 1 (LTP) internal cavity. Two different palmitate binding modes (1 and 2), with similar protein–ligand interaction energies, have been identified using a variety of simulation strategies. These strategies include applying experimental protein–ligand atom–atom distance restraints during the simulation, or protonating the palmitate ligand, or using the vacuum GROMOS 54B7 force‐field parameter set for the ligand during the initial stages of the simulations. In both the binding modes identified the palmitate carboxylate head group hydrogen bonds with main chain amide groups in helix A, residues 4 to 19, of the protein. In binding mode 1 the hydrogen bonds are to Lys 11, Cys 13, and Leu 14 and in binding mode 2 to Thr 15, Tyr 16, Val 17, Ser 24 and also to the OH of Thr 15. In both cases palmitate binding exploits irregularity of the intrahelical hydrogen‐bonding pattern in helix A of barley LTP due to the presence of Pro 12. Simulations of two variants of barley LTP, namely the single mutant Pro12Val and the double mutant Pro12Val Pro70Val, show that Pro 12 is required for persistent palmitate binding in the LTP cavity. Overall, the work identifies key MD simulation approaches for characterizing the details of protein–ligand interactions in complexes where NMR data provide insufficient restraints.     

Probing the interaction between cHAVc3 peptide and the EC1 domain of E-cadherin using NMR and molecular dynamics simulations

The goal of this work is to probe the interaction between cyclic cHAVc3 peptide and the EC1 domain of human E-cadherin protein. Cyclic cHAVc3 peptide (cyclo(1,6)Ac-CSHAVC-NH₂) binds to the EC1 domain as shown by chemical shift perturbations in the 2D ¹H,-¹⁵N-HSQC NMR spectrum. The molecular dynamics (MD) simulations of the EC1 domain showed folding of the C-terminal tail region into the main head region of the EC1 domain. For cHAVc3 peptide, replica exchange molecular dynamics (REMD) simulations generated five structural clusters of cHAVc3 peptide. Representative structures of cHAVc3 and the EC1 structure from MD simulations were used in molecular docking experiments with NMR constraints to determine the binding site of the peptide on EC1. The results suggest that cHAVc3 binds to EC1 around residues Y36, S37, I38, I53, F77, S78, H79, and I94. The dissociation constants (K_d values) of cHAVc3 peptide to EC1 were estimated using the NMR chemical shifts data and the estimated K_ds are in the range of .5 × 10^?5–7.0 × 10^?5 M.     

MD and NMR studies of alpha-bungarotoxin surface accessibility

Protein surface accessibility represents a dimension of structural biology which has not been discussed in details so far, in spite of its fundamental role in controlling the molecular recognition process. In the present report the surface accessibility of alpha-bungarotoxin, a small and well characterized protein, has been investigated by analyzing its interaction with solvent and paramagnetic molecules in an integrated way. The presence of strong hydration sites, identified by a combined analysis of MD simulation and NMR results, seems to prevent the access of Gd(III)DTPA-BMA to the protein surface. On the contrary, the limited hydration of the alpha-bungarotoxin active site favors frequent encounters between the paramagnetic probe and the protein in the latter region. All the data obtained here for alpha-bungarotoxin suggest that shape and stability of the solvation shell control its surface accessibility and, hence, intermolecular interactions in a way which could be common to many other proteins.     

Characterization of the structure and dynamics of amyloidogenic variants of human lysozyme by NMR spectroscopy

The structures and dynamics of the native states of two mutational variants of human lysozyme, I56T and D67H, both associated with non-neuropathic systemic amyloidosis, have been investigated by NMR spectroscopy. The (1)H and (15)N main-chain amide chemical shifts of the I56T variant are very similar to those of the wild-type protein, but those of the D67H variant are greatly altered for 28 residues in the beta-domain. This finding is consistent with the X-ray crystallographic analysis, which shows that the structure of this variant is significantly altered from that of the wild-type protein in this region. The (1)H-(15)N heteronuclear NOE values show that, with the exception of V121, every residue in the wild-type and I56T proteins is located in tightly packed structures characteristic of the native states of most proteins. In contrast, D67H has a region of substantially increased mobility as shown by a dramatic decrease in heteronuclear NOE values of residues near the site of mutation. Despite this unusual flexibility, the D67H variant has no greater propensity to form amyloid fibrils in vivo or in vitro than has I56T. This finding indicates that it is the increased ability of the variants to access partially folded conformations, rather than intrinsic changes in their native state properties, that is the origin of their amyloidogenicity.