Lipocalin Q83 reveals a dual ligand binding mode with potential implications for the functions of siderocalins

Siderocalins are particular lipocalins that participate in the innate immune response by interfering with bacterial siderophore-mediated iron uptake. Additionally, siderocalins are involved in several physiological and pathological processes such as inflammation, iron delivery, tissue differentiation, and cancer progression. Here we show that siderocalin Q83 displays an unexpected dual ligand binding mode as it can bind enterobactin and unsaturated fatty acids simultaneously. The solution structure of the siderocalin Q83 in complex with arachidonic acid and enterobactin reveals molecular details of this novel dual binding mode and the determinants of fatty acid binding specificity. Our results suggest that Q83 is a metabolic hub linking iron and fatty acid pathways. This unexpected coupling might contribute to the pleiotropic functions of siderocalins.     

The v-myc-induced Q83 Lipocalin Is a Siderocalin

Siderocalins are atypical lipocalins able to capture siderophores with high affinity. They contribute to the innate immune response by interfering with bacterial siderophore-mediated iron uptake but are also involved in numerous physiological processes such as inflammation, iron delivery, tissue differentiation, and cancer progression. The Q83 lipocalin was originally identified based on its overexpression in quail embryo fibroblasts transformed by the v-myc oncogene. We show here that Q83 is a siderocalin, binding the siderophore enterobactin with an affinity and mode of binding nearly identical to that of neutrophil gelatinase-associated lipocalin (NGAL), the prototypical siderocalin. This strengthens the role of siderocalins in cancer progression and inflammation. In addition, we also present the solution structure of Q83 in complex with intact enterobactin and a detailed analysis of the Q83 binding mode, including mutagenesis of the critical residues involved in enterobactin binding. These data provide a first insight into the molecular details of siderophore binding and delineate the common molecular properties defining the siderocalin protein family.     

Protein dynamics and motions in relation to their functions: several case studies and the underlying mechanisms

Proteins are dynamic entities in cellular solution with functions governed essentially by their dynamic personalities. We review several dynamics studies on serine protease proteinase K and HIV-1 gp120 envelope glycoprotein to demonstrate the importance of investigating the dynamic behaviors and molecular motions for a complete understanding of their structure–function relationships. Using computer simulations and essential dynamic (ED) analysis approaches, the dynamics data obtained revealed that: (i) proteinase K has highly flexible substrate-binding site, thus supporting the induced-fit or conformational selection mechanism of substrate binding; (ii) Ca²⁺ removal from proteinase K increases the global conformational flexibility, decreases the local flexibility of substrate-binding region, and does not influence the thermal motion of catalytic triad, thus explaining the experimentally determined decreased thermal stability, reduced substrate affinity, and almost unchanged catalytic activity upon Ca²⁺ removal; (iii) substrate binding affects the large concerted motions of proteinase K, and the resulting dynamic pocket can be connected to substrate binding, orientation, and product release; (iv) amino acid mutations 375 S/W and 423 I/P of HIV-1 gp120 have distinct effects on molecular motions of gp120, facilitating 375 S/W mutant to assume the CD4-bound conformation, while 423 I/P mutant to prefer for CD4-unliganded state. The mechanisms underlying protein dynamics and protein–ligand binding, including the concept of the free energy landscape (FEL) of the protein–solvent system, how the ruggedness and variability of FEL determine protein’s dynamics, and how the three ligand-binding models, the lock-and-key, induced-fit, and conformational selection are rationalized based on the FEL theory are discussed in depth.     

Investigation of Domain Motions in Bacteriophage T4 Lysozyme

Abstract

Hinge-bending in T4 lysozyme has been inferred from single amino acid mutant crystalline allomorphs by Matthews and coworkers. This raises an important question: are the different conformers in the unit cell artifacts of crystal packing forces, or do they represent different solution state structures? The objective of this theoretical study is to determine whether domain motions and hinge-bending could be simulated in T4 lysozyme using molecular dynamics. An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. Molecular dynamics calculations were computed using the Discover software package (Biosym Technologies). All hydrogen atoms were modeled explicitly with the inclusion of all 152 crystallographic waters at a temperature of 300 K. The native T4 lysozyme molecular dynamics simulation demonstrated hinge-bending in the protein. Relative domain motions between the N-terminal and C-terminal domains were evident. The enzyme hinge bending sites resulted from small changes in backbone atom conformations over several residues rather than rotation about a single bound. Two hinge loci were found in the simulation. One locus comprises residues 8–14 near the C-terminal of the A helix; the other site, residues 77–83 near the C-terminal of the C helix. Comparison of several snapshot structures from the dynamics trajectory clearly illustrates domain motions between the two lysozyme lobes. Time correlated atomic motions in the protein were analyzed using a dynamical cross-correlation map. We found a high degree of correlated atomic motions in each of the domains and, to a lesser extent, anticorrelated motions between the two domains. We also found that the hairpin loop in the N-terminal lobe (residues 19–24) acted as a mobile ‘flap’ and exhibited highly correlated dynamic motions across the cleft of the active site, especially with residue 142.     

