Closing the side-chain gap in protein loop modeling

The success of structure-based drug design relies on accurate protein modeling where one of the key issues is the modeling and refinement of loops. This study takes a critical look at modeled loops, determining the effect of re-sampling side-chains after the loop conformation has been generated. The results are evaluated in terms of backbone and side-chain conformations with respect to the native loop. While models can contain loops with high quality backbone conformations, the side-chain orientations could be poor, and therefore unsuitable for ligand docking and structure-based design. In this study, we report on the ability to model loop side-chains accurately using a variety of commercially available algorithms that include rotamer libraries, systematic torsion scans and knowledge-based methods. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Intrinsic energy landscapes of amino acid side-chains

Amino acid side-chain conformational properties influence the overall structural and dynamic properties of proteins and, therefore, their biological functions. In this study, quantum mechanical (QM) potential energy surfaces for the rotation of side-chain χ(1) and χ(2) torsions in dipeptides in the alphaR, beta, and alphaL backbone conformations were calculated. The QM energy surfaces provide a broad view of the intrinsic conformational properties of each amino acid side-chain. The extent to which intrinsic energetics dictates side-chain orientation was studied through comparisons of the QM energy surfaces with χ(1) and χ(2) free energy surfaces from probability distributions obtained from a survey of high resolution crystal structures. In general, the survey probability maxima are centered in minima of the QM surfaces as expected for sp(3) (or sp(2) for χ(2) of Asn, Phe, Trp, and Tyr) atom centers with strong variations between amino acids occurring in the energies of the minima indicating intrinsic differences in rotamer preferences. High correlations between the QM and survey data were found for hydrophobic side-chains except Met, suggesting minimal influence of the protein and solution environments on their conformational distributions. Conversely, low correlations for polar or charged side-chains indicate a dominant role of the environment in stabilizing conformations that are not intrinsically favored. Data also link the presence of off-rotamers in His and Trp to favorable interactions with the backbone. Results also suggest that the intrinsic energetics of the side-chains of Phe and Tyr may play important roles in protein folding and stability. Analyses on whether intrinsic side-chain energetics can influence backbone preference identified a strong correlation for residues in the alphaL backbone conformation. It is suggested that this correlation reflects the intrinsic instability of the alphaL backbone such that assumption of this backbone conformation is facilitated by intrinsically favorable side-chain conformations. Together our results offer a broad overview of the conformational properties of amino acid side-chains and the QM data may be used as target data for force field optimization.     

The Influence of the Amino Acid Side Chains on the Raman Optical Activity Spectra of Proteins

The Raman optical activity (ROA) spectra of proteins show distinct patterns arising from the secondary structure. It is generally believed that the spectral contributions of the side-chains largely cancel out because of their flexibility and the occurrence of many side-chains with different conformations. Yet, the influence of the side-chains on the ROA patterns assigned to different secondary structures is unknown. Here, the first systematic study of the influence of all amino acid side-chains on the ROA patterns is presented based on density functional theory (DFT) calculations of an extensive collection of peptide models that include many different side-chain and secondary structure conformations. It was shown that the contributions of the side-chains to a large extent average out with conformational flexibility. However, specific side-chain conformations can have significant contributions to the ROA patterns. It was also shown that α-helical structure is very sensitive to both the exact backbone conformation and the side-chain conformation. Side-chains with χ₁≈−60° generate ROA patterns alike those in experiment. Aromatic side-chains strongly influence the amide III ROA patterns. Because of the huge structural sensitivity of ROA, the spectral patterns of proteins arise from extensive conformational averaging of both the backbone and the side-chains. The averaging results in the fine spectral details and relative intensity differences observed in experimental spectra.     

Synthesis and conformational analysis of a cyclic peptide obtained via i to i+4 intramolecular side-chain to side-chain azide-alkyne 1,3-dipolar cycloaddition

