Solution structures of Mycobacterium tuberculosis thioredoxin C and models of intact thioredoxin system suggest new approaches to inhibitor and drug design

Here, we report the NMR solution structures of Mycobacterium tuberculosis (M. tuberculosis) thioredoxin C in both oxidized and reduced states, with discussion of structural changes that occur in going between redox states. The NMR solution structure of the oxidized TrxC corresponds closely to that of the crystal structure, except in the C‐terminal region. It appears that crystal packing effects have caused an artifactual shift in the α4 helix in the previously reported crystal structure, compared with the solution structure. On the basis of these TrxC structures, chemical shift mapping, a previously reported crystal structure of the M. tuberculosis thioredoxin reductase (not bound to a Trx) and structures for intermediates in the E. coli thioredoxin catalytic cycle, we have modeled the complete M. tuberculosis thioredoxin system for the various steps in the catalytic cycle. These structures and models reveal pockets at the TrxR/TrxC interface in various steps in the catalytic cycle, which can be targeted in the design of uncompetitive inhibitors as potential anti‐mycobacterial agents, or as chemical genetic probes of function. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Impact of deglycosylation and thermal stress on conformational stability of a full length murine igG2a monoclonal antibody: Observations from molecular dynamics simulations

With the rise of antibody based therapeutics as successful medicines, there is an emerging need to understand the fundamental antibody conformational dynamics and its implications towards stability of these medicines. Both deglycosylation and thermal stress have been shown to cause conformational destabilization and aggregation in monoclonal antibodies. Here, we study instabilities caused by deglycosylation and by elevated temperature (400 K) by performing molecular dynamic simulations on a full length murine IgG2a mAb whose crystal structure is available in the Protein Data bank. C^α‐atom root mean square deviation and backbone root mean square fluctuation calculations show that deglycosylation perturbs quaternary and tertiary structures in the C_H2 domains. In contrast, thermal stress pervades throughout the antibody structure and both Fabs and Fc regions are destabilized. The thermal stress applied in this study was not sufficient to cause large scale unfolding within the simulation time and most amino acid residues showed similar average solvent accessible surface area and secondary structural conformations in all trajectories. C_H3 domains were the most successful at resisting the conformational destabilization. The simulations helped identify aggregation prone regions, which may initiate cross‐β motif formation upon deglycosylation and upon applying thermal stress. Deglycosylation leads to increased backbone fluctuations and solvent exposure of a highly conserved APR located in the edge β‐strand A of the C_H2 domains. Aggregation upon thermal stress is most likely initiated by two APRs that overlap with the complementarity determining regions. This study has important implications for rational design of antibody based therapeutics that are resistant towards aggregation. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Prediction of protein domain boundaries from inverse covariances

It has been known even since relatively few structures had been solved that longer protein chains often contain multiple domains, which may fold separately and play the role of reusable functional modules found in many contexts. In many structural biology tasks, in particular structure prediction, it is of great use to be able to identify domains within the structure and analyze these regions separately. However, when using sequence data alone this task has proven exceptionally difficult, with relatively little improvement over the naive method of choosing boundaries based on size distributions of observed domains. The recent significant improvement in contact prediction provides a new source of information for domain prediction. We test several methods for using this information including a kernel smoothing‐based approach and methods based on building alpha‐carbon models and compare performance with a length‐based predictor, a homology search method and four published sequence‐based predictors: DOMCUT, DomPRO, DLP‐SVM, and SCOOBY‐DOmain. We show that the kernel‐smoothing method is significantly better than the other ab initio predictors when both single‐domain and multidomain targets are considered and is not significantly different to the homology‐based method. Considering only multidomain targets the kernel‐smoothing method outperforms all of the published methods except DLP‐SVM. The kernel smoothing method therefore represents a potentially useful improvement to ab initio domain prediction. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

A probabilistic model for secondary structure prediction from protein chemical shifts

Protein chemical shifts encode detailed structural information that is difficult and computationally costly to describe at a fundamental level. Statistical and machine learning approaches have been used to infer correlations between chemical shifts and secondary structure from experimental chemical shifts. These methods range from simple statistics such as the chemical shift index to complex methods using neural networks. Notwithstanding their higher accuracy, more complex approaches tend to obscure the relationship between secondary structure and chemical shift and often involve many parameters that need to be trained. We present hidden Markov models (HMMs) with Gaussian emission probabilities to model the dependence between protein chemical shifts and secondary structure. The continuous emission probabilities are modeled as conditional probabilities for a given amino acid and secondary structure type. Using these distributions as outputs of first‐ and second‐order HMMs, we achieve a prediction accuracy of 82.3%, which is competitive with existing methods for predicting secondary structure from protein chemical shifts. Incorporation of sequence‐based secondary structure prediction into our HMM improves the prediction accuracy to 84.0%. Our findings suggest that an HMM with correlated Gaussian distributions conditioned on the secondary structure provides an adequate generative model of chemical shifts. Proteins 2013; © 2012 Wiley Periodicals, Inc.     

A note on pterosaur nesting behavior

Based on examination of eggshell structure and predicted vapor conductances in eggshells in recently described material from Argentina and China we conclude that pterosaurs buried their eggs. Egg-burying imposes theoretical restrictions on the distribution of pterosaurs, both geographically and spatially, raises the possibility of thermal sex determination and supports previous suggestions that they exhibited nesting fidelity. Some features associated with egg-burying, such as weight savings, are likely to have been fortuitous pre-adaptations for these flying reptiles, but others may have disadvantaged them relative to avian competitors or increased their vulnerability to extinction in a cooling climate.     

