Application of circular dichroism to study RNA folding transitions

Examining DNA Structures with In-droplet Hydrogen/Deuterium Exchange Mass Spectrometry

Capillary vibrating sharp-edge spray ionization (cVSSI) combined with hydrogen/deuterium exchange-mass spectrometry (HDX-MS) has been utilized to characterize different solution-phase DNA conformers including DNA G-quadruplex topologies as well as triplex DNA and duplex DNA. In general, G-quadruplex DNA shows a wide range of protection of hydrogens extending from ~12% to ~21% deuterium incorporation. Additionally, the DNA sequences selected to represent parallel, antiparallel, and hybrid G-quadruplex topologies exhibit slight differences in deuterium uptake levels which appear to loosely relate to overall conformer stability. Notably, the exchange level for one of the hybrid sequence sub topologies of G-quadruplex DNA (24 TTG) is significantly different (compared with the others studied here) despite the DNA sequences being highly comparable. For the quadruplex-forming sequences, correlation analysis suggests protection of base hydrogens involved in tetrad hydrogen bonding. For duplex DNA ~19% deuterium incorporation is observed while only ~16% is observed for triplex DNA. This increased protection of hydrogens may be due to the added backbone scaffolding and Hoogsteen base pairing of the latter species. These experiments lay the groundwork for future studies aimed at determining the structural source of this protection as well as the applicability of the approach for ascertaining different oligonucleotide folds, co-existing conformations, and/or overall conformer flexibility.

Se-Glargine: Chemical Synthesis of a Basal Insulin Analogue Stabilized by an Internal Diselenide Bridge

Insulin has long provided a model for studies of protein folding and stability, enabling enhanced treatment of diabetes mellitus via analogue design. We describe the chemical synthesis of a basal insulin analogue stabilized by substitution of an internal cystine (A6-A11) by a diselenide bridge. The studies focused on insulin glargine (formulated as Lantus® and Toujeo®; Sanofi). Prepared at pH 4 in the presence of zinc ions, glargine exhibits a shifted isoelectric point due to a basic B chain extension (ArgB31 -ArgB32 ). Subcutaneous injection leads to pH-dependent precipitation of a long-lived depot. Pairwise substitution of CysA6 and CysA11 by selenocysteine was effected by solid-phase peptide synthesis; the modified A chain also contained substitution of AsnA21 by Gly, circumventing acid-catalyzed deamidation. Although chain combination of native glargine yielded negligible product, in accordance with previous synthetic studies, the pairwise selenocysteine substitution partially rescued this reaction: substantial product was obtained through repeated combination, yielding a stabilized insulin analogue. This strategy thus exploited both (a) the unique redox properties of selenocysteine in protein folding and (b) favorable packing of an internal diselenide bridge in the native state, once achieved. Such rational optimization of protein folding and stability may be generalizable to diverse disulfide-stabilized proteins of therapeutic interest.

Rapid Kinetics of Pistol Ribozyme: Insights into Limits to RNA Catalysis

Pistol ribozyme (Psr) is a distinct class of small endonucleolytic ribozymes, which are important experimental systems for defining fundamental principles of RNA catalysis and designing valuable tools in biotechnology. High-resolution structures of Psr, extensive structure-function studies, and computation support a mechanism involving one or more catalytic guanosine nucleobases acting as a general base and divalent metal ion-bound water acting as an acid to catalyze RNA 2'-O-transphosphorylation. Yet, for a wide range of pH and metal ion concentrations, the rate of Psr catalysis is too fast to measure manually and the reaction steps that limit catalysis are not well understood. Here, we use stopped-flow fluorescence spectroscopy to evaluate Psr temperature dependence, solvent H/D isotope effects, and divalent metal ion affinity and specificity unconstrained by limitations due to fast kinetics. The results show that Psr catalysis is characterized by small apparent activation enthalpy and entropy changes and minimal transition state H/D fractionation, suggesting that one or more pre-equilibrium steps rather than chemistry is rate limiting. Quantitative analyses of divalent ion dependence confirm that metal aquo ion pK_a correlates with higher rates of catalysis independent of differences in ion binding affinity. However, ambiguity regarding the rate-limiting step and similar correlation with related attributes such as ionic radius and hydration free energy complicate a definitive mechanistic interpretation. These new data provide a framework for further interrogation of Psr transition state stabilization and show how thermal instability, metal ion insolubility at optimal pH, and pre-equilibrium steps such as ion binding and folding limit the catalytic power of Psr suggesting potential strategies for further optimization.

The structural features of the ligand-free moaA riboswitch and its ion-dependent folding

Riboswitches are structural elements of mRNA involved in the regulation of gene expression by responding to specific cellular metabolites. To fulfil their regulatory function, riboswitches prefold into an active state, the so-called binding competent form, that guarantees metabolite binding and allows a consecutive refolding of the RNA. Here, we describe the folding pathway to the binding competent form as well as the ligand free structure of the moaA riboswitch of E. coli. This RNA proposedly responds to the molybdenum cofactor (Moco), a highly oxygen-sensitive metabolite, essential in the carbon and sulfur cycles of eukaryotes. K⁺- and Mg²⁺-dependent footprinting assays and spectroscopic investigations show a high degree of structure formation of this RNA already at very low ion-concentrations. Mg²⁺ facilitates additionally a general compaction of the riboswitch towards its proposed active structure. We show that this fold agrees with the earlier suggested secondary structure which included also a long-range tetraloop/tetraloop-receptor like interaction. Metal ion cleavage assays revealed specific Mg²⁺-binding pockets within the moaA riboswitch. These Mg²⁺ binding pockets are good indicators for the potential Moco binding site, since in riboswitches, Mg²⁺ was shown to be necessary to bind phosphate-carrying metabolites. The importance of the phosphate and of other functional groups of Moco is highlighted by binding assays with tetrahydrobiopterin, the reduced and oxygen-sensitive core moiety of Moco. We demonstrate that the general molecular shape of pterin by its own is insufficient for the recognition by the riboswitch.

pH-responsive nanoparticles based on POEOMA-b-PDPA block copolymers for RNA encapsulation, protection and cell delivery

The use of gene-based products, such as DNA or RNA, is increasingly being explored for various innovative therapies. However, the success of these strategies is highly dependent on the effective delivery of these biomolecules to target cells. Therefore, the development of pH-responsive nanoparticles comprises the creation of intelligent delivery systems with high therapeutic efficiency. In this work, the pH-responsiveness of the poly(2-(diisopropylamino)ethyl methacrylate)) (PDPA) block was investigated for the encapsulation and delivery of small RNAs (sRNA) to cancer cells. The pH responsiveness was dependent on the protonation profile of the tertiary amines of PDPA, which directly affected the electrostatic interactions established with RNA. Thus, block copolymers based on poly(oligo(ethylene oxide) methyl ether methacrylate) (POEOMA) and PDPA, POEOMA-b-PDPA, were synthesized by supplemental activator and reducing agent atom transfer radical polymerization (SARA ATRP). The structure of the block copolymers was characterized by size exclusion chromatography and ¹H NMR spectroscopy. The copolymers allowed effective complexation of model sRNAs and a pre-miRNA with efficiencies of about 89 % and 91 %, respectively. The characterization by dynamic light scattering revealed that these systems had sizes between 76 and 1375 nm. It was also found that the morphology of the polyplexes depended on the pH, since the preparation at a pH lower than the pKa of the copolymers resulted in spherical but polydisperse particles, while higher pH values resulted in nanoparticles with more homogeneous size, but altered morphology. Moreover, due to pH-responsiveness, it was achieved the release of RNA at pH higher than the pKa of the copolymers, while maintaining its integrity. The polyplexes also showed a high potential to protect RNA from RNases. The transfection of a lung cancer model (A549) and fibroblast cell lines showed that these polyplexes did not cause cell toxicity. In addition, the polyplexes enabled the effective transfection of the A549 cell line with pre-miRNA-29b and miRNA-29b, resulting in a decrease of expression levels of the target DNMT3B gene by approximately 51 % and 47 %, respectively. Overall, the POEOMA-b-PDPA copolymers proved to be a promising strategy for developing responsive delivery systems, that can play a critical role in some diseases, such as cancer, where pH varies between the intra and extracellular environments.

Challenges with Simulating Modified RNA: Insights into Role and Reciprocity of Experimental and Computational Approaches

