A novel long distance base-pairing interaction in human immunodeficiency virus type 1 RNA occludes the Gag start codon

The 5'-untranslated region (5'-UTR) is the most conserved part of the HIV-1 RNA genome, and it contains regulatory motifs that mediate various steps in the viral life cycle. Previous work showed that the 5'-terminal 290 nucleotides of HIV-1 RNA adopt two mutually exclusive secondary structures, long distance interaction (LDI) and branched multiple hairpin (BMH). BMH has multiple hairpins, including the dimer initiation signal (DIS) hairpin that mediates RNA dimerization. LDI contains a long distance base-pairing interaction that occludes the DIS region. Consequently, the two conformations differ in their ability to form RNA dimers. In this study, we have presented evidence that the full-length 5'-UTR also adopts the LDI and BMH conformations. The downstream 290-352 region, including the Gag start codon, folds differently in the context of the LDI and BMH structures. These nucleotides form an extended hairpin structure in the LDI conformation, but the same sequences create a novel long distance interaction with upstream U5 sequences in the BMH conformation. The presence of this U5-AUG duplex was confirmed by computer-assisted RNA structure prediction, biochemical analyses, and a phylogenetic survey of different virus isolates. The U5-AUG duplex may influence translation of the Gag protein because it occludes the start codon of the Gag open reading frame.     

Inhibition of 5′-UTR RNA Conformational Switching in HIV-1 Using Antisense PNAs

Background

The genome of retroviruses, including HIV-1, is packaged as two homologous (+) strand RNA molecules, noncovalently associated close to their 5′-end in a region called dimer linkage structure (DLS). Retroviral HIV-1 genomic RNAs dimerize through complex interactions between dimerization initiation sites (DIS) within the (5′-UTR). Dimer formation is prevented by so calledLong Distance Interaction (LDI) conformation, whereas Branched Multiple Hairpin (BMH) conformation leads to spontaneous dimerization.

Methods and Results

We evaluated the role of SL1 (DIS), PolyA Hairpin signal and a long distance U5-AUG interaction by in-vitro dimerization, conformer assay and coupled dimerization and template-switching assays using antisense PNAs. Our data suggests evidence that PNAs targeted against SL1 produced severe inhibitory effect on dimerization and template-switching processes while PNAs targeted against U5 region do not show significant effect on dimerization and template switching, while PNAs targeted against AUG region showed strong inhibition of dimerization and template switching processes.

Conclusions

Our results demonstrate that PNA can be used successfully as an antisense to inhibit dimerization and template switching process in HIV -1 and both of the processes are closely linked to each other. Different PNA oligomers have ability of switching between two thermodynamically stable forms. PNA targeted against DIS and SL1 switch, LDI conformer to more dimerization friendly BMH form. PNAs targeted against PolyA haipin configuration did not show a significant change in dimerization and template switching process. The PNA oligomer directed against the AUG strand of U5-AUG duplex structure also showed a significant reduction in RNA dimerization as well as template- switching efficiency.The antisense PNA oligomers can be used to regulate the shift in the LDI/BMH equilibrium.     

The HIV-1 leader RNA conformational switch regulates RNA dimerization but does not regulate mRNA translation

The untranslated leader RNA is the most conserved part of the human immunodeficiency virus type I (HIV-1) genome. It contains many regulatory motifs that mediate a variety of steps in the viral life cycle. Previous work showed that the full-length leader RNA can adopt two alternative structures: a long distance interaction (LDI) and a branched multiple-hairpin (BMH) structure. The BMH structure exposes the dimer initiation site (DIS) hairpin, whereas this motif is occluded in the LDI structure. Consequently, these structures differ in their capacity to form RNA dimers in vitro. The BMH structure is dimerization-competent, due to DIS hairpin formation, but also presents the splice donor (SD) and RNA packaging (Psi) hairpins. In the LDI structure, an extended RNA packaging (Psi(E)) hairpin is folded, which includes the splice donor site and gag coding sequences. The gag initiation codon is engaged in a long distance base pairing interaction with sequences in the upstream U5 region in the BMH structure, thus forming the evolutionarily conserved U5-AUG duplex. Therefore, the LDI-BMH equilibrium may affect not only the process of RNA dimer formation but also translation initiation. In this study, we designed mutations in the 3'-terminal region of the leader RNA that alter the equilibrium of the LDI-BMH structures. The mutant leader RNAs are affected in RNA dimer formation, but not in their translation efficiency. These results indicate that the LDI-BMH status does not regulate HIV-1 RNA translation, despite the differential presentation of the gag initiation codon in both leader RNA structures.     

