5β-reduced neuroactive steroids as modulators of growth and viability of postnatal neurons and glia

Endogenous neurosteroids (NS) and their synthetic analogs, neuroactive steroids (NAS), are potentially useful drug-like compounds affecting the pathophysiology of miscellaneous central nervous system disorders (e.g. Alzheimer´s disease, epilepsy, depression, etc.). Additionally, NS have been shown to promote neuron viability and neurite outgrowth upon injury. The molecular, structural and physicochemical basis of the NS effect on neurons is so far not fully understood, and the development of new, biologically relevant assays is essential for their comparative analysis and for assessment of their mechanism of action. Here, we report the development of a novel, plate-based, high-content in vitro assay for screening of NS and newly synthesized, 5β-reduced NAS for the promotion of postnatal neuron survival and neurite growth using fluorescent, postnatal mixed cortical neuron cultures isolated from thy1-YFP transgenic mice. The screen allows a detailed time course analysis of different parameters, such as the number of neurons or neurite lengths of 7-day, in vitro neuron cultures. Using the screen, we identify a new NAS, compound 42, that promotes the survival and growth of postnatal neurons significantly better than several endogenous NS (dehydroepiandrosterone, progesterone, and allopregnanolone). Interestingly, we demonstrate that compound 42 also promotes the proliferation of glia (in particular oligodendrocytes) and that the glial function is critical for its neuron growth support. Computational analysis of the biological data and calculated physicochemical properties of tested NS and NAS demonstrated that their biological activity is proportional to their lipophilicity. Together, the screen proves useful for the selection of neuron-active NAS and the comparative evaluation of their biologically relevant structural and physicochemical features.

Proteolytic profiles of two isoforms of human AMBN expressed in E. coli by MMP-20 and KLK-4 proteases

Ameloblastin is a protein in biomineralization of tooth enamel. However recent results indicate that this is probably not its only role in an organism. Enamel matrix formation represents a complex process enabled via specific crosslinking of two proteins - the most abundant amelogenin and the ameloblastin (AMBN). The human AMBN (hAMBN) gene possesses 13 protein coding exons with alternatively spliced transcripts and the longest isoform about 447 amino acid residues. It has been described that AMBN molecules in vitro assemble into oligomers via a sequence encoded by exon 5. Enamel is formed by the processing of enamel proteins by two specific proteases - enamelysin (MMP-20) and kallikrein 4 (KLK-4). The scaffold made of AMEL and non-amelogenin proteins is cleaved and removed from the developed tooth enamel. The hAMBN is expressed in two isoforms (ISO I and II), which could lead to their different utilization determined by distinct proteolytic profiles. In this study, we compared proteolytic profiles of both isoforms of hAMBN expressed in E. coli after proteolysis by MMP-20, KLK-4, and their 1:2 mixture. Proteolysis products were analysed and cleavage sites were identified by mass spectrometry. The proteolytic profiles of two AMBN isoforms showed different results, although we have to determine that the analysed AMBN was not post-translationally modified as expressed in prokaryotic cells. These results may lead to the suggestion of potentially divergent roles of AMBN isoforms cleavage products in various cell signalling pathways such as calcium buffering or signalling cascades.

Fusion of two unrelated protein domains in a chimera protein and its 3D prediction: Justification of the x-ray reference structures as a prediction benchmark

Proteins are naturally formed by domains edging their functional and structural properties. A domain out of the context of an entire protein can retain its structure and to some extent also function on its own. These properties rationalize construction of artificial fusion multidomain proteins with unique combination of various functions. Information on the specific functional and structural characteristics of individual domains in the context of new artificial fusion proteins is inevitably encoded in sequential order of composing domains defining their mutual spatial positions. So the challenges in designing new proteins with new domain combinations lie dominantly in structure/function prediction and its context dependency. Despite the enormous body of publications on artificial fusion proteins, the task of their structure/function prediction is complex and nontrivial. The degree of spatial freedom facilitated by a linker between domains and their mutual orientation driven by noncovalent interactions is beyond a simple and straightforward methodology to predict their structure with reasonable accuracy. In the presented manuscript, we tested methodology using available modeling tools and computational methods. We show that the process and methodology of such prediction are not straightforward and must be done with care even when recently introduced AlphaFold II is used. We also addressed a question of benchmarking standards for prediction of multidomain protein structures-x-ray or Nuclear Magnetic Resonance experiments. On the study of six two-domain protein chimeras as well as their composing domains and their x-ray structures selected from PDB, we conclude that the major obstacle for justified prediction is inappropriate sampling of the conformational space by the explored methods. On the other hands, we can still address particular steps of the methodology and improve the process of chimera proteins prediction.

