Helix bundles and coiled coils in alpha-spectrin and tropomyosin: a theoretical CD study

The dipole interaction model is used to investigate the effects of interactions between helices and supertwisting of helices by determining whether the predicted UV absorption and CD spectra for the three-helix bundle and coiled coil are significantly different from spectra for the single straight alpha-helix. Crystallographic data by Yan et al. for alpha-spectrin are used to construct a three-helix bundle of poly(L-alanine) modeling the protein. Backbone torsion angles represented by Fourier series are used to generate supertwisted helices and coiled coil models of poly(L-alanine) that have pitch, radius, and residue repeat similar to experimental crystallographic data on tropomyosin. Calculated CD spectra are compared with available experimental data. Theoretical spectra for the three-helix bundle and the supertwisted structures are quite similar to predictions for the straight alpha-helix of the same length with similar torsion angles, suggesting that CD is primarily dependent on the average backbone conformation and would not be a sensitive tool for distinguishing between single straight helices and closely packed or twisted alpha-helices.     

Polyproline II helix conformation in a proline-rich environment: a theoretical study

Interest centers here on whether a polyproline II helix can propagate through adjacent non-proline residues, and on shedding light on recent experimental observations suggesting the presence of significant PP(II) structure in a short alanine-based peptide with no proline in the sequence. For this purpose, we explored the formation of polyproline II helices in proline-rich peptides with the sequences Ac-(Pro)(3)-X-(Pro)(3)-Gly-Tyr-NH(2), with X = Pro (PPP), Ala (PAP), Gln (PQP), Gly (PGP), and Val (PVP), and Ac-(Pro)(3)-Ala-Ala-(Pro)(3)-Gly-Tyr-NH(2) (PAAP), by using a theoretical approach that includes a solvent effect as well as cis <--> trans isomerization of the peptide groups and puckering conformations of the pyrrolidine ring of the proline residues. Since (13)C chemical shifts have proven to be useful for identifying secondary-structure preferences in proteins and peptides, and because values of the dihedral angles (phi,psi) are the main determinants of their magnitudes, we have, therefore, computed the Boltzmann-averaged (13)C chemical shifts for the guest residues in the PXP peptide (X = Pro, Ala, Gln, Gly, and Val) with a combination of approaches, involving molecular mechanics, statistical mechanics, and quantum mechanics. In addition, an improved procedure was used to carry out the conformational searches and to compute the solvent polarization effects faster and more accurately than in previous work. The current theoretical work and additional experimental evidence show that, in short proline-rich peptides, alanine decreases the polyproline II helix content. In particular, the theoretical evidence accumulated in this work calls into question the proposal that alanine has a strong preference to adopt conformations in the polyproline II region of the Ramachandran map.     

The effect of salt concentration on DNA conformation transition: a molecular-dynamics study

We performed three 3-ns molecular dynamics simulations of d(CGCGAATTCGCG)₂ using the AMBER 8 package to determine the effect of salt concentration on DNA conformational transitions. All the simulations were started with A-DNA, with different salt concentrations, and converged with B-DNA with similar conformational parameters. However, the dynamic processes of the three MD simulations were very different. We found that the conformation transition was slow in the solution with higher salt concentration. To determine the cause of this retardation, we performed three additional 1.5-ns simulations starting with B-DNA and with the salt concentrations corresponding to the simulations mentioned above. However, astonishingly, there was no delayed conformation evolution found in any of the three simulations. Thus, our simulation conclusion is that higher salt concentrations slows the A → B conformation transition, but have no effect on the final stable structure. Figure A-DNA and B-DNA. (a) is the canonical A-DNA, and (b) is the canonical B-DNA. Looking from the central major groove     

The influence of the peptide bond on the conformation of amino acids: a theoretical and FT-IR matrix-isolation study of N-acetylproline

A combined experimental matrix-isolation FT-IR and theoretical study has been performed to investigate the conformational behavior of N-acetylproline. The conformational landscape of N-acetylproline was explored using successively higher computational methods, i.e. HF, DFT(B3LYP) and finally MP2. The exploration resulted in 10 conformations with a relative energy difference smaller than 22 kJ.mol(-1) at the HF/3-21G level of theory. These conformations led to six different conformations after DFT(B3LYP) optimizations. Further optimization at the MP2/6-31++G** level of theory resulted in the same six conformations, all of them with an energy difference smaller than 11.5kJ.mol(-1). One conformation with an intramolecular H-bond was found which was energetically the most favorable conformation. The vibrational and thermodynamical features were calculated using the DFT and MP2 methodologies. In the experimental matrix-isolation FT-IR spectrum, the most stable conformation was dominant and at least two non-H-bonded conformations could be identified. An experimental rotamerization constant between the H-bonded and the other non-H-bonded conformations was estimated and appeared to agree reasonably well with the theoretical MP2 predictions. Some new spectral features of N-acetylproline compared to proline were discovered which might be used to discriminate between the acetylated and non-acetylated form.     

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.     

