Comprehensive analysis of intramolecular G-quadruplex structures: furthering the understanding of their formalism

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.

ASC-G4, an algorithm to calculate advanced structural characteristics of G-quadruplexes

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.

Optimization of G-Quadruplex Ligands through a SAR Study Combining Parallel Synthesis and Screening of Cationic Bis(acylhydrazones)

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.

Deciphering the Mechanism of Inhibition of SERCA1a by Sarcolipin Using Molecular Simulations

SERCA1a is an ATPase calcium pump that transports Ca²⁺ from the cytoplasm to the sarco/endoplasmic reticulum lumen. Sarcolipin (SLN), a transmembrane peptide, regulates the activity of SERCA1a by decreasing its Ca²⁺ 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²⁺ binding sites, which would decrease the affinity for Ca²⁺. 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 ³⁵¹DKTG³⁵⁴ 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²⁺ low-affinity E2-like state in the presence of SLN and toward a Ca²⁺ 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.

Identification of New Potent Acetylcholinesterase Inhibitors Using Virtual Screening and in vitro Approaches

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.

Flagellar microtubule doublet assembly in vitro reveals a regulatory role of tubulin C-terminal tails

Microtubule doublets (MTDs), consisting of an incomplete B-microtubule at the surface of a complete A-microtubule, provide a structural scaffold mediating intraflagellar transport and ciliary beating. Despite the fundamental role of MTDs, the molecular mechanism governing their formation is unknown. We used a cell-free assay to demonstrate a crucial inhibitory role of the carboxyl-terminal (C-terminal) tail of tubulin in MTD assembly. Removal of the C-terminal tail of an assembled A-microtubule allowed for the nucleation of a B-microtubule on its surface. C-terminal tails of only one A-microtubule protofilament inhibited this side-to-surface tubulin interaction, which would be overcome in vivo with binding protein partners. The dynamics of B-microtubule nucleation and its distinctive isotropic elongation was elucidated by using live imaging. Thus, inherent interaction properties of tubulin provide a structural basis driving flagellar MTD assembly.

Efficient conformational sampling and weak scoring in docking programs? Strategy of the wisdom of crowds

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 Å.

Structure-based mutational analysis of ICAT residues mediating negative regulation of β-catenin co-transcriptional activity

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.

Benchmark of four popular virtual screening programs: construction of the active/decoy dataset remains a major determinant of measured performance

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 Abstract

vSDC: a method to improve early recognition in virtual screening when limited experimental resources are available

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 abstract

An atomistic view of microtubule stabilization by GTP

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.

Expanding the chemical diversity of CK2 inhibitors

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.

A DFT study on the relative affinity for oxygen of the alpha and beta subunits of hemoglobin

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.

Nucleosome chiral transition under positive torsional stress in single chromatin fibers

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.

Molecular dynamics assignment of NMR correlation times to specific motions in a "basket-handle porphyrin" heme

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.

CO migration pathways in cytochrome P450cam studied by molecular dynamics simulations

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.

Ligand migration and escape pathways in haem proteins

Biophysical techniques developed during the last three decades have provided an increasingly detailed description of the internal processes associated with ligand capture and release by haem proteins. Myoglobin has long been the paradigm for these studies. More recently, cytochrome P450cam (the camphor-metabolizing cytochrome P450 from Pseudomonas putida) has also received considerable interest. In spite of sharing the same prosthetic group, the Fe(II)-haem, these proteins are structurally unrelated and they perform different functions. Recent works show that both proteins exhibit a common feature: a series of permanent or fluctuating, mostly hydrophobic, cavities of the protein matrix are providing transient docking sites as well as migration, escape and possibly entry pathways for the ligand. Remarkably, these systems of cavities connect the distal and the proximal regions of the haem, a disposition that may contribute to ligand capture enhancement.

Flexibility and Plasticity of Human Centrin 2 Binding to the Xeroderma Pigmentosum Group C Protein (XPC) from Nuclear Excision Repair,

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.

Internal cavities and ligand passageways in human hemoglobin characterized by molecular dynamics simulations

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.

Disentangling ligand migration and heme pocket relaxation in cytochrome P450cam

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.

Integrating three views of Arf1 activation dynamics

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.

