Molecular modeling of phytochrome

Application of molecular dynamics simulation for exploring the roles of plant biomolecules in promoting environmental health

Understanding the dynamic changes of plant biomolecules is vital for exploring their mechanisms in the environment. Molecular dynamics (MD) simulation has been widely used to study structural evolution and corresponding properties of plant biomolecules at the microscopic scale. Here, this review (i) outlines structural properties of plant biomolecules, and the crucial role of MD simulation in advancing studies of the biomolecules; (ii) describes the development of MD simulation in plant biomolecules, determinants of simulation, and analysis parameters; (iii) introduces the applications of MD simulation in plant biomolecules, including the response of the biomolecules to multiple stresses, their roles in corrosive environments, and their contributions in improving environmental health; (iv) reviews techniques integrated with MD simulation, such as molecular biology, quantum mechanics, molecular docking, and machine learning modeling, which bridge gaps in MD simulation. Finally, we make suggestions on determination of force field types, investigation of plant biomolecule mechanisms, and use of MD simulation in combination with other techniques. This review provides comprehensive summaries of the mechanisms of plant biomolecules in the environment revealed by MD simulation and validates it as an applicable tool for bridging gaps between macroscopic and microscopic behavior, providing insights into the wide application of MD simulation in plant biomolecules.

Arginine of retinoic acid receptor beta which coordinates with the carboxyl group of retinoic acid functions independent of the amino acid residues responsible for retinoic acid receptor subtype ligand specificity

The biological actions of retinoic acid (RA) are mediated by retinoic acid receptors (RARalpha, RARbeta, and RARgamma) and retinoid X receptors (RXRalpha, RXRbeta, and RXRgamma). Consistent with the X-ray crystal structures of RARalpha and RARgamma, site-directed mutagenesis studies have demonstrated the importance of a conserved Arg residue (alphaArg(276), betaArg(269), and gammaArg(278)) for coordination with the carboxyl group of RA. However, mutation of Arg(269) to Ala in RARbeta causes only a 3- to 6-fold increase in the K(d) for RA and EC(50) in RA-dependent transcriptional transactivation assays while the homologous mutation in either RARalpha or RARgamma causes a 110-fold and a 45-fold increase in EC(50) value, respectively. To further investigate the nature of this difference, we prepared mutant RARs to determine the effect of conversion of betaR269A to a mutant which mimics either RARalpha ligand selectivity (betaA225S/R269A) or RARgamma ligand selectivity (betaI263M/R269A/V338A). Our results demonstrate that in RARbeta mutants that acquire either RARalpha or RARgamma ligand specificity the Arg(269) position responsible for coordination with the carboxyl group of retinoids continued to function like that of RARbeta. Furthermore, three mutant receptors (betaA225S/R269A, betaA225S/F279, and alphaF286A) were found to have a greater than wild-type affinity for the RARalpha-selective ligand Am580. Finally, a homology-based computer model of the ligand binding domain (LBD) of RARbeta and the X-ray crystal structures of the LBD of both RARalpha and RARgamma are used to describe potential mechanisms responsible for the increased affinity of some mutants for Am580 and for the difference in the effect of mutation of Arg(269) in RARbeta compared to its homologous Arg in RARalpha and RARgamma.

Molecular modeling of archaebacterial bipolar tetraether lipid membranes

Membranes composed of glycerol dialkylnonitol tetraether (GDNT) lipids from the thermoacidophilic archaebacterium Sulfolobus acidocaldarius have been studied by molecular modeling. GDNT membranes containing eight cyclopentane rings in the molecule are packed much tighter than those without rings. When containing eight cyclopentane rings, the beta-D-galactosyl-D-glucose head-group of GDNT runs almost parallel to the membrane surface. However, when containing no rings, the head-group is oriented perpendicular to the membrane surface. Using molecular dynamics calculations, we have also conducted comparative studies of membrane packing between GDNT and various non-archaebacterial membranes. Compared to gel state dipalmitoylphosphatidylcholine (DPPC) and gel state distearoylphosphatidylcholine (DSPC) bilayers, the GDNT membrane with eight cyclopentane rings has a more negative interaction energy, thus a tighter membrane packing, while the GDNT without rings is less tightly packed than gel state DSPC. Based on the calculated interaction energies, the GDNT membranes (with and without rings) are much more tightly packed than DPhPC (an ester-linked diphytanyl PC) and DPhyPC (an ether-linked diphytanyl PC) bilayers. This suggests that the branched methyl group in the phytanyl chain is not the major contributor of the tight packing of GDNT membranes. The biological implication of this study is that the cyclopentane ring could increase GDNT membrane thermal stability. This explains why the number of cyclopentane rings in archaebacterial lipid increases with increasing growth temperature. Perhaps, through the ring-temperature compensation mechanism the plasma membrane of thermoacidophilic archaebacteria is able to maintain a tight and rigid structure, consequently, a constant proton gradient between the extracellular (pH 2.5) and intracellular compartment (pH 6.5), over a wide range of growth temperatures.

