Exploring the interaction energies for the binding of hydroxydiphenyl ethers to enoyl-acyl carrier protein reductases

It is now well established that the potent anti-microbial compound, triclosan, interrupts the type II fatty acid synthesis by inhibiting the enzyme enoyl-ACP reductase in a number of organisms. Existence of a high degree of similarity between the recently discovered enoyl-ACP reductase from R falciparum and B. napus enzyme permitted building of a satisfactory model for the former enzyme that explained some of the key aspects of the enzyme such as its specificity for binding to the cofactor and the inhibitor. We now report the interaction energies between triclosan and other hydroxydiphenyl ethers with the enzymes from B. napus, E. coli and R falciparum. Examination of the triclosan-enzyme interactions revealed that subtle differences exist in the ligand binding sites of the enzymes from different sources i.e., B. napus, E. coli and P falciparum. A comparison of their binding propensities thus determined should aid in the design of effective inhibitors for the respective enzymes.

Effect of Substrate Binding Loop Mutations on the Structure, Kinetics, and Inhibition of Enoyl Acyl Carrier Protein Reductase from Plasmodium falciparum

Enoyl acyl carrier protein reductase (ENR), which catalyzes the final and rate limiting step of fatty acid elongation, has been validated as a potential drug target. Triclosan is known to be an effective inhibitor for this enzyme. We mutated the substrate binding site residue Ala372 of the ENR of Plasmodium falciparum (PfENR) to Methionine and Valine which increased the affinity of the enzyme towards triclosan to almost double, close to that of Escherichia coli ENR (EcENR) which has a Methionine at the structurally similar position of Ala372 of PfENR. Kinetic studies of the mutants of PfENR and the crystal structure analysis of the A372M mutant revealed that a more hydrophobic environment enhances the affinity of the enzyme for the inhibitor. A triclosan derivative showed a threefold increase in the affinity towards the mutants compared to the wild type, due to additional interactions with the A372M mutant as revealed by the crystal structure. The enzyme has a conserved salt bridge which stabilizes the substrate binding loop and appears to be important for the active conformation of the enzyme. We generated a second set of mutants to check this hypothesis. These mutants showed loss of function, except in one case, where the crystal structure showed that the substrate binding loop is stabilized by a water bridge network. (C) 2011 IUBMB mum Life, 63(1): 30-41,2011

A soluble preparation from developing groundnut seeds (Arachis hypogaea) catalyzes de novo synthesis of long chain fatty acids

A 100,000 x g supernatant fraction prepared from developing groundnut seeds (30-35 days after flowering) catalyzed the synthesis of fatty acids from [l-14C]acetate at a rate of 120nmoles of acetate incorporated per hr per gram fresh weight of tissue. 90% of this incorporated label was associated with fatty acids. The major fatty acids formed were stearic- (77%) and palmitic acids (14%) with 4% of oleic acid. The fatty acid synthetase activity was stable when stored at 0-4 degrees C for at least fifteen days. It is concluded from these results that acetyl-coA carboxylase and all the enzymes of fatty acid synthetase from developing groundnut seeds are soluble.

