How hormone receptor-DNA binding affects nucleosomal DNA: the role of symmetry

Nuclear Receptors in Liver Fibrosis

Nuclear receptors are ligand-activated transcription factors that regulate gene expression of a variety of key molecular signals involved in liver fibrosis. The primary cellular driver of liver fibrogenesis are activated hepatic stellate cells. Different NRs regulate the hepatic expression of pro-inflammatory and pro-fibrogenic cytokines that promote the transformation of hepatic stellate cells into fibrogenic myofibroblasts. Importantly, nuclear receptors regulate gene expression circuits that promote hepatic fibrogenesis and/or allow liver fibrosis regression. In this review, we highlight the direct and indirect influence of nuclear receptors on liver fibrosis, with a focus on hepatic stellate cells, and discuss potential therapeutic effects of nuclear receptor modulation in regard to anti-fibrotic and anti-inflammatory effects. Further research on nuclear receptors-related signaling may lead to the clinical development of effective anti-fibrotic therapies for patients with liver disease.

Genome Dashboards: Framework and Examples

Genomics is a sequence-based informatics science and a three-dimensional-structure-based material science. However, in practice, most genomics researchers utilize sequence-based informatics approaches or three-dimensional-structure-based material science techniques, not both. This division is, at least in part, the result of historical developments rather than a fundamental necessity. The underlying computational tools, experimental techniques, and theoretical models were developed independently. The primary result presented here is a framework for the unification of informatics- and physics-based data associated with DNA, nucleosomes, and chromatin. The framework is based on the mathematical representation of geometrically exact rods and the generalization of DNA basepair step parameters. Data unification enables researchers to integrate computational, experimental, and theoretical approaches for the study of chromatin biology. The framework can be implemented using model-view-controller design principles, existing genome browsers, and existing molecular visualization tools. We developed a minimal, web-based genome dashboard, G-Dash-min, and applied it to two simple examples to demonstrate the usefulness of data unification and proof of concept. Genome dashboards developed using the framework and design principles presented here are extensible and customizable and are therefore more broadly applicable than the examples presented. We expect a number of purpose-specific genome dashboards to emerge as a novel means of investigating structure-function relationships for genomes that range from basepairs to entire chromosomes and for generating, validating, and testing mechanistic hypotheses.

Modeling, molecular dynamics, and docking assessment of transcription factor rho: a potential drug target in Brucella melitensis 16M

The zoonotic disease brucellosis, a chronic condition in humans affecting renal and cardiac systems and causing osteoarthritis, is caused by Brucella, a genus of Gram-negative, facultative, intracellular pathogens. The mode of transmission and the virulence of the pathogens are still enigmatic. Transcription regulatory elements, such as rho proteins, play an important role in the termination of transcription and/or the selection of genes in Brucella. Adverse effects of the transcription inhibitors play a key role in the non-successive transcription challenges faced by the pathogens. In the investigation presented here, we computationally predicted the transcription termination factor rho (TtFRho) inhibitors against Brucella melitensis 16M via a structure-based method. In view the unknown nature of its crystal structure, we constructed a robust three-dimensional homology model of TtFRho’s structure by comparative modeling with the crystal structure of the Escherichia coli TtFRho (Protein Data Bank ID: 1PVO) as a template in MODELLER (v 9.10). The modeled structure was optimized by applying a molecular dynamics simulation for 2 ns with the CHARMM (Chemistry at HARvard Macromolecular Mechanics) 27 force field in NAMD (NAnoscale Molecular Dynamics program; v 2.9) and then evaluated by calculating the stereochemical quality of the protein. The flexible docking for the interaction phenomenon of the template consists of ligand-related inhibitor molecules from the ZINC (ZINC Is Not Commercial) database using a structure-based virtual screening strategy against minimized TtFRho. Docking simulations revealed two inhibitors compounds – ZINC24934545 and ZINC72319544 – that showed high binding affinity among 2,829 drug analogs that bind with key active-site residues; these residues are considered for protein-ligand binding and unbinding pathways via steered molecular dynamics simulations. Arg215 in the model plays an important role in the stability of the protein-ligand complex via a hydrogen bonding interaction by aromatic-π contacts, and the ADMET (absorption, distribution, metabolism, and excretion) analysis of best leads indicate nontoxic in nature with good potential for drug development.

A zinc-finger like metal binding site in the nucleosome

UV spectroscopy demonstrated that chicken mononucleosomes bind Co(II) and Zn(II) ions at submicromolar concentrations in a tetrahedral mode, at a conserved zinc finger-like site, composed of Cys110 and His113 residues of both H3 molecules. Neither of these metal ions substituted for another, indicating a limited binding reversibility. Molecular modeling indicated that the tetrahedral site is formed by unhindered rotations around Calpha-Cbeta bonds in the side chains of the zinc binding residues. The resulting local rearrangement of the protein structure shields the bound metal ion from the solvent, explaining the observed lack of reversibility of the binding. Consequences of these findings for zinc homeostasis, metal toxicology and nucleosomal regulation are discussed.

