Molecular modelling of miraculin: Structural analyses and functional hypotheses

Miraculin is a plant protein that displays the peculiar property of modifying taste by swiching sour into a sweet taste. Its monomer is flavourless at all pH as well as at high concentration; the dimer form elicits its taste-modifying activity at acidic pH; a tetrameric form is also reported as active. Two histidine residues, located in exposed regions, are the main responsible of miraculin activity, as demonstrated by mutagenesis studies. Since structural data of miraculin are not available, we have predicted its three-dimensional structure and simulated both its dimer and tetramer forms by comparative modelling and molecular docking techniques. Finally, molecular dynamics simulations at different pH conditions have indicated that at acidic pH the dimer assumes a widely open conformation, in agreement with the hypotheses coming from other studies.     

Computational analyses of curcuminoid analogs against kinase domain of HER2

Background

Human epidermal growth factor receptor 2 (HER2) has an important role in cancer aggressiveness and poor prognosis. HER2 has been used as a drug target for cancers. In particular, to effectively treat HER2-positive cancer, small molecule inhibitors were developed to target HER2 kinase. Knowing that curcumin has been used as food to inhibit cancer activity, this study evaluated the efficacy of natural curcumins and curcumin analogs as HER2 inhibitors using in vitro and in silico studies. The curcumin analogs considered in this study composed of 4 groups classified by their core structure, β-diketone, monoketone, pyrazole, and isoxazole.

Results

In the present study, both computational and experimental studies were performed. The specificity of curcumin analogs selected from the docked results was examined against human breast cancer cell lines. The screened curcumin compounds were then subjected to molecular dynamics simulation study. By modifying curcumin analogs, we found that protein-ligand affinity increases. The benzene ring with a hydroxyl group could enhance affinity by forming hydrophobic interactions and the hydrogen bond with the hydrophobic pocket. Hydroxyl, carbonyl or methoxy group also formed hydrogen bonds with residues in the adenine pocket and sugar pocket of HER2-TK. These modifications could suggest the new drug design for potentially effective HER2-TK inhibitors. Two outstanding compounds, bisdemethylcurcumin (AS-KTC006) and 3,5-bis((E)-3,4-dimethoxystyryl)isoxazole (AS-KTC021 ),were well oriented in the binding pocket almost in the simulation time, 30 ns. This evidence confirmed the results of cell-based assays and the docking studies. They possessed more distinguished interactions than known HER2-TK inhibitors, considering them as a promising drug in the near future.

Conclusions

The series of curcumin compounds were screened using a computational molecular docking and followed by human breast cancer cell lines assay. Both AS-KTC006 and AS-KTC021 could inhibit breast cancer cell lines though inhibiting of HER2-TK. The intermolecular interactions were confirmed by molecular dynamics simulation studies. This information would explore more understanding of curcuminoid structures and HER2-TK.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-261) contains supplementary material, which is available to authorized users.     

肝癌特异性鼠源及人源化单链抗体基因的构建及大肠杆菌中的表达

为探讨一株肝细胞癌特异性鼠源及其人源化单链抗体基因在大肠杆菌中的可溶性表达策略并比较二者对抗原的结构能力,在三种载体中分别以融合、分泌及胞内表达的方式进行了研究,表达产物均以包涵体形式存在;对复性后的单链抗体以细胞ＥＬＩＳＡ及竞争抑制流式细胞仪法进行检测,表明人源化单链抗体和鼠源单链抗体有相近的抗原结合能力。结论是：大肠杆菌中表达的基因工程单链抗体的可溶性可能主要由自身氨基酸一级序列决定;先前的设     

Evolutionarily evolved discriminators in the 3-TPR domain of the Toc64 family involved in protein translocation at the outer membrane of chloroplasts and mitochondria