Analysis of proteins with the `hot dog' fold: Prediction of function and identification of catalytic residues of hypothetical proteins

Background: The hot dog fold has been found in more than sixty proteins since the first report of its existence about a decade ago. The fold appears to have a strong association with fatty acid biosynthesis, its regulation and metabolism, as the proteins with this fold are predominantly coenzyme A-binding enzymes with a variety of substrates located at their active sites. Results: We have analyzed the structural features and sequences of proteins having the hot dog fold. This study reveals that though the basic architecture of the fold is well conserved in these proteins, significant differences exist in their sequence, nature of substrate and oligomerization. Segments with certain conserved sequence motifs seem to play crucial structural and functional roles in various classes of these proteins. Conclusion: The analysis led to predictions regarding the functional classification and identification of possible catalytic residues of a number of hot dog fold-containing hypothetical proteins whose structures were determined in high throughput structural genomics projects.

Hydrolysis of sphingomyelin to ceramide with hydrofluoric acid

Sphingomyelin was hydrolyzed to ceramide near quantitatively on treatment with 40% HF at 40°C for 72 hr. The reaction of sphingomyelin with HF is much slower than phosphoglycerides. HF treatment did not alter either the fatty acid composition or the stereochemical configuration of the sphingosine moiety of ceramide formed.

Ligand dependent intra and inter subunit communication in human tryptophanyl tRNA synthetase as deduced from the dynamics of structure networks

Homodimeric protein tryptophanyl tRNA synthetase (TrpRS) has a Rossmann fold domain and belongs to the 1c subclass of aminoacyl tRNA synthetases. This enzyme performs the function of acylating the cognate tRNA. This process involves a number of molecules (2 protein subunits, 2 tRNAs and 2 activated Trps) and thus it is difficult to follow the complex steps in this process. Structures of human TrpRS complexed with certain ligands are available. Based on structural and biochemical data, mechanism of activation of Trp has been speculated. However, no structure has yet been solved in the presence of both the tRNA(Trp) and the activated Trp (TrpAMP). In this study, we have modeled the structure of human TrpRS bound to the activated ligand and the cognate tRNA. In addition, we have performed molecular dynamics (MD) simulations on these models as well as other complexes to capture the dynamical process of ligand induced conformational changes. We have analyzed both the local and global changes in the protein conformation from the protein structure network (PSN) of MD snapshots, by a method which was recently developed in our laboratory in the context of the functionally monomeric protein, methionyl tRNA synthetase. From these investigations, we obtain important information such as the ligand induced correlation between different residues of this protein, asymmetric binding of the ligands to the two subunits of the protein as seen in the crystal structure analysis, and the path of communication between the anticodon region and the aminoacylation site. Here we are able to elucidate the role of dimer interface at a level of detail, which has not been captured so far.

Ligand Specificity of Group I Biotin Protein Ligase of Mycobacterium tuberculosis

Background: Fatty acids are indispensable constituents of mycolic acids that impart toughness & permeability barrier to the cell envelope of M. tuberculosis. Biotin is an essential co-factor for acetyl-CoA carboxylase (ACC) the enzyme involved in the synthesis of malonyl-CoA, a committed precursor, needed for fatty acid synthesis. Biotin carboxyl carrier protein (BCCP) provides the co-factor for catalytic activity of ACC. Methodology/Principal Findings: BPL/BirA (Biotin Protein Ligase), and its substrate, biotin carboxyl carrier protein (BCCP) of Mycobacterium tuberculosis (Mt) were cloned and expressed in E. coli BL21. In contrast to EcBirA and PhBPL, the similar to 29.5 kDa MtBPL exists as a monomer in native, biotin and bio-5'AMP liganded forms. This was confirmed by molecular weigt profiling by gel filtration on Superdex S-200 and Dynamic Light Scattering (DLS). Computational docking of biotin and bio-5'AMP to MtBPL show that adenylation alters the contact residues for biotin. MtBPL forms 11 H-bonds with biotin, relative to 35 with bio-5'AMP. Docking simulations also suggest that bio-5'AMP hydrogen bonds to the conserved `GRGRRG' sequence but not biotin. The enzyme catalyzed transfer of biotin to BCCP was confirmed by incorporation of radioactive biotin and by Avidin blot. The K-m for BCCP was similar to 5.2 mu M and similar to 420 nM for biotin. MtBPL has low affinity (K-b = 1.06 x 10(-6) M) for biotin relative to EcBirA but their K-m are almost comparable suggesting that while the major function of MtBPL is biotinylation of BCCP, tight binding of biotin/bio-5'AMP by EcBirA is channeled for its repressor activity. Conclusions/Significance: These studies thus open up avenues for understanding the unique features of MtBPL and the role it plays in biotin utilization in M. tuberculosis.

Structures and molecular-dynamics studies of three active-site mutants of bovine pancreatic phospholipase A(2)

Phospholipase A(2) hydrolyzes phospholipids at the sn-2 position to cleave the fatty-acid ester bond of L-glycerophospholipids. The catalytic dyad (Asp99 and His48) along with a nucleophilic water molecule is responsible for enzyme hydrolysis. Furthermore, the residue Asp49 in the calcium-binding loop is essential for controlling the binding of the calcium ion and the catalytic action of phospholipase A2. To elucidate the structural role of His48 and Asp49, the crystal structures of three active-site single mutants H48N, D49N and D49K have been determined at 1.9 angstrom resolution. Although the catalytically important calcium ion is present in the H48N mutant, the crystal structure shows that proton transfer is not possible from the catalytic water to the mutated residue. In the case of the Asp49 mutants, no calcium ion was found in the active site. However, the tertiary structures of the three active-site mutants are similar to that of the trigonal recombinant enzyme. Molecular-dynamics simulation studies provide a good explanation for the crystallographic results.

Crystal structure of dimeric FabZ of Plasmodium falciparum reveals conformational switching to active hexamers by peptide flips

