The application of 1H NMR chemical shift calculations to diastereotopic groups in proteins.

We have calculated chemical shifts for a range of diastereotopic protons in proteins (i.e. methylene protons, and the methyl groups of valine and leucine residues), using a recently optimised method for chemical shift calculation. The calculations are based on crystal structure coordinates, and have been compared with experimental stereospecific assignments. The results indicate that chemical shifts can be used to suggest stereospecific assignments with about 80% probability of being correct, in cases where both the experimental and the calculated chemical shift differences between a pair of diastereotopic protons are greater than 0.3 ppm. Inaccurate calculations are shown to be caused in most cases by differences between crystal and solution structures. Furthermore, chemical shift calculations based on NMR structures are shown to be capable of acting as a further constraint on structure, by limiting the range of side-chain conformations adopted in structures calculated from NMR data.     

Analysis of side-chain conformational distributions in neutrophil peptide-5 NMR structures

The side-chain conformations have been analyzed in the antimicrobial peptide, Neutrophil Peptide-5 (NP-5), whose structure was independently generated from nmr-derived distance constraints using a distance geometry algorithm. The side-chain and peptide dihedral angle distributions in the nmr structures were compared with those constructed from a data base of high-resolution protein crystal structures. The side-chain conformational preferences for NP-5 in solution are significantly different from those observed in the crystal structure data base. These results indicate that the side-chain conformations are quite disordered for many of the residues of NP-5. The absence of a correlation between the width of the conformational distribution and surface accessibility suggests that the disorder may be due to limitations in the structural information extracted from the nmr data rather than to molecular motion. However, it is also observed that the degree of conformational disorder is only weakly correlated with the number of nuclear Overhauser enhancements to a given side chain. Possible reasons for this are discussed. Molecular mechanics refinement of these structures did not significantly change the side-chain populations. Anomolously wide distributions are observed for rotations about the peptide bonds and the disulfide bonds in the NP-5 distance geometry structures, which are improved by the refinement. The very high degree of order observed for the central dihedral angle of the disulfide bond in the high-resolution crystal data base suggests that the rotation about this bond in proteins is determined by the local potential.     

Determination of interspin distances between spin labels attached to insulin: comparison of electron paramagnetic resonance data with the X-ray structure

A method was developed to determine the interspin distances of two or more nitroxide spin labels attached to specific sites in proteins. This method was applied to different conformations of spin-labeled insulins. The electron paramagnetic resonance (EPR) line broadening due to dipolar interaction is determined by fitting simulated EPR powder spectra to experimental data, measured at temperatures below 200 K to freeze the protein motion. The experimental spectra are composed of species with different relative nitroxide orientations and interspin distances because of the flexibility of the spin label side chain and the variety of conformational substates of proteins in frozen solution. Values for the average interspin distance and for the distance distribution width can be determined from the characteristics of the dipolar broadened line shape. The resulting interspin distances determined for crystallized insulins in the R6 and T6 structure agree nicely with structural data obtained by x-ray crystallography and by modeling of the spin-labeled samples. The EPR experiments reveal slight differences between crystal and frozen solution structures of the B-chain amino termini in the R6 and T6 states of hexameric insulins. The study of interspin distances between attached spin labels can be applied to obtain structural information on proteins under conditions where other methods like two-dimensional nuclear magnetic resonance spectroscopy or x-ray crystallography are not applicable.     

Reconciling the solution and X-ray structures of the villin headpiece helical subdomain: molecular dynamics simulations and double mutant cycles reveal a stabilizing cation-pi interaction

