Intrinsically disordered electronegative clusters improve stability and binding specificity of RNA-binding proteins.

RNA-binding proteins play crucial roles in various cellular functions and contain abundant disordered protein regions. The disordered regions in RNA-binding proteins are rich in repetitive sequences, such as poly-K/R, poly-N/Q, poly-A, and poly-G residues. Our bioinformatic analysis identified a largely neglected repetitive sequence family we define as electronegative clusters (ENCs) that contain acidic residues and/or phosphorylation sites. The abundance and length of ENCs exceed other known repetitive sequences. Despite their abundance, the functions of ENCs in RNA-binding proteins are still elusive. To investigate the impacts of ENCs on protein stability, RNA-binding affinity, and specificity, we selected one RNA-binding protein, the ribosomal biogenesis factor 15 (Nop15), as a model. We found that the Nop15 ENC increases protein stability and inhibits nonspecific RNA binding, but minimally interferes with specific RNA binding. To investigate the effect of ENCs on sequence specificity of RNA binding, we grafted an ENC to another RNA-binding protein, Ser/Arg-rich splicing factor 3. Using RNA Bind-n-Seq, we found that the engineered ENC inhibits disparate RNA motifs differently, instead of weakening all RNA motifs to the same extent. The motif site directly involved in electrostatic interaction is more susceptible to the ENC inhibition. These results suggest that one of functions of ENCs is to regulate RNA binding via electrostatic interaction. This is consistent with our finding that ENCs are also overrepresented in DNA-binding proteins, whereas underrepresented in halophiles, in which nonspecific nucleic acid binding is inhibited by high concentrations of salts.

Structural and Metal Ion Effects on Human Topoisomerase IIα Inhibition by α-(N)-Heterocyclic Thiosemicarbazones.

Our previous research has shown that α-(N)-heterocyclic thiosemicarbazone (TSC) metal complexes inhibit human topoisomerase IIα (TopoIIα), while the ligands without metals do not. To find out the structural elements of TSC that are important for inhibiting TopoIIα, we have synthesized two series of α-(N)-heterocyclic TSCs with various substrate ring segments, side chain substitutions, and metal ions, and we have examined their activities in TopoIIα-mediated plasmid DNA relaxation and cleavage assays. Our goal is to explore the structure-activity relationship of α-(N)-heterocyclic TSCs and their effect on TopoIIα. Our data suggest that, similar to Cu(II)-TSCs, Pd(II)-TSC complexes inhibit plasmid DNA relaxation mediated by TopoIIα. In TopoIIα-mediated plasmid DNA cleavage assays, the Cu(II)-TSC complexes induce higher levels of DNA cleavage than their Pd(II) counterparts. The Cu(II)-TSC complexes with methyl, ethyl, and tert-butyl substitutions are slightly more effective than those with benzyl and phenyl groups. The α-(N)-heterocyclic ring substrates of the TSCs, including benzoylpyridine, acetylpyridine, and acetylthiazole, do not exhibit a significant difference in TopoIIα-mediated DNA cleavage. Our data suggest that the metal ion of TSC complexes plays a predominant role in inhibition of TopoIIα, the side chain substitution of the terminal nitrogen plays a secondary role, while the substrate ring segment has the least effect. Our molecular modeling data support the biochemical data, which together provide a mechanism by which Cu(II)-TSC complexes stabilize TopoIIα-mediated cleavage complexes.

Molecular interactions between photosystem I and ferredoxin: an integrated energy frustration and experimental model.

The stromal domain (PsaC, PsaD, and PsaE) of photosystem I (PSI) reduces transiently bound ferredoxin (Fd) or flavodoxin. Experimental structures exist for all of these protein partners individually, but no experimental structure of the PSI/Fd or PSI/flavodoxin complexes is presently available. Molecular models of Fd docked onto the stromal domain of the cyanobacterial PSI site are constructed here utilizing X-ray and NMR structures of PSI and Fd, respectively. Predictions of potential protein-protein interaction regions are based on experimental site-directed mutagenesis and cross-linking studies to guide rigid body docking calculations of Fd into PSI, complemented by energy landscape theory to bring together regions of high energetic frustration on each of the interacting proteins. The results identify two regions of high localized frustration on the surface of Fd that contain negatively charged Asp and Glu residues. This study predicts that these regions interact predominantly with regions of high localized frustration on the PsaC, PsaD, and PsaE chains of PSI, which include several residues predicted by previous experimental studies.

