Chemical shifts and three-dimensional protein structures

Summary During the past three years it has become possible to compute ab initio the ¹³C, ¹⁵N and ¹⁹F NMR chemical shifts of many sites in native proteins. Chemical shifts are beginning to become a useful supplement to more established methods of solution structure determination, and may find utility in solid-state analysis as well. From ¹³C NMR, information on , and torsions can be obtained, permitting both assignment verification, and structure refinement and prediction. For ¹⁵N, both torsional and hydrogen-bonding effects are important, while for ¹⁹F, chemical shifts are primarily indicators of the local charge field. Chemical shift calculations are still slow, but shielding hypersurfaces — the shift as a function of the dihedral angles that define the molecular conformation — are becoming accessible. Over the next few years, theoretical and computer hardware improvements will enable more routine use of chemical shifts in structural studies, including the study of metal-ligand interactions, the analysis of drug and substrate binding and catalysis, the study of folding/unfolding pathways, as well as the characterization of conformational substates. Rather than simply being a necessary prerequisite for multidimensional NMR, chemical shifts and chemical shift non-equivalence due to folding are now beginning to be useful for structural characterization.     

Computation of relaxation matrix elements from incomplete NOESY data sets

Summary The structural determination of biological molecules in solution by NMR relies on the determination of a set of interatomic distances obtained by measurement of intramolecular nuclear Overhauser effects (NOE). It is shown in this paper that it is possible to obtain the accurate relaxation rate (and hence the interatomic distance) from the direct measurement of a single NOE signal. The precise analysis of a NOESY peak evolution with respect to the mixing time allows the evaluation of the relaxation parameters for the pair of spins under consideration. This is done without any assumption on the relaxation of unmeasured spins, or on the movement of the molecule. The theoretical basis of this method is presented. In order to evaluate the proposed method, a simulated case on the protein BPTI is studied, which shows that the method performs very well even in the case of noisy data sets.     

Hox2 Genes Are Required for Tonotopic Map Precision and Sound Discrimination in the Mouse Auditory Brainstem

1. Download high-res image (303KB)
2. Download full-size image

     

Refined structure of bovine carbonic anhydrase III at 2.0 Å resolution

The three-dimensional structure of bovine carbonic anhydrase III (BCA III) from red skeletal muscle cells has been determined by molecular replacement methods. The structure has been refined at 2.0 Å resolution by both constrained and restrained structure-factor least squares refinement. The current crystallographic R-value is 19.2% and 121 solvent molecules have so far been found associated with the protein. The structure is highly similar to the refined structure of human carbonic anhydrase II. Some differences in amino acid sequence and structure between the two isoenzymes are discussed. In BCA III, Lys 64 and Arg 91 (His 64 and Ile 91 in HCA II) are both pointing out from the active site cavity forming salt bridges with Glu 4 and Asp 72 (His 4 and Asp 72 in HCA II), respectively. However, Arg 67 and Phe 198 (Asn 67 and Leu 198 in HCA II) are oriented towards the zinc ion and significantly reduce the volume of the active site cavity. Phe 198 particularly reduces the size of the substrate binding region at the “deep water” position at the bottom of the cavity and we sugest that this is one of the major reasons for the differences in catalytic properties of isoenzyme III as compared to isozyme II. © 1993 Wiley-Liss, Inc.     

Refinement of d(GCGAAGC) hairpin structure using one- and two-bond residual dipolar couplings

The structure of the ¹³C,¹⁵N-labeled d(GCGAAGC) hairpin, as determined by NMR spectroscopy and refined using molecular dynamics with NOE-derived distances, torsion angles, and residual dipolar couplings (RDCs), is presented. Although the studied molecule is of small size, it is demonstrated that the incorporation of diminutive RDCs can significantly improve local structure determination of regions undefined by the conventional restraints. Very good correlation between the experimental and back-calculated small one- and two-bond ¹H-¹³C, ¹H-¹⁵N, ¹³C-¹³C and ¹³C-¹⁵N coupling constants has been attained. The final structures clearly show typical features of the miniloop architecture. The structure is discussed in context of the extraordinary stability of the d(GCGAAGC) hairpin, which originates from a complex interplay between the aromatic base stacking and hydrogen bonding interactions.     