The 36-residue helical subdomain of the villin headpiece, HP36, is one of the smallest cooperatively folded proteins, folding on the microsecond time scale. The domain is an extraordinarily popular model system for both experimental and computational studies of protein folding. The structure of HP36 has been determined using X-ray crystallography and NMR spectroscopy, with the resulting structures exhibiting differences in helix packing, van der Waals contacts, and hydrogen bonding. It is important to determine the solution structure of HP36 with as much accuracy as possible since this structure is widely used as a reference for simulations and experiments. We complement the existing data by using all-atom molecular dynamics simulations with explicit solvent to evaluate which of the experimental models is the better representation of HP36 in solution. After simulation for 50 ns initiated with the NMR structure, we observed that the protein spontaneously adopts structures with a backbone conformation, core packing, and C-capping motif on the third helix that are more consistent with the crystal structure. We also examined hydrogen bonding and side chain packing interactions between D44 and R55 and between F47 and R55, respectively, which were observed in the crystal structure but not in the NMR-based solution structure. Simulations showed large fluctuations in the distance between D44 and R55, while the distance between F47 and R55 remained stable, suggesting the formation of a cation-pi interaction between those residues. Experimental double mutant cycles confirmed that the F47-R55 pair has a larger energetic coupling than the D44-R55 interaction. Overall, these combined experimental and computational studies show that the X-ray crystal structure is the better reference structure for HP36 in solution at neutral pH. Our analysis also shows how detailed molecular dynamics simulations combined with experimental validation can help bridge the gap between NMR and crystallographic methods.     

Two new crystal forms of the choline-binding domain of the major pneumococcal autolysin: insights into the dynamics of the active homodimer

Very little is known about the in vivo regulation of the catalytic activity of the major pneumococcal autolysin (LytA), a surface-exposed enzyme that rules the self-destruction of pneumococcal cells through degradation of their peptidoglycan backbone. Two new crystal forms of the cell wall anchoring domain of LytA were obtained, and their structures were solved and refined to 2.4A and 2.8A resolution. The domain is a homodimer with a boomerang-like shape in which the tertiary structure of each monomer is comprised by six independent beta hairpins arranged in a superhelical fashion. Choline molecules at the hydrophobic interface of consecutive hairpins maintain this unique structure. The C-terminal hairpin (last 13 residues of LytA) in the solenoid is responsible for the formation of the catalytically active homodimer. Although the general fold in the structures derived from both crystal forms is essentially the same, two different conformations of the basic homodimer are observed. Biochemical approaches have demonstrated the fundamental role of the 11 C-terminal residues in the catalytic activity of LytA. The studies reported here reveal the importance of some amino acid residues at the C terminus in the determination of the relative distance of the active dimeric form of the autolysin, which appears to be essential for the catalytic activity of this enzyme.     

Reduced representation model of protein structure prediction: statistical potential and genetic algorithms.

A reduced representation model, which has been described in previous reports, was used to predict the folded structures of proteins from their primary sequences and random starting conformations. The molecular structure of each protein has been reduced to its backbone atoms (with ideal fixed bond lengths and valence angles) and each side chain approximated by a single virtual united-atom. The coordinate variables were the backbone dihedral angles phi and psi. A statistical potential function, which included local and nonlocal interactions and was computed from known protein structures, was used in the structure minimization. A novel approach, employing the concepts of genetic algorithms, has been developed to simultaneously optimize a population of conformations. With the information of primary sequence and the radius of gyration of the crystal structure only, and starting from randomly generated initial conformations, I have been able to fold melittin, a protein of 26 residues, with high computational convergence. The computed structures have a root mean square error of 1.66 A (distance matrix error = 0.99 A) on average to the crystal structure. Similar results for avian pancreatic polypeptide inhibitor, a protein of 36 residues, are obtained. Application of the method to apamin, an 18-residue polypeptide with two disulfide bonds, shows that it folds apamin to native-like conformations with the correct disulfide bonds formed.     

Structural analysis of flexible proteins in solution by small angle X-ray scattering combined with crystallography

In the last few years, SAXS of biological materials has been rapidly evolving and promises to move structural analysis to a new level. Recent innovations in SAXS data analysis allow ab initio shape predictions of proteins in solution. Furthermore, experimental scattering data can be compared to calculated scattering curves from the growing data base of solved structures and also identify aggregation and unfolded proteins. Combining SAXS results with atomic resolution structures enables detailed characterizations in solution of mass, radius, conformations, assembly, and shape changes associated with protein folding and functions. SAXS can efficiently reveal the spatial organization of protein domains, including domains missing from or disordered in known crystal structures, and establish cofactor or substrate-induced conformational changes. For flexible domains or unstructured regions that are not amenable for study by many other structural techniques, SAXS provides a unique technology. Here, we present SAXS shape predictions for PCNA that accurately predict a trimeric ring assembly and for a full-length DNA repair glycosylase with a large unstructured region. These new results in combination with illustrative published data show how SAXS combined with high resolution crystal structures efficiently establishes architectures, assemblies, conformations, and unstructured regions for proteins and protein complexes in solution.     

