Why the Energy Landscape of Barnase Is Hierarchical

We have used NMR and computational methods to characterize the dynamics of the ribonuclease barnase over a wide range of timescales in free and inhibitor-bound states. Using temperature- and denaturant-dependent measurements of chemical shift, we show that barnase undergoes frequent and highly populated hinge bending. Using relaxation dispersion, we characterize a slower and less populated motion with a rate of 750 ± 200 s⁻¹, involving residues around the lip of the active site, which occurs in both free and bound states and therefore suggests conformational selection. Normal mode calculations characterize correlated hinge bending motions on a very rapid timescale. These three measurements are combined with previous measurements and molecular dynamics calculations on barnase to characterize its dynamic landscape on timescales from picoseconds to milliseconds and length scales from 0.1 to 2.5 nm. We show that barnase has two different large-scale fluctuations: one on a timescale of 10⁻⁹−10⁻⁶ s that has no free energy barrier and is a hinge bending that is determined by the architecture of the protein; and one on a timescale of milliseconds (i.e., 750 s⁻¹) that has a significant free energy barrier and starts from a partially hinge-bent conformation. These two motions can be described as hierarchical, in that the more highly populated faster motion provides a platform for the slower (less probable) motion. The implications are discussed. The use of temperature and denaturant is suggested as a simple and general way to characterize motions on the intermediate ns-μs timescale.

Development and Validation of 2D Difference Intensity Analysis for Chemical Library Screening by Protein-Detected NMR Spectroscopy

An academic chemical screening approach was developed by using 2D protein-detected NMR, and a 352-chemical fragment library was screened against three different protein targets. The approach was optimized against two protein targets with known ligands: CXCL12 and BRD4. Principal component analysis reliably identified compounds that induced nonspecific NMR crosspeak broadening but did not unambiguously identify ligands with specific affinity (hits). For improved hit detection, a novel scoring metric-difference intensity analysis (DIA)-was devised that sums all positive and negative intensities from 2D difference spectra. Applying DIA quickly discriminated potential ligands from compounds inducing nonspecific NMR crosspeak broadening and other nonspecific effects. Subsequent NMR titrations validated chemotypes important for binding to CXCL12 and BRD4. A novel target, mitochondrial fission protein Fis1, was screened, and six hits were identified by using DIA. Screening these diverse protein targets identified quinones and catechols that induced nonspecific NMR crosspeak broadening, hampering NMR analyses, but are currently not computationally identified as pan-assay interference compounds. The results established a streamlined screening workflow that can easily be scaled and adapted as part of a larger screening pipeline to identify fragment hits and assess relative binding affinities in the range of 0.3-1.6 mm. DIA could prove useful in library screening and other applications in which NMR chemical shift perturbations are measured.

Multi-probe relaxation dispersion measurements increase sensitivity to protein dynamics

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion measurements are a valuable tool for the characterization of structural transitions on the micro-millisecond timescale. While the measurement of (15)N relaxation dispersion is now routine, the measurements with alternative nuclei remain limited. Here we report (15)N as well as (1)H R2 relaxation dispersion measurements of the N23PP/S148A "dynamic knockout" mutant of dihydrofolate reductase. The (1)H dispersion measurements are complementary to (15)N data as many additional residues are observed to have dispersive behavior for the (1)H nucleus. Simultaneous fitting of the dispersion profiles for the two nuclei increases the accuracy of exchange parameters determined for individual residues and clustered groups of residues. The different sensitivity of the two nuclei to changes in backbone torsional angles, ring currents, and hydrogen bonding effects provides important insights into the nature of the structural changes that take place during the exchange process. We observe clear evidence of direct and indirect hydrogen bond effects for the (15)N and (1)H chemical shift changes in the active-site, modulation of ring current shielding in the CD-loop and backbone torsional changes in a cluster of residues associated with the C-terminus. This work demonstrates the power of combined (1)H and (15)N probes for the study of backbone dynamics on the micro-millisecond timescale though the analysis of chemical shift changes.

Correction of erroneously packed protein's side chains in the NMR structure based on ab initio chemical shift calculations

In this work, protein side chain (1)H chemical shifts are used as probes to detect and correct side-chain packing errors in protein's NMR structures through structural refinement. By applying the automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) method for ab initio calculation of chemical shifts, incorrect side chain packing was detected in the NMR structures of the Pin1 WW domain. The NMR structure is then refined by using molecular dynamics simulation and the polarized protein-specific charge (PPC) model. The computationally refined structure of the Pin1 WW domain is in excellent agreement with the corresponding X-ray structure. In particular, the use of the PPC model yields a more accurate structure than that using the standard (nonpolarizable) force field. For comparison, some of the widely used empirical models for chemical shift calculations are unable to correctly describe the relationship between the particular proton chemical shift and protein structures. The AF-QM/MM method can be used as a powerful tool for protein NMR structure validation and structural flaw detection.

1H NMR spectra part 31: 1H chemical shifts of amides in DMSO solvent

