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.

Dynamic Water-Mediated Hydrogen Bonding in a Collagen Model Peptide

In the canonical (G-X-Y)(n) sequence of the fibrillar collagen triple helix, stabilizing direct interchain hydrogen bonding connects neighboring chains. Mutations of G can disrupt these interactions and are linked to connective tissue diseases. Here we integrate computational approaches with nuclear magnetic resonance (NMR) to obtain a dynamic view of hydrogen bonding distributions in the (POG)(4)(-)(POA)-(POG)(5) peptide, showing that the solution conformation, dynamics, and hydrogen bonding deviate from the reported X-ray crystal structure in many aspects. The simulations and NMR data provide clear evidence of inequivalent environments in the three chains. Molecular dynamics (MD) simulations indicate direct interchain hydrogen bonds in the leading chain, water bridges in the middle chain, and nonbridging waters in the trailing chain at the G → A substitution site. Theoretical calculations of NMR chemical shifts using a quantum fragmentation procedure can account for the unusual downfield NMR chemical shifts at the substitution sites and are used to assign the resonances to the individual chains. The NMR and MD data highlight the sensitivity of amide shifts to changes in the acceptor group from peptide carbonyls to water. The results are used to interpret solution NMR data for a variety of glycine substitutions and other sequence triplet interruptions to provide new connections between collagen sequences, their associated structures, dynamical behavior, and their ability to recognize collagen receptors.

Bayesian inference of conformational state populations from computational models and sparse experimental observables

We present a Bayesian inference approach to estimating conformational state populations from a combination of molecular modeling and sparse experimental data. Unlike alternative approaches, our method is designed for use with small molecules and emphasizes high-resolution structural models, using inferential structure determination with reference potentials, and Markov Chain Monte Carlo to sample the posterior distribution of conformational states. As an application of the method, we determine solution-state conformational populations of the 14-membered macrocycle cineromycin B, using a combination of previously published sparse Nuclear Magnetic Resonance (NMR) observables and replica-exchange molecular dynamic/Quantum Mechanical (QM)-refined conformational ensembles. Our results agree better with experimental data compared to previous modeling efforts. Bayes factors are calculated to quantify the consistency of computational modeling with experiment, and the relative importance of reference potentials and other model parameters.

Motional timescale predictions by molecular dynamics simulations: case study using proline and hydroxyproline sidechain dynamics

We propose a new approach for force field optimizations which aims at reproducing dynamics characteristics using biomolecular MD simulations, in addition to improved prediction of motionally averaged structural properties available from experiment. As the source of experimental data for dynamics fittings, we use ¹³C NMR spin-lattice relaxation times T₁ of backbone and sidechain carbons, which allow to determine correlation times of both overall molecular and intramolecular motions. For structural fittings, we use motionally averaged experimental values of NMR J couplings. The proline residue and its derivative 4-hydroxyproline with relatively simple cyclic structure and sidechain dynamics were chosen for the assessment of the new approach in this work. Initially, grid search and simplexed MD simulations identified large number of parameter sets which fit equally well experimental J couplings. Using the Arrhenius-type relationship between the force constant and the correlation time, the available MD data for a series of parameter sets were analyzed to predict the value of the force constant that best reproduces experimental timescale of the sidechain dynamics. Verification of the new force-field (termed as AMBER99SB-ILDNP) against NMR J couplings and correlation times showed consistent and significant improvements compared to the original force field in reproducing both structural and dynamics properties. The results suggest that matching experimental timescales of motions together with motionally averaged characteristics is the valid approach for force field parameter optimization. Such a comprehensive approach is not restricted to cyclic residues and can be extended to other amino acid residues, as well as to the backbone. Proteins 2014; 82:195–215. © 2013 Wiley Periodicals, Inc.

Contribution of high-energy conformations to NMR chemical shifts, a DFT-BOMD study

This paper highlights the relevance of including the high-energy conformational states sampled by Born-Oppenheimer molecular dynamics (BOMD) in the calculation of time-averaged NMR chemical shifts. Our case study is the very flexible glycerol molecule that undergoes interconversion between conformers in a nonrandom way. Along the sequence of structures from one backbone conformer to another, transition states have been identified. The three (13)C NMR chemical shifts of the molecule were estimated by averaging their calculated values over a large set of BOMD snapshots. The simulation time needed to obtain a good agreement with the two signals present in the experimental spectrum is shown to be dependent on the atomic orbital basis set used for the dynamics, with a necessary longer trajectory for the most extended basis sets. The large structural deformations with respect to the optimized conformer geometries that occur along the dynamics are related to a kinetically driven conformer distribution. Calculated conformer type populations are in good agreement with experimental gas phase microwave results.

Comprehensive multiphase NMR spectroscopy: basic experimental approaches to differentiate phases in heterogeneous samples

Heterogeneous samples, such as soils, sediments, plants, tissues, foods and organisms, often contain liquid-, gel- and solid-like phases and it is the synergism between these phases that determine their environmental and biological properties. Studying each phase separately can perturb the sample, removing important structural information such as chemical interactions at the gel-solid interface, kinetics across boundaries and conformation in the natural state. In order to overcome these limitations a Comprehensive Multiphase-Nuclear Magnetic Resonance (CMP-NMR) probe has been developed, and is introduced here, that permits all bonds in all phases to be studied and differentiated in whole unaltered natural samples. The CMP-NMR probe is built with high power circuitry, Magic Angle Spinning (MAS), is fitted with a lock channel, pulse field gradients, and is fully susceptibility matched. Consequently, this novel NMR probe has to cover all HR-MAS aspects without compromising power handling to permit the full range of solution-, gel- and solid-state experiments available today. Using this technology, both structures and interactions can be studied independently in each phase as well as transfer/interactions between phases within a heterogeneous sample. This paper outlines some basic experimental approaches using a model heterogeneous multiphase sample containing liquid-, gel- and solid-like components in water, yielding separate (1)H and (13)C spectra for the different phases. In addition, (19)F performance is also addressed. To illustrate the capability of (19)F NMR soil samples, containing two different contaminants, are used, demonstrating a preliminary, but real-world application of this technology. This novel NMR approach possesses a great potential for the in situ study of natural samples in their native state.

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.