Structure of guanosine-3',5'-cytidine monophosphate. I. Semi-empirical potential energy calculations and model-building

The conformation and packing scheme for guanosine-3′, 5′-cytidine monophosphate, GpC, were computed by minimizing the classical potential energy. The lowest energy conformation of the isolated molecule had dihedral angles in the range of helical RNA's and the sugar pucker was C3′ endo. This was used as the starting conformation in a packing search over orientation space, the dihedral angles being flexible in this step also. The packing search was restricted by constraints from our x-ray data, namely, (1) the dimensions of the monoclinic unit cell and its pseudo-C2 symmetry (the real space group is P2₁), (2) the location of the phosphorous atom, and (3) the orientation of the bases. In addition, a geometric function was devised to impose Watson-Crick base pairing. Thus, a trial structure could be sought without explicit inclusion of intermolecular potentials. An interactive computer graphics system was used for visualizing the calculated structures. The packing searches yielded two lowest energy schemes in which the molecules had the same conformation (similar to double-helical RNA) but different orientations within the unit cell. One of these was refined by standard x-ray methods to a discrepancy index of 14.4% in the C2 pseudocell. This served as the starting structure for the subsequent refinement in the real P2₁ cell.⁵     

Structure and function of bone collagen fibrils

The intermolecular volume of fully hydrated collagen fibrils from a number of mineralized and non-mineralized tissues of adult rats has been determined both by an exclusion technique and by a method which involves the monitoring of specific X-ray diffraction parameters. The intermolecular volume of either bone or dentinal fibrils is approximately twice that of either tail or achilles tendon, and the most frequent intermolecular distance in bone or dentine fibrils is approximately 3 Å larger than of the tendons.A number of fibrillar structures are most compatible with the intermolecular volume of rat tail tendon. These include hexagonal molecular packing and orthogonal arrays of microfibrils comprising seven parallel molecular strands. The intermolecular volume of bone or dentinal collagen fibrils, on the other hand, appears to arise from structures having a disordered or pseudo-hexagonal molecular packing, in which the most frequent intermolecular distance is about 19 Å.The space associated with collagen fibrils in adult bone is such that 70 to 80% of the mineral is located within the intermolecular space of the fibrils—approximately equal amounts of mineral being in spaces having lateral dimensions of 25 to 75 Å and 6 to 12 Å, respectively. Particles located in the latter kind of intermolecular space probably constitute, to a large extent, the non-crystalline mineral phase of adult bone.The stereo-chemical constraints on the transport of mineral ions into and within collagen fibrils of bone and tendon support the postulate that bone collagen is an in vivo catalyst for mineral deposition and further suggests that its catalytic activity may be partially regulated through its molecular packing.     

Systematic solution to homo‐oligomeric structures determined by NMR

Protein structure determination by NMR has predominantly relied on simulated annealing‐based conformational search for a converged fold using primarily distance constraints, including constraints derived from nuclear Overhauser effects, paramagnetic relaxation enhancement, and cysteine crosslinkings. Although there is no guarantee that the converged fold represents the global minimum of the conformational space, it is generally accepted that good convergence is synonymous to the global minimum. Here, we show such a criterion breaks down in the presence of large numbers of ambiguous constraints from NMR experiments on homo‐oligomeric protein complexes. A systematic evaluation of the conformational solutions that satisfy the NMR constraints of a trimeric membrane protein, DAGK, reveals 9 distinct folds, including the reported NMR and crystal structures. This result highlights the fundamental limitation of global fold determination for homo‐oligomeric proteins using ambiguous distance constraints and provides a systematic solution for exhaustive enumeration of all satisfying solutions. Proteins 2015; 83:651–661. © 2015 Wiley Periodicals, Inc.     

A fast coarse filtering method for peptide identification by mass spectrometry

MOTIVATION: We reformulate the problem of comparing mass-spectra by mapping spectra to a vector space model. Our search method leverages a metric space indexing algorithm to produce an initial candidate set, which can be followed by any fine ranking scheme. RESULTS: We consider three distance measures integrated into a multi-vantage point index structure. Of these, a semi-metric fuzzy-cosine distance using peptide precursor mass constraints performs the best. The index acts as a coarse, lossless filter with respect to the SEQUEST and ProFound scoring schemes, reducing the number of distance computations and returned candidates for fine filtering to about 0.5% and 0.02% of the database respectively. The fuzzy cosine distance term improves specificity over a peptide precursor mass filter, reducing the number of returned candidates by an order of magnitude. Run time measurements suggest proportional speedups in overall search times. Using an implementation of ProFound's Bayesian score as an example of a fine filter on a test set of Escherichia coli protein fragmentation spectra, the top results of our sample system are consistent with that of SEQUEST.     

