Thermal adaptation of mRNA secondary structure: stability versus lability

2021 ◽  
Vol 118 (45) ◽  
pp. e2113324118
Author(s):  
Ming-ling Liao ◽  
Yun-wei Dong ◽  
George N. Somero

Macromolecular function commonly involves rapidly reversible alterations in three-dimensional structure (conformation). To allow these essential conformational changes, macromolecules must possess higher order structures that are appropriately balanced between rigidity and flexibility. Because of the low stabilization free energies (marginal stabilities) of macromolecule conformations, temperature changes have strong effects on conformation and, thereby, on function. As is well known for proteins, during evolution, temperature-adaptive changes in sequence foster retention of optimal marginal stability at a species’ normal physiological temperatures. Here, we extend this type of analysis to messenger RNAs (mRNAs), a class of macromolecules for which the stability–lability balance has not been elucidated. We employ in silico methods to determine secondary structures and estimate changes in free energy of folding (ΔGfold) for 25 orthologous mRNAs that encode the enzyme cytosolic malate dehydrogenase in marine mollusks with adaptation temperatures spanning an almost 60 °C range. The change in free energy that occurs during formation of the ensemble of mRNA secondary structures is significantly correlated with adaptation temperature: ΔGfold values are all negative and their absolute values increase with adaptation temperature. A principal mechanism underlying these adaptations is a significant increase in synonymous guanine + cytosine substitutions with increasing temperature. These findings open up an avenue of exploration in molecular evolution and raise interesting questions about the interaction between temperature-adaptive changes in mRNA sequence and in the proteins they encode.

2006 ◽  
Vol 7 (1) ◽  
pp. 37-43 ◽  
Author(s):  
T. A. Hughes ◽  
J. N. McElwaine

Secondary structures within the 5′ untranslated regions of messenger RNAs can have profound effects on the efficiency of translation of their messages and thereby on gene expression. Consequently they can act as important regulatory motifs in both physiological and pathological settings. Current approaches to predicting the secondary structure of these RNA sequences find the structure with the global-minimum free energy. However, since RNA folds progressively from the 5′ end when synthesised or released from the translational machinery, this may not be the most probable structure. We discuss secondary structure prediction based on local-minimisation of free energy with thermodynamic fluctuations as nucleotides are added to the 3′ end and show that these can result in different secondary structures. We also discuss approaches for studying the extent of the translational inhibition specified by structures within the 5′ untranslated region.


Author(s):  
Balaji Selvam ◽  
Ya-Chi Yu ◽  
Liqing Chen ◽  
Diwakar Shukla

<p>The SWEET family belongs to a class of transporters in plants that undergoes large conformational changes to facilitate transport of sugar molecules across the cell membrane. However, the structures of their functionally relevant conformational states in the transport cycle have not been reported. In this study, we have characterized the conformational dynamics and complete transport cycle of glucose in OsSWEET2b transporter using extensive molecular dynamics simulations. Using Markov state models, we estimated the free energy barrier associated with different states as well as 1 for the glucose the transport mechanism. SWEETs undergoes structural transition to outward-facing (OF), Occluded (OC) and inward-facing (IF) and strongly support alternate access transport mechanism. The glucose diffuses freely from outside to inside the cell without causing major conformational changes which means that the conformations of glucose unbound and bound snapshots are exactly same for OF, OC and IF states. We identified a network of hydrophobic core residues at the center of the transporter that restricts the glucose entry to the cytoplasmic side and act as an intracellular hydrophobic gate. The mechanistic predictions from molecular dynamics simulations are validated using site-directed mutagenesis experiments. Our simulation also revealed hourglass like intermediate states making the pore radius narrower at the center. This work provides new fundamental insights into how substrate-transporter interactions actively change the free energy landscape of the transport cycle to facilitate enhanced transport activity.</p>


2019 ◽  
Vol 26 (10) ◽  
pp. 743-750 ◽  
Author(s):  
Remya Radha ◽  
Sathyanarayana N. Gummadi

