A structural model of the E. coli PhoB Dimer in the transcription initiation complex

<p>Abstract</p> <p>Background</p> <p>There exist > 78,000 proteins and/or nucleic acids structures that were determined experimentally. Only a small portion of these structures corresponds to those of protein complexes. While homology modeling is able to exploit knowledge-based potentials of side-chain rotomers and backbone motifs to infer structures for new proteins, no such general method exists to extend our understanding of protein interaction motifs to novel protein complexes.</p> <p>Results</p> <p>We use a Motif Binding Geometries (MBG) approach, to infer the structure of a protein complex from the database of complexes of homologous proteins taken from other contexts (such as the helix-turn-helix motif binding double stranded DNA), and demonstrate its utility on one of the more important regulatory complexes in biology, that of the RNA polymerase initiating transcription under conditions of phosphate starvation. The modeled PhoB/RNAP/σ-factor/DNA complex is stereo-chemically reasonable, has sufficient interfacial Solvent Excluded Surface Areas (SESAs) to provide adequate binding strength, is physically meaningful for transcription regulation, and is consistent with a variety of known experimental constraints.</p> <p>Conclusions</p> <p>Based on a straightforward and easy to comprehend concept, "proteins and protein domains that fold similarly could interact similarly", a structural model of the PhoB dimer in the transcription initiation complex has been developed. This approach could be extended to enable structural modeling and prediction of other bio-molecular complexes. Just as models of individual proteins provide insight into molecular recognition, catalytic mechanism, and substrate specificity, models of protein complexes will provide understanding into the combinatorial rules of cellular regulation and signaling.</p

Structural modeling of the catalytic subunit-regulatory subunit dimeric complex of the camp-dependent protein kinase.

The cAMP-dependent protein kinase (PKA) is a multifunctional kinase that serves as a prototype for understanding second messenger signaling and protein phosphorylation. In the absence of a cAMP signal, PKA exists as a dimer of dimers, consisting of two regulatory (R) and two catalystic (C) subunits. Based on experimentally derived data (i.e., crystal structures of the R and C subunits, mutagenesis data identifying points of subunit-subunit contacts), the neutron scattering derived model for the heterodimer (Zhao et al., 1998) and using a set of computational approaches (homology modeling, Monte Carlo simulation), they have developed a high-resolution model of the RII{alpha}-C{alpha} dimer. The nature of the subunit-subunit interface was studied. The model reveals an averaged size dimer interface (2100 Angstrom{sup 2}) that is distant from the pseudo-substrate binding site on the C subunit. The additional contacts made by the pseudosubstrate increases the stability of the dimeric complex. Based on a set of R-C dimer structures derived using a simulated annealing approach, specific interactions (hydrogen bonds) between the two subunits and were identified

Mismatched dNTP incorporation by DNA polymerase β does not proceed via globally different conformational pathways†

Understanding how DNA polymerases control fidelity requires elucidation of the mechanisms of matched and mismatched dNTP incorporations. Little is known about the latter because mismatched complexes do not crystallize readily. In this report, we employed small-angle X-ray scattering (SAXS) and structural modeling to probe the conformations of different intermediate states of mammalian DNA polymerase β (Pol β) in its wild-type and an error-prone variant, I260Q. Our structural results indicate that the mismatched ternary complex lies in-between the open and the closed forms, but more closely resembles the open form for WT and the closed form for I260Q. On the basis of molecular modeling, this over-stabilization of mismatched ternary complex of I260Q is likely caused by formation of a hydrogen bonding network between the side chains of Gln260, Tyr296, Glu295 and Arg258, freeing up Asp192 to coordinate MgdNTP. These results argue against recent reports suggesting that mismatched dNTP incorporations follow a conformational path distinctly different from that of matched dNTP incorporation, or that its conformational closing is a major contributor to fidelity

Design and synthesis of new 2-arylnaphthyridin-4-ones as potent antitumor agents targeting tumorigenic cell lines

