Trp-Cage Folding on Organic Surfaces

Molecular View on the iRGD Peptide Binding Mechanism: Implications for Integrin Activity and Selectivity Profiles

Receptor-selective peptides are widely used as smart carriers for specific tumor-targeted delivery. A remarkable example is the cyclic nonapeptide iRGD (CRGDKPGDC, 1) that couples intrinsic cytotoxic effects with striking tumor-homing properties. These peculiar features are based on a rather complex multistep mechanism of action, where the primary event is the recognition of RGD integrins. Despite the high number of preclinical studies and the recent success of a phase I trial for the treatment of pancreatic ductal adenocarcinoma (PDAC), there is little information available about the iRGD three-dimensional (3D) structure and integrin binding properties. Here, we re-evaluate the peptide's affinity for cancer-related integrins including not only the previously known targets αvβ3 and αvβ5 but also the αvβ6 isoform, which is known to drive cell growth, migration, and invasion in many malignancies including PDAC. Furthermore, we use parallel tempering in the well-tempered ensemble (PT-WTE) metadynamics simulations to characterize the in-solution conformation of iRGD and extensive molecular dynamics calculations to fully investigate its binding mechanism to integrin partners. Finally, we provide clues for fine-tuning the peptide's potency and selectivity profile, which, in turn, may further improve its tumor-homing properties.

A Transfer Free Energy Based Implicit Solvent Model for Protein Simulations in Solvent Mixtures: Urea-Induced Denaturation as a Case Study

We developed a method for implicit solvent molecular dynamics simulations of proteins in solvent mixtures (model with implicit solvation thermodynamics, MIST). The MIST method introduces experimental group transfer free energies to the generalized Born formulation for generating molecular trajectories without the need for developing rigorous explicit-solvent force fields for multicomponent solutions. As a test case, we studied the urea-induced denaturation of the Trp-cage miniprotein in water. We demonstrate that our method allows efficient exploration of the conformational space of the protein in only a few hundreds of nanoseconds of all-atom unbiased simulations. Furthermore, selective implementation of the transfer free energies of specific peptide groups, backbone, and side chains enables us to decouple their specific energetic contributions to the conformational changes of the protein. The approach herein developed can readily be extended to the investigation of complex matrices as well as to the characterization of protein aggregation. The MIST method is implemented in Plumed (ver. 2.8) as a separate module called SASA.

The effect of charged residue substitutions on the thermodynamics of protein-surface interactions

The interactions of proteins with surfaces are important in both biological processes and biotechnologies. In contrast to decades of study regarding the biophysics of proteins in bulk solution, however, our mechanistic understanding of the biophysics of proteins interacting with surfaces remains largely qualitative. In response, we have set to explore quantitatively the thermodynamics of protein-surface interactions. In this work, we explore systematically the role of electrostatics in modulating the interaction between proteins and charged surfaces. In particular, we use electrochemistry to explore the extent to which a macroscopic, hydroxyl-coated surface held at a slightly negative potential affects the folding thermodynamics of surface-attached protein variants with different composition of charged amino acids. Doing so, we find that attachment to the surface generally leads to a net stabilization, presumably due to excluded volume effects that reduce the entropy of the unfolded state. The magnitude of this stabilization, however, is strongly correlated with the charged-residue content of the protein. In particular, we find statistically significant correlations with both the net charge of the protein, with greater negative charge leading to less stabilization by the surface, and with the number of arginines, with more arginines leading to greater stabilization. Such findings refine our understanding of protein-surface interactions, providing in turn a guiding rationale to achieve the functional deposition of proteins on artificial surfaces for implementation in, for example, protein-based biotechnologies.

Distinct and Nonadditive Effects of Urea and Guanidinium Chloride on Peptide Solvation

