Correlation of protein structure and dynamics to scalar couplings across hydrogen bonds

NMR-observable scalar couplings across hydrogen bonds in nucleic acids and proteins present a quantitative measure for the geometry and--by the implicit experimental time averaging--dynamics of hydrogen bonds. We have carried out in-depth molecular dynamics (MD) simulations with various force fields on three proteins: ubiquitin, the GB1 domain of protein G, and the SMN Tudor domain, for which experimental h3JNC' scalar couplings of backbone hydrogen bonds and various high-resolution X-ray structures are available. Theoretical average values for h3JNC' were calculated from the snapshots of these MD simulations either by density functional theory or by a geometric parametrization (Barfield, M. J. Am. Chem. Soc. 2002, 124, 4158-4168). No significant difference was found between the two methods. The results indicate that time-averaging using explicit water solvation in the MD simulations improves significantly the agreement between experimental and theoretical values for the lower resolution structures ubiquitin (1.8 A), Tudor domain (1.8 A), and protein G (2.1 A). Only marginal improvement is found for the high-resolution structure (1.1 A) of protein G. Hence, experimental h3JNC' values are compatible with a static, high-resolution structural model. The MD averaging of the low-resolution structures moves the averages of the rHO distance and the H...O=C angle theta closer to their respective values in the high-resolution structures, thereby improving the agreement using experimental h3JNC' data. In contrast, MD averaging with implicit water models deteriorates the agreement with experiment for all proteins. The differing behavior can be explained by an artifactual lengthening of H-bonds caused by the implicit water models.     

Long-range magnetization transfer between uncoupled nuclei by dipole-dipole cross-correlated relaxation: a precise probe of beta-sheet geometry in proteins

Interference between dipolar interactions in covalently linked (13)C-(1)H and nonlinked (1)H-(1)H pairs can be used to generate antiphase magnetization between noncoupled spins. The buildup rate of such antiphase terms is highly sensitive to local geometry, in particular the interproton distance and the (13)C-(1)H-(1)H internuclear angle. These rates have been measured for opposing C(alpha)H(alpha) pairs in antiparallel beta-sheets in the third Igg-binding domain of protein G (GB3) and in HIV protease, complexed with the inhibitor DMP323. For GB3, good agreement with the 1.1-A crystal structure is found. However, this agreement rapidly deteriorates with decreasing resolution of the corresponding X-ray structure. For HIV protease, two separate crystal structures that differ by less than 0.2 A from one another exhibit lower agreement in their predicted cross-correlated relaxation rates relative to one another than is found between experimental rates and the average of the rates predicted for the two structures. These data indicate that quantitative measurement of these cross-correlated relaxation rates can provide highly accurate structural information in macromolecules.     

Protein structure determination by high-resolution solid-state NMR spectroscopy: application to microcrystalline ubiquitin

High-resolution solid-state NMR spectroscopy has become a promising method for the determination of three-dimensional protein structures for systems which are difficult to crystallize or exhibit low solubility. Here we describe the structure determination of microcrystalline ubiquitin using 2D (13)C-(13)C correlation spectroscopy under magic angle spinning conditions. High-resolution (13)C spectra have been acquired from hydrated microcrystals of site-directed (13)C-enriched ubiquitin. Inter-residue carbon-carbon distance constraints defining the global protein structure have been evaluated from 'dipolar-assisted rotational resonance' experiments recorded at various mixing times. Additional constraints on the backbone torsion angles have been derived from chemical shift analysis. Using both distance and dihedral angle constraints, the structure of microcrystalline ubiquitin has been refined to a root-mean-square deviation of about 1 A. The structure determination strategies for solid samples described herein are likely to be generally applicable to many proteins that cannot be studied by X-ray crystallography or solution NMR spectroscopy.     

Crystal structure of a genomically encoded fosfomycin resistance protein (FosA) at 1.19 A resolution by MAD phasing off the L-III edge of Tl(+)

