Backbone dynamics in collagen

Water scaffolding in collagen: Implications on protein dynamics as revealed by solid-state NMR

Solid-state NMR studies of collagen samples of various origins confirm that the amplitude of collagen backbone and sidechain motions increases significantly on increasing the water content. This conclusion is supported by the changes observed in three different NMR observables: (i) the linewidth dependence on the 1H decoupling frequency; (ii) 13C CSA changes for the peptide carbonyl groups, and (iii) dephasing rates of 1H-13C dipolar couplings. In particular, a nearly threefold increase in motional amplitudes of the backbone librations about C-Cα or N-Cα bonds was found on increasing the added water content up to 47 wt%D2 O. On the basis of the frequencies of NMR observables involved, the timescale of the protein motions dependent on the added water content is estimated to be of the order of microseconds. This estimate agrees with that from wideline T2(1)H NMR measurements. Also, our wideline 1H NMR measurements revealed that the timescale of the microsecond motions in proteins reduces significantly on increasing the added water content, i.e., an ∼15-fold increase in protein motional frequencies is observed on increasing the added water content to 45 wt% D2 O. The observed changes in collagen dynamics is attributed to the increase in water translational diffusion on increasing the amount of added water, which leads to more frequent "bound water/free water" exchange on the protein surface, accompanied by the breakage and formation of new hydrogen bonds with polar functionalities of protein.

Vanishing amplitude of backbone dynamics causes a true protein dynamical transition: 2H NMR studies on perdeuterated C-phycocyanin

Using a combination of H2 nuclear magnetic resonance (NMR) methods, we study internal rotational dynamics of the perdeuterated protein C-phycocyanin (CPC) in dry and hydrated states over broad temperature and dynamic ranges with high angular resolution. Separating H2 NMR signals from methyl deuterons, we show that basically all backbone deuterons exhibit highly restricted motion occurring on time scales faster than microseconds. The amplitude of this motion increases when a hydration shell exists, while it decreases upon cooling and vanishes near 175 K. We conclude that the vanishing of the highly restricted motion marks a dynamical transition, which is independent of the time window and of a fundamental importance. This conclusion is supported by results from experimental and computational studies of the proteins myoglobin and elastin. In particular, we argue based on findings in molecular dynamics simulations that the behavior of the highly restricted motion of proteins at the dynamical transition resembles that of a characteristic secondary relaxation of liquids at the glass transition, namely the nearly constant loss. Furthermore, H2 NMR studies on perdeuterated CPC reveal that, in addition to highly restricted motion, small fractions of backbone segments exhibit weakly restricted dynamics when temperature and hydration are sufficiently high.

Concise NMR approach for molecular dynamics characterizations in organic solids

Molecular dynamics characterisations in solids can be carried out selectively using dipolar-dephasing experiments. Here we show that the introduction of a sum of Lorentzian and Gaussian functions greatly improve fittings of the "intensity versus time" data for protonated carbons in dipolar-dephasing experiments. The Lorentzian term accounts for remote intra- and intermolecular (1)H-(13)C dipole-dipole interactions, which vary from one molecule to another or for different carbons within the same molecule. Thus, by separating contributions from weak remote interactions, more accurate Gaussian decay constants, T(dd), can be extracted for directly bonded (1)H-(13)C dipole-dipole interactions. Reorientations of the (1)H-(13)C bonds lead to the increase of T(dd), and by measuring dipolar-dephasing constants, insight can be gained into dynamics in solids. We have demonstrated advantages of the method using comparative dynamics studies in the α and γ polymorphs of glycine, cyclic amino acids L-proline, DL-proline and trans-4-hydroxy-L-proline, the Ala residue in different dipeptides, as well as adamantane and hexamethylenetetramine. It was possible to distinguish subtle differences in dynamics of different carbon sites within a molecule in polymorphs and in L- and DL-forms. The presence of overall molecular motions is shown to lead to particularly large differences in dipolar-dephasing experiments. The differences in dynamics can be attributed to differences in noncovalent interactions. In the case of hexamethylenetetramine, for example, the presence of C-H···N interactions leads to nearly rigid molecules. Overall, the method allows one to gain insight into the role of noncovalent interactions in solids and their influence on the molecular dynamics.

Experimental verification of force fields for molecular dynamics simulations using Gly-Pro-Gly-Gly

