Revisiting “reverse hydrophobic effect”: Applicable only to coil mutations at the surface

In a seminal paper, Pakula and Sauer (Nature, 1990, 344, 363–364) demonstrated that the increase in side‐chain hydrophobicity has a reverse relationship with protein stability. We have addressed this problem with several examples of mutants that span at different locations in protein structure based on secondary structure and solvent accessibility. We confirmed that the stability change upon single coil mutation at exposed region is reversely correlated with hydrophobicity with a single exception. In addition, we found the existence of such relationship in partially buried coil mutants. The stability of exposed helical mutants is governed by conformational properties. In buried and partially buried helical and strand mutants properties reflecting hydrophobicity have direct relationship with stability, whereas an opposite relationship was obtained with entropy and flexibility. The structural analysis of partially buried/exposed mutants showed that the surrounding residues are important for the stability change upon mutation. These results provide insights to understand the general behavior for the stability of proteins upon amino acid substitutions. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 591–599, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Global dynamics of two-compartment models for cell production systems with regulatory mechanisms

We present a global stability analysis of two-compartment models of a hierarchical cell production system with a nonlinear regulatory feedback loop. The models describe cell differentiation processes with the stem cell division rate or the self-renewal fraction regulated by the number of mature cells. The two-compartment systems constitute a basic version of the multicompartment models proposed recently by Marciniak-Czochra and collaborators [25] to investigate the dynamics of the hematopoietic system. Using global stability analysis, we compare different regulatory mechanisms. For both models, we show that there exists a unique positive equilibrium that is globally asymptotically stable if and only if the respective reproduction numbers exceed one. The proof is based on constructing Lyapunov functions, which are appropriate to handle the specific nonlinearities of the model. Additionally, we propose a new model to test biological hypothesis on the regulation of the fraction of differentiating cells. We show that such regulatory mechanism is incapable of maintaining homeostasis and leads to unbounded cell growth. Potential biological implications are discussed.     

The role of transmissible diseases in the Holling-Tanner predator-prey model

Transmissible diseases are known to induce remarkable major behavioral changes in predator-prey systems. However, little attention has been paid to model such situations. The latter would allow to predict useful applications in both dynamics and control. Here the Holling-Tanner model is revisited to account for the influence of a transmissible disease, under the assumption that it spreads among the prey species only. We have found the equilibria and analyzed the behavior of the system around each one of them. A threshold result determining when the disease dies out has been identified. We also investigated the parametric space under which the system enters into Hopf and transcritical bifurcations, around the disease free equilibrium. The system is shown to experience neither saddle-node nor pitch-fork bifurcation. Global stability results are obtained by constructing suitable Lyapunov functions.     

Insight into the stability of cross‐β amyloid fibril from molecular dynamics simulation

Amyloid fibrils are considered to play causal roles in the pathogenesis of amyloid‐related degenerative diseases such as Alzheimer's disease, type II diabetes mellitus, the transmissible spongiform encephalopathies, and prion disease. The mechanism of fibril formation is still hotly debated and remains an important open question. In this study, we utilized molecular dynamics (MD) simulation to analyze the stability of hexamer for eight class peptides. The MD results suggest that VEALYL and MVGGVV‐1 are the most stable ones, then SNQNNY, followed by LYQLEN, MVGGVV‐2, VQIVYK, SSTSAA, and GGVVIA. The statistics result indicates that hydrophobic residues play a key role in stabilizing the zipper interface. Single point and two linkage mutants of MVGGVV‐1 confirmed that both Met1 and Val2 are key hydrophobic residues. This is consistent with the statistics analysis. The stability results of oligomer for MVGGVV‐1 suggest that the intermediate state should be trimer (3‐0) and tetramer (2‐2). These methods can be used in stabilization study of other amyloid fibril. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 578–586, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

具有多时变分布时滞的Lurie控制系统的绝对稳定性

研究了一类具有多时变分布时滞的Lurie控制系统的绝对稳定性.通过构造合适的Lyapunov函数,并利用矩阵不等式技术,获得了该控制系统绝对稳定性的充分条件.     

Conformational stability of the full‐atom hexameric model of the ClpB chaperone from Escherichia coli

