Homology model of human retinoic acid metabolising enzyme cytochrome P450 26A1 (CYP26A1): active site architecture and ligand binding

Homology models of cytochrome P450 RA1 (CYP26A1) were constructed using three human P450 structures, CYP2C8, CYP2C9 and CYP3A4 as templates for the model building. Using MOE software the lowest energy CYP26A1 model was then assessed for stereochemical quality and side chain environment. Further active site optimisation of the CYP26A1 model built using the CYP3A4 template was performed by molecular dynamics to generate a final CYP26A1 model. The natural substrate, all-trans-retinoic acid (atRA), and inhibitor R 15866, were docked into the model allowing further validation of the active site architecture. Using the docking studies structurally and functionally important residues were identified with subsequent characterisation of secondary structure. Multiple hydrophobic interactions, including the side chains of TRP112, PHE299, PHE222, PHE84, PHE374 and PRO371, are important for binding of atRA and R115866. Additional hydrogen bonding interactions were noted as follows: atRA-- C==O of the atRA carboxylate group and ARG86; R115866--benzothiazole nitrogen and the backbone NH of SER115.     

Targeting methionyl tRNA synthetase: design,synthesis and antibacterial activity against Clostridium difficile of novel 3-biaryl-N-benzylpropan-1-amine derivatives

The synthesis of a series of benzimidazole-N-benzylpropan-1-amines and adenine-N-benzylpropan-1-amines is described. Subsequent evaluation against two strains of the anaerobic bacterium Clostridium difficile was performed with three amine derivatives displaying MIC values of 16?μg/mL. Molecular docking studies of the described amines determined that the amines interact within two active site pockets of C. difficile methionyl tRNA synthetase with methoxy substituents in the benzyl ring and an adenine biaryl moiety resulting in optimal binding interactions.     

Homology modeling,molecular docking and spectra assay studies of sterol 14α-demethylase from <Emphasis Type="Italic">Penicillium digitatum</Emphasis>

Sterol 14α-demethylase from Penicillium digitatum (PdCYP51) is a prime target of antifungal drugs for citrus disease in plants. To design novel antifungal compounds, a homology model of PdCYP51 was constructed using the recently reported crystal structure of human CYP51 as the template. Molecular docking was performed to investigate the interaction of four commercial fungicides with the modeled enzyme. The side chain of these compounds interplayed with PdCYP51 mainly through hydrophobic and van der Waals interactions. Biochemical spectra analysis of inhibitors combined with PdCYP51 are also compatible with the docking results. This is the first molecular modeling for PdCYP51 based on the eukaryotic crystal structure of CYP51. The structural information and binding site mapping of PdCYP51 for different inhibitors obtained from this study could aid in screening and designing new antifungal compounds targeting this enzyme.     

Structural and functional conservation profiles of novel cathepsin L-like proteins identified in the Drosophila melanogaster genome

Cathepsin L is a cysteine protease which degrades connective tissue proteins including collagen, elastin, and fibronectin. In this study, five well-characterized cathepsin L proteins from different arthropods were used as query sequences for the Drosophila genome database. The search yielded 10 cathepsin L-like sequences, of which eight putatively represent novel cathepsin L-like proteins. To understand the phylogenetic relationship among these cathepsin L-like proteins, a phylogenetic tree was constructed based on their sequences. In addition, models of the tertiary structures of cathepsin L were constructed using homology modeling methods and subjected to molecular dynamics simulations to obtain reasonable structure to understand its dynamical behavior. Our findings demonstrate that all of the potential Drosophila cathepsin L-like proteins contain at least one cathepsin propeptide inhibitor domain. Multiple sequence alignment and homology models clearly highlight the conservation of active site residues, disulfide bonds, and amino acid residues critical for inhibitor binding. Furthermore, comparative modeling indicates that the sequence/structure/function profiles and active site architectures are conserved.     

Fragment screening and structural analyses highlight the ATP-assisted ligand binding for inhibitor discovery against type 1 methionyl-tRNA synthetase