Biophysical characterization of an ensemble of intramolecular i-motifs formed by the human c-MYC NHE III1 P1 promoter mutant sequence.

i-Motif-forming sequences are present in or near the regulatory regions of >40% of all genes, including known oncogenes. We report here the results of a biophysical characterization and computational study of an ensemble of intramolecular i-motifs that model the polypyrimidine sequence in the human c-MYC P1 promoter. Circular dichroism results demonstrate that the mutant sequence (5'-CTT TCC TAC CCTCCC TAC CCT AA-3') can adopt multiple "i-motif-like," classical i-motif, and single-stranded structures as a function of pH. The classical i-motif structures are predominant in the pH range 4.2-5.2. The "i-motif-like" and single-stranded structures are the most significant species in solution at pH higher and lower, respectively, than that range. Differential scanning calorimetry results demonstrate an equilibrium mixture of at least three i-motif folded conformations with Tm values of 38.1, 46.6, and 49.5 degrees C at pH 5.0. The proposed ensemble of three folded conformations includes the three lowest-energy conformations obtained by computational modeling and two folded conformers that were proposed in a previous NMR study. The NMR study did not report the most stable conformer found in this study.

General library-based Monte Carlo technique enables equilibrium sampling of semi-atomistic protein models.

We introduce "library-based Monte Carlo" (LBMC) simulation, which performs Boltzmann sampling of molecular systems based on precalculated statistical libraries of molecular-fragment configurations, energies, and interactions. The library for each fragment can be Boltzmann distributed and thus account for all correlations internal to the fragment. LBMC can be applied to both atomistic and coarse-grained models, as we demonstrate in this "proof-of-principle" report. We first verify the approach in a toy model and in implicitly solvated all-atom polyalanine systems. We next study five proteins, up to 309 residues in size. On the basis of atomistic equilibrium libraries of peptide-plane configurations, the proteins are modeled with fully atomistic backbones and simplified Go-like interactions among residues. We show that full equilibrium sampling can be obtained in days to weeks on a single processor, suggesting that more accurate models are well within reach. For the future, LBMC provides a convenient platform for constructing adjustable or mixed-resolution models: the configurations of all atoms can be stored at no run-time cost, while an arbitrary subset of interactions is "turned on".

Binding Mechanisms of Electron Transport Proteins with Cyanobacterial Photosystem I: An Integrated Computational and Experimental Model.

The stromal domain (PsaC, D, and E) of photosystem I (PSI) in cyanobacteria accepts electrons from PsaA and PsaB of photosystem I (PSI). These electrons are then used in the reduction of transiently bound ferredoxin (Fd) or flavodoxin. Experimental X-ray and NMR structures are known for all of these protein partners separately, yet to date, there is no known experimental structure of the PSI/Fd complexes in the published literature. Computational models of Fd docked with the stromal domain of cyanobacterial PSI were assembled here starting from X-ray and NMR structures of PSI and Fd. Predicted models of specific regions of protein-protein interactions were built on the basis of energetic frustration, residue conservation and evolutionary importance, as well as from experimental site-directed mutagenesis and cross-linking studies. Microsecond time-scale molecular dynamics simulations of the PSI/Fd complexes suggest, rather than a single complex structure, the possible existence of multiple transient complexes of Fd bound to PSI.

Hydropathic analysis of the free energy differences in anthracycline antibiotic binding to DNA.