Anti-barnacle activity of novel simple alkyl isocyanides derived from citronellol

Twenty novel simple alkyl isocyanides derived from citronellol were synthesized and evaluated for their antifouling activity and toxicity against cypris larvae of the barnacle, Balanus amphitrite. The anti-barnacle activity of the synthesized isocyanides was in the EC₅₀ range of 0.08–1.49 μg ml^?1. Simple isocyanides containing a benzoate and chloro group showed the most potent anti-barnacle activity. In addition, none of the synthesized compounds showed significant toxicity and LC₅₀ values were <10 μg ml^?1. The LC₅₀/EC₅₀ ratios of almost all of the synthesized compounds were >10². The results indicate that these simple isocyanides are promising low-toxicity antifouling agents.     

Larval settlement and metamorphosis of the mussel Mytilus coruscus in response to natural biofilms

Settlement and metamorphosis of pediveliger larvae of Mytilus coruscus in response to natural biofilms was investigated in the laboratory. Pediveliger larvae settled and metamorphosed in response to biofilms and post-larval settlement and metamorphosis increased with biofilm age. The activity of the biofilm was positively correlated with biofilm age, dry weight, bacterial density and diatom density, but had no apparent relationship with chlorophyll a concentration. The change in bacterial community composition corresponding to biofilm age may explain differences in the age-dependent inducing activities of biofilms, which in turn may play an important role in larval settlement in this species.     

Cyanobacteria-containing biofilms from a Mayan monument in Palenque,Mexico

Surfaces of buildings at the archaeological site of Palenque, Mexico, are colonized by cyanobacteria that form biofilms, which in turn cause aesthetic and structural damage. The structural characterization and species composition of biofilms from the walls of one of these buildings, El Palacio, are reported. The distribution of photosynthetic microorganisms in the biofilms, their relationship with the colonized substratum, and the three-dimensional structure of the biofilms were studied by image analysis. The differences between local seasonal microenvironments at the Palenque site, the bioreceptivity of stone and the relationship between biofilms and their substrata are described. The implications for the development and permanence of species capable of withstanding temporal heterogeneity in and on El Palacio, mainly due to alternating wet and dry seasons, are discussed. Knowledge on how different biofilms contribute to biodegradation or bioprotection of the substratum can be used to develop maintenance and conservation protocols for cultural heritage.     

Molecular Dynamics Study on Protein and it's Water Structure at High Pressure

Abstract

We have performed NPT molecular dynamics simulations (Langevin Piston Method) on two types of solvated proteins-‘denaturation-unfavorable’ protein (insulin) and ‘denaturation-favorable protein’ (ribonuclease A) at high pressure (from 1 bar up to 20 kbar). The method is based on the extended system formalism introduced by Andersen, where the deterministic equations of motion for the piston degree of freedom are replaced by Langevin equation. We report the structural changes of proteins (ribonuclease A and insulin) and water molecules through radius of gyration, solvent accessible surface area, hydrogen bond pattern, and the topology of water clusters connected by the hydrogen bonded circular network. The solvent accessibility of ribonuclease A is mainly decreased by hydrophilic residues rather than hydrophobic residues under high pressure. From the results of hydrogen bond analysis, we have found that α-helix is more stable than β-sheet under high pressure. In addition, from the analysis of the water cluster, we have observed that for ribonuclease A, 5-membered ring structure is more favorable than 6-membered ring at higher pressure. However, for insulin, the ratio of 5 to 6-ring is constant over the pressure ranges for which we have performed MD simulation. This indicates that the water structure around insulin does not change under high pressure.     

Study of PE and iPP orientations on the surface of carbon nanotubes by using molecular dynamic simulations

In this study, the early stage of interfacial crystallisation behaviour of low molecular weight polyethylene (PE) and isotactic polypropylene (iPP) oligomer on the surface of carbon nanotubes (CNTs) with different diameters, chiralities and topography structures was studied using MD simulations. We started to simulate the effect of CNTs chirality and diameter on PE molecular chain orientation, and then the effect of CNTs topography structure on PE and iPP molecular chain orientation was investigated. Finally, some experiments were carried out to prove the simulated results. Our study shows that for CNTs with a diameter comparable with the radius of gyration (R_g) of a polymer chain, an easy orientation of PE chains along CNTs axis is observed for all the systems of the CNTs with different chiralities due to a geometric confinement effect. For CNTs with a much larger diameter, multiple orientation of PE chains is induced on its surface due to the lattice matching between graphite lattice and PE molecular chains. In this case, the chirality of CNTs dominates the orientation of graphite lattice, which determines the orientation of PE chains arrangement on CNTs surface. More importantly, it was found that the groove structure formed by CNT bundles is very useful for the stabilisation of polymer chain, and thus facilitates the orientation of molecular chain along the long axis of CNTs. As a result, a novel nanohybrid shish–kebab (NHSK) structure with CNTs acting as central shish while polymer lamellae as kebab can be successfully obtained for both PE with zigzag conformation and iPP with helical conformation. This simulation result was well supported by the experimental observation. Our study could provide not only a deep understanding of the origin of the polymer chain orientation on CNTs surface but also the guidance for the preparation of polymer/CNTs nanocomposites with novel NHSK structure.