Coiled coils at the edge of configurational heterogeneity. Structural analyses of parallel and antiparallel homotetrameric coiled coils reveal configurational sensitivity to a single solvent-exposed amino acid substitution

A detailed understanding of the mechanisms by which particular amino acid sequences can give rise to more than one folded structure, such as for proteins that undergo large conformational changes or misfolding, is a long-standing objective of protein chemistry. Here, we describe the crystal structures of a single coiled-coil peptide in distinct parallel and antiparallel tetrameric configurations and further describe the parallel or antiparallel crystal structures of several related peptide sequences; the antiparallel tetrameric assemblies represent the first crystal structures of GCN4-derived peptides exhibiting such a configuration. Intriguingly, substitution of a single solvent-exposed residue enabled the parallel coiled-coil tetramer GCN4-pLI to populate the antiparallel configuration, suggesting that the two configurations are close enough in energy for subtle sequence changes to have important structural consequences. We present a structural analysis of the small changes to the helix register and side-chain conformations that accommodate the two configurations and have supplemented these results using solution studies and a molecular dynamics energetic analysis using a replica exchange methodology. Considering the previous examples of structural nonspecificity in coiled-coil peptides, the findings reported here not only emphasize the predisposition of the coiled-coil motif to adopt multiple configurations but also call attention to the associated risk that observed crytstal structures may not represent the only (or even the major) species present in solution.     

Polynucleotides. 13. Stoichiometric and thermodynamic studies of polynucleotide helices with non-complementary residues

Stoichiometry and thermodynamic properties of polyadenylate–polyuridylate helices containing varying proportions of near-randomly distributed non-complementary uridine residues were charactrized from an analysis of their mixing curves and melting profiles measured at 259 nm and at appropriate longer wavelength isochromic points. The noncomplementary residues in this homopolymer–copolymer system (in which the homopolymer has the capacity to readjust with respect to the residues with which it is in opposition) show absolute preference for an extrahelical conformation even when situated in … AAUAA … sequences and must occur therefore as single loops. As the frequency of extrahelical residues in creases, the electrostatic energy of these complexes becomes greater, and is particularly severe for the three-stranded helices. Thus, an adenyl-ate-uridylate copolymer containing 35.2 mole percent uridine residues does not form a three-stranded complex with polyuridylate even in 0.7M Na⁺at O°C. The imperfections introduced into the helix lattice by extrahelical residues decrease the cooperativity of thermal denaturation as well as T_m. However, for the helices with extrahelical residues in low frequency (～1 per helix turn) only small increases in concentration of charge-neutralizing ions are required to bring T_m to the level of their perfect analogs. Two-stranded helices with a higher density of extra helical residues (～5 per helix turn) show [Na⁺] dependence of T_m characteristic of perfect three-stranded helices. These findings together with the absence of an effect of these imperfections on the hypochromicity per base-pair suggest only minimal disruption of helix continuity or distortion of stacking interactions that normally in volve the base pairs or triplets.     

Chain flexibility and configurational dimensions of left handed Z- DNA and right handed B-DNA helices

The configurational behaviour of flexible helices of right handed B- and left handed Z-types have been analysed using statistical mechanical procedures. The configuration-dependent parameter, most importantly, the persistence length has been computed, using the heminucleotide scheme of treating polynucleotide chains under the approximation that perturbations in the backbone torsions produce sufficient flexibility in these helices. The values of persistence lengths obtained for Z-helices are very much higher than that of B-helices indicating that former is less flexible compared to the latter. These are in accordance with the results obtained recently on B- and Z-forms of poly(dG-dC). (dG-dC) using light scattering studies. Also the persistence lengths of B_II-DNA helices characterised by a skew 3'-hemiucleotide (ε-270°), and also when they coexist with B-DNA have been computed and the values lie within the range of experimentally reported values on B-helices. It is argued that the decrease in the persistence length values of B-DNA at higher salt concentration is due to additional small fluctuations in sugar residue torsions induced due to neutralisation of electrostatic repulsions between adjacent phosphates of the nucleotide. Noteworthy is that these are correlated to winding angle variations and the consequent bending of the helix. Contribution No. 659.     

Effect of O(2′)-methylation on the stability of polynucleotide helices