Thyroid hormone receptor DNA binding is required for both positive and negative gene regulation

The beta isoform of thyroid hormone receptor (TR-beta) has a key role in the feedback regulation of the hypothalamic-pituitary-thyroid (H-P-T) axis. The mechanism of trans-repression of the hypothalamic thyrotropin-releasing hormone (TRH) and pituitary thyroid-stimulating hormone (TSH) subunit genes, however, remains poorly understood. A number of distinct mechanisms for TR-beta-mediated negative regulation by thyroid hormone have been proposed, including those that require and do not require DNA binding. To clarify the importance of DNA binding in negative regulation, we constructed a DNA-binding mutant of TR-beta in which two amino acids within the P box were altered (GSG for EGG) to resemble that found in the glucocorticoid receptor (GR). We termed this mutant GS125, and as expected, it displayed low binding affinities for positive and negative thyroid hormone-response element (pTRE and nTRE, respectively) in gel-mobility shift assays. In transient transfection assays, the GS125 mutant abolished transactivation on three classic pTREs (DR+4, LAP, and PAL) and all negatively regulated promoters in the H-P-T axis (TRH, TSH-beta, and TSH-alpha). However, GS125 TR-beta bound to a composite TR/GR-response element and was fully functional on this hybrid TR/GR-response element. Moreover, the GS125 TR-beta mutant displayed normal interactions with transcriptional cofactors in mammalian two-hybrid assays. These data do not support a DNA-binding independent mechanism for thyroid hormone negative regulation in the H-P-T axis.

The investigation of the effects of counterions in protein dynamics simulations

Molecular simulations able to exactly represent solvated charged proteins are helpful in understanding protein dynamics, structure and function. In the present study we have used two different starting structures of papain (a typical, stable, globular protein of intermediate net charge) and different modeling procedures to evaluate some effects of counterions in simulations. A number of configurations have been generated and relaxed for each system by various combinations of constrained simulated annealing and molecular dynamics procedures, using the AMBER force field. The analysis of trajectories shows that the simulations of solvated proteins are moderately sensitive to the presence of counterions. However, this sensitivity is highly dependent on the starting model and different procedures of equilibration used. The neutralized systems tend to evince smaller root mean square deviations regardless of the system investigated and the simulation procedure used. The results of parameterized fitting of the simulated structures to the crystallographic data, giving quantitative measure of the total charge influence on the stability of various elements of the secondary structure, revealed a clear scatter of different reactions of various systems' secondary structures to counterions addition: some systems apparently were stabilized when neutralized, while the others were not. Thus, one cannot unequivocally state, despite consideration of specific simulation conditions, whether protein secondary structures are more stable when they have neutralized charges. This suggests that caution should be taken when claiming the stabilizing effect of counterions in simulations other than those involving small, unstable polypeptides or highly charged proteins.

Structure, interaction, dynamics and solvent effects on the DNA-EcoRI complex in aqueous solution from molecular dynamics simulation

A 0.7-ns molecular dynamics simulation of the DNA-EcoRI complex in a 7.0-A solvent shell indicated a stable behavior of the system. No significant evaporation or smearing of the solvent's outer boundary occurred. The structure and the intermolecular interactions were found to be well maintained during the simulation. The interaction pattern in the simulation was found to be very similar to that in the crystal structure. Most of the specific interactions between the DNA and the protein were found to be enhanced in the simulation compared to that in the crystal structure as a result of improved interaction geometry. The nonspecific interactions were found to be stronger than the specific ones. The specific interactions between the N7 atoms of Gua(4) or Ade(5) or Ade(6) and the protein were found to be present over almost the entire time of the simulation, whereas hydrogen bonds involving the amino groups of the Ade(5) and Ade(6) with the protein were found to be relatively weaker, with lower probability and shorter lifetime. The time evolution of the root mean square deviations of the DNA and the protein were highly correlated even at the later part of the simulation, showing the tight binding between them. Several long-lived water bridges were found between the DNA backbone atoms and the protein and also between the two protein monomers, which increased the overall stability of the complex. The two protein monomers were found to interact strongly with each other. The energy of the DNA kink deformation was estimated as approximately 31 kcal/mol.

Steered molecular dynamics simulation of the Rieske subunit motion in the cytochrome bc(1) complex

Crystallographic structures of the mitochondrial ubiquinol/cytochrome c oxidoreductase (cytochrome bc(1) complex) suggest that the mechanism of quinol oxidation by the bc(1) complex involves a substantial movement of the soluble head of the Rieske iron-sulfur protein (ISP) between reaction domains in cytochrome b and cytochrome c(1) subunits. In this paper we report the results of steered molecular dynamics simulations inducing, through an applied torque within 1 ns, a 56 degrees rotation of the soluble domain of ISP. For this purpose, a solvated structure of the bc(1) complex in a phospholipid bilayer (a total of 206,720 atoms) was constructed. A subset of 91,061 atoms was actually simulated with 45,131 moving atoms. Point charge distributions for the force field parametrization of heme groups and the Fe(2)S(2) cluster of the Rieske protein included in the simulated complex were determined. The simulations showed that rotation of the soluble domain of ISP is actually feasible. Several metastable conformations of the ISP during its rotation were identified and the interactions stabilizing the initial, final, and intermediate positions of the soluble head of the ISP domain were characterized. A pathway for proton conduction from the Q(o) site to the solvent via a water channel has been identified.

