Molecular docking and dynamics simulation analyses unraveling the differential enzymatic catalysis by plant and fungal laccases with respect to lignin biosynthesis and degradation

Laccase, widely distributed in bacteria, fungi, and plants, catalyzes the oxidation of wide range of compounds. With regards to one of the important physiological functions, plant laccases are considered to catalyze lignin biosynthesis while fungal laccases are considered for lignin degradation. The present study was undertaken to explain this dual function of laccases using in-silico molecular docking and dynamics simulation approaches. Modeling and superimposition analyses of one each representative of plant and fungal laccases, namely, Populus trichocarpa and Trametes versicolor, respectively, revealed low level of similarity in the folding of two laccases at 3D levels. Docking analyses revealed significantly higher binding efficiency for lignin model compounds, in proportion to their size, for fungal laccase as compared to that of plant laccase. Residues interacting with the model compounds at the respective enzyme active sites were found to be in conformity with their role in lignin biosynthesis and degradation. Molecular dynamics simulation analyses for the stability of docked complexes of plant and fungal laccases with lignin model compounds revealed that tetrameric lignin model compound remains attached to the active site of fungal laccase throughout the simulation period, while it protrudes outwards from the active site of plant laccase. Stability of these complexes was further analyzed on the basis of binding energy which revealed significantly higher stability of fungal laccase with tetrameric compound than that of plant. The overall data suggested a situation favorable for the degradation of lignin polymer by fungal laccase while its synthesis by plant laccase.     

Pharmacophore-directed Homology Modeling and Molecular Dynamics Simulation of G Protein-coupled Receptor： Study of Possible Binding Modes of 5-HT2c Receptor Agonists

A new pharmacophore-based modeling procedure, including homology modeling, pharmacophore study, flexible molecular docking, and long-time molecular dynamics (MD) simulations, was employed to construct the structure of the human 5-HT_(2C) receptor and determine the characteristics of binding modes of 5-HT_(2C) receptor agonists. An agonist-receptor complex has been constructed based on homology modeling and a pharmacophore hypothesis model based on some high active compounds. Then MD simulations of the ligand-receptor complex in an explicit membrane environment were carried out. The conformation of the 5- HT_(2C) receptor during MD simulation was explored, and the stable binding modes of the studied agonist were determined. Flexible molecular docking of several structurally diverse agonists of the human 5-HT_(2C) receptor was carried out, and the general binding modes of these agonists were investigated. According to the models presented in this work and the results of Flexi-Dock, the involvement of the amino acid residues Asp134, Ser138, Ash210, Asn331, Tyr358, Ile131, Ser132, Val135, Thr139, Ile189, Val202, Val208, Leu209, Phe214, Val215, Gly218, Ser219, Phe223, Trp324, Phe327, and Phe328 in agonist recognition was studied. The obtained binding modes of the human 5-HT_(2C) receptor agonists have good agreement with the site-directed mutagenesis data and other studies.     

抗癌晶体蛋白Parasporin-2空间结构的研究

通过同源建模得到了抗癌晶体蛋白Parasporin-2(Cry46Aa1)活性区的初始三维结构,利用分子动力学方法对初始三维结构进行优化。为了评估模型的好坏,利用Ramachandran plot和结构匹配等方法对模型进行评价。结果显示,所得的Parasporin-2空间模型良好。比较了Cry46Aa1与Cry46Ab1这两种高度同源的抗癌晶体蛋白的空间结构,找出了其空间结构上的不同之处。通过大分子对接程序Hex4.5,模拟了Parasporin-2与受体GPI-瞄固蛋白(CD59)的相互作用。结果显示,Parasporin-2中参与对接的活性位点位于两个α-螺旋下方的凹槽内,而CD59中的四个loop倾向于插入该凹槽内。研究结果为Parasporin-2的抗癌机制研究及其理性改造奠定了基础。     

Homology modeling of the Abl-SH3 domain

A tertiary structure model of the Abl-SH3 domain is predicted by using homology modeling techniques coupled to molecular dynamics simulations. Two template proteins were used, Fyn-SH3 and Spc-SH3. The refined model was extensively checked for errors using criteria based on stereochemistry, packing, solvation free-energy, accessible surface areas, and contact analyses. The different checking methods do not totally agree, as each one evaluates a different characteristic of protein structures. Several zones of the protein are more susceptible to incorporating errors. These include residues 13, 15, 35, 39, 45, 46, 50, and 60. An interesting finding is that the measurement of the Cα chirality correlated well with the rest of the criteria, suggesting that this parameter might be a good indicator of correct local conformation. Deviations of more than 4 degrees may be indicative of poor local structure. © 1994 Wiley-Liss, Inc.     

Isozymes inhibited by active site blocking: versatility of calcium indifferent hesperidin binding to phospholipase A2 and its significance

sPLA₂ is released under inflammatory conditions from neutrophils, basophils and T-cells. They cleave the cellular phospholipids leading to the release of arachidonic acid and there by provide intermediates for biosynthesis of inflammatory mediators. The focus of this study is on the interaction of hesperidin, a natural flavonoid with Group IB, IIA, and V and X isozymes of sPLA₂. Affinity of hesperidin towards PLA₂ isozymes was analyzed through enzymatic studies and molecular modeling. The experiments showed that hesperidin competitively inhibited PLA₂ with IC₅₀ of 5.1?µM. Molecular modeling studies revealed the association of hesperidin with the docking scores ?6.90, ?9.53, ?5.63 and ?8.29?kcal for isozymes Group IB, IIA, V and X of PLA₂ respectively. Their binding energy values were calculated as ?20.25, ?21.63, ?21.66 and ?33.43?kcal for the Group IB, IIA, V and X respectively. Structural model for Group V was made by homology modeling since no structural coordinates were available. Molecular dynamics studies were carried out to evaluate the structural stability of protein ligand complex. The analyses showed that hesperidin blocked the entry of the substrate to the active site of PLA₂ and it was indifferent to the differences of the isozymes. Hence, hesperidin might serve as lead for designing highly specific anti-inflammatory drugs directed to the PLA₂ isozyme specific to various diseases, with IC₅₀ value of therapeutic significance.     

