Structural model for an oligonucleotide containing a bulged guanosine by NMR and energy minimization

We present three-dimensional structural models for a DNA oligomer containing a bulged guanosine based on proton NMR data and energy minimization computations. The nonexchangeable proton resonances of the duplex 5'd(GATGGGCAG).d(CTGCGCCATC) are assigned by nuclear Overhauser effect spectroscopy (NOESY) and correlated spectroscopy connectivities, and the NMR spectrum is compared with that of a regular 8-mer of similar sequence, 5'd(GATGGCAG).d(CTGCCATC). Experimental proton-proton distances are obtained from NOESY spectra acquired with mixing times of 100, 150, and 200 ms. A refined three-dimensional structure for the bulge-containing duplex is calculated from regular B DNA starting coordinates by using the AMBER molecular mechanics program [Weiner, S. J., Kollman, P. A., Case, D. A., Singh, U. C., Ghio, C., Alagona, G., Profeta, S., & Weiner, P. (1984) J. Am. Chem. Soc. 106, 765-784]. We compare structures obtained by building the helix in three and four base pair increments with structures obtained by direct minimization of the entire nine base sequence, with and without experimental distance constraints. The general features of all the calculated structures are very similar. The helix is of the B family, with the extra guanine stacked into the helix, and the helix axis is bent by 18-23 degrees, in agreement with gel mobility data for bulge-containing sequences [Rice, J. A. (1987) Ph.D. Thesis, Yale University].     

Computation of mixed phosphatidylcholine-cholesterol bilayer structures by energy minimization.

The energetically preferred structures of dimyristoylphosphatidylcholine (DMPC)-cholesterol bilayers were determined at a 1:1 mole ratio. Crystallographic symmetry operations were used to generate planar bilayers of cholesterol and DMPC. Energy minimization was carried out with respect to bond rotations, rigid body motions, and the two-dimensional lattice constants. The lowest energy structures had a hydrogen bond between the cholesterol hydroxyl and the carbonyl oxygen of the sn-2 acyl chain, but the largest contribution to the intermolecular energy was from the nonbonded interactions between the flat alpha surface of cholesterol and the acyl chains of DMPC. Two modes of packing in the bilayer were found; in structure A (the global minimum), unlike molecules are nearest neighbors, whereas in structure B (second lowest energy) like-like intermolecular interactions predominate. Crystallographic close packing of the molecules in the bilayer was achieved, as judged from the molecular areas and the bilayer thickness. These energy-minimized structures are consistent with the available experimental data on mixed bilayers of lecithin and cholesterol, and may be used as starting points for molecular dynamics or other calculations on bilayers.     

Accurate flexible fitting of high-resolution protein structures into cryo-electron microscopy maps using coarse-grained pseudo-energy minimization

Cryo-electron microscopy (cryo-EM) has been widely used to explore conformational states of large biomolecular assemblies. The detailed interpretation of cryo-EM data requires the flexible fitting of a known high-resolution protein structure into a low-resolution cryo-EM map. To this end, we have developed what we believe is a new method based on a two-bead-per-residue protein representation, and a modified form of the elastic network model that allows large-scale conformational changes while maintaining pseudobonds and secondary structures. Our method minimizes a pseudo-energy which linearly combines various terms of the modified elastic network model energy with a cryo-EM-fitting score and a collision energy that penalizes steric collisions. Unlike previous flexible fitting efforts using the lowest few normal modes, our method effectively utilizes all normal modes so that both global and local structural changes can be fully modeled. We have validated our method for a diverse set of 10 pairs of protein structures using simulated cryo-EM maps with a range of resolutions and in the absence/presence of random noise. We have shown that our method is both accurate and efficient compared with alternative techniques, and its performance is robust to the addition of random noise. Our method is also shown to be useful for the flexible fitting of three experimental cryo-EM maps.     

Conformation of a bulge-containing oligomer from a hot-spot sequence by NMR and energy minimization

Two-dimensional nmr data on a bulge-containing oligodeoxyribonucleotide, 5'dGATGGGCAG.dCTGACCCATC, and a regular oligomer of similar sequence, 5'dGATGGCAG.dCTGCCATC, are presented. The nonexchangeable protons are assigned from sequential nuclear Overhauser effect spectroscopy (NOESY) connectivities. The two-dimensional NOE (NOESY) and correlated (COSY) spectra of the bulge-containing oligomer are compared to those of the perfect 8-mer. Experimental proton-proton distances are determined from NOESY spectra acquired with mixing times of 100, 150, and 200 ms, using comparable distances in the B-DNA region of the molecule as a calibration. With this approach, measured distances do not depend systematically on mixing time. Energy minimization techniques are used to calculate a three-dimensional structure for the bulge-containing oligomer in agreement with the nmr data. The helix is of the B family, with the extra adenine stacked into the helix, and the helix axis is bent by 20 degrees.     

