Nature of the charged-group effect on the stability of the C-peptide helix

Toward a Monte Carlo simulation of protein systems in amino-acid sequence space

In this article, we present our strategy for studying amino-acid sequence dependences on protein structures. For this purpose, performing Metropolis Monte Carlo simulations in the amino-acid sequence space is necessary. We want to use a coarse-grained protein model with an accurate potential energy function. We introduce a method for optimizing potential-energy parameters based on the native protein structure database, Protein Data Bank.

A functional helix shuffled variant of the B domain of Staphylococcus aureus

The B domain of protein A is a biotechnologically important three-helix bundle protein. It binds the Fc fragment of antibodies with helix 1/2 and the Fab region with helix 2/3. Here we designed a helix shuffled variant by changing the connectivity of the helices, in order to redesign the helix bundle, yielding altered helix-loop-helix properties. The new loops that generate the new connectivity were created in several protein libraries, and Fc binding variants were selected for a detailed biochemical characterization. We were able to create variants with Fc binding affinity at the same level as the wild type B but with significantly reduced thermal stability. The NMR structure proved that the overall three-dimensional structure was maintained not only in the helix shuffled variant but also points to some potential local differences to wild-type B, which could be the reason for the reduced thermal stability. Therefore, protein A is an example of an optimized structure being more important for stability than for function. Using the helix shuffled variant as a ligand on an affinity column facilitates a robust and straightforward purification of antibodies, but allows for a milder elution at less extreme pH. Therefore, the helix shuffled variant is a suitable ligand to purify more pH-sensitive antibodies.

Statistical analysis of the unique characteristics of secondary structures in proteins

Protein folding is a complex process influenced by the primary sequence of amino acids. Early studies focused on understanding whether the specificity or the conservation of properties of amino acids was crucial for folding into secondary structures such as α-helices, β-sheets, turns, and coils. However, with the advent of artificial intelligence (AI) and machine learning (ML), the emphasis has shifted towards the precise nature and occurrence of specific amino acids. In our study, we analyzed a large set of proteins from diverse organisms to identify unique features of secondary structures, particularly in terms of the distribution of polar, non-polar, and charged amino acid residues. We found that α-helices tend to have a higher proportion of charged and non-polar groups compared to other secondary structures and that the presence of oppositely charged amino acid residues in helices stabilizes them, facilitating the formation of longer helices. These characteristics are distinct to α-helices. This study offers valuable insights for researchers in the field of protein design, enabling the de-novo creation of short helical peptides for a range of applications. We have also developed a web server for extensive analysis of proteins from different databases. The web server is housed at https://proseqanalyser.iitgn.ac.in/.

Short Oligourea Foldamers as N- or C-Caps for Promoting α-Helix Formation in Water

While foldamers have been extensively studied as protein mimics and especially as α-helix mimics, their use as capping motif to enhance α-helix propensity remains comparatively much limited. In this study, we leverage the structural similarities between urea-based helical foldamers and α-helix to investigate the efficacy of oligoureas as N- or C-caps for reinforcing α-helical structures in water. Short oligoureas, comprising 3 to 4 residues, were strategically introduced at the N- or C-terminus of two peptide sequences (S-peptide and an Ala-rich model sequence). The impact of these foldamer insertions on peptide conformation was examined using electronic circular dichroism (ECD) and solution NMR. This research identifies specific foldamer sequences capable of promoting α-helicity when incorporated at either terminus of the peptides. Not only does this work broaden the application scope of foldamers, but it also provides valuable insights into novel strategies for modulating peptide conformation in aqueous environments. The findings presented in this study may have implications for peptide design and the development of bioactive foldamer-based peptide mimics.

Helix-Coil Transition at a Glycine Following a Nascent α-Helix: A Synergetic Guidance Mechanism for Helix Growth

A detailed understanding of forces guiding the rapid folding of a polypeptide from an apparently random coil state to an ordered α-helical structure following the rate-limiting preorganization of the initial three residue backbones into helical conformation is imperative to comprehending and regulating protein folding and for the rational design of biological mimetics. However, several details of this process are still unknown. First, although the helix-coil transition was proposed to originate at the residue level (J. Chem. Phys. 1959, 31, 526-535; J. Chem. Phys. 1961, 34, 1963-1974), all helix-folding studies have only established it between time-averaged bulk states of a long-lived helix and several transiently populated random coils, along the whole helix model sequence. Second, the predominant thermodynamic forces driving either this two-state transition or the faster helix growth following helix nucleation are still unclear. Third, the conformational space of the random coil state is not well-defined unlike its corresponding α-helix. Here we investigate the restrictions placed on the conformational space of a Gly residue backbone, as a result of it immediately succeeding a nascent α-helical turn. Analyses of the temperature-dependent 1D-, 2D-NMR, FT-IR, and CD spectra and GROMACS MD simulation trajectory of a Gly residue backbone following a model α-helical turn, which is artificially rigidified by a covalent hydrogen bond surrogate, reveal that: (i) the α-helical turn guides the ϕ torsion of the Gly exclusively into either a predominantly populated entropically favored α-helical (α-ϕ) state or a scarcely populated random coil (RC-ϕ) state; (ii) the α-ϕ state of Gly in turn favors the stability of the preceding α-helical turn, while the RC-ϕ state disrupts it, revealing an entropy-driven synergetic guidance for helix growth in the residue following helix nucleation. The applicability of a current synergetic guidance mechanism to explain rapid helix growth in folded and unfolded states of proteins and helical peptides is discussed.

Further Optimization and Validation of the Classical Drude Polarizable Protein Force Field