Using enhanced-sampling replica exchange fully atomistic molecular dynamics simulations, we show that, individually, urea and guanidinium chloride (GdmCl) denature the Trpcage protein, but remarkably, the helical segment ¹NLYIQWL⁷ of the protein is stabilized in mixed denaturant solutions. GdmCl induces protein denaturation via a combination of direct and indirect effects involving dehydration of the protein and destabilization of stabilizing salt bridges. In contrast, urea denatures the protein through favorable protein-urea preferential interactions, with peptide-specific indirect effects of urea on the water structure around the protein. In the case of the helical segment of Trpcage, urea "oversolvates" the peptide backbone by reorganizing water molecules from the peptide side chains to the peptide backbone. An intricate nonadditive thermodynamic balance between GdmCl-induced dehydration of the peptide and the urea-induced changes in solvation structure triggers partial counteraction to urea denaturation and stabilization of the helix.

Surface Attachment Enhances the Thermodynamic Stability of Protein L

Despite the importance of protein-surface interactions in both biology and biotechnology, our understanding of their origins is limited due to a paucity of experimental studies of the thermodynamics behind such interactions. In response, we have characterized the extent to which interaction with a chemically well-defined macroscopic surface alters the stability of protein L. To do so, we site-specifically attached a redox-reporter-modified protein variant to a hydroxy-terminated monolayer on a gold surface and then used electrochemistry to monitor its guanidine denaturation and determine its folding free energy. Comparison with the free energy seen in solution indicates that interaction with this surface stabilizes the protein by 6 kJ mol^-1 , a value that is in good agreement with theoretical estimates of the entropic consequences of surface-induced excluded volume effects, thus suggesting that chemically specific interactions with this surface (e.g., electrostatic interactions) are limited in magnitude.

Molecular Dynamics Simulation of Protein Biosurfactants

Surfaces and interfaces are ubiquitous in nature and are involved in many biological processes. Due to this, natural organisms have evolved a number of methods to control interfacial and surface properties. Many of these methods involve the use of specialised protein biosurfactants, which due to the competing demands of high surface activity, biocompatibility, and low solution aggregation may take structures that differ from the traditional head–tail structure of small molecule surfactants. As well as their biological functions, these proteins have also attracted interest for industrial applications, in areas including food technology, surface modification, and drug delivery. To understand the biological functions and technological applications of protein biosurfactants, it is necessary to have a molecular level description of their behaviour, in particular at surfaces and interfaces, for which molecular simulation is well suited to investigate. In this review, we will give an overview of simulation studies of a number of examples of protein biosurfactants (hydrophobins, surfactin, and ranaspumin). We will also outline some of the key challenges and future directions for molecular simulation in the investigation of protein biosurfactants and how this can help guide future developments.

Quantitative measurements of protein-surface interaction thermodynamics

Whereas proteins generally remain stable upon interaction with biological surfaces, they frequently unfold on and adhere to artificial surfaces. Understanding the physicochemical origins of this discrepancy would facilitate development of protein-based sensors and other technologies that require surfaces that do not compromise protein structure and function. To date, however, only a small number of such artificial surfaces have been reported, and the physics of why these surfaces support functional biomolecules while others do not has not been established. Thus motivated, we have developed an electrochemical approach to determining the folding free energy of proteins site-specifically attached to chemically well-defined, macroscopic surfaces. Comparison with the folding free energies seen in bulk solution then provides a quantitative measure of the extent to which surface interactions alter protein stability. As proof-of-principle, we have characterized the FynSH3 domain site-specifically attached to a hydroxyl-coated surface. Upon guanidinium chloride denaturation, the protein unfolds in a reversible, two-state manner with a free energy within 2 kJ/mol of the value seen in bulk solution. Assuming that excluded volume effects stabilize surface-attached proteins, this observation suggests there are countervening destabilizing interactions with the surface that, under these conditions, are similar in magnitude. Our technique constitutes an unprecedented experimental tool with which to answer long-standing questions regarding the molecular-scale origins of protein-surface interactions and to facilitate rational optimization of surface biocompatibility.

