MM/PBSA analysis of molecular dynamics simulations of bovine beta-lactoglobulin: free energy gradients in conformational transitions?

A 200 nanoseconds all-atom simulation of the pH-dependent EF loop transition in bovine β-lactoglobulin. The role of the orientation of the E89 side chain

In silico molecular dynamics (MD) using crystallographic and NMR data was used to simulate the effects of the protonation state of E89 on the pH-dependent conformational rearrangement of the EF loop, also known as the Tanford transition, in a series of apo-β-lactoglobulin (BLG) structures. Compared to existing studies these simulations were carried out over a much longer time scale (200 ns where the stability of the transition can be evaluated) and used an explicit water model. We considered eight different entries from the Brookhaven Protein Data Bank (PDB) separated into two groups. We observed that fixing the protonation state of E89 prompts the transition of the EF loop only when its side chain is oriented under the loop and into the entrance of the interior cavity. The motion of the EF loop occurs mostly as a step-function and its timing varies greatly from ∼ 20 ns to ∼170 ns from the beginning of the simulation. Once the transition is completed, the protein appears to reach a stable conformation as in a true two-state transition. We also observed novel findings. When the transition occurs, the hydrogen bond between E89 and S116 is replaced with a salt bridge with Lys residues in the βC-CD loop-βD motif. This electrostatic interaction causes the distortion of this motif as well as the protrusion of the GH loop into the aperture of the cavity with the result of limiting the increase of its contour area despite the rotation of the EF loop.Communicated by Ramaswamy H. Sarma.

Exploring the molecular basis of UG-rich RNA recognition by the human splicing factor TDP-43 using molecular dynamics simulation and free energy calculation

Transactivation response element RNA/DNA-binding protein 43 (TDP-43) is involved in the regulation of alternative splicing of human neurodegenerative disease-related genes through binding to long UG-rich RNA sequences. Mutations in TDP-43, most in the homeodomain, cause neurological disorders such as amyotrophic lateral sclerosis and fronto temporal lobar degeneration. Several mutants destabilize the structure and disrupt RNA-binding activity. The biological functions of these mutants have been characterized, but the structural basis behind the loss of RNA-binding activity is unclear. Focused on the specific TDP-43-ssRNA complex (PDB code 4BS2), we applied molecular dynamics simulations and the molecular mechanics Poisson-Boltzmann surface area free energy calculation to characterize and explore the structural and dynamic effects between ssRNA and TDP-43. The energetic analysis indicated that the intermolecular van der Waals interaction and nonpolar solvation energy play an important role in the binding process of TDP-43 and ssRNA. Compared with the wild-type TDP-43, the reduction of the polar or non-polar interaction between all the mutants F149A, D105A/S254A, R171A/D174A, F147L/F149L/F229L/F231L and ssRNA is the main reason for the reduction of its binding free energy. Decomposing energies suggested that the extensive interactions between TDP-43 and the nitrogenous bases of ssRNA are responsible for the specific ssRNA recognition by TDP-43. These results elucidated the TDP-43-ssRNA interaction comprehensively and further extended our understanding of the previous experimental data. The uncovering of TDP-43-ssRNA recognition mechanism will provide us useful insights and new chances for the development of anti-neurodegenerative drugs.

Structural mechanism of the Tanford transition of bovine β-lactoglobulin through microsecond molecular dynamics simulations