TRPM5 Channel Binds Calcium-Binding Proteins Calmodulin and S100A1

Melastatin transient receptor potential (TRPM) channels belong to one of the most significant subgroups of the transient receptor potential (TRP) channel family. Here, we studied the TRPM5 member, the receptor exposed to calcium-mediated activation, resulting in taste transduction. It is known that most TRP channels are highly modulated through interactions with extracellular and intracellular agents. The binding sites for these ligands are usually located at the intracellular N- and C-termini of the TRP channels, and they can demonstrate the character of an intrinsically disordered protein (IDP), which allows such a region to bind various types of molecules. We explored the N-termini of TRPM5 and found the intracellular regions for calcium-binding proteins (CBPs) the calmodulin (CaM) and calcium-binding protein S1 (S100A1) by in vitro binding assays. Furthermore, molecular docking and molecular dynamics simulations (MDs) of the discovered complexes confirmed their known common binding interface patterns and the uniqueness of the basic residues present in the TRPM binding regions for CaM/S100A1.

TRPM7 N-terminal region forms complexes with calcium binding proteins CaM and S100A1

Transient receptor potential melastatin 7 (TRPM7) represents melastatin TRP channel with two significant functions, cation permeability and kinase activity. TRPM7 is widely expressed among tissues and is therefore involved in a variety of cellular functions representing mainly Mg²⁺ homeostasis, cellular Ca²⁺ flickering, and the regulation of DNA transcription by a cleaved kinase domain translocated to the nucleus. TRPM7 participates in several important biological processes in the nervous and cardiovascular systems. Together with the necessary function of the TRPM7 in these tissues and its recently analyzed overall structure, this channel requires further studies leading to the development of potential therapeutic targets. Here we present the first study investigating the N-termini of TRPM7 with binding regions for important intracellular modulators calmodulin (CaM) and calcium-binding protein S1 (S100A1) using in vitro and in silico approaches. Molecular simulations of the discovered complexes reveal their potential binding interfaces with common interaction patterns and the important role of basic residues present in the N-terminal binding region of TRPM.

The order of PDZ3 and TrpCage in fusion chimeras determines their properties-a biophysical characterization

Most of the structural proteins known today are composed of domains that carry their own functions while keeping their structural properties. It is supposed that such domains, when taken out of the context of the whole protein, can retain their original structure and function to a certain extent. Information on the specific functional and structural characteristics of individual domains in a new context of artificial fusion proteins may help to reveal the rules of internal and external domain communication. Moreover, this could also help explain the mechanism of such communication and address how the mutual allosteric effect plays a role in a such multi-domain protein system. The simple model system of the two-domain fusion protein investigated in this work consisted of a well-folded PDZ3 domain and an artificially designed small protein domain called Tryptophan Cage (TrpCage). Two fusion proteins with swapped domain order were designed to study their structural and functional features as well as their biophysical properties. The proteins composed of PDZ3 and TrpCage, both identical in amino acid sequence but different in composition (PDZ3-TrpCage, TrpCage-PDZ3), were studied using circualr dichroism (CD) spectrometry, analytical ultracentrifugation, and molecular dynamic simulations. The biophysical analysis uncovered different structural and denaturation properties of both studied proteins, revealing their different unfolding pathways and dynamics.

Mapping of CaM, S100A1 and PIP2-Binding Epitopes in the Intracellular N- and C-Termini of TRPM4

