Calculations of the free energy of interaction of the c-Fos−c-Jun coiled coil: Effects of the solvation model and the inclusion of polarization effects

The leucine zipper region of activator protein-1 (AP-1) comprises the c-Jun and c-Fos proteins and constitutes a well-known coiled coil protein−protein interaction motif. We have used molecular dynamics (MD) simulations in conjunction with the molecular mechanics/Poisson−Boltzmann generalized-Born surface area [MM/PB(GB)SA] methods to predict the free energy of interaction of these proteins. In particular, the influence of the choice of solvation model, protein force field, and water potential on the stability and dynamic properties of the c-Fos−c-Jun complex were investigated. Use of the AMBER polarizable force field ff02 in combination with the polarizable POL3 water potential was found to result in increased stability of the c-Fos−c-Jun complex. MM/PB(GB)SA calculations revealed that MD simulations using the POL3 water potential give the lowest predicted free energies of interaction compared to other nonpolarizable water potentials. In addition, the calculated absolute free energy of binding was predicted to be closest to the experimental value using the MM/GBSA method with independent MD simulation trajectories using the POL3 water potential and the polarizable ff02 force field, while all other binding affinities were overestimated.

Use of protein A gene fusions for the analysis of structure-function relationship of the transactivator protein C of bacteriophage Mu

A sensitive dimerization assay for DNA binding proteins has been developed using gene fusion technology. For this purpose, we have engineered a gene fusion using protein A gene of Staphylococcus aureus and C gene, the late gene transactivator of bacteriophage Mu. The C gene was fused to the 3' end of the gene for protein A to generate an A- C fusion. The overexpressed fusion protein was purified in a single step using immunoglobulin affinity chromatography. Purified fusion protein exhibits DNA binding activity as demonstrated by electrophoretic mobility shift assays. When the fusion protein A-C was mixed with C and analyzed for DNA binding, in addition to C and A-C specific complexes, a single intermediate complex comprising of a heterodimer of C and A-C fusion proteins was observed. Further, the protein A moiety in the fusion protein A-C does not contribute to DNA binding as demonstrated by proteolytic cleavage and circular dichroism (CD) analysis. The assay has also been applied to analyze the DNA binding domain of C protein by generating fusions between protein A and N- and C-terminal deletion mutants of C. The results indicate a role for the region towards the carboxy terminal of the protein in DNA binding. The general applicability of this method is discussed.

Observation of a mixed antiparallel and parallel beta-sheet motif in the crystal structure of Boc-Ala-Ile-Aib-OMe

A diastereomeric mixture of the tripeptide Boc-Ala-Ile-Aib-OMe crystallized in the space group P1 from CH3OH/H2O. The unit cell parameters are a = 10.593(2) A, b = 14.377(3) A, c = 17.872(4) A, alpha = 104.41(2) degrees, beta = 90.55(2) degrees, gamma = 106.91(2) degrees, V = 2512.4 A3, Z = 4. X-Ray crystallographic studies show the presence of four molecules in the asymmetric unit consisting of two pairs of diastereomeric peptides, Boc-L-Ala-L-Ile-Aib-OMe and Boc-L-Ala-D-Ile-Aib-OMe. The four molecules in the asymmetric unit form a rarely found mixed antiparallel and parallel beta-sheet hydrogen bond motif. The Ala and (L,D)-Ile residues in all the four molecules adopt the extended conformations, while the phi, psi values of the Aib residues are in the right-handed helical region. In one of the molecules the Ile sidechain adopts the unusual gauche conformation about the C beta-C gamma bond.

NMR and X-ray crystallographic studies on cyclic tetrapeptide, cyclo (D-Phe-Pro-Sar-Gly)