The fosfomycin resistance protein (FosA) catalyzes the Mn(II)- and K+-dependent addition of glutathione to the oxirane of the antibiotic fosfomycin. The crystal structure of FosA from Pseudomonas aeruginosa was solved at a resolution of 1.19 A by multiwavelength anomalous diffraction at the L-III edge of a Tl+ derivative. The structure solution took advantage of the ability of Tl+ to substitute for K+. The existence of multiple Tl sites in the asymmetric unit suggests that this may be a generally useful technique for phasing protein crystal structures. A 1.35 A resolution structure with phosphate bound in the active site shows that the Mn(II) center has a rare four-coordinate geometry. The structure of the fosfomycin complex at 1.19 A resolution indicates that the Mn(II) center is close to five-coordinate with trigonal bipyramidal geometry and a ligand set consisting of two histidines (H7 and H64) and one phosphonate oxygen occupying the equatorial sites and the carboxylate of E110 at one of the apical sites. The oxirane oxygen of the substrate is located at the other apical site but is 0.2 A beyond the average Mn-O distance for five-coordinate Mn(II). The Mn(II) center is proposed to stabilize the alkoxide in the transition state, while the nearby hydroxyl group of T9 acts as a proton donor in the reaction. The K+ ion located 6.5 A from the Mn(II) appears to help orient the substrate for nucleophilic attack.     

Evaluation of backbone proton positions and dynamics in a small protein by liquid crystal NMR spectroscopy

NMR measurements of a large set of protein backbone one-bond dipolar couplings have been carried out to refine the structure of the third IgG-binding domain of Protein G (GB3), previously solved by X-ray crystallography at a resolution of 1.1 A. Besides the commonly used bicelle, poly(ethylene glycol), and filamentous phage liquid crystalline media, dipolar couplings were also measured when the protein was aligned inside either positively or negatively charged stretched acrylamide gels. Refinement of the GB3 crystal structure against the (13)C(alpha)-(13)C' and (13)C'-(15)N dipolar couplings improves the agreement between experimental and predicted (15)N-(1)H(N) as well as (13)C(alpha)-(1)H(alpha) dipolar couplings. Evaluation of the peptide bond N-H orientations shows a weak anticorrelation between the deviation of the peptide bond torsion angle omega from 180 degrees and the angle between the N-H vector and the C'-N-C(alpha) plane. The slope of this correlation is -1, indicating that, on average, pyramidalization of the peptide N contributes to small deviations from peptide bond planarity ( = 179.3 +/- 3.1 degrees ) to the same degree as true twisting around the C'-N bond. Although hydrogens are commonly built onto crystal structures assuming the N-H vector orientation falls on the line bisecting the C'-N-C(alpha) angle, a better approximation adjusts the C(alpha)-C'-N-H torsion angle to -2 degrees. The (15)N-(1)H(N) dipolar data do not contradict the commonly accepted motional model where angular fluctuations of the N-H bond orthogonal to the peptide plane are larger than in-plane motions, but the amplitude of angular fluctuations orthogonal the C(alpha)(i-1)-N(i)-C(alpha)(i) plane exceeds that of in-plane motions by at most 10-15 degrees. Dipolar coupling analysis indicates that for most of the GB3 backbone, the amide order parameters, S, are highly homogeneous and vary by less than +/-7%. Evaluation of the H(alpha) proton positions indicates that the average C(alpha)-H(alpha) vector orientation deviates by less than 1 degrees from the direction that makes ideal tetrahedral angles with the C(alpha)-C(beta) and C(alpha)-N vectors.     

The first solution structure of a paramagnetic copper(II) protein: the case of oxidized plastocyanin from the cyanobacterium Synechocystis PCC6803.