The CHARMM Drude-2013 polarizable force field (FF) was developed to include the explicit treatment of induced electronic polarizability, resulting in a more accurate description of the electrostatic interactions in molecular dynamics (MD) simulations. While the Drude-2013 protein FF has shown success in improving the folding properties of α-helical peptides and to reproduce experimental observables in simulations up to 1 μs, some limitations were noted regarding the stability of β-sheet structures in simulations longer than 100 ns as well as larger deviations from crystal structures in simulations of a number of proteins compared to the additive CHARMM36 protein FF. The origin of the instability has been identified and appears to be primarily due to overestimated atomic polarizabilities and induced dipole-dipole interactions on the Cβ, Cγ, and Cδ side chain atoms. To resolve this and other issues, a number of aspects of the model were revisited, resulting in Drude-2019 protein FF. Backbone parameters were optimized targeting the conformational properties of the (Ala)5 peptide in solution along with gas phase properties of the alanine dipeptide. Dipeptides that contain N-acetylated and N'-methylamidated termini, excluding Gly, Pro, and Ala, were used as models to optimize the atomic polarizabilities and Thole screening factors on selected Cβ, Cγ, and Cδ carbons by targeting quantum mechanical (QM) dipole moments and molecular polarizabilities. In addition, to obtain better conformational properties, side chain χ1 and χ2 dihedral parameters were optimized targeting QM data for the respective side chain dipeptide conformations as well as Protein Data Bank survey data based on the χ1, χ2 sampling from Hamiltonian replica-exchange MD simulations of (Ala)4-X-(Ala)4 in solution, where X is the amino acid of interest. Further improvements include optimizing nonbonded interactions between charged residues to reproduce QM interaction energies of the charged-protein model compounds and experimental osmotic pressures. Validation of the optimized Drude protein FF includes MD simulations of a collection of peptides and proteins including β-sheet structures, as well as transmembrane ion channels. Results showed that the updated Drude-2019 protein FF yields smaller overall root-mean-square differences of proteins as compared to the additive CHARMM36m and Drude-2013 FFs as well as similar or improved agreement with experimental NMR properties, allowing for long time scale simulation studies of proteins and more complex biomolecular systems in conjunction with the remainder of the Drude polarizable FF.

MALDI-TOF mass spectrometry-based quantification of C-peptide in diabetes patients

BACKGROUND

Serum C-peptide concentrations reflect insulin secretion and beta cell function and can be used to diagnose and distinguish type-1 and type-2 diabetes. C-peptide is a more accurate indicator of insulin status than direct insulin measurement for monitoring patients with diabetes. However, the current methods available for C-peptide quantification exhibit poor reproducibility, are costly, and require highly trained laboratory personnel. Here, we have developed and evaluated a matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS)-based assay to standardize C-peptide measurements, providing highly accurate and comparable results across testing systems and laboratories.

METHODS

C-peptide from human serum was enriched using antibody-conjugated magnetic beads. The eluted isolates were further modified with 6-aminoquinolyl-N-hydroxysuccinimidyl carbamate (AQC) to enhance the ionization of naturally acidic C-peptide. After desalting with ZipTips, the samples were subjected to MALDI-TOF MS analysis. Recombinant human C-peptide was used to develop the assay, and a heavy isotope labeled human C-peptide was used as an internal standard for quantification.

RESULTS

The MALDI-TOF MS method was validated in accordance with the restrictions of the device, with a limit of quantitation of 25 pmol/L. A correlation between the MAL-DI-TOF MS assay and a reference method was conducted using patient samples. The resulting regression revealed good agreement.

CONCLUSIONS

A simple, high-throughput, cost effective and quantitative MALDI-TOF MS C-peptide assay has been successfully developed and validated in clinical serum samples.

Thermodynamics of Helix formation in small peptides of varying lengthin vacuo, implicit solvent and explicit solvent: Comparison between AMBER force fields

Helix formation is of great significance in protein folding. The helix-forming tendencies of amino acids are accumulated along the sequence to determine the helix-forming tendency of peptides. Computer simulation can be used to model this process in atomic details and give structural insights. In the current work, we employ equilibrate-state free energy simulation to systematically study the folding/unfolding thermodynamics of a series of mutated peptides. Two AMBER force fields including AMBER99SB and AMBER14SB are compared. The new 14SB force field uses refitted torsion parameters compared with 99SB and they share the same atomic charge scheme. We find that in vacuo the helix formation is mutation dependent, which reflects the different helix propensities of different amino acids. In general, there are helix formers, helix indifferent groups and helix breakers. The helical structure becomes more favored when the length of the sequence becomes longer, which arises from the formation of additional backbone hydrogen bonds in the lengthened sequence. Therefore, the helix indifferent groups and helix breakers will become helix formers in long sequences. Also, protonation-dependent helix formation is observed for ionizable groups. In 14SB, the helical structures are more stable than in 99SB and differences can be observed in their grouping schemes, especially in the helix indifferent group. In solvents, all mutations are helix indifferent due to protein–solvent interactions. The decrease in the number of backbone hydrogen bonds is the same with the increase in the number of protein–water hydrogen bonds. The 14SB in explicit solvent is able to capture the free energy minima in the helical state while 14SB in implicit solvent, 99SB in explicit solvent and 99SB in implicit solvent cannot. The helix propensities calculated under 14SB agree with the corresponding experimental values, while the 99SB results obviously deviate from the references. Hence, implicit solvent models are unable to correctly describe the thermodynamics even for the simple helix formation in isolated peptides. Well-developed force fields and explicit solvents are needed to correctly describe the protein dynamics. Aside from the free energy, differences in conformational ensemble under different force fields in different solvent models are observed. The numbers of hydrogen bonds formed under different force fields agree and they are mostly determined by the solvent model.

One Peptide Reveals the Two Faces of α-Helix Unfolding-Folding Dynamics

