Effect of urea on the β-hairpin conformational ensemble and protein denaturation mechanism

Molecular mechanism of the common and opposing cosolvent effects of fluorinated alcohol and urea on a coiled coil protein

Alcohols and urea are widely used as effective protein denaturants. Among monohydric alcohols, 2,2,2-trifluoroethanol (TFE) has large cosolvent effects as a helix stabilizer in proteins. In contrast, urea efficiently denatures ordered native structures, including helices, into coils. These opposing cosolvent effects of TFE and urea are well known, even though both preferentially bind to proteins; however, the underlying molecular mechanism remains controversial. Cosolvent-dependent relative stability between native and denatured states is rigorously related to the difference in preferential binding parameters (PBPs) between these states. In this study, GCN4-p1 with two-stranded coiled coil helices was employed as a model protein, and molecular dynamics simulations for the helix dimer and isolated coil were conducted in aqueous solutions with 2 M TFE and urea. As 2 M cosolvent aqueous solutions did not exhibit clustering of cosolvent molecules, we were able to directly investigate the molecular origin of the excess PBP without considering the enhancement effect of PBPs arising from the concentration fluctuations. The calculated excess PBPs of TFE for the helices and those of urea for the coils were consistent with experimentally observed stabilization of helix by TFE and that of coil by urea. The former was caused by electrostatic interactions between TFE and side chains of the helices, while the latter was attributed to both electrostatic and dispersion interactions between urea and the main chains. Unexpectedly, reverse-micelle-like orientations of TFE molecules strengthened the electrostatic interactions between TFE and the side chains, resulting in strengthening of TFE solvation.

Nanopore single-molecule biosensor in protein denaturation analysis

Unclear issues in protein studies include but not limited to the stability and denaturation mechanism in the presence of denaturants. Herein, we report a dynamic monitoring approach based on nanopore single-molecule biosensor, which can detect the protein's folding and unfolding transitions by recording a nanopore ionic current. When gradually increasing the concentration of denaturant guanidine hydrochloride (GdmCl), sensitive responses were observed with lysozyme unfolding. The emergence of the featured biphasic-pulse demonstrated the existence of a stable intermediate. It was the first time to experimentally confirm the dynamic equilibrium between the intermediate and the native states at single molecule level, therefore consolidating the standpoint of lysozyme denaturation process following the three-state model. Additionally, we got more insights into the conformation about the intermediate as globular-like structure, larger gyration radius, and enhanced positive charge density. We considered that the manner of denaturant toward lysozyme adopts the "direct" model based on stronger electrostatic and van der Waals forces. Nanopore biosensor exhibited excellent sensitivity with a low detection concentration of 280 pM and reproducibility in analysing the folding intermediate of lysozyme.

Does Urea Preferentially Interact with Amide Moieties or Nonpolar Sidechains? A Question Answered Through a Judicious Selection of Model Systems

The transfer model suggests that urea unfolds proteins mainly by increasing the solubility of the amide backbone, probably through urea-induced increase in hydrogen bonding. Other studies suggest that urea addition increases the magnitude of solvent-solute van der Waals interactions, which increases the solubility of nonpolar sidechains. More recent analyses hypothesize that urea has a similar effect in increasing the solubility of backbone and sidechain groups. In this work, we compare the effects of urea addition on the solvation of amides and alkyl groups. At first, we study the effects of urea addition upon solvent hydrogen bonding acidity and basicity through the perturbation in the fluorescence spectrum of probes 1-AN and 1-DMAN. Our results demonstrate that the solvent's hydrogen bonding properties are minimally affected by urea addition. Subsequently, we show that urea addition does not perturb the intra-molecular hydrogen bonding in salicylic acid significantly. Finally, we investigate how urea preferentially interacts with amide and alkyl groups moieties in water by comparing the effects of urea addition upon the solubility of acetaminophen and 4-tertbutylphenol. We show that urea affects amide and t-butyl solubility (lowers the transfer free energy of both amide (backbone) and alkyl (sidechain) groups) in a similar fashion. In other words, preferential interaction of urea with both moieties contributes to protein denaturation.