Methionyl-tRNA synthetase (MetRS) charges tRNA^Met with l-methionine (L-Met) to decode the ATG codon for protein translation, making it indispensable for all cellular lives. Many gram-positive bacteria use a type 1 MetRS (MetRS1), which is considered a promising antimicrobial drug target due to its low sequence identity with human cytosolic MetRS (HcMetRS, which belongs to MetRS2). Here, we report crystal structures of a representative MetRS1 from Staphylococcus aureus (SaMetRS) in its apo and substrate-binding forms. The connecting peptide (CP) domain of SaMetRS differs from HcMetRS in structural organization and dynamic movement. We screened 1049 chemical fragments against SaMetRS preincubated with or without substrate ATP, and ten hits were identified. Four cocrystal structures revealed that the fragments bound to either the L-Met binding site or an auxiliary pocket near the tRNA CCA end binding site of SaMetRS. Interestingly, fragment binding was enhanced by ATP in most cases, suggesting a potential ATP-assisted ligand binding mechanism in MetRS1. Moreover, co-binding with ATP was also observed in our cocrystal structure of SaMetRS with a class of newly reported inhibitors that simultaneously occupied the auxiliary pocket, tRNA site and L-Met site. Our findings will inspire the development of new MetRS1 inhibitors for fighting microbial infections.     

Caseoperoxidase,mixed β-casein–SDS–hemin–imidazole complex: a nano artificial enzyme

A novel peroxidase-like artificial enzyme, named “caseoperoxidase”, was biomimetically designed using a nano artificial amino acid apo-protein hydrophobic pocket. This four-component nano artificial enzyme containing heme–imidazole–β-casein–SDS exhibited high activity growth and k_cat performance toward the native horseradish peroxidase demonstrated by the steady state kinetics using UV–vis spectrophotometry. The hydrophobicity and secondary structure of the caseoperoxidase were studied by ANS fluorescence and circular dichroism spectroscopy. Camel β-casein (Cβ-casein) was selected as an appropriate apo-protein for the heme active site because of its innate flexibility and exalted hydrophobicity. This selection was confirmed by homology modeling method. Heme docking into the newly obtained Cβ-casein structure indicated one heme was mainly incorporated with Cβ-casein. The presence of a main electrostatic site for the active site in the Cβ-casein was also confirmed by experimental methods through Wyman binding potential and isothermal titration calorimetry. The existence of Cβ-casein protein in this biocatalyst lowered the suicide inactivation and provided a suitable protective role for the heme active-site. Additional experiments confirmed the retention of caseoperoxidase structure and function as an artificial enzyme.     

Involvement of cysteine 306 and alanine 63 in the thermostability and oligomeric organization of glucose isomerase from <Emphasis Type="Italic">Streptomyces</Emphasis> sp. SK

The implication of the original alanine 63 (Ala63) and the unique cysteine 306 (Cys306) residues in the thermostability of the Streptomyces sp. SK glucose isomerase (SKGI) were investigated by site-directed mutagenesis and homology modelling. The Cys306 to Ala mutation within SKGI dramatically affected its thermal stability by decreasing the half-life from 80 to 15 min at 90°C while the Ala63 to Ser replacement shifted this half-life to 65 min. The electrophoretic analysis proves that the residue Cys306 participates in oligomerization of the SKGI. Its stabilizing role is materialized by hydrogen bonds established with arginines at positions 284 and 259, as deduced from the constructed three-dimensional model. We have also shown that the presence of an Ala63 instead of Ser63 seems to be more suitable for enzyme thermostability by maintaining hydrophobic pocket that contributes to the protection of the enzyme active site.     

Computer-aided molecular design of 1H-imidazole-2,4-diamine derivatives as potential inhibitors of <Emphasis Type="Italic">Plasmodium falciparum</Emphasis> DHFR enzyme

Design and discovery of new potential inhibitors of Plasmodium falciparum dihydrofolate reductase (PfDHFR), equally active against both the wild-type and mutant strains, is urgently needed. In this study, a computer-aided molecular design approach that involved ab initio molecular orbital and density functional theory calculations, along with molecular electrostatic potential analysis, and molecular docking studies was employed to design 15 1H-imidazole-2,4-diamine derivatives as potential inhibitors of PfDHFR enzyme. Visual inspection of the binding modes of the compounds demonstrated that they all interact, via H-bond interactions, with key amino acid residues (Asp54, Ileu/Leu164, Asn/Ser108 and Ile14) similar to those of WR99210 (3) in the active site of the enzymes used in the study. These interactions are known to be essential for enzyme inhibition. These compounds showed better or comparable binding affinities to that of the bound ligand (WR99210). In silico toxicity predictions, carried out using TOPKAT software, also indicated that the compounds are non-toxic.     

A Three Dimensional Model of Human Thiopurine Methyltransferase; Ligand Interactions and Structural Consequences of Naturally Occurring Mutations

