Molecular Docking and Interaction Studies of Identified Abscisic Acid Receptors in Oryza sativa: An In-Silico Perspective on Comprehending Stress Tolerance Mechanisms

Background: Abiotic stresses affect plants in several ways and as such, phytohormones such as abscisic acid (ABA) play an important role in conferring tolerance towards these stresses. Hence, to comprehend the role of ABA and its interaction with receptors of the plants, a thorough investigation is essential. Aim: The current study aimed to identify the ABA receptors in Oryza sativa, to find the receptor that binds best with ABA and to examine the mutations present to help predict better binding of the receptors with ABA. Methods: Protein sequences of twelve PYL (Pyrabactin resistance 1) and seven PP2C (type 2C protein phosphatase) receptors were retrieved from the Rice Annotation Project database and their 3D structures were predicted using RaptorX. Protein-ligand molecular docking studies between PYL and ABA were performed using AutoDock 1.5.6, followed by 100ns molecular dynamic simulation studies using Desmond to determine the acceptable conformational changes after docking via root mean square deviation RMSD plot analysis. Protein-protein docking was then carried out in three sets: PYL-PP2Cs, PYL-ABA-PP2C and PYL(mut)-ABA-PP2C to scrutinize changes in structural conformations and binding energies between complexes. The amino acids of interest were mapped at their respective genomic coordinates using SNP-seek database to ascertain if there were any naturally occurring single nucleotide polymorphisms (SNPs) responsible for triggering rice PYLs mutations. Results: Initial protein-ligand docking studies revealed good binding between the complexes, wherein PYL6-ABA complex showed the best energy of -8.15 kcal/mol. The 100ns simulation studies revealed changes in the RMSD values after docking, indicating acceptable conformational changes. Furthermore, mutagenesis study performed at specific PYL-ABA interacting residues followed by downstream PYL(mut)-ABA-PP2C protein-protein docking results after induction of mutations demonstrated binding energy of -8.17 kcal/mol for PP2C79-PYL11-ABA complex. No naturally occurring SNPs that were responsible for triggering rice PYL mutations were identified when specific amino acid coordinates were mapped at respective genomic coordinates. Conclusion: Thus, the present study provides valuable insights on the interactions of ABA receptors in rice and induced mutations in PYL11 that can enhance the downstream interaction with PP2C.