X-ray analysis has provided structural data about a pH-driven conformational change in β-lactoglobulin (BLG) known as the Tanford transition, which occurs at around pH 7 and involves the EF loop, which acts as a lid closing the internal cavity of the protein below pH 7 and opening it above pH 7. NMR studies using wild-type BLG have encountered problems trying to explain the Tanford transition, however, they have provided important insight using a dimeric BLG mutant, revealing that the opening and closure of the EF loop consists of two types of motions in the microsecond and milliseconds timescales. This provides valuable information indicating that the dimeric state is a good model to study the Tanford transition, although the understanding of this structural change is still lacking at the atomic level. We performed microsecond molecular dynamics (MD) simulations starting from different conformations of BLG in the monomeric and dimeric state, with protonated and deprotonated E89, in order to explore the Tanford transition. Our results provide structural information for the transition from the closed to the open conformation in BLG and show it occurs in the dimeric state in the microsecond timescale, in line with the fast motion observed through NMR experiments. In addition, MD simulations coupled to MMGBSA approach indicated that the most populated conformer of BLG in the open state is able to bind ligands with similar affinity to that of BLG at neutral pH obtained through crystallographic experiments.Communicated by Ramaswamy H. Sarma.

Effective prediction model and determination of binding residues influential for inhibitors targeting HIV-1 integrase-LEDGF/p75 interface by employing solvent accessible surface area energy as key determinant

Development of a highly accurate prediction model for protein-ligand inhibition has been a major challenge in drug discovery. Herein, we describe a novel predictive model for the inhibition of HIV-1 integrase (IN)-LEDGF/p75 protein-protein interaction. The model was constructed using energy parameters approximated from molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations. Chemometric analysis using partial least squares (PLS) regression revealed that solvent accessible surface area energy (ΔG_SASA) is the major determinant parameter contributing greatly to the prediction accuracy. PLS prediction model on the ΔG_SASA values collected from 41 complexes yielded a strong correlation between the predicted and the actual inhibitory activities (R² = 0.9666, RMSEC of pIC₅₀ values = 0.0890). Additionally, for the test set of 14 complexes, the model performed satisfactorily with very low pIC₅₀ errors (Q² = 0.5168, RMSEP = 0.3325). A strong correlation between the buried surface areas on the IN protein, when bound with IN-LEDGF/p75 inhibitors, and the respective ΔG_SASA values was also obtained. Furthermore, the current method could identify 'hot spots'of amino acid residues highly influential to the inhibitory activity prediction. This could present fruitful implications in binding site determination and future inhibitor developments targeting protein-protein interactions.Communicated by Ramaswamy H. Sarma.

β-lactoglobulin's conformational requirements for ligand binding at the calyx and the dimer interphase: a flexible docking study

β-lactoglobulin (BLG) is an abundant milk protein relevant for industry and biotechnology, due significantly to its ability to bind a wide range of polar and apolar ligands. While hydrophobic ligand sites are known, sites for hydrophilic ligands such as the prevalent milk sugar, lactose, remain undetermined. Through the use of molecular docking we first, analyzed the known fatty acid binding sites in order to dissect their atomistic determinants and second, predicted the interaction sites for lactose with monomeric and dimeric BLG. We validated our approach against BLG structures co-crystallized with ligands and report a computational setup with a reduced number of flexible residues that is able to reproduce experimental results with high precision. Blind dockings with and without flexible side chains on BLG showed that: i) 13 experimentally-determined ligands fit the calyx requiring minimal movement of up to 7 residues out of the 23 that constitute this binding site. ii) Lactose does not bind the calyx despite conformational flexibility, but binds the dimer interface and an alternate Site C. iii) Results point to a probable lactolation site in the BLG dimer interface, at K141, consistent with previous biochemical findings. In contrast, no accessible lysines are found near Site C. iv) lactose forms hydrogen bonds with residues from both monomers stabilizing the dimer through a claw-like structure. Overall, these results improve our understanding of BLG's binding sites, importantly narrowing down the calyx residues that control ligand binding. Moreover, our results emphasize the importance of the dimer interface as an insufficiently explored, biologically relevant binding site of particular importance for hydrophilic ligands. Furthermore our analyses suggest that BLG is a robust scaffold for multiple ligand-binding, suitable for protein design, and advance our molecular understanding of its ligand sites to a point that allows manipulation to control binding.