The (1)H chemical shifts of 48 amides in DMSO solvent are assigned and presented. The solvent shifts Δδ (DMSO-CDCl3 ) are large (1-2 ppm) for the NH protons but smaller and negative (-0.1 to -0.2 ppm) for close range protons. A selection of the observed solvent shifts is compared with calculated shifts from the present model and from GIAO calculations. Those for the NH protons agree with both calculations, but other solvent shifts such as Δδ(CHO) are not well reproduced by the GIAO calculations. The (1)H chemical shifts of the amides in DMSO were analysed using a functional approach for near ( ≤ 3 bonds removed) protons and the electric field, magnetic anisotropy and steric effect of the amide group for more distant protons. The chemical shifts of the NH protons of acetanilide and benzamide vary linearly with the π density on the αN and βC atoms, respectively. The C=O anisotropy and steric effect are in general little changed from the values in CDCl3. The effects of substituents F, Cl, Me on the NH proton shifts are reproduced. The electric field coefficient for the protons in DMSO is 90% of that in CDCl3. There is no steric effect of the C=O oxygen on the NH proton in an NH…O=C hydrogen bond. The observed deshielding is due to the electric field effect. The calculated chemical shifts agree well with the observed shifts (RMS error of 0.106 ppm for the data set of 257 entries).

Charge invariant protein-water relaxation in GB1 via ultrafast tryptophan fluorescence

The protein–water interface is a critical determinant of protein structure and function, yet the precise nature of dynamics in this complex system remains elusive. Tryptophan fluorescence has become the probe of choice for such dynamics on the picosecond time scale (especially via fluorescence “upconversion”). In the absence of ultrafast (“quasi-static”) quenching from nearby groups, the TDFSS (time-dependent fluorescence Stokes shift) for exposed Trp directly reports on dipolar relaxation near the interface (both water and polypeptide). The small protein GB1 contains a single Trp (W43) of this type, and its structure is refractory to pH above 3. Thus, it can be used to examine the dependence of dipolar relaxation upon charge reconfiguration with titration. Somewhat surprisingly, the dipolar dynamics in the 100 fs to 100 ps range were unchanged with pH, although nanosecond yield, rates, and access all changed. These results were rationalized with the help of molecular dynamics (including QM-MM) simulations that reveal a balancing, sometimes even countervailing influence of protein and water dipoles. Interestingly, these simulations also showed the dominant influence of water molecules which are associated with the protein interface for up to 30 ps yet free to rotate at approximately “bulk” water rates.

Protein dielectric constants determined from NMR chemical shift perturbations

Understanding the connection between protein structure and function requires a quantitative understanding of electrostatic effects. Structure-based electrostatic calculations are essential for this purpose, but their use has been limited by a long-standing discussion on which value to use for the dielectric constants (ε(eff) and ε(p)) required in Coulombic and Poisson-Boltzmann models. The currently used values for ε(eff) and ε(p) are essentially empirical parameters calibrated against thermodynamic properties that are indirect measurements of protein electric fields. We determine optimal values for ε(eff) and ε(p) by measuring protein electric fields in solution using direct detection of NMR chemical shift perturbations (CSPs). We measured CSPs in 14 proteins to get a broad and general characterization of electric fields. Coulomb's law reproduces the measured CSPs optimally with a protein dielectric constant (ε(eff)) from 3 to 13, with an optimal value across all proteins of 6.5. However, when the water-protein interface is treated with finite difference Poisson-Boltzmann calculations, the optimal protein dielectric constant (ε(p)) ranged from 2 to 5 with an optimum of 3. It is striking how similar this value is to the dielectric constant of 2-4 measured for protein powders and how different it is from the ε(p) of 6-20 used in models based on the Poisson-Boltzmann equation when calculating thermodynamic parameters. Because the value of ε(p) = 3 is obtained by analysis of NMR chemical shift perturbations instead of thermodynamic parameters such as pK(a) values, it is likely to describe only the electric field and thus represent a more general, intrinsic, and transferable ε(p) common to most folded proteins.

Using chemical shift perturbation to characterise ligand binding

Chemical shift perturbation (CSP, chemical shift mapping or complexation-induced changes in chemical shift, CIS) follows changes in the chemical shifts of a protein when a ligand is added, and uses these to determine the location of the binding site, the affinity of the ligand, and/or possibly the structure of the complex. A key factor in determining the appearance of spectra during a titration is the exchange rate between free and bound, or more specifically the off-rate koff. When koff is greater than the chemical shift difference between free and bound, which typically equates to an affinity Kd weaker than about 3μM, then exchange is fast on the chemical shift timescale. Under these circumstances, the observed shift is the population-weighted average of free and bound, which allows Kd to be determined from measurement of peak positions, provided the measurements are made appropriately. (1)H shifts are influenced to a large extent by through-space interactions, whereas (13)Cα and (13)Cβ shifts are influenced more by through-bond effects. (15)N and (13)C' shifts are influenced both by through-bond and by through-space (hydrogen bonding) interactions. For determining the location of a bound ligand on the basis of shift change, the most appropriate method is therefore usually to measure (15)N HSQC spectra, calculate the geometrical distance moved by the peak, weighting (15)N shifts by a factor of about 0.14 compared to (1)H shifts, and select those residues for which the weighted shift change is larger than the standard deviation of the shift for all residues. Other methods are discussed, in particular the measurement of (13)CH3 signals. Slow to intermediate exchange rates lead to line broadening, and make Kd values very difficult to obtain. There is no good way to distinguish changes in chemical shift due to direct binding of the ligand from changes in chemical shift due to allosteric change. Ligand binding at multiple sites can often be characterised, by simultaneous fitting of many measured shift changes, or more simply by adding substoichiometric amounts of ligand. The chemical shift changes can be used as restraints for docking ligand onto protein. By use of quantitative calculations of ligand-induced chemical shift changes, it is becoming possible to determine not just the position but also the orientation of ligands.

