Destruction of long-range interactions by a single mutation in lysozyme

[Modulation of Aggregation and Immunogenicity of a Protein: Based on the Study of Hen Lysozyme]

This study focuses on the modulation of protein aggregation and immunogenicity. As a starting point for investigating long-range interactions within a non-native protein, the effects of perturbing denatured protein states on their aggregation, including the formation of amyloid fibrils, were evaluated. The effects of adducts, sugar modifications, and stabilization on protein aggregation were then examined. We also investigated how protein immunogenicity was affected by enhancing protein conformational stability and other factors.

Roles of Tryptophan and Charged Residues on the Polymorphisms of Amyloids Formed by K-Peptides of Hen Egg White Lysozyme Investigated through Molecular Dynamics Simulations

Atomistic molecular dynamics simulations of amyloid models, consisting of the previously reported STDY-K-peptides and K-peptides from the hen egg white lysozyme (HEWL), were performed to address the effects of charged residues and pH observed in an in vitro study. Simulation results showed that amyloid models with antiparallel configurations possessed greater stability and compactness than those with parallel configurations. Then, peptide chain stretching and ordering were measured through the end-to-end distance and the order parameter, for which the amyloid models consisting of K-peptides and the STDY-K-peptides at pH 2 displayed a higher level of chain stretching and ordering. After that, the molecular mechanics energy decomposition and the radial distribution function (RDF) clearly displayed the importance of Trp62 to the K-peptide and the STDY-K-peptide models at pH 2. Moreover, the results also displayed how the negatively charged Asp52 disrupted the interaction networks and prevented the amyloid formation from STDY-K-peptide at pH 7. Finally, this study provided an insight into the interplay between pH conditions and molecular interactions underlying the formation of amyloid fibrils from short peptides contained within the HEWL. This served as a basis of understanding towards the design of other amyloids for biomaterial applications.

Analyses of Mutation Displacements from Homology Models

Evaluation of the structural perturbations introduced by a single amino acid mutation is the main issue for protein structural biology. We propose here to present some recent advances in methods, allowing the splitting of distortion between the actual substitution effect and the contribution of the local flexibility of the position where the mutation occurs. Its main drawback is the need of many structures with a single mutation in each of them. To bypass this difficulty, we propose to use molecular modeling tools, with several software enabling us to build a model from a template, given the sequence. As a proof of concept, we rely on a gold standard, the human lysozyme. Both wild-type and three mutant structures are available in the PDB. Two of these mutations result in amyloid fibril formation, and the last one is neutral. As a conclusion, irrespective of the algorithm used for modeling, side chain conformations at the site of mutation are reliable, although long-range effects are out of reach of these tools.

The effects of mutation and modification on the structure and stability of human lysozyme: A molecular link between carbamylation and atherosclerosis

Amino acid mutations in some proteins such as lysozyme lead to genetically disorder variants and adverse pathogenic consequences. Recently, amino acid modifications were known as a risk factor in many related diseases such as uremia and atherosclerosis, showing the importance of these surface-structure changes. Although the structural consequences of the hereditary proteins have been examined extensively, such effects for the protein modifications are known to a lesser extent. One drawback in the examination of protein modifications is hardness in experimental detection of modifications by techniques such as NMR and crystallography. Molecular modeling and simulation can help to understand such phenomena at the molecular levels. It is more rational that the effects of both mutation and modification can be compared in a single protein model. Here, molecular dynamics simulation is used to compare the effects of a disease-related carbamylation modification and an amyloidogenic mutation (D67H) in human lysozyme as a model protein. The results show that the carbamylation adversely effects on the tertiary structure, leading to the similar unfolding pathway to the hereditary amyloidogenic form. The carbamylation leads to the instability of the overall protein conformation, especially on the β-domain, which is a characteristic of hereditary amyloidosis in human lysozymes. The aggregation behaviors of both modified and mutant lysozyme were examined by molecular docking calculations. The results showed that the partially unfolded lysozyme might form tight protein aggregates upon carbamylation similar to the amyloidogenic variant. Both single and all-residues carbamylations impose serious conformational changes to the tertiary structure of lysozyme. It was obtained that carbamylation of lysozyme strongly effects on the stability of N-terminal β-sheet, which can produce a highly unstable conformation. The results of this study not only show the adverse structural consequences of a disease-associated post-translational modification, but it also may be very helpful to understand the molecular basis for many carbamylation-related diseases such as atherosclerosis in ESRD patients. The results show that non-native post-translational modifications may be as structurally important as hereditary mutations.

Analyses of displacements resulting from a point mutation in proteins

The effects of a single residue substitution on the protein backbone are frequently quite small and there are many other potential sources of structural variation for protein. We present here a methodology considering different sources of distortions in order to isolate the very effect of the mutation. To validate our methodology, we consider a well-studied family with many single mutants: the human lysozyme. Most of the perturbations are expected to be at the very localisation of the mutation, but in many cases the effects are propagated at long range. We show that the distances between the mutated residue and the 5% most disturbed residues exponentially decreases. One third of the affected residues are in direct contact with the mutated position; the remaining two thirds are potential allosteric effects. We confirm the reliability of the residues identified as significantly perturbed by comparing our results to experimental studies. We confirm with the present method all the previously identified perturbations. This study shows that mutations have long-range impact on protein backbone that can be detected, although the displacement of the affected atoms is small.

