Abstract: Antimicrobial peptides are drugs used against a wide range of pathogens which present a great advantage: in contrast with antibiotics they do not develop resistance. The wide spectrum of antimicrobial peptides advertises them in the research and pharmaceutical industry as attractive starting points for obtaining new, more effective analogs. Here we predict the antimicrobial activity against Bacillus subtilis (expressed as minimal inhibitory concentration values) for 33 mastoparan analogs and their new derivatives by a non-aligned 3D-QSAR (quantitative structure-activity relationship) method. We establish the contribution to antimicrobial activity of molecular descriptors (hydrophobicity, hydrogen bond donor and steric), correlated with contributions from the membrane environment (sodium, potassium, chloride). Our best QSAR models show significant cross-validated correlation q(2) (0.55-0.75), fitted correlation r(2) (greater than 0.90) coefficients and standard error of prediction SDEP (less than 0.250). Moreover, based on our most accurate 3D-QSAR models, we propose nine new mastoparan analogs, obtained by computational mutagenesis, some of them predicted to have significantly improved antimicrobial activity compared to the parent compound.
Abstract: Antidepressants and antipsychotics are psychiatric agents used for the treatment of various types of psychiatric diseases. Although currently among the most commonly prescribed drugs, their effectiveness and adverse effects are the topic of many studies and controversial claims. Here we generate QSAR models based on compounds series including 20 drugs recommended for two critical psychiatric diseases: depression and schizophrenia and we use these QSAR models to predict the biological activity of these 20 antidepressants and antipsychotics. We establish the membrane ions' contributions (sodium, potassium, calcium and iron) mediated by water to the antagonism of these drugs at the 5-HT1A receptor. The reliability of our QSAR models in predicting compounds activity is indicated by significant values for cross-validated correlation q2 (0.60–0.76) and fitted correlation r2 (0.96–0.98) coefficients. Our results indicate that potassium, calcium and iron play a key role for the antagonistic activity of drugs at the 5-HT1A receptor. Moreover, based on the established QSAR equations, we analysed 24 new escitalopram derivatives as possibly improved antidepressants targeting the 5-HT1A receptor. We identified that the presence of methyl groups and hydrogen atoms improves antidepressant activity while the simultaneous presence of ethyl, propyl or halogens decreased drastically antidepressant activity at the 5-HT1A site.
Abstract: TREK-1, a mechanosensitive K channel from the two-pore family (K(2)P), is involved in protective regulation of the resting potential in CNS neurons and other tissues. The structure of TREK-1 and the basis of its sensitivity to stretch and variety of lipid-soluble factors remain unknown. Using existing K channel structures as modeling templates, TREK-1 was envisioned as a two-fold symmetrical complex with the gate formed primarily by the centrally positioned TM2b helices of the second homologous repeat. Opening was modeled as a conical expansion of the barrel separating TM2b's accompanied by extension of TM2a helices with the cytoplasmic TM2a-TM1b connector. Seeking first experimental support to the models we have accomplished thermodynamic analysis of mouse TREK-1 gating and functional testing of several deletion mutants. The predicted increase of the channel in-plane area (ΔA) of ~5 nm(2) in models was supported by the experimental ΔA of ~4 nm(2) derived from the slope of open probability versus membrane tension in HEK-293T cells and their cytoskeleton-depleted blebs. In response to steps of suction, wild-type channel produced transient currents in cell-attached patches and mostly sustained currents upon patch excision. TREK-1 motifs not present in canonical K channels include divergent cytoplasmic N- and C-termini, and a characteristic 50-residue extracellular loop in the first homologous repeat. Deletion of the extracellular loop (Δ76-124) reduced the average current density in patches, increased spontaneous activity and generated a larger sub-population of high-conductance channels, while activation by tension augmented by arachidonic acid was fully retained. Further deletion of the C-terminal end (Δ76-124/Δ334-411) removed voltage dependency but otherwise produced no additional effect. In an attempt to generate a cysteine-free version of the channel, we mutated two remaining cysteines 159 and 219 in the transmembrane region. C219A did not compromise channel activity, whereas the C159A/S mutants were essentially inactive. Treatment with β-mercaptoethanol suggested that none of these cysteines form functionally-important disulfides.
Abstract: Mechanosensitive TREK channels belong to the family of K2P channels, a family of widely distributed, well modulated channels that uniquely have two similar or identical subunits, each with two TM1-P-TM2 motifs. Our goal is to build viable structural models of TREK channels, as representatives of K2P channels family. The structures available to be used as templates belong to the 2TM channels superfamily. These have low sequence similarity and different structural features: four symmetrically arranged subunits, each having one TM1-P-TM2 motif. Our model building strategy used two subunits of the template (KcsA) to build one subunit of the target (TREK-1). Our models of the Closed channel were adjusted to differ substantially from those of the template, e.g., TM2 of the 2nd repeat is near the axis of the pore whereas TM2 of the 1st repeat is far from the axis. Segments linking the two repeats and immediately following the last TM segment were modeled ab initio as α-helices based on helical periodicities of hydrophobic and hydrophilic residues, highly conserved and poorly conserved residues, and statistically related positions from multiple sequence alignments. The models were further refined by two-fold symmetry-constrained MD simulations using a protocol we developed previously. We also built models of the Open state and suggest a possible tension-activated gating mechanism characterized by helical motion with two-fold symmetry. Our models are consistent with deletion/truncation mutagenesis and thermodynamic analysis of gating described in the accompanying paper.
Abstract: There are presented six new QSAR models, which are correlating antipsychotic activity ( pK i values at dopamine D1–D4 and serotonine (5-HT2C, 5-HT2A) receptors) with physicochemical parameters. A large data set of typical and atypical antipsychotics already approved for clinical treatment including as well representatives with new chemical structures which are exhibiting antipsychotic activity (tetrahydrofuran-, benzamide-, 3-aminoethyl-1-tetralones-, piperazine-, benzothiazepine- and pyrrolobenzazepine-derivatives) were incorporated within this study. The appropriate 2D and internal-3D molecular descriptors could be generated by use of the computational software MOE (Molecular Operating Environment). Significant q 2 (0.63–0.76) and r 2 (0.70–0.78) correlation coefficients were obtained, indicating that the established equations can be used within a wide range of binding constants ( pK i = 5 to 9.76). By use of these linear models an assembly of new aripiprazole structures could be established. Some of them are showing significantly improved antipsychotic activity in comparison with the parent compound.
Abstract: SUMMARY The cell wall, a strong extraprotoplasmic layer surrounding plant cells that mainly consists of a variety of polysaccharides, constitutes a major barrier for potential parasites. Plant-parasitic nematodes are well equipped to overcome this barrier as they produce and secrete cell-wall-degrading enzymes. Expression profiling of various life stages of the potato cyst nematode Globodera rostochiensis revealed a novel pectate lyase gene (Gr-pel2, 759 bp). The Gr-PEL2 protein showed highest similarity to pectate lyases from the facultative plant-parasitic nematodes Bursaphelenchus mucronatus and B. xylophilus and the soil-inhabiting saprophytic Streptomyces and Frankia species (i.e. 40-42% identity and 58-60% similarity), whereas only a remote relatedness to the previously identified Gr-PEL1 was observed (i.e. 28% identity and 43% similarity). Transient expression of Gr-pel2 in leaves of Nicotiana benthamiana resulted in severe malformations of the infiltrated tissues, not relating to maceration and soft rot symptoms. Ca(2+) is known to be essential for pectate lyase activity, and the most likely calcium-binding site was identified in the Gr-PEL2 protein by combining homology modelling of the three-dimensional structure, site-directed mutagenesis and transient expression in leaves. A highly charged cleft in Gr-PEL2, which is likely to be involved in substrate binding and which is also significantly more hydrophobic in Gr-PEL1, was shown to be essential for protein activity. Our results underline the broad spectrum of pectate lyases and cell-wall-degrading enzymes necessary for successful parasitism by cyst nematodes.
Abstract: O-Glycan biosynthesis is initiated by the transfer of N-acetylgalactosamine (GalNAc) from a nucleotide sugar donor (UDP-GalNAc) to Ser/Thr residues of an acceptor substrate. The detailed transfer mechanism, catalyzed by the UDP-GalNAc polypeptide:N-acetyl-alpha-galactosaminyltransferases (ppGalNAcTs), remains unclear despite structural information available for several isoforms in complex with substrates at various stages along the catalytic pathway. We used all-atom molecular dynamics simulations with explicit solvent and counterions to study the conformational dynamics of ppGalNAcT-2 in several enzymatic states along the catalytic pathway. ppGalNAcT-2 is simulated both in the presence and in the absence of substrates and reaction products to examine the role of conformational changes in ligand binding. In multiple 40-ns-long simulations of more than 600 ns total run time, we studied systems ranging from 45,000 to 95,000 atoms. Our simulations accurately identified dynamically active regions of the protein, as previously revealed by the X-ray structures, and permitted a detailed, atomistic description of the conformational changes of loops near the active site and the characterization of the ensemble of structures adopted by the transferase complex on the transition pathway between the ligand-bound and ligand-free states. In particular, the conformational transition of a functional loop adjacent to the active site from closed (active) to open (inactive) is correlated with the rotameric state of the conserved residue W331. Analysis of water dynamics in the active site revealed that internal water molecules have an important role in enhancing the enzyme flexibility. We also found evidence that charged side chains in the active site rearrange during site opening to facilitate ligand binding. Our results are consistent with the single-displacement transfer mechanism previously proposed for ppGalNAcTs based on X-ray structures and mutagenesis data and provide new evidence for possible functional roles of certain amino acids conserved across several isoforms.
Abstract: We present here a neural networks method designed to predict biological activity based on a local representation of the ligand. The compounds of the series are represented by a vector mapping for each of four substituent properties: volume, log P, dipole moment and a simple 'steric' parameter relating to its shape. This ligand representation was tested using neural networks on a set of 42 cyclic-urea derivatives, inhibiting HIV-1 protease. The leave-one-out cross-validation using all descriptors in the input gave a correlation factor between prediction and experiment of 0.76 for the overall set and 0.88 when three outliers were left out. To rank the significance of the four descriptors, we further tested all combinations of two and three parameters for each substituent, using two disjunctive testing sets of five inhibitors. In these sets, vectors with extreme descriptor values were used either in the training or the testing set (sets A and B, respectively). The method is a very good interpolator (set A, 95+/-2% accuracy) but a less effective extrapolator (set B, 85+/-2% accuracy). Generally, the combinations including the 'steric' parameter predict better than average, while those containing the volume are less effective. The best prediction, 98.8+/-1.2%, was obtained when log P, the dipole and the steric parameter were used on set A. At the opposite end, the lowest ranked descriptor set was obtained when replacing log P with the volume, giving 92.3+/-6.7% accuracy over the set A.
Abstract: Esophageal secretions from endoparasitic sedentary nematodes are thought to play key roles throughout plant parasitism, in particular during the invasion of the root tissue and the initiation and maintenance of the nematode feeding site (NFS) essential for nematode development. The secretion in planta of esophageal cell-wall-degrading enzymes by migratory juveniles has been shown, suggesting a role for these enzymes in the invasion phase. Nevertheless, the secretion of an esophageal gland protein into the NFS by nematode sedentary stages has never been demonstrated. The calreticulin Mi-CRT is a protein synthesized in the esophageal glands of the root-knot nematode Meloidogyne incognita. After three-dimensional modeling of the Mi-CRT protein, a surface peptide was selected to raise specific antibodies. In planta immunolocalization showed that Mi-CRT is secreted by migratory and sedentary stage nematodes, suggesting a role for Mi-CRT throughout parasitism. During the maintenance of the NFS, the secreted Mi-CRT was localized outside the nematode at the tip of the stylet. In addition, Mi-CRT accumulation was observed along the cell wall of the giant cells that compose the feeding site, providing evidence for a nematode esophageal protein secretion into the NFS.
Abstract: Dopachrome tautomerase (Dct) is a type I membrane protein and an important regulatory enzyme that plays a pivotal role in the biosynthesis of melanin and in the rapid metabolism of its toxic intermediates. Dct-mutant melanocytes carrying the slaty or slaty light mutations were derived from the skin of newborn congenic C57BL/6J non-agouti black mice and were used to study the effect(s) of these mutations on the intracellular trafficking of Dct and on the pigmentation of the cells. Dct activity is 3-fold lower in slaty cells compared with non-agouti black melanocytes, whereas slaty light melanocytes have a surprisingly 28-fold lower Dct activity. Homology modelling of the active site of Dct suggests that the slaty mutation [R194Q (Arg194-->Gln)] is located in the active site and may alter the ability of the enzyme to transform the substrate. Transmembrane prediction methods indicate that the slaty light mutation [G486R (Gly486-->Arg)] may result in the sliding of the transmembrane domain towards the N-terminus, thus interfering with Dct function. Chemical analysis showed that both Dct mutations increase pheomelanin and reduce eumelanin produced by melanocytes in culture. Thus the enzymatic activity of Dct may play a role in determining whether the eumelanin or pheomelanin pathway is preferred for pigment biosynthesis.
Abstract: Southern analysis showed that Gr-EXPB1, a functional expansin from the potato cyst nematode Globodera rostochiensis, is member of a multigene family, and EST data suggest expansins to be present in other plant parasitic nematodes as well. Homology modeling predicted that Gr-EXPB1 domain 1 (D1) has a flat beta-barrel structure with surface-exposed aromatic rings, whereas the 3D structure of Gr-EXPB1-D2 was remarkably similar to plant expansins. Gr-EXPB1 shows highest sequence similarity to two extracellular proteins from saprophytic soil-inhabiting Actinobacteria, and includes a bacterial type II carbohydrate-binding module. These results support the hypothesis that a number of pathogenicity factors of cyst nematodes is of procaryotic origin and were acquired by horizontal gene transfer.
Abstract: We recently reported statistical analysis of structural data on glycosidic linkages. Here we extend this analysis to the glycan-protein linkage, and the peptide primary, secondary, and tertiary structures around N-glycosylation sites. We surveyed 506 glycoproteins in the Protein Data Bank crystallographic database, giving 2592 glycosylation sequons (1683 occupied) and generated a database of 626 nonredundant sequons with 386 occupied. Deviations in the expected amino acid composition were seen around occupied asparagines, particularly an increased occurrence of aromatic residues before the asparagine and threonine at position +2. Glycosylation alters the asparagine side chain torsion angle distribution and reduces its flexibility. There is an elevated probability of finding glycosylation sites in which secondary structure changes. An 11-class taxonomy was developed to describe protein surface geometry around glycosylation sites. Thirty-three percent of the occupied sites are on exposed convex surfaces, 10% in deep recesses and 20% on the edge of grooves with the glycan filling the cleft. A surprisingly large number of glycosylated asparagine residues have a low accessibility. The incidence of aromatic amino acids brought into close contact with the glycan by the folding process is higher than their normal levels on the surface or in the protein core. These data have significant implications for control of sequon occupancy and evolutionary selection of glycosylation sites and are discussed in relation to mechanisms of protein fold stabilization and regional quality control of protein folding. Hydrophobic protein-glycan interactions and the low accessibility of glycosylation sites in folded proteins are common features and may be critical in mediating these functions.
