Views of helical peptides: a proposal for the position of 3(10)-helix along the thermodynamic folding pathway

Helical twists and β-turns in structures at serine-proline sequences: Stabilization of cis-proline and type VI β-turns via C-H/O interactions

Structures at serine-proline sites in proteins were analyzed using a combination of peptide synthesis with structural methods and bioinformatics analysis of the PDB. Dipeptides were synthesized with the proline derivative (2S,4S)-(4-iodophenyl)hydroxyproline [hyp(4-I-Ph)]. The crystal structure of Boc-Ser-hyp(4-I-Ph)-OMe had two molecules in the unit cell. One molecule exhibited cis-proline and a type VIa2 β-turn (BcisD). The cis-proline conformation was stabilized by a C-H/O interaction between Pro C-Hα and the Ser side-chain oxygen. NMR data were consistent with stabilization of cis-proline by a C-H/O interaction in solution. The other crystallographically observed molecule had trans-Pro and both residues in the PPII conformation. Two conformations were observed in the crystal structure of Ac-Ser-hyp(4-I-Ph)-OMe, with Ser adopting PPII in one and the β conformation in the other, each with Pro in the δ conformation and trans-Pro. Structures at Ser-Pro sequences were further examined via bioinformatics analysis of the PDB and via DFT calculations. Ser-Pro versus Ala-Pro sequences were compared to identify bases for Ser stabilization of local structures. C-H/O interactions between the Ser side-chain Oγ and Pro C-Hα were observed in 45% of structures with Ser-cis-Pro in the PDB, with nearly all Ser-cis-Pro structures adopting a type VI β-turn. 53% of Ser-trans-Pro sequences exhibited main-chain COi•••HNi+3 or COi•••HNi+4 hydrogen bonds, with Ser as the i residue and Pro as the i + 1 residue. These structures were overwhelmingly either type I β-turns or N-terminal capping motifs on α-helices or 310-helices. These results indicate that Ser-Pro sequences are particularly potent in favoring these structures. In each, Ser is in either the PPII or β conformation, with the Ser Oγ capable of engaging in a hydrogen bond with the amide N-H of the i + 2 (type I β-turn or 310-helix; Ser χ1 t) or i + 3 (α-helix; Ser χ1 g+) residue. Non-proline cis amide bonds can also be stabilized by C-H/O interactions.

Effects of an external static EF on the conformational transition of 5-HT1A receptor: A molecular dynamics simulation study

The serotonin receptor subtype 1A (5-HT1AR), one of the G-protein-coupled receptor (GPCR) family, has been implicated in several neurological conditions. Understanding the activation and inactivation mechanism of 5-HT1AR at the molecular level is critical for discovering novel therapeutics in many diseases. Recently there has been a growing appreciation for the role of external electric fields (EFs) in influencing the structure and activity of biomolecules. In this study, we used molecular dynamics (MD) simulations to examine conformational features of active states of 5-HT1AR and investigate the effect of an external static EF with 0.02 V/nm applied on the active state of 5-HT1AR. Our results showed that the active state of 5-HT1AR maintained the native structure, while the EF led to structural modifications in 5-HT1AR, particularly inducing the inward movement of transmembrane helix 6 (TM6). Furthermore, it disturbed the conformational switches associated with activation in the CWxP, DRY, PIF, and NPxxY motifs, consequently predisposing an inclination towards the inactive-like conformation. We also found that the EF led to an overall increase in the dipole moment of 5-HT1AR, encompassing TM6 and pivotal amino acids. The analyses of conformational properties of TM6 showed that the changed secondary structure and decreased solvent exposure occurred upon the EF condition. The interaction of 5-HT1AR with the membrane lipid bilayer was also altered under the EF. Our findings reveal the molecular mechanism underlying the transition of 5-HT1AR conformation induced by external EFs, which offer potential novel insights into the prospect of employing structure-based EF applications for GPCRs.

Stabilization of 310-Helices in Macrocycles Using Dominant Rotor Methodology

Stabilization of biologically relevant structural motifs has been a long-standing challenge. Here we show that atropisomeric dominant rotors can stabilize rare 310-helices in macrocycles. The target molecules were prepared using solid-phase peptide synthesis and subjected to extensive structural analysis. Molecular dynamics (MD) simulations enabled us to acquire solution structures for the target molecules, which offered evidence for stable 310-helix formation, ordinarily a metastable state. The 310-helices were shown to retain helicity after heating to 100 °C for 72 h. Moreover, the crude atropisomeric mixtures could be thermally enriched toward 310-helical macrocycles with selectivities of >20:1.

Unveiling potential PET degrading eukaryotes through in silico bioprospecting of PETases