Energetic, Structural and Dynamic Properties of Nucleobase-Urea Interactions that Aid in Urea Assisted RNA Unfolding

Understanding the structure-function relationships of RNA has become increasingly important given the realization of its functional role in various cellular processes. Chemical denaturation of RNA by urea has been shown to be beneficial in investigating RNA stability and folding. Elucidation of the mechanism of unfolding of RNA by urea is important for understanding the folding pathways. In addition to studying denaturation of RNA in aqueous urea, it is important to understand the nature and strength of interactions of the building blocks of RNA. In this study, a systematic examination of the structural features and energetic factors involving interactions between nucleobases and urea is presented. Results from molecular dynamics (MD) simulations on each of the five DNA/RNA bases in water and eight different concentrations of aqueous urea, and free energy calculations using the thermodynamic integration method are presented. The interaction energies between all the nucleobases with the solvent environment and the transfer free energies become more favorable with respect to increase in the concentration of urea. Preferential interactions of urea versus water molecules with all model systems determined using Kirkwood-Buff integrals and two-domain models indicate preference of urea by nucleobases in comparison to water. The modes of interaction between urea and the nucleobases were analyzed in detail. In addition to the previously identified hydrogen bonding and stacking interactions between urea and nucleobases that stabilize the unfolded states of RNA in aqueous solution, NH-π interactions are proposed to be important. Dynamic properties of each of these three modes of interactions have been presented. The study provides fundamental insights into the nature of interaction of urea molecules with nucleobases and how it disrupts nucleic acids.

Effect of Urea, Arginine, and Ethanol Concentration on Aggregation of 179CVNITV184 Fragment of Sheep Prion Protein

Understanding protein aggregation is of utmost importance as it is responsible for causing several neurodegenerative diseases and one of the serious impediments in large-scale biopharmaceutical production. The prion protein is responsible for pathological states in fatal transmissible spongiform conditions, such as Creutzfeldt-Jakob disease and bovine spongiform encephalopathy. The peptide fragment 178-191 of Syrian hamster prion protein is known to be amyloidogenic. Here, we identified the fragment ₁₇₉CVNITV₁₈₄ as an aggregation-prone fragment in sheep prion protein. This fragment is conserved sequence among sheep and Syrian hamster prion protein and also falls in the previously identified amyloidogenic sequence. The mechanistic details of the aggregation behavior are analyzed in three different concentrations of urea, arginine, and ethanol. Urea and arginine are found to be aggregation suppressors, but ethanol enhances the protein aggregation through β-sheet formation. We have also analyzed the influence of these osmolyte on water dynamics in the presence of the octamer of this aggregation-prone fragment and correlated this water dynamics with the aggregation behavior of the octamer.

Assessment and optimization of collective variables for protein conformational landscape: GB1 β-hairpin as a case study

Collective variables (CVs), when chosen judiciously, can play an important role in recognizing rate-limiting processes and rare events in any biomolecular systems. However, high dimensionality and inherent complexities associated with such biochemical systems render the identification of an optimal CV a challenging task, which in turn precludes the elucidation of an underlying conformational landscape in sufficient details. In this context, a relevant model system is presented by a 16-residue β-hairpin of GB1 protein. Despite being the target of numerous theoretical and computational studies for understanding the protein folding, the set of CVs optimally characterizing the conformational landscape of the β-hairpin of GB1 protein has remained elusive, resulting in a lack of consensus on its folding mechanism. Here we address this by proposing a pair of optimal CVs which can resolve the underlying free energy landscape of the GB1 hairpin quite efficiently. Expressed as a linear combination of a number of traditional CVs, the optimal CV for this system is derived by employing the recently introduced time-structured independent component analysis approach on a large number of independent unbiased simulations. By projecting the replica-exchange simulated trajectories along these pair of optimized CVs, the resulting free energy landscape of this system is able to resolve four distinct well-separated metastable states encompassing the extensive ensembles of folded, unfolded, and molten globule states. Importantly, the optimized CVs were found to be capable of automatically recovering a novel partial helical state of this protein, without needing to explicitly invoke helicity as a constituent CV. Furthermore, a quantitative sensitivity analysis of each constituent in the optimized CV provided key insights on the relative contributions of the constituent CVs in the overall free energy landscapes. Finally, the kinetic pathways connecting these metastable states, constructed using a Markov state model, provide an optimum description of the underlying folding mechanism of the peptide. Taken together, this work offers a quantitatively robust approach toward comprehensive mapping of the underlying folding landscape of a quintessential model system along its optimized CV.

