Molecular Dynamics Simulation of Ligand Dissociation from Liver Fatty Acid Binding Protein

The mechanisms of how ligands enter and leave the binding cavity of fatty acid binding proteins (FABPs) have been a puzzling question over decades. Liver fatty acid binding protein (LFABP) is a unique family member which accommodates two molecules of fatty acids in its cavity and exhibits the capability of interacting with a variety of ligands with different chemical structures and properties. Investigating the ligand dissociation processes of LFABP is thus a quite interesting topic, which however is rather difficult for both experimental approaches and ordinary simulation strategies. In the current study, random expulsion molecular dynamics simulation, which accelerates ligand motions for rapid dissociation, was used to explore the potential egress routes of ligands from LFABP. The results showed that the previously hypothesized “portal region” could be readily used for the dissociation of ligands at both the low affinity site and the high affinity site. Besides, one alternative portal was shown to be highly favorable for ligand egress from the high affinity site and be related to the unique structural feature of LFABP. This result lends strong support to the hypothesis from the previous NMR exchange studies, which in turn indicates an important role for this alternative portal. Another less favored potential portal located near the N-terminal end was also identified. Identification of the dissociation pathways will allow further mechanistic understanding of fatty acid uptake and release by computational and/or experimental techniques.     

A single amino acid mutation in zebrafish (Danio rerio) liver bile acid-binding protein can change the stoichiometry of ligand binding

In all of the liver bile acid-binding proteins (L-BABPs) studied so far, it has been found that the stoichiometry of binding is of two cholate molecules per internal binding site. In this paper, we describe the expression, purification, crystallization, and three-dimensional structure determination of zebrafish (Danio rerio) L-BABP to 1.5A resolution, which is currently the highest available for a protein of this family. Since we have found that in zebrafish, the stoichiometry of binding in the protein cavity is of only one cholate molecule per wild type L-BABP, we examined the role of two crucial amino acids present in the binding site. Using site-directed mutagenesis, we have prepared, crystallized, and determined the three-dimensional structure of co-crystals of two mutants. The mutant G55R has the same stoichiometry of binding as the wild type protein, whereas the C91T mutant changes the stoichiometry of binding from one to two ligand molecules in the cavity and therefore appears to be more similar to the other members of the L-BABP family. Based on the presence or absence of a single disulfide bridge, it can be postulated that fish should bind a single cholate molecule, whereas amphibians and higher vertebrates should bind two. Isothermal titration calorimetry has also revealed the presence in the wild type protein and the G55R mutant of an additional binding site, different from the first and probably located on the surface of the molecule.     

Solution structure of the highly acidic protein HI1450 from Haemophilus influenzae, a putative double-stranded DNA mimic

The solution structure of the acidic protein HI1450 from Haemophilus influenzae has been determined by NMR spectroscopy. HI1450 has homologues in ten other bacterial species including Escherichia coli, Vibrio cholerae, and Yersinia pestis but there are no functional assignments for any members of the family. Thirty-one of the amino acids in this 107-residue protein are aspartates or glutamates, contributing to an unusually low pI of 3.72. The secondary structure elements are arranged in an alpha-alpha-beta-beta-beta-beta order with the two alpha helices packed against the same side of an anti-parallel four-stranded beta meander. Two large loops, one between beta1 and beta2 and the other between beta2 and beta3 bend almost perpendicularly across the beta-strands in opposite directions on the non-helical side of the beta-sheet to form a conserved hydrophobic cavity. The HI1450 structure has some similarities to the structure of the double-stranded DNA (dsDNA) mimic uracil DNA glycosylase inhibitor (Ugi) including the distribution of surface charges and the position of the hydrophobic cavity. Based on these similarities, as well as having a comparable molecular surface to dsDNA, we propose that HI1450 may function as a dsDNA mimic in order to inhibit or regulate an as yet unidentified dsDNA binding protein.     

Insight into the interaction sites between fatty acid binding proteins and their ligands