Molecular dynamics of the FixJ receiver domain: movement of the beta4-alpha4 loop correlates with the in and out flip of Phe101

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.

New insights into the allosteric mechanism of human hemoglobin from molecular dynamics simulations

It is still difficult to obtain a precise structural description of the transition between the deoxy T-state and oxy R-state conformations of human hemoglobin, despite a large number of experimental studies. We used molecular dynamics with the Path Exploration with Distance Constraints (PEDC) method to provide new insights into the allosteric mechanism at the atomic level, by simulating the T-to-R transition. The T-state molecule in the absence of ligands was seen to have a natural propensity for dimer rotation, which nevertheless would be hampered by steric hindrance in the "joint" region. The binding of a ligand to the alpha subunit would prevent such hindrance due to the coupling between this region and the alpha proximal histidine, and thus facilitate completion of the dimer rotation. Near the end of this quaternary transition, the "switch" region adopts the R conformation, resulting in a shift of the beta proximal histidine. This leads to a sliding of the beta-heme, the effect of which is to open the beta-heme's distal side, increasing the accessibility of the Fe atom and thereby the affinity of the protein. Our simulations are globally consistent with the Perutz strereochemical mechanism.

Motions in hemoglobin studied by normal mode analysis and energy minimization: evidence for the existence of tertiary T-like, quaternary R-like intermediate structures

The normal mode analysis of human hemoglobin showed the presence in the deoxy T-state of one main preferential direction that brings the structure close to the R-state, with a low-energy variation, while in the oxy R-state there are several modes that point towards the T-state, but with higher energy variations and less contribution to the transition. The displacement along a combination of normal modes, followed by energy minimization, starting from the R-state, did not allow one to obtain a structure significantly different from that of R, showing that the fully oxygenated hemoglobin is trapped in a deep and narrow potential energy minimum. On the contrary, starting from the deoxy T-state, the displacement along a combination of normal modes, followed by energy minimization, yielded an intermediate structure, that we designate Tmin(d1), which is closer to R; the normal modes of Tmin(d1) indicated that the potential energy minimum in the vicinity of this structure is as narrow as that of R but less deep. The procedure of displacement along the modes, followed by energy minimization, was applied to Tmin(d1), yielding Tmin(d2); then the procedure was repeated, yielding the intermediate structures Tmin(d3) and Tmin(d4). The structures Tmin(d2), Tmin(d3) and Tmin(d4) are not significantly different from each other, indicating that they are trapped in a narrow, deep energy minimum. This procedure revealed the existence of at least two intermediate sets of structures between T and R: the first one, Tmin(d1), is different from the T and R structures, while the second set, Tmin(d2), Tmin(d3) and Tmin(d4), is quaternary R-like and tertiary T-like, where the contacts at the interfaces alpha1 beta1 and alpha1 beta2 are R-like, and the alpha and beta heme environments are still T-like.

Analysis of the low frequency normal modes of the T-state of aspartate transcarbamylase

Aspartate transcarbamylase (ATCase) is an important control enzyme in the pyrimidine biosynthetic pathway in Escherichia coli. It is a classic example of an allosteric protein and has been extensively studied biochemically, kinetically and structurally. As yet, however, a detailed model for the cooperative transition between the tensed (T) and relaxed (R) forms of the protein does not exist. In this work we have calculated the low frequency normal modes of the CTP-ligated T-state of ATCase with the aim of identifying some of the motions that could be important in initiating the transition. The calculated modes, of frequencies lower than 5 per cm, produce root-mean-square coordinate deviations for the atoms which are a substantial fraction of those derived from the crystallographic B-factors. Some of the modes result in displacements which change the quaternary structure of the protein (in particular the elongation of the protein and the relative rotation of the subunits) in such a way that the R-state structure is approached. The implication of these mode motions for the overall T-->R transition process is discussed.

Computation of low-frequency normal modes in macromolecules: improvements to the method of diagonalization in a mixed basis and application to hemoglobin

The method of diagonalization in a mixed basis (DIMB) that was published previously (Mouawad.), L. and Perahia D., Biopolymers 33, 599, 1993), allows the computation of the low-frequency vibrational modes for large macromolecules. Improvements to this method are presented here, namely the single and double truncation window techniques. The best convergence rate is obtained with the double truncation windows, which couple most efficiently the parts of the macromolecule which are far in sequence but close in space. Both methods were applied to the T-state of hemoglobin, to compare their efficiency. The resulting modes are analyzed in order to study the pathways between T- and R-states of this protein. They show that the quaternary conformational are mainly due to one mode at 2 cm-1.