Molecular determinants of the binding of various endogenous modulators to transient receptor potential (TRP) channels are crucial for the understanding of necessary cellular pathways, as well as new paths for rational drug designs. The aim of this study was to characterise interactions between the TRP cation channel subfamily melastatin member 4 (TRPM4) and endogenous intracellular modulators-calcium-binding proteins (calmodulin (CaM) and S100A1) and phosphatidylinositol 4, 5-bisphosphate (PIP2). We have found binding epitopes at the N- and C-termini of TRPM4 shared by CaM, S100A1 and PIP2. The binding affinities of short peptides representing the binding epitopes of N- and C-termini were measured by means of fluorescence anisotropy (FA). The importance of representative basic amino acids and their combinations from both peptides for the binding of endogenous TRPM4 modulators was proved using point alanine-scanning mutagenesis. In silico protein-protein docking of both peptides to CaM and S100A1 and extensive molecular dynamics (MD) simulations enabled the description of key stabilising interactions at the atomic level. Recently solved cryo-Electron Microscopy (EM) structures made it possible to put our findings into the context of the entire TRPM4 channel and to deduce how the binding of these endogenous modulators could allosterically affect the gating of TRPM4. Moreover, both identified binding epitopes seem to be ideally positioned to mediate the involvement of TRPM4 in higher-order hetero-multimeric complexes with important physiological functions.

An intrinsically disordered proteins community for ELIXIR

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) are now recognised as major determinants in cellular regulation. This white paper presents a roadmap for future e-infrastructure developments in the field of IDP research within the ELIXIR framework. The goal of these developments is to drive the creation of high-quality tools and resources to support the identification, analysis and functional characterisation of IDPs. The roadmap is the result of a workshop titled “An intrinsically disordered protein user community proposal for ELIXIR” held at the University of Padua. The workshop, and further consultation with the members of the wider IDP community, identified the key priority areas for the roadmap including the development of standards for data annotation, storage and dissemination; integration of IDP data into the ELIXIR Core Data Resources; and the creation of benchmarking criteria for IDP-related software. Here, we discuss these areas of priority, how they can be implemented in cooperation with the ELIXIR platforms, and their connections to existing ELIXIR Communities and international consortia. The article provides a preliminary blueprint for an IDP Community in ELIXIR and is an appeal to identify and involve new stakeholders.

An ameloblastin C-terminus variant is present in human adipose tissue

Objective

Transcriptional regulatory elements in the ameloblastin (AMBN) promoter indicate that adipogenesis may influence its expression. The objective here was to investigate if AMBN is expressed in adipose tissue, and have a role during differentiation of adipocytes.

Design

AMBN expression was examined in adipose tissue and adipocytes by real-time PCR and ELISA. Distribution of ameloblastin was investigated by immunofluorescence in sections of human subcutaneous adipose tissue. The effect of recombinant proteins resembling AMBN and its processed products on proliferation of primary human pre-adipocytes and murine 3T3-L1 cell lines was measured by [³H]-thymidine incorporation. The effect on adipocyte differentiation was evaluated by the expression profile of the adipogenic markers PPARγ and leptin, and the content of lipids droplets (Oil-Red-O staining).

Results

AMBN was found to be expressed in human adipose tissue, human primary adipocytes, and in 3T3-L1 cells. The C-terminus of the AMBN protein and a 45 bp shorter splice variant was identified in human subcutaneous adipose tissue. The expression of AMBN was found to increase four-fold during differentiation of 3T3-L1 cells. Administration of recombinant AMBN reduced the proliferation, and enhanced the expression of PPARγ and leptin in 3T3-L1 and human pre-adipocytes, respectively.

Conclusions

The AMBN C-terminus variant was identified in adipocytes. This variant may be encoded from a short splice variant. Increased expression of AMBN during adipogenesis and its effect on adipogenic factors suggests that AMBN also has a role in adipocyte development.

Exquisite ligand stereoselectivity of a Drosophila juvenile hormone receptor contrasts with its broad agonist repertoire

The sesquiterpenoid juvenile hormone (JH) is vital to insect development and reproduction. Intracellular JH receptors have recently been established as basic helix-loop-helix transcription factor (bHLH)/PAS proteins in Drosophila melanogaster known as germ cell–expressed (Gce) and its duplicate paralog, methoprene-tolerant (Met). Upon binding JH, Gce/Met activates its target genes. Insects possess multiple native JH homologs whose molecular activities remain unexplored, and diverse synthetic compounds including insecticides exert JH-like effects. How the JH receptor recognizes its ligands is unknown. To determine which structural features define an active JH receptor agonist, we tested several native JHs and their nonnative geometric and optical isomers for the ability to bind the Drosophila JH receptor Gce, to induce Gce-dependent transcription, and to affect the development of the fly. Our results revealed high ligand stereoselectivity of the receptor. The geometry of the JH skeleton, dictated by two stereogenic double bonds, was the most critical feature followed by the presence of an epoxide moiety at a terminal position. The optical isomerism at carbon C11 proved less important even though Gce preferentially bound a natural JH enantiomer. The results of receptor-ligand–binding and cell-based gene activation assays tightly correlated with the ability of different geometric JH isomers to induce gene expression and morphogenetic effects in the developing insects. Molecular modeling supported the requirement for the proper double-bond geometry of JH, which appears to be its major selective mechanism. The strict stereoselectivity of Gce toward the natural hormone contrasts with the high potency of synthetic Gce agonists of disparate chemistries.

