Side-chain Packing Research Articles

Peptide flat-rod formation by precise arrangement among enantiomeric hydrophobic helices

Precise control of molecular arrangement is essential for functional molecular assemblies. A linear (I-shaped) amphiphilic block copolypeptide, polysarcosine-b-(l-Leu-Aib)6 (I-SL12), which has a hydrophilic polysarcosine (PSar) chain and a hydrophobic helical block, was reported to self-assemble into nanotubes by regular packing of the Leu side chains. Here, we have synthesized a T-shaped amphiphilic block copolypeptide, (l-Leu-Aib)3-AzF(PSar)-Aib-(l-Leu-Aib)2 (T-SL12), to investigate the effect of molecular geometry on the morphology of molecular assemblies. Unlike conventional I-SL12, T-SL12 self-assembles into helical nanotubes. A mixture of T-SL12 (a right-handed helix) and polysarcosine-b-(d-Leu-Aib)6 (I-SdL12, a left-handed helix) formed flat rod-shaped structures, while the mixture of T-SL12 and I-SL12 (both right-handed) forms nanotubes with an 80-nm diameter. This result indicates that stereo-complexes was formed between T-SL12 and I-SdL12. Peptidic flat-rod were obtained at ratios of T-SL12 and I-SdL12 from 1:1 to 1:3 (wt/wt), although their width (ca. 12 nm) and length (50–200 nm) did not change with stoichiometry. The thickness (6 nm) of the flat rod was measured by AFM. From these dimensions, we propose that the minor axis of peptidic flat-rod is composed of two stereo-complexed heterodimers of T-SL12 and I-SdL12 by orienting the I-SdL12s facing each other, and that this four-peptide unit is repeated side-by-side along the long axis.

Heterochirality-Mediated Cross-Strand Nested Hydrophobic Interaction Effects Manifested in Surface-Bound Peptide Assembly Structures.

Amino acid chirality has been envisioned as an important strategy to regulate structure and function of peptide self-assembled architectures. However, the molecular mechanism of chirality effects in peptide assemblies remains largely elusive. Here, the assembly structures of l-peptide polyphenylalanine F10 (FFFFFFFFFF) and block heterochiral peptide F5f5 (FFFFFfffff) composed of two FFFFF repeat blocks with opposite chirality were characterized at the single-molecule level by using scanning tunneling microscopy. Each peptide formed two distinctively different assembly structures on the HOPG surface, in which peptide chains took parallel and antiparallel β-sheet conformations, respectively. The molecular-level observations revealed that the staggered arrangement of cross-strand side chains achieved in the antiparallel β-sheet structure of the block heterochiral peptide facilitated intimate packing of side chains and maximized inter-residue van der Waals interactions, which led to more residues participating in assembly and greatly stabilized the β-sheet structure of the surface-bound peptide assembly, but such cross-strand nested interactions were not accessible in the heterochiral parallel β-sheet structure and the enantiomerically pure assembly structures. This work could contribute to the molecular insights of stereochemical interactions in peptide assemblies and feasibility of extending this block heterochirality pattern to other peptides with various lengths and amino acid compositions for structural regulations.

Comparison of side-chain dispersion in protein structures determined by cryo-EM and X-ray crystallography.

An evaluation of systematic differences in local structure and conformation in the interior of protein tertiary structures determined by crystallography and by cryo-electron microscopy (cryo-EM) is reported. The expectation is that any consistent differences between the derived atomic models could provide insights into variations in side-chain packing that result from differences in specimens prepared for analysis between these two methods. By computing an atomic packing score, which provides a quantitative measure of clustering of side-chain atoms in the core of the tertiary structures, it is found that, in general, for structures determined by cryo-EM, side chains are more dispersed than in structures determined by X-ray crystallography over a similar resolution range. This trend is also observed in the packing comparison at subunit interfaces. Similar trends were observed in the packing comparison at the core of tertiary structures of the same proteins determined by both X-ray and cryo-EM methods. It is proposed here that the reduced dispersion of side chains in protein crystals could be due to some level of dehydration in 3D crystals prepared for X-ray crystallography and also because the higher rate of freezing of protein samples for cryo-EM may enable preservation of a more native conformation.