Molecular models of six anthracycline antibiotics and their complexes with 32 distinct DNA octamer sequences were created and analyzed using HINT (Hydropathic INTeractions) to describe binding. The averaged binding scores were then used to calculate the free energies of binding for comparison with experimentally determined values. In parsing our results based on specific functional groups of doxorubicin, our calculations predict a free energy contribution of -3.6 +/- 1.1 kcal x mol(-1) (experimental -2.5 +/- 0.5 kcal x mol(-1)) from the groove binding daunosamine sugar. The net energetic contribution of removing the hydroxyl at position C9 is -0.7 +/- 0.7 kcal x mol(-1) (-1.1 +/- 0.5 kcal x mol(-1)). The energetic contribution of the 3' amino group in the daunosamine sugar (when replaced with a hydroxyl group) is -3.7 +/- 1.1 kcal x mol(-1) (-0.7 +/- 0.5 kcal x mol(-1)). We propose that this large discrepancy may be due to uncertainty in the exact protonation state of the amine. The energetic contribution of the hydroxyl group at C14 is +0.4 +/- 0.6 kcal x mol(-1) (-0.9 +/- 0.5 kcal x mol(-1)), largely due to unfavorable hydrophobic interactions between the hydroxyl oxygen and the methylene groups of the phosphate backbone of the DNA. Also, there appears to be considerable conformational uncertainty in this region. This computational procedure calibrates our methodology for future analyses where experimental data are unavailable.

A new bisintercalating anthracycline with picomolar DNA binding affinity.

A new bisintercalating anthracycline (WP762) has been designed, in which monomeric units of daunorubicin have been linked through their amino groups on the daunosamine moieties using an m-xylenyl linker. Differential scanning calorimetry and UV melting experiments were used to measure the ultratight binding of WP762 to DNA. The binding constant for the interaction of WP762 with herring sperm DNA was determined to be 7.3 (+/-0.2) x 10(12) M(-1) at 20 degrees C. The large favorable binding free energy of -17.3 kcal mol(-1) was found to result from a large negative enthalpic contribution of -33.8 kcal mol(-1) and an opposing entropic term (-TDeltaS = +16.5 kcal mol(-1)). A comparative molecular modeling study rationalized the increased binding by the m-xylenyl linker of WP762 positioning in the DNA minor groove compared to the p-xylenyl linker found in WP631, the first bis-anthracycline of this type. The cytotoxicity of WP762 was compared to that of other anthracyclines in Jurkat T lymphocytes. These studies, together with an analysis of the cell-cycle traverse in the presence of WP762, suggest that in these cells the new drug is more cytotoxic than the structurally related WP631.

A computational model for anthracycline binding to DNA: tuning groove-binding intercalators for specific sequences.

Anthracycline antibiotics such as doxorubicin and its analogues have been in common use as anticancer drugs for almost half a century. There has been intense interest in the DNA binding sequence specificity of these compounds in recent years, with the hope that a compound could be identified that could possibly modulate gene expression or exhibit reduced toxicity. To computationally analyze this phenomenon, we have constructed molecular models of 65 doxorubicin analogues and their complexes with eight distinct DNA octamer sequences. The HINT (Hydropathic INTeractions) program was utilized to describe binding, including differences in the functional group contributions as well as sequence selectivity. Of these 65 compounds, two compounds were calculated to have a selectivity (the calculated DeltaDeltaG(sel) between the sequence with the strongest binding and the second strongest binding sequence) greater than -0.75 kcal mol(-1) for one sequence over all others, 10 compounds were specific between -0.50 and -0.74 kcal mol(-1), 18 compounds were specific between -0.25 and -0.49 kcal mol(-1), and 35 compounds were virtually nonspecific with a DeltaDeltaG below -0.24 kcal mol(-1). Several compounds have been identified from this study that include features which may enhance sequence selectivity, including several with a halogen in lieu of the 4'-OH in the daunosamine sugar, one compound with a nonaromatic six-membered ring (pirarubicin) in place of the 4'-OH, and a compound with an aromatic ring in the vicinity of the C(14) region (zorubicin). Removal of the methoxy group at the C(4) position on the aglycone portion also appears to add potency and selectivity (idarubicin). Overall, efficient computational methods are presented that can be utilized to analyze the free energy of binding and sequence selectivity of both known and designed analogues of doxorubicin to identify future lead compounds for further experimental research.

Homology modeling of the CheW coupling protein of the chemotaxis signaling complex.