A riboswitch regulates RNA dimerization and packaging in human immunodeficiency virus type 1 virions

The genome of retroviruses, including human immunodeficiency virus type 1 (HIV-1), consists of two identical RNA strands that are packaged as noncovalently linked dimers. The core packaging and dimerization signals are located in the downstream part of the untranslated leader of HIV-1 RNA-the Psi and the dimerization initiation site (DIS) hairpins. The HIV-1 leader can adopt two alternative conformations that differ in the presentation of the DIS hairpin and consequently in their ability to dimerize in vitro. The branched multiple-hairpin (BMH) structure folds the poly(A) and DIS hairpins, but these domains are base paired in a long distance interaction (LDI) in the most stable LDI conformation. This LDI-BMH riboswitch regulates RNA dimerization in vitro. It was recently shown that the Psi hairpin structure is also presented differently in the LDI and BMH structures. Several detailed in vivo studies have indicated that sequences throughout the leader RNA contribute to RNA packaging, but how these diverse mutations affect the packaging mechanism is not known. We reasoned that these effects may be due to a change in the LDI-BMH equilibrium, and we therefore reanalyzed the structural effects of a large set of leader RNA mutations that were presented in three previous studies (J. L. Clever, D. Mirandar, Jr., and T. G. Parslow, J. Virol. 76:12381-12387, 2002; C. Helga-Maria, M. L. Hammarskjold, and D. Rekosh, J. Virol. 73:4127-4135, 1999; R. S. Russell, J. Hu, V. Beriault, A. J. Mouland, M. Laughrea, L. Kleiman, M. A. Wainberg, and C. Liang, J. Virol. 77:84-96, 2003). This analysis revealed a strict correlation between the status of the LDI-BMH equilibrium and RNA packaging. Furthermore, a correlation is apparent between RNA dimerization and RNA packaging, and these processes may be coordinated by the same LDI-BMH riboswitch mechanism.     

A model study of protein nascent chain and cotranslational folding using hydrophobic-polar residues

To study protein nascent chain folding during biosynthesis, we investigate the folding behavior of models of hydrophobic and polar (HP) chains at growing length using both two-dimensional square lattice model and an optimized three-dimensional 4-state discrete off-lattice model. After enumerating all possible sequences and conformations of HP heteropolymers up to length N = 18 and N = 15 in two and three-dimensional space, respectively, we examine changes in adopted structure, stability, and tolerance to single point mutation as the nascent chain grows. In both models, we find that stable model proteins have fewer folded nascent chains during growth, and often will only fold after reaching full length. For the few occasions where partial chains of stable proteins fold, these partial conformations on average are very similar to the corresponding parts of the final conformations at full length. Conversely, we find that sequences with fewer stable nascent chains and sequences with native-like folded nascent chains are more stable. In addition, these stable sequences in general can have many more point mutations and still fold into the same conformation as the wild type sequence. Our results suggest that stable proteins are less likely to be trapped in metastable conformations during biosynthesis, and are more resistant to point-mutations. Our results also imply that less stable proteins will require the assistance of chaperone and other factors during nascent chain folding. Taken together with other reported studies, it seems that cotranslational folding may not be a general mechanism of in vivo protein folding for small proteins, and in vitro folding studies are still relevant for understanding how proteins fold biologically.     

The efficiency of folding of some proteins is increased by controlled rates of translation in vivo. A hypothesis

We propose that the way in which some proteins fold is affected by the rates at which regions of their polypeptide chains are translated in vivo. Furthermore, we suggest that their gene sequences have evolved to control the rate of translational elongation such that the synthesis of defined portions of their polypeptide chains is separated temporally. We stress that many proteins are capable of folding efficiently into their native conformations without the help of differential translation rates. For these proteins the amino acid sequence does indeed contain all the information needed for the polypeptide chain to fold correctly (even in vitro, after denaturation). However, other proteins clearly do not fold efficiently into their native conformation in vitro. We argue that the efficiency of folding of these problematic proteins in vivo may be improved by controlled synthesis of the nascent polypeptide.     

Stability and conformation of the dimeric HIV-1 genomic RNA 5′UTR

