Conformational Propensities Research Articles

Cyclic arylopeptoid oligomers: synthesis and conformational propensities of peptide-mimetic aromatic macrocycles.

Macrocyclic peptide-mimetic molecules are attracting renewed attention and have found widespread application in research fields ranging from biochemical science to materials science. Herein, we describe the synthesis and structural elucidation of cyclo[n]-p-arylopeptoids (classified into cyclic aromatic ε-amino acids) bearing various side chains, namely, C[n]pAP(Rn) (where n inside brackets denotes the number of main chain units and R inside parentheses represents side chains). We investigate the influence of n and R on the macrocyclization efficiency of linear p-arylopeptoid oligomers (n = 3, 4, 5) under high-dilution conditions with or without slow addition. The structures of the cyclo[4]-p-arylopeptoids (C[4]pAP(Rn)) and their conformational dynamics are disclosed on the basis of single-crystal X-ray analyses, viable-temperature (VT) 1H NMR studies, and density functional theory (DFT) calculations. We found two representative conformations (open and closed) of cyclo[4]-p-arylopeptoids (C[4]pAP(Rn)) in the solid state and whose preference in the solution state was most likely dependent on solvent polarity. We believe that this simple but dynamic macrocyclic peptide-mimetic molecular scaffold would be attractive for developing new functional molecular tools based on rational molecular design as well as molecular library screening strategies.

Structural insights into a thermostable variant of human carbonic anhydrase II.

Carbonic anhydrase is an enzyme of interest for many biotechnological developments including carbon sequestration. These applications often require harsh conditions, so there is a need for the development of thermostable variants. One of the most thermostable human carbonic anhydrase II (HCAIIts) variants was patented in 2006. Here, we report the ultra-high resolution crystal structure of HCAIIts. The structural changes seen are consistent with each of the six mutations involved acting largely independently and variously resulting in increased H-bonding, improved packing, and reduced side chain entropy loss on folding to yield the increased stability. We further suggest that for four of the mutations, improvements in backbone conformational energetics is also a contributor and that considerations of such conformational propensities of individual amino acids are often overlooked.

The influence of solvent on conformational properties of peptides with Aib residue\u2014a DFT study

The conformational propensities of the Aib residue on the example of two model peptides Ac-Aib-NHMe (1) and Ac-Aib-NMe2 (2), were studied by B3LYP and M06-2X functionals, in the gas phase and in the polar solvents. To verify the reliability of selected functionals, we also performed MP2 calculations for the tested molecules in vacuum. Polarizable continuum models (PCM and SMD) were used to estimate the solvent effect. Ramachandran maps were calculated to find all energy minima. Noncovalent intramolecular interactions due to hydrogen-bonds and dipole attractions between carbonyl groups are responsible for the relative stabilities of the conformers. In order to verify the theoretical results, the available conformations of similar X-ray structures from the Cambridge Crystallographic Data Center (CCDC) were analyzed. The results of the calculations show that both derivatives with the Aib residue in the gas phase prefer structures stabilized by intramolecular N–H⋯O hydrogen bonds, i.e., C5 and C7 conformations, while polar solvent promotes helical conformation with φ, ψ values equal to +/−60°, +/−40°. In addition, in the case of molecule 2, the helical conformation is the only one available in the polar environment. This result is fully consistent with the X-ray data.Graphical abstractEffect of solvent on the Ramachandran maps of the model peptides with Aib residue

Plant Polypeptide Hormone Systemin Prefers Polyproline II Conformation in Solution.

Systemin, an 18 amino-acid-signaling peptide, was the first plant polypeptide hormone to be discovered. Earlier structural studies involving NMR spectroscopy indicated a lack of definite structure in solution while circular dichroism spectroscopy suggested the presence of left-handed polyproline II (PPII) conformation. Here, we report the results of molecular dynamics simulations of the peptide in explicit solvent with two different force fields, namely, ff99SBildn and ff99IDPs, both of which showed a large propensity for PPII-like conformations in spite of showing differing features for other conformational characteristics. More remarkably, the conformations with predicted chemical shifts that agreed better with the NMR observations had a larger than average PPII content, especially for the ff99IDPs force field. An independent docking calculation of the molecule with the putative receptor SR160 also retained this conformational preference for PPII structure. The results suggest PPII to be an important class of conformation for systemin which may have a role in its bioactivity.

