Anton O. Chugunov

PLATINUM: a web tool for analysis of hydrophobic/hydrophilic organization of biomolecular complexes

UNLABELLED The PLATINUM (Protein-Ligand ATtractions Investigation NUMerically) web service is designed for analysis and visualization of hydrophobic/hydrophilic properties of biomolecules supplied as 3D-structures. Furthermore, PLATINUM provides a number of tools for quantitative characterization of the hydrophobic/hydrophilic match in biomolecular complexes e.g. in docking poses. These complement standard scoring functions. The calculations are based on the concept of empirical Molecular Hydrophobicity Potential (MHP). AVAILABILITY The PLATINUM web tool as well as detailed documentation and tutorial are available free of charge for academic users at http://model.nmr.ru/platinum/. PLATINUM requires Java 5 or higher and Adobe Flash Player 9. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Evidence that interaction between conserved residues in transmembrane helices 2, 3 and 7 are crucial for human VPAC1 receptor activation.

The VPAC1 receptor belongs to family B of G protein-coupled receptors (GPCR-B) and is activated upon binding of the vasoactive intestinal peptide (VIP). Despite the recent determination of the structure of the N terminus of several members of this receptor family, little is known about the structure of the transmembrane (TM) region and about the molecular mechanisms leading to activation. In the present study, we designed a new structural model of the TM domain and combined it with experimental mutagenesis experiments to investigate the interaction network that governs ligand binding and receptor activation. Our results suggest that this network involves the cluster of residues Arg188 in TM2, Gln380 in TM7, and Asn229 in TM3. This cluster is expected to be altered upon VIP binding, because Arg188 has been shown previously to interact with Asp3 of VIP. Several point mutations at positions 188, 229, and 380 were experimentally characterized and were shown to severely affect VIP binding and/or VIP-mediated cAMP production. Double mutants built from reciprocal residue exchanges exhibit strong cooperative or anticooperative effects, thereby indicating the spatial proximity of residues Arg188, Gln380, and Asn229. Because these residues are highly conserved in the GPCR-B family, they can moreover be expected to have a general role in mediating function.

Homology Modeling of MT1 and MT2 Receptors

Melatonin is a neurohormone synthesized and secreted mainly during the dark period of the circadian cycle by the pineal gland. It has already been proved to be involved in a number of chronobiological processes, most of them being mediated by its membranar receptors MT1 and MT2. Both are members of the GPCR class and, despite the interest they elicit, their 3D structure is still to be described. Models for both human MT1 and MT2 receptors have been constructed by homology modeling, using the X-ray structure of bovine rhodopsin as template. These models have been evaluated in terms of hydrophobic properties of the helices and refined to take into account the rearrangement of GPCRs necessary for their activation, thus leading to a putative activated model for each subtype.

Recent advances in computational modeling of α-helical membrane-active peptides.

Membrane-active peptides (MAPs) represent a broad variety of molecules, and biological functions of most are directly associated with their ability to interact with membranes. Taking into account the effect of MAPs on living cells they can be nominally divided into three major groups - fusion (FPs), antimicrobial/cytolytic (AMPs/CPs) and cell-penetrating (CPPs) peptides. Although spatial structure of different MAPs varies to a great extent, linear α-helical peptides represent the most studied class. These peptides possess relatively simple structural organization and share a set of similar molecular features, which make them very attractive to both experimental and computational studies. Here, we review different molecular modeling methods in prospective of their applications to study of α-helical MAPs. The most sophisticated of them, such as molecular dynamics simulations, give atomistic information about molecular interactions driving peptide binding to the water-lipid interface, cooperative mechanisms of membrane destabilization and thermodynamics of these processes. Significant progress has been achieved in this field during the last few years, resulting in a possibility to observe computationally MAPs action in realistic peptide-to-lipid ratios and over the microsecond timescale. Other relatively simple but powerful approaches allow assessment of important characteristics of MAPs such as α-helical propensity, amphiphilicity, total hydrophobicity, and spatial distribution of charge and hydrophobic/hydrophilic properties, etc. Altogether, computational methods provide efficient basis for rational design of MAPs with predefined properties and a spectrum of biological activities.

