Molecular dynamics simulation of interaction between phospholipid membrane and pyrazine and its derivatives

We have performed a comparative molecular dynamics simulation of the diffusion process of the heterocyclic compound pyrazine and its methylated derivatives into the model membrane phospholipid bilayer. Several structural and dynamical bilayer parameters were measured, and qualitative interrelations between parameter changes and the substituted pyrazine structure were studied. The simulation results support the hypothesis that molecular mechanisms of biological effects of substituted pyrazines involve dissolution of the effector molecule in the membrane bilayer and subsequent changes in bilayer properties. This stage can provide the means for pyrazine molecules to interact with integral membrane proteins, directly or indirectly through the changed lipid environment of the protein.     

Ubericin A, a class IIa bacteriocin produced by Streptococcus uberis

Streptococcus uberis, a causal agent of bovine mastitis, produces ubericin A, a 5.3-kDa class IIa (pediocin-like) bacteriocin, which was purified and characterized. The uba locus comprises two overlapping genes: ubaA (ubericin A precursor peptide) and ubaI (putative immunity protein). Ubericin A is the first streptococcal class IIa bacteriocin to be characterized.     

Molecular dynamics simulation study of tubulin dimer interaction with cytostatics

Colchicine, podophylotoxin, and indibulin are natural cytostatics that are used in the treatment of neoplasms. However, application of the compounds is restricted due to their high toxicity and low specificity. Computational experiments modeling tubulin interactions with the cytostatics seem a promising approach to design new analogues of the above-mentioned drugs with higher cytostatic activity and lower toxicity. Therefore, the CHARMM software was used to examine the macromolecules using molecular dynamics and mechanics methods. Particularly, a procedure was applied according to which molecules of each studied cytostatics were placed at several various random positions around the predicted binding site on tubulin. As a result, cytostatic binding regions were identified on the tubulin molecule. It was shown that, during the interaction, structural alterations occurred in these regions that may be responsible for tubulin polymerization. Thus, alterations have been revealed for the first time in the structure of tubulin in the regions of cytostatic binding that can substantially affect its function.     

Development of a PCR-based assay for rapid detection of class IIa bacteriocin genes

Aims: We have developed a PCR‐based assay using custom designed panel of primers which allows rapid detection of class IIa bacteriocin‐coding genes. To demonstrate the applicability of the developed assay, the method was applied on 40 metagenomic DNA preparations isolated from native microbiota of Polish artisanal cheeses produced in the Tatra Mountains. Methods and Results: The developed assay was designed on the basis of a large scale alignment of class IIa bacteriocin‐coding genes. A panel of seven primer pairs with confirmed ability to detect class IIa bacteriocin‐coding sequences was obtained. The following study has revealed a superb bacteriocinogenic potential of all forty analysed cheese samples. Conclusions: The majority of obtained sequences were lactic acid bacteria (LAB) related, although some sequences showed significant similarity to bacteriocin‐coding sequences present in non‐LAB bacteriocin producers. The results suggest that several potentially new bacteriocin‐coding sequences were found. Significance and Impact of the Study: The developed assay can be extremely helpful in establishing whether isolates from the environment of interest have a potential of synthesizing antilisterial class IIa bacteriocins. Application of the approach may represent a useful tool contributing to ecological studies looking for valuable probiotic, bacteriocinogenic microbiota developing in foods.     

NMR study on the interaction between MHC class I protein and its antigen peptide

A major histcompatibility complex (MHC) class I protein H-2K(b) was expressed in a large scale as a fusion protein with thioredoxin and hexahistidine at the N-terminus to analyze the interaction with the antigen peptide SIYRYYGL. NMR spectra of the peptide in the mixture solution with the protein showed very broad signals, indicating the obviously clear existence of the dynamic interaction between the class I protein and the antigen peptide. The interaction of the protein and peptide was discussed as well as the surrounding atmosphere of the peptide in the complex.     

Molecular dynamics simulation of class 3 aldehyde dehydrogenase