Intramolecular side-chain to side-chain cyclization is an established approach to achieve stabilization of specific conformations and a recognized strategy to improve resistance toward proteolytic degradation. To this end, cyclizations, which are bioisosteric to the lactam-type side-chain to side-chain modification and do not require orthogonal protection schemes, are of great interest. Herein, we report the employment of Cu(I)-catalyzed 1,3-dipolar cycloaddition of side chains modified with azido and alkynyl functions and explore alternative synthetic routes to efficiently generate 1,4-disubstituted [1,2,3]triazolyl-containing cyclopeptides. The solid-phase assembly of the linear precursor including epsilon-azido norleucine and the propargylglycine (Pra) in positions i and i+4, respectively, was accomplished by either subjecting the resin-bound peptide to selective on-resin diazo transformation of a Lys into the Nle(epsilon-N3) or the incorporation of Fmoc-Nle(epsilon-N3)-OH during the stepwise build-up of the resin-bound peptide 1b. Solution-phase Cu(I)-catalyzed 1,3-dipolar cycloaddition converts the linear precursor Ac-Lys-Gly-Nle(epsilon-N3)-Ser-Ile-Gln-Pra-Leu-Arg-NH2 (2) into the 1,4-disubstituted [1,2,3]triazolyl-containing cyclopeptide [Ac-Lys-Gly-Xaa(&(1))-Ser-Ile-Gln-Yaa(&(2))-Leu-Arg-NH2][(&(1)(CH2)4-1,4-[1,2,3]triazolyl-CH2&(2))] (3). The conformational preferences of the model cyclopeptide 3 (III), which is derived from the sequence of a highly helical and potent i to i+4 side-chain to side-chain lactam-containing antagonist of parathyroid hormone-related peptide (PTHrP), are compared to the corresponding lactam analogue Ac[Lys(13)(&(1)),Asp(17)(&(2))]hPTHrP(11-19)NH2 (II). CD and NMR studies of 3 and II in water/hexafluoroacetone (HFA) (50:50, v/v) revealed a high prevalence of turn-helical structures involving in particular the cyclic regions of the molecule. Despite a slight difference of the backbone arrangement, the side-chains of Ser, Gln, and Ile located at the i+1 to i+3 of the ring-forming sequences share the same spatial orientation. Both cyclopeptides differ regarding the location of the turn-helical segment, which in II involves noncyclized residues while in 3 it overlaps with residues involved in the cyclic structure. Therefore, the synthetic accessibility and conformational similarity of i to i+4 side-chain to side-chain cyclopeptide containing the 1,4-disubstituted [1,2,3]triazolyl moiety to the lactam-type one may result in similar bioactivities.     

A new method for side-chain conformation prediction using a Hopfield network and reproduced rotamers

We present a new side-chain prediction method based on energy minimization using a Hopfield network, focusing on the buried residues of proteins. In this method, the network is composed of automata assigned to each rotamer to restrict side-chain conformational space. We reproduced a rotamer library that enabled us to more widely cover the space for side-chain conformations than those previously produced. The accuracy of the side-chain modeling was estimated by three standards: root mean square deviations (rmsds) between the modeled and the crystal structures, the percentages of correctly predicted side-chain torsion angles, and the percentages of correctly predicted hydrogen bonds. The average rmsd for buried side chains of 21 proteins was 1.10 Å. The value was almost always improved relative to the previous works. The percentage of side-chain X₁ angles for buried residues was 87.3%. By considering the hydrogen bond energy, the average percentage of correctly predicted hydrogen bonds rose from 33% without hydrogen bond energy to 52% with the bond energy. We applied this method to homology modeling, where the protein backbone used to predict side-chain conformations deviates from the correct conformation, and could predict side-chain conformations as correctly as those using the correct backbones. © 1996 by John Wiley & Sons, Inc.     

How side-chain hydrophilicity modulates morphology and charge transport in mixed conducting polymers

