Mapping the druggable allosteric space of G-protein coupled receptors: a fragment-based molecular dynamics approach.

To address the problem of specificity in G-protein coupled receptor (GPCR) drug discovery, there has been tremendous recent interest in allosteric drugs that bind at sites topographically distinct from the orthosteric site. Unfortunately, structure-based drug design of allosteric GPCR ligands has been frustrated by the paucity of structural data for allosteric binding sites, making a strong case for predictive computational methods. In this work, we map the surfaces of the beta1 (beta1AR) and beta2 (beta2AR) adrenergic receptor structures to detect a series of five potentially druggable allosteric sites. We employ the FTMAP algorithm to identify 'hot spots' with affinity for a variety of organic probe molecules corresponding to drug fragments. Our work is distinguished by an ensemble-based approach, whereby we map diverse receptor conformations taken from molecular dynamics (MD) simulations totaling approximately 0.5 micros. Our results reveal distinct pockets formed at both solvent-exposed and lipid-exposed cavities, which we interpret in light of experimental data and which may constitute novel targets for GPCR drug discovery. This mapping data can now serve to drive a combination of fragment-based and virtual screening approaches for the discovery of small molecules that bind at these sites and which may offer highly selective therapies

Efficient and Accurate Modeling of Conformational Transitions in Proteins: The Case of c-Src Kinase

The theoretical computational modeling of large conformational transitions occurring in biomolecules still represents a challenge. Here, we present an accurate "in silico" description of the activation and deactivation mechanisms of human c-Src kinases, a fundamental process regulating several crucial cell functions. Our results clearly show that by applying an efficient and automated algorithm able to drive the molecular dynamics (MD) sampling along the pathway between the two c-Src conformational states - the active state and the inactive state - it is possible to accurately describe, at reduced computational costs, the molecular mechanism underlying these large conformational rearrangements. This procedure, combining the MD simulations with the sampling along the well-defined principal motions connecting the two conformational states, allows to provide a description well beyond the present computational limits, and it is easily applicable to different systems where the structures of both the initial and final states are known

Label-free optical detection of single enzyme-reactant reactions and associated conformational changes

Monitoring the kinetics and conformational dynamics of single enzymes is crucial in order to better understand their biological functions as these motions and structural dynamics are usually unsynchronized among the molecules. Detecting the enzyme-reactant interactions and associated conformational changes of the enzyme on a single molecule basis, however, remain as a challenge with established optical techniques due to the commonly required labeling of the reactants or the enzyme itself. The labeling process is usually non-trivial and the labels themselves might skew the physical properties of the enzyme. Here we demonstrate an optical, label-free method capable of observing enzymatic interactions and the associated conformational changes on the single molecule level. We monitor polymerase/DNA interactions via the strong near-field enhancement provided by plasmonic nanorods resonantly coupled to whispering gallery modes in microcavities. Specifically, we employ two different recognition schemes: one in which the kinetics of polymerase/DNA interactions are probed in the vicinity of DNA-functionalized nanorods, and the other in which these interactions are probed via the magnitude of conformational changes in the polymerase molecules immobilized on nanorods. In both approaches we find that low and high polymerase activities can be clearly discerned via their characteristic signal amplitude and signal length distributions. Furthermore, the thermodynamic study of the monitored interactions suggests the occurrence of DNA polymerization. This work constitutes a proof-of-concept study of enzymatic activities via plasmonically enhanced microcavities and establishes an alternative and label-free method capable of investigating structural changes in single molecules

Morphine activates neuroinflammation in a manner parallel to endotoxin

Opioids create a neuroinflammatory response within the CNS, compromising opioid-induced analgesia and contributing to various unwanted actions. How this occurs is unknown but has been assumed to be via classic opioid receptors. Herein, we provide direct evidence that morphine creates neuroinflammation via the activation of an innate immune receptor and not via classic opioid receptors. We demonstrate that morphine binds to an accessory protein of Toll-like receptor 4 (TLR4), myeloid differentiation protein 2 (MD-2), thereby inducing TLR4 oligomerization and triggering proinflammation. Small-molecule inhibitors, RNA interference, and genetic knockout validate the TLR4/MD-2 complex as a feasible target for beneficially modifying morphine actions. Disrupting TLR4/MD-2 protein–protein association potentiated morphine analgesia in vivo and abolished morphine-induced proinflammation in vitro, the latter demonstrating that morphine-induced proinflammation only depends on TLR4, despite the presence of opioid receptors. These results provide an exciting, nonconventional avenue to improving the clinical efficacy of opioids.Xiaohui Wang, Lisa C. Loram, Khara Ramos, Armando J. de Jesus, Jacob Thomas, Kui Cheng, Anireddy Reddy, Andrew A. Somogyi, Mark R. Hutchinson, Linda R. Watkins and Hang Yi