Detailed computational study of the active site of the hepatitis C viral RNA polymerase to aid novel drug design

The hepatitis C virus (HCV) RNA polymerase, NS5B, is a leading target for novel and selective HCV drug design. The enzyme has been the subject of intensive drug discovery aimed at developing direct acting antiviral (DAA) agents that inhibit its activity and hence prevent the virus from replicating its genome. In this study, we focus on one class of NS5B inhibitors, namely nucleos(t)ide mimetics. Forty-one distinct nucleotide structures have been modeled within the active site of NS5B for the six major HCV genotypes. Our comprehensive modeling protocol employed 287 different molecular dynamics simulations combined with the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) methodology to rank and analyze these structures for all genotypes. The binding interactions of the individual compounds have been investigated and reduced to the atomic level. The present study significantly refines our understanding of the mode of action of NS5B-nucleotide-inhibitors, identifies the key structural elements necessary for their activity, and implements the tools for ranking the potential of additional much needed novel inhibitors of NS5B.

Wards in the keyway: amino acids with anomalous pK(a)s in calycins

As a follow-up to our recent analysis of the electrostatics of bovine β-lactoglobulin (Eberini et al. in Amino Acids 42:2019-2030, 2011), we investigated whether the occurrence in the native structure of calycins-the superfamily to which β-lactoglobulin belongs-of amino acids with anomalous pK (a)s is an infrequent or, on the contrary, a common occurrence, and whether or not a general pattern may be recognized. To this aim, we randomly selected four calycins we had either purified from natural sources or prepared with recombinant DNA technologies during our previous and current structural and functional studies on this family. Their pIs vary over several pH units and their known functions are as diverse as carriers, enzymes, immunomodulators and/or extracellular chaperones. In our survey, we used both in silico prediction methods and in vitro procedures, such as isoelectric focusing, electrophoretic titration curves and spectroscopic techniques. By comparing the results under native conditions (no exposure of the proteins to chaotropic agents) to those after protein unfolding (in the presence of 8 M urea), a shift is observed in the pK (a) of at least one amino acid per protein, which results in a measurable change in pI. Three types of amino acids are involved: Cys, Glu, and His, their position varies along the calycin sequence. Although no common mechanism may thus be recognized, we hypothesize that the 'normalization' of anomalous pK (a)s may be the phenomenon that accompanies, and favors, structural rearrangements such as those involved in ligand binding by these proteins. An interesting, if anecdotal, validation to this view comes from the behavior of human retinol binding protein, for which the pI of the folded and liganded protein is intermediate between those of the folded and unliganded and of the unfolded protein forms. Likewise, both solid (from crystallography) and solution state (from CD spectroscopy) data confirm that the protein undergoes structural rearrangement upon retinol binding.

Effects of temperature on the p53-DNA binding interactions and their dynamical behavior: comparing the wild type to the R248Q mutant

Background

The protein p53 plays an active role in the regulation of cell cycle. In about half of human cancers, the protein is inactivated by mutations located primarily in its DNA-binding domain. Interestingly, a number of these mutations possess temperature-induced DNA-binding characteristics. A striking example is the mutation of Arg248 into glutamine or tryptophan. These mutants are defective for binding to DNA at 310 K although they have been shown to bind specifically to several p53 response elements at sub-physiological temperatures (298–306 K).

Methodology/Principal Findings

This important experimental finding motivated us to examine the effects of temperature on the structure and configuration of R248Q mutant and compare it to the wild type protein. Our aim is to determine how and where structural changes of mutant variants take place due to temperature changes. To answer these questions, we compared the mutant to the wild-type proteins from two different aspects. First, we investigated the systems at the atomistic level through their DNA-binding affinity, hydrogen bond networks and spatial distribution of water molecules. Next, we assessed changes in their long-lived conformational motions at the coarse-grained level through the collective dynamics of their side-chain and backbone atoms separately.

Conclusions