Fatty acid binding proteins (FABPs), are evolutionarily conserved small cytoplasmic proteins that occur in many tissue-specific types. One of their primary functions is to facilitate the clearance of the cytoplasmic matrix from free fatty acids and of other detergent-like compounds. Crystallographic studies of FABP proteins have revealed a well defined binding site located deep inside their β-clam structure that is hardly exposed to the bulk solution. However, NMR measurements revealed that, when the protein is equilibrated with its ligands, residues that are clearly located on the outer surface of the protein do interact with the ligand. To clarify this apparent contradiction we applied molecular dynamics simulations to follow the initial steps associated with the FABP–fatty acid interaction using, as a model, the interaction of toad liver basic FABP, or chicken liver bile acid binding protein, with a physiological concentration of palmitate ions. The simulations (~200 ns of accumulated time) show that fatty acid molecules interact, unevenly, with various loci on the protein surface, with the favored regions being the portal and the anti-portal domains. Random encounters with palmitate at these regions led to lasting adsorption to the surface, while encounters at the outer surface of the β-clam were transient. Therefore, we suggest that the protein surface is capable of sequestering free fatty acids from solution, where brief encounters evolve into adsorbed states, which later mature by migration of the ligand into a more specific binding site.     

Discovery of Fragment Molecules That Bind the Human Peroxiredoxin 5 Active Site

The search for protein ligands is a crucial step in the inhibitor design process. Fragment screening represents an interesting method to rapidly find lead molecules, as it enables the exploration of a larger portion of the chemical space with a smaller number of compounds as compared to screening based on drug-sized molecules. Moreover, fragment screening usually leads to hit molecules that form few but optimal interactions with the target, thus displaying high ligand efficiencies. Here we report the screening of a homemade library composed of 200 highly diverse fragments against the human Peroxiredoxin 5 protein. Peroxiredoxins compose a family of peroxidases that share the ability to reduce peroxides through a conserved cysteine. The three-dimensional structures of these enzymes ubiquitously found throughout evolution have been extensively studied, however, their biological functions are still not well understood and to date few inhibitors have been discovered against these enzymes. Six fragments from the library were shown to bind to the Peroxiredoxin 5 active site and ligand-induced chemical shift changes were used to drive the docking of these small molecules into the protein structure. The orientation of the fragments in the binding pocket was confirmed by the study of fragment homologues, highlighting the role of hydroxyl functions that hang the ligands to the Peroxiredoxin 5 protein. Among the hit fragments, the small catechol molecule was shown to significantly inhibit Peroxiredoxin 5 activity in a thioredoxin peroxidase assay. This study reports novel data about the ligand-Peroxiredoxin interactions that will help considerably the development of potential Peroxiredoxin inhibitors.     

Solution structure of ileal lipid binding protein in complex with glycocholate.

Ileal lipid binding protein (ILBP) is a cytosolic lipid-binding protein that binds both bile acids and fatty acids. We have determined the solution structure of porcine ILBP in complex with glycocholate by homonuclear and heteronuclear two-dimensional NMR spectroscopy. The conformation of the protein-ligand complex was determined by restrained energy minimization and simulated annealing calculations after docking the glycocholate ligand into the protein structure. The overall tertiary structure of ILBP is highly analogous to the three-dimensional structures of several other intracellular lipid binding proteins (LBPs). Like the apo-structure, the bile-acid complex of ILBP is composed of 10 anti-parallel beta-strands that form a water-filled clam-shell structure, and two short alpha-helices. Chemical shift data indicated that the bile acid ligand is bound inside the protein cavity. Furthermore, 13C-edited heteronuclear single-quantum correlation-NOESY experiments showed NOE contacts between several aromatic residues located in the proposed bile acid portal region and the 13C-labeled ligand. A single bile acid molecule is bound inside the protein, with the steroid moiety penetrating deep into the water-accessible internal cavity, such that ring A is located right above the plane of the Trp49 indole ring. The carboxylate tail of the ligand is protruding from the proposed bile acid portal into the surrounding aqueous solution. The body of the steroid moiety is oriented with the nonpolar face in contact with the mostly hydrophobic residues of beta-strands C, D and E, while the polar face shows contacts with the side-chains of Tyr97, His99, Glu110 and Arg121 in beta-strands H, I and J. Thus, the conformational arrangement of the ligand complex suggests that the binding affinity of ILBP for bile acid molecules is based mainly on strong hydrophobic interactions inside the protein cavity. Furthermore, this binding mode explains how ILBP can transport unconjugated and conjugated bile acids.     

Determining molecular binding sites on human serum albumin by displacement of oleic acid

