The role of alpha-, 3(10)-, and pi-helix in helix-->coil transitions

Mutations in LIFR rewire the JAK/STAT signaling pathway: A study unveiling mechanistic details of Stüve-Wiedemann syndrome

Stüve-Wiedemann syndrome (SWS), a rare autosomal recessive disorder, characterized by diminutive size, curvature of the elongated bones, bent fingers, episodes of heightened body temperature, respiratory distress or periods of breath-holding, and challenges with feeding, especially causes fatality in infants. SWS is an outcome of potential missense mutations in the leukemia inhibitory factor receptor gene reflected as numerous amino acid mutations at protein level. Employing in silico tools and techniques like mutational screening with Pred_MutHTP, I-Mutant2.0, PANTHER.db, PolyPhen, to classify mutations as deleterious/destabilizing, in conjunction with experimental data analysis, P136A and S279P emerged as 'effect'-causing mutations. Pre-existing knowledge suggests, SWS progression is effectuated conformationally altered and dysfunctional LIFR, unable to bind to LIF and further form the LIF/LIFR/gp130 signalling complex. To gain functional insights into the effect of the said mutations on the wild type protein, an all-atom, explicit, solvent molecular dynamics simulation was performed following docking approaches. Consequently, referring to the RMSD, RMSF, protein dynamic network analysis, energy landscape plots and domain motion analysis, it was revealed that unbound LIFR_WT was more prone to LIF binding as usual whereas the mutants exhibited considerable domain closure to inhibit LIF binding. We conducted binding affinity analysis via MM/GBSA and dissociation constant estimation after LIFR-LIF docking and found the WT_complex to be more stable and compact as a whole when compared to the flexible mutant complexes thus being associated with SWS. Our study offers a route for understanding molecular level implications upon LIFR mutations which opens an avenue for therapeutic interventions.

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.

An in-silico receptor-pharmacophore based multistep molecular docking and simulation study to evaluate the inhibitory potentials against NS1 of DENV-2

DENV-2 strain is the most fatal and infectious of the five dengue virus serotypes. The non-structural protein NS1 encoded by its genome is the most significant protein required for viral pathogenesis and replication inside the host body. Thus, targeting the NS1 protein and designing an inhibitor to limit its stability and secretion is a propitious attempt in our fight against dengue. Four novel inhibitors are designed to target the conserved cysteine residues (C55, C313, C316, and C329) and glycosylation sites (N130 and N207) of the NS1 protein in an attempt to halt the spread of the dengue infection in the host body altogether. Numerous computer-aided drug designing techniques including molecular docking, molecular dynamics simulation, virtual screening, principal component analysis, and dynamic cross-correlation matrix were employed to determine the structural and functional activity of the NS1-inhibitor complexes. From our analysis, it was evident that the extent of structural and atomic level fluctuations of the ligand-bound protein exhibited a declining trend in contrast to unbound protein which was prominently noticeable through the RMSD, RMSF, Rg, and SASA graphs. The ADMET analysis of the four ligands revealed a promising pharmacokinetics and pharmacodynamic profile, along with good bioavailability and toxicity properties. The proposed drugs when bound to the targeted cavities resulted in stable conformations in comparison to their unbound state, implying they have good affinity promising effective drug action. Thus, they can be tested in vitro and used as potential anti-dengue drugs.Communicated by Ramaswamy H. Sarma.

Identification of amentoflavone as a potent SARS-CoV-2 Mpro inhibitor: a combination of computational studies and in vitro biological evaluation