RNA is critical to a broad spectrum of biological and viral processes. This functional diversity is a result of their dynamic nature; the variety of three-dimensional structures that they can fold into; and a host of post-transcriptional chemical modifications. While there are many experimental techniques to study the structural dynamics of biomolecules, molecular dynamics simulations (MDS) play a significant role in complementing experimental data and providing mechanistic insights. The accuracy of the results obtained from MDS is determined by the underlying physical models i.e., the force-fields, that steer the simulations. Though RNA force-fields have received a lot of attention in the last decade, they still lag compared to their protein counterparts. The chemical diversity imparted by the RNA modifications adds another layer of complexity to an already challenging problem. Insight into the effect of RNA modifications upon RNA folding and dynamics is lacking due to the insufficiency or absence of relevant experimental data. This review provides an overview of the state of MDS of modified RNA, focusing on the challenges in parameterization of RNA modifications as well as insights into relevant reference experiments necessary for their calibration.

Peptide Model of the Mutant Proinsulin Syndrome. I. Design and Clinical Correlation

The mutant proinsulin syndrome is a monogenic cause of diabetes mellitus due to toxic misfolding of insulin's biosynthetic precursor. Also designated mutant INS-gene induced diabetes of the young (MIDY), this syndrome defines molecular determinants of foldability in the endoplasmic reticulum (ER) of β-cells. Here, we describe a peptide model of a key proinsulin folding intermediate and variants containing representative clinical mutations; the latter perturb invariant core sites in native proinsulin (LeuB15→Pro, LeuA16→Pro, and PheB24→Ser). The studies exploited a 49-residue single-chain synthetic precursor (designated DesDi), previously shown to optimize in vitro efficiency of disulfide pairing. Parent and variant peptides contain a single disulfide bridge (cystine B19-A20) to provide a model of proinsulin's first oxidative folding intermediate. The peptides were characterized by circular dichroism and redox stability in relation to effects of the mutations on (a) in vitro foldability of the corresponding insulin analogs and (b) ER stress induced in cell culture on expression of the corresponding variant proinsulins. Striking correlations were observed between peptide biophysical properties, degree of ER stress and age of diabetes onset (neonatal or adolescent). Our findings suggest that age of onset reflects the extent to which nascent structure is destabilized in proinsulin's putative folding nucleus. We envisage that such peptide models will enable high-resolution structural studies of key folding determinants and in turn permit molecular dissection of phenotype-genotype relationships in this monogenic diabetes syndrome. Our companion study (next article in this issue) employs two-dimensional heteronuclear NMR spectroscopy to define site-specific perturbations in the variant peptides.

Tenuous Transcriptional Threshold of Human Sex Determination. I. SRY and Swyer Syndrome at the Edge of Ambiguity

Male sex determination in mammals is initiated by SRY, a Y-encoded transcription factor. The protein contains a high-mobility-group (HMG) box mediating sequence-specific DNA bending. Mutations causing XY gonadal dysgenesis (Swyer syndrome) cluster in the box and ordinarily arise de novo. Rare inherited variants lead to male development in one genetic background (the father) but not another (his sterile XY daughter). De novo and inherited mutations occur at an invariant Tyr adjoining the motif's basic tail (box position 72; Y127 in SRY). In SRY-responsive cell lines CH34 and LNCaP, de novo mutations Y127H and Y127C reduced SRY activity (as assessed by transcriptional activation of principal target gene Sox9) by 5- and 8-fold, respectively. Whereas Y127H impaired testis-specific enhancer assembly, Y127C caused accelerated proteasomal proteolysis; activity was in part rescued by proteasome inhibition. Inherited variant Y127F was better tolerated: its expression was unperturbed, and activity was reduced by only twofold, a threshold similar to other inherited variants. Biochemical studies of wild-type (WT) and variant HMG boxes demonstrated similar specific DNA affinities (within a twofold range), with only subtle differences in sharp DNA bending as probed by permutation gel electrophoresis and fluorescence resonance-energy transfer (FRET); thermodynamic stabilities of the free boxes were essentially identical. Such modest perturbations are within the range of species variation. Whereas our cell-based findings rationalize the de novo genotype-phenotype relationships, a molecular understanding of inherited mutation Y127F remains elusive. Our companion study uncovers cryptic biophysical perturbations suggesting that the para-OH group of Y127 anchors a novel water-mediated DNA clamp.

Triplex-Forming Peptide Nucleic Acids with Extended Backbones

Peptide nucleic acid (PNA) forms a triple helix with double-stranded RNA (dsRNA) stabilized by a hydrogen-bonding zipper formed by PNA's backbone amides (N-H) interacting with RNA phosphate oxygens. This hydrogen-bonding pattern is enabled by the matching ∼5.7 Å spacing (typical for A-form dsRNA) between PNA's backbone amides and RNA phosphate oxygens. We hypothesized that extending the PNA's backbone by one -CH2 - group might bring the distance between PNA amide groups closer to 7 Å, which is favourable for hydrogen bonding to the B-form dsDNA phosphate oxygens. Extension of the PNA backbone was expected to selectively stabilize PNA-DNA triplexes compared to PNA-RNA. To test this hypothesis, we synthesized triplex-forming PNAs that had the pseudopeptide backbones extended by an additional -CH2 - group in three different positions. Isothermal titration calorimetry measurements of the binding affinity of these extended PNA analogues for the matched dsDNA and dsRNA showed that, contrary to our structural reasoning, extending the PNA backbone at any position had a strong negative effect on triplex stability. Our results suggest that PNAs might have an inherent preference for A-form-like conformations when binding double-stranded nucleic acids. It appears that the original six-atom-long PNA backbone is an almost perfect fit for binding to A-form nucleic acids.

Reassessment of an Innovative Insulin Analogue Excludes Protracted Action yet Highlights the Distinction between External and Internal Diselenide Bridges

Long-acting insulin analogues represent the most prescribed class of therapeutic proteins. An innovative design strategy was recently proposed: diselenide substitution of an external disulfide bridge. This approach exploited the distinctive physicochemical properties of selenocysteine (U). Relative to wild type (WT), Se-insulin[C7U^A , C7U^B ] was reported to be protected from proteolysis by insulin-degrading enzyme (IDE), predicting prolonged activity. Because of this strategy's novelty and potential clinical importance, we sought to validate these findings and test their therapeutic utility in an animal model of diabetes mellitus. Surprisingly, the analogue did not exhibit enhanced stability, and its susceptibility to cleavage by either IDE or a canonical serine protease (glutamyl endopeptidase Glu-C) was similar to WT. Moreover, the analogue's pharmacodynamic profile in rats was not prolonged relative to a rapid-acting clinical analogue (insulin lispro). Although [C7U^A , C7U^B ] does not confer protracted action, nonetheless its comparison to internal diselenide bridges promises to provide broad biophysical insight.

"Register-shift" insulin analogs uncover constraints of proteotoxicity in protein evolution

Globular protein sequences encode not only functional structures (the native state) but also protein foldability, i.e. a conformational search that is both efficient and robustly minimizes misfolding. Studies of mutations associated with toxic misfolding have yielded insights into molecular determinants of protein foldability. Of particular interest are residues that are conserved yet dispensable in the native state. Here, we exploited the mutant proinsulin syndrome (a major cause of permanent neonatal-onset diabetes mellitus) to investigate whether toxic misfolding poses an evolutionary constraint. Our experiments focused on an invariant aromatic motif (Phe^B24-Phe^B25-Tyr^B26) with complementary roles in native self-assembly and receptor binding. A novel class of mutations provided evidence that insulin can bind to the insulin receptor (IR) in two different modes, distinguished by a "register shift" in this motif, as visualized by molecular dynamics (MD) simulations. Register-shift variants are active but defective in cellular foldability and exquisitely susceptible to fibrillation in vitro Indeed, expression of the corresponding proinsulin variant induced endoplasmic reticulum stress, a general feature of the mutant proinsulin syndrome. Although not present among vertebrate insulin and insulin-like sequences, a prototypical variant ([Gly^B24]insulin) was as potent as WT insulin in a rat model of diabetes. Although in MD simulations the shifted register of receptor engagement is compatible with the structure and allosteric reorganization of the IR-signaling complex, our results suggest that this binding mode is associated with toxic misfolding and so is disallowed in evolution. The implicit threat of proteotoxicity limits sequence variation among vertebrate insulins and insulin-like growth factors.

Investigations on the interface of nucleic acid aptamers and binding targets