The amyloid concentric β-barrel hypothesis: Models of synuclein oligomers, annular protofibrils, lipoproteins, and transmembrane channels.

Amyloid beta (Aβ of Alzheimer's disease) and α-synuclein (α-Syn of Parkinson's disease) form large fibrils. Evidence is increasing however that much smaller oligomers are more toxic and that these oligomers can form transmembrane ion channels. We have proposed previously that Aβ42 oligomers, annular protofibrils, and ion channels adopt concentric β-barrel molecular structures. Here we extend that hypothesis to the superfamily of α, β, and γ-synucleins. Our models of numerous synuclein oligomers, annular protofibrils, tubular protofibrils, lipoproteins, and ion channels were developed to be consistent with sizes, shapes, molecular weights, and secondary structures of assemblies as determined by electron microscopy and other studies. The models have the following features: (1) all subunits have identical structures and interactions; (2) they are consistent with conventional β-barrel theory; (3) the distance between walls of adjacent β-barrels is between 0.6 and 1.2 nm; (4) hydrogen bonds, salt bridges, interactions among aromatic side-chains, burial and tight packing of hydrophobic side-chains, and aqueous solvent exposure of hydrophilic side-chains are relatively optimal; and (5) residues that are identical among distantly related homologous proteins cluster in the interior of most oligomers whereas residues that are hypervariable are exposed on protein surfaces. Atomic scale models of some assemblies were developed.

AlphaFold 2: Why It Works and Its Implications for Understanding the Relationships of Protein Sequence, Structure, and Function.

AlphaFold 2 (AF2) was the star of CASP14, the last biannual structure prediction experiment. Using novel deep learning, AF2 predicted the structures of many difficult protein targets at or near experimental resolution. Here, we present our perspective of why AF2 works and show that it is a very sophisticated fold recognition algorithm that exploits the completeness of the library of single domain PDB structures. It has also learned local side chain packing rearrangements that enable it to refine proteins to high resolution. The benefits and limitations of its ability to predict the structures of many more proteins at or close to atomic detail are discussed.

Temperature induced helical contraction and expansion in branched polycarbodiimides and their solvent vapor sensing properties

Helical polycarbodiimides with either naphthalene or anthracene side chains display chiroptical switching properties in solution and the solid state. This interesting chiroptical switching occurs due to the expansion and contraction of the helical backbone and the simultaneous reorientation of the aromatic side chains. Herein, we modified the aliphatic side-chain structure to induce the switching behavior in simple phenyl substituted polycarbodiimides. The polymerization of 1-hexynyl-5-phenylcarbodiimide introduced a terminal alkyne appended to the polymer backbone. By using the CLICK reaction, we were able to attach steric bulky branched alkyl groups to poly(N-phenyl)-N′-(n-branched alkyl) carbodiimide (P-1B). To increase the solubility and modulate the aliphatic sidechain packing of the polymer, 1-hexynyl-5-phenylcarbodiimide was copolymerized with phenyl octadecyl carbodiimide in different mole ratios. IR spectra and Specific Optical Rotation (SOR) studies showed these simple phenyl polymers show temperature driven solution state switching between two and sometimes three conformational states akin to what was only previously observed in polymers with polyaromatic side chains. In addition, one of the copolymers, P-2B, synthesized from the 50/50 mole feed ratio of the branched monomer and phenyl octadecyl carbodiimide showed solid-state sensing of various volatile organic compounds (VOCs) and solvent vapor.

Dimerization of Cadherin-11 involves multi-site coupled unfolding and strand swapping

