Elucidation of the conformational free energy landscape in H.pylori LuxS and its implications to catalysis

Background: One of the major challenges in understanding enzyme catalysis is to identify the different conformations and their populations at detailed molecular level in response to ligand binding/environment. A detail description of the ligand induced conformational changes provides meaningful insights into the mechanism of action of enzymes and thus its function. Results: In this study, we have explored the ligand induced conformational changes in H. pylori LuxS and the associated mechanistic features. LuxS, a dimeric protein, produces the precursor (4,5-dihydroxy-2,3-pentanedione) for autoinducer-2 production which is a signalling molecule for bacterial quorum sensing. We have performed molecular dynamics simulations on H. pylori LuxS in its various ligand bound forms and analyzed the simulation trajectories using various techniques including the structure network analysis, free energy evaluation and water dynamics at the active site. The results bring out the mechanistic details such as co operativity and asymmetry between the two subunits, subtle changes in the conformation as a response to the binding of active and inactive forms of ligands and the population distribution of different conformations in equilibrium. These investigations have enabled us to probe the free energy landscape and identify the corresponding conformations in terms of network parameters. In addition, we have also elucidated the variations in the dynamics of water co-ordination to the Zn2+ ion in LuxS and its relation to the rigidity at the active sites. Conclusions: In this article, we provide details of a novel method for the identification of conformational changes in the different ligand bound states of the protein, evaluation of ligand-induced free energy changes and the biological relevance of our results in the context of LuxS structure-function. The methodology outlined here is highly generalized to illuminate the linkage between structure and function in any protein of known structure.

Network approach for capturing ligand-induced subtle global changes in protein structures

Ligand-induced conformational changes in proteins are of immense functional relevance. It is a major challenge to elucidate the network of amino acids that are responsible for the percolation of ligand-induced conformational changes to distal regions in the protein from a global perspective. Functionally important subtle conformational changes (at the level of side-chain noncovalent interactions) upon ligand binding or as a result of environmental variations are also elusive in conventional studies such as those using root-mean-square deviations (r.m.s.d.s). In this article, the network representation of protein structures and their analyses provides an efficient tool to capture these variations (both drastic and subtle) in atomistic detail in a global milieu. A generalized graph theoretical metric, using network parameters such as cliques and/or communities, is used to determine similarities or differences between structures in a rigorous manner. The ligand-induced global rewiring in the protein structures is also quantified in terms of network parameters. Thus, a judicious use of graph theory in the context of protein structures can provide meaningful insights into global structural reorganizations upon perturbation and can also be helpful for rigorous structural comparison. Data sets for the present study include high-resolution crystal structures of serine proteases from the S1A family and are probed to quantify the ligand-induced subtle structural variations.

Dramatic Structural Changes Resulting from the Loss of a Crucial Hydrogen Bond in the Hinge Region Involved in C-Terminal Helix Swapping in SurE: A Survival Protein from Salmonella typhimurium

Domain swapping is an interesting feature of some oligomeric proteins in which each protomer of the oligomer provides an identical surface for exclusive interaction with a segment or domain belonging to another protomer. Here we report results of mutagenesis experiments on the structure of C-terminal helix swapped dimer of a stationary phase survival protein from Salmonella typhimurium (StSurE). Wild type StSurE is a dimer in which a large helical segment at the C-terminus and a tetramerization loop comprising two beta strands are swapped between the protomers. Key residues in StSurE that might promote C-terminal helix swapping were identified by sequence and structural comparisons. Three mutants in which the helix swapping is likely to be avoided were constructed and expressed in E. coli. Three-dimensional X-ray crystal structures of the mutants H234A and D230A/H234A could be determined at 2.1 angstrom and 2.35 angstrom resolutions, respectively. Contrary to expectations, helix swapping was mostly retained in both the mutants. The loss of the crucial D230 OD2- H234 NE2 hydrogen bond (2.89 angstrom in the wild type structure) in the hinge region was compensated by new inter and intra-chain interactions. However, the two fold molecular symmetry was lost and there were large conformational changes throughout the polypeptide. In spite of these changes, the dimeric structure and an approximate tetrameric organization were retained, probably due to the interactions involving the tetramerization loop. Mutants were mostly functionally inactive, highlighting the importance of precise inter-subunit interactions for the symmetry and function of StSurE.

Switch Region for Pathogenic Structural Change in Conformational Disease and Its Prediction