Nucleic acid aptamers are single-stranded DNA or RNA of 20-100 nucleotides in length that have attracted substantial scientific interest due to their ability to specifically bind to target molecules via the formation of three-dimensional structures. Compared to traditional protein antibodies, aptamers have several advantages, such as their small size, high binding affinity, specificity, flexible structure, being chemical synthesizable and modifiable, good biocompatibility, high stability and low immunogenicity, which all contribute to their widely applications in the biomedical field. To date, much progress has been made in the study and applications of aptamers, however, detailed information on how aptamers bind to their targets is still scarce. Over the past few decades, many methods have been introduced to investigate the aptamer-target binding process, such as measuring the main kinetic or thermodynamic parameters, detecting the structural changes of the binding complexes, etc. Apart from traditional physicochemical methods, various types of molecular docking programs have been applied to simulate the aptamer-target interactions, while these simulations also have limitations. To facilitate the further research on the interactions, herein, we provide a brief review to illustrate the recent advances in the study of aptamer-target interactions. We summarize the binding targets of aptamers, such as small molecules, macromolecules, and even cells. Their binding constants (KD) are also summarized. Methods to probe the aptamer-target binding process, such as surface plasmon resonance (SPR), circular dichroism spectroscopy (CD), isothermal titration calorimetry (ITC), footprinting assay, truncation and mutation assay, nuclear magnetic resonance spectroscopy (NMR), X-ray crystallography and molecular docking simulation are indicated. The binding forces mediating the aptamer-target interactions, such as hydrogen bonding, electrostatic interaction, the hydrophobic effect, π-π stacking and van der Waals forces are summarized. The challenges and future perspectives are also discussed.

Structure-based stabilization of insulin as a therapeutic protein assembly via enhanced aromatic-aromatic interactions

Key contributions to protein structure and stability are provided by weakly polar interactions, which arise from asymmetric electronic distributions within amino acids and peptide bonds. Of particular interest are aromatic side chains whose directional π-systems commonly stabilize protein interiors and interfaces. Here, we consider aromatic-aromatic interactions within a model protein assembly: the dimer interface of insulin. Semi-classical simulations of aromatic-aromatic interactions at this interface suggested that substitution of residue TyrB26 by Trp would preserve native structure while enhancing dimerization (and hence hexamer stability). The crystal structure of a [TrpB26]insulin analog (determined as a T3Rf3 zinc hexamer at a resolution of 2.25 Å) was observed to be essentially identical to that of WT insulin. Remarkably and yet in general accordance with theoretical expectations, spectroscopic studies demonstrated a 150-fold increase in the in vitro lifetime of the variant hexamer, a critical pharmacokinetic parameter influencing design of long-acting formulations. Functional studies in diabetic rats indeed revealed prolonged action following subcutaneous injection. The potency of the TrpB26-modified analog was equal to or greater than an unmodified control. Thus, exploiting a general quantum-chemical feature of protein structure and stability, our results exemplify a mechanism-based approach to the optimization of a therapeutic protein assembly.

Solution structure of an ultra-stable single-chain insulin analog connects protein dynamics to a novel mechanism of receptor binding

Domain-minimized insulin receptors (IRs) have enabled crystallographic analysis of insulin-bound "micro-receptors." In such structures, the C-terminal segment of the insulin B chain inserts between conserved IR domains, unmasking an invariant receptor-binding surface that spans both insulin A and B chains. This "open" conformation not only rationalizes the inactivity of single-chain insulin (SCI) analogs (in which the A and B chains are directly linked), but also suggests that connecting (C) domains of sufficient length will bind the IR. Here, we report the high-resolution solution structure and dynamics of such an active SCI. The hormone's closed-to-open transition is foreshadowed by segmental flexibility in the native state as probed by heteronuclear NMR spectroscopy and multiple conformer simulations of crystallographic protomers as described in the companion article. We propose a model of the SCI's IR-bound state based on molecular-dynamics simulations of a micro-receptor complex. In this model, a loop defined by the SCI's B and C domains encircles the C-terminal segment of the IR α-subunit. This binding mode predicts a conformational transition between an ultra-stable closed state (in the free hormone) and an active open state (on receptor binding). Optimization of this switch within an ultra-stable SCI promises to circumvent insulin's complex global cold chain. The analog's biphasic activity, which serendipitously resembles current premixed formulations of soluble insulin and microcrystalline suspension, may be of particular utility in the developing world.

Rewriting nature's assembly manual for a ssRNA virus

Satellite tobacco necrosis virus (STNV) is one of the smallest viruses known. Its genome encodes only its coat protein (CP) subunit, relying on the polymerase of its helper virus TNV for replication. The genome has been shown to contain a cryptic set of dispersed assembly signals in the form of stem-loops that each present a minimal CP-binding motif AXXA in the loops. The genomic fragment encompassing nucleotides 1-127 is predicted to contain five such packaging signals (PSs). We have used mutagenesis to determine the critical assembly features in this region. These include the CP-binding motif, the relative placement of PS stem-loops, their number, and their folding propensity. CP binding has an electrostatic contribution, but assembly nucleation is dominated by the recognition of the folded PSs in the RNA fragment. Mutation to remove all AXXA motifs in PSs throughout the genome yields an RNA that is unable to assemble efficiently. In contrast, when a synthetic 127-nt fragment encompassing improved PSs is swapped onto the RNA otherwise lacking CP recognition motifs, assembly is partially restored, although the virus-like particles created are incomplete, implying that PSs outside this region are required for correct assembly. Swapping this improved region into the wild-type STNV1 sequence results in a better assembly substrate than the viral RNA, producing complete capsids and outcompeting the wild-type genome in head-to-head competition. These data confirm details of the PS-mediated assembly mechanism for STNV and identify an efficient approach for production of stable virus-like particles encapsidating nonnative RNAs or other cargoes.

An ultra-stable single-chain insulin analog resists thermal inactivation and exhibits biological signaling duration equivalent to the native protein

Thermal degradation of insulin complicates its delivery and use. Previous efforts to engineer ultra-stable analogs were confounded by prolonged cellular signaling in vivo, of unclear safety and complicating mealtime therapy. We therefore sought an ultra-stable analog whose potency and duration of action on intravenous bolus injection in diabetic rats are indistinguishable from wild-type (WT) insulin. Here, we describe the structure, function, and stability of such an analog, a 57-residue single-chain insulin (SCI) with multiple acidic substitutions. Cell-based studies revealed native-like signaling properties with negligible mitogenic activity. Its crystal structure, determined as a novel zinc-free hexamer at 2.8 Å, revealed a native insulin fold with incomplete or absent electron density in the C domain; complementary NMR studies are described in the accompanying article. The stability of the analog (ΔG_U 5.0(±0.1) kcal/mol at 25 °C) was greater than that of WT insulin (3.3(±0.1) kcal/mol). On gentle agitation, the SCI retained full activity for >140 days at 45 °C and >48 h at 75 °C. These findings indicate that marked resistance to thermal inactivation in vitro is compatible with native duration of activity in vivo Further, whereas WT insulin forms large and heterogeneous aggregates above the standard 0.6 mm pharmaceutical strength, perturbing the pharmacokinetic properties of concentrated formulations, dynamic light scattering, and size-exclusion chromatography revealed only limited SCI self-assembly and aggregation in the concentration range 1-7 mm Such a combination of favorable biophysical and biological properties suggests that SCIs could provide a global therapeutic platform without a cold chain.

Identification of functional tetramolecular RNA G-quadruplexes derived from transfer RNAs

RNA G-quadruplex (RG4) structures are involved in multiple biological processes. Recent genome-wide analyses of human mRNA transcriptome identified thousands of putative intramolecular RG4s that readily assemble in vitro but shown to be unfolded in vivo. Previously, we have shown that mature cytoplasmic tRNAs are cleaved during stress response to produce tRNA fragments that function to repress translation in vivo. Here we report that these bioactive tRNA fragments assemble into intermolecular RG4s. We provide evidence for the formation of uniquely stable tetramolecular RG4 structures consisting of five tetrad layers formed by 5ʹ-terminal oligoguanine motifs of an individual tRNA fragment. RG4 is required for functions of tRNA fragments in the regulation of mRNA translation, a critical component of cellular stress response. RG4 disruption abrogates tRNA fragments ability to trigger the formation of Stress Granules in vivo. Collectively, our data rationalize the existence of naturally occurring RG4-assembling tRNA fragments and emphasize their regulatory roles.

RNA G-quadruplexes (RG4) occur in vivo and have regulatory roles in mRNA metabolism. Here the authors show that the guanine residue stretches at the 5’ end of tRNA-derived stress-induced RNAs mediate the formation of tetramolecular RG4 structures, which play a role in the post-transcriptional regulation of gene expression.

Extending Halogen-based Medicinal Chemistry to Proteins: IODO-INSULIN AS A CASE STUDY