The conformation of the synthetic cyclic tetrapeptide cyclo(D-Phe-Pro-Sar-Gly) has been determined in solution using the nuclear magnetic resonance technique and in the crystal state by X-ray crystallography. Results showed that the peptide exhibited two different conformations in solution, conformer 1 having cis-trans-cis-trans peptide bonds and conformer 2 having trans-cis-trans-cis peptide bonds. No intramolecular hydrogen bonds were observed in the structures. The X-ray diffraction studies showed the crystals to be orthorhombic with space group P2(1)2(1)2(1) with unit-cell dimensions, a = 5.790, b = 10.344, c = 31.446 A, Z = 4, R = 0.104 for 2301 observed reflections. The crystal structure showed only one type of conformer having cis-trans-cis-trans peptide bonds similar to the conformer 1 in solution.

Theoretical Studies on Coupled Electron and Proton Transfer in Cytochrome c Oxidase

The complexity of life is based on an effective energy transduction machinery, which has evolved during the last 3.5 billion years. In aerobic life, the utilization of the high oxidizing potential of molecular oxygen powers this machinery. Oxygen is safely reduced by a membrane bound enzyme, cytochrome c oxidase (CcO), to produce an electrochemical proton gradient over the mitochondrial or bacterial membrane. This gradient is used for energy-requiring reactions such as synthesis of ATP by F0F1-ATPase and active transport. In this thesis, the molecular mechanism by which CcO couples the oxygen reduction chemistry to proton-pumping has been studied by theoretical computer simulations. By building both classical and quantum mechanical model systems based on the X-ray structure of CcO from Bos taurus, the dynamics and energetics of the system were studied in different intermediate states of the enzyme. As a result of this work, a mechanism was suggested by which CcO can prevent protons from leaking backwards in proton-pumping. The use and activation of two proton conducting channels were also enlightened together with a mechanism by which CcO sorts the chemical protons from pumped protons. The latter problem is referred to as the gating mechanism of CcO, and has remained a challenge in the bioenergetics field for more than three decades. Furthermore, a new method for deriving charge parameters for classical simulations of complex metalloenzymes was developed.

Diversity and Distribution of Conidae from the TamilNadu Coast of India (Mollusca: Caenogastropoda: Conidae)

A survey of the marine gastropod genus Conus Linnaeus was conducted along the TamilNadu Coast of India to explore the regional geographic distribution and diversity. The 60 species observed increased the number of Indian Conidae from 77 to 81. Conus imperialis Linne, C. mitratus Hwass in Bruguiere, C. striolatus Kiener and C. violaceus Gmelin are newly recorded from the study area. Conus amadis Gmelin was the most widely distributed species. The highest diversity (48 species) occurred in the Gulf of Mannar, followed by 22 species from northern, six from southern, and five from the Palk Bay regions. We suggest that the rich diversity recorded in the Gulf of Mannar reflects the physical conditions, microhabitats and required resources such as food and shelter that favour the occurrence of the large number of Conus species.

Revision of the family Alycidae (Acariformes, Acari), with special reference to European species

The thesis provides a proposal to divide Alycidae G. Canestrini & Fanzago into two subfamilies and four tribes. This new hierarchy is based on a reassessment and reranking of new and previously known synapomorphies of the clusters concerned by cladistic analysis, using 60 morphological characters for 48 ingroup species. The basic characters of the taxa are illustrated either by SEM micrographs (Scanning Electron Microscopy) or by outline drawings. The presented classification includes the definitions of Alycini G. Canestrini & Fanzago new rank; Bimichaeliini Womersley new rank; Petralycini new rank; and the (re)descriptions of Alycus C.L. Koch, Pachygnathus Dugès, Amphialycus Zachvatkin, Bimichaelia Thor and Laminamichaelia gen. nov. The species described or redescribed are: Pachygnathus wasastjernae sp. nov. from Kvarken (Merenkurkku), Finland; Pachygnathus villosus Dugès (in Oken); Alycus roseus C.L. Koch; Alycus denasutus (Grandjean) comb. and stat. nov.; Alycus trichotus (Grandjean) comb. nov.; Alycus marinus (Schuster) comb. nov.; Amphialycus (Amphialycus) pentophthalmus Zachvatkin; Amphialycus (Amphialycus) leucogaster (Grandjean); and Amphialycus (Orthacarus) oblongus (Halbert) comb. nov.; Bimichaelia augustana (Berlese); Bimichaelia sarekensis Trägårdh; Laminamichaelia setigera (Berlese) comb. nov.; Laminamichelia arbusculosa (Grandjean) comb. nov.; Laminamichelia subnuda (Berlese) comb. nov. and Petralycus unicornis Grandjean. Fourteen nominal species were found to be junior synonymies. The importance of sensory organs in taxonomy is well recognized, but inclusion of the elaborate skin pattern seemed to improve essentially the usefulness of the prodorsal sensory area. The detailed pictures of the prodorsa of the European alycids could be used like passport photographs for the species. A database like this of prodorsa of other mite taxa as well might be an answer to future needs of species identification in soil zoology, ecology and conservation.