Protein folding, misfolding and aggregation: The importance of two-electron stabilizing interactions.

Proteins associated with neurodegenerative diseases are highly pleiomorphic and may adopt an all-α-helical fold in one environment, assemble into all-β-sheet or collapse into a coil in another, and rapidly polymerize in yet another one via divergent aggregation pathways that yield broad diversity of aggregates’ morphology. A thorough understanding of this behaviour may be necessary to develop a treatment for Alzheimer’s and related disorders. Unfortunately, our present comprehension of folding and misfolding is limited for want of a physicochemical theory of protein secondary and tertiary structure. Here we demonstrate that electronic configuration and hyperconjugation of the peptide amide bonds ought to be taken into account to advance such a theory. To capture the effect of polarization of peptide linkages on conformational and H-bonding propensity of the polypeptide backbone, we introduce a function of shielding tensors of the Cα atoms. Carrying no information about side chain-side chain interactions, this function nonetheless identifies basic features of the secondary and tertiary structure, establishes sequence correlates of the metamorphic and pH-driven equilibria, relates binding affinities and folding rate constants to secondary structure preferences, and manifests common patterns of backbone density distribution in amyloidogenic regions of Alzheimer’s amyloid β and tau, Parkinson’s α-synuclein and prions. Based on those findings, a split-intein like mechanism of molecular recognition is proposed to underlie dimerization of Aβ, tau, αS and PrPC, and divergent pathways for subsequent association of dimers are outlined; a related mechanism is proposed to underlie formation of PrPSc fibrils. The model does account for: (i) structural features of paranuclei, off-pathway oligomers, non-fibrillar aggregates and fibrils; (ii) effects of incubation conditions, point mutations, isoform lengths, small-molecule assembly modulators and chirality of solid-liquid interface on the rate and morphology of aggregation; (iii) fibril-surface catalysis of secondary nucleation; and (iv) self-propagation of infectious strains of mammalian prions.

Using Cooperatively Folded Peptides To Measure Interaction Energies and Conformational Propensities.

The rates and equilibria of the folding of biopolymers are determined by the conformational preferences of the subunits that make up the sequence of the biopolymer and by the interactions that are formed in the folded state in aqueous solution. Because of the centrality of these processes to life, quantifying conformational propensities and interaction strengths is vitally important to understanding biology. In this Account, we describe our use of peptide model systems that fold cooperatively yet are small enough to be chemically synthesized to measure such quantities. The necessary measurements are made by perturbing an interaction or conformation of interest by mutation and measuring the difference between the folding free energies of the wild type (in which the interaction or conformation is undisturbed) and the mutant model peptides (in which the interaction has been eliminated or the conformational propensities modified). With the proper controls and provided that the peptide model system in question folds via a two-state process, these folding free energy differences can be accurate measures of interaction strengths or conformational propensities. This method has the advantage of having high sensitivity and high dynamic range because the energies of interest are coupled to folding free energies, which can be measured with precisions on the order of a few tenths of a kilocalorie by well-established biophysical methods, like chaotrope or thermal denaturation studies monitored by fluorescence or circular dichroism. In addition, because the model peptides can be chemically synthesized, the full arsenal of natural and unnatural amino acids can be used to tune perturbations to be as drastic or subtle as desired. This feature is particularly noteworthy because it enables the use of analytical tools developed for physical organic chemistry, especially linear free energy relationships, to decompose interaction energies into their component parts to obtain a deeper understanding of the forces that drive interactions in biopolymers. We have used this approach, primarily with the WW domain derived from the human Pin1 protein as our model system, to assess hydrogen bond strengths (especially those formed by backbone amides), the dependence of hydrogen bond strengths on the environment in which they form, β-turn propensities of both natural sequences and small molecule β-turn mimics, and the energetics of carbohydrate-protein interactions. In each case, the combination of synthetic accessibility, the ease of measuring folding energies, and the robustness of the structure of the Pin1 WW domain to mutation enabled us to obtain incisive measurements of quantities that have been challenging to measure by other methods.

Metamorphic Proteins: Emergence of Dual Protein Folds from One Primary Sequence.