Protein backbone and sidechain torsion angles predicted from NMR chemical shifts using artificial neural networks

A new program, TALOS-N, is introduced for predicting protein backbone torsion angles from NMR chemical shifts. The program relies far more extensively on the use of trained artificial neural networks than its predecessor, TALOS+. Validation on an independent set of proteins indicates that backbone torsion angles can be predicted for a larger, ≥90 % fraction of the residues, with an error rate smaller than ca 3.5 %, using an acceptance criterion that is nearly two-fold tighter than that used previously, and a root mean square difference between predicted and crystallographically observed (ϕ, ψ) torsion angles of ca 12º. TALOS-N also reports sidechain χ(1) rotameric states for about 50 % of the residues, and a consistency with reference structures of 89 %. The program includes a neural network trained to identify secondary structure from residue sequence and chemical shifts.

NightShift: NMR shift inference by general hybrid model training--a framework for NMR chemical shift prediction

Background

NMR chemical shift prediction plays an important role in various applications in computational biology. Among others, structure determination, structure optimization, and the scoring of docking results can profit from efficient and accurate chemical shift estimation from a three-dimensional model.

A variety of NMR chemical shift prediction approaches have been presented in the past, but nearly all of these rely on laborious manual data set preparation and the training itself is not automatized, making retraining the model, e.g., if new data is made available, or testing new models a time-consuming manual chore.

Results

In this work, we present the framework NightShift (NMR Shift Inference by General Hybrid Model Training), which enables automated data set generation as well as model training and evaluation of protein NMR chemical shift prediction.

In addition to this main result – the NightShift framework itself – we describe the resulting, automatically generated, data set and, as a proof-of-concept, a random forest model called Spinster that was built using the pipeline.

Conclusion

By demonstrating that the performance of the automatically generated predictors is at least en par with the state of the art, we conclude that automated data set and predictor generation is well-suited for the design of NMR chemical shift estimators.

The framework can be downloaded from https://bitbucket.org/akdehof/nightshift. It requires the open source Biochemical Algorithms Library (BALL), and is available under the conditions of the GNU Lesser General Public License (LGPL). We additionally offer a browser-based user interface to our NightShift instance employing the Galaxy framework via https://ballaxy.bioinf.uni-sb.de/.

1H NMR spectra. Part 30(+): 1H chemical shifts in amides and the magnetic anisotropy, electric field and steric effects of the amide group

The (1)H spectra of 37 amides in CDCl(3) solvent were analysed and the chemical shifts obtained. The molecular geometries and conformational analysis of these amides were considered in detail. The NMR spectral assignments are of interest, e.g. the assignments of the formamide NH(2) protons reverse in going from CDCl(3) to more polar solvents. The substituent chemical shifts of the amide group in both aliphatic and aromatic amides were analysed using an approach based on neural network data for near (≤3 bonds removed) protons and the electric field, magnetic anisotropy, steric and for aromatic systems π effects of the amide group for more distant protons. The electric field is calculated from the partial atomic charges on the N.C═O atoms of the amide group. The magnetic anisotropy of the carbonyl group was reproduced with the asymmetric magnetic anisotropy acting at the midpoint of the carbonyl bond. The values of the anisotropies Δχ(parl) and Δχ(perp) were for the aliphatic amides 10.53 and -23.67 (×10(-6) Å(3)/molecule) and for the aromatic amides 2.12 and -10.43 (×10(-6) Å(3)/molecule). The nitrogen anisotropy was 7.62 (×10(-6) Å(3)/molecule). These values are compared with previous literature values. The (1)H chemical shifts were calculated from the semi-empirical approach and also by gauge-independent atomic orbital calculations with the density functional theory method and B3LYP/6-31G(++) (d,p) basis set. The semi-empirical approach gave good agreement with root mean square error of 0.081 ppm for the data set of 280 entries. The gauge-independent atomic orbital approach was generally acceptable, but significant errors (ca. 1 ppm) were found for the NH and CHO protons and also for some other protons.

Assessing and refining molecular dynamics simulations of proteins with nuclear magnetic resonance data

The sophistication of the force fields, algorithms and hardware used for molecular dynamics (MD) simulations of proteins is continuously increasing. No matter how advanced the methodology, however, it is essential to evaluate the appropriateness of the structures sampled in a simulation by comparison with quantitative experimental data. Solution nuclear magnetic resonance (NMR) data are particularly useful for checking the quality of protein simulations, as they provide both structural and dynamic information on a variety of temporal and spatial scales. Here, various features and implications of using NMR data to validate and bias MD simulations are outlined, including an overview of the different types of NMR data that report directly on structural properties and of relevant simulation techniques. The focus throughout is on how to properly account for conformational averaging, particularly within the context of the assumptions inherent in the relationships that link NMR data to structural properties.

Interpreting protein structural dynamics from NMR chemical shifts