The NMR solution structure of oxidized plastocyanin from the cyanobacterium Synechocystis PCC6803 is here reported. The protein contains paramagnetic copper(II), whose electronic relaxation times are quite unfavorable for NMR solution studies. The structure has been solved on the basis of 1041 meaningful NOESY cross-peaks, 18 1D NOEs, 26 T(1) values, 96 dihedral angle constraints, and 18 H-bonds. The detection of broad hyperfine-shifted signals and their full assignment allowed the identification of the copper(II) ligands and the determination of the Cu-S-C-H dihedral angle for the coordinated cysteine. The global root-mean-square deviation from the mean structure for the solution structure family is 0.72 +/- 0.14 and 1.16 +/- 0.17 A for backbone and heavy atoms, respectively. The structure is overall quite satisfactory and represents a breakthrough, in that it includes paramagnetic copper proteins among the metalloproteins for which solution structures can be afforded. The comparison with the available X-ray structure of a triple mutant is also performed.     

X-ray absorption spectroscopy of hemes and hemeproteins in solution: multiple scattering analysis

A full quantitative analysis of Fe K-edge X-ray absorption spectra has been performed for hemes in two porphynato complexes, that is, iron(III) tetraphenylporphyrin chloride (Fe(III)TPPCl) and iron(III) tetraphenylporphyrin bis(imidazole) (Fe(III)TPP(Imid)2), in two protein complexes whose X-ray structure is known at atomic resolution (1.0 A), that is, ferrous deoxy-myoglobin (Fe(II)Mb) and ferric aquo-myoglobin (Fe(III)MbH2O), and in ferric cyano-myoglobin (Fe(III)MbCN), whose X-ray structure is known at lower resolution (1.4 A). The analysis has been performed via the multiple scattering approach, starting from a muffin tin approximation of the molecular potential. The Fe-heme structure has been obtained by analyzing independently the Extended X-ray Absorption Fine Structure (EXAFS) region and the X-ray Absorption Near Edge Structure (XANES) region. The EXAFS structural results are in full agreement with the crystallographic values of the models, with an accuracy of +/- 0.02 A for Fe-ligand distances, and +/-6 degrees for angular parameters. All the XANES features above the theoretical zero energy (in the lower rising edge) are well accounted for by single-channel calculations, for both Fe(II) and Fe(III) hemes, and the Fe-N p distance is determined with the same accuracy as EXAFS. XANES evaluations of Fe-5th and Fe-6th ligand distances are determined with 0.04-0.07 A accuracy; a small discrepancy with EXAFS (0.01 to 0.05 A beyond the statistical error), is found for protein compounds. Concerns from statistical correlation among parameters and multiple minima in the parameter space are discussed. As expected, the XANES accuracy is slightly lower than what was found for polarized XANES on Fe(III)MbCN single crystal (0.03-0.04 A), and states the actual state-of-the-art of XANES analysis when used to extract heme-normal parameters in a solution spectrum dominated by heme-plane scattering.     

Automated protein structure determination from NMR spectra

Fully automated structure determination of proteins in solution (FLYA) yields, without human intervention, three-dimensional protein structures starting from a set of multidimensional NMR spectra. Integrating existing and new software, automated peak picking over all spectra is followed by peak list filtering, the generation of an ensemble of initial chemical shift assignments, the determination of consensus chemical shift assignments for all (1)H, (13)C, and (15)N nuclei, the assignment of NOESY cross-peaks, the generation of distance restraints, and the calculation of the three-dimensional structure by torsion angle dynamics. The resulting, preliminary structure serves as additional input to the second stage of the procedure, in which a new ensemble of chemical shift assignments and a refined structure are calculated. The three-dimensional structures of three 12-16 kDa proteins computed with the FLYA algorithm coincided closely with the conventionally determined structures. Deviations were below 0.95 A for the backbone atom positions, excluding the flexible chain termini. 96-97% of all backbone and side-chain chemical shifts in the structured regions were assigned to the correct residues. The purely computational FLYA method is suitable for substituting all manual spectra analysis and thus overcomes a main efficiency limitation of the NMR method for protein structure determination.     

Close contacts between carbonyl oxygen atoms and aromatic centers in protein structures: pi...pi or lone-pair...pi interactions?