Background:pH is one of the decisive macromolecular properties of proteins that significantly affects enzyme structure, stability and reaction rate. Change in pH may protonate or deprotonate the side group of aminoacid residues in the protein, thereby resulting in changes in chemical and structural features. Hence studies on the kinetics of enzyme deactivation by pH are important for assessing the bio-functionality of industrial enzymes. L-asparaginase is one such important enzyme that has potent applications in cancer therapy and food industry.Objective:The objective of the study is to understand and analyze the influence of pH on deactivation and stability of Vibrio cholerae L-asparaginase.Methods:Kinetic studies were conducted to analyze the effect of pH on stability and deactivation of Vibrio cholerae L-asparaginase. Circular Dichroism (CD) and Differential Scanning Calorimetry (DSC) studies have been carried out to understand the pH-dependent conformational changes in the secondary structure of V. cholerae L-asparaginase.Results:The enzyme was found to be least stable at extreme acidic conditions (pH< 4.5) and exhibited a gradual increase in melting temperature from 40 to 81 °C within pH range of 4.0 to 7.0. Thermodynamic properties of protein were estimated and at pH 7.0 the protein exhibited ΔG37of 26.31 kcal mole-1, ΔH of 204.27 kcal mole-1 and ΔS of 574.06 cal mole-1 K-1.Conclusion:The stability and thermodynamic analysis revealed that V. cholerae L-asparaginase was highly stable over a wide range of pH, with the highest stability in the pH range of 5.0–7.0.


2021 ◽  
Vol 18 (1) ◽  
pp. 13-19
Author(s):  
A.I. Rabadanova ◽  

The steady growth in the number of drug addicts, especially among young people, dictates the need to find ways to prevent and treat this disease. In this regard, there is a need for a more detailed study of the mechanisms of the course of this disease using modern research methods, such as atomic force microscopy and fluorescence analysis of amino acid residues. Purpose of the work: to reveal the structural and functional state of erythrocyte membranes in drug addiction. Materials and methods. The studies were carried out on the erythrocyte membranes of 60 subjects suffering from heroin addiction. The shape and topography of the erythrocyte surface were studied, and spectral analysis of the proteins of the erythrocyte membranes was carried out. Results. The conducted AFM studies of erythrocyte membranes indicate the heterogeneity of the surface mechanical properties of the erythrocyte membranes of drug addicts. The data obtained indicate an acceleration of the aging process of erythrocytes in drug addiction, which goes in two ways: the formation of outgrowths on the plasmolemma, which subsequently die off (echinocytes) and invagination of the plasmolemma of erythrocytes (spherocytes). The fluorescence spectrum of amino acids in erythrocytes of drug addicts is characterized by a significant decrease in the intensity of almost all peaks and a shift of the fluorescence peak to the short-wave region. Findings. With drug addiction, changes in the structural integrity of red blood cells are noted. In people with drug addiction, in comparison with healthy people, there is a higher variability of the morphology of erythrocytes, which is expressed in a significant increase in the proportion of echinocytes and spherocytes against the background of a significant decrease in the number of discocytes. For the membrane proteins of erythrocytes of drug addicts, conformational changes are characteristic, manifested in a decrease in the intensity of fluorescence of aromatic amino acids, which indicates their structural modification and significant vulnerability of the hematopoietic system. They are largely determined by changes in the fluorescence intensity of tryptophan and, to a lesser extent, tyrosine, which indicates the preservation of the three-dimensional structure of the protein.


2017 ◽  
Vol 114 (21) ◽  
pp. E4158-E4167 ◽  
Author(s):  
Bogdan Lev ◽  
Samuel Murail ◽  
Frédéric Poitevin ◽  
Brett A. Cromer ◽  
Marc Baaden ◽  
...  

Pentameric ligand-gated ion channels control synaptic neurotransmission by converting chemical signals into electrical signals. Agonist binding leads to rapid signal transduction via an allosteric mechanism, where global protein conformational changes open a pore across the nerve cell membrane. We use all-atom molecular dynamics with a swarm-based string method to solve for the minimum free-energy gating pathways of the proton-activated bacterial GLIC channel. We describe stable wetted/open and dewetted/closed states, and uncover conformational changes in the agonist-binding extracellular domain, ion-conducting transmembrane domain, and gating interface that control communication between these domains. Transition analysis is used to compute free-energy surfaces that suggest allosteric pathways; stabilization with pH; and intermediates, including states that facilitate channel closing in the presence of an agonist. We describe a switching mechanism that senses proton binding by marked reorganization of subunit interface, altering the packing of β-sheets to induce changes that lead to asynchronous pore-lining M2 helix movements. These results provide molecular details of GLIC gating and insight into the allosteric mechanisms for the superfamily of pentameric ligand-gated channels.


2017 ◽  
Author(s):  
Mizuki Takemoto ◽  
Yongchan Lee ◽  
Ryuichiro Ishitani ◽  
Osamu Nureki