An NMR method was developed for determining binding sites of small molecules on human serum albumin (HSA) by competitive displacement of (13)C-labeled oleic acid. This method is based on the observation that in the crystal structure of HSA complexed with oleic acid, two principal drug-binding sites, Sudlow's sites I (warfarin) and II (ibuprofen), are also occupied by fatty acids. In two-dimensional [(1)H,(13)C]heteronuclear single quantum coherence NMR spectra, seven distinct resonances were observed for the (13)C-methyl-labeled oleic acid as a result of its binding to HSA. Resonances corresponding to the major drug-binding sites were identified through competitive displacement of molecules that bind specifically to each site. Thus, binding of molecules to these sites can be followed by their displacement of oleic acids. Furthermore, the amount of bound ligand at each site can be determined from changes in resonance intensities. For molecules containing fluorine, binding results were further validated by direct observations of the bound ligands using (19)F NMR. Identifying the binding sites for drug molecules on HSA can aid in determining the structure-activity relationship of albumin binding and assist in the design of molecules with altered albumin binding.     

Studies on ligand binding to histidine triad nucleotide binding protein 1

Histidine triad nucleotide binding protein (HINT1) is an intracellular protein that binds purine mononucleotides. Strong sequence conservation suggests that these proteins play a fundamental role in cell biology, however its exact cellular function continues to remain elusive. nuclear magnetic resonance (NMR) studies using STD and HSQC were conducted to observe ligand binding to HINT1. These studies were confirmed using fluorescence spectroscopy titrations. We found that AICAR, the first non-phosphate containing ligand, binds to mouse histidine triad nucleotide binding protein 1 (HINT1). Chemical shift perturbations are mapped onto the X-ray structure showing AICAR binds at the same site as GMP. The NMR results demonstrated that this method will be valuable for the future screening of small molecules that can be used to modulate the function of HINT1.     

Crystal structure of YecO from Haemophilus influenzae (HI0319) reveals a methyltransferase fold and a bound S-adenosylhomocysteine.

The crystal structure of YecO from Haemophilus influenzae (HI0319), a protein annotated in the sequence databases as hypothetical, and that has not been assigned a function, has been determined at 2.2-A resolution. The structure reveals a fold typical of S-adenosyl-L-methionine-dependent (AdoMet) methyltransferase enzymes. Moreover, a processed cofactor, S-adenosyl-L-homocysteine (AdoHcy), is bound to the enzyme, further confirming the biochemical function of HI0319 and its sequence family members. An active site arginine, shielded from bulk solvent, interacts with an anion, possibly a chloride ion, which in turn interacts with the sulfur atom of AdoHcy. The AdoHcy and nearby protein residues delineate a small solvent-excluded substrate binding cavity of 162 A(3) in volume. The environment surrounding the cavity indicates that the substrate molecule contains a hydrophobic moiety and an anionic group. Many of the residues that define the cavity are invariant in the HI0319 sequence family but are not conserved in other methyltransferases. Therefore, the substrate specificity of YecO enzymes is unique and differs from the substrate specificity of all other methyltransferases sequenced to date. Examination of the Enzyme Commission list of methyltransferases prompted a manual inspection of 10 possible substrates using computer graphics and suggested that the ortho-substituted benzoic acids fit best in the active site.     

Ligand Binding Modulates the Structural Dynamics and Compactness of the Major Birch Pollen Allergen

Pathogenesis-related plant proteins of class-10 (PR-10) are essential for storage and transport of small molecules. A prominent member of the PR-10 family, the major birch pollen allergen Bet v 1, is the main cause of spring pollinosis in the temperate climate zone of the northern hemisphere. Bet v 1 binds various ligand molecules to its internal cavity, and immunologic effects of the presence of ligand have been discussed. However, the mechanism of binding has remained elusive. In this study, we show that in solution Bet v 1.0101 is conformationally heterogeneous and cannot be represented by a single structure. NMR relaxation data suggest that structural dynamics are fundamental for ligand access to the protein interior. Complex formation then leads to significant rigidification of the protein along with a compaction of its 3D structure. The data presented herein provide a structural basis for understanding the immunogenic and allergenic potential of ligand binding to Bet v 1 allergens.     

Double threading through DNA: NMR structural study of a bis-naphthalene macrocycle bound to a thymine-thymine mismatch

