Analysis of Hydrogen Bonding and Stability of Protein Secondary Structures in Molecular Dynamics Simulation

Molecular Dynamics (MD) simulations provide an atomic level account of the molecular motions and have proven to be immensely useful in the investigation of the dynamical structure of proteins. Once an MD trajectory is obtained, specific interactions at the molecular level can be directly studied by setting up appropriate combinations of distance and angle monitors. However, if a study of the dynamical behavior of secondary structures in proteins becomes important, this approach can become unwieldy. We present herein a method to study the dynamical stability of secondary structures in proteins, based on a relatively simple analysis of backbone hydrogen bonds. The method was developed for studying the thermal unfolding of beta-lactamases, but can be extended to other systems and adapted to study relevant properties.

Conjugacy invariant pseudo-norms, representability and RNA secondary structures

To a reasonable approximation, a secondary structures of RNA is determined by Watson-Crick pairing without pseudo-knots in such a way as to minimise the number of unpaired bases: We show that this minimal number is determined by the maximal conjugacy-invariant pseudo-norm on the free group on two generators subject to bounds on the generators. This allows us to construct lower bounds on the minimal number of unpaired bases by constructing conjugacy invariant pseudo-norms. We show that one such construction, based on isometric actions on metric spaces, gives a sharp lower bound. A major goal here is to formulate a purely mathematical question, based on considering orthogonal representations, which we believe is of some interest independent of its biological roots.

Role of pH controlled DNA secondary structures in the reversible dispersion/precipitation and separation of metallic and semiconducting single-walled carbon nanotubes

Single-stranded DNA (ss-DNA) oligomers (dA(20), d(C(3)TA(2))(3)C-3] or dT(20)) are able to disperse single-walled carbon nanotubes (SWNTs) in water at pH 7 through non-covalent wrapping on the nanotube surface. At lower pH, an alteration of the DNA secondary structure leads to precipitation of the SWNTs from the dispersion. The structural change of dA(20) takes place from the single-stranded to the A-motif form at pH 3.5 while in case of d(C(3)TA(2))(3)C-3] the change occurs from the single-stranded to the i-motif form at pH 5. Due to this structural change, the DNA is no longer able to bind the nanotube and hence the SWNT precipitates from its well-dispersed state. However, this could be reversed on restoring the pH to 7, where the DNA again relaxes in the single-stranded form. In this way the dispersion and precipitation process could be repeated over and over again. Variable temperature UV-Vis-NIR and CD spectroscopy studies showed that the DNA-SWNT complexes were thermally stable even at similar to 90 degrees C at pH 7. Broadband NIR laser (1064 nm) irradiation also demonstrated the stability of the DNA-SWNT complex against local heating introduced through excitation of the carbon nanotubes. Electrophoretic mobility shift assay confirmed the formation of a stable DNA-SWNT complex at pH 7 and also the generation of DNA secondary structures (A/i-motif) upon acidification. The interactions of ss-DNA with SWNTs cause debundling of the nanotubes from its assembly. Selective affinity of the semiconducting SWNTs towards DNA than the metallic ones enables separation of the two as evident from spectroscopic as well as electrical conductivity studies.

Identification of local variations within secondary structures of proteins

Secondary-structure elements (SSEs) play an important role in the folding of proteins. Identification of SSEs in proteins is a common problem in structural biology. A new method, ASSP (Assignment of Secondary Structure in Proteins), using only the path traversed by the C atoms has been developed. The algorithm is based on the premise that the protein structure can be divided into continuous or uniform stretches, which can be defined in terms of helical parameters, and depending on their values the stretches can be classified into different SSEs, namely -helices, 3(10)-helices, -helices, extended -strands and polyproline II (PPII) and other left-handed helices. The methodology was validated using an unbiased clustering of these parameters for a protein data set consisting of 1008 protein chains, which suggested that there are seven well defined clusters associated with different SSEs. Apart from -helices and extended -strands, 3(10)-helices and -helices were also found to occur in substantial numbers. ASSP was able to discriminate non--helical segments from flanking -helices, which were often identified as part of -helices by other algorithms. ASSP can also lead to the identification of novel SSEs. It is believed that ASSP could provide a better understanding of the finer nuances of protein secondary structure and could make an important contribution to the better understanding of comparatively less frequently occurring structural motifs. At the same time, it can contribute to the identification of novel SSEs. A standalone version of the program for the Linux as well as the Windows operating systems is freely downloadable and a web-server version is also available at .

