Homology Modeling Based Solution Structure of Hoxc8-DNA Complex: Role of Context Bases Outside TAAT Stretch

Abstract

The 3D structure of neither Hoxc8 nor Hoxc8-DNA complex is known. The repressor protein Hoxc8 binds to the TAAT stretch of the promoter of the osteopontin gene and modulates its expression. Over expression of the osteopontin gene is related to diseases like osteoporosis, multiple sclerosis, cancer et cetera. In this paper we have proposed a 3D structure of Hoxc8-DNA complex obtained by Homology modeling and molecular dynamics (MD) simulation in explicit water. The crystal structure (9ant.pdb) of Antennapedia homeodomain in complex with its DNA sequence was chosen as the template based on (i) high sequence identity (85% for the protein and 60% for the DNA) and (ii) the presence of the TAAT stretch in interaction with the protein. The resulting model was refined by MD simulation for 2.0ns in explicit water. This refined model was then characterized in terms of the structural and the interactional features to improve our understanding of the mechanism of Hoxc8-DNA recognition. The interaction pattern shows that the residues Ile¹⁹⁵, Gln¹⁹⁸, and Asn¹⁹⁹, and the bases S2-⁴TAATG⁸ are most important for recognition suggesting the stretch TAATG as the ‘true recognition element’ in the present case. A strong and long-lived water bridge connecting Gln¹⁹⁸ and the base of S1-C⁷ complementary to S2-G⁸ was observed. Our predicted model of Hoxc8-DNA complex provides us with features that are consistent with the available experimental data on Hoxc8 and the general features of other homeodomain-DNA complexes. The predictions based on the model are also amenable to experimental verification.     

Molecular dynamics analysis of the engrailed homeodomain-DNA recognition

Molecular dynamics (MD) simulations were performed for investigating the role of Gln50 in the engrailed homeodomain-DNA recognition. Employing the crystal structure of free engrailed homeodomain and homeodomain-DNA complex as a starting structure, we carried out MD simulations of: (i) the complex between engrailed homeodomain and a 20 base-pair DNA containing TAATTA core sequence; (ii) the free engrailed homeodomain. The simulations show that homeodomain flexibility does not depend on its ligation state. The engrailed homeodomain shows similar flexibility, and the recognition helix-3 shows very similar characteristic of high rigidity and limited conformational space in two complexation states. At the same time, DNA structure has also no obvious conformational fluctuations. These results preclude the possibility of the side chain of Gln50 forming direct hydrogen bonds to the core DNA bases. MD simulations confirm a few well-conserved sites for water-mediated hydrogen bonds from protein to DNA are occupied by water molecules, and Gln50 interacts with corresponding core DNA bases through water-mediated hydrogen bonds. So Gln50 plays a relatively modest role in determining the affinity and specificity of the engrailed homeodomain. In addition, the electrostatic interaction between homeodomain and phosphate backbone of the DNA is a main factor for N- and C-terminal arm becoming ordered upon DNA binding.     

Role of N-terminal helix in interaction of ribosomal protein S15 with 16S rRNA

The position and conformation of the N-terminal helix of free ribosomal protein S15 was earlier found to be modified under various conditions. This variability was supposed to provide the recognition by the protein of its specific site on 16S rRNA. To test this hypothesis, we substituted some amino acid residues in this helix and assessed effects of these substitutions on the affinity of the protein for 16S rRNA. The crystal structure of the complex of one of these mutants (Thr3Cys S15) with the 16S rRNA fragment was determined, and a computer model of the complex containing another mutant (Gln8Met S15) was designed. The available and new information was analyzed in detail, and the N-terminal helix was concluded to play no significant role in the specific binding of the S15 protein to its target on 16S rRNA.     

Crystal structure of an engrailed homeodomain-DNA complex at 2.8 A resolution: a framework for understanding homeodomain-DNA interactions

