The Molecular Recognition Mechanism for Superoxide Dismutase Presequence Binding to the Mitochondrial Protein Import Receptor Tom20 from Oryza sativa Involves an LRTLA Motif

Discovery of the new alpha-glucosidase inhibitor with therapeutic potential in type 2 diabetes mellitus by a novel high-throughput virtual screening and free energy evaluation

Type 2 diabetes can cause a variety of complications, significantly affecting people's health. Given their ability to suppress carbohydrate digestion, alpha-glucosidase inhibitors are effective treatments for diabetes. However, the current approved glucosidase inhibitors' side effects of abdominal discomfort limit their use. We used the compound Pg3R from the natural fruit berry as a reference, screening against a large database of 22 million compounds to identify potential health-friendly alpha-glucosidase inhibitors. Ligand-based screening enables us to identify 3968 ligands that exhibit structural similarity compared to the natural compound. These lead hits were used for LeDock, and their binding free energies were evaluated by MM/GBSA. Among the top-scoring candidates, ZINC263584304 exhibited the strongest binding affinity to alpha-glucosidase, with a "low-fat" structural characteristic. Its recognition mechanism was further investigated by microsecond MD simulations and free energy landscapes, exhibiting novel conformational changes during the binding process. Our study provided a novel alpha-glucosidase inhibitor with the potential to treat type 2 diabetes.

Molecular Insights into Substrate Binding of the Outer Membrane Enzyme OmpT

The enzyme OmpT of the outer membrane of Escherichia coli shows proteolytic activity and cleaves peptides and proteins. Using molecular dynamics simulations in a fully hydrated lipid bilayer on a time scale of hundreds of nanoseconds, we draw a detailed atomic picture of substrate recognition in the OmpT-holo enzyme complex. Hydrogen bonds and salt bridges are essential for maintaining the integrity of the active site and play a central role for OmpT in recognizing its substrate. Electrostatic interactions are critical at all stages from approaching the substrate to docking at the active site. Computational alanine scanning based on the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) approach confirms the importance of multiple residues in the active site that form salt bridges. The substrate fluctuates along the axis of the β-barrel, which is associated with oscillations of the binding cleft formed by the residue pairs D210-H212 and D83-D85. Principal component analysis suggests that substrate and protein movements are correlated. We observe the transient presence of putative catalytic water molecules near the active site, which may be involved in the nucleophilic attack on the cleavable peptide bond of the substrate.

SHP-1 knockdown suppresses mitochondrial biogenesis and aggravates mitochondria-dependent apoptosis induced by all trans retinal through the STING/AMPK pathways

BACKGROUND

Oxidative stress-caused damage to the retinal pigment epithelium (RPE) underlies the onset and progression of age-related macular degeneration (AMD). Impaired mitochondrial biogenesis sensitizes RPE cells to mitochondrial dysfunction, energy insufficiency and death. Src-homology 2 domain-containing phosphatase (SHP)-1 is important in regulating immune responses and cell survival. However, its roles in cell survival are not always consistent. Until now, the effects of SHP-1 on RPE dysfunction, especially mitochondrial homeostasis, remain to be elucidated. We sought to clarify the effects of SHP-1 in RPE cells in response to atRAL-induced oxidative stress and determine the regulatory mechanisms involved.

METHODS

In the all trans retinal (atRAL)-induced oxidative stress model, we used the vector of lentivirus to knockdown the expression of SHP-1 in ARPE-19 cells. CCK-8 assay, Annexin V/PI staining and JC-1 staining were utilized to determine the cell viability, cell apoptosis and mitochondrial membrane potential. We also used immunoprecipitation to examine the ubiquitination modification of stimulator of interferon genes (STING) and its interaction with SHP-1. The expression levels of mitochondrial marker, proteins related to mitochondrial biogenesis, and signaling molecules involved were examined by western blotting analysis.

RESULTS

We found that SHP-1 knockdown predisposed RPE cells to apoptosis, aggravated mitochondrial damage, and repressed mitochondrial biogenesis after treatment with atRAL. Immunofluoresent staining and immunoprecipitation analysis confirmed that SHP-1 interacted with the endoplasmic reticulum-resident STING and suppressed K63-linked ubiquitination and activation of STING. Inhibition of STING with the specific antagonist H151 attenuated the effects of SHP-1 knockdown on mitochondrial biogenesis and oxidative damage. The adenosine monophosphate-activated protein kinase (AMPK) pathway acted as the crucial downstream target of STING and was involved in the regulatory processes.

CONCLUSIONS