Investigation of ligand selectivity in CYP3A7 by molecular dynamics simulations

Cytochrome P450 (CYP) 3A7 plays a crucial role in the biotransformation of the metabolized endogenous and exogenous steroids. To compare the metabolic capabilities of CYP3A7–ligands complexes, three endogenous ligands were selected, namely dehydroepiandrosterone (DHEA), estrone, and estradiol. In this study, a three-dimensional model of CYP3A7 was constructed by homology modeling using the crystal structure of CYP3A4 as the template and refined by molecular dynamics simulation (MD). The docking method was adopted, combined with MD simulation and the molecular mechanics generalized born surface area method, to probe the ligand selectivity of CYP3A7. These results demonstrate that DHEA has the highest binding affinity, and the results of the binding free energy were in accordance with the experimental conclusion that estrone is better than estradiol. Moreover, several key residues responsible for substrate specificity were identified on the enzyme. Arg372 may be the most important residue due to the low interaction energies and the existence of hydrogen bond with DHEA throughout simulation. In addition, a cluster of Phe residues provides a hydrophobic environment to stabilize ligands. This study provides insights into the structural features of CYP3A7, which could contribute to further understanding of related protein structures and dynamics.     

Cell‐free protein synthesis and purification of human dopamine D2 receptor long isoform

The human dopamine D2 receptor long isoform (D2L) has significant implications in neurological and neuropsychiatric disorders such as Parkinson's disease and schizophrenia. Detailed structural knowledge of this receptor is limited owing to its highly hydrophobic nature, which leads to protein aggregation and host toxicity when expressed in cellular systems. The newly emerging field of cell‐free protein expression presents numerous advantages to overcome these challenges. This system utilizes protein synthesis machinery and exogenous DNA to synthesize functional proteins outside of intact cells. This study utilizes two different cell‐free systems for the synthesis of human dopamine D2L receptor. These include the Escherichia coli lysate‐based system and the wheat‐germ lysate‐based system. The bacterial cell‐free method used pET 100/D‐TOPO vector to synthesize hexa‐histidine‐tagged D2L receptor using a dialysis bag system; the resulting protein was purified using nickel‐nitrilotriacetic acid affinity resin. The wheat germ system used pEU–glutathione‐S‐transferase (GST) vector to synthesize GST‐tagged D2L receptor using a bilayer translation method; the resulting protein was purified using a GST affinity resin. The presence and binding capacity of the synthesized D2L receptor was confirmed by immunoblotting and radioligand competition assays, respectively. Additionally, in‐gel protein sequencing via Nano LC‐MS/MS was used to confirm protein synthesis via the wheat germ system. The results showed both systems to synthesize microgram quantities of the receptor. Improved expression of this highly challenging protein can improve research and understanding of the human dopamine D2L receptor. © 2013 American Institute of Chemical Engineers Biotechnol. Prog., 29:601–608, 2013     

A coupling of homology modeling with multiple molecular dynamics simulation for identifying representative conformation of GPCR structures: a case study on human bombesin receptor subtype-3

In this study, a computational pipeline was therefore devised to overcome homology modeling (HM) bottlenecks. The coupling of HM with molecular dynamics (MD) simulation is useful in that it tackles the sampling deficiency of dynamics simulations by providing good-quality initial guesses for the native structure. Indeed, HM also relaxes the severe requirement of force fields to explore the huge conformational space of protein structures. In this study, the interaction between the human bombesin receptor subtype-3 and MK-5046 was investigated integrating HM, molecular docking, and MD simulations. To improve conformational sampling in typical MD simulations of GPCRs, as in other biomolecules, multiple trajectories with different initial conditions can be employed rather than a single long trajectory. Multiple MD simulations of human bombesin receptor subtype-3 with different initial atomic velocities are applied to sample conformations in the vicinity of the structure generated by HM. The backbone atom conformational space distribution of replicates is analyzed employing principal components analysis. As a result, the averages of structural and dynamic properties over the twenty-one trajectories differ significantly from those obtained from individual trajectories.     

How good are comparative models in the understanding of protein dynamics?

The 3D structure of a protein is essential to understand protein dynamics. If experimentally determined structure is unavailable, comparative models could be used to infer dynamics. However, the effectiveness of comparative models, compared to experimental structures, in inferring dynamics is not clear. To address this, we compared dynamics features of ~800 comparative models with their crystal structures using normal mode analysis. Average similarity in magnitude, direction, and correlation of residue motions is >0.8 (where value 1 is identical) indicating that the dynamics of models and crystal structures are highly similar. Accuracy of 3D structure and dynamics is significantly higher for models built on multiple and/or high sequence identity templates (>40%). Three-dimensional (3D) structure and residue fluctuations of models are closer to that of crystal structures than to templates (TM score 0.9 vs 0.7 and square inner product 0.92 vs 0.88). Furthermore, long-range molecular dynamics simulations on comparative models of RNase 1 and Angiogenin showed significant differences in the conformational sampling of conserved active-site residues that characterize differences in their activity levels. Similar analyses on two EGFR kinase variant models highlight the effect of mutations on the functional state-specific αC helix motions and these results corroborate with the previous experimental observations. Thus, our study adds confidence to the use of comparative models in understanding protein dynamics.