Three-dimensional modeling of and ligand docking to vitamin D receptor ligand binding domain

The ligand binding domain of the human vitamin D receptor (VDR) was modeled based on the crystal structure of the retinoic acid receptor. The ligand binding pocket of our VDR model is spacious at the helix 11 site and confined at the β-turn site. The ligand 1α,25-dihydroxyvitamin D3 was assumed to be anchored in the ligand binding pocket with its side chain heading to helix 11 (site 2) and the A-ring toward the β-turn (site 1). Three residues forming hydrogen bonds with the functionally important 1α- and 25-hydroxyl groups of 1α,25-dihydroxyvitamin D3 were identified and confirmed by mutational analysis: the 1α-hydroxyl group is forming pincer-type hydrogen bonds with S237 and R274 and the 25-hydroxyl group is interacting with H397. Docking potential for various ligands to the VDR model was examined, and the results are in good agreement with our previous three-dimensional structure-function theory.

Three-dimensional structure of poliovirus receptor bound to poliovirus

Poliovirus initiates infection by binding to its cellular receptor (Pvr). We have studied this interaction by using cryoelectron microscopy to determine the structure, at 21-Å resolution, of poliovirus complexed with a soluble form of its receptor (sPvr). This density map aided construction of a homology-based model of sPvr and, in conjunction with the known crystal structure of the virus, allowed delineation of the binding site. The virion does not change significantly in structure on binding sPvr in short incubations at 4°C. We infer that the binding configuration visualized represents the initial interaction that is followed by structural changes in the virion as infection proceeds. sPvr is segmented into three well-defined Ig-like domains. The two domains closest to the virion (domains 1 and 2) are aligned and rigidly connected, whereas domain 3 diverges at an angle of ≈60°. Two nodules of density on domain 2 are identified as glycosylation sites. Domain 1 penetrates the “canyon” that surrounds the 5-fold protrusion on the capsid surface, and its binding site involves all three major capsid proteins. The inferred pattern of virus–sPvr interactions accounts for most mutations that affect the binding of Pvr to poliovirus.

The two-dimensional IR nonlinear spectroscopy of a cyclic penta-peptide in relation to its three-dimensional structure

A form of two-dimensional (2D) vibrational spectroscopy, which uses two ultrafast IR laser pulses, is used to examine the structure of a cyclic penta-peptide in solution. Spectrally resolved cross peaks occur in the off-diagonal region of the 2D IR spectrum of the amide I region, analogous to those in 2D NMR spectroscopy. These cross peaks measure the coupling between the different amide groups in the structure. Their intensities and polarizations relate directly to the three-dimensional structure of the peptide. With the help of a model coupling Hamiltonian, supplemented by density functional calculations, the spectra of this penta-peptide can be regenerated from the known solution phase structure. This 2D-IR measurement, with an intrinsic time resolution of less than 1 ps, could be used in all time regimes of interest in biology.

Photoactive yellow protein: A structural prototype for the three-dimensional fold of the PAS domain superfamily

PAS domains are found in diverse proteins throughout all three kingdoms of life, where they apparently function in sensing and signal transduction. Although a wealth of useful sequence and functional information has become recently available, these data have not been integrated into a three-dimensional (3D) framework. The very early evolutionary development and diverse functions of PAS domains have made sequence analysis and modeling of this protein superfamily challenging. Limited sequence similarities between the ∼50-residue PAS repeats and one region of the bacterial blue-light photosensor photoactive yellow protein (PYP), for which ground-state and light-activated crystallographic structures have been determined to high resolution, originally were identified in sequence searches using consensus sequence probes from PAS-containing proteins. Here, we found that by changing a few residues particular to PYP function, the modified PYP sequence probe also could select PAS protein sequences. By mapping a typical ∼150-residue PAS domain sequence onto the entire crystallographic structure of PYP, we show that the PAS sequence similarities and differences are consistent with a shared 3D fold (the PAS/PYP module) with obvious potential for a ligand-binding cavity. Thus, PYP appears to prototypically exhibit all the major structural and functional features characteristic of the PAS domain superfamily: the shared PAS/PYP modular domain fold of ∼125–150 residues, a sensor function often linked to ligand or cofactor (chromophore) binding, and signal transduction capability governed by heterodimeric assembly (to the downstream partner of PYP). This 3D PAS/PYP module provides a structural model to guide experimental testing of hypotheses regarding ligand-binding, dimerization, and signal transduction.

