Iterative Model Refinement Research Articles

Structure and Mechanism of Iron Translocation by a Dps Protein from Microbacterium arborescens

Dps (DNA protection during starvation) enzymes are a major class of dodecameric proteins that bacteria use to detoxify their cytosol through the uptake of reactive iron species. In the stationary growth phase of bacteria, Dps enzymes are primarily used to protect DNA by biocrystallization. To characterize the wild type Dps protein from Microbacterium arborescens that displays additional catalytic functions (amide hydrolysis and synthesis), we determined the crystal structure to a resolution of 2.05 Å at low iron content. The structure shows a single iron at the ferroxidase center coordinated by an oxo atom, one water molecule, and three ligating residues. An iron-enriched protein structure was obtained at 2 Å and shows the stepwise uptake of two hexahydrated iron atoms moving along channels at the 3-fold axis before a restriction site inside the channels requires removal of the hydration sphere. Supporting biochemical data provide insight into the regulation of this acylamino acid hydrolase. Moreover, the peroxidase activity of the protein was determined. The influence of iron and siderophores on the expression of acylamino acid hydrolase was monitored during several stages of cell growth. Altogether our data provide an interesting view of an unusual Dps-like enzyme evolutionarily located apart from the large Dps sequence clusters.

A workflow to facilitate three-dimensional geometrical modelling of complex poly-deformed geological units

The three-dimensional (3D) geometry of complexly deformed regions is often beyond the scope of simple 2D and 2.5D representation in cross-sections and block diagrams. Methods must be developed for fully three-dimensional representation. A workflow for such three-dimensional modelling of structurally complex areas is presented. Data requirements are tailored to typical results of structural field work in strongly deformed rocks from mid-crustal levels. The workflow is based on data evaluation and data export using standard geographical information system and database management systems. Three-dimensional modelling in polyphase deformed areas is highly dependent on the correct interpretation of variations in the orientation of the dominant planar fabrics. The computer-aided earth-modelling software 3D GeoModeller is well adapted for managing this problem. It calculates the geometry of geological interfaces taking into account simultaneously the foliation data and location of lithological contacts. The workflow is based on iterative model refinement based on interactive editing of the geometry in section and map views with data assessment and pre-processing using GIS software. This approach assures internal consistency of the resulting three-dimensional models. Modelling of the Lower Lepontine Nappes in the Central Alps is used as an illustrative example from a complexly deformed terrain.

Image-Guided Modeling of Virus Growth and Spread

Although many tools of cellular and molecular biology have been used to characterize single intracellular cycles of virus growth, few culture methods exist to study the dynamics of spatially spreading viruses over multiple generations. We have previously developed a method that addresses this need by tracking the spread of focal infections using immunocytochemical labeling and digital imaging. Here, we build reaction-diffusion models to account for spatio-temporal patterns formed by the spreading viral infection front as well as data from a single cycle of virus growth (one-step growth). Systems with and without the interferon-mediated antiviral response of the host cells are considered. Dynamic images of the spreading infections guide iterative model refinement steps that lead to reproduction of all of the salient features contained in the images, not just the velocity of the infection front. The optimal fits provide estimates for key parameters such as virus-host binding and the production rate of interferon. For the examined data, highly-lumped infection models that ignore the one-step growth dynamics provide a comparable fit to models that more accurately account for these dynamics, highlighting the fact that increased model complexity does not necessarily translate to improved fit. This work demonstrates how model building can facilitate the interpretation of experiments by highlighting contributions from both biological and methodological factors.

Iterative-build OMIT maps: map improvement by iterative model building and refinement without model bias

A procedure for carrying out iterative model building, density modification and refinement is presented in which the density in an OMIT region is essentially unbiased by an atomic model. Density from a set of overlapping OMIT regions can be combined to create a composite 'iterative-build' OMIT map that is everywhere unbiased by an atomic model but also everywhere benefiting from the model-based information present elsewhere in the unit cell. The procedure may have applications in the validation of specific features in atomic models as well as in overall model validation. The procedure is demonstrated with a molecular-replacement structure and with an experimentally phased structure and a variation on the method is demonstrated by removing model bias from a structure from the Protein Data Bank.

