Volumetric properties underlying ligand binding in a monomeric hemoglobin: A high-pressure NMR study

The 2/2 hemoglobin of the cyanobacterium Synechococcus sp. PCC 7002, GlbN, coordinates the heme iron with two histidines and exists either with a b heme or with a covalently attached heme. The binding of exogenous ligands displaces the distal histidine and induces a conformational rearrangement involving the reorganization of internal void volumes. The formation of passageways within the resulting conformation is thought to facilitate ligand exchange and play a functional role. Here we monitored the perturbation induced by pressure on the ferric bis-histidine and cyanide-bound states of GlbN using ¹H–¹⁵N HSQC NMR spectroscopy. We inspected the outcome with a statistical analysis of 170 homologous 2/2 hemoglobin sequences. We found that the compression landscape of GlbN, as represented by the variation of an average chemical shift parameter, was highly sensitive to ligand swapping and heme covalent attachment. Stabilization of rare conformers was observed at high pressures and consistent with cavity redistribution upon ligand binding. In all states, the EF loop was found to be exceptionally labile to pressure, suggesting a functional role as a semi-flexible hinge between the adjacent helices. Finally, coevolved clusters presented a common pattern of compensating pressure responses. The high-pressure dissection combined with protein sequence analysis established locations with volumetric signatures relevant to residual communication of 2/2 hemoglobins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.     

Enhanced sampling of protein conformational states for dynamic cross-docking within the protein-protein docking server SwarmDock

The formation of specific protein-protein interactions is often a key to a protein's function. During complex formation, each protein component will undergo a change in the conformational state, for some these changes are relatively small and reside primarily at the sidechain level; however, others may display notable backbone adjustments. One of the classic problems in the protein-docking field is to be able to a priori predict the extent of such conformational changes. In this work, we investigated three protocols to find the most suitable input structure conformations for cross-docking, including a robust sampling approach in normal mode space. Counterintuitively, knowledge of the theoretically best combination of normal modes for unbound-bound transitions does not always lead to the best results. We used a novel spatial partitioning library, Aether Engine (see Supplementary Materials ), to efficiently search the conformational states of 56 receptor/ligand pairs, including a recent CAPRI target, in a systematic manner and selected diverse conformations as input to our automated docking server, SwarmDock, a server that allows moderate conformational adjustments during the docking process. In essence, here we present a dynamic cross-docking protocol, which when benchmarked against the simpler approach of just docking the unbound components shows a 10% uplift in the quality of the top docking pose.     

Variations in genotype–phenotype correlations in phenylalanine hydroxylase deficiency in Chinese Han population

Background

The value of genotyping to predict variant phenotypes in patients with phenylalanine hydroxylase (Pah) deficiency is a matter of debate. However, there exists no comprehensive population relationship study focused on the Han Chinese.

Methods

We analyzed genotype–phenotype correlation for 186 different genotypes in 338 unrelated Chinese patients harboring 109 different Pah mutations. Two systems were used in this process. The first was a phenotype prediction system based on arbitrary values (AV) attributed to each mutation. The second was a pair-wise correlation analysis. The observed phenotype for AV analysis was the corresponding metabolic phenotype stratified according to the pretreatment phenylalanine (Phe) value.

Results

We found that the observed phenotype matched the predicted phenotype in 54.41% of 272 patients for whom AV information was available; the highest degree of concordance (61.83%) was found in patients with null/null genotypes, whereas the lowest “concordance rate” (32.69%) was observed for patients with expected mild-PKU phenotype. There are repeated inconsistencies for such mutations as R241C, R243Q, R261Q, V388M, V399V, R408Q, A434D and EX6-96A>G which are associated with variable phenotypes in patients with identical genotype. Significant correlations were disclosed between pretreatment Phe values and predicted residual activity (r = − 0.45643, P < 0.0001) or AV sum (r = − 0.59523, P < 0.0001).

Conclusion

Our study supports the notion that the Pah mutation genotype is the main determinant of metabolic phenotype in most patients in a particular population, and provided novel insights into the values that underpin the subsequent treatment and the prognosis of PKU in Chinese.     

Molecular dynamics simulations of heat conduction in multi-walled carbon nanotubes

Heat conduction in multi-walled carbon nanotubes (MWNTs) was studied using non-equilibrium molecular dynamics simulations. This research focuses on the effects of the multi-wall structure of the MWNTs on the heat conduction. The results show that the thermal conductivity of a MWNT is almost the same as that of the corresponding single-walled carbon nanotubes (SWNTs) rather than much smaller as has been suggested. Thus, the multi-wall structure does not significantly affect the thermal conduction in the MWNTs. Analysis of the temperature profiles and the phonon density of states confirms that there is almost no heat transport between the MWNT layers and that each layer conducts heat nearly independently along parallel channels. This is physically reasonable since the weak inter-wall interactions and large interfacial thermal resistances make the MWNT layers behave like parallel thermal circuits.     