Improving Protein Structure Prediction by New Strategies: Experimental Insights and the Genetic Algorithm

Three different approaches to improve tertiary fold prediction using the genetic algorithm are discussed: (i) Refinement of the search strategy, (ii) combination of prediction and experiment and (iii) inclusion of experimental data as selection criteria into the genetic algorithm. Examples from our current work are presented for refined strategies against crowding in solution space, definition of domain boundaries and secondary structure in combination with experiment, and direct incorporation of experimentally known distance constraints into the fitness function.Electronic Supplementary Material available.     

Three-dimensional crystals of CaATPase from sarcoplasmic reticulum. Symmetry and molecular packing.

Structural studies of CaATPase from sarcoplasmic reticulum have so far been restricted to low resolution due to the poor order of two-dimensional crystal forms. However, we report that three-dimensional microcrystals of detergent-solubilized CaATPase diffract to 7.2 A in x-ray powder patterns and may therefore provide an opportunity to study CaATPase structure at higher resolutions. In the present study, we have characterized the symmetry and molecular packing of negatively stained crystals by electron microscopy (em). By altering the detergent-to-lipid ratio, different sized crystals were produced, which adhere to an em grid in different orientations. Thus, we obtained micrographs of three different projections and from these determined unit cell dimensions to be 151 X 51 X 158 A and the three-dimensional space group to be C2 with an angle beta very close to 90 degrees; x-ray powder patterns of hydrated, unstained crystals yielded dimensions of 166 X 58 X 164 A. Micrographs from each of two principal projections were averaged to produce two-dimensional density maps. Based on these maps and on the previously determined low-resolution structure of CaATPase, a packing diagram for these three-dimensional crystals is presented and major intermolecular contacts are proposed.     

Transverse relaxation optimized 3D and 4D 15N/15N separated NOESY experiments of 15N labeled proteins

Distance constraints are an important complement to orientational constraints. While a high-resolution monomer structure of the ion channel forming polypeptide, gramicidin A, has been solved with 120 orientational constraints, the precise geometry of the dimer interface has not been characterized. Here, using both ¹³C and ¹⁵N labeled gramicidin A samples in hydrated phospholipid bilayers, both inter- and intramolecular distances have been measured with a recently developed simultaneous frequency and amplitude modulation (SFAM) solid-state NMR scheme. Using this approach ¹⁵N-¹³C₁ residual dipolar couplings across a hydrogen bond as small as 20 ± 2 Hz have been characterized. While such distances are on the order of 4.2 ± 0.2 Å, the spectroscopy is complicated by rapid global motion of the molecular structure about the bilayer normal and channel axis. Consequently, the nominal 40 Hz dipolar coupling is averaged depending on the orientation of the internuclear vector with respect to the motional axis. The intermolecular distance confirmed the previously described monomeric structure, while the intramolecular distance across the monomer–monomer interface defined this junction and confirmed the previous model of this interface.     

A genetic algorithm based molecular modeling technique for RNA stem-loop structures.

A new modeling technique for arriving at the three dimensional (3-D) structure of an RNA stem-loop has been developed based on a conformational search by a genetic algorithm and the following refinement by energy minimization. The genetic algorithm simultaneously optimizes a population of conformations in the predefined conformational space and generates 3-D models of RNA. The fitness function to be optimized by the algorithm has been defined to reflect the satisfaction of known conformational constraints. In addition to a term for distance constraints, the fitness function contains a term to constrain each local conformation near to a prepared template conformation. The technique has been applied to the two loops of tRNA, the anticodon loop and the T-loop, and has found good models with small root mean square deviations from the crystal structure. Slightly different models have also been found for the anticodon loop. The analysis of a collection of alternative models obtained has revealed statistical features of local variations at each base position.     

Structured motifs search.

In this paper, we describe an algorithm for the localization of structured models, i.e. sequences of (simple) motifs and distance constraints. It basically combines standard pattern matching procedures with a constraint satisfaction solver, and it has the ability, not present in similar tools, to search for partial matches. A significant feature of our approach, especially in terms of efficiency for the application context, is that the (potentially) exponentially many solutions to the considered problem are represented in compact form as a graph. Moreover, the time and space necessary to build the graph are linear in the number of occurrences of the component patterns.     

Purification, crystallisation and preliminary crystallographic studies of succinate:ubiquinone oxidoreductase from Escherichia coli.