Many diseases are believed to be related to abnormal protein folding. In the first step of such pathogenic structural changes, misfolding occurs in regions important for the stability of the native structure. This destabilizes the normal protein conformation, while exposing the previously hidden aggregation-prone regions, leading to subsequent errors in the folding pathway. Sites involved in this first stage can be deemed switch regions of the protein, and can represent perfect binding targets for drugs to block the abnormal folding pathway and prevent pathogenic conformational changes. In this study, a prediction algorithm for the switch regions responsible for the start of pathogenic structural changes is introduced. With an accuracy of 94%, this algorithm can successfully find short segments covering sites significant in triggering conformational diseases (CDs) and is the first that can predict switch regions for various CDs. To illustrate its effectiveness in dealing with urgent public health problems, the reason of the increased pathogenicity of H5N1 influenza virus is analyzed; the mechanisms of the pandemic swine-origin 2009 A(H1N1) influenza virus in overcoming species barriers and in infecting large number of potential patients are also suggested. It is shown that the algorithm is a potential tool useful in the study of the pathology of CDs because: (1) it can identify the origin of pathogenic structural conversion with high sensitivity and specificity, and (2) it provides an ideal target for clinical treatment.

Real-Time Study of Genomic DNA Structural Changes upon Interaction with Small Molecules Using Dual-Polarization Interferometry

The structural changes of genomic DNA upon interaction with small molecules have been studied in real time using dual-polarization interferometry (DPI). Native or thermally denatured DNA was immobilized on the silicon oxynitride surface via a preadsorbed poly(ethylenimine) (PEI) layer. The mass loading was similar for both types of DNA, however, native DNA formed a looser and thicker layer due to its rigidity, unlike the more flexible denatured DNA, which mixed with PEI to form a denser and thinner layer. Ethidium bromide (EtBr), a classical intercalator, induced the large thickness decrease and density increase of native DNA (double-stranded), but a slight increase in both the thickness and density of denatured DNA (partial single-stranded).

Conformational transition of DNA induced by cationic lipid vesicle in acidic solution: spectroscopy investigation

The conformational transition of DNA induced by the interaction between DNA and a cationic lipid vesicle, didodecyidimethylammonium bromide (DDAB), had been investigated by circular dichroism (CD) and UV spectroscopy methods. We used singular value decomposition least squares method (SVDLS) to analyze the experimental CD spectra. Although pH value influenced the conformation of DNA in solution, the results showed that upon binding to double helical DNA, positively charged liposomes induced a conformational transition of DNA molecules from the native B-form to more compact conformations. At the same time, no obvious conformational changes occurred at single-strand DNA (ssDNA). While the cationic lipid vesicles and double-strand DNA (dsDNA) were mixed at a high molar ratio of DDAB vesicles to dsDNA, the conformation of dsDNA transformed from the B-form to the C-form resulting in an increase in duplex stability (DeltaT(m) = 8 +/- 0.4 degreesC). An increasing in T-m was also observed while the cationic lipid vesicles interacted with ssDNA.

Characterization of the structural and functional changes of hemoglobin in dimethyl sulfoxide by spectroscopic techniques

Circular dichroism (CD), fourier transform infrared (FTIR), and fluorescence spectroscopy were used to explore the effect of dimethyl sulfoxide (DMSO) on the structure and function of hemoglobin (Hb). The native tertiary structure was disrupted completely when the concentration of DMSO reached 50% (v/v), which was determined by loss of the characteristic Soret CD spectrum. Loss of the native tertiary structure could be mainly caused by breaking the hydrogen bonds, between the heme propionate groups and nearby surface amino acid residues, and by disorganizing the hydrophobic interior of this protein. Upon exposure of Hb to 52% DMSO for ca. 12 h in a D2O medium no significant change in 1652 cm(-1) band of the FTIR spectrum was produced, which demonstrated that alpha-helical structure predominated. When the concentration of DMSO increased to 57%: (1) the band at 1652 cm(-1) disappeared with the appearance of two new bands located at 1661 and 1648 cm(-1); (2) another new band at 1623 cm(-1) was attributed to the formation of intermolecular beta-sheet or aggregation, which was the direct consequence of breaking of the polypeptide chain by the competition of S=O groups in DMSO with C=O groups in amide bonds. Further increasing the DMSO concentration to 80%, the intensity at 1623 cm(-1) increased, and the bands at 1684, 1661 and 1648 cm(-1) shifted to 1688, 1664 and 1644 cm(-1), respectively. These changes showed that the native secondary structure of Hb was last and led to further aggregation and increase of the content of 'free' amide C=O groups. In pure DMSO solvent, the major band at 1664 cm(-1) indicated that almost all of both the intermolecular beta-sheet and any residual secondary structure were completely disrupted. The red shift of the fluorescence emission maxima showed that the tryptophan residues were exposed to a greater hydrophilic environment as the DMSO content increased. GO-binding experiment suggested that the biological function of Hb was disrupted seriously even if the content of DMSO was 20%. (C) 1998 Elsevier Science B.V. All rights reserved.