The three-dimensional model of human thiopurine methyltransferase (hTPMT) was constructed by molecular modeling. A multiple alignment of AdoMet dependent methyltransferases based on a structural superposition of the AdoMet binding domain of Hhai, TaqI and rCOMT was used in the modeling procedure. The reliability of the model was examined by comparing its conformation and packing properties with those of Hhai, TaqI and rCOMT and structures in the PDB-database. The examined criteria indicated a reliable model structure. The model gave insight into the structural effects of naturally occurring mutations of the hTPMT allele, and was used to characterize the ligand interactions of the protein. The residues Gln42 and Glu91 were predicted to participate in AdoMet binding through H-bond interactions whereas Phe146 participates through Van der Waal interaction. The cationic methyl-sulphonium group of AdoMet was located close to the aromatic residue Phe40. The model also indicated that substrates interact with hTPMT situated in a pocket consisting of the hydrophobic residues Phe40, Met148, Val184, Val220 and the charged residues Lys145, Glu218, Lys219. These residues were also included in a predictive explanation for the inhibitor/substrate preference of the enzyme. The most frequent of naturally occurring mutations was predicted to cause alterations on the surface of the protein with minor/none structural consequences. The mutation Ala80-Pro seemed directly to cause an inactive enzyme by disrupting the structure of the binding site of AdoMet.Electronic Supplementary Material available.     

3D-structure of human estrogenic 17β-HSD1: binding with various steroids

Human estrogenic dehydrogenase (17β-HSD1) catalyses the last step in the biosynthesis of the active estrogens that stimulate the proliferation of breast cancer cells. While the primary substrate for the enzyme is estrone, the enzyme has some activity for the non-estrogenic substrates. To better understand the structure–function relationships of 17β-HSD1 and to provide a better ground for the design of inhibitors, we have determined the crystal structures of 17β-HSD1 in complex with different steroids.The structure of the complex of estradiol with the enzyme determined previously (Azzi et al., Nature Structural Biology 3, 665–668) showed that the narrow active site was highly complementary to the substrate. The substrate specificity is due to a combination of hydrogen bonding and hydrophobic interactions between the steroid and the enzyme binding pocket. We have now determined structures of 17β-HSD1 in complex with dihydrotestosterone and 20α-OH-progesterone. In the case of the C19 androgen, several residues within the enzyme active site make some small adjustments to accommodate the increased bulk of the substrate. In addition, the C19 steroids bind in a slightly different position from estradiol with shifts in positions of up to 1.4 Å. The altered binding position avoids unfavorable steric interactions between Leu 149 and the C19 methyl group (Han et al., unpublished). The known kinetic parameters for these substrates can be rationalized in light of the structures presented. These results give evidence for the structural basis of steroid recognition by 17β-HSD1 and throw light on the design of new inhibitors for this pivotal steroid enzyme.     

β‐Ionone as putative semiochemical suggested by ligand binding on an odorant‐binding protein of Hylamorpha elegans and electroantennographic recordings

Currently, odorant‐binding proteins (OBPs) are considered the first filter for olfactory information for insects and constitute an interesting target for pest control. Thus, an OBP (HeleOBP) from the scarab beetle Hylamorpha elegans (Burmeister) was identified, and ligand‐binding assays based on fluorescence and in silico approaches were performed, followed by a simulated binding assay. Fluorescence binding assays showed slight binding for most of the ligands tested, including host‐plant volatiles. A high binding affinity was obtained for β‐ionone, a scarab beetle‐related compound. However, the binding of its analogue α‐ionone was weaker, although it is still considered good. On the other hand, through a three‐dimensional model of HeleOBP constructed by homology, molecular docking was carried out with 29 related ligands to the beetle. Results expressed as free binding energy and fit quality (FQ) indicated strong interactions of sesquiterpenes and terpenoids (α‐ and β‐ionone) with HeleOBP as well as some aromatic compounds. Residues such as His102, Tyr105 and Tyr113 seemed to participate in the interactions previously mentioned. Both in silico scores supported the experimental affinity for the strongest ligands. Therefore, the activity of α‐ionone, β‐ionone and 2‐phenyl acetaldehyde at antennal level was studied using electroantenography (EAG). Results showed that the three ligands are electrophysiologically active. However, an aliquot of β‐ionone (represented by 3.0 ng) elicited stronger EAG responses in antennae of males than of females. Finally, the role of these ligands as potential semiochemicals for H. elegans is discussed.     

Homology model of 1alpha,25-dihydroxyvitamin D3 24-hydroxylase cytochrome P450 24A1 (CYP24A1): active site architecture and ligand binding