Conservation of residue interactions in a family of Ca-binding proteins

In the TNC family of Ca-binding proteins (calmodulin, parvalbumin, intestinal calcium binding protein and troponin C) approximately 70 well-conserved amino acid sequences and six crystal structures are known. We find a clear correlation between residue contacts in the structures and residue conservation in the sequences: residues with strong sidechain-sidechain contacts in the three-dimenesional structure tend to be the more conserved in the sequence. This is one way to quantify the intuitive notion of the importance of sidechain interactions for maintaining protein three-dimensional structure in evolution and may usefully be taken into account in planning point mutations in protein engineering.     

High-resolution solution structure of reduced French bean plastocyanin and comparison with the crystal structure of poplar plastocyanin

The three-dimensional solution structure of reduced (CuI) plastocyanin from French bean leaves has been determined by distance geometry and restrained molecular dynamics methods using constraints obtained from 1H n.m.r. (nuclear magnetic resonance) spectroscopy. A total of 1244 experimental constraints were used, including 1120 distance constraints, 103 dihedral angle constraints and 21 hydrogen bond constraints. Stereospecific assignments were made for 26 methylene groups and the methyls of 11 valines. Additional constraints on copper co-ordination were included in the restrained dynamics calculations. The structures are well defined with average atomic root-mean-square deviations from the mean of 0.45 A for all backbone heavy atoms and 1.08 A for side-chain heavy atoms. French bean plastocyanin adopts a beta-sandwich structure in solution that is similar to the X-ray structure of reduced poplar plastocyanin; the average atomic root-mean-square difference between 16 n.m.r. structures and the X-ray structure is 0.76 A for all backbone heavy atoms. The conformations of the side-chains that constitute the hydrophobic core of French bean plastocyanin are very well defined. Of 47 conserved residues that populate a single chi 1 angle in solution, 43 have the same rotamer in the X-ray structure. Many surface side-chains adopt highly preferred conformations in solution, although the 3J alpha beta coupling constants often indicate some degree of conformational averaging. Some surface side-chains are disordered in both the solution and crystal structures of plastocyanin. There is a striking correlation between measures of side-chain disorder in solution and side-chain temperature factors in the X-ray structure. Side-chains that form a distinctive acidic surface region, believed to be important in binding other electron transfer proteins, appear to be disordered. Fifty backbone amide protons form hydrogen bonds to carbonyls in more than 60% of the n.m.r. structures; 45 of these amide protons exchange slowly with solvent deuterons. Ten hydrogen bonds are formed between side-chain and backbone atoms, eight of which are correlated with decreased proton exchange. Of the 60 hydrogen bonds formed in French bean plastocyanin, 56 occur in the X-ray structure of the poplar protein; two of the missing hydrogen bonds are absent as a result of mutations. It appears that molecular dynamics refinement of highly constrained n.m.r. structures allows accurate prediction of the pattern of hydrogen bonding.     

Catalytic mechanism of cyclophilin as observed in molecular dynamics simulations: pathway prediction and reconciliation of X-ray crystallographic and NMR solution data

Cyclophilins are proteins that catalyze X-proline cis-trans interconversion, where X represents any amino acid. Its mechanism of action has been investigated over the past years but still generates discussion, especially because until recently structures of the ligand in the cis and trans conformations for the same system were lacking. X-ray crystallographic structures for the complex cyclophilin A and HIV-1 capsid mutants with ligands in the cis and trans conformations suggest a mechanism where the N-terminal portion of the ligand rotates during the cis-trans isomerization. However, a few years before, a C-terminal rotating ligand was proposed to explain NMR solution data. In the present study we use molecular dynamics (MD) simulations to generate a trans structure starting from the cis structure. From simulations starting from the cis and trans structures obtained through the rotational pathways, the seeming contradiction between the two sets of experimental data could be resolved. The simulated N-terminal rotated trans structure shows good agreement with the equivalent crystal structure and, moreover, is consistent with the NMR data. These results illustrate the use of MD simulation at atomic resolution to model structural transitions and to interpret experimental data.     