The Three-Dimensional Structure of Aspergillus niger Pectin Lyase B at 1.7-Å Resolution1

The three-dimensional structure of Aspergillus niger pectin lyase B (PLB) has been determined by crystallographic techniques at a resolution of 1.7 Å. The model, with all 359 amino acids and 339 water molecules, refines to a final crystallographic R factor of 16.5%. The polypeptide backbone folds into a large right-handed cylinder, termed a parallel β helix. Loops of various sizes and conformations protrude from the central helix and probably confer function. The largest loop of 53 residues folds into a small domain consisting of three antiparallel β strands, one turn of an α helix, and one turn of a 310 helix. By comparison with the structure of Erwinia chrysanthemi pectate lyase C (PelC), the primary sequence alignment between the pectate and pectin lyase subfamilies has been corrected and the active site region for the pectin lyases deduced. The substrate-binding site in PLB is considerably less hydrophilic than the comparable PelC region and consists of an extensive network of highly conserved Trp and His residues. The PLB structure provides an atomic explanation for the lack of a catalytic requirement for Ca2+ in the pectin lyase family, in contrast to that found in the pectate lyase enzymes. Surprisingly, however, the PLB site analogous to the Ca2+ site in PelC is filled with a positive charge provided by a conserved Arg in the pectin lyases. The significance of the finding with regard to the enzymatic mechanism is discussed.

Faithful representation of similarities among three-dimensional shapes in human vision.

Efficient and reliable classification of visual stimuli requires that their representations reside a low-dimensional and, therefore, computationally manageable feature space. We investigated the ability of the human visual system to derive such representations from the sensory input-a highly nontrivial task, given the million or so dimensions of the visual signal at its entry point to the cortex. In a series of experiments, subjects were presented with sets of parametrically defined shapes; the points in the common high-dimensional parameter space corresponding to the individual shapes formed regular planar (two-dimensional) patterns such as a triangle, a square, etc. We then used multidimensional scaling to arrange the shapes in planar configurations, dictated by their experimentally determined perceived similarities. The resulting configurations closely resembled the original arrangements of the stimuli in the parameter space. This achievement of the human visual system was replicated by a computational model derived from a theory of object representation in the brain, according to which similarities between objects, and not the geometry of each object, need to be faithfully represented.

Three-dimensional structure of human protein kinase C interacting protein 1, a member of the HIT family of proteins.

The three-dimensional structure of protein kinase C interacting protein 1 (PKCI-1) has been solved to high resolution by x-ray crystallography using single isomorphous replacement with anomalous scattering. The gene encoding human PKCI-1 was cloned from a cDNA library by using a partial sequence obtained from interactions identified in the yeast two-hybrid system between PKCI-1 and the regulatory domain of protein kinase C-beta. The PKCI-1 protein was expressed in Pichia pastoris as a dimer of two 13.7-kDa polypeptides. PKCI-1 is a member of the HIT family of proteins, shown by sequence identity to be conserved in a broad range of organisms including mycoplasma, plants, and humans. Despite the ubiquity of this protein sequence in nature, no distinct function has been shown for the protein product in vitro or in vivo. The PKCI-1 protomer has an alpha+beta meander fold containing a five-stranded antiparallel sheet and two helices. Two protomers come together to form a 10-stranded antiparallel sheet with extensive contacts between a helix and carboxy terminal amino acids of a protomer with the corresponding amino acids in the other protomer. PKCI-1 has been shown to interact specifically with zinc. The three-dimensional structure has been solved in the presence and absence of zinc and in two crystal forms. The structure of human PKCI-1 provides a model of this family of proteins which suggests a stable fold conserved throughout nature.

