Fusion Peptide Interaction with Lipid Bilayer: Modeling with Monte Carlo Simulation and Continuum Electrostatics Calculation

Abstract

Conformation of 20-residue peptide E5, an analog of the fusion peptide of influenza virus hemagglutinin, was explored by Monte-Carlo technique starting with the fully buried in the membrane ideal α-helix. The lipid bilayer (of 30 Å width) together with surrounding water were modeled by the atomic solvation parameters. During the simulation, residues 2–18 of the peptide retained α-helical conformation, and the peptide was found to be partially immersed into the bilayer. In the resulting low-energy conformers, the N-terminus was buried inside the membrane, its position with respect to the bilayer surface (Z_NT) being varied from 2.5 to 7.5 Å, and the orientation of the helical axis relative to the membrane plane (Θ) – from 10 to 35°. The low-energy conformers (below -200kcal/mol) were clustered in the space (Z_NT, Θ) into 4 groups. To select low-energy states of the peptide compatible with NMR data, we calculated pKa values of E5 ionizable groups and compared them with the experimental values. It was shown that the best correlation coefficient (0.87) and rmsd (0.68 in pH units) were obtained for the group of states which is characterized by Θ = 15–19° and Z_NT = 3.5–4.5Å.     

Zebra: a web server for bioinformatic analysis of diverse protein families

During evolution of proteins from a common ancestor, one functional property can be preserved while others can vary leading to functional diversity. A systematic study of the corresponding adaptive mutations provides a key to one of the most challenging problems of modern structural biology – understanding the impact of amino acid substitutions on protein function. The subfamily-specific positions (SSPs) are conserved within functional subfamilies but are different between them and, therefore, seem to be responsible for functional diversity in protein superfamilies. Consequently, a corresponding method to perform the bioinformatic analysis of sequence and structural data has to be implemented in the common laboratory practice to study the structure–function relationship in proteins and develop novel protein engineering strategies. This paper describes Zebra web server – a powerful remote platform that implements a novel bioinformatic analysis algorithm to study diverse protein families. It is the first application that provides specificity determinants at different levels of functional classification, therefore addressing complex functional diversity of large superfamilies. Statistical analysis is implemented to automatically select a set of highly significant SSPs to be used as hotspots for directed evolution or rational design experiments and analyzed studying the structure–function relationship. Zebra results are provided in two ways – (1) as a single all-in-one parsable text file and (2) as PyMol sessions with structural representation of SSPs. Zebra web server is available at http://biokinet.belozersky.msu.ru/zebra.     

146 Amyloid fibril polymorphism probed by advanced vibrational spectroscopy

Amyloid fibrils are associated with many neurodegenerative diseases. All known amyloids including pathogenic and nonpathogenic forms display functional and structural heterogeneity (polymorphism) which determines the level of their toxicity. Despite a significant biological and biomedical importance, the nature of the amyloid fibril polymorphism remains elusive. We utilized for the first time three most advanced vibrational techniques to probe the core, the surface, and supramolecular chirality of fibril polymorphs. A new type of folding, aggregation phenomenon, spontaneous refolding from one polymorph to another, was discovered (Kurouski, Lauro et al., 2010). Hydrogen–deuterium exchange deep UV resonance Raman spectroscopy (Oladepo, Xiong et al., 2012) combined with advanced statistical analysis (Shashilov & Lednev, 2010) allowed for structural characterization of the highly ordered cross-β core of amyloid fibrils. We reported several examples showing significant variations in the core structure for fibril polymorphs. Amyloid fibrils are generally composed of several protofibrils and may adopt variable morphologies, such as twisted ribbons or flat-like sheets. We discovered the existence of another level of amyloid polymorphism, namely, that associated with fibril supramolecular chirality. Two chiral polymorphs of insulin, which can be controllably grown by means of small pH variations, exhibit opposite signs of vibrational circular dichroism (VCD) spectra (Kurouski, Dukor et al. 2012). VCD supramolecular chirality is correlated not only by the apparent fibril handedness but also by the sense of supramolecular chirality from a deeper level of chiral organization at the protofilament level of fibril structure. A small pH change initiates spontaneous transformation of insulin fibrils from one polymorph to another. As a result, fibril supramolecular chirality overturns both accompanying morphological and structural changes (Kurouski, Dukor et al. 2012). No conventional methods could probe the fibril surface despite its significant role in the biological activity. We utilized tip-enhanced Raman spectroscopy (TERS) to characterize the surface structure of an individual fibril due to a high depth and lateral spatial resolution of the method in the nanometer range (Kurouski, Deckert-Gaudig et al. 2012). It was found that the surface is strongly heterogeneous and consists of clusters with various protein conformations and amino acid composition.     

Ethidium Probing of the Parallel Double- and Four-stranded Structures Formed by the Telomeric DNA Sequences dG(GT)4G and d(GT)5

Abstract

