Plasmon tunability in metallodielectric metamaterials

The dielectric properties of metamaterials consisting of periodically arranged metallic nanoparticles of spherical shape are calculated by rigorously solving Maxwell's equations. Effective dielectric functions are obtained by comparing the reflectivity of planar surfaces limiting these materials with Fresnel's formulas for equivalent homogeneous media, showing mixing and splitting of individual-particle modes due to inter-particle interaction. Detailed results for simple cubic and fcc crystals of aluminum spheres in vacuum, silver spheres in vacuum, and silver spheres in a silicon matrix are presented. The filling fraction of the metal f is shown to determine the position of the plasmon modes of these metamaterials. Significant deviations are observed with respect to Maxwell-Garnett effective medium theory for large f, and multiple plasmons are predicted to exist in contrast to Maxwell-Garnett theory.Comment: 6 pages, 4 figure

On the origin of the Boson peak in globular proteins

We study the Boson Peak phenomenology experimentally observed in globular proteins by means of elastic network models. These models are suitable for an analytic treatment in the framework of Euclidean Random Matrix theory, whose predictions can be numerically tested on real proteins structures. We find that the emergence of the Boson Peak is strictly related to an intrinsic mechanical instability of the protein, in close similarity to what is thought to happen in glasses. The biological implications of this conclusion are also discussed by focusing on a representative case study.Comment: Proceedings of the X International Workshop on Disordered Systems, Molveno (2006

Nonlinearity of Mechanochemical Motions in Motor Proteins

The assumption of linear response of protein molecules to thermal noise or structural perturbations, such as ligand binding or detachment, is broadly used in the studies of protein dynamics. Conformational motions in proteins are traditionally analyzed in terms of normal modes and experimental data on thermal fluctuations in such macromolecules is also usually interpreted in terms of the excitation of normal modes. We have chosen two important protein motors - myosin V and kinesin KIF1A - and performed numerical investigations of their conformational relaxation properties within the coarse-grained elastic network approximation. We have found that the linearity assumption is deficient for ligand-induced conformational motions and can even be violated for characteristic thermal fluctuations. The deficiency is particularly pronounced in KIF1A where the normal mode description fails completely in describing functional mechanochemical motions. These results indicate that important assumptions of the theory of protein dynamics may need to be reconsidered. Neither a single normal mode, nor a superposition of such modes yield an approximation of strongly nonlinear dynamics.Comment: 10 pages, 6 figure

Using Entropy Maximization to Understand the Determinants of Structural Dynamics beyond Native Contact Topology

Comparison of elastic network model predictions with experimental data has provided important insights on the dominant role of the network of inter-residue contacts in defining the global dynamics of proteins. Most of these studies have focused on interpreting the mean-square fluctuations of residues, or deriving the most collective, or softest, modes of motions that are known to be insensitive to structural and energetic details. However, with increasing structural data, we are in a position to perform a more critical assessment of the structure-dynamics relations in proteins, and gain a deeper understanding of the major determinants of not only the mean-square fluctuations and lowest frequency modes, but the covariance or the cross-correlations between residue fluctuations and the shapes of higher modes. A systematic study of a large set of NMR-determined proteins is analyzed using a novel method based on entropy maximization to demonstrate that the next level of refinement in the elastic network model description of proteins ought to take into consideration properties such as contact order (or sequential separation between contacting residues) and the secondary structure types of the interacting residues, whereas the types of amino acids do not play a critical role. Most importantly, an optimal description of observed cross-correlations requires the inclusion of destabilizing, as opposed to exclusively stabilizing, interactions, stipulating the functional significance of local frustration in imparting native-like dynamics. This study provides us with a deeper understanding of the structural basis of experimentally observed behavior, and opens the way to the development of more accurate models for exploring protein dynamics

On the conservation of the slow conformational dynamics within the amino acid kinase family: NAGK the paradigm

N-Acetyl-L-Glutamate Kinase (NAGK) is the structural paradigm for examining the catalytic mechanisms and dynamics of amino acid kinase family members. Given that the slow conformational dynamics of the NAGK (at the microseconds time scale or slower) may be rate-limiting, it is of importance to assess the mechanisms of the most cooperative modes of motion intrinsically accessible to this enzyme. Here, we present the results from normal mode analysis using an elastic network model representation, which shows that the conformational mechanisms for substrate binding by NAGK strongly correlate with the intrinsic dynamics of the enzyme in the unbound form. We further analyzed the potential mechanisms of allosteric signalling within NAGK using a Markov model for network communication. Comparative analysis of the dynamics of family members strongly suggests that the low-frequency modes of motion and the associated intramolecular couplings that establish signal transduction are highly conserved among family members, in support of the paradigm sequence→structure→dynamics→function © 2010 Marcos et al

Molecular Mechanics of the α-Actinin Rod Domain: Bending, Torsional, and Extensional Behavior

α-Actinin is an actin crosslinking molecule that can serve as a scaffold and maintain dynamic actin filament networks. As a crosslinker in the stressed cytoskeleton, α-actinin can retain conformation, function, and strength. α-Actinin has an actin binding domain and a calmodulin homology domain separated by a long rod domain. Using molecular dynamics and normal mode analysis, we suggest that the α-actinin rod domain has flexible terminal regions which can twist and extend under mechanical stress, yet has a highly rigid interior region stabilized by aromatic packing within each spectrin repeat, by electrostatic interactions between the spectrin repeats, and by strong salt bridges between its two anti-parallel monomers. By exploring the natural vibrations of the α-actinin rod domain and by conducting bending molecular dynamics simulations we also predict that bending of the rod domain is possible with minimal force. We introduce computational methods for analyzing the torsional strain of molecules using rotating constraints. Molecular dynamics extension of the α-actinin rod is also performed, demonstrating transduction of the unfolding forces across salt bridges to the associated monomer of the α-actinin rod domain

Variation in cartilage T2 and T2* mapping of the wrist: a comparison between 3- and 7-T MRI

Abstract Background To analyze regional variations in T2 and T2* relaxation times in wrist joint cartilage and the triangular fibrocartilage complex (TFCC) at 3 and 7 T and to compare values between field strengths. Methods Twenty-five healthy controls and 25 patients with chronic wrist pain were examined at 3 and 7 T on the same day using T2- and T2*-weighted sequences. Six different regions of interest (ROIs) were evaluated for cartilage and 3 ROIs were evaluated at the TFCC based on manual segmentation. Paired t-tests were used to compare T2 and T2* values between field strengths and between different ROIs. Spearman’s rank correlation was calculated to assess correlations between T2 and T2* time values at 3 and 7 T. Results T2 and T2* time values of the cartilage differed significantly between 3 and 7 T for all ROIs (p ≤ 0.045), with one exception: at the distal lunate, no significant differences in T2 values were observed between field strengths. T2* values differed significantly between 3 and 7 T for all ROIs of the TFCC (p ≤ 0.001). Spearman’s rank correlation between 3 and 7 T ranged from 0.03 to 0.62 for T2 values and from 0.01 to 0.48 for T2* values. T2 and T2* values for cartilage varied across anatomic locations in healthy controls at both 3 and 7 T. Conclusion Quantitative results of T2 and T2* mapping at the wrist differ between field strengths, with poor correlation between 3 and 7 T. Local variations in cartilage T2 and T2* values are observed in healthy individuals. Relevance statement T2 and T2* mapping are feasible for compositional imaging of the TFCC and the cartilage at the wrist at both 3 and 7 T, but the clinical interpretation remains challenging due to differences between field strengths and variations between anatomic locations. Key points •Field strength and anatomic locations influence T2 and T2* values at the wrist. •T2 and T2* values have a poor correlation between 3 and 7 T. •Local reference values are needed for each anatomic location for reliable interpretation. Graphical Abstract </jats:sec

A lightweight, flow-based toolkit for parallel and distributed bioinformatics pipelines

<p>Abstract</p> <p>Background</p> <p>Bioinformatic analyses typically proceed as chains of data-processing tasks. A pipeline, or 'workflow', is a well-defined protocol, with a specific structure defined by the topology of data-flow interdependencies, and a particular functionality arising from the data transformations applied at each step. In computer science, the dataflow programming (DFP) paradigm defines software systems constructed in this manner, as networks of message-passing components. Thus, bioinformatic workflows can be naturally mapped onto DFP concepts.</p> <p>Results</p> <p>To enable the flexible creation and execution of bioinformatics dataflows, we have written a modular framework for parallel pipelines in Python ('PaPy'). A PaPy workflow is created from re-usable components connected by data-pipes into a directed acyclic graph, which together define nested higher-order map functions. The successive functional transformations of input data are evaluated on flexibly pooled compute resources, either local or remote. Input items are processed in batches of adjustable size, all flowing one to tune the trade-off between parallelism and lazy-evaluation (memory consumption). An add-on module ('NuBio') facilitates the creation of bioinformatics workflows by providing domain specific data-containers (<it>e.g</it>., for biomolecular sequences, alignments, structures) and functionality (<it>e.g</it>., to parse/write standard file formats).</p> <p>Conclusions</p> <p>PaPy offers a modular framework for the creation and deployment of parallel and distributed data-processing workflows. Pipelines derive their functionality from user-written, data-coupled components, so PaPy also can be viewed as a lightweight toolkit for extensible, flow-based bioinformatics data-processing. The simplicity and flexibility of distributed PaPy pipelines may help users bridge the gap between traditional desktop/workstation and grid computing. PaPy is freely distributed as open-source Python code at <url>http://muralab.org/PaPy</url>, and includes extensive documentation and annotated usage examples.</p

On the Communication of Scientific Results: The Full-Metadata Format

In this paper, we introduce a scientific format for text-based data files, which facilitates storing and communicating tabular data sets. The so-called Full-Metadata Format builds on the widely used INI-standard and is based on four principles: readable self-documentation, flexible structure, fail-safe compatibility, and searchability. As a consequence, all metadata required to interpret the tabular data are stored in the same file, allowing for the automated generation of publication-ready tables and graphs and the semantic searchability of data file collections. The Full-Metadata Format is introduced on the basis of three comprehensive examples. The complete format and syntax is given in the appendix

Change in Allosteric Network Affects Binding Affinities of PDZ Domains: Analysis through Perturbation Response Scanning

The allosteric mechanism plays a key role in cellular functions of several PDZ domain proteins (PDZs) and is directly linked to pharmaceutical applications; however, it is a challenge to elaborate the nature and extent of these allosteric interactions. One solution to this problem is to explore the dynamics of PDZs, which may provide insights about how intramolecular communication occurs within a single domain. Here, we develop an advancement of perturbation response scanning (PRS) that couples elastic network models with linear response theory (LRT) to predict key residues in allosteric transitions of the two most studied PDZs (PSD-95 PDZ3 domain and hPTP1E PDZ2 domain). With PRS, we first identify the residues that give the highest mean square fluctuation response upon perturbing the binding sites. Strikingly, we observe that the residues with the highest mean square fluctuation response agree with experimentally determined residues involved in allosteric transitions. Second, we construct the allosteric pathways by linking the residues giving the same directional response upon perturbation of the binding sites. The predicted intramolecular communication pathways reveal that PSD-95 and hPTP1E have different pathways through the dynamic coupling of different residue pairs. Moreover, our analysis provides a molecular understanding of experimentally observed hidden allostery of PSD-95. We show that removing the distal third alpha helix from the binding site alters the allosteric pathway and decreases the binding affinity. Overall, these results indicate that (i) dynamics plays a key role in allosteric regulations of PDZs, (ii) the local changes in the residue interactions can lead to significant changes in the dynamics of allosteric regulations, and (iii) this might be the mechanism that each PDZ uses to tailor their binding specificities regulation