Proton NMR studies of the structural and dynamical effect of chemical modification of a single aromatic side-chain in a snake cardiotoxin. Relation to the structure of the putative binding site and the cytolytic activity of the toxin

This paper presents the comparative comprehensive analysis of NMR structural parameters (NOEs, scalar coupling, chemical shifts) of toxin gamma, a cardiotoxin isolated from the venom of Naja nigricollis, and three chemical derivatives, i.e. the 2-nitrophenylsulphonyl (NPS)-Trp11, 3-nitro-Tyr22 and 3-nitro-Tyr51 derivatives. In previous work, the chemical modifications of single side chains have suggested that these aromatic residues, in association with several lysine residues, contributed to the cytotoxicity of toxin gamma. Analysis of these results based on the refined solution structure of the toxin has resulted in the proposal of a conserved phospholipid binding site through which cardiotoxins are likely to interact with the membrane of target cells. The present work shows that modifications of either the tryptophan residue or the tyrosine residues, which are within or near the proposed binding site, have no influence on the three-dimensional structure of the protein. On the other hand, the proton exchange study of the backbone amides indicates that the structural core of the protein is destabilized in the three derivatives. This corresponds to a decrease of the overall stability of the protein as indicated by the comparative solvent denaturation study of the unmodified toxin gamma and the Trp11 derivative. More specifically, the dynamics of the three-stranded beta sheet, a part of the structural core, are highly perturbed by the chemical modifications. This sheet was previously proposed as a part of the phospholipid binding site of cardiotoxins. The dynamical perturbation of this site appears to be correlated with the decrease in toxicity of the chemical derivatives.

Do cardiotoxins possess a functional site? Structural and chemical modification studies reveal the functional site of the cardiotoxin from Naja nigricollis

Examination of the literature has revealed that regarding the amino acid sequences, cardiotoxins constitute a family of homogeneous compounds. In contrast, cardiotoxins appear heterogeneous as far as their biological and spectroscopic properties are concerned. As a result, comparison between these molecules with a view to establishing structure-activity correlations is complicated. We have therefore reviewed recent works aiming at identifying the functional site of a defined cardiotoxin, ie toxin gamma from the venom of the spitting cobra Naja nigricollis. The biological and structural properties of toxin gamma are first described. In particular, a model depicting the 3-dimensional structure of the toxin studied by NMR spectroscopy is proposed. The toxin polypeptide chain is folded into 3 adjacent loops rich in beta-sheet structure connected to a small globular core containing the 4 disulfide bonds. A number of derivatives chemically modified at a single aromatic or amino group have been prepared. The structure of each derivative was probed by emission fluorescence, circular dichroism and NMR spectroscopy. Also tested was the ability of the derivatives to kill mice, depolarize excitable cell membranes and lyse epithelial cells. Modification of some residues in the first loop, in particular Lys-12 and at the base of the second loop substantially affected biological properties, with no sign of concomitant structural modifications other than local changes. Modifications in other regions much less affected the biological properties of the toxin. A plausible functional site for toxin gamma involving loop I and the base of loop II is presented. It is stressed that the functional site of other cardiotoxins may be different.

Flexibility and folding of phosphoglycerate kinase

Flexibility and folding of phosphoglycerate kinase, a two-domain monomeric enzyme, have been studied using a wide variety of methods including theoretical approaches. Mutants of yeast phosphoglycerate kinase have been prepared in order to introduce cysteinyl residues as local probes throughout the molecule without perturbating significantly the structural or the functional properties of the enzyme. The apparent reactivity of a unique cysteine in each mutant has been used to study the flexibility of PGK. The regions of larger mobility have been found around residue 183 on segment beta F in the N-domain and residue 376 on helix XII in the C-domain. These regions are also parts of the molecule which unfold first. Ligand binding induces conformational motions in the molecule, especially in the regions located in the cleft. Moreover, the results obtained by introducing a fluorescent probe covalently linked to a cysteine are in agreement with the helix scissor motion of helices 7 and 14 assumed by Blake to direct the hinge bending motion of the domains during the catalytic cycle. The folding process of both horse muscle and yeast phosphoglycerate kinases involves intermediates. These intermediates are more stable in the horse muscle than in the yeast enzyme. In both enzymes, domains behave as structural modules capable of folding and stabilizing independently, but in the horse muscle enzyme the C-domain is more stable and refolds prior to the N-domain, contrary to that which has been observed in the yeast enzyme. A direct demonstration of the independence of domains in yeast phosphoglycerate kinase has been provided following the obtention of separated domains by site-directed mutagenesis. These domains have a native-like structure and refold spontaneously after denaturation by guanidine hydrochloride.