A comparison of optimal and suboptimal RNA secondary structures predicted by free energy minimization with structures determined by phylogenetic comparison.

This article describes the latest version of an RNA folding algorithm that predicts both optimal and suboptimal solutions based on free energy minimization. A number of RNA's with known structures deduced from comparative sequence analysis are folded to test program performance. The group of solutions obtained for each molecule is analysed to determine how many of the known helixes occur in the optimal solution and in the best suboptimal solution. In most cases, a structure about 80% correct is found with a free energy within 2% of the predicted lowest free energy structure.     

Global minimization of an off-lattice potential energy function using a chaperone-based refolding method

A global energy minimization method based on what is known about the mechanisms of the GroEL/GroES chaperonin system is applied to two 22-mers of an off-lattice protein model whose native states are beta-hairpins and which have structural similarity to short peptides known to interact strongly with the GroEL substrate binding domain. These model substrates have been used by other workers to test the effectiveness of a number of global minimization techniques, and are regarded as providing a significant challenge. The minimization method developed here is progressively elaborated from an initial simple form that targets exposed hydrophobic regions for unfolding to include a refolding phase that encourages the later recompactification of partly unfolded substrate; this refolding phase is seen to be crucial in the successful application of the method. The optimal handling of hydrophilic monomers within the model is also systematically explored, and it is seen that the best interpretation of their role is one that allows the chaperonin model to operate in "proofreading" mode whereby misfolded substrates are recognized by their surface exposure of a large proportion of hydrophobic monomers. The final version of the model allows native-like structures to be found quickly, on average for the two 22-mer substrates after 6 or 7 chaperone contacts. These results compare very favorably with those that have been obtained elsewhere using generic global minimization methods such as those based on thermal annealing. The paper concludes with a discussion of the place of the technique within the general category of hypersurface deformation methods for global minimization, and with suggestions as to how the chaperone-based method developed here could be elaborated so as to be effective on longer substrate chains that give rise to more complex tertiary structures in their native states.     

Protein contributions to redox potentials of homologous rubredoxins: an energy minimization study.

The energetic contributions of the protein to the redox potential in an iron-sulfur protein are studied via energy minimization, comparing homologous rubredoxins from Clostridium pasteurianum, Desulfovibrio gigas, Desulfovibrio vulgaris, and Pyrococcus furiosus. The reduction reaction was divided into 1) the change in the redox site charge without allowing the protein to respond and 2) the relaxation of the protein in response to the new charge state, focusing on the latter. The energy minimizations predict structural relaxation near the redox site that agrees well with that in crystal structures of oxidized and reduced P. furiosus rubredoxin, but underpredicts it far from the redox site. However, the relaxation energies from the energy-minimized structures agree well with those from the crystal structures, because the polar groups near the redox site are the main determinants and the charged groups are all located at the surface and thus are screened dielectrically. Relaxation energies are necessary for good agreement with experimentally observed differences in reduction energies between C. pasteurianum and the other three rubredoxins. Overall, the relaxation energy is large (over 500 mV) from both the energy-minimized and the crystal structures. In addition, the range in the relaxation energy for the different rubredoxins is large (300 mV), because even though the structural perturbations of the polar groups are small, they are very near the redox site. Thus the relaxation energy is an important factor to consider in reduction energetics.     

Translocation of a tRNA with an extended anticodon through the ribosome

Coordinated translocation of the tRNA-mRNA complex by the ribosome occurs in a precise, stepwise movement corresponding to a distance of three nucleotides along the mRNA. Frameshift suppressor tRNAs generally contain an extra nucleotide in the anticodon loop and they subvert the normal mechanisms used by the ribosome for frame maintenance. The mechanism by which suppressor tRNAs traverse the ribosome during translocation is poorly understood. Here, we demonstrate translocation of a tRNA by four nucleotides from the A site to the P site, and from the P site to the E site. We show that translocation of a punctuated mRNA is possible with an extra, unpaired nucleotide between codons. Interestingly, the NMR structure of the four nucleotide anticodon stem-loop reveals a conformation different from the canonical tRNA structure. Flexibility within the loop may allow conformational adjustment upon A site binding and for interacting with the four nucleotide codon in order to shift the mRNA reading frame.     

NMR structures of paramagnetic metalloproteins