Molecular dynamics (MD) simulation of the rat class 3 aldehyde dehydrogenase (ALDH) with nicotinamide dinucleotide (NAD) cofactors and explicit water molecules are reported. Our results demonstrate that MD simulation using the latest methodologies can maintain the crystal structure of the enzyme, as well as closely reproduce the short timescale dynamics of the enzyme. Furthermore, the examination of the distance between the nucleophilic Cys-243 and the NAD cofactor reveal important fluctuations that could be linked to ALDH catalysis. Finally, our quantum mechanical model of benzaldehyde in the active site of ALDH demonstrates that the enzyme requires only minor conformational changes to be poised for nucleophilic attack on the substrate.     

Action of divergicin M35, a class IIa bacteriocin, on liposomes and Listeria

AIMS: The mode of action of divergicin M35, a class IIa bacteriocin, was studied against Listeria monocytogenes with sensitive (DivS) and resistant (DivM) phenotypes, as well as on synthetic phospholipid liposomes. METHODS AND RESULTS: Divergicin-induced release of 1,6-diphenyl-1,3,5-hexatriene (DPH) from zwitterionic (DMPC) and anionic (DMPC/DMPG, 4:1) liposomes, divergicin binding to liposomes, intracellular ATP concentration, cation efflux, cell affinity for hydrocarbons and cell lysis were measured and cell damage was visualized by fluorescence imaging and transmission electron microscopy. Divergicin M35 at 5 microg ml(-1) induced DPH efflux from anionic and zwitterionic liposomes at rates of about 2.58% and 1.61% per minute, respectively. DPH efflux rate from anionic liposomes was reduced by about 1.83% and 2.1% per minute in the presence of Li+ and Ca2+, respectively. Binding affinity of divergicin M35 to anionic and zwitterionic liposomes was about 86% and 63%, respectively. Intracellular ATP decreased in the sensitive and the resistant strains by 96.7% and 72.8%, respectively after 20 min of exposure to 5 microg ml(-1) divergicin M35. Lysis of the sensitive strain reached 57% in 18 h at a concentration of 5 microg ml(-1) when compared with the lysis of the divergicin-resistant strain (38.8%). The K+ and Na+ efflux from the divergicin-sensitive strain reached 87% and 80% of the total ion content within 5 min of exposure. This strain also showed higher affinity for hydrocarbons. CONCLUSIONS: The cell death of listerial strains upon addition of divergicin M35 could result from ATP depletion, K+ and Na+ efflux, and bacteriolysis. This triple biological effect was attenuated in the DivM strain. SIGNIFICANCE AND IMPACT OF THE STUDY: This study contributed to the understanding of the mode of action of divergicin M35, a pediocin-like bacteriocin.     

Cloning, expression, and identification of a novel class IIa bacteriocin in the Escherichia coli cell-free protein expression system

The NB-C1 gene, acquired from the result of data mining of the lactic acid bacteria genome, is a novel potential class IIa bacteriocin gene with the characteristic YGNGVxC cluster. To produce soluble NB-C1 efficiently and overcome issues of protein toxicity, we adopted a GFP fusion strategy using an Escherichia coli cell-free protein expression system. We constructed the expression vector pIVEX2.4d-GFP-NB-C1, which was expressed in both the batch mode and the continuous exchange cell-free (CECF) systems. The amount of soluble fusion protein achieved from the CECF system (2.2 mg/ml) was approximately three times higher than that in the batch mode (0.73 mg/ml). The soluble fusion protein was purified via one-step Ni–NTA affinity chromatography, with a concentration of 0.26 mg/ml and a purity of 95%. The purified NB-C1 showed strong antimicrobial activity against the indicator bacteria Listeria monocytogenes.     

Enterocin T, a novel class IIa bacteriocin produced by Enterococcus sp. 812

Enterococcus sp. 812, isolated from fresh broccoli, was previously found to produce a bacteriocin active against a number of Gram-positive bacteria, including Listeria monocytogenes. Bacteriocin activity decreased slightly after autoclaving (121 °C for 15 min), but was inactivated by protease K. Mass spectrometry analysis revealed the bacteriocin mass to be approximately 4,521.34 Da. N-terminal amino acid sequencing yielded a partial sequence, NH₂-ATYYGNGVYXDKKKXWVEWGQA, by Edman degradation, which contained the consensus class IIa bacteriocin motif YGNGV in the N-terminal region. The obtained partial sequence showed high homology with some enterococcal bacteriocins; however, no identical peptide or protein was found. This peptide was therefore considered to be a novel bacteriocin produced by Enterococcus sp. 812 and was termed enterocin T.     