Polyamines and its analogue modulates amyloid fibrillation in lysozyme: A comparative investigation

BACKGROUND

Polyamines can induce protein aggregation that can be related to the physiology of the cellular function. Polyamines have been implicated in protein aggregation which may lead to neuropathic and non neuropathic amyloidosis.

SCOPE OF REVIEW

Change in the level of polyamine concentration has been associated with ageing and neurodegeneration such as Parkinson's disease, Alzheimer's disease. Lysozyme aggregation in the presence of polyamines leads to non neuropathic amyloidosis. Polyamine analogues can suppress or inhibit protein aggregation suggesting their efficacy against amyloidogenic protein aggregates.

MAJOR CONCLUSIONS

In this study we report the comparative interactions of lysozyme with the polyamine analogue, 1-naphthyl acetyl spermine in comparison with the biogenic polyamines through spectroscopy, calorimetry, imaging and docking techniques. The findings revealed that the affinity of binding varied as spermidine > 1-naphthyl acetyl spermine > spermine. The biogenic polyamines accelerated the rate of fibrillation significantly, whereas the analogue inhibited the rate of fibrillation to a considerable extent. The polyamines bind near the catalytic diad residues viz. Glu35 and Asp52, and in close proximity of Trp62 residue. However, the analogue showed dual nature of interaction where its alkyl amine region bind in same way as the biogenic polyamines bind to the catalytic site, while the naphthyl group makes hydrophobic contacts with Trp62 and Trp63, thereby suggesting its direct influence on fibrillation.

GENERAL SIGNIFICANCE

This study, thus, potentiates, the development of a polyamine analogue that can perform as an effective inhibitor targeted towards aggregation of amyloidogenic proteins.

Dual mechanism of ionic liquid-induced protein unfolding

Ionic liquids (ILs) are gaining attention as protein stabilizers and refolding additives.

Molecular Origin of the Stability Difference in Four Shark IgNAR Constant Domains

Improving the stability of antibodies for manufacture and shelf life is one of the main focuses of antibody engineering. One stabilization strategy is to perform specific mutations in human antibodies based on highly stable antibodies in other species. To identify the key residues for mutagenesis, it is necessary to understand the roles of these residues in stabilizing the antibody. Here, we use molecular dynamics simulations to study the molecular origin of the four shark immunoglobulin new antigen receptors constant domains (C1-C4). According to the unfolding pathways and the conformational free energy surfaces in 8 M urea at 380 K, the C2 domain is the most stable, followed by C4, C1, and C3, which agrees with the experimental findings. The C1 and C3 domains follow a common unfolding pathway in which the unfolding starts from the edge strands, particularly strand g, and then gradually progresses to the inner strands. Detailed structural analysis of the C2 domain reveals a "sandwich-like" R339-E322-R341 salt-bridge cluster on strand g, which grants ultrahigh stability to the C2 domain. We further design two sets of mutations by mutating E322 to alanine or setting all atomic charges in E322 to zero to break the salt-bridge cluster in the C2 domain, which confirms the importance of the salt bridges in stability. In the C4 domain, the D80-K104 salt bridge on strand g also strengthens the stability. On the other hand, in the C1 and C3 domains, there is no salt bridge on strand g. In addition to the salt bridges, the overall hydrophobicity score of the hydrophobic core is also positively correlated with the domain stability. Our findings provide a detailed microscopic picture of the molecular origin of the four shark immunoglobulin new antigen receptors constant domains that not only explains the differences in their structural stability but also provides important insights into future antibody design.

In silico screening of cancer-associated mutations in the HSA domain of BRG1 and its role in affecting the Arp-HSA sub-complex of SWI/SNF

SWI/SNF (SWItch/Sucrose Non-Fermentable) complexes regulate the gene expression programs by remodeling the nucleosome architecture of the chromatin functional elements. These large multi-component complexes comprise eight or more subunits and are conserved from yeast to human. Noticeably, nuclear actin and actin-related proteins (Arps) are an integral part of these complexes and are known to directly interact with the helicase-SANT-associated (HSA) domain of ATPase subunit. Recently, SWI/SNF subunits are gaining importance because of the prevalence of cancer-causing mutations associated with them. The functional characterization of the mutations in the SWI/SNF subunits is important for understanding their role in tumorigenesis and identifying potential therapeutic strategies. To study the actin-related complex of human SWI/SNF and the cancer-associated mutations interfering Arp assembly with the ATPase subunit, we modelled the structure of the β-actin-BAF53A-HSA complex based on the yeast Arp-HSA complex (PDB ID: 4I6M). Seven deleterious mutations in the HSA domain of BRG1 were identified based on the functional screening of cancer-associated mutations in the COSMIC database. Detailed structural analysis of the six mutations (R466H, R469W, Y489C, K502N, R513Q and R521P) based on molecular dynamics (MD) simulations reveal the distinct effect of each mutation in destabilizing the structure of the Arp-HSA complex. Predominantly we could notice the long-range effect of the HSA mutations in influencing the dynamics of the Arp subunits.

Pathogenic mutation R959W alters recognition dynamics of dysferlin inner DysF domain

