The Role of Electrostatic Interactions in Folding of β-Proteins

A Hydroxylamine-Mediated Amidination of Lysine Residues That Retains the Protein's Positive Charge

Lysine-specific peptide and protein modification strategies are widely used to study charge-related functions and applications. However, these strategies often result in the loss of the positive charge on lysine, significantly impacting the charge-related properties of proteins. Herein, we report a strategy to preserve the positive charge and selectively convert amines in lysine side chains to amidines using nitriles and hydroxylamine under aqueous conditions. Various unprotected peptides and proteins were successfully modified with a high conversion rate. Moreover, the reactive amidine moiety and derived modification site enable subsequent secondary modifications. Notably, positive charges were retained during the modification. Therefore, positive charge-related protein properties, such as liquid-liquid phase separation behaviour of α-synuclein, were not affected. This strategy was subsequently applied to a lysine rich protein to develop an amidine-containing coacervate DNA complex with outstanding mechanical properties. Overall, our innovative strategy provides a new avenue to explore the characteristics of positively charged proteins.

Improved Random Batch Ewald Method in Molecular Dynamics Simulations

The random batch Ewald (RBE) is an efficient and accurate method for molecular dynamics (MD) simulations of physical systems at the nano/microscale. The method shows great potential to solve the computational bottleneck of long-range interactions, motivating a necessity to accelerate short-range components of the nonbonded interactions for a further speedup of MD simulations. In this work, we present an improved RBE method for the nonbonding interactions by introducing the random batch idea to constructing neighbor lists for the treatment of both the short-range part of the Ewald splitting and the Lennard-Jones potential. The efficiency that the novel neighbor list algorithm owes to the stochastic minibatch strategy can significantly reduce the total number of neighbors. We obtain the error estimate and convergence by theoretical analysis and implement the improved RBE method in the LAMMPS package. Benchmark simulations are performed to demonstrate the accuracy and stability of the algorithm. Numerical tests on computer performance by conducting large-scaled MD simulations for systems including up to 0.1 billion water molecules are run on massive clusters with up to 50000 CPU cores, demonstrating the attractive features such as the high parallel scalability and memory-saving of the method in comparison to the existing methods.

Protein dynamics of human serum albumin at hypothermic temperatures investigated by temperature jump

Human serum albumin (HSA) is the most abundant protein in human plasma. Most protein dynamics studies of HSA have been performed above the hyperthermia temperature (>42 °C), so information on the dynamics under hypothermic conditions (<35 °C) is lacking. In this work, a tryptophan-based fluorescence temperature jump system was employed to investigate the thermally-induced dynamic process of HSA at a physiological concentration of ca. 45 mg mL^-1 and pH = ca. 7 upon an instantaneous temperature increase from 25 °C to 30-43 °C. The observed kinetics manifested a three-state consecutive feature, . Upon analysis with the Arrhenius model, the rate coefficients k₁ and k₂ manifested piecewise temperature dependence, and the turning-point temperature was found to be ca. 34 °C, coinciding with the upper bound of hypothermic temperature. Meanwhile, the corresponding activation energies of the transitions at 34-43 °C were lower than those at 30-34 °C, suggesting that protein conformational adjustments at 34-43 °C were more feasible than those at hypothermic temperatures. These observations provided a fresh viewpoint on the relationship between the energetics of protein dynamics and the apparent functioning of a given protein at the molecular level.

Sampling of the conformational landscape of small proteins with Monte Carlo methods