Oligonucleotides 3′-d(GT)₅-(CH₂CH₂O)₃-d(GT)₅-3′ (parGT), containing GT repeats present in the telomeric DNA from Saccharomyces cerevisiae, had been demonstrated to form bimolecular structure, GT-quadruplex (qGT) [O. F. Borisova et al. FEBS Letters 306, 140–142 (1992)]. Four d(GT)₅ strands of the GT-quadruplex are parallel and form five G-quartets while thymines are bulged out. The four GT repeats when flanked by guanines, 3′-dG(TG)₄G-(CH₂CH₂O)₃- dG(GT)₄G-3′ (hp-GT), had been shown to form a novel parallel-stranded (ps) double helix with G·G and T·T base pairs (hp-GT ps-DNA) [A. K. Shchyolkina et al. J. Biomol. Struct. Dynam. 18, 493–503 (2001)]. In the present study the intercalator ethidium bromide (Et) was used for probing the two structures. The mode of Et binding and its effect on thermostability of qGT and hp-GT were compared. The quantum yield (q) and the fluorescence lifetime (τ) of Et:qGT (q = 0.15 ±0.01 and τ = 24 ±1 ns) and Et:hp-GT (q = 0.10 ± 0.01 and τ = 16.5 ± 1 ns) indicative of intercalation mode of Et binding were determined. Et binding to qGT was found to be cooperative with corresponding coefficient ω = 3.9 ± 0.1 and the binding constant K= (6.4 ± 0.1)·10M^?1. The maximum number of Et molecules intercalating into GT-quadruplex is as high as twice the number of inerspaces between G-quartets (eight in our case). The data conform to the model of Et association with GT-quadruplex suggested earlier [O. F. Borisova et al. Mol. Biol. (Russ) 35, 732–739 (2001)]. The anticooperative type of Et binding was observed in case of hp- GT ps-DNA, with the maximum number of bound Et molecules, N = 4 ÷ 5, and the association constant K = (1.5 ± 0.1)·10⁵ M^?1. Thermodynamic parameters of formation of Et:qGT and EtBr:hp-GT complexes were calculated from UV thermal denaturation profiles.     

Melting of Cross-Linked DNA IV. Methods for Computer Modeling of Total Influence on DNA Melting of Monofunctional Adducts,Intrastrand and Interstrand Cross-Links Formed by Molecules of an Antitumor Drug

Abstract

A theoretical method is developed for calculation of melting curves of covalent complexes of DNA with antitumor drugs. The method takes into account all the types of chemical modifications of the double helix caused by platinum compounds and DNA alkylating agents: 1) monofunctional adducts bound to one nucleotide; 2) intrastrand cross-links which appear due to bidentate binding of a drug molecule to two nucleotides that are included into the same DNA strand; 3) interstrand cross-links caused by bidentate binding of a molecule to two nucleotides of different strands. The developed calculation method takes into account the following double helix alterations at sites of chemical modifications: 1) a change in stability of chemically modified base pairs and neighboring ones, that is caused by all the types of chemical modifications; 2) a change in the energy of boundaries between helical and melted regions at sites of chemical modification (local alteration of the factor of cooperativity of DNA melting), that is caused by all the types of chemical modifications, too; 3) a change in the loop entropy factor of melted regions that include interstrand cross-links; 4) the prohibition of divergence of DNA strands in completely melted DNA molecules, which is caused by interstrand cross-links only. General equations are derived, and three calculation methods are proposed to calculate DNA melting curves and the parameters that characterize the helix-coil transition.     

Nucleosome-like Complex of the Histone from the Hyperthermophile Methanopyrus Kandleri (MkaH) with Linear DNA

Abstract

The MkaH protein from the archaeon Methanopyrus kandleri, an unusual assembly of two histone-fold domains in a single polypeptide chain, demonstrates high structural similarity to eukaryal histones. We studied the DNA binding and self-association properties of MkaH by means of the electrophoretic mobility shift assay (EMSA), electron microscopy (EM), chemical cross-linking, and analytical gel filtration. EMSA showed an increased mobility of linear DNA complexed with MkaH protein with a maximum at a protein-DNA weight ratio (Rw) of ≈3; the mobility decreased at higher protein concentration. EM of the complexes formed at Rw ≤ 3 revealed formation of isometric loops encompassing 71 +/- 7 bp of DNA duplex. At high values of Rw (≥9) thickened compact nucleoprotein structures were observed; no individual loops were seen within the complexes. Gel filtration chromatography and chemical fixation indicated that in the absence of DNA the dominant form of the MkaH in solution, unlike other archaeal histones, is a stable dimer (pseudo-tetramer of the histone-fold domain) apparently resembling the eukaryal (H3-H4)₂ tetramer. Similarly, dimers are the dominant form of the protein interacting with DNA. The properties of MkaH supporting the assignment of its intermediate position between other archaeal and eukaryal histones are discussed.     

Melting of Cross-Linked DNA V. Cross-Linking Effect Caused by Local Stabilization of the Double Helix

Abstract

