Identification and structural analysis of a zebrafish apo and holo cellular retinol-binding protein

Cellular retinol-binding proteins (CRBPs) are cytoplasmic retinol-specific binding proteins. Mammalian CRBPs have been thoroughly characterised previously. Here we report on the identification and X-ray structural analysis of the apo (1.7A resolution) and holo (1.4A resolution) forms of a zebrafish CRBP. According to amino acid sequence and structure analyses, the zebrafish CRBP that we have identified resembles closely mammalian CRBP II, suggesting that it is the zebrafish orthologue of this mammalian CRBP type. Zebrafish CRBP forms a tight complex with all-trans retinol, producing an absorption spectrum similar to those of mammalian holo-CRBPs, albeit slightly blue-shifted. The superposition of the alpha-carbon atoms of the liganded (complexed with retinol) and unliganded forms of zebrafish CRBP shows significant differences in correspondence of the betaC-betaD (residues 55-58) and betaE-betaF (residues 74-77) turns, providing evidence for the occurrence of conformational changes accompanying retinol binding/release. Remarkable and well-defined ligand-dependent conformational changes in the protein region comprising the two beta-turns affect both the main chain and the side-chains of several residues. The two beta-turns project towards the interior of the cavity devoid of ligand of the apoprotein. The side-chains of F57, Y60 and L77 change substantially their orientation and position in the apoprotein relative to the holoprotein. In the beta-barrel internal cavity of apo-CRBP they occupy some of the space that is otherwise occupied by bound retinol in holo-CRBP, and are displaced from these positions on ligand binding. These results indicate that a flexible area encompassing the betaC-betaD and betaE-betaF turns may serve as the ligand portal and that these turns undergo conformational changes associated with the not yet clarified mechanism of retinol binding and release in CRBPs.     

Exploration of the conformational landscape in pregnane X receptor reveals a new binding pocket

Ligand‐regulated pregnane X receptor (PXR), a member of the nuclear receptor superfamily, plays a central role in xenobiotic metabolism. Despite its critical role in drug metabolism, PXR activation can lead to adverse drug‐drug interactions and early stage metabolism of drugs. Activated PXR can induce cancer drug resistance and enhance the onset of malignancy. Since promiscuity in ligand binding makes it difficult to develop competitive inhibitors targeting PXR ligand binding pocket (LBP), it is essential to identify allosteric sites for effective PXR antagonism. Here, molecular dynamics (MD) simulation studies unravelled the existence of two different conformational states, namely “expanded” and “contracted”, in apo PXR ligand binding domain (LBD). Ligand binding events shifted this conformational equilibrium and locked the LBD in a single “ligand‐adaptable” conformational state. Ensemble‐based computational solvent mapping identified a transiently open potential small molecule binding pocket between α5 and α8 helices, named “α8 pocket”, whose opening‐closing mechanism directly correlated with the conformational shift in LBD. A virtual hit identified through structure‐based virtual screening against α8 pocket locks the pocket in its open conformation. MD simulations further revealed that the presence of small molecule at allosteric site disrupts the LBD dynamics and locks the LBD in a “tightly‐contracted” conformation. The molecular details provided here could guide new structural studies to understand PXR activation and antagonism.     

Probing flexibility and "induced-fit" phenomena in aldose reductase by comparative crystal structure analysis and molecular dynamics simulations