Computer simulation provides an increasingly realistic picture of large-scale conformational change of proteins, but investigations remain fundamentally constrained by the femtosecond timestep of molecular dynamics simulations. For this reason, many biologically interesting questions cannot be addressed using accessible state-of-the-art computational resources. Here, we report the development of an all-atom Monte Carlo approach that permits the modelling of the large-scale conformational change of proteins using standard off-the-shelf computational hardware and standard all-atom force fields. We demonstrate extensive thermodynamic characterization of the folding process of the α-helical Trp-cage, the Villin headpiece and the β-sheet WW-domain. We fully characterize the free energy landscape, transition states, energy barriers between different states, and the per-residue stability of individual amino acids over a wide temperature range. We demonstrate that a state-of-the-art intramolecular force field can be combined with an implicit solvent model to obtain a high quality of the folded structures and also discuss limitations that still remain.

Nanostructured Ni-Cu Electrocatalysts for the Oxygen Evolution Reaction

Ni-based materials are promising electrocatalysts for the oxygen evolution reaction (OER) for water splitting in alkaline media. We report the synthesis and OER electrocatalysis of both Ni-Cu nanoparticles (20-50 nm in diameter) and Ni-Cu nanoclusters (<20 metal atoms). Analysis of mass spectral data from matrix-assisted laser desorption/ionization and electrospray ionization techniques demonstrates that discrete heterobimetallic Ni-Cu nanoclusters capped with glutathione ligands were successfully synthesized. Ni-Cu nanoclusters with a 52:48 mol % Ni:Cu metal composition display an OER onset overpotential of 50 mV and an overpotential of 150 mV at 10 mA cm-2, which makes this catalyst one of the most efficient nonprecious metal OER catalysts. The durability of the nanocluster catalysts on carbon electrodes can be extended by appending them to electrodes modified with TiO2 nanoparticles. Infrared spectroscopy results indicate that the aggregation dynamics of the glutathione ligands change during catalysis. Taken together, these results help explain the reactivity of a novel class of nanostructured Ni-Cu OER catalysts, which are underexplored alternatives to more commonly studied Ni-Fe, Ni-Co, and Ni-Mn materials.

Folding and Unfolding of the Short Light-Triggered β-Hairpin Peptide AzoChignolin Occurs within 100 ns

To map the underlying molecular mechanisms of folding dynamics in proteins, light-operated peptides have emerged as promising tools. In this study, we reveal the complete sequence of light-induced structural changes of AzoChignolin, a short β-hairpin peptide containing an azobenzene photoswitch in its loop region. Light-triggered structural changes were monitored by time-resolved IR spectroscopy. Formation and destruction of the hairpin structure is very fast and occurs within 100 ns for AzoChignolin in methanol. Atomistic molecular dynamics simulations using two explicit solvents, methanol and water, revealed the underlying molecular processes and allowed us to gain further insight into the reaction mechanism. Despite its rapid reaction time, hairpin formation in these solvents is not force-driven by the molecular switch but proceeded via formation of interstrand hydrogen bonds and contacts between aromatic residues. Moreover, the combined experimental and theoretical study demonstrates that the solvent (methanol vs water) does not dictate the velocity of β-hairpin formation in the AzoChignolin peptide comprising only a few hydrophobic residues in the strands.

Peptide Side-COOH Groups Have Two Distinct Conformations under Biorelevant Conditions

The carboxyl (COOH) side chain groups of amino acids, such as aspartic acid, play an important role in biochemical processes, including enzymatic proton transport. In many theoretical studies, it was found that the (bio)chemical reactivity of the carboxyl group strongly depends on the conformation of this group. Interestingly, up to now there has been no experimental investigation of the geometry and the stability of different COOH conformers under biorelevant conditions. Here, we investigate the conformational isomerism of the side chain COOH group of N-acetyl aspartic acid amide using polarization-resolved two-dimensional infrared spectroscopy. We find that the carboxyl group shows two distinct near-planar conformers (syn and anti) when dissolved in water at room temperature. Both conformers are significantly populated in aqueous solution (75 ± 10% and 25 ± 10% for syn and anti, respectively). Molecular dynamics simulations show that the anti conformer interacts more strongly with water molecules than the syn conformer, explaining why this conformer is significantly present in aqueous solution.