Phosphorylation Modulates Ameloblastin Self-assembly and Ca 2+ Binding

Ameloblastin (AMBN), an important component of the self-assembled enamel extra cellular matrix, contains several in silico predicted phosphorylation sites. However, to what extent these sites actually are phosphorylated and the possible effects of such post-translational modifications are still largely unknown. Here we report on in vitro experiments aimed at investigating what sites in AMBN are phosphorylated by casein kinase 2 (CK2) and protein kinase A (PKA) and the impact such phosphorylation has on self-assembly and calcium binding. All predicted sites in AMBN can be phosphorylated by CK2 and/or PKA. The experiments show that phosphorylation, especially in the exon 5 derived part of the molecule, is inversely correlated with AMBN self-assembly. These results support earlier findings suggesting that AMBN self-assembly is mostly dependent on the exon 5 encoded region of the AMBN gene. Phosphorylation was significantly more efficient when the AMBN molecules were in solution and not present as supramolecular assemblies, suggesting that post-translational modification of AMBN must take place before the enamel matrix molecules self-assemble inside the ameloblast cell. Moreover, phosphorylation of exon 5, and the consequent reduction in self-assembly, seem to reduce the calcium binding capacity of AMBN suggesting that post-translational modification of AMBN also can be involved in control of free Ca²⁺ during enamel extra cellular matrix biomineralization. Finally, it is speculated that phosphorylation can provide a functional crossroad for AMBN either to be phosphorylated and act as monomeric signal molecule during early odontogenesis and bone formation, or escape phosphorylation to be subsequently secreted as supramolecular assemblies that partake in enamel matrix structure and mineralization.

A community proposal to integrate proteomics activities in ELIXIR

Computational approaches have been major drivers behind the progress of proteomics in recent years. The aim of this white paper is to provide a framework for integrating computational proteomics into ELIXIR in the near future, and thus to broaden the portfolio of omics technologies supported by this European distributed infrastructure. This white paper is the direct result of a strategy meeting on ‘The Future of Proteomics in ELIXIR’ that took place in March 2017 in Tübingen (Germany), and involved representatives of eleven ELIXIR nodes.

These discussions led to a list of priority areas in computational proteomics that would complement existing activities and close gaps in the portfolio of tools and services offered by ELIXIR so far. We provide some suggestions on how these activities could be integrated into ELIXIR’s existing platforms, and how it could lead to a new ELIXIR use case in proteomics. We also highlight connections to the related field of metabolomics, where similar activities are ongoing. This white paper could thus serve as a starting point for the integration of computational proteomics into ELIXIR. Over the next few months we will be working closely with all stakeholders involved, and in particular with other representatives of the proteomics community, to further refine this paper.

Ameloblastin Peptides Modulates the Osteogenic Capacity of Human Mesenchymal Stem Cells

During amelogenesis the extracellular enamel matrix protein AMBN is quickly processed into 17 kDa (N-terminus) and 23 kDa (C-terminus) fragments. In particular, alternatively spliced regions derived by exon 5/6 within the N-terminus region are known to be critical in biomineralization. Human mesenchymal stem cells (hMSC) also express and secrete AMBN, but it is unclear if this expression has effects on the hMSC themselves. If, as suggested from previous findings, AMBN act as a signaling molecule, such effects could influence hMSC growth and differentiation, as well as promoting the secretion of other signaling proteins like cytokines and chemokines. If AMBN is found to modulate stem cell behavior and fate, it will impact our understanding on how extracellular matrix molecules can have multiple roles during development ontogenesis, mineralization and healing of mesenchymal tissues. Here we show that synthetic peptides representing exon 5 promote hMSC proliferation. Interestingly, this effect is inhibited by the application of a 15 aa peptide representing the alternatively spliced start of exon 6. Both peptides also influence gene expression of RUNX2 and osteocalcin, and promote calcium deposition in cultures, indicating a positive influence on the osteogenic capacity of hMSC. We also show that the full-length AMBN-WT and N-terminus region enhance the secretion of RANTES, IP-10, and IL-8. In contrast, the AMBN C-terminus fragment and the exon 5 deleted AMBN (DelEx5) have no detectable effects on any of the parameters investigated. These findings suggest the signaling effect of AMBN is conveyed by processed products, whereas the effect on proliferation is differentially modulated through alternative splicing during gene expression.