AbstractSecondary active transporters translocate their substrates using the electrochemical potentials of other chemicals, undergoing large-scale conformational changes. Despite extensive structural studies, the atomic details of the transport mechanism still remain elusive. Here we performed a series of all-atom molecular dynamics simulations of the triose-phosphate/phosphate translocator (TPT), which exports organic phosphates in the chloroplast stroma in strict counter exchange with inorganic phosphate (Pi). Biased sampling methods, including string method and umbrella sampling, successfully reproduced the conformational changes between the inward– and outward-facing states, along with the substrate binding. The free energy landscape of this entire TPT transition pathway demonstrated the alternating access and substrate translocation mechanisms, which revealed Pi is relayed by positively charged residues along the transition pathway. Furthermore, the conserved Glu207 functions as a “molecular switch”, linking the local substrate binding and the global conformational transition. Our results provide atomic-detailed insights into the energy coupling mechanism of antiporter.


Author(s):  
Michela Zuffo ◽  
Aurélie Gandolfini ◽  
Brahim Heddi ◽  
Anton Granzhan

ABSTRACTDNA is polymorphic since, despite its ubiquitous presence as a double-stranded helix, it is able to fold into a plethora of other secondary structures both in vitro and in cells. Despite the considerable advances that have been made in understanding this structural diversity, its high-throughput investigation still faces severe limitations. This mainly stems from the lack of suitable label-free methods, combining a fast and cheap experimental workflow with high information content. Here, we explore the use of intrinsic fluorescence emitted by nucleic acids for this scope. After a preliminary assessment of the suitability of this phenomenon for tracking the conformational changes of DNA, we examined the intrinsic steady-state emission spectra of an 89-membered set of synthetic oligonucleotides with reported conformation (G-quadruplexes, i-motifs, single- and double stranded DNA) by means of multivariate analysis. Specifically, principal component analysis of emission spectra resulted in successful clustering of oligonucleotides into three corresponding conformational groups, albeit without discrimination between single- and double-stranded structures. Linear discriminant analysis of the same training set was exploited for the assessment of new sequences, allowing the evaluation of their G4-forming propensity. Our method does not require any labelling agent or dye, avoiding the related intrinsic bias, and can be utilized to screen novel sequences of interest in a high-throughput and cost-effective manner. In addition, we observed that left-handed (Z-) G4 structures were systematically more fluorescent than most other G4 structures, almost reaching the quantum yield of 5′-d[(G3T)3G3]-3′ (G3T), the most fluorescent G4 structure reported to date. This property is likely to arise from the similar base-stacking geometry in both types of structures.


2019 ◽  
Vol 20 (1) ◽  
Author(s):  
Guangyao Zhou ◽  
Jackson Loper ◽  
Stuart Geman

Abstract Background A folding RNA molecule encounters multiple opportunities to form non-native yet energetically favorable pairings of nucleotide sequences. Given this forbidding free-energy landscape, mechanisms have evolved that contribute to a directed and efficient folding process, including catalytic proteins and error-detecting chaperones. Among structural RNA molecules we make a distinction between “bound” molecules, which are active as part of ribonucleoprotein (RNP) complexes, and “unbound,” with physiological functions performed without necessarily being bound in RNP complexes. We hypothesized that unbound molecules, lacking the partnering structure of a protein, would be more vulnerable than bound molecules to kinetic traps that compete with native stem structures. We defined an “ambiguity index”—a normalized function of the primary and secondary structure of an individual molecule that measures the number of kinetic traps available to nucleotide sequences that are paired in the native structure, presuming that unbound molecules would have lower indexes. The ambiguity index depends on the purported secondary structure, and was computed under both the comparative (“gold standard”) and an equilibrium-based prediction which approximates the minimum free energy (MFE) structure. Arguing that kinetically accessible metastable structures might be more biologically relevant than thermodynamic equilibrium structures, we also hypothesized that MFE-derived ambiguities would be less effective in separating bound and unbound molecules. Results We have introduced an intuitive and easily computed function of primary and secondary structures that measures the availability of complementary sequences that could disrupt the formation of native stems on a given molecule—an ambiguity index. Using comparative secondary structures, the ambiguity index is systematically smaller among unbound than bound molecules, as expected. Furthermore, the effect is lost when the presumably more accurate comparative structure is replaced instead by the MFE structure. Conclusions A statistical analysis of the relationship between the primary and secondary structures of non-coding RNA molecules suggests that stem-disrupting kinetic traps are substantially less prevalent in molecules not participating in RNP complexes. In that this distinction is apparent under the comparative but not the MFE secondary structure, the results highlight a possible deficiency in structure predictions when based upon assumptions of thermodynamic equilibrium.


Sign in / Sign up

Export Citation Format

Share Document