The Escherichia coli heat shock protein ClpB, a member of the Hsp100 family, plays a crucial role in cellular thermotolerance. In co‐operation with the Hsp70 chaperone system, it is able to solubilize proteins aggregated by heat shock conditions and refold them into the native state in an ATP‐dependent way. It was established that the mechanism of ClpB action depends on the formation of a ring‐shaped hexameric structure and the translocation of a protein substrate through an axial channel. The structural aspects of this process are not fully known. By means of homology modeling and protein–protein docking, we obtained a model of the hexameric arrangement of the full‐length ClpB protein complexed with ATP. A molecular dynamics simulation of this model was performed to assess its flexibility and conformational stability. The high mobility of the “linker” M‐domain, essential for the renaturing activity of ClpB, was demonstrated, and the size and shape of central channel were analyzed. In this model, we propose the coordinates for a loop between b4 and B6 structural elements, not defined in previous structural research, which faces the inside of the channel and may therefore play a role in substrate translocation. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 47–60, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Synthesis of quadruplex‐forming tetra‐end‐linked oligonucleotides: Effects of the linker size on quadruplex topology and stability

G‐quadruplexes are characteristic structural arrangements of guanine‐rich DNA sequences that abound in regions with relevant biological significance. These structures are highly polymorphic differing in the number and polarity of the strands, loop composition, and conformation. Furthermore, the cation species present in solution strongly influence the topology of the G‐quadruplexes. Recently, we reported the synthesis and structural studies of new G‐quadruplex forming oligodeoxynucleotides (ODNs) in which the 3′‐ and/or the 5′‐ends of four ODN strands are linked together by a non‐nucleotidic tetra‐end‐linker (TEL). These TEL‐ODN analogs having the sequence TGGGGT are able to form parallel G‐quadruplexes characterized by a remarkable high thermal stability. We report here an investigation about the influence of the reduction of the TEL size on the molecularity, topology, and stability of the resulting TEL‐G‐quadruplexes using a combination of circular dichroism (CD), CD melting, ¹H NMR spectroscopy, gel electrophoresis, and molecular modeling data. We found that all TEL‐(TGGGGT)₄ analogs, regardless the TEL size and the structural orientation of the ODN branches, formed parallel TEL‐G‐quadruplexes. The molecular modeling studies appear to be consistent with the experimental CD and NMR data revealing that the G‐quadruplexes formed by TEL‐ODNs having the longer TEL (L 1 ‐ 4 ) are more stable than the corresponding G‐quadruplexes having the shorter TEL (S 1 ‐ 4 ). The relative stability of S 1 ‐ 4 was also reported. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 466–477, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Structure and conformational stability of a tetrameric thermostable N‐succinylamino acid racemase

The N‐succinylamino acid racemases (NSAAR) belong to the enolase superfamily and they are large homooctameric/hexameric species that require a divalent metal ion for activity. We describe the structure and stability of NSAAR from Geobacillus kaustophilus (GkNSAAR) in the absence and in the presence of Co²⁺ by using hydrodynamic and spectroscopic techniques. The Co²⁺, among other assayed divalent ions, provides the maximal enzymatic activity at physiological pH. The protein seems to be a tetramer with a rather elongated shape, as shown by AU experiments; this is further supported by the modeled structure, which keeps intact the largest tetrameric oligomerization interfaces observed in other homooctameric members of the family, but it does not maintain the octameric oligomerization interfaces. The native functional structure is mainly formed by α‐helix, as suggested by FTIR and CD deconvoluted spectra, with similar percentages of structure to those observed in other protomers of the enolase superfamily. At low pH, the protein populates a molten‐globule‐like conformation. The GdmCl denaturation occurs through a monomeric intermediate, and thermal denaturation experiments indicate a high thermostability. The presence of the cofactor Co²⁺ did alter slightly the secondary structure, but it did not modify substantially the stability of the protein. Thus, GkNSAAR is one of the few members of the enolase family whose conformational propensities and stability have been extensively characterized. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 757–772, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Calculating biological behaviors of epigenetic states in the phage λ life cycle