Facet-Specific Adsorption of Tripeptides at Aqueous Au Interfaces: Open Questions in Reconciling Experiment and Simulation

The adsorption of three homo-tripeptides, HHH, YYY, and SSS, at the aqueous Au interface is investigated, using molecular dynamics simulations. We find that consideration of surface facet effects, relevant to experimental conditions, opens up new questions regarding interpretations of current experimental findings. Our well-tempered metadynamics simulations predict the rank ordering of the tripeptide binding affinities at aqueous Au(111) to be YYY > HHH > SSS. This ranking differs with that obtained from existing experimental data which used surface-immobilized Au nanoparticles as the target substrate. The influence of Au facet on these experimental findings is then considered, via our binding strength predictions of the relevant amino acids at aqueous Au(111) and Au(100)(1 × 1). The Au(111) interface supports an amino acid ranking of Tyr > HisA ≃ HisH > Ser, matching that of the tripeptides on Au(111), while the ranking on Au(100) is HisA > Ser ≃ Tyr ≃ HisH, with only HisA showing non-negligible binding. The substantial reduction in Tyr amino acid affinity for Au(100) vs Au(111) offers one possible explanation for the experimentally observed weaker adsorption of YYY on the nanoparticle-immobilized substrate compared with HHH. In a separate set of simulations, we predict the structures of the adsorbed tripeptides at the two aqueous Au facets, revealing facet-dependent differences in the adsorbed conformations. Our findings suggest that Au facet effects, where relevant, may influence the adsorption structures and energetics of biomolecules, highlighting the possible influence of the structural model used to interpret experimental binding data.

Simulations of disordered proteins and systems with conformational heterogeneity

Intrinsically disordered proteins (IDPs) and protein regions can facilitate a wide variety of complex physiological processes such as binding, signaling, and formation of membraneless organelles. They can however also play pathological roles by aggregating into cytotoxic oligomers and fibrils. Characterizing the structure and function of disordered proteins is an onerous task, primarily because these proteins adopt transient structures, which are difficult to capture in experiments. Simulations have emerged as a powerful tool for interpreting and augmenting experimental measurements of IDPs. In this review we focus on computer simulations of disordered protein structures, functions, assemblies, and emerging questions that, taken together, give an overview of the field as it exists today.

Protein States as Symmetry Transitions in the Correlation Matrices

Over the last few years, there has been significant progress in the knowledge on protein folding. However, some aspects of protein folding still need further attention. One of these is the exact relationship between the folded and unfolded states and the differences between them. Whereas the folded state is well known, at least from a structural point of view (just think of the thousands of structures in online databases), the unfolded state is more elusive. Also, these are dynamic states of matter, and this aspect cannot be overlooked. Molecular dynamics-derived correlation matrices are an invaluable source of information on the protein dynamics. Here, bulk eigenvalue spectra of the correlation matrices obtained from the Trp-cage dynamics in the folded and unfolded states have been analyzed. The associated modes represent localized vibrations and are significantly affected by the fine details of the structure and interactions. Therefore, these bulk modes can be used as probes of the protein local dynamics in different states. The results of these analyses show that the correlation matrices describing the folded and unfolded dynamics belong to different symmetry classes. This finding provides new support to the phase-transition models of protein folding.

Strong Electrostatic Interactions Lead to Entropically Favorable Binding of Peptides to Charged Surfaces

