Deciphering fine molecular details of proteins' structure and function with a Protein Surface Topography (PST) method

Broad-spectrum activity of membranolytic cationic macrocyclic peptides against multi-drug resistant bacteria and fungi

The emergence of multidrug-resistant (MDR) strains causes severe problems in the treatment of microbial infections owing to limited treatment options. Antimicrobial peptides (AMPs) are drawing considerable attention as promising antibiotic alternative candidates to combat MDR bacterial and fungal infections. Herein, we present a series of small amphiphilic membrane-active cyclic peptides composed, in part, of various nongenetically encoded hydrophilic and hydrophobic amino acids. Notably, lead cyclic peptides 3b and 4b showed broad-spectrum activity against drug-resistant Gram-positive (MIC = 1.5-6.2 µg/mL) and Gram-negative (MIC = 12.5-25 µg/mL) bacteria, and fungi (MIC = 3.1-12.5 µg/mL). Furthermore, lead peptides displayed substantial antibiofilm action comparable to standard antibiotics. Hemolysis (HC50 = 230 µg/mL) and cytotoxicity (>70 % cell viability against four different mammalian cells at 100 µg/mL) assay results demonstrated the selective lethal action of 3b against microbes over mammalian cells. A calcein dye leakage experiment substantiated the membranolytic effect of 3b and 4b, which was further confirmed by scanning electron microscopy. The behavior of 3b and 4b in aqueous solution and interaction with phospholipid bilayers were assessed by employing nuclear magnetic resonance (NMR) spectroscopy in conjunction with molecular dynamics (MD) simulations, providing a solid structural basis for understanding their membranolytic action. Moreover, 3b exhibited stability in human blood plasma (t1/2 = 13 h) and demonstrated no signs of resistance development against antibiotic-resistant S. aureus and E. coli. These findings underscore the potential of these newly designed amphiphilic cyclic peptides as promising anti-infective agents, especially against Gram-positive bacteria.

Dynamic molecular portraits of ion-conducting pores characterize functional states of TRPV channels

Structural biology is solving an ever-increasing number of snapshots of ion channel conformational ensembles. Deciphering ion channel mechanisms, however, requires understanding the ensemble dynamics beyond the static structures. Here, we present a molecular modeling-based approach characterizing the ion channel structural intermediates, or their "dynamic molecular portraits", by assessing water and ion conductivity along with the detailed evaluation of pore hydrophobicity and residue packing. We illustrate the power of this approach by analyzing structures of few vanilloid-subfamily transient receptor potential (TRPV) channels. Based on the pore architecture, there are three major states that are common for TRPVs, which we call α-closed, π-closed, and π-open. We show that the pore hydrophobicity and residue packing for the open state is most favorable for the pore conductance. On the contrary, the α-closed state is the most hydrophobic and always non-conducting. Our approach can also be used for structural and functional classification of ion channels.

The scorpion toxin BeKm-1 blocks hERG cardiac potassium channels using an indispensable arginine residue

BeKm-1 is a peptide toxin from scorpion venom that blocks the pore of the potassium channel hERG (Kv11.1) in the human heart. Although individual protein structures have been resolved, the structure of the complex between hERG and BeKm-1 is unknown. Here, we used molecular dynamics and ensemble docking, guided by previous double-mutant cycle analysis data, to obtain an in silico model of the hERG-BeKm-1 complex. Adding to the previous mutagenesis study of BeKm-1, our model uncovers the key role of residue Arg20, which forms three interactions (a salt bridge and hydrogen bonds) with the channel vestibule simultaneously. Replacement of this residue even by lysine weakens the interactions significantly. In accordance, the recombinantly produced BeKm-1R20K mutant exhibited dramatically decreased activity on hERG. Our model may be useful for future drug design attempts.

Structural Analysis and Activity Correlation of Amphiphilic Cyclic Antimicrobial Peptides Derived from the [W4R4] Scaffold

