Changes in Lysozyme Flexibility upon Mutation Are Frequent, Large and Long-Ranged

We investigate changes in human c-type lysozyme flexibility upon mutation via a Distance Constraint Model, which gives a statistical mechanical treatment of network rigidity. Specifically, two dynamical metrics are tracked. Changes in flexibility index quantify differences within backbone flexibility, whereas changes in the cooperativity correlation quantify differences within pairwise mechanical couplings. Regardless of metric, the same general conclusions are drawn. That is, small structural perturbations introduced by single point mutations have a frequent and pronounced affect on lysozyme flexibility that can extend over long distances. Specifically, an appreciable change occurs in backbone flexibility for 48% of the residues, and a change in cooperativity occurs in 42% of residue pairs. The average distance from mutation to a site with a change in flexibility is 17–20 Å. Interestingly, the frequency and scale of the changes within single point mutant structures are generally larger than those observed in the hen egg white lysozyme (HEWL) ortholog, which shares 61% sequence identity with human lysozyme. For example, point mutations often lead to substantial flexibility increases within the β-subdomain, which is consistent with experimental results indicating that it is the nucleation site for amyloid formation. However, β-subdomain flexibility within the human and HEWL orthologs is more similar despite the lowered sequence identity. These results suggest compensating mutations in HEWL reestablish desired properties

On Conduction in a Bacterial Sodium Channel

Voltage-gated Na+-channels are transmembrane proteins that are responsible for the fast depolarizing phase of the action potential in nerve and muscular cells. Selective permeability of Na+ over Ca2+ or K+ ions is essential for the biological function of Na+-channels. After the emergence of the first high-resolution structure of a Na+-channel, an anionic coordination site was proposed to confer Na+ selectivity through partial dehydration of Na+ via its direct interaction with conserved glutamate side chains. By combining molecular dynamics simulations and free-energy calculations, a low-energy permeation pathway for Na+ ion translocation through the selectivity filter of the recently determined crystal structure of a prokaryotic sodium channel from Arcobacter butzleri is characterised. The picture that emerges is that of a pore preferentially occupied by two ions, which can switch between different configurations by crossing low free-energy barriers. In contrast to K+-channels, the movements of the ions appear to be weakly coupled in Na+-channels. When the free-energy maps for Na+ and K+ ions are compared, a selective site is characterised in the narrowest region of the filter, where a hydrated Na+ ion, and not a hydrated K+ ion, is energetically stable

Replication Pauses of the Wild-Type and Mutant Mitochondrial DNA Polymerase Gamma: A Simulation Study

The activity of polymerase γ is complicated, involving both correct and incorrect DNA polymerization events, exonuclease activity, and the disassociation of the polymerase:DNA complex. Pausing of pol-γ might increase the chance of deletion and depletion of mitochondrial DNA. We have developed a stochastic simulation of pol-γ that models its activities on the level of individual nucleotides for the replication of mtDNA. This method gives us insights into the pausing of two pol-γ variants: the A467T substitution that causes PEO and Alpers syndrome, and the exonuclease deficient pol-γ (exo−) in premature aging mouse models. To measure the pausing, we analyzed simulation results for the longest time for the polymerase to move forward one nucleotide along the DNA strand. Our model of the exo− polymerase had extremely long pauses, with a 30 to 300-fold increase in the time required for the longest single forward step compared to the wild-type, while the naturally occurring A467T variant showed at most a doubling in the length of the pauses compared to the wild-type. We identified the cause of these differences in the polymerase pausing time to be the number of disassociations occurring in each forward step of the polymerase

Binding Free Energy Landscape of Domain-Peptide Interactions