A quantitative connection of experimental and simulated folding landscapes by vibrational spectroscopy

We break the barrier between simulation and experiment by comparing identical computed and experimental infrared observables.

For small molecule reaction kinetics, computed reaction coordinates often mimic experimentally measured observables quite accurately. Although nowadays simulated and measured biomolecule kinetics can be compared on the same time scale, a gap between computed and experimental observables remains. Here we directly compared temperature-jump experiments and molecular dynamics simulations of protein folding dynamics using the same observable: the time-dependent infrared spectrum. We first measured the stability and folding kinetics of the fastest-folding β-protein, the GTT35 WW domain, using its structurally specific infrared spectrum. The relaxation dynamics of the peptide backbone, β-sheets, turn, and random coil were measured independently by probing the amide I′ region at different frequencies. Next, the amide I′ spectra along folding/unfolding molecular dynamics trajectories were simulated by accurate mixed quantum/classical calculations. The simulated time dependence and spectral amplitudes at the exact experimental probe frequencies provided relaxation and folding rates in agreement with experimental observations. The calculations validated by experiment yield direct structural evidence for a rate-limiting reaction step where an intermediate state with either the first or second hairpin is formed. We show how folding switches from a more homogeneous (apparent two-state) process at high temperature to a more heterogeneous process at low temperature, where different parts of the WW domain fold at different rates.

Folding free energy landscapes of β-sheets with non-polarizable and polarizable CHARMM force fields

Molecular dynamics (MD) simulations of peptides and proteins offer atomic-level detail into many biological processes, although the degree of insight depends on the accuracy of the force fields used to represent them. Protein folding is a key example in which the accurate reproduction of folded-state conformations of proteins and kinetics of the folding processes in simulation is a longstanding goal. Although there have been a number of recent successes, challenges remain in capturing the full complexity of folding for even secondary-structure elements. In the present work, we have used all-atom MD simulations to study the folding properties of one such element, the C-terminal β-hairpin of the B1 domain of streptococcal protein G (GB1). Using replica-exchange umbrella sampling simulations, we examined the folding free energy of two fixed-charge CHARMM force fields, CHARMM36 and CHARMM22^*, as well as a polarizable force field, the CHARMM Drude-2013 model, which has previously been shown to improve the folding properties of α-helical peptides. The CHARMM22^* and Drude-2013 models are in rough agreement with experimental studies of GB1 folding, while CHARMM36 overstabilizes the β-hairpin. Additional free-energy calculations show that small adjustments to the atomic polarizabilities in the Drude-2013 model can improve both the backbone solubility and folding properties of GB1 without significantly affecting the model's ability to properly fold α-helices. We also identify a non-native salt bridge in the β-turn region that overstabilizes the β-hairpin in the C36 model. Finally, we demonstrate that tryptophan fluorescence is insufficient for capturing the full β-hairpin folding pathway.

Soluble Zwitterionic Poly(sulfobetaine) Destabilizes Proteins

The widespread interest in neutral, water-soluble polymers such as poly(ethylene glycol) (PEG) and poly(zwitterions) such as poly(sulfobetaine) (pSB) for biomedical applications is due to their widely assumed low protein binding. Here we demonstrate that pSB chains in solution can interact with proteins directly. Moreover, pSB can reduce the thermal stability and increase the protein folding cooperativity relative to proteins in buffer or in PEG solutions. Polymer-dependent changes in the tryptophan fluorescence spectra of three structurally-distinct proteins reveal that soluble, 100 kDa pSB interacts directly with all three proteins and changes both the local polarity near tryptophan residues and the protein conformation. Thermal denaturation studies show that the protein melting temperatures decrease by as much as ∼1.9 °C per weight percent of polymer and that protein folding cooperativity increases by as much as ∼130 J mol^-1 K^-1 per weight percent of polymer. The exact extent of the changes is protein-dependent, as some proteins exhibit increased stability, whereas others experience decreased stability at high soluble pSB concentrations. These results suggest that pSB is not universally protein-repellent and that its efficacy in biotechnological applications will depend on the specific proteins used.