In this investigation, semiempirical NMR chemical shift prediction methods are used to evaluate the dynamically averaged values of backbone chemical shifts obtained from unbiased molecular dynamics (MD) simulations of proteins. MD-averaged chemical shift predictions generally improve agreement with experimental values when compared to predictions made from static X-ray structures. Improved chemical shift predictions result from population-weighted sampling of multiple conformational states and from sampling smaller fluctuations within conformational basins. Improved chemical shift predictions also result from discrete changes to conformations observed in X-ray structures, which may result from crystal contacts, and are not always reflective of conformational dynamics in solution. Chemical shifts are sensitive reporters of fluctuations in backbone and side chain torsional angles, and averaged (1)H chemical shifts are particularly sensitive reporters of fluctuations in aromatic ring positions and geometries of hydrogen bonds. In addition, poor predictions of MD-averaged chemical shifts can identify spurious conformations and motions observed in MD simulations that may result from force field deficiencies or insufficient sampling and can also suggest subsets of conformational space that are more consistent with experimental data. These results suggest that the analysis of dynamically averaged NMR chemical shifts from MD simulations can serve as a powerful approach for characterizing protein motions in atomistic detail.

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 impact of hydrogen bonding on amide 1H chemical shift anisotropy studied by cross-correlated relaxation and liquid crystal NMR spectroscopy

Site-specific ¹H chemical shift anisotropy (CSA) tensors have been derived for the well-ordered backbone amide moieties in the B3 domain of protein G (GB3). Experimental input data include residual chemical shift anisotropy (RCSA), measured in six mutants that align differently relative to the static magnetic field when dissolved in a liquid crystalline Pf1 suspension, and cross-correlated relaxation rates between the ¹H^N CSA tensor and either the ¹H−¹⁵N, the ¹H−¹³C′, or the ¹H−¹³C^α dipolar interactions. Analyses with the assumption that the ¹H^N CSA tensor is symmetric with respect to the peptide plane (three-parameter fit) or without this premise (five-parameter fit) yield very similar results, confirming the robustness of the experimental input data, and that, to a good approximation, one of the principal components orients orthogonal to the peptide plane. ¹H^N CSA tensors are found to deviate strongly from axial symmetry, with the most shielded tensor component roughly parallel to the N−H vector, and the least shielded component orthogonal to the peptide plane. DFT calculations on pairs of N-methyl acetamide and acetamide in H-bonded geometries taken from the GB3 X-ray structure correlate with experimental data and indicate that H-bonding effects dominate variations in the ¹H^N CSA. Using experimentally derived ¹H^N CSA tensors, the optimal relaxation interference effect needed for narrowest ¹H^N TROSY line widths is found at ∼1200 MHz.

Site-specific backbone amide (15)N chemical shift anisotropy tensors in a small protein from liquid crystal and cross-correlated relaxation measurements

Site-specific (15)N chemical shift anisotropy (CSA) tensors have been derived for the well-ordered backbone amide (15)N nuclei in the B3 domain of protein G (GB3) from residual chemical shift anisotropy (RCSA) measured in six different mutants that retain the native structure but align differently relative to the static magnetic field when dissolved in a liquid crystalline Pf1 suspension. This information is complemented by measurement of cross-correlated relaxation rates between the (15)N CSA tensor and either the (15)N-(1)H or (15)N-(13)C' dipolar interaction. In agreement with recent solid state NMR measurements, the (15)N CSA tensors exhibit only a moderate degree of variation from averaged values, but have larger magnitudes in alpha-helical (-173 +/- 7 ppm) than in beta-sheet (-162 +/- 6 ppm) residues, a finding also confirmed by quantum computations. The orientations of the least shielded tensor component cluster tightly around an in-peptide-plane vector that makes an angle of 19.6 +/- 2.5 degrees with the N-H bond, with the asymmetry of the (15)N CSA tensor being slightly smaller in alpha-helix (eta = 0.23 +/- 0.17) than in beta-sheet (eta = 0.31 +/- 0.11). The residue-specific (15)N CSA values are validated by improved agreement between computed and experimental (15)N R(1rho) relaxation rates measured for (15)N-{(2)H} sites in GB3, which are dominated by the CSA mechanism. Use of residue-specific (15)N CSA values also results in more uniform generalized order parameters, S(2), and predicts considerable residue-by-residue variations in the magnetic field strengths where TROSY line narrowing is most effective.

Structural analysis of the binding of the diquaternary pyridophenazine derivative dqdppn to B-DNA oligonucleotides

The interaction of the ethylene-bipyridyldiylium-naphthaphenazine dication, dqdppn, with several hexa- and octanucleotide duplexes has been studied using CD and NMR. Taken together, these studies reveal that with the hexanucleotide, dqdppn intercalates into the terminal base pair, and causes a large twisting of the terminal base pair. In contrast, with all three octanucleotides, dqdppn intercalates more centrally within the sequence. The NMR-derived structures of two of the binding complexes demonstrate that dqdppn intercalates from the major groove in an unusual 'side-on' geometry, rather than threading through the helix. An analysis of these results indicates that the preferred binding site is not sequence-specific, but primarily at the most conformationally flexible DNA step.

Pressure-dependent structure changes in barnase on ligand binding reveal intermediate rate fluctuations

