Fractal and complex network analyses of protein molecular dynamics

Based on protein molecular dynamics, we investigate the fractal properties of energy, pressure and volume time series using the multifractal detrended fluctuation analysis (MF-DFA) and the topological and fractal properties of their converted horizontal visibility graphs (HVGs). The energy parameters of protein dynamics we considered are bonded potential, angle potential, dihedral potential, improper potential, kinetic energy, Van der Waals potential, electrostatic potential, total energy and potential energy. The shape of the h(q)h(q) curves from MF-DFA indicates that these time series are multifractal. The numerical values of the exponent h(2)h(2) of MF-DFA show that the series of total energy and potential energy are non-stationary and anti-persistent; the other time series are stationary and persistent apart from series of pressure (with H≈0.5H≈0.5 indicating the absence of long-range correlation). The degree distributions of their converted HVGs show that these networks are exponential. The results of fractal analysis show that fractality exists in these converted HVGs. For each energy, pressure or volume parameter, it is found that the values of h(2)h(2) of MF-DFA on the time series, exponent λλ of the exponential degree distribution and fractal dimension dBdB of their converted HVGs do not change much for different proteins (indicating some universality). We also found that after taking average over all proteins, there is a linear relationship between 〈h(2)〉〈h(2)〉 (from MF-DFA on time series) and 〈dB〉〈dB〉 of the converted HVGs for different energy, pressure and volume.

Nutritional status and the foraging behaviour of Bactrocera tryoni with particular reference to protein bait spray

Poisoned protein baits comprise a recognized method for controlling tephritid fruit flies in the form of a ‘lure-and-kill’ technique. However, little is known about how a fly's internal protein and carbohydrate levels (i.e. nutritional status) might influence the efficacy of this control. In the present study, the relationships between the internal levels of protein (as measured by total body nitrogen) and carbohydrate (as measured by total body carbon) of the fruit fly Bactrocera tryoni (Froggatt) (Diptera: Tephritidae) are investigated, as well as its foraging behaviours in response to protein, fruit and cue-lure (a male-specific attractant) baits. Small cage behavioural experiments are conducted using flies from cultures of different nutritional status and wild flies sampled from the field during the fruiting cycle of a guava crop. For female flies, increasing total body nitrogen is correlated with decreased protein foraging and increased oviposition activity; increasing total body carbon levels generate the same behavioural changes except that the oviposition response is not significant. For males, there are no significant correlations between changes in total body nitrogen and total body carbon and protein or cue-lure foraging. For wild flies from the guava orchard, almost all of them are sexually mature when entering the crop and, over the entire season, total body nitrogen and total body carbon levels are such that protein hunger is unlikely for most flies. The results infer strongly that the requirements of wild, sexually mature flies for protein are minimal and that flies can readily gain sufficient nutrients from wild sources for their physiological needs. The results offer a mechanistic explanation for the poor response of male and mature female fruit flies to protein bait spray.

Secondary structure element alignment kernel method for prediction of protein structural classes

In this paper, we aim at predicting protein structural classes for low-homology data sets based on predicted secondary structures. We propose a new and simple kernel method, named as SSEAKSVM, to predict protein structural classes. The secondary structures of all protein sequences are obtained by using the tool PSIPRED and then a linear kernel on the basis of secondary structure element alignment scores is constructed for training a support vector machine classifier without parameter adjusting. Our method SSEAKSVM was evaluated on two low-homology datasets 25PDB and 1189 with sequence homology being 25% and 40%, respectively. The jackknife test is used to test and compare our method with other existing methods. The overall accuracies on these two data sets are 86.3% and 84.5%, respectively, which are higher than those obtained by other existing methods. Especially, our method achieves higher accuracies (88.1% and 88.5%) for differentiating the α + β class and the α/β class compared to other methods. This suggests that our method is valuable to predict protein structural classes particularly for low-homology protein sequences. The source code of the method in this paper can be downloaded at http://math.xtu.edu.cn/myphp/math/research/source/SSEAK_source_code.rar.

Cardiometabolic effects of genetic upregulation of the interleukin 1 receptor antagonist: A Mendelian randomisation analysis

