Search

Disinfectant biocides are chemicals that are heavily used for disinfection purposes in households, hospitals, and agrifood industry. The most common type of biocides are quaternary ammonium compounds (QAs), notably benzalkonium chloride (BAC) and didecyldimethylammonium chloride (DDAC), which have been shown to inhibit cholinesterases. This study aims to evaluate the effect of these biocides towards different cholinesterases using both enzyme inhibition and molecular docking experiments. Acetylcholinesterase (AChE) from Drosophila melanogaster (DM-AChE), Electrophorus electricus (EE-AChE), bovine erythrocytes (BE-AChE) and butyrylcholinesterase from horse serum (BChE) were selected for this study. Using a colorimetric assay, all these enzymes were shown to be inhibited in a competitive form by both biocides, BAC and DDAC, with the exception of DM-AChE, which was inhibited in a non-competitive manner by BAC. Molecular docking experiments enabled to identify structural determinants involved in the different modes of inhibition observed. More particularly, our results suggest that non-competitive inhibition of DM-AChE by BAC could be related to the binding of the inhibitor into a more extended active site compared to other cholinesterases., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

G-quadruplexes (G4) are helical structures found in guanine-rich DNA or RNA sequences. Generally, their formalism is based on a few dozen structures, which can produce some inconsistencies or incompleteness. Using the website ASC-G4, we analyzed the structures of 333 intramolecular G4s, of all types, which allowed us to clarify some key concepts and present new information. To each of the eight distinguishable topologies corresponds a groove-width signature and a predominant glycosidic configuration (gc) pattern governed by the directions of the strands. The relative orientations of the stacking guanines within the strands, which we quantified and related to their vertical gc successions, determine the twist and tilt of the helices. The latter impact the minimum groove widths, which represent the space available for lateral ligand binding. The G4 four helices have similar twists, even when these twists are irregular, meaning that they have various angles along the strands. Despite its importance, the vertical gc succession has no strict one-to-one relationship with the topology, which explains the discrepancy between some topologies and their corresponding circular dichroism spectra. This study allowed us to introduce the new concept of platypus G4s, which are structures with properties corresponding to several topologies., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

A sensitive and robust electrochemical cholinesterase-based sensor was developed to detect the quaternary ammonium (QAs) biocides most frequently found in agri-food industry wash waters: benzalkonium chloride (BAC) and didecyldimethylammonium chloride (DDAC). To reach the maximum residue limit of 28 nM imposed by the European Union (EU), two types of cholinesterases were tested, acetylcholinesterase (AChE, from Drosophila melanogaster ) and butyrylcholinesterase (BChE, from horse serum). The sensors were designed by entrapping AChE or BChE on cobalt phthalocyanine-modified screen-printed carbon electrodes. The limits of detection (LOD) of the resulting biosensors were 38 nM for DDAC and 320 nM for BAC, using, respectively, AChE and BChE. A simple solid-phase extraction step was used to concentrate the samples before biosensor analysis, allowing for the accurate determination of DDAC and BAC in tap water with limits of quantification (LOQ) as low as 2.7 nM and 5.3 nM, respectively. Additional assays demonstrated that the use of a phosphotriesterase enzyme allows for the total removal of interferences due to the possible presence of organophosphate insecticides in the sample. The developed biosensors were shown to be stable during 3 months storage at 4 °C.