New software for finding transition states by probing accessible,or ergodic,regions

An automated, iterative approach to finding the lowest energy, ionic diffusion paths through a periodic structure has been developed within our new code (written in FORTRAN 77 and named Bubble). The approach is quite general in that it can be applied to find, at a chosen temperature, the accessible (ergodic) regions of a hyper-surface, which is defined across a uniform grid [1 Schön, J.C., Putz, H. and Jansen, M. 1996. Studying the energy hypersurface of continuous system—The threshold algorithm. J. Phys.-Conden. Matt., 8: 143[Crossref] , [Google Scholar]]. We describe both our implementation within the Bubble code and its application to locating the approximate transition states for Mg interstitial diffusion in forsterite, which can then be refined using standard transition state searching [2 Banerjee, A., Adams, N., Simons, J. and Shepard, R. 1985. Search for stationary points on surfaces. J. Phys. Chem., 89: 52[Crossref], [Web of Science ®] , [Google Scholar]].     

Subunit Secondary Structure in Filamentous Viruses: Predictions and Observations

Abstract

The algorithm of Gamier, Osguthorpe and Robson (J. Mol. Biol. 120, 97–120, 1978) for prediction of protein secondary structure has been applied to the coat protein sequences of six filamentous bacteriophages: fd, Ifl, IKe, Pfl, Xf and Pf3. For subunits of Class I virions (fd, Ifl, IKe), the algorithm predicts a very high percentage of helix in comparison to other structure types, which is in accord with the results of laser Raman and circular dichroism measurements. For subunits of the Class II virions (Pfl, Xf, Pf3), the algorithm consistently predicts a predominance of β structure, which is compatible with the demonstrated facility for conversion of Class II subunits from α-helix to β-strand under appropriate experimental conditions (Thomas, Prescott and Day, J. Mol. Biol. 165, 321–356, 1983). Even when the algorithm is biased to favor helix, the Class II virion subunits are predicted to contain considerably more strand than helix. Qualitatively similar results are obtained using the algorithm of Chou and Fasman {Adv. Enzym. 47, 45–148,45-148). Therefore, both predictive and experimental methods indicate a distinction between Gass I and II subunits, which is reflected in a greater tendency of the latter to adopt other than uniform β-helical conformation. The results suggest a possible model for the disassembly of filamentous viruses which may involve the unraveling of coat protein helices at the N terminus.     

Theophylline Steady-State Pharmacokinetics: Recent Concepts and Their Application in Chronotherapy of Reactive Airway Diseases

With the introduction of sustained-release theophyllíne formulations for once-daily dosing or for unequally divided twice-daily dosing, comparison with conventional equally divided twice-daily dosing has been focused on nocturnal serum theophylline concentrations (STCs), plateau properties and peak-trough fluctuation. The merits of various steady-state characteristics such as nocturnal excess, plateau time, residual concentration, peak-trough fluctuation, swing and AUC fluctuation are illustrated by 15 data sets from 7 multiple-dose studies, each including either 10–12 healthy volunteers or 12–20 COPD-patients. In all of the studies, STCs were determined at least every 2 hr over a 24-hr period in steady-state. Included in the studies were 7 sustained-release theophylline formulations which were administered either once daily (in the morning or in the evening), or twice daily (either equally divided, or unequally divided with one-third of the dose being given in the morning and two-thirds in the evening).     

Ocean energy and the law of the sea: The need for a protocol

Although fossil fuels are the overwhelming source of energy for the world, and will continue to be so for the foreseeable future, demographic, environmental, political, and economic factors indicate that interest in alternative, renewable sources of energy will grow. There is a need for both global and national policies on ocean energy management. In particular, coastal states and the energy industry would benefit from guidelines that helped to create a predictable, stable environment in which long‐term, high‐cost research, development, and investment decisions could be made with confidence. Coastal states have jurisdiction over the maritime zones most relevant to energy production, but many lack the expertise and funds to develop this potential source. Industry must operate within the control of coastal states and will not be served well by a plethora of differing legal interpretations and unilaterally imposed restrictions and obligations from state to state. An Ocean Energy Protocol to the 1982 UN Convention on the Law of the Sea would afford governments and industry the opportunity to clarify their respective obligations and address particular interests for mutual benefit.