Background To investigate potential cardiovascular and other effects of long-term pharmacological interleukin 1 (IL-1) inhibition, we studied genetic variants that produce inhibition of IL-1, a master regulator of inflammation. Methods We created a genetic score combining the effects of alleles of two common variants (rs6743376 and rs1542176) that are located upstream of IL1RN, the gene encoding the IL-1 receptor antagonist (IL-1Ra; an endogenous inhibitor of both IL-1α and IL-1β); both alleles increase soluble IL-1Ra protein concentration. We compared effects on inflammation biomarkers of this genetic score with those of anakinra, the recombinant form of IL-1Ra, which has previously been studied in randomised trials of rheumatoid arthritis and other inflammatory disorders. In primary analyses, we investigated the score in relation to rheumatoid arthritis and four cardiometabolic diseases (type 2 diabetes, coronary heart disease, ischaemic stroke, and abdominal aortic aneurysm; 453 411 total participants). In exploratory analyses, we studied the relation of the score to many disease traits and to 24 other disorders of proposed relevance to IL-1 signalling (746 171 total participants). Findings For each IL1RN minor allele inherited, serum concentrations of IL-1Ra increased by 0·22 SD (95% CI 0·18–0·25; 12·5%; p=9·3 × 10−33), concentrations of interleukin 6 decreased by 0·02 SD (−0·04 to −0·01; −1·7%; p=3·5 × 10−3), and concentrations of C-reactive protein decreased by 0·03 SD (−0·04 to −0·02; −3·4%; p=7·7 × 10−14). We noted the effects of the genetic score on these inflammation biomarkers to be directionally concordant with those of anakinra. The allele count of the genetic score had roughly log-linear, dose-dependent associations with both IL-1Ra concentration and risk of coronary heart disease. For people who carried four IL-1Ra-raising alleles, the odds ratio for coronary heart disease was 1·15 (1·08–1·22; p=1·8 × 10−6) compared with people who carried no IL-1Ra-raising alleles; the per-allele odds ratio for coronary heart disease was 1·03 (1·02–1·04; p=3·9 × 10−10). Per-allele odds ratios were 0·97 (0·95–0·99; p=9·9 × 10−4) for rheumatoid arthritis, 0·99 (0·97–1·01; p=0·47) for type 2 diabetes, 1·00 (0·98–1·02; p=0·92) for ischaemic stroke, and 1·08 (1·04–1·12; p=1·8 × 10−5) for abdominal aortic aneurysm. In exploratory analyses, we observed per-allele increases in concentrations of proatherogenic lipids, including LDL-cholesterol, but no clear evidence of association for blood pressure, glycaemic traits, or any of the 24 other disorders studied. Modelling suggested that the observed increase in LDL-cholesterol could account for about a third of the association observed between the genetic score and increased coronary risk. Interpretation Human genetic data suggest that long-term dual IL-1α/β inhibition could increase cardiovascular risk and, conversely, reduce the risk of development of rheumatoid arthritis. The cardiovascular risk might, in part, be mediated through an increase in proatherogenic lipid concentrations. Funding UK Medical Research Council, British Heart Foundation, UK National Institute for Health Research, National Institute for Health Research Cambridge Biomedical Research Centre, European Research Council, and European Commission Framework Programme 7.

The function of the RNA-binding protein TEL1 in moss reveals ancient regulatory mechanisms of shoot development

The shoot represents the basic body plan in land plants. It consists of a repeated structure composed of stems and leaves. Whereas vascular plants generate a shoot in their diploid phase, non-vascular plants such as mosses form a shoot (called the gametophore) in their haploid generation. The evolution of regulatory mechanisms or genetic networks used in the development of these two kinds of shoots is unclear. TERMINAL EAR1-like genes have been involved in diploid shoot development in vascular plants. Here, we show that disruption of PpTEL1 from the moss Physcomitrella patens, causes reduced protonema growth and gametophore initiation, as well as defects in gametophore development. Leafy shoots formed on ΔTEL1 mutants exhibit shorter stems with more leaves per shoot, suggesting an accelerated leaf initiation (shortened plastochron), a phenotype shared with the Poaceae vascular plants TE1 and PLA2/LHD2 mutants. Moreover, the positive correlation between plastochron length and leaf size observed in ΔTEL1 mutants suggests a conserved compensatory mechanism correlating leaf growth and leaf initiation rate that would minimize overall changes in plant biomass. The RNA-binding protein encoded by PpTEL1 contains two N-terminus RNA-recognition motifs, and a third C-terminus non-canonical RRM, specific to TEL proteins. Removal of the PpTEL1 C-terminus (including this third RRM) or only 16–18 amino acids within it seriously impairs PpTEL1 function, suggesting a critical role for this third RRM. These results show a conserved function of the RNA-binding PpTEL1 protein in the regulation of shoot development, from early ancestors to vascular plants, that depends on the third TEL-specific RRM.