Insulin, a protein critical for metabolic homeostasis, provides a classical model for protein design with application to human health. Recent efforts to improve its pharmaceutical formulation demonstrated that iodination of a conserved tyrosine (TyrB26) enhances key properties of a rapid-acting clinical analog. Moreover, the broad utility of halogens in medicinal chemistry has motivated the use of hybrid quantum- and molecular-mechanical methods to study proteins. Here, we (i) undertook quantitative atomistic simulations of 3-[iodo-TyrB26]insulin to predict its structural features, and (ii) tested these predictions by X-ray crystallography. Using an electrostatic model of the modified aromatic ring based on quantum chemistry, the calculations suggested that the analog, as a dimer and hexamer, exhibits subtle differences in aromatic-aromatic interactions at the dimer interface. Aromatic rings (TyrB16, PheB24, PheB25, 3-I-TyrB26, and their symmetry-related mates) at this interface adjust to enable packing of the hydrophobic iodine atoms within the core of each monomer. Strikingly, these features were observed in the crystal structure of a 3-[iodo-TyrB26]insulin analog (determined as an R6 zinc hexamer). Given that residues B24-B30 detach from the core on receptor binding, the environment of 3-I-TyrB26 in a receptor complex must differ from that in the free hormone. Based on the recent structure of a "micro-receptor" complex, we predict that 3-I-TyrB26 engages the receptor via directional halogen bonding and halogen-directed hydrogen bonding as follows: favorable electrostatic interactions exploiting, respectively, the halogen's electron-deficient σ-hole and electronegative equatorial band. Inspired by quantum chemistry and molecular dynamics, such "halogen engineering" promises to extend principles of medicinal chemistry to proteins.

Tetracycline determines the conformation of its aptamer at physiological magnesium concentrations

Synthetic riboswitches are versatile tools for the study and manipulation of biological systems. Yet, the underlying mechanisms governing its structural properties and regulation under physiological conditions are poorly studied. We performed spectroscopic and calorimetric experiments to explore the folding kinetics and thermodynamics of the tetracycline-binding aptamer, which can be employed as synthetic riboswitch, in the range of physiological magnesium concentrations. The dissociation constant of the ligand-aptamer complex was found to strongly depend on the magnesium concentration. At physiological magnesium concentrations, tetracycline induces a significant conformational shift from a compact, but heterogeneous intermediate state toward the completely formed set of tertiary interactions defining the regulation-competent structure. Thus, the switching functionality of the tetracycline-binding aptamer appears to include both a conformational rearrangement toward the regulation-competent structure and its thermodynamic stabilization.

Contribution of TyrB26 to the Function and Stability of Insulin: STRUCTURE-ACTIVITY RELATIONSHIPS AT A CONSERVED HORMONE-RECEPTOR INTERFACE

Crystallographic studies of insulin bound to receptor domains have defined the primary hormone-receptor interface. We investigated the role of Tyr(B26), a conserved aromatic residue at this interface. To probe the evolutionary basis for such conservation, we constructed 18 variants at B26. Surprisingly, non-aromatic polar or charged side chains (such as Glu, Ser, or ornithine (Orn)) conferred high activity, whereas the weakest-binding analogs contained Val, Ile, and Leu substitutions. Modeling of variant complexes suggested that the B26 side chains pack within a shallow depression at the solvent-exposed periphery of the interface. This interface would disfavor large aliphatic side chains. The analogs with highest activity exhibited reduced thermodynamic stability and heightened susceptibility to fibrillation. Perturbed self-assembly was also demonstrated in studies of the charged variants (Orn and Glu); indeed, the Glu(B26) analog exhibited aberrant aggregation in either the presence or absence of zinc ions. Thus, although Tyr(B26) is part of insulin's receptor-binding surface, our results suggest that its conservation has been enjoined by the aromatic ring's contributions to native stability and self-assembly. We envisage that such classical structural relationships reflect the implicit threat of toxic misfolding (rather than hormonal function at the receptor level) as a general evolutionary determinant of extant protein sequences.

Interaction of the tetracyclines with double-stranded RNAs of random base sequence: new perspectives on the target and mechanism of action

The 16S rRNA binding mechanism proposed for the antibacterial action of the tetracyclines does not explain their mechanism of action against non-bacterial pathogens. In addition, several contradictory base pairs have been proposed as their binding sites on the 16S rRNA. This study investigated the binding of minocycline and doxycycline to short double-stranded RNAs (dsRNAs) of random base sequences. These tetracyclines caused a dose-dependent decrease in the fluorescence intensities of 6-carboxyfluorescein (FAM)-labelled dsRNA and ethidium bromide (EtBr)-stained dsRNA, indicating that both drugs bind to dsRNA of random base sequence in a manner that is competitive with the binding of EtBr and other nucleic acid ligands often used as stains. This effect was observable in the presence of Mg(2+). The binding of the tetracyclines to dsRNA changed features of the fluorescence emission spectra of the drugs and the CD spectra of the RNA, and inhibited RNase III cleavage of the dsRNA. These results indicate that the double-stranded structures of RNAs may have a more important role in their interaction with the tetracyclines than the specific base pairs, which had hitherto been the subject of much investigation. Given the diverse functions of cellular RNAs, the binding of the tetracyclines to their double-stranded helixes may alter the normal processing and functioning of the various biological processes they regulate. This could help to explain the wide range of action of the tetracyclines against various pathogens and disease conditions.

Full characterization and cytotoxic activity of new silver(i) and copper(i) helicates with quaterpyridine

DNA binding and cells' cycle disorders caused by four new dinuclear Ag(i) and Cu(i) double helicates with quaterpyridine ligands.

Using sm-FRET and denaturants to reveal folding landscapes

RNA folding studies aim to clarify the relationship among sequence, tridimensional structure, and biological function. In the last decade, the application of single-molecule fluorescence resonance energy transfer (sm-FRET) techniques to investigate RNA structure and folding has revealed the details of conformational changes and timescale of the process leading to the formation of biologically active RNA structures with subnanometer resolution on millisecond timescales. In this review, we initially summarize the first wave of single-molecule FRET-based RNA techniques that focused on analyzing the influence of mono- and divalent metal ions on RNA function, and how these studies have provided very valuable information about folding pathways and the presence of intermediate and low-populated states. Next, we describe a second generation of single-molecule techniques that combine sm-FRET with the use of chemical denaturants as an emerging powerful approach to reveal information about the dynamics and energetics of RNA folding that remains hidden using conventional sm-FRET approaches. The main advantages of using the competing interplay between folding agents such as metal ions and denaturants to observe and manipulate the dynamics of RNA folding and RNA-ligand interactions is discussed in the context of the adenine riboswitch aptamer.

Distinctive structural motifs of RNA G-quadruplexes composed of AGG, CGG and UGG trinucleotide repeats

Trinucleotide repeats are microsatellite sequences that are polymorphic in length. Their expansion in specific genes underlies a number of neurodegenerative disorders. Using ultraviolet-visible, circular dichroism, nuclear magnetic resonance (NMR) spectroscopies and electrospray ionization mass spectrometry, the structural preferences of RNA molecules composed of two and four repeats of AGG, CGG and UGG in the presence of K⁺, Na⁺ and NH₄⁺ were analysed. (AGG)₂A, (AGG)₄A, p(UGG)₂U and p(UGG)₄U strongly prefer folding into G-quadruplexes, whereas CGG-containing sequences can adopt different types of structure depending on the cation and on the number of repeats. In particular, the two-repeat CGG sequence folds into a G-quadruplex in potassium buffer. We also found that each G-quadruplex fold is different: A:(G:G:G:G)A hexads were found for (AGG)₂A, whereas mixed G:C:G:C tetrads and U-tetrads were observed in the NMR spectra of G(CGG)₂C and p(UGG)₂U, respectively. Finally, our NMR study highlights the influence of the strand sequence on the structure formed, and the influence of the intracellular environment on the folding. Importantly, we highlight that although potassium ions are prevalent in cells, the structures observed in the HeLa cell extract are not always the same as those prevailing in biophysical studies in the presence of K⁺ ions.

Biophysical optimization of a therapeutic protein by nonstandard mutagenesis: studies of an iodo-insulin derivative

Insulin provides a model for the therapeutic application of protein engineering. A paradigm in molecular pharmacology was defined by design of rapid-acting insulin analogs for the prandial control of glycemia. Such analogs, a cornerstone of current diabetes regimens, exhibit accelerated subcutaneous absorption due to more rapid disassembly of oligomeric species relative to wild-type insulin. This strategy is limited by a molecular trade-off between accelerated disassembly and enhanced susceptibility to degradation. Here, we demonstrate that this trade-off may be circumvented by nonstandard mutagenesis. Our studies employed Lys(B28), Pro(B29)-insulin ("lispro") as a model prandial analog that is less thermodynamically stable and more susceptible to fibrillation than is wild-type insulin. We have discovered that substitution of an invariant tyrosine adjoining the engineered sites in lispro (Tyr(B26)) by 3-iodo-Tyr (i) augments its thermodynamic stability (ΔΔGu 0.5 ± 0.2 kcal/mol), (ii) delays onset of fibrillation (lag time on gentle agitation at 37 °C was prolonged by 4-fold), (iii) enhances affinity for the insulin receptor (1.5 ± 0.1-fold), and (iv) preserves biological activity in a rat model of diabetes mellitus. (1)H NMR studies suggest that the bulky iodo-substituent packs within a nonpolar interchain crevice. Remarkably, the 3-iodo-Tyr(B26) modification stabilizes an oligomeric form of insulin pertinent to pharmaceutical formulation (the R6 zinc hexamer) but preserves rapid disassembly of the oligomeric form pertinent to subcutaneous absorption (T6 hexamer). By exploiting this allosteric switch, 3-iodo-Tyr(B26)-lispro thus illustrates how a nonstandard amino acid substitution can mitigate the unfavorable biophysical properties of an engineered protein while retaining its advantages.