Metalloproteins represent a large share of the proteome and many of them contain paramagnetic metal ions. The knowledge, at atomic resolution, of their structure in solution is important to understand processes in which they are involved, such as electron transfer mechanisms, enzymatic reactions, metal homeostasis and metal trafficking, as well as interactions with their partners. Formerly considered as unfeasible, the first structure in solution by nuclear magnetic resonance (NMR) of a paramagnetic protein was obtained in 1994. Methodological and instrumental advancements pursued over the last decade are such that NMR structure of paramagnetic proteins may be now routinely obtained. We focus here on approaches and problems related to the structure determination of paramagnetic proteins in solution through NMR spectroscopy. After a survey of the background theory, we show how the effects produced by the presence of a paramagnetic metal ion on the NMR parameters, which are in many cases deleterious for the detection of NMR spectra, can be overcome and turned into an additional source of structural restraints. We also briefly address features and perspectives given by the use of 13C-detected protonless NMR spectroscopy for proteins in solution. The structural information obtained through the exploitation of a paramagnetic center are discussed for some Cu2+ -binding proteins and for Ca2+ -binding proteins, where the replacement of a diamagnetic metal ion with suitable paramagnetic metal ions suggests novel approaches to the structural characterization of proteins containing diamagnetic and NMR-silent metal ions.     

Conformational studies of cyclic peptide structures in solution from 1H-Nmr data by distance geometry calculation and restrained energy minimization

The three-dimensional structure of a cyclic bouvardin analogue, cyclo (-Pro-MeTyr-Ala-MeTyr-MeTyr-D-Ala-) has been determined by distance geometry calculation and restrained energy minimization from nmr data. The preparation of the input for the distance geometry calculations, the modification of the amino acid library, and the analysis of the structures were done with the aid of a recently developed software package, GEOM. A great variety of different initial structures were explored to check the uniqueness of the determined solution structure. Calculations with 500 different initial structures and two different strategies led to a uniquely determined backbone conformation with a root mean square deviations value of 0.4 A. The backbone structure consists of two beta-turns, a beta-II turn at Pro1-MeTyr2, and a beta-VI turn at MeTyr4-MeTyr5. The efficiency of the two calculation strategies were compared in order to propose an optimal means for performing distance geometry calculations with cyclic structures.     

Conformational changes in the selectivity filter of the open-state KcsA channel: an energy minimization study

Potassium channels switch between closed and open conformations and selectively conduct K⁺ ions. There are at least two gates. The TM2 bundle at the intracellular site is the primary gate of KcsA, and rearrangements at the selectivity filter (SF) act as the second gate. The SF blocks ion flow via an inactivation process similar to C-type inactivation of voltage-gated K⁺ channels. We recently generated the open-state conformation of the KcsA channel. We found no major, possibly inactivating, structural changes in the SF associated with this massive inner-pore rearrangement, which suggests that the gates might act independently. Here we energy-minimize the open state of wild-type and mutant KcsA, validating in silico structures of energy-minimized SFs by comparison with crystallographic structures, and use these data to gain insight into how mutation, ion depletion, and K⁺ to Na⁺ substitution influence SF conformation. Both E71 or D80 protonations/mutations and the presence/absence of protein-buried water molecule(s) modify the H-bonding network stabilizing the P-loops, spawning numerous SF conformations. We find that the inactivated state corresponds to conformations with a partially unoccupied or an entirely empty SF. These structures, involving modifications in all four P-loops, are stabilized by H-bonds between amide H and carbonyl O atoms from adjacent P-loops, which block ion passage. The inner portions of the P-loops are more rigid than the outer parts. Changes are localized to the outer binding sites, with innermost site S4 persisting in the inactivated state. Strong binding by Na⁺ locally contracts the SF around Na⁺, releasing ligands that do not participate in Na⁺ coordination, and occluding the permeation pathway. K⁺ selectivity primarily appears to arise from the inability of the SF to completely dehydrate Na⁺ ions due to basic structural differences between liquid water and the “quasi-liquid” SF matrix.     

Multibilayer structure of dimyristoylphosphatidylcholine dihydrate as determined by energy minimization

Complete energy minimization was carried out on the multibilayer crystal structure of 1,2-dimyristoyl-sn-glycero-3-phosphocholine dihydrate (DMPC.2H2O), starting from the X-ray structure determination reported by Pearson and Pascher (1979) Nature 281, 499-501. The asymmetric unit contains two nonidentical DMPC molecules and four water molecules. Minimization removed the acyl chain disorder present in the X-ray structure and caused the carbon planes of the acyl chains to become mutually parallel. Two energy-minimized structures (structures I and II) were found which mainly differed in the hydrogen-bonding arrangement of the waters of hydration. In structure I as in the X-ray structure, one of the water molecules forms a hydrogen-bonded bridge between successive bilayers; but in structure II, all hydrogen bonds are satisfied on the same bilayer. Structure II corresponds to the global energy minimum and is also a suitable structure for single bilayers. The lattice constants and cell volume of the minimized structures are close to the experimental values. The electrostatic force between DMPC bilayers is attractive. The mean hydration energy of the water is -14.2 kcal/mol, which is 2.5 kcal/mol lower than the binding energy of ice.     