Every amino acid exhibits a different propensity for distinct structural conformations. Hence, decoding how the primary amino acid sequence undergoes the transition to a defined secondary structure and its final three-dimensional fold is presently considered predictable with reasonable certainty. However, protein sequences that defy the first principles of secondary structure prediction (they attain two different folds) have recently been discovered. Such proteins, aptly named metamorphic proteins, decrease the conformational constraint by increasing flexibility in the secondary structure and thereby result in efficient functionality. In this review, we discuss the major factors driving the conformational switch related both to protein sequence and to structure using illustrative examples. We discuss the concept of an evolutionary transition in sequence and structure, the functional impact of the tertiary fold, and the pressure of intrinsic and external factors that give rise to metamorphic proteins. We mainly focus on the major components of protein architecture, namely, the α-helix and β-sheet segments, which are involved in conformational switching within the same or highly similar sequences. These chameleonic sequences are widespread in both cytosolic and membrane proteins, and these folds are equally important for protein structure and function. We discuss the implications of metamorphic proteins and chameleonic peptide sequences in de novo peptide design.

Probing the Conformation-Dependent Preferential Binding of Ethanol to Cationic Glycylalanylglycine in Water/Ethanol by Vibrational and NMR Spectroscopy.

The conformational propensity of amino acid residues is determined by an intricate balance of peptide-solvent and solvent-solvent interactions. To explore how the systematic replacement of water by a cosolvent affects the solvation of both the amino acid backbone and side chains, we performed a combined vibrational spectroscopy and NMR study of cationic glycylalanylglycine (GAG) in different ethanol/water mixtures of between 0 and 42 mol percent ethanol. Classical model peptide N'-methylacetamide was used as a reference system to probe solvent-induced spectroscopic changes. The alanine residue of GAG in water is known to exhibit a very high propensity for polyproline II (pPII). Adding up to 30 mol % ethanol at room temperature leads only to minor changes in the Ramachandran distribution of alanine, which mostly changes within the individual conformational subspaces. A further increase in the ethanol fractions leads to a destabilization of pPII and a stabilization of β-strand conformations. At higher temperatures, different degrees of enthalpy-entropy compensations lead to a much stronger influence of ethanol on the peptide's conformational distribution. Ethanol-induced changes in chemical shifts and amide I wavenumbers strongly suggest that ethanol replaces water preferentially in the solvation shell of the polar C-terminal peptide group and of the alanine side chain, whereas the N-terminal group remains mostly hydrated. Furthermore, we found that ethanol interacts more strongly with the peptide if the latter adopts β-strand conformations. This leads to an unusual positive temperature coefficient for the chemical shift of the C-terminal amide proton. Our data suggests a picture in which GAG eventually accumulates at water-ethanol interfaces if the ethanol fractions exceed 0.3, which explains why the further addition of ethanol eventually causes self-aggregation and the subsequent formation of a hydrogel.

The influence of N-terminal acetylation on micelle-induced conformational changes and aggregation of α-Synuclein.

The biological function of α-Synuclein has been related to binding to lipids and membranes but these interactions can also mediate α-Synuclein aggregation, which is associated to Parkinson’s disease and other neuropathologies. In brain tissue α-Synuclein is constitutively N-acetylated, a modification that plays an important role in its conformational propensity, lipid and membrane binding, and aggregation propensity. We studied the interactions of the lipid-mimetic SDS with N-acetylated and non-acetylated α-Synuclein, as well as their early-onset Parkinson’s disease variants A30P, E46K and A53T. At low SDS/protein ratios α-Synuclein forms oligomeric complexes with SDS micelles with relatively low α-helical structure. These micellar oligomers can efficiently nucleate aggregation of monomeric α-Synuclein, with successive formation of oligomers, protofibrils, curly fibrils and mature amyloid fibrils. N-acetylation reduces considerably the rate of aggregation of WT α-Synuclein. However, in presence of any of the early-onset Parkinson’s disease mutations the protective effect of N-acetylation against micelle-induced aggregation becomes impaired. At higher SDS/protein ratios, N-acetylation favors another conformational transition, in which a second type of α-helix-rich, non-aggregating oligomers become stabilized. Once again, the Parkinson’s disease mutations disconnect the influence of N-acetylation in promoting this transition. These results suggest a cooperative link between the N-terminus and the region of the mutations that may be important for α-Synuclein function.