These findings suggest that SHP-1 knockdown potentiates STING overactivation and represses mitochondrial biogenesis and cell survival, at least in part by blocking the AMPK pathway in RPE cells. Therefore, restoring mitochondrial health by regulating SHP-1 in RPE cells may be a potential therapeutic strategy for degenerative retinal diseases including AMD.

Mechanism of zinc ejection by disulfiram in nonstructural protein 5A

Hepatitis C virus (HCV) is a notorious member of the Flaviviridae family of enveloped, positive-strand RNA viruses. Non-structural protein 5A (NS5A) plays a key role in HCV replication and assembly. NS5A is a multi-domain protein which includes an N-terminal amphipathic membrane anchoring alpha helix, a highly structured domain-1, and two intrinsically disordered domains 2-3. The highly structured domain-1 contains a zinc finger (Zf)-site, and binding of zinc stabilizes the overall structure, while ejection of this zinc from the Zf-site destabilizes the overall structure. Therefore, NS5A is an attractive target for anti-HCV therapy by disulfiram, through ejection of zinc from the Zf-site. However, the zinc ejection mechanism is poorly understood. To disclose this mechanism based on three different states, A-state (NS5A protein), B-state (NS5A + Zn), and C-state (NS5A + Zn + disulfiram), we have performed molecular dynamics (MD) simulation in tandem with DFT calculations in the current study. The MD results indicate that disulfiram triggers Zn ejection from the Zf-site predominantly through altering the overall conformation ensemble. On the other hand, the DFT assessment demonstrates that the Zn adopts a tetrahedral configuration at the Zf-site with four Cys residues, which indicates a stable protein structure morphology. Disulfiram binding induces major conformational changes at the Zf-site, introduces new interactions of Cys39 with disulfiram, and further weakens the interaction of this residue with Zn, causing ejection of zinc from the Zf-site. The proposed mechanism elucidates the therapeutic potential of disulfiram and offers theoretical guidance for the advancement of drug candidates.

High-throughput virtual screening and microsecond MD simulations to identify potential sugar mimic of the solute-binding protein BlAXBP of the ABC transporter from Bifidobacterium animalis subsp. Lactis

Xylotetraose is a prebiotic oligosaccharide can be utilized by the ABC transporter of the gut microbiota Bifidobacteria. BlAXBP is the solute binding protein of the ABC transporter, and its complex with xylotetraose has been solved by X-ray crystallography. Here, we have identified novel sugar mimic of BlAXBP by applying a high-throughput virtual screening of ZINC database containing a huge library with ∼22 M compounds. To begin with, we identified 18,571 ligands by a ligand-based virtual screening. Further, a total of 3968 compounds were selected for molecular docking due to their Tanimoto coefficient's value were larger than a cutoff of 0.08. The molecular mechanics-generalized born surface area was used to evaluate the binding free energies, and the top 10 ligands with free energies below an energy threshold of -35.22 kcal/mol were selected. ZINC13783511 formed the most stable complex with BlAXBP and its recognition mechanism were further explored by microsecond MD simulations in explicit solvent. Free energy landscapes were used to evaluate conformational changes of BlAXBP in its ligand free and binding states. Collectively, this work identified potential novel sugar mimics to BlAXBP, providing novel atomic-level understanding of the binding mechanism.

Recognition motifs rather than phylogenetic origin influence the ability of targeting peptides to import nuclear-encoded recombinant proteins into rice mitochondria

Mitochondria fulfil essential functions in respiration and metabolism as well as regulating stress responses and apoptosis. Most native mitochondrial proteins are encoded by nuclear genes and are imported into mitochondria via one of several receptors that recognize N-terminal signal peptides. The targeting of recombinant proteins to mitochondria therefore requires the presence of an appropriate N-terminal peptide, but little is known about mitochondrial import in monocotyledonous plants such as rice (Oryza sativa). To gain insight into this phenomenon, we targeted nuclear-encoded enhanced green fluorescent protein (eGFP) to rice mitochondria using six mitochondrial pre-sequences with diverse phylogenetic origins, and investigated their effectiveness by immunoblot analysis as well as confocal and electron microscopy. We found that the ATPA and COX4 (Saccharomyces cerevisiae), SU9 (Neurospora crassa), pFA (Arabidopsis thaliana) and OsSCSb (Oryza sativa) peptides successfully directed most of the eGFP to the mitochondria, whereas the MTS2 peptide (Nicotiana plumbaginifolia) showed little or no evidence of targeting ability even though it is a native plant sequence. Our data therefore indicate that the presence of particular recognition motifs may be required for mitochondrial targeting, whereas the phylogenetic origin of the pre-sequences probably does not play a key role in the success of mitochondrial targeting in dedifferentiated rice callus and plants.

EAE-induced upregulation of mitochondrial MnSOD is associated with increases of mitochondrial SGK1 and Tom20 protein in the mouse kidney cortex