The effect of ribose(O2′)-methylation on the stability of (O2′)-methylated polyribonucleotide helices has been studied by conformational energy calculations. The preferred orientation of the methyl group is found to further stabilize the helical phosphodiester conformation (g^?,g^?) due to the enhanced short-range interactions arising between the methyl groups and the adjacent ribose moieties. The experimentally observed increase in melting temperature of (O2′)-alkylated polyribonucleotides is thus attributable to the enhanced stability of the helical backbone conformation.     

Polynucleotides. XIV. Photochemical evidence for an extrahelical, solvent-accessible environment of non-complementary residues in polynucleotide helices

Ultraviolet photohydration and photodimerization reactions have been used to characterize the local environment of non-complementary uridine residues that are either isolated or clustered within the copolymer strand of poly(A,U) · poly(U) double helices. The initial photoreaction kinetics, measured from absorbance changes, are independent of the wavelength of measurement, and contain a fast-reacting component, absent in the analogous perfect helix, which corresponds quantitatively to the non-complementary residues. These results are consistent with an unstacked conformation of all non-complementary residues in these complexes. The initial kinetics demonstrate further that the bases are extrahelical and fully accessible to solvent, since their rate of photohydration is the same as for UMP. Photodimerization of these extrahelical U residues does not appear to occur within clusters of less than five residues, suggesting that neighboring extrahelical bases within smaller loops are stereochemically constrained from cyclobutane bond formation. Modification of extrahelical residues does not significantly disrupt adjacent helical secondary structure. Hence, the stoichiometric proportions of U residues that are intrahelical and extrahelical could be confirmed from the relation of hypochromicity to percentage photomodification. The resistance to photomodification of helical segments within the homopolymer-copolymer complexes (in sufficiently high counterion) is equivalent to that of the perfect analogue; the local environment of helical U residues within helices-with-loops and the perfect helix must, therefore, be essentially the same.     

A computer method for finding common base paired helices in aligned sequences: application to the analysis of random sequences.

We describe a new computer program that identifies conserved secondary structures in aligned nucleotide sequences of related single-stranded RNAs. The program employs a series of hash tables to identify and sort common base paired helices that are located in identical positions in more than one sequence. The program gives information on the total number of base paired helices that are conserved between related sequences and provides detailed information about common helices that have a minimum of one or more compensating base changes. The program is useful in the analysis of large biological sequences. We have used it to examine the number and type of complementary segments (potential base paired helices) that can be found in common among related random sequences similar in base composition to 16S rRNA from Escherichia coli. Two types of random sequences were analyzed. One set consisted of sequences that were independent but they had the same mononucleotide composition as the 16S rRNA. The second set contained sequences that were 80% similar to one another. Different results were obtained in the analysis of these two types of random sequences. When 5 sequences that were 80% similar to one another were analyzed, significant numbers of potential helices with two or more independent base changes were observed. When 5 independent sequences were analyzed, no potential helices were found in common. The results of the analyses with random sequences were compared with the number and type of helices found in the phylogenetic model of the secondary structure of 16S ribosomal RNA. Many more helices are conserved among the ribosomal sequences than are found in common among similar random sequences. In addition, conserved helices in the 16S rRNAs are, on the average, longer than the complementary segments that are found in comparable random sequences. The significance of these results and their application in the analysis of long non-ribosomal nucleotide sequences is discussed.     

The application of polyelectrolyte limiting laws to the helix–coil transition of DNA. VI. The numerical value of the axial phosphate spacing for the coil form

It is concluded on the basis of comparison of polyelectrolyte theory with published data that the mean phosphate spacing b along the contour axis of an unfolded polynucleotide single strand is in the range 3–4 Å (polyelectrolyte parameter ξ ≈ 2), regardless of temperature, base composition, or extent of stacking. This result is consistent with the low-angle X-ray scattering measurements of Gulik, Inoue, and Luzzati on poly(C). No conclusion may be drawn from this value of b concerning the structure of the chain skeleton or the spatial arrangement of the bases other than that the chain is far from an all-trans local conformation (for which b would be about 6–7 Å, the length of a nucleotide unit). The structural implications, or lack thereof, are discussed in detail in the following paper.     

Homodimeric mitochondrial phosphate transport protein. Transient subunit/subunit contact site between the transport relevant transmembrane helices A