Dysferlin, a 220 kD protein, plays a major role in regulating plasma membrane repair in muscle cells. Mutations in the dysferlin inner DysF domain are known to cause different types of muscular dystrophy, including limb-girdle muscular dystrophy type 2B (LGMD2B) and Miyoshi myopathy (MM). Replacement of arginine in position 959 by tryptophan has been frequently associated with both LGMD2B and MM, but the molecular mechanisms by which this mutation alters dysferlin function remain unknown. In this study, we have used protein binding site predictions and microsecond molecular dynamics (MD) simulations to determine the effects pathogenic mutation R959W on the structural dynamics of dysferlin inner DysF domain. Analysis of 2 μs long MD trajectories revealed that mutation R959W does not induce local destabilization, unfolding or misfolding of the domain. We used a binding site predictor to discover a protein-binding site (residues T958-I966 and E1031-H1037) that resembles pincers in shape. Cartesian principal component analysis and interresidue distance distributions of the wild-type domain showed that the predicted protein-binding site undergoes a pincer motion, and populates two structural states, open and closed. We found that mutation R959W inhibits the pincer motion of the protein-binding site and completely shifts the equilibrium toward the open state. These differences in the structural dynamics of the predicted binding site suggest that mutation R959W alters recognition dynamics of the inner DysF domain. Based on these findings and on previous experimental studies, we propose a novel role for the inner DysF domain in muscle membrane repair through recruitment of dysferlin to plasma membrane. In conclusion, these findings have important implications for our understanding of the structural aspects of muscular dystrophies in atomic-level resolution.

Nanomechanics of Protein Unfolding Outside a Generic Nanopore

Protein folding and unfolding have been the subject of active research for decades. Most of previous studies in protein unfolding were focused on temperature, chemical, and/or force (such as in atomic force microscopy (AFM)) induced denaturations. Recent studies on the functional roles of proteasomes (such as ClpXP) revealed a different unfolding process in cell, during which a target protein is mechanically unfolded and pulled into a confined, pore-like geometry for degradation. While the proteasome nanomachine has been extensively studied, the mechanism for unfolding proteins with the proteasome pore is still poorly understood. Here, we investigate the mechanical unfolding process of ubiquitin with (or really outside) a generic nanopore, and compare such process with that in the AFM pulling experiment. Unexpectedly, the required force for protein unfolding through a pore can be much smaller than that by the AFM. Simulation results also unveiled different nanomechanics, tearing fracture vs shearing friction, in these two distinct types of mechanical unfoldings.

Empirical Optimization of Interactions between Proteins and Chemical Denaturants in Molecular Simulations

Chemical denaturants are the most commonly used perturbation applied to study protein stability and folding kinetics as well as the properties of unfolded polypeptides. We build on recent work balancing the interactions of proteins and water, and accurate models for the solution properties of urea and guanidinium chloride, to develop a combined force field that is able to capture the strength of interactions between proteins and denaturants. We use solubility data for a model tetraglycine peptide in each denaturant to tune the protein-denaturant interaction by a novel simulation methodology. We validate the results against data for more complex sequences: single-molecule Förster resonance energy transfer data for a 34-residue fragment of the globular protein CspTm and photoinduced electron transfer quenching data for the disordered peptides C(AGQ)nW in denaturant solution as well as the chemical denaturation of the mini-protein Trp cage. The combined force field model should aid our understanding of denaturation mechanisms and the interpretation of experiment.

Chelerythrine–lysozyme interaction: spectroscopic studies, thermodynamics and molecular modeling exploration

The binding of the iminium and alkanolamine forms of chelerythrine to lysozyme (Lyz) was investigated by spectroscopy and molecular modeling studies.

Bio-mimicking of proline-rich motif applied to carbon nanotube reveals unexpected subtleties underlying nanoparticle functionalization

Here, we report computational studies of the SH3 protein domain interacting with various single-walled carbon nanotubes (SWCNT) either bare or functionalized by mimicking the proline-rich motif (PRM) ligand (PPPVPPRR) and compare it to the SH3-PRM complex binding. With prolines or a single arginine attached, the SWCNT gained slightly on specificity when compared with the bare control, whereas with multi-arginine systems the specificity dropped dramatically to our surprise. Although the electrostatic interaction provided by arginines is crucial in the recognition between PRM and SH3 domain, our results suggest that attaching multiple arginines to the SWCNT has a detrimental effect on the binding affinity. Detailed analysis of the MD trajectories found two main factors that modulate the specificity of the binding: the existence of competing acidic patches at the surface of SH3 that leads to "trapping and clamping" by the arginines, and the rigidity of the SWCNT introducing entropic penalties in the proper binding. Further investigation revealed that the same "clamping" phenomenon exits in the PRM-SH3 system, which has not been reported in previous literature. The competing effects between nanoparticle and its functionalization components revealed by our model system should be of value to current and future nanomedicine designs.

Cytotoxicity of graphene: recent advances and future perspective

Graphene, a unique two-dimensional single-atom-thin nanomaterial with exceptional structural, mechanical, and electronic properties, has spurred an enormous interest in many fields, including biomedical applications, which at the same time ignites a growing concern on its biosafety and potential cytotoxicity to human and animal cells. In this review, we present a summary of some very recent studies on this important subject with both experimental and theoretical approaches. The molecular interactions of graphene with proteins, DNAs, and cell membranes (both bacteria and mammalian cells) are discussed in detail. Severe distortions in structures and functions of these biomacromolecules by graphene are identified and characterized. For example, the graphene is shown to disrupt bacteria cell membranes by insertion/cutting as well as destructive extraction of lipid molecules directly. More interestingly, this cytotoxicity has been shown to have implications in de novo design of nanomedicine, such as graphene-based band-aid, a potential 'green' antibiotics due to its strong physical-based (instead of chemical-based) antibacterial capability. These studies have provided a better understanding of graphene nanotoxicity at both cellular and molecular levels, and also suggested therapeutic potential by using graphene's cytotoxicity against bacteria cells.