Transport of polypeptides across membranes is a general and essential cellular process utilised by molecular machines. At least one component of these complexes contains a domain composed of three tetratricopeptide repeat (3-TPR) motifs. We have focussed on the receptor Toc64 to elucidate the evolved functional specifications of its 3-TPR domain. Toc64 is a component of the Toc core complex and functionally replaces Tom70 at the outer membrane of mitochondria in plants. Its 3-TPR domain recognises the conserved C-terminus of precursor-bound chaperones. We built homology models of the 3-TPR domain of chloroplastic Toc64 from different species and of the mitochondrial isoform from Arabidopsis. Guided by modelling, we identified residues essential for functional discrimination of the differently located isoforms to be located almost exclusively on the convex surface of the 3-TPR domain. The only exception is at568Ser/ps557Met, which is positioned in the ligand-binding groove. The functional implications of the homology models are discussed. Figure Motion contained within the 2nd eigenvector of the Calpha covariance matrix of the 3-TPR domain of atToc64-V indicated by a porcupine plot Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Homology modeling and dynamics study of aureusidin synthase--an important enzyme in aurone biosynthesis of snapdragon flower

Aurones, a class of plant flavonoids, provide bright yellow color on some important ornamental flowers, such as cosmos, coreopsis, and snapdragon (Antirrhinum majus). Recently, it has been elucidated that aureusidin synthase (AUS), a homolog of plant polyphenol oxidase (PPO), plays a key role in the yellow coloration of snapdragon flowers. In addition, it has been shown that AUS is a chalcone-specific PPO specialized for aurone biosynthesis. AUS gene has been successfully demonstrated as an attractive tool to engineer yellow flowers in blue flowers. Despite these biological studies, the structural basis for the specificity of substrate interactions of AUS remains elusive. In this study, we performed homology modeling of AUS using Grenache PPO and Sweet potato catechol oxidase (CO). An AUS-inhibitor was then developed from the initial homology model based on the CO and subsequently validated. We performed a thorough study between AUS and PTU inhibitor by means of interaction energy, which indicated the most important residues in the active site that are highly conserved. Analysis of the molecular dynamics simulations of the apo enzyme and ligand-bound complex showed that complex is relatively stable than apo and the active sites of both systems are flexible. The results from this study provide very helpful information to understand the structure-function relationships of AUS.     

Computational development of an alpha1A-adrenoceptor model in a membrane mimic

A number of subtypes of the alpha-adrenoceptor have been identified; however, knowledge of the three-dimensional structures of such membrane proteins is limited, and no crystal structure of an alpha-adrenoceptor is available to date. We have developed and analysed homology models of the alpha1A-adrenoceptor subtype based on the bovine rhodopsin crystal structure (1l9 h). Subsequent structural refinement was performed through molecular dynamics simulations using the Amber 7 suite of programs with a biphasic H2O/CHCl3/H2O cell utilised to mimic the receptor's natural membrane environment.     

Setting up and running molecular dynamics simulations of membrane proteins

Molecular dynamics simulations have become a popular and powerful technique to study lipids and membrane proteins. We present some general questions and issues that should be considered prior to embarking on molecular dynamics simulation studies of membrane proteins and review common simulation methods. We suggest a practical approach to setting up and running simulations of membrane proteins, and introduce two new (related) methods to embed a protein in a lipid bilayer. Both methods rely on placing lipids and the protein(s) on a widely spaced grid and then 'shrinking' the grid until the bilayer with the protein has the desired density, with lipids neatly packed around the protein. When starting from a grid based on a single lipid structure, or several potentially different lipid structures (method 1), the bilayer will start well-packed but requires more equilibration. When starting from a pre-equilibrated bilayer, either pure or mixed, most of the structure of the bilayer stays intact, reducing equilibration time (method 2). The main advantages of these methods are that they minimize equilibration time and can be almost completely automated, nearly eliminating one time consuming step in MD simulations of membrane proteins.     

Molecular modeling of the von Willebrand factor A2 Domain and the effects of associated type 2A von Willebrand disease mutations

A homology model for the A2 domain of von Willebrand factor (VWF) is presented. A large number of target–template alignments were combined into a consensus alignment and used for constructing the model from the structures of six template proteins. Molecular dynamics simulation was used to study the structural and dynamic effects of eight mutations introduced into the model, all associated with type 2A von Willebrand disease. It was found that the group I mutations G1505R, L1540P and S1506L cause significant deviations over multiple regions of the protein, coupled to significant thermal fluctuations for G1505R and L1540P. This suggests that protein instability may be responsible for their intracellular retention. The group II mutations R1597W, E1638K and G1505E caused single loop displacements near the physiologic VWF proteolysis site between Y1605–M1606. These modest structural changes may affect interactions between VWF and the ADAMTS13 protease. The group II mutations I1628T and L1503Q caused no significant structural change in the protein, suggesting that inclusion of the protease in this model is necessary for understanding their effect.Figure Homology model of the von Willebrand factor A2 domainElectronic Supplementary Material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s00894-004-0194-9     