This study addresses the environmental problem of PET plastic through in silico bioprospecting for the identification and experimental validation of novel PET degrading eukaryotes through the in silico bioprospectingI of PETases, employing a methodology that combines Hidden Markov Models (HMMs), clustering techniques, molecular docking, and dynamic simulations. A total of 424 putative PETase sequences were identified from 219 eukaryotic organisms, highlighting six sequences with low affinity energies. The Aspergillus luchuensis sequence showed the lowest Gibbs free energy and exhibited stability at different temperatures in molecular dynamics assays. Experimental validation, through a plate clearance assay and HPLC, confirmed PETase activity in three wild-type fungal strains, with A. luchuensis showing the highest efficiency. The results obtained demonstrate the effectiveness of combining computational and experimental approaches as proof of concept to discover and validate eukaryotes with PET-degrading capabilities opening new perspectives for the sustainable management of this type of waste and contributing to its environmental mitigation.

Dissecting the geometric and hydrophobic constraints of stapled peptides

Stapled peptides are a promising class of molecules with potential as highly specific probes of protein-protein interactions and as therapeutics. Hydrocarbon stapling affects the peptide properties through the interplay of two factors: enhancing the overall hydrophobicity and constraining the conformational flexibility. By constructing a series of virtual peptides, we study the role of each factor in modulating the structural properties of a hydrocarbon-stapled peptide PM2, which has been shown to enter cells, engage its target Mouse Double Minute 2 (MDM2), and activate p53. Hamiltonian replica exchange molecular dynamics (HREMD) simulations suggest that hydrocarbon stapling favors helical populations of PM2 through a combination of the geometric constraints and the enhanced hydrophobicity of the peptide. To further understand the conformational landscape of the stapled peptides along the binding pathway, we performed HREMD simulations by restraining the peptide at different distances from MDM2. When the peptide approaches MDM2, the binding pocket undergoes dehydration which appears to be greater in the presence of the stapled peptide compared with the linear peptide. In the binding pocket, the helicity of the stapled peptide is increased due to the favorable interactions between the peptide residues as well as the staple and the microenvironment of the binding pocket, contributing to enhanced affinity. The dissection of the multifaceted mechanism of hydrocarbon stapling into individual factors not only deepens fundamental understanding of peptide stapling, but also provides guidelines for the design of new stapled peptides.

Structural and biological characterization of shortened derivatives of the cathelicidin PMAP-36

Cathelicidins, a family of host defence peptides in vertebrates, play an important role in the innate immune response, exhibiting antimicrobial activity against many bacteria, as well as viruses and fungi. This work describes the design and synthesis of shortened analogues of porcine cathelicidin PMAP-36, which contain structural changes to improve the pharmacokinetic properties. In particular, 20-mers based on PMAP-36 (residues 12-31) and 13-mers (residues 12-24) with modification of amino acid residues at critical positions and introduction of lipid moieties of different lengths were studied to identify the physical parameters, including hydrophobicity, charge, and helical structure, required to optimise their antibacterial activity. Extensive conformational analysis, performed by CD and NMR, revealed that the substitution of Pro25-Pro26 with Ala25-Lys26 increased the α-helix content of the 20-mer peptides, resulting in broad-spectrum antibacterial activity against Escherichia coli, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa and Staphylococcus epidermidis strains. Interestingly, shortening to just 13 residues resulted in only a slight decrease in antibacterial activity. Furthermore, two sequences, a 13-mer and a 20-mer, did not show cytotoxicity against HaCat cells up to 64 µM, indicating that both derivatives are not only effective but also selective antimicrobial peptides. In the short peptide, the introduction of the helicogenic α-aminoisobutyric acid forced the helix toward a prevailing 310 structure, allowing the antimicrobial activity to be maintained. Preliminary tests of resistance to Ser protease chymotrypsin indicated that this modification resulted in a peptide with an increased in vivo lifespan. Thus, some of the PMAP-36 derivatives studied in this work show a good balance between chain length, antibacterial activity, and selectivity, so they represent a good starting point for the development of even more effective and proteolysis-resistant active peptides.

A class of peptides designed to replicate and enhance the Receptor for Hyaluronic Acid Mediated Motility binding domain

The extra-cellular matrix (ECM) is a complex and rich microenvironment that is exposed and over-expressed across several injury or disease pathologies. Biomaterial therapeutics are often enriched with peptide binders to target the ECM with greater specificity. Hyaluronic acid (HA) is a major component of the ECM, yet to date, few HA adherent peptides have been discovered. A class of HA binding peptides was designed using B(X₇)B hyaluronic acid binding domains inspired from the helical face of the Receptor for Hyaluronic Acid Mediated Motility (RHAMM). These peptides were bioengineered using a custom alpha helical net method, allowing for the enrichment of multiple B(X₇)B domains and the optimisation of contiguous and non-contiguous domain orientations. Unexpectedly, the molecules also exhibited the behaviour of nanofiber forming self-assembling peptides and were investigated for this characteristic. Ten 23-27 amino acid residue peptides were assessed. Simple molecular modelling was used to depict helical secondary structures. Binding assays were performed with varying concentrations (1-10 mg/mL) and extra-cellular matrices (HA, collagens I-IV, elastin, and Geltrex). Concentration mediated secondary structures were assessed using circular dichroism (CD), and higher order nanostructures were visualized using transmission electron microscopy (TEM). All peptides formed the initial apparent 3₁₀/alpha-helices, yet peptides 17x-3, 4, BHP3 and BHP4 were HA specific and potent (i.e., a significant effect) binders at increasing concentrations. These peptides shifted from apparent 3₁₀/alpha-helical structures at low concentration to beta-sheets at increasing concentration and also formed nanofibers which are noted as self-assembling structures. Several of the HA binding peptides outperformed our positive control (mPEP35) at 3-4 times higher concentrations, and were enhanced by self-assembly as each of these groups had observable nanofibers. STATEMENT OF SIGNIFICANCE: Specific biomolecules or peptides have played a crucial role in developing materials or systems to deliver key drugs and therapeutics to a broad spectrum of diseases and disorders. In these diseased tissues, cells build protein/sugar networks, which are uniquely exposed and great targets to deliver drugs to. Hyaluronic acid (HA) is involved in every stage of injury and is abundant in cancer. To date, only two HA specific peptides have been discovered. In our work, we have designed a way to model and trace binding regions as they appear on the face of a helical peptide. Using this method we have created a family of peptides enriched with HA binding domains that stick with 3-4 higher affinity than those previously discovered.