The complex and specific pMHC interactions with diverse HIV-1 TCR clonotypes reveal a structural basis for alterations in CTL function

Immune control of viral infections is modulated by diverse T cell receptor (TCR) clonotypes engaging peptide-MHC class I complexes on infected cells, but the relationship between TCR structure and antiviral function is unclear. Here we apply in silico molecular modeling with in vivo mutagenesis studies to investigate TCR-pMHC interactions from multiple CTL clonotypes specific for a well-defined HIV-1 epitope. Our molecular dynamics simulations of viral peptide-HLA-TCR complexes, based on two independent co-crystal structure templates, reveal that effective and ineffective clonotypes bind to the terminal portions of the peptide-MHC through similar salt bridges, but their hydrophobic side-chain packings can be very different, which accounts for the major part of the differences among these clonotypes. Non-specific hydrogen bonding to viral peptide also accommodates greater epitope variants. Furthermore, free energy perturbation calculations for point mutations on the viral peptide KK10 show excellent agreement with in vivo mutagenesis assays, with new predictions confirmed by additional experiments. These findings indicate a direct structural basis for heterogeneous CTL antiviral function.

Dissecting the contributions of β-hairpin tyrosine pairs to the folding and stability of long-lived human γD-crystallins

Ultraviolet-radiation-induced damage to and aggregation of human lens crystallin proteins are thought to be a significant pathway to age-related cataract. The aromatic residues within the duplicated Greek key domains of γ- and β-crystallins are the main ultraviolet absorbers and are susceptible to direct and indirect ultraviolet damage. The previous site-directed mutagenesis studies have revealed a striking difference for two highly conserved homologous β-hairpin Tyr pairs, at the N-terminal domain (N-td) and C-terminal domain (C-td), respectively, in their contribution to the overall stability of HγD-Crys, but why they behave so differently still remains a mystery. In this paper, we systematically investigated the underlying molecular mechanism and detailed contributions of these two Tyr pairs with large scale molecular dynamics simulations. A series of different tyrosine-to-alanine pair(s) substitutions were performed in either the N-td, the C-td, or both. Our results suggest that the Y45A/Y50A pair substitution in the N-td mainly affects the stability of the N-td itself, while the Y133A/Y138A pair substitution in the C-td leads to a more cooperative unfolding of both N-td and C-td. The stability of motif 2 in the N-td is mainly determined by the interdomain interface, while motif 1 in the N-td or motifs 3 and 4 in the C-td are mainly stabilized by the intradomain hydrophobic core. The damage to any tyrosine pair(s) can directly introduce some apparent water leakage to the hydrophobic core at the interface, which in turn causes a serious loss in the stability of the N-td. However, for the C-td substitutions, it may further impair the stable "sandwich-like" Y133-R167-Y138 cluster (through cation-π interactions) in the wild-type, thus causing the loop regions near the residue A138 to undergo large fluctuations, which in turn results in the intrusion of water into the hydrophobic core of the C-td and induces the C-td to lose its stability. These findings help resolve the "mystery" on why these two Tyr pairs display such a striking difference in their contributions to the overall protein stability despite their highly homologous nature.

UV-radiation induced disruption of dry-cavities in human γD-crystallin results in decreased stability and faster unfolding

Age-onset cataracts are believed to be expedited by the accumulation of UV-damaged human γD-crystallins in the eye lens. Here we show with molecular dynamics simulations that the stability of γD-crystallin is greatly reduced by the conversion of tryptophan to kynurenine due to UV-radiation, consistent with previous experimental evidences. Furthermore, our atomic-detailed results reveal that kynurenine attracts more waters and other polar sidechains due to its additional amino and carbonyl groups on the damaged tryptophan sidechain, thus breaching the integrity of nearby dry center regions formed by the two Greek key motifs in each domain. The damaged tryptophan residues cause large fluctuations in the Tyr-Trp-Tyr sandwich-like hydrophobic clusters, which in turn break crucial hydrogen-bonds bridging two β-strands in the Greek key motifs at the “tyrosine corner”. Our findings may provide new insights for understanding of the molecular mechanism of the initial stages of UV-induced cataractogenesis.

Interactions between proteins and carbon-based nanoparticles: exploring the origin of nanotoxicity at the molecular level

The widespread application of nanomaterials has spurred an interest in the study of interactions between nanoparticles and proteins due to the biosafety concerns of these nanomaterials. In this review, a summary is presented of some of the recent studies on this important subject, especially on the interactions of proteins with carbon nanotubes (CNTs) and metallofullerenols. Two potential molecular mechanisms have been proposed for CNTs' inhibition of protein functions. The driving forces of CNTs' adsorption onto proteins are found to be mainly hydrophobic interactions and the so-called π-π stacking between CNTs' carbon rings and proteins' aromatic residues. However, there is also recent evidence showing that endohedral metallofullerenol Gd@C82 (OH)22 can be used to inhibit tumor growth, thus acting as a potential nanomedicine. These recent findings have provided a better understanding of nanotoxicity at the molecular level and also suggested therapeutic potential by using nanoparticles' cytotoxicity against cancer cells.

