Intramolecular autoinhibition regulates the selectivity of PRPF40A tandem WW domains for proline-rich motifs.

PRPF40A plays an important role in the regulation of pre-mRNA splicing by mediating protein-protein interactions in the early steps of spliceosome assembly. By binding to proteins at the 5´ and 3´ splice sites, PRPF40A promotes spliceosome assembly by bridging the recognition of the splices. The PRPF40A WW domains are expected to recognize proline-rich sequences in SF1 and SF3A1 in the early spliceosome complexes E and A, respectively. Here, we combine NMR, SAXS and ITC to determine the structure of the PRPF40A tandem WW domains in solution and characterize the binding specificity and mechanism for proline-rich motifs recognition. Our structure of the PRPF40A WW tandem in complex with a high-affinity SF1 peptide reveals contributions of both WW domains, which also enables tryptophan sandwiching by two proline residues in the ligand. Unexpectedly, a proline-rich motif in the N-terminal region of PRPF40A mediates intramolecular interactions with the WW tandem. Using NMR, ITC, mutational analysis in vitro, and immunoprecipitation experiments in cells, we show that the intramolecular interaction acts as an autoinhibitory filter for proof-reading of high-affinity proline-rich motifs in bona fide PRPF40A binding partners. We propose that similar autoinhibitory mechanisms are present in most WW tandem-containing proteins to enhance binding selectivity and regulation of WW/proline-rich peptide interaction networks.

Structural basis for the distinct roles of non-conserved Pro116 and conserved Tyr124 of BCH domain of yeast p50RhoGAP.

p50RhoGAP is a key protein that interacts with and downregulates the small GTPase RhoA. p50RhoGAP is a multifunctional protein containing the BNIP-2 and Cdc42GAP Homology (BCH) domain that facilitates protein-protein interactions and lipid binding and the GAP domain that regulates active RhoA population. We recently solved the structure of the BCH domain from yeast p50RhoGAP (YBCH) and showed that it maintains the adjacent GAP domain in an auto-inhibited state through the ß5 strand. Our previous WT YBCH structure shows that a unique kink at position 116 thought to be made by a proline residue between alpha helices α6 and α7 is essential for the formation of intertwined dimer from asymmetric monomers. Here we sought to establish the role and impact of this Pro116. However, the kink persists in the structure of P116A mutant YBCH domain, suggesting that the scaffold is not dictated by the proline residue at this position. We further identified Tyr124 (or Tyr188 in HBCH) as a conserved residue in the crucial ß5 strand. Extending to the human ortholog, when substituted to acidic residues, Tyr188D or Tyr188E, we observed an increase in RhoA binding and self-dimerization, indicative of a loss of inhibition of the GAP domain by the BCH domain. These results point to distinct roles and impact of the non-conserved and conserved amino acid positions in regulating the structural and functional complexity of the BCH domain.

4,4-Difluoroproline as a Unique 19F NMR Probe of Proline Conformation.

Despite the importance of proline conformational equilibria (trans versus cis amide and exo versus endo ring pucker) on protein structure and function, there is a lack of convenient ways to probe proline conformation. 4,4-Difluoroproline (Dfp) was identified to be a sensitive 19F NMR-based probe of proline conformational biases and cis-trans isomerism. Within model compounds and disordered peptides, the diastereotopic fluorines of Dfp exhibit similar chemical shifts (ΔδFF = 0-3 ppm) when a trans X-Dfp amide bond is present. In contrast, the diastereotopic fluorines exhibit a large (ΔδFF = 5-12 ppm) difference in chemical shift in a cis X-Dfp prolyl amide bond. DFT calculations, X-ray crystallography, and solid-state NMR spectroscopy indicated that ΔδFF directly reports on the relative preference of one proline ring pucker over the other: a fluorine which is pseudo-axial (i.e., the pro-4R-F in an exo ring pucker, or the pro-4S-F in an endo ring pucker) is downfield, while a fluorine which is pseudo-equatorial (i.e., pro-4S-F when exo, or pro-4R-F when endo) is upfield. Thus, when a proline is disordered (a mixture of exo and endo ring puckers, as at trans-Pro in peptides in water), it exhibits a small Δδ. In contrast, when the Pro is ordered (i.e., when one ring pucker is strongly preferred, as in cis-Pro amide bonds, where the endo ring pucker is strongly favored), a large Δδ is observed. Dfp can be used to identify inherent induced order in peptides and to quantify proline cis-trans isomerism. Using Dfp, we discovered that the stable polyproline II helix (PPII) formed in the denatured state (8 M urea) exhibits essentially equal populations of the exo and endo proline ring puckers. In addition, the data with Dfp suggested the specific stabilization of PPII by water over other polar solvents. These data strongly support the importance of carbonyl solvation and n → π* interactions for the stabilization of PPII. Dfp was also employed to quantify proline cis-trans isomerism as a function of phosphorylation and the R406W mutation in peptides derived from the intrinsically disordered protein tau. Dfp is minimally sterically disruptive and can be incorporated in expressed proteins, suggesting its broad application in understanding proline cis-trans isomerization, protein folding, and local order in intrinsically disordered proteins.

