Molecular modeling, mutational analysis and steroid specificity of the ligand binding pocket of mPRα (PAQR7): Shared ligand binding with AdipoR1 and its structural basis.

The 7-transmembrane architecture of adiponectin receptors (AdipoRs), determined from their X-ray crystal structures, was used for homology modeling of another progesterone and adipoQ receptor (PAQR) family member, membrane progesterone receptor alpha (mPRα). The mPRα model identified excess positively charged residues on the cytosolic side, suggesting it has the same membrane orientation as AdipoRs with an intracellular N-terminus. The homology model showed identical amino acid residues to those forming the zinc binding pocket in AdipoRs, which strongly implies that zinc is also present in mPRα. The homology model showed a critical H-bond interaction between the glutamine (Q) residue at 206 in the binding pocket and the 20-carbonyl of progesterone. Mutational analysis showed no progesterone binding to the arginine (R) 206 mutant and modeling predicted this was due to the strong positive charge of arginine stabilizing the presence of an oleic acid (C18:1) molecule in the binding pocket, as observed in the X-rays of AdipoRs. High Zn2+ concentrations are predicted to form a salt with the carboxylate group of the oleic acid, thereby eliminating its binding to the free fatty acid (FFA) binding pocket, and allowing progesterone to bind. This is supported by experiments showing 100 µM Zn2+ addition restored [3H]-progesterone binding of the Q206R mutant to levels in WT mPRα and increased [3H]-progesterone binding to mPRγ and AdipoR1 which have arginine residues in this region. The model predicts hydrophobic interactions of progesterone with amino acid residues surrounding the binding pocket, including valine 146 in TM3, which when mutated into a polar serine resulted in a complete loss of [3H]-progesterone binding. The mPRα model showed there is no hydrogen bond donor in the vicinity of the 3-keto group of progesterone and ligand structure-activity studies with 3-deoxy steroids revealed that, unlike the nuclear progesterone receptor, the 3-carbonyl oxygen is not essential for binding to mPRα. Interestingly, the small synthetic AdipoR agonist, AdipoRon, displayed binding affinity for mPRα and mimicked progesterone signaling, whereas D-e-MAPP, a ceramidase inhibitor, blocked progesterone signaling. Thus, critical residues around the binding pocket and steroid structures that bind mPRα, as well as similarities with AdipoRs, can be predicted from the homology model.

KNIME-based Analysis of Off-Target Effect of Drugs Related to The Molecular 2D Fingerprint.

PURPOSE: Structurally similar molecules are likely to have similar biological activity. In this study, similarity searching based on molecular 2D fingerprint was performed to analyze off-target effects of drugs. The purpose of this study is to determine the correlation between the adverse effects and drug off-targets. METHODS: A workflow was built using KNIME to run dataset preparation of twenty-nine targets from ChEMBL, generate molecular 2D fingerprints of the ligands, calculate the similarity between ligand sets, and compute the statistical significance using similarity ensemble approach (SEA). Tanimoto coefficients (Tc) are used as a measure of chemical similarity in which the values between 0.2 and 0.4 are the most common for the majority of ligand pairs and considered to be insignificant similar. RESULT: The majority of ligand sets are unrelated, as is evidenced by the intrinsic chemical differences and the classification of statistical significance based on expectation value. The rank-ordered expectation value of inter-target similarity showed a correlation with off-target effects of the known drugs. CONCLUSION: Similarity-searching using molecular 2D fingerprint can be applied to predict off-targets and correlate them to the adverse effects of the drugs. KNIME as an open-source data analytic platform is applicable to build a workflow for data mining of ChEMBL database and generating SEA statistical model.

QSAR-based molecular signatures of prenylated (iso)flavonoids underlying antimicrobial potency against and membrane-disruption in Gram positive and Gram negative bacteria.