Computational methods for the description of pharmacologically relevant platinum complexes--molecular structure and bond dissociation

Cancer is after cardiovascular disease the most frequent cause of death in Europe. In 28 of 53 countries considered in this area it is already the leading cause of death and expected to gain even more importance until the year 2020. Amongst the large arsenal of different anti-cancer drugs, platinum drugs belong to the first developed anticancer drugs and still have a large impact on cancer therapy. Nevertheless therapy with platinum-anticancer drugs is accompanied by severe adverse effects caused by frequent interactions with the amino acids of different human proteins. Computational chemistry offers methods to study such interactions and even those of not yet synthesized drugs in silico. For such studies a profound knowledge of the prediction quality of various computational methods towards platinum-drug-like complexes is necessary. By this article we are aiming on delivering important accuracy information of the frequently used computational methods. Most important findings are the high performance of the double hybrid functional B2PLYP for the calculation of geometries, even in small basis sets, followed by BP86 and PBE and the still acceptable performance of the semi-empirical Method PM6-D3H4X for extremely large systems. To follow absolute energies of the dissociation process, LPNO-CEPA and B3LYP-D3 can be suggested while SCS-MP2 shows an extremely narrow standard deviation and a low maximum error, which make it an ideal candidate for relative energy calculations in the exploration of reaction mechanisms.

The characterization of a novel S100A1 binding site in the N-terminus of TRPM1

Transient receptor potential melastatin-1 channel (TRPM1) is an important mediator of calcium influx into the cell that is expressed in melanoma and ON-bipolar cells. Similar to other members of the TRP channel family, the intracellular N- and C- terminal domains of TRPM1 are expected to play important roles in the modulation of TRPM1 receptor function. Among the most commonly occurring modulators of TRP channels are the cytoplasmically expressed calcium binding proteins calmodulin and S100 calcium-binding protein A1 (S100A1), but the interaction of TRPM1 with S100A1 has not been described yet. Here, using a combination of biophysical and bioinformatics methods, we have determined that the N-terminal L242-E344 region of TRPM1 is a S100A1 binding domain. We show that formation of the TRPM1/S100A1 complex is calcium-dependent. Moreover, our structural model of the complex explained data obtained from fluorescence spectroscopy measurements revealing that the complex formation is facilitated through interactions of clusters positively charged (K271A, R273A, R274A) and hydrophobic (L263A, V270A, L276A) residues at the N-terminus of TRPM1. Taken together, our data suggest a molecular mechanism for the potential regulation of TRPM1 by S100A1.

Sequence-Specific Recognition of DNA by Proteins: Binding Motifs Discovered Using a Novel Statistical/Computational Analysis

Decades of intensive experimental studies of the recognition of DNA sequences by proteins have provided us with a view of a diverse and complicated world in which few to no features are shared between individual DNA-binding protein families. The originally conceived direct readout of DNA residue sequences by amino acid side chains offers very limited capacity for sequence recognition, while the effects of the dynamic properties of the interacting partners remain difficult to quantify and almost impossible to generalise. In this work we investigated the energetic characteristics of all DNA residue—amino acid side chain combinations in the conformations found at the interaction interface in a very large set of protein—DNA complexes by the means of empirical potential-based calculations. General specificity-defining criteria were derived and utilised to look beyond the binding motifs considered in previous studies. Linking energetic favourability to the observed geometrical preferences, our approach reveals several additional amino acid motifs which can distinguish between individual DNA bases. Our results remained valid in environments with various dielectric properties.

PIP2 and PIP3 interact with N-terminus region of TRPM4 channel