Protein Folding Cooperativity and Thermodynamic Barriers of the Simplest β-Sheet Fold: A Survey of WW Domains

Theory and experiments have shown that microsecond folding proteins exhibit characteristic thermodynamic properties that reflect the limited cooperativity of folding over marginal barriers (downhill folding). Those studies have mostly focused on proteins with large α-helical contents and small size, which tend to be the fastest folders. A key open question is whether such properties are also present in the fastest all-β proteins. We address this issue by investigating the unfolding thermodynamics of a collection of WW domains as representatives of the simplest β-sheet fold. WW domains are small microsecond folders, although they do not fold as fast as their α-helical counterparts. In previous work on the NEDD4-WW4 domain, we reported deviations from two-state thermodynamics that were less apparent and thus suggestive of an incipient downhill scenario. Here we investigate the unfolding thermodynamics of four other WW domains (NEDD4-WW3, YAP65-WW1(L30K), FBP11-WW1, and FBP11-WW2) by performing all of the thermodynamic tests for downhill folding that have been previously developed on α-helical proteins. This set of five WW domains shares low sequence identity and include examples from two specificity classes, thus providing a comprehensive survey. Thermodynamic analysis of the four new WW domains consistently reveals all of the properties of downhill folding equilibria, which are in all cases more marked than what we found before in NEDD4-WW4. Our results show that fast-folding all-β proteins do share limited cooperativity and gradual unfolding thermodynamics with fast α-helical proteins and suggest that the free energy barrier to folding of natural proteins is mostly determined by size and fold topology and much less by the specific amino acid sequence.

Parallel folding pathways of Fip35 WW domain explained by infrared spectra and their computer simulation

We present a calculation of the amide I' infrared (IR) spectra of the folded, unfolded, and intermediate states of the WW domain Fip35, a model system for β-sheet folding. Using an all-atom molecular dynamics simulation in which multiple folding and unfolding events take place we identify six conformational states and then apply perturbed matrix method quantum-mechanical calculations to determine their amide I' IR spectra. Our analysis focuses on two states previously identified as Fip35 folding intermediates and suggests that a three-stranded core similar to the folded state core is the main source of the spectroscopic differences between the two intermediates. In particular, we propose a hypothesis for why folding via one of these intermediates was not experimentally observed by IR T-jump.

Meandering Down the Energy Landscape of Protein Folding: Are We There Yet?

As judged by a single publication metric, the activity in the protein folding field has been declining over the past 5 years, after enjoying a decade-long growth. Does this development indicate that the field is sunsetting or is this decline only temporary? Upon surveying a small territory of its landscape, we find that the protein folding field is still quite active and many important findings have emerged from recent experimental studies. However, it is also clear that only continued development of new techniques and methods, especially those enabling dissection of the fine details and features of the protein folding energy landscape, will fuel this old field to move forward.

Effects of naturally occurring charged mutations on the structure, stability, and binding of the Pin1 WW domain