Three-dimensional tracking of motile bacteria near a solid planar surface.

Knowing how motile bacteria move near and along a solid surface is crucial to understanding such diverse phenomena as the migration of infectious bacteria along a catheter, biofilm growth, and the movement of bacteria through the pore spaces of saturated soil, a critical step in the in situ bioremediation of contaminated aquifers. In this study, a tracking microscope is used to record the three-dimensional motion of Escherichia coli near a planar glass surface. Data from the tracking microscope are analyzed to quantify the effects of bacteria-surface interactions on the swimming behavior of bacteria. The speed of cells approaching the surface is found to decrease in agreement with the mathematical model of Ramia et al. [Ramia, M., Tullock, D. L. & Phan-Tien, N. (1993) Biophys J. 65,755-778], which represents the bacteria as spheres with a single polar flagellum rotating at a constant rate. The tendency of cells to swim adjacent to the surface is shown in computer-generated reproductions of cell traces. The attractive interaction potential between the cells and the solid surface is offered as one of several possible explanations for this tendency.

Three-dimensional scroll waves of cAMP could direct cell movement and gene expression in Dictyostelium slugs.

Complex three-dimensional waves of excitation can explain the observed cell movement pattern in Dictyostelium slugs. Here we show that these three-dimensional waves can be produced by a realistic model for the cAMP relay system [Martiel, J. L. & Goldbeter, A. (1987) Biophys J. 52, 807-828]. The conversion of scroll waves in the prestalk zone of the slug into planar wave fronts in the prespore zone can result from a smaller fraction of relaying cells in the prespore zone. Further, we show that the cAMP concentrations to which cells in a slug are exposed over time display a simple pattern, despite the complex spatial geometry of the waves. This cAMP distribution agrees well with observed patterns of cAMP-regulated cell type-specific gene expression. The core of the spiral, which is a region of low cAMP concentration, might direct expression of stalk-specific genes during culmination.

Three-dimensional inverse modelling of magnetic anomaly sources based on a genetic algorithm.

We present a modelling method to estimate the 3-D geometry and location of homogeneously magnetized sources from magnetic anomaly data. As input information, the procedure needs the parameters defining the magnetization vector (intensity, inclination and declination) and the Earth's magnetic field direction. When these two vectors are expected to be different in direction, we propose to estimate the magnetization direction from the magnetic map. Then, using this information, we apply an inversion approach based on a genetic algorithm which finds the geometry of the sources by seeking the optimum solution from an initial population of models in successive iterations through an evolutionary process. The evolution consists of three genetic operators (selection, crossover and mutation), which act on each generation, and a smoothing operator, which looks for the best fit to the observed data and a solution consisting of plausible compact sources. The method allows the use of non-gridded, non-planar and inaccurate anomaly data and non-regular subsurface partitions. In addition, neither constraints for the depth to the top of the sources nor an initial model are necessary, although previous models can be incorporated into the process. We show the results of a test using two complex synthetic anomalies to demonstrate the efficiency of our inversion method. The application to real data is illustrated with aeromagnetic data of the volcanic island of Gran Canaria (Canary Islands).

Critical behavior of three-dimensional disordered Potts models with many states

We study the 3D Disordered Potts Model with p = 5 and p = 6. Our numerical simulations (that severely slow down for increasing p) detect a very clear spin glass phase transition. We evaluate the critical exponents and the critical value of the temperature, and we use known results at lower p values to discuss how they evolve for increasing p. We do not find any sign of the presence of a transition to a ferromagnetic regime.

First-order transition in a three-dimensional disordered system

We present the first detailed numerical study in three dimensions of a first-order phase transition that remains first order in the presence of quenched disorder (specifically, the ferromagnetic-paramagnetic transition of the site-diluted four states Potts model). A tricritical point, which lies surprisingly near the pure-system limit and is studied by means of finite-size scaling, separates the first-order and second-order parts of the critical line. This investigation has been made possible by a new definition of the disorder average that avoids the diverging-variance probability distributions that plague the standard approach. Entropy, rather than free energy, is the basic object in this approach that exploits a recently introduced microcanonical Monte Carlo method.