Our previous demonstration that severe experimental autoimmune encephalomyelitis (EAE) increases MnSOD protein abundance in the mouse kidney cortex led this study to elucidate the underlying mechanism with monensin-treated HEK293 cells as a model. Severe EAE increases mitochondrial protein abundance of SGK1 kinase and Tom20, a critical subunit of mitochondrial translocase in the renal cortex. In HEK293 cells, catalase inhibits monensin-induced increases of mitochondrial SGK1 and Tom20 protein levels. Further, GSK650394, a specific inhibitor of SGK1 reduces monensin-induced increase of mitochondrial protein abundance of Tom20 and MnSOD. Finally, RNAi of Tom20 reduces the effect of monensin on MnSOD. MnSOD and Tom20 physically associate with each other. In conclusion, in HEK293 cells, mitochondrial reactive oxygen species increase protein abundance of mitochondrial SGK1, which leads to a rise of mitochondrial Tom20, resulting in importing MnSOD protein into the mitochondria. This could be a mechanism by which severe EAE up-regulates mitochondrial MnSOD in the kidney cortex.

Identification of AtHsp90.6 involved in early embryogenesis and its structure prediction by molecular dynamics simulations

Heat-shock protein of 90 kDa (Hsp90) is a key molecular chaperone involved in folding the synthesized protein and controlling protein quality. Conformational dynamics coupled to ATPase activity in N-terminal domain is essential for Hsp90's function. However, the relevant process is still largely unknown in plant Hsp90s, especially those required for plant embryogenesis which is inextricably tied up with human survival. Here, AtHsp90.6, a member of Hsp90 family in Arabidopsis, was firstly identified as a protein essential for embryogenesis. Thus we modelled AtHsp90.6 in its functionally closed 'lid-down' and open 'lid-up' states, exploring the nucleotide binding mechanism in these two states. Free energy landscape and electrostatic potential analysis revealed the switching mechanism between these two states. Collectively, this study quantitatively analysed the conformational changes of AtHsp90.6 bound to ATP or ADP. This result may help us understand the mechanism of action of AtHsp90.6 in future.

Proteomic analysis of stachyose contribution to the growth of Lactobacillus acidophilus CICC22162

Stachyose is a functional oligosaccharide, acting as a potential prebiotic for colonic fermentation.

Understanding the impact of Fc glycosylation on its conformational changes by molecular dynamics simulations and bioinformatics

N-linked glycosylation of Fc at N297 plays an important role in its effector function, aberrance of which would cause disease pathogenesis. Here, we performed all-atom molecular dynamics simulations to explore the effects of Fc glycosylation on its dynamics behaviors. Firstly, equilibrium simulations suggested that Fc deglycosylation was able to induce residual flexibility in its CH2 domain. Besides, the free energy landscape revealed three minimum energy wells in deglycosylated Fc, representing its "open", "semi-closed" and "closed" states. However, we could only observe the "open" state of glycosylated Fc. Supportively, principal component analysis emphasized the prominent motion of delyclosylated Fc and dynamically depicted how it changed from the "open" state to its "closed" state. Secondly, we studied the recognition mechanism of the Fc binding to its partners. Energy decomposition analysis identified key residues of Fc to recognize its two partners P13 and P34. Evidently, electrostatic potential surfaces showed that electrostatic attraction helped to stabilize the interaction between Fc and its partners. Also, relative binding free energies explained different binding affinities in Fc-P13 and Fc-P34. Collectively, these results together provided the structural basis for understanding conformational changes of deglycosylated Fc and the recognition mechanism of the Fc binding to its partners.

Molecular dynamics simulation and bioinformatics study on chloroplast stromal ridge complex from rice (Oryza sativa L.)

Background

Rice (Oryza sativa L.) is one of the most important cereal crops in the world and its yield is closely related to the photosynthesis efficiency. The chloroplast stromal ridge complex consisting of PsaC-PsaD-PsaE plays an important role in plant photosynthesis, which has been a subject of many studies. Till now, the recognition mechanism between PsaC and PsaD in rice is still not fully understood.

Results

Here, we present the interaction features of OsPsaC and OsPsaD by molecular dynamics simulations and bioinformatics. Firstly, we identified interacting residues in the OsPsaC-OsPsaD complex during simulations. Significantly, important hydrogen bonds were observed in residue pairs R19-E103, D47-K62, R53-E63, Y81-R20, Y81-R61 and L26-V105. Free energy calculations suggested two salt bridges R19-E103 and D47-K62 were essential to maintain the OsPsaC-OsPsaD interaction. Supportively, electrostatic potentials surfaces of OsPsaD exhibited electrostatic attraction helped to stabilize the residue pairs R19-E103 and D47-K62. In particular, the importance of R19 was further verified by two 500 ns CG-MD simulations. Secondly, this study compared the stromal ridge complex in rice with that in other organisms. Notably, alignments of amino acids showed these two salt bridges R19-E103 and D47-K62 also existed in other organisms. Electrostatic potentials surfaces and X-ray structural analysis strongly suggested the stromal ridge complex in other organisms adopted a similar and general recognition mechanism.