The transient receptor potential melastatin 4 (TRPM4) is a calcium-activated non-selective ion channel broadly expressed in a variety of tissues. Receptor has been identified as a crucial modulator of numerous calcium dependent mechanisms in the cell such as immune response, cardiac conduction, neurotransmission and insulin secretion. It is known that phosphoinositide lipids (PIPs) play a unique role in the regulation of TRP channel function. However the molecular mechanism of this process is still unknown. We characterized the binding site of PIP2 and its structural analogue PIP3 in the E733-W772 proximal region of the TRPM4 N-terminus via biophysical and molecular modeling methods. The specific positions R755 and R767 in this domain were identified as being important for interactions with PIP2/PIP3 ligands. Their mutations caused a partial loss of PIP2/PIP3 binding specificity. The interaction of PIP3 with TRPM4 channels has never been described before. These findings provide new insight into the ligand binding domains of the TRPM4 channel.

Intrinsically disordered enamel matrix protein ameloblastin forms ribbon-like supramolecular structures via an N-terminal segment encoded by exon 5

Tooth enamel, the hardest tissue in the body, is formed by the evolutionarily highly conserved biomineralization process that is controlled by extracellular matrix proteins. The intrinsically disordered matrix protein ameloblastin (AMBN) is the most abundant nonamelogenin protein of the developing enamel and a key element for correct enamel formation. AMBN was suggested to be a cell adhesion molecule that regulates proliferation and differentiation of ameloblasts. Nevertheless, detailed structural and functional studies on AMBN have been substantially limited by the paucity of the purified nondegraded protein. With this study, we have developed a procedure for production of a highly purified form of recombinant human AMBN in quantities that allowed its structural characterization. Using size exclusion chromatography, analytical ultracentrifugation, transmission electron, and atomic force microscopy techniques, we show that AMBN self-associates into ribbon-like supramolecular structures with average widths and thicknesses of 18 and 0.34 nm, respectively. The AMBN ribbons exhibited lengths ranging from tens to hundreds of nm. Deletion analysis and NMR spectroscopy revealed that an N-terminal segment encoded by exon 5 comprises two short independently structured regions and plays a key role in self-assembly of AMBN.

Hypertrophic cardiomyopathy: from mutation to functional analysis of defective protein

Aim

To analyze the genesis of hypertrophic cardiomyopathy on a large cohort of patients from molecular genetics point of view and perform the functional analysis of the 3D molecular model of defective myosin-7 protein in silico.

Methods

The study enrolled 153 patients with diagnosed hypertrophic cardiomyopathy from different parts of the Czech Republic. DNA samples were analyzed for mutations in exons 21 and 22 of the MYH7 gene, which have been associated with high mutation clustering. The 3D model of human myosin-7 was built using the x-ray structure of nucleotide-free scallop myosin S1 as the structural template. We performed de novo structure prediction of mutant and wild type peptides spanning the 769-788 amino acids region of the myosin-7 protein.

Results

The Arg⁸⁷⁰His and Asp⁷⁷⁸Val amino acid alterations were found in 2 unrelated patients with a severe form of hypertrophic cardiomyopathy. The Asp⁷⁷⁸Val variation was chosen for subsequent 3D molecular modeling in silico. The mutation of the Asp by Val not only changes the character of the interaction pattern with other amino acids or ions but Val, being a small hydrophobic amino acid, can also completely change the stability of the region.

Conclusion

Mutation location in the MYH7 gene and changes in amino acid composition may have a crucial negative impact on the outcome of the disease in patients with hypertrophic cardiomyopathy. In addition, a mutation that changes the charge of the amino acid is more likely to affect protein function than a conservative mutation.

Another role of proline: stabilization interactions in proteins and protein complexes concerning proline and tryptophane