Anticancer peptides mechanisms, simple and complex

Peptide therapy has started since 1920s with the advent of insulin application, and now it has emerged as a new approach in treatment of diseases including cancer. Using anti-cancer peptides (ACPs) is a promising way of cancer therapy as ACPs are continuing to be approved and arrived at major pharmaceutical markets. Traditional cancer treatments face different problems like intensive adverse effects to patient's body, cell resistance to conventional chemical drugs and in some worse cases the occurrence of cell multidrug resistance (MDR) of cancerous tissues against chemotherapy. On the other hand, there are some benefits conceived for peptides usage in treatment of diseases specifically cancer, as these compounds present favorable characteristics such as smaller size, high activity, low immunogenicity, good biocompatibility in vivo, convenient and rapid way of synthesis, amenable to sequence modification and revision and there is no limitation for the type of cargo they carry. It is possible to achieve an optimum molecular and functional structure of peptides based on previous experience and bank of peptide motif data which may result in novel peptide design. Bioactive peptides are able to form pores in cell membrane and induce necrosis or apoptosis of abnormal cells. Moreover, recent researches have focused on the tumor recognizing peptide motifs with the ability to permeate to cancerous cells with the aim of cancer treatment at earlier stages. In this strategy the most important factors for addressing cancer are choosing peptides with easy accessibility to tumor cell without cytotoxicity effect towards normal cells. The peptides must also meet acceptable pharmacokinetic requirements. In this review, the characteristics of peptides and cancer cells are discussed. The various mechanisms of peptides' action proposed against cancer cells make the next part of discussion. It will be followed by giving information on peptides application, various methods of peptide designing along with introducing various databases. Future aspects of peptides for employing in area of cancer treatment come as conclusion at the end.

Free-energy landscapes and insertion pathways for peptides in membrane environment

Free-energy landscapes for short peptides-specifically for variants of the pH low insertion peptide (pHLIP)-in the heterogeneous environment of a lipid bilayer or cell membrane are constructed, taking into account a set of dominant interactions and the conformational preferences of the peptide backbone. Our methodology interprets broken internal H-bonds along the backbone of a polypeptide as statistically interacting quasiparticles, activated from the helix reference state. The favored conformation depends on the local environment (ranging from polar to nonpolar), specifically on the availability of external H-bonds (with H_{2}O molecules or lipid headgroups) to replace internal H-bonds. The dominant side-chain contribution is accounted for by residue-specific transfer free energies between polar and nonpolar environments. The free-energy landscape is sensitive to the level of pH in the aqueous environment surrounding the membrane. For high pH, we identify pathways of descending free energy that suggest a coexistence of membrane-adsorbed peptides with peptides in solution. A drop in pH raises the degree of protonation of negatively charged residues and thus increases the hydrophobicity of peptide segments near the C terminus. For low pH, we identify insertion pathways between the membrane-adsorbed state and a stable trans-membrane state with the C terminus having crossed the membrane.

From Folding to Assembly: Functional Supramolecular Architectures of Peptides Comprised of Non-Canonical Amino Acids

The engineering of biological molecules is the fundamental concept behind the design of complex materials with desirable functions. Over the last few decades, peptides and proteins have emerged as useful building blocks for well-defined nanostructures with controlled size and dimensions. Short peptides in particular have received much attention due to their inherent biocompatibility, lower synthetic cost, and ease of tunability. In addition to the diverse self-assembling properties of short peptides comprising coded amino acids and their emerging applications in nanotechnology, there is now growing interest in the properties of peptides composed of non-canonical amino acids. Such non-natural oligomers have been shown in recent years to form well-defined secondary structures similar to natural proteins, with the ability to self-assemble to generate a wide variety of nanostructures with excellent biostability. This review describes recent events in the development of supramolecular assemblies of peptides composed completely of non-coded amino acids and their hybrid analogues. Special attention is paid to understanding the supramolecular assemblies at the atomic level and to considering their potential applications in nanotechnology.

A highly conserved 310 helix within the kinesin motor domain is critical for kinesin function and human health