The C-terminal domain is sufficient for endonuclease activity of Neisseria gonorrhoeae MutL

The mutL gene of Neisseria gonorrhoeae has been cloned and the gene product purified. We have found that the homodimeric N. gonorrhoeae MutL (NgoL) protein displays an endonuclease activity that incises covalently closed circular DNA in the presence of Mn2+, Mg2+ or Ca2+ ions, unlike human MutL alpha which shows endonuclease activity only in the presence of Mn2+. We report in the present paper that the C-terminal domain of N. gonorrhoeae MutL (NgoL-CTD) consisting of amino acids 460-658 exhibits Mn2+-dependent endonuclease activity. Sedimentation velocity, sedimentation equilibrium and dynamic light scattering experiments show NgoL-CTD to be a dimer. The probable endonucleolytic active site is localized to a metal-binding motif, DMHAX(2)EX(4)E, and the nicking endonuclease activity is dependent on the integrity of this motif. By in vitro comparison of wild-type and it mutant NgoL-CTD protein, we show that the latter protein exhibits highly reduced endonuclease activity. We therefore suggest that the mode of excision initiation in DNA mismatch repair may be different in organisms that lack MutH protein, but have MutL proteins that harbour the D[M/Q]HAX(2)EX(4)E motif.

A semi-analytical particle filter for identification of nonlinear oscillators

Particle filters find important applications in the problems of state and parameter estimations of dynamical systems of engineering interest. Since a typical filtering algorithm involves Monte Carlo simulations of the process equations, sample variance of the estimator is inversely proportional to the number of particles. The sample variance may be reduced if one uses a Rao-Blackwell marginalization of states and performs analytical computations as much as possible. In this work, we propose a semi-analytical particle filter, requiring no Rao-Blackwell marginalization, for state and parameter estimations of nonlinear dynamical systems with additively Gaussian process/observation noises. Through local linearizations of the nonlinear drift fields in the process/observation equations via explicit Ito-Taylor expansions, the given nonlinear system is transformed into an ensemble of locally linearized systems. Using the most recent observation, conditionally Gaussian posterior density functions of the linearized systems are analytically obtained through the Kalman filter. This information is further exploited within the particle filter algorithm for obtaining samples from the optimal posterior density of the states. The potential of the method in state/parameter estimations is demonstrated through numerical illustrations for a few nonlinear oscillators. The proposed filter is found to yield estimates with reduced sample variance and improved accuracy vis-a-vis results from a form of sequential importance sampling filter.

DNA- and chromatin-condensing properties of rat testes H1a and H1t compared to those of rat liver H1bdec: H1t is a poor condenser of chromatin

Histones H1a and H1t are two major linker histone variants present at the pachytene interval of mammalian spermatogenesis. The DNA- and chromatin-condensing properties of these two variants isolated from rat testes were studied and compared with those from rat liver. For this purpose, the histone H1 subtypes were purified from the respective tissues using bath acid and salt extraction procedures, Circular dichroism studies revealed that acid exposure during isolation affects the alpha-helical structure of both the globular domain (in the presence of 1 M NaCl) and the C-terminal lambda-tail (in the presence of 60% trifluoroethanol). The condensation of rat oligonucleosomal DNA, as measured by circular dichroism spectroscopy, by the salt-extracted histone H1 was at least 10 times more efficient than condensation by the acid-extracted histone H1. A site size of 16-20 base pairs was calculated for the salt-extracted histone H1. Among the different histone H1 subtypes, somatic histone H1bdec had the highest DNA-condensing property, followed by histone H1a and histone H1t. All the salt-extracted histones condensed rat oligonucleosomal DNA more efficiently than linear pBR-322 DNA, Histones H1bdec and H1a condensed histone H1-depleted chromatin, prepared from rat liver nuclei, with relatively equal efficiency. On the other hand, there was no condensation of histone H1-depleted chromatin with the testes specific histone H1t. A comparison of the amino acid sequences of histone H1d (rat) and histone H1t (rat) revealed several interesting differences in the occurrence of DNA-binding motifs at the C-terminus. A striking observation is the presence of a direct repeat of an octapeptide motif K(A)T(S)PKKA(S)K(T)K(A) in histone H1d that is absent in histone H1t.