Thermodynamic and structural effect of urea and guanidine chloride on the helical and on a hairpin fragment of GB1 from molecular simulations

With the help of molecular-dynamics simulations, we studied the effect of urea and guanidine chloride on the thermodynamic and structural properties of the helical fragment of protein GB1, comparing them with those of its second beta hairpin. We showed that the helical fragment in different solvents populates an ensemble of states that is more complex than that of the hairpin, and thus the associated experimental observables (circular-dichroism spectra, secondary chemical shifts, m values), that we back-calculated from the simulations and compared with the actual data, are more difficult to interpret. We observed that in the case of both peptides, urea binds tightly to their backbone, while guanidine exerts its denaturing effect in a more subtle way, strongly affecting the electrostatic properties of the solution. This difference can have consequences in the way denaturation experiments are interpreted. Proteins 2017; 85:753-763. © 2016 Wiley Periodicals, Inc.

Comprehensive structural and dynamical view of an unfolded protein from the combination of single-molecule FRET, NMR, and SAXS

The properties of unfolded proteins are essential both for the mechanisms of protein folding and for the function of the large group of intrinsically disordered proteins. However, the detailed structural and dynamical characterization of these highly dynamic and conformationally heterogeneous ensembles has remained challenging. Here we combine and compare three of the leading techniques for the investigation of unfolded proteins, NMR spectroscopy (NMR), small-angle X-ray scattering (SAXS), and single-molecule Förster resonance energy transfer (FRET), with the goal of quantitatively testing their consistency and complementarity and for obtaining a comprehensive view of the unfolded-state ensemble. Using unfolded ubiquitin as a test case, we find that its average dimensions derived from FRET and from structural ensembles calculated using the program X-PLOR-NIH based on NMR and SAXS restraints agree remarkably well; even the shapes of the underlying intramolecular distance distributions are in good agreement, attesting to the reliability of the approaches. The NMR-based results provide a highly sensitive way of quantifying residual structure in the unfolded state. FRET-based nanosecond fluorescence correlation spectroscopy allows long-range distances and chain dynamics to be probed in a time range inaccessible by NMR. The combined techniques thus provide a way of optimally using the complementarity of the available methods for a quantitative structural and dynamical description of unfolded proteins both at the global and the local level.

Contribution of Electrostatics in the Fibril Stability of a Model Ionic-Complementary Peptide

In this work we quantified the role of electrostatic interactions in the self-assembly of a model amphiphilic peptide (RADA 16-I) into fibrillar structures by a combination of size exclusion chromatography and molecular simulations. For the peptide under investigation, it is found that a net charge of +0.75 represents the ideal condition to promote the formation of regular amyloid fibrils. Lower net charges favor the formation of amorphous precipitates, while larger net charges destabilize the fibrillar aggregates and promote a reversible dissociation of monomers from the ends of the fibrils. By quantifying the dependence of the equilibrium constant of this reversible reaction on the pH value and the peptide net charge, we show that electrostatic interactions contribute largely to the free energy of fibril formation. The addition of both salt and a charged destabilizer (guanidinium hydrochloride) at moderate concentration (0.3-1 M) shifts the monomer-fibril equilibrium toward the fibrillar state. Whereas the first effect can be explained by charge screening of electrostatic repulsion only, the promotion of fibril formation in the presence of guanidinium hydrochloride is also attributed to modifications of the peptide conformation. The results of this work indicate that the global peptide net charge is a key property that correlates well with the fibril stability, although the peptide conformation and the surface charge distribution also contribute to the aggregation propensity.

Dispersion interactions between urea and nucleobases contribute to the destabilization of RNA by urea in aqueous solution