KIF1A is a critical cargo transport motor within neurons. More than 100 known mutations result in KIF1A-associated neurological disorder (KAND), a degenerative condition for which there is no cure. A missense mutation, P305L, was identified in children diagnosed with KAND, but the molecular basis for the disease is unknown. We find that this conserved residue is part of an unusual 3₁₀ helix immediately adjacent to the family-specific K-loop, which facilitates a high microtubule-association rate. We find that the mutation negatively affects several biophysical parameters of the motor. However, the microtubule-association rate of the motor is most markedly affected, revealing that the presence of an intact K-loop is not sufficient for its function. We hypothesize that the 3₁₀ helix facilitates a specific K-loop conformation that is critical for its function. We find that the function of this proline is conserved in kinesin-1, revealing a fundamental principle of the kinesin motor mechanism.

Time-dependent conformational analysis of ALK5-lumican complex in presence of graphene and graphene oxide employing molecular dynamics and MMPBSA calculation

Lumican, an extracellular matrix protein avails wound healing by binding to ALK5 membrane receptor (TGF-beta receptor I). Their interaction enables epithelialization and substantiates rejuvenation of injured tissue. To enrich permanence of ALK5-lumican interaction, we employed graphene and graphene oxide co-factors. Herein, this study explicates concomitancy of graphene and graphene oxide with ALK5-lumican. We performed an in silico approach involving molecular modelling, molecular docking, molecular dynamics for 200 ns, DSSP analysis and MMPBSA calculations. Results of molecular dynamics indicate cofactors influential in altering bioactive site of lumican than ALK5. Similarly, MMPBSA calculations unveiled binding energy of apoenzyme as -108.09 kcal/mol, holoenzyme (G) as -79.20 kcal/mol and holoenzyme (GO) as -114.33 kcal/mol. This concludes graphene oxide lucrative in enhancing binding energy of ALK5-lumican in holoenzyme (GO) via coil formation of Lum C₁₃ domain. In contrast, graphene reduced binding energy of ALK5-lumican in holoenzyme (G) modifying Lum C₁₃ into beta sheets. MMPBSA residual contribution analysis of Lum C₁₃ residues revealed binding energy of -13.9 kcal/mol for apoenzyme, -6.8 kcal/mol for holoenzyme (G) and -19.5 kcal/mol for holoenzyme (GO). This supports coil formation propitious for better ALK5-Lum interaction. Highest SASA energy of -21.05 kcal/mol of holoenzyme (G) assures graphene reasonable for improved ALK5-lumican hydrophobicity. As per the motive of the study, graphene oxide enriches permanence of ALK5-lumican. This provides counsel for plausible exploitation of lumican and graphene oxide as targeted/nano drug delivery system to reinstate acute wounds, chronic wounds, corneal wounds, hypertrophic scars and keloids in near future. Communicated by Ramaswamy H. Sarma.

Non-canonical Shedding of TNFα by SPPL2a Is Determined by the Conformational Flexibility of Its Transmembrane Helix

Ectodomain (EC) shedding defines the proteolytic removal of a membrane protein EC and acts as an important molecular switch in signaling and other cellular processes. Using tumor necrosis factor (TNF)α as a model substrate, we identify a non-canonical shedding activity of SPPL2a, an intramembrane cleaving aspartyl protease of the GxGD type. Proline insertions in the TNFα transmembrane (TM) helix strongly increased SPPL2a non-canonical shedding, while leucine mutations decreased this cleavage. Using biophysical and structural analysis, as well as molecular dynamic simulations, we identified a flexible region in the center of the TNFα wildtype TM domain, which plays an important role in the processing of TNFα by SPPL2a. This study combines molecular biology, biochemistry, and biophysics to provide insights into the dynamic architecture of a substrate's TM helix and its impact on non-canonical shedding. Thus, these data will provide the basis to identify further physiological substrates of non-canonical shedding in the future.

Atomistic Simulations of Heme Dissociation Pathways in Human Methemoglobins Reveal Hidden Intermediates

Heme dissociations disrupt function and structural integrity of human hemoglobin and trigger various cardiovascular complications. These events become significant in methemoglobins that have undergone autoxidation of ferrous into ferric heme. We have structurally characterized the heme disassociation pathways for adult tetrameric methemoglobins using all-atom molecular dynamics simulations. These reveal that bis-histidine hemichromes, characterized here by the coordination of heme iron to both the F8 (proximal) and E7 (distal) histidines, are seen as intermediates following dissociation of the water molecule distally bound to each heme iron. Later, the breaking of coordination between heme iron and proximal histidine disrupts the F helix and pushes it away from the heme cavity, enabling both bulk solvent penetration and disruption of tetramer interface interactions. The interactions inhibiting heme dissociation were then seen to be (i) either a direct or a water-molecule-mediated interaction between distal histidine and heme iron and (ii) stacking between heme and the αCE1/βCD1 phenylalanine residue. These interactions are less important in the β than in α subunits due to a more flexible β subunit CE loop region. The absence of a distal histidine interaction in the H(E7)L mutant and increased heme cavity volume in the V(E11)A mutant both promoted heme escape from the protein interior. Adult and fetal hemoglobins were seen to share a general heme disassociation pathway and intermediates due to the conservation of key heme pocket residues. The intermediates seen here are analyzed in light of experimental studies of heme dissociation and pathways of certain hemoglobinopathies.

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.