In our ongoing quest to design effective antimicrobial peptides (AMPs), this study aimed to elucidate the mechanisms governing cyclic amphiphilic AMPs and their interactions with membranes. The objective was to discern the nature of these interactions and understand how peptide sequence and structure influence antimicrobial activity. We introduced modifications into the established cyclic AMP peptide, [W4R4], incorporating an extra aromatic hydrophobic residue (W), a positively charged residue (R), or the unique 2,5-diketopiperazine (DKP). This study systematically explored the structure-activity relationships (SARs) of a series of cyclic peptides derived from the [W4R4] scaffold, including the first synthesis and evaluation of [W4R4(DKP)]. Structural, dynamic, hydrophobic, and membrane-binding properties of four cyclic peptides ([W4R4], [W5R4], [W4R5], [W4R4(DKP)]) were explored using molecular dynamics simulations within a DOPC/DOPG lipid bilayer that mimics the bacterial membrane. The results revealed distinct SARs linking antimicrobial activity to parameters such as conformational plasticity, immersion depth in the bilayer, and population of the membrane binding mode. Notably, [W4R5] exhibited an optimal "activity/binding to the bacterial membrane" pattern. This multidisciplinary approach efficiently decoded finely regulated SAR profiles, laying a foundation for the rational design of novel antimicrobial peptides.

Molecular pathway and structural mechanism of human oncochannel TRPV6 inhibition by the phytocannabinoid tetrahydrocannabivarin

The calcium-selective oncochannel TRPV6 is an important driver of cell proliferation in human cancers. Despite increasing interest of pharmacological research in developing synthetic inhibitors of TRPV6, natural compounds acting at this channel have been largely neglected. On the other hand, pharmacokinetics of natural small-molecule antagonists optimized by nature throughout evolution endows these compounds with a medicinal potential to serve as potent and safe next-generation anti-cancer drugs. Here we report the structure of human TRPV6 in complex with tetrahydrocannabivarin (THCV), a natural cannabinoid inhibitor extracted from Cannabis sativa. We use cryo-electron microscopy combined with electrophysiology, calcium imaging, mutagenesis, and molecular dynamics simulations to identify THCV binding sites in the portals that connect the membrane environment surrounding the protein to the central cavity of the channel pore and to characterize the allosteric mechanism of TRPV6 inhibition. We also propose the molecular pathway taken by THCV to reach its binding site. Our study provides a foundation for the development of new TRPV6-targeting drugs.

Mixed mechanism of conformational selection and induced fit as a molecular recognition process in the calreticulin family of proteins

The fundamental question on the mechanism of molecular recognition during ligand binding has attracted a lot of scientific scrutiny. The two competing theories of ligand binding-"induced fit" and "conformational selection" have been proposed to explain biomolecular recognition. Since exploring a family of proteins with similar structural architectures and conserved functional roles can provide valuable insight into the significance of molecular structure and function, we performed molecular dynamics simulations on the calreticulin family of proteins, which specifically recognize monoglucosylated N-glycan during the protein folding process. Atomistic simulations of lectins in free and bound forms demonstrated that they exist in several conformations spanning from favorable to unfavorable for glycan binding. Our analysis was confined to the carbohydrate recognition domain (CRD) of these lectins to demonstrate the degree of conservation in protein sequence and structure and relate them with their function. Furthermore, we computed the lectin-glycan binding affinity using the mmPBSA approach to identify the most favorable lectin conformation for glycan binding and compared the molecular interaction fields in terms of noncovalent bond interactions. We also demonstrated the involvement of Tyr and Trp residues in the CRD with the non-reducing end glucose and central mannose residues, which contribute to some of the specific interactions. Furthermore, we analyzed the conformational changes in the CRD through SASA, RMSFs and protein surface topography mapping of electrostatic and hydrophobic potentials. Our findings demonstrate a hybrid mechanism of molecular recognition, initially driven by conformational selection followed by glycan-induced fluctuations in the key residues to strengthen the glycan binding interactions.

PST-PRNA: prediction of RNA-binding sites using protein surface topography and deep learning

MOTIVATION

Protein-RNA interactions play essential roles in many biological processes, including pre-mRNA processing, post-transcriptional gene regulation and RNA degradation. Accurate identification of binding sites on RNA-binding proteins (RBPs) is important for functional annotation and site-directed mutagenesis. Experimental assays to sparse RBPs are precise and convincing but also costly and time consuming. Therefore, flexible and reliable computational methods are required to recognize RNA-binding residues.

RESULTS