Small-molecule inhibitors of SARS-CoV-2 Mpro that block the active site pocket of the viral main protease have been considered potential therapeutics for the development of drugs against SARS-CoV-2. Here, we report the identification of amentoflavone (a biflavonoid) through docking-based virtual screening of a library comprised of 231 compounds consisting of flavonoids and isoflavonoids. The docking results were further substantiated through extensive analysis of the data obtained from all-atom 150 ns MD simulation. End-state effective free energy calculations using MM-PBSA calculations further suggested that (Ra)-amentoflavone (C3'-C8''-atropisomer) may show a greater binding affinity towards the Mpro than (Sa)-amentoflavone. In vitro cytotoxicity assay established that amentoflavone showed a high CC50 value indicating much lower toxicity. Further, potent inhibition of the Mpro by amentoflavone was established by studying the effect on HEK293T cells treated with SARS-CoV-2 Mpro expressing plasmid.Communicated by Ramaswamy H. Sarma.

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.

Non-canonical helical transitions and conformational switching are associated with characteristic flexibility and disorder indices in TRP and Kv channels

Structural evidence and much experimental data have demonstrated the presence of non-canonical helical substructures (π and 3₁₀) in regions of great functional relevance both in TRP as in Kv channels. Through an exhaustive compositional analysis of the sequences underlying these substructures, we find that each of them is associated with characteristic local flexibility profiles, which in turn are implicated in significant conformational rearrangements and interactions with specific ligands. We found that α-to-π helical transitions are associated with patterns of local rigidity whereas α-to-3₁₀ transitions are mainly leagued with high local flexibility profiles. We also study the relationship between flexibility and protein disorder in the transmembrane domain of these proteins. By contrasting these two parameters, we located regions showing a sort of structural discrepancy between these similar but not identical protein attributes. Notably, these regions are presumably implicated in important conformational rearrangements during the gating in those channels. In that sense, finding these regions where flexibility and disorder are not proportional allows us to detect regions with potential functional dynamism. From this point of view, we highlighted some conformational rearrangements that occur during ligand binding events, the compaction, and refolding of the outer pore loops in several TRP channels, as well as the well-known S4 motion in Kv channels.

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.

Influence of pH on indole-dependent heterodimeric interactions between Anopheles gambiae odorant-binding proteins OBP1 and OBP4

Insect Odorant Binding Proteins (OBPs) constitute important components of their olfactory apparatus, as they are essential for odor recognition. OBPs undergo conformational changes upon pH change, altering their interactions with odorants. Moreover, they can form heterodimers with novel binding characteristics. Anopheles gambiae OBP1 and OBP4 were found capable of forming heterodimers possibly involved in the specific perception of the attractant indole. In order to understand how these OBPs interact in the presence of indole and to investigate the likelihood of a pH-dependent heterodimerization mechanism, the crystal structures of OBP4 at pH 4.6 and 8.5 were determined. Structural comparison to each other and with the OBP4-indole complex (3Q8I, pH 6.85) revealed a flexible N-terminus and conformational changes in the α4-loop-α5 region at acidic pH. Fluorescence competition assays showed a weak binding of indole to OBP4 that become further impaired at acidic pH. Additional Molecular Dynamic and Differential Scanning Calorimetry studies displayed that the influence of pH on OBP4 stability is significant compared to the modest effect of indole. Furthermore, OBP1-OBP4 heterodimeric models were generated at pH 4.5, 6.5, and 8.5, and compared concerning their interface energy and cross-correlated motions in the absence and presence of indole. The results indicate that the increase in pH may induce the stabilization of OBP4 by increasing its helicity, thereby enabling indole binding at neutral pH that further stabilizes the protein and possibly promotes the creation of a binding site for OBP1. A decrease in interface stability and loss of correlated motions upon transition to acidic pH may provoke the heterodimeric dissociation allowing indole release. Finally, we propose a potential OBP1-OBP4 heterodimer formation/disruption mechanism induced by pH change and indole binding.

Despite their structural similarities, the cytosolic isoforms of human Hsp90 show different behaviour in thermal unfolding due to their conformation: An FTIR study

Heat shock proteins 90 (Hsp90) are chaperones that promote the proper folding of other proteins under high temperature stress situations. Hsp90s are highly conserved and ubiquitous proteins, and in mammalian cells, they are localized in the cytoplasm, endoplasmic reticulum, and mitochondria. Cytoplasmic Hsp90 are named Hsp90α and Hsp90β and differ mainly in their expression pattern: Hsp90α is expressed under stress conditions, while Hsp90β is a constitutive protein. Structurally, both share the same characteristics by presenting three well-conserved domains, one of which, the N-terminal domain, has a binding site for ATP to which various drugs targeting this protein, including radicicol, can bind. The protein is mainly found in dimeric form and adopts different conformations depending on the presence of ligands, co-chaperones and client proteins. In this study, some aspects of structure and thermal unfolding of cytoplasmic human Hsp90 were analysed by infrared spectroscopy. The effect on Hsp90β of binding with a non-hydrolysable ATP analogue and radicicol was also examined. The results obtained showed that despite the high similarity in secondary structure the two isoforms exhibit substantial differences in their behaviour during thermal unfolding, as Hsp90α exhibits higher thermal stability, slower denaturation process and different event sequence during unfolding. Ligand binding strongly stabilizes Hsp90β and slightly modifies the secondary structure of the protein as well. Most likely, these structural and thermostability characteristics are closely related to the conformational cycling of the chaperone and its propensity to exist in monomer or dimer form.

Redβ177 annealase structure reveals details of oligomerization and λ Red-mediated homologous DNA recombination

The Redβ protein of the bacteriophage λ red recombination system is a model annealase which catalyzes single-strand annealing homologous DNA recombination. Here we present the structure of a helical oligomeric annealing intermediate of Redβ, consisting of N-terminal residues 1-177 bound to two complementary 27mer oligonucleotides, determined via cryogenic electron microscopy (cryo-EM) to a final resolution of 3.3 Å. The structure reveals a continuous binding groove which positions and stabilizes complementary DNA strands in a planar orientation to facilitate base pairing via a network of hydrogen bonding. Definition of the inter-subunit interface provides a structural basis for the propensity of Redβ to oligomerize into functionally significant long helical filaments, a trait shared by most annealases. Our cryo-EM structure and molecular dynamics simulations suggest that residues 133-138 form a flexible loop which modulates access to the binding groove. More than half a century after its discovery, this combination of structural and computational observations has allowed us to propose molecular mechanisms for the actions of the model annealase Redβ, a defining member of the Redβ/RecT protein family.

Redβ annealase catalyses single-strand annealing homologous DNA recombination. Here, the authors present a cryo-EM structure of a Redβ annealing intermediate bound to two complementary 27mer oligonucleotides.

The neuropeptide galanin adopts an irregular secondary structure

Human galanin is a 30-residue neuropeptide targeted for development of analgesics, antidepressants, and anticonvulsants. While previous work from our group and others has already produced significant insights into galanin's N-terminal region, no extant structures of galanin in databases include its full-length sequence and the function of its C-terminus remains ambiguous. We report the NMR solution structure of full-length human galanin C-terminal amide, determined from 2D ¹H-¹H COSY, TOCSY, and ROESY NMR data. Galanin adopts an irregular helical structure across its N-terminus, likely the average of several coiling states. We present the NMR structure of a peptide encompassing the C-terminus of galanin as a stand-alone fragment. The C-terminus of full-length galanin appears to indirectly assist the intramolecular association of hydrophobic sidechains within its N-terminus, remotely rigidifying their position when compared to previously studied N-terminal galanin fragments. By contrast, there is flexibility in the C-terminus of galanin, characterized by two i to i + 2 hydrogen-bonded turns within an otherwise dynamic backbone. The C-terminal portion of the peptide renders it soluble, and plays a hitherto undescribed biophysical role in pre-organizing the galanin receptor binding epitope. We speculate that hydrophilic microdomains of signaling peptides, hormones, and perhaps intrinsically disordered proteins may also function similarly.

Hiding in plain sight: three chemically distinct α-helix types

Linus Pauling in 1950 published a three-dimensional model for a universal protein secondary structure motif which he initially called the alpha-spiral. Jack Dunitz, then a postdoc in Pauling's lab suggested to Pauling that the term helix is more accurate than spiral when describing the right-handed peptide and protein coiled structures. Pauling agreed, hence the rise of the alpha-helix, and, by extension, the ‘double helix’ structure of DNA. Although structural biologists and protein chemists are familiar with varying polar and apolar characters of amino acids in alpha-helices, to non-experts the three chemically distinct alpha-helix types classified here may hide in plain sight.

Peptidic Scaffolds Enable Rapid and Multivariate Secondary Sphere Evolution for an Abiotic Metallocatalyst

Metalloenzymes have benefited from the iterative process of evolution to achieve the precise arrangements of secondary sphere non-covalent interactions that enhance metal-centered catalysis. Iterative synthesis of scaffolds that display complex secondary sphere elements in abiotic systems can be highly challenging and time-intensive. To overcome this synthetic bottleneck, we developed a highly modular and rapid synthetic strategy, leveraging the efficiency of solid-phase peptide synthesis and conformational control afforded by non-canonical residues to construct a ligand platform displaying up to four unique residues of varying electronics and sterics in the secondary coordination sphere. As a proof-of-concept that peptidic secondary sphere can cooperate with the metal complex, we applied this scaffold to a well-known, modestly active C-H oxidizing Fe catalyst to evolve specific non-covalent interactions that is more than double its catalytic activity. Solution-state NMR structures of several catalyst variants suggest that higher activity is correlated with a hydrophobic pocket above the Fe center that may enhance the formation of a catalyst-substrate complex. Above all, we show that peptides are a convenient, highly modular, and structurally defined ligand platform for creating secondary coordination spheres that comprise multiple, diverse functional groups.

Shaking the β-Bulges

β-bulges are irregularities inside the β-sheets. They represent more than 3 percent of the protein residues, i.e., they are as frequent as 3.₁₀ helices. In terms of evolution, β-bulges are not more conserved than any other local protein conformations within homologous protein structures. In a first of its kind study, we have investigated the dynamical behaviour of β-bulges using the largest known set of protein molecular dynamics simulations. We observed that more than 50 percent of the existing β-bulges in protein crystal structures remained stable during dynamics while more than1/6^th were not stable at all and disappeared entirely. Surprisingly, 1.1 percent of β-bulges that appeared remained stable. β-bulges have been categorized in different subtypes. The most common β-bulges' types are the smallest insertion in β-strands (namely AC and AG); they are found as stable as the whole β-bulges dataset. Low occurring types (namely PC and AS), that have the largest insertions, are significantly more stable than expected. Thus, this pioneer study allowed to precisely quantify the stability of the β-bulges, demonstrating their structural robustness, with few unexpected cases raising structural questions.

The Ca2+ response of a smart forisome protein is dependent on polymerization

Forisomes are giant self-assembling mechanoproteins that undergo reversible structural changes in response to Ca²⁺ and various other stimuli. Artificial forisomes assembled from the monomer MtSEO-F1 can be used as smart biomaterials, but the molecular basis of their functionality is not understood. To determine the role of protein polymerization in forisome activity, we tested the Ca²⁺ association of MtSEO-F1 dimers (the basic polymerization unit) by circular dichroism spectroscopy and microscale thermophoresis. We found that soluble MtSEO-F1 dimers neither associate with Ca²⁺ nor undergo structural changes. However, polarization modulation infrared reflection absorption spectroscopy revealed that aggregated MtSEO-F1 dimers and fully-assembled forisomes associate with Ca²⁺ , allowing the hydration of poorly-hydrated protein areas. A change in the signal profile of complete forisomes indicated that Ca²⁺ interacts with negatively-charged regions in the protein complexes that only become available during aggregation. We conclude that aggregation is required to establish the Ca²⁺ response of forisome polymers.

Quantitative Assessment of Chirality of Protein Secondary Structures and Phenylalanine Peptide Nanotubes

In this study we consider the features of spatial-structure formation in proteins and their application in bioengineering. Methods for the quantitative assessment of the chirality of regular helical and irregular structures of proteins are presented. The features of self-assembly of phenylalanine (F) into peptide nanotubes (PNT), which form helices of different chirality, are also analyzed. A method is proposed for calculating the magnitude and sign of the chirality of helix-like peptide nanotubes using a sequence of vectors for the dipole moments of individual peptides.

Adsorption of lysozyme on gold surfaces in the presence of an external electric potential

Adsorbed protein films consist of essential building blocks of many biotechnological and biomedical devices. The electrostatic potential may significantly modulate the protein behaviour on surfaces, affecting their structure and biological activity. In this study, lysozyme was used to investigate the effects of applied electric potentials on adsorption and the protein structure. The pH and the surface charge determine the amount and secondary structure of adsorbed lysozyme on a gold surface. In-situ measurements using polarization modulation infrared reflection absorption spectroscopy indicated that the concentration of both the adsorbed anions and the lysozyme led to conformational changes in the protein film, which was demonstrated by a greater amount of aggregated β-sheets in films fabricated at net positive charges of the Au electrode (E_ads > E_pzc). The changes in secondary structure involved two parallel processes. One comprised changes in the hydration/hydrogen-bond network at helices, leading to diverse helical structures: α-, 3₁₀- and/or π-helices. In the second process β-turns, β-sheets, and random coils displayed an ability to form aggregated β-sheet structures. The study illuminates the understanding of electrical potential-dependent changes involved in the protein misfolding process.

Energetics and structure of alanine-rich α-helices via adaptive steered molecular dynamics

The energetics and hydrogen bonding profiles of the helix-to-coil transition were found to be an additive property and to increase linearly with chain length, respectively, in alanine-rich α-helical peptides. A model system of polyalanine repeats was used to establish this hypothesis for the energetic trends and hydrogen bonding profiles. Numerical measurements of a synthesized polypeptide Ac-Y(AEAAKA)kF-NH2 and a natural α-helical peptide a2N (1-17) provide evidence of the hypothesis's generality. Adaptive steered molecular dynamics was employed to investigate the mechanical unfolding of all of these alanine-rich polypeptides. We found that the helix-to-coil transition is primarily dependent on the breaking of the intramolecular backbone hydrogen bonds and independent of specific side-chain interactions and chain length. The mechanical unfolding of the α-helical peptides results in a turnover mechanism in which a 310-helical structure forms during the unfolding, remaining at a near constant population and thereby maintaining additivity in the free energy. The intermediate partially unfolded structures exhibited polyproline II helical structure as previously seen by others. In summary, we found that the average force required to pull alanine-rich α-helical peptides in between the endpoints-namely the native structure and free coil-is nearly independent of the length or the specific primary structure.

Role of intramolecular hydrogen bonds in promoting electron flow through amino acid and oligopeptide conjugates

Elucidating the factors that control charge transfer rates in relatively flexible conjugates is of importance for understanding energy flows in biology as well as assisting the design and construction of electronic devices. Here, we report ultrafast electron transfer (ET) and hole transfer (HT) between a corrole (Cor) donor linked to a perylene-diimide (PDI) acceptor by a tetrameric alanine (Ala)₄ Selective photoexcitation of the donor and acceptor triggers subpicosecond and picosecond ET and HT. Replacement of the (Ala)₄ linker with either a single alanine or phenylalanine does not substantially affect the ET and HT kinetics. We infer that electronic coupling in these reactions is not mediated by tetrapeptide backbone nor by direct donor-acceptor interactions. Employing a combination of NMR, circular dichroism, and computational studies, we show that intramolecular hydrogen bonding brings the donor and the acceptor into proximity in a "scorpion-shaped" molecular architecture, thereby accounting for the unusually high ET and HT rates. Photoinduced charge transfer relies on a (Cor)NH^…O=C-NH^…O=C(PDI) electronic-coupling pathway involving two pivotal hydrogen bonds and a central amide group as a mediator. Our work provides guidelines for construction of effective donor-acceptor assemblies linked by long flexible bridges as well as insights into structural motifs for mediating ET and HT in proteins.

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.

A Short Peptide Hydrogel with High Stiffness Induced by 310-Helices to β-Sheet Transition in Water

Biological gels generally require polymeric chains that produce long-lived physical entanglements. Low molecular weight colloids offer an alternative to macromolecular gels, but often require ad-hoc synthetic procedures. Here, a short biomimetic peptide composed of eight amino acid residues derived from squid sucker ring teeth proteins is demonstrated to form hydrogel in water without any cross-linking agent or chemical modification and exhibits a stiffness on par with the stiffest peptide hydrogels. Combining solution and solid-state NMR, circular dichroism, infrared spectroscopy, and X-ray scattering, the peptide is shown to form a supramolecular, semiflexible gel assembled from unusual right-handed 310-helices stabilized in solution by π-π stacking. During gelation, the 310-helices undergo conformational transition into antiparallel β-sheets with formation of new interpeptide hydrophobic interactions, and molecular dynamic simulations corroborate stabilization by cross β-sheet oligomerization. The current study broadens the range of secondary structures available to create supramolecular hydrogels, and introduces 310-helices as transient building blocks for gelation via a 310-to-β-sheet conformational transition.

Polyproline chains destabilize the Alzheimer's amyloid-β protofibrils: A molecular dynamics simulation study

Alzheimer's is a fatal neurodegenerative disease for which there is no cure at present. The disease is characterized by the presence of plaques, principally comprising the amyloid-β peptide (viz., β-sheet) in the brains of a patient. In our present work, we study the interaction of these β-sheets with a different number of repeating units of proline (β-sheet breaker) by docking and all atom molecular dynamics simulations. Our results indicate that proline can break the amyloid protofibrils apart, cause them to break their β-sheet structure, and in some cases even induce the formation of 3₁₀ helices, which may be intermediates in the unfolding of these β-sheets. We have also observed that some of the important hydrogen bonds and salt bridges between chains were disrupted by proline and the tight interatomic packing of atoms in the fibrils was made relatively loose. Proline chains had a tendency to make several contacts with charged residues. Proline chains binded well to the fibrils by strong electrostatic interactions while hydrophobic interactions played a less important role. This leads to the conclusion that proline can break the amyloid fibrils apart and can be considered in the design of novel peptide-based drugs to treat Alzheimer's disease and potentially other diseases caused by the misfolding of proteins into β-sheets.

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.

Molecular interactions and mutational impact upon rhodopsin (G90→D90) for hindering dark adaptation of eye: A comparative structural level outlook for signaling mechanism in night blindness

For night blindness, a detailed structural exploration of the interactions among G-protein receptor rhodopsin, transducin and arrestin was performed. Rhodopsin is responsible for dim light vision while a point mutation (G90→D90) results in an adverse change in its photo-transduction. The validated 3D models of the three proteins were utilized, and upon mutation and interactions, rhodopsin attained higher stability (evaluated through thermodynamic energy calculations, electrostatic surface potential and solvent accessible area), thereby participating strongly with transducin. Conformational switches in mutated rhodopsin also depicted a firm conformation with few 3₁₀ helices accompanied by increased percentage of pure α-helices and sheets. All evaluations were corroborated through paired T-tests. Glu33 (glycosylated unit in the N-terminal zone) of rhodopsin plays a chief role in the overall interaction pattern. Arg69 and Glu33 from wild-type rhodopsin participated in ionic interactions, while the latter set of ionic interaction remained preserved even after mutation. Cys323 (C-terminal residue) and Arg69 formed H-bonds from the wild-type rhodopsin. Cys323 exceptionally supports cellular signaling pattern in the non-mutated situation and for the non-sufferers of night-blindness. Ser297 and Tyr43 from mutated rhodopsin reside in helices and interact with Thr32 of transducin, preserving the steady conformation in activated interacted state, even in the dark. Ser297 lies adjoined to Lys296 (retinal attachment site), which resides in NPXXY motif (an "activation switch" for signal transduction). Thus, the molecular facet for involvement of photo-transduction, which holds a paramount zone in ophthalmology, was dealt with. This might instigate the future prospect for drug discovery to prevent such mutations.

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.

Long Distance Modulation of Disorder-to-Order Transitions in Protein Allostery

Elucidation of the molecular details of allosteric communication between distant sites in a protein is key to understanding and manipulating many biological regulatory processes. Although protein disorder is acknowledged to play an important thermodynamic role in allostery, the molecular mechanisms by which this disorder is harnessed for long distance communication are known for a limited number of systems. Transcription repression by the Escherichia coli biotin repressor, BirA, is allosterically activated by binding of the small molecule effector biotinoyl-5'-AMP. The effector acts by promoting BirA dimerization, which is a prerequisite for sequence-specific binding to the biotin biosynthetic operon operator sequence. A 30 Å distance separates the effector binding and dimerization surfaces in BirA, and previous studies indicate that allostery is mediated, in part, by disorder-to-order transitions on the two coupled sites. In this work, combined experimental and computational methods have been applied to investigate the molecular basis of allosteric communication in BirA. Double-mutant cycle analysis coupled with thermodynamic measurements indicates functional coupling between residues in disordered loops on the two distant surfaces. All atom molecular dynamics simulations reveal that this coupling occurs through long distance reciprocal modulation of the structure and dynamics of disorder-to-order transitions on the two surfaces.

Characterization of Biomolecular Helices and Their Complementarity Using Geometric Analysis

A general method is presented to characterize the helical properties of potentially irregular helices, such as those found in protein secondary and tertiary structures and nucleic acids. The method was validated using artificial helices with varying numbers of points, points per helical turn, pitch, and radius. The sensitivity of the method was validated by applying increasing amounts of random perturbation to the coordinates of these helices; 399 360 helices in total were evaluated. In addition, the helical parameters of protein secondary structure elements and nucleic acid helices were analyzed. Generally, at least seven points were required to recapitulate the parameters of a helix using our method. The method can also be used to calculate the helical parameters of nucleic acid-binding proteins, like TALE, enabling direct analysis of their helix complementarity to sequence-dependent DNA distortions.

A tendril perversion in a helical oligomer: trapping and characterizing a mobile screw-sense reversal

Helical oligomers of achiral monomers adopt domains of uniform screw sense, which are occasionally interrupted by screw-sense reversals. These rare, elusive, and fast-moving features have eluded detailed characterization. We now describe the structure and habits of a screw-sense reversal trapped within a fragment of a helical oligoamide foldamer of the achiral quaternary amino acid 2-aminoisobutyric acid (Aib). The reversal was enforced by compelling the amide oligomer to adopt a right-handed screw sense at one end and a left-handed screw sense at the other. The trapped reversal was characterized by X-ray crystallography, and its dynamic properties were monitored by NMR and circular dichroism, and modelled computationally. Raman spectroscopy indicated that a predominantly helical architecture was maintained despite the reversal. NMR and computational results indicated a stepwise shift from one screw sense to another on moving along the helical chain, indicating that in solution the reversal is not localised at a specific location, but is free to migrate across a number of residues. Analogous unconstrained screw-sense reversals that are free to move within a helical structure are likely to provide the mechanism by which comparable helical polymers and foldamers undergo screw-sense inversion.

Utilizing Electron Spin Echo Envelope Modulation To Distinguish between the Local Secondary Structures of an α-Helix and an Amphipathic 310-Helical Peptide

Electron spin echo envelope modulation (ESEEM) spectroscopy was used to distinguish between the local secondary structures of an α-helix and a 3₁₀-helix. Previously, we have shown that ESEEM spectroscopy in combination with site-directed spin labeling (SDSL) and ²H-labeled amino acids (i) can probe the local secondary structure of α-helices, resulting in an obvious deuterium modulation pattern, where i+4 positions generally show larger ²H ESEEM peak intensities than i+3 positions. Here, we have hypothesized that due to the unique turn periodicities of an α-helix (3.6 residues per turn with a pitch of 5.4 Å) and a 3₁₀-helix (3.1 residues per turn with a pitch of 5.8-6.0 Å), the opposite deuterium modulation pattern would be observed for a 3₁₀-helix. In this study, ²H-labeled d₁₀-leucine (Leu) was substituted at a specific Leu residue (i) and a nitroxide spin label was positioned 2, 3, and 4 residues away (denoted i+2 to i+4) on an amphipathic model peptide, LRL₈. When LRL₈ is solubilized in trifluoroethanol (TFE), the peptide adopts an α-helical structure, and alternatively, forms a 3₁₀-helical secondary structure when incorporated into liposomes. Larger ²H ESEEM peaks in the FT frequency domain data were observed for the i+4 samples when compared to the i+3 samples for the α-helix whereas the opposite pattern was revealed for the 3₁₀-helix. These unique patterns provide pertinent local secondary structural information to distinguish between the α-helical and 3₁₀-helical structural motifs for the first time using this ESEEM spectroscopic approach with short data acquisition times (∼30 min) and small sample concentrations (∼100 μM) as well as providing more site-specific secondary structural information compared to other common biophysical approaches, such as CD.

Temperature-dependent molecular dynamics study reveals an ionic liquid induced 310 - to α-helical switch in a neurotoxin

Thermal melting and recooling of AuIB, a neurotoxic conopeptide and a highly potent nonaddictive pain reliever is investigated thoroughly in water and an ionic liquid (IL) 1-butyl-3-methylimidazolium Chloride, [Im₄₁ ][Cl] by classical molecular dynamics simulations. Structural evolution of AuIB in water and the IL is observed at different temperatures between 305 and 400 K, to explore how highly viscous ionic solvents affect the peptide structure as compared to conventional solvent water. At 305 K, unlike water, the coercive effect of IL frustrates AuIB secondary structural motifs significantly. As the temperature is raised, a very interesting IL induced conformational transition from 3₁₀ - to α-helix is noticed in the peptide, presumably triggered by a significant restructuring of the peptide H-bond network. The backbone length distributions of the peptide indicate that the IL induced conformational switching is accompanied by a reduction of the axial rise of the helical region, encompassing the residues Pro-6 to Ala-10. Further, we estimated the void space available to the peptide for its structural relaxation within the first solvation shell of ∼5 Å in water as well as in IL. A temperature increase by 100 K, opens up an estimated void volume of ∼70 ų , equivalent to the volume of approximately six water molecules, around the peptide in IL. Cooling simulations of AuIB point to the crucial interplay between thermodynamically favored AuIB conformers and their kinetic control. This study provides a comprehensive understanding of the ionic solvation of biomolecules reinforcing previous experimental findings.

How To Minimize Artifacts in Atomistic Simulations of Membrane Proteins, Whose Crystal Structure Is Heavily Engineered: β₂-Adrenergic Receptor in the Spotlight

Atomistic molecular dynamics (MD) simulations are used extensively to elucidate membrane protein properties. These simulations are based on three-dimensional protein structures that in turn are often based on crystallography. The protein structures resolved in crystallographic studies typically do not correspond to pristine proteins, however. Instead the crystallized proteins are commonly engineered, including structural modifications (mutations, replacement of protein sequences by antibodies, bound ligands, etc.) whose impact on protein structure and dynamics is largely unknown. Here we explore this issue through atomistic MD simulations (∼5 μs in total), focusing on the β2-adrenergic receptor (β2AR) that is one of the most studied members of the G-protein coupled receptor superfamily. Starting from an inactive-state crystal structure of β2AR, we remove the many modifications in β2AR systematically one at a time, in six consecutive steps. After each step, we equilibrate the system and simulate it quite extensively. The results of this step-by-step approach highlight that the structural modifications used in crystallization can affect ligand and G-protein binding sites, packing at the transmembrane-helix interface region, and the dynamics of connecting loops in β2AR. When the results of the systematic step-by-step approach are compared to an all-at-once technique where all modifications done on β2AR are removed instantaneously at the same time, it turns out that the step-by-step method provides results that are superior in terms of maintaining protein structural stability. The results provide compelling evidence that for membrane proteins whose 3D structure is based on structural engineering, the preparation of protein structure for atomistic MD simulations is a delicate and sensitive process. The results show that most valid results are found when the structural modifications are reverted slowly, one at a time.

Removal of the Side Chain at the Active-Site Serine by a Glycine Substitution Increases the Stability of a Wide Range of Serine β-Lactamases by Relieving Steric Strain

Serine β-lactamases are bacterial enzymes that hydrolyze β-lactam antibiotics. They utilize an active-site serine residue as a nucleophile, forming an acyl-enzyme intermediate during hydrolysis. In this study, thermal denaturation experiments as well as X-ray crystallography were performed to test the effect of substitution of the catalytic serine with glycine on protein stability in serine β-lactamases. Six different enzymes comprising representatives from each of the three classes of serine β-lactamases were examined, including TEM-1, CTX-M-14, and KPC-2 of class A, P99 of class C, and OXA-48 and OXA-163 of class D. For each enzyme, the wild type and a serine-to-glycine mutant were evaluated for stability. The glycine mutants all exhibited enhanced thermostability compared to that of the wild type. In contrast, alanine substitutions of the catalytic serine in TEM-1, OXA-48, and OXA-163 did not alter stability, suggesting removal of the Cβ atom is key to the stability increase associated with the glycine mutants. The X-ray crystal structures of P99 S64G, OXA-48 S70G and S70A, and OXA-163 S70G suggest that removal of the side chain of the catalytic serine releases steric strain to improve enzyme stability. Additionally, analysis of the torsion angles at the nucleophile position indicates that the glycine mutants exhibit improved distance and angular parameters of the intrahelical hydrogen bond network compared to those of the wild-type enzymes, which is also consistent with increased stability. The increased stability of the mutants indicates that the enzyme pays a price in stability for the presence of a side chain at the catalytic serine position but that the cost is necessary in that removal of the serine drastically impairs function. These findings support the stability-function hypothesis, which states that active-site residues are optimized for substrate binding and catalysis but that the requirements for catalysis are often not consistent with the requirements for optimal stability.

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.

The Physiological Molecular Shape of Spectrin: A Compact Supercoil Resembling a Chinese Finger Trap

The primary, secondary, and tertiary structures of spectrin are reasonably well defined, but the structural basis for the known dramatic molecular shape change, whereby the molecular length can increase three-fold, is not understood. In this study, we combine previously reported biochemical and high-resolution crystallographic data with structural mass spectroscopy and electron microscopic data to derive a detailed, experimentally-supported quaternary structure of the spectrin heterotetramer. In addition to explaining spectrin’s physiological resting length of ~55-65 nm, our model provides a mechanism by which spectrin is able to undergo a seamless three-fold extension while remaining a linear filament, an experimentally observed property. According to the proposed model, spectrin’s quaternary structure and mechanism of extension is similar to a Chinese Finger Trap: at shorter molecular lengths spectrin is a hollow cylinder that extends by increasing the pitch of each spectrin repeat, which decreases the internal diameter. We validated our model with electron microscopy, which demonstrated that, as predicted, spectrin is hollow at its biological resting length of ~55-65 nm. The model is further supported by zero-length chemical crosslink data indicative of an approximately 90 degree bend between adjacent spectrin repeats. The domain-domain interactions in our model are entirely consistent with those present in the prototypical linear antiparallel heterotetramer as well as recently reported inter-strand chemical crosslinks. The model is consistent with all known physical properties of spectrin, and upon full extension our Chinese Finger Trap Model reduces to the ~180-200 nm molecular model currently in common use.

Spectrins are cytoskeletal and scaffolding proteins ubiquitously expressed in essentially all cell-types. Despite unequivocal evidence for a short physiological length of ~55–65 nm at rest, spectrin is typically represented as an extended ~200 nm molecule that is implied based on crystallographic structures of a number of tandem repeats. Here, we incorporate previously reported biochemical and crystallographic data with structural mass spectroscopy and electron microscopic data to derive a detailed, experimentally-supported quaternary structure of the physiological compact form of spectrin. In addition to explaining spectrin’s physiological resting length (~55–65 nm), our model provides a mechanism by which spectrin can undergo a seamless three-fold extension, which is an experimentally observed property that is responsible for restoration of cell shape after mechanical deformation.

Native like helices in a specially designed β peptide in the gas phase

First principles simulations and gas phase spectroscopy suggest equilibrium of helices for an oligomer of open chain β amino acids.

The influence of flanking secondary structures on amino Acid content and typical lengths of 3/10 helices

We used 3D structures of a highly redundant set of bacterial proteins encoded by genes of high, average, and low GC-content. Four types of connecting bridges—regions situated between any of two major elements of secondary structure (alpha helices and beta strands)—containing a pure random coil were compared with connecting bridges containing 3/10 helices. We included discovered trends in the original “VVTAK Connecting Bridges” algorithm, which is able to predict more probable conformation for a given connecting bridge. The highest number of significant differences in amino acid usage was found between 3/10 helices containing bridges connecting two beta strands (they have increased Phe, Tyr, Met, Ile, Leu, Val, and His usages but decreased usages of Asp, Asn, Gly, and Pro) and those without 3/10 helices. The typical (most common) length of 3/10 helices situated between two beta strands and between beta strand and alpha helix is equal to 5 amino acid residues. The preferred length of 3/10 helices situated between alpha helix and beta strand is equal to 3 residues. For 3/10 helices situated between two alpha helices, both lengths (3 and 5 amino acid residues) are typical.

Structural diversity of Alzheimer's disease amyloid-β dimers and their role in oligomerization and fibril formation

Alzheimer's disease (AD) is associated with the formation of toxic amyloid-β (Aβ)42 oligomers, and recent evidence supports a role for Aβ dimers as building blocks for oligomers. Molecular dynamics simulation studies have identified clans for the dominant conformations of Aβ42 forming dimers; however, it is unclear if a larger spectrum of dimers is involved and which set(s) of dimers might evolve to oligomers verse fibrils. Therefore, for this study we generated multiple structural conformations of Aβ42, using explicit all-atom molecular dynamics, and then clustering the different structures based on key conformational similarities. Those matching a selection threshold were then used to model a process of oligomerization. Remarkably, we showed a greater diversity in Aβ dimers than previously described. Depending on the clan family, different types of Aβ dimers were obtained. While some had the tendency to evolve into oligomeric rings, others formed fibrils of diverse characteristics. Then we selected the dimers that would evolve to membranephilic annular oligomers. Nearly one third of the 28 evaluated annular oligomers had the dimer interfaces between the neighboring Aβ42 monomers with possible salt bridges between the residue K28 from one side and either residue E22 or D23 on the other. Based on these results, key amino acids were identified for point mutations that either enhanced or suppressed the formation and toxicity of oligomer rings. Our studies suggest a greater diversity of Aβ dimers. Understanding the structure of Aβ dimers might be important for the rationale design of small molecules that block formation of toxic oligomers.

Fluctuations of an exposed π-helix involved in lipoxygenase substrate recognition

The second helix in lipoxygenases adapts to permit substrate access to the active site, but details of this process are varied and poorly understood. We therefore examined the dynamics of helix 2 in solutions of spin-labeled soybean lipoxygenase-1 and spin relaxation at 60 K of the spin-labels by catalytic iron. Helix 2 in soybean lipoxygenase structures is surface-exposed and contains one turn of π-helix, centrally located. A site-directed spin-label scan of 18 of the 21 helix 2 residues, and electron paramagnetic resonance, showed that the π-helical segment became unusually mobile, on a nanosecond time scale, under conditions favoring substrate binding (pH 9 and lipid addition), while segments before and after had relatively unchanged dynamics. Backbone dynamics of residues in the π-helical segment appeared to be correlated, at pH 9. Samples also were frozen to examine the polarity and proticity of the local environments, the effect of the local environment on intrinsic relaxation, and dipolar relaxation by two symmetries of catalytic iron. The average hyperfine tensor component, Azz, of four π-helix residues decreased by 1.75 G, with an increase in pH from 7 to 9, while it remained unaffected for nearby buried residues. Power saturation data suggested the change in polarity specific to the π-helix altered the intrinsic relaxation rates. Different symmetries of iron contributed to distance-dependent magnetic relaxation. We interpret these data to mean that a π-helix in the second helix of plant lipoxygenases is highly dynamic and is the site where lipid chains penetrate to inner helices that outline the substrate pocket.

Molecular simulations of a dynamic protein complex: role of salt-bridges and polar interactions in configurational transitions

Ion charge pairs and hydrogen bonds have been extensively studied for their roles in stabilizing protein complexes and in steering the process of protein association. Recently, it has become clear that some protein complexes are dynamic in that they interconvert between several alternate configurations. We have previously characterized one such system: the EphA2:SHIP2 SAM-SAM heterodimer by solution NMR. Here we carried out extensive all-atom molecular-dynamics simulations on a microsecond time-scale starting with different NMR-derived structures for the complex. Transitions are observed between several discernible configurations at average time intervals of 50-100 ns. The domains reorient relative to one another by substantial rotation and a slight shifting of the interfaces. Bifurcated and intermediary salt-bridge and hydrogen-bond interactions play a role in the transitions in a process that can be described as moving along a "monkey-bar". We notice an increased density of salt bridges near protein interaction surfaces that appear to enable these transitions, also suggesting why the trajectories can become kinetically hindered in regions where fewer of such interactions are possible. In this context, even microsecond molecular-dynamics simulations are not sufficient to sample the energy landscape unless the structures remain close to their experimentally derived low-energy configurations.