Organic mixed ionic-electronic conductors (OMIECs) are a developing class of organic electronic materials distinguished by their dual modes of conduction. The side-chains of OMIEC polymers are responsible for forming a percolating electrolyte phase that mediates doping and ionic conduction. Despite this critical role, design rules for OMIEC side-chains are still nascent and their effects on OMIEC morphology and charge transport have yet to be systematically studied. Here we perform the first dedicated coarse-grained molecular dynamics study of OMIECs where the side-chain identity and distribution are systematically varied using a random copolymer architecture. The simulations recapitulate the nonlinear progression of the morphology from an interfacially gated electrolyte when large fractions of hydrophobic side-chains are incorporated, to an electrolyte swelled morphology after crossing a threshold of approximately 40% polar side-chains. Kinetic Monte Carlo simulations were used to characterize the charge transport behaviors in these systems, revealing two interesting maxima in the mobility at 40% and 100% polar side-chain fractions, respectively. With respect to maximizing the charge mobility and conductivity, these simulations suggest that a uniform hydrophilic side-chain distribution is optimal and that there are few advantages to using mixed side-chains in a random copolymer architecture. These results also suggest several alternative side-chain engineering strategies for optimizing OMIEC performance.     

Fine grained sampling of residue characteristics using molecular dynamics simulation

In a fine-grained computational analysis of protein structure, we investigated the relationships between a residue's backbone conformations and its side-chain packing as well as conformations. To produce continuous distributions in high resolution, we ran molecular dynamics simulations over a set of protein folds (dynameome). In effect, the dynameome dataset samples not only the states well represented in the PDB but also the known states that are not well represented in the structural database. In our analysis, we characterized the mutual influence among the backbone ?,ψ angles with the first side-chain torsion angles (χ₁) and the volumes occupied by the side-chains. The dependencies of these relationships on side-chain environment and amino acids are further explored. We found that residue volumes exhibit dependency on backbone 2° structure conformation: side-chains pack more densely in extended β-sheet than in α-helical structures. As expected, residue volumes on the protein surface were larger than those in the interior. The first side-chain torsion angles are found to be dependent on the backbone conformations in agreement with previous studies, but the dynameome dataset provides higher resolution of rotamer preferences based on the backbone conformation. All three gauche^?, gauche⁺, and trans rotamers show different patterns of ?,ψ dependency, and variations in χ₁ value are skewed from their canonical values to relieve the steric strains. By demonstrating the utility of dynameomic modeling on the native state ensemble, this study reveals details of the interplay among backbone conformations, residue volumes and side-chain conformations.     

T-Analyst: a program for efficient analysis of protein conformational changes by torsion angles

T-Analyst is a user-friendly computer program for analyzing trajectories from molecular modeling. Instead of using Cartesian coordinates for protein conformational analysis, T-Analyst is based on internal bond-angle-torsion coordinates in which internal torsion angle movements, such as side-chain rotations, can be easily detected. The program computes entropy and automatically detects and corrects angle periodicity to produce accurate rotameric states of dihedrals. It also clusters multiple conformations and detects dihedral rotations that contribute hinge-like motions. Correlated motions between selected dihedrals can also be observed from the correlation map. T-Analyst focuses on showing changes in protein flexibility between different states and selecting representative protein conformations for molecular docking studies. The program is provided with instructions and full source code in Perl.     

An efficient, path-independent method for free-energy calculations

Classical free-energy methods depend on the definition of physical or nonphysical integration paths to calculate free-energy differences between states. This procedure can be problematic and computationally expensive when the states of interest do not overlap and are far apart in phase space. Here we introduce a novel method to calculate free-energy differences that is path-independent by transforming each end state into a reference state in which the vibrational entropy is the sole component of the total entropy, thus allowing direct computation of the relative free energy. We apply the method to calculate side-chain entropies of a beta-hairpin-forming peptide in a variety of backbone conformations, demonstrating its importance in determining structural propensities. We find that low-free-energy conformations achieve their stability through optimal trade off between enthalpic gains due to favorable interatomic interactions and entropic losses incurred by the same.     

NMR spectra of carotenoporphyrins. Computer-assisted conformational analysis