Structural Impact of Phosphorylation and Dielectric Constant Variation on Synaptotagmin's IDR

We used time-resolved Förster resonance energy transfer, circular dichroism, and molecular dynamics simulation to investigate the structural dependence of synaptotagmin 1's intrinsically disordered region (IDR) on phosphorylation and dielectric constant. We found that a peptide corresponding to the full-length IDR sequence, a ∼60-residue strong polyampholyte, can sample structurally collapsed states in aqueous solution, consistent with its κ-predicted behavior, where κ is a sequence-dependent parameter that is used to predict IDR compaction. In implicit solvent simulations of this same sequence, lowering the dielectric constant to more closely mimic the environment near a lipid bilayer surface promoted further sampling of collapsed structures. We then examined the structural tendencies of central region residues of the IDR in isolation. We found that the exocytosis-modulating phosphorylation of Thr¹¹² disrupts a local disorder-to-order transition induced by trifluoroethanol/water mixtures that decrease the solution dielectric constant and stabilize helical structure. Implicit solvent simulations on these same central region residues testing the impact of dielectric constant alone converge on a similar result, showing that helical structure is formed with higher probability at a reduced dielectric. In these helical conformers, lysine-aspartic acid salt bridges contribute to stabilization of transient secondary structure. In contrast, phosphorylation results in formation of salt bridges unsuitable for helix formation. Collectively, these results suggest a model in which phosphorylation and compaction of the IDR sequence regulate structural transitions that in turn modulate neuronal exocytosis.

Interaction of Halictine-Related Antimicrobial Peptides with Membrane Models

We have investigated structural changes of peptides related to antimicrobial peptide Halictine-1 (HAL-1) induced by interaction with various membrane-mimicking models with the aim to identify a mechanism of the peptide mode of action and to find a correlation between changes of primary/secondary structure and biological activity. Modifications in the HAL-1 amino acid sequence at particular positions, causing an increase of amphipathicity (Arg/Lys exchange), restricted mobility (insertion of Pro) and consequent changes in antimicrobial and hemolytic activity, led to different behavior towards model membranes. Secondary structure changes induced by peptide-membrane interaction were studied by circular dichroism, infrared spectroscopy, and fluorescence spectroscopy. The experimental results were complemented by molecular dynamics calculations. An α-helical structure has been found to be necessary but not completely sufficient for the HAL-1 peptides antimicrobial action. The role of alternative conformations (such as β-sheet, PPII or 3₁₀-helix) also seems to be important. A mechanism of the peptide mode of action probably involves formation of peptide assemblies (possibly membrane pores), which disrupt bacterial membrane and, consequently, allow membrane penetration.

Helical Structure of Recombinant Melittin

Melittin is an extensively studied, 26-residue toxic peptide from honey bee venom. Because of its versatility in adopting a variety of secondary (helix or coil) and quaternary (monomer or tetramer) structures in various environments, melittin has been the focus of numerous investigations as a model peptide in protein folding studies as well as in studies involving binding to proteins, lipids, and polysaccharides. A significant body of evidence supports the view that melittin binds to these macromolecules in a predominantly helical conformation, but detailed structural knowledge of this conformation is lacking. In this report, we present nuclear magnetic resonance (NMR)-based structural insights into the helix formation of recombinant melittin in the presence of trifluoroethanol (TFE): a known secondary structure inducer in peptides. These studies were performed at neutral pH, with micromolar amounts of the peptide. Using nuclear Overhauser effect (NOE)-derived distance restraints from three-dimensional NMR spectra, we determined the atomic resolution solution NMR structure of recombinant melittin bearing a TFE-stabilized helix. To circumvent the complications with structure determination of small peptides with high conformational flexibility, we developed a workflow for enhancing proton NOEs by increasing the viscosity of the medium. In the TFE-containing medium, recombinant monomeric melittin forms a long, continuous helical structure, which consists of the N- and C-terminal α-helices and the noncanonical 3₁₀-helix in the middle. The noncanonical 3₁₀-helix is missing in the previously solved X-ray structure of tetrameric melittin and the NMR structure of melittin in methanol. Melittin's structure in TFE-containing medium provides insights into melittin's conformational transitions, which are relevant to the peptide's interactions with its biological targets.

Stabilized coronavirus spikes are resistant to conformational changes induced by receptor recognition or proteolysis

Severe acute respiratory syndrome coronavirus (SARS-CoV) emerged in 2002 as a highly transmissible pathogenic human betacoronavirus. The viral spike glycoprotein (S) utilizes angiotensin-converting enzyme 2 (ACE2) as a host protein receptor and mediates fusion of the viral and host membranes, making S essential to viral entry into host cells and host species tropism. As SARS-CoV enters host cells, the viral S is believed to undergo a number of conformational transitions as it is cleaved by host proteases and binds to host receptors. We recently developed stabilizing mutations for coronavirus spikes that prevent the transition from the pre-fusion to post-fusion states. Here, we present cryo-EM analyses of a stabilized trimeric SARS-CoV S, as well as the trypsin-cleaved, stabilized S, and its interactions with ACE2. Neither binding to ACE2 nor cleavage by trypsin at the S1/S2 cleavage site impart large conformational changes within stabilized SARS-CoV S or expose the secondary cleavage site, S2′.