The understanding of fast folding dynamics of single α-helices comes mostly from studies on rationally designed peptides displaying sequences with high helical propensity. The folding/unfolding dynamics and energetics of α-helix conformations in naturally occurring peptides remains largely unexplored. Here we report the study of a protein fragment analogue of the C-peptide from bovine pancreatic ribonuclease-A, RN80, a 13-amino acid residue peptide that adopts a highly populated helical conformation in aqueous solution. ¹H NMR and CD structural studies of RN80 showed that α-helix formation displays a pH-dependent bell-shaped curve, with a maximum near pH 5, and a large decrease in helical content in alkaline pH. The main forces stabilizing this short α-helix were identified as a salt bridge formed between Glu-2 and Arg-10 and the cation-π interaction involving Tyr-8 and His-12. Thus, deprotonation of Glu-2 or protonation of His-12 are essential for the RN80 α-helix stability. In the present study, RN80 folding and unfolding were triggered by laser-induced pH jumps and detected by time-resolved photoacoustic calorimetry (PAC). The photoacid proton release, amino acid residue protonation, and unfolding/folding events occur at different time scales and were clearly distinguished using time-resolved PAC. The partial unfolding of the RN80 α-helix, due to protonation of Glu-2 and consequent breaking of the stabilizing salt bridge between Glu-2 and Arg-10, is characterized by a concentration-independent volume expansion in the sub-microsecond time range (0.8 mL mol^-1, 369 ns). This small volume expansion reports the cost of peptide backbone rehydration upon disruption of a solvent-exposed salt bridge, as well as backbone intrinsic expansion. On the other hand, RN80 α-helix folding triggered by His-12 protonation and subsequent formation of a cation-π interaction leads to a microsecond volume contraction (-6.0 mL mol^-1, ∼1.7 μs). The essential role of two discrete side chain interactions, a salt bridge, and in particular a single cation-π interaction in the folding dynamics of a naturally occurring α-helix peptide is uniquely revealed by these data.

Exendin-4-encapsulated dissolving microneedle arrays for efficient treatment of type 2 diabetes

Dissolving microneedles (DMNs) are microscopic needles capable of delivering encapsulated compounds and releasing them into the skin in a minimally invasive manner. Most studies indicate that encapsulating therapeutics in DMNs is an efficacious approach; however, the importance of evaluating the activity of encapsulated compounds, during the fabrication process, has not been examined in detail. Conducting an analysis of thermal, chemical, and physical stress factors, including temperature, pH, and the interaction of the polymer and therapeutics mixture during preparation, is essential for retaining the activity of encapsulated therapeutics during and after fabrication. Here, we optimised the thermal, chemical, and physical parameters for the fabrication of exendin-4 (Ex-4)-encapsulated DMNs (Ex-4 DMNs). Ex-4, a peptide agonist of glucagon-like peptide (GLP) receptor, is used for glycaemic control in patients with type 2 diabetes. Our findings indicate that optimising the parameters involved in DMN fabrication retained the activity of Ex-4 by up to 98.3 ± 1.5%. Ex-4 DMNs reduced the blood-glucose level in diabetic mice with efficiency similar to that of a subcutaneous injection. We believe that this study paves way for the commercialisation of an efficient and minimally invasive treatment for patients with type 2 diabetes.

Comparative Molecular Dynamics Analysis of RNase-S Complex Formation

Limited proteolysis of RNase-A yields a short N-terminal S-peptide segment and the larger S-protein. Binding of S-peptide to S-protein results in the formation of an enzymatically active RNase-S protein. S-peptide undergoes a transition from intrinsic disorder to an ordered helical state upon association with S-protein to form RNase-S and is an excellent model system to study coupled folding and binding. To better understand the dynamics of the RNases-S complex and its isolated partners, comparative molecular dynamics simulations have been performed. In agreement with experiment, we find significant conformational fluctuations of the isolated S-peptide compatible with a disordered regime and only little residual helical structure. In the RNase-S complex, the N-terminal helix of S-peptide unfolds and refolds repeatedly on the microsecond timescale, indicating that the α-helical structure is only part of the equilibrium regime for these residues whereas the C-terminal residues are confined to the helical conformation that is found in the x-ray structure. This is also in line with systematic, in silico Alanine scanning free-energy simulations, which indicate that the major contribution to complex stability emerges from the C-terminal helical turn, consisting of residues 8-13 in S-peptide whereas the N-terminal S-peptide residues 1-7 make only minor contributions. Comparative simulations of S-protein in the presence and absence of S-peptide reveal that the isolated S-protein is significantly more flexible than in the complex, and undergoes a global pincerlike conformational change that narrows the S-peptide binding cleft. The narrowed binding cleft adds a barrier for complex formation likely influencing the binding kinetics. This conformational change is reversed by S-peptide association, which also stabilizes conformational fluctuations in S-protein. Such global motions associated with binding are also likely to play a role for other coupled peptide folding and binding processes at peptide binding regions on protein surfaces.

Design of Stable α-Helical Peptides and Thermostable Proteins in Biotechnology and Biomedicine