A computer-assisted method of conformational analysis for porphyrin molecules bearing flexible side-chains has been developed. The method utilizes the ring current-induced chemical shift changes of the side-chain protons which arise from the porphyrin macrocycle and any attached aryl rings. The treatment has been applied to a series of carotenoporphyrin molecules, which are important as models for a variety of photophysical processes in biological systems. Chemical shift data of sufficient accuracy for the conformational analysis were obtained from 500 MHz NMR experiments. The conformations of the carotenoporphyrins varied from extended ones with the carotenoid well away from the porphyrin ring to tightly folded species, depending on molecular constitution. The analytical method can be extended to other porphyrin-based systems.     

Gas-phase models of gamma turns: effect of side-chain/backbone interactions investigated by IR/UV spectroscopy and quantum chemistry

The conformations of laser-desorbed jet-cooled short peptide chains Ac-Phe-Xxx-NH2 (Xxx=Gly, Ala, Val, and Pro) have been investigated by IR/UV double resonance spectroscopy and density-functional-theory (DFT) quantum chemistry calculations. Singly gamma-folded backbone conformations (betaL-gamma) are systematically observed as the most stable conformers, showing that in these two-residue peptide chains, the local conformational preference of each residue is retained (betaL for Phe and gamma turn for Xxx). Besides, beta turns are also spontaneously formed but appear as minor conformers. The theoretical analysis suggests negligible inter-residue interactions of the main conformers, which enables us to consider these species as good models of gamma turns. In the case of valine, two similar types of gamma turns, differing by the strength of their hydrogen bond, have been found both experimentally and theoretically. This observation provides evidence for a strong flexibility of the peptide chain, whose minimum-energy structures are controlled by side-chain/backbone interactions. The qualitative conformational difference between the present species and the reversed sequence Ac-Xxx-Phe-NH2 is also discussed.     

Backbone and side-chain ordering in a small protein

We investigate the relation between backbone and side-chain ordering in a small protein. For this purpose, we have performed multicanonical simulations of the villin headpiece subdomain HP-36, an often used toy model in protein studies. Concepts of circular statistics are introduced to analyze side-chain fluctuations. In contrast to earlier studies on homopolypeptides [Wei et al., J. Phys. Chem. B 111, 4244 (2007)], we do not find collective effects leading to a separate transition. Rather, side-chain ordering is spread over a wide temperature range. Our results indicate a thermal hierarchy of ordering events, with side-chain ordering appearing at temperatures below the helix-coil transition but above the folding transition. We conjecture that this thermal hierarchy reflects an underlying temporal order, and that side-chain ordering facilitates the search for the correct backbone topology.     

Statistical mechanics of protein allostery: roles of backbone and side-chain structural fluctuations

A statistical mechanical model of allosteric transition of proteins is developed by extending the structure-based model of protein folding to cases that a protein has two different native conformations. Partition function is calculated exactly within the model and free-energy surfaces associated with allostery are derived. In this paper, the model of allosteric transition proposed in a previous paper [Proc. Natl. Acad. Sci. U.S.A 134, 7775 (2010)] is reformulated to describe both fluctuation in side-chain configurations and that in backbone structures in a balanced way. The model is applied to example proteins, Ras, calmodulin, and CheY: Ras undergoes the allosteric transition between guanosine diphosphate (GDP)-bound and guanosine triphosphate (GTP)-bound forms, and the model results show that the GDP-bound form is stabilized enough to prevent unnecessary signal transmission, but the conformation in the GTP-bound state bears large fluctuation in side-chain configurations, which may help to bind multiple target proteins for multiple pathways of signaling. The calculated results of calmodulin show the scenario of sequential ordering in Ca(2+) binding and the associated allosteric conformational change, which are realized though the sequential appearing of pre-existing structural fluctuations, i.e., fluctuations to show structures suitable to bind Ca(2+) before its binding. Here, the pre-existing fluctuations to accept the second and third Ca(2+) ions are dominated by the side-chain fluctuation. In CheY, the calculated side-chain fluctuation of Tyr106 is coordinated with the backbone structural change in the β4-α4 loop, which explains the pre-existing Y-T coupling process in this protein. Ability of the model to explain allosteric transitions of example proteins supports the view that the large entropic effects lower the free-energy barrier of allosteric transition.     