Kalium: a database of potassium channel toxins from scorpion venom

Kalium (http://kaliumdb.org/) is a manually curated database that accumulates data on potassium channel toxins purified from scorpion venom (KTx). This database is an open-access resource, and provides easy access to pages of other databases of interest, such as UniProt, PDB, NCBI Taxonomy Browser, and PubMed. General achievements of Kalium are a strict and easy regulation of KTx classification based on the unified nomenclature supported by researchers in the field, removal of peptides with partial sequence and entries supported by transcriptomic information only, classification of β-family toxins, and addition of a novel λ-family. Molecules presented in the database can be processed by the Clustal Omega server using a one-click option. Molecular masses of mature peptides are calculated and available activity data are compiled for all KTx. We believe that Kalium is not only of high interest to professional toxinologists, but also of general utility to the scientific community. Database URL: http://kaliumdb.org/

Differences in binding sites of two melatonin receptors help to explain their selectivity to some melatonin analogs : A molecular modeling study

Abstract Numerous diseases have been linked to the malfunction of G-protein coupled receptors (GP-CRs). Their adequate treatment requires rational design of new high-affinity and high-selectivity drugs targeting these receptors. In this work, we report three-dimensional models of the human MT1 and MT2 melatonin receptors, members of the GPCR family. The models are based on the X-ray structure of bovine rhodopsin. The computational approach employs an original procedure for optimization of receptor-ligand structures. It includes rotation of one of the transmembrane α-helices around its axis with simultaneous assessment of quality of the resulting complexes according to a number of criteria we have developed for this purpose. The optimal geometry of the receptor-ligand binding is selected based on the analysis of complementarity of hydrophobic/hydrophilic properties between the ligand and its protein environment in the binding site. The elaborated “optimized” models are employed to explore the details of protein-ligand interactions for melatonin and a number of its analogs with known affinity to MT1 and MT2 receptors. The models permit rationalization of experimental data, including those that were not used in model building. The perspectives opened by the constructed models and by the optimization procedure in the design of new drugs are discussed.

Modular Organization of α-Toxins from Scorpion Venom Mirrors Domain Structure of Their Targets - Sodium Channels

Background: Scorpion α-toxins affect voltage-gated sodium channels in both mammals and insects. Results: We perform thorough computational analyses of α-toxin molecular architecture and structure-function relationship. Conclusion: Taxon specificity of “orphan” toxins can be predicted from a structural perspective. Significance: The proposed surface mapping technique is a new tool to analyze protein-protein complexes. To gain success in the evolutionary “arms race,” venomous animals such as scorpions produce diverse neurotoxins selected to hit targets in the nervous system of prey. Scorpion α-toxins affect insect and/or mammalian voltage-gated sodium channels (Navs) and thereby modify the excitability of muscle and nerve cells. Although more than 100 α-toxins are known and a number of them have been studied into detail, the molecular mechanism of their interaction with Navs is still poorly understood. Here, we employ extensive molecular dynamics simulations and spatial mapping of hydrophobic/hydrophilic properties distributed over the molecular surface of α-toxins. It is revealed that despite the small size and relatively rigid structure, these toxins possess modular organization from structural, functional, and evolutionary perspectives. The more conserved and rigid “core module” is supplemented with the “specificity module” (SM) that is comparatively flexible and variable and determines the taxon (mammal versus insect) specificity of α-toxin activity. We further show that SMs in mammal toxins are more flexible and hydrophilic than in insect toxins. Concomitant sequence-based analysis of the extracellular loops of Navs suggests that α-toxins recognize the channels using both modules. We propose that the core module binds to the voltage-sensing domain IV, whereas the more versatile SM interacts with the pore domain in repeat I of Navs. These findings corroborate and expand the hypothesis on different functional epitopes of toxins that has been reported previously. In effect, we propose that the modular structure in toxins evolved to match the domain architecture of Navs.

Liquid but Durable: Molecular Dynamics Simulations Explain the Unique Properties of Archaeal-Like Membranes

Archaeal plasma membranes appear to be extremely durable and almost impermeable to water and ions, in contrast to the membranes of Bacteria and Eucaryota. Additionally, they remain liquid within a temperature range of 0–100°C. These are the properties that have most likely determined the evolutionary fate of Archaea, and it may be possible for bionanotechnology to adopt these from nature. In this work, we use molecular dynamics simulations to assess at the atomistic level the structure and dynamics of a series of model archaeal membranes with lipids that have tetraether chemical nature and “branched” hydrophobic tails. We conclude that the branched structure defines dense packing and low water permeability of archaeal-like membranes, while at the same time ensuring a liquid-crystalline state, which is vital for living cells. This makes tetraether lipid systems promising in bionanotechnology and material science, namely for design of new and unique membrane nanosystems.

Temperature-sensitive gating of TRPV1 channel as probed by atomistic simulations of its trans- and juxtamembrane domains

Heat-activated transient receptor potential channel TRPV1 is one of the most studied eukaryotic proteins involved in temperature sensation. Upon heating, it exhibits rapid reversible pore gating, which depolarizes neurons and generates action potentials. Underlying molecular details of such effects in the pore region of TRPV1 is of a crucial importance to control temperature responses of the organism. Despite the spatial structure of the channel in both open (O) and closed (C) states is known, microscopic nature of channel gating and mechanism of thermal sensitivity are still poorly understood. In this work, we used unrestrained atomistic molecular dynamics simulations of TRPV1 (without N- and C-terminal cytoplasmic domains) embedded into explicit lipid bilayer in its O- and C-states. We found that the pore domain with its neighboring loops undergoes large temperature-dependent conformational transitions in an asymmetric way, when fragments of only one monomer move with large amplitude, freeing the pore upon heating. Such an asymmetrical gating looks rather biologically relevant because it is faster and more reliable than traditionally proposed “iris-like” symmetric scheme of channel opening. Analysis of structural, dynamic, and hydrophobic organization of the pore domain revealed entropy growth upon TRPV1 gating, which is in line with current concepts of thermal sensitivity.

Structural Insight into Specificity of Interactions between Nonconventional Three-finger Weak Toxin from Naja kaouthia (WTX) and Muscarinic Acetylcholine Receptors.

Background: Cobras “three-finger” nonconventional toxin WTX allosterically modulates muscarinic receptors (mAChRs). Results: Activity of several WTX mutants was analyzed; toxin spatial structure and dynamics were determined; and complexes of toxin with M1 and M3 mAChRs were modeled. Conclusion: Flexible loop II is the major determinant for toxin binding to different mAChRs. Significance: Structural framework for rationalization of target-specific positive/negative allosteric regulation of mAChRs is provided. Weak toxin from Naja kaouthia (WTX) belongs to the group of nonconventional “three-finger” snake neurotoxins. It irreversibly inhibits nicotinic acetylcholine receptors and allosterically interacts with muscarinic acetylcholine receptors (mAChRs). Using site-directed mutagenesis, NMR spectroscopy, and computer modeling, we investigated the recombinant mutant WTX analogue (rWTX) which, compared with the native toxin, has an additional N-terminal methionine residue. In comparison with the wild-type toxin, rWTX demonstrated an altered pharmacological profile, decreased binding of orthosteric antagonist N-methylscopolamine to human M1- and M2-mAChRs, and increased antagonist binding to M3-mAChR. Positively charged arginine residues located in the flexible loop II were found to be crucial for rWTX interactions with all types of mAChR. Computer modeling suggested that the rWTX loop II protrudes to the M1-mAChR allosteric ligand-binding site blocking the entrance to the orthosteric site. In contrast, toxin interacts with M3-mAChR by loop II without penetration into the allosteric site. Data obtained provide new structural insight into the target-specific allosteric regulation of mAChRs by “three-finger” snake neurotoxins.