Evolution of conserved secondary structures and their function in transcriptional regulation networks

Background: Many conserved secondary structures have been identified within conserved elements in the human genome, but only a small fraction of them are known to be functional RNAs. The evolutionary variations of these conserved secondary structures in h

Secondary structures of mRNAs and splicing of introns in eukaryotes.

After analyzing the secondary structures of 68 exon-intron-exon and the corresponding exon-exon sequence segments, it is found that about 90% of 5' and 3' terminal bases G (splicing sites) of introns are situated in the loops of secondary structures or at the ends of stems near the loops, and most of "G" s in loops are closed to the ends of loops. Approximately 92% of the connecting sites of the adjoining exons also show the similar features. About 82% of the branch point "A" s are situated in loops or at the ends of stems near the loops. Splicing sites and branch points approach each other in space because of the folding.

The rapid assignment of ruminal fungi to presumptive genera using ITS1 and ITS2 RNA secondary structures to produce group-specific fingerprints

Danny S. Tuckwell, Matthew J. Nicholson, Christopher S. McSweeney, Michael K. Theodorou and Jayne L. Brookman (2005). The rapid assignment of ruminal fungi to presumptive genera using ITS1 and ITS2 RNA secondary structures to produce group-specific fingerprints. Microbiology, 151 (5) pp.1557-1567 Sponsorship: BBSRC / Stapledon Memorial Trust RAE2008

Modeling protein evolution using secondary structures

L’évolution des protéines est un domaine important de la recherche en bioinformatique et catalyse l'intérêt de trouver des outils d'alignement qui peuvent être utilisés de manière fiable et modéliser avec précision l'évolution d'une famille de protéines. TM-Align (Zhang and Skolnick, 2005) est considéré comme l'outil idéal pour une telle tâche, en termes de rapidité et de précision. Par conséquent, dans cette étude, TM-Align a été utilisé comme point de référence pour faciliter la détection des autres outils d'alignement qui sont en mesure de préciser l'évolution des protéines. En parallèle, nous avons élargi l'actuel outil d'exploration de structures secondaires de protéines, Helix Explorer (Marrakchi, 2006), afin qu'il puisse également être utilisé comme un outil pour la modélisation de l'évolution des protéines.

Origin and evolution of 3'UTR of flaviviruses: long direct repeats as a basis for the formation of secondary structures and their significance for virus transmission

The 3' untranslated regions (3'UTRs) of flaviviruses are reviewed and analyzed in relation to short sequences conserved as direct repeats (DRs). Previously, alignments of the 3'UTRs have been constructed for three of the four recognized flavivirus groups, namely mosquito-borne, tick-borne, and nonclassified flaviviruses (MBFV, TBFV, and NCFV, respectively). This revealed (1) six long repeat sequences (LRSs) in the 3'UTR and open-reading frame (ORF) of the TBFV, (2) duplication of the 3'UTR of the NCFV by intramolecular recombination, and (3) the possibility of a common origin for all DRs within the MBFV. We have now extended this analysis and review it in the context of all previous published analyses. This has been achieved by constructing a robust alignment between all flaviviruses using the published DRs and secondary RNA structures as "anchors" to reveal additional homologies along the 3'UTR. This approach identified nucleotide regions within the MBFV, NKV (no-known vector viruses), and NCFV 3'UTRs that are homologous to different LRSs in the TBFV 3'UTR and ORF. The analysis revealed that some of the DRs and secondary RNA structures described individually within each flavivirus group share common evolutionary origins. The 3'UTR of flaviviruses, and possibly the ORF, therefore probably evolved through multiple duplication of an RNA domain, homologous to the LRS previously identified only in the TBFV. The short DRs in all virus groups appear to represent the evolutionary remnants of these domains rather than resulting from new duplications. The relevance of these flavivirus DRs to evolution, diversity, 3'UTR enhancer function, and virus transmission is reviewed.

