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

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

Monosialogangliosides and Their Interaction with Cholera Toxin—Investigation by Molecular Modeling and Molecular Mechanics

Abstract

Molecular mechanics and molecular dynamics studies are performed to investigate the conformational preference of cell surface monosialogangliosides (GM3, GM2 and GM1) in aqueous environment. Water mediated hydrogen bonding network plays a significant role in the structural stabilization of GM3, GM2 and GM1. The spatial flexibility of NeuNAc of gangliosides at the binding site of cholera toxin reveals a limited allowed eulerian space of 2.4% with a much less allowed eulerian space (1.4%) for external galactose of GM1. The molecular mechanics of monosialoganglioside-cholera toxin complex reveals that cholera toxin can accommodate the monosialogangliosides in three different modes. Direct and water mediated hydrogen bonding interactions stabilize these binding modes and play an essential role in defining the order of specificity for different monosialogangliosides towards cholera toxin. This study identifies the NeuNAc binding site as a site for design of inhibitors that can restrict the pathogenic activity of cholera toxin.     

A comparative modeling and molecular docking study on Mycobacterium tuberculosis targets involved in peptidoglycan biosynthesis

An alarming rise of multidrug-resistant Mycobacterium tuberculosis strains and the continuous high global morbidity of tuberculosis have reinvigorated the need to identify novel targets to combat the disease. The enzymes that catalyze the biosynthesis of peptidoglycan in M. tuberculosis are essential and noteworthy therapeutic targets. In this study, the biochemical function and homology modeling of MurI, MurG, MraY, DapE, DapA, Alr, and Ddl enzymes of the CDC1551 M. tuberculosis strain involved in the biosynthesis of peptidoglycan cell wall are reported. Generation of the 3D structures was achieved with Modeller 9.13. To assess the structural quality of the obtained homology modeled targets, the models were validated using PROCHECK, PDBsum, QMEAN, and ERRAT scores. Molecular dynamics simulations were performed to calculate root mean square deviation (RMSD) and radius of gyration (Rg) of MurI and MurG target proteins and their corresponding templates. For further model validation, RMSD and Rg for selected targets/templates were investigated to compare the close proximity of their dynamic behavior in terms of protein stability and average distances. To identify the potential binding mode required for molecular docking, binding site information of all modeled targets was obtained using two prediction algorithms. A docking study was performed for MurI to determine the potential mode of interaction between the inhibitor and the active site residues. This study presents the first accounts of the 3D structural information for the selected M. tuberculosis targets involved in peptidoglycan biosynthesis.     

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.     

Conformational Similarity Indices Between Different Residues in Proteins and α-Helix Propensities

Abstract

Various amino acid similarity matrices have been derived using data on physicochemical properties and molecular evolution. Conformational similarity indices, CS_XX′, between different residues are computed here using the distribution of the main-chain and side-chain torsion angles and the values have been used to cluster amino acids in proteins. A subset of these parameters, CS_AX′ indicates the extent of similarity in the main-chain and side-chain conformations (φ ψ and χ₁) of different residues (X) with Ala (A) and is found to have strong correlation with α-helix propensities. However, no subset of CS_XX′ provides any linear relationship with β-sheet propensities, suggesting that the conformational feature favouring the location of a residue in an a-helix is different from the one favouring the β-sheet. Conformationally similar residues (close CS_AX values) have similar steric framework of the side-chain (linear/branched, aliphatic/aromatic), irrespective of the polarity or hydrophobicity. Cooperative nucleation of helix may be facile for a contiguous stretch of residues with high overall CS_AX values.     

In silico prediction of drug resistance due to S247R mutation of Influenza H1N1 neuraminidase protein

We present here in silico studies on antiviral drug resistance due to a novel mutation of influenza A/H1N1 neuraminidase (NA) protein. Influenza A/H1N1 virus was responsible for a recent pandemic and is currently circulating among the seasonal influenza strains. M2 and NA are the two major viral proteins related to pathogenesis in humans and have been targeted for drug designing. Among them, NA is preferred because the ligand-binding site of NA is highly conserved between different strains of influenza virus. Different mutations of the NA active site residues leading to drug resistance or susceptibility of the virus were studied earlier. We report here a novel mutation (S247R) in the NA protein that was sequenced earlier from the nasopharyngeal swab from Sri Lanka and Thailand in the year 2009 and 2011, respectively. Another mutation (S247N) was already known to confer resistance to oseltamivir. We did a comparative study of these two mutations vis-a-vis the drug-sensitive wild type NA to understand the mechanism of drug resistance of S247N and to predict the probability of the novel S247R mutation to become resistant to the currently available drugs, oseltamivir and zanamivir. We performed molecular docking- and molecular dynamics-based analysis of both the mutant proteins and showed that mutation of S247R affects drug binding to the protein by positional displacement due to altered active site cavity architecture, which in turn reduces the affinity of the drug molecules to the NA active site. Our analysis shows that S247R may have high probability of being resistant.     