Proline-tryptophan complexes derived from experimental structures are investigated by quantum chemical procedures known to properly describe the London dispersion energy. We study two geometrical arrangements: the "L-shaped", stabilized by an H-bond, and the "stacked-like", where the two residues are in parallel orientation without any H-bond. Interestingly, the interaction energies in both cases are comparable and very large ( approximately 7 kcal mol(-1)). The strength of stabilization in the stacked arrangement is rather surprising considering the fact that only one partner has an aromatic character. The interaction energy decomposition using the SAPT method further demonstrates the very important role of dispersion energy in such arrangement. To elucidate the structural features responsible for this unexpectedly large stabilization we examined the role of the nitrogen heteroatom and the importance of the cyclicity of the proline residue. We show that the electrostatic interaction due to the presence of the dipole, caused by the nitrogen heteroatom, contributes largely to the strength of the interaction. Nevertheless, the cyclic arrangement of proline, which allows for the largest amount of dispersive contact with the aromatic partner, also has a notable-effect. Geometry optimizations carried out for the "stacked-like" complexes show that the arrangements derived from protein structure are close to their gas phase optimum geometry, suggesting that the environment has only a minor effect on the geometry of the interaction. We conclude that the strength of proline non-covalent interactions, combined with this residue's rigidity, might be the explanation for its prominent role in protein stabilization and recognition processes.

MedusaScore: an accurate force field-based scoring function for virtual drug screening

Virtual screening is becoming an important tool for drug discovery. However, the application of virtual screening has been limited by the lack of accurate scoring functions. Here, we present a novel scoring function, MedusaScore, for evaluating protein-ligand binding. MedusaScore is based on models of physical interactions that include van der Waals, solvation, and hydrogen bonding energies. To ensure the best transferability of the scoring function, we do not use any protein-ligand experimental data for parameter training. We then test the MedusaScore for docking decoy recognition and binding affinity prediction and find superior performance compared to other widely used scoring functions. Statistical analysis indicates that one source of inaccuracy of MedusaScore may arise from the unaccounted entropic loss upon ligand binding, which suggests avenues of approach for further MedusaScore improvement.

Chemical compound navigator: A web-based chem-BLAST, chemical taxonomy-based search engine for browsing compounds

A novel technique to annotate, query, and analyze chemical compounds has been developed and is illustrated by using the inhibitor data on HIV protease-inhibitor complexes. In this method, all chemical compounds are annotated in terms of standard chemical structural fragments. These standard fragments are defined by using criteria, such as chemical classification; structural, chemical, or functional groups; and commercial, scientific or common names or synonyms. These fragments are then organized into a data tree based on their chemical substructures. Search engines have been developed to use this data tree to enable query on inhibitors of HIV protease (http://xpdb.nist.gov/hivsdb/hivsdb.html). These search engines use a new novel technique, Chemical Block Layered Alignment of Substructure Technique (Chem-BLAST) to search on the fragments of an inhibitor to look for its chemical structural neighbors. This novel technique to annotate and query compounds lays the foundation for the use of the Semantic Web concept on chemical compounds to allow end users to group, sort, and search structural neighbors accurately and efficiently. During annotation, it enables the attachment of "meaning" (i.e., semantics) to data in a manner that far exceeds the current practice of associating "metadata" with data by creating a knowledge base (or ontology) associated with compounds. Intended users of the technique are the research community and pharmaceutical industry, for which it will provide a new tool to better identify novel chemical structural neighbors to aid drug discovery.

Application of InChI to curate, index, and query 3-D structures

The HIV structural database (HIVSDB) is a comprehensive collection of the structures of HIV protease, both of unliganded enzyme and of its inhibitor complexes. It contains abstracts and crystallographic data such as inhibitor and protein coordinates for 248 data sets, of which only 141 are from the Protein Data Bank (PDB). Efficient annotation, indexing, and querying of the inhibitor data is crucial for their effective use for technological and industrial applications. The application of IUPAC International Chemical Identifier (InChI) to index, curate, and query inhibitor structures HIVSDB is described.

Inhibitors of HIV-1 protease: a major success of structure-assisted drug design

Retroviral protease (PR) from the human immunodeficiency virus type 1 (HIV-1) was identified over a decade ago as a potential target for structure-based drug design. This effort was very successful. Four drugs are already approved, and others are undergoing clinical trials. The techniques utilized in this remarkable example of structure-assisted drug design included crystallography, NMR, computational studies, and advanced chemical synthesis. The development of these drugs is discussed in detail. Other approaches to designing HIV-1 PR inhibitors, based on the concepts of symmetry and on the replacement of a water molecule that had been found tetrahedrally coordinated between the enzyme and the inhibitors, are also discussed. The emergence of drug-induced mutations of HIV-1 PR leads to rapid loss of potency of the existing drugs and to the need to continue the development process. The structural basis of drug resistance and the ways of overcoming this phenomenon are mentioned.