Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

Development and evaluation of a video-based modeling strategy for the dissemination phase of a school health tobacco -use prevention program

Even the best school health education programs will be unsuccessful if they are not disseminated effectively in a manner that encourages classroom adoption and implementation. This study involved two components: (1) the development of a videotape intervention to be used in the dissemination phase of a 4-year, NCI-funded diffusion study and (2) the evaluation of that videotape intervention strategy in comparison with a print (information transfer) strategy. Conceptualization has been guided by Social Learning Theory, Diffusion Theory, and communication theory. Additionally, the PRECEDE Framework has been used. Seventh and 8th grade classroom teachers from Spring Branch Independent School District in west Houston participated in the evaluation of the videotape and print interventions using a 57-item preadoption survey instrument developed by the UT Center for Health Promotion Research and Development. Two-way ANOVA was used to study individual score differences for five outcome variables: Total Scale Score (comprised of 57 predisposing, enabling, and reinforcing items), Adoption Characteristics Subscale, Attitude Toward Innovation Subscale, Receptivity Toward Innovation, and Reinforcement Subscale. The aim of the study is to compare the effect upon score differences of video and print interventions alone and in combination. Seventy-three 7th and 8th grade classroom teachers completed the study providing baseline and post-intervention measures on factors related to the adoption and implementation of tobacco-use prevention programs. Two-way ANOVA, in relation to the study questions, found significant scoring differences for those exposed to the videotape intervention alone for both the Attitude Toward Innovation Subscale and the Receptivity to Adopt Subscale. No significant results were found to suggest that print alone influences favorable scoring differences between baseline and post-intervention testing. One interaction effect was found suggesting video and print combined are more effective for influencing favorable scoring differences for the Reinforcement for the Adoption Subscale.^ This research is unique in that it represents a newly emerging field in health promotion communications research with implications for Social Learning Theory, Diffusion Theory, and communication science that are applicable to the development of improved school health interventions. ^

Strontium isotopes in pore waters of east equatorial Pacific sediments

Pore-water samples from the equatorial sedimentary bulge area show reversals in depth profiles of 87Sr/86Sr ratios at the sediment/basement interface. Results of this work support inferences made from previous pore-water data (from DSDP drilling in the area) that large-scale horizontal advection of seawater has occurred through the basement underlying the thick sedimentary sequence in this region. The area of apparent advection includes the eastern part of the equatorial high-productivity zone and part of the Guatemala Basin. We attempted to find links between the observed near-basement reversals in pore-water chemistry and sedimentary thickness, age, and topography of the area. Most of the sites that show horizontal advection have disturbed basement topography or outcrops within 10 to 20 km, suggesting that the cooling effects of outcrops may extend for at least 20 km horizontally. Heat-flow data from the area were compared to determine whether sites showing near-bottom chemistry reversals were consistent with areas of low conductive heat flow. This was generally true for the area of the sedimentary bulge and Guatemala Basin. Not enough pore-water data from the Nazca Plate were available to establish any reliable systematics. Because the high-productivity area is well-sealed from hydrothermal circulation, the missing heat must be lost by horizontal advective heat transport. From profiles of strontium isotopes and other elements that show departure from seawater values with increasing depth in the sediments, but return to seawater values near the basement, it appears that water flows relatively freely through much of the oceanic crust, even when sealed by considerable sedimentary cover.

Chemical and Sr isotopic compositions of pore fluids and marine carbonates from ODP Hole 130-807A