The macrocyclic bis-naphthalene macrocycle (2,7-BisNP), belonging to the cyclobisintercalator family of DNA ligands, recognizes T-T mismatch sites in duplex DNA with high affinity and selectivity, as evidenced by thermal denaturation experiments and NMR titrations. The binding of this macrocycle to an 11-mer DNA oligonucleotide containing a T-T mismatch was studied using NMR spectroscopy and NMR-restrained molecular modeling. The ligand forms a single type of complex with the DNA, in which one of the naphthalene rings of the ligand occupies the place of one of the mismatched thymines, which is flipped out of the duplex. The second naphthalene unit of the ligand intercalates at the A-T base pair flanking the mismatch site, leading to encapsulation of its thymine residue via double stacking. The polyammonium linking chains of the macrocycle are located in the minor and the major grooves of the oligonucleotide and participate in the stabilization of the complex by formation of hydrogen bonds with the encapsulated thymine base and the mismatched thymine remaining inside the helix. The study highlights the uniqueness of this cyclobisintercalation binding mode and its importance for recognition of DNA lesion sites by small molecules.     

The crystal structure of Rv0813c from Mycobacterium tuberculosis reveals a new family of fatty acid-binding protein-like proteins in bacteria

The gene Rv0813c from Mycobacterium tuberculosis, which codes for a hypothetical protein of unknown function, is conserved within the order Actinomycetales but absent elsewhere. The crystal structure of Rv0813c reveals a new family of proteins that resemble the fatty acid-binding proteins (FABPs) found in eukaryotes. Rv0813c adopts the 10-stranded beta-barrel fold typical of FABPs but lacks the double-helix insert that covers the entry to the binding site in the eukaryotic proteins. The barrel encloses a deep cavity, at the bottom of which a small cyclic ligand was found to bind to the hydroxyl group of Tyr192. This residue is part of a conserved Arg-X-Tyr motif much like the triad that binds the carboxylate group of fatty acids in FABPs. Most of the residues forming the internal surface of the cavity are conserved in homologous protein sequences found in CG-rich prokaryotes, strongly suggesting that Rv0813c is a member of a new family of bacterial FABP-like proteins that may have roles in the recognition, transport, and/or storage of small molecules in the bacterial cytosol.     

Molecular dynamics simulations of the adipocyte lipid binding protein reveal a novel entry site for the ligand

The adipocyte lipid binding protein (ALBP) binds fatty acids (FA) in a cavity that is inaccessible from the bulk. Therefore, the penetration of the FA necessitates conformational changes whose nature is still unknown. It was suggested that the lipid first enters through a "portal region" which consists of the alphaII helix and the adjacent tight turns. The initial event in the ligand binding must be the interaction of the lipid with the protein surface. To analyze this interaction, we have carried out three molecular dynamics simulations of the apo-ALBP, with a palmitate ion initially located at different positions near the protein surface. The simulation indicated that the ligand could adsorb to the protein in more than one location. Yet, in one case, the ligand managed to penetrate the protein by entering a newly formed cavity some 10 A deep. The entry site is located near the N-terminus, at the junction between the loops connecting the beta-strands. Further progression of the penetration seems to be arrested by the need for desolvation of the COOH end of the headgroup. Evolutionary analysis showed that amino acids in this entry site are well conserved. On the basis of these observations, we suggest that the ligand may enter the protein from a site other than the portal region. Furthermore, the rate-limiting step is proposed to be the desolvation of the FA polar headgroup.     

High Selectivity of the γ-Aminobutyric Acid Transporter 2 (GAT-2, SLC6A13) Revealed by Structure-based Approach

The solute carrier 6 (SLC6) is a family of ion-dependent transporters that mediate uptake into the cell of osmolytes such as neurotransmitters and amino acids. Four SLC6 members transport GABA, a key neurotransmitter that triggers inhibitory signaling pathways via various receptors (e.g., GABA_A). The GABA transporters (GATs) regulate the concentration of GABA available for signaling and are thus targeted by a variety of anticonvulsant and relaxant drugs. Here, we characterize GAT-2, a transporter that plays a role in peripheral GABAergic mechanisms, by constructing comparative structural models based on crystallographic structures of the leucine transporter LeuT. Models of GAT-2 in two different conformations were constructed and experimentally validated, using site-directed mutagenesis. Computational screening of 594,166 compounds including drugs, metabolites, and fragment-like molecules from the ZINC database revealed distinct ligands for the two GAT-2 models. 31 small molecules, including high scoring compounds and molecules chemically related to known and predicted GAT-2 ligands, were experimentally tested in inhibition assays. Twelve ligands were found, six of which were chemically novel (e.g., homotaurine). Our results suggest that GAT-2 is a high selectivity/low affinity transporter that is resistant to inhibition by typical GABAergic inhibitors. Finally, we compared the binding site of GAT-2 with those of other SLC6 members, including the norepinephrine transporter and other GATs, to identify ligand specificity determinants for this family. Our combined approach may be useful for characterizing interactions between small molecules and other membrane proteins, as well as for describing substrate specificities in other protein families.     