Conformational kinetics reveals affinities of protein conformational states.

Most biological reactions rely on interplay between binding and changes in both macromolecular structure and dynamics. Practical understanding of this interplay requires detection of critical intermediates and determination of their binding and conformational characteristics. However, many of these species are only transiently present and they have often been overlooked in mechanistic studies of reactions that couple binding to conformational change. We monitored the kinetics of ligand-induced conformational changes in a small protein using six different ligands. We analyzed the kinetic data to simultaneously determine both binding affinities for the conformational states and the rate constants of conformational change. The approach we used is sufficiently robust to determine the affinities of three conformational states and detect even modest differences in the protein's affinities for relatively similar ligands. Ligand binding favors higher-affinity conformational states by increasing forward conformational rate constants and/or decreasing reverse conformational rate constants. The amounts by which forward rate constants increase and reverse rate constants decrease are proportional to the ratio of affinities of the conformational states. We also show that both the affinity ratio and another parameter, which quantifies the changes in conformational rate constants upon ligand binding, are strong determinants of the mechanism (conformational selection and/or induced fit) of molecular recognition. Our results highlight the utility of analyzing the kinetics of conformational changes to determine affinities that cannot be determined from equilibrium experiments. Most importantly, they demonstrate an inextricable link between conformational dynamics and the binding affinities of conformational states.

Metal ions and protein folding: conformational and functional interplay

Dissertation presented to obtain a PhD degree in Biochemistry at Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa

Hydration, conformational states and kinetics of yeast hexokinase PII /

The kinetic study of the coupled enzymatic reaction involving monomeric yeast hexokinase PII (HK) and yeast glucose-6-phosphate dehydrogenase (G-6-PDH) yields a Michaelis constant of 0.15 ± 0.01 mM for D-glucose. At pH 8.7 HK is present in monomeric form. The addition of polyethylene glycol (PEG), to the reaction mixture increased the affinity of HK for glucose, independent ofMW of the PEG from 2000 to 10000. The osmotic stress exerted by PEG can be used to measure the change in number of water molecules that accompany enzyme conformational changes (Rand, et al., 1993). Results indicate that the G-6-PDH is not osmotically sensitive and thus, the change in the number of PEG-inaccessible water molecules (ANw) measured in the coupled reaction is only the difference between the glucose-bound and glucosefree conformations of HK. ANw ~ 450 with PEGs of MW > 2000 under conditions for both binding (Reid and Rand, 1997) and kinetic assays. The contribution water may play in the binding of ATP (Km = 0.24 + 0.02 mM) has also been examined. It was found that in this case ANw = (for osmotic pressures < 2.8x10* dynes/cm^), suggesting no additional numbers of waters are displaced when ATP binds to HK. Osmotic pressure experiments were also performed with dimeric HK. It was determined that both the monomeric and dimeric forms of HK give the same ANw under low pressures. If this large ANw is due to conformational flexibility, it would appear that the flexibility is not reduced upon dimerization ofthe enzyme.

Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

Conformational Properties of Angiotensin II and Its Active and Inactive TOAC-Labeled Analogs in the Presence of Micelles. Electron Paramagnetic Resonance, Fluorescence, and Circular Dichroism Studies

The interaction between angiotensin II (AII, DRVYIHPF) and its analogs carrying 2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4-carboxylic acid (TOAC) and detergents-negatively charged sodium dodecyl sulfate (SDS) and zwitterionic N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS)-was examined by means of EPR, CD, and fluorescence. EPR spectra of partially active TOAC(1)-AII and inactive TOAC(3)-AII in aqueous solution indicated fast tumbling, the freedom of motion being greater at the N-terminus. Line broadening occurred upon interaction with micelles. Below SDS critical micelle concentration, broader lines indicated complex formation with tighter molecular packing than in micelles. Small changes in hyperfine splittings evinced TOAC location at the micelle-water interface. The interaction with anionic micelles was more effective than with zwitterionic micelles. Peptide-micelle interaction caused fluorescence increase. The TOAC-promoted intramolecular fluorescence quenching was more, pronounced for TOAC(3)-AII because of the proximity between the nitroxide and Tyr(4). CD spectra showed that although both AII and TOAC(1)-AII presented flexible conformations in water, TOAC(3)-AII displayed conformational restriction because of the TOAC-imposed bend (Schreier et al., Biopolymers 2004, 74, 389). In HPS, conformational changes were observed for the labeled peptides at neutral and basic pH. In SDS, all peptides underwent pH-dependent conformational changes. Although the spectra suggested similar folds for All and TOAC(1)-AII, different conformations were acquired by TOAC(3)-AII. The membrane environment has been hypothesized to shift conformational equilibria so as to stabilize the receptor-bound conformation of ligands. The fact that TOAC(3)-AII is unable to acquire conformations similar to those of native AII and partially active TOAC(1)-AII is probably the explanation for its lack of biological activity. (C) 2009 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 92: 525-537, 2009.