To develop new anticancer drug candidates from 2-arylnaphthyridin-4-one (AN), we have designed and synthesized a series of 3′-hydroxy and 6-hydroxy derivatives of AN. The results of cytotoxicity screening indicated that the replacement of the 3′-methoxy moiety on the C-ring phenyl group of AN (6a–e) with 3′-hydroxy (7a–e) made no significant effect on the inhibitory activity against HL-60, Hep3B and NCI-H460 cancer cell lines. On the other hand, replacing the 6-methoxy group on the A-ring of AN (6g–i) with a 6-hydroxy group (7g–i) resulted in reduced inhibitory activity against the above three cancer cell lines. Among the above-mentioned target compounds, 2-(3-hydroxyphenyl)-5-methyl-1,8-naphthyridin-4(1H)-one (7a) demonstrated the greatest potency and the best selectivity toward tumorigenic cancer cell lines. In a 7a preliminary mechanism of action study in Hep3B hepatoma cells, 7a showed the effects on microtubules followed by cell cycle arrest and sequentially led to apoptosis

A reduced-coordinate approach to modeling RNA 3-D structures

With the realization of RNA molecules capable of performing very specific functions (e.g., catalytic RNAs and RNAs that bind ligand with affinity and specificity of an anti-body) and contrary to the traditional view that structure of RNA molecules being functionally passive, it has become clear that studying the 3-dimensional (3-D) folding of RNA molecules is a very important task. In the absence of sufficient number of experimentally determined RNA structures available up-to-date, folding of RNA structures computationally provides an alternative approach in studying the 3-D structure of RNA molecules. We have developed a computational approach for folding RNA 3-D structures. The method is conceptually simple and general. It consists of two major components. The first being the arrangement of all helices in space. Once the helices are positioned and oriented in space, structures of the connecting loops are modeled and inserted between the helices. Any number of structural constraints derived either experimentally or theoretically can be used to guide the folding processes. A conformational sampling approach is developed with structural equilibration using the Metropolis Monte Carlo simulation. The lengths of various loop sizes (ranging from 1 base to 7 bases) are calculated based on a set of RNA structures deposited in PDB as well as a set of loop structures constructed using our method. The validity of using the averaged loop lengths of the connecting loops as distance constraints for arranging the helices in space is studied

A Computational Approach to Modeling Nucleic Acid Hairpin Structures

Hairpin is a structural motif frequently observed in both RNA and DNA molecules. This motif is involved specifically in various biological functions (e.g., gene expression and regulation). To understand how these hairpin motifs perform their functions, it is important to study their structures. Compared to protein structural motifs, structures of nucleic acid hairpins are less known. Based on a set of reduced coordinates for describing nucleic acid structures and a sampling algorithm that equilibrates structures using Metropolis Monte Carlo simulation, we developed a method to model nucleic acid hairpin structures. This method was used to predict the structure of a DNA hairpin with a single-guanosine loop. The lowest energy structure from the ensemble of 200 sampled structures has a RMSD of <1.5 Å, from the structure determined using NMR. Additional constraints for the loop bases were introduced for modeling an RNA hairpin with two nucleotides in the loop. The modeled structure of this RNA hairpin has extensive base stacking and an extra hydrogen bond (between the CYT in the loop and a phosphate oxygen), as observed in the NMR structure

Three-dimensional model of a selective theophylline-binding RNA molecule

We propose a three-dimensional (3D) model for an RNA molecule that selectively binds theophylline but not caffeine. This RNA, which was found using SELEX [Jenison, R.D., et al., Science (1994) 263:1425] is 10,000 times more specific for theophylline (Kd=320 nM) than for caffeine (Kd=3.5 mM), although the two ligands are identical except for a methyl group substituted at N7 (present only in caffeine). The binding affinity for ten xanthine-based ligands was used to derive a Comparative Molecular Field Analysis (CoMFA) model (R{sup 2} = 0.93 for 3 components, with cross-validated R{sup 2} of 0.73), using the SYBYL and GOLPE programs. A pharmacophoric map was generated to locate steric and electrostatic interactions between theophylline and the RNA binding site. This information was used to identify putative functional groups of the binding pocket and to generate distance constraints. Based on a model for the secondary structure (Jenison et al., idem), the 3D structure of this RNA was then generated using the following method: each helical region of the RNA molecule was treated as a rigid body; single-stranded loops with specific end-to-end distances were generated. The structures of RNA-xanthine complexes were studied using a modified Monte Carlo algorithm. The detailed structure of an RNA-ligand complex model, as well as possible explanations for the theophylline selectivity will be discussed