Asymmetric structural motions of the homomeric alpha7 nicotinic receptor ligand binding domain revealed by molecular dynamics simulation

A homology model of the ligand binding domain of the alpha7 nicotinic receptor is constructed based on the acetylcholine-binding protein crystal structure. This structure is refined in a 10 ns molecular dynamics simulation. The modeled structure proves fairly resilient, with no significant changes at the secondary or tertiary structural levels. The hypothesis that the acetylcholine-binding protein template is in the activated or desensitized state, and the absence of a bound agonist in the simulation suggests that the structure may also be relaxing from this state to the activatable state. Candidate motions that take place involve not only the side chains of residues lining the binding sites, but also the subunit positions that determine the overall shape of the receptor. In particular, two nonadjacent subunits move outward, whereas their partners counterclockwise to them move inward, leading to a marginally wider interface between themselves and an overall asymmetric structure. This in turn affects the binding sites, producing two that are more open and characterized by distinct side-chain conformations of W54 and L118, although motions of the side chains of all residues in every binding site still contribute to a reduction in binding site size, especially the outward motion of W148, which hinders acetylcholine binding. The Cys loop at the membrane interface also displays some flexibility. Although the short simulation timescale is unlikely to sample adequately all the conformational states, the pattern of observed motions suggests how ligand binding may correlate with larger-scale subunit motions that would connect with the transmembrane region that controls the passage of ions. Furthermore, the shape of the asymmetry with binding sites of differing affinity for acetylcholine, characteristic of other nicotinic receptors, may be a natural property of the relaxed, activatable state of alpha7.     

Dynamical differences of hemoglobin and the ionotropic glutamate receptor in different states revealed by a new dynamics alignment method

A new algorithm for comparison of protein dynamics is presented. Compared protein structures are superposed and their modes of motions are calculated using the anisotropic network model. The obtained modes are aligned using the dynamic programming algorithm of Needleman and Wunsch, commonly used for sequence alignment. Dynamical comparison of hemoglobin in the T and R2 states reveals that the dynamics of the allosteric effector 2,3‐bisphosphoglycerate binding site is different in the two states. These differences can contribute to the selectivity of the effector to the T state. Similar comparison of the ionotropic glutamate receptor in the kainate+(R,R)‐2b and ZK bound states reveals that the kainate+(R,R)‐2b bound states slow modes describe upward motions of ligand binding domain and the transmembrane domain regions. Such motions may lead to the opening of the receptor. The upper lobes of the LBDs of the ZK bound state have a smaller interface with the amino terminal domains above them and have a better ability to move together. The present study exemplifies the use of dynamics comparison as a tool to study protein function. Proteins 2017; 85:1507–1517. © 2014 Wiley Periodicals, Inc.     

Conformational and dynamics changes induced by bile acids binding to chicken liver bile acid binding protein

The correlation between protein motions and function is a central problem in protein science. Several studies have demonstrated that ligand binding and protein dynamics are strongly correlated in intracellular lipid binding proteins (iLBPs), in which the high degree of flexibility, principally occurring at the level of helix-II, CD, and EF loops (the so-called portal area), is significantly reduced upon ligand binding. We have recently investigated by NMR the dynamic properties of a member of the iLBP family, chicken liver bile acid binding protein (cL-BABP), in its apo and holo form, as a complex with two bile salts molecules. Binding was found to be regulated by a dynamic process and a conformational rearrangement was associated with this event. We report here the results of molecular dynamics (MD) simulations performed on apo and holo cL-BABP with the aim of further characterizing the protein regions involved in motion propagation and of evaluating the main molecular interactions stabilizing bound ligands. Upon binding, the root mean square fluctuation values substantially decrease for CD and EF loops while increase for the helix-loop-helix region, thus indicating that the portal area is the region mostly affected by complex formation. These results nicely correlate with backbone dynamics data derived from NMR experiments. Essential dynamics analysis of the MD trajectories indicates that the major concerted motions involve the three contiguous structural elements of the portal area, which however are dynamically coupled in different ways whether in the presence or in the absence of the ligands. Motions of the EF loop and of the helical region are part of the essential space of both apo and holo-BABP and sample a much wider conformational space in the apo form. Together with NMR results, these data support the view that, in the apo protein, the flexible EF loop visits many conformational states including those typical of the holo state and that the ligand acts stabilizing one of these pre-existing conformations. The present results, in agreement with data reported for other iLBPs, sharpen our knowledge on the binding mechanism for this protein family.     

