Dependence of Protein Adsorption on Wetting Behavior of UHMWPE-HA-Al2O3-CNT Hybrid Biocomposites

Ultrahigh-molecular-weight polyethylene (UHMWPE) is used as an articulating surface in total hip and knee joint replacement. In order to enhance long-term durability/wear resistance properties, UHMWPE-based polymer-ceramic hybrid composites are being developed. Surface properties such as wettability and protein adsorption alter with reinforcement or with change in surface chemistry. From this perspective, the wettability and protein adsorption behavior of compression-molded UHMWPE-hydroxyapatite (HA)-aluminum oxide (Al2O3)-carbon nanotube (CNT) composites were analyzed in conjunction with surface roughness. The combined effect of Al2O3 and CNT shows enhancement of the contact angle by similar to 37A degrees compared with the surface of the UHMWPE matrix reinforced with HA. In reference to unreinforced UHMWPE, protein adsorption density also increased by similar to 230% for 2 wt.%HA-5 wt.%Al2O3-2 wt.%CNT addition to UHMWPE. An important conclusion is that the polar and dispersion components of the surface free energy play a significant role in wetting and protein adsorption than do the total free energy or chemistry of the surface. The results of this study have major implications for the biocompatibility of these newly developed biocomposites.

Insights into the Fold Organization of TIM Barrel from Interaction Energy Based Structure Networks

There are many well-known examples of proteins with low sequence similarity, adopting the same structural fold. This aspect of sequence-structure relationship has been extensively studied both experimentally and theoretically, however with limited success. Most of the studies consider remote homology or ``sequence conservation'' as the basis for their understanding. Recently ``interaction energy'' based network formalism (Protein Energy Networks (PENs)) was developed to understand the determinants of protein structures. In this paper we have used these PENs to investigate the common non-covalent interactions and their collective features which stabilize the TIM barrel fold. We have also developed a method of aligning PENs in order to understand the spatial conservation of interactions in the fold. We have identified key common interactions responsible for the conservation of the TIM fold, despite high sequence dissimilarity. For instance, the central beta barrel of the TIM fold is stabilized by long-range high energy electrostatic interactions and low-energy contiguous vdW interactions in certain families. The other interfaces like the helix-sheet or the helix-helix seem to be devoid of any high energy conserved interactions. Conserved interactions in the loop regions around the catalytic site of the TIM fold have also been identified, pointing out their significance in both structural and functional evolution. Based on these investigations, we have developed a novel network based phylogenetic analysis for remote homologues, which can perform better than sequence based phylogeny. Such an analysis is more meaningful from both structural and functional evolutionary perspective. We believe that the information obtained through the ``interaction conservation'' viewpoint and the subsequently developed method of structure network alignment, can shed new light in the fields of fold organization and de novo computational protein design.

Cyclic AMP-induced Conformational Changes in Mycobacterial Protein Acetyltransferases

The activities of a number of proteins are regulated by the binding of cAMP and cGMP to cyclic nucleotide binding (CNB) domains that are found associated with one or more effector domains with diverse functions. Although the conserved architecture of CNB domains has been extensively studied by x-ray crystallography, the key to unraveling the mechanisms of cAMP action has been protein dynamics analyses. Recently, we have identified a novel cAMP-binding protein from mycobacteria, where cAMP regulates the activity of an associated protein acetyltransferase domain. In the current study, we have monitored the conformational changes that occur upon cAMP binding to the CNB domain in these proteins, using a combination of bioluminescence resonance energy transfer and amide hydrogen/deuterium exchange mass spectrometry. Coupled with mutational analyses, our studies reveal the critical role of the linker region (positioned between the CNB domain and the acetyltransferase domain) in allosteric coupling of cAMP binding to activation of acetyltransferase catalysis. Importantly, major differences in conformational change upon cAMP binding were accompanied by stabilization of the CNB and linker domain alone. This is in contrast to other cAMP-binding proteins, where cyclic nucleotide binding has been shown to involve intricate and parallel allosteric relays. Finally, this powerful convergence of results from bioluminescence resonance energy transfer and hydrogen/deuterium exchange mass spectrometry reaffirms the power of solution biophysical tools in unraveling mechanistic bases of regulation of proteins in the absence of high resolution structural information.