Influence of nucleosome structure on the three-dimensional folding of idealized minichromosomes

BACKGROUND

The closed circular, multinucleosome-bound DNA comprising a minichromosome provides one of the best known examples of chromatin organization beyond the wrapping of the double helix around the core of histone proteins. This higher level of chain folding is governed by the topology of the constituent nucleosomes and the spatial disposition of the intervening protein-free DNA linkers.

RESULTS

By simplifying the protein-DNA assembly to an alternating sequence of virtual bonds, the organization of a string of nucleosomes on the minichromosome can be treated by analogy to conventional chemical depictions of macromolecular folding in terms of the bond lengths, valence angles, and torsions of the chain. If the nucleosomes are evenly spaced and the linkers are sufficiently short, regular minichromosome structures can be identified from analytical expressions that relate the lengths and angles formed by the virtual bonds spanning the nucleosome-linker repeating units to the pitch and radius of the organized quaternary structures that they produce.

CONCLUSIONS

The resulting models with 4-24 bound nucleosomes illustrate how a minichromosome can adopt the low-writhe folding motifs deduced from biochemical studies, and account for published images of the 30 nm chromatin fiber and the simian virus 40 (SV40) nucleohistone core. The marked sensitivity of global folding to the degree of protein-DNA interactions and the assumed nucleosomal shape suggest potential mechanisms for chromosome rearrangements upon histone modification.

Unbinding of retinoic acid from its receptor studied by steered molecular dynamics

Retinoic acid receptor (RAR) is a ligand-dependent transcription factor that regulates the expression of genes involved in cell growth, differentiation, and development. Binding of the retinoic acid hormone to RAR is accompanied by conformational changes in the protein which induce transactivation or transrepression of the target genes. In this paper we present a study of the hormone binding/unbinding process in order to clarify the role of some of the amino acid contacts and identify possible pathways of the all-trans retinoic acid binding/unbinding to/from human retinoic acid receptor (hRAR)-gamma. Three possible pathways were explored using steered molecular dynamics simulations. Unbinding was induced on a time scale of 1 ns by applying external forces to the hormone. The simulations suggest that the hormone may employ one pathway for binding and an alternative "back door" pathway for unbinding.

Structure prediction of a complex between the chromosomal protein HMG-D and DNA

Non-histone chromosomal proteins are an important part of nuclear structure and function due to their ability to interact with DNA to form and modulate chromatin structure and regulate gene expression. However, the understanding of the function of chromosomal proteins at the molecular level has been hampered by the lack of structures of chromosomal protein-DNA complexes. We have carried out a molecular dynamics modeling study to provide insight into the mode of DNA binding to the chromosomal HMG-domain protein, HMG-D. Three models of a complex of HMG-D bound to DNA were derived through docking the protein to two different DNA fragments of known structure. Molecular dynamics simulations of the complexes provided data indicating the most favorable model. This model was further refined by molecular dynamics simulation and extensively analyzed. The structure of the corresponding HMG-D-DNA complex exhibits many features seen in the NMR structures of the sequence-specific HMG-domain-DNA complexes, lymphoid enhancer factor 1 (LEF-1) and testis determining factor (SRY). The model reveals differences from these known structures that suggest how chromosomal proteins bind to many different DNA sequences with comparable affinity.

Binding of the estrogen receptor to DNA. The role of waters

Molecular dynamics simulations are carried out to investigate the binding of the estrogen receptor, a member of the nuclear hormone receptor family, to specific and non-specific DNA. Two systems have been simulated, each based on the crystallographic structure of a complex of a dimer of the estrogen receptor DNA binding domain with DNA. One structure includes the dimer and a consensus segment of DNA, ds(CCAGGTCACAGTGACCTGG); the other structure includes the dimer and a nonconsensus segment of DNA, ds(CCAGAACACAGTGACCTGG). The simulations involve an atomic model of the protein-DNA complex, counterions, and a sphere of explicit water with a radius of 45 A. The molecular dynamics package NAMD was used to obtain 100 ps of dynamics for each system with complete long-range electrostatic interactions. Analysis of the simulations revealed differences in the protein-DNA interactions for consensus and nonconsensus sequences, a bending and unwinding of the DNA, a slight rearrangement of several amino acid side chains, and inclusion of water molecules at the protein-DNA interface region. Our results indicate that binding specificity and stability is conferred by a network of direct and water mediated protein-DNA hydrogen bonds. For the consensus sequence, the network involves three water molecules, residues Glu-25, Lys-28, Lys-32, Arg-33, and bases of the DNA. The binding differs for the nonconsensus DNA sequence in which case the fluctuating network of hydrogen bonds allows water molecules to enter the protein-DNA interface. We conclude that water plays a role in furnishing DNA binding specificity to nuclear hormone receptors.