In this work we measured 1H NMR chemical shifts for the ribonuclease barnase at pressures from 3 MPa to 200 MPa, both free and bound to d(CGAC). Shift changes with pressure were used as restraints to determine the change in structure with pressure. Free barnase is compressed by approximately 0.7%. The largest changes are on the ligand-binding face close to Lys-27, which is the recognition site for the cleaved phosphate bond. This part of the protein also contains the buried water molecules. In the presence of d(CGAC), the compressibility is reduced by approximately 70% and the region of structural change is altered: the ligand-binding face is now almost incompressible, whereas changes occur at the opposite face. Because compressibility is proportional to mean square volume fluctuation, we conclude that in free barnase, volume fluctuation is largest close to the active site, but when the inhibitor is bound, the fluctuations become much smaller and are located mainly on the opposite face. The timescale of the fluctuations is nanoseconds to microseconds, consistent with the degree of ordering required for the fluctuations, which are intermediate between rapid uncorrelated side-chain dynamics and slow conformational transitions. The high-pressure technique is therefore useful for characterizing motions on this relatively inaccessible timescale.

TALOS+: a hybrid method for predicting protein backbone torsion angles from NMR chemical shifts

NMR chemical shifts in proteins depend strongly on local structure. The program TALOS establishes an empirical relation between 13C, 15N and 1H chemical shifts and backbone torsion angles phi and psi (Cornilescu et al. J Biomol NMR 13 289-302, 1999). Extension of the original 20-protein database to 200 proteins increased the fraction of residues for which backbone angles could be predicted from 65 to 74%, while reducing the error rate from 3 to 2.5%. Addition of a two-layer neural network filter to the database fragment selection process forms the basis for a new program, TALOS+, which further enhances the prediction rate to 88.5%, without increasing the error rate. Excluding the 2.5% of residues for which TALOS+ makes predictions that strongly differ from those observed in the crystalline state, the accuracy of predicted phi and psi angles, equals +/-13 degrees . Large discrepancies between predictions and crystal structures are primarily limited to loop regions, and for the few cases where multiple X-ray structures are available such residues are often found in different states in the different structures. The TALOS+ output includes predictions for individual residues with missing chemical shifts, and the neural network component of the program also predicts secondary structure with good accuracy.

Improving NMR protein structure quality by Rosetta refinement: a molecular replacement study

The structure of human protein HSPC034 has been determined by both solution nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography. Refinement of the NMR structure ensemble, using a Rosetta protocol in the absence of NMR restraints, resulted in significant improvements not only in structure quality, but also in molecular replacement (MR) performance with the raw X-ray diffraction data using MOLREP and Phaser. This method has recently been shown to be generally applicable with improved MR performance demonstrated for eight NMR structures refined using Rosetta (Qian et al., Nature 2007;450:259-264). Additionally, NMR structures of HSPC034 calculated by standard methods that include NMR restraints have improvements in the RMSD to the crystal structure and MR performance in the order DYANA, CYANA, XPLOR-NIH, and CNS with explicit water refinement (CNSw). Further Rosetta refinement of the CNSw structures, perhaps due to more thorough conformational sampling and/or a superior force field, was capable of finding alternative low energy protein conformations that were equally consistent with the NMR data according to the Recall, Precision, and F-measure (RPF) scores. On further examination, the additional MR-performance shortfall for NMR refined structures as compared with the X-ray structure were attributed, in part, to crystal-packing effects, real structural differences, and inferior hydrogen bonding in the NMR structures. A good correlation between a decrease in the number of buried unsatisfied hydrogen-bond donors and improved MR performance demonstrates the importance of hydrogen-bond terms in the force field for improving NMR structures. The superior hydrogen-bond network in Rosetta-refined structures demonstrates that correct identification of hydrogen bonds should be a critical goal of NMR structure refinement. Inclusion of nonbivalent hydrogen bonds identified from Rosetta structures as additional restraints in the structure calculation results in NMR structures with improved MR performance.

Determination of protein structures in the solid state from NMR chemical shifts

Solid-state NMR spectroscopy does not require proteins to form crystalline or soluble samples and can thus be applied under a variety of conditions, including precipitates, gels, and microcrystals. It has recently been shown that NMR chemical shifts can be used to determine the structures of the native states of proteins in solution. By considering the cases of two proteins, GB1 and SH3, we provide an initial demonstration here that this type of approach can be extended to the use of solid-state NMR chemical shifts to obtain protein structures in the solid state without the need for measuring interatomic distances.

De novo protein structure generation from incomplete chemical shift assignments

NMR chemical shifts provide important local structural information for proteins. Consistent structure generation from NMR chemical shift data has recently become feasible for proteins with sizes of up to 130 residues, and such structures are of a quality comparable to those obtained with the standard NMR protocol. This study investigates the influence of the completeness of chemical shift assignments on structures generated from chemical shifts. The Chemical-Shift-Rosetta (CS-Rosetta) protocol was used for de novo protein structure generation with various degrees of completeness of the chemical shift assignment, simulated by omission of entries in the experimental chemical shift data previously used for the initial demonstration of the CS-Rosetta approach. In addition, a new CS-Rosetta protocol is described that improves robustness of the method for proteins with missing or erroneous NMR chemical shift input data. This strategy, which uses traditional Rosetta for pre-filtering of the fragment selection process, is demonstrated for two paramagnetic proteins and also for two proteins with solid-state NMR chemical shift assignments.

Pressure-induced changes in the solution structure of the GB1 domain of protein G