Lone-pair...pi and, more recently, pi...pi interactions have been studied in small molecule crystal structures, and they are the focus of attention in some biomolecules. In this study, we have systematically analyzed 500 high-resolution protein structures (resolution < or =1.8 A) and identified 286 examples in which carbonyl oxygen atoms approach the aromatic centers within a distance of 3.5 A. Contacts involving backbone carbonyl oxygens are frequently observed in helices and, to some extent, in strands. Geometrical characterization indicates that these contacts have geometry in between that of an ideal pi...pi and a lone-pair...pi interaction. Quantum mechanical calculations using 6-311++G** basis sets reveal that these contacts give rise to energetically favorable interactions and, along with MD simulations, indicate that such interactions could stabilize secondary structures.     

GRID: A high‐resolution protein structure refinement algorithm

The energy‐based refinement of protein structures generated by fold prediction algorithms to atomic‐level accuracy remains a major challenge in structural biology. Energy‐based refinement is mainly dependent on two components: (1) sufficiently accurate force fields, and (2) efficient conformational space search algorithms. Focusing on the latter, we developed a high‐resolution refinement algorithm called GRID. It takes a three‐dimensional protein structure as input and, using an all‐atom force field, attempts to improve the energy of the structure by systematically perturbing backbone dihedrals and side‐chain rotamer conformations. We compare GRID to Backrub, a stochastic algorithm that has been shown to predict a significant fraction of the conformational changes that occur with point mutations. We applied GRID and Backrub to 10 high‐resolution (≤ 2.8 Å) crystal structures from the Protein Data Bank and measured the energy improvements obtained and the computation times required to achieve them. GRID resulted in energy improvements that were significantly better than those attained by Backrub while expending about the same amount of computational resources. GRID resulted in relaxed structures that had slightly higher backbone RMSDs compared to Backrub relative to the starting crystal structures. The average RMSD was 0.25 ± 0.02 Å for GRID versus 0.14 ± 0.04 Å for Backrub. These relatively minor deviations indicate that both algorithms generate structures that retain their original topologies, as expected given the nature of the algorithms. © 2012 Wiley Periodicals, Inc.     

Differential ordering of the protein backbone and side chains during protein folding revealed by site-specific recombinant infrared probes

The time scale for ordering of the polypeptide backbone relative to the side chains is a critical issue in protein folding. The interplay between ordering of the backbone and ordering of the side chains is particularly important for the formation of β-sheet structures, as the polypeptide chain searches for the native stabilizing cross-strand interactions. We have studied these issues in the N-terminal domain of protein L9 (NTL9), a model protein with mixed α/β structure. We have developed a general approach for introducing site-specific IR probes for the side chains (azide) and backbone ((13)C═(18)O) using recombinant protein expression. Temperature-jump time-resolved IR spectroscopy combined with site-specific labeling enables independent measurement of the respective backbone and side-chain dynamics with single residue resolution. We have found that side-chain ordering in a key region of the β-sheet structure occurs on a slower time scale than ordering of the backbone during the folding of NTL9, likely as a result of the transient formation of non-native side-chain interactions.     

Energy-based reconstruction of a protein backbone from its alpha-carbon trace by a Monte-Carlo method