The experimentally observed effect of temperature on the DNA-binding properties of p53 is reproduced. Analysis of atomistic and coarse-grained data reveal that changes in binding are determined by a few key residues and provide a rationale for the mutant-loss of binding at physiological temperatures. The findings can potentially enable a rescue strategy for the mutant structure.

How to obtain statistically converged MM/GBSA results

The molecular mechanics/generalized Born surface area (MM/GBSA) method has been investigated with the aim of achieving a statistical precision of 1 kJ/mol for the results. We studied the binding of seven biotin analogues to avidin, taking advantage of the fact that the protein is a tetramer with four independent binding sites, which should give the same estimated binding affinities. We show that it is not enough to use a single long simulation (10 ns), because the standard error of such a calculation underestimates the difference between the four binding sites. Instead, it is better to run several independent simulations and average the results. With such an approach, we obtain the same results for the four binding sites, and any desired precision can be obtained by running a proper number of simulations. We discuss how the simulations should be performed to optimize the use of computer time. The correlation time between the MM/GBSA energies is approximately 5 ps and an equilibration time of 100 ps is needed. For MM/GBSA, we recommend a sampling time of 20-200 ps for each separate simulation, depending on the protein. With 200 ps production time, 5-50 separate simulations are required to reach a statistical precision of 1 kJ/mol (800-8000 energy calculations or 1.5-15 ns total simulation time per ligand) for the seven avidin ligands. This is an order of magnitude more than what is normally used, but such a number of simulations is needed to obtain statistically valid results for the MM/GBSA method.

Excited protein states of human tear lipocalin for low- and high-affinity ligand binding revealed by functional AB loop motion

Human tear lipocalin (TL), a prominent member of lipocalin family, exhibits functional and structural promiscuity. The plasticity of loop regions modulates entry to the ligand pocket at the "open" end of the eight-stranded beta-barrel. Site-directed multi-distance measurements using fluorescence resonance energy transfer between functional loops register two excited protein states for low- and high-affinity ligand binding. At low pH, the longest loop AB adopts the conformation of the low-affinity excited protein state that matches the crystal structure of holo-TL at pH 8. A "crankshaft" like movement is detected for the loop AB in a low pH transition. At pH 7.3 the holo-protein assumes a high-affinity excited protein state, in which the loop AB is more compact (RMS=3.1A). In the apo-holo transition, the reporter Trp 28 moves about 4.5A that reflects a decrease in distance between Glu27 and Lys108. This interaction fixes the loop AB conformation for the high-affinity mode. No such movement is detected at low pH, where Glu27 is protonated. Data strongly indicate that the protonation state of Glu27 modulates the conformation of the loop AB for high- and low-affinity binding.

pH-Dependent conformational changes in tear lipocalin by site-directed tryptophan fluorescence

Tear lipocalin (TL), a major protein of human tears, binds a broad array of endogenous ligands. pH-dependent ligand binding in TL may have functional implications in tears. Previously, conformational selections of the AB and GH loops have been implicated in ligand binding by site-directed tryptophan fluorescence (SDTF). In this study, SDTF was applied to the AB and GH loops to investigate pH-driven conformational changes relevant to ligand binding. Both loops demonstrate significant but distinct conformational rearrangements over a wide pH range. In the low-pH transition, from 7.3 to 3.0, residues of the GH loop exhibit decreased solvent accessibilities. In acrylamide quenching experiments, the average quenching rate constant (k(q), accessibility parameter) of the residues in the GH loop is decreased approximately 38%, from 2.1 x 10(9) to 1.3 x 10(9) M(-1) s(-1). However, despite the significant changes in accessibilities for some residues in the AB loop, the average accessibility per residue remained unchanged (average k(q) = 1.2 M(-1) s(-1)). Accordingly, the low-pH transition induces conformational changes that reshuffle the accessibility profiles of the residues in the AB loop. A significant difference in the titration curves between the holo and apo forms of the W28 mutant suggests that the protonation states of the residues around position 28 modulate conformational switches of the AB loop relevant to ligand binding.