Homology modeling of human histone deacetylase 10 and design of potential selective inhibitors

Abstract

Histone deacetylases (HDACs) are implicated in the pathology of various cancers, and their pharmacological blockade has proven to be promising in reversing the malignant phenotypes. However, lack of crystal structures of some of the human HDAC isoforms (e.g., HDAC10) hinders the design of the isoform-selective inhibitor. Here, the recently solved X-ray crystal structure of Danio rerio (zebrafish) HDAC10 (Protein Data Bank (PDB) ID; 5TD7, released on 24 May 2017) was retrieved from the PDB and used as a template structure to model the three-dimensional structure of human HDAC10. The overall quality of the best model (M0017) was assessed by computing its z-score—a measure of the deviation of the total energy of the structure with respect to an energy distribution derived from random conformations and by docking of known HDAC10 inhibitors to its catalytic cavity. Furthermore, to identify potential HDAC10-selective inhibitor ligand-based virtual screening was carried out against the ZINC database. The free modeled structure of HDAC10 and its complexes with quisinostat and the highest-ranked compound ZINC19749069 were submitted to molecular dynamics simulation. The comparative analysis of root-mean-squared deviation, root-mean-squared fluctuation, radius of gyration (Rg), and potential energy of these systems showed that HDAC10-ZINC19749069 complex remained the most stable over time. Thus, M0017 could be potentially used for structure-based inhibitor against HDAC10, and ZINC19749069 may provide a scaffold for further optimization.

Communicated by Ramaswamy H. Sarma     

Homology modeling and docking studies of FabH (β-ketoacyl-ACP synthase III) enzyme involved in type II fatty acid biosynthesis of Chlorella variabilis: a potential algal feedstock for biofuel production

The concept of using microalgae as an alternative renewable source of biofuel has gained much importance in recent years. However, its commercial feasibility is still an area of concern for researchers. Unraveling the fatty acid metabolic pathway and understanding structural features of various key enzymes regulating the process will provide valuable insights to target microalgae for augmented oil content. FabH (β-ketoacyl-acyl carrier protein synthase; KAS III) is a condensing enzyme catalyzing the initial elongation step of type II fatty acid biosynthetic process and acyl carrier protein (ACP) facilitates the shuttling of the fatty acyl intermediates to the active site of the respective enzymes in the pathway. In the present study, a reliable three-dimensional structure of FabH from Chlorella variabilis, an oleaginous green microalga was modeled and subsequently the key residues involved in substrate binding were determined by employing protein–protein docking and molecular dynamics (MD) simulation protocols. The FabH-ACP complex having the lowest docking energy score showed the binding of ACP to the electropositive FabH surface with strong hydrogen bond interactions. The MD simulation results indicated that the substrate-complexed FabH adopted a more stable conformation than the free enzyme. Further, the FabH structure retained its stability throughout the simulation although noticeable displacements were observed in the loop regions. Molecular simulation studies suggested the importance of crucial hydrogen bonding of the conserved Arg⁹¹ of FabH with Glu⁵³ and Asp⁵⁶ of ACP for exhibiting high affinity between the enzyme and substrate. The molecular modeling results are consistent with available experimental results on the flexibility of FabH and the present study provides first in silico insights into the structural and dynamical aspect of catalytic mechanism of FabH, which could be used for further site-specific mutagenic experiments to develop engineered high oil-yielding microalgal strains for biofuel production.     

Molecular modeling studies on nucleoside hydrolase from the biological warfare agent Brucella suis