An automatic procedure is proposed for reconstruction of a protein backbone from its C(alpha)-trace; it is based on optimization of a simplified energy function of a peptide backbone, given its alpha-carbon trace. The energy is expressed as a sum of the energies of interaction between backbone peptide groups that are not neighbors in the sequence, the energies of local interactions within all amino acid residues, and a harmonic penalty function accounting for the conservation of standard bond angles. The energy of peptide group interactions is calculated using the assumption that each peptide group acts as a point dipole. For local interaction energy, use is made of a two-dimensional Fourier series expansion of the energies of model terminally blocked amino acid residues, calculated with the Empirical Conformational Energy Program for Peptides (ECEPP/3) force field in the angles lambda((1)) and lambda((2)) defining the rotation of peptide groups adjacent to a C(alpha) carbon atom about the corresponding C(alpha) em leader C(alpha) virtual-bond axes. To explore all possible rotations of peptide groups within a fixed C(alpha)-trace, a Monte Carlo search is carried out. The initial lambda angles are calculated by aligning the dipoles of the peptide groups that are close in space, subject to the condition of favorable local interactions. After the Monte Carlo search is accomplished with the simplified energy function, the energy of the structure is minimized with the ECEPP/3 force field, with imposition of distance constraints corresponding to the initial C(alpha)-trace geometry. The procedure was tested on model alpha-helices and beta-sheets, as well as on the crystal structure of the immunoglobulin binding protein (PDB code: 1IGD, an alpha/beta protein). In all cases, complete backbone geometry was reconstructed with a root-mean-square (rms) deviation of 0.5 A from the all-atom target structure.     

Transplant-insert-constrain-relax-assemble (TICRA): protein-ligand complex structure modeling and application to kinases

We introduce TICRA (transplant-insert-constrain-relax-assemble), a method for modeling the structure of unknown protein-ligand complexes using the X-ray crystal structures of homologous proteins and ligands with known activity. We present results from modeling the structures of protein kinase-inhibitor complexes using p38 and Lck as examples. These examples show that the TICRA method may be used prospectively to create and refine models for protein kinase-inhibitor complexes with an overall backbone rmsd of less than 0.75 ? for the kinase domain, when compared to published X-ray crystal structures. Further refinement of the models of the kinase domains of p38 and Lck in complex with their cognate ligands from the published crystal structures was able to improve the rmsd's of the model complexes to below 0.5 ?. Our results show that TICRA is a useful approach to the problem of structure-based drug design in cases where little structural information is available for the target proteins and the binding mode of active compounds is unknown.     

3D structure determination of the Crh protein from highly ambiguous solid-state NMR restraints

In a wide variety of proteins, insolubility presents a challenge to structural biology, as X-ray crystallography and liquid-state NMR are unsuitable. Indeed, no general approach is available as of today for studying the three-dimensional structures of membrane proteins and protein fibrils. We here demonstrate, at the example of the microcrystalline model protein Crh, how high-resolution 3D structures can be derived from magic-angle spinning solid-state NMR distance restraints for fully labeled protein samples. First, we show that proton-mediated rare-spin correlation spectra, as well as carbon-13 spin diffusion experiments, provide enough short, medium, and long-range structural restraints to obtain high-resolution structures of this 2 x 10.4 kDa dimeric protein. Nevertheless, the large number of 13C/15N spins present in this protein, combined with solid-state NMR line widths of about 0.5-1 ppm, induces substantial ambiguities in resonance assignments, preventing 3D structure determination by using distance restraints uniquely assigned on the basis of their chemical shifts. In the second part, we thus demonstrate that an automated iterative assignment algorithm implemented in a dedicated solid-state NMR version of the program ARIA permits to resolve the majority of ambiguities and to calculate a de novo 3D structure from highly ambiguous solid-state NMR data, using a unique fully labeled protein sample. We present, using distance restraints obtained through the iterative assignment process, as well as dihedral angle restraints predicted from chemical shifts, the 3D structure of the fully labeled Crh dimer refined at a root-mean-square deviation of 1.33 A.     

The Double‐Histidine Cu2+‐Binding Motif: A Highly Rigid,Site‐Specific Spin Probe for Electron Spin Resonance Distance Measurements

The development of ESR methods that measure long‐range distance distributions has advanced biophysical research. However, the spin labels commonly employed are highly flexible, which leads to ambiguity in relating ESR measurements to protein‐backbone structure. Herein we present the double‐histidine (dHis) Cu²⁺‐binding motif as a rigid spin probe for double electron–electron resonance (DEER) distance measurements. The spin label is assembled in situ from natural amino acid residues and a metal salt, requires no postexpression synthetic modification, and provides distance distributions that are dramatically narrower than those found with the commonly used protein spin label. Simple molecular modeling based on an X‐ray crystal structure of an unlabeled protein led to a predicted most probable distance within 0.5 Å of the experimental value. Cu²⁺ DEER with the dHis motif shows great promise for the resolution of precise, unambiguous distance constraints that relate directly to protein‐backbone structure and flexibility.     