NMR study of hybrid hemoglobins containing unnatural heme: effect of heme modification on their tertiary and quaternary structures

The effect of heme modification on the tertiary and quaternary structures of hemoglobins was examined by utilizing the NMR spectra of the reconstituted [mesohemoglobin (mesoHb), deuterohemoglobin (deuteroHb)] and hybrid heme (meso-proto, deutero-proto) hemoglobins (Hbs). The heme peripheral modification resulted in the preferential downfield shift of the proximal histidine N1H signal for the beta subunit, indicating nonequivalence of the structural change induced by the heme modification in the alpha and beta subunits of Hb. In the reconstituted and hybrid heme Hbs, the exchangeable proton resonances due to the intra- and intersubunit hydrogen bonds, which have been used as the oxy and deoxy quaternary structural probes, were shifted by 0.2-0.3 ppm from that of native Hb upon the beta-heme substitution. This suggests that, in the fully deoxygenated form, the quaternary structure of the reconstituted Hbs is in an "imperfect" T state in which the hydrogen bonds located at the subunit interface are slightly distorted by the conformational change of the beta subunit. Moreover, the two heme orientations are found in the alpha subunit of deuteroHb, but not in the beta subunit of deuteroHb, and in both the alpha and beta subunits of mesoHb. The tertiary and quaternary structural changes in the Hb molecule induced by the heme peripheral modification were also discussed in relation to their functional properties.     

Ruthenium-iron hybrid hemoglobins as a model for partially liganded hemoglobin: NMR studies of their tertiary and quaternary structures

Diruthenium-substituted Ru-Fe hybrid hemoglobins (Hb) were synthesized by heme substitution from protoheme to ruthenium (II) carbonyldeuteroporphyrin in the alpha or beta subunits. As the carbon monoxide coordinated to ruthenium (II) is not released under physiological conditions, deoxygenated Ru-Fe hybrid derivatives [alpha(Fe)2 beta(Ru-CO)2 and alpha(Ru-CO)2 beta(Fe)2] can serve as models for half-liganded Hbs. On the basis of proton NMR spectra of hyperfine-shifted proton resonances, these Ru-Fe hybrid Hbs have only small structural changes in the heme environment of the partner subunits at low pH. The proton NMR spectra of the intersubunit hydrogen-bonded protons also showed that the quaternary structures of the two complementary hybrids both remain in the "T-like state" at low pH, suggesting that the T to R structural conversion is induced by ligation of the third ligand molecule. Marked conformational changes in the heme vicinity are observed at high pH only for alpha(Ru-CO)2 beta(Fe)2, and its quaternary structure is converted into the "R state"; the alpha(Fe)2 beta(Ru-CO)2 hybrid does not undergo this change. This implies that the free-energy difference between the two quaternary states is smaller in the alpha-liganded hybrid than in the beta-liganded one.     

NMR evidence for independent domain structures in zoocin A, an antibacterial exoenzyme

NMR was used to obtain spectroscopic evidence supporting a two domain model for zoocin A in which an N-terminal catalytic domain is linked by a threonine-proline rich linker to a target recognition domain responsible for recognizing the cell wall of bacteria susceptible to the bacteriolytic action of the enzyme. When cloned and separately expressed, each domain retains the folding found in the whole enzyme. Additionally, spectroscopy suggests that the target recognition domain has a conformation typical of a soluble globular protein, while the catalytic domain aggregates at low millimolar concentrations.     

NMR structures of polytopic integral membrane proteins

Membrane proteins represent up to 30% of the proteins in all organisms, they are involved in many biological processes and are the molecular targets for around 50% of validated drugs. Despite this, membrane proteins represent less than 1% of all high-resolution protein structures due to various challenges associated with applying the main biophysical techniques used for protein structure determination. Recent years have seen an explosion in the number of high-resolution structures of membrane proteins determined by NMR spectroscopy, especially for those with multiple transmembrane-spanning segments. This is a review of the structures of polytopic integral membrane proteins determined by NMR spectroscopy up to the end of the year 2010, which includes both β-barrel and α-helical proteins from a number of different organisms and with a range in types of function. It also considers the challenges associated with performing structural studies by NMR spectroscopy on membrane proteins and how some of these have been overcome, along with its exciting potential for contributing new knowledge about the molecular mechanisms of membrane proteins, their roles in human disease, and for assisting drug design.