The CD of two-chain coiled coils: experiments on tropomyosin and tropomyosin segments in the tyrosine/disulfide spectral region

CD experiments are reported for several coiled-coil species in the tyrosine/disulfide (approximately 250-350-nm) region. Intact noncross-linked tropomyosin (approximately 3 degrees C) shows a negative nonsymmetric band maximal at 280 nm. This spectrum is the sum over six tyrosines/chain, and has conformational significance, since it disappears on denaturation. Experiments on an excised coiled-coil segment, each of whose chains comprise residues 11-127 of the tropomyosin sequence and only one tyrosine (Y60), reveal that not all tyrosines are alike. The spectrum at 3 degrees C shows a small negative maximum at approximately 285 nm and a substantial, hitherto unknown, positive band at approximately 270 nm, the latter masked in the parent protein by the negative contribution from the other tyrosines. A noncross-linked coiled-coil segment comprising residues 142-281, in which Y60 is absent, shows no such positive band. This peculiarity of Y60 is confirmed by absorbance spectra, with the extinction coefficient of Y60 larger in benign media than the average of the other tyrosines. Intact (3 degrees C) C190 cross-linked tropomyosin is known to yield, besides tyrosine contributions, a positive maximum at approximately 300 nm. Subtracting the corresponding data for noncross-linked tropomyosin shows that the disulfide spectrum itself actually has two equal, partly resolved bands at, respectively, 250 and 280 nm. The existence of a chiral disulfide argues for a relatively rigid, perhaps strained, local coiled coil. A C190 cross-linked segment comprising residues 142-281 shows a chiral disulfide spectrum like tropomyosin's, but another segment, comprising residues 168-284, shows none; thus removal of residues 142-167 causes loss of chirality at C190, over 20 residues away. These spectra thus contain important information on the subtle local differences in coiled-coil structures.     

Role of the membrane interface on the conformation of the caveolin scaffolding domain: a CD and NMR study

Circular dichroism (CD) and NMR spectroscopy were used to study the conformational properties of two synthetic peptides, D82-R101 and D82-I109, encompassing the caveolin scaffolding domain (D82-R101), in the presence of dodecylphosphocholine (DPC) micelles. Our data show that a stable helical conformation of the caveolin scaffolding domain in a membrane mimicking system is only obtained for the peptide including the L102-I109 hydrophobic stretch, a part of the caveolin intra-membrane domain. Through chemical shift variations, an ensemble of six residues of the D82-L109 peptide, mainly located in the V(94)TKYWFYR(101) motif were found to detect the presence of phosphatidylserine solubilized in DPC micelles. Our results constitute a first step for elucidating at a residue level the conformational properties of the central region of the caveolin-1 protein.     

Different effects of trifluoroethanol and glycerol on the stability of tropomyosin helices and the head-to-tail complex

Tropomyosin (Tm) is a dimeric coiled-coil protein, composed of 284 amino acids (410 A), that forms linear homopolymers through head-to-tail interactions at low ionic strength. The head-to-tail complex involves the overlap of approximately nine N-terminal residues of one molecule with nine C-terminal residues of another Tm molecule. In this study, we investigate the influence of 2,2,2-trifluoroethanol (TFE) and glycerol on the stability of recombinant Tm fragments (ASTm1-142, Tm143-284(5OHW269)) and of the dimeric head-to-tail complex formed by the association of these two fragments. The C-terminal fragment (Tm143-284(5OHW269)) contains a 5-hydroxytryptophan (5OHW) probe at position 269 whose fluorescence is sensitive to the head-to-tail interaction and allows us to accompany titrations of Tm143-284(5OHW269) with ASTm1-142 to calculate the dissociation constant (Kd) and the interaction energy at TFE and glycerol concentrations between 0% and 15%. We observe that TFE, but not glycerol, reduces the stability of the head-to-tail complex. Thermal denaturation experiments also showed that the head-to-tail complex increases the overall conformational stability of the Tm fragments. Urea and thermal denaturation assays demonstrated that both TFE and glycerol increase the stability of the isolated N- and C-terminal fragments; however, only TFE caused a significant reduction in the cooperativity of unfolding these fragments. Our results show that these two cosolvents stabilize the structures of individual Tm fragments in different manners and that these differences may be related to their opposing effects on head-to-tail complex formation.     