Analysis on conservation of disulphide bonds and their structural features in homologous protein domain families

Background: Disulphide bridges are well known to play key roles in stability, folding and functions of proteins. Introduction or deletion of disulphides by site-directed mutagenesis have produced varying effects on stability and folding depending upon the protein and location of disulphide in the 3-D structure. Given the lack of complete understanding it is worthwhile to learn from an analysis of extent of conservation of disulphides in homologous proteins. We have also addressed the question of what structural interactions replaces a disulphide in a homologue in another homologue. Results: Using a dataset involving 34,752 pairwise comparisons of homologous protein domains corresponding to 300 protein domain families of known 3-D structures, we provide a comprehensive analysis of extent of conservation of disulphide bridges and their structural features. We report that only 54% of all the disulphide bonds compared between the homologous pairs are conserved, even if, a small fraction of the non-conserved disulphides do include cytoplasmic proteins. Also, only about one fourth of the distinct disulphides are conserved in all the members in protein families. We note that while conservation of disulphide is common in many families, disulphide bond mutations are quite prevalent. Interestingly, we note that there is no clear relationship between sequence identity between two homologous proteins and disulphide bond conservation. Our analysis on structural features at the sites where cysteines forming disulphide in one homologue are replaced by non-Cys residues show that the elimination of a disulphide in a homologue need not always result in stabilizing interactions between equivalent residues. Conclusion: We observe that in the homologous proteins, disulphide bonds are conserved only to a modest extent. Very interestingly, we note that extent of conservation of disulphide in homologous proteins is unrelated to the overall sequence identity between homologues. The non-conserved disulphides are often associated with variable structural features that were recruited to be associated with differentiation or specialisation of protein function.

Identification of actin as an estradiol 17-beta-stimulated protein in the human placenta

Addition of estradiol 17-beta to first trimester human placental minces resulted in an increased synthesis of a protein of apparent molecular weight 45 kDa. The specific involvement of estrogen in the stimulation of this protein was established by demonstrating a reduction in the level of this protein by the addition of CCS 16949 A, an inhibitor of aromatase, a key enzyme in the biosynthesis of estradiol 17-beta and ICI 182,780, an estrogen receptor antagonist. The protein was purified to homogeneity and N-terminal sequencing of two of the internal peptides obtained by enzymatic digestion of the protein, as well as the absence of a free N-terminal indicated that it could be actin. This was confirmed by Western blotting using commercially available actin antiserum. The role of estradiol 17-beta in the stimulation of actin synthesis in human placenta was also established by monitoring the quantitative inhibition of DNase I by actin.

The effect of pH on the unfolding pathway and stability of ribosome-inactivating protein abrin-II

The effect of pH on the unfolding pathway acid the stability of the toxic protein abrin-II have been studied by increasing denaturant concentrations of guanidine hydrochloride and by monitoring the change in 8,1-anilino naphthalene sulfonic acid (ANS) fluorescence upon binding to the hydrophobic sites of the protein. Intrinsic protein fluorescence, far and near UV-circular dichroism (CD) spectroscopy and ANS binding studies reveal that the unfolding of abrin-II occurs through two intermediates at pH 7.2 and one intermediate at pH 4.5. At pH 7.2, the two subunits A and B of abrin-II unfold sequentially. The native protein is more stable at pH 4.5 than at pH 7.2. However, the stability of the abrin-II A-subunit is not affected by a change in pH. These observations may assist in an understanding of the physiologically relevant transmembrane translocation of the toxin.