α-Helices are the most frequently occurring elements of the secondary structure in water-soluble globular proteins. Their increased conformational stability is among the main reasons for the high thermal stability of proteins in thermophilic bacteria. In addition, α-helices are often involved in protein interactions with other proteins, nucleic acids, and the lipids of cell membranes. That is why the highly stable α-helical peptides used as highly active and specific inhibitors of protein-protein and other interactions have recently found more applications in medicine. Several different approaches have been developed in recent years to improve the conformational stability of α-helical peptides and thermostable proteins, which will be discussed in this review. We also discuss the methods for improving the permeability of peptides and proteins across cellular membranes and their resistance to intracellular protease activity. Special attention is given to the SEQOPT method (http://mml.spbstu.ru/services/seqopt/), which is used to design conformationally stable short α-helices.

Grid-based backbone correction to the ff12SB protein force field for implicit-solvent simulations

Force fields, such as Amber's ff12SB, can be fairly accurate models of the physical forces in proteins and other biomolecules. When coupled with accurate solvation models, force fields are able to bring insight into the conformational preferences, transitions, pathways, and free energies for these biomolecules. When computational speed/cost matters, implicit solvent is often used but at the cost of accuracy. We present an empirical grid-like correction term, in the spirit of cMAPs, to the combination of the ff12SB protein force field and the GBneck2 implicit-solvent model. Ff12SB-cMAP is parametrized on experimental helicity data. We provide validation on a set of peptides and proteins. Ff12SB-cMAP successfully improves the secondary structure biases observed in ff12SB + Gbneck2. Ff12SB-cMAP can be downloaded ( https://github.com/laufercenter/Amap.git ) and used within the Amber package. It can improve the agreement of force fields + implicit solvent with experiments.

The complex interplay between ligand binding and conformational structure of the folate binding protein (folate receptor): Biological perspectives

This review analyzes how interplay between folate binding and changes in folate binding protein (FBP) conformation/self-association affects the biological function of FBP. Concentration-dependent, reversible self-association of hydrophobic apo-FBP at pI=7.4 is associated with decreased affinity for folate, probably due to shielding of binding sites between interacting hydrophobic patches. Titration with folate removes apo-monomers, favoring dissociation of self-associated apo-FBP into apo-monomers. Folate anchors to FBP through a network of hydrogen bonds and hydrophobic interactions, and the binding induces a conformational change with formation of hydrophilic and stable holo-FBP. Holo-FBP exhibits a ligand-mediated concentration-dependent self-association into multimers of great thermal and chemical stability due to strong intermolecular forces. Both ligand and FBP are thus protected against biological/physicochemical decomposition. In biological fluids with low FBP concentrations, e.g., saliva, semen and plasma, hydrophobic apo-monomers and hydrophilic holo-monomers associate into stable asymmetrical complexes with aberrant binding kinetics unless detergents, e.g., cholesterol or phospholipids are present.

High-resolution crystal structures of protein helices reconciled with three-centered hydrogen bonds and multipole electrostatics

Theoretical and experimental evidence for non-linear hydrogen bonds in protein helices is ubiquitous. In particular, amide three-centered hydrogen bonds are common features of helices in high-resolution crystal structures of proteins. These high-resolution structures (1.0 to 1.5 Å nominal crystallographic resolution) position backbone atoms without significant bias from modeling constraints and identify Φ = -62°, ψ = -43 as the consensus backbone torsional angles of protein helices. These torsional angles preserve the atomic positions of α-β carbons of the classic Pauling α-helix while allowing the amide carbonyls to form bifurcated hydrogen bonds as first suggested by Némethy et al. in 1967. Molecular dynamics simulations of a capped 12-residue oligoalanine in water with AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications), a second-generation force field that includes multipole electrostatics and polarizability, reproduces the experimentally observed high-resolution helical conformation and correctly reorients the amide-bond carbonyls into bifurcated hydrogen bonds. This simple modification of backbone torsional angles reconciles experimental and theoretical views to provide a unified view of amide three-centered hydrogen bonds as crucial components of protein helices. The reason why they have been overlooked by structural biologists depends on the small crankshaft-like changes in orientation of the amide bond that allows maintenance of the overall helical parameters (helix pitch (p) and residues per turn (n)). The Pauling 3.6₁₃ α-helix fits the high-resolution experimental data with the minor exception of the amide-carbonyl electron density, but the previously associated backbone torsional angles (Φ, Ψ) needed slight modification to be reconciled with three-atom centered H-bonds and multipole electrostatics. Thus, a new standard helix, the 3.6_13/10-, Némethy- or N-helix, is proposed. Due to the use of constraints from monopole force fields and assumed secondary structures used in low-resolution refinement of electron density of proteins, such structures in the PDB often show linear hydrogen bonding.

Current status of protein force fields for molecular dynamics simulations

The current status of classical force fields for proteins is reviewed. These include additive force fields as well as the latest developments in the Drude and AMOEBA polarizable force fields. Parametrization strategies developed specifically for the Drude force field are described and compared with the additive CHARMM36 force field. Results from molecular simulations of proteins and small peptides are summarized to illustrate the performance of the Drude and AMOEBA force fields.

Effects of charge states, charge sites and side chain interactions on conformational preferences of a series of model peptide ions

The factors affecting conformational preference of gas phase peptide ions are investigated by IM-MS and molecular dynamics simulation.

Effect of external pulling forces on the length distribution of peptides

BACKGROUND

The distribution of the length of a polypeptide, or that of the distance between any two of its atoms, is an important property as it can be analytically or numerically estimated for a number of polymer models. Importantly, it is directly measurable through a number of different experimental techniques. Length distributions can be straightforwardly assessed from molecular dynamics simulation; however, true convergence through full accurate coverage of the length range is difficult to achieve.

METHODS

The application of external constant force combined with the weighted-histogram analysis method (WHAM) is used to enhance sampling of unlikely 'long' or 'short' conformations and obtain the potential of mean force, while also collecting dynamic properties of the chain under variable tension.

RESULTS

We demonstrate the utility of constant force to enhance the sampling efficiency and obtain experimentally measurable quantities on a series of short peptides, including charge-rich sequences that are known to be highly helical but whose properties are distinct from those of helical peptides undergoing helix-coil transitions.

CONCLUSIONS

Force-enhanced sampling enhances the range and accuracy of the length-based potential of mean force of the peptide, in particular those sequences that contain increased numbers of charged residues.

GENERAL SIGNIFICANCE

This approach allows users to simultaneously probe the force-dependent behaviour of peptides directly, enhance the range and accuracy of the length-based PMF of the peptide and also test the convergence of simulations by comparing the overlap of PMF profiles from different constant forces. This article is part of a special issue entitled Recent developments of molecular dynamics.

Translating Koch's postulates to identify matrix metalloproteinase roles in postmyocardial infarction remodeling: cardiac metalloproteinase actions (CarMA) postulates

The first matrix metalloproteinase (MMP) was described in 1962; and since the 1990s, cardiovascular research has focused on understanding how MMPs regulate many aspects of cardiovascular pathology from atherosclerosis formation to myocardial infarction and stroke. Although much information has been gleaned by these past reports, to a large degree MMP cardiovascular biology remains observational, with few studies homing in on cause and effect relationships. Koch's postulates were first developed in the 19th century as a way to establish microorganism function and were modified in the 20th century to include methods to establish molecular causality. In this review, we outline the concept for establishing a similar approach to determine causality in terms of MMP functions. We use left ventricular remodeling postmyocardial infarction as an example, but this approach will have broad applicability across both the cardiovascular and the MMP fields.