In this work, we propose PST-PRNA, a novel model for predicting RNA-binding sites (PRNA) based on protein surface topography (PST). Taking full advantage of the 3D structural information of protein, PST-PRNA creates representative topography images of the entire protein surface by mapping it onto a unit spherical surface. Four kinds of descriptors are encoded to represent residues on the surface. Then, the potential features are integrated and optimized by using deep learning models. We compile a comprehensive non-redundant RBP dataset to train and test PST-PRNA using 10-fold cross-validation. Numerous experiments demonstrate PST-PRNA learns successfully the latent structural information of protein surface. On the non-redundant dataset with sequence identity of 0.3, PST-PRNA achieves area under the receiver operating characteristic curves (AUC) value of 0.860 and Matthew's correlation coefficient value of 0.420. Furthermore, we construct a completely independent test dataset for justification and comparison. PST-PRNA achieves AUC value of 0.913 on the independent dataset, which is superior to the other state-of-the-art methods.

AVAILABILITY AND IMPLEMENTATION

The code and data are available at https://www.github.com/zpliulab/PST-PRNA. A web server is freely available at http://www.zpliulab.cn/PSTPRNA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Dynamic "Molecular Portraits" of Biomembranes Drawn by Their Lateral Nanoscale Inhomogeneities

To date, it has been reliably shown that the lipid bilayer/water interface can be thoroughly characterized by a sophisticated so-called "dynamic molecular portrait". The latter reflects a combination of time-dependent surface distributions of various physicochemical properties, inherent in both model lipid bilayers and natural multi-component cell membranes. One of the most important features of biomembranes is their mosaicity, which is expressed in the constant presence of lateral inhomogeneities, the sizes and lifetimes of which vary in a wide range-from 1 to 10³ nm and from 0.1 ns to milliseconds. In addition to the relatively well-studied macroscopic domains (so-called "rafts"), the analysis of micro- and nanoclusters (or domains) that form an instantaneous picture of the distribution of structural, dynamic, hydrophobic, electrical, etc., properties at the membrane-water interface is attracting increasing interest. This is because such nanodomains (NDs) have been proven to be crucial for the proper membrane functioning in cells. Therefore, an understanding with atomistic details the phenomena associated with NDs is required. The present mini-review describes the recent results of experimental and in silico studies of spontaneously formed NDs in lipid membranes. The main attention is paid to the methods of ND detection, characterization of their spatiotemporal parameters, the elucidation of the molecular mechanisms of their formation. Biological role of NDs in cell membranes is briefly discussed. Understanding such effects creates the basis for rational design of new prospective drugs, therapeutic approaches, and artificial membrane materials with specified properties.

Potassium channel blocker crafted by α-hairpinin scaffold engineering

The α-Hairpinins are a family of plant defense peptides with a common fold presenting two short α-helices stabilized by two invariant S-S-bridges. We have shown previously that substitution of just two amino acid residues in a wheat α-hairpinin Tk-AMP-X2 leads to Tk-hefu-2 that features specific affinity to voltage-gated potassium channels K_V1.3. Here, we utilize a combined molecular modeling approach based on molecular dynamics simulations and protein surface topography technique to improve the affinity of Tk-hefu-2 to K_V1.3 while preserving its specificity. An important advance of this work compared with our previous studies is transition from the analysis of various physicochemical properties of an isolated toxin molecule to its consideration in complex with its target, a membrane-bound ion channel. As a result, a panel of computationally designed Tk-hefu-2 derivatives was synthesized and tested against K_V1.3. The most active mutant Tk-hefu-10 showed a half-maximal inhibitory concentration of ∼150 nM being >10 times more active than Tk-hefu-2 and >200 times more active than the original Tk-hefu. We conclude that α-hairpinins provide an attractive disulfide-stabilized scaffold for the rational design of ion channel inhibitors. Furthermore, the success rate can be considerably increased by the proposed "target-based" iterative strategy of molecular design.

Mapping the underlying mechanisms of fibrinogen benzothiazole drug interactions using computational and experimental approaches