Urea has long been used to investigate protein folding and, more recently, RNA folding. Studies have proposed that urea denatures RNA by participating in stacking interactions and hydrogen bonds with nucleic acid bases. In this study, the ability of urea to form unconventional stacking interactions with RNA bases is investigated using ab initio calculations (RI-MP2 and CCSD(T) methods with the aug-cc-pVDZ basis set). A total of 29 stable nucleobase-urea stacked complexes are identified in which the intermolecular interaction energies (up to -14 kcal/mol) are dominated by dispersion effects. Natural bond orbital (NBO) and atoms in molecules (AIM) calculations further confirm strong interactions between urea and nucleobases. Calculations on model systems with multiple urea and water molecules interacting with a guanine base lead to a hypothesis that urea molecules along with water are able to form cage-like structures capable of trapping nucleic acid bases in extrahelical states by forming both hydrogen-bonded and dispersion interactions, thereby contributing to the unfolding of RNA in the presence of urea in aqueous solution.

Computational structural biology of opioid receptors

The publication of high-resolution structures for all of the opioid receptor subfamilies has unveiled exciting opportunities for mechanistic insight into the molecular mechanisms underlying the biology of nociception, reward, and higher cognitive functions, as well as promises for progress in several clinical areas such as pain management, physiological dependence, addiction, and mood disorders. To turn this promise into novel and improved therapeutic entities, however, this information needs to be supplemented with research strategies that explore the dynamic behavior of the proteins and their interactions with other receptors and ligands in their physiological environment.Here we describe state-of-the-art molecular dynamics computational protocols, based on all-atom and coarse-grained modeling techniques, designed to estimate crucial thermodynamic and kinetic parameters describing the binding of small-molecule ligands and the formation of supramolecular complexes.

Tackling sampling challenges in biomolecular simulations

Molecular dynamics (MD) simulations are a powerful tool to give an atomistic insight into the structure and dynamics of proteins. However, the time scales accessible in standard simulations, which often do not match those in which interesting biological processes occur, limit their predictive capabilities. Many advanced sampling techniques have been proposed over the years to overcome this limitation. This chapter focuses on metadynamics, a method based on the introduction of a time-dependent bias potential to accelerate sampling and recover equilibrium properties of a few descriptors that are able to capture the complexity of a process at a coarse-grained level. The theory of metadynamics and its combination with other popular sampling techniques such as the replica exchange method is briefly presented. Practical applications of these techniques to the study of the Trp-Cage miniprotein folding are also illustrated. The examples contain a guide for performing these calculations with PLUMED, a plugin to perform enhanced sampling simulations in combination with many popular MD codes.

P-loop conformation governed crizotinib resistance in G2032R-mutated ROS1 tyrosine kinase: clues from free energy landscape

Tyrosine kinases are regarded as excellent targets for chemical drug therapy of carcinomas. However, under strong purifying selection, drug resistance usually occurs in the cancer cells within a short term. Many cases of drug resistance have been found to be associated with secondary mutations in drug target, which lead to the attenuated drug-target interactions. For example, recently, an acquired secondary mutation, G2032R, has been detected in the drug target, ROS1 tyrosine kinase, from a crizotinib-resistant patient, who responded poorly to crizotinib within a very short therapeutic term. It was supposed that the mutation was located at the solvent front and might hinder the drug binding. However, a different fact could be uncovered by the simulations reported in this study. Here, free energy surfaces were characterized by the drug-target distance and the phosphate-binding loop (P-loop) conformational change of the crizotinib-ROS1 complex through advanced molecular dynamics techniques, and it was revealed that the more rigid P-loop region in the G2032R-mutated ROS1 was primarily responsible for the crizotinib resistance, which on one hand, impaired the binding of crizotinib directly, and on the other hand, shortened the residence time induced by the flattened free energy surface. Therefore, both of the binding affinity and the drug residence time should be emphasized in rational drug design to overcome the kinase resistance.