How force unfolding differs from chemical denaturation

Single-molecule force spectroscopies are remarkable tools for studying protein folding and unfolding, but force unfolding explores protein configurations that are potentially very different from the ones traditionally explored in chemical or thermal denaturation. Understanding these differences is crucial because such configurations serve as starting points of folding studies, and thus can affect both the folding mechanism and the kinetics. Here we provide a detailed comparison of both chemically induced and force-induced unfolded state ensembles of ubiquitin based on extensive, all-atom simulations of the protein either extended by force or denatured by urea. As expected, the respective unfolded states are very different on a macromolecular scale, being fully extended under force with no contacts and partially extended in urea with many nonnative contacts. The amount of residual secondary structure also differs: A significant population of α-helices is found in chemically denatured configurations but such helices are absent under force, except at the lowest applied force of 30 pN where short helices form transiently. We see that typical-size helices are unstable above this force, and β-sheets cannot form. More surprisingly, we observe striking differences in the backbone dihedral angle distributions for the protein unfolded under force and the one unfolded by denaturant. A simple model based on the dialanine peptide is shown to not only provide an explanation for these striking differences but also illustrates how the force dependence of the protein dihedral angle distributions give rise to the worm-like chain behavior of the chain upon force.

Constant pH molecular dynamics of proteins in explicit solvent with proton tautomerism

pH is a ubiquitous regulator of biological activity, including protein-folding, protein-protein interactions, and enzymatic activity. Existing constant pH molecular dynamics (CPHMD) models that were developed to address questions related to the pH-dependent properties of proteins are largely based on implicit solvent models. However, implicit solvent models are known to underestimate the desolvation energy of buried charged residues, increasing the error associated with predictions that involve internal ionizable residue that are important in processes like hydrogen transport and electron transfer. Furthermore, discrete water and ions cannot be modeled in implicit solvent, which are important in systems like membrane proteins and ion channels. We report on an explicit solvent constant pH molecular dynamics framework based on multi-site λ-dynamics (CPHMD(MSλD)). In the CPHMD(MSλD) framework, we performed seamless alchemical transitions between protonation and tautomeric states using multi-site λ-dynamics, and designed novel biasing potentials to ensure that the physical end-states are predominantly sampled. We show that explicit solvent CPHMD(MSλD) simulations model realistic pH-dependent properties of proteins such as the Hen-Egg White Lysozyme (HEWL), binding domain of 2-oxoglutarate dehydrogenase (BBL) and N-terminal domain of ribosomal protein L9 (NTL9), and the pKa predictions are in excellent agreement with experimental values, with a RMSE ranging from 0.72 to 0.84 pKa units. With the recent development of the explicit solvent CPHMD(MSλD) framework for nucleic acids, accurate modeling of pH-dependent properties of both major class of biomolecules-proteins and nucleic acids is now possible.

The interrelationship between ligand binding and thermal unfolding of the folate binding protein. The role of self-association and pH

The present study utilized a combination of DLS (dynamic light scattering) and DSC (differential scanning calorimetry) to address thermostability of high-affinity folate binding protein (FBP), a transport protein and cellular receptor for the vitamin folate. At pH7.4 (pI=7-8) ligand binding increased concentration-dependent self-association of FBP into stable multimers of holo-FBP. DSC of 3.3μM holo-FBP showed Tm (76°C) and molar enthalpy (146kcalM(-1)) values increasing to 78°C and 163kcalM(-1) at 10μM holo-FBP, while those of apo-FBP were 55°C and 105kcalM(-1). Besides ligand binding, intermolecular forces involved in concentration-dependent multimerization thus contribute to the thermostability of holo-FBP. Hence, thermal unfolding and dissociation of holo-FBP multimers occur simultaneously consistent with a gradual decrease from octameric to monomeric holo-FBP (10μM) in DLS after a step-wise rise in temperature to 78°C≈Tm. Stable holo-FBP multimers may protect naturally occurring labile folates against decomposition or bacterial utilization. DSC established an interrelationship between diminished folate binding at pH5, especially in NaCl-free buffers, and low thermostability. Positively charged apo-FBP was almost completely unfolded and aggregated at pH5 (Tm 38°C) and holo-FBP, albeit more thermostable, was labile with aggregation tendency. Addition of 0.15M NaCl increased thermostability of apo-FBP drastically, and even more so that of holo-FBP. Electrostatic forces thus seem to contribute to a diminished thermostability at low pH. Fluorescence spectroscopy after irreversible thermal unfolding of FBP revealed a weak-affinity folate binding.

Force Field for Peptides and Proteins based on the Classical Drude Oscillator

Presented is a polarizable force field based on a classical Drude oscillator framework, currently implemented in the programs CHARMM and NAMD, for modeling and molecular dynamics (MD) simulation studies of peptides and proteins. Building upon parameters for model compounds representative of the functional groups in proteins, the development of the force field focused on the optimization of the parameters for the polypeptide backbone and the connectivity between the backbone and side chains. Optimization of the backbone electrostatic parameters targeted quantum mechanical conformational energies, interactions with water, molecular dipole moments and polarizabilities and experimental condensed phase data for short polypeptides such as (Ala)₅. Additional optimization of the backbone φ, ψ conformational preferences included adjustments of the tabulated two-dimensional spline function through the CMAP term. Validation of the model included simulations of a collection of peptides and proteins. This 1^st generation polarizable model is shown to maintain the folded state of the studied systems on the 100 ns timescale in explicit solvent MD simulations. The Drude model typically yields larger RMS differences as compared to the additive CHARMM36 force field (C36) and shows additional flexibility as compared to the additive model. Comparison with NMR chemical shift data shows a small degradation of the polarizable model with respect to the additive, though the level of agreement may be considered satisfactory, while for residues shown to have significantly underestimated S² order parameters in the additive model, improvements are calculated with the polarizable model. Analysis of dipole moments associated with the peptide backbone and tryptophan side chains show the Drude model to have significantly larger values than those present in C36, with the dipole moments of the peptide backbone enhanced to a greater extent in sheets versus helices and the dipoles of individual moieties observed to undergo significant variations during the MD simulations. Although there are still some limitations, the presented model, termed Drude-2013, is anticipated to yield a molecular picture of peptide and protein structure and function that will be of increased physical validity and internal consistency in a computationally accessible fashion.