Structure of crambin in solution, crystal and in the trajectories of molecular dynamics simulations

The mechanisms of the three-dimensional crambin structure alterations in the crystalline environments and in the trajectories of the molecular dynamics simulations in the vacuum and crystal surroundings have been analyzed. In the crystalline state and in the solution the partial regrouping of remote intramolecular packing contacts, involved in the formation and stabilization of the tertiary structure of the crambin molecule, occurs in NMR structures. In the crystalline state it is initiated by the formation of the intermolecular contacts, the conformational influence of its appearance is distributed over the structure. The changes of the conformations and positions of the residues of the loop segments, where the intermolecular contacts of the crystal surroundings are preferably concentrated, are most observable. Under the influence of these contacts the principal change of the regular secondary structure of crambin is taking place: extension of the two-strand β structure to the three-strand structure with the participation of the single last residue N46 of the C-terminal loop. In comparison with the C-terminal loop the more profound changes are observed in the conformation and the atomic positions of the backbone atoms and in the solvent accessibility of the residues of the interhelical loop. In the solution of the ensemble of the 8 NMR structures relative accessibility to the solvent differs more noticeably also in the region of the loop segments and rather markedly in the interhelical loop. In the crambin cryogenic crystal structures the positions of the atoms of the backbone and/or side chain of 14–18 of 46 residues are discretely disordered. The disorganizations of at least 8 of 14 residues occur directly in the regions of the intermolecular contacts and another 5 residues are disordered indirectly through the intramolecular contacts with the residues of the intermolecular contacts. Upon the molecular dynamics simulation in the vacuum surrounding as in the solution of the crystalline structure of crambin the essential changes of the backbone conformation are caused by the intermolecular contacts absence, but partly masked by the structure changes owing to the nonpolar H atoms absence on the simulated structure. The intermolecular contact absence is partly manifested upon the molecular dynamics simulation of the crambin crystal with one protein molecule. Compared to the crystal structure the lengths of the interpeptide hydrogen bonds and other interresidue contacts in an average solution NMR structure are somewhat shorter and accordingly the energy of the interpeptide hydrogen bonds is better. This length shortening can occur at the stage of the refinement of the NMR structures of the crambin and other proteins by its energy minimizations in the vacuum surroundings and not exist in the solution protein structures.     

The linker of calmodulin lacking Glu84 is elongated in solution,although it is bent in the crystal

The solution structure of a mutant calmodulin (des84) lacking Glu84 in the central helix linking the two calmodulin lobes is substantially different from its crystal structure. As determined by small-angle X-ray scattering, the radius of gyration and the maximum linear dimension of des84 in the presence of 0.1 mM calcium are 20.8 Å and 62.5 Å, respectively. These respective dimensions are larger than those expected from the crystal structure of des84, 18.5 Å and 55.0 Å, and smaller than those expected from the crystal structure of wild type, 22.8 Å and 67.5 Å. The distance distribution function of des84 indicates that it assumes an elongated, dumbbell shape in solution. The solution scattering profile of des84 is indistinguishable from that of wild-type calmodulin. The calcium-dependent binding of melittin to des84 causes a change in its shape from elongated to spherical, as seen with other calmodulins. In comparison with calcium-saturated des84, calcium-free des84 is slightly elongated; a slight compaction is observed with native calmodulin. The observed differences between the averaged solution structure and the crystal structure of des84 suggests that an ensemble of structures is available to calmodulin in solution and that its target need not induce a change in its conformation. These results support the theory that the linker region of the central helix of calmodulin functions as a flexible tether. © 1996 Wiley-Liss, Inc.     

Closed conformation of the active site loop of rabbit muscle triosephosphate isomerase in the absence of substrate: evidence of conformational heterogeneity