The solution structure of the GB1 domain of protein G at a pressure of 2 kbar is presented. The structure was calculated as a change from an energy-minimised low-pressure structure using (1)H chemical shifts. Two separate changes can be characterised: a compression/distortion, which is linear with pressure; and a stabilisation of an alternative folded state. On application of pressure, linear chemical shift changes reveal that the backbone structure changes by about 0.2 A root mean square, and is compressed by about 1% overall. The alpha-helix compresses, particularly at the C-terminal end, and moves toward the beta-sheet, while the beta-sheet is twisted, with the corners closest to the alpha-helix curling up towards it. The largest changes in structure are along the second beta-strand, which becomes more twisted. This strand is where the protein binds to IgG. Curved chemical shift changes with pressure indicate that high pressure also populates an alternative structure with a distortion towards the C-terminal end of the helix, which is likely to be caused by insertion of a water molecule.

Structural change in a B-DNA helix with hydrostatic pressure

Study of the effects of pressure on macromolecular structure improves our understanding of the forces governing structure, provides details on the relevance of cavities and packing in structure, increases our understanding of hydration and provides a basis to understand the biology of high-pressure organisms. A study of DNA, in particular, helps us to understand how pressure can affect gene activity. Here we present the first high-resolution experimental study of B-DNA structure at high pressure, using NMR data acquired at pressures up to 200 MPa (2 kbar). The structure of DNA compresses very little, but is distorted so as to widen the minor groove, and to compress hydrogen bonds, with AT pairs compressing more than GC pairs. The minor groove changes are suggested to lead to a compression of the hydration water in the minor groove.

Consistent blind protein structure generation from NMR chemical shift data

Protein NMR chemical shifts are highly sensitive to local structure. A robust protocol is described that exploits this relation for de novo protein structure generation, using as input experimental parameters the (13)C(alpha), (13)C(beta), (13)C', (15)N, (1)H(alpha) and (1)H(N) NMR chemical shifts. These shifts are generally available at the early stage of the traditional NMR structure determination process, before the collection and analysis of structural restraints. The chemical shift based structure determination protocol uses an empirically optimized procedure to select protein fragments from the Protein Data Bank, in conjunction with the standard ROSETTA Monte Carlo assembly and relaxation methods. Evaluation of 16 proteins, varying in size from 56 to 129 residues, yielded full-atom models that have 0.7-1.8 A root mean square deviations for the backbone atoms relative to the experimentally determined x-ray or NMR structures. The strategy also has been successfully applied in a blind manner to nine protein targets with molecular masses up to 15.4 kDa, whose conventional NMR structure determination was conducted in parallel by the Northeast Structural Genomics Consortium. This protocol potentially provides a new direction for high-throughput NMR structure determination.

Cooperative Hydrogen Bonding Effects Are Key Determinants of Backbone Amide Proton Chemical Shifts in Proteins

A computational methodology for backbone amide proton chemical shift (delta(H)) predictions based on ab initio quantum mechanical treatment of part of the protein is presented. The method is used to predict and interpret 13 delta(H) values in protein G and ubiquitin. The predicted amide-amide delta(H) values are within 0.6 ppm of experiment, with a root-mean-square deviation (RMSD) of 0.3 ppm. We show that while the hydrogen bond geometry is the most important delta(H)-determinant, longer-range cooperative effects of extended hydrogen networks make significant contributions to delta(H). We present a simple model that accurately relates the protein structure to delta(H).

Multiple binding modes of the camptothecin family to DNA oligomers

The binding constants of camptothecin, topotecan and its lactone ring-opened carboxylate derivative to DNA octamers were measured by UV and NMR spectroscopy. The self-association of topotecan (TPT) was also measured. The carboxylate form of TPT binds in the same way as the lactone, but more weakly. Titration of TPT into d(GCGATCGC)2 shows a preferred location stacked onto the terminal G1 base. However, the intermolecular NOEs cannot be reconciled with a single conformation of the complex, and suggest a model of a limited number of conformations in fast exchange. MD calculations on four pairs of starting structures with TPT stacked onto the G1-C8 base pair in different orientations were therefore performed. The use of selected experimental "docking" restraints yielded ten MD trajectories covering a wide conformational space. From a combination of calculated free energies, NOEs and chemical shifts, some of the structures produced could be eliminated, and it is concluded that the data are consistent with two major families of conformations in fast exchange. One of these is the conformation found in a crystal of a TPT/DNA/topoisomerase I ternary complex [Proc. Natl. Acad. Sci. USA 2002, 99, 15 387-15 392].

Flap opening mechanism of HIV-1 protease

The active site of aspartic proteases, such as HIV-1 protease (PR), is covered by one or more flaps, which restrict access to the active site. For HIV-1 PR, X-ray diffraction studies suggested that in the free enzyme the two flaps are packed onto each other loosely in a semi-open conformation, while molecular dynamics (MD) studies observed that the flaps can also separate into open conformations. In this study, the mechanism of flap opening and the structure and dynamics of HIV-1 PR with semi-open and open flap conformations were investigated using molecular dynamics simulations. The flaps showed complex dynamic behavior as two distinct mechanisms of flap opening and various stable flap conformations (semi-open, open and curled) were observed during the simulations. A network of weakly polar interactions between the flaps were proposed to be responsible for stabilizing the semi-open flap conformation. It is hypothesized that such interactions could be responsible for making flap opening a highly sensitive gating mechanism which control access to the active site.