Mn-SOD与抗癌胚抗原单链抗体基因的融合及高效表达

将Ｍｎ－ＳＯＤ与抗癌胚抗原（ＣＥＡ）单链抗体基因（Ｓｃ－Ｆｖ ｇｅｎｅ）融合,重组到含Ｔ７启动子的表达载体ｐＥＴ－２２ｂ（＋）中,构建表达质粒ｐＥＴＭｎ－ＳＯＤ－ＳｃＦｖ,并转化大肠杆菌ＢＬ２１（ＤＥ３）,进行高效表达,表达物占菌体可溶性总蛋白的２４％。ＳＤＳ－ＰＡＧＥ和蛋白质和迹图谱显示表达物分子量为４５ｋＤ与融合基因编码蛋白质的理论值相符。该蛋白质在大肠杆菌中为泌型表达有利于纯化。ＲＩＡ测定表     

Electron Microscopic and Physico-Chemical Studies of DNA Complexes with Synthetic Oligopeptides: Binding Specificity and DNA Compact Structures

Abstract

Binding to DNA of two synthetic peptides, Val-Thr-Thr-Val-Val-NH-NH-Dns and Thr-Val- Thr-Lys-Val-Gly-Thr-Lsy-Val-Gly-Thr-Val-Val-NH-NH-Dns (where Dns is a residue of 5- dimethylaminonaphthalene-l-sulfonic acid), has been studied by circular dichroism, electron microscopy and fluorescence methods. It has been found that these two peptides can self- associate in aqueous solution as follows from the fact that concentration-dependent changes are observed in the UV absorbance and fluorescence spectra. The two peptides can bind to DNA both in self-associated and monomeric forms. The pentapeptide in the β-associated form binds more strongly to poly(dG) · poly(dC) than to poly[d(A-C)] · poly[d(G-T)] and poly(dA) · poly(dT) whereas the tridecapeptide exhibits an opposite order of preferences binding more strongly to poly[d(A-C)] · poly[d(G-T)] and poly(dA) · poly(dT) than to poly(dG) · poly(dC).

Binding is a cooperative process which is accompanied by the DNA compaction at peptide/DNA base pair ratios greater than l. At the initial stage of the compaction process, the coalescence of DNA segments covered by bound peptide molecules leads to the formation of DNA loops stabilized by the interaction between peptide molecules bound to different DNA segments. Further increase in the peptide/DNA ratio leads to the formation of rod-like structures each consisting of two or more double-stranded DNA segments. The final stage of the compaction process involves folding of fibrillar macromolecular complexes into a globular structure containing only one DNA molecule.     

Molecular Dynamics Study of the Structure,Flexibility and Dynamics of Thermostable L1 Lipase at High Temperatures

Molecular Dynamics (MD) simulations have been used to understand how protein structure, dynamics, and flexibility are affected by adaptation to high temperature for several years. We report here the results of the high temperature MD simulations of Bacillus stearothermophilus L1 (L1 lipase). We found that the N-terminal moiety of the enzyme showed a high flexibility and dynamics during high temperature simulations which preceded and followed by clear structural changes in two specific regions; the small domain and the main catalytic domain or core domain of the enzyme. These two domains interact with each other through a Zn²⁺-binding coordination with Asp-61 and Asp-238 from the core domain and His-81 and His-87 from the small domain. Interestingly, the His-81 and His-87 were among the highly fluctuated and mobile residues at high temperatures. The results appear to suggest that tight interactions of Zn²⁺-binding coordination with specified residues became weak at high temperature which suggests the contribution of this region to the thermostability of the enzyme. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Localisation and dynamics of sodium counterions around DNA in solution from molecular dynamics simulation