Effects of Proline on Internal Friction in Simulated Folding Dynamics of Several Alanine-Based α-Helical Peptides.

We have studied in silico the effect of proline, a model cosolvent, on local and global friction coefficients in (un)folding of several typical alanine-based α-helical peptides. Local friction is related to dwell times of a single, ensemble-averaged hydrogen bond (HB) within each peptide. Global friction is related to energy dissipated in a series of configurational changes of each peptide experienced by increasing the number of HBs during folding. Both of these approaches are important in relation to future atomic force microscopic-based measurements of internal friction via force-clamp single-molecule force spectroscopy. Molecular dynamics (MD) simulations for six peptides, namely, ALA5, ALA8, ALA15, ALA21, (AAQAA)3, and H2N-GN(AAQAA)2G-COONH2, have been conducted at 2 and 5 M proline solutions in water. Using previously obtained MD data for these peptides in pure water as well as upgraded theoretical models, we obtained variations of local and global internal friction coefficients as a function of solution viscosity. The results showed the substantial role of proline in stabilizing the folded state and slowing the overall folding dynamics. Consequently, larger friction coefficients were obtained at larger viscosities. The local and global internal friction, i.e., respective, friction coefficients approximated to zero viscosity, was also obtained. The evolution of friction coefficients with viscosity was weakly dependent on the number of concurrent folding pathways but was rather dominated by a stabilizing effect of proline on the folded states. Obtained values of local and global internal friction showed qualitatively similar results and a clear dependency on the structure of the studied peptide.

Misacylation of tRNA with Ser-Pro Dipeptide for In Vitro Transcription-Translation.

Serine-proline (Ser-Pro) backbone-modified dipeptide analogues are powerful tools to investigate the role of cis-trans isomerization in the regulation of the cell cycle and transcription. These studies have previously been limited to synthetic peptides, whose synthesis is a challenge for larger peptides due to the compounding yield loss incurred in each step. We now introduce a method for the aminoacylation of tRNA with dipeptides and dipeptide analogs to permit the installation of cis- and trans-locked Ser-Pro analogues into full-length proteins. To that end, we synthesized the 3,5-dinitrobenzyl (DNB)-activated esters of a native Ser-Pro dipeptide and its cis- and trans-locked alkene analogs. Murakami et al. created the DNB flexizyme (dFx), a ribozyme that acylates tRNA with DNB esters of amino acids to permit unnatural amino acids to be incorporated into proteins. A tRNA from yeast that recognizes the amber stop codon, along with the dFx flexizyme, were generated by in vitro transcription with T7 RNA polymerase. dFx was used to successfully catalyze the chemical misacylation of truncated amber tRNA with the Ser-Pro-DNB activated dipeptide. This method allows the introduction of non-native Ser-Pro dipeptide mimics into full-length proteins by in vitro transcription-translation. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Synthesis of 3,5-dinitrobenzyl activated esters of Ser-Pro Basic Protocol 2: Preparation of truncated amber tRNA Basic Protocol 3: Acylation of amber-tRNA by the dFx flexizyme Basic Protocol 4: PAGE electrophoresis of tRNASerPro.