1H chemical shifts in NMR: Part 22-Prediction of the 1H chemical shifts of alcohols, diols and inositols in solution, a conformational and solvation investigation

The (1)H NMR spectra of a number of alcohols, diols and inositols are reported and assigned in CDCl(3), D(2)O and DMSO-d(6) (henceforth DMSO) solutions. These data were used to investigate the effects of the OH group on the (1)H chemical shifts in these molecules and also the effect of changing the solvent. Inspection of the (1)H chemical shifts of those alcohols which were soluble in both CDCl(3) and D(2)O shows that there is no difference in the chemical shifts in the two solvents, provided that the molecules exist in the same conformation in the two solvents. In contrast, DMSO gives rise to significant and specific solvation shifts. The (1)H chemical shifts of these compounds in the three solvents were analysed using the CHARGE model. This model incorporates the electric field, magnetic anisotropy and steric effects of the functional group for long-range protons together with functions for the calculation of the two- and three-bond effects. The long-range effect of the OH group was quantitatively explained without the inclusion of either the C--O bond anisotropy or the C--OH electric field. Differential beta and gamma effects for the 1,2-diol group needed to be included to obtain accurate chemical shift predictions. For DMSO solution the differential solvent shifts were calculated in CHARGE on the basis of a similar model, incorporating two-bond, three-bond and long-range effects. The analyses of the (1)H spectra of the inositols and their derivatives in D(2)O and DMSO solution also gave the ring (1)H,(1)H coupling constants and for DMSO solution the CH--OH couplings and OH chemical shifts. The (1)H,(1)H coupling constants were calculated in the CHARGE program by an extension of the cos(2)phi equation to include the orientation effects of electronegative atoms and the CH--OH couplings by a simple cos(2)phi equation. Comparison of the observed and calculated couplings confirmed the proposed conformations of myo-inositol, chiro-inositol, quebrachitol and allo-inositol. The OH chemical shifts were also calculated in the CHARGE program. Comparison of the observed and calculated OH chemical shifts and CH.OH couplings suggested the existence of intramolecular hydrogen bonding in a myo-inositol derivative.

Computational assignment of protein backbone NMR peaks by efficient bounding and filtering

NMR resonance assignment is one of the key steps in solving an NMR protein structure. The assignment process links resonance peaks to individual residues of the target protein sequence, providing the prerequisite for establishing intra- and inter-residue spatial relationships between atoms. The assignment process is tedious and time-consuming, which could take many weeks. Though there exist a number of computer programs to assist the assignment process, many NMR labs are still doing the assignments manually to ensure quality. This paper presents a new computational method based on the combination of a suite of algorithms for automating the assignment process, particularly the process of backbone resonance peak assignment. We formulate the assignment problem as a constrained weighted bipartite matching problem. While the problem, in the most general situation, is NP-hard, we present an efficient solution based on a branch-and-bound algorithm with effective bounding techniques using two recently introduced approximation algorithms. We also devise a greedy filtering algorithm for reducing the search space. Our experimental results on 70 instances of (pseudo) real NMR data derived from 14 proteins demonstrate that the new solution runs much faster than a recently introduced (exhaustive) two-layer algorithm and recovers more correct peak assignments than the two-layer algorithm. Our result demonstrates that integrating different algorithms can achieve a good tradeoff between backbone assignment accuracy and computation time.

1H chemical shifts in NMR. Part 21--prediction of the 1H chemical shifts of molecules containing the ester group: a modelling and ab initio investigation

The 1H NMR spectra of 24 compounds containing the ester group are given and assigned. These data were used to investigate the effect of the ester group on the 1H chemical shifts in these molecules. These effects were analysed using the CHARGE model, which incorporates the electric field, magnetic anisotropy and steric effects of the functional group for long-range protons together with functions for the calculation of the two- and three-bond effects. The effect of the ester electric field was given by considering the partial atomic charges on the three atoms of the ester group. The anisotropy of the carbonyl group was reproduced with an asymmetric magnetic anisotropy acting at the midpoint of the carbonyl bond with values of Deltachi(parl) and Deltachi(perp) of 10.1 x 10(-30) and -17.1 x 10(-30) cm3 molecule(-1). An aromatic ring current (=0.3 times the benzene ring current) was found to be necessary for pyrone but none for maleic anhydride. This result was confirmed by GIAO calculations. The observed 1H chemical shifts in the above compounds were compared with those calculated by CHARGE and the ab initio GIAO method (B3LYP/6-31G**). For the 24 compounds investigated with 150 1H chemical shifts spanning a range of ca 10 ppm, the CHARGE model gave an excellent r.m.s. error (obs - calc) of <0.1 ppm. The GIAO calculations gave a very reasonable r.m.s. error of ca 0.2 ppm although larger deviations of ca 0.5 ppm were observed for protons near to the electronegative atoms. The accurate predictions of the 1H chemical shifts given by the CHARGE model were used in the conformational analysis of the vinyl esters methyl acrylate and methyl crotonate. An illustration of the use of the CHARGE model in the prediction of the 1H spectrum of a complex organic molecule (benzochromen-6-one) is also given.

GROMACS: Fast, flexible, and free