Structure and folding of disulfide-rich miniproteins: insights from molecular dynamics simulations and MM-PBSA free energy calculations

The fold of small disulfide-rich proteins largely relies on two or more disulfide bridges that are main components of the hydrophobic core. Because of the small size of these proteins and their high cystine content, the cysteine connectivity has been difficult to ascertain in some cases, leading to uncertainties and debates in the literature. Here, we use molecular dynamics simulations and MM-PBSA free energy calculations to compare similar folds with different disulfide pairings in two disulfide-rich miniprotein families, namely the knottins and the short-chain scorpion toxins, for which the connectivity has been discussed. We first show that the MM-PBSA approach is able to discriminate the correct knotted topology of knottins from the laddered one. Interestingly, a comparison of the free energy components for kalata B1 and MCoTI-II suggests that cyclotides and squash inhibitors, although sharing the same scaffold, are stabilized through different interactions. Application to short-chain scorpion toxins suggests that the conventional cysteine pairing found in many homologous toxins is significantly more stable than the unconventional pairing reported for maurotoxin and for spinoxin. This would mean that native maurotoxin and spinoxin are not at the lowest free energy minimum and might result from kinetically rather than thermodynamically driven oxidative folding processes. For both knottins and toxins, the correct or conventional disulfide connectivities provide lower flexibilities and smaller deviations from the initial conformations. Overall, our work suggests that molecular dynamics simulations and the MM-PBSA approach to estimate free energies are useful tools to analyze and compare disulfide bridge connectivities in miniproteins.

Structure and dynamics of Cu(I) binding in copper chaperones Atox1 and CopZ: a computer simulation study

Copper chaperones deliver reduced copper (i.e., Cu(I)) to metal-binding domains of P-type ATPases in the cytoplasm of a range of organisms. Both chaperones and target domains have a ferredoxin-like fold and metal-binding motifs involving two Cys residues. Here, we investigated the Cu-binding geometry and structural dynamics of two homologous Cu(I) chaperones, Homo sapiens Atox1 and Bacillus subtilis CopZ, using a combination of quantum mechanical-molecular mechanics (QM-MM) and classical molecular dynamics (MD) methods. Our QM-MM optimized geometries for the holo- proteins suggested that Cu(I) in Atox1 favors a linear Cys(S)-Cu-Cys(S) arrangement but that this angle is close to 150 degrees in CopZ. Classical MD simulations suggest that both Atox1 and CopZ apo- forms have an increased conformational flexibility as compared to the respective holo- forms. This difference is most pronounced in CopZ and correlates with a lower in vitro thermal stability. Both average fluctuation (i.e., rmsd) and radius of gyration data demonstrate that the effects of Cu(I) coordination extend throughout the proteins. Distinct deviations between the two homologues were found in protein-solvent interactions, entropy of Cu(I) binding, and apo-protein Cys-Cys distance distributions. Our in silico results provide new insights into copper chaperone behavior with direct implications for copper transport mechanisms in vivo.

Protein-Ligand Binding Affinity Predictions by Implicit Solvent Simulations: A Tool for Lead Optimization?

Continuum electrostatics is combined with rigorous free-energy calculations in an effort to deliver a reliable and efficient method for in silico lead optimization. The methodology is tested by calculation of the relative binding free energies of a set of inhibitors of neuraminidase, cyclooxygenase2, and cyclin-dependent kinase 2. The calculated free energies are compared to the results obtained with explicit solvent simulations and empirical scoring functions. For cyclooxygenase2, deficiencies in the continuum electrostatics theory are identified and corrected with a modified simulation protocol. For neuraminidase, it is shown that a continuum representation of the solvent leads to markedly different protein-ligand interactions compared to the explicit solvent simulations, and a reconciliation of the two protocols is problematic. Cyclin-dependent kinase 2 proves more challenging, and none of the methods employed in this study yield high quality predictions. Despite the differences observed, for these systems, the use of an implicit solvent framework to predict the ranking of congeneric inhibitors to a protein is shown to be faster, as accurate or more accurate than the explicit solvent protocol, and superior to empirical scoring schemes.