Global changes in local protein dynamics reduce the entropic cost of carbohydrate binding in the arabinose-binding protein

Protein dynamics make important but poorly understood contributions to molecular recognition phenomena. To address this, we measure changes in fast protein dynamics that accompany the interaction of the arabinose-binding protein (ABP) with its ligand, d-galactose, using NMR relaxation and molecular dynamics simulation. These two approaches present an entirely consistent view of the dynamic changes that occur in the protein backbone upon ligand binding. Increases in the amplitude of motions are observed throughout the protein, with the exception of a few residues in the binding site, which show restriction of dynamics. These counter-intuitive results imply that a localised binding event causes a global increase in the extent of protein dynamics on the pico- to nanosecond timescale. This global dynamic change constitutes a substantial favourable entropic contribution to the free energy of ligand binding. These results suggest that the structure and dynamics of ABP may be adapted to exploit dynamic changes to reduce the entropic costs of binding.     

Molecular simulations of β‐lactoglobulin complexed with fatty acids reveal the structural basis of ligand affinity to internal and possible external binding sites

The interaction of saturated fatty acids of different length (C8:0 to C18:0) with β‐lactoglobulin (βLG) was investigated by molecular dynamics simulation and docking approaches. The results show that the presence of such ligands in the hydrophobic central cavity of βLG, known as the protein calyx, determines an enhancement of atomic fluctuations compared with the unliganded form, especially for loops at the entrance of the binding site. Concerted motions are evidenced for protein regions that could favor the binding of ligands. The mechanism of anchoring of fatty acids of different length is similar for the carboxylate head‐group, through electrostatic interactions with the side chains of Lys60/Lys69. The key protein residues to secure the hydrocarbon chain are Phe105/Met107, which adapt their conformation upon ligand binding. In particular, Phe105 provides an additional hydrophobic clamp only for the tail of the two fatty acids with the longest chains, palmitic, and stearic acid, which are known to bind βLG with a high affinity. The search of additional external binding sites for fatty acids, distinct from the calyx, was also carried out for palmitic acid. Two external sites with a lower affinity were identified as secondary sites, one consisting in a hydrophobic cavity allowing two distinct binding modes for the fatty acid, and the other corresponding to a surface crevice close to the protein α‐helix. The overall results provide a comprehensive picture of the dynamical behavior of βLG in complex with fatty acids, and elucidate the structural basis of the binding of these physiological ligands. Proteins 2014; 82:2609–2619. © 2014 Wiley Periodicals, Inc.     

Relations between vitamin D and fatty acid binding properties of vitamin D-binding protein

The interaction of fatty acids with bovine vitamin D-binding protein (DBP) was studied using a partition equilibrium method. This protein has one high affinity site for binding of fatty acids with an association constant Ka = 7 x 10(5) M-1 for palmitic acid and Ka = 6 x 10(5) M-1 for arachidonic acid. Competition experiments showed that palmitic acid hardly competes with 25-hydroxycholecalciferol for binding to DBP. However, arachidonic acid showed comparatively a stronger competition for binding to this protein. The great difference in competition of palmitic and arachidonic acids with 25-hydroxycholecalciferol may be related to changes in DBP conformation promoted by the binding of different ligands.     

MptpA Kinetics Enhanced by Allosteric Control of an Active Conformation

Understanding allostery in the Mycobacterium tuberculosis low molecular weight protein tyrosine phosphatase (MptpA) is a subject of great interest since MptpA is one of two protein tyrosine phosphatases (PTPs) from the pathogenic organism Mycobacterium tuberculosis expressed during host cell infection. Here, we combine computational modeling with solution NMR spectroscopy and we find that Q75 is an allosteric site. Removal of the polar side chain of Q75 by mutation to leucine results in a cascade of events that reposition the acid loop over the active site and relocates the catalytic aspartic acid (D126) at an optimal position for proton donation to the leaving aryl group of the substrate and for subsequent hydrolysis of the thiophosphoryl intermediate. The computational analysis is consistent with kinetic data, and NMR spectroscopy, showing that the Q75L mutant exhibits enhanced reaction kinetics with similar substrate binding affinity. We anticipate that our findings will motivate further studies on the possibility that MptpA remains passivated during the chronic state of infection and increases its activity as part of the pathogenic life cycle of M. tuberculosis possibly via allosteric means.     

