Functional dynamics of the folded ankyrin repeats of I kappa B alpha revealed by nuclear magnetic resonance.

Inhibition of nuclear factor kappaB (NF-kappaB) is mainly accomplished by IkappaB alpha, which consists of a signal response sequence at the N-terminus, a six-ankyrin repeat domain (ARD) that binds NF-kappaB, and a C-terminal PEST sequence. Previous studies with the ARD revealed that the fifth and sixth repeats are only partially folded in the absence of NF-kappaB. Here we report NMR studies of a truncated version of IkappaB alpha, containing only the first four ankyrin repeats, IkappaB alpha(67-206). This four-repeat segment is well-structured in the free state, enabling full resonance assignments to be made. H-D exchange, backbone dynamics, and residual dipolar coupling (RDC) experiments reveal regions of flexibility. In addition, regions consistent with the presence of micro- to millisecond motions occur periodically throughout the repeat structure. Comparison of the RDCs with the crystal structure gave only moderate agreement, but an ensemble of structures generated by accelerated molecular dynamics gave much better agreement with the measured RDCs. The regions showing flexibility correspond to those implicated in entropic compensation for the loss of flexibility in ankyrin repeats 5 and 6 upon binding to NF-kappaB. The regions showing micro- to millisecond motions in the free protein are the ends of the beta-hairpins that directly interact with NF-kappaB in the complex

Crystal structures of trypanosoma brucei oligopeptidase B broaden the paradigm of catalytic regulation in prolyl oligopeptidase family enzymes

Oligopeptidase B cleaves after basic amino acids in peptides up to 30 residues. As a virulence factor in bacteria and trypanosomatid pathogens that is absent in higher eukaryotes, this is a promising drug target. Here we present ligand-free open state and inhibitor-bound closed state crystal structures of oligopeptidase B from Trypanosoma brucei, the causative agent of African sleeping sickness. These (and related) structures show the importance of structural dynamics, governed by a fine enthalpic and entropic balance, in substrate size selectivity and catalysis. Peptides over 30 residues cannot fit the enzyme cavity, preventing the complete domain closure required for a key propeller Asp/Glu to fix the catalytic His and Arg in the catalytically competent conformation. This size exclusion mechanism protects larger peptides and proteins from degradation. Similar bacterial prolyl endopeptidase and archael acylaminoacyl peptidase structures demonstrate this mechanism is conserved among oligopeptidase family enzymes across all three domains of life

Experimental Validation of the Predicted Binding Site of Escherichia coli K1 Outer Membrane Protein A to Human Brain Microvascular Endothelial Cells: Identification of Critical Mutations That Prevent E. coli Meningitis

Escherichia coli K1, the most common cause of meningitis in neonates, has been shown to interact with GlcNAc1–4GlcNAc epitopes of Ecgp96 on human brain microvascular endothelial cells (HBMECs) via OmpA (outer membrane protein A). However, the precise domains of extracellular loops of OmpA interacting with the chitobiose epitopes have not been elucidated. We report the loop-barrel model of these OmpA interactions with the carbohydrate moieties of Ecgp96 predicted from molecular modeling. To test this model experimentally, we generated E. coli K1 strains expressing OmpA with mutations of residues predicted to be critical for interaction with the HBMEC and tested E. coli invasion efficiency. For these same mutations, we predicted the interaction free energies (including explicit calculation of the entropy) from molecular dynamics (MD), finding excellent correlation (R^2 = 90%) with experimental invasion efficiency. Particularly important is that mutating specific residues in loops 1, 2, and 4 to alanines resulted in significant inhibition of E. coli K1 invasion in HBMECs, which is consistent with the complete lack of binding found in the MD simulations for these two cases. These studies suggest that inhibition of the interactions of these residues of Loop 1, 2, and 4 with Ecgp96 could provide a therapeutic strategy to prevent neonatal meningitis due to E. coli K1

Mining electron density for functionally relevant protein polysterism in crystal structures.

This review focuses on conceptual and methodological advances in our understanding and characterization of the conformational heterogeneity of proteins. Focusing on X-ray crystallography, we describe how polysterism, the interconversion of pre-existing conformational substates, has traditionally been analyzed by comparing independent crystal structures or multiple chains within a single crystal asymmetric unit. In contrast, recent studies have focused on mining electron density maps to reveal previously 'hidden' minor conformational substates. Functional tests of the importance of minor states suggest that evolutionary selection shapes the entire conformational landscape, including uniquely configured conformational substates, the relative distribution of these substates, and the speed at which the protein can interconvert between them. An increased focus on polysterism may shape the way protein structure and function is studied in the coming years

Conformational Entropy in PEPCK Catalysis: Dynamic Motions Critical to Function