A molecular mechanism for Lys(49)-phospholipase A(2) activity based on ligand-induced conformational change

Agkistrodon contortrix laticinctus myotoxin is a Lys(49)- phospholipase A(2) (EC 3.1.1.4) isolated from the venom of the serpent A contortrix laticinctus (broad-banded copperhead). We present here three monomeric crystal structures of the myotoxin, obtained under different crystallization conditions. The three forms present notable structural differences and reveal that the presence of a ligand in the active site (naturally presumed to be a fatty acid) induces the exposure of a hydrophobic surface (the hydrophobic knuckle) toward the C terminus. The knuckle in A contortrix laticinctus myotoxin involves the side chains of Phe(121) and Phe(124) and is a consequence of the formation of a canonical structure for the main chain within the region of residues 118-125. Comparison with other Lys(49)-phospholipase A(2) myotoxins shows that although the knuckle is a generic structural motif common to all members of the family, it is not readily recognizable by simple sequence analyses. An activation mechanism is proposed that relates fatty acid retention at the active site to conformational changes within the C-terminal region, a part of the molecule that has long been associated with Ca2+-independent membrane damaging activity and myotoxicity. This provides, for the first time, a direct structural connection between the phospholipase active site and the C-terminal myotoxic site, justifying the otherwise enigmatic conservation of the residues of the former in supposedly catalytically inactive molecules.

Whydration effects on DNA double helix stability modulates ligand binding to natural DNA in response to changes in water activity

In this work we present evidence that water molecules are actively involved on the control of binding affinity and binding site discrimination of a drug to natural DNA. In a previous study, the effect of water activity (a(w)) on the energetic parameters of actinomycin-D intercalation to natural DNA was determined using the osmotic stress method (39). This earlier study has shown evidence that water molecules act as an allosteric regulator of ligand binding to DNA via the effect of water activity on the long-range stability of the DNA secondary structure. In this work we have carried out DNA circularization experiments using the plasmid pUC18 in the absence of drugs and in the presence of different neutral solutes to evaluate the contribution of water activity to the energetics of DNA helix unwinding. The contribution of water to these independent reactions were made explicit by the description of how the changes in the free energy of ligand binding to DNA and in the free energy associated with DNA helix torsional deformation are linked to a(w) via changes in structural hydration. Taken together, the results of these studies reveal an extensive linkage between ligand binding affinity and site binding discrimination, and long range helix conformational changes and DNA hydration, This is strong evidence that water molecules work as a classical allosteric regulator of ligand binding to the DNA via its contribution to the stability of the double helix secondary structure, suggesting a possible mechanism by which the biochemical machinery of DNA processing takes advantage of the low activity of water into the cellular milieu.

Effect of Solvent-Induced Coil to Helix Conformational Change on the Two-Photon Absorption Spectrum of Poly(3,6-phenanthrene)

This paper investigates the effect of solvent-induced conformational changes of poly(3,6-phenanthrene) on their two-photon absorption (2PA). Such effect was studied employing the wavelength-tunable femtosecond Z-scan technique and modeled using the sum-over-essential states approach. We observed a strong reduction of the 2PA cross-section when the sample was prepared in hexane (poor solvent) in comparison to chloroform (good solvent), which is related to the conformation adopted by the polymer in each case. In chloroform it adopts a random coil conformation, as opposed to the one-handed helix conformation in hexane. Our results pointed out that the coil to helix conformation change decreases the degree of molecular planarity of the polymer pi-conjugated backbone, which is primarily responsible for their optical nonlinearity, contributing to diminishing the effective transition dipole moments and, consequently, the 2PA cross-section. Moreover, by studying the nonlinear response with different light polarization, we showed that, although the solvent-induced conformational change does not alter the molecular symmetry of the polymer, it modifies considerably the direction of the transition dipole moments between the excited states.