Dynamics and Mechanism of the Tanford Transition of Bovine β-Lactoglobulin Studied using Heteronuclear NMR Spectroscopy

The Tanford transition is a conformational change of bovine beta-lactoglobulin (betaLG) occurring at around pH 7, identified originally on the basis of optical rotatory dispersion and the accessibility of a thiol group. X-ray analysis has suggested that a conformational change to the EF-loop is responsible for the Tanford transition, with the loop closing the hydrophobic cavity of the beta-barrel of the betaLG molecule below pH 7 and flipping to open the cavity above pH 7. To clarify the dynamics of this conformational change, NMR measurements were made at neutral pH. Since severe signal broadening due to monomer-dimer equilibrium prevented NMR measurements of wild-type betaLG at neutral pH, we searched for optimal sample conditions, finding that a disulfide bond-linked dimer of the mutant A34C gives an HSQC spectrum without signal broadening. The HSQC and CD spectra indicated that in overall conformation A34C is similar to wild-type betaLG, suggesting that the A34C dimer is a good model with which to study the structure and dynamics of the wild-type at neutral pH. The pH-dependent HSQC signal changes and Lipari-Szabo type relaxation analyses of the A34C dimer revealed that the conformational change to the EF-loop occurs above pH 7. We observed two types of motions in the EF-loop region; relatively fast (micro- to milliseconds) and slow (milliseconds or slower) conformational exchanges of the residues located in the hinge and top of the EF-loop regions, respectively. Furthermore, the GH-loop adjacent to the EF-loop exhibited conformational change at a pH slightly lower than that at which the EF-loop motions occurred. From these observations, we propose a three-step mechanism of conformational change in the EF-loop leading to the Tanford transition, in which the GH-loop conformational change, the cleavage of the hydrogen bonds at the hinge, and the flip of the EF-loop occur sequentially.

NMR dynamic studies suggest that allosteric activation regulates ligand binding in chicken liver bile acid-binding protein

Apo chicken liver bile acid-binding protein has been structurally characterized by NMR. The dynamic behavior of the protein in its apo- and holo-forms, complexed with chenodeoxycholate, has been determined via (15)N relaxation and steady state heteronuclear (15)N((1)H) nuclear Overhauser effect measurements. The dynamic parameters were obtained at two pH values (5.6 and 7.0) for the apoprotein and at pH 7.0 for the holoprotein, using the model free approach. Relaxation studies, performed at three different magnetic fields, revealed a substantial conformational flexibility on the microsecond to millisecond time scales, mainly localized in the C-terminal face of the beta-barrel. The observed dynamics are primarily caused by the protonation/deprotonation of a buried histidine residue, His(98), located on this flexible face. A network of polar buried side chains, defining a spine going from the E to J strand, is likely to provide the long range connectivity needed to communicate motion from His(98) to the EF loop region. NMR data are accompanied by molecular dynamics simulations, suggesting that His(98) protonation equilibrium is the triggering event for the modulation of a functionally important motion, i.e. the opening/closing at the protein open end, whereas ligand binding stabilizes one of the preexisting conformations (the open form). The results presented here, complemented with an analysis of proteins belonging to the intracellular lipid-binding protein family, are consistent with a model of allosteric activation governing the binding mechanism. The functional role of this mechanism is thoroughly discussed within the framework of the mechanism for the enterohepatic circulation of bile acids.

A decoy set for the thermostable subdomain from chicken villin headpiece, comparison of different free energy estimators

BACKGROUND