Salts drive controllable multilayered upright assembly of amyloid-like peptides at mica/water interface

Surface-assisted self-assembly of amyloid-like peptides has received considerable interest in both amyloidosis research and nanotechnology in recent years. Despite extensive studies, some controlling factors, such as salts, are still not well understood, even though it is known that some salts can promote peptide self-assemblies through the so-called "salting-out" effect. However, they are usually noncontrollable, disordered, amorphous aggregates. Here, we show via a combined experimental and theoretical approach that a conserved consensus peptide NH2-VGGAVVAGV-CONH2 (GAV-9) (from representative amyloidogenic proteins) can self-assemble into highly ordered, multilayered nanofilaments, with surprising all-upright conformations, under high-salt concentrations. Our atomic force microscopy images also demonstrate that the vertical stacking of multiple layers is highly controllable by tuning the ionic strength, such as from 0 mM (monolayer) to 100 mM (mainly double layer), and to 250 mM MgCl2 (double, triple, quadruple, and quintuple layers). Our atomistic molecular dynamics simulations then reveal that these individual layers have very different internal nanostructures, with parallel β-sheets in the first monolayer but antiparallel β-sheets in the subsequent upper layers due to their different microenvironment. Further studies show that the growth of multilayered, all-upright nanostructures is a common phenomenon for GAV-9 at the mica/water interface, under a variety of salt types and a wide range of salt concentrations.

Hydrophobic interaction drives surface-assisted epitaxial assembly of amyloid-like peptides

The molecular mechanism of epitaxial fibril formation has been investigated for GAV-9 (NH(3)(+)-VGGAVVAGV-CONH(2)), an amyloid-like peptide extracted from a consensus sequence of amyloidogenic proteins, which assembles with very different morphologies, "upright" on mica and "flat" on the highly oriented pyrolytic graphite (HOPG). Our all-atom molecular dynamics simulations reveal that the strong electrostatic interaction induces the "upright" conformation on mica, whereas the hydrophobic interaction favors the "flat" conformation on HOPG. We also show that the epitaxial pattern on mica is ensured by the lattice matching between the anisotropic binding sites of the basal substrate and the molecular dimension of GAV-9, accompanied with a long-range order of well-defined β-strands. Furthermore, the binding free energy surfaces indicate that the longitudinal assembly growth is predominantly driven by the hydrophobic interaction along the longer crystallographic unit cell direction of mica. These findings provide a molecular basis for the surface-assisted molecular assembly, which might also be useful for the design of de novo nanodevices.

Analysis of core region from egg white lysozyme forming amyloid fibrils

Some of the lysozyme mutants in humans cause systemic amyloidosis. Hen egg white lysozyme (HEWL) has been well studied as a model protein of amyloid fibrils formation. We previously identified an amyloid core region consisting of nine amino acids (designated as the K peptide), which is present at 54-62 in HEWL. The K peptide, with tryptophan at its C- terminus, has the ability of self-aggregation. In the present work we focused on its structural properties in relation to the formation of fibrils. The K peptide alone formed definite fibrils having β-sheet structures by incubation of 7 days under acidic conditions at 37°C. A substantial number of fibrils were generated under this pH condition and incubation period. Deletion and substitution of tryptophan in the K peptide resulted in no formation of fibrils. Tryptophan 62 in lysozyme was suggested to be especially crucial to forming amyloid fibrils. We also show that amyloid fibrils formation of the K peptide requires not only tryptophan 62 but also a certain length containing hydrophobic amino acids. A core region is involved in the significant formation of amyloid fibrils of lysozyme.

Collapse of unfolded proteins in a mixture of denaturants

Both urea and guanidinium chloride (GdmCl) are frequently used as protein denaturants. Given that proteins generally adopt extended or unfolded conformations in either aqueous urea or GdmCl, one might expect that the unfolded protein chains will remain or become further extended due to the addition of another denaturant. However, a collapse of denatured proteins is revealed using atomistic molecular dynamics simulations when a mixture of denaturants is used. Both hen egg-white lysozyme and protein L are found to undergo collapse in the denaturant mixture. The collapse of the protein conformational ensembles is accompanied by a decreased solubility and increased non-native self-interactions of hydrophobic residues in the urea/GdmCl mixture. The increase of non-native interactions rather than the native contacts indicates that the proteins experience a simple collapse transition from the fully denatured states. During the protein collapse, the relatively stronger denaturant GdmCl displays a higher tendency to be absorbed onto the protein surface due to their stronger electrostatic interactions with proteins. At the same time, urea molecules also accumulate near the protein surface, resulting in an enhanced "local crowding" for the protein near its first solvation shell. This rearrangement of denaturants near the protein surface and crowded local environment induce the protein collapse, mainly by burying their hydrophobic residues. These findings from molecular simulations are then further explained by a simple analytical model based on statistical mechanics.

A fast parallel clustering algorithm for molecular simulation trajectories