New insights into an old protein: the functional diversity of mammalian glyceraldehyde-3-phosphate dehydrogenase

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was considered a classical glycolytic protein examined for its pivotal role in energy production. It was also used as a model protein for analysis of protein structure and enzyme mechanisms. The GAPDH gene was utilized as a prototype for studies of genetic organization, expression and regulation. However, recent evidence demonstrates that mammalian GAPDH displays a number of diverse activities unrelated to its glycolytic function. These include its role in membrane fusion, microtubule bundling, phosphotransferase activity, nuclear RNA export, DNA replication and DNA repair. These new activities may be related to the subcellular localization and oligomeric structure of GAPDH in vivo. Furthermore, other investigations suggest that GAPDH is involved in apoptosis, age-related neurodegenerative disease, prostate cancer and viral pathogenesis. Intriguingly, GAPDH is also a unique target of nitric oxide. This review discusses the functional diversity of GAPDH in relation to its protein structure. The mechanisms through which mammalian cells may utilize GAPDH amino acid sequences to provide these new functions and to determine its intracellular localization are considered. The interrelationship between new GAPDH activities and its role in cell pathologies is addressed.

Computer analysis of phytochrome sequences and reevaluation of the phytochrome secondary structure by Fourier transform infrared spectroscopy

A repertoire of various methods of computer sequence analysis was applied to phytochromes in order to gain new insights into their structure and function. A statistical analysis of 23 complete phytochrome sequences revealed regions of non-random amino acid composition, which are supposed to be of particular structural or functional importance. All phytochromes other than phyD and phyE from Arabidopsis have at least one such region at the N-terminus between residues 2 and 35. A sequence similarity search of current databases indicated striking homologies between all phytochromes and a hypothetical 84.2-kDa protein from the cyanobacterium Synechocystis. Furthermore, scanning the phytochrome sequences for the occurrence of patterns defined in the PROSITE database detected the signature of the WD repeats of the beta-transducin family within the functionally important 623-779 region (sequence numbering of phyA from Avena) in a number of phytochromes. A multiple sequence alignment performed with 23 complete phytochrome sequences is made available via the IMB Jena World-Wide Web server (http://www.imb-jena.de/PHYTO.html). It can be used as a working tool for future theoretical and experimental studies. Based on the multiple alignment striking sequence differences between phytochromes A and B were detected directly at the N-terminal end, where all phytochromes B have an additional stretch of 15-42 amino acids. There is also a variety of positions with totally conserved but different amino acids in phytochromes A and B. Most of these changes are found in the sequence segment 150-200. It is, therefore, suggested that this region might be of importance in determining the photosensory specificity of the two phytochromes. The secondary structure prediction based on the multiple alignment resulted in a small but significant beta-sheet content. This finding is confirmed by a reevaluation of the secondary structure using FTIR spectroscopy.

Biogenesis of thylakoid membranes with emphasis on the process in Chlamydomonas

Recent results obtained by electron microscopic and biochemical analyses of greening Chlamydomonas reinhardtii y1 suggest that localized expansion of the plastid envelope is involved in thylakoid biogenesis. Kinetic analyses of the assembly of light-harvesting complexes and development of photosynthetic function when degreened cells of the alga are exposed to light suggest that proteins integrate into membrane at the level of the envelope. Current information, therefore, supports the earlier conclussion that the chloroplast envelope is a major biogenic structure, from which thylakoid membranes emerge. Chloroplast development in Chlamydomonas provides unique opportunities to examine in detail the biogenesis of thylakoids.