Cadherin extracellular domain 1 (EC1) mediates homophilic dimerization in adherens junctions. Conserved Trp2 and Trp4 residues in type II cadherins anchor the EC1 A strand intermolecularly in strand-swapped dimers. Herein, NMR spectroscopy is used to elucidate the roles of Trp2 and Trp4 in Cadherin-11 dimerization. The monomeric state, with the A strand and Trp side chains packed intramolecularly, is in equilibrium with sparsely populated partially and fully A-strand-exposed states, in which Trp2 (and Trp4, respectively) side-chain packing is disrupted. Exchange kinetics between the major state and the partially (fully) A-strand-exposed state is slow-intermediate (intermediate-fast). A separate very fast process exchanges ordered and random-coil BC-loop conformations with populations dependent on A-strand exposure and dimerization status. In addition, very slow processes connect the folded A-strand-exposed conformation to partially unfolded states, which may represent additional domain-swapping intermediates. The dimerization mechanism of type II cadherins is revealed as coupled folding and strand swapping.

Optimal Protein Sequence Design Mitigates Mechanical Failure in Silk β-Sheet Nanocrystals.

The excellent mechanical strength and toughness of spider silk are well characterized experimentally and understood atomistically using computational simulations. However, little attention has been focused on understanding whether the amino acid sequence of β-sheet nanocrystals, which is the key to rendering strength to silk fiber, is optimally chosen to mitigate molecular-scale failure mechanisms. To investigate this, we modeled β-sheet nanocrystals of various representative small/polar/hydrophobic amino acid repeats for determining the sequence motif having superior nanomechanical tensile strength and toughness. The constant velocity pulling of the central β-strand in the nanocrystal, using steered molecular dynamics, showed that homopolymers of small amino acid (alanine/alanine-glycine) sequence motifs, occurring in natural silk fibroin, have better nanomechanical properties than other modeled structures. Further, we analyzed the hydrogen bond (HB) and β-strand pull dynamics of modeled nanocrystals to understand the variation in their rupture mechanisms and explore sequence-dependent mitigating factors contributing to their superior mechanical properties. Surprisingly, the enhanced side-chain interactions in homopoly-polar/hydrophobic amino acid models are unable to augment backbone HB cooperativity to increase mechanical strength. Our analyses suggest that nanocrystals of pristine silk sequences most likely achieve superior mechanical strength by optimizing side-chain interaction, packing, and main-chain HB interactions. Thus, this study suggests that the nanocrystal β-sheet sequence plays a crucial role in determining the nanomechanical properties of silk, and the evolutionary process has optimized it in natural silk. This study provides insight into the molecular design principle of silk with implications in the genetically modified artificial synthesis of silk-like biomaterials.

Computational IR Spectroscopy of Insulin Dimer Structure and Conformational Heterogeneity.

We have investigated the structure and conformational dynamics of insulin dimer using a Markov state model (MSM) built from extensive unbiased atomistic molecular dynamics simulations and performed infrared spectral simulations of the insulin MSM to describe how structural variation within the dimer can be experimentally resolved. Our model reveals two significant conformations to the dimer: a dominant native state consistent with other experimental structures of the dimer and a twisted state with a structure that appears to reflect a ∼55° clockwise rotation of the native dimer interface. The twisted state primarily influences the contacts involving the C-terminus of insulin's B chain, shifting the registry of its intermolecular hydrogen bonds and reorganizing its side-chain packing. The MSM kinetics predict that these configurations exchange on a 14 μs time scale, largely passing through two Markov states with a solvated dimer interface. Computational amide I spectroscopy of site-specifically 13C18O labeled amides indicates that the native and twisted conformation can be distinguished through a series of single and dual labels involving the B24F, B25F, and B26Y residues. Additional structural heterogeneity and disorder is observed within the native and twisted states, and amide I spectroscopy can also be used to gain insight into this variation. This study will provide important interpretive tools for IR spectroscopic investigations of insulin structure and transient IR kinetics experiments studying the conformational dynamics of insulin dimer.

Enzyme catalysis prior to aromatic residues: Reverse engineering of a dephospho-CoA kinase.