Enzymes increase the rate at which chemical reactions occur. How they achieve this rate increase has been an area of intense research for many years. Multiple studies have shown that dynamic motions of individual protein segments can play a role in enzyme function. Our enzyme of interest, phosphoenolpyruvate carboxykinase (PEPCK), has multiple dynamic motions that work in concert during the catalytic cycle. One of these mobile elements, termed the omega-loop due to its architecture, acts like a lid that closes over the active site and has been hypothesized to serve multiple functions during the catalytic cycle. Furthermore, the act of opening and closing the lid domain is thought to be a delicate free energy balance, meaning the lid pays an energetic penalty to close, which is offset by the energetic contribution from substrate binding. To investigate the role of the loop, specifically the open/closed transition, and to probe for unappreciated roles, multiple mutant enzymes were created. In the first set of mutations the lid region was removed to investigate other potential roles for the omega-loop. In the second set two different amino acid residues (Ala-467 and Glu-89) were mutated to either increase the energetic penalty of lid closure or decrease the energetic contribution from substrate binding. The resultant mutants were characterized via structure/function experiments. The data revealed two new roles for the omega-loop lid, confirmed the delicate free energy balance, and shed light on the energetic pathway between the active site and the lid. Taking into account the data from the WT and all seven mutant PEPCK enzymes, a revised model for the role of the omega-loop lid during catalysis is proposed: The role of the lid during the catalytic cycle is to correctly position the substrates in the active site, stabilize the R-loop and N/C-terminal lobe closure, and sequester and protect the reaction intermediate

Energetic and Dynamic Analysis of Inhibitor Binding to Drug-Resistant HIV-1 Proteases: A Dissertation

HIV-1 protease is a very important drug target for AIDS therapy. Nine protease inhibitors have been proved by FDA and used in AIDS treatment. Due to the high replication rate and the lack of fidelity of the HIV-1 reverse transcriptase, HIV-1 virus developed various drug-resistant variants. Although experimental methods such as crystallography and isothermal titration calorimetry provide structural and thermodynamic data on drug-resistant variants, they are unable to discern the mechanism by which the mutations confer resistance to inhibitors. Understanding the drug-resistance mechanism is crucial for developing new inhibitors more tolerant to the drug-resistant mutations. Computational methods such as free energy calculations and molecular dynamic simulations can provide insights to the drug resistance mechanism at an atomic level. In this thesis, I have focused on the elucidation of the energetic and dynamics of key drug-resistant variants of HIV-1 protease. Two multi-drug resistant variants, in comparison with wild-type HIV-1 protease were used for the comparisons: Flap+ (L10I, G48V, I54V, and V82A) which contains a combination of flap and active site mutations and ACT (V82T, I84V) that only contains active site mutations. In Chapter II, I applied free energy simulations and decomposition methods to study the differential mechanism of resistance to the two variants, Flap+ and ACT, to the recently FDA-approved protease inhibitor darunavir (DRV). In this study, the absolute and relative binding free energies of DRV with wild-type protease and the two protease variants were calculated with MM-PB/GBSA and thermodynamic integration methods, respectively. And the predicted results are in good agreement with the ITC experimental results. Free energy decomposition elucidates the mutations alter not only its own interaction with DRV but also other residues by changing the geometry of binding pocket. And the VdW interactions between the bis-THF group of DRV is predominant even in the drug-resistant variants. At the end of this chapter, I offer suggestions on developing new inhibitors that are based on DRV but might be less susceptible to drug-resistant mutations. In Chapter III, 20-ns MD simulations of the apo wildtype protease and the apo drug-resistant protease variant Flap+ are analyzed and compared. In these studies, these mutations have been found to decrease the protease flexibility in the apo form but increase the mobility when the protease is binding with inhibitor. In Chapter IV, more details of the free energy simulation and decomposition are discussed. NMR relaxation experiments were set up as a control for the MD simulation study of the dynamics of the Flap+ variant. The difficulty of finishing the NMR experiment is discussed and the solution and some preliminary results are shown. In summary, the scope of this thesis was to use computational methods to study drug-resistant protease variants’ thermodynamic and dynamic properties to illuminate the mechanism of protease drug resistance. This knowledge will contribute to rational design of new protease inhibitors which bind more tightly to the protease and hinder the development of drug-resistant mutations

Modulating Ligand Dissociation through Methyl Isomerism in Accessory Sites: Binding of Retinol to Cellular Carriers

Due to the poor aqueous solubility of retinoids, evolution has tuned their binding to cellular proteins to address specialized physiological roles by modulating uptake, storage, and delivery to specific targets. With the aim to disentangle the structure-function relationships in these proteins and disclose clues for engineering selective carriers, the binding mechanism of the two most abundant retinol-binding isoforms was explored by using enhanced sampling molecular dynamics simulations and surface plasmon resonance. The distinctive dynamics of the entry portal site in the holo species was crucial to modulate retinol dissociation. Remarkably, this process is controlled at large extent by the replacement of Ile by Leu in the two isoforms, thus suggesting that a fine control of ligand release can be achieved through a rigorous selection of conservative mutations in accessory sites

A physical model for PDZ-domain/peptide interactions

The PDZ domain is an interaction motif that recognizes and binds the C-terminal peptides of target proteins. PDZ domains are ubiquitous in nature and help assemble multiprotein complexes that control cellular organization and signaling cascades. We present an optimized energy function to predict the binding free energy (ΔΔG) of PDZ domain/peptide interactions computationally. Geometry-optimized models of PDZ domain/peptide interfaces were built using Rosetta, and protein and peptide side chain and backbone degrees of freedom are minimized simultaneously. Using leave-one-out cross-validation, Rosetta’s energy function is adjusted to reproduce experimentally determined ΔΔG values with a correlation coefficient of 0.66 and a standard deviation of 0.79 kcal mol−1. The energy function places an increased weight on hydrogen bonding interactions when compared to a previously developed method to analyze protein/protein interactions. Binding free enthalpies (ΔΔH) and entropies (ΔS) are predicted with reduced accuracies of R = 0.60 and R = 0.17, respectively. The computational method improves prediction of PDZ domain specificity from sequence and allows design of novel PDZ domain/peptide interactions