DNA-condensing and chromatin-condensing properties of rat testes hla and hit compared to those of rat-liver hlbdec - hlt is a poor condenser of chromatin

Histones H1a and H1t are two major linker histone variants present at the pachytene interval of mammalian spermatogenesis. The DNA- and chromatin-condensing properties of these two variants isolated from rat testes were studied and compared with those from rat liver. For this purpose, the histone H1 subtypes were purified from the respective tissues using bath acid and salt extraction procedures, Circular dichroism studies revealed that acid exposure during isolation affects the alpha-helical structure of both the globular domain (in the presence of 1 M NaCl) and the C-terminal lambda-tail (in the presence of 60% trifluoroethanol). The condensation of rat oligonucleosomal DNA, as measured by circular dichroism spectroscopy, by the salt-extracted histone H1 was at least 10 times more efficient than condensation by the acid-extracted histone H1. A site size of 16-20 base pairs was calculated for the salt-extracted histone H1. Among the different histone H1 subtypes, somatic histone H1bdec had the highest DNA-condensing property, followed by histone H1a and histone H1t. All the salt-extracted histones condensed rat oligonucleosomal DNA more efficiently than linear pBR-322 DNA, Histones H1bdec and H1a condensed histone H1-depleted chromatin, prepared from rat liver nuclei, with relatively equal efficiency. On the other hand, there was no condensation of histone H1-depleted chromatin with the testes specific histone H1t. A comparison of the amino acid sequences of histone H1d (rat) and histone H1t (rat) revealed several interesting differences in the occurrence of DNA-binding motifs at the C-terminus. A striking observation is the presence of a direct repeat of an octapeptide motif K(A)T(S)PKKA(S)K(T)K(A) in histone H1d that is absent in histone H1t.

Molecular epidemiology of human rhinoviruses