Genetic characterisation and heterologous expression of leucocin C, a class IIa bacteriocin from Leuconostoc carnosum 4010

Leuconostoc carnosum 4010 is a protective culture for meat products. It kills the foodborne pathogen Listeria monocytogenes by producing two class IIa (pediocin-like) bacteriocins, leucocin A and leucocin C. The genes for leucocin A production have previously been characterised from Leuconostoc gelidum UAL 187, whereas no genetic studies about leucocin C has been published. Here, we characterised the genes for the production of leucocins A and C in L. carnosum 4010. In this strain, leucocin A and leucocin C operons were localised in different plasmids. Unlike in L. gelidum, leucocin A operon in L. carnosum 4010 only contained the structural and the immunity genes lcaAB without transporter genes lcaECD. On the contrary, leucocin C cluster included two intact operons. Novel genes lecCI encode the leucocin C precursor and the 97-aa immunity protein LecI, respectively. LecI shares 48 % homology with the immunity proteins of sakacin P and listeriocin. Another leucocin C operon lecXTS, encoding an ABC transporter and an accessory protein, was 97 % identical with the leucocin A transporter operon lcaECD of L. gelidum. For heterologous expression of leucocin C in Lactococcus lactis, the mature part of the lecC gene was fused with the signal sequence of usp45 in the secretion vector pLEB690. L. lactis secreted leucocin C efficiently, as shown by large halos on lawns of L. monocytogenes and Leuconostoc mesenteroides indicators. The function of LecI was then demonstrated by expressing the gene lecI in L. monocytogenes. LecI-producing Listeria was less sensitive to leucocin C than the vector strain, thus corroborating the immunity function of LecI.     

NMR solution structure of ImB2, a protein conferring immunity to antimicrobial activity of the type IIa bacteriocin, carnobacteriocin B2

Bacteriocins produced by lactic acid bacteria are potent antimicrobial compounds which are active against closely related bacteria. Producer strains are protected against the effects of their cognate bacteriocins by immunity proteins that are located on the same genetic locus and are coexpressed with the gene encoding the bacteriocin. Several structures are available for class IIa bacteriocins; however, to date, no structures are available for the corresponding immunity proteins. We report here the NMR solution structure of the 111-amino acid immunity protein for carnobacteriocin B2 (ImB2). ImB2 folds into a globular domain in aqueous solution which contains an antiparallel four-helix bundle. Extensive packing by hydrophobic side chains in adjacent helices forms the core of the protein. The C-terminus, containing a fifth helix and an extended strand, is held against the four-helix bundle by hydrophobic interactions with helices 3 and 4. Most of the charged and polar residues in the protein face the solvent. Helix 3 is well-defined to residue 55, and a stretch of nascent helix followed by an unstructured loop joins it to helix 4. No interaction is observed between ImB2 and either carnobacteriocin B2 (CbnB2) or its precursor. Protection from the action of CbnB2 is only observed when ImB2 is expressed within the cell. The loop between helices 3 and 4, and a hydrophobic pocket which it partially masks, may be important for interaction with membrane receptors responsible for sensitivity to class IIa bacteriocins.     

Crystal structure and mutagenic analysis of a bacteriocin immunity protein,Mun-im

Bacteriocin-producing lactic acid bacteria (LAB) possess a self-protection factor, which is generally called an immunity protein. In this study, we determine the crystal structure of an immunity protein, designated Mun-im, which was classified into subgroup B immunity proteins for class IIa bacteriocins. The Mun-im protein takes a left-turning antiparallel four-helix bundle structure with the flexible N- and C-terminal parts. Although the amino acid sequences of the subgroup B immunity proteins are distinguished from those of the subgroup A, the core structure of Mun-im is well-superimposed with that of the subgroup A immunity protein, EntA-im, and the C-terminus of both proteins is flexible. However, the C-terminus of Mun-im is obviously shorter than that of the subgroup A. We found through mutagenic study of Mun-im that the C-terminus and the K86 residue on the helix 4 in the immunity protein molecule are important for expression of the immunity activity on the subgroup B immunity proteins.     