Mechanism of glycan receptor recognition and specificity switch for avian, swine, and human adapted influenza virus hemagglutinins: a molecular dynamics perspective.

Hemagglutinins (HA's) from duck, swine, and human influenza viruses have previously been shown to prefer avian and human glycan receptor analogues with distinct topological profiles, pentasaccharides LSTa (alpha-2,3 linkage) and LSTc (alpha-2,6 linkage), in comparative molecular dynamics studies. On the basis of detailed analyses of the dynamic motions of the receptor binding domains (RBDs) and interaction energy profiles with individual glycan residues, we have identified approximately 30 residue positions in the RBD that present distinct profiles with the receptor analogues. Glycan binding constrained the conformational space sampling by the HA. Electrostatic steering appeared to play a key role in glycan binding specificity. The complex dynamic behaviors of the major SSE and trimeric interfaces with or without bound glycans suggested that networks of interactions might account for species specificity in these low affinity and high avidity (multivalent) interactions between different HA and glycans. Contact frequency, energetic decomposition, and H-bond analyses revealed species-specific differences in HA-glycan interaction profiles, not readily discernible from crystal structures alone. Interaction energy profiles indicated that mutation events at the set of residues such as 145, 156, 158, and 222 would favor human or avian receptor analogues, often through interactions with distal asialo-residues. These results correlate well with existing experimental evidence, and suggest new opportunities for simulation-based vaccine and drug development

Molecular modeling of an antigenic complex between a viral peptide and a class I major histocompatibility glycoprotein

Computer simulation of the conformations of short antigenic peptides (&lo residues) either free or bound to their receptor, the major histocompatibility complex (MHC)- encoded glycoprotein H-2 Ld, was employed to explain experimentally determined differences in the antigenic activities within a set of related peptides. Starting for each sequence from the most probable conformations disclosed by a pattern-recognition technique, several energyminimized structures were subjected to molecular dynamics simulations (MD) either in vacuo or solvated by water molecules. Notably, antigenic potencies were found to correlate to the peptides propensity to form and maintain an overall a-helical conformation through regular i,i + 4 hydrogen bonds. Accordingly, less active or inactive peptides showed a strong tendency to form i,i+3 hydrogen bonds at their Nterminal end. Experimental data documented that the C-terminal residue is critical for interaction of the peptide with H-2 Ld. This finding could be satisfactorily explained by a 3-D Q.S.A.R. analysis postulating interactions between ligand and receptor by hydrophobic forces. A 3-D model is proposed for the complex between a high-affinity nonapeptide and the H- 2 Ld receptor. First, the H-2 Ld molecule was built from X-ray coordinates of two homologous proteins: HLA-A2 and HLA-Aw68, energyminimized and studied by MD simulations. With HLA-A2 as template, the only realistic simulation was achieved for a solvated model with minor deviations of the MD mean structure from the X-ray conformation. Water simulation of the H-2 Ld protein in complex with the antigenic nonapeptide was then achieved with the template- derived optimal parameters. The bound peptide retains mainly its a-helical conformation and binds to hydrophobic residues of H-2 Ld that correspond to highly polymorphic positions of MHC proteins. The orientation of the nonapeptide in the binding cleft is in accordance with the experimentally determined distribution of its MHC receptor-binding residues (agretope residues). Thus, computer simulation was successfully employed to explain functional data and predicts a-helical conformation for the bound peptid

Deciphering Chemomechanical Couplings in Proteins Using Microsecond-level Molecular Dynamics Simulations

All-atom molecular dynamics (MD) simulations combine the high temporal resolution of experimental methods like smFRET and spatial resolution of methods like x-ray crystallography, to provide a detailed dynamic picture of biomolecular processes. Here, microsecond-level atomistic MD simulations have been used to characterize chemomechanical couplings in human fibroblast growth factor 1 (hFGF1) and the spike proteins of SARS CoV-1 and SARS-CoV-2. hFGF1 is a globular signaling protein that is involved in several physiological processes ranging from cell proliferation to wound healing. Experimental studies have previously described the low proteolytic and thermal stability of hFGF1, in addition to the stabilizing role of heparin. Here, a conformational change in the hFGF1 heparin-binding pocket that occurs only when heparin is absent, is described for the first time. Comparisons with experimental data indicate that this conformational transition is implicated in the low thermal stability of hFGF1. Unique electrostatic interactions that contribute to heparin-mediated stabilization are also described. This work also describes a novel binding affinity estimation approach involving restrained umbrella sampling simulations. The absolute binding affinity for the hFGF1-heparin interaction determined using this approach is in very good agreement with data from isothermal titration calorimetry (ITC) experiments. This binding affinity study revealed that restraining ligand orientation is essential for effective sampling along a protein-ligand distance collective variable.The differential dynamic behavior of the SARS-CoV-1 and CoV-2 spike proteins is also described in this work. Spike protein activation is the first step in the “effective binding” process leading to interaction with the human ACE2 receptor. This study shows that the active form of the CoV-1 spike protein is less stable than that of the CoV-2 spike protein and that the energy barriers associated with activation and inactivation are higher in CoV-2. A “pseudo-inactive” state of the CoV-1 spike protein is described for the first time, wherein the N-terminal domain (NTD) interacts with the receptor-binding domain (RBD). This highlights the potential role of the NTD in spike protein inactivation. The relatively slower kinetics of spike protein activation and inactivation in CoV-2 indicate that it might spend more time bound to the ACE2 receptor than CoV-1, which in turn might provide an explanation for the higher transmissibility of CoV-2

Polycation-π Interactions Are a Driving Force for Molecular Recognition by an Intrinsically Disordered Oncoprotein Family

Molecular recognition by intrinsically disordered proteins (IDPs) commonly involves specific localized contacts and target-induced disorder to order transitions. However, some IDPs remain disordered in the bound state, a phenomenon coined "fuzziness", often characterized by IDP polyvalency, sequence-insensitivity and a dynamic ensemble of disordered bound-state conformations. Besides the above general features, specific biophysical models for fuzzy interactions are mostly lacking. The transcriptional activation domain of the Ewing's Sarcoma oncoprotein family (EAD) is an IDP that exhibits many features of fuzziness, with multiple EAD aromatic side chains driving molecular recognition. Considering the prevalent role of cation-π interactions at various protein-protein interfaces, we hypothesized that EAD-target binding involves polycation- π contacts between a disordered EAD and basic residues on the target. Herein we evaluated the polycation-π hypothesis via functional and theoretical interrogation of EAD variants. The experimental effects of a range of EAD sequence variations, including aromatic number, aromatic density and charge perturbations, all support the cation-π model. Moreover, the activity trends observed are well captured by a coarse-grained EAD chain model and a corresponding analytical model based on interaction between EAD aromatics and surface cations of a generic globular target. EAD-target binding, in the context of pathological Ewing's Sarcoma oncoproteins, is thus seen to be driven by a balance between EAD conformational entropy and favorable EAD-target cation-π contacts. Such a highly versatile mode of molecular recognition offers a general conceptual framework for promiscuous target recognition by polyvalent IDPs. © 2013 Song et al

Functional Diversity and Structural Disorder in the Human Ubiquitination Pathway

The ubiquitin-proteasome system plays a central role in cellular regulation and protein quality control (PQC). The system is built as a pyramid of increasing complexity, with two E1 (ubiquitin activating), few dozen E2 (ubiquitin conjugating) and several hundred E3 (ubiquitin ligase) enzymes. By collecting and analyzing E3 sequences from the KEGG BRITE database and literature, we assembled a coherent dataset of 563 human E3s and analyzed their various physical features. We found an increase in structural disorder of the system with multiple disorder predictors (IUPred - E1: 5.97%, E2: 17.74%, E3: 20.03%). E3s that can bind E2 and substrate simultaneously (single subunit E3, ssE3) have significantly higher disorder (22.98%) than E3s in which E2 binding (multi RING-finger, mRF, 0.62%), scaffolding (6.01%) and substrate binding (adaptor/substrate recognition subunits, 17.33%) functions are separated. In ssE3s, the disorder was localized in the substrate/adaptor binding domains, whereas the E2-binding RING/HECT-domains were structured. To demonstrate the involvement of disorder in E3 function, we applied normal modes and molecular dynamics analyses to show how a disordered and highly flexible linker in human CBL (an E3 that acts as a regulator of several tyrosine kinase-mediated signalling pathways) facilitates long-range conformational changes bringing substrate and E2-binding domains towards each other and thus assisting in ubiquitin transfer. E3s with multiple interaction partners (as evidenced by data in STRING) also possess elevated levels of disorder (hubs, 22.90% vs. non-hubs, 18.36%). Furthermore, a search in PDB uncovered 21 distinct human E3 interactions, in 7 of which the disordered region of E3s undergoes induced folding (or mutual induced folding) in the presence of the partner. In conclusion, our data highlights the primary role of structural disorder in the functions of E3 ligases that manifests itself in the substrate/adaptor binding functions as well as the mechanism of ubiquitin transfer by long-range conformational transitions. © 2013 Bhowmick et al