The 87Sr/86Sr ratios and Sr concentrations in sediment and pore fluids are used to evaluate the rates of calcite recrystallization at ODP Site 807A on the Ontong Java Plateau, an 800-meter thick section of carbonate ooze and chalk. A numerical model is used to evaluate the pore fluid chemistry and Sr isotopes in an accumulating section. The deduced calcite recrystallization rate is 2% per million years (%/Myr) near the top of the section and decreases systematically in older parts of the section such that the rate is close to 0.1/age (in years). The deduced recrystallization rates have important implications for the interpretation of Ca and Mg concentration profiles in the pore fluids. The effect of calcite recrystallization on pore fluid chemistry is described by the reaction length, L, which varies by element, and depends on the concentration in pore fluid and solid. When L is small compared to the thickness of the sedimentary section, the pore fluid concentration is controlled by equilibrium or steady-state exchange with the solid phase, except within a distance L of the sediment-water interface. When L is large relative to the thickness of sediment, the pore fluid concentration is mostly controlled by the boundary conditions and diffusion. The values of L for Ca, Sr, and Mg are of order 15, 150, and 1500 meters, respectively. L_Sr is derived from isotopic data and modeling, and allows us to infer the values of L_Ca and L_Mg. The small value for L_Ca indicates that pore fluid Ca concentrations, which gradually increase down section, must be equilibrium values that are maintained by solution-precipitation exchange with calcite and do not reflect Ca sources within or below the sediment column. The pore fluid Ca measurements and measured alkalinity allow us to calculate the in situ pH in the pore fluids, which decreases from 7.6 near the sediment-water interface to 7.1+/-0.1 at 400-800 mbsf. While the calculated pH values are in agreement with some of the values measured during ODP Leg 130, most of the measurements are artifacts. The large value for L_Mg indicates that the pore fluid Mg concentrations at 807A are not controlled by calcite-fluid equilibrium but instead are determined by the changing Mg concentration of seawater during deposition, modified by aqueous diffusion in the pore fluids. We use the pore fluid Mg concentration profile at Site 807A to retrieve a global record for seawater Mg over the past 35 Myr, which shows that seawater Mg has increased rapidly over the past 10 Myr, rather than gradually over the past 60 Myr. This observation suggests that the Cenozoic rise in seawater Mg is controlled by continental weathering inputs rather than by exchange with oceanic crust. The relationship determined between reaction rate and age in silicates and carbonates is strikingly similar, which suggests that reaction affinity is not the primary determinant of silicate dissolution rates in nature.

Pore water isotopic ratios of CO2 and SO4, and sediment CH4 content of ODP Leg 164 sites

A unique set of geochemical pore-water data, characterizing the sulfate reduction and uppermost methanogenic zones, has been collected at the Blake Ridge (offshore southeastern North America) from Ocean Drilling Program (ODP) Leg 164 cores and piston cores. The d13C values of dissolved CO2 (sum CO2) are as 13C-depleted as -37.7 per mil PDB (Site 995) at the sulfate-methane interface, reflecting a substantial contribution of isotopically light carbon from methane. Although the geochemical system is complex and difficult to fully quantify, we use two methods to constrain and illustrate the intensity of anaerobic methane oxidation in Blake Ridge sediments. An estimate using a two-component mixing model suggests that ~24% of the carbon residing in the sum CO2 pool is derived from biogenic methane. Independent diagenetic modeling of a methane concentration profile (Site 995) indicates that peak methane oxidation rates approach 0.005 µmol/cm**3/yr, and that anaerobic methane oxidation is responsible for consuming ~35% of the total sulfate flux into the sediments. Thus, anaerobic methane oxidation is a significant biogeochemical sink for sulfate, and must affect interstitial sulfate concentrations and sulfate gradients. Such high proportions of sulfate depletion because of anaerobic methane oxidation are largely undocumented in continental rise sediments with overlying oxic bottom waters. We infer that the additional amount of sulfate depleted through anaerobic methane oxidation, fueled by methane flux from below, causes steeper sulfate gradients above methane-rich sediments. Similar pore water chemistries should occur at other methane-rich, continental-rise settings associated with gas hydrates.

Depth profiles of pore-water and solid-phase constituents, respiration rates and grain size distribution in sediments of the Clarion-Clipperton Fracture Zone