Aldose reductase is a promising target for the treatment of diabetic complications, and as such, has become the focus of various drug design projects. As revealed by a survey of available crystal structures, the protein shows pronounced induced-fit effects upon ligand binding. Although helping to explain the enzyme's substrate promiscuity, phenomena of this kind are still responsible for significant complications in structure-based design efforts directed to aldose reductase. Accordingly, a deeper understanding of the principles governing conformational alterations in this enzyme would be of utmost practical importance. As a first step in addressing this issue, molecular dynamics (MD) simulations have been carried out. The ultrahigh resolution crystal structure of aldose reductase complexed with inhibitor IDD594 served as ideal starting point for a set of different simulations of nanosecond time scale: the native complexed state with bound inhibitor, the uncomplexed state (after removal of the inhibitor) at standard temperature, and the uncomplexed state at elevated temperature. The reference simulation of the complex exhibits extraordinary stability of the overall fold, whereas two distinct conformational substates are found for the binding-site region. In contrast, already at standard temperature pronounced changes are observed in the binding region during the simulation of the uncomplexed state. Leu300, for example, closes the access to the pocket opened by IDD594. On the other hand, conformations around the catalytic site are highly conserved, with the His110-Tyr48-NADP+ orientation being stabilized by a water molecule. Detailed analysis of the trajectories allows to reveal a set of distinct conformational substates that may prove useful as alternative structural templates in virtual screening for new aldose reductase inhibitors.     

Access Path to the Ligand Binding Pocket May Play a Role in Xenobiotics Selection by AhR

Understanding of multidrug binding at the atomic level would facilitate drug design and strategies to modulate drug metabolism, including drug transport, oxidation, and conjugation. Therefore we explored the mechanism of promiscuous binding of small molecules by studying the ligand binding domain, the PAS-B domain of the aryl hydrocarbon receptor (AhR). Because of the low sequence identities of PAS domains to be used for homology modeling, structural features of the widely employed HIF-2α and a more recent suitable template, CLOCK were compared. These structures were used to build AhR PAS-B homology models. We performed molecular dynamics simulations to characterize dynamic properties of the PAS-B domain and the generated conformational ensembles were employed in in silico docking. In order to understand structural and ligand binding features we compared the stability and dynamics of the promiscuous AhR PAS-B to other PAS domains exhibiting specific interactions or no ligand binding function. Our exhaustive in silico binding studies, in which we dock a wide spectrum of ligand molecules to the conformational ensembles, suggest that ligand specificity and selection may be determined not only by the PAS-B domain itself, but also by other parts of AhR and its protein interacting partners. We propose that ligand binding pocket and access channels leading to the pocket play equally important roles in discrimination of endogenous molecules and xenobiotics.     

Modeling of Arylamide Helix Mimetics in the p53 Peptide Binding Site of hDM2 Suggests Parallel and Anti-Parallel Conformations Are Both Stable

The design of novel α-helix mimetic inhibitors of protein-protein interactions is of interest to pharmaceuticals and chemical genetics researchers as these inhibitors provide a chemical scaffold presenting side chains in the same geometry as an α-helix. This conformational arrangement allows the design of high affinity inhibitors mimicking known peptide sequences binding specific protein substrates. We show that GAFF and AutoDock potentials do not properly capture the conformational preferences of α-helix mimetics based on arylamide oligomers and identify alternate parameters matching solution NMR data and suitable for molecular dynamics simulation of arylamide compounds. Results from both docking and molecular dynamics simulations are consistent with the arylamides binding in the p53 peptide binding pocket. Simulations of arylamides in the p53 binding pocket of hDM2 are consistent with binding, exhibiting similar structural dynamics in the pocket as simulations of known hDM2 binders Nutlin-2 and a benzodiazepinedione compound. Arylamide conformations converge towards the same region of the binding pocket on the 20 ns time scale, and most, though not all dihedrals in the binding pocket are well sampled on this timescale. We show that there are two putative classes of binding modes for arylamide compounds supported equally by the modeling evidence. In the first, the arylamide compound lies parallel to the observed p53 helix. In the second class, not previously identified or proposed, the arylamide compound lies anti-parallel to the p53 helix.     

Binding of retinoids and beta-carotene to beta-lactoglobulin. Influence of protein modifications