Structure of the Human Lung Cytochrome P450 2A13

The human lung cytochrome P450 2A13 (CYP2A13) activates the nicotine-derived procarcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) into DNA-altering compounds that cause lung cancer. Another cytochrome P450, CYP2A6, is also present in human lung, but at much lower levels. Although these two enzymes are 93.5% identical, CYP2A13 metabolizes NNK with much lower K(m) values than does CYP2A6. To investigate the structural differences between these two enzymes the structure of CYP2A13 was determined to 2.35A by x-ray crystallography and compared with structures of CYP2A6. As expected, the overall CYP2A13 and CYP2A6 structures are very similar with an average root mean square deviation of 0.5A for the Calpha atoms. Like CYP2A6, the CYP2A13 active site cavity is small and highly hydrophobic with a cluster of Phe residues composing the active site roof. Active site residue Asn(297) is positioned to hydrogen bond with an adventitious ligand, identified as indole. Amino acid differences between CYP2A6 and CYP2A13 at positions 117, 300, 301, and 208 relate to different orientations of the ligand plane in the two protein structures and may underlie the significant variations observed in binding and catalysis of many CYP2A ligands. In addition, docking studies suggest that residues 365 and 366 may also contribute to differences in NNK metabolism.

Interpretation of ensembles created by multiple iterative rebuilding of macromolecular models.

Automation of iterative model building, density modification and refinement in macromolecular crystallography has made it feasible to carry out this entire process multiple times. By using different random seeds in the process, a number of different models compatible with experimental data can be created. Sets of models were generated in this way using real data for ten protein structures from the Protein Data Bank and using synthetic data generated at various resolutions. Most of the heterogeneity among models produced in this way is in the side chains and loops on the protein surface. Possible interpretations of the variation among models created by repetitive rebuilding were investigated. Synthetic data were created in which a crystal structure was modelled as the average of a set of ;perfect' structures and the range of models obtained by rebuilding a single starting model was examined. The standard deviations of coordinates in models obtained by repetitive rebuilding at high resolution are small, while those obtained for the same synthetic crystal structure at low resolution are large, so that the diversity within a group of models cannot generally be a quantitative reflection of the actual structures in a crystal. Instead, the group of structures obtained by repetitive rebuilding reflects the precision of the models, and the standard deviation of coordinates of these structures is a lower bound estimate of the uncertainty in coordinates of the individual models.

REFINEMENT: A search framework for the identification of interferon-responsive elements in DNA sequences – a case study with ISRE and GAS

Interferons (IFN) are a family of pleiotropic secreted proteins that play a key role in mediating antiviral and apoptotic responses, and in immune modulation. Interferons induce a large number of genes through activating the janus tyrosine kinase (JAK)–signal transducers and activators of transcription proteins (STAT) pathway, and the binding of transcription factors to upstream regions of the inducible genes (interferon-stimulated gene, ISG) at specific DNA regulatory elements known as interferon-stimulated response element (ISRE) and gamma-activated sequence (GAS). We have previously performed DNA micro-arrays on peripheral blood mononuclear cells (PBMC) treated with interferon-α in culture and showed that approximately 700 genes are significantly modulated ( P ≤ 0.001). In order to search for ISRE and GAS we have developed a framework called regulatory element finding with iteration and effective model refinement (REFINEMENT) using an existing program (HMMER) and a standard discriminating scoring technique. Although REFINEMENT uses existing programs, our framework itself is novel as it effectively discriminates occurrences using an iterative model refinement technique. REFINEMENT has detected either ISRE or GAS sequence in all of the genes shown to be induced at a P-value ≤ 0.001. There were far more functional occurrences in ISRE than in GAS, suggesting that ISRE plays a greater role in response to interferon-α than GAS sequences. This method can be used to identify such sequences in any set of genes. REFINEMENT is non-commercial and is accessible at http://cancer.informatics.indiana.edu/ttsukaha/miltonlab/refinement/.