Molecular dynamics simulation of Pf1 coat protein.

The results of molecular dynamics simulations of Pf1 coat protein are described and compared to experimental NMR data on both the membrane bound and structural forms of this viral coat protein. Molecular dynamics simulations of the 46 residue coat protein and related model sequences were performed according to a simple protocol. The simulations were initiated with the polypeptides in a completely uniform alpha helical conformation in a dielectric continuum (epsilon = 2) and the motions of individual residues were followed as a function of time by monitoring the angular fluctuations of amide NH bond vectors. The simulations of Pf1 coat protein were able to identify the same mobile and structured segments found in experimental NMR studies of the membrane bound form of the protein (Shon, K.-J., Y. Kim, L. A. Colnago, and S. J. Opella. 1991. Science (Wash. DC). 252:1303-1305). Significantly, in addition to mobile amino and carboxyl terminal regions, a mobile internal loop was found that connects the rigid hydrophobic and amphipathic helices in the protein. NMR experiments show that this mobile loop is present in both the viral and membrane bound forms of the protein and that it plays a role in viral assembly (Nambudripad, R., W. Stark, S. J. Opella, and L. Makowski. 1991. Science (Wash. DC) 252:1305-1308). The results of simulations of several alanine based 46 residue polypeptides with some of the charged residues present in the Pf1 coat protein sequence suggest that interactions between the Asp 14 and Asp 18 sidechains and the peptide backbone are responsible for the formation of the mobile loop.(ABSTRACT TRUNCATED AT 250 WORDS)     

Molecular dynamics simulation and nmr study of a blood group H trisaccharide

Molecular dynamics simulations in vacuum and solution have been carried out on 2′-α-L -fucosyllactitol, a model for blood group H in conjunction with two-dimensional nmr measurements on the same compound. Three independent starting conformations for the dynamics were chosen from low energy conformations obtained by a ?/ψ grid search. Nine 5 ns vacuum simulations of the trisaccharide were performed, employing three different ways to treat electrostatic interactions for each starting conformation: distance-dependent dielectric with ε = r, constant dielectric with ε = 1, or constant dielectric with ε = 80. In vacuum, transitions of ? and ψ for the α-L -Fuc-(1 → 2)-β-D -Gal element occur in a cooperative manner. The virtual distance obtained for H1 in fucose to H2 in galactose from nuclear Overhauser effect spectroscopy experiments agree with one of the conformations of the trisaccharide in one of the three 100 ps aqueous simulations (?/ψ ca. ?100°/150°), indicating this may be a dominant solution conformation. The rms fluctuations of the ?- and ψ-dihedral angles were ～ 10° for a conformational state, both in the vacuum and the aqueous simulations. For the simulations in vacuum, the agreement with experimental NOE data is reasonable when a constant dielectric of 1 is used (major conformers having ?/ψ ca. ?100°/150° and ?140°/100°), whereas the agreement was poor with a constant dielectric of 80. Translational diffusion coefficients calculated from the simulation of the oligosaccharides were 0.12–0.18 × 10^?5 cm²/s and from nmr measurements 0.27 × 10^?5 cm²/s. © 1994 John Wiley & Sons, Inc.     

肿瘤坏死因子与其受体相互作用的计算机模拟研究

利用同源模建方法,以肿瘤坏死因子(TNF)R55受体胞外区的晶体结构为参考模板,预测了TNFR75受体胞外区Cys18-Phe147片段的三维结构。根据R55受体胞外区与LT相结合的复合物的晶体结构,预测了TNF与R55及R75胞外区的复合物的三维结构,模拟了TNF与受体之间的相互作用。由于TNF与受体的作用形式是三聚体对三聚体,因此在模拟TNF与受体相互作用时选择了包括一个非对称的TNF三聚体和一个受体(R55或R75)单体的模拟系统。结合已有的突变体实验结果,利用计算机分析手段,发现了一些TNF突变体之所以具有受体选择性的三维结构基础和发挥了关键作用的氨基酸残基以及这些残基之间的主要作用形式。研究深化了对已有的突变体实验结果的认识,建立了不同的实验结果之间的内在关联,为以后有目的的新型突变体设计和实验研究打下了基础。     

Molecular dynamics simulations of interaction between sub-bituminous coal and water