The binding of retinol, retinyl acetate, retinoic acid and beta-carotene to native, esterified and alkylated beta-lactoglobulin was followed by quenching of tryptophan fluorescence. Three studied retinoids bind to native or modified beta-lactoglobulin in 1:1 molar ratios, with apparent dissociation constants in the range of 10(-8) M. The maximum tryptophan fluorescence quenching of unmodified beta-lactoglobulin by beta-carotene is observed at the ligand/protein ratio of 1:2. Esterification and alkylation of beta-lactoglobulin shift the ratio of beta-carotene/protein to 1:1. In all the cases, except for retinoic acid binding to N-ethyllysyl-BLG, the performed chemical modifications of beta-lactoglobulin enhance protein binding affinity. Measured apparent dissociation constants of beta-carotene complexes with native and modified beta-lactoglobulin are an order of magnitude lower from binding constants of other studied retinoids.     

Role of proximal His93 in nitric oxide binding to metmyoglobin. Application of continuum solvation in Monte Carlo protein simulations.

Monte Carlo protein simulations with continuum solvation were used to explore the conformational mobility of NO within the active site of metmyoglobin. To the best of our knowledge this is the first application of a continuum solvation model for exploring protein binding sites. The usefulness of the Monte Carlo conformational analysis was demonstrated in comparative molecular dynamics simulations. Analysis of conformer populations revealed that Monte Carlo conformational analysis is more effective in mapping the relevant region of the potential surface. Distribution of low-energy conformers obtained for the metmyoglobin-NO complex was found to depend on the orientation of proximal His93. Different charge distributions corresponding to the two experimentally verified possible torsions of this proximal residue result in strong binding of NO or its release to a nearby hydrophobic trap. Conformer populations obtained by Monte Carlo conformational analysis were grouped into three main families: one, with the NO directly bound to the iron, appears when the CA-CB-CG-CD2 torsion of His93 is at its ligand binding value (-113 degrees); and two conformers exist where NO is trapped in a nearby hydrophobic pocket, the same cavity that was determined to be the geminate trap of CO in ferrous Mb as a result of the torsional flip of His93 to its ligand releasing state (-125 degrees). Based on this analysis, we suggest that the electrostatic rearrangement coupled to the conformational fluctuation of the proximal His leads to the geminate trapping of the ligand. Conformational rearrangement of the proximal side would provide the possibility of rebinding of the ligand to Fe.     

Expect the unexpected or caveat for drug designers: multiple structure determinations using aldose reductase crystals treated under varying soaking and co-crystallisation conditions

In structure-based drug design, accurate crystal structure determination of protein-ligand complexes is of utmost importance in order to elucidate the binding characteristics of a putative lead to a given target. It is the starting point for further design hypotheses to predict novel leads with improved properties. Often, crystal structure determination is regarded as ultimate proof for ligand binding providing detailed insight into the specific binding mode of the ligand to the protein. This widely accepted practise relies on the assumption that the crystal structure of a given protein-ligand complex is unique and independent of the protocol applied to produce the crystals. We present two examples indicating that this assumption is not generally given, even though the composition of the mother liquid for crystallisation was kept unchanged: Multiple crystal structure determinations of aldose reductase complexes obtained under varying crystallisation protocols concerning soaking and crystallisation exposure times were performed resulting in a total of 17 complete data sets and ten refined crystal structures, eight in complex with zopolrestat and two complexed with tolrestat. In the first example, a flip of a peptide bond is observed, obviously depending on the crystallisation protocol with respect to soaking and co-crystallisation conditions. This peptide flip is accompanied by a rupture of an H-bond formed to the bound ligand zopolrestat. The indicated enhanced local mobility of the complex is in agreement with the results of molecular dynamics simulations. As a second example, the aldose reductase-tolrestat complex is studied. Unexpectedly, two structures could be obtained: one with one, and a second with four inhibitor molecules bound to the protein. They are located in and near the binding pocket facilitated by crystal packing effects. Accommodation of the four ligand molecules is accompanied by pronounced shifts concerning two helices interacting with the additional ligands.     

Complete sequence and model for the A2 subunit of the carotenoid pigment complex, crustacyanin