Comprehensive Analysis of Coupled Proline Cis-Trans States in Bradykinin Using ωBP-REMD Simulations.

It is well-known that proline (Pro) cis-trans isomerization plays a decisive role in the folding and stabilization of proteins. The conformational coupling between isomerization states of different Pro residues in proteins during conformational adaptation processes is not well understood. In the present work, we investigate the coupled cis-trans isomerization of three Pro residues using bradykinin (BK), a partially unstructured nonapeptide hormone, as a model system. We use a recently developed enhanced-sampling molecular dynamics method (ω-bias potential replica exchange molecular dynamics; ωBP-REMD) that allows us to exhaustively sample all combinations of Pro isomer states and obtain converged probability densities of all eight state combinations within 885 ns ωBP-REMD simulations. In agreement with experiment, the all-trans state is seen to be the preferred isomer of zwitterionic aqueous BK. In about a third of its structures, this state presents the characteristic C-terminal ß-turn conformation; however, other isomer combinations also contribute significantly to the structural ensemble. Unbiased probabilities can be projected onto the peptide bond dihedral angles of the three Pro residues. This unveils the interdependence of the individual Pro isomerization states, i.e., a possible coupling of the different Pro isomers. The cis/trans equilibrium of a Pro residue can change by up to 2.5 kcal·mol-1, depending on the isomerization state of other Pro residues. For example, for Pro7, the simulations indicate that its cis state becomes favored compared to its trans state when Pro2 is switched from the trans state to the cis state. Our findings demonstrate the efficiency of the ωBP-REMD methodology and suggest that the coupling of Pro isomerization states may play an even more decisive role in larger folded proteins subject to more conformational restraints.

Conserved proline residues prevent dimerization and aggregation in the ß-lactamase BlaC.

Evolution leads to conservation of amino acid residues in protein families. Conserved proline residues are usually considered to ensure the correct folding and to stabilize the three-dimensional structure. Surprisingly, proline residues that are highly conserved in class A ß-lactamases were found to tolerate various substitutions without large losses in enzyme activity. We investigated the roles of three conserved prolines at positions 107, 226, and 258 in the ß-lactamase BlaC from Mycobacterium tuberculosis and found that mutations can lead to dimerization of the enzyme and an overall less stable protein that is prone to aggregate over time. For the variant Pro107Thr, the crystal structure shows dimer formation resembling domain swapping. It is concluded that the proline substitutions loosen the structure, enhancing multimerization. Even though the enzyme does not lose its properties without the conserved proline residues, the prolines ensure the long-term structural integrity of the enzyme.

Polymorphic Structure Determination of the Macrocyclic Drug Paritaprevir by MicroED.

Paritaprevir is an orally bioavailable, macrocyclic drug used for treating chronic Hepatitis C virus (HCV) infection. Its structures have been elusive to the public until recently when one of the crystal forms is solved by microcrystal electron diffraction (MicroED). In this work, the MicroED structures of two distinct polymorphic crystal forms of paritaprevir are reported from the same experiment. The different polymorphs show conformational changes in the macrocyclic core, as well as the cyclopropyl sulfonamide and methyl pyrazinamide substituents. Molecular docking shows that one of the conformations fits well into the active site pocket of the HCV non-structural 3/4A (NS3/4A) serine protease target, and can interact with the pocket and catalytic triad via hydrophobic interactions and hydrogen bonds. These results can provide further insight for optimization of the binding of acyl sulfonamide inhibitors to the HCV NS3/4A serine protease. In addition, this also demonstrates the opportunity to derive different polymorphs and distinct macrocycle conformations from the same experiments using MicroED.

Bioinspired Herbicides-Ionic Liquids or Liquid Cocrystals?