An automated approach to network features of protein structure ensembles

Network theory applied to protein structures provides insights into numerous problems of biological relevance. The explosion in structural data available from PDB and simulations establishes a need to introduce a standalone-efficient program that assembles network concepts/parameters under one hood in an automated manner. Herein, we discuss the development/application of an exhaustive, user-friendly, standalone program package named PSN-Ensemble, which can handle structural ensembles generated through molecular dynamics (MD) simulation/NMR studies or from multiple X-ray structures. The novelty in network construction lies in the explicit consideration of side-chain interactions among amino acids. The program evaluates network parameters dealing with topological organization and long-range allosteric communication. The introduction of a flexible weighing scheme in terms of residue pairwise cross-correlation/interaction energy in PSN-Ensemble brings in dynamical/chemical knowledge into the network representation. Also, the results are mapped on a graphical display of the structure, allowing an easy access of network analysis to a general biological community. The potential of PSN-Ensemble toward examining structural ensemble is exemplified using MD trajectories of an ubiquitin-conjugating enzyme (UbcH5b). Furthermore, insights derived from network parameters evaluated using PSN-Ensemble for single-static structures of active/inactive states of 2-adrenergic receptor and the ternary tRNA complexes of tyrosyl tRNA synthetases (from organisms across kingdoms) are discussed. PSN-Ensemble is freely available from http://vishgraph.mbu.iisc.ernet.in/PSN-Ensemble/psn_index.html.

Solvent Sensitivity of Protein Unfolding: Dynamical Study of Chicken Villin Headpiece Subdomain in Water Ethanol Binary Mixture

We carry out a series of long atomistic molecular dynamics simulations to study the unfolding of a small protein, chicken villin headpiece (HP-36), in water-ethanol (EtOH) binary mixture. The prime objective of this work is to explore the sensitivity of protein unfolding dynamics toward increasing concentration of the cosolvent and unravel essential features of intermediates formed in search of a dynamical pathway toward unfolding. In water ethanol binary mixtures, HP-36 is found to unfold partially, under ambient conditions, that otherwise requires temperature as high as similar to 600 K to denature in pure aqueous solvent. However, an interesting course of pathway is observed to be followed in the process, guided by the formation of unique intermediates. The first step of unfolding is essentially the separation of the cluster formed by three hydrophobic (phenylalanine) residues, namely, Phe-7, Phe-11, and Phe-18, which constitute the hydrophobic core, thereby initiating melting of helix-2 of the protein. The initial steps are similar to temperature-induced unfolding as well as chemical unfolding using DMSO as cosolvent. Subsequent unfolding steps follow a unique path. As water-ethanol shows composition-dependent anomalies, so do the details of unfolding dynamics. With an increase in cosolvent concentration, different partially unfolded intermediates are found to be formed. This is reflected in a remarkable nonmonotonic composition dependence of several order parameters, including fraction of native contacts and protein-solvent interaction energy. The emergence of such partially unfolded states can be attributed to the preferential solvation of the hydrophobic residues by the ethyl groups of ethanol. We further quantify the local dynamics of unfolding by using a Marcus-type theory.

The Presence of Multiple Cellular Defects Associated with a Novel G50E Iron-Sulfur Cluster Scaffold Protein ( ISCU) Mutation Leads to Development of Mitochondrial Myopathy