Three-dimensional conformational crystallographic binding-modes are of paramount importance to understand the docking mechanism of protein-ligand interactions and to identify potential "leading drugs" conformers towards rational drugs-design. Herein, we present an integrated computational-experimental study tackling the problem of multiple binding modes among the ligand 3-(2-Benzothiazolylthio)-propane sulfonic acid (BTS) and the fibrinogen receptor (E-region). Based on molecular docking simulations, we found that the free energy of binding values for nine of different BTS-docking complexes (i.e., BTS-pose_1-9) were very close. We have also identified a docking-mechanism of BTS-interaction mainly based on non-covalent hydrophobic interactions with H-bond contacts stabilizing the fibrinogen-BTS docking complexes. Interestingly, the different BTS-poses_1-9 were found to be able to block the fibrinogen binding site (E-region) by inducing local perturbations in effector and allosteric residues, reducing the degree of collectivity in its flexibility normal modes. As such, we theoretically suggest that the BTS-binding modes can significantly affect the physiological condition of the unoccupied fibrinogen protein structure by bringing global and local perturbations in the frequency domain spectra. The proposed theoretical mechanisms, the interactions involved and the conformational changes suggested, were further corroborated by different experimental techniques such as isothermal titration calorimetry (ITC), zeta potential, UV-vis, fluorescence and small angle X-ray scattering (SAXS). The combined results shall open new avenues towards the application of complex supra-molecular information in rational drugs-design.

Ribbon α-Conotoxin KTM Exhibits Potent Inhibition of Nicotinic Acetylcholine Receptors

KTM is a 16 amino acid peptide with the sequence WCCSYPGCYWSSSKWC. Here, we present the nuclear magnetic resonance (NMR) structure and bioactivity of this rationally designed α-conotoxin (α-CTx) that demonstrates potent inhibition of rat α3β2-nicotinic acetylcholine receptors (rα3β2-nAChRs). Two bioassays were used to test the efficacy of KTM. First, a qualitative PC12 cell-based assay confirmed that KTM acts as a nAChR antagonist. Second, bioactivity evaluation by two-electrode voltage clamp electrophysiology was used to measure the inhibition of rα3β2-nAChRs by KTM (IC50 = 0.19 ± 0.02 nM), and inhibition of the same nAChR isoform by α-CTx MII (IC50 = 0.35 ± 0.8 nM). The three-dimensional structure of KTM was determined by NMR spectroscopy, and the final set of 20 structures derived from 32 distance restraints, four dihedral angle constraints, and two disulfide bond constraints overlapped with a mean global backbone root-mean-square deviation (RMSD) of 1.7 ± 0.5 Å. The structure of KTM did not adopt the disulfide fold of α-CTx MII for which it was designed, but instead adopted a flexible ribbon backbone and disulfide connectivity of C2-C16 and C3-C8 with an estimated 12.5% α-helical content. In contrast, α-CTx MII, which has a native fold of C2-C8 and C3-C16, has an estimated 38.1% α-helical secondary structure. KTM is the first reported instance of a Framework I (CC-C-C) α-CTx with ribbon connectivity to display sub-nanomolar inhibitory potency of rα3β2-nAChR subtypes.

Qualitative Assay to Detect Dopamine Release by Ligand Action on Nicotinic Acetylcholine Receptors

A pheochromocytoma of the rat adrenal medulla derived (a.k.a. PC12) cell-based assay for dopamine measurement by luminescence detection was customized for the qualitative evaluation of agonists and antagonists of nicotinic acetylcholine receptors (nAChRs). The assay mechanism begins with ligand binding to transmembrane nAChRs, altering ion flow into the cell and inducing dopamine release from the cell. Following release, dopamine is oxidized by monoamine oxidase generating hydrogen peroxide that catalyzes a chemiluminescence reaction involving luminol and horseradish peroxidase, thus producing a detectable response. Results are presented for the action of nAChR agonists (acetylcholine, nicotine, and cytisine), and antagonists (α-conotoxins (α-CTxs) MII, ImI, LvIA, and PeIA) that demonstrate a luminescence response correlating to the increase or decrease of dopamine release. A survey of cell growth and treatment conditions, including nerve growth factor, nicotine, ethanol, and temperature, led to optimal assay requirements to achieve maximal signal intensity and consistent response to ligand treatment. It was determined that PC12 cells treated with a combination of nerve growth factor and nicotine, and incubated at 37 °C, provided favorable results for a reduction in luminescence signal upon treatment of cells with α-CTxs. The PC12 assay is intended for use as a fast, efficient, and economic qualitative method to assess the bioactivity of molecules that act on nAChRs, in which testing of ligand-nAChR binding hypotheses and computational predictions can be validated. As a screening method for nAChR bioactivity, lead compounds can be assessed for their likelihood of exhibiting desired bioactivity prior to being subjected to more complex quantitative methods, such as electrophysiology or live animal studies.