The Double‐Histidine Cu2+‐Binding Motif: A Highly Rigid,Site‐Specific Spin Probe for Electron Spin Resonance Distance Measurements

The development of ESR methods that measure long‐range distance distributions has advanced biophysical research. However, the spin labels commonly employed are highly flexible, which leads to ambiguity in relating ESR measurements to protein‐backbone structure. Herein we present the double‐histidine (dHis) Cu²⁺‐binding motif as a rigid spin probe for double electron–electron resonance (DEER) distance measurements. The spin label is assembled in situ from natural amino acid residues and a metal salt, requires no postexpression synthetic modification, and provides distance distributions that are dramatically narrower than those found with the commonly used protein spin label. Simple molecular modeling based on an X‐ray crystal structure of an unlabeled protein led to a predicted most probable distance within 0.5 Å of the experimental value. Cu²⁺ DEER with the dHis motif shows great promise for the resolution of precise, unambiguous distance constraints that relate directly to protein‐backbone structure and flexibility.     

Comparative structural analysis of a histone-like protein from Spiroplasma melliferum in the crystalline state and in solution

A solution of a histone-like protein from Spiroplasma melliferum (HUSpm) was examined by small-angle X-ray scattering (SAXS). The experimental SAXS curve was compared with those calculated for the HUSpm structures from the PDB databank obtained by both X-ray diffraction analysis and nuclear magnetic resonance spectroscopy. The model of the HUSpm structure in solution, which best agrees with the experimental SAXS data, has a shorter distance between the centers of mass of the HUSpm monomers compared to the crystal structure, indicating that the HUSpm monomers can be located closer to each other in solution than in the crystalline state.     

Methyl proton contacts obtained using heteronuclear through-bond transfers in solid-state NMR spectroscopy

A two-dimensional proton-mediated carbon-carbon correlation experiment that relies on through-bond heteronuclear magnetization transfers is demonstrated in the context of solid-state NMR of proteins. This new experiment, dubbed J-CHHC by analogy to the previously developed dipolar CHHC techniques, is shown to provide selective and sensitive correlations in the methyl region of 2D spectra of crystalline organic compounds. The method is then demonstrated on a microcrystalline sample of the dimeric protein Crh (2 x 10.4 kDa). A total of 34 new proton-proton contacts involving side-chain methyl groups were observed in the J-CHHC spectrum, which had not been observed with the conventional experiment. The contacts were then used as additional distance restraints for the 3D structure determination of this microcrystalline protein. Upon addition of these new distance restraints, which are in large part located in the hydrophobic core of the protein, the root-mean-square deviation with respect to the X-ray structure of the backbone atom coordinates of the 10 best conformers of the new ensemble of structures is reduced from 1.8 to 1.1 A.     

Increasing the accuracy of solution NMR structures of membrane proteins by application of residual dipolar couplings. High-resolution structure of outer membrane protein A

The structure determination of membrane proteins is one of the most challenging applications of solution NMR spectroscopy. The paucity of distance information available from the highly deuterated proteins employed requires new approaches in structure determination. Here we demonstrate that significant improvement in the structure accuracy of the membrane protein OmpA can be achieved by refinement with residual dipolar couplings (RDCs). The application of charged polyacrylamide gels allowed us to obtain two alignments and accurately measure numerous heteronuclear dipolar couplings. Furthermore, we have demonstrated that using a large set of RDCs in the refinement can yield a structure with 1 A rms deviation to the backbone of the high-resolution crystal structure. Our simulations with various data sets indicate that dipolar couplings will be critical for obtaining accurate structures of membrane proteins.