Background: Muscle-specific deficiency of iron-sulfur (Fe-S) cluster scaffold protein (ISCU) leads to myopathy. Results: Cells carrying the myopathy-associated G50E ISCU mutation demonstrate impaired Fe-S cluster biogenesis and mitochondrial dysfunction. Conclusion: Reduced mitochondrial respiration as a result of diminished Fe-S cluster synthesis results in muscle weakness in myopathy patients. Significance: The molecular mechanism behind disease progression should provide invaluable information to combat ISCU myopathy. Iron-sulfur (Fe-S) clusters are versatile cofactors involved in regulating multiple physiological activities, including energy generation through cellular respiration. Initially, the Fe-S clusters are assembled on a conserved scaffold protein, iron-sulfur cluster scaffold protein (ISCU), in coordination with iron and sulfur donor proteins in human mitochondria. Loss of ISCU function leads to myopathy, characterized by muscle wasting and cardiac hypertrophy. In addition to the homozygous ISCU mutation (g.7044GC), compound heterozygous patients with severe myopathy have been identified to carry the c.149GA missense mutation converting the glycine 50 residue to glutamate. However, the physiological defects and molecular mechanism associated with G50E mutation have not been elucidated. In this report, we uncover mechanistic insights concerning how the G50E ISCU mutation in humans leads to the development of severe ISCU myopathy, using a human cell line and yeast as the model systems. The biochemical results highlight that the G50E mutation results in compromised interaction with the sulfur donor NFS1 and the J-protein HSCB, thus impairing the rate of Fe-S cluster synthesis. As a result, electron transport chain complexes show significant reduction in their redox properties, leading to loss of cellular respiration. Furthermore, the G50E mutant mitochondria display enhancement in iron level and reactive oxygen species, thereby causing oxidative stress leading to impairment in the mitochondrial functions. Thus, our findings provide compelling evidence that the respiration defect due to impaired biogenesis of Fe-S clusters in myopathy patients leads to manifestation of complex clinical symptoms.

Propensity to Form Amyloid Fibrils Is Encoded as Excitations in the Free Energy Landscape of Monomeric Proteins

Protein aggregation, linked to many of diseases, is initiated when monomers access rogue conformations that are poised to form amyloid fibrils. We show, using simulations of src SH3 domain, that mechanical force enhances the population of the aggregation-prone (N*) states, which are rarely populated under force free native conditions but are encoded in the spectrum of native fluctuations. The folding phase diagrams of SH3 as a function of denaturant concentration (C]), mechanical force (f), and temperature exhibit an apparent two-state behavior, without revealing the presence of the elusive N* states. Interestingly, the phase boundaries separating the folded and unfolded states at all C] and f fall on a master curve, which can be quantitatively described using an analogy to superconductors in a magnetic field. The free energy profiles as a function of the molecular extension (R), which are accessible in pulling experiments, (R), reveal the presence of a native-like N* with a disordered solvent-exposed amino-terminal beta-strand. The structure of the N* state is identical with that found in Fyn SH3 by NMR dispersion experiments. We show that the timescale for fibril formation can be estimated from the population of the N* state, determined by the free energy gap separating the native structure and the N* state, a finding that can be used to assess fibril forming tendencies of proteins. The structures of the N* state are used to show that oligomer formation and likely route to fibrils occur by a domain-swap mechanism in SH3 domain. (C) 2014 Elsevier Ltd. All rights reserved.

Diffusion on a rugged energy landscape with spatial correlations