This research provides information about combinations of several amino acids, including l-proline (Pro), l-arginine (Arg), and l-histidine (His), with phenoxyacetic acid herbicides (MCPA and 2,4-D). Five amino acid ionic liquids (AAILs), one amino acid higher-melting salt (AAHMS), and two amino acid liquid cocrystals (AALCs) were obtained in high yields (>90%). The ionization of the six new structures was confirmed by NMR, IR, and molecular modeling. X-ray crystallography was used to definitively confirm the binding location of the mobile hydrogen. Furthermore, we propose a computational method for estimating the energy of specific hydrogen bond(s) in AAIL crystals based on the NBO and QTAIM hydrogen bond parameters obtained by model calculations. An in-depth analysis of the structures allowed to answer the question posed in the title, ionic liquids or liquid cocrystals? AAILs based on arginine and histidine were obtained. In contrast, combining proline with MCPA and 2,4-D led to AALCs. Finally, the compounds were analyzed to measure their herbicidal activity. These studies proved that the novel form of MCPA or 2,4-D improved its ability to control weeds compared to commercial formulations containing the same active ingredients.

Isomerization of proline-46 in the N-terminal tail of galectin-3 enhances T cell apoptosis via the ROS-ERK pathway.

Galectin-3 (Gal-3) is unique in the galectin family, due to the presence of a long N-terminal tail (NT) arising from its conserved carbohydrate recognition domain (CRD). Although functional significance of the NT has remained elusive, our previous studies demonstrated the importance of NT prolines to Gal-3 function. Here, we show that during the time Gal-3 stands in solution for three or more days, Gal-3 NT undergoes a slow, intra-molecular, time-dependent conformational/dynamical change associated with proline cis-trans isomerization. From initial dissolution of Gal-3 in buffer to three days in solution, Gal-3-mediated T cell apoptosis is enhanced from 23 % to 37 %. Western blotting and flow cytometry show that the enhancement occurs via the ROS-ERK pathway, and not by the PKC-ERK pathway. To assess which proline(s) is (are) responsible for this effect, we individually mutated all 14 NT prolines within the first 68 residues to alanines, and assessed their effect on ROS production. Our study shows that isomerization of P46 alone is responsible for the upregulation of ROS and T cell apoptosis. NMR studies show that this unique effect is mediated by a change in dynamic interactions between the NT and CRD F-face, which in turn leads to this change in Gal-3 function.

Dissection of the role of a Src homology 3 domain in the evolution of binding preference of paralogous proteins.

Protein-protein interactions (PPIs) drive many cellular processes. Some interactions are directed by Src homology 3 (SH3) domains that bind proline-rich motifs on other proteins. The evolution of the binding specificity of SH3 domains is not completely understood, particularly following gene duplication. Paralogous genes accumulate mutations that can modify protein functions and, for SH3 domains, their binding preferences. Here, we examined how the binding of the SH3 domains of 2 paralogous yeast type I myosins, Myo3 and Myo5, evolved following duplication. We found that the paralogs have subtly different SH3-dependent interaction profiles. However, by swapping SH3 domains between the paralogs and characterizing the SH3 domains freed from their protein context, we find that very few of the differences in interactions, if any, depend on the SH3 domains themselves. We used ancestral sequence reconstruction to resurrect the preduplication SH3 domains and examined, moving back in time, how the binding preference changed. Although the most recent ancestor of the 2 domains had a very similar binding preference as the extant ones, older ancestral domains displayed a gradual loss of interaction with the modern interaction partners when inserted in the extant paralogs. Molecular docking and experimental characterization of the free ancestral domains showed that their affinity with the proline motifs is likely not the cause for this loss of binding. Taken together, our results suggest that a SH3 and its host protein could create intramolecular or allosteric interactions essential for the SH3-dependent PPIs, making domains not functionally equivalent even when they have the same binding specificity.

4-Thiaproline accelerates the slow folding phase of proteins containing cis prolines in the native state by two orders of magnitude.