Homology models of cytochrome P450 24A1 (CYP24A1) were constructed using three human P450 structures, CYP2C8, CYP2C9 and CYP3A4 as templates for the model building. Using molecular operating environment (MOE) software the lowest energy CYP24A1 model was then assessed for stereochemical quality and side chain environment. Further active site optimisation of the CYP24A1 model built using the CYP3A4 template was performed by molecular dynamics to generate a final CYP24A1 model. The natural substrate, 1,25-dihydroxyvitamin D(3) (calcitriol) and the CYP24 inhibitor (R)-N-(2-(1H-imidazol-1-yl)-2-phenylethyl)-4'-chlorobiphenyl-4-carboxamide ((R)-VID-400) were docked into the model allowing further validation of the active site architecture. Using the docking studies structurally and functionally important residues were identified with subsequent characterisation of secondary structure.     

Molecular modeling of <Emphasis Type="Italic">Trypanosoma cruzi</Emphasis> glutamate cysteine ligase and investigation of its interactions with glutathione

Trypanosoma cruzi glutamate cysteine ligase (TcGCL) is considered a potential drug target to develop novel antichagasic drugs. We have used a variety of computational methods to investigate the interactions between TcGCL with Glutathione (GSH). The three-dimensional structure of TcGCL was constructed by comparative modeling methods using the Saccharomyces cerevisiae glutamate cysteine ligase as template. Molecular dynamics simulations were used to validate the TcGCL model and to analyze the molecular interactions with GSH. Using RMSD clustering, the most prevalent GSH binding modes were identified paying attention to the residues involved in the molecular interactions. The GSH binding modes were used to propose pharmacophore models that can be exploited in further studies to identify novel antichagasic compounds.     

Molecular Modeling of Cytochrome P450 1A1: Enzyme-Substrate Interactions and Substrate Binding Affinities

Abstract

Human cytochrome P450 1A1, which is present in lungs, plays an important role in the metabolic activation of chemical carcinogens, and in particular, is thought to be linked to lung cancer. The mechanism of carcinogenesis is related to the enzyme's ability to oxidize highly toxic compounds, such as polycyclic aromatic hydrocarbons (PAHs), to their carcinogenic derivatives.

In order to better understand P450 1A1 function, a homology model of this enzyme has been constructed. The model has been based on the structure of P450 2C5, the first mammalian P450 to be crystallized. The coordinates of the model have been calculated using a consensus strategy, and the resulting structure has been evaluated with the ProStat and Profiles-3D programs. P450 1A1 substrates, such as benzo[a]pyrene, ethoxyresorufin and methoxyre- sorufin, were then docked into the active site of the model, and key amino acid residues able to interact with the substrate, have been identified. The analysis of enzyme-substrate interactions indicated that hydrophobic interactions are mainly responsible for binding of these substrates in the active site. Moreover, the non-bond enzyme-substrate interaction energy for ethoxyresorufin was lower than that for methoxyresorufin, which is consistent with higher activity of 1A1 towards the former substrate. Key residue Val-382 may play an important role in these interactions. Additionally, we performed binding free energy calculations for the three substrates. The obtained values were similar to those observed experimentally, which suggests that this approach might be useful for prediction of binding constants.     

Molecular Interaction Study to Explore the <i>Nigella sativa</i> Bioactive Components as an Inhibitor of Peptide Deformylase to Inhibit the <i>Xanthomonas oryzae</i> pv. oryzae via Applying Computational Approach

Bacterial leaf blight (BLB) caused by Xanthomonas oryzae pv. oryzae (Xoo) is one of the most damaging diseases to rice across the world. Various chemicals have been employed so far for the management of bacterial leaf blight. On the other hand, these compounds are damaging to the ecosystem and have an impact on non-target species such as humans and animals. As a result, there is a need to create a new natural inhibitor for BLB management. Deformylase (PDF) enzyme is present in all eubacteria and its necessity in bacterial protein synthesis reveals it as an attractive target for drug development. In this study, the active components of Nigella sativa have been selected based on their previously reported antimicrobial activity and screened on the active site of bacterial PDF by the in silico art of techniques. Among these compounds, dithymoquinone and p-cymene strongly bind with the PDF enzyme with binding energy values of 7.77 kcal/mol and 7.26 kcal/mol, respectively, which is comparatively higher than the control compound (−6.73 kcal/mol). Hence, the “dithymoquinone-PDF” and “p-cymene-PDF” complexes were selected for further study, and their stability was assessed by molecular dynamic (MD) simulation. In MD simulation, both selected compounds exhibited steady-state interaction with PDF for 20 ns. It has been hypothesized that p-cymene and dithymoquinone inhibit peptide deformylase and could be used as antibacterials or pesticides against Xoo against the BLB disease.     

Purification and Identification of a Novel Leucine Aminopeptidase from <Emphasis Type="Italic">Bacillus thuringiensis israelensis</Emphasis>

