Structure of the St. Louis encephalitis virus postfusion envelope trimer

St. Louis encephalitis virus (SLEV) is a mosquito-borne flavivirus responsible for several human encephalitis outbreaks over the last 80 years. Mature flavivirus virions are coated with dimeric envelope (E) proteins that mediate attachment and fusion with host cells. E is a class II fusion protein, the hallmark of which is a distinct dimer-to-trimer rearrangement that occurs upon endosomal acidification and insertion of hydrophobic fusion peptides into the endosomal membrane. Herein, we report the crystal structure of SLEV E in the posfusion trimer conformation. The structure revealed specific features that differentiate SLEV E from trimers of related flavi- and alphaviruses. SLEV E fusion loops have distinct intermediate spacing such that they are positioned further apart than previously observed in flaviviruses but closer together than Semliki Forest virus, an alphavirus. Domains II and III (DII and DIII) of SLEV E also adopt different angles relative to DI, which suggests that the DI-DII joint may accommodate spheroidal motions. However, trimer interfaces are well conserved among flaviviruses, so it is likely the differences observed represent structural features specific to SLEV function. Analysis of surface potentials revealed a basic platform underneath flavivirus fusion loops that may interact with the anionic lipid head groups found in membranes. Taken together, these results highlight variations in E structure and assembly that may direct virus-specific interactions with host determinants to influence pathogenesis

Identification of a better Homo sapiens Class II HDAC inhibitor through binding energy calculations and descriptor analysis

Human papillomaviruses (HPVs) are the most common on sexually transmitted viruses in the world. HPVs are responsible for a large spectrum of deseases, both benign and malignant. The certain types of HPV are involved in the development of cervical cancer. In attemps to find additional drugs in the treatment of cervical cancer, inhibitors of the histone deacetylases (HDAC) have received much attention due to their low cytotoxic profiles and the E6/E7 oncogene function of human papilomavirus can be completely by passed by HDAC inhibition. The histone deacetylase inhibitors can induce growth arrest, differentiation and apoptosis of cancer cells. HDAC class I and class II are considered the main targets for cancer. Therefore, the six HDACs class II was modeled and about two inhibitors (SAHA and TSA) were docked using AutoDock4.2, to each of the inhibitor in order to identify the pharmacological properties. Based on the results of docking, SAHA and TSA were able to bind with zinc ion in HDACs models as a drug target. SAHA was satisfied almost all the properties i.e., binding affinity, the Drug-Likeness value and Drug Score with 70% oral bioavailability and the carbonyl group of these compound fits well into the active site of the target where the zinc is present. Hence, SAHA could be developed as potential inhibitors of class II HDACs and valuable cervical cancer drug candidate

The Isoelectric Region of Proteins: A Systematic Analysis

Background: Binding of proteins in ion exchange chromatography is dominated by electrostatic interactions and can be tuned by adjusting pH and ionic strength of the solvent. Therefore, the isoelectric region (IER), the pH region of almost zero charge near the pI, has been used to predict the binding properties of proteins. Principal findings: Usually the IER is small and binding and elution is carried out at pH values near to the pI. However, some proteins with an extended IER have been shown to bind and elute far away from its pI. To analyze factors that mediate the size of the IER and to identify proteins with an extended IER, two protein families consisting of more than 7000 proteins were systematically investigated. Most proteins were found to have a small IER and thus are expected to bind or elute near to their pI, while only a small fraction of less than 2 % had a large IER. Conclusions: Only four factors, the number of histidines, the pI, the number of titratable amino acids and the ratio of acidic to basic residues, are sufficient to reliably classify proteins by their IER based on their sequence only, and thus to predict their binding and elution behaviour in ion exchange chromatography

Computational Identification of Uncharacterized Cruzain Binding Sites

Chagas disease, caused by the unicellular parasite Trypanosoma cruzi, claims 50,000 lives annually and is the leading cause of infectious myocarditis in the world. As current antichagastic therapies like nifurtimox and benznidazole are highly toxic, ineffective at parasite eradication, and subject to increasing resistance, novel therapeutics are urgently needed. Cruzain, the major cysteine protease of Trypanosoma cruzi, is one attractive drug target. In the current work, molecular dynamics simulations and a sequence alignment of a non-redundant, unbiased set of peptidase C1 family members are used to identify uncharacterized cruzain binding sites. The two sites identified may serve as targets for future pharmacological intervention

Minimum Free Energy Path of Ligand-Induced Transition in Adenylate Kinase

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme

Hydrophilicity Matching – A Potential Prerequisite for the Formation of Protein-Protein Complexes in the Cell

A binding event between two proteins typically consists of a diffusional search of binding partners for one another, followed by a specific recognition of the compatible binding sites resulting in the formation of the complex. However, it is unclear how binding partners find each other in the context of the crowded, constantly fluctuating, and interaction-rich cellular environment. Here we examine the non-specific component of protein-protein interactions, which refers to those physicochemical properties of the binding partners that are independent of the exact details of their binding sites, but which can affect their localization or diffusional search for one another. We show that, for a large set of high-resolution experimental 3D structures of binary, transient protein complexes taken from the DOCKGROUND database, the binding partners display a surprising, statistically significant similarity in terms of their total hydration free energies normalized by a size-dependent variable. We hypothesize that colocalization of binding partners, even within individual cellular compartments such as the cytoplasm, may be influenced by their relative hydrophilicity, potentially in response to local hydrophilic gradients

