Coarse‐grained and all‐atom modeling of structural states and transitions in hemoglobin

Hemoglobin (Hb), an oxygen‐binding protein composed of four subunits (α1, α2, β1, and β2), is a well‐known example of allosteric proteins that are capable of cooperative ligand binding. Despite decades of studies, the structural basis of its cooperativity remains controversial. In this study, we have integrated coarse‐grained (CG) modeling, all‐atom simulation, and structural data from X‐ray crystallography and wide‐angle X‐ray scattering (WAXS), aiming to probe dynamic properties of the two structural states of Hb (T and R state) and the transitions between them. First, by analyzing the WAXS data of unliganded and liganded Hb, we have found that the structural ensemble of T or R state is dominated by one crystal structure of Hb with small contributions from other crystal structures of Hb. Second, we have used normal mode analysis to identify two distinct quaternary rotations between the α1β1 and α2β2 dimer, which drive the transitions between T and R state. We have also identified the hot‐spot residues whose mutations are predicted to greatly change these quaternary motions. Third, we have generated a CG transition pathway between T and R state, which predicts a clear order of quaternary and tertiary changes involving α and β subunits in Hb. Fourth, we have used the accelerated molecular dynamics to perform an all‐atom simulation starting from the T state of Hb, and we have observed a transition toward the R state of Hb. Further analysis of crystal structural data and the all‐atom simulation trajectory has corroborated the order of quaternary and tertiary changes predicted by CG modeling. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Remarkable disparity in mechanical response among the extracellular domains of type I and II cadherins

Cadherins, a large family of calcium-dependent adhesion molecules, are critical for intercellular adhesion. While crystallographic structures for several cadherins show clear structural similarities, their relevant adhesive strengths vary and their mechanisms of adhesion between types I and II cadherin subfamilies are still unclear. Here, stretching of cadherins was explored experimentally by atomic force microscopy and computationally by steered molecular dynamics (SMD) simulations, where partial unfolding of the E-cadherin ectodomains was observed. The SMD simulations on strand-swapping cadherin dimers displayed similarity in binding strength, suggesting contributions of other mechanisms to explain the strength differences of cell adhesion in vivo. Systematic simulations on the unfolding of the extracellular domains of type I and II cadherins revealed diverse pathways. However, at the earliest stage, a remarkable similarity in unfolding was observed for the various type I cadherins that was distinct from that for type II cadherins. This likely correlates positively with their distinct adhesive properties, suggesting that the initial forced deformation in type I cadherins may be involved in cadherin-mediated adhesion.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:25     

Comparative study of the binding mode between cytochrome P450 17A1 and prostate cancer drugs in the absence of haem iron

Abstract

According to the X-ray crystal structures of CYP17A1 (including its complexes with inhibitors), it is shown that a hydrogen bond exists between CYP17A1 and its inhibitors (such as abiraterone and TOK-001). Previous short MD simulations (50?ns) suggested that the binding of abiraterone to CYP17A1 is stronger than that of TOK-001. In this work, by carrying out long atomistic MD simulations (200?ns) of CYP17A1 and its complexes with abiraterone and TOK-001, we observed a binding mode between CYP17A1 and abiraterone, which is different from the binding mode between CYP17A1 and TOK-001. In the case of abiraterone binding, the unfilled volume in the active site cavity increases the freedom of movement of abiraterone within CYP17A1, leading to the collective motions of the helices G and B′ as well as the breaking of hydrogen bond existing between the 3β-OH group of abiraterone and N202 of CYP17A1. However, the unfilled volume in the active site cavity can be occupied by the benzimidazole ring of TOK-001, restraining the motion of TOK-001. By pulling the two inhibitors (abiraterone and TOK-001) out of the binding pocket in CYP17A1, we discovered that abiraterone and TOK-001 were moved from their binding sites to the surface of protein similarly through the channels formed by the helices G and B′. In addition, based on the free energy calculations, one can see that it is energetically favorable for the two inhibitors (abiraterone and TOK-001) to enter into the binding pocket in CYP17A1.     

Binding orientation and specificity of calmodulin to rat olfactory cyclic nucleotide-gated ion channel

Calmodulin (CaM), the primary intracellular Ca²⁺ receptor, regulates a large number of key enzymes and controls a wide spectrum of important biological responses. Recognition between CaM and its target sequence in rat olfactory cyclic nucleotide-gated ion channel (OLFp) was investigated by circular dichroism (CD), fluorescence, and NMR spectroscopy. Fluorescence data showed the OLFp tightly bound to CaM with a dissociation constant of 12?nM in a 1:1 stoichiometry. Far-UV CD data showed that approximately 60% of OLFp residues formed α-helical structures when associated with CaM. NMR data showed that most of the ¹⁵N–¹H HSQC cross-peaks of the ¹⁵N-labeled CaM not only shifted but also split into two sets of peaks upon association with the OLFp. Our data indicated that the two distinct CaM/OLFp complexes existed simultaneously with stable structures that were not interexchangeable within the NMR time scale. In light of the palindromic sequence of OLFp (FQRIVRLVGVIRDW) for CaM targeting, we proposed that the helical OLFp with C₂ symmetry may bind to CaM in two orientations. This hypothesis is supported by the observation that only one set of ¹⁵N–¹H HSQC cross-peaks of the ¹⁵N-labeled CaM was detected upon association with OLFp-M13 chimeric peptide (OLFMp), a mutated OLFp lacking the palindromic feature. The binding specificity of OLFMp to CaM was restored when the palindromic feature was destroyed. Binding modes of CaM/OLFp and CaM/OLFMp simulated by molecular docking were in accord with their distinct patterns observed in HSQC spectra. Our studies suggest that the palindromic residues in OLFp are crucial for the orientation-specific recognition by CaM.     