Protein surface topography as a tool to enhance the selective activity of a potassium channel blocker

Tk-hefu is an artificial peptide designed based on the α-hairpinin scaffold, which selectively blocks voltage-gated potassium channels K_v1.3. Here we present its spatial structure resolved by NMR spectroscopy and analyze its interaction with channels using computer modeling. We apply protein surface topography to suggest mutations and increase Tk-hefu affinity to the K_v1.3 channel isoform. We redesign the functional surface of Tk-hefu to better match the respective surface of the channel pore vestibule. The resulting peptide Tk-hefu-2 retains K_v1.3 selectivity and displays ∼15 times greater activity compared with Tk-hefu. We verify the mode of Tk-hefu-2 binding to the channel outer vestibule experimentally by site-directed mutagenesis. We argue that scaffold engineering aided by protein surface topography represents a reliable tool for design and optimization of specific ion channel ligands.

Mutagenesis of α-Conotoxins for Enhancing Activity and Selectivity for Nicotinic Acetylcholine Receptors

Nicotinic acetylcholine receptors (nAChRs) are found throughout the mammalian body and have been studied extensively because of their implication in a myriad of diseases. α-Conotoxins (α-CTxs) are peptide neurotoxins found in the venom of marine snails of genus Conus. α-CTxs are potent and selective antagonists for a variety of nAChR isoforms. Over the past 40 years, α-CTxs have proven to be valuable molecular probes capable of differentiating between closely related nAChR subtypes and have contributed greatly to understanding the physiological role of nAChRs in the mammalian nervous system. Here, we review the amino acid composition and structure of several α-CTxs that selectively target nAChR isoforms and explore strategies and outcomes for introducing mutations in native α-CTxs to direct selectivity and enhance binding affinity for specific nAChRs. This review will focus on structure-activity relationship studies involving native α-CTxs that have been rationally mutated and molecular interactions that underlie binding between ligand and nAChR isoform.

Vasoinhibin comprises a three-helix bundle and its antiangiogenic domain is located within the first 79 residues

Vasoinhibin belongs to a family of angiogenesis inhibitors generated when the fourth α-helix (H4) of the hormone prolactin (PRL) is removed by specific proteolytic cleavage. The antiangiogenic properties are absent in uncleaved PRL, indicating that conformational changes create a new bioactive domain. However, the solution structure of vasoinhibin and the location of its bioactive domain are unknown. Molecular dynamic simulation (MD) showed that the loss of H4 exposes the hydrophobic nucleus of PRL and leads to the compression of the molecule into a three-helix bundle that buries the hydrophobic nucleus again. Compression occurs by the movement of loop 1 (L1) and its interaction with α-helix 1 (H1) generating a new L1 conformation with electrostatic and hydrophobic surfaces distinct from those of PRL, that may correspond to a bioactive domain. Consistent with this model, a recombinant protein containing the first 79 amino acids comprising H1 and L1 of human PRL inhibited the proliferation and migration of endothelial cells and upregulated the vasoinhibin target genes, IL1A and ICAM1. This bioactivity was comparable to that of a conventional vasoinhibin having the 123 residues encompassing H1, L1, Η2, L2, and Η3 of human PRL. These findings extend the vasoinhibin family to smaller proteins and provide important structural information, which will aid in antiangiogenic drug development.

The role of hydrophobic /hydrophilic balance in the activity of structurally flexible vs. rigid cytolytic polypeptides and analogs developed on their basis

INTRODUCTION