Tyrosine Kinase Activation and Conformational Flexibility: Lessons from Src-Family Tyrosine Kinases.

Protein kinases are enzymes that catalyze the covalent transfer of the γ-phosphate of an adenosine triphosphate (ATP) molecule onto a tyrosine, serine, threonine, or histidine residue in the substrate and thus send a chemical signal to networks of downstream proteins. They are important cellular signaling enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Unregulated protein kinase activity is often associated with a wide range of diseases, therefore making protein kinases major therapeutic targets. A prototypical system of central interest to understand the regulation of kinase activity is provided by tyrosine kinase c-Src, which belongs to the family of Src-related non-receptor tyrosine kinases (SFKs). Although the broad picture of autoinhibition via the regulatory domains and via the phosphorylation of the C-terminal tail is well characterized from a structural point of view, a detailed mechanistic understanding at the atomic-level is lacking. Advanced computational methods based on all-atom molecular dynamics (MD) simulations are employed to advance our understanding of tyrosine kinase activation. The computational studies suggest that the isolated kinase domain (KD) is energetically most favorable in the inactive conformation when the activation loop (A-loop) of the KD is not phosphorylated. The KD makes transient visits to a catalytically competent active-like conformation. The process of bimolecular trans-autophosphorylation of the A-loop eventually locks the KD in the active state. Activating point mutations may act by slightly increasing the population of the active-like conformation, enhancing the availability of the A-loop to be phosphorylated. The Src-homology 2 (SH2) and Src-homology 3 (SH3) regulatory domains, depending upon their configuration, either promote the inactive or the active state of the kinase domain. In addition to the roles played by the SH3, SH2, and KD, the Src-homology 4-Unique domain (SH4-U) region also serves as a key moderator of substrate specificity and kinase function. Thus, a fundamental understanding of the conformational propensity of the SH4-U region and how this affects the association to the membrane surface are likely to lead to the discovery of new intermediate states and alternate strategies for inhibition of kinase activity for drug discovery. The existence of a multitude of KD conformations poses a great challenge aimed at the design of specific inhibitors. One promising computational strategy to explore the conformational flexibility of the KD is to construct Markov state models from aggregated MD data.

Perfluoro-tert-butyl Homoserine Is a Helix-Promoting, Highly Fluorinated, NMR-Sensitive Aliphatic Amino Acid: Detection of the Estrogen Receptor·Coactivator Protein–Protein Interaction by 19F NMR

Highly fluorinated amino acids can stabilize proteins and complexes with proteins, via enhanced hydrophobicity, and provide novel methods for identification of specific molecular events in complex solutions, via selective detection by 19F NMR and the absence of native 19F signals in biological contexts. However, the potential applications of 19F NMR in probing biological processes are limited both by the strong propensities of most highly fluorinated amino acids for the extended conformation and by the relatively modest sensitivity of NMR spectroscopy, which typically constrains measurements to mid-micromolar concentrations. Herein, we demonstrate that perfluoro-tert-butyl homoserine exhibits a propensity for compact conformations, including α-helix and polyproline helix (PPII), that is similar to that of methionine. Perfluoro-tert-butyl homoserine has nine equivalent fluorines that do not couple to any other nuclei, resulting in a sharp singlet that can be sensitively detected rapidly at low micromolar concentrations. Perfluoro-tert-butyl homoserine was incorporated at sites of leucine residues within the α-helical LXXLL short linear motif of estrogen receptor (ER) coactivator peptides. A peptide containing perfluoro-tert-butyl homoserine at position i + 3 of the ER coactivator LXXLL motif exhibited a Kd of 2.2 μM for the estradiol-bound estrogen receptor, similar to that of the native ligand. 19F NMR spectroscopy demonstrated the sensitive detection (5 μM concentration, 128 scans) of binding of the peptide to the ER and of inhibition of protein-protein interaction by the native ligand or by the ER antagonist tamoxifen. These results suggest diverse potential applications of perfluoro-tert-butyl homoserine in probing protein function and protein-protein interfaces in complex solutions.

Exploring the Effects on the Conformational Propensity of Alanine in the Unblocked Tripeptide Glycyl-Analyl-Glycine in Water/Ethanol Mixtures