We implemented a GPU-powered parallel k-centers algorithm to perform clustering on the conformations of molecular dynamics (MD) simulations. The algorithm is up to two orders of magnitude faster than the CPU implementation. We tested our algorithm on four protein MD simulation datasets ranging from the small Alanine Dipeptide to a 370-residue Maltose Binding Protein (MBP). It is capable of grouping 250,000 conformations of the MBP into 4000 clusters within 40 seconds. To achieve this, we effectively parallelized the code on the GPU and utilize the triangle inequality of metric spaces. Furthermore, the algorithm's running time is linear with respect to the number of cluster centers. In addition, we found the triangle inequality to be less effective in higher dimensions and provide a mathematical rationale. Finally, using Alanine Dipeptide as an example, we show a strong correlation between cluster populations resulting from the k-centers algorithm and the underlying density. © 2012 Wiley Periodicals, Inc.

Molecular dynamics simulations of Ago silencing complexes reveal a large repertoire of admissible 'seed-less' targets

To better understand the recognition mechanism of RISC and the repertoire of guide-target interactions we introduced G:U wobbles and mismatches at various positions of the microRNA (miRNA) 'seed' region and performed all-atom molecular dynamics simulations of the resulting Ago-miRNA:mRNA ternary complexes. Our simulations reveal that many modifications, including combinations of multiple G:U wobbles and mismatches in the seed region, are admissible and result in only minor structural fluctuations that do not affect overall complex stability. These results are further supported by analyses of HITS-CLIP data. Lastly, introduction of disruptive mutations revealed a bending motion of the PAZ domain along the L1/L2 'hinge' and a subsequent opening of the nucleic-acid-binding channel. Our findings suggest that the spectrum of a miRNA's admissible targets is different from what is currently anticipated by the canonical seed-model. Moreover, they provide a likely explanation for the previously reported sequence-dependent regulation of unintended targeting by siRNAs.

Free-energy simulations reveal that both hydrophobic and polar interactions are important for influenza hemagglutinin antibody binding

Antibodies binding to conserved epitopes can provide a broad range of neutralization to existing influenza subtypes and may also prevent the propagation of potential pandemic viruses by fighting against emerging strands. Here we propose a computational framework to study structural binding patterns and detailed molecular mechanisms of viral surface glycoprotein hemagglutinin (HA) binding with a broad spectrum of neutralizing monoclonal antibody fragments (Fab). We used rigorous free-energy perturbation (FEP) methods to calculate the antigen-antibody binding affinities, with an aggregate underlying molecular-dynamics simulation time of several microseconds (∼2 μs) using all-atom, explicit-solvent models. We achieved a high accuracy in the validation of our FEP protocol against a series of known binding affinities for this complex system, with <0.5 kcal/mol errors on average. We then introduced what to our knowledge are novel mutations into the interfacial region to further study the binding mechanism. We found that the stacking interaction between Trp-21 in HA2 and Phe-55 in the CDR-H2 of Fab is crucial to the antibody-antigen association. A single mutation of either W21A or F55A can cause a binding affinity decrease of ΔΔG > 4.0 kcal/mol (equivalent to an ∼1000-fold increase in the dissociation constant K(d)). Moreover, for group 1 HA subtypes (which include both the H1N1 swine flu and the H5N1 bird flu), the relative binding affinities change only slightly (< ±1 kcal/mol) when nonpolar residues at the αA helix of HA mutate to conservative amino acids of similar size, which explains the broad neutralization capability of antibodies such as F10 and CR6261. Finally, we found that the hydrogen-bonding network between His-38 (in HA1) and Ser-30/Gln-64 (in Fab) is important for preserving the strong binding of Fab against group 1 HAs, whereas the lack of such hydrogen bonds with Asn-38 in most group 2 HAs may be responsible for the escape of antibody neutralization. These large-scale simulations may provide new insight into the antigen-antibody binding mechanism at the atomic level, which could be essential for designing more-effective vaccines for influenza.

Correction of both NBD1 energetics and domain interface is required to restore ΔF508 CFTR folding and function

The folding and misfolding mechanism of multidomain proteins remains poorly understood. Although thermodynamic instability of the first nucleotide-binding domain (NBD1) of ΔF508 CFTR (cystic fibrosis transmembrane conductance regulator) partly accounts for the mutant channel degradation in the endoplasmic reticulum and is considered as a drug target in cystic fibrosis, the link between NBD1 and CFTR misfolding remains unclear. Here, we show that ΔF508 destabilizes NBD1 both thermodynamically and kinetically, but correction of either defect alone is insufficient to restore ΔF508 CFTR biogenesis. Instead, both ΔF508-NBD1 energetic and the NBD1-MSD2 (membrane-spanning domain 2) interface stabilization are required for wild-type-like folding, processing, and transport function, suggesting a synergistic role of NBD1 energetics and topology in CFTR-coupled domain assembly. Identification of distinct structural deficiencies may explain the limited success of ΔF508 CFTR corrector molecules and suggests structure-based combination corrector therapies. These results may serve as a framework for understanding the mechanism of interface mutation in multidomain membrane proteins.

Aggregates with lysozyme and ovalbumin show features of amyloid-like fibrils