The wide variety of protein structures and functions results from the diverse properties of the 20 canonical amino acids. The generally accepted hypothesis is that early protein evolution was associated with enrichment of a primordial alphabet, thereby enabling increased protein catalytic efficiencies and functional diversification. Aromatic amino acids were likely among the last additions to genetic code. The main objective of this study was to test whether enzyme catalysis can occur without the aromatic residues (aromatics) by studying the structure and function of dephospho-CoA kinase (DPCK) following aromatic residue depletion. We designed two variants of a putative DPCK from Aquifex aeolicus by substituting (a) Tyr, Phe and Trp or (b) all aromatics (including His). Their structural characterization indicates that substituting the aromatics does not markedly alter their secondary structures but does significantly loosen their side chain packing and increase their sizes. Both variants still possess ATPase activity, although with 150-300 times lower efficiency in comparison with the wild-type phosphotransferase activity. The transfer of the phosphate group to the dephospho-CoA substrate becomes heavily uncoupled and only the His-containing variant is still able to perform the phosphotransferase reaction. These data support the hypothesis that proteins in the early stages of life could support catalytic activities, albeit with low efficiencies. An observed significant contraction upon ligand binding is likely important for appropriate organization of the active site. Formation of firm hydrophobic cores, which enable the assembly of stably structured active sites, is suggested to provide a selective advantage for adding the aromatic residues.

Protein sequence optimization with a pairwise decomposable penalty for buried unsatisfied hydrogen bonds

In aqueous solution, polar groups make hydrogen bonds with water, and hence burial of such groups in the interior of a protein is unfavorable unless the loss of hydrogen bonds with water is compensated by formation of new ones with other protein groups. For this reason, buried "unsatisfied" polar groups making no hydrogen bonds are very rare in proteins. Efficiently representing the energetic cost of unsatisfied hydrogen bonds with a pairwise-decomposable energy term during protein design is challenging since whether or not a group is satisfied depends on all of its neighbors. Here we describe a method for assigning a pairwise-decomposable energy to sidechain rotamers such that following combinatorial sidechain packing, buried unsaturated polar atoms are penalized. The penalty can be any quadratic function of the number of unsatisfied polar groups, and can be computed very rapidly. We show that inclusion of this term in Rosetta sidechain packing calculations substantially reduces the number of buried unsatisfied polar groups.

Altering CLC stoichiometry by reducing non-polar side-chains at the dimerization interface

CLC-ec1 is a Cl−/H+ antiporter that forms stable homodimers in lipid bilayers, with a free energy of −10.9 kcal/mol in 2:1 POPE/POPG lipid bilayers. The dimerization interface is formed by four transmembrane helices: H, I, P and Q, that are lined by non-polar side-chains that come in close contact, yet it is unclear as to whether their interactions drive dimerization. To investigate whether non-polar side-chains are required for dimer assembly, we designed a series of constructs where side-chain packing in the dimer state is significantly reduced by making 4–5 alanine substitutions along each helix (H-ala, I-ala, P-ala, Q-ala). All constructs are functional and three purify as stable dimers in detergent micelles despite the removal of significant side-chain interactions. On the other hand, H-ala shows the unique behavior of purifying as a mixture of monomers and dimers, followed by a rapid and complete conversion to monomers. In lipid bilayers, all four constructs are monomeric as examined by single-molecule photobleaching analysis. Further study of the H-helix shows that the single mutation L194A is sufficient to yield monomeric CLC-ec1 in detergent micelles and lipid bilayers. X-ray crystal structures of L194A reveal the protein re-assembles to form dimers, with a structure that is identical to wild-type. Altogether, these results demonstrate that non-polar membrane embedded side-chains play an important role in defining dimer stability, but the stoichiometry is highly contextual to the solvent environment. Furthermore, we discovered that L194 is a molecular hot-spot for defining dimerization of CLC-ec1.

OPUS-Rota3: Improving Protein Side-Chain Modeling by Deep Neural Networks and Ensemble Methods.