Short peptides have been shown to be model systems for obtaining the conformational propensities of individual amino acids. Recent studies on cationic GxG (x: different guest amino acid residues) peptides in water revealed the Ramachandran plots for 16 amino acid residues. Out of the natural amino acids, alanine in GAG has been shown to stand out owing to its unusually high propensity for polyproline II (pPII) conformations (0.72). Thus, its conformational distribution differs substantially from the classical ones of Ramachandran and Flory. Peptide-water interactions have been proposed as the major determinant of alanine's preference for pPII. To elucidate how peptide hydration affects conformational preferences of alanine we used FTIR, vibrational circular dichroism, polarized Raman and NMR spectroscopy to determine Ramachandran plots of GAG in aqueous solutions with 3, 12, and 42mol% ethanol. We observed a modest decrease of the pPII fraction from 0.72 to 0.65 and an increase in β-strand population from 0.15 to 0.20 from pure water to 12 mol% ethanol. Interestingly, the corresponding changes of the enthalpic and entropic differences between pPII and β -strand by far exceed changes in Gibbs energy at room temperature owing to enthalpy-entropy compensation effects. Our thermodynamic data reveal that the cosolvent stabilizes pPII enthalpically (−10 kJ/mol in water to −18 kJ/mol in 12 mol% ethanol) but destabilizes it even more entropically (−23 to −50 J/mol∗K). Comparisons of IR and NMR spectra of GAG and N- methylacetamide in the same ethanol/water mixtures reveal changes of peptide hydration which causes the observed destabilization of pPII of alanine and the concomitant increase of its conformational entropy.

Intrinsic α helix propensities compact hydrodynamic radii in intrinsically disordered proteins.

Proteins that lack tertiary stability under normal conditions, known as intrinsically disordered, exhibit a wide range of biological activities. Molecular descriptions for the biology of intrinsically disordered proteins (IDPs) consequently rely on disordered structural models, which in turn require experiments that assess the origins to structural features observed. For example, while hydrodynamic size is mostly insensitive to sequence composition in chemically denatured proteins, IDPs show strong sequence-specific effects in the hydrodynamic radius (Rh ) when measured under normal conditions. To investigate sequence-modulation of IDP Rh , disordered ensembles generated by a hard sphere collision model modified with a structure-based parameterization of the solution energetics were used to parse the contributions of net charge, main chain dihedral angle bias, and excluded volume on hydrodynamic size. Ensembles for polypeptides 10-35 residues in length were then used to establish power-law scaling relationships for comparison to experimental Rh from 26 IDPs. Results showed the expected outcomes of increased hydrodynamic size from increases in excluded volume and net charge, and compaction from chain-solvent interactions. Chain bias representing intrinsic preferences for α helix and polyproline II (PPII ), however, modulated Rh with intricate dependence on the simulated propensities. PPII propensities at levels expected in IDPs correlated with heightened Rh sensitivity to even weak α helix propensities, indicating bias for common (φ, ψ) are important determinants of hydrodynamic size. Moreover, data show that IDP Rh can be predicted from sequence with good accuracy from a small set of physicochemical properties, namely intrinsic conformational propensities and net charge. Proteins 2017; 85:296-311. © 2016 Wiley Periodicals, Inc.

PH-Responsive Mechanistic Switch Regulates the Formation of Dendritic and Fibrillar Nanostructures of a Functional Amyloid.

In contrast to pathological amyloids, functional amyloids are involved in crucial physiological functions. For instance, the melanosomal protein comprising a highly amyloidogenic polypeptide repeat domain assembles into amyloid fibrils that act as templates for melanin biosynthesis within acidic melanosomes. However, the mechanism-morphology-function relationship of functional amyloids is poorly understood. Here, we demonstrate that the repeat domain of the melanosomal protein exhibits two distinct types of aggregation pathways that display nanoscale polymorphism in acidic pH. In the pH range of 4.5-6, aggregation proceeds via a typical nucleation-dependent mechanism, resulting in the formation of highly ordered β-rich curvy thread-like fibrils. On the contrary, at pH < 4.5, aggregation occurs through a rapid nucleation-independent isodesmic polymerization process that yields dendritic aggregates having lower degree of internal packing. These dendritic nanostructures can be converted into more stable fibrils by switching the pH. The nanoscale polymorphism associated with the mechanistic switch is likely to be mediated by the altered conformational propensities and intermolecular interactions due to the protonation/deprotonation of critical glutamate residues. We propose that this striking shift in the mechanism that dictates the nanoscale morphology regulates the melanosomal maturation.