The first part of this work investigates the molecular epidemiology of a human enterovirus (HEV), echovirus 30 (E-30). This project is part of a series of studies performed in our research team analyzing the molecular epidemiology of HEV-B viruses. A total of 129 virus strains had been isolated in different parts of Europe. The sequence analysis was performed in three different genomic regions: 420 nucleotides (nt) in the VP4/VP2 capsid protein coding region, the entire VP1 capsid protein coding gene of 876 nt, and 150 nt in the VP1/2A junction region. The analysis revealed a succession of dominant sublineages within a major genotype. The temporally earlier genotypes had been replaced by a genetically homogenous lineage that has been circulating in Europe since the late 1970s. The same genotype was found by other research groups in North America and Australia. Globally, other cocirculating genetic lineages also exist. The prevalence of a dominant genotype makes E-30 different from other previously studied HEVs, such as polioviruses and coxsackieviruses B4 and B5, for which several coexisting genetic lineages have been reported. The second part of this work deals with molecular epidemiology of human rhinoviruses (HRVs). A total of 61 field isolates were studied in the 420-nt stretch in the capsid coding region of VP4/VP2. The isolates were collected from children under two years of age in Tampere, Finland. Sequences from the clinical isolates clustered in the two previously known phylogenetic clades. Seasonal clustering was found. Also, several distinct serotype-like clusters were found to co-circulate during the same epidemic season. Reappearance of a cluster after disappearing for a season was observed. The molecular epidemiology of the analyzed strains turned out to be complex, and we decided to continue our studies of HRV. Only five previously published complete genome sequences of HRV prototype strains were available for analysis. Therefore, all designated HRV prototype strains (n=102) were sequenced in the VP4/VP2 region, and the possibility of genetic typing of HRV was evaluated. Seventy-six of the 102 prototype strains clustered in HRV genetic group A (HRV-A) and 25 in group B (HRV-B). Serotype 87 clustered separately from other HRVs with HEV species D. The field strains of HRV represented as many as 19 different genotypes, as judged with an approximate demarcation of a 20% nt difference in the VP4/VP2 region. The interserotypic differences of HRV were generally similar to those reported between different HEV serotypes (i.e. about 20%), but smaller differences, less than 10%, were also observed. Because some HRV serotypes are genetically so closely related, we suggest that the genetic typing be performed using the criterion "the closest prototype strain". This study is the first systematic genetic characterization of all known HRV prototype strains, providing a further taxonomic proposal for classification of HRV. We proposed to divide the genus Human rhinoviruses into HRV-A and HRV-B. The final part of the work comprises a phylogenetic analysis of a subset (48) of HRV prototype strains and field isolates (12) in the nonstructural part of the genome coding for the RNA-dependent RNA polymerase (3D). The proposed division of the HRV strains in the species HRV-A and HRV-B was also supported by 3D region. HRV-B clustered closer to HEV species B, C, and also to polioviruses than to HRV-A. Intraspecies variation within both HRV-A and HRV-B was greater in the 3D coding region than in the VP4/VP2 coding region, in contrast to HEV. Moreover, the diversity of HRV in 3D exceeded that of HEV. One group of HRV-A, designated HRV-A', formed a separate cluster outside other HRV-A in the 3D region. It formed a cluster also in the capsid region, but located within HRV-A. This may reflect a different evolutionary history of distinct genomic regions among HRV-A. Furthermore, the tree topology within HRV-A in the 3D region differed from that in the VP4/VP2, suggesting possible recombination events in the evolution of the strains. No conflicting phylogenies were observed in any of the 12 field isolates. Possible recombination was further studied using the Similarity and Bootscanning analyses of the complete genome sequences of HRV available in public databases. Evidence for recombination among HRV-A was found, as HRV2 and HRV39 showed higher similarity in the nonstructural part of the genome. Whether HRV2 and HRV39 strains - and perhaps also some other HRV-A strains not yet completely sequenced - are recombinants remains to be determined.

How many inpatients in our hospitals have foot complications? The foot disease in inpatients study

Background Foot complications have been found to affect large proportions of hospital in patients with diabetes. However, no studies have investigated the proportion of foot complications affecting all people in general inpatient populations. The aims of this cross-sectional study were to investigate the point-prevalence of different foot complications in general inpatient populations, analyse differences in diabetes and non-diabetes sub-groups, and examine characteristics of people primarily admitted for a foot complication. Methods Eligible participants were all adults admitted overnight, for any reason, into five diverse hospitals on one day; excluding maternity, mental health and cognitively impaired patients. All participants underwent a physical foot examination, by trained podiatrists using validated measures, to clinically diagnose different foot complications; including foot wounds, infections, deformity, peripheral arterial disease (PAD) and peripheral neuropathy (PN). Data were also collected on participants' primary reason for admission and a range of demographic, social determinant, medical history, foot complication history, self-care and footwear risk factors. Results Overall, 733 participants consented (83% of eligible participants); mean(±SD) age 62(±19) years, 480 (55.8%) male and 172 (23.5%) had diabetes. Foot complication prevalence included: wounds 9.0% (95% CI) (5.1-8.7), infections 3.3% (2.2-4.9), deformity 22.4% (19.5-26.7), PAD 21.0% (18.2-24.1) and PN 22.0% (19.1-25.1). Diabetes populations had significantly more foot complications than non-diabetes (p < 0.01); wounds (15.7% vs 7.0%), infections (7.1% vs 2.2%), deformity (30.5% vs 19.9%), PAD (35.1% vs 16.7%) and PN (43.3% vs 15.4%). Foot complications were the primary reason for admission in 7.4% (95% CI) (5.7-9.5) of all participants. In a backwards stepwise multivariate analysis having a foot complication as the primary reason for admission was independently associated (OR (95% CI) with foot wounds (18.9 (7.3-48.7)), foot infections (6.0 (1.6-22.4)), history of amputation (4.7 (1.3-17.0) and PAD (2.9 (1.3-6.6)). Conclusions Findings of this study indicate one in every ten hospital inpatients had an active foot wound or infection. In patients with diabetes had significantly higher proportions of foot complications than non-diabetes inpatients. Remarkably one in every thirteen inpatients in this study were primarily hospitalised for a foot complication. Further research and policy is required to tackle this seemingly large inpatient foot complication burden.