Mining characteristic relations bind to RNA secondary structures

The identification of RNA secondary structures has been among the most exciting recent developments in biology and medical science. It has been recognized that there is an abundance of functional structures with frameshifting, regulation of translation, and splicing functions. However, the inherent signal for secondary structures is weak and generally not straightforward due to complex interleaving substrings. This makes it difficult to explore their potential functions from various structure data. Our approach, based on a collection of predicted RNA secondary structures, allows us to efficiently capture interesting characteristic relations in RNA and bring out the top-ranked rules for specified association groups. Our results not only point to a number of interesting associations and include a brief biological interpretation to them. It assists biologists in sorting out the most significant characteristic structure patterns and predicting structurefunction relationships in RNA.

Investigating the effect of different positioning of lysine residues along the peptide chain of mastoparans for their secondary structures and biological activities

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

De Novo Design of secondary structures: Synthesis and conformational studies

Chemists have long sought to extrapolate the power of biological catalysis and recognition to synthetic systems. These efforts have focused largely on low molecular weight catalysts and receptors; however, biological systems themselves rely almost exclusively on polymers, proteins and RNA, to perform complex chemical functions. Proteins and RNA are unique in their ability to adopt compact, well-ordered conformations, and specific folding provides precise spatial orientation of the functional groups that comprise the “active site”. These features suggest that identification of new polymer backbones with discrete and predictable folding propensities (“foldamers”) will provide a basis for design of molecular machines with unique capabilities. The foldamer approach complements current efforts to design unnatural properties into polypeptides and polynucleotides. The aim of this thesis is the synthesis and conformational studies of new classes of foldamers, using a peptidomimetic approach. Moreover their attitude to be utilized as ionophores, catalysts, and nanobiomaterials were analyzed in solution and in the solid state. This thesis is divided in thematically chapters that are reported below. It begins with a very general introduction (page 4) which is useful, but not strictly necessary, to the expert reader. It is worth mentioning that paragraph I.3 (page 22) is the starting point of this work and paragraph I.5 (page 32) isrequired to better understand the results of chapters 4 and 5. In chapter 1 (page 39) is reported the synthesis and conformational analysis of a novel class of foldamers containing (S)-β3-homophenylglycine [(S)-β3-hPhg] and D- 4-carboxy-oxazolidin-2-one (D-Oxd) residues in alternate order is reported. The experimental conformational analysis performed in solution by IR, 1HNMR, and CD spectroscopy unambiguously proved that these oligomers fold into ordered structures with increasing sequence length. Theoretical calculations employing ab initio MO theory suggest a helix with 11-membered hydrogenbonded rings as the preferred secondary structure type. The novel structures enrich the field of peptidic foldamers and might be useful in the mimicry of native peptides. In chapter 2 cyclo-(L-Ala-D-Oxd)3 and cyclo-(L-Ala-DOxd) 4 were prepared in the liquid phase with good overall yields and were utilized for bivalent ions chelation (Ca2+, Mg2+, Cu2+, Zn2+ and Hg2+); their chelation skill was analyzed with ESI-MS, CD and 1HNMR techniques and the best results were obtained with cyclo-(L-Ala-D-Oxd)3 and Mg2+ or Ca2+. Chapter 3 describes an application of oligopeptides as catalysts for aldol reactions. Paragraph 3.1 concerns the use of prolinamides as catalysts of the cross aldol addition of hydroxyacetone to aromatic aldeydes, whereas paragraphs 3.2 and 3.3 are about the catalyzed aldol addition of acetone to isatins. By means of DFT and AIM calculations, the steric and stereoelectronic effects that control the enantioselectivity in the cross-aldol addition of acetone to isatin catalysed by L-proline have been studied, also in the presence of small quantities of water. In chapter 4 is reported the synthesis and the analysis of a new fiber-like material, obtained from the selfaggregation of the dipeptide Boc-L-Phe-D-Oxd-OBn, which spontaneously forms uniform fibers consisting of parallel infinite linear chains arising from singleintermolecular N-H···O=C hydrogen bonds. This is the absolute borderline case of a parallel β-sheet structure. Longer oligomers of the same series with general formula Boc-(L-Phe-D-Oxd)n-OBn (where n = 2-5), are described in chapter 5. Their properties in solution and in the solid state were analyzed, in correlation with their attitude to form intramolecular hydrogen bond. In chapter 6 is reported the synthesis of imidazolidin-2- one-4-carboxylate and (tetrahydro)-pyrimidin-2-one-5- carboxylate, via an efficient modification of the Hofmann rearrangement. The reaction affords the desired compounds from protected asparagine or glutamine in good to high yield, using PhI(OAc)2 as source of iodine(III).