Rugged energy landscapes find wide applications in diverse fields ranging from astrophysics to protein folding. We study the dependence of diffusion coefficient (D) of a Brownian particle on the distribution width (epsilon) of randomness in a Gaussian random landscape by simulations and theoretical analysis. We first show that the elegant expression of Zwanzig Proc. Natl. Acad. Sci. U.S.A. 85, 2029 (1988)] for D(epsilon) can be reproduced exactly by using the Rosenfeld diffusion-entropy scaling relation. Our simulations show that Zwanzig's expression overestimates D in an uncorrelated Gaussian random lattice - differing by almost an order of magnitude at moderately high ruggedness. The disparity originates from the presence of ``three-site traps'' (TST) on the landscape - which are formed by the presence of deep minima flanked by high barriers on either side. Using mean first passage time formalism, we derive a general expression for the effective diffusion coefficient in the presence of TST, that quantitatively reproduces the simulation results and which reduces to Zwanzig's form only in the limit of infinite spatial correlation. We construct a continuous Gaussian field with inherent correlation to establish the effect of spatial correlation on random walk. The presence of TSTs at large ruggedness (epsilon >> k(B)T) gives rise to an apparent breakdown of ergodicity of the type often encountered in glassy liquids. (C) 2014 AIP Publishing LLC.

Multidimensional free energy surface of unfolding of HP-36: Microscopic origin of ruggedness

The protein folding funnel paradigm suggests that folding and unfolding proceed as directed diffusion in a multidimensional free energy surface where a multitude of pathways can be traversed during the protein's sojourn from initial to final state. However, finding even a single pathway, with the detail chronicling of intermediates, is an arduous task. In this work we explore the free energy surface of unfolding pathway through umbrella sampling, for a small globular a-helical protein chicken-villin headpiece (HP-36) when the melting of secondary structures is induced by adding DMSO in aqueous solution. We find that the unfolding proceeds through the initial separation or melting of aggregated hydrophobic core that comprises of three phenylalanine residues (Phe7, Phe11, and Phe18). This separation is accompanied by simultaneous melting of the second helix. Unfolding is found to be a multistage process involving crossing of three consecutive minima and two barriers at the initial stage. At a molecular level, Phe18 is observed to reorient itself towards other hydrophobic grooves to stabilize the intermediate states. We identify the configuration of the intermediates and correlate the intermediates with those obtained in our previous works. We also give an estimate of the barriers for different transition states and observe the softening of the barriers with increasing DMSO concentration. We show that higher concentration of DMSO tunes the unfolding pathway by destabilizing the third minimum and stabilizing the second one, indicating the development of a solvent modified, less rugged pathway. The prime outcome of this work is the demonstration that mixed solvents can profoundly transform the nature of the energy landscape and induce unfolding via a modified route. A successful application of Kramer's rate equation correlating the free energy simulation results shows faster rate of unfolding with increasing DMSO concentration. This work perhaps presents the first systematic theoretical study of the effect of a chemical denaturant on the microscopic free energy surface and rates of unfolding of HP-36. (C) 2014 AIP Publishing LLC.

Folding of Protein L with Implications for Collapse in the Denatured State Ensemble

A fundamental question in protein folding is whether the coil to globule collapse transition occurs during the initial stages of folding (burst phase) or simultaneously with the protein folding transition. Single molecule fluorescence resonance energy transfer (FRET) and small-angle X-ray scattering (SAXS) experiments disagree on whether Protein L collapse transition occurs during the burst phase of folding. We study Protein L folding using a coarse-grained model and molecular dynamics simulations. The collapse transition in Protein L is found to be concomitant with the folding transition. In the burst phase of folding, we find that FRET experiments overestimate radius of gyration, R-g, of the protein due to the application of Gaussian polymer chain end-to-end distribution to extract R-g from the FRET efficiency. FRET experiments estimate approximate to 6 angstrom decrease in R-g when the actual decrease is approximate to 3 angstrom on guanidinium chloride denaturant dilution from 7.5 to 1 M, thereby suggesting pronounced compaction in the protein dimensions in the burst phase. The approximate to 3 angstrom decrease is close to the statistical uncertainties of the R-g data measured from SAXS experiments, which suggest no compaction, leading to a disagreement with the FRET experiments. The transition-state ensemble (TSE) structures in Protein L folding are globular and extensive in agreement with the Psi-analysis experiments. The results support the hypothesis that the TSE of single domain proteins depends on protein topology and is not stabilized by local interactions alone.