Amino-acid-dependent main-chain torsion-energy terms for protein systems

Many commonly used force fields for protein systems such as AMBER, CHARMM, GROMACS, OPLS, and ECEPP have amino-acid-independent force-field parameters for main-chain torsion-energy terms. Here, we propose a new type of amino-acid-dependent torsion-energy terms in the force fields. As an example, we applied this approach to AMBER ff03 force field and determined new amino-acid-dependent parameters for ψ (N-C(α)-C-N) and ζ (C(β)-C(α)-C-N) angles for each amino acid by using our optimization method, which is one of the knowledge-based approach. In order to test the validity of the new force-field parameters, we then performed folding simulations of α-helical and β-hairpin peptides, using the optimized force field. The results showed that the new force-field parameters gave structures more consistent with the experimental implications than the original AMBER ff03 force field.

Testing the diffusing boundary model for the helix-coil transition in peptides

The dynamics of peptide α-helices have been studied extensively for many years, and the kinetic mechanism of the helix-coil dynamics has been discussed controversially. Recent experimental results have suggested that equilibrium helix-coil dynamics are governed by movement of the helix/coil boundary along the peptide chain, which leads to slower unfolding kinetics in the helix center compared with the helix ends and position-independent helix formation kinetics. We tested this diffusion of boundary model in helical peptides of different lengths by triplet-triplet energy transfer measurements and compared the data with simulations based on a kinetic linear Ising model. The results show that boundary diffusion in helical peptides can be described by a classical, Einstein-type, 1D diffusion process with a diffusion coefficient of 2.7⋅10(7) (amino acids)(2)/s or 6.1⋅10(-9) cm(2)/s. In helices with a length longer than about 40 aa, helix unfolding by coil nucleation in a helical region occurs frequently in addition to boundary diffusion. Boundary diffusion is slowed down by helix-stabilizing capping motifs at the helix ends in agreement with predictions from the kinetic linear Ising model. We further tested local and nonlocal effects of amino acid replacements on helix-coil dynamics. Single amino acid replacements locally affect folding and unfolding dynamics with a ϕf-value of 0.35, which shows that interactions leading to different helix propensities for different amino acids are already partially present in the transition state for helix formation. Nonlocal effects of amino acid replacements only influence helix unfolding (ϕf = 0) in agreement with a diffusing boundary mechanism.

The N-terminal helix controls the transition between the soluble and amyloid states of an FF domain

Background

Protein aggregation is linked to the onset of an increasing number of human nonneuropathic (either localized or systemic) and neurodegenerative disorders. In particular, misfolding of native α-helical structures and their self-assembly into nonnative intermolecular β-sheets has been proposed to trigger amyloid fibril formation in Alzheimer’s and Parkinson’s diseases.

Methods

Here, we use a battery of biophysical techniques to elucidate the conformational conversion of native α-helices into amyloid fibrils using an all-α FF domain as a model system.

Results

We show that under mild denaturing conditions at low pH this FF domain self-assembles into amyloid fibrils. Theoretical and experimental dissection of the secondary structure elements in this domain indicates that the helix 1 at the N-terminus has both the highest α-helical and amyloid propensities, controlling the transition between soluble and aggregated states of the protein.

Conclusions

The data illustrates the overlap between the propensity to form native α-helices and amyloid structures in protein segments.

Significance

The results presented contribute to explain why proteins cannot avoid the presence of aggregation-prone regions and indeed use stable α-helices as a strategy to neutralize such potentially deleterious stretches.

Replica-Permutation Method with the Suwa-Todo Algorithm beyond the Replica-Exchange Method

We propose a new method for molecular dynamics and Monte Carlo simulations, which is referred to as the replica-permutation method (RPM), to realize more efficient sampling than the replica-exchange method (REM). In RPM, not only exchanges between two replicas but also permutations among more than two replicas are performed. Furthermore, instead of the Metropolis algorithm, the Suwa-Todo algorithm is employed for replica-permutation trials to minimize its rejection ratio. We applied RPM to particles in a double-well potential energy, Met-enkephalin in a vacuum, and a C-peptide analog of ribonuclease A in explicit water. For comparison purposes, replica-exchange molecular dynamics simulations were also performed. As a result, RPM sampled not only the temperature space but also the conformational space more efficiently than REM for all systems. From our simulations of C-peptide, we obtained the α-helix structure with salt bridges between Gly2 and Arg10, which is known in experiments. Calculating its free-energy landscape, the folding pathway was revealed from an extended structure to the α-helix structure with the salt bridges. We found that the folding pathway consists of the two steps: The first step is the "salt-bridge formation step," and the second step is the "α-helix formation step."

pH-replica exchange molecular dynamics in proteins using a discrete protonation method

Protonation equilibria in biological molecules modulates structure, dynamics, and function. A pH-replica exchange molecular dynamics (pH-REMD) method is described here to improve the coupling between conformational and protonation sampling. Under a Hamiltonian replica exchange setup, conformations are swapped between two neighboring replicas, which themselves are at different pHs. The method has been validated on a series of biological systems. We applied pH-REMD to a series of model compounds, to an terminally charged ADFDA pentapeptide, and to a heptapeptide derived from the ovomucoid third domain (OMTKY3). In all of those systems, the predicted pK(a) by pH-REMD is very close to the experimental value and almost identical to the ones obtained by constant pH molecular dynamics (CpH MD). The method presented here, pH-REMD, has the advantage of faster convergence properties due to enhanced sampling of both conformation and protonation spaces.

The amyloid precursor protein copper binding domain histidine residues 149 and 151 mediate APP stability and metabolism