A Three-Dimensional Stereotaxic MRI Brain Atlas of the Cichlid Fish Oreochromis mossambicus

The African cichlid Oreochromis mossambicus (Mozambique tilapia) has been used as a model system in a wide range of behavioural and neurobiological studies. The increasing number of genetic tools available for this species, together with the emerging interest in its use for neurobiological studies, increased the need for an accurate hodological mapping of the tilapia brain to supplement the available histological data. The goal of our study was to elaborate a three-dimensional, high-resolution digital atlas using magnetic resonance imaging, supported by Nissl staining. Resulting images were viewed and analysed in all orientations (transverse, sagittal, and horizontal) and manually labelled to reveal structures in the olfactory bulb, telencephalon, diencephalon, optic tectum, and cerebellum. This high resolution tilapia brain atlas is expected to become a very useful tool for neuroscientists using this fish model and will certainly expand their use in future studies regarding the central nervous system.

Characterising the three-dimensional ozone distribution of a tidally locked Earth-like planet

We simulate the 3D ozone distribution of a tidally locked Earth-like exoplanet using the high-resolution, 3D chemistry climate model CESM1(WACCM) and study how the ozone layer of a tidally locked Earth (TLE) (ΩTLE = 1/365 days) differs from that of our present-day Earth (PDE) (ΩPDE = 1/1 day). The middle atmosphere reaches a steady state a symptotically within the first 80 days of the simulation. An upwelling, centred on the subsolar point, is present on the day side while a downwelling, centred on the antisolar point, is present on the night side. In the mesosphere, we find similar global ozone distributions for the TLE and the PDE, with decreased ozone on the day side and enhanced ozone on the night side. In the lower mesosphere, a jet stream transitions into a large-scale vortex around a low-pressure system, located at low latitudes of the TLE night side. In the middle stratosphere, the concentration of odd oxygen is approximately equal to that of the ozone [(Ox) ≈ (O3)]. At these altitudes, the lifetime of odd oxygen is ~16 h and the transport processes significantly contribute to the global distribution of stratospheric ozone. Compared to the PDE, where the strong Coriolis force acts as a mixing barrier between low and high latitudes, the transport processes of the TLE are governed by jet streams variable in the zonal and meridional directions. In the middle stratosphere of the TLE, we find high ozone values on the day side, due to the increased production of atomic oxygen on the day side, where it immediately recombines with molecular oxygen to form ozone. In contrast, the ozone is depleted on the night side, due to changes in the solar radiation distribution and the presence of a downwelling. As a result of the reduced Coriolis force, the tropical and extratropical air masses are well mixed and the global temperature distribution of the TLE stratosphere has smaller horizontal gradients than the PDE. Compared to the PDE, the total ozone column global mean is reduced by ~19.3 %. The day side and the night side total ozone column means are reduced by 23.21 and 15.52 %, respectively. Finally, we present the total ozone column (TOC) maps as viewed by a remote observer for four phases of the TLE during its revolution around the star. The mean TOC values of the four phases of the TLE vary by up to 23 %.

Methods of visualization and analysis of cardiac depolarization in the three dimensional space

The master thesis presents methods for intellectual analysis and visualization 3D EKG in order to increase the efficiency of ECG analysis by extracting additional data. Visualization is presented as part of the signal analysis tasks considered imaging techniques and their mathematical description. Have been developed algorithms for calculating and visualizing the signal attributes are described using mathematical methods and tools for mining signal. The model of patterns searching for comparison purposes of accuracy of methods was constructed, problems of a clustering and classification of data are solved, the program of visualization of data is also developed. This approach gives the largest accuracy in a task of the intellectual analysis that is confirmed in this work. Considered visualization and analysis techniques are also applicable to the multi-dimensional signals of a different kind.