Manganese nodules of the Clarion-Clipperton Fracture Zone (CCFZ) in the NE Pacific Ocean are highly enriched in Ni, Cu, Co, Mo and rare-earth elements, and thus may be the subject of future mining operations. Elucidating the depositional and biogeochemical processes that contribute to nodule formation, as well as the respective redox environment in both, water column and sediment, supports our ability to locate future nodule deposits and evaluates the potential ecological and environmental effects of future deep-sea mining. For these purposes we evaluated the local hydrodynamics and pore-water geochemistry with respect to the nodule coverage at four sites in the eastern CCFZ. Furthermore, we carried out selective leaching experiments at these sites in order to assess the potential mobility of Mn in the solid phase, and compared them with the spatial variations in sedimentation rates. We found that the oxygen penetration depth is 180 - 300 cm at all four sites, while reduction of Mn and NO3- is only significant below the oxygen penetration depth at sites with small or no nodules on the sediment surface. At the site without nodules, potential microbial respiration rates, determined by incubation experiments using 14C-labelled acetate, are slightly higher than at sites with nodules. Leaching experiments showed that surface sediments covered with big or medium-sized nodules are enriched in mobilizable Mn. Our deep oxygen measurements and pore-water data suggest that hydrogenetic and oxic-diagenetic processes control the present-day nodule growth at these sites, since free manganese from deeper sediments is unable to reach the sediment surface. We propose that the observed strong lateral contrasts in nodule size and abundance are sensitive to sedimentation rates, which in turn, are controlled by small-scale variations in seafloor topography and bottom-water current intensity.

Pore water chemistry of sediment core GeoB15101-7 and water chemistry of CTD casts GeoB15101-1 and GeoB15101-6

Submarine brine lakes feature sharp and persistent concentration gradients between seawater and brine, though these should be smoothed out by free diffusion in open ocean settings. The anoxic Urania basin of the Eastern Mediterranean contains an ultra sulfidic, hypersaline brine of Messinian origin above a thick layer of suspended sediments. With a dual modeling approach we reconstruct its contemporary stratification by geochemical solute transport fundamentals, and show that thermal convection is required to maintain mixing in the brine and mud layer. The origin of the Urania basin stratification was dated to 1650 years before present, which may be linked to a major earthquake in the region. The persistence of the chemoclines may be key to the development of diverse and specialized microbial communities. Ongoing thermal convection in the fluid mud layer may have important, yet unresolved consequences for sedimentological and geochemical processes, also in similar environments.

Geochemical data from pore water and solids from cores GS08-153-91GC12 and GS08-153-53-GC6, Arctic Mid-Ocean Ridge

Microbial communities and their associated metabolic activity in marine sediments have a profound impact on global biogeochemical cycles. Their composition and structure are attributed to geochemical and physical factors, but finding direct correlations has remained a challenge. Here we show a significant statistical relationship between variation in geochemical composition and prokaryotic community structure within deep-sea sediments. We obtained comprehensive geochemical data from two gravity cores near the hydrothermal vent field Loki's Castle at the Arctic Mid-Ocean Ridge, in the Norwegian-Greenland Sea. Geochemical properties in the rift valley sediments exhibited strong centimeter-scale stratigraphic variability. Microbial populations were profiled by pyrosequencing from 15 sediment horizons (59,364 16S rRNA gene tags), quantitatively assessed by qPCR, and phylogenetically analyzed. Although the same taxa were generally present in all samples, their relative abundances varied substantially among horizons and fluctuated between Bacteria- and Archaea-dominated communities. By independently summarizing covariance structures of the relative abundance data and geochemical data, using principal components analysis, we found a significant correlation between changes in geochemical composition and changes in community structure. Differences in organic carbon and mineralogy shaped the relative abundance of microbial taxa. We used correlations to build hypotheses about energy metabolisms, particularly of the Deep Sea Archaeal Group, specific Deltaproteobacteria, and sediment lineages of potentially anaerobic Marine Group I Archaea. We demonstrate that total prokaryotic community structure can be directly correlated to geochemistry within these sediments, thus enhancing our understanding of biogeochemical cycling and our ability to predict metabolisms of uncultured microbes in deep-sea sediments.