One of the key pathological hallmarks of Alzheimer disease (AD) is the accumulation of the APP-derived amyloid β peptide (Aβ) in the brain. Altered copper homeostasis has also been reported in AD patients and is thought to increase oxidative stress and to contribute to toxic Aβ accumulation and regulate APP metabolism. The potential involvement of the N-terminal APP copper binding domain (CuBD) in these events has not been investigated. Based on the tertiary structure of the APP CuBD, we examined the histidine residues of the copper binding site (His(147), His(149), and His(151)). We report that histidines 149 and 151 are crucial for CuBD stability and APP metabolism. Co-mutation of the APP CuBD His(149) and His(151) to asparagine decreased APP proteolytic processing, impaired APP endoplasmic reticulum-to-Golgi trafficking, and promoted aberrant APP oligomerization in HEK293 cells. Expression of the triple H147N/H149N/H151N-APP mutant led to up-regulation of the unfolded protein response. Using recombinant protein encompassing the APP CuBD, we found that insertion of asparagines at positions 149 and 151 altered the secondary structure of the domain. This study identifies two APP CuBD residues that are crucial for APP metabolism and suggests an additional role of this domain in APP folding and stability besides its previously identified copper binding activity. These findings are of major significance for the design of novel AD therapeutic drugs targeting this APP domain.

PROTEIN FORCE-FIELD PARAMETERS OPTIMIZED WITH THE PROTEIN DATA BANK II: COMPARISONS OF FORCE FIELDS BY FOLDING SIMULATIONS OF SHORT PEPTIDES

In Paper I of this series, the formulations of the optimization method of existing force-field parameters for protein systems have been presented. We then applied it to five sets of force-field parameters, namely, AMBER parm94, AMBER parm96, AMBER parm99, CHARMM version 22, and OPLS-AA. In order to test the validity of these force fields, the folding simulations of α-helical and β-hairpin peptides have been performed with each of the original and optimized force-field parameters. We found that all five modified force-field parameters gave both α-helical and β-hairpin structures more consistent with the experimental implications than the original force fields.

Novel lynx spider toxin shares common molecular architecture with defense peptides from frog skin

A unique 30-residue cationic peptide oxyopinin 4a (Oxt 4a) was identified in the venom of the lynx spider Oxyopes takobius (Oxyopidae). Oxt 4a contains a single N-terminally located disulfide bond, Cys4-Cys10, and is structurally different from any spider toxin studied so far. According to NMR findings, the peptide is disordered in water, but assumes a peculiar torpedo-like structure in detergent micelles. It features a C-terminal amphipathic α-helical segment (body; residues 12-25) and an N-terminal disulfide-stabilized loop (head; residues 1-11), and has an unusually high density of positive charge in the head region. Synthetic Oxt 4a was produced and shown to possess strong and broad-spectrum cytolytic and antimicrobial activity. cDNA cloning showed that the peptide is synthesized in the form of a conventional prepropeptide with an acidic prosequence. Unlike other arachnid toxins, Oxt 4a exhibits striking similarity with defense peptides from the skin of ranid frogs that contain the so-called Rana-box motif (a C-terminal disulfide-enclosed loop). Parallelism or convergence is apparent on several levels: the structure, function and biosynthesis of a lynx spider toxin are mirrored by those of Rana-box peptides from frogs.

Venom peptides from solitary hunting wasps induce feeding disorder in lepidopteran larvae

The cell lytic activity and toxicity against lepidopteran larvae of 13 venom peptides (4 OdVPs and 9 EpVPs) from two solitary hunting wasps, Orancistrocerus drewseni and Eumenes pomiformis, were examined with mastoparan as a reference peptide. Of the 13 peptides, 7 were predicted to have α-helical structures that exhibit the typical character of amphipathic α-helical antimicrobial peptides. The remaining peptides exhibited coil structures; among these, EpVP5 possesses two Cys residues that form an internal disulfide bridge. All the helical peptides including mastoparan showed antimicrobial and insect cell lytic activities, whereas only two of them were hemolytic against human erythrocytes. The helical peptides induced a feeding disorder when injected into the vicinity of the head and thorax of Spodoptera exigua larvae, perhaps because their non-specific neurotoxic or myotoxic action induced cell lysis. At low concentrations, however, these helical peptides increased cell permeability without inducing cell lysis. These findings suggest that the helical venom peptides may function as non-specific neurotoxins or myotoxins and venom-spreading factors at low concentrations, as well as preservatives for long-term storage of the prey via antimicrobial, particularly antifungal, activities.

Short-lived alpha-helical intermediates in the folding of beta-sheet proteins

Several lines of evidence point strongly toward the importance of highly alpha-helical intermediates in the folding of all globular proteins, regardless of their native structure. However, experimental refolding studies demonstrate no observable alpha-helical intermediate during refolding of some beta-sheet proteins and have dampened enthusiasm for this model of protein folding. In this study, beta-sheet proteins were hypothesized to have potential to form amphiphilic helices at a period of <3.6 residues/turn that matches or exceeds the potential at 3.6 residues/turn. Hypothetically, such potential is the basis for an effective and unidirectional mechanism by which highly alpha-helical intermediates might be rapidly disassembled during folding and potentially accounts for the difficulty in detecting highly alpha-helical intermediates during the folding of some proteins. The presence of this potential was confirmed, indicating that a model entailing ubiquitous formation of alpha-helical intermediates during the folding of globular proteins predicts previously unrecognized features of primary structure. Further, the folding of fatty acid binding protein, a predominantly beta-sheet protein that exhibits no apparent highly alpha-helical intermediate during folding, was dramatically accelerated by 2,2,2-trifluoroethanol, a solvent that stabilizes alpha-helical structure. This observation suggests that formation of an alpha-helix can be a rate-limiting step during folding of a predominantly beta-sheet protein and further supports the role of highly alpha-helical intermediates in the folding of all globular proteins.

Constant pH replica exchange molecular dynamics in biomolecules using a discrete protonation model