Cancers can eventually confer drug resistance to the continued medication. In most cases, mutations occurred in a drug target can attenuate the binding affinity of the drugs. Here, we studied the drug resistance mechanisms of the mutations G2032R in the ROS1 tyrosine kinase in fusion-type NSCLC. It is well known that the phosphate-binding loop (P-loop) plays a vital role in the binding of competitive inhibitors in tyrosine kinases, and numerous mutations have been found occurred around the P-loop, which may affect the binding/unbinding process of a drug. Free energy surfaces were constructed to characterize the impact of the mutation to the binding/unbinding process of a well-known NSCLC drug, crizotinib. Two advanced free energy calculation methods, namely funnel based well-tempered metadynamics and umbrella sampling based absolute binding free energy calculation achieved consistent results with the experimental data, suggesting that the rigid P-loop of the mutated target was mainly responsible for the crizotinib resistance to ROS1 tyrosine kinase.

The Parkinson disease-linked LRRK2 protein mutation I2020T stabilizes an active state conformation leading to increased kinase activity

The effect of leucine-rich repeat kinase 2 (LRRK2) mutation I2020T on its kinase activity has been controversial, with both increased and decreased effects being reported. We conducted steady-state and pre-steady-state kinetic studies on LRRKtide and its analog LRRKtide(S). Their phosphorylation differs by the rate-limiting steps: product release is rate-limiting for LRRKtide and phosphoryl transfer is rate-limiting for LRRKtide(S). As a result, we observed that the I2020T mutant is more active than wild type (WT) LRRK2 for LRRKtide(S) phosphorylation, whereas it is less active than WT for LRRKtide phosphorylation. Our pre-steady-state kinetic data suggest that (i) the I2020T mutant accelerates the rates of phosphoryl transfer of both reactions by 3-7-fold; (ii) this increase is masked by a rate-limiting product release step for LRRKtide phosphorylation; and (iii) the observed lower activity of the mutant for LRRKtide phosphorylation is a consequence of its instability: the concentration of the active form of the mutant is 3-fold lower than WT. The I2020T mutant has a dramatically low KATP and therefore leads to resistance to ATP competitive inhibitors. Two well known DFG-out or type II inhibitors are also weaker toward the mutant because they inhibit the mutant in an unexpected ATP competitive mechanism. The I2020 residue lies next to the DYG motif of the activation loop of the LRRK2 kinase domain. Our modeling and metadynamic simulations suggest that the I2020T mutant stabilizes the DYG-in active conformation and creates an unusual allosteric pocket that can bind type II inhibitors but in an ATP competitive fashion.

Effect of urea concentration on aggregation of amyloidogenic hexapeptides (NFGAIL)

We have performed large-scale all-atom molecular dynamics (MD) simulations to study the aggregation behavior of four NFGAIL hexapeptides in the aqueous urea solution, with a urea concentration ranging from 0 to 5 M. We find that urea in general suppresses the peptide aggregation, but suppression slows down in the intermediation concentration regime around 3 M. Two competing mechanisms of urea are determined: urea molecules accumulated near the first solvation shell (FSS) tend to unfold the hexapeptide, which favors aggregation; on the other hand, the tight hydrogen bonds formed between urea and peptide mainchains hinder the association of peptides which disfavors the formation of the β-sheet. Furthermore, the different nonlinear urea concentration dependences of the urea-peptide and peptide-peptide hydrogen bonds lead to a nonmonotonic behavior, with a weak enhancement in the peptide aggregation around 3 M.

Exploring early stages of the chemical unfolding of proteins at the proteome scale

After decades of using urea as denaturant, the kinetic role of this molecule in the unfolding process is still undefined: does urea actively induce protein unfolding or passively stabilize the unfolded state? By analyzing a set of 30 proteins (representative of all native folds) through extensive molecular dynamics simulations in denaturant (using a range of force-fields), we derived robust rules for urea unfolding that are valid at the proteome level. Irrespective of the protein fold, presence or absence of disulphide bridges, and secondary structure composition, urea concentrates in the first solvation shell of quasi-native proteins, but with a density lower than that of the fully unfolded state. The presence of urea does not alter the spontaneous vibration pattern of proteins. In fact, it reduces the magnitude of such vibrations, leading to a counterintuitive slow down of the atomic-motions that opposes unfolding. Urea stickiness and slow diffusion is, however, crucial for unfolding. Long residence urea molecules placed around the hydrophobic core are crucial to stabilize partially open structures generated by thermal fluctuations. Our simulations indicate that although urea does not favor the formation of partially open microstates, it is not a mere spectator of unfolding that simply displaces to the right of the folded←→unfolded equilibrium. On the contrary, urea actively favors unfolding: it selects and stabilizes partially unfolded microstates, slowly driving the protein conformational ensemble far from the native one and also from the conformations sampled during thermal unfolding.