The cis/trans isomerization of peptidyl-prolyl peptide bonds is often the bottleneck of the refolding reaction for proteins containing cis proline residues in the native state. Proline (Pro) analogues, especially C4-substituted fluoroprolines, have been widely used in protein engineering to enhance the thermodynamic stability of peptides and proteins and to investigate folding kinetics. 4-thiaproline (Thp) has been shown to bias the ring pucker of Pro, to increase the cis population percentage of model peptides in comparison to Pro, and to diminish the activation energy barrier for the cis/trans isomerization reaction. Despite its intriguing properties, Thp has been seldom incorporated into proteins. Moreover, the impact of Thp on the folding kinetics of globular proteins has never been reported. In this study, we show that upon incorporation of Thp at cisPro76 into the thioredoxin variant Trx1P the half-life of the refolding reaction decreased from ~2 h to ~35 s. A dramatic acceleration of the refolding rate could be observed also for the protein pseudo wild-type barstar upon replacement of cisPro48 with Thp. Quantum chemical calculations suggested that the replacement of the Cγ H2 group by a sulfur atom in the pyrrolidine ring, might lower the barrier for cis/trans rotation due to a weakened peptide bond. The protein variants retained their thermodynamic stability upon incorporation of Thp, while the catalytic and enzymatic activities of the modified Trx1P remained unchanged. Our results show that the Pro isostere Thp might accelerate the rate of the slow refolding reaction for proteins containing cis proline residues in the native state, independent from the local structural environment.

Chiral Recognition of Diketopiperazine Containing Proline Residues by (-)-Epigallocatechin-3-O-gallate in Water.

The stoichiometry and precipitate yield of a complex of (-)-epigallocatechin-3-O-gallate (EGCg) and cyclo(Pro-Xxx) (Xxx = phenylalanine (Phe), tyrosine (Tyr)) were evaluated using integrated values of their proton signals by quantitative 1H-NMR (q NMR). It was determined to be a 1 : 1 complex of EGCg and cyclo(Pro-Xxx). The change in the chemical shift value of proton signals of cyclo(Pro-Xxx) in 1H-NMR spectra by adding standard amounts of EGCg was investigated. Differences in chemical shift values of H8α, H7αß, H8ß, H10, H9, and H3 proton signals between cyclo(L-Pro-L-Phe) and cyclo(D-Pro-D-Phe), and those of H8α, H7αß, H8ß, H10, H9, H3, and H13 proton signals between cyclo(L-Pro-L-Tyr) and cyclo(D-Pro-D-Tyr) were observed as a significant difference at 54 mmol/L of EGCg. It was found that their chirality was clearly recognized by EGCg. The significant difference in the change of the chemical shift value of H8α proton signals between cyclo(L-Pro-L-Xxx) and cyclo(D-Pro-D-Xxx) was the largest, and the difference was considered to have resulted from the difference in the ratio of extended conformer in equilibrium between folded and extended conformers. Such a significant difference in change values between cyclo(L-Pro-D-Xxx) and cyclo(D-Pro-L-Xxx) was not observed due to a rigid intramolecular CH-π interaction. EGCg did not clearly recognize the chirality of cyclo(L-Pro-D-Xxx) and cyclo(D-Pro-L-Xxx).

Effects of L-ascorbic acid (C6H8O6: Vit-C) on collagen amino acids: DFT study.

Vitamin C plays a very important role in the repair of connective tissue, especially for sports whose training causes the most damage to this tissue. Therefore, many people believe that L-ascorbic acid (C6H8O6: vitamin C) reduces the recovery time between sports exercises. The most abundant form of structural protein in the body is collagen. Collagen is characterized by a high concentration of the three amino acids glycine (Gly), proline (Pro), and hydroxyproline (Hyp), which creates its characteristic triple helix structure. Therefore, in this study, the effect of vitamin C presence on the sequence, interaction, and orientation of amino acids for collagen formation is investigated using computational simulation. This study aimed to investigate the mechanism of action of vitamin C in terms of thermodynamics and structure of the reaction. The calculations are performed using density function theory (DFT) by the base set of B3LYP/6-311++G (p,d). The results show that the presence of vitamin C is effective in the formation of collagen protein for this interaction and the mechanism of amino acid sequence (Gly-Hyp-Pro) is better in the formation of collagen protein in the presence of vitamin C. The presence of Vit-C in the formation and direction of hydroxyproline (Hyp) causes its separation from the prolyl 5-hydroxylase enzyme. In the absence of vitamin C, the reaction stops at this stage and proline cannot be converted into hydroxyproline. The computational data shows vitamin C prevents unwanted interactions and directs amino acid reactions to repair connective tissue (collagen). Therefore, vitamin C acts as a cofactor in the Prolyl 5-Hydroxylase enzyme and causes it to convert proline to hydroxyl.