Chain length, temperature and solvent effects on the structural properties of α-aminoisobutyric acid homooligopeptides

Non-coded α-amino acids, originally exploited by nature, have been successfully reproduced by recent synthetic strategies to confer special structural and functional properties to small peptides. The most known and well-studied atypical residue is α-aminoisobutyric acid (Aib), which is contained in a fairly large number of peptides with known antibiotic effects. Here, we report on a molecular dynamics (MD) study of a series of homooligopeptides based on α-aminoisobutyric acid (Aib) with increasing length (Ac-(Aib)n-NMe, n = 5, 6, 7 and 10) and at various temperatures, employing a recent extension of the AMBER force field tailored for the Aib residue. Solvent effects have been analyzed by comparative MD simulations of a heptapeptide in water and dimethylsulfoxide at different temperatures. Our results show that the preference for the 310- and/or α-helix structures, which typically characterize Aib based peptides, is finely tuned by several factors including the chain length, temperature and solvent nature. While the transitions between intra-molecular i → i + 3 and i → i + 4 hydrogen bonds characterizing 310 and α-helices, respectively, are rather fast in small peptides (in the picosecond timescale), our analysis shows that the above physical and chemical factors modulate the relative equilibrium populations of the two helical structures. The obtained results nicely agree with available experimental data and support the use of the new force field for modeling Aib containing peptides.

Structural and mechanistic insights into an archaeal DNA-guided Argonaute protein

Argonaute (Ago) proteins in eukaryotes are known as key players in post-transcriptional gene silencing¹, while recent studies on prokaryotic Agos hint at their role in the protection against invading DNA^2,3. Here, we present crystal structures of the apo enzyme and a binary Ago-guide complex of the archaeal Methanocaldococcus jannaschii (Mj) Ago. Binding of a guide DNA leads to large structural rearrangements. This includes the structural transformation of a hinge region containing a switch helix, which has been shown for human Ago2 to be critical for the dynamic target search process^4-6. To identify key residues crucial for MjAgo function, we analysed the effect of several MjAgo mutants. We observe that the nature of the 3' and 5' nucleotides in particular, as well as the switch helix, appear to impact MjAgo cleavage activity. In summary, we provide insights into the molecular mechanisms that drive DNA-guided DNA silencing by an archaeal Ago.

The activity of prolactin releasing peptide correlates with its helicity

The prolactin releasing peptide (PrRP) is involved in regulating food intake and body weight homeostasis, but molecular details on the activation of the PrRP receptor remain unclear. C-terminal segments of PrRP with 20 (PrRP20) and 13 (PrRP8-20) amino acids, respectively, have been suggested to be fully active. The data presented herein indicate this is true for the wildtype receptor only; a 5-10-fold loss of activity was found for PrRP8-20 compared to PrRP20 at two extracellular loop mutants of the receptor. To gain insight into the secondary structure of PrRP, we used CD spectroscopy performed in TFE and SDS. Additionally, previously reported NMR data, combined with ROSETTANMR, were employed to determine the structure of amidated PrRP20. The structural ensemble agrees with the spectroscopic data for the full-length peptide, which exists in an equilibrium between α- and 3(10)-helix. We demonstrate that PrRP8-20's reduced propensity to form an α-helix correlates with its reduced biological activity on mutant receptors. Further, distinct amino acid replacements in PrRP significantly decrease affinity and activity but have no influence on the secondary structure of the peptide. We conclude that formation of a primarily α-helical C-terminal region of PrRP is critical for receptor activation.

High-resolution structural characterization of Noxa, an intrinsically disordered protein, by microsecond molecular dynamics simulations

High-resolution characterization of the structure and dynamics of intrinsically disordered proteins (IDPs) remains a challenging task. Consequently, a detailed understanding of the structural and functional features of IDPs remains limited, as very few full-length disordered proteins have been structurally characterized. We have performed microsecond-long molecular dynamics (MD) simulations of Noxa, the smallest member of the large Bcl-2 family of apoptosis regulating proteins, to characterize in atomic-level detail the structural features of a disordered protein. A 2.5 μs MD simulation starting from an unfolded state of the protein revealed the formation of a central antiparallel β-sheet structure flanked by two disordered segments at the N- and C-terminal ends. This topology is in reasonable agreement with protein disorder predictions and available experimental data. We show that this fold plays an essential role in the intracellular function and regulation of Noxa. We demonstrate that unbiased MD simulations in combination with a modern force field reveal structural and functional features of disordered proteins at atomic-level resolution.

Aqueous self-assembly of short hydrophobic peptides containing norbornene amino acid into supramolecular structures with spherical shape

The preparation and self-assembly of short hydrophobic peptides containing the non-coded norbornene amino acid is reported. The formation of a supramolecular assembly in water was assessed by TEM and DLS.