A constant pH replica exchange molecular dynamics (REMD) method is proposed and implemented to improve coupled protonation and conformational state sampling. By mixing conformational sampling at constant pH (with discrete protonation states) with a temperature ladder, this method avoids conformational trapping. Our method was tested and applied to seven different biological systems. The constant pH REMD not only predicted pKa correctly for small, model compounds but also converged faster than constant pH molecular dynamics (MD). We further tested our constant pH REMD on a heptapeptide from ovomucoid third domain (OMTKY3). Although constant pH REMD and MD produced very close pKa values, the constant pH REMD showed its advantage in the efficiency of conformational and protonation state samplings.

Gas-phase acidities of cysteine-polyglycine peptides: the effect of the cysteine position

The sequence and conformational effects on the gas-phase acidities of peptides have been studied by using two pairs of isomeric cysteine-polyglycine peptides, CysGly(3,4)NH(2) and Gly(3,4)CysNH(2). The extended Cooks kinetic method was employed to determine the gas-phase acidities using a triple quadrupole mass spectrometer with an electrospray ionization source. The ion activation was achieved via collision-induced dissociation experiments. The deprotonation enthalpies (Delta(acid)H) were determined to be 323.9 +/- 2.5 kcal/mol (CysGly(3)NH(2)), 319.2 +/- 2.3 kcal/mol (CysGly(4)NH(2)), 333.8 +/- 2.1 kcal/mol (Gly(3)CysNH(2)), and 321.9 +/- 2.8 kcal/mol (Gly(4)CysNH(2)), respectively. The corresponding deprotonation entropies (Delta(acid)S) of the peptides were estimated. The gas-phase acidities (Delta(acid)G) were derived to be 318.4 +/- 2.5 kcal/mol (CysGly(3)NH(2)), 314.9 +/- 2.3 kcal/mol (CysGly(4)NH(2)), 327.5 +/- 2.1 kcal/mol (Gly(3)CysNH(2)), and 317.4 +/- 2.8 kcal/mol (Gly(4)CysNH(2)), respectively. Conformations and energetic information of the neutral and anionic peptides were calculated through simulated annealing (Tripos), geometry optimization (AM1), and single point energy calculations (B3LYP/6-31+G(d)), respectively. Both neutral and deprotonated peptides adopt many possible conformations of similar energies. All neutral peptides are mainly random coils. The two C-cysteine anionic peptides, Gly(3,4)(Cys-H)(-)NH(2), are also random coils. The two N-cysteine anionic peptides, (Cys-H)(-)Gly(3,4)NH(2), may exist in both random coils and stretched helices. The two N-cysteine peptides, CysGly(3)NH(2) and CysGly(4)NH(2), are significantly more acidic than the corresponding C-terminal cysteine ones, Gly(3)CysNH(2) and Gly(4)CysNH(2). The stronger acidities of the former may come from the greater stability of the thiolate anion resulting from the interaction with the helix-macrodipole, in addition to the hydrogen bonding interactions.

Intrinsic propensities of amino acid residues in GxG peptides inferred from amide I' band profiles and NMR scalar coupling constants

A reliable intrinsic propensity scale of amino acid residues is indispensable for an assessment of how local conformational distributions in the unfolded state can affect the folding of peptides and proteins. Short host-guest peptides, such as GxG tripeptides, are suitable tools for probing such propensities. To explore the conformational distributions sampled by the central amino acid residue in these motifs, we combined vibrational (IR, Raman, and VCD) with NMR spectroscopy. The data were analyzed in terms of a superposition of two-dimensional Gaussian distribution functions in the Ramachandran space pertaining to subensembles of polyproline II, beta-strand, right- and left-handed helical, and gamma-turn-like conformations. The intrinsic propensities of eight amino acid residues (x = A, V, F, L, S, E, K, and M) in GxG peptides were determined as mole fractions of these subensembles. Our results show that alanine adopts primarily (approximately 80%) a PPII-like conformation, while valine and phenylalanine were found to sample PPII and beta-strand-like conformations equally. The centers of the respective beta-strand distributions generally do not coincide with canonical values of dihedral angles of residues in parallel or antiparallel beta-strands. In fact, the distributions for most residues found in the beta-region significantly overlap the PPII-region. A comparison with earlier reported results for trivaline reveals that the terminal valines increase the beta-strand propensity of the central valine residue even further. Of the remaining investigated amino acids, methionine preferred PPII the most (0.64), and E, S, L, and K exhibit moderate (0.56-0.45) PPII propensities. Residues V, F, S, E, and L sample, to a significant extent, a region between the canonical PPII and (antiparallel) beta-strand conformations. This region coincides with the sampling reported for L and V using theoretical predictions (Tran et al. Biochemistry 2005, 44, 11369). The distributions of all investigated residues differ from coil library and computationally predicted distributions in that they do not exhibit a substantial sampling of helical conformations. We conclude that this sampling of helical conformations arises from the context dependence, for example, neighboring residues, in proteins and longer peptides, some of which is long-range.

The effect of acidic residues and amphipathicity on the lytic activities of mastoparan peptides studied by fluorescence and CD spectroscopy

Some mastoparan peptides extracted from social wasps display antimicrobial activity and some are hemolytic and cytotoxic. Although the cell specificity of these peptides is complex and poorly understood, it is believed that their net charges and their hydrophobicity contribute to modulate their biological activities. We report a study, using fluorescence and circular dichroism spectroscopies, evaluating the influence of these two parameters on the lytic activities of five mastoparans in zwitterionic and anionic phospholipid vesicles. Four of these peptides, extracted from the venom of the social wasp Polybia paulista, present both acidic and basic residues with net charges ranging from +1 to +3 which were compared to Mastoparan-X with three basic residues and net charge +4. Previous studies revealed that these peptides have moderate-to-strong antibacterial activity against Gram-positive and Gram-negative microorganisms and some of them are hemolytic. Their affinity and lytic activity in zwitterionic vesicles decrease with the net electrical charges and the dose response curves are more cooperative for the less charged peptides. Higher charged peptides display higher affinity and lytic activity in anionic vesicles. The present study shows that the acidic residues play an important role in modulating the peptides' lytic and biological activities and influence differently when the peptide is hydrophobic or when the acidic residue is in a hydrophilic peptide.