Chemical composition of sediment core GeoB12309-5 from the northern Arabian Sea, its pore water profile and age depth relationship as well as CTD data obtained during cruise M74/3

The Arabian Sea off the Pakistan continental margin is characterized by one of the world's largest oxygen minimum zones (OMZ). The lithology and geochemistry of a 5.3 m long gravity core retrieved from the lower boundary of the modern OMZ (956 m water depth) were used to identify late Holocene changes in oceanographic conditions and the vertical extent of the OMZ. While the lower part of the core (535 - 465 cm, 5.04 - 4.45 cal kyr BP, Unit 3) is strongly bioturbated indicating oxic bottom water conditions, the upper part of the core (284 - 0 cm, 2.87 cal kyr BP to present, Unit 1) shows distinct and well-preserved lamination, suggesting anoxic bottom waters. The transitional interval from 465 to 284 cm (4.45 - 2.87 cal kyr BP, Unit 2) contains relicts of lamination which are in part intensely bioturbated. These fluctuations in bioturbation intensity suggest repetitive changes between anoxic and oxic/suboxic bottom-water conditions between 4.45 - 2.87 cal kyr BP. Barium excess (Baex) and total organic carbon (TOC) contents do not explain whether the increased TOC contents found in Unit 1 are the result of better preservation due to low BWO concentrations or if the decreased BWO concentration is a result of increased productivity. Changes in salinity and temperature of the outflowing water from the Red Sea during the Holocene influenced the water column stratification and probably affected the depth of the lower boundary of the OMZ in the northern Arabian Sea. Even if we cannot prove certain scenarios, we propose that the observed downward shift of the lower boundary of the OMZ was also impacted by a weakened Somali Current and a reduced transport of oxygen-rich Indian Central Water into the Arabian Sea, both as a response to decreased summer insolation and the continuous southward shift of the Intertropical Convergence Zone during the late Holocene.

Uranium contents and isotope ratios of uranium and strontium in sediments, leachates and pore fluids from ODP Hole 162-984A

Bulk dissolution rates for sediment from ODP Site 984A in the North Atlantic are determined using the 234U/238U activity ratios of pore water, bulk sediment, and leachates. Site 984A is one of only several sites where closely spaced pore water samples were obtained from the upper 60 meters of the core; the sedimentation rate is high (11-15 cm/ka), hence the sediments in the upper 60 meters are less than 500 ka old. The sediment is clayey silt and composed mostly of detritus derived from Iceland with a significant component of biogenic carbonate (up to 30%). The pore water 234U/238U activity ratios are higher than seawater values, in the range of 1.2 to 1.6, while the bulk sediment 234U/238U activity ratios are close to 1.0. The 234U/238U of the pore water reflects a balance between the mineral dissolution rate and the supply rate of excess 234U to the pore fluid by a-recoil injection of 234Th. The fraction of 238U decays that result in a-recoil injection of 234U to pore fluid is estimated to be 0.10 to 0.20 based on the 234U/238U of insoluble residue fractions. The calculated bulk dissolution rates, in units of g/g/yr are in the range of 0.0000004 to 0.000002 1/yr. There is significant down-hole variability in pore water 234U/238U activity ratios (and hence dissolution rates) on a scale of ca. 10 m. The inferred bulk dissolution rate constants are 100 to 1000 times slower than laboratory-determined rates, 100 times faster than rates inferred for older sediments based on Sr isotopes, and similar to weathering rates determined for terrestrial soils of similar age. The results of this study suggest that U isotopes can be used to measure in situ dissolution rates in fine-grained clastic materials. The rate estimates for sediments from ODP Site 984 confirm the strong dependence of reactivity on the age of the solid material: the bulk dissolution rate (R_d) of soils and deep-sea sediments can be approximately described by the expression R_d ~ 0.1 1/age for ages spanning 1000 to 500,000,000 yr. The age of the material, which encompasses the grain size, surface area, and other chemical factors that contribute to the rate of dissolution, appears to be a much stronger determinant of dissolution rate than any single physical or chemical property of the system.