ASC-G4 is an algorithm for the calculation of the advanced structural characteristics of G-quadruplexes (G4). It allows the unambiguous determination of the intramolecular G4 topology, based on the oriented strand numbering. It also resolves the ambiguity in the determination of the guanine glycosidic configuration. With this algorithm, we showed that the use of the C3' or C5' atoms to calculate the groove width in G4 is more appropriate than the P atoms and that the groove width does not always reflect the space available within the groove. For the latter, the minimum groove width is more appropriate. The application of ASC-G4 to 207 G4 structures guided the choices made for the calculations. A website based on ASC-G4 (http://tiny.cc/ASC-G4) was created, where the user uploads his G4 structure and gets its topology, the types of its loops and their lengths, the presence of snapbacks and bulges, the distribution of guanines in the tetrads and strands, the glycosidic configuration of these guanines, their rise, the groove widths, the minimum groove widths, the tilt and twist angles, the backbone dihedral angles, etc. It also provides a large number of atom-atom and atom-plane distances that are relevant to evaluating the quality of the structure., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

G-quadruplexes (G4s), secondary structures adopted by guanine-rich DNA and RNA sequences, are implicated in numerous biological processes and have been suggested as potential drug targets. Accordingly, there is an increasing interest in developing high-throughput methods that allow the generation of congeneric series of G4-targeting molecules ("ligands") and investigating their interactions with the targets. We have developed an operationally simple method of parallel synthesis to generate "ready-to-screen" libraries of cationic acylhydrazones, a motif that we have previously identified as a promising scaffold for potent, biologically active G4 ligands. Combined with well-established screening techniques, such as fluorescence melting, this method enables the rapid synthesis and screening of combinatorial libraries of potential G4 ligands. Following this protocol, we synthesized a combinatorial library of 90 bis(acylhydrazones) and screened it against five different nucleic acid structures. This way, we were able to analyze the structure-activity relationships within this series of G4 ligands, and identified three novel promising ligands whose interactions with G4-DNAs of different topologies were studied in detail by a combination of several biophysical techniques, including native mass spectrometry, and molecular modeling., (© 2022 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

SERCA1a is an ATPase calcium pump that transports Ca 2+ from the cytoplasm to the sarco/endoplasmic reticulum lumen. Sarcolipin (SLN), a transmembrane peptide, regulates the activity of SERCA1a by decreasing its Ca 2+ transport rate, but its mechanism of action is still not well-understood. To decipher this mechanism, we have performed normal mode analysis in the all-atom model, with the SERCA1a-SLN complex, or the isolated SERCA1a, embedded in an explicit membrane. The comparison of the results allowed us to provide an explanation at the atomic level for the action of SLN that is in good agreement with experimental observations. In our analyses, the presence of SLN locally perturbs the TM6 transmembrane helix and as a consequence modifies the position of D800, one of the key metal-chelating residues. Additionally, it reduces the flexibility of the gating residues, V304, and E309 in TM4, at the entrance of the Ca 2+ binding sites, which would decrease the affinity for Ca 2+ . Unexpectedly, SLN has also an effect on the ATP binding site more than 35 Å away, due to the straightening of TM5, a long helix considered as the spine of the protein. The straightening of TM5 modifies the structure of the P-N linker that sits above it, and which comprises the 351 DKTG 354 conserved motif, resulting in an increase of the distance between ATP and the phosphorylation site. As a consequence, the turn-over rate could be affected. All this gives SERCA1a the propensity to go toward a Ca 2+ low-affinity E2-like state in the presence of SLN and toward a Ca 2+ high-affinity E1-like state in the absence of SLN. In addition to a general mechanism of inhibition of SERCA1a regulatory peptides, this study also provides an insight into the conformational transition between the E2 and E1 states., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Barbot, Beswick, Montigny, Quiniou, Jamin and Mouawad.)

Acetylcholinesterase (AChE) is currently the most favorable target for the symptomatic treatment and reduction of Alzheimer's disease (AD). In order to identify new potent inhibitors of this enzyme, we describe herein a new structure-based virtual screening (SBVS) using the Institut Curie-CNRS chemical library (ICCL), which contained at the screening date 14307 compounds. The strategy undertaken in this work consisted of the use of several docking programs in SBVS calculations followed by the application of a consensus method (vSDC) and a scrupulous visual analysis. It allowed us to obtain a high degree of success, with a yield of almost 86 %, since 12 hits were identified among only 14 molecules tested in vitro. Still more remarkably, 6 of these hits were more active than galantamine, the reference inhibitor. These hits were predicted to have good ADMET properties. The two most promising compounds can serve as leads for AD treatment., (© 2019 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Background: In drug design, an efficient structure-based optimization of a ligand needs the precise knowledge of the protein-ligand interactions. In the absence of experimental information, docking programs are necessary for ligand positioning, and the choice of a reliable program is essential for the success of such an optimization. The performances of four popular docking programs, Gold, Glide, Surflex and FlexX, were investigated using 100 crystal structures of complexes taken from the Directory of Useful Decoys-Enhanced database., Results: The ligand conformational sampling was rather efficient, with a correct pose found for a maximum of 84 complexes, obtained by Surflex. However, the ranking of the correct poses was not as efficient, with a maximum of 68 top-rank or 75 top-4 rank correct poses given by Glidescore. No relationship was found between either the sampling or the scoring performance of the four programs and the properties of either the targets or the small molecules, except for the number of ligand rotatable bonds. As well, no exploitable relationship was found between each program performance in docking and in virtual screening; a wrong top-rank pose may obtain a good score that allows it to be ranked among the most active compounds and vice versa. Also, to improve the results of docking, the strengths of the programs were combined either by using a rescoring procedure or the United Subset Consensus (USC). Oddly, positioning with Surflex and rescoring with Glidescore did not improve the results. However, USC based on docking allowed us to obtain a correct pose in the top-4 rank for 87 complexes. Finally, nine complexes were scrutinized, because a correct pose was found by at least one program but poorly ranked by all four programs. Contrarily to what was expected, except for one case, this was not due to weaknesses of the scoring functions., Conclusions: We conclude that the scoring functions should be improved to detect the correct poses, but sometimes their failure may be due to other varied considerations. To increase the chances of success, we recommend to use several programs and combine their results. Graphical abstract Summary of the results obtained by semi-rigid docking of crystallographic ligands. The docking was done on 100 protein-ligand X-ray structures, taken from the DUD-E database, and using four programs, Glide, Gold, Surflex and FlexX. Based on the docking results, we applied our United Subset Consensus method (USC), for which only the top4-rank poses are relevant. The number of complexes for which the best pose is correct, is represented by the gray boxes, the blue and red boxes correspond to the number of complexes with a correct pose ranked as the top 1 or within the top 4. A pose is considered correct when its root-mean-square deviation from the crystal structure is less than 2 Å.

ICAT (Inhibitor of β-CAtenin and TCF) is a small acidic protein that negatively regulates β-catenin co-transcriptional activity by competing with TCF/LEF factors in their binding to β-catenin superhelical core. In melanoma cells, ICAT competes with LEF1 to negatively regulate the M-MITF and NEDD9 target genes. The structure of ICAT consists of two domains: the 3-helix bundle N-terminal domain binds to β-catenin Armadillo (Arm) repeats 10-12 and the C-terminal tail binds to Arm repeats 5-9. To elucidate the structural mechanisms governing ICAT/β-catenin interactions in melanoma cells, three ICAT residues Y15, K19 and V22 in the N-terminal domain, contacting hydrophobic β-catenin residue F660, were mutated and interaction was assessed by immunoprecipitation. Despite the moderate hydrophobicity of the contact, its removal completely abolished the interaction. In the ICAT C-terminal tail consensus sequence, neutralization of the electrostatic interactions between residues D66, E75 and β-catenin residues K435, K312, coupled to deletion of the hydrophobic contact between F71 and β-catenin R386, markedly reduced, but failed to abolish the ICAT-mediated negative regulation of M-MITF and NEDD9 promoters. We conclude that in melanoma cells, anchoring of ICAT N-terminal domain to β-catenin through the hook made by residue F660, trapped in the pincers formed by ICAT residues Y15 and V22, is crucial for stabilizing the ICAT/β-catenin complex. This is a prerequisite for binding of the consensus peptide to Arm repeats 5-9 and competition with LEF1. Differences between ICAT and LEF1 in their affinity for β-catenin may rely on the absence in ICAT of hydrophilic residues between D66 and F71.

Background: In a structure-based virtual screening, the choice of the docking program is essential for the success of a hit identification. Benchmarks are meant to help in guiding this choice, especially when undertaken on a large variety of protein targets. Here, the performance of four popular virtual screening programs, Gold, Glide, Surflex and FlexX, is compared using the Directory of Useful Decoys-Enhanced database (DUD-E), which includes 102 targets with an average of 224 ligands per target and 50 decoys per ligand, generated to avoid biases in the benchmarking. Then, a relationship between these program performances and the properties of the targets or the small molecules was investigated., Results: The comparison was based on two metrics, with three different parameters each. The BEDROC scores with α = 80.5, indicated that, on the overall database, Glide succeeded (score > 0.5) for 30 targets, Gold for 27, FlexX for 14 and Surflex for 11. The performance did not depend on the hydrophobicity nor the openness of the protein cavities, neither on the families to which the proteins belong. However, despite the care in the construction of the DUD-E database, the small differences that remain between the actives and the decoys likely explain the successes of Gold, Surflex and FlexX. Moreover, the similarity between the actives of a target and its crystal structure ligand seems to be at the basis of the good performance of Glide. When all targets with significant biases are removed from the benchmarking, a subset of 47 targets remains, for which Glide succeeded for only 5 targets, Gold for 4 and FlexX and Surflex for 2., Conclusion: The performance dramatic drop of all four programs when the biases are removed shows that we should beware of virtual screening benchmarks, because good performances may be due to wrong reasons. Therefore, benchmarking would hardly provide guidelines for virtual screening experiments, despite the tendency that is maintained, i.e., Glide and Gold display better performance than FlexX and Surflex. We recommend to always use several programs and combine their results. Graphical AbstractSummary of the results obtained by virtual screening with the four programs, Glide, Gold, Surflex and FlexX, on the 102 targets of the DUD-E database. The percentage of targets with successful results, i.e., with BDEROC(α = 80.5) > 0.5, when the entire database is considered are in Blue, and when targets with biased chemical libraries are removed are in Red.

Background: In drug design, one may be confronted to the problem of finding hits for targets for which no small inhibiting molecules are known and only low-throughput experiments are available (like ITC or NMR studies), two common difficulties encountered in a typical academic setting. Using a virtual screening strategy like docking can alleviate some of the problems and save a considerable amount of time by selecting only top-ranking molecules, but only if the method is very efficient, i.e. when a good proportion of actives are found in the 1-10 % best ranked molecules., Results: The use of several programs (in our study, Gold, Surflex, FlexX and Glide were considered) shows a divergence of the results, which presents a difficulty in guiding the experiments. To overcome this divergence and increase the yield of the virtual screening, we created the standard deviation consensus (SDC) and variable SDC (vSDC) methods, consisting of the intersection of molecule sets from several virtual screening programs, based on the standard deviations of their ranking distributions., Conclusions: SDC allowed us to find hits for two new protein targets by testing only 9 and 11 small molecules from a chemical library of circa 15,000 compounds. Furthermore, vSDC, when applied to the 102 proteins of the DUD-E benchmarking database, succeeded in finding more hits than any of the four isolated programs for 13-60 % of the targets. In addition, when only 10 molecules of each of the 102 chemical libraries were considered, vSDC performed better in the number of hits found, with an improvement of 6-24 % over the 10 best-ranked molecules given by the individual docking programs.Graphical abstractIn drug design, for a given target and a given chemical library, the results obtained with different virtual screening programs are divergent. So how to rationally guide the experimental tests, especially when only a few number of experiments can be made? The variable Standard Deviation Consensus (vSDC) method was developed to answer this issue. Left panel the vSDC principle consists of intersecting molecule sets, chosen on the basis of the standard deviations of their ranking distributions, obtained from various virtual screening programs. In this study Glide, Gold, FlexX and Surflex were used and tested on the 102 targets of the DUD-E database. Right panel Comparison of the average percentage of hits found with vSDC and each of the four programs, when only 10 molecules from each of the 102 chemical libraries of the DUD-E database were considered. On average, vSDC was capable of finding 38 % of the findable hits, against 34 % for Glide, 32 % for Gold, 16 % for FlexX and 14 % for Surflex, showing that with vSDC, it was possible to overcome the unpredictability of the virtual screening results and to improve them.

A microtubule is a dynamic system formed of αβ-tubulins. The presence of nonhydrolyzable guanosine-5'-triphosphate (GTP)/guanosine diphosphate (GDP) on the β-tubulins provokes microtubule polymerization/depolymerization. Despite the large number of experimental studies of this dynamical process, its mechanism is still unclear. To provide insights into this mechanism we studied the first depolymerization steps of GDP/GTP-bound microtubules by normal-mode analysis with the all-atom model. We also constructed a depolymerizing microtubule and compared it to cryo-electron microscopy tomograms (cyro-ET). The results show that during depolymerization, the protofilaments not only curve but twist to weaken their lateral interactions. These interactions are stabilized by GTP, but not evenly. Not all of the interface residues are of equal importance: five of them, belonging to the H2-S3 loop, play a special role; acting as a lock whose key is the γ-phosphate of GTP. Sequence alignments of several tubulins confirm the importance of these residues., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

Protein kinase casein kinase 2 (CK2) is a serine/threonine kinase with evidence of implication in growth dysregulation and apoptosis resistance, making it a relevant target for cancer therapy. Several CK2 inhibitors have been developed showing variable efficiency, emphasizing the need to expand the chemical diversity of those inhibitors. We report the identification and characterization of 2,8-difurandicarboxylic acid derivatives as a new class of nanomolar ATP-competitive inhibitors. Selectivity profiling pointed out proviral insertion Moloney virus kinases (Pim kinases) as the only other kinases that are significantly inhibited. By combining structure-activity relationship analysis with structural determination, we were able to determine the binding mode of these inhibitors for both kinases and to explain their strong inhibitory potency. Essential chemical features necessary for activity on both kinases were then identified. The described compounds are not cell permeable: however, they could provide a lead for developing novel inhibitors usable also in vivo. Given the similar but not redundant pathophysiological functions of CK2 and Pim family members, such inhibitors would provide new attractive leads for targeted cancer therapy. This work highlights that 2 functionally related kinases from different kinome branches display exquisite sensitivity to a common inhibitor.

A series of 8-arylated purine derivatives bearing either an aniline or an alkyl amide at position 6 were found to inhibit glycogen synthase kinase-3, with good selectivity over ten kinases. Molecular modeling studies indicated that the most active compounds (8a and 8e), adopt a planar conformation, close to the shape of AMPPNP in the crystal structure of GSK-3. These compounds are stabilized by hydrophobic contacts between the 8-aromatic group and the protein adenine pocket and by electrostatic contacts., (Copyright (c) 2010 Elsevier Masson SAS. All rights reserved.)

Centrin, an EF-hand calcium-binding protein, has been shown to be involved in the duplication of centrosomes, and Sfi1 (Suppressor of fermentation-induced loss of stress resistance protein 1) is one of its centrosomal targets. There are three isoforms of human centrin, but here we only considered centrin 2 (HsCen2). This protein has the ability to bind to any of the approximately 25 repeats of human Sfi1 (hSfi1) with more or less affinity. In this study, we mainly focused on the 17th repeat (R17-hSfi1-20), which presents the highest level of similarity with a well-studied 17-residue peptide (P17-XPC) from human xeroderma pigmentosum complementation group C protein, another centrin target for DNA repair. The only known structure of HsCen2 was resolved in complex with P17-XPC. The 20-residue peptide R17-hSfi1-20 exhibits the motif L8L4W1, which is the reverse of the XPC motif, W1L4L8. Consequently, the dipole of the helix formed by this motif has a reverse orientation. We wished to ascertain the impact of this reversal on the structure, dynamics and affinity of centrin. To address this question, we determined the structure of C-HsCen2 [the C-terminal domain of HsCen2 (T94-Y172)] in complex with R17-hSfi1-20 and monitored its dynamics by NMR, after having verified that the N-terminal domain of HsCen2 does not interact with the peptide. The structure shows that the binding mode is similar to that of P17-XPC. However, we observed a 2 -A translation of the R17-hSfi1-20 helix along its axis, inducing less anchorage in the protein and the disruption of a hydrogen bond between a tryptophan residue in the peptide and a well-conserved nearby glutamate in C-HsCen2. NMR dynamic studies of the complex strongly suggested the existence of an unusual calcium secondary binding mode in calcium-binding loop III, made possible by the uncommon residue composition of this loop. The secondary metal site is only populated at high calcium concentration and depends on the type of bound ligand.

The EF-hand calcium-binding proteins may exist either in an extended or a compact conformation. This conformation is sometimes correlated with the function of the calcium-binding protein. For those proteins whose structure and function are known, calcium sensors are usually extended and calcium buffers compact; hence, there is interest in predicting the form of the protein starting from its sequence. In the present study, we used two different procedures: one that already exists in the literature, the sosuidumbbell algorithm, mainly based on the charges of the two EF-hand domains, and the other comprising a novel procedure that is based on linker average hydrophilicity. The linker consists of the residues that connect the domains. The two procedures were tested on 17 known-structure calcium-binding proteins and then applied to 59 unknown-structure centrins. The sosuidumbbell algorithm yielded the correct conformations for only 15 of the known-structure proteins and predicted that all centrins should be in a closed form. The linker average hydrophilicity procedure discriminated well between all the extended and non-extended forms of the known-structure calcium-binding proteins, and its prediction concerning centrins reflected well their phylogenetic classification. The linker average hydrophilicity criterion is a simple and powerful means to discriminate between extended and non-extended forms of calcium-binding proteins. What is remarkable is that only a few residues that constitute the linker (between 2 and 20 in our tested sample of proteins) are responsible for the form of the calcium-binding protein, showing that this form is mainly governed by short-range interactions.

Protein kinase CK2 is a Ser/Thr kinase, with a constitutive activity, that is considered as a promising target for cancer therapy. The currently available CK2 inhibitors lack the potency and the pharmacological properties necessary to be suitable and successful in clinical settings. We report the development of new potent CK2 inhibitors from salicylaldehyde derivatives identified by automated screening of a proprietary small-molecule library. Docking simulations and analysis of the structure-activity relationship for the hits allowed to determine their binding modes on CK2, and to carry out the optimization of their structures. This strategy led to the discovery of potent CK2 inhibitors with novel structures, one of which was able to inhibit CK2 activity in living cells and promote tumor cell death. The essential features required for potent CK2 inhibitory activity of this class of compounds are discussed.

None of the already described CK2 inhibitors did fulfill the requirements for successful clinical settings. In order to find innovative CK2 inhibitors based on new scaffolds, we have performed a high-throughput screening of diverse chemical libraries. We report here the identification and characterization of several classes of new inhibitors. Whereas some share characteristics of previously known CK2 inhibitors, others are chemically unrelated and may represent new opportunities for the development of better CK2 inhibitors. By combining structure-activity relationships with a docking procedure, we were able to determine the binding mode of these inhibitors. Interestingly, beside the identification of several nanomolar ATP-competitive inhibitors, one class of chemical inhibitors displays a non-ATP competitive mode of inhibition, a feature that suggests that CK2 possess distinct druggable binding sites. For the most promising inhibitors, selectivity profiling was performed. We also provide evidence that some chemical compounds are inhibiting CK2 in living cells. Finally, the collected data allowed us to draw the rules about the chemical requirements for CK2 inhibition both in vitro and in a cellular context.

Using magnetic tweezers to investigate the mechanical response of single chromatin fibers, we show that fibers submitted to large positive torsion transiently trap positive turns at a rate of one turn per nucleosome. A comparison with the response of fibers of tetrasomes (the [H3-H4](2) tetramer bound with approximately 50 bp of DNA) obtained by depletion of H2A-H2B dimers suggests that the trapping reflects a nucleosome chiral transition to a metastable form built on the previously documented right-handed tetrasome. In view of its low energy, <8 kT, we propose that this transition is physiologically relevant and serves to break the docking of the dimers on the tetramer that in the absence of other factors exerts a strong block against elongation of transcription by the main RNA polymerase.

Previous laser flash photolysis investigations between 100 and 300 K have shown that the kinetics of CO rebinding with cytochrome P450(cam)(camphor) consist of up to four different processes revealing a complex internal dynamics after ligand dissociation. In the present work, molecular dynamics simulations were undertaken on the ternary complex P450(cam)(cam)(CO) to explore the CO migration pathways, monitor the internal cavities of the protein, and localize the CO docking sites. One trajectory of 1 nsec with the protein in a water box and 36 trajectories of 1 nsec in the vacuum were calculated. In each trajectory, the protein contained only one CO ligand on which no constraints were applied. The simulations were performed at 200, 300, and 320 K. The results indicate the presence of seven CO docking sites, mainly hydrophobic, located in the same moiety of the protein. Two of them coincide with xenon binding sites identified by crystallography. The protein matrix exhibits eight persistent internal cavities, four of which corresponding to the ligand docking sites. In addition, it was observed that water molecules entering the protein were mainly attracted into the polar pockets, far away from the CO docking sites. Finally, the identified CO migration pathways provide a consistent interpretation of the experimental rebinding kinetics.

DFT calculations are carried out on computational models of the active center of the alpha and beta subunits of hemoglobin in both its oxygenated (R) and deoxygenated (T) states. The computational models are defined by the full heme group, including all porphyrin substituents, and the four amino acids closer to it. The role of the protein environment is introduced by freezing the position of the alpha carbon atom of each of the four amino acids to the positions they have in the available PDB structures. Oxygen affinity is then evaluated by computing the energy difference between the optimized structures of the oxygenated and deoxygenated forms of each model. The results indicate a higher affinity of the alpha subunits over the beta ones. Analysis of the computed structures points out to the strength of the hydrogen bond between the distal histidine and the oxygen molecule as a key factor in discriminating the different systems.

Human centrin 2 is a component of the nucleotide excision repair system, as a subunit of the heterotrimer including xeroderma pigmentosum group C protein (XPC) and hHR23B. The C-terminal domain of centrin (C-HsCen2) binds strongly a peptide from the XPC protein (P1-XPC: N(847)-R(863)). Here, we characterize the solution Ca(2+)-dependent structural and molecular features of the C-HsCen2 in complex with P1-XPC, mainly using NMR spectroscopy and molecular modeling. The N-terminal half of the peptide, organized as an alpha helix is anchored into a deep hydrophobic cavity of the protein, because of three bulky hydrophobic residues in position 1-4-8 and electrostatic contacts with the centrin helix E. Investigation of the whole centrin interactions shows that the N-terminal domain of the protein is not involved in the complex formation and is structurally independent from the peptide-bound C-terminal domain. The complex may exist in three different binding conformations corresponding to zero, one, and two Ca(2+)-bound states, which may exchange with various rates and have distinct structural stability. The various features of the intermolecular interaction presented here constitute a centrin-specific mode for the target binding.

Molecular dynamics simulations of the unliganded T state of human hemoglobin showed the existence of a spontaneous, very wide cavity on the distal side of the alpha subunit. This cavity consists of three tunnels spreading from the vicinity of the iron atom (the ligand binding site) to the surface of the subunit, constituting possible passageways for the entrance of the ligand. A fourth passageway was characterized due to the trajectory of water molecules entering or leaving the heme pocket. Analogous passages were observed in the beta subunits. They all appear and disappear dynamically, although some parts of them are more persistent along the trajectories. The most persistent regions within these tunnels correspond to all the xenon docking sites of human cytoglobin and to some of those of sperm whale and horse heart myoglobins and group I truncated hemoglobins.

Iron (II) basket-handle porphyrins (BHP) are a series of encumbered heme models designed several years ago to mimic the ligand binding site of hemoproteins. Contrary to expectations, kinetic investigations have revealed that the k(on) rates for CO and/or O2 binding were only marginally affected by the assumed central steric hindrance of the iron atom. Thus, it was hypothesized that the internal dynamics of the molecule might be at the origin of the poor steric protection. To address this issue, measurements of nuclear magnetic resonance relaxation rates, fluorescence anisotropy experiments, and molecular dynamics simulations were undertaken. The size of BHP is small enough to allow the simulation in explicit chloroform with an almost complete sampling of the conformational space. The order parameters calculated from the MD trajectory compare well with the NMR experimental data and the predicted rotational correlation time corresponding to the Brownian motion of the molecule is in good agreement with the fluorescence measurements. Moreover, combining the results obtained using the three techniques allows the attribution of each internal NMR correlation time to a particular internal motion, revealing that even such medium-sized molecules are able to display quite complex internal dynamics. In particular, the handle phenyls that were assumed to sandwich the porphyrin have in fact a vanishing probability to be found in the proximity of the iron atom. They are therefore unable to reduce ligand accessibility significantly, which may explain the behavior of the k(on) rates.

In this work we show that ligand migration and active site conformational relaxation can occur independently of each other in hemoproteins. The complicated kinetics of carbon monoxide rebinding with cytochrome P450cam display up to five distinct processes between 77 K and 300 K. They were disentangled by using a combination of three approaches: 1), the competition of the ligand with xenon for the occupation of internal protein cavities; 2), the modulation of the amount of distal steric hindrance within the heme pocket by varying the nature of the substrate; and 3), molecular mechanics calculations to support the proposed heme-substrate relaxation mechanism and to seek internal cavities. In cytochrome P450cam, active site conformational relaxation results from the displacement of the substrate toward the heme center upon photodissociation of the ligand. It is responsible for the long, puzzling bimodal nature of the rebinding kinetics observed down to 77 K. The relaxation rate is strongly substrate-dependent. Ligand migration is slower and is observed only above 135 K. Migration and return rates are independent of the substrate.