Solution synthesis, conformational analysis, and antimicrobial activity of three alamethicin F50/5 analogs bearing a trifluoroacetyl label

We prepared, by solution-phase methods, and fully characterized three analogs of the membrane-active peptaibiotic alamethicin F50/5, bearing a single trifluoroacetyl (Tfa) label at the N-terminus, at position 9 (central region) or at position 19 (C-terminus), and with the three Gln at positions 7, 18, and 19 replaced by Glu(OMe) residues. To add the Tfa label at position 9 or 19, a γ-trifluoroacetylated α,γ-diaminobutyric acid (Dab) residue was incorporated as a replacement for the original Val(9) or Glu(OMe)(19) amino acid. We performed a detailed conformational analysis of the three analogs (using FT-IR absorption, CD, 2D-NMR, and X-ray diffraction), which clearly showed that Tfa labeling does not introduce any dramatic backbone modification in the predominantly α-helical structure of the parent peptaibiotic. The results of an initial solid-state (19)F-NMR study on one of the analogs favor the conclusion that the Tfa group is a very promising reporter for the analysis of peptaibioticmembrane interactions. Finally, we found that the antimicrobial activities of the three newly synthesized analogs depend on the position of the Tfa label in the peptide sequence.

SAXS data based global shape analysis of trigger factor (TF) proteins from E. coli, V. cholerae, and P. frigidicola: resolving the debate on the nature of monomeric and dimeric forms

Dimerization of bacterial chaperone trigger factor (TF) is an inherent protein concentration based property which available biophysical characterization and crystal structures have kept debatable. We acquired small-angle X-ray scattering (SAXS) intensity data from different TF homologues from Escherichia coli (ECTF), Vibrio cholerae (VCTF), and Psychrobacter frigidicola (PFTF) while varying each protein concentration. We found that ECTF and VCTF adopt a compact dimeric shape at higher concentrations which did not resemble the "back-to-back" conformation reported earlier for ECTF from crystallography (PDB ID: 1W26 ). In contrast, PFTF remained monomeric throughout the concentration range 2-90 μM displaying a multimodal open extended conformation. OLIGOMER analysis showed that both the ECTF and VCTF remained completely monomeric at lower concentrations (2-11 μM), while, at higher concentrations (60-90 μM), they adopted a dimeric form. Interestingly, the equilibrium existed in the medium concentration range (>11 and <60 μM), which correlates with the physiological concentration (40-50 μM) of TF in cell cytoplasm. Additionally, circular dichroism data revealed that solution structures of ECTF and VCTF contain predominantly α-helical content, while PFTF contains 310-helical content.

1H-Azepine-4-amino-4-carboxylic acid: a new α,α-disubstituted ornithine analogue capable of inducing helix conformations in short Ala-Aib pentapeptides

A very efficient synthesis of orthogonally protected 1H-azepine-4-amino-4-carboxylic acid, abbreviated as Azn, a conformationally restricted analogue of ornithine, was realized. It was obtained on a gram scale in good overall yield in five steps, three of which did not require isolation of the intermediates, starting from the readily available 1-amino-4-oxo-cyclohexane-4-carboxylic acid. Both enantiomers were used for the preparation of pentapeptide models containing Ala, Aib, and Azn. Conformational studies using both spectroscopic techniques (NMR, CD) and molecular dynamics on model 5-mer peptides showed that the (R)-Azn isomer possesses a marked helicogenic effect.

The critical main-chain length for helix formation in water: determined in a peptide series with alternating Aib and Ala residues exclusively and detected with ECD spectroscopy

Critical main-chain length for peptide helix formation in the crystal (solid) state and in organic solvents has been already reported. In this short communication, we describe our results aiming at assessing the aforementioned parameter in water solution. To this goal, we synthesized step-by-step by solution procedures a complete series of N-terminally acetylated, C-terminally methoxylated oligopeptides, characterized only by alternating Aib and Ala residues, from the dimer to the nonamer level. All these compounds were investigated by electronic circular dichroism in the far-UV region in water solution as a function of chemical structure, namely presence/absence of an ester moiety or a negative charge at the C-terminus, and temperature. We find that the critical main-chain lengths for 3(10)- and α-helices, although still formed to a limited extent, in aqueous solution are six and eight residues, respectively.

Lateral sorting in model membranes by cholesterol-mediated hydrophobic matching

Theoretical studies predict hydrophobic matching between transmembrane domains of proteins and bilayer lipids to be a physical mechanism by which membranes laterally self-organize. We now experimentally study the direct consequences of mismatching of transmembrane peptides of different length with bilayers of different thicknesses at the molecular level. In both model membranes and simulations we show that cholesterol critically constrains structural adaptations at the peptide-lipid interface under mismatch. These constraints translate into a sorting potential and lead to selective lateral segregation of peptides and lipids according to their hydrophobic length.

Structural insights into the oligomerization and architecture of eukaryotic membrane pore-forming toxins