Conclusions

These results together provided structure basis and dynamics behavior to understand recognition and assembly of the stromal ridge complex in rice.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0877-0) contains supplementary material, which is available to authorized users.

Mitochondrial targeting of the Arabidopsis F1-ATPase γ-subunit via multiple compensatory and synergistic presequence motifs

The majority of mitochondrial proteins are encoded in the nuclear genome and imported into mitochondria posttranslationally from the cytosol. An N-terminal presequence functions as the signal for the import of mitochondrial proteins. However, the functional information in the presequence remains elusive. This study reports the identification of critical sequence motifs from the presequence of Arabidopsis thaliana F1-ATPase γ-subunit (pFAγ). pFAγ was divided into six 10-amino acid segments, designated P1 to P6 from the N to the C terminus, each of which was further divided into two 5-amino acid subdivisions. These P segments and their subdivisions were substituted with Ala residues and fused to green fluorescent protein (GFP). Protoplast targeting experiments using these GFP constructs revealed that pFAγ contains several functional sequence motifs that are dispersed throughout the presequence. The sequence motifs DQEEG (P4a) and VVRNR (P5b) were involved in translocation across the mitochondrial membranes. The sequence motifs IAARP (P2b) and IAAIR (P3a) participated in binding to mitochondria. The sequence motifs RLLPS (P2a) and SISTQ (P5a) assisted in pulling proteins into the matrix, and the sequence motif IAARP (P2b) functioned in Tom20-dependent import. In addition, these sequence motifs exhibit complex relationships, including synergistic functions. Thus, multiple sequence motifs dispersed throughout the presequence are proposed to function cooperatively during protein import into mitochondria.

Binding mechanism between Hsp90 and Sgt1 explored by homology modeling and molecular dynamics simulations in rice

The Hsp90 (for heat shock protein90) and the Sgt1 (for suppressor of the G2 allele of skp1) are widely distributed in animals, yeast, and plants. The former functions as molecular chaperon activating a series of client proteins, the latter functions as an adaptor protein participating in multiple biological processes such as immunity response through interactions with different protein complexes. In the present study, we have constructed a homology model of Hsp90-Sgt1 complex in rice based on a recently resolved structure from barley and Arabidopsis to explore its binding mechanisms and to understand the detailed interaction profile. A total of 20 ns explicit solvent molecular dynamics simulations combined with MM-GBSA computations and virtual alanine scanning were performed for the modeled complex. In the final structure, three strong salt bridges were found between OsHsp90 and OsSgt1, D217(OsHsp90)-K186(OsSgt1), D218(OsHsp90)-K237(OsSgt1) and K161(OsHsp90)-E239(OsSgt1). Besides, residue Y173 of OsSgt1 played a vital role in the interactions with OsHsp90, the detailed interactions were discussed. These results would help us understand the critical features determining the Hsp90-Sgt1 binding process.

Isolation of thylakoid membrane complexes from rice by a new double-strips BN/SDS-PAGE and bioinformatics prediction of stromal ridge subunits interaction

Thylakoid membrane complexes of rice (Oryza sativa L.) play crucial roles in growth and crop production. Understanding of protein interactions within the complex would provide new insights into photosynthesis. Here, a new "Double-Strips BN/SDS-PAGE" method was employed to separate thylakoid membrane complexes in order to increase the protein abundance on 2D-gels and to facilitate the identification of hydrophobic transmembrane proteins. A total of 58 protein spots could be observed and subunit constitution of these complexes exhibited on 2D-gels. The generality of this new approach was confirmed using thylakoid membrane from spinach (Spinacia oleracea) and pumpkin (Cucurita spp). Furthermore, the proteins separated from rice thylakoid membrane were identified by the mass spectrometry (MS). The stromal ridge proteins PsaD and PsaE were identified both in the holo- and core- PSI complexes of rice. Using molecular dynamics simulation to explore the recognition mechanism of these subunits, we showed that salt bridge interactions between residues R19 of PsaC and E168 of PasD as well as R75 of PsaC and E91 of PsaD played important roles in the stability of the complex. This stromal ridge subunits interaction was also supported by the subsequent analysis of the binding free energy, the intramolecular distances and the intramolecular energy.