The three Cys of the yeast (Saccharomyces cerevisiae) mitochondrial phosphate transport protein (PTP) subunit were replaced with Ser. The seven mutants (single, double, and complete Cys replacements) were expressed in yeast, and the homodimeric mutant PTPs were purified from the mitochondria and reconstituted. The pH gradient-dependent net phosphate (Pi) transport uptake rates (initial conditions: 1 mM [Pi]e, pHe 6.80; 0 mM [Pi]i, pHi 8.07) catalyzed by these reconstituted mutants are similar to those of the wild-type protein and range from 15 to 80 micromol Pi/min mg PTP protein. Aerobic media inhibit only the Pi uptake rates catalyzed by PTPs with the conserved (yeast and bovine) Cys28. This inhibition in the proteoliposomes is 84-95% and can be completely reversed by dithiothreitol. Transport by the wild type as well as by all mutant proteins with Cys28 is more than 90% inhibited by mersalyl. Transport catalyzed by mutant proteins with only Cys300 or only Cys134 is less sensitive, and that catalyzed by the no Cys mutant shows 40% inhibition by mersalyl. When dithiothreitol is removed from purified single Cys mutant proteins, only the mutant protein with Cys28 appears as a homodimer in a nonreducing SDS polyacrylamide gel. Thus, the function relevant transmembrane helix A, with Cys 28 about equidistant from the two inner membrane surfaces, is in close contact with parts of transmembrane helix A of the other subunit in the functional homodimeric PTP. The results identify for the first time not only a transmembrane helix contact site between the two subunits of a homodimeric mitochondrial transport protein but also a contact site that if locked into position blocks transport. The results are related to two available secondary transporter structures (lactose permease, glycerol-3-phosphate transporter) as well as to a low resolution projection structure and a high resolution structure of monomers of inhibitor ADP/ATP carrier complexes.     

The stability of transmembrane helices: A molecular dynamics study on the isolated helices of bacteriorhodopsin

Bacteriorhodopsin (bR) continues to be a proven testing ground for the study of integral membrane proteins (IMPs). It is important to study the stability of the individual helices of bR, as they are postulated to exist as independently stable transmembrane helices (TMHs) and also for their utility as templates for modeling other IMPs with the postulated seven-helix bundle topology. Toward this purpose, the seven helices of bR have been studied by molecular dynamics simulation in this study. The suitability of using the backbone-dependent rotamer library of side-chain conformations arrived at from the data base of globular protein structures in the case TMHs has been tested by another set of 7 helix simulations with the side-chain orientations taken from this library. The influence of the residue's net charge on the helix stability was examined by simulating the helices III, IV, and VI (from both of the above sets of helices) with zero net charge on the side chains. The results of these 20 simulations demonstrate in general the stability of the isolated helices of bR in conformity with the two-stage hypothesis of IMP folding. However, the helices I, II, V, and VII are more stable than the other three helices. The helical nature of certain regions of III, IV, and VI are influenced by factors such as the net charge and orientation of several residues. It is seen that the residues Arg, Lys, Asp, and Glu (charged residues), and Ser, Thr, Gly, and Pro, play a crucial role in the stability of the helices of bR. The backbone-dependent rotamer library for the side chains is found to be suitable for the study of TMHs in IMP. © 1996 John Wiley & Sons, Inc.     

Mutual interactions of the phosphate groups in locally deformed backbones of various DNA double helices at high salt concentrations

Changes in the free energy of mutual phosphate group interactions are calculated that accompany bending of the A-, B- and Z-DNA backbones in 0.7, 2.1 and 4.2 mol/l NaCl aqueous solutions. The bending is often found to be favoured in the direction of the double helix grooves; B-DNA prefers bending into the major groove while minor groove is the preferred bending direction of A-DNA in the presence of 0.7 mol/l NaCl. Interestingly, the preferences are reversed in 4.2 mol/l NaCl. Further stabilization of A-DNA and B-DNA backbones is achieved in some cases if bending is combined with suitable local double helix twist alterations. Bending tendencies of Z-DNA backbone are generally weaker and they decrease, in contrast to B-DNA and A-DNA, with the increasing ionic strength.     

Sequence studies of nonradioactive Mycoplasma tRNAPhe with the aid of polynucleotide phosphorylase and polynucleotide kinase

The known methods of enzymatic phosphorylation with [(32)P]phosphate of the 3'- or 5'-hydroxyl group of an oligonucleotide have been applied to oligonucleotides derived from Mycoplasma tRNA(Phe). The fingerprints obtained by both methods are very similar to each other and to that of uniformly labelled tRNA. The sequence of some oligonucleotides was determined by partial digestion of the 3'-phosphorylated fragment with spleen phosphodiesterase and of the corresponding 5'-phosphorylated fragment with venom phosphodiesterase.     