Being important representatives of various proteomes, membrane-active cationic peptides (CPs) are attractive objects as lead compounds in the design of new antibacterial, anticancer, antifungal, and antiviral molecules. Numerous CPs are found in insect and snake venoms, where many of them reveal cytolytic properties. Due to advances in omics technologies, the number of such peptides is growing dramatically. Areas covered: To understand structure-function relationships for CPs in a living cell, detailed analysis of their hydrophobic/hydrophilic properties is indispensable. We consider two structural classes of membrane-active CPs: latarcins (Ltc) from spider and cardiotoxins (CTXs) from snake venoms. While the former are void off disulfide bonds and conformationally flexible, the latter are structurally rigid and cross-linked with disulfide bonds. In order to elucidate structure-activity relationships behind their antibacterial, anticancer, and hemolytic effects, the properties of these polypeptides are considered on a side-by-side basis. Expert commentary: An ever-increasing number of venom-derived membrane-active polypeptides require new methods for identification of their functional propensities and sequence-based design of novel pharmacological substances. We address these issues considering a number of the designed peptides, based either on Ltc or CTX sequences. Experimental and computer modeling techniques required for these purposes are delineated.

CARD domain of rat RIP2 kinase: Refolding, solution structure, pH-dependent behavior and protein-protein interactions

RIP2, one of the RIP kinases, interacts with p75 neurotrophin receptor, regulating the neuron survival, and with NOD1 and NOD2 proteins, causing the innate immune response against gram-negative and gram-positive bacteria via its caspase recruitment domain (CARD). This makes RIP2 a prospective target for novel therapies, aimed to modulate the inflammatory diseases and neurogenesis/neurodegeneration. Several studies report the problems with the stability of human RIP2 CARD and its production in bacterial hosts, which is a prerequisite for the structural investigation with solution NMR spectroscopy. In the present work, we report the high yield production and refolding protocols and resolve the structure of rat RIP2 CARD. The structure reveals the important differences to the previously published conformation of the homologous human protein. Using solution NMR, we characterized the intramolecular mobility and pH-dependent behavior of RIP2 CARD, and found the propensity of the protein to form high-order oligomers at physiological pH while being monomeric under acidic conditions. The oligomerization of protein may be explained, based on the electrostatic properties of its surface. Analysis of the structure and sequences of homologous proteins reveals the residues which are significant for the unusual fold of RIP2 CARD domains from different species. The high-throughput protein production/refolding protocols and proposed explanation for the protein oligomerization, provide an opportunity to design the stabilized variants of RIP2 CARD, which could be used to study the structural details of RIP2/NOD1/NOD2 interaction and perform the rational drug design.

Mixed-Probe Simulation and Probe-Derived Surface Topography Map Analysis for Ligand Binding Site Identification

Membrane proteins represent a considerable fraction of pharmaceutical drug targets. A computational technique to identify ligand binding pockets in these proteins is therefore of great importance. We recently reported such a technique called pMD-membrane that utilizes small molecule probes to detect ligand binding sites and surface hotspots on membrane proteins based on probe-based molecular dynamics simulation. The current work extends pMD-membrane to a diverse set of small organic molecular species that can be used as cosolvents during simulation of membrane proteins. We also describe a projection technique for globally quantifying probe densities on the protein surface and introduce a technique to construct surface topography maps directly from the probe-binding propensity of surface residues. The map reveals surface patterns and geometric features that aid in filtering out high probe density hotspots lacking pocketlike characteristics. We demonstrate the applicability of the extended pMD-membrane and the new analysis tool by exploring the druggability of full-length G12D, G12V, and G13D oncogenic K-Ras mutants bound to a negatively charged lipid bilayer. Using data from 30 pMD-membrane runs conducted in the presence of a 2.8 M cosolvent made up of an equal proportion of seven small organic molecules, we show that our approach robustly identifies known allosteric ligand binding sites and other reactive regions on K-Ras. Our results also show that accessibility of some pockets is modulated by differential membrane interactions.

High-Affinity α-Conotoxin PnIA Analogs Designed on the Basis of the Protein Surface Topography Method

Despite some success for small molecules, elucidating structure-function relationships for biologically active peptides - the ligands for various targets in the organism - remains a great challenge and calls for the development of novel approaches. Some of us recently proposed the Protein Surface Topography (PST) approach, which benefits from a simplified representation of biomolecules' surface as projection maps, which enables the exposure of the structure-function dependencies. Here, we use PST to uncover the "activity pattern" in α-conotoxins - neuroactive peptides that effectively target nicotinic acetylcholine receptors (nAChRs). PST was applied in order to design several variants of the α-conotoxin PnIA, which were synthesized and thoroughly studied. Among the best was PnIA[R9, L10], which exhibits nanomolar affinity for the α7 nAChR, selectivity and a slow wash-out from this target. Importantly, these mutations could hardly be delineated by "standard" structure-based drug design. The proposed combination of PST with a set of experiments proved very efficient for the rational construction of new bioactive molecules.