Resveratrol Binding to Human Serum Albumin

Abstract

Resveratrol (Res), a polyphenolic compound found largely in the skin of red grape and wine, exhibits a wide range of pharmaceutical properties and plays a role in prevention of human cardiovascular diseases [Pendurthi et al., Arterioscler. Thromb. Vasc. Biol. 19, 419–426 (1999)]. It shows a strong affinity towards protein binding and used as inhibitor for cyclo- oxygenase and ribonuclease reductase. The aim of this study was to examine the interaction of resveratrol with human serum albumin (HSA) in aqueous solution at physiological conditions, using a constant protein concentration (0.3 mM) and various pigment contents μM to mM). FTIR, UV-Visible, CD, and fluorescence spectroscopic methods were used to determine the resveratrol binding mode, the binding constant and the effects of pigment complexation on protein secondary structure.

Structural analysis showed that resveratrol bind non-specifically (H-bonding) via polypeptide polar groups with overall binding constant of K_Res = 2.56× 10⁵ M^?1. The protein secondary structure, analysed by CD spectroscopy, showed no major alterations at low resveratrol concentrations (0.125 mM), whereas at high pigment content (1 mM), major increase of α-helix from 57% (free HSA) to 62% and a decrease of β-sheet from 10% (free HSA) to 7% occurred in the resveratrol-HSA complexes. The results indicate a partial stabilization of protein secondary structure at high resveratrol content.     

Flux analysis and metabolomics for systematic metabolic engineering of microorganisms

Rational engineering of metabolism is important for bio-production using microorganisms. Metabolic design based on in silico simulations and experimental validation of the metabolic state in the engineered strain helps in accomplishing systematic metabolic engineering. Flux balance analysis (FBA) is a method for the prediction of metabolic phenotype, and many applications have been developed using FBA to design metabolic networks. Elementary mode analysis (EMA) and ensemble modeling techniques are also useful tools for in silico strain design. The metabolome and flux distribution of the metabolic pathways enable us to evaluate the metabolic state and provide useful clues to improve target productivity. Here, we reviewed several computational applications for metabolic engineering by using genome-scale metabolic models of microorganisms. We also discussed the recent progress made in the field of metabolomics and ¹³C-metabolic flux analysis techniques, and reviewed these applications pertaining to bio-production development. Because these in silico or experimental approaches have their respective advantages and disadvantages, the combined usage of these methods is complementary and effective for metabolic engineering.     

Effects of Age,Sex, Index Admission,and Predominant Polarity on the Seasonality of Acute Admissions For Bipolar Disorder: A Population-Based Study

Bipolar disorder seasonality has been documented previously, though information on the effect of demographic and clinical variables on seasonal patterns is scant. This study examined effects of age, sex, index admission, and predominant polarity on bipolar disorder seasonality in a nationwide population. An inpatient cohort admitted to hospital exclusively for mental illness was derived from the Taiwan National Health Insurance Research Database for 2002–2007. The authors identified 9619 inpatients with bipolar disorder, who had generated 15 078 acute admission records. An empirical mode decomposition method was used to identify seasonal oscillations in bipolar admission data, and regression and cross-correlation analyses were used to quantify the degree and timing of bipolar admission seasonality. Results for seasonality timing found that manic or mixed episodes peak in spring or summer, and depressive episodes peak in winter. Analysis for degree of seasonality revealed that (1) the polarity of patients' index admission predicted the seasonality of relapse admissions; (2) seasonality was significant in female admissions for depressive episodes and in male admissions for manic episodes; (3) young adults displayed a higher degree of seasonality for acute admissions than middle-aged adults; and (4) patients with predominantly depressive admissions displayed a higher degree of seasonality than patients with predominantly manic admissions. Demographic and clinical variables were found to affect the seasonality of acute admissions for bipolar disorders. These findings highlight the need for research on identification and management of seasonal features in bipolar patients. (Author correspondence: ccyang@physionet.org)     

湖南工程学院生物工程专业实践 教学模式创新与实践

分析生物工程专业人才需求的现状与发展前景,构建了以社会需求为导向,以实践课程体系为载体,以能力培养为主线,以学生科技创新活动和教师科研为驱动,项目研究与能力培养相结合的创新实践教学模式.实践表明,创新实践教学模式极大提高了学生的学习兴趣和教师的工作热情.     

PATERNITY PROTECTION CAN PROVIDE A KICK‐START FOR THE EVOLUTION OF MALE‐ONLY PARENTAL CARE

Sperm competition and uncertainty of paternity hamper the evolution of male parental care. Thus, maternal care predominates in most taxa. What if males can, however, limit cuckoldry by guarding the eggs postmating? Here, we show that this provides a reason to reconsider an old and nowadays rather discredited hypothesis: that external fertilization is associated with male care because the parent who releases its gametes first can depart leaving the other in a “cruel bind,” having to care for the offspring. In our model, protection of paternity provides an additional incentive for the male to stay associated with its young. When we then assume that offspring survive better if guarded, paternity protection proves enough to kick‐start the evolution of male‐only parental care from a scenario with no care. This fits with data from fishes, where male‐only care is associated with external fertilization, whereas female‐only care almost always evolves after an initial transition to internal fertilization. Our model unifies disparate hypotheses regarding parental care roles and provides support for the idea that care roles can be influenced by sex differences in selection to be physically close to the offspring, including selection that is initially not based on offspring survival.     

广西爬树蕨属两种植物的细胞分类学研究

首次报道了中国广西产两种爬树蕨属植物的染色体数目及繁殖方式。爬树蕨和桂南爬树蕨的染色体数目都是n=40,2n=80(2x); 每个孢子囊的孢子数为64,是有性生殖二倍体。结合已有资料对爬树蕨属中的系统位置进行了讨论。