Prenylated flavonoids and isoflavonoids are phytochemicals with remarkable antibacterial activity. In this study, 30 prenylated (iso)flavonoids were tested against Listeria monocytogenes and Escherichia coli (the latter in combination with an efflux pump inhibitor). Minimum inhibitory concentrations of the most active compounds ranged between 6.3-15.0 µg/mL. Quantitative structure-activity relationships (QSAR) analysis was performed and linear regression models were proposed with R2 between 0.77-0.80, average R2m between 0.70-0.75, Q2LOO between 0.66-0.69, and relatively low amount of descriptors. Shape descriptors (related to flexibility and globularity), together with hydrophilic/hydrophobic volume and surface area descriptors, were identified as important molecular characteristics related to activity. A 3D pharmacophore model explaining the effect of the prenyl position on the activity of compounds was developed for each bacterium. These models predicted active compounds with an accuracy of 71-88%. With regard to the mode of action, good antibacterial prenylated (iso)flavonoids with low relative hydrophobic surface area caused remarkable membrane permeabilization, whereas those with higher relative hydrophobic surface area did not. Based on the QSAR and membrane permeabilization studies, the mode of action of antibacterial prenylated (iso)flavonoids was putatively rationalized.

Molden 2.0: quantum chemistry meets proteins.

Since the first distribution of Molden in 1995 and the publication of the first article about this software in 2000 work on Molden has continued relentlessly. A few of the many improved or fully novel features such as improved and broadened support for quantum chemistry calculations, preparation of ligands for use in drug design related softwares, and working with proteins for the purpose of ligand docking.

Structural basis for non-genuine phenolic acceptor substrate specificity of Streptomyces roseochromogenes prenyltransferase CloQ from the ABBA/PT-barrel superfamily.

Acceptor substrate specificity of Streptomyces roseochromogenes prenyltransferase SrCloQ was investigated using different non-genuine phenolic compounds. RP-UHPLC-UV-MSn was used for the tentative annotation and quantification of the prenylated products. Flavonoids, isoflavonoids and stilbenoids with different types of substitution were prenylated by SrCloQ, although with less efficiency than the genuine substrate 4-hydroxyphenylpyruvate. The isoflavan equol, followed by the flavone 7,4'-dihydroxyflavone, were the best non-genuine acceptor substrates. B-ring C-prenylation was in general preferred over A-ring C-prenylation (ratio 5:1). Docking studies of non-genuine acceptor substrates with the B-ring oriented towards the donor substrate dimethylallyl pyrophosphate, showed that the carbonyl group of the C-ring was able to make stabilizing interactions with the residue Arg160, which might determine the preference observed for B-ring prenylation. No reaction products were formed when the acceptor substrate had no phenolic hydroxyl groups. This preference can be explained by the essential hydrogen bond needed between a phenolic hydroxyl group and the residue Glu281. Acceptor substrates with an additional hydroxyl group at the C3' position (B-ring), were mainly O3'-prenylated (> 80% of the reaction products). This can be explained by the proximity of the C3' hydroxyl group to the donor substrate at the catalytic site. Flavones were preferred over isoflavones by SrCloQ. Docking studies suggested that the orientation of the B-ring and of the phenolic hydroxyl group at position C7 (A-ring) of flavones towards the residue Tyr233 plays an important role in this observed preference. Finally, the insights obtained on acceptor substrate specificity and regioselectivity for SrCloQ were extended to other prenyltransferases from the CloQ/NhpB family.

Visually impaired researchers get their hands on quantum chemistry: application to a computational study on the isomerization of a sterol.

In molecular sciences, articles tend to revolve around 2D representations of 3D molecules, and sighted scientists often resort to 3D virtual reality software to study these molecules in detail. Blind and visually impaired (BVI) molecular scientists have access to a series of audio devices that can help them read the text in articles and work with computers. Reading articles published in this journal, though, is nearly impossible for them because they need to generate mental 3D images of molecules, but the article-reading software cannot do that for them. We have previously designed AsteriX, a web server that fully automatically decomposes articles, detects 2D plots of low molecular weight molecules, removes meta data and annotations from these plots, and converts them into 3D atomic coordinates. AsteriX-BVI goes one step further and converts the 3D representation into a 3D printable, haptic-enhanced format that includes Braille annotations. These Braille-annotated physical 3D models allow BVI scientists to generate a complete mental model of the molecule. AsteriX-BVI uses Molden to convert the meta data of quantum chemistry experiments into BVI friendly formats so that the entire line of scientific information that sighted people take for granted-from published articles, via printed results of computational chemistry experiments, to 3D models-is now available to BVI scientists too. The possibilities offered by AsteriX-BVI are illustrated by a project on the isomerization of a sterol, executed by the blind co-author of this article (HBW).