Comparative analyses of the secondary structures of synthetic and intracellular yeast MFA2 mRNAs

The overall folded (global) structure of mRNA may be critical to translation and turnover control mechanisms, but it has received little experimental attention. Presented here is a comparative analysis of the basic features of the global secondary structure of a synthetic mRNA and the same intracellular eukaryotic mRNA by dimethyl sulfate (DMS) structure probing. Synthetic MFA2 mRNA of Saccharomyces cerevisiae first was examined by using both enzymes and chemical reagents to determine single-stranded and hybridized regions; RNAs with and without a poly(A) tail were compared. A folding pattern was obtained with the aid of the mfold program package that identified the model that best satisfied the probing data. A long-range structural interaction involving the 5′ and 3′ untranslated regions and causing a juxtaposition of the ends of the RNA, was examined further by a useful technique involving oligo(dT)-cellulose chromatography and antisense oligonucleotides. DMS chemical probing of A and C nucleotides of intracellular MFA2 mRNA was then done. The modification data support a very similar intracellular structure. When low reactivity of A and C residues is found in the synthetic RNA, ≈70% of the same sites are relatively more resistant to DMS modification in vivo. A slightly higher sensitivity to DMS is found in vivo for some of the A and C nucleotides predicted to be hybridized from the synthetic structural model. With this small mRNA, the translation process and mRNA-binding proteins do not block DMS modifications, and all A and C nucleotides are modified the same or more strongly than with the synthetic RNA.

Function of RNA secondary structures in transcriptional attenuation of the Bacillus subtilis pyr operon

The Bacillus subtilis pyr operon is regulated by exogenous pyrimidines by a transcriptional attenuation mechanism. Transcription in vitro from pyr DNA templates specifying attenuation regions yielded terminated and read-through transcripts of the expected lengths. Addition of the PyrR regulatory protein plus UMP led to greatly increased termination. Synthetic antisense deoxyoligonucleotides were used to probe possible secondary structures in the pyr mRNA that were proposed to play roles in controlling attenuation. Oligonucleotides predicted to disrupt terminator structures suppressed termination, whereas oligonucleotides predicted to disrupt the stem of antiterminator stem-loops strongly promoted termination at the usual termination site. Oligonucleotides that disrupt a previously unrecognized stem-loop structure, called the anti-antiterminator, the formation of which interferes with formation of the downstream antiterminator, suppressed termination. We propose that transcriptional attenuation of the pyr operon is governed by switching between alternative antiterminator versus anti-antiterminator plus terminator structures, and that PyrR acts by UMP-dependent binding to and stabilization of the anti-antiterminator.