Pore-forming toxins (PFTs) are proteins that are secreted as soluble molecules and are inserted into membranes to form oligomeric transmembrane pores. In this paper, we report the crystal structure of Fragaceatoxin C (FraC), a PFT isolated from the sea anemone Actinia fragacea, at 1.8 Å resolution. It consists of a crown-shaped nonamer with an external diameter of about 11.0 nm and an internal diameter of approximately 5.0 nm. Cryoelectron microscopy studies of FraC in lipid bilayers reveal the pore structure that traverses the membrane. The shape and dimensions of the crystallographic oligomer are fully consistent with the membrane pore. The FraC structure provides insight into the interactions governing the assembly process and suggests the structural changes that allow for membrane insertion. We propose a nonameric pore model that spans the membrane by forming a lipid-free α-helical bundle pore.

Structural and thermodynamic investigations on the aggregation and folding of acylphosphatase by molecular dynamics simulations and solvation free energy analysis

Protein engineering method to study the mutation effects on muscle acylphosphatase (AcP) has been actively applied to describe kinetics and thermodynamics associated with AcP aggregation as well as folding processes. Despite the extensive mutation experiments, the molecular origin and the structural motifs for aggregation and folding kinetics as well as thermodynamics of AcP have not been rationalized at the atomic resolution. To this end, we have investigated the mutation effects on the structures and thermodynamics for the aggregation and folding of AcP by using the combination of fully atomistic, explicit-water molecular dynamics simulations, and three-dimensional reference interaction site model theory. The results indicate that the A30G mutant with the fastest experimental aggregation rate displays considerably decreased α1-helical contents as well as disrupted hydrophobic core compared to the wild-type AcP. Increased solvation free energy as well as hydrophobicity upon A30G mutation is achieved due to the dehydration of hydrophilic side chains in the disrupted α1-helix region of A30G. In contrast, the Y91Q mutant with the slowest aggregation rate shows a non-native H-bonding network spanning the mutation site to hydrophobic core and α1-helix region, which rigidifies the native state protein conformation with the enhanced α1-helicity. Furthermore, Y91Q exhibits decreased solvation free energy and hydrophobicity compared to wild type due to more exposed and solvated hydrophilic side chains in the α1-region. On the other hand, the experimentally observed slower folding rates in both mutants are accompanied by decreased helicity in α2-helix upon mutation. We here provide the atomic-level structures and thermodynamic quantities of AcP mutants and rationalize the structural origin for the changes that occur upon introduction of those mutations along the AcP aggregation and folding processes.

What can solid state NMR contribute to our understanding of protein folding?

Complete understanding of the folding process that connects a structurally disordered state of a protein to an ordered, biochemically functional state requires detailed characterization of intermediate structural states with high resolution and site specificity. While the intrinsically inhomogeneous and dynamic nature of unfolded and partially folded states limits the efficacy of traditional X-ray diffraction and solution NMR in structural studies, solid state NMR methods applied to frozen solutions can circumvent the complications due to molecular motions and conformational exchange encountered in unfolded and partially folded states. Moreover, solid state NMR methods can provide both qualitative and quantitative structural information at the site-specific level, even in the presence of structural inhomogeneity. This article reviews relevant solid state NMR methods and their initial applications to protein folding studies. Using either chemical denaturation to prepare unfolded states at equilibrium or a rapid freezing apparatus to trap non-equilibrium, transient structural states on a sub-millisecond time scale, recent results demonstrate that solid state NMR can contribute essential information about folding processes that is not available from more familiar biophysical methods.

Toward detecting the formation of a single helical turn by 2D IR cross peaks between the amide-I and -II modes

We have combined two-dimensional infrared (2D IR) spectroscopy and isotope substitutions to reveal the vibrational couplings between a pair of amide-I and -II modes that are several residues away but directly connected through a hydrogen bond in a helical peptide. This strategy is demonstrated on a 3(10)-helical hexapeptide, Z-Aib-L-Leu-(Aib)2-Gly-Aib-OtBu, and its 13C=18O-Leu monolabeled and 13C=18O-Leu/15N-Gly bis-labeled isotopomers in CDCl3. The isotope-dependent amide-I/II cross peaks clearly show that the second and fourth peptide linkages are vibrationally coupled as they are in proximity, forming a 3(10)-helical turn. The experimental spectra are compared to simulations based on a vibrational exciton Hamiltonian model that fully takes into account the amide-I and -II modes. The amide-II local mode frequency is evaluated by a new model based on the effects of hydrogen-bond geometry and sites. Ab initio nearest-neighbor coupling maps of the amide-I/I, -I/II, -II/I and -II/II modes are generated by isotopically isolating the local modes of N-acetyl-glycine N'-methylamide (AcGlyNHMe). Longer range couplings are modeled by transition charge interactions. The effects of the capping groups are incorporated and isotope effects are analyzed based on ab initio calculations of six model compounds. The main features of the 2D IR spectra are reproduced by this modeling. The conformational sensitivity of the isotope-dependent amide-I/II cross peaks is discussed in comparison with the calculated spectra for a semiextended structure. Our experimental and theoretical study demonstrates that the combination of 2D IR and 13C=18O/15N labeling is a useful structural method for detecting helical turn formation with residue-level specificity.