Thermodynamic characterization of the conformational stability of the homodimeric protein, pea lectin

The conformational stability of the homodimeric pea lectin was determined by both isothermal urea-induced and thermal denaturation in the absence and presence of urea. The denaturation profiles were analyzed to obtain the thermodynamic parameters associated with the unfolding of the protein. The data not only conform to the simple A(2) double left right arrow 2U model of unfolding but also are well described by the linear extrapolation model for the nature of denaturant-protein interactions. In addition, both the conformational stability (Delta G(s)) and the Delta C-p for the protein unfolding is quite high, at about 18.79 kcal/ mol and 5.32 kcal/(mol K), respectively, which may be a reflection of the relatively larger size of the dimeric molecule (M-r 49 000) and, perhaps, a consequent larger buried hydrophobic core in the folded protein. The simple two-state (A(2) double left right arrow 2U) nature of the unfolding process, with the absence of any monomeric intermediate, suggests that the quaternary interactions alone may contribute significantly to the conformational stability of the oligomer-a point that may be general to many oligomeric proteins.

Fluctuations in protein synthesis from a single RNA template: Stochastic kinetics of ribosomes

Proteins are polymerized by cyclic machines called ribosomes, which use their messenger RNA (mRNA) track also as the corresponding template, and the process is called translation. We explore, in depth and detail, the stochastic nature of the translation. We compute various distributions associated with the translation process; one of them-namely, the dwell time distribution-has been measured in recent single-ribosome experiments. The form of the distribution, which fits best with our simulation data, is consistent with that extracted from the experimental data. For our computations, we use a model that captures both the mechanochemistry of each individual ribosome and their steric interactions. We also demonstrate the effects of the sequence inhomogeneities of real genes on the fluctuations and noise in translation. Finally, inspired by recent advances in the experimental techniques of manipulating single ribosomes, we make theoretical predictions on the force-velocity relation for individual ribosomes. In principle, all our predictions can be tested by carrying out in vitro experiments.

The critical role of water in spider silk and its consequence for protein mechanics

Due to its remarkable mechanical and biological properties, there is considerable interest in understanding, and replicating, spider silk's stress-processing mechanisms and structure-function relationships. Here, we investigate the role of water in the nanoscale mechanics of the different regions in the spider silk fibre, and their relative contributions to stress processing. We propose that the inner core region, rich in spidroin II, retains water due to its inherent disorder, thereby providing a mechanism to dissipate energy as it breaks a sacrificial amide-water bond and gains order under strain, forming a stronger amide-amide bond. The spidroin I-rich outer core is more ordered under ambient conditions and is inherently stiffer and stronger, yet does not on its own provide high toughness. The markedly different interactions of the two proteins with water, and their distribution across the fibre, produce a stiffness differential and provide a balance between stiffness, strength and toughness under ambient conditions. Under wet conditions, this balance is destroyed as the stiff outer core material reverts to the behaviour of the inner core.

Immunogenic and protective properties of haemagglutinin protein (H) of Rinderpest virus expressed by a recombinant baculovirus

The hemagglutinin (H) protein of Rinderpest virus expressed by a recombinant buculovirus used as a vaccine produced high titres of neutralizing antibody to Rinderpest virus in the vaccinated cattle, comparable to the levels produced by live attenuated vaccine. The immunized cattle were protected against a vaccine-virus challenge, as demonstrated by the failure of development of antibodies to N protein of the vaccine virus. The lack of replication of vaccine virus in the immunized cattle indicated that they are capable of showing a protective response if challenged with a virulent virus.

Photoinduced DNA and Protein Cleavage Activity of Ferrocene-Conjugated Ternary Copper(II) Complexes

