Most mitochondrial proteins are synthesized in the cytosol as precursor and imported into the mitochondria by Tom complexes (translocase of outer membrane complexes). Knowledge of the binding mechanism between precursor and Tom20 in plants is very limited. Here, computational methods are employed to improve our understanding of the interactions between both molecules. To this end, we model mitochondrial superoxide dismutase precursor (pSOD) in complex with Tom20 in Oryza sativa (OsTom20). In a first stage, five main binding modes were generated using clustering analysis, energy minimization, and expert knowledge. In a second stage, the quality and validity of the resulting complexes is assessed by molecular dynamics (MD) simulations with a generalized Born solvation model. The change in binding free energies is estimated using a computational alanine scanning technique. We identified a particularly favorable complex between pSOD and OsTom20, exhibiting the lowest binding free energy among all candidates and correlating well with experimental data. Furthermore, three independent explicit solvent MD simulations of this structure, each of 100 ns duration, reveal that hydrophobic interactions occur between pSOD and OsTom20, in particular between L158 of pSOD and W81 of OsTom20, as evidenced by analysis of intermolecular distances and corresponding relative free energy landscapes. L158 is part of an interacting LRTLA motif. These results provide new insight into the structural basis and dynamics of precursor recognition by Tom20 in plant, and their generality is supported by sequence alignments with seven other plants.

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.