The biology and behavior of bacteriophage regulation have been the focus of classical investigations of molecular control of gene expression. Both qualitative and quantitative aspects of this behavior have been systematically characterized experimentally. Complete understanding of the robustness and stability of the genetic circuitry for the lysis-lysogeny switch remains an unsolved puzzle. It is an excellent test case for our understanding of biological behavior of an integrated network based on its physical, chemical, DNA, protein, and functional properties. We have used a new approach to non-linear dynamics to formulate a new mathematical model, performed a theoretical study on the phage life cycle, and solved the crucial part of this puzzle. We find a good quantitative agreement between the theoretical calculation and published experimental observations in the protein number levels, the lysis frequency in the lysogen culture, and the lysogenization frequency for mutants of O_R. We also predict the desired robustness for the genetic switch. We believe that this is the first successful example in the quantitative calculation of robustness and stability of the phage regulatory network, one of the simplest and most well-studied regulatory systems.     

Characterization of the mutant lytic state in lambda expression systems

The two propagative phases of bacteriophage lambda, lysogeny and lysis, can be used in concert to enhance productivity of recombinant expression systems. Lambda vectors carrying mutations to prevent both cell lysis and lambda DNA packaging in the lytic state have been shown to yield 100% stability of the product gene in lysogeny and to produce up to 15% of total cell protein as product beta-galactosidase in a mutant lytic state.(14) Despite these mutations, partial lysis of the culture was observed following induction of the cells from a lysogenic phase into the lytic state. To understand better the phage-host cell interactions and to investigate the possible cause(s) of lysis in these highly productive expression systems, we have made a detailed study of the suppressor-free system JM105(NM1070). We have found high levels of product (15% of total cell protein as beta-glactosidase) to be due chiefly to a high-copy number of lambda DNA in the mutant lytic state. There is partial lysis of the culture even in this suppressor-free system caused by a low-level natural suppression of the amber mutation in gene S of NM1070, resulting in accumulation of lambda endolysin. We have also monitored changes in cell growth and morphology upon induction of the lysogen. There is a slight increase in cell number that follows a linear relationship with time and a 25-fold increase in cell volume during recombinat protein production in the mutant lytic state.     

Stability in a multi-species assemblage of large herbivores in East Africa: an alternative interpretation

Summary Prins and Douglas-Hamilton (1990) analyzed data based on nine census counts of large herbivore species in Lake Manyara National Park, northern Tanzania, over the period 1959–1984. Their major conclusion was that even if individual species-showed large fluctuations in numbers, the different species compensated the fluctuations of the other species in a way resulting in a constancy of total herbivore biomass, constancy of plant biomass consumption, and overall stability of the system under natural conditions. The authors believed that they had found a support for this view by calculating the stability index based on correlations between numbers of large herbivores. In this paper I show that Prins and Douglas-Hamilton's calculation of the stability index was not justified. Grazing and browsing pressure by large herbivores in Lake Manyara National Park seems to be remarkably constant. However, available information does not allow any rigorous conclusions about the stability of this community. We need more data from other systems dominated by large herbivores to be able to make comparisons and to be able to say which systems are more stable and in what sense.     

Sequence organization of the origins of DNA replication in lambdoid coliphages

We have determined the sequences of the ori region DNA of several phage lambda mutants and hybrids, which shed light on the mechanism of DNA replication in the lambdoid phages. These include the heterologous substitution hybrids lambda rep82:lambda and lambda rep80:lambda, a pseudorevertant of the ori-r93 mutant lambda r93hot5, and the insertion mutant lambda pk35. The ori regions of the three lambdoid phages, lambda, phi 80 and 82, all have repeated sequences, termed iterons, and A . T-rich zones. We note that a similar arrangement of DNA is also found in several other prokaryotic origins of replication. lambda and phi 80 have four iterons, and 82 has five. The origin of lambda r93hot5 is unusual in that contains only three iterons, yet the phage grows normally. Analysis of this mutant indicates that the spacing of iterons is crucial to ori function, whereas their number is not. This argues against the cloverleaf model for lambda ori structure (Hobom et al., 1979). In lambda pk35 the drug resistance element Tn903 is inserted into the "inceptor" (ice) site, proposed to be crucial for lambda replication initiation (Hobom et al., 1979); yet this phage grows normally.     

Mg2+ ions bind at the C‐terminal region of skeletal muscle α‐tropomyosin