The crystal structure of beta-hydroxyacyl acyl carrier protein dehydratase of Plasmodium falciparum (PfFabZ) has been determined at a resolution of 2.4 angstrom. PfFabZ has been found to exist as a homodimer (d-PfFabZ) in the crystals of the present study in contrast to the reported hexameric form (h-PfFabZ) which is a trimer of dimers crystallized in a different condition. The catalytic sites of this enzyme are located in deep narrow tunnel-shaped pockets formed at the dimer interface. A histidine residue from one subunit of the dimer and a glutamate residue from the other subunit lining the tunnel form the catalytic dyad in the reported crystal structures. While the position of glutamate remains unaltered in the crystal structure of d-PffabZ compared to that in b-PfFabZ, the histidine residue takes up an entirely different conformation and moves away from the tunnel leading to a His-Phe cis-trans peptide flip at the histidine residue. In addition, a loop in the vicinity has been observed to undergo a similar flip at a Tyr-Pro peptide bond. These alterations not only prevent the formation of a hexamer but also distort the active site geometry resulting in a dimeric form of FabZ that is incapable of substrate binding. The dimeric state and an altered catalytic site architecture make d-PfFabZ distinctly different from the FabZ structures described so far. Dynamic light scattering and size exclusion chromatographic studies clearly indicate a pH-related switching of the dimers to active hexamers. (c) 2006 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserv.

Interfacial microrheology as a tool to study viscoelastic transitions in nanoconfined soft matter

We present a method to perform in situ microrheological measurements on monolayers of soft materials undergoing viscoelastic transitions under compression. Using the combination of a Langmuir trough mounted on the inverted microscope stage of a laser scanning confocal microscope we track the motion of individual fluorescent quantum dots partly dispersed in monolayers spread at the air-water interface. From the calculated mean square displacement of the probe particles and extending a well established scheme of the generalized Stokes-Einstein relation in bulk to the interface we arrive at the viscoelastic modulus for the respective monolayers as a function of surface density. Measurements on monolayers of glassy as well as nonglassy polymers and a standard fatty acid clearly show sensitivity of our technique to subtle variations, in the viscoelastic properties of the highly confined materials under compression. Evidence for possible spatial variations of such viscoelastic properties at a given surface density for the fatty acid monolayer is also provided.

Interfacial microrheology as a tool to study viscoelastic transitions in nanoconfined soft matter

We present a method to perform in situ microrheological measurements on monolayers of soft materials undergoing viscoelastic transitions under compression. Using the combination of a Langmuir trough mounted on the inverted microscope stage of a laser scanning confocal microscope we track the motion of individual fluorescent quantum dots partly dispersed in monolayers spread at the air-water interface. From the calculated mean square displacement of the probe particles and extending a well established scheme of the generalized Stokes-Einstein relation in bulk to the interface we arrive at the viscoelastic modulus for the respective monolayers as a function of surface density. Measurements on monolayers of glassy as well as nonglassy polymers and a standard fatty acid clearly show sensitivity of our technique to subtle variations, in the viscoelastic properties of the highly confined materials under compression. Evidence for possible spatial variations of such viscoelastic properties at a given surface density for the fatty acid monolayer is also provided.

Defective in Cuticular Ridges (DCR) of Arabidopsis thaliana, a Gene Associated with Surface Cutin Formation, Encodes a Soluble Diacylglycerol Acyltransferase

A key step in the triacylglycerol (TAG) biosynthetic pathway is the final acylation of diacylglycerol (DAG) by DAG acyltransferase. In silico analysis has revealed that the DCR (defective in cuticular ridges) (At5g23940) gene has a typical HX4D acyltransferase motif at the N-terminal end and a lipid binding motif VX(2)GF at the middle of the sequence. To understand the biochemical function, the gene was overexpressed in Escherichia coli, and the purified recombinant protein was found to acylate DAG specifically in an acyl-CoA-dependent manner. Overexpression of At5g23940 in a Saccharomyces cerevisiae quadruple mutant deficient in DAG acyltransferases resulted in TAG accumulation. At5g23940 rescued the growth of this quadruple mutant in the oleate-containing medium, whereas empty vector control did not. Lipid particles were localized in the cytosol of At5g23940-transformed quadruple mutant cells, as observed by oil red O staining. There was an incorporation of 16-hydroxyhexadecanoic acid into TAG in At5g23940-transformed cells of quadruple mutant. Here we report a soluble acyl-CoA-dependent DAG acyltransferase from Arabidopsis thaliana. Taken together, these data suggest that a broad specific DAG acyltransferase may be involved in the cutin as well as in the TAG biosynthesis by supplying hydroxy fatty acid.