On the interpretation of the observed linear free energy relationship in phosphate hydrolysis: a thorough computational study of phosphate diester hydrolysis in solution

The hydrolysis of phosphate esters is crucially important to biological systems, being involved in, among other things, signaling, energy transduction, biosynthesis, and the regulation of protein function. Despite this, there are many questions that remain unanswered in this important field, particularly with regard to the preferred mechanism of hydrolysis of phosphate esters, which can proceed through any of multiple pathways that are either associative or dissociative in nature. Previous comparisons of calculated and observed linear free energy relationships (LFERs) for phosphate monoester dianions with different leaving groups showed that the TS character gradually changes from associative to dissociative with the increasing acidity of the leaving group, while reproducing the experimental LFER. Here, we have generated ab initio potential energy surfaces for the hydrolysis of phosphate diesters in solution, with a variety of leaving groups. Once again, the reaction changes from a compact concerted pathway to one that is more expansive in character when the acidity of the leaving group increases. When such systems are examined in solution, it is essential to take into consideration the contribution of solute to the overall activation entropy, which remains a major computational challenge. The popular method of calculating the entropy using a quasi-harmonic approximation appears to markedly overestimate the configurational entropy for systems with multiple occupied energy wells. We introduce an improved restraint release approach for evaluating configurational entropies and apply this approach to our systems. We demonstrate that when this factor is taken into account, it is possible to reproduce the experimental LFER for this system with reasonable accuracy.     

A mechanical explanation for cytoskeletal rings and helices in bacteria

Several bacterial proteins have been shown to polymerize into coils or rings on cell membranes. These include the cytoskeletal proteins MreB, FtsZ, and MinD, which together with other cell components make up what is being called the bacterial cytoskeleton. We believe that these shapes arise, at least in part, from the interaction of the inherent mechanical properties of the protein polymers and the constraints imposed by the curved cell membrane. This hypothesis, presented as a simple mechanical model, was tested with numerical energy-minimization methods from which we found that there are five low-energy polymer morphologies on a rod-shaped membrane: rings, lines, helices, loops, and polar-targeted circles. Analytic theory was used to understand the possible structures and to create phase diagrams that show which parameter combinations lead to which structures. Inverting the results, it is possible to infer the effective mechanical bending parameters of protein polymers from fluorescence images of their shapes. This theory also provides a plausible explanation for the morphological changes exhibited by the Z ring in a sporulating Bacillus subtilis; is used to calculate the mechanical force exerted on a cell membrane by a polymer; and allows predictions of polymer shapes in mutant cells.     

`Single addition' and `transnucleotidation' reactions catalyzed by polynucleotide phosphorylase. Effect of enzymatic removal of inorganic phosphate during reaction

The reaction of the tetranucleotide, pA-A(2)-A, with 2'(3')-0-(alpha-methoxyethyl)uridine 5'-diphosphate, Mg(2+) ions, and M. luteus polynucleotide phosphorylase followed by mild acid treatment to remove the blocking groups results in a 49% yield of the desired single addition product, pA-A(3)-U, together with smaller amounts of pA-A-U, pA-A-A, pA-A(2)-U, pA-A(2)-A, pA-A(3)-A, pA-A(4)-U, and pA-A(4)-A. The side products are thought to arise from the phosphorolysis of the acceptor molecule by the inorganic phosphate formed in the reaction mixture and from subsequent additions to the various oligonucleotide species by the resulting adenosine 5'-diphosphate. A system developed for the removal of inorganic phosphate as it is formed in the synthesis involves the addition to the reaction mixture of calf spleen nucleoside phosphorylase and nicotinamide riboside and, under these conditions, pA-A(3)-U can be prepared in 90% yield with essentially no side products. Under similar conditions, pA-A(3)-A, pA-A(3)-G, and pA-A(3)-C may be prepared from pA-A(2)-A and the appropriate blocked nucleoside diphosphate in yields of 85-94%. The incubation of pA-A(2)-A alone with polynucleotide phosphorylase exhibits the phenomenon of "transnucleotidation" in that the molecule is partially converted to oligonucleotides of smaller and larger chain lengths. In the presence of the phosphate removal system, however, the tetranucleotide is not attacked by the enzyme, and thus, "transnucleotidation" appears to be simply a combination of phosphorolytic and addition reactions catalyzed by trace amounts of inorganic phosphate contaminating the enzyme and/or the substrate.