The effects of ligands on the conformation of phosphoglycerate kinase: fluorescence anisotropy decay and theoretical interpretation

Horse muscle phosphoglycerate kinase (PGK) is a monomer folded into two widely distant domains. In the glycolytic pathway, this enzyme catalyzes the first reaction that produces ATP. It was suggested, by analogy with yeast hexokinase, that a hinge-bending motion may be induced by the binding of specific substrates to the protein. To analyze such a motion, or any structural changes induced by ligand binding, fluorescence anisotropy decay of tryptophan residues in free and liganded PGK was studied. At 293 K, for the free protein and the binary complex with 3-phosphoglycerate, a single correlation time of 26 ns was observed, corresponding to the rotation of the overall protein, whereas upon addition of MgADP, this correlation time decreased to 10 ns. Such a decrease cannot be merely due to a change of the protein's shape and volume. To explain this, it was suggested that the fluorescence anisotropy decay of the PGK-MgADP complex corresponded to the rotation of the only buried tryptophan (Trp 335). The rotational paths of this tryptophan, in the presence and absence of the nucleotide, were established by potential energy minimization calculations. The results indicated that MgADP induces a displacement of helix alpha-13 that decreases the rotational energy barrier of Trp 335 from 16 kcal/mol in the free protein to 8 kcal/mol in the complex.

Introduction of internal cysteines as conformational probes in yeast phosphoglycerate kinase

Several mutants of yeast phosphoglycerate kinase, each containing only one internal cysteine residue, were constructed from a single mutant devoid of cysteine. These cysteines were introduced as local conformational probes in selected buried positions. The enzyme activity, conformational characteristics and stability indicated that the mutations introduced only small perturbations in the molecule. The folding-unfolding process mediated by guanidine hydrochloride under equilibrium conditions was studied by following the variations in ellipticity and the reactivity of the cysteine residue towards 5,5'-dithiobis(nitrobenzoate). The process was found to be reversible except for mutant C97A, V49C, suggesting that this region located in helix I might be crucial in determining an intermediate on the folding pathway. The transitions obtained by the two signals did not coincide, indicating that the local structures, in several parts inside the molecule, are more sensitive to the denaturant than the overall conformation.

Study of the fast-reacting cysteines in phosphoglycerate kinase using chemical modification and site-directed mutagenesis

Horse muscle phosphoglycerate kinase, like other mammalian phosphoglycerate kinases, contains seven cysteine residues of which two react rapidly with 5,5'-dithio-bis(2-nitrobenzoate) (Nbs2) following second-order kinetics (k = 640 M-1.s-1). Selective cyanylation of the fast-reacting cysteines, followed by chemical cleavage and subsequent sodium dodecyl sulfate/polyacrylamide gel electrophoresis analysis of the resulting polypeptides, suggested that these cysteines are at positions 378 and 379. Cysteine residues were introduced into yeast phosphoglycerate kinase by site-directed mutagenesis. Mutant enzymes, each containing only one cysteine residue at position 364, 376, or 377, were constructed from a mutant devoid of cysteine (Cys97----Ala). In the last two mutants, the cysteines were at positions corresponding to Cys378 and Cys379, respectively, in the horse muscle enzyme. The chemical reactivity of the cysteine groups in these latter two yeast mutant enzymes was similar to that of the fast-reacting cysteines in the horse muscle enzyme. Furthermore, they were similarly modified upon substrate binding. All these data demonstrate unambiguously that the fast-reacting cysteines in the horse muscle enzyme are Cys378 and Cys379.