Probing the role of the fully conserved Cys126 in triosephosphate isomerase by site-specific mutagenesis - distal effects on dimer stability

Cys126 is a completely conserved residue in triosephosphate isomerase that is proximal to the active site but has been ascribed no specific role in catalysis. A previous study of the C126S and C126A mutants of yeast TIM reported substantial catalytic activity for the mutant enzymes, leading to the suggestion that this residue is implicated in folding and stability [Gonzalez-Mondragon E et al. (2004) Biochemistry43, 3255–3263]. We re-examined the role of Cys126 with the Plasmodium falciparum enzyme as a model. Five mutants, C126S, C126A, C126V, C126M, and C126T, were characterized. Crystal structures of the 3-phosphoglycolate-bound C126S mutant and the unliganded forms of the C126S and C126A mutants were determined at a resolution of 1.7–2.1 Å. Kinetic studies revealed an approximately five-fold drop in kcat for the C126S and C126A mutants, whereas an approximately 10-fold drop was observed for the other three mutants. At ambient temperature, the wild-type enzyme and all five mutants showed no concentration dependence of activity. At higher temperatures (> 40 °C), the mutants showed a significant concentration dependence, with a dramatic loss in activity below 15 μm. The mutants also had diminished thermal stability at low concentration, as monitored by far-UV CD. These results suggest that Cys126 contributes to the stability of the dimer interface through a network of interactions involving His95, Glu97, and Arg98, which form direct contacts across the dimer interface.

Structural basis for the functional and inhibitory mechanisms of beta-hydroxyacyl-acyl carrier protein dehydratase (FabZ) of Plasmodium falciparum

The beta-hydroxyacyl-acyl carrier protein dehydratase of Plasmodium falciparum (PfFabZ) catalyzes the third and important reaction of the fatty acid elongation cycle. The crystal structure of PfFabZ is available in hexameric (active) and dimeric (inactive) forms. However, PfFabZ has not been crystallized with any bound inhibitors until now. We have designed a new condition to crystallize PfFabZ with its inhibitors bound in the active site, and determined the crystal structures of four of these complexes. This is the first report on any FabZ enzyme with active site inhibitors that interact directly with the catalytic residues. Inhibitor binding not only stabilized the substrate binding loop but also revealed that the substrate binding tunnel has an overall shape of ``U''. In the crystal structures, residue Phe169 located in the middle of the tunnel was found to be in two different conformations, open and closed. Thus, Phe169, merely by changing its side chain conformation, appears to be controlling the length of the tunnel to make it suitable for accommodating longer substrates. The volume of the substrate binding tunnel is determined by the sequence as well as by the conformation of the substrate binding loop region and varies between organisms for accommodating fatty acids of different chain lengths. This report on the crystal structures of the complexes of PfFabZ provides the structural basis of the inhibitory mechanism of the enzyme that could be used to improve the potency of inhibitors against an important component of fatty acid synthesis common to many infectious organisms. (C) 2011 Elsevier Inc. All rights reserved.

Deorphanization of Malonyl CoA:ACP Transacylase Drug Target in Plasmodium falciparum (PfFabD) Using Bacterial Antagonists: A Piggyback' Approach for Antimalarial Drug Discovery

Quest for new drug targets in Plasmodium sp. has underscored malonyl CoA:ACP transacylase (PfFabD) of fatty acid biosynthetic pathway in apicoplast. In this study, a piggyback approach was employed for the receptor deorphanization using inhibitors of bacterial FabD enzymes. Due to the lack of crystal structure, theoretical model was constructed using the structural details of homologous enzymes. Sequence and structure analysis has localized the presence of two conserved pentapeptide motifs: GQGXG and GXSXG and five key invariant residues viz., Gln109, Ser193, Arg218, His305 and Gln354 characteristic of FabD enzyme. Active site mapping of PfFabD using substrate molecules has disclosed the spatial arrangement of key residues in the cavity. As structurally similar molecules exhibit similar biological activities, signature pharmacophore fingerprints of FabD antagonists were generated using 0D-3D descriptors for molecular similarity-based cluster analysis and to correlate with their binding profiles. It was observed that antagonists showing good geometrical fitness score were grouped in cluster-1, whereas those exhibiting high binding affinities in cluster-2. This study proves important to shed light on the active site environment to reveal the hotspot for binding with higher affinity and to narrow down the virtual screening process by searching for close neighbors of the active compounds.