Comparative study of two box H/ACA ribonucleoprotein pseudouridine-synthases: relation between conformational dynamics of the guide RNA, enzyme assembly and activity

Multiple RNA-guided pseudouridine synthases, H/ACA ribonucleoprotein particles (RNPs) which contain a guide RNA and four proteins, catalyze site-specific post-transcriptional isomerization of uridines into pseudouridines in substrate RNAs. In archaeal particles, the guide small RNA (sRNA) is anchored by the pseudouridine synthase aCBF5 and the ribosomal protein L7Ae. Protein aNOP10 interacts with both aCBF5 and L7Ae. The fourth protein, aGAR1, interacts with aCBF5 and enhances catalytic efficiency. Here, we compared the features of two H/ACA sRNAs, Pab21 and Pab91, from Pyrococcus abyssi. We found that aCBF5 binds much more weakly to Pab91 than to Pab21. Surprisingly, the Pab91 sRNP exhibits a higher catalytic efficiency than the Pab21 sRNP. We thus investigated the molecular basis of the differential efficiencies observed for the assembly and catalytic activity of the two enzymes. For this, we compared profiles of the extent of lead-induced cleavages in these sRNAs during a stepwise reconstitution of the sRNPs, and analyzed the impact of the absence of the aNOP10–L7Ae interaction. Such probing experiments indicated that the sRNAs undergo a series of conformational changes upon RNP assembly. These changes were also evaluated directly by circular dichroism (CD) spectroscopy, a tool highly adapted to analyzing RNA conformational dynamics. In addition, our results reveal that the conformation of helix P1 formed at the base of the H/ACA sRNAs is optimized in Pab21 for efficient aCBF5 binding and RNP assembly. Moreover, P1 swapping improved the assembly of the Pab91 sRNP. Nonetheless, efficient aCBF5 binding probably also relies on the pseudouridylation pocket which is not optimized for high activity in the case of Pab21.

Intermolecular domain docking in the hairpin ribozyme: metal dependence, binding kinetics and catalysis

The hairpin ribozyme is a prototype small, self-cleaving RNA motif. It exists naturally as a four-way RNA junction containing two internal loops on adjoining arms. These two loops interact in a cation-driven docking step prior to chemical catalysis to form a tightly integrated structure, with dramatic changes occurring in the conformation of each loop upon docking. We investigate the thermodynamics and kinetics of the docking process using constructs in which loop A and loop B reside on separate molecules. Using a novel CD difference assay to isolate the effects of metal ions linked to domain docking, we find the intermolecular docking process to be driven by sub-millimolar concentrations of the exchange-inert Co(NH₃)₆³⁺. RNA self-cleavage requires binding of lower-affinity ions with greater apparent cooperativity than the docking process itself, implying that, even in the absence of direct coordination to RNA, metal ions play a catalytic role in hairpin ribozyme function beyond simply driving loop-loop docking. Surface plasmon resonance assays reveal remarkably slow molecular association, given the relatively tight loop-loop interaction. This observation is consistent with a “double conformational capture” model in which only collisions between loop A and loop B molecules that are simultaneously in minor, docking-competent conformations are productive for binding.

Insulin fibrillation and protein design: topological resistance of single-chain analogs to thermal degradation with application to a pump reservoir

Insulin is susceptible to thermal fibrillation, a misfolding process that leads to nonnative cross-β assembly analogous to pathological amyloid deposition. Pharmaceutical formulations are ordinarily protected from such degradation by sequestration of the susceptible monomer within native protein assemblies. With respect to the safety and efficacy of insulin pumps, however, this strategy imposes an intrinsic trade-off between pharmacokinetic goals (rapid absorption and clearance) and the requisite physical properties of a formulation (prolonged shelf life and stability within the reservoir). Available rapid-acting formulations are suboptimal in both respects; susceptibility to fibrillation is exacerbated even as absorption is delayed relative to the ideal specifications of a closed-loop system. To circumvent this molecular trade-off, we exploited structural models of insulin fibrils and amyloidogenic intermediates to define an alternative protective mechanism. Single-chain insulin (SCI) analogs were shown to be refractory to thermal fibrillation with maintenance of biological activity for more than 3 months under conditions that promote the rapid fibrillation and inactivation of insulin. The essential idea exploits an intrinsic incompatibility between SCI topology and the geometry of cross-β assembly. A peptide tether was thus interposed between the A- and B-chains whose length was (a) sufficiently long to provide the "play" needed for induced fit of the hormone on receptor binding and yet (b) sufficiently short to impose a topological barrier to fibrillation. Our findings suggest that ultrastable monomeric SCI analogs may be formulated without protective self-assembly and so permit simultaneous optimization of pharmacokinetics and reservoir life.

Conformational dynamics of insulin

We have exploited a prandial insulin analog to elucidate the underlying structure and dynamics of insulin as a monomer in solution. A model was provided by insulin lispro (the active component of Humalog(®); Eli Lilly and Co.). Whereas NMR-based modeling recapitulated structural relationships of insulin crystals (T-state protomers), dynamic anomalies were revealed by amide-proton exchange kinetics in D(2)O. Surprisingly, the majority of hydrogen bonds observed in crystal structures are only transiently maintained in solution, including key T-state-specific inter-chain contacts. Long-lived hydrogen bonds (as defined by global exchange kinetics) exist only at a subset of four α-helical sites (two per chain) flanking an internal disulfide bridge (cystine A20-B19); these sites map within the proposed folding nucleus of proinsulin. The anomalous flexibility of insulin otherwise spans its active surface and may facilitate receptor binding. Because conformational fluctuations promote the degradation of pharmaceutical formulations, we envisage that "dynamic re-engineering" of insulin may enable design of ultra-stable formulations for humanitarian use in the developing world.

Deciphering the hidden informational content of protein sequences: foldability of proinsulin hinges on a flexible arm that is dispensable in the mature hormone

Protein sequences encode both structure and foldability. Whereas the interrelationship of sequence and structure has been extensively investigated, the origins of folding efficiency are enigmatic. We demonstrate that the folding of proinsulin requires a flexible N-terminal hydrophobic residue that is dispensable for the structure, activity, and stability of the mature hormone. This residue (Phe(B1) in placental mammals) is variably positioned within crystal structures and exhibits (1)H NMR motional narrowing in solution. Despite such flexibility, its deletion impaired insulin chain combination and led in cell culture to formation of non-native disulfide isomers with impaired secretion of the variant proinsulin. Cellular folding and secretion were maintained by hydrophobic substitutions at B1 but markedly perturbed by polar or charged side chains. We propose that, during folding, a hydrophobic side chain at B1 anchors transient long-range interactions by a flexible N-terminal arm (residues B1-B8) to mediate kinetic or thermodynamic partitioning among disulfide intermediates. Evidence for the overall contribution of the arm to folding was obtained by alanine scanning mutagenesis. Together, our findings demonstrate that efficient folding of proinsulin requires N-terminal sequences that are dispensable in the native state. Such arm-dependent folding can be abrogated by mutations associated with β-cell dysfunction and neonatal diabetes mellitus.

A synonymous single nucleotide polymorphism in DeltaF508 CFTR alters the secondary structure of the mRNA and the expression of the mutant protein

Recent advances in our understanding of translational dynamics indicate that codon usage and mRNA secondary structure influence translation and protein folding. The most frequent cause of cystic fibrosis (CF) is the deletion of three nucleotides (CTT) from the cystic fibrosis transmembrane conductance regulator (CFTR) gene that includes the last cytosine (C) of isoleucine 507 (Ile507ATC) and the two thymidines (T) of phenylalanine 508 (Phe508TTT) codons. The consequences of the deletion are the loss of phenylalanine at the 508 position of the CFTR protein (DeltaF508), a synonymous codon change for isoleucine 507 (Ile507ATT), and protein misfolding. Here we demonstrate that the DeltaF508 mutation alters the secondary structure of the CFTR mRNA. Molecular modeling predicts and RNase assays support the presence of two enlarged single stranded loops in the DeltaF508 CFTR mRNA in the vicinity of the mutation. The consequence of DeltaF508 CFTR mRNA "misfolding" is decreased translational rate. A synonymous single nucleotide variant of the DeltaF508 CFTR (Ile507ATC), that could exist naturally if Phe-508 was encoded by TTC, has wild type-like mRNA structure, and enhanced expression levels when compared with native DeltaF508 CFTR. Because CFTR folding is predominantly cotranslational, changes in translational dynamics may promote DeltaF508 CFTR misfolding. Therefore, we propose that mRNA "misfolding" contributes to DeltaF508 CFTR protein misfolding and consequently to the severity of the human DeltaF508 phenotype. Our studies suggest that in addition to modifier genes, SNPs may also contribute to the differences observed in the symptoms of various DeltaF508 homozygous CF patients.