The crystal structure of a complex containing the engrailed homeodomain and a duplex DNA site has been determined at 2.8 A resolution and refined to a crystallographic R factor of 24.4%. In this complex, two separate regions of the 61 amino acid polypeptide contact a TAAT subsite. An N-terminal arm fits into the minor groove, and the side chains of Arg-3 and Arg-5 make contacts near the 5' end of this "core consensus" binding site. An alpha helix fits into the major groove, and the side chains of IIe-47 and Asn-51 contact base pairs near the 3' end of the TAAT site. This "recognition helix" is part of a structurally conserved helix-turn-helix unit, but these helices are longer than the corresponding helices in the lambda repressor, and the relationship between the helix-turn-helix unit and the DNA is significantly different.     

The third helix of the murine Hoxc8 homeodomain facilitates protein transduction in mammalian cells

Previously, we have demonstrated that purified Hoxc8 homeoprotein has the ability to penetrate the cellular membrane and can be transduced efficiently into COS-7 cells. Moreover, the Hoxc8 protein is able to form a complex with DNA molecules in vitro and helps the DNA be delivered intracellularly, serving as a gene delivery vehicle. Here, we further analyzed the membrane transduction activity of Hoxc8 protein and provide the evidence that the 16 amino acid (a.a.191-206, 2.23 kDa) third helix of murine Hoxc8 protein is an efficient protein transduction domain (PTD). When the 16 amino acid peptide was fused at the carboxyl terminal of enhanced green fluorescence protein (EGFP), the fusion proteins were transduced efficiently into the primary pig fetal fibroblast cells. The transduction efficiency increased in a concentration-dependent manner up to 1 μM, and appeared to plateau above a concentration of 1 μM. When tandem multimers of PTD, EGFP-PTD(2), EGFP-PTD(3), EGFP-PTD(4), and EGFP-PTD(5), were analyzed at 500 nM of concentration, the penetrating efficiency increased in a dose-dependent manner. As the number of PTDs increased, the EGFP signal also increased, although the signal maintained plateau after EGFP-PTD(3). These results indicate that the 16 amino acid third helix is the key element responsible for the membrane transduction activity of Hoxc8 proteins, and further suggest that the small peptide could serve as a therapeutic delivery vehicle for large cargo proteins.     

Comparative molecular dynamics simulations of HIV-1 integrase and the T66I/M154I mutant: binding modes and drug resistance to a diketo acid inhibitor

HIV-1 IN is an essential enzyme for viral replication and an interesting target for the design of new pharmaceuticals for use in multidrug therapy of AIDS. L-731,988 is one of the most active molecules of the class of beta-diketo acids. Individual and combined mutations of HIV-1 IN at residues T66, S153, and M154 confer important degrees of resistance to one or more inhibitors belonging to this class. In an effort to understand the molecular mechanism of the resistance of T66I/M154I IN to the inhibitor L-731,988 and its specific binding modes, we have carried out docking studies, explicit solvent MD simulations, and binding free energy calculations. The inhibitor was docked against different protein conformations chosen from prior MD trajectories, resulting in 2 major orientations within the active site. MD simulations have been carried out for the T66I/M154I DM IN, DM IN in complex with L-731,988 in 2 different orientations, and 1QS4 IN in complex with L-731,988. The results of these simulations show a similar dynamical behavior between T66I/M154I IN alone and in complex with L-731,988, while significant differences are observed in the mobility of the IN catalytic loop (residues 138-149). Water molecules bridging the inhibitor to residues from the active site have been identified, and residue Gln62 has been found to play an important role in the interactions between the inhibitor and the protein. This work provides information about the binding modes of L-731,988, as well as insight into the mechanism of inhibitor-resistance in HIV-1 integrase.     

Comparison of protein solution structures refined by molecular dynamics simulation in vacuum,with a generalized Born model,and with explicit water

