Structural dynamics of wild type and mutated forms of human L1 endonuclease and insights into its sequence specific nucleic acid binding mechanism: A molecular dynamics study

Insights into the RNA binding mechanism of human L1-ORF1p: a molecular dynamics study

The recognition and binding of nucleic acids by ORF1p, an L1 retrotransposon protein, have not yet been clearly understood due to the lack of structural knowledge. The present study attempts to identify the probable single-stranded RNA binding pathway of trimeric ORF1p using computational methods like ligand mapping methodology combined with molecular dynamics simulations. Using the ligand mapping methodology, the possible RNA interacting sites on the surface of the trimeric ORF1p were identified. The crystal structure of the ORF1p timer and an RNA molecule of 29 nucleotide bases in length were used to generate the structure of the ORF1p complex based on information on predicted binding sites as well as the functional states of the CTD. The various complexes of ORF1p-RNA were generated using polyU, polyA and L1RNA sequences and were simulated for a period of 75 ns. The observed stable interaction pattern was used to propose the possible binding pathway. Based on the binding free energy for complex formation, both polyU and L1RNA complexes were identified as stable complexes, while the complex formed with polyA was the least stable one. Furthermore, the importance of the residues in the CC domain (Lys137 and Arg141), the RRM loop (Arg206, Arg210 and Arg211) and the CTD (Arg 261 and Arg262) of all three chains in stabilizing the wrapped RNA has been highlighted in this study. The presence of several electrostatic interactions including H-bond interactions increases the affinity towards RNA and hence plays a vital role in retaining the wrapped position of RNA around ORF1p. Altogether, this study presents one of the possible RNA binding pathways of ORF1p and clearly highlights the functional state of ORF1p visited during RNA binding.

Molecular Modeling and Structural Stability of Wild-Type and Mutant CYP51 from Leishmania major: In Vitro and In Silico Analysis of a Laboratory Strain

Cutaneous leishmaniasis is a neglected tropical disease and a major public health in the most countries. Leishmania major is the most common cause of cutaneous leishmaniasis. In the Leishmania parasites, sterol 14α-demethylase (CYP51), which is involved in the biosynthesis of sterols, has been identified as an attractive target for development of new therapeutic agents. In this study, the sequence and structure of CYP51 in a laboratory strain (MRHO/IR/75/ER) of L. major were determined and compared to the wild-type strain. The results showed 19 mutations including seven non-synonymous and 12 synonymous ones in the CYP51 sequence of strain MRHO/IR/75/ER. Importantly, an arginine to lysine substitution at position of 474 resulted in destabilization of CYP51 (ΔΔG = 1.17 kcal/mol) in the laboratory strain; however, when the overall effects of all substitutions were evaluated by 100 ns molecular dynamics simulation, the final structure did not show any significant changes (p-value < 0.05) in stability parameter of the strain MRHO/IR/75/ER compared to the wild-type protein. The energy level for the CYP51 of wild-type and MRHO/IR/75/ER strain were −40,027.1 and −39,706.48 Kcal/mol respectively. The overall Root-mean-square deviation (RMSD) deviation between two proteins was less than 1 Å throughout the simulation and Root-mean-square fluctuation (RMSF) plot also showed no substantial differences between amino acids fluctuation of the both protein. The results also showed that, these mutations were located on the protein periphery that neither interferes with protein folding nor with substrate/inhibitor binding. Therefore, L. major strain MRHO/IR/75/ER is suggested as a suitable laboratory model for studying biological role of CYP51 and inhibitory effects of sterol 14α-demethylase inhibitors.