Tropomyosin (Tm) is a dimeric coiled‐coil protein that polymerizes through head‐to‐tail interactions. These polymers bind along actin filaments and play an important role in the regulation of muscle contraction. Analysis of its primary structure shows that Tm is rich in acidic residues, which are clustered along the molecule and may form sites for divalent cation binding. In a previous study, we showed that the Mg²⁺‐induced increase in stability of the C‐terminal half of Tm is sensitive to mutations near the C‐terminus. In the present report, we study the interaction between Mg²⁺ and full‐length Tm and smaller fragments corresponding to the last 65 and 26 Tm residues. Although the smaller Tm peptide (Tm_{259‐284(W269)}) is flexible and to large extent unstructured, the larger Tm_{220‐284(W269)} fragment forms a coiled coil in solution whose stability increases significantly in the presence of Mg²⁺. NMR analysis shows that Mg²⁺ induces chemical shift perturbations in both Tm_{220‐284(W269)} and Tm_{259‐284(W269)} in the vicinity of His276, in which are located several negatively charged residues. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 583–590, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

一类离散双线性生态系统的正规化控制

双线性系统能对非线性系统很好的近似,可对生态、生物等过程中的许多现象进行描述,对双线性系统研究具有一定的实际价值与理论意义.本文研究了一类离散双线性生态系统的全局渐近稳定问题.针对该系统,给出了一种简单的正规化反馈控制律.运用Lyapunov稳定性理论证明了在此控制律下的闭环系统是全局渐近稳定的.     

Impact of α‐hydroxy‐propanodeoxyguanine adducts on DNA duplex energetics: Opposite base modulation and implications for mutagenicity and genotoxicity

Acrolein is an α,β‐unsaturated aldehyde that is a major environmental pollutant, as well as a product of cellular metabolism. DNA bases react with acrolein to form two regioisomeric exocyclic guanine adducts, namely γ‐hydroxy‐propanodeoxyguanosine (γ‐OH‐PdG) and its positional isomer α‐hydroxy‐propanodeoxyguanosine (α‐OH‐PdG). The γ‐OH‐PdG isomer adopts a ring‐opened conformation with minimal structural perturbation of the DNA host duplex. Conversely, the α‐OH‐PdG isomer assumes a ring‐closed conformation that significantly disrupts Watson‐Crick base‐pair alignments within the immediate vicinity of the damaged site. We have employed a combination of calorimetric and spectroscopic techniques to characterize the thermodynamic origins of these lesion‐induced structural alterations. Specifically, we have assessed the energetic impact of α‐OH‐PdG centered within an 11‐mer duplex by hybridizing the adduct‐containing oligonucleotide with its complementary strand harboring a central base N [where N = C or A], yielding a pair of duplexes containing the nascent lesion (α‐OH‐PdG·C) or mismatched adduct (α‐OH‐PdG·A), respectively. Our data reveal that the nascent lesion is highly destabilizing, whereas its mismatched counterpart partially ameliorates α‐OH‐PdG‐induced destabilization. Collectively, our data provide energetic characterizations of the driving forces that modulate error‐free versus error‐prone DNA translesion synthesis. The biological implications of our findings are discussed in terms of energetically probing acrolein‐mediated mutagenicity versus adduct‐induced genotoxicity. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 370–382, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

The diversity of FtsY‐lipid interactions

The bacterial signal recognition particle (SRP) receptor FtsY forms a complex with the SRP Ffh to target nascent polypeptide chains to the bacterial inner membrane. How FtsY interacts with lipids and associates to the membrane is unclear. Here, we show that vesicle binding leads to partial protection against proteolytic degradation and a change in secondary structure, which differs depending on whether the lipids are simple mixtures of zwitterionic and anionic lipids, mimics of Escherichia coli lipids, or lysolipids. Lipid binding alters the stability of FtsY. Thermal unfolding of FtsY in buffer shows two transitions, one occurring at ～60°C and the other at ～90°C. The thermal intermediate accumulating between 60 and 90°C has structural features in common with the state induced by binding to E. coli lipids. E. coli lipid extract induces a single transition around 70°C, anionic lipids have no effect while cooperative unfolding is completely removed in lysolipids. Thus, the lipid environment profoundly influences the dynamic properties of FtsY, leading to three different kinds of FtsY‐lipid interactions with different effects on structure, proteolytic protection, and stability, and is driven both by hydrophobic and electrostatic interactions. Trypsin digestion experiments highlight the central role of the N‐domain in lipid contacts, whereas the A‐ and G‐domains appear to play a more minor part. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 595–606, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Computational analysis of glycoside hydrolase family 1 specificities