Homology models of the E. coli and T. maritima chemotaxis protein CheW were constructed to assess the quality of structural predictions and their applicability in chemotaxis research: i) a model of E. coli CheW was constructed using the T. maritima CheW NMR structure as a template, and ii) a model of T. maritima CheW was constructed using the E. coli CheW NMR structure as a template. The conformational space accessible to the homology models and to the NMR structures was investigated using molecular dynamics and Monte Carlo simulations. The results show that even though static homology models of CheW may be partially structurally different from their corresponding experimentally determined structures, the conformational space they can access through their dynamic variations can be similar, for specific regions of the protein, to that of the experimental NMR structures. When CheW homology models are allowed to explore their local accessible conformational space, modeling can provide a rational path to predicting CheW interactions with the MCP and CheA proteins of the chemotaxis complex. Homology models of CheW (and potentially, of other chemotaxis proteins) should be seen as snapshots of an otherwise larger ensemble of accessible conformational space.

Molecular modeling and biophysical analysis of the c-MYC NHE-III1 silencer element.

G-Quadruplex and i-Motif-forming sequences in the promoter regions of several oncogenes show promise as targets for the regulation of oncogenes. In this study, molecular models were created for the c-MYC NHE-III(1) (nuclease hypersensitivity element III(1)) from two 39-base complementary sequences. The NHE modeled here consists of single folded conformers of the polypurine intramolecular G-Quadruplex and the polypyrimidine intramolecular i-Motif structures, flanked by short duplex DNA sequences. The G-Quadruplex was based on published NMR structural data for the c-MYC 1:2:1 loop isomer. The i-Motif structure is theoretical (with five cytosine-cytosine pairs), where the central intercalated cytosine core interactions are based on NMR structural data obtained for a tetramolecular [d(A(2)C(4))(4)] model i-Motif. The loop structures are in silico predictions of the c-MYC i-motif loops. The porphyrin meso-tetra(N-methyl-4-pyridyl)porphine (TMPyP4), as well as the ortho and meta analogs TMPyP2 and TMPyP3, were docked to six different locations in the complete c-MYC NHE. Comparisons are made for drug binding to the NHE and the isolated G-Quadruplex and i-Motif structures. NHE models both with and without bound cationic porphyrin were simulated for 100 ps using molecular dynamics techniques, and the non-bonded interaction energies between the DNA and porphyrins calculated for all of the docking interactions.

Biophysical studies of the c-MYC NHE III1 promoter: model quadruplex interactions with a cationic porphyrin.

Regulation of the structural equilibrium of G-quadruplex-forming sequences located in the promoter regions of oncogenes by the binding of small molecules has shown potential as a new avenue for cancer chemotherapy. In this study, microcalorimetry (isothermal titration calorimetry and differential scanning calorimetry), electronic spectroscopy (ultraviolet-visible and circular dichroism), and molecular modeling were used to probe the complex interactions between a cationic porphryin mesotetra (N-methyl-4-pyridyl) porphine (TMPyP4) and the c-MYC PU 27-mer quadruplex. The stoichiometry at saturation is 4:1 mol of TMPyP4/c-MYC PU 27-mer G-quadruplex as determined by isothermal titration calorimetry, circular dichroism, and ultraviolet-visible spectroscopy. The four independent TMPyP4 binding sites fall into one of two modes. The two binding modes are different with respect to affinity, enthalpy change, and entropy change for formation of the 1:1 and 2:1, or 3:1 and 4:1 complexes. Binding of TMPyP4, at or near physiologic ionic strength ([K(+)] = 0.13 M), is described by a "two-independent-sites model." The two highest-affinity sites exhibit a K(1) of 1.6 x 10(7) M(-1) and the two lowest-affinity sites exhibit a K(2) of 4.2 x 10(5) M(-1). Dissection of the free-energy change into the enthalpy- and entropy-change contributions for the two modes is consistent with both "intercalative" and "exterior" binding mechanisms. An additional complexity is that there may be as many as six possible conformational quadruplex isomers based on the sequence. Differential scanning calorimetry experiments demonstrated two distinct melting events (T(m)1 = 74.7 degrees C and T(m)2 = 91.2 degrees C) resulting from a mixture of at least two conformers for the c-MYC PU 27-mer in solution.