The active site loop of triosephosphate isomerase (TIM) exhibits a hinged-lid motion, alternating between the two well defined "open" and "closed" conformations. Until now the closed conformation had only been observed in protein complexes with substrate analogues. Here, we present the first rabbit muscle apo TIM structure, refined to 1.5A resolution, in which the active site loop is either in the open or in the closed conformation in different subunits of the enzyme. In the closed conformation described here, the lid loop residues participate in stabilizing hydrogen bonds characteristic of holo TIM structures, whereas chemical interactions observed in the open loop conformation are similar to those found in the apo structures of TIM. In the closed conformation, a number of water molecules are observed at the projected ligand atom positions that are hydrogen bonded to the active site residues. Additives used during crystallization (DMSO and Tris molecules and magnesium atoms) were modeled in the electron density maps. However, no specific binding of these molecules is observed at, or close to, the active site and the lid loop. To further investigate this unusual closed conformation of the apo enzyme, two more rabbit muscle TIM structures, one in the same and another in a different crystal form, were determined. These structures present the open lid conformation only, indicating that the closed conformation cannot be explained by crystal contact effects. To rationalize why the active site loop is closed in the absence of ligand in one of the subunits, extensive comparison with previously solved TIM structures was carried out, supported by the bulk of available experimental information about enzyme kinetics and reaction mechanism of TIM. The observation of both open and closed lid conformations in TIM crystals might be related to a persistent conformational heterogeneity of this protein in solution.     

Molecular Morphology of Eukaryotic Class-1 Translation Termination Factor eRF1 in Solution

The integral structural parameters and the shape of the molecule of human translation termination factor eRF1 were determined from the small-angle X-ray scattering in solution. The molecular shapes were found by bead modeling with nonlinear minimization of the root-mean-square deviation of the calculated from the experimental scattering curve. Comparisons of the small-angle scattering curves computed for atomic-resolution structures of eRF1 with the experimental data on scattering from solution testified that the crystal and the solution conformations are close. In the ribosome, the distance between the eRF1 motifs GGQ and NIKS must be shorter than in crystal or solution (75 versus 100–107 Å). Therefore, like its bacterial counterpart RF2, the eukaryotic eRF1 must change its conformation as it binds to the ribosome. The conformational mobility of eukaryotic and prokaryotic class-1 release factors is another feature making them functionally akin to tRNA.     

Molecular morphology of eukaryotic class I translation termination factor eRF1 in solution

The integral structural parameters and the shape of the molecule of human translation termination factor eRF1 were determined from the small-angle X-ray scattering in solution. The molecular shapes were found by bead modeling with nonlinear minimization of the root-mean-square deviation of the calculated from the experimental scattering curve. Comparisons of the small-angle scattering curves computed for atomic-resolution structures of eRF1 with the experimental data on scattering from solution testified that the crystal and the solution conformations are close. In the ribosome, the distance between the eRF1 motifs GGQ and NIKS must be shorter than in crystal or solution (75 versus 107-112 A). Therefore, like its bacterial counterpart RF2, the eukaryotic eRF1 must change its conformation as it binds to the ribosome. The conformational mobility of eukaryotic and prokaryotic class-1 release factors is another feature making them functionally akin to tRNA.     

Crystal and solution structures of the B-DNA dodecamer d(CGCAAATTTGCG) probed by Raman spectroscopy: heterogeneity in the crystal structure does not persist in the solution structure

The self-complementary dodecamer d(CGCAAATTTGCG) crystallizes as a double helix of the B form and manifests a Raman spectrum with features not observed in Raman spectra of either DNA solutions or wet DNA fibers. A number of Raman bands are assigned to specific nucleoside sugar and phosphodiester conformations associated with this model B-DNA crystal structure. The Raman bands proposed as markers of the crystalline B-DNA structure are compared and contrasted with previously proposed markers of Z-DNA and A-DNA crystals. The results indicate that the three canonical forms of DNA can be readily distinguished by Raman spectroscopy. However, unlike Z-DNA and A-DNA, which retain their characteristic Raman fingerprints in aqueous solution, the B-DNA Raman spectrum is not completely conserved between crystal and solution states. The Raman spectra reveal greater heterogeneity of nucleoside conformations (sugar puckers) in the DNA molecules of the crystal structure than in those of the solution structure. The results are consistent with conversion of one-third of the dG residues from the C2'-endo/anti conformation in the solution structure to another conformation, deduced to be C1'-exo/anti, in the crystal. The dodecamer crystal also exhibits unusually broad Raman bands at 790 and 820 cm-1, associated with the geometry of the phosphodiester backbone and indicating a wider range of (alpha, zeta) backbone torsion angles in the crystal than in the solution structure. The results suggest that backbone torsion angles in the CGC and GCG sequences, which flank the central AAATTT sequence, are significantly different for crystal and solution structures, the former containing the greater diversity.(ABSTRACT TRUNCATED AT 250 WORDS)     