A membrane protein complex, succinate dehydrogenase (SQR) from Escherichia coli has been purified and crystallised. This enzyme is composed of four subunits containing FAD, three iron-sulphur clusters and one haem b as prosthetic groups. The obtained crystals belong to the hexagonal space group P6(3) with the unit-cell dimensions of a=b=123.8 A and c=214.6 A. An asymmetric unit of the crystals contains one SQR monomer (M(r) 120 kDa). A data set is now available at 4.0 A resolution with 88.1% completeness and 0.106 R(merge). We have obtained a molecular replacement solution that shows sensible molecular packing, using the soluble domain of E. coli QFR (fumarate reductase) as a search model. The packing suggests that E. coli SQR is a crystallographic trimer rather than a dimer as observed for the E. coli QFR.     

Heterodimeric structure of superoxide dismutase in complex with its metallochaperone

The copper chaperone for superoxide dismutase (CCS) activates the eukaryotic antioxidant enzyme copper, zinc superoxide dismutase (SOD1). The 2.9 A resolution structure of yeast SOD1 complexed with yeast CCS (yCCS) reveals that SOD1 interacts with its metallochaperone to form a complex comprising one monomer of each protein. The heterodimer interface is remarkably similar to the SOD1 and yCCS homodimer interfaces. Striking conformational rearrangements are observed in both the chaperone and target enzyme upon complex formation, and the functionally essential C-terminal domain of yCCS is well positioned to play a key role in the metal ion transfer mechanism. This domain is linked to SOD1 by an intermolecular disulfide bond that may facilitate or regulate copper delivery.     

Structure determination of symmetric homo-oligomers by a complete search of symmetry configuration space, using NMR restraints and van der Waals packing

Structural studies of symmetric homo-oligomers provide mechanistic insights into their roles in essential biological processes, including cell signaling and cellular regulation. This paper presents a novel algorithm for homo-oligomeric structure determination, given the subunit structure, that is both complete, in that it evaluates all possible conformations, and data-driven, in that it evaluates conformations separately for consistency with experimental data and for quality of packing. Completeness ensures that the algorithm does not miss the native conformation, and being data-driven enables it to assess the structural precision possible from data alone. Our algorithm performs a branch-and-bound search in the symmetry configuration space, the space of symmetry axis parameters (positions and orientations) defining all possible C(n) homo-oligomeric complexes for a given subunit structure. It eliminates those symmetry axes inconsistent with intersubunit nuclear Overhauser effect (NOE) distance restraints and then identifies conformations representing any consistent, well-packed structure to within a user-defined similarity level. For the human phospholamban pentamer in dodecylphosphocholine micelles, using the structure of one subunit determined from a subset of the experimental NMR data, our algorithm identifies a diverse set of complex structures consistent with the nine intersubunit NOE restraints. The distribution of determined structures provides an objective characterization of structural uncertainty: backbone RMSD to the previously determined structure ranges from 1.07 to 8.85 A, and variance in backbone atomic coordinates is an average of 12.32 A(2). Incorporating vdW packing reduces structural diversity to a maximum backbone RMSD of 6.24 A and an average backbone variance of 6.80 A(2). By comparing data consistency and packing quality under different assumptions of oligomeric number, our algorithm identifies the pentamer as the most likely oligomeric state of phospholamban, demonstrating that it is possible to determine the oligomeric number directly from NMR data. Additional tests on a number of homo-oligomers, from dimer to heptamer, similarly demonstrate the power of our method to provide unbiased determination and evaluation of homo-oligomeric complex structures.     

Structure of bonito heart ferricytochrome c and some remarks on molecular interaction in its crystalline state.

The structure analysis of bonito heart ferricytochrome c was carried out at 2.8 A resolution by X-ray diffraction. The overall features of the molecule are virtually identical with those of bonito ferrocytochromes c and other cytochromes c. In the present work, the modes of molecular packing among cytochromes c were also compared by means of intermolecular distance maps. Some differences in the structures of ferro- and ferricytochrome c may exist on the surface of the molecules.     

1H NMR assignments of apo calcyclin and comparative structural analysis with calbindin D9k and S100 beta.