The complete sequence has been determined for the A2 subunit of crustacyanin, an astaxanthin-binding protein from the carapace of the lobster Homarus gammarus. The polypeptide chain is 174 residues long and is similar to proteins of the retinol-binding protein superfamily. Some regions of the sequence are most similar to the retinol-binding protein, beta-lactoglobulin subgroup, while the disulphide bonding pattern is more akin to that seen in the porphyrin binding proteins insecticyanin and bilin-binding protein. It is beginning to appear as though this superfamily of proteins, characterized by a similar gross structural framework, may be further subdivided into interrelated subclasses. Model building based on the coordinates of the known structure of human plasma retinol-binding protein and on empirical prediction algorithms has allowed the putative identification of side-chains which line the binding cavity. This pocket is larger than in retinol binding protein and beta-lactoglobulin but does not allow the carotenoid to adopt a folded conformation. The amino acid composition of the pocket does not support a 'charge-shift'-type hypothesis to support the bathochromic shift phenomenon which takes place on interaction of the chromophore with the protein. Instead aromatic side-chains may play a prominent role.     

On the binding determinants of the glutamate agonist with the glutamate receptor ligand binding domain

Ionotropic glutamate receptors (GluRs) are ligand-gated membrane channel proteins found in the central neural system that mediate a fast excitatory response of neurons. In this paper, we report theoretical analysis of the ligand-protein interactions in the binding pocket of the S1S2 (ligand binding) domain of the GluR2 receptor in the closed conformation. By utilizing several theoretical methods ranging from continuum electrostatics to all-atom molecular dynamics simulations and quantum chemical calculations, we were able to characterize in detail glutamate agonist binding to the wild-type and E705D mutant proteins. A theoretical model of the protein-ligand interactions is validated via direct comparison of theoretical and Fourier transform infrared spectroscopy (FTIR) measured frequency shifts of the ligand's carboxylate group vibrations [Jayaraman et al. (2000) Biochemistry 39, 8693-8697; Cheng et al. (2002) Biochemistry 41, 1602-1608]. A detailed picture of the interactions in the binding site is inferred by analyzing contributions to vibrational frequencies produced by protein residues forming the ligand-binding pocket. The role of mobility and hydrogen-bonding network of water in the ligand-binding pocket and the contribution of protein residues exposed in the binding pocket to the binding and selectivity of the ligand are discussed. It is demonstrated that the molecular surface of the protein in the ligand-free state has mainly positive electrostatic potential attractive to the negatively charged ligand, and the potential produced by the protein in the ligand-binding pocket in the closed state is complementary to the distribution of the electrostatic potential produced by the ligand itself. Such charge complementarity ensures specificity to the unique charge distribution of the ligand.     

Dynamics of cholesterol exchange in the oxysterol binding protein family

The oxysterol-binding protein-related protein (ORP) family is essential to sterol transfer and sterol-dependent signal transduction in eukaryotes. The crystal structure of one ORP family member, yeast Osh4, is known in apo and sterol-bound states. In the bound state, a 29 residue N-terminal lid region covers the opening of the cholesterol-binding tunnel, preventing cholesterol exchange. Equilibrium and steered molecular dynamics (MD) simulations of Osh4 were carried out to characterize the mechanism of cholesterol exchange. While most of the structural core was stable during the simulations, the lid was partly opened in the apo equilibrium MD simulation. Helix α7, which undergoes the largest conformational change in the crystallized bound and apo states, is conformationally coupled to the opening of the lid. The movement of α7 helps create a docking site for donor or acceptor membranes in the open state. In the steered MD simulations of cholesterol dissociation, we observed complete opening of the lid covering the cholesterol-binding tunnel. Cholesterol was found to exit the binding pocket in a step-wise process involving (i) the breaking of water-mediated hydrogen bonds and van der Waals contacts within the binding pocket, (ii) opening of the lid covering the binding pocket, and (iii) breakage of transient cholesterol contacts with the rim of the pocket and hydrophobic residues on the interior face of the lid.     

EF loop conformational change triggers ligand binding in beta-lactoglobulins