Side-chain modeling is critical for protein structure prediction since the uniqueness of the protein structure is largely determined by its side-chain packing conformation. In this paper, differing from most approaches that rely on rotamer library sampling, we first propose a novel side-chain rotamer prediction method based on deep neural networks, named OPUS-RotaNN. Then, on the basis of our previous work OPUS-Rota2, we propose an open-source side-chain modeling framework, OPUS-Rota3, which integrates the results of different methods into its rotamer library as the sampling candidates. By including OPUS-RotaNN into OPUS-Rota3, we conduct our experiments on three native backbone test sets and one non-native backbone test set. On the native backbone test set, CAMEO-Hard61 for example, OPUS-Rota3 successfully predicts 51.14% of all side-chain dihedral angles with a tolerance criterion of 20° and outperforms OSCAR-star (50.87%), SCWRL4 (50.40%), and FASPR (49.85%). On the non-native backbone test set DB379-ITASSER, the accuracy of OPUS-Rota3 is 52.49%, better than OSCAR-star (48.95%), FASPR (48.69%), and SCWRL4 (48.29%). All the source codes including the training codes and the data we used are available at https://github.com/thuxugang/opus_rota3.

Microstructure of bottlebrush poly(n-alkyl methacrylate)s beyond side chain packing

Based on results of small and wide angle X-ray scattering measurements, we report a new microstructural feature in bottlebrush poly (n-alkyl methacrylate)s larger than side chain packing. Temperature dependent X-ray experiments confirmed that such large scale aggregation structure is originated from the crystallization of the side chains. Two-dimensional X-ray scattering patterns revealed that the structure is due to the formation of disordered regions locate on average two layers of well crystallized side chains apart. Crystallization of the side chains rejected backbone entanglements to the surrounding area forming regions of high entanglement density so that restricting the crystallization of side chains. As such a density contrast is created between the crystallized regions and the disordered ones which results in strong small angle X-ray scattering signals.

FASPR: an open-source tool for fast and accurate protein side-chain packing.

Protein structure and function are essentially determined by how the side-chain atoms interact with each other. Thus, accurate protein side-chain packing (PSCP) is a critical step toward protein structure prediction and protein design. Despite the importance of the problem, however, the accuracy and speed of current PSCP programs are still not satisfactory. We present FASPR for fast and accurate PSCP by using an optimized scoring function in combination with a deterministic searching algorithm. The performance of FASPR was compared with four state-of-the-art PSCP methods (CISRR, RASP, SCATD and SCWRL4) on both native and non-native protein backbones. For the assessment on native backbones, FASPR achieved a good performance by correctly predicting 69.1% of all the side-chain dihedral angles using a stringent tolerance criterion of 20°, compared favorably with SCWRL4, CISRR, RASP and SCATD which successfully predicted 68.8%, 68.6%, 67.8% and 61.7%, respectively. Additionally, FASPR achieved the highest speed for packing the 379 test protein structures in only 34.3 s, which was significantly faster than the control methods. For the assessment on non-native backbones, FASPR showed an equivalent or better performance on I-TASSER predicted backbones and the backbones perturbed from experimental structures. Detailed analyses showed that the major advantage of FASPR lies in the optimal combination of the dead-end elimination and tree decomposition with a well optimized scoring function, which makes FASPR of practical use for both protein structure modeling and protein design studies. The web server, source code and datasets are freely available at https://zhanglab.ccmb.med.umich.edu/FASPR and https://github.com/tommyhuangthu/FASPR. Supplementary data are available at Bioinformatics online.

Polymorphic A\u03b242 fibrils adopt similar secondary structure but differ in cross-strand side chain stacking interactions within the same \u03b2-sheet

Formation of polymorphic amyloid fibrils is a common feature in neurodegenerative diseases involving protein aggregation. In Alzheimer’s disease, different fibril structures may be associated with different clinical sub-types. Structural basis of fibril polymorphism is thus important for understanding the role of amyloid fibrils in the pathogenesis and progression of these diseases. Here we studied two types of Aβ42 fibrils prepared under quiescent and agitated conditions. Quiescent Aβ42 fibrils adopt a long and twisted morphology, while agitated fibrils are short and straight, forming large bundles via lateral association. EPR studies of these two types of Aβ42 fibrils show that the secondary structure is similar in both fibril polymorphs. At the same time, agitated Aβ42 fibrils show stronger interactions between spin labels across the full range of the Aβ42 sequence, suggesting a more tightly packed structure. Our data suggest that cross-strand side chain packing interactions within the same β-sheet may play a critical role in the formation of polymorphic fibrils.