Experimental NMR verification of MD simulations using 12 different force fields (AMBER, CHARMM, GROMOS, and OPLS-AA) and 5 different water models has been undertaken to identify reliable MD protocols for structure and dynamics elucidations of small open chain peptides containing Gly and Pro. A conformationally flexible tetrapeptide Gly-Pro-Gly-Gly was selected for NMR (3)J-coupling, chemical shift, and internuclear distance measurements, followed by their calculations using 2 μs long MD simulations in water. In addition, Ramachandran population maps for Pro-2 and Gly-3 residues of GPGG obtained from MD simulations were used for detailed comparisons with similar maps from the protein data bank (PDB) for large number of Gly and Pro residues in proteins. The MD simulations revealed strong dependence of the populations and geometries of preferred backbone and side chain conformations, as well as the time scales of the peptide torsional transitions on the force field used. On the basis of the analysis of the measured and calculated data, AMBER99SB is identified as the most reliable force field for reproducing NMR measured parameters, which are dependent on the peptide backbone and the Pro side chain geometries and dynamics. Ramachandran maps showing the dependence of conformational populations as a function of backbone ϕ/ψ angles for Pro-2 and Gly-3 residues of GPGG from MD simulations using AMBER99SB, AMBER03, and CHARMM were found to resemble similar maps for Gly and Pro residues from the PDB survey. Three force fields (AMBER99, AMBER99ϕ, and AMBER94) showed the least satisfactory agreement with both the solution NMR and the PDB survey data. The poor performance of these force fields is attributed to their propensity to overstabilize helical peptide backbone conformations at the Pro-2 and Gly-3 residues. On the basis of the similarity of the MD and PDB Ramachandran plots, the following sequence of transitions is suggested for the Gly backbone conformation: α(L) ⇆ β(PR) ⇆ β(S) ⇆ β(P) ⇆ α, where backbone secondary structures α(L) and α are associated with helices and turns, β(P) and β(PR) correspond to the left- and right-handed polyproline II structures and β(S) denotes the fully stretched backbone conformation. Compared to the force field dependence, less significant, but noteworthy, variations in the populations of the peptide backbone conformations were observed. For different solvent models considered, a correlation was noted between the number of torsional transitions in GPGG and the water self-diffusion coefficient on using TIP3P, TIP4P, and TIP5P models. In addition to MD results, we also report DFT derived Karplus relationships for Gly and Pro residues using B972 and B3LYP functionals.

The effects of anticalcification treatments and hydration on the molecular dynamics of bovine pericardium collagen as revealed by 13C solid-state NMR

This article describes a solid-state NMR (SSNMR) investigation of the influence of hydration and chemical cross-linking on the molecular dynamics of the constituents of the bovine pericardium (BP) tissues and its relation to the mechanical properties of the tissue. Samples of natural phenethylamine-diepoxide (DE)- and glutaraldehyde (GL)-fixed BP were investigated by (13)C cross-polarization SSNMR to probe the dynamics of the collagen, and the results were correlated to the mechanical properties of the tissues, probed by dynamical mechanical analysis. For samples of natural BP, the NMR results show that the higher the hydration level the more pronounced the molecular dynamics of the collagen backbone and sidechains, decreasing the tissue's elastic modulus. In contrast, in DE- and GL-treated samples, the collagen molecules are more rigid, and the hydration seems to be less effective in increasing the collagen molecular dynamics and reducing the mechanical strength of the samples. This is mostly attributed to the presence of cross-links between the collagen plates, which renders the collagen mobility less dependent on the water absorption in chemically treated samples.

Structural and Dynamic Aspects of Hydrogen-Bonded Complexes and Inclusion Compounds Containing α,ω-Dicyanoalkanes and Urea, Investigated by Solid-State 13C and 2H NMR Techniques

Solid-state 13C NMR and 2H NMR techniques have been used to investigate structural and dynamic properties of the 1,4-dicyanobutane/urea and 1,5-dicyanopentane/urea 1:1 hydrogen-bonded complexes and the 1,6-dicyanohexane/urea inclusion compound. The pure crystalline phase of urea has also been investigated. The 13C NMR studies have focused on 13C chemical shift anisotropy and second-order quadrupolar effects (arising from 13C-14N interaction) for the urea molecules and the cyano groups of the alpha,omega-dicyanoalkanes. Parameters describing these interactions are derived and are discussed in relation to the known structural properties of these materials. Comparison of 13C chemical shift anisotropies of the cyano carbons and rates of 13C dipolar dephasing suggest that 1,4-dicyanobutane and 1,5-dicyanopentane are effectively static, whereas 1,6-dicyanohexane has greater mobility. 2H NMR line shape analysis for the 1,4-dicyanobutane/urea-d4 and 1,5-dicyanopentane/urea-d4 complexes indicates that the only motion of the urea molecules that is effective on the 2H NMR time scale is a rapid libration about the C=O bond over an angular range of about 26 degrees . For the 1,6-dicyanohexane/urea-d4 inclusion compound, the 2H NMR line shape is consistent with a motion comprising 180 degrees jumps about the C=O bond at rates that are intermediate on the 2H NMR time scale. In addition, rapid libration about the C=O bond also occurs over an angular range of about 20 degrees . The dynamic properties of the urea molecules in these materials are compared with those of urea molecules in other crystalline environments.