Homodimeric chicken galectin CG-1B (C-14): Crystal structure and detection of unique redox-dependent shape changes involving inter- and intrasubunit disulfide bridges by gel filtration, ultracentrifugation, site-directed mutagenesis, and peptide mass fingerprinting

Intrafamily gene diversification has led to three prototype galectins in chicken [i.e., chicken galectin (CG)-1A, CG-1B, and CG-2] that show distinct expression profiles and developmental regulation. In order to pinpoint structural disparities among them, we determined the crystal structure of CG-1B. Alteration of the position of the Trp ring in the lectin site and the presence of only two ordered water molecules therein, as well as changes in the interface region between the two subunits, set the structure of CG-1B clearly apart from that of CG-1A. Intriguingly, the unique presence of two Cys residues at positions 2 and 7 in the N-terminal region translated into formation of an intersubunit disulfide bridge between the Cys7 residues of the homodimer in the crystal. In solution, oxidation is associated with significant shape changes in the dimeric protein and the additional occurrence of a compacted form with an intrasubunit disulfide bridge between Cys2 and Cys7. The single-site mutant C7S/C7V was not subjected to such changes, supporting the crucial role of Cys7 in redox-dependent shape changes. These results point to the functional significance of the distinctive presence of the two Cys residues in the N-terminal region of CG-1B.     

The solution structure of a monomeric insulin. A two-dimensional 1H-NMR study of des-(B26-B30)-insulin in combination with distance geometry and restrained molecular dynamics.

The solution conformation of des-(B26-B30)-insulin (DPI) has been investigated by 1H-NMR spectroscopy. A set of 250 approximate interproton distance restraints, derived from two-dimensional nuclear Overhauser enhancement spectra, were used as the basis of a structure determination using distance geometry (DG) and distance-bound driven dynamics (DDD). Sixteen DG structures were optimized using energy minimization (EM) and submitted to short 5-ps restrained molecular dynamics (RMD) simulations. A further refinement of the DDD structure with the lowest distance errors was done by energy minimization, a prolonged RMD simulation in vacuo and a time-averaged RMD simulation. An average structure was obtained from a trajectory generated during 20-ps RMD. The final structure was compared with the des-(B26-B30)-insulin crystal structure refined by molecular dynamics and the 2-Zn crystal structure of porcine insulin. This comparison shows that the overall structure of des-(B26-B30)-insulin is retained in solution with respect to the crystal structures with a high flexibility at the N-terminal part of the A chain and at the N-terminal and C-terminal parts of the B chain. In the RMD run a high mobility of Gly A1, Asn A21 and of the side chain of Phe B25 is noticed. One of the conformations adopted by des-(B26-B30)-insulin in solution is similar to that of molecule 1 (Chinese nomenclature) in the crystal structure of porcine insulin.     

Computer-based screening of functional conformers of proteins

A long-standing goal in biology is to establish the link between function, structure, and dynamics of proteins. Considering that protein function at the molecular level is understood by the ability of proteins to bind to other molecules, the limited structural data of proteins in association with other bio-molecules represents a major hurdle to understanding protein function at the structural level. Recent reports show that protein function can be linked to protein structure and dynamics through network centrality analysis, suggesting that the structures of proteins bound to natural ligands may be inferred computationally. In the present work, a new method is described to discriminate protein conformations relevant to the specific recognition of a ligand. The method relies on a scoring system that matches critical residues with central residues in different structures of a given protein. Central residues are the most traversed residues with the same frequency in networks derived from protein structures. We tested our method in a set of 24 different proteins and more than 260,000 structures of these in the absence of a ligand or bound to it. To illustrate the usefulness of our method in the study of the structure/dynamics/function relationship of proteins, we analyzed mutants of the yeast TATA-binding protein with impaired DNA binding. Our results indicate that critical residues for an interaction are preferentially found as central residues of protein structures in complex with a ligand. Thus, our scoring system effectively distinguishes protein conformations relevant to the function of interest.