Stability and aggregation propensity do not fully account for the association of various germline variable domain gene segments with light chain amyloidosis.

Variable domain (VL) gene segments exhibit variable tendencies to be associated with light chain amyloidosis (AL). While few of them are very frequent in AL and give rise to most of the amyloidogenic light chains compiled at the sequence databases, other are rarely found among the AL cases. To analyze to which extent these tendencies depend on folding stability and aggregation propensity of the germline VL protein, we characterized VL proteins encoded by four AL-associated germline gene segments and one not associated to AL. We found that the AL-associated germline rVL proteins differ widely in conformational stability and propensity to in vitro amyloid aggregation. While in vitro the amyloid formation kinetics of these proteins correlate well with their folding stabilities, the folding stability does not clearly correlate with their germline's frequencies in AL. We conclude that the association of the VL genes segments to amyloidosis is not determined solely by the folding stability and aggregation propensity of the germline VL protein. Other factors, such as the frequencies of destabilizing mutations and susceptibility to proteolysis, must play a role in determining the light chain amyloidogenicity.

Representing environment-induced helix-coil transitions in a coarse grained peptide model

Coarse grained (CG) models are widely used in studying peptide self-assembly and nanostructure formation. One of the recurrent challenges in CG modeling is the problem of limited transferability, for example to different thermodynamic state points and system compositions. Understanding transferability is generally a prerequisite to knowing for which problems a model can be reliably used and predictive. For peptides, one crucial transferability question is whether a model reproduces the molecule's conformational response to a change in its molecular environment. This is of particular importance since CG peptide models often have to resort to auxiliary interactions that aid secondary structure formation. Such interactions take care of properties of the real system that are per se lost in the coarse graining process such as dihedral-angle correlations along the backbone or backbone hydrogen bonding. These auxiliary interactions may then easily overstabilize certain conformational propensities and therefore destroy the ability of the model to respond to stimuli and environment changes, i.e. they impede transferability. In the present paper we have investigated a short peptide with amphiphilic EALA repeats which undergoes conformational transitions between a disordered and a helical state upon a change in pH value or due to the presence of a soft apolar/polar interface. We designed a base CG peptide model that does not carry a specific (backbone) bias towards a secondary structure. This base model was combined with two typical approaches of ensuring secondary structure formation, namely a C α -C α -C α -C α pseudodihedral angle potential or a virtual site interaction that mimics hydrogen bonding. We have investigated the ability of the two resulting CG models to represent the environment-induced conformational changes in the helix-coil equilibrium of EALA. We show that with both approaches a CG peptide model can be obtained that is environment-transferable and that correctly represents the peptide's conformational response to different stimuli compared to atomistic reference simulations. The two types of auxiliary interactions lead to different kinetic behavior as well as to different structural properties for fully formed helices and folding intermediates, and we discuss the advantages and disadvantages of these approaches.

Molecular Determinants of Tubulin's C-Terminal Tail Conformational Ensemble.

Tubulin is important for a wide variety of cellular processes including cell division, ciliogenesis, and intracellular trafficking. To perform these diverse functions, tubulin is regulated by post-translational modifications (PTM), primarily at the C-terminal tails of both the α- and β-tubulin heterodimer subunits. The tubulin C-terminal tails are disordered segments that are predicted to extend from the ordered tubulin body and may regulate both intrinsic properties of microtubules and the binding of microtubule associated proteins (MAP). It is not understood how either interactions with the ordered tubulin body or PTM affect tubulin's C-terminal tails. To probe these questions, we developed a method to isotopically label tubulin for C-terminal tail structural studies by NMR. The conformational changes of the tubulin tails as a result of both proximity to the ordered tubulin body and modification by mono- and polyglycine PTM were determined. The C-terminal tails of the tubulin dimer are fully disordered and, in contrast with prior simulation predictions, exhibit a propensity for β-sheet conformations. The C-terminal tails display significant chemical shift differences as compared to isolated peptides of the same sequence, indicating that the tubulin C-terminal tails interact with the ordered tubulin body. Although mono- and polyglycylation affect the chemical shift of adjacent residues, the conformation of the C-terminal tail appears insensitive to the length of polyglycine chains. Our studies provide important insights into how the essential disordered domains of tubulin function.