Beta-lactoglobulins, belonging to the lipocalin family, are a widely studied group of proteins, characterized by the ability to solubilize and transport hydrophobic ligands, especially fatty acids. Despite many reports, the mechanism of ligand binding and the functional role of these proteins is still unclear, and many contradicting concepts are often encountered in the literature. In the present paper the comparative analysis of the binding properties of beta-lactoglobulins has been performed using sequence-derived information, structure-based electrostatic calculations, docking simulations, and NMR experiments. Our results reveal for the first time the mechanism of beta-lactoglobulin ligand binding, which is completely determined by the opening-closing of EF loop, triggered by Glu89 protonation. The alkaline shift observed for Glu89 pKa in porcine beta-lactoglobulin (pKa 9.7) with respect to the bovine species (pKa 5.5) depends upon the interplay of electrostatic effects of few nearby key residues. Porcine protein is therefore able to bind fatty acids provided that the appropriate pH solution conditions are met (pH > 8.6), where the EF loop conformational change can take place. The unusually high pH of binding detected for porcine beta-lactoglobulin seems to be functional to lipases activity. Theoretical pKa calculations extended to representative beta-lactoglobulins allowed the identification of key residues involved in structurally and functionally important electrostatic interactions. The results presented here provide a strong indication that the described conformational change is a common feature of all beta-lactoglobulins.     

ATP dependent NS3 helicase interaction with RNA: insights from molecular simulations

Non-structural protein 3 (NS3) helicase from hepatitis C virus is an enzyme that unwinds and translocates along nucleic acids with an ATP-dependent mechanism and has a key role in the replication of the viral RNA. An inchworm-like mechanism for translocation has been proposed based on crystal structures and single molecule experiments. We here perform atomistic molecular dynamics in explicit solvent on the microsecond time scale of the available experimental structures. We also construct and simulate putative intermediates for the translocation process, and we perform non-equilibrium targeted simulations to estimate their relative stability. For each of the simulated structures we carefully characterize the available conformational space, the ligand binding pocket, and the RNA binding cleft. The analysis of the hydrogen bond network and of the non-equilibrium trajectories indicates an ATP-dependent stabilization of one of the protein conformers. Additionally, enthalpy calculations suggest that entropic effects might be crucial for the stabilization of the experimentally observed structures.     

Complexes of porcine odorant binding protein with odorant molecules belonging to different chemical classes

Porcine odorant binding protein (pOBP) is a monomer of 157 amino acid residues, purified in abundance from pig nasal mucosa. In contrast to the observation on lipocalins as retinol binding protein (RBP), major urinary protein (MUP) or bovine odorant binding protein (bOBP), no naturally occurring ligand was found in the beta-barrel cavity of pOBP. Porcine OBP was therefore chosen as a simple model for structure/function studies with odorant molecules. In competition experiments with tritiated pyrazine, the affinity of pOBP towards several odorant molecules belonging to different chemical classes has been found to be of the micromolar order, with a 1:1 stoichiometry. The X-ray structures of pOBP complexed to these molecules were determined at resolution between 2.15 and 1.4 A. As expected, the electron density of the odorant molecules was observed into the hydrophobic beta-barrel of the lipocalin. Inside this cavity, very few specific interactions were established between the odorant molecule and the amino acid side-chains, which did not undergo significant conformational change. The high B-factors observed for the odorant molecules as well as the existence of alternative conformations reveal a non-specific mode of binding of the odorant molecules in the cavity.     

Interplay between structural rigidity and electrostatic interactions in the ligand binding domain of GluR2

Using molecular dynamics (MD) simulations, computational protein modifications, and a novel theoretical methodology that determines structural rigidity/flexibility (the FIRST algorithm), we investigate how molecular structure and dynamics of the glutamate receptor ligand binding domain (GluR2 S1S2) facilitate its conformational transition. S1S2 is a two-lobe protein, which undergoes a cleft closure conformational transition upon binding an agonist in the cleft between the two lobes; hence it is expected that the mechanism of this conformational transition can be characterized as a hinge-type. However, in the rigidity analysis one lobe of the protein is identified as a single rigid cluster while the other one is structurally flexible, inconsistent with a presumed mechanical hinge mechanism. Instead, we characterize the cleft-closing transition as a load and lock mechanism. We find that when two cross-cleft hydrogen bonds are disrupted the protein undergoes a rapid cleft opening transition. At the same time, the dynamical behavior of the cleft in the presence of the glutamate ligand is only weakly affected by the S652 peptide bond in its flipped conformation observed in the crystal structure. The residue E705 plays significant role in stabilization of the closed conformation via electrostatic interactions. The presence of the E705-K730 salt bridge seems to correlate strongly withthe cleft opening transition in the MD simulations.     