Thermodynamic analyses can provide key insights into the origins of protein self-assembly on surfaces, protein function, and protein stability. However, obtaining quantitative measurements of thermodynamic observables from unbiased classical simulations of peptide or protein adsorption is challenging because of sampling limitations brought on by strong biomolecule/surface binding forces as well as time scale limitations. We used the parallel tempering metadynamics in the well-tempered ensemble (PTMetaD-WTE) enhanced sampling method to study the adsorption behavior and thermodynamics of several explicitly solvated model peptide adsorption systems, providing new molecular-level insight into the biomolecule adsorption process. Specifically studied were peptides LKα14 and LKβ15 and trpcage miniprotein adsorbing onto a charged, hydrophilic self-assembled monolayer surface functionalized with a carboxylic acid/carboxylate headgroup and a neutral, hydrophobic methyl-terminated self-assembled monolayer surface. Binding free energies were calculated as a function of temperature for each system and decomposed into their respective energetic and entropic contributions. We investigated how specific interfacial features such as peptide/surface electrostatic interactions and surface-bound ion content affect the thermodynamic landscape of adsorption and lead to differences in surface-bound conformations of the peptides. Results show that upon adsorption to the charged surface, configurational entropy gains of the released solvent molecules dominate the configurational entropy losses of the bound peptide. This behavior leads to an apparent increase in overall system entropy upon binding and therefore to the surprising and seemingly nonphysical result of an apparent increased binding free energy at elevated temperatures. Opposite effects and conclusions are found for the neutral surface. Additional simulations demonstrate that by adjusting the ionic strength of the solution, results that show the expected physical behavior, i.e., peptide binding strength that decreases with increasing temperature or is independent of temperature altogether, can be recovered on the charged surface. On the basis of this analysis, an overall free energy for the entire thermodynamic cycle for peptide adsorption on charged surfaces is constructed and validated with independent simulations.

Conformations of Myoglobin-Derived Peptides at the Air-Water Interface

The conformational change exhibited by proteins at liquid interfaces, such as the air-water and oil-water interfaces, has long been of interest both for understanding protein structure outside of native environments and for applications in areas including food technology and pharmaceuticals. Using molecular simulation, this article studies the conformations of two peptides derived from myoglobin, for which the emulsification behavior has been studied. Both peptides were found to readily adsorb onto the air-water interface, with one of these (experimentally, the more effective stabilizer) adopting a flat, extended conformation and the other peptide remaining close to its solution conformation.

Surface force measurements and simulations of mussel-derived peptide adhesives on wet organic surfaces

Translating sticky biological molecules-such as mussel foot proteins (MFPs)-into synthetic, cost-effective underwater adhesives with adjustable nano- and macroscale characteristics requires an intimate understanding of the glue's molecular interactions. To help facilitate the next generation of aqueous adhesives, we performed a combination of surface forces apparatus (SFA) measurements and replica-exchange molecular dynamics (REMD) simulations on a synthetic, easy to prepare, Dopa-containing peptide (MFP-3s peptide), which adheres to organic surfaces just as effectively as its wild-type protein analog. Experiments and simulations both show significant differences in peptide adsorption on CH3-terminated (hydrophobic) and OH-terminated (hydrophilic) self-assembled monolayers (SAMs), where adsorption is strongest on hydrophobic SAMs because of orientationally specific interactions with Dopa. Additional umbrella-sampling simulations yield free-energy profiles that quantitatively agree with SFA measurements and are used to extract the adhesive properties of individual amino acids within the context of MFP-3s peptide adhesion, revealing a delicate balance between van der Waals, hydrophobic, and electrostatic forces.

Correlation Analysis of Trp-Cage Dynamics in Folded and Unfolded States

A fundamental and still debated problem is how folded structures of proteins are related to their unfolded state. Besides the classical view, in which a large number of conformations characterize the unfolded state while the folded one is dominated by a single structure, recently a reassessment of the denatured state has been suggested. A growing amount of evidence indicates that not only the folded but also the unfolded state is at least partially organized. Here, we try to answer the question of how different protein dynamics is in folded and unfolded states by performing all-atom molecular dynamics simulations on the model protein Trp-cage. Random matrix theory inspired analysis of the correlation matrices has been carried out. The spectra of these correlation matrices show that the low rank modes of Trp-cage dynamics are outside of the limit expected for a random system both in folded and in unfolded conditions. These findings shed light on the nature of the unfolded state of the proteins, suggesting that it is much less random than previously thought.