The localisation and dynamics of sodium counterions around the DNA duplex d(AGCGTACTAGTACGCT)₂ corresponding to the trp operator fragment used in the crystal structure of the half site complex has been studied by a 1.4 ns molecular dynamics simulation in explicit solvent. A continuous and well-defined counterion density is shown to be present around the minor groove, while density patches are found in the major groove in regions where DNA bending is observed. A residence time analysis reveals the dynamic nature of these distributions. The resulting picture agrees with previous theoretical and experimental studies of A-tract DNA sequences, and is consistent with the polyelectrolyte condensation model. Received: 12 August 1999 / Revised version: 11 November 1999 / Accepted: 7 December 1999     

Discovery of a potent and selective adenylyl cyclase type 8 agonist by docking-based virtual screening

Adenylyl cyclases (ACs), which are responsible for catalyzing the conversion of adenosine triphosphate (ATP) into the second messenger cyclic adenosine monophosphate (cAMP), play a critical role in cell signal transduction. In this study, a combined approach involving docking-based virtual screening, with the combination of homology modeling followed by an in-vitro, and cell-based biological assay have been performed for discovering a class of novel potent and selective isoform adenylyl cyclase type 8 (AC8) agonist. The computer-aided virtual screening was used to identify fourteen virtual cluster compounds as potential hits which were further subjected to rigorous bioassays. A novel hit compound VHC-7 (ethyl 3-(2,4-dichlorobenzyl)-2-oxoindoline-3-carboxylate) was identified as a highly potent selective AC8 agonist with EC₅₀ value of 0.1052 ± 0.038 µM. Remarkably, the molecule herein reported can be explored further to discover greater number of hit compounds with better pharmacokinetic properties as well as to serve as a promising novel hit agonist of AC8 for the treatment of various central nervous system disorders and its associated diseases.     

猴血SV40抗体和SV40 DNA携带关系分析

目的　检测 5 7份恒河猴血清及对应的猴血中抗SV4 0抗体、SV4 0DNA的携带情况 ,找出抗体滴度与DNA携带的相关关系。方法　用间接免疫荧光法检测猴血中SV4 0中和抗体 ,聚合酶链反应 (PCR)法检测SV4 0st抗原DNA的携带情况。结果　 5 7份猴血清中 ,SV4 0抗体阳性率为 93% (5 3 5 7) ,抗体滴度最高为 1∶12 80 ,最低为 1∶10 ,GMT为 15 9.5 6 ,PCR检测猴血淋巴细胞SV4 0DNA阳性率为 2 4 .5 6 % (14 5 7) ,当抗体为 1 80以下时 ,病毒整合于血细胞 ,1∶80以上时 ,对st抗原基因的整合有抑制作用。结论　恒河猴SV4 0DNA的携带情况与抗体阳性和滴度呈反相关     

Expression of anti human IL-4 and IL-6 scFvs in transgenic tobacco plants

The two murine single-chain Fv (scFv) genes against human interleukin IL-4 and IL-6 cytokines were cloned in a plant expression vector (pGEJAE1) and mobilized to Agrobacterium tumefaciens. Tobacco leaf discs were co-cultured with Agrobacterium and transferred to selective media for regeneration. The tobacco in vitro plants produced scFvs against human IL-4 and IL-6. Only 8% of transformed plants expressing anti-IL-4 scFv were obtained versus 76% of transformed plants expressing anti-IL-6 scFv. In addition, some plants producing anti-IL-4 and anti-IL-6 scFvs aged more rapidly in in vitro conditions and in greenhouse pots than did control plants. Western blot analysis showed that the transformed Nicotiana tabacum plants contained proteins with an apparent molecular mass on electrophoresis of ca. 32 kDa, corresponding to the predicted size of the scFvs. As entire plant root seemed to accumulate more scFv than did leaves, we decided to continue working with isolated roots. Anti-IL-6 scFvs were detected in cultivated roots and their culture media. Functional anti-IL-6 scFv accounted for 0.16–0.18% of total soluble proteins. The affinity of the anti-IL-6 scFv produced in plants and measured by Biacore was similar to that of scFv produced in Escherichia coli. The high levels of antibody accumulation in isolated roots and secretion into the medium demonstrate the potential for producing recombinant protein in bioreactor systems.these authors contributed equally to this workthese authors contributed equally to this work     