The interaction of egg-white lysozyme with N-ovalbumin, the native form of egg-white ovalbumin with the denaturation temperature, Tm, of 78 °C, was investigated by the inhibition of lysozyme muramidase activity, differential scanning calorimetry, and circular dichroism assay as indicators. Signals for the interaction were the most prominent when the mixture of lysozyme and N-ovalbumin was co-heated at 72 °C, slightly lower than the Tm of N-ovalbumin. The interaction was also marked when unheated lysozyme was mixed with N-ovalbumin preheated at 72 °C. Moreover, the mixture rapidly formed fibrous precipitates, which were positive for thioflavin T fluorescent emission, a marker for the amyloid fibril formation. Also electron microscopic observation exhibited features of fibrils. The interaction potency of ovalbumin was ascribed to the tryptic fragment ILELPFASGT MSMLVLLPDE VSGLEQLESIINFEK (residues 229–263), derived from the 2B strands 2 and 3 of ovalbumin. From lysozyme, on the other hand, the chymotryptic peptide RNRCKGTDVQAW (residues 112–123), including cluster 6, and the chymotryptic/tryptic peptide GILQINSRW (residues 54–62), including cluster 3, were responsible for the interaction with N-ovalbumin. Interestingly, this nonamer peptide was found to have the ability to self-aggregate. To the authors knowledge, this may be the first report to document the possible involvement of dual proteins in the formation of amyloid-like fibrils.

Plugging into proteins: poisoning protein function by a hydrophobic nanoparticle

Nanoscale particles have become promising materials in many fields, such as cancer therapeutics, diagnosis, imaging, drug delivery, catalysis, as well as biosensors. In order to stimulate and facilitate these applications, there is an urgent need for the understanding of the nanoparticle toxicity and other risks involved with these nanoparticles to human health. In this study, we use large-scale molecular dynamics simulations to study the interaction between several proteins (WW domains) and carbon nanotubes (one form of hydrophobic nanoparticles). We have found that the carbon nanotube can plug into the hydrophobic core of proteins to form stable complexes. This plugging of nanotubes disrupts and blocks the active sites of WW domains from binding to the corresponding ligands, thus leading to the loss of the original function of the proteins. The key to this observation is the hydrophobic interaction between the nanoparticle and the hydrophobic residues, particularly tryptophans, in the core of the domain. We believe that these findings might provide a novel route to the nanoparticle toxicity on the molecular level for the hydrophobic nanoparticles.

Infrared study of the folding mechanism of a helical hairpin: porcine PYY

The helical hairpin motif plays a key role as a receptor site in DNA binding and protein-protein interactions. Thus, various helical hairpins have recently been developed to assess the factors that control the DNA and/or protein binding affinities of this structural motif and to form synthetic templates for protein and drug design. In addition, several lines of evidence suggest that rapid acquisition of a helical hairpin structure from the unfolded ensemble may guide the rapid formation of helical proteins. Despite its importance as a crucial structural element in protein folding and binding, the folding mechanism of the helical hairpin motif has not been thoroughly studied. Herein, we investigate the structural determinants of the folding kinetics of a naturally occurring helical hairpin (porcine PYY) that is free of disulfide bonds and metal ion-induced cross-links using an infrared temperature-jump technique. It is found that mutations in the turn region predominantly increase the barrier of folding irrespective of the temperature, whereas the effect of mutations that perturb the hydrophobic interactions between the two helices is temperature-dependent. At low temperatures, deletion of hydrophobic side chains is found to predominantly affect the unfolding rate, while the opposite is observed at high temperatures. These results are interpreted in terms of a folding mechanism in which the turn is formed in the transition state and also based on the assumption that cross-strand hydrophobic contacts exist in the thermally unfolded state of PYY.

Single mutation effects on conformational change and membrane deformation of influenza hemagglutinin fusion peptides

The single mutation effect on the conformational change and membrane permeation of influenza hemagglutinin fusion peptides has been studied with molecular dynamics simulations. A total of seven peptides, including wild-type fusion peptide and its six single point mutants (G1E, G1S, G1V, G4V, E11A, and W14A, all with no fusion or hemifusion activity) are examined systematically, which covers a wide range of mutation sites as well as mutant residue types (both hydrophobic and hydrophilic). The wild-type shows a kink structure (inversed V-shape), which facilitates the interaction between the fusion peptide and the lipid bilayer, as well as the interaction between the two arms of the fusion peptide. All mutants show a strong tendency toward a linear alpha-helix conformation, with the initial kink structure in the wild-type broken. More interestingly, one of the key hydrophobic residues around the initial kink region, Phe-9, is found to flip away from the membrane surface in most of these mutants. This conformational change causes a loss of key interactions between the original two arms of the inversed V-shape of the wild-type, thus disabling the kink structure, which results in the stabilization of the linear alpha-helix structure. The fusion peptides also display significant impact on the membrane structure deformation. The thickness of the lipid bilayer surrounding the wild-type fusion peptide decreases significantly, which induces a positive curvature of lipid bilayer. All the single mutations examined here reduce this membrane structural deformation, supporting the fusion activity data from experiments.

Chemistry of Hofmeister anions and osmolytes

The study of the interactions of salts and osmolytes with macromolecules in aqueous solution originated with experiments concerning protein precipitation more than 100 years ago. Today, these solutes are known to display recurring behavior for myriad biological and chemical processes. Such behavior depends both on the nature and concentration of the species in solution. Despite the generality of these effects, our understanding of the molecular-level details of ion and osmolyte specificity is still quite limited. Here, we review recent studies of the interactions between anions and urea with model macromolecular systems. A mechanism for specific ion effects is elucidated for aqueous systems containing charged and uncharged polymers, polypeptides, and proteins. The results clearly show that the effects of the anions are local and involve direct interactions with macromolecules and their first hydration shell. Also, a hydrogen-bonding mechanism is tested for the urea denaturation of proteins with some of these same systems. In that case, direct hydrogen bonding can be largely discounted as the key mechanism for urea stabilization of uncollapsed and/or unfolded structures.