A novel leucine aminopeptidase was purified from a Bacillus thuringiensis israelensis (Bti) culture. The purification stages included heating the concentrated supernatant to 65°C for 90 min, anion-exchange chromatography by DEAE cellulose, and hydrophobic chromatography by phenyl Sepharose. The specific activity of leucine aminopeptidase after the hydrophobic chromatography increased by 215.5-fold and the yield was 16%. The molecular weight of the active enzyme was 59 kDa. Mass spectrometry analysis of the 59-kDa leucine aminopeptidase revealed that this protein has at least 41% homology with the cytosol leucine aminopeptidase produced by Bacillus cereus. Maximal leucine aminopeptidase activity occurred at 65°C, pH 10 toward leucine as the amino acid terminus. The enzyme was strongly inhibited by bestatin, dithiothreitol, and 1,10-phenanthroline, indicating that the enzyme might be considered as a metallo-aminopeptidase that has disulfide bonds at the catalytic site or at a region that influences its configuration. Examination of the purified leucine aminopeptidase’s effect on the activation of the protoxin Cyt1Aa from Bti revealed that when it acts synergistically with Bti endogenous proteases, it has only a minor role in the processing of Cyt1Aa into an active toxin.     

Predicting the functionally distinct residues in the heme,cation, and substrate-binding sites of peroxidase from stress-tolerant mangrove specie, <Emphasis Type="Italic">Avicennia marina</Emphasis>

Recent work was conducted to predict the structure of functionally distinct regions of Avicennia marina peroxidase (AP) by using the structural coordinates of barley grains peroxidase as the template. This enzyme is utilized by all living organisms in many biosynthetic or degradable processes and in defense against oxidative stress. The homology model showed some distinct structural changes in the heme, calcium, and substrate-binding regions. Val53 was found to be an important coordinating residue between distal calcium ion and the distal heme site while Ser176 is coordinated to the proximal histidine through Ala174 and Leu172. Different ionic and hydrogen-bonded interactions were also observed in AP. Analyses of various substrate–enzyme interactions revealed that the substrate-binding pocket is provided by the residues, His41, Phe70, Gly71, Asp138, His139, and Lys176; the later three residues are not conserved in the peroxidase family. We have also performed structural comparison of the A. marina peroxidase with that of two class III salt-sensitive species, peanut and soybean. Four loop regions were found to have largest structural deviation. The overall protein sequence was also analyzed for the presence of probable post-translational modification sites and the functional significance of these sites were outlined.     

Novel thermostable lipase from <Emphasis Type="Italic">Bacillus circulans</Emphasis> IIIB153: comparison with the mesostable homologue at sequence and structure level

Thermophilic Bacillus circulans IIIB153 isolated from hot springs of North West Himalayas, India, produced an extracellular lipase, which exhibited significant biofilm disruption property on the static biofilm disruption model with a single species of Actinomyces viscosous. The gene encoding the lipase was cloned and overexpressed in Escherichia coli. Recombinant Bacillus circulans lipase (BCL), a monomer with molecular mass of 43 kDa also exhibited significant biofilm disruption activity. The enzyme was optimally active at 60°C, pH 8.5 and retained >70% of its original activity after 1 h incubation at 60°C. 3D structure of BCL developed by homology modeling showed a typical α/β hydrolase fold, a characteristic feature of lipolytic enzymes. Comparison of thermostable BCL with mesostable lipase from Chromobacterium viscosum at the sequence and structure level showed distinct variations in the structural features, with the presence of a high content of proline residues, aromatic amino acids and salt bridges. These features along with the presence of zinc-binding site observed in BCL structure could have a potential role in thermal stability of the enzyme.     

Structural stability and functional analysis of L-asparaginase from <Emphasis Type="Italic">Pyrococcus furiosus</Emphasis>

We report studies on an L-asparaginase from Pyrococcus furiosus, cloned and expressed in Escherichia coli and purified to homogeneity. Protein stability and enzyme kinetic parameters were determined. The enzyme was found to be thermostable, natively dimeric, and glutaminase-free, with optimum activity at pH 9.0. It showed a K _m of 12 mM and a substrate inhibition profile above 20 mM L-asparagine. Urea could not induce unfolding and enzyme inactivation; however, with guanidine hydrochloride (GdnCl) a two-state unfolding pattern was observed. Reduced activity and an altered near-UV-CD signal for protein at low GdnCl concentration (1 M) suggested tertiary structural changes at the enzyme active site. A homology three-dimensional model was developed and the structural information was combined with activity and stability data to give functional clues about the asparaginase.