Side chain packing states of chitosan‐based supramolecular derivatives containing long alkyl side chains

Side chain packing states of chitosan‐based supramolecular derivatives containing long alkyl side chains

PPerturb: A Server for Predicting Long-Distance Energetic Couplings and Mutation-Induced Stability Changes in Proteins via Perturbations.

The strength of intraprotein interactions or contact network is one of the dominant factors determining the thermodynamic stabilities of proteins. The nature and the extent of connectivity of this network also play a role in allosteric signal propagation characteristics upon ligand binding to a protein domain. Here, we develop a server for rapid quantification of the strength of an interaction network by employing an experimentally consistent perturbation approach previously validated against a large data set of 375 mutations in 19 different proteins. The web server can be employed to predict the extent of destabilization of proteins arising from mutations in the protein interior in experimentally relevant units. Moreover, coupling distances—a measure of the extent of percolation on perturbation—and overall perturbation magnitudes are predicted in a residue-specific manner, enabling a first look at the distribution of energetic couplings in a protein or its changes upon ligand binding. We show specific examples of how the server can be employed to probe for the distribution of local stabilities in a protein, to examine changes in side chain orientations or packing before and after ligand binding, and to predict changes in stabilities of proteins upon mutations of buried residues. The web server is freely available at http://pbl.biotech.iitm.ac.in/pPerturb and supports recent versions of all major browsers.

Helix-helix interaction viewed in an angle frame indicates a role of the size of sidechains in packing

Three-dimensional packing is essential for protein folding. Homology alignment and de nova methods have been used to develop algorithms that would predict protein packing in proteins. Even though hydrophobic cores of proteins are formed through the packing of side chains of hydrophobic residues, the stereological arrangement of individual amino acids and in three dimensional hydrophobic cores is rather difficult to determine. In order to simplify the description of a protein core packing, in this manuscript, we have determined the relationship between angles, distances, and residue usage between two helices. Data from this analysis demonstrate that the contact residues determine the angle and the distance between two helical elements. This approach provides a means to predict the three-dimensional packing of helices in proteins and allows for an understanding of the biological interaction within proteins and among proteins based on surface contact residue parameters.

Toward the Accuracy and Speed of Protein Side-Chain Packing: ASystematic Study on Rotamer Libraries.

Protein rotamers refer to the conformational isomers taken by the side-chains of amino acids to accommodate specific structural folding environments. Since accurate modeling of atomic interactions is difficult, rotamer information collected from experimentally solved protein structures is often used to guide side-chain packing in protein folding and sequence design studies. Many rotamer libraries have been built in the literature but there is little quantitative guidance on which libraries should be chosen for different structural modeling studies. Here, we performed a comparative study of six widely used rotamer libraries and systematically examined their suitability for protein folding and sequence design in four aspects: (1) side-chain match accuracy, (2) side-chain conformation prediction, (3) de novo protein sequence design, and (4) computational time cost. We demonstrated that, compared to the backbone-dependent rotamer libraries (BBDRLs), the backbone-independent rotamer libraries (BBIRLs) generated conformations that more closely matched the native conformations due to the larger number of rotamers in the local rotamer search spaces. However, more practically, using an optimized physical energy function incorporated into a simulated annealing Monte Carlo searching scheme, we showed that utilization of the BBDRLs could result in higher accuracies in side-chain prediction and higher sequence recapitulation rates in protein design experiments. Detailed data analyses showed that the major advantage of BBDRLs lies in the energy term derived from the rotamer probabilities that are associated with the individual backbone torsion angle subspaces. This term is important for distinguishing between amino acid identities as well as the rotamer conformations of an amino acid. Meanwhile, the backbone torsion angle subspace-specific rotamer search drastically speeds up the searching time, despite the significantly larger number of total rotamers in the BBDRLs. These results should provide important guidance for the development and selection of rotamer libraries for practical protein design and structure prediction studies.