Further along the Road Less Traveled: AMBER ff15ipq, an Original Protein Force Field Built on a Self-Consistent Physical Model

We present the AMBER ff15ipq force field for proteins, the second-generation force field developed using the Implicitly Polarized Q (IPolQ) scheme for deriving implicitly polarized atomic charges in the presence of explicit solvent. The ff15ipq force field is a complete rederivation including more than 300 unique atomic charges, 900 unique torsion terms, 60 new angle parameters, and new atomic radii for polar hydrogens. The atomic charges were derived in the context of the SPC/Eb water model, which yields more-accurate rotational diffusion of proteins and enables direct calculation of nuclear magnetic resonance (NMR) relaxation parameters from molecular dynamics simulations. The atomic radii improve the accuracy of modeling salt bridge interactions relative to contemporary fixed-charge force fields, rectifying a limitation of ff14ipq that resulted from its use of pair-specific Lennard-Jones radii. In addition, ff15ipq reproduces penta-alanine J-coupling constants exceptionally well, gives reasonable agreement with NMR relaxation rates, and maintains the expected conformational propensities of structured proteins/peptides, as well as disordered peptides—all on the microsecond (μs) time scale, which is a critical regime for drug design applications. These encouraging results demonstrate the power and robustness of our automated methods for deriving new force fields. All parameters described here and the mdgx program used to fit them are included in the AmberTools16 distribution.

Structural properties of the linkers connecting the N- and C- terminal domains in the MocR bacterial transcriptional regulators

Peptide inter-domain linkers are peptide segments covalently linking two adjacent domains within a protein. Linkers play a variety of structural and functional roles in naturally occurring proteins. In this work we analyze the sequence properties of the predicted linker regions of the bacterial transcriptional regulators belonging to the recently discovered MocR subfamily of the GntR regulators. Analyses were carried out on the MocR sequences taken from the phyla Actinobacteria, Firmicutes, Alpha-, Beta- and Gammaproteobacteria. The results suggest that MocR linkers display phylum-specific characteristics and unique features different from those already described for other classes of inter-domain linkers. They show an average length significantly higher: 31.8 ± 14.3 residues reaching a maximum of about 150 residues. Compositional propensities displayed general and phylum-specific trends. Pro is dominating in all linkers. Dyad propensity analysis indicate Pro–Pro as the most frequent amino acid pair in all linkers. Physicochemical properties of the linker regions were assessed using amino acid indices relative to different features: in general, MocR linkers are flexible, hydrophilic and display propensity for β-turn or coil conformations. Linker sequences are hypervariable: only similarities between MocR linkers from organisms related at the level of species or genus could be found with sequence searches. The results shed light on the properties of the linker regions of the new MocR subfamily of bacterial regulators and may provide knowledge-based rules for designing artificial linkers with desired properties.

Structural propensities of kinase family proteins from a Potts model of residue co-variation.

Understanding the conformational propensities of proteins is key to solving many problems in structural biology and biophysics. The co-variation of pairs of mutations contained in multiple sequence alignments of protein families can be used to build a Potts Hamiltonian model of the sequence patterns which accurately predicts structural contacts. This observation paves the way to develop deeper connections between evolutionary fitness landscapes of entire protein families and the corresponding free energy landscapes which determine the conformational propensities of individual proteins. Using statistical energies determined from the Potts model and an alignment of 2896 PDB structures, we predict the propensity for particular kinase family proteins to assume a "DFG-out" conformation implicated in the susceptibility of some kinases to type-II inhibitors, and validate the predictions by comparison with the observed structural propensities of the corresponding proteins and experimental binding affinity data. We decompose the statistical energies to investigate which interactions contribute the most to the conformational preference for particular sequences and the corresponding proteins. We find that interactions involving the activation loop and the C-helix and HRD motif are primarily responsible for stabilizing the DFG-in state. This work illustrates how structural free energy landscapes and fitness landscapes of proteins can be used in an integrated way, and in the context of kinase family proteins, can potentially impact therapeutic design strategies.