The inclusion of explicit solvent water in molecular dynamics refinement of NMR structures ought to provide the most physically meaningful accounting for the effects of solvent on structure, but is computationally expensive. In order to evaluate the validity of commonly used vacuum refinements and of recently developed continuum solvent model methods, we have used three different methods to refine a set of NMR solution structures of a medium sized protein, Escherichia coliglutaredoxin 2, from starting structures calculated using the program DYANA. The three different refinement protocols used molecular dynamics simulated annealing with the program AMBER in vacuum (VAC), including a generalized Born (GB) solvent model, and a full calculation including explicit solvent water (WAT). The structures obtained using the three methods of refinements were very similar, a reflection of their generally well-determined nature. However, the structures refined with the generalized Born model were more similar to those from explicit water refinement than those refined in vacuum. Significant improvement was seen in the percentage of backbone dihedral angles in the most favored regions of , space and in hydrogen bond pattern for structures refined with the GB and WAT models, compared with the structures refined in vacuum. The explicit water calculation took an average of 200 h of CPU time per structure on an SGI cluster, compared to 15–90 h for the GB calculation (depending on the parameters used) and 2 h for the vacuum calculation. The generalized Born solvent model proved to be an excellent compromise between the vacuum and explicit water refinements, giving results comparable to those of the explicit water calculation. Some improvement for and angle distribution and hydrogen bond pattern can also be achieved by energy minimizing the vacuum structures with the GB model, which takes a much shorter time than MD simulations with the GB model.     

Hoxc8 early enhancer of the Indonesian coelacanth, Latimeria menadoensis

Hoxc8 early enhancer controls the initiation and establishment phase of Hoxc8 expression in the mouse. Comparative studies indicate the presence of Hoxc8 early enhancer sequences in different vertebrate clades including mammals, birds and fish. Previous studies have shown differences between teleost and mammalian Hoxc8 early enhancers with respect to sequence and organization of protein binding elements. This raises the question of when the Hoxc8 early enhancer arose and how it has become modified in different vertebrate lineages. Here, we describe Hoxc8 early enhancer from the Indonesian coelacanth, Latimeria menadoensis. Coelacanths are the only extant lobefinned fish whose genome is tractable to genome analysis. The Latimeria Hoxc8 early enhancer sequence more closely resembles that of the mouse than that of Fugu or zebrafish. When assayed for enhancer activity by reporter gene analysis in transgenic mouse embryos, Latimeria Hoxc8 early enhancer directs expression to the posterior neural tube and mesoderm similar to that of the mouse enhancer. These observations support a close relationship between coelacanths and tetrapods and place the origin of a common Hoxc8 early enhancer sequence within the sarcopterygian lineage. The divergence of teleost (actinopterygii) Hoxc8 early enhancer may reflect a case of relaxed selection or other forms of instability induced by genome duplication events.     

Collective properties of hydration: long range and specificity of hydrophobic interactions.

We report results of molecular dynamics (MD) simulations of composite model solutes in explicit molecular water solvent, eliciting novel aspects of the recently demonstrated, strong many-body character of hydration. Our solutes consist of identical apolar (hydrophobic) elements in fixed configurations. Results show that the many-body character of PMF is sufficiently strong to cause 1) a remarkable extension of the range of hydrophobic interactions between pairs of solute elements, up to distances large enough to rule out pairwise interactions of any type, and 2) a SIF that drives one of the hydrophobic solute elements toward the solvent rather than away from it. These findings complement recent data concerning SIFs on a protein at single-residue resolution and on model systems. They illustrate new important consequences of the collective character of hydration and of PMF and reveal new aspects of hydrophobic interactions and, in general, of SIFs. Their relevance to protein recognition, conformation, function, and folding and to the observed slight yet significant nonadditivity of functional effects of distant point mutations in proteins is discussed. These results point out the functional role of the configurational and dynamical states (and related statistical weights) corresponding to the complex configurational energy landscape of the two interacting systems: biomolecule + water.     

Conformation of a trimannoside bound to mannose-binding protein by nuclear magnetic resonance and molecular dynamics simulations