DNA interstrand cross-links are usually formed due to bidentate covalent or coordination binding of a cross-linking agent to nucleotides of different strands. However interstrand linkages can be also caused by any type of chemical modification that gives rise to a strong local stabilization of the double helix. These stabilized sites conserve their helical structure and prevent local and total strand separation at temperatures above the melting of ordinary AT and GC base pairs. This local stabilization makes DNA melting fully reversible and independent of strand concentration like ordinary covalent interstrand cross-links. The stabilization can be caused by all the types of chemical modifications (interstrand cross-links, intrastrand cross-links or monofunctional adducts) if they give rise to a strong enough local stabilization of the double helix. Our calculation demonstrates that an increase in stability by 25 to 30 kcal in the free energy of a single base pair of the double helix is sufficient for this “cross-linking effect” (i.e. conserving the helicity of this base pair and preventing strand separation after melting of ordinary base pairs). For the situation where there is more then one stabilized site in a DNA duplex (e.g., 1 stabilized site per 1000 bp), a lower stabilization per site is sufficient for the “cross-linking effect” (18–20 kcal). A substantial increase in DNA stability was found in various experimental studies for some metal-based anti-tumor compounds. These compounds may give rise to the effect described above. If ligand induced stabilization is distributed among several neighboring base pairs, a much lower minimum increase per stabilized base pair is sufficient to produce the cross-linking effect (1 bp- 24.4 kcal; 5 bp- 5.3 kcal; 10 bp- 2.9 kcal, 25 bp- 1.4 kcal; 50 bp- 1.0 kcal). The relatively weak non-covalent binding of histones or protamines that cover long regions of DNA (20–40 bp) can also cause this effect if the salt concentration of the solution is sufficiently low to cause strong local stabilization of the double helix. Stretches of GC pairs more than 25 bp in length inserted into poly(AT) DNA also exhibit properties of stabilizing interstrand cross-links.     

Structural basis for cytokinin receptor signaling: an evolutionary approach

Cytokinins are ubiquitous plant hormones; their signal is perceived by sensor histidine kinases—cytokinin receptors. This review focuses on recent advances on cytokinin receptor structure, in particular sensing module and adjacent domains which play an important role in hormone recognition, signal transduction and receptor subcellular localization. Principles of cytokinin binding site organization and point mutations affecting signaling are discussed. To date, more than 100 putative cytokinin receptor genes from different plant species were revealed due to the total genome sequencing. This allowed us to employ an evolutionary and bioinformatics approaches to clarify some new aspects of receptor structure and function. Non-transmembrane areas adjacent to the ligand-binding CHASE domain were characterized in detail and new conserved protein motifs were recovered. Putative mechanisms for cytokinin-triggered receptor activation were suggested.     

Morphology of the jaw apparatus in 8 species of Patellogastropoda (Mollusca, Gastropoda) with special reference to Testudinalia tesulata (Lottiidae)

The fine structure of the jaw apparatus was studied by scanning electron microscopy in eight species of Patellogastropoda. The jaw apparatus is an unpaired two-layered dorsolateral structure with anterior and posterior wings attached to the odontophore by muscles. The jaw of Testudinalia tesulata (O.F. Müller, 1776) is a derivative of the cuticle typical for the foregut. The tissue forming the jaw is a specialized foregut epithelium (gnathoepithelium), consisting of a special type of cells called gnathoblasts. The jaw grows in areas of the epithelium characterized by high concentration of electron-dense vesicles, ER and long microvilli that penetrate deep into the jaw plate. This indicates that the gnathoblasts take an active part in jaw growth. In most cases, these areas of the gnathoepithelium are highly folded. The main differences between the species studied are form and thickness of the frontal edge of the jaw. These differences do not correlate with the systematic position of the species studied but likely depend more on the feeding mode. The transmission electron microscopy studies yielded new morphological criteria for comparison between various gastropod species and other members of Trochozoa, in particular, Annelida. The jaws of Annelida are cuticular structures formed on the surface of specialized epithelial cells, often also called gnathoblasts. The jaw of Patellogastropoda can be attributed to the first type of annelid jaw formation characterized by an epithelium with long microvilli and continuous growth.     

Modeling of Peptides in Implicit Membrane-Mimetic Media

Abstract

Lipid bilayer plays a crucial role in folding of membrane peptides and their stabilization in the membrane-bound state. Correct treatment of the media effects is thus essential for realistic simulations of peptides in bilayers. Previously (Volynsky et al., 1999), we proposed an efficient solvation model which mimics heterogeneous membrane-water system. The model is based on combined employment of atomic solvation parameters for water and hydrocarbon, which approximate hydrated headgroups and acyl chains of lipids, respectively. In this study, the model is employed in non-restrained Monte Carlo simulations of several peptides: totally apolar 20-residue poly-L-Leu, hydrophobic peptide with polar edges, and strongly amphiphilic pep-tide. The principal goals are: to explore energy landscape of these peptides in membrane; to characterize the structures of low-energy states and their orientations with respect to the bilayer. Simulations were performed starting from different structures (unordered or helical) and orientations. It was found that the membrane environment significantly promotes an α-helical conformation for all the peptides, while their energetically favourable orientations are quite different. Thus, poly-Leu was immobilized inside the membrane, the hydrophobic peptide with polar termini adapted transbilayer orientation, whereas the amphiphilic peptide stayed on the lipid-water interface in peripherial orientation. Energy barriers between different states were characterized. The computational results were compared with the experimental structural data.