The homodimeric S100 protein calcyclin has been studied in the apo state by two-dimensional 1H NMR spectroscopy. Using a combination of scalar correlation and NOE experiments, sequence-specific 1H NMR assignments were obtained for all but one backbone and > 90% of the side-chain resonances. To our knowledge, the 2 x 90 residue (20 kDa) calcyclin dimer is the largest protein system for which such complete assignments have been made by purely homonuclear methods. Sequential and medium-range NOEs and slowly exchanging backbone amide protons identified directly the four helices and the short antiparallel beta-type interaction between the two binding loops that comprise each subunit of the dimer. Further analysis of NOEs enabled the unambiguous assignment of 556 intrasubunit distance constraints, 24 intrasubunit hydrogen bonding constraints, and 2 x 26 intersubunit distance constraints. The conformation of the monomer subunit was refined by distance geometry and restrained molecular dynamics calculations using the intrasubunit constraints only. Calculation of the dimer structure starting from this conformational ensemble has been reported elsewhere. The extent of structural homology among the apo calcyclin subunit, the monomer subunit of apo S100 beta, and monomeric apo calbindin D9k has been examined in detail by comparing 1H NMR chemical shifts and secondary structures. This analysis was extended to a comprehensive comparison of the three-dimensional structures of the calcyclin monomer subunit and calbindin D9k, which revealed greater similarity in the packing of their hydrophobic cores than was anticipated previously. Together, these results support the hypothesis that all members of the S100 family have similar core structures and similar modes of dimerization. Analysis of the amphiphilicity of Helix IV is used to explain why calbindin D9k is monomeric, but full-length S100 proteins form homodimers.     

Determination of the NMR solution structure of a specific DNA complex of the Myb DNA-binding domain

Summary The solution structure of a specific DNA complex of the minimum DNA-binding domain of the mouse c-Myb protein was determined by distance geometry calculations using a set of 1732 nuclear Overhauser enhancement (NOE) distance restraints. In order to determine the complex structure independent of the initial guess, we have developed two different procedures for the docking calculation using simulated annealing in four-dimensional space (4D-SA). One is a multiple-step procedure, where the protein and the DNA were first constructed independently by 4D-SA using only the individual intramolecular NOE distance restraints. Here, the initial structure of the protein was a random coil and that of the DNA was a typical B-form duplex. Then, as the starting structure for the next docking procedure, the converged protein and DNA structures were placed in random molecular orientations, separated by 50 Å. The two molecules were docked by 4D-SA utilizing all the restraints, including the additional 66 intermolecular distance restraints. The second procedure comprised a single step, in which a random-coil protein and a typical B-form DNA duplex were first placed 70 Å from each other. Then, using all the intramolecular and intermolecular NOE distance restraints, the complex structure was constructed by 4D-SA. Both procedures yielded the converged complex structures with similar quality and structural divergence, but the multiple-step procedure has much better convergence power than the single-step procedure. A model study of the two procedures was performed to confirm the structural quality, depending upon the number of intermolecular distance restraints, using the X-ray structure of the engrailed homeodomain-DNA complex.Abbreviations rmsd root-mean-square deviation - NOE nuclear Overhauser enhancement - 4D-SA simulated annealing in four-dimensional space - Myb-R2R3 repeats 2 and 3 of the DNA-binding domain of the c-Myb protein - DNA 16 Myb-specific binding DNA duplex with 16 base pairs - IHDD-C residues 3 to 59 of the C-chain of the engrailed homeodomain-DNA complex - DNA11 DNA duplex with base pairs 9 to 19 of the engrailed homeodomain-DNA complex     

Change in the crystal packing of soybean beta-amylase mutants substituted at a few surface amino acid residues

In spite of the high similarity of amino acid sequence and three-dimensional structure between soybean beta-amylase (SBA) and sweet potato beta-amylase (SPB), their quaternary structure is quite different, being a monomer in SBA and a tetramer in SPB. Because most of the differences in amino acid sequences are found in the surface region, we tested the tetramerization of SBA by examining mutations of residues located at the surface. We designed the SBA tetramer using the SPB tetramer structure as a model and calculating the change of accessible surface area (DeltaASA) for each residue in order to select sites for the mutation. Two different mutant genes encoding SB3 (D374Y/L481R/P487D) and SB4 (K462S added to SB3), were constructed for expression in Escherichia coli and the recombinant proteins were purified. They existed as a monomer in solution, but gave completely different crystals from the native SBA. The asymmetric unit of the mutants contains four molecules, while that of native SBA contains one. The interactions of the created interfaces revealed that there were more intermolecular interactions in the SB3 than in the SB4 tetramer. The substituted charged residues on the surface are involved in interactions with adjacent molecules in a different way, forming a new crystal packing pattern.     

Low-resolution structures of proteins in solution retrieved from X-ray scattering with a genetic algorithm.