The proteins Arno and Gea2 of the Sec7 family can promote GDP-GTP exchange on Arf1, a small GTP-binding protein, which coordinates coated vesicle formation for protein transport within the cell. Crystal structures of the essential Sec7 domain (Sec7d) of Gea2 in the free and Arf1-bound forms suggest that conformational dynamics of the Sec7d as well as those of the G-protein play a role in nucleotide exchange. Starting from a set of complementary crystal structures, we compared the collective movements of unbound Gea2 and Arno Sec7 domains, Arf1-GDP, and the Arf1-Gea2(Sec7d) nucleotide-free complex using normal modes analyses. In all unbound Sec7d analyses, significant low-energy movements were found to lead to closure of the hydrophobic groove towards the form seen in the Arf1-Gea2(Sec7d) complex, suggesting that groove closure is a general feature of the Sec7 family. Low-energy movements in Arf1-GDP implicate critical switch 1 and 2 residues which are coupled to modifications in the myristoylated N-terminal-helix binding site at the other end of the "interswitch" beta hairpin. It is suggested that Sec7d groove closure upon docking of the two molecules may permit extraction of switch 1 from Arf1-GDP and prepare the complex for movement of the interswitch, which is central to the membrane-linked exchange activity. Large-scale collective movements in the Arf1-Sec7d complex appear to participate in the insertion of the Sec7d Glu finger into the GDP binding site to promote actual nucleotide release.

FixJ is a two-domain response regulator involved in nitrogen fixation in Sinorhizobium meliloti. Recent X-ray characterization of both the native (unphosphorylated) and the active (phosphorylated) states of the protein identify conformational changes of the beta4-alpha4 loop and the conserved residue Phe101 as the key switches in activation. These structures also allowed investigation of the transition between conformations of this two-component regulatory receiver domain by molecular dynamics simulations. The path for the conformational change was studied with a distance constraint directing the system from one state to the other. The simulations provide evidence for a correlation between the conformation of the beta4-alpha4 loop and the orientation of the residue Phe101. A model presenting the sequence of events during the activation/deactivation process is discussed.