Probing of Exosites Leads to Novel Inhibitor Scaffolds of HCV NS3/4A Proteinase

Hepatitis C is a treatment-resistant disease affecting millions of people worldwide. The hepatitis C virus (HCV) genome is a single-stranded RNA molecule. After infection of the host cell, viral RNA is translated into a polyprotein that is cleaved by host and viral proteinases into functional, structural and non-structural, viral proteins. Cleavage of the polyprotein involves the viral NS3/4A proteinase, a proven drug target. HCV mutates as it replicates and, as a result, multiple emerging quasispecies become rapidly resistant to anti-virals, including NS3/4A inhibitors.To circumvent drug resistance and complement the existing anti-virals, NS3/4A inhibitors, which are additional and distinct from the FDA-approved telaprevir and boceprevir α-ketoamide inhibitors, are required. To test potential new avenues for inhibitor development, we have probed several distinct exosites of NS3/4A which are either outside of or partially overlapping with the active site groove of the proteinase. For this purpose, we employed virtual ligand screening using the 275,000 compound library of the Developmental Therapeutics Program (NCI/NIH) and the X-ray crystal structure of NS3/4A as a ligand source and a target, respectively. As a result, we identified several novel, previously uncharacterized, nanomolar range inhibitory scaffolds, which suppressed of the NS3/4A activity in vitro and replication of a sub-genomic HCV RNA replicon with a luciferase reporter in human hepatocarcinoma cells. The binding sites of these novel inhibitors do not significantly overlap with those of α-ketoamides. As a result, the most common resistant mutations, including V36M, R155K, A156T, D168A and V170A, did not considerably diminish the inhibitory potency of certain novel inhibitor scaffolds we identified.Overall, the further optimization of both the in silico strategy and software platform we developed and lead compounds we identified may lead to advances in novel anti-virals

Synergy between trastuzumab and pertuzumab for human epidermal growth factor 2 (Her2) from colocalization: an in silico based mechanism

10.1186/bcr2888Breast Cancer Research133-BCRR

Prediction of Co-Receptor Usage of HIV-1 from Genotype

Human Immunodeficiency Virus 1 uses for entry into host cells a receptor (CD4) and one of two co-receptors (CCR5 or CXCR4). Recently, a new class of antiretroviral drugs has entered clinical practice that specifically bind to the co-receptor CCR5, and thus inhibit virus entry. Accurate prediction of the co-receptor used by the virus in the patient is important as it allows for personalized selection of effective drugs and prognosis of disease progression. We have investigated whether it is possible to predict co-receptor usage accurately by analyzing the amino acid sequence of the main determinant of co-receptor usage, i.e., the third variable loop V3 of the gp120 protein. We developed a two-level machine learning approach that in the first level considers two different properties important for protein-protein binding derived from structural models of V3 and V3 sequences. The second level combines the two predictions of the first level. The two-level method predicts usage of CXCR4 co-receptor for new V3 sequences within seconds, with an area under the ROC curve of 0.937±0.004. Moreover, it is relatively robust against insertions and deletions, which frequently occur in V3. The approach could help clinicians to find optimal personalized treatments, and it offers new insights into the molecular basis of co-receptor usage. For instance, it quantifies the importance for co-receptor usage of a pocket that probably is responsible for binding sulfated tyrosine

Cytoskeletal Signaling: Is Memory Encoded in Microtubule Lattices by CaMKII Phosphorylation?

Memory is attributed to strengthened synaptic connections among particular brain neurons, yet synaptic membrane components are transient, whereas memories can endure. This suggests synaptic information is encoded and ‘hard-wired’ elsewhere, e.g. at molecular levels within the post-synaptic neuron. In long-term potentiation (LTP), a cellular and molecular model for memory, post-synaptic calcium ion (Ca2+) flux activates the hexagonal Ca2+-calmodulin dependent kinase II (CaMKII), a dodacameric holoenzyme containing 2 hexagonal sets of 6 kinase domains. Each kinase domain can either phosphorylate substrate proteins, or not (i.e. encoding one bit). Thus each set of extended CaMKII kinases can potentially encode synaptic Ca2+ information via phosphorylation as ordered arrays of binary ‘bits’. Candidate sites for CaMKII phosphorylation-encoded molecular memory include microtubules (MTs), cylindrical organelles whose surfaces represent a regular lattice with a pattern of hexagonal polymers of the protein tubulin. Using molecular mechanics modeling and electrostatic profiling, we find that spatial dimensions and geometry of the extended CaMKII kinase domains precisely match those of MT hexagonal lattices. This suggests sets of six CaMKII kinase domains phosphorylate hexagonal MT lattice neighborhoods collectively, e.g. conveying synaptic information as ordered arrays of six “bits”, and thus “bytes”, with 64 to 5,281 possible bit states per CaMKII-MT byte. Signaling and encoding in MTs and other cytoskeletal structures offer rapid, robust solid-state information processing which may reflect a general code for MT-based memory and information processing within neurons and other eukaryotic cells