A novel haem-binding interface in the 22 kDa haem-binding protein p22HBP

The 22 kDa haem-binding protein, p22HBP, is highly expressed in erythropoietic tissues and binds to a range of metallo- and non-metalloporphyrin molecules with similar affinities, suggesting a role in haem regulation or synthesis. We have determined the three-dimensional solution structure of p22HBP and mapped the porphyrin-binding site, which comprises a number of loops and a alpha-helix all located on a single face of the molecule. The structure of p22HBP is related to the bacterial multi-drug resistance protein BmrR, and is the first protein with this fold to be identified in eukaryotes. Strikingly, the porphyrin-binding site in p22HBP is located in a similar position to the drug-binding site of BmrR. These similarities suggest that the broad ligand specificity observed for both BmrR and p22HBP may result from a conserved ligand interaction mechanism. Taken together, these data suggest that the both the fold and its associated function, that of binding to a broad range of small hydrophobic molecules, are ancient, and have been adapted throughout evolution for a variety of purposes.     

Toward novel inhibitors against KdsB: a highly specific and selective broad-spectrum bacterial enzyme

KdsB (3-deoxy-manno-octulosonate cytidylyltransferase) is a highly specific and selective bacterial enzyme that catalyzes KDO (3-Deoxy-D-mano-oct-2-ulosonic acid) activation in KDO biosynthesis pathway. Failure in KDO biosynthesis causes accumulation of lipid A in the bacterial outer membrane that leads to cell growth arrest. This study reports a combinatorial approach comprising virtual screening of natural drugs library, molecular docking, computational pharmacokinetics, molecular dynamics simulation, and binding free energy calculations for the identification of potent lead compounds against the said enzyme. Virtual screening demonstrated 1460 druglike compounds in a total of 4800, while molecular docking illustrated Ser13, Arg14, and Asp236 as the anchor amino acids for recognizing and binding the inhibitors. Functional details of the enzyme in complex with the best characterized compound-226 were explored through two hundred nanoseconds of MD simulation. The ligand after initial adjustments jumps into the active cavity, followed by the deep cavity, and ultimately backward rotating movement toward the initial docked site of the pocket. During the entire simulation period, Asp236 remained in contact with the ligand and can be considered as a major catalytic residue of the enzyme. Radial distribution function confirmed that toward the end of the simulation, strengthening of ligand-receptor occurred with ligand and enzyme active residues in close proximity. Binding free energy calculations via MM(PB/GB)SA and Waterswap reaction coordinates, demonstrated the high affinity of the compound for enzyme active site residues. These findings can provide new avenues for designing potent compounds against notorious bacterial pathogens.     

Flexibility and ligand exchange in a buried cavity mutant of T4 lysozyme studied by multinuclear NMR

The Leu99-->Ala mutant of T4 lysozyme contains a large internal cavity in the core of its C-terminal domain that is capable of reversibly binding small hydrophobic compounds. Although the cavity is completely buried, molecules such as benzene or xenon can exchange rapidly in and out. The dynamics of the unliganded protein have been compared to the wild-type protein by measuring the NMR spin relaxation rates of backbone amide and side chain methyl nuclei. Many residues surrounding the cavity were found to be affected by a chemical exchange process with a rate of 1500 +/- 200 s(-1), which is quenched upon addition of saturating amounts of the ligand xenon. The relationship between the structure, dynamics, and energetics of the T4 lysozyme mutant is discussed.     

Water movement during ligand unbinding from receptor site

An 1-ns unbinding trajectory of retinol from the bovine serum retinol-binding protein has been obtained from molecular dynamics simulations. The behavior of water during ligand unbinding has never been studied in detail. I described a new method for defining a binding site, located the water molecules involved in the binding site, and examined their movements during unbinding. I found that there were only small changes in the binding site. During unbinding, the number of water molecules inside the binding site decreased, with some water molecules exhibiting movements similar in magnitude to bulk water, and there were rearrangements of the hydrogen bonds. This work represents the first detailed study of the behavior of water during an unbinding process.