Enzyme-specific activation versus leaving group ability.

Enzyme-specific activation and the substrate mimetics strategy are effective ways to circumvent the limited substrate recognition often encountered in protease-catalyzed peptide synthesis. A key structural element in both approaches is the guanidinophenyl (OGp) ester, which enables important interactions for affinity and recognition by the enzyme--at least, this is usually the explanation given for its successful application. In this study we show that leaving group ability is of equal or even greater importance. To this end we used both experimental and computational methods: 1) synthesis of close analogues of OGp, and their evaluation in a dipeptide synthesis assay with trypsin, 2) molecular docking studies to provide insights into the binding mode, and 3) ab initio calculations to evaluate their electronic properties.

Phenylalanine 368 of multidrug resistance-associated protein 4 (MRP4/ABCC4) plays a crucial role in substrate-specific transport activity.

Multidrug resistance-associated protein 4 (MRP4) is a membrane transporter that mediates the cellular efflux of a wide range of anionic drugs and endogenous molecules. MRP4 transport can influence the pharmacokinetics of drugs and their metabolites, therefore more knowledge about the molecular determinants important for its transport function would be of relevance. Here, we substituted amino acids Phe(368), Trp(995), and Arg(998) with conservative or non-conservative residues, and determined the effect on transport of the model substrates estradiol 17-ß-d-glucuronide (E(2)17ßG), cyclic guanosine monophosphate (cGMP), methotrexate (MTX), and folic acid into membrane vesicles isolated from baculovirus transduced HEK293 cells overexpressing the mutant MRP4 proteins. This revealed that all Arg(998) mutations appeared to be deleterious, whereas the effect of a Phe(368) or Trp(995) replacement was dependent on the amino acid introduced and the substrate studied. Substitution of Phe(368) with Trp (F368W) induced a gain-of-function of E(2)17ßG transport and a loss-of-function of MTX transport, which could not be attributed to an altered substrate binding. Moreover, we did not observe any modification in ATP or ADP handling for F368W. These results, in combination with docking of substrates in a homology model of MRP4 in the inward- and outward-facing conformation, suggest that Phe(368) and Trp(995) do not play an important role in the initial binding of substrates. They, however, might interact with the substrates during rearrangement of helixes for substrate translocation, funneling the substrates to the exit site in the outward-facing conformation.

Triazole fungicides can induce cross-resistance to medical triazoles in Aspergillus fumigatus.

BACKGROUND: Azoles play an important role in the management of Aspergillus diseases. Azole resistance is an emerging global problem in Aspergillus fumigatus, and may develop through patient therapy. In addition, an environmental route of resistance development has been suggested through exposure to 14α-demethylase inhibitors (DMIs). The main resistance mechanism associated with this putative fungicide-driven route is a combination of alterations in the Cyp51A-gene (TR(34)/L98H). We investigated if TR(34)/L98H could have developed through exposure to DMIs. METHODS AND FINDINGS: Thirty-one compounds that have been authorized for use as fungicides, herbicides, herbicide safeners and plant growth regulators in The Netherlands between 1970 and 2005, were investigated for cross-resistance to medical triazoles. Furthermore, CYP51-protein homology modeling and molecule alignment studies were performed to identify similarity in molecule structure and docking modes. Five triazole DMIs, propiconazole, bromuconazole, tebuconazole, epoxiconazole and difenoconazole, showed very similar molecule structures to the medical triazoles and adopted similar poses while docking the protein. These DMIs also showed the greatest cross-resistance and, importantly, were authorized for use between 1990 and 1996, directly preceding the recovery of the first clinical TR(34)/L98H isolate in 1998. Through microsatellite genotyping of TR(34)/L98H isolates we were able to calculate that the first isolate would have arisen in 1997, confirming the results of the abovementioned experiments. Finally, we performed induction experiments to investigate if TR(34)/L98H could be induced under laboratory conditions. One isolate evolved from two copies of the tandem repeat to three, indicating that fungicide pressure can indeed result in these genomic changes. CONCLUSIONS: Our findings support a fungicide-driven route of TR(34)/L98H development in A. fumigatus. Similar molecule structure characteristics of five triazole DMIs and the three medical triazoles appear the underlying mechanism of cross resistance development. Our findings have major implications for the assessment of health risks associated with the use of triazole DMIs.