An Achilles' heel in an amyloidogenic protein and its repair: insulin fibrillation and therapeutic design

Insulin fibrillation provides a model for a broad class of amyloidogenic diseases. Conformational distortion of the native monomer leads to aggregation-coupled misfolding. Whereas beta-cells are protected from proteotoxicity by hexamer assembly, fibrillation limits the storage and use of insulin at elevated temperatures. Here, we have investigated conformational distortions of an engineered insulin monomer in relation to the structure of an insulin fibril. Anomalous (13)C NMR chemical shifts and rapid (15)N-detected (1)H-(2)H amide-proton exchange were observed in one of the three classical alpha-helices (residues A1-A8) of the hormone, suggesting a conformational equilibrium between locally folded and unfolded A-chain segments. Whereas hexamer assembly resolves these anomalies in accordance with its protective role, solid-state (13)C NMR studies suggest that the A-chain segment participates in a fibril-specific beta-sheet. Accordingly, we investigated whether helicogenic substitutions in the A1-A8 segment might delay fibrillation. Simultaneous substitution of three beta-branched residues (Ile(A2) --> Leu, Val(A3) --> Leu, and Thr(A8) --> His) yielded an analog with reduced thermodynamic stability but marked resistance to fibrillation. Whereas amide-proton exchange in the A1-A8 segment remained rapid, (13)Calpha chemical shifts exhibited a more helical pattern. This analog is essentially without activity, however, as Ile(A2) and Val(A3) define conserved receptor contacts. To obtain active analogs, substitutions were restricted to A8. These analogs exhibit high receptor-binding affinity; representative potency in a rodent model of diabetes mellitus was similar to wild-type insulin. Although (13)Calpha chemical shifts remain anomalous, significant protection from fibrillation is retained. Together, our studies define an "Achilles' heel" in a globular protein whose repair may enhance the stability of pharmaceutical formulations and broaden their therapeutic deployment in the developing world.

Dissecting the influence of Mg2+ on 3D architecture and ligand-binding of the guanine-sensing riboswitch aptamer domain

Long-range tertiary interactions determine the three-dimensional structure of a number of metabolite-binding riboswitch RNA elements and were found to be important for their regulatory function. For the guanine-sensing riboswitch of the Bacillus subtilis xpt-pbuX operon, our previous NMR-spectroscopic studies indicated pre-formation of long-range tertiary contacts in the ligand-free state of its aptamer domain. Loss of the structural pre-organization in a mutant of this RNA (G37A/C61U) resulted in the requirement of Mg²⁺ for ligand binding. Here, we investigate structural and stability aspects of the wild-type aptamer domain (Gsw) and the G37A/C61U-mutant (Gsw^loop) of the guanine-sensing riboswitch and their Mg²⁺-induced folding characteristics to dissect the role of long-range tertiary interactions, the link between pre-formation of structural elements and ligand-binding properties and the functional stability. Destabilization of the long-range interactions as a result of the introduced mutations for Gsw^loop or the increase in temperature for both Gsw and Gsw^loop involves pronounced alterations of the conformational ensemble characteristics of the ligand-free state of the riboswitch. The increased flexibility of the conformational ensemble can, however, be compensated by Mg²⁺. We propose that reduction of conformational dynamics in remote regions of the riboswitch aptamer domain is the minimal pre-requisite to pre-organize the core region for specific ligand binding.

Extended structures in RNA folding intermediates are due to nonnative interactions rather than electrostatic repulsion

RNA folding occurs via a series of transitions between metastable intermediate states for Mg(2+) concentrations below those needed to fold the native structure. In general, these folding intermediates are considerably less compact than their respective native states. Our previous work demonstrates that the major equilibrium intermediate of the 154-residue specificity domain (S-domain) of the Bacillus subtilis RNase P RNA is more extended than its native structure. We now investigate two models with falsifiable predictions regarding the origins of the extended intermediate structures in the S-domains of the B. subtilis and the Escherichia coli RNase P RNA that belong to different classes of P RNA and have distinct native structures. The first model explores the contribution of electrostatic repulsion, while the second model probes specific interactions in the core of the folding intermediate. Using small-angle X-ray scattering and Langevin dynamics simulations, we show that electrostatics plays only a minor role, whereas specific interactions largely account for the extended nature of the intermediate. Structural contacts in the core, including a nonnative base pair, help to stabilize the intermediate conformation. We conclude that RNA folding intermediates adopt extended conformations due to short-range, nonnative interactions rather than generic electrostatic repulsion of helical domains. These principles apply to other ribozymes and riboswitches that undergo functionally relevant conformational changes.

Contribution of residue B5 to the folding and function of insulin and IGF-I: constraints and fine-tuning in the evolution of a protein family

Proinsulin exhibits a single structure, whereas insulin-like growth factors refold as two disulfide isomers in equilibrium. Native insulin-related growth factor (IGF)-I has canonical cystines (A6-A11, A7-B7, and A20-B19) maintained by IGF-binding proteins; IGF-swap has alternative pairing (A7-A11, A6-B7, and A20-B19) and impaired activity. Studies of mini-domain models suggest that residue B5 (His in insulin and Thr in IGFs) governs the ambiguity or uniqueness of disulfide pairing. Residue B5, a site of mutation in proinsulin causing neonatal diabetes, is thus of broad biophysical interest. Here, we characterize reciprocal B5 substitutions in the two proteins. In insulin, His(B5) --> Thr markedly destabilizes the hormone (DeltaDeltaG(u) 2.0 +/- 0.2 kcal/mol), impairs chain combination, and blocks cellular secretion of proinsulin. The reciprocal IGF-I substitution Thr(B5) --> His (residue 4) specifies a unique structure with native (1)H NMR signature. Chemical shifts and nuclear Overhauser effects are similar to those of native IGF-I. Whereas wild-type IGF-I undergoes thiol-catalyzed disulfide exchange to yield IGF-swap, His(B5)-IGF-I retains canonical pairing. Chemical denaturation studies indicate that His(B5) does not significantly enhance thermodynamic stability (DeltaDeltaG(u) 0.2 +/- 0.2 kcal/mol), implying that the substitution favors canonical pairing by destabilizing competing folds. Whereas the activity of Thr(B5)-insulin is decreased 5-fold, His(B5)-IGF-I exhibits 2-fold increased affinity for the IGF receptor and augmented post-receptor signaling. We propose that conservation of Thr(B5) in IGF-I, rescued from structural ambiguity by IGF-binding proteins, reflects fine-tuning of signal transduction. In contrast, the conservation of His(B5) in insulin highlights its critical role in insulin biosynthesis.

Crystal structure of a "nonfoldable" insulin: impaired folding efficiency despite native activity

Protein evolution is constrained by folding efficiency ("foldability") and the implicit threat of toxic misfolding. A model is provided by proinsulin, whose misfolding is associated with beta-cell dysfunction and diabetes mellitus. An insulin analogue containing a subtle core substitution (Leu(A16) --> Val) is biologically active, and its crystal structure recapitulates that of the wild-type protein. As a seeming paradox, however, Val(A16) blocks both insulin chain combination and the in vitro refolding of proinsulin. Disulfide pairing in mammalian cell culture is likewise inefficient, leading to misfolding, endoplasmic reticular stress, and proteosome-mediated degradation. Val(A16) destabilizes the native state and so presumably perturbs a partial fold that directs initial disulfide pairing. Substitutions elsewhere in the core similarly destabilize the native state but, unlike Val(A16), preserve folding efficiency. We propose that Leu(A16) stabilizes nonlocal interactions between nascent alpha-helices in the A- and B-domains to facilitate initial pairing of Cys(A20) and Cys(B19), thus surmounting their wide separation in sequence. Although Val(A16) is likely to destabilize this proto-core, its structural effects are mitigated once folding is achieved. Classical studies of insulin chain combination in vitro have illuminated the impact of off-pathway reactions on the efficiency of native disulfide pairing. The capability of a polypeptide sequence to fold within the endoplasmic reticulum may likewise be influenced by kinetic or thermodynamic partitioning among on- and off-pathway disulfide intermediates. The properties of [Val(A16)]insulin and [Val(A16)]proinsulin demonstrate that essential contributions of conserved residues to folding may be inapparent once the native state is achieved.