Phase Separation and Fibrillization of Human Annexin A7 Are Mediated by Its Proline-Rich Domain.

Human annexin A7, a calcium- and phospholipid-binding protein, governs calcium homeostasis, plasma membrane repair, apoptosis, and tumor progression. A7 contains an N-terminal proline-rich domain (PRD; 180 residues, ∼24% prolines) that determines its functional specificity. Using microscopy and dye-binding assays, we show that recombinant A7 and its isolated PRD spontaneously phase separate into spherical condensates, which subsequently transform into ß-sheet-rich fibrils. We demonstrate that fibrillization of A7-PRD proceeds via primary nucleation and fibril-catalyzed secondary nucleation processes, as determined by chemical kinetics, providing a mechanistic basis for its amyloid assembly. This study confirms and highlights a subclass of eukaryotic PRDs prone to forming aggregates with important physiological and pathological implications.

Assignment of the disordered, proline-rich N-terminal domain of the tumour suppressor p53 protein using 1HN and 1Hα-detected NMR measurements.

Protein p53 is mostly known for playing a key role in tumour suppression, and mutations in the p53 gene are amongst the most frequent genomic events accompanying oncogenic transformation. Continuous research is conducted to target disordered proteins/protein regions for cancer therapy, for which atomic level information is also necessary. The disordered N-terminal part of p53 contains the transactivation and the proline-rich domains-which besides being abundant in proline residues-contains repetitive Pro-Ala motifs. NMR assignment of such repetitive, proline-rich regions is challenging due to the lack of amide protons in the 1HN-detected approaches, as well as due to the small chemical shift dispersion. In the present study we perform the full assignment of the p531-100 region by applying a combination of 1HN- and 1Hα-detected NMR experiments. We also show the increased information content when using real-time homo- and heteronuclear decoupled acquisition schemes. On the other hand, we highlight the presence of minor proline species, and using Pro-selective experiments we determine the corresponding cis or trans conformation. Secondary chemical shifts for (Cα-Cß) atoms indicate the disordered nature of this region, with expected helical tendency for the TAD1 region. As the role of the proline-rich domain is yet not well understood our results can contribute to further successful investigations.

Effects of structural changes on antibacterial activity and cytotoxicity due to proline substitutions in chimeric peptide HnMc.

Owing to the rapidly increasing emergence of multidrug-resistant pathogens, antimicrobial peptides (AMPs) are being explored as next-generation antibiotics. However, AMPs present in nature are highly toxic and exhibit low antibacterial activity. Simple modifications, such as amino acid substitution, can enhance antimicrobial activity and cell selectivity. Herein, we show that HnMc-W, substituted by the Phe1Trp analog of HnMc, a chimeric peptide, resulted in membranolytic antibacterial action and enhanced salt tolerance, whereas HnMc-WP1 with one Ser9Pro substitution resulted in a proline-kink helical structure that increased salt-tolerant antibacterial effects and reduced cytotoxicity. In addition, the HnMc-WP2 peptide, designed with a PXXP motif, had a flexible central hinge in its α-helical structure due to the introduction of two Pro and two Gln (X positions, by deletion of two Gln at positions 16 and 17) residues instead of Ser at position. HnMc-WP2 exhibited excellent antibacterial effects without cytotoxicity in vitro. Moreover, its potent antibacterial activity was demonstrated in a drug-resistant Pseudomonas aeruginosa-infected mouse model in vivo. Our findings provide valuable information for the design of peptides with high antibacterial activity and cell selectivity.