Modeling phase transitions in mixtures of β-γ lens crystallins

We analyze the experimentally determined phase diagram of a γD-βB1 crystallin mixture. Proteins are described as dumbbells decorated with attractive sites to allow inter-particle interaction. We use thermodynamic perturbation theory to calculate the free energy of such mixtures and, by applying equilibrium conditions, also the compositions and concentrations of the co-existing phases. Initially we fit the Tcloudversus packing fraction η measurements for a pure (x2 = 0) γD solution in 0.1 M phosphate buffer at pH = 7.0. Another piece of experimental data, used to fix the model parameters, is the isotherm x2vs. η at T = 268.5 K, at the same pH and salt content. We use the conventional Lorentz-Berthelot mixing rules to describe cross interactions. This enables us to determine: (i) model parameters for pure βB1 crystallin protein and to calculate; (ii) complete equilibrium surface (Tcloud-x2-η) for the crystallin mixtures. (iii) We present the results for several isotherms, including the tie-lines, as also the temperature-packing fraction curves. Good agreement with the available experimental data is obtained. An interesting result of these calculations is evidence of the coexistence of three phases. This domain appears for the region of temperatures just out of the experimental range studied so far. The input parameters, leading good description of experimental data, revealed a large difference between the numbers of the attractive sites for γD and βB1 proteins. This interesting result may be related to the fact that γD has a more than nine times smaller quadrupole moment than its partner in the mixture.

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.

Structuprint: a scalable and extensible tool for two-dimensional representation of protein surfaces

BACKGROUND

The term 'molecular cartography' encompasses a family of computational methods for two-dimensional transformation of protein structures and analysis of their physicochemical properties. The underlying algorithms comprise multiple manual steps, whereas the few existing implementations typically restrict the user to a very limited set of molecular descriptors.

RESULTS

We present Structuprint, a free standalone software that fully automates the rendering of protein surface maps, given - at the very least - a directory with a PDB file and an amino acid property. The tool comes with a default database of 328 descriptors, which can be extended or substituted by user-provided ones. The core algorithm comprises the generation of a mould of the protein surface, which is subsequently converted to a sphere and mapped to two dimensions, using the Miller cylindrical projection. Structuprint is partly optimized for multicore computers, making the rendering of animations of entire molecular dynamics simulations feasible.

CONCLUSIONS

Structuprint is an efficient application, implementing a molecular cartography algorithm for protein surfaces. According to the results of a benchmark, its memory requirements and execution time are reasonable, allowing it to run even on low-end personal computers. We believe that it will be of use - primarily but not exclusively - to structural biologists and computational biochemists.

Natural compounds interacting with nicotinic acetylcholine receptors: from low-molecular weight ones to peptides and proteins

Nicotinic acetylcholine receptors (nAChRs) fulfill a variety of functions making identification and analysis of nAChR subtypes a challenging task. Traditional instruments for nAChR research are d-tubocurarine, snake venom protein α-bungarotoxin (α-Bgt), and α-conotoxins, neurotoxic peptides from Conus snails. Various new compounds of different structural classes also interacting with nAChRs have been recently identified. Among the low-molecular weight compounds are alkaloids pibocin, varacin and makaluvamines C and G. 6-Bromohypaphorine from the mollusk Hermissenda crassicornis does not bind to Torpedo nAChR but behaves as an agonist on human α7 nAChR. To get more selective α-conotoxins, computer modeling of their complexes with acetylcholine-binding proteins and distinct nAChRs was used. Several novel three-finger neurotoxins targeting nAChRs were described and α-Bgt inhibition of GABA-A receptors was discovered. Information on the mechanisms of nAChR interactions with the three-finger proteins of the Ly6 family was found. Snake venom phospholipases A₂ were recently found to inhibit different nAChR subtypes. Blocking of nAChRs in Lymnaea stagnalis neurons was shown for venom C-type lectin-like proteins, appearing to be the largest molecules capable to interact with the receptor. A huge nAChR molecule sensible to conformational rearrangements accommodates diverse binding sites recognizable by structurally very different compounds.