The delicate equilibrium between the folded and functional structure of a protein and its unfolded state is highly dependent on environmental variables such as the solvent. For example the co-solvent urea is a well-known protein denaturant that displaces the equilibrium towards unstructured and non-functional conformations of proteins. However the molecular mechanism behind its ability remains an enigma and the interpretation of the experimental data is still ambiguous. By analyzing a set of representative proteins through extensive molecular dynamics simulations in urea, we provide a robust and consensus picture of the first stages of urea-driven protein unfolding and elucidate the role of urea in accelerating protein unfolding. Our results suggest that urea, thanks to its stickiness and slow diffusion, benefits from the intrinsic flexibility of proteins and stabilizes partially open-states, slowly driving the protein toward unfolding.

Free-energy landscape of protein oligomerization from atomistic simulations

In the realm of protein-protein interactions, the assembly process of homooligomers plays a fundamental role because the majority of proteins fall into this category. A comprehensive understanding of this multistep process requires the characterization of the driving molecular interactions and the transient intermediate species. The latter are often short-lived and thus remain elusive to most experimental investigations. Molecular simulations provide a unique tool to shed light onto these complex processes complementing experimental data. Here we combine advanced sampling techniques, such as metadynamics and parallel tempering, to characterize the oligomerization landscape of fibritin foldon domain. This system is an evolutionarily optimized trimerization motif that represents an ideal model for experimental and computational mechanistic studies. Our results are fully consistent with previous experimental nuclear magnetic resonance and kinetic data, but they provide a unique insight into fibritin foldon assembly. In particular, our simulations unveil the role of nonspecific interactions and suggest that an interplay between thermodynamic bias toward native structure and residual conformational disorder may provide a kinetic advantage.

Structure and dynamics of urea/water mixtures investigated by vibrational spectroscopy and molecular dynamics simulation

Urea/water is an archetypical "biological" mixture and is especially well-known for its relevance to protein thermodynamics as urea acts as a protein denaturant at high concentration. This behavior has given rise to an extended debate concerning urea's influence on water structure. On the basis of a variety of methods and of definitions of the water structure, urea has been variously described as a structure-breaker, a structure-maker, or as remarkably neutral toward water. Because of its sensitivity to microscopic structure and dynamics, vibrational spectroscopy can help resolve these debates. We report experimental and theoretical spectroscopic results for the OD stretch of HOD/H2O/urea mixtures (linear IR, 2DIR, and pump-probe anisotropy decay) and for the CO stretch of urea-D4/D2O mixtures (linear IR only). Theoretical results are obtained using existing approaches for water and a modification of a frequency map developed for acetamide. All absorption spectra are remarkably insensitive to urea concentration, consistent with the idea that urea only very weakly perturbs the water structure. Both this work and experiments by Rezus and Bakker, however, show that water's rotational dynamics are slowed down by urea. Analysis of the simulations casts doubt on the suggestion that urea immobilizes particular doubly hydrogen bonded water molecules.

Molecular Insight into Different Denaturing Efficiency of Urea, Guanidinium, and Methanol: A Comparative Simulation Study

We have designed various nanoslit systems, whose opposing surfaces can be either hydrophobic, hydrophilic, or simply a water-vapor interface, for the molecular dynamics simulation of confined water with three different protein denaturants, i.e., urea, guanidinium chloride (GdmCl), and methanol, respectively. Particular attention is placed on the preferential adsorption of the denaturant molecules onto the opposing surfaces and associated resident time in the vicinal layer next to the surfaces, as well as their implication in the denaturing efficiency of different denaturant molecules. Our simulation results show that among the three denaturants, the occupancy of methanol in the vicinal layer is the highest while the residence time of Gdm is the longest. Although the occupancy and the residence time of urea in the vicinal layer is less than those of the other two denaturant molecules, urea entails "all-around" properties for being a highly effective denaturant. The distinct characteristics of three denaturants may suggest a different molecular mechanism for the protein denaturation. This comparative simulation by design allows us to gain additional insights, on the molecular level, into the denaturation effect and related hydrophobic effect.

