Thermostability of the PYL–PP2C Heterodimer Is Dependent on Magnesium: In Silico Insights into the Link between Heat Stress Response and Magnesium Deficiency in Plants

Modeling and dynamical analysis of the full-length structure of factor XII with zinc

Zinc (II), the second most abundant transition metal in blood, binds to the initiator of the contact pathway, factor XII (FXII). This binding induces conformational changes in the structure of FXII eventually leading to its activation. Despite many in vitro and in vivo studies on zinc-mediated activation of FXII, its molecular mechanism remains elusive mainly due to absence of a full-length structural model of FXII. To this end, this study investigated the role of zinc in the structure and dynamics of the full-length structure FXII that was obtained through molecular modeling. We have used four structural templates covering more than 70% of the FXII sequence and the remaining interconnecting regions were built by loop modeling. The resulting full-length structure of FXII contained disordered regions, but in comparison to the AlphaFold (AF) prediction, our full-length model represented a more realistic structure because of the disordered regions which were modeled to yield a more compact full-length structure in our model than the AF structure. Other than the disordered regions, our model and AF prediction were highly similar. The resulting full-length FXII structure was used to generate different systems representing the zinc-bound form (holo). Further to assess the contribution of the disulfide bridges, we also analyzed the apo and holo FXII structures with oxidized or reduced cysteine side-chains. Simulations suggested zinc binding conferred rigidity to the structure, particularly to the light chain of FXII. Zinc binding alone was sufficient to limit the backbone flexibility while 15 disulfide bonds, which were scattered throughout the structure, made a less significant contribution to the backbone rigidity. Altogether our results provide insights into the first realistic full-length structure of FXII focusing on the impact of structural zinc and disulfide bridges in the dynamics of this structure.

A comprehensive transcriptome analysis of contrasting rice cultivars highlights the role of auxin and ABA responsive genes in heat stress response

Sensing a change in ambient temperature is key to survival among all living organisms. Temperature fluctuations due to climate change are a matter of grave concern since it adversely affects growth and eventually the yield of crop plants, including two of the major cereals, i.e., rice and wheat. Thus, to understand the response of rice seedlings to elevated temperatures, we performed microarray-based transcriptome analysis of two contrasting rice cultivars, Annapurna (heat tolerant) and IR64 (heat susceptible), by subjecting their seedlings to 37 °C and 42 °C, sequentially. The transcriptome analyses revealed a set of uniquely regulated genes and related pathways in red rice cultivar Annapurna, particularly associated with auxin and ABA as a part of heat stress response in rice. The changes in expression of few auxin and ABA associated genes, such as OsIAA13, OsIAA20, ILL8, OsbZIP12, OsPP2C51, OsDi19-1 and OsHOX24, among others, were validated under high-temperature conditions using RT-qPCR. In particular, the expression of auxin-inducible SAUR genes was enhanced considerably at both elevated temperatures. Further, using genes that expressed inversely under heat vs. cold temperature conditions, we built a regulatory network between transcription factors (TF) such as HSFs, NAC, WRKYs, bHLHs or bZIPs and their target gene pairs and determined regulatory coordination in their expression under varying temperature conditions. Our work thus provides useful insights into temperature-responsive genes, particularly under elevated temperature conditions, and could serve as a resource of candidate genes associated with thermotolerance or downstream components of temperature sensors in rice.

Going to the roots of reduced magnesium dietary intake: A tradeoff between climate changes and sources

Magnesium is essential in plants where it is associated with chlorophyll pigments and serves as a cofactor of enzymes implicated in photosynthesis and metabolism. It is an essential nutrient for animals, involved in hundreds metabolic reaction and crucial for the biological activity of ATP. Not surprisingly, magnesium deficiency is detrimental for the health of plants and animals. In humans, subclinical magnesium deficiency is common and generates chronic inflammation, which is the common denominator of a wide range of mental and physical health problems from metabolic diseases to cognitive impairment, from osteopenia and sarcopenia to depression. It is ascertained that magnesium content in fruits and vegetables dropped in the last fifty years, and about 80% of this metal is lost during food processing. As a consequence, a large percentage of people all over the world does not meet the minimum daily magnesium requirement. In this scoping review, we summarize how agronomic and environmental factors, including global warming, affect magnesium content and availability in the soil and, consequently, in the food chain, with the aim of attracting the interest of botanists, agronomists, animal and human nutritionists and physicians to work on a strategy that grants adequate magnesium intake for everybody.

Genome-wide phylogenetic and structural analysis reveals the molecular evolution of the ABA receptor gene family

The plant hormone abscisic acid (ABA) plays a crucial role during the plant life cycle as well as in adaptive responses to environmental stresses. The core regulatory components of ABA signaling in plants are the pyrabactin resistance1/PYR1-like/regulatory component of ABA receptor family (PYLs), which comprise the largest plant hormone receptor family known. They act as negative regulators of members of the protein phosphatase type 2C family. Due to the biological importance of PYLs, many researchers have focused on their genetic redundancy and consequent functional divergence. However, little is understood of their evolution and its impact on the generation of regulatory diversity. In this study, we identify positive selection and functional divergence in PYLs through phylogenetic reconstruction, gene structure and expression pattern analysis, positive selection analysis, functional divergence analysis, and structure comparison. We found the correlation of desensitization of PYLs under specific modifications in the molecular recognition domain with functional diversification. Hence, an interesting antagonistic co-evolutionary mechanism is proposed for the functional diversification of ABA receptor family proteins. We believe a compensatory evolutionary pathway may have occurred.

Molecular Dynamics Simulations Reveal Differentiated Context-Dependent Conformational Dynamics of Two Proteins of the Same Family

The Arabidopsis pyrabactin resistant 1 (PYR1)-like family of proteins (PYLs) are receptors of abscisic acid (ABA), an essential small signaling molecule in plants. Here, we report a comparative molecular dynamics (MD) study on two PYL members, PYR1 and PYL10, which, despite their highly similar sequences and structures, have been suggested to belong to two different subclasses of PYLs, one being dimeric and relying on binding to ABA to inhibit downstream type 2C protein phosphatases (PP2Cs) and the other being monomeric and able to constitutively inhibit downstream PP2Cs without ABA. MD simulations have been carried out on these proteins in various monomeric or complexation states. Analyses of the simulations unambiguously confirm that ABA has large effects on the conformational dynamics of PYR1 but not PYL10, whereas a downstream PP2C has much larger effects on PYL10 than on PYR1. The differentiated effects are consistent with the functional differences between the two proteins. Potential of mean forces (PMFs) calculated by umbrella sampling showed that binding to ABA strengthens the PYR1-PP2C complex, increasing the PMF change for dissociation from 7.5 to 12.0 kcal mol^-1. On the other hand, the same PMF change for an apo-PYL10-PP2C complex was computed to be 9.5 kcal mol^-1, suggesting stronger binding in apo-PYL10-PP2C than in apo-PYR1-PP2C. Several specific sequence features that may contribute to the functional differentiation between PYR1 and PYL10 are suggested based on the intersubunit residue-residue contacts occurred in the simulations.