Quantum mechanical polar surface area.

A correlation has been established between the absorbed fraction of training-set molecules after oral administration in humans and the Quantum Mechanical Polar Surface Area (QMPSA). This correlation holds for the QMPSA calculated with structures where carboxyl groups are deprotonated. The correlation of the absorbed fraction and the QMPSA calculated on the neutral gas phase optimized structures is much less pronounced. This suggests that the absorption process is mainly determined by polar interactions of the drug molecules in water solution. Rules are given to derive the optimal polar/apolar ranges of the electrostatic potential.

The structure-function relationship of the Aspergillus fumigatuscyp51A L98H conversion by site-directed mutagenesis: the mechanism of L98H azole resistance.

Since 1998, the rapid emergence of multi-azole-resistance (MAR) was observed in Aspergillus fumigatus in the Netherlands. Two dominant mutations were found in the cyp51A gene, a 34bp tandem repeat (TR) in the promoter region combined with a leucine to histidine substitution at codon 98 (L98H). In this study, we show that molecular dynamics simulations combined with site-directed mutagenesis of amino acid substitutions in the cyp51A gene, correlate to the structure-function relationship of the L98H substitution conferring to MAR in A. fumigatus. Because of a L98H directed change in the flexibility of the loops, that comprise a gate-like structure in the protein, the capacity of the two ligand entry channels is modified by narrowing the diameter and thereby binding of azoles is obstructed. Moreover, the L98H induced relocation of tyrosine 121 and tyrosine 107 seems to be related to the MAR phenotype, without affecting the biological activity of the CYP51A protein. Site-directed mutagenesis showed that both the 34bp TR and the L98H mutation are required to obtain the MAR phenotype. Furthermore, the amino acid leucine in codon 98 in A. fumigatus is highly conserved and important for maintaining the structure of the CYP51A protein that is essential for azole docking.

Azole resistance profile of amino acid changes in Aspergillus fumigatus CYP51A based on protein homology modeling.

Molecular studies have shown that the majority of azole resistance in Aspergillus fumigatus is associated with amino acid substitutions in the cyp51A gene. To obtain insight into azole resistance mutations, the cyp51A gene of 130 resistant and 76 susceptible A. fumigatus isolates was sequenced. Out of 130 azole-resistant isolates, 105 contained a tandem repeat of 34 bp in the promoter region and a leucine-to-histidine substitution in codon 98 (designated TR/L98H). Additionally, in 12 of these TR/L98H resistant isolates, the mutations S297T and F495I were found, and in 1 isolate, the mutation F495I was found. In eight azole-resistant isolates, known azole resistance mutations were detected in codon G54, G138, or M220. In three azole-susceptible isolates, the mutation E130D, L252L, or S400I was found and in 13 azole-susceptible isolates but also in 1 azole-resistant isolate, the mutations F46Y, G98G, M172V, N248T, D255E, L358L, E427K, and C454C were found. All of the nonsynonymous mutations, apart from the mutations in codons G54, G138, and M220 and L98H, were located at the periphery of the protein, as determined by a structural model of the A. fumigatus Cyp51A protein, and were predicted neither to interact with azole compounds nor to affect structural integrity. Therefore, this wide diversity of mutations in the cyp51A gene in azole-susceptible A. fumigatus isolates is not correlated with azole resistance. Based on the Cyp51A protein homology model, the potential correlation of a mutation to azole resistance can be predicted.

The non-gastric H,K-ATPase as a tool to study the ouabain-binding site in Na,K-ATPase.