Alternation of DNA and Solvent Layers in the A Form of d(GGCGCC) Obtainhted by Ethanol Crystallization

Abstract

We have determined by X-ray crystallography the structure of the hexamer duplex d(GGCGCC)₂ in the A-form using ethanol as a precipitant. The same sequence had previously been crystallized in the B-form, but with 2-methyl-2, 4-pentanediol as a precipitant. It appears that ethanol precipitation is a useful method to induce the formation of A-form crystals of DNA. Packing of the molecules in the crystal has unique features: the known interaction of A-DNA duplexes between terminal base-pairs and the minor groove of neighbor molecules is combined with a superstructure consisting in an alternation of DNA layers and solvent layers (water/ions). This organization in layers has been observed before, also with hexamers in the A conformation which crystallize in the same space group (C222 ₁). The solvent layer has a precise thickness, although very few ordered water molecules can be detected. Another feature of this crystal is its large unit cell, which gives rise to an asymmetric unit with three hexamer duplexes. One of the three duplexes is quite different from the other two in several aspects: the number of base pairs per turn, the twist pattern, the mean value of the twist angle and the fact that one terminal base-pair is not stacked as part of the duplex and appears to be disordered. So the variability in conformation of this sequence is remarkable.     

The structure of ColE1 rop in solution

Summary The structure of the ColE1 repressor of primer (rop) protein in solution was determined from the proton nuclear magnetic resonance data by a combined use of distance geometry and restrained molecular dynamics calculations. A set of structures was determined with low internal energy and virtually no violations of the experimental distance restraints. Rop forms homodimers: Two helical hairpins are arranged as an antiparallel four helix bundle with a left-handed rope-like twist of the helix axes and with left-handed bundle topology. The very compact packing of the side chains in the helix interfaces of the rop coiled-coil structure may well account for its high stability. Overall, the solution structure is highly similar to the recently determined X-ray structure (Banner, D.W., Kokkinidis, M. and Tsernoglou, D. (1987)J. Mol. Biol.,196, 657–675), although there are minor differences in regions where packing forces appear to influence the crystal structure.Abbreviations rop repressor of primer - NMR nuclear magnetic resonance - NOE nuclear Overhauser enhancement - NOESY NOE spectroscopy - RAN Set Structures generated from random choice of the dihedrai angles - HEL Set Structures generated from random choice of the dihedral angles restricted to ranges allowed for helices - MD molecular dynamics - EM energy minimization - RMSD root-mean-square deviation of atomic positions     

Knowledge-based Homology Modeling and Experimental Determination of Amino Acid Side Chain Accessibility by the Laser Photo CIDNP (Chemically Induced Dynamic Nuclear Polarization) Approach in Solution: Lessons from the Small Sialidase of Clostridium perfringens

The success of knowledge-based homology modelling is critically dependent on the predictive potency of the program structure-based calculations, which attempt to translate homologous sequences into three-dimensional structures, and on the actual relevance of the crystal structure for the protein topology. As quality control, experimental data for selected parameters of the proteins conformation are required. Using the crystal structure of the sialidase of Salmonella typhimurium as framework for model building of the homologous enzyme from Clostridium perfringens, a set of energy-minimised conformers is derived. These proteins present e.g. Tyr, Trp and His residues with an assessable area on the surface, since the side chains of these amino acid residues are responsive to chemically induced dynamic nuclear polarization (CIDNP), monitored by NMR. Hence, as first lesson, a comparative analysis for model-derived and experimentally determined values can be performed. The second lesson of this study concerns the notable impact of single amino acid substitutions (Tyr/Phe, Cys/Ser) on the surface accessibility of the CIDNP-reactive amino acid side chains in mutant forms of the sialidase. Corroborating the predictions from the theoretical calculations, the spectra of the engineered mutants reveal marked and non-uniform alterations. Thus, the effect of apparently rather conservative amino acid substitutions on a distinct conformational aspect of this protein, even at distant sites, should not be underestimated.