Estimators of free energies are routinely used to judge the quality of protein structural models. As these estimators still present inaccuracies, they are frequently evaluated by discriminating native or native-like conformations from large ensembles of so-called decoy structures.

RESULTS

A decoy set is obtained from snapshots taken from 5 long (100 ns) molecular dynamics (MD) simulations of the thermostable subdomain from chicken villin headpiece. An evaluation of the energy of the decoys is given using: i) a residue based contact potential supplemented by a term for the quality of dihedral angles; ii) a recently introduced combination of four statistical scoring functions for model quality estimation (FRST); iii) molecular mechanics with solvation energy estimated either according to the generalized Born surface area (GBSA) or iv) the Poisson-Boltzmann surface area (PBSA) method.

CONCLUSION

The decoy set presented here has the following features which make it attractive for testing energy scoring functions:1) it covers a broad range of RMSD values (from less than 2.0 A to more than 12 A);2) it has been obtained from molecular dynamics trajectories, starting from different non-native-like conformations which have diverse behaviour, with secondary structure elements correctly or incorrectly formed, and in one case folding to a native-like structure. This allows not only for scoring of static structures, but also for studying, using free energy estimators, the kinetics of folding;3) all structures have been obtained from accurate MD simulations in explicit solvent and after molecular mechanics (MM) energy minimization using an implicit solvent method. The quality of the covalent structure therefore does not suffer from steric or covalent problems. The statistical and physical effective energy functions tested on the set behave differently when native simulation snapshots are included or not in the set and when averaging over the trajectory is performed.

Induction of unique structural changes in guanine-rich DNA regions by the triazoloacridone C-1305, a topoisomerase II inhibitor with antitumor activities

We recently reported that the antitumor triazoloacridone, compound C-1305, is a topoisomerase II poison with unusual properties. In this study we characterize the DNA interactions of C-1305 in vitro, in comparison with other topoisomerase II inhibitors. Our results show that C-1305 binds to DNA by intercalation and possesses higher affinity for GC- than AT-DNA as revealed by surface plasmon resonance studies. Chemical probing with DEPC indicated that C-1305 induces structural perturbations in DNA regions with three adjacent guanine residues. Importantly, this effect was highly specific for C-1305 since none of the other 22 DNA interacting drugs tested was able to induce similar structural changes in DNA. Compound C-1305 induced stronger structural changes in guanine triplets at higher pH which suggested that protonation/deprotonation of the drug is important for this drug-specific effect. Molecular modeling analysis predicts that the zwitterionic form of C-1305 intercalates within the guanine triplet, resulting in widening of both DNA grooves and aligning of the triazole ring with the N7 atoms of guanines. Our results show that C-1305 binds to DNA and induces very specific and unusual structural changes in guanine triplets which likely plays an important role in the cytotoxic and antitumor activity of this unique compound.

Functional analysis of a novel RUNX2 missense mutation found in a family with cleidocranial dysplasia

Mutations of the RUNX2 gene result in dominantly inherited cleidocranial dysplasia (CCD). RUNX2 encodes for an osteoblast-specific transcription factor, which recognizes specific DNA sequences by the runt domain. DNA binding is stabilized by the interaction with the protein CBFbeta, which induces structural modifications of the runt domain. A novel 574G > A RUNX2 missense mutation has been found in members of a family clinically diagnosed with CCD. This mutation causes the glycine at position 192 to change to arginine (G192R), in loop 9 of the runt domain. Unlike other residues of loop 9, G192 does not establish DNA contacts. Accordingly, the G192R mutant showed a 50% reduction in binding activity compared to the wild-type runt domain. However, the mutation completely abolished the activating properties of the protein on osteocalcin promoter. Moreover, the G192R mutant exerts a dominant-negative effect when overexpressed. Computer modeling indicated that the G192R mutation perturbs not only loop 9, but also other parts of the runt domain, suggesting impairment of the interaction with CBFbeta.