On CO Egress from,and Re‐Uptake by,the Enzyme MauG,as a Mimic of the Acquisition of Oxidizing Agents by the pre‐MADH_MauG System. A Molecular Mechanics Approach

In this work, by applying a non‐deterministic, randomly‐oriented minimal force to the dissociated CO ligand of the MauG‐CO system, the molecular‐dynamics (MD) behavior of this system could be quickly unraveled. It turned out that CO has no marked directional egress from the high‐spin c‐heme iron distal pocket. Rather, CO is able to exploit all interstices created during the protein fluctuations. Nonetheless, no steady route toward the surrounding solvent was ever observed: CO jumped first into other binding pockets before being able to escape the protein. In a few cases, on hitting the surrounding H₂O molecules, CO was observed to reverse direction, re‐entering the protein. A contention that conformational inversion of the P107 ring provides a gate to the iron ion is not supported by the present simulations.     

Isolation of a new ligand-carrying casein fragment from bovine mammary gland microsomes.

Whilst looking for components involved in retinol metabolism in secreting mammary gland cells, a 12 kDa protein was isolated. This protein had bound a ligand with characteristics of retinol. N-Terminal sequencing and amino acid analysis showed that this protein is highly homologous with an alpha-s1-casein fragment. No ligand was found for beta-lactoglobulin, previously thought to be involved in retinol metabolism.     

Evidence of chemical exchange in recombinant Major Urinary Protein and quenching thereof upon pheromone binding

The internal dynamics of recombinant Major Urinary Protein (rMUP) have been investigated by monitoring transverse nitrogen-15 relaxation using multiple-echo Carr–Purcell–Meiboom–Gill (CPMG) experiments. While the ligand-free protein (APO-rMUP) features extensive evidence of motions on the milliseconds time scale, the complex with 2-methoxy-3-isobutylpyrazine (HOLO-rMUP) appears to be much less mobile on this time scale. At 308 K, exchange rates k _ex = 500–2000 s⁻¹ were typically observed in APO-rMUP for residues located adjacent to a β-turn comprising residues 83–87. These residues occlude an entry to the binding pocket and have been proposed to be a portal for ligand entry in other members of the lipocalin family, such as the retinol binding protein and the human fatty-acid binding protein. Exchange rates and populations are largely uncorrelated, suggesting local ‘breathing’ motions rather than a concerted global conformational change. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Guanidine-II aptamer conformations and ligand binding modes through the lens of molecular simulation

Regulation of gene expression via riboswitches is a widespread mechanism in bacteria. Here, we investigate ligand binding of a member of the guanidine sensing riboswitch family, the guanidine-II riboswitch (Gd-II). It consists of two stem–loops forming a dimer upon ligand binding. Using extensive molecular dynamics simulations we have identified conformational states corresponding to ligand-bound and unbound states in a monomeric stem–loop of Gd-II and studied the selectivity of this binding. To characterize these states and ligand-dependent conformational changes we applied a combination of dimensionality reduction, clustering, and feature selection methods. In absence of a ligand, the shape of the binding pocket alternates between the conformation observed in presence of guanidinium and a collapsed conformation, which is associated with a deformation of the dimerization interface. Furthermore, the structural features responsible for the ability to discriminate against closely related analogs of guanidine are resolved. Based on these insights, we propose a mechanism that couples ligand binding to aptamer dimerization in the Gd-II system, demonstrating the value of computational methods in the field of nucleic acids research.