beta-Strand interactions at the domain interface critical for the stability of human lens gammaD-crystallin

Human age-onset cataracts are believed to be caused by the aggregation of partially unfolded or covalently damaged lens crystallin proteins; however, the exact molecular mechanism remains largely unknown. We have used microseconds of molecular dynamics simulations with explicit solvent to investigate the unfolding process of human lens gammaD-crystallin protein and its isolated domains. A partially unfolded folding intermediate of gammaD-crystallin is detected in simulations with its C-terminal domain (C-td) folded and N-terminal domain (N-td) unstructured, in excellent agreement with biochemical experiments. Our simulations strongly indicate that the stability and the folding mechanism of the N-td are regulated by the interdomain interactions, consistent with experimental observations. A hydrophobic folding core was identified within the C-td that is comprised of a and b strands from the Greek key motif 4, the one near the domain interface. Detailed analyses reveal a surprising non-native surface salt-bridge between Glu135 and Arg142 located at the end of the ab folded hairpin turn playing a critical role in stabilizing the folding core. On the other hand, an in silico single E135A substitution that disrupts this non-native Glu135-Arg142 salt-bridge causes significant destabilization to the folding core of the isolated C-td, which, in turn, induces unfolding of the N-td interface. These findings indicate that certain highly conserved charged residues, that is, Glu135 and Arg142, of gammaD-crystallin are crucial for stabilizing its hydrophobic domain interface in native conformation, and disruption of charges on the gammaD-crystallin surface might lead to unfolding and subsequent aggregation.

Free energy simulations reveal a double mutant avian H5N1 virus hemagglutinin with altered receptor binding specificity

Historically, influenza pandemics have been triggered when an avian influenza virus or a human/avian reassorted virus acquires the ability to replicate efficiently and become transmissible in the human population. Most critically, the major surface glycoprotein hemagglutinin (HA) must adapt to the usage of human-like (alpha-2,6-linked) sialylated glycan receptors. Therefore, identification of mutations that can switch the currently circulating H5N1 HA receptor binding specificity from avian to human might provide leads to the emergence of pandemic H5N1 viruses. To define such mutations in the H5 subtype, here we provide a computational framework that combines molecular modeling with extensive free energy simulations. Our results show that the simulated binding affinities are in good agreement with currently available experimental data. Moreover, we predict that one double mutation (V135S and A138S) in HA significantly enhances alpha-2,6-linked receptor recognition by the H5 subtype. Our simulations indicate that this double mutation in H5N1 HA increases the binding affinity to alpha-2,6-linked sialic acid receptors by 2.6 +/- 0.7 kcal/mol per HA monomer that primarily arises from the electrostatic interactions. Further analyses reveal that introduction of this double mutation results in a conformational change in the receptor binding pocket of H5N1 HA. As a result, a major rearrangement occurs in the hydrogen-bonding network of HA with the human receptor, making the human receptor binding pattern of double mutant H5N1 HA surprisingly similar to that observed in human H1N1 HA. These large scale molecular simulations on single and double mutants thus provide new insights into our understanding toward human adaptation of the avian H5N1 virus.

Single mutation induced H3N2 hemagglutinin antibody neutralization: a free energy perturbation study

The single mutation effect on the binding affinity of H3N2 viral protein hemagglutinin (HA) with the monoclonical antibody fragment (Fab) is studied in this paper using the free energy perturbation (FEP) simulations. An all-atom protein model with explicit solvents is used to perform an aggregate of several microsecond FEP molecular dynamics simulations. A recent experiment shows that a single mutation in H3N2 HA, T131I, increases the antibody-antigen dissociation constant Kd by a factor of approximately 4000 (equivalent to a binding affinity decrease of approximately 5 kcal/mol), thus introducing an escape of the antibody (Ab) neutralization. Our FEP result confirms this experimental finding by estimating the HA-Ab binding affinity decrease of 5.2 +/- 0.9 kcal/mol but with a somewhat different molecular mechanism from the experimental findings. Detailed analysis reveals that this large binding affinity decrease in the T131I mutant is mainly due to the displacement of two bridge water molecules otherwise present in the wild-type HA/Ab interface. The decomposition of the binding free energy supports this observation, as the major contribution to the binding affinity is from the electrostatic interactions. In addition, we find that the loss of the binding affinity is also related to the large conformational distortion of one loop (loop 155-161) in the unbound state of the mutant. We then simulate all other possible mutations for this specific mutation site T131, and predict a few more mutations with even larger decreases in the binding affinity (i.e., better candidates for antibody neutralization), such as T131W, T131Y, and T131F. As for further validation, we have also modeled another mutation, S157L, with experimental binding affinity available (Kd increasing approximately 500 times), and found a binding affinity decrease of 4.1 +/- 1.0 kcal/mol, which is again in excellent agreement with experiment. These large scale simulations might provide new insights into the detailed physical interaction, possible future escape mutation, and antibody-antigen coevolution relationship between influenza virus and human antibodies.