Peptide recognition domains (PRDs) are ubiquitous protein domains which mediate large numbers of protein interactions in the cell. How these PRDs are able to recognize peptide sequences in a rapid and specific manner is incompletely understood. We explore the peptide binding process of PDZ domains, a large PRD family, from an equilibrium perspective using an all-atom Monte Carlo (MC) approach. Our focus is two different PDZ domains representing two major PDZ classes, I and II. For both domains, a binding free energy surface with a strong bias toward the native bound state is found. Moreover, both domains exhibit a binding process in which the peptides are mostly either bound at the PDZ binding pocket or else interact little with the domain surface. Consistent with this, various binding observables show a temperature dependence well described by a simple two-state model. We also find important differences in the details between the two domains. While both domains exhibit well-defined binding free energy barriers, the class I barrier is significantly weaker than the one for class II. To probe this issue further, we apply our method to a PDZ domain with dual specificity for class I and II peptides, and find an analogous difference in their binding free energy barriers. Lastly, we perform a large number of fixed-temperature MC kinetics trajectories under binding conditions. These trajectories reveal significantly slower binding dynamics for the class II domain relative to class I. Our combined results are consistent with a binding mechanism in which the peptide C terminal residue binds in an initial, rate-limiting step

The cerebrospinal fluid proteome in HIV infection: change associated with disease severity