Small-angle x-ray solution scattering (SAXS) is analyzed with a new method to retrieve convergent model structures that fit the scattering profiles. An arbitrary hexagonal packing of several hundred beads containing the problem object is defined. Instead of attempting to compute the Debye formula for all of the possible mass distributions, a genetic algorithm is employed that efficiently searches the configurational space and evolves best-fit bead models. Models from different runs of the algorithm have similar or identical structures. The modeling resolution is increased by reducing the bead radius together with the search space in successive cycles of refinement. The method has been tested with protein SAXS (0.001 < S < 0.06 A(-1)) calculated from x-ray crystal structures, adding noise to the profiles. The models obtained closely approach the volumes and radii of gyration of the known structures, and faithfully reproduce the dimensions and shape of each of them. This includes finding the active site cavity of lysozyme, the bilobed structure of gamma-crystallin, two domains connected by a stalk in betab2-crystallin, and the horseshoe shape of pancreatic ribonuclease inhibitor. The low-resolution solution structure of lysozyme has been directly modeled from its experimental SAXS profile (0.003 < S < 0.03 A(-1)). The model describes lysozyme size and shape to the resolution of the measurement. The method may be applied to other proteins, to the analysis of domain movements, to the comparison of solution and crystal structures, as well as to large macromolecular assemblies.     

Structural homology of spinach acyl carrier protein and Escherichia coli acyl carrier protein based on NMR data

Acyl carrier proteins (ACPs) from spinach and from Escherichia coli have been used to demonstrate the utility of proton NMR for comparison of homologous structures. The structure of E. coli ACP had been previously determined and modeled as a rapid equilibrium among multiple conformational forms (Kim and Prestegard, Biochemistry 28:8792–8797, 1989). Spinach ACP showed two slowly exchanging forms and could be manipulated into one form for structural study. Here we compare this single form to postulated multiple forms of E. coli ACP using the limited amount of NOE data available for the spinach protein. A number of long-range NOE contacts were present between homologous residues in both spinach and E. coli ACP, suggesting tertiary structural homology. To allow a more definitive structural comparison, a method was developed to use spinach ACP NOE constraints to search for regions of structural divergence from two postulated forms of E. coli ACP. The homologous regions of the two protein sequences were aligned, additional distance constraints were extracted from the E. coli structure, and these were mapped onto the spinach sequence. These distance constraints were combined with experimental NOE constraints and a distance geometry simulated annealing protocol was used to test for compatibility of the constraints. All of the experimental spinach NOE constraints could be successfully combined with the E. coli data, confirming the general hypothesis of structural homology. A better fit was obtained with one form, suggesting a preferential stabilization of that form in the spinach case. Proteins 27:131–143 © 1997 Wiley-Liss, Inc.     

Solution structures of proteins from NMR data and modeling: alternative folds for neutrophil peptide 5

The structure of neutrophil peptide 5 in solution has recently been reported (Pardi et al., 1988). The structure determination was accomplished by using a distance geometry algorithm and 107 interproton distance constraints obtained from 2D NMR data. In each of the eight independent solutions to the distance geometry equations, the overall fold of the polypeptide backbone was identical and the root mean square (rms) deviation between backbone atoms of the superimposed structures was small (approximately 2.4 A). In this paper we report additional NP-5 structures obtained by using a new structure generation algorithm: a Monte Carlo search in torsion angle space. These structures have a large rms backbone deviation from the distance geometry structures (approximately 5.0 A). The backbone topologies differ in significant respects from the distance geometry structures and from each other. Structures are found that are pseudo mirror images of part or all of the fold corresponding to that first obtained with the distance geometry procedure. For small proteins, the problem of distinguishing the correct structure among pseudo mirror images is likely to be greater than previously recognized. When a set of test distance constraints constructed from a novel Monte Carlo structure is used as input in the distance geometry algorithm, the fold of the resulting structure does not correspond to that of the target. The results also demonstrate that the previously accepted criteria (the magnitude of the rms deviation between multiple solutions of the distance geometry equations) for defining the accuracy and precision of a peptide structure generated from NMR data are inadequate. An energetic analysis of structures corresponding to the different folding topologies has been carried out. The molecular mechanics energies obtained by minimization and molecular dynamics refinement provide sufficient information to eliminate certain alternative structures. On the basis of a careful comparison of the different trial structures with the experimental data, it is concluded that the NP-5 peptide fold which was originally reported is most consistent with the data. An alternative fold corresponding to structures with low energies and small total distance violations is ruled out because for this fold predicted NOEs are not observed experimentally.