The effect of sonication on the hydrocarbon chain conformation in model membrane systems: a Raman spectroscopic study.

Raman spectra are presented for egg lecithin above and below the gelliquid crystal phase transition, and several regions of the Raman spectrum are shown to be sensitive to conformational changes in the hydrocarbon chains. These regions are used to investigate the effect of sonication on the structure of egg lecithin and dipalmitoyl lecithin vesicles. Sonication of both egg lecithin above Tm, and dipalmitoyl lecithin above and below Tm produces no change in the relative population of trans and gauche isomers in any of the systems studied. Sonication does however appear to effect interchain interactions, a possible consequence of imperfect packing towards the center of the bilayers in vesicle systems.     

Protein-film voltammetry: a theoretical study of the temperature effect using square-wave voltammetry

Square-wave voltammetry of surface redox reactions is considered as an adequate model for a protein-film voltammetric setup. Here we develop a theoretical approach to analyze the effects of temperature on square-wave voltammograms. The performed simulations address the surface redox reactions featuring slow, modest and fast electron transfer. The theoretical calculations show that the temperature affects the square-wave voltammetric responses in a complex way resulting in a variety of peak shapes. Temperature effects on the phenomena known as "quasireversible maximum" and "split SW peaks" are also analyzed. The simulated results can be used to analyze the redox mechanisms and kinetic parameters of electron transfer reactions in protein-film criovoltammetry and other surface-confined redox systems. Our analysis also shows how "abnormal" features present in some square-wave voltammetric studies can easily be misinterpreted by postulating "multiple species", "stable radicals", or additional processes. Finally we provide a simple algorithm to use the "quasireversible maximum" to determine the activation energy of electron transfer reactions by surface redox systems.     

The effect of dominance on food hoarding: a game theoretical model

Many food hoarding animals live in small groups structured by rank. The presence of conspecifics in the hoarding area increases the risk of losing stored supplies. The possibility of stealing from others depends on a forager's rank in the group. Highly ranked individuals can steal from subordinates and also protect their own caches. Since storing incurs both costs and benefits, the optimal hoarding investment will differ between individuals of different rank. In a game theoretical model, we investigate how dominant and subordinate individuals should optimize their hoarding effort. Our model imagines animals that are large-scale hoarders in autumn and dependent on stored supplies for winter survival. Many examples can be found in the bird families Paridae and Corvidae, but the model can be used for any hoarding species that forage in groups. Predictions from the model are as follows: First, subordinates should store more than dominants, but in a predictable environment, this difference will decrease as the environment gets harsher. Under harsh conditions, dominants should store almost as much as subordinates and, later, spend almost as much time retrieving their own caches as subordinates. Second, if on the other hand, bad winter conditions were not expected when storing, dominants should spend more time pilfering caches from subordinates. Third, in populations that are highly dependent on stored supplies, dominants should store relatively more than in populations that are less dependent on stored supplies. Fourth, harsher environments will favor hoarding. And finally, if dominant individuals store, it implies that hoarders have a selfish recovery advantage over conspecific pilferers.     

Influence of the 2'-hydroxyl group conformation on the stability of A-form helices in RNA

The 2'-hydroxyl group plays fundamental roles in both the structure and the function of RNA, and is the major determinant of the conformational and thermodynamic differences between RNA and DNA. Here, we report a conformational analysis of 2'-OH groups of the HIV-2 TAR RNA by means of NMR scalar coupling measurements in solution. Our analysis supports the existence of a network of water molecules spanning the minor groove of an RNA A-form helix, as has been suggested on the basis of a high-resolution X-ray study of an RNA duplex. The 2'-OH protons of the lower stem nucleotides of the TAR RNA project either towards the O3' or towards the base, where the 2'-OH group can favorably participate in H-bonding interactions with a water molecule situated in the nucleotide base plane. We observe that the k(ex) rate of the 2'-OH proton with the bulk solvent anti-correlates with the base-pair stability, confirming the involvement of the 2'-OH group in a collective network of H-bonds, which requires the presence of canonical helical secondary structure. The methodology and conformational analysis presented here are broadly applicable and facilitate future studies aimed to correlate the conformation of the 2'-OH group with both the structure and the function of RNA and RNA-ligand complexes.     

Conformational studies of proteoglycans: theoretical studies on the conformation of heparin.