<p>Abstract</p> <p>Background</p> <p>Central nervous system (CNS) infection is a nearly universal feature of untreated systemic HIV infection with a clinical spectrum that ranges from chronic asymptomatic infection to severe cognitive and motor dysfunction. Analysis of cerebrospinal fluid (CSF) has played an important part in defining the character of this evolving infection and response to treatment. To further characterize CNS HIV infection and its effects, we applied advanced high-throughput proteomic methods to CSF to identify novel proteins and their changes with disease progression and treatment.</p> <p>Results</p> <p>After establishing an <it>accurate mass and time </it>(AMT) tag database containing 23,141 AMT tags for CSF peptides, we analyzed 91 CSF samples by LC-MS from 12 HIV-uninfected and 14 HIV-infected subjects studied in the context of initiation of antiretroviral therapy and correlated abundances of identified proteins a) within and between subjects, b) with all other proteins across the entire sample set, and c) with "external" CSF biomarkers of infection (HIV RNA), immune activation (neopterin) and neural injury (neurofilament light chain protein, NFL). We identified a mean of 2,333 +/- 328 (SD) peptides covering 307 +/-16 proteins in the 91 CSF sample set. Protein abundances differed both between and within subjects sampled at different time points and readily separated those with and without HIV infection. Proteins also showed inter-correlations across the sample set that were associated with biologically relevant dynamic processes. One-hundred and fifty proteins showed correlations with the external biomarkers. For example, using a threshold of cross correlation coefficient (Pearson's) ≤ -0.3 and ≥0.3 for potentially meaningful relationships, a total of 99 proteins correlated with CSF neopterin (43 negative and 56 positive correlations) and related principally to neuronal plasticity and survival and to innate immunity. Pathway analysis defined several networks connecting the identified proteins, including one with amyloid precursor protein as a central node.</p> <p>Conclusions</p> <p>Advanced CSF proteomic analysis enabled the identification of an array of novel protein changes across the spectrum of CNS HIV infection and disease. This initial analysis clearly demonstrated the value of contemporary state-of-the-art proteomic CSF analysis as a discovery tool in HIV infection with likely similar application to other neurological inflammatory and degenerative diseases.</p

Corresponding Functional Dynamics across the Hsp90 Chaperone Family: Insights from a Multiscale Analysis of MD Simulations

Understanding how local protein modifications, such as binding small-molecule ligands, can trigger and regulate large-scale motions of large protein domains is a major open issue in molecular biology. We address various aspects of this problem by analyzing and comparing atomistic simulations of Hsp90 family representatives for which crystal structures of the full length protein are available: mammalian Grp94, yeast Hsp90 and E.coli HtpG. These chaperones are studied in complex with the natural ligands ATP, ADP and in the Apo state. Common key aspects of their functional dynamics are elucidated with a novel multi-scale comparison of their internal dynamics. Starting from the atomic resolution investigation of internal fluctuations and geometric strain patterns, a novel analysis of domain dynamics is developed. The results reveal that the ligand-dependent structural modulations mostly consist of relative rigid-like movements of a limited number of quasi-rigid domains, shared by the three proteins. Two common primary hinges for such movements are identified. The first hinge, whose functional role has been demonstrated by several experimental approaches, is located at the boundary between the N-terminal and Middle-domains. The second hinge is located at the end of a three-helix bundle in the Middle-domain and unfolds/unpacks going from the ATP- to the ADP-state. This latter site could represent a promising novel druggable allosteric site common to all chaperones

Comparative Omics-Driven Genome Annotation Refinement: Application across Yersiniae

Genome sequencing continues to be a rapidly evolving technology, yet most downstream aspects of genome annotation pipelines remain relatively stable or are even being abandoned. The annotation process is now performed almost exclusively in an automated fashion to balance the large number of sequences generated. One possible way of reducing errors inherent to automated computational annotations is to apply data from omics measurements (i.e. transcriptional and proteomic) to the un-annotated genome with a proteogenomic-based approach. Here, the concept of annotation refinement has been extended to include a comparative assessment of genomes across closely related species. Transcriptomic and proteomic data derived from highly similar pathogenic Yersiniae (Y. pestis CO92, Y. pestis Pestoides F, and Y. pseudotuberculosis PB1/+) was used to demonstrate a comprehensive comparative omic-based annotation methodology. Peptide and oligo measurements experimentally validated the expression of nearly 40% of each strain's predicted proteome and revealed the identification of 28 novel and 68 incorrect (i.e., observed frameshifts, extended start sites, and translated pseudogenes) protein-coding sequences within the three current genome annotations. Gene loss is presumed to play a major role in Y. pestis acquiring its niche as a virulent pathogen, thus the discovery of many translated pseudogenes, including the insertion-ablated argD, underscores a need for functional analyses to investigate hypotheses related to divergence. Refinements included the discovery of a seemingly essential ribosomal protein, several virulence-associated factors, a transcriptional regulator, and many hypothetical proteins that were missed during annotation

Ligand-Dependent Conformations and Dynamics of the Serotonin 5-HT2A Receptor Determine Its Activation and Membrane-Driven Oligomerization Properties

From computational simulations of a serotonin 2A receptor (5-HT2AR) model complexed with pharmacologically and structurally diverse ligands we identify different conformational states and dynamics adopted by the receptor bound to the full agonist 5-HT, the partial agonist LSD, and the inverse agonist Ketanserin. The results from the unbiased all-atom molecular dynamics (MD) simulations show that the three ligands affect differently the known GPCR activation elements including the toggle switch at W6.48, the changes in the ionic lock between E6.30 and R3.50 of the DRY motif in TM3, and the dynamics of the NPxxY motif in TM7. The computational results uncover a sequence of steps connecting these experimentally-identified elements of GPCR activation. The differences among the properties of the receptor molecule interacting with the ligands correlate with their distinct pharmacological properties. Combining these results with quantitative analysis of membrane deformation obtained with our new method (Mondal et al, Biophysical Journal 2011), we show that distinct conformational rearrangements produced by the three ligands also elicit different responses in the surrounding membrane. The differential reorganization of the receptor environment is reflected in (i)-the involvement of cholesterol in the activation of the 5-HT2AR, and (ii)-different extents and patterns of membrane deformations. These findings are discussed in the context of their likely functional consequences and a predicted mechanism of ligand-specific GPCR oligomerization

Role of biomechanics in the understanding of normal, injured, and healing ligaments and tendons

Ligaments and tendons are soft connective tissues which serve essential roles for biomechanical function of the musculoskeletal system by stabilizing and guiding the motion of diarthrodial joints. Nevertheless, these tissues are frequently injured due to repetition and overuse as well as quick cutting motions that involve acceleration and deceleration. These injuries often upset this balance between mobility and stability of the joint which causes damage to other soft tissues manifested as pain and other morbidity, such as osteoarthritis