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.     

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.     

Coagulation Factor Binding Orientation and Dimerization May Influence Infectivity of Adenovirus-Coagulation Factor Complexes

Adenoviruses (Ads) are promising vectors for therapeutic interventions in humans. When injected into the bloodstream, Ad vectors can bind several vitamin K-dependent blood coagulation factors, which contributes to virus sequestration in the liver by facilitating transduction of hepatocytes. Although both coagulation factors FVII and FX bind the hexon protein of human Ad serotype 5 (HAdv5) with a very high affinity, only FX appears to play a role in mediating Ad-hepatocyte transduction in vivo. To understand the discrepancy between efficacy of FVII binding to hexon and its apparently poor capacity for supporting virus cell entry, we analyzed the HAdv5-FVII complex by using high-resolution cryo-electron microscopy (cryo-EM) followed by molecular dynamic flexible fitting (MDFF) simulations. The results indicate that although hexon amino acids T423, E424, and T425, identified earlier as critical for FX binding, are also involved in mediating binding of FVII, the FVII GLA domain sits within the surface-exposed hexon trimer depression in a different orientation from that found for FX. Furthermore, we found that when bound to hexon, two proximal FVII molecules interact via their serine protease (SP) domains and bury potential heparan sulfate proteoglycan (HSPG) receptor binding residues within the dimer interface. In contrast, earlier cryo-EM studies of the Ad-FX interaction showed no evidence of dimer formation. Dimerization of FVII bound to Ad may be a contributing mechanistic factor for the differential infectivity of Ad-FX and Ad-FVII complexes, despite high-affinity binding of both these coagulation factors to the virus.     

Discovery of subtype selective muscarinic receptor antagonists as alternatives to atropine using in silico pharmacophore modeling and virtual screening methods

Muscarinic acetylcholine receptors (mAChRs) have five known subtypes which are widely distributed in both the peripheral and central nervous system for regulation of a variety of cholinergic functions. Atropine is a well known muscarinic subtype non-specific antagonist that competitively inhibits acetylcholine (ACh) at postganglionic muscarinic sites. Atropine is used to treat organophosphate (OP) poisoning and resulting seizures in the warfighter because it competitively inhibits acetylcholine (ACh) at the muscarinic cholinergic receptors. ACh accumulates due to OP inhibition of acetylcholinesterase (AChE), the enzyme that hydrolyzes ACh. However, atropine produces several unwanted side-effects including dilated pupils, blurred vision, light sensitivity, and dry mouth. To overcome these side-effects, our goal was to find an alternative to atropine that emphasizes M1 (seizure prevention) antagonism but has minimum M2 (cardiac) and M3 (e.g., eye) antagonism so that an effective less toxic medical countermeasure may be developed to protect the warfighter against OP and other chemical warfare agents (CWAs). We adopted an in silico pharmacophore modeling strategy to develop features that are characteristics of known M1 subtype-selective compounds and used the model to identify several antagonists by screening an in-house (WRAIR-CIS) compound database. The generated model for the M1 selectivity was found to contain two hydrogen bond acceptors, one aliphatic hydrophobic, and one ring aromatic feature distributed in a 3D space. From an initial identification of about five hundred compounds, 173 compounds were selected through principal component and cluster analyses and in silico ADME/Toxicity evaluations. Next, these selected compounds were evaluated in a subtype-selective in vitro radioligand binding assay. Twenty eight of the compounds showed antimuscarinic activity. Nine compounds showed specificity for M1 receptors and low specificity for M3 receptors. The pK_i values of the compounds range from 4.5 to 8.5 nM in comparison to a value of 8.7 nM for atropine. 2-(diethylamino)ethyl 2,2-diphenylpropanoate (ZW62841) was found have the best desired selectivity. None of the newly found compounds were previously reported to exhibit antimuscarinic specificity. Both theoretical and experimental results are presented.