The new program OPAL for molecular dynamics simulations and energy refinements of biological macromolecules

Summary A new program for molecular dynamics (MD) simulation and energy refinement of biological macromolecules, OPAL, is introduced. Combined with the supporting program TRAJEC for the analysis of MD trajectories, OPAL affords high efficiency and flexibility for work with diferent force fields, and offers a user-friendly interface and extensive trajectory analysis capabilities. Salient features are computational speeds of up to 1.5 GFlops on vector supercomputers such as the NEC SX-3, ellipsoidal boundaries to reduce the system size for studies in explicit solvents, and natural treatment of the hydrostatic pressure. Practical applications of OPAL are illustrated with MD simulations of pure water, energy minimization of the NMR structure of the mixed disulfide of a mutant E. coli glutaredoxin with glutathione in different solvent models, and MD simulations of a small protein, pheromone Er-2, using either instantaneous or time-averaged NMR restraints, or no restraints.Abbreviations D diffusion constant in cm²/s - Er-2 pheromone 2 from Euplotes raikovi - GFlop one billion floating point operations per second - Grx(C14S)-SG mixed disulfide between a mutant E. coli glutaredoxin, with Cys¹⁴ replaced by Ser, and glutathione - MD molecular dynamics - NOE nuclear Overhauser enhancement - rmsd root-mean-square deviation - density in g/cm³     

Protein structure determination using a database of interatomic distance probabilities

The accelerated pace of genomic sequencing has increased the demand for structural models of gene products. Improved quantitative methods are needed to study the many systems (e.g., macromolecular assemblies) for which data are scarce. Here, we describe a new molecular dynamics method for protein structure determination and molecular modeling. An energy function, or database potential, is derived from distributions of interatomic distances obtained from a database of known structures. X-ray crystal structures are refined by molecular dynamics with the new energy function replacing the Van der Waals potential. Compared to standard methods, this method improved the atomic positions, interatomic distances, and side-chain dihedral angles of structures randomized to mimic the early stages of refinement. The greatest enhancement in side-chain placement was observed for groups that are characteristically buried. More accurate calculated model phases will follow from improved interatomic distances. Details usually seen only in high-resolution refinements were improved, as is shown by an R-factor analysis. The improvements were greatest when refinements were carried out using X-ray data truncated at 3.5 A. The database potential should therefore be a valuable tool for determining X-ray structures, especially when only low-resolution data are available.     

Cryo‐EM model validation using independent map reconstructions

An increasing number of cryo‐electron microscopy (cryo‐EM) density maps are being generated with suitable resolution to trace the protein backbone and guide sidechain placement. Generating and evaluating atomic models based on such maps would be greatly facilitated by independent validation metrics for assessing the fit of the models to the data. We describe such a metric based on the fit of atomic models with independent test maps from single particle reconstructions not used in model refinement. The metric provides a means to determine the proper balance between the fit to the density and model energy and stereochemistry during refinement, and is likely to be useful in determining values of model building and refinement metaparameters quite generally.     

Good student questions in inquiry learning

Acquisition of scientific reasoning is one of the big challenges in education. A popular educational strategy advocated for acquiring deep knowledge is inquiry-based learning, which is driven by emerging ‘good questions’. This study will address the question: ‘Which design features allow learners to refine questions while preserving student ownership of the inquiry process?’ This design-based research has been conducted over several years with advanced high-school biology classes. The results confirm the central role of question elaboration as an interactive process that leads from vague to complex and adequate. To make this happen, the inquiry process must extend over a long time, learners and teachers should share a knowledge improvement goal, and text produced by students should be structured by question–answer pairs addressing a single concept using authentic resources. Further features are discussion with peers, teacher feedback with respect to answer elaboration and conceptual differentiation, and finally, teacher guidance that should fade out in successive inquiry cycles to ensure student responsibility.