X-ray absorption and molecular dynamics study of cation binding sites in the purple membrane

The present work describes the results of a study aimed at identifying candidate cation binding sites on the extracellular region of bacteriorhodopsin, including a site near the retinal pocket. The approach used is a combined effort involving computational chemistry methods (computation of cation affinity maps and molecular dynamics) together with the Extended X-Ray Absorption Fine Structure (EXAFS) technique to obtain relevant information about the local structure of the protein in the neighborhood of Mn(2+) ions in different affinity binding sites. The results permit the identification of a high-affinity binding site where the ion is coordinated simultaneously to Asp212(-) and Asp85(-). Comparison of EXAFS data of the wild type protein with the quadruple mutant E9Q/E74Q/E194Q/E204Q at pH 7.0 and 10.0 demonstrate that extracellular glutamic acid residues are involved in cation binding.     

Global and local motions in ribonuclease A: a molecular dynamics study

The understanding of protein dynamics is one of the major goals of structural biology. A direct link between protein dynamics and function has been provided by x-ray studies performed on ribonuclease A (RNase A) (B. F. Rasmussen et al., Nature, 1992, Vol. 357, pp. 423-424; L. Vitagliano et al., Proteins: Structure, Function, and Genetics, 2002, Vol. 46, pp. 97-104). Here we report a 3 ns molecular dynamics simulation of RNase A in water aimed at characterizing the dynamical behavior of the enzyme. The analysis of local and global motions provides interesting insight on the dynamics/function relationship of RNase A. In agreement with previous crystallographic reports, the present study confirms that the RNase A active site is constituted by rigid (His12, Asn44, Thr45) and flexible (Lys41, Asp83, His119, Asp121) residues. The analysis of the global motions, performed using essential dynamics, shows that the two beta-sheet regions of RNase A move coherently in opposite directions, thus modifying solvent accessibility of the active site, and that the mixed alpha/3(10)-helix (residues 50-60) behaves as a mechanical hinge during the breathing motion of the protein. These data demonstrate that this motion, essential for RNase A substrate binding and release, is an intrinsic dynamical property of the ligand-free enzyme.     