Proton Transport Across Collagen Fibrils and Scaffolds: The Role of Hydroxyproline.

Collagen is one of the most studied proteins due to its fundamental role in creating fibrillar structures and supporting tissues in our bodies. Accordingly, collagen is also one of the most used proteins for making tissue-engineered scaffolds for various types of tissues. To date, the high abundance of hydroxyproline (Hyp) within collagen is commonly ascribed to the structure and stability of collagen. Here, we hypothesize a new role for the presence of Hyp within collagen, which is to support proton transport (PT) across collagen fibrils. For this purpose, we explore here three different collagen-based hydrogels: the first is prepared by the self-assembly of natural collagen fibrils, and the second and third are based on covalently linking between collagen via either a self-coupling method or with an additional cross-linker. Following the formation of the hydrogel, we introduce here a two-step reaction, involving (1) attaching methanesulfonyl to the -OH group of Hyp, followed by (2) removing the methanesulfonyl, thus reverting Hyp to proline (Pro). We explore the PT efficiency at each step of the reaction using electrical measurements and show that adding the methanesulfonyl group vastly enhances PT, while reverting Hyp to Pro significantly reduces PT efficiency (compared with the initial point) with different efficiencies for the various collagen-based hydrogels. The role of Hyp in supporting the PT can assist in our understanding of the physiological roles of collagen. Furthermore, the capacity to modulate conductivity across collagen is very important to the use of collagen in regenerative medicine.

Efficient and accurate calculation of proline cis/trans isomerization free energies from Hamiltonian replica exchange molecular dynamics simulations.

Proline cis/trans isomerization plays an important role in many biological processes but occurs on time scales not accessible to brute-force molecular dynamics (MD) simulations. We have designed a new Hamiltonian replica exchange scheme, ω-bias potential replica exchange molecular dynamics (ωBP-REMD), to efficiently and accurately calculate proline cis/trans isomerization free energies. ωBP-REMD is applied to various proline-containing tripeptides and a biologically important proline residue in the N2-domain of the gene-3-protein of phage fd in the wildtype and mutant variants of the protein. Excellent cis/trans transition rates are obtained. Reweighting of the sampled probability distribution along the peptide bond dihedral angle allows construction of the corresponding free-energy profile and calculation of the cis/trans isomerization free energy with high statistical precision. Very good agreement with experimental data is obtained. ωBP-REMD outperforms standard umbrella sampling in terms of convergence and agreement with experiment and strongly reduces perturbation of the local structure near the proline residue.

Experimental NOE, Chemical Shift, and Proline Isomerization Data Provide Detailed Insights into Amelotin Oligomerization.

Amelotin is an intrinsically disordered protein (IDP) rich in Pro residues and is involved in hydroxyapatite mineralization. It rapidly oligomerizes under physiological conditions of pH and pressure but reverts to its monomeric IDP state at elevated pressure. We identified a 105-residue segment of the protein that becomes ordered upon oligomerization, and we used pressure-jump NMR spectroscopy to measure long-range NOE contacts that exist exclusively in the oligomeric NMR-invisible state. The kinetics of oligomerization and dissociation were probed at the residue-specific level, revealing that the oligomerization process is initiated in the C-terminal half of the segment. Using pressure-jump NMR, the degree of order in the oligomer at the sites of Pro residues was probed by monitoring changes in cis/trans equilibria relative to the IDP state after long-term equilibration under oligomerizing conditions. Whereas most Pro residues revert to trans in the oligomeric state, Pro-49 favors a cis configuration and three Pro residues retain an unchanged cis fraction, pointing to their local lack of order in the oligomeric state. NOE contacts and secondary 13C chemical shifts in the oligomeric state indicate the presence of an 11-residue α-helix, preceded by a small intramolecular antiparallel ß-sheet, with slower formation of long-range intermolecular interactions to N-terminal residues. Although none of the models generated by AlphaFold2 for the amelotin monomer was consistent with experimental data, subunits of a hexamer generated by AlphaFold-Multimer satisfied intramolecular NOE and chemical shift data and may provide a starting point for developing atomic models for the oligomeric state.