During HIV-1 assembly, the viral Gag polyprotein specifically selects the dimeric RNA genome for packaging into new virions. The 5′ untranslated region (5′UTR) of the dimeric genome may adopt a conformation that is optimal for recognition by Gag. Further conformational rearrangement of the 5′UTR, promoted by the nucleocapsid (NC) domain of Gag, is predicted during virus maturation. Two 5′UTR dimer conformations, the kissing dimer (KD) and the extended dimer (ED), have been identified in vitro, which differ in the extent of intermolecular basepairing. Whether 5′UTRs from different HIV-1 strains with distinct sequences have access to the same dimer conformations has not been determined. Here, we applied fluorescence cross-correlation spectroscopy and single-molecule Förster resonance energy transfer imaging to demonstrate that 5′UTRs from two different HIV-1 subtypes form (KDs) with divergent stabilities. We further show that both 5′UTRs convert to a stable dimer in the presence of the viral NC protein, adopting a conformation consistent with extensive intermolecular contacts. These results support a unified model in which the genomes of diverse HIV-1 strains adopt an ED conformation.     

Folding of alpha(r)beta and epsilonbeta reverse turns; a nanosecond molecular dynamics simulation of the hexapeptide MSALNT and the octapeptide NMSALNTL in water.

Folding of the hexapeptide MSALNT and the octapeptide NMSALNTL were investigated using 2.8 ns molecular dynamics (MD) simulations in aqueous solution. In the simulation, the central sequence SALN of the hexapeptide folded rapidly within 200 ps into an alpha(r)beta turn conformation (type VIII conformation) and remained in this conformation for the rest of the trajectory. The sequence SALN of the octapeptide needed 2 ns to fold via epsilonbeta conformations into a similar conformation. The results join the sequences into a growing group of sequences which have a tendency to form secondary structures and thereby to direct protein folding. The structures of the reverse turn conformations were in accordance with the experimental results (Hakalehto et al., Eur J. Biochem. 250, 19-29 (1997)). The main driving force of folding seems to be the hydrophobic interaction between the side chains of Ala and Leu at the i+1 and i+2 positions of the beta-turn.     

Optimal HP configurations of proteins by combining local search with elastic net algorithm

The prediction of protein conformation from its amino-acid sequence is one of the most prominent problems in computational biology. But it is NP-hard. Here, we focus on an abstraction widely studied of this problem, the two-dimensional hydrophobic-polar protein folding problem (2D HP PFP). Mathematical optimal model of free energy of protein is established. Native conformations are often sought using stochastic sampling methods, but which are slow. The elastic net (EN) algorithm is one of fast deterministic methods as travelling salesman problem (TSP) strategies. However, it cannot be applied directly to protein folding problem, because of fundamental differences in the two types of problems. In this paper, how the 2D HP protein folding problem can be framed in terms of TSP is shown. Combination of the modified elastic net algorithm and novel local search method is adopted to solve this problem. To our knowledge, this is the first application of EN algorithm to 2D HP model. The results indicate that our approach can find more optimal conformations and is simple to implement, computationally efficient and fast.     

Failures of inverse folding and threading with gapped alignment

To calculate the tertiary structure of a protein from its amino acid sequence, the thermodynamic approach requires a potential function of sequence and conformation that has its global minimum at the native conformation for many different proteins. Here we study the behavior of such functions for the simplest model system that still has some of the features of the protein folding problem, namely two-dimensional square lattice chain configurations involving two residue types. First we show that even the given contact potential, which by definition is used to identify the folding sequences and their unique native conformations, cannot always correctly select which sequences will fold to a given structure. Second, we demonstrate that the given contact potential is not always able to favor the native alignment of a native sequence on its own native conformation over other gapped alignments of different folding sequences onto that same conformation. Because of these shortcomings, even in this simple model system in which all conformations and all native sequences are known and determined directly by the given potential, we must reexamine our expectations for empirical potentials used for inverse folding and gapped alignment on more realistic representations of proteins. © 1996 Wiley-Liss, Inc.     

Cotranslational Folding of Proteins

We suppose that folding of proteins occurs cotranslationally by the following scheme. The polypeptide chains enter the folding sites from protein translocation complexes (ribosome, translocation machinery incorporated in membranes) directionally with the N-terminus and gradually. The chain starts to fold as soon as its N-terminal residue enters the folding site from the translocation complex. The folding process accompanies the translocation of the chain to its folding site and is completed after the C-terminal residue leaves the translocation complex. Proteins fold in sequential stages, by translocation of their polypeptide into folding compartments. At each stage a particular conformation of the N-terminal part of the chain that has emerged from the translocation complex is formed. The formation of both the particular conformations of the N-terminal chain segment at each folding stage and the final native protein conformation at the last stage occurs in a time that does not exceed the duration of the fastest elongation cycle on the ribosome.     