Brucella suis is a dangerous biological warfare agent already used for military purposes. This bacteria cause brucellosis, a zoonosis highly infective and difficult to fight. An important selective target for chemotherapy against this disease is nucleoside hydrolase (NH), an enzyme still not found in mammals. We present here the first three-dimensional structure of B. suis NH (BsNH) and propose this enzyme as a molecular target to the drug design in the fight against brucellosis. In addition, we performed molecular docking studies, aiming to analyze the three-dimensional positioning of nine known inhibitors of Chritidia fasciculata NH (CfNH) in the active sites of BsNH and CfNH. We also analyzed the main interactions of some of these compounds inside the active site of BsNH and the relevant factors to biological activity. These results, together with further molecular dynamics (MD) simulations, pointed out to the most promising compound as lead for the design of potential inhibitors of BsNH. Most of the docking and MD results corroborated to each other and the docking results also suggested a good correlation with experimental data.     

A novel parameterization scheme for energy equations and its use to calculate the structure of protein molecules

A novel scheme for the parameterization of a type of “potential energy” function for protein molecules is introduced. The function is parameterized based on the known conformations of previously determined protein structures and their sequence similarity to a molecule whose conformation is to be calculated. Once parameterized, minima of the potential energy function can be located using a version of simulated annealing which has been previously shown to locate global and near-global minima with the given functional form. As a test problem, the potential was parameterized based on the known structures of the rubredoxins from Desulfovibrio vulgaris, Desulfovibrio desulfuricans, and Clostridium pasteurianum, which vary from 45 to 54 amino acids in length, and the sequence alignments of these molecules with the rubredoxin sequence from Desulfovibrio gigas. Since the Desulfovibrio gigas rubredeoxin conformation has also been determined, it is possible to check the accuracy of the results. Ten simulated-annealing runs from random starting conformations were performed. Seven of the 10 resultant conformations have an all-C_α rms deviation from the crystallographically determined conformation of less than 1.7 Å. For five of the structures, the rms deviation is less than 0.8 Å. Four of the structures have conformations which are virtually identical to each other except for the position of the carboxy-terminal residue. This is also the conformation which is achieved if the determined crystal structure is minimized with the same potential. The all-C_α rms difference between the crystal and minimized crystal structures is 0.6 Å. It is further observed that the “energies” of the structures according to the potential function exhibit a strong correlation with rms deviation from the native structure. The conformations of the individual model structures and the computational aspects of the modeling procedure are discussed. © 1993 Wiley-Liss, Inc.     

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.     

Mechanism of the plant cytochrome P450 for herbicide resistance: a modelling study

Plant cytochrome P450 is a key enzyme responsible for the herbicide resistance but the molecular basis of the mechanism is unclear. To understand this, four typical plant P450s and a widely resistant herbicide chlortoluron were analysed by carrying out homology modelling, molecular docking, molecular dynamics simulations and binding free energy analysis. Our results demonstrate that: (i) the putative hydrophobic residues located in the F-helix and polar residues in I-helix are critical in the herbicide resistance; (ii) the binding mode analysis and binding free energy calculation indicate that the distance between catalytic site of chlortoluron and heme of P450, as well as the binding affinity are key elements affecting the resistance for plants. In conclusion, this work provides a new insight into the interactions of plant P450s with herbicide from a molecular level, offering valuable information for the future design of novel effective herbicides which also escape from the P450 metabolism.     

A profound computational study to prioritize the natural compound inhibitors against the P. falciparum orotidine-5-monophosphate decarboxylase enzyme

Abstract

In the current contribution, a multicomplex-based pharmacophore modeling approach was employed on the structural proteome of Plasmodium falciparum orotidine-5-monophosphate decarboxylase enzyme (PfOMPDC). Among the constructed pharmacophore models, the representative hypotheses were selected as the primary filter to screen the molecules with the complementary features responsible for showing inhibition. Thereafter, auxiliary evaluations were performed on the screened candidates via drug-likeness and molecular docking studies. Subsequently, the stability of the docked protein-ligand complexes was scrutinized by employing molecular dynamics simulations and molecular mechanics-Poisson Boltzmann surface area based free binding energy calculations. The stability the docked candidates was compared with the highly active crystallized inhibitor (3S9Y-FNU) to seek more potential candidates. All the docked molecules displayed stable dynamic behavior and high binding free energy in comparison to 3S9Y-FNU. The employed workflow resulted in the retrieval of five drug-like candidates with diverse scaffolds that may show inhibitory activity against PfOMPDC and could be further used as the novel scaffold to develop novel antimalarials.

Communicated by Ramaswamy H. Sarma