Based on studies with chimeras between (non-)gastric H,K-ATPase and Na,K-ATPase, a model for the ouabain binding site has recently been presented (Qiu et al. J.Biol.Chem. 280 (2005) 32349). In this model, hydrogen bonds between specific amino acid residues of Na,K-ATPase and hydroxyl groups of ouabain play a crucial role. In the present study, a series of ouabain analogues were tested on baculovirus-expressed Na,K-ATPase and an ouabain-sensitive mutant of non-gastric H,K-ATPase (D312E/ S319G/ A778P/ I795L/ F802C). For each analogue, the results obtained by measuring ATPase inhibition and [(3)H]ouabain replacement agreed rather well. In Na,K-ATPase, strophanthidin had a 7-10 times higher and digoxin a 4-12 times lower affinity than ouabain. The results of the non-gastric H,K-ATPase mutant were rather similar to that of Na,K-ATPase with exception of dihydro-ouabain that showed a much lower affinity with the non-gastric H,K-ATPase mutant. Docking studies showed that all analogues bind to the same pocket in Na,K-ATPase. However, the amino acids to which hydrogen bonds were formed differed and depended on the availability of hydroxyl or keto groups in the ouabain analogues.

QCRNA 1.0: a database of quantum calculations for RNA catalysis.

This work outlines a new on-line database of quantum calculations for RNA catalysis (QCRNA) available via the worldwide web at http://theory.chem.umn.edu/QCRNA. The database contains high-level density functional calculations for a large range of molecules, complexes and chemical mechanisms important to phosphoryl transfer reactions and RNA catalysis. Calculations are performed using a strict, consistent protocol such that a wealth of cross-comparisons can be made to elucidate meaningful trends in biological phosphate reactivity. Currently, around 2000 molecules have been collected in varying charge states in the gas phase and in solution. Solvation was treated with both the PCM and COSMO continuum solvation models. The data can be used to study important trends in reactivity of biological phosphates, or used as benchmark data for the design of new semiempirical quantum models for hybrid quantum mechanical/molecular mechanical simulations.

Reconstruction of the complete ouabain-binding pocket of Na,K-ATPase in gastric H,K-ATPase by substitution of only seven amino acids.

Although cardiac glycosides have been used as drugs for more than 2 centuries and their primary target, the sodium pump (Na,K-ATPase), has already been known for 4 decades, their exact binding site is still elusive. In our efforts to define the molecular basis of digitalis glycosides binding we started from the fact that a closely related enzyme, the gastric H,K-ATPase, does not bind glycosides like ouabain. Previously, we showed that a chimera of these two enzymes, in which only the M3-M4 and M5-M6 hairpins were of Na,K-ATPase, bound ouabain with high affinity (Koenderink, J. B., Hermsen, H. P. H., Swarts, H. G. P., Willems, P. H. G. M., and De Pont, J. J. H. H. M. (2000) Proc. Natl. Acad. Sci. U. S. A. 97, 11209-11214). We also demonstrated that only three amino acids (Phe(783), Thr(797), and Asp(804)) present in the M5-M6 hairpin of Na,K-ATPase were sufficient to confer high affinity ouabain binding to a chimera which contained in addition the M3-M4 hairpin of Na,K-ATPase (Qiu, L. Y., Koenderink, J. B., Swarts, H. G., Willems, P. H., and De Pont, J. J. H. H. M. (2003) J. Biol. Chem. 278, 47240-47244). To further pinpoint the ouabain-binding site here we used a chimera-based loss-of-function strategy and identified four amino acids (Glu(312), Val(314), Ile(315), Gly(319)), all present in M4, as being important for ouabain binding. In a final gain-of-function study we showed that a gastric H,K-ATPase that contained Glu(312), Val(314), Ile(315), Gly(319), Phe(783), Thr(797), and Asp(804) of Na,K-ATPase bound ouabain with the same affinity as the native enzyme. Based on the E(2)P crystal structure of Ca(2+)-ATPase we constructed a homology model for the ouabain-binding site of Na,K-ATPase involving all seven amino acids as well as several earlier postulated amino acids.