Optimal HP configurations of proteins by combining local search with elastic net algorithm

The prediction of protein conformation from its amino-acid sequence is one of the most prominent problems in computational biology. But it is NP-hard. Here, we focus on an abstraction widely studied of this problem, the two-dimensional hydrophobic-polar protein folding problem (2D HP PFP). Mathematical optimal model of free energy of protein is established. Native conformations are often sought using stochastic sampling methods, but which are slow. The elastic net (EN) algorithm is one of fast deterministic methods as travelling salesman problem (TSP) strategies. However, it cannot be applied directly to protein folding problem, because of fundamental differences in the two types of problems. In this paper, how the 2D HP protein folding problem can be framed in terms of TSP is shown. Combination of the modified elastic net algorithm and novel local search method is adopted to solve this problem. To our knowledge, this is the first application of EN algorithm to 2D HP model. The results indicate that our approach can find more optimal conformations and is simple to implement, computationally efficient and fast.     

Positive selection dictates the choice between kinetic and thermodynamic protein folding and stability in subtilases

Subtilisin E (SbtE) is a member of the ubiquitous superfamily of serine proteases called subtilases and serves as a model for understanding propeptide-mediated protein folding mechanisms. Unlike most proteins that adopt thermodynamically stable conformations, the native state of SbtE is trapped into a kinetically stable conformation. While kinetic stability offers distinct functional advantages to the native state, the constraints that dictate the selection between kinetic and thermodynamic folding and stability remain unknown. Using highly conserved subtilases, we demonstrate that adaptive evolution of sequence dictates selection of folding pathways. Intracellular and extracellular serine proteases (ISPs and ESPs, respectively) constitute two subfamilies within the family of subtilases that have highly conserved sequences, structures, and catalytic activities. Our studies on the folding pathways of subtilisin E (SbtE), an ESP, and its homologue intracellular serine protease 1 (ISP1), an ISP, show that although topology, contact order, and hydrophobicity that drive protein folding reactions are conserved, ISP1 and SbtE fold through significantly different pathways and kinetics. While SbtE absolutely requires the propeptide to fold into a kinetically trapped conformer, ISP1 folds to a thermodynamically stable state more than 1 million times faster and independent of a propeptide. Furthermore, kinetics establish that ISP1 and SbtE fold through different intermediate states. An evolutionary analysis of folding constraints in subtilases suggests that observed differences in folding pathways may be mediated through positive selection of specific residues that map mostly onto the protein surface. Together, our results demonstrate that closely related subtilases can fold through distinct pathways and mechanisms, and suggest that fine sequence details can dictate the choice between kinetic and thermodynamic folding and stability.     

Reconstruction of the src-SH3 protein domain transition state ensemble using multiscale molecular dynamics simulations

We use an integrated computational approach to reconstruct accurately the transition state ensemble (TSE) for folding of the src-SH3 protein domain. We first identify putative TSE conformations from free energy surfaces generated by importance sampling molecular dynamics for a fully atomic, solvated model of the src-SH3 protein domain. These putative TSE conformations are then subjected to a folding analysis using a coarse-grained representation of the protein and rapid discrete molecular dynamics simulations. Those conformations that fold to the native conformation with a probability (P(fold)) of approximately 0.5, constitute the true transition state. Approximately 20% of the putative TSE structures were found to have a P(fold) near 0.5, indicating that, although correct TSE conformations are populated at the free energy barrier, there is a critical need to refine this ensemble. Our simulations indicate that the true TSE conformations are compact, with a well-defined central beta sheet, in good agreement with previous experimental and theoretical studies. A structured central beta sheet was found to be present in a number of pre-TSE conformations, however, indicating that this element, although required in the transition state, does not define it uniquely. An additional tight cluster of contacts between highly conserved residues belonging to the diverging turn and second beta-sheet of the protein emerged as being critical elements of the folding nucleus. A number of commonly used order parameters to identify the transition state for folding were investigated, with the number of native Cbeta contacts displaying the most satisfactory correlation with P(fold) values.