A model of the carbohydrate recognition domain of the serum form of mannose-binding protein (MBP) from rat complexed with methyl 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside is presented. Allowed conformations for the bound sugar were derived from simulated annealing protocols incorporating distance restraints computed from transferred NOESY spectra. The resulting sugar conformations were then modeled into the MBP binding site, and these models of the complex were refined using molecular dynamics (MD) simulations in the presence of solvent water. These studies indicate that only one of the two major conformations of the alpha(1-->6) linkage found in solution is significantly populated in the bound state (omega = 60 degrees ), whereas the alpha(1-->3) linkage samples at least two states, similar to its behavior in free solution. The bound conformation allows direct hydrogen bonds to form between the sugar and K182 of MBP, in addition to other water-mediated hydrogen bonds. Estimates of binding constants of candidate complexes based on changes in solvent-accessible surface areas upon binding support the NMR and MD results. These estimates further suggest that the enthalpic gains of the additional sugar-MBP interactions in a trisaccharide as opposed to a monosaccharide are offset by entropic penalties, offering an explanation for previous binding data.     

A novel protocol of energy optimisation for predicted protein structures built by homology modelling

Homology modelling was applied to predict the three-dimensional (3D) structures of six sets of lipase proteins. Sequence identities between the target and template were 34.6, 44.9, 57.4, 69.9, 79.0 and 86.2%, respectively. Then, eight different protocols including three optimising factors [periodically bounded cell (PBC) water, molecular dynamics (MD) simulation, ‘grade-unpacking’ strategy or ‘combinatorial’ strategy] were used to refine the initial model of each system. By comparing the energy-optimised models with the true 3D structure of the target protein in terms of all backbone atoms' root mean square deviation, we determined a novel but all-purpose protocol for model refinement. The protocol refined a homology model by adopting the ‘grade-unpacking’ strategy for energy minimisation while the model was solvated in PBC water. Furthermore, by comparing the influence of each single optimising factor on the accuracy of the refined structure, we found that introducing the MD simulation into the model refinement method would decrease the accuracy of the final protein structure while methods with either PBC water or the ‘grade-unpacking’ strategy would increase the accuracy of the final model.     

Structure at 2.3 A resolution of the cytochrome bc(1) complex from the yeast Saccharomyces cerevisiae co-crystallized with an antibody Fv fragment

BACKGROUND: The cytochrome bc(1) complex is part of the energy conversion machinery of the respiratory and photosynthetic electron transfer chains. This integral membrane protein complex catalyzes electron transfer from ubiquinol to cytochrome c. It couples the electron transfer to the electrogenic translocation of protons across the membrane via a so-called Q cycle mechanism. RESULTS: The cytochrome bc(1) complex from the yeast Saccharomyces cerevisiae was crystallized together with a bound antibody Fv fragment. The structure was determined at 2.3 A resolution using multiple isomorphous replacement, and refined to a crystallographic R factor of 22.2% (R(free) = 25.4%). The complex is present as a homodimer. Each 'monomer' of the refined model includes 2178 amino acid residues of subunits COR1, QCR2, COB, CYT1, RIP1, QCR6, QCR7, QCR8 and QCR9 of the cytochrome bc(1) complex and of the polypeptides V(H) and V(L) of the Fv fragment, the cofactors heme b(H), heme b(L), heme c(1), the [2Fe-2S] cluster and 346 water molecules. The Fv fragment binds to the extrinsic domain of the [2Fe-2S] Rieske protein and is essential for formation of the crystal lattice. CONCLUSIONS: The approach to crystallize membrane proteins as complexes with specific antibody fragments appears to be of general importance. The structure of the yeast cytochrome bc(1) complex reveals in detail the binding sites of the natural substrate coenzyme Q6 and the inhibitor stigmatellin. Buried water molecules close to the binding sites suggest possible pathways for proton uptake and release. A comparison with other cytochrome bc(1) complexes shows features that are specific to yeast.