Toward an atomistic description of the urea-denatured state of proteins

We present here the characterization of the structural, dynamics, and energetics of properties of the urea-denatured state of ubiquitin, a small prototypical soluble protein. By combining state-of-the-art molecular dynamics simulations with NMR and small-angle X-ray scattering data, we were able to: (i) define the unfolded state ensemble, (ii) understand the energetics stabilizing unfolded structures in urea, (iii) describe the dedifferential nature of the interactions of the fully unfolded proteins with urea and water, and (iv) characterize the early stages of protein refolding when chemically denatured proteins are transferred to native conditions. The results presented herein are unique in providing a complete picture of the chemically unfolded state of proteins and contribute to deciphering the mechanisms that stabilize the native state of proteins, as well as those that maintain them unfolded in the presence of urea.

Type II kinase inhibitors show an unexpected inhibition mode against Parkinson's disease-linked LRRK2 mutant G2019S

A number of well-known type II inhibitors (ATP-noncompetitive) that bind kinases in their DFG-out conformation were tested against wild-type LRRK2 and the most common Parkinson's disease-linked mutation, G2019S. We found that traditional type II inhibitors exhibit surprising variability in their inhibition mechanism between the wild type (WT) and the G2019S mutant of LRRK2. The type II kinase inhibitors were found to work in an ATP-competitive fashion against the G2019S mutant, whereas they appear to follow the expected noncompetitive mechanism against WT. Because the G2019S mutation lies in the DXG motif (DYG in LRRK2 but DFG in most other kinases) of the activation loop, we explored the structural consequence of the mutation on loop dynamics using an enhanced sampling method called metadynamics. The simulations suggest that the G2019S mutation stabilizes the DYG-in state of LRRK2 through a series of hydrogen bonds, leading to an increase in the conformational barrier between the active and inactive forms of the enzyme and a relative stabilization of the active form. The conformational bias toward the active form of LRRK2 mutants has two primary consequences. (1) The mutant enzyme becomes hyperactive, a known contributor to the Parkinsonian phenotype, as a consequence of being "locked" into the activated state, and (2) the mutation creates an unusual allosteric pocket that can bind type II inhibitors but in an ATP-competitive fashion. Our results suggest that developing type II inhibitors, which are generally considered superior to type I inhibitors because of desirable selectivity profiles, might be especially challenging for the G2019S LRRK2 mutant.

Quantitative characterization of local protein solvation to predict solvent effects on protein structure

Characterization of solvent preferences of proteins is essential to the understanding of solvent effects on protein structure and stability. Although it is generally believed that solvent preferences at distinct loci of a protein surface may differ, quantitative characterization of local protein solvation has remained elusive. In this study, we show that local solvation preferences can be quantified over the entire protein surface from extended molecular dynamics simulations. By subjecting microsecond trajectories of two proteins (lysozyme and antibody fragment D1.3) in 4 M glycerol to rigorous statistical analyses, solvent preferences of individual protein residues are quantified by local preferential interaction coefficients. Local solvent preferences for glycerol vary widely from residue to residue and may change as a result of protein side-chain motions that are slower than the longest intrinsic solvation timescale of ∼10 ns. Differences of local solvent preferences between distinct protein side-chain conformations predict solvent effects on local protein structure in good agreement with experiment. This study extends the application scope of preferential interaction theory and enables molecular understanding of solvent effects on protein structure through comprehensive characterization of local protein solvation.

Folding Kinetics and Unfolded State Dynamics of the GB1 Hairpin from Molecular Simulation