Glycoside hydrolase family 1 consists of beta-glucosidases, beta-galactosidases, 6-phospho-beta-galactosidases, myrosinases, and other enzymes having similar primary and tertiary structures but diverse specificities. Among these enzymes, beta-glucosidases hydrolyze cellobiose to glucose, and therefore they are key players in any cellulose to glucose process. All family members attack beta-glycosidic bonds between a pyranosyl glycon and an aglycon, but most have little specificity for the aglycon or for the bond configuration. Furthermore, glycon specificity is not absolute. Sixteen family members (six beta-glucosidases, two cyanogenic beta-glucosidases, one 6-phospho-beta-galactosidase, two myrosinases, and five beta-glycosidases) have known tertiary structures. We have used automated docking to computationally bind disaccharides with allopyranosyl, galactopyranosyl, glucopyranosyl, mannopyranosyl, 6-phosphogalactopyranosyl, and 6-phosphoglucopyranosyl glycons, all linked by beta-(1,2), beta-(1,3), beta-(1,4), and beta-(1,6)-glycosidic bonds to beta-glucopyranoside aglycons, along with beta-(1,1-thio)-allopyranosyl, -galactopyranosyl, -glucopyranosyl, and -mannopyranosyl) beta-glucopyranosides, into all of these structures to investigate the structural determinants of their enzyme specificities. The following are the eight active-site residues: Glu191, Thr194, Phe205, Asn285, Arg336, Asn376, Trp378, and Trp465 (Zea mays beta-glucosidase numbering), that control a significant amount of glycon specificity. (c) 2008 Wiley Periodicals, Inc. Biopolymers 89: 1021-1031, 2008.This article was originally published online as an accepted preprint. The "Published Online" date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com.     

时滞大系统鲁棒稳定性补充结果

利用代数Ｒｉｃｃａｔｉ方法和Ｌｙａｐｕｎｏｖ稳定性,本文分别得到带有匹配条件和非匹配条件的时滞大系统的鲁棒稳定性判别准则,通过仿真实例进一步说明本文所得稳定性的容许界限优于文献中的结论。     

Stability engineering of anti-EGFR scFv antibodies by rational design of a lambda-to-kappa swap of the VL framework using a structure-guided approach

Phage-display technology facilitates rapid selection of antigen-specific single-chain variable fragment (scFv) antibodies from large recombinant libraries. ScFv antibodies, composed of a V_H and V_L domain, are readily engineered into multimeric formats for the development of diagnostics and targeted therapies. However, the recombinant nature of the selection strategy can result in V_H and V_L domains with sub-optimal biophysical properties, such as reduced thermodynamic stability and enhanced aggregation propensity, which lead to poor production and limited application. We found that the C10 anti-epidermal growth factor receptor (EGFR) scFv, and its affinity mutant, P2224, exhibit weak production from E. coli. Interestingly, these scFv contain a fusion of lambda3 and lambda1 V-region (LV3 and LV1) genes, most likely the result of a PCR aberration during library construction. To enhance the biophysical properties of these scFvs, we utilized a structure-based approach to replace and redesign the pre-existing framework of the V_L domain to one that best pairs with the existing V_H. We describe a method to exchange lambda sequences with a more stable kappa3 framework (KV3) within the V_L domain that incorporates the original lambda DE-loop. The resulting scFvs, C10KV3_LV1DE and P2224KV3_LV1DE, are more thermodynamically stable and easier to produce from bacterial culture. Additionally, C10KV3_LV1DE and P2224KV3_LV1DE retain binding affinity to EGFR, suggesting that such a dramatic framework swap does not significantly affect scFv binding. We provide here a novel strategy for redesigning the light chain of problematic scFvs to enhance their stability and therapeutic applicability.