The high moisture content of sub-bituminous coal is associated with the interactions between coal and water. Because of complex composition and structure, the graphite surface modified by hydroxyl, carboxyl and carbonyl groups was used to represent the surface model of sub-bituminous coal according to XPS results. Density profiles for oxygen atoms and hydrogen atoms indicate that the coal surface properties affect the structural and dynamic characteristics of the interfacial water molecules. The interfacial water exhibits much more ordering than bulk water. The results of radial distribution functions, mean square displacement and local self-diffusion coefficient for water molecule related to three oxygen moieties confirmed that the water molecules prefer to absorb with carboxylic groups, and adsorption of water molecules at the hydroxy and carbonyl is similar.     

Molecular dynamics investigation of the interaction between DNA and distamycin

The complex of the minor groove binding drug distamycin and the B-DNA oligomer d-(CGCAAATTTGCG) was investigated by molecular dynamics simulations. For this purpose, accurate atomic partial charges of distamycin were determined by extended quantum chemical calculations. The complex was simulated without water but with hydrated counterions. The oligomer without the drug was simulated in the same fashion and also with 1713 water molecules and sodium counterions. The simulations revealed that the binding of distamycin in the minor groove induces a stiffening of the DNA helix. The drug also prevents a transition from B-DNA to A-DNA that was found to occur rapidly (30 ps) in the segment without bound distamycin in a water-free environment but not in simulations including water. In other simulations, we investigated the relaxation processes after distamycin was moved from its preferred binding site, either radially or along the minor groove. Binding in the major groove was simulated as well and resulted in a bound configuration with the guanidinium end of distamycin close to two phosphate groups. We suggest that, in an aqueous environment, tight hydration shells covering the DNA backbone prevent such an arrangement and thus lead to distamycin's propensity for minor groove binding.     

Molecular dynamics simulation of a phospholipid membrane

We present the results of molecular dynamics (MD) simulations of a phospholipid membrane in water, including full atomic detail. The goal of the simulations was twofold: first we wanted to set up a simulation system which is able to reproduce experimental results and can serve as a model membrane in future simulations. This goal being reached it is then further possible to gain insight in to those properties that are experimentally more difficult to access. The system studied is dipalmitoylphosphatidylcholine/water, consisting of 5408 atoms. Using original force field parameters the membrane turned out to approach a gel-like state. With slight changes of the parameters, the system adopted a liquid-crystalline state. Separate 80 ps runs were performed on both the gel and liquid-crystalline systems. Comparison of MD results with reliable experimental data (bilayer repeat distance, surface area per lipid, tail order parameters, atom distributions) showed that our simulations, especially the one in the liquid-crystalline phase, can serve as a realistic model for a phospholipid membrane. Further analysis of the trajectories revealed valuable information on various properties. In the liquid-crystalline phase, the interface turns out to be quite diffuse, with water molecules penetrating into the bilayer to the position of the carbonyl groups. The 10–90% width of the interface turns out to be 1.3 nm and the width of the hydrocarbon interior 3.0 nm. The headgroup dipoles are oriented at a small angle with respect to the bilayer plane. The resulting charge distribution is almost completely cancelled by the water molecules. The electron density distribution shows a large dip in the middle of the membrane. In this part the tails are more flexible. The mean life time between dihedral transitions is 20 ps. The average number of gauche angles per tail is 3.5. The occurrence of kinks is not a significant feature.Abbreviations MD molecular dynamics - DPPC dipalmitoylphosphatidylcholine - SPC simple point charges - DPPE dipalmitoylphosphatidylethanolamine Correspondence to: H. J. C. Berendsen     

Molecular dynamics simulation of a phosphatidylglycerol membrane

Although molecular dynamics simulations are an important tool for studying membrane systems, relatively few simulations have used anionic lipids. This paper reports the first simulation of a pure phosphatidylglycerol (PG) bilayer. The properties of this equilibrated palmitoyloleoylphosphatidylglycerol membrane agree with experimental observations of PG membranes and with previous simulations of monolayers and mixed bilayers containing PG lipids. These simulations also provide interesting insights into hydrogen bonding interactions in PG membranes. This equilibrated membrane will be a useful starting point for simulations of membrane proteins interacting with PG lipids.