The C-terminal β-hairpin of protein G is a 16-residue peptide that folds in a two-state fashion akin to many larger proteins. However, with an experimental folding time of ∼6 μs, it remains a challenging system for all-atom, explicitly solvated, molecular dynamics simulations. Here, we use a large simulation data set (0.7 ms total) of the hairpin at 300 and 350 K to interpret its folding via a master equation approach. We find a separation of over an order of magnitude between the longest and second longest relaxation times, with the slowest relaxation corresponding to folding. However, in spite of this apparent two-state dynamics, the folding rate determined based on a first-passage time analysis depends on the initial conditions chosen, with a nonexponential distribution of first passage times being obtained in some cases. Using the master equation model, we are now able to account quantitatively for the observed distribution of first passage times. The deviation from the expected exponential distribution for a two-state system arises from slow dynamics in the unfolded state, associated with formation and melting of helical structures. Our results help to reconcile recent findings of slow dynamics in unfolded proteins with observed two-state folding kinetics. At the same time, they indicate that care is required in estimating folding kinetics from many short folding simulations. Last, we are able to use the master equation model to obtain details of the folding mechanism and folding transition state, which appear consistent with the "zipper" mechanism inferred from the experiment.

Cosolvent effects on protein stability

Proteins are marginally stable, and the folding/unfolding equilibrium of proteins in aqueous solution can easily be altered by the addition of small organic molecules known as cosolvents. Cosolvents that shift the equilibrium toward the unfolded ensemble are termed denaturants, whereas those that favor the folded ensemble are known as protecting osmolytes. Urea is a widely used denaturant in protein folding studies, and the molecular mechanism of its action has been vigorously debated in the literature. Here we review recent experimental as well as computational studies that show an emerging consensus in this problem. Urea has been shown to denature proteins through a direct mechanism, by interacting favorably with the peptide backbone as well as the amino acid side chains. In contrast, the molecular mechanism by which the naturally occurring protecting osmolyte trimethylamine N-oxide (TMAO) stabilizes proteins is not clear. Recent studies have established the strong interaction of TMAO with water. Detailed molecular simulations, when used with force fields that incorporate these interactions, can provide insight into this problem. We present the development of a model for TMAO that is consistent with experimental observations and that provides physical insight into the role of cosolvent-cosolvent interaction in determining its preferential interaction with proteins.

Urea-Water Solvation Forces on Prion Structures

Solvation forces are crucial determinants in the equilibrium between the folded and unfolded state of proteins. Particularly interesting are the solvent forces of denaturing solvent mixtures on folded and misfolded states of proteins involved in neurodegeneration. The C-terminal globular domain of the ovine prion protein (1UW3) and its analogue H2H3 in the α-rich and β-rich conformation were used as model structures to study the solvation forces in 4 M aqueous urea using molecular dynamics. The model structures display very different secondary structures and solvent exposures. Most protein atoms favor interactions with urea over interactions with water. The force difference between protein-urea and protein-water interactions correlates with hydrophobicity; i.e., urea interacts preferentially with hydrophobic atoms, in agreement with results from solvent transfer experiments. Solvent Shannon entropy maps illustrate the mobility gradient of the urea-water mixture from the first solvation shell to the bulk. Single urea molecules replace water in the first solvation shell preferably at locations of relatively high solvent entropy.

Assessing the relative stability of dimer interfaces in g protein-coupled receptors

Considerable evidence has accumulated in recent years suggesting that G protein-coupled receptors (GPCRs) associate in the plasma membrane to form homo- and/or heteromers. Nevertheless, the stoichiometry, fraction and lifetime of such receptor complexes in living cells remain topics of intense debate. Motivated by experimental data suggesting differing stabilities for homomers of the cognate human β1- and β2-adrenergic receptors, we have carried out approximately 160 microseconds of biased molecular dynamics simulations to calculate the dimerization free energy of crystal structure-based models of these receptors, interacting at two interfaces that have often been implicated in GPCR association under physiological conditions. Specifically, results are presented for simulations of coarse-grained (MARTINI-based) and atomistic representations of each receptor, in homodimeric configurations with either transmembrane helices TM1/H8 or TM4/3 at the interface, in an explicit lipid bilayer. Our results support a definite contribution to the relative stability of GPCR dimers from both interface sequence and configuration. We conclude that β1- and β2-adrenergic receptor homodimers with TM1/H8 at the interface are more stable than those involving TM4/3, and that this might be reconciled with experimental studies by considering a model of oligomerization in which more stable TM1 homodimers diffuse through the membrane, transiently interacting with other protomers at interfaces involving other TM helices.