Pin1 is a peptidyl-prolyl cis-trans isomerase, whose WW domain specifically recognizes the pSer/Thr-Pro motif. Pin1 is involved in multiple phosphorylation events that regulate the activities of various substrates, and Pin1 deregulation has been reported in various diseases, including cancer and Alzheimer's disease. The WW domain of Pin1 has been used as a small model protein to investigate the folding mechanisms of the β-sheet structure by studying the effect of mutations or its naturally occurring variants. However, only a few studies have investigated the structure and binding of Pin1 WW mutants. In the present work, two naturally occurring Pin1 WW variants, namely, G20D and S16R, derived from the cynomolgus monkey and African green monkey, respectively, were selected to investigate the influence of charge mutation on the structure, stability, and binding properties of the Pin1 WW domain. Analysis using a combination of nuclear magnetic resonance (NMR) and chemical shift-based calculations revealed that the G20D and S16R mutants had high structural similarity to the wild-type Pin1 WW domain. However, the presence of a charge mutation significantly decreased the stability of the Pin1 WW domain. Both the wild-type and G20D forms of the Pin1 WW domain utilized a three-site mode to bind to a phosphorylated Tau peptide, pT231, whereas the S16R mutant binds to the pT231 peptide either in a non-specific manner or through a totally different binding mechanism. Correspondingly, the wild-type and two mutant Pin1 WW domains showed different binding affinities to the Tau phosphopeptide. Considering that the WW domain participates in the catalytic activity of the Pin1 isomerase, our study represents a novel approach for studying Pin1 function through the analysis of its naturally occurring mutants.

Dual time-resolved temperature-jump fluorescence and infrared spectroscopy for the study of fast protein dynamics

Time-resolved temperature-jump (T-jump) coupled with fluorescence and infrared (IR) spectroscopy is a powerful technique for monitoring protein dynamics. Although IR spectroscopy of the polypeptide amide I mode is more technically challenging, it offers complementary information because it directly probes changes in the protein backbone, whereas, fluorescence spectroscopy is sensitive to the environment of specific side chains. With the advent of widely tunable quantum cascade lasers (QCL) it is possible to efficiently probe multiple IR frequencies with high sensitivity and reproducibility. Here we describe a dual time-resolved T-jump fluorescence and IR spectrometer and its application to study protein folding dynamics. A Q-switched Ho:YAG laser provides the T-jump source for both time-resolved IR and fluorescence spectroscopy, which are probed by a QCL and Ti:Sapphire laser, respectively. The Ho:YAG laser simultaneously pumps the time-resolved IR and fluorescence spectrometers. The instrument has high sensitivity, with an IR absorbance detection limit of <0.2mOD and a fluorescence sensitivity of 2% of the overall fluorescence intensity. Using a computer controlled QCL to rapidly tune the IR frequency it is possible to create a T-jump induced difference spectrum from 50ns to 0.5ms. This study demonstrates the power of the dual time-resolved T-jump fluorescence and IR spectroscopy to resolve complex folding mechanisms by complementary IR absorbance and fluorescence measurements of protein dynamics.

The Dependence of Carbohydrate-Aromatic Interaction Strengths on the Structure of the Carbohydrate

Interactions between proteins and carbohydrates are ubiquitous in biology. Therefore, understanding the factors that determine their affinity and selectivity are correspondingly important. Herein, we have determined the relative strengths of intramolecular interactions between a series of monosaccharides and an aromatic ring close to the glycosylation site in an N-glycoprotein host. We employed the enhanced aromatic sequon, a structural motif found in the reverse turns of some N-glycoproteins, to facilitate face-to-face monosaccharide-aromatic interactions. A protein host was used because the dependence of the folding energetics on the identity of the monosaccharide can be accurately measured to assess the strength of the carbohydrate-aromatic interaction. Our data demonstrate that the carbohydrate-aromatic interaction strengths are moderately affected by changes in the stereochemistry and identity of the substituents on the pyranose rings of the sugars. Galactose seems to make the weakest and allose the strongest sugar-aromatic interactions, with glucose, N-acetylglucosamine (GlcNAc) and mannose in between. The NMR solution structures of several of the monosaccharide-containing N-glycoproteins were solved to further understand the origins of the similarities and differences between the monosaccharide-aromatic interaction energies. Peracetylation of the monosaccharides substantially increases the strength of the sugar-aromatic interaction in the context of our N-glycoprotein host. Finally, we discuss our results in light of recent literature regarding the contribution of electrostatics to CH-π interactions and speculate on what our observations imply about the absolute conservation of GlcNAc as the monosaccharide through which N-linked glycans are attached to glycoproteins in eukaryotes.