Mycobacterium tuberculosis FtsZ requires at least one arginine residue at the C-terminal end for polymerization in vitro

We examined whether C-terminal residues of soluble recombinant FtsZ of Mycobacterium tuberculosis (MtFtsZ) have any role in MtFtsZ polymerization in vitro. MtFtsZ-delta C1, which lacks C-terminal extreme Arg residue (underlined in the C-terminal extreme stretch of 13 residues, DDDDVDVPPFMRR), but retaining the penultimate Arg residue (DDDDVDVPPFMR), polymerizes like full-length MtFtsZ in vitro. However, MtFtsZ-delta C2 that lacks both the Arg residues at the C-terminus (DDDDVDVPPFM), neither polymerizes at pH 6.5 nor forms even single- or double-stranded filaments at pH 7.7 in the presence of 10 mM CaCl2. Neither replacement of the penultimate Arg residue, in the C-terminal Arg deletion mutant DDDDVDVPPFMR, with Lys or His or Ala or Asp (DDDDVDVPPFMK/H/A/D) enabled polymerization. Although MtFtsZ-delta C2 showed secondary and tertiary structural changes, which might have affected polymerization, GTPase activity of MtFtsZ-delta C2 was comparable to that of MtFtsZ. These data suggest that MtFtsZ requires an Arg residue as the extreme C-terminal residue for polymerization in vitro. The polypeptide segment containing C-terminal 67 residues, whose coordinates were absent from MtFtsZ crystal structure, was modeled on tubulin and MtFtsZ dimers. Possibilities for the influence of the C-terminal Arg residues on the stability of the dimer and thereby on MtFtsZ polymerization have been discussed.

Recognition of Interaction Interface Residues in Low-Resolution Structures of Protein Assemblies Solely from the Positions of C alpha Atoms

Background: The number of available structures of large multi-protein assemblies is quite small. Such structures provide phenomenal insights on the organization, mechanism of formation and functional properties of the assembly. Hence detailed analysis of such structures is highly rewarding. However, the common problem in such analyses is the low resolution of these structures. In the recent times a number of attempts that combine low resolution cryo-EM data with higher resolution structures determined using X-ray analysis or NMR or generated using comparative modeling have been reported. Even in such attempts the best result one arrives at is the very course idea about the assembly structure in terms of trace of the C alpha atoms which are modeled with modest accuracy. Methodology/Principal Findings: In this paper first we present an objective approach to identify potentially solvent exposed and buried residues solely from the position of C alpha atoms and amino acid sequence using residue type-dependent thresholds for accessible surface areas of C alpha. We extend the method further to recognize potential protein-protein interface residues. Conclusion/Significance: Our approach to identify buried and exposed residues solely from the positions of C alpha atoms resulted in an accuracy of 84%, sensitivity of 83-89% and specificity of 67-94% while recognition of interfacial residues corresponded to an accuracy of 94%, sensitivity of 70-96% and specificity of 58-94%. Interestingly, detailed analysis of cases of mismatch between recognition of interface residues from C alpha positions and all-atom models suggested that, recognition of interfacial residues using C alpha atoms only correspond better with intuitive notion of what is an interfacial residue. Our method should be useful in the objective analysis of structures of protein assemblies when positions of only C alpha positions are available as, for example, in the cases of integration of cryo-EM data and high resolution structures of the components of the assembly.