Fast side-chain losses in keV ion-induced dissociation of protonated peptides

In this paper, we compare ionization and dissociation of a series of singly and doubly protonated peptides, namely leucine enkephalin, bradykinin, LHRH and substance P as induced by collisions with keV H⁺, He⁺ and He²⁺. For all peptides under study, the fragmentation pattern depends strongly on the electronic structure of the projectile ions. Immonium ions, side-chains and their fragments dominate the spectrum whereas fragments due to peptide backbone cleavage are weak or even almost absent for He⁺. Here, resonant electron capture from the peptide is ruled out and only interaction channels accompanied by much higher excitation contribute. Cleavage of the side-chain linkage appears to be a process alternative to backbone fragmentation occurring after internal vibrational redistribution of excitation energy. Depending on the peptide, this process can lead to the loss of a side-chain cation (leucine enkephalin, LHRH) or a neutral side-chain (substance P).     

A Comparison of Cysteine-Conjugated Nitroxide Spin Labels for Pulse Dipolar EPR Spectroscopy

The structure-function and materials paradigms drive research on the understanding of structures and structural heterogeneity of molecules and solids from materials science to structural biology. Functional insights into complex architectures are often gained from a suite of complementary physicochemical methods. In the context of biomacromolecular structures, the use of pulse dipolar electron paramagnetic resonance spectroscopy (PDS) has become increasingly popular. The main interest in PDS is providing long-range nanometre distance distributions that allow for identifying macromolecular topologies, validating structural models and conformational transitions as well as docking of quaternary complexes. Most commonly, cysteines are introduced into protein structures by site-directed mutagenesis and modified site-specifically to a spin-labelled side-chain such as a stable nitroxide radical. In this contribution, we investigate labelling by four different commercial labelling agents that react through different sulfur-specific reactions. Further, the distance distributions obtained are between spin-bearing moieties and need to be related to the protein structure via modelling approaches. Here, we compare two different approaches to modelling these distributions for all four side-chains. The results indicate that there are significant differences in the optimum labelling procedure. All four spin-labels show differences in the ease of labelling and purification. Further challenges arise from the different tether lengths and rotamers of spin-labelled side-chains; both influence the modelling and translation into structures. Our comparison indicates that the spin-label with the shortest tether in the spin-labelled side-group, (bis-(2,2,5,5-Tetramethyl-3-imidazoline-1-oxyl-4-yl) disulfide, may be underappreciated and could increase the resolution of structural studies by PDS if labelling conditions are optimised accordingly.     

An unconventional ligand-binding mechanism of substrate-binding proteins: MD simulation and Markov state model analysis of BtuF

In conventional “Venus Flytrap” mechanism, substrate-binding proteins (SBPs) interconvert between the open and closed conformations. Upon ligand binding, SBPs form a tightly closed conformation with the ligand bound at the interface of two domains. This mechanism was later challenged by many type III SBPs, such as the vitamin B₁₂-binding protein BtuF, in which the apo- and holo-state proteins adopt very similar conformations. Here, we combined molecular dynamics simulation and Markov state model analysis to study the conformational dynamics of apo- and B₁₂-bound BtuF. The results indicate that the crystal structures represent the only stable conformation of BtuF. Meanwhile, both apo- and holo-BtuF undergo large-scale interdomain motions with little energy cost. B₁₂ binding casts little restraints on the interdomain motions, suggesting that ligand binding affinity is enhanced by the remaining conformational entropy of holo-BtuF. These results reveal a new paradigm of ligand recognition mechanism of SBPs. © 2019 Wiley Periodicals, Inc.     

Collision‐free poisson motion planning in ultra high‐dimensional molecular conformation spaces

The function of protein, RNA, and DNA is modulated by fast, dynamic exchanges between three‐dimensional conformations. Conformational sampling of biomolecules with exact and nullspace inverse kinematics, using rotatable bonds as revolute joints and noncovalent interactions as holonomic constraints, can accurately characterize these native ensembles. However, sampling biomolecules remains challenging owing to their ultra‐high dimensional configuration spaces, and the requirement to avoid (self‐) collisions, which results in low acceptance rates. Here, we present two novel mechanisms to overcome these limitations. First, we introduce temporary constraints between near‐colliding links. The resulting constraint varieties instantaneously redirect the search for collision‐free conformations, and couple motions between distant parts of the linkage. Second, we adapt a randomized Poisson‐disk motion planner, which prevents local oversampling and widens the search, to ultra‐high dimensions. Tests on several model systems show that the sampling acceptance rate can increase from 16% to 70%, and that the conformational coverage in loop modeling measured as average closeness to existing loop conformations doubled. Correlated protein motions identified with our algorithm agree with those from MD simulations. © 2018 Wiley Periodicals, Inc.     

On the molecular origin of cold denaturation of globular proteins

A polypeptide chain can adopt very different conformations, a fundamental distinguishing feature of which is the water accessible surface area, WASA, that is a measure of the layer around the polypeptide chain where the center of water molecules cannot physically enter, generating a solvent-excluded volume effect. The large WASA decrease associated with the folding of a globular protein leads to a large decrease in the solvent-excluded volume, and so to a large increase in the configurational/translational freedom of water molecules. The latter is a quantity that depends upon temperature. Simple calculations over the -30 to 150 °C temperature range, where liquid water can exist at 1 atm, show that such a gain decreases significantly on lowering the temperature below 0 °C, paralleling the decrease in liquid water density. There will be a temperature where the destabilizing contribution of the polypeptide chain conformational entropy exactly matches the stabilizing contribution of the water configurational/translational entropy, leading to cold denaturation.     

Recent developments in peptide ligation independent of amino acid side-chain functional group

Chemical synthesis of peptides and proteins has evolved into an indispensable tool for chemical biology. Peptide ligation is a straightforward technique for joining two short peptide fragments together via a native peptide bond to afford a larger natural peptide or protein. However, the junction sites are limited to several specific amino acids because most peptide ligations involve participation of the side-chain functional groups of the junction-site amino acids. To overcome such intrinsic limitations, “general” peptide ligations which do not rely on the side-chain functional group have been developed. This review summarized the recent developments in peptide ligations that are independent of side-chain functional group of ligation-junction-site amino acid.     

Importance of chirality and reduced flexibility of protein side chains: a study with square and tetrahedral lattice models

Side chains of amino acid residues are the determining factor that distinguishes proteins from other unstable chain polymers. In simple models they are often represented implicitly (e.g., by spin states) or simplified as one atom. Here we study side chain effects using two-dimensional square lattice and three-dimensional tetrahedral lattice models, with explicitly constructed side chains formed by two atoms of different chirality and flexibility. We distinguish effects due to chirality and effects due to side chain flexibilities, since residues in proteins are L residues, and their side chains adopt different rotameric states. For short chains, we enumerate exhaustively all possible conformations. For long chains, we sample effectively rare events such as compact conformations and obtain complete pictures of ensemble properties of conformations of these models at all compactness region. This is made possible by using sequential Monte Carlo techniques based on chain growth method. Our results show that both chirality and reduced side chain flexibility lower the folding entropy significantly for globally compact conformations, suggesting that they are important properties of residues to ensure fast folding and stable native structure. This corresponds well with our finding that natural amino acid residues have reduced effective flexibility, as evidenced by statistical analysis of rotamer libraries and side chain rotatable bonds. We further develop a method calculating the exact side chain entropy for a given backbone structure. We show that simple rotamer counting underestimates side chain entropy significantly for both extended and near maximally compact conformations. We find that side chain entropy does not always correlate well with main chain packing. With explicit side chains, extended backbones do not have the largest side chain entropy. Among compact backbones with maximum side chain entropy, helical structures emerge as the dominating configurations. Our results suggest that side chain entropy may be an important factor contributing to the formation of alpha helices for compact conformations.