ATP binding enables broad antibiotic selectivity of aminoglycoside phosphotransferase(3')-IIIa: an elastic network analysis

The bacterial enzyme aminoglycoside phosphotransferase(3′)-IIIa (APH) confers resistance against a wide range of aminoglycoside antibiotics. In this study, we use the Gaussian network model to investigate how the binding of nucleotides and antibiotics influences the dynamics and thereby the ligand binding properties of APH. Interestingly, in NMR experiments, the dynamics differ significantly in various APH complexes, although crystallographic studies indicate that no larger conformational changes occur upon ligand binding. Isothermal titration calorimetry also shows different thermodynamic contributions to ligand binding. Formation of aminoglycoside-APH complexes is enthalpically driven, while the enthalpic change upon aminoglycoside binding to the nucleotide-APH complex is much smaller. The differential effects of nucleotide binding and antibiotic binding to APH can be explained theoretically by single-residue fluctuations and correlated motions of the enzyme. The surprising destabilization of β-sheet residues upon nucleotide binding, as seen in hydrogen/deuterium exchange experiments, shows that the number of closest neighbors does not fully explain residue flexibility. Additionally, we must consider correlated motions of dynamic protein domains, which show that not only connectivity but also the overall protein architecture is important for protein dynamics.     

Dynamics based alignment of proteins: an alternative approach to quantify dynamic similarity

Background  

The dynamic motions of many proteins are central to their function. It therefore follows that the dynamic requirements of a protein are evolutionary constrained. In order to assess and quantify this, one needs to compare the dynamic motions of different proteins. Comparing the dynamics of distinct proteins may also provide insight into how protein motions are modified by variations in sequence and, consequently, by structure. The optimal way of comparing complex molecular motions is, however, far from trivial. The majority of comparative molecular dynamics studies performed to date relied upon prior sequence or structural alignment to define which residues were equivalent in 3-dimensional space.     

Oncogenic mutations on Rac1 affect global intrinsic dynamics underlying GTP and PAK1 binding

Rac1 is a small member of the Rho GTPase family. One of the most important downstream effectors of Rac1 is a serine/threonine kinase, p21-activated kinase 1 (PAK1). Mutational activation of PAK1 by Rac1 has oncogenic signaling effects. Here, although we focus on Rac1-PAK1 interaction by atomic-force-microscopy-based single-molecule force spectroscopy experiments, we explore the effect of active mutations on the intrinsic dynamics and binding interactions of Rac1 by Gaussian network model analysis and molecular dynamics simulations. We observe that Rac1 oncogenic mutations are at the hinges of three global modes of motion, suggesting the mechanical changes as potential markers of oncogenicity. Indeed, the dissociation of wild-type Rac1-PAK1 complex shows two distinct unbinding dynamic states that are reduced to one with constitutively active Q61L and oncogenic Y72C mutant Rac1, as revealed by single-molecule force spectroscopy experiments. Q61L and Y72C mutations change the mechanics of the Rac1-PAK1 complex by increasing the elasticity of the protein and slowing down the transition to the unbound state. On the other hand, Rac1’s intrinsic dynamics reveal more flexible GTP and PAK1-binding residues on switches I and II with Q61L, Y72C, oncogenic P29S and Q61R, and negative T17N mutations. The cooperativity in the fluctuations of GTP-binding sites around the p-loop and switch I decreases in all mutants, mostly in Q61L, whereas some PAK1-binding residues display enhanced coupling with GTP-binding sites in Q61L and Y72C and within each other in P29S. The predicted binding free energies of the modeled Rac1-PAK1 complexes show that the change in the dynamic behavior likely means a more favorable PAK1 interaction. Overall, these findings suggest that the active mutations affect intrinsic functional dynamic events and alter the mechanics underlying the binding of Rac1 to GTP and upstream and downstream partners including PAK1.     

Analysis of the relationship between enzyme activity and its internal motion using nuclear magnetic resonance: 15N relaxation studies of wild-type and mutant lysozyme

A mutant lysozyme where R14 and H15 are deleted together has higher activity and a similar binding ability to an inhibitor, trimer of N-acetylglucosamine ((NAG)3), compared with wild-type lysozyme. Since this has been attributed to intrinsic protein dynamic properties, we investigated the relationship between the activity and the internal motions of proteins. Backbone dynamics of the free and the complex forms with the (NAG)3 have been studied by measurement of the 15N T1 and T2 relaxation rates and NOE determinations at 600 MHz. Analysis of the data using the model-free formalism showed that the generalized order parameters (S2) were almost the same in wild-type and mutant lysozyme in unbound state, indicating that the mutation had little effect on the global internal motions. On the other hand, in the presence of (NAG)3, although some signals located around the active site were broadened or decreased in intensity because of strong perturbation by (NAG)3, there were several residues that showed increased or decreased backbone S2 in the complexed lysozymes. A comparison of the internal motions of the wild-type and mutant complexes showed a number of distinct dynamic differences between them. In particular, many residues located at or near active-site regions (turn 1, strand 2, turn 2 and long loop), displayed greater backbone dynamics reflecting the order parameter in mutant complex relative to mutant free. Furthermore, the Rex values at the loop C-D region, which was considered to be important for enzymatic activity, significantly increased. From these results, it was suggested that variations in the dynamics of these regions may play an important role in the enzyme activity.     

Molecular dynamics studies of caspase-3

Caspase-3 is a fundamental target for pharmaceutical interventions against a variety of diseases involving disregulated apoptosis. The enzyme is active as a dimer with two symmetry-related active sites, each featuring a Cys-His catalytic dyad and a selectivity loop, which recognizes the characteristic DEVD pattern of the substrate. Here, a molecular dynamics study of the enzyme in complex with two pentapeptide substrates DEVDG is presented, which provides a characterization of the dynamic properties of the active form in aqueous solution. The mobility of the substrate and that of the catalytic residues are rather low indicating a distinct preorganization effect of the Michaelis complex. An essential mode analysis permits us to identify coupled motions between the two monomers. In particular, it is found that the motions of the two active site loops are correlated and tend to steer the substrate toward the reactive center, suggesting that dimerization has a distinct effect on the dynamic properties of the active site regions. The selectivity loop of one monomer turns out to be correlated with the N-terminal region of the p12 subunit of the other monomer, an interaction that is also found to play a fundamental role in the electrostatic stabilization of the quaternary structure. To further characterize the specific influence of dimerization on the enzyme essential motions, a molecular dynamics analysis is also performed on the isolated monomer.