This article describes the software suite GROMACS (Groningen MAchine for Chemical Simulation) that was developed at the University of Groningen, The Netherlands, in the early 1990s. The software, written in ANSI C, originates from a parallel hardware project, and is well suited for parallelization on processor clusters. By careful optimization of neighbor searching and of inner loop performance, GROMACS is a very fast program for molecular dynamics simulation. It does not have a force field of its own, but is compatible with GROMOS, OPLS, AMBER, and ENCAD force fields. In addition, it can handle polarizable shell models and flexible constraints. The program is versatile, as force routines can be added by the user, tabulated functions can be specified, and analyses can be easily customized. Nonequilibrium dynamics and free energy determinations are incorporated. Interfaces with popular quantum-chemical packages (MOPAC, GAMES-UK, GAUSSIAN) are provided to perform mixed MM/QM simulations. The package includes about 100 utility and analysis programs. GROMACS is in the public domain and distributed (with source code and documentation) under the GNU General Public License. It is maintained by a group of developers from the Universities of Groningen, Uppsala, and Stockholm, and the Max Planck Institute for Polymer Research in Mainz. Its Web site is http://www.gromacs.org.

Protein chemical shifts arising from alpha-helices and beta-sheets depend on solvent exposure

The NMR chemical shifts of certain atomic nuclei in proteins ((1)H(alpha),(13)C(alpha), and (13)C(beta)) depend sensitively on whether or not the amino acid residue is part of a secondary structure (alpha-helix, beta-sheet), and if so, whether it is helix or sheet. The physical origin of the different chemical shifts of atomic nuclei in alpha-helices versus beta-sheets is a problem of long standing. We report that the chemical shift contributions arising from secondary structure (secondary structure shifts) depend strongly on the extent of exposure to solvent. This behavior is observed for (1)H(alpha), (13)C(alpha), and (13)C(beta) (sheet), but not for(13)C(beta) (helix), whose secondary structure shifts are small. When random coil values are subtracted from the chemical shifts of all(1)H(alpha) nuclei (Pro residues excluded) and the residual chemical shifts are summed to plot the mean values against solvent exposure, the results give a funnel-shaped curve that approaches a small value at full-solvent exposure. When chemical shifts are plotted instead against E(local), the electrostatic contribution to conformational energy produced by local dipole-dipole interactions, a well characterized dependence of (1)H(alpha) chemical shifts on E(local) is found. The slope of this plot varies with both the type of amino acid and the extent of solvent exposure. These results indicate that secondary structure shifts are produced chiefly by the electric field of the protein, which is screened by water dipoles at residues in contact with solvent.

NMR structure determination of a modified DNA oligonucleotide containing a new intercalating nucleic acid

The intercalating nucleic acid (INA) presented in this paper is a novel 1-O-(1-pyrenylmethyl)glycerol DNA intercalator that induces high thermal affinity for complementary DNA. The duplex examined contained two INA intercalators, denoted X, inserted directly opposite each other: d(C(1)T(2)C(3)A(4)A(5)C(6)X(7)C(8)A(9)A(10)G(11)C(12)T(13)):d(A(14)G(15)C(16)T(17)-T(18)G(19)X(20)G(21)T(22)T(23)G(24)A(25)G(26)). Unlike most other nucleotide analogues, DNA with INA inserted has a lower affinity for hybridizing to complementary DNA with an INA inserted directly opposite than to complementary unmodified DNA. In this study we used two-dimensional (1)H NMR spectroscopy to determine a high-resolution solution structure of the weak INA-INA duplex. A modified ISPA approach was used to obtain interproton distance bounds from NOESY cross-peak intensities. These distance bounds were used as restraints in molecular dynamics (rMD) calculations. Twenty final structures were generated for the duplex from a B-type DNA starting structure. The root-mean-square deviation (RMSD) of the coordinates for the 20 structures of the complex was 1.95 A. This rather large value, together with broad lines in the area of insertion, reflect the high degree of internal motion in the complex. The determination of the structure revealed that both intercalators were situated in the center of the helix, stacking with each other and the neighboring nucleobases. The intercalation of the INAs caused an unwinding of the helix in the insertion area, creating a ladderlike structure. The structural changes observed upon intercalation were mainly of local character; however, a broadening of the minor groove was found throughout the helix.

Molecular dynamics simulations of HPr under hydrostatic pressure

The histidine-containing protein (HPr) plays an important role in the phosphotransferase system (PTS). The deformations induced on the protein structure at high hydrostatic pressure values (4, 50, 100, 150, and 200 MPa) were previously (H. Kalbitzer, A. Görler, H. Li, P. Dubovskii, A. Hengstenberg, C. Kowolik, H. Yamada, and K. Akasaka, Protein Science 2000, Vol. 9, pp. 693-703) analyzed by NMR experiments: the nonlinear variations of the amide chemical shifts at high pressure values were supposed to arise from induced shifts in the protein conformational equilibrium. Molecular dynamics (MD) simulations are here performed, to analyze the protein internal mobility at 0.1 MPa, and to relate the nonlinear variations of chemical shifts observed at high pressure, to variations in conformational equilibrium. The global features of the protein structure are only slightly modified along the pressure. Nevertheless, the values of the Voronoi residues volumes show that the residues of alpha-helices are more compressed that those belonging to the beta-sheet. The alpha-helices are also displaying the largest internal mobility and deformation in the simulations. The nonlinearity of the 1H chemical shifts, computed from the MD simulation snapshots, is in qualitative agreement with the nonlinearity of the experimentally observed chemical shifts.