Ferrocene-conjugated ternary copper(II) complexes [Cu(L)(B)](ClO4)(2), where L is FcCH(2)N(CH2Py)(2) (Fc = (eta(5)-C5H4)Fe-II(eta(5)-C5H5)) and B is a phenanthroline base, viz., 2,2'-bipyridine (bpy, 1), 1, 10-phenanthroline (phen, 2), dipyrido[3,2-d:2',3'-f]quinoxaline (dpq, 3), and dipyrido[3,2-a:2',3'-c]phenazine (dppz, 4), have been synthesized and characterized by various spectroscopic and analytical techniques. The bpy complex 1, as its hexafluorophosphate salt, has been structurally characterized by X-ray crystallography. The molecular structure shows the copper(II) center having an essentially square-pyramidal coordination geometry in which L with a pendant ferrocenyl (Fc) moiety and bpy show respective tridentate and bidentate modes of binding to the metal center. The complexes are redox active, showing a reversible cyclic voltammetric response of the Fc(+)-Fc couple near 0.5 V vs SCE and a quasi-reversible Cu(II)-Cu(I) couple near 0.0 V. Complexes 2-4 show binding affinity to calf thymus (CT) DNA, giving binding constant (K-b) values in the range of 4.2 x 10(4) to 2.5 x 10(5) M-1. Thermal denaturation and viscometric titration data suggest groove binding and/or a partial intercalative mode of binding of the complexes to CT DNA. The complexes show good binding propensity to the bovine serum albumin (BSA) protein, giving K-BSA values of similar to 10(4) M-1 for the bpy and phen complexes and similar to 10(5) M-1 for the dpq and dppz complexes. Complexes 2-4 exhibit efficient chemical nuclease activity in the presence of 3-mercapto-propionic acid (MPA) as a reducing agent or hydrogen peroxide (H2O2) as an oxidizing agent. Mechanistic studies reveal formation of hydroxyl radicals as the reactive species. The dpq and dppz complexes are active in cleaving supercoiled (SC) pUC19 DNA on photoexposure to visible light of different wavelengths including red light using an argon-krypton mixed gas ion laser. Mechanistic investigations using various inhibitors reveal the fort-nation of hydroxyl radicals in the DNA photocleavage reactions. The dppz complex 4, which shows efficient photoioduced BSA cleavage activity, is a potent multifunctional model nuclease and protease in the chemistry of photodynamic therapy (PDT) of cancer.

Mutations in STIL, Encoding a Pericentriolar and Centrosomal Protein, Cause Primary Microcephaly

Primary microcephaly (MCPH) is an autosomal-recessive congenital disorder characterized by smaller-than-normal brain size and mental retardation. MCPH is genetically heterogeneous with six known loci: MCPH1-MCPH6. We report mapping of a novel locus, MCPH7, to chromosome 1p32.3-p33 between markers D1S2797 and D1S417, corresponding to a physical distance of 8.39 Mb. Heterogeneity analysis of 24 families previously excluded from linkage to the six known MCPH loci suggested linkage of five families (20.83%) to the MCPH7 locus. In addition, four families were excluded from linkage to the MCPH7 locus as well as all of the six previously known loci, whereas the remaining 15 families could not be conclusively excluded or included. The combined maximum two-point LOD score for the linked families was 5.96 at marker D1S386 at theta = 0.0. The combined multipoint LOD score was 6.97 between markers D1S2797 and D1S417. Previously, mutations in four genes, MCPH1, CDK5RAP2, ASPM, and CENPJ, that code for centrosomal proteins have been shown to cause this disorder. Three different homozygous mutations in STIL, which codes for a pericentriolar and centrosomal protein, were identified in patients from three of the five families linked to the MCPH7 locus; all are predicted to truncate the STIL protein. Further, another recently ascertained family was homozygous for the same mutation as one of the original families. There was no evidence for a common haplotype. These results suggest that the centrosome and its associated structures are important in the control of neurogenesis in the developing human brain.

A soluble preparation from developing groundnut seeds (Arachis hypogaea) catalyzes de novo synthesis of long chain fatty acids

A 100,000 x g supernatant fraction prepared from developing groundnut seeds (30-35 days after flowering) catalyzed the synthesis of fatty acids from [l-14C]acetate at a rate of 120nmoles of acetate incorporated per hr per gram fresh weight of tissue. 90% of this incorporated label was associated with fatty acids. The major fatty acids formed were stearic- (77%) and palmitic acids (14%) with 4% of oleic acid. The fatty acid synthetase activity was stable when stored at 0-4 degrees C for at least fifteen days. It is concluded from these results that acetyl-coA carboxylase and all the enzymes of fatty acid synthetase from developing groundnut seeds are soluble.