Interplay of structure, hydration and thermal stability in formacetal modified oligonucleotides: RNA may tolerate nonionic modifications better than DNA

DNA and RNA oligonucleotides having formacetal internucleoside linkages between uridine and adenosine nucleosides have been prepared and studied using UV thermal melting, osmotic stress, and X-ray crystallography. Formacetal modifications have remarkably different effects on double helical RNA and DNA-the formacetal stabilizes the RNA helix by +0.7 degrees C but destabilizes the DNA helix by -1.6 degrees C per modification. The apparently hydrophobic formacetal has little effect on hydration of RNA but decreases the hydration of DNA, which suggests that at least part of the difference in thermal stability may be related to differences in hydration. A crystal structure of modified DNA shows that two isolated formacetal linkages fit almost perfectly in an A-type helix (decamer). Taken together, the data suggest that RNA may tolerate nonionic backbone modifications better than DNA. Overall, formacetal appears to be an excellent mimic of phosphate linkage in RNA and an interesting modification for potential applications in fundamental studies and RNA-based gene control strategies, such as RNA interference.

Design of an insulin analog with enhanced receptor binding selectivity: rationale, structure, and therapeutic implications

Insulin binds with high affinity to the insulin receptor (IR) and with low affinity to the type 1 insulin-like growth factor (IGF) receptor (IGFR). Such cross-binding, which reflects homologies within the insulin-IGF signaling system, is of clinical interest in relation to the association between hyperinsulinemia and colorectal cancer. Here, we employ nonstandard mutagenesis to design an insulin analog with enhanced affinity for the IR but reduced affinity for the IGFR. Unnatural amino acids were introduced by chemical synthesis at the N- and C-capping positions of a recognition alpha-helix (residues A1 and A8). These sites adjoin the hormone-receptor interface as indicated by photocross-linking studies. Specificity is enhanced more than 3-fold on the following: (i) substitution of Gly(A1) by D-Ala or D-Leu, and (ii) substitution of Thr(A8) by diaminobutyric acid (Dab). The crystal structure of [D-Ala(A1),Dab(A8)]insulin, as determined within a T(6) zinc hexamer to a resolution of 1.35 A, is essentially identical to that of human insulin. The nonstandard side chains project into solvent at the edge of a conserved receptor-binding surface shared by insulin and IGF-I. Our results demonstrate that modifications at this edge discriminate between IR and IGFR. Because hyperinsulinemia is typically characterized by a 3-fold increase in integrated postprandial insulin concentrations, we envisage that such insulin analogs may facilitate studies of the initiation and progression of cancer in animal models. Future development of clinical analogs lacking significant IGFR cross-binding may enhance the safety of insulin replacement therapy in patients with type 2 diabetes mellitus at increased risk of colorectal cancer.

Enhancing the activity of a protein by stereospecific unfolding: conformational life cycle of insulin and its evolutionary origins

A central tenet of molecular biology holds that the function of a protein is mediated by its structure. An inactive ground-state conformation may nonetheless be enjoined by the interplay of competing biological constraints. A model is provided by insulin, well characterized at atomic resolution by x-ray crystallography. Here, we demonstrate that the activity of the hormone is enhanced by stereospecific unfolding of a conserved structural element. A bifunctional beta-strand mediates both self-assembly (within beta-cell storage vesicles) and receptor binding (in the bloodstream). This strand is anchored by an invariant side chain (Phe(B24)); its substitution by Ala leads to an unstable but native-like analog of low activity. Substitution by d-Ala is equally destabilizing, and yet the protein diastereomer exhibits enhanced activity with segmental unfolding of the beta-strand. Corresponding photoactivable derivatives (containing l- or d-para-azido-Phe) cross-link to the insulin receptor with higher d-specific efficiency. Aberrant exposure of hydrophobic surfaces in the analogs is associated with accelerated fibrillation, a form of aggregation-coupled misfolding associated with cellular toxicity. Conservation of Phe(B24), enforced by its dual role in native self-assembly and induced fit, thus highlights the implicit role of misfolding as an evolutionary constraint. Whereas classical crystal structures of insulin depict its storage form, signaling requires engagement of a detachable arm at an extended receptor interface. Because this active conformation resembles an amyloidogenic intermediate, we envisage that induced fit and self-assembly represent complementary molecular adaptations to potential proteotoxicity. The cryptic threat of misfolding poses a universal constraint in the evolution of polypeptide sequences.

Analysis of natural variants of the hepatitis C virus internal ribosome entry site reveals that primary sequence plays a key role in cap-independent translation

The HCV internal ribosome entry site (IRES) spans a region of ∼340 nt that encompasses most of the 5′ untranslated region (5′UTR) of the viral mRNA and the first 24–40 nt of the core-coding region. To investigate the implication of altering the primary sequence of the 5′UTR on IRES activity, naturally occurring variants of the 5′UTR were isolated from clinical samples and analyzed. The impact of the identified mutations on translation was evaluated in the context of RLuc/FLuc bicistronic RNAs. Results show that depending on their location within the RNA structure, these naturally occurring mutations cause a range of effects on IRES activity. However, mutations within subdomain IIId hinder HCV IRES-mediated translation. In an attempt to explain these data, the dynamic behavior of the subdomain IIId was analyzed by means of molecular dynamics (MD) simulations. Despite the loss of function, MD simulations predicted that mutant G266A/G268U possesses a structure similar to the wt-RNA. This prediction was validated by analyzing the secondary structure of the isolated IIId RNAs by circular dichroism spectroscopy in the presence or absence of Mg²⁺ ions. These data strongly suggest that the primary sequence of subdomain IIId plays a key role in HCV IRES-mediated translation.

The structure of a mutant insulin uncouples receptor binding from protein allostery. An electrostatic block to the TR transition

The zinc insulin hexamer undergoes allosteric reorganization among three conformational states, designated T(6), T(3)R(3)(f), and R(6). Although the free monomer in solution (the active species) resembles the classical T-state, an R-like conformational change is proposed to occur upon receptor binding. Here, we distinguish between the conformational requirements of receptor binding and the crystallographic TR transition by design of an active variant refractory to such reorganization. Our strategy exploits the contrasting environments of His(B5) in wild-type structures: on the T(6) surface but within an intersubunit crevice in R-containing hexamers. The TR transition is associated with a marked reduction in His(B5) pK(a), in turn predicting that a positive charge at this site would destabilize the R-specific crevice. Remarkably, substitution of His(B5) (conserved among eutherian mammals) by Arg (occasionally observed among other vertebrates) blocks the TR transition, as probed in solution by optical spectroscopy. Similarly, crystallization of Arg(B5)-insulin in the presence of phenol (ordinarily a potent inducer of the TR transition) yields T(6) hexamers rather than R(6) as obtained in control studies of wild-type insulin. The variant structure, determined at a resolution of 1.3A, closely resembles the wild-type T(6) hexamer. Whereas Arg(B5) is exposed on the protein surface, its side chain participates in a solvent-stabilized network of contacts similar to those involving His(B5) in wild-type T-states. The substantial receptor-binding activity of Arg(B5)-insulin (40% relative to wild type) demonstrates that the function of an insulin monomer can be uncoupled from its allosteric reorganization within zinc-stabilized hexamers.

Characterizing riboswitch function: identification of Mg2+ binding site in T box antiterminator RNA

T box bacterial genes utilize a riboswitch mechanism to regulate gene expression at the transcriptional level. Complementary base pairing of the 5'-untranslated mRNA with uncharged cognate tRNA stabilizes formation of an antiterminator element and permits complete transcription. In the absence of tRNA, a mutually exclusive RNA terminator element forms and results in transcription termination. This regulatory mechanism requires divalent metal ions at the antitermination event. The structural effects of Mg(2+) binding to antiterminator model RNA were investigated to ascertain if this requirement is due to the presence of a specific metal ion binding site in the antiterminator. Spectroscopic analysis identified the presence of a hydrated, diffuse Mg(2+) binding site. The results indicate that the mechanistic requirement for divalent metal ions is not due to Mg(2+)-induced pre-formation of a functional antiterminator receptor; rather, Mg(2+) binds in a helical region of high phylogenetic sequence conservation adjacent to the tRNA binding site.

Design of an active ultrastable single-chain insulin analog: synthesis, structure, and therapeutic implications