Various models proposed for heparin have been examined by a stereochemical approach involving contact distance criteria and potential energy calculations. The present study suggests that the model favored by Atkins and coworkers [Biochem. J. (1973) 135 , 729–733 and (1974) 143 , 251–252] is not stereochemically satisfactory. An alternative model has been proposed with N-acetyl-D -glucosamine and one of the uronides in the ⁴C₁ conformation and the other uronide (probably sulfated) in the ¹C₄ conformation. The observed variations in the tetrasaccharide periodicities (16.5–17.3 Å) in different crystalline modifications of heparin have been attributed to possible changes in the rotational angles about the interunit glycosidic bonds rather than a change in the pyranose ring conformation. The proposed model is also independent of the observed variation in the relative composition of uronic acid residues in heparin. These conclusions are in disagreement with those of earlier workers.     

The conformation and stability of ribonucleic acids: modeling base sequence effects in double stranded helices

Base sequence effects within double stranded RNA oligomers of A and Z conformations have been studied by molecular modeling using a methodological approach specifically adapted to nucleic acids. Calculations on symmetric oligomers having homonucleotide or dinucleotide repeating base sequences show that sequence changes can produce modifications in overall conformation, influence the degree of internal hydrogen bonding and strongly affect stability.     

The effect of sonication on the hydrocarbon chain conformation in model membrane systems: A Raman spectroscopic study

Raman spectra are presented for egg lecithin above and below the gel-liquid crystal phase transition, and several regions of the Raman spectrum are shown to be sensitive to conformational changes in the hydrocarbon chains. These regions are used to investigate the effect of sonication on the structure of egg lecithin and dipalmitoyl lecithin vesicles.Sonication of both egg lecithin above T_m, and dipalmitoyl lecithin above and below T_m produces no change in the relative population of trans and gauche isomers in any of the systems studied. Sonication does however appear to effect interchain interactions, a possible consequence of imperfect packing towards the center of the bilayers in vesicle systems.     

Circular dichroism study of NaSCN effect on conformation of poly-L-lysine

Effects of NaSCN, urea and KCl on alpha, beta and random conformations of poly-L-lysine (PLL) in water at room temperature were examined and compared quantitatively on the basis of the rotational strength of maximal peak by means of circular dichroism (CD) measurement. Alpha and beta helical conformational change of PLL was markedly concentration dependent in both the cases of NaSCN and urea, but not KCl. Among these salts, the distortion potency of millimolar concentrations of NaSCN on both alpha and beta conformations was undoubtedly several hundred times stronger than for the other salts, showing a slightly lesser effect on the alpha conformation as compared with that on the beta helical one, while there was no significant effect on random conformations even in maximal salt concentrations. The concentration required to alter the peptide conformation was substantially smaller for urea than for KCl, but both urea and KCl exhibited more effectiveness on alpha than beta conformation in contrast to NaSCN throughout the respective concentrations.     

Fluorescence studies of the conformation of pyrene-labeled tropomyosin: effects of F-actin and myosin subfragment 1

The fluorescence of pyrene-TM [rabbit skeletal tropomyosin (TM) labeled at Cys with N-(1-pyrenyl)maleimide] consists of monomer and excimer bands [Betcher-Lange, S., & Lehrer, S.S. (1978) J. Biol. Chem. 253, 3757-3760]; an increase in excimer fluorescence with temperature is due to a shift in equilibrium from a chain-closed state (N) to a chain-open state (X) associated with a helix pretransition [Graceffa, P., & Lehrer, S.S. (1980) J. Biol. Chem. 255, 11296-11300]. In this study, we show that the presence of appreciable excimer fluorescence at temperatures below the N----X pretransition (initial excimer) is due to perturbation of the TM chain-chain interaction by the pyrenes at Cys-190. Fluorescence and ATPase titrations indicated that the label caused a decrease in TM binding to F-actin primarily due to reduced end to end TM interactions on the actin filament. Under conditions where pyrene-TM was bound to F-actin, however, the excimer fluorescence did not increase with temperature, indicating that F-actin stabilizes tropomyosin by inhibiting the N----X transition. The binding of myosin subfragment 1 (S1) to pyrene-TM-F-actin at low ratios to actin caused time-dependent changes in fluorescence. After equilibrium was reached, the initial excimer fluorescence was markedly reduced and remained constant over the pretransition temperature range. Further stabilization of tropomyosin conformation on F-actin is therefore associated with S1 binding. Effects of the binding of S1 to the F-actin-tropomyosin thin filament on the state of tropomyosin were studied by monitoring the monomer fluorescence of pyrene-TM.(ABSTRACT TRUNCATED AT 250 WORDS)