Single-chain insulin (SCI) analogs provide insight into the inter-relation of hormone structure, function, and dynamics. Although compatible with wild-type structure, short connecting segments (<3 residues) prevent induced fit upon receptor binding and so are essentially without biological activity. Substantial but incomplete activity can be regained with increasing linker length. Here, we describe the design, structure, and function of a single-chain insulin analog (SCI-57) containing a 6-residue linker (GGGPRR). Native receptor-binding affinity (130 +/- 8% relative to the wild type) is achieved as hindrance by the linker is offset by favorable substitutions in the insulin moiety. The thermodynamic stability of SCI-57 is markedly increased (DeltaDeltaG(u) = 0.7 +/- 0.1 kcal/mol relative to the corresponding two-chain analog and 1.9 +/- 0.1 kcal/mol relative to wild-type insulin). Analysis of inter-residue nuclear Overhauser effects demonstrates that a native-like fold is maintained in solution. Surprisingly, the glycine-rich connecting segment folds against the insulin moiety: its central Pro contacts Val(A3) at the edge of the hydrophobic core, whereas the final Arg extends the A1-A8 alpha-helix. Comparison between SCI-57 and its parent two-chain analog reveals striking enhancement of multiple native-like nuclear Overhauser effects within the tethered protein. These contacts are consistent with wild-type crystal structures but are ordinarily attenuated in NMR spectra of two-chain analogs, presumably due to conformational fluctuations. Linker-specific damping of fluctuations provides evidence for the intrinsic flexibility of an insulin monomer. In addition to their biophysical interest, ultrastable SCIs may enhance the safety and efficacy of insulin replacement therapy in the developing world.

A large collapsed-state RNA can exhibit simple exponential single-molecule dynamics

The process of large RNA folding is believed to proceed from many collapsed structures to a unique functional structure requiring precise organization of nucleotides. The diversity of possible structures and stabilities of large RNAs could result in non-exponential folding kinetics (e.g. stretched exponential) under conditions where the molecules have not achieved their native state. We describe a single-molecule fluorescence resonance energy transfer (FRET) study of the collapsed-state region of the free energy landscape of the catalytic domain of RNase P RNA from Bacillus stearothermophilus (C(thermo)). Ensemble measurements have shown that this 260 residue RNA folds cooperatively to its native state at >or=1 mM Mg(2+), but little is known about the conformational dynamics at lower ionic strength. Our measurements of equilibrium conformational fluctuations reveal simple exponential kinetics that reflect a small number of discrete states instead of the expected inhomogeneous dynamics. The distribution of discrete dwell times, collected from an "ensemble" of 300 single molecules at each of a series of Mg(2+) concentrations, fit well to a double exponential, which indicates that the RNA conformational changes can be described as a four-state system. This finding is somewhat unexpected under [Mg(2+)] conditions in which this RNA does not achieve its native state. Observation of discrete well-defined conformations in this large RNA that are stable on the seconds timescale at low [Mg(2+)] (<0.1 mM) suggests that even at low ionic strength, with a tremendous number of possible (weak) interactions, a few critical interactions may produce deep energy wells that allow for rapid averaging of motions within each well, and yield kinetics that are relatively simple.

The A-chain of Insulin Contacts the Insert Domain of the Insulin Receptor

The contribution of the insulin A-chain to receptor binding is investigated by photo-cross-linking and nonstandard mutagenesis. Studies focus on the role of Val(A3), which projects within a crevice between the A- and B-chains. Engineered receptor alpha-subunits containing specific protease sites ("midi-receptors") are employed to map the site of photo-cross-linking by an analog containing a photoactivable A3 side chain (para-azido-Phe (Pap)). The probe cross-links to a C-terminal peptide (residues 703-719 of the receptor A isoform, KTFEDYLHNVVFVPRPS) containing side chains critical for hormone binding (underlined); the corresponding segment of the holoreceptor was shown previously to cross-link to a Pap(B25)-insulin analog. Because Pap is larger than Val and so may protrude beyond the A3-associated crevice, we investigated analogs containing A3 substitutions comparable in size to Val as follows: Thr, allo-Thr, and alpha-aminobutyric acid (Aba). Substitutions were introduced within an engineered monomer. Whereas previous studies of smaller substitutions (Gly(A3) and Ser(A3)) encountered nonlocal conformational perturbations, NMR structures of the present analogs are similar to wild-type insulin; the variant side chains are accommodated within a native-like crevice with minimal distortion. Receptor binding activities of Aba(A3) and allo-Thr(A3) analogs are reduced at least 10-fold; the activity of Thr(A3)-DKP-insulin is reduced 5-fold. The hormone-receptor interface is presumably destabilized either by a packing defect (Aba(A3)) or by altered polarity (allo-Thr(A3) and Thr(A3)). Our results provide evidence that Val(A3), a site of mutation causing diabetes mellitus, contacts the insert domain-derived tail of the alpha-subunit in a hormone-receptor complex.

Folding of a universal ribozyme: the ribonuclease P RNA

Ribonuclease P is among the first ribozymes discovered, and is the only ubiquitously occurring ribozyme besides the ribosome. The bacterial RNase P RNA is catalytically active without its protein subunit and has been studied for over two decades as a model system for RNA catalysis, structure and folding. This review focuses on the thermodynamic, kinetic and structural frameworks derived from the folding studies of bacterial RNase P RNA.

Characterization of the mechanical unfolding of RNA pseudoknots

The pseudoknot is an important RNA structural element that provides an excellent model system for studying the contributions of tertiary interactions to RNA stability and to folding kinetics. RNA pseudoknots are also of interest because of their key role in the control of ribosomal frameshifting by viral RNAs. Their mechanical properties are directly relevant to their unfolding by ribosomes during translation. We have used optical tweezers to study the kinetics and thermodynamics of mechanical unfolding and refolding of single RNA molecules. Here we describe the unfolding of the frameshifting pseudoknot from infectious bronchitis virus (IBV), three constituent hairpins, and three mutants of the IBV pseudoknot. All four pseudoknots cause −1 programmed ribosomal frameshifting. We have measured the free energies and rates of mechanical unfolding and refolding of the four frameshifting pseudoknots. Our results show that the IBV pseudoknot requires a higher force than its corresponding hairpins to unfold. Furthermore, its rate of unfolding changes little with increasing force, in contrast with the rate of hairpin unfolding. The presence of Mg²⁺ significantly increases the kinetic barriers to unfolding the IBV pseudoknot, but has only a minor effect on the hairpin unfolding. The greater mechanical stability of pseudoknots compared to hairpins, and their kinetic insensitivity to force supports the hypothesis that −1 frameshifting depends on the difficulty of unfolding the mRNA.

Investigation ofde novo Totally Random Biosequences, Part III

Fold is essential to RNA properties, and, in particular, its thermodynamic stability can be used to monitor RNA-protein or RNA-ligand interactions, and to engineer RNA with novel or improved properties. While clearly valuable, experimental determination of RNA folding stability by traditional biophysical techniques requires substantial amounts of pure sample and rather expensive equipment. In this paper, we report a new, simple approach to the determination of RNA folding stability by coupling enzymatic digestion and temperature denaturation. The assay, named RNA folding stability Test (RNA Foster), is designed to probe the fraction of folded RNA (f(fold)) in an equilibrium mixture of folded and unfolded ones as a function of temperature. The simplicity of RNA Foster suggests that it can easily be scaled up for high-throughput studies of RNA folding stability both in basic and applied research.

A counterintuitive Mg2+-dependent and modification-assisted functional folding of mitochondrial tRNAs

Mitochondrial tRNAs (mtRNAs) often lack domains and posttranscriptional modifications that are found in cytoplasmic tRNAs. These structural and chemical elements normally stabilize the folding of cytoplasmic tRNAs into canonical structures that are competent for aminoacylation and translation. For example, the dihydrouridine (D) stem and loop domain is involved in the tertiary structure of cytoplasmic tRNAs through hydrogen bonds and a Mg2+ bridge to the ribothymidine (T) stem and loop domain. These interactions are often absent in mtRNA because the D-domain is truncated or missing. Using gel mobility shift analyses, UV, circular dichroism and NMR spectroscopies and aminoacylation assays, we have investigated the functional folding interactions of chemically synthesized and site-specifically modified mitochondrial and cytoplasmic tRNAs. We found that Mg2+ is critical for folding of the truncated D-domain of bovine mtRNAMet with the tRNA's T-domain. Contrary to the expectation that Mg2+ stabilizes RNA folding, the mtRNAMet D-domain structure was unfolded and relaxed, rather than stabilized in the presence of Mg2+. Because the D-domain is transcribed prior to the T-domain, we conclude that Mg2+ prevents misfolding of the 5'-half of bovine mtRNAMet facilitating its correct interaction with the T-domain. The interaction of the mtRNAMet D-domain with the T-domain was enhanced by a pseudouridine located in either the D or T-domains compared to that of the unmodified RNAs (Kd=25.3, 24.6 and 44.4 microM, respectively). Mg2+ also affected the folding interaction of a yeast mtRNALeu1, but had minimal effect on the folding of an Escherichia coli cytoplasmic tRNALeu. The D-domain